EN221 Summary 1 Tensor Stuff • Divergence : T ∇ · u = ∂iui ∇ · u = u nda, ZR Z∂R T ∇ · T = ∂iTije j ∇ · T = T nda. ZR Z∂R ∇ ⊗ u = Jacobian ∇ ⊗ u = u ⊗ nda, ZR Z∂R (matrix divergence: columns stay separate) • Box product : [ a , b, c] = a · ( b ∧ c) • Levi-Civita tensor: δi, 1 δj, 1 δk, 1 1 ( ijk) an even permut. of ( 123) , εijk = det δi, 2 δj, 2 δk, 2 = [ ei , e j , ek ] = − 1 ( ijk) an odd permut. of ( 123) , δ δ δ 0 if not. i, 3 j, 3 k, 3 e j ∧ ek = εijkei. det( abc) = εijkaibjck εijkεilm = δjl δkm − δjmδkl , εijkεijl = 2δkl , εijkεijk = 6 • Principal Invariants: IA = λ 1 + λ2 + λ 3 = ([ Aa, b, c] + [ a , Ab, c] + [ a , b, Ac])/[ a , b, c] = tr A, 1 II = λ λ + λ λ + λ λ = ([ Aa , Ab, c] + [ Aa , b, Ac] + [ a , Ab, Ac])/[ a , b, c] = [ tr2 A − tr A2 ] , A 1 2 2 3 1 3 2 IIIA = λ 1 λ 2 λ 3 = [ Aa , Ab, Ac]/[ a , b, c] = det A. T • Adjugate/Cofactor of a Tensor: A∗ ( a ∧ b) = ( Aa) ∧ ( Ab) ⇒ A∗ = det A( A− ) . 1 ∂t det A( t) = det A tr(( ∂tA) A− ) • Tensor Product : TO ⊗ FROM T ei ⊗ e j = eie j ( u ⊗ v) a = u( v · a) ( u ⊗ v)( w ⊗ x) = v · w( u ⊗ x) ( u ⊗ v) A = u ⊗ ( ATv) • Skewsymmetric matrices: Rotation around axis Ω given by orthogonal matrix Q( t) . T x˙ = Q˙ x ⇒ ∂t( Q Q) = 0, W = Q˙ Q T, W = − WT. Wx = Ω ∧ x. 2 Kinematics 2.1 Static • Reference and deformed configurations. 1 2 Section 2 • Deformation gradient : assumed regular. J = det F 0. x( X) = X + u( X) , F = ∇X ⊗ x( X) , 1 F− = ∂x jXβEβ ⊗ ei • Isochoric : J = 1 . • Polar decomposition : ◦ F = RU, T 2 1 F F = U , R = FU − . ◦ F = V R. Features: ◦ Is unique. ◦ R is rotation of principal axes. ◦ R average of all rotations. ◦ Principal axes of V are Rui. ◦ σ( V) = σ( U) . ◦ R = vk ⊗ uk . ◦ F = λ kvk ⊗ uk . • Left/Right Cauchy-Green Deformation Tensor: FFT/FTF SPD. • Strain : 1 E = ( FTF − Id) ( Lagrangean: | dx | 2 − | dX | 2 = 2dX · EdX) , 2 1 E = ( Id − F TF 1 ) ( Eulerian: | dx | 2 − | dX | 2 = 2dx · E dx) . ′ 2 − − ′ • Stretch: λ( M) = ( M · FTFM) 1/2 = | UM | . Has local maxima and minima when M is an eigenvector of U. • Transformation of area elements: nda = F∗ NdA • Deformation gradient in cylindrical coordinates: Given r θ = f( R, Θ,Z) , z we have 1 F = ∂ x ⊗ E + ∂ x ⊗ E + ∂ x ⊗ E . R R R Θ Θ Z Z Also expressible as mixed tensor from E( R, Θ ,Z ) to E( r,θ,z) : ∂r 1 ∂r ∂r ∂R R ∂Θ ∂z r ∂θ r ∂θ r ∂θ F = . 1 ∂R R ∂Θ 1 ∂z ∂z 1 ∂z ∂z ∂R R ∂Θ ∂Z Caveat for mixed tensors: tr( F) Fii. However det, V, U as usual. Also works for spherical basis, but more complicated. Kinematics 3 2.1.1 Static Examples 1 • Pure shear: F = λe1 ⊗ e1 + λ − e2 ⊗ e2 . • Simple shear: F = Id + λe1 ⊗ e 2 . • Pure bending: x ( R − Y) sin α( x) y = R − ( R − Y) cos α( x) , J = ( R − Y) α ′. z Z • Tension and torsion : X cos α λZ − Y sin α λZ x √λ l √λ l X α Y α y = sin λZ + cos λZ . √λ l √λ l z λZ • Turning a cylinder inside out . 2.2 Dynamic • Steady motion : ∂/∂t v( x , t) = 0. • Material/Lagrangean POV : focus on particle, expressions in terms of X and t → Solids. • Spatial/Eulerian POV : focus on point in space, expressions in terms of x and t → Fluids. • Lines: ◦ Path line : Curve traced by a fixed particle. ◦ Streamlines: Field lines of velocity in Eulerian POV. Both coincide under steady motion. • Material derivative : ∂ϕ ϕ˙ = + ∇xϕ · v , ∂t ∂w w˙ = + ( ∇x ⊗ w) v , ∂t ∂T T˙ = + ( ∇x ⊗ T) v. ∂T • Acceleration : a = v˙ . • Velocity gradient : L = ∇x ⊗ v ⇒ F˙ = LF (chain rule). F requires a “reference state”, L does not. • dx˙ = F˙ dX = LFdX = Ldx. Assume dx = m| dx | . | dx | Strain rate : • = m · Lm = m · Dm | dx | m˙ = Lm − m( m · Lm) • Stretch and Spin : L = D + W, D = DT, W = − WT. D11 : stretching rate of a line element along the 1-direction D12: (roughly) change in angle between the 1- and 2-direction. Principal axes pi of D are rigidly rotating about 1 ω = curl v 2 with Wpi = ω × pi. • Vorticity : curl v = 2 · angular velocity. (Letter here is also ω.) 4 Section 3 • J˙ = J tr L = J div v. • Integrals over moving contours: v · dx = v( x , t) · FdX ICt ICR d v · dx = ( v˙ ( x , t) · F + v( x , t) · LF) dX dt ICt ICR = v˙ ( x , t) + LTv( x , t)dx ICt T • Integrals over moving surfaces: Similar, taking into account that F∗ = JF− . d u · nds = ( u˙ + u tr( L) − Lu) · nds. dt ZSt ZSt • Integrals over moving volumes/Reynolds’ Transport Theorem : d ϕ( x)dv = ϕ˙ + ϕ tr( L)dv. dt ZR t ZRt Observe that tr( L) = div v, which is zero in the incompressible case. • Circulation : v · dx = curl v · ds = ω · ds ICt ZSt ZSt 1 LTv = ∇v2 2 if circulation-preserving d 0 = v · dx = v˙ ( x , t) + LTv( x , t)dx dt ICt ICt 1 = v˙ ( x , t)dx + ∇v2 dx ICt ICt 2 = a( x , t)dx ICt = curl a · ds ZSt If a = ∇ψ, then the motion is circulation-preserving. If circulation-preserving, then curl a = ω˙ + ωtr( L) − Lω = 0. Then consider the product rule on d 1 ( JF ω) = = 0 dt − to find Cauchy’s vorticity formula : 1 ω = Fω . J ref Field lines of vorticity are vortex lines. If the motion is circulation-preserving, these are material curves. 3 Balance Laws and Field Equations • Conservation of Mass: Assumption: Jρ = ρref ( referential) . Balance Laws and Field Equations 5 Therefore, ρ˙ J + ρJ˙ = 0 ρ˙ J + ρJ div v = 0 ρ˙ + ρ div v = 0 ∂tρ + div( ρv) = 0. • Transport Theorem with density : d ρΦdv = ρΦ˙ dv. dt ZRt ZRt • Linear Momentum : M = ρv ◦ Stress vector: t( n) is force/unit area. ◦ Balance law : d ρv = ρv˙ = ρbdv + t( n) da. dt ZRt ZRt ZR t Z∂Rt T ◦ Stress tensor/Cauchy’s Theorem : σ n = t( n) . Derivation: − t( n) = − t( n) by pillbox and balance law. − − Tetrahedron argument: n the general normal of the coordinate-system-boxed tetrahe- dron. Then other faces ai = ani, where a is the area of the complicated face. Volume h a/3. Apply 1/a · balance law, let h → 0. Assume continuity, derive linear depen- dence by assuming values are locally constant. ◦ Updated balance law: ρa = ρb + ∇ · σ ZR t ZR t ◦ Field equations: ρrefx¨ = ∇X · s + ρrefb ( referential) ρa = ∇x · σ + ρb ( spatial) 1 T T • Nominal Stress/Conjugate stress: s = JF− σ. ( σ nda = s NdA can be directly verified.) Also called Piola-Kirchoff stress. s T is 2nd Piola-Kirchoff stress. • Angular Momentum : H = ρx ∧ v ◦ Non-polar material : no contact torques. ◦ Balance law : ( ) d ∗ ρx ∧ v = ρx ∧ v˙ = ρ( x ∧ b + c)dv + x ∧ t( n) da dt ZRt ZRt ZRt Z∂Rt c is body torque. Equality ( ∗ ) follows because x-derivatives vanish once ∧ v is applied. ◦ Subsituting Cauchy’s Theorem into the balance law gives x ∧ ( ∇x · σ + ρb) = ρ( x ∧ b)dv + x ∧ σnda ZRt ZRt Z∂Rt x ∧ ( ∇x · σ) x ∧ σnda ZR t Z∂R t T View in component form, apply Gauß, derive εijkσji = 0 ⇒ σ = σ . ◦ Field equations: s TFT = F s ( referential) σT = σ ( spatial) 6 Section 3 • Vector identities: 1 ( A · ∇) A = ∇| A | 2 + ( ∇ × A) × A 2 ( A · ∇) A = ( ∇ ⊗ A) A. Use these identities to rewrite the v˙ as ∇( v 2 ) for irrotational flow. • Types of fluid flow : ◦ Inviscid : σ = − p Id ⇒ div σ = − ∇p. ◦ Incompressible : ρ˙ = 0 or div v = 0. ◦ Steady : ∂tv = 0. ρ˙ = v · ∇ ρ. ◦ Irrotational : ω = 0 or v = ∇ϕ. ◦ Elastic : p( ρ) ◦ Ideal=incompressible : div v = 0, J = 1 . • Rayleigh-Plesset equation : Begin with deformation of spherical shell (with extent!), assume J ≡ 1 . Derive ODE. • Conservative potentials: b = − ∇β ◦ Elastic or ideal flow here is circulation preserving, i.e. a = − ∇something. − Have 1 a = − ∇p( ρ) + b ρ 1 = − p ( ρ) ∇ ρ − ∇ β. ρ ′ − Define ρ 1 ε( ρ) = p′( ρ′)dρ ′ Z0 ρ ′ ⇒ ∇ε( ρ) = ε ′( ρ) ∇ρ − Therefore a = − ∇( ε( ρ) + β) . − For ideal fluid substitute p/ ρ0 for ε. ◦ Bernoulli’s Theorem : − Flow irrotational (i.e. v = − ∇ ϕ): v2 ∇ ∂tϕ + + ε( ρ) + β = 0. 2 Proof: Just rewrite, obtaining v2 /2 from second term of material derivative. − Flow steady: v2 • + ε( ρ) + β = 0, 2 i.e. this quantity is constant along streamlines. Proof: Exploit v · v˙ = v · ∇( v2 ) − Flow both irrotational and steady: v2 ∇ + ε( ρ) + β = 0. 2 Balance Laws and Field Equations 7 • Acoustic wave equation : ◦ Assume ρ = ρ0 + δρ, | v | ≪ 1 , |∇v | ≪ 1 . ◦ Start with ∂t ( ∇ · v) , use cons. of. momentum without second order term, cons. of mass as ∂tρ + ρ0div( v) = 0. 2 2 ◦ δρtt = c ∇ δρ, with c = ∂ρp. p • Mach number: assume steadiness b = 0, use v · v˙ in terms of c2 . v2 v · ( ρv) = v · v˙ ρ 1 − • c2 m 4 ◦ Supersonic nozzle m < 1 , m > 1 . • Conservation of Energy : ◦ Balance law : d K( R ) = − S( R ) + P( R ) dt t t t d 1 ρv · vdv = − tr( σD)dv + ρb · v + σn · vda dt 2 Z ZR t ZRt Z∂Rt Kinetic energy K( t) Stress power S( Rt ) Power supplied P ( Rt ) Proof: Multiply Equation of Motion by v, integrate by parts in the σ term.