稳定与屈曲
一直觉得这些写稳定的书籍有一个事情没有讲清楚，就是什么是稳定。一上来就讲几种失稳（稳定分叉、不稳定分叉等）会让人觉得有些不知所措。这里就谈谈自己的理解。……
这里只讲杆单元刚度矩阵的推导，基于前提是已经根据势能最小原理得到了局部坐标系下单元刚度矩阵表达式，即对于弹性杆单元，通用单元刚度矩阵表达式为……
这里只讲杆单元刚度矩阵的推导，基于前提是已经根据势能最小原理得到了局部坐标系下单元刚度矩阵表达式，即对于弹性杆单元，通用单元刚度矩阵表达式为：
——沿杆单元方向建立局部坐标，起点为0，端点为杆单元长度
——沿杆单元方向位移
对于其他符号意义可参见于另一文：有限元推导基本符号
参考书籍：辛科维奇《有限元方法》第五版曾攀等译
采用两节点拉格朗日多项式形式的插值函数（关于插值函数可参见于《有限元方法》第五版第八章）：
即这里采用的这里形函数为：
由于形状规则，对于上述坐标选取下积分依旧方便计算，故理论上可以不采用单元映射方法（更多可见于可参见于《有限元方法》第五版第九章）。不过这里通用起见，采用等参单元，即另取局部坐标，满足：
相应形函数为：
雅可比矩阵为：
这里： 则： 故： 由于积分容易得到，故不采用高斯积分方法而直接求解定积分： 如果扩展到y向，则有：
存在转换矩阵
(这里符号与前述雅可比矩阵相同，但矩阵并不相同)使得：
故： 则： 相应有整体坐标系下单元刚度矩阵为：
每一个单元通过节点转换矩阵集成入整体刚度矩阵中去，即： 相应的如果有节点荷载向量为，则得到最终有限元方程为： 上述杆单元刚度矩阵的推导可以认为是最简单的了。对于实际的复杂单元，上述求解单元刚度矩阵过程中需要进行积分。由于一般情形下无法求出不定积分，故通常采用高斯积分形式。
对于梁单元亦可采用上述方法进行推导，只不过形函数、应变微分算子、物理方程等不相同。
上问推导可能存在的一个困惑是怎么会出现两个局部坐标系。其原因在于本文推导时为讲究方法通用性应用了等参变换，实际上对于本例而言是没有必要的，不进行等参变换依旧能得到相同结果。由于采用了等参变换，第一个局部坐标就相当于第二个局部坐标的整体坐标了。