焊接ansys-生死单元实例及解释

/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.0010,217.8525472

TBPT,,0.019768838,255.5291284

TBPT,,0.039221021,284.7003

TBPT,,0.076978671,310.3430985

TBPT,,0.113353542,321.3356538

TBPT,,0.163817474,328.2458777

TBPT,,0.215140347,329.937

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,1.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.0010,217.8525472

TBPT,,0.019768838,255.5291284

TBPT,,0.039221021,284.7003

TBPT,,0.076978671,310.3430985

TBPT,,0.113353542,321.3356538

TBPT,,0.163817474,328.2458777

TBPT,,0.215140347,329.937

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,1.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