资源描述：
《生死单元技术.doc》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。
1、最近在做焊接模拟分析，做了一点东西和大家一起分享一下，还有很多问题是需要和大家一起讨论的！欢迎高手指正，谢谢〜题在两个物体Z间焊接一圈悍缝，分析冷却麻的结构变形和应力！模熨见apdl（为了计算速度我简化了）solution1：热结构耦合分析可以采用右•接方法，即选用耦合单元solid5,它同时包括了温度和位移白由度，同时附加了电磁特性，这里我们不关注。我们采用瞬态分析的方法，一开始把焊接单元全部杀死，这里的杀死意味着单元的刚度等属性被赋予一个小量（默认的是le-6）0然麻随着焊接过程依次激活单元,加载温度,认为悍料温度为1500度，也是其材料的参考温度。焊完一圈麻冷却降温！但遗憾的
2、是这样的分析是基于线弹性理论的，solid5单元屮无法进行報性分析。以下是直接法的apdl程序：/title,WeldAnalysisby^ElementBirthandDeathM/FILNAME,welding_direct,1/prep7et,l,5!1号材料是钢!2号材料是焊料MG-51T!假设他们的刚度随温度是变化的。MPTEMP,1,0MPTEMP,2,500MPTEMP,3,1000MPTEMP,4,1500MPDATA,EX,l,,2E+（）05MPDATA,EX,l,,1.5e5MPDATA,EX,l„7e4MPDATA,EX,l,,le4MPDATA,PRXY,
3、l,,0.3MPDATA,PRXY,l,,0.3MPDATA,PRXY,l,,0.3MPDATA,PRXY,l,,0.3MPDATA,EX,2,,2E+005MPDATA,EX,2„1.5e5MPDATA,EX,2„7e4MPDATA,EX,2„le4MPDATA,PRXY,2,,0.3MPDATA,PRXY,2,,0.3MPDATA,PRXY,2,,0.3MPDATA,PRXY,2,,0.3!材料密度（假设为常值）mp,dens,l,8e-3rnp,dens,2,8e-3!热膨胀系数（假设为常值）mp,alpxj,l.2e-5mp,alpx,2J.2e-5!热传导系数（假设为常值
4、）mp,kxx,1,0.03mp,kxx,2,0.03!比热（假设为常值）mpc1,0.6mp,c,2,0.6!由于该5号单元还有磁H由度，此处假设一磁特性，但并不影响我们所关心的结果mpjnurxJJe-10mp,murx,2,le-10mp,rsvx,l,le-10mp,rsvx,2,le-10!假设焊料焊上去麻的初始温度是1500°Cmp,reft,1,25mp,reft,2,1500!下面建立儿何模型csys,0k,1,7,0,0k,2,7,0,4k,3,56,0,4k,4,56,0,11k,5,60,0,11k,6,60,0,8k,7,60,0,0k,8,63,0,8k,
5、9,63,0,11k,12,62,0,8kJ3,62,0,-12k,14,60,0,-12k,10,0,0,0k,ll,0,0,la,1,2,3,4,5,6,7a,5,6,&9a,7,6,12,13」4ET,2,SHELL63!划分单元esize,2type,2matJamesh,l,3,2mat,2amesh,2TYPE,1EXTOPT,ESIZE,25,0,EXTOPT,ACLEAR,1EXTOPT,ATTR,1,1,1vrotat,aHM,„JO,l1„2cswplaj1,1aclear,alletdel,2asel,s,,,18,33,15da,all,ux,0da,all
6、,uy,0da,all,uz,0/soluantype,4!瞬态分析trnopt,fullOUTRES,ALL,1I!假设模型表面施加对流，冷却时也是！allselsfa.allj，conv,1.1e-5,25!由于第2,5体是焊料所在区域，因此，首先将该区域的单元“死”掉vsel,s”,2,5,3eslv*get,wnum,elem„count*dim,eorder”wnum*dim,ne„wnum*dim,ney^wnum*get,emax,elem„numjYiax*get,emin,elemMnum,minmine=0!下面的DO循环用于将焊料区的单元按柱坐标csyslly
7、方向排序，以便后面!模拟焊料逐步“生长“的过程ii=0*do,i,emin,emax*if,esel(i),eq丄then••••■11=11+1ne(ii)=i*endif*enddooJJ.wnum*get,ney(i),elem,ne(i),cent,y水enddo*dojl,l,wnumminy=1000*doj丄wnum*if,ney(i)Jt,miny>thenminy二ney(i)mine=ne(i)*•J=»*else*if5ney(i),eq,min
