Theory Background for the Membrane Interface
A 3D membrane is similar to a shell but it has only translational degrees of freedom and the results are constant in the thickness direction.
The thickness of the membrane is d, which can vary over the element. The displacements are interpolated by Lagrange or Serendipity shape functions.
A 2D axisymmetric membrane is similar to the 3D membrane and it has a nonzero circumferential strain in the out-of-plane direction.
Local Coordinate Systems
Boundaries
Many quantities for a membrane can best be interpreted in a local coordinate system aligned to the membrane surface. Material data, initial stresses-strains, and constitutive laws are always represented in the local coordinate system.
This local membrane surface coordinate system is defined by the boundary coordinate system (t1, t2, n).
The quantities like stresses and strains are also available as results in the global coordinate system after a transformation from a local (boundary) system.
Local Edge System
Many features, such as an edge load, allow input in an edge local coordinate system. The orthogonal local edge coordinate system directions xl, yl, and zl are defined so that:
The first direction (xl) is along the edge. This direction can be visualized by selecting the Show edge directions arrows check box in the View settings.
The third direction (zl) is the same as the membrane normal direction of the adjacent boundary.
The second direction (yl) is in the plane of the shell and orthogonal to the edge. It is formed by the cross product of zl and xl; yl = zl × xl.
The local edge system can be visualized by plotting the components of the local edge transformation matrix with an Arrow Line plot. The matrix components are defined per feature. For instance, the variable name for the xx-component is <interface>.<feature_tag>.TleXX.
When an edge is shared between two or more boundaries, the directions may not always be unique. It is then possible to use the control Face Defining the Local Orientations to select from which boundary the normal direction zl should be picked. The default is Use face with lowest number.
If the geometry selection contains several edges, the only available option is Use face with lowest number, since the list of adjacent boundaries would then be different for each edge. For each edge in the selection, the face with the lowest number attached to that edge is then used for the definition of the normal orientation.
Strain-Displacement Relation
The kinematic relations of the membrane element are first expressed along the global coordinate axes. The strains are then transformed to the element local direction. Since the membrane is defined only on a boundary, derivatives in all spatial directions are not directly available. This makes the derivation of the strain tensor somewhat different from what is used in solid mechanics.
The deformation gradient F is in general defined as the gradient of the current coordinates with respect to the original coordinates:
In the Membrane interface, a tangential deformation gradient is computed as
Here is a displacement gradient computed using the tangential derivative operator, and N is the normal vector to the undeformed membrane. FT now contains information about the stretching in the plane of the membrane.
Since the tangential deformation gradient does not contain any information about the transversal stretch λn, it must be augmented by the normal deformation gradient FN to define the full deformation gradient tensor. It is given by
where n is the normal vector to the deformed membrane. For anisotropic materials, the shear deformation gradient FS is also necessary to define the full deformation gradient tensor. It is given by
where t1 and t2 are the tangent vectors to the deformed membrane. The full deformation gradient tensor F is the sum of tangential, shear and normal deformation gradient tensors
Note that FS is only nonzero for anisotropic materials, otherwise FS = 0.
The Right Cauchy-Green tensor C is generally defined as
From C, the Green-Lagrange strains are computed using the standard expression
The local tangential strains in the membrane are calculated by transformation of this strain tensor into the local coordinate system.
The Jacobian J is the ratio between the current volume and the initial volume. In full 3D it is defined as
In the membrane, only the C tensor is available, so instead the following expression is used:
The area scale factor is also computed as
In the case of geometrically linear analysis, a linearized version of the strain tensor is used.
Constitutive Relation and Weak Contributions
The constitutive relations for the membrane on the reference surface are similar to those used in the Solid Mechanics interface.
The thermal strains and initial stresses-strains (only for the in-plane directions of the membrane) are added in the constitutive relation in a similar manner as it is done in Solid Mechanics.
The weak expressions in the Membrane interface are similar to that of linear elastic continuum mechanics.
See also Analysis of Deformation in the documentation of the Solid Mechanics interface.
External Loads
Contributions to the virtual work from the external load are of the form
where the forces (F) can be distributed over a boundary or an edge or be concentrated in a point. In the special case of a follower load, defined by its pressure p, the force intensity is where n is the normal in the deformed configuration.
For a follower load, the change in midsurface area is taken into account, and integration of the load is done in the spatial frame.
Stress Calculations
The stresses are computed by applying the constitutive law to the computed strains.
The membrane does not support transverse and bending forces so the only section forces it support is the membrane force defined as:
where is the local stress tensor and contains only in-plane stress components.