双椭球热源的函数表达式

双椭球热源的函数表达式

～双椭球的函数表达式 pi=3.141592654 finish step=0.002

/clear v=0.006

/filname,reverse-shuang tuo qiu,0 !count=ly/step

/title,hanhanahanhan count=10

/units,si !***********************************

/prep7 !定义单元类型

!***********************************

ET,1,SOLID70 !建立几何模型

keyopt,1,2,1 !***********************************

ET,2,SOLID90 wpstyle,,,,,,,,0 MPTEMP,1,20,100,200,300,400,500, block,0,0.006,0,ly,0,-lz MPTEMP,7,600,700,800,900,1000,1500

block,0,0.012,0,ly,0,-lz MPDATA,KXX,1,1,152,159,164,165,168,196,

block,0,0.024,0,ly,0,-lz MPDATA,KXX,1,7,209,162,162,162,162,162

block,0,lx,0,ly,0,-lz MPDATA,DENS,1,1,2680,2662,2646,2620,26vovlap,all 00,2580, /view,1,1,1,1 MPDATA,DENS,1,7,2436,2384,2384,2384,23!定义焊缝区单元尺寸0.002 84,2384 lesize,7,0.001

MPDATA,C,1,1,956,963,1047,1130,1224,130lesize,2,0.001 8, lesize,12,0.001 MPDATA,C,1,7,1412,1084,1084,1084,1084,1lesize,11,0.001 084

lesize,6,0.001 MPDATA,HF,1,1,8.22,11.0,13.7,23.2,33.4,46.lesize,3,0.001 8, lesize,8,0.001

MPDATA,HF,1,7,58.0,68.5,68.5,68.5,68.5,68.lesize,1,0.001 5

lesize,4,0.001

lesize,9,0.001 lx=0.08 lesize,5,0.001 ly=0.2 lesize,10,0.001

lz=0.005 !定义过渡区的单元尺寸

r=0.002 lesize,15,r

D=0.005～焊缝宽度 lesize,18,r

ah=D/2～双椭球热源参数定义 !定义远离焊缝区的单元尺寸 chf=2*ah lesize,34,1.25*0.001 chb=8*ah lesize,35,1.25*0.001 bh=0.005

lesize,47,1.25*0.001 rf=2*chf/(chf+chb) lesize,46,1.25*0.001

rb=2*chb/(chf+chb) lesize,39,4*r K=0.7 功率有效系数 lesize,42,4*r I=120 lesize,57,4*r U=24 lesize,60,4*r q0=I*U lesize,27,4*r q=q0*K !有效功率lesize,30,4*r

lesize,58,4*r !*********************************** lesize,59,4*r nsel,s,ext,

!*********************************** nsel,u,loc,x,0 !------改

!划分网格，生生有限元模型 sf,all,conv,200,293 !------改

!*********************************** allsel,all

type,1 !***********************************

mshkey,1 !*********************************** mshape,0 !双椭球热源的加载vmesh,1 !*********************************** vmesh,7 !前半球节点的选取

type,2 !*********************************** mshkey,0

LOCAL,11,2,0,disy,0.000,,,,chf/ah,bh/ah, mshape,1 NSEL,S,LOC,X,0,ah vmesh,5 NSEL,R,LOC,Y,0,90

vmesh,6 NSEL,R,LOC,Z,-90,0 numcmp,area CM,FRONT,NODE

save

LOCAL,13,0,0,disy,0.000,,,,1,1, !*********************************** !!!前半球函数、表格定义

!循环加载、求解 *DEL,_FNCNAME

!*********************************** *DEL,_FNCMTID

!施加载荷 *DEL,_FNC_C1

/solu *DEL,_FNC_C2

*do,i,0,count *DEL,_FNC_C3

disy=i*step *DEL,_FNC_C4

t=(disy+step)/v *DEL,_FNC_C5

antype,4 *DEL,_FNCCSYS

trnopt,full *DEL,'TFRONT'

tunif,293 *SET,_FNCNAME,'TFRONT' outres,basic,last *DIM,_FNC_C1,,1 time,t *DIM,_FNC_C2,,1

nropt,full,,on *DIM,_FNC_C3,,1

autots,1 *DIM,_FNC_C4,,1

nsubst,1,,,1 *DIM,_FNC_C5,,1

kbc,1 *SET,_FNC_C1(1),q

lnsrch,1 *SET,_FNC_C2(1),rf

tsres,erase *SET,_FNC_C3(1),ah

!TINTP,,,,1,,,,, *SET,_FNC_C4(1),chf

timint,1,therm *SET,_FNC_C5(1),bh

!*********************************** *SET,_FNCCSYS,13

!焊接过程开始