设为首页收藏本站联系我们
切换到宽版

声振论坛

 找回密码
 我要加入

QQ登录

只需一步，快速开始

声振论坛»论坛 › 工程软件 › Ansys讨论区 › 一个完整的三维移动高斯分布表面热源焊接例子
12下一页
返回列表 发新帖
查看: 6549|回复: 15

[综合讨论] 一个完整的三维移动高斯分布表面热源焊接例子

[复制链接]
AaronSpark
发表于 2006-5-16 07:33 | 显示全部楼层 |阅读模式

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

您需要 登录 才可以下载或查看，没有账号？我要加入

x
  1. /prep7
  2. /UNITS,si

  4. et,1,solid70
  5. mp,c,1,544.28
  6. mp,kxx,1,15.06
  7. mp,dens,1,4500

  9. block,0,30,0,15,0,3

  11. esize,1
  12. vmesh,all
  13. FINISH

  15. /SOLU
  16. !PARAMETER TABLE TO INPUT HEAT FLUX AS A MOVING SOURCE
  17. !SPEED OF HEAT SOURCE IN X DIRECTION v=(200/60)m/s
  18. !GAUSS VARIATION ASSUMED IN X & Y DIRECTION AS q=qmax/EXP(K*R**2)
  19. !TOTAL HEAT INPUT Q=3.14*qmax/K
  20. !DIMENSIONING TABLE NAMED "FLUX2"
  21. v=1
  22. L=30
  23. NUM_TIME=L/v+1
  24. NUM_X=L/v+1

  26. *DIM,FLUX2,TABLE,NUM_TIME,NUM_X,7,TIME,X,Y

  28. *DO,K,1,7,1

  30. *DO,I,1,NUM_TIME,1
  31. FLUX2(I,0,K)=(I-1)
  32. *ENDDO

  34. *DO,J,1,NUM_X,1
  35. FLUX2(0,J,K)=(J*v-v)
  36. *ENDDO

  38. Q=1000
  39. KAPA=1/4
  40. qmax=((Q*KAPA)/3.14159)
  41. PERCENT=0.98
  42. RADIUS=ABS(SQRT(1/KAPA*LOG(1/(1-PERCENT))))

  44. !N=[RADIUS/v]

  46. *DO,ANTIME,1,NUM_TIME,1
  47. *DO,N,0,15,1
  48. XDIST=(N*v)
  49. *IF,(ANTIME-N),GT,0,THEN
  50. *IF,SQRT((XDIST)**2+((RADIUS/6)*(K-1))**2),LE,RADIUS,THEN
  51. !
  52. GAUSS=qmax/EXP(KAPA*((XDIST)**2+((RADIUS/6)*(K-1))**2))
  53. FLUX2(ANTIME,ANTIME-N,K)=GAUSS*10e5
  54. !
  55. *ELSE
  56. *ENDIF
  57. *ENDIF
  58. !
  59. *IF,(ANTIME+N),LE,NUM_TIME,THEN
  60. *IF,SQRT((XDIST)**2+((RADIUS/6)*(K-1))**2),LE,RADIUS,THEN
  61. !
  62. GAUSS=qmax/EXP(KAPA*((XDIST)**2+((RADIUS/6)*(K-1))**2))
  63. FLUX2(ANTIME,ANTIME+N,K)=GAUSS*10e5
  64. !
  65. *ELSE
  66. *ENDIF
  67. *ENDIF
  68. *ENDDO
  69. *ENDDO
  70. !
  71. !INPUT Y DIRECTION DIVISIONS
  72. !---------------------------
  73. FLUX2(0,0,K)=(K-1)*(RADIUS/6)
  74. *ENDDO

  76. antype,trans
  77. tunif,20
  78. kbc,1

  80. *SET,delt,1
  81. timint,on
  82. autots,off

  84. *do,tm,1e-6,30,0.5
  85. /VIEW, 1 ,1,1,1
  86. /ANG, 1
  87. /REP,FAST
  88. FLST,2,1,5,ORDE,1
  89. FITEM,2,2
  90. /GO
  91. !*
  92. TIME,tm
  93. AUTOTS,0
  94. DELTIM, , , ,1
  95. KBC,0
  96. !*
  97. TSRES,ERASE
  98. FLST,2,1,5,ORDE,1
  99. FITEM,2,2
  100. /GO
  101. !*
  102. !*
  103. SFA,P51X,1,HFLUX, %FLUX2%
  104. /STATUS,SOLU
  105. SOLVE
  106. OUTRES,ALL,ALL,
  107. /PSF,hflux,,2

  109. /replot

  111. *enddo
  112. OUTRES,ALL,ALL,

  114. FINISH
复制代码

回复
分享到：

使用道具 举报

AaronSpark
 楼主| 发表于 2006-5-16 07:34 | 显示全部楼层
这个例子的看点:

1.建模.

2.死活单元.

3.相变.

4.局部坐标系.


  1. /PREP7
  2. !*********************
  3. ET,1,SOLID70 !thermal solid

  5. MP,KXX,1,.5e-3 !weld on right
  6. MP,C,1,.2
  7. MP,DENS,1,.2833
  8. MPTEMP,1,0,2643,2750,2875,3000 !enthalpy as function of time
  9. MPDATA,ENTH,1,1,0,128.1,163.8,174.2,184.6

  11. MP,KXX,2,.5e-3 !solid material
  12. MP,C,2,.2
  13. MP,DENS,2,.2833

  15. MP,KXX,3,0.5e-3 !weld on left
  16. MP,C,3,.2 !for t=[0,1000]: only right weld, left weld has no therm conduct
  17. MP,DENS,3,.2833 !for t>1000: left weld is activated

  19. W=1.20
  20. TL=.34
  21. TWELD=.2 ! weld thickness
  22. TPIVOT=TL ! thickness of support bar
  23. TPHALF=TPIVOT/2
  24. WIDTH=W ! width of support bar
  25. WHALF=WIDTH/2
  26. LCYL=2.5 ! length of cylinder
  27. LHALF=LCYL/2
  28. RID=.50 ! inside radius of cylinder
  29. ROD=.75 ! outside radius of cylinder
  30. *AFUN,DEG
  31. ANG1=50-ASIN(TPHALF/ROD)-2*ASIN(TWELD/(2*ROD))
  32. ANG2=50-ASIN(TPHALF/ROD)
  33. ANG3=50+ASIN(TPHALF/ROD)
  34. ANG4=50+ASIN(TPHALF/ROD)+2*ASIN(TWELD/(2*ROD))
  35. WPST,-1 ,1 ,-10 ,10 ,.003 ,3

  37. /PNUM,VOLU,1
  38. /VIEW,1,-.2,.2,1
  39. CYLIND,RID,ROD,-LHALF,-WHALF,-90,ANG1
  40. CYLIND,RID,ROD,-WHALF,+WHALF,-90,ANG1
  41. CYLIND,RID,ROD,+WHALF,+LHALF,-90,ANG1
  42. CYLIND,RID,ROD,-LHALF,-WHALF,ANG1,ANG2
  43. CYLIND,RID,ROD,-WHALF,+WHALF,ANG1,ANG2
  44. CYLIND,RID,ROD,+WHALF,+LHALF,ANG1,ANG2
  45. CYLIND,RID,ROD,-LHALF,-WHALF,ANG2,ANG3
  46. CYLIND,RID,ROD,-WHALF,+WHALF,ANG2,ANG3
  47. CYLIND,RID,ROD,+WHALF,+LHALF,ANG2,ANG3
  48. CYLIND,RID,ROD,-LHALF,-WHALF,ANG3,ANG4
  49. CYLIND,RID,ROD,-WHALF,+WHALF,ANG3,ANG4
  50. CYLIND,RID,ROD,+WHALF,+LHALF,ANG3,ANG4
  51. CYLIND,RID,ROD,-LHALF,-WHALF,ANG4,90
  52. CYLIND,RID,ROD,-WHALF,+WHALF,ANG4,90
  53. CYLIND,RID,ROD,+WHALF,+LHALF,ANG4,90

  55. /TYPE,1,4
  56. WPROT,50 ! vol 16 (used to trim welds)
  57. BLOCK,.6,rod*cos(asin(tphalf/rod))+tweld,-tphalf,tphalf,-whalf,whalf
  58. CSWPLA,11,0
  59. KSEL,S,LOC,X,.75*COS(ASIN(TPHALF/ROD))
  60. KSEL,R,LOC,Y,-TPHALF
  61. KSEL,R,LOC,Z,LHALF
  62. *GET,KNUM1,KP,,NUM,MAX
  63. KSEL,S,LOC,X,.75*COS(ASIN(TPHALF/ROD))
  64. KSEL,R,LOC,Y,TPHALF
  65. KSEL,R,LOC,Z,LHALF
  66. *GET,KNUM2,KP,,NUM,MAX
  67. KSEL,ALL
  68. WPRO,-50
  69. CSYS,0
  70. LOCAL,12,0,KX(KNUM1),KY(KNUM1)
  71. CSYS,0
  72. LOCAL,13,0,KX(KNUM2),KY(KNUM2)
  73. WPCSYS,-1,12
  74. VSEL,U,,,1,16
  75. CYLIND,,TWELD,-WHALF,WHALF,(ANG1+ANG2)/2-90,70
  76. /DIST,1,.6
  77. /FOCUS,1,.47,.59
  78. WPCSYS,-1,13
  79. CYLIND,,TWELD,-WHALF,WHALF,35,90+(ANG3+ANG4)/2
  80. VSEL,ALL
  81. BOPT,keep,yes,
  82. bopt,version,rv52
  83. VSEL,S,,,5,17,12
  84. VSBV,17,5
  85. VSEL,A,,,16,19,3
  86. VSBV,19,16
  87. VSEL,A,,,11,18,7
  88. VSBV,18,11
  89. VSEL,A,,,16,21,5
  90. VSBV,21,16
  91. VSEL,ALL
  92. WPSTYLE,default
  93. /FOCUS,,.5,.5
  94. /DIST,1,.6
  95. CYLIN,.6,SQRT((ROD+.01)**2+TPHALF**2),-2,+2,10,90
  96. BOPT,keep,no
  97. VSEL,S,,,20,22,2
  98. VSEL,A,,,23
  99. VSEL,ALL
  100. VOVLAP,20,22,23
  101. VDELE,16,21,1,1
  102. VDELE,28,,,1
  103. /TYPE,1,4

  105. CSYS,0
  106. WPCSYS,-1
  107. WPROT,50

  109. NUMSTR,VOLU,29
  110. block,ROD+.01,ROD*COS(ASIN(TPHALF/ROD))+TWELD,-TPHALF,TPHALF,-WHALF,WHALF
  111. NUMSTR,VOLU,30
  112. BLOCK,ROD*COS(ASIN(TPHALF/ROD))+TWELD,3,-TPHALF,TPHALF,-WHALF,WHALF
  113. NUMMRG,KP
  114. VSEL,S,,,24,26,2 ! right hand weld (make mat 1)
  115. ESIZE,.1
  116. VMESH,ALL
  117. MAT,3
  118. vsel,s,,,25,27,2 ! left hand weld (make mat 3 - dead)
  119. vmesh,all
  120. ALLSEL,ALL
  121. VSEL,S,,,24,27,1
  122. VSEL,A,,,30
  123. VSEL,INVE
  124. ESIZE,.1
  125. MAT,2
  126. VMESH,ALL
  127. ESIZE,.2
  128. LSEL,S,,,198,199
  129. LESIZE,ALL,,,8,5,1
  130. LSEL,S,,,196,201,5
  131. LESIZE,ALL,,,8,.2,1
  132. ALLSEL,ALL
  133. VMESH,30
  134. NUMCMP,all
  135. WPSTYLE,,,,,,,,0
  136. esel,s,mat,,3
  137. ekill,all
  138. allsel,all
  139. SAVE,mesh,DB
  140. FINI
  141. !
  142. !
  143. !*********************
  144. /SOLU
  145. !*********************
  146. ANTYPE,TRANS

  148. !Initialization of temperatures
  149. !----------------------------
  150. TIMINT,OFF
  151. ESEL,S,MAT,,1 !weld on right
  152. NSLE
  153. D,ALL,TEMP,3000 !To(weld)=3000
  154. NSEL,INVE
  155. D,ALL,TEMP,70 !To(steel)=70
  156. TIME,1
  157. KBC,0
  158. NSUBST,2,2,2
  159. OUTRES,BASIC,2
  160. ALLSEL,all
  161. SOLVE
  162. SAVE,therm,db

  164. !Phase change for weld on right (time varies from 1 to 100 and from 100 to 1000)
  165. !----------------------------------------------------------------------------------
  166. OUTRES,NSOL,1
  167. DDELE,ALL,TEMP
  168. TIMINT,ON
  169. TINTP,,,,1
  170. TIME,100 !phase change
  171. DELTIME,1,.5,10
  172. AUTOTS,ON
  173. ASEL,S,EXT
  174. ASEL,U,LOC,X,0
  175. SFA,ALL,,CONV,5E-5,70
  176. ALLSEL,ALL
  177. SOLVE

  179. TIME,1000 !solidification
  180. DELTIME,50,10,500
  181. AUTOTS,ON
  182. SOLVE
  183. FINI

  185. !Phase change for weld on right (time varies from 1e3 to 1100 and from 1100 to 2000)
  186. !---------------------------------------------------------------------------------------
  187. /PREP7
  188. EALIVE,ALL
  189. ALLSEL,all
  190. FINI

  192. /SOLU
  193. ANTYPE,TRANS,REST
  194. vsel,s,,,17
  195. vsel,a,,,19
  196. nslv,,1
  197. d,all,temp,3000
  198. time,1001
  199. nsubst,2,2,2
  200. allsel,all
  201. solve

  203. ANTYPE,TRANS,REST
  204. DDELE,ALL,TEMP
  205. TIME,1100
  206. DELTIME,1,.5,10
  207. SOLVE
  208. TIME,2000
  209. DELTIME,100,10,200
  210. SOLVE
  211. FINI
复制代码

评分

1

查看全部评分

回复 支持 反对

使用道具 举报

lorry
发表于 2006-9-7 21:23 | 显示全部楼层
楼主的例子从何而来？望能提供更多的例子。
回复 支持 反对

使用道具 举报

hxb944210
发表于 2006-9-26 10:12 | 显示全部楼层
例子不错，就是少了一点注释，跟编程一样，没有注释的代码别人总是很难看懂的。建议今后再上传类似例子时能有相关的注释。
回复 支持 反对

使用道具 举报

新元素
发表于 2006-12-19 10:04 | 显示全部楼层
也是
我是新手
现在看命令流有点困难
建议
加上注释
毕竟还有广大的初级爱好者
谢谢
回复 支持 反对

使用道具 举报

datace100
发表于 2008-4-25 17:15 | 显示全部楼层
感谢LZ ！这个例子很好啊
回复 支持 反对

使用道具 举报

xsdcaifu
发表于 2008-4-26 12:47 | 显示全部楼层
例子很好，这是两个例子吗？:@)
回复 支持 反对

使用道具 举报

bavil
发表于 2008-5-9 14:17 | 显示全部楼层
zhi支持，就是看的不是很懂
回复 支持 反对

使用道具 举报

onlyu
发表于 2008-5-9 17:18 | 显示全部楼层
不久我也要做这方面的东西，感谢了先，以后还得多多请教
回复 支持 反对

使用道具 举报

lixiaoyu0323
发表于 2009-12-9 15:46 | 显示全部楼层
非常感谢，好好看咯:@)
回复 支持 反对

使用道具 举报

yaojinnan2009
发表于 2009-12-18 16:23 | 显示全部楼层
谢谢啦。。。呵呵
回复 支持 反对

使用道具 举报

守护蓝天
发表于 2011-1-3 15:23 | 显示全部楼层
下来研究下，现在正在做这方面的工作，急需这样的例子。还有这样的例子吗?
回复 支持 反对

使用道具 举报

lzwlws
发表于 2011-3-9 11:44 | 显示全部楼层
正需要这个，只是得一行一行看，有点麻烦啊
回复 支持 反对

使用道具 举报

lzwlws
发表于 2011-3-9 11:50 | 显示全部楼层
回复 1 # AaronSpark 的帖子

问下楼主，第28行的7，和47行的15怎么确定的？
回复 支持 反对

使用道具 举报

gension
发表于 2011-3-30 13:10 | 显示全部楼层
能讲解一下否？
回复 支持 反对

使用道具 举报

下一页 »
12下一页
返回列表 发新帖
高级模式
B Color Image Link Quote Code Smilies
您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

QQ|小黑屋|Archiver|手机版|联系我们|声振论坛

GMT+8, 2024-9-22 07:08 , Processed in 0.072242 second(s), 19 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

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