/PREP7
ET,1,SOLID70
MP,DENS,1,7930
MP,C,1,502
mptemp,1,24,40,60,80,110,150
mpdata,kxx,1,1,9.33,10.42,11.58,12.7,14,18
/VIEW,1,1,1,1
!*************定义基板尺寸***********
lx=0.080
lz=0.060
ly=0.005
!*************定义焊接参数***********
weld=0.012 !每道焊道宽度
wheight=0.003 !焊道的高度
LSIZE=0.002
V=0.004 !速度
tinc=LSIZE/V !每小段的时间间隔
pi=3.1415926
U=25
I=180
Q=U*I !电源有效功率
effect=0.8
f1=0.6
f2=1.4
ah=0.003
chf=2*ah
chb=8*ah
bh=0.003
Qmax11=6*effect*3**(1/2)*f1*q/(pi*pi**(1/2)*ah*bh*chf)
Qmax12=6*effect*3**(1/2)*f2*q/(pi*pi**(1/2)*ah*bh*chb)
!*****************建立模型,生成网格**************
block,0,lx,0,ly,0,lz
wpoff,0,ly,lz/2
block,0,lx,0,wheight*5,-weld/2,weld/2
wpoff,0,0,-weld/2
vsbw,all,,delete
wpoff,0,0,weld
vsbw,all,,delete
wpoff,0,-ly,-weld/2
block,-0.02,0,0,ly,-lz/2,lz/2
block,lx,lx+0.02,0,ly,-lz/2,lz/2
wpoff,0,ly,weld/2
vglue,all
FLST,5,4,4,ORDE,4
FITEM,5,14
FITEM,5,16
FITEM,5,19
FITEM,5,39
CM,_Y,LINE
LSEL, , , ,P51X
CM,_Y1,LINE
CMSEL,,_Y
!*
LESIZE,_Y1, , ,40, , , , ,1
!*
FLST,5,4,4,ORDE,4
FITEM,5,23
FITEM,5,-24
FITEM,5,69
FITEM,5,71
CM,_Y,LINE
LSEL, , , ,P51X
CM,_Y1,LINE
CMSEL,,_Y
!*
LESIZE,_Y1, , ,6, , , , ,1
!*
FLST,5,4,4,ORDE,4
FITEM,5,13
FITEM,5,15
FITEM,5,70
FITEM,5,72
CM,_Y,LINE
LSEL, , , ,P51X
CM,_Y1,LINE
CMSEL,,_Y
!*
LESIZE,_Y1, , ,10, , , , ,1
!*
MSHAPE,0,3D
MSHKEY,1
!*
CM,_Y,VOLU
VSEL, , , , 8
CM,_Y1,VOLU
CHKMSH,'VOLU'
CMSEL,S,_Y
!*
VMESH,_Y1
!*
CMDELE,_Y
CMDELE,_Y1
CMDELE,_Y2
!*
ESIZE,2e-3,0,
/DIST,1,1.08222638492,1
/REP,FAST
/DIST,1,1.08222638492,1
/REP,FAST
/DIST,1,1.08222638492,1
/REP,FAST
/DIST,1,1.08222638492,1
/REP,FAST
/DIST,1,1.08222638492,1
/REP,FAST
/DIST,1,1.08222638492,1
/REP,FAST
/DIST,1,1.08222638492,1
/REP,FAST
/DIST,1,1.08222638492,1
/REP,FAST
/DIST,1,1.08222638492,1
/REP,FAST
/DIST,1,1.08222638492,1
/REP,FAST
/DIST,1,1.08222638492,1
/REP,FAST
/DIST,1,1.08222638492,1
/REP,FAST
FLST,5,3,6,ORDE,3
FITEM,5,1
FITEM,5,7
FITEM,5,9
CM,_Y,VOLU
VSEL, , , ,P51X
CM,_Y1,VOLU
CHKMSH,'VOLU'
CMSEL,S,_Y
!*
VMESH,_Y1
!*
CMDELE,_Y
CMDELE,_Y1
CMDELE,_Y2
!*
FLST,5,2,6,ORDE,2
FITEM,5,10
FITEM,5,-11
CM,_Y,VOLU
VSEL, , , ,P51X
CM,_Y1,VOLU
CHKMSH,'VOLU'
CMSEL,S,_Y
!*
VSWEEP,_Y1
!*
CMDELE,_Y
CMDELE,_Y1
CMDELE,_Y2
!*
/UI,MESH,OFF
/solu
antype,trans
trnopt,full
nropt,full,,on
AUTOTS,on
lnsrch,on
tunif,25
!**************杀死焊缝单元********************
nsel,s,loc,y,ly,ly+wheight*5
esln,s,1
CM,weld_elem,ELEM
ekill,all
esel,s,elem,,weld_elem
nsle,s
cm,weld_node,node
ALLSEL,ALL
esel,s,live
eplot
!**********施加对流载荷*************
! 熔积层暴露空气,辐射和对流
ALLSEL,ALL
NSEL,S,LOC,y,0,0.020
NSEL,R,EXT
SF,ALL,CONV,rad,25 !空气
!**********约束杀死单元的温度约束*************
nsel,s,node,,weld_node
D,ALL,temp,298
!************定义数组维数**********
*do,j,1,5,1
MAX_X=1+lx/LSIZE
!*************定义table数组******************************************
*do,i,1,MAX_X,1
XX=(I-1)*LSIZE
CSYS,0
nsel,s,loc,y,ly+wheight*(j-1),ly+wheight*j
nsel,r,loc,z,lz/2-weld/2,lz/2+weld/2
esln,s,1
DDELE,ALL
allsel,all
!!!!!!!生死单元的处理
csys,0
nsel,s,loc,y,ly+wheight*(j-1),ly+wheight*j
nsel,r,loc,z,lz/2-weld/2,lz/2+weld/2
nsel,r,loc,x,XX-chf/2,XX+chf/2
esln,s,1
ealive,all
eplot
!!!!!!双椭球体热源的加载
ALLSEL,ALL
!!!双椭球热源的加载
!!!前半球节点的选取
local,11,2,XX,ly+wheight*j,lz/2,,,,bh/chf,ah/chf,
NSEL,S,LOC,X,0,chf
NSEL,R,LOC,Y,270,360
NSEL,R,LOC,Z,0,180
CM,FRONT,NODE
ALLSEL,ALL
!!!!!!!!!!!!!!!!!!!!!!!!!!!
LOCAL,13,0,XX,ly+wheight*j,lz/2, , , ,1,1,
*DEL,_FNCNAME
*DEL,_FNCMTID
*DEL,_FNC_C1
*DEL,_FNC_C2
*DEL,_FNC_C3
*DEL,_FNC_C4
*DEL,_FNCCSYS
*DEL,'front'
*SET,_FNCNAME,'front'
*DIM,_FNC_C1,,1
*DIM,_FNC_C2,,1
*DIM,_FNC_C3,,1
*DIM,_FNC_C4,,1
*SET,_FNC_C1(1),qmax11
*SET,_FNC_C2(1),chf
*SET,_FNC_C3(1),bh
*SET,_FNC_C4(1),ah
*SET,_FNCCSYS,13
! /INPUT,front.func,,,1
*DIM,%_FNCNAME%,TABLE,6,20,1,,,,%_FNCCSYS%
!
! Begin of equation: Qmax11*exp(-3(({x}/chf)^2+({y}/bh)^2+({z}/ah)^2))
*SET,%_FNCNAME%(0,0,1), 0.0, -999
*SET,%_FNCNAME%(2,0,1), 0.0
*SET,%_FNCNAME%(3,0,1), %_FNC_C1(1)%
*SET,%_FNCNAME%(4,0,1), %_FNC_C2(1)%
*SET,%_FNCNAME%(5,0,1), %_FNC_C3(1)%
*SET,%_FNCNAME%(6,0,1), %_FNC_C4(1)%
*SET,%_FNCNAME%(0,1,1), 1.0, -1, 0, 0, 0, 0, 0
*SET,%_FNCNAME%(0,2,1), 0.0, -2, 0, 1, 0, 0, -1
*SET,%_FNCNAME%(0,3,1), 0, -3, 0, 1, -1, 2, -2
*SET,%_FNCNAME%(0,4,1), 0.0, -1, 0, 1, 2, 4, 18
*SET,%_FNCNAME%(0,5,1), 0.0, -2, 0, 2, 0, 0, -1
*SET,%_FNCNAME%(0,6,1), 0.0, -4, 0, 1, -1, 17, -2
*SET,%_FNCNAME%(0,7,1), 0.0, -1, 0, 1, 3, 4, 19
*SET,%_FNCNAME%(0,8,1), 0.0, -2, 0, 2, 0, 0, -1
*SET,%_FNCNAME%(0,9,1), 0.0, -5, 0, 1, -1, 17, -2
*SET,%_FNCNAME%(0,10,1), 0.0, -1, 0, 1, -4, 1, -5
*SET,%_FNCNAME%(0,11,1), 0.0, -2, 0, 1, 4, 4, 20
*SET,%_FNCNAME%(0,12,1), 0.0, -4, 0, 2, 0, 0, -2
*SET,%_FNCNAME%(0,13,1), 0.0, -5, 0, 1, -2, 17, -4
*SET,%_FNCNAME%(0,14,1), 0.0, -2, 0, 1, -1, 1, -5
*SET,%_FNCNAME%(0,15,1), 0.0, -1, 0, 3, 0, 0, -2
*SET,%_FNCNAME%(0,16,1), 0.0, -4, 0, 1, -1, 3, -2
*SET,%_FNCNAME%(0,17,1), 0.0, -1, 0, 1, -3, 3, -4
*SET,%_FNCNAME%(0,18,1), 0.0, -1, 7, 1, -1, 0, 0
*SET,%_FNCNAME%(0,19,1), 0.0, -2, 0, 1, 17, 3, -1
*SET,%_FNCNAME%(0,20,1), 0.0, 99, 0, 1, -2, 0, 0
! End of equation: Qmax11*exp(-3(({x}/chf)^2+({y}/bh)^2+({z}/ah)^2))
!-->
BF,FRONT,HGEN,%front%
CMDELE,FRONT
ALLSEL,ALL
!!!后半球节点的选取
local,12,2,XX,ly+wheight*j,lz/2,,,,bh/chb,ah/chb,
NSEL,S,LOC,X,0,chb
NSEL,R,LOC,Y,180,270
NSEL,R,LOC,Z,0,180
CM,BEHIND,NODE
ALLSEL,ALL
!!!!!!!!!!!!!!!!!!!!!!!!!!!
LOCAL,14,0,XX,ly+wheight*j,lz/2, , , ,1,1,
*DEL,_FNCNAME
*DEL,_FNCMTID
*DEL,_FNC_C1
*DEL,_FNC_C2
*DEL,_FNC_C3
*DEL,_FNC_C4
*DEL,_FNCCSYS
*DEL,'back'
*SET,_FNCNAME,'back'
*DIM,_FNC_C1,,1
*DIM,_FNC_C2,,1
*DIM,_FNC_C3,,1
*DIM,_FNC_C4,,1
*SET,_FNC_C1(1),qmax12
*SET,_FNC_C2(1),chb
*SET,_FNC_C3(1),bh
*SET,_FNC_C4(1),ah
*SET,_FNCCSYS,0
! /INPUT,back.func,,,1
*DIM,%_FNCNAME%,TABLE,6,20,1,,,,%_FNCCSYS%
!
! Begin of equation: Qmax12*exp(-3(({x}/chb)^2+({y}/bh)^2+({z}/ah)^2))
*SET,%_FNCNAME%(0,0,1), 0.0, -999
*SET,%_FNCNAME%(2,0,1), 0.0
*SET,%_FNCNAME%(3,0,1), %_FNC_C1(1)%
*SET,%_FNCNAME%(4,0,1), %_FNC_C2(1)%
*SET,%_FNCNAME%(5,0,1), %_FNC_C3(1)%
*SET,%_FNCNAME%(6,0,1), %_FNC_C4(1)%
*SET,%_FNCNAME%(0,1,1), 1.0, -1, 0, 0, 0, 0, 0
*SET,%_FNCNAME%(0,2,1), 0.0, -2, 0, 1, 0, 0, -1
*SET,%_FNCNAME%(0,3,1), 0, -3, 0, 1, -1, 2, -2
*SET,%_FNCNAME%(0,4,1), 0.0, -1, 0, 1, 2, 4, 18
*SET,%_FNCNAME%(0,5,1), 0.0, -2, 0, 2, 0, 0, -1
*SET,%_FNCNAME%(0,6,1), 0.0, -4, 0, 1, -1, 17, -2
*SET,%_FNCNAME%(0,7,1), 0.0, -1, 0, 1, 3, 4, 19
*SET,%_FNCNAME%(0,8,1), 0.0, -2, 0, 2, 0, 0, -1
*SET,%_FNCNAME%(0,9,1), 0.0, -5, 0, 1, -1, 17, -2
*SET,%_FNCNAME%(0,10,1), 0.0, -1, 0, 1, -4, 1, -5
*SET,%_FNCNAME%(0,11,1), 0.0, -2, 0, 1, 4, 4, 20
*SET,%_FNCNAME%(0,12,1), 0.0, -4, 0, 2, 0, 0, -2
*SET,%_FNCNAME%(0,13,1), 0.0, -5, 0, 1, -2, 17, -4
*SET,%_FNCNAME%(0,14,1), 0.0, -2, 0, 1, -1, 1, -5
*SET,%_FNCNAME%(0,15,1), 0.0, -1, 0, 3, 0, 0, -2
*SET,%_FNCNAME%(0,16,1), 0.0, -4, 0, 1, -1, 3, -2
*SET,%_FNCNAME%(0,17,1), 0.0, -1, 0, 1, -3, 3, -4
*SET,%_FNCNAME%(0,18,1), 0.0, -1, 7, 1, -1, 0, 0
*SET,%_FNCNAME%(0,19,1), 0.0, -2, 0, 1, 17, 3, -1
*SET,%_FNCNAME%(0,20,1), 0.0, 99, 0, 1, -2, 0, 0
! End of equation: Qmax12*exp(-3(({x}/chb)^2+({y}/bh)^2+({z}/ah)^2))
!-->
BF,BEHIND,HGEN,%back%
CMDELE,BEHIND
ALLSEL,ALL
time,i*tinc+50*(j-1)
nsubst,10
SOLVE
BFDELE,all,HGEN
esel,s,live
eplot
OUTRES,ALL,ALL,
*ENDDO
*enddo
*DO,ii,1,120
csys,0
nsel,s,ext
sf,all,conv,30,25
ALLSEL,ALL
TIME,5*ii+250
solve
*ENDDO |