本帖最后由 水幽寒 于 2009-3-15 13:54 编辑
只做过简单一点的研究,不知道对你是否有用处,对于梁单元好像没有办法施加温度梯度,所以研究过实体单元的温度梯度问题。
是做了一个悬臂梁,截面10*10m,长度100米
一端固结,在顶面向下6米范围内,施加从14度到0度的温度梯度,分别采用了ansys模型和midas模型进行了比较分析,结果还是比较理想的。
这是ansys的源文件
finish
/clear
/prep7
et,1,solid45
mp,ex,1,3.45e10
mp,dens,1,2400
mp,prxy,1,0.17
k,1,
k,2,10
k,3,,-10
k,4,10,-10
a,1,2,4,3
voffst,1,100
esize,2
vmesh,all
nsel,all
MP,alpx,1,1.0e-5
!MP,kxx,1,1.33
nsel,all
a=ndinqr(0,13)
*do,i,1,a
*if,abs(ny(i)),le,6,then
bf,i,temp,14/6*(6-abs(ny(i)))
*else
bf,i,temp,0
*endif
*enddo
/PBF,TEMP, ,1
/replot
nsel,s,loc,z,0
d,all,all
/solu
allsel
solve
附件中是midas的源文件
结果:单位MPa
ansys上缘0.81的压应力,下缘,1.26压应力
midas上缘0.79压应力,下缘,1.18压应力 |