第六章-用有限单元法解解平面问题.ppt

§6-11 应用变分原理导出有限单元法基本方程 ·最小势能原理 一个体系在势能最小时, 系统方能处在稳定平衡状态. 最小势能原理： 总势能： EP是以位移函数d为自变量的函数, 故称为泛函(函数的函数); 一个泛函的最小值所对应的自变量函数, 就是泛函偏导数(变分)为0时, 所对应的自变量函数. §6-11 应用变分原理导出有限单元法基本方程 ·最小势能原理 根据最小势能原理: 可推导出和虚功原理相同的有限元基本方程. 例 题 深梁在跨中受集中力F的作用, 若取 μ = 0、t = 1试用有限单元法求解跨中的位移。 F 1m 1m 例 题 深梁在跨中受集中力F的作用, 若取 μ = 0、t = 1试用有限单元法求解跨中的位移。 解: ① 根据对称性此问题可取1/2分析 F F/2 例 题 深梁在跨中受集中力F的作用, 若取 μ = 0、t = 1试用有限单元法求解跨中的位移。 解: ② 对单元和结点进行编号 F/2 1 I II 2 3 4 I II 1 i 2 4 j m 1 4 3 i j m 单元号 I II 局部编码 整体编码 i 4 1 j 1 4 m 2 3 例 题 深梁在跨中受集中力F的作用, 若取 μ = 0、t = 1试用有限单元法求解跨中的位移。 解: ③ 叠加整体刚度矩阵 整体刚度矩阵(8×8)： 例 题 深梁在跨中受集中力F的作用, 若取 μ = 0、t = 1试用有限单元法求解跨中的位移。 解: ③ 叠加整体刚度矩阵 K11： kjj I + kii II i j m i j m kjj I kii II 例 题 深梁在跨中受集中力F的作用, 若取 μ = 0、t = 1试用有限单元法求解跨中的位移。 解: ③ 叠加整体刚度矩阵 K12： kjm I i j m i j m kjm I 例 题 深梁在跨中受集中力F的作用, 若取 μ = 0、t = 1试用有限单元法求解跨中的位移。 解: ③ 叠加整体刚度矩阵 K23： 例 题 深梁在跨中受集中力F的作用, 若取 μ = 0、t = 1试用有限单元法求解跨中的位移。 解: ③ 叠加整体刚度矩阵 整体刚度矩阵(8×8)： 例 题 深梁在跨中受集中力F的作用, 若取 μ = 0、t = 1试用有限单元法求解跨中的位移。 解: ④ 整体位移列阵 F/2 1 I II 2 3 4 例 题 深梁在跨中受集中力F的作用, 若取 μ = 0、t = 1试用有限单元法求解跨中的位移。 解: ⑤ 整体结点荷载列阵 F/2 1 I II 2 3 4 例 题 深梁在跨中受集中力F的作用, 若取 μ = 0、t = 1试用有限单元法求解跨中的位移。 解: ⑥ 求解线性方程组 例 题 深梁在跨中受集中力F的作用, 若取 μ = 0、t = 1试用有限单元法求解跨中的位移。 解: ⑥ 求解线性方程组 F/2 1 I II 2 3 4 显然位移 v3 v4 例 题 深梁在跨中受集中力F的作用, 若取 μ = 0、t = 1试用有限单元法求解跨中的位移。 解: 方法2 F/2 1 I II 2 3 4 大量元素被删除, 浪费计算量 例 题 深梁在跨中受集中力F的作用, 若取 μ = 0、t = 1试用有限单元法求解跨中的位移。 解: 方法2 F/2 1 I II 2 3 4 根据要计算的位移, 缩减整体刚度矩阵,可大幅降低计算量 只需计算3、4结点位移 求10个Krs 求3个Krs 例 题 深梁在跨中受集中力F的作用, 若取 μ = 0、t = 1试用有限单元法求解跨中的位移。 解: 方法2 F/2 1 I II 2 3 4 ⑥ 求解线性方程组 P142 6-5 6-6 本 章 作 业 * ·求解步骤 (1)划分单元网格, 对单元和结点编号 1 2 4 3 5 6 8 7 9 10 12 11 13 14 16 15 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 §6-8 解题的具体步骤 单元的划分 ·求解步骤 (2)选定直角坐标系, 按程序要求填写和输入有关信息 主要信息： 每个结点的坐标值; 每个单元的信息(单元内部结点的整体编码); 材料的弹性常数值; 荷载信息; 约束信息; §6-8 解题的具体步骤 单元的划分 ·求解步骤 (3)使用计算机程序进行计算 程序流程： 读入信息 形成整体刚度矩阵K 形成整体结点荷载列阵FL