节点力的求取及debug命令讨论帖

zhuweihong2004

为了求取某节点在载荷作用下的内力，若用组装后的整体刚度矩阵则比较麻烦，因为通过K*U=F反推节点内力前提就是得列出所有节点的各个方向的唯一，然后代入公式中求得F,这对集中载荷的静态分析是可行的，但是当载荷为动载荷，特别是在进行谐波分析时，得到得往往是很多组关于频率变化的位移，那么再回代方程的话就比较麻烦，个人采用的一下方法,不知道合适与否，仅分享经验，若有不同意见，大家可以讨论：
1，首先用debug命令获得单元的刚度，质量及阻尼矩阵，
2，根据所求节点信息，找出所有共享该节点的单元，
3，组装局部单元的刚度矩阵，组装局部的单元的所有节点的位移，这么做的依据是在模型和刚度矩阵的组装过程中，不相关的两个节点，彼此的单元在刚度矩阵中是没有贡献的，这么一来，刚度矩阵的维数就大大的降阶了，求解也比较简单，
4，将所有数据代入公式，很方便的可求得节点力
    本人采用的是幅值相位的方式加载的，对于按实虚部加载的方式没有验证，同时，如果有阻尼力存在的话，按照广义刚度矩阵的定义，则必须按实虚部方式加载。
最后还有几个待解决的问题：怎么定义某个面上的力（不是单元应力）和速度，相位的因素在上述过程中如何考虑
    以上是本人新近学习的几点心得，有什么错误的地方还请指教，谢谢