有限元（梁系）.ppt

第三章 梁系结构的有限元法 梁的基本概念 梁--承受弯曲变形的杆形构件； 梁模型--平面假设将梁简化成一维问题，可由截面上两中性轴交点的连线（梁轴线）代表； 梁的变形特点—既与挠度关联，也与弯角关联；（分别展示四种基本情况）； 挠度与弯角的关联性：弯角为挠度沿轴线方向的一阶导数（斜率，小变形）。 梁单元的基本概念 节点位移--平面梁：挠度和弯角，2自由度；广义平面梁：增加一个轴向位移（拉压），3自由度；空间梁：两向挠度和两向转角，4自由度；广义空间梁（略）； 节点力—按做功方式与节点位移一一对应，如：挠度对应剪力，弯角对应弯矩，等等。 梁结构有限元法的基本步骤 结构离散--梁结构离散为若干个单元，之间以节点相连，节点为关联单元共有（连续性）； 单元分析--建立单元节点力与节点位移的关系；包括：单元位移插值、应变与应力分析，节点力合成； 节点平衡--对于每一节点，建立所有相关节点力与外力的平衡关系；得到一组线性方程； 约束处理—引入已知节点位移（约束条件），使上述方程组成为可解方程组； 求解节点位移，进一步可计算单元应变和应力以及约束节点力。 梁单元分析 一、平面梁单元的位移函数 1.梁单元的局部坐标系 设单元的位移函数为： 二、建立节点位移与节点力关系 1、 轴向节点力 2、弯矩节点力 3、剪切节点力 由梁微元的平衡关系（简单推导）有： 约束处理 一、消行消列法 把总刚度方程进行换行换列（原方程的解仍保持不变）后重新进行排列。注意：换行的同时必须换列，才能使变换后的方程组与原方程组完全等数 这样就可以把总刚度方程分为两组。其中： ----已知结点荷载列向量 ----未知结点荷载列向量 ----未知结点位移列向量 ----已知结点位移列向量 展开（3-13）式，我们可得到如下两个方程。 为简单起见，令 （固定约束） 则由式（3-15a） 得： 再把（3-16）代入式（3-15b）可直接求得 因此，只要引入已知的边界约束条件 ，未知的结点位移 及未知的支反力 都是可求的。 二、置大数法（乘大数法） 设结构的约束条件为： 结构的总刚方程为： 设： 为已知位移 ，求解步骤为： 1．修改[K]：将主元素 → 。其中N为充分大数，使得 一般取N= ；N太大易溢出，N太小精度不高。 2．修改{F}：将 → ，经过这样的处理后，总刚度矩阵的第r个方程为： 由于N ，故可得： 对S个约束条件进行同样的处理。经过上述修改后，整个刚度方程完全满足了结构在支座处的位移边界条件，如果边界位移条件足以消除结构的刚性位移，则修改后的刚度矩阵已是非奇异的对称矩阵，通过求解结构刚度方程，即可求得未知的结点位移。 三、划零置一法 仍设约束 ， 为已知的支承处的位移。 结构的总刚度方程为： 按下面的做法修改总刚度方程： 1. 修改[K]：令 2. 修改{F}：令 （ 为给定数据） 修改后的第r个方程为： 满足给定边界条件 修改后的第i个方程为： 即： 仍然满足原来的第i个方程。 注：消行消列法一般用于手算。 置大数法与划零置一法用于电算。不同的是置大数法在约束处理上是近似方法；而后者在约束处理上为精确法。代价就是置大数法处理比置一法更方便、简单。 载荷处理（load process） 杆系结构的有限元方法就其理论体系本身，并不涉及跨间载荷，而假设结构所作用的载荷均为结点载荷。对于结构上作用的非结点载荷，对其进行某种变换（conversion），把其转换为等效的结点载荷。 就载荷等效变换的方法而言，在杆系结构中，一种是用虚功原理（比较严谨），另一种是用位移法的概念来进行，比较简单。 位移法的处理方法：按结构力学位移法解题思想求该单元的固端力和固端力矩，将这些固端力和固端力矩反方向作业于结构的对应结点，即该单元非结点的等效结点载荷。 * * l (e) i j 有： 其中 为待定系数 确定待定系数： l (e) i j 把 待定系数代入位移插值函数，有： 容易验证 ： (3-1a)，(3-1b)或(3-2a)，(3-2b)称为平面梁单元的位移插值函数 在局部坐标中： 由梁截面正应力合成（简单推导）可得 ： 对节点i(X=0)有： 对节点j(X=L)有： 在节点i有： 三、梁单元的单元刚度矩阵 1、单元刚度矩阵 用