*HEADING
COUPLED TEMP-DISP ANALYSIS TO DRIVE MASS DIFFUSION       
*** ALL OUTPUT REQUESTS ARE REQUIRED FOR SUBSEQUENT MASS 
*** DIFFUSION ANALYSIS
*PREPRINT,ECHO=NO,MODEL=YES,HISTORY=NO
*NODE, NSET=A
1,0.0,-0.005
*NODE, NSET=B
161,0.1,-0.005
*NODE, NSET=C
2737,0.1,0.005
*NODE, NSET=D
2577,0.0,0.005
*NFILL,NSET=LEFT
A,D,16,161
*NFILL,NSET=RIGHT
B,C,16,161
*NFILL,NSET=PLATE
LEFT,RIGHT,160,1
*NSET, NSET=TOP,GENERATE
2577,2737
*NSET,NSET=BOTTOM,GENERATE
1,161
*NSET,NSET=TEMP,GENERATE
1,161,2
323,483,2
645,805,2
967,1127,2
1289,1449,2
1611,1771,2
1933,2093,2
2255,2415,2
2577,2737,2
** GENERATE A NODE SET FOR OBSERVATION AT X=7.5cm
*NSET,NSET=RPOINT,GENERATE
121,2697,161
*NSET,NSET=LPOINT,GENERATE
17,2593,161
*ELEMENT,TYPE=CPE8T,ELSET=SC2D
1,1,3,325,323,2,164,324,162
*ELGEN, ELSET=SC2D
1,80,2,1,8,322,80,1
*ELSET,ELSET=LPOINT,GENERATE
9,569,80
*ELSET,ELSET=RPOINT,GENERATE
61,621,80
*ELSET,ELSET=SINVOUT,GENERATE
1,240
401,640
*SOLID SECTION,ELSET=SC2D,MATERIAL=STEEL
1.0D0, 
*MATERIAL,NAME=STEEL
*ELASTIC
2e11,0.3
*CONDUCTIVITY
50., 
*PHYSICAL CONSTANTS,ABSOLUTEZ=0
**
*BOUNDARY
** ZERO X DISPLACEMENT AT LHS
LEFT,1,1
1289,2,2
**
**
*FILE FORMAT, ZERO INCREMENT
*INITIAL CONDITIONS, TYPE=TEMPERATURE
PLATE,273.0D-0
**TEMP,273.0D-0
**
**
*STEP,INC=100,AMPLITUDE=RAMP
*COUPLED TEMPERATURE-DISPLACEMENT,STEADY STATE
0.2,2.0,0.2,0.2
**
*CLOAD, OP=NEW
161,   1,  -1.5625e+03
322,1,  -5.4688e+03
483,1,  -2.3438e+03
644,1,  -3.9062e+03
805,1,  -1.5625e+03
966,1,  -2.3438e+03
1127,1,  -7.8125e+02
1288,1,  -7.8125e+02
1449,1,0
1610,1,7.8125e+02
1771,1,7.8125e+02
1932,1,2.3438e+03
2093,1,1.5625e+03
2254,1,3.9062e+03
2415,1,2.3438e+03
2576,1,5.4688e+03
2737,1,1.5625e+03
**
** TEMPERATURES ENTERED IN KELVIN
**
*BOUNDARY
TOP,11,11,4.23e2
BOTTOM,11,11,1.23E2
**
*RESTART,WRITE,FREQ=2
**
** THE FILE OUTPUT REQUESTS FOR THIS PROBLEM ARE REQUIRED
** TO DRIVE THE SUBSEQUENT MASS DIFFUSION ANALYSIS.
**
*NODE FILE
NT, 
*OUTPUT,FIELD
*NODE OUTPUT
POR,
*OUTPUT,FIELD
*ELEMENT OUTPUT
FLVEL,
*OUTPUT,FIELD
*NODE OUTPUT
NT,
*EL FILE,POSITION=AVERAGED AT NODES,ELSET=SINVOUT,FREQ=1
SINV, 
*EL PRINT,POSITION=AVERAGED AT NODES,ELSET=SINVOUT,FREQ=0
SINV, 
*NODE PRINT,FREQ=0
*END STEP
*STEP,INC=100,AMPLITUDE=RAMP
HOLD TEMP AND PRESS CONSTANT FOR MASS DIFF STEADY STATE
*COUPLED TEMPERATURE-DISPLACEMENT,STEADY STATE
100.0D0,100.0
**
*END STEP