- 积分
- 0
- 注册时间
- 2014-2-26
- 仿真币
-
- 最后登录
- 1970-1-1
|
发表于 2018-7-15 11:24:42
|
显示全部楼层
来自 中国
可以通过在地形面加nstress(就是给起伏面加垂直于面的应力),这样就起伏面上有应力了,然后判定哪些节点力不为零就提取哪些节点的编号和坐标就可以了啊
5.2.2 搜索临空面
临空面是指岩土体滑动时自由空间的边界面,如斜坡的坡面,硐室的边墙、坝基下游河床面等。在FLAC3D建模中一般的水平或竖直的边界面很容易通过range限定,但对于起伏不规则的边界面比如山地、不规则斜坡的坡面就要通过FLAC3D内置FISH语言搜寻。寻找临空面的目的有三:
(1)、检查网格、判断模型的对错。
(2)、在临空面上施加载荷。
(3)、便于前后处理。
一般搜寻的方法也有三个:
(1)、寻找面在单元体遍历中出现的次数,边界面只能出现一次,超过一次即为岩体内部单元表面而非临空面。
(2)、利用interface接触面单元或shell、liner结构单元会自动捕捉临空面的特性搜寻临空面。
(3)、给模型施加静水压力的apply nstress命令也可以自动搜寻临空面,施加的静水压力除以单元体表面积的应力被平均分配到单元体节点上,通过只在所需的临空面施加静水压力,然后寻找应力不为零的节点即找到临空面所有节点编号,临空面也就找到了。
上面三种方法中,第一种适用于比较麻烦,需要比较清楚FLAC3D模型网格信息组成规律。第二种和第三种方法都比较适合于寻找底面和四周侧面全是平面而上表面为曲面的上临空面搜寻。第二种方法有时在边界处有些单元体节点因为命令不熟悉原因可能没有包括进来,建立好interface后仍需要fish函数提取相应的节点编号。
第三种方法具体代码如下,由于静水压力是被平均分配到节点上的,只要靠近边界处单元表面被施加了静水压力那么边界处单元节点就会被分配到应力,通过检测应力不为零的节点即可找到临空面所有节点编号。
|
|