第三章 平面线弹性问题有限元.pptVIP

第三章 平面线弹性问题的有限元 ②单元刚度矩阵得物理意义 由虚功原理得节点力与节点位移的关系 根据矩阵乘法可得 由此可知：单元刚度矩阵的2×2阶子矩阵[Krs]表示使节点S（s＝i，j，m）产生单位位移时，在节点r（r＝ i，j，m ）上所需要施加的节点力得大小。 第三章 平面线弹性问题的有限元 如果将Krs展开 就是使节点s产生一个水平单位的位移，在r节点上所需要施加得垂直力 r s 第三章 平面线弹性问题的有限元 ③单元刚度矩阵的特性 a. 对称性 [K]e是对称矩阵，即Kij＝Kji，可以用功的互等原理来证明 Fr 由Fr引起的位移Usr Fs 由Fs引起的位移Urs 对弹性体有Fr·Urs＝ Fs·Usr 对于单元e，在节点i，j上 FjxUj＝FiyVi 假定Uj＝Vi＝1 则Fjx＝Kji Fiy＝Kij ∴Kij＝Kji i j Fiy Vi Uj Fjx b. 奇异性 不存在 是由于没有引入位移边界条件的结果。 c. 分块性 第三章 平面线弹性问题的有限元 §3.3 总刚矩阵的形成及其修正 1. 利用节点平衡组成总刚矩阵 外力＝节点力→未知的内力 节点荷载＝节点力 （单刚矩阵描述了节点力与节点位移的关系） 第三章 平面线弹性问题的有限元 第三章 平面线弹性问题的有限元 根据节点受力分析可得 R1x，R1y，R2y为约束反力，R2x＝R3x＝R4x＝0 第三章 平面线弹性问题的有限元 以上可以将它们按节点分单元写成矩阵形式，即 第三章 平面线弹性问题的有限元 写成分块矩阵的形式为： 那么对于节点力与节点位移的关系可表示为： 第三章 平面线弹性问题的有限元 将单元刚度矩阵、节点位移、代入节点力与节点荷载的平衡方程得： 第三章 平面线弹性问题的有限元 一般地 如果一个单元组合有n个节点，m个单元，那么总体刚度矩阵就是m个单刚矩阵由6×6阶升阶2n×2n阶以后叠加的总和。 第三章 平面线弹性问题的有限元 2. 总刚矩阵的形成规律 ①[Krs]当r＝s时，就是共用节点的所有单元的单刚矩阵子块的叠加结果。 ② [Krs]当r≠s时，若rs是结构体的内边，就是共用rs这条边的单元的对应的子块的叠加，如53边，若rs是外边，就是使用这个单元的对应子块，如14边。 ③若r、s不同属于一个单元时，则：Krs＝0。 总刚矩阵[Krs]中的两个脚标，r表示节点力作用的节点编号，s表示产生位移的节点编号，当r＝s时，该节点位移与所有共用单元在该节点的节点力有关，当r≠s时，节点s的位移就与rs所在单元的节点r的节点力有关，当r、s不属于同一单元，则节点s的位移与节点r的节点力没有直接关系。 第三章 平面线弹性问题的有限元 3. 总刚矩阵的特征 a. 对称性 Krs＝Ksr c. 奇异性 不存在 b. 稀疏性 证明 第三章 平面线弹性问题的有限元 4. 位移边界条件的引入及总刚矩阵的修正 ①零位移边界条件，删除对应的列和行 ②非零位移边界条件 两种方法： a. 一般方法 i）令已知位移u0i对应的荷载：Ri＝u0i； ii）修改总刚矩阵中相应行、列元素，使： Kij＝Kji＝0， Kii＝1； iii）修改荷载矢量{R}，令： b. 主对角线元素优先法 i）将与已知节点位移ui＝u0对应的总刚矩阵中的主对角线元素Kii乘以一个绝对大的数，如A＝1012； ii）将对应荷载项改为AKiiu0； iii）其余各项保持不变。 岩土工程地质问题分析中常用的方法 也是一种常用方法 第三章 平面线弹性问题的有限元 §3.4 总刚矩阵存贮 1. 总刚矩阵的存贮方法 ①根据对称性，只存贮其下三角部分 ②根据稀疏性，只存贮非零元素 为了便于研究存贮方法，计算存贮地址 我们引入两个概念 半带宽：每一行第一个非零元素到对角线元素之间元素个数。 行宽：不包括对角线元素的带宽。 按不同行的不同带宽的存贮方法，我们称为变带宽存贮。 将二维数组变为一维数组存贮。 第三章 平面线弹性问题的有限元 2. 变带宽一维数组存贮的实现 用一维数组AK（J）来存贮总刚矩阵的下三角部分的元素。 设每一行对角线在一维数组的地址为D2 则任一元素KD1D3在一维数组的地址为： J＝D2－（D1－D3）＝D2－D1＋D3 第三章 平面线弹性问题的有限元 现在要计算主对角线元素的地址D2： 某行主对角线存贮地址＝该行行宽＋上行主对角线元素存贮地址＋1 LA（I）＝LA（I）＋LA（I－1）＋1 某行主对角线元素存贮地址 该行行宽 上行主对角线元素存贮地址 下面计算行宽 第三章 平面线弹性问题的有限元 第D1行的行宽： LA（D1）＝（D1－D3）max 行号 列号 D10，D30 D1≥D3 （只考虑下三角部分