水灌内流固耦合问题Ansys实现

TNC

该题实现了一水罐在减速状态下，罐壁的应力分布和流体（水）的压强与速度分布

1. /PREP7
   
2. K,1
   
3. K,2,,1,,
   
4. K,3,,-1,,
   
5. K,4,,,1,
   
6. K,5,,,-1,
   
7. LARC,2,5,1,1,
   
8. LARC,2,4,1,1,
   
9. LARC,4,3,1,1,
   
10. LARC,3,5,1,1,
    
11. FLST,2,4,4
    
12. FITEM,2,2
    
13. FITEM,2,1
    
14. FITEM,2,4
    
15. FITEM,2,3
    
16. AL,P51X
    
17. FLST,2,1,5,ORDE,1
    
18. FITEM,2,1
    
19. ARSCALE,P51X, , ,1,0.5,0.95, ,0,1
    
20. FLST,2,1,5,ORDE,1
    
21. FITEM,2,1
    
22. VEXT,P51X, , ,5.1,0,0,,,,
    
23. SPHERE,1, ,0,360,
    
24. wpro,,,90.000000
    
25. VSBW, 2
    
26. wpro,,,-90.000000
    
27. FLST,3,1,6,ORDE,1
    
28. FITEM,3,4
    
29. VGEN, ,P51X, , , , ,5.1, , ,1
    
30. VPLOT
    
31. VDELE, 4, , ,1
    
32. FLST,2,1,6,ORDE,1
    
33. FITEM,2,3
    
34. VLSCAL,P51X, , ,0.3,0.5,0.95, ,0,1
    
35. VPLOT
    
36. FLST,3,1,6,ORDE,1
    
37. FITEM,3,3
    
38. VSYMM,X,P51X, , , ,0,0
    
39. FLST,3,1,6,ORDE,1
    
40. FITEM,3,3
    
41. VGEN, ,P51X, , , 5.1, , , , ,1
    
42. SAVE
    
43. WPOFF,1.55,,,
    
44. WPRO,,,90.000000
    
45. ASBW,3
    
46. WPRO,,,-90.000000
    
47. WPOFF,2,,,
    
48. wpro,,,90.000000
    
49. ASBW,14
    
50. CSYS,0
    
51. WPAVE,0,0,0
    
52. CSYS,0
    
53. wpro,,,-90.000000
    
54. wpoff,,,0.4
    
55. ASBW,15
    
56. WPOFF,,,0.3
    
57. ASBW,16
    
58. wpoff,,,0.1
    
59. ASBW,15
    
60. CSYS,0
    
61. WPAVE,0,0,0
    
62. CSYS,0
    
63. wpoff,1.55,,
    
64. wpro,,,90.000000
    
65. ASBW,6
    
66. wpoff,,,2
    
67. ASBW,19
    
68. CSYS,0
    
69. WPAVE,0,0,0
    
70. CSYS,0
    
71. wpro,,,-90.000000
    
72. wpoff,,,-0.5
    
73. ASBW,20
    
74. CSYS,0
    
75. WPAVE,0,0,0
    
76. CSYS,0
    
77. wpoff,,,0.4
    
78. ASBW,4
    
79. wpoff,,,-0.8
    
80. ASBW,5
    
81. CSYS,0
    
82. WPAVE,0,0,0
    
83. CSYS,0
    
84. /VIEW, 1 ,1,1,1
    
85. /ANG, 1
    
86. /REP,FAST
    
87. FINISH
    
88. 
    
89. /PREP7
    
90. et,1,142,
    
91. et,2,93
    
92. BTOL,0.10e-9
    
93. vadd,1,2
    
94. vadd,3,4
    
95. BTOL,0.10e-4
    
96. VPLOT
    
97. 
    
98. asel,s,,,1
    
99. asel,a,,,4
    
100. asel,a,,,7,8
     
101. asel,a,,,12,21
     
102. asel,a,,,25,27
     
103. asel,a,,,30
     
104. aatt,2,2,2 ! (material 2)
     
105. alls
     
106. vsel,s,,,1
     
107. vatt,1,1,1 ! (material 1)
     
108. alls
     
109. !
     
110. 
     
111. !划分网格(存在疑问)
     
112. eshape,1
     
113. asel,s,,,1
     
114. asel,a,,,4
     
115. asel,a,,,7,8
     
116. asel,a,,,12,21
     
117. asel,a,,,25,27
     
118. asel,a,,,30
     
119. amesh,all
     
120. CM,_Y,VOLU
     
121. VSEL, , , , 1
     
122. CM,_Y1,VOLU
     
123. CHKMSH,'VOLU'
     
124. CMSEL,S,_Y
     
125. !*
     
126. VMESH,_Y1
     
127. !*
     
128. CMDELE,_Y
     
129. CMDELE,_Y1
     
130. CMDELE,_Y2
     
131. !*
     
132. 
     
133. !设置流体环境
     
134. et,1,142
     
135. et,2,0 ! becomes the Null Element(流固耦合场计算处理方法)
     
136. !! CFD Solution Control
     
137. flda,solu,flow,1 !Solves the momentum and pressure equations
     
138. flda,solu,turb,1 !Activates the turbulence model
     
139. flda,iter,exec,400 !Sets iteration and output controls for steady state analyses
     
140. flda,outp,sumf,10 !set to be defout
     
141. !! CFD Property Information
     
142. flda,prot,dens,constant
     
143. flda,prot,visc,constant
     
144. flda,nomi,dens,1000. ! 1000 kg/m3 for density - water
     
145. flda,nomi,visc,4.6E-4 ! 4.6E-4 kg-s/m (viscosity of water)
     
146. flda,conv,pres,1.E-8 ! Tighten pressure equation convergence
     
147. FLDATA18,METH,PRES,1 ! TDMA
     
148. FLDATA19,TDMA,PRES,200,
     
149. FLDA,RELX,PRES,0.5 !
     
150. FLDA,RELX,Vx,0.5
     
151. FLDATA24,TURB,MODL,5
     
152. !! CFD Boundary Conditions (Applied to Solid Model)
     
153. asel,s,,,7
     
154. asel,a,,,12
     
155. asel,a,,,25
     
156. asel,a,,,27
     
157. da,all,vx,0.,1
     
158. asel,s,,,16,21
     
159. asel,a,,,8
     
160. asel,a,,,14
     
161. asel,a,,,30
     
162. asel,a,,,26
     
163. da,all,vy,0.,1
     
164. da,all,vz,0.,1
     
165. alls
     
166. /title,Fluid Analysis
     
167. physics,write,fluid,fluid
     
168. !
     
169. 
     
170. !设置固体环境
     
171. physics,clear
     
172. et,1,0 ! The Null element for the fluid region
     
173. et,2,93 ! Gasket element - material 2
     
174. R,2,0.001, , , , , ,
     
175. mp,ex,2,3.2947E+9 ! Young's modulus for rubber
     
176. mp,nuxy,2,0.4270 ! Poisson's ratio for the rubber
     
177. asel,s,,,1
     
178. asel,a,,,4
     
179. nsla,,1 !Select all nodes associated with the selected areas
     
180. d,all,ux,0.
     
181. d,all,uy,0.
     
182. d,all,uz,0. ! Fix the bottom of the gasket
     
183. alls
     
184. /title,structural analysis
     
185. finish
     
186. 
     
187. /solu
     
188. antype,static
     
189. nlgeom,on
     
190. cnvtol,f,,,,-1 !Sets convergence(一致性) values for nonlinear analyses;here,no minimum is enforced
     
191. physics,write,struc,struc
     
192. physics,clear
     
193. SAVE
     
194. /solu
     
195. physics,read,fluid
     
196. flda,iter,exec,50
     
197. vsel,s,,,1
     
198. nslv,s,0
     
199. BF,all,FORC,3.68,,, !Defines a nodal body force load
     
200. alls
     
201. solve
     
202. !vsel,s,,,1
     
203. !nslv,s,0
     
204. !BF,all,FORC,3.68,,, !Defines a nodal body force load(另外两个加速度值:5.21,6.25)
     
205. !alls
     
206. !solve
     
207. vsel,s,,,1
     
208. nslv,s,0
     
209. BFDELE,all,forc
     
210. fini
     
211. 
     
212. physics,read,struc ! Read in structures environment
     
213. /assign,esave,struc,esav ! Files for restarting nonlinear structure
     
214. /assign,emat,struc,emat
     
215. /solu
     
216. asel,s,,,1
     
217. asel,a,,,4
     
218. asel,a,,,7,8
     
219. asel,a,,,12,21
     
220. asel,a,,,25,27
     
221. asel,a,,,30
     
222. nsla,,1
     
223. esel,s,type,,2
     
224. ldread,pres,last,,,2,,rfl ! Apply pressure surface load from Flotran
     
225. alls
     
226. rescontrol,,none ! Do not use multiframe restart for nonlinear
     
227. solve
     
228. FINISH

复制代码

FSI

设置非线性迭代收敛残差， 那个CONVERGENCE是收敛的意思，不是一致性。晕

TNC

搞错了，不好意思