在山区地方修建一个高14m的土石坝,土石坝拦住的是14m高的稀泥巴。我想得到的是坝体是否安全，坝体及地基土的整体应力和应变。坝体坡度是按照1：1.75的比例 \wqi_[A
`Xvrf
我建的模型是： U~ QIO O
New +gh6eY8
;title FP =%e]vJ
;Analysis of embankment dam =m6;]16 D
;set logfile Analysisofembankmentdam H 'Q4 IRT
;set log on ~#+ Hhc(
;1 几何建模=================================================
+T02AS
Generate zone brick & >AI65g
p0=(0,0,0) & rAq2
p1=(28,0,0) & 'Y{ fah
p2=(0,1,0) & |GJSAs"L@
p3=(0,0,14) & +]P? ?`,R;
P4=(28,1,0) & @DG $
P5=(0,1,14) & $Kn{x!,"(
P6=(3,0,14) & Q[ J%
P7=(3,1,14) & tb#. Y
size 14 1 7 ratio 1 1 1 ;矩形 xG i,\K\:
plot gp surface n!e4"|4~z
gen zone brick & "HSAwe`5jU
p0 0 0 -10 p1 add 50 0 0 p2 add 0 1 0 p3 add 0 0 10 & -@ f5d
size 25 1 10 ratio 1 1 1 ;矩形 oZTKG'
attach face range z -0.1 0.1 ;联接网格 6jT+kq)
;镜像网格 'N#,,d/G
Generate zone reflect normal (1,0,0) origin=(0,0,0) dq[X:3i
group soildam range z 0 14 [ h20y
group groundsill1 range z -2 0 n 5h4]u
group groundsill2 range z -6 -2 + HvEiY
group groundsill3 range z -10 -6 + >?"P ^
; 2 指定地基土为Mohr Coulomb 模型============================================= *^+8_%;1
model null range z 0 15 ;对坝体设置空单元 U 1!6%x
model moh range z -11 0 ;对地基土网格进行赋值 BMMWP
pro coh 1e10 fric 30 di 2 bulk 6.7e6 she 3.1e6 ten 1e10 range group groundsill1 @ *'$QD,
ini density 1800 range group groundsill1 nw|ls2
pro coh 1e10 fric 10 di 0 bulk 1.11e7 she 3.70e6 ten 1e10 range group groundsill2 b C{1LY0
ini density 1700 range group groundsill2 70IBE[T&
pro coh 1e10 fric 30 di 3 bulk 3.33e8 she 1.11e8 ten 1e10 range group groundsill3 Oo<L~7B
ini density 1700 range group groundsill3 uaN0X"
; 3 边界条件================================================================= Obwj=_+upd
fix x range x 49.9 50.1 dV(61C0wn
fix x range x -50.1 -49.9 Ux5pw
fix y range y -0.1 0.1 9D,`9L5-=
fix y range y 0.9 1.1 *mH++3h
fix x y z range z -10.1 -9.9 R]hilb'a
; 4 初始应力状态========================================================= $2 N)m:X0
set grav 0 0 -9.81 J_) .H d
set large SJ4[n.tPI
ini szz -300e3 grad 0 0 18000 jinDKJ,n;
ini sxx -150e3 grad 0 0 9000 ;s
m )f
ini syy -150e3 grad 0 0 9000 ZdQt!
plot cont szz outline on *V}T}nK7
solve ;第一次求解 /~RY{ c@#L
pro coh 1e4 fric 30 di 2 bulk 1.07e7 she 6.4e6 ten 1e5 range group groundsill1 <2 Q+? L{
pro coh 0.8e4 fric 20 di 0 bulk 5.14e6 she 1.97e6 ten 1e5 range group groundsill2 @}sxA9 a
pro coh 2.5e4 fric 30 di 3 bulk 3.33e8 she 1.11e8 ten 14433.7 range group groundsill3 ;Gf,$dbWn
solve ;第二次求解 }cUq1r-bW
ini xdisp 0 ydisp 0 zdisp 0 ;将节点位移清零 6 2{ (i'K
ini xvel 0 yvel 0 zvel 0 ;将节点速度清零 ;
C/:$l
;plot con xdis kOi@QLdN
save 1地基初始应力场.sav GcXh
V
;5 土石坝施工过程模拟================================================= >p"ytRu^
model elastic range z 0 3 ;1 对坝体激活单元 im Gg3'
pro bulk 1.07e7 she 6.4e6 range z 0 3 ~wtl\-cY
ini density 1800 range z 0 3 ft$@':F
solve @,b tQ_'X
;plot con xdis BWRM
gN'.
model elastic range z 3 6 ;2 对坝体激活单元
vV 5dW
pro bulk 1.07e7 she 6.4e6 range z 3 6 n1D,0+N=
ini density 1800 range z 3 6 j?a^fcXB
solve J${'?!N
model elastic range z 6 9 ;3 对坝体激活单元 3^ NHV g
pro bulk 1.07e7 she 6.4e6 range z 6 9 ^lA Dq']
ini density 1800 range z 6 9 tS|gQUF17
solve ,ru2C_LQ
model elastic range z 9 14 ;4 对坝体激活单元 Z{<& 2*
pro bulk 1.07e7 she 6.4e6 range z 9 14
5cY([4,
ini density 1800 range z 9 14 =z _.R E
solve ,T:Uk*Bj
ini xdisp 0 ydisp 0 zdisp 0 ;将节点位移清零 NrQGoAOw
ini xvel 0 yvel 0 zvel 0 ;将节点速度清零 3.[ fTrzJ
save 2填土石坝.sav =QG@{?JTl
pause &z'N Q !uV
model moh range group soildam ;对坝体激活实际模型 !$}:4}5 6F
pro coh 1e4 fric 30 di 2 bulk 1.07e8 she 6.4e7 ten 1e5 range group soildam ； \v|nRn,`-
solve >/W
save 3填土石坝.sav hsYS<]
pause IY$v%%2WZ
,ddoII
; 2稀泥浆应力状态如何模拟？？？ K ?R*
)_
apply nstress -3.08e5 gradient 0,0,2.2e4 range z -0.5 14.1 x 2.9 28 954!ED|F(
apply szz -3.08e5 range z -2 2 x 28.1 50 wrtJ8O(
step 1 t}$WP&XRG<
Plot surface yellow fap red \53(D7 +
pause B 4fMD]
; 监测变量 D}.Pk>5
hist unbal ;1 ^g^R[8
hist gp zdisp 3 0 14 ;2 Y 8$Y]2
hist gp xdisp 3 0 14 ;3 e+`LtEve0
hist gp zdisp -3 0 14 ;4 T'W)RYnwl
hist gp xdisp -3 0 14 ;5 =E E>QM
hist gp zdisp 28 0 0 ;6 [WZGu6$SU
hist gp xdisp 28 0 0 ;7 Wd9y8z;
hist gp zdisp -28 0 0 ;8 &A}@@d
hist gp xdisp -28 0 0 ;9 gXrXVv<)yw
solve d2cslD d
save 4填土石坝求解完.sav /#WRd}IjK
H > Y0R
1 我这里的内聚力和摩擦角很小的时候，变形特别大？？？请问如何考虑 TU-aL
2稀泥浆应力状态如何模拟？？？ ?%)G %2
4gyC?#Ede
