1
|
Princelle D, Davico G, Viceconti M. Comparative validation of two patient-specific modelling pipelines for predicting knee joint forces during level walking. J Biomech 2023; 159:111758. [PMID: 37659354 DOI: 10.1016/j.jbiomech.2023.111758] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/29/2022] [Revised: 07/28/2023] [Accepted: 08/07/2023] [Indexed: 09/04/2023]
Abstract
Over the past few years, the use of computer models and simulations tailored to the patient's physiology to assist clinical decision-making has increased enormously.While several pipelines to develop personalized models exist, their adoption on a large scale is still limited due to the required niche computational skillset and the lengthy operations required. Novel toolboxes, such as STAPLE, promise to streamline and expedite the development of image-based skeletal lower limb models. STAPLE-generated models can be rapidly generated, with minimal user input, and present similar joint kinematics and kinetics compared to models developed employing the established INSIGNEO pipeline. Yet, it is unclear how much the observed discrepancies scale up and affect joint contact force predictions. In this study, we compared image-based musculoskeletal models developed (i) with the INSIGNEO pipeline and (ii) with a semi-automated pipeline that combines STAPLE and nmsBuilder, and assessed their accuracy against experimental implant data.Our results showed that both pipelines predicted similar total knee joint contact forces between one another in terms of profiles and average values, characterized by a moderately high level of agreement with the experimental data. Nonetheless, the Student t-test revealed statistically significant differences between both pipelines. Of note, the STAPLE-based pipeline required considerably less time than the INSIGNEO pipeline to generate a musculoskeletal model (i.e., 60 vs 160 min). This is likely to open up opportunities for the use of personalized musculoskeletal models in clinical practice, where time is of the essence.
Collapse
Affiliation(s)
- Domitille Princelle
- Medical Technology Lab, IRCCS Istituto Ortopedico Rizzoli, Bologna, Italy; Department of Industrial Engineering, Alma Mater Studiorum - University of Bologna, Italy.
| | - Giorgio Davico
- Medical Technology Lab, IRCCS Istituto Ortopedico Rizzoli, Bologna, Italy; Department of Industrial Engineering, Alma Mater Studiorum - University of Bologna, Italy.
| | - Marco Viceconti
- Medical Technology Lab, IRCCS Istituto Ortopedico Rizzoli, Bologna, Italy; Department of Industrial Engineering, Alma Mater Studiorum - University of Bologna, Italy
| |
Collapse
|
2
|
Valente G, Benedetti MG, De Paolis M, Donati DM, Taddei F. Differences in hip musculoskeletal loads between limbs during daily activities in patients with 3D-printed hemipelvic reconstructions following tumor surgery. Gait Posture 2023; 102:56-63. [PMID: 36924596 DOI: 10.1016/j.gaitpost.2023.03.005] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 08/05/2022] [Revised: 03/07/2023] [Accepted: 03/09/2023] [Indexed: 03/18/2023]
Abstract
BACKGROUND Anatomical custom-made prostheses, thanks to computer-aided design and 3D-printing technology, help improve osseointegration and reduce mechanical complications in bone reconstructions following bone tumors. A recent quantitative analysis of long-term recovery in patients with 3D-printed reconstructions following pelvic tumor surgery showed asymmetries in ground reaction forces between limbs during different motor activities, while standing very good motor performance and quality of life. RESEARCH QUESTION We analyzed hip contact forces and muscle forces in that cohort of six patients with an innovative custom-made reconstruction of the hemipelvis, and we tested the hypothesis that asymmetries in ground reaction forces would result in more marked differences in musculoskeletal forces. METHODS State-of-the-art musculoskeletal modeling in an optimization-based inverse-dynamics workflow was used to calculate hip contact forces and muscle forces during five motor activities, and the differences between limbs were statistically evaluated across the motor activity cycles and on the force peaks. RESULTS The musculoskeletal loads were found to be not symmetric, as hip loads were generally higher in the contralateral limb. We found significant differences in considerable portions of the motor activities cycles except squat, load symmetry indices indicating a load increase (median up to 25%) on the contralateral limb, especially during stair descent and chair rise/sit, and significantly higher values in the contralateral limb at force peaks. SIGNIFICANCE We confirmed the hypothesis that residual asymmetries found in ground reaction forces were amplified when hip musculoskeletal loads were investigated, reflecting a shift of the loads toward the intact limb. Despite the general trend of higher loads found in the contralateral hip, this cannot be considered a risk of overloading, as both hips supported loads in a physiological range or lower, indicating a likely optimal recovery.
Collapse
Affiliation(s)
- Giordano Valente
- Bioengineering and Computing Laboratory, IRCCS Istituto Ortopedico Rizzoli, Bologna, Italy.
| | - Maria Grazia Benedetti
- Physical Medicine and Rehabilitation Unit, IRCCS Istituto Ortopedico Rizzoli, Bologna, Italy
| | - Massimiliano De Paolis
- Department of Orthopaedics, IRCCS Azienda Ospedaliero-Universitaria di Bologna, Bologna, Italy
| | | | - Fulvia Taddei
- Bioengineering and Computing Laboratory, IRCCS Istituto Ortopedico Rizzoli, Bologna, Italy
| |
Collapse
|
3
|
Real-time replication of three-dimensional and time-varying physiological loading cycles for bone and implant testing: A novel protocol demonstrated for the proximal human femur while walking. J Mech Behav Biomed Mater 2021; 124:104817. [PMID: 34536802 DOI: 10.1016/j.jmbbm.2021.104817] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/04/2021] [Revised: 08/26/2021] [Accepted: 09/03/2021] [Indexed: 11/24/2022]
Abstract
In vitro real-time replication of three-dimensional, time-varying load profiles acting on human bones during physical activity can advance bone and implant testing protocols. This study aimed to develop a novel protocol for applying the three-dimensional, time-varying hip contact force while walking to a human femur specimen. The target force profile was obtained from the literature. A proximal femur from an elderly female donor was instrumented using ten rosette strain gages and tested using a custom-made hexapod robot. A load-control algorithm determined the robot position generating the target force at low frequency (0.0004 Hz). Five cycles of the robot position were played back at five intermediate frequencies up to real-time (0.04, 0.08, 0.16, 0.4, and 0.8 Hz). The hip reaction force, the length of the actuators (position), and cortical strains were compared. The error in the load-control force was 0.3 ± 4.2 N (mean ± SD). The last three force, position, and strain cycles varied by less than 1.1% for every frequency analyzed. Across frequencies, the force increased by 28% at 0.8 Hz as a logarithmic function of frequency (R2 = 0.98). The position and strain error linearly increased with frequency up to 0.4 Hz. The median position error and the interquartile range of the strain error reached 15% and 13% at 0.8 Hz. Changes of force and cortical strain at increasing frequencies were linearly related (R2 = 0.99). Therefore, the protocol developed can provide repeatable three-dimensional time-varying load profiles, although the comparison of the specimen deformation obtained across frequencies should be considered with care, particularly in the higher frequency range. This information supports the design of dynamic tests of bone and implants.
Collapse
|
4
|
Melzner M, Süß F, Dendorfer S. The impact of anatomical uncertainties on the predictions of a musculoskeletal hand model - a sensitivity study. Comput Methods Biomech Biomed Engin 2021; 25:156-164. [PMID: 34180730 DOI: 10.1080/10255842.2021.1940974] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/21/2022]
Abstract
Outputs of musculoskeletal models should be considered probabilistic rather than deterministic as they are affected by inaccuracies and estimations associated with the development of the model. One of these uncertainties being critical for modeling arises from the determination of the muscles' line of action and the physiological cross-sectional area. Therefore, the aim of this study was to evaluate the outcome sensitivity of model predictions from a musculoskeletal hand model in comparison to the uncertainty of these input parameters. For this purpose, the kinematics and muscle activities of different hand movements (abduction of the fingers, abduction of the thumb, and flexion of the thumb) were recorded. One thousand simulations were calculated for each movement using the Latin hypercube sampling method with a corresponding variation of the muscle origin/insertion points and the cross-sectional area. Comparing the standard hand to simulations incorporating uncertainties of input parameters shows no major deviations in on- and off-set time point of muscle activities. About 60% of simulations are located within a ± 30% interval around the standard model concerning joint reaction forces. The comparison with the variation of the input data leads to the conclusion that the standard hand model is able to provide not over-scattered outcomes and, therefore, can be considered relatively stable. These results are of practical importance to the personalization of a musculoskeletal model with subject-specific bone geometries and hence changed muscle line of action.
Collapse
Affiliation(s)
- Maximilian Melzner
- Laboratory for Biomechanics, Ostbayerische Technische Hochschule (OTH) Regensburg, Regensburg, Germany.,Regensburg Center of Biomedical Engineering, OTH and University Regensburg, Regensburg, Germany
| | - Franz Süß
- Laboratory for Biomechanics, Ostbayerische Technische Hochschule (OTH) Regensburg, Regensburg, Germany.,Regensburg Center of Biomedical Engineering, OTH and University Regensburg, Regensburg, Germany
| | - Sebastian Dendorfer
- Laboratory for Biomechanics, Ostbayerische Technische Hochschule (OTH) Regensburg, Regensburg, Germany.,Regensburg Center of Biomedical Engineering, OTH and University Regensburg, Regensburg, Germany
| |
Collapse
|
5
|
A marker registration method to improve joint angles computed by constrained inverse kinematics. PLoS One 2021; 16:e0252425. [PMID: 34048476 PMCID: PMC8162579 DOI: 10.1371/journal.pone.0252425] [Citation(s) in RCA: 7] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/09/2021] [Accepted: 05/15/2021] [Indexed: 11/21/2022] Open
Abstract
Accurate computation of joint angles from optical marker data using inverse kinematics methods requires that the locations of markers on a model match the locations of experimental markers on participants. Marker registration is the process of positioning the model markers so that they match the locations of the experimental markers. Markers are typically registered using a graphical user interface (GUI), but this method is subjective and may introduce errors and uncertainty to the calculated joint angles and moments. In this investigation, we use OpenSim to isolate and quantify marker registration–based error from other sources of error by analyzing the gait of a bipedal humanoid robot for which segment geometry, mass properties, and joint angles are known. We then propose a marker registration method that is informed by the orientation of anatomical reference frames derived from surface-mounted optical markers as an alternative to user registration using a GUI. The proposed orientation registration method reduced the average root-mean-square error in both joint angles and joint moments by 67% compared to the user registration method, and eliminated variability among users. Our results show that a systematic method for marker registration that reduces subjective user input can make marker registration more accurate and repeatable.
Collapse
|
6
|
Automatic generation of personalised skeletal models of the lower limb from three-dimensional bone geometries. J Biomech 2020; 116:110186. [PMID: 33515872 DOI: 10.1016/j.jbiomech.2020.110186] [Citation(s) in RCA: 18] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/24/2020] [Revised: 10/06/2020] [Accepted: 12/11/2020] [Indexed: 02/07/2023]
Abstract
The generation of personalised and patient-specific musculoskeletal models is currently a cumbersome and time-consuming task that normally requires several processing hours and trained operators. We believe that this aspect discourages the use of computational models even when appropriate data are available and personalised biomechanical analysis would be beneficial. In this paper we present a computational tool that enables the fully automatic generation of skeletal models of the lower limb from three-dimensional bone geometries, normally obtained by segmentation of medical images. This tool was evaluated against four manually created lower limb models finding remarkable agreement in the computed joint parameters, well within human operator repeatability. The coordinate systems origins were identified with maximum differences between 0.5 mm (hip joint) and 5.9 mm (subtalar joint), while the joint axes presented discrepancies between 1° (knee joint) to 11° (subtalar joint). To prove the robustness of the methodology, the models were built from four datasets including both genders, anatomies ranging from juvenile to elderly and bone geometries reconstructed from high-quality computed tomography as well as lower-quality magnetic resonance imaging scans. The entire workflow, implemented in MATLAB scripting language, executed in seconds and required no operator intervention, creating lower extremity models ready to use for kinematic and kinetic analysis or as baselines for more advanced musculoskeletal modelling approaches, of which we provide some practical examples. We auspicate that this technical advancement, together with upcoming progress in medical image segmentation techniques, will promote the use of personalised models in larger-scale studies than those hitherto undertaken.
Collapse
|
7
|
A multi-scale modelling framework combining musculoskeletal rigid-body simulations with adaptive finite element analyses, to evaluate the impact of femoral geometry on hip joint contact forces and femoral bone growth. PLoS One 2020; 15:e0235966. [PMID: 32702015 PMCID: PMC7377390 DOI: 10.1371/journal.pone.0235966] [Citation(s) in RCA: 33] [Impact Index Per Article: 8.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/03/2020] [Accepted: 06/25/2020] [Indexed: 11/23/2022] Open
Abstract
Multi-scale simulations, combining muscle and joint contact force (JCF) from musculoskeletal simulations with adaptive mechanobiological finite element analysis, allow to estimate musculoskeletal loading and predict femoral growth in children. Generic linearly scaled musculoskeletal models are commonly used. This approach, however, neglects subject- and age-specific musculoskeletal geometry, e.g. femoral neck-shaft angle (NSA) and anteversion angle (AVA). This study aimed to evaluate the impact of proximal femoral geometry, i.e. altered NSA and AVA, on hip JCF and femoral growth simulations. Musculoskeletal models with NSA ranging from 120° to 150° and AVA ranging from 20° to 50° were created and used to calculate muscle and hip JCF based on the gait analysis data of a typically developing child. A finite element model of a paediatric femur was created from magnetic resonance images. The finite element model was morphed to the geometries of the different musculoskeletal models and used for mechanobiological finite element analysis to predict femoral growth trends. Our findings showed that hip JCF increase with increasing NSA and AVA. Furthermore, the orientation of the hip JCF followed the orientation of the femoral neck axis. Consequently, the osteogenic index, which is a function of cartilage stresses and defines the growth rate, barely changed with altered NSA and AVA. Nevertheless, growth predictions were sensitive to the femoral geometry due to changes in the predicted growth directions. Altered NSA had a bigger impact on the growth results than altered AVA. Growth simulations based on mechanobiological principles were in agreement with reported changes in paediatric populations.
Collapse
|
8
|
Machine learning methods to support personalized neuromusculoskeletal modelling. Biomech Model Mechanobiol 2020; 19:1169-1185. [DOI: 10.1007/s10237-020-01367-8] [Citation(s) in RCA: 15] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/22/2020] [Accepted: 07/08/2020] [Indexed: 12/19/2022]
|
9
|
Hiasa Y, Otake Y, Takao M, Ogawa T, Sugano N, Sato Y. Automated Muscle Segmentation from Clinical CT Using Bayesian U-Net for Personalized Musculoskeletal Modeling. IEEE TRANSACTIONS ON MEDICAL IMAGING 2020; 39:1030-1040. [PMID: 31514128 DOI: 10.1109/tmi.2019.2940555] [Citation(s) in RCA: 52] [Impact Index Per Article: 13.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/10/2023]
Abstract
We propose a method for automatic segmentation of individual muscles from a clinical CT. The method uses Bayesian convolutional neural networks with the U-Net architecture, using Monte Carlo dropout that infers an uncertainty metric in addition to the segmentation label. We evaluated the performance of the proposed method using two data sets: 20 fully annotated CTs of the hip and thigh regions and 18 partially annotated CTs that are publicly available from The Cancer Imaging Archive (TCIA) database. The experiments showed a Dice coefficient (DC) of 0.891±0.016 (mean±std) and an average symmetric surface distance (ASD) of 0.994±0.230 mm over 19 muscles in the set of 20 CTs. These results were statistically significant improvements compared to the state-of-the-art hierarchical multi-atlas method which resulted in 0.845 ± 0.031 DC and 1.556 ± 0.444 mm ASD. We evaluated validity of the uncertainty metric in the multi-class organ segmentation problem and demonstrated a correlation between the pixels with high uncertainty and the segmentation failure. One application of the uncertainty metric in active-learning is demonstrated, and the proposed query pixel selection method considerably reduced the manual annotation cost for expanding the training data set. The proposed method allows an accurate patient-specific analysis of individual muscle shapes in a clinical routine. This would open up various applications including personalization of biomechanical simulation and quantitative evaluation of muscle atrophy.
Collapse
|
10
|
An Anatomical-Based Subject-Specific Model of In-Vivo Knee Joint 3D Kinematics From Medical Imaging. APPLIED SCIENCES-BASEL 2020. [DOI: 10.3390/app10062100] [Citation(s) in RCA: 18] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/11/2022]
Abstract
Biomechanical models of the knee joint allow the development of accurate procedures as well as novel devices to restore the joint natural motion. They are also used within musculoskeletal models to perform clinical gait analysis on patients. Among relevant knee models in the literature, the anatomy-based spatial parallel mechanisms represent the joint motion using rigid links for the ligaments’ isometric fibres and point contacts for the articular surfaces. To customize analyses, therapies and devices, there is the need to define subject-specific models, but relevant procedures and their accuracy are still questioned. A procedure is here proposed and validated to define a customized knee model based on a spatial parallel mechanism. Computed tomography, magnetic resonance and 3D-video-fluoroscopy were performed on a healthy volunteer to define the personalized model geometry. The model was then validated by comparing the measured and the replicated joint motion. The model showed mean absolute difference and standard deviations in translations and rotations, respectively of 0.98 ± 0.40 mm and 0.68 ± 0.29 ° for the tibia–femur motion, and of 0.77 ± 0.15 mm and 2.09 ± 0.69 ° for the patella–femur motion. These results show that accurate personalized spatial models of knee kinematics can be obtained from in-vivo imaging.
Collapse
|
11
|
Zuk M, Syczewska M, Pezowicz C. Sensitivity analysis of the estimated muscle forces during gait with respect to the musculoskeletal model parameters and dynamic simulation techniques. J Biomech Eng 2018; 140:2694845. [PMID: 30098142 DOI: 10.1115/1.4040943] [Citation(s) in RCA: 13] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/22/2017] [Indexed: 11/08/2022]
Abstract
The purpose of the current study was to investigate the robustness of dynamic simulation results in the presence of uncertainties resulting from application of a scaled-generic musculoskeletal model instead of a subject-specific model as well as the effect of the choice of simulation method on the obtained muscle forces. The performed sensitivity analysis consisted of the following multibody parameter modifications: maximum isometric muscle forces, number of muscles, the hip joint centre location, segment masses as well as different dynamic simulation methods, namely static optimization with three different criteria and a computed muscle control algorithm (hybrid approach combining forward and inverse dynamics). Twenty-four different models and fifty-five resultant dynamic simulation data sets were analysed. The effects of model perturbation on the magnitude and profile of muscle forces were compared. It has been shown that estimated muscle forces are very sensitive to model parameters. The greatest impact was observed in the case of the force magnitude of the muscles generating high forces during gait (regardless of the modification introduced). However, the force profiles of those muscles were preserved. Relatively large differences in muscle forces were observed for different simulation techniques, which included both magnitude and profile of muscle forces. Personalization of model parameters would affect the resultant muscle forces and seems to be necessary to improve general accuracy of the estimated parameters. However, personalization alone will not ensure high accuracy due to the still unresolved muscle force sharing problem.
Collapse
Affiliation(s)
- Magdalena Zuk
- Faculty of Mechanical Engineering, Wrocław University of Science and Technology, Wrocław, Poland
| | - Malgorzata Syczewska
- Department of Paediatric Rehabilitation, The Children's Memorial Health Institute, Warsaw, Poland
| | - Celina Pezowicz
- Faculty of Mechanical Engineering, Wrocław University of Science and Technology, Wrocław, Poland
| |
Collapse
|
12
|
Modenese L, Montefiori E, Wang A, Wesarg S, Viceconti M, Mazzà C. Investigation of the dependence of joint contact forces on musculotendon parameters using a codified workflow for image-based modelling. J Biomech 2018; 73:108-118. [DOI: 10.1016/j.jbiomech.2018.03.039] [Citation(s) in RCA: 50] [Impact Index Per Article: 8.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/04/2017] [Revised: 02/09/2018] [Accepted: 03/21/2018] [Indexed: 11/24/2022]
|
13
|
Valente G, Crimi G, Vanella N, Schileo E, Taddei F. nmsBuilder: Freeware to create subject-specific musculoskeletal models for OpenSim. COMPUTER METHODS AND PROGRAMS IN BIOMEDICINE 2017; 152:85-92. [PMID: 29054263 DOI: 10.1016/j.cmpb.2017.09.012] [Citation(s) in RCA: 70] [Impact Index Per Article: 10.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 04/05/2017] [Revised: 08/05/2017] [Accepted: 09/15/2017] [Indexed: 06/07/2023]
Abstract
BACKGROUND AND OBJECTIVE Musculoskeletal modeling and simulations of movement have been increasingly used in orthopedic and neurological scenarios, with increased attention to subject-specific applications. In general, musculoskeletal modeling applications have been facilitated by the development of dedicated software tools; however, subject-specific studies have been limited also by time-consuming modeling workflows and high skilled expertise required. In addition, no reference tools exist to standardize the process of musculoskeletal model creation and make it more efficient. Here we present a freely available software application, nmsBuilder 2.0, to create musculoskeletal models in the file format of OpenSim, a widely-used open-source platform for musculoskeletal modeling and simulation. nmsBuilder 2.0 is the result of a major refactoring of a previous implementation that moved a first step toward an efficient workflow for subject-specific model creation. METHODS nmsBuilder includes a graphical user interface that provides access to all functionalities, based on a framework for computer-aided medicine written in C++. The operations implemented can be used in a workflow to create OpenSim musculoskeletal models from 3D surfaces. A first step includes data processing to create supporting objects necessary to create models, e.g. surfaces, anatomical landmarks, reference systems; and a second step includes the creation of OpenSim objects, e.g. bodies, joints, muscles, and the corresponding model. RESULTS We present a case study using nmsBuilder 2.0: the creation of an MRI-based musculoskeletal model of the lower limb. The model included four rigid bodies, five degrees of freedom and 43 musculotendon actuators, and was created from 3D surfaces of the segmented images of a healthy subject through the modeling workflow implemented in the software application. CONCLUSIONS We have presented nmsBuilder 2.0 for the creation of musculoskeletal OpenSim models from image-based data, and made it freely available via nmsbuilder.org. This application provides an efficient workflow for model creation and helps standardize the process. We hope this would help promote personalized applications in musculoskeletal biomechanics, including larger sample size studies, and might also represent a basis for future developments for specific applications.
Collapse
Affiliation(s)
- Giordano Valente
- Medical Technology Laboratory, Rizzoli Orthopaedic Institute, via di Barbiano 1/10, 40136 Bologna, Italy.
| | - Gianluigi Crimi
- Medical Technology Laboratory, Rizzoli Orthopaedic Institute, via di Barbiano 1/10, 40136 Bologna, Italy
| | - Nicola Vanella
- Medical Technology Laboratory, Rizzoli Orthopaedic Institute, via di Barbiano 1/10, 40136 Bologna, Italy
| | - Enrico Schileo
- Computational Bioengineering Laboratory, Rizzoli Orthopaedic Institute, via di Barbiano 1/10, 40136 Bologna, Italy
| | - Fulvia Taddei
- Medical Technology Laboratory, Rizzoli Orthopaedic Institute, via di Barbiano 1/10, 40136 Bologna, Italy
| |
Collapse
|
14
|
Valente G, Pitto L, Schileo E, Piroddi S, Leardini A, Manfrini M, Taddei F. Relationship between bone adaptation and in-vivo mechanical stimulus in biological reconstructions after bone tumor: A biomechanical modeling analysis. Clin Biomech (Bristol, Avon) 2017; 42:99-107. [PMID: 28131017 DOI: 10.1016/j.clinbiomech.2017.01.017] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 07/07/2016] [Revised: 01/17/2017] [Accepted: 01/19/2017] [Indexed: 02/07/2023]
Abstract
BACKGROUND Biomechanical interpretations of bone adaptation in biological reconstructions following bone tumors would be crucial for orthopedic oncologists, particularly if based on quantitative observations. This would help plan for surgical treatments, rehabilitative programs and communication with the patients. We aimed to analyze the biomechanical adaptation of a femoral reconstruction after Ewing sarcoma according to an increasingly-used surgical technique, and to relate in-progress bone resorption to the mechanical stimulus induced by different motor activities. METHODS We created a multiscale musculoskeletal and finite element model from CT scans and motion analysis data at a 76-month follow-up of a patient, to analyze muscle and joint loads, and to compare the mechanical competence of the reconstructed bone with the contralateral limb, in the current real condition and in a possible revision surgery that removed proximal screws. FINDINGS Our results showed strategies of muscle coordination that led to differences in joint loads between limbs more marked in more demanding motor activities, and generally larger in the contralateral limb. The operated femur presented a markedly low ratio of physiological strain due to load-sharing with the metal implant, particularly in the lateral aspect. The possible revision surgery would help restore a physiological strain configuration, while the safety of the reconstruction would not be threatened. INTERPRETATION We suggest that bone resorption is related to load-sharing and to the internal forces exerted during movement, and the mechanical stimulus should be improved by adopting modifications in the surgical treatment and by promoting physical therapy aimed at specific muscle strengthening.
Collapse
Affiliation(s)
- Giordano Valente
- Medical Technology Laboratory, Rizzoli Orthopaedic Institute, Bologna, Italy.
| | - Lorenzo Pitto
- Medical Technology Laboratory, Rizzoli Orthopaedic Institute, Bologna, Italy
| | - Enrico Schileo
- Computational Bioengineering Laboratory, Rizzoli Orthopaedic Institute, Bologna, Italy
| | - Sabina Piroddi
- Medical Technology Laboratory, Rizzoli Orthopaedic Institute, Bologna, Italy
| | - Alberto Leardini
- Movement Analysis Laboratory, Rizzoli Orthopaedic Institute, Bologna, Italy
| | - Marco Manfrini
- Orthopedic and Traumatologic Clinic for Musculoskeletal Tumors, Rizzoli Orthopaedic Institute, Bologna, Italy
| | - Fulvia Taddei
- Medical Technology Laboratory, Rizzoli Orthopaedic Institute, Bologna, Italy
| |
Collapse
|
15
|
Weese J, Lorenz C. Four challenges in medical image analysis from an industrial perspective. Med Image Anal 2016; 33:44-49. [PMID: 27344939 DOI: 10.1016/j.media.2016.06.023] [Citation(s) in RCA: 38] [Impact Index Per Article: 4.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/28/2016] [Revised: 06/14/2016] [Accepted: 06/15/2016] [Indexed: 12/12/2022]
Abstract
Today's medical imaging systems produce a huge amount of images containing a wealth of information. However, the information is hidden in the data and image analysis algorithms are needed to extract it, to make it readily available for medical decisions and to enable an efficient work flow. Advances in medical image analysis over the past 20 years mean there are now many algorithms and ideas available that allow to address medical image analysis tasks in commercial solutions with sufficient performance in terms of accuracy, reliability and speed. At the same time new challenges have arisen. Firstly, there is a need for more generic image analysis technologies that can be efficiently adapted for a specific clinical task. Secondly, efficient approaches for ground truth generation are needed to match the increasing demands regarding validation and machine learning. Thirdly, algorithms for analyzing heterogeneous image data are needed. Finally, anatomical and organ models play a crucial role in many applications, and algorithms to construct patient-specific models from medical images with a minimum of user interaction are needed. These challenges are complementary to the on-going need for more accurate, more reliable and faster algorithms, and dedicated algorithmic solutions for specific applications.
Collapse
Affiliation(s)
- Jürgen Weese
- Philips Research Hamburg, Röntgenstrasse 22 - 24, D-22335 Hamburg, Germany.
| | - Cristian Lorenz
- Philips Research Hamburg, Röntgenstrasse 22 - 24, D-22335 Hamburg, Germany.
| |
Collapse
|
16
|
Carbone V, van der Krogt M, Koopman H, Verdonschot N. Sensitivity of subject-specific models to Hill muscle–tendon model parameters in simulations of gait. J Biomech 2016; 49:1953-1960. [DOI: 10.1016/j.jbiomech.2016.04.008] [Citation(s) in RCA: 31] [Impact Index Per Article: 3.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/29/2015] [Revised: 03/31/2016] [Accepted: 04/02/2016] [Indexed: 11/16/2022]
|
17
|
Effect of lower-limb joint models on subject-specific musculoskeletal models and simulations of daily motor activities. J Biomech 2015; 48:4198-205. [PMID: 26506255 DOI: 10.1016/j.jbiomech.2015.09.042] [Citation(s) in RCA: 30] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/27/2015] [Revised: 09/23/2015] [Accepted: 09/26/2015] [Indexed: 11/21/2022]
Abstract
Understanding the validity of using musculoskeletal models is critical, making important to assess how model parameters affect predictions. In particular, assumptions on joint models can affect predictions from simulations of movement, and the identification of image-based joints is unavoidably affected by uncertainty that can decrease the benefits of increasing model complexity. We evaluated the effect of different lower-limb joint models on muscle and joint contact forces during four motor tasks, and assessed the sensitivity to the uncertainties in the identification of anatomical four-bar-linkage joints. Three MRI-based musculoskeletal models having different knee and ankle joint models were created and used for the purpose. Model predictions were compared against a baseline model including simpler and widely-adopted joints. In addition, a probabilistic analysis was performed by perturbing four-bar-linkage joint parameters according to their uncertainty. The differences between models depended on the motor task analyzed, and there could be marked differences at peak loading (up to 2.40 BW at the knee and 1.54 BW at the ankle), although they were rather small over the motor task cycles (up to 0.59 BW at the knee and 0.31 BW at the ankle). The model including more degrees of freedom showed more discrepancies in predicted muscle activations compared to measured muscle activity. Further, including image-based four-bar-linkages was robust to simulate walking, chair rise and stair ascent, but not stair descent (peak standard deviation of 2.66 BW), suggesting that joint model complexity should be set according to the imaging dataset available and the intended application, performing sensitivity analyses.
Collapse
|
18
|
Sensitivity of femoral strain calculations to anatomical scaling errors in musculoskeletal models of movement. J Biomech 2015; 48:3606-15. [DOI: 10.1016/j.jbiomech.2015.08.001] [Citation(s) in RCA: 19] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/21/2014] [Revised: 08/03/2015] [Accepted: 08/04/2015] [Indexed: 11/20/2022]
|
19
|
Martelli S, Calvetti D, Somersalo E, Viceconti M. Stochastic modelling of muscle recruitment during activity. Interface Focus 2015; 5:20140094. [PMID: 25844155 DOI: 10.1098/rsfs.2014.0094] [Citation(s) in RCA: 37] [Impact Index Per Article: 4.1] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/12/2022] Open
Abstract
Muscle forces can be selected from a space of muscle recruitment strategies that produce stable motion and variable muscle and joint forces. However, current optimization methods provide only a single muscle recruitment strategy. We modelled the spectrum of muscle recruitment strategies while walking. The equilibrium equations at the joints, muscle constraints, static optimization solutions and 15-channel electromyography (EMG) recordings for seven walking cycles were taken from earlier studies. The spectrum of muscle forces was calculated using Bayesian statistics and Markov chain Monte Carlo (MCMC) methods, whereas EMG-driven muscle forces were calculated using EMG-driven modelling. We calculated the differences between the spectrum and EMG-driven muscle force for 1-15 input EMGs, and we identified the muscle strategy that best matched the recorded EMG pattern. The best-fit strategy, static optimization solution and EMG-driven force data were compared using correlation analysis. Possible and plausible muscle forces were defined as within physiological boundaries and within EMG boundaries. Possible muscle and joint forces were calculated by constraining the muscle forces between zero and the peak muscle force. Plausible muscle forces were constrained within six selected EMG boundaries. The spectrum to EMG-driven force difference increased from 40 to 108 N for 1-15 EMG inputs. The best-fit muscle strategy better described the EMG-driven pattern (R (2) = 0.94; RMSE = 19 N) than the static optimization solution (R (2) = 0.38; RMSE = 61 N). Possible forces for 27 of 34 muscles varied between zero and the peak muscle force, inducing a peak hip force of 11.3 body-weights. Plausible muscle forces closely matched the selected EMG patterns; no effect of the EMG constraint was observed on the remaining muscle force ranges. The model can be used to study alternative muscle recruitment strategies in both physiological and pathophysiological neuromotor conditions.
Collapse
Affiliation(s)
- Saulo Martelli
- Medical Device Research Institute, School of Computer Science, Engineering and Mathematics , Flinders University , Tonsley, South Australia 5042 , Australia ; North West Academic Centre , The University of Melbourne , St Albans, Victoria 3021 , Australia
| | - Daniela Calvetti
- Department of Mathematics , Applied Mathematics, and Statistics, Case Western Reserve University , Cleveland, OH 44106-7058 , USA
| | - Erkki Somersalo
- Department of Mathematics , Applied Mathematics, and Statistics, Case Western Reserve University , Cleveland, OH 44106-7058 , USA
| | - Marco Viceconti
- Department of Mechanical Engineering and INSIGNEO Institute for in silico Medicine , University of Sheffield , Sheffield S1 3JD , UK
| |
Collapse
|
20
|
TLEM 2.0 – A comprehensive musculoskeletal geometry dataset for subject-specific modeling of lower extremity. J Biomech 2015; 48:734-41. [DOI: 10.1016/j.jbiomech.2014.12.034] [Citation(s) in RCA: 97] [Impact Index Per Article: 10.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Accepted: 11/27/2014] [Indexed: 11/20/2022]
|
21
|
Valente G, Pitto L, Testi D, Seth A, Delp SL, Stagni R, Viceconti M, Taddei F. Are subject-specific musculoskeletal models robust to the uncertainties in parameter identification? PLoS One 2014; 9:e112625. [PMID: 25390896 PMCID: PMC4229232 DOI: 10.1371/journal.pone.0112625] [Citation(s) in RCA: 123] [Impact Index Per Article: 12.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/21/2014] [Accepted: 10/20/2014] [Indexed: 11/22/2022] Open
Abstract
Subject-specific musculoskeletal modeling can be applied to study musculoskeletal disorders, allowing inclusion of personalized anatomy and properties. Independent of the tools used for model creation, there are unavoidable uncertainties associated with parameter identification, whose effect on model predictions is still not fully understood. The aim of the present study was to analyze the sensitivity of subject-specific model predictions (i.e., joint angles, joint moments, muscle and joint contact forces) during walking to the uncertainties in the identification of body landmark positions, maximum muscle tension and musculotendon geometry. To this aim, we created an MRI-based musculoskeletal model of the lower limbs, defined as a 7-segment, 10-degree-of-freedom articulated linkage, actuated by 84 musculotendon units. We then performed a Monte-Carlo probabilistic analysis perturbing model parameters according to their uncertainty, and solving a typical inverse dynamics and static optimization problem using 500 models that included the different sets of perturbed variable values. Model creation and gait simulations were performed by using freely available software that we developed to standardize the process of model creation, integrate with OpenSim and create probabilistic simulations of movement. The uncertainties in input variables had a moderate effect on model predictions, as muscle and joint contact forces showed maximum standard deviation of 0.3 times body-weight and maximum range of 2.1 times body-weight. In addition, the output variables significantly correlated with few input variables (up to 7 out of 312) across the gait cycle, including the geometry definition of larger muscles and the maximum muscle tension in limited gait portions. Although we found subject-specific models not markedly sensitive to parameter identification, researchers should be aware of the model precision in relation to the intended application. In fact, force predictions could be affected by an uncertainty in the same order of magnitude of its value, although this condition has low probability to occur.
Collapse
Affiliation(s)
- Giordano Valente
- Medical Technology Laboratory, Rizzoli Orthopaedic Institute, Bologna, Italy
| | - Lorenzo Pitto
- Medical Technology Laboratory, Rizzoli Orthopaedic Institute, Bologna, Italy
| | - Debora Testi
- BioComputing Competence Centre, SCS s.r.l., Bologna, Italy
| | - Ajay Seth
- Department of Bioengineering, Stanford University, Stanford, California, United States of America
| | - Scott L. Delp
- Department of Bioengineering, Stanford University, Stanford, California, United States of America
- Department of Mechanical Engineering, Stanford University, Stanford, California, United States of America
| | - Rita Stagni
- Department of Electrical, Electronic and Information Engineering, University of Bologna, Bologna, Italy
| | - Marco Viceconti
- Department of Mechanical Engineering and INSIGNEO Institute for In Silico Medicine, University of Sheffield, Sheffield, United Kingdom
| | - Fulvia Taddei
- Medical Technology Laboratory, Rizzoli Orthopaedic Institute, Bologna, Italy
| |
Collapse
|
22
|
Martelli S, Valente G, Viceconti M, Taddei F. Sensitivity of a subject-specific musculoskeletal model to the uncertainties on the joint axes location. Comput Methods Biomech Biomed Engin 2014; 18:1555-63. [DOI: 10.1080/10255842.2014.930134] [Citation(s) in RCA: 43] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/25/2022]
|
23
|
Compensation by nonoperated joints in the lower limbs during walking after endoprosthetic knee replacement following bone tumor resection. Clin Biomech (Bristol, Avon) 2013; 28:898-903. [PMID: 24008057 DOI: 10.1016/j.clinbiomech.2013.08.005] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 04/12/2013] [Revised: 08/14/2013] [Accepted: 08/20/2013] [Indexed: 02/07/2023]
Abstract
BACKGROUND Endoprosthetic knee replacement is often used to preserve joint function in patients with bone tumors of the distal femur or proximal tibia. Recently, because of improved oncologic outcome, surgeons are focusing more on the functional outcome of patients with musculoskeletal tumors. We hypothesized that patients who have undergone endoprosthetic knee replacement are forced to compensate for deficiency in their operated joint during walking. In this study, we investigated differences in gait kinematics, kinetics, and energetics between patients with endoprosthetic knee replacement and healthy subjects. METHODS We performed gait analysis for 8 patients who underwent endoprosthetic knee replacement after bone tumor resection and 8 matched healthy subjects. Gait kinematics, kinetics, and energetics of patients' ipsilateral and contralateral limbs were compared with those of healthy subjects by using Dunnett's test. FINDINGS Compared with healthy subjects, patients showed increased negative joint power around the ipsilateral ankle, greater second peak in the contralateral vertical ground reaction forces, and abnormal hip movement on both sides after initial contact. INTERPRETATION Patients tended to compensate for dysfunction of the reconstructed knee by muscles around the ipsilateral ankle and contralateral hip, with increased load on the contralateral limb during walking. These differences could lead to secondary impairments. Further analysis, including musculoskeletal simulation and assessment of long-term functional outcome with regard to secondary musculoskeletal impairment, is needed to verify the significance of the change in gait and to determine the need for special care for secondary musculoskeletal dysfunction in these patients.
Collapse
|
24
|
Influence of weak hip abductor muscles on joint contact forces during normal walking: probabilistic modeling analysis. J Biomech 2013; 46:2186-93. [PMID: 23891175 DOI: 10.1016/j.jbiomech.2013.06.030] [Citation(s) in RCA: 61] [Impact Index Per Article: 5.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/19/2012] [Revised: 06/28/2013] [Accepted: 06/30/2013] [Indexed: 11/22/2022]
Abstract
The weakness of hip abductor muscles is related to lower-limb joint osteoarthritis, and joint overloading may increase the risk for disease progression. The relationship between muscle strength, structural joint deterioration and joint loading makes the latter an important parameter in the study of onset and follow-up of the disease. Since the relationship between hip abductor weakness and joint loading still remains an open question, the purpose of this study was to adopt a probabilistic modeling approach to give insights into how the weakness of hip abductor muscles, in the extent to which normal gait could be unaltered, affects ipsilateral joint contact forces. A generic musculoskeletal model was scaled to each healthy subject included in the study, and the maximum force-generating capacity of each hip abductor muscle in the model was perturbed to evaluate how all physiologically possible configurations of hip abductor weakness affected the joint contact forces during walking. In general, the muscular system was able to compensate for abductor weakness. The reduced force-generating capacity of the abductor muscles affected joint contact forces to a mild extent, with 50th percentile mean differences up to 0.5 BW (maximum 1.7 BW). There were greater increases in the peak knee joint loads than in loads at the hip or ankle. Gluteus medius, particularly the anterior compartment, was the abductor muscle with the most influence on hip and knee loads. Further studies should assess if these increases in joint loading may affect initiation and progression of osteoarthritis.
Collapse
|
25
|
Chèze L, Moissenet F, Dumas R. State of the art and current limits of musculo-skeletal models for clinical applications. ACTA ACUST UNITED AC 2012. [DOI: 10.1051/sm/2012026] [Citation(s) in RCA: 11] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022]
|
26
|
Fregly BJ, Boninger ML, Reinkensmeyer DJ. Personalized neuromusculoskeletal modeling to improve treatment of mobility impairments: a perspective from European research sites. J Neuroeng Rehabil 2012; 9:18. [PMID: 22463378 PMCID: PMC3342221 DOI: 10.1186/1743-0003-9-18] [Citation(s) in RCA: 48] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/29/2011] [Accepted: 03/30/2012] [Indexed: 05/26/2023] Open
Abstract
Mobility impairments due to injury or disease have a significant impact on quality of life. Consequently, development of effective treatments to restore or replace lost function is an important societal challenge. In current clinical practice, a treatment plan is often selected from a standard menu of options rather than customized to the unique characteristics of the patient. Furthermore, the treatment selection process is normally based on subjective clinical experience rather than objective prediction of post-treatment function. The net result is treatment methods that are less effective than desired at restoring lost function. This paper discusses the possible use of personalized neuromusculoskeletal computer models to improve customization, objectivity, and ultimately effectiveness of treatments for mobility impairments. The discussion is based on information gathered from academic and industrial research sites throughout Europe, and both clinical and technical aspects of personalized neuromusculoskeletal modeling are explored. On the clinical front, we discuss the purpose and process of personalized neuromusculoskeletal modeling, the application of personalized models to clinical problems, and gaps in clinical application. On the technical front, we discuss current capabilities of personalized neuromusculoskeletal models along with technical gaps that limit future clinical application. We conclude by summarizing recommendations for future research efforts that would allow personalized neuromusculoskeletal models to make the greatest impact possible on treatment design for mobility impairments.
Collapse
Affiliation(s)
- Benjamin J Fregly
- Department of Mechanical & Aerospace Engineering, University of Florida, 231 MAE-A Building, PO Box 116250, Gainesville, FL 32611-6250, USA.
| | | | | |
Collapse
|