焊接过程生死单元程序

finish
/nopr
/CONFIG,nres,3000 !设定结果集中载荷步的步数
/CONFIG,norstgm,0 !结果集中包含几何结构的信息
/CONFIG,nproc,1 !设定处理器个数
/FILNAME, welding_on_plate, 1
/prep7

width=0.05
wwidth=0.0046
wdepth=0.008
wheight=0.001
length=0.05
K,1,wwidth,0,0,
K,2,width,0,0,
K,3,width,0,-0.01,
K,4,-width,0,-0.01,
K,5,-width,0,0,
K,6,-wwidth,0,0,
K,7,0,0,-wdepth,
K,8,0,0,0.001,
LSTR,1,2
LSTR,2,3
LSTR,3,4
LSTR,4,5
LSTR,6,5
LSTR,6,7
LSTR,7,1
LARC,6,1,8
AL,5,4,3,2,1,7,6
AL,6,7,8
VOFFST,1,-length,
VOFFST,2,-length,
wprot,0,0,90
VSBW,1
VSBW,2
FLST,2,2,6,ORDE,2
FITEM,2,1
FITEM,2,3
VDELE,P51X, , ,1
NUMCMP,ALL
wpoff,0,0,0.015
VSBW,1
wpoff,0,0,0.015
VSBW,3
VSEL,ALL
VGLUE,ALL
NUMCMP,ALL
ET,1,SOLID70
MPTEMP,,,,,,,,
MPTEMP,1,298
MPTEMP,2,398
MPTEMP,3,498
MPTEMP,4,598
MPTEMP,5,698
MPTEMP,6,798
MPTEMP,7,898
MPTEMP,8,998
MPTEMP,9,1098
MPDATA,KXX,1,,31
MPDATA,KXX,1,,32.5
MPDATA,KXX,1,,33.5
MPDATA,KXX,1,,34
MPDATA,KXX,1,,33
MPDATA,KXX,1,,33.5
MPDATA,KXX,1,,33.8
MPDATA,KXX,1,,34
MPDATA,KXX,1,,30
MPTEMP,,,,,,,,
MPTEMP,1,298
MPTEMP,2,398
MPTEMP,3,498
MPTEMP,4,598
MPTEMP,5,698
MPTEMP,6,798
MPTEMP,7,898
MPTEMP,8,998
MPTEMP,9,1098
MPDATA,C,1,,410
MPDATA,C,1,,500
MPDATA,C,1,,530
MPDATA,C,1,,570
MPDATA,C,1,,600
MPDATA,C,1,,720
MPDATA,C,1,,810
MPDATA,C,1,,1080
MPDATA,C,1,,800
MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,DENS,1,,7890
MPTEMP,,,,,,,,
MPTEMP,1,298
MPTEMP,2,398
MPTEMP,3,498
MPTEMP,4,598
MPTEMP,5,698
MPTEMP,6,798
MPTEMP,7,898
MPTEMP,8,998
MPTEMP,9,1098
MPDATA,HF,1,,100
MPDATA,HF,1,,350
MPDATA,HF,1,,480
MPDATA,HF,1,,520
MPDATA,HF,1,,1000
MPDATA,HF,1,,1050
MPDATA,HF,1,,1500
MPDATA,HF,1,,2000
MPDATA,HF,1,,3000
MPTEMP,,,,,,,,
MPTEMP,1,298
MPTEMP,2,398
MPTEMP,3,498
MPTEMP,4,598
MPTEMP,5,698
MPTEMP,6,798
MPTEMP,7,898
MPTEMP,8,998
MPTEMP,9,1098
MPDATA,KXX,2,,31
MPDATA,KXX,2,,32.5
MPDATA,KXX,2,,33.5
MPDATA,KXX,2,,34
MPDATA,KXX,2,,33
MPDATA,KXX,2,,33.5
MPDATA,KXX,2,,33.8
MPDATA,KXX,2,,34
MPDATA,KXX,2,,30
MPTEMP,,,,,,,,
MPTEMP,1,298
MPTEMP,2,398
MPTEMP,3,498
MPTEMP,4,598
MPTEMP,5,698
MPTEMP,6,798
MPTEMP,7,898
MPTEMP,8,998
MPTEMP,9,1098
MPDATA,C,2,,410
MPDATA,C,2,,500
MPDATA,C,2,,530
MPDATA,C,2,,570
MPDATA,C,2,,600
MPDATA,C,2,,720
MPDATA,C,2,,810
MPDATA,C,2,,1080
MPDATA,C,2,,800
MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,DENS,2,,7890
MPTEMP,,,,,,,,
MPTEMP,1,298
MPTEMP,2,398
MPTEMP,3,498
MPTEMP,4,598
MPTEMP,5,698
MPTEMP,6,798
MPTEMP,7,898
MPTEMP,8,998
MPTEMP,9,1098
MPDATA,HF,2,,100
MPDATA,HF,2,,350
MPDATA,HF,2,,480
MPDATA,HF,2,,520
MPDATA,HF,2,,1000
MPDATA,HF,2,,1050
MPDATA,HF,2,,1500
MPDATA,HF,2,,2000
MPDATA,HF,2,,3000
SAVE
MAT,

1
LESIZE,2,,,10,
LESIZE,14,,,10,
LESIZE,1,,,6,
LESIZE,6,,,6,
EXTOPT,ACLEAR,1
EXTOPT,VSWE,TETS,1
EXTOPT,VSWE,AUTO,0
EXTOPT,ESIZE,24,0
VSWEEP,1,2,4
!!!!!!!!
MAT,1
LESIZE,4,,,10,
LESIZE,3,,,10,
LESIZE,18,,,10,
EXTOPT,ACLEAR,1
EXTOPT,VSWE,TETS,1
EXTOPT,VSWE,AUTO,0
EXTOPT,ESIZE,50,0
VSWEEP,2,8,9
!!!!!!!!!!
MAT,1
LESIZE,26,,,10,
LESIZE,22,,,10,
LESIZE,27,,,10,
LESIZE,13,,,4,
EXTOPT,ESIZE,50,0
VSWEEP,3,16,13
!!!!!!!!!!!
MAT,1
LESIZE,33,,,6,
LESIZE,34,,,6,
LESIZE,35,,,6,
LESIZE,36,,,6,
mshkey,1
mshape,0,3d
vmesh,4
!!!!!!!!!!
wpcsys,-1,0
ALLSEL,ALL
!!!!!!!!!!!!杀死单元
VSEL,S, , , 2
ESLV,S
CM,weld_elem,ELEM
ekill,all
esel,s,elem,,weld_elem
nsle,s
cm,weld_node,node
ALLSEL,ALL
!!!!!!!!!!!!!焊接物理参数的设置
D=0.005
ah=D/2
chf=2*ah
chb=8*ah
bh=0.010
rf=2*chf/(chf+chb)
rb=2*chb/(chf+chb)
M=0.73
W0=3500
W=W0*M
STEP=0.002
N=length/STEP !!!!!!25
V=0.002
SAVE
!!!!!!!!!!求解阶段



/SOLU
ANTYPE,4
TRNOPT,FULL
TUNIF,298,
OUTRES,NSOL,all

NROPT,FULL,,ON
AUTOTS,1
NSUBST,5, , ,1 !!!
KBC,1
lnsrch,1 !!!
TSRES,ERASE
TINTP,,,,1,,,,,!!!
TIMINT,0,struct
TIMINT,1,THERM
TIMINT,0,MAG
!!!!!!!!!!!!!边界条件的施加
NSEL,S,EXT !选择外部节点
NSEL,U,LOC,X,0 !不选x=0的节点
SF,ALL,CONV,200,293 !在节点上施加对流换热系数
ALLSEL,ALL !选择全部图元
!!!!!!!!!!!!!!!!!!!!!!
*DO,I,0,N
YY=I*STEP
T=(YY+STEP)/V !!!!!!!每步焊接时间，时间跟踪
TIME,T
!!!!!!!生死单元的处理
wpcsys,-1,0
nsel,s,node,,weld_node
nsel,r,loc,y,YY,YY+chf
cm,cur_node,node
esln,s
ealive,all
nsel,s,node,,weld_node
nsel,u,node,,cur_node
cm,weld_node,node

!!!!!!双椭球体热源的加载
ALLSEL,ALL
!!!双椭球热源的加载
!!!前半球节点的选取

LOCAL,11,2,0,YY,0, , , ,chf/ah,bh/ah,
NSEL,S,LOC,X,0,ah
NSEL,R,LOC,Y,,90
NSEL,R,LOC,Z,-90,90
CM,FRONT,NODE
ALLSEL,ALL

!!!!!!!!!!!!!!!!!!!!!!!!!!!
LOCAL,13,0,0,YY,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,'TFRONT' !!!!!!!!
*SET,_FNCNAME,'TFRONT'
*DIM,_FNC_C1,,1
*DIM,_FNC_C2,,1
*DIM,_FNC_C3,,1
*DIM,_FNC_C4,,1
*DIM,_FNC_C5,,1
*SET,_FNC_C1(1),rf !!!rf
*SET,_FNC_C2(1),w !!!Q
*SET,_FNC_C3(1),chf !!!chf
*SET,_FNC_C4(1),ah !!!ah
*SET,_FNC_C5(1),bh !!!bh
*SET,_FNCCSYS,13
*DIM,%_FNCNAME%,TABLE,7,38,1,,,,%_FNCCSYS%
*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%(7,0,1), %_F

NC_C5(1)%
*SET,%_FNCNAME%(0,1,1), 1.0, -1, 16, 3, 0, 0, 0
*SET,%_FNCNAME%(0,2,1), 0.0, -2, 0, 6, 0, 0, -1
*SET,%_FNCNAME%(0,3,1), 0, -3, 0, 1, -2, 3, -1
*SET,%_FNCNAME%(0,4,1), 0.0, -1, 0, 1, -3, 3, 17
*SET,%_FNCNAME%(0,5,1), 0.0, -2, 0, 1, -1, 3, 18
*SET,%_FNCNAME%(0,6,1), 0.0, -1, 0, 0, 0, 0, 0
*SET,%_FNCNAME%(0,7,1), 0.0, -3, 0, 1, 0, 0, -1
*SET,%_FNCNAME%(0,8,1), 0.0, -4, 0, 1, -1, 2, -3
*SET,%_FNCNAME%(0,9,1), 0.0, -1, 0, 2, 0, 0, 2
*SET,%_FNCNAME%(0,10,1), 0.0, -3, 0, 1, 2, 17, -1
*SET,%_FNCNAME%(0,11,1), 0.0, -1, 0, 2, 0, 0, 19
*SET,%_FNCNAME%(0,12,1), 0.0, -5, 0, 1, 19, 17, -1
*SET,%_FNCNAME%(0,13,1), 0.0, -1, 0, 1, -3, 4, -5
*SET,%_FNCNAME%(0,14,1), 0.0, -3, 0, 2, 0, 0, 3
*SET,%_FNCNAME%(0,15,1), 0.0, -5, 0, 1, 3, 17, -3
*SET,%_FNCNAME%(0,16,1), 0.0, -3, 0, 2, 0, 0, 20
*SET,%_FNCNAME%(0,17,1), 0.0, -6, 0, 1, 20, 17, -3
*SET,%_FNCNAME%(0,18,1), 0.0, -3, 0, 1, -5, 4, -6
*SET,%_FNCNAME%(0,19,1), 0.0, -5, 0, 1, -1, 1, -3
*SET,%_FNCNAME%(0,20,1), 0.0, -1, 0, 2, 0, 0, 4
*SET,%_FNCNAME%(0,21,1), 0.0, -3, 0, 1, 4, 17, -1
*SET,%_FNCNAME%(0,22,1), 0.0, -1, 0, 2, 0, 0, 21
*SET,%_FNCNAME%(0,23,1), 0.0, -6, 0, 1, 21, 17, -1
*SET,%_FNCNAME%(0,24,1), 0.0, -1, 0, 1, -3, 4, -6
*SET,%_FNCNAME%(0,25,1), 0.0, -3, 0, 1, -5, 1, -1
*SET,%_FNCNAME%(0,26,1), 0.0, -1, 0, 3, 0, 0, -3
*SET,%_FNCNAME%(0,27,1), 0.0, -5, 0, 1, -1, 3, -3
*SET,%_FNCNAME%(0,28,1), 0.0, -1, 0, 1, -4, 3, -5
*SET,%_FNCNAME%(0,29,1), 0.0, -1, 7, 1, -1, 0, 0
*SET,%_FNCNAME%(0,30,1), 0.0, -3, 0, 1, -2, 3, -1
*SET,%_FNCNAME%(0,31,1), 0.0, -1, 0, 1, 20, 3, 21
*SET,%_FNCNAME%(0,32,1), 0.0, -2, 0, 1, -1, 3, 19
*SET,%_FNCNAME%(0,33,1), 0.0, -1, 0, 3.14159265358979310, 0, 0, -2
*SET,%_FNCNAME%(0,34,1), 0.0, -4, 0, 1, -2, 3, -1
*SET,%_FNCNAME%(0,35,1), 0.0, -1, 16, 3.14159265358979310, 0, 0, 0
*SET,%_FNCNAME%(0,36,1), 0.0, -2, 0, 1, -4, 3, -1
*SET,%_FNCNAME%(0,37,1), 0.0, -1, 0, 1, -3, 4, -2
*SET,%_FNCNAME%(0,38,1), 0.0, 99, 0, 1, -1, 0, 0
!!!!!!!!!!!!!!!!!!
BF,FRONT,HGEN, %TFRONT%
CMDELE,FRONT
ALLSEL,ALL
!!!后半球节点的选取
LOCAL,12,2,0,YY,wheight, , , ,chb/ah,bh/ah,
NSEL,S,LOC,X,0,ah
NSEL,R,LOC,Y,270,360
NSEL,R,LOC,Z,-90,90
CM,BEHIND,NODE

!!!后半球加载
!!!!!!!!!!!!!!!!!!!
LOCAL,14,0,0,YY,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),rb !!!rb
*SET,_FNC_C2(1),w !!!Q
*SET,_FNC_C3(1),chb !!!chb
*SET,_FNC_C4(1),ah !!!ah
*SET,_FNC_C5(1),bh !!!bh
*SET,_FNCCSYS,14 !!! 0
*DIM,%_FNCNAME%,TABLE,7,38,1,,,,%_FNCCSYS%
*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%(7,0,1), %_FNC_C5(1)%
*SET,%_FNCNAME%(0,1,1), 1.0, -1, 16, 3, 0, 0, 0
*SET,%_FNCNAME%(0,2,1), 0.0, -2, 0, 6, 0, 0, -1
*SET,%_FNCNAME%(0,3,1), 0, -3, 0, 1, -2, 3, -1
*SET,%_FNCNAME%(0,4,1), 0.0, -1, 0, 1, -3, 3, 17
*SET,%_FNCNAME%(0,5,1), 0.0, -2, 0, 1, -1, 3, 18
*SET,%_FNCNAME%(0,6,1), 0.0, -1, 0, 0, 0, 0, 0
*SET,%_FNCNAME%(0,7,1), 0.0, -3, 0, 1, 0, 0, -1
*SET,%_FNCNAME%(0,8,1), 0.0, -4, 0, 1, -1, 2, -3
*SET,%_FNCNAME%(0,9,1), 0.0, -1, 0, 2, 0, 0, 2
*SET,%_FNCNAME%(0,10,1), 0.0, -3, 0, 1, 2, 17, -1
*SET,%_FNCNAME%(0,11,1), 0.0, -1, 0, 2, 0, 0, 19
*SET,%_FNCNAME%(0,12,1), 0.0, -5, 0, 1, 19, 17, -1
*SET,%_FNCNAME%(0,13,1), 0.0, -1, 0, 1, -3, 4, -5
*SET,%_FNCNAME%(0,14,1), 0.0, -3, 0, 2, 0, 0, 3
*SET,%_FNCNAME%(0,15,1), 0.0, -5, 0, 1, 3, 17, -3
*SET,%_FNCNAME%(0,16,1), 0.0, -3, 0, 2, 0, 0, 20
*SET,%_FNCNAME%(0,17,1), 0.0, -6, 0, 1, 20, 17, -3
*SET,%_FNCNAME%(0,18,1), 0.0, -3, 0, 1, -5, 4, -6
*SET,%_FNCNAME%(0,19,1), 0.0, -5, 0, 1, -1, 1, -3
*SET,%_FNCNAME%(0,20,1), 0.0, -1, 0, 2, 0, 0, 4
*SET,%_FNCNAME%(0,21,1), 0.0, -3, 0, 1, 4, 17, -1
*SET,%_FNCNAME%(0,22,1), 0.0, -1, 0, 2, 0, 0, 21
*SET,%_FNCNAME%(0,23,1), 0.0, -6, 0, 1, 21, 17, -1
*SET,%_FNCNAME%(0,24,1), 0.0, -1, 0, 1, -3, 4, -6
*SET,%_FNCNAME%(0,25,1), 0.0, -3, 0, 1, -5, 1, -1
*SET,%_FNCNAME%(0,26,1), 0.0, -1, 0, 3, 0, 0, -3
*SET,%_FNCNAME%(0,27,1), 0.0, -5, 0, 1, -1, 3, -3
*SET,%_FNCNAME%(0,28,1), 0.0, -1, 0, 1, -4, 3, -5
*SET,%_FNCNAME%(0,29,1), 0.0, -1, 7, 1, -1, 0, 0
*SET,%_FNCNAME%(0,30,1), 0.0, -3, 0, 1, -2, 3, -1
*SET,%_FNCNAME%(0,31,1), 0.0, -1, 0, 1, 20, 3, 21
*SET,%_FNCNAME%(0,32,1), 0.0, -2, 0, 1, -1, 3, 19
*SET,%_FNCNAME%(0,33,1), 0.0, -1, 0, 3.14159265358979310, 0, 0, -2
*SET,%_FNCNAME%(0,34,1), 0.0, -4, 0, 1, -2, 3, -1
*SET,%_FNCNAME%(0,35,1), 0.0, -1, 16, 3.14159265358979310, 0, 0, 0
*SET,%_FNCNAME%(0,36,1), 0.0, -2, 0, 1, -4, 3, -1
*SET,%_FNCNAME%(0,37,1), 0.0, -1, 0, 1, -3, 4, -2
*SET,%_FNCNAME%(0,38,1), 0.0, 99, 0, 1, -1, 0, 0
BF,BEHIND,HGEN, %TBEHIND%
CMDELE,BEHIND
ALLSEL,ALL

SOLVE
BFDELE,all,HGEN
*ENDDO
SAVE
!!!添加冷却载荷第二步
TIME,length/V+5
NSUBST,15, , ,1
!!!DELTIM,0.05,0.01,0.1,1
SOLVE
!!!添加冷却载荷第二步
CSYS,0
TIME,length/V+20
NSUBST,15, , ,1
!!!DELTIM,0.2,0.1,0.4,1
SOLVE