为什么沉降不对称啊我做的模型仿真币全给你

laoyin2

FINISH  ! Make sure we are at BEGIN level   
/CLEAR,NOSTART  ! Clear model since no SAVE found   
/NOPR   
/PMETH,OFF,0
KEYW,PR_SET,1   
KEYW,PR_STRUC,1



/PREP7  
!*********************面
ET,1,plane42


!**********************混凝土
ET,2,SOLID45

!********************8围岩
ET,3,SOLID65




!******************************混凝土mp1
!18.96mpa(兆帕)
MP,EX,1,20000e06
MP,PRXY,1,0.2  


!**************************************岩土mp2


MP,EX,2,300e6
MP,NUXY,2,0.4
mp,dens,2,2e3



!*******************************点  
K, ,0,0,,   
K, ,-4,-4,,
K, ,-4,4,,  
K, ,4,-4,,  
K, ,4,4,,   


k,,2,2,
k,,2,-2
k,,-2,-2
k,,-2,2


!************************线***************************************


LSTR,       2,       4  
LSTR,       4,       5  
LSTR,       5,       3  
LSTR,       3,       2


LSTR,       3,       9  
LSTR,       6,       5  
LSTR,       7,       4  
LSTR,       8,       2  

larc,6,7,1,2.828
larc,7,8,1,2.828
larc,8,9,1,2.828
larc,9,6,1,2.828


al,9,10,11,12


al,3,5,6,12
al,4,5,8,11
al,1,7,8,10
al,2,6,7,9

lsel,s,line,,1,12,1

lesize,all,,,4


asel,s,area,,2,5,1
mshkey,3
amesh,all

asel,s,area,,1
mshkey,3
amesh,all




!******************************************拉伸**********************************
mat,1                    
! 混凝土
type,2
!
esize,,5
! Anzahl der Elemente in L鋘gsrichtung
asel,s,area,,1           
! 选面从几到几
vext,all,,,,,20         
! 拉伸50m
allsel


!拉伸
mat,2                    
! 岩土
type,3
!砌衬
esize,,5
! Anzahl der Elemente in L鋘gsrichtung
asel,s,area,,2,5,1
! 选面从几到几
vext,all,,,,,20         
! 拉伸50m
ALLSEL,ALL  


!**********************清除面
asel,s,type,,1
aclear,all

!重力荷载
acel,,10 !加速度


!边界条件

nsel,s,loc,y,-4
! Markiert alle Punkte in y-Richtung auf H鰄e null
D,All,UY,0
! Freiheitsgrad UY festgehalten ...

allsel
nsel,s,loc,x,4
d,all,ux,0

allsel
nsel,s,loc,x,-4
d,all,ux,0

allsel
nsel,s,loc,z,0
d,all,uz,0

allsel
nsel,s,loc,z,20
d,all,uz,0



!last********************************************************last
/solu
!under face!可以不考虑  以后加荷载







!solution*********************************


/SOL  
!nsel,all
alls
eplo
ALLSEL



/solu                !进入求解器
antype,static            !指定为静力学求解器
pred,on                !打开时间步长预测器
lnsrch,on            !打开线性搜索
nlgeom,on            !打开大变形求解
!pcg 法
!EQSLV,PCG, ,1,
!nropt,full            !设定全N-R求解
OUTRES,ALL,all,      !最后一步，all，所有项
!NSUBST,subs,100,10
!Substepps
NEQIT,1000   
!迭代次数
!CNVTOL,F, ,0.01,2, ,   !精度是0.01

laoyin2

模型很简单 就是2种材料 一个圆形 一个方形 在重力下变形为什么沉降不一样呢 报酬50元 帮我看出来

laoyin2

我的目的是想让他们变形不一样 可是有个面两种材料黏在一起了呢

wuchg999

两种模型材料弹摸相差太大了，你单独显示每一种材料就对了

nufeng

1. FINISH  ! Make sure we are at BEGIN level   
   
2. /CLEAR,NOSTART  ! Clear model since no SAVE found   
   
3. /COM,ANSYS RELEASE 11.0    UP20070125       14:41:44    06/24/2010
   
4. /NOPR   
   
5. /PMETH,OFF,0
   
6. KEYW,PR_SET,1   
   
7. KEYW,PR_STRUC,1
   
8. /PREP7  
   
9. !*********************  
   
10. ET,1,plane42
    
11. !**********************
    
12. ET,2,SOLID45
    
13. !********************8  
    
14. ET,3,SOLID65
    
15. !******************************混凝土mp1
    
16. !18.96mpa(兆帕)   
    
17. MP,EX,1,20000e06
    
18. MP,PRXY,1,0.2   
    
19. !**************************************岩土mp2
    
20. MP,EX,2,300e6   
    
21. MP,NUXY,2,0.4   
    
22. mp,dens,2,2e3   
    
23. !*******************************
    
24. K, ,0,0,,   
    
25. K, ,-4,-4,,
    
26. K, ,-4,4,,  
    
27. K, ,4,-4,,  
    
28. K, ,4,4,,   
    
29. k,,2,2,
    
30. k,,2,-2
    
31. k,,-2,-2
    
32. k,,-2,2
    
33. !************************线***************************************
    
34. LSTR,       2,       4  
    
35. LSTR,       4,       5  
    
36. LSTR,       5,       3  
    
37. LSTR,       3,       2  
    
38. LSTR,       3,       9  
    
39. LSTR,       6,       5  
    
40. LSTR,       7,       4  
    
41. LSTR,       8,       2  
    
42. larc,6,7,1,2.828
    
43. larc,7,8,1,2.828
    
44. larc,8,9,1,2.828
    
45. larc,9,6,1,2.828
    
46. al,9,10,11,12   
    
47. al,3,5,6,12
    
48. al,4,5,8,11
    
49. al,1,7,8,10
    
50. al,2,6,7,9  
    
51. lsel,s,line,,1,12,1
    
52. 
    
53. FLST,2,5,5,ORDE,2   
    
54. FITEM,2,1   
    
55. FITEM,2,-5  
    
56. VEXT,all, , ,0,0,20,,,,
    
57. !************************mesh
    
58. vatt,2,,3
    
59. vsel,s,,,1
    
60. vatt,1,,2
    
61. alls
    
62. vmesh,all
    
63. acel,,10 !  
    
64. !   
    
65. nsel,s,loc,y,-4
    
66. ! Markiert alle Punkte in y-Richtung auf H鰄e null
    
67. D,All,UY,0  
    
68. ! Freiheitsgrad UY festgehalten ...
    
69. allsel  
    
70. nsel,s,loc,x,4  
    
71. d,all,ux,0  
    
72. allsel  
    
73. nsel,s,loc,x,-4
    
74. d,all,ux,0  
    
75. allsel  
    
76. nsel,s,loc,z,0  
    
77. d,all,uz,0  
    
78. allsel  
    
79. nsel,s,loc,z,20
    
80. d,all,uz,0  
    
81. !last********************************************************last   
    
82. FINISH  
    
83. /SOL
    
84. !under face!
    
85. !solution*********************************  
    
86. /SOL
    
87. !nsel,all   
    
88. alls
    
89. eplo
    
90. ALLSEL  
    
91. /solu                !  
    
92. antype,static            !  
    
93. pred,on                !
    
94. lnsrch,on            !  
    
95. nlgeom,on            !  
    
96. !pcg
    
97. !EQSLV,PCG, ,1,
    
98. !nropt,full            !设定全N-R
    
99. OUTRES,ALL,all,      !最后一步，all
    
100. !NSUBST,subs,100,10
     
101. !Substepps  
     
102. NEQIT,1000  
     
103. !   
     
104. !CNVTOL,F, ,0.01,2, ,   !精度是0.01  
     
105. /STATUS,SOLU
     
106. SOLVE   
     
107. FINISH  
     
108. /POST1  
     
109. plnsol,u,y  

复制代码

laoyin2

不对你那样这两个材料就连在一起了

nufeng

两个材料不连， 要用接触分析吧！？

nufeng

求解之前定义解除对，否则不收敛

1. 
   
2. /prep7
   
3. !**********************混凝土
   
4. ET,1,SOLID45
   
5. !********************8围岩
   
6. ET,2,SOLID65
   
7. 
   
8. 
   
9. !******************************混凝土mp1
   
10. !18.96mpa(兆帕)
    
11. MP,EX,1,20000e06
    
12. MP,PRXY,1,0.2  
    
13. mp,dens,1,2e3
    
14. !**************************************岩土mp2
    
15. 
    
16. MP,EX,2,300e6
    
17. MP,NUXY,2,0.4
    
18. mp,dens,2,2e3
    
19. BLOCK,-4,4,-4,4,0,20,   
    
20. CYL4, , ,2*sqrt(2), , , ,20
    
21. VOVLAP,all
    
22. VDELE,2, , ,1
    
23. CYL4, , ,2*sqrt(2), , , ,20
    
24. LESIZE,all,1
    
25. VSWEEP,all
    
26. vsel,s,,,3
    
27. eslv,s
    
28. EMODIF,all,MAT,2,  
    
29. EMODIF,all,TYPE,2,
    
30. 
    
31. /sol
    
32. allsel
    
33. !重力荷载
    
34. acel,,10 !加速度
    
35. 
    
36. !边界条件
    
37. nsel,s,loc,y,-4
    
38. ! Markiert alle Punkte in y-Richtung auf H鰄e null
    
39. D,All,UY,0
    
40. ! Freiheitsgrad UY festgehalten ...
    
41. allsel
    
42. nsel,s,loc,x,4
    
43. d,all,ux,0
    
44. allsel
    
45. nsel,s,loc,x,-4
    
46. d,all,ux,0
    
47. allsel
    
48. nsel,s,loc,z,0
    
49. d,all,uz,0
    
50. allsel
    
51. nsel,s,loc,z,20
    
52. d,all,uz,0
    
53. 
    
54. !last********************************************************last
    
55. /solu
    
56. !under face!可以不考虑  以后加荷载
    
57. !solution*********************************
    
58. 
    
59. /SOL  
    
60. !nsel,all
    
61. alls
    
62. eplo
    
63. ALLSEL
    
64. 
    
65. /solu                !进入求解器
    
66. antype,static            !指定为静力学求解器
    
67. pred,on                !打开时间步长预测器
    
68. lnsrch,on            !打开线性搜索
    
69. nlgeom,on            !打开大变形求解
    
70. !pcg 法
    
71. !EQSLV,PCG, ,1,
    
72. !nropt,full            !设定全N-R求解
    
73. OUTRES,ALL,all,      !最后一步，all，所有项
    
74. !NSUBST,subs,100,10
    
75. !Substepps
    
76. NEQIT,1000   
    
77. !迭代次数
    
78. !CNVTOL,F, ,0.01,2, ,   !精度是0.01
    

复制代码

ETsy

建立楼主把复合材料的建模研究一下