在山区地方修建一个高14m的土石坝,土石坝拦住的是14m高的稀泥巴。我想得到的是坝体是否安全，坝体及地基土的整体应力和应变。坝体坡度是按照1：1.75的比例 hh$V[/iK
B>!OW2q0D
我建的模型是： pod=|(c
New GCEcg&s=\S
;title T`?{Is['(
;Analysis of embankment dam @ 7WWoy
;set logfile Analysisofembankmentdam 5VXI/Lw#
;set log on ,{tK{XpS
;1 几何建模================================================= :hA=(iz
Generate zone brick & ^#Ruw?D
p0=(0,0,0) & <NVSF6`
p1=(28,0,0) & a#_=c>h;
p2=(0,1,0) & '%} k"&t$i
p3=(0,0,14) & O `pqS\H
P4=(28,1,0) & /%P|<[<
[
P5=(0,1,14) & -#e3aXe
P6=(3,0,14) & %( tu<
P7=(3,1,14) & @+,J^[ y
size 14 1 7 ratio 1 1 1 ;矩形 K< ;I*cAX
plot gp surface ;]/emw=a
gen zone brick & +wwb+aG6{
p0 0 0 -10 p1 add 50 0 0 p2 add 0 1 0 p3 add 0 0 10 & M zWVsV
size 25 1 10 ratio 1 1 1 ;矩形 Vp5i i]B4
attach face range z -0.1 0.1 ;联接网格 Tw{}Ht_Qq
;镜像网格 \*(A1Vk
Generate zone reflect normal (1,0,0) origin=(0,0,0) `wzb}"gLsM
group soildam range z 0 14 x bsk
group groundsill1 range z -2 0 Kpx(x0^2
group groundsill2 range z -6 -2 ,
Ac
gsC
group groundsill3 range z -10 -6 S T *\ Q
; 2 指定地基土为Mohr Coulomb 模型============================================= 42{\u 08Z
model null range z 0 15 ;对坝体设置空单元 U< aT%^_
model moh range z -11 0 ;对地基土网格进行赋值 .h0@Vs
pro coh 1e10 fric 30 di 2 bulk 6.7e6 she 3.1e6 ten 1e10 range group groundsill1 t\/i9CBn
ini density 1800 range group groundsill1 &ml7368@
pro coh 1e10 fric 10 di 0 bulk 1.11e7 she 3.70e6 ten 1e10 range group groundsill2 -F+
)N$CW
ini density 1700 range group groundsill2 {4%B^+}T
pro coh 1e10 fric 30 di 3 bulk 3.33e8 she 1.11e8 ten 1e10 range group groundsill3 )e1&[0
ini density 1700 range group groundsill3 bJ$6[H-:
; 3 边界条件================================================================= M$O*@])
fix x range x 49.9 50.1
,'KS:`m!
fix x range x -50.1 -49.9 a(qij&>
fix y range y -0.1 0.1 L93PDp4v
fix y range y 0.9 1.1 2Gw2k8g&
fix x y z range z -10.1 -9.9 q3CcXYY
; 4 初始应力状态========================================================= -{^I T`
set grav 0 0 -9.81 ry2ZVIFa
set large G1vg2'A
ini szz -300e3 grad 0 0 18000 VX;br1$X
ini sxx -150e3 grad 0 0 9000 2EU((Q`>=(
ini syy -150e3 grad 0 0 9000 lh N2xg5x
plot cont szz outline on Hw(_l,Xf
solve ;第一次求解 U i^~A
pro coh 1e4 fric 30 di 2 bulk 1.07e7 she 6.4e6 ten 1e5 range group groundsill1 D>U b)i
pro coh 0.8e4 fric 20 di 0 bulk 5.14e6 she 1.97e6 ten 1e5 range group groundsill2 e)#O -y
pro coh 2.5e4 fric 30 di 3 bulk 3.33e8 she 1.11e8 ten 14433.7 range group groundsill3 A$H;2T5N
solve ;第二次求解 mw[T[
ini xdisp 0 ydisp 0 zdisp 0 ;将节点位移清零 91Fx0(
ini xvel 0 yvel 0 zvel 0 ;将节点速度清零 ! b=jD;<
;plot con xdis !P_8D*^9
save 1地基初始应力场.sav _ yu d
;5 土石坝施工过程模拟================================================= TVVr<r
model elastic range z 0 3 ;1 对坝体激活单元 W$I^Ej}>$
pro bulk 1.07e7 she 6.4e6 range z 0 3 4IT`8n~
ini density 1800 range z 0 3 `[&) X
solve ;$\?o
;plot con xdis 08+cN T
model elastic range z 3 6 ;2 对坝体激活单元 aR2 Vvo
pro bulk 1.07e7 she 6.4e6 range z 3 6 ;a Y.CgX
ini density 1800 range z 3 6 b%TS37`^[
solve E6=JL$"
model elastic range z 6 9 ;3 对坝体激活单元 I"4Lma
pro bulk 1.07e7 she 6.4e6 range z 6 9 ;VL v2J*
ini density 1800 range z 6 9 =}$YZ uzmU
solve <fZ?F=
model elastic range z 9 14 ;4 对坝体激活单元 Y%i=u:}fm
pro bulk 1.07e7 she 6.4e6 range z 9 14 |\ Zs oA
ini density 1800 range z 9 14 Hju7gP=y}
solve =Fd !wkB'{
ini xdisp 0 ydisp 0 zdisp 0 ;将节点位移清零 T[z}^"
ini xvel 0 yvel 0 zvel 0 ;将节点速度清零 B*zR/?U^
save 2填土石坝.sav t=K;/ 1
pause (D?%(f
model moh range group soildam ;对坝体激活实际模型 h0=Q .Yz6
pro coh 1e4 fric 30 di 2 bulk 1.07e8 she 6.4e7 ten 1e5 range group soildam ； 0c$0<2D%
solve <#[_S$ 54
save 3填土石坝.sav y,i ~w |4
pause BkH- d z
J#```cB
; 2稀泥浆应力状态如何模拟？？？ BJ
c'4>
apply nstress -3.08e5 gradient 0,0,2.2e4 range z -0.5 14.1 x 2.9 28
sp/l-a
apply szz -3.08e5 range z -2 2 x 28.1 50 @AvD V$F
step 1 qG=>eRR
Plot surface yellow fap red {x-iBg9#l2
pause /9K,W)h_
; 监测变量 o9j*Yz
hist unbal ;1 0L
^WTq
hist gp zdisp 3 0 14 ;2 =)2sehU/
hist gp xdisp 3 0 14 ;3 4)cQU.(*k
hist gp zdisp -3 0 14 ;4 XMR$I&;G8
hist gp xdisp -3 0 14 ;5 RW?F{Jy{
hist gp zdisp 28 0 0 ;6 iq25|{1$
hist gp xdisp 28 0 0 ;7 =u*\P!$
hist gp zdisp -28 0 0 ;8 iEA$`LhO\A
hist gp xdisp -28 0 0 ;9 DR d|m<Z
solve t,YnweH
save 4填土石坝求解完.sav 6 a E:v R2
ir72fSe
1 我这里的内聚力和摩擦角很小的时候，变形特别大？？？请问如何考虑 /u h?F
2稀泥浆应力状态如何模拟？？？ {c7ZA %T~R
]ZkR~?
