Matching Items (14)
Filtering by
- All Subjects: Fatigue
- All Subjects: Finite element method
- Creators: Liu, Yongming
Description
Aluminum alloys and their composites are attractive materials for applications requiring high strength-to-weight ratios and reasonable cost. Many of these applications, such as those in the aerospace industry, undergo fatigue loading. An understanding of the microstructural damage that occurs in these materials is critical in assessing their fatigue resistance. Two distinct experimental studies were performed to further the understanding of fatigue damage mechanisms in aluminum alloys and their composites, specifically fracture and plasticity. Fatigue resistance of metal matrix composites (MMCs) depends on many aspects of composite microstructure. Fatigue crack growth behavior is particularly dependent on the reinforcement characteristics and matrix microstructure. The goal of this work was to obtain a fundamental understanding of fatigue crack growth behavior in SiC particle-reinforced 2080 Al alloy composites. In situ X-ray synchrotron tomography was performed on two samples at low (R=0.1) and at high (R=0.6) R-ratios. The resulting reconstructed images were used to obtain three-dimensional (3D) rendering of the particles and fatigue crack. Behaviors of the particles and crack, as well as their interaction, were analyzed and quantified. Four-dimensional (4D) visual representations were constructed to aid in the overall understanding of damage evolution. During fatigue crack growth in ductile materials, a plastic zone is created in the region surrounding the crack tip. Knowledge of the plastic zone is important for the understanding of fatigue crack formation as well as subsequent growth behavior. The goal of this work was to quantify the 3D size and shape of the plastic zone in 7075 Al alloys. X-ray synchrotron tomography and Laue microdiffraction were used to non-destructively characterize the volume surrounding a fatigue crack tip. The precise 3D crack profile was segmented from the reconstructed tomography data. Depth-resolved Laue patterns were obtained using differential-aperture X-ray structural microscopy (DAXM), from which peak-broadening characteristics were quantified. Plasticity, as determined by the broadening of diffracted peaks, was mapped in 3D. Two-dimensional (2D) maps of plasticity were directly compared to the corresponding tomography slices. A 3D representation of the plastic zone surrounding the fatigue crack was generated by superimposing the mapped plasticity on the 3D crack profile.
ContributorsHruby, Peter (Author) / Chawla, Nikhilesh (Thesis advisor) / Solanki, Kiran (Committee member) / Liu, Yongming (Committee member) / Arizona State University (Publisher)
Created2014
Description
The objective of this research is to develop robust, accurate, and adaptive algorithms in the framework of the extended finite element method (XFEM) for fracture analysis of highly heterogeneous materials with complex internal geometries. A key contribution of this work is the creation of novel methods designed to automate the incorporation of high-resolution data, e.g. from X-ray tomography, that can be used to better interpret the enormous volume of data generated in modern in-situ experimental testing. Thus new algorithms were developed for automating analysis of complex microstructures characterized by segmented tomographic images.
A centrality-based geometry segmentation algorithm was developed to accurately identify discrete inclusions and particles in composite materials where limitations in imaging resolution leads to spurious connections between particles in close contact.To allow for this algorithm to successfully segment geometry independently of particle size and shape, a relative centrality metric was defined to allow for a threshold centrality criterion for removal of voxels that spuriously connect distinct geometries.
To automate incorporation of microstructural information from high-resolution images, two methods were developed that initialize signed distance fields on adaptively-refined finite element meshes. The first method utilizes a level set evolution equation that is directly solved on the finite element mesh through Galerkins method. The evolution equation is formulated to produce a signed distance field that matches geometry defined by a set of voxels segmented from tomographic images. The method achieves optimal convergence for the order of elements used. In a second approach, the fast marching method is employed to initialize a distance field on a uniform grid which is then projected by least squares onto a finite element mesh. This latter approach is shown to be superior in speed and accuracy.
Lastly, extended finite element method simulations are performed for the analysis of particle fracture in metal matrix composites with realistic particle geometries initialized from X-ray tomographic data. In the simulations, particles fracture probabilistically through a Weibull strength distribution. The model is verified through comparisons with the experimentally-measured stress-strain response of the material as well as analysis of the fracture. Further, simulations are then performed to analyze the effect of mesh sensitivity, the effect of fracture of particles on their neighbors, and the role of a particles shape on its fracture probability.
A centrality-based geometry segmentation algorithm was developed to accurately identify discrete inclusions and particles in composite materials where limitations in imaging resolution leads to spurious connections between particles in close contact.To allow for this algorithm to successfully segment geometry independently of particle size and shape, a relative centrality metric was defined to allow for a threshold centrality criterion for removal of voxels that spuriously connect distinct geometries.
To automate incorporation of microstructural information from high-resolution images, two methods were developed that initialize signed distance fields on adaptively-refined finite element meshes. The first method utilizes a level set evolution equation that is directly solved on the finite element mesh through Galerkins method. The evolution equation is formulated to produce a signed distance field that matches geometry defined by a set of voxels segmented from tomographic images. The method achieves optimal convergence for the order of elements used. In a second approach, the fast marching method is employed to initialize a distance field on a uniform grid which is then projected by least squares onto a finite element mesh. This latter approach is shown to be superior in speed and accuracy.
Lastly, extended finite element method simulations are performed for the analysis of particle fracture in metal matrix composites with realistic particle geometries initialized from X-ray tomographic data. In the simulations, particles fracture probabilistically through a Weibull strength distribution. The model is verified through comparisons with the experimentally-measured stress-strain response of the material as well as analysis of the fracture. Further, simulations are then performed to analyze the effect of mesh sensitivity, the effect of fracture of particles on their neighbors, and the role of a particles shape on its fracture probability.
ContributorsYuan, Rui (Author) / Oswald, Jay (Thesis advisor) / Chawla, Nikhilesh (Committee member) / Liu, Yongming (Committee member) / Solanki, Kiran (Committee member) / Chen, Kangping (Committee member) / Arizona State University (Publisher)
Created2015
Description
Interstitial impurity atoms can significantly alter the chemical and physical properties of the host material. Oxygen impurity in HCP titanium is known to have a considerable strengthening effect mainly through interactions with dislocations. To better understand such an effect, first the role of oxygen on various slip planes in titanium is examined using generalized stacking fault energies (GSFE) computed by the first principles calculations. It is shown that oxygen can significantly increase the energy barrier to dislocation motion on most of the studied slip planes. Then the Peierls-Nabbaro model is utilized in conjunction with the GSFE to estimate the Peierls stress ratios for different slip systems. Using such information along with a set of tension and compression experiments, the parameters of a continuum scale crystal plasticity model, namely CRSS values, are calibrated. Effect of oxygen content on the macroscopic stress-strain response is further investigated through experiments on oxygen-boosted samples at room temperature. It is demonstrated that the crystal plasticity model can very well capture the effect of oxygen content on the global response of the samples. It is also revealed that oxygen promotes the slip activity on the pyramidal planes.
The effect of oxygen impurity on titanium is further investigated under high cycle fatigue loading. For that purpose, a two-step hierarchical crystal plasticity for fatigue predictions is presented. Fatigue indicator parameter is used as the main driving force in an energy-based crack nucleation model. To calculate the FIPs, high-resolution full-field crystal plasticity simulations are carried out using a spectral solver. A nucleation model is proposed and calibrated by the fatigue experimental data for notched titanium samples with different oxygen contents and under two load ratios. Overall, it is shown that the presented approach is capable of predicting the high cycle fatigue nucleation time. Moreover, qualitative predictions of microstructurally small crack growth rates are provided. The multi-scale methodology presented here can be extended to other material systems to facilitate a better understanding of the fundamental deformation mechanisms, and to effectively implement such knowledge in mesoscale-macroscale investigations.
The effect of oxygen impurity on titanium is further investigated under high cycle fatigue loading. For that purpose, a two-step hierarchical crystal plasticity for fatigue predictions is presented. Fatigue indicator parameter is used as the main driving force in an energy-based crack nucleation model. To calculate the FIPs, high-resolution full-field crystal plasticity simulations are carried out using a spectral solver. A nucleation model is proposed and calibrated by the fatigue experimental data for notched titanium samples with different oxygen contents and under two load ratios. Overall, it is shown that the presented approach is capable of predicting the high cycle fatigue nucleation time. Moreover, qualitative predictions of microstructurally small crack growth rates are provided. The multi-scale methodology presented here can be extended to other material systems to facilitate a better understanding of the fundamental deformation mechanisms, and to effectively implement such knowledge in mesoscale-macroscale investigations.
ContributorsGholami Bazehhour, Benyamin (Author) / Solanki, Kiran N (Thesis advisor) / Liu, Yongming (Committee member) / Oswald, Jay J (Committee member) / Rajagopalan, Jagannathan (Committee member) / Jiao, Yang (Committee member) / Arizona State University (Publisher)
Created2018
Description
Fatigue is a degradation process of materials that would lead to failure when materials are subjected to cyclic loadings. During past centuries, various of approaches have been proposed and utilized to help researchers understand the underlying theories of fatigue behavior of materials, as well as design engineering structures so that catastrophic disasters that arise from fatigue failure could be avoided. The stress-life approach is the most classical way that academia applies to analyze fatigue data, which correlates the fatigue lifetime with stress amplitudes during cyclic loadings. Fracture mechanics approach is another well-established way, by which people regard the cyclic stress intensity factor as the driving force during fatigue crack nucleation and propagation, and numerous models (such as the well-known Paris’ law) are developed by researchers.
The significant drawback of currently widely-used fatigue analysis approaches, nevertheless, is that they are all cycle-based, limiting researchers from digging into sub-cycle regime and acquiring real-time fatigue behavior data. The missing of such data further impedes academia from validating hypotheses that are related to real-time observations of fatigue crack nucleation and growth, thus the existence of various phenomena, such as crack closure, remains controversial.
In this thesis, both classical stress-life approach and fracture-mechanics-based approach are utilized to study the fatigue behavior of alloys. Distinctive material characterization instruments are harnessed to help collect and interpret key data during fatigue crack growth. Specifically, an investigation on the sub-cycle fatigue crack growth behavior is enabled by in-situ SEM mechanical testing, and a non-uniform growth mechanism within one loading cycle is confirmed by direct observation as well as image interpretation. Predictions based on proposed experimental procedure and observations show good match with cycle-based data from references, which indicates the credibility of proposed methodology and model, as well as their capability of being applied to a wide range of materials.
The significant drawback of currently widely-used fatigue analysis approaches, nevertheless, is that they are all cycle-based, limiting researchers from digging into sub-cycle regime and acquiring real-time fatigue behavior data. The missing of such data further impedes academia from validating hypotheses that are related to real-time observations of fatigue crack nucleation and growth, thus the existence of various phenomena, such as crack closure, remains controversial.
In this thesis, both classical stress-life approach and fracture-mechanics-based approach are utilized to study the fatigue behavior of alloys. Distinctive material characterization instruments are harnessed to help collect and interpret key data during fatigue crack growth. Specifically, an investigation on the sub-cycle fatigue crack growth behavior is enabled by in-situ SEM mechanical testing, and a non-uniform growth mechanism within one loading cycle is confirmed by direct observation as well as image interpretation. Predictions based on proposed experimental procedure and observations show good match with cycle-based data from references, which indicates the credibility of proposed methodology and model, as well as their capability of being applied to a wide range of materials.
ContributorsLiu, Siying (Author) / Liu, Yongming (Thesis advisor) / Jiao, Yang (Committee member) / Nian, Qiong (Committee member) / Arizona State University (Publisher)
Created2018
Description
Widespread knowledge of fracture mechanics is mostly based on previous models that generalize crack growth in materials over several loading cycles. The objective of this project is to characterize crack growth that occurs in titanium alloys, specifically Grade 5 Ti-6Al-4V, at the sub-cycle scale, or within a single loading cycle. Using scanning electron microscopy (SEM), imaging analysis is performed to observe crack behavior at ten loading steps throughout the loading and unloading paths. Analysis involves measuring the incremental crack growth and crack tip opening displacement (CTOD) of specimens at loading ratios of 0.1, 0.3, and 0.5. This report defines the relationship between crack growth and the stress intensity factor, K, of the specimens, as well as the relationship between the R-ratio and stress opening level. The crack closure phenomena and effect of microcracks are discussed as they influence the crack growth behavior. This method has previously been used to characterize crack growth in Al 7075-T6. The results for Ti-6Al-4V are compared to these previous findings in order to strengthen conclusions about crack growth behavior.
ContributorsNazareno, Alyssa Noelle (Author) / Liu, Yongming (Thesis director) / Jiao, Yang (Committee member) / Mechanical and Aerospace Engineering Program (Contributor) / Barrett, The Honors College (Contributor)
Created2018-05
Description
In-situ fatigue damage diagnosis and prognosis is a challenging problem for both metallic and composite materials and structures. There are various uncertainties arising from material properties, component geometries, measurement noise, feature extraction techniques, and modeling errors. It is essential to manage and incorporate these uncertainties in order to achieve accurate damage detection and remaining useful life (RUL) prediction.
The aim of this study is to develop an integrated fatigue damage diagnosis and prognosis framework for both metallic and composite materials. First, Lamb waves are used as the in-situ damage detection technique to interrogate the damaged structures. Both experimental and numerical analysis for the Lamb wave propagation within aluminum are conducted. The RUL of lap joints under variable and constant fatigue loading is predicted using the Bayesian updating by incorporating damage detection information and various sources of uncertainties. Following this, the effect of matrix cracking and delamination in composite laminates on the Lamb wave propagation is investigated and a generalized probabilistic delamination size and location detection framework using Bayesian imaging method (BIM) is proposed and validated using the composite fatigue testing data. The RUL of the open-hole specimen is predicted using the overall stiffness degradation under fatigue loading. Next, the adjoint method-based damage detection framework is proposed considering the physics of heat conduction or elastic wave propagation. Different from the classical wave propagation-based method, the received signal under pristine condition is not necessary for estimating the damage information. This method can be successfully used for arbitrary damage location and shape profiling for any materials with higher accuracy and resolution. Finally, some conclusions and future work are generated based on the current investigation.
The aim of this study is to develop an integrated fatigue damage diagnosis and prognosis framework for both metallic and composite materials. First, Lamb waves are used as the in-situ damage detection technique to interrogate the damaged structures. Both experimental and numerical analysis for the Lamb wave propagation within aluminum are conducted. The RUL of lap joints under variable and constant fatigue loading is predicted using the Bayesian updating by incorporating damage detection information and various sources of uncertainties. Following this, the effect of matrix cracking and delamination in composite laminates on the Lamb wave propagation is investigated and a generalized probabilistic delamination size and location detection framework using Bayesian imaging method (BIM) is proposed and validated using the composite fatigue testing data. The RUL of the open-hole specimen is predicted using the overall stiffness degradation under fatigue loading. Next, the adjoint method-based damage detection framework is proposed considering the physics of heat conduction or elastic wave propagation. Different from the classical wave propagation-based method, the received signal under pristine condition is not necessary for estimating the damage information. This method can be successfully used for arbitrary damage location and shape profiling for any materials with higher accuracy and resolution. Finally, some conclusions and future work are generated based on the current investigation.
ContributorsPeng, Tishun (Author) / Liu, Yongming (Thesis advisor) / Chattopadhyay, Aditi (Committee member) / Mignolet, Marc (Committee member) / Papandreou-Suppappola, Antonia (Committee member) / Tang, Pingbo (Committee member) / Arizona State University (Publisher)
Created2016
Description
An integrated experimental and numerical investigation for laser-generated optoacoustic wave propagation in structural materials is performed. First, a multi-physics simulation model is proposed to simulate the pulsed laser as a point heat source which hits the surface of an aluminum sheet. The pulsed laser source can generate a localized heating on the surface of the plate and induce an in-plane stress wave. ANSYS – a finite element analysis software – is used to build the 3D model and a coupled thermal-mechanical simulation is performed in which the heat flux is determined by an empirical laser-heat conversion relationship. The displacement and stress field-histories are obtained to get the time of arrival and wave propagation speed of the stress wave. The effect of an added point mass is investigated in detail to observe the local material perturbation and remote wave signals. Following this, the experimental investigation of optoacoustic wave is also performed. A new experimental setup and control is developed and assembled in-house. Various laser firing parameters are investigated experimentally and the optimal combination is used for the experimental testing. Matrix design for different testing conditions is also proposed to include the effect of wave path, sampling procedure, and local point mass on the optoacoustic wave propagation. The developed numerical simulation results are validated with experimental observations. It is shown that the proposed experimental setup can offer a potential fast scanning method for damage detection (local property change) for plate-like structural component.
ContributorsLiu, Chen (Author) / Liu, Yongming (Thesis advisor) / Wang, Liping (Committee member) / Jiao, Yang (Committee member) / Arizona State University (Publisher)
Created2016
Description
In this paper, at first, analytical formulation of J-integral for a non-local particle model (VCPM) using atomic scale finite element method is proposed for fracture analysis of 2D solids. A brief review of classical continuum-based J-integral and anon-local lattice particle method is given first. Following this, detailed derivation for the J-integral in discrete particle system is given using the energy equivalence and stress-tensor mapping between the continuum mechanics and lattice-particle system.With the help of atomistic finite element method, the J-integral is expressed as a summation of the corresponding terms in the particle system.
Secondly, a coupling algorithm between a non-local particle method (VCPM) and the classical finite element method (FEM) is discussed to gain the advantages of both methods for fracture analysis in large structures. In this algorithm, the discrete VCPM particle and the continuum FEM domains are solved within a unified theoretical framework. A transitional element technology is developed to smoothly link the 10-particles element with the traditional FEM elements to guaranty the continuity and consistency at the coupling interface. An explicit algorithm for static simulation is developed.
Finally, numerical examples are illustrated for the accuracy, convergence, and path-independence of the derived J-integral formulation. Discussions on the comparison with alternative estimation methods and potential application for fracture simulation are given. The accuracy and efficiency of the coupling algorithm are tested by several benchmark problems such as static crack simulation.
Secondly, a coupling algorithm between a non-local particle method (VCPM) and the classical finite element method (FEM) is discussed to gain the advantages of both methods for fracture analysis in large structures. In this algorithm, the discrete VCPM particle and the continuum FEM domains are solved within a unified theoretical framework. A transitional element technology is developed to smoothly link the 10-particles element with the traditional FEM elements to guaranty the continuity and consistency at the coupling interface. An explicit algorithm for static simulation is developed.
Finally, numerical examples are illustrated for the accuracy, convergence, and path-independence of the derived J-integral formulation. Discussions on the comparison with alternative estimation methods and potential application for fracture simulation are given. The accuracy and efficiency of the coupling algorithm are tested by several benchmark problems such as static crack simulation.
ContributorsZope, Jayesh (Author) / Liu, Yongming (Thesis advisor) / Oswald, Jay (Committee member) / Jiang, Hanqing (Committee member) / Arizona State University (Publisher)
Created2016
Description
Composite materials are now beginning to provide uses hitherto reserved for metals in structural systems such as airframes and engine containment systems, wraps for repair and rehabilitation, and ballistic/blast mitigation systems. These structural systems are often subjected to impact loads and there is a pressing need for accurate prediction of deformation, damage and failure. There are numerous material models that have been developed to analyze the dynamic impact response of polymer matrix composites. However, there are key features that are missing in those models that prevent them from providing accurate predictive capabilities. In this dissertation, a general purpose orthotropic elasto-plastic computational constitutive material model has been developed to predict the response of composites subjected to high velocity impacts. The constitutive model is divided into three components – deformation model, damage model and failure model, with failure to be added at a later date. The deformation model generalizes the Tsai-Wu failure criteria and extends it using a strain-hardening-based orthotropic yield function with a non-associative flow rule. A strain equivalent formulation is utilized in the damage model that permits plastic and damage calculations to be uncoupled and capture the nonlinear unloading and local softening of the stress-strain response. A diagonal damage tensor is defined to account for the directionally dependent variation of damage. However, in composites it has been found that loading in one direction can lead to damage in multiple coordinate directions. To account for this phenomena, the terms in the damage matrix are semi-coupled such that the damage in a particular coordinate direction is a function of the stresses and plastic strains in all of the coordinate directions. The overall framework is driven by experimental tabulated temperature and rate-dependent stress-strain data as well as data that characterizes the damage matrix and failure. The developed theory has been implemented in a commercial explicit finite element analysis code, LS-DYNA®, as MAT213. Several verification and validation tests using a commonly available carbon-fiber composite, Toyobo’s T800/F3900, have been carried and the results show that the theory and implementation are efficient, robust and accurate.
ContributorsHoffarth, Canio (Author) / Rajan, Subramaniam D. (Thesis advisor) / Goldberg, Robert (Committee member) / Neithalath, Narayanan (Committee member) / Mobasher, Barzin (Committee member) / Liu, Yongming (Committee member) / Arizona State University (Publisher)
Created2016
Description
Cohesive zone model is one of the most widely used model for fracture analysis, but still remains open ended field for research. The earlier works using the cohesive zone model and Extended finite element analysis (XFEM) have been briefly introduced followed by an elaborate elucidation of the same concepts. Cohesive zone model in conjugation with XFEM is used for analysis in static condition in order to check its applicability in failure analysis. A real time setup of pipeline failure due to impingement is analyzed along with a detailed parametric study to understand the influence of the prominent design variable. After verifying its good applicability, a creep model is built for analysis where the cohesive zone model with XFEM is used for a time dependent creep loading. The challenge in this simulation was to achieve coupled behavior of cracks initiation and propagation along with creep loading. By using Design of Experiment, the results from numerical simulation were used to build an equation for life prediction for creep loading condition. The work was further extended to account for fatigue damage accumulation for high cycle fatigue loading in cohesive elements. A model was conceived to account for damage due to fatigue loading along within cohesive zone model for cohesive elements in ABAQUS simulation software. The model was verified by comparing numerical modelling of Double cantilever beam under high cycle fatigue loading and experiment results from literature. The model was also applied to a major industrial problem of blistering in Cured-In-Plane liner pipelines and a demonstration of its failure is shown. In conclusion, various models built on cohesive zone to address static and time dependent loading with real time scenarios and future scope of work in this field is discussed.
ContributorsChandrasekhar, Vishal (Author) / Liu, Yongming (Thesis advisor) / Oswald, Jay (Committee member) / Jiang, Hanqing (Committee member) / Arizona State University (Publisher)
Created2016