求教：不同单元大小空隙水压力斑点图的处理

durongjun2003

采用论坛中海棠和君之风的程序做了一个小实验。
实验目的：不同和单元大小对空隙水压力的影响、
实验步骤：（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）