请教大家在激光成形中，如何加载一个恒速移动的热光源？

sunclark

首先定义一个载荷步
把热源作为高斯热源或椭球分布，把它离散化，生成表面热源或体内热源
运用*DO和*ENDDO循环，同时控制工作坐标系跟随热源一起移动，一起加载，即可实现横速移动的热源加载，我这里给你一个例子，是我作硕士论文的时候一个激光－YAG复合焊的很成功的小例子
/solu
TUNIF,20
ANTYPE,4
LOCAL,11,0,0,0,thick, , , ,1,1   
NSEL,U,LOC,X,width/2   
NSEL,R,EXT
SF,ALL,CONV,8.22,20
LOCAL,11,0,0,0,thick, , , -90,1,1
NSEL,U,LOC,X,-thick
NSEL,R,EXT
SF,ALL,CONV,8.22,20         
ALLSEL,ALL
remain=r/5
n=length/2/(r1/12)
Y1=0

*do,i,1,6
step1=length/2/n  
Y1=Y1+step1
allsel,all
LOCAL,11,1,0,Y1,thick, , , ,1,1  
NSEL,S,LOC,X,0,r1/6+ remain
NSEL,R,LOC,Z,0
NPLOT
SF,ALL,HFLUX,47151991.61*1.2
allsel,all

NSEL,S,LOC,X,r1/6-remain,2*r1/6+remain   
NSEL,R,LOC,Z,0
NPLOT
SF,ALL,HFLUX,40005744.44*1.2
allsel,all

NSEL,S,LOC,X,2*r1/6-remain,3*r1/6+remain   
NSEL,R,LOC,Z,0
NPLOT
SF,ALL,HFLUX,28798124.84*1.2
allsel,all

NSEL,S,LOC,X,3*r1/6-remain ,4*r1/6+remain  
NSEL,R,LOC,Z,0
NPLOT
SF,ALL,HFLUX,17588146.23*1.2
allsel,all

NSEL,S,LOC,X,4*r1/6-remain ,5*r1/6+remain  
NSEL,R,LOC,Z,0
NPLOT
SF,ALL,HFLUX,9113427.331*1.2
allsel,all

NSEL,S,LOC,X,5*r1/6-remain ,6*r1/6+remain  
NSEL,R,LOC,Z,0
NPLOT
SF,ALL,HFLUX,4006253.636*1.2
allsel,all

LOCAL,11,0,0,Y1,thick, , ,-90 ,1,1
nsel,s,loc,Y,-r-remain,r+remain
nsel,R,loc,X,0,-0.001+remain   
nsel,r,loc,Z,0
nplot
SF,ALL,HFLUX,200115088.1
allsel,all

nsel,s,loc,Y,-r-remain,r+remain
nsel,R,loc,X,-0.001-remain,-0.002+remain
nsel,r,loc,Z,0
nplot                       
SF,ALL,HFLUX,155849787.3
allsel,all

nsel,s,loc,Y,-r-remain,r+remain
nsel,R,loc,X,-0.002-remain,-0.003+remain
nsel,r,loc,Z,0
nplot
SF,ALL,HFLUX,121375936.4
allsel,all

nsel,s,loc,Y,-r-remain,r+remain
nsel,R,loc,X,-0.003-remain，-0.004+remain
nsel,r,loc,Z,0
nplot
SF,ALL,HFLUX,94527674.33
allsel,all

nsel,s,loc,Y,-r-remain,r+remain
nsel,R,loc,X,-0.004-remain,-0.005+remain
nsel,r,loc,Z,0
nplot
SF,ALL,HFLUX,73618226.79
allsel,all

nsel,s,loc,Y,-r-remain,r+remain
nsel,R,loc,X,-0.005-remain,-0.006+remain
nsel,r,loc,Z,0
nplot
SF,ALL,HFLUX,57333932.66
allsel,all

nsel,s,loc,Y,-r-remain,r+remain
nsel,R,loc,X,-0.006-remain,-0.007+remain
nsel,r,loc,Z,0
nplot
SF,ALL,HFLUX,44651711.65
allsel,all

TIME,(step1/v)*i
NROPT,FULL,,ON
AUTOTS,1  
NSUBST,25, , ,1
KBC,1   
lnsrch,1
TSRES,ERASE
TINTP,,,,1,,,,,
SOLVE

LOCAL,11,1,0,Y1,thick, , , ,1,1  
NSEL,S,LOC,X,0,r1+remain   
NSEL,R,LOC,Z,0
SFDELE,ALL,HFLUX
SF,ALL,CONV,8.22,20
ALLSEL,ALL

LOCAL,11,0,0,Y1,thick, , ,-90 ,1,1
nsel,s,loc,Y,-r-remain,r+remain
nsel,R,loc,X,0,-thick
nsel,r,loc,Z,0
SFDELE,ALL,HFLUX
SF,ALL,CONV,8.22,20
allsel,all
*enddo
save

sunclark

当然没有问题，QQ上聊

sunclark

你好，zhangjg，我把你给我的命令流给你加些注释，针对你最后出现的问题，你的热源可能定义错了，应该是功率除以面积，而不是周长，并相关修改如下
       另外，你既然定义单位为国际单位，就不要再重新定义你的计算单位为K，这也可能是计算中出现警告的原因之一，
       既然加载了热源，移动过后要把它卸载，所以你的结果会出现异常，你程序问题不小，建议你还是对照ansys的help文件认真的理解

/prep7  ！进入前处理
/units,si   !设置单位制为国际单位制
*set,x0,0  ！赋值，类似于C++中的“＝”，下同
*set,y0,0
*set,rb,0.002  !定义激光光斑半径
*set,Q,8000   !定义激光的功率密度
*set,PI,3.1415926
*set,v,0.05   !定义激光的移动速度
*set,tim,0
QC=Q/(PI* rb*rb)        
et,1,solid70  ! 定义单元类型
mp,dens,1,7833  ！定义密度
mp,c,1,465  ！定义比热
mp,kxx,1,54  ！定义热传导系数
block,0,0.052,0,0.025,-0.025,0   ！建长方体模型   
/view,1,1,1,1  !转变视图方向为ISO视角
esize,0.002,    ！单元大小
mshkey,0   ！网格设置
mshape,1,3D
vmesh,all   ！划分网格
nsel,s,loc,x,0,0.05   ！选择激光光斑区域
nsel,r,loc,y,0,0.005
nsel,r,loc,z,-0.01,0
nrefine,all,3      !对激光光斑作用的区域进行网格细化
save

x_center=x0+rb  !设置光斑的初始位置
y_center=y0
然后通过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
光斑的移动采用循环编制,每一次移动一个光斑的位置.在加载区域的选择上,由于圆形区域
不好选择,所以采用圆形区域的外接正方形区域来近似.
再输入第二部分命令流
/sol
antype,4,new  ！定义分析类型
trnopt,full
toffst,0 !设置采用的温度单位为K
tunif,300  !设置初始温度为300K
!---------------------------
*do,i,1,13
tim=tim+0.04
time,tim !设置载荷步作用时间
x_center=x0+(i-1)*2*rb !设置当前作用区域的中心点坐标        (每次x向前移动一个直径的距离）
y_center=y0
*set,laser(4,0,1),x_center !改变定义的函数中响应的中心点坐标值

nsel,s,loc,x,x_center-0.002,x_center+0.002
nsel,r,loc,y,y_center,y_center+0.002
nsel,r,loc,z,0

sf,all,hflux,%laser% !在选择的节点上面施加定义好的热流密度函数
allsel,all
nsubst,5
solve
*enddo
     你的求解设置也是需要重新考虑的
     祝你好运，以上只是我的一点小意见，仅供参考

gzgshang

"

block,0,0.052,0,0.025,-0.025,0 ！建长方体模型
.................................................
nsel,s,loc,x,x_center-0.002,x_center+0.002
nsel,r,loc,y,y_center,y_center+0.002
nsel,r,loc,z,0
sf,all,hflux,%laser% !在选择的节点上面施加定义好的热流密度函数
"

高斯方程是一个关于x,y的二维形式，沿x轴移动的热源y值选的为什么选在厚度当中，我不太明白，请您赐教！

trs20042004

我的毕业设计是激光打孔的 是不是也能和sunclark学习那 我的QQ8403345 希望你能加我 以后可以请教