# Optimization of precharge placement in sheet molding compound process

Fariba Ebrahimian, Sebastian Rodriguez, Daniele Di Lorenzo, Francisco Chinesta Soria

► To cite this version:

Fariba Ebrahimian, Sebastian Rodriguez, Daniele Di Lorenzo, Francisco Chinesta Soria. Optimization of precharge placement in sheet molding compound process. International Journal of Material Forming, 2024, 17 (4), pp.34. 10.1007/s12289-024-01836-w. http://dx.doi.org/10.1007/s12289-024-01836-w

HAL Id: hal-04768694

https://hal.science/hal-04768694v1

Submitted on 6 Nov 2024

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
# Optimization of precharge placement in sheet molding compound process

Fariba Ebrahimian<fnref n="1" />, ID · Sebastian Rodriguez<fnref n="1" /> · Daniele Di Lorenzo<fnref n="1" />,<fnref n="2" /> · Francisco Chinesta<fnref n="1" />,<fnref n="2" />

Received: 6 March 2024 / Accepted: 20 May 2024 / Published online: 1 June 2024
© The Author(s) 2024

## Abstract

This study aims to provide precise predictions for the compression of reinforced polymers during the sheet Molding Compound (SMC) process, ensuring the attainment of a predefined structure while preventing material overflow during the process. The primary challenge revolves around identifying the optimal initial shape to prevent material rebound during the process. To confront this issue, a numerical model is utilized, faithfully simulating the SMC process and forming the foundation for our investigations. Furthermore, to optimize the pre-fill stage, a surrogate model is proposed to enhance modeling efficiency, and then an inverse analysis method is applied. This approach of minimizing material rebound during the SMC process results in a reliable metamodel to predict an initial mass shape accurately and at a low computational cost, thus ensuring the squeezed material fits the mold shape.

**Keywords** Surrogate model · SMC process · Feedforward neural networks · Model Order Reduction(MOR) · Inverse analysis

## Introduction

Over the last two decades, the benefits of utilizing polymer and composite structures have become increasingly evident. These structures offer numerous advantages, including lightweight strength, corrosion resistance, and design flexibility [1]. These features reduce costs and increase value in industries where performance and durability matter. Recognizing the advantages of polymer and composite struc-

tures, various industries have adopted these materials to revolutionize their manufacturing processes and product designs. Among these, the automotive sector stands out as a prime example of innovation in material application. Between the various methods, the sheet molding compound (SMC) process has emerged as the preferred choice for manufacturing composite structures in the automotive industry, primarily due to its exceptional high-volume production capabilities.

In the SMC, polyester resin and glass or carbon fibers are combined and heated to produce a hot substance that is then compressed in a mold [2]. The process is illustrated in Fig. 1.

However, the process has its challenges because critical process parameters, such as curing time, mold closing speed, molding pressure, and precharge specifications (including geometry, placement, and size), significantly influence the final product's quality in the compression molding process [4]. The mechanical properties of the final structure are considerably influenced by the incorporation of glass or carbon fibers within the resin [5]. This influence stems from the alteration of the SMC material's concentration and flow behavior within the mold due to the presence of fibers. Consequently,

All authors contributed equally to this work.

<figure><img src="image_2.png" /></figure>

Fariba Ebrahimian
fariba.ebrahimian@ensam.eu

Sebastian Rodriguez
sebastian.rodriguez_iturra@ensam.eu

Daniele Di Lorenzo
daniele.di.lorenzo@ensam.eu

Francisco Chinesta
francisco.chinesta@ensam.eu

<footnote n="1">¹ PIMM Lab, ENSAM Institute of Technology 151 Boulevard de l'Hôpital, Paris 75013, France</footnote>

<footnote n="2">² ESI Group, 3bis Rue Saarinen, 94528 Rungis CEDEX, France</footnote>
<figure><img src="image_3.png" /><figcaption>Fig. 1 Squeeze flow phenomenon observed in a confined narrow gap [3]</figcaption></figure>

numerous researchers have focused on understanding the SMC process of reinforced polymers [6–8].

Simulating and optimizing the manufacturing process is crucial to produce high-quality polymer and composite structures, especially reinforced polymers. Since the SMC cycle test costs are not cheap, computer simulations are used as a reference [9]. However, in contemporary scientific research, the consensus is that numerical modeling is a resource-intensive approach, mainly when dealing with complex phenomena. Therefore, surrogate models emerge as an attractive solution to speed up optimization processes [10, 11]. Surrogate models can be built employing different techniques, among the most popular one finds Genetic Algorithms(GA) [12], Proper Generalized Decomposition(PGD) [13, 14], sparse-PGD [15], and Neural Networks [16].

Neural networks possess significant advantages in the domain of data analysis and optimization. Their notable adaptability allows them to effectively tackle diverse tasks, including pattern recognition, image classification, and natural language processing. This versatility positions neural networks as a powerful tool for complex data analysis [17]. Additionally, the capability of neural networks to comprehend intricate non-linear relationships make them quite adept at handling complex data patterns and achieving high predictive accuracy across various real-world applications [18–20].

In this paper, a surrogate for the SMC process is proposed based on the use of the POD model reduction method and the use of neural networks. The objective is to optimize the SMC parameters using a rapid simulation model in a reverse engineering context [21]. Therefore, the initial shape of the mass is optimized using this surrogate so that at the end of the process the mass coincides with the edge of the mold, thus minimizing the rebound of the resin in the mold. In this regard, established optimization methodologies such as gradient descent [22], Newton’s method, Levenberg-Marquardt [23, 24] algorithm could be employed.

The paper is organized into several sections. Section 2 introduces the mathematical formulation and assumptions used for simulating compression and material flow. In Section 3, the numerical implementation and methodology introduced in the previous section are presented in detail. Additionally, Section 4 provides an explanation of the surrogate model construction, specifically the Feedforward neural network that is used to optimize the SMC process. Section 5 presents a numerical example to demonstrate the application of the surrogate for compression molding behavior and optimization of precharge. Finally, Section 6 provides conclusions and perspectives.

## Modeling of SMC process

In this section, we outline the mathematical modeling employed that will be used to numerically simulate a Sheet Molding Compound (SMC) process necessary to build a surrogate model. Our focus lies in capturing the evolution of the flow-front for each initial shape during the process. Given the complexity of the process and the need to prioritize optimizing precharge dimensions, it is practical to simplify the SMC process at this stage. Consequently, we have chosen to neglect the effect of fiber orientation in the material flow kinematics. In order to describe this process, it’s essential to identify and explain the key relationships and equations that govern the compression process. These will be detailed in the following section.

## Squeeze flow approximation

In modeling the compression process, we assume they behave as a Newtonian fluid. Also, the motion of fibers does not influence the main flow field of resin. The initial equation for modeling molding processes such as Resin Transfer
Molding (RTM) [25] or Sheet Molding Compound (SMC) [7] is derived from the Navier-Stokes equations. To simplify further the model, we make use of the narrow gaps hypothesis ($\frac{2h}{L} \ll 1$), as illustrated in Fig. 2. Additionally, by assuming the validity of the narrow gaps hypothesis, the lubrication theory becomes a valuable tool for significantly reducing computational complexity in the modeling process [3].

So, the governing equation for flow becomes:

$$\nabla P = \eta \Delta \mathbf{v}, \qquad (1)$$

where $P$ is the pressure, $\mathbf{v} = (u, v, w)$ is the velocity vector and $\eta$ is the dynamic viscosity while $\nabla$ and $\Delta$ represent respectively the gradient and Laplacian operator. In what follows, the velocity gradient in the thickness direction is much larger than in the $x$ and $y$ (in-plain gradients) direction and the component of velocity in $z$ direction ($w$) is neglected with respect to the in-plain component:

$$\left\{ \begin{array}{l} \frac{\partial u}{\partial z} \gg \frac{\partial u}{\partial x}, \quad \frac{\partial u}{\partial z} \gg \frac{\partial u}{\partial y} \\ \frac{\partial v}{\partial z} \gg \frac{\partial v}{\partial x}, \quad \frac{\partial v}{\partial z} \gg \frac{\partial v}{\partial y} \\ w \approx 0 \end{array} \right. . \qquad (2)$$

Therefore, in the Eq. 1, $P = P(x, y)$. By performing the integration with the $z$-variable and enforcing no-slip boundary conditions at both gap walls, we obtain:

$$\mathbf{v}(x, y, z) = \frac{1}{2\eta} \nabla P (z^2 - h^2), \qquad (3)$$

representing the velocity of the flow. The velocity averaged over the gap thickness is then:

$$\tilde{\mathbf{v}}(x, y) = \frac{1}{2h} \int_{-h}^{+h} \mathbf{v}(x, y) dz = -\frac{h^2}{3\eta} \nabla P, \qquad (4)$$

<figure><img src="image_4.png" /><figcaption>Fig. 2 Squeeze flow phenomenon observed in a confined narrow gap</figcaption></figure>

while the flow-rate ($\mathbf{q}(x, y)$) is given by:

$$\mathbf{q}(x, y) = 2h\tilde{\mathbf{v}}(x, y). \qquad (5)$$

Finally, the mass conservation is given by:

$$\nabla \cdot \mathbf{q}(x, y) = U, \qquad (6)$$

where the compression rate $U$ is assumed to be constant. Now, by replacing expression Eq. 5 into Eq. 6, and considering $\tilde{U} = \frac{U}{2h}$, we have:

$$-\frac{h^2}{3\eta} \Delta P = \tilde{U}. \qquad (7)$$

It is important to note that as the fluid is subjected to compression, its thickness reduces progressively with time and can be described by the following relation:

$$h(t) = h_0 - \frac{U \times t}{2}, \qquad (8)$$

where $h_0$ is the initial thickness. Ultimately, Eq. 7 can be reformulated in weak form, as:

$$\frac{h^2}{3\eta} \int_{\Omega_f(t)} \nabla \hat{P} \cdot \nabla P \, d\Omega = \int_{\Omega_f(t)} \hat{P} \tilde{U} \, d\Omega. \qquad (9)$$

where the fluid domain at time $t$ is defined by $\Omega_f(t)$ and $\hat{P}$ is a test function. By discretizing the weak formulation, we can calculate the pressure, and then according to Eq. 4, we have the average velocity for the fluid domain.

## Flow front tracking

To consider the temporal evolution of the fluid boundary ($\Omega_f(t)$) in Eq. 9, two main approaches can be employed, namely, the fixed or modified mesh techniques [26]. Since the modified mesh involves re-meshing at each time step and it can be highly costly [27], to lower the computational burden, we adopted the fixed mesh approach [25, 28, 29].

This approach consists of using an auxiliary function $I(\mathbf{x}, t)$, defined across the entire domain ($\Omega$), to account for the continuous updating of the boundary. In particular, within the fluid domain $\Omega_f(t)$, $I$ is equal to one, while in the empty domain ($\Omega_e(t) = \Omega - \Omega_f(t)$), it is zero:

$$I(\mathbf{x}, t) = \begin{cases} 1 & \text{if } \mathbf{x} \in \Omega_f(t) \\ 0 & \text{if } \mathbf{x} \in \Omega_e(t) \end{cases}. \qquad (10)$$
At time equals 0, I variable is determined by the initial shape of precharge which is defined in advance, this means:

$$I(\mathbf{x}, t = 0) = I_{\text{init}}. \qquad (11)$$

The time evolution of the fluid domain is governed by a scalar and linear advection equation that can be expressed as follows:

$$\frac{\partial I(\mathbf{x}, t)}{\partial t} + \mathbf{v}(\mathbf{x}, t) \cdot \nabla I(\mathbf{x}, t) = 0. \qquad (12)$$

To solve the transport (12), we used the volume of fluid technique along with the Total Variation Diminishing (TVD) method [30].

In order to use the fixed mesh approach we need to extend the Eq. 9 to the whole domain. Thus considering:

* Null pression in the the empty region of the mold:

$$P = 0 \text{ in } \Omega_e(t). \qquad (13)$$

* Flow model in the fluid region:

$$\nabla \cdot \left( -\frac{h^2}{3\eta} \nabla P \right) = \tilde{U} \text{ in } \Omega_f(t). \qquad (14)$$

Using these assumptions, combining the Eqs. 13 and 14, and introducing the I function, the weak formulation defined in the whole domain reads:

$$\int_{\Omega} (I \nabla \hat{P} (\frac{h^2}{3\eta} \nabla P) + (1-I) \hat{P} P) d\Omega = \int_{\Omega} I \hat{P} \tilde{U} d\Omega. \quad (15)$$

A schematic illustration of the model is presented in Fig. 3. To achieve time evolution of SMC process, Eqs. 15 and 12 are solved.

## Numerical implementation

In order to build a workflow that can be easily adapted to simulate a diverse spectrum of shapes, we parameterize the shape of the initially precharge with the help of the Non-Uniform Rational B-Splines (NURBS) method [31]. In particular, a set of m control points is essential to create the curve shape for the initial fluid (precharge). The first and m^th control points are intentionally set to be identical, ensuring the creation of a closed curve. Additionally, we strategically position the second and the (m-1)^th along the same line to eliminate sharp

<figure><img src="image_5.png" /><figcaption>Fig. 3 Mass variations across compression moulding</figcaption></figure>

corners, promoting smoothness and continuity. Therefore, (m-2) independent points are considered to parameterize the shape. Figure 4 shows a shape generated with this methodology.

After fixing the independent points and generating the initial shape configuration, the SMC process has been simulated, employing Finite Element Method (FEM) and Control Volume (CV) techniques in MATLAB. Figure 5 illustrates the approach used to create the numerical simulation model for Sheet Molding Compound (SMC) process.

We consider the elements of the mesh that are inside this curve as full elements ($I(x, 0) = I_{\text{init}} = 1, x \in \Omega_f$). Sub-

<figure><img src="image_6.png" /><figcaption>Fig. 4 Illustration of the shape generated by NURBS</figcaption></figure>
<figure><img src="image_7.png" /><figcaption>Fig. 5 Numerical simulation of compression molding process</figcaption></figure>

sequently, Eq. 15 is solved to determine the pressure and velocity profile along the fluid boundary. Based on the flow pattern evolution, the fluid domain boundary is updated based on Eq. 12. Figure 6 shows the snapshot of the flow front capturing the boundary evolution.

<figure><img src="image_8.png" /><figcaption>Fig. 6 Recovery of the boundary at each time-step</figcaption></figure>

At each time step, the boundary is determined by considering the center of mass of each element falling within the boundary (where the value $I$ ranges between 0 and 1). Then, by interpolating between these points and sampling with a constant number of points ($n$), following a uniform discretization of the angular polar coordinate, ensuring that the number of points used to describe the boundary remains fixed for each time step. This process ensures consistency in the number of points on the boundary, providing stability in the representation of the flow front tracking.

The boundaries describing all the points are ultimately saved in the format of ($X_t, Y_t$):

$$X_t = \begin{pmatrix} x_1 \\ \vdots \\ x_n \end{pmatrix}_t, \quad Y_t = \begin{pmatrix} y_1 \\ \vdots \\ y_n \end{pmatrix}_t, \qquad (16)$$

where $t = \{1, ..., N_t\}$ define the time-step of the simulation, and $n$ is the number of chosen points that discretize the boundary curve. Or in a more compact expression for each time step:

$$R_t = \begin{bmatrix} X_t \\ Y_t \end{bmatrix}. \qquad (17)$$
<figure><img src="image_9.png" /><figcaption>**Fig. 7** Feedforward neural network algorithm</figcaption></figure>

To study different mold shapes, we process (by compressing) various primitive shapes. To create initial shapes, we used a Latin Hypercube Sampling method to produce different control points for the NURBS curves. This sampling method is particularly advantageous because it ensures that the obtained samples accurately reflect the underlying distribution of the data compared to other random sampling methods, and it also reduces the sample size.

Finally, K sampled initial shapes are simulated, and the evolution of the boundary for each case and time step is stored in a matrix.

$$D = \left[ \underbrace{[\mathbf{R}_1^1, \dots, \mathbf{R}_{N_i}^1]}_{\text{1st simulation}}, \dots, \underbrace{[\mathbf{R}_1^K, \dots, \mathbf{R}_{N_i}^K]}_{\text{Kth simulation}} \right]. \quad (18)$$

In the following section, we will utilize the stored boundary to construct the surrogate model.

## Surrogate model

Our main goal is to figure out the best initial shape for each unique mold shape. Therefore, we need to apply an inverse analysis to the SMC model. In order to intelligently process, analyze data, and apply reverse analysis, a data-driven model is proposed. One data-driven model that can significantly reduce computational costs consists of Neural Networks (NN) along with model-order reduction techniques. Neural Networks are capable of learning non-linear relationships and identifying complex patterns in data, which is crucial for solving the SMC problem due to the intricate geometry to which the deformed mass can evolve.

In order to create a Neural Network (NN) model, data of different initial shapes is required. The previous section introduces a numerical model that defines the changing boundaries of the shape during the compression process.

A Feedforward Neural Network (FFNN) is a machine learning tool that links input data to output data. In this case, each vector of Control Points (CP) (that represents the initial shape by (m - 2) independent points) and time connected to

<figure><img src="image_10.png" /><figcaption>**Fig. 8** Updated feedforward neural network algorithm by SVD</figcaption></figure>
<figure><img src="image_11.png" /><figcaption>**Fig. 9** Reverse analysis strategy to identify the precharge state</figcaption></figure>

a specific boundary in that time ($R_t$) as the output (Fig. 7). Through this method, we can develop a neural network model that has the capability to forecast the compressed mass at any given time. Utilizing this approach enables us to have various solutions for every initial mass geometry. Therefore, we can select the appropriate time for compressing the mass based on the mold geometry.

This study trained the FFNN from K samples and their associated flow front evolution stored in matrix D. To establish an analytical correlation between input and output, i hidden layers consisting of j neurons are implied, and the activation function implemented in the hidden layer could depend on the complex patterns and relationships within the data.

<figure><img src="image_12.png" /><figcaption>**Fig. 10** Results at various time steps of a mass compression using the FEM solver</figcaption></figure>

However, training a neural network model can be challenging due to its complex nature and high computational demands. To mitigate these challenges and prevent over fitting, and computing time savings, we can reduce the complexity of neural network models by decreasing the output size. This approach conserves computational resources and leads to more efficient training of FFNN models. Therefore, the Model Order Reduction technique will be applied to the output dataset (D), after which we can train the model using the reduced data.

<figure><img src="image_13.png" /><figcaption>**Fig. 11** In black, the initial shape used to compare the performance of the trained NN with the FEM solver. The mass has been compressed for 10 seconds, and the results from the two approaches are compared. In particular, in blue, the approximation of the neural network, and in red, the solution of the FEM solver</figcaption></figure>
<figure><img src="image_14.png" /><figcaption>Fig. 12 Mold shape chosen to verify the inverse procedure</figcaption></figure>

One of the methods to reduce the dimension is the Singular Value Decomposition (SVD) [32, 33], which is applied to the matrix D obtaining:

$$D = U \Sigma \hat{V}, \quad D \in \mathbb{R}^{2n \times (K \cdot N_t)}. \qquad (19)$$

Where $U \in \mathbb{R}^{2n \times 2n}$ and $V \in \mathbb{R}^{(K \cdot N_t) \times (K \cdot N_t)}$ are unitary and such that $\Sigma \in \mathbb{R}^{2n \times (K \cdot N_t)}$ is a matrix with real, non-negative entries on the diagonal and zeros out the diagonal. By using the SVD and pre-multiplying (19) by $\hat{U}$, we have:

$$\hat{U} D = \Sigma \hat{V} = \lambda. \qquad (20)$$

Where $\lambda$ matrix is essentially a projection of the D matrix onto a lower dimensional to simplify computations and improve analysis. Therefore, the FFNN is trained by a new

<figure><img src="image_15.png" /><figcaption>Fig. 13 Results of the inverse procedure: the shape of the mass identified in as the initial shape to use for the chosen mold (dashed curve) is depicted in violet. In the inverse analysis, the output compression time for this shape is 9.9 seconds</figcaption></figure>

<figure><img src="image_16.png" /><figcaption>Fig. 14 Comparison of the compressed (9.9 seconds) initial shape identified by the inverse procedure with the shape of the mold used as input for the inverse procedure</figcaption></figure>

dataset $\lambda$ that consists of the z first decomposition models. Figure 8 shows a schematic representation of the neural network architecture used in this paper.

## Inverse analysis for precharge reconstruction

Given our objective of determining the optimal initial shape based on the mold shape, we must employ an inverse analysis approach to the neural network model developed to generate boundary profiles in response to various initial shapes. In particular, to minimize the objective function:

$$f = \|\lambda^{\text{mold}} - \lambda^{\text{NN}}\|_2^2. \qquad (21)$$

where $\lambda^{\text{mold}}$ is the arbitrary mold shape, $\lambda^{\text{NN}}$ is generate for a random initial shape parameters and $\| \cdot \|_2$ represents the classic $L_2$ norm. We adopted gradient descent techniques, thus the initial shape parameters $X_k$ will be updated in each iterate accordingly:

$$X_{k+1} = X_k - \alpha \nabla f, \qquad (22)$$

with $\alpha$ being the learning rate. Figure 9 shows the inverse analyses applied to reach the initial shape for particular mold shapes.

## Case study analysis

In this section, we analyze the methodology outlined in the previous section through a numerical case study. Following the rationale presented in Section 4, we construct a surrogate model by simulating 200 molding processes using a FEM solver. In these simulations, the key parameter variations involve the shape of the initial mass, which is represented by nine independent control points of a NURBS. For all shapes
<figure><img src="image_17.png" /><figcaption>Fig. 15 (a) Shape of the mold (dash line) used as input for the inverse analysis, initial shape identified by the inverse analysis (violet line) ; (b) Comparison of the compressed (7.2 seconds) initial shape identified by the inverse procedure with the shape of the mold</figcaption></figure>

considered, a body thickness of 20 mm and a dynamic viscosity of 0.1 Pa/s are imposed. The molding process assumes a compression rate of 1 mm/s and a compression time of 10 seconds. Figure 10 illustrates one of the simulations, showing the boundary shape propagation at different time step during the compression. Throughout all numerical simulations, mass conservation is ensured by monitoring the relative error given by :

$$ Error = 100 \times \frac{|M_{end} - M_{int}|}{M_{int}}, \qquad (23) $$

where $M_{int}$ denotes the mass of the fluid before the compression process starts, and $M_{end}$ represents the compressed mass. The error in mass conservation observed at the end of all simulations performed was no greater than 2.5%, which is acceptable, considering the size of the mesh elements used for the simulation.

The simulated shapes are divided into training and test sets with an 80-20% split ratio. The shapes in the training set are then used to construct the D matrix Eq. 18 and then the SVD method is applied to identify the first 7 modes that best represent these shapes. These modes, combined with time, serve as inputs for a NN designed to simulate the SMC process (Section 4). In this configuration, the NN consists of a single hidden layer containing 50 neurons, with the hyperbolic tangent (tanh) activation function employed. Once the NN is trained, it is tested over the test set. Figure 11 shows the output of the NN of one of the shapes from the test set compressed for 10 seconds wrt to the numerical reference (FEM). The relative error (RE) between the two results (FEM and NN) is calculated as follows:

$$ RE = 100 \times \frac{\|\mathbf{x}_i - \mathbf{x}'_i\|_2}{\|\mathbf{x}_i\|_2}, \qquad (24) $$

where, $\mathbf{x} = (x, y)$ and $\mathbf{x}' = (x', y')$ represent the coordinates of the boundary points of the reference and the NN approximation, respectively. The relative error (24) for all the shapes in the test set was no greater than 1.3%.

After training the Neural Network (NN), we can utilize inverse analysis to deduce the geometry of the initial shape based on the mold's geometry (Section 4.1). Suppose we have a new shape of the mold (Fig. 12) that is outside the train-test set.

Applying the inverse procedure we obtain as initial of mass the one showed in Fig. 13. According to the inversion algorithm, the identified shape will reach the mold's boundary within 9.9 seconds.

To validate the outcome of inverse procedure, the identified initial shape is compressed for 10 seconds and compared with the mold in Fig. 14. The relative error between the compressed mass and mold shapes is 0.06%.

To demonstrate the effectiveness of the presented method, different shapes of the mold have been considered. Figure 15 displays the results of the inverse procedure in the case of a mold with star shape. For this case the relative error between the compressed mass and mold shapes is 0.3%.

## Conclusions and perspectives

The primary focus of our investigation was to pinpoint the optimal initial shape, a critical factor in preventing material return during compression and ensuring the successful attainment of the intended molded shape. To achieve this objective, we employed a Feedforward neural network trained using data from the high fidelity model. In particular, the inverse analysis was conducted within the neural network framework using gradient descent (GD) techniques. The results suggest the ability to predict the optimal shape corresponding to a specific compression duration. As a result, this research can
create a background for investigating the behavior of fibers in order to improve the physical properties of the final product in the SMC process.

**Acknowledgements** We acknowledge the financial support from the European Union Horizon 2020 research and innovation program under the Marie Skłodowska-Curie grant agreement No. 956401 (XS-Meta). We confirm that this paper is not being submitted to any other journal, and there is no conflict of interest between the authors.

<figure><img src="image_18.png" /></figure>

## Declarations

**Conflicts of interest** The authors declare no conflict of interest.

**Open Access** This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.

## References

1. Hsissou R, Seghiri R, Benzekri Z, Hilali M, Rafik M, Elharfi A (2021) Polymer composite materials: a comprehensive review. Compos Struct 262:113640

2. Park C, Lee W (2012) Compression molding in polymer matrix composites. In: Manufacturing techniques for polymer matrix composites (PMCs), pp 47–94. Elsevier. Amsterdam

3. Barasinski A, Chinesta F, Ghnatios C, De Luca P (2023) Process simulation: moulding processes. In: Design and manufacture of structural composites, pp 401–426. Elsevier. Amsterdam

4. Orgéas L, Dumont PJ (2011) Sheet molding compounds. Wiley encyclopedia of composites, pp 1–36

5. Binetruy C, Chinesta F, Keunings R (2015) Flows in Polymers, Reinforced Polymers and Composites: A multi-Scale Approach. Springer, Berlin

6. Perez M, Guevelou S, Abisset-Chavanne E, Chinesta F, Keunings R (2017) From dilute to entangled fibre suspensions involved in the flow of reinforced polymers: a unified framework. J Non-Newtonian Fluid Mech 250:8–17

7. Perez M, Prono D, Ghnatios C, Abisset E, Duval JL, Chinesta F (2020) Advanced modeling and simulation of sheet moulding compound (smc) processes. Int J Mater Form 13:675–685

8. Babu BJ, Waigaonkar S, Rajput A (2011) Experimental investigations and optimization of processibility of sheet moulding compound. Journal of Polymer Engineering

9. Lee L, Fan J, Kim J, Im Y-T (1991) Flow analysis of sheet molding compounds in compression molding. Int Polym Process 6(1):61–72

10. Pfrommer J, Zimmerling C, Liu J, Kärger L, Henning F, Beyerer J (2018) Optimisation of manufacturing process parameters using

deep neural networks as surrogate models. Procedia CIRP 72:426–431

11. Kudela J, Matousek R (2022) Recent advances and applications of surrogate models for finite element method computations: a review. Soft Comput 26(24):13709–13733

12. Kim M-S, Lee WI, Han W-S, Vautrin A (2011) Optimisation of location and dimension of smc precharge in compression moulding process. Comput & Struct 89(15–16):1523–1534

13. Ghnatios C, Masson F, Huerta A, Leygue A, Cueto E, Chinesta F (2012) Proper generalized decomposition based dynamic data-driven control of thermal processes. Comput Methods Appl Mech Eng 213:29–41

14. Chinesta F, Cueto E, Huerta A (2014) Pgd for solving multidimensional and parametric models. Separated representations and pgd-based model reduction: fundamentals and applications. Springer. Berlin, Heidelberg, pp 27–89

15. Ibáñez R, Abisset-Chavanne E, Ammar A, González D, Cueto E, Huerta A, Duval JL, Chinesta Fet al (2018) A multidimensional data-driven sparse identification technique: the sparse proper generalized decomposition. Complexity 2018

16. Heaton J (2018) Ian goodfellow, yoshua bengio, and aaron courville: deep learning: The mit press, 2016, 800 pp, doi: 0262035618. Genet Program Evolvable Mach 19(1-2):305–307

17. Hesthaven JS, Rozza G (2008) Proper generalized decomposition of multiscale problems. SIAM J Numer Anal 46(1):348–377

18. Bengio Y, Courville A, Vincent P (2013) Representation learning: a review and new perspectives. IEEE Trans Pattern Anal Mach Intell 35(8):1798–1828

19. Nielsen MA (2015) Neural Networks and Deep Learning, vol 25. Determination press, San Francisco, CA, USA

20. Hernández Q, Badías A, González D, Chinesta F, Cueto E (2021) Structure-preserving neural networks. J Comput Phys 426:109950

21. Weichert D, Link P, Stoll A, Rüping S, Ihlenfeldt S, Wrobel S (2019) A review of machine learning for the optimization of production processes. J Adv Manuf Technol 104(5–8):1889–1902

22. Kindermann J, Linden A (1990) Inversion of neural networks by gradient descent. Parallel Comput 14(3):277–286

23. Levenberg K (1944) A method for the solution of certain non-linear problems in least squares. Q Appl Math 2(2):164–168

24. Marquardt DW (1963) An algorithm for least-squares estimation of nonlinear parameters. J Soc Ind Appl Math 11(2):431–441

25. Rodriguez S, Monteiro E, Mechbal N, Rebillat M, Chinesta F (2023) Hybrid twin of rtm process at the scarce data limit. Int J Mater Form 16(4):40

26. Tezduyar TE (2001) Finite element methods for flow problems with moving boundaries and interfaces. Arch Comput Methods Eng 8:83–130

27. Loon R, Anderson PD, Vosse FN, Sherwin SJ (2007) Comparison of various fluid-structure interaction methods for deformable bodies. Comput & Struct 85(11–14):833–843

28. Chinesta F, Mabrouki T, Ramon A (2002) Some difficulties in the flow front treatment in fixed mesh simulations of composites forming processes

29. Chinesta F, Cueto E, Quintela P, Paredes J (2004) Induced anisotropy in foams forming processes: modelling and simulation. J Mater Process Technol - J MATER PROCESS TECHNOL 155:1482–1488. https://doi.org/10.1016/j.jmatprotec.2004.04.378

30. Garcia J, Gascón L, Chinesta F (2003) A fixed mesh numerical method for modelling the flow in liquid composites moulding processes using a volume of fluid technique. Comput Methods Appl Mech Eng 192(7–8):877–893

31. Piegl L, Tiller W (1996) The NURBS Book. Springer, Berlin, Germany

32. Lange K, Lange K (2010) Singular value decomposition. Numerical analysis for statisticians, pp 129–142
33. Brunton SL, Kutz JN (2022) Data-driven Science and Engineering: Machine Learning, Dynamical Systems, and Control. Cambridge University Press, Cambridge, UK

**Publisher's Note** Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.