*HEADING
C3D10MT  TET ELEMENTS: MODEL CHANGE VERIFICATION, ELASTIC MATERIAL
*RESTART,WRITE
*PREPRINT,MODEL=YES,HISTORY=NO

** TEST 3: C3D10
**
*NODE,NSET=NALL
101,0.,5.,0.5
102,1.,5.,1.
103,1.,5.,0.
104,1.,6.,0.5
111,.5,5.,0.75
112,1,5,0.5
113,.5,5,.25
114,.5,5.5,.5
115,1,5.5,.75
116,1,5.5,.25
105,2,5,.5
117,1.5,5,.75
118,1.5,5,.25
119,1.5,5.5,0.5
**
121,3.,5.,0.5
122,4.,5.,1.
123,4.,5.,0.
124,4.,6.03,0.5
131,3.5,5.,0.75
132,4,5,0.5
133,3.5,5,.25
134,3.5,5.515,.5
135,4,5.515,.75
136,4,5.515,.25
125,5,5,.5
137,4.5,5,.75
138,4.5,5.,.25
139,4.5,5.515,0.5
*NSET,NSET=BOT
101,102,103,105,111,112,113,117,118,121,122,123
125,131,132,133,137,138
*NSET,NSET=TIP
104,124
*NSET,NSET=MIDSIDE
134,135,136,139
*NSET,NSET=MIDSIDEW
134,135,136,139,114,115,116,119
*NSET,NSET=DEFBC
124, 
*ELEMENT,TYPE=C3D10MT,ELSET=TEST3
101,101,102,103,104,111,112,113,114,115,116
102,102,105,103,104,117,118,112,115,119,116
103,121,122,123,124,131,132,133,134,135,136
104,122,125,123,124,137,138,132,135,139,136
*ELSET,ELSET=REFER
101,102
*ELSET,ELSET=TEST
103,104
**
**
*ELSET,ELSET=EALL
TEST,REFER
*SOLID SECTION,MATERIAL=ELAS,ELSET=EALL
*MATERIAL,NAME=ELAS
*ELASTIC
1.E6, .25
*EXPANSION
1.E-6, 
*DENSITY
.2829, 
*CONDUCTIVITY
 7.872e-4, 
*INITIAL CONDITIONS,TYPE=TEMPERATURE
NALL,70.,0.
*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 
*COUPLED TEMPERATURE-DISPLACEMENT, STEADY STATE
0.25,1.
*PRINT,MODEL CHANGE=YES
*BOUNDARY,OP=NEW
BOT,PINNED
DEFBC,2,2,-0.03
DEFBC,1,1,0.0
DEFBC,3,3,0.0
MIDSIDE,2,2,-0.015
MIDSIDE,1,1,0.0
MIDSIDE,3,3,0.0
101,11,,70.
102,11,,70.
121,11,,70.
122,11,,70.
115,11,,70.
135,11,,70.
*NODE PRINT,FREQ=999
RF, 
*EL PRINT,POSITION=CENTROID,FREQ=999
S, 
E, 
*END STEP
**
**
*STEP,NLGEOM
 STEP #2: REMOVE TEST ELEMENTS
*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
*COUPLED TEMPERATURE-DISPLACEMENT, STEADY STATE
1.,1.
*MODEL CHANGE,ADD=STRAINFREE,TYPE=ELEMENT
TEST, 
*BOUNDARY,OP=NEW
BOT,PINNED
101,11,,70.
102,11,,70.
121,11,,70.
122,11,,70.
115,11,,70.
135,11,,70.
*END STEP
**
**
*STEP,NLGEOM
 STEP #4: APPLY DLOAD TYPE GRAV 
*COUPLED TEMPERATURE-DISPLACEMENT, STEADY STATE
0.25,1.
*DLOAD,OP=NEW
TEST, GRAV,150.,0.,1.,0.
REFER,GRAV,150.,0.,1.,0.
*END STEP
**
**
*STEP,NLGEOM
 STEP #5: APPLY DLOAD TYPE BX
*COUPLED TEMPERATURE-DISPLACEMENT, STEADY STATE
.25,1.0
*DLOAD, OP=NEW
TEST, BX, 150.
REFER, BX, 150. 
*END STEP
**
**
*STEP,NLGEOM
 STEP #6: APPLY DLOAD TYPE BY
*COUPLED TEMPERATURE-DISPLACEMENT, STEADY STATE
.25,1.0
*DLOAD, OP=NEW
TEST, BY, 150.
REFER,BY, 150.
*END STEP
**
**
*STEP,NLGEOM
 STEP #7: APPLY DLOAD TYPE BZ
*COUPLED TEMPERATURE-DISPLACEMENT, STEADY STATE
.25,1.0
*DLOAD, OP=NEW
TEST, BZ, 150.
REFER,BZ, 150.
*END STEP
**
**
*STEP,NLGEOM
 STEP #8: APPLY DLOAD TYPE HP1
*COUPLED TEMPERATURE-DISPLACEMENT, STEADY STATE
.25,1.0
*DLOAD, OP=NEW
TEST, HP1, 1.E3, 2., 0. 
REFER, HP1, 1.E3, 2., 0.
*END STEP
**
**
*STEP,NLGEOM
 STEP #9: APPLY DLOAD TYPE HP2
*COUPLED TEMPERATURE-DISPLACEMENT, STEADY STATE
.25,1.0
*DLOAD, OP=NEW
TEST, HP2, 1.E3, 2., 0. 
REFER, HP2, 1.E3, 2., 0.
*END STEP
**
**
*STEP,NLGEOM
 STEP #10: APPLY DLOAD TYPE HP3
*COUPLED TEMPERATURE-DISPLACEMENT, STEADY STATE
.25,1.0
*DLOAD, OP=NEW
TEST, HP3, 1.E3, 2., 0. 
REFER, HP3, 1.E3, 2., 0.
*END STEP
**
**
**
*STEP,NLGEOM
 STEP #11: APPLY DLOAD TYPE HP4
*COUPLED TEMPERATURE-DISPLACEMENT, STEADY STATE
.25,1.0
*DLOAD, OP=NEW
TEST, HP4, 1.E3, 2., 0. 
REFER, HP4, 1.E3, 2., 0.
*END STEP
**
**
*STEP,NLGEOM
 STEP #12: APPLY DLOAD TYPE P1
*COUPLED TEMPERATURE-DISPLACEMENT, STEADY STATE
.25,1.0
*DLOAD, OP=NEW
TEST, P1, 1.E3, 2., 0. 
REFER, P1, 1.E3, 2., 0.
*END STEP
**
**
*STEP,NLGEOM
 STEP #13: APPLY DLOAD TYPE P2
*COUPLED TEMPERATURE-DISPLACEMENT, STEADY STATE
.25,1.0
*DLOAD, OP=NEW
TEST, P2, 1.E3, 2., 0. 
REFER, P2, 1.E3, 2., 0.
*END STEP
**
**
*STEP,NLGEOM
 STEP #14: APPLY DLOAD TYPE P3
*COUPLED TEMPERATURE-DISPLACEMENT, STEADY STATE
.25,1.0
*DLOAD, OP=NEW
TEST, P3, 1.E3, 2., 0. 
REFER, P3, 1.E3, 2., 0.
*END STEP
**
**
**
*STEP,NLGEOM
 STEP #15: APPLY DLOAD TYPE P4
*COUPLED TEMPERATURE-DISPLACEMENT, STEADY STATE
.25,1.0
*DLOAD, OP=NEW
TEST, P4, 1.E3, 2., 0. 
REFER, P4, 1.E3, 2., 0.
*END STEP