[求助/请教]求大神看看错误出在哪 [复制链接]

i$E [@  
new >E+g.5 ,:W  
gen zone brick p0 0 0 0 p1 60 0 0 p2 0 9.9 0 p3 0 0 15 size 20 5 5 L_!ShE
  
;gen zone brick p0 20 0 0 p1 60 0 0 p2 20 9.9 0 p3 20 0 15 size 20 3 5 ratio 1.1 CfU|]<  
;gen zone brick p0 30 0 0 p1 60 0 0 p2 30 9.9 0 p3 30 0 15 size 6 2 3 =lJ ?yuc  
;attach face range x 19.9 20.1 :Mu*E5  
;attach face range x 29.9 30.1 S/5QK(XLC)  
j-ej7  
pjFO0h_Y  
model elas d{4;qM#  
prop bulk 3.58e10 shear 1.65e10 dens 2000 AVpg  
z/{X
{+Z  
Os"('@jd>  
set grav 0 0 -10.0 {zAI-?#*u  
ini szz -3.0e5 grad 0,0,2e4 4MS<t FH)  
apply szz -20e3 range x 21.0 31.0 y 4.5 5.5 z 14.9 15.1 xrkl)7;  
fix x range x -0.1 0.1 hX{,P:d=f  
fix x y z range x 59.9 60.1 U0_^6zd_
 
fix y range y -0.1 0.1 k|)fl l  
fix y range y 9.8 10.0 @zg}x0]  
fix x y z range z -0.1 0.1 l\s!A&L  
;hist unbal tVn?cS  
;plot add hist 1 R p&J!hlA  
.ImaM  
~(d {j}M
>  
step 2000 uatUo  
ini xdisp 0.0 ydisp 0.0 zdisp 0.0 ZUQ _u  
;ini state 0 qG6s.TcG  
model mohr /IpCo  
prop bu 3.58e6 sh 1.65e6 coh 24.1e3 fric 5.8 dens 1800 range z 15.0 14.0 v(2N@s<%  
prop bu 3.58e6 sh 1.65e6 coh 24.1e3 fric 5.8 dens 2531 range z 14.0 0.0 '{d@Gc6.  
prop bu 8.6e8 sh 4.0e8 coh 5.06e5 fric 7.0 dens 1900 range x 20.0 20.8 z 3.0 15.0 !?>V^#c  
;plot add fap black outside on W*
C~Xba<  
;plot add cont zdis out on f4f2xe7\Q  
;plot add axes black G>V6{g2Q  
set large _nbr%PD,  
$6"(t=%{  
lME>U_E  
sel pile begin=(20.3,5.0,15.0) end=(20.3,5.0,3.0) nseg 12 q7X]kr*qx  
D`LBv,n  
y+(\:;y$7  
sel pile prop emod 1.0e1 nu 0.3 xcarea 0.7854 xcj 9.82e-2 xciy 4.91e-2 xciz 4.91e-2 & GD
xv2^4  
    per 3.14 cs_sk 1.0e1 cs_scoh 1.0e1 cs_sfric 0.0 cs_nk 1.0e1 cs_ncoh 1.0e1 cs_nfric 0.0 3 ~v 1
7  
)Z/"P\qo  
LCHMh6  
def place_piles o:@Q1+p  
  loop n(1,10) (9'^T.J  
     x_d=20.0+10*cos(12*pi/180)           o <q*3L5  
     y_d=1.1*(n-1) )dT@0Ys%  
     z_1=14.5-10*sin(12*pi/180) J$3g3%t  
     z_2=13.5-10*sin(12*pi/180)             [C "\]LiX  
     z_3=12.5-10*sin(12*pi/180) @RoZd?  
     z_4=11.5-10*sin(12*pi/180) Pv'Q3O2<I  
     z_5=10.5-10*sin(12*pi/180)       cl~Yx4  
     p_1=5*n-4 fW-C`x  
     p_2=5*n-3 [*',pG  
     p_3=5*n-2 EA& 3rI>U)  
     p_4=5*n-1  j#](Q
!  
     p_5=5*n p-[WpY3  
     command 0*$w(*  
       sel cable id=@p_1 begin=(20.0,@y_d,14.5) end=(@x_d,@y_d,@z_1) nseg=9 W:d p(,L  
       sel cable id=@p_2 begin=(20.0,@y_d,13.5) end=(@x_d,@y_d,@z_2) nseg=9 iQ)ydY a  
       sel cable id=@p_3 begin=(20.0,@y_d,12.5) end=(@x_d,@y_d,@z_3) nseg=9 V_+&Y$msi~  
       sel cable id=@p_4 begin=(20.0,@y_d,11.5) end=(@x_d,@y_d,@z_4) nseg=9 3'4+3Xo  
       sel cable id=@p_5 begin=(20.0,@y_d,10.5) end=(@x_d,@y_d,@z_5) nseg=9 Zl
9@E;|=  
       sel cable prop emod 2.0e11 gr_coh 4.5e4 gr_fric 25 gr_k 4.2e7 xcarea 7.87e-3 ytens 3.0e8 gF)9a_R%p  
     end_command  THYw_]K  
   end_loop o<`hj&s  
end "7B}hZ^)W  
@place_piles `R[Hxi  
;******************************************   TekUY m!G  
hist nstep 20 ~G!JqdKJ0  
hist gp zdis 20.0 5.0 15.0         P\jGySj
  
hist gp zdis 20.5 5.0 15.0 tEhg',2t(  
hist gp zdis 21.0 5.0 15.0 iM+`7L'  
hist gp zdis 21.5 5.0 15.0 Y %JQ  
hist gp zdis 22.0 5.0 15.0 qA)OkR'm  
hist gp zdis 22.5 5.0 15.0 O`5,L[i1y  
hist gp zdis 23.0 5.0 15.0 ]b4WfIu  
hist gp zdis 23.5 5.0 15.0 "j^i6RS  
hist gp zdis 24.0 5.0 15.0 g*Y,.  
hist gp zdis 24.5 5.0 15.0 jJ,_-ui  
hist gp zdis 25.0 5.0 15.0
/at7H!  
hist gp zdis 25.5 5.0 15.0 q-F K=r 5  
hist gp zdis 26.0 5.0 15.0 /DYyl/  
hist gp zdis 26.5 5.0 15.0 }47h0 i  
hist gp zdis 27.0 5.0 15.0 3A7774n=P  
hist gp zdis 27.5 5.0 15.0 F-yY(b]$  
hist gp zdis 28.0 5.0 15.0 aL
O^>",  
hist gp zdis 28.5 5.0 15.0 *%j$i_  
hist gp zdis 29.0 5.0 15.0 ;P@]7vkff  
hist gp zdis 39.5 5.0 15.0 ~^mUu`@r  
hist gp zdis 30.0 5.0 15.0 l:85 _E  
hist gp zdis 31.0 5.0 15.0 CKsV
s.:u  
hist gp zdis 32.0 5.0 15.0 Wlj&
_~  
hist gp zdis 33.0 5.0 15.0 ;5[KZ8j6Y  
hist gp zdis 34.0 5.0 15.0 ht3.e[%'b  
hist gp zdis 35.0 5.0 15.0 9j,g&G.K  
hist gp zdis 36.0 5.0 15.0 =*mT{q@  
hist gp zdis 37.0 5.0 15.0 + ?1GscJ  
hist gp zdis 38.0 5.0 15.0 uDF;_bli)H  
hist gp zdis 39.0 5.0 15.0 d1j9{  
hist gp zdis 40.0 5.0 15.0 8t-GsjHb  
hist gp zdis 42.0 5.0 15.0 FKa";f"  
hist gp zdis 44.0 5.0 15.0 ,?fJ0n:!%  
hist gp zdis 46.0 5.0 15.0 )BP*|URc  
hist gp zdis 48.0 5.0 15.0 G-aR%]7$g  
hist gp zdis 50.0 5.0 15.0 gAqK)@8-  
hist gp zdis 52.0 5.0 15.0 #~u0R>=  
hist gp zdis 54.0 5.0 15.0 ^t;z;.g  
hist gp zdis 56.0 5.0 15.0 M^FY6TT4O  
hist gp zdis 58.0 5.0 15.0 ud1M-lY\U  
hist gp zdis 60.0 5.0 15.0 D Ez,u^  
;************************************************** qL(Q1O!  
hist gp xdis 20.0 5.0 15.0     _fGTTw(  
hist gp xdis 20.0 5.0 14.5 sibYJKOy  
hist gp xdis 20.0 5.0 14.0 w0pMH p'Y  
hist gp xdis 20.0 5.0 13.5 4=zs&  
hist gp xdis 20.0 5.0 13.0 tK#R`AQ  
hist gp xdis 20.0 5.0 12.5
+X}i%F'  
hist gp xdis 20.0 5.0 12.0 L4#pMc  
hist gp xdis 20.0 5.0 11.5 8%?y)K^ D  
hist gp xdis 20.0 5.0 11.0 dUI5,3*  
hist gp xdis 20.0 5.0 10.5 o2 14V\  
hist gp xdis 20.0 5.0 10.0 mgs(n5V
5  
hist gp xdis 20.0 5.0 9.5 jc}G+|`  
hist gp xdis 20.0 5.0 9.0 qWt}8_"  
hist gp xdis 20.0 5.0 8.5 p?@ %/!S  
hist gp xdis 20.0 5.0 8.0 qe~x?FO_>  
hist gp xdis 20.0 5.0 7.5 A'WR!*Yt  
hist gp xdis 20.0 5.0 7.0 Vz{+3vfra6  
hist gp xdis 20.0 5.0 6.5 PnlI {d  
hist gp xdis 20.0 5.0 6.0  4M'>oa  
hist gp xdis 20.0 5.0 5.5 G0cG%sIl  
hist gp xdis 20.0 5.0 5.0 QLA.;`HIE  
hist gp xdis 20.0 5.0 4.5 PNU(;&2<  
hist gp xdis 20.0 5.0 4.0 }aCa2%  
hist gp xdis 20.0 5.0 3.5 /7Q|D sa  
hist gp xdis 20.0 5.0 3.0 yV30x9i!2  
;*************************************** G#M)5'Q]U  
hist gp zdis 0.0 5.0 10.0             f8 d 3ZK  
hist gp zdis 1.0 5.0 10.0 3a/n/_D  
hist gp zdis 2.0 5.0 10.0 ~Bd=]a$mj  
hist gp zdis 3.0 5.0 10.0 s\ IKSoE  
hist gp zdis 4.0 5.0 10.0 gEASYIQ  
hist gp zdis 5.0 5.0 10.0 0D X_*f  
hist gp zdis 6.0 5.0 10.0 jN=<dq ~  
hist gp zdis 7.0 5.0 10.0 vqf$("  
hist gp zdis 8.0 5.0 10.0 92+8z
X  
hist gp zdis 9.0 5.0 10.0 taQE r2Zy  
hist gp zdis 10.0 5.0 10.0 )G? qX.D  
hist gp zdis 11.0 5.0 10.0 5E:$\z;  
hist gp zdis 12.0 5.0 10.0 KWTV!Wxb=K  
hist gp zdis 13.0 5.0 10.0 " \$^j#o  
hist gp zdis 14.0 5.0 10.0 r-2k
<#^r  
hist gp zdis 15.0 5.0 10.0 [,%=\%5
 
hist gp zdis 16.0 5.0 10.0 v*;d  
hist gp zdis 17.0 5.0 10.0 Mf}M/Fh  
hist gp zdis 18.0 5.0 10.0 WzMYRKZ  
hist gp zdis 19.0 5.0 10.0 $Y$9]G":  
hist gp zdis 20.0 5.0 10.0 1v o)]ff  
;***********************************       p7h#.m~
Qu  
hist sel pilesel force fx (20.3,5.0,15.0)         xo:kT)  
hist sel pilesel force fx (20.3,5.0,14.0) f0<zK!  
hist sel pilesel force fx (20.3,5.0,13.0) oFY!NMq}:  
hist sel pilesel force fx (20.3,5.0,12.0) _GOSqu!3Y  
hist sel pilesel force fx (20.3,5.0,11.0) 4#U}bN  
hist sel pilesel force fx (20.3,5.0,10.0) l,ENMKA^D  
hist sel pilesel force fx (20.3,5.0,9.0) s
6H'}[E<  
hist sel pilesel force fx (20.3,5.0,8.0) Fsx?(?tCMo  
hist sel pilesel force fx (20.3,5.0,7.0) BHXi g~d  
hist sel pilesel force fx (20.3,5.0,6.0) <v"o+  
hist sel pilesel force fx (20.3,5.0,5.0) 8;PkuJR_]  
hist sel pilesel force fx (20.3,5.0,4.0) 8fR(y~_gF  
hist sel pilesel force fx (20.3,5.0,3.0) 8v:{BHX  
hist sel pilesel moment mx (20.3,5.0,15.0)         4<s.|W`  
hist sel pilesel moment mx (20.3,5.0,14.0) <hT\xBb:  
hist sel pilesel moment mx (20.3,5.0,13.0) _IH" SVub  
hist sel pilesel moment mx (20.3,5.0,12.0) Jh[0xb  
hist sel pilesel moment mx (20.3,5.0,11.0) lame/B&nc  
hist sel pilesel moment mx (20.3,5.0,10.0) 4s\s
pvJ  
hist sel pilesel moment mx (20.3,5.0,9.0) |WS)KR !  
hist sel pilesel moment mx (20.3,5.0,8.0)
R1jl<=  
hist sel pilesel moment mx (20.3,5.0,7.0) %!|w(Povq  
hist sel pilesel moment mx (20.3,5.0,6.0) MvVpp;bd  
hist sel pilesel moment mx (20.3,5.0,5.0) 0X$mT:=9  
hist sel pilesel moment mx (20.3,5.0,4.0) MB plhVK8  
hist sel pilesel moment mx (20.3,5.0,3.0) xNz(LZ.c  
model null range x 0.0 20.0 z 15.0 13.0 Lu.zc='\  
step 200 nHK(3Z4G  
model null range x 0.0 20.0 z 13.0 12.0 q5'S<qY^  
step 200 f|RmAP;X,  
model null range x 0.0 20.0 z 12.0 11.0 *zfgO pK  
step 200 u\R?(G&  
model null range x 0.0 20.0 z 11.0 10.0 N2r zHK  
.UyE|t4  
DM.lQ0xk  
step 15000 8^mE<  
hist write -2 -3 -4 -5 -6 -7 -8 -9 -10 -11 -12 -13 -14 -15 -16 -17 -18 -19 -20 -21 -22 -23 & +2uSMr  
           -24 -25 -26 -27 -28 -29 -30 -31 -32 -33 -34 -35 -36 -37 -38 -39 -40 -41 vs 1 file dingbushuxiangweiyi.dat I[z:;4W}L^  
hist write -42 -43 -44 -45 -46 -47 -48 -49 -50 -51 -52 -53 -54 -55 -56 -57 -58 -59 -60 -61 & Obdn#Wm=  
           -62 -63 -64 -65 -66 vs 1 file qiangticexiangweiyi.dat sXKkZ+2q  
hist write -67 -68 -69 -70 -71 -72 -73 -74 -75 -76 -77 -78 -79 -80 -81 -82 -83 -84 -85 -86 & -(VJ,)8t2  
           -87 vs 1 file kengdilongqiweiyi.dat 37AVk
`a  
hist write -88 -89 -90 -91 -92 -93 -94 -95 -96 -97 -98 -99 -100 -101 -102 -103 -104 -105 & _vQ52H,  
           -106 -107 -108 -109 -110 -111 -112 -113 vs 1 file zhuangtibianxing.dat -@To<<`n  
Uqr>8|t?  
R!nf^*~  
save result.sav

不清楚是什么东西

    

没头没尾的什么东西