*HEADING
 ELBOW ELEMENTS: MODEL CHANGE VERIFICATION
*PREPRINT,MODEL=YES,HISTORY=NO
**
** TEST 1: ELBOW31
**
*NODE,NSET=NALL
1,  0.0, 1.
2,  1.0, 1.
3,  2.0, 1.
11, 0.0, 0.0, .0
12, 1.5, -.1, .1
13, 3.0, -.3, .2
*NSET,NSET=LEFT
1,11
*NSET,NSET=MID
12, 
*NSET,NSET=RIGHT
13, 
*NSET,NSET=TIP
3,13
*ELEMENT,TYPE=ELBOW31,ELSET=TEST1
1, 1,2
2, 2,3
11,11,12
12,12,13
*ELSET,ELSET=REFER
1,2
*ELSET,ELSET=TEST
11,12
*ELSET,ELSET=EOUT
1,11
**
** TEST 2: ELBOW32
**
*NODE,NSET=NALL
101,  5.0, 1.
102,  6.0, 1.
103,  7.0, 1.
111, 5.0, 0.0, .0
112, 6.5, -.1, .1
113, 8.0, -.3, .2
*NSET,NSET=LEFT
101,111
*NSET,NSET=MID
112, 
*NSET,NSET=RIGHT
113, 
*NSET,NSET=TIP
103,113
*ELEMENT,TYPE=ELBOW32,ELSET=TEST1
101, 101,102,103
111, 111,112,113
*ELSET,ELSET=REFER
101, 
*ELSET,ELSET=TEST
111, 
*ELSET,ELSET=EOUT
101,111
**
** TEST 3: ELBOW31B
**
*NODE,NSET=NALL
201,  0.0, 4.
202,  1.0, 4.
203,  2.0, 4.
211,  0.0, 3.0, .0
212,  1.5, 2.9, .1
213,  3.0, 2.7, .2
*NSET,NSET=LEFT
201,211
*NSET,NSET=MID
212, 
*NSET,NSET=RIGHT
213, 
*NSET,NSET=TIP
203,213
*ELEMENT,TYPE=ELBOW31B,ELSET=TEST1
201, 201,202
202, 202,203
211, 211,212
212, 212,213
*ELSET,ELSET=REFER
201,202
*ELSET,ELSET=TEST
211,212
*ELSET,ELSET=EOUT
201,211
**
** TEST 4: ELBOW31C
**
*NODE,NSET=NALL
301,  5.0, 4.
302,  6.0, 4.
303,  7.0, 4.
311,  5.0, 3.0, .0
312,  6.5, 2.9, .1
313,  8.0, 2.7, .2
*NSET,NSET=LEFT
301,311
*NSET,NSET=MID
312, 
*NSET,NSET=RIGHT
313, 
*NSET,NSET=TIP
303,313
*ELEMENT,TYPE=ELBOW31C,ELSET=TEST1
301, 301,302
302, 302,303
311, 311,312
312, 312,313
*ELSET,ELSET=REFER
301,302
*ELSET,ELSET=TEST
311,312
*ELSET,ELSET=EOUT
301,311
**
**
**  REMAINDER OF INPUT DECK SAME FOR ALL TESTS
**
*ELSET,ELSET=EALL
TEST,REFER
*BEAM SECTION,SECTION=ELBOW,TEMPERATURE=GRADIENTS,MATERIAL=ELAS,ELSET=TEST1
1.,.1,10.
0.,10.,0.
3,6,1
*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.
*NORMAL
TEST,LEFT,0.,1.,0.
TEST,MID,.086,.9963,0.
TEST,RIGHT,.11723,.9931,0.
**
**
*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,7
RIGHT,1,1,-1.0
RIGHT,2,2, 0.3
RIGHT,3,3,-0.2
RIGHT,4,5, 0.0
RIGHT,6,6,.11750075
MID,1,1,-0.5
MID,2,2, 0.1
MID,3,3,-0.1
MID,4,5, 0.0
MID,6,6,.086105945
*NODE PRINT,FREQ=999
U, 
RF, 
*EL PRINT,ELSET=EOUT,FREQ=999
1,
S,E
*OUTPUT,FIELD,FREQ=999
*NODE OUTPUT
RF, 
*ELEMENT OUTPUT,POSITION=CENTROID
E, 
S, 
*OUTPUT,HISTORY,FREQ=999
*NODE OUTPUT, NSET=NALL
RF, 
*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,7
*NODE PRINT,FREQ=999
RF, 
*END STEP
**
**
*STEP,NLGEOM
 STEP #4: APPLY LOAD TO BOTH ELEMENTS 
*STATIC
0.25,1.
*DLOAD
TEST, GRAV,50.,0.,1.,0.
REFER,GRAV,50.,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,2
LEFT,1,7
*TEMPERATURE
NALL,80.,10.,0.
*EL PRINT,ELSET=EOUT,FREQ=999
1,
THE,E,
S,TEMP
*OUTPUT,FIELD,FREQ=999
*ELEMENT OUTPUT,ELSET=EOUT
S,THE,E,TEMP,SF
*OUTPUT,HISTORY
*ELEMENT OUTPUT,ELSET=EOUT
SF
*END STEP
**
**
*STEP,NLGEOM
 STEP #8: FREQUENCY EXTRACTION
*FREQUENCY, eigensolver=lanczos, SIM=NO
66, , , 
*BOUNDARY,OP=NEW
LEFT,1,7
*EL PRINT,FREQ=0
*NODE PRINT,FREQ=0
*OUTPUT,FIELD,FREQ=0
*END STEP
**
**
*STEP,PERTURBATION
 STEP #9: MODAL DYNAMIC
*MODAL DYNAMIC
.002,.01
*CLOAD
TIP,2,1.5E2
*NODE PRINT,FREQ=999
RF, 
*EL PRINT,ELSET=EOUT,FREQ=999
1,
S,E
*OUTPUT,FIELD,FREQ=999
*NODE OUTPUT
RF, 
*OUTPUT,HISTORY,FREQ=999
*NODE OUTPUT, NSET=NALL
RF, 
*OUTPUT,FIELD,FREQ=999
*ELEMENT OUTPUT,POSITION=CENTROID
E, 
S, 
*END STEP
**
**
*STEP,PERTURBATION
 STEP #10: STEADY-STATE DYNAMICS, MODAL BASED
*STEADY STATE DYNAMICS,INTERVAL=RANGE
10., 20., 5
*CLOAD
TIP,2,15.
*CLOAD,IMAG
TIP,2,2.0
*NODE PRINT,FREQ=999
RF, 
*EL PRINT,ELSET=EOUT,FREQ=999
1,
S,E
*OUTPUT,FIELD,FREQ=999
*NODE OUTPUT
RF, 
*OUTPUT,HISTORY,FREQ=999
*NODE OUTPUT, NSET=NALL
RF, 
*OUTPUT,FIELD,FREQ=999
*ELEMENT OUTPUT,ELSET=EOUT
1
S,E
*END STEP
**
**
*STEP,PERTURBATION
 STEP #11: STEADY-STATE DYNAMICS, DIRECT
*STEADY STATE DYNAMICS,DIRECT
10., 20., 5
*BOUNDARY,OP=NEW
LEFT,1,7
*CLOAD
TIP,2,15.
*NODE PRINT,FREQ=999
RF, 
*EL PRINT,ELSET=EOUT,FREQ=999
1,
S,E
*OUTPUT,FIELD,FREQ=999
*NODE OUTPUT
RF, 
*OUTPUT,HISTORY,FREQ=999
*NODE OUTPUT,NSET=NALL
RF, 
*OUTPUT,FIELD, FREQ=0
*ELEMENT OUTPUT,POSITION=CENTROID
*END STEP
**
**
*STEP,PERTURBATION
 STEP #11: STEADY-STATE DYNAMICS, SUBSPACE PROJECTION
*STEADY STATE DYNAMICS,SUBSPACE PROJECTION,INTERVAL=RANGE
10., 20., 5
*CLOAD
TIP,2,15.
*NODE PRINT,FREQ=999
RF, 
*EL PRINT,ELSET=EOUT,FREQ=999
1,
S,E
*OUTPUT,FIELD,FREQ=999
*NODE OUTPUT
RF, 
*OUTPUT,HISTORY,FREQ=999
*NODE OUTPUT,NSET=NALL
RF, 
*OUTPUT,FIELD, FREQ=0
*ELEMENT OUTPUT,POSITION=CENTROID
*END STEP
**
**
*STEP,NLGEOM
 STEP #12: DYNAMICS
*DYNAMIC
.002, .01
*BOUNDARY,OP=NEW
LEFT,1,7
*EL PRINT,ELSET=EOUT,FREQ=999
1,
S,E
*EL PRINT,ELSET=EOUT,FREQ=999
1,2,3
MISES
*NODE PRINT,FREQ=999
RF, 
*CLOAD,AMP=RAMP
TIP,2,15.
*OUTPUT,OP=ADD,FIELD,FREQ=10
*ELEMENT OUTPUT,ELSET=EOUT
MISESMAX
*END STEP