Matching Items (77)
Filtering by
- All Subjects: Mechanical Engineering
- All Subjects: Aerospace Engineering
- Creators: Liu, Yongming
- Member of: Theses and Dissertations
Description
There is a concerted effort in developing robust systems health monitoring/management (SHM) technology as a means to reduce the life cycle costs, improve availability, extend life and minimize downtime of various platforms including aerospace and civil infrastructure. The implementation of a robust SHM system requires a collaborative effort in a variety of areas such as sensor development, damage detection and localization, physics based models, and prognosis models for residual useful life (RUL) estimation. Damage localization and prediction is further complicated by geometric, material, loading, and environmental variabilities. Therefore, it is essential to develop robust SHM methodologies by taking into account such uncertainties. In this research, damage localization and RUL estimation of two different physical systems are addressed: (i) fatigue crack propagation in metallic materials under complex multiaxial loading and (ii) temporal scour prediction near bridge piers. With little modifications, the methodologies developed can be applied to other systems.
Current practice in fatigue life prediction is based on either physics based modeling or data-driven methods, and is limited to predicting RUL for simple geometries under uniaxial loading conditions. In this research, crack initiation and propagation behavior under uniaxial and complex biaxial fatigue loading is addressed. The crack propagation behavior is studied by performing extensive material characterization and fatigue testing under in-plane biaxial loading, both in-phase and out-of-phase, with different biaxiality ratios. A hybrid prognosis model, which combines machine learning with physics based modeling, is developed to account for the uncertainties in crack propagation and fatigue life prediction due to variabilities in material microstructural characteristics, crack localization information and environmental changes. The methodology iteratively combines localization information with hybrid prognosis models using sequential Bayesian techniques. The results show significant improvements in the localization and prediction accuracy under varying temperature.
For civil infrastructure, especially bridges, pier scour is a major failure mechanism. Currently available techniques are developed from a design perspective and provide highly conservative scour estimates. In this research, a fully probabilistic scour prediction methodology is developed using machine learning to accurately predict scour in real-time under varying flow conditions.
Current practice in fatigue life prediction is based on either physics based modeling or data-driven methods, and is limited to predicting RUL for simple geometries under uniaxial loading conditions. In this research, crack initiation and propagation behavior under uniaxial and complex biaxial fatigue loading is addressed. The crack propagation behavior is studied by performing extensive material characterization and fatigue testing under in-plane biaxial loading, both in-phase and out-of-phase, with different biaxiality ratios. A hybrid prognosis model, which combines machine learning with physics based modeling, is developed to account for the uncertainties in crack propagation and fatigue life prediction due to variabilities in material microstructural characteristics, crack localization information and environmental changes. The methodology iteratively combines localization information with hybrid prognosis models using sequential Bayesian techniques. The results show significant improvements in the localization and prediction accuracy under varying temperature.
For civil infrastructure, especially bridges, pier scour is a major failure mechanism. Currently available techniques are developed from a design perspective and provide highly conservative scour estimates. In this research, a fully probabilistic scour prediction methodology is developed using machine learning to accurately predict scour in real-time under varying flow conditions.
ContributorsNeerukatti, Rajesh Kumar (Author) / Chattopadhyay, Aditi (Thesis advisor) / Jiang, Hanqing (Committee member) / Liu, Yongming (Committee member) / Rajadas, John (Committee member) / Yekani Fard, Masoud (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
Finite element simulations modeling the hydrodynamic impact loads subjected to an elastomeric coating were performed to develop an understanding of the performance and failure mechanisms of protective coatings for cavitating environments.
In this work, two major accomplishments were achieved: 1) scaling laws were developed from hydrodynamic principles and numerical simulations to allow conversion of measured distributions of pressure peaks in a cavitating flow to distributions of microscopic impact loadings modeling individual bubble collapse events, and 2) a finite strain, thermo-mechanical material model for polyurea-based elastomers was developed using a logarithmic rate formulation and implemented into an explicit finite element code.
Combining the distribution of microscopic impact loads and finite element modeling, a semi-quantitative predictive framework is created to calculate the energy dissipation within the coating which can further the understanding of temperature induced coating failures.
The influence of coating thickness and elastomer rheology on the dissipation of impact energies experienced in cavitating flows has also been explored.
The logarithmic formulation has many desired features for the polyurea constitutive model, such as objectivity, integrability, and additive decomposition compatibility.
A review and discussion on the kinematics in large deformation, including a comparison between Lagrangian and Eulerian descriptions, are presented to explain the issues in building rate-dependent constitutive models in finite strains.
When comparing the logarithmic rate with other conventional rates in test examples, the logarithmic rate shows a better conservation of objectivity and integrability.
The modeling framework was validated by comparing predictions against temperatures measured within coatings subjected to a cavitating jet.
Both the experiments and models show that the temperatures generated, even under mild flow conditions, raise the coating temperature by a significant amount, suggesting that the failure of these coatings under more aggressive flows is thermally induced.
The models show that thin polyurea coatings synthesized with shorter molecular weight soft segments dissipate significantly less energy per impact and conduct heat more efficiently.
This work represents an important step toward understanding thermally induced failure in elastomers subjected to cavitating flows, which provides a foundation for design and optimization of coatings with enhanced erosion resistance.
In this work, two major accomplishments were achieved: 1) scaling laws were developed from hydrodynamic principles and numerical simulations to allow conversion of measured distributions of pressure peaks in a cavitating flow to distributions of microscopic impact loadings modeling individual bubble collapse events, and 2) a finite strain, thermo-mechanical material model for polyurea-based elastomers was developed using a logarithmic rate formulation and implemented into an explicit finite element code.
Combining the distribution of microscopic impact loads and finite element modeling, a semi-quantitative predictive framework is created to calculate the energy dissipation within the coating which can further the understanding of temperature induced coating failures.
The influence of coating thickness and elastomer rheology on the dissipation of impact energies experienced in cavitating flows has also been explored.
The logarithmic formulation has many desired features for the polyurea constitutive model, such as objectivity, integrability, and additive decomposition compatibility.
A review and discussion on the kinematics in large deformation, including a comparison between Lagrangian and Eulerian descriptions, are presented to explain the issues in building rate-dependent constitutive models in finite strains.
When comparing the logarithmic rate with other conventional rates in test examples, the logarithmic rate shows a better conservation of objectivity and integrability.
The modeling framework was validated by comparing predictions against temperatures measured within coatings subjected to a cavitating jet.
Both the experiments and models show that the temperatures generated, even under mild flow conditions, raise the coating temperature by a significant amount, suggesting that the failure of these coatings under more aggressive flows is thermally induced.
The models show that thin polyurea coatings synthesized with shorter molecular weight soft segments dissipate significantly less energy per impact and conduct heat more efficiently.
This work represents an important step toward understanding thermally induced failure in elastomers subjected to cavitating flows, which provides a foundation for design and optimization of coatings with enhanced erosion resistance.
ContributorsLiao, Xiao (Author) / Oswald, Jay (Thesis advisor) / Liu, Yongming (Committee member) / Jiang, Hanqing (Committee member) / Rajan, Subramaniam D. (Committee member) / Platte, Rodrigo (Committee member) / Arizona State University (Publisher)
Created2016
Description
A previously developed small time scale fatigue crack growth model is improved, modified and extended with an emphasis on creating the simplest models that maintain the desired level of accuracy for a variety of materials. The model provides a means of estimating load sequence effects by continuously updating the crack opening stress every cycle, in a simplified manner. One of the significant phenomena of the crack opening stress under negative stress ratio is the residual tensile stress induced by the applied compressive stress. A modified coefficient is introduced to determine the extent to which residual stress impact the crack closure and is observed to vary for different materials. Several other literature models for crack closure under constant loading are also reviewed and compared with the proposed model. The modified model is then shown to predict several sets of published test results under constant loading for a variety of materials.
The crack opening stress is formalized as a function of the plastic zone sizes at the crack tip and the current crack length, which provided a means of approximation, accounting for both acceleration and retardation effects in a simplified manner. A sensitivity parameter is introduced to modify the enlarged plastic zone due to overload, to better fit the delay cycles with the test data and is observed to vary for different materials. Furthermore, the interaction effect induced by the combination of overload and underload sequence is modeled by depleting the compressive plastic zone due to an overload with the tensile plastic zone due to an underload. A qualitative analysis showed the simulation capacity of the small time scale model under different load types. A good agreement between prediction and test data for several irregular load types proved the applicability of the small time scale model under variable amplitude loading.
The crack opening stress is formalized as a function of the plastic zone sizes at the crack tip and the current crack length, which provided a means of approximation, accounting for both acceleration and retardation effects in a simplified manner. A sensitivity parameter is introduced to modify the enlarged plastic zone due to overload, to better fit the delay cycles with the test data and is observed to vary for different materials. Furthermore, the interaction effect induced by the combination of overload and underload sequence is modeled by depleting the compressive plastic zone due to an overload with the tensile plastic zone due to an underload. A qualitative analysis showed the simulation capacity of the small time scale model under different load types. A good agreement between prediction and test data for several irregular load types proved the applicability of the small time scale model under variable amplitude loading.
ContributorsVenkatesan, Karthik Rajan (Author) / Liu, Yongming (Thesis advisor) / Oswald, Jay (Committee member) / Jiang, Hanqing (Committee member) / Arizona State University (Publisher)
Created2016
Stochastic multiscale modeling and statistical characterization of complex polymer matrix composites
Description
There are many applications for polymer matrix composite materials in a variety of different industries, but designing and modeling with these materials remains a challenge due to the intricate architecture and damage modes. Multiscale modeling techniques of composite structures subjected to complex loadings are needed in order to address the scale-dependent behavior and failure. The rate dependency and nonlinearity of polymer matrix composite materials further complicates the modeling. Additionally, variability in the material constituents plays an important role in the material behavior and damage. The systematic consideration of uncertainties is as important as having the appropriate structural model, especially during model validation where the total error between physical observation and model prediction must be characterized. It is necessary to quantify the effects of uncertainties at every length scale in order to fully understand their impact on the structural response. Material variability may include variations in fiber volume fraction, fiber dimensions, fiber waviness, pure resin pockets, and void distributions. Therefore, a stochastic modeling framework with scale dependent constitutive laws and an appropriate failure theory is required to simulate the behavior and failure of polymer matrix composite structures subjected to complex loadings. Additionally, the variations in environmental conditions for aerospace applications and the effect of these conditions on the polymer matrix composite material need to be considered. The research presented in this dissertation provides the framework for stochastic multiscale modeling of composites and the characterization data needed to determine the effect of different environmental conditions on the material properties. The developed models extend sectional micromechanics techniques by incorporating 3D progressive damage theories and multiscale failure criteria. The mechanical testing of composites under various environmental conditions demonstrates the degrading effect these conditions have on the elastic and failure properties of the material. The methodologies presented in this research represent substantial progress toward understanding the failure and effect of variability for complex polymer matrix composites.
ContributorsJohnston, Joel Philip (Author) / Chattopadhyay, Aditi (Thesis advisor) / Liu, Yongming (Committee member) / Jiang, Hanqing (Committee member) / Dai, Lenore (Committee member) / Rajadas, John (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
Description
This investigation develops small-size reduced order models (ROMs) that provide an accurate prediction of the response of only part of a structure, referred to as component-centric ROMs. Four strategies to construct such ROMs are presented, the first two of which are based on the Craig-Bampton Method and start with a set of modes for the component of interest (the β component). The response in the rest of the structure (the α component) induced by these modes is then determined and optimally represented by applying a Proper Orthogonal Decomposition strategy using Singular Value Decomposition. These first two methods are effectively basis reductions techniques of the CB basis. An approach based on the “Global - Local” Method generates the “global” modes by “averaging” the mass property over α and β comp., respectively (to extract a “coarse” model of α and β) and the “local” modes orthogonal to the “global” modes to add back necessary “information” for β. The last approach adopts as basis for the entire structure its linear modes which are dominant in the β component response. Then, the contributions of other modes in this part of the structure are approximated in terms of those of the dominant modes with close natural frequencies and similar mode shapes in the β component. In this manner, the non-dominant modal contributions are “lumped” onto the dominant ones, to reduce the number of modes for a prescribed accuracy. The four approaches are critically assessed on the structural finite element model of a 9-bay panel with the modal lumping-based method leading to the smallest sized ROMs. Therefore, it is extended to the nonlinear geometric situation and first recast as a rotation of the modal basis to achieve unobservable modes. In the linear case, these modes completely disappear from the formulation owing to orthogonality. In the nonlinear case, however, the generalized coordinates of these modes are still present in the nonlinear terms of the observable modes. A closure-type algorithm is then proposed to eliminate the unobserved generalized coordinates. This approach, its accuracy and computational savings, was demonstrated on a simple beam model and the 9-bay panel model.
ContributorsWang, Yuting (Author) / Mignolet, Marc P (Thesis advisor) / Jiang, Hanqing (Committee member) / Liu, Yongming (Committee member) / Oswald, Jay (Committee member) / Rajan, Subramaniam D. (Committee member) / Spottswood, Stephen M (Committee member) / Arizona State University (Publisher)
Created2017
Description
This study uses Computational Fluid Dynamics (CFD) modeling to analyze the
dependence of wind power potential and turbulence intensity on aerodynamic design of a
special type of building with a nuzzle-like gap at its rooftop. Numerical simulations using
ANSYS Fluent are carried out to quantify the above-mentioned dependency due to three
major geometric parameters of the building: (i) the height of the building, (ii) the depth of
the roof-top gap, and (iii) the width of the roof-top gap. The height of the building is varied
from 8 m to 24 m. Likewise, the gap depth is varied from 3 m to 5 m and the gap width
from 2 m to 4 m. The aim of this entire research is to relate these geometric parameters of
the building to the maximum value and the spatial pattern of wind power potential across
the roof-top gap. These outcomes help guide the design of the roof-top geometry for wind
power applications and determine the ideal position for mounting a micro wind turbine.
From these outcomes, it is suggested that the wind power potential is greatly affected by
the increasing gap width or gap depth. It, however, remains insensitive to the increasing
building height, unlike turbulence intensity which increases with increasing building
height. After performing a set of simulations with varying building geometry to quantify
the wind power potential before the installation of a turbine, another set of simulations is
conducted by installing a static turbine within the roof-top gap. The results from the latter
are used to further adjust the estimate of wind power potential. Recommendations are made
for future applications based on the findings from the numerical simulations.
dependence of wind power potential and turbulence intensity on aerodynamic design of a
special type of building with a nuzzle-like gap at its rooftop. Numerical simulations using
ANSYS Fluent are carried out to quantify the above-mentioned dependency due to three
major geometric parameters of the building: (i) the height of the building, (ii) the depth of
the roof-top gap, and (iii) the width of the roof-top gap. The height of the building is varied
from 8 m to 24 m. Likewise, the gap depth is varied from 3 m to 5 m and the gap width
from 2 m to 4 m. The aim of this entire research is to relate these geometric parameters of
the building to the maximum value and the spatial pattern of wind power potential across
the roof-top gap. These outcomes help guide the design of the roof-top geometry for wind
power applications and determine the ideal position for mounting a micro wind turbine.
From these outcomes, it is suggested that the wind power potential is greatly affected by
the increasing gap width or gap depth. It, however, remains insensitive to the increasing
building height, unlike turbulence intensity which increases with increasing building
height. After performing a set of simulations with varying building geometry to quantify
the wind power potential before the installation of a turbine, another set of simulations is
conducted by installing a static turbine within the roof-top gap. The results from the latter
are used to further adjust the estimate of wind power potential. Recommendations are made
for future applications based on the findings from the numerical simulations.
ContributorsKailkhura, Gargi (Author) / Huang, Huei-Ping (Thesis advisor) / Rajagopalan, Jagannathan (Committee member) / Forzani, Erica (Committee member) / Arizona State University (Publisher)
Created2017
Description
The study of response of various materials to intense dynamic loading events,
such as shock loading due to high-velocity impacts, is extremely important in a wide
variety of military and industrial applications. Shock loading triggers extreme states,
leading to high pressures and strain rates, and neglecting strength is a typical
approximation under such conditions. However, recent results have shown that strength
effects are larger than expected, so they must be taken into account. Recently,
hydrodynamic instabilities, the most common being the Rayleigh-Taylor (RTI) and
Richtmyer-Meshkov (RMI) instabilities, have been used to infer the dynamic strength of
materials at high pressure conditions. In our experiments and simulations, a novel RMI
approach is used, in which periodic surface perturbations are made on high purity
aluminium target, which was laser ablated to create a rippled shock front. Due to the
slow linear growth rate of RMI, the evolution of the perturbations on the back surface of
the sample as a result of the rippled shock can be measured via Transient Imaging
Displacement Interferometry (TIDI). The velocity history at the free surface was
recorded by spatially resolved laser velocimetry. These measurements were compared
with the results from the simulations, which were implemented using rate independent
and rate dependent material models, to characterize the dynamic strength of the
material. Simulations using the elastic-perfectly plastic model, which is rate
independent, failed to provide a value of dynamic yield strength that would match
experimental measurements of perturbation amplitudes. The Preston-Tonks-Wallace
(PTW) model, which is rate dependent model, worked well for aluminium. This model
was, in turn, used as a reference for calibrating the rate dependent Steinberg-Lund
model and the results from simulations using the calibration models were also compared
to experimental measurements.
such as shock loading due to high-velocity impacts, is extremely important in a wide
variety of military and industrial applications. Shock loading triggers extreme states,
leading to high pressures and strain rates, and neglecting strength is a typical
approximation under such conditions. However, recent results have shown that strength
effects are larger than expected, so they must be taken into account. Recently,
hydrodynamic instabilities, the most common being the Rayleigh-Taylor (RTI) and
Richtmyer-Meshkov (RMI) instabilities, have been used to infer the dynamic strength of
materials at high pressure conditions. In our experiments and simulations, a novel RMI
approach is used, in which periodic surface perturbations are made on high purity
aluminium target, which was laser ablated to create a rippled shock front. Due to the
slow linear growth rate of RMI, the evolution of the perturbations on the back surface of
the sample as a result of the rippled shock can be measured via Transient Imaging
Displacement Interferometry (TIDI). The velocity history at the free surface was
recorded by spatially resolved laser velocimetry. These measurements were compared
with the results from the simulations, which were implemented using rate independent
and rate dependent material models, to characterize the dynamic strength of the
material. Simulations using the elastic-perfectly plastic model, which is rate
independent, failed to provide a value of dynamic yield strength that would match
experimental measurements of perturbation amplitudes. The Preston-Tonks-Wallace
(PTW) model, which is rate dependent model, worked well for aluminium. This model
was, in turn, used as a reference for calibrating the rate dependent Steinberg-Lund
model and the results from simulations using the calibration models were also compared
to experimental measurements.
ContributorsGopalakrishnan, Ashish (Author) / Peralta, Pedro (Thesis advisor) / Rajagopalan, Jagannathan (Committee member) / Solanki, Kiran (Committee member) / Arizona State University (Publisher)
Created2017