ProductsAbaqus/StandardAbaqus/Explicit Problem descriptionA schematic of the model is shown in Figure 1. The structure is a composite plate composed of orthotropic layers of equal thickness. It is simply supported at its ends and bounded along its edges to impose plane strain conditions in the y-direction. Each layer models a fiber/matrix composite with the following properties:
where L signifies the direction parallel to the fibers and T signifies the transverse direction. In Abaqus/Standard two methods are used to specify the lay-up definition for the conventional shell element model. First, a composite shell section is defined to specify the thickness, number of integration points, material name, and orientation of each layer. Second, a composite general shell section is defined to specify the thickness, material name, and angle of orientation relative to the section orientation (the default shell directions in this case) for each layer. In Abaqus/Explicit only the former method is used. The material properties are specified using the orthotropic elastic in plane stress definition. The orientation of the fibers in each layer is defined by an in-plane rotation angle measured relative to the local shell directions or relative to an orientation definition given for the general shell section. In addition to the methods outlined above, a third method of stacking continuum shell elements is used to specify the lay-up definition for a composite model. This method can be used effectively to study localized behavior, since continuum shell elements handle high aspect ratios between the in-plane dimension and the thickness dimension well. The lay-up definition for the continuum (solid) element model in Abaqus/Standard is specified using a composite solid section definition. The thickness, material name, and orientation definition are specified for each layer. A distributed load with a sinusoidal distribution in space, , is applied to the top of the composite plate. In Abaqus/Standard the load is applied using user subroutine DLOAD in a static linear analysis step. In addition, an Abaqus/Standard input file is included that demonstrates the use of the DCOUP3D element to apply this distributed load. In Abaqus/Explicit the load is applied instantaneously at time 0. Two composite plates are analyzed in this example. The first is a two-layer plate with the fibers oriented parallel and orthogonal to the x-axis in the bottom and top layer, respectively. In the second plate, which has three layers of equal thickness, the fibers in the outer layers are oriented parallel to the x-axis, while the fibers in the middle layer are orthogonal to the x-axis. The span-to-thickness ratio of the plates, , is varied from 4 to 30 in the Abaqus/Standard analysis; in Abaqus/Explicit this ratio is 4 throughout the analysis. A 1 × 10 mesh of second-order S8R shell elements is used to model the plates in Abaqus/Standard. A 2 × 10 mesh of first-order S4R shell elements is used to model the plates in Abaqus/Explicit. The S4R, S8R, and S8RT shell elements are well-suited for modeling thick composite shells since they account for transverse shear flexibility. Five integration points are specified through the thickness of each layer with the models that use the shell section. This provides sufficient data to describe the stress distributions through the thickness of each layer. For the models that use the general shell section, only three points are available for output. (Since the analysis is linear elastic, three points are sufficient to determine all fields through the thickness.) The plate with the lowest span-to-thickness ratio is also analyzed with Abaqus/Standard using a 1 × 10 mesh of second-order C3D20R composite solid elements. To illustrate the stacking capability of continuum shell elements, several meshes are provided for the two- and three-layer plates with a span-to-thickness ratio of 4. The two-layer plate is modeled with a 2 × 10 mesh of SC8R elements, each element representing a single layer of the 90/0 composite plate. One model of the three-layer plate uses a 1 × 10 mesh of SC8R elements using a single element through the thickness with a composite section definition. Another model of the three-layer plate uses a 3 × 10 mesh of SC8R elements, each element representing a single layer of the 0/90/0 composite plate. Additional models of the three-layer plate with 6, 12, and 24 elements through the thickness are provided. In these models each composite layer is modeled with 2, 4, and 8 elements through the thickness, respectively. Additional input files using SC8R elements are included to illustrate defining the stacking and thickness direction independent of the element nodal connectivity. Failure measuresThe plane stress orthotropic failure measures are defined in Plane stress orthotropic failure measures. To demonstrate their use, let the limit stresses and limit strains be given as follows:
The scaling factor for the Tsai-Wu coefficient is 0.0. These values are chosen such that failure occurs under the stress-based failure criteria for the given loading in the two-layer case with 4. Results and discussionThe results for each of the analyses are discussed in the following sections. Abaqus/Standard resultsFigure 2 shows the maximum z-displacement as a function of the span-to-thickness ratio of the two- and three-layer plates in a normalized form as As seen in the figure, the finite element displacements for both the two- and three-layer plates agree well with the prediction from elasticity theory for a wide range of s values. The CPT results are stiff at low values of s since shear flexibility is neglected. For 4, Figure 3 and Figure 4 show the transverse shear stress (TSHR13) and the axial stress (S11) distributions through the plate thickness for the two-layer plate normalized as and Figure 5 and Figure 6 show the corresponding results for the three-layer plate. It is seen that the shell element results are much closer to the predictions of CPT than to elasticity theory because of the assumption of linear stress variation through the thickness in the first-order shear flexible theory used for elements such as S8R and S4R. Figure 7 compares the elasticity solution of the transverse shear distribution for the three-layer plate to an approximate solution using the output variable SSAVG4. SSAVG4 is the average transverse shear stress in the local 1-direction. Since SSAVG4 is constant over an element, mesh refinement (in this case 24 continuum shell elements through the thickness) is typically required to capture the variation of shear stress through the thickness of the plate. The output variables CTSHR13 and CTSHR23 offer a more economical alternative to SSAVG4 and SSAVG5 for estimating shear stress in stacked continuum shells. Figure 8 and Figure 9 show very good agreement between the elasticity solution of the transverse shear distribution for the three- and two-layer plates to the solution using the output variable CTSHR13 for a 3 × 10 and 2 × 10 mesh of continuum shell elements, respectively. The shear stress computed using CTSHR13 is continuous across the continuum shell element interfaces. In addition, while the estimates of the transverse shear distributions using SSAVG4 and CTSHR13 (shown in Figure 7 and Figure 8) are both good, using CTSHR13 requires a mesh of only 3 continuum shell elements through the thickness, as compared to 24 elements for SSAVG4. Figure 10 compares the transverse shear stress distribution obtained with the solid element model with the shell element result. The figure shows that the transverse shear stresses predicted by solid elements do not vanish at the free surfaces of the structure. It also shows that the stress is discontinuous at layer interfaces. The reason for this is that in the composite solid element, the transverse shear stresses are obtained directly from the displacement field in contrast to the shell element, where the transverse shear stresses are obtained from an equilibrium calculation. These deficiencies decrease if the number of solid elements used in the discretization through the section thickness is increased. Although the transverse shear stresses are inaccurate, the displacement field and components of stress in the plane of the layer (not shown here) are in much better agreement with the analytical result. In fact, these results are somewhat better than the results obtained with the S8R elements. The composite solid elements were not used to analyze the thinner plates since the solid elements would not have any advantage over plate elements in that case. For 10, Figure 11 and Figure 12 show that the transverse shear and axial stress distributions of the finite element results—along with the CPT predictions—agree with elasticity theory. The stress distributions become more accurate with increasing span-to-thickness ratio (as the plate becomes thinner in comparison to the span). In Figure 13 and Figure 14 the maximum stress theory and Tsai-Wu theory failure indices are plotted as a function of the normalized distance from the midsurface for the two- and three-layer cases, respectively. The indices are calculated at the center of the plate for S8R elements with 4. Values of the failure index greater than or equal to 1.0 indicate failure. Discontinuous jumps in the failure index occur at layer boundaries as a result of the orientation of the material. The strain levels are well below those required for failure, so no strain-based failure indices are plotted. Abaqus/Explicit resultsThe explicit dynamic analysis is run for a sufficiently long time so that a quasi-static state is reached—that is, the plates are in steady-state vibration. Since step loadings are applied, static solutions of stresses can be obtained as half of their vibration amplitudes. Figure 15 and Figure 16 show the transverse shear stress (TSHR13) and the axial stress (S11) distributions through the plate thickness for the two-layer S4R model normalized as: and compared with classical plate theory (CPT) and linear elasticity theory. Figure 17 and Figure 18 show the corresponding results for the three-layer plate. In Figure 19 and Figure 20, the maximum stress theory and Tsai-Wu theory failure indices are plotted as a function of the normalized distance from the midsurface for the two- and three-layer cases, respectively. The indices are calculated at the center of the plate. Values of the failure index greater than or equal to 1.0 indicate failure. Discontinuous jumps in the failure index occur at layer boundaries due to the orientation of the material. The strain levels are well below those required for failure, so no strain-based failure indices are plotted. Input filesAbaqus/Standard input files
Abaqus/Explicit input files
References
FiguresFigure 1. Composite plate subject to distributed loading.
Figure 2. Maximum deflection of two- and three-layer plates with various
span-to-thickness ratios;
Abaqus/Standard
analysis.
Figure 3. Transverse shear stress distribution through the thickness of a
two-layer plate (4);
Abaqus/Standard
analysis.
Figure 4. Axial stress distribution through the thickness of a two-layer plate
(4);
Abaqus/Standard
analysis.
Figure 5. Transverse shear stress distribution through the thickness of a
three-layer plate (4);
Abaqus/Standard
analysis.
Figure 6. Axial stress distribution through the thickness of a three-layer plate
(4);
Abaqus/Standard
analysis.
Figure 7. Comparison of the elasticity solution of the transverse shear stress
distribution in a three-layer plate to the output variable SSAVG4 with 24 SC8R elements stacked through the thickness;
Abaqus/Standard
analysis.
Figure 8. Comparison of the elasticity solution of the transverse shear stress
distribution in a three-layer plate to the output variable CTSHR13 with 3 SC8R elements stacked through the thickness;
Abaqus/Standard
analysis.
Figure 9. Comparison of the elasticity solution of the transverse shear stress
distribution in a two-layer plate to the output variable CTSHR13 with 2 SC8R elements stacked through the thickness;
Abaqus/Standard
analysis.
Figure 10. Transverse shear stress distribution through the thickness of a
three-layer plate (4):
shells versus solid elements;
Abaqus/Standard
analysis.
Figure 11. Transverse shear stress distribution through the thickness of a
three-layer plate (10);
Abaqus/Standard
analysis.
Figure 12. Axial stress distribution through the thickness of a three-layer plate
(10);
Abaqus/Standard
analysis.
Figure 13. Maximum stress theory and Tsai-Wu theory (0.0)
failure indices as a function of normalized distance from the midsurface.
Two-layer plate, 4;
Abaqus/Standard
analysis.
Figure 14. Maximum stress theory and Tsai-Wu theory (0.0)
failure indices as a function of normalized distance from the midsurface.
Three-layer plate, 4;
Abaqus/Standard
analysis.
Figure 15. Transverse shear stress distribution through the thickness of a
two-layer plate;
Abaqus/Explicit
analysis.
Figure 16. Axial stress distribution through the thickness of a two-layer plate;
Abaqus/Explicit
analysis.
Figure 17. Transverse shear stress distribution through the thickness of a
three-layer plate;
Abaqus/Explicit
analysis.
Figure 18. Axial stress distribution through the thickness of a three-layer
plate;
Abaqus/Explicit
analysis.
Figure 19. Maximum stress theory and Tsai-Wu theory failure indices as a function
of normalized distance from the midsurface. Two-layer plate;
Abaqus/Explicit
analysis.
Figure 20. Maximum stress theory and Tsai-Wu theory failure indices as a function
of normalized distance from the midsurface. Three-layer plate;
Abaqus/Explicit
analysis.
|