平板焊接之后极限拉伸载荷的计算问题，求教

ouyangyu3

1. FINISH
   
2. /CLEAR
   
3. /PREP7
   
4. /TITLE,YLJS
   
5. ET,1,SOLID70
   
6. 
   
7. MPTEMP,1,25,500,1000,1500
   
8. MPDATA,EX,1,1,1.96E11,1.4275E11,0.7575E11,0.0425E11
   
9. MPDATA,EX,2,1,1.96E11,1.4275E11,0.7575E11,0.0425E11
   
10. TB,MISO,1,4
    
11. TBTEMP,25
    
12. TBPT,,0.004,7.84E+8
    
13. TBPT,,0.0045,8.685E+8
    
14. TBPT,,0.005,9.0E+8
    
15. TBPT,,0.009,9.18E+8
    
16. TBPT,,0.1,9.18E+8
    
17. TBTEMP,500
    
18. TBPT,,0.004,5.71E+8
    
19. TBPT,,0.005,5.75E+8
    
20. TBPT,,0.009,5.88E+8
    
21. TBPT,,0.01,5.88E+8
    
22. TBTEMP,1000
    
23. TBPT,,0.004,3.03E+8
    
24. TBPT,,0.005,3.06E+8
    
25. TBPT,,0.009,3.13E+8
    
26. TBPT,,0.01,3.13E+8
    
27. TBTEMP,1500
    
28. TBPT,,0.004,1.7E+007
    
29. TBPT,,0.005,1.7E+007
    
30. TBPT,,0.009,1.8E+007
    
31. TBPT,,0.01,1.8E+007
    
32. TB,MISO,2,4,
    
33. TBTEMP,25
    
34. TBPT,,0.004,7.84E+8
    
35. TBPT,,0.0045,8.685E+8
    
36. TBPT,,0.005,9.0E+8
    
37. TBPT,,0.009,9.18E+8
    
38. TBPT,,0.1,9.18E+8
    
39. TBTEMP,500
    
40. TBPT,,0.004,5.71E+8
    
41. TBPT,,0.005,5.75E+8
    
42. TBPT,,0.009,5.88E+8
    
43. TBPT,,0.01,5.88E+8
    
44. TBTEMP,1000
    
45. TBPT,,0.004,3.03E+8
    
46. TBPT,,0.005,3.06E+8
    
47. TBPT,,0.009,3.13E+8
    
48. TBPT,,0.01,3.13E+8
    
49. TBTEMP,1500
    
50. TBPT,,0.004,1.7E+007
    
51. TBPT,,0.005,1.7E+007
    
52. TBPT,,0.009,1.8E+007
    
53. TBPT,,0.01,1.8E+007
    
54. MP,PRXY,1,0.3
    
55. MP,PRXY,2,0.3
    
56. MP,ALPX,1,1.06e-5
    
57. MP,ALPX,2,1.06e-5
    
58. MP,KXX,1,16.3
    
59. MP,KXX,2,16.3
    
60. MP,C,1,502
    
61. MP,C,2,502
    
62. MP,MURX,1,1
    
63. MP,MURX,2,1
    
64. MP,DENS,1,7800
    
65. MP,DENS,2,7800
    
66. UIMP,1,REFT,25
    
67. UIMP,2,REFT,1300
    
68. 
    
69. hflengh=0.003
    
70. leftlengh=0.4
    
71. rightlengh=0.4
    
72. width=0.25
    
73. depth=0.008
    
74. k,1,0,0,0
    
75. k,2,leftlengh,0,0
    
76. k,3,leftlengh,width,0
    
77. k,4,0,width,0
    
78. k,5,0,0,depth
    
79. k,6,leftlengh,0,depth
    
80. k,7,leftlengh,width,depth
    
81. k,8,0,width,depth
    
82. k,9,leftlengh+hflengh,0,0
    
83. k,10,leftlengh+hflengh+rightlengh,0,0
    
84. k,11,leftlengh+hflengh+rightlengh,width,0
    
85. k,12,leftlengh+hflengh,width,0
    
86. k,13,leftlengh+hflengh,0,depth
    
87. k,14,leftlengh+hflengh+rightlengh,0,depth
    
88. k,15,leftlengh+hflengh+rightlengh,width,depth
    
89. k,16,leftlengh+hflengh,width,depth
    
90. V,5,1,4,8,6,2,3,7
    
91. V,6,2,3,7,13,9,12,16
    
92. V,13,9,12,16,14,10,11,15
    
93. VGLUE,ALL
    
94. NUMCMP,ALL
    
95. 
    
96. LESIZE,4,0.005,,,
    
97. LESIZE,11,0.005,,,
    
98. LESIZE,3,0.004,,,
    
99. TYPE,1
    
100. MAT,1
     
101. MSHAPE,0,3D
     
102. MSHKEY,1
     
103. VMESH,1
     
104. 
     
105. LESIZE,28,0.005,,,
     
106. LESIZE,27,0.005,,,
     
107. LESIZE,26,0.004,,,
     
108. TYPE,1
     
109. MAT,1
     
110. MSHAPE,0,3D
     
111. MSHKEY,1
     
112. VMESH,3
     
113. 
     
114. LESIZE,19,0.003,,,
     
115. TYPE,1
     
116. MAT,2
     
117. MSHAPE,0,3D
     
118. MSHKEY,1
     
119. VMESH,2
     
120. NUMCMP,ALL
     
121. EPLOT
     
122. 
     
123. SFA,1,1,CONV,62.5,25
     
124. SFA,2,1,CONV,62.5,25
     
125. SFA,3,1,CONV,62.5,25
     
126. SFA,4,1,CONV,62.5,25
     
127. SFA,5,1,CONV,62.5,25
     
128. SFA,14,1,CONV,62.5,25
     
129. SFA,16,1,CONV,62.5,25
     
130. SFA,12,1,CONV,62.5,25
     
131. SFA,13,1,CONV,62.5,25
     
132. SFA,15,1,CONV,62.5,25
     
133. OUTRES,ALL,ALL
     
134. 
     
135. VSEL,S,,,2
     
136. ESLV,S
     
137. CM,TOTAL,ELEM
     
138. CMSEL,S,TOTAL
     
139. EKILL,ALL
     
140. ESEL,S,LIVE
     
141. EPLOT
     
142. 
     
143. MXT=1300
     
144. DT09=2.399
     
145. DT01=0.001
     
146. T=0
     
147. tt=0
     
148. 
     
149. /SOL
     
150. ANTYPE,4
     
151. TRNOPT,FULL
     
152. LUMPM,0
     
153. ESEL,S,LIVE
     
154. OUTRES,ALL,ALL
     
155. IC,ALL,TEMP,25
     
156. ALLSEL,ALL
     
157. 
     
158. *DO,I,1,50
     
159. ESEL,S,,,16000+I*2
     
160. EALIVE,ALL
     
161. ESEL,S,LIVE
     
162. EPLOT
     
163. T=T+DT09
     
164. TIME,T
     
165. NSUBST,1
     
166. KBC=1
     
167. ESEL,S,,,16000+I*2
     
168. NSLE,S,ALL
     
169. D,ALL,TEMP,MXT
     
170. ALLSEL,ALL
     
171. SOLVE
     
172. T=T+DT01
     
173. TIME,T
     
174. NSUBST,1
     
175. KBC=1
     
176. ESEL,S,,,16000+I*2
     
177. NSLE,S,ALL
     
178. DDELE,ALL,TEMP,MAXT
     
179. ALLSEL,ALL
     
180. SOLVE
     
181. *ENDDO
     
182. 
     
183. *DO,I,1,50
     
184. ESEL,S,,,15999+I*2
     
185. EALIVE,ALL
     
186. ESEL,S,LIVE
     
187. EPLOT
     
188. T=T+DT09
     
189. TIME,T
     
190. NSUBST,1
     
191. KBC=1
     
192. ESEL,S,,,15999+I*2
     
193. NSLE,S,ALL
     
194. D,ALL,TEMP,MXT
     
195. ALLSEL,ALL
     
196. SOLVE
     
197. T=T+DT01
     
198. TIME,T
     
199. NSUBST,1
     
200. KBC=1
     
201. ESEL,S,,,15999+I*2
     
202. NSLE,S,ALL
     
203. DDELE,ALL,TEMP,MAXT
     
204. ALLSEL,ALL
     
205. SOLVE
     
206. *ENDDO
     
207. 
     
208. TIME,1200
     
209. DELTIM,20
     
210. SOLVE
     
211. TIME,3600
     
212. DELTIM,40
     
213. SOLVE
     
214. TIME,7200
     
215. DELTIM,90
     
216. SOLVE
     
217. SAVE,'TEMP','db','F:\ANSYS\REAL'
     
218. 
     
219. /PREP7
     
220. ETCHG,TTS
     
221. EPLOT
     
222. ALLSEL,ALL
     
223. DL,1, ,ALL
     
224. DL,3, ,ALL
     
225. DL,22, ,ALL
     
226. DL,26, ,ALL
     
227. DT24=2.4
     
228. DT20=20
     
229. DT40=40
     
230. DT90=90
     
231. TT=0
     
232. 
     
233. /SOL
     
234. ANTYPE,0
     
235. ISWRITE,ON
     
236. ALLSEL,ALL
     
237. EPLOT
     
238. CMSEL,S,TOTAL
     
239. EKILL,ALL
     
240. ESEL,S,LIVE
     
241. EPLOT
     
242. 
     
243. *DO,I,1,50
     
244. ESEL,S,,,16000+I*2
     
245. EALIVE,ALL
     
246. ESEL,S,LIVE
     
247. EPLOT
     
248. TT=TT+DT24
     
249. ESEL,ALL
     
250. TIME,TT
     
251. NSUBST,1
     
252. LDREAD,TEMP,,, TT, ,'temp','rth','F:\ANSYS\'
     
253. SOLVE
     
254. *ENDDO
     
255. 
     
256. *DO,I,1,50
     
257. ESEL,S,,,15999+I*2
     
258. EALIVE,ALL
     
259. ESEL,S,LIVE
     
260. EPLOT
     
261. TT=TT+DT24
     
262. ESEL,ALL
     
263. TIME,TT
     
264. NSUBST,1
     
265. LDREAD,TEMP,,, TT, ,'temp','rth','F:\ANSYS\'
     
266. SOLVE
     
267. *ENDDO
     
268. 
     
269. *DO,I,1,60
     
270. TT=TT+DT20
     
271. TIME,TT
     
272. NSUBST,1
     
273. LDREAD,TEMP,,, TT, ,'temp','rth','F:\ANSYS\'
     
274. SOLVE
     
275. *ENDDO
     
276. *DO,I,1,90
     
277. TT=TT+DT40
     
278. TIME,TT
     
279. NSUBST,1
     
280. LDREAD,TEMP,,, TT, ,'temp','rth','F:\ANSYS\'
     
281. SOLVE
     
282. *ENDDO
     
283. *DO,I,1,80
     
284. TT=TT+DT90
     
285. TIME,TT
     
286. NSUBST,1
     
287. LDREAD,TEMP,,, TT, ,'temp','rth','F:\ANSYS\'
     
288. SOLVE
     
289. *ENDDO
     
290. SAVE,'result2','db','F:\ansys\real\'
     
291. 
     
292. 
     
293. 
     
294. *************************************************************极限强度计算
     
295. 
     
296. 
     
297. 
     
298. /SOLU
     
299. ANTYPE,0
     
300. ALLSEL,ALL
     
301. DDELE,ALL,ALL
     
302. alls
     
303. nsel,s,loc,X,0
     
304. d,all,all
     
305. alls
     
306. nsel,s,loc,X,0.803
     
307. cp,1,ux,all
     
308. *get,nodenum,node,0,count
     
309. f,all,fX,5000000/nodenum
     
310. D,ALL,UY,,,,,UZ
     
311. alls
     
312. RESCONTROL,DEFINE
     
313. NLGEOM,on
     
314. arclen,on
     
315. NSUBST,100
     
316. outres,all,all
     
317. SOLVE

复制代码

在平板对接焊接，生死单元模拟之后，运用RESCONTROL,进行重启动分析，计算焊接残余应力时候的边界条件，删除与不删除这两种情况下，计算的极限拉伸载荷。都很大，跟在直接计算极限拉伸载荷（是1840KN），相差了几倍的关系，想请教下高手，这个如何解决，是我的命令错了么？
如果运用初始应力施加，ISWRITE输出应力文件，然后ISFILE读入应力值加入模型中计算的极限载荷值跟1840KN接近，但是运用LIST看了下施加的应力场，跟焊接得到的残余应力场有出入，这是怎么回事？计算时候说的是：由于材料定义的是非线性，但是应力场加在了弹性材料上，这个是怎么回事？是不是计算的值不准确了？帮忙解答下。谢谢各位了！

ouyangyu3

自己先顶下，各位大侠，多说说自己的想法啊，小弟叩谢了啊！急着呢