*HEADING
 CONTINUUM ELEMENTS: MODEL CHANGE VERIFICATION, ELASTIC MATERIAL
 WITH ORIENTATION
*RESTART,WRITE
*PREPRINT,MODEL=YES,HISTORY=NO
**
** TEST 24: CPS8
**
*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=CPS8,ELSET=TEST24
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, 
**
**  REMAINDER OF INPUT DECK SAME FOR ALL TESTS
**
*SOLID SECTION,MATERIAL=ELASR,ELSET=REFER,ORIENTATION=ELASR
1.0, 
*ORIENTATION,NAME=ELASR,SYSTEM=CYLINDRICAL
 4.5,-1.5,0,4.5,-1.5,1.0
*MATERIAL,NAME=ELASR
*ELASTIC,TYPE=ENGINEERING CONSTANTS
1000.,1100,1200,0.21,0.22,0.25,1500,1600
1700, 
*EXPANSION
1.E-3, 
*DENSITY
1.E-3, 
*SOLID SECTION,MATERIAL=ELAST,ELSET=TEST,ORIENTATION=ELAST
1.0, 
*ORIENTATION,NAME=ELAST,SYSTEM=CYLINDRICAL
4.5,-4.5,0,4.5,-4.5,1.0
*MATERIAL,NAME=ELAST
*ELASTIC,TYPE=ENGINEERING CONSTANTS
1000.,1100,1200,0.21,0.22,0.25,1500,1600
1700, 
*EXPANSION
1.E-3, 
*DENSITY
1.E-3, 
*INITIAL CONDITIONS,TYPE=TEMPERATURE
NALL,70.
*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,2
RIGHT,1,1,-1.0
RIGHT,2,2, 0.2
MID,1,1,-0.5
MID,2,2, 0.1
*NODE PRINT,FREQ=999
RF, 
*EL PRINT,FREQ=999
S,E
*OUTPUT,HISTORY,VARIABLE=PRESELECT
*OUTPUT,FIELD,FREQ=999
*ELEMENT OUTPUT
E,S
*NODE OUTPUT
U,
RF,
*END STEP
**
**
*STEP
 STEP #2: REMOVE TEST ELEMENTS
*STATIC
1.,1.
*MODEL CHANGE,REMOVE,TYPE=ELEMENT
TEST,
*END STEP
**
**
*STEP
 STEP #3: ADD TEST ELEMENTS BACK IN STRAIN-FREE
*STATIC
1.,1.
*MODEL CHANGE,ADD=STRAINFREE,TYPE=ELEMENT
TEST, 
*BOUNDARY,OP=NEW
LEFT,1,2
*END STEP
**
**
*STEP
 STEP #4: APPLY LOAD TO BOTH ELEMENTS 
*STATIC
0.1,1.
*DLOAD
TEST,P4,-0.25
REFER,P4,-0.25
*END STEP
**
**
*STEP
 STEP #5: REMOVE TEST ELEMENT
*STATIC
.1,1.0
*MODEL CHANGE,REMOVE,TYPE=ELEMENT
TEST,
*END STEP
**
**
*STEP
 STEP #6: ADD TEST ELEMENTS BACK IN WITH-STRAIN
*STATIC
.1,1.0
*MODEL CHANGE,ADD=WITHSTRAIN,TYPE=ELEMENT
TEST,
*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
*OUTPUT,HISTORY,VARIABLE=PRESELECT,FREQ=1
*OUTPUT,FIELD,FREQ=999
*ELEMENT OUTPUT
E,S,THE,TEMP
*NODE OUTPUT
U,
RF,
*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
*OUTPUT,FIELD,FREQ=0
*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
*OUTPUT,FIELD,FREQ=999
*NODE OUTPUT
RF, 
*ELEMENT OUTPUT
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
*OUTPUT,HISTORY,VARIABLE=PRESELECT
*OUTPUT,FIELD,FREQ=1
*ELEMENT OUTPUT
E,S,
*NODE OUTPUT
U,
RF,
*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
*END STEP
**
**
*STEP,PERTURBATION
 STEP #12: STEADY-STATE DYNAMICS, SUBSPACE PROJECTION
*STEADY STATE DYNAMICS,SUBSPACE PROJECTION,INTERVAL=RANGE
1000., 10000., 5
*CLOAD
TIP,2,0.2
*NODE PRINT,FREQ=999
RF, 
*EL PRINT,FREQ=999
S,E
*END STEP
**
**
*STEP
 STEP #13: DYNAMICS
*DYNAMIC
.0004, .004
*BOUNDARY,OP=NEW
LEFT,1,2
*EL PRINT,FREQ=999
S,E
*NODE PRINT,FREQ=999
RF, 
*CLOAD,AMP=RAMP
TIP,2,0.2
*END STEP