在山区地方修建一个高14m的土石坝,土石坝拦住的是14m高的稀泥巴。我想得到的是坝体是否安全，坝体及地基土的整体应力和应变。坝体坡度是按照1：1.75的比例 ixJUq o
-_jV.`t
我建的模型是： V:42\b7x
New >l$vu-k)~4
;title $XS0:C0
;Analysis of embankment dam ~L(_q]
;set logfile Analysisofembankmentdam 0w >DU^+
;set log on c ;3bX6RD*
;1 几何建模================================================= $ ,k SR}
Generate zone brick & -X@;"0v
p0=(0,0,0) & O$
i6r]j_
p1=(28,0,0) & oeXNb4; 4
p2=(0,1,0) & nf< <]iHf
p3=(0,0,14) & 5iola}6
P4=(28,1,0) & CiP-Zh[gZ
P5=(0,1,14) & < %Qw
dEO
P6=(3,0,14) & SwQ.tK1p
P7=(3,1,14) & S(A0),
size 14 1 7 ratio 1 1 1 ;矩形 X.0/F6U
plot gp surface d9/E^)TT
gen zone brick & dE5DH~ldV
p0 0 0 -10 p1 add 50 0 0 p2 add 0 1 0 p3 add 0 0 10 & 1A23 G$D
size 25 1 10 ratio 1 1 1 ;矩形 ;{|a~e?Y
attach face range z -0.1 0.1 ;联接网格 V mQ7M4j*
;镜像网格 H+zn:j@ ~L
Generate zone reflect normal (1,0,0) origin=(0,0,0) #SY 8Zv
group soildam range z 0 14 \Rn.ug
group groundsill1 range z -2 0 $"{V],:T
|
group groundsill2 range z -6 -2 AK<ZP ?0
group groundsill3 range z -10 -6 ADX}
; 2 指定地基土为Mohr Coulomb 模型============================================= v)JQb-<
model null range z 0 15 ;对坝体设置空单元 XA])<dZ
model moh range z -11 0 ;对地基土网格进行赋值 \h^bOxh
pro coh 1e10 fric 30 di 2 bulk 6.7e6 she 3.1e6 ten 1e10 range group groundsill1 3 &*0n^g
ini density 1800 range group groundsill1 aMTY{
pro coh 1e10 fric 10 di 0 bulk 1.11e7 she 3.70e6 ten 1e10 range group groundsill2 ,8u u,,c
ini density 1700 range group groundsill2 OU=9fw
pro coh 1e10 fric 30 di 3 bulk 3.33e8 she 1.11e8 ten 1e10 range group groundsill3 ;U<) $5
ini density 1700 range group groundsill3 $52Te3n
; 3 边界条件=================================================================
Y6A]dk
fix x range x 49.9 50.1 RCt)qh+
fix x range x -50.1 -49.9 Ja- D}|;
fix y range y -0.1 0.1 .Jdw:
fix y range y 0.9 1.1 DT&[W<oN
fix x y z range z -10.1 -9.9 ?Di, '
; 4 初始应力状态========================================================= [Hdk=p
set grav 0 0 -9.81 Fga9
set large K.
G#[
ini szz -300e3 grad 0 0 18000 @{_PO {=\C
ini sxx -150e3 grad 0 0 9000 NrvS/ cI!t
ini syy -150e3 grad 0 0 9000 Pyi PhOJe
plot cont szz outline on '4sT+q
solve ;第一次求解 \3q{E",\>@
pro coh 1e4 fric 30 di 2 bulk 1.07e7 she 6.4e6 ten 1e5 range group groundsill1 ilK8V4k<T)
pro coh 0.8e4 fric 20 di 0 bulk 5.14e6 she 1.97e6 ten 1e5 range group groundsill2 m@JU).NKCS
pro coh 2.5e4 fric 30 di 3 bulk 3.33e8 she 1.11e8 ten 14433.7 range group groundsill3 |PN-,f{ -
solve ;第二次求解 ;hs:wLVa"
ini xdisp 0 ydisp 0 zdisp 0 ;将节点位移清零 |xz qYu?o
ini xvel 0 yvel 0 zvel 0 ;将节点速度清零 6\86E$f=h
;plot con xdis _}] o~
save 1地基初始应力场.sav 'OGOT0(
;5 土石坝施工过程模拟================================================= 4 \(;}M-R{
model elastic range z 0 3 ;1 对坝体激活单元 I#l9
pro bulk 1.07e7 she 6.4e6 range z 0 3 dh/:H/k kR
ini density 1800 range z 0 3 %9mCgHQ9
solve (Cp :NS
;plot con xdis :0T]p"y4
model elastic range z 3 6 ;2 对坝体激活单元 qn @Qd9Sf
pro bulk 1.07e7 she 6.4e6 range z 3 6 ?HIc=
ini density 1800 range z 3 6 7kn=j6I
solve *DkA$Eu3u
model elastic range z 6 9 ;3 对坝体激活单元 {CH\TmSz
pro bulk 1.07e7 she 6.4e6 range z 6 9 ,WOF)
ini density 1800 range z 6 9 lGB7(
solve 9[N' HpQ3
model elastic range z 9 14 ;4 对坝体激活单元 X_
>B7(k
pro bulk 1.07e7 she 6.4e6 range z 9 14 :kZ2N67
ini density 1800 range z 9 14 ^OG^ %
x"
solve p!'wOThO`
ini xdisp 0 ydisp 0 zdisp 0 ;将节点位移清零 p)ZlQ .d#Y
ini xvel 0 yvel 0 zvel 0 ;将节点速度清零 z@y*
jT
save 2填土石坝.sav ?l,i(I
pause oW]&]*>J
model moh range group soildam ;对坝体激活实际模型 +bm2vIh$
pro coh 1e4 fric 30 di 2 bulk 1.07e8 she 6.4e7 ten 1e5 range group soildam ； =Ak>2
solve |
7>1)
save 3填土石坝.sav iy
tSC
pause RA[` Cp"
MbnV5 b:X
; 2稀泥浆应力状态如何模拟？？？ z!Kadqns
apply nstress -3.08e5 gradient 0,0,2.2e4 range z -0.5 14.1 x 2.9 28 z i>f436-
apply szz -3.08e5 range z -2 2 x 28.1 50 hl~(& D1^
step 1 ?p5RSt
Plot surface yellow fap red ;$i9gP[|m
pause u\qyh9 s
; 监测变量 @
x*#7Y
hist unbal ;1 -lL*WA `
hist gp zdisp 3 0 14 ;2 tBBN62^ X
hist gp xdisp 3 0 14 ;3 F__>`Do l
hist gp zdisp -3 0 14 ;4 (Xq eX (s
hist gp xdisp -3 0 14 ;5 mS~3 QV
hist gp zdisp 28 0 0 ;6 pq5 )Ug
hist gp xdisp 28 0 0 ;7 o\]e}+1[o
hist gp zdisp -28 0 0 ;8 e;3$7$n Pv
hist gp xdisp -28 0 0 ;9 l O,
2
solve Lu:!vTRmw
save 4填土石坝求解完.sav j<deTK;.
|0f\>X I
1 我这里的内聚力和摩擦角很小的时候，变形特别大？？？请问如何考虑 glHag"(
2稀泥浆应力状态如何模拟？？？ any\}
wX 41R]pF
