用重启动修改材料参数的问题

dragonwen

模型初始温度为175度，初始状态无应力（Reference temp=175），然后，求解的第一步模型降温到100，在此过程中材料的CTE为20E-6；在求解的第二步继续降温到25度，在第二步中材料的CTE变为5E-6。
下图中的蓝色曲线是某一点Z方向位移的的时间历程曲线。从图中可以明显看出，在第二步开始时，出现了一个明显的突变，这个突变产生的原因是材料CTE的改变。但是同样可以看出，在程序计算第二步的时候，他并没有累加第一步的计算结果，而是仍然以175作为参考温度，根据此时材料的CTE和温度来计算物体的变形（Deformation=CTE*ΔT）。
从图中的初始变形量和最终变形量的比较可以说明这一点。图中降温到100度后材料的变形量竟然比继续降温到25度的变形量还大，这十分不合理。真实的情况应该是模型第二步的变形是在第一步变形的基础上进行，也就是说随着时间的增加，总的变形量应该逐渐增大，虽然第二步的变形量不如第一步的变形量大，就像图中的红色曲线所示的那样。

这个问题我在版内搜索过，但是都没有一个解决办法。。。再次提出。。希望高手出招。。。

[ 本帖最后由 dragonwen 于 2007-2-9 17:38 编辑 ]

dragonwen

下面是命令流文件：   

/PREP7

/VIEW,1,1,1,1   
/ANG,1  
/REP,FAST

ET,1,SOLID45
BLC4,,,100E-3,100E-3,1000E-3

MP,EX,1,160E9
MP,PRXY,1,0.3
MP,ALPX,1,20E-6

MP,EX,2,160E9
MP,PRXY,2,0.3
MP,ALPX,2,5E-6

TYPE,1
MAT,1
ESIZE,10E-3,0,
MSHAPE,0,3D
MSHKEY,1

VMESH,ALL

NSEL,S,LOC,Z,0
D,ALL,,,,,,ALL,,,,,
ALLSEL,ALL
EPLOT
SAVE

/TITLE,STEP1: CTE=20E-6
/SOLU
TREF,175
NSEL,S,,,ALL
BF,ALL,TEMP,100
ALLSEL,ALL
EPLOT

TIME,10
AUTOTS,1
NSUBST,20,40,10
OUTRES,ALL,ALL
OUTPR,BASIC,ALL

SOLVE

/TITLE,STEP2: CTE=5E-5
ESEL,S,MAT,,1
MPCHG,2,ALL
ALLSEL,ALL

ANTYPE,,REST

NSEL,S,,,ALL
BF,ALL,TEMP,25
ALLSEL,ALL
EPLOT

TIME,20
AUTOTS,1
NSUBST,20,40,10
OUTRES,ALL,ALL
OUTPR,BASIC,ALL

SOLVE


/POST26
NSOL,2,122,U,Z, UZ_2
STORE,MERGE
XVAR,1  
PLVAR,2,
/XRANGE,DEFAULT
/YRANGE,-3E-3,2E-3,1
/REPLOT