不规则区域平面弹性问题的正则区域重心插值配点法.pdf

卷 期 棳 棾椀 棽 计算 力 学 学 报 斨旓旍棶棾椀 斘旓棶棽 年 月 旇旈旑斿旙斿斒旓旛旘旑斸旍旓旀 旓旐旛旚斸旚旈旓旑斸旍斖斿斻旇斸旑旈斻旙 棽棸棻椄 棿 斄 旘旈旍棽棸棻椄 旔 旔 椇 王兆清 斈斚斏棻棸棶椃椀棻棻 旙旍旞棽棸棻椂棻棽棸椃棸棸棻 旉 不规则区域平面弹性问题的 正则区域重心插值配点法 灣棻 棻 棻 棽 王兆清 棳 纪思源 棳 徐子康 棳 金 棬 棳 椈 棳 棭 山东建筑 学 工程力学研究所 济南 山东建筑 学 理学院 济南 棻棶 棽椀棸棻棸棻 棽棶 棽椀棸棻棸棻 椇 棳 摘 要 将不规则区域嵌入到规则的矩形区域 在矩形区域上将弹性平面问题的控制方程采用重心斕斸旘斸旑斿插 旂 旂 棳 暎 棳 值离散 得到控制方程矩阵形式的离散表达式 在边界节点上利用重心插值离散边界条件 规则区域采用置换法 棳 棳 棳 施加边界条件 不规则区域采用附加法施加边界条件 得到求解平面弹性问题的过约束线性代数方程组 采用最 棳 暎 小二乘法进行求解 得到整个规则区域上的位移数值解 利用重心插值计算得到不规则区域内任意节点的位移 棴棻棿 棳 暎 暎 值 计算精度可到棻棸 以上 数值算例验证了所建立方法的有效性和计算精度 椇 椈 椈 椈 椈 椈 关键词 不规则区域 平面弹性问题 正则区域法 重心斕斸旘斸旑斿插值 微分矩阵 配点法 旂 旂 中图分类号椇 文献标志码椇 文章编号椇 棬 棭 斚棽棿棻 斄 棻棸棸椃灢棿椃棸椄棽棸棻椄棸棽灢棸棻椆椀灢棸椃 棳 暎 似展开 数值求解微分方程的配点法 通过微分矩 棻 引 言 棳 阵来表达未知函数的各阶导数 微分矩阵通过对未 求解弹性力学平面问题的主要方法有解析解