Nonlinear Finite Element Method

Nonlinear Finite Element Method
13/12/2004
Nonlinear Finite Element Method

Lectures include discussion of the nonlinear finite element method. It is preferable to have completed Introduction to Nonlinear Finite Element Analysis available in summer session. If not, students are required to study on their own before participating this course. Reference:Toshiaki.,Kubo. Introduction: Tensor Analysis For Nonlinear Finite Element Method (Hisennkei Yugen Yoso no tameno Tensor Kaiseki no Kiso),Maruzen. Lecture references are available and downloadable at http://www.sml.k.utokyo.ac.jp/members/nabe/lecture2004 They should be posted on the website by the day before scheduled meeting, and each students are expected to come in with a copy of the reference. Lecture notes from previous year are available and downloadable, also at http://www.sml.k.u.tokyo.ac.jp/members/nabe/lecture2003 You may find the course title, Advanced Finite Element Method but the contents covered are the same I will cover this year. I will assign the exercises from this year, and expect the students to hand them in during the following lecture. They are not the requirements and they will not be graded, however it is important to actually practice calculate in deeper understanding the finite element method. For any questions, contact me at [email protected]
Nonlinear Finite Element Method Lecture Schedule

1. 10/ 4 Finite element analysis in boundary value problems and the differential equations 2. 10/18 Finite element analysis in linear elastic body 3. 10/25 Isoparametric solid element (program) 4. 11/ 1 Numerical solution and boundary condition processing for system of linear equations with exercises 5. 11/ 8 Basic program structure of the linear finite element method(program) 6. 11/15 Finite element formulation in geometric nonlinear problems(program) 7. 11/22 Static analysis techniquehyperelastic body and elastic-plastic material for nonlinear equations (program) 8. 11/29 Exercises for Lecture7 9. 12/ 6 Dynamic analysis technique and eigenvalue analysis in the nonlinear equations 10. 12/13 Constitutive element 11. 12/20 Numerical solution skyline methoditerative method for the system of linear equations 12. 1/17 ALE finite element fluid analysis 13. 1/24 ALE finite element fluid analysis
Dynamic Analysis Technique

Equation of motion at time t+t ,
M represents a mass matrix and C represents a damping matrix. Here, we state the direct time integration method, which integrate the equation as it is by direction of time. In linear problems, it is sometimes analyzed by superposition of natural mode, and it is called the mode superposition method. Given a solution at time t , a method of evaluating an accurate solution of t+t is called an implicit method.linear acceleration method and Newmark- method, etc.), while a method of making a prediction based on a solution of t is called an explicit method.central difference Normally, M and C stay constant regardless of time, yet Q changes. Linearize the inner vectors then adopt the Newton-Raphson method.
In this way, the above includes the unknown quantities thus, usually the three quantities are related under some supposition, and put them all together to make single variable in analyzation.
Linear Acceleration Method

Linear acceleration method is hardly used in actual analysis. Literally, it supposes the linear change in acceleration in time t and t+t .
Having
as a variable, and integrate the equation above then substitute
may yield,
Here we have,
Newmark - Method 1
Newmark- method is first introduced by Newmark in 1959 and since then it has contributed various benefits in applications to the area of dynamics. This technique is signified as unconditional stability of time integration scheme, thus the dynamic analysis can be conducted with a great time increment step . Moreover, introduction of numerical damping not only to the higher order modes but also to the lower order modes is possible. The equation of motion at time t+t is given by,
Newmark- method assumes the relation of deformation and velocity at time t+t from t as in the following equation.
Newmark - Method 2
If we set = method. If we set , = equations.
,= ,=
then the equation may coincide with the equation in linear acceleration in above equation, then and are given by the following
Apparently, we assume this relation to take a constant average velocity at t to t+t.
The combinations for =
and =
are idiomatically called the trapezoidal rule.
Newmark - Method 3
In transforming the relational expression of Newmark- method with a deformation at t+t in addition to the deformation, acceleration and velocity at t, we can express the acceleration and the velocity at t+t can be given by,
In linear problems we can use the linear stiffness matrix K and express the inner force
Executing equation of Newmark- method in linear system,
Newmark - Method 4
Executing equation of Newmark- method in linear system,
Linear analysis programming can be implemented if
is adopted in stead of
the usual stiffness matrix.

Under linear system, given
and =
, the solution is stable regardless of the value in
time stept. While , = and = are given, if the value t takes below the certain value, we may obtain a guarantee in solution stability, in other words, which implies a conditional stability. However, where > , the numerical damping effect may take place.
Newmark - Method 5
In the same way we handled stationary problems, iteration method based on the following equation can be implemented. on Newton-Raphson
Where signs in upper right side (k) indicate the k th iteration. In conducting iteration by Newmark method, if deformation is given unknown, take the difference of acceleration, deformation and velocity relational expression at the k th and k 1th iteration.
System of linear equation with unknown variable
Newmark - Method 6
Newmark- method implementation procedures are indicated in the following. 1. Initial value calculation (a) Obtain the mass matrix M and the damping matrix C. (b) Calculate the following invariables.
2. Calculations at each step (a)Calculate the significant stiffness matrix at time t+t.
(b) Calculate the significant load at time t+t.
(c) Evaluate the deformation of time t+t: t+
(d) For nonlinear problems, iteration can be calculated as,
iii. Calculate the significant load for the i 1 th iteration.
iv.Calculate the deformation increment for the i th iteration.
v. Calculate the deformation increment of single increment step.
vi. If convergence exists, then the deformation, velocity and acceleration at time t+t can be evaluated in the following calculation.

Nonlinear Finite Element Method

Uploaded by

Copyright:

Available Formats

Nonlinear Finite Element Method

Uploaded by

Document Information

Original Title

Copyright

Available Formats

Share this document

Share or Embed Document

Sharing Options

Did you find this document useful?

Is this content inappropriate?

Copyright:

Available Formats

Nonlinear Finite Element Method

Uploaded by

Copyright:

Available Formats

Nonlinear Finite Element Method

Nonlinear Finite Element Method

Nonlinear Finite Element Method Lecture Schedule

Dynamic Analysis Technique

Linear Acceleration Method

as a variable, and integrate the equation above then substitute

If we set = method. If we set , = equations.

Apparently, we assume this relation to take a constant average velocity at t to t+t.

The combinations for =

are idiomatically called the trapezoidal rule.

Executing equation of Newmark- method in linear system,

Linear analysis programming can be implemented if

the usual stiffness matrix.

, the solution is stable regardless of the value in

System of linear equation with unknown variable

(b) Calculate the significant load at time t+t.

(c) Evaluate the deformation of time t+t: t+

(d) For nonlinear problems, iteration can be calculated as,

iii. Calculate the significant load for the i 1 th iteration.

iv.Calculate the deformation increment for the i th iteration.

v. Calculate the deformation increment of single increment step.

You might also like