在山区地方修建一个高14m的土石坝,土石坝拦住的是14m高的稀泥巴。我想得到的是坝体是否安全，坝体及地基土的整体应力和应变。坝体坡度是按照1：1.75的比例 6\l F
#_\**%,<
我建的模型是： 6MOwn*%5k
New J4]"@0 ?6
;title Rm3W&hQ
;Analysis of embankment dam [e,xC!2
;set logfile Analysisofembankmentdam kPOk.F%)
;set log on 8]&Fu 3M^
;1 几何建模================================================= X,@nD@
Generate zone brick & +;!^aNJ ,
p0=(0,0,0) & }vZfp5Y
p1=(28,0,0) & +Q"s!\5
p2=(0,1,0) & s]yZ<uA
p3=(0,0,14) & 3B[tbU(
P4=(28,1,0) & @}jg5}
P5=(0,1,14) & fBTNI`#
P6=(3,0,14) & ]l }v
P7=(3,1,14) & ^Zs ^
size 14 1 7 ratio 1 1 1 ;矩形 @< P;F
plot gp surface )%/ Ni^
gen zone brick & h_X'O3r
p0 0 0 -10 p1 add 50 0 0 p2 add 0 1 0 p3 add 0 0 10 & ral=`/p
size 25 1 10 ratio 1 1 1 ;矩形 L m$K R!z
attach face range z -0.1 0.1 ;联接网格 }#Up:o]A!
;镜像网格 Up-^km
Generate zone reflect normal (1,0,0) origin=(0,0,0) "/O07l1Q<
group soildam range z 0 14 0<A*I{,4L
group groundsill1 range z -2 0 H^*[TX=#[
group groundsill2 range z -6 -2 DB"z93Mr<K
group groundsill3 range z -10 -6 a@ |`!<5
; 2 指定地基土为Mohr Coulomb 模型============================================= 3I'7+?@@l
model null range z 0 15 ;对坝体设置空单元 k=t{o
model moh range z -11 0 ;对地基土网格进行赋值 W S vhC
pro coh 1e10 fric 30 di 2 bulk 6.7e6 she 3.1e6 ten 1e10 range group groundsill1 aIT0t0.
ini density 1800 range group groundsill1 0dx%b677d
pro coh 1e10 fric 10 di 0 bulk 1.11e7 she 3.70e6 ten 1e10 range group groundsill2 ;U0w<>4L
ini density 1700 range group groundsill2 11S{XbU
pro coh 1e10 fric 30 di 3 bulk 3.33e8 she 1.11e8 ten 1e10 range group groundsill3 Z>`frL
ini density 1700 range group groundsill3 5 3+C;]J
; 3 边界条件================================================================= vzAY+EEx
fix x range x 49.9 50.1 PtPx (R3
fix x range x -50.1 -49.9 _ y'g11 \
fix y range y -0.1 0.1 e}gGl<((g
fix y range y 0.9 1.1 {'P?wv
fix x y z range z -10.1 -9.9 ko!]vHB9`
; 4 初始应力状态========================================================= M$v\7vBgO!
set grav 0 0 -9.81 NxH%%>o>
set large xE_~.EoB
ini szz -300e3 grad 0 0 18000 1"A"AMZf
ini sxx -150e3 grad 0 0 9000 $XyDw|z[
ini syy -150e3 grad 0 0 9000 ty;a!yjC
plot cont szz outline on ,o2x,I
solve ;第一次求解 ).Z
U0fV
pro coh 1e4 fric 30 di 2 bulk 1.07e7 she 6.4e6 ten 1e5 range group groundsill1 =J-&usX
pro coh 0.8e4 fric 20 di 0 bulk 5.14e6 she 1.97e6 ten 1e5 range group groundsill2 M;g"rpM
pro coh 2.5e4 fric 30 di 3 bulk 3.33e8 she 1.11e8 ten 14433.7 range group groundsill3 X.e4pLwGK
solve ;第二次求解 abe5 As r
ini xdisp 0 ydisp 0 zdisp 0 ;将节点位移清零 N#OO{`":Z`
ini xvel 0 yvel 0 zvel 0 ;将节点速度清零 WYvcN8F
;plot con xdis f#38QP-T
save 1地基初始应力场.sav DBfq9%J _
;5 土石坝施工过程模拟================================================= 4\U"e *
model elastic range z 0 3 ;1 对坝体激活单元 gW G>}M@
pro bulk 1.07e7 she 6.4e6 range z 0 3 .$&vSOgd(
ini density 1800 range z 0 3 GOCe&?
solve M%13b$i~f
;plot con xdis u~t% GIg
model elastic range z 3 6 ;2 对坝体激活单元 YFcMU5_F
pro bulk 1.07e7 she 6.4e6 range z 3 6 8N<2R T8W
ini density 1800 range z 3 6 o
-x=/b
solve 9-Qu5 L~
model elastic range z 6 9 ;3 对坝体激活单元 d<_IC7$u>
pro bulk 1.07e7 she 6.4e6 range z 6 9 06af{FXsGb
ini density 1800 range z 6 9 fn //j7 j
solve sEb*GF*.V
model elastic range z 9 14 ;4 对坝体激活单元 V#&S&dn
pro bulk 1.07e7 she 6.4e6 range z 9 14 skLr6Cs|
ini density 1800 range z 9 14 nqY arHi
solve =dQ46@
ini xdisp 0 ydisp 0 zdisp 0 ;将节点位移清零 K?S5C8
ini xvel 0 yvel 0 zvel 0 ;将节点速度清零 V'N]u (^
save 2填土石坝.sav #FGj)pu
pause -!j5j:RR
model moh range group soildam ;对坝体激活实际模型 |@D%y &
pro coh 1e4 fric 30 di 2 bulk 1.07e8 she 6.4e7 ten 1e5 range group soildam ； E_xCRfw_i]
solve GKXd"8z]
save 3填土石坝.sav Y0Rk:Njc
pause GQ9H>Ssz
jJkc vC8d
; 2稀泥浆应力状态如何模拟？？？ 78'3&,+si
apply nstress -3.08e5 gradient 0,0,2.2e4 range z -0.5 14.1 x 2.9 28 I&&[ ':
apply szz -3.08e5 range z -2 2 x 28.1 50 ?Ql<s8
step 1 AbMf8$$3SH
Plot surface yellow fap red *kGk.a =
pause . P! pC
; 监测变量 p ^I#9(PT
hist unbal ;1 DoE N`K\U
hist gp zdisp 3 0 14 ;2 ynq^ztBVe
hist gp xdisp 3 0 14 ;3 K{G\=yJ((
hist gp zdisp -3 0 14 ;4 7@!ne&8Z?
hist gp xdisp -3 0 14 ;5 Ql nI &o
hist gp zdisp 28 0 0 ;6 }(UU~V
hist gp xdisp 28 0 0 ;7 ;`Wh^Qgi
hist gp zdisp -28 0 0 ;8 %\ !3tN
hist gp xdisp -28 0 0 ;9 H3 |x
solve ]#O~lq
save 4填土石坝求解完.sav h"$], =
{@67'jL
1 我这里的内聚力和摩擦角很小的时候，变形特别大？？？请问如何考虑 I|P#|0< 2
2稀泥浆应力状态如何模拟？？？ ,h9N, bIQg
Tml>>O
