!杀灭闸墩单元 nsel,s,loc,y,5,10 esln,s,1 ekill,all
!激活第一层单元 nsel,s,loc,y,5,8 esln,s,1 ealive,all
nsel,s,loc,y,0,8 nsel,inve d,all,temp,0 allsel
!!!!!!!!!!!!FINISH
!!!!!!/PREP7 !!(如果不加此句，循环后所出的图形时间是从第六天开始的，相当于Main Menu>Preprocessor ） delt=.2 tim=0.001
timint,on !打开时间积分 TINTP,,,,1.0,,
*do,day,1,10 !!!!!/prep7 tim=tim+delt
time,tim !指定荷载步的结束时间
deltim,0.2 !设置荷载步的子步长(即一天内再分成1/0.5段） AUTOS,-1 !打开自动时间步长
DIFHgen=20498472.32*exp(-0.36*((tim-0.1)**0.74))*((tim-0.1)**(-0.26)) !VSEL,ALL
BFV,2,HGEN,DIFHgen !加生热荷载 KBC,0 !荷载线性变化（ramped) SFA,8,1,CONV,14,10 SFA,9,1,CONV,14,10 SFA,10,1,CONV,14,10 SFA,11,1,CONV,14,10 ALLSEL /solu
SOLVE *enddo
!!!!!!开始应力场的求解 /prep7 etchg,tts
mp,alpx,1,0.7e-5 !!!热膨胀系数
mp,ex,1,29.5e9 !!!90天的弹性模量
mp,prxy,1,0.167 !!!泊松比 mp,dens,1,2441 !!!!密度
!!!!!!!!利用循环实现弹性模量随时间变化，按顺序激活单元，并导入每个步长时间段的温度场，得出变化的应力场 time,0.001 nropt,full
!!!!!!!!!!!!solve tim=0.001 delt=0.2
!!!!!!!!先把2,3层单元全部杀死，再激活第一层 !杀灭闸墩单元 nsel,s,loc,y,5,10 esln,s,1 ekill,all
!激活第一层单元 nsel,s,loc,y,5,8 esln,s,1 ealive,all
nsel,s,loc,y,0,8 nsel,inve d,all,all,0 allsel
*do,tm,1,10 tim=tim+delt time,tim deltim,0.2
!改变第1层弹性模量 /prep7
ecx=1.41*27.68*10**9*(1-exp(-0.4*tim**0.34)) mp,ex,1,ecx nsel,s,loc,y,5,8 esln,s,1 emodi,all finish /solu
!!!!!!结构计算的边界条件 nsel,s,loc,y,5,8 nsel,r,loc,x,0 d,all,ux,0 alls
nsel,s,loc,y,5,8 nsel,r,loc,x,12 d,all,ux,0 alls
nsel,s,loc,y,5,8 nsel,r,loc,z,0 d,all,uz,0 alls
nsel,s,loc,y,5,8 nsel,r,loc,z,7 d,all,uz,0 alls
!!!!!!!!!!考虑重力因素 grav=9.81
acel 0,grav,0 !!!!!!!!!!施加y向重力加速度
!!!!!!!读入热分析结果 !allsel,all !nsel,all
!BFVDELE,all,TEMP
!SFADELE,all,TEMP !这四句好像没有什么用 !!!!!ldread,temp,,,tim,,,rth ldread,temp,0,0,tim,,,rth tref,0 solve
*enddo !!!结束56天的循环 fini
【利用ANSYS单元生死和热-结构耦合分析】
!下面的命令流进行的是一个简单的二维焊接分析, 利用ANSYS单元生死和热-结构耦合分析功能进
!行焊接过程仿真, 计算焊接过程中的温度分布和应力分布以及冷却后的焊缝残余应力。
/title,Weld Analysis by \ /prep7
/unit,si ! 采用国际单位制 !
et,1,13,4 ! 13号二维耦合单元, 同时具有温度和位移自由度 et,2,13,4 !
! 1号材料是钢
! 2号材料是铝 ! 3号材料是铜
! 铝是本次分析中的焊料, 它将钢结构部分和铜结构部分焊接起来 ! 下面是在几个温度点下, 各材料的弹性模量 mptemp,1,20,500,1000,1500,2000
mpdata,ex,1,1,1.93e11,1.50e11,0.70e11,0.10e11,0.01e11 mpdata,ex,2,1,1.02e11,0.50e11,0.08e11,0.001e11,0.0001e11 mpdata,ex,3,1,1.17e11,0.90e11,0.30e11,0.05e11,0.005e11 !
! 假设各材料都是双线性随动硬化弹塑性本构关系
! 下面是各材料在各温度点下的屈服应力和屈服后的弹性模量 tb,bkin,1,5 tbtemp,20,1
tbdata,1,1200e6,0.193e11 tbtemp,500,2
tbdata,1, 933e6,0.150e11 tbtemp,1000,3
tbdata,1, 435e6,0.070e11 tbtemp,1500,4
tbdata,1, 70e6,0.010e11 tbtemp,2000,5
tbdata,1, 7e6,0.001e11 !
tb,bkin,2,5 tbtemp,20,1
tbdata,1,800e6,0.102e11 tbtemp,500,2
tbdata,1,400e6,0.050e11 tbtemp,1000,3
tbdata,1, 70e6,0.008e11 tbtemp,1500,4
tbdata,1, 1e6,0.0001e11 tbtemp,2000,5
tbdata,1,0.1e6,0.00001e11 !
tb,bkin,3,5 tbtemp,20,1
tbdata,1,900e6,0.117e11 tbtemp,500,2
tbdata,1,700e6,0.090e11 tbtemp,1000,3
tbdata,1,230e6,0.030e11 tbtemp,1500,4
tbdata,1, 40e6,0.005e11
tbtemp,2000,5
tbdata,1, 4e6,0.0005e11 !
! 材料密度(假设为常值) mp,dens,1,8030 mp,dens,2,4850 mp,dens,3,8900
! 热膨胀系数(假设为常值) mp,alpx,1,1.78e-5 mp,alpx,2,9.36e-6 mp,alpx,3,1.66e-5
! 泊松比(假设为常值) mp,nuxy,1,0.29 mp,nuxy,2,0.30 mp,nuxy,3,0.30
! 热传导系数(假设为常值) mp,kxx,1,16.3 mp,kxx,2,7.44 mp,kxx,3,393
! 比热(假设为常值) mp,c,1,502 mp,c,2,544 mp,c,3,385
! 热膨胀系数(假设为常值)
! 由于该13号单元还有磁自由度, 此处假设一磁特性, 但并不影响我们所关心的结果 mp,murx,1,1 mp,murx,2,1 mp,murx,3,1
! 假设焊料(铝)焊上去后的初始温度是1500℃ mp,reft,1,20 mp,reft,2,1500 mp,reft,3,20 !
! 下面建立几何模型 csys,0 k,1,0,0,0 k,2,0.5,0,0 k,3,1,0,0 k,4,0,0.3,0 k,5,0.35,0.3,0 k,6,0.65,0.3,0 k,7,1,0.3,0 a,1,2,5,4 a,2,6,5