学位论文 > 优秀研究生学位论文题录展示

基于化学流变学的树脂传递模塑工艺过程的数值模拟及灵敏度分析

作　者: 丁妍羽
导　师: 孙胜
学　校: 山东大学
专　业: 材料加工工程
关键词: 树脂传递模塑边缘效应化学流变数值模拟灵敏度分析
分类号: TQ320.66
类　型: 博士论文
年　份: 2011年
下　载: 117次
引　用: 0次
阅　读: 论文下载

内容摘要

近些年来,树脂传递模塑工艺(Resin Transfer Molding, RTM)是一种发展很快的复合材料制备工艺。它是指先把按性能和结构要求设计好的增强纤维材料预制体制件置放在模具型腔内,然后使用注射设备在一定压力下,将低粘度树脂体系注入模具型腔中,树脂固化后成型为有纤维预制件支撑的产品工件。由于具有低成本和高生产效率等优势,RTM工艺在很多领域的应用在日益扩大,如航空航天、舰船、汽车、建筑等。树脂充填过程是RTM工艺成型过程中的关键环节,模具结构、纤维与树脂的材料性质、加热方式、边缘效应、树脂固化反应现象、工艺参数等都会影响产品的性能和质量。因此,综合考虑各种影响因素,深入研究树脂浸润机理以及动力学过程,将有助于建立复合材料制品制造的缺陷控制体系,从而提高产品质量、降低生产成本。本文以RTM工艺树脂流体的充填过程为研究对象,应用多门类学科知识,建立并分析在复杂外场条件下材料体系、工艺条件、材料结构形成与充填行为之间的数学函数关系；耦合求解不同的化学与物理场量方程；自主开发等温和非等温条件下树脂充流动过程的数值模拟核心程序；引入灵敏度分析方法,定量分析工艺及材料性能参数对树脂流场的影响程度及规律。主要工作与结论如下：针对纤维预制体具有双级多孔介质的特点,采用局部体积平均化理论推导了预制体内树脂的流动控制方程,此动量方程考虑了树脂流体惯性力与粘性力的作用,比一般学者使用的Darcy定律的应用范围更广；基于矩形截面通道的完全发展流动控制方程以及等效渗透率的数学表达式,改进了标准Navier-Stokes方程,获得了边缘区域内树脂的流动控制方程,此模型考虑了模腔厚度对树脂流场的影响；根据大分子链的递归性质,应用全概率理论,推导了逐步聚合反应和链式聚合反应时树脂体系重均分子量的数值计算式。上述方程结合树脂固化反应动力学模型以及化学流变场的基本控制方程,便可实现树脂浸润过程的数学建模。为了减少计算复杂性,采用单区域方法,树脂在纤维预制体内和边缘区域内的流动采用一组通用控制方程描述,此方法可避免显式列出这两个区域界面处的边界条件；考虑树脂相与空气相之间的相互作用,流体流动采用两相(树脂相和空气相)流动模型来处理,此方法可避免显式列出树脂流动前沿处的边界条件；应用有限体积法推导了相关场量控制方程的离散表达式；引入SIMPLE算法以及交错网格技术,便可实现压力场与速度场的耦合求解；采用直接积分法或向后差分方法,实现了固化反应动力学方程的离散；采用VOF/PLIC方法实时追踪充填过程中的树脂流动前沿的推进过程。上述方法为完整实施RTM工艺充填阶段的数值模拟提供关键技术支持。编写了RTM工艺树脂充填过程的数值模拟核心程序；采用唯象模型和分子模型分别描述两种树脂体系材料(双马来酰亚胺树脂体系和环氧树脂体系)在充填过程中的化学流变行为；模拟分析了树脂固化反应程度、相对分子量、流体粘度、树脂流场随时间的演变规律,以及固化反应、流体温度、充填速度和渗透率对树脂流动前沿、模腔内化学—物理场量分布的影响,获得许多有意义的认识。为了定性和定量地研究工艺及材料性能参数对RTM充填过程的影响,引入灵敏度分析方法,推导了相关流场物理量的灵敏度方程以及数学计算式,建立了树脂流动前沿形状的灵敏度表达式,确定了灵敏度方程与树脂流场控制方程的耦合求解方法；数值分析了在恒压和恒流速入口边界条件下,树脂温度、流体速度、注射压力、纤维预制体渗透率及边缘区域渗透率对充填时间和树脂流动前沿影响规律和程度,为提高充填效率与产品质量提出了有价值的建议。分析了RTM工艺充模过程中温度场与树脂流场之间的关系特点,考虑厚度方向上的热传导现象,建立了二维流场/三维温度场的耦合求解方法,推导了三维能量方程在空间域和时间域上的离散形式,编写了非等温充填过程数值模拟核心程序；模拟分析了空间节点温度、树脂固化反应程度及粘度随时间的演变规律,探讨了边缘效应以及热传导系数对相关物理量空间分布的影响规律。结果表明,边缘效应会加剧模腔内相关物理场量的不均匀分布。

全文目录

摘要  13-15
ABSTRACT  15-18
符号说明  18-21
第一章绪论  21-35
  1.1 引言  21
  1.2 树脂传递模塑成型(RTM)工艺概述  21-24
    1.2.1 RTM工艺的特点  21-22
    1.2.2 RTM工艺的应用  22-24
    1.2.3 RTM工艺过程缺陷产生的机理分析  24
  1.3 RTM工艺相关工作的研究进展  24-31
    1.3.1 RTM工艺充模过程的研究现状  24-28
    1.3.2 边缘效应研究现状  28-30
    1.3.3 树脂化学流变行为的研究进展  30
    1.3.4 RTM工艺充模过程灵敏度分析的研究进展  30-31
  1.4 课题的提出及研究方案  31-33
    1.4.1 课题的提出  31-32
    1.4.2 研究方案  32-33
  1.5 本文主要研究内容  33-35
第二章基于化学流变学的RTM工艺充模过程的数学建模  35-52
  2.1 引言  35-36
  2.2 树脂流动的基本控制方程  36-42
    2.2.1 纤维预制体内树脂流动控制方程的建立  36-40
    2.2.2 边缘区域内树脂流动控制方程的建立  40-42
  2.3 固化反应动力学模型  42-43
  2.4 树脂结构场的基本控制方程  43-48
    2.4.1 逐步聚合反应聚合物体系的重均分子量  44-46
    2.4.2 链式聚合反应聚合物体系的重均分子量  46-48
  2.5 化学流变模型的建立  48-51
    2.5.1 唯象模型  49-50
    2.5.2 分子理论模型  50-51
  2.6 本章小结  51-52
第三章 RTM工艺充模过程数值模拟的关键技术研究  52-72
  3.1 引言  52
  3.2 有限体积法的区域离散化  52-53
  3.3 动量方程在交错网格中的离散  53-59
    3.3.1 动量方程的通用形式  53-56
    3.3.2 动量方程在时间域与空间域上的离散  56-59
  3.4 SIMPLE算法:压力场—速度场的耦合求解  59-61
    3.4.1 压力修正方程  59-61
    3.4.2 SIMPLE方法的具体实施步骤  61
  3.5 固化反应动力学方程的数值计算方法  61-62
  3.6 代数方程组的求解  62-64
    3.6.1 三对角阵算法(TDMA)  63
    3.6.2 块修正技术  63-64
  3.7 自由表面追踪技术—流体体积函数法(VOF法)  64-68
    3.7.1 VOF方法原理  64-65
    3.7.2 运动界面重构  65-66
    3.7.3 运动界面推进  66-68
  3.8 RTM工艺充模过程数值模拟流程  68-71
  3.9 本章小结  71-72
第四章双马来酰亚胺树脂体系与环氧树脂体系RTM工艺等温充模过程的数值模拟  72-94
  4.1 引言  72
  4.2 控制方程的建立  72-77
    4.2.1 树脂流动控制方程的建立  72-74
    4.2.2 固化反应动力学模型及化学流变模型的建立  74-77
  4.3 初边界条件的设定及基本假设  77-78
    4.3.1 初边界条件  77-78
    4.3.2 基本假设  78
  4.4 数值模拟流程  78-79
  4.5 模拟验证  79-81
  4.6 RTM工艺充模过程的模拟分析  81-92
    4.6.1 BMI树脂体系的模拟结果与分析  81-87
    4.6.2 EPOXY树脂体系的模拟结果与分析  87-92
  4.7 本章小结  92-94
第五章 RTM工艺充模过程的灵敏度分析  94-118
  5.1 引言  94
  5.2 灵敏度分析的理论基础  94-96
  5.3 树脂流动的控制方程  96-97
  5.4 灵敏度方程的建立  97-99
    5.4.1 基本物理量灵敏度方程的推导  97-99
    5.4.2 树脂流动前沿形状灵敏度方程的建立  99
  5.5 边界条件  99-100
  5.6 控制方程的求解  100-102
  5.7 灵敏度分析的实施步骤  102
  5.8 算例验证  102-110
    5.8.1 恒压入口边界条件算例  104-107
    5.8.2 恒流速入口边界条件算例  107-110
  5.9 灵敏度分析的模拟算例与结果讨论  110-116
    5.9.1 恒压入口边界条件算例  111-114
    5.9.2 恒流速入口边界条件算例  114-116
  5.10 本章小结  116-118
第六章温度场-流场耦合的RTM工艺非等温充模过程数值模拟  118-138
  6.1 引言  118
  6.2 控制方程的建立  118-121
  6.3 控制方程的求解  121-124
    6.3.1 网格划分及方程离散  121-123
    6.3.2 边界条件  123-124
    6.3.3 代数方程的求解  124
  6.4 RTM工艺非等温充模过程的数值模拟流程  124
  6.5 模拟结果与实验结果的比较  124-130
  6.6 算例数值分析与结果讨论  130-137
    6.6.1 边缘效应对温度场的影响  130-135
    6.6.2 热传导系数对温度场的影响  135-137
  6.7 本章小结  137-138
第七章结论与展望  138-142
  7.1 结论  138-141
  7.2 展望  141-142
参考文献  142-155
致谢  155-156
攻读博士学位期间完成的论文  156-157
攻读博士学位期间参与的科研项目  157-159
附录  159-176
学位论文评阅及答辩情况表  176

基于化学流变学的树脂传递模塑工艺过程的数值模拟及灵敏度分析

内容摘要

全文目录

相似论文