请教：基坑桩锚支护模拟过程

888438

请教高手，如何用flac3d模拟桩锚支护的基坑开挖过程，本人编了一些命令流，模拟结果很差，想请高手指点如何模拟这个过程，谢谢。

dreamyhw

期待高手来解惑！！

zouzzz

期待！

月三更

请问楼主在模拟过程中喷射的混凝土是用的什么单元啊？有应力吗？我模拟的总是没有应力出现。急死了……

gaoxin66

new
# t# e; }6 \& ?7 o/ `gen zone brick size 10 10 10 &
- E1 W5 l! {& F5 |2 jp0 0 0 0 p1 10 0 0 p3 0 0 10
' Z" M0 _% X' L& a/ x3 p0 _model mohr
9 Y# K- N- ]3 S0 N0 C  q0 oprop bulk 1e8 shear 0.3e8 fric 35 coh 1e10 ten 1e10
; i) Z. J# p: I' {7 Wset grav 0 0 -9.81$ ~) ^" S& @; U4 V
ini dens 1000  X( {) v1 H' n9 P
fix x range x 0
& o  {: s; f9 Z( }fix x range x 10
: W  Q, R+ s0 mfix y range y 0
: C% z: {7 c6 y, q3 ^- I0 Mfix y range y 10
3 r1 s7 B. a) ]: H. Hfix z range z 0
! o7 k- |- K- whist unbal
  }6 i$ s9 c9 V2 `, ]8 z+ @# X6 wset hist_rep 10: T: S! `! x; B- B
solve/ j) H. n/ i; w0 e
save 初始平衡.sav
* T+ t: N& w, c& c1 y: s
8 v% x. w+ T0 g4 ^, v% [ini xdisp 0 ydisp 0 zdisp 02 M0 }  O6 b8 U0 v
ini xvel 0 yvel 0 zvel 0
( C, D9 u3 v7 d; s; Emodel null range x 0 2 y 0 2 z 0 24 f' M5 ?2 B8 f5 M0 c% d4 Q
solve5 i1 z& ]: K5 ]1 s3 b+ b
save 1次开挖.sav/ G3 @. l0 L, [  _& ~* ?  V
) c% C# |0 _& v
restore 初始平衡.sav
* U# p7 a" u. I6 A+ t0 tini xdisp 0 ydisp 0 zdisp 0( A, O5 T  [$ }2 g0 m
ini xvel 0 yvel 0 zvel 0
8 _* a% i5 {% b# O- g3 [* q7 Fmodel null range x 0 2 y 0 2 z 0 29 ]" s2 x/ `+ m* s" A2 [! E
def place_cables
* q6 X! i9 y2 t, J6 D, C" \    loop m(1,2)  N9 ^9 ]7 I1 }1 Q
    caseof m1 U- b" w: Q1 d8 u6 H7 K8 I- [
    case 1# F0 F) S$ W, E' m& F0 W5 t% U
    x_d1=0.5
! l) U, P; C3 m5 U    x_d2=0.5
* M* G  h( u. i; ]    x_d3=0.5
- f. C/ f; _' `$ V) X8 ?    y_d1=0.5
% q' ^: s. F' K    y_d2=0.5
, U  Z- i: Q6 k) B4 k' f: N    y_d3=0.5: D0 D7 j; p4 y4 K$ i
    z_d1=1.5
0 |. @3 A5 h+ O" }" D4 g    z_d2=2.0
1 k: g' w: V8 R1 ^: F    z_d3=3.51 r" [) t  U/ _( H
       a=1
7 E! ^. `( f  ^9 O0 U' h! w       b=2
/ D6 P& I4 A: z# n% @       c=3
6 E* \& M, q+ ]       d=6" t; {+ ^3 H1 f0 ?: ~
       id_11=1" p$ z# K6 }1 @8 C3 ]2 c* k
       id_111=1! \  L0 ^) k: o
    case 2
4 u3 J6 X0 m6 N% [    x_d1=0.5
  B# v. L2 I7 l; T# S& N    x_d2=0.5
( v/ o' w$ @  `    x_d3=0.5
2 B! d; j) G+ c+ X) V0 Y  c    y_d1=1.5
" L/ V, s- ?& |9 B2 p5 o) G    y_d2=1.5
1 }8 w$ j8 Q+ {    y_d3=1.5
" f, A, a7 g& c) \8 O9 e. x    z_d1=1.5) ^$ ~" ]7 x  m" b5 Z
    z_d2=2.0, x7 H  |! M: T. f7 \  t
    z_d3=3.5' N( q; x4 g6 x- Z) t  V
       a=15$ n! b' }0 F5 W$ B/ d
       b=16
2 g* i* F  @' \& c# F$ |3 m       c=17  I& b# s$ S3 T9 {( I6 h
       d=20
4 i- W; T0 ^# o0 x* e       id_11=3
* [& K0 X2 w' V' F( [       id_111=15* ]* H/ `' r, v; a# d6 x3 b* [
   endcase
% D; K; ?; [2 B& _% j8 \  J   loop n(1,2)
% z0 F& ^3 v, Z6 f. v9 M   x_d1=x_d1+1
9 c" B1 b% Y( L/ a- H   x_d2=x_d2+1" p- _- K$ T2 h
   x_d3=x_d3+1
# T  N. e+ O! l' m$ D) S   y_d1=y_d1% Q" c9 \- Q7 l6 J' l7 q$ K
   y_d2=y_d2! C$ e' M5 ^# `8 T. s8 s
   y_d3=y_d38 F) ~# R4 M1 w
   z_d1=z_d1, G- |6 m7 G5 }6 q+ |
   z_d2=z_d2; n) Q0 r7 q. R. b( x
   z_d3=z_d36 c0 `5 S  i/ E. X- F" J
       a=a+(n-1)*70 ?& ]7 {& M' c, i9 P& C7 i( u: [
       b=b+(n-1)*7
: X' Z' |% M  O7 k6 ~1 B& m. o       c=c+(n-1)*78 j7 t% |8 h8 }1 ]
       d=d+(n-1)*7
! z% A4 m; Z- O" ]        id_11=id_11+1
1 C+ u7 w9 f9 c; Y       id_111=id_111+(n-1)*7
# {" b5 S( d/ B; o( B! c) C  command
/ s$ u/ \. a% F  sel cable id=id_11 beg x_d1 y_d1 z_d1 end x_d2 y_d2 z_d2 nseg 2' u0 Q9 ~- A& A& z7 m. x2 N
  sel cable id=id_11 beg x_d2 y_d2 z_d2 end x_d3 y_d3 z_d3 nseg 4" a3 U6 f- f& p/ p5 F
  sel cable id=id_11 prop emod 2e10 ytension 3.1e5 xcarea 0.0004906 &
8 o! R  P1 y( j6 J7 Y+ N& z   gr_coh 1 gr_k 1 gr_per 0.0785 range cid a,b
" ~; N$ d5 ]  b* S& O8 e5 \  sel cable id= id_11 prop emod 2e10 ytension 3.1e5 xcarea 0.0004906 &
8 V( i# h0 l. }) s7 h: s3 f   gr_coh 10e5 gr_k 2e7 gr_per 0.0785 range cid c,d
2 Q/ [% R" Y6 Z# y) [9 L9 z  sel delete link range id id_1112 e4 E# w4 |3 U% R2 V2 g4 m
  sel link id=id_111 id_111 target zone
. B3 X$ c+ K0 b1 W4 N  sel link attach xdir=rigid ydir=rigid zdir=rigid xrdir=rigid yrdir=rigid zrdir=rigid range id id_111! O. r- v: B+ f. B
  sel cable id=id_11 pretension 60e3 range cid a,d
$ G3 C' k" Z- D9 f6 s7 g  endcommand! h" }  \& e7 X3 I
  end_loop) w  X+ G- W' S- P! L% U
  end_loop
% V" R- [7 `6 x& ~  end1 ]" g( E& |# ?" c# d- o: }) P
  place_cables

ycding

楼主可参此贴：
http://forum.simwe.com/thread-872629-1-1.html