双椭球热源的函数表达式

双椭球热源的函数表达式
～双椭球的函数表达式 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
!焊接过程开始
*DIM,%_FNCNAME%,TABLE,7,38,1,,,,%_FNCCSYS% %_FNCNAME%(0,0,1)= 0.0, -999
%_FNCNAME%(2,0,1)= 0.0
%_FNCNAME%(3,0,1)= %_FNC_C1(1)%
%_FNCNAME%(4,0,1)= %_FNC_C2(1)%
%_FNCNAME%(5,0,1)= %_FNC_C3(1)%
%_FNCNAME%(6,0,1)= %_FNC_C4(1)%
%_FNCNAME%(7,0,1)= %_FNC_C5(1)%
%_FNCNAME%(0,1,1)= 1.0, -1, 16, 3, 0, 0, 0 %_FNCNAME%(0,2,1)= 0.0, -2, 0, 6, 0, 0, -1 %_FNCNAME%(0,3,1)= 0, -3, 0, 1, -2, 3, -1
%_FNCNAME%(0,4,1)= 0.0, -1, 0, 1, -3, 3, 17 %_FNCNAME%(0,5,1)= 0.0, -2, 0, 1, -1, 3, 18 %_FNCNAME%(0,6,1)= 0.0, -1, 0, 0, 0, 0, 0
%_FNCNAME%(0,7,1)= 0.0, -3, 0, 1, 0, 0, -1 %_FNCNAME%(0,8,1)= 0.0, -4, 0, 1, -1, 2, -3 %_FNCNAME%(0,9,1)= 0.0, -1, 0, 3, 0, 0, -4
%_FNCNAME%(0,10,1)= 0.0, -3, 0, 1, -4, 3, -1 %_FNCNAME%(0,11,1)= 0.0, -1, 0, 2, 0, 0, 2 %_FNCNAME%(0,12,1)= 0.0, -4, 0, 1, 2, 17, -1
%_FNCNAME%(0,13,1)= 0.0, -1, 0, 2, 0, 0, 19 %_FNCNAME%(0,14,1)= 0.0, -5, 0, 1, 19, 17, -1 %_FNCNAME%(0,15,1)= 0.0, -1, 0, 1, -4, 4, -5
%_FNCNAME%(0,16,1)= 0.0, -4, 0, 2, 0, 0, 3 %_FNCNAME%(0,17,1)= 0.0, -5, 0, 1, 3, 17, -4 %_FNCNAME%(0,18,1)= 0.0, -4, 0, 2, 0, 0, 20
%_FNCNAME%(0,19,1)= 0.0, -6, 0, 1, 20, 17, -4 %_FNCNAME%(0,20,1)= 0.0, -4, 0, 1, -5, 4, -6 %_FNCNAME%(0,21,1)= 0.0, -5, 0, 1, -1, 1, -4
%_FNCNAME%(0,22,1)= 0.0, -1, 0, 2, 0, 0, 4 %_FNCNAME%(0,23,1)= 0.0, -4, 0, 1, 4, 17, -1 %_FNCNAME%(0,24,1)= 0.0, -1, 0, 2, 0, 0, 21
%_FNCNAME%(0,25,1)= 0.0, -6, 0, 1, 21, 17, -1 %_FNCNAME%(0,26,1)= 0.0, -1, 0, 1, -4, 4, -6 %_FNCNAME%(0,27,1)= 0.0, -4, 0, 1, -5, 1, -1
%_FNCNAME%(0,28,1)= 0.0, -1, 0, 1, -3, 3, -4 %_FNCNAME%(0,29,1)= 0.0, -1, 7, 1, -1, 0, 0 %_FNCNAME%(0,30,1)= 0.0, -3, 0, 1, -2, 3, -1
%_FNCNAME%(0,31,1)= 0.0, -1, 16, 3.14159265358979310, 0, 0, 0 %_FNCNAME%(0,32,1)= 0.0, -2, 0, 3.14159265358979310, 0, 0, -1
%_FNCNAME%(0,33,1)= 0.0, -4, 0, 1, -2, 3, -1 %_FNCNAME%(0,34,1)= 0.0, -1, 0, 1, -4, 3, 20 %_FNCNAME%(0,35,1)= 0.0, -2, 0, 1, -1, 3, 19
%_FNCNAME%(0,36,1)= 0.0, -1, 0, 1, -2, 3, 21 %_FNCNAME%(0,37,1)= 0.0, -2, 0, 1, -3, 4, -1 %_FNCNAME%(0,38,1)= 0.0, 99, 0, 1, -2, 0, 0
BF,FRONT,HGEN,%TFRONT%
CMDELE,FRONT
ALLSEL,ALL
!!!后半球节点的选取
LOCAL,12,2,0,disy,0,,,,chb/ah,bh/ah, NSEL,S,LOC,X,0,ah
NSEL,R,LOC,Y,-90,0
nsel,r,loc,y,-90,0
CM,BEHIND,NODE
!!!后半球加载
LOCAL,14,0,0,disy,0,,,,1,1, !!!后半球函数、表格定义
!!!后半球函数、表格定义
*DEL,_FNCNAME
*DEL,_FNCMTID
*DEL,_FNC_C1
*DEL,_FNC_C2
*DEL,_FNC_C3
*DEL,_FNC_C4
*DEL,_FNC_C5
*DEL,_FNCCSYS
*DEL,'TBEHIND'
*SET,_FNCNAME,'TBEHIND'
*DIM,_FNC_C1,,1
*DIM,_FNC_C2,,1
*DIM,_FNC_C3,,1
*DIM,_FNC_C4,,1
*DIM,_FNC_C5,,1
*SET,_FNC_C1(1),q
*SET,_FNC_C2(1),rb
*SET,_FNC_C3(1),ah
*SET,_FNC_C4(1),chb
*SET,_FNC_C5(1),bh
*SET,_FNCCSYS,14
*DIM,%_FNCNAME%,TABLE,7,38,1,,,,%_FNCCSYS%
%_FNCNAME%(0,0,1)= 0.0, -999 %_FNCNAME%(2,0,1)= 0.0
%_FNCNAME%(3,0,1)= %_FNC_C1(1)% %_FNCNAME%(4,0,1)= %_FNC_C2(1)%
%_FNCNAME%(5,0,1)= %_FNC_C3(1)% %_FNCNAME%(6,0,1)= %_FNC_C4(1)%
%_FNCNAME%(7,0,1)= %_FNC_C5(1)% %_FNCNAME%(0,1,1)= 1.0, -1, 16, 3, 0, 0, 0
%_FNCNAME%(0,2,1)= 0.0, -2, 0, 6, 0, 0, -1
%_FNCNAME%(0,3,1)= 0, -3, 0, 1, -2, 3, -1
%_FNCNAME%(0,4,1)= 0.0, -1, 0, 1, -3, 3, 17
%_FNCNAME%(0,5,1)= 0.0, -2, 0, 1, -1, 3, 18 %_FNCNAME%(0,6,1)= 0.0,
-1, 0, 0, 0, 0, 0 %_FNCNAME%(0,7,1)= 0.0, -3, 0, 1, 0, 0, -1
%_FNCNAME%(0,8,1)= 0.0, -4, 0, 1, -1, 2, -3 %_FNCNAME%(0,9,1)= 0.0, -1, 0, 3, 0, 0, -4 %_FNCNAME%(0,10,1)= 0.0, -3, 0, 1, -4, 3, -1
%_FNCNAME%(0,11,1)= 0.0, -1, 0, 2, 0, 0, 2 %_FNCNAME%(0,12,1)= 0.0, -4, 0, 1, 2, 17, -1 %_FNCNAME%(0,13,1)= 0.0, -1, 0, 2, 0, 0, 19
%_FNCNAME%(0,14,1)= 0.0, -5, 0, 1, 19, 17, -1 %_FNCNAME%(0,15,1)= 0.0, -1, 0, 1, -4, 4, -5 %_FNCNAME%(0,16,1)= 0.0, -4, 0, 2, 0, 0, 3
%_FNCNAME%(0,17,1)= 0.0, -5, 0, 1, 3, 17, -4 %_FNCNAME%(0,18,1)= 0.0, -4, 0, 2, 0, 0, 20 %_FNCNAME%(0,19,1)= 0.0, -6, 0, 1, 20, 17, -4
%_FNCNAME%(0,20,1)= 0.0, -4, 0, 1, -5, 4, -6 %_FNCNAME%(0,21,1)= 0.0, -5, 0, 1, -1, 1, -4 %_FNCNAME%(0,22,1)= 0.0, -1, 0, 2, 0, 0, 4
%_FNCNAME%(0,23,1)= 0.0, -4, 0, 1, 4, 17, -1 %_FNCNAME%(0,24,1)= 0.0, -1, 0, 2, 0, 0, 21 %_FNCNAME%(0,25,1)= 0.0, -6, 0, 1, 21, 17, -1
%_FNCNAME%(0,26,1)= 0.0, -1, 0, 1, -4, 4, -6 %_FNCNAME%(0,27,1)= 0.0, -4, 0, 1, -5, 1, -1 %_FNCNAME%(0,28,1)= 0.0, -1, 0, 1, -3, 3, -4
%_FNCNAME%(0,29,1)= 0.0, -1, 7, 1, -1, 0, 0 %_FNCNAME%(0,30,1)= 0.0, -3, 0, 1, -2, 3, -1 %_FNCNAME%(0,31,1)= 0.0, -1, 16, 3.14159265358979310, 0, 0, 0
%_FNCNAME%(0,32,1)= 0.0, -2, 0, 3.14159265358979310, 0, 0, -1
%_FNCNAME%(0,33,1)= 0.0, -4, 0, 1, -2, 3, -1 %_FNCNAME%(0,34,1)= 0.0, -1, 0, 1, -4, 3, 20 %_FNCNAME%(0,35,1)= 0.0, -2, 0, 1, -1, 3, 19
%_FNCNAME%(0,36,1)= 0.0, -1, 0, 1, -2, 3, 21 %_FNCNAME%(0,37,1)= 0.0, -2, 0, 1, -3, 4, -1 %_FNCNAME%(0,38,1)= 0.0, 99, 0, 1, -2, 0, 0
BF,BEHIND,HGEN,%TBEHIND%
CMDELE,BEHIND
ALLSEL,ALL
SOLVE
BFDELE,all,HGEN
*ENDDO
SAVE
!*********************************** !冷却阶段
!***********************************
time,40
nsubst,25
solve
time,640,
nsubst,20
solve
time,2500
nsubst,50
solve
outres,all,all finish
～下面是双椭球的函数式: ～6*sqrt(3)*q*rf*exp(-
3*({X}^2/ah^2+{Y}^2/chf^2+{Z}^2/bh^2))/({PI}*sqrt({PI})*chf*ah*bh) ～6*sqrt(3)*q*rb*exp(-
3*({X}^2/ah^2+{Y}^2/chb^2+{Z}^2/bh^2))/({PI}*sqrt({PI})*chb*ah*bh)。