焊接

clplpsb1985

一个激光焊点的传热模拟，目的是得到其温度场。里面的参数随意设定的
SOLID70单元
感觉建模不行，想搞成1/4球的形式，但划分网格时出问题。恳请高手指点
finish
/clear,all
/prep7
/units,si
*set,rb,0.001
*set,length,0.05
*set,width,0.05
*set,height,0.05
et,1,solid70
mptemp,1,20,250,500,750,1000,1500,1700,2500
mpdata,c,1,1,460,480,530,675,670,660,780,820
mpdata,kxx,1,1,50,47,40,27,30,35,140,142
mpdata,dens,1,1,7820,7700,7610,7550,7490,7350,7300,7090
mpdata,hf,1,1,100,350,520,1000,1500,3000,3100,3500
mpdata,c,2,1,460,480,530,675,670,660,780,820
mpdata,kxx,2,1,50,47,40,27,30,35,140,142
mpdata,dens,2,1,7820,7700,7610,7550,7490,7350,7300,7090
mpdata,hf,2,1,100,350,520,1000,1500,3000,3100,3500
block,0,length,0,height,0,-width
block,length,2*length,0,height,0,-width
block,0,length,height,2*height,0,-width
block,0,length,0,height,-width,-2*width
block,length,2*length,height,2*height,0,-width
block,length,2*length,0,height,-width,-2*width
block,0,length,height,2*height,-width,-2*width
block,length,2*length,height,2*height,-width,-2*width
allsel,all
vglue,all
numcmp,all
/view,1,1,1,1
vsel,s,volu,,1
LESIZE,1,,,height/rb
LESIZE,4,,,length/rb
LESIZE,9,,,width/rb
MSHKEY,1
VATT,1,1,1
mat,1
type,1
VMESH,1
allsel,all
vsel,u,volu,,1
mat,2
lesize,38,,,10,5
lesize,39,,,10,5
lesize,29,,,10,5
lesize,31,,,10,5
lesize,35,,,10,5
lesize,34,,,10,5
vmesh,all
allsel,all
/solu
antype,trans,new
TIMINT,OFF
KBC,1
time,0.001
trnopt,full
tunif,27
KK=3
Qmax=1234567
R=0.01
MAX_TIME=1
MAX_X= length/rb+1
MAX_Y=height/rb+1
*DIM,FLUX2,TABLE,MAX_X,MAX_Y,MAX_TIME,X,Y,TIME  
k=max_time
*DO,I,1,MAX_X,1
FLUX2(I,0,K)=(I-1)*rb
*ENDDO
*DO,J,1,MAX_Y,1
FLUX2(0,J,K)=(J-1)*rb
*ENDDO
*DO,I,1,MAX_X,1
*DO,J,1,MAX_Y,1  
DISTANCE=SQRT(((I-1)*rb)**2+((J-1)*rb)**2)
*IF,DISTANCE,LE,R,THEN
FLUX2(I,J,K)=Qmax/EXP(kk*DISTANCE**2/(0.01**2))
*ELSE
FLUX2(I,J,K)=0
*ENDIF
*ENDDO
*ENDDO
FLUX2(0,0,K)=(K-1)
asel,s,,,2
nsla,r,1
sf,all,hflux,%flux2%
allsel,all
asel,s,loc,x,0
asel,a,loc,y,0
asel,a,loc,z,0
nsla,r,1
nsel,inve
sf,all,conv,27,27
allsel,all
nsubst,10,20,5
lswrite,1
Save

sfdele,all,hflux
ddele,all,temp
OUTRES,all,1  
TIMINT,ON
TINTP,,,,1
TIME,100
kbc,0
DELTIME,1,.5,10
neqit,25
AUTOTS,ON
ALLSEL,ALL
lswrite,2
lssolve,1,2