Fatigue crack propagation analysis of plaque rupture

Rupture of atheromatous plaque is the major cause of stroke or heart attack. Considering that the cardiovascular system is a classic fatigue environment, plaque rupture was treated as a chronic fatigue crack growth process in this study. Fracture mechanics theory was introduced to describe the stress status at the crack tip and Paris' law was used to calculate the crack growth rate. The effect of anatomical variation of an idealized plaque cross-section model was investigated. The crack initiation was considered to be either at the maximum circumferential stress location or at any other possible locations around the lumen. Although the crack automatically initialized at the maximum circumferential stress location usually propagated faster than others, it was not necessarily the most critical location where the fatigue life reached its minimum. We found that the fatigue life was minimum for cracks initialized in the following three regions: the midcap zone, the shoulder zone, and the backside zone. The anatomical variation has a significant influence on the fatigue life. Either a decrease in cap thickness or an increase in lipid pool size resulted in a significant decrease in fatigue life. Comparing to the previously used stress analysis, this fatigue model provides some possible explanations of plaque rupture at a low stress level in a pulsatile cardiovascular environment, and the method proposed here may be useful for further investigation of the mechanism of plaque rupture based on in vivo patient data.

MRI-based fatigue analysis predicts plaque vulnerability: A study comparing symptomatic and asymptomatic individuals

BACKGROUND: Rupture of atheromatous plaque in the carotid artery often leads to thrombosis and subsequent stroke. The mechanism of plaque rupture is not entirely clear but is thought to be a multi-factorial process involving thinning and weakening of the fibrous cap and biomechanical stress as the trigger leading to plaque rupture. As the cardiovascular system is a classic fatigue environment, the weakening of plaque leading to rupture may be a fatigue process, which is a symptomatically quiescent but potentially progressive failure process. In this study, we used a fatigue analysis based on in vivo magnetic resonance imaging (MRI) to investigate the rupture initiation location, crack propagation path and fatigue life within plaques of asymptomatic and symptomatic individuals. METHODS: Forty non-consecutive subjects (20 symptomatic and 20 asymptomatic) underwent high-resolution multi-sequence in vivo MRI of the carotid bifurcation. Fatigue analysis was performed based on the plaque geometry derived from in vivo MRI of the carotid artery at the point of maximum stenosis. Paris’ Law in fracture mechanics is adopted to determine the fatigue crack growth rate. Incremental crack propagation was dynamically simulated based on stress distributions. Plaque initiation location, crack propagation path and fatigue cycle of symptomatic and asymptomatic individuals were compared. RESULTS: Cracks were often found to begin at the lumen wall at areas of stress concentration. The preferred rupture direction was radial from the lumen center. The crack initially advanced slowly but accelerated as it developed, depending on plaque morphology. The fatigue cycles of symptomatic plaques were significantly less than those in the asymptomatic group (2.3 ± 0.9 vs 3.1 ± 0.7 (x106); p = 0.003). CONCLUSIONS: The number of cycles to rupture in symptomatic patients was higher than those predicted in asymptomatic patients by fatigue analysis, suggesting the possibility that plaques with a less fatigue life may be more prone to be symptomatic and rupture. If further validated by large-scale longitudinal studies, fatigue analysis based on high resolution in vivo MRI could potentially act as a useful tool for risk assessment of carotid atheroma.

Scale-dependent fracture mode in Cu-Ni laminate composites

Fracture behavior of Cu-Ni laminate composites has been investigated by tensile testing. It was found that as the individual layer thickness decreases from 100 to 20nm, the resultant fracture angle of the Cu-Ni laminate changes from 72 degrees to 50 degrees. Cross-sectional observations reveal that the fracture of the Ni layers transforms from opening to shear mode as the layer thickness decreases while that of the Cu layers keeps shear mode. Competition mechanisms were proposed to understand the variation in fracture mode of the metallic laminate composites associated with length scale.

Development of a multi-scale finite element model of the osteoporotic lumbar vertebral body for the investigation of apparent level vertebra mechanics and micro-level trabecular mechanics.

Osteoporotic spinal fractures are a major concern in ageing Western societies. This study develops a multi-scale finite element (FE) model of the osteoporotic lumbar vertebral body to study the mechanics of vertebral compression fracture at both the apparent (whole vertebral body) and micro-structural (internal trabecular bone core)levels. Model predictions were verified against experimental data, and found to provide a reasonably good representation of the mechanics of the osteoporotic vertebral body. This novel modelling methodology will allow detailed investigation of how trabecular bone loss in osteoporosis affects vertebral stiffness and strength in the lumbar spine.

Anisotropic damage mechanics as a novel approach to improve pre-and post-failure borehole stability analysis

Anisotropic damage distribution and evolution have a profound effect on borehole stress concentrations. Damage evolution is an irreversible process that is not adequately described within classical equilibrium thermodynamics. Therefore, we propose a constitutive model, based on non-equilibrium thermodynamics, that accounts for anisotropic damage distribution, anisotropic damage threshold and anisotropic damage evolution. We implemented this constitutive model numerically, using the finite element method, to calculate stress–strain curves and borehole stresses. The resulting stress–strain curves are distinctively different from linear elastic-brittle and linear elastic-ideal plastic constitutive models and realistically model experimental responses of brittle rocks. We show that the onset of damage evolution leads to an inhomogeneous redistribution of material properties and stresses along the borehole wall. The classical linear elastic-brittle approach to borehole stability analysis systematically overestimates the stress concentrations on the borehole wall, because dissipative strain-softening is underestimated. The proposed damage mechanics approach explicitly models dissipative behaviour and leads to non-conservative mud window estimations. Furthermore, anisotropic rocks with preferential planes of failure, like shales, can be addressed with our model.

Effect of crack depth and specimen width on fracture toughness of a carbon steel in the ductile-brittle transition region

The effects of crack depth (a/W) and specimen width W on the fracture toughness and ductile±brittle transition have been investigated using three-point bend specimens. Finite element analysis is employed to obtain the stress-strain fields ahead of the crack tip. The results show that both normalized crack depth (a/W) and specimen width (W) affect the fracture toughness and ductile±brittle fracture transition. The measured crack tip opening displacement decreases and ductile±brittle transition occurs with increasing crack depth (a/W) from 0.1 to 0.2 and 0.3. At a fixed a/W (0.2 or 0.3), all specimens fail by cleavage prior to ductile tearing when specimen width W increases from 25 to 40 and 50 mm. The lower bound fracture toughness is not sensitive to crack depth and specimen width. Finite element analysis shows that the opening stress in the remaining ligament is elevated with increasing crack depth or specimen width due to the increase of in-plane constraint. The average local cleavage stress is dependent on both crack depth and specimen width but its lower bound value is not sensitive to constraint level. No fixed distance can be found from the cleavage initiation site to the crack tip and this distance increases gradually with decreasing inplane constraint.

Mechanics of Polyropylene-Seed-Coat-Fibres Composites AndPolyropylene –Wood Fibres Composites-A Comparative Study

Forestry by-products have potential applications as components of wood composites. Replacement of conventional pine radiata wood-fibres by the fibres from the seeds (SCF) of the by-products, require determining and optimizing the mechanical properties to producing highest quality products. Response to mechanical stress is an important aspect to consider towards partial or full replacement of the wood-fibres by SCFs. In the present study the critical strain energy release rate, and the fracture toughness are derived from the published data. The present work uses rules of mixture to derive the mechanical and the physical properties of the SCF and relates the performance of the composites of the wood-fibres and the SCF to chemical composition, dispersion, weight and Vf of the fibres. We have also derived the Gc, the critical strain energy release rate, KIC, the fracture toughness of the composites.

Development of a novel experimental model to investigate the influence of mechanics on bone healing

This dissertation proposed a novel experimental model combining a defect configuration with an active instrumented fixation device to investigate the influence of mechanics on bone healing. The proposed defect configuration aimed to minimise physiological loading within an experimental fracture gap and the instrumented fixator was used for the application of controlled displacements and in vivo stiffness monitoring of the healing process. This thesis has provided a novel approach to advance current knowledge and understanding of mechanobiology, which has been limited in previous experimental models.

The mechanics of disc degeneration from an engineering and computational modelling viewpoint

The anatomy and microstructure of the spine and in particular the intervertebral disc are intimately linked to how they operate in vivo and how they distribute loads to the adjacent musculature and bony anatomy. The degeneration of the intervertebral discs may be characterised by a loss of hydration, loss of disc height, a granular texture and the presence of annular lesions. As such, degeneration of the intervertebral discs compromises the mechanical integrity of their components and results in adaption and modification in the mechanical means by which loads are distributed between adjacent spinal motion segments.

Simulation of the nutrient supply in fracture healing

The healing process for bone fractures is sensitive to mechanical stability and blood supply at the fracture site. Most currently available mechanobiological algorithms of bone healing are based solely on mechanical stimuli, while the explicit analysis of revascularization and its influences on the healing process have not been thoroughly investigated in the literature. In this paper, revascularization was described by two separate processes: angiogenesis and nutrition supply. The mathematical models for angiogenesis and nutrition supply have been proposed and integrated into an existing fuzzy algorithm of fracture healing. The computational algorithm of fracture healing, consisting of stress analysis, analyses of angiogenesis and nutrient supply, and tissue differentiation, has been tested on and compared with animal experimental results published previously. The simulation results showed that, for a small and medium-sized fracture gap, the nutrient supply is sufficient for bone healing, for a large fracture gap, non-union may be induced either by deficient nutrient supply or inadequate mechanical conditions. The comparisons with experimental results demonstrated that the improved computational algorithm is able to simulate a broad spectrum of fracture healing cases and to predict and explain delayed unions and non-union induced by large gap sizes and different mechanical conditions. The new algorithm will allow the simulation of more realistic clinical fracture healing cases with various fracture gaps and geometries and may be helpful to optimise implants and methods for fracture fixation.

Modelling the effects of bone fragment contact in fracture healing

The fracture healing process is modulated by the mechanical environment created by imposed loads and motion between the bone fragments. Contact between the fragments obviously results in a significantly different stress and strain environment to a uniform fracture gap containing only soft tissue (e.g. haematoma). The assumption of the latter in existing computational models of the healing process will hence exaggerate the inter-fragmentary strain in many clinically-relevant cases. To address this issue, we introduce the concept of a contact zone that represents a variable degree of contact between cortices by the relative proportions of bone and soft tissue present. This is introduced as an initial condition in a two-dimensional iterative finite element model of a healing tibial fracture, in which material properties are defined by the volume fractions of each tissue present. The algorithm governing the formation of cartilage and bone in the fracture callus uses fuzzy logic rules based on strain energy density resulting from axial compression. The model predicts that increasing the degree of initial bone contact reduces the amount of callus formed (periosteal callus thickness 3.1mm without contact, down to 0.5mm with 10% bone in contact zone). This is consistent with the greater effective stiffness in the contact zone and hence, a smaller inter-fragmentary strain. These results demonstrate that the contact zone strategy reasonably simulates the differences in the healing sequence resulting from the closeness of reduction.