*HEADING
 CONTINUUM ELEMENTS: MODEL CHANGE VERIFICATION, ELASTIC MATERIAL
 Coupled temperature-displacement elements,
*RESTART,WRITE
*PREPRINT,MODEL=YES,HISTORY=NO
** TEST 30: CPE4T
**
*NODE,NSET=NALL
40001,  55.0,  0.0
40002,  56.0,  0.0
40003,  56.0, -1.0
40004,  55.0, -1.0
40005,  57.0,  0.0
40006,  57.0, -1.0
40007,  55.0, -3.0
40008,  56.5, -3.1
40009,  56.5, -4.1
40010,  55.0, -4.0
40011,  58.0, -3.2
40012,  58.0, -4.2
*NSET,NSET=LEFT
40001,40004,40007,40010
*NSET,NSET=MID
40008,40009
*NSET,NSET=RIGHT
40011,40012
*NSET,NSET=TIP
40005,40006,40011,40012
*ELEMENT,TYPE=CPE4T,ELSET=TEST30
40001, 40001,40004,40003,40002
40002, 40002,40003,40006,40005
40003, 40007,40010,40009,40008
40004, 40008,40009,40012,40011
*ELSET,ELSET=REFER
40001,40002
*ELSET,ELSET=TEST
40003,40004
**
** TEST 31: CPS4T
**
*NODE,NSET=NALL
50001,  55.0, -9.0
50002,  56.0, -9.0
50003,  56.0, -10.0
50004,  55.0, -10.0
50005,  57.0, -9.0
50006,  57.0, -10.0
50007,  55.0, -12.0
50008,  56.5, -12.1
50009,  56.5, -13.1
50010,  55.0, -13.0
50011,  58.0, -12.2
50012,  58.0, -13.2
*NSET,NSET=LEFT
50001,50004,50007,50010
*NSET,NSET=MID
50008,50009
*NSET,NSET=RIGHT
50011,50012
*NSET,NSET=TIP
50005,50006,50011,50012
*ELEMENT,TYPE=CPS4T,ELSET=TEST31
50001, 50001,50004,50003,50002
50002, 50002,50003,50006,50005
50003, 50007,50010,50009,50008
50004, 50008,50009,50012,50011
*ELSET,ELSET=REFER
50001,50002
*ELSET,ELSET=TEST
50003,50004
**
** TEST 32: CAX8T
**
*NODE,NSET=NALL
60001,  55.0, -18.0
60002,  56.0, -18.0
60003,  56.0, -19.0
60004,  55.0, -19.0
60005,  57.0, -18.0
60006,  57.0, -19.0
60007,  57.0, -18.5
60008,  55.0, -18.5
60009,  55.0, -21.0
60010,  56.5, -21.1
60011,  56.5, -22.1
60012,  55.0, -22.0
60013,  58.0, -21.2
60014,  58.0, -22.2
60015,  58.0, -21.7
60016,  55.0, -21.5
*NSET,NSET=LEFT
60001,60004,60008,60009,60012,60016
*NSET,NSET=MID
60010,60011
*NSET,NSET=RIGHT
60013,60014,60015
*NSET,NSET=TIP
60005,60006,60007,60013,60014,60015
*ELEMENT,TYPE=CAX8T,ELSET=TEST32
60001, 60001,60004,60006,60005,60008,60003,60007,60002
60002, 60009,60012,60014,60013,60016,60011,60015,60010
*ELSET,ELSET=REFER
60001, 
*ELSET,ELSET=TEST
60002, 
****************************beginmods
** TEST 33: CPE4RT
**
*NODE,NSET=NALL
70001,  55.0, -27.0
70002,  56.0, -27.0
70003,  56.0, -28.0
70004,  55.0, -28.0
70005,  57.0, -27.0
70006,  57.0, -28.0
70007,  55.0, -30.0
70008,  56.5, -30.1
70009,  56.5, -31.1
70010,  55.0, -31.0
70011,  58.0, -30.2
70012,  58.0, -31.2
*NSET,NSET=LEFT
70001,70004,70007,70010
*NSET,NSET=MID
70008,70009
*NSET,NSET=RIGHT
70011,70012
*NSET,NSET=TIP
70005,70006,70011,70012
*ELEMENT,TYPE=CPE4RT,ELSET=TEST33
70001, 70001,70004,70003,70002
70002, 70002,70003,70006,70005
70003, 70007,70010,70009,70008
70004, 70008,70009,70012,70011
*ELSET,ELSET=REFER
70001,70002
*ELSET,ELSET=TEST
70003,70004
**
** TEST 34: CPS4RT
**
*NODE,NSET=NALL
80001,  55.0, -36.0
80002,  56.0, -36.0
80003,  56.0, -37.0
80004,  55.0, -37.0
80005,  57.0, -36.0
80006,  57.0, -37.0
80007,  55.0, -39.0
80008,  56.5, -39.1
80009,  56.5, -40.1
80010,  55.0, -40.0
80011,  58.0, -39.2
80012,  58.0, -40.2
*NSET,NSET=LEFT
80001,80004,80007,80010
*NSET,NSET=MID
80008,80009
*NSET,NSET=RIGHT
80011,80012
*NSET,NSET=TIP
80005,80006,80011,80012
*ELEMENT,TYPE=CPS4RT,ELSET=TEST34
80001, 80001,80004,80003,80002
80002, 80002,80003,80006,80005
80003, 80007,80010,80009,80008
80004, 80008,80009,80012,80011
*ELSET,ELSET=REFER
80001,80002
*ELSET,ELSET=TEST
80003,80004
*****************************endmods
**
*NSET, NSET=FIXEDTEMP
40001,
50001,
60001,
70001,
80001,
40009,
50009,
60009,
70009,
80009,
**  REMAINDER OF INPUT DECK SAME FOR ALL TESTS
**
*ELSET,ELSET=EALL
TEST,REFER
*SOLID SECTION,MATERIAL=ELAS,ELSET=TEST30
1.0, 
*SOLID SECTION,MATERIAL=ELAS,ELSET=TEST31
1.0, 
*SOLID SECTION,MATERIAL=ELAS,ELSET=TEST32
1.0, 
*SOLID SECTION,MATERIAL=ELAS,ELSET=TEST33
1.0, 
*SOLID SECTION,MATERIAL=ELAS,ELSET=TEST34
1.0, 
*MATERIAL,NAME=ELAS
*ELASTIC
1000., .25
*EXPANSION
1.E-3, 
*DENSITY
1.E-3, 
*CONDUCTIVITY
1.0, 
*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
*COUPLED TEMPERATURE-DISPLACEMENT,STEADY STATE
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
FIXEDTEMP, 11,,70.
*NODE PRINT,FREQ=999
RF, 
*EL PRINT,FREQ=999
S,E
*END STEP
**
**
*STEP,NLGEOM
 STEP #2: REMOVE TEST ELEMENTS
**STATIC
*COUPLED TEMPERATURE-DISPLACEMENT,STEADY STATE
1,1.
*MODEL CHANGE,REMOVE,TYPE=ELEMENT
TEST,
*END STEP
**
**
*STEP, NLGEOM
 STEP #3: ADD TEST ELEMENTS BACK IN STRAIN-FREE
**STATIC
*COUPLED TEMPERATURE-DISPLACEMENT,STEADY STATE
1.,1.
*MODEL CHANGE,ADD=STRAINFREE,TYPE=ELEMENT
TEST, 
*BOUNDARY,OP=NEW
LEFT,1,2
FIXEDTEMP, 11,,70.
*END STEP
**
**
*STEP,NLGEOM
 STEP #4: APPLY LOAD TO BOTH ELEMENTS 
**STATIC
*COUPLED TEMPERATURE-DISPLACEMENT,STEADY STATE
0.1,1.
*DLOAD
TEST,P4,-0.25
REFER,P4,-0.25
*END STEP
**
**
*STEP,NLGEOM
 STEP #5: REMOVE TEST ELEMENT
**STATIC
*COUPLED TEMPERATURE-DISPLACEMENT,STEADY STATE
.1,1.0
*MODEL CHANGE,REMOVE,TYPE=ELEMENT
TEST,
*END STEP
**
**
*STEP,NLGEOM
 STEP #6: ADD TEST ELEMENTS BACK IN WITH-STRAIN
**STATIC
*COUPLED TEMPERATURE-DISPLACEMENT,STEADY STATE
.1,1.0
*MODEL CHANGE,ADD=WITHSTRAIN,TYPE=ELEMENT
TEST,
*END STEP
**
**
*STEP,NLGEOM
 STEP #7: APPLY THERMAL LOAD
*COUPLED TEMPERATURE-DISPLACEMENT,STEADY STATE
.1,1.0
*BOUNDARY,OP=NEW,FIXED
TIP,1
LEFT,1,2
*BOUNDARY,OP=NEW
NALL,11,11,100.
*EL PRINT,FREQ=999
S,THE,TEMP
*END STEP