《第三章 杆梁结构的有限元分析原理》.ppt

第3章 杆系结构的有限元法;杆、梁单元概述;本章主要内容;3.1有限元分析的完整过程;1）用标准化的分段小单元来逼近原结构 2）寻找能够满足位移边界条件的许可位移场 3）基于位移场的最小势能原理来求解 基本变量为： ;完整的求解过程; 单元位移模式：u(x)=a0+a1x 单元节点条件：u(0)=u1， u(l)=u2 将式（b）代入式（a）， 从而得 ;回代得 写成矩阵形式为;根据几何方程可得单元应变的表达 单元应变写成矩阵形式为 简记为;根据物理方程可得单元应力的表达 单元应力写成矩阵形式为 简记为;;写成矩阵形式为;在得到各个单元的势能表达式后，需要进行离散单元的装配，以求出整个系统的总势能，对于该系统，总势能包括两个单元部分;处理边界条件是获取可能位移场，将左端的约束条件，即u1=0代入 上式可以得到简化的势能表达式;由于上式是基于许可位移场的表达的系统势能，这是由全部节点位 移分段所插值出的位移场为全场许位移场，且基本未知量为节点位 移，根据最小势能原理（即针对未知位移求一阶变分）有 ;将结构参数和外载荷代入上式有 求解得（单位m）;7）计算单元应变;8）计算单元应力;对于单元势能的表达，对其取极值有 具体地对于单元1，有 其中R1是节点1的支反力，P2是单元1的节点2所受的力，即单元2对该节 点的作用力，将前面求得的节点位移代入上式可得支反力大小。;以上是一个简单结构有限元方法求解得完整过程，对于复杂结构，其 求解过程完全相同，由于每一个步骤都具备标准化和规范性的特征， 所以可以在计算机上编程而自动实现。 讨论1：对于一个单元的势能取极值，所得到的方程为节点的位移和节 点力之间的关系，也称为单元的平衡关系，由此可以求出每一个单元 所受的节点力。;讨论2：由前面的步骤，我们也可以直接将各个单元的刚度矩阵按照节 点编号的对应位置来进行装配，即在未处理边界条件之前，先形成整 体刚度矩阵。 其物理意义是，表示在未处理边界条件前的基于节点描述的总体平衡 关系。在对该方程进行位移边界条件的处理后就可以求解，这样与先 处理边界条件再求系统势能的最小值所获得的方程完全相同。; 小结 有限元分析的基本步骤及表达式;有限元分析的基本步骤及表达式;一　拉压杆单元;　　用结点位移表示单元上任意截面的位移。对拉压杆单元，可以取其位移为一次多项式，即 　　　　　 由位移的边界条件： 可得系数 、 为：　　 这样，截面任意一点的位移 为： 　　　　　 用矩阵表示为： 　　　　　 其中 ; 根据材料力学中应变的定义，有 　　　　　 这里 为应变矩阵。由虎克定律，其应力为： 　　　　　 ; 利用虚位移原理求单元刚度矩阵，设杆端i、j分别产生虚位移 、 ，则由此引起的杆轴截面任意位置的虚位移为：　　　　　 对应的虚应变为： 　　 根据虚位移原理虚功方程（力乘以虚位移得虚功、外力虚功等于变形虚功），有： 　　　　　 将上式整理得： ;式中 ：为局部坐标系下单元结点荷载矩阵。 　 　　　记　 　　　　　 则可以得到拉压杆单元的单元刚度方程为： 　　　　　 这里 为局部坐标系下的单元刚度矩阵 ， 为局部坐标系下等效结点荷载矩阵，但值得指出的是：分布荷载 中可以包含集中荷载。根据定义，可以进一步求得单元刚度矩阵为： 　　　　　 ;二　扭转杆单元; 由材料力学可知，截面扭矩为： 　　　　 　 式中： 　　我们利用极小势能原理来进行单元分析，杆单元的势能用泛函表示为： 　　　　　 ; 由极小势能原理，取上述泛函的变分 ，可得： 　　 或者写为： 　 设： 　　　　 可得扭转杆单元的单元刚度方程为： 　　　　　 可以看到，其形式与拉压杆单元的单元刚度方程完全一致。同样，由上式可以进一步求得其局部坐标系下得单元刚度矩阵为： 　　　　　 ;三　只计弯曲的杆单元;;　　　　　 　　取挠曲线方程为 的三次多项式，即单元上任意一点的挠度为： 　　　　　 根据单元的位移边界条件： 　　　　时：