Skip to main content

Biomechanical effect of intertrochanteric curved varus osteotomy on stress reduction in femoral head osteonecrosis: a finite element analysis

Abstract

Background

Intertrochanteric curved varus osteotomy (CVO) has been widely used to remove the necrotic bone away from the weight-bearing portion in the treatment of osteonecrosis of the femoral head (ONFH). However, whether all types of necrosis will benefit from CVO, in terms of the stress level, the effect of different center-edge (CE) angles of acetabulum on stress distribution of necrosis after CVO, and the relationship between the intact ratio and the stress of necrosis, has never been addressed. The purpose of the study was to evaluate the influence of CVO on the stress reduction in necrotic bone using a finite element analysis (FEA) with different CE angles.

Methods

CVO finite element models of the hip joint were simulated with a lesion of 60°. The osteotomy angles were divided into four configurations (15°, 20°, 25°, and 30°), and three types (A, B, and C1) of lesions were established based on the Japanese Investigation Committee (JIC) classification. In addition, two CE angles (18° and 33°) of acetabulum were considered. The maximum and mean von Mises stress were analyzed in terms of the necrotic bone by a physiological loading condition. Moreover, the correlation of the intact ratio measured in 3D and the stress distribution after CVO was analyzed.

Results

Stress reduction was obtained after CVO. For type B, the CVO angle was 20° (0.61 MPa), and for type C1, the CVO angle was 30° (0.77 MPa), if the mean stress level was close to type A (0.61 MPa), as a standard. The maximum and mean von Mises stress were higher in the CE angle of 18°models, respectively. The intact ratio measured in 3D had a good negative correlation with stress after CVO and had more influence on stress distribution in comparison to other geometric parameters.

Conclusions

For making decisions about the biomechanics of CVO, a CVO angle of > 20° was recommended for type B and > 30° was safe for type C1. The risk of progressive collapse was increased in the insufficient situation of the weight-bearing portion after CVO. The intact ratio could provide information about clinical outcomes and stress distribution after CVO.

Background

Osteonecrosis of the femoral head (ONFH) is a progressive pathological process that occurs due to multiple factors affecting the blood supply of femoral head [1]. Among the preservation procedures for the treatment of ONFH, intertrochanteric curved varus osteotomy (CVO) is a proximal femoral osteotomy technique that is widely used in Japan [2,3,4,5,6,7,8,9,10,11,12]. It is performed with the aim of removing necrotic regions of the femoral head from the weight-bearing portion, decreasing the load subjected to infarction in order to prevent collapse in the early stage. The key point of CVO is to create a curved varus osteotomy line though the tip of great trochanter and the middle of the lesser trochanter without exceeding the intertrochanteric crest. The osteotomy center is not necessarily located in the femoral head center [13].

Clinically, to estimate the collapse risk after CVO, the intact ratio is used and defined as the intact surface area divided by the weight-bearing area in the femoral head. A post-operative intact ratio of 33.0% was necessary if a satisfactory result was to be achieved [14], whereas it was prone to collapse if the necrotic lesion was located in the anterolateral weight-bearing portion or with insufficient femoral head coverage [15]. Moreover, the post-operative intact ratio of CVO <33.3% and lateral center-edge (LCE) angle <25° were identified as independent factors leading to radiographic failure and conversion to total hip arthroplasty from clinical reports [16]. However, from a biomechanical viewpoint, few studies have focused on the stress changes of necrotic bone and the influence of different LCE angles on the stress distribution of necrotic bone after CVO. In addition, the correlation between the intact ratio and stress in the necrotic bone has not been discussed.

A finite element analysis (FEA) is a popular computer simulation for evaluating the stress and strain state in the human body, including studies on ONFH [9, 17,18,19,20]. Thus, the purpose of the present study was to evaluate the influence of CVO on the stress in the necrotic femoral head using an FEA. We also evaluated the effect of the LCE angle on the stress distribution in the segment of necrotic bone for better decision-making in relation to CVO.

Methods

This biomechanical study was approved by our ethical review committee (NO. 0-0672). A healthy volunteer (sex: male, age: 27 years, height: 164 cm, body weight: 66kg) without any musculoskeletal disease or history of hip joint operations was recruited. Quantitative computed tomography (QCT) was scanned in combination with a calibration phantom (B-MAS200, Kyoto-kagaku, Kyoto, Japan). The resolution of each CT image was 512 by 512 pixels with a slice thickness of 1.0 mm, and the pixel size was 0.782 mm/pixel under 120 KV and 102.50 mA conditions.

CVO model with ONFH

The intact hip joint model was constructed from the QCT data by segmenting the bony structure of the ilium and femur using a medical image processing software program (Mimics 22, Materialise, Belgium) (Fig. 1a). The cartilage geometry was created by filling the clearance gap between the acetabulum and the femoral head. Then, the cartilage model was divided into acetabular and femoral head cartilage with a spherical surface [21].

Fig. 1
figure 1

Construction of the CVO model with femoral head necrosis. a The hip joint model extracted from CT data. b The three locations with a cone angle of 60°: type A, type B, and type C1. c Intertrochanteric cured varus osteotomy with a rotation angle of 15°, 20°, 25°, and 30°

The necrotic bone model was created based on the Japanese Investigation Committee (JIC) classification [22]. The size of the necrotic lesion was set as a conoid shape (cone angle: 60°) to represent an early stage, low risk, and pre-collapse of the lesion [18]. The typical three locations of the necrotic bone (types A, B, C1) were simulated. (Fig. 1b). In brief, the weight-bearing portion of the acetabulum was trisected; the vertical line through lateral node of each aliquot reached the point over the femoral head that was defined as the boundary of the lateral locations of types A, B, C1 in the central coronal section of the femoral head of CT images. The locations matched with the changes of the weight-bearing portion, which were decided by different LCE angles in this study.

The CVO models were created using a CAD system (SolidWorks 2016, SolidWorks Corp., USA) (Fig. 1c). The osteotomy plane was set from the frontal view, and the femur model was cut along the plane by Boolean operations. The osteomized proximal part was rotated and fixed to the distal femur with four different angles of 15°, 20°, 25°, and 30° without implants.

3D Measurement of geometric parameters

Femoral head coverage area was measured in each hip joint model [23]. The spline curve was fit to the rim of the acetabulum (Fig. 2a, black line, lateral to medial) and projected on the surface of the femoral head to create an area of femoral head coverage (Fig. 2b). The weight-bearing area was defined as the lunate surface in this study. The border curve between the lunate surface and the acetabular fossa (Fig. 2a, red line) was projected on the surface of the femoral head to create the weight-bearing area (Fig. 2c). The intact area was defined as the area with the necrotic region removed from the weight-bearing area (Fig. 2d). The intact ratio was calculated as the intact area/the weight-bearing area. The results of the measurement are summarized in (Table 1).

Fig. 2
figure 2

Three-dimensional measurement of the femoral head coverage, weight-bearing area, and intact area. a The acetabular rim line (in black) for determining the acetabular coverage and acetabular arch line (in red) for determining the weight-bearing area. b The acetabular coverage projected on the femoral head in yellow. c The weight-bearing area projected on the femoral head in yellow. d1–d3 The intact area in yellow of types A, B, and C1 was determined by removing the necrotic area in orange (top view)

Table 1 3D-Measurement of geometrical parameters

Material properties assignment and boundary configurations

The hip joint models were meshed 928,127 four-node tetrahedral elements with an element size of approximately 1 mm for bone and cartilage, which are assumed to be isotropic and elastic materials. For the femur model, heterogeneous Young’s modulus was assigned to each element based on QCT data (Fig. 3a) [24]. Briefly, the parameters used for converting Hounsfield Units (HU) to radiographic CT density (ρQCT(g/cm3) (Eq. (1)) were calculated from the B-MAS200 phantom [25] and from ρQCT to Ash density (ρash(g/cm3) (Eq. (2)) [26], then the apparent density that was calculated from the Ash density with a ratio of 0.6 [27] was converted to the elastic modulus (Eq. (3)) [28].

$$ {\rho}_{QCT}\left(\mathrm{g}/{\mathrm{cm}}^3\right)=0.9863\mathrm{HU}\hbox{-} 2.0804 $$
(1)
$$ {\rho}_{ash}\left(\mathrm{g}/{\mathrm{cm}}^3\right)=0.877\times {\rho}_{QCT}+0.0789 $$
(2)
$$ E=6850{\rho}_{app}^{1.49} $$
(3)
Fig. 3
figure 3

The finite element model of CVO. a The femur was mapped with isotropic heterogeneous material properties based on grayscale values from the QCT data. b Lateral center-edge angle = 33° to represent the normal weight-bearing portion and lateral center-edge angle = 18° to represent the insufficient weight-bearing portion. c Ground reaction force was imposed at the distal end of the femur (red arrow); several muscles (blue dotted lines) and the ligaments (purple lines) around the hip joint were attached; the interfaces of the CVO were bounded without fixation

The necrotic lesion, cartilage and ilium were assigned with homogeneous material properties based on the literature (Table 2) [18, 29, 30].

Table 2 Material properties for FE model

To simulate the standing loading, a ground reaction force of 700 N was applied to a rigid plate fixed to the distal end of the femur. Seven muscles around the hip joint were modeled using connector elements based on a previous study (Table 3) [29]. The forces can minimize the bending moment in every cross-section of the femur [31]. The hip capsular ligaments were modeled as a 1D springs element. The stiffness of the superior iliofemoral, inferior iliofemoral, ischiofemoral, and pubofemoral ligaments were set at 96, 102, 40, and 36 N/mm, respectively. The width of ligament was simulated by several springs; the numbers of springs was adopted from a previous study [32]. The stiffness of each spring was calculated as the stiffness of ligament divided by the number of springs. The position of origin and insertion of springs was determined at the surface nodes of bone and mimicked the anatomical position of the ligament as much as possible [33]. The pubic symphysis and sacroiliac joint were fully fixed to prevent translation and rotation. The osteotomy interface was set as a bond, and the contact area of cartilage was modeled as frictionless.

Table 3 Loading for finite element model

We also modified the acetabulum shape to investigate the influence of the weight-bearing portion on the stress state in the necrotic lesion after CVO (Fig. 3b). The model with an LCE angle of 33° represented the normal weight-bearing portion, and the model with an LCE angle of 18° represented the insufficient portion [34, 35]. Finally, a total of 30 different finite element models were created simulating three locations of necrosis combined with four CVO angles in two different LCE angle conditions (Fig. 3c). The FE analysis was performed using a general-purpose FEA software program (ABAQUS 2019, Dassault Systems, Providence, RI).

Sensitivity studies and mesh convergence tests

Changes in the assumed material properties were performed to investigate how such properties affect the stress and strain predictions for sensitivity studies in physiological conditions using the intact hip joint model. We analyzed the variations in the FE model prediction, when the elastic modulus of the pelvis bone and articular cartilage were changed to ±10 % of the initial values [30]. The results of the sensitivity studies are presented in Table 4. In our selected models, the variation of the sensitivity studies was confirmed to be <5% by evaluating the maximum von Mises stress of the femur. This information demonstrated that the alteration in material parameters did not result in significant changes. According to the modulus and Poisson’s ratios, the FE-predicted stress was relatively insensitive to change [36]. This indicates that the material coefficient used in the study was stable and effective.

Table 4 Sensitivity studies of FE models

Mesh convergence tests were performed to confirm the discretization for the analysis. A mesh was generally considered to be sufficiently refined when an increase in mesh resolution yielded less than a 5% change in the result [37]. The number of tetrahedral elements in the intact femur model was increased to analyze the stress variation by evaluating the maximum von Mises stress of the femur (Table 5). Based on the convergence test, the elements of 928,127 with an average element size of 1 mm (percentage variation within 5%) was chosen as the final FE model to minimize the discretization error for the verification of the calculation [38].

Table 5 Convergence tests of FE models

Results

The effectiveness of CVO on stress reduction in necrotic bone

The maximum and mean von Mises stress in the necrotic bone were reduced with an increase in the CVO angle (Fig. 4). In particular, in the insufficient weight-bearing models (LCE angle of 18°), the maximum stress was reduced from 15.67 to 6.13 MPa for type C1 and from 13.73 to 3.76 MPa for type B (Fig. 4a). Correspondingly, the stress concentration at the interface of healthy and necrotic bone was weakened, accompanied by an increase in the CVO angle (Figs. 5 and 6). For type A, with both LCE angles, slight stress reduction was observed after CVO in comparison to types B and C1 (Fig. 4b, d).

Fig. 4
figure 4

Von Mises stress in different models. a Maximum von Mises stress of the models. b The maximum von Mises stress reduction after CVO. c The mean von Mises stress of the models. d The mean von Mises stress reduction after CVO

Fig. 5
figure 5

Top view of the von Mises stress distribution of types A, B, and C1 in two different weight-bearing portions with several CVO angles

Fig. 6
figure 6

Coronal central section view of the von Mises stress distribution of types A, B, C1 in two different weight-bearing portions with several CVO angles

Influence of the weight-bearing portion

The LCE angles affected the maximum and mean stress in the necrotic bone, with the exception of type A. The stress with an LCE angle of 18° were higher in comparison to the stress with an LCE angle of 33°, for types B (maximum stress: 13.73 vs 6.83 MPa) and C1 (maximum stress: 15.67 vs 8.14 MPa) in CVO 0° models (Fig. 4a). With an LCE angle of 18°, extensive stress distribution was observed in the necrotic region in comparison to an LCE angle of 33° (Figs. 5 and 6).

Relationship between geometric measurements and stress of necrosis

The correlation of these geometric measurements with the mean von Mises stress of the necrotic region was compared (Fig. 7). From the results of this linear regression analysis, the geometric measurements and stress had a negative relationship with good fitness, as determined by the R2 value. The stress of necrosis decreased by increasing the intact ratio, intact area, weight-bearing area, and femoral head coverage. However, the intact ratio had more influence on the stress of the necrotic bone in comparison to other geometric measurements with a large slope value.

Fig. 7
figure 7

The correlation of femoral head coverage, weight-bearing area, intact area, and intact ratio with the mean von Mises stress of necrosis of types A, B, and C1

Discussion

Proximal femoral osteotomy for the treatment of ONFH is performed with the aim of eliminating the mechanical factors suffered by necrotic bone. This is why the mechanical simulation of such procedures is necessary for assisting preoperative decision-making [9]. In this study, the biomechanical effectiveness of CVO was addressed by evaluating the stress of necrotic bone after CVO (Fig. 4). All types of necrotic lesions would mechanically benefit from CVO with different levels of stress reduction. As expected, in the case of types B and C1, the higher degree of CVO performed, the greater the degree of stress reduction that was achieved, because the necrotic bone is removed from the weight-bearing portion to a medial site by varus osteotomy. For example, the maximum stress of type C1 from 8.14 MPa in 0° situation was reduced to 4.53 MPa after a CVO of 30° with a maximum stress reduction of 44.3% (Fig. 4a, b). These results suggested that CVO was an effective way to reduce the stress level in necrotic bone.

Our simulation demonstrated that the appropriate CVO angle depended on the type of necrotic lesion. For type A, there was little change in the stress after CVO, and the stress level was relatively lower in 0° models. This result can support the clinical repots [7, 14, 16, 39,40,41,42,43,44]. Therefore, it was not effective to perform the CVO with type A. For type B, satisfactory results could be achieved from 20° (0.61 MPa), and for type C1, an angle of at least 30° (0.77 MPa) was considered safe, if the stress level of type A (0.61 MPa) was regarded as a standard, based on the results predicted in this study.

Unfavorable biomechanical situation of the hip such as acetabular dysplasia, which can increase contact stress as well as has been identified as a risk factor for the collapse of necrotic bone, leads to highly nonuniform contact stress distribution in comparison to the normal hip joint [16]. In this study, the stress in the insufficient weight-bearing portion models was higher in comparison to normal weight-bearing models. In the case of type C1, the mean stress value was 2.85 MPa of LCE angle 18° vs. 1.57 MPa of LCE angle 33° in the no osteotomy model, and it was higher (0.90 MPa vs 0.77 MPa) in the LCE angle 18° model even after CVO 30°. Type A had received the least influence from the weight-bearing portion (Fig. 4c). From these results, the influence of the weight-bearing portion on stress distribution of lesion was location-dependent, the more lateral site located, the more stress suffered. Thus, further procedures, such as periacetabular osteotomy, should be considered, when CVO is performed for a hip with an unfavorable biomechanical situation.

Stress concentrations were observed at not only the lateral interface region but also anterior region for types B and C1 (Figs. 5 and 6), which indicated that the anterior side was also at risk of collapse. The finding was supported by a clinical report on CVO [39]. This report suggested that the anterior necrotic angle was associated with the progressive collapse of the anterior lesion.

The intact ratio is used as the most important indicator for evaluating the risk of collapse after CVO in the clinical setting [14]. It is calculated from the transposed intact articular surface of the femoral head and the weight-bearing area of the acetabulum on anteroposterior radiographs [39]. However, the relationship between the intact ratio measured in 3D and the stress of necrotic bone after CVO was firstly discussed in this study. We found that a strong negative correction between the intact ratio and the stress level of necrotic bone after CVO had more influence on stress distribution with a larger slope value in the linear regression model, in comparison to the other parameters of femoral head coverage, weight-bearing area, and intact area (Fig. 7). Thus, it was confirmed that the intact ratio could be a bridge between the clinical results and the stress reduction of necrosis after CVO.

The present study was associated with some limitations. (1) Only one healthy subject was evaluated; thus, the deviation derived from patient-specific differences could not be considered [45, 46]. The study should have investigated additional types of necrosis from patient data. (2) Finite element models of physiological loading conditions should generally be validated from electromyography data [47,48,49]. For the present study, validation of the predicted results depended on similarity to the physical phenomenon of the real structure of femoral head with an intact model subjected to loading. (3) The interfaces of the CVO were assumed to be fully bonded, which more likely simulates the healed stage after CVO without taking the influence of the fixation device into account when simulating the initial postoperative stage. (4) The mechanical property of human articular cartilage displayed a viscoelastic and nonlinear mechanical response to loading [46, 50]; however, an isotropic homogeneous material property was assigned in the study [29, 35]. (5) This simulation only addressed mechanical factors that would affect the stress distribution of necrotic bone; however, biological reaction that may be different after each CVO also needs to be considered.

Conclusions

Understanding the complicated interdependence of the size and location of necrotic lesions and the configuration of the hip and pelvis is important in making decisions regarding optimal treatment. This finite element study suggested all types of lesions would mechanically benefit from CVO. However, for type B, an angle of >20° was recommended, and for type C1, an angle of at least 30° was safe, decreasing to the stress level of type A. An insufficient weight-bearing portion at acetabulum such as dysplasia can lead to the stress concentration at the necrotic bone rather than the normal hip joint. The intact ratio could be a bridge between the clinical outcomes and stress level in CVO.

Availability of data and materials

The data and materials are available from the corresponding author.

Abbreviations

CVO:

Curved varus osteotomy

ONFH:

Osteonecrosis of the femoral head

LCE:

Lateral center-edge angle

FEA:

Finite element analysis

QCT:

Quantitative computed tomography

JIC:

Japanese Investigation Committee

HU:

Hounsfield Unit

Eq:

Equation

CAD:

Computer-assisted design

vs:

Versus

3D:

Three-dimensional

References

  1. Moya-Angeler J, Gianakos AL, Villa JC, Ni A, Lane JM. Current concepts on osteonecrosis of the femoral head. World J Orthop. 2015;6(8):590–601. https://doi.org/10.5312/wjo.v6.i8.590.

    Article  PubMed  PubMed Central  Google Scholar 

  2. Mont MA, Carbone JJ, Fairbank AC. Core decompression versus nonoperative management for osteonecrosis of the hip. Clin Orthop Relat Res. 1996;324:169–78. https://doi.org/10.1097/00003086-199603000-00020.

    Article  Google Scholar 

  3. Veillette CJH, Mehdian H, Schemitsch EH, Mckee MD. Survivorship analysis and radiographic outcome following tantalum rod insertion for osteonecrosis of the femoral head. J Bone Joint Surg Am. 2006;88:48–55.

    PubMed  Google Scholar 

  4. Keizer SB, Kock NB, Dijkstra PDS, Taminiau AHM, Nelissen RGHH. Treatment of avascular necrosis of the hip by a non-vascularised cortical graft. J Bone Joint Surg Br. 2006;88(4):460–6.

    Article  CAS  Google Scholar 

  5. Chen WP, Tai CL, Shih CH, Hsieh PH, Leou MC, Lee MS. selection of fixation devices in proximal femur rotational osteotomy: clinical complications and finite element analysis. Clin Biomech. 2004;19(3):255–62. https://doi.org/10.1016/j.clinbiomech.2003.12.003.

    Article  Google Scholar 

  6. Lee YK, Park CH, Ha YC, Kim DY, Lyu SH, Koo KH. Comparison of surgical parameters and results between curved varus osteotomy and rotational osteotomy for osteonecrosis of the femoral head. Clin Orthop Surg. 2017;9(2):160–8. https://doi.org/10.4055/cios.2017.9.2.160.

    Article  PubMed  PubMed Central  Google Scholar 

  7. Sakano S, Hasegawa Y, Torii Y, Kawasaki M, Ishiguro N. Curved intertrochanteric varus osteotomy for osteonecrosis of the femoral head. J Bone Joint Surg Br. 2004;86(3):359–65.

    Article  CAS  Google Scholar 

  8. Nozawa M, Enomoto F, Shitoto K, Matsuda K, Maezawa K, Kurosawa H. Rotational acetabular osteotomy for osteonecrosis with collapse of the femoral head in young patients. J Bone Joint Surg. 2005;87(3):514–20. https://doi.org/10.2106/00004623-200503000-00006.

    Article  PubMed  Google Scholar 

  9. Baker KJ, Brown TD, Brand RA. A finite-element analysis of the effects of intertrochanteric osteotomy on stresses in femoral head osteonecrosis. Clin Orthop Relat Res. 1989;249:183–98.

    Google Scholar 

  10. Quaranta M, Miranda L, Oliva F, Aletto C, Maffulli N. Osteotomies for avascular necrosis of the femoral head. Br Med Bull. 2021;137(1):98–111. https://doi.org/10.1093/bmb/ldaa044.

    Article  PubMed  Google Scholar 

  11. Migliorini F, Maffulli N, Eschweiler J, Tingart M, Baroncini A. Core decompression isolated or combined with bone marrow-derived cell therapies for femoral head osteonecrosis. Expert Opin Biol Ther. 2021;21(3):423–30. https://doi.org/10.1080/14712598.2021.1862790.

    Article  CAS  PubMed  Google Scholar 

  12. Sadile F, Bernasconi A, Russo S, Maffulli N. Core decompression versus other joint preserving treatments for osteonecrosis of the femoral head: a meta-analysis. Br Med Bull. 2016;118(1):33–49. https://doi.org/10.1093/bmb/ldw010.

    Article  PubMed  PubMed Central  Google Scholar 

  13. Nishio A, Sugioka Y. A new technique of the varus osteotomy at the upper end of the femur. Orthop Trauma. 1971;20(3):381–6. https://doi.org/10.5035/nishiseisai.20.381.

    Article  Google Scholar 

  14. Zhao G, Yamamoto T, Ikemura S, Motomura G, Mawatari T, Nakashima Y, et al. Radiological outcome analysis of transtrochanteric curved varus osteotomy for osteonecrosis of the femoral head at a mean follow-up of 12.4 years. J Bone Joint Surg Br. 2010;92(6):781–6.

    Article  CAS  Google Scholar 

  15. Roush TF, Olson SA, Pietrobon R, Braga L, Urbaniak JR. Influence of acetabular coverage on hip survival after free vascularized fibular grafting for femoral head osteonecrosis. J Bone Joint Surg Am. 2006;88(10):2152–8. https://doi.org/10.2106/JBJS.E.00469.

    Article  PubMed  Google Scholar 

  16. Okura T, Hasegawa Y, Morita D, Osawa Y, Ishiguro N. What factors predict the failure of curved intertrochanteric varus osteotomy for the osteonecrosis of the femoral head? Arch Orthop Trauma Surg. 2016;136(12):1647–55. https://doi.org/10.1007/s00402-016-2563-0.

    Article  PubMed  Google Scholar 

  17. Brown TD, Mutschler TA, Ferguson AB. A non-linear finite element analysis of some early collapse processes in femoral head osteonecrosis. J Biomech. 1982;15(9):705–15. https://doi.org/10.1016/0021-9290(82)90024-0.

    Article  CAS  PubMed  Google Scholar 

  18. Lee MS, Tai CL, Senan V, Shih CH, Lo SW, Chen WP. The effect of necrotic lesion size and rotational degree on the stress reduction in transtrochanteric rotational osteotomy for femoral head osteonecrosis - a three-dimensional finite-element simulation. Clin Biomech. 2006;21(9):969–76. https://doi.org/10.1016/j.clinbiomech.2006.05.005.

    Article  Google Scholar 

  19. Li T, Huang Z, Li Y, Xue Z, Sun J, Gao H, et al. Prediction of collapse using patient-specific finite element analysis of osteonecrosis of the femoral head. Orthop Surg. 2019;11(5):794–800. https://doi.org/10.1111/os.12520.

    Article  PubMed  PubMed Central  Google Scholar 

  20. Yang JW, Koo KH, Lee MC, Yang P. D.Noh M, Kim SY, Kim K, Ha YC, Joun MS. Mechanics of femoral head osteonecrosis using three-dimensional finite element method. Arch Orthop Trauma Surg. 2002;122(2):88–92. https://doi.org/10.1007/s004020100324.

    Article  PubMed  Google Scholar 

  21. Henak CR, Carruth ED, Anderson AE, Harris MD, Ellis BJ, Peters CL, et al. Finite element predictions of cartilage contact mechanics in hips with retroverted acetabula. Osteoarthr Cartil. 2013;21(10):1522–9. https://doi.org/10.1016/j.joca.2013.06.008.

    Article  CAS  Google Scholar 

  22. Sugano N, Atsumi T, Ohzono K, Kubo T, Hotokebuchi T, Takaoka K. The 2001 revised criteria for diagnosis, classification, and staging of idiopathic osteonecrosis of the femoral head. J Orthop Sci. 2002;7(5):601–5. https://doi.org/10.1007/s007760200108.

    Article  PubMed  Google Scholar 

  23. Hansen BJ, Harris MD, Anderson LA, Peters CL, Weiss JA, Anderson AE. Correlation between radiographic measures of acetabular morphology with 3D femoral head coverage in patients with acetabular retroversion. Acta Orthop. 2012;83(3):233–9. https://doi.org/10.3109/17453674.2012.684138.

    Article  PubMed  PubMed Central  Google Scholar 

  24. Tanaka H, Yamako G, Kurishima H, Yamashita S, Mori Y, Chiba D, et al. Biomechanical analysis of supra-acetabular insufficiency fracture using finite element analysis. J Orthop Sci. 2018;23(5):825–33. https://doi.org/10.1016/j.jos.2018.04.005.

    Article  PubMed  Google Scholar 

  25. Knowles NK, Reeves JM, Ferreira LM. Quantitative Computed Tomography (QCT) derived Bone Mineral Density (BMD) in finite element studies: a review of the literature. J Exp Orthop. 2016;3(1):36. https://doi.org/10.1186/s40634-016-0072-2.

  26. Schileo E, Dall'ara E, Taddei F, Malandrino A, Schotkamp T, Baleani M, et al. An accurate estimation of bone density improves the accuracy of subject-specific finite element models. J Biomech. 2008;41(11):2483–91. https://doi.org/10.1016/j.jbiomech.2008.05.017.

    Article  PubMed  Google Scholar 

  27. Ali AA, Cristofolini L, Schileo E, Hu H, Taddei F, Kim RH, et al. Specimen-specific modeling of hip fracture pattern and repair. J Biomech. 2014;47(2):536–43. https://doi.org/10.1016/j.jbiomech.2013.10.033.

    Article  PubMed  Google Scholar 

  28. Morgan EF, Bayraktar HH, Keaveny TM. Trabecular bone modulus-density relationships depend on anatomic site. J Biomech. 2003;36(7):897–904. https://doi.org/10.1016/S0021-9290(03)00071-X.

    Article  PubMed  Google Scholar 

  29. Zhou GQ, Pang ZH, Chen QQ, He W, Chen ZQ, Chen LL, et al. Reconstruction of the biomechanical transfer path of femoral head necrosis: a subject-specific finite element investigation. Comput Biol Med. 2014;52:96–101. https://doi.org/10.1016/j.compbiomed.2014.04.002.

    Article  PubMed  Google Scholar 

  30. Ravera EP, Crespo MJ, Guarnieri FA, Braidot AA. Stress in human pelvis throughout the gait cycle: development, evaluation and sensitivity studies of a finite element model. IFMBE Proceeding. 2015;49:246–9. https://doi.org/10.1007/978-3-319-13117-7_64.

    Article  Google Scholar 

  31. Sverdlova NS, Witzel U. Principles of determination and verification of muscle forces in the human musculoskeletal system: muscle forces to minimise bending stress. J Biomech. 2010;43(3):387–96. https://doi.org/10.1016/j.jbiomech.2009.09.049.

    Article  PubMed  Google Scholar 

  32. Zou Z, Chávez-Arreola A, Mandal P, Board TN, Alonso-Rasgado T. Optimization of the position of the acetabulum in a ganz periacetabular osteotomy by finite element analysis. J Orthop Res. 2013;31(3):472–9. https://doi.org/10.1002/jor.22245.

    Article  PubMed  Google Scholar 

  33. Hewitt J, Guilak F, Glisson R, Vail TP. Regional material properties of the human hip joint capsule ligaments. J Orthop Res. 2001;19(3):359–64. https://doi.org/10.1016/S0736-0266(00)00035-8.

    Article  CAS  PubMed  Google Scholar 

  34. Wylie JD, Kapron AL, Peters CL, Aoki SK, Maak TG. Relationship between the lateral center-edge angle and 3-dimensional acetabular coverage. Orthop J Sport Med. 2017;5(4):2325967117700589.

    Article  Google Scholar 

  35. Zhao X, Chosa E, Totoribe K, Deng G. Effect of periacetabular osteotomy for acetabular dysplasia clarified by three-dimensional finite element analysis. J Orthop Sci. 2010;15(5):632–40. https://doi.org/10.1007/s00776-010-1511-z.

    Article  PubMed  Google Scholar 

  36. Anderson AE, Peters CL, Tuttle BD, Weiss JA. Subject-specific finite element model of the pelvis: development, validation and sensitivity studies. J Biomech Eng. 2005;127(3):364–73. https://doi.org/10.1115/1.1894148.

    Article  PubMed  Google Scholar 

  37. Anderson AE, Ellis BJ, Weiss JA. Verification, validation and sensitivity studies in computational biomechanics. Comput Methods Biomech Biomed Engin. 2007;10(3):171–84. https://doi.org/10.1080/10255840601160484.

    Article  PubMed  PubMed Central  Google Scholar 

  38. Oreskes N, Shrader-Frechette K, Belitz K. Verification, validation, and confirmation of numerical models in the earth sciences. Science. 1994;263(5147):641–6. https://doi.org/10.1126/science.263.5147.641.

    Article  CAS  PubMed  Google Scholar 

  39. Kubo Y, Motomura G, Ikemura S, Hatanaka H, Utsunomiya T, Hamai S, et al. Effects of anterior boundary of the necrotic lesion on the progressive collapse after varus osteotomy for osteonecrosis of the femoral head. J Orthop Sci. 2020;25(10):145–51. https://doi.org/10.1016/j.jos.2019.02.014.

    Article  PubMed  Google Scholar 

  40. Shi L, Fan Y, Zhang C, Shen J. Value of virtual reality technology in preoperative planning of transtrochanteric curved varus osteotomy for avascular necrosis of femoral head in adults. Zhongguo Xiu Fu Chong Jian Wai Ke Za Zhi. 2019;33(8):923–8. https://doi.org/10.7507/1002-1892.201903083.

    Article  PubMed  Google Scholar 

  41. Osawa Y, Seki T, Okura T, Takegami Y, Ishiguro N, Hasegawa Y. Curved intertrochanteric varus osteotomy vs total hip arthroplasty for osteonecrosis of the femoral head in patients under 50 years old. J Arthroplasty. 2020;35(6):1600–5. https://doi.org/10.1016/j.arth.2020.01.026.

    Article  PubMed  Google Scholar 

  42. Takegami Y, Komatsu D, Seki T, Ishiguro N, Hasegawa Y. Total hip arthroplasty after failed curved intertrochanteric varus osteotomy for avascular necrosis of the femoral head. Nagoya J Med Sci. 2016;78(1):89–97.

    PubMed  PubMed Central  Google Scholar 

  43. Ikemura S, Yamamoto T, Jingushi S, Nakashima Y, Mawatari T, Iwamoto Y. Leg-length discrepancy after transtrochanteric curved varus osteotomy for osteonecrosis of the femoral head. J Bone Joint Surg Br. 2007;89(6):725–9.

    Article  CAS  Google Scholar 

  44. Hamanishi M, Yasunaga Y, Yamasaki T, Mori R, Shoji T, Ochi M. The clinical and radiographic results of intertrochanteric curved varus osteotomy for idiopathic osteonecrosis of the femoral head. Arch Orthop Trauma Surg. 2014;134(3):305–10. https://doi.org/10.1007/s00402-013-1919-y.

    Article  PubMed  Google Scholar 

  45. Shi J, Chen J, Wu J, Chen F, Huang G, Wang Z, et al. Evaluation of the 3D finite element method using a tantalum rod for osteonecrosis of the femoral head. Med Sci Monit. 2014;20:2556–64. https://doi.org/10.12659/MSM.890920.

    Article  PubMed  PubMed Central  Google Scholar 

  46. Bae JY, Kwak DS, Park KS, Jeon I. Finite element analysis of the multiple drilling technique for early osteonecrosis of the femoral head. Ann Biomed Eng. 2013;41(12):2528–37. https://doi.org/10.1007/s10439-013-0851-1.

    Article  PubMed  Google Scholar 

  47. Kaze AD, Maas S, Arnoux PJ, Wolf C, Pape D. A finite element model of the lower limb during stance phase of gait cycle including the muscle forces. BioMed Eng OnLine. 2017;16(1):138. https://doi.org/10.1186/s12938-017-0428-6.

    Article  Google Scholar 

  48. Seo JW, Kang DW, Kim JY, Yang ST, Kim DH, Choi JS, et al. Finite element analysis of the femur during stance phase of gait based on musculoskeletal model simulation. Biomed Mater Eng. 2014;24(6):2485–93. https://doi.org/10.3233/BME-141062.

    Article  PubMed  Google Scholar 

  49. Hicks JL, Uchida TK, Seth A, Rajagopal A, Delp SL. Is my model good enough? Best practices for verification and validation of musculoskeletal models and simulations of movement. J Biomech Eng. 2015;137(2):020905. https://doi.org/10.1115/1.4029304.

    Article  PubMed  Google Scholar 

  50. Robinson DL, Kersh ME, Walsh NC, Ackland DC, de Steiger RN, Pandy MG. Mechanical properties of normal and osteoarthritic human articular cartilage. J Mech Behav Biomed Mater. 2016;61:96–109. https://doi.org/10.1016/j.jmbbm.2016.01.015.

    Article  CAS  PubMed  Google Scholar 

Download references

Acknowledgements

I would like to thank Kenshin Matsue for using his valuable time to check the results. I would also like to thank the Kyoto-kagaku company for providing the calibration phantom (B-MAS200) for the subject-specific bone mineral density measurement.

Funding

This work was supported by JSPS KAKENHI Grant Number JP17K01362. The funder had no role in design of the study, the collection, analysis, and interpretation of the data, or in writing the manuscript.

Author information

Authors and Affiliations

Authors

Contributions

YW, GY, YN, and EC designed the study and collected the materials and data. YW, TO, and GY conducted the finite element simulation. YW and GY played a major role in writing the manuscript. HA, YN, and EC developed the research question and guided the study. The authors read and approved the final manuscript before submission.

Corresponding author

Correspondence to Go Yamako.

Ethics declarations

Ethics approval and consent to participate

The study was approved by the Ethical Review Committee of Miyazaki University (NO. 0-0672), and informed consent was obtained from subjects.

Consent for publication

Not applicable

Competing interests

The authors declare that they have no competing interest.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Wang, Y., Yamako, G., Okada, T. et al. Biomechanical effect of intertrochanteric curved varus osteotomy on stress reduction in femoral head osteonecrosis: a finite element analysis. J Orthop Surg Res 16, 465 (2021). https://doi.org/10.1186/s13018-021-02614-z

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s13018-021-02614-z