Plate elements are three- or four-node elements formulated in three-dimensional space. These elements are used to model and analyze objects such as pressure vessels, or structures such as automobile body parts.
The out-of-plane rotational DOF is not considered for plate elements. You can apply the other rotational DOFs and all the translational DOFs as needed.
Nodal forces, nodal moments (except when about an axis normal to the element face), pressures (normal to the element face), acceleration/gravity, centrifugal and thermal loads are supported.
Surface-based loads (pressure, surface force, and so on, but not constraints) and element properties (thickness, element normal coordinate, and so on) are applied to an entire plate element. Since these items are based on the surface number of the lines forming the element, and since each element could be composed of lines on four different surface numbers, how these items are applied depend on whether the mesh is created automatically (by either the mesher from a CAD model or the 2D mesh generation), or whether the mesh is created by hand. The surface number of the individual lines that form an element are combined as indicated in Table 1 to create a surface number for the whole element. Loads and properties are then applied to the entire element based on the element's surface number.
How Mesh Was Created |
Definition of Surface Number of Element |
Midplane Mesh from CAD Model |
All elements coincident with the collapsed surface of the CAD model are on the CAD's surface number regardless of the surface number of the lines. |
Plate/Shell Mesh from CAD Model |
All elements coincident with the surface of the CAD model are on the CAD's surface number regardless of the surface number of the lines. |
2D Mesh from Sketches |
All elements are assigned to surface number 1 regardless of the surface number of the lines. |
Hand-built Mesh and Modified Automatic Meshes |
The highest surface number of any line on the element determines the Surface Number of the entire element. This is the voting rule. |
Table 1: Definition of Plate Element Surface Number |
To enter the element parameters, select the Element Definition entry in the tree view for the plate element part or parts, right-click, and choose Edit Element Definition. Alternatively, select the part or parts in the display area or tree view, right-click, and choose Edit Element Data.
Material Model: Specify the material model for this part in the Material Model drop-down Menu. If the material properties in all directions are identical, select the Isotropic option. If the material properties vary along two orthogonal axes, select the Orthotropic option. (The orientation of the orthotropic axes is then defined using the Nodal Order Method option. See below.)
Element Formulation: Specify which type of element formulation is used for this part in the Element Formulation drop-down menu. The Veubeke option uses the theory by B. Fraeijs de Veubeke for plate formulation for displaced and equilibrium models. This option is recommended for plate elements that have little or no warpage. The Reduced Shear option uses the constant linear strain triangle (CLST) with reduced shear integration and Hsieh, Clough and Tocher (HCT) plate bending element theories. This option is recommended for plate elements that contain significant warpage. The Linear Strain option uses the CLST without reduced shear integration and HCT plate bending element theories. The Constant Strain option uses the constant strain triangle (CST) and HCT plate bending element theories.
Temperature Method: There are three options for performing a thermal stress analysis with plate elements. These are selected in the Temperature Method drop-down menu. If the Stress Free option is selected, the thermal strain (ε) is calculated as the product of the difference of the nodal temperatures (Tnode) applied to the model and the Stress Free Reference Temperature (Tref), and the thermal coefficient of expansion (α): ε = α(Tnode-Tref). The Stress Free Reference Temperature is entered in the appropriate field of the Element Definition dialog box. If the Mean option is selected, the thermal strain is calculated as the product of the Mean Temperature Difference (entered in the spreadsheet) and the thermal coefficient of expansion: ε = α(Mean Temperature Difference). If the Nodal dT option is selected, the thermal strain is calculated as the product of the difference of the nodal temperatures applied to the model and 0 degrees and the thermal coefficient of expansion: ε = α(Tnode-0). (Also see delta T thru thickness below.)
Twisting coefficient ratio: The undefined rotational degree of freedom (the direction perpendicular to the element) for a plate element is assigned an artificial stiffness to help stabilize the solution. The magnitude of the artificial stiffness equals the Twisting coefficient ratio times the smallest bending stiffness of the element.
The linear plate element is a combination of planar plate and membrane elements. The rotational degree of freedom perpendicular to the plate element is undefined on a local basis. When combined with other plate elements at an angle, the global rotational degree of freedom is defined. (Visualize this as the in-plane rotation in one element having a component in the out-of-plane direction for the adjacent element.) To avoid a singularity (unknown solution) in the solution of the global stiffness matrix, the twisting coefficient is used to create an artificial stiffness on a local basis. This local stiffness is added to the global stiffness matrix. If this artificial stiffness is too large, the solution behaves as if the model is partially tied down in the twisting direction.
Values for the twisting coefficient ratio that are too large may cause a significant artificial constraint, especially where plates meet at an angle. Values that are too small can increase the maximum/minimum stiffness ratio. A large maximum/minimum stiffness ratio may cause a warning and can make the matrix harder to solve, increasing the chance of an inaccurate solution. (The warning is output during the assembly of the stiffness matrix and before the solving operation. It may be followed by solution warnings which are a much more serious indicator of problems.)
The maximum/minimum stiffness ratio is not always independent of the units. If the maximum and minimum stiffnesses were due to tension, then the units of each (such as N/mm) are canceled. With plate elements, the maximum stiffness is often a tension (units of force/length) and the minimum stiffness is often the out-of-plane rotation (units like force*length/radian), so the maximum stiffness divided by the minimum stiffness does have units. The Twisting coefficient ratio may need to be adjusted depending on the units in use.
Properties: The majority of the Element Definition input is entered in a spreadsheet. The specifics of the input depend on the selection in the Properties drop-down menu and the Use mid-plane mesh thickness check box. The options are as follows:
Figure 1: Thickness of a Plate Element
![]() |
Figure 2: Determining the Element Normal The edge-on view of the plate element is shown. |
|
Figure 3: Local 1 and 2 axes for Plate Elements The dots along the side of the element are at the midpoint of the side. |
![]() |
|
Delta T thru thickness |
= (Ttop - Tbottom) / thickness = (100 - 80 °F) / (0.1 inch) = 200 °F/ inch |
Figure 4: Temperature Gradient Through a Plate Element |