Efficacy and efficiency of multivariate linear regression for rapid prediction of femoral strain fields during activity.

Multivariate Linear Regression-based (MLR) surrogate models were explored to reduce the computational cost of predicting femoral strains during normal activity in comparison with finite element analysis. The musculoskeletal model of one individual, the finite-element model of the right femur, and experimental force and motion data for normal walking, fast walking, stair ascent, stair descent, and rising from a chair were obtained from a previous study. Equivalent Von Mises strain was calculated for 1000 frames uniformly distributed across activities. MLR surrogate models were generated using training sets of 50, 100, 200 and 300 samples. The finite-element and MLR analyses were compared using linear regression. The Root Mean Square Error (RMSE) and the 95th percentile of the strain error distribution were used as indicators of average and peak error. The MLR model trained using 200 samples (RMSE < 108 µε; peak error < 228 µε) was used as a reference. The finite-element method required 66 s per frame on a standard desktop computer. The MLR model required 0.1 s per frame plus 1848 s of training time. RMSE ranged from 1.2% to 1.3% while peak error ranged from 2.2% to 3.6% of the maximum micro-strain (5020 µε). Performance within an activity was lower during early and late stance, with RMSE of 4.1% and peak error of 8.6% of the maximum computed micro-strain. These results show that MLR surrogate models may be used to rapidly and accurately estimate strain fields in long bones during daily physical activity.

[1]  Tim W Dorn,et al.  Muscular strategy shift in human running: dependence of running speed on hip and ankle muscle performance , 2012, Journal of Experimental Biology.

[2]  M. Pandy,et al.  Mechanical Loading of the Femoral Neck in Human Locomotion , 2018, Journal of bone and mineral research : the official journal of the American Society for Bone and Mineral Research.

[3]  Marcus G Pandy,et al.  Dynamic simulation of knee-joint loading during gait using force-feedback control and surrogate contact modelling. , 2017, Medical engineering & physics.

[4]  Ayman Habib,et al.  OpenSim: Open-Source Software to Create and Analyze Dynamic Simulations of Movement , 2007, IEEE Transactions on Biomedical Engineering.

[5]  Mark Taylor,et al.  Computationally efficient prediction of bone-implant interface micromotion of a cementless tibial tray during gait. , 2014, Journal of biomechanics.

[6]  Kozo Nakamura,et al.  Prediction of proximal femur strength using a CT-based nonlinear finite element method: differences in predicted fracture load and site with changing load and boundary conditions. , 2009, Bone.

[7]  Diogo M. Geraldes,et al.  Consideration of multiple load cases is critical in modelling orthotropic bone adaptation in the femur , 2015, Biomechanics and Modeling in Mechanobiology.

[8]  Barbara Reggiani,et al.  Finite-Element Modeling of Bones From CT Data: Sensitivity to Geometry and Material Uncertainties , 2006, IEEE Transactions on Biomedical Engineering.

[9]  Mark Taylor,et al.  A Computational Efficient Method to Assess the Sensitivity of Finite-Element Models: An Illustration With the Hemipelvis. , 2016, Journal of biomechanical engineering.

[10]  G. Niebur,et al.  Comparison of the elastic and yield properties of human femoral trabecular and cortical bone tissue. , 2004, Journal of biomechanics.

[11]  Ahmet Erdemir,et al.  Adaptive surrogate modeling for efficient coupling of musculoskeletal control and tissue deformation models. , 2009, Journal of biomechanical engineering.

[12]  Kozo Nakamura,et al.  Prediction of strength and strain of the proximal femur by a CT-based finite element method. , 2007, Journal of biomechanics.

[13]  Massimo Sartori,et al.  MOtoNMS: A MATLAB toolbox to process motion data for neuromusculoskeletal modeling and simulation , 2015, Source Code for Biology and Medicine.

[14]  Prasanth B. Nair,et al.  Use of a statistical model of the whole femur in a large scale, multi-model study of femoral neck fracture risk. , 2009, Journal of biomechanics.

[15]  M. Pandy,et al.  Sensitivity of model predictions of muscle function to changes in moment arms and muscle-tendon properties: a Monte-Carlo analysis. , 2012, Journal of biomechanics.

[16]  E. Schneider,et al.  Influence of muscle forces on femoral strain distribution. , 1998, Journal of biomechanics.

[17]  Mark Taylor,et al.  Four decades of finite element analysis of orthopaedic devices: where are we now and what are the opportunities? , 2015, Journal of biomechanics.

[18]  G. Faust,et al.  Finite element analysis of four different implants inserted in different positions to stabilize an idealized trochanteric femoral fracture. , 2009, Injury.

[19]  Wei Sun,et al.  A deep learning approach to estimate stress distribution: a fast and accurate surrogate of finite-element analysis , 2018, Journal of The Royal Society Interface.

[20]  Benjamin J Fregly,et al.  Surrogate modeling of deformable joint contact using artificial neural networks. , 2015, Medical engineering & physics.

[21]  Georg N Duda,et al.  Machine learning techniques for the optimization of joint replacements: Application to a short-stem hip implant , 2017, PloS one.

[22]  Peter Pivonka,et al.  Femoral shaft strains during daily activities: Implications for atypical femoral fractures. , 2014, Clinical biomechanics.

[23]  M. Pandy,et al.  Sensitivity of femoral strain calculations to anatomical scaling errors in musculoskeletal models of movement. , 2015, Journal of biomechanics.

[24]  A. Dickinson,et al.  Full-field in vitro measurements and in silico predictions of strain shielding in the implanted femur after total hip arthroplasty , 2015, Proceedings of the Institution of Mechanical Engineers. Part H, Journal of engineering in medicine.

[25]  W. Hayes,et al.  Fracture prediction for the proximal femur using finite element models: Part I--Linear analysis. , 1991, Journal of biomechanical engineering.

[26]  Ju Zhang,et al.  Accuracy of femur reconstruction from sparse geometric data using a statistical shape model , 2017, Computer methods in biomechanics and biomedical engineering.

[27]  Saulo Martelli,et al.  Femoral Neck Strain during Maximal Contraction of Isolated Hip-Spanning Muscle Groups , 2016, Comput. Math. Methods Medicine.

[28]  M. T. Bah,et al.  Efficient computational method for assessing the effects of implant positioning in cementless total hip replacements. , 2011, Journal of biomechanics.

[29]  James C Coburn,et al.  Prediction of contact mechanics in metal-on-metal Total Hip Replacement for parametrically comprehensive designs and loads. , 2015, Journal of biomechanics.

[30]  T. Besier,et al.  Influence of collars on the primary stability of cementless femoral stems: A finite element study using a diverse patient cohort , 2017, Journal of orthopaedic research : official publication of the Orthopaedic Research Society.

[31]  Raphael T Haftka,et al.  Surrogate articular contact models for computationally efficient multibody dynamic simulations. , 2010, Medical engineering & physics.

[32]  W. Hayes,et al.  Fracture prediction for the proximal femur using finite element models: Part II--Nonlinear analysis. , 1991, Journal of biomechanical engineering.

[33]  J. Favre,et al.  A neural network model to predict knee adduction moment during walking based on ground reaction force and anthropometric measurements. , 2012, Journal of biomechanics.

[34]  E. Perilli,et al.  Development of a surrogate model based on patient weight, bone mass and geometry to predict femoral neck strains and fracture loads. , 2017, Journal of biomechanics.

[35]  Raphael T Haftka,et al.  Two-dimensional surrogate contact modeling for computationally efficient dynamic simulation of total knee replacements. , 2009, Journal of biomechanical engineering.

[36]  R. Brand,et al.  Testing the daily stress stimulus theory of bone adaptation with natural and experimentally controlled strain histories. , 1997, Journal of biomechanics.

[37]  Monica Reggiani,et al.  Biofeedback for Gait Retraining Based on Real-Time Estimation of Tibiofemoral Joint Contact Forces , 2017, IEEE Transactions on Neural Systems and Rehabilitation Engineering.

[38]  Marco Viceconti,et al.  Subject-specific finite element models of long bones: An in vitro evaluation of the overall accuracy. , 2006, Journal of biomechanics.

[39]  Marco Viceconti,et al.  Sensitivity of a subject-specific musculoskeletal model to the uncertainties on the joint axes location , 2015, Computer methods in biomechanics and biomedical engineering.