Analysis of the Edge Effect on Flow Patterns in Liquid Composites Molding

In liquid composites molding, the cutting of the fiber preform is often not sufficiently precise leaving a small clearance between the reinforcement and the mold edges. This clearance creates a preferential flow path for the resin which may disrupt the filling of the mold cavity. Experience has shown that even a small clearance of 1 or 2 mm could have a significant effect, especially if the preform has a high fiber content. A model is thus needed to predict this channeling effect in order to take it into account in computer simulations of the mold filling process.This paper presents a model to describe this phenomenon. The idea is to characterize simultaneously the flow in the channel and through the reinforcement. The model is derived using Navier–Stokes equations in the open channel and Darcy's law in the porous preform. From this model, an equivalent porous medium was defined for which an equivalent permeability tensor can be computed as a function of the channel geometry. Numerical simulations performed with the computer software RTMFLOT developed in our laboratory have shown that in some limit cases, i.e., a large channel (5–6 mm) or for very low porosity reinforcements, the transverse flow from the channel to the preform can be neglected in the model while still obtaining quite a good prediction of the flow edge effect. In other cases, however, namely for a clearance of intermediate size (2–3 mm) which is the most common case in RTM, or for a higher porosity of the reinforcement, the transverse flow from the channel to the rest of the preform must be taken into account. Experimental data to validate the proposed model are also presented.     

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.     

Numerical simulation of injection/compression liquid composite molding. Part 1. Mesh generation

This paper presents a numerical simulation of injection/compression liquid composite molding, where the fiber preform is compressed to a desired degree after an initial charge of resin has been injected into the mold. Due to the possibility of an initial gap at the top of the preform and out-of-plane heterogeneity in the multi-layered fiber preform, a full three-dimensional (3D) flow simulation is essential. We propose an algorithm to generate a suitable 3D finite element mesh, starting from a two-dimensional shell mesh representing the geometry of the mold cavity. Since different layers of the preform have different compressibilities, and since properties such as permeability are a strong function of the degree of compression, a simultaneous prediction of preform compression along with the resin flow is necessary for accurate mold-filling simulation. The algorithm creates a coarser mechanical mesh to simulate compression of the preform, and a finer flow mesh to simulate the motion of the resin in the preform and gap. Lines connected to the top and bottom plates of the mold, called spines, are used as conduits for the nodes. A method to generate a surface parallel to a given surface, thereby maintaining the thickness of the intermediate space, is used to construct the layers of the preform in the mechanical mesh. The mechanical mesh is further subdivided along the spines to create the flow mesh. Examples of the three-dimensional meshes generated by the algorithm are presented.     

Preform permeability variation with porosity of fiberglass and carbon mats

An experimental investigation on fiber bed permeability variation with porosity for three types of reinforcement mats is performed. The reinforcements consist of plain-weave carbon, plain-weave fiberglass, and chopped fiberglass mats. Resin flow experiments are performed in a rectangular cavity with different fiber volume fractions. RL 440 epoxy resin is used as the working fluid in the experiments. Several layers of mats are laid inside the mold in each experiment and resin is injected at a constant pressure. The effects of reinforcement type and porosity on fiber bed permeability are investigated. Fiber mat permeability of woven mats show large degrees of anisotropy. Resin flow in chopped fiberglass mats is circular, suggesting an isotropic permeability tensor. In all the three cases, preform permeability increases with fiber bed porosity in a non-linear fashion. The results of this investigation could be employed in optimization of liquid composite molding manufacturing processes.     

Modeling forces generated within rigid liquid composite molding tools. Part A: Experimental study

The term liquid composite molding (LCM) encompasses a growing list of processes, including resin transfer molding (RTM), injection/compression molding (I/CM), and resin infusion (a.k.a. VARTM). All LCM techniques involve compressive deformation of the fiber reinforcement prior to, and in many cases during mold filling. Forces acting on molds are primarily due to the requirement to compact the reinforcement, and pressure generated due to resin flow through these fibrous structures. An experimental study of the forces exerted on a mold during the RTM and I/CM processes is presented here. Two reinforcing materials have been considered, exhibiting significantly different resistance to compaction. The evolution of mold clamping force has been shown to be strongly influenced by the complex, non-elastic compaction behaviour of fiber reinforcements. The important effects include stress relaxation, an apparent lubrication by the injected fluid, and permanent deformation. Efforts to simulate the experiments will be presented in Part B of this study.     

LCM 充模过程中的边缘效应

边缘效应是复合材料液体模塑成型技术(Liquid composites molding , LCM) 中常见的纤维预成型体铺敷缺陷之一。采用单向流动法研究了边缘效应对纤维预成型体渗透率及充模过程的影响, 结合其等效渗透率的理论预测模型对不同纤维体积含量、不同缝隙宽度条件下的边缘效应进行了模拟与分析, 提出了一边缘效应强弱的表征因子, 并以一较复杂的模腔的充模过程为实例提出了对边缘效应的在线监控策略及处理方案。      

Numerical simulation of injection/compression liquid composite molding. Part 2: preform compression

In the injection/compression liquid composite molding process (I/C-LCM), a liquid polymer resin is injected into a partially open mold, which contains a preform of reinforcing fibers. After some or all of the resin has been injected, the mold is closed, compressing the preform and causing additional resin flow. This paper addresses compression of the preform, with particular emphasis on modeling three-dimensional mold geometries and multi-layer preforms in which the layers have different mechanical responses. First, a new constitutive relation is developed to model the mechanical response of fiber mats during compression. We introduce a new form of nonlinear elasticity for transversely isotropic materials. A special case of this form is chosen that includes the compressive stress generated by changes in mat thickness, but suppresses all other responses. This avoids the need to model slip of the preform along the mold surface. Second, a finite element method, based on the principle of virtual displacement, is developed to solve for the deformation of the preform at any stage of mold closing. The formulation includes both geometric and material nonlinearities, and uses a full Newton–Raphson iteration in the solution. An open gap above the preform can be incorporated by treating the gap as a distinct material layer with a very small stiffness. Examples show that this approach successfully predicts compression in dry preforms for three-dimensional I/C-LCM molds.     

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.     

注射条件对LCM工艺非饱和流动特性影响

双尺度多孔纤维预制体填充过程中延迟浸润的非饱和流动现象,对基于树脂流过区域为完全饱和区域的充模理论及模拟方法提出了挑战。通过控制体/有限单元（CV/FE）法结合沉浸函数实现了液体模塑成型工艺（LCM）中非饱和填充浸润的数值模拟,并对比了恒压下的实验结果,验证了其可靠性。分析讨论了注射口压力、流量和液体黏度对双尺度多孔纤维织物非饱和填充浸润特性的影响。结果表明:在允许误差内,该数值模拟结果可靠,可用于分析讨论各因素对双尺度多孔织物非饱和流动特性的影响;填充浸润过程中,纤维织物内部非饱和区域长度并非保持不变,而是随着填充浸润的进行经历了4个变化过程;不同注射条件下,压力、流量及黏度对非饱和流动特性影响不同。研究结果对合理控制注射条件及流体特性实现双尺度多孔纤维预制件的完全浸润具有指导意义。      

Analysis and minimization of void formation during resin transfer molding process

In the resin transfer molding process, residual air in the pores of fiber preform results in dry spots and microvoids in the finished product. The dry spots are usually formed due to irregular permeability of fiber mat and improper injection locations. The microvoids result from non-uniform microarchitecture of the fiber preform, and they are transported through the gap between fiber tows during infiltration of the resin. In this study, a real-time simulation/control method was proposed to actively control the formation and the transport of air voids during the mold filling. The flow equations were solved in real time to predict the change of the flow front shape. The flow front was detected by optical sensors and the control actions were taken based on the sensor signals. Through this automated simulation/control scheme, a real-time control of resin flow could effectively avoid the dry spots and minimize the formation of microvoids by modulating the injection pressure.     

The effect of fabric and fiber tow shear on dual scale flow and fiber bundle saturation during liquid molding of textile composites

In Liquid Composite Molding (LCM) processes, a fibrous reinforcement preform is placed or draped over a mold surface, the mold is closed and a resin is either injected under pressure or infused under vacuum to cover all the spaces in between the fibers of the preform to create a composite part. LCM is used in a variety of manufacturing applications, from the aerospace to the medical industries. In this manufacturing process, the properties of the fibrous reinforcement inside the closed mold is of great concern. Preform structure, volume fraction, and permeability all influence the processing characteristics and final part integrity. When preform fabrics are draped over a mold surface, the geometry and characteristics of both the bulk fabric and fiber tow bundles change as the fabric shears to conform to the mold curvature. Numerical simulations can predict resin flow in dual scale fabrics in which one can separately track the filling of the fiber tows in addition to flow of resin within the bulk fabric. The effect of the deformation of the bulk fabric due to draping over the tool surface has been previously addressed by accounting for the change in fiber volume fraction and permeability during the filling of a mold. In this work, we investigate the effect of shearing of the fiber tows in addition to bulk deformation during the dual scale filling. We model the influence of change in fiber tow characteristics due to draping and deformation on mold filling and compare it with the results when the fiber tow deformation effect is ignored. Model experiments are designed and conducted with a dual scale fabric to characterize the change in permeability of fiber tow with deformation angle. Simulations which account for dual scale shear demonstrate that the tow saturation rate is affected, requiring longer fill times, or higher pressures to completely saturate fiber tows in areas of a mold with high local shear. This should prove useful in design of components for applications in which it is imperative to ensure that there are no unfilled fiber tows in the final fabricated component.     

Fabric structure and mold curvature effects on preform permeability and mold filling in the RTM process. Part II. Predictions and comparisons with experiments

In Part I, an experimental study was completed in a series of five molds, each having corners of different radii (from 0.06 to 8.0 in.). The primary goal of this work has been to determine whether corners in LCM molds significantly affect the filling process, by altering the structure of the preform locally in such regions. Consistent trends were found for each series of experiments completed in the mold, for which the same preform type, and number of layers were used. For constant flow rate injection, the required injection pressures to fill the two molds with tighter radii was significantly increased as compared to the other molds. Composite parts were manufactured in these molds, measurements made on these parts revealing design flaws, the cavity thicknesses not being equal in all sections of the molds. Several numerical simulations are presented in this paper, the goal being to separate any effects due to the varying thicknesses from effects due to the corners present. Careful simulations have been completed, taking into account the actual thicknesses in each mold, and the resulting preform volume fraction. Experimentally measured permeability data was employed, and the predicted injection pressures match very well for all four molds studied, seeming to indicate that the corners present have not affected the filling. A simple model for preform compression in corners has been developed, which predicts local permeability modifications due to in-plane compression of the fabric layers. These predictions have been employed in conjunction with an existing tool to analytically predict the permeability components of a preform in a flat cavity. This code requires from the user a geometrical model of a preform unit cell, this data being measured from samples cut from the parts manufactured. The resulting predictions for injection pressure are good for an entirely predictive approach, underpredicting experimental values by only 30–60%. Sensitivity analyses have demonstrated the strong relationship between permeability and the details of the preform cell. Two numerical studies were completed to determine how sensitive the injection pressure curves are to reduced permeabilities in the corner regions. For the two injection schemes having two different gate locations, pressures were not significantly affected, while the permeabilities in this region were reduced up to 100 times. Though the molds used were not ideal for isolating effects on mold filling due to corner radii, the evidence presented does not show the existence of any strong behavior related to mold radii.     

基于混合网格方法的VARTM工艺充模仿真与实验验证

针对VARTM工艺的特点,建立了充模过程树脂流动和预成型体变形行为数学模型。提出了基于混合网格方法的VARTM充模仿真算法,在该算法中,模具型腔几何模型进行二维或三维网格划分,在每个真空袋表面单元上增加一个一维附属单元,用于在仿真过程中实时地吸收或挤出因真空袋变形而产生的局部树脂体积变化,形成混合网格仿真模型;求解过程中,对树脂流动和预成型体变形分别进行求解后,基于上述混合网格模型进行两者耦合操作,实现了仿真精度和速度的统一。搭建了VARTM充模实验平台,进行了一维充模实验,通过仿真结果与实验测量结果对比,验证了本文算法的正确性。最后,通过三维仿真算例,验证了算法对三维复杂结构和顺序浇口策略仿真的可行性。     

Approximate numerical method for prediction of temperature distribution in flow through narrow gaps containing porous media

The paper examines an approximate solution to the problem of temperature distribution in the flow through a narrow gap filled with porous medium. This problem is relevant for composite material manufacturing processes, such as the resin transfer molding and its derivatives. In such processes, cold resin is sometimes injected into a closed heated cavity containing a network of stationary fibers. The cavity can be long and wide, but it is usually only a few millimeters thick. The approximate solution is based on an estimated through-the thickness temperature profile that allows one to reduce the complexity from fully three-dimensional to two-dimensional and significantly improve the computational efficiency. The implementation of this solution within an existing simulation code for non-isothermal filling of closed mold is briefly described. Finally, several examples are presented to illustrate the accuracy of the proposed method to predict the midplane temperature of the resin as it impregnates the preform to fill the mold cavity.The support for this work has been provided by the Office of Naval Research (ONR) under Grant #N00014-97-C-0415 for the Advanced Materials Intelligent Processing Center at the University of Delaware.     

Development of adaptive injection flow rate and pressure control algorithms for resin transfer molding

To prevent dry spot formation during fabrication of composite parts by Resin Transfer Molding (RTM), a control interface and four different adaptive control algorithms have been developed and tested with numerical simulations. The interface is capable of controlling the flow pattern of resin as it fills a mold containing a preform of fiber reinforcement, provided that the mold is equipped with multiple inlet gates, a single vent and a spinal sensor system that continuously feeds the interface with the resin flow front locations along the spine lines connecting the inlet gates to the vent. Four different adaptive control algorithms targeting on injection flow rate control, injection pressure control, linearly-corrected pressure control, and the combined flow rate and linearly-corrected pressure control have been proposed and incorporated with the control interface. To provide desirable controllability of the filling process and effective utilization of the resin dispensing equipment, the final formulations were optimized by means of numerical simulations of a rectangular RTM part containing different permeability distributions. The results were compared to investigate the strengths and weaknesses of the spinal adaptive control algorithms in terms of dry spot size, filling speed, and the minimum responding speed of injection pump. Finally, a complex geometry case study was conducted to validate and highlight the spinal adaptive control algorithms’ capability in handling flow disturbance for a complex RTM mold filling process which involves irregular mold geometry, multiple inserts, significant permeability and racetracking variations, and non-straight spinal sensors.     

Application of the level set method to the simulation of resin transfer molding

A new methodology is presented to simulate mold filling in resin transfer molding (RTM) using a combination of the level set and boundary element methods (BEMs). RTM is a composite manufacturing process where a liquid resin is injected in a closed rigid mold containing a dry fibrous reinforcement. Process simulation is motivated by the importance of tracking accurately the motion of the flow front during the mold filling stage. The BEM solves the equation governing the resin flow and the level set method is implemented to track the resin front in the mold. This formulation opens up new opportunities to improve RTM flow simulations and optimize injection molds. The present paper focuses on isothermal resin flow in undeformable porous medium. The implementation of the numerical algorithm is described and several examples of two-dimensional filling with single or multiple injection gates are presented. The robustness of the coupling and the ability to predict accurately the position of the front by this new model are discussed. It is also shown how dry spot formation can be tracked precisely during the simulation and how a generalization of this approach allows predicting resin flow across obstacles.     

A Two-layer Model for the Simulation of the VARTM Process with Resin Distribution Layer

Vacuum assisted resin transfer molding (VARTM) is one of the important processes to fabricate high performance composites. In this process, resin is drawn into the mold to impregnate the fiber reinforcement to a form composite. A resin distribution layer with high permeability was often introduced on top of the fiber reinforcement to accelerate the filling speed. Due to the difference of the flow resistance in the resin distribution layer and the reinforcement as well as the resulting through thickness transverse flow, the filling flow field is intrinsically three-dimensional. This study developed a two-layer model with two-dimensional formulation to simulate the filling flow of the VARTM process with a resin distribution layer. Two-dimensional flow was considered in each layer and a transverse flow in the thickness direction was estimated between the two layers. Thermal analysis including the transverse convection was also performed to better simulate the temperature distribution.     

An optically-based inverse method to measure in-plane permeability fields of fibrous reinforcements

Structural composite manufacturing relying on Liquid Composite Molding technologies is strongly affected by local variability of the fibrous reinforcement. Optical techniques using light transmission are used and allow field measurements of areal weight (and fibre volume fraction) of glass fibre reinforcement. The coupling of obtained areal weight mappings along with injection flow fronts is used to extract in-plane permeability fields. The current work presents results with a focus on glass random mats, but the method can be adapted to any glass fibrous medium. A study of convergence and error due to discretization is performed. Also the influence of the stacking of fibrous layers on the preform variability is analyzed. The major advantage of the proposed technique is a relatively fast acquisition of statistical data on reinforcement variability, which can be later utilized in stochastic based process simulations.     

Numerical prediction of saturation in dual scale fibrous reinforcements during Liquid Composite Molding

This paper presents a fractional flow model based on two-phase flow, resin and air, through a porous medium to simulate numerically Liquid Composites Molding (LCM) processes. It allows predicting the formation, transport and compression of voids in the modeling of LCM. The equations are derived by combining Darcy’s law and mass conservation for each phase (resin/air). In the model, the relative permeability and capillary pressure depend on saturation. The resin is incompressible and the air slightly compressible. Introducing some simplifications, the fractional flow model consists of a saturation equation coupled with a pressure/velocity equation including the effects of air solubility and compressibility. The introduction of air compressibility in the pressure equation allows for the numerical prediction of the experimental behavior at low constant resin injection flow rate. A good agreement was obtained between the numerical prediction of saturation in a glass fiber reinforcement and the experimental observations during the filling of a test mold by Resin Transfer Molding (RTM).