- Research article
- Open Access
Biomechanical effect of intertrochanteric curved varus osteotomy on stress reduction in femoral head osteonecrosis: a finite element analysis
Journal of Orthopaedic Surgery and Research volume 16, Article number: 465 (2021)
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.
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.
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.
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.
Osteonecrosis of the femoral head (ONFH) is a progressive pathological process that occurs due to multiple factors affecting the blood supply of femoral head . 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 .
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 , whereas it was prone to collapse if the necrotic lesion was located in the anterolateral weight-bearing portion or with insufficient femoral head coverage . 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 . 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.
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 .
The necrotic bone model was created based on the Japanese Investigation Committee (JIC) classification . 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 . 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 . 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).
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) . 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  and from ρQCT to Ash density (ρash(g/cm3) (Eq. (2)) , then the apparent density that was calculated from the Ash density with a ratio of 0.6  was converted to the elastic modulus (Eq. (3)) .
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) . The forces can minimize the bending moment in every cross-section of the femur . 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 . 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 . 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.
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 . 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 . This indicates that the material coefficient used in the study was stable and effective.
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 . 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 .
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).
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.
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 . 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 . 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 . 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 . It is calculated from the transposed intact articular surface of the femoral head and the weight-bearing area of the acetabulum on anteroposterior radiographs . 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.
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.
Curved varus osteotomy
Osteonecrosis of the femoral head
Lateral center-edge angle
Finite element analysis
Quantitative computed tomography
Japanese Investigation Committee
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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
The authors declare that they have no competing interest.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
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