Stable Orthotropic Materials

Total Page:16

File Type:pdf, Size:1020Kb

Stable Orthotropic Materials Eurographics/ ACM SIGGRAPH Symposium on Computer Animation (2014) Vladlen Koltun and Eftychios Sifakis (Editors) Stable Orthotropic Materials Yijing Li and Jernej Barbicˇ University of Southern California, United States Abstract Isotropic Finite Element Method (FEM) deformable object simulations are widely used in computer graphics. Several applications (wood, plants, muscles) require modeling the directional dependence of the material elas- tic properties in three orthogonal directions. We investigate orthotropic materials, a special class of anisotropic materials where the shear stresses are decoupled from normal stresses. Orthotropic materials generalize trans- versely isotropic materials, by exhibiting different stiffnesses in three orthogonal directions. Orthotropic materials are, however, parameterized by nine values that are difficult to tune in practice, as poorly adjusted settings easily lead to simulation instabilities. We present a user-friendly approach to setting these parameters that is guaran- teed to be stable. Our approach is intuitive as it extends the familiar intuition known from isotropic materials. We demonstrate our technique by augmenting linear corotational FEM implementations with orthotropic materials. Categories and Subject Descriptors (according to ACM CCS): I.3.6 [Computer Graphics]: Methodology and Techniques—Interaction Techniques, I.6.8 [Simulation and Modeling]: Types of Simulation—Animation 1. Introduction parameters, by parameterizing the six Poisson’s ratios using a stable one-dimensional parameter family, similar to the in- Simulation of three-dimensional solid deformable models tuition from isotropic simulation. This makes it possible to is important in many applications in computer graphics, easily augment existing simulation solvers with stable and robotics, special effects and virtual reality. Most applications intuitive anisotropic effects. We support large deformations in these fields have been limited to isotropic materials, i.e., using corotational Finite Element Method simulation. materials that are equally elastic in all directions. Many real materials are, however, stiffer in some directions than oth- 2. Related Work ers. The space of such anisotropic materials is vast and not easy to navigate, tune or control. In this paper, we study or- Anisotropic materials are discussed in many references, see, thotropic materials. Orthotropic materials exhibit different e.g. [Bow11]. Transversely isotropic hyperelastic materials stiffness in three orthogonal directions; formally, they pos- were presented by Bonet and Burton [BB98]. Picinbono sess three orthogonal planes of rotational symmetry. They et al. [PDA01] proposed a non-linear FEM model to sim- form an intuitive subset of all anisotropic materials, as they ulate soft tissues with large deformations and transversely generalize the familiar isotropic, and transversely isotropic, isotropic behavior. Thije et al. [TTAH07] addressed the in- materials to materials with three different stiffness values stabilities that occur under strong anisotropy, and provided a in some three orthogonal directions. Although simpler than simple updated Lagrangian FEM scheme to handle the prob- fully general anisotropic materials, orthotropic materials still lem. Picinbono et al. [PLDA00] described a surgery simu- require tuning nine independent parameter values. In prac- lator that can model linear transversely isotropic materials tice, this task is difficult due to the large number of param- at haptic rates and also presented a nonlinear transversely eters and because many of the settings lead to unstable sim- isotropic model for medical simulation [PDA03]. Sermesant ulations in a non-obvious way. In this paper, we study or- et al. [SCD∗01,SDA06] and Talbot et al. [TMD∗13] adopted thotropic materials from the point of view of practical sim- a transversely isotropic material in constructing an electro- ulation in computer graphics and related fields. We demon- mechanical model of the heart. Allard et al. [AMC∗09] used strate how to intuitively and stably tune orthotropic material a 2D anisotropic material to simulate thin soft tissue tearing, c The Eurographics Association 2014. Yijing Li & Jernej Barbiˇc/ Stable Orthotropic Materials Comas et al. [CTA∗08] implemented a transversely isotropic volume of tet e; and Be is a 6 × 12 matrix determined by the visco-hyperelastic model on the GPU, and Teran [TSB∗05] initial shape of tet e (see [Sha90] or [MG04]). and Sifakis [SNF05, Sif07] simulated human muscles with a transversely isotropic, quasi-incompressible model. Irv- Unlike isotropic materials that are parameterized by a sin- ing et al. [ITF04] proposed a robust, large-deformation in- gle Young’s modulus and Poisson’s ratio, orthotropic mate- vertible simulation method and demonstrated it with trans- rials have three different Young’s moduli E1;E2;E3, one for versely isotropic models. Previous anisotropic applications each orthogonal direction, and six Poisson’s ratios ni j; for in computer graphics focused on transversely isotropic ma- i 6= j; only three of which are independent. Young’s mod- terials where two directions have equal stiffness, leading ulus Ei gives the stiffness of the material when loaded in to five tunable parameters. We generalize materials to or- orthogonal direction i: Poisson’s ratio ni j gives the contrac- thotropic materials with three distinct stiffnesses in three tion in direction j when the extension is applied in direc- orthogonal directions, and present an intuitive approach to tion i: In a general anisotropic material, both the normal and tune the resulting nine parameters. To the best of our knowl- shear components of strain affect both the normal and shear edge, we are first work in computer graphics to analyze components of stress, i.e., matrix C is dense. In orthotropic orthotropic materials in substantial detail. Previous papers materials, however, the normal and shear components are on orthotropic materials in engineering assumed that the decoupled: normal stresses only cause normal strains, and nine orthotropic parameters are given or measured from shear stresses only cause shear strains. Furthermore, individ- real materials [VBCW81], whereas we provide an intuitive ual shear stresses in the 12;23;31 planes are decoupled from way for the users to tune them and ensure they are stable. each other: strain ei j (i 6= j) only depends on stress si j via a Our work uses corotational linear FEM materials introduced scalar parameter (shear modulus) mi j: Under these assump- in [MG04]. Construction of the stiffness matrix for linear tions, the elasticity tensor has 9 free parameters and takes a FEM materials can also be found, for example, in [Sha90]. block-diagonal form. It is easiest to first state its inverse 2 1 − n21 − n31 0 0 0 3 E1 E2 E3 3. Orthotropic Materials 6− n12 1 − n32 0 0 0 7 6 E1 E2 E3 7 6− n13 − n23 1 0 0 0 7 We now introduce orthotropic materials. Given the defor- 6 E1 E2 E3 7 Sortho = 6 7: (4) mation gradient F; the Green-Lagrange strain is defined as 6 0 0 0 1 0 0 7 6 m12 7 e3×3 = (FT F − I)=2; and the Cauchy stress s 3×3 gives 6 0 0 0 0 1 0 7 4 m23 5 the elastic forces per surface area in a unit direction n; 0 0 0 0 0 1 as s 3×3n [Sha90]. Note that we can operate with Cauchy m31 stresses here as they are equivalent to other forms of stresses The orthotropic elasticity tensor is then (Piola) due to the small-deformation analysis; we achieve large deformations via co-rotational linear FEM [MG04]. A 0 The 6 × 6 elasticity tensor S relates strain e to stress s via C = S −1 = ; for (5) ortho ortho 0 B e = S s; where we have unrolled the 3×3 symmetric matri- 3×3 3×3 2 3 ces e = [ei j]i j and s = [si j]i j into 6-vectors, using the E1(1 − n23n32) E2(n12 + n32n13) E3(n13 + n12n23) 12, 23, 31 ordering of the shear components as in [Sha90]: A = ¡4E1(n21 + n31n23) E2(1 − n13n31) E3(n23 + n21n13)5; E ( + ) E ( + ) E ( − ) T 1 n31 n21n32 2 n32 n12n31 3 1 n12n21 e = [e11 e22 e33 2e12 2e23 2e31] ; (1) (6) T s = [s11 s22 s33 s12 s23 s31] : (2) 2 3 m12 0 0 Components 11;22;33 are called normal components, B = 4 0 m23 0 5; and (7) whereas 12;23;31 are referred to as shear components. The 0 0 m31 inverse elasticity tensor C = S −1 relates s to e; via s = 1 ¡ = : (8) C e: The elasticity tensor must be symmetric and therefore it 1 − n12n21 − n23n32 − n31n13 − 2n21n32n13 has 21 independent entries for a general anisotropic material, Equations4 and5 give elasticity tensors with respect to the world coordinate axes. A general orthotropic material, 2C C C C C C 3 11 12 13 14 15 16 however, assumes the block-diagonal form given in Equa- 6 C22 C23 C24 C25 C267 6 7 tions4 and5 only in a special orthogonal basis, given by 6 C33 C34 C35 C367 C = 6 7: (3) the three principal axes where the stiffnesses are E1;E2;E3: 6 C44 C C 7 6 45 467 In other bases (including world-coordinate axes), its form 4 Sym: C55 C565 looks generic, as in Equation3. Therefore, to model or- C 66 thotropic materials whose principal axes are not aligned with Once C is known, the stiffness matrix for a linear tetrahe- the world axes, we need to convert elasticity tensors from dral element is computed as Ke = V eBeT C eBe; where V e is one basis to another. For a basis given by a rotation Q; the c The Eurographics Association 2014. Yijing Li & Jernej Barbiˇc/ Stable Orthotropic Materials elasticity tensor C transforms as follows: practice, it is very tedious to tune these parameters, as sub- (1) (2) optimal values easily cause the simulation to explode, or in- T K 2K troduce undue stiffness or other poor simulation behavior.
Recommended publications
  • 10-1 CHAPTER 10 DEFORMATION 10.1 Stress-Strain Diagrams And
    EN380 Naval Materials Science and Engineering Course Notes, U.S. Naval Academy CHAPTER 10 DEFORMATION 10.1 Stress-Strain Diagrams and Material Behavior 10.2 Material Characteristics 10.3 Elastic-Plastic Response of Metals 10.4 True stress and strain measures 10.5 Yielding of a Ductile Metal under a General Stress State - Mises Yield Condition. 10.6 Maximum shear stress condition 10.7 Creep Consider the bar in figure 1 subjected to a simple tension loading F. Figure 1: Bar in Tension Engineering Stress () is the quotient of load (F) and area (A). The units of stress are normally pounds per square inch (psi). = F A where: is the stress (psi) F is the force that is loading the object (lb) A is the cross sectional area of the object (in2) When stress is applied to a material, the material will deform. Elongation is defined as the difference between loaded and unloaded length ∆푙 = L - Lo where: ∆푙 is the elongation (ft) L is the loaded length of the cable (ft) Lo is the unloaded (original) length of the cable (ft) 10-1 EN380 Naval Materials Science and Engineering Course Notes, U.S. Naval Academy Strain is the concept used to compare the elongation of a material to its original, undeformed length. Strain () is the quotient of elongation (e) and original length (L0). Engineering Strain has no units but is often given the units of in/in or ft/ft. ∆푙 휀 = 퐿 where: is the strain in the cable (ft/ft) ∆푙 is the elongation (ft) Lo is the unloaded (original) length of the cable (ft) Example Find the strain in a 75 foot cable experiencing an elongation of one inch.
    [Show full text]
  • Static Analysis of Isotropic, Orthotropic and Functionally Graded Material Beams
    Journal of Multidisciplinary Engineering Science and Technology (JMEST) ISSN: 2458-9403 Vol. 3 Issue 5, May - 2016 Static analysis of isotropic, orthotropic and functionally graded material beams Waleed M. Soliman M. Adnan Elshafei M. A. Kamel Dep. of Aeronautical Engineering Dep. of Aeronautical Engineering Dep. of Aeronautical Engineering Military Technical College Military Technical College Military Technical College Cairo, Egypt Cairo, Egypt Cairo, Egypt [email protected] [email protected] [email protected] Abstract—This paper presents static analysis of degrees of freedom for each lamina, and it can be isotropic, orthotropic and Functionally Graded used for long and short beams, this laminated finite Materials (FGMs) beams using a Finite Element element model gives good results for both stresses and deflections when compared with other solutions. Method (FEM). Ansys Workbench15 has been used to build up several models to simulate In 1993 Lidstrom [2] have used the total potential different types of beams with different boundary energy formulation to analyze equilibrium for a conditions, all beams have been subjected to both moderate deflection 3-D beam element, the condensed two-node element reduced the size of the of uniformly distributed and transversal point problem, compared with the three-node element, but loads within the experience of Timoshenko Beam increased the computing time. The condensed two- Theory and First order Shear Deformation Theory. node system was less numerically stable than the The material properties are assumed to be three-node system. Because of this fact, it was not temperature-independent, and are graded in the possible to evaluate the third and fourth-order thickness direction according to a simple power differentials of the strain energy function, and thus not law distribution of the volume fractions of the possible to determine the types of criticality constituents.
    [Show full text]
  • Constitutive Relations: Transverse Isotropy and Isotropy
    Objectives_template Module 3: 3D Constitutive Equations Lecture 11: Constitutive Relations: Transverse Isotropy and Isotropy The Lecture Contains: Transverse Isotropy Isotropic Bodies Homework References file:///D|/Web%20Course%20(Ganesh%20Rana)/Dr.%20Mohite/CompositeMaterials/lecture11/11_1.htm[8/18/2014 12:10:09 PM] Objectives_template Module 3: 3D Constitutive Equations Lecture 11: Constitutive relations: Transverse isotropy and isotropy Transverse Isotropy: Introduction: In this lecture, we are going to see some more simplifications of constitutive equation and develop the relation for isotropic materials. First we will see the development of transverse isotropy and then we will reduce from it to isotropy. First Approach: Invariance Approach This is obtained from an orthotropic material. Here, we develop the constitutive relation for a material with transverse isotropy in x2-x3 plane (this is used in lamina/laminae/laminate modeling). This is obtained with the following form of the change of axes. (3.30) Now, we have Figure 3.6: State of stress (a) in x1, x2, x3 system (b) with x1-x2 and x1-x3 planes of symmetry From this, the strains in transformed coordinate system are given as: file:///D|/Web%20Course%20(Ganesh%20Rana)/Dr.%20Mohite/CompositeMaterials/lecture11/11_2.htm[8/18/2014 12:10:09 PM] Objectives_template (3.31) Here, it is to be noted that the shear strains are the tensorial shear strain terms. For any angle α, (3.32) and therefore, W must reduce to the form (3.33) Then, for W to be invariant we must have Now, let us write the left hand side of the above equation using the matrix as given in Equation (3.26) and engineering shear strains.
    [Show full text]
  • Orthotropic Material Homogenization of Composite Materials Including Damping
    View metadata, citation and similar papers at core.ac.uk brought to you by CORE provided by Lirias Orthotropic material homogenization of composite materials including damping A. Nateghia,e, A. Rezaeia,e,, E. Deckersa,d, S. Jonckheerea,b,d, C. Claeysa,d, B. Pluymersa,d, W. Van Paepegemc, W. Desmeta,d a KU Leuven, Department of Mechanical Engineering, Celestijnenlaan 300B box 2420, 3001 Heverlee, Belgium. E-mail: [email protected]. b Siemens Industry Software NV, Digital Factory, Product Lifecycle Management – Simulation and Test Solutions, Interleuvenlaan 68, B-3001 Leuven, Belgium. c Ghent University, Department of Materials Science & Engineering, Technologiepark-Zwijnaarde 903, 9052 Zwijnaarde, Belgium. d Member of Flanders Make. e SIM vzw, Technologiepark Zwijnaarde 935, B-9052 Ghent, Belgium. Introduction In the multiscale analysis of composite materials obtaining the effective equivalent material properties of the composite from the properties of its constituents is an important bridge between different scales. Although, a variety of homogenization methods are already in use, there is still a need for further development of novel approaches which can deal with complexities such as frequency dependency of material properties; which is specially of importance when dealing with damped structures. Time-domain methods Frequency-domain method Material homogenization is performed in two separate steps: Wave dispersion based homogenization is performed to a static step and a dynamic one. provide homogenized damping and material
    [Show full text]
  • 2.1.5 Stress-Strain Relation for Anisotropic Materials
    ,/... + -:\ 2.1.5Stress-strain Relation for AnisotropicMaterials The generalizedHook's law relatingthe stressesto strainscan be writtenin the form U2, 13,30-32]: o,=C,,€, (i,i =1,2,.........,6) Where: o, - Stress('ompttnents C,, - Elasticity matrix e,-Straincomponents The material properties (elasticity matrix, Cr) has thirty-six constants,stress-strain relations in the materialcoordinates are: [30] O1 LI CT6 a, l O2 03_ a, -t ....(2-6) T23 1/ )1 T32 1/1l 'Vt) Tt2 C16 c66 I The relations in equation (3-6) are referred to as characterizing anisotropicmaterials since there are no planes of symmetry for the materialproperties. An alternativename for suchanisotropic material is a triclinic material. If there is one plane of material property symmetry,the stress- strainrelations reduce to: -'+1-/ ) o1 c| ct2 ct3 0 0 crc €l o2 c12 c22 c23 0 0 c26 a: ct3 C23 ci3 0 0 c36 o-\ o3 ....(2-7) 723 000cqqc450 y)1 T3l 000c4sc5s0 Y3l T1a ct6 c26 c36 0 0 c66 1/t) Where the plane of symmetryis z:0. Such a material is termed monoclinicwhich containstfirfen independentelastic constants. If there are two brthogonalplanes of material property symmetry for a material,symmetry will existrelative to a third mutuallyorthogonal plane. The stress-strainrelations in coordinatesaligned with principal material directionswhich are parallel to the intersectionsof the three orthogonalplanes of materialsymmetry, are'. Ol ctl ct2ctl000 tl O2 Ct2 c22c23000 a./ o3 ct3 c23c33000 * o-\ ....(2-8) T23 0 00c4400 y)1 yll 731 0 000css0 712 0 0000cr'e 1/1) 'lt"' -' , :t i' ,/ Suchmaterial is calledorthotropic material.,There is no interaction betweennormal StreSSeS ot r 02, 03) andanisotropic materials.
    [Show full text]
  • Analysis of Deformation
    Chapter 7: Constitutive Equations Definition: In the previous chapters we’ve learned about the definition and meaning of the concepts of stress and strain. One is an objective measure of load and the other is an objective measure of deformation. In fluids, one talks about the rate-of-deformation as opposed to simply strain (i.e. deformation alone or by itself). We all know though that deformation is caused by loads (i.e. there must be a relationship between stress and strain). A relationship between stress and strain (or rate-of-deformation tensor) is simply called a “constitutive equation”. Below we will describe how such equations are formulated. Constitutive equations between stress and strain are normally written based on phenomenological (i.e. experimental) observations and some assumption(s) on the physical behavior or response of a material to loading. Such equations can and should always be tested against experimental observations. Although there is almost an infinite amount of different materials, leading one to conclude that there is an equivalently infinite amount of constitutive equations or relations that describe such materials behavior, it turns out that there are really three major equations that cover the behavior of a wide range of materials of applied interest. One equation describes stress and small strain in solids and called “Hooke’s law”. The other two equations describe the behavior of fluidic materials. Hookean Elastic Solid: We will start explaining these equations by considering Hooke’s law first. Hooke’s law simply states that the stress tensor is assumed to be linearly related to the strain tensor.
    [Show full text]
  • Crack Tip Elements and the J Integral
    EN234: Computational methods in Structural and Solid Mechanics Homework 3: Crack tip elements and the J-integral Due Wed Oct 7, 2015 School of Engineering Brown University The purpose of this homework is to help understand how to handle element interpolation functions and integration schemes in more detail, as well as to explore some applications of FEA to fracture mechanics. In this homework you will solve a simple linear elastic fracture mechanics problem. You might find it helpful to review some of the basic ideas and terminology associated with linear elastic fracture mechanics here (in particular, recall the definitions of stress intensity factor and the nature of crack-tip fields in elastic solids). Also check the relations between energy release rate and stress intensities, and the background on the J integral here. 1. One of the challenges in using finite elements to solve a problem with cracks is that the stress field at a crack tip is singular. Standard finite element interpolation functions are designed so that stresses remain finite a everywhere in the element. Various types of special b c ‘crack tip’ elements have been designed that 3L/4 incorporate the singularity. One way to produce a L/4 singularity (the method used in ABAQUS) is to mesh L the region just near the crack tip with 8 noded elements, with a special arrangement of nodal points: (i) Three of the nodes (nodes 1,4 and 8 in the figure) are connected together, and (ii) the mid-side nodes 2 and 7 are moved to the quarter-point location on the element side.
    [Show full text]
  • Multidisciplinary Design Project Engineering Dictionary Version 0.0.2
    Multidisciplinary Design Project Engineering Dictionary Version 0.0.2 February 15, 2006 . DRAFT Cambridge-MIT Institute Multidisciplinary Design Project This Dictionary/Glossary of Engineering terms has been compiled to compliment the work developed as part of the Multi-disciplinary Design Project (MDP), which is a programme to develop teaching material and kits to aid the running of mechtronics projects in Universities and Schools. The project is being carried out with support from the Cambridge-MIT Institute undergraduate teaching programe. For more information about the project please visit the MDP website at http://www-mdp.eng.cam.ac.uk or contact Dr. Peter Long Prof. Alex Slocum Cambridge University Engineering Department Massachusetts Institute of Technology Trumpington Street, 77 Massachusetts Ave. Cambridge. Cambridge MA 02139-4307 CB2 1PZ. USA e-mail: [email protected] e-mail: [email protected] tel: +44 (0) 1223 332779 tel: +1 617 253 0012 For information about the CMI initiative please see Cambridge-MIT Institute website :- http://www.cambridge-mit.org CMI CMI, University of Cambridge Massachusetts Institute of Technology 10 Miller’s Yard, 77 Massachusetts Ave. Mill Lane, Cambridge MA 02139-4307 Cambridge. CB2 1RQ. USA tel: +44 (0) 1223 327207 tel. +1 617 253 7732 fax: +44 (0) 1223 765891 fax. +1 617 258 8539 . DRAFT 2 CMI-MDP Programme 1 Introduction This dictionary/glossary has not been developed as a definative work but as a useful reference book for engi- neering students to search when looking for the meaning of a word/phrase. It has been compiled from a number of existing glossaries together with a number of local additions.
    [Show full text]
  • Lecture II Polarization Anisotropy Spectrum
    Lecture II Polarization Anisotropy Spectrum Wayne Hu Crete, July 2008 Damping Tail Dissipation / Diffusion Damping • Imperfections in the coupled fluid → mean free path λC in the baryons • Random walk over diffusion scale: geometric mean of mfp & horizon λD ~ λC√N ~ √λCη >> λC • Overtake wavelength: λD ~ λ ; second order in λC/λ • Viscous damping for R<1; heat conduction damping for R>1 1.0 N=η / λC Power λD ~ λC√N 0.1 perfect fluid instant decoupling λ 500 1000 1500 Silk (1968); Hu & Sugiyama (1995); Hu & White (1996) l Dissipation / Diffusion Damping • Rapid increase at recombination as mfp ↑ • Independent of (robust to changes in) perturbation spectrum • Robust physical scale for angular diameter distance test (ΩK, ΩΛ) Recombination 1.0 Power 0.1 perfect fluid instant decoupling recombination 500 1000 1500 Silk (1968); Hu & Sugiyama (1995); Hu & White (1996) l Damping • Tight coupling equations assume a perfect fluid: no viscosity, no heat conduction • Fluid imperfections are related to the mean free path of the photons in the baryons −1 λC = τ_ where τ_ = neσT a is the conformal opacity to Thomson scattering • Dissipation is related to the diffusion length: random walk approximation p p p λD = NλC = η=λC λC = ηλC the geometric mean between the horizon and mean free path • λD/η∗ ∼ few %, so expect the peaks :> 3 to be affected by dissipation Equations of Motion • Continuity k Θ_ = − v − Φ_ ; δ_ = −kv − 3Φ_ 3 γ b b where the photon equation remains unchanged and the baryons follow number conservation with ρb = mbnb • Euler k v_ = k(Θ + Ψ) − π − τ_(v − v ) γ 6 γ γ b a_ v_ = − v + kΨ +τ _(v − v )=R b a b γ b where the photons gain an anisotropic stress term πγ from radiation viscosity and a momentum exchange term with the baryons and are compensated by the opposite term in the baryon Euler equation Viscosity • Viscosity is generated from radiation streaming from hot to cold regions • Expect k π ∼ v γ γ τ_ generated by streaming, suppressed by scattering in a wavelength of the fluctuation.
    [Show full text]
  • Cosmic Microwave Background
    1 29. Cosmic Microwave Background 29. Cosmic Microwave Background Revised August 2019 by D. Scott (U. of British Columbia) and G.F. Smoot (HKUST; Paris U.; UC Berkeley; LBNL). 29.1 Introduction The energy content in electromagnetic radiation from beyond our Galaxy is dominated by the cosmic microwave background (CMB), discovered in 1965 [1]. The spectrum of the CMB is well described by a blackbody function with T = 2.7255 K. This spectral form is a main supporting pillar of the hot Big Bang model for the Universe. The lack of any observed deviations from a 7 blackbody spectrum constrains physical processes over cosmic history at redshifts z ∼< 10 (see earlier versions of this review). Currently the key CMB observable is the angular variation in temperature (or intensity) corre- lations, and to a growing extent polarization [2–4]. Since the first detection of these anisotropies by the Cosmic Background Explorer (COBE) satellite [5], there has been intense activity to map the sky at increasing levels of sensitivity and angular resolution by ground-based and balloon-borne measurements. These were joined in 2003 by the first results from NASA’s Wilkinson Microwave Anisotropy Probe (WMAP)[6], which were improved upon by analyses of data added every 2 years, culminating in the 9-year results [7]. In 2013 we had the first results [8] from the third generation CMB satellite, ESA’s Planck mission [9,10], which were enhanced by results from the 2015 Planck data release [11, 12], and then the final 2018 Planck data release [13, 14]. Additionally, CMB an- isotropies have been extended to smaller angular scales by ground-based experiments, particularly the Atacama Cosmology Telescope (ACT) [15] and the South Pole Telescope (SPT) [16].
    [Show full text]
  • Cross-Linked Polymers and Rubber Elasticity Chapter 9 (Sperling)
    Cross-linked Polymers and Rubber Elasticity Chapter 9 (Sperling) • Definition of Rubber Elasticity and Requirements • Cross-links, Networks, Classes of Elastomers (sections 1-3, 16) • Simple Theory of Rubber Elasticity (sections 4-8) – Entropic Origin of Elastic Retractive Forces – The Ideal Rubber Behavior • Departures from the Ideal Rubber Behavior (sections 9-11) – Non-zero Energy Contribution to the Elastic Retractive Forces – Stress-induced Crystallization and Limited Extensibility of Chains (How to make better elastomers: High Strength and High Modulus) – Network Defects (dangling chains, loops, trapped entanglements, etc..) – Semi-empirical Mooney-Rivlin Treatment (Affine vs Non-Affine Deformation) Definition of Rubber Elasticity and Requirements • Definition of Rubber Elasticity: Very large deformability with complete recoverability. • Molecular Requirements: – Material must consist of polymer chains. Need to change conformation and extension under stress. – Polymer chains must be highly flexible. Need to access conformational changes (not w/ glassy, crystalline, stiff mat.) – Polymer chains must be joined in a network structure. Need to avoid irreversible chain slippage (permanent strain). One out of 100 monomers must connect two different chains. Connections (covalent bond, crystallite, glassy domain in block copolymer) Cross-links, Networks and Classes of Elastomers • Chemical Cross-linking Process: Sol-Gel or Percolation Transition • Gel Characteristics: – Infinite Viscosity – Non-zero Modulus – One giant Molecule – Solid
    [Show full text]
  • Chapter 10: Elasticity and Oscillations
    Chapter 10 Lecture Outline 1 Copyright © The McGraw-Hill Companies, Inc. Permission required for reproduction or display. Chapter 10: Elasticity and Oscillations •Elastic Deformations •Hooke’s Law •Stress and Strain •Shear Deformations •Volume Deformations •Simple Harmonic Motion •The Pendulum •Damped Oscillations, Forced Oscillations, and Resonance 2 §10.1 Elastic Deformation of Solids A deformation is the change in size or shape of an object. An elastic object is one that returns to its original size and shape after contact forces have been removed. If the forces acting on the object are too large, the object can be permanently distorted. 3 §10.2 Hooke’s Law F F Apply a force to both ends of a long wire. These forces will stretch the wire from length L to L+L. 4 Define: L The fractional strain L change in length F Force per unit cross- stress A sectional area 5 Hooke’s Law (Fx) can be written in terms of stress and strain (stress strain). F L Y A L YA The spring constant k is now k L Y is called Young’s modulus and is a measure of an object’s stiffness. Hooke’s Law holds for an object to a point called the proportional limit. 6 Example (text problem 10.1): A steel beam is placed vertically in the basement of a building to keep the floor above from sagging. The load on the beam is 5.8104 N and the length of the beam is 2.5 m, and the cross-sectional area of the beam is 7.5103 m2.
    [Show full text]