- 积分
- 1
- 注册时间
- 2007-5-28
- 仿真币
-
- 最后登录
- 1970-1-1
|
采用论坛中海棠和君之风的程序做了一个小实验。
实验目的:不同和单元大小对空隙水压力的影响、
实验步骤:(1)利用海棠提供的程序,在借用fish和surfer程序的基础上,生成具有起伏的地表层及地表以下的基岩层。(2)为了生成曲形水面,在生成起伏地形的模型时,共生成两套模型,其中一套为实际模型zhengchang2,另外一套的边界比实际模型稍大,暂定为kuodazu5。在实际模型zhengchang2中,基岩层定义为组 1 ,地表定义为组 2 ,对应的文件为“zhengchang2.flac3d”。在kuodazu5模型中,基岩层定义为组 1_da,地表定义为组 2_da。(3)调用“小水面-13.txt”文本文档,完成水面的生成和水压的生成。
出现的问题:(1)由于在基岩组的网格划分上存在不同大小的单元,在实际模型zhengchang2和扩大模型kuodazu5中,垂直方向 z 的变化范围在-20m到-10m的单元大小为2m×1m,而在-10m到基岩组顶面的范围内单元大小则为4m×2m。由此造成了图17出现的问题,在单元大小不同的地方水压出现了斑点状。
以下是模型的调用代码:n
;================
;导入网格数据以生成水面
impgrid kuodazu5.flac3d
;为建立水面设置“虚拟界面”而分组为水上部分和水下部分
group water_above range group 2_da any
group water_below range group water_above not
gen separate water_below
interface 1 wrap water_below water_above
set grav 0 0 -10
water den 1000
;==================================
;利用遍历界面单元节点生成水面
def water_table
p_i=i_head
p_ie=i_elem_head(p_i) ;界面1单元指针赋予p_ie
loop while p_ie # null
;返回界面单元的三个顶点的地址
p_gp1=ie_vert(p_ie,1)
p_gp2=ie_vert(p_ie,2)
p_gp3=ie_vert(p_ie,3)
;以这三个界面单元顶点为水面的三个顶点,生成水面
;========interface的节点和网格节点的地址是分别存放的
x1=in_pos(p_gp1,1)-0.001
y1=in_pos(p_gp1,2)-0.001
z1=in_pos(p_gp1,3)-0.001
x2=in_pos(p_gp2,1)+0.001
y2=in_pos(p_gp2,2)+0.001
z2=in_pos(p_gp2,3)+0.001
x3=in_pos(p_gp3,1)+0.001
y3=in_pos(p_gp3,2)+0.001
z3=in_pos(p_gp3,3)+0.001
command
water table face x1,y1,z1 x2,y2,z2 x3,y3,z3
endcommand
p_ie=ie_next(p_ie)
endloop
end
;删除单元组water_below和组water_above,为导入实际的材料分组做准备
model null range group water_below
model null range group water_above
;===========================
;建立实际材料分组部分
impgrid zhengchang2.flac3d
attach face
;建立真实的“界面”
gen separate 1
interface 2 wrap 1 2
;生成水面
water_table
;删除虚拟界面单元
int 1 dele
;彻底删除原先的分组,减少存储空间
dele range group water_below
dele range group water_above
;施加直边约束
fix x y z rang z -19.9 -20.1
fix x range x 3.9 4.1
fix x range x 63.9 64.1
fix y range y 1.9 2.1
fix y range y 47.9 48.1
plot set rotation 130 290 0 ang 22.5 mag 1.25
plot add axes red scal 0.5
plo add water lgre
plo add sketch red
plo sho
save 小水面-13.sav
调用文件在附件中,请各位不吝赐教,在此先谢过了。
生成的效果图如下所示。
图1 水面位置图
图2 整个模型水压效果图1
图3 整个模型水压效果图2
图4 垂直方向-10m及以上高程水压效果图(单元大小4m×2m)
图5 垂直方向-20m到-10m顶面水压效果图(单元大小2m×1m)
图6 垂直方向-20m到-10m底面水压效果图(单元大小2m×1m)
|
本帖子中包含更多资源
您需要 登录 才可以下载或查看,没有账号?注册
×
|