为什么我的高斯热源加不上去啊?

yxb345

我要在面A7(从节点205开始向下移动)上施加移动的高斯热源,为什么计算完成后温度云图总是20度?





!工作环境
fini
/cle
/filn,welding          !工作名
!进入前处理
/prep7
/units,si
*afun,rad
et,1,shell57
et,2,surf152
keyopt,2,8,4
TOFFST,273                  

!定义模型的各几何参数
boxb=0.4    !箱形梁宽度
boxh=0.7    !箱形梁高度
Ib=0.4        !工字形斜撑宽度
Ih=0.4      !工字形斜撑腹板高度
boxl=6.0    !箱形梁长度
Il=3.95     !工字形斜撑长度
pi=3.1415926
Ith=pi/4      !工字形斜撑角度
Ibt=0.021        !工字形斜撑翼缘厚度
Iht=0.021        !工字形斜撑腹板厚度
boxbt=0.02    !箱形梁宽度方向厚度
boxht=0.02    !箱形梁高度方向厚度
Jc=0.6
weldb=0.02  
ns1=0.1     !网格尺寸
ns2=0.05
ns3=0.005
!设置材料参数
r,1,0.021                                                        !实常数
r,2,0.021                                                           !焊条
MPTEMP,1,20,100,300,500,600,650                    
MP,DENS,1,2.8E3                           
MP,C,1,770.045,1.671,-0.000152
MP,KXX,1,121
                                                    !双线性随动硬化弹塑性本构关系


!!!!!!!!!!!!!!!箱形梁关键点!!!!!!!
k,1,0,0,0
k,2,0,-boxh,0
k,3,-boxl/4,0,0
k,4,-boxl/4,-boxh,0
k,5,-boxl/6,0,0
k,6,-boxl/6,-boxh,0
k,7,-Jc*cos(Ith)-(Ib+0.2)*sin(Ith)/2,0,0
k,8,-Jc*cos(Ith)-(Ib+0.2)*sin(Ith)/2,-boxh,0
k,9,-boxl/2+Il*cos(Ith)-Ib*sin(Ith)/2,0,0
k,10,-boxl/2+Il*cos(Ith)-Ib*sin(Ith)/2,-boxh,0
!!!!!!工字形关键点!!!!!!!!!!!
k,11,-boxl/5+Ib*sin(Ith)/2,boxl*tan(Ith)/5+Ib*cos(Ith)/2,0
k,12,-boxl/5,boxl*tan(Ith)/5,0
k,13,-boxl/5-Ib*sin(Ith)/2,boxl*tan(Ith)/5-Ib*cos(Ith)/2,0
k,14,0.25*(-boxl/2)+Ib*sin(Ith)/2,0.25*boxl*tan(Ith)/2+Ib*cos(Ith)/2,0
k,15,0.25*(-boxl/2),0.25*(boxl*tan(Ith)/2),0
k,16,0.25*(-boxl/2)-Ib*sin(Ith)/2,0.25*boxl*tan(Ith)/2-Ib*cos(Ith)/2,0
k,17,-Jc*cos(Ith)+Ib*sin(Ith)/2,Jc*sin(Ith)+Ib*cos(Ith)/2,0
k,18,-Jc*cos(Ith),Jc*sin(Ith),0
k,19,-Jc*cos(Ith)-Ib*sin(Ith)/2,Jc*sin(Ith)-Ib*cos(Ith)/2,0
k,20,-boxl/2+Il*cos(Ith)+Ib*sin(Ith)/2,boxl*tan(Ith)/2+Ib*cos(Ith)/2-Il*sin(Ith),0
k,21,-boxl/2+Il*cos(Ith),boxl*tan(Ith)/2-Il*sin(Ith),0
k,22,-boxl/2+Il*cos(Ith)-Ib*sin(Ith)/2,boxl*tan(Ith)/2-Ib*cos(Ith)/2-Il*sin(Ith),0
!!!!!!!节点板关键点!!!!!!!!!!!
k,23,-boxl/2+Il*cos(Ith)-(Ib/2+Weldb)*sin(Ith),boxl*tan(Ith)/2-(Ib/2+Weldb)*cos(Ith)-Il*sin(Ith),0
k,24,0,(Ib+0.2)*cos(Ith)/2+Jc*sin(Ith),0
k,25,-Jc*cos(Ith)+(Ib+0.2)*sin(Ith)/2,Jc*sin(Ith)+(Ib+0.2)*cos(Ith)/2,0
k,26,-Jc*cos(Ith)+(Ib/2+Weldb)*sin(Ith),Jc*sin(Ith)+(Ib/2+Weldb)*cos(Ith),0
k,27,-Jc*cos(Ith)-(Ib/2+Weldb)*sin(Ith),Jc*sin(Ith)-(Ib/2+Weldb)*cos(Ith),0
k,28,-Jc*cos(Ith)-(Ib+0.2)*sin(Ith)/2,Jc*sin(Ith)-(Ib+0.2)*cos(Ith)/2,0
k,29,-Jc*cos(Ith)-(Ib+0.2)*sin(Ith)/2,boxl*tan(Ith)/2-(Ib/2+Weldb)*cos(Ith)-Il*sin(Ith),0
k,30,-Jc*cos(Ith)-(Ib+0.2)*sin(Ith)/2,boxl*tan(Ith)/2-(Ib/2+Weldb)*cos(Ith)-(Il+Weldb)*sin(Ith),0
k,31,0,boxl*tan(Ith)/2+(Ib/2+Weldb)*cos(Ith)-Il*sin(Ith),0
k,32,0,boxl*tan(Ith)/2+Ib*cos(Ith)/2-Il*sin(Ith),0
k,33,0,boxl*tan(Ith)/2+Ib*cos(Ith)/2-(Il+Weldb)*sin(Ith),0
k,34,0,boxl*tan(Ith)/2-Ib*cos(Ith)/2-(Il+Weldb)*sin(Ith),0
k,35,0,boxl*tan(Ith)/2-(Ib/2+Weldb)*cos(Ith)-(Il+Weldb)*sin(Ith),0
k,36,-boxl/2+Il*cos(Ith)+(Ib/2+Weldb)*sin(Ith),boxl*tan(Ith)/2+(Ib/2+Weldb)*cos(Ith)-Il*sin(Ith),0
k,37,-boxl/2+Il*cos(Ith)+Ib*sin(Ith)/2+Weldb/cos(Ith),boxl*tan(Ith)/2+Ib*cos(Ith)/2-Il*sin(Ith),0
k,38,-boxl/2+(Il+Weldb)*cos(Ith)+Ib*sin(Ith)/2,boxl*tan(Ith)/2+Ib*cos(Ith)/2-(Il+Weldb)*sin(Ith),0
k,39,-boxl/2+(Il+Weldb)*cos(Ith)-Ib*sin(Ith)/2,boxl*tan(Ith)/2-Ib*cos(Ith)/2-(Il+Weldb)*sin(Ith),0
k,40,-boxl/2+(Il+Weldb)*cos(Ith)-(Ib/2+Weldb)*sin(Ith),boxl*tan(Ith)/2-(Ib/2+Weldb)*cos(Ith)-(Il+Weldb)*sin(Ith),0
k,41,-boxl/2+(Il+Weldb)*cos(Ith),boxl*tan(Ith)/2-(Il+Weldb)*sin(Ith),0
k,42,0,boxl*tan(Ith)/2-(Il+Weldb)*sin(Ith),0
!!!网格密度为ns1的面
a,5,6,4,3
a,11,12,15,14
a,12,13,16,15
!!!网格密度为ns2的面
a,7,8,6,5
a,14,17,18,15
a,15,16,19,18
!!!网格密度为ns3的面
!!!!!!焊条!!!!
a,19,27,23,22
a,22,23,40,39
a,22,21,41,39
a,21,20,38,41
a,20,36,37,38
a,17,26,36,20

!!!!!!以下为自由网格面!!!!!!!
!!!!箱形梁竖面!!!!
a,1,2,10,9
a,9,10,8,7
a,17,18,21,20
a,19,18,21,22


!!!!!!节点板!!!!!
a,24,25,26,36,31
a,36,31,32,37
a,37,32,33,38
a,33,38,41,42
a,41,42,34,39
a,34,35,40,39
a,40,35,1,9
a,40,30,7,9
a,29,30,40,23
a,27,28,29,23
!划分网格1-3
*do,i1,1,3
alls
asel,s,,,i1
aesize,i1,ns1
type,1
mat,1
mshape,0,2d
mshkey,1
amesh,all
*enddo
!划分网格4-6
*do,i1,4,6
alls
asel,s,,,i1
aesize,i1,ns2
type,1
mat,1
mshape,0,2d
mshkey,1
amesh,all
*enddo
!划分网格7-12!!!!!        !!焊条分网
*do,i1,7,12
alls
asel,s,,,i1
lsla,s,all
lesize,all,ns3
type,1
mat,1
mshape,0,2d
mshkey,1
amesh,all
*enddo
!划分网格13-26!!!!!      
*do,i1,13,26
alls
asel,s,,,i1
!aesize,i1,ns2
type,1
mat,1
mshape,1,2d
mshkey,0
amesh,all
*enddo
nkpt,,all
nummrg,all

local,11,1,-boxl/2+(Il+Weldb)*cos(Ith),boxl*tan(Ith)/2-(Il+Weldb)*sin(Ith),0
csys,11
alls
agen,2,all,,,,-45,,,,1
csys,0
dsys,0
!205  -0.40279     0.52015      0.0000
!271  -0.40279     0.21279      0.0000
!519  -0.41279     0.20279      0.0000
!943   0.27214E-01 0.20279      0.0000
!939   0.17214E-01 0.21279      0.0000
!959   0.17214E-01 0.52015      0.0000
!*******************
I=100 !焊接电流
U=25 !焊接电弧电压
effcient=0.85 !焊接热效率η
r=0.005 !电弧有效加热半径
Qm=3*I*U*effcient/(pi*r**2) !最大热流高斯密度
kk=0.85 !能量集中系数
Qmax=Qm*kk   !!!!!!!!!!!!!!      
ly=0.3073593697   !!!!由KDIST得到        
lx=Weldb
ns3=0.005    !焊条单元网格长度
v=0.010      !焊接速度
tm=ly/v      !焊接总时间
tinc=ns3/v   !单步时间长度
max_time=1+tm/tinc
max_x=1+lx/ns3
max_y=1+ly/ns3

type,2
esurf,0

!以下为输入高斯移动热源
*dim,flux2,table,max_y,max_x,max_time,y,x,time   
*do,k,1,max_time,1  
  *do,i,1,max_y,1
   flux2(i,0,k)=(i-1)*ns3
  *enddo   
  *do,j,1,max_x,1   
   flux2(0,j,k)=(j-1)*ns3
  *enddo     
  *do,i,1,max_y,1  
    *do,j,1,max_x,1  
     xcenter=-0.40279   !节点205的X坐标
     ycenter=0.52015-v*(k-1)*tinc
     distance=sqrt(((i-1)*ns3-ycenter)**2+((j-1)*ns3-xcenter)**2)  
    *if,distance,le,r,then
        flux2(i,j,k)=1e05*Qmax/exp(3*distance**2/(r**2))      !高斯移动热源公式
    *else
        flux2(i,j,k)=0
    *endif  
    *enddo
  *enddo
  flux2(0,0,k)=(k-1)*tinc
*enddo
!!!!!!!!!!!!!高斯移动热源完毕!!!!!!!!!!!
/solu
/solu
antype,trans
tintp,0.005,,,1,,,   !设置向后积分
kbc,1
timint,off
autots,on
ESEL,S,TYPE,,2
SFE,ALL,,CONV,,100
SFE,ALL,,CONV,2,20    !表面效应单元施加对流系数
ALLSEL,ALL
tref,20
time,1e-6
d,all,temp,20
NSUBST,1
solve               !进行稳态求解，得到初始条件
ddele,all,temp

*do,t,tinc,tm,tinc
/solu
antype,,rest
time,t
deltim,tinc,tinc,tinc            
kbc,1
outres,all,all
esel,s,type,,1
nsle,s
sf,all,hflux,%flux2%        !
autots,off
                 
solcontrol,on            
timint,on
allsel,all
solve
finish
*enddo