热-应力问题（各位大哥，版主帮帮忙！）

wbc4817

我是模拟一长方体铝合金材料在水中淬火过程，就是把均匀温度的工件投入一定温度的水中，差不多冷却下来时，看工件的残余应力分布情况。开始我选用solid5直接法做，结果不对，后来看到帖子说是这个单元不支持塑性变形。接着我用间接法做，热solid70单元，转为结构solid45，定义了材料随温度变化的屈服强度，把温度场的结果文件作为载荷施加给体，但为什么读取应力值的时候，每一时刻的结果都一样的啊？好像就是最后一个时刻的结（但我温度场每一时刻的结果都保存下来了），读取等效应力值比较小，plastic equivalent stress和材料在这一温度的屈服强度差不多。但我需要的是随时间变化的应力值。
各位大侠，版主们帮帮小弟，是不是ansys这样的三维模型无法做啊？还是我哪边出错了呀？非常感谢！


下面是一个二维的命令流，问题和三维的一样，应力值不变化。还有几个警告，望大哥们能帮我一下，问题到底出在哪里？谢谢！
警告：*** WARNING ***                         CP =      27.000   TIME= 10:48:16
The ETCHG command has changed some element types.  Please examine not   
only the KEYOPTs, but also the real constants, material properties,     
boundary conditions, and loadings, and then modify as needed.           
*** WARNING ***                         CP =      27.938   TIME= 10:49:20
The current solution may have been produced using different model or   
boundary condition data than is currently stored.  POST1 results may   
be erroneous unless you perform a new solution using the stored data.   
*** WARNING ***                         CP =      29.203   TIME= 10:50:08
No constraints have been defined using the D command.                    


命令流：
/BATCH  
/COM,ANSYS RELEASE 10.0    UP20050718       10:43:17    07/21/2008
/input,menust,tmp,'',,,,,,,,,,,,,,,,1   
/GRA,POWER
/GST,ON
/PLO,INFO,3
/GRO,CURL,ON
/CPLANE,1   
/REPLOT,RESIZE  
WPSTYLE,,,,,,,,0
/FILNAME,7.21,0
/PREP7  
!*  
ET,1,PLANE77
!*  
!*  
MPTEMP,,,,,,,,  
MPTEMP,1,0  
MPDATA,KXX,1,,130   
MPTEMP,,,,,,,,  
MPTEMP,1,0  
MPDATA,C,1,,841.91  
MPTEMP,,,,,,,,  
MPTEMP,1,0  
MPDATA,DENS,1,,2810
RECTNG,0,0.08,0,0.04,   
LPLOT   
FLST,5,2,4,ORDE,2   
FITEM,5,1   
FITEM,5,3   
CM,_Y,LINE  
LSEL, , , ,P51X
CM,_Y1,LINE
CMSEL,,_Y   
!*  
LESIZE,_Y1, , ,40, , , , ,1
!*  
FLST,5,2,4,ORDE,2   
FITEM,5,2   
FITEM,5,4   
CM,_Y,LINE  
LSEL, , , ,P51X
CM,_Y1,LINE
CMSEL,,_Y   
!*  
LESIZE,_Y1, , ,20, , , , ,1
!*  
CM,_Y,AREA  
ASEL, , , ,       1
CM,_Y1,AREA
CHKMSH,'AREA'   
CMSEL,S,_Y  
!*  
MSHKEY,1
AMESH,_Y1   
MSHKEY,0
!*  
CMDELE,_Y   
CMDELE,_Y1  
CMDELE,_Y2  
!*  
EPLOT   
FINISH  
/SOL
!*  
ANTYPE,4
!*  
TRNOPT,FULL
LUMPM,0
!*  
DELTIM,1,0,0
OUTRES,ERASE
OUTRES,ALL,ALL  
AUTOTS,1
TIME,120
!*  
TUNIF,470,  
ALLSEL,ALL  
FLST,5,4,4,ORDE,2   
FITEM,5,1   
FITEM,5,-4  
LSEL,S, , ,P51X
FLST,2,4,4,ORDE,2   
FITEM,2,1   
FITEM,2,-4  
  
!*  
SFL,P51X,CONV,1000, ,20,
!*  
ALLSEL,ALL  
ALLSEL,ALL  
/STATUS,SOLU
SOLVE   
FINISH  
/POST1  
SET,FIRST   
!*  
/EFACET,1   
PLNSOL, TEMP,, 0
SET,LIST,999
SET,,, ,,, ,60  
!*  
/EFACET,1   
PLNSOL, TEMP,, 0
SET,LAST
!*  
/EFACET,1   
PLNSOL, TEMP,, 0
FINISH  
/PREP7  
ETCHG,TTS   
!*  
MPTEMP,,,,,,,,  
MPTEMP,1,0  
MPDATA,EX,1,,7.1e10
MPDATA,PRXY,1,,0.33
MPTEMP,,,,,,,,  
MPTEMP,1,0  
UIMP,1,REFT,,,  
MPDATA,ALPX,1,,6.8e-5   
FINISH  
/SOL
!*  
!*  
NLGEOM,0
SSTIF,0
NROPT,AUTO, ,OFF
EQSLV,ICCG, ,0,
PRECISION,0
MSAVE,0
TOFFST,0,   
!*  
FINISH  
/POST1  
FINISH  
/SOL
TREF,20,
LDREAD,TEMP,,, , ,'7.21','rth',' '  
ALLSEL,ALL  
/STATUS,SOLU
SOLVE   
FINISH  
/POST1  
!*  
/EFACET,1   
PLNSOL, S,EQV, 0,1.0
SET,FIRST   
!*  
/EFACET,1   
PLNSOL, S,EQV, 0,1.0

[ 本帖最后由 wbc4817 于 2008-7-28 10:15 编辑 ]

下天雄

帖命令流……

hjma123

间接法要对应好时间点，否则难以得到正确结果。

woodcopper

你每一时刻的温度值保留下来不等于读进去了，仔细检查你的命令流使用是否争取，时间节点是否准确

chichest

我也和你一样的问题，大侠们请指点一下怎么才能对应好时间点？

wbc4817

帖子重新编辑了一下

[ 本帖最后由 wbc4817 于 2008-7-28 10:14 编辑 ]

woodcopper

其实呢你的问题描述已经很清楚了，但是不知道你的操作是如何的，所以我们也不能确定问题在什么地方也不能给你指出问题在哪里，如果你不是很介意，可以把命令流贴出来给大家看看，共同找到问题的处在什么地方，否则这样的求助真的是让想帮助的你的兄弟们也无从下手

wbc4817

我已经重新编辑过帖子了，有命令流，兄弟你帮我看看，可以吗？谢谢！

woodcopper

看了你的命令流，我计算过了，确实不变
不过我更改了一下LDREAD,TEMP,,, , ,'7.21','rth',' '  你这个命令
为LDREAD,TEMP,,, 1, ,'7.21','rth',' '  和LDREAD,TEMP,,, 100, ,'7.21','rth',' '  分别计算得出的结论是变化的，现在遇到的问题是每次计算只带入其中一个时间点的温度值进行进算，所以要循环的增加时间节点，然后输出每次时间节点下计算出来的结果，自己处理所要得到的值，我现在是这样理解的

所以那位兄弟解决一下这个问题

woodcopper

对了怎么贴图片了？

woodcopper

不知道用循环的命令应该怎么解决这个问题

sweetboy9063

这个问题搞明白了，能学不少东西啊。