*HEADING
 CONTINUUM ELEMENTS: MODEL CHANGE VERIFICATION, ELASTIC MATERIAL
*RESTART,WRITE
*PREPRINT,MODEL=YES,HISTORY=NO
**
** TEST 1: 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=TEST1
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 2: CPE4H
**
*NODE,NSET=NALL
101, -5.0, -9.0
102, -4.0, -9.0
103, -4.0, -10.0
104, -5.0, -10.0
105, -3.0, -9.0
106, -3.0, -10.0
107, -5.0, -12.0
108, -3.5, -12.1
109, -3.5, -13.1
110, -5.0, -13.0
111, -2.0, -12.2
112, -2.0, -13.2
*NSET,NSET=LEFT
101,104,107,110
*NSET,NSET=MID
108,109
*NSET,NSET=RIGHT
111,112
*NSET,NSET=TIP
105,106,111,112
*ELEMENT,TYPE=CPE4H,ELSET=TEST2
101, 101,104,103,102
102, 102,103,106,105
103, 107,110,109,108
104, 108,109,112,111
*ELSET,ELSET=REFER
101,102
*ELSET,ELSET=TEST
103,104
**
** TEST 3: CPE4I
**
*NODE,NSET=NALL
201, -5.0, -18.0
202, -4.0, -18.0
203, -4.0, -19.0
204, -5.0, -19.0
205, -3.0, -18.0
206, -3.0, -19.0
207, -5.0, -21.0
208, -3.5, -21.1
209, -3.5, -22.1
210, -5.0, -22.0
211, -2.0, -21.2
212, -2.0, -22.2
*NSET,NSET=LEFT
201,204,207,210
*NSET,NSET=MID
208,209
*NSET,NSET=RIGHT
211,212
*NSET,NSET=TIP
205,206,211,212
*ELEMENT,TYPE=CPE4I,ELSET=TEST3
201, 201,204,203,202
202, 202,203,206,205
203, 207,210,209,208
204, 208,209,212,211
*ELSET,ELSET=REFER
201,202
*ELSET,ELSET=TEST
203,204
**
** TEST 4: CPE4R
**
*NODE,NSET=NALL
301, -5.0, -27.0
302, -4.0, -27.0
303, -4.0, -28.0
304, -5.0, -28.0
305, -3.0, -27.0
306, -3.0, -28.0
307, -5.0, -30.0
308, -3.5, -30.1
309, -3.5, -31.1
310, -5.0, -31.0
311, -2.0, -30.2
312, -2.0, -31.2
*NSET,NSET=LEFT
301,304,307,310
*NSET,NSET=MID
308,309
*NSET,NSET=RIGHT
311,312
*NSET,NSET=TIP
305,306,311,312
*ELEMENT,TYPE=CPE4R,ELSET=TEST4
301, 301,304,303,302
302, 302,303,306,305
303, 307,310,309,308
304, 308,309,312,311
*ELSET,ELSET=REFER
301,302
*ELSET,ELSET=TEST
303,304
**
** TEST 5: 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=TEST5
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 6: CPE8H
**
*NODE,NSET=NALL
501,  5.0, -9.0
502,  6.0, -9.0
503,  6.0, -10.0
504,  5.0, -10.0
505,  7.0, -9.0
506,  7.0, -10.0
507,  7.0, -9.5
508,  5.0, -9.5
509,  5.0, -12.0
510,  6.5, -12.1
511,  6.5, -13.1
512,  5.0, -13.0
513,  8.0, -12.2
514,  8.0, -13.2
515,  8.0, -12.7
516,  5.0, -12.5
*NSET,NSET=LEFT
501,504,508,509,512,516
*NSET,NSET=MID
510,511
*NSET,NSET=RIGHT
513,514,515
*NSET,NSET=TIP
505,506,507,513,514,515
*ELEMENT,TYPE=CPE8H,ELSET=TEST6
501, 501,504,506,505,508,503,507,502
502, 509,512,514,513,516,511,515,510
*ELSET,ELSET=REFER
501, 
*ELSET,ELSET=TEST
502, 
**
** TEST 7: CPE8R
**
*NODE,NSET=NALL
601,  5.0, -18.0
602,  6.0, -18.0
603,  6.0, -19.0
604,  5.0, -19.0
605,  7.0, -18.0
606,  7.0, -19.0
607,  7.0, -18.5
608,  5.0, -18.5
609,  5.0, -21.0
610,  6.5, -21.1
611,  6.5, -22.1
612,  5.0, -22.0
613,  8.0, -21.2
614,  8.0, -22.2
615,  8.0, -21.7
616,  5.0, -21.5
*NSET,NSET=LEFT
601,604,608,609,612,616
*NSET,NSET=MID
610,611
*NSET,NSET=RIGHT
613,614,615
*NSET,NSET=TIP
605,606,607,613,614,615
*ELEMENT,TYPE=CPE8R,ELSET=TEST7
601, 601,604,606,605,608,603,607,602
602, 609,612,614,613,616,611,615,610
*ELSET,ELSET=REFER
601, 
*ELSET,ELSET=TEST
602, 
**
**  REMAINDER OF INPUT DECK SAME FOR ALL TESTS
**
*ELSET,ELSET=EALL
TEST,REFER
*SOLID SECTION,MATERIAL=ELAS,ELSET=EALL
1.0, 
*MATERIAL,NAME=ELAS
*ELASTIC
1000., .25
*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
.5,1.
*MODEL CHANGE,REMOVE,TYPE=ELEMENT
TEST,
*BOUNDARY,OP=NEW
LEFT,1,2
*END STEP
**
**
*STEP
 STEP #3: ADD TEST ELEMENTS BACK IN STRAIN-FREE
*STATIC
1.,1.
*MODEL CHANGE,ADD=STRAINFREE,TYPE=ELEMENT
TEST, 
*END STEP
**
**
*STEP
 STEP #4: APPLY LOAD TO BOTH ELEMENTS 
*STATIC
0.1,1.
*DLOAD
TEST,TRVEC4,-0.25,0,-1,0
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
 STEP 11: FREQUENCY EXTRACTION FOR SUBSPACE PROJECTION
*FREQUENCY, eigensolver=lanczos, SIM=NO
124, , 
*BOUNDARY,OP=NEW
LEFT,1,2
*EL PRINT,FREQ=0
*NODE PRINT,FREQ=0
*OUTPUT,FIELD,FREQ=0
*END STEP
**
**
*STEP,PERTURBATION
 STEP #12: 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
*OUTPUT,FIELD,FREQ=999
*NODE OUTPUT
RF, 
*END STEP
**
**
*STEP,PERTURBATION
 STEP #13: STEADY-STATE DYNAMICS, SUBSPACE PROJECTION
*STEADY STATE DYNAMICS,SUBSPACE PROJECTION, INTERVAL=RANGE
1000., 10000., 5
**if 124 modes were used in the *FREQUENCY step, results would be , SIM=NO
**identical with *SSD,direct
*CLOAD
TIP,2,0.2
*NODE PRINT,FREQ=999
RF, 
*EL PRINT,FREQ=999
S,E
*END STEP
**
**
*STEP
 STEP #14: 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