首页> 中文学位 >成型充填过程中非等温非牛顿粘性流动的ALE有限元与无网格自适应耦合模拟
【6h】

成型充填过程中非等温非牛顿粘性流动的ALE有限元与无网格自适应耦合模拟

代理获取

目录

文摘

英文文摘

声明

1绪论

1.1课题的工程背景及研究意义

1.2控制方程的数值求解

1.2.1克服LBB条件的数值方法

1.2.2对流占优问题的数值方法

1.3自由面追踪技术研究现状

1.3.1基于欧拉描述的方法

1.3.2基于拉格朗日描述的方法

1.3.3基于任意拉格朗日欧拉描述的方法

1.4无网格法简介

1.4.1无网格法的发展历程

1.4.2无网格法评价

1.5本文的主要工作

参考文献

2伽辽金型无网格法介绍

2.1引言

2.2线弹性力学的控制方程及其Galerkin弱形式

2.3移动最小二乘法(Moving Least Square,MLS)

2.3.1 MLS基本原理

2.3.2节点的权函数和影响域

2.3.3 MLS形函数及其导数的形成过程

2.3.4形成过程的改进

2.4本质边界条件的施加

2.4.1拉格朗日乘子法

2.4.2修正变分原理法

2.4.3罚函数法

2.4.4 Nitsche法

2.4.5位移约束方程法

2.4.6与有限元耦合法

2.5区域积分的实现

2.5.1背景格子结构积分(Integration with background cell structure)

2.5.2背景网格积分(Integration with background mesh)

2.5.3节点积分

2.5.4移动最小二乘积分

2.6本章小结

参考文献

3成型充填过程模拟的数学模型

3.1引言

3.2 ALE运动学描述

3.3 ALE描述下的控制方程组

3.4本构关系

3.5充填问题的几何描述和边界条件

3.6 ALE自由面追踪及网格生成方法

3.6.1自由面更新及碰壁控制

3.6.2实时计算网格的生成

3.7成型充填过程模拟的实施步骤

3.8本章小结

参考文献

4不可压缩N-S方程的压力稳定型分步算法

4.1引言

4.2稳态不可压缩Stokes流动问题的PSPG稳定化技术

4.2.1 PSPG稳定化求解格式的导出

4.2.2数值验证

4.3经典的分步算法(C-FSA)

4.4有限增量微积分(Finite Increment Calculus,FIC)理论

4.5基于FIC的压力稳定化的不可压缩N-S方程的导出

4.6压力稳定型分步算法(PS-FSA)的导出

4.7压力稳定性分析

4.8压力稳定型CBS分步算法

4.9数值算例

4.9.1方腔流问题

4.9.2平面Poisseuille流动问题

4.10本章小结

参考文献

5有限元与无网格的自适应耦合空间离散及其在充填模拟中的应用

5.1引言

5.2连续掺混法(Continuous Blending Method,CBM)

5.3基于CBM的区域划分方案

5.4有限元与无网格的自适应耦合空间离散方法

5.5耦合方案在充填模拟中的应用

5.6数值算例

5.6.1平面Poisseuille流动问题

5.6.2典型型腔充填过程模拟

5.7本章小结

参考文献

6非等温非牛顿粘性流的基于ALE描述的广义特征线Galerkin法

6.1引言

6.2经典的特征线Galerkin法

6.3 ALE描述下的广义特征线Galerkin法

6.4非等温非牛顿流的广义特征线法基分裂算法(PS-GCBS)

6.5数值算例

6.5.1平面Poisseuille流动问题

6.5.2 4∶1收缩流动问题

6.5.3典型型腔充填过程模拟

6.6本章小结

参考文献

7充填过程有限元与无网格自适应耦合模拟的程序实现

7.1引言

7.2数据结构说明

7.3主程序流程说明

7.4各模块流程说明

7.4.1自由面追踪和网格生成模块

7.4.2有限元与无网格自适应耦合插值模块

7.4.3控制方程数值求解模块

8总结与展望

8.1总结

8.2展望

参考文献

论文创新点摘要

攻读博士学位期间发表学术论文情况

致谢

展开▼

摘要

注塑成型是目前塑料成型的主要方法之一,生产的制品如磁带盒、计算机键盘等。该工艺的流程为,首先将熔融的高聚物熔体以较大的流速注射入模具型腔中,直至充满整个腔体,然后待塑料熔体充分固化后开启模具脱出塑料制品,最后再循环进入下一个制品的生产。 为提高模具的生产效率和水平,人们迫切需要发展针对该工艺过程的高效的数值分析工具来取代传统的费时费力的试模修模的过程。因而,在过去的数十年中,注塑成型过程的数值模拟引起了越来越广泛的关注。尽管如此,仍然还有很多问题有待进一步的研究,如自由面的准确追踪以及稳定、高效、精确地求解控制方程等。本论文的研究工作正是以解决这些问题为目标而展开的。 移动自由面的准确追踪对于注塑成型过程的数值模拟具有十分重要的意义。根据建立方程所采用参考系的不同,目前已有的处理方法大体上可归类为拉格朗日方法、欧拉方法以及任意拉格朗日欧拉方法(ALE)。在ALE方法中,计算网格的运动独立于物质构型的运动,结合了拉格朗日方法和欧拉方法的优点。本文第三章将给出一个基于ALE方法的自由面追踪和网格生成.重生成模型,可在准确确定移动自由面空间位置的同时尽量避免发生网格扭曲。在所发展的模型中,充填域变质量体系的实时网格生成被简化为移动自由面附近区域的多边形三角化过程,大大减少了网格处理所需的计算量。此外,还引入了局部Laplacian光顺方案来提高计算网格的质量,并发展了多种机制来处理自由面与模具边界的接触问题。 实现注塑成型模拟的另外一个关键问题是须发展能稳定求解控制方程初边值问题的数值方法。众所周知,采用标准Galerkin方法模拟不可压缩流动问题通常会导致两种类型的虚假数值振荡。其一来源于控制方程组的可约化混合特性,该特性限制了速度.压力(u-p)插值空间的自由选取。不满足所谓LBB条件的非协调的插值模式,如u-p的等低阶插值模式,将导致压力场的虚假空间振荡。第二种类型的数值振荡是由控制方程的对流特性引起的,尤其易发生在对流占优的情况下,这是由于标准的空间 Galerkin离散方法只对自伴随的算子方程有效。 对于第一个问题,本文通过在有限增量微积分 (FIC)的理论框架下重建质量守恒方程,并理性地引入一个辅助变量来避免空间高阶导数的计算,提出了压力稳定型分步算法(PS-FSA),有效消除了压力场的虚假数值振荡。与经典的分步算法(C-FSA) 相比,该算法具有更好的压力稳定性。此外,本文还采用这种压力稳定机制改善了经典的成型充填过程中非等温非牛顿粘性流动的ALE有限元与无网格自适应耦合模拟特征线基分裂算法(CBS)的压力稳定性。PS-FSA算法具体的推导过程以及它的数值验证将在第四章中给出。对于第二个问题,本文发展了基于ALE描述的广义特征线Galerkin(CG)方法,有效处理了对流占优问题控制方程的空间离散。而且,与经典的CG方法相比,该方法能使用更大的时间步长。实际上,经典的CG方法以及经典的 Crank-Nicolson (CN)方法能分别作为两个特例而被纳入到本文所提出的广义CG方法的框架中。此外,本文方法还将经典CG方法由欧拉描述下推广到ALE描述下,以便于和前述ALE自由面追踪技术相结合。应指出的是,若将参考系固定于空间不动,ALE描述可退化为欧拉描述,即相对于经典CG方法来说,本文方法更具一般性。将所发展的广义CG方法和分步算法以及第四章给出的压力稳定技术相结合,本文还提出了一个迭代型压力稳定的广义CBS算法,并将其应用于非等温非牛顿流的数值模拟中。该算法具体的推导过程以及它的数值验证将在第六章中给出。 此外,本文还提出了一个有限元与无网格离散区域剖分的自适应算法,可自动地将整个计算域划分成分别应用有限元法(FE)、无网格法(MF)以及它们的耦合插值进行空间离散的三个子域。基于该自适应算法以及连续掺混法(CBM),本文提出了有限元与无网格的自适应耦合空间离散方法,并将其应用于注塑成型的数值模拟中。所发展的自适应耦合方法可在充分结合有限元法与无网格法各自优点的同时避免它们各自的缺点。该方法具体的实现细节以及展现它相对于传统有限元法和无网格法优势的数值验证将在第五章中给出。 为了论文的完整性,本文第二章还简单介绍了伽辽金型无网格法的基本概念和它所面临的问题以及相应的处理方法。实现本文算法的计算机程序及其数据结构的说明放在第七章。第八章总结全文并展望进一步的研究内容。

著录项

相似文献

  • 中文文献
  • 外文文献
  • 专利
代理获取

客服邮箱:kefu@zhangqiaokeyan.com

京公网安备:11010802029741号 ICP备案号:京ICP备15016152号-6 六维联合信息科技 (北京) 有限公司©版权所有
  • 客服微信

  • 服务号