- •Copyright Notice
- •Table of Contents
- •Chapter 1. Introduction
- •I. Basic Concepts
- •Examples:
- •Why Finite Element Method?
- •Applications of FEM in Engineering
- •Examples:
- •A Brief History of the FEM
- •FEM in Structural Analysis (The Procedure)
- •Example:
- •Available Commercial FEM Software Packages
- •Objectives of This FEM Course
- •II. Review of Matrix Algebra
- •Linear System of Algebraic Equations
- •Matrix Addition and Subtraction
- •Scalar Multiplication
- •Matrix Multiplication
- •Transpose of a Matrix
- •Symmetric Matrix
- •Unit (Identity) Matrix
- •Determinant of a Matrix
- •Singular Matrix
- •Matrix Inversion
- •Examples:
- •Solution Techniques for Linear Systems of Equations
- •Positive Definite Matrix
- •Differentiation and Integration of a Matrix
- •Types of Finite Elements
- •III. Spring Element
- •One Spring Element
- •Spring System
- •Checking the Results
- •Notes About the Spring Elements
- •Example 1.1
- •Chapter 2. Bar and Beam Elements
- •I. Linear Static Analysis
- •II. Bar Element
- •Stiffness Matrix --- Direct Method
- •Stiffness Matrix --- A Formal Approach
- •Example 2.1
- •Example 2.2
- •Distributed Load
- •Bar Elements in 2-D and 3-D Space
- •2-D Case
- •Transformation
- •Stiffness Matrix in the 2-D Space
- •Element Stress
- •Example 2.3
- •Example 2.4 (Multipoint Constraint)
- •3-D Case
- •III. Beam Element
- •Simple Plane Beam Element
- •Direct Method
- •Formal Approach
- •3-D Beam Element
- •Example 2.5
- •Equivalent Nodal Loads of Distributed Transverse Load
- •Example 2.6
- •Example 2.7
- •FE Analysis of Frame Structures
- •Example 2.8
- •Chapter 3. Two-Dimensional Problems
- •I. Review of the Basic Theory
- •Plane (2-D) Problems
- •Stress-Strain-Temperature (Constitutive) Relations
- •Strain and Displacement Relations
- •Equilibrium Equations
- •Exact Elasticity Solution
- •Example 3.1
- •II. Finite Elements for 2-D Problems
- •A General Formula for the Stiffness Matrix
- •Constant Strain Triangle (CST or T3)
- •Linear Strain Triangle (LST or T6)
- •Linear Quadrilateral Element (Q4)
- •Quadratic Quadrilateral Element (Q8)
- •Example 3.2
- •Transformation of Loads
- •Stress Calculation
- •I. Symmetry
- •Types of Symmetry:
- •Examples:
- •Applications of the symmetry properties:
- •Examples:
- •Cautions:
- •II. Substructures (Superelements)
- •Physical Meaning:
- •Mathematical Meaning:
- •Advantages of Using Substructures/Superelements:
- •Disadvantages:
- •III. Equation Solving
- •Direct Methods (Gauss Elimination):
- •Iterative Methods:
- •Gauss Elimination - Example:
- •Iterative Method - Example:
- •IV. Nature of Finite Element Solutions
- •Stiffening Effect:
- •V. Numerical Error
- •VI. Convergence of FE Solutions
- •Type of Refinements:
- •Examples:
- •VII. Adaptivity (h-, p-, and hp-Methods)
- •Error Indicators:
- •Examples:
- •Chapter 5. Plate and Shell Elements
- •Applications:
- •Forces and Moments Acting on the Plate:
- •Stresses:
- •Relations Between Forces and Stresses
- •Thin Plate Theory ( Kirchhoff Plate Theory)
- •Examples:
- •Under uniform load q
- •Thick Plate Theory (Mindlin Plate Theory)
- •II. Plate Elements
- •Kirchhoff Plate Elements:
- •Mindlin Plate Elements:
- •Discrete Kirchhoff Element:
- •Test Problem:
- •Mesh
- •III. Shells and Shell Elements
- •Example: A Cylindrical Container.
- •Shell Theory:
- •Shell Elements:
- •Curved shell elements:
- •Test Cases:
- •Chapter 6. Solid Elements for 3-D Problems
- •I. 3-D Elasticity Theory
- •Stress State:
- •Strains:
- •Stress-strain relation:
- •Displacement:
- •Strain-Displacement Relation:
- •Equilibrium Equations:
- •Stress Analysis:
- •II. Finite Element Formulation
- •Displacement Field:
- •Stiffness Matrix:
- •III. Typical 3-D Solid Elements
- •Tetrahedron:
- •Hexahedron (brick):
- •Penta:
- •Element Formulation:
- •Solids of Revolution (Axisymmetric Solids)
- •Axisymmetric Elements
- •Applications
- •Chapter 7. Structural Vibration and Dynamics
- •I. Basic Equations
- •A. Single DOF System
- •B. Multiple DOF System
- •Example
- •II. Free Vibration
- •III. Damping
- •IV. Modal Equations
- •V. Frequency Response Analysis
- •VI. Transient Response Analysis
- •B. Modal Method
- •Cautions in Dynamic Analysis
- •Examples
- •Chapter 8. Thermal Analysis
- •Further Reading
Lecture Notes: Introduction to Finite Element Method |
Chapter 6. Solid Elements for 3-D Problems |
Chapter 6. Solid Elements for 3-D Problems
I. 3-D Elasticity Theory
Stress State:
y
z
y , v
σy
τ yx
τ yz
τzy
τzx
σz
z, w
F
x
τxy
σ x
τxz |
x, u |
|
© 1997-2002 Yijun Liu, University of Cincinnati |
138 |
Lecture Notes: Introduction to Finite Element Method |
Chapter 6. Solid Elements for 3-D Problems |
|
|
|
|
|
|
|
|
σ x |
|
|
|
||||
|
|
|
|
|
|
|
|
|
σ |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
σ |
|
y |
|
|
|
|
|
|
|
σ={σ }= |
|
|
z |
|
|
, |
||||||
|
|
|
|
|
|
|
|
|
|||||||
|
|
|
|
|
|
|
|
|
τ xy |
|
|
|
|||
|
|
|
|
|
|
|
|
τ |
yz |
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
τ zx |
|
|
|||||
Strains: |
|
|
|
|
|
|
|
|
|
||||||
|
|
|
|
|
|
|
ε x |
|
|
|
|
||||
|
|
|
|
|
|
|
ε y |
|
|
|
|
||||
|
|
|
|
|
|
|
ε |
|
|
|
|
|
|||
|
ε= {ε }= |
γ |
|
z |
|
|
, |
||||||||
|
|
|
|
|
|
|
|
|
xy |
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
γ yz |
|
|
|
||||
|
|
|
|
|
|
|
γ zx |
|
|
|
|
||||
Stress-strain relation: |
|||||||||||||||
|
σ |
|
|
|
|
|
|
|
|
|
|
|
1 − v |
||
|
x |
|
|
|
|
|
|
|
|
|
|
|
v |
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
σ y |
|
|
|
|
|
|
|
|
|
|
v |
|||
|
σz |
|
= |
|
E |
|
|
|
|
|
|
0 |
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
(1 |
|
|
|
|
|
|
|
|||||
|
τxy |
|
+ v)(1 − 2v) |
0 |
|||||||||||
|
τ |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
yz |
|
|
|
|
|
|
|
|
|
|
|
||
|
τzx |
|
|
|
|
|
|
|
|
|
|
|
0 |
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
or [σij ]
or [εij ]
v |
v |
|
0 |
1 − v |
v |
|
0 |
v |
1 − v |
1 |
0 |
0 |
0 |
− 2v |
|
|
2 |
||
|
|
|
|
0 |
0 |
|
0 |
0 |
0 |
|
0 |
0
0
0
0
1 − 2v
2
0
(1)
(2)
0 |
|
ε |
|
|
0 |
|
|
||
|
|
εx |
||
0 |
|
|
|
y |
0 |
|
|
ε |
|
|
|
|
z |
|
0 |
|
γxy |
||
|
γ |
|
||
|
|
|
|
yz |
1 − 2v |
γzx |
|||
2 |
|
|
|
|
or |
σ = Eε |
(3) |
© 1997-2002 Yijun Liu, University of Cincinnati |
139 |
Lecture Notes: Introduction to Finite Element Method |
Chapter 6. Solid Elements for 3-D Problems |
|||
Displacement: |
|
|
|
|
u(x, |
y, z) |
u |
|
|
|
|
1 |
|
(4) |
u = v(x, y, z) |
= u2 |
|
||
|
|
|
|
|
w(x, y, z) |
u3 |
|
|
|
Strain-Displacement Relation:
ε |
|
|
= |
∂u , |
|
ε |
|
= |
|
∂v |
|
, ε |
|
= |
|
∂w |
, |
|
|
|
|
|
||||||||
x |
|
|
y |
|
∂y |
z |
|
∂z |
|
|
|
|
|
|||||||||||||||||
|
|
|
|
∂x |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||
γ |
|
|
= |
∂v |
+ |
∂u |
, |
γ |
|
|
|
= |
∂w |
+ |
∂v |
, γ |
|
= |
∂u + |
∂w |
(5) |
|||||||||
xy |
∂x |
∂y |
yz |
∂y |
∂z |
xz |
∂x |
|||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
∂z |
|
||||||||||||||||
or |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1 |
∂ui |
|
|
∂u j |
|
|
|
|
|
(i, j =1,2,3) |
|
|
||||||||||||
εij |
|
= |
|
2 |
|
|
|
+ |
|
|
|
|
|
, |
|
|
|
|
|
|
||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||
|
|
|
|
|
|
∂x j |
|
|
∂xi |
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||
or simply,
εij |
= |
1 |
(ui, j +u j,i ) |
(tensor notation) |
|
|
2 |
|
|
© 1997-2002 Yijun Liu, University of Cincinnati |
140 |
Lecture Notes: Introduction to Finite Element Method |
Chapter 6. Solid Elements for 3-D Problems |
Equilibrium Equations: |
|
|
||||||||||||||||||
|
|
∂σx |
+ |
|
|
∂τxy |
|
+ |
|
|
∂τxz |
|
+ f x = 0 , |
|||||||
|
|
∂x |
|
|
∂y |
|
∂z |
|
||||||||||||
|
|
|
|
|
|
|
|
|
|
|
||||||||||
|
∂τ yx |
|
+ |
|
∂σ y |
+ |
|
∂τ yz |
|
+ f y = 0 , |
||||||||||
|
|
∂x |
|
|
∂y |
|
|
∂z |
|
|||||||||||
|
|
|
|
|
|
|
|
|
|
|
||||||||||
|
∂τzx |
|
+ |
∂τzy |
|
+ |
∂σz |
|
|
+ fz = 0 , |
||||||||||
|
|
|
|
|
|
∂y |
|
|
||||||||||||
|
|
∂x |
|
|
|
|
|
|
|
∂z |
|
|
||||||||
or |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
σij, j + fi |
= 0 |
|||||||||
Boundary Conditions (BC’s): |
||||||||||||||||||||
ui = ui , |
|
|
|
on |
|
|
Γu ( specified displacement ) |
|||||||||||||
|
|
|
|
|
on |
|
Γσ |
( specified traction ) |
||||||||||||
ti = ti , |
|
|
|
|||||||||||||||||
(traction ti =σij |
|
n j |
) |
|
|
|
||||||||||||||
p
n
Γσ
Γ( =Γu +Γσ )
Γu
Stress Analysis:
(6)
(7)
Solving equations in (3), (5) and (6) under the BC’s in (7) provides the stress, strain and displacement fields (15 equations for 15 unknowns for 3-D problems). Analytical solutions are difficult to find!
© 1997-2002 Yijun Liu, University of Cincinnati |
141 |
