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

焊接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实例
焊接模拟 ansys 实例 !下面的命令流进行的是一个简单的二维焊接分析, 利用 ANSYS 单元生死和热-结构 耦合分析功能进 !行焊接过程仿真, 计算焊接过程中的温度分布...
Ansys 单元生死一例
Ansys 单元生死一例_机械/仪表_工程科技_专业资料。Ansys 单元生死一例默认分类 2010-04-07 11:32:02 阅读 258 评论 0 finish /clear 字号:大中小 订阅 /...
ANSYS生死单元应用总结
ANSYS 生死单元应用总结 摘要:随着超高层建筑、大跨度结构、地下工程等大量的涌现...比如 超高层建筑与大跨度结构的施工过程模拟分析,钢结构焊接与退火 的模拟,隧道...
Ansys单元生死功能
Ansys单元生死功能_IT/计算机_专业资料。Ansys 单元...本过 程对于由相变引起的模型效应 (如焊接过程中原...下面的例子是杀死总应变超过许用值的单元: /SOLU ...
Ansys单元生死功能
OPT,ESTIF,ESEL,EKILL,NSLE,NSEL,D,F,SF 和 BF 命 令得到更详细的解释。...焊接ansys-生死单元实例... 2页 1下载券 ansys单元生死解析 22页 1下载券 ...
ANSYS Example05单元生死算例 (ANSYS)
例子将利用 ANSYS 提供的单元生死”功能 来进行一个门式框架的施工模拟 施工分为三步 1: 建立立柱和临时支撑 2: 安装横梁 3: 去掉临时支撑 知识要点 (1...
ANSYS生死单元应用总结
ANSYS生死单元应用总结_机械/仪表_工程科技_专业资料。龙源期刊网 http://www....比如超高层建筑与大跨度结构的施工过程模拟分析,钢结 构焊接与退火的模拟,隧道...
ANSYS单元生死总结2
ANSYS单元生死总结2_机械/仪表_工程科技_专业资料。...本过程对于由相变引起的模型效 应 (如焊接过程中原...这 样会导致不收敛 求解命令流实例: NLGEOM,ON !...
焊缝焊接收缩量的ANSYS仿真分析
本文运用 ANSYS 的热分析功能对焊接的收缩进行仿 真。该仿真存在的难点是热结构耦合、单元生死、材料的弹塑性、AP 今日推荐 89份文档 ...
更多相关标签:
ansys生死单元实例 | ansys焊接生死单元 | ansys焊接模拟实例 | ansys生死单元 | ansys中生死单元 | ansys生死单元命令流 | ansys静力学分析实例 | ansys谐响应分析实例 |