*2.5. Model Task*

A rectangular beam of 140 mm × 28 mm × 14 mm was used for the algorithm implementation. Eight-node hexahedral finite elements were used for the calculations. The kinematic loading of 1 mm was used in the numerical simulation. The length of the kinematic loading region was 20 mm. In Figure 2, the loading scheme is presented; the *Vcon* region is marked by a green color. The end faces of the beam were fixed.

**Figure 2.** Loading scheme; *U* is applied displacements; the green region is *Vcon* region.

The mechanical properties of acrylonitrile butadiene styrene were used for further production by additive manufacturing. So, Young's modulus was equal to 200 GPa, the shear modulus was equal to 71.5 GPa, and the Poisson ratio was 0.4. For [*σ*]inf, 10% of maximum von Mises stress in the construction was used. The stop condition was as follows:

$$\max(|\lambda\_i - \lambda\_{i-\prime}|, |\beta\_i - \beta\_{i-1}|) < \varepsilon \tag{16}$$

where *ε* was equal to 10−<sup>3</sup> .
