在山区地方修建一个高14m的土石坝,土石坝拦住的是14m高的稀泥巴。我想得到的是坝体是否安全，坝体及地基土的整体应力和应变。坝体坡度是按照1：1.75的比例 TD*AFR3Oz
3? <A]"X .
我建的模型是： }6pr.-J
New qc.TY p
;title 5(\/ b<#
;Analysis of embankment dam 'AWWdz
;set logfile Analysisofembankmentdam ]b >XN8y.
;set log on 8dLmsk ^
;1 几何建模================================================= !gV{[j?~zr
Generate zone brick & :-U& _%#w
p0=(0,0,0) & @:B}QxC
p1=(28,0,0) & Y@q9
p2=(0,1,0) & Im-qGB0C
p3=(0,0,14) & Z_dL@\#|
P4=(28,1,0) & ~`Vo0Z*S
P5=(0,1,14) & yMM2us#*+q
P6=(3,0,14) & >EA\KrjW
P7=(3,1,14) & b?cO+PY01
size 14 1 7 ratio 1 1 1 ;矩形 G9xO>Xp^Al
plot gp surface I(kEvfxc"
gen zone brick & +A1*e+/b\
p0 0 0 -10 p1 add 50 0 0 p2 add 0 1 0 p3 add 0 0 10 & y xT} hMa
size 25 1 10 ratio 1 1 1 ;矩形 RTbV!I
attach face range z -0.1 0.1 ;联接网格 >dgq2ok!u
;镜像网格 ar
7.O;e
Generate zone reflect normal (1,0,0) origin=(0,0,0) _qk&W_u
group soildam range z 0 14 \(=xc2
group groundsill1 range z -2 0 [-"ZuU G
group groundsill2 range z -6 -2 :6%ivS
group groundsill3 range z -10 -6 8[t*VIXI
; 2 指定地基土为Mohr Coulomb 模型============================================= hT_Q_1,
model null range z 0 15 ;对坝体设置空单元 nO'C2)bBSG
model moh range z -11 0 ;对地基土网格进行赋值 *' es(]W
pro coh 1e10 fric 30 di 2 bulk 6.7e6 she 3.1e6 ten 1e10 range group groundsill1 ;XyryCo
ini density 1800 range group groundsill1 DzA'M X
pro coh 1e10 fric 10 di 0 bulk 1.11e7 she 3.70e6 ten 1e10 range group groundsill2 htrtiJ1
ini density 1700 range group groundsill2 i"Hc( lg
pro coh 1e10 fric 30 di 3 bulk 3.33e8 she 1.11e8 ten 1e10 range group groundsill3 3G 5xIr6
ini density 1700 range group groundsill3 (Rr C<5"
; 3 边界条件================================================================= o(> #}[N}
fix x range x 49.9 50.1 Z
eY *5m
fix x range x -50.1 -49.9 Ktt(l-e +
fix y range y -0.1 0.1 )+Z.J]$O-
fix y range y 0.9 1.1 #H
|p)2k
fix x y z range z -10.1 -9.9 z19%!k
; 4 初始应力状态========================================================= )h>dD
set grav 0 0 -9.81 ]oz >/\!
set large ^jb;4nf
ini szz -300e3 grad 0 0 18000 xL<c/B`-:
ini sxx -150e3 grad 0 0 9000 ^?\|2H
ini syy -150e3 grad 0 0 9000 9An \uH)mL
plot cont szz outline on U6wy^!_X9
solve ;第一次求解 UUbO\_&y
pro coh 1e4 fric 30 di 2 bulk 1.07e7 she 6.4e6 ten 1e5 range group groundsill1 D%?9[Qb
pro coh 0.8e4 fric 20 di 0 bulk 5.14e6 she 1.97e6 ten 1e5 range group groundsill2 nQ|($V1?W
pro coh 2.5e4 fric 30 di 3 bulk 3.33e8 she 1.11e8 ten 14433.7 range group groundsill3 Y `$\o
solve ;第二次求解 LfU? 1:Du
ini xdisp 0 ydisp 0 zdisp 0 ;将节点位移清零 xe(7q1
ini xvel 0 yvel 0 zvel 0 ;将节点速度清零 I` j G
;plot con xdis l O*
save 1地基初始应力场.sav tQxxm= >
;5 土石坝施工过程模拟================================================= l_9Z zN
model elastic range z 0 3 ;1 对坝体激活单元 &Qj1uf92.
pro bulk 1.07e7 she 6.4e6 range z 0 3 9C Ki$L
ini density 1800 range z 0 3 ,JbP~2M~%
solve m :~y:.
;plot con xdis c68y\
model elastic range z 3 6 ;2 对坝体激活单元 5 A 5t
pro bulk 1.07e7 she 6.4e6 range z 3 6 @e\
@EW
ini density 1800 range z 3 6 "Y5 :{Kj
solve cD!E.2[
model elastic range z 6 9 ;3 对坝体激活单元 c05-1
pro bulk 1.07e7 she 6.4e6 range z 6 9 u0)9IZxc
ini density 1800 range z 6 9 vr?u= _%Z
solve ./.aL Th
model elastic range z 9 14 ;4 对坝体激活单元 P|lDW|}D@
pro bulk 1.07e7 she 6.4e6 range z 9 14 G; pmR^
ini density 1800 range z 9 14 n) D
solve =;Co0Q`
ini xdisp 0 ydisp 0 zdisp 0 ;将节点位移清零 XhWo~zh"
ini xvel 0 yvel 0 zvel 0 ;将节点速度清零 lk81IhI
save 2填土石坝.sav y0?HZ Xq
pause qe
e_wx
model moh range group soildam ;对坝体激活实际模型 cH :&S=>h
pro coh 1e4 fric 30 di 2 bulk 1.07e8 she 6.4e7 ten 1e5 range group soildam ； r|
\""
solve y] O&w{m$
save 3填土石坝.sav O}2/w2n
pause e0ni
eLgq
)
; 2稀泥浆应力状态如何模拟？？？ XDyo=A]
apply nstress -3.08e5 gradient 0,0,2.2e4 range z -0.5 14.1 x 2.9 28 v_v>gPl,
apply szz -3.08e5 range z -2 2 x 28.1 50 &
@_PY
step 1 X&rsWk
Plot surface yellow fap red <4@8T7
pause N'l2$8
; 监测变量 7)2 Q
hist unbal ;1 9B83H V4 J
hist gp zdisp 3 0 14 ;2 (Jj xrZ+L
hist gp xdisp 3 0 14 ;3 #D?w,<_8,
hist gp zdisp -3 0 14 ;4 #uC}IX2n
hist gp xdisp -3 0 14 ;5 %z-s o?gF
hist gp zdisp 28 0 0 ;6 -byaV;T?"
hist gp xdisp 28 0 0 ;7 n ;vZY
hist gp zdisp -28 0 0 ;8 >o& %via}
hist gp xdisp -28 0 0 ;9 6CGk*s
solve 
