97-欧拉有限元及无网格自适应耦合模拟.pdfVIP

No．7 第43卷第7期 机械工程学报 v01．43 2007年7月 Jul． 2007 CHINESEJoURNALOFMECHANICALENGINEERING 成形充填过程的任意拉格朗日一欧拉 有限元与无网格自适应耦合模拟· 段庆林李锡夔 (大连理工大学工业装备结构分析国家重点实验室大连116023) 摘要：在任意拉格朗日一欧拉(劬i仃a叫La鲫ngi锄．Eulcri姐，√虹凰)方法框架中发展用于成形充填过程数值模拟的 有限元与无网格自适应耦合方案。该方案自适应地在网格发生扭曲的区域采用无网格法空间离散以保证求解的精 度和稳定性，而在网格质量较好的区域以及本质边界上保留使用有限元法空间离散以提高计算效率和便于施加本 质边界条件，从而在充分发掘有限元法和无网格法各自优点的同时避开它们各自的缺点，在一定程度上保证该方 法同时具有较快的计算速度和较好的健壮性。采用ALE参考构形描写充填过程中的熔体质量运动及守恒，可在 不频繁变动网格的前提下准确跟踪自由面。应用压力稳定型分步算法求解控制方程，可方便采用不满足 该方法相对于单一的有限元法和无网格法的优越性，以及对含自由面变质量体系流动问题数值模拟的有效性。 关键词：ALE方法成形充填无网格自由面追踪 自适应耦合方法 中图分类号：0357．1 和物质点，则可在不频繁变动网格的前提下精确跟 O前言 踪自由面。 问题②的困难在于不可压缩模型要求速度、压 a_buSka． 成形充填过程代表一类含自由界面的变质量 力的插值空间须满足所谓Lady吐％skaja-B 体系流动问题，对此类过程实施有效数值模拟是发 Bre盟i(LBB)条件，否则将导致压力场的虚假数值振 展注塑模计算机辅助制造技术亟待解决的核心问题 荡。尽管有若干种插值模式能够满足该条件，但更 之一。为此须解决以下两个关键问题：①精确跟踪 便于应用的速度、压力的同低阶插值模式并不满足 熔体自由面并形成相应计算网格。②在给定计算网 了各自的稳定化方案。值得注意的是，LI等【6】在有 格上准确求解控制方程。 限增量微积分的理论框架下，通过引入一个附加变 ’问题①的困难在于任何时刻要保证自由面是 量，发展了压力稳定型分步算法，有效消除了压力 一个包含熔体阵面上粒子的物质面，且要满足自由 场的虚假振荡。 ． 面上的法向和切向应力恒为零(如假定不存在表面 应着重指出的是，虽然ALE方法能够在一定程 张力1。自由面跟踪方法可分类为Euleri蛆方法、 度上解决困扰Lag啪gi髓方法的网格扭曲问题，但 Lag啪gi锄方法和Eul耐an-La酎强gi趾方法。 在流体运动发生剧烈变化的区域仍将不可避免地出 Eul嘶瓶方法尽管具有域内网格固定的优点，但易