1.6.1 Plate and Shell Element: Geometry Modeling Considerations
1.6.1 Plate and Shell Element: Geometry Modeling Considerations
1.6.1 Plate and Shell Element: Geometry Modeling Considerations
Page 1 of 12
2. While assigning nodes to an element in the input data, it is essential that the nodes be specified either clockwise or counter clockwise (see the following figure). For better efficiency, similar elements should be numbered sequentially.
Figure 1-15: Examples of correct and incorrect numbering sequences
mk:@MSITStore:C:\SProV8i\STAAD\Help\Technical_Reference_2007.chm::/General_Des... 9/4/2013
Page 2 of 12
3. Element aspect ratio should not be excessive. They should be on the order of 1:1, and preferably less than 4:1. 4. Individual elements should not be distorted. Angles between two adjacent element sides should not be much larger than 90 and never larger than 180.
Figure 1-16: Some examples of good and bad elements in terms of the angles
mk:@MSITStore:C:\SProV8i\STAAD\Help\Technical_Reference_2007.chm::/General_Des... 9/4/2013
Page 3 of 12
7. Temperature load due to difference in temperature between top and bottom surfaces of the element.
Theoretical Basis
The STAAD plate finite element is based on hybrid finite element formulations. An incomplete quadratic stress distribution is assumed. For plane stress action, the assumed stress distribution is as follows.
Figure 1-17: Assumed stress distribution
a1 through a10 = constants of stress polynomials The following quadratic stress distribution is assumed for plate bending action:
Figure 1-18: Quadratic stress distribution assumed for bending
mk:@MSITStore:C:\SProV8i\STAAD\Help\Technical_Reference_2007.chm::/General_Des... 9/4/2013
Page 4 of 12
a1 through a13 = constants of stress polynomials The distinguishing features of this finite element are: 1. Displacement compatibility between the plane stress component of one element and the plate bending component of an adjacent element which is at an angle to the first (see the following figure) is achieved by the elements. This compatibility requirement is usually ignored in most flat shell/plate elements.
Figure 1-19: Adjacent elements at some angle
2. The out of plane rotational stiffness from the plane stress portion of each element is usefully incorporated and not treated as a dummy as is usually done in most commonly available commercial software. 3. Despite the incorporation of the rotational stiffness mentioned previously, the elements satisfy the patch test absolutely. 4. These elements are available as triangles and quadrilaterals, with corner nodes only, with each node having six degrees of freedom. 5. These elements are the simplest forms of flat shell/plate elements possible with corner nodes only and six degrees of freedom per node. Yet solutions to sample problems converge rapidly to accurate answers even with a large mesh size. 6. These elements may be connected to plane/space frame members with full displacement compatibility. No additional restraints/releases are required. 7. Out of plane shear strain energy is incorporated in the formulation of the plate bending component. As a result, the elements respond to Poisson boundary conditions which are considered to be more accurate than the customary Kirchoff boundary conditions. 8. The plate bending portion can handle thick and thin plates, thus extending the usefulness of the plate elements into a multiplicity of problems. In addition, the thickness of the plate is taken into consideration in calculating the out of plane
mk:@MSITStore:C:\SProV8i\STAAD\Help\Technical_Reference_2007.chm::/General_Des... 9/4/2013
Page 5 of 12
shear. 9. The plane stress triangle behaves almost on par with the well known linear stress triangle. The triangles of most similar flat shell elements incorporate the constant stress triangle which has very slow rates of convergence. Thus the triangular shell element is very useful in problems with double curvature where the quadrilateral element may not be suitable. 10. Stress retrieval at nodes and at any point within the element.
mk:@MSITStore:C:\SProV8i\STAAD\Help\Technical_Reference_2007.chm::/General_Des... 9/4/2013
Page 6 of 12
Title
Description
SQX, SQY Shear stresses (Force/ unit length/ thickness) SX, SY SXY Membrane stresses (Force/unit length/ thickness) Inplane Shear Stress (Force/unit length/ thickness) Moments per unit width (Force x Length/length) (For Mx, the unit width is a unit distance parallel to the local Y axis. For My, the unit width is a unit distance parallel to the local X axis. Mx and My cause bending, while Mxy causes the element to twist out-of-plane.) SMAX, SMIN Principal stresses in the plane of the element (Force/unit area). The 3rd principal stress is 0.0 Maximum 2D shear stress in the plane of the element (Force/unit area) 3D Von Mises stress at the top and bottom surfaces, where: VM = 0.707[(SMAX - SMIN)2 + SMAX2 + SMIN2]1/2 TRESCAT, Tresca stress, where TRESCA = MAX[ |(SmaxTRESCAB Smin)| , |(Smax)| , |(Smin)| ] Notes 1. All element stress output is in the local coordinate system. The direction and sense of the element stresses are explained in the following section. 2. To obtain element stresses at a specified point within the element, you must provide the location (local X, local Y) in the coordinate system for the element. The origin of the local coordinate system coincides with the center of the element. 3. The 2 nonzero Principal stresses at the surface (SMAX & SMIN), the maximum 2D shear stress (TMAX), the 2D orientation of the principal plane (ANGLE), the 3D Von Mises stress (VONT & VONB), and the 3D Tresca stress (TRESCAT &
TMAX
VONT, VONB
mk:@MSITStore:C:\SProV8i\STAAD\Help\Technical_Reference_2007.chm::/General_Des... 9/4/2013
Page 7 of 12
TRESCAB) are also printed for the top and bottom surfaces of the elements. The top and the bottom surfaces are determined on the basis of the direction of the local zaxis. 4. The third principal stress is assumed to be zero at the surfaces for use in Von Mises and Tresca stress calculations. However, the TMAX and ANGLE are based only on the 2D inplane stresses (SMAX & SMIN) at the surface. The 3D maximum shear stress at the surface is not calculated but would be equal to the 3D Tresca stress divided by 2.0.
mk:@MSITStore:C:\SProV8i\STAAD\Help\Technical_Reference_2007.chm::/General_Des... 9/4/2013
Page 8 of 12
Mx is the Bending Moment on the local x face and the local x-face is the face perpendicular to the local x-axis. My is the Bending Moment on the local y face and the local y-face is the face perpendicular to the local y-axis.
Figure 1-23: Stress caused by Mx
mk:@MSITStore:C:\SProV8i\STAAD\Help\Technical_Reference_2007.chm::/General_Des... 9/4/2013
Page 9 of 12
mk:@MSITStore:C:\SProV8i\STAAD\Help\Technical_Reference_2007.chm::/General_Des... 9/4/2013
Page 10 of 12
mk:@MSITStore:C:\SProV8i\STAAD\Help\Technical_Reference_2007.chm::/General_Des... 9/4/2013
Page 11 of 12
Members, plate elements, solid elements and surface elements can all be part of a single STAAD model. The MEMBER INCIDENCES input must precede the INCIDENCE input for plates, solids or surfaces. All INCIDENCEs must precede other input such as properties, constants, releases, loads, etc. The selfweight of the finite elements is converted to joint loads at the connected nodes and is not used as an element pressure load.
mk:@MSITStore:C:\SProV8i\STAAD\Help\Technical_Reference_2007.chm::/General_Des... 9/4/2013
Page 12 of 12
Reference Materials
5.11 Joint Coordinates Specification 5.13 Elements and Surfaces 5.13.1 Plate and Shell Element Incidence Specification 5.13.2 Solid Element Incidences Specification 5.14 Plate Element Mesh Generation 5.14.1 Element Mesh Generation 5.21 Element/Surface Property Specification 5.21.1 Element Property Specification 5.24 Element Plane Stress and Ignore Inplane Rotation Specification 5.32.3 Element Load Specifications 5.32.3.1 Element Load Specification - Plates 5.22.2 Element Release Specification
mk:@MSITStore:C:\SProV8i\STAAD\Help\Technical_Reference_2007.chm::/General_Des... 9/4/2013