当前位置:首页 >> 机械/仪表 >>

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


/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

相关文章:
Ansys单元生死功能
Ansys单元生死功能_IT/计算机_专业资料。Ansys 单元...本过 程对于由相变引起的模型效应 (如焊接过程中原...下面的例子是杀死总应变超过许用值的单元: /SOLU ...
单元生死法
本过程对于由相变引起的模型 效应(如焊接过程中原不...使用 ANSYS 结果控制单元生死 在许多时候,用户并不...下面的例子是杀死总应变超过许用值的单元: /SOLU ...
ANSYS Example05单元生死算例 (ANSYS)
ANSYS Example05单元生死算例 (ANSYS)_工学_高等教育_教育专区。ANSYS Example05单元生死算例 (ANSYS)05 单元生死算例 (ANSYS) 土木工程中经常需要对施工过程进行...
ANSYS单元生死总结2
ANSYS单元生死总结2_机械/仪表_工程科技_专业资料。...本过程对于由相变引起的模型效 应 (如焊接过程中原...这 样会导致不收敛 求解命令流实例: NLGEOM,ON !...
ANSYS生死单元应用总结
比如超高层建筑与大跨度结构的施工过程模拟分析,钢结 构焊接与退火的模拟,隧道...ANSYS生死单元 2387人阅读 22页 2下载券 单元生死算例 (ANSYS) 文... 1385...
ANSYS单元的生和死教程
单元生死选项ANSYS 高级分析技术指南 SIMU-ONLINE BBS...本过 程对于由相变引起的模型效应 (如焊接过程中原...下面的例子是杀死总应变超过许用值的单元: /SOLU ...
ansys利用单元生死模拟拱桥浇筑过程
ANSYS单元死活应用探讨 5页 免费 生死单元技术 5页 免费a​n​s​y​...a​n​s​y​s​应​用​实​例​,​很​有​借​鉴...
ANSYS单元生死功能模拟门式刚架施工
Ansys单元生死功能 7页 免费 基于ANSYS单元生死技术的... 4页 1财富值 ANSYS...Ansys单元生死功能 7页 免费 单元生死算例 (ANSYS) 3页 2财富值 3第四讲 ...
基于ANSYS的焊接残余应力的有限元分析
降低了焊接件的刚度和稳定性,采用有限元分析软件 ANSYS 的热 —结构耦合和单元生死分析功能对钢板焊接过程进行模拟仿真,分析焊接残余应力分布的规 律,为复杂焊接结构...
钎焊解释
一些焊接名词的解释 4页 2下载券 焊接ansys-生死单元实例... 2页 1下载券 ...钎焊简介 钎焊简介 钎焊是采用比母材熔点低的金属材料作钎料,将焊件(母材)与...
更多相关标签: