ansys焊接温度场模拟

上课睡觉

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

复制代码

俠菩提

請問在workbench中是否可以使用以上的方法來做焊接題目

wjdcumt

俠菩提 发表于 2018-11-25 16:23
請問在workbench中是否可以使用以上的方法來做焊接題目

workbench可以直接使用移动热源进行操作