施加径向约束问题

18596764675

问题：怎么样对一个竖向放置（轴线为Y轴，端面在XZ平面），圆心不在坐标原点的空心圆柱体，类似于图1所示，目的是在外表面施加垂直于面的约束，即径向约束。


用下面方法操作后显示的约束，如图2所示，不知道问题出现在哪里，请问下正确的操作应该是怎么样，望解答，十分感谢！











1.激活柱坐标系：
CSYS,12.选择要施加径向约束曲面上的所有节点 ：
FLST,5，  …           
FITEM,5, …
FITEM,5, …
……
ASEL,S, , ,P51X
NSLA,S,1
NPLOT
3.旋转节点坐标系，使之与激活的坐标系平行：
FLST,2  …
FITEM,2, …
FITEM,2, …
……
NROTAT,P51X
4.对所选节点施加X方向约束：
FLST,  …
FITEM,2, …
FITEM,2, …
/GO
D,P51X, , , , , ,UX, , , , ,

htbbzzg

  需要先定义一个局部圆柱坐标系，其 Z 轴平行于总体 Y 轴，且通过圆柱体的的圆心。然后将要施加径向约束曲面上的所有节点的节点坐标系转换为该局部圆柱坐标系，再在这些节点上施加 Ux 的约束。
  你的命令流似乎是这样做的，只是缺少一些内容，不知其中有什么问题。
  另外，也可能是显示问题，在施加约束前，增加一个命令：
    DSYS,12
  试试。

njweiwei2

htbbzzg 发表于 2018-12-6 08:48
需要先定义一个局部圆柱坐标系，其 Z 轴平行于总体 Y 轴，且通过圆柱体的的圆心。然后将要施加径向约束曲 ...

DSYS是显示坐标系命令，只影响显示结果；比如直角坐标系里面的弧线在圆柱坐标系里成了直线；

建议检查一下第12坐标系的类型，确定是圆柱坐标系，而不是直角；另外，可以试试先把边界条件转换到节点上（虽然从命令上看是直接加载到节点上的，但是就当预防万一吧）。具体命令为：
sbctran
/replo

18596764675

htbbzzg 发表于 2018-12-6 08:48
需要先定义一个局部圆柱坐标系，其 Z 轴平行于总体 Y 轴，且通过圆柱体的的圆心。然后将要施加径向约束曲 ...

按照您说的操作方法

local,11,1,4.8,0,-20,0,90,0
创建一个局部柱坐标系11，其 Z 轴平行于总体 Y 轴，且通过圆柱体的的圆心，圆心在总体坐标系中为4.8,0,-20
出来的结果如图片所示，看着有点奇怪，不知道怎么回事

18596764675

njweiwei2 发表于 2018-12-6 09:40
DSYS是显示坐标系命令，只影响显示结果；比如直角坐标系里面的弧线在圆柱坐标系里成了直线；

建议检查一 ...

好的，谢谢您的解答

18596764675

htbbzzg 发表于 2018-12-6 08:48
需要先定义一个局部圆柱坐标系，其 Z 轴平行于总体 Y 轴，且通过圆柱体的的圆心。然后将要施加径向约束曲 ...

模型如图，上表面自由，下表面固定，其他四周的面约束其各自的法向。
命令流如下，不知道正不正确，麻烦您帮忙看一下，谢谢了
/PREP7  
et,1,45
et,2,200  
KEYOPT,2,1,6
KEYOPT,2,2,0
MPTEMP,,,,,,,,  
MPTEMP,1,0  
MPDATA,EX,1,,1200000000
MPDATA,PRXY,1,,0.35
MPTEMP,,,,,,,,  
MPTEMP,1,0  
MPDATA,DENS,1,,2500
RECTNG,4,-4,2.4,-3.6,   
K,100,0,0,0,
K,101,0,0,-20,  
K,102,4.8,0,-24.8,  
K,103,24.8,0,-24.8,
K,104,4.8,0,-20,
LSTR,     100,     101  
LSTR,     102,     103  
LSTR,     101,     104  
FLST,2,1,8  
FITEM,2,4.8,0,-20   
WPAVE,P51X  
wpro,,90.000000,
!*  
LARC,102,101,104,4.8,   
LDELE,      7   
TYPE,   2   
MAT,       1
REAL,   
ESYS,       0   
SECNUM,
!*  
FLST,5,4,4,ORDE,2   
FITEM,5,1   
FITEM,5,-4  
CM,_Y,LINE  
LSEL, , , ,P51X
CM,_Y1,LINE
CMSEL,,_Y   
!*  
LESIZE,_Y1, , ,32, , , , ,1
!*  
FLST,5,3,4,ORDE,3   
FITEM,5,5   
FITEM,5,-6  
FITEM,5,8   
CM,_Y,LINE  
LSEL, , , ,P51X
CM,_Y1,LINE
CMSEL,,_Y   
!*  
LESIZE,_Y1, , ,80, , , , ,1
!*  
CM,_Y,AREA  
ASEL, , , ,       1
CM,_Y1,AREA
CHKMSH,'AREA'   
CMSEL,S,_Y  
!*  
MSHKEY,1
AMESH,_Y1   
MSHKEY,0
!*  
CMDELE,_Y   
CMDELE,_Y1  
CMDELE,_Y2  
!*  
LPLOT   
TYPE,   1   
MAT,       1
REAL,   
ESYS,       0   
SECNUM,
!*  
FLST,8,3,4  
FITEM,8,5   
FITEM,8,8   
FITEM,8,6   
VDRAG,       1, , , , , ,P51X   
WPCSYS,-1,0
WPSTYLE,,,,,,,,0
/REPLOT
FLST,2,3,5,ORDE,3   
FITEM,2,1   
FITEM,2,13  
FITEM,2,15  
!*  
/GO
DA,P51X,UZ,
FLST,2,3,5,ORDE,3   
FITEM,2,3   
FITEM,2,5   
FITEM,2,16  
!*  
/GO
DA,P51X,UX,
FLST,2,3,5,ORDE,3   
FITEM,2,2   
FITEM,2,7   
FITEM,2,12  
!*  
/GO
DA,P51X,ALL,
local,11,1,4.8,0,-20,0,90,0
FLST,5,2,5,ORDE,2   
FITEM,5,8   
FITEM,5,10  
ASEL,S, , ,P51X
NSLA,S,1
NPLOT   
FLST,2,5346,1,ORDE,18   
FITEM,2,1090
FITEM,2,-1091   
FITEM,2,1123
FITEM,2,-1155   
FITEM,2,1187
FITEM,2,-1217   
FITEM,2,88210   
FITEM,2,-88211  
FITEM,2,88243   
FITEM,2,-88275  
FITEM,2,88307   
FITEM,2,-88337  
FITEM,2,89299   
FITEM,2,-89614  
FITEM,2,92064   
FITEM,2,-94512  
FITEM,2,96962   
FITEM,2,-99410  
NROTAT,P51X
FLST,2,5346,1,ORDE,18   
FITEM,2,1090
FITEM,2,-1091   
FITEM,2,1123
FITEM,2,-1155   
FITEM,2,1187
FITEM,2,-1217   
FITEM,2,88210   
FITEM,2,-88211  
FITEM,2,88243   
FITEM,2,-88275  
FITEM,2,88307   
FITEM,2,-88337  
FITEM,2,89299   
FITEM,2,-89614  
FITEM,2,92064   
FITEM,2,-94512  
FITEM,2,96962   
FITEM,2,-99410  
!*  
/GO
D,P51X, , , , , ,UX, , , , ,
ALLSEL,ALL  
VPLOT

njweiwei2

18596764675 发表于 2018-12-6 21:32
模型如图，上表面自由，下表面固定，其他四周的面约束其各自的法向。
命令流如下，不知道正不正确，麻烦 ...

正确。
去除底面节点以后看到的约束（即确认侧面是否添加的是法向约束）如图1：
底面节点的约束如图2（此处为了显示清楚，只显示了一部分）：

18596764675

njweiwei2 发表于 2018-12-7 08:27
正确。
去除底面节点以后看到的约束（即确认侧面是否添加的是法向约束）如图1：
底面节点的约束如图2（此 ...

非常感谢您的回复，您的意思是说我发的命令流和操作都没有问题么，那我怎么挑取任意一条弧线（比如上表面的弧线）显示径向约束都如同我上面发的图片里面显示的那样，不是完全指向径向的按您说的也查看了下，如下图所示，图中红色框框部分显示很奇怪

njweiwei2

18596764675 发表于 2018-12-7 08:33
非常感谢您的回复，您的意思是说我发的命令流和操作都没有问题么，那我怎么挑取任意一条弧线（比如上表面 ...

那就奇怪了，我用的就是你的命令流呢。
试试看：

sbctran
/replo