请教梁单元自由度的释放问题

wazon

版主，大侠，你们好！我有一个问题想请教，问题如下：现在想用刚性单元和梁单元来模拟铰接， 蜘蛛网用RBE3单元，中间的用梁单元CBEAM，想释放掉转动自由度。不过释放了转动的话，就报错了，不释放就不报错，在hypermesh里面设置的。


释放之后的文件：（注：本想释放转动，见第二行的4，4：不知道这样是不是错了？）
                    1-        CBEAM   507782  13      450030  450018  0.0     1.0     0.0             ++00000B
                    2-        ++00000B4       4                                                        +00000C
                    3-        CBEAM   507783  13      450018  450029  0.0     1.0     0.0             ++000009
                    4-        ++0000094       4                                                        +00000A
                    5-        CBEAM   507784  13      450027  450015  0.0     1.0     0.0             ++00000F
                    6-        ++00000F4       4                                                        +00000G
                    7-        CBEAM   507785  13      450015  450026  0.0     1.0     0.0             ++00000D
                    8-        ++00000D4       4                                                        +00000E
                    9-        CBEAM   507787  13      450028  450024  0.0     1.0     0.0             ++000003
                   10-        ++0000034       4                                                        +000004
                   11-        CBEAM   507788  13      450024  450025  0.0     1.0     0.0             ++000001
                   12-        ++0000014       4                                                        +000002
                   13-        CBEAM   507789  13      450031  450021  0.0     1.0     0.0             ++000007
                   14-        ++0000074       4                                                        +000008
                   15-        CBEAM   507792  13      450021  450032  0.0     1.0     0.0             ++000005
                   16-        ++0000054       4                                                        +000006
报错： *** USER FATAL MESSAGE 2040 (EMG)
     SINGULAR MATRIX FOR ELEMENT   507782
*** USER FATAL MESSAGE 2040 (EMG)
     SINGULAR MATRIX FOR ELEMENT   507783
*** USER FATAL MESSAGE 2040 (EMG)
     SINGULAR MATRIX FOR ELEMENT   507784
*** USER FATAL MESSAGE 2040 (EMG)
     SINGULAR MATRIX FOR ELEMENT   507785
*** USER FATAL MESSAGE 2040 (EMG)
     SINGULAR MATRIX FOR ELEMENT   507787
*** USER FATAL MESSAGE 2040 (EMG)
     SINGULAR MATRIX FOR ELEMENT   507788
*** USER FATAL MESSAGE 2040 (EMG)
     SINGULAR MATRIX FOR ELEMENT   507789
^^^ USER   WARNING MESSAGE 9031 (ERRPH1)   
^^^ NOGO ENCOUNTERED IN SUBDMAP    SEMG     
*** USER INFORMATION MESSAGE 4110 (OUTPX2)

如果换成不释放的话，就是在pins a  和pins b 位置输入0，而不是4的话就不报错了。
如下：
CBEAM     507794      13  450018  4500300.0     1.0     0.0            
CBEAM     507795      13  450018  4500290.0     1.0     0.0            
CBEAM     507796      13  450027  4500150.0     1.0     0.0            
CBEAM     507797      13  450015  4500260.0     1.0     0.0            
CBEAM     507798      13  450032  4500210.0     1.0     0.0            
CBEAM     507800      13  450021  4500310.0     1.0     0.0            
CBEAM     507801      13  450028  4500240.0     1.0     0.0            
CBEAM     507802      13  450024  4500250.0     1.0     0.0

请问我的问题出在哪里？实在想不明白，帮助文件也看了。
PA,PB in flags for beam ends A and B, respectively; used to remove
connections between the grid point and selected degrees-offreedom
of the beam. The degrees-of-freedom are defined in the
element’s coordinate system and the pin flags are applied at the
offset ends of the beam (see Figure 8-12). The beam must have
stiffness associated with the PA and PB degrees-of-freedom to be
released by the pin flags. For example, if PA = 4, the PBEAM
entry must have a nonzero value for J, the torsional stiffness. (Up
to five of the unique Integers 1 through 6 with no embedded
blanks.) Pin flags are not allowed for beam p-elements.
网上找了好多，都没找到关于梁单元CBEAM单元自由度释放的问题，想借版主大侠一些宝贵的时间帮忙看看，谢谢！

ywh362568

问题描述挺清楚的
除了两个mpc，其他约束怎么加的呢

wazon

感谢版主！是这样的，释放自由度和不释放自由度时其它约束及边界条件情况是一样的。。。自己尝试了一下，
1-        CBEAM   507782  13      450030  450018  0.0     1.0     0.0             ++00000B
2-        ++00000B123       456                                                        +00000C
或者
1-        CBEAM   507782  13      450030  450018  0.0     1.0     0.0             ++00000B
2-        ++00000B123       123                                                        +00000C
这样的计算都是可以的，不报错，尽管计算的结果跟我要的大相径庭
......还是搞不懂为啥换成4  4 就不对了

ywh362568

表面上看似乎有些奇怪，释放一个反而计算出错，2040错误，没见过
释放多了反而不出错
是过约束造成的矩阵奇异，出错？
你的rbe3定义了全部的自由度？

wazon

是的，rbe3单元定义了全部的自由度。。。我感觉这个可能跟rbe3没有关系，因为我自己找了一个小例子，让这个CBEAM连接壳单元也有这样的问题

*** USER FATAL MESSAGE 2040 (EMG)
     SINGULAR MATRIX FOR ELEMENT      219
^^^ USER   WARNING MESSAGE 9031 (ERRPH1)   
^^^ NOGO ENCOUNTERED IN SUBDMAP    SEMG

yema82811

这个BEAM元的4自由度是 壳元的6自由度，然后RBE3的主点的自由度最好只约束3个平动自由度（考虑主点的刚体运动），现在楼主约束了RBE3主点的所有自由度，释放了BEAM两边的4自由度，不知道这边是不是有问题。

yema82811

然后你用BEAM元直接连接壳元，只有一个beam元，没有技术处理，自由度不匹配，当然会有问题

wazon

版主和楼上兄弟你们好，谢谢你们的指点。我参考了一下资料，貌似不能同时释放梁单元两端的4自由度，只释放一边就可以达到释放自由度的目的。即改成：
1-        CBEAM   507782  13      450030  450018  0.0     1.0     0.0             ++00000B
2-        ++00000B4                                                                +00000C
是可以计算的，而且计算的结果也貌似可以。。。具体什么原因还是不得而知，呵呵，至少当前的问题算是解决了。

yema82811

昨天分析的太复杂了，其实一个很简单的问题，就是beam元通过RBE3连接壳元的变体而已。

而且也不是所说的不能同时释放单元两边的自由度。可以同时释放。不知楼主是要分析一个什么样的连接，和怎么样去判断计算的结果没有问题。