1.1 Patient selection criteria
Inclusion criteria: (1) patients who can be consistently followed up post-surgery; (2) individuals without a history of hip preservation failure within 6 months after the surgical procedure; (3) participants who have not undergone any additional treatments during the follow-up period; (4) subjects who have undergone pre- and post-operative CT examinations, as well as periodic imaging assessments throughout the follow-up duration.
Exclusion criteria: (1) patients who did not comply with the follow-up; (2) patients experiencing severe adverse reactions, allergies, or accidents; (3) individuals with incomplete imaging data.
This study included a total of 33 patients (33 hips) who underwent rotational osteotomy of the base of the femoral neck at our hospital between January 2017 and December 2021.
Ethical approval was obtained from the Medical Ethics Committee of the Affiliated Hospital of Nanjing University of Traditional Chinese Medicine (2023NL-058-02).After providing all the enrolled patients with detailed information about the surgical procedure, their written informed consent was obtained. All procedures are carried out in accordance with the “Helsinki Declaration” of the World Medical Association.
1.2 General information.
The study included a total of 33 participants, comprising 30 males (representing 30 hips) and 3 females (representing 3 hips), with an age range of 18-48 years (mean age: 31.8). Among them, there were 13 cases of left hip involvement, 12 cases of right hip involvement, and 4 cases of bipartite hips. The types of osteonecrosis of the femoral head (ONFH) observed in this study were as follows: idiopathic ONFH accounted for 15 cases (15 hips), hormonal ONFH accounted for 9 cases (9 hips), alcoholic ONFH accounted for7 cases(7hips), and traumatic ONFH accounted for2cases(2hips).
1.3 Surgical methods
After preoperative planning and administration of general anesthesia, the patient is positioned laterally for surgery. A lateral incision is made on the affected hip joint to expose and protect the insertion points of gluteus medius muscle and rectus femoris muscle. Osteotomy is performed at the greater trochanter with careful preservation of both upper and lower attachments. The anterior side of the femoral neck joint capsule is dissected to assess any damage to the glenoid labrum. Subsequently, dislocation of the femoral head allows observation of its anterior and posterior collapse status. To evaluate blood supply, a hole is drilled using a Cricket's pin in order to investigate perfusion within the femoral head. On the posterior side of greater trochanter, soft tissue surrounding external rotator muscle group is meticulously peeled off as a flap.
1.4 Observation index and follow-up
The TAD, Cal-TAD, and NSA measurements were conducted on immediate postoperative frontal and lateral radiographs by two independent observers who were attending orthopaedic surgeons. Both observers received a detailed briefing regarding the TAD, Cal-TAD, and NSA measurements[ 11-12 ].
1.5 Finite element analysis: an example of FNS
The CT data of a male patient, aged 30 and weighing 80kg, with left-sided alcoholic femoral head necrosis ARCO type III was collected. The patient had a history of heavy drinking for 5 years.
CT scan parameters included a layer thickness of 1mm and scanning range from the upper edge of the pelvis to 15cm below the lesser trochanter. The CT data was saved in DICOM format, along with an original data copy.
1.5.1 Main software
Mimics 21.0 software, Geomagic Studio 2017 software, Solidworks 2021 software, ANSYS 17.0 software.
1.5.2 Modelling process
The CT image in Mimics Research 21.0 was segmented based on a specific threshold to obtain an initial model of the femoral head. Subsequently, the initial model was imported into Geomagic 2017 for establishing a preliminary geometric representation. Various processes including noise reduction, removal of regional features, and frosting were performed to optimize the femoral head model. After surface fitting, the optimized model was further assembled using SolidWorks 2021 software to establish the final preoperative geometric representation (Figure 1). This representation was then imported into ANASYS 17.0 software for material assignment, body mesh and surface mesh generation, as well as finite element analysis.
1.5.3 Mesh Division
The meshing settings of the model: the number of cells and nodes after meshing are shown in Table 1.
Table 1: Number of nodes and cells in the preoperative finite element model
Sports event Rotational osteotomy model of the base of the femoral neck
Number of nodes 160456
Number of units 89005
1.5.4 Establishment of different NSA and FNS models after rotational osteotomy via the base of the femoral neck
The femur model was imported into the 3D engineering software SolidWorks 2021 as a STEP file. A cross-section was created at the base of the femoral neck, with the osteotomy surface perpendicular to the axis of the femoral neck. The necrotic area was repositioned outside of the weight-bearing region by rotating the osteotomy block forward along the axis of the femoral neck, resulting in merging of both proximal and distal osteotomy surfaces. Six distinct postoperative NSA models were established based on variations in osteotomy surfaces.
The internal fixation device model was developed using Johnson & Johnson's FNS data, comprising a femoral neck splice plate, a femoral neck power bar, and a femoral neck anti-rotation screw. The FNS model was assembled with the femoral osteotomy model and positioned according to clinical fixation methods to create the final model (Figure 2). Subsequently, ANASYS17.0 software was employed for material assignment, body mesh formation, surface mesh generation, and finite element analysis.
1.5.5 Material property assignment
The material parameters for cortical bone, cancellous bone, necrotic zone, and internal fixation were inputted into ANSYS 17.0(Detailed information is presented in Table 2), with corresponding material properties assigned to each component of the analysis. The FNS internal fixation device was fabricated using titanium alloy material [7-8].
1.5.6 Loading load and setting constraints
According to relevant biomechanical studies, the unilateral loading pressure is calculated as 2.5 times the product of body weight (M), gravity coefficient (G), and a constant value of 0.5N. In this study, loads of 920N, 1500N, and 2000N were applied to the weight-bearing area above the femoral head to investigate the biomechanical changes following rotational osteotomy under different stress conditions. The weight-bearing area was defined as an arc spanning 40° in both internal and external directions from the center of the femoral head, as well as an arc spanning 80° in the anterior-posterior direction relative to its center (Figure 3) [9]. The distal end of the femur was immobilized as a fixed interface.
1.5.7 Evaluation index
Examine the variations in maximum displacement values and stress concentration phenomena at the osteotomy surface across different models of rotational osteotomy via the base of the femoral neck.
1.6 Statistical methods
The data were statistically analyzed using SPSS 27.0 software and GraphPad Prism 8. The normality of the distribution of continuous variables was assessed using the Shapiro-Wilk test. For variables with a symmetrical distribution, mean (mean) and standard deviation (SD) were reported, while median and interquartile range (IQR) were used for non-normal distributions. Categorical data were presented as absolute numbers and percentages.
Statistical comparisons of continuous variables were assessed using the Student's t-test for normally distributed variables and the Mann-Whitney U test for non-normally distributed variables. Binary logistic regression analyses were conducted with TAD, Cal-TAD, and NSA as independent variables, while the occurrence of total hip arthroplasty for failed hip preservation was considered as the dependent variable. Intraclass correlation coefficients (ICC) along with 95% confidence intervals (CI) were calculated for interval data. The results of the logistic regression analysis were visualized using forest plots and Nomogram plots. The Nomogram plots were utilized to develop a scale based on the magnitude of regression coefficients from all independent variables, assigning scores to each level of every value for each independent variable. Subsequently, a total score was computed for each patient, followed by calculating the probability of outcome occurrence at that specific time through a transformation function between score and outcome probability.
The area under the ROC curve (AUC) was calculated to assess the performance of TAD, Cal-TAD, and NSA models and determine the optimal threshold predictive value. AUC values were interpreted as follows: AUC = 0.5 indicated no discrimination; 0.5 < AUC ≤ 0.7 represented low precision; 0.7 < AUC ≤ 0.9 indicated moderate precision; AUC > 0.9 suggested high precision; and an AUC of 1 denoted perfect discrimination. The optimal thresholds for TAD, Cal-TAD, and NSA were determined based on Youden's J statistic by identifying the point on the ROC curve with maximum distance from the diagonal line.
LASSO regression effectively mitigated overfitting by compressing the coefficients of variables in the regression model, leading to the generation of a Risk Score. Utilizing risk factor analysis, we examined the association between the Risk Score and both survival status and survival time, subsequently generating a heatmap illustrating the relationship between Risk Group and various risk factors.