桩土接触的命令流，运行后不收敛，请指点

car2008

这是一段桩土接触的命令流，运行后不收敛，不知怎样改，请高手指点以下，非常感谢
桩承受竖向荷载,桩土间建立了接触面

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. 
    
23. LSEL,S,,,1,3,2
    
24. LESIZE,ALL,,,2
    
25. LSEL,S,,,2,4,2
    
26. LESIZE,ALL,,,16
    
27. AMESH,1
    
28. LSEL,S,,,6,8,2
    
29. LESIZE,ALL,0.5
    
30. LSEL,S,,,5
    
31. LESIZE,ALL,0.5
    
32. LSEL,S,,,7
    
33. LESIZE,ALL,0.5
    
34. AMESH,2
    
35. LSEL,S,,,11
    
36. LESIZE,ALL,0.5
    
37. LSEL,S,,,12
    
38. LESIZE,ALL,0.5
    
39. LSEL,S,,,10
    
40. LESIZE,ALL,0.5
    
41. LSEL,S,,,9
    
42. LESIZE,ALL,0.5
    
43. AMESH,3
    
44. 
    
45. EXTOPT,ESIZE,2,0,
    
46. EXTOPT,ACLEAR,1
    
47. TYPE,2
    
48. MAT,1
    
49. VEXT,1,,,,,-1
    
50. EXTOPT,ESIZE,2,0,
    
51. EXTOPT,ACLEAR,1
    
52. MAT,2
    
53. VEXT,2,3,1,,,-1
    
54. ALLSEL
    
55. /VIEW,1,1,1,1
    
56. EPLOT
    
57. 
    
58. K,1001,5,6,-5
    
59. K,1002,5,0,-5
    
60. K,1003,5,-16,-5
    
61. L,18,1001
    
62. L,1001,1002
    
63. L,1002,17
    
64. L,1002,1003
    
65. L,1003,21
    
66. AL,22,36,37,38
    
67. AL,30,38,39,40
    
68. VEXT,19,20,1,-5
    
69. VPLOT
    
70. 
    
71. VSEL,S,,,4,5,1
    
72. VPLOT
    
73. LSEL,S,,,36,52,1
    
74. LESIZE,ALL,0.5
    
75. TYPE,2
    
76. MAT,2
    
77. VMESH,ALL
    
78. ALLSEL
    
79. VSEL,S,,,2,5,1
    
80. VPLOT
    
81. ESLV,S
    
82. NSLV,S,1
    
83. EPLOT
    
84. NUMMRG,ALL
    
85. NUMCMP,ALL
    
86. ALLSEL
    
87. /PNUM,AREA,1
    
88. /VIEW,1,1,1,1
    
89. /ANG,1
    
90. /PNUM,MAT,1,2
    
91. EPLOT
    
92. allsel
    
93. ET,3,TARGE170   
    
94. ET,4,CONTA173   
    
95. KEYOPT,4,9,0  
    
96. KEYOPT,4,12,2
    
97. r,10
    
98. real,10
    
99. mp,mu,2,0.2
    
100. vsel,s,,,1
     
101. vplot
     
102. asel,s,,,4,6,1     
     
103. aplot
     
104. cm,target,area  
     
105. type,3  
     
106. nsla,s,1
     
107. esurf,all   
     
108. esel,s,type,,3  
     
109. eplot   
     
110. /psymb,esys,1   
     
111. allsel
     
112. 
     
113. vsel,s,,,2,5,1  
     
114. vplot
     
115. asel,s,,,13
     
116. asel,a,,,17
     
117. asel,a,,,22
     
118. aplot
     
119. cm,contact,area
     
120. type,4  
     
121. nsla,s,1
     
122. esurf,all   
     
123. esel,s,type,,4  
     
124. eplot   
     
125. /psymb,esys,1   
     
126. allsel  
     
127. gplot   
     
128. SAVE
     
129. finish      
     
130. 
     
131. /solu
     
132. asel,s,loc,x,0
     
133. aplot
     
134. da,all,symm
     
135. asel,s,loc,x,5
     
136. da,all,all
     
137. aplot
     
138. asel,s,loc,z,0
     
139. da,all,symm
     
140. asel,s,loc,z,-5
     
141. da,all,all
     
142. asel,s,loc,y,-16
     
143. da,all,all
     
144. allsel
     
145. acel,0,9.8,0,
     
146. nsel,s,loc,y,8
     
147. esln,s
     
148. eplot
     
149. sfe,all,3,pres,,500
     
150. antype,static
     
151. nlgeom,on
     
152. time,1
     
153. nsub,100
     
154. outres,all,all
     
155. allsel
     
156. solv
     
157. finish
     
158. 总是出现不收敛问题，请高手指点
     

复制代码

[ 本帖最后由 16443 于 2008-9-9 08:54 编辑 ]

car2008

今天得到启发，可能是D-P材料参数设置出问题，正在调整中，

car2008

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

原因 出在单位 要统一  
见 Ansys 帮助文件
内聚力与弹性模量的单位 要一致