求助关于裂纹梁的裂纹刚度矩阵如何添加

mihuaguai

用matlab编写裂纹梁的有限元程序。
需要将裂纹刚度矩阵添加到总体刚度矩阵当中，看了很多论文，但自己在实现的时候都出问题。

很多论文里都说“采用线弹簧模型，裂纹单元无长度，所以不考虑裂纹对总质量矩阵的共享。裂纹的存在对两侧节点横向位移没有影响，只考虑转角即可。裂纹单元刚度矩阵为[Kt -Kt;-Kt Kt]。获得了裂纹等效刚度矩阵后，可以很方便的加到总刚度阵中。”

可是我就不明白这个很方便是怎么加的。试了好几种，出来的结果感觉也不对。
不知道谁做过裂纹方面的研究，求达人指教。

另外，我参考的论文有 余志刚《基于高阶有限元方法的裂纹斜梁振动特性分析》
                                    陈雪峰 《基于小波有限元的悬臂梁裂纹识别》
                                    S. Chinchalkar 《determination of crack location in beams using natural frequencies》


谢谢。