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

焊接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生死单元应用总结
ANSYS 生死单元应用总结 摘要:随着超高层建筑、大跨度结构、地下工程等大量的涌现...比如 超高层建筑与大跨度结构的施工过程模拟分析,钢结构焊接与退火 的模拟,隧道...
ANSYS单元生死总结2
ANSYS单元生死总结2_机械/仪表_工程科技_专业资料。...本过程对于由相变引起的模型效 应 (如焊接过程中原...这 样会导致不收敛 求解命令流实例: NLGEOM,ON !...
ANSYS Example05单元生死算例 (ANSYS)
ANSYS Example05单元生死算例 (ANSYS)_工学_高等教育_教育专区。ANSYS Example05单元生死算例 (ANSYS)05 单元生死算例 (ANSYS) 土木工程中经常需要对施工过程进行...
ANSYS生死单元应用总结
ANSYS生死单元应用总结_机械/仪表_工程科技_专业资料。龙源期刊网 http://www....比如超高层建筑与大跨度结构的施工过程模拟分析,钢结 构焊接与退火的模拟,隧道...
单元生死算例 (ANSYS) 文档
单元生死算例 (ANSYS) 文档_建筑/土木_工程科技_专业资料。土木工程中经常需要对施工过程进行模拟。 很多复杂工程构件的最不利 受力状态往往未必是在结构完工以后,...
ANSYS单元的生和死教程
单元生死选项ANSYS 高级分析技术指南 SIMU-ONLINE BBS...本过 程对于由相变引起的模型效应 (如焊接过程中原...下面的例子是杀死总应变超过许用值的单元: /SOLU ...
ansys分析实例集
简单的二维焊接分析-单元生死实例下面的命令流进行的是一个简单的二维焊接分析, 接 过程仿真, 利用 ANSYS 单元生死和热-结构耦合分析功能进行焊 计算焊接过程中的...
ANSYS单元生死功能模拟门式刚架施工
Ansys单元生死功能 7页 免费 基于ANSYS单元生死技术的... 4页 1财富值 ANSYS...Ansys单元生死功能 7页 免费 单元生死算例 (ANSYS) 3页 2财富值 3第四讲 ...
更多相关标签:
ansys生死单元实例 | ansys焊接模拟实例 | ansys焊接实例 | ansys生死单元 | ansys生死单元教程 | ansys单元生死技术 | 生死单元 模态 ansys | ansys workbench实例 |