*HEADING
 5 DOF SHELL ELEMENTS: MODEL CHANGE VERIFICATION, ELASTIC MATERIAL
*PREPRINT,MODEL=YES,HISTORY=NO
***********************
** TEST 1: S4R5
***********************
*NODE,NSET=NALL
101, 0.0, 1.,0.0
102, 1.0, 1.,0.0
103, 2.0, 1.,0.0
104, 0.0, 1.,1.0
105, 1.0, 1.,1.0
106, 2.0, 1.,1.0
111, 0.0, 0.,0.0
112, 1.5, 0.,0.0
113, 3.0, 0.,0.0
114, 0.0, 0.,1.0
115, 1.5, 0.,1.0
116, 3.0, 0.,1.0
*NSET,NSET=LEFT
101,111,104,114
*NSET,NSET=MID
112,115
*NSET,NSET=RIGHT
113,116
*NSET,NSET=TIP
103,113,106,116
*NSET,NSET=SIDE
102,103,105,106,112,113,115,116
*ELEMENT,TYPE=S4R5,ELSET=TEST1
101,101,102,105,104
102,102,103,106,105
111,111,112,115,114
112,112,113,116,115
*ELSET,ELSET=REFER
101,102
*ELSET,ELSET=TEST
111,112
***********************
** TEST 3: STRI65
***********************
*NODE,NSET=NALL
301, 10., 4.,0.0
302, 11., 4.,0.0
303, 12., 4.,0.0
304, 10.0, 4.,1.0
305, 11.0, 4.,1.0
306, 12.0, 4.,1.0
307, 10.5, 4.,0.0
308, 11.5, 4.,0.0
309, 10.0, 4.,0.5
310, 11.0, 4.,0.5
311, 12.0, 4.,0.5
312, 10.5, 4.,1.0
313, 11.5, 4.,1.0
314, 10.5, 4.,0.5
315, 11.5, 4.,0.5
321, 10.0, 3.,0.0
322, 11.5, 3.,0.0
323, 13.0, 3.,0.0
324, 10.0, 3.,1.0
325, 11.5, 3.,1.0
326, 13.0, 3.,1.0
327, 10.75,3.,0.0
328, 12.25,3.,0.0
329, 10.0, 3.,0.5
330, 11.5, 3.,0.5
331, 13.0, 3.,0.5
332, 10.75,3.,1.0
333, 12.25,3.,1.0
334, 10.75,3.,0.5
335, 12.25,3.,0.5
*NSET,NSET=LEFT
301,321,304,309,324,329
*NSET,NSET=MID
322,330,325
*NSET,NSET=MIDR
328,333,335
*NSET,NSET=MIDL
327,332,334
*NSET,NSET=RIGHT
323,331,326
*NSET,NSET=TIP
303,323,311,306,331,326
*NSET,NSET=SIDE
302,303,305,306,307,308,312
313,322,323,325,326,327,328,332,333
*ELEMENT,TYPE=STRI65,ELSET=TEST3
301, 301,302,304,307,314,309
302, 302,305,304,310,312,314
303, 302,303,305,308,315,310
304, 303,306,305,311,313,315
311, 321,322,324,327,334,329
312, 322,325,324,330,332,334
313, 322,323,325,328,335,330
314, 323,326,325,331,333,335
*ELSET,ELSET=REFER
301,302,303,304
*ELSET,ELSET=TEST
311,312,313,314
*********************
** TEST 4: S8R5
*********************
*NODE,NSET=NALL
401, 0.0, 4.,0.0
402, 1.0, 4.,0.0
403, 2.0, 4.,0.0
404, 0.0, 4.,1.0
405, 1.0, 4.,1.0
406, 2.0, 4.,1.0
407, 0.5, 4.,0.0
408, 1.5, 4.,0.0
409, 0.0, 4.,0.5
410, 1.0, 4.,0.5
411, 2.0, 4.,0.5
412, 0.5, 4.,1.0
413, 1.5, 4.,1.0
421, 0.0, 3.,0.0
422, 1.5, 3.,0.0
423, 3.0, 3.,0.0
424, 0.0, 3.,1.0
425, 1.5, 3.,1.0
426, 3.0, 3.,1.0
427, 0.75,3.,0.0
428, 2.25,3.,0.0
429, 0.0, 3.,0.5
430, 1.5, 3.,0.5
431, 3.0, 3.,0.5
432, 0.75,3.,1.0
433, 2.25,3.,1.0
*NSET,NSET=LEFT
401,421,404,409,424,429
*NSET,NSET=MID
422,430,425
*NSET,NSET=MIDR
428,433
*NSET,NSET=MIDL
427,432
*NSET,NSET=RIGHT
423,431,426
*NSET,NSET=TIP
403,423,411,406,431,426
*NSET,NSET=SIDE
402,403,405,406,407,408,412
413,422,423,425,426,427,428,432,433
*ELEMENT,TYPE=S8R5,ELSET=TEST4
401, 401,402,405,404,407,410,412,409
402, 402,403,406,405,408,411,413,410
411, 421,422,425,424,427,430,432,429
412, 422,423,426,425,428,431,433,430
*ELSET,ELSET=REFER
401,402
*ELSET,ELSET=TEST
411,412
***********************
** TEST 5: S9R5
***********************
*NODE,NSET=NALL
501, 5., 4.,0.0
502, 6., 4.,0.0
503, 7., 4.,0.0
504, 5.0, 4.,1.0
505, 6.0, 4.,1.0
506, 7.0, 4.,1.0
507, 5.5, 4.,0.0
508, 6.5, 4.,0.0
509, 5.0, 4.,0.5
510, 6.0, 4.,0.5
511, 7.0, 4.,0.5
512, 5.5, 4.,1.0
513, 6.5, 4.,1.0
514, 5.5, 4.,0.5
515, 6.5, 4.,0.5
521, 5.0, 3.,0.0
522, 6.5, 3.,0.0
523, 8.0, 3.,0.0
524, 5.0, 3.,1.0
525, 6.5, 3.,1.0
526, 8.0, 3.,1.0
527, 5.75,3.,0.0
528, 7.25,3.,0.0
529, 5.0, 3.,0.5
530, 6.5, 3.,0.5
531, 8.0, 3.,0.5
532, 5.75,3.,1.0
533, 7.25,3.,1.0
534, 5.75,3.,0.5
535, 7.25,3.,0.5
*NSET,NSET=LEFT
501,521,504,509,524,529
*NSET,NSET=MID
522,530,525
*NSET,NSET=MIDR
528,533,535
*NSET,NSET=MIDL
527,532,534
*NSET,NSET=RIGHT
523,531,526
*NSET,NSET=TIP
503,523,511,506,531,526
*NSET,NSET=SIDE
502,503,505,506,507,508,512
513,522,523,525,526,527,528,532,533
*ELEMENT,TYPE=S9R5,ELSET=TEST5
501, 501,502,505,504,507,510,512,509,514
502, 502,503,506,505,508,511,513,510,515
511, 521,522,525,524,527,530,532,529,534
512, 522,523,526,525,528,531,533,530,535
*ELSET,ELSET=REFER
501,502
*ELSET,ELSET=TEST
511,512
************************************************
**  REMAINDER OF INPUT DECK SAME FOR ALL TESTS
************************************************
*ELSET,ELSET=EALL
TEST,REFER
*SHELL SECTION,MATERIAL=ELAS,ELSET=EALL
.1, 
*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,INC=100,NLGEOM
 STEP #1: DISPLACE TEST SPECIMEN SO FINAL SHAPE IS
          IDENTICAL TO INITIAL SHAPE OF REFERENCE ELEMENT 
*STATIC
0.05,1.
*PRINT,MODEL CHANGE=YES
*BOUNDARY,OP=NEW
LEFT,1,6,0.0
MID,4,5,0.0
MIDL,4,5,0.0
MIDR,4,5,0.0
RIGHT,4,5,0.0
SIDE,3,5,0.0
RIGHT,1,1,-1.0
RIGHT,6,6,.001
MID,1,1,-0.5
MIDR,1,1,-.75
MIDL,1,1,-.25
*NODE PRINT,FREQ=999
RF, 
*EL PRINT,POSITION=CENTROID,FREQ=999
1,
S,E
*NODE FILE,FREQ=999
RF, 
*EL FILE,POSITION=CENTROID,FREQ=999
E, S
*OUTPUT,FIELD,FREQ=999
*ELEMENT OUTPUT
S
*OUTPUT,FIELD,FREQ=999
*ELEMENT OUTPUT,POSITION=CENTROID
E,S
MISESONLY,PRESSONLY
*NODE OUTPUT
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,6
*END STEP
**
**
*STEP,INC=100,NLGEOM
 STEP #4: APPLY LOAD TO BOTH ELEMENTS 
*STATIC
0.025,1.
*DLOAD
TEST, GRAV,50.,1.,1.,.1
REFER,GRAV,50.,1.,1.,.1
*END STEP
**
**
**
*STEP,NLGEOM
 STEP #5: APPLY THERMAL LOAD
*STATIC
.25,1.0
*BOUNDARY,OP=NEW,FIXED
LEFT,1,6
*TEMPERATURE
NALL,100.
*EL PRINT,POSITION=CENTROID,FREQ=999
1,
S,THE,
E,TEMP
*EL FILE,POSITION=CENTROID,FREQ=999
S,THE
E,TEMP
*END STEP
**
**
*STEP,NLGEOM
 STEP #6: FREQUENCY EXTRACTION
*FREQUENCY, eigensolver=lanczos, SIM=NO
30, , 
*BOUNDARY,OP=NEW
LEFT,1,6
*EL PRINT,FREQ=0
*NODE PRINT,FREQ=0
*OUTPUT,FIELD,FREQ=0
*MODAL FILE 
*END STEP
**
**
*STEP,PERTURBATION
 STEP #7: MODAL DYNAMIC
*MODAL DYNAMIC
.0075,.075
*CLOAD
TIP,2,150.
*NODE PRINT,FREQ=999
RF, 
*EL PRINT,POSITION=CENTROID,FREQ=999
1,
S,E
*NODE FILE,FREQ=999
RF, 
*EL FILE,POSITION=CENTROID,FREQ=999
E,S
*OUTPUT,FIELD,FREQ=999
*ELEMENT OUTPUT
S
*OUTPUT,FIELD,FREQ=999
*ELEMENT OUTPUT,POSITION=CENTROID
E,S
*NODE OUTPUT
RF,
*END STEP
**
**
*STEP,PERTURBATION
 STEP #8: STEADY-STATE DYNAMICS, MODAL BASED
*STEADY STATE DYNAMICS
3.,30., 3
*CLOAD
TIP,2,150.
*NODE PRINT,FREQ=999
RF, 
*EL PRINT,POSITION=CENTROID,FREQ=999
1,
S,E
*NODE FILE,FREQ=999
RF, 
*EL FILE,POSITION=CENTROID,FREQ=0
*OUTPUT,FIELD,FREQ=999
*ELEMENT OUTPUT
S
*OUTPUT,FIELD,FREQ=999
*ELEMENT OUTPUT,POSITION=CENTROID
E,S
*NODE OUTPUT
RF,
*END STEP
**
**
*STEP,PERTURBATION
 STEP #9: STEADY-STATE DYNAMICS, DIRECT
*STEADY STATE DYNAMICS,DIRECT
3., 30., 3
*BOUNDARY,OP=NEW
LEFT,1,6
*CLOAD
TIP,2,150.
*NODE PRINT,FREQ=999
RF, 
*EL PRINT,POSITION=CENTROID,FREQ=999
1,
S,E
*NODE FILE,FREQ=999
RF, 
*EL FILE,POSITION=CENTROID,FREQ=0
*END STEP
**
**
*STEP,NLGEOM
 STEP #10: DYNAMICS
*DYNAMIC
.0075,.075
*BOUNDARY,OP=NEW
LEFT,1,6
*EL PRINT,POSITION=CENTROID,FREQ=999
1,
S,E
*NODE PRINT,FREQ=999
RF, 
*DLOAD,OP=NEW
*OUTPUT,FIELD,FREQ=999
*ELEMENT OUTPUT
S
*OUTPUT,FIELD,FREQ=999
*ELEMENT OUTPUT,POSITION=CENTROID
E,S
*NODE OUTPUT
RF,
*END STEP