Skip to main content

Dynamic hip kinematics during squatting before and after total hip arthroplasty

Abstract

Background

The difference in in vivo kinematics before and after total hip arthroplasty (THA) for the same subjects and the clearance between the liner and neck during squatting have been unclear. The purpose of the present study was to clarify (1) the changes in the in vivo kinematics between prosthetic hips and osteoarthritis hips of the same subjects and (2) the extent of the liner-to-neck clearance during squatting under weight-bearing conditions.

Methods

This study consisted of 10 patients who underwent unilateral THA for symptomatic osteoarthritis. Using a flat-panel X-ray detector, we obtained continuous radiographs during squatting. We analyzed the hip joint’s movements using three-dimensional-to-two-dimensional model-to-image registration techniques. We also quantified the minimum distance at maximum flexion and extension, and the minimum angle at maximum flexion between the liner and stem neck.

Results

The maximum hip flexion angles post-THA (80.7° [range, 69.4–98.6°]) changed significantly compared with the pre-THA values (71.7° [range, 55.2°–91.2°]). The pelvic tilt angle (posterior +, anterior−) at the maximum hip flexion post-THA (10.4° [range, − 6.7° to 26.9°]) was significantly smaller than that at pre-THA (16.6° [range, − 3° to 40.3°]). The minimum anterior and posterior liner-to-neck distances averaged 10.9 and 8.0 mm, respectively, which was a significant difference. The minimum liner-to-neck angle at maximum flexion averaged 34.7° (range, 20.7°–46.3°). No liner-to-neck contact occurred in any of the hips.

Conclusion

THA increased the range of hip joint motion and the pelvis tilted anteriorly more after than before THA, with sufficient liner-to-neck clearance during squatting. These data may be beneficial for advising patients after THA regarding postoperative activity restrictions in daily life.

Background

Total hip arthroplasty (THA) is the best surgical procedure for patients with end-stage osteoarthritis (OA) of the hip joint and is highly effective in relieving pain and improving function [1,2,3]. The clinical success of THA allows patients to resume their activities of daily living easily, and this includes squatting motions [4]. In Non-Western cultures, squatting is one of the fundamental activities of daily living and requires deep ranges of motion in flexion [5, 6]. Due to such requirements for deep flexibility of the prosthetic hip, there are concerns for impingement and dislocation after THA. Therefore, the knowledge of the in vivo kinematics of the prosthetic hip associated with squatting and the clearance between the liner and neck (liner-to-neck clearance) could be beneficial for advising for patients after THA regarding postoperative activity restrictions in daily life.

Accurate evaluations of kinematics under weight-bearing conditions have been achieved using three-dimensional (3D)-to-two-dimensional (2D) model-to-image registration techniques [7,8,9]. Recently, these techniques have been applied to the kinematic analyses of hips affected by OA and subsequent THA procedures [10, 11]. We previously reported that, with respect to squatting, patients with OA were unable to deeply flex their femurs due to limited range of motion (ROM) of the hip joints. Additionally, patients with OA tilted their pelvis more posteriorly to maintain a deeply flexed posture than did healthy subjects [10]. In prosthetic hips, Koyanagi et al. revealed that the mean maximum hip flexion ROM was 86.2 °, which is smaller than the mean maximum hip flexion angle of 95.4 ° in normal native hips [5, 12]. However, to the best of our knowledge, no previous report has demonstrated the changes in in vivo 3D kinematics during squatting pre- and post-THA for the same subjects using the accurate 3D-to-2D model-to-image registration techniques.

The purpose of the present study was to clarify (1) the changes in the in vivo kinematics between prosthetic hips and OA hips of the same subjects and (2) the extent of the liner-to-neck clearance during squatting under weight-bearing conditions.

Methods

Patients

The protocol in the current study was approved by our institutional review board. All patients provided informed consent to participate in this study. Between August 2012 and October 2015, 177 patients with 200 hips underwent primary cementless THA by two senior surgeons (Y.N. and S.H.). Among the original 177 patients, 10 of them (10 hips) met all the inclusion criteria, which were (1) unilateral THA for symptomatic OA (2) no previous surgery of the ipsilateral hip, (3) no previous surgery of the spine or other joints, (4) obtaining the informed consent prior to THA, and (5) ability to squat safely without assistance before and after THA. Demographic data of the patients are shown in Table 1. There were 6 women and 4 men with a mean age at the time of THA of 65 ± 8 years (range, 55–84 years). The mean height was 156 ± 7 cm (range, 148–168 cm), and the mean body mass index (BMI) was 23 ± 4 kg/m2 (range, 18–31 kg/m2). According to the Kellgren-Lawrence scale, 2 hips were classified as grade III and 8 hips were classified as grade IV. [13] The pre- and post-THA mean Harris hip scores were 48 ± 7 (range, 40–57) and 95 ± 3 (range, 91–99), respectively. The mean post-THA follow-up was 43 ± 11 months (range, 23–59 months). No patient had a history of any complication after THA.

Table 1 Demographic data

Surgical procedures and implants

All operations were performed through a posterolateral approach with repair of the posterior soft tissue [14]. The combined anteversion technique was adopted to cope with the wide range of femoral anteversion of hip dysplasia at our institution. The cup was placed according to stem anteversion so that combined anteversion ranged from 40° to 60° [15]. A cementless hemispherical press-fit cup, straight metaphyseal fit stem, and high cross-linked ultra-high molecular weight polyethylene liner (AMS and PerFix HA; Aeonian; Kyocera, Kyoto, Japan) were used [16, 17]. All femoral heads were alumina ceramic, and the head sizes were 32 mm in 9 cases and 26 mm in 1 case. There are two types of rims in the Aeonian liner, one is a 15 ° elevated rim (elevated rim liner) to prevent dislocation of the femoral head, and the other is not elevated (flat liner). Four hips used the elevated rim liners, and 6 hips used the flat liner, respectively. When using the elevated rim liner, we recorded the location on the cup where the top of the rim was placed.

The orientations of the acetabular cup and stem were measured using the postoperative computed tomography (CT) data. The cup inclination was measured as the angle of abduction using the inter-tear-drop line as the baseline (radiographic inclination). The cup anteversion was measured as the angle of anteversion in the sagittal plane (operative anteversion). Femoral anteversion was measured as the angle of anteversion between the prosthetic femoral neck and transepicondylar axis [18].

Kinematic analysis

We essentially followed the method in accordance with previous reports, and we partially referenced the kinematics data regarding pre-THA [10, 19]. The 3D positions and orientations of the pelvis, acetabular cup, femur, and femoral stem during squatting were determined using 3D-to-2D model-to-image registration techniques. Continuous radiographic images were used to survey squatting movements using a flat-panel X-ray detector (Ultimax-I, Toshiba, Tochigi, Japan) with the following parameters: image area of 420 mm × 420 mm, resolution of 0.274 mm × 0.274 mm/pixel, and frame rate of 3.5 frames/s (Fig. 1). Each subject routinely underwent CT (Aquilion, Toshiba, Tochigi, Japan) with a 512 × 512 image matrix, a 0.35 × 0.35 pixel dim, and 1-mm-slice thicknesses from the superior edge of the pelvis to just below the knee joint line. Anatomical coordinate systems of the pelvis and femur were embedded in each bone model derived from CT data according to our previous study [10, 19]. Computer simulation was performed to generate virtual, digitally reconstructed radiographs in which the light source and the projected plane parameters were set to be identical to the actual radiographic imaging conditions. Each model silhouette was matched with the actual silhouette by translating and rotating the 3D model to minimize the number of unmatched pixels between the silhouettes. The orientation of the femur relative to the pelvis: hip movements, was determined using the Cardan/Euler angle system in x-y-z order (flexion/extension, adduction/abduction, internal/external rotation). The maximum errors associated with tracking the position of the femur/stem relative to the pelvis/acetabular cup were 0.36/0.43 mm, 0.37/0.48 mm, and 0.48 °/0.52 °, respectively, for in-plane translation, out-of-plane translation, and rotation, respectively [10, 19].

Fig. 1
figure 1

The hip motions during squatting were captured as continuous X-ray images using a flat panel X-ray detector. Patients before and after total hip arthroplasty stood from a squat position with their heel down

Regarding the liner-to-neck clearance, we quantified the minimum distance at maximum flexion and extension and the angle at maximum flexion between the liner and stem neck using a computer-aided design software program ([CATIA V5]; Dassault Systemes, Vélizy-Villacoublay, France) (Fig. 2) [19, 20].

Fig. 2
figure 2

The minimum distance at maximum flexion (a) at maximum extension (b), and the angle at maximum flexion (c) between the liner and stem neck

Statistical analysis

All data are expressed as mean ± standard deviation (SD) and were tested for normality using the Shapiro-Wilk test. To compare hip kinematics before and after THA in the same patient, normally distributed variables were evaluated using the paired t test. Non-normally distributed variables were evaluated using the independent Wilcoxon signed-rank test. Statistical significance was defined as a P value < .05. The statistical analyses were performed using JMP Software (Version 11; SAS Institute, Cary, NC, USA).

Results

Orientation of the components

The mean cup inclination and mean cup anteversion were 39.1 ± 5.4° (range, 30.1°–48.0°) and 18.5 ± 12.0° (range, 4.3°–40.0°), respectively. The mean stem anteversion was 31.7 ± 6.8° (range, 22.0°–42.1°), and the mean combined anteversion was 50.3 ± 9.3° (range, 36.9°–67.5°) (Table 2).

Table 2 Component data

Kinematics of the hip joint

The maximum hip flexion angles were determined during squatting at 10% of the squat ascent cycle pre-THA, and at 15% of the squat ascent cycle post-THA, respectively (Fig. 3). The kinematics data of each patient are shown in Table 3. The maximum femoral and hip flexion angles post-THA (94.9 ± 6.8° [range, 85.4°–104.6°] and 80.7 ± 10.0° [range, 69.4°–98.6°], respectively) changed significantly (P = .013 and P = .005, respectively) compared with the pre-THA values (86.2 ± 7.3° [range, 70.9°–94.3°] and 71.7 ± 11.9° [range, 55.2°–91.2°], respectively) (Figs. 3 and 4). In the squatting position, the pelvis tilted backward and gradually tilted forward with standing. The pelvic tilt angle (posterior +, anterior−) at the maximum hip flexion post-THA (10.4 ± 10.4° [range, − 6.7° to 26.9°] was significantly (P = .0463) smaller than that at pre-THA (16.6 ± 13.3° [range, − 3° to 40.3°]) (Fig. 5).

Fig. 3
figure 3

The hip flexion/extension angles during squatting in patients before total hip arthroplasty (THA) (gray lines) and patients after THA (black lines). Error bars show standard deviation

Table 3 Patient kinematic data
Fig. 4
figure 4

The femoral flexion/extension angles during squatting in patients before total hip arthroplasty (THA) (gray lines) and patients after THA (black lines). Error bars show standard deviation

Fig. 5
figure 5

Posterior/anterior pelvic tilt angles [posterior +, anterior−] during squatting in patients before total hip arthroplasty (THA) (gray lines) and patients after THA (black lines). Error bars show standard deviation

Liner-to-neck clearance

The liner-to-neck clearance data for the each patient are shown in Table 3. The minimum anterior and posterior liner-to-neck distances averaged 10.9 ± 2.2 mm (range, 7.3–14.8 mm) and 8.0 ± 2.4 mm (range, 3.6–11.2 mm), respectively, and the difference was significant (P = .0143). The minimum liner-to-neck angle at maximum flexion averaged 34.7 ± 7.3° (range, 20.7–46.3°). No liner-to-neck contact occurred in any of the hips during squatting.

Discussion

To our knowledge, no previous reports have demonstrated the changes in in vivo hip kinematics during squatting pre- and post-THA for the same subjects using accurate 3D-to-2D model-to-image registration techniques. The current study elucidated the dynamic hip kinematics during squatting before and after THA and quantified the liner-to-neck clearance. The maximum hip flexion significantly changed from 72° pre-THA to 81° post-THA, and the pelvic tilt angle at the maximum hip flexion post-THA was 10° , which was significantly smaller than the 17° measurement pre-THA. The anterior and posterior liner-to-neck distances averaged 11 and 8 mm, respectively, and the differences were significant. The mean minimum liner-to-neck angle was 35°, and there was no liner-to-neck contact during squatting in any of the hips.

A previous study found that patients with OA were not able to flex their femurs deeply due to limited ROM of the hip joints during squatting. They also tend to tilt their pelvises more posteriorly to maintain a deeply flexed posture than did healthy subjects [10]. The present study revealed that the maximum hip flexion was improved by approximately 9°. The pelvic tilt at the maximum hip flexion post-THA was significantly inclined and measured approximately 7° more anteriorly than pre-THA for the same subjects. We found that THA increased the range of the femoral and hip joint motions during squatting, and the compensation to maintain deeply flexed postures was reduced. Sagittal pelvic mobility is integral in flexing the torso to maintain balance and to allow the large hip flexion angles that are essential for deep squats [21].

It has been reported that approximately 95 to 102° of maximum hip flexion occurs during squatting in healthy subjects [5, 22]. Our study showed that the maximum hip flexion significantly changed from 72° pre-THA to 81° post-THA. Although THA provided increased ROM in patients with OA, prosthetic hips were not able to recover the kinematics to the level of healthy hips. Catelli et al. reported that joint kinematics of the pelvis and hip do not return to the level of healthy hips after THA during squatting and dual-mobility implant combined with the poorer functional scores [23]. A limited range of hip flexion persisted even after THA, which might have affected postoperative functional outcomes during activities of daily living, especially those requiring deeply flexed postures.

Specific postures might include potential risks of prosthetic impingement and could cause postoperative dislocation, and polyethylene wear after THA [24, 25]. However, we found only a few reports which quantified the liner-to-neck distances while squatting. In our study, the anterior and posterior liner-to-neck distances averaged 11 and 8 mm, respectively, and the differences were significant. Furthermore, the minimum liner-to-neck angle at maximum flexion was 21°. Koyanagi et al. reported that the minimum angle leading up to theoretical prosthetic impingement was more than 10° [12]. Based on the analysis, we concluded that the liner-to-neck clearance during squatting after THA was sufficient in these patients.

There are several limitations in the present study. First, the number of patients in this study was small. A larger sample size with a wider variation in component position might result in greater statistical reliability and reveal correlation between implant positioning and minimum liner-to-neck clearance; this topic needs further clarification. Second, kinematic processing of radiographic measurements imparts a risk of radiation exposure. However, we believe that this represents an important data-driven approach to provide feedback on the activities of daily living specific advice for each patient. Third, this study examined only a single component design. Although the design is similar to that of many others that are currently available, the results could differ. Fourth, the sequential motions during squatting were collected into twice, because even the large flat-panel X-ray detector that was used in this study provided a limited field of view (FOV). It might be necessary that the development of fluoroscopy is expected to achieve quite a large FOV. Finally, the patients included in this study were all Japanese with lower BMI compared to the Caucasian average. The patients with a BMI ≥ 30 kg/m2 could show different kinematic data with limited maximum hip flexion angle.

Conclusion

We quantified the change in hip kinematics before and after THA and the liner-to-neck clearance while squatting using the 3D-to-2D model-to-image registration techniques. THA increased the ranges of femoral and hip joint motion and the pelvis tilted anteriorly more after than before THA, with sufficient liner-to-neck clearance during squatting. These data may be beneficial for advising patients after THA regarding postoperative activity restrictions in daily life.

Abbreviations

2D:

Two-dimensional

3D:

Three-dimensional

BMI:

Body mass index

CT:

Computed tomography

DDH:

Developmental dysplasia of the hip

FOV:

Field of view

HHS:

Harris hip score

OA:

Osteoarthritis

ROM:

Range of motion

SD:

Standard deviation

THA:

Total hip arthroplasty

References

  1. Learmonth ID, Young C, Rorabeck C. The operation of the century: total hip replacement. Lancet. 2007;370(9597):1508–19.

    Article  PubMed  Google Scholar 

  2. Chang RW, Pellisier JM, Hazen GB. A cost-effectiveness analysis of total hip arthroplasty for osteoarthritis of the hip. JAMA. 1996;275(11):858–65.

    Article  PubMed  CAS  Google Scholar 

  3. Jacobsen S, Sonne-Holm S, Soballe K, Gebuhr P, Lund B. Hip dysplasia and osteoarthrosis: a survey of 4151 subjects from the osteoarthrosis substudy of the Copenhagen City Heart Study. Acta Orthop. 2005;76(2):149–58.

    Article  PubMed  Google Scholar 

  4. Rissanen P, Aro S, Slatis P, Sintonen H, Paavolainen P. Health and quality of life before and after hip or knee arthroplasty. J Arthroplast. 1995;10(2):169–75.

    Article  CAS  Google Scholar 

  5. Hemmerich A, Brown H, Smith S, Marthandam SS, Wyss UP. Hip, knee, and ankle kinematics of high range of motion activities of daily living. J Orthop Res. 2006;24(4):770–81.

    Article  PubMed  CAS  Google Scholar 

  6. Tang H, Du H, Tang Q, Yang D, Shao H, Zhou Y. Chinese patients’ satisfaction with total hip arthroplasty: what is important and dissatisfactory? J Arthroplast. 2014;29(12):2245–50.

    Article  Google Scholar 

  7. Komistek RD, Dennis DA, Mahfouz M. In vivo fluoroscopic analysis of the normal human knee. Clin Orthop Relat Res. 2003;410:69–81.

    Article  Google Scholar 

  8. Hamai S, Moro-oka TA, Miura H, Shimoto T, Higaki H, Fregly BJ, Iwamoto Y, Banks SA. Knee kinematics in medial osteoarthritis during in vivo weight-bearing activities. J Orthop Res. 2009;27(12):1555–61.

    Article  PubMed  Google Scholar 

  9. Ishimaru M, Shiraishi Y, Ikebe S, Higaki H, Hino K, Onishi Y, Miura H. Three-dimensional motion analysis of the patellar component in total knee arthroplasty by the image matching method using image correlations. J Orthop Res. 2014;32(5):619–26.

    Article  PubMed  Google Scholar 

  10. Hara D, Nakashima Y, Hamai S, Higaki H, Ikebe S, Shimoto T, Yoshimoto K, Iwamoto Y. Dynamic hip kinematics in patients with hip osteoarthritis during weight-bearing activities. Clin Biomech. 2016;32:150–6.

    Article  Google Scholar 

  11. Tsai T-Y, Li J-S, Wang S, Lin H, Malchau H, Li G, Rubash H, Kwon Y-M. A novel dual fluoroscopic imaging method for determination of THA kinematics: in-vitro and in-vivo study. J Biomech. 2013;46(7):1300–4.

    Article  PubMed  Google Scholar 

  12. Koyanagi J, Sakai T, Yamazaki T, Watanabe T, Akiyama K, Sugano N, Yoshikawa H, Sugamoto K. In vivo kinematic analysis of squatting after total hip arthroplasty. Clin Biomech. 2011;26(5):477–83.

    Article  Google Scholar 

  13. Kellgren JH, Lawrence JS. Radiological assessment of osteo-arthrosis. Ann Rheum Dis. 1957;16(4):494–502.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  14. Pellicci PM, Bostrom M, Poss R. Posterior approach to total hip replacement using enhanced posterior soft tissue repair. Clin Orthop Relat Res. 1998;355:224–8.

    Article  Google Scholar 

  15. Nakashima Y, Hirata M, Akiyama M, Itokawa T, Yamamoto T, Motomura G, Ohishi M, Hamai S, Iwamoto Y. Combined anteversion technique reduced the dislocation in cementless total hip arthroplasty. Int Orthop. 2014;38(1):27–32.

    Article  PubMed  Google Scholar 

  16. Sato T, Nakashima Y, Akiyama M, Yamamoto T, Mawatari T, Itokawa T, Ohishi M, Motomura G, Hirata M, Iwamoto Y. Wear resistant performance of highly cross-linked and annealed ultra-high molecular weight polyethylene against ceramic heads in total hip arthroplasty. J Orthop Res. 2012;30(12):2031–7.

    Article  PubMed  CAS  Google Scholar 

  17. Nakashima Y, Sato T, Yamamoto T, Motomura G, Ohishi M, Hamai S, Akiyama M, Hirata M, Hara D, Iwamoto Y. Results at a minimum of 10 years of follow-up for AMS and PerFix HA-coated cementless total hip arthroplasty: impact of cross-linked polyethylene on implant longevity. J Orthop Sci. 2013;18(6):962–8.

    Article  PubMed  CAS  Google Scholar 

  18. Yoshioka Y, Siu D, Cooke TD. The anatomy and functional axes of the femur. J Bone Joint Surg Am. 1987;69(6):873–80.

    Article  PubMed  CAS  Google Scholar 

  19. Hara D, Nakashima Y, Hamai S, Higaki H, Ikebe S, Shimoto T, Yoshimoto K, Iwamoto Y. Dynamic hip kinematics during the golf swing after total hip arthroplasty. Am J Sports Med. 2016;44(7):1801–9.

    Article  PubMed  Google Scholar 

  20. Yoshimoto K, Hamai S, Higaki H, Gondoh H, Nakashima Y. Visualization of a cam-type femoroacetabular impingement while squatting using image-matching techniques: a case report. Skelet Radiol. 2017;46(9):1277–82.

    Article  Google Scholar 

  21. Lamontagne M, Kennedy MJ, Beaule PE. The effect of cam FAI on hip and pelvic motion during maximum squat. Clin Orthop Relat Res. 2009;467(3):645–50.

    Article  PubMed  Google Scholar 

  22. Hara D, Nakashima Y, Hamai S, Higaki H, Ikebe S, Shimoto T, Hirata M, Kanazawa M, Kohno Y, Iwamoto Y. Kinematic analysis of healthy hips during weight-bearing activities by 3D-to-2D model-to-image registration technique. Biomed Res Int. 2014;2014:457573.

    Article  PubMed  PubMed Central  Google Scholar 

  23. Catelli DS, Kowalski E, Beaule PE, Lamontagne M. Does the dual-mobility hip prosthesis produce better joint kinematics during extreme hip flexion task? J Arthroplast. 2017;32(10):3206–12.

    Article  Google Scholar 

  24. Shon WY, Baldini T, Peterson MG, Wright TM, Salvati EA. Impingement in total hip arthroplasty a study of retrieved acetabular components. J Arthroplasty. 2005;20(4):427–35.

    Article  PubMed  Google Scholar 

  25. Marchetti E, Krantz N, Berton C, Bocquet D, Fouilleron N, Migaud H, Girard J. Component impingement in total hip arthroplasty: frequency and risk factors. A continuous retrieval analysis series of 416 cup. Orthop Traumatol Surg Res. 2011;97(2):127–33.

    Article  PubMed  CAS  Google Scholar 

Download references

Funding

This study was supported by JSPS KAKENHI Grant No. 15K10450 and 25870499, grant from the Japan Orthopaedics and Traumatology Foundation, Inc. (No. 263), and grant from the Nakatomi Foundation.

Availability of data and materials

The datasets supporting the conclusions of the present study are available from the corresponding author on reasonable request.

Author information

Authors and Affiliations

Authors

Contributions

KK, SH, DH, SI, HH, and YN contributed to the conception and design of the study. KK, SH, DH, KY, and KS performed acquisition of data. KK, SH, DH, SI, HH, HG, and YW conducted data analysis, and KK, SH, and YN contributed to data interpretation and preparation of the manuscript. All authors read and approved the final version of the manuscript.

Corresponding author

Correspondence to Satoshi Hamai.

Ethics declarations

Ethics approval and consent to participate

All procedures performed in studies involving human participants were in accordance with the ethical standards of our institution and the 1964 Helsinki Declaration and its later amendments or comparable ethical standards. This study was approved by the Institutional Review Board (IRB) of Kyushu University (IRB number 24–55). The subject was given informed consent to participate in this IRB approved study and was informed of the risk of radiation exposure required.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

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 distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. 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.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Komiyama, K., Hamai, S., Hara, D. et al. Dynamic hip kinematics during squatting before and after total hip arthroplasty. J Orthop Surg Res 13, 162 (2018). https://doi.org/10.1186/s13018-018-0873-3

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s13018-018-0873-3

Keywords