Search

The aim of this study was to determine whether specific three-dimensional aortic shape features, extracted via statistical shape analysis (SSA), correlate with the development of thoracic ascending aortic dissection (TAAD) risk and associated aortic hemodynamics. Thirty-one patients followed prospectively with ascending thoracic aortic aneurysm (ATAA), who either did (12 patients) or did not (19 patients) develop TAAD, were included in the study, with aortic arch geometries extracted from computed tomographic angiography (CTA) imaging. Arch geometries were analyzed with SSA, and unsupervised and supervised (linked to dissection outcome) shape features were extracted with principal component analysis (PCA) and partial least squares discriminant analysis (PLS-DA), respectively. We determined PLS-DA to be effective at separating dissection and no-dissection patients (p=0.0010), with decreased tortuosity and more equal ascending and descending aortic diameters associated with higher dissection risk. In contrast, neither PCA nor traditional morphometric parameters (maximum diameter, tortuosity, or arch volume) were effective at separating dissection and no-dissection patients. The arch shapes associated with higher dissection probability were supported with hemodynamic insight. Computational fluid dynamics (CFD) simulations revealed a correlation between the PLS-DA shape features and wall shear stress (WSS), with higher maximum WSS in the ascending aorta associated with increased risk of dissection occurrence. Our work highlights the potential importance of incorporating higher dimensional geometric assessment of aortic arch anatomy in TAAD risk assessment, and in considering the interdependent influences of arch shape and hemodynamics as mechanistic contributors to TAAD occurrence.

BackgroundBleb presence in intracranial aneurysms (IAs) is a known indication of instability and vulnerability.ObjectiveTo develop and evaluate predictive models of bleb development in IAs based on hemodynamics, geometry, anatomical location, and patient population.MethodsCross-sectional data (one time point) of 2395 IAs were used for training bleb formation models using machine learning (random forest, support vector machine, logistic regression, k-nearest neighbor, and bagging). Aneurysm hemodynamics and geometry were characterized using image-based computational fluid dynamics. A separate dataset with 266 aneurysms was used for model evaluation. Model performance was quantified by the area under the receiving operating characteristic curve (AUC), true positive rate (TPR), false positive rate (FPR), precision, and balanced accuracy.ResultsThe final model retained 18 variables, including hemodynamic, geometrical, location, multiplicity, and morphology parameters, and patient population. Generally, strong and concentrated inflow jets, high speed, complex and unstable flow patterns, and concentrated, oscillatory, and heterogeneous wall shear stress patterns together with larger, more elongated, and more distorted shapes were associated with bleb formation. The best performance on the validation set was achieved by the random forest model (AUC=0.82, TPR=91%, FPR=36%, misclassification error=27%).ConclusionsBased on the premise that aneurysm characteristics prior to bleb formation resemble those derived from vascular reconstructions with their blebs virtually removed, machine learning models can identify aneurysms prone to bleb development with good accuracy. Pending further validation with longitudinal data, these models may prove valuable for assessing the propensity of IAs to progress to vulnerable states and potentially rupturing.

BackgroundBlebs are important secondary structures of intracranial aneurysms associated with increased rupture risk and can affect local wall stress and hemodynamics. Mechanisms of bleb development and evolution are not clearly understood. We investigate the relationship between blebs with different wall characteristics and local hemodynamics and rupture sites.MethodsBlebs with different wall appearances in intra-operative videos were analyzed with image-based computational fluid dynamics. Thin red blebs were compared against thick atherosclerotic/hyperplastic white/yellow blebs. Rupture points were identified in videos of ruptured aneurysms harboring blebs.ResultsThin blebs tended to be closer to the inflow than atherosclerotic blebs of the same aneurysm (P=0.0234). Blebs near the inflow had higher velocity (P=0.0213), vorticity (P=0.0057), shear strain rate (P=0.0084), wall shear stress (WSS) (P=0.0085), and WSS gradient (P=0.0151) than blebs far from the inflow. In a subset of 12 ruptured aneurysms harboring blebs, rupture points were associated with thin blebs in 42% of aneurysms, atherosclerotic blebs in 25%, and were away from blebs in the remaining 33%.ConclusionsNot all blebs are equal; some have thin translucent walls while others have thick atherosclerotic walls. Thin blebs tend to be located closer to the inflow than atherosclerotic blebs. Blebs near the inflow are exposed to stronger flows with higher and spatially variable WSS than blebs far from the inflow which tend to have uniformly lower WSS. Aneurysms can rupture at thin blebs, atherosclerotic blebs, and even away from blebs. Further study of wall failure in aneurysms with different bleb types is needed.

Amorphous silicon (a-Si) in the form of various nanoscale architectures (nanoparticles, nanoflakes, nanorods, nanowires, hollow nanotubes, and thin films) is a promising alternative anode to traditional carbon with considerable potential for high capacity lithium-ion batteries. However, silicon-based systems suffer from colossal volume expansion-related stresses leading to electrode decrepitation causing rapid capacity fade, and consequently, poor cycling performance. Loss of mechanical integrity of the anode due to lithium alloying-induced stresses is the primary reason for this capacity fade. We hypothesized that disintegration of the a-Si thin film occurs due to the mechanical constraints against lithium intercalation/alloying-induced volume expansion imposed by the current collector, and alteration of this mechanical constraint using an intervening soft layer will considerably improve the stability of the anode system. To test this hypothesis, we utilized an innovative custom-developed multiphysics finite element framework to model the anode system. Our simulations revealed that modulus mismatch between the silicon thin film and the elastic layer sandwich between this layer and the current collector is a key design parameter critical for improving the mechanical integrity of the silicon thin film. We also determined that the adhesion between the different anode components is important in prolonging the anode life.

The aim of this study was to determine whether specific three-dimensional aortic shape features, extracted via statistical shape analysis (SSA), correlate with the development of thoracic ascending aortic dissection (TAAD) risk and associated aortic hemodynamics. Thirty-one patients followed prospectively with ascending thoracic aortic aneurysm (ATAA), who either did (12 patients) or did not (19 patients) develop TAAD, were included in the study, with aortic arch geometries extracted from computed tomographic angiography (CTA) imaging. Arch geometries were analyzed with SSA, and unsupervised and supervised (linked to dissection outcome) shape features were extracted with principal component analysis (PCA) and partial least squares discriminant analysis (PLS-DA), respectively. We determined PLS-DA to be effective at separating dissection and no-dissection patients ([Formula: see text]), with decreased tortuosity and more equal ascending and descending aortic diameters associated with higher dissection risk. In contrast, neither PCA nor traditional morphometric parameters (maximum diameter, tortuosity, or arch volume) were effective at separating dissection and no-dissection patients. The arch shapes associated with higher dissection probability were supported with hemodynamic insight. Computational fluid dynamics (CFD) simulations revealed a correlation between the PLS-DA shape features and wall shear stress (WSS), with higher maximum WSS in the ascending aorta associated with increased risk of dissection occurrence. Our work highlights the potential importance of incorporating higher dimensional geometric assessment of aortic arch anatomy in TAAD risk assessment, and in considering the interdependent influences of arch shape and hemodynamics as mechanistic contributors to TAAD occurrence.

OBJECTIVE: To assess ascending aortic distensibility and build geometry and distensibility-based patient-specific stress distribution maps in patients sustaining type A aortic dissection (TAAD) using pre-dissection non-invasive imaging. METHODS: Review of charts from patients undergoing surgical repair of TAAD (n=351) led to the selection of a subset population (n=7) with two or more pre-dissection CTAs and echocardiograms at least one year prior to dissection. Ascending aortic wall biomechanical properties (aortic strain, distensibility and stiffness) were compared to age and size matched non-dissected non-aneurysmal controls. Patient-specific aortic strain served as an input in aortic geometry-based simulated three-dimensional reconstructions to generate longitudinal and circumferential wall stress maps. Inspection of peri-operative dissection scans and intraoperative visual examination confirmed primary tear locations. RESULTS: Pre-dissection echocardiography revealed ascending aortas of patients sustaining TAAD to exhibit decreased aortic wall strain (14.50 ± 1.13% vs 8.49 ± 1.08%; p < 0.01), decreased distensibility (4.26 ± 0.44 vs 2.39 ± 0.33 10(−6)cm(2)dyne(−1); p < 0.01), increased stiffness (3.84 ± 0.24 vs 7.48 ± 1.05; p < 0.001), and increased longitudinal wall stress (246 ± 22 vs 172 ± 37 kPa; p < 0.01). There was no significant difference in circumferential wall stress. Pre-dissection CTA models revealed overlap between regions of increased longitudinal wall stress and primary tear sites. CONCLUSIONS: Using pre-dissection imaging, we identified increased stiffness and longitudinal wall stress in ascending aortas of dissection patients. Patient-specific imaging-derived biomechanical property maps like these may be instrumental toward designing better prediction models of aortic dissection potential.

Biomechanical failure of the artery wall can lead to rupture, a catastrophic event with a high rate of mortality. Thus, there is a pressing need to understand failure behavior of the arterial wall. Uniaxial testing remains the most common experimental technique to assess tissue failure properties. However, the relationship between intrinsic failure parameters of the tissue and measured uniaxial failure properties is not fully established. Furthermore, the effect of the experimental variables, such as specimen shape and boundary conditions, on the measured failure properties is not well understood. We developed a finite element model capable of recapitulating pre-failure and post-failure uniaxial biomechanical response of the arterial tissue specimen. Intrinsic stiffness, strength and fracture toughness of the vessel wall tissue were used as the input material parameters to the model. Two uniaxial testing protocols were considered: a conventional setup with a rectangular specimen held at the grips by cardboard inserts, and the other used a dogbone specimen with soft foam inserts at the grips. Our computational study indicated negligible differences in the peak stress and post-peak mechanical behavior between these two testing protocols. It was also found that the tissue experienced only modest localized failure until higher levels of applied stretch beyond the peak stress. A robust cohesive model was capable of modeling the post-peak biomechanical response, which was primarily governed by tissue fracture toughness. Our results suggest that the post-peak region, in conjunction with the peak stress, must be considered to evaluate the complete biomechanical failure behavior of the soft tissue.

BACKGROUND: Calcification was recently found to be present in the majority of cerebral aneurysms, though how calcification and the presence or absence of co-localized lipid pools affect failure properties is still unknown. OBJECTIVE: The primary objective is to quantify the biomechanical effect of a macro-calcification with surrounding Near-Calcification Region (NCR) of varying mechanical properties on tissue failure behavior. METHODS: We utilized a structurally informed finite element model to simulate pre-failure and failure behavior of a human cerebral tissue specimen modeled as a composite containing a macro-calcification and surrounding NCR, embedded in a fiber matrix composite. Data from multiple imaging modalities was combined to quantify the collagen organization and calcification geometry. An idealized parametric model utilizing the calibrated model was used to explore the impact of NCR properties on tissue failure. RESULTS: Compared to tissue without calcification, peak stress was reduced by 82% and 49% for low modulus (representing lipid pool) and high modulus (simulating increase in calcification size) of the NCR, respectively. Failure process strongly depended on NCR properties with lipid pools blunting the onset of complete failure. When the NCR was calcified, the sample was able to sustain larger overall stress, however the failure process was abrupt with nearly simultaneous failure of the loaded fibers. CONCLUSIONS: Failure of calcified vascular tissue is strongly influenced by the ultrastructure in the vicinity of the calcification. Computational modeling of failure in fibrous soft tissues can be used to understand how pathological changes impact the tissue failure process, with potentially important clinical implications.

Right ventricular (RV) remodeling and longitudinal fiber reorientation in the setting of pulmonary hypertension (PH) affects ventricular structure and function, eventually leading to RV failure. Characterizing the kinematics of myocardial fibers helps better understanding the underlying mechanisms of fiber realignment in PH. In the current work, high-frequency ultrasound imaging and structurally-informed finite element (FE) models were employed for an exploratory evaluation of the stretch-induced kinematics of RV fibers. Image-based experimental evaluation of fiber kinematics in porcine myocardium revealed the capability of affine assumptions to effectively approximate myofiber realignment in the RV free wall. The developed imaging framework provides a noninvasive modality to quantify transmural RV myofiber kinematics in large animal models. FE modeling results demonstrated that chronic pressure overload, but not solely an acute rise in pressures, results in kinematic shift of RV fibers towards the longitudinal direction. Additionally, FE simulations suggest a potential protective role for concentric hypertrophy (increased wall thickness) against fiber reorientation, while eccentric hypertrophy (RV dilation) resulted in longitudinal fiber realignment. Our study improves the current understanding of the role of different remodeling events involved in transmural myofiber reorientation in PH. Future experimentations are warranted to test the model-generated hypotheses.

BackgroundBlebs are rupture risk factors in intracranial aneurysms (IAs), but their prevalence, distribution, and associations with clinical factors as well as their causes and effects on aneurysm vulnerability remain unclear.MethodsA total of 122 blebs in 270 IAs selected for surgery were studied using patient-specific vascular reconstructions from 3D angiographic images. Bleb geometry, location on the aneurysm, and frequency of occurrence in aneurysms at different locations were analyzed. Associations between gender, age, smoking, hypertension, hormone therapy, dental infection, and presence of blebs were investigated.ResultsOf all aneurysms with blebs, 77% had a single bleb and 23% had multiple blebs. Only 6% of blebs were at the neck, while 46% were in the body and 48% in the dome. Aneurysms with blebs were larger (pConclusionsBlebs are common in IAs, and most aneurysms harboring blebs have a single bleb. Blebs in the aneurysm neck are rare, but they are equally common in the body and dome. The presence of blebs in IAs was associated with dental infection, and negatively associated with hormone replacement therapy.

Rotator cuff tear propagation is multifactorial and may be due to localized changes in mechanical properties from tendon remodeling based on the inhomogeneous stresses experienced by a tendon with a tear. The objective of this study was to investigate the effect of localized tendon remodeling on tear propagation for simulated supraspinatus tendon tears. A validated computational model of a supraspinatus tendon using subject-specific geometry and material properties with a 1 cm wide anterior tear was used. The medial edge of the supraspinatus tendon was displaced 5 mm to induce tear propagation and cohesive elements were used to model tear propagation. Four remodeling scenarios were investigated: (1) Baseline (no remodeling), (2) Positive remodeling (increased fiber stiffness) and (3) Negative remodeling (decreased fiber stiffness) at tear tips, and (4) Negative remodeling along the medial-lateral tear edge. Output parameters included the amount of tear propagation, critical load to propagate the tear, and maximum principal stress at the tear tips. Positive remodeling at the tear tips resulted in the largest amount of tear propagation (18.4 mm), highest peak maximum principal stress (25.2 MPa), and lowest critical load to propagate the tear (249N). Conversely, negative remodeling at the tear tips resulted in the least amount of tear propagation (16 mm), lowest peak maximum principal stress (17.6 MPa) and highest critical load to propagate the tear (278N). Overall, remodeling at the tear tips has the greatest effect on tear propagation. Therefore, a better method for clinicians to measure tendon stiffness at the tear tips would be helpful to improve outcome of patients.

Rotator cuff tendons undergo degeneration with age, which could have an impact on tear propagation. The objective of this study was to predict tear propagation for different levels of tissue degeneration using an experimentally validated finite element model of a supraspinatus tendon. It was hypothesized that greater amounts of degeneration will result in tear propagation at lower loads than tendons with less degeneration. Using a previously-validated computational model of supraspinatus tendon, 1-cm tears were introduced in the anterior, middle, and posterior thirds of the tendon. Cohesive elements were assigned subject-specific failure properties to model tear propagation, and tendon degeneration ranging from "minimal" to "severe" was modeled by modifying its mechanical properties. Tears in tendons with severe degeneration required the smallest loads to propagate (122-207 N). Posterior tears required greater loads compared to middle and anterior tears at all levels of degeneration. Stress and strain required for tear propagation decreased substantially with degeneration, ranging from 8.5 MPa and 32.6% strain for minimal degeneration and 0.6 MPa and 4.5% strain for severe degeneration. Overall, this work indicates that greater amounts of tendon degeneration lead to greater risk of tear propagation, supporting the need for early detection and treatment of rotator cuff tears.

High lethality of aortic dissection necessitates accurate predictive metrics for dissection risk assessment. The not infrequent incidence of dissection at aortic diameters

Rotator cuff tears remain a significant clinical problem with a high incidence rate and severe clinical burden. Previous computational models developed to study rotator cuff tears have not modeled tissue damage and tear propagation. The objective of this study was to predict tear propagation for various combinations of tear size and location using an experimentally validated finite element model of supraspinatus tendon. It was hypothesized that larger rotator cuff tears propagate at lower loads than smaller tears, and that posterior tears require higher loads to propagate than anterior tears. Using a previously validated computational model of supraspinatus tendon, tears of size 0.5–1.5 cm were introduced to the tendon geometry in the anterior, middle, and posterior tendon thirds. Cohesive elements were assigned subject-specific failure properties and used to model tissue damage and tear propagation. A displacement of 5 mm was applied to the medial tendon edge to induce tear propagation. Model outputs included critical load required to propagate the tear, and principal stress and maximum principal strain at the anterior and posterior tear tips. For all tear sizes, posterior tears required the highest loads to propagate (247–567 N). Anterior tears generally required the least load to propagate (171–280 N). Stress and strain were larger on the articular side (maximum 33.9% articular strain vs 27.8% bursal strain). Overall, larger tears located in the anterior supraspinatus tendon that interrupt the rotator cable are most at risk for tear propagation, and should be carefully followed by clinicians when considering treatment options.

We examine the stretching behavior of rubber-plastic composites composed of a layer of styrene-ethylene/propylene-styrene (SEPS) rubber, bonded to a layer of linear low density polyethylene (LLDPE) plastic. Dog-bone shaped samples of rubber, plastic, and rubber-plastic bilayers with rubber : plastic thickness ratio in the range of 1.2-9 were subjected to uniaxial tension tests. The degree of inhomogeneity of deformation was quantified by digital image correlation analysis of video recordings of these tests. In tension, the SEPS layer showed homogeneous deformation, whereas the LLDPE layer showed necking followed by stable drawing owing to its elastoplastic deformation behavior and post-yield strain hardening. Bilayer laminates showed behavior intermediate between the plastic and the rubber, with the degree of necking and drawing reducing as the rubber : plastic ratio increased. A simple model was developed in which the force in the bilayer was taken as the sum of forces in the plastic and the rubber layers measured independently. By applying a mechanical energy balance to this model, the changes in bilayer necking behavior with rubber thickness could be predicted qualitatively.

Polymeric plastics deform irreversibly (i.e., inelastically) whereas rubbers deform reversibly, i.e., elastically. Thus, uniaxially stretching a rubber-plastic bilayer composite beyond its yield point can create an elastic strain mismatch between the two layers. Upon release, the bilayer may then bend out-of-plane. We quantify the mechanics of such stretch-release-induced shape changes in rectangular specimens of rubber-plastic bilayers. We show a remarkable dependence of the final shape upon the stretch applied prior to release. At small stretch, all bilayers bend into arch or roll shapes with the plastic on the convex face. At large stretch, the bilayers bend into half-tubes with the plastic, now heavily-wrinkled, becoming the concave face. Thus, the sign and direction of the curvature both flip as applied stretch increases. Between these two extremes, saddle shapes appear which have characteristics of both arches as well as half-tubes. Sufficiently narrow samples show different behavior: they transition from arches to helices as stretch increases. All these shapes are mono-stable. We document numerous ways in which the mechanics of rubber-plastic bilayers differs from that of fully-elastic bilayers. Most importantly, yielding of the plastic layer during the shape change strongly affects the mechanics of the elastic–plastic bilayers, and yielding accompanied by plastic wrinkling has an especially large effect. A strain energy model illustrates how the shape change is dictated by the change in the ratio of elastic strain mismatch in the two directions due to the formation of wrinkles at the rubber-plastic interface.

Rotator cuff tears are a significant clinical problem previously investigated by unvalidated computational models that either use simplified geometry or isotropic elastic material properties to represent the tendon. The objective of this study was to develop an experimentally validated, finite element model of supraspinatus tendon using specimen-specific geometry and inhomogeneous material properties to predict strains in intact supraspinatus tendon at multiple abduction angles. Three-dimensional tendon surface strains were determined at 60 deg, 70 deg, and 90 deg of glenohumeral abduction for articular and bursal surfaces of supraspinatus tendon during cyclic loading (5–200 N, 50 cycles, 20 mm/min) to serve as validation data for computational model predictions. A finite element model was developed using the tendon geometry and inhomogeneous material properties to predict surface strains for loading conditions mimicking experimental loading conditions. Experimental strains were directly compared with computational model predictions to validate the model. Overall, the model successfully predicted magnitudes of strains that were within the experimental repeatability of 3% strain of experimental measures on both surfaces of the tendon. Model predictions and experiments showed the largest strains to be located on the articular surface (∼8% strain) between the middle and the anterior edge of the tendon. Importantly, the reference configuration chosen to calculate strains had a significant effect on strain calculations, and therefore, must be defined with an innovative optimization algorithm. This study establishes a rigorously validated specimen-specific (both geometry and material properties) computational model using novel surface strain measurements for the use in investigating the function of the supraspinatus tendon and to ultimately predict the propagation of supraspinatus tendon tears based on the tendon's mechanical environment.

Vaginal tears are very common and can lead to severe complications such as hemorrhaging, fecal incontinence, urinary incontinence, and dyspareunia. Despite the implications of vaginal tears on women's health, there are currently no experimental studies on the tear behavior of vaginal tissue. In this study, planar equi-biaxial tests on square specimens of vaginal tissue, with sides oriented along the longitudinal direction (LD) and circumferential direction (CD), were conducted using swine as animal model. Three groups of specimens were mechanically tested: the NT group (n = 9), which had no pre-imposed tear, the longitudinal tear (LT) group (n = 9), and the circumferential tear (CT) group (n = 9), which had central pre-imposed elliptically shaped tears with major axes oriented in the LD and the CD, respectively. Through video recording during testing, axial strains were measured for the NT group using the digital image correlation (DIC) technique and axial displacements of hook clamps were measured for the NT, LT, and CT groups in the LD and CD. The swine vaginal tissue was found to be highly nonlinear and somewhat anisotropic. Up to normalized axial hook displacements of 1.15, no tears were observed to propagate, suggesting that the vagina has a high resistance to further tearing once a tear has occurred. However, in response to biaxial loading, the size of the tears for the CT group increased significantly more than the size of the tears for the LT group (p = 0.003). The microstructural organization of the vagina is likely the culprit for its tear resistance and orientation-dependent tear behavior. Further knowledge on the structure–function relationship of the vagina is needed to guide the development of new methods for preventing the severe complications of tearing.

Background: This study's aim was to determine whether the inclusion of superficial perineal structures in a finite-element simulation of vaginal delivery impacts the pubovisceral muscle and perineal body, two common sites of birth-related injury. The hypothesis, inferred from prevailing literature, was that these structures would have minimal influence (differences less than ±10%). Methods: Two models were made using the Visible Human Project's female cadaver to create a rigid, fixed pelvis, musculature held by spring attachments to that pelvis, and a rigid, ellipsoidal fetal head prescribed with an inferior displacement to simulate delivery. Injury site stretch ratios and fetal head and perineal body displacements and angles of progression were compared between the Omitted Model (which excluded the superficial perineal structures as is common practice) and the Included Model (which included them). Results: Included Model stretch ratios were +107%, −9.84% and −14.6% compared to Omitted Model perineal body and right and left pubovisceral muscles, respectively. Included Model peak perineal body inferior displacement was +72.5% greater while similar anterior–posterior displacements took longer to reach. Conclusion: These results refute our hypothesis, suggesting superficial perineal structures impact simulations of vaginal delivery by inhibiting perineal body anterior–posterior displacement, which stretches and inferiorly displaces the perineal body.

Incorporation of collagen structural information into the study of biomechanical behavior of ascending thoracic aortic (ATA) wall tissue should provide better insight into the pathophysiology of ATA. Structurally motivated constitutive models that include fiber dispersion and recruitment can successfully capture overall mechanical response of the arterial wall tissue. However, these models cannot examine local microarchitectural features of the collagen network, such as the effect of fiber disruptions and interaction between fibrous and non-fibrous components, which may influence emergent biomechanical properties of the tissue. Motivated by this need, we developed a finite element based three-dimensional structural model of the lamellar units of the ATA media that directly incorporates the collagen fiber microarchitecture. The fiber architecture was computer generated utilizing network features, namely fiber orientation distribution, intersection density and areal concentration, obtained from image analysis of multiphoton microscopy images taken from human aneurysmal ascending thoracic aortic media specimens with bicuspid aortic valve (BAV) phenotype. Our model reproduces the typical J-shaped constitutive response of the aortic wall tissue. We found that the stress state in the non-fibrous matrix was homogeneous until the collagen fibers were recruited, but became highly heterogeneous after that event. The degree of heterogeneity was dependent upon local network architecture with high stresses observed near disrupted fibers. The magnitude of non-fibrous matrix stress at higher stretch levels was negatively correlated with local fiber density. The localized stress concentrations, elucidated by this model, may be a factor in the degenerative changes in aneurysmal ATA tissue.

Heterostructures of silicon and carbon nanotubes (CNT) have been widely studied as Li-ion battery anodes. The focus of the current study is to investigate the role of silicon configurations on the mechanical integrity of the Si–CNT heterostructured anodes during electrochemical cycling. We hypothesize that void nucleation and growth in silicon during electrochemical cycling of Li can induce fracture and eventual failure. To test this hypothesis, we utilized a custom developed multiphysics finite element modeling framework considering the lithium diffusion induced elasto-plastic deformation of silicon. We systematically varied the silicon component configuration and enumerated the stress field within it for one complete electrochemical cycle. Resulting evolution of stress state reveals that reducing the mechanical constraints on Si reduces the plastic flow of the material, and thus possibility of void nucleation and growth. We find that the Si droplet configuration is mechanically stable while the continuous Si coating configuration is prone to void growth induced mechanical failure. Present analysis provides a mechanistic understanding of the effect of Si configurations in heterostructured electrodes on its mechanical integrity, which can help in design of next-generation hetersostructured electrodes with improved capacity retention.

Many plastics show necking and drawing behavior in tension, sometimes called “cold drawing”. In contrast, elastomers stretch homogeneously in tension. We examine the tensile behavior of rubber–plastic laminate composites using 3D finite element simulations and an analytical model. A rate-independent constitutive behavior was adopted in which the modulus at small-strain, strain hardening at large strain, and yield stress (only for the plastic) can all be varied independently. For sufficiently small rubber/plastic thickness ratio, layered composites show necking and drawing wherein a tensile bar coexists in two strain states, one with a large stretch (necked region) and the other with a modest stretch (unnecked region). With increasing rubber/plastic thickness ratio, the two strain states approach each other in a manner resembling a second order phase transition culminating in a critical point. Above this critical rubber/plastic thickness ratio, the layered composites stretch homogeneously. An analytical model based on adding the First Piola–Kirchoff stresses of the rubber and plastic layers, along with a modification for inelastic deformation, is shown to capture most of the results of 3D simulations accurately. We comment on the practical relevance of these results to toughening relatively brittle plastics, and more specifically, the critical importance of strain hardening of the rubber.

Type A Aortic Dissection (TAAD) is a life-threatening condition involving delamination of ascending aortic media layers. While current clinical guidelines recommend surgical intervention for aneurysm diameter > 5.5 cm, high incidence of TAAD in patients below this diameter threshold indicates the pressing need for improved evidence-based risk prediction metrics. Construction of such metrics will require the knowledge of the biomechanical failure properties of the aortic wall tissue under biaxial loading conditions. We utilized a fiber-level finite element based structural model of the aortic tissue to quantify the relationship between aortic tissue strength and physiologically relevant biaxial stress state for nonaneurysmal and aneurysmal patient cohorts with tricuspid aortic valve phenotype. We found that the model predicted strength of the aortic tissue under physiologic biaxial loading conditions depends on the stress biaxiality ratio, defined by the ratio of the longitudinal and circumferential components of the tissue stress. We determined that predicted biaxial tissue strength is statistically similar to its uniaxial circumferential strength below biaxiality ratios of 0.68 and 0.69 for nonaneurysmal and aneurysmal cohorts, respectively. Beyond this biaxiality ratio, predicted biaxial strength for both cohorts reduced drastically to a magnitude statistically similar to its longitudinal strength. We identified fiber-level failure mechanisms operative under biaxial stress state governing aforementioned tissue failure behavior. These findings are an important first step towards the development of mechanism-based TAAD risk assessment metrics for early identification of high-risk patients.

Although uniaxial tensile testing is commonly used to evaluate failure properties of vascular tissue, there is no established protocol for specimen shape or gripping method. Large percentages of specimens are reported to fail near the clamp and can potentially confound the studies, or, if discarded will result in sample waste. The objective of this study is to identify sample geometry and clamping conditions that can achieve consistent failure in the midregion of small arterial specimens, even for vessels from older individuals. Failure location was assessed in 17 dogbone specimens from human cerebral and sheep carotid arteries using soft inserts. For comparison with commonly used protocols, an additional 22 rectangular samples were tested using either sandpaper or foam tape inserts. Midsample failure was achieved in 94% of the dogbone specimens, while only 14% of the rectangular samples failed in the midregion, the other 86% failing close to the clamps. Additionally, we found midregion failure was more likely to be abrupt, caused by cracking or necking. In contrast, clamp failure was more likely to be gradual and included a delamination mode not seen in midregion failure. Hence, this work provides an approach that can be used to obtain consistent midspecimen failure, avoiding confounding clamp-related artifacts. Furthermore, with consistent midregion failure, studies can be designed to image the failure process in small vascular samples providing valuable quantitative information about changes to collagen and elastin structure during the failure process.

Exposure following pelvic organ prolapse repair has been observationally associated with wrinkling of the implanted mesh. The purpose of this study was to quantify the impact of variable boundary conditions on the out-of-plane deformations of mesh subjected to tensile loading. Using photogrammetry and surface curvature analyses, deformed geometries were accessed for two commercially available products. Relative to standard clamping methods, the amount of out-of-plane deformation significantly increased when point loads were introduced to simulate suture fixation in-vivo. These data support the hypothesis that regional increases in the concentration of mesh potentially enhance the host’s foreign body response, leading to exposure.

Each year, millions of Americans suffer bone fractures, often requiring internal fixation. Current devices, like plates and screws, are made with permanent metals or resorbable polymers. Permanent metals provide strength and biocompatibility, but cause long-term complications and may require removal. Resorbable polymers reduce long-term complications, but are unsuitable for many load-bearing applications. To mitigate complications, degradable magnesium (Mg) alloys are being developed for craniofacial and orthopedic applications. Their combination of strength and degradation make them ideal for bone fixation. Previously, we conducted a pilot study comparing Mg and titanium devices with a rabbit ulna fracture model. We observed Mg device degradation, with uninhibited healing. Interestingly, we observed bone formation around degrading Mg, but not titanium, devices. These results highlighted the potential for these fixation devices. To better assess their efficacy, we conducted a more thorough study assessing 99.9% Mg devices in a similar rabbit ulna fracture model. Device degradation, fracture healing, and bone formation were evaluated using microcomputed tomography, histology and biomechanical tests. We observed device degradation throughout, and calculated a corrosion rate of 0.40 ± 0.04 mm/year after 8 weeks. In addition, we observed fracture healing by 8 weeks, and maturation after 16 weeks. In accordance with our pilot study, we observed bone formation surrounding Mg devices, with complete overgrowth by 16 weeks. Bend tests revealed no difference in flexural load of healed ulnae with Mg devices compared to intact ulnae. These data suggest that Mg devices provide stabilization to facilitate healing, while degrading and stimulating new bone formation.

Nerve blocks are frequently performed by anesthesiologists to control pain. For sciatic nerve blocks, the optimal placement of the needle tip between its paraneural sheath and epineurial covering is challenging, even under ultrasound guidance, and frequently results in nerve puncture. We performed needle penetration tests on cadaveric isolated paraneural sheath (IPS), isolated nerve (IN), and the nerve with overlying paraneural sheath (NPS), and quantified puncture force requirement and fracture toughness of these specimens to assess their role in determining the clinical risk of nerve puncture. We found that puncture force (123 ± 17 mN) and fracture toughness (45.48 ± 9.72 J m−2) of IPS was significantly lower than those for NPS (1440 ± 161 mN and 1317.46 ± 212.45 Jm−2, respectively), suggesting that it is not possible to push the tip of the block needle through the paraneural sheath only, without pushing it into the nerve directly, when the sheath is lying directly over the nerve. Results of this study provide a physical basis for tangential placement of the needle as the ideal situation for local anesthetic deposition, as it allows for the penetration of the sheath along the edge of the nerve without entering the epineurium.

We examine the mechanics of three-layer composite films composed of an elastomeric layer sandwiched between two thin surface layers of plastic. Upon stretching and releasing such composite films, they develop a highly wrinkled surface texture. The mechanism for this texturing is that during stretching, the plastic layers yield and stretch irreversibly whereas the elastomer stretches reversibly. Thus upon releasing, the plastic layers buckle due to compressive stress imposed by the elastomer. Experiments are conducted using SEPS elastomer and 50 micron thick LLDPE plastic films. Stretching and releasing the composites to 2-5 times their original length induces buckles with wavelength on the order of 200 microns, and the wavelength decreases as the stretching increases. FEM simulations reveal that plastic deformation is involved at all stages during this process: (1) during stretching, the plastic layer yields in tension; (2) during recovery, the plastic layer first yields in-plane in compression and then buckles; (3) post-buckling, plastic hinges are formed at high-curvature regions. Homogeneous wrinkles are predicted only within a finite window of material properties: if the yield stress is too low, the plastic layers yield in-plane, without wrinkling, whereas if the yield stress is too high, non-homogeneous wrinkles are predicted. This approach to realizing highly wrinkled textures offers several advantages, most importantly the fact that high aspect ratio wrinkles (amplitude to wavelength ratios exceeding 0.4) can be realized.

Sputter-deposited amorphous silicon thin films on metallic copper current collectors are widely studied as lithium-ion anode systems. Electrochemical results indicate these electrodes exhibit near theoretical capacity for first few cycles; however delamination at the thin film–current collector interface causes rapid capacity fade leading to poor cycling performance. Primary reason for this interfacial delamination is the mechanical stress generated due to colossal volume expansion of silicon during lithiation. The focus of the current study is to present a mechanistic understanding of the role of mechanical properties of the current collector on this characteristic delamination behavior during electrochemical cycling. Toward this end, we have developed a computational framework that accounts for the coupled diffusion induced large deformation in silicon, elasto-plastic deformation of the current collector, as well as the nucleation and propagation of interfacial delamination. We have also performed a detailed parametric study to investigate the effect of mechanical properties of the current collector on the delamination of the thin film–current collector interface. We have accordingly determined that current collectors with low elastic modulus such as graphite can completely suppress interfacial delamination. Our analysis thus provides a sound mechanistic approach for designing next generation Si thin film anodes with improved capacity retention.

We recently reported a mechanistic model to link micro-architectural information to the delamination strength (Sd) of human ascending thoracic aorta (ATA). That analysis demonstrated that the number density (N) and failure energy (Uf) of the radially-oriented collagen fibers contribute to the Sd of both aneurysmal (ATAA) and non-aneurysmal (CTRL-ATA) aortic tissue. Among the set of ATAA samples, we studied specimens from patients displaying bicuspid (BAV) and tricuspid aortic valve (TAV) morphologic phenotypes. Results from our prior work were based on the assumption that the Uf was independent of dissection direction. In the current study, we excluded that assumption and hypothesized that Uf correlates with the Sd of ATAA. To test the hypothesis, we used previously-reported experimentally-determined Sd measurements and N of radially-oriented collagen fibers as input in our validated mechanistic model to calculate Uf for BAV-ATAA, TAV-ATAA and CTRL-ATA tissue specimens. The results of our analysis revealed that Uf is significantly lower for both BAV-ATAA and TAV-ATAA compared to CTRL-ATA cases, and does not differ between BAV-ATAA and TAV-ATAA. Furthermore, we found that Uf is consistent between circumferential-radial and longitudinal-radial planes in either of BAV-ATAA, TAV-ATAA or CTRL-ATA specimens. These findings employ a novel mechanistic model to increase our understanding of the putative interrelationship between biomechanical properties, extracellular matrix biology, and failure energy of aortic dissection.

This study's aim was to determine whether the inclusion of superficial perineal structures in a finite-element simulation of vaginal delivery impacts the pubovisceral muscle and perineal body, two common sites of birth-related injury. The hypothesis, inferred from prevailing literature, was that these structures would have minimal influence (differences less than ±10%). Two models were made using the Visible Human Project's female cadaver to create a rigid, fixed pelvis, musculature held by spring attachments to that pelvis, and a rigid, ellipsoidal fetal head prescribed with an inferior displacement to simulate delivery. Injury site stretch ratios and fetal head and perineal body displacements and angles of progression were compared between the Omitted Model (which excluded the superficial perineal structures as is common practice) and the Included Model (which included them). Included Model stretch ratios were +107%, −9.84% and −14.6% compared to Omitted Model perineal body and right and left pubovisceral muscles, respectively. Included Model peak perineal body inferior displacement was +72.5% greater while similar anterior–posterior displacements took longer to reach. These results refute our hypothesis, suggesting superficial perineal structures impact simulations of vaginal delivery by inhibiting perineal body anterior–posterior displacement, which stretches and inferiorly displaces the perineal body.

Amorphous silicon thin film deposited on copper current collector is a promising candidate for the high capacity lithium-ion battery (LIB) anode. However these systems exhibit a rapid capacity fade, and as a result, poor cyclic performance. Interfacial delamination due to lithium intercalation induced stress is the primary reason behind this capacity fade. For the case of a clean interface, crack blunting resulting from the plastic flow of the metallic substrate will eventually arrest the delamination. However, experimental observations suggest a complete delamination of the thin film from the substrate indicating the existence of interface embrittling mechanisms. Further experiments have revealed the evidence of segregation of lithium into the interfacial region that can potentially embrittle the interface. The focus of the current study is to investigate the role of such irreversible mechanisms at the interface on its delamination response during electrochemical cycling. Towards this end, we have developed a computational framework that accounts for coupled diffusion induced large deformation in silicon thin film, elasto-plastic deformation of the copper current collector, as well as the nucleation and propagation of interfacial delamination. A detailed parametric study has been presented to investigate the effect of segregation induced embrittlement on the delamination growth. Present analysis provides a mechanistic understanding of the delamination behavior of the interface between the silicon thin film and copper current collector that may help in the design of future thin film based LIB anodes with improved cyclic performance.

Large volume change and associated stress generation is known to cause failure of the silicon thin film anode used for Lithium-ion batteries after a few cycles. Experimental observations suggest that plastic deformation of the underlying Cu substrate and degradation of the active/inactive interface are the primary reasons responsible for the capacity fade. The goal of the present study is to examine the interplay between these mechanisms using a computational mechanics approach. In the present study, a novel multi-physics finite element framework has been developed to simulate the lithiation and de-lithiation induced failure of amorphous Si (a-Si) thin film on Cu foil. The numerical framework is based on the finite deformation of the active silicon wherein diffusion of lithium occurs, plastic deformation of the Cu foil, and debonding of the active/inactive interface. The effect of substrate property, interfacial energy and kinetics of interface degradation has been examined.

Computational cost is a major issue in modern large-scale simulations used across different disciplines of science and engineering. Computationally efficient surrogate models that can represent the original model with desired accuracy have been explored in the recent past. However, with the exception of few efforts, most of these techniques rely on a reduced order representation of the original complex model, resulting in a loss of information. In this paper we demonstrate the applicability of high dimensional model representation (HDMR) technique in addressing this issue while preserving the original model dimension. We will discuss the applicability of this surrogate modeling technique in the field of feasibility analysis drawing examples from process systems and materials design. It will be shown that the original physical models can be essentially considered as a black box, and same methodology can be applied across all the examples studied. It is found that the accuracy of the surrogate models depends on the order of the approximation and number of sampling points employed. While first-order approximation is largely inadequate, second-order approximation is sufficient for the model systems studied. Sampling requirement is also dramatically low for the construction of these surrogate models.

Silicon dioxide (thickness 350 nm and 969 nm) and silicon nitride (thickness 218 nm) films deposited on silicon substrate using plasma enhanced chemical vapor deposition process were investigated using a Berkovich nanoindenter. The load-depth measurements revealed that the oxide films have lower modulus and hardness compared to the silicon substrate, where as the nitride film has a higher hardness and slightly lower modulus than the substrate. To delineate the substrate effect, a phenomenological model, that captures most of the ‘continuous stiffness measurement’ data, was proposed and then extended on both sides to determine the film and substrate properties. The modulus and hardness of the oxide film were around 53 GPa and 4–8 GPa where as those of the nitride film were around 150 GPa and 19 GPa, respectively. These values compare well with the measurements reported elsewhere in the literature.

A computational method for arbitrary crack motion through a finite element mesh, termed as the generalized cohesive element technique, is presented. In this method, an element with an internal discontinuity is replaced by two superimposed elements with a combination of original and imaginary nodes. Conventional cohesive zone modeling, limited to crack propagation along the edges of the elements, is extended to incorporate the intra-element mixed-mode crack propagation. Proposed numerical technique has been shown to be quite accurate, robust and mesh insensitive provided the cohesive zone ahead of the crack tip is resolved adequately. A series of numerical examples is presented to demonstrate the validity and applicability of the proposed method.

Single-walled carbon nanotube (SWNT)/poly(methyl methacrylate) (PMMA) composites were prepared using coagulation method. The electrical conductivity and the electromagnetic interference (EMI) shielding of SWNT/PMMA composites over the X-band (8–12 GHz) and the microwave (200–2000 MHz) frequency range have been investigated. The electrical conductivity of composites increases with SWNT loading by 13 orders of magnitude, from 10−15 to 10−2 Ω−1 cm−1 with a percolation threshold of about 3 wt% SWNTs. The effect of the sample thickness on the shielding effectiveness has been studied, and correlated to the electrical conductivity of composites. The data suggest that SWNT/PMMA composites containing higher SWNT loading (above 10 wt%) be useful for EMI shielding and those with lower SWNT loading be useful for electrostatic charge dissipation. The dominant shielding mechanism of SWNT/PMMA composites was also discussed. POLYM. ENG. SCI., 2009. © 2009 Society of Plastics Engineers

In this feature article, dynamic indentation-induced deformation and fracture phenomena in brittle solids are related to the deformation behavior and fracture patterns evolved during dynamic impact on structural ceramics. Examples of indentation-induced localized amorphization in boron carbide (B4C) and differences in fracture characteristics due to static and dynamic indentations on B4C and zirconium diboride–silicon carbide (ZrB2–SiC) ceramic composite are presented. These features are compared with the fracture patterns evolved during the high-velocity ball impact on SiC. The influences of processing-induced (inherent) flaws and deformation-induced flaws on fracture initiation in a brittle material are discussed. Recent developments in analytical modeling of indentation- and impact-induced damage in ceramics and computational efforts on fracture and fragmentation response of ceramics are thoroughly reviewed. Finally, some future research directions are identified.

Shear band spacing in Zr-based bulk metallic glasses (BMGs) under dynamic loads is found to vary with position and local strain rate in the indented region. To investigate the dependence of shear band evolution characteristics on local strain rate and normal stress, a micromechanical model based on momentum diffusion is proposed. The thermo-mechanical model takes into account the normal stress dependence of yield stress, the free volume theory and the associated viscosity change within the shear band region. Temperature rise is obtained from the balance between the heat diffusion to the adjacent regions from a shear band and the heat generation due to the accumulated plastic work in a shear band. The parametric study has revealed that thermal effects play a minor role when the critical shear displacement is below 10 nm (as in nanoindentation) but become significant when the shear displacement accumulated in a shear band is of the order of hundreds of nanometers (as in uniaxial compression and in dynamic indentations). Finally, it is found that the normal stress plays a crucial role in the deformation behavior of BMGs by not only decreasing the time for shear band formation but also increasing the temperature rise significantly.

Single-walled carbon nanotube (SWCNT) and ethylene vinyl acetate (EVA) composites were synthesized in an internal mixer by melt mixing. The electrical conductivity as well as electromagnetic interference (EMI) shielding effectiveness (SE) over the X-band (8–12 GHz) and microwave (200–2,000 MHz) frequency ranges of these composites were investigated. It was observed that the electrical conductivity of composites increases with increasing SWCNT loading. A percolation threshold of about 3.5 wt.% was obtained and the electrical conductivity of EVA was increased by ten orders of magnitude, from 10−14 to 10−4 Ω−1 cm−1. The effect of sample thickness on SE was investigated. The correlation between SE and conductivity of the composites is discussed. The experimental data showed that the SE of the composites containing higher carbon nanotube loadings (above 10 wt.%) could be used as an EMI shielding material and lower SWCNT loadings could be used for the dissipation of electrostatic charge.

Diatom frustules have been suggested for numerous nanotechnological applications. Experimental studies using nanoindenter have shown that the hardness and the stiffness of the frustules vary with location of indentation. To gain further insight, a computational framework has been developed where the Berkovich nanoindentation experiments were simulated by a rigid-deformable contact process. Three different approaches that provide progressively increasing level of understanding of the deformation behavior of frustules were adopted. The differences in the mechanical responses of the frustule due to variation of indentation location, size of pores, and distribution of pores were analyzed. It has been found that the effective stiffness of the frustule is linearly related to the porosity level and does not depend on the frustule size or its pore architecture. It has been shown that a 3D porous shell computational model is more appropriate to simulate the experimentally obtained mechanical response of diatom frustules.