Objective
To verify the behavior of the MTC4 thick plate bending element in a skew shape
Problem Description
A skewed, simply supported plate is loaded with a uniform pressure load of 1 psi.
Material properties: E = 1e5 psi, ν = 0.3
Geometric properties: L = 100 in, h = 1 in
Finite Element Model
16, 64, 256, 1024 shell elements
Model type: 2D Plate Bending (MITC4 thick plate formulation)
Results
The displacement at the plate center (C) is given by [Ref 1].
Unit: displacement - in
Displacement Dz |
ENERCALC 3D |
[Ref 1] |
4 x 4 mesh |
3.9182 |
3.9182 |
8 x 8 mesh |
3.8991 |
3.8991 |
16 x 16 mesh |
4.1875 |
4.1875 |
32 x 32 mesh |
4.4098 |
4.4098 |
Comments
The displacements given by ENERCALC 3D are identical to the referenced values. The correct theoretical displacement is given as 4.640 in.
Reference
[1]. Sa, Jorge, Valente and Areias “Development of shear locking-free shell elements using an enhanced assumed strain formulation”, International Journal of Numerical Methods in Engineering, 2002; 53: 1721-1750