*HEADING
 CONTINUUM ELEMENTS: MODEL CHANGE VERIFICATION, ELASTIC MATERIAL
 WITH REBAR
*RESTART,WRITE
*PREPRINT,MODEL=YES,HISTORY=NO
**
** TEST 39: CPE4
**
*NODE,NSET=NALL
1,  -5.0,  0.0
2,  -4.0,  0.0
3,  -4.0, -1.0
4,  -5.0, -1.0
5,  -3.0,  0.0
6,  -3.0, -1.0
7,  -5.0, -3.0
8,  -3.5, -3.1
9,  -3.5, -4.1
10, -5.0, -4.0
11, -2.0, -3.2
12, -2.0, -4.2
*NSET,NSET=LEFT
1,4,7,10
*NSET,NSET=MID
8,9
*NSET,NSET=RIGHT
11,12
*NSET,NSET=TIP
5,6,11,12
*ELEMENT,TYPE=CPE4,ELSET=TEST39
1, 1,4,3,2
2, 2,3,6,5
3, 7,10,9,8
4, 8,9,12,11
*ELSET,ELSET=REFER
1,2
*ELSET,ELSET=TEST
3,4
**
** TEST 40: CPE8
**
*NODE,NSET=NALL
401,  5.0,  0.0
402,  6.0,  0.0
403,  6.0, -1.0
404,  5.0, -1.0
405,  7.0,  0.0
406,  7.0, -1.0
407,  7.0, -0.5
408,  5.0, -0.5
409,  5.0, -3.0
410,  6.5, -3.1
411,  6.5, -4.1
412,  5.0, -4.0
413,  8.0, -3.2
414,  8.0, -4.2
415,  8.0, -3.7
416,  5.0, -3.5
*NSET,NSET=LEFT
401,404,408,409,412,416
*NSET,NSET=MID
410,411
*NSET,NSET=RIGHT
413,414,415
*NSET,NSET=TIP
405,406,407,413,414,415
*ELEMENT,TYPE=CPE8,ELSET=TEST40
401, 401,404,406,405,408,403,407,402
402, 409,412,414,413,416,411,415,410
*ELSET,ELSET=REFER
401, 
*ELSET,ELSET=TEST
402, 
**
** TEST 41: CPS4
**
*NODE,NSET=NALL
1001,  15.0,  0.0
1002,  16.0,  0.0
1003,  16.0, -1.0
1004,  15.0, -1.0
1005,  17.0,  0.0
1006,  17.0, -1.0
1007,  15.0, -3.0
1008,  16.5, -3.1
1009,  16.5, -4.1
1010,  15.0, -4.0
1011,  18.0, -3.2
1012,  18.0, -4.2
*NSET,NSET=LEFT
1001,1004,1007,1010
*NSET,NSET=MID
1008,1009
*NSET,NSET=RIGHT
1011,1012
*NSET,NSET=TIP
1005,1006,1011,1012
*ELEMENT,TYPE=CPS4,ELSET=TEST41
1001, 1001,1004,1003,1002
1002, 1002,1003,1006,1005
1003, 1007,1010,1009,1008
1004, 1008,1009,1012,1011
*ELSET,ELSET=REFER
1001,1002
*ELSET,ELSET=TEST
1003,1004
**
** TEST 42: CPS8
**
*NODE,NSET=NALL
4001,  25.0,  0.0
4002,  26.0,  0.0
4003,  26.0, -1.0
4004,  25.0, -1.0
4005,  27.0,  0.0
4006,  27.0, -1.0
4007,  27.0, -0.5
4008,  25.0, -0.5
4009,  25.0, -3.0
4010,  26.5, -3.1
4011,  26.5, -4.1
4012,  25.0, -4.0
4013,  28.0, -3.2
4014,  28.0, -4.2
4015,  28.0, -3.7
4016,  25.0, -3.5
*NSET,NSET=LEFT
4001,4004,4008,4009,4012,4016
*NSET,NSET=MID
4010,4011
*NSET,NSET=RIGHT
4013,4014,4015
*NSET,NSET=TIP
4005,4006,4007,4013,4014,4015
*ELEMENT,TYPE=CPS8,ELSET=TEST42
4001, 4001,4004,4006,4005,4008,4003,4007,4002
4002, 4009,4012,4014,4013,4016,4011,4015,4010
*ELSET,ELSET=REFER
4001, 
*ELSET,ELSET=TEST
4002, 
**
** Axisymmetric Elements 
**
*Node, Nset=NALL
   6001,          35.,           0.
   6002,          36.,           0.
   6003,          36.,          -1.
   6004,          35.,          -1.
   6005,          37.,           0.
   6006,          37.,          -1.
   6007,          35.,          -3.
   6008,         36.5,         -3.1
   6009,         36.5,         -4.1
   6010,          35.,          -4.
   6011,          38.,         -3.2
   6012,          38.,         -4.2
  10001,          45.,           0.
  10002,          46.,           0.
  10003,          46.,          -1.
  10004,          45.,          -1.
  10005,          47.,           0.
  10006,          47.,          -1.
  10007,          47.,         -0.5
  10008,          45.,         -0.5
  10009,          45.,          -3.
  10010,         46.5,         -3.1
  10011,         46.5,         -4.1
  10012,          45.,          -4.
  10013,          48.,         -3.2
  10014,          48.,         -4.2
  10015,          48.,         -3.7
  10016,          45.,         -3.5
  11016,         35.4,          -1.
  11017,         35.4,           0.
  11018,         36.4,          -1.
  11019,         36.4,           0.
  12016,         45.8,          -1.
  12017,         45.8,         -0.5
  12018,         45.8,           0.
  13016,         35.6,        -4.04
  13017,         35.6,        -3.04
  13018,         37.1,        -4.14
  13019,         37.1,        -3.14
  14016,         46.2,        -4.08
  14017,         46.2,        -3.58
  14018,         46.2,        -3.08
*Element, type=CAX4
6001, 6001, 6004, 6003, 6002
6002, 6002, 6003, 6006, 6005
6003, 6007, 6010, 6009, 6008
6004, 6008, 6009, 6012, 6011
*Element, type=CAX8
10001, 10001, 10004, 10006, 10005, 10008, 10003, 10007, 10002
10002, 10009, 10012, 10014, 10013, 10016, 10011, 10015, 10010
**
*Elset,Elset=ONE
6001, 6002
*Elset,Elset=TWO
10001,
*Elset,Elset=THREE
6003, 6004
*Elset,Elset=FOUR
10002,
**
*ELEMENT,TYPE=SFMAX1,ELSET=SFMAX1-1
11002,11016,11017
11003,11018,11019
*SURFACE SECTION,ELSET=SFMAX1-1
*REBAR LAYER
REBAR1,1.000000,0.300000,0.000000,RBMAT1,-30.000000,1
*EMBEDDED ELEMENT, HOST ELSET=ONE
11002,
11003,
*ELEMENT,TYPE=SFMAX2,ELSET=SFMAX2-2
12002,12016,12017,12018
*SURFACE SECTION,ELSET=SFMAX2-2
*REBAR LAYER
REBAR1,1.000000,0.300000,0.000000,RBMAT1,-30.000000,1
*EMBEDDED ELEMENT, HOST ELSET=TWO
12002,
*ELEMENT,TYPE=SFMAX1,ELSET=SFMAX1-3
13002,13016,13017
13003,13018,13019
*SURFACE SECTION,ELSET=SFMAX1-3
*REBAR LAYER
REBAR2,1.000000,0.300000,0.000000,RBMAT2,-30.000000,1
*EMBEDDED ELEMENT, HOST ELSET=THREE
13002,
13003,
*ELEMENT,TYPE=SFMAX2,ELSET=SFMAX2-4
14002,14016,14017,14018
*SURFACE SECTION,ELSET=SFMAX2-4
*REBAR LAYER
REBAR2,1.000000,0.300000,0.000000,RBMAT2,-30.000000,1
*EMBEDDED ELEMENT, HOST ELSET=FOUR
14002,
**
**
*NSET,NSET=LEFT
6001,6004,6007,6010
*NSET,NSET=MID
6008,6009
*NSET,NSET=RIGHT
6011,6012
*NSET,NSET=TIP
6005,6006,6011,6012
*NSET,NSET=LEFT
10001,10004,10008,10009,10012,10016
*NSET,NSET=MID
10010,10011
*NSET,NSET=RIGHT
10013,10014,10015
*NSET,NSET=TIP
10005,10006,10007,10013,10014,10015
**
**
*Elset, elset=EALL_CAX
 6001,  6002,  6003,  6004, 10001, 10002
*Solid Section, elset=EALL_CAX, material=ELAS
1.,
**
**  REMAINDER OF INPUT DECK SAME FOR ALL TESTS
**
*ELSET,ELSET=EALL
TEST,REFER
*SOLID SECTION,MATERIAL=ELAS,ELSET=EALL
1.0, 
*ELSET,ELSET=RTEST
TEST,6003,6004,10002,SFMAX1-3,SFMAX2-4
*MATERIAL,NAME=ELAS
*ELASTIC
10.e6, .25
*EXPANSION
1.E-6, 
*DENSITY
1.E-4, 
*INITIAL CONDITIONS,TYPE=TEMPERATURE
NALL,70.
*AMPLITUDE,NAME=RAMP
0.,0.,  1.,1.
*REBAR,ELEMENT=CONTINUUM,MATERIAL=RBMAT1,NAME=REBAR1
REFER,1.0,0.3,-30,0.4,1
*REBAR,ELEMENT=CONTINUUM,MATERIAL=RBMAT2,NAME=REBAR2
TEST,1.0,0.3,-30,0.4,1
*MATERIAL,NAME=RBMAT1
*ELASTIC
30.E6, 
*MATERIAL,NAME=RBMAT2
*ELASTIC
30.E6, 
**
**
*FILE FORMAT,ZERO INCREMENT
*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,2
RIGHT,1,1,-1.0
RIGHT,2,2, 0.2
MID,1,1,-0.5
MID,2,2, 0.1
*OUTPUT,FIELD,FREQ=999
*ELEMENT OUTPUT
 S,E
*ELEMENT OUTPUT, REBAR
 S,E
*NODE PRINT,FREQ=999
RF, 
*EL PRINT,FREQ=999
S,E
*EL PRINT,REBAR,FREQ=999
S,E
*NODE FILE,FREQ=999
RF, 
*EL FILE,FREQ=999
E, S
*EL FILE, REBAR, FREQ=999
E, S
*END STEP
**
**
*STEP
 STEP #2: REMOVE TEST ELEMENTS
*STATIC
1.,1.
*MODEL CHANGE,REMOVE,TYPE=ELEMENT
RTEST,
*END STEP
**
**
*STEP
 STEP #3: ADD TEST ELEMENTS BACK IN STRAIN-FREE
*STATIC
1.,1.
*MODEL CHANGE,ADD=STRAINFREE,TYPE=ELEMENT
RTEST, 
*BOUNDARY,OP=NEW
LEFT,1,2
*EL PRINT,REBAR,FREQ=999
S,E
*END STEP
**
**
*STEP
 STEP #4: APPLY LOAD TO BOTH ELEMENTS 
*STATIC
0.1,1.
*DLOAD
TEST,P4,-0.25
REFER,P4,-0.25
6003,P4,-0.25
6004,P4,-0.25
6001,P4,-0.25
6002,P4,-0.25
10001,P4,-0.25
10002,P4,-0.25
*EL PRINT,REBAR,FREQ=999
S,E
*EL PRINT,FREQ=999
S,E
*END STEP
**
**
*STEP
 STEP #5: REMOVE TEST ELEMENT
*STATIC
.1,1.0
*MODEL CHANGE,REMOVE,TYPE=ELEMENT
RTEST,
*END STEP
**
**
*STEP
 STEP #6: ADD TEST ELEMENTS BACK IN WITH-STRAIN
*STATIC
.1,1.0
*MODEL CHANGE,ADD=WITHSTRAIN,TYPE=ELEMENT
RTEST,
*EL PRINT,REBAR,FREQ=999
S,E
*EL PRINT,FREQ=999
S,E
*END STEP
**
**
*STEP
 STEP #7: APPLY THERMAL LOAD
*STATIC
.1,1.0
*BOUNDARY,OP=NEW,FIXED
TIP,1
LEFT,1,2
*TEMPERATURE
NALL,100.
*EL PRINT,FREQ=999
S,THE,TEMP
*EL PRINT,REBAR,FREQ=999
S,THE,TEMP
*EL FILE,FREQ=999
S,THE,TEMP
*EL FILE,REBAR,FREQ=999
S,E
*ENERGY FILE,FREQ=999
*END STEP
**
**
*STEP
 STEP #8: FREQUENCY EXTRACTION
*FREQUENCY, eigensolver=lanczos, SIM=NO
32, , 
*BOUNDARY,OP=NEW
LEFT,1,2
*EL PRINT,FREQ=0
*NODE PRINT,FREQ=0
*MODAL FILE
*OUTPUT,FIELD,VAR=PRESELECT
*OUTPUT,HISTORY
*END STEP
**
**
*STEP,PERTURBATION
 STEP #9: MODAL DYNAMIC
*MODAL DYNAMIC
.0004,.004
*CLOAD
TIP,2,0.2
*NODE PRINT,FREQ=999
RF, 
*EL PRINT,FREQ=999
S,E
*EL PRINT,REBAR,FREQ=999
S,E
*NODE FILE,FREQ=999
RF, 
*EL FILE,FREQ=999
E,S
*EL FILE,REBAR,FREQ=999
E,S
*END STEP
**
**
*STEP,PERTURBATION
 STEP #10: STEADY-STATE DYNAMICS, MODAL BASED
*STEADY STATE DYNAMICS,INTERVAL=RANGE
1000., 10000., 5
*CLOAD
TIP,2,0.2
*NODE PRINT,FREQ=999
RF, 
*EL PRINT,FREQ=999
S,E
*EL PRINT,REBAR,FREQ=999
S,E
*NODE FILE,FREQ=999
RF, 
*EL FILE,FREQ=0
*END STEP
**
**
*STEP,PERTURBATION
 STEP #11: STEADY-STATE DYNAMICS, DIRECT
*STEADY STATE DYNAMICS,DIRECT
1000., 10000., 5
*BOUNDARY,OP=NEW
LEFT,1,2
*CLOAD
TIP,2,0.2
*NODE PRINT,FREQ=999
RF, 
*EL PRINT,FREQ=999
S,E
*EL PRINT,REBAR,FREQ=999
S,E
*NODE FILE,FREQ=999
RF, 
*EL FILE,FREQ=0
*output,field
*END STEP
**
**
*STEP
 STEP #12: DYNAMICS
*DYNAMIC,DIRECT
.0004, .004
*BOUNDARY,OP=NEW
LEFT,1,2
*EL PRINT,FREQ=999
S,E
*EL PRINT,REBAR,FREQ=999
S,E
*NODE PRINT,FREQ=999
RF, 
*CLOAD,AMP=RAMP
TIP,2,0.2
*ENERGY FILE,FREQ=999
*END STEP