声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 5732|回复: 34

[前后处理] 大家看看这个网格是如何划分出来的

[复制链接]
发表于 2006-5-31 13:53 | 显示全部楼层 |阅读模式

马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。

您需要 登录 才可以下载或查看,没有账号?我要加入

x
第一个图是一本博士学位论文上有限元划分的图
第二个图是我划分的(计算不收敛)
想请教下大家
第一个图的网格是如何划分出来的?
我认为的难点就是圆柱不能Mapped划分,要划分出这种六面体只能用sweep
但是圆柱和长方体不能一起sweep,只能各自sweep
但是各自sweep的效果,明显不好,见图三
所以只有在接触处用free 然后进行单元退化
然而生成的四面体较多计算时不收敛
如何能划分出如图1这样规则而且耦合的又好的网格出来呢?
希望大家畅所欲言 ~~
小弟不胜感激
1.jpg
2.jpg
3.jpg
回复
分享到:

使用道具 举报

发表于 2006-5-31 17:02 | 显示全部楼层
试着把上面的小圆柱分割下来
规则的体就可以用mapped了
 楼主| 发表于 2006-5-31 17:16 | 显示全部楼层
但是做起来真的很难啊 ~~
怎么分割呢?
 楼主| 发表于 2006-6-1 11:08 | 显示全部楼层
长方体里面还有钢筋
不能divide
哎 ~~
用四面体单元吧
不收敛
快疯了
 楼主| 发表于 2006-6-1 16:35 | 显示全部楼层
  1. /PREP7
  2. et,1,link8 !钢筋单元
  3. et,2,solid65 !承台体单元
  4. et,3,solid95 !过渡单元
  5. et,4,solid45 !加载区分散荷载用
  6. r,1,113 !12号钢筋
  7. r,2
  8. KEYOPT,2,1,1
  9. KEYOPT,2,3,0
  10. KEYOPT,2,5,0
  11. KEYOPT,2,6,0
  12. KEYOPT,2,7,1
  13. KEYOPT,2,8,0
  14. mp,ex,1,2.1e5 !定义钢筋属性
  15. mp,prxy,1,0.3
  16. tb,bkin,1
  17. tbdata,,335
  18. mp,ex,2,3e4 !定义混凝土承台体属性
  19. mp,prxy,2,0.2
  20. *SET,fc,20.97
  21. *SET,ft,2.2
  22. TB,CONC,2,
  23. TBTEMP,0
  24. TBDATA,,0.5,0.9,ft,-1,,
  25. TBDATA,,,,,,,
  26. tb,miso,2
  27. tbtemp,0
  28. tbpt,,0,0
  29. tbpt,,0.0002,0.0002*3e4
  30. tbpt,,0.0004,0.36*fc
  31. tbpt,,0.0006,0.51*fc
  32. tbpt,,0.0008,0.64*fc
  33. tbpt,,0.001,0.75*fc
  34. tbpt,,0.0012,0.84*fc
  35. tbpt,,0.0014,0.91*fc
  36. tbpt,,0.0016,0.96*fc
  37. tbpt,,0.0018,0.99*fc
  38. tbpt,,0.002,fc
  39. tbpt,,0.0025,0.95*fc
  40. tbpt,,0.003,0.9*fc
  41. tbpt,,0.0035,0.85*fc
  42. tbpt,,0.004,0.8*fc
  43. mp,ex,3,3.65e4 ! 定义桩、柱混凝土属性
  44. mp,prxy,3,0.2
  45. mp,ex,4,3.3e4 !定义加载区单元属性
  46. mp,prxy,4,0.2

  47. n,1,50,50,50 !建钢筋单元
  48. NGEN,15,1,1 , , ,50, , ,1,
  49. NGEN,9,15,1, 15, 1, ,50, ,1,
  50. type,1
  51. mat,1
  52. r,1
  53. *do,ii,1,14,1
  54. e,ii,ii+1
  55. *enddo
  56. FLST,4,14,2,ORDE,2
  57. FITEM,4,1
  58. FITEM,4,-14
  59. EGEN,9,15,P51X, , , , , , , , ,50, ,
  60. *do,ii,1,106,15
  61. e,ii,ii+15
  62. *enddo
  63. FLST,5,8,2,ORDE,2
  64. FITEM,5,127
  65. FITEM,5,-134
  66. ENGEN,8,15,1,P51X, , , , , , , ,50, , ,

  67. blc4,,,800,500,400 !建混凝土承台体模型
  68. wpoff,,,50
  69. vsbw,all
  70. vglue,all
  71. wpcsys,-1,0
  72. wpoff,300,150,400
  73. blc4,,,200,200,150 !建柱体
  74. wpcsys,-1,0
  75. wpoff,100,100
  76. CYL4, , ,50, , , ,-80 !建桩体
  77. VGEN,2,4, , , ,300, , ,0
  78. VGEN,3,4, 5, ,300, , , ,0
  79. allsel
  80. vsel,s,loc,z,500,550 !加载区划分单元
  81. aslv,a
  82. lsla,s
  83. lesize,all,50
  84. vatt,4,2,4
  85. vmesh,all
  86. allsel
  87. vsel,s,loc,z,400,500 !柱体划分单元
  88. aslv,a
  89. lsla,s
  90. lesize,all,50
  91. vatt,3,2,2
  92. vmesh,all
  93. allsel
  94. vsel,s,loc,z,50,400 !承台体划分单元
  95. aslv,a
  96. lsla,s
  97. lesize,all,50
  98. vatt,2,2,2
  99. vmesh,all
  100. allsel
  101. vsel,s,loc,z,-80,50 !划分桩及其桩承台接触层
  102. SMRT,7
  103. mshape,1,3d
  104. MSHKEY,0
  105. vatt,3,2,3
  106. vmesh,all
  107. allsel,all
  108. tchg,95,92,0
  109. allsel

  110. asel,s,loc,z,-80 !桩底施加约束
  111. nsla,s
  112. d,all,all
  113. FINISH
  114. /solu
  115. ANTYPE,0 !非线性设置
  116. CNVTOL,F, ,0.05,2,0.01,
  117. outres,all,all
  118. autots,on
  119. lnsrch,on
  120. arclen,off
  121. ncnv,2
  122. neqit,20
  123. pred,on

  124. allsel !施加荷载
  125. nsubst,20
  126. allsel
  127. nsel,s,loc,z,550
  128. f,all,fz,-16000
  129. kbc,0
  130. nsel,all
  131. lswrite,1
  132. allsel
  133. fcum,add,1
  134. nsubst,20
  135. allsel
  136. nsel,s,loc,z,550
  137. f,all,fz,-16000
  138. kbc,0
  139. nsel,all
  140. lswrite,2
  141. allsel
  142. fcum,add,1
  143. nsubst,20
  144. allsel
  145. nsel,s,loc,z,550
  146. f,all,fz,-8000
  147. kbc,0
  148. nsel,all
  149. lswrite,3
  150. allsel
  151. fcum,add,1
  152. nsubst,20
  153. allsel
  154. nsel,s,loc,z,550
  155. f,all,fz,-8000
  156. kbc,0
  157. nsel,all
  158. lswrite,4
  159. allsel
  160. fcum,add,1
  161. nsubst,20
  162. allsel
  163. nsel,s,loc,z,550
  164. f,all,fz,-4000
  165. kbc,0
  166. nsel,all
  167. lswrite,5
  168. save
复制代码
发表于 2006-6-1 17:20 | 显示全部楼层

看看这个网格是如何划分出来

这个是用abq做的,
picture1:
1.png
发表于 2006-6-1 17:22 | 显示全部楼层

回复:(linqus)看看这个网格是如何划分出来

picture2:
一半模型
1.png
发表于 2006-6-1 17:38 | 显示全部楼层

回复:(linqus)回复:(linqus)看看这个网格是如何...

picture3
底部:

共32976个单元
1.png
 楼主| 发表于 2006-6-1 19:30 | 显示全部楼层
看起来比较不错哦
不过好多单元数
尺寸取的比我的要小很多 ~~
发表于 2006-6-1 20:28 | 显示全部楼层
这个网格划的真不错,希望作者能把划分的方法奉献一下,让大家也学习一下。谢谢!
 楼主| 发表于 2006-6-2 10:15 | 显示全部楼层
。。。。
他是用abq做的
发表于 2006-6-3 09:16 | 显示全部楼层

用adina自由划分做的。

尝试着用adina做一个,自由划分,
效果还行:

pic1:9792elements
1.jpg
发表于 2006-6-3 09:17 | 显示全部楼层

用adina自由划分做的。

pic2:
22464elements
1.jpg
发表于 2006-6-3 09:20 | 显示全部楼层

用adina自由划分做的。

由于没有加以控制,
所以划分时间较abq长一些,
还有,底部网格与abq划分不同。
手头没有ansys,没做。呵呵

底部:
1.jpg
 楼主| 发表于 2006-6-4 09:47 | 显示全部楼层

回复:(渊源1983)大家看看这个网格是如何划分出来的...

强人!!
不过这两种软件都不熟悉
也没时间去熟悉了
就想用ansys做下 ~~
您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

QQ|小黑屋|Archiver|手机版|联系我们|声振论坛

GMT+8, 2024-11-11 19:06 , Processed in 0.078853 second(s), 21 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

快速回复 返回顶部 返回列表