[Numerical simulation in musculoskeletal biomechanics : Application perspectives and possibilities]. / Numerische Simulation in der muskuloskelettalen Biomechanik : Anwendungsperspektiven und Möglichkeiten.

BACKGROUND: Computational research methods, such as finite element analysis (FEA) and musculoskeletal multi-body simulation (MBS), are important in musculoskeletal biomechanics because they enable a better understanding of the mechanics of the musculoskeletal system, as well as the development and evaluation of orthopaedic implants. These methods are used to analyze clinically relevant issues in various anatomical regions, such as the hip, knee, shoulder joints and spine. Preoperative simulation can improve surgical planning in orthopaedics and predict individual results. EXAMPLES FROM PRACTICE: In this article, the methods of FE analysis and MBS are explained using two practical examples, and the activities of the "Numerical Simulation" cluster of the "Musculoskeletal Biomechanics Research Network (MSB-NET)" are presented in more detail. An outlook classifies numerical simulation in the age of artificial intelligence and draws attention to the relevance of simulation in the (re)approval of implants.

Impact of Structural Compliance of a Six Degree of Freedom Joint Simulator on Virtual Ligament Force Calculation in Total Knee Endoprosthesis Testing.

The AMTI VIVO™ six degree of freedom joint simulator allows reproducible preclinical testing of joint endoprostheses under specific kinematic and loading conditions. When testing total knee endoprosthesis, the articulating femoral and tibial components are each mounted on an actuator with two and four degrees of freedom, respectively. To approximate realistic physiological conditions with respect to soft tissues, the joint simulator features an integrated virtual ligament model that calculates the restoring forces of the ligament apparatus to be applied by the actuators. During joint motion, the locations of the ligament insertion points are calculated depending on both actuators' coordinates. In the present study, we demonstrate that unintended elastic deformations of the actuators due to the specifically high contact forces in the artificial knee joint have a considerable impact on the calculated ligament forces. This study aims to investigate the effect of this structural compliance on experimental results. While the built-in algorithm for calculating the ligament forces cannot be altered by the user, a reduction of the ligament force deviations due to the elastic deformations could be achieved by preloading the articulating implant components in the reference configuration. As a proof of concept, a knee flexion motion with varying ligament conditions was simulated on the VIVO simulator and compared to data derived from a musculoskeletal multibody model of a total knee endoprosthesis.

Influence of posterior cruciate ligament tension on tibiofemoral and patellofemoral joint contact mechanics in cruciate-retaining total knee replacement: a combined musculoskeletal multibody and finite-element simulation.

The influence of posterior cruciate ligament (PCL) tension on the clinical outcome of cruciate-retaining total knee replacement (CR-TKR) remains controversial. Various numerical approaches have been used to study this influence systematically, but the models used are limited by certain assumptions and simplifications. Therefore, the objective of this computational study was to develop a combined musculoskeletal multibody and finite-element simulation during a squat motion to 90° knee flexion with a CR-TKR design to overcome previous limitations regarding model inputs. In addition, different PCL tensions (tight, lax, resected) were modeled and the influence on tibiofemoral and resurfaced patellofemoral joint dynamics and contact stresses was evaluated. The effect of the PCL on knee joint dynamics and contact stresses was more pronounced at higher flexion angles. Tibiofemoral joint dynamics were influenced and a tight PCL induced increased posterior femoral translation during flexion. The maximum contact stress in the tibial insert increased from 20.6 MPa to 22.5 MPa for the resected and tightest PCL at 90° knee flexion. Patellofemoral joint dynamics were only slightly affected by PCL tension. However, the maximum contact stress in the patellar component decreased from 58.0 MPa to 53.7 MPa for the resected and tightest PCL at 90° knee flexion. The combination of musculoskeletal multibody and finite-element simulation is a sufficient method to comprehensively investigate knee joint dynamics and contact stresses in CR-TKR. The PCL tension after CR-TKR affects joint dynamics and contact stresses at the articulating implant surfaces.

Exploration of the Advanced VIVOTM Joint Simulator: An In-Depth Analysis of Opportunities and Limitations Demonstrated by the Artificial Knee Joint.

In biomechanical research, advanced joint simulators such as VIVOTM offer the ability to test artificial joints under realistic kinematics and load conditions. Furthermore, it promises to simplify testing with advanced control approaches and the ability to include virtual ligaments. However, the overall functionality concerning specific test setup conditions, such as the joint lubrication or control algorithm, has not been investigated in-depth so far. Therefore, the aim of this study was to analyse the basic functionality of the VIVOTM joint simulator with six degrees of freedom in order to highlight its capabilities and limitations when testing a total knee endoprostheses using a passive flexion-extension movement. For this, different test setup conditions were investigated, e.g., the control method, repeatability and kinematic reproducibility, waveform frequency, lubrication, and implant embedding. The features offered by the VIVOTM joint simulator are useful for testing joint endoprostheses under realistic loading scenarios. It was found that the results were highly influenced by the varying test setup conditions, although the same mechanical load case was analysed. This study highlights the difficulties encountered when using six degrees of freedom joint simulators, contributes to their understanding, and supports users of advanced joint simulators through functional and tribological analysis of joint endoprostheses.

Do hip resurfacing and short hip stem arthroplasties differ from conventional hip stem replacement regarding impingement-free range of motion?

Total hip joint replacement (THR) is clinically well-established. In this context, the resulting range of motion (ROM) is crucial for patient satisfaction when performing joint movements. However, the ROM for THR with different bone preserving strategies (short hip stem and hip resurfacing) raises the question of whether the ROM is comparable with conventional hip stems. Therefore, this computer-based study aimed to investigate the ROM and type of impingement for different implant systems. An established framework with computer-aided design 3D models based on magnetic resonance imaging data of 19 patients with hip osteoarthritis was used to analyse the ROM for three different implant systems (conventional hip stem vs. short hip stem vs. hip resurfacing) during typical joint movements. Our results revealed that all three designs led to mean maximum flexion higher than 110°. However, hip resurfacing showed less ROM (-5% against conventional and -6% against short hip stem). No significant differences were observed between the conventional and short hip stem during maximum flexion and internal rotation. Contrarily, a significant difference was detected between the conventional hip stem and hip resurfacing during internal rotation (p = 0.003). The ROM of the hip resurfacing was lower than the conventional and short hip stem during all three movements. Furthermore, hip resurfacing shifted the impingement type to implant-to-bone impingement compared with the other implant designs. The calculated ROMs of the implant systems achieved physiological levels during maximum flexion and internal rotation. However, bone impingement was more likely during internal rotation with increasing bone preservation. Despite the larger head diameter of hip resurfacing, the ROM examined was substantially lower than that of conventional and short hip stem.

Patella height influences patellofemoral contact and kinematics following cruciate-retaining total knee replacement.

The role of patella height is discussed controversially in total knee arthroplasty (TKA). Therefore, this computational study aims to systematically analyze the biomechanical effect of different patella heights on patellofemoral (PF) forces and kinematics after cruciate-retaining (CR) TKA. We implemented a CR bicondylar TKA with a dome patellar button in a validated dynamic musculoskeletal multibody model of a male human knee joint. Retropatellar dynamics (contact force [N], shear force [N], patellar shift [mm], tilt [°], and rotation [°]) were evaluated during dual-limb squat motion (flexion from 0° to 90°) with simulated active muscle forces and the effects of different patella heights (Blackburne-Peel [BP] ratio of 0.39, 0.49, 0.65, 0.85, 1.01, and 1.1 were systematically examined). As active knee flexion increased, PF contact force also increased. Patella alta (BP = 1.1) resulted in higher PF contact forces compared to normal patella height (BP = 0.65) by up to 16%. Contrarily, patella baja was associated with decreased PF forces by 7%. Compared to patella baja (BP = 0.39), patella alta (BP = 1.1) considerably increased the contact force by up to 25%. Different patellar heights mainly affected PF shear forces during early knee flexion. Concerning PF kinematics, patella alta (BP = 1.1) yielded a greater lateral tilt of more than 4° and higher patellar rotation by up to 3° during deep knee flexion, compared to normal patella height (BP = 0.65). Our computational study indicates that patella alta is associated with the highest PF contact and shear force after the implantation of a CR bicondylar TKA. This should be considered in PF disorders following TKA.

Computer-based analysis of different component positions and insert thicknesses on tibio-femoral and patello-femoral joint dynamics after cruciate-retaining total knee replacement.

BACKGROUND: Positioning of the implant components and tibial insert thickness constitute critical aspects of total knee replacement (TKR) that influence the postoperative knee joint dynamics. This study aimed to investigate the impact of implant component positioning (anterior-posterior and medio-lateral shift) and varying tibial insert thickness on the tibio-femoral (TF) and patello-femoral (PF) joint kinematics and contact forces after cruciate-retaining (CR)-TKR. METHOD: A validated musculoskeletal multibody simulation (MMBS) model with a fixed-bearing CR-TKR during a squat motion up to 90° knee flexion was deployed to calculate PF and TF joint dynamics for varied implant component positions and tibial insert thicknesses. Evaluation was performed consecutively by comparing the respective knee joint parameters (e.g. contact force, quadriceps muscle force, joint kinematics) to a reference implant position. RESULTS: The PF contact forces were mostly affected by the anterior-posterior as well as medio-lateral positioning of the femoral component (by 3 mm anterior up to 31 % and by 6 mm lateral up to 14 %). TF contact forces were considerably altered by tibial insert thickness (24 % in case of + 4 mm increase) and by the anterior-posterior position of the femoral component (by 3 mm posterior up to 16 %). Concerning PF kinematics, a medialised femoral component by 6 mm increased the lateral patellar tilt by more than 5°. CONCLUSIONS: Our results indicate that regarding PF kinematics and contact forces the positioning of the femoral component was more critical than the tibial component. The positioning of the femoral component in anterior-posterior direction on and PF contact force was evident. Orthopaedic surgeons should strictly monitor the anterior-posterior as well as the medio-lateral position of the femoral component and the insert thickness.

Barefoot walking is more stable in the gait of balance recovery in older adults.

BACKGROUND: Perturbation-based balance training on a treadmill is an emerging method of gait stability training with a characteristic task nature that has had positive and sustained effects on balance recovery strategies and fall reduction. Little is known about the effects produced by shod and barefoot walking. We aimed to investigate which is more appropriate, shod or barefoot walking, for perturbation-based balance training in older adults. METHODS: Fourteen healthy older adults (age: 68.29 ± 3.41 years; body height: 1.76 ± 0.10 m; body mass: 81.14 ± 14.52 kg) performed normal and trip-like perturbed walking trials, shod and barefoot, on a treadmill of the Gait Real-time Analysis Interactive Lab. The marker trajectories data were processed by Human Body Model software embedded in the Gait Offline Analysis Tool. The outcomes of stride length variability, stride time variability, step width variability, and swing time variability were computed and statistically analyzed by a two-way repeated-measures analysis of variance (ANOVA) based on gait pattern (normal gait versus perturbed recovery gait) and footwear condition (shod versus barefoot). RESULTS: Footwear condition effect (p = 0.0310) and gait pattern by footwear condition interaction effect (p = 0.0055) were only observed in swing time variability. Gait pattern effects were detected in all four outcomes of gait variability. CONCLUSIONS: Swing time variability, independent of gait speed, could be a valid indicator to differentiate between footwear conditions. The lower swing time variability in perturbed recovery gait suggests that barefoot walking may be superior to shod walking for perturbation-based balance training in older adults.

Lower extremity joint compensatory effects during the first recovery step following slipping and stumbling perturbations in young and older subjects.

BACKGROUND: The lower extremity may play a crucial role in compensating for gait perturbations. The study aimed to explore the mechanism of perturbation compensation by investigating the gait characteristics and lower extremity joint moment effects in young (YS) and older subjects (OS) during the first recovery gait following slipping (slipping_Rec1) and stumbling (stumbling_Rec1). METHOD: An automatic perturbation-triggered program was developed using D-Flow software based on the Gait Real-time Analysis Interactive Lab to induce the two aforementioned perturbations. Marker trajectories and ground reaction forces were recorded from 15 healthy YS (age: 26.53 ± 3.04 years; body height: 1.73 ± 0.07 m; body mass: 66.81 ± 11.44 kg) and 15 healthy OS (age: 68.33 ± 3.29 years; body height: 1.76 ± 0.10 m; body mass: 81.13 ± 13.99 kg). The Human Body Model was used to compute the variables of interest. One-way analysis of variance and independent samples t-test statistical analyses were performed. RESULTS: In slipping_Rec1 and stumbling_Rec1, the change in gait pattern was mainly reflected in a significant increase in step width, no alterations in step length and stance/swing ratio were revealed. Based on perturbed task specificity, lower extremity joint moments increased or decreased at specific phases of the gait cycle in both YS and OS in slipping_Rec1 and stumbling_Rec1 compared to normal gait. The two perturbed gaits reflected the respective compensatory requirements for the lower extremity joints, with both sagittal and frontal joint moments producing compensatory effects. The aging effect was not reflected in the gait pattern, but rather in the hip extension moment during the initial stance of slipping_Rec1. CONCLUSIONS: Slipping appears to be more demanding for gait recovery than stumbling. Gait perturbation compensatory mechanisms for OS should concentrate on ankle strategy in the frontal plane and counter-rotation strategy around the hip.

Compensatory Responses During Slip-Induced Perturbation in Patients With Knee Osteoarthritis Compared With Healthy Older Adults: An Increased Risk of Falls?

Background: Functional impairment of the knee joint affected by osteoarthritis and loss of muscle strength leads to a significant increase in the number of falls. Nevertheless, little is known about strategies for coping with gait perturbations in patients with knee osteoarthritis (KOA). Thus, this study aimed to examine the compensatory strategies of patients with KOA in response to a backward slip perturbation compared with healthy older adults. Methods: An automated perturbation program was developed by using D-Flow software based on the Gait Real-time Analysis Interactive Lab, and an induced backward slip perturbation was implemented on nine patients with severe KOA (68.89 ± 3.59 years) and 15 age-matched healthy older adults (68.33 ± 3.29 years). Step length, gait speed, range of motion, vertical ground reaction forces, lower extremity joint angles, and joint moments were computed and analyzed. Results: Compared with older adults, patients with KOA had significantly lower step length, gait speed, and vertical ground reaction forces in both normal walking and the first recovery step following backward slip perturbations. Inadequate flexion and extension of joint angles and insufficient generation of joint moments predispose patients with KOA to fall. Hip extension angle and flexion moment, knee range of motion, and vertical ground reaction forces are key monitoring variables. Conclusion: The risk of falls for patients with KOA in response to backward slip perturbations is higher. Patients with KOA should focus not only on quadriceps muscle strength related to knee range of motion but also on improving hip extensor strength and activation through specific exercises. Targeted resistance training and perturbation-based gait training could be better options.

Intraoperative analysis of patellofemoral joint morphology before and after total knee arthroplasty : A feasibility study.

BACKGROUND: Patellofemoral problems after total knee arthroplast (TKA) are frequent and often associated with a change in the geometry of the trochlear groove. OBJECTIVE: The present study aimed to analyze the feasibility of intraoperative examination of the patellofemoral joint geometry before and after the implantation of bicondylar total knee replacements without exposing the patient to radiation. MATERIAL AND METHODS: The patellofemoral joint morphology geometries of 33 patients before and after implantation of a bicondylar total knee replacement was intraoperatively analyzed using a digital scanning method. Femoral surface geometries were extracted from intraoperatively recorded tactile scanning data using an in-house algorithm. The geometries were then characterized by sulcus height, medial femoral condyle height and lateral femoral condyle height. RESULTS: Our preliminary results show that these key geometric parameters are only partially restored after TKA leading to a distinctly different shaped profile of the anterior distal femur. Maximum and minimum mean differences in sulcus height, medial femoral condyle height, and lateral femoral condyle height before and after surgery were 2.00 mm/-3.06 mm, 2.51 mm/-6.25 mm and 2.74 mm/-3.58 mm, respectively. CONCLUSION: A new method for intraoperative analysis of the patellofemoral joint morphology before and after TKA without radiation exposure was developed and utilized. Even with the use of modern total knee designs, the patellofemoral articulation is considerably changed postoperatively as quantified by the key parameters of sulcus height, medial femoral condyle height, and lateral femoral condyle height. This may result in altered knee biomechanics and might explain persistent patellofemoral disorders, which are often reported by patients after TKA.

Kinematics and kinetics comparison of ultra-congruent versus medial-pivot designs for total knee arthroplasty by multibody analysis.

Nowadays, several configurations of total knee arthroplasty (TKA) implants are commercially available whose designs resulted from clinical and biomechanical considerations. Previous research activities led to the development of the so-called medial-pivot (MP) design. However, the actual benefits of the MP, with respect to other prosthesis designs, are still not well understood. The present work compares the impact of two insert geometries, namely the ultra-congruent (UC) and medial-pivot (MP), on the biomechanical behaviour of a bicondylar total knee endoprosthesis. For this purpose, a multibody model of a lower limb was created alternatively integrating the two implants having the insert geometry discretized. Joint dynamics and contact pressure distributions were evaluated by simulating a squat motion. Results showed a similar tibial internal rotation range of about 3.5°, but an early rotation occurs for the MP design. Furthermore, the discretization of the insert geometry allowed to efficiently derive the contact pressure distributions, directly within the multibody simulation framework, reporting peak pressure values of 33 MPa and 20 MPa for the UC and MP, respectively. Clinically, the presented findings confirm the possibility, through a MP design, to achieve a more natural joint kinematics, consequently improving the post-operative patient satisfaction and potentially reducing the occurrence of phenomena leading to the insert loosening.

The Influence of Mathematical Definitions on Patellar Kinematics Representations.

A correlation between patellar kinematics and anterior knee pain is widely accepted. However, there is no consensus on how they are connected or what profile of patellar kinematics would minimize anterior knee pain. Nevertheless, answering this question by merging existing studies is further complicated by the variety of ways to describe patellar kinematics. Therefore, this study describes the most frequently used conventions for defining patellar kinematics, focusing on the rotations. The similarities and differences between the Cardan sequences and angles calculated by projecting axes are analyzed. Additionally, a tool is provided to enable the conversion of kinematic data between definitions in different studies. The choice of convention has a considerable impact on the absolute values and the clinical characteristics of the patello-femoral angles. In fact, the angles that result from using different mathematical conventions to describe a given patello-femoral rotation from our analyses differ up to a Root Mean Squared Error of 111.49° for patellar flexion, 55.72° for patellar spin and 35.39° for patellar tilt. To compare clinical kinematic patello-femoral results, every dataset must follow the same convention. Furthermore, researchers should be aware of the used convention's implications to ensure reproducibility when interpreting and comparing such data.

Computational Analysis of Bone Remodeling in the Proximal Tibia Under Electrical Stimulation Considering the Piezoelectric Properties.

The piezoelectricity of bone is known to play a crucial role in bone adaptation and remodeling. The application of an external stimulus such as mechanical strain or electric field has the potential to enhance bone formation and implant osseointegration. Therefore, in the present study, the objective is to investigate bone remodeling under electromechanical stimulation as a step towards establishing therapeutic strategies. For the first time, piezoelectric bone remodeling in the human proximal tibia under electro-mechanical loads was analyzed using the finite element method in an open-source framework. The predicted bone density distributions were qualitatively and quantitatively assessed by comparing with the computed tomography (CT) scan and the bone mineral density (BMD) calculated from the CT, respectively. The effect of model parameters such as uniform initial bone density and reference stimulus on the final density distribution was investigated. Results of the parametric study showed that for different values of initial bone density the model predicted similar but not identical final density distribution. It was also shown that higher reference stimulus value yielded lower average bone density at the final time. The present study demonstrates an increase in bone density as a result of electrical stimulation. Thus, to minimize bone loss, for example, due to physical impairment or osteoporosis, mechanical loads during daily physical activities could be partially replaced by therapeutic electrical stimulation.

Reporting checklist for verification and validation of finite element analysis in orthopedic and trauma biomechanics.

Finite element analysis (FEA) has become a fundamental tool for biomechanical investigations in the last decades. Despite several existing initiatives and guidelines for reporting on research methods and results, there are still numerous issues that arise when using computational models in biomechanical investigations. According to our knowledge, these problems and controversies lie mainly in the verification and validation (V&V) process as well as in the set-up and evaluation of FEA. This work aims to introduce a checklist including a report form defining recommendations for FEA in the field of Orthopedic and Trauma (O&T) biomechanics. Therefore, a checklist was elaborated which summarizes and explains the crucial methodologies for the V&V process. In addition, a report form has been developed which contains the most important steps for reporting future FEA. An example of the report form is shown, and a template is provided, which can be used as a uniform basis for future documentation. The future application of the presented report form will show whether serious errors in biomechanical investigations using FEA can be minimized by this checklist. Finally, the credibility of the FEA in the clinical area and the scientific exchange in the community regarding reproducibility and exchangeability can be improved.

Finite element analysis of bone remodelling with piezoelectric effects using an open-source framework.

Bone tissue exhibits piezoelectric properties and thus is capable of transforming mechanical stress into electrical potential. Piezoelectricity has been shown to play a vital role in bone adaptation and remodelling processes. Therefore, to better understand the interplay between mechanical and electrical stimulation during these processes, strain-adaptive bone remodelling models without and with considering the piezoelectric effect were simulated using the Python-based open-source software framework. To discretise numerical attributes, the finite element method (FEM) was used for the spatial variables and an explicit Euler scheme for the temporal derivatives. The predicted bone apparent density distributions were qualitatively and quantitatively evaluated against the radiographic scan of a human proximal femur and the bone apparent density calculated using a bone mineral density (BMD) calibration phantom, respectively. Additionally, the effect of the initial bone density on the resulting predicted density distribution was investigated globally and locally. The simulation results showed that the electrically stimulated bone surface enhanced bone deposition and these are in good agreement with previous findings from the literature. Moreover, mechanical stimuli due to daily physical activities could be supported by therapeutic electrical stimulation to reduce bone loss in case of physical impairment or osteoporosis. The bone remodelling algorithm implemented using an open-source software framework facilitates easy accessibility and reproducibility of finite element analysis made.

Stress and strain distribution in femoral heads for hip resurfacing arthroplasty with different materials: A finite element analysis.

Femoral bone loss due to stress and strain shielding is a common problem in hip resurfacing arthroplasty (HRA), which arises from the different stiffness of implant materials and the adjacent bone. Usually, the implants used in HRA are made of cobalt-chromium alloy (CoCr). As a novel concept, implants may also be made of ceramics, whose stiffness exceeds that of the adjacent bone by a multiple. Therefore, this computational study aimed to evaluate whether poly (ether-ether-ketone) (PEEK) or a hybrid material with a PEEK body and ceramic surface made of alumina toughened zirconia (ATZ) might be more suitable implant alternatives for HRA, as they can avoid stress and strain shielding. A reconstructed model of a human femur with an HRA implant was simulated, whereby the material of the HRA was varied between CoCr, ATZ, zirconia toughened alumina (ZTA), PEEK, and a hybrid PEEK-ATZ material. The implant fixation method also varied (cemented or cementless). The simulated models were compared with an intact model to analyze stress and strain distribution in the femoral head and neck. The strain distribution was evaluated at a total of 30,344 (cemented HRA) and 63,531 (uncemented HRA) nodes in the femoral head and neck region and divided into different strain regions (<400 µm/m: atrophy; 400-3000 µm/m: bone preserving and building; 3000-20,000 µm/m: yielding and >20,000 µm/m fracture). In addition, the mechanical stability of the implants was evaluated. When the material of the HRA implant was simulated as metal or ceramic while evaluating the strains, it was seen that around 22-26% of the analyzed nodes in the femoral head and neck were in an atrophic region, 47-51% were in a preserving or building region, and 27-28% were in a yielding region. In the case of PEEK implant, less than 0.5% of the analyzed nodes were in an atrophic region, 66-69% in a preserving or building region, and 31-34% in a yielding region. The fixation technique also had a small influence. When a hybrid HRA was simulated, the strains at the analyzed nodes depended on the thickness of the ceramic material. In conclusion, the material of the HRA implant was crucial in terms of stress and strain distribution in the adjacent bone. HRA made of PEEK or a hybrid material leads to decisively reduced stress and strain alteration compared to stiffer materials such as CoCr, ATZ, and ZTA. This confirms the potential for reduction in stress and strain shielding in the femoral head with the use of a hybrid material with a PEEK body for HRA.

Musculoskeletal Multibody Simulation Analysis on the Impact of Patellar Component Design and Positioning on Joint Dynamics after Unconstrained Total Knee Arthroplasty.

Patellofemoral (PF) disorders are considered a major clinical complication after total knee replacement (TKR). Malpositioning and design of the patellar component impacts knee joint dynamics, implant fixation and wear propagation. However, only a limited number of studies have addressed the biomechanical impact of the patellar component on PF dynamics and their results have been discussed controversially. To address these issues, we implemented a musculoskeletal multibody simulation (MMBS) study for the systematical analysis of the patellar component's thickness and positioning on PF contact forces and kinematics during dynamic squat motion with virtually implanted unconstrained cruciate-retaining (CR)-TKR. The patellar button thickness clearly increased the contact forces in the PF joint (up to 27%). Similarly, the PF contact forces were affected by superior-inferior positioning (up to 16%) and mediolateral positioning (up to 8%) of the patellar button. PF kinematics was mostly affected by the mediolateral positioning and the thickness of the patellar component. A medialization of 3 mm caused a lateral patellar shift by up to 2.7 mm and lateral patellar tilt by up to 1.6°. However, deviations in the rotational positioning of the patellar button had minor effects on PF dynamics. Aiming at an optimal intraoperative patellar component alignment, the orthopedic surgeon should pay close attention to the patellar component thickness in combination with its mediolateral and superior-inferior positioning on the retropatellar surface. Our generated MMBS model provides systematic and reproducible insight into the effects of patellar component positioning and design on PF dynamics and has the potential to serve as a preoperative analysis tool.

A round-robin finite element analysis of human femur mechanics between seven participating laboratories with experimental validation.

Finite element analysis is a common tool that has been used for the past few decades to predict the mechanical behavior of bone. However, to our knowledge, there are no round-robin finite element analyses of long human bones with more than two participating biomechanics laboratories published yet, where the results of the experimental tests were not known in advance. We prepared a fresh-frozen human femur for a compression test in a universal testing machine measuring the strains at 10 bone locations as well as the deformation of the bone in terms of the displacement of the loading point at a load of 2 kN. The computed tomography data of the bone with a calibration phantom as well as the orientation of the bone in the testing machine with the according boundary conditions were delivered to seven participating laboratories. These were asked to perform a finite element analysis simulating the experimental setup and deliver their results to the coordinator without knowing the experimental results. Resultantly, four laboratories had deviations from the experimentally measured strains of less than 40%, and three laboratories had deviations of their numerically determined values compared to the experimental data of more than 120%. These deviations are thought to be based on different material laws and material data, as well as different material mapping methods. Investigations will be conducted to clarify and assess the reasons for the large deviations in the numerical data. It was shown that the precision of finite element models of the human femur is not yet as developed as desired by the biomechanics community.

Neuro-musculoskeletal flexible multibody simulation yields a framework for efficient bone failure risk assessment.

Fragility fractures are a major socioeconomic problem. A non-invasive, computationally-efficient method for the identification of fracture risk scenarios under the representation of neuro-musculoskeletal dynamics does not exist. We introduce a computational workflow that integrates modally-reduced, quantitative CT-based finite-element models into neuro-musculoskeletal flexible multibody simulation (NfMBS) for early bone fracture risk assessment. Our workflow quantifies the bone strength via the osteogenic stresses and strains that arise due to the physiological-like loading of the bone under the representation of patient-specific neuro-musculoskeletal dynamics. This allows for non-invasive, computationally-efficient dynamic analysis over the enormous parameter space of fracture risk scenarios, while requiring only sparse clinical data. Experimental validation on a fresh human femur specimen together with femur strength computations that were consistent with literature findings provide confidence in the workflow: The simulation of an entire squat took only 38 s CPU-time. Owing to the loss (16% cortical, 33% trabecular) of bone mineral density (BMD), the strain measure that is associated with bone fracture increased by 31.4%; and yielded an elevated risk of a femoral hip fracture. Our novel workflow could offer clinicians with decision-making guidance by enabling the first combined in-silico analysis tool using NfMBS and BMD measurements for optimized bone fracture risk assessment.