经典的双锥密封分析，不过还是有点问题，很急，望各位斑竹大侠看看到底哪有问题

cosmoszzz

由于还是有点问题，所以这次只贴命令流，等我问题解决了，我再好好整理一下帖子，下调权限分享给论坛的兄弟。
双锥密封结构分析（没有问题）见附录1
由于我需要加温度场做温度和结构耦合分析，采用间接法，稳态耦合也已经没有问题，见附录2
无奈现在还要分析瞬态温度应力，耦合时候发现瞬态时间0.1s没有问题，时间大于1E5没有问题，时间(50,1E5)怎么都不收敛，命令流基本一样，除了求解瞬态温度的时间不一样，哪位好心的高手能帮助一下，真的很急。
有问题的瞬态耦合命令流见附录3。

apprent

楼主的附录1：

Large-sized Double-cone model_1.rar (2.47 KB)
下载次数:15
2008-11-9 13:13




阅读权限: 10
没有问题的双锥密封命令流

1. Finish
   
2. /clear
   
3. /filn,Large-sized Double-cone model_1
   
4. /TITLE,Large-sized Double-cone model_1
   
5. 
   
6. 
   
7. !_____________________参数设定
   
8. 
   
9. !!!!!!!!!!!!!!!!!!!!双锥环1尺寸
   
10. szh1_A=153
    
11. szh1_B=41
    
12. szh1_C=84
    
13. szh1_Di=3200
    
14. szh1_alpha=30
    
15. szh1_g=3.2
    
16. szh1_E=10 !双锥环上部间隙
    
17. 
    
18. !!!!!!!!!!!!!!!!!!!!模型结构尺寸
    
19. deltap=600 !平盖总高
    
20. delta=200
    
21. L=1420
    
22. theta=360/32
    
23. 
    
24. !!!!!!!!!!!!!!!!!!!!螺栓尺寸
    
25. LS_D=130 !螺栓直径
    
26. LSD_D=140 !平盖上螺栓洞直径
    
27. 
    
28. 
    
29. nut_h=80 !螺母高
    
30. nut_r=95 !螺母外切半径
    
31. 
    
32. forc_bolt=526760
    
33. Pi=12.6
    
34. 
    
35. 
    
36. !_____________________定义单元及材料属性
    
37. /prep7
    
38. et,1,solid45
    
39. mp,ex,1,1.98E5 !双锥环35钢材料参数
    
40. mp,nuxy,1,0.3
    
41. 
    
42. mp,ex,2,1.98E5 !顶底盖和筒体材料16Mn参数
    
43. mp,nuxy,2,0.3
    
44. 
    
45. mp,ex,3,1.98E5 !螺栓材料40CrNiMo参数
    
46. mp,nuxy,3,0.3
    
47. 
    
48. mp,mu,4,0.3 !定义摩擦系数
    
49. 
    
50. 
    
51. 
    
52. 
    
53. 
    
54. !________________________________________创建模型
    
55. !!!!!!!!!!!!!!!!!!!!创建双锥环
    
56. k,1001,0,0,0
    
57. k,1002,0,1,0
    
58. *afun,deg !必须把弧度换成角度
    
59. k,1,szh1_Di/2,0
    
60. k,2,kx(1)+(szh1_A-szh1_c)/2*tan(szh1_alpha),ky(1)+(szh1_A-szh1_c)/2
    
61. k,3,kx(2),ky(2)+szh1_c
    
62. k,4,kx(1),ky(1)+szh1_A
    
63. k,5,kx(3)-szh1_B,ky(4)
    
64. k,6,kx(5),ky(1)
    
65. a,1,2,3,4,5,6
    
66. 
    
67. !!!!!!!!!!!!!!!!!!!!创建双锥环间隙
    
68. k,7,kx(6)-szh1_g,ky(6)
    
69. k,8,kx(7),ky(5)+szh1_E
    
70. k,9,kx(4),ky(8)
    
71. a,4,5,6,7,8,9
    
72. 
    
73. blc4,0,0,szh1_Di/2+delta,deltap
    
74. blc4,kx(1),ky(1),delta,-L+ky(2)
    
75. 
    
76. aptn,all
    
77. 
    
78. kwpave,8
    
79. wprota,,90
    
80. asbw,6
    
81. kwpave,5
    
82. asbw,3
    
83. asbw,7
    
84. kwpave,3
    
85. asbw,9
    
86. kwpave,2
    
87. asbw,10
    
88. 
    
89. wprota,,,90
    
90. kwpave,9
    
91. asbw,4
    
92. kwpave,8
    
93. asbw,12
    
94. 
    
95. alls
    
96. adele,1,2,,1
    
97. adele,11,,,1
    
98. 
    
99. kl,33,0.5 !生成kp5
    
100. kwpave,5
     
101. asbw,10
     
102. asbw,3
     
103. asbw,7
     
104. 
     
105. lglue,all
     
106. 
     
107. alls
     
108. vrotat,all,,,,,,1001,1002,theta !平盖与上锥面接触号22，筒体与下锥面接触面号43
     
109. 
     
110. k,2001,szh1_Di/2,0
     
111. k,2002,kx(2001)+(szh1_A-szh1_c)/2*tan(szh1_alpha),ky(2001)+(szh1_A-szh1_c)/2
     
112. k,2003,kx(2002),ky(2002)+szh1_c
     
113. k,2004,kx(2001),ky(2001)+szh1_A
     
114. k,2005,kx(2003)-szh1_B,ky(2004)
     
115. k,2006,kx(2005),ky(2001)
     
116. a,2001,2002,2003,2004,2005,2006 !形成a58
     
117. 
     
118. vrotat,58,,,,,,1001,1002,theta !双锥环v13,上锥面61,下锥面59
     
119. 
     
120. 
     
121. !!!!!!!!!!!!!!!!!!!!创建平盖螺栓洞
     
122. kl,44,0.5 !生成kp56
     
123. wpstyl,defa
     
124. kwpave,56
     
125. wprota,,90
     
126. 
     
127. cyl4,,,LSD_D/2,,,,ky(56)-szh1_A/2 !ky(56)是active cs下的坐标
     
128. alls
     
129. bopt,numb,off
     
130. vsel,s,volu,,1,3
     
131. vsel,a,volu,,9,11
     
132. vsel,a,volu,,14
     
133. vsbv,all,14 !挖出平盖螺栓孔洞
     
134. 
     
135. alls
     
136. cyl4,,,LS_D/2,,,,ky(56)+L-ky(2)
     
137. vsel,s,volu,,1,5,4
     
138. vsel,a,volu,,8
     
139. vsbv,all,1 !挖出筒体螺栓洞
     
140. 
     
141. 
     
142. cyl4,,,LS_D/2,,,,ky(56)+L-ky(2)+120 !创建螺栓,120是螺栓超出平盖部分长度
     
143. vgen,2,1,,,,120,,,,1
     
144. 
     
145. rpr4,6,,,nut_r,,-nut_h !螺母
     
146. 
     
147. 
     
148. !!!!!!!!!!!!!!!!!!!创建cm
     
149. alls
     
150. vsel,s,volu,,1,3
     
151. vsel,a,volu,,5
     
152. vsel,a,volu,,19,20 !在平盖上映出螺母相关线
     
153. vptn,all
     
154. vsel,s,,,8,10
     
155. vadd,all !螺栓体号v1
     
156. 
     
157. vsel,s,,,23
     
158. cm,luomu,volu
     
159. 
     
160. alls
     
161. vsel,s,,,6,7
     
162. vsel,a,,,15,18
     
163. vsel,a,,,21,22
     
164. vsel,a,,,4,12,8
     
165. cm,pinggai,volu
     
166. 
     
167. alls
     
168. vsel,s,,,11,14,3
     
169. cm,tongti,volu
     
170. 
     
171. alls
     
172. vsel,s,,,13
     
173. cm,szh,volu
     
174. 
     
175. !_____________________划分网格
     
176. wpstyl,defa !wp回原点，很可能在后面做预紧单元时有用
     
177. alls
     
178. mat,1 !双锥环网格
     
179. esize,10
     
180. lesize,108,,,16,,,,,1
     
181. cmsel,s,szh
     
182. vsweep,all
     
183. 
     
184. mat,2
     
185. cmsel,s,pinggai
     
186. esize,20
     
187. lesize,55,,,16,,,,,1
     
188. lesize,3,,,4,,,,,1
     
189. lesize,90,,,16,,,,,1
     
190. lesize,26,,,17,,,,,1
     
191. lesize,74,,,16,,,,,1
     
192. lesize,13,,,10,,,,,1
     
193. lesize,22,,,35,,,,,1
     
194. lesize,75,,,35,,,,,1
     
195. vsweep,15
     
196. vsweep,18
     
197. vsweep,17
     
198. vsweep,16
     
199. vsweep,7
     
200. vsweep,6
     
201. vsweep,12
     
202. vsweep,4
     
203. 
     
204. lesize,5,,,10,,,,,1
     
205. lesize,18,,,10,,,,,1
     
206. lesize,40,,,10,,,,,1
     
207. lesize,38,,,10,,,,,1
     
208. lesize,6,,,3,,,,,1
     
209. lesize,2,,,3,,,,,1
     
210. lesize,45,,,3,,,,,1
     
211. lesize,39,,,3,,,,,1
     
212. lesize,204,,,2,,,,,1
     
213. lesize,202,,,2,,,,,1
     
214. lesize,203,,,2,,,,,1
     
215. lesize,205,,,2,,,,,1
     
216. lesize,168,,,3,,,,,1
     
217. lesize,169,,,3,,,,,1
     
218. lesize,166,,,3,,,,,1
     
219. lesize,171,,,3,,,,,1
     
220. lesize,198,,,3,,,,,1
     
221. lesize,200,,,3,,,,,1
     
222. lesize,150,,,20,,,,,1
     
223. lesize,149,,,20,,,,,1
     
224. lesize,153,,,20,,,,,1
     
225. lesize,147,,,20,,,,,1
     
226. MSHKEY,0
     
227. MSHAPE,1,3d
     
228. esize,12
     
229. vmesh,21
     
230. vmesh,22
     
231. 
     
232. cmsel,s,tongti
     
233. esize,15
     
234. lesize,31,,,4,,,,,1
     
235. lesize,21,,,20,,,,,1
     
236. lesize,82,,,16,,,,,1
     
237. lesize,70,,,16,,,,,1
     
238. vsweep,all
     
239. 
     
240. mat,3
     
241. vsel,s,,,1 !螺栓
     
242. vmesh,all
     
243. 
     
244. cmsel,s,luomu
     
245. vmesh,all
     
246. !!!!!!!!!!!!!!!!!!网格划分完毕
     
247. 
     
248. !___________________________将各部分模型节点压缩，接触面部分节点不能压缩，否则将是刚性连接，不是接触
     
249. !!!!!!!!!!!!!平盖，筒体，螺栓，螺母部分
     
250. alls
     
251. cmsel,s,pinggai
     
252. cmsel,a,tongti
     
253. cmsel,a,luomu
     
254. vsel,a,,,1
     
255. nslv,s,1
     
256. nummrg,node
     
257. 
     
258. 
     
259. !____________________________创建预紧单元
     
260. alls
     
261. psmesh,1,yujin,,volu,1,0,y,szh1_A/2,,,,lsyj
     
262. 
     
263. 
     
264. !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
     
265. !____________________________创建接触
     
266. 
     
267. !!!!!!!!!!!!!!!!上锥面和平盖接触
     
268. alls
     
269. asel,s,area,,61 !上锥面
     
270. cm,szm_up,area
     
271. alls
     
272. asel,s,area,,22 !平盖下
     
273. cm,pg_down,area
     
274. 
     
275. et,3,170
     
276. et,4,174
     
277. keyopt,4,9,0
     
278. keyopt,4,10,2
     
279. 
     
280. mat,4
     
281. r,2
     
282. real,2
     
283. 
     
284. cmsel,s,pg_down !目标面
     
285. type,3
     
286. nsla,s,1
     
287. esln,s,0
     
288. esurf,all
     
289. 
     
290. alls
     
291. cmsel,s,szm_up !接触面
     
292. type,4
     
293. nsla,s,1
     
294. esln,s,0
     
295. esurf,all
     
296. 
     
297. 
     
298. !!!!!!!!!!!!!!!!下锥面和筒体接触
     
299. alls
     
300. asel,s,area,,59 !下锥面
     
301. cm,szm_down,area
     
302. alls
     
303. asel,s,area,,43 !底座
     
304. cm,dz,area
     
305. 
     
306. et,5,170
     
307. et,6,174
     
308. keyopt,6,9,0
     
309. keyopt,6,10,2
     
310. 
     
311. r,3
     
312. real,3
     
313. mat,4
     
314. cmsel,s,dz !目标面
     
315. type,5
     
316. nsla,s,1
     
317. esln,s,0
     
318. esurf,all
     
319. 
     
320. alls
     
321. cmsel,s,szm_down !接触面
     
322. type,6
     
323. nsla,s,1
     
324. esln,s,0
     
325. esurf,all
     
326. 
     
327. !!!!!!!!!!!!!!!!双锥环内表面和平盖接触
     
328. alls
     
329. asel,s,area,,63 !双锥环内表面
     
330. cm,szh_in,area
     
331. alls
     
332. asel,s,area,,36,41,5 !平盖与双锥环接触部分
     
333. cm,pg_szh,area
     
334. 
     
335. et,7,170
     
336. et,8,174
     
337. keyopt,8,10,2
     
338. 
     
339. r,4,,,10,,2
     
340. real,4
     
341. mat,4
     
342. cmsel,s,szh_in !接触面
     
343. type,8
     
344. nsla,s,1
     
345. esln,s,0
     
346. esurf,all
     
347. 
     
348. alls
     
349. cmsel,s,pg_szh !目标面
     
350. type,7
     
351. nsla,s,1
     
352. esln,s,0
     
353. esurf,all
     
354. 
     
355. fini
     
356. 
     
357. 
     
358. !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!接触创建完毕
     
359. /solu
     
360. alls
     
361. lsclear,all
     
362. nropt,full
     
363. pred,on
     
364. lnsrch,on
     
365. antype,0
     
366. autots,1
     
367. nsubst,50,400,1
     
368. outres,all,all
     
369. EQSLV,pcg,1E-5
     
370. NCNV,0,0,0,0,0 !发生错误不中断
     
371. 
     
372. wpstyl,defa
     
373. csys,5
     
374. asel,s,loc,y,0
     
375. asel,a,loc,y,theta
     
376. da,all,symm
     
377. 
     
378. csys,0
     
379. asel,s,loc,y,-L+ky(2)
     
380. nsla,s,1
     
381. d,all,uy
     
382. 
     
383. 
     
384. !_________________________________________加载
     
385. 
     
386. 
     
387. !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!预紧力
     
388. alls
     
389. sload,1,9,lock,forc,forc_bolt,1,2
     
390. time,10
     
391. alls
     
392. lswrite,1
     
393. 
     
394. time,20
     
395. lswrite,2
     
396. 
     
397. asel,s,,,40
     
398. asel,a,,,64
     
399. asel,a,,,34
     
400. asel,a,,,27
     
401. asel,a,,,51
     
402. asel,a,,,62
     
403. asel,a,,,63
     
404. asel,a,,,36
     
405. asel,a,,,41
     
406. sfa,all,,pres,Pi
     
407. alls
     
408. time,30
     
409. lswrite,3
     
410. 
     
411. !lssolve,1,3,1
     
412. !fini

复制代码

cosmoszzz

其实我接触的参数定义的是有点问题，但是结果符合实验数据。如果常规的在双锥内侧把接触对间隙关闭(开始有一点间隙)，并定义ICON值，结果是有，收敛性很好，但是是不对的，切记切记。

zhuredhigh

附录1在哪里呀？

cosmoszzz

附录1我的权限调的是10，等我问题解决后，下调权限。

huanghe198507

请问有阅读要权限为10的高手吗、帮我下一下附录一？
谢谢

huanghe198507

非常感谢
6# apprent

huanghe198507

为什么楼主的程序在ＡＮＳＹＳ10.0中调试不通呢？
还有楼主为什么确定关闭间隙的方法一定不对？
2# cosmoszzz

yi1xiang

看看我的权限是多少

whj0396

你好，我想问下你能上传下正确的命令流吗，我做这个模拟好久了一直出现错误，看到您的命令流试了下发现也调试不成功，不胜感激

huojunjie

非常感谢

blues-et

:):):):):):)

狂拳之血

apprent 发表于 2009-11-18 23:54
楼主的附录1：

Large-sized Double-cone model_1.rar (2.47 KB)

学习了，谢谢！