*HEADING
: *TEMPERATURE AND *PRESSURE WITH *MASS DIFFUSION
*RESTART,WRITE
*PREPRINT,ECHO=YES,HISTORY=NO,MODEL=NO
*NODE,NSET=ALL
 1, 0.
 2, 7.
 3, 7., 7.
 4, 0., 7.
 11, 10.
 12, 17.
 13, 17., 7.
 14, 10., 7.
*NSET,NSET=FIX1
 1,2
 11,12
*NSET,NSET=FIX2
 2,3
 12,13
*NSET,NSET=FIX3
 3,4
 13,14
*NSET,NSET=FIX4
 1,4
 11,14
*ELEMENT,TYPE=DC2D4, ELSET=EALL
 1,1,2,3,4
 2,11,12,13,14
*SOLID SECTION,MATERIAL=A1, ELSET=EALL
*MATERIAL,NAME=A1
*DIFFUSIVITY
 3.77E-5, 
*KAPPA,TYPE=TEMP
-273, ,0.
-283, ,10.
-293, ,20.
*KAPPA,TYPE=PRESS
10., 
*SOLUBILITY
 1., 
*PHYSICAL CONSTANTS,ABSOLUTEZ=-273
*AMPLITUDE,NAME=TEST
0.0,0.0,0.5,0.1,1.0,0.75,1.0001,1.0
*INITIAL CONDITION,TYPE=TEMPERATURE,FILE=xpressic.odb,STEP=1,INC=0
*INITIAL CONDITION,TYPE=PRESSURE,FILE=xpressic,STEP=2,INC=6
*STEP,AMPLITUDE=RAMP
*MASS DIFFUSION,STEADY STATE
0.1,1.0
*FILE FORMAT,ZERO INCREMENT,ASCII
*BOUNDARY
3,11,11,100.0
13,11,11,100.0
*EL FILE
conc, 
isol, 
esol, 
*EL FILE, POSITION=AVERAGED AT NODES
conc, 
isol, 
esol, 
*EL FILE
SOL, 
*END STEP
**
**
**
*STEP,AMPLITUDE=RAMP
*MASS DIFFUSION,STEADY STATE
0.1,1.0
*PRESSURE STRESS,AMPLITUDE=TEST
1,-1
2,0
3,-1
4,-2
*PRESSURE STRESS,AMPLITUDE=TEST,INPUT=xpressinp.inp
*END STEP
**
**
**
*STEP
*MASS DIFFUSION,STEADY STATE
0.1,1.0
*PRESSURE STRESS,OP=NEW
*END STEP
**
**
**
*STEP
*MASS DIFFUSION,STEADY STATE
0.1,1.0
*TEMPERATURE,AMPLITUDE=TEST
1,10
2,20
3,10
4,0
11,10
12,20
13,10
14,0
*END STEP
**
**
**
*STEP
*MASS DIFFUSION,STEADY STATE
0.1,1.0
*TEMPERATURE,OP=NEW
*END STEP
**
**
**
*STEP
making pressure stress for element 1 different from element 2
*MASS DIFFUSION,STEADY STATE
0.1,1.0
*PRESSURE STRESS,AMPLITUDE=TEST,OP=MOD
1,-2
2,0
3,-2
4,-4
*PRESSURE STRESS,AMPLITUDE=TEST,input=xpressinp.inp
*TEMPERATURE,AMPLITUDE=TEST
1,10
2,20
3,10
4,0
11,10
12,20
13,10
14,0
*END STEP
*STEP
Pressure stress on element 1 is made equal to that in element 2.
The results for both elements will be identical at the end of
this step.
*MASS DIFFUSION,STEADY STATE
0.1,1.0
*PRESSURE STRESS,AMPLITUDE=TEST,OP=MOD
1,-1
2,0
3,-1
4,-2
*PRESSURE STRESS,AMPLITUDE=TEST,input=xpressinp.inp
*TEMPERATURE,AMPLITUDE=TEST
1,10
2,20
3,10
4,0
11,10
12,20
13,10
14,0
*END STEP