Thin Walled Cylinder Buckling

A cylinder is loaded with a uniformly distributed line load along the top edge. Critical buckling load is determined. Different meshes and differences between Belytschko-Tsay shell (elform 2), S/R Hughes-Liu shell (elform 6), Belytschko-Wong-Chiang shell (elform 10) and fully integrated shell (elform 16) can be studied.

*BOUNDARY_SPC_SET
*CONTROL_IMPLICIT_BUCKLE
*CONTROL_IMPLICIT_GENERAL
*CONTROL_IMPLICIT_SOLUTION
*CONTROL_IMPLICIT_EIGENVALUE
*CONTROL_SHELL
*CONTROL_TERMINATION
*DATABASE_BINARY_D3PLOT
*DEFINE_CURVE
*ELEMENT_SHELL
*END
*HOURGLASS
*KEYWORD
*LOAD_NODE_SET
*MAT_ELASTIC
*NODE
*PART
*SECTION_SHELL
*SET_NODE_LIST_TITLE
*TITLE

*TITLE
Buckling of a Thin Walled Cylinder Under Compression
*CONTROL_IMPLICIT_GENERAL
         1  0.100000         2         1         2
*CONTROL_IMPLICIT_SOLUTION
         2        11        15  0.001000  0.010000       0.0  0.900000  1.000000
         2         1         1         2
         0         1       0.0         1         2
*CONTROL_IMPLICIT_BUCKLE
         3
*CONTROL_TERMINATION
  1.000000         0       0.0       0.0       0.0
*CONTROL_SHELL
  20.00000         0         0         0         2         1         1         1
       0.0         0
*DATABASE_BINARY_D3PLOT
  0.010000
*PART
         1         1         1         0         1
*SECTION_SHELL
         1         2       0.0         0         1       0.0         0         1
  0.100000  0.100000  0.100000  0.100000         0       0.0
*MAT_ELASTIC
         1  0.0100001.0000e+07  0.300000       0.0       0.0       0.0
*HOURGLASS
         1         4       0.0         0       0.0       0.0       0.0       0.0
*SET_NODE_LIST_TITLE
bottom nodes
         1       0.0       0.0       0.0       0.0
         1        31        61        91       121       151       181       211
       241       271       301       331       361       391       421       451
       481       511       541       571       601       631       661       691
       721       751       781       811       841       871       901       931
       961       991      1021      1051      1081      1111      1141      1171
      1201      1231      1261      1291      1321      1351      1381      1411
      1441      1471      1501      1531      1561      1591      1621      1651
      1681      1711      1741      1771      1801      1831      1861      1891
      1921      1951      1981      2011      2041      2071      2101      2131
      2161      2191      2221      2251
*SET_NODE_LIST_TITLE
top nodes
         2       0.0       0.0       0.0       0.0
        30        60        90       120       150       180       210       240
       270       300       330       360       390       420       450       480
       510       540       570       600       630       660       690       720
       750       780       810      1200      1230      1260      1290      1320
      1350      1380      1410      1440      1470      1500      1530      1560
      1590      1620      1650      1680      1710      1740      1920      1950
      1980      2010      2040      2070      2100      2130      2160      2190
      2220      2250      2280       840       870       900       930       960
       990      1020      1050      1080      1110      1140      1170      1770
      1800      1830      1860      1890
*BOUNDARY_SPC_SET
         1         0         1         1         1         1         1         1
*BOUNDARY_SPC_SET
         2         0         1         1         0         0         0         0
*DEFINE_CURVE
         1         0       0.0       0.0       0.0       0.0
                 0.0                 0.0
          1.00000000         13.15789474
*LOAD_NODE_SET
         2         3         1 -1.000000
*ELEMENT_SHELL
*NODE
*END