Protocol for constructing subject-specific biomechanical models of knee joint

A robust protocol for building subject-specific biomechanical models of the human knee joint is proposed which uses magnetic resonance imaging, motion analysis and force platform data in conjunction with detailed 3D finite element models. The proposed protocol can be used for determining stress and strain distributions and contact kinetics in different knee elements at different body postures during various physical activities. Several examples are provided to highlight the capabilities and potential applications of the proposed protocol. This includes preliminary results on the role of body weight on the stresses and strains induced in the knee articular cartilages and meniscus during single-leg stance and calculations of the induced stresses and ligament forces during the gait cycle.     

A biomechanical model to assess the contribution of pelvic musculature weakness to the development of stress urinary incontinence

A biomechanical model of the female pelvic support system was developed to explore the contribution of pelvic floor muscle defect to the development of stress urinary incontinence (SUI). From a pool of 135 patients, clinical data of 26 patients with pelvic muscular defect were used in modelling. The model was employed to estimate the parameters that describe the stiffness properties of the vaginal wall and ligament tissues for individual patients. The parameters were then implemented into the model to evaluate for each patient the impact of pelvic muscular defect on the vaginal apex support and the bladder neck support, a factor that relates to the onset of SUI. For the modelling analysis, the compromise of pelvic muscular support was demonstrated to contribute to vaginal apex prolapse and bladder neck prolapse, a condition commonly seen in SUI patients, while simulated conditions of restored muscular support were shown to help re-establish both vaginal apex and bladder neck supports. The findings illustrate the significance of pelvic muscle strength to vaginal support and urinary continence; therefore, the clinical recommendation of pelvic muscle strengthening, such as Kegel exercises, has been shown to be an effective treatment for patients with SUI symptoms.     

Image-guided simulation of tissue deformation using a mechanical model on a surgical application

This paper presents a method for localizing the position of a liver and a tumor within the tissue during a minimally invasive liver operation. From pre-operative CT scans, the liver volume and its internal structures are segmented using image-processing techniques. Based on these segmentations, a three-dimensional mechanical model is built to compute the liver volume and internal structure displacement under boundary conditions such as external forces from the surgical instrument. This can help the surgeon understand the motion of internal structures when manipulating the liver. To validate our method, an experiment on a porcine liver explant was performed to assess the difference between actual tissue motion and the mechanical model.     

Influence of connection type on the biomechanical behavior of distal extension mandibular removable partial dentures supported by implants and natural teeth

Few studies are performed to evaluate the influence of connection type on the stress distribution of distal extension mandibular removable partial dentures (RPDs) supported by both implants and natural teeth. In this study, five three-dimensional finite element models were prepared to simulate mandibular bilateral partially edentulous arches. Four were RPDs supported by both implants and natural teeth, and the other one was RPDs supported only by natural teeth. The maximum equivalent (EQV) stress values of bone around implants, the abutments, and the mucosa displacements of the related supporting structures were measured. It was found that a non-rigid telescopic coping was more favorable to protect the implant than a rigid telescopic coping. Compared with other connection types, the easy resilient attachment (ERA) system seemed to be effective to associate implant without complications. However, the results obtained in the present study should be cautiously interpreted in the clinic.     

The effect of hydrogel injection on cardiac function and myocardial mechanics in a computational post-infarction model

An emerging therapy to limit adverse heart remodelling following myocardial infarction (MI) is the injection of polymers into the infarcted left ventricle (LV). In the few numerical studies carried out in this field, the definition and distribution of the hydrogel in the infarcted myocardium were simplified. In this computational study, a more realistic biomaterial distribution was simulated after which the effect on cardiac function and mechanics was studied. A validated finite element heart model was used in which an antero-apical infarct was defined. Four infarct models were created representing different temporal phases in the progression of a MI. Hydrogel layers were simulated in the infarcted myocardium in each model. Biomechanical and functional improvement of the LV was found after hydrogel inclusion in the ischaemic models representing the early phases of MI. In contrast, only functional but no mechanical restitution was shown in the scar model due to hydrogel presence.     

Patient-specific finite element model of the spine and spinal cord to assess the neurological impact of scoliosis correction: preliminary application on two cases with and without intraoperative neurological complications

Scoliosis is a 3D deformation of the spine and rib cage. For severe cases, surgery with spine instrumentation is required to restore a balanced spine curvature. This surgical procedure may represent a neurological risk for the patient, especially during corrective maneuvers. This study aimed to computationally simulate the surgical instrumentation maneuvers on a patient-specific biomechanical model of the spine and spinal cord to assess and predict potential damage to the spinal cord and spinal nerves. A detailed finite element model (FEM) of the spine and spinal cord of a healthy subject was used as reference geometry. The FEM was personalized to the geometry of the patient using a 3D biplanar radiographic reconstruction technique and 3D dual kriging. Step by step surgical instrumentation maneuvers were simulated in order to assess the neurological risk associated to each maneuver. The surgical simulation methodology implemented was divided into two parts. First, a global multi-body simulation was used to extract the 3D displacement of six vertebral landmarks, which were then introduced as boundary conditions into the personalized FEM in order to reproduce the surgical procedure. The results of the FEM simulation for two cases were compared to published values on spinal cord neurological functional threshold. The efficiency of the reported method was checked considering one patient with neurological complications detected during surgery and one control patient. This comparison study showed that the patient-specific hybrid model reproduced successfully the biomechanics of neurological injury during scoliosis correction maneuvers.     

Influence of custom-made implant designs on the biomechanical performance for the case of immediate post-extraction placement in the maxillary esthetic zone: a finite element analysis

Due to the increasing adoption of immediate implantation strategies and the rapid development of the computer aided design/computer aided manufacturing technology, a therapeutic concept based on patient-specific implant dentistry has recently been reintroduced by many researchers. However, little information is available on the designs of custom-made dental implant systems, especially their biomechanical behavior. The influence of the custom-made implant designs on the biomechanical performance for both an immediate and a delayed loading protocol in the maxillary esthetic zone was evaluated by means of the finite element (FE) method. FE models of three dental implants were considered: a state of the art cylindrical implant and two custom-made implants designed by reverse engineering technology, namely a root-analogue implant and a root-analogue threaded implant. The von Mises stress distributions and micro-motions around the bone-implant interfaces were calculated using ANSYS software. In a comparison of the three implant designs for both loading protocols, a favorable biomechanical performance was observed for the use of root-analogue threaded implant which approximated the geometry of natural anterior tooth and maintained the original long-axis. The results indicated that bone-implant interfacial micro-motion was reduced and a favorable stress distribution after osseointegration was achieved.     

Biomechanical effects of posterior pedicle fixation techniques on the adjacent segment for the treatment of thoracolumbar burst fractures: a biomechanical analysis

Abstract

Posterior pedicle fixation technique is a common method for treating thoracolumbar burst fractures, but the effect of different fixation techniques on the postoperative spinal mechanical properties has not been clearly defined, especially on adjacent segments. A finite element model of T10-L2 with moderate T12 vertebra burst fracture was constructed to investigate biomechanical behavior of three posterior pedicle screw fixation techniques. Compared with traditional short-segment 4 pedicle screw fixation (TS-4) and intermediate long-segment 6 pedicle screw fixation (IL-6), mono-segment 4 pedicle screw fixation (MS-4) provides a safer surgical selection to prevent the secondary degeneration of adjacent segments in the long-term.     

Numerical analysis of the influence of nucleus pulposus removal on the biomechanical behavior of a lumbar motion segment

Nucleus replacement was deemed to have therapeutic potential for patients with intervertebral disc herniation. However, whether a patient would benefit from nucleus replacement is technically unclear. This study aimed to investigate the influence of nucleus pulposus (NP) removal on the biomechanical behavior of a lumbar motion segment and to further explore a computational method of biomechanical characteristics of NP removal, which can evaluate the mechanical stability of pulposus replacement. We, respectively, reconstructed three types of models for a mildly herniated disc and three types of models for a severely herniated disc based on a L4–L5 segment finite element model with computed tomography image data from a healthy adult. First, the NP was removed from the herniated disc models, and the biomechanical behavior of NP removal was simulated. Second, the NP cavities were filled with an experimental material (Poisson's ratio = 0.3; elastic modulus = 3 MPa), and the biomechanical behavior of pulposus replacement was simulated. The simulations were carried out under the five loadings of axial compression, flexion, lateral bending, extension, and axial rotation. The changes of the four biomechanical characteristics, i.e. the rotation degree, the maximum stress in the annulus fibrosus (AF), joint facet contact forces, and the maximum disc deformation, were computed for all models. Experimental results showed that the rotation range, the maximum AF stress, and joint facet contact forces increased, and the maximum disc deformation decreased after NP removal, while they changed in the opposite way after the nucleus cavities were filled with the experimental material.     

Influence of anisotropic bone properties on the biomechanical behavior of the acetabular cup implant: a multiscale finite element study

Although the biomechanical behavior of the acetabular cup (AC) implant is determinant for the surgical success, it remains difficult to be assessed due to the multiscale and anisotropic nature of bone tissue. The aim of the present study was to investigate the influence of the anisotropic properties of peri-implant trabecular bone tissue on the biomechanical behavior of the AC implant at the macroscopic scale. Thirteen bovine trabecular bone samples were imaged using micro-computed tomography (μCT) with a resolution of 18 μm. The anisotropic biomechanical properties of each sample were determined at the scale of the centimeter based on a dedicated method using asymptotic homogenization. The material properties obtained with this multiscale approach were used as input data in a 3D finite element model to simulate the macroscopic mechanical behavior of the AC implant under different loading conditions. The largest stress and strain magnitudes were found around the equatorial rim and in the polar area of the AC implant. All macroscopic stiffness quantities were significantly correlated (R² > 0.85, p < 6.5 e-6) with BV/TV (bone volume/total volume). Moreover, the maximum value of the von Mises stress field was significantly correlated with BV/TV (R² > 0.61, p < 1.6 e-3) and was always found at the bone-implant interface. However, the mean value of the microscopic stress (at the scale of the trabeculae) decrease as a function of BV/TV for vertical and torsional loading and do not depend on BV/TV for horizontal loading. These results highlight the importance of the anisotropic properties of bone tissue.     

The potential role of variations in juvenile hip geometry on the development of Legg-Calvé-Perthes disease: a biomechanical investigation

Legg-Calvé-Perthes disease (LCP) is one of the most poorly understood diseases in paediatric orthopaedics. One common trait of LCP is the marked morphological difference between healthy and pathological hips, early deviations of which (i.e. prior to disease onset) have been suggested to lead to the overload and collapse of the epiphysis. Here, the impact of common variations in geometry is investigated with a finite element model of a juvenile femur under single leg standing and landing. Here, the impact of typical variations in geometry is investigated with a finite element model of a juvenile femur under single leg standing and landing. The variations appear to have only a limited effect on the stress distribution in the femoral epiphysis even during high impact activities. This suggests that, for this individual at least, they would be unlikely to cause epiphyseal overload and collapse, even in the presence of a skeletally immature epiphysis.     

Development and validation of a computational model of the knee joint for the evaluation of surgical treatments for osteoarthritis

A three-dimensional (3D) knee joint computational model was developed and validated to predict knee joint contact forces and pressures for different degrees of malalignment. A 3D computational knee model was created from high-resolution radiological images to emulate passive sagittal rotation (full-extension to 65°-flexion) and weight acceptance. A cadaveric knee mounted on a six-degree-of-freedom robot was subjected to matching boundary and loading conditions. A ligament-tuning process minimised kinematic differences between the robotically loaded cadaver specimen and the finite element (FE) model. The model was validated by measured intra-articular force and pressure measurements. Percent full scale error between FE-predicted and in vitro-measured values in the medial and lateral compartments were 6.67% and 5.94%, respectively, for normalised peak pressure values, and 7.56% and 4.48%, respectively, for normalised force values. The knee model can accurately predict normalised intra-articular pressure and forces for different loading conditions and could be further developed for subject-specific surgical planning.     

Effects of sitting postures on risks for deep tissue injury in the residuum of a transtibial prosthetic-user: a biomechanical case study

Transtibial amputation prosthetic-users are at risk of developing deep tissue injury (DTI) while donning their prosthesis for prolonged periods; however, no study addresses the mechanical loading of the residuum during sitting with a prosthesis. We combined MRI-based 3D finite element modelling of a residuum with an injury threshold and a muscle damage law to study risks for DTI in one sitting subject in two postures: 30°-knee-flexion vs. 90°-knee-flexion. We recorded skin-socket pressures, used as model boundary conditions. During the 90°-knee-flexion simulations, major internal muscle injuries were predicted (>1000 mm³). In contrast, the 30°-knee-flexion simulations only produced minor injury ( < 14 mm³). Predicted injury rates at 90°-knee-flexion were over one order of magnitude higher than those at 30°-knee-flexion. We concluded that in this particular subject, prolonged 90°-knee-flexion sitting theoretically endangers muscle viability in the residuum. By expanding the studies to large subject groups, this research approach can support development of guidelines for DTI prevention in prosthetic-users.     

Patient-specific simulation of a stentless aortic valve implant: the impact of fibres on leaflet performance

In some cases of aortic valve leaflet disease, the implant of a stentless biological prosthesis represents an excellent option for aortic valve replacement (AVR). In particular, if compared with the implant of mechanical valves, it provides a more physiological haemodynamic performance and a reduced thrombogeneticity, avoiding the use of anticoagulants. The clinical outcomes of AVR are strongly dependent on an appropriate choice of both prosthesis size and replacement technique, which is, at present, strictly related to surgeon's experience and skill. This represents the motivation for patient-specific finite element analysis able to virtually reproduce stentless valve implantation. With the aim of performing reliable patient-specific simulations, we remark that, on the one hand, it is not well established in the literature whether bioprosthetic leaflet tissue is isotropic or anisotropic; on the other hand, it is of fundamental importance to incorporate an accurate material model to realistically predict post-operative performance. Within this framework, using a novel computational methodology to simulate stentless valve implantation, we test the impact of using different material models on both the stress pattern and post-operative coaptation parameters (i.e. coaptation area, length and height). As expected, the simulation results suggest that the material properties of the valve leaflets affect significantly the post-operative prosthesis performance.     

Development of a finite element musculoskeletal model with the ability to predict contractions of three-dimensional muscles

Representation of realistic muscle geometries is needed for systematic biomechanical simulation of musculoskeletal systems. Most of the previous musculoskeletal models are based on multibody dynamics simulation with muscles simplified as one-dimensional (1D) line-segments without accounting for the large muscle attachment areas, spatial fibre alignment within muscles and contact and wrapping between muscles and surrounding tissues. In previous musculoskeletal models with three-dimensional (3D) muscles, contractions of muscles were among the inputs rather than calculated, which hampers the predictive capability of these models. To address these issues, a finite element musculoskeletal model with the ability to predict contractions of 3D muscles was developed. Muscles with realistic 3D geometry, spatial muscle fibre alignment and muscle-muscle and muscle-bone interactions were accounted for. Active contractile stresses of the 3D muscles were determined through an efficient optimization approach based on the measured kinematics of the lower extremity and ground force during gait. This model also provided stresses and strains of muscles and contact mechanics of the muscle-muscle and muscle-bone interactions. The total contact force of the knee predicted by the model corresponded well to the in vivo measurement. Contact and wrapping between muscles and surrounding tissues were evident, demonstrating the need to consider 3D contact models of muscles. This modelling framework serves as the methodological basis for developing musculoskeletal modelling systems in finite element method incorporating 3D deformable contact models of muscles, joints, ligaments and bones.     

Validation of a full body finite element model (THUMS) for running-type impacts to the lower extremity

The purpose of this study was to determine whether modifying an existing, highly biofidelic full body finite element model [total human model for safety (THUMS)] would produce valid amplitude and temporal shock wave characteristics as it travels proximally through the lower extremity. Modifying an existing model may be more feasible than developing a new model, in terms of cost, labour and expertise. The THUMS shoe was modified to more closely simulate the material properties of a heel pad. Relative errors in force and acceleration data from experimental human pendulum impacts and simulated THUMS impacts were 22% and 54%, respectively, across the time history studied. The THUMS peak acceleration was attenuated by 57.5% and took 19.7 ms to travel proximally along the lower extremity. Although refinements may be necessary to improve force and acceleration timing, the modified THUMS represented, to a certain extent, shock wave propagation and attenuation demonstrated by living humans under controlled impact conditions.     

A biomechanical model to simulate the effect of a high vertical loading on trunk flexural stiffness

A human trunk model was developed to simulate the effect of a high vertical loading on trunk flexural stiffness. A force–length relationship is attributed to each muscle of the multi-body model. Trunk stiffness and muscle forces were evaluated experimentally and numerically for various applied loads. Experimental evaluation of trunk stiffness was carried out by measuring changes in reaction force following a sudden horizontal displacement at the T10 level prior to paraspinal reflexes induction. Results showed that the trunk stiffness increases under small applied loads, peaks when the loads were further increased and decreases when higher loads are applied. A sensitivity analysis to muscle force–length relationship is provided to determine the model's limitations. This model pointed out the importance of taking into account the changes in muscle length to evaluate the effect of spinal loads beyond the safe limit that cannot be evaluated experimentally and to predict the trunk instability under vertical load.     

The effect of meniscal tears and resultant partial meniscectomies on the knee contact stresses: a finite element analysis

Knee osteoarthritis (OA) is believed to result from high levels of contact stresses on the articular cartilage and meniscus after meniscal damage. This study investigated the effect of meniscal tears and partial meniscectomies on the peak compressive and shear stresses in the human knee joint. An elaborate three-dimensional finite element model of knee joint including bones, articular cartilages, menisci and main ligaments was developed from computed tomography and magnetic resonance imaging images. This model was used to model four types of meniscal tears and their resultant partial meniscectomies and analysed under an axial 1150 N load at 0° flexion. Three different conditions were compared: a healthy knee joint, a knee joint with medial meniscal tears and a knee joint following partial meniscectomies. The numerical results showed that each meniscal tear and its resultant partial meniscectomy led to an increase in the peak compressive and shear stresses on the articular cartilages and meniscus in the medial knee compartment, especially for partial meniscectomy. Among the four types of meniscal tears, the oblique tear resulted in the highest values of the peak compressive and shear stresses. For the four partial meniscectomies, longitudinal meniscectomy led to the largest increase in these two stresses. The lateral compartment was minimally affected by all the simulations. The results of this study demonstrate meniscal tear and its resultant partial meniscectomy has a positive impact on the maintenance of high levels of contact stresses, which may improve the progression of knee OA, especially for partial meniscectomy. Surgeons should adopt a prudent strategy to preserve the greatest amount of meniscus possible.     

A partial linear model in the outcome-dependent sampling setting to evaluate the effect of prenatal PCB exposure on cognitive function in children

Outcome-dependent sampling (ODS) has been widely used in biomedical studies because it is a cost-effective way to improve study efficiency. However, in the setting of a continuous outcome, the representation of the exposure variable has been limited to the framework of linear models, due to the challenge in terms of both theory and computation. Partial linear models (PLM) are a powerful inference tool to nonparametrically model the relation between an outcome and the exposure variable. In this article, we consider a case study of a PLM for data from an ODS design. We propose a semiparametric maximum likelihood method to make inferences with a PLM. We develop the asymptotic properties and conduct simulation studies to show that the proposed ODS estimator can produce a more efficient estimate than that from a traditional simple random sampling design with the same sample size. Using this newly developed method, we were able to explore an open question in epidemiology: whether in utero exposure to background levels of polychlorinated biphenyls (PCBs) is associated with children's intellectual impairment. Our model provides further insights into the relation between low-level PCB exposure and children's cognitive function. The results shed new light on a body of inconsistent epidemiologic findings.