我来出几个题目，看看有没有高手

怎么没人做啊？我做了下高斯热源加载圆柱的，不过有很多问题啊，根本就不实高斯加载，而变成了一个圆环加载，命令流如下，还请版主和各位大虾看看问题出在那阿？
/prep7
/units,si
*afun,deg
*set,element,0.005
*set,inradius,0.05
*set,outradius,0.06
*set,degree,90
*set,length,0.15
et,1,solid70
mp,dens,1,2640
mptemp,1,100,200,400,600,800,1000
mptemp,7,1200
mpdata,c,1,1,660.4,836,973.9,1061.7,1086.8,1103.5
mpdata,c,1,7,1128.6
mpdata,enth,1,1,0,7.583e4,2.592e5,4.619e5,6.875e5,9.528e5
mpdata,enth,1,7,1.278e6
mptemp,1,20,100,200,400,600,800
mptemp,7,1000,1200
mpdata,kxx,1,1,1.38,1.53,1.65,1.86,2.02,2.15
mpdata,kxx,1,7,2.29,2.43
csys,1
k,1,inradius,0,0
k,2,inradius,degree,0
k,3,inradius,degree,length
k,4,inradius,0,length
k,5,outradius,0,0
k,6,outradius,degree,0
k,7,outradius,degree,length
k,8,outradius,0,length
v,1,2,3,4,5,6,7,8
/view,1,1,1,1
lesize,5,,,10
lesize,6,,,90
lesize,8,,,60
mshkey,1
vatt,1,1,1
vmesh,1
save
FINISH
/SOL
kbc,1
antype,trans,new
timint,on
autos,on
outres,,all
tunif,15
*set,delt,1
*SET,qmax,10000000
*SET,r,0.005
*SET,pi,3.14159
*SET,lsize,2
*SET,v,90
*SET,tinc,lsize/v
*dim,flux2,table,44,60,44,y,z,time
*do,k,1,2,1
*do,i,1,44,1
*SET,flux2(i,0,k),(i-1)*lsize
*enddo
*do,j,1,60,1
*SET,flux2(0,j,k),(j-1)*0.005
*enddo
*do,i,1,44,1
*do,j,1,60,1
*SET,ycenter,v*(k-1)*tinc
*SET,zcenter,length/2
*afun,deg
*SET,distance,sqrt((((i-1)*lsize-ycenter)/(180)*pi*0.06)**2+((j-1)*0.005-zcenter)**2)
*if,distance,le,2*r,then
*SET,flux2(i,j,k),qmax/EXP(DISTANCE**2/(r**2))
*else
*SET,flux2(i,j,k),0
*endif
*enddo
*enddo
*SET,Y0,0
*SET,flux2(0,0,k),k*tinc !
Nsel,s,loc,x,0.06
!Nsel,r,loc,y,y0-1-1e-6,y0+1+1e-6
Sf,all,hflux,%flux2%
*SET,Y0,y0+2
Allsel,all
Time,(k-1)*tinc
Solve
*enddo
fini

是耦合问题，对曲面来说，关键是高斯热源的加载，热源在平板上服从高斯分布，在曲面上也可按高斯分布来计算，这些好办。但是如果焊缝是角焊缝，一部分是曲面，一部分是平面，这个时候移动热源的加载就有点麻烦了

应该是个耦合问题，关于平板的问题多，关于曲面的还没有见过。我也想知道

期待高手
论坛刚刚建立
应该广邀专家级高手