Design Data and Finite Element Analysis of 3D Printed Poly( ε -Caprolactone)-Based Lattice Scaffolds: Inﬂuence of Type of Unit Cell, Porosity, and Nozzle Diameter on the Mechanical Behavior

: Material extrusion additive manufacturing (MEAM) is an advanced manufacturing method that produces parts via layer-wise addition of material. The potential of MEAM to prototype lattice structures is remarkable, but restrictions imposed by manufacturing processes lead to practical limits on the form and dimension of structures that can be produced. For this reason, such structures are mainly manufactured by selective laser melting. Here, the capabilities of fused ﬁlament fabrication (FFF) to produce custom-made lattice structures are explored by combining the 3D printing process, including computer-aided design (CAD), with the ﬁnite element method (FEM). First, we generated four types of 3D CAD scaffold models with different geometries (reticular, triangular, hexagonal, and wavy microstructures) and tunable unit cell sizes (1–5 mm), and then, we printed them using two nozzle diameters (i.e., 0.4 and 0.8 mm) in order to assess the printability limitation. The mechanical behavior of the above-mentioned lattice scaffolds was studied using FEM, combining compressive modulus (linear and nonlinear) and shear modulus. Using this approach, it was possible to print functional 3D polymer lattice structures with some discrepancies between nozzle diameters, which allowed us to elucidate critical parameters of printing in order to obtain printed that lattices (1) fully comply with FFF guidelines, (2) are capable of bearing different compressive loads, (3) possess tunable porosity, and (3) overcome surface quality and accuracy issues. In addition, these ﬁndings allowed us to develop 3D printed wrist brace orthosis made up of lattice structures, minimally invasive (4 mm of thick), lightweight (<20 g), and breathable (porosity >80%), to be used for the rehabilitation of patients with neuromuscular disease, rheumatoid arthritis, and beyond. Altogether, our ﬁndings addressed multiple challenges associated with the development of polymeric lattice scaffolds with FFF, offering a new tool for designing speciﬁc devices with tunable mechanical behavior and porosity.


Introduction
Highly controlled lattice architectures are seeing growing interest in different application fields, ranging from the manufacturing of lightweight orthosis [1][2][3], scaffolds for tissue engineering [4][5][6][7][8][9][10][11], implantable medical devices [12][13][14], stimuli-responsive softrobotics [15][16][17], and topology optimization for architectural structures [18]. Such lattice structures comprise a base form, namely unit cell, which is spatially repeated across the entire structure either in a fixed x, y, z geometry or with varying degrees of asymmetry [19][20][21][22][23]. From a design perspective, the key advantages of lattice structures include dimensional accuracy, surface quality, residual stresses, and microstructures, dictating not only the mechanical performance of the overall structure but also other critical features such as surface area, porosity, fluidic dynamics, and materials consumption [24,25]. Therefore, Eng 2022, 3 control over these design parameters is of fundamental interest for the optimization and prediction of the behavior of these structures for specific applications.
Material extrusion additive manufacturing (MEAM), also known as 3D printing, is an advanced manufacturing method able to produce lattice structures (and beyond) via precise layer-wise addition of material through its extrusion from the millimeter-or micronscaled nozzle [26]. However, the potential of MEAM to prototype lattice structures is hampered by restrictions imposed by manufacturing processes that lead to practical limits on the form and dimension of structures that can be produced. Indeed, as previously reported by Panetta and co-workers [27], such structures are mainly manufactured by selective laser sintering (SLS) or stereolithography (SLA) [28] since MEAM does not have the adequate resolution to precisely print microstructures. Nevertheless, limitations such as minimum cell size, thickness of walls/edges/faces, cell geometry, nozzle diameter, printing parameters, or porosity have not been fully evaluated. On the other hand, as highlighted by Silva et al., [28], albeit MEAM could seem inadequate compared to SLS and SLA for the production of lattice structures, it allows to use biocompatible and bioabsorbable polymers more suitable for biomedical applications, the feedstock is safer and easier to handle and does not require additional post-processing.
For these reasons, in this study, we aimed to investigate the capability of MEAM to produce custom-made lattice structures, combining fused filament fabrication (FFF), including computer-aided design (CAD), with the finite element method (FEM) to assess their printability, tune the porosity, and predict their mechanical behaviors, useful for tissue engineering and biomedical implants.
We focused our attention on generating and examining four types of 3D CAD scaffold models with different geometries (reticular, triangular, hexagonal, and wavy) as the unit cells. In addition, paying attention to the available polymers for 3D printing purposes, we used poly(ε-caprolactone) (PCL) because of its low melting point of 60 • C, high printability, soft-like properties, biocompatibility, and ease of multi-functionalizations with bioactive cues (i.e., biomolecules, proteins/peptides, growth factors, and nano-composites) [29].
Each PCL-based lattice structure was 3D printed with different unit cells (ranging from 1 to 5 mm) using two nozzle diameters (i.e., 0.4 and 0.8 mm) in order to assess the printability limitation and the scaffolds manufacturability in terms of minimum unit cell size, and pore size compared to the nozzle diameter. Further, the mechanical behavior of such lattice scaffolds was studied using FEM, combining compressive modulus and shear modulus either in linear or nonlinear hyperelastic models. The CAD-based FEM method allowed us to predict and understand the mechanical properties of the engineered scaffolds vs. the nozzle diameter, thus providing tools for the design and control of the geometric configuration of the lattice structures.
Ultimately, this study aims to address the design complexity issues for lattice structures creation by MEAM, realizing a CAD/FEM-driven methodology that could be a useful tool for designers, materials scientists, and MEAM engineers who wish to apply complex structures for advanced designs, where the techniques described could be readily applied for specific applications, such as light-weighting orthosis, scaffolds for targeted tissues or 3D cell cultures, controlled drug-release, medical devices, and soft robotics.
We hypothesize that the controlled biomechanics and porosity of lattice scaffolds fabricated by 3D printing could be a useful tool obtained through designing different inner geometries via the FEM method.

Materials
Poly-ε-Caprolactone (PCL) (Filoalfa, Ozzero, Italy) was used as polymer filament without further modifications. All materials were handled with gloved hands and following standard surface analysis laboratory practices to minimize any possible contamination.

CAD-Based Lattice Structures Design
Four different lattice structures CAD designs were created using Autodesk Inventor ® (Autodesk 2020, McInnis Parkway San Rafael, CA, USA), simulating the superimposition of different paths during the 3D printing procedure. The reticular, triangular, hexagonal, and wavy geometries were designed making the angle between different paths variables. In order to create the structures to be tested, each pattern was considered as a linear repetition of a unit cell size; five cell sizes were designed and tested (ranging from 1 to 5 mm). All designed scaffolds featured 5 × 5 cell units in the x and y directions, with 15 layers of thickness for triangular and hexagonal geometries and 16 layers of thickness for reticular and wavy ones. The porosity of the CAD models was obtained by measuring the occupied volume of each scaffold, and the designed porosity (P CAD ) was calculated as follows (Equation (1)): in which V f is the occupied volume of the scaffold, and V tot is the bulk volume of the CAD model.

Finite Element Method (FEM) Modeling
In order to predict the mechanical behavior of designed lattice scaffolds and to assess the dependency of the mechanical properties with respect to the printing path width, the compressive and shear moduli were simulated by finite element analysis using ABAQUS/standard ® (Dassault Systemes Simulia Corporation., Johnston, RI, USA). Each CAD model was imported into the FE software as a continuous unique part. The PCL-based scaffolds properties were simulated by both linear and nonlinear elastic behavior. In the linear elastic condition, the elastic modulus and Poisson's ratio of PCL were assumed to be 350 MPa and 0.3, respectively, as reported previously [30,31]. The elastic modulus was then used in the FE simulation with the assumption of linear elastic deformation and considering the PCL material as isotropic. In the STEP section of the model, the Nlgeom option was set to on in order to allow results with nonlinear behavior. An automatic incrementation was selected for the simulation, with a maximum number of increments equal to 100 and an increment size increasing from an initial size of 0.01 up to 0.1 to speed up the process. No particular problems in the iterations were registered. The minimum increment size was left to a default 1.0 × 10 −5 value. The output requested from the simulations contains stresses, strains, displacement, and Force/Reaction to allow us to calculate the related compressive modulus. These variables were calculated every n number of increments. Then, compression simulation was performed on the CAD models constraining the lower face of the scaffold, denying any possible movement or rotation with respect to the reference system, and 1.5% compressive strain was exerted to the top of the model. The compressive elastic modulus of each scaffold was obtained from Equation (2): in which σ comp is the equivalent maximum stress reached for the imposed displacement, and ε is the imposed strain. Since: where R is the reaction force of the constrained lower side of the scaffold, and A eq is the resistant equivalent area of the scaffold (see Figure S1 for further details). The imposed strain was calculated from Equation (4) as follows: where δ is the imposed displacement, and t is the thickness of the scaffold (see Figure S1 for further details). Hence, the final formulation of the compressive elastic modulus, taking into account all the boundary conditions already mentioned, is: Afterward, shear moduli of scaffolds with four types of microstructures were estimated. The directions of the scaffold's length were called x-direction, while the in-plane orthogonal direction was called y-direction. A displacement parallel to the xand ydirection was imposed to the upper surface of each scaffold, and the linear shear modulus was then calculated with Equation (6): where τ is the shear stress of the homogenized lattice structure for the maximum imposed displacement, and γ is the shear strain of the simulation. The shear stress can be then obtained using: where R x and R y are the reaction forces in the xand y-direction of the constrained portion of the tested scaffold, while A eq is the resistant equivalent area of the scaffold. The shear strain instead is calculated as follows: where δ is the applied displacement of the upper part of the scaffold and t is the total thickness of the scaffold as calculated before. Hence, by combining the previous equations, it is possible to obtain the final shear modulus formulation (Equation (9)): For nonlinear elastic behavior, the compressive mechanical behavior of each scaffold was predicted by using the nonlinear hyperelastic van der Waals material model [31,32]. For this model, the strain energy function is: where U is the strain energy potential, J el is the elastic volume ratio, I 1 and I 2 are the first and second invariants of the deviatoric strain, µ, λ m , α, and β describe the deviatoric behavior of the material, and D introduces the compressibility [33].
Since in van der Waals hyperelastic material model, µ 0 (initial shear modulus) is equal to µ (shear behavior of the material in van der Waals), D value can be calculated as follow: The utilized parameters for the van der Waals hyperelastic model for describing the nonlinear mechanical behavior of PCL were µ = 140, λ m = 1.9, α = 6, β = 0, and D = 0.00659, as reported previously [31].

Fused Deposition (FDM) Printing
All lattice structures were 3D printed using fused deposition modeling (FDM) technology through QXXL Full-Metal (Sharebot, Nibionno, Italy). PCL filaments (1.75 mm diameter) were extruded through a 0.4 mm and 0.8 mm diameter nozzle at the temperature of 130 • C. The bed temperature was set at 40 • C. The layer height was 0.4 mm, and the nozzle travel rate was 3600 mm min −1 . Each object was modeled using Autodesk Inventor ® software, then exported in stereolithography (.stl) format and sliced using Simplify3D ® software to be printed. To ensure the adhesion of the scaffolds to the surface of the printing platform, a 3D glue stick (Magigoo™, Swieqi, Malta) was used.

Spectroscopic Characterizations of 3D Printed Lattice Structures
Raman analysis was conducted using Progeny™ spectrometer (Rigaku Corporation, Tokio, Japan) with a 1064 nm laser source and selectable laser output set at 490 mW. The spectral range is 200-2500 cm −1 with transmission-type volume phase grating. The spectral resolution is 15-18 cm −1 , and the detector is a thermoelectrically cooled indium gallium arsenide (InGaAs) [34]. The samples were analyzed with 60 cumulative scans with optimized laser power, aperture size, and duration (7 s) per exposure in order to achieve the best signal-to-noise ratio. All of the obtained spectra data were plotted using Origin™ 8 software (One Roundhouse Plaza, Northampton, MA, USA).

CAD-Based Scaffolds Design and Printability
Once unit cell configurations have been generated as a surface, the first procedure was to convert this into a three-dimensional form (see Figure S1 and Methods for further details). For triangular and hexagonal geometries, the scaffold consisted of 5 × 5 unit cell repetitions of 15 layers, while the reticular and wavy scaffolds consisted of 5 × 5 cell unit repetitions of 16 layers (Figure 1). Then such CAD models were 3D printed using Poly-ε-Caprolactone (PCL) as feedstock material in order to evaluate the printability of each scaffold by changing the nozzle diameter as a parameter of interest.
Starting from the 0.4 mm nozzle diameter, the minimum unit cell sizes allowed for printing were 2 mm for reticular and triangular geometries and 1 mm for hexagonal and wavy ones. Instead, the maximum value of unit cell achievable from such nozzle diameter was equal to 4 mm for reticular and triangular geometries and 3 mm for hexagonal and wavy ( Figure S2A-D). In these unit cell ranges, the overall accuracy and quality of printed structures were optimal. On the contrary, a smaller unit cell size compared to the previously mentioned values led to a fusion of such scaffolds where the geometry can hardly be distinguished, and the shape and thickness of struts were inconsistent. Larger unit cell size strongly compromised the inner scaffold structure, leading to large voids generated from lack of material due to uncontrolled filament deposition during printing: extruded filament could not secure adhesion to the previous layers, and the printing failed.
For the lattice structures, 3D printed using a 0.8 mm nozzle, only larger unit cell sizes for each pattern can be printed successfully ( Figure S2E-H). In particular, we were able to print unit cell sizes ranging from 3 mm up to 5 mm for reticular and triangular geometries and 2-4 mm for wavy and hexagonal structures. Because of a wider printed path, the manufacturability of smaller unit cell sizes was hampered, as this resulted in the scaffold structures melting. Overall, with the use of either 0.4 mm or 0.8 mm nozzle diameter, it was possible to tune the unit cell sizes in order to obtain stable or well-defined lattice scaffolds that fully comply with FFF guidelines, avoiding stringing inside unit cells and lumps of filaments.
However, for the latter, we must keep in mind the extruder type (i.e., direct drive extruder and Bowden drive extruder) that could influence the outcomes of the overall structure. Briefly, the direct drive extruder has a direct drive gear close to the heating system and nozzle; hence the extrusion and retraction commands are more immediate, creating a more stable material deposition control and better resolution [35]. Contrary, in the Bowden drive extruder, the drive gear is installed outside the extruder head [35].
In order to handle such lattice structures based on the type of extruder available, the retraction settings must be adjusted to control the flow and the correct filament deposition. In this case, adopting the Sharebot QXXL 3D printer (see Materials and Methods for further details), we had a direct drive extruder, and the optimal retraction distance and retraction speed were 0.5 mm and 4200 mm min −1 , respectively. In the case of using Bowden drive extruder, with PCL as feedstock material, for obtaining the same accuracy in the deposition of filament (and therefore the same final lattice scaffolds), high values of retraction were needed to prevent unnecessary extrusion while the extruder travels between printed sections. Specifically, such settings must be of 6 mm and 4800 mm min −1 for retraction distance and retraction speed, respectively (data not shown).
Next, we examined the porosity of CAD-based lattice scaffolds by varying the unit cell types as it can affect the topology of the entire structure, the fluid flow, cell ingrowth, vascularization, nutrition/oxygen supply, and energy absorption. Analyzing the data, it can be clearly seen for all scaffold architectures that the porosity increases linearly as the unit cell size increases, as would be expected (Figure 2). In addition, for the same unit cell size, the porosity of the 3D printed scaffolds with the 0.8 mm nozzle is always lower than those with the 0.4 mm. It should be emphasized that all scaffolds are highly porous, ranging from 46% (for wavy scaffold with 1 mm unit cell) to 87.50% (for triangular scaffold with 4 mm unit cell). Overall, these data demonstrate a relationship between unit cell size and lattice structures, which is useful for fine-tuning the properties of scaffolds for tissue engineering, 3D cell cultures, bio-robotics, and orthopedics, but also to select which structure and which unit cell is most suitable for a specific use. It is further believed that by fine-tuning the lattice CAD models and/or adjusting the base unit cell size, a greater number of effective lattice gradients can be realized, allowing increased flexibility while adjusting the porosity for a given scaffold thickness.

FEM Analysis
For linear FEM analysis, the elastic modulus of PCL filament was considered to be equal to 350 MPa, and the Poisson's ratio was set to 0.3 (see Materials and Methods).
A simple mesh sensitivity analysis was performed in order to obtain a plateau value of the compression and shear modulus for each tested scaffold. The FE model was generated with a mesh size of 0.5 mm, thus obtaining a number of finite elements of the order of 10 5 . The elastic modulus was used in the FEM models with the assumption of linear elastic deformation, the FEM analysis was performed for each scaffold, and the mechanical behavior under compression was simulated. The distribution of von Mises stress for comparing the scaffolds with different geometries and nozzle diameters is shown in Figure 3. It is interesting to note that in the reticular, triangular, and wavy scaffolds, when subjected to stress, we have noticed the increase in the compressive elastic modulus due to the existence of thicker load columns. Conversely, the amount of maximum stress for the hexagonal scaffolds was lower than those of other scaffolds due to the absence of proper supporting columns. We assumed that the design parameters could influence the bonding strength between layers and filaments in the printed scaffolds, and the layer bonding of the hexagonal scaffold could be significantly weaker compared to other tested scaffolds because in every three layers, their contact points shifted, which could result in a decrease of the bonding strength. Using the FEM simulation, the amount of force after completion of analysis was obtained, and the equivalent compressive elastic modulus was calculated by the ratio of compressive stress to compressive strain (see Equation (2) in the Materials and Methods).
The elastic modulus calculation showed how the mechanical properties of each scaffold could be strongly influenced by nozzle diameter and unit cell sizes (Figure 4). The compressive elastic moduli of the scaffolds with different geometries were predicted to have a broad range from 0.07 to 116.48 MPa. This range can be controlled by the design parameters such as unit cell size, geometries, porosity, nozzle diameter, and feedstock material. Considering a similar pattern-with the same unit cell size-it is possible to notice an increasing value of compressive modulus for the nozzle 0.8 mm compared to that of 0.4 mm. Instead, for the same path width, by increasing the unit cell sizes, the compressive features tend to decrease, thus weakening the scaffolds.
Leaving out the highest compressive moduli of wavy and hexagonal geometries at 1 mm of unit cell size (due to the lower void content in the scaffolds), taking into consideration the unit cell of 2 mm, it is highlighted that the more rigid structure is guaranteed by the reticular geometry (33.92 MPa) while the triangular and wavy patterns possess similar moduli (27.61 MPa and 28.8 MPa, respectively); this difference between patterns is maintained even with larger unit cells and larger print paths (Table 1). Rather, the hexagonal pattern displayed compression moduli of an order of magnitude lower than the other geometries, and this, as mentioned above, is due to the absence of thicker load columns.  In addition, the shear modulus of each microstructure was estimated by simulating shear tests on each scaffold in two lateral directions (Figures 5 and S3). We observed that the values of shear moduli for reticular and wavy scaffolds have no directional dependence, probably due to the symmetry in their geometries. On the contrary, the triangular and hexagonal scaffolds have an uneven elastic modulus in lateral directions, with the elastic moduli slightly different from xto y-direction (Table S1). Altogether, shear moduli were dependent on both the type of geometries and unit cell sizes. Moreover, shear moduli of reticular and wavy scaffolds were similar in xand y-direction, showing the highest shear moduli to compared other structures due to their supporting load columns; on the contrary, triangular and hexagonal scaffolds have the highest shear modulus in y-direction; on the contrary, the decreased value in x-direction was probably because of the highest amount of inter-layer contact points in these directions.
Lastly, we sought to estimate the compressive behavior of such lattice scaffolds via nonlinear FEM analysis ( Figure 6A-D), using the Van der Waals hyperelastic model (see Materials and Methods for further details), as this model would better describe the material behavior of PCL filament [31]. Here, the compressive moduli of different scaffolds were calculated considering the initial and final point of the simulation, for a strain of 5%, since for all models, the compressive moduli were similar in small strain ranges but, by increasing the strain value, the stress-strain curves tend to diverge, showing a softening ( Figure 6E). The results showed that minimum variations of compressive modulus are registered for hexagonal geometries, where a difference close to 0 has been observed for all unit cell sizes. Higher softening is registered for reticular, wavy, and triangular geometries, with a variation ranging from 10 to 18%. Additionally, in this case, the reticular scaffolds resulted the stiffer, followed by wavy and triangular, whereas the hexagonal scaffolds remained the softer structures, maintaining almost the same compressive moduli with respect to the linear model (Table 2). Overall, the nonlinear results showed that all PCL-based lattice structures with different unit cell sizes demonstrated hyperelastic behavior under compression. Nevertheless, we are aware that to better understand the nonlinearities in such structures would be useful to carry out predictive FEM analysis both with various load models and with higher strain ranges, coupled with experimental data.

Spectroscopy Characterization of 3D Printed Scaffolds
To explore if printing parameters could influence the chemical structures of PCL-based lattice scaffolds, which play an important role in the scaffolds overall behavior (such as biodegradation rate, crystallinity, and morphology), we carried out Raman spectroscopy analysis ( Figure 6F). Raman spectra of lattice scaffolds printed with 0.4 and 0.8 mm nozzle diameter exhibited a similar peak pattern with minimal intensity changes at 915 and 961 cm −1 (in blue) associated with C-COO stretching, and three narrow bands at 1033 cm −1 , 1066 cm −1 , and 1110 cm −1 (in light green) associated with C-C stretching. Some other bands are visible approximately at 1284 and 1301 cm −1 (in magenta) due to the symmetric CH 2 stretching, and at 1419 cm −1 , 1443 cm −1 , and 1468 cm −1 (in purple) ascribable to CH 2 bending; lastly, a peak was observed at 1722 cm −1 (in orange) due to C=O stretching vibration.
In summary, this spectroscopic characterization seems to suggest that the printing parameters do not alter the chemical structures of such lattice scaffolds.

3D Printed Wrist Brace Orthosis with Different Lattice Geometries
Finally, we wondered if these results allowed us to go beyond the manufacturability of scaffolds for tissue engineering. The results obtained enabled us to develop 3D printed wrist and hand brace orthosis made up of lattice structures that might be minimally invasive (but stiffer enough to support the limb), lightweight, and breathable (due to the possibility of fine-tuning the porosity of such structures). In recent years, 3D printing has offered an innovative method for rapid and customizable prototyping in the biomedical field, making these advantages possible. However, it is seldom discussed in the fields of hand rehabilitation, and many clinicians and hand therapists are still unaware of how this technology could be a useful and effective tool in the rehabilitation of patients with neuromuscular disease, rheumatoid arthritis, and beyond [36]. Therefore, to shed some light on the grey areas, we envisioned how to engineer novel orthosis with lattice geometries (Figure 7). Since reticular and triangular geometries provided better results both in terms of mechanical properties and porosity, we have chosen them to design the wrist and hand orthosis. Briefly, the 3D printed lattice structure orthosis has been created using Autodesk Inventor ® . The procedure to obtain the final device was divided into two different phases: after the measurement of the anatomic region of interest, all data are collected to draw the first representation of the orthosis. Then, via extrusion, the edges of the orthosis and the reinforced portions are created. Subsequently, linear repetition of lattice geometry was created and cut to fit the inner boundary of the drawn orthosis. Finally, combining the boundary with the selected pattern, it is possible to export it in stereolithography (STL) format and slice it using Simplify3D ® software. The final 3D printed orthosis was 4 mm thick, weighed <20 g, and possessed breathability >80%.
Surprisingly, the combination of PCL as feedstock material and lattice geometries with tunable mechanical features allowed us to obtain shape-memory orthosis. Such property is of particular interest as by using these structures coupled to PCL, it is no longer necessary to thermoform the polymer (as required for PLA, for instance), but on the contrary, PCL-based orthosis immediately takes the shape of the limb and gets back to its original shape once removed from the limb. This proof of concept gives us the impetus to test and validate these orthoses in future works.

Conclusions
In conclusion, we reported the design data and printing parameters for 3D printing of functional PCL-based lattice structures that fully comply with FFF guidelines are capable of bearing different compressive and shear loads, and with tunable porosity.
Based on the 3D CAD models, the scaffolds with four types of designed geometries (reticular, triangular, hexagonal, and wavy) were fabricated through FDM, elucidating the key value of the nozzle diameter and unit cell sizes in order to obtain stable or well-defined lattice scaffolds avoiding stringing inside unit cells and lumps of filament.
Further, we demonstrated that the FEM combined with 3D printing has the potential to predict the mechanical features of the scaffolds either in linear or nonlinear models, useful to select the most appropriate geometry, unit cell size, and nozzle diameter for a specific need without any fabrication, thus offering a new tool for designing specific devices with tunable mechanical behavior and porosity for tissue engineering, 3D cell cultures, controlled drug-release, medical devices, soft robotics, and beyond the usefulness of the PCL-based 3D printed lattice scaffolds for the fabrication of stiffer, lightweight, and breathable wrist and hand brace orthosis is demonstrated.
Supplementary Materials: The following are available online at https://www.mdpi.com/article/10 .3390/eng3010002/s1, Figure S1: Explanation picture about the quantities reported in Equations (2) and (3), Figure S2: 3D printed PCL-based geometries with different nozzle diameters, Figure S3: The von Mises stress-strain distributions in the simulated shear test by FEM in the y-directions, Table S1: Comparison of shear moduli of four types of scaffolds with different nozzle diameter and unit cell sizes.
Funding: This research received no external funding.

Data Availability Statement:
The authors confirm that the data supporting the findings of this study are available within the article [and/or] its Supplementary Materials.