Material Preparation and 3D Reconstruction: Our study commenced with the acquisition of high-precision CT images of the femur, serving as the foundational data for constructing subsequent three-dimensional finite element models. Employing MIMICS, a specialized medical image processing software, we meticulously extracted the geometric information of the femur from these images, culminating in a detailed 3D reconstruction.
Model Optimization: Post-reconstruction, models often exhibit imperfections such as sharp edges, overlapping surfaces, and voids. To rectify these issues, we utilized GEOMAGIC software to optimize the model, focusing on smoothing sharp edges, eliminating overlapping surfaces, and filling voids. This process ensured the geometric integrity of the model, closely mimicking the authentic femoral geometry and laying a robust foundation for subsequent finite element analysis.
Finite Element Analysis: Following optimization, the enhanced IGES file was imported into ANSYS for three-dimensional finite element analysis. Within the simulation, two distinct material properties were defined: cortical bone and trabecular bone. Cortical bone was assigned a Young's modulus of 17,000 MPa and a Poisson's ratio of 0.3; trabecular bone had a Young's modulus of 50 MPa and a Poisson's ratio of 0.4. For element selection, we consistently used Quadratic Tetrahedral elements, with an element size set at 5 mm.
Experimental Design: We devised multiple experiments to evaluate the impact of three independent variables—traction duration, traction force[3,4], and traction angle[5]—on the stress characteristics and deformation of the femur. Our experimental plan was as follows:
Investigating Traction Force and Duration:
Established three traction force levels (360N, 400N, 500N) representing insufficient, adequate, and excessive traction forces, respectively.
Analyzed the variation in stress values across different traction forces and determined the maximum tolerable traction duration for each force level.
This exploration aimed to elucidate the relationship between traction force and duration, guiding the dynamic adjustment of traction force over time.
Utilized least squares linear regression analysis to validate the correlation between traction force, maximum tolerable duration, and stress values.
Traction Angle Variability:
Maintaining constant traction force and duration, we adjusted the traction angle in a two-dimensional space.
Examined the stress features and deformation across the femur under varying traction angles.
Employed least squares linear regression analysis to confirm the association between angle and stress values.
We also do theoretical analysis and derivation for offset in three-dimensional space, and the schematic diagram is shown in Fig. 2:
To begin with, let L1 represent the line of normal traction force, and L1′ denote the line of traction force after deviation. Next, select point P on L1′, and construct plane n parallel to L1 passing through P, along with plane m perpendicular to plane n. Then, project L1 onto plane m to obtain L2 within the same plane, where the angle between L2 and L1′ is denoted as θ. Following that, project L2 onto plane n to yield L3, and the angle between L3 and L2 is designated as ϕ. The angle between L1′ and L3 is noted as ψ. Our goal is to solve for angles θ and ϕ so that compensations can be made in the corresponding directions.
The detailed steps for solving these angles are as follows:
Choose point Q on L1, measure the distance between points P and Q, denoted as d′.
Measure the distance from point Q to plane n, designated as d, and the angle between d and d′ is recorded as α.
Determine the distance between the perpendicular from point P to L3 and point V, referred to as d′′.
Using the Law of Cosines, d2 + d′2 − d′′2 = 2dd′cos(α), calculate d′′.
Let W be the intersection point of L1′ ,L2, and L3. Calculate angle θ using the distances dwv (from W to V), dwp (from W to P), and d′′.
Applying the Law of Cosines again, dwv2 + dwp2 − d′′2 = 2dwvdwpcos(θ), find θ.
Utilizing the Triple Cosine Formula, determine the product of cos(θ) and cos(ϕ).
Measure the distance from point P to plane m. Since L2 is the projection of L1′ onto plane m, the foot of the perpendicular from L1′ to plane m will fall on L2, marked as X.
Measure dpx and dwx at points P and W, respectively.
In triangle pwx, applying the Law of Cosines, calculate β.
Finally, derive ϕ from the calculated values.
Compensate for the angles θ and ϕ to re-establish parallelism between the lines of traction force and fracture force.
Integrating the findings from the aforementioned two sets of experiments, we explored the interrelationship between traction weight, traction time, and traction angle. Our focus was on investigating the impact of dynamic adjustments in traction weight and time, along with the effect of maintaining parallel traction angles on stress characteristics. The results substantiated that dynamic parallel traction yields the most favorable outcomes. And Table 1 shows the attribute setting of 3D finite element analysis.
Table 1
Settings for Material Properties
Part | Young’s modulus(E)[MPa] | Poisson’s ratio(v) | Element type | Element Size |
Cortical bone | 17000 | 0.3 | Quadractic tetrahedral | 5mm |
Cartilage | 50 | 0.4 | Quadractic tetrahedral | 5mm |