Nonlinear Finite Element Analysis: Mema/Cses 648 Spring 2002
Nonlinear Finite Element Analysis: Mema/Cses 648 Spring 2002
Nonlinear Finite Element Analysis: Mema/Cses 648 Spring 2002
You are required to develop a geometrically nonlinear finite element analysis computer program
based on Newton-Raphson iterative technique to analyze static bending of linear elastic plates
using the first-order shear deformation plate theory (see Chapter 6 of the Lecture Notes). Use
the following definitions of plate extensional (Aij ) and bending (Dij ) stiffnesses in terms of the
engineering constants:
E1 h
A11 = , A12 = ν21 A11 , A22 = A11 (E2 /E1 ) , ν21 E1 = ν12 E2
1 − ν12 ν21
E1 h3
D11 = , D12 = ν21 D11 , D22 = D11 (E2 /E1 )
12(1 − ν12 ν21 )
G12 h3
A66 = G12 h, A44 = Ks G23 h, A55 = Ks G13 h, D66 = , Ks = 5/6
12
where Ks is the shear correction factor and h is the total thickness of the plate.
Use a simply-supported square plate with a = 10 in., h = 1 in., E1 = 7.8×106 psi, E2 = 2.6×106
psi, ν12 =0.25, G12 = G13 = G23 = 1.3 × 106 psi, and subjected to uniformly distributed transverse
load (i.e., q = q0 =constant) to validate your program. Use 8 × 8 and 4 × 4 uniform meshes of linear
and nine-node quadratic rectangular elements in a quadrant of the plate, respectively. Use load
steps such that the load parameter P ≡ q0 a4 /E2 h4 is equal to 5, and use twelve load steps. Plot
(a) load vs. center deflection (w0 /h versus P ), and (b) load vs. maximum stress (σxx and σxy ).
Use a converegnce tolerance of ε = 10−2 and a reasonable value of IT M AX.
You must submit: (1) a hard copy of the parts of the computer program that you have de-
veloped, (2) input for the problems, and (3) plots/tables of the numerical results obtained for the
problems.
ANNOUNCEMENT
Test 1 is scheduled for March 7th. You may pick up the test on Thursday (March 7th)
and return by noon on March 8th. The test will require you to solve a given set of nonlinear
problems using the progarms you have developed. You must also include all formulative steps.