紧急求助：关于铰接点处理的问题

fuguoyang

我用fortran编写的空间刚架程序计算两个铰接点的门式刚架，对较接点我采用ansys中的方法，一个铰接点定义两个节点，然后将这两个节点的x、y、z平动自由度相等，而转动自由度分别独立，不知怎么回事算不出结果，通过调试我发现LDLT分解刚度距阵时D的一个元素为0，应该是约束不过导致刚度距阵的奇异。但我不知道怎么会出现这种情况？是我对铰接点的处理不对还是？请大家指点！

cae_ren

小弟觉得 铰接应该是只有一个旋转的自由度吧
应该耦合掉其他5个自由度吧  因为除了旋转的那个方向 另外两个旋转方向时无法相对转动的。

刚开始耦合  不确定思路正确与否  仅交流下