*Heading
Composite Cylinder
*Include, input=dsacentrifuge_node_elem.inp
*Elset, elset=_I1,  generate
 841,  936,    1
*Shell Section, elset=_I1, material=Aluminum
3., 3
*Elset, elset=_I2,  generate
 721,  840,    1
*Shell Section, elset=_I2, material=Aluminum
3., 3
*Elset, elset=_I3,  generate
 937,  960,    1
*Shell Section, elset=_I3, Composite, Offset=0.155 
    3.00, 3, Aluminum, 0.0
 <TT800>, 3,     T800, 00.0
<THM400>, 3,    HM400, <THETA1>
<THM400>, 3,    HM400, <THETA2>
 <TT800>, 3,     T800, 00.0
<THM400>, 3,    HM400, <THETA2>
<THM400>, 3,    HM400, <THETA1>
 <TT800>, 3,     T800, 00.0
    1.00, 3,     T800, 90.0
*Elset, elset=_I4,  generate
 145,  720,    1
*Shell Section, elset=_I4, Composite 
 <TT800>, 3,     T800, 00.0
<THM400>, 3,    HM400, <THETA1>
<THM400>, 3,    HM400, <THETA2>
 <TT800>, 3,     T800, 00.0
<THM400>, 3,    HM400, <THETA2>
<THM400>, 3,    HM400, <THETA1>
 <TT800>, 3,     T800, 00.0
*Elset, elset=_I5,  generate
 961,  984,    1
*Shell Section, elset=_I5, Composite, Offset=0.092
    2.00, 3, Aluminum, 0.0
 <TT800>, 3,     T800, 00.0
<THM400>, 3,    HM400, <THETA1>
<THM400>, 3,    HM400, <THETA2>
 <TT800>, 3,     T800, 00.0
<THM400>, 3,    HM400, <THETA2>
<THM400>, 3,    HM400, <THETA1>
 <TT800>, 3,     T800, 00.0
    1.00, 3,     T800, 90.0
*Elset, elset=_I6,  generate
  49,  144,    1
*Shell Section, elset=_I6, material=Aluminum
2., 3
*Elset, elset=_I7,  generate
  1,  48,   1
*Shell Section, elset=_I7, material=Aluminum
2., 3
*Nset, nset=_G4
   1,  33,  34,  35,  36,  37,  38,  39,  40,  
  41,  42,  43,  44,  45,  46,  47,  48,  49,  
  50,  51,  52,  53,  54,  55 
*Nset, nset=_G5
   8, 204, 205, 206, 207, 208, 209, 210, 211, 
 212, 213, 214, 215, 216, 217, 218, 219, 220, 
 221, 222, 223, 224, 225, 226
*Elset, elset=_G4  
   1,   3,   5,   7,   9,  11,  13,  15,  17,  
  19,  21,  23,  25,  27,  29,  31,  33,  35,  
  37,  39,  41,  43,  45,  47, 844, 848, 852, 
 856, 860, 864, 868, 872, 876, 880, 884, 888, 
 892, 896, 900, 904, 908, 912, 916, 920, 924, 
 928, 932, 936
*Nset, nset=_G9,  generate
    1,  1008,     1
*Elset, elset=_G9,   generate
   1,  984,    1
** 
** MATERIALS
** 
*Material, name=Aluminum
*Density
 2.854e-09,
*Elastic
70000., 0.33
*Material, name=HM400
*Density
 1.63e-09,
*Elastic, type=LAMINA
233967.,14778., 0.032,10191., 5777., 5634.
*Material, name=T800
*Density
 1.63e-09,
*Elastic, type=LAMINA
177000.,14920.,  0.21, 5700., 5700., 5630.
*Parameter
THETA = 45.0
THETA1= THETA
THETA2=-THETA
TT800=0.15
THM400=0.50
ALPHA=0.0
*Design Parameter
THM400,THETA,TT800,ALPHA
*Parameter Shape Variation,Parameter=ALPHA,Step=1,
File=dsacentrifuge_freq,mode=4
*Dsa controls
TT800, CD, 0.22500E-03
*Elset, elset=slice, generate
  5, 6, 1
 57,60, 1
193,216,1
731,735,1
849,852,1
938,938,1
982,982,1
*Elset, elset=QA_TEST_ESET
877
*Nset, nset=cut, generate
 76,  76, 1
145, 145, 1
434, 456, 1
110, 110, 1
195, 195, 1
*Nset, nset=QA_TEST_NSET
4,
*Step,  nlgeom, dsa
Revolve The cylinder
*Static
0.2, 1., 1e-05, 1.
** 
** BOUNDARY CONDITIONS
** 
** Name: BC-1 Type: Displacement/Rotation
*Boundary
_G4, 1, 3
_G5, 2, 3
** 
** LOADS
** 
** Name: Centrif   Type: CENTRIFUGAL
*Dload
_G9,CENTRIF,1096620.00,0.0,0.0,0.0,1.0,0.0,0.0
** 
** OUTPUT REQUESTS
** 
*Output, field, Freq=5
*Node Output
U,
*Element Output
1,2,3,4,5,6,7,8,9,10,11
S,
SF,
*Design Response, Freq=5
*Node Response
U,
*Element Response
1,2,3,4,5,6,7,8,9,10,11
S,
SF,
*Design Response, Freq=5
*Node Response, Nset=QA_TEST_NSET
U,
*Element Response, Elset=QA_TEST_ESET
1,2,3,4,5,6,7,8,9,10,11
S,
*El Print, freq=0
*Node Print, freq=0
*End Step