*HEADING
: *TEMPERATURE AND *PRESSURE WITH *MASS DIFFUSION
*RESTART,WRITE
*PREPRINT,ECHO=YES,HISTORY=NO,MODEL=NO
*NODE,NSET=ALL
*NSET,NSET=FIX1
*NSET,NSET=FIX2
*NSET,NSET=FIX3
*NSET,NSET=FIX4
*ELEMENT,TYPE=DC2D4, ELSET=EALL
*SOLID SECTION,MATERIAL=A1, ELSET=EALL
*MATERIAL,NAME=A1
*DIFFUSIVITY
*KAPPA,TYPE=TEMP
*KAPPA,TYPE=PRESS
*SOLUBILITY
*PHYSICAL CONSTANTS,ABSOLUTEZ=-273
*AMPLITUDE,NAME=TEST
0.0,0.0,0.5,0.1,1.0,0.75,1.0,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
*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
*PRESSURE STRESS,AMPLITUDE=TEST
*PRESSURE STRESS,AMPLITUDE=TEST,INPUT=xpressinp.inp
*END STEP
**
**
**
*STEP
*MASS DIFFUSION,STEADY STATE
*PRESSURE STRESS,OP=NEW
*END STEP
**
**
**
*STEP
*MASS DIFFUSION,STEADY STATE
*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
*TEMPERATURE,OP=NEW
*END STEP
**
**
**
*STEP
making pressure stress for element 1 different from element 2
*MASS DIFFUSION,STEADY STATE
*PRESSURE STRESS,AMPLITUDE=TEST,OP=MOD
*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
*PRESSURE STRESS,AMPLITUDE=TEST,OP=MOD
*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