在山区地方修建一个高14m的土石坝,土石坝拦住的是14m高的稀泥巴。我想得到的是坝体是否安全，坝体及地基土的整体应力和应变。坝体坡度是按照1：1.75的比例 /&Vgo ~.J
,`nl";Zc
我建的模型是： O,A}p:Pgs
New 7J$5dFV2
;title ,Z 1W3;O
;Analysis of embankment dam |&xjuBC
;set logfile Analysisofembankmentdam |= U(8t
;set log on "RG #e +
;1 几何建模================================================= /*2sg>e'QF
Generate zone brick & @[] A&)B
p0=(0,0,0) & R$xk cg2(
p1=(28,0,0) & nbM7 >tnsk
p2=(0,1,0) & f8X /kz
p3=(0,0,14) & YTo^Q&
P4=(28,1,0) & 5q>u]n9]
P5=(0,1,14) & t
9Dr%#
P6=(3,0,14) & JJn+H&[B
P7=(3,1,14) & AG==A&d>$
size 14 1 7 ratio 1 1 1 ;矩形 z,#3YC{'
plot gp surface 9e xHR&>{
gen zone brick & ;th]/ G
p0 0 0 -10 p1 add 50 0 0 p2 add 0 1 0 p3 add 0 0 10 & 8OW504AD
size 25 1 10 ratio 1 1 1 ;矩形 9b}AZ]$
attach face range z -0.1 0.1 ;联接网格 8G0DuMI5
;镜像网格 A[.5B i
Generate zone reflect normal (1,0,0) origin=(0,0,0) T.j&UEsd
group soildam range z 0 14 TPeBb8v 8D
group groundsill1 range z -2 0 I-`qo7dQ_S
group groundsill2 range z -6 -2 O&c~7tM%
group groundsill3 range z -10 -6 a vI
; 2 指定地基土为Mohr Coulomb 模型============================================= Y
=BXV7\
model null range z 0 15 ;对坝体设置空单元 B%e#u.'6
model moh range z -11 0 ;对地基土网格进行赋值 7h3JH
pro coh 1e10 fric 30 di 2 bulk 6.7e6 she 3.1e6 ten 1e10 range group groundsill1 fpK`
ini density 1800 range group groundsill1 eS{ xma
pro coh 1e10 fric 10 di 0 bulk 1.11e7 she 3.70e6 ten 1e10 range group groundsill2 3ZKaqwK
ini density 1700 range group groundsill2 T1}9^3T?{
pro coh 1e10 fric 30 di 3 bulk 3.33e8 she 1.11e8 ten 1e10 range group groundsill3 F I I>6c
ini density 1700 range group groundsill3 1>1i i
; 3 边界条件================================================================= raPUx _$PH
fix x range x 49.9 50.1 m K2M1r
fix x range x -50.1 -49.9 \p&a c&]
fix y range y -0.1 0.1 NXyuv7%5=
fix y range y 0.9 1.1 > =Z@)PAe
fix x y z range z -10.1 -9.9 8*V8B=q}K
; 4 初始应力状态========================================================= D;l)&"| r?
set grav 0 0 -9.81 G r
a(DGX
set large v8L&F9
o
ini szz -300e3 grad 0 0 18000 d{LQr}_o$$
ini sxx -150e3 grad 0 0 9000 >&kb|)
ini syy -150e3 grad 0 0 9000 +^6v%z
plot cont szz outline on b[I;6HW
solve ;第一次求解 $ *u{i4b
pro coh 1e4 fric 30 di 2 bulk 1.07e7 she 6.4e6 ten 1e5 range group groundsill1 /`Wd+
pro coh 0.8e4 fric 20 di 0 bulk 5.14e6 she 1.97e6 ten 1e5 range group groundsill2 9ywPWT[^
pro coh 2.5e4 fric 30 di 3 bulk 3.33e8 she 1.11e8 ten 14433.7 range group groundsill3 |NiW r1&i0
solve ;第二次求解 6$ JRV
ini xdisp 0 ydisp 0 zdisp 0 ;将节点位移清零 b@&ydgmaQ
ini xvel 0 yvel 0 zvel 0 ;将节点速度清零 J& IFn/JK$
;plot con xdis =>7\s}QZ
save 1地基初始应力场.sav @Fl&@ $
;5 土石坝施工过程模拟================================================= 4gNF;
model elastic range z 0 3 ;1 对坝体激活单元 .C2.j[>
pro bulk 1.07e7 she 6.4e6 range z 0 3 pd rF/U+
ini density 1800 range z 0 3 #(`@D7S"
solve sN
`NZyG
;plot con xdis jSj
(ZU6
model elastic range z 3 6 ;2 对坝体激活单元 .7nr :P
pro bulk 1.07e7 she 6.4e6 range z 3 6 W2a 9P_
ini density 1800 range z 3 6 G *f5B
solve #VwA?$4g`
model elastic range z 6 9 ;3 对坝体激活单元 {=NHidi ~
pro bulk 1.07e7 she 6.4e6 range z 6 9 $K,6!FyBa
ini density 1800 range z 6 9 oV c
l (
solve ) n[ oP%
model elastic range z 9 14 ;4 对坝体激活单元 Kb# }f/
pro bulk 1.07e7 she 6.4e6 range z 9 14 sO!YM5v8
ini density 1800 range z 9 14 Yjr6/&ML
solve QzQTE-SQ
ini xdisp 0 ydisp 0 zdisp 0 ;将节点位移清零 *<**rY*
ini xvel 0 yvel 0 zvel 0 ;将节点速度清零 B !hrr
save 2填土石坝.sav &dV|~xA6N
pause U*'
YGv
model moh range group soildam ;对坝体激活实际模型 t*? CD.S
pro coh 1e4 fric 30 di 2 bulk 1.07e8 she 6.4e7 ten 1e5 range group soildam ； 62A b4!
solve F<UEipe/N
save 3填土石坝.sav ~!,Q<?
pause R":nG7o
k_=~ObA$g
; 2稀泥浆应力状态如何模拟？？？ ~la=rh 3
apply nstress -3.08e5 gradient 0,0,2.2e4 range z -0.5 14.1 x 2.9 28 Ass8 c]H@
apply szz -3.08e5 range z -2 2 x 28.1 50 tN5b rf
step 1 3d{v5. C#X
Plot surface yellow fap red \sIRV}Tk}N
pause 1MV\
^l_
; 监测变量 w:z@!<
hist unbal ;1 s1!_zf_
hist gp zdisp 3 0 14 ;2 .bm#|X)RO
hist gp xdisp 3 0 14 ;3 =9
TAs? =
hist gp zdisp -3 0 14 ;4 \kpk-[W*x{
hist gp xdisp -3 0 14 ;5 !8|] R
hist gp zdisp 28 0 0 ;6 at/v.U |F
hist gp xdisp 28 0 0 ;7 ^;2dZgJ 4^
hist gp zdisp -28 0 0 ;8 G~{xTpL
hist gp xdisp -28 0 0 ;9 PUUB n"U-
solve 9 GdrJ~h
save 4填土石坝求解完.sav %^l77 :O
bT @7&
1 我这里的内聚力和摩擦角很小的时候，变形特别大？？？请问如何考虑 #_+T@|r
2稀泥浆应力状态如何模拟？？？ |f^/((: D
=#i#IF42?
