Long-term implant fixation and stress-shielding in total hip replacement

Implant fixation implies a strong and durable mechanical bond between the prosthetic component and host skeleton. Assuming the short-term impediments to implant fixation are successfully addressed and that longer-term issues such as late infection and mechanical failure of the components are avoided, the biological response of the host tissue to the presence of the implant is critical to long-term success. In particular, maintenance of adequate peri-prosthetic bone stock is a key factor. Two major causes of bone loss in the supporting bone are adverse bone remodeling in response to debris shed from the implant and stress-shielding. Here, I review some of the major lessons learned from studying stress-shielding-induced bone loss. It is well known that stress-shielding can be manipulated by altering implant design, but less well appreciated that the development of bone anabolic agents may make it possible to reduce the severity of stress-shielding and the associated bone loss by augmenting the host skeleton through the use of locally or systemically delivered agents. In most cases, mechanical, material and biological factors do not act in isolation, emphasizing that it is often not possible to optimize all boundary conditions.

Comparative modeling of transcranial magnetic and electric stimulation in mouse, monkey, and human

Transcranial magnetic stimulation (TMS) and transcranial electric stimulation (TES) are increasingly popular methods to noninvasively affect brain activity. However, their mechanism of action and dose-response characteristics remain under active investigation. Translational studies in animals play a pivotal role in these efforts due to a larger neuroscientific toolset enabled by invasive recordings. In order to translate knowledge gained in animal studies to humans, it is crucial to generate comparable stimulation conditions with respect to the induced electric field in the brain. Here, we conduct a finite element method (FEM) modeling study of TMS and TES electric fields in a mouse, capuchin and macaque monkeys, and a human model. We systematically evaluate the induced electric fields and analyze their relationship to head and brain anatomy. We find that with increasing head size, TMS-induced electric field strength first increases and then decreases according to a two-term exponential function. TES-induced electric field strength strongly decreases from smaller to larger specimen with up to 100x fold differences across species. Our results can serve as a basis to compare and match stimulation parameters across studies in animals and humans.

Development of a Metric for Predicting Brain Strain Responses Using Head Kinematics

Diffuse brain injuries are caused by excessive brain deformation generated primarily by rapid rotational head motion. Metrics that describe the severity of brain injury based on head motion often do not represent the governing physics of brain deformation, rendering them ineffective over a broad range of head impact conditions. This study develops a brain injury metric based on the response of a second-order mechanical system, and relates rotational head kinematics to strain-based brain injury metrics: maximum principal strain (MPS) and cumulative strain damage measure (CSDM). This new metric, universal brain injury criterion (UBrIC), is applicable over a broad range of kinematics encountered in automotive crash and sports. Efficacy of UBrIC was demonstrated by comparing it to MPS and CSDM predicted in 1600 head impacts using two different finite element (FE) brain models. Relative to existing metrics, UBrIC had the highest correlation with the FE models, and performed better in most impact conditions. While UBrIC provides a reliable measurement for brain injury assessment in a broad range of head impact conditions, and can inform helmet and countermeasure design, an injury risk function was not incorporated into its current formulation until validated strain-based risk functions can be developed and verified against human injury data.

Skeletal maturity leads to a reduction in the strain magnitudes induced within the bone: a murine tibia study

Bone adapts to changes in the local mechanical environment (e.g. strains) through formation and resorption processes. However, the bone adaptation response is significantly reduced with increasing age. The mechanical strains induced within the bone by external loading are determined by bone morphology and tissue material properties. Although it is known that changes in bone mass, architecture and bone tissue quality occur with age, to what extent they contribute to the altered bone adaptation response remains to be determined. This study investigated alterations in strains induced in the tibia of different aged female C57Bl/6J mice (young, 10-week-old; adult, 26-week-old; and elderly, 78-week-old) subjected to in vivo compressive loading. Using a combined in vivo/in silico approach, the strains in the bones were assessed by both strain gauging and finite element modeling experiments. In cortical bone, strain magnitudes induced at the mid-diaphysis decreased by 20% from young to adult mice and by 15% from adult to elderly mice. In the cancellous bone (at the proximal metaphysis), induced strains were 70% higher in young compared with adult and elderly mice. Taking into account previous studies showing a reduced bone adaptation response to mechanical loading in adulthood, these results suggest that the diminished adaptive response is in part due to a reduction in the strains induced within the bone.

The Shrinking Brain: Cerebral Atrophy Following Traumatic Brain Injury

Cerebral atrophy in response to traumatic brain injury is a well-documented phenomenon in both primary investigations and review articles. Recent atrophy studies focus on exploring the region-specific patterns of cerebral atrophy; yet, there is no study that analyzes and synthesizes the emerging atrophy patterns in a single comprehensive review. Here we attempt to fill this gap in our current knowledge by integrating the current literature into a cohesive theory of preferential brain tissue loss and by identifying common risk factors for accelerated atrophy progression. Our review reveals that observations for mild traumatic brain injury remain inconclusive, whereas observations for moderate-to-severe traumatic brain injury converge towards robust patterns: brain tissue loss is on the order of 5% per year, and occurs in the form of generalized atrophy, across the entire brain, or focal atrophy, in specific brain regions. The most common regions of focal atrophy are the thalamus, hippocampus, and cerebellum in gray matter and the corpus callosum, corona radiata, and brainstem in white matter. We illustrate the differences of generalized and focal gray and white matter atrophy on emerging deformation and stress profiles across the whole brain using computational simulation. The characteristic features of our atrophy simulations-a widening of the cortical sulci, a gradual enlargement of the ventricles, and a pronounced cortical thinning-agree well with clinical observations. Understanding region-specific atrophy patterns in response to traumatic brain injury has significant implications in modeling, simulating, and predicting injury outcomes. Computational modeling of brain atrophy could open new strategies for physicians to make informed decisions for whom, how, and when to administer pharmaceutical treatment to manage the chronic loss of brain structure and function.

Resistive exercise in astronauts on prolonged spaceflights provides partial protection against spaceflight-induced bone loss

Bone loss in astronauts during spaceflight may be a risk factor for osteoporosis, fractures and renal stone formation. We previously reported that the bisphosphonate alendronate, combined with exercise that included an Advanced Resistive Exercise Device (ARED), can prevent or attenuate group mean declines in areal bone mineral density (aBMD) measured soon after ~ 6-month spaceflights aboard the International Space Station (ISS). It is unclear however if the beneficial effects on postflight aBMD were due to individual or combined effects of alendronate and ARED. Hence, 10 additional ISS astronauts were recruited who used the ARED (ARED group) without drug administration using similar measurements in the previous study, i.e., densitometry, biochemical assays and analysis of finite element (FE) models. In addition densitometry data (DXA and QCT only) were compared to published data from crewmembers (n = 14-18) flown prior to in-flight access to the ARED (Pre-ARED). Group mean changes from preflight (± SD %) were used to evaluate effects of countermeasures as sequentially modified on the ISS (i.e., Pre-ARED vs. ARED; ARED vs. Bis+ARED). Spaceflight durations were not significantly different between groups. Postflight bone density measurements were significantly reduced from preflight in the Pre-ARED group. As previously reported, combined Bis+ARED prevented declines in all DXA and QCT hip densitometry and in estimates of FE hip strengths; increased the aBMD of lumbar spine; and prevented elevations in urinary markers for bone resorption during spaceflight. ARED without alendronate partially attenuated declines in bone mass but did not suppress biomarkers for bone resorption or prevent trabecular bone loss. Resistive exercise in the ARED group did not prevent declines in hip trabecular vBMD, but prevented reductions in cortical vBMD of the femoral neck, in FE estimate of hip strength for non-linear stance (NLS) and in aBMD of the femoral neck. We conclude that a bisphosphonate, when combined with resistive exercise, enhances the preservation of bone mass because of the added suppression of bone resorption in trabecular bone compartment not evident with ARED alone.

Ventricular wall biomaterial injection therapy after myocardial infarction: Advances in material design, mechanistic insight and early clinical experiences

Intramyocardial biomaterial injection therapy for myocardial infarction has made significant progress since concept initiation more than 10 years ago. The interim successes and progress in the first 5 years have been extensively reviewed. During the last 5 years, two phase II clinical trials have reported their long term follow up results and many additional biomaterial candidates have reached preclinical and clinical testing. Also in recent years deeper investigations into the mechanisms behind the beneficial effects associated with biomaterial injection therapy have been pursued, and a variety of process and material parameters have been evaluated for their impact on therapeutic outcomes. This review explores the advances made in this biomaterial-centered approach to ischemic cardiomyopathy and discusses potential future research directions as this therapy seeks to positively impact patients suffering from one of the world's most common sources of mortality.

Finite element modeling of impulsive excitation and shear wave propagation in an incompressible, transversely isotropic medium

Elastic properties of materials can be measured by observing shear wave propagation following localized, impulsive excitations and relating the propagation velocity to a model of the material. However, characterization of anisotropic materials is difficult because of the number of elasticity constants in the material model and the complex dependence of propagation velocity relative to the excitation axis, material symmetries, and propagation directions. In this study, we develop a model of wave propagation following impulsive excitation in an incompressible, transversely isotropic (TI) material such as muscle. Wave motion is described in terms of three propagation modes identified by their polarization relative to the material symmetry axis and propagation direction. Phase velocities for these propagation modes are expressed in terms of five elasticity constants needed to describe a general TI material, and also in terms of three constants after the application of two constraints that hold in the limit of an incompressible material. Group propagation velocities are derived from the phase velocities to describe the propagation of wave packets away from the excitation region following localized excitation. The theoretical model is compared to the results of finite element (FE) simulations performed using a nearly incompressible material model with the five elasticity constants chosen to preserve the essential properties of the material in the incompressible limit. Propagation velocities calculated from the FE displacement data show complex structure that agrees quantitatively with the theoretical model and demonstrates the possibility of measuring all three elasticity constants needed to characterize an incompressible, TI material.

Machine learning and individual variability in electric field characteristics predict tDCS treatment response

BACKGROUND

Transcranial direct current stimulation (tDCS) is widely investigated as a therapeutic tool to enhance cognitive function in older adults with and without neurodegenerative disease. Prior research demonstrates that electric current delivery to the brain can vary significantly across individuals. Quantification of this variability could enable person-specific optimization of tDCS outcomes. This pilot study used machine learning and MRI-derived electric field models to predict working memory improvements as a proof of concept for precision cognitive intervention.

METHODS

Fourteen healthy older adults received 20 minutes of 2 mA tDCS stimulation (F3/F4) during a two-week cognitive training intervention. Participants performed an N-back working memory task pre-/post-intervention. MRI-derived current models were passed through a linear Support Vector Machine (SVM) learning algorithm to characterize crucial tDCS current components (intensity and direction) that induced working memory improvements in tDCS responders versus non-responders.

MAIN RESULTS

SVM models of tDCS current components had 86% overall accuracy in classifying treatment responders vs. non-responders, with current intensity producing the best overall model differentiating changes in working memory performance. Median current intensity and direction in brain regions near the electrodes were positively related to intervention responses (r=0.811,p<0.001 and r=0.774,p=0.001).

CONCLUSIONS

This study provides the first evidence that pattern recognition analyses of MRI-derived tDCS current models can provide individual prognostic classification of tDCS treatment response with 86% accuracy. Individual differences in current intensity and direction play important roles in determining treatment response to tDCS. These findings provide important insights into mechanisms of tDCS response as well as proof of concept for future precision dosing models of tDCS intervention.

Image-based computational assessment of vascular wall mechanics and hemodynamics in pulmonary arterial hypertension patients

Pulmonary arterial hypertension (PAH) is a disease characterized by an elevated pulmonary arterial (PA) pressure. While several computational hemodynamic models of the pulmonary vasculature have been developed to understand PAH, they are lacking in some aspects, such as the vessel wall deformation and its lack of calibration against measurements in humans. Here, we describe a computational modeling framework that addresses these limitations. Specifically, computational models describing the coupling of hemodynamics and vessel wall mechanics in the pulmonary vasculature of a PAH patient and a normal subject were developed. Model parameters, consisting of linearized stiffness E of the large vessels and Windkessel parameters for each outflow branch, were calibrated against in vivo measurements of pressure, flow and vessel wall deformation obtained, respectively, from right-heart catheterization, phase-contrast and cine magnetic resonance images. Calibrated stiffness E of the proximal PA was 2.0 and 0.5 MPa for the PAH and normal models, respectively. Calibrated total compliance C_T and resistance R_T of the distal vessels were, respectively, 0.32 ml/mmHg and 11.3 mmHg∗min/l for the PAH model, and 2.93 ml/mmHg and 2.6 mmHg∗min/l for the normal model. These results were consistent with previous findings that the pulmonary vasculature is stiffer with more constricted distal vessels in PAH patients. Individual effects on PA pressure due to remodeling of the distal and proximal compartments of the pulmonary vasculature were also investigated in a sensitivity analysis. The analysis suggests that the remodeling of distal vasculature contributes more to the increase in PA pressure than the remodeling of proximal vasculature.

An integrated electromechanical-growth heart model for simulating cardiac therapies

An emerging class of models has been developed in recent years to predict cardiac growth and remodeling (G&R). We recently developed a cardiac G&R constitutive model that predicts remodeling in response to elevated hemodynamics loading, and a subsequent reversal of the remodeling process when the loading is reduced. Here, we describe the integration of this G&R model to an existing strongly coupled electromechanical model of the heart. A separation of timescale between growth deformation and elastic deformation was invoked in this integrated electromechanical-growth heart model. To test our model, we applied the G&R scheme to simulate the effects of myocardial infarction in a realistic left ventricular (LV) geometry using the finite element method. We also simulate the effects of a novel therapy that is based on alteration of the infarct mechanical properties. We show that our proposed model is able to predict key features that are consistent with experiments. Specifically, we show that the presence of a non-contractile infarct leads to a dilation of the left ventricle that results in a rightward shift of the pressure volume loop. Our model also predicts that G&R is attenuated by a reduction in LV dilation when the infarct stiffness is increased.

Effects of collagen microstructure and material properties on the deformation of the neural tissues of the lamina cribrosa

It is widely considered that intraocular pressure (IOP)-induced deformation within the neural tissue pores of the lamina cribrosa (LC) contributes to neurodegeneration and glaucoma. Our goal was to study how the LC microstructure and mechanical properties determine the mechanical insult to the neural tissues within the pores of the LC. Polarized light microscopy was used to measure the collagen density and orientation in histology sections of three sheep optic nerve heads (ONH) at both mesoscale (4.4μm) and microscale (0.73μm) resolutions. Mesoscale fiber-aware FE models were first used to calculate ONH deformations at an IOP of 30mmHg. The results were then used as boundary conditions for microscale models of LC regions. Models predicted large insult to the LC neural tissues, with 95th percentile 1st principal strains ranging from 7 to 12%. Pores near the scleral boundary suffered significantly higher stretch compared to pores in more central regions (10.0±1.4% vs. 7.2±0.4%; p=0.014; mean±SD). Variations in material properties altered the minimum, median, and maximum levels of neural tissue insult but largely did not alter the patterns of pore-to-pore variation, suggesting these patterns are determined by the underlying structure and geometry of the LC beams and pores. To the best of our knowledge, this is the first computational model that reproduces the highly heterogeneous neural tissue strain fields observed experimentally.

STATEMENT OF SIGNIFICANCE

The loss of visual function associated with glaucoma has been attributed to sustained mechanical insult to the neural tissues of the lamina cribrosa due to elevated intraocular pressure. Our study is the first computational model built from specimen-specific tissue microstructure to consider the mechanics of the neural tissues of the lamina separately from the connective tissue. We found that the deformation of the neural tissue was much larger than that predicted by any recent microstructure-aware models of the lamina. These results are consistent with recent experimental data and the highest deformations were found in the region of the lamina where glaucomatous damage first occurs. This study provides new insight into the complex biomechanical environment within the lamina.

Micromechanical poroelastic finite element and shear-lag models of tendon predict large strain dependent Poisson's ratios and fluid expulsion under tensile loading

As tendons are loaded, they reduce in volume and exude fluid to the surrounding medium. Experimental studies have shown that tendon stretching results in a Poisson's ratio greater than 0.5, with a maximum value at small strains followed by a nonlinear decay. Here we present a computational model that attributes this macroscopic observation to the microscopic mechanism of the load transfer between fibrils under stretch. We develop a finite element model based on the mechanical role of the interfibrillar-linking elements, such as thin fibrils that bridge the aligned fibrils or macromolecules such as glycosaminoglycans (GAGs) in the interfibrillar sliding and verify it with a theoretical shear-lag model. We showed the existence of a previously unappreciated structure-function mechanism whereby the Poisson's ratio in tendon is affected by the strain applied and interfibrillar-linker properties, and together these features predict tendon volume shrinkage under tensile loading. During loading, the interfibrillar-linkers pulled fibrils toward each other and squeezed the matrix, leading to the Poisson's ratio larger than 0.5 and fluid expulsion. In addition, the rotation of the interfibrillar-linkers with respect to the fibrils at large strains caused a reduction in the volume shrinkage and eventual nonlinear decay in Poisson's ratio at large strains. Our model also predicts a fluid flow that has a radial pattern toward the surrounding medium, with the larger fluid velocities in proportion to the interfibrillar sliding.

A general framework for application of prestrain to computational models of biological materials

It is often important to include prestress in computational models of biological tissues. The prestress can represent residual stresses (stresses that exist after the tissue is excised from the body) or in situ stresses (stresses that exist in vivo, in the absence of loading). A prestressed reference configuration may also be needed when modeling the reference geometry of biological tissues in vivo. This research developed a general framework for representing prestress in finite element models of biological materials. It is assumed that the material is elastic, allowing the prestress to be represented via a prestrain. For prestrain fields that are not compatible with the reference geometry, the computational framework provides an iterative algorithm for updating the prestrain until equilibrium is satisfied. The iterative framework allows for enforcement of two different constraints: elimination of distortion in order to address the incompatibility issue, and enforcing a specified in situ fiber strain field while allowing for distortion. The framework was implemented as a plugin in FEBio (www.febio.org), making it easy to maintain the software and to extend the framework if needed. Several examples illustrate the application and effectiveness of the approach, including the application of in situ strains to ligaments in the Open Knee model (simtk.org/home/openknee). A novel method for recovering the stress-free configuration from the prestrain deformation gradient is also presented. This general purpose theoretical and computational framework for applying prestrain will allow analysts to overcome the challenges in modeling this important aspect of biological tissue mechanics.

Structural and mechanical characterisation of bridging veins: A review

Bridging veins drain the venous blood from the cerebral cortex into the superior sagittal sinus (SSS) and doing so they bridge the subdural space. Despite their importance in head impact biomechanics, little is known about their properties with respect to histology, morphology and mechanical behaviour. Knowledge of these characteristics is essential for creating a biofidelic finite element model to study the biomechanics of head impact, ultimately leading to the improved design of protective devices by setting up tolerance criteria. This paper presents a comprehensive review of the state-of-the-art knowledge on bridging veins. Tolerance criteria to prevent head injury through impact have been set by a number of research groups, either directly through impact experiments or by means of finite element (FE) simulations. Current state-of-the-art FE head models still lack a biofidelic representation of the bridging veins. To achieve this, a thorough insight into their nature and behaviour is required. Therefore, an overview of the general morphology and histology is provided here, showing the clearly heterogeneous nature of the bridging vein complex, with its three different layers and distinct morphological and histological changes at the region of outflow into the superior sagittal sinus. Apart from a complex morphology, bridging veins also exhibit complex mechanical behaviour, being nonlinear, viscoelastic and prone to damage. Existing material models capable of capturing these properties, as well as methods for experimental characterisation, are discussed. Future work required in bridging vein research is firstly to achieve consensus on aspects regarding morphology and histology, especially in the outflow cuff segment. Secondly, the advised material models need to be populated with realistic parameters through biaxial mechanical experiments adapted to the dimensions of the bridging vein samples. Finally, updating the existing finite element head models with these parameters will render them truly biofidelic, allowing the establishment of accurate tolerance criteria and, ultimately, better head protection devices.

Imaging of current flow in the human head during transcranial electrical therapy

BACKGROUND

It has been assumed that effects caused by tDCS or tACS neuromodulation are due to electric current flow within brain structures. However, to date, direct current density distributions in the brains of human subjects have not been measured. Instead computational models of tDCS or tACS have been used to predict electric current and field distributions for dosimetry and mechanism analysis purposes.

OBJECTIVE/HYPOTHESIS

We present the first in vivo images of electric current density distributions within the brain in four subjects undergoing transcranial electrical stimulation.

METHODS

Magnetic resonance electrical impedance tomography (MREIT) techniques encode current flow in phase images. In four human subjects, we used MREIT to measure magnetic flux density distributions caused by tACS currents, and then calculated current density distributions from these data. Computational models of magnetic flux and current distribution, constructed using contemporaneously collected T₁-weighted structural MRI images, were co-registered to compare predicted and experimental results.

RESULTS

We found consistency between experimental and simulated magnetic flux and current density distributions using transtemporal (T7-T8) and anterior-posterior (Fpz-Oz) electrode montages, and also differences that may indicate a need to improve models to better interpret experimental results. While human subject data agreed with computational model predictions in overall scale, differences may result from factors such as effective electrode surface area and conductivities assumed in models.

CONCLUSIONS

We believe this method may be useful in improving reproducibility, assessing safety, and ultimately aiding understanding of mechanisms of action in electrical and magnetic neuromodulation modalities.

A preliminary modeling investigation into the safe correction zone for high tibial osteotomy

BACKGROUND

High tibial osteotomy (HTO) re-aligns the weight-bearing axis (WBA) of the lower limb. The surgery reduces medial load (reducing pain and slowing progression of cartilage damage) while avoiding overloading the lateral compartment. The optimal correction has not been established. This study investigated how different WBA re-alignments affected load distribution in the knee, to consider the optimal post-surgery re-alignment.

METHODS

We collected motion analysis and seven Tesla MRI data from three healthy subjects, and combined this data to create sets of subject-specific finite element models (total=45 models). Each set of models simulated a range of potential post-HTO knee re-alignments. We shifted the WBA from its native alignment to between 40% and 80% medial-lateral tibial width (corresponding to 2.8°-3.1° varus and 8.5°-9.3° valgus), in three percent increments. We then compared stress/pressure distributions in the models.

RESULTS

Correcting the WBA to 50% tibial width (0° varus-valgus) approximately halved medial compartment stresses, with minimal changes to lateral stress levels, but provided little margin for error in undercorrection. Correcting the WBA to a more commonly-used 62%-65% tibial width (3.4°-4.6° valgus) further reduced medial stresses but introduced the danger of damaging lateral compartment tissues. To balance optimal loading environment with that of the historical risk of under-correction, we propose a new target: WBA correction to 55% tibial width (1.7°-1.9° valgus), which anatomically represented the apex of the lateral tibial spine.

CONCLUSIONS

Finite element models can successfully simulate a variety of HTO re-alignments. Correcting the WBA to 55% tibial width (1.7°-1.9° valgus) optimally distributes medial and lateral stresses/pressures.

Fatigue lifetime prediction of a reduced-diameter dental implant system: Numerical and experimental study

OBJECTIVE

To validate the fatigue lifetime of a reduced-diameter dental implant system predicted by three-dimensional finite element analysis (FEA) by testing physical implant specimens using an accelerated lifetime testing (ALT) strategy with the apparatus specified by ISO 14801.

METHODS

A commercially-available reduced-diameter titanium dental implant system (Straumann Standard Plus NN) was digitized using a micro-CT scanner. Axial slices were processed using an interactive medical image processing software (Mimics) to create 3D models. FEA analysis was performed in ABAQUS, and fatigue lifetime was predicted using fe-safe^® software. The same implant specimens (n=15) were tested at a frequency of 2Hz on load frames using apparatus specified by ISO 14801 and ALT. Multiple step-stress load profiles with various aggressiveness were used to improve testing efficiency. Fatigue lifetime statistics of physical specimens were estimated in a reliability analysis software (ALTA PRO). Fractured specimens were examined using SEM with fractographic technique to determine the failure mode.

RESULTS

FEA predicted lifetime was within the 95% confidence interval of lifetime estimated by experimental results, which suggested that FEA prediction was accurate for this implant system. The highest probability of failure was located at the root of the implant body screw thread adjacent to the simulated bone level, which also agreed with the failure origin in physical specimens.

SIGNIFICANCE

Fatigue lifetime predictions based on finite element modeling could yield similar results in lieu of physical testing, allowing the use of virtual testing in the early stages of future research projects on implant fatigue.

Patient-specific bone modeling and analysis: the role of integration and automation in clinical adoption

Patient-specific analysis of bones is considered an important tool for diagnosis and treatment of skeletal diseases and for clinical research aimed at understanding the etiology of skeletal diseases and the effects of different types of treatment on their progress. In this article, we discuss how integration of several important components enables accurate and cost-effective patient-specific bone analysis, focusing primarily on patient-specific finite element (FE) modeling of bones. First, the different components are briefly reviewed. Then, two important aspects of patient-specific FE modeling, namely integration of modeling components and automation of modeling approaches, are discussed. We conclude with a section on validation of patient-specific modeling results, possible applications of patient-specific modeling procedures, current limitations of the modeling approaches, and possible areas for future research.

Robust augmented reality registration method for localization of solid organs' tumors using CT-derived virtual biomechanical model and fluorescent fiducials

BACKGROUND

Augmented reality (AR) is the fusion of computer-generated and real-time images. AR can be used in surgery as a navigation tool, by creating a patient-specific virtual model through 3D software manipulation of DICOM imaging (e.g., CT scan). The virtual model can be superimposed to real-time images enabling transparency visualization of internal anatomy and accurate localization of tumors. However, the 3D model is rigid and does not take into account inner structures' deformations. We present a concept of automated AR registration, while the organs undergo deformation during surgical manipulation, based on finite element modeling (FEM) coupled with optical imaging of fluorescent surface fiducials.

METHODS

Two 10 × 1 mm wires (pseudo-tumors) and six 10 × 0.9 mm fluorescent fiducials were placed in ex vivo porcine kidneys (n = 10). Biomechanical FEM-based models were generated from CT scan. Kidneys were deformed and the shape changes were identified by tracking the fiducials, using a near-infrared optical system. The changes were registered automatically with the virtual model, which was deformed accordingly. Accuracy of prediction of pseudo-tumors' location was evaluated with a CT scan in the deformed status (ground truth). In vivo: fluorescent fiducials were inserted under ultrasound guidance in the kidney of one pig, followed by a CT scan. The FEM-based virtual model was superimposed on laparoscopic images by automatic registration of the fiducials.

RESULTS

Biomechanical models were successfully generated and accurately superimposed on optical images. The mean measured distance between the estimated tumor by biomechanical propagation and the scanned tumor (ground truth) was 0.84 ± 0.42 mm. All fiducials were successfully placed in in vivo kidney and well visualized in near-infrared mode enabling accurate automatic registration of the virtual model on the laparoscopic images.

CONCLUSIONS

Our preliminary experiments showed the potential of a biomechanical model with fluorescent fiducials to propagate the deformation of solid organs' surface to their inner structures including tumors with good accuracy and automatized robust tracking.

Micropatterned cell sheets as structural building blocks for biomimetic vascular patches

To successfully develop a functional tissue-engineered vascular patch, recapitulating the hierarchical structure of vessel is critical to mimic mechanical properties. Here, we use a cell sheet engineering strategy with micropatterning technique to control structural organization of bovine aortic vascular smooth muscle cell (VSMC) sheets. Actin filament staining and image analysis showed clear cellular alignment of VSMC sheets cultured on patterned substrates. Viability of harvested VSMC sheets was confirmed by Live/Dead^® cell viability assay after 24 and 48 h of transfer. VSMC sheets stacked to generate bilayer VSMC patches exhibited strong inter-layer bonding as shown by lap shear test. Uniaxial tensile testing of monolayer VSMC sheets and bilayer VSMC patches displayed nonlinear, anisotropic stress-stretch response similar to the biomechanical characteristic of a native arterial wall. Collagen content and structure were characterized to determine the effects of patterning and stacking on extracellular matrix of VSMC sheets. Using finite-element modeling to simulate uniaxial tensile testing of bilayer VSMC patches, we found the stress-stretch response of bilayer patterned VSMC patches under uniaxial tension to be predicted using an anisotropic hyperelastic constitutive model. Thus, our cell sheet harvesting system combined with biomechanical modeling is a promising approach to generate building blocks for tissue-engineered vascular patches with structure and mechanical behavior mimicking native tissue.

Collagen fiber interweaving is central to sclera stiffness

The mechanical properties of the microstructural components of sclera are central to eye physiology and pathology. Because these parameters are extremely difficult to measure directly, they are often estimated using inverse-modeling matching deformations of macroscopic samples measured experimentally. Although studies of sclera microstructure show collagen fiber interweaving, current models do not account for this interweaving or the resulting fiber-fiber interactions, which might affect parameter estimates. Our goal was to test the hypothesis that constitutive parameters estimated using inverse modeling differ if models account for fiber interweaving and interactions. We developed models with non-interweaving or interweaving fibers over a wide range of volume fractions (36-91%). For each model, we estimated fiber stiffness using inverse modeling matching biaxial experimental data of human sclera. We found that interweaving increased the estimated fiber stiffness. When the collagen volume fraction was 64% or less, the stiffness of interweaving fibers was about 1.25 times that of non-interweaving fibers. For higher volume fractions, the ratio increased substantially, reaching 1.88 for a collagen volume fraction of 91%. Simulating a model (interweaving/non-interweaving) using the fiber stiffness estimated from the other model produced substantially different behavior, far from that observed experimentally. These results show that estimating microstructural component mechanical properties is highly sensitive to the assumed interwoven/non-interwoven architecture. Moreover, the results suggest that interweaving plays an important role in determining the structural stiffness of sclera, and potentially of other soft tissues in which the collagen fibers interweave. STATEMENT OF SIGNIFICANCE: The collagen fibers of sclera are interwoven, but numerical models do not account for this interweaving or the resulting fiber-fiber interactions. To determine if interweaving matters, we examined the differences in the constitutive model parameters estimated using inverse modeling between models with interweaving and non-interweaving fibers. We found that the estimated stiffness of the interweaving fibers was up to 1.88 times that of non-interweaving fibers, and that the estimate increased with collagen volume fraction. Our results suggest that fiber interweaving is a fundamental characteristic of connective tissues, additional to anisotropy, density and orientation. Better characterization of interweaving, and of its mechanical effects is likely central to understanding microstructure and biomechanics of sclera and other soft tissues.

The influence of running shoe with different carbon-fiber plate designs on internal foot mechanics: A pilot computational analysis

A carbon-fiber plate (CFP) embedded into running shoes is a commonly applied method to improve running economy, but little is known in regard the effects of CFP design features on internal foot mechanics. This study aimed to explore how systematic changes in CFP geometrical variations (i.e., thickness and location) can alter plantar pressure and strain under the forefoot as well as metatarsal stress state through computational simulations. A foot-shoe finite element (FE) model was built and different CFP features including three thicknesses (1 mm, 2 mm, and 3 mm) and three placements (high-loaded (just below the insole), mid-loaded (in between the midsole), and low-loaded (just above the outsole)) were further modulated within the shoe sole. Simulations were conducted at the impact peak instant during forefoot strike running. Compared with the no-CFP shoe, peak plantar pressure and compressive strain under the forefoot consistently decreased when the CFP thickness increased, and the low-loaded conditions were found more effective (peak pressure decreased up to 31.91% and compressive strain decreased up to 18.61%). In terms of metatarsal stress, CFP designs resulted in varied effects and were dependent on their locations. Specifically, high-loaded CFP led to relatively higher peak metatarsal stress without the reduction trend as thickness increased (peak stress increased up to 12.91%), while low-loaded conditions showed a gradual reduction in peak stress, decreasing by 0.74%. Therefore, a low-loaded thicker CFP should be considered to achieve the pressure-relief effects of running shoes without the expense of increased metatarsal stress.

3D correction over 2years with anterior vertebral body growth modulation: A finite element analysis of screw positioning, cable tensioning and postoperative functional activities

BACKGROUND

Anterior vertebral body growth modulation is a fusionless instrumentation to correct scoliosis using growth modulation. The objective was to biomechanically assess effects of cable tensioning, screw positioning and post-operative position on tridimensional correction.

METHODS

The design of experiments included two variables: cable tensioning (150/200N) and screw positioning (lateral/anterior/triangulated), computationally tested on 10 scoliotic cases using a personalized finite element model to simulate spinal instrumentation, and 2years growth modulation with the device. Dependent variables were: computed Cobb angles, kyphosis, lordosis, axial rotation and stresses exerted on growth plates. Supine functional post-operative position was simulated in addition to the reference standing position to evaluate corresponding growth plate's stresses.

FINDINGS

Simulated cable tensioning and screw positioning had a significant impact on immediate and after 2years Cobb angle (between 5°-11°, p<0.01). Anterior screw positioning significantly increased kyphosis after 2years (6°-8°, p=0.02). Triangulated screw positioning did not significantly impact axial rotation but significantly reduced kyphosis (8°-10°, p=0.001). Growth plates' stresses were increased by 23% on the curve's convex side with cable tensioning, while screw positioning rather affected anterior/posterior distributions. Supine position significantly affected stress distributions on the apical vertebra compared to standing position (respectively 72% of compressive stresses on convex side vs 55%).

INTERPRETATION

This comparative numerical study showed the biomechanical possibility to adjust the fusionless instrumentation parameters to improve correction in frontal and sagittal planes, but not in the transverse plane. The convex side stresses increase in the supine position may suggest that growth modulation could be accentuated during nighttime.

Mechanical damage of tympanic membrane in relation to impulse pressure waveform - A study in chinchillas

Mechanical damage to middle ear components in blast exposure directly causes hearing loss, and the rupture of the tympanic membrane (TM) is the most frequent injury of the ear. However, it is unclear how the severity of injury graded by different patterns of TM rupture is related to the overpressure waveforms induced by blast waves. In the present study, the relationship between the TM rupture threshold and the impulse or overpressure waveform has been investigated in chinchillas. Two groups of animals were exposed to blast overpressure simulated in our lab under two conditions: open field and shielded with a stainless steel cup covering the animal head. Auditory brainstem response (ABR) and wideband tympanometry were measured before and after exposure to check the hearing threshold and middle ear function. Results show that waveforms recorded in the shielded case were different from those in the open field and the TM rupture threshold in the shielded case was lower than that in the open field (3.4 ± 0.7 vs. 9.1 ± 1.7 psi or 181 ± 1.6 vs. 190 ± 1.9 dB SPL). The impulse pressure energy spectra analysis of waveforms demonstrates that the shielded waveforms include greater energy at high frequencies than that of the open field waves. Finally, a 3D finite element (FE) model of the chinchilla ear was used to compute the distributions of stress in the TM and the TM displacement with impulse pressure waves. The FE model-derived change of stress in response to pressure loading in the shielded case was substantially faster than that in the open case. This finding provides the biomechanical mechanisms for blast induced TM damage in relation to overpressure waveforms. The TM rupture threshold difference between the open and shielded cases suggests that an acoustic role of helmets may exist, intensifying ear injury during blast exposure.