平面弹性问题位移—应力混合重心插值配点法

作     者：张磊 

作者单位：山东建筑大学 

学位级别：硕士

导师姓名：王兆清

授予年度：2017年

学科分类：07[理学] 08[工学] 070102[理学-计算数学] 0701[理学-数学] 080102[工学-固体力学] 0801[工学-力学（可授工学、理学学位）] 

主      题：弹性力学问题 重心Lagrange插值 位移-应力混合公式 配点法 微分矩阵 

摘      要：弹性力学问题可归结为二阶耦合椭圆形偏微分方程边值问题。工程中遇到的大部分问题都难以得到其解析解。为求解弹性力学方程,工程实际中广泛采用数值求解技术。本文提出数值分析平面弹性问题的位移-应力混合重心插值配点法。将弹性力学控制方程表达为位移和应力的耦合偏微分方程组,采用重心插值近似未知量,利用重心插值微分矩阵得到平面问题控制方程的矩阵形式离散表达式。使用重心插值离散位移和应力边界条件,采用附加法施加边界条件,得到求解平面弹性问题的过约束线性代数方程组,应用最小二乘法求解过约束方程组,得到平面弹性问题位移和应力数值解。对于不规则区域的弹性力学问题,采用重心Lagrange插值正则区域法,将不规则区域嵌入规则区域,在规则区域上采用重心Lagrange插值近似未知函数。利用配点法强迫微分方程在离散节点处精确成立,得到规则区域位移-应力混合方程组。在不规则区域的边界上取若干节点,由规则区域内的重心插值插值节点的未知函数,得到一个边界条件的约束代数方程。将位移-应力混合方程的离散方程和边界条件的约束方程组合成一个新的过约束代数方程组,应用最小二乘法求解过约束方程组,得到平面弹性问题位移和应力数值解。本文提供的5个规则区域的数值算例和4个不规则区域的数值算例结果表明:重心Lagrange插值配点法和重心插值正则区域法的运用,可以有效的解决规则区域和不规则区域的平面弹性问题。重心Lagrange插值配点法不仅计算公式简单、节点适应性好、程序通用性强、而且计算精度非常高。