[物理]有限元法基础及Ansys应用10.pptVIP

Event Name Here BEA Confidential 教师：李鸿秋 lihongqiu@jit.edu.cn 第二章 平面问题有限元法 第3章 杆单元和梁单元 3.1 杆件系统的有限元分析方法 3.1 杆件系统的有限元分析方法 3.1 杆件系统的有限元分析方法 3.1 杆件系统的有限元分析方法 3.1 杆件系统的有限元分析方法 3.1 杆件系统的有限元分析方法 3.1 杆件系统的有限元分析方法 3.1 杆件系统的有限元分析方法 3.1 杆件系统的有限元分析方法 3.1 杆件系统的有限元分析方法 3.1 杆件系统的有限元分析方法 3.1 杆件系统的有限元分析方法 3.1 杆件系统的有限元分析方法 单元应变场 　　根据单元的位移场函数式，由几何方程可以得到单元的应变场表达式， 记为： 　　这里，B矩阵称为几何矩阵。B矩阵可以表示为分块矩阵的形式 其中 注意：矩形单元的应变场为一次线性函数。 单元应力场 　　由物理方程及应变矩阵，可以得到单元的应力场表达式， 其中　　　 为应力矩阵，D称为弹性矩阵，对于平面应力问题， 注意：矩形单元的应力场为一次线性函数。 将应力矩阵表示为分块矩阵的形式 其中： 对于平面应变问题，只需将E换为　　　，　换为　　　。 单元刚度矩阵 　　和三角形单元一样，可以根据虚功理导出节点位移向量和节点力向量之间关系，即单元的刚度矩阵　　，可以将其写成分块的形式。 其中 　　对于平面应力问题，如果单元厚度t为常数，则可得刚度矩阵的显式形式： 积分得： 单元刚度矩阵 单元的体积力和表面力引起的节点力和三角形三节点单元的计算公式相同。由于位移分量在x为常数及y为常数的直线上是线性变化的，因此节点的分配也符合静力等效原则。 （1）对于单元的自重W，载荷列阵为： 单元等效节点力 （2）如果单元在一个边界上受到三角形分布的表面力，在该边界上一个节点处为0，另一个节点处为最大，在将总表面力的1/3移植到前一个节点，2/3移植到后一个节点。 矩形单元和三角形单元的比较：矩形单元中的应力分量不是常量，若弹性体中采用相同数目的节点时，矩形单元比比三角形单元能更好的反映应力急剧变化的情况。 但是，矩形单元比三角形单元的适应性差，不能适应斜线及曲线边界，且不便于对不同部位采用大小不等的单元。 单元等效节点力 整体平衡方程：根据单元的刚度矩阵K，等效节点力矩阵，按对号入座的方式叠加组装，整体刚度矩阵和节点载荷矩阵，从而得到整体平衡方程。 八、计算结果处理 ?1=(??+??)/2 ?2=(??+??+??+??+??+??)/6 ?1=(??+??)/2 ?2=(??+??)/2 连续体有限元分析中的几次近似 1. 逼近性离散（网格剖分） 提高精度的办法： 减小单元尺寸 改变单元形状 本章主要介绍利用杆单元及梁单元进行结构静力学的有限元分析原理。首先介绍了杆单元的分析方法，详细给出了采用杆单元进行有限元分析的整个过程,以加深对有限元法的理解。 杆件只承受轴向力，可以视为一种特殊的梁单元，本节将采用有限元法来分析杆件系统，以下给出规范的有限元法中关于杆单元的推导过程，以及整个杆系的求解过程。 如图3-1所示的杆件结构，左端铰支，右端作用一个集中力，相关参数如图。具体求解过程如下： 图 3-1 杆件结构 （1）确定坐标系、单元离散，确定位移变量, 外载荷及边界条件。 要建立两种坐标系：单元坐标系（局部坐标系）、整体坐标系。根据自然离散, 坐标系建立成一维, 单元划分为两个, 给出相应的节点1、2、3以及相应的坐标值（见图3-1）。在局部坐标系中，取杆单元的左端点为坐标原点，图3-2为任取的一个杆单元。 图 3-2 杆单元 对于两个节点的杆单元，存在如下节点力和节点位移的关系式 (3.1) 其中， 称为单元刚度矩阵 （2）确定位移模式 假设单元位移场： 取其线性部分，系数 、 可由节点位移 、 确定 (3.2) （3）形函数矩阵的推导 由单元的节点条件, 两个节点坐标为x1、x2,两个节点位移为 ， ，代入上式插值模式公式得： 求解得到 这样， 可以写成如下矩阵形式 导出 (3.3) 得到形函数矩阵 (3.4) 记节点位移矢量是 (3.5) 因此，用形函数矩阵表达的单元内任一点的