Product: ABAQUS/Standard
This example illustrates the use of the *ORIENTATION option (Orientations, Section 2.2.5 of the ABAQUS Analysis User's Manual) in the analysis of multilayered, laminated, composite shells. The problem considered is the linear analysis of a flat plate made from two layers oriented at 45°, subjected to a uniform pressure loading. The example verifies simple laminated composite plate analysis. The ABAQUS results are compared with the analytical solution given in Spilker et al. (1976). The cross-section is not balanced, so the response includes membrane-bending coupling. Composite failure measures are defined for the plane stress orthotropic material.
The structure is a two-layer, composite, orthotropic, square plate that is simply supported on its edges. The layers are oriented at 45° with respect to the plate edges. Figure 1.1.21 shows the loading and the plate dimensions. Each layer has the following material properties:
These properties are specified directly in the *ELASTIC, TYPE=LAMINA option (Linear elastic behavior, Section 17.2.1 of the ABAQUS Analysis User's Manual), which is provided for defining linear elastic behavior for a lamina under plane stress conditions. More general orthotropic properties (for solid continuum elements) can be specified with the *ELASTIC, TYPE=ORTHO option.In this example the plate is considered to be at an arbitrary angle to the global axis system to make use of the *ORIENTATION option for illustration purposes. The plate is shown in Figure 1.1.22.
The boundary conditions require that displacements that are transverse and normal to the shell edges are fixed, but motions that are parallel to the edges are permitted. The *TRANSFORM option (Transformed coordinate systems, Section 2.1.5 of the ABAQUS Analysis User's Manual) has been used to define a convenient set of local displacement degrees of freedom so that the boundary conditions and the output of nodal variables can be interpreted more easily.
The *ORIENTATION option is used to define the direction of the layers. The rotation of the material axes of the layers with respect to the standard directions used by ABAQUS for stress and strain components in shells is defined on data lines in four of the models used and, again for illustration purposes, by means of user subroutine ORIENT in four other models. The section is not balanced since it has only two layers in different orientations, which results in membrane-bending coupling. The motion does not exhibit symmetry for the same reason, and the entire shell must be modeled.
An alternative means of defining the layer orientation is to use the *ORIENTATION option to define the orientation of the section and then to define the in-plane angle of rotation relative to the section orientation directly with the layer data following the *SHELL SECTION or *SHELL GENERAL SECTION option. In this case the section force and section strain are calculated in the section orientation directions (rather than the default shell directions).
Three types of models are used. One is an 8 × 8 mesh of S9R5 elements, which are shell elements that allow transverse shear along lines in the element. However, the analytical solution of Spilker et al. uses thin shell theory, which neglects transverse shear effects. We have, therefore, introduced an artificially high transverse shear stiffness in this model by using the *TRANSVERSE SHEAR STIFFNESS option.
The second type of model is a 16 × 16 mesh of triangular shells; models for both S3R and SC6R elements are provided. These elements are general-purpose shell elements that allow transverse shear deformation. An artificially high transverse shear stiffness is introduced by using the *TRANSVERSE SHEAR STIFFNESS option. No mesh convergence studies have been performed, but finer meshes should improve accuracy since these elements use a constant bending strain approximation.
The third type of model is made up of STRI65 shell elements, which are also based on the discrete Kirchhoff theory. An 8 × 8 mesh is used.
To demonstrate the use of composite failure measures (Plane stress orthotropic failure measures, Section 17.2.3 of the ABAQUS Analysis User's Manual), limit stresses are defined with the *FAIL STRESS option. The stress-based failure criteria are defined as follows:
Printed failure indices are requested for maximum stress theory (MSTRS) and Tsai-Hill theory (TSAIH). All failure measures are written to the results file (CFAILURE).Table 1.1.21 summarizes the results by comparing displacement and moment values to the analytical solution. It is clear by the results presented in the table that all models give good results, with the second-order models providing higher accuracy than the first-order S3R model, as would be expected.
Figure 1.1.23 shows the failure surface for Tsai-Hill theory (i.e., those stress values that, for a given , yield a failure index 1.0), along with the stress state at each section point in the center of the plate. Only section point 6 has a stress state outside the failure surface (1.0).
S3R element model with the orientation for the material defined with *ORIENTATION.
S3R element model with the orientation for the material defined in user subroutine ORIENT.
User subroutine ORIENT used in anisoplate_s3r_usr_orient.inp.
SC6R element model with the orientation for the material defined with *ORIENTATION.
SC6R element model with the orientation for the material defined in user subroutine ORIENT.
SC6R model with the orientation for the shell section defined with *ORIENTATION and the orientation for the material defined by an angle on the data lines for *SHELL GENERAL SECTION.
User subroutine ORIENT used in anisoplate_sc6r_usr_orient.inp.
S9R5 model with the orientation for the material defined with *ORIENTATION.
S9R5 model with the orientation for the material defined in user subroutine ORIENT.
User subroutine ORIENT used in anisoplate_s9r5_usr_orient.inp.
S9R5 model with the orientation for the shell section defined with *ORIENTATION and the orientation for the material defined by an angle on the data lines for *SHELL SECTION.
S9R5 model with the orientation for the shell section defined with *ORIENTATION and the orientation for the material defined by an angle on the data lines for *SHELL GENERAL SECTION.
STRI65 element model with the orientation for the material defined with *ORIENTATION.
STRI65 element model with the orientation for the material defined in user subroutine ORIENT.
User subroutine ORIENT used in anisoplate_stri65_usr_orient.inp.
Spilker, R. L., S. Verbiese, O. Orringer, S. E. French, E. A. Witmer, and A. Harris, Use of the Hybrid-Stress Finite-Element Model for the Static and Dynamic Analysis of Multilayer Composite Plates and Shells, Report for the Army Materials and Mechanics Research Center, Watertown, MA, 1976.