- 积分
- 0
- 注册时间
- 2011-3-9
- 仿真币
-
- 最后登录
- 1970-1-1
|
楼主 |
发表于 2011-5-24 20:35:50
|
显示全部楼层
来自 北京
new
;1 几何建模============================
generate zone brick & ;建立区域1
p0=(0,0,0)&
p1=(55,0,0)&
p2=(0,1,0)&
p3=(0,0,20.5)&
p4=(55,1,0)&
p5=(0,1,20.5)&
p6=(3.5,0,20.5)&
p7=(3.5,1,20.5)&
size 30 1 20 ratio 1 1 1
plot gp surface
gen zone brick & ;建立区域2
p0 0 0 -55 p1 55 0 -55 p2 0 1 -55 p3 0 0 0 &
size 30 1 30 ratio 1 1 0.98
gen zone brick & ;建立区域3
p0 55 0 -55 p1 110 0 -55 p2 55 1 -55 p3 55 0 0 &
size 20 1 30 ratio 1.04 1 0.98
generate zone reflect normal(1,0,0) origin=(0,0,0) ;镜像网格
attach face ;联接网格的面
group soildam range z 0 20.5
group groundsill1 range z -20 0
group groundsill2 range z -55 -20
;2边界条件===========================================
fix x range x 109.9 110.1 ;固定模型右侧边界的X方向速度
fix x range x -109.9 -110.1 ;固定模型左侧边界的X方向速度
fix x y z range z -54.9 -55.1 ;固定模型底部边界的三个方向速度
fix y ;固定模型所有节点Y方向速度
;3指定地基土为Mohr Coulomb模型======================
model null range z 0 20.5 ;对坝体设置空单元
model moh range z -55 0 ;对地基土网格进行赋值
pro bulk 30e6 shear 6.5e6 coh 10e10 ten 1e10 range group groundsill1
ini density 1850 range group groundsill1
pro bulk 4.67e9 shear 2.8e9 coh 10e10 ten 1e10 range group groundsill2
ini density 2540 range group groundsill2
;4初始应力状态=======================================
set grav 0 0 -9.81
solve ;第一次求解
pro bulk 30e6 shear 6.5e6 coh 2.0e4 fric 18 range group groundsill1 ;赋真实参数
pro bulk 4.67e9 shear 2.8e9 coh 1.0e6 fric 32 range group groundsill2
solve ;第二次求解
save 地基初始应力.sav |
|