# Minimizing Deformations during HP MJF 3D Printing

Karel Ráž, Zdeněk Chval, Sacha Thomann

► To cite this version:

Karel Ráž, Zdeněk Chval, Sacha Thomann. Minimizing Deformations during HP MJF 3D Printing. Materials, 2023, 16 (23), pp.7389. 10.3390/ma16237389. hal-05245312

HAL Id: hal-05245312

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

Submitted on 8 Sep 2025

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
Article

# Minimizing Deformations during HP MJF 3D Printing

Karel Ráž <fnref n="1" />,* , Zdeněk Chval <fnref n="1" /> and Sacha Thomann <fnref n="2" /> ID

**Abstract:** (1) Background: The purpose of this study was to investigate deformations that occur during additive manufacturing by the HP (Hewlett-Packard) Multi Jet Fusion (MJF) process. These deformations affect the final properties of 3D-printed parts, and proper compensating technology has to be developed in order to minimize these deformations. (2) Methods: Parts were printed with powder composed of nylon plastic infused with glass beads (PA12GB). The HP MJF technology was used during investigations. All parts (specimens) were measured at different points over an extended period to follow the deformations at each point. Different finite element simulations were performed to compare them with real results and assess the viability of using simulations to save time. Various modules of the Digimat software, such as additive manufacturing (AM), material focused (MF), finite element (FE), and computer-aided engineering (CAE), were used to run the simulations. (3) Results: It was found that the printing position of the part in the printer had an impact on deformations. When the part was simulated in a tilted position but alone (deformation: 7.19 mm), the value of the deformation was 1.49 mm greater than when the other parts (two comparable parts) were simulated at the same time (deformation: 5.7 mm). The difference between the simulation with the three parts together (deformation: 5.7 mm) and reality (deformation: 3.44 mm) was 2.26 mm. Finally, the difference between the simulated single part (deformation: 7.19 mm) and the real part (deformation: 3.44) was 3.75 mm. (4) Conclusions: The results of this study will contribute to a better understanding of deformation mechanisms and will suggest solutions for improving the quality of printed parts. Three-dimensional printing is a rapidly growing technology that offers numerous possibilities across various fields. However, one commonly encountered issue is the deformation of printed parts. Methods for minimizing deformations were studied during the 3D printing process using HP MJF technology. Various factors contributing to deformation were investigated, and different techniques for reducing them were explored.

**Keywords:** 3D printing; additive manufacturing; PA12GB; Digimat; deformation; MJF

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

check for
updates

**Citation:** Ráž, K.; Chval, Z.; Thomann, S. Minimizing Deformations during HP MJF 3D Printing. *Materials* **2023**, *16*, 7389.
https://doi.org/10.3390/ma16237389

Academic Editors: Young-Hag Koh and Jiapeng Sun

Received: 4 October 2023
Revised: 30 October 2023
Accepted: 25 November 2023
Published: 28 November 2023

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

**Copyright:** © 2023 by the authors.
Licensee MDPI, Basel, Switzerland.
This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/).

## 1. Introduction

Three-dimensional printing of industrial parts is more and more used to create complex and cheaper products. Additive manufacturing has advantages like design freedom and cost reduction, but there are also disadvantages like the printing time and the appearance of warping deformation in the part due to high temperature [1,2]. There exist different types of additive manufacturing technologies, such as SLS, FDM, FFF, or MJF. For each of these printing techniques, there are different deformations after the printing process [3,4]. This research is focused on deformation and warpage prediction using FEM simulation tools. Deformation after the production process can occur due to various effects (high temperature, wrong positioning in printing area, etc.).

Deformation in 3D-printed parts is a common issue and can occur for various reasons.

Understanding this mechanism of deformation is crucial for improving the quality and reliability of 3D-printed components. Common problems associated with deformation

<footnote n="1">¹ Faculty of Mechanical Engineering, Regional Technological Institute, University of West Bohemia, Univerzitni 8, 306 14 Plzen, Czech Republic; zdchval@fst.zcu.cz</footnote>

<footnote n="2">² Polytech Montpellier, MI4, Université de Montpellier, Pl. Eugène Bataillon, 34090 Montpellier, France; sacha.thomann@etu.umontpellier.fr</footnote>

<footnote>* Correspondence: kraz@fst.zcu.cz; Tel.: +420-377-638-751</footnote>
in 3D-printed parts are warping (caused mainly by different cooling rates across the part), cracking (caused by internal stresses), shrinkage (caused by thermal contractions), and layer misalignments. This research is focused on shrinkage and warping [5].

Additive manufacturing (AM) has revolutionized the way we design, produce, and interact with objects, opening up new horizons in various industries. One of the most promising techniques in the field of AM is the high-performance Multi Jet Fusion (HP MJF) process [6,7], which has gained significant attention for its ability to fabricate complex, functional 3D parts with exceptional speed and precision. However, like any manufacturing process, HP MJF is not without its challenges and limits. One critical aspect that demands a comprehensive understanding is the potential deformation of 3D parts during the printing process.

In this article, the deformation in 3D parts produced through the HP MJF process (see Figure 1) is deeply studied and described. The underlying mechanisms that lead to distortion and warping, such as residual stresses, are explored. Understanding these factors is crucial for optimizing printing parameters and ensuring the manufacturability of complex geometries [8–10].

<figure><img src="image_4.png" /><figcaption>**Figure 1.** A diagram of the MJF process.</figcaption></figure>

The influence of the production method on part deformation and strength is investi-
gated so that engineers can leverage this knowledge to design parts that are more resilient
to deformation while maintaining their intended functionalities. The interplay between
part orientation, support structures, and build layout plays a pivotal role in mitigating
deformation and achieving successful outcomes.

In conclusion, this article aims to shed light on the phenomenon of deformation
in 3D parts fabricated using the HP MJF process. By comprehending the underlying
mechanisms, designers, engineers, and manufacturers can make informed decisions to
minimize deformation and optimize the overall manufacturing process [11–13]. Through
this exploration, this research aspires to unlock the full potential of the HP MJF technology
in diverse industries worldwide, where complex, high-performance 3D polymer parts are
at the forefront of innovation.

**2. Materials and Methods**

The HP Multi Jet Fusion (MJF) process involves powdered materials being selectively
fused layer by layer to build complex 3D structures (Figure 1) [14]. This advanced 3D
printing technology uses a powder-based additive manufacturing process [15,16]. A layer
of powdered unfused material is spread over the built platform, and inkjet arrays selectively
apply a fusing agent and detailing agent onto the powder according to the digital design.
A heating element then fuses the treated areas together to form a solid layer. This layer-by-
layer process continues until the 3D print is complete [17]. The final properties of the part
are highly dependent on the layer quality and thickness [18–20].

The material considered in this research was PA12GB. The thickness of one layer
of PA12GB is around 0.08 mm, and this technology results in actual productivity of
4115 cm³/h [21,22]. The processing temperature for PA12GB is around 160 to 180 °C [23].
This temperature has to be understood with respect to the powder melting point of PA12GB,
which is 186 °C. The material PA12GB has a tensile strength of 30 MPa (according to the of-
ficial datasheet). The guaranteed value of elongation at break, as provided by the producer of the material, is 6.5% [24,25].

The HP MJF technique offers impressive benefits, such as high speed, precision, and the ability to produce complex parts with excellent detail resolution and mechanical properties. It also introduces complexities related to thermal, mechanical, and material properties. These factors can contribute to varying degrees of deformation, impacting the final accuracy, dimensional stability, and overall quality of the printed parts.

HP MJF has applications in various industries, including aerospace, automotive, healthcare, and consumer goods [26]. Its versatility and efficiency make it a preferred choice for rapid prototyping, customized manufacturing, and production of end-use parts with high reliability. However, careful consideration and optimization of the printing parameters are necessary to mitigate potential deformations and ensure optimal results in the final printed parts [27].

Furthermore, the impact of the material properties on deformation behavior was examined, recognizing that different powder materials react differently under the heat and pressure of the MJF process. Identifying the ideal material and manufacturing settings and considering correct mechanical characteristics is vital for obtaining highly accurate and dimensionally stable parts [28,29].

### 3. Microstructure

To carry out this study, the digital definition of the characteristics of the used material had to be the first step. This approach would enable us to obtain a virtual representation that is as faithful as possible to reality, allowing us to anticipate the behavior expected under real conditions.

#### 3.1. Material Behavior

To begin with, FEM analysis was used to numerically define the material, i.e., PA12GB. This is a composite material made by incorporating glass beads into polyamide 12 (PA12), which is a type of nylon, to reinforce the material. Material properties for the combination of PA12 and glass beads were obtained using Digimat MF (see Table 1) and the material datasheet (see Table 2).

<table><caption>**Table 1.** Engineering constants available from the Digimat MF library.</caption><thead><tr><th>Parameters</th><th>PA12</th><th>GB</th><th>PA12GB</th></tr></thead><tbody><tr><td>Density (g/cm³)</td><td>1.01</td><td>2.54</td><td>1.33</td></tr><tr><td>E (GPa)</td><td>1.735</td><td>72.000</td><td>2.685</td></tr><tr><td>G (GPa)</td><td>N/A</td><td>N/A</td><td>1.03</td></tr><tr><td>N</td><td>0.39</td><td>0.22</td><td>0.37</td></tr><tr><td>Yield stress</td><td>22</td><td>N/A</td><td>N/A</td></tr><tr><td>Hardening modulus</td><td>18</td><td>N/A</td><td>N/A</td></tr><tr><td>Hardening exponent</td><td>70</td><td>N/A</td><td>N/A</td></tr><tr><td>Linear hardening modulus</td><td>15</td><td>N/A</td><td>N/A</td></tr></tbody></table>

<table><caption>**Table 2.** Engineering constants from the manufacturer's datasheet.</caption><thead><tr><th>Parameters</th><th>PA12GB</th></tr></thead><tbody><tr><td>Density (g/cm³)</td><td>1.30</td></tr><tr><td>E (GPa)</td><td>2.850</td></tr><tr><td>G (GPa)</td><td>1.03</td></tr><tr><td>N</td><td>0.21</td></tr></tbody></table>

The constitutive law used to model PA12 utilized elastoplastic law for the J2 plasticity model and exponential and linear laws for the isotropic hardening model [30]. An isotropic elastic model was used to model the glass beads [31].
Now that the material properties had been defined, the microstructure properties had to be described. The PA12GB microstructure is divided into two phases:

* Phase one: PA12 phase, which is a matrix phase;

* Phase two: Glass beads phase, which is an inclusion phase with many spherical inclusions, the total mass fraction of which is 0.4 of phase one.

The elastic domain simulation with Hooke’s law is possible according to the following equation:

$$E = \sigma \times \epsilon \qquad (1)$$

A failure model had to be defined within Digimat in order to generate plastic domain behavior and run the analysis. The 3D transversely isotropic Tsai–Hill model was chosen because it incorporates interactions between stress components, even though it does not take into account the tension–compression asymmetry [32–34], which is negligible in our case. The 3D transversely isotropic Tsai–Hill model under deformation is written in Digimat as follows [35]:

$$f_A = \sqrt{F_A(\epsilon)}, \qquad (2)$$

With

$$F_A(\epsilon) = \frac{\epsilon_{11}^2}{X^2} - \frac{\epsilon_{11}(\epsilon_{22} + \epsilon_{33})}{X^2} + \frac{(\epsilon_{22}^2 + \epsilon_{33}^2)}{Y^2} + \left(\frac{1}{X^2} - \frac{2}{Y^2}\right)\epsilon_{22}\epsilon_{33} + \frac{((2\epsilon_{12})^2 + (2\epsilon_{13})^2)}{S^2} + \left(\frac{1}{Y^2} - \frac{1}{4X^2}\right)(2\epsilon_{23})^2 \qquad (3)$$

where

* $\epsilon_{11}$: Principal strain in the fiber direction.

* $\epsilon_{22}$: Principal strain perpendicular to the fiber.

* $\epsilon_{33}$: Principal strain perpendicular to the fiber.

* $\epsilon_{23}$: Strain perpendicular to the fiber.

* $f_A$: Failure indicator; failure occurs when indicator is more or equal to 1.

* X, Y, S are unitless maximum strain parameters.

* In our case of rounded glass beads, this equation can be simplified because the fiber is globular and its properties are the same in all three directions. This complex equation is used because of its possible use for other (nonglobular) phases [36,37].

When everything is defined, Digimat can be used to obtain the stress versus strain curve, but this is only a simulated result. Tensile tests on test specimens with various orientations during printing were performed in order to determine the material’s experimental curves. This made it possible to compare the experimental and simulated results.

The highest difference was for orientation 0° and 90° [38,39], while orientation 45° had an average value between the two. The example of specimen orientation is shown in Figure 2.

<figure><img src="image_5.png" /><figcaption>Figure 2. Different orientations for 3D printing.</figcaption></figure>
Each result is plotted onto the same graph in Figure 3. The gray dashed curve represents the stress-strain curve of PA12GB at an angle of 0 degrees, while the red curve represents the curve for PA12GB at an angle of 90 degrees. The difference was observed in the plastic domain, where failure occurred much earlier for the 90° test compared to the 0° test.

<figure><img src="image_6.png" /><figcaption>**Figure 3.** Stress-strain curves for PA12GB.</figcaption></figure>

The Young's modulus was higher for the 90° orientation, and the rupture occurred at lower strain (compared to 0°) because loading is affected more by the cohesivity between layers [40]. The linear behavior (Young's modulus from the datasheet) of the material is shown in Figure 3 using the green dotted line. The last curve in yellow represents the curve generated by Digimat, which closely matches the tensile curve at 0°. The material created via Digimat therefore seems consistent with reality and can be used to predict deformation. The obtained values have to be comparable to the official value of tensile strength from the HP datasheet, which is 30 MPa. The difference is caused by the not absolutely homogenous printing process in reality compared to the ideal printing process achieved when obtaining the value for the datasheet. The influence of other parts within the printing area on the strength (heat transfer and temperature distribution) is significant (tensile strength decreasing by 10–15%) [41,42].

The orientation of parts during the printing process also affects the accuracy of the print. This geometric accuracy is negligible with respect to the thermal deformation caused by the design and thermal expansion of the material.

## 3.2. *Microstructure Simulation*

An examination of the material's microstructure, represented by the representative volume element (as depicted in Figure 4), was conducted in order to simulate part deflection and behavior during the additive manufacturing process. This analysis aimed to visualize stress distribution within the material, thereby enhancing our internal comprehension of its properties. The stiffness and thermal properties of PA12GB results from this material modeling. These crucial values were internally used within Digimat for further simulations of the additive manufacturing process. This step was necessary for correct material modeling within FEM analysis.

Digimat FE was used to generate the microstructure composed of the PA12 matrix and numerous spherical inclusions of glass beads. When a mechanical strain loading of
0.15 mm/mm was applied, von Mises stresses were obtained, as shown in Figure 5. The stress/contraction in the main direction are shown in Figure 6.

<figure><img src="image_7.png" /><figcaption>Figure 4. Representative volume element (RVE).</figcaption></figure>

<figure><img src="image_8.png" /><figcaption>Figure 5. von Mises stress, $\sigma_{Vm}$, in microstructure (MPa).</figcaption></figure>

<figure><img src="image_9.png" /><figcaption>Figure 6. Stress in direction 11 (main direction) $\sigma_{11}$ (deformation amplified 14 times) (MPa).</figcaption></figure>
As can be seen in Figures 5 and 6, the stresses were mainly distributed in the glass beads, while they were lower in the PA12 matrix.

## 4. Deflection Simulation

Simulations of 3D-printed polymer parts were carried out in order to observe the different deformations that can occur in these parts as a function of different parameters. The software Digimat 2023.3 AM was used for this task.

The first step of simulation is definition of the material characterization, followed by importation to a CAD model to properly define orientation and positioning. The simulation of the building process depends on the quality of the mesh (number of voxels). This step results in coupled thermal-structural FEM simulation, which can be used to obtain results of deflection and prediction of defects. Elements of each layer are turned on with respect to the movement of the printing head.

### 4.1. Part Design

A rectangular plate with a thick edge on half the perimeter and a rib running diagonally across the piece was chosen. The part was designed in Ansys and measured 150 mm × 200 mm × 20 mm (Figure 7). It is obvious that the design of this part will lead to deformation after the printing process. This deformation is caused by the CAD design (one-sided rib). The thermal behavior of the polymer PA12GB during cooling causes shrinkage. This shrinkage is significant in areas where more materials are concentrated. Therefore, this will be the part that deforms upwards in the direction of the rib location. The correct design should be with ribs on both sides, but this design was not possible to use because of the function of the part.

<figure><img src="image_10.png" /><figcaption>Figure 7. Design of the test part.</figcaption></figure>

### 4.2. Digimat AM Setup

To effectively simulate the deformation of the component, it is imperative to select Digimat AM parameters that correspond with real-world conditions. The printing chamber's dimensions aligned with those employed in the HP printer (380 × 284 × 380 mm). The simulation encompassed a comprehensive full-build simulation, accounting not only for the part itself but also for the surrounding powder material [43].

As outlined in the Materials and Methods selection, we opted for PA12GB based on its technical characteristics outlined in the previous section, which closely approximate the actual material. We maintained the printing process parameters without alteration with a consistent printing temperature of 330 degrees Celsius.
However, one parameter was modified. The thickness of each layer of powder, normally 0.08 mm, was increased to 1 mm for reasons relating to computer processing power and storage. For this reason, at the Digimat meshing stage of the part, the voxel size was set to 2 mm, generating a mesh with 410,688 voxels in full build. The process took 963 min, or 16 h 3 min, on a computer with an Intel(R) Core(TM) i5-4590 CPU @3.30 GHz. Given that the number of voxels increases quadratically with decreasing voxel size, this was the best compromise between time and accuracy, as shown in Figures 8 and 9.

<figure><img src="image_11.png" /><figcaption>**Figure 8.** Quantity of voxels depending on mesh size.</figcaption></figure>

<figure><img src="image_12.png" /><figcaption>**Figure 9.** Simulation time depending on voxel number.</figcaption></figure>

Meshing was then carried out on Digimat layer by layer. The mesh quality is shown in Figure 10.

### 4.3. *Results from Digimat FEM Analysis*

Once the part had been meshed, three simulations were carried out to determine the differences generated by the orientation of the print. The first one was a simulation of the single part oriented horizontally (flat position) in the printing chamber. The simulation considered the coupled structural-thermal solution with element birth function. Heat transfer coefficients were set up by the solution process with respect to the selected material [44,45].
<figure><img src="image_13.png" /><figcaption>**Figure 10.** Part meshed with Digimat—detail of a corner.</figcaption></figure>

The deformations for the flat orientation occurred in a circular pattern originating from the center of the component, as illustrated in Figure 11. This occurrence was due to the cooling process, where the part initially cooled from its outer edges towards the center, causing it to contract inward to its midpoint. It is important to note that the most significant deformations occurred at the corners, where there was a lower amount of material. Nevertheless, it is worth noting that the overall deformation of the part appeared to be minimal, with a maximum deflection of only 1.03 mm.

<figure><img src="image_14.png" /><figcaption>**Figure 11.** Simulated part in flat position in the printer with PA12GB (mm).</figcaption></figure>

The second simulation was of the piece tilted as it would be when printed in real life (approximately 45° orientation in the printing chamber). Here, the corner with the most material was at the lowest point (the corner where all three ribs were connected), meaning that it would be printed first. Figure 12 shows that the maximum deformation was located at the end of the rib, i.e., at the last printed area.

The deformation was caused by one main factor: the fact that the part had already begun to cool on the lower part, thus gradually creating greater deformation at the upper end of the part. The maximum deflection in this simulation was 7.72 mm. The difference between this and the previous result was caused by completely different heat transfer during the printing process.

The last simulation was of the part at the same angle as before but this time as in the real world, with the part being printed at the same time as the two others (slightly different parts). Here, the simulation took place with three parts, as shown in Figure 13, but this research was focused only on the deformations developing on the initial part (the same part as in the previous simulations).
<figure><img src="image_15.png" /><figcaption>**Figure 12.** Simulated part in tilted position in the printer (mm).</figcaption></figure>

<figure><img src="image_16.png" /><figcaption>**Figure 13.** Positioning of parts for the simulation as in reality.</figcaption></figure>

The part simulated in the previous section was printed as shown with the red arrow in Figure 14, which represents a different orientation than flat. It is obvious that the testing parts were produced within the print job with other parts. The thermal effect of these parts was neglected. These smaller parts were not part of this research; they were added to the print job by the printer operator, and the correct position of all parts was set up according to the official HP software Netfabb, which was used to communicate with the HP 4200 printer.

The results were almost the same as those from the previous simulation; however, the deformations were less significant: 5.7 mm for this one (Figure 15) compared to 7.72 mm for the previous one. This is due to the fact that when a part is finished being printed, it remains warmer, and the printing of the other parts continues to heat the powder generally, thus increasing the cooling time of the part while homogenizing the temperature.

It can be concluded that the part deforms differently depending on the position in which it is printed. If minimal deformation is desired, it should be printed flat with other parts around it, even though it may tend to deform at multiple points. On the other hand, if a larger deformation at a single point is negligible, it might be preferable to print it inclined. Printing flat parts horizontally is generally minimally used because it causes problems with surface quality, such as sinks or 'elephant skin'. This flat orientation is shown just for comparing. It will not be used in reality because it will cause the detailing and fusing agent to remain on the upper surface during printing.
<figure><img src="image_17.png" /><figcaption>Figure 14. Parts in the HP MJF 4200 print box.</figcaption></figure>

<figure><img src="image_18.png" /><figcaption>Figure 15. Simulated part with other parts (which are hidden) as in reality (mm).</figcaption></figure>

**5. Real Deflection**

The measurement of deformations after the printing process has to be performed in order to validate the simulations from the virtual prediction. The difference between the simulation and reality has to be quantified.

In order to check whether the simulated results are correct, it is necessary to observe the deformations that occur in the real part. Figure 16 shows the deformed printed part. The dimension points are shown in order to measure deflection in time.

<figure><img src="image_19.png" /><figcaption>Figure 16. Printed part with measured parameters.</figcaption></figure>
## 5.1. Measurements

Measurements were taken at different times after the printing process. A total of 15 measurements (15 specimens) were performed at each time in order to obtain statistically correct results (by the 3D measuring station). The first measurement was just after the printing had finished, then two and five days later, and then every week for three weeks. To obtain long-term information, measurements were also taken six and eight weeks after printing.

From Table 3, it is clear that almost all deformations were negligible in view of the values achieved.

<table><caption>**Table 3.** Deformations over time for the initial part in millimeters.</caption><thead><tr><th>Value in Time (mm)</th><th>δh</th><th>h</th><th>H</th><th>l</th><th>e</th><th>b</th></tr></thead><tbody><tr><td>Designed in CAD</td><td>0.00</td><td>10.00</td><td>20.00</td><td>150.00</td><td>2.00</td><td>16.45</td></tr><tr><td>Immediately after print</td><td>3.10</td><td>10.15</td><td>20.20</td><td>151.00</td><td>2.00</td><td>16.37</td></tr><tr><td>2 days after print</td><td>3.82</td><td>10.13</td><td>20.11</td><td>149.95</td><td>2.00</td><td>16.41</td></tr><tr><td>5 days after print</td><td>3.35</td><td>10.03</td><td>20.18</td><td>149.66</td><td>1.93</td><td>16.38</td></tr><tr><td>1 week after print</td><td>3.52</td><td>10.3</td><td>20.20</td><td>149.79</td><td>1.89</td><td>16.35</td></tr><tr><td>2 weeks after print</td><td>3.68</td><td>10.05</td><td>20.21</td><td>149.81</td><td>1.86</td><td>16.31</td></tr><tr><td>3 weeks after print</td><td>4.20</td><td>9.99</td><td>20.24</td><td>149.81</td><td>1.90</td><td>16.37</td></tr><tr><td>6 weeks after print</td><td>3.47</td><td>9.97</td><td>20.19</td><td>149.80</td><td>1.89</td><td>16.36</td></tr><tr><td>8 weeks after print</td><td>3.44</td><td>10.01</td><td>20.22</td><td>149.80</td><td>1.89</td><td>16.34</td></tr></tbody></table>

The only non-negligible value was that at the end of the rib (δh), which should be zero. This value increased then decreased and seemed to stabilize at around 3.42 mm of deformation (Figure 17). It has to be mentioned that the part was designed specially without an opposite rib (nonbalanced design) in order to achieve some deformation to compare.

<figure><img src="image_20.png" /><figcaption>**Figure 17.** Measurement of the deflection at the end of the rib over time (δh).</figcaption></figure>

## 5.2. Discussion and Comparison of Results

By comparing the results of advanced simulations with the data obtained from real experiments, a correlation between the virtual world and the physical world can be established. It is obvious that the simulations and reality appeared to be consistent as the deformation occurred at the same point. However, there were a few differences between the values of the maximum deflection in the corner (δh).
When the part was simulated in a tilted position but alone (deformation: 7.19 mm), the value of the deformation was 1.49 mm greater than when the other parts were simulated at the same time (deformation: 5.7 mm). The difference between the simulation with the three parts together (deformation: 5.7 mm) and reality (deformation: 3.44 mm) was 2.26 mm. Finally, the difference between the simulated single part (deformation: 7.19 mm) and the real part (deformation: 3.44) was 3.75 mm. This is not an insignificant difference. Further research and development in this area must be carried out to bring virtual and real printing closer together and thereby design parts with minimal deformation. All results are shown in Table 4.

<table><caption>**Table 4.** Comparison of results.</caption><thead><tr><th></th><th>Desired Value in CAD</th><th>FEM—Flat Position</th><th>FEM—Tilted Position One Part</th><th>FEM—Tilted Position Three Parts</th><th>Steady State Value—Real Measurement</th></tr></thead><tbody><tr><th>δh</th><td>0.00</td><td>1.03</td><td>7.19</td><td>5.70</td><td>3.44</td></tr></tbody></table>

## 6. Conclusions

This study delved into the investigation of deformations that occur during the Hewlett-Packard (HP) Multi Jet Fusion (MJF) 3D printing process. The study focused on minimizing deformations by exploring different printing orientations and highlighting the factors that influence the deformation of printed parts.

Simulation results were compared with experimental data, demonstrating good consistency between the two, by considering the printing position, the material, the size, and even the quantity of pieces during printing.

The difference between reality (after stabilization of dimension in time) and simulation was 2.26 mm. This is the lowest achieved error with maximal accuracy of simulation for a full printing job. Further approximation of results is possible only by FEM mesh refinement. This is limited by the actual available computing power.

The insights gained into the deformation mechanisms and simulation methods provide valuable information for optimizing the design of parts printed using the HP MJF technology. This research contributes to enhancing the quality of printed parts and harnessing the full potential of the HP MJF technology across various industrial applications. In summary, this study underscores the significance of comprehending deformation mechanisms and employing effective strategies in order to minimize deformations during HP MJF 3D printing. Future research will continue with the description of material models of different composite polymers and refining results in order to minimize differences between simulations and the parts produced by the MJF (and other) technology. The main limit of this study is the limited computing power. The actual considered size of one element was 2 mm. This value should be lower than the layer size during the actual printing process (0.08 mm) in order to obtain more precise results. This means that elements that are 25 times smaller should be used for FEM simulations. This size of elements was not possible to simulate with the computational power that was used.

**Author Contributions:** Conceptualization, K.R.; methodology, K.R. and S.T.; numerical simulations, K.R. and S.T.; experimental data processing, S.T. and Z.C.; writing—original draft preparation, K.R. and Z.C. All authors have read and agreed to the published version of the manuscript.

**Funding:** This work was supported by the Grant Agency of the University of West Bohemia, grant "Long-term conceptual development of a research organization".

**Institutional Review Board Statement:** Not applicable.

**Informed Consent Statement:** Not applicable.

**Data Availability Statement:** Data are contained within the article.

**Acknowledgments:** The article was prepared under project "Long-term conceptual development of a research organization".
**Conflicts of Interest:** The authors declare no conflict of interest.

**References**

1. Chua, C.K.; Wong, C.H.; Yeong, W.Y. Benchmarking for Additive Manufacturing. In *Standards, Quality Control and Measurement Sciences in 3D Printing and Additive Manufacturing*, 1st ed.; Academic Press: Cambridge, MA, USA, 2017; pp. 181–212.

2. Shahrubudin, N.; Lee, T.C.; Ramlan, R. An Overview on 3D Printing Technology: Technological, Materials, and Applications. *Procedia Manuf.* **2019**, *35*, 1286–1296. [CrossRef]

3. Farkas, A.Z.; Galatanu, S.-V.; Nagib, R. The Influence of Printing Layer Thickness and Orientation on the Mechanical Properties of DLP 3D-Printed Dental Resin. *Polymers* **2023**, *15*, 1113. [CrossRef] [PubMed]

4. Liu, X.; Tey, W.S.; Choo, J.Y.C.; Chen, J.; Tan, P.; Cai, C.; Ong, A.; Zhao, L.; Zhou, K. Enhancing the mechanical strength of multi jet fusion-printed polyamide 12 and its glass fiber-reinforced composite via high-temperature annealing. *Addit. Manuf.* **2021**, *46*, 102205. [CrossRef]

5. Archenti, A.; Maffei, A. Proceedings of the International Conference on Advanced Manufacturing Engineering and Technologies NEWTECH. 2013. Available online: https://www.diva-portal.org/smash/get/diva2:660817/FULLTEXT09.pdf#page=163 (accessed on 15 September 2023).

6. Geng, Z.; Bidanda, B. Geometric precision analysis for Additive Manufacturing processes: A comparative study. *Precis. Eng.* **2021**, *69*, 68–76. [CrossRef]

7. Kechagias, J.; Chaidas, D.; Vidakis, N.; Salonitis, K.; Vaxevanidis, N.M. Key Parameters Controlling Surface Quality and Dimensional Accuracy: A Critical Review of FFF Process. *Mater. Manuf. Process.* **2022**, *37*, 963–984. [CrossRef]

8. Yelamanchi, B.; Mummareddy, B.; Santiago, C.C.; Ojoawo, B.; Metsger, K.; Helfferich, B.; Zapka, J.; Silani, F.; MacDonald, E.; Cortes, P. Mechanical and fatigue performance of pressurized vessels fabricated with multi jet fusion™ for automotive applications. *Addit. Manuf.* **2021**, *44*, 102048. [CrossRef]

9. Alghamdi, S.; John, S.; Choudhury, N.; Dutta, N. Additive Manufacturing of Polymer Materials: Progress, Promise and Challenges. *Polymers* **2021**, *13*, 753. [CrossRef] [PubMed]

10. Richards, R.W.; Sims, D. Glass-Filled Thermoplastics—Effects of Glass Variables and Processing on Properties. *Composites* **1971**, *2*, 214–220. [CrossRef]

11. Li, Y.; Long, S.; Liu, Q.; Lv, H.; Liu, M. Resistive Switching Performance Improvement via Modulating Nanoscale Conductive Filament, Involving the Application of Two-Dimensional Layered Materials. *Small* **2017**, *13*, 1604306. [CrossRef]

12. Cai, C.; Tey, W.S.; Chen, J.; Zhu, W.; Liu, X.; Liu, T.; Zhao, L.; Zhou, K. Comparative study on 3d printing of polyamide 12 by selective laser sintering and multi jet fusion. *J. Mater. Process. Technol.* **2021**, *228*, 116882. [CrossRef]

13. Atakok, G.; Kam, M.; Koc, H.B. Tensile, three-point bending and impact strength of 3d printed parts using pla and recycled pla filaments: A statistical investigation. *J. Mater. Res. Technol.* **2022**, *18*, 1542–1554. [CrossRef]

14. Lee, P.-H.; Chung, H.; Lee, S.W.; Yoo, J.; Ko, J. Review: Dimensional Accuracy in Additive Manufacturing Processes. In Proceedings of the ASME 2014 International Manufacturing Science and Engineering Conference, Detroit, MI, USA, 9–13 June 2014.

15. Alomarah, A.; Ruan, D.; Masood, S.; Gao, Z. Compressive properties of a novel additively manufactured 3d auxetic structure. *Smart Mater. Struct.* **2019**, *28*, 085019. [CrossRef]

16. 3DPrint.com Company. HP Reveals More Info about Their Multi Jet Fusion 3D Printing Technology, Plans for Second 3D Printer. 2016. Available online: https://3dprint.com/113630/hp-multi-jet-fusion-plans-info/ (accessed on 15 September 2023).

17. Shim, J.S.; Kim, J.-E.; Jeong, S.H.; Choi, Y.J.; Ryu, J.J. Printing accuracy, mechanical properties, surface characteristics, and microbial adhesion of 3D-printed resins with various printing orientations. *J. Prosthet. Dent.* **2020**, *124*, 468–475. [CrossRef] [PubMed]

18. Nowacki, J.; Sieczkiewicz, N. Problems of determination of MultiJet 3D printing distortions using a 3D scanner. *Arch. Mater. Sci. Eng.* **2020**, *103*, 30–41. [CrossRef]

19. Fradl, D.; Panditaratne, J.; Bi, J.; Fu, R.; Oancea, V. Finite element simulation of the multi jet fusion (mjf™) process using abaqus. In Proceedings of the Scince of the Age of Experience Conference, Boston, MA, USA, 23–27 May 2017.

20. HP Development Company, L.P. HP 3D High Reusability PA 12 Glass Beads. 2019. Available online: https://static1.sw-cdn.net/files/cms/materials/data-sheets/HP-MJF-PA12GB-datasheet.pdf (accessed on 15 September 2023).

21. Mazzanti, V.; Malagutti, L.; Mollica, F. FDM 3D Printing of Polymers Containing Natural Fillers: A Review of Their Mechanical Properties. *Polymers* **2019**, *11*, 1094. [CrossRef] [PubMed]

22. Chan, H.K.; Griffin, J.; Lim, J.J.; Zeng, F.; Chiu, A.S.F. The Impact of 3D Printing Technology on the Supply Chain: Manufacturing and Legal Perspectives. *Int. J. Prod. Econ.* **2018**, *205*, 156–162. [CrossRef]

23. Das, A.K.; Agar, D.A.; Rudolfsson, M.; Larsson, S.H. A Review on Wood Powders in 3D Printing: Processes, Properties and Potential Applications. *J. Mater. Res. Technol.* **2021**, *15*, 241–255. [CrossRef]

24. Avanzini, A.; Battini, D.; Pandini, S. Static and fatigue behavior in presence of notches for polyamide 12 (pa12) additively manufactured via multi jet fusion™ process. *Int. J. Fatigue* **2022**, *161*, 106912. [CrossRef]

25. Osswald, P.V.; Obst, P.; Mazzei Capote, G.A.; Friedrich, M.; Rietzel, D.; Witt, G. Failure criterion for pa 12 multi-jet fusion additive manufactured parts. *Addit. Manuf.* **2021**, *37*, 101668. [CrossRef]

26. Raz, K.; Chval, Z.; Milsimerova, A. Thermal specification of 3d printed injection moulds made from pa12gb. *IOP Conf. Ser. Mater. Sci. Eng.* **2021**, *1199*, 012009. [CrossRef]
27. Rosso, S.; Meneghello, R.; Biasetto, L.; Grigolato, L.; Concheri, G.; Savio, G. In-depth comparison of polyamide 12 parts manufactured by multi jet fusion and selective laser sintering. *Addit. Manuf.* **2020**, *36*, 101713. [CrossRef]

28. O'Connor, H.J.; Dowling, D.P. Comparison between the properties of polyamide 12 and glass bead filled polyamide 12 using the multi jet fusion printing process. *Addit. Manuf.* **2020**, *31*, 100961. [CrossRef]

29. Yang, F.; Zobeiry, N.; Mamidala, R.; Chen, X. A review of aging, degradation, and reusability of pa12 powders in selective laser sintering additive manufacturing. *Mater. Today Commun.* **2023**, *34*, 105279. [CrossRef]

30. Suder, J.; Bobovsky, Z.; Mlotek, J.; Vocetka, M.; Zeman, Z.; Safar, M. Experimental analysis of temperature resistance of 3d printed pla components. *MM Sci. J.* **2021**, *1*, 4322–4327. [CrossRef]

31. Mehrpouya, M.; Tuma, D.; Vaneker, T.; Afrasiabi, M.; Bambach, M.; Gibson, I. Multimaterial powder bed fusion techniques. *Rapid Prototyp. J.* **2022**, *28*, 1–19. [CrossRef]

32. Belter, J.T.; Dollar, A.M. Strengthening of 3d printed fused deposition manufactured parts using the fill compositing technique. *PLoS ONE* **2015**, *10*, e0122915. [CrossRef] [PubMed]

33. Heitkamp, T.; Girnth, S.; Kuschmitz, S.; Waldt, N.; Klawitter, G.; Vietor, T. Experimental and Numerical Investigation of the Mechanical Properties of 3D-Printed Hybrid and Non-Hybrid Composites. *Polymers* **2023**, *15*, 1164. [CrossRef]

34. Gadelmoula, A.; Aldahash, S.A. Tribological Properties of Glass Bead-Filled Polyamide 12 Composite Manufactured by Selective Laser Sintering. *Polymers* **2023**, *15*, 1268. [CrossRef]

35. Tiwari, A.S.; Yang, S. Energy Consumption Modeling of 3D-Printed Carbon-Fiber-Reinforced Polymer Parts. *Polymers* **2023**, *15*, 1290. [CrossRef]

36. Issabayeva, Z.; Shishkovsky, I. Prediction of The Mechanical Behavior of Polylactic Acid Parts with Shape Memory Effect Fabricated by FDM. *Polymers* **2023**, *15*, 1162. [CrossRef]

37. Duvoisin, C.; Horst, D. Additive Manufacturing at Industry 4.0: A Review. *Int. J. Eng. Tech. Res.* **2018**, *8*, 3–8.

38. Bandyopadhyay, A.; Gualtieri, T.; Heer, B.; Bose, S. Introduction to Additive Manufacturing. In *Additive Manufacturing*, 2nd ed.; Bandyopadhyay, A., Bose, S., Eds.; CRC Press: Boca Raton, FL, USA, 2019; pp. 1–23.

39. Tawlik, M.; Nemat-Alla, M.; Dewidar, M. Enhancing the properties of aluminum alloys fabricated using wire β arc additive manufacturing technique—A review. *J. Mater. Res. Technol.* **2021**, *13*, 754–768.

40. Mendricky, R.; Fris, D. Analysis of the Accuracy and the Surface Roughness of FDM/FFF Technology and Optimisation of Process Parameters. *Tech. Gaz.* **2020**, *4*, 1166–1173.

41. Multi Jet Fusion printing tips and tricks. In *HP Guide*; HP: Barcelona, Spain, 2018; pp. 1–21.

42. Delfs, P.; Töws, M.; Schmid, H.-J. Optimized build orientation of additive manufactured parts for improved surface quality and build time. *Addit. Manuf.* **2016**, *12*, 314–320. [CrossRef]

43. Choren, J.A.; Heinrich, S.M.; Silver-Thorn, M.B. Young's modulus and porosity relationships for additive manufacturing applications. *J. Mater. Sci.* **2013**, *48*, 5103–5112. [CrossRef]

44. Lee, K.P.M.; Pandelidi, C.; Kajtaz, M. Build orientation effects on mechanical properties and porosity of polyamide-11 fabricated via multi jet fusion. *Addit. Manuf.* **2020**, *36*, 101533. [CrossRef]

45. Poszvek, G.; Wiedermann, C.; Markl, E.; Bauer, J.M.; Seemann, R.; Durakbasa, N.M.; Lackner, M. Fused Filament Fabrication of Ceramic Components for Home Use. In *Lecture Notes in Mechanical Engineering*; Springer International Publishing: Cham, Switzerland, 2021; pp. 121–139. ISBN 9783030627836.

**Disclaimer/Publisher's Note:** The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.