计力第2章-2.pptVIP

有限元方程： 1. 每一个方程式表示: 一个结点在一个自由度方向的平衡条件， 上式是弹性体所有结点平衡方程的汇集。 2. 所选择的位移函数保证了变形连续性， 协调条件自然满足。 3. 已考虑了外载荷势能，故方程包含外力边界条件 4. 方程尚不满足几何边界条件，目前还不能解出全部结点位移。 如： 其中： 为给定的结点方向位移。 2. 引入位移边界条件的方法 （1） 直接代入法 将向量 a、P 分成两部分： —— 待定的结点位移； —— 已知的结点位移。 （2.2.32） 将刚度阵 K 分成4个分块矩阵： 方程（2.2.32）变为： （2.2.52） 由 K 的对称性，显然有： 将方程（2.2.52）展开： （2.2.53） 令： （2.2.54） —— 给出 2n－l个方程，求解 2n－l个结点位移。 （1） 直接代入法 （2.2.54） （2.2.55） 讨论： （1）需重新组合方程，新的方程阶数降低了。 （2）重新组合方程时，已将原方程的顺序改变，不易编程实现。 —— 给出 2n－l个方程，求解 2n－l个结点位移。 缺点：矩阵要重新排列，引起不必要的麻烦。 （2） 对角线元素改“1”法 设： （如：无移动的铰支座、链杆支座等。） 当： 时， 优点： （1）不改变原方程的阶数； （2）不改变结点未知量的顺序； —— 给程序编制带来方便。 缺点： 不能用于非零约束 （3） 对角线元素乘大数法 设： 则将 其中： 第 j 个方程 修改后的 第 j 个方程成为： 第 j 个方程近似成为： 优点： （1）方程的阶数不变； （2）结点未知量的顺序不变； —— 给程序编制带来方便。 （3）不易引起修改后的方程成为病态方程，提高计算精度与稳定性。 对角线元素乘大数法 —— 为现有有限元程序中常用的方法。 引入位移边界条件后，得到的有限元方程表示成： 上述方程中，刚度矩阵 不再是奇异矩阵， 可以具体求解。 2.2.8 应用举例 * * * 中国矿业大学 王培荣 2.2.5 结构刚度矩阵 K 和 结构等效结点载荷列阵 P 的集成 总刚度阵 整体等效结点载荷列阵 1. 刚度矩阵 K 的集成 设系统有 n 个结点，n e 个单元： 式中： （2.2.50） 的作用是： 将每个单元刚度矩阵的子块，按该单元结点的总体编号送入总刚中合适的位置。 这一公式是为了表示方便，实际程序中，只须将单刚根据单元结点编号送入总刚即可。 3 10 9 e i j m 说明： （1） —— 只有9个子矩阵不为零，其余子矩阵全为零。 （2） —— 即将 扩大到与总刚度阵 中的各子矩阵按照 （3） —— 即将 单元结点的实际编码放入总刚度阵中。 反应单元刚度阵 Ke 对总刚度阵的贡献。 K 同阶矩阵，以便矩阵相加。 2. 结构等效结点载荷列阵 P 的集成 （2.2.51） 3 10 9 e i j m 对每一个单元 的作用是将单元载荷 列阵作用在每个单元结点上的等效载荷， 送到总载荷列阵的相应位置。 例： 总刚度矩阵形成示例 1 2 3 4 5 6 ① ② ③ ④ 单元号： 对应的整体节点号： ① ② ③ ④ （1）各单元刚度阵： （2）将单元刚度阵扩大成与总刚度阵相同的阶数 1 2 3 1 2 3 1 2 3 4 5 6 1 2 3 4 5 6 1 2 3 4 5 6 ① ② ③ ④ 1 2 3 4 5 6 1 2 3 4 5 6 2 5 3 2 5 3 1 2 3 4 5 6 ① ② ③ ④ 5 6 3 5 6 3 1 2 3 4 5 6 1 2 3 4 5 6 1 2 3 4 5 6 ① ② ③ ④ 2 4 5 2 4 5 1 2 3 4 5 6 1 2 3 4 5 6 1 2 3 4 5 6 ① ② ③ ④ 1 2 3 4 5 6 ① ② ③ ④ 1 2 3 4 5 6 1 2 3 4 5 6 2.2.6 结构总刚度矩阵 K 的特点 1. K 中元素 kij 的物理意义 当结构的第 j 个结点位移方向上发生单位位移，而其它结点位移方向上位移为零时，需在第 i 个结点位移方向上施加的结点力大小。 2. K 的性质（特点） （1）对称性 （2）奇异性 每行（或每列）的所元素之和等于零， K 的