Author Contributions
Conceptualization, N.D.L., C.J.G. and K.S.; methodology, C.J.G. and N.D.L.; finite element analyses and parametric studies, I.R.; conversion of scanning file, G.K. and S.V.; orthopedic expertise, K.S. and D.G.; finite element expertise, C.J.G. and N.D.L.; optimization expertise, N.D.L. and G.K.; writing—original draft preparation, I.R. and G.K.; writing—review and editing, C.J.G. and N.D.L.; visualization, F.P., G.K., S.V. and I.R.; supervision, C.J.G. and N.D.L.; and project administration, N.D.L. All authors have read and agreed to the published version of the manuscript.
Figure 1.
(a) Original scanned spine base. (b) Final smoothen model.
Figure 1.
(a) Original scanned spine base. (b) Final smoothen model.
Figure 2.
(a) CAD representation of spinal brace. (b) Division of the brace into individual surfaces to facilitate assignment of loads and boundary conditions.
Figure 2.
(a) CAD representation of spinal brace. (b) Division of the brace into individual surfaces to facilitate assignment of loads and boundary conditions.
Figure 3.
(a) Schematic illustration of the forces of the tethering straps on the rear face of the brace. (b) Imposed displacements and boundary conditions in the numerical simulation of the brace.
Figure 3.
(a) Schematic illustration of the forces of the tethering straps on the rear face of the brace. (b) Imposed displacements and boundary conditions in the numerical simulation of the brace.
Figure 4.
(a) Surfaces of attachment points of the straps at the rear face of the spinal brace. (b) Surfaces of supporting pads of the spinal brace.
Figure 4.
(a) Surfaces of attachment points of the straps at the rear face of the spinal brace. (b) Surfaces of supporting pads of the spinal brace.
Figure 5.
Description of the interaction surfaces between the brace and the human body. (a) Inner sides of the surfaces of supporting pads. (b) Outer sides of the surfaces of the human body.
Figure 5.
Description of the interaction surfaces between the brace and the human body. (a) Inner sides of the surfaces of supporting pads. (b) Outer sides of the surfaces of the human body.
Figure 6.
Pressure–overclosure chart describing the law of the brace–body interface in the pressure direction.
Figure 6.
Pressure–overclosure chart describing the law of the brace–body interface in the pressure direction.
Figure 7.
Contact pressures on the surfaces representing the human body (First numerical model).
Figure 7.
Contact pressures on the surfaces representing the human body (First numerical model).
Figure 8.
Displacements of the brace: (a) in the X-axis, (b) in the Y-axis, (c) in the Ζ-axis (First numerical model).
Figure 8.
Displacements of the brace: (a) in the X-axis, (b) in the Y-axis, (c) in the Ζ-axis (First numerical model).
Figure 9.
Von Mises stress results (MPa) for brace with thickness: (a) t = 5 mm, (b) t = 4 mm, (c) t = 3 mm, (d) t = 2 mm.
Figure 9.
Von Mises stress results (MPa) for brace with thickness: (a) t = 5 mm, (b) t = 4 mm, (c) t = 3 mm, (d) t = 2 mm.
Figure 10.
Von Mises stress results (MPa) for material with modulus of elasticity: (a) E = 2340 MPa, (b) E = 1520 MPa, (c) E = 1000 MPa, and (d) E = 320 MPa.
Figure 10.
Von Mises stress results (MPa) for material with modulus of elasticity: (a) E = 2340 MPa, (b) E = 1520 MPa, (c) E = 1000 MPa, and (d) E = 320 MPa.
Figure 11.
Section thickness—Displacements chart describing the change in deformation with increasing the thickness of the section of the brace.
Figure 11.
Section thickness—Displacements chart describing the change in deformation with increasing the thickness of the section of the brace.
Figure 12.
Modulus of elasticity–displacements chart describing the change in deformation with increasing modulus of elasticity.
Figure 12.
Modulus of elasticity–displacements chart describing the change in deformation with increasing modulus of elasticity.
Figure 13.
Geometry and dimensions of tensile strength specimen in mm.
Figure 13.
Geometry and dimensions of tensile strength specimen in mm.
Figure 14.
(a) Specimen orientation and (b) inner geometry of 3D-printed specimens. The light and dark blue color presents the outer walls and the orange color the infill.
Figure 14.
(a) Specimen orientation and (b) inner geometry of 3D-printed specimens. The light and dark blue color presents the outer walls and the orange color the infill.
Figure 15.
(a) Specimen placed in the testing machine and (b) specimen placed in in the testing machine after the end of the experiment.
Figure 15.
(a) Specimen placed in the testing machine and (b) specimen placed in in the testing machine after the end of the experiment.
Figure 16.
True stress–true strain curves for test specimens of direction 1 [0°].
Figure 16.
True stress–true strain curves for test specimens of direction 1 [0°].
Figure 17.
True stress–true strain curves for test specimens of direction 2 [90°].
Figure 17.
True stress–true strain curves for test specimens of direction 2 [90°].
Figure 18.
Results of primary stresses on the spinal brace: (a) S11; (b) S22.
Figure 18.
Results of primary stresses on the spinal brace: (a) S11; (b) S22.
Figure 19.
Displacements of the brace: (a) in the X-axis, (b) in the Y-axis, (c) in the Ζ-axis (Second numerical model).
Figure 19.
Displacements of the brace: (a) in the X-axis, (b) in the Y-axis, (c) in the Ζ-axis (Second numerical model).
Figure 20.
Contact pressures on the surfaces representing the human body (Second Numerical model).
Figure 20.
Contact pressures on the surfaces representing the human body (Second Numerical model).
Table 1.
Max von Mises stresses for various values of the brace cross-section thickness.
Table 1.
Max von Mises stresses for various values of the brace cross-section thickness.
Max von Mises Stress | Section Thickness |
---|
Smax = 20.29 MPa | t = 5 mm |
Smax = 18.15 MPa | t = 4 mm |
Smax = 16.60 MPa | t = 3 mm |
Smax = 14.06 MPa | t = 2 mm |
Table 2.
Max displacements in the three axes for various values of the brace cross-section thickness.
Table 2.
Max displacements in the three axes for various values of the brace cross-section thickness.
Section Thickness | Displacements |
---|
X-Axis | Y-Axis | Z-Axis |
---|
t = 5 mm | 13.65mm | 25.81 mm | 14.74 mm |
t = 4 mm | 14.28 mm | 26.53 mm | 16.21 mm |
t = 3 mm | 15.11 mm | 27.41 mm | 18.01 mm |
t = 2 mm | 16.32 mm | 27.64 mm | 20.35 mm |
Table 3.
Max von Mises stresses for various values of modulus of elasticity of the brace material.
Table 3.
Max von Mises stresses for various values of modulus of elasticity of the brace material.
Max von Mises Stress | Modulus of Elasticity E |
---|
Smax = 40.36 MPa | E = 2340 MPa |
Smax = 26.91 MPa | E = 1520 MPa |
Smax =18.15 MPa | E = 1000 MPa |
Smax = 6.17 MPa | E = 320 MPa |
Table 4.
Max displacements in the three axes for various values of the modulus of elasticity of the brace material.
Table 4.
Max displacements in the three axes for various values of the modulus of elasticity of the brace material.
Modulus of Elasticity, E | Displacements |
---|
X-Axis | Y-Axis | Z-Axis |
---|
E = 2340 Mpa | 14.13 mm | 25.00 mm | 15.40 mm |
E = 1520 Mpa | 14.21 mm | 25.83 mm | 15.82 mm |
E = 1000 Mpa | 14.28 mm | 26.53 mm | 16.21 mm |
E = 320 Mpa | 14.44 mm | 28.25 mm | 17.23 mm |
Table 5.
Summary of tensile test results.
Table 5.
Summary of tensile test results.
Specimen | Modulus of Elasticity [MPa] | Ultimate Tensile Stress [MPa] | Ultimate Strain |
---|
1.1 | 3127.1 | 47.06 | 0.0351 |
1.2 | 2857.4 | 48.49 | 0.0394 |
1.3 | 3004.0 | 48.36 | 0.0469 |
1.4 | 2988.4 | 48.97 | 0.0505 |
1.5 | 3241.4 | 50.42 | 0.0223 |
3.1 | 3121.2 | 35.76 | 0.0217 |
3.2 | 2871.3 | 32.82 | 0.0131 |
3.3 | 2887.7 | 37.46 | 0.0200 |
3.4 | 3089.4 | 33.10 | 0.0297 |
3.5 | 2518.7 | 34.17 | 0.0152 |
Table 6.
Descriptive statistics of 3D-printed PLA specimens.
Table 6.
Descriptive statistics of 3D-printed PLA specimens.
Orientation | Measure | Range | Min | Max | Mean | St. Dev | Confidence Interval (Two Sided) |
---|
95% | 99% |
---|
Xmean> | Xmean< | Xmean> | Xmean< |
---|
0° | E [MPa] | 384.00 | 2857.40 | 3241.40 | 3043.66 | 146.10 | 2848.97 | 3238.35 | 2719.95 | 3367.37 |
σult [MPa] | 3.35 | 47.06 | 50.42 | 48.66 | 1.21 | 46.96 | 50.36 | 45.83 | 51.49 |
90° | E [MPa] | 602.50 | 2518.70 | 3121.20 | 2897.66 | 240.38 | 2592.19 | 3203.13 | 2389.75 | 3405.57 |
σult [MPa] | 4.63 | 32.82 | 37.46 | 34.66 | 1.94 | 32.31 | 37.01 | 30.76 | 38.57 |
Table 7.
Material properties of 3D-printed PLA simulation.
Table 7.
Material properties of 3D-printed PLA simulation.
Material | E1 | E2 | v12 | G12 | G13 | G23 |
---|
| MPa | MPa | - | MPa | MPa | MPa |
PLA | 3043.66 | 2898 | 0.330 | 1140 | 1140 | 1140 |