梁单元的有限元实现——有限元实践笔记（3）

讲完了弹簧单元和杆单元，接下来讲一讲稍微复杂一点的梁单元的有限元实现。

什么是梁单元

梁单元为细长结构构件，通常只承受横向载荷。无论是杆单元还是梁单元都是结构力学中基本的承力构件之一，二者都是一维单元，物理特征基本相同，二者之间的区别主要是承力形式不同，规定杆单元只承受轴向力，梁单元为可承受横向力和弯曲力矩。下图中规定了梁单元节点位移、转角、力、力矩的正方向，力和位移的正方向指向y轴的正向，转角和力矩为逆时针为正。

引言

采用直接刚度法推导线弹性梁单元的刚度矩阵，应用有限元分析的一般步骤实现梁单元的结构的静力学问题求解。
本文由以下几个部分组成：

• Euler-Bernoulli 梁理论简介
• 局部坐标系下梁单元的刚度矩阵和方程的推导
• 刚度矩阵和方程的组装
• 坐标变换
• 平面刚架结构的有限元分析实例及matlab代码

1 Euler-Bernoulli 梁理论简介

Euler-Bernoulli 梁理论是由Euler和Bernoulli两个人推导得到的，并由此命名。Euler-Bernoulli 梁理论是基于平断面假定建立的，即梁弯曲前后，垂直于纵向主轴的平面保持为平面并垂直于主轴，只考虑弯曲变形，该理论适用于线弹性梁的力学行为分析。
![在这里插入图片描述](https://img-blog.csdnimg.cn/20210124153423513.png?x-oss-process=image/watermark,type_ZmF
梁单元的控制方程为：

E为弹性模量，I为截面惯性矩，w(x)为作用在梁上的分布力，v为梁的横向位移，即挠度，x为纵向坐标。

2 局部坐标系下梁单元的刚度矩阵和方程的推导

有限元分析的一般步骤如下：
step1 单元类型的选择
确定单元自由度，本部分仅对承受横向载荷的梁单元进行讨论，因此每个单元节点有两个自由度，单元节点位移为$v_1$和 ϕ 1 \phi_1