双尺度纤维织物二维非饱和流动的数值模拟与实验

液体模塑成型工艺(LCM)中非饱和流动的填充模拟对于在虚拟空间中快速、高效地优化工艺参数具有重要意义。采用了一种模拟双尺度纤维织物在等温条件下非饱和流动的双尺度计算模型,通过引入沉浸函数求解宏观-微观流动控制方程组,同时考虑了在微观浸渍中毛细压力的影响,在有限元/控制体积算法中实现了对非饱和流动的数值模拟。随后对三向缝合纤维织物进行了二维径向填充实验,将实验结果与数值模拟的预测值对比。结果表明,该计算模型可以较精确地模拟双尺度纤维织物中的非饱和流动。在此计算模型的基础上,讨论了流体黏度、注射流量及纤维束孔隙率对非饱和填充浸润的影响。结果表明,不同流体黏度、注射流量及纤维束孔隙率对纤维织物填充过程中非饱和区域长度、入口压力曲线及填充时间影响不同。研究结果可以对合理预测纤维织物的浸润及树脂填充过程中入口压力提供指导。     

液体模塑成型工艺二维径向非饱和流动数值模拟

通过引入沉浸函数建立了双尺度多孔介质非饱和流动模型,并采用有限元/控制体积法实现了恒压及恒流注射条件下液体模塑成型（LCM）工艺二维径向非饱和流动的数值模拟,得到了不同注射条件下纤维织物内的压力场分布及半饱和区域长度随时间的变化规律,并将双尺度非饱和理论结果与单尺度饱和理论结果进行对比。结果表明：非饱和流动过程中,半饱和区域内的压力和压力梯度明显下降;半饱和区域长度随时间逐渐增加随后保持稳定,当流动前沿到达出口后半饱和区域长度开始逐渐减小;当两个主方向渗透率不同时,沿主方向半饱和区域长度也不同,渗透率越大该方向的半饱和区域长度也越大,纤维织物完全浸润时间取决于较小的渗透率。研究结果对合理预测树脂填充过程中压力分布及纤维预制件的浸润具有指导意义。     

恒压下树脂在双尺度多孔介质中的流动特性

基于复合材料液态模塑(LCM)工艺过程中存在半饱和区域的实验现象以及对预制体双尺度效应的逐步认识, 一些学者提出用沉浸模型来研究双尺度多孔介质的不饱和流动。通过体积均匀化方法描述了双尺度多孔介质复合材料液态模塑工艺模型的特征, 得到含有沉浸项的双尺度多孔介质的质量守恒方程, 并采用有限元法对方程进行数值求解, 通过具体算例计算了考虑双尺度效应时恒压树脂注射下不同时段的压力分布状态, 得到树脂在填充过程中流动前沿半饱和区域从出现到消失的过程, 采用不同注射压力进行模拟并比较。结果表明, 与单尺度多孔介质模型不同, 双尺度多孔介质模型更能反映实际树脂填充过程中出现的半饱和区域现象。     

基于“离位”增韧技术Z向注射RTM成型的浸润研究

针对"离位"增韧技术和Z-RTM成型技术,引入饱和度参数修正Darcy定律,建立描述树脂在纤维预制件中非稳态流动的偏微分方程,研究恒流注射过程中体积流量、树脂黏度和纤维预制件渗透率等工艺参数对非稳态浸润过程注入压力的影响,模拟树脂在层间未增韧和增韧纤维预制件束内和束间的流动。结果表明:数值模拟结果具有可靠性;随着注射时间的增加,纤维预制件内部各点的压力增加;随着体积流量、树脂黏度的增加,注入压力线性增加,而随着纤维渗透率的增加,注入压力减少,符合Darcy定律;实现了树脂在纤维预制件细微观层次浸润的可视化,这种可视化结果为预测树脂在预制件中的宏观流动提供了重要补充,并为实际工艺提供了一定指导作用。     

层间纳米纤维膜对玻纤预制体渗流特性的影响

为了研究层间纳米纤维膜对玻纤织物渗流特性的影响,使用超景深三维显微镜表征了纳米纤维含量对玻纤织物微观结构的影响,采用径向法测量了纳米纤维膜夹层玻纤织物预制体的渗透率,重点分析了纳米纤维含量对玻纤织物预制体渗流模式的影响。结果表明:玻璃纤维束间的毫米尺度区域被纳米纤维膜填充而离散成微米尺度区域;预制体孔隙率及渗透率值均随着纳米纤维含量的增加而减小;随着纳米纤维含量的增加,复合预制体表现出的各向异性程度逐渐减小;树脂宏观流动前沿内部分饱和区域面积比例随纳米纤维含量的增加而增大;相同纳米纤维含量预制体的部分饱和区域面积比例随注入时间的增加呈先增大后减小趋势。     

树脂传递模塑成型工艺中嵌套效应引起渗透率变异的实验与数值模拟

树脂传递模塑成型工艺(RTM)中最重要的变形模式之一是厚度方向压缩。厚度方向压缩减小了织物预成型体的厚度,使织物预成型体局部结构形式发生改变从而引起嵌套效应。嵌套效应不仅会减少织物预成型体的厚度,增加纤维的体积分数并改变孔隙率,而且相邻织物层嵌套效应具有一定的空间分散性,从而使得织物预成型体渗透率具有变异性。本文针对低黏度树脂设计了一种实验装置用以测量局部渗透率的空间分散性,随后建立了随机嵌套单胞模型,利用ANSYS/CFX有限元软件实现了单胞填充浸润的数值模拟,通过流量分析获得局部渗透率,并研究了渗透率的统计分布。通过实验结果与数值模拟结果相对比,验证数值模拟结果的可靠性。最后,基于渗透率的统计分布建立了随机渗透率场,并进行填充浸润的数值模拟,通过与传统恒定渗透率的方法进行比较,证明该方法具有更高的先进性。研究结果可以对未来RTM工艺的稳健性优化提供依据。      

非饱和多孔介质动力应变局部变化数值模拟

在已有研究工作基础上对非饱和多孔介质应变局部化问题进行研究,给出非饱和多孔介质的分析控制方程,其中饱和度与毛细压力关系由实验给出。采取适用于非饱和砂土的改进的广义塑性本构模型对应局部变化过程进行数值模拟,给出了试件应变局部化发展过程以及孔隙压力的变化规律。对初始饱和土中所产生非饱和剪切带进行计算的结果表明,采用非饱和模型较饱和模型将获得更为合理的结果。     

测试条件对集束纤维单向渗透率的影响

为掌握渗透率测试的影响因素,自行研制了单向纤维集束的浸渗特性测试系统,对比了液体沿单向纤维集束的饱和与非饱和渗透特性,分析了不同压力作用形式、纤维体积分数、液体种类对非饱和渗透率测试结果的影响。结果表明：相同条件下纤维集束的饱和渗透率大于非饱和渗透率测试结果;当外压较小时毛细作用对非饱和渗透率测试结果的影响显著;当纤维体积分数增加时集束纤维的渗透率呈减小趋势;相同压力差下真空驱动测得的渗透率比正压驱动测试结果大。     

树脂传递模塑工艺中的非饱和流动过程模拟与实验研究

针对编织类纤维增强体的纤维束之间与纤维束内孔隙的双尺度特点,建立了平纹织物的细观结构模型,并推导了汇函数的数学表达式。建立了局部细观流动特征的非饱和流动控制方程,利用有限元/控制体积方法求解,得到了局部饱和度分布。与实验进行比较,吻合较好。      

基于壁面滑移的薄壁注塑流动三维数值模拟

基于壁面滑移速度模型,且考虑黏度的压力依赖性,运用三维有限元数值模拟技术,对薄壁注塑流动过程进行数值模拟,研究壁面滑移对厚度方向熔体速度分布、注射时间、注射压力和锁模力的影响。结果表明,壁面滑移对薄壁注塑流动过程具有重要的影响。通过模拟结果与实验结果进行对比,以验证数值模拟结果的可靠性。结果还表明,只有同时考虑壁面滑移效应和黏度的压力依赖性,才能更准确地描述薄壁注塑流动过程。     

The Sink Term to Different Kinds of Fibrous Mat During the Unsaturated Filling Process

The fibrous pre-form of resin transfer molding is a dual-scale porous medium with two distinct scales of pores, i.e., pores in intra- and inter-tow, which produce an unsaturated infiltration phenomenon during filling. A sink term representing the delayed flow rate from the inter-tow gap into the intra-tow one is introduced to establish governing equations. This study mainly analyzes the sink term by tow saturation during the microscopic flow. First, fiber-tow permeability is calculated by FLOTRAN of ANSYS, Second, periodic unit cells are built according to different structures, and the concrete expression of the sink term is indirectly obtained through the numerical simulation and date fitting of tow saturation under different pressure and viscosity conditions. Results indicate that: the FLOTRAN module can be used to calculate the permeability of fiber tow in two directions; Moreover, the filling time and infiltration process for diverse unit cells with the same volume fraction are different; under the same injection condition, different unit cells have different parameters for the sink term.     

Application of calcium carbonate in resin transfer molding process: An experimental investigation

The resin transfer molding (RTM) process is used to manufacture advanced composite materials made of continuous glass or carbon fibers embedded in a thermoset polymer matrix. In this process, a fabric preform is prepared, and is then placed into a mold cavity. After the preform is compacted between the mold parts, thermoset polymer is transferred from an injection machine to the mold cavity through injection gate(s). Resin flows through the porous fabric, and eventually flows out through the ventilation port(s). After the resin cure process (cross‐linking of the polymer), the mold is opened and the part is removed. The objective of this study is to verify the application of calcium carbonate mixed in resin in the RTM process. Several rectilinear infiltration experiments were conducted using glass fiber mat molded in a RTM system with cavity dimensions of 320 × 150 × 3.6 mm, room temperature, maximum injection pressure 0.202 bar and different content of CaCO₃ (10 and 40%) and particle size (mesh opening 38 and 75 µm). The results show that the use of filled resin with CaCO₃ influences the preform impregnation during the RTM molding, changing the filling time and flow front position, however it is possible to make composite with a good quality and low cost.     

Analysis of non-isothermal mold filling process in resin transfer molding (RTM) and structural reaction injection molding (SRIM)

In this paper, we present a modeling and numerical simulation of a mold filling process in resin transfer molding/structural reaction injection molding utilizing the homogenization method. Conventionally, most of the mold filling analyses have been based on a macroscopic flow model utilizing Darcy's law. While Darcy's law is successful in describing the averaged flow field within the mold cavity packed with a porous fiber preform, it requires experiments to obtain the permeability tensor and is limited to the case of porous fiber preform-it can not be used to model the resin flow through a double porous fiber preform. In the current approach, the actual flow field is considered, to which the homogenization method is applied to obtain the averaged flow model. The advantages of the current approach are: parameters such as the permeability and effective heat conductivity of the impregnanted fiber preform can be calculated; the actual flow field as well as averaged flow field can be obtained; and the resin flow through a double porous fiber preform can be modelled. In the presentation, we first derive the averaged flow model for the resin flow through a porous fiber preform and compare it with that of other methods. Next, we extend the result to the case of double porous fiber preform. An averaged flow model for the resin flow through a double porous fiber preform is derived, and a simulation program is developed which is capable of predicting the flow pattern and temperature distribution in the mold filling process. Finally, an example of a three dimensional part is provided.     

Modeling of heat transfer and unsaturated flow in woven fiber reinforcements during direct injection-pultrusion process of thermoplastic composites

This paper provides a methodology for the modeling of heat transfer and polymer flow during direct thermoplastic injection pultrusion process. Pultrusion was initially developed with thermosets which have low viscosity. But the impregnation becomes a critical point with thermoplastics which exhibit higher viscosity. There are very few reported works on direct thermoplastic impregnation with injection within the die. In addition, the rare studies have not adequately addressed the issue of unsaturated flow in woven fiber reinforcements. The solution proposed here, models the polymer flow through dual-scale porous media. A heat transfer model is coupled to a flow model enriched with a sink term. Specific changes of variables are made so as to model the steady state solution of unsaturation along a continuous process. The sink term, added to the continuity equation, represents the absorption rate of polymer by the bundles. Data were measured on a pultrusion line and micrographs confirmed the modeling strategy with an unsaturated flow approach. The flow modeling coupled to heat transfer of the thermoplastic pultrusion process aims at determining the saturation evolution through the die so as to manufacture pultruded profiles with the lowest residual porosity.     

A numerical study of filling process through multilayer preforms in resin injection/compression molding

In resin injection/compression molding (RI/CM), a preform often comprises layers of different fiber reinforcements. Each fiber reinforcement has unique through thickness and in-plane permeabilities as well as compressibility, creating a heterogeneous porous medium in the mold cavity. In the present article, numerical simulation is utilized to investigate the filling process of RI/CM in such a heterogeneous porous medium. The filling stage is simulated in a full three-dimensional space by using control volume/finite element method and based upon an appropriate filling algorithm. The flow in the open gap which may be present in the mold cavity is modeled by Darcy’s law using an equivalent permeability. Numerical simulations of filling process for preforms containing two and three layers of different reinforcements in various stacking sequences are conducted with the aid of computer code developed in this study. Results show that the injection time as well as flow front progression depends on fiber types in the whole preform, fiber stacking sequence and open gap provided in the mold cavity. Simulated results also suggest that the presence of open gap at top of reinforcement can lead to both low injection time and uniform flow pattern.     

基于多相流的LCM工艺气泡生成仿真

分析了液体模塑成型工艺（LCM）下织物预成型体中双尺度流动以及由此造成的空气裹入,进而产生细观及微观气泡的现象。基于多相流（VOF）方法建立了树脂空气两相流体在单胞内部流动的数学模型,并确定了该模型中多孔介质阻力源项和毛细力源项的具体形式。基于Fluent软件的UDF功能实现了上述两相流模型的数值求解,研究了平纹织物单胞中的两相流动以及2种气泡的生成过程。对Rovcloth 2454织物的气泡生成仿真结果显示,毛细数 Ca 对气泡的产生有决定性作用：当毛细数接近临界毛细数 Ca^c时,气泡产生量最低,而当Ca小于Ca^c时,容易产生细观气泡,反之容易产生微观气泡。通过与文献中的理论预测和实验数据对比,验证了本文算法的正确性。