杭州锐达数字技术有限公司
查看: 557|回复: 2

[CFD及热分析] ansys焊接温度场模拟

[复制链接]
发表于 2015-9-29 18:17 | 显示全部楼层 |阅读模式

马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。

您需要 登录 才可以下载或查看,没有帐号?我要加入

x
  1. /UNITS,SI
  2. /CONFIG,NRES,10000
  3. /PREP7
  4. ET,1,SOLID70    !单元类型,1表示此单元编号为1
  5. MP,DENS,1,7930   !定义材料密度   
  6. MP,C,1,502       !定义材料比热
  7. mptemp,1,20,100,500
  8. mpdata,kxx,1,1,12.1,16.3,21.4   !kxx热导系数
  9. /VIEW,1,1,1,1
  10. !*************定义基板尺寸***********
  11. lx=0.1
  12. lz=0.1
  13. ly=0.006    !基板大小
  14. !*************定义焊接参数***********
  15. layer=7                              !总的焊接道数
  16. lap=1/3                          !两道之间搭接量
  17. weld=0.01                        !每道焊道宽度
  18. wweld=layer*weld-(layer-1)*lap*weld      !焊道的总宽度
  19. wheight=0.004                    !焊道的高度
  20. remain=0.01
  21. LSIZE=0.002   
  22. V=0.008          !速度
  23. pi=3.1415926
  24. U=25
  25. I=180
  26. Q=U*I         !电源有效功率(在此,可以对Q通过UI来调整)
  27. R=0.006         !电源有效热半径
  28. effect=0.8
  29. Qmax=effect*Q/(pi*R*R) !中心处最大热流
  30. tm=(lx-2*remain)/V         !一道的总时间
  31. tinc=LSIZE/V    !每小段的时间间隔
  32. !*****************建立模型,生成网格**************
  33. block,0,lx,0,ly,0,lz   !创建一个块体积,可以根据基板尺寸调整坐标
  34. wpoff,remain,ly,(lz-wweld)/2  !随着坐标系轴的变化平移工作平面
  35. block,0,lx-2*remain,0,wheight,0,wweld
  36. vsbw,all,,delete   !对block布尔运算
  37. wpoff,,,wweld
  38. vsbw,all,,delete
  39. wprota,90,90    !转动工作面
  40. vsbw,all,,delete
  41. wpoff,,,lx-2*remain
  42. vsbw,all,,delete
  43. save
  44. vglue,all    !粘合所有的几何体
  45. lsel,s,,,9,12,1
  46. lsel,a,,,25,26,1
  47. lsel,a,,,30
  48. lsel,a,,,32
  49. lsel,a,,,38
  50. lsel,a,,,45
  51. lsel,a,,,59
  52. lsel,a,,,73
  53. lesize,all,0.005,,,0.2,,,,1
  54. allsel,all
  55. lsel,s,,,13,16,1
  56. lsel,a,,,19,21,1
  57. lsel,a,,,23,24,1
  58. lsel,a,,,34
  59. lsel,a,,,41,44,1
  60. lsel,a,,,63,64,1
  61. lsel,a,,,77,78,1
  62. lsel,a,,,81,82,1
  63. lesize,all,0.001,,,1,,,,1
  64. allsel,all
  65. lsel,s,,,1
  66. lsel,a,,,3,6,1
  67. lsel,a,,,8
  68. lsel,a,,,27
  69. lsel,a,,,29
  70. lsel,a,,,31
  71. lsel,a,,,33
  72. lsel,a,,,35
  73. lsel,a,,,37
  74. lsel,a,,,40
  75. lsel,a,,,50
  76. lsel,a,,,53
  77. lsel,a,,,55
  78. lsel,a,,,62
  79. lsel,a,,,67
  80. lsel,a,,,69
  81. lsel,a,,,85
  82. lesize,all,0.002,,,1,,,,1
  83. allsel,all
  84. lsel,s,,,46
  85. lsel,a,,,86
  86. lsel,a,,,88
  87. lesize,all,0.001,,,5,,,,1
  88. allsel,all
  89. vsweep,all
  90. /PNUM,DEFA
  91. EPLOT
  92. fini
  93. !*******************进入后处理器*********************
  94. /solu
  95. antype,trans,new
  96. !nlgeom,on
  97. timint,0,struct
  98. timint,1,therm
  99. timint,0,mag
  100. timint,0,elect
  101. tref,25
  102. nropt,auto
  103. autos,on
  104. kbc,0
  105. pred,on
  106. lnsrch,on
  107. !************杀死焊缝去单元***********
  108. nsel,s,loc,y,ly,ly+wheight
  109. esln,s,1
  110. ekill,all
  111. allsel,all
  112. esel,s,live
  113. eplot
  114. !**********施加对流载荷*************
  115. esel,s,live
  116. nsel,s,loc,z,0
  117. nsel,a,loc,z,lz
  118. nsel,a,loc,x,0
  119. nsel,a,loc,x,lx
  120. sf,all,conv,10,25
  121. !************定义数组维数**********
  122. MAX_TIME=1+tm/tinc
  123. !MAX_X=1+(lx-2*remain)/LSIZE
  124. MAX_X=1+lx/LSIZE
  125. MAX_Z=1+lz/LSIZE
  126. !*************定义table数组
  127. ******************************************
  128. *DIM,heatflux,TABLE,MAX_X,MAX_Z,MAX_TIME,X,Z,TIME
  129. *do,n,1,layer,1
  130. m=mod(n,2)
  131. *if,m,eq,1,then
  132. !n=1
  133.         *DO,K,1,MAX_TIME,1
  134.         *DO,I,1,MAX_X,1
  135.         heatflux(I,0,K)=(I-1)*LSIZE
  136.         *ENDDO
  137.         *DO,J,1,MAX_Z,1
  138.         heatflux(0,J,K)=(J-1)*LSIZE
  139.         *ENDDO
  140.         *DO,I,1,MAX_X,1
  141.         *DO,J,1,MAX_Z,1
  142.         XCENTER=remain+V*(K-1)*tinc
  143.         ZCENTER=2*(n-1)*weld/3+(lz-wweld)/2+weld/2
  144.         DISTANCE=SQRT(((I-1)*LSIZE-XCENTER)**2+((J-1)*LSIZE-
  145. ZCENTER)**2)
  146.         *IF,DISTANCE,LE,R,THEN
  147.         heatflux(I,J,K)=Qmax/EXP(2*DISTANCE**2/(R**2))
  148.         *ELSE
  149.         heatflux(I,J,K)=0
  150.         *ENDIF
  151.         *ENDDO
  152.         *ENDDO
  153.         heatflux(0,0,K)=(K-1)*tinc
  154. !**********激活单元*********
  155.         esel,s,live
  156.         nsel,s,loc,x,xcenter-0.002,xcenter+0.001
  157.         nsel,r,loc,z,zcenter-weld/2,zcenter+weld/2
  158.         nsel,r,loc,y,ly,ly+wheight
  159.         esln,s,0
  160.         ealive,all
  161.         allsel,all
  162. !**********施加热流载荷**********
  163.         esel,s,live
  164.         nsel,s,loc,Y,ly
  165.         sf,all,hflux,%heatflux%
  166.         allsel,all
  167. !************求解*****************
  168.         time,(n-1)*(tm+tinc)+k*lsize/v
  169.         nsubst,2
  170.         SOLVE
  171. !*************删除热流载荷***********
  172.         nsel,s,loc,y,ly
  173.         sfdele,all,hflux
  174.         allsel,all
  175.         esel,s,live
  176.         eplot
  177.         OUTRES,ALL,ALL,
  178.         /PSF,HFLUX,,2
  179.         /REPLOT
  180.         OUTRES,ALL,ALL,
  181.         *ENDDO
  182.         
  183. *else
  184. !n=2
  185.         *DO,K,1,MAX_TIME,1
  186.         *DO,I,1,MAX_X,1
  187.         heatflux(I,0,K)=(I-1)*LSIZE
  188.         *ENDDO
  189.         *DO,J,1,MAX_Z,1
  190.         heatflux(0,J,K)=(J-1)*LSIZE
  191.         *ENDDO
  192.         *DO,I,1,MAX_X,1
  193.         *DO,J,1,MAX_Z,1
  194.         XCENTER=lx-remain-V*(K-1)*tinc
  195.         ZCENTER=2*(n-1)*weld/3+(lz-wweld)/2+weld/2
  196.         DISTANCE=SQRT(((I-1)*LSIZE-XCENTER)**2+((J-1)*LSIZE-
  197. ZCENTER)**2)
  198.         *IF,DISTANCE,LE,R,THEN
  199.         heatflux(I,J,K)=Qmax/EXP(2*DISTANCE**2/(R**2))
  200.         *ELSE
  201.         heatflux(I,J,K)=0
  202.         *ENDIF
  203.         *ENDDO
  204.         *ENDDO
  205.         heatflux(0,0,K)=(K-1)*tinc
  206. !**********激活单元*********
  207.         esel,s,live
  208.         nsel,s,loc,x,xcenter-0.001,xcenter+0.002
  209.         nsel,r,loc,z,zcenter-weld/2,zcenter+weld/2
  210.         nsel,r,loc,y,ly,ly+wheight
  211.         esln,s,0
  212.         ealive,all
  213.         allsel,all
  214. !**********施加热流载荷**********
  215.         esel,s,live
  216.         nsel,s,loc,Y,ly
  217.         sf,all,hflux,%heatflux%
  218.         allsel,all
  219. !************求解*****************
  220.         time,(n-1)*(tm+tinc)+k*lsize/v
  221.         nsubst,2
  222.         SOLVE
  223. !*************删除热流载荷***********
  224.         nsel,s,loc,y,ly
  225.         sfdele,all,hflux
  226.         allsel,all
  227.         esel,s,live
  228.         eplot
  229.         OUTRES,ALL,ALL,
  230.         /PSF,HFLUX,,2
  231.         /REPLOT
  232.         OUTRES,ALL,ALL,
  233.         *ENDDO
  234. *endif
  235. *enddo
  236. FINISH
复制代码


回复
分享到:

使用道具 举报

发表于 2018-11-25 16:23 | 显示全部楼层
請問在workbench中是否可以使用以上的方法來做焊接題目

点评

workbench可以直接使用移动热源进行操作  详情 回复 发表于 2018-12-3 15:44
发表于 2018-12-3 15:44 | 显示全部楼层
俠菩提 发表于 2018-11-25 16:23
請問在workbench中是否可以使用以上的方法來做焊接題目

workbench可以直接使用移动热源进行操作
您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

快速回复 返回顶部 返回列表