给相交的两面墙分别施加法向约束，遇到问题

aj2000

各位好，

我想给下图所示三棱柱（无90度内角）施加法向的约束，包括一个底面及两个侧面。
我用的如下语句：
apply nvel 0.  range z -0.0001 0.0001   ; 底面
apply nvel 0.  range y -0.0001 0.0001   ; 将X轴含在内的那个侧面
macro face3 'plane dip 90 dd 300 origin 0., -0.0001,0 above'  ; 另一个侧面
apply nvel 0.  range face3

结果发现棱边上的节点并没有被约束住，而可以产生非零位移。

我在手册上找不到解决办法。请问哪位能指点一下，多谢！

[ 本帖最后由 aj2000 于 2007-8-2 19:27 编辑 ]

aj2000

自己改了一下命令，麻烦多了，但是解决问题了，
就是为了避免同一个gp 出现两个local coordinate system，如下：

; fix the gridpoint at origin
fix x y z range x -0.0001 0.0001 y -0.0001 0.0001 z -0.0001 0.0001

; fix the grid point on z axis (not including the gridpoint at origin)
fix x y range x -0.0001 0.0001 y -0.0001 0.0001 z 0.0001 100

; fix the edge on global x axis (not including the gridpoint at origin)
fix y z range x 0.0001 10 y -0.0001 0.00001 z -0.0001 0.0001

; fix the other edge
macro endface 'plane dip 90 dd 300 origin 0., -0.0001,0 above'
macro zheight 'plane dip 0 dd 300 origin 0.,0.,0.0001 below'
apply nvel 0. plane dip 90 dd 300 range endface zheight x 0.0001 100
apply dvel 0. plane dip 90 dd 300 range endface zheight x 0.0001 100

; fix the bottom face (not including the two edges)
macro endface1 'plane dip 90 dd 300 origin 0., -0.0001,0 below'
fix z range z -0.0001 0.0001 y 0.00001 100 x 0.0001 100 endface1

; fix the face on y=0 (not including the edges)
macro zheight1 'plane dip 0 dd 300 origin 0.,0.,0.0001 above'
fix y range x 0.0001 100 y -0.0001 0.0001 zheight1

; fix the other face (not including the edges)
apply nvel 0. plane dip 90 dd 300 range endface x 0.0001 100


真是复杂啊，但是自己想不到简洁的命令了，希望有更好想法的大侠指教。
谢谢

[ 本帖最后由 aj2000 于 2007-8-2 19:26 编辑 ]

benjackxu

就得用你这种局部坐标系的办法才行。而且想给象你的网格模型上的某个侧面施加法向速度边界条件（FIX）在FLAC3D中好像是实现不了的 。

aj2000

是呀， Fix 语句只适用于Global 坐标系

www

看看，学些一下