*HEADING
C3D10MT TET ELEMENTS: MODEL CHANGE VERIFICATION, ELASTIC MATERIAL
**RESTART,WRITE
*PREPRINT,MODEL=YES,HISTORY=NO
** TEST 1: C3D10MT
**
*SYSTEM
0.,5.,0.
*NODE,NSET=NALL
201,0.,5.,0.5
202,1.,5.,1.
203,1.,5.,0.
204,1.,6.,0.5
211,.5,5.,0.75
212,1,5,0.5
213,.5,5,.25
214,.5,5.5,.5
215,1,5.5,.75
216,1,5.5,.25
205,2,5,.5
217,1.5,5,.75
218,1.5,5,.25
219,1.5,5.5,0.5
**
221,3.,5.,0.5
222,4.,5.,1.
223,4.,5.,0.
224,4.,6.03,0.5
231,3.5,5.,0.75
232,4,5,0.5
233,3.5,5,.25
234,3.5,5.515,.5
235,4,5.515,.75
236,4,5.515,.25
225,5,5,.5
237,4.5,5,.75
238,4.5,5.,.25
239,4.5,5.515,0.5
*NSET,NSET=BOT
201,202,203,205,211,212,213,217,218,221,222,223
225,231,232,233,237,238
*NSET,NSET=TIP
204,224
*NSET,NSET=MIDSIDE
234,235,236,239
*NSET,NSET=MIDSIDEW
234,235,236,239,214,215,216,219
*NSET,NSET=DEFBC
224, 
*ELEMENT,TYPE=C3D10MT,ELSET=TEST5
201,201,202,203,204,211,212,213,214,215,216
202,202,205,203,204,217,218,212,215,219,216
203,221,222,223,224,231,232,233,234,235,236
204,222,225,223,224,237,238,232,235,239,236
*ELSET,ELSET=REFER
201,202
*ELSET,ELSET=TEST
203,204
**
**
**
**
*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
201,11,,70.
215,11,,70.
221,11,,70.
235,11,,70.
*NODE PRINT,FREQ=999
RF, 
*EL PRINT,POSITION=CENTROID,FREQ=999
S, 
LE, 
*OUTPUT,FIELD,FREQ=999
*NODE OUTPUT
RF, 
*ELEMENT OUTPUT
LE, 
S, 
*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
201,11,,70.
215,11,,70.
221,11,,70.
235,11,,70.
*END STEP
**
**
*STEP,NLGEOM
 STEP #4: APPLY LOAD TO BOTH ELEMENTS 
*COUPLED TEMPERATURE-DISPLACEMENT,STEADY STATE
0.25,1.
*DLOAD
TEST, GRAV,150.,0.,1.,0.
REFER,GRAV,150.,0.,1.,0.
*END STEP
**
**
*STEP,NLGEOM
 STEP #5: REMOVE TEST ELEMENT
*COUPLED TEMPERATURE-DISPLACEMENT,STEADY STATE
.25,1.0
*MODEL CHANGE,REMOVE,TYPE=ELEMENT
TEST,
*END STEP
**
**
*STEP,NLGEOM
 STEP #6: ADD TEST ELEMENTS BACK IN WITH-STRAIN
*COUPLED TEMPERATURE-DISPLACEMENT,STEADY STATE
.25,1.0
*MODEL CHANGE,ADD=WITHSTRAIN,TYPE=ELEMENT
TEST,
*END STEP
**
**
*STEP,NLGEOM
 STEP #7: APPLY THERMAL LOAD
*COUPLED TEMPERATURE-DISPLACEMENT,STEADY STATE
.25,1.0
*BOUNDARY,OP=NEW,FIXED
BOT,1,3
TIP,1,3
MIDSIDEW,1,3
*BOUNDARY, OP=NEW
NALL,11,,100.
*EL PRINT,POSITION=CENTROID,FREQ=999
S, 
THE, 
LE, 
TEMP, 
*OUTPUT,FIELD,FREQ=999
*ELEMENT OUTPUT
S, 
THE, 
LE, 
TEMP, 
*END STEP
**
**