{"id":3221,"date":"2024-01-25T15:33:44","date_gmt":"2024-01-25T21:33:44","guid":{"rendered":"https:\/\/kimeca.com.mx\/?p=3221"},"modified":"2026-01-04T18:24:47","modified_gmt":"2026-01-05T00:24:47","slug":"lamin-comp-shells-buck-of-a-cylindrical","status":"publish","type":"post","link":"https:\/\/kimeca.com.mx\/index.php\/lamin-comp-shells-buck-of-a-cylindrical\/","title":{"rendered":"Laminated composite shells: buckling of a cylindrical panel with a circular hole with Abaqus"},"content":{"rendered":"

<\/span>Laminated composite shells: buckling of a cylindrical panel with a circular hole with Abaqus<\/a><\/span><\/h2>\n

This example illustrates modeling a thin, laminated composite shell in the presence of buckling.<\/a><\/p>\n

SIMULIA Community<\/a><\/p>\n

Training Course: Analysis of Composite Materials with Abaqus<\/a><\/p>\n

\"Laminated<\/p>\n

Animation-1.<\/em><\/strong><\/p>\n

<\/h2>\n

<\/span>Geometry and model<\/strong><\/span><\/h2>\n
\n

The structure analyzed is shown in\u00a0Figure 1\u00a0and was originally studied experimentally by Knight and Starnes (1984). The test specimen is a cylindrical panel with a 355.6 mm (14 in) square platform and a 381 mm (15 in) radius of curvature, so that the panel covers a 55.6\u00b0 arc of the cylinder. The panel contains a centrally located hole of 50.8 mm (2 in) diameter. The shell consists of 16 layers of unidirectional graphite fibers in an epoxy resin. Each layer is 0.142 mm (.0056 in) thick. The layers are arranged in the symmetric stacking sequence {\u00b1<\/span><\/span><\/span><\/span><\/span><\/span>45\/90\/0\/0\/90\/\u2213<\/span><\/span><\/span><\/span><\/span><\/span>45} degrees repeated twice. The nominal orthotropic elastic material properties as defined by Stanley (1985) are<\/p>\n

\"\"<\/p>\n

where the 1-direction is along the fibers, the 2-direction is transverse to the fibers in the surface of the lamina, and the 3-direction is normal to the lamina.<\/p>\n

The panel is fully clamped on the bottom edge, clamped except for axial motion on the top edge and simply supported along its vertical edges. Three analyses are considered. The first is a linear (prebuckling) analysis in which the panel is subjected to a uniform end shortening of 0.8 mm (.0316 in). The total axial force and the distribution of axial force along the midsection are used to compare the results with those obtained by Stanley (1985). The second analysis consists of an eigenvalue extraction of the first five buckling modes. The buckling loads and mode shapes are also compared with those presented by Stanley (1985). Finally, a nonlinear load-deflection analysis is done to predict the postbuckling behavior, using the modified Riks algorithm. For this analysis an initial imperfection is introduced. The imperfection is based on the fourth buckling mode extracted during the second analysis. These results are compared with those of Stanley (1985) and with the experimental measurements of Knight and Starnes (1984).<\/p>\n

The mesh used in\u00a0Abaqus<\/span>\u00a0is shown in\u00a0Figure 2. The anisotropic material behavior precludes any symmetry assumptions, hence the entire panel is modeled. The same mesh is used with the 4-node shell element (type\u00a0S4R5<\/span>) and also with the 9-node shell element (type\u00a0S9R5<\/span>); the 9-node element mesh, thus, has about four times the number of degrees of freedom as the 4-node element mesh. The 6-node triangular shell element\u00a0STRI65<\/span>\u00a0is also used; it employs two triangles for each quadrilateral element of the second-order mesh. Mesh generation is facilitated by specifying node fill and node mapping, as shown in the input data. In this model specification of the relative angle of orientation to define the material orientation within each layer, along with orthotropic elasticity in plane stress, makes the definition of the laminae properties straightforward.<\/p>\n

The shell elements used in this example use an approximation to thin shell theory, based on a numerical penalty applied to the transverse shear strain along the element edges. These elements are not universally applicable to the analysis of composites since transverse shear effects can be significant in such cases and these elements are not designed to model them accurately. Here, however, the geometry of the panel is that of a thin shell; and the symmetrical lay-up, along with the relatively large number of laminae, tends to diminish the importance of transverse shear deformation on the response.<\/p>\n<\/div>\n

<\/h2>\n

 <\/p>\n

<\/span>Relation between stress resultants and generalized strains<\/strong><\/span><\/h2>\n
\n

The shell section is most easily defined by giving the layer thickness, material, and orientation, in which case\u00a0Abaqus<\/span>\u00a0preintegrates to obtain the section stiffness properties. However, the user can choose to input the section stiffness properties directly instead, as follows.<\/p>\n

In\u00a0Abaqus<\/span>\u00a0a lamina is considered as an orthotropic sheet in plane stress. The principal material axes of the lamina (see\u00a0Figure 3) are longitudinal, denoted by\u00a0L<\/span>; transverse to the fiber direction in the surface of the lamina, denoted by\u00a0T<\/span>; and normal to the lamina surface, denoted by\u00a0N<\/span><\/span><\/span><\/span><\/span><\/em>.<\/span><\/span><\/span>\u00a0The constitutive relations for a general orthotropic material in the principal directions (L<\/span><\/span>,<\/span><\/span>T<\/span><\/span>,<\/span><\/span>N<\/span><\/span><\/span><\/span><\/span><\/span><\/span><\/em>) are<\/p>\n

\"\"<\/p>\n

In terms of the data required to define orthotropic elasticity by specifying terms in the elastic stiffness matrix in\u00a0Abaqus<\/span>\u00a0these are<\/p>\n

\"\"<\/p>\n

This matrix is symmetric and has nine independent constants. If we assume a state of plane stress, then\u00a0\u03c3<\/span><\/span><\/span>N<\/span><\/span><\/span><\/span><\/span><\/span><\/span><\/span><\/em>\u00a0is taken to be zero. This yields<\/p>\n

\"\"<\/p>\n

where<\/p>\n

\"\"<\/p>\n

The correspondence between these terms and the usual engineering constants that might be given for a simple orthotropic layer in a laminate is<\/p>\n

\"\"<\/p>\n

The parameters used on the right-hand side of the above equation are those that must be provided as part of the definition of orthotropic elasticity in plane stress.<\/p>\n

If the (1<\/span><\/span>,<\/span><\/span>2<\/span><\/span>,<\/span><\/span>N<\/span><\/span><\/span><\/span><\/span><\/span><\/span><\/em>) system denotes the standard shell basis directions that\u00a0Abaqus<\/span>\u00a0chooses by default, the local stiffness components must be rotated to this system to construct the lamina’s contribution to the general shell section stiffness. Since\u00a0Q<\/span><\/span><\/span>i<\/span><\/span>j<\/span><\/span><\/span><\/span><\/span><\/span><\/span><\/span><\/span><\/em>\u00a0represent fourth-order tensors, in the case of a lamina they are oriented at an angle\u00a0\u03b8<\/span><\/span><\/span><\/span><\/span><\/span><\/em>\u00a0to the standard shell basis directions used in\u00a0Abaqus<\/span>. Hence, the transformation is<\/p>\n

\"\"<\/p>\n

where Qij<\/i><\/span>\u00a0are the stiffness coefficients in the standard shell basis directions used by\u00a0Abaqus<\/span>.<\/p>\n

Abaqus<\/span>\u00a0assumes that a laminate is a stack of laminae arranged with the principal directions of each layer in different orientations. The various layers are assumed to be rigidly bonded together. The section force and moment resultants per unit length in the normal basis directions in a given layer can be defined on this basis as<\/p>\n

\"\"<\/p>\n

where\u00a0h<\/span><\/em>\u00a0is the thickness of the layer.<\/p>\n

This leads to the relations<\/p>\n

\"\"<\/p>\n

where the components of this section stiffness matrix are given by<\/p>\n

\"\"<\/p>\n

\"\"<\/p>\n

\"\"<\/p>\n

\"\"<\/p>\n

\"\"<\/p>\n<\/div>\n

<\/h2>\n

 <\/p>\n

\n

<\/span>Results and discussion<\/strong><\/span><\/h2>\n
\n

The total axial force necessary to compress the panel 0.803 mm (0.0316 in) is 100.2 kN (22529 lb) for the mesh of\u00a0S9R5<\/span>\u00a0elements, 99.5 kN (22359 lb) for the mesh of\u00a0S4R5<\/span>\u00a0elements, and 100.3 kN (22547 lb) for the mesh of\u00a0STRI65<\/span>\u00a0elements. These values match closely with the result of 100 kN (22480 lb) reported by Stanley (1985).\u00a0Figure 5\u00a0shows the displaced configuration and a profile of axial force along the midsection of the panel (at\u00a0z<\/span><\/span>=<\/span><\/span><\/em>L<\/span><\/span>\/<\/span><\/span><\/em>2<\/em>)<\/span><\/span><\/span><\/span><\/span><\/span><\/span><\/span>. It is interesting to note that the axial load is distributed almost evenly across the entire panel, with only a very localized area near the hole subjected to an amplified stress level. This suggests that adequate results for this linear analysis could also be obtained with a coarser mesh that has a bias toward the hole.<\/p>\n

The second stage of the analysis is the eigenvalue buckling prediction. To obtain the buckling predictions with\u00a0Abaqus<\/span>, an eigenvalue buckling prediction step is run. In this step nominal values of load are applied. The magnitude that is used is not of any significance, since eigenvalue buckling is a linear perturbation procedure: the stiffness matrix and the stress stiffening matrix are evaluated at the beginning of the step without any of this load applied. The eigenvalue buckling prediction step calculates the eigenvalues that, multiplied with the applied load and added to any \u201cbase state\u201d loading, are the predicted buckling loads. The eigenvectors associated with the eigenvalues are also obtained.<\/p>\n

The buckling predictions are summarized in\u00a0Table 1\u00a0and\u00a0Figure 6. The buckling load predictions from\u00a0Abaqus<\/span>\u00a0are higher than those reported by Stanley. The eigenmode predictions given by the mesh using element types\u00a0S4R5<\/span>,\u00a0S9R5<\/span>, and\u00a0STRI65<\/span>\u00a0are all the same and agree well with those reported by Stanley. Stanley makes several important observations that remain valid for the\u00a0Abaqus<\/span>\u00a0results: (1) the eigenvalues are closely spaced; (2) nevertheless, the mode shapes vary significantly in character; (3) the first buckling mode bears the most similarity to the linear prebuckling solution; (4) there is no symmetry available that can be utilized for computational efficiency.<\/p>\n

Following the eigenvalue buckling analyses, nonlinear postbuckling analysis is carried out by imposing an imperfection based on the fourth buckling mode. The maximum initial perturbation is 10% of the thickness of the shell. The load versus normalized displacement plots for the\u00a0S9R5<\/span>\u00a0mesh, the\u00a0S4R5<\/span>\u00a0mesh, and the\u00a0STRI65<\/span>\u00a0mesh are compared with the experimental results and those given by Stanley in\u00a0Figure 7. The overall response prediction is quite similar for the\u00a0Abaqus<\/span>\u00a0elements, although the general behavior predicted by Stanley is somewhat different. The\u00a0Abaqus<\/span>\u00a0results show a peak load slightly above the buckling load predicted by the eigenvalue extraction, while Stanley’s results show a significantly lower peak load. In addition, the\u00a0Abaqus<\/span>\u00a0results show rather less loss of strength after the initial peak, followed quite soon by positive stiffness again. Neither the\u00a0Abaqus<\/span>\u00a0results nor Stanley’s results agree closely with the experimentally observed dramatic loss of strength after peak load. Stanley ascribes this to material failure (presumably delamination), which is not modeled in his analyses or in these.<\/p>\n

Figure 8\u00a0shows the deformed configurations for the panel during its postbuckling response. The plots show the results for\u00a0S4R5<\/span>, but the pattern is similar for\u00a0S9R5<\/span>\u00a0and\u00a0STRI65<\/span>. The response is quite symmetric initially; but, as the critical load is approached, a nonsymmetric dimple develops and grows, presumably accounting for the panel’s loss of strength. Later in the postbuckling response another wrinkle can be seen to be developing.<\/p>\n

\n<\/div>\n<\/div>\n

\"\"Animation-2.<\/strong><\/em><\/p>\n

<\/h2>\n

 <\/p>\n

\n

<\/span>Tables<\/strong><\/span><\/h2>\n
\n\n\n\n\n\n\n\n\n\n\n\n\n\n\n\n\n\n\n\n\n\n\n\n\n\n\n\n
Table 1.\u00a0<\/span>Summary of buckling load predictions.<\/span><\/caption>\n
Mode 1<\/td>\nStanley<\/td>\n107.0 kN (24054 lb)<\/td>\n<\/tr>\n
S9R5<\/td>\n113.4 kN (25501 lb)<\/td>\n<\/tr>\n
S4R5<\/td>\n115.5 kN (25964 lb)<\/td>\n<\/tr>\n
S4<\/td>\n114.3 kN (25696 lb)<\/td>\n<\/tr>\n
STRI65<\/td>\n113.8 kN (25579 lb)<\/td>\n<\/tr>\n
Mode 2<\/td>\nStanley<\/td>\n109.6 kN (24638 lb)<\/td>\n<\/tr>\n
S9R5<\/td>\n117.6 kN (26429 lb)<\/td>\n<\/tr>\n
S4R5<\/td>\n121.2 kN (27244 lb)<\/td>\n<\/tr>\n
S4<\/td>\n116.5 kN (26196 lb)<\/td>\n<\/tr>\n
STRI65<\/td>\n117.8 kN (26492 lb)<\/td>\n<\/tr>\n
Mode 3<\/td>\nStanley<\/td>\n116.2 kN (26122 lb)<\/td>\n<\/tr>\n
S9R5<\/td>\n120.3 kN (27049 lb)<\/td>\n<\/tr>\n
S4R5<\/td>\n124.7 kN (28042 lb)<\/td>\n<\/tr>\n
S4<\/td>\n124.1 kN (27889 lb)<\/td>\n<\/tr>\n
STRI65<\/td>\n121.1 kN (27217 lb)<\/td>\n<\/tr>\n
Mode 4<\/td>\nStanley<\/td>\n140.1 kN (31494 lb)<\/td>\n<\/tr>\n
S9R5<\/td>\n147.5 kN (33161 lb)<\/td>\n<\/tr>\n
S4R5<\/td>\n156.1 kN (35092 lb)<\/td>\n<\/tr>\n
S4<\/td>\n152.3 kN (34247 lb)<\/td>\n<\/tr>\n
STRI65<\/td>\n146.9 kN (33015 lb)<\/td>\n<\/tr>\n
Mode 5<\/td>\nStanley<\/td>\n151.3 kN (34012 lb)<\/td>\n<\/tr>\n
S9R5<\/td>\n171.3 kN (38512 lb)<\/td>\n<\/tr>\n
S4R5<\/td>\n181.5 kN (40800 lb)<\/td>\n<\/tr>\n
S4<\/td>\n184.2 kN (41413 lb)<\/td>\n<\/tr>\n
STRI65<\/td>\n172.8 kN (38843 lb)<\/td>\n<\/tr>\n<\/tbody>\n<\/table>\n<\/div>\n<\/article>\n
\n

\n\n<\/map>\n<\/p>\n<\/article>\n

<\/h2>\n

 <\/p>\n

\n

<\/span>Figures<\/strong><\/span><\/h2>\n
\n
Figure 1. Geometry for cylindrical panel with hole.<\/span>
\"\"<\/figure>\n
Figure 2. Mesh for cylindrical panel with hole.<\/span>
\"\"<\/figure>\n
Figure 3. Typical lamina.<\/span>
\"\"<\/figure>\n
Figure 4. Typical laminate.<\/span>
\"\"<\/figure>\n
Figure 5. Displaced shape and axial force distribution.<\/span>
\"\"<\/figure>\n
Figure 6. Buckling modes, element types\u00a0S4R5<\/span>,\u00a0S9R5<\/span>, and\u00a0STRI65<\/span>.<\/span>
\"\"<\/figure>\n
Figure 7. Load-displacement response.<\/span>
\"\"<\/figure>\n
Figure 8. Postbuckling deformations: 10% h imperfection with\u00a0S4R5<\/span>.<\/span>
\"\"<\/figure>\n<\/div>\n<\/article>\n","protected":false},"excerpt":{"rendered":"

Laminated composite shells: buckling of a cylindrical panel with a circular hole with Abaqus<\/p>\n","protected":false},"author":11,"featured_media":3218,"comment_status":"open","ping_status":"open","sticky":false,"template":"","format":"standard","meta":{"_et_pb_use_builder":"","_et_pb_old_content":"","_et_gb_content_width":"","footnotes":"","jetpack_publicize_message":"","jetpack_publicize_feature_enabled":true,"jetpack_social_post_already_shared":true,"jetpack_social_options":{"image_generator_settings":{"template":"highway","default_image_id":0,"font":"","enabled":false},"version":2}},"categories":[46,34,44,130,36,132],"tags":[84,126,128,106,68,124,58],"class_list":["post-3221","post","type-post","status-publish","format-standard","has-post-thumbnail","hentry","category-abaqus","category-dassault-systemes","category-kimeca","category-service","category-simulia","category-training","tag-simulia-abaqus-fea","tag-buckling-failure","tag-composite-materials","tag-dassault-systemes","tag-fea","tag-kimeca","tag-simulia"],"jetpack_publicize_connections":[],"_links":{"self":[{"href":"https:\/\/kimeca.com.mx\/index.php\/wp-json\/wp\/v2\/posts\/3221","targetHints":{"allow":["GET"]}}],"collection":[{"href":"https:\/\/kimeca.com.mx\/index.php\/wp-json\/wp\/v2\/posts"}],"about":[{"href":"https:\/\/kimeca.com.mx\/index.php\/wp-json\/wp\/v2\/types\/post"}],"author":[{"embeddable":true,"href":"https:\/\/kimeca.com.mx\/index.php\/wp-json\/wp\/v2\/users\/11"}],"replies":[{"embeddable":true,"href":"https:\/\/kimeca.com.mx\/index.php\/wp-json\/wp\/v2\/comments?post=3221"}],"version-history":[{"count":0,"href":"https:\/\/kimeca.com.mx\/index.php\/wp-json\/wp\/v2\/posts\/3221\/revisions"}],"wp:featuredmedia":[{"embeddable":true,"href":"https:\/\/kimeca.com.mx\/index.php\/wp-json\/wp\/v2\/media\/3218"}],"wp:attachment":[{"href":"https:\/\/kimeca.com.mx\/index.php\/wp-json\/wp\/v2\/media?parent=3221"}],"wp:term":[{"taxonomy":"category","embeddable":true,"href":"https:\/\/kimeca.com.mx\/index.php\/wp-json\/wp\/v2\/categories?post=3221"},{"taxonomy":"post_tag","embeddable":true,"href":"https:\/\/kimeca.com.mx\/index.php\/wp-json\/wp\/v2\/tags?post=3221"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}