*HEADING
: SHELL OFFSET WITH REBARS FOR LINEAR DYNAMICS [S8R]
*RESTART,WRITE,FREQ=999
*NODE,NSET=ALL
1,0,0,-0.1
2,10.,0,-0.1
3,10.,1.0,-0.1
4,0.,1.0,-0.1
5,5.,0,-0.1
6,10.,0.5,-0.1
7,5.,1.0,-0.1
8,0.,0.5,-0.1
9,0,0,0.1
10,10.,0,0.1
11,10.,1.0,0.1
12,0.,1.0,0.1
13,5.,0,0.1
14,10.,0.5,0.1
15,5.,1.0,0.1
16,0.,0.5,0.1
*NSET,NSET=SHELLB
1,2,3,4,5,6,7,8
*NSET,NSET=SHELLT
9,10,11,12,13,14,15,16
*NSET,NSET=LEFT
1,4,8,9,12,16
*NSET,NSET=BOT
1,2,5,9,10,13
*ELEMENT,TYPE=S8R,ELSET=SHELLB
1, 1,2,3,4,5,6,7,8
*ELEMENT,TYPE=S8R,ELSET=SHELLT
2,9,10,11,12,13,14,15,16
*SHELL SECTION,MATERIAL=ELMAT, ELSET=SHELLB,OFFSET=SNEG
0.2,3
*REBAR LAYER
REBARB,0.1,0.025,-0.05,RBMAT1,0,1
*SHELL SECTION,MATERIAL=ELMAT, ELSET=SHELLT,OFFSET=0.5
0.2,3
*REBAR LAYER
REBART,0.1,0.025,0.05,RBMAT1,0,1
*MATERIAL,NAME=ELMAT
*ELASTIC
3.E6,0.25
*DENSITY
7.5E-4, 
*MATERIAL,NAME=RBMAT1
*ELASTIC
30.E6, 
*DENSITY
7.5E-4, 
*NSET,NSET=RIGHT
2,3,6,10,11,14
*BOUNDARY
LEFT,1,6
RIGHT,1,4
RIGHT,6,6
*SPECTRUM,NAME=SPEC,TYPE=ACCELERATION
 0., 0., 0.
 1.E5, 1000., 0.
 1.E5, 10000., 0.
*STEP
*FREQUENCY, eigensolver=lanczos
2, , 
*MODAL FILE
*END STEP
******************************************
*STEP
*STEADY STATE DYNAMICS,DIRECT
0.01, 
0.1, 
1., 
*CLOAD
RIGHT,5,-100.
*END STEP
*STEP
*MODAL DYNAMIC
.001,.003
*CLOAD
 RIGHT,5,-100.
*OUTPUT,FIELD
*ELEMENT OUTPUT,ELSET=SHELLB
 1,2,3
 E,Inten
*ELEMENT OUTPUT,ELSET=SHELLT
 1,2,3
 E,Inten
*ELEMENTOUTPUT,REBAR
Inten,
*NODE OUTPUT,NSET=SHELLB
 U, 
*NODE OUTPUT,NSET=SHELLT
 U,
*END STEP
******************************************
*STEP
*RESPONSE SPECTRUM,SUM=ABS,COMP=ALG
 SPEC,0.,0.,1.,386.4
*END STEP
*STEP
*FREQUENCY, eigensolver=lanczos
30, , 
*NODE FILE,F=0
*EL FILE,F=0
*NODE PRINT,F=0
*EL PRINT,F=0
*END STEP
**
*STEP
*STEADY STATE DYNAMICS,SUBSPACE PROJECTION
0.01, 
0.1, 
1., 
*CLOAD
RIGHT,5,-100.
*EL FILE,ELSET=SHELLB
 1,2,3
 S,
 E,
*EL FILE,REBAR=REBARB,ELSET=SHELLB
 S,E
*EL PRINT,ELSET=SHELLB
 1,2,3
 S,
 E,
*EL PRINT,REBAR=REBARB,ELSET=SHELLB
 S,E
*NODE PRINT,NSET=SHELLB
 U,
 RF,
*NODE FILE,NSET=SHELLB
 U, RF,
*EL FILE,ELSET=SHELLT
 1,2,3
 S,
 E,
*EL FILE,REBAR=REBART,ELSET=SHELLT
S,E
*EL PRINT,ELSET=SHELLT
 1,2,3
 S,
 E,
*EL PRINT,REBAR=REBART,ELSET=SHELLT
S,E
*NODE PRINT,NSET=SHELLT
 U,
 RF,
*NODE FILE,NSET=SHELLT
 U, RF,
*OUTPUT,FIELD
*ELEMENTOUTPUT,REBAR
Inten,
*ELEMENT OUTPUT,ELSET=SHELLB
 1,2,3
 S,
 E,Inten
*ELEMENT OUTPUT,ELSET=SHELLT
 1,2,3
 S,
 E,Inten
*NODE OUTPUT,NSET=SHELLB
 U, RF,
*NODE OUTPUT,NSET=SHELLT
 U, RF,
*OUTPUT,HISTORY
*ELEMENT OUTPUT,ELSET=SHELLB
 1,2,3
 S,
 E,
*ELEMENT OUTPUT,ELSET=SHELLT
 1,2,3
 S,
 E,
*NODE OUTPUT,NSET=SHELLB
 U, RF,
*NODE OUTPUT,NSET=SHELLT
 U, RF,
*END STEP