声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 1254|回复: 2

[结构分析] 桩土接触的命令流,运行后不收敛,请指点

[复制链接]
发表于 2008-9-8 15:16 | 显示全部楼层 |阅读模式

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

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

x
这是一段桩土接触的命令流,运行后不收敛,不知怎样改,请高手指点以下,非常感谢
桩承受竖向荷载,桩土间建立了接触面
  1. /PREP7   
  2. ET,1,PLANE42
  3. ET,2,SOLID45
  4. MP,EX,1,2.5e10
  5. MP,NUXY,1,0.2  
  6. MP,DENS,1,2500
  7. MP,EX,2,2.5E8  
  8. MP,NUXY,2,0.4  
  9. MP,DENS,2,2000
  10. TB,DP,2,,,  
  11. TBDATA,1,19,32,30  
  12. RECTNG,0,1,0,8,
  13. RECTNG,1,5,0,6,
  14. RECTNG,0,5,0,-16,   
  15. /PNUM,AREA,1
  16. /PNUM,LINE,1
  17. ASEL,S,,,2,3,1  
  18. AGLUE,ALL   
  19. NUMCMP,ALL  
  20. ALLSEL
  21. APLOT

  22. LSEL,S,,,1,3,2
  23. LESIZE,ALL,,,2
  24. LSEL,S,,,2,4,2
  25. LESIZE,ALL,,,16
  26. AMESH,1
  27. LSEL,S,,,6,8,2
  28. LESIZE,ALL,0.5
  29. LSEL,S,,,5
  30. LESIZE,ALL,0.5
  31. LSEL,S,,,7
  32. LESIZE,ALL,0.5
  33. AMESH,2
  34. LSEL,S,,,11
  35. LESIZE,ALL,0.5
  36. LSEL,S,,,12
  37. LESIZE,ALL,0.5
  38. LSEL,S,,,10
  39. LESIZE,ALL,0.5
  40. LSEL,S,,,9
  41. LESIZE,ALL,0.5
  42. AMESH,3

  43. EXTOPT,ESIZE,2,0,
  44. EXTOPT,ACLEAR,1
  45. TYPE,2
  46. MAT,1
  47. VEXT,1,,,,,-1
  48. EXTOPT,ESIZE,2,0,
  49. EXTOPT,ACLEAR,1
  50. MAT,2
  51. VEXT,2,3,1,,,-1
  52. ALLSEL
  53. /VIEW,1,1,1,1
  54. EPLOT

  55. K,1001,5,6,-5
  56. K,1002,5,0,-5
  57. K,1003,5,-16,-5
  58. L,18,1001
  59. L,1001,1002
  60. L,1002,17
  61. L,1002,1003
  62. L,1003,21
  63. AL,22,36,37,38
  64. AL,30,38,39,40
  65. VEXT,19,20,1,-5
  66. VPLOT

  67. VSEL,S,,,4,5,1
  68. VPLOT
  69. LSEL,S,,,36,52,1
  70. LESIZE,ALL,0.5
  71. TYPE,2
  72. MAT,2
  73. VMESH,ALL
  74. ALLSEL
  75. VSEL,S,,,2,5,1
  76. VPLOT
  77. ESLV,S
  78. NSLV,S,1
  79. EPLOT
  80. NUMMRG,ALL
  81. NUMCMP,ALL
  82. ALLSEL
  83. /PNUM,AREA,1
  84. /VIEW,1,1,1,1
  85. /ANG,1
  86. /PNUM,MAT,1,2
  87. EPLOT
  88. allsel
  89. ET,3,TARGE170   
  90. ET,4,CONTA173   
  91. KEYOPT,4,9,0  
  92. KEYOPT,4,12,2
  93. r,10
  94. real,10
  95. mp,mu,2,0.2
  96. vsel,s,,,1
  97. vplot
  98. asel,s,,,4,6,1     
  99. aplot
  100. cm,target,area  
  101. type,3  
  102. nsla,s,1
  103. esurf,all   
  104. esel,s,type,,3  
  105. eplot   
  106. /psymb,esys,1   
  107. allsel

  108. vsel,s,,,2,5,1  
  109. vplot
  110. asel,s,,,13
  111. asel,a,,,17
  112. asel,a,,,22
  113. aplot
  114. cm,contact,area
  115. type,4  
  116. nsla,s,1
  117. esurf,all   
  118. esel,s,type,,4  
  119. eplot   
  120. /psymb,esys,1   
  121. allsel  
  122. gplot   
  123. SAVE
  124. finish      

  125. /solu
  126. asel,s,loc,x,0
  127. aplot
  128. da,all,symm
  129. asel,s,loc,x,5
  130. da,all,all
  131. aplot
  132. asel,s,loc,z,0
  133. da,all,symm
  134. asel,s,loc,z,-5
  135. da,all,all
  136. asel,s,loc,y,-16
  137. da,all,all
  138. allsel
  139. acel,0,9.8,0,
  140. nsel,s,loc,y,8
  141. esln,s
  142. eplot
  143. sfe,all,3,pres,,500
  144. antype,static
  145. nlgeom,on
  146. time,1
  147. nsub,100
  148. outres,all,all
  149. allsel
  150. solv
  151. finish
  152. 总是出现不收敛问题,请高手指点
复制代码

[ 本帖最后由 16443 于 2008-9-9 08:54 编辑 ]
回复
分享到:

使用道具 举报

 楼主| 发表于 2008-9-23 14:59 | 显示全部楼层
今天得到启发,可能是D-P材料参数设置出问题,正在调整中,
 楼主| 发表于 2008-9-26 19:33 | 显示全部楼层

收敛了


TB,DP,2,,,  
TBDATA,1,19,32,30
改为
TB,DP,2,,,  
TBDATA,1,19000,32,30  
即可收敛

原因 出在单位 要统一  
见 Ansys 帮助文件
内聚力与弹性模量的单位 要一致
您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

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

GMT+8, 2024-11-14 02:00 , Processed in 0.075826 second(s), 17 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

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