关于移动热源引起的热屈曲问题求助

wangml

各位斑竹大侠，问题整了好久也没有结果，只好到这里求助了。希望能得到帮助。
问题描述：沿直线移动的热源扫描金属薄板，由于局部受热受到周围的约束，会产生很大的压应力，造成薄板受热区域屈曲，并最终沿扫描线弯曲。我想用ansys对薄板进行热屈曲分析，不知道可行不？下面是移动热源扫描薄板的相关代码，需要定义了laser,convect函数才能运行。由于薄板热屈曲跟热源的扫描过程有关系，不知道该怎么进行屈曲分析，请各位斑竹给看一下，

/com
/CWD,'D:\123'
/prep7
/units,si !设置单位制为国际单位制

width = 0.04
length  = 0.04
height = 0.00084

n=1 !定义扫描的行数和列数
x0=width/(n+1)
y0=length/(n+1)
rb=0.001*2 !定义定义等离子电弧半径
Q=800 !定义等离子电弧的功率
PI=3.1415926
v=0.004 !定义等离子电弧的移动速度

x_center=-2*rb !设置等离子电弧的初始位置
y_center=y0

QC=0.25*Q/(PI*rb*rb)

h = 0
et,1,solid70
mp,dens,1,7900
mp,c,1,512

mp,kxx,1,16.3
mp,kyy,1,16.3
mp,kzz,1,16.3

block,0,width,0,length,-height,0 !创建一个长方体作为半无限大平板的近似

allsel
lsel,s,loc,y,0
lsel,r,loc,z,0
LESIZE,all,rb
allsel
lsel,s,loc,x,0
lsel,r,loc,z,0  
LESIZE,all,rb
allsel
lsel,s,loc,x,0
lsel,r,loc,y,0  
LESIZE,all,height/3

MSHAPE,0,3D
MSHKEY,1
VSEL, , , ,1
VMESH,all  
allsel

!然后通过gui方式调用函数编辑器定义要加载的热流密度函数,函数的形式为
!QC*exp(-2*(({X}-x_center)^2+({Y}-y_center)^2)/(rb^2)),保存后退出,利用函数加载器加载这个函数,在弹出的对话框中定义Table Parameter Name为
!laser,然后将函数式中的参数分别定义为QC=QC,x_center=x_center,
!y_center=y_center,rb=rb

fini
/sol
antype,trans,new
trnopt,full
toffst,273            !设置采用的温度单位为摄氏度.
tunif,27              !设置初始温度为27摄氏度
tref,27
!---------------------------

nsel,s,loc,x,x_center-rb,x_center+rb
nsel,r,loc,y,y_center-rb,y_center+rb
nsel,r,loc,z,0
sf,all,hflux,%laser%  !在选择的节点上面施加定义好的热流密度函数
allsel,all

nsubst,5,20,

physics,write,thermal
physics,clear

fini
/units,si
/prep7
et,1,45

mptemp
mptemp,1,27,500,1000,1500,2000
mpdata,ex,1,1,2.0e11,1.7e11,0.9e11,0.2e11,0.01e11

mptemp
mptemp,1,27,500,1000,1500,2000
mpdata,prxy,1,1,0.3,0.3,0.3,0.3,0.3

!mp,alpx,1,1.9e-5
mptemp
mptemp,1,27,100,500,1000,1500,2000
mpdata,alpx,1,1,1.6e-5,1.66e-5,1.8e-5,2e-5,2.2e-5,2.2e-5

tb,bkin,1,6,,,
tbtemp,27
tbdata,1,200e6,0.8e10
tbtemp,500
tbdata,1,110e6,0.5e10
tbtemp,1000
tbdata,1,60e6,0.3e10
tbtemp,1500
tbdata,1,1e6,0.06e10
tbtemp,2000
tbdata,1,0.1e6,0.01e10

antype,static,new
pstres,1
SSTIF,ON
eqslv,saprse

toffst,273                !设置采用的温度单位为摄氏度.
tunif,27                  !设置初始温度为27摄氏度
tref,27

nsel,s,loc,x,width/2
dsym,symm,x
nsel,all

nsel,s,loc,y,length/2
dsym,symm,y
allsel

nsel,s,loc,x,width/2
nsel,r,loc,y,length/2
nsel,r,loc,z,-height
d,all,uz,0
allsel
physics,write,struct
physics,clear

j=1
x_center = -2*rb
y_center = y0
*set,laser(4,0,1),x_center
*set,laser(5,0,1),y_center
*set,convect(4,0,1),x_center
*set,convect(5,0,1),y_center
tim = rb/v

*do,i,1,width/rb+2
!*do,i,1,3
fini
/filename,temp3
physics,read,thermal
/sol
*if,i,eq,1,then
antype,trans,new                  
*else
/prep7
antype,trans,rest     
*endif

x_center=x_center+rb         !设置当前作用区域的中心点坐标

*set,laser(4,0,1),x_center   !改变定义的函数中响应的中心点坐标值
nsel,s,loc,x,x_center-rb,x_center+rb
nsel,r,loc,y,y_center-rb,y_center+rb
nsel,r,loc,z,0
sf,all,hflux,%laser%         !在选择的节点上面施加定义好的热流密度函数
allsel
/solu
*set,convect(4,0,1),x_center !改变定义的函数中响应的中心点坐标值
nsel,s,loc,x,x_center-4*rb,x_center+4*rb
nsel,r,loc,y,y_center-4*rb,y_center+4*rb
nsel,r,loc,z,-height
sf,all,conv,%convect%,27     !在选择的节点上面施加定义好的热流密度函数
allsel
time,tim                     !设置载荷步作用时间
rescontrol,,none
solve
save
sfdele,all,hflux
sfdele,all,conv
allsel,all
fini

physics,read,struct
fini
/filename,stru3
/assign,esav,stru3,esav                                                        
/assign,emat,stru3,emat
/sol
*if,i,eq,1,then
antype,static,new                  
*else
parsav,all                   ! Save parameters for convergence check
resume                       ! Resume DB - to return original node positions
parresume
/prep7
antype,static,rest  
fini
*endif
/solu
ldread,temp,,,tim,,temp3,rth !选择温度场计算结果作为结构分析的载荷
time,tim                     !设置载荷步作用时间
rescontrol,,none
solve
save
fini

/filename,buck3  
                                 
/solu
antype,buckle,new
PSTRES,on
bucopt,lanb,1,0
mxpand,1,0,0,yes
time,tim
solve
fini

physics,read,struct
/assign,esav,buck3,esav                                                        
/assign,emat,buck3,emat
/solu
upcoord,0.001,off
OUTRES,ALL,ALL
NLGEOM,1
PSTRES,1
SSTIF,ON
/sol
*if,i,eq,1,then
antype,static,new            
*else
parsav,all                   ! Save parameters for convergence check
resume                       ! Resume DB - to return original node positions
parresume
/prep7
antype,static,rest
fini
*endif
/solu
ldread,temp,,,tim,,temp3,rth !选择温度场计算结果作为结构分析的载荷
time,tim                     !设置载荷步作用时间
rescontrol,,none
solve
save
fini
tim = tim+rb/v
NLGEOM,0
/assign,esav                 !redirect to initial database
/assign,emat
*enddo

wangml

好冷清，问了这么久都没有人回答！

wangml

斑竹你在那里？？？我需要你和大家的帮助！！！！

iomega

我不明白的是为什么大家都在 solve循环中改变节点上的热流密度来实现移动热源？为什么不在函数定义中用时间变量？比如：

QC*exp(-2*(({X}-x_center- V*{Time})^2+({Y}-y_center)^2)/(rb^2)) 来实现热源在x方向以V速度移动。

liujin2000161

关键不是总在一条直线上运动！比如走螺旋线形等！

wangml

有道理，之前我到没有考虑过。

wangml

有没有用ansys做移动热源引起的薄板热屈曲分析的，给提供各联系方式我有问题，已经半年了没有解决！maoluwang@sina.com

ytswq

俺也是做类似研究的，不过是菜鸟，希望高手指点！QQ：343054045

csuwumin

我也在做类似的课题，好难哦

csuwumin

好方法，支持一下