Mit2 092 f09_lec07
0 likes382 views
This document outlines the finite element analysis process for solids and fluids in a lecture format. It discusses the derivation of key formulas, the application of interpolation functions, and the criteria for achieving accurate analytical solutions within finite element methods. The importance of equilibrium, compatibility, and the refinement of mesh or interpolation order to improve accuracy is emphasized.
Mit2 092 f09_lec07
- 1. • • • • z }| { Z 2.092/2.093 — Finite Element Analysis of Solids & Fluids I Fall ‘09 Lecture 7 -Finite Element Solution Process, cont’d Prof. K. J. Bathe MIT OpenCourseWare Recall that we have established for the general system some expressions: KU = R (1) K = Σ K(m) ; R = RB + RS mZ K(m) = B(m)T C(m)B(m)dV (m) (2) V (m) Then, since C(m)T = C(m), Z T K(m)T = B(m)T C(m)T B(m)T dV (m) = K(m) V (m) Therefore, K is symmetric! Z R(m) B = H(m)T fB(m)dV (m) (3) V (m) For statics, we use f B(m), but for linear dynamics we must modify it: u¨(m) fB(m) = f˜B(m) − ρ(m) H(m)U ¨ (4) where f˜B(m) is the body force vector excluding inertial loads, and accelerations are interpolated in the same way as displacements. Now, Eq. (1) becomes MU ¨+ KU = R M =Σ M (m) ; M(m) = H(m)T ρ(m)H(m)dV (m) (5) mV (m) We need the initial conditions for U (t = 0) and U˙ (t = 0). We will calculate only the element matrices corresponding to the element nodes. Example Reading assignment: Example 4.5 ⎡ ⎤ ⎡ ⎤ ⎡ ⎤ 0 000 0 K = ⎣ 0 ⎦ + ⎣ 0 ⎦ = ⎣⎦ •• ×ו• + ×× 0000 0 ×× ×× The nonzero elements in the first term correspond to element 1, and those in the second term correspond to element 2. 1
- 2. # | {z } |{z} | {z } # Lecture 7 Finite Element Solution Process, cont’d 2.092/2.093, Fall ‘09 Example: Two-Dimensional Analysis Reading assignment: Example 4.6 This is a thin plate, hence use plane stress assumptions: ⎤ ⎡ 1 ν 0 C = E 1 − ν2 ⎢⎢⎣ ⎥⎥⎦ ν 10 00 1−ν 2 ⎞⎛⎞ ⎛⎞ ⎛ ∂u εxx τxx∂x ∂v ⎜⎜⎝ ⎟⎟⎠ = ⎜⎜⎝ ⎟⎟⎠ , ⎜⎜⎝ ⎟⎟⎠ εyy τyy ∂y ∂u + ∂v ∂x γxy τxy∂y ⎤⎡ u1 u2 u3 u4 v1 v2 v3 v4 8×1 ⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎣ ⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎦ ⎞ ⎛ 4 Σhiui i 4 u(x, y) v(x, y) ⎜⎝ ⎟⎠ = H= Σhivi i 2×8 2×1 h1 h2 h3 h4 0000 H = 0000 h1 h2 h3 h4 hi are called “interpolation functions,” and are defined as 1 at node i and zero at all other nodes. 2
- 3. Lecture 7 Finite Element Solution Process, cont’d 2.092/2.093, Fall ‘09 h1(x, y)= 14 (1 + x)(1 + y) h2(x, y)= 1 (1 − x)(1 + y) 4 h3(x, y)= 1 (1 − x)(1 − y) 4 h4(x, y)= 1 (1 + x)(1 − y) 4 3
- 4. Z Z Lecture 7 Finite Element Solution Process, cont’d 2.092/2.093, Fall ‘09 Note that h1 + h2 + h3 + h4 = 1. The interpolation function hi fulfills the following conditions: Constant strain • • Rigid body motion • Compatibility ⎡⎤ h1,x h2,x h3,x h4,x 0000 ⎢⎥ B = ⎣ 0000 h1,y h2,y h3,y h4,y ⎦ h1,y h2,y h3,y h4,y h1,x h2,x h3,x h4,x Next, we define F (m), the element nodal forces: F (m) = B(m)T τ (m)dV (m) (A) V (m) We assume U has already been calculated, so τ (m) = C(m)B(m)U (B) Two Important Properties I. At any node, the sum of the element nodal point forces is in equilibrium with the externally applied nodal loads; that is: ΣF (m) = R ;ΣF (m) = F mm Why do we have this property? Use (A) and (B): Σ B(m)T C(m)B(m)dV (m) U = KU = R mV (m) II. Each element (m) is in equilibrium under its nodal forces. Proof: ε(m)T =0 Zz r }| { U T F (m) = U T B(m)T τ (m)dV (m) rr V (m) In this expression, UrT represents a rigid body motion. Now we use rigid body motions corresponding to translations and rotations. 4
- 5. Lecture 7 Finite Element Solution Process, cont’d 2.092/2.093, Fall ‘09 In general, to obtain the exact analytical solution, the following should be satisfied: • Differential equilibrium both in the volume, and on the surface • Compatibility Stress-strain law • For a properly formulated finite element procedure, the displacement interpolation functions ensure that compatibility is satisfied. The stress-strain law is also enforced; however, differential equilibrium is not satisfied at every point in the continuum. Rather, equilibrium is only satisfied at every node and on each element (see properties I and II above). As such, the finite element solution is approximate, and this approximation becomes more accurate as the mesh is refined, or the interpolation order of the element is increased. 5
- 6. MIT OpenCourseWare http://ocw.mit.edu 2.092 / 2.093 Finite Element Analysis of Solids and Fluids I Fall 2009 For information about citing these materials or our Terms of Use, visit: http://ocw.mit.edu/terms.