不平行于x,y,z轴的边界平面，边界条件怎么设置？？

xigang

边界平面平行于x,y,z轴，边界条件比较容易设置：fix x y range x 9.9 10.1

最近做一个挡土墙阳角的应力计算，遇到了不平行于x,y,z轴的边界平面，挡土墙也不平行于x,y,z轴

我用下面命令设置：

fix x y range plane normal x1,y1,z1 origin a,b,c above  &
                  plane normal x1,y1,z1 origin d,e,f  below

这样可以吗？？　　

第二个问题：


想在挡土墙和填土间设置接触单元，挡土墙平行x,y,z轴，比较容易设置

interface 1 face range x 9.9 10.1 y 1 10  z 1 10

但挡土墙不平行与x,y,z轴，怎么确定挡土墙的范围？？
命令改怎么写？？

下面两幅图片，红色部分为挡土墙

xigang

请版主为小弟解答下困惑

benjackxu

你是说固定面法线方向的速度？可以采用局部坐标系，但在两面交线处会与全局坐标系有冲突，似乎没有解决的办法。个人觉得这个采用应力边界条件等效要好处理得多。

ssw1216

1)斜边如果是在一个平面上采用面控制固定，如果为非平面只有采用节点控制
约束面
fix( )   range( )  plane norm ( )  origin ( )
好像也可以用
apply nvelocity 0 range plane
约束点
fix x or y or z ran (id)
用用 gpfix查看下
2) int可以这样生成
interface 1 face range plane norm ( , , ) origin (, , ) dist 0.1
                               x (, ) y ( , ) z ( , )
用 plot int查看下
试试看!

xigang

边界面是一个平面

挡墙接触面设置，主要是挡墙这个平面范围不知道如何确定
normal 和origin 倒是比较好定

xigang

边界面条件设置，FLAC借助一很薄厚度的平面来设置，
在边界面不平行于x,y,z轴的时候，法向和原点都好确定，
但这个很薄的厚度怎么控制？？

ssw1216

原帖由 xigang 于 2007-10-15 15:45 发表
边界面是一个平面

挡墙接触面设置，主要是挡墙这个平面范围不知道如何确定
normal 和origin 倒是比较好定

x (, ) y ( , ) z ( , )
是在X,Y,Z轴上投影的上下限
坐标总知道吧

benjackxu

对FLAC3D中不平行x、y、z的边界平面施加速度边界条件，我一直有疑问：
通常我们是需要该平面的法线方向速度为一定值，这时：
fix x y z ran ()
fix x ran()或者fix x yran（）就不能等效了，因为我们需要的只是约束法线方向的速度（只约束一个方向），那么就只有用局部坐标系来进行
比如
apply nvelocity 0 range plane
不过这个时候麻烦就来了，这个面的约束条件和平行x、y、z的边界平面约束条件在两面的交线处是会冲突的。如果将所有面（包括平行x、y、z的边界平面）
都采用局部坐标系
即采用：
apply nvelocity ( ) range plane
apply dvelocity 0 range plane
apply svelocity 0 range plane
的命令组合，还是会出现问题（我没有解决过 ）。不知道大家如何解决，那样岂不是所有分析模型都要把模型边界切得四四方方才行 ？

ssw1216

法线方向的速度就是两个基向量方向的合速度
故仅约束法向量分解后的两个基向量方向即可
用这个命令可以
fix x y range plane norm (-4.86,-8.56,0) origin (924.21,492.23,45) dis 0.1 &
                       x (889.25,1266.77) y (297.57,511.7) z (5,55)
dist 0.1 的概念就相当于
range x -0.1 0.1
如果怕交面处不好定义
可以设定X,Y,X的范围为小些
:D

ssw1216

偶的意思是分解处理

xigang

还有一个问题，就是，挡墙的底部设置接触面，
挡墙底部是一个直角梯形的范围（镜像一半），如果用x(,) y(,) z(,) 好像不行，
这个挡墙的底部接触面怎么定范围 ？？

FLAC3D前处理不好设置了，如果能开发成ANSYS设置边界条件就好了

feiyulx

个人认为坐标系统是整体的Ｘ－Ｙ－Ｚ，还是局部的Ｎ－Ｄ－Ｓ都没有关系，只要保证系统是静定的（甚至可以是超静定的，只是变位要产生很大的应力），我们关心部位的的作用效应的位置离边界足够远，以至于边界条件对其影响可以忽略不计，模型就是合理的。个人认为定义位移边界条件还是用整体坐标系统为好。

xigang

按照局部坐标思路，设置边界条件：

两侧使用局部坐标: apply svel 0.0 plane norm origin dist 0.1
                                apply dvel 0.0 plane norm origin dist 0.1
模型底部和背侧使用整体坐标：fix x y z range z -13.9 -14.1

结果有冲突提示：cannot have applied velocities in global and local axes at same gp
出现在两侧与底部拐角处的网格点同时使用不同坐标系

后来，我把全部边界面，背侧，两侧，底部都改成局部坐标系,使用：
两侧 apply svel/dvel0.0 plane norm (,,) orig (,,) dist 0.1
背侧:apply svel/del 0.0 plane norm -1 0 0 origin 0 0 0 range x -0.1 0.1
底部：apply svel/dvel/nvel plane norm 0 0 -1 origin (,,) range z -0.1 0.1

运行后，还是出现问题了：
gridpoint alreadly has applied velocities in another local axes system

问题还是出现在边界面交汇拐角处，同一侧 不能都使用apply svel/dvel  
否则就与另一侧的 apply svel/dvel 发生冲突，只能一侧使用apply svel
另一侧使用apply dvel ,
这很显然，另一个问题又出现了：

只要当每侧边界面不平行于x,y,z轴，每侧只能设置滚轴边界了！！！！
真的如benjackxu版主所说，模型只能切成四四方方的模型了！！！


这看来是FLAC3D的一个大的缺憾，不知道当时开发人员有没有想到这一点哦，有条件的兄弟可以像总部反应这一问题哦

我在想，如果能设置命令 fix x y z  同样也适用于fix svel dvel nvel  就好了！！！

[ 本帖最后由 xigang 于 2007-10-17 10:00 编辑 ]

benjackxu

将模型选得比关心的分析区域大得多，而且确知已测得模型边界的地应力的情况下，我觉得可以用应力边界条件替代。
就是用apply nstress ran ()
不知道这样是否可行？

bob410

这个问题最终由定论不？

kjj0501

学习了一下
以后碰到类似的问题可以参考一下了

hxm19920629

mark一下，学习了