Papers by Author: Fabrice Schmidt

Abstract: Carbon Fiber Reinforced Polymers (CFRPs) are essential to the aerospace industry, offering superior strength-to-weight ratios. Currently, the manufacturing of primary structures via standard autoclave curing is a robust, mastered process that successfully minimizes defects, keeping porosity levels below critical thresholds (typically < 1 %). Consequently, porosity is generally not considered as an issue in standard, optimized production lines.However, this stability may be affected by emerging industrial paradigms aimed at increasing production rates and reducing costs. The shift toward accelerated manufacturing – characterized by rapid heating rates, shortened cure cycles and by new manufacturing processes – and the introduction of complex material architectures risk re-introducing significant porosity. In parallel, there is currently no numerical model capable of accurately predicting porosity formation and evolution under these complex conditions. Existing simulation approaches are typically macroscopic and rely on homogenized porous media assumptions, failing to capture the essential micro-scale interactions between bubbles and fibres.To address this gap, this study presents an extended, custom multi-physics Computational Fluid Dynamics (CFD) solver built upon an existing OpenFOAM framework. The goal is to provide the first predictive tool for void evolution within realistic microstructures. The numerical framework couples a two-phase compressible flow model with the complete thermo-chemo-rheological physics of thermoset curing.The solver is applied to 2D Representative Volume Elements (RVEs) of a prepreg ply. Simulations of a standard autoclave cycle demonstrated the solver's ability to capture micro-scale dynamics, showing how voids are compressed and transported during the resin viscosity drop before being frozen at gelation. A parametric study comparing 3-bars and 7-bars pressures confirmed the model's physical ability in predicting void volume reduction.While currently focused on mechanical compression, the tool is designed to support the development of future manufacturing cycles. Future work will incorporate moisture diffusion physics and includes experimental validation via X-ray micro-tomography and in-situ synchrotron monitoring.

Abstract: Nowadays, injection stretch blow molding (ISBM) process represents the most employed technology to produce plastic bottles. An important step of this process is the heat conditioning stage which is performed within infrared ovens by the use of powerful halogen lamps. Homogenizing the temperature distribution along and inside the preform at the end of this conditioning stage is one of the key parameters to determine the final quality of the bottle (thickness, mechanical properties, transparency...). In this research work, a numerical software has been developed to simulate recycled PET (rPET) preforms infrared heating inside the industrial ISBM ovens, where both rotation and translation of the preform across different heating modules occur. In addition, the presence of a fan system involving a forced convective condition inside the ovens is also considered using a Computational Fluid Dynamics (CFD) approach instead of using a conventional heat transfer coefficient.

Abstract: In previous studies [1, , we have presented a detailed formulation of a macroscopic analytical model of the optical propagation of laser beams in the case of unidirectional thermoplastic composites materials. This analytical model presented a first step which concerns the estimation of the laser beam intensity at the welding interface. It describes the laser light path in scattering transparent composites (first component) by introducing light scattering ratio and scattering standard deviation. The absorption was assumed to be negligible in regard to the scattering effect. In this current paper, in order to describe completely the laser welding process in composite materials, we introduce the absorption phenomenon in the model, in the absorbing material (second component), in order to determine the radiative heat source generated at the welding interface. Finally, we will be able to perform a three dimensional temperature field calculation using a commercial FEM software. In laser welding process, the temperature distribution inside the irradiated materials is essential in order to optimize the process. Experimental measurements will be performed in order to valid the analytical model.

Abstract: Composite stamping is a two steps process that includes an infrared heating oven in order to melt the composite sheets before forming. This study deals with the numerical simulation of the heating step of the process. The numerical model has been validated using three woven glass and carbon / PA6.6 composites provided by Solvay Rhodia. This type of simulation consists in solving the heat equation with a radiative flux that characterizes the interaction of the material with the IR heating. The model thus considers the IR properties of the material (emission and reflexion). Considering a homogeneous composite, the optical and thermal properties of sheets have been first measured. The material’s emissivity has been measured using a FTIR spectrometer from the reflective and transmitive spectra, by using the Kirchoff law and considering a Lambertian material. Three complementary measurement techniques were used to determine the thermal properties of the composites. Differential Scanning Calorimetry (DSC) measurements have been performed to identify the heat capacity of the composites. On another hand, a hot disc system (measurements performed at the LTN, France) has been used in transient conditions to determine the heat capacity and the thermal conductivity of the composites is all three directions. Finally, the in-plane thermal matrix of conductivity has also been measured by thermography by using an inverse method. The simulation of composites heating has been performed with Comsol MultiphysicsTM and the simulation procedure was validated by comparison with experimental results. The simulated IR oven is composed of 9 IR emittors provided by Toshiba Lighting Company that emit mainly in short IR wavelength (0.75-2µm). The emission properties of the tungsten filament were implemented in order to simulate the IR heating. Free convective heat transfer was also taken into account in the oven. In order to validate the model, an experimental set-up was instrumented with a calibrated IR pyrometer that measured the back side of the heated composite sheets. The experimental results confirm a low thermal gradient through composite thickness, in particular for carbon-reinforced composite. This result is consistent with the low Biot number of the composites. Moreover, experimental and simulated temperatures are in good agreement with an error lower than 15% in the entire heating stage from room to melt temperature.

Abstract: Epoxy resins have several applications in the aerospace and automobile industry. Because of their good adhesive properties, superior mechanical, chemical and thermal properties, and resistance to fatigue and micro cracking, they produce high performance composites. In the technology presented here, the composite is cured in an IR oven which includes halogen lamps. The liquid resin infusion (LRI) process is used to manufacture the composite, whereby liquid resin is infused through a fiber reinforcement previously laid up in a one-sided mold. These epoxy resins release an exothermic heat flux during the curing process, which can possibly cause an excessive temperature in the thickness. Consequently, for the production of high performance composites, it is necessary to know the thermal behavior of the composite during curing. Therefore, IR interactions with the graphite/epoxy system were modeled as a surface radiation transport. In our work, we have studied IR interactions with the composite, which is placed in an IR oven. Using an IR spectrometer Bruker Vertex 70 (1-27 μm), we measured radiative properties and determined the fraction of IR rays absorbed by the composite. Since it is necessary to optimize the manufacturing time and costs and to determine the performance of these composites, the purpose of this study is to model the IR curing of a composite part (carbon fiber reinforced epoxy matrix) in the infrared oven. The work consists in two parts. In the first part, a FE thermal model based on radiosity method was developed, for the prediction of the infrared incident heat flux on the top surface of the composite during the curing process. This model was validated using a reference solution based on ray tracing algorithms developed in Matlab® (In-lab software called Rayheat based on ray tracing algorithms is used to compute the radiative heat flux that impacts the composite). Through the FE thermal model, an optimization study on the percentage power of each infrared heater is performed in order to optimize the incident IR heat flux uniformity on the composite. This optimization is performed using the Matlab® optimization algorithms based on Sequential Quadratic Programming method. In a second part, the optimized parameters set is used in a three-dimensional numerical model which is developed in the finite element commercial software Comsol Multiphysics ™, where the heat balance equation is coupled with the cure kinetic model of the resin. This numerical model allows calculation of the temperature distribution in the composite during curing, which is a key parameter that affects its mechanical properties. In this model, we can predict also the evolution of the degree of cure as function of time. Experimental measurements were used to validate simulations of the whole infrared composite curing process. Keywords: Curing composite, infrared oven, Radiation, Optimization, Epoxy resin, Carbon fibers.

Abstract: In the industrial process, the moisture of the clay sheet obtained by extrusion and pressed to form a tile varies in time. It depends on the nature and the mixing of the raw materials during the production. In order to model and undersand the influence of the moisture on the pressing step, it is necessary to determine the parameters of the rheological and tribological laws. A study of the rheological behaviour, based on free compression tests on cylinder samples, allowed to use an elasto-visco-plastic behaviour for the extruded clay paste. The different constitutive parameters were estimated by an inverse analysis based on the experimental force/displacement curves. The identification was performed with the optimisation algorithm implemented in the commercial software Forge® 2009. The influence of the water content in the paste on the rheological parameters was identified and fitted using linear models. The friction factor was measured from tests on a rectilign tribometer. To understand the influence of the moisture, we simulated a compression test, using Forge® involving the shaping of a tile lug. This geometry is representative of the state of stress during the pressing of the tile, in an area currently sujected to defects. The numerical model show that an increase of eighteen percent of the moisture allows to decrease by half the pressing force.

Abstract: The main purpose of this study is to cure a 3D geometry composite part (carbon fiber reinforced epoxy matrix) using an infrared oven. The work consists of two parts. In the first part, a FE thermal model was developed, for the prediction of the infrared incident heat flux on the top surface of the composite during the curing process. This model was validated using a reference solution based on ray tracing algorithms developed in Matlab®. Through the FE thermal model, an optimization study on the percentage power of each infrared heater is performed in order to optimize the incident IR heat flux uniformity on the composite. This optimization is performed using the Matlab® optimization algorithms based on Sequential Quadratic Programming and dynamically linked with the FE software COMSOL Multiphysics®. In a second part, the optimized parameters set is used in a model developed for the thermo-kinetic simulations of the composite IR curing process and the predictions of the degree of cure and temperature distribution in the composite part during the curing process.

Abstract: Liquid Composite Molding (LCM) is a popular manufacturing process used in many industries. In Resin Transfer Molding (RTM), the liquid resin flows through the fibrous preform placed in a mold. Numerical simulation of the filling stage is a useful tool in mold design. In this paper the implemented method is based on coupling a Boundary Element Method (BEM) with a Level Set tracking. The present contribution is a two-dimensional approach, decoupled from kinetics, thermal analysis and reinforcement deformation occurring during the flow. Applications are presented and tested, including a flow close to industrial conditions.

Abstract: The use of composite materials in large structures has increased in recent years. Aircraft industry has recently begun to investigate the field of Liquid Composite Molding (LCM) through research programs because of its ability to produce large parts at low cost. The present paper focuses on modeling a 3D radial impregnation through an anisotropic fibrous perform. As a preliminary work, it is assumed an isothermal flow and no hydro-mechanical coupling. Governing equations are Darcy’s law and mass conservation. Simulation is performed combining Boundary Element Method (BEM) with a lagrangian moving mesh method. An analytical solution is developed to assess the numerical model.