*HEADING
COUPLED TEMP-DISP ANALYSIS TO DRIVE MASS DIFFUSION
*** ALL OUTPUT REQUESTS ARE REQUIRED FOR SUBSEQUENT MASS
*** DIFFUSION ANALYSIS
*PREPRINT,ECHO=NO,MODEL=NO,HISTORY=NO
*NODE, NSET=A
*NODE, NSET=B
*NODE, NSET=C
*NODE, NSET=D
*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
*NSET,NSET=BOTTOM,GENERATE
*NSET,NSET=TEMP,GENERATE
** GENERATE A NODE SET FOR OBSERVATION AT X=7.5cm
*NSET,NSET=RPOINT,GENERATE
*NSET,NSET=LPOINT,GENERATE
*ELEMENT,TYPE=CPE8T,ELSET=SC2D
*ELGEN, ELSET=SC2D
*ELSET,ELSET=LPOINT,GENERATE
*ELSET,ELSET=RPOINT,GENERATE
*ELSET,ELSET=SINVOUT,GENERATE
*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
TEMP,273.0D-0
**
**
*STEP,INC=100,AMPLITUDE=RAMP
*COUPLED TEMPERATURE-DISPLACEMENT,DELTMX=20
0.2D0,2.0
**
*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
**
** 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,NSET=TEMP
NT,
*OUTPUT,FIELD
*NODE OUTPUT
POR,
*OUTPUT,FIELD
*ELEMENT OUTPUT
FLVEL,
*OUTPUT,FIELD
*NODE OUTPUT,NSET=TEMP
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
100.0D0,100.0
**
*END STEP