请教一个热应力问题

Forthuth

各位大侠好，恳请出手帮帮忙，问题简述如下：
两个大物体夹住一个小物体，两者均从175℃均匀降温到25℃，问系统的内部应力分布？具体问题如下：
fini
/clea
/prep7
ET,1,solid185
mp,ex,1,1.5e5
mp,prxy,1,0.3
mp,alpx,1,1.5e-5
mptemp,1,20,50,80,100,140,180
mpdata,ex,2,1,18800,18000,16000,14000,10000,6000
mpdata,prxy,2,1,0.3,0.3,0.3,0.3,0.3,0.3
mpdata,alpx,2,1,0.7e-5,0.7e-5,0.7e-5,3.4e-5,3.4e-5,3.4e-5
block,0,5,0,2,0,2
block,5,7,0,2,0,2
block,7,12,0,2,0,2
vglue,all
numcmp,all
vsel,s,loc,x,5,7
vatt,2,,1
allsel,all
esize,0.2
vmesh,all
vplot,all
fini
/solv
tref,175
da,5,all
da,8,all
allsel,all
然后呢，怎样设置载荷步？这问题应该很简单吧

Forthuth

没人帮我看看吗? 我重新改写一下
fini
/clea
/prep7
ET,1,solid185
mp,ex,1,1.5e5
mp,prxy,1,0.3
mp,alpx,1,1.5e-5
mptemp,1,20,40,60,80,100,120
mpdata,ex,2,1,18800,18000,16000,14000,10000,6000
mpdata,prxy,2,1,0.3,0.3,0.3,0.3,0.3,0.3
mpdata,alpx,2,1,0.7e-5,0.7e-5,0.7e-5,3.4e-5,3.4e-5,3.4e-5
block,0,5,0,2,0,2
block,5,7,0,2,0,2
block,7,12,0,2,0,2
vglue,all
numcmp,all
vsel,s,loc,x,5,7
vatt,2,,1
allsel,all
esize,0.2
vmesh,all
vplot,all
da,5,all
da,8,all
fini
/solv
tref,120
antype,4
trnopt,full
nlgeom,1
nropt,full
outres,all,all
tunif,100
time,1
solv

nsubst,2
kbc,0
*do,n,0,5,1
tunif,100-20*n
time,5*60*n+5*60
solve
*enddo
fini

/post1
.............

但是怎么算都不对,计算机跟本就不管120-40范围内材料常数,也就是计算机只直接计算最后一步,按ΔT=120-20计算,如果在循环中加入tref,120-20*n,则根本没热应力,问怎样把程序上一步计算的应力结果作为下一步计算的初始条件进行计算,难道这不是程序默认的计算方式吗?
先谢过了!

Forthuth

另外加入petres,1命令也是不管用的

Forthuth

没人吗?模型很简单,算起来也很快,不会超过半数2分钟的.

Forthuth

改了改，有人能帮我确认一下看看我的有没有什么错误吗？

fini
/clea
/filname,example-1
/prep7
ET,1,solid5
mp,ex,1,1.5e11
mp,prxy,1,0.3
mp,kxx,1,400
mp,alpx,1,1.5e-5
mp,c,1,500
mp,dens,1,5000

mptemp,1,20,50,80,100,140,180
mpdata,ex,2,1,1.88e10,1.8e10,1.6e10,1.4e10,1e10,0.6e10
mpdata,prxy,2,1,0.3,0.3,0.3,0.3,0.3,0.3
mpdata,alpx,2,1,0.7e-5,0.7e-5,0.7e-5,3.4e-5,3.4e-5,3.4e-5
mpdata,dens,2,1,2000,2000,2000,2000,2000,2000
mpdata,c,2,1,500,500,500,500,500,500
mpdata,kxx,2,1,300,300,300,300,300,300

block,0,0.05,0,0.02,0,0.02
block,0.05,0.07,0,0.02,0,0.02
block,0.07,0.12,0,0.02,0,0.02

vglue,all
numcmp,all

vsel,s,loc,x,0.05,0.07
vatt,2,,1

allsel,all
esize,0.004
vmesh,all
vplot,all

da,5,ux,0
da,5,uy,0
da,8,ux,0
da,8,uy,0
fini

/solv
antype,4
trnopt,full
nlgeom,1
nropt,full
Toffst,273
tref,175
outres,all,all

tunif,175
timint,off
time,0.01
allsel,all
solve

time,600
timint,on
nsubst,20,20,10
autots,on
kbc,0
asel,u,,,6,7,1
sfa,all,1,conv,1000,25
allsel,all
lswrite,1

kbc,0
time,620
timint,on
asel,u,,,6,7,1
sfa,all,1,conv,1,0
vsel,s,mat,,2
bfv,all,hgen,8e6
allsel,all
nsubst,20,20,10
lswrite,2

kbc,0
time,1500
timint,on
asel,u,,,6,7,1
sfa,all,1,conv,1000,25
vsel,s,mat,,2
bfv,all,hgen,8e6
allsel,all
nsubst,20,20,10
lswrite,3
lssolve,1,3,1,

/post1
SET,first
PLNSOL, TEMP,, 0,1.0
SET,NEXT
PLNSOL, TEMP,, 0,1.0

Forthuth

都摆了快一个月了就没人来看看这个问题吗？

Forthuth

结果: 找到 “热应力” 相关主题 500 个
一个个的查，找得我吐血啊！
高手在哪里？

Forthuth

我放了一个多月啊！总算是弄明白了吧
    在弹性范围内应变=εe+εt，只有弹性应变εe产生热应力，弹性应力只与最终的热变形有关所以中间温度的杨氏模量不改变最终的应力场分布，以前我看过一篇文章上的一个公式没怎么细想就认定σ=∫E(T)*α(T)dT了。
    至于那个alpx不管用是因为“alpx是割线线胀系数，由它得到的热应变只与最后的温度差有关。ANSYS里还有一个即时线胀系数（Instantaneous coefficients of thermal expansion ），CTEX（CTEY、CTEZ），……”

lovehuijin

你对于热应力分析都搞明白了么?咱们可不可以交流下。qq：375859283

Forthuth

什么叫都搞明白了？碰到不会的就想办法弄会吧，我不太上QQ的

liuzengxun1966

(1)所有物体都均匀降温，没有温度梯度和约束，应力应该为零。
(2)antype,4 是瞬态分析，您的问题可能不是瞬态问题，应该是静态antype,0
(3)您的问题总在变。
（4）传热分析应为antype，4；但是应力问题是antype，0
自己再琢磨一下吧。论坛只是大家讨论和学习的地方，没有义务可言。说话更应该客气。

Forthuth

(1)，边界上有位移束缚，所以有热应力
(2)，我用的是solid5，本来是想做热力直接耦合进行瞬态分析的，所以设成了antype,4，没有改过来，不过这似乎不影响结果，因为我没有设置瞬态分析所需的必要参数
(3)，问题在变是为了表述更清晰一些，其实质并没有什么变化
放了一个来月没人搭理，有点小小的抱怨也应该可以理解，我也确实很少上QQ(一个月合起来可能有一两个小时吧)，我自己都不知道是不是都搞明白了，所以只好反问了。我只是说实话，没什么别的意思，见谅