平面强化有限单元法的h型网格自适应.pdf

第 45卷第 12期 浙 江 大 学 学 报 (工学版) JournalofZhejiangUniversity(EngineeringScience) Vo1．45NO．12 2011年 12月 Dec．201l D()I：10．3785／j．issn．1008—973X．2011．12．013 平面强化有限单元法的h型网格 自适应 凌道盛，l、令方 ，黄根清，黄 博 (浙江大学 软弱土与环境土工教育部重点实验室，岩土工程研究所 ，浙江 杭州 310058) 摘 要 ：自适应有限元 中的单元细划会不可避免地弓i入非规则节点，针对非规则节点的多点约束处理方法存在 自 由度多、数值稳定性差等缺点 ，根据强化有限单元法 (FEM i)中数学 网格和物理网格分离的特点，构造一种能够 统一处理任意阶非规则节点的位移模式关联法则——物理模式重构法 ，并导}H相应 的有限元列式．修正 了常规有 限元 中的超收敛应力恢复法和 ZZ(Zienkiwicz—Zhu)后验误差估 计理论 ，提 出强化有限单元法 的后验误差估计方法 及相应 的 h型 网格 自适应策 略．物理模式重构法在单元层次上直接将非规则节点的位移用具有独立 自由度 的数学 节点位移线性表达。既保证了单元间的位移连续性，又不需要约束处理的多余 自由度，保证了数值稳定性．数值算 例表 明：强化有 限单元法 的后验误差估计方法及 自适应策略能给 出满足精度要求 的网格． 关键词 ：强化有限单元法 ；位移模式关联法则；后验误差估计 ；h型网格 自适应 ；非规则节点 中图分类号 ：0 343．1；TU 43 文献标志码 ：A 文章编号 ：1008—973X(2O11)12—215O—O9 h—adaptiveenhanced finiteelementmethod forplaneproblems IING Dao—sheng，BU Ling～fang，HUANG Gen—qing，HUANG Bo (MOEKeyLaboratoryofSoJtSoilsandGeoenvironmentalEngineering，InstituteoJGeotechnialEngineering， ZhejiangUniversity，Hangzhou310058，China) Abstract：The multipointconstraintmethod，which iscurrently applied to dealwith thehanging node introduced by the h—refinementin standard finite element method，requiresmore degrees offreedom (DOFs)and degradesthenumericalstability，therefore，aconvenientmethodto treatthehanging node waspresentedintheframeworkoftheenhancedfiniteelementmethod (FEM 。‘)．Accordingtoseparated mathematicalandphysicalmeshes，anew correlationruleforFEM ，namelyreconstruction ofphysical displacementmode，wasproposedtohandlethehangingnodewithanyirregularityindexinauniform way， andthecorrespondingfiniteelementformulaswerethenderived．Thesuperconvergentpatchrecoveryand ZZ (Zienkiwicz—Zhu)errorestimatemethodwereim