《有限元分析基础教程》（曾攀）笔记二-梁单元有限元方程推导

不得不说，Mathematica真是个好东西，以前学习有限元的时候，对于书中的方程推导，看到了就看过去了，从没有想过要自己推导一遍，原因是手工推导太复杂。有了MM，原来很复杂的东西突然变得简单了。

1.单元几何描述

上图是纯弯梁单元，长度l，弹模E，面积A，惯性矩I。两个节点1和2的位移列阵为

\[
q^{e}=[v_{1},\theta_{1},v_{2},\theta_{2}]^{T}
\]

$v$是挠度(defection)，或者叫位移；$\theta$是转角(slope)。需注意的是$v$和$\theta$的方向，一个是向上，一个是逆时针。

两个节点的节点力矩阵为

\[
P^{e}=[P_{v1},M_{1},P_{v2},M_{2}]^{T}
\]

当然实际情况往往是在梁的长度方向上作用有荷载，而不是只在节点处有，这时就要进行荷载等效，后面会有说明。注意这两个矩阵都是列矩阵。

需要注意的是，节点力矩阵表示的的是节点上的所有的力，不仅包括荷载引起的等效节点力，还包括节点的反力，反力矩等。

2.单元位移场表达

由于有4个位移节点的已知条件，那么假设纯弯曲梁单元的位移挠度函数具有四个待定系数，如下形式

\begin{equation}
v(x)=a_{0}+a_{1}x+a_{2}x^{2}+a_{3}x^{3}
\end{equation}
对于两端节点，位移和转角分别为$v_{1},\theta_{1},v_{2},\theta_{2}$，注意挠曲线方程在一点出的导