/title,Weld Analysis by "Element Birth and Death"
/FILNAME,welding_indirect,1
/prep7
et,1,90
! 1号材料是钢
! 2号材料是焊料MG-51T
！假设他们的刚度随温度是变化的。
MPTEMP,1,0
MPTEMP,2,500
MPTEMP,3,1000
MPTEMP,4,1500
MPDATA,EX,1,,2E+005
MPDATA,EX,1,,1.5e5
MPDATA,EX,1,,7e4
MPDATA,EX,1,,1e4
MPDATA,PRXY,1,,0.3
MPDATA,PRXY,1,,0.3
MPDATA,PRXY,1,,0.3
MPDATA,PRXY,1,,0.3
TB,MISO,1,4,7,0
TBTEMP,0
TBPT,,0.001089,217.8525472
TBPT,,0.019768838,255.5291284
TBPT,,0.039221021,284.700893
TBPT,,0.076978671,310.3430985
TBPT,,0.113353542,321.3356538
TBPT,,0.163817474,328.2458777
TBPT,,0.215140347,329.9893647
TBTEMP,500
TBPT,,0.001133333,170
TBPT,,0.019768838,207.6765812
TBPT,,0.039221021,236.8483458
TBPT,,0.076978671,262.4905513
TBPT,,0.113353542,273.4831066
TBPT,,0.163817474,280.3933305
TBPT,,0.215140347,282.1368175
TBTEMP,1000
TBPT,,0.001128571,79
TBPT,,0.019768838,116.6765812
TBPT,,0.039221021,145.8483458
TBPT,,0.076978671,171.4905513
TBPT,,0.113353542,182.4831066
TBPT,,0.163817474,189.3933305
TBPT,,0.215140347,191.1368175
TBTEMP,1500
TBPT,,0.0011,11
TBPT,,0.019768838,48.67658123
TBPT,,0.039221021,77.84834581
TBPT,,0.076978671,103.4905513
TBPT,,0.113353542,114.4831066
TBPT,,0.163817474,121.3933305
TBPT,,0.215140347,123.1368175
MPDATA,EX,2,,2E+005
MPDATA,EX,2,,1.5e5
MPDATA,EX,2,,7e4
MPDATA,EX,2,,1e4
MPDATA,PRXY,2,,0.3
MPDATA,PRXY,2,,0.3
MPDATA,PRXY,2,,0.3
MPDATA,PRXY,2,,0.3
TB,MISO,2,4,7,0
TBTEMP,0
TBPT,,0.001089,217.8525472
TBPT,,0.019768838,255.5291284
TBPT,,0.039221021,284.700893
TBPT,,0.076978671,310.3430985
TBPT,,0.113353542,321.3356538
TBPT,,0.163817474,328.2458777
TBPT,,0.215140347,329.9893647
TBTEMP,500
TBPT,,0.001133333,170
TBPT,,0.019768838,207.6765812
TBPT,,0.039221021,236.8483458
TBPT,,0.076978671,262.4905513
TBPT,,0.113353542,273.4831066
TBPT,,0.163817474,280.3933305
TBPT,,0.215140347,282.1368175
TBTEMP,1000
TBPT,,0.001128571,79
TBPT,,0.019768838,116.6765812
TBPT,,0.039221021,145.8483458
TBPT,,0.076978671,171.4905513
TBPT,,0.113353542,182.4831066
TBPT,,0.163817474,189.3933305
TBPT,,0.215140347,191.1368175
TBTEMP,1500
TBPT,,0.0011,11
TBPT,,0.019768838,48.67658123
TBPT,,0.039221021,77.84834581
TBPT,,0.076978671,103.4905513
TBPT,,0.113353542,114.4831066
TBPT,,0.163817474,121.3933305
TBPT,,0.215140347,123.1368175
! 材料密度(假设为常值)
mp,dens,1,8e-3
mp,dens,2,8e-3
! 热膨胀系数(假设为常值)
mp,alpx,1,1.2e-5
mp,alpx,2,1.2e-5
! 热传导系数(假设为常值)
mp,kxx,1,0.03
mp,kxx,2,0.03
! 比热(假设为常值)
mp,c,1,0.6
mp,c,2,0.6
! 假设焊料焊上去后的初始温度是1500℃
mp,reft,1,25
mp,reft,2,1500
!
! 下面建立几何模型
csys,0
k,1,7,0,0
k,2,7,0,4
k,3,56,0,4
k,4,56,0,11
k,5,60,0,11
k,6,60,0,8
k,7,60,0,0
k,8,63,0,8
k,9,63,0,11
k,12,62,0,8
k,13,62,0,-12
k,14,60,0,-12
k,10,0,0,0
k,11,0,0,1
a,1,2,3,4,5,6,7
a,5,6,8,9
a,7,6,12,13,14
ET,2,SHELL63
! 划分单元
esize,2
type,2
mat,1
amesh,1,3,2
mat,2
amesh,2
TYPE, 1
EXTOPT,ESIZE,25,0,
EXTOPT,ACLEAR,1
EXTOPT,ATTR,1,1,1
vrotat,all,,,,,,10,11,,2
cswpla,11,1
aclear,all
etdel,2
/solu
antype,4 ! 瞬态分析
trnopt,full
!
! 假设模型表面施加对流，冷却时也是！
allsel
sfa,all,1,conv,1.1e-5,25
! 由于第2，5体是焊料所在区域，因此，首先将该区域的单元"死"掉
vsel,s,,,2,5,3
eslv
*get,wnum,elem,,count
*dim,eorder,,wnum
*dim,ne,,wnum
*dim,ney,,wnum
*get,emax,elem,,num,max
*get,emin,elem,,num,min
mine=0
! 下面的DO循环用于将焊料区的单元按柱坐标csys11y方向排序，以便后面
! 模拟焊料逐步"生长"的过程
ii=0
*do,i,emin,emax
*if,esel(i),eq,1,then
ii=ii+1
ne(ii)=i
*endif
*enddo
*do,i,1,wnum
*get,ney(i),elem,ne(i),cent,y
*enddo
*do,i1,1,wnum
miny=1000
*do,i,1,wnum
*if,ney(i),lt,miny,then
miny=ney(i)
mine=ne(i)
j=i
*else
*if,ney(i),eq,miny,then
miny=ney(i)
mine=ne(i)
j=i
*endif
*endif
*enddo
eorder(i1)=mine
ney(j)=1001
*enddo
max_tem=1500 ! 按照前面假设，焊料的初始温度为1500℃
dt=1! 焊接8个单元所需的时间
t=0 ! 起始时间
allsel
save,welding_indirect,db
nsubst,1
*do,i,1,wnum,8
ddele,all,temp
*do,j,1,20
！此步单元的温度加载,在每个出生单元的20个节点上施加温度
d,nelem(eorder(i),j),temp,max_tem
d,nelem(eorder(i+1),j),temp,max_tem
d,nelem(eorder(i+2),j),temp,max_tem
d,nelem(eorder(i+3),j),temp,max_tem
d,nelem(eorder(i+4),j),temp,max_tem
d,nelem(eorder(i+5),j),temp,max_tem
d,nelem(eorder(i+6),j),temp,max_tem
d,nelem(eorder(i+7),j),temp,max_tem
*enddo
t=t+dt
time,t
nsubst,1
solve
*enddo
!冷却过程分析
time,5000
ddele,all,temp
nsubst,1
solve
finish
!热分析结束，开始结构分析。
/prep7
etchg,tts
tunif,25
! 在模型加位移约束
asel,s,,,18,33,15
da,all,ux,0
da,all,uy,0
da,all,uz,0
allsel
finish
/solu
antype,0
！先将焊接单元杀死（相当于杀死单元提供原来刚度等参数都乘以一个小量，默认是1e-6）
vsel,s,,,2,5,3
eslv
ekill,all
allsel
LDREAD,TEMP,25,1,, ,'welding_indirect','rth',
nsubst,5,100,1
time,1
solve
！激活焊接单元。
ealive,all
LDREAD,TEMP,26,1,, ,'welding_indirect','rth',
time,2
nsubst,100,10000,1
solve