杆与梁结构的有限元法.ppt

§2.3.1 简单梁单元 再设： 同理，由梁的变形公式和平衡条件可求得刚度矩阵的第二列元素： 梁单元变形 由刚度方程可得： §2.3.1 简单梁单元 同样的方法可以求出其余2列元素，从而求出单元刚度矩阵： 显然，与弹簧和杆单元一样，该梁单元的刚度矩阵具有如下性质： 1）对称性； 2）奇异性； 3）主对角元素恒正。 刚度矩阵求得后，单元特性就完全确定： §2.3.1 简单梁单元 采用矩阵分块方法和运算规则，对梁单元的刚度方程按节点进行分块。 单元节点力列阵分块： 单元节点位移列阵分块： 分块形式的单元刚度矩阵： 上面每一子块均为2×1子列阵。 每一子块均为2×2子矩阵 3、单元刚度方程的分块 §2.3.1 简单梁单元 将上式按分块矩阵乘法展开，得两个矢量方程（共4个代数方程）： 因此，单元刚度方程分块形式表示为： 从上面方程可以看出梁单元刚度矩阵子块的物理意义：相关节点位移对对应节点力的贡献。 上面按分块形式表示的单元刚度方程——节点力～节点位移关系在整体分析中集成单元特性时更加简洁，在有限元分析中广泛采用。 §2.3.1 简单梁单元 §2.3.1 简单梁单元 三、离散结构的整体分析 设已知分块形式的各单元特性方程： §2.3.1 简单梁单元 以离散结构的各节点作为隔离体，以节点2为例，建立其平衡方程。 单元节点力的反作用力 外载荷 单元节点力 单元节点力 节点2的受力分为两类： 1）外载荷： 2）单元（1）、（2）上节点力的反作用力： §2.3.1 简单梁单元 由节点2的静力平衡条件得： 单元节点力的反作用力 外载荷 单元节点力 单元节点力 节点2的外载荷=节点2对其所有相连单元的节点力之和（节点总内力） 也就是节点2所受外载荷 要分配到相连的单元上。 §2.3.1 简单梁单元 由前面给出的单元（1）、（2）分块形式单元刚度方程代入节点2的平衡方程： §2.3.1 简单梁单元 同理，由节点3的平衡可得： 由节点1、4的平衡得： 将上面4个节点的平衡方程合并，写成矩阵形式得： §2.3.1 简单梁单元 上式简写为： —— 结构节点位移列阵（ 8 ×1） —— 结构节点载荷列阵（ 8 ×1） —— 结构总刚度矩阵（8 ×8） —— 结构（系统）有限元平衡方程 §2.3.1 简单梁单元 结构总刚度矩阵也可以由各单元刚度矩阵扩大到整体规模后叠加而成，方法同前面的弹簧单元和杆单元。 由于单元刚度矩阵在扩大和叠加过程中，其具有的性质（对称、奇异、主对角元恒正）不变，因此结构总刚度矩阵仍然保持这些性质。 总刚度矩阵中有大量元素为0，因此矩阵具有稀疏性 非零元素沿主对角线呈带状分布（节点编号满足一定条件）。 结构总刚度矩阵的讨论： §2.3.1 简单梁单元 总之，从弹簧、直杆和梁结构有限元总刚度矩阵的特点可以归纳出结构有限元总刚度矩阵的性质如下： 1）对称性； 2）奇异性； 3）稀疏性； 4）非零元素带状分布 §2.3.1 简单梁单元 结构有限元平衡方程的讨论： 平衡方程左边总刚度矩阵与位移列阵之积等于结构中各节点的总节点力（各节点对相关单元作用力之叠加）；因此，总刚每行各子块表征相应节点位移对该行对应总节点力的贡献——总刚子块的物理意义。 §2.3.1 简单梁单元 平衡方程右端是各节点外载荷，左端是由节点位移和单元刚度矩阵子块叠加计算得到的总节点力。因此，有限元平衡方程表征了系统各节点所受外载荷与所受所有相关单元反作用总力（总节点力）之间的平衡。 结构有限元平衡方程可以叙述为： 总节点力（内力） = 节点外载荷。 §2.3.1 简单梁单元 对于特定结构，方程中必存在已知位移和相应的未知载荷（支反力），因此，平衡方程求解前必须进行约束处理，分离出关于未知位移的方程进行求解。然后再用求出的位移，通过剩余方程求出支反力。 §2.3.2 平面内一般梁单元 拉伸、弯曲组合 单元变形特征 节点位移分量 节点载荷分量 平面梁单元 整体节点位移 整体节点载荷 节点自由度：3 一、单元与节点 平面刚架 模拟 单元有2个节点：i，j 局部坐标系下节点位移分量： 轴向位移： 横向挠度： 转角： 局部坐标系下节点力分量： 轴向力： 横向剪力： 弯矩： §2.3.2 平面内一般梁单元 单元描述： 二、局部坐标系下平面梁单元刚度方程 单元有6个位移分量 6个自由度 单元节点位移列阵： 单元节点力列阵： §2.3.2 平面内一般梁单元 单元描述：