Finite element modeling with patient-specific geometry to assess clinical risks of percutaneous pulmonary valve implantation.

BACKGROUND: Percutaneous pulmonary valve implantation (PPVI) is a non-surgical treatment for right ventricular outflow tract (RVOT) dysfunction. During PPVI, a stented valve, delivered via catheter, replaces the dysfunctional pulmonary valve. Stent oversizing allows valve anchoring within the RVOT, but overexpansion can intrude on the surrounding structures. Potentially dangerous outcomes include aortic valve insufficiency (AVI) from aortic root (AR) distortion and myocardial ischemia from coronary artery (CA) compression. Currently, risks are evaluated via balloon angioplasty/sizing before stent deployment. Patient-specific finite element (FE) analysis frameworks can improve pre-procedural risk assessment, but current methods require hundreds of hours of high-performance computation. METHODS: We created a simplified method to simulate the procedure using patient-specific FE models for accurate, efficient pre-procedural PPVI (using balloon expandable valves) risk assessment. The methodology was tested by retrospectively evaluating the clinical outcome of 12 PPVI candidates. RESULTS: Of 12 patients (median age 14.5 years) with dysfunctional RVOT, 7 had native RVOT and 5 had RV-PA conduits. Seven patients had undergone successful RVOT stent/valve placement, three had significant AVI on balloon testing, one had left CA compression, and one had both AVI and left CA compression. A model-calculated change of more than 20% in lumen diameter of the AR or coronary arteries correctly predicted aortic valve sufficiency and/or CA compression in all the patients. CONCLUSION: Agreement between FE results and clinical outcomes is excellent. Additionally, these models run in 2-6 min on a desktop computer, demonstrating potential use of FE analysis for pre-procedural risk assessment of PPVI in a clinically relevant timeframe.

An Approach to Quantify Anisotropic Multiaxial Failure of the Annulus Fibrosus.

Tears in the annulus fibrosus (AF) of the intervertebral disk (IVD) occur due to multiaxial loading on the spine. However, most existing AF failure studies measure uniaxial stress, not the multiaxial stress at failure. Delamination theory, which requires advanced structural knowledge and knowledge about the interactions between the AF fibers and matrix, has historically been used to understand and predict AF failure. Alternatively, a simple method, the Tsai-Hill yield criteria, could describe multiaxial failure of the AF. This yield criteria uses the known tissue fiber orientation and an equation to establish the multiaxial failure stresses that cause failure. This paper presents a method to test the multiaxial failure stress of the AF experimentally and evaluate the potential for the Tsai-Hill model to predict these failure stresses. Porcine AF was cut into a dogbone shape at three distinct angles relative to the primary lamella direction (parallel, transverse, and oblique). Then, each dogbone was pulled to complete rupture. The Cauchy stress in the material's fiber coordinates was calculated. These multiaxial stress parameters were used to optimize the coefficients of the Tsai-Hill yield. The coefficients obtained for the Tsai-Hill model vary by an order of magnitude between the fiber and transverse directions, and these coefficients provide a good description of the AF multiaxial failure stress. These results establish both an experimental approach and the use of the Tsai-Hill model to explain the anisotropic failure behavior of the tissue.

A Graphical Approach to Visualize and Interpret Biochemically Coupled Biomechanical Models.

The last decade has seen the emergence of progressively more complex mechanobiological models, often coupling biochemical and biomechanical components. The complexity of these models makes interpretation difficult, and although computational tools can solve model equations, there is considerable potential value in a simple method to explore the interplay between different model components. Pump and system performance curves, long utilized in centrifugal pump selection and design, inspire the development of a graphical technique to depict visually the performance of biochemically-coupled mechanical models. Our approach is based on a biochemical performance curve (analogous to the classical pump curve) and a biomechanical performance curve (analogous to the system curve). Upon construction of the two curves, their intersection, or lack thereof, describes the coupled model's equilibrium state(s). One can also observe graphically how an applied perturbation shifts one or both curves, and thus how the other component will respond, without rerunning the full model. While the upfront cost of generating the performance curve graphic varies with the efficiency of the model components, the easily interpretable visual depiction of what would otherwise be nonintuitive model behavior is valuable. Herein, we outline how performance curves can be constructed and interpreted for biochemically-coupled biomechanical models and apply the technique to two independent models in the cardiovascular space. The performance curve approach can illustrate and help identify weaknesses in model construction, inform user-applied perturbations and fitting procedures to generate intended behaviors, and improve the efficiency of the model generation and application process.

Atherosclerotic Calcifications Have a Local Effect on the Peel Behavior of Human Aortic Media.

Aortic dissections, characterized by the propagation of a tear through the layers of the vessel wall, are critical, life-threatening events. Aortic calcifications are a common comorbidity in both acute and chronic dissections, yet their impact on dissection mechanics remains unclear. Using micro-computed tomography (CT) imaging, peel testing, and finite element modeling, this study examines the interplay between atherosclerotic calcifications and dissection mechanics. Samples cut from cadaveric human thoracic aortas were micro-CT imaged and subsequently peel-tested to map peel tension curves to the location of aortic calcifications. Empirical mode decomposition separated peel tension curves into high and low-frequency components, with high-frequency effects corresponding to interlamellar bonding mechanics and low-frequency effects to peel tension fluctuations. Finally, we used an idealized finite element model to examine how stiff calcifications affect aortic failure mechanics. Results showed that atherosclerosis influences dissection behavior on multiple length scales. Experimentally, atherosclerotic samples exhibited higher peel tensions and greater variance in the axial direction. The variation was driven by increased amplitudes of low-frequency tension fluctuations in diseased samples, indicating that more catastrophic propagations occur near calcifications. The simulations corroborated this finding, suggesting that the low-frequency changes resulted from the presence of a stiff calcification in the vessel wall. There were also modifications to the high-frequency peel mechanics, a response likely attributable to alterations in the microstructure and interlamellar bonding within the media. Considered collectively, these findings demonstrate that dissection mechanics are modified in aortic media nearby and adjacent to aortic calcifications.

The Lumbar Facet Capsular Ligament Becomes More Anisotropic and the Fibers Become Stiffer With Intervertebral Disc and Facet Joint Degeneration.

Degeneration of the lumbar spine, and especially how that degeneration may lead to pain, remains poorly understood. In particular, the mechanics of the facet capsular ligament may contribute to low back pain, but the mechanical changes that occur in this ligament with spinal degeneration are unknown. Additionally, the highly nonlinear, heterogeneous, and anisotropic nature of the facet capsular ligament makes understanding mechanical changes more difficult. Clinically, magnetic resonance imaging (MRI)-based signs of degeneration in the facet joint and the intervertebral disc (IVD) correlate. Therefore, this study examined how the nonlinear, heterogeneous mechanics of the facet capsular ligament change with degeneration of the lumbar spine as characterized using MRI. Cadaveric human spines were imaged via MRI, and the L2-L5 facet joints and IVDs were scored using the Fujiwara and Pfirrmann grading systems. Then, the facet capsular ligament was isolated and biaxially loaded. The nonlinear mechanical properties of the ligament were obtained using a nonlinear generalized anisotropic inverse mechanics analysis (nGAIM). Then a Holzapfel-Gasser-Ogden (HGO) model was fit to the stress-strain data obtained from nGAIM. The facet capsular ligament is stiffer and more anisotropic at larger Pfirrmann grades and higher Fujiwara scores than at lower grades and scores. Analysis of ligament heterogeneity showed all tissues are highly heterogeneous, but no distinct spatial patterns of heterogeneity were found. These results show that degeneration of the lumbar spine including the facet capsular ligament appears to be occurring as a whole joint phenomenon and advance our understanding of lumbar spine degeneration.

Characterizing Tissue Remodeling and Mechanical Heterogeneity in Cerebral Aneurysms.

Accurately assessing the complex tissue mechanics of cerebral aneurysms (CAs) is critical for elucidating how CAs grow and whether that growth will lead to rupture. The factors that have been implicated in CA progression - blood flow dynamics, immune infiltration, and extracellular matrix remodeling - all occur heterogeneously throughout the CA. Thus, it stands to reason that the mechanical properties of CAs are also spatially heterogeneous. Here, we present a new method for characterizing the mechanical heterogeneity of human CAs using generalized anisotropic inverse mechanics, which uses biaxial stretching experiments and inverse analyses to determine the local Kelvin moduli and principal alignments within the tissue. Using this approach, we find that there is significant mechanical heterogeneity within a single acquired human CA. These results were confirmed using second harmonic generation imaging of the CA's fiber architecture and a correlation was observed. This approach provides a single-step method for determining the complex heterogeneous mechanics of CAs, which has important implications for future identification of metrics that can improve accuracy in prediction risk of rupture.

In Situ Lumbar Facet Capsular Ligament Strains Due to Joint Pressure and Residual Strain.

The lumbar facet capsular ligament, which surrounds and limits the motion of each facet joint in the lumbar spine, has been recognized as being mechanically significant and has been the subject of multiple mechanical characterization studies in the past. Those studies, however, were performed on isolated tissue samples and thus could not assess the mechanical state of the ligament in vivo, where the constraints of attachment to rigid bone and the force of the joint pressure lead to nonzero strain even when the spine is not loaded. In this work, we quantified these two effects using cadaveric lumbar spines (five spines, 20 total facet joints harvested from L2 to L5). The effect of joint pressure was measured by injecting saline into the joint space and tracking the 3D capsule surface motion via digital image correlation, and the prestrain due to attachment was measured by dissecting a large section of the tissue from the bone and by tracking the motion between the on-bone and free states. We measured joint pressures of roughly 15-40 kPa and local first principal strains of up to 25-50% when 0.3 mL of saline was injected into the joint space; the subsequent increase in pressure and strain were more modest for further increases in injection volume, possibly due to leakage of fluid from the joint. The largest stretches were in the bone-to-bone direction in the portions of the ligament spanning the joint space. When the ligament was released from the vertebrae, it shrank by an average of 4-5%, with local maximum (negative) principal strain values of up to 30%, on average. Based on these measurements and previous tests on isolated lumbar facet capsular ligaments, we conclude that the normal in vivo state of the facet capsular ligament is in tension, and that the collagen in the ligament is likely uncrimped even when the spine is not loaded.

An Experimental-Computational Approach to Quantify Blood Rheology in Sickle Cell Disease.

In sickle cell disease, aberrant blood flow due to oxygen-dependent changes in red cell biomechanics is a key driver of pathology. Most studies to date have focused on the potential role of altered red cell deformability and blood rheology in precipitating vaso-occlusive crises. Numerous studies, however, have shown that sickle blood flow is affected even at high oxygen tensions, suggesting a potentially systemic role for altered blood flow in driving pathologies, including endothelial dysfunction, ischemia, and stroke. In this study, we applied a combined experimental-computation approach that leveraged an experimental platform that quantifies sickle blood velocity fields under a range of oxygen tensions and shear rates. We computationally fitted a continuum model to our experimental data to generate physics-based parameters that capture patient-specific rheological alterations. Our results suggest that sickle blood flow is altered systemically, from the arterial to the venous circulation. We also demonstrated the application of this approach as a tool to design patient-specific transfusion regimens. Finally, we demonstrated that patient-specific rheological parameters can be combined with patient-derived vascular models to identify patients who are at higher risk for cerebrovascular complications such as aneurysm and stroke. Overall, this study highlights that sickle blood flow is altered systemically, which can drive numerous pathologies, and this study demonstrates the potential utility of an experimentally parameterized continuum model as a predictive tool for patient-specific care.

Investigation of Pathophysiological Aspects of Aortic Growth, Remodeling, and Failure Using a Discrete-Fiber Microstructural Model.

Aortic aneurysms are inherently unpredictable. One can never be sure whether any given aneurysm may rupture or dissect. Clinically, the criteria for surgical intervention are based on size and growth rate, but it remains difficult to identify a high-risk aneurysm, which may require intervention before the cutoff criteria, versus an aneurysm than can be treated safely by more conservative measures. In this work, we created a computational microstructural model of a medial lamellar unit (MLU) incorporating (1) growth and remodeling laws applied directly to discrete, individual fibers, (2) separate but interacting fiber networks for collagen, elastin, and smooth muscle, (3) active and passive smooth-muscle cell mechanics, and (4) failure mechanics for all three fiber types. The MLU model was then used to study different pathologies and microstructural anomalies that may play a role in vascular growth and failure. Our model recapitulated many aspects of arterial remodeling under hypertension with no underlying genetic syndrome including remodeling dynamics, tissue mechanics, and failure. Syndromic effects (smooth muscle cell (SMC) dysfunction or elastin fragmentation) drastically changed the simulated remodeling process, tissue behavior, and tissue strength. Different underlying pathologies were able to produce similarly dilatated vessels with different failure properties, providing a partial explanation for the imperfect nature of aneurysm size as a predictor of outcome.

Dicer1 Deficiency in the Idiopathic Pulmonary Fibrosis Fibroblastic Focus Promotes Fibrosis by Suppressing MicroRNA Biogenesis.

RATIONALE: The lung extracellular matrix (ECM) in idiopathic pulmonary fibrosis (IPF) mediates progression of fibrosis by decreasing fibroblast expression of miR-29 (microRNA-29), a master negative regulator of ECM production. The molecular mechanism is undefined. IPF-ECM is stiffer than normal. Stiffness drives fibroblast ECM production in a YAP (yes-associated protein)-dependent manner, and YAP is a known regulator of miR-29. Therefore, we tested the hypothesis that negative regulation of miR-29 by IPF-ECM was mediated by mechanotransduction of stiffness. OBJECTIVES: To determine how IPF-ECM negatively regulates miR-29. METHODS: We decellularized lung ECM using detergents and prepared polyacrylamide hydrogels of defined stiffness by varying acrylamide concentrations. Mechanistic studies were guided by immunohistochemistry of IPF lung and used cell culture, RNA-binding protein assays, and xenograft models. MEASUREMENTS AND MAIN RESULTS: Contrary to our hypothesis, we excluded fibroblast mechanotransduction of ECM stiffness as the primary mechanism deregulating miR-29. Instead, systematic examination of miR-29 biogenesis revealed a microRNA processing defect that impeded processing of miR-29 into its mature bioactive forms. Immunohistochemical analysis of the microRNA processing machinery in IPF lung specimens revealed decreased Dicer1 expression in the procollagen-rich myofibroblastic core of fibroblastic foci compared with the focus perimeter and adjacent alveolar walls. Mechanistically, IPF-ECM increased association of the Dicer1 transcript with RNA binding protein AUF1 (AU-binding factor 1), and Dicer1 knockdown conferred primary human lung fibroblasts with cell-autonomous fibrogenicity in zebrafish and mouse lung xenograft models. CONCLUSIONS: Our data identify suppression of fibroblast Dicer1 expression in the myofibroblast-rich IPF fibroblastic focus core as a central step in the mechanism by which the ECM sustains fibrosis progression in IPF.

Editors' Choice Papers for 2018.

No abstract.

Effects of Collagen Heterogeneity on Myocardial Infarct Mechanics in a Multiscale Fiber Network Model.

The scar that forms after a myocardial infarction is often characterized by a highly disordered architecture but generally exhibits some degree of collagen fiber orientation, with a resulting mechanical anisotropy. When viewed in finer detail, however, the heterogeneity of the sample is clear, with different subregions exhibiting different fiber orientations. In this work, we used a multiscale finite element model to explore the consequences of the heterogeneity in terms of mechanical behavior. To do so, we used previously obtained fiber alignment maps of rat myocardial scar slices (n = 15) to generate scar-specific finite element meshes that were populated with fiber models based on the local alignment state. These models were then compared to isotropic models with the same sample shape and fiber density, and to homogeneous models with the same sample shape, fiber density, and average fiber alignment as the scar-specific models. All simulations involved equibiaxial extension of the sample with free motion in the third dimension. We found that heterogeneity led to a lower degree of mechanical anisotropy and a higher level of local stress concentration than the corresponding homogeneous model, and also that fibers failed in the heterogeneous model at much lower macroscopic strains than in the isotropic and homogeneous models. Taken together, these results suggest that scar heterogeneity may impair myocardial mechanical function both in terms of anisotropy and strength, and that individual variations in scar heterogeneity could be an important consideration for understanding scar remodeling and designing therapeutic interventions for patients after myocardial infarction.

Ex Vivo Mechanical Tests and Multiscale Computational Modeling Highlight the Importance of Intramural Shear Stress in Ascending Thoracic Aortic Aneurysms.

Ascending thoracic aortic aneurysms (ATAAs) are anatomically complex in terms of architecture and geometry, and both complexities contribute to unpredictability of ATAA dissection and rupture in vivo. The goal of this work was to examine the mechanism of ATAA failure using a combination of detailed mechanical tests on human tissue and a multiscale computational model. We used (1) multiple, geometrically diverse, mechanical tests to characterize tissue properties; (2) a multiscale computational model to translate those results into a broadly usable form; and (3) a model-based computer simulation of the response of an ATAA to the stresses generated by the blood pressure. Mechanical tests were performed in uniaxial extension, biaxial extension, shear lap, and peel geometries. ATAA tissue was strongest in circumferential extension and weakest in shear, presumably because of the collagen and elastin in the arterial lamellae. A multiscale, fiber-based model using different fiber properties for collagen, elastin, and interlamellar connections was specified to match all of the experimental data with one parameter set. Finally, this model was used to simulate ATAA inflation using a realistic geometry. The predicted tissue failure occurred in regions of high stress, as expected; initial failure events involved almost entirely interlamellar connections, consistent with arterial dissection-the elastic lamellae remain intact, but the connections between them fail. The failure of the interlamellar connections, paired with the weakness of the tissue under shear loading, is suggestive that shear stress within the tissue may contribute to ATAA dissection.

Perspectives on Sharing Models and Related Resources in Computational Biomechanics Research.

The role of computational modeling for biomechanics research and related clinical care will be increasingly prominent. The biomechanics community has been developing computational models routinely for exploration of the mechanics and mechanobiology of diverse biological structures. As a result, a large array of models, data, and discipline-specific simulation software has emerged to support endeavors in computational biomechanics. Sharing computational models and related data and simulation software has first become a utilitarian interest, and now, it is a necessity. Exchange of models, in support of knowledge exchange provided by scholarly publishing, has important implications. Specifically, model sharing can facilitate assessment of reproducibility in computational biomechanics and can provide an opportunity for repurposing and reuse, and a venue for medical training. The community's desire to investigate biological and biomechanical phenomena crossing multiple systems, scales, and physical domains, also motivates sharing of modeling resources as blending of models developed by domain experts will be a required step for comprehensive simulation studies as well as the enhancement of their rigor and reproducibility. The goal of this paper is to understand current perspectives in the biomechanics community for the sharing of computational models and related resources. Opinions on opportunities, challenges, and pathways to model sharing, particularly as part of the scholarly publishing workflow, were sought. A group of journal editors and a handful of investigators active in computational biomechanics were approached to collect short opinion pieces as a part of a larger effort of the IEEE EMBS Computational Biology and the Physiome Technical Committee to address model reproducibility through publications. A synthesis of these opinion pieces indicates that the community recognizes the necessity and usefulness of model sharing. There is a strong will to facilitate model sharing, and there are corresponding initiatives by the scientific journals. Outside the publishing enterprise, infrastructure to facilitate model sharing in biomechanics exists, and simulation software developers are interested in accommodating the community's needs for sharing of modeling resources. Encouragement for the use of standardized markups, concerns related to quality assurance, acknowledgement of increased burden, and importance of stewardship of resources are noted. In the short-term, it is advisable that the community builds upon recent strategies and experiments with new pathways for continued demonstration of model sharing, its promotion, and its utility. Nonetheless, the need for a long-term strategy to unify approaches in sharing computational models and related resources is acknowledged. Development of a sustainable platform supported by a culture of open model sharing will likely evolve through continued and inclusive discussions bringing all stakeholders at the table, e.g., by possibly establishing a consortium.

Computational Parametric Analysis of the Mechanical Response of Structurally Varying Pacinian Corpuscles.

The Pacinian corpuscle (PC) is a cutaneous mechanoreceptor that senses low-amplitude, high-frequency vibrations. The PC contains a nerve fiber surrounded by alternating layers of solid lamellae and interlamellar fluid, and this structure is hypothesized to contribute to the PC's role as a band-pass filter for vibrations. In this study, we sought to evaluate the relationship between the PC's material and geometric parameters and its response to vibration. We used a spherical finite element mechanical model based on shell theory and lubrication theory to model the PC's outer core. Specifically, we analyzed the effect of the following structural properties on the PC's frequency sensitivity: lamellar modulus (E), lamellar thickness (h), fluid viscosity (µ), PC outer radius (Ro), and number of lamellae (N). The frequency of peak strain amplification (henceforth "peak frequency") and frequency range over which strain amplification occurred (henceforth "bandwidth") increased with lamellar modulus or lamellar thickness and decreased with an increase in fluid viscosity or radius. All five structural parameters were combined into expressions for the relationship between the parameters and peak frequency, ωpeak=1.605×10-6N3.475(Eh/µRo), or bandwidth, B=1.747×10-6N3.951(Eh/µRo). Although further work is needed to understand how mechanical variability contributes to functional variability in PCs and how factors such as PC eccentricity also affect PC behavior, this study provides two simple expressions that can be used to predict the impact of structural or material changes with aging or disease on the frequency response of the PC.

Collagen Organization in Facet Capsular Ligaments Varies With Spinal Region and With Ligament Deformation.

The spinal facet capsular ligament (FCL) is primarily comprised of heterogeneous arrangements of collagen fibers. This complex fibrous structure and its evolution under loading play a critical role in determining the mechanical behavior of the FCL. A lack of analytical tools to characterize the spatial anisotropy and heterogeneity of the FCL's microstructure has limited the current understanding of its structure-function relationships. Here, the collagen organization was characterized using spatial correlation analysis of the FCL's optically obtained fiber orientation field. FCLs from the cervical and lumbar spinal regions were characterized in terms of their structure, as was the reorganization of collagen in stretched cervical FCLs. Higher degrees of intra- and intersample heterogeneity were found in cervical FCLs than in lumbar specimens. In the cervical FCLs, heterogeneity was manifested in the form of curvy patterns formed by collections of collagen fibers or fiber bundles. Tensile stretch, a common injury mechanism for the cervical FCL, significantly increased the spatial correlation length in the stretch direction, indicating an elongation of the observed structural features. Finally, an affine estimation for the change of correlation length under loading was performed which gave predictions very similar to the actual values. These findings provide structural insights for multiscale mechanical analyses of the FCLs from various spinal regions and also suggest methods for quantitative characterization of complex tissue patterns.

Cellular Microbiaxial Stretching to Measure a Single-Cell Strain Energy Density Function.

The stress in a cell due to extracellular mechanical stimulus is determined by its mechanical properties, and the structural organization of many adherent cells suggests that their properties are anisotropic. This anisotropy may significantly influence the cells' mechanotransductive response to complex loads, and has important implications for development of accurate models of tissue biomechanics. Standard methods for measuring cellular mechanics report linear moduli that cannot capture large-deformation anisotropic properties, which in a continuum mechanics framework are best described by a strain energy density function (SED). In tissues, the SED is most robustly measured using biaxial testing. Here, we describe a cellular microbiaxial stretching (CµBS) method that modifies this tissue-scale approach to measure the anisotropic elastic behavior of individual vascular smooth muscle cells (VSMCs) with nativelike cytoarchitecture. Using CµBS, we reveal that VSMCs are highly anisotropic under large deformations. We then characterize a Holzapfel-Gasser-Ogden type SED for individual VSMCs and find that architecture-dependent properties of the cells can be robustly described using a formulation solely based on the organization of their actin cytoskeleton. These results suggest that cellular anisotropy should be considered when developing biomechanical models, and could play an important role in cellular mechano-adaptation.

Isotropic Failure Criteria Are Not Appropriate for Anisotropic Fibrous Biological Tissues.

The von Mises (VM) stress is a common stress measure for finite element models of tissue mechanics. The VM failure criterion, however, is inherently isotropic, and therefore may yield incorrect results for anisotropic tissues, and the relevance of the VM stress to anisotropic materials is not clear. We explored the application of a well-studied anisotropic failure criterion, the Tsai­Hill (TH) theory, to the mechanically anisotropic porcine aorta. Uniaxial dogbones were cut at different angles and stretched to failure. The tissue was anisotropic, with the circumferential failure stress nearly twice the axial (2.67 ± 0.67 MPa compared to 1.46 ± 0.59 MPa). The VM failure criterion did not capture the anisotropic tissue response, but the TH criterion fit the data well (R2 = 0.986). Shear lap samples were also tested to study the efficacy of each criterion in predicting tissue failure. Two-dimensional failure propagation simulations showed that the VM failure criterion did not capture the failure type, location, or propagation direction nearly as well as the TH criterion. Over the range of loading conditions and tissue geometries studied, we found that problematic results that arise when applying the VM failure criterion to an anisotropic tissue. In contrast, the TH failure criterion, though simplistic and clearly unable to capture all aspects of tissue failure, performed much better. Ultimately, isotropic failure criteria are not appropriate for anisotropic tissues, and the use of the VM stress as a metric of mechanical state should be reconsidered when dealing with anisotropic tissues.

Failure of the Porcine Ascending Aorta: Multidirectional Experiments and a Unifying Microstructural Model.

The ascending thoracic aorta is poorly understood mechanically, especially its risk of dissection. To make better predictions of dissection risk, more information about the multidimensional failure behavior of the tissue is needed, and this information must be incorporated into an appropriate theoretical/computational model. Toward the creation of such a model, uniaxial, equibiaxial, peel, and shear lap tests were performed on healthy porcine ascending aorta samples. Uniaxial and equibiaxial tests showed anisotropy with greater stiffness and strength in the circumferential direction. Shear lap tests showed catastrophic failure at shear stresses (150-200 kPa) much lower than uniaxial tests (750-2500 kPa), consistent with the low peel tension (∼60 mN/mm). A novel multiscale computational model, including both prefailure and failure mechanics of the aorta, was developed. The microstructural part of the model included contributions from a collagen-reinforced elastin sheet and interlamellar connections representing fibrillin and smooth muscle. Components were represented as nonlinear fibers that failed at a critical stretch. Multiscale simulations of the different experiments were performed, and the model, appropriately specified, agreed well with all experimental data, representing a uniquely complete structure-based description of aorta mechanics. In addition, our experiments and model demonstrate the very low strength of the aorta in radial shear, suggesting an important possible mechanism for aortic dissection.

A computational model of flow and species transport in the mesangium.

A variety of macromolecules accumulate in the glomerular mesangium in many different diseases, but the physics of the transport of these molecules within the mesangial matrix has not been extensively studied. We present a computational model of convection and diffusion within the porous mesangial matrix and apply this model to the specific instance of immunoglobulin A (IgA) transport in IgA nephropathy. We examine the influence of physiological factors including glomerular basement membrane (GBM) thickness and mesangial matrix density on the total accumulation of IgA. Our results suggest that IgA accumulation can be understood by relating convection and diffusion, thus demonstrating the importance of intrinsic glomerular factors.