求助边坡加锚 [复制链接]

***********建立边坡模型，并进行网格化分************ $Zn>W@\  
;**********第1块*********************** :Qu.CvYF  
gen zone wedge p0(325.371509775074,0,376.635571694099) p1 m`&6[[)6~  
(237.509927879625,0,270.925904145511) p2 LS/ZZAN u  
(325.371509775074,10,376.635571694099) p3 nmw#4yHYy:  
(325.363927118656,0,298.609017598765) p4 m$`4.>J  
(237.509927879625,10,270.925904145511) p5 S
oHw9FtS  
(325.363927118656,10,298.609017598765 size 5 1 5 wBCBZs$H  
;**********第2块*********************** J3 xi5S  
gen zone brick p0(325.363927118656,0,298.609017598765) p1 ^tL]QE?|  
(393.318698206292,0,321.843169023481) p2 Um!LF"Z  
(325.363927118656,10,298.609017598765) p3 MjW{JR)I  
(325.371509775074,0,376.635571694099) p4 D\Fu4Eg  
(393.318698206292,10,321.843169023481) p5 NQx>u  
(325.371509775074,10,376.635571694099) p6 t vp kc;  
(393.326805700341,0,415.229734921956) p7 eI
cIl2  
(393.326805700341,10,415.229734921956) size 5 1 5 )1/J5DI @8  
;**********第3块*********************** \SooIEl@  
gen zone brick p0(393.318698206292,0,321.843169023481) p1 uoq|l  
(475.498640213354,0,354.365597897843) p2 1}q(Pn2  
(393.318698206292,10,321.843169023481) p3 %Vb~}sT:  
(393.326805700341,0,415.229734921956) p4 Iuve~ugO  
(475.498640213354,10,354.365597897843) p5 x-T7 tr&(  
(393.326805700341,10,415.229734921956) p6 J9&#);(  
(475.50683818756,0,438.723854768785) p7 ]#;u]  
(475.50683818756,10,438.723854768785) size 5 1 5 8Xa{.y"  
;**********第4块*********************** kS62]v]  
gen zone brick p0(475.498640213354,0,354.365597897843) p1 \7WZFh%:  
(538.970620034247,0,384.218214400487) p2 2m,t<Y;  
(475.498640213354,10,354.365597897843) p3 zB8J|uG  
(475.50683818756,0,438.723854768785) p4 uCjbb  
(538.970620034247,10,384.218214400487) p5 .Fx-$Yqy  
(475.50683818756,10,438.723854768785) p6 \Icd>>)*  
(538.978541658805,0,475.463833495854) p7 ~.Er  
(538.978541658805,10,475.463833495854) size 5 1 5 :!w;Y;L:+  
;**********第5块*********************** ~pw%p77)  
gen zone brick p0(538.970620034247,0,384.218214400487) p1 H,(4a2zx  
(633.394334999949,0,434.810340192916) p2 {#N,&?[  
(538.970620034247,10,384.218214400487) p3 ._p^0UxT  
(538.978541658805,0,475.463833495854) p4 x75 3o\u!  
(633.394334999949,10,434.810340192916) p5 9gFfbvd  
(538.978541658805,10,475.463833495854) p6 :M$8<03>F  
(633.4060080761,0,554.92786334483) p7 2,rjy|R`  
(633.4060080761,10,554.92786334483) size 5 1 5 3oC^"723  
;**********第6块*********************** xJ^pqb  
gen zone brick p0(633.394334999949,0,434.810340192916) p1
s+9b.  
(702.197699772515,0,474.027380981266) p2 RY}:&vWDk  
(633.394334999949,10,434.810340192916) p3 0Wb3M"#9<  
(633.4060080761,0,554.92786334483) p4 obK6GG?ZE  
(702.197699772515,10,474.027380981266) p5 YR"IPyj  
(633.4060080761,10,554.92786334483) p6 IAd[_<9D  
(702.208936770931,0,603.461295685074) p7 Of;$ VK'  
(702.208936770931,10,603.461295685074) size 5 1 5 _SrkR7  
;**********第7块*********************** u vyvy  
gen zone brick p0(702.197699772515,0,474.027380981266) p1 +7Qj%x\  
(764.976557373504,0,510.487720438224) p2 uA%cie  
(702.197699772515,10,474.027380981266) p3 u>.>hQ  
(702.208936770931,0,603.461295685074) p4 & &:ZY4`  
(764.976557373504,10,510.487720438224) p5 Ai~d  
(702.208936770931,10,603.461295685074) p6 }
LwKi-G?  
(764.989858326288,0,663.695394258972) p7 e@DVf  
(764.989858326288,10,663.695394258972) size 5 1 5 /Z2 g>  
;**********第8块*********************** \|7Y"WEQ  
gen zone brick p0(764.976557373504,0,510.487720438224) p1 snVeOe#'S  
(802.260882016456,0,533.452982234561) p2 3uuB/8  
(764.976557373504,10,510.487720438224) p3 oz'^.+uvE  
(764.989858326288,0,663.695394258972) p4 6'|NALW  
(802.260882016456,10,533.452982234561) p5 ^.>XDUO F  
(764.989858326288,10,663.695394258972) p6 `L @`l  
(802.27850543261,0,714.800304082947) p7 S[y?>  
(802.27850543261,10,714.800304082947) size 5 1 5 FTh/1"a  
;**********第9块*********************** TUi<  
gen zone brick p0(802.260882016456,0,533.452982234561) p1 /t04}+,e^  
(867.934972023531,0,579.397279648651) p2 `;
Xwv)  
(802.260882016456,10,533.452982234561) p3 o(YF`;OhvS  
(802.27850543261,0,714.800304082947) p4 K 5AArI  
(867.934972023531,10,579.397279648651) p5 Lf+3nN  
(802.27850543261,10,714.800304082947) p6 vPsf{[Kr  
(867.950570236155,0,759.066053053997) p7 'ALe>\WO  
(867.950570236155,10,759.066053053997) size 5 1 5 -:Jn|=  
;**********第10块*********************** r5Xi2!  
gen zone brick p0(867.934972023531,0,579.397279648651) p1 ]m\:XhI*<  
(938.839302411703,0,641.808924386638) p2 nXW]9zC"/  
(867.934972023531,10,579.397279648651) p3 @rlL'|&X*  
(867.950570236155,0,759.066053053997) p4 n==+NL  
(938.839302411703,10,641.808924386638) p5 \GCT3$  
(867.950570236155,10,759.066053053997) p6 gcF V$  
(938.85393110367,0,810.310226793018) p7 72sBx3 ;  
(938.85393110367,10,810.310226793018) size 5 1 5 .~%,eF;l$  
;**********第11块*********************** t+aE*Q  
gen zone brick p0(938.839302411703,0,641.808924386638) p1 *40Z}1ng  
(1006.05516602482,0,712.294023994194) p2 sDS0cc6e  
(938.839302411703,10,641.808924386638) p3 15cgmZsS  
(938.85393110367,0,810.310226793018) p4 sf,9Ym  
(1006.05516602482,10,712.294023994194) p5 d[&Ah~,  
(938.85393110367,10,810.310226793018) p6 pW5PF)([  
(1006.06859579434,0,866.985479665564) p7 kOV6O?h  
(1006.06859579434,10,866.985479665564) size 5 1 5 lT+N{[kLt*  
;**********第12块*********************** ;'oi7b  
gen zone brick p0(1006.05516602482,0,712.294023994194) p1 6AKT-r.  
(1081.21471243618,0,798.618156480781) p2 D n}TO*  
(1006.05516602482,10,712.294023994194) p3 d=a$Gd_$  
(1006.06859579434,0,866.985479665564) p4 q;<Q-jr&O  
(1081.21471243618,10,798.618156480781) p5 +pjU
4>)  
(1006.06859579434,10,866.985479665564) p6 l)qGG$7$  
(1081.22524240954,0,906.972998236051) p7 M(.]?+  
(1081.22524240954,10,906.972998236051) size 5 1 5 jO5Wemqf  
;**********第13块*********************** ;f[@zo><r  
gen zone wedge p0(1223.59495360001,0,974.235305806856) p1 {%8=qJ3@  
(1081.22524240954,0,906.972998236051) p2 H8$";T(I  
(1223.59495360001,10,974.235305806856) p3 9=G dj!L  
(1081.21471243618,0,798.618156480781) p4 |"Fm<  
(1081.22524240954,10,906.972998236051) p5 *cc|(EM  
(1081.21471243618,10,798.618156480781 size 5 1 5 QD^"cPC)mM  
group 1 3&
Fqd  
;**********第14块*********************** 70E@h=oQ  
gen zone brick p0(0,0,0) p1(237.483596932835,0,-2.30787651873532E-02) pJ_>^i=  
p2(0,10,0) p3(0,0,277.520524977887) p4(237.483596932835,10,- W C3b_ia  
2.30787651873532E-02) p5(0,10,277.520524977887) p6 ]Czq A c  
(237.509927879625,0,270.925904145511) p7 sx][X itR+  
(237.509927879625,10,270.925904145511) size 5 1 5 S^  JUQx7  
;**********第15块*********************** ZIJTGa}B q  
gen zone brick p0(237.483597172023,0,-2.06174928067938E-02) p1 +zzS  
(325.334905084496,0,-.031616195723271) p2(237.483597172023,10,- *Edr\P
 
2.06174928067938E-02) p3(237.509927879625,0,270.925904145511) p4 8_uh2`+Bvb  
(325.334905084496,10,-.031616195723271) p5 9S{?@*V  
(237.509927879625,10,270.925904145511) p6 7WN$ rl5/  
(325.363927118656,0,298.609017598765) p7 z1LY|8$G  
(325.363927118656,10,298.609017598765) size 5 1 5 vW03nt86  
;**********第16块*********************** O;RNmiVoq  
gen zone brick p0(325.334905412165,0,-2.82444352031419E-02) p1 .KxE>lJbqM  
(393.29075394712,0,-3.41441235817683E-02) p2(325.334905412165,10,- ;Rd\yAG  
2.82444352031419E-02) p3(325.363927118657,0,298.609017598765) p4 :s
nn-e0l  
(393.29075394712,10,-3.41441235817683E-02) p5 6gD|QC~;  
(325.363927118657,10,298.609017598765) p6 }>m3V2>[  
(393.318698206292,0,321.843169023481) p7 l`vr({A  
(393.318698206292,10,321.843169023481) size 5 1 5 N4wMAT:h  
;**********第17块*********************** "yPKdwP  
gen zone brick p0(393.29075394712,0,-3.41441235817634E-02) p1 C<hb{
$@  
(475.464198790122,0,-4.12781337960892E-02) p2(393.29075394712,10,- du^r EMb%  
3.41441235817634E-02) p3(393.318698206292,0,321.843169023481) p4 \2AXW@xE  
(475.464198790122,10,-4.12781337960892E-02) p5 (Ek=0;Cr  
(393.318698206292,10,321.843169023481) p6 (=w ff5U  
(475.498640213354,0,354.365597897843) p7 @v=A)L  
(475.498640213354,10,354.365597897843) size 5 1 5 ,CjJO -  
;**********第18块*********************** 33w(Pw  
gen zone brick p0(475.464198790122,0,-4.12781337960892E-02) p1 O
p ;){JT  
(538.937259494696,0,-.046788642261033) p2(475.464198790122,10,- vA#?\j2  
4.12781337960892E-02) p3(475.498640213355,0,354.365597897843) p4 F>rf cW2  
(538.937259494696,10,-.046788642261033) p5 Kvh6D"  
(475.498640213355,10,354.365597897843) p6 5U.,iQ(d  
(538.970620034247,0,384.218214400487) p7 YL@d+ -\  
(538.970620034247,10,384.218214400487) size 5 1 5 )q'~<QxI\  
;**********第19块*********************** >}B~~C;  
gen zone brick p0(538.937259494696,0,-.046788642261033) p1 LJBDB6  
(633.352073949682,0,-6.15494458749579E-02) p2(538.937259494696,10, z<s4-GJ)?  
-.046788642261033) p3(538.970620034247,0,384.218214400487) p4 q^+Z>
  
(633.352073949682,10,-6.15494458749579E-02) p5 $:<G=  
(538.970620034247,10,384.218214400487) p6 bn8?-  
(633.394334999947,0,434.810340192915) p7 cP[]\r+Kj  
(633.394334999947,10,434.810340192915) size 5 1 5 g#<M/qn  
;**********第20块*********************** et :v4^*f  
gen zone brick p0(633.352074587578,0,-5.49854052973544E-02) p1 dWhF[q"  
(702.156541085974,0,-6.09587677106447E-02) p2(633.352074587578,10,- -O>*` O>M  
5.49854052973544E-02) p3(633.394334999949,0,434.810340192916) p4 P8;f^3V(+/  
(702.156541085974,10,-6.09587677106447E-02) p5 2O)2#N  
(633.394334999949,10,434.810340192916) p6 ot.R Gpg%  
(702.197699772506,0,474.027380981261) p7 'W|@d8}h  
(702.197699772506,10,474.027380981261) size 5 1 5 :]-? l4(%  
;**********第21块*********************** -I{J]L$S#  
gen zone brick p0(702.156541085974,0,-6.09587677105841E-02) p1 AV?<D.<  
(764.932232855049,0,-6.64087330509346E-02) p2(702.156541085974,10,- U4,hEnJBT  
6.09587677105841E-02) p3(702.197699772515,0,474.027380981266) p4 0%5x&vx'S  
(764.932232855049,10,-6.64087330509346E-02) p5 ?mW;%d~]  
(702.197699772515,10,474.027380981266) p6 jY5BVTWnV  
(764.976557373504,0,510.487720438224) p7 -cnlj  
(764.976557373504,10,510.487720438224) size 5 1 5 K7$Vl"l  
;**********第22块*********************** g bwg3$!9  
gen zone brick p0(764.932232855049,0,-6.64087330509346E-02) p1 T(6B,  
(802.209034042135,0,-6.96449741618499E-02) p2(764.932232855049,10,- bHH=MLZR:  
6.64087330509346E-02) p3(764.976557373516,0,510.487720438231) p4 8Zvh"Z?  
(802.209034042135,10,-6.96449741618499E-02) p5 .@;,'Xw1~  
(764.976557373516,10,510.487720438231) p6 f>C|qDmT  
(802.260882016456,0,533.452982234561) p7 >jBnNA@  
(802.260882016456,10,533.452982234561) size 5 1 5 Y$6W~j  
;**********第23块*********************** o!M*cyq  
gen zone brick p0(802.209034042135,0,-6.96449741618499E-02) p1 O7\)C]A  
(867.884664242919,0,-7.53467019837899E-02) p2(802.209034042135,10,- AZ
adNuL/  
6.96449741618499E-02) p3(802.260882016456,0,533.452982234561) p4
Z|a\rNv  
(867.884664242919,10,-7.53467019837899E-02) p5 (r )fx  
(802.260882016456,10,533.452982234561) p6 parC~)b_  
(867.934972023531,0,579.397279648651) p7 -~ycr[}x  
(867.934972023531,10,579.397279648651) size 5 1 5 9{5 c}bX  
;**********第24块*********************** g63?(+Fz  
gen zone brick p0(867.884664242919,0,-7.53467019837739E-02) p1 8v']>5S]#  
(938.783575736664,0,-8.15018967641663E-02) p2(867.884664242919,10,- {>=#7e-]  
7.53467019837739E-02) p3(867.934972023531,0,579.397279648651) p4 m7~[f7U  
(938.783575736664,10,-8.15018967641663E-02) p5 c}g:vh  
(867.934972023531,10,579.397279648651) p6 1w|V'e?kb  
(938.839302411703,0,641.808924386638) p7 M)!:o/!cS  
(938.839302411703,10,641.808924386638) size 5 1 5 &)|3OJ'o  
;**********第25块*********************** s\i.pd:Q  
gen zone brick p0(938.783575736664,0,-8.15018967641232E-02) p1 G/C5o=cY  
(1005.99331957413,0,-8.73368109502664E-02) p2(938.783575736664,10,- .
G ~,h  
8.15018967641232E-02) p3(938.839302411703,0,641.808924386638) p4 [E
V}P&U  
(1005.99331957413,10,-8.73368109502664E-02) p5 O"x/O#66  
(938.839302411703,10,641.808924386638) p6 N0
G-/  
(1006.05516602482,0,712.294023994194) p7 |A@Gch fd  
(1006.05516602482,10,712.294023994194) size 5 1 5 Eb=;D1)y]  
;**********第26块*********************** =v]
eQIp  
gen zone brick p0(1005.99331957413,0,-8.73368109502664E-02) p1  \l8$1p  
(1081.13709323497,0,-9.38605297728576E-02) p2(1005.99331957413,10,- YP vg(T  
8.73368109502664E-02) p3(1006.05516602482,0,712.294023994194) p4 d<l
-Ldle  
(1081.13709323497,10,-9.38605297728576E-02) p5 Y&_1U/}h  
(1006.05516602482,10,712.294023994194) p6 {cBLm/C  
(1081.21471243616,0,798.618156480762) p7 9=Rj9%  
(1081.21471243616,10,798.618156480762) size 5 1 5 G.c@
4Wz+  
;**********第27块*********************** N^8 lfc$a  
gen zone brick p0(1081.13709323497,0,-9.38605297728576E-02) p1 UT7"
.1H  
(1223.50026541007,0,-.118900318576736) p2(1081.13709323497,10,- AAt<{  
9.38605297728576E-02) p3(1081.21471243618,0,798.618156480781) p4 8yc?9&/|  
(1223.50026541007,10,-.118900318576736) p5 ld*RL:G  
(1081.21471243618,10,798.618156480781) p6 zVs|go>F  
(1223.5949536,0,974.235305806856) p7(1223.5949536,10,974.235305806856) me`(J y<  
size 5 1 5 aXefi'!6  
;**********第28块*********************** $[P>nRhW  
gen zone brick p0(1223.50026664235,0,-.106220001073744) p1 WS)u{ or  
(1369.5591226882,0,-.118900318576736) p2(1223.50026664235,10, JTg0T+  
-.106220001073744) p3(1223.59495360001,0,974.235305806858) p4 O@bDMg  
(1369.5591226882,10,-.118900318576736) p5 6'N!)b^-  
(1223.59495360001,10,974.235305806858) p6 CmPix]YMQ  
(1369.5591226882,0,1127.69137015435) p7 )04lf*ti  
(1369.5591226882,10,1127.69137015435) size 5 1 5 ICgyCsZ,  
group 2 range group 1 not ';?b99  
***********定义岩土参数********** @7*Ag~MRb  
model elas range group 2 /A) v$Bv=  
model mohr range group 2 not er0ClvB  
property bulk 20e9 shear 10e9 density 2400 range group 2 |d1%N'Ll  
property bulk 2e9 shear 1e9 cohesion 12e3 friction 21 density 2000 n"{oj7E0a  
range group 2 not ?OPAf4h  
*************施加锚索************* {bvm83{T  
sel cable id=1 begin(325.3175,0.5,376.6356) end(516.1374,0.5,287.6800) e/h7x\Z  
nseg=10 $W;IW$  
sel cable id=2 begin(460.8365,0.5,434.5298) end(662.0450,0.5,340.7048) 8M;G@ Q80  
nseg=10 id.W"5+  
sel cable id=3 begin(572.3805,0.5,503.5765) end(771.5718,0.5,410.6942) |_;Vb  
nseg=10 J8yi#A>+  
sel cable id=4 begin(668.6169,0.5,579.7655) end(874.9548,0.5,483.5486) D;Jb'Be  
nseg=10 !j.jvI%e;  
sel cable id=5 begin(759.5143,0.5,658.4420) end(958.6288,0.5,565.5934) Zm@ O[:~  
nseg=10 ;.r >  
sel cable id=6 begin(841.5323,0.5,741.2590) end WtTwY8HC  
(1058.6469,0.5,640.0168) nseg=10 #Rdq^TGMi;  
sel cable id=7 begin(943.9166,0.5,814.5790) end P'6(HT>F?  
(1159.8808,0.5,713.8732) nseg=10 nNff~u)I  
sel cable id=8 begin(1037.8491,0.5,883.8945) end(1242.9164,0.5,788.27) !S',V&Yb  
nseg=10 K*Tvo`  
*************定义锚索参数************* l7r!fAV-f  
sel cable property xcarea 9e-3 emod 2e11 ytension 1295e3 ycompression (FAd'$lhX}  
-15e4 & IK-E{,iK
c  
gr_k= 3e9 gr_per=1.13e-2 gr_coh=600e3 gr_fric=0 slide off 6\9 9WQ  
sel cable pretension=600e3 `-N&cc  
*************约束模型加自重*************** uz+b  
plot add axes red ?$^qcpJCp  
fix y p }bTI5  
fix z range z 0 hrRX=  
fix x range x 0 any x 1369.5591 any fE/8;v!
=  
set gravity 0 0 -9.8 =}~hbPJM  
history unbal -j_J1P0,  
plot hist 1 kM?p>V6  
***********求解************** 8}W06k>)%  
solve y]
`@%V2P  
cuP5cL/Y  
l<$c.GgFd  
B1A5b=6G<  
  上面是我编的命令流，可是计算结果反映加锚后边坡的塑性区和加锚前是一样的，位移变化也没有任何区别，锚索没有起到应有的作用，不知为何，烦请各位高手帮小弟一把呀，谢谢！ V ;)q?ZHg  

自己顶下！