*HEADING
 BEAM ELEMENTS: MODEL CHANGE VERIFICATION, ELASTIC MATERIAL
*PREPRINT,MODEL=YES,HISTORY=NO
**
** TEST 1: B21
**
*NODE,NSET=NALL
1,  0.0, 1.
2,  1.0, 1.
3,  2.0, 1.
11, 0.0, 0.
12, 1.5, -.1
13, 3.0, -.3
*NSET,NSET=LEFT
1,11
*NSET,NSET=MID
12, 
*NSET,NSET=RIGHT
13, 
*NSET,NSET=TIP
3,13
*ELEMENT,TYPE=B21,ELSET=TEST1
1, 1,2
2, 2,3
11,11,12
12,12,13
*ELSET,ELSET=REFER
1,2
*ELSET,ELSET=TEST
11,12
**
** TEST 2: B22
**
*NODE,NSET=NALL
101, 5.0, 1.
102, 6.0, 1.
103, 7.0, 1.
111, 5.0, 0.
112, 6.5, -.1
113, 8.0, -.3
*NSET,NSET=LEFT
101,111
*NSET,NSET=MID
112, 
*NSET,NSET=RIGHT
113, 
*NSET,NSET=TIP
103,113
*ELEMENT,TYPE=B22,ELSET=TEST2
101,101,102,103
111,111,112,113
*ELSET,ELSET=REFER
101,102
*ELSET,ELSET=TEST
111,112
**
** TEST 3: B23
**
*NODE,NSET=NALL
201, 10.0, 1.
202, 11.0, 1.
203, 12.0, 1.
211, 10.0, 0.
212, 11.5, -.1
213, 13.0, -.3
*NSET,NSET=LEFT
201,211
*NSET,NSET=MID
212, 
*NSET,NSET=RIGHT
213, 
*NSET,NSET=TIP
203,213
*ELEMENT,TYPE=B23,ELSET=TEST3
201,201,202
202,202,203
211,211,212
212,212,213
*ELSET,ELSET=REFER
201,202
*ELSET,ELSET=TEST
211,212
**
** TEST 4: B21H
**
*NODE,NSET=NALL
301, 0.0, 4.
302, 1.0, 4.
303, 2.0, 4.
311, 0.0, 3.
312, 1.5, 2.9
313, 3.0, 2.7
*NSET,NSET=LEFT
301,311
*NSET,NSET=MID
312, 
*NSET,NSET=RIGHT
313, 
*NSET,NSET=TIP
303,313
*ELEMENT,TYPE=B21H,ELSET=TEST4
301, 301,302
302, 302,303
311, 311,312
312, 312,313
*ELSET,ELSET=REFER
301,302
*ELSET,ELSET=TEST
311,312
**
** TEST 5: B22H
**
*NODE,NSET=NALL
401, 5., 4.
402, 6., 4.
403, 7., 4.
411, 5.0, 3.
412, 6.5, 2.9
413, 8.0, 2.7
*NSET,NSET=LEFT
401,411
*NSET,NSET=MID
412, 
*NSET,NSET=RIGHT
413, 
*NSET,NSET=TIP
403,413
*ELEMENT,TYPE=B22H,ELSET=TEST5
401, 401,402,403
411, 411,412,413
*ELSET,ELSET=REFER
401,402
*ELSET,ELSET=TEST
411,412
**
** TEST 6: B23H
**
*NODE,NSET=NALL
501, 10.0,  4.0
502, 11.0,  4.0
503, 12.0,  4.0
511, 10.0,  3.0
512, 11.5,  2.9
513, 13.0,  2.7
*NSET,NSET=LEFT
501,511
*NSET,NSET=MID
512, 
*NSET,NSET=RIGHT
513, 
*NSET,NSET=TIP
503,513
*ELEMENT,TYPE=B23H,ELSET=TEST6
501, 501,502
502, 502,503
511, 511,512
512, 512,513
*ELSET,ELSET=REFER
501,502
*ELSET,ELSET=TEST
511,512
**
** TEST 7: B21 WITH REBAR
**
*NODE,NSET=NALL
601, 0.0,  7.0
602, 1.0,  7.0
603, 2.0,  7.0
611, 0.0,  6.0
612, 1.5,  5.9
613, 3.0,  5.7
*NSET,NSET=LEFT
601,611
*NSET,NSET=MID
612, 
*NSET,NSET=RIGHT
613, 
*NSET,NSET=TIP
603,613
*ELEMENT,TYPE=B21,ELSET=TEST7
601, 601,602
602, 602,603
611, 611,612
612, 612,613
*ELSET,ELSET=REFER
601,602
*ELSET,ELSET=TEST
611,612
*REBAR,ELEMENT=BEAM,NAME=REBAR7,MATERIAL=ELAS
TEST7,.5,0.,.75
**
** TEST 8: B23 WITH REBAR
**
*NODE,NSET=NALL
701, 5.0,  7.0
702, 6.0,  7.0
703, 7.0,  7.0
711, 5.0,  6.0
712, 6.5,  5.9
713, 8.0,  5.7
*NSET,NSET=LEFT
701,711
*NSET,NSET=MID
712, 
*NSET,NSET=RIGHT
713, 
*NSET,NSET=TIP
703,713
*ELEMENT,TYPE=B23,ELSET=TEST8
701, 701,702
702, 702,703
711, 711,712
712, 712,713
*ELSET,ELSET=REFER
701,702
*ELSET,ELSET=TEST
711,712
*REBAR,ELEMENT=BEAM,NAME=REBAR7,MATERIAL=ELAS
TEST8,.5,0.,.75
**
**  REMAINDER OF INPUT DECK SAME FOR ALL TESTS
**
*ELSET,ELSET=EALL
TEST,REFER
*BEAM SECTION,SECTION=RECT,TEMPERATURE=GRADIENTS,MATERIAL=ELAS,ELSET=EALL
1.,2.
0.,0.,-1.
*NORMAL
TEST,LEFT,0.,1.,0.
TEST,MID,.086,.9963,0.
TEST,RIGHT,.11723,.9931,0.
*MATERIAL,NAME=ELAS
*ELASTIC
1000., .25
*EXPANSION
1.E-3, 
*DENSITY
1.E-3, 
*INITIAL CONDITIONS,TYPE=TEMPERATURE
NALL,70.,0.
*AMPLITUDE,NAME=RAMP
0.,0.,  1.,1.
**
**
*STEP,NLGEOM
 STEP #1: DISPLACE TEST SPECIMEN SO FINAL SHAPE IS
          IDENTICAL TO INITIAL SHAPE OF REFERENCE ELEMENT 
*STATIC
0.25,1.
*PRINT,MODEL CHANGE=YES
*BOUNDARY,OP=NEW
LEFT,1,6
RIGHT,1,1,-1.0
RIGHT,2,2, 0.3
RIGHT,6,6,.11750075
MID,1,1,-0.5
MID,2,2, 0.1
MID,6,6,.086105945
*NODE PRINT,NSET=LEFT,FREQ=999
RF, 
*EL PRINT,POSITION=CENTROID,FREQ=999
1,
S,E
*OUTPUT,FIELD,FREQ=999
*NODE OUTPUT,NSET=LEFT
RF, 
*OUTPUT,HISTORY,FREQ=999
*NODE OUTPUT,NSET=LEFT
RF, 
*OUTPUT,FIELD,FREQ=999
*ELEMENT OUTPUT,ELSET=TEST
 S
*OUTPUT,FIELD,FREQ=999
*ELEMENT OUTPUT,ELSET=TEST,POSITION=CENTROID
E, S
*END STEP
**
**
*STEP,NLGEOM
 STEP #2: REMOVE TEST ELEMENTS
*STATIC
1.,1.
*MODEL CHANGE,REMOVE,TYPE=ELEMENT
TEST,
*END STEP
**
**
*STEP,NLGEOM
 STEP #3: ADD TEST ELEMENTS BACK IN STRAIN-FREE
*STATIC
1.,1.
*MODEL CHANGE,ADD=STRAINFREE,TYPE=ELEMENT
TEST, 
*BOUNDARY,OP=NEW
LEFT,1,6
*END STEP
**
**
*STEP,NLGEOM
 STEP #4: APPLY LOAD TO BOTH ELEMENTS 
*STATIC
0.25,1.
*DLOAD
TEST, GRAV,5000.,0.,1.,0.
REFER,GRAV,5000.,0.,1.,0.
*END STEP
**
**
*STEP,NLGEOM
 STEP #5: REMOVE TEST ELEMENT
*STATIC
.25,1.0
*MODEL CHANGE,REMOVE,TYPE=ELEMENT
TEST,
*END STEP
**
**
*STEP,NLGEOM
 STEP #6: ADD TEST ELEMENTS BACK IN WITH-STRAIN
*STATIC
.25,1.0
*MODEL CHANGE,ADD=WITHSTRAIN,TYPE=ELEMENT
TEST,
*END STEP
**
**
*STEP,NLGEOM
 STEP #7: APPLY THERMAL LOAD
*STATIC
.25,1.0
*BOUNDARY,OP=NEW,FIXED
TIP,1
LEFT,1,6
*TEMPERATURE
NALL,100.,20.
*EL PRINT,POSITION=CENTROID,FREQ=999
1,
S,THE,E,TEMP
*OUTPUT,FIELD,FREQ=999
*ELEMENT OUTPUT,ELSET=TEST
S,THE,E,TEMP
*END STEP
**
**
*STEP,NLGEOM
 STEP #8: FREQUENCY EXTRACTION
*FREQUENCY, eigensolver=lanczos, SIM=NO
72, , 
*BOUNDARY,OP=NEW
LEFT,1,6
*EL PRINT,FREQ=0
*NODE PRINT,FREQ=0
*OUTPUT,FIELD,FREQ=0
*END STEP
**
**
*STEP,PERTURBATION
 STEP #9: MODAL DYNAMIC
*MODAL DYNAMIC
.0004,.004
*CLOAD
TIP,2,150.
*NODE PRINT,NSET=LEFT,FREQ=999
RF, 
*EL PRINT,POSITION=CENTROID,FREQ=999
1,
S,E
*OUTPUT,FIELD,FREQ=999
*NODE OUTPUT,NSET=LEFT
RF, 
*OUTPUT,HISTORY,FREQ=999
*NODE OUTPUT,NSET=LEFT
RF, 
*OUTPUT,FIELD,FREQ=999
*ELEMENT OUTPUT,ELSET=TEST
S
*OUTPUT,FIELD,FREQ=999
*ELEMENT OUTPUT,ELSET=TEST,POSITION=CENTROID
E,S
*END STEP
**
**
*STEP,PERTURBATION
 STEP #10: STEADY-STATE DYNAMICS, MODAL BASED
*STEADY STATE DYNAMICS,INTERVAL=RANGE
10., 100., 5
*CLOAD
TIP,2,150.
*NODE PRINT,NSET=LEFT,FREQ=999
RF, 
*EL PRINT,POSITION=CENTROID,FREQ=999
1,
S,E
*OUTPUT,FIELD,FREQ=999
*NODE OUTPUT,NSET=LEFT
RF, 
*OUTPUT,HISTORY,FREQ=999
*NODE OUTPUT,NSET=LEFT
RF, 
*END STEP
**
**
*STEP,PERTURBATION
 STEP #11: STEADY-STATE DYNAMICS, DIRECT
*STEADY STATE DYNAMICS,DIRECT
10., 100., 5
*BOUNDARY,OP=NEW
LEFT,1,6
*CLOAD
TIP,2,150.
*NODE PRINT,NSET=LEFT,FREQ=999
RF, 
*EL PRINT,POSITION=CENTROID,FREQ=999
1,
S,E
*OUTPUT,FIELD,FREQ=999
*NODE OUTPUT,NSET=LEFT
RF, 
*OUTPUT,HISTORY,FREQ=999
*NODE OUTPUT,NSET=LEFT
RF, 
*OUTPUT,FIELD,FREQ=999
*ELEMENT OUTPUT,ELSET=TEST
S,
*OUTPUT,FIELD,FREQ=999
*ELEMENT OUTPUT,ELSET=TEST,POSITION=CENTROID
S,E
*********************
*output,field
*END STEP
**
**
*STEP,PERTURBATION
 STEP #12: STEADY-STATE DYNAMICS, SUBSPACE PROJECTION
*STEADY STATE DYNAMICS,SUBSPACE PROJECTION,INTERVAL=RANGE
10., 100., 5
*CLOAD
TIP,2,150.
*NODE PRINT,NSET=LEFT,FREQ=999
RF, 
*EL PRINT,POSITION=CENTROID,FREQ=999
1,
S,E

*END STEP
***********************************
*STEP,PERTURBATION
 STEP #13: STEADY-STATE DYNAMICS, DIRECT REALONLY
*STEADY STATE DYNAMICS,DIRECT,REALONLY
10., 100., 5
*BOUNDARY,OP=NEW
LEFT,1,6
*CLOAD
TIP,2,150.
*NODE PRINT,NSET=LEFT,FREQ=999
RF, 
*EL PRINT,POSITION=CENTROID,FREQ=999
1,
S,E
*OUTPUT,FIELD,FREQ=999
*NODE OUTPUT,NSET=LEFT
RF, 
*OUTPUT,HISTORY,FREQ=999
*NODE OUTPUT,NSET=LEFT
RF, 
*OUTPUT,FIELD,FREQ=999
*ELEMENT OUTPUT,ELSET=TEST
S,
*OUTPUT,FIELD,FREQ=999
*ELEMENT OUTPUT,ELSET=TEST,POSITION=CENTROID
S,E
*********************
*output,field
*END STEP
**
**
*STEP,PERTURBATION
 STEP #14: STEADY-STATE DYNAMICS, SUBSPACE PROJECTION, REALONLY
*STEADY STATE DYNAMICS,SUBSPACE PROJECTION,INTERVAL=RANGE,REALONLY
10., 100., 5
*CLOAD
TIP,2,150.
*NODE PRINT,NSET=LEFT,FREQ=999
RF, 
*EL PRINT,POSITION=CENTROID,FREQ=999
1,
S,E

*END STEP
**
*STEP,NLGEOM
 STEP #15: DYNAMICS
*DYNAMIC
.0004, .004
*BOUNDARY,OP=NEW
LEFT,1,6
*EL PRINT,POSITION=CENTROID,FREQ=999
1,
S,E
*NODE PRINT,NSET=LEFT,FREQ=999
RF, 
*CLOAD,AMP=RAMP
TIP,2,150.
*END STEP