Estimation of load conditions and strain distribution for in vivo murine tibia compression loading using experimentally informed finite element models.
Ontology highlight
ABSTRACT:
SUBMITTER: Pickering E
PROVIDER: S-EPMC7856106 | biostudies-literature |
REPOSITORIES: biostudies-literature
Similar Datasets
Project description:High-resolution peripheral quantitative computed tomography (HR-pQCT) based micro-finite element (μFE) analysis allows accurate prediction of stiffness and ultimate load of standardised (∼1 cm) distal radius and tibia sections. An alternative homogenized finite element method (hFE) was recently validated to compute the ultimate load of larger (∼2 cm) distal radius sections that include Colles' fracture sites. Since the mechanical integrity of the weight-bearing distal tibia is gaining clinical interest, it has been shown that the same properties can be used to predict the strength of both distal segments of the radius and the tibia. Despite the capacity of hFE to predict structural properties of distal segments of the radius and the tibia, the limitations of such homogenization scheme remain unclear. Therefore, the objective of this study is to build a complete mechanical data set of the compressive behavior of distal segments of the tibia and to compare quantitatively the structural properties with the hFE predictions. As a further aim, it is intended to verify whether hFE is also able to capture the post-yield strain localisation or fracture zones in such a bone section, despite the absence of strain softening in the constitutive model. Twenty-five fresh-frozen distal parts of tibias of human donors were used in this study. Sections were cut corresponding to an in-house triple-stack protocol HR-pQCT scan, lapped, and scanned using micro computed tomography (μCT). The sections were tested in compression until failure, unloaded and scanned again in μCT. Volumetric bone mineral density (vBMD) and bone mineral content (BMC) were correlated to compression test results. hFE analysis was performed in order to compare computational predictions (stiffness, yield load and plastic deformation field pattern) with the compressive experiment. Namely, strain localization was assessed based on digital volume correlation (DVC) results and qualitatively compared to hFE predictions by comparing mid-slices patterns. Bone mineral content (BMC) showed a good correlation with stiffness (R2 = 0.92) and yield (R2 = 0.88). Structural parameters also showed good agreement between the experiment and hFE for both stiffness (R2 = 0.96, slope = 1.05 with 95 % CI [0.97, 1.14]) and yield (R2 = 0.95, slope = 1.04 [0.94, 1.13]). The qualitative comparison between hFE and DVC strain localization patterns allowed the classification of the samples into 3 categories: bad (15 sections), semi (8), and good agreement (2). The good correlations between BMC or hFE and experiment for structural parameters were similar to those obtained previously for the distal part of the radius. The failure zones determined by hFE corresponded to registration only in 8 % of the cases. We attribute these discrepancies to local elastic/plastic buckling effects that are not captured by the continuum-based FE approach exempt from strain softening. A way to improve strain localization hFE prediction would be to use longer distal segments with intact cortical shells, as done for the radius. To conclude, the used hFE scheme captures the elastic and yield response of the tibia sections reliably but not the subsequent failure process.
| S-EPMC10999809 | biostudies-literature
Project description:High-resolution finite element models derived from micro-computed tomography images are often used to study the effects of trabecular microarchitecture and loading mode on tissue stress, but the degree to which existing finite element methods correctly predict the location of tissue failure is not well characterized. In the current study, we determined the relationship between the location of highly strained tissue, as determined from high-resolution finite element models, and the location of tissue microdamage, as determined from three-dimensional fluoroscopy imaging, which was performed after the microdamage was generated in-vitro by mechanical testing. Fourteen specimens of human vertebral cancellous bone were assessed (8 male donors, 2 female donors, 47-78 years of age). Regions of stained microdamage, were 50-75% more likely to form in highly strained tissue (principal strains exceeding 0.4%) than elsewhere, and generally the locations of the regions of microdamage were significantly correlated (p<0.05) with the locations of highly strained tissue. This spatial correlation was stronger for the largest regions of microdamage (?1,000,000?m(3) in volume); 87% of large regions of microdamage were located near highly strained tissue. Together, these findings demonstrate that there is a strong correlation between regions of microdamage and regions of high strain in human cancellous bone, particularly for the biomechanically more important large instances of microdamage.
| S-EPMC4663152 | biostudies-literature
Project description:Our objective was to identify precise mechanical metrics of the proximal tibia which differentiated OA and normal knees. We developed subject-specific FE models for 14 participants (7 OA, 7 normal) who were imaged three times each for assessing precision (repeatability). We assessed various mechanical metrics (minimum principal and von Mises stress and strain as well as structural stiffness) across the proximal tibia for each subject. In vivo precision of these mechanical metrics was assessed using CV%RMS. We performed parametric and non-parametric statistical analyses and determined Cohen's d effect sizes to explore differences between OA and normal knees. For all FE-based mechanical metrics, average CV%RMS was less than 6%. Minimum principal stress was, on average, 75% higher in OA versus normal knees while minimum principal strain values did not differ. No difference was observed in structural stiffness. FE modeling could precisely quantify and differentiate mechanical metrics variations in normal and OA knees, in vivo. This study suggests that bone stress patterns may be important for understanding OA pathogenesis at the knee.
| S-EPMC6068127 | biostudies-literature
Project description:New treatments for bone diseases require testing in animal models before clinical translation, and the mouse tibia is among the most common models. In vivo micro-Computed Tomography (microCT)-based micro-Finite Element (microFE) models can be used for predicting the bone strength non-invasively, after proper validation against experimental data. Different modelling techniques can be used to estimate the bone properties, and the accuracy associated with each is unclear. The aim of this study was to evaluate the ability of different microCT-based microFE models to predict the mechanical properties of the mouse tibia under compressive load. Twenty tibiae were microCT scanned at 10.4 µm voxel size and subsequently compressed at 0.03 mm/s until failure. Stiffness and failure load were measured from the load-displacement curves. Different microFE models were generated from each microCT image, with hexahedral or tetrahedral mesh, and homogeneous or heterogeneous material properties. Prediction accuracy was comparable among models. The best correlations between experimental and predicted mechanical properties, as well as lower errors, were obtained for hexahedral models with homogeneous material properties. Experimental stiffness and predicted stiffness were reasonably well correlated (R2 = 0.53-0.65, average error of 13-17%). A lower correlation was found for failure load (R2 = 0.21-0.48, average error of 9-15%). Experimental and predicted mechanical properties normalized by the total bone mass were strongly correlated (R2 = 0.75-0.80 for stiffness, R2 = 0.55-0.81 for failure load). In conclusion, hexahedral models with homogeneous material properties based on in vivo microCT images were shown to best predict the mechanical properties of the mouse tibia.
| S-EPMC8154847 | biostudies-literature
Project description:Although previous studies suggested that the interfacial tension gamma(cc) acting along cell-cell boundaries and the effective viscosity mu of the cell cytoplasm could be measured by compressing a spherical aggregate of cells between parallel plates, the mechanical understanding necessary to extract this information from these tests-tests that have provided the surface tension sigma(cm) acting along cell-medium interfaces-has been lacking. These tensions can produce net forces at the subcellular level and give rise to cell motions and tissue reorganization, the rates of which are regulated by mu. Here, a three-dimensional (3D) cell-based finite element model provides insight into the mechanics of the compression test, where these same forces are at work, and leads to quantitative relationships from which the effective viscosity mu of the cell cytoplasm, the tension gamma(cc) that acts along internal cell-cell interfaces and the surface tension sigma(cp) along the cell-platen boundaries can be determined from force-time curves and aggregate profiles. Tests on 5-day embryonic chick mesencephalon, neural retina, liver, and heart aggregates show that all of these properties vary significantly with cell type, except gamma(cc), which is remarkably constant. These properties are crucial for understanding cell rearrangement and tissue self-organization in contexts that include embryogenesis, cancer metastases, and tissue engineering.
| S-EPMC2799989 | biostudies-literature
Project description:Introduction: Three-dimensional (3D)-printed custom pelvic implants have become a clinically viable option for patients undergoing pelvic cancer surgery with resection of the hip joint. However, increased clinical utilization has also necessitated improved implant durability, especially with regard to the compression screws used to secure the implant to remaining pelvic bone. This study evaluated six different finite element (FE) screw modeling methods for predicting compression screw pullout and fatigue failure in a custom pelvic implant secured to bone using nine compression screws. Methods: Three modeling methods (tied constraints (TIE), bolt load with constant force (BL-CF), and bolt load with constant length (BL-CL)) generated screw axial forces using functionality built into Abaqus FE software; while the remaining three modeling methods (isotropic pseudo-thermal field (ISO), orthotropic pseudo-thermal field (ORT), and equal-and-opposite force field (FOR)) generated screw axial forces using iterative physics-based relationships that can be implemented in any FE software. The ability of all six modeling methods to match specified screw pretension forces and predict screw pullout and fatigue failure was evaluated using an FE model of a custom pelvic implant with total hip replacement. The applied hip contact forces in the FE model were estimated at two locations in a gait cycle. For each of the nine screws in the custom implant FE model, likelihood of screw pullout failure was predicted using maximum screw axial force, while likelihood of screw fatigue failure was predicted using maximum von Mises stress. Results: The three iterative physics-based modeling methods and the non-iterative Abaqus BL-CL method produced nearly identical predictions for likelihood of screw pullout and fatigue failure, while the other two built-in Abaqus modeling methods yielded vastly different predictions. However, the Abaqus BL-CL method required the least computation time, largely because an iterative process was not needed to induce specified screw pretension forces. Of the three iterative methods, FOR required the fewest iterations and thus the least computation time. Discussion: These findings suggest that the BL-CL screw modeling method is the best option when Abaqus is used for predicting screw pullout and fatigue failure in custom pelvis prostheses, while the iterative physics-based FOR method is the best option if FE software other than Abaqus is used.
| S-EPMC11372789 | biostudies-literature
Project description:Radiofrequency ablation is the most common minimally invasive therapy used in the United States to treat hepatocellular carcinoma and liver metastases. The ability to perform real-time temperature imaging while a patient is undergoing ablation therapy may help reduce the high recurrence rates following ablation therapy. Ultrasound echo signals undergo time shifts with increasing temperature due to sound speed and thermal expansion, which are tracked using both 1D cross correlation and 2D block matching based speckle tracking methods. In this paper, we present a quantitative evaluation of the accuracy and precision of temperature estimation using the above algorithms on both simulated and experimental data. A finite element analysis simulation of radiofrequency ablation of hepatic tissue was developed. Finite element analysis provides a method to obtain the exact temperature distribution along with a mapping of the tissue displacement due to thermal expansion. These local displacement maps were combined with the displacement due to speed of sound changes and utilized to generate ultrasound radiofrequency frames at specified time increments over the entire ablation procedure. These echo signals provide an ideal test-bed to evaluate the performance of both speckle tracking methods, since the estimated temperature results can be compared directly to the exact finite element solution. Our results indicate that the 1D cross-correlation (CC) method underestimates the cumulative displacement by 0.20mm, while the underestimation with 2D block matching (BM) is about 0.14 mm after 360 s of ablation. The 1D method also overestimates the size of the ablated region by 5.4% when compared to 2.4% with the 2D method after 720 s of ablation. Hence 2D block matching provides better tracking of temperature variations when compared to the 1D cross-correlation method over the entire duration of the ablation procedure. In addition, results obtained using 1D cross-correlation diverge from the ideal finite element results after 7 min of ablation and for temperatures greater than 65 degrees C. In a similar manner, experimental results presented using a tissue-mimicking phantom also demonstrate that the maximum percent difference with 2D block matching was 5%, when compared to 31% with the 1D method over the 700 s heating duration on the phantom.
| S-EPMC2777513 | biostudies-literature
Project description:Tendon's viscoelastic behaviors are important to the tissue mechanical function and cellular mechanobiology. When loaded in longitudinal tension, tendons often have a large Poisson's ratio (ν>2) that exceeds the limit of incompressibility for isotropic material (ν=0.5), indicating that tendon experiences volume loss, inducing poroelastic fluid exudation in the transverse direction. Therefore, transverse poroelasticity is an important contributor to tendon material behavior. Tendon hydraulic permeability which is required to evaluate the fluid flow contribution to viscoelasticity, is mostly unavailable, and where available, varies by several orders of magnitude. In this manuscript, we quantified the transverse poroelastic material parameters of rat tail tendon fascicles by conducting transverse osmotic loading experiments, in both tension and compression. We used a multi-start optimization method to evaluate the parameters using biphasic finite element modeling. Our tendon samples had a transverse hydraulic permeability of 10-4 to 10-5 mm4. (Ns)-1 and showed a significant tension-compression nonlinearity in the transverse direction. Further, using these results, we predict hydraulic permeability during longitudinal (fiber-aligned) tensile loading, and the spatial distribution of fluid flow during osmotic loading. These results reveal novel aspects of tendon mechanics and can be used to study the physiomechanical response of tendon in response to mechanical loading.
| S-EPMC7438606 | biostudies-literature
Project description:The aim of this study was to investigate the role of the superficial zone on the mechanical behavior of articular cartilage. Confined compression of articular cartilage was modeled using a biphasic finite element analysis to calculate the one-dimensional deformation of the extracellular matrix (ECM) and movement of the interstitial fluid through the ECM and articular surface. The articular cartilage was modeled as an inhomogeneous, nonlinear hyperelastic biphasic material with depth and strain-dependent material properties. Two loading conditions were simulated, one where the superficial zone was loaded with a porous platen (normal test) and the other where the deep zone was loaded with the porous platen (upside down test). Compressing the intact articular cartilage with 0.2 MPa stress reduced the surface permeability by 88%. Removing the superficial zone increased the rate of change for all mechanical parameters and decreased the fluid support ratio of the tissue, resulting in increased tissue deformation. Apparent permeability linearly increased after superficial removal in the normal test, yet it did not change in the upside down test. Orientation of the specimen affected the time-dependent biomechanical behavior of the articular cartilage, but not equilibrium behavior. The two tests with different specimen orientations resulted in very different apparent permeabilities, suggesting that in an experimental study which quantifies material properties of an inhomogeneous material, the specimen orientation should be stated along with the permeability result. The current study provides new insights into the role of the superficial zone on mechanical behavior of the articular cartilage.
| S-EPMC4274183 | biostudies-literature
Project description:The mechanical behavior of aluminum hexagonal honeycombs subjected to out-of-plane dynamic indentation and compression loads has been investigated numerically using ANSYS/LS-DYNA in this paper. The finite element (FE) models have been verified by previous experimental results in terms of deformation pattern, stress-strain curve, and energy dissipation. The verified FE models have then been used in comprehensive numerical analysis of different aluminum honeycombs. Plateau stress, σpl, and dissipated energy (EI for indentation and EC for compression) have been calculated at different strain rates ranging from 10² to 10⁴ s-1. The effects of strain rate and t/l ratio on the plateau stress, dissipated energy, and tearing energy have been discussed. An empirical formula is proposed to describe the relationship between the tearing energy per unit fracture area, relative density, and strain rate for honeycombs. Moreover, it has been found that a generic formula can be used to describe the relationship between tearing energy per unit fracture area and relative density for both aluminum honeycombs and foams.
| S-EPMC5456648 | biostudies-other