- •Copyright Notice
- •Table of Contents
- •FEM in Structural Analysis (The Procedure)
- •Example:
- •Chapter 1. Introduction
- •I. Basic Concepts
- •Examples:
- •Why Finite Element Method?
- •Applications of FEM in Engineering
- •Examples:
- •A Brief History of the FEM
- •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
- •Type of Refinements:
- •Examples:
- •VI. Convergence of FE Solutions
- •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 7. Structural Vibration and Dynamics |
II. Free Vibration
Study of the dynamic characteristics of a structure:
•natural frequencies
•normal modes (shapes)
Let f(t) = 0 and C = 0 (ignore damping) in the dynamic equation (8) and obtain
Mu&&+Ku = 0 |
(12) |
Assume that displacements vary harmonically with time, that is,
u (t ) = u sin( ωt ), u& (t ) = ω u cos( ωt ),
u&&(t ) = −ω 2 u sin( ωt ),
where u is the vector of nodal displacement amplitudes. Eq. (12) yields,
[K − ω 2 M ]u = 0 |
(13) |
This is a generalized eigenvalue problem (EVP).
Solutions?
© 1997-2003 Yijun Liu, University of Cincinnati |
163 |
Lecture Notes: Introduction to Finite Element Method |
Chapter 7. Structural Vibration and Dynamics |
Trivial solution: u = 0 for any values of ω (not interesting).
Nontrivial solutions: u ≠ 0 only if
K −ω2M |
|
= 0 |
(14) |
|
This is an n-th order polynomial of ω2, from which we can find n solutions (roots) or eigenvalues ωi.
ωi (i = 1, 2, …, n) are the natural frequencies (or characteristic frequencies) of the structure.
ω1 (the smallest one) is called the fundamental frequency.
For each ωi , Eq. (13) gives one solution (or eigen) vector
[K − ω i 2 M ]u i = 0 .
u i (i=1,2,…,n) are the normal modes (or natural modes, mode shapes, etc.).
Properties of Normal Modes
u iT K u j |
= 0 , |
|
u iT M u j = 0 , |
for i ≠ j, |
(15) |
if ωi ≠ωj . That is, modes are orthogonal (or independent) to each other with respect to K and M matrices.
© 1997-2003 Yijun Liu, University of Cincinnati |
164 |
Lecture Notes: Introduction to Finite Element Method |
Chapter 7. Structural Vibration and Dynamics |
Normalize the modes:
u iT M u i = 1, |
|
u iT K u i = ω i2 . |
(16) |
Note:
•Magnitudes of displacements (modes) or stresses in normal mode analysis have no physical meaning.
•For normal mode analysis, no support of the structure is necessary.
ωi = 0 there are rigid body motions of the whole or a part of the structure.
apply this to check the FEA model (check for mechanism or free elements in the models).
•Lower modes are more accurate than higher modes in the FE calculations (less spatial variations in the lower modesfewer elements/wave length are needed).
© 1997-2003 Yijun Liu, University of Cincinnati |
165 |
Lecture Notes: Introduction to Finite Element Method |
Chapter 7. Structural Vibration and Dynamics |
Example:
y |
v2 |
|
|
|
|
ρ, A, EI |
|
θ2 |
1 |
2 |
x |
L |
|
|
|
|
|
v |
|
|
|
0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
[K −ω2M] |
2 |
|
= |
|
, |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||
|
|
|
θ2 |
|
|
0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
K = |
EI |
12 |
|
−6L |
|
|
M = |
|
ρAL |
156 |
−22L |
|
|||||||||||||||||||
|
|
|
|
|
|
|
|
2 |
, |
|
|
|
|
|
|
|
|
|
|
|
|
2 |
. |
|
|||||||
3 |
|
|
|
|
|
|
|
|
|
|
420 |
|
|
|
|
||||||||||||||||
|
L |
−6L |
|
|
4L |
|
|
|
|
|
|
|
−22L |
4L |
|
|
|||||||||||||||
EVP: |
|
12−156λ |
|
−6L +22Lλ |
|
=0, |
|
|
|
|
|
|
|
||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||
|
|
|
−6L + |
22Lλ |
2 |
2 |
λ |
|
|
|
|
|
|
|
|
|
|
||||||||||||||
|
|
|
4L |
−4L |
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||
in which λ = ω 2 ρAL4 / 420 EI . |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||
Solving the EVP, we obtain, |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
EI |
|
|
12 |
v2 |
|
|
1 |
|
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
, |
|
|
|
= 1.38 |
|
|
, |
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
4 |
|
|
|
|
|
|
|||||||||||
|
|
|
|
|
|
|
|
ω1 = 3.533 |
|
|
|
|
|
|
|
||||||||||||||||
|
|
#2 |
|
|
|
|
|
|
|
|
|
|
ρAL |
|
|
|
|
|
|
θ2 |
|
|
|
|
|
|
|||||
#3 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
12 |
|
|
|
1 |
|
|
L |
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1 |
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
EI |
|
v2 |
|
|
|
||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||
|
|
|
|
|
|
|
|
ω2 |
= |
|
|
|
|
|
|
, |
|
|
|
= |
|
|
|
. |
|||||||
#1 |
|
|
|
|
|
|
34.81 |
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
ρAL4 |
|
|
|
|
|
θ |
2 |
7.62 |
|
|
|||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
2 |
|
|
|
L |
||||||
Exact solutions: |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
EI |
|
|
12 |
|
|
|
|
|
|
|
|
EI |
|
|
12 |
|
|
|
|
|
|
|||||
ω1 = 3.516 |
|
|
|
|
|
|
|
, |
|
ω2 = 22.03 |
|
|
|
|
|
. |
|
|
|
|
|
|
|||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||
|
|
|
|
ρAL4 |
|
|
|
|
|
|
|
ρAL4 |
|
|
|
|
|
|
|
||||||||||||
We can see that mode 1 is calculated much more accurately than mode 2, with one beam element.
© 1997-2003 Yijun Liu, University of Cincinnati |
166 |
