Multiscale modeling of unsaturated flow in dual-scale fiber preforms of liquid composite molding III: reactive flows

The woven, stitched or braided fabrics used in liquid composite molding (LCM) display partial saturation behind moving flow-front in an LCM mold which is caused by delayed impregnation of fiber tows. In this part 3 of the present series of three papers, a novel multiscale approach proposed in parts 1 and 2 [1] and [2] is adapted for modeling the unsaturated flow observed in the dual-scale fabrics of LCM under non-isothermal, reactive conditions. The volume-averaged species or resin cure equation, in conjunction with volume-averaged mass, momentum and energy (temperature) equations, is employed to model the reactive resin flow in the inter-tow (gap) and intra-tow (tow) regions with coupling expressed through several sink and source terms in the governing equations. A coarse global-mesh is used to solve the global (gap) flow over the entire domain, and a fine local mesh in form of the unit-cell of periodic fabrics is employed to solve the local (tow) flows. The multiscale algorithm based on the hierarchical computational grids is then extended to solve the dual-scale flow under reactive conditions. The simulation is compared with a two-color experiment and a previously published two-layer model. Significant differences between the temperatures and cures of the gap and tow regions of the dual-scale porous medium are observed. The ratio of pore volumes in the tow and gap regions, the effective thermal conductivity in the tows, and the reaction rate are identified as the important parameters for temperature and cure distributions in the gap and tow regions.     

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

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

Variations in unsaturated flow with flow direction in resin transfer molding: An experimental investigation

The dual-scale nature of fiber preforms due to the presence of large continuous gaps between fiber tows gives rise to the unsaturated flow in resin transfer molding (RTM) process which is characterized by a droop in the injection pressure history due to the delayed absorption of fiber tows (the ‘sink’ effect). In this study, we experimentally investigate the effect of change in flow direction on the unsaturated flow in three anisotropic dual-scale fiber mats. A series of 1-D mold-filling experiments involving a constant flow rate were conducted for a unidirectional woven fiber-mat, a biaxial stitched mat, and a triaxial stitched fiber-mat along with a reference single-scale random mat. In the case of the unidirectional mats, the droop in the inlet-pressure history, signifying the strength of the sink effect, is found to be strongest for flow along the micro-channels aligned with fiber tows. The droop, and hence the sink effect, is observed to weaken progressively for flow-directions at 45° and 90° to this principal direction. In the case of the biaxial and triaxial mats, the situation is more complex due to the multi-layer construction of such mats: maximum droop is found when mats are oriented at a 45° angle with respect to the fiber-mat coordinate, and it weakens in the 0° and 90° directions. The unsaturated flow effect is also quantified by measuring percentage deviation in the area under the experimental curve from that of the predicted curve. A clear correlation between the droop (through the percentage deviation) and the permeability along a flow direction in the unidirectional mats is observable, though such a relationship eludes the triaxial mat. The effect of unsaturated flow on liquid-front progress during the 1-D experiment was also studied. In contrast to the reference single-scale random mat where the observed front progress closely follow the prediction based on the single-scale physics, a small difference was observed between the observed and predicted front progress for the three dual-scale mats considered. However the difference was too small to yield any significant correlation with the flow direction.     

Multiscale modeling of unsaturated flow in dual-scale fiber preforms of liquid composite molding I: Isothermal flows

Mold-filling simulation of unsaturated flows in LCM is important for optimizing mold design quickly and cost-effectively in the virtual space. For the first time, a true multiscale approach is developed for simulating the unsaturated flow under isothermal conditions in the dual-scale fiber-mats of RTM. To solve the coupled macro-micro equation-set, a coarse global mesh is used to solve the global flow equations over the entire domain while fine local meshes in form of the periodic unit-cells of fabrics are employed to solve the local tow-impregnation process. A multiscale algorithm based on hierarchical computational grids has been proposed to simulate the unsaturated flow in the dual-scale fiber mats under isothermal conditions. The predictions are compared with measurements for a 1-D flow experiment which indicates that the proposed approach can be used to simulate the unsaturated flow accurately through dual-scale fiber mats in LCM without the use of any fitting parameters.     

Modelling microscopic flow in woven fabric reinforcements and its application in dual-scale resin infusion modelling

A physical unit cell impregnation model is proposed for the micro-scale flow in plain woven reinforcements. The modelling results show a characteristic relationship between tow impregnation speed, the surrounding local macro-scale resin pressure and the tow saturation within the unit cell. This relationship has been formulated into a mathematical algorithm which can be directly incorporated into a continuum dual-scale model to predict the ‘sink’ term. The results using the dual-scale model show a sharp resin front in inter-tow-pore spaces and a partially saturated front region in intra-tow-pore spaces. This demonstrates that the impregnation of fibre tows lags behind the resin front in the macro pore spaces. The modelling results are in agreement with two reported experimental observations. It has been shown that the unsaturated region at the flow front could increase or have a fixed length under different circumstances. These differences are due to the variation in tow impregnation speed (or the time required for the tow to become fully impregnated), the weave architecture and the nesting and packing of plies. The modelling results have also demonstrated the drooping of the inlet pressure when flow is carried out under constant injection rates. The implementation of the algorithm into a dual-scale model shows coherence with a single-scale unsaturated model, but demonstrates an advantage in flexibility, precision and convenience in application.     

Multiscale modeling of unsaturated flow of dual-scale fiber preform in liquid composite molding II: Non-isothermal flows

A novel multiscale approach is developed for modeling non-isothermal flows under unsaturated conditions in the dual-scale fabrics of liquid composite molding (LCM). The flow and temperature governing equations at the global or gap or inter-tow (∼m) level and the local or intra-tow (∼mm) levels are based on a previous dual-scale volume averaging method. To solve the coupled equations at two length-scales, a coarse global mesh is used to solve the global flow over the entire domain, and a fine local mesh in form of the unit-cell of periodic fabrics is employed to solve the local tow-impregnation process. (The latter is used to compute sink terms required for solving the former.) A multiscale algorithm based on the hierarchical computational grids is then proposed to solve the dual-scale flow under non-isothermal (but non-reactive) conditions. To test the proposed multiscale model, we first carry out a validation study in which the temperature histories predicted by the multiscale method are compared with experimental data available in a publication for a simple 1-D flow. Despite the lack of information about various model parameters, a reasonably good comparison with the experimental results is achieved. Then, the non-isothermal flow through a simple 1-D flow domain is carried out and the predictions of the multiscale simulation are compared with those of a previously published two-layer model. The multiscale predictions are found to be very similar to the two-layer predictions. A significant difference between the gap and tow temperatures is observed. The ratio of pore volumes in the tow and gap regions, thermal conductivity of the tows, and fiber types are identified as the important parameters for temperature distributions in the gap and tow regions. A further comparison with the single-scale flow simulation highlights significant differences between the conventional single-scale and the proposed dual-scale modeling approaches.     

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.     

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

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

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

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

Longitudinal permeability determination of dual-scale fibrous materials

In this work, the longitudinal permeability of squarely packed dual-scale fiber preforms is studied theoretically. These fiber preforms are composed of aligned porous tows and the tows are tightly packed. The effective permeability is calculated as a parallel-like network of intra-tow permeability and inter-tow permeability, which are quantified by Darcy’s law and the inscribed radius between tows, respectively. The jump velocity at the interface between inter-tow fluids and porous tows is considered, as derived by substituting Beavers and Joseph’s correlation into Brinkman’s equation. We further examine the effects of intra-tow permeability on the effective permeability of the fibrous system with three interface conditions: (1) interface velocity = 0, (2) interface velocity = mean intra-tow velocity, and (3) interface velocity = jump velocity. The jump-velocity-based model is found to be closest to numerical data. The influence of the fiber volume fraction of tows on the effective permeability is also analyzed.     

Correlation of void distribution to VARTM manufacturing techniques

Vacuum assisted resin transfer molding (VARTM) is a liquid composite molding (LCM) technique used to manufacture large scale composite structures. Fiber preforms are placed on a tool surface and covered by a flow enhancement layer and a plastic bag. A vacuum is drawn on the system to infuse the resin. When the resin does not fully saturate the empty regions in between the fibers, voids are created. The fiber tows in woven and stitched preforms have a much lower permeability as compared to the bulk permeability of the fabric. Hence, fiber tows saturate with resin later than the pores between fiber tows and are more prone to voids.This study explores the impact of extended resin bleeding time and additional flow resistance at the vent on the void content within fiber tows both experimentally and by numerical simulation. Samples were machined from each of the manufactured panels and analyzed using image analysis techniques to obtain a relative void content. The experimental results were compared to results obtained by numerical simulation.The experimental void distribution showed that if resin is not allowed to bleed or if no external resistance is attached at the vent, the void content over the length of the part is not uniform. All void levels reduced when resistance was added or bleeding was allowed. The discrepancy between experimental and numerical results was addressed by including deformable distribution media in numerical model to capture the continuation of resin flow after the injection gate is closed.     

Simulation of void formation in liquid composite molding processes

Air entrapment within and between fiber tows during preform permeation in liquid composite molding (LCM) processes leads to undesirable quality in the resulting composite material with defects such as discontinuous material properties, failure zones, and visual flaws. Essential to designing processing conditions for void-free filling is the development of an accurate prediction of local air entrapment locations as the resin permeates the preform. To this end, the study presents a numerical simulation of the infiltrating dual-scale resin flow through the actual architecture of plain weave fibrous preforms accounting for the capillary effects within the fiber bundles. The numerical simulations consider two-dimensional cross sections and full three-dimensional representations of the preform to investigate the relative size and location of entrapped voids for a wide range of flow, preform geometry, and resin material properties. Based on the studies, a generalized paradigm is presented for predicting the void content as a function of the Capillary and Reynolds numbers governing the materials and processing. Optimum conditions for minimizing air entrapment during processing are also presented and discussed.     

Prediction of optimal flow front velocity to minimize void formation in dual scale fibrous reinforcements

Liquid Composite Molding (LCM) is an increasingly used class of processes to manufacture high performance composites. Engineering fabrics commonly used in LCM generally have a dual scale architecture in terms of porosity: microscopic pores exist between the filaments in the fiber tows, while macroscopic pores appear between the tows. Capillary flows in fiber tows play a major role on the quality of composites made by resin injection through fibrous reinforcements. This paper reports on an investigation on fabric imbibition characterization and subsequent evaluation of the optimal flow front velocity during resin injection through fibrous reinforcements. The goal is to devise more robust LCM processes and improve part quality. In order to evaluate a priori the injection conditions that minimize void formation, an impregnation model is developed based on imbibition characterization. This approach allows predicting the optimal front velocity without having to model complex dual scale flows through fibrous reinforcements and without performing expensive and time-consuming fabrication tests. After a summary of previous imbibition results obtained with a probe fluid, the optimal modified capillary numbers are computed by the new predictive model and the values are compared with results reported in the literature on void formation in LCM processes. Afterwards, capillary rise measurements are carried out with four infiltration fluids in order to evaluate the range of optimal flow front velocity that minimizes void formation. This characterization is implemented with vinyl ester resin, epoxy anhydride resin, styrene and anhydride. Finally, the optimal flow front velocity is evaluated for several fabric configurations.     

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.     

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

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

Modelling the effect of material properties and process parameters on tow impregnation in out-of-autoclave prepregs

Out-of-autoclave prepregs based on woven fabrics initially consist of dry tows and resin-rich areas. The tows allow air evacuation in the initial stages of processing and are subsequently infiltrated by surrounding matrix. The following study analyzes the relationship between material properties, process parameters and tow impregnation for three OOA prepregs. First, a representative model for tow impregnation is developed. Then, the model parameters are determined and the model predictions are correlated to impregnation data measured by X-ray microtomography. Finally, the model is used in a parametric study to investigate the effect of fibre architecture, cure cycle temperature and resin initial degree of cure on tow impregnation rate and to predict the possible occurrence of flow-induced micro-voids.     

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.     

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.     

Thermoplastic Polyurethane Reinforced with Continuous Carbon Fiber Tows: Manufacturing Technique and Fabric Property Evaluation

Co-extrusion technique is used to combine carbon fiber tows and thermoplastic polyurethane (TPU) blends into composite carbon fibers (CCFs). The TPU blends serve as the shell and the carbon fiber tows serve as the core. CCFs are fabricated into woven composites and knitted composites with electromagnetic interference shielding effectiveness (EMI SE). The effects of the content of carbon fiber tows on different fabrics are examined. The failure mechanism of CCFs is improved as the fiber slippage along the axial direction replaces the fiber fracture. Fabrics made of CCFs are soft. The laminated woven composites have an average EMI SE of -27.20 dB at 0–3000 MHz and -44.99 dB at 2450 MHz, suggesting a good EMI attenuation capability. This study proposes the use of co-extrusion technique followed coating and fabrication of carbon fibers. This brand-new concept would be vital in the development of carbon fiber composites.

     

Influence of Scale Specific Features on the Progressive Damage of Woven Ceramic Matrix Composites (CMCs)

It is well known that failure of a material is a locally driven event. In the case of ceramic matrix composites (CMCs), significant variations in the microstructure of the composite exist and their significance on both deformation and life response need to be assessed. Examples of these variations include changes in the fiber tow shape, tow shifting/nesting and voids within and between tows. In the present work, the influence of many of these scale specific architectural features of woven ceramic composite are examined stochastically at both the macroscale (woven repeating unit cell (RUC)) and structural scale (idealized using multiple RUCs). The recently developed MultiScale Generalized Method of Cells (MSGMC) methodology is used to determine the overall deformation response, proportional elastic limit (first matrix cracking), and failure under tensile loading conditions and associated probability distribution functions. Prior results showed that the most critical architectural parameter to account for is weave void shape and content with other parameters being less in severity. Current results show that statistically only the post-elastic limit region (secondary hardening modulus and ultimate tensile strength) is impacted by local uncertainties both at the macro and structural level.