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

***********建立边坡模型，并进行网格化分************ Qi&!IG  
;**********第1块*********************** Ih;6(5z  
gen zone wedge p0(325.371509775074,0,376.635571694099) p1 MQ$[jOAqP  
(237.509927879625,0,270.925904145511) p2 W|NzdxCY  
(325.371509775074,10,376.635571694099) p3 ^cZF#%k  
(325.363927118656,0,298.609017598765) p4 U)=StpTT  
(237.509927879625,10,270.925904145511) p5 -n?|,cO  
(325.363927118656,10,298.609017598765 size 5 1 5 Cl7IP<.  
;**********第2块*********************** rqxoqcZ  
gen zone brick p0(325.363927118656,0,298.609017598765) p1 }W#Gf.$6C  
(393.318698206292,0,321.843169023481) p2 @g }r*U?  
(325.363927118656,10,298.609017598765) p3 vd|P
THV_  
(325.371509775074,0,376.635571694099) p4 c+f~>AaI  
(393.318698206292,10,321.843169023481) p5 RgHPYf{  
(325.371509775074,10,376.635571694099) p6 ]!d #2(  
(393.326805700341,0,415.229734921956) p7 Tsxl4ZK  
(393.326805700341,10,415.229734921956) size 5 1 5 Lp5LRw  
;**********第3块*********************** -;a}'1HOE  
gen zone brick p0(393.318698206292,0,321.843169023481) p1 Uq<c+4)5  
(475.498640213354,0,354.365597897843) p2 2o] V q  
(393.318698206292,10,321.843169023481) p3 RI!!?hYm  
(393.326805700341,0,415.229734921956) p4 EFW'D=&h8  
(475.498640213354,10,354.365597897843) p5 M9.jJf  
(393.326805700341,10,415.229734921956) p6 t4GG@`  
(475.50683818756,0,438.723854768785) p7 Y#{ L}  
(475.50683818756,10,438.723854768785) size 5 1 5 $iUK, ?  
;**********第4块*********************** DxdiXf[j  
gen zone brick p0(475.498640213354,0,354.365597897843) p1 b2aPo M=  
(538.970620034247,0,384.218214400487) p2 3!_y@sWx
 
(475.498640213354,10,354.365597897843) p3 X~*1  
(475.50683818756,0,438.723854768785) p4 BKPXXR  
(538.970620034247,10,384.218214400487) p5 b/:9^&z  
(475.50683818756,10,438.723854768785) p6 :;cKns0OA  
(538.978541658805,0,475.463833495854) p7 y#F( xm+L  
(538.978541658805,10,475.463833495854) size 5 1 5 cgNK67"(  
;**********第5块*********************** eGHxiC  
gen zone brick p0(538.970620034247,0,384.218214400487) p1 ^ b{0|:  
(633.394334999949,0,434.810340192916) p2 6%sX<)n%]  
(538.970620034247,10,384.218214400487) p3 Z*tB=  
(538.978541658805,0,475.463833495854) p4 Rb Jl;  
(633.394334999949,10,434.810340192916) p5 oTS*k: C'  
(538.978541658805,10,475.463833495854) p6 Z `sM(?m  
(633.4060080761,0,554.92786334483) p7 ts`c_hH,1'  
(633.4060080761,10,554.92786334483) size 5 1 5 26B]b{Iz{  
;**********第6块*********************** jAB~XaT,  
gen zone brick p0(633.394334999949,0,434.810340192916) p1 dT4?8:  
(702.197699772515,0,474.027380981266) p2 )s5Q4m!  
(633.394334999949,10,434.810340192916) p3 QFY1@2EC  
(633.4060080761,0,554.92786334483) p4 \L>3E#R-Q  
(702.197699772515,10,474.027380981266) p5 5%,J@&5G s  
(633.4060080761,10,554.92786334483) p6 KhIg  
(702.208936770931,0,603.461295685074) p7 }b)7gd=  
(702.208936770931,10,603.461295685074) size 5 1 5 yDg`9q.ckm  
;**********第7块*********************** =e=sK'NvD  
gen zone brick p0(702.197699772515,0,474.027380981266) p1 /oKa?iT  
(764.976557373504,0,510.487720438224) p2 ~OD}`  
(702.197699772515,10,474.027380981266) p3 Azvj(j  
(702.208936770931,0,603.461295685074) p4 =+%QfuK  
(764.976557373504,10,510.487720438224) p5 -`sK?*[{J  
(702.208936770931,10,603.461295685074) p6 cK%Sty'8+  
(764.989858326288,0,663.695394258972) p7
xok8  
(764.989858326288,10,663.695394258972) size 5 1 5 as=Z_a:0N  
;**********第8块*********************** rV*9=  
gen zone brick p0(764.976557373504,0,510.487720438224) p1 _CN5,mLNRk  
(802.260882016456,0,533.452982234561) p2 Q'/v-bd?o  
(764.976557373504,10,510.487720438224) p3 Jz(wXp  
(764.989858326288,0,663.695394258972) p4 `qn
Sq(tNq  
(802.260882016456,10,533.452982234561) p5 lK;|ciq"c7  
(764.989858326288,10,663.695394258972) p6 89e<,f`h  
(802.27850543261,0,714.800304082947) p7 -L%tiz`_  
(802.27850543261,10,714.800304082947) size 5 1 5 :c[iS~ ~Y  
;**********第9块*********************** T}d%XMXq  
gen zone brick p0(802.260882016456,0,533.452982234561) p1 z2*>5c%  
(867.934972023531,0,579.397279648651) p2 aMKi`EW  
(802.260882016456,10,533.452982234561) p3 Opx"'HC@G  
(802.27850543261,0,714.800304082947) p4  G`8i{3:  
(867.934972023531,10,579.397279648651) p5 ! eXDN  
(802.27850543261,10,714.800304082947) p6 O W`yv  
(867.950570236155,0,759.066053053997) p7 g
\^7
Q  
(867.950570236155,10,759.066053053997) size 5 1 5 6..G/,T
B  
;**********第10块*********************** K,(37Id'  
gen zone brick p0(867.934972023531,0,579.397279648651) p1 *;@V5[^3I?  
(938.839302411703,0,641.808924386638) p2 $L%gQkz_  
(867.934972023531,10,579.397279648651) p3 Ox
9WH4E  
(867.950570236155,0,759.066053053997) p4 DdDwMq  
(938.839302411703,10,641.808924386638) p5 bwN>E+  
(867.950570236155,10,759.066053053997) p6 O(q1R#n-}+  
(938.85393110367,0,810.310226793018) p7 S bI7<_  
(938.85393110367,10,810.310226793018) size 5 1 5 m%au* 0p  
;**********第11块*********************** kEiWE|  
gen zone brick p0(938.839302411703,0,641.808924386638) p1 MG6taOO!  
(1006.05516602482,0,712.294023994194) p2 ;%wQ
nhg  
(938.839302411703,10,641.808924386638) p3 'b(V8x  
(938.85393110367,0,810.310226793018) p4 J:WO%P=Q  
(1006.05516602482,10,712.294023994194) p5 {[&$W8Li  
(938.85393110367,10,810.310226793018) p6 xh|NmZg  
(1006.06859579434,0,866.985479665564) p7 Ot,eAiaX  
(1006.06859579434,10,866.985479665564) size 5 1 5 <4UF/G)  
;**********第12块*********************** 1@" L  
gen zone brick p0(1006.05516602482,0,712.294023994194) p1 x[L/d"Wf  
(1081.21471243618,0,798.618156480781) p2 K@d,8[  
(1006.05516602482,10,712.294023994194) p3 Zw6UH;5  
(1006.06859579434,0,866.985479665564) p4 [C_Dv-d  
(1081.21471243618,10,798.618156480781) p5 H|@R+  
(1006.06859579434,10,866.985479665564) p6 bz$)@gLc  
(1081.22524240954,0,906.972998236051) p7 wP:ab  
(1081.22524240954,10,906.972998236051) size 5 1 5 93I.Wp_{  
;**********第13块*********************** e+416 ~X v  
gen zone wedge p0(1223.59495360001,0,974.235305806856) p1 ~Psv[b=]  
(1081.22524240954,0,906.972998236051) p2 ,Hq*zc
c  
(1223.59495360001,10,974.235305806856) p3 f-]5Z
hM'  
(1081.21471243618,0,798.618156480781) p4 !| ObNS  
(1081.22524240954,10,906.972998236051) p5 V% -wZL/  
(1081.21471243618,10,798.618156480781 size 5 1 5 &\Kp_AR  
group 1 A_nu:K-  
;**********第14块*********************** *F[;D7sZ~  
gen zone brick p0(0,0,0) p1(237.483596932835,0,-2.30787651873532E-02) +7%}SV 2)  
p2(0,10,0) p3(0,0,277.520524977887) p4(237.483596932835,10,- R/5@*mv{  
2.30787651873532E-02) p5(0,10,277.520524977887) p6 N-C=O  
(237.509927879625,0,270.925904145511) p7 ;<^t)8E  
(237.509927879625,10,270.925904145511) size 5 1 5 ? }t[  
;**********第15块*********************** &4]~s:F  
gen zone brick p0(237.483597172023,0,-2.06174928067938E-02) p1 s,~)5nL  
(325.334905084496,0,-.031616195723271) p2(237.483597172023,10,- )Pj8{.t4  
2.06174928067938E-02) p3(237.509927879625,0,270.925904145511) p4 lxSCN6  
(325.334905084496,10,-.031616195723271) p5 gxa@da  
(237.509927879625,10,270.925904145511) p6 V'n4iM  
(325.363927118656,0,298.609017598765) p7 ..u2IdEu  
(325.363927118656,10,298.609017598765) size 5 1 5 slPLc  
;**********第16块*********************** 1uz7E  
gen zone brick p0(325.334905412165,0,-2.82444352031419E-02) p1 &v}c3wL]  
(393.29075394712,0,-3.41441235817683E-02) p2(325.334905412165,10,- 4<c#3]  
2.82444352031419E-02) p3(325.363927118657,0,298.609017598765) p4 (>.+tq}  
(393.29075394712,10,-3.41441235817683E-02) p5 Fk
q;Q  
(325.363927118657,10,298.609017598765) p6 `)Z+]5:  
(393.318698206292,0,321.843169023481) p7 Cx,)$!1  
(393.318698206292,10,321.843169023481) size 5 1 5 Pgdv)i3  
;**********第17块*********************** Y\cQ"9  
gen zone brick p0(393.29075394712,0,-3.41441235817634E-02) p1 @-d0~.S  
(475.464198790122,0,-4.12781337960892E-02) p2(393.29075394712,10,- uIeD.I'@{5  
3.41441235817634E-02) p3(393.318698206292,0,321.843169023481) p4 Y#[Wv1hi  
(475.464198790122,10,-4.12781337960892E-02) p5 T1Gy_
G/  
(393.318698206292,10,321.843169023481) p6 B`I9  
(475.498640213354,0,354.365597897843) p7
T~Z7kc'  
(475.498640213354,10,354.365597897843) size 5 1 5 vBRW5@  
;**********第18块*********************** `y%1K|Y=  
gen zone brick p0(475.464198790122,0,-4.12781337960892E-02) p1 i7
*cpNPO  
(538.937259494696,0,-.046788642261033) p2(475.464198790122,10,- &uX|Ksq  
4.12781337960892E-02) p3(475.498640213355,0,354.365597897843) p4 WJXQM[  
(538.937259494696,10,-.046788642261033) p5 HW,2x}[  
(475.498640213355,10,354.365597897843) p6 335\0~;3  
(538.970620034247,0,384.218214400487) p7 Zcq4?-&  
(538.970620034247,10,384.218214400487) size 5 1 5 ZA;VA=)\8  
;**********第19块*********************** XwerQwO=  
gen zone brick p0(538.937259494696,0,-.046788642261033) p1 E]opA$JQ  
(633.352073949682,0,-6.15494458749579E-02) p2(538.937259494696,10, Z{#;my*X|  
-.046788642261033) p3(538.970620034247,0,384.218214400487) p4 ]E8S`[Vn  
(633.352073949682,10,-6.15494458749579E-02) p5 aC
9PlKI  
(538.970620034247,10,384.218214400487) p6 /d/Quro  
(633.394334999947,0,434.810340192915) p7 #"3az8u  
(633.394334999947,10,434.810340192915) size 5 1 5 <;T$?J9  
;**********第20块*********************** @-g'BvS  
gen zone brick p0(633.352074587578,0,-5.49854052973544E-02) p1 hp:8e@  
(702.156541085974,0,-6.09587677106447E-02) p2(633.352074587578,10,- W5#5RK"uX  
5.49854052973544E-02) p3(633.394334999949,0,434.810340192916) p4 %\u
EV  
(702.156541085974,10,-6.09587677106447E-02) p5 \Gzo^w  
(633.394334999949,10,434.810340192916) p6 <kY||  
(702.197699772506,0,474.027380981261) p7 ]t'bd<O  
(702.197699772506,10,474.027380981261) size 5 1 5 Xfq]vQ/{  
;**********第21块*********************** *CnrzrKtQ  
gen zone brick p0(702.156541085974,0,-6.09587677105841E-02) p1 P#8]m(  
(764.932232855049,0,-6.64087330509346E-02) p2(702.156541085974,10,- BUp,bJpO  
6.09587677105841E-02) p3(702.197699772515,0,474.027380981266) p4 slV7,4S&!  
(764.932232855049,10,-6.64087330509346E-02) p5 HJV8P2f8`  
(702.197699772515,10,474.027380981266) p6 74hQ?Atw:  
(764.976557373504,0,510.487720438224) p7 tWo MUp  
(764.976557373504,10,510.487720438224) size 5 1 5 L.]mC!  
;**********第22块*********************** Xy;!Q`h(  
gen zone brick p0(764.932232855049,0,-6.64087330509346E-02) p1 %uVbI'n)  
(802.209034042135,0,-6.96449741618499E-02) p2(764.932232855049,10,- HDTdOG)  
6.64087330509346E-02) p3(764.976557373516,0,510.487720438231) p4 4h[S`;D0Vf  
(802.209034042135,10,-6.96449741618499E-02) p5 9YtdE*,k  
(764.976557373516,10,510.487720438231) p6 0xNlO9b/  
(802.260882016456,0,533.452982234561) p7 p]7Gj&a  
(802.260882016456,10,533.452982234561) size 5 1 5 4Pm
+0=E
 
;**********第23块*********************** x:Nd>Fb  
gen zone brick p0(802.209034042135,0,-6.96449741618499E-02) p1 [;]@PKW?w  
(867.884664242919,0,-7.53467019837899E-02) p2(802.209034042135,10,- $ZQ?E^> B  
6.96449741618499E-02) p3(802.260882016456,0,533.452982234561) p4 YflotlT}  
(867.884664242919,10,-7.53467019837899E-02) p5 %5zztReI  
(802.260882016456,10,533.452982234561) p6 2k^dxk~$V;  
(867.934972023531,0,579.397279648651) p7 k*;2QED  
(867.934972023531,10,579.397279648651) size 5 1 5 @3UVl^T  
;**********第24块*********************** q|8p4X}/]  
gen zone brick p0(867.884664242919,0,-7.53467019837739E-02) p1 /_,} o7@t~  
(938.783575736664,0,-8.15018967641663E-02) p2(867.884664242919,10,- ~6hG"t]:  
7.53467019837739E-02) p3(867.934972023531,0,579.397279648651) p4 w
|1Gb[  
(938.783575736664,10,-8.15018967641663E-02) p5 .W>8bg'u9  
(867.934972023531,10,579.397279648651) p6 9hG+?   
(938.839302411703,0,641.808924386638) p7 fC$~3v  
(938.839302411703,10,641.808924386638) size 5 1 5 }PR^Dj.  
;**********第25块*********************** ,0~/ Cn  
gen zone brick p0(938.783575736664,0,-8.15018967641232E-02) p1 @c8s<9I]  
(1005.99331957413,0,-8.73368109502664E-02) p2(938.783575736664,10,- !Qg%d&q.Sx  
8.15018967641232E-02) p3(938.839302411703,0,641.808924386638) p4 >QusXD"L>  
(1005.99331957413,10,-8.73368109502664E-02) p5 \sy;ca)[6g  
(938.839302411703,10,641.808924386638) p6 B1&H5gxgN  
(1006.05516602482,0,712.294023994194) p7 -CTsB)=\,  
(1006.05516602482,10,712.294023994194) size 5 1 5 ^z[-pTY  
;**********第26块*********************** $LPu_FJ  
gen zone brick p0(1005.99331957413,0,-8.73368109502664E-02) p1 <N{pMz  
(1081.13709323497,0,-9.38605297728576E-02) p2(1005.99331957413,10,- mndUQN_Gb  
8.73368109502664E-02) p3(1006.05516602482,0,712.294023994194) p4 EmDA\9~@R  
(1081.13709323497,10,-9.38605297728576E-02) p5 b=#3p  
(1006.05516602482,10,712.294023994194) p6 FC= %_y  
(1081.21471243616,0,798.618156480762) p7
J6!t"eB+  
(1081.21471243616,10,798.618156480762) size 5 1 5 biFN]D  
;**********第27块*********************** p,?8s%  
gen zone brick p0(1081.13709323497,0,-9.38605297728576E-02) p1 GMw|@?:{  
(1223.50026541007,0,-.118900318576736) p2(1081.13709323497,10,- v!;E1  
9.38605297728576E-02) p3(1081.21471243618,0,798.618156480781) p4 ,]N!I%SI  
(1223.50026541007,10,-.118900318576736) p5 !n`ogzOh  
(1081.21471243618,10,798.618156480781) p6 q5$z:'zE  
(1223.5949536,0,974.235305806856) p7(1223.5949536,10,974.235305806856) YLlw:jN  
size 5 1 5 *5i~N}  
;**********第28块*********************** zr0_SCh;2  
gen zone brick p0(1223.50026664235,0,-.106220001073744) p1 /WDz;,X  
(1369.5591226882,0,-.118900318576736) p2(1223.50026664235,10, h>Z$ n`T  
-.106220001073744) p3(1223.59495360001,0,974.235305806858) p4 |U_48
 
(1369.5591226882,10,-.118900318576736) p5 N_vVEIO9  
(1223.59495360001,10,974.235305806858) p6 u+]v.Mt  
(1369.5591226882,0,1127.69137015435) p7 ^}2 ie|  
(1369.5591226882,10,1127.69137015435) size 5 1 5 dk
V%Pyj  
group 2 range group 1 not 3I_
"vk  
***********定义岩土参数********** J3sO%4sYR  
model elas range group 2 y*M
,&,$  
model mohr range group 2 not z]k=sk  
property bulk 20e9 shear 10e9 density 2400 range group 2 >#@1 I  
property bulk 2e9 shear 1e9 cohesion 12e3 friction 21 density 2000 -(n[^48K  
range group 2 not qj71 rj  
*************施加锚索************* }lgqRg)F9[  
sel cable id=1 begin(325.3175,0.5,376.6356) end(516.1374,0.5,287.6800) b8J\Lm|J  
nseg=10 YbC6&_  
sel cable id=2 begin(460.8365,0.5,434.5298) end(662.0450,0.5,340.7048) jl<rxO?-F  
nseg=10 ^i{B8]2,  
sel cable id=3 begin(572.3805,0.5,503.5765) end(771.5718,0.5,410.6942) 6e@ O88=  
nseg=10 blZiz2F  
sel cable id=4 begin(668.6169,0.5,579.7655) end(874.9548,0.5,483.5486) SnbH`\U"  
nseg=10 ~'WvIA (  
sel cable id=5 begin(759.5143,0.5,658.4420) end(958.6288,0.5,565.5934) 1Tn!.E *  
nseg=10 ?ng?>!  
sel cable id=6 begin(841.5323,0.5,741.2590) end ;G.m;5A
  
(1058.6469,0.5,640.0168) nseg=10 %PG::b  
sel cable id=7 begin(943.9166,0.5,814.5790) end R]%ZqT{PS  
(1159.8808,0.5,713.8732) nseg=10 DPnKr/  
sel cable id=8 begin(1037.8491,0.5,883.8945) end(1242.9164,0.5,788.27) ?6   
nseg=10 !oU$(,#9  
*************定义锚索参数*************
x5m .MQ J  
sel cable property xcarea 9e-3 emod 2e11 ytension 1295e3 ycompression ]lF'o&v]  
-15e4 & %E7+W{?*1  
gr_k= 3e9 gr_per=1.13e-2 gr_coh=600e3 gr_fric=0 slide off )o[ O%b  
sel cable pretension=600e3 G'( %8\  
*************约束模型加自重*************** .2hQ!)+  
plot add axes red Y*``C):K%  
fix y ,{c9Lv%@J  
fix z range z 0 {KL<Hx2M  
fix x range x 0 any x 1369.5591 any =JyYU*G4  
set gravity 0 0 -9.8 qy@gW@IU  
history unbal uNnwz%w  
plot hist 1 1vu4}%nD  
***********求解************** h*hV  
solve 9AJ!7J#v"  
>TsJ0E?3x  
w\;=3C`  
uLb- NxQ-  
  上面是我编的命令流，可是计算结果反映加锚后边坡的塑性区和加锚前是一样的，位移变化也没有任何区别，锚索没有起到应有的作用，不知为何，烦请各位高手帮小弟一把呀，谢谢！ Be2@9  

自己顶下！