关于梁挠曲线方程解析解与有限元解的简单对比及其他

一、对比

对于受均布荷载的简支梁，假设梁的长度L=1，均布荷载大小为1，弹模E为1，惯性矩I为1，那么梁的挠曲线方程的解析解为

$$
v(x)=-\left(\frac{x^4}{24}-\frac{x^3}{12}+\frac{x}{24}\right)
$$

梁单元的形函数如下

根据形函数的定义，只要知道了单元两端节点的位移，即可用形函数求出单元上任意位置的挠度。即

$$
v(x)=N(x)\cdot q^{e}
$$

其中$q^{e}$为节点的位移列阵。

下面对比一下，由形函数得到的挠度曲线与挠曲线的解析解之间的差别，为了对比方便，这里采用形函数时，节点的位移均由挠曲线解析解求出，换言之&