博士论文序言
貌似在大论文里序言部分并不是必须的，然而我却以为这就像块自留地一样，可以相对自由的漫无目的说开去。先讲两个小故事。博士论文说到头其实就是求特征值实部的问题，事实上开始从这方面着手是相对较晚的。某日（在大论文已经写作很久之时）翻看本科时常微分方程教材时，说实话有些吃惊。这门课是本科生涯中考得最高分的一门课，然而当时教学范围只有前三章，而全书共有四章。那么最后一章是什么呢？其大标题是“稳定性与……
这里只讲杆单元刚度矩阵的推导，基于前提是已经根据势能最小原理得到了局部坐标系下单元刚度矩阵表达式，即对于弹性杆单元，通用单元刚度矩阵表达式为……
这里只讲杆单元刚度矩阵的推导，基于前提是已经根据势能最小原理得到了局部坐标系下单元刚度矩阵表达式，即对于弹性杆单元，通用单元刚度矩阵表达式为：
——沿杆单元方向建立局部坐标，起点为0，端点为杆单元长度
——沿杆单元方向位移
对于其他符号意义可参见于另一文：有限元推导基本符号
参考书籍：辛科维奇《有限元方法》第五版曾攀等译
采用两节点拉格朗日多项式形式的插值函数（关于插值函数可参见于《有限元方法》第五版第八章）：
即这里采用的这里形函数为：
由于形状规则，对于上述坐标选取下积分依旧方便计算，故理论上可以不采用单元映射方法（更多可见于可参见于《有限元方法》第五版第九章）。不过这里通用起见，采用等参单元，即另取局部坐标，满足：
相应形函数为：
雅可比矩阵为：
这里： 则： 故： 由于积分容易得到，故不采用高斯积分方法而直接求解定积分： 如果扩展到y向，则有：
存在转换矩阵
(这里符号与前述雅可比矩阵相同，但矩阵并不相同)使得：
故： 则： 相应有整体坐标系下单元刚度矩阵为：
每一个单元通过节点转换矩阵集成入整体刚度矩阵中去，即： 相应的如果有节点荷载向量为，则得到最终有限元方程为： 上述杆单元刚度矩阵的推导可以认为是最简单的了。对于实际的复杂单元，上述求解单元刚度矩阵过程中需要进行积分。由于一般情形下无法求出不定积分，故通常采用高斯积分形式。
对于梁单元亦可采用上述方法进行推导，只不过形函数、应变微分算子、物理方程等不相同。
上问推导可能存在的一个困惑是怎么会出现两个局部坐标系。其原因在于本文推导时为讲究方法通用性应用了等参变换，实际上对于本例而言是没有必要的，不进行等参变换依旧能得到相同结果。由于采用了等参变换，第一个局部坐标就相当于第二个局部坐标的整体坐标了。