【求助】《ANSYS14.0超级学习手册》三辊卷板机仿真错误

liujj15

手头有这本书，这个算例是17章第一个，我是初学者，跟着提示做下来感觉问题很多。

命令流如下，我自己加了一些注释，大家忙的话可以帮我跑一下看看是不是我自己计算机的问题（做到第一个solve就行了），稍微有点时间的话可以详细看看里面哪里有问题，多谢多谢！
我自己做到 0.332858795s 处就判定不收敛了，再ANSYS14.0和ANSYS15.0都做过了

1. <p><p>/BATCH  
   
2. ! /COM,ANSYS RELEASE 14.0    UP20111024       11:37:55    09/10/2012              
   
3. /input,start140,ans,'D:\Program Files\ANSYS Inc\v140\ANSYS\apdl\',,,,,,,,,,,,,,,,1  
   
4. /PREP7  
   
5. ! /PNUM,VOLU,ON   
   
6. ET,1,SHELL181
   
7. ET,2,SOLID186   
   
8. MP,EX,1,2E11
   
9. MP,DENS,1,7800  
   
10. MP,NUXY,1,0.3   
    
11. MP,EX,2,2E11
    
12. MP,DENS,2,7800  
    
13. MP,NUXY,2,0.3
    
14. ! 材料2双线性随动强化  
    
15. TB,BKIN,2,1
    
16. TBTEMP,0
    
17. TBDATA,,240E6,2E8   
    
18. R,1,0.02 ! 疑问1：这个实数的定义似乎是后面定义shell181的厚度，不过似乎这种方法是没有用的？现在应该是都用section了？
    
19. ! 下面是画出三个轧辊
    
20. CYLIND,0.38/2,0,0.3,2.3,0,360   
    
21. CYLIND,0.38/4,0,0.3,0,0,360
    
22. CYLIND,0.38/4,0,2.6,2.3,0,360   
    
23. VGLUE,ALL   
    
24. VSEL,NONE   
    
25. WPOFF,0.54/2,-0.4   
    
26. CYLIND,0.3/2,0,0.3,2.3,0,360
    
27. CYLIND,0.3/4,0,0.3,0,0,360  
    
28. CYLIND,0.3/4,0,2.6,2.3,0,360
    
29. VGLUE,ALL   
    
30. VSEL,NONE   
    
31. WPOFF,-0.54
    
32. CYLIND,0.3/2,0,0.3,2.3,0,360
    
33. CYLIND,0.3/4,0,0.3,0,0,360  
    
34. CYLIND,0.3/4,0,2.6,2.3,0,360
    
35. VGLUE,ALL
    
36. ! 因为轧辊需要转动，我看到的一个文献上解释说：solid没有办法转动，只能用shell，所以这里删掉了volume只保留了area  
    
37. ALLS
    
38. VDELE,ALL
    
39. ! 定义硬点  
    
40. HPTCREATE,AREA,5,,COORD
    
41. HPTCREATE,AREA,10,,COORD,0,0,2.6
    
42. HPTCREATE,AREA,34,,COORD,-0.54/2,-0.4,2.6   
    
43. HPTCREATE,AREA,29,,COORD,-0.54/2,-0.4,0
    
44. HPTCREATE,AREA,22,,COORD,0.54/2,-0.4,2.6
    
45. HPTCREATE,AREA,17,,COORD,0.54/2,-0.4,0
    
46. ! 定义网格划分参数并划分网格，用了三角网格
    
47. MAT,1   
    
48. TYPE,1  
    
49. REAL,1  
    
50. ESIZE,0.05  
    
51. SMRTSIZE,8  
    
52. MSHAPE,1
    
53. AMESH,ALL
    
54. ! 画出被卷的钢板并进行网格划分  
    
55. BLOCK,0.7,-1.0,0.3/2+.02,0.3/2+.02+0.03,0.5,2.1
    
56. LESIZE,96,,,1   
    
57. MAT,2   
    
58. TYPE,2  
    
59. MSHKEY,1
    
60. MSHAPE,0
    
61. ESIZE,0.05  
    
62. VMESH,ALL
    
63. ! 定义4，5号单元用来接触，疑问2：漏了编号3，不过似乎没有什么影响  
    
64. ET,4,TARGE170   
    
65. ET,5,CONTA174
    
66. ! 定义conta174的12号关键字为rough  
    
67. KEYOPT,5,12,1
    
68. ! 定义三个real，主要为了设置刚度系数为0.1  
    
69. R,10,,,0.1  
    
70. R,11,,,0.1  
    
71. R,12,,,0.1
    
72. ! 选择40号面，应该是钢板的上表面，并进一步选择了这个上表面除开z坐标0.8~1.8范围内的节点
    
73. ! 相当于只是选择了钢板的两侧范围内的节点，下面还要选择轧辊上面两侧的节点，创建接触对
    
74. ASEL,S,,,40
    
75. NSLA,S,1
    
76. NSEL,U,LOC,Z,0.8,1.8
    
77. MAT,2   
    
78. REAL,10
    
79. TYPE,4  
    
80. ESURF   
    
81. ALLS
    
82. ASEL,S,,,3,4,1  
    
83. NSLA,S,1
    
84. NSEL,U,LOC,Z,0.9,1.7
    
85. NSEL,U,LOC,Z,0.3,0.4
    
86. NSEL,U,LOC,Z,2.2,2.3
    
87. MAT,1   
    
88. REAL,10
    
89. TYPE,5  
    
90. ESURF   
    
91. ALLS
    
92. ASEL,S,,,39
    
93. NSLA,S,1
    
94. NSEL,U,LOC,Z,0.8,1.8
    
95. MAT,2   
    
96. REAL,11
    
97. TYPE,4  
    
98. ESURF   
    
99. ALLS
    
100. ASEL,S,,,27,28,1
     
101. NSLA,S,1
     
102. NSEL,U,LOC,Z,0.9,1.7
     
103. NSEL,U,LOC,Z,0.3,0.4
     
104. NSEL,U,LOC,Z,2.2,2.3
     
105. MAT,1   
     
106. REAL,11
     
107. TYPE,5  
     
108. ESURF   
     
109. ALLS
     
110. ASEL,S,,,39
     
111. NSLA,S,1
     
112. NSEL,U,LOC,Z,0.8,1.8
     
113. MAT,2   
     
114. REAL,12
     
115. TYPE,4  
     
116. ESURF   
     
117. ALLS
     
118. ASEL,S,,,15,16,1
     
119. NSLA,S,1
     
120. NSEL,U,LOC,Z,0.9,1.7
     
121. NSEL,U,LOC,Z,0.3,0.4
     
122. NSEL,U,LOC,Z,2.2,2.3
     
123. MAT,1   
     
124. REAL,12
     
125. TYPE,5  
     
126. ESURF   
     
127. ALLS
     
128. FINI
     
129. ! 设置求解载荷及各项参数
     
130. /SOLU   
     
131. D,ALL,UZ
     
132. D,ALL,ROTX  
     
133. D,ALL,ROTY  
     
134. NSEL,S,LOC,X
     
135. NSEL,R,LOC,Y
     
136. D,ALL,UY,-0.08  
     
137. D,ALL,UX
     
138. D,ALL,UZ
     
139. ALLS
     
140. CSYS,4  
     
141. ASEL,S,LOC,Y,-2,0.3/2+0.01  
     
142. NSLA,S,1
     
143. D,ALL,ALL   
     
144. ALLS
     
145. ANTYPE,TRANS
     
146. NEQIT,100   
     
147. LNSRCH,ON   
     
148. NLGEOM,ON   
     
149. TIME,1  
     
150. AUTOT,ON
     
151. NSUBST,30,80,20
     
152. KBC,0   
     
153. OUTRES,ALL,ALL  
     
154. SOLVE   
     
155. !以下模拟下辊转动，驱动钢板前进，但是因为我第一步都做不过，所以以下就根本没有做了
     
156. CSYS,4  
     
157. ASEL,S,LOC,Y,-2,0.3/2+0.01  
     
158. NSLA,S,1
     
159. DDELE,ALL,ALL   
     
160. ALLS
     
161. D,ALL,UZ
     
162. D,ALL,ROTX  
     
163. D,ALL,ROTY  
     
164. NSEL,S,LOC,X,0  
     
165. NSEL,R,LOC,Y
     
166. D,ALL,ROTZ,-3.1415926   
     
167. D,ALL,UX
     
168. D,ALL,UY
     
169. ALLS
     
170. WPOFF,0.54  
     
171. NSEL,S,LOC,X,0  
     
172. NSEL,R,LOC,Y
     
173. D,ALL,UX
     
174. D,ALL,UY
     
175. D,ALL,ROTZ,-3.1415926   
     
176. TIME,2  
     
177. ALLS
     
178. KBC,0   
     
179. NSUBST,200,1000,100
     
180. SOLVE   
     
181. ! /COM,ANSYS RELEASE 15.0    UP20131014       08:33:41    03/04/2016              
     
182. ! /GRA,POWER  
     
183. ! /GST,ON
     
184. ! /PLO,INFO,3
     
185. ! /GRO,CURL,ON
     
186. ! /CPLANE,1   
     
187. ! /REPLOT,RESIZE  
     
188. WPSTYLE,,,,,,,,0
     
189. ! LGWRITE,'lizi-17-1','lgw','C:\Users\CeShi\Desktop\',COMMENT   
     

复制代码

liujj15

有什么其他信息需要的可以在下面留言~我继续添加和解释

liujj15

更新：查出了一个问题，钢板和轧辊没有接触上，留有的空隙导致了刚体位移，把CONTA174的关键字5设为1（自动闭合），情况好了一点，但是运算到一多半仍然出现不收敛的情况