Users:General FEM Analysis/Elements Reference/Membrane1

Revision as of 08:01, 2 September 2010

General Description

Element Type

• This membrane element is a plane stress element with 3 degrees of freedom per node
• The membrane element neglects the bending stiffness
• The ratio of the thickness and the lengths in the other two directions is much smaller than one ( t/Lx << 1 and t/Ly << 1).
• The thickness is constant over the element

Degrees of Freedom

For the membrane element use the 3 translatoric degrees of freedom. DISP_X, DISP_Y, DISP_Z

Input Parameters

Parameter Description

Example of a Complete Input Block

EL-PROP 1 : MEMBRANE1
MAT= EL-MAT 1        THICKNESS=1.0
PRESTRESS      SIG11=1.0     SIG22=1.0      SIG12=0.0
A_X=1.0     A_Y=0.0    A_Z=0.0
B_X=0.0    B_Y=1.0    B_Z=0.0  
LAGRANGE=UPDATED

Element Loading

Pressure

Dead Load

Snow Load

Theory

The theory and finite element formulation is described in detail in the work of Dieringer ^[1] and Linhard ^[2]

For the correct use of the membrane element and the interpretation of related results, the following aspects should be considered:

Material parameter

With the parameter MAT the material, that would be used in the calculation, for the membrane element is defined. A very popular mistake is that a material number would be used, that isn't defined already. Please ensure that the material that you would like to use at this point are already exits. The following materials are testet for the membrane element:

• linear elastic isotropic
• linear elastic orthotropic (Münsch-Reinhardt)
• multilinear elastic isotropic
• elastoplastic isotropic
• material on the basis of response functions

Thickness parameter

The parameter THICKNESS defines the thickness of the membrane. The thickness should be constant over the element.

Prestress directions on the surface

The parameters A_X, A_Y, A_Z and B_X, B_Y, B_Z are used to define the prestress directions on the surface. In this approach the principle directions of the prestress are defined in a plane area (see figure below). The definition of the area is given by the two vectors f₁ and f₂. The normal vector of the area can be calculated with the cross product of the in plane vectors f₃=f₁ x f₂. Afterwards the line of intersection T₁ of the area which is given by f₁ and f₃ and the curved surface can be calculated. In this approach T₁ is interpreted as the first principle direction of the prestress on the curved surface. With the assumption that T₃ is equal to the surface normal vector G₃ (not uniformed), the second direction of the prestress is calculated as T₂=T₁ x T₃. W.r.t. the parameters for the input file, only the plane area with the vectors f₁ and f₂ have to be defined. Referring to the depict approach the vector A defines the vector f₁ and the vector B defines the vector f₂.

Prestress state &sigma

SIG11, SIG22, SIG12 describes the prestress of the membrane element. The element is based on the plane stress assuption. Due to that only normal and shear stresses in the midplane have to be defined. SIGG11 is the stress acting in T₁, SIG22 acting in T₂ and SIG12 is the in plane shear, whereas SIG12=SIG21 (see figure below).

Lagrange type

With the LARANGE parameter it is possible to switch between form finding and statical/dynamical analysis. For the value UPDATED the element is for form finding and for the value TOTAL the element is for statical/dynamical analysis. It is important that the LANGRANGE parameter match to the type of analysis.

References

1. ↑ Dieringer, F.: Implementierung eines geometrisch nichtlinearen Membranelements in einer objektorientierten Programmierumgebung, Masterthesis, Chair of structural analysis, Technische Universität München, 2009
2. ↑ Linhard, J.: Numerisch-mechanische Betrachtung des Entwurfsprozesses von Membrantragwerken, PhD-Thesis, Chair of structural analysis, Technische Universität München, 2009