基于矩形独立覆盖初步实现结构静力分析的自动计算-无机材料学报.pdfVIP

　 第３３卷第２期 长　 江　 科　 学　 院　 院　 报 Ｖｏｌ．３３　 Ｎｏ．２　 　 　 　 ２ ０ １６ 年 ２ 月 Ｊｏｕｒｎａｌ ｏｆ Ｙａｎｇｔｚｅ Ｒｉｖｅｒ Ｓｃｉｅｎｔｉｆｉｃ Ｒｅｓｅａｒｃｈ Ｉｎｓｔｉｔｕｔｅ Ｆｅｂ． ２ ０ １ ６　 　 － 　 　 ｄｏｉ：１０．１１９８８／ ｃｋｙｙｂ ２０１６，３３（０２）：１４４ １５０ 基于矩形独立覆盖初步实现 结构静力分析的自动计算 苏海东，龚亚琦，颉志强，祁勇峰 （长江科学院材料与结构研究所，武汉　 ４３００１０ ） 摘　 要：采用前期研究的基于矩形独立覆盖的新型数值流形法，提出结构线弹性静力分析的自动计算方法，包括自 动的前处理、自适应分析等。 根据独立覆盖的特点提出几个后验误差指标：独立覆盖之间条形连接区域的应变连 续性指标；边界应力指标和独立覆盖的高阶误差指标。 利用新方法的ｈ型网格加密及ｐ型升阶的方便性，选择一 － 种路径尝试ｈ ｐ型的混合自适应，其中，对于矩形独立覆盖采用简单的二分法实现覆盖加密。 通过几个二维算例 验证了新方法实现自动计算的可行性，只需人工输入结构外形、材料参数和边界条件，其它工作完全交由计算机完 成，最终得到满足一定精度的计算结果。 － 关键词：数值流形方法；独立覆盖；自动计算；误差估计；ｈ ｐ混合自适应；静力分析 － － － 中图分类号：ＴＢ１１５；ＴＶ３１１　 　 　 文献标志码：Ａ　 　 　 文章编号：１００１ ５４８５（２０１６）０２ ０１４４ ０７ 通过结点的正确连接，而实际工程结构的复杂性使 １　 研究背景 得网格剖分的难度较大。 同时，由于求解域内物理 场分布的复杂程度不同，需要网格密度的变化，比如 计算机辅助工程 （Ｃｏｍｐｕｔｅｒ Ａｉｄｅｄ Ｅｎｇｉｎｅｅｒｉｎｇ， 目前广泛使用的ｈ型自适应需要对精度不佳的区域 ＣＡＥ），主要是指用计算机辅助求解复杂工程和产品 进行网格细分，涉及到的大小网格过渡通常不太方 的结构强度、刚度、动力响应、弹塑性等力学性能的 便，要保持网格合适形状较为不易。 在计算收敛效 ［１］ 分析计算 。 提高分析效率对ＣＡＥ意义重大，未来 率方面，对于存在强奇异性的裂纹尖端（以下简称 的发展趋势是让计算机自动、高效地完成分析，即 裂尖）等复杂应力场分布区域，用常规的有限单元 ＣＡＥ分析技术的自动化。 插值函数逼近往往效果不佳。 本文提出自动计算的思想，其目标是基本做到 广义有限元在结点上引入级数提高近似函数 无需人工参与就能获得理想结果的分析，当然，必要 逼近能力，但其主要问题是高阶近似函数情况下的 － 的结构外形、材料参数和边界条件的人工输入除外。 ［３ ４］ 方程组线性相关 。 扩展有限元在单元结点上增 要达到此目标，必须解决 ＣＡＥ前处理的自动化、计