Matching Items (98)
Filtering by
- All Subjects: Mechanical Engineering
- Creators: Liu, Yongming
- Creators: Oswald, Jay
- Member of: Theses and Dissertations
Description
Recently, the use of zinc oxide (ZnO) nanowires as an interphase in composite materials has been demonstrated to increase the interfacial shear strength between carbon fiber and an epoxy matrix. In this research work, the strong adhesion between ZnO and carbon fiber is investigated to elucidate the interactions at the interface that result in high interfacial strength. First, molecular dynamics (MD) simulations are performed to calculate the adhesive energy between bare carbon and ZnO. Since the carbon fiber surface has oxygen functional groups, these were modeled and MD simulations showed the preference of ketones to strongly interact with ZnO, however, this was not observed in the case of hydroxyls and carboxylic acid. It was also found that the ketone molecules ability to change orientation facilitated the interactions with the ZnO surface. Experimentally, the atomic force microscope (AFM) was used to measure the adhesive energy between ZnO and carbon through a liftoff test by employing highly oriented pyrolytic graphite (HOPG) substrate and a ZnO covered AFM tip. Oxygen functionalization of the HOPG surface shows the increase of adhesive energy. Additionally, the surface of ZnO was modified to hold a negative charge, which demonstrated an increase in the adhesive energy. This increase in adhesion resulted from increased induction forces given the relatively high polarizability of HOPG and the preservation of the charge on ZnO surface. It was found that the additional negative charge can be preserved on the ZnO surface because there is an energy barrier since carbon and ZnO form a Schottky contact. Other materials with the same ionic properties of ZnO but with higher polarizability also demonstrated good adhesion to carbon. This result substantiates that their induced interaction can be facilitated not only by the polarizability of carbon but by any of the materials at the interface. The versatility to modify the magnitude of the induced interaction between carbon and an ionic material provides a new route to create interfaces with controlled interfacial strength.
ContributorsGalan Vera, Magdian Ulises (Author) / Sodano, Henry A (Thesis advisor) / Jiang, Hanqing (Committee member) / Solanki, Kiran (Committee member) / Oswald, Jay (Committee member) / Speyer, Gil (Committee member) / Arizona State University (Publisher)
Created2013
Description
This research examines the current challenges of using Lamb wave interrogation methods to localize fatigue crack damage in a complex metallic structural component subjected to unknown temperatures. The goal of this work is to improve damage localization results for a structural component interrogated at an unknown temperature, by developing a probabilistic and reference-free framework for estimating Lamb wave velocities and the damage location. The methodology for damage localization at unknown temperatures includes the following key elements: i) a model that can describe the change in Lamb wave velocities with temperature; ii) the extension of an advanced time-frequency based signal processing technique for enhanced time-of-flight feature extraction from a dispersive signal; iii) the development of a Bayesian damage localization framework incorporating data association and sensor fusion. The technique requires no additional transducers to be installed on a structure, and allows for the estimation of both the temperature and the wave velocity in the component. Additionally, the framework of the algorithm allows it to function completely in an unsupervised manner by probabilistically accounting for all measurement origin uncertainty. The novel algorithm was experimentally validated using an aluminum lug joint with a growing fatigue crack. The lug joint was interrogated using piezoelectric transducers at multiple fatigue crack lengths, and at temperatures between 20°C and 80°C. The results showed that the algorithm could accurately predict the temperature and wave speed of the lug joint. The localization results for the fatigue damage were found to correlate well with the true locations at long crack lengths, but loss of accuracy was observed in localizing small cracks due to time-of-flight measurement errors. To validate the algorithm across a wider range of temperatures the electromechanically coupled LISA/SIM model was used to simulate the effects of temperatures. The numerical results showed that this approach would be capable of experimentally estimating the temperature and velocity in the lug joint for temperatures from -60°C to 150°C. The velocity estimation algorithm was found to significantly increase the accuracy of localization at temperatures above 120°C when error due to incorrect velocity selection begins to outweigh the error due to time-of-flight measurements.
ContributorsHensberry, Kevin (Author) / Chattopadhyay, Aditi (Thesis advisor) / Liu, Yongming (Committee member) / Papandreou-Suppappola, Antonia (Committee member) / Arizona State University (Publisher)
Created2013
Description
Shock loading is a complex phenomenon that can lead to failure mechanisms such as strain localization, void nucleation and growth, and eventually spall fracture. Studying incipient stages of spall damage is of paramount importance to accurately determine initiation sites in the material microstructure where damage will nucleate and grow and to formulate continuum models that account for the variability of the damage process due to microstructural heterogeneity. The length scale of damage with respect to that of the surrounding microstructure has proven to be a key aspect in determining sites of failure initiation. Correlations have been found between the damage sites and the surrounding microstructure to determine the preferred sites of spall damage, since it tends to localize at and around the regions of intrinsic defects such as grain boundaries and triple points. However, considerable amount of work still has to be done in this regard to determine the physics driving the damage at these intrinsic weak sites in the microstructure. The main focus of this research work is to understand the physical mechanisms behind the damage localization at these preferred sites. A crystal plasticity constitutive model is implemented with different damage criteria to study the effects of stress concentration and strain localization at the grain boundaries. A cohesive zone modeling technique is used to include the intrinsic strength of the grain boundaries in the simulations. The constitutive model is verified using single elements tests, calibrated using single crystal impact experiments and validated using bicrystal and multicrystal impact experiments. The results indicate that strain localization is the predominant driving force for damage initiation and evolution. The microstructural effects on theses damage sites are studied to attribute the extent of damage to microstructural features such as grain orientation, misorientation, Taylor factor and the grain boundary planes. The finite element simulations show good correlation with the experimental results and can be used as the preliminary step in developing accurate probabilistic models for damage nucleation.
ContributorsKrishnan, Kapil (Author) / Peralta, Pedro (Thesis advisor) / Mignolet, Marc (Committee member) / Sieradzki, Karl (Committee member) / Jiang, Hanqing (Committee member) / Oswald, Jay (Committee member) / Arizona State University (Publisher)
Created2013
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
Cavitation erosion is a significant cause of wear in marine components, such as impellers, propellers or rudders. While the erosion process has been widely studied on metals, the effect of cavitation on polymers is not well-understood. The stress response in metals differs greatly from that of polymers, e.g. rate and temperature effects are far more important, thus damage and wear mechanisms of polymers under cavitating flows are significantly different. In this work, heat-driven failure caused by viscous dissipation and void nucleation resulting from tensile stresses arising from stress wave reflections are investigated as two possible material failure mechanisms.
As a first step in developing a fundamental understanding of the cavitation erosion process on polymer surfaces, simulations are performed of the collapse of individual bubbles against a compliant surface e.g. metallic substrates with polyurea coatings. The surface response of collapse-driven impact loads is represented by a idealized, time-dependent, Gaussian pressure distribution on the surface. A two-dimensional distribution of load radii and durations is considered corresponding to characteristic of cavitating flows accelerated erosion experiments. Finite element simulations are performed to fit a response curve that relates the loading parameters to the energy dissipated in the coating and integrated with collapse statistics to generate an expected heat input into the coating.
The impulsive pressure, which is generated due to bubble collapse, impacts the material and generates intense shock waves. The stress waves within the material reflects by interaction with the substrate. A transient region of high tensile stress is produced by the interaction of these waves. Simulations suggests that maximum hydrostatic tension which cause failure of polyurea layer is observed in thick coating. Also, the dissipated viscous energy and corresponding temperature rise in a polyurea is calculated, and it is concluded that temperature has influence on deformation.
As a first step in developing a fundamental understanding of the cavitation erosion process on polymer surfaces, simulations are performed of the collapse of individual bubbles against a compliant surface e.g. metallic substrates with polyurea coatings. The surface response of collapse-driven impact loads is represented by a idealized, time-dependent, Gaussian pressure distribution on the surface. A two-dimensional distribution of load radii and durations is considered corresponding to characteristic of cavitating flows accelerated erosion experiments. Finite element simulations are performed to fit a response curve that relates the loading parameters to the energy dissipated in the coating and integrated with collapse statistics to generate an expected heat input into the coating.
The impulsive pressure, which is generated due to bubble collapse, impacts the material and generates intense shock waves. The stress waves within the material reflects by interaction with the substrate. A transient region of high tensile stress is produced by the interaction of these waves. Simulations suggests that maximum hydrostatic tension which cause failure of polyurea layer is observed in thick coating. Also, the dissipated viscous energy and corresponding temperature rise in a polyurea is calculated, and it is concluded that temperature has influence on deformation.
ContributorsPanwar, Ajay (Author) / Oswald, Jay (Thesis advisor) / Dooley, Kevin (Committee member) / Chen, Kangping (Committee member) / Arizona State University (Publisher)
Created2015
Description
Recent studies of the occurrence of post-flutter limit cycle oscillations (LCO) of the F-16 have provided good support to the long-standing hypothesis that this phenomenon involves a nonlinear structural damping. A potential mechanism for the appearance of nonlinearity in the damping are the nonlinear geometric effects that arise when the deformations become large enough to exceed the linear regime. In this light, the focus of this investigation is first on extending nonlinear reduced order modeling (ROM) methods to include viscoelasticity which is introduced here through a linear Kelvin-Voigt model in the undeformed configuration. Proceeding with a Galerkin approach, the ROM governing equations of motion are obtained and are found to be of a generalized van der Pol-Duffing form with parameters depending on the structure and the chosen basis functions. An identification approach of the nonlinear damping parameters is next proposed which is applicable to structures modeled within commercial finite element software.
The effects of this nonlinear damping mechanism on the post-flutter response is next analyzed on the Goland wing through time-marching of the aeroelastic equations comprising a rational fraction approximation of the linear aerodynamic forces. It is indeed found that the nonlinearity in the damping can stabilize the unstable aerodynamics and lead to finite amplitude limit cycle oscillations even when the stiffness related nonlinear geometric effects are neglected. The incorporation of these latter effects in the model is found to further decrease the amplitude of LCO even though the dominant bending motions do not seem to stiffen as the level of displacements is increased in static analyses.
The effects of this nonlinear damping mechanism on the post-flutter response is next analyzed on the Goland wing through time-marching of the aeroelastic equations comprising a rational fraction approximation of the linear aerodynamic forces. It is indeed found that the nonlinearity in the damping can stabilize the unstable aerodynamics and lead to finite amplitude limit cycle oscillations even when the stiffness related nonlinear geometric effects are neglected. The incorporation of these latter effects in the model is found to further decrease the amplitude of LCO even though the dominant bending motions do not seem to stiffen as the level of displacements is increased in static analyses.
ContributorsSong, Pengchao (Author) / Mignolet, Marc P (Thesis advisor) / Chattopadhyay, Aditi (Committee member) / Oswald, Jay (Committee member) / Arizona State University (Publisher)
Created2015
Description
The football helmet is a device used to help mitigate the occurrence of impact-related traumatic (TBI) and minor traumatic brain injuries (mTBI) in the game of American football. The current design methodology of using a hard shell with an energy absorbing liner may be adequate for minimizing TBI, however it has had less effect in minimizing mTBI. The latest research in brain injury mechanisms has established that the current design methodology has produced a helmet to reduce linear acceleration of the head. However, angular accelerations also have an adverse effect on the brain response, and must be investigated as a contributor of brain injury.
To help better understand how the football helmet design features effect the brain response during impact, this research develops a validated football helmet model and couples it with a full LS-DYNA human body model developed by the Global Human Body Modeling Consortium (v4.1.1). The human body model is a conglomeration of several validated models of different sections of the body. Of particular interest for this research is the Wayne State University Head Injury Model for modeling the brain. These human body models were validated using a combination of cadaveric and animal studies. In this study, the football helmet was validated by laboratory testing using drop tests on the crown of the helmet. By coupling the two models into one finite element model, the brain response to impact loads caused by helmet design features can be investigated. In the present research, LS-DYNA is used to study a helmet crown impact with a rigid steel plate so as to obtain the strain-rate, strain, and stress experienced in the corpus callosum, midbrain, and brain stem as these anatomical regions are areas of concern with respect to mTBI.
To help better understand how the football helmet design features effect the brain response during impact, this research develops a validated football helmet model and couples it with a full LS-DYNA human body model developed by the Global Human Body Modeling Consortium (v4.1.1). The human body model is a conglomeration of several validated models of different sections of the body. Of particular interest for this research is the Wayne State University Head Injury Model for modeling the brain. These human body models were validated using a combination of cadaveric and animal studies. In this study, the football helmet was validated by laboratory testing using drop tests on the crown of the helmet. By coupling the two models into one finite element model, the brain response to impact loads caused by helmet design features can be investigated. In the present research, LS-DYNA is used to study a helmet crown impact with a rigid steel plate so as to obtain the strain-rate, strain, and stress experienced in the corpus callosum, midbrain, and brain stem as these anatomical regions are areas of concern with respect to mTBI.
ContributorsDarling, Timothy (Author) / Rajan, Subramaniam D. (Thesis advisor) / Muthuswamy, Jitendran (Thesis advisor) / Oswald, Jay (Committee member) / Mignolet, Marc (Committee member) / Arizona State University (Publisher)
Created2014
Description
Monte Carlo simulations are traditionally carried out for the determination of the amplification of forced vibration response of turbomachine/jet engine blades to mistuning. However, this effort can be computationally time consuming even when using the various reduced order modeling techniques. Accordingly, some investigations in the past have focused on obtaining simple approximate estimates for this amplification. In particular, two of these have proposed the use of harmonic patterns of the blade properties around the disk as an approximate alternative to the many random patterns of Monte Carlo analyses. These investigations, while quite encouraging, have relied solely on single degree of freedom per sector models of the rotor.
In this light, the overall focus of the present effort is a revisit of harmonic
mistuning of rotors focusing first the confirmation of the previously obtained findings with a more detailed model of the blisk in both conditions of an isolated blade-dominated resonance and of a veering between blade and disk dominated modes. The latter condition cannot be simulated by a single degree of freedom per sector model. Further, the analysis will consider the distinct cases of mistuning due to variations of material properties (Young's modulus) and geometric properties (geometric mistuning). In the single degree of freedom model, both mistuning types are equivalent but they are not, as demonstrated here, in more realistic models. The difference arises because changes in geometry induce not only changes in natural frequencies of the blades alone but of their modes and the importance of these two sources of variability is discussed with both Monte Carlo simulation and harmonic mistuning results.
The present investigation focuses also on the possible extension of the harmonic mistuning concept and of its quantitative information that can be derived from such analyses. From it, a novel measure of blade-disk coupling is introduced and assessed in comparison with the coupling index introduced in the past. In conclusions, the low cost of harmonic mistuning computations in comparison with full Monte Carlo simulations is
demonstrated to be worthwhile to elucidate the basic behavior of the mistuned rotor in a random setting.
In this light, the overall focus of the present effort is a revisit of harmonic
mistuning of rotors focusing first the confirmation of the previously obtained findings with a more detailed model of the blisk in both conditions of an isolated blade-dominated resonance and of a veering between blade and disk dominated modes. The latter condition cannot be simulated by a single degree of freedom per sector model. Further, the analysis will consider the distinct cases of mistuning due to variations of material properties (Young's modulus) and geometric properties (geometric mistuning). In the single degree of freedom model, both mistuning types are equivalent but they are not, as demonstrated here, in more realistic models. The difference arises because changes in geometry induce not only changes in natural frequencies of the blades alone but of their modes and the importance of these two sources of variability is discussed with both Monte Carlo simulation and harmonic mistuning results.
The present investigation focuses also on the possible extension of the harmonic mistuning concept and of its quantitative information that can be derived from such analyses. From it, a novel measure of blade-disk coupling is introduced and assessed in comparison with the coupling index introduced in the past. In conclusions, the low cost of harmonic mistuning computations in comparison with full Monte Carlo simulations is
demonstrated to be worthwhile to elucidate the basic behavior of the mistuned rotor in a random setting.
ContributorsSahoo, Saurav (Author) / Mignolet, Marc Paul (Thesis advisor) / Chattopadhyay, Aditi (Committee member) / Oswald, Jay (Committee member) / Arizona State University (Publisher)
Created2014
Description
The focus of this investigation is on the formulation and a validation of reduced order models (ROMs) for the prediction of the response of structures with embedded piezoelectric actuators. The ROMs considered here are those constructed in a nonintrusive manner from a commercial finite element software, NASTRAN is adopted here. Notwithstanding the popularity of piezoelectric materials in structural dynamics related applications such as structural health monitoring and energy harvesting, not all commercial finite element software allow directly their modeling. In such cases, e.g., with NASTRAN, one can proceed with an analogy and replace the electric actuation in the piezoelectric material by a fictitious thermal effect producing the same strain. This process recasts the determination of a ROM for a structure with embedded piezoelectric actuator into a similar ROM but for a heated structure, the framework of which has recently been developed. Yet, the temperature field resulting from the analogy would be quite different from the one considered in past effort and would excite a broad array of structural modes. Accordingly, as a preamble to considering a beam with a piezoelectric layer, a simpler plate model is considered that is subjected to a uniform temperature but a complex pressure loading that excites the entire set of modes of the plate in the broad frequency band considered. The very good match of the predictions obtained by this ROM in comparison to their full finite element counterparts provides the necessary confidence to next address a beam with embedded piezoelectic actuator. The test model considered for this validation is a built-up nano beam analyzed recently in nonlinear geometric conditions by full finite elements and by a non-intrusive ROM procedure under harmonic variations of the piezoelectic voltage. This structural model and its loading conditions are very different from those considered in past applications of nonintrusive ROMs, thus the excellent results obtained here provide further support of the broad generality of the nonintrusive ROM methodology, including of the appropriateness of the "dual modes" basis functions.
ContributorsVyas, Varun (Author) / Mignolet, Marc (Thesis advisor) / Hollkamp, Joseph (Committee member) / Liu, Yongming (Committee member) / Arizona State University (Publisher)
Created2014
Description
Fission products in nuclear fuel pellets can affect fuel performance as they change the fuel chemistry and structure. The behavior of the fission products and their release mechanisms are important to the operation of a power reactor. Research has shown that fission product release can occur through grain boundary (GB) at low burnups. Early fission gas release models, which assumed spherical grains with no effect of GB diffusion, did not capture the early stage of the release behavior well. In order to understand the phenomenon at low burnup and how it leads to the later release mechanism, a microstructurally explicit model is needed. This dissertation conducted finite element simulations of the transport behavior using 3-D microstructurally explicit models. It looks into the effects of GB character, with emphases on conditions that can lead to enhanced effective diffusion. Moreover, the relationship between temperature and fission product transport is coupled to reflect the high temperature environment.
The modeling work began with 3-D microstructure reconstruction for three uranium oxide samples with different oxygen stoichiometry: UO2.00 UO2.06 and UO2.14. The 3-D models were created based on the real microstructure of depleted UO2 samples characterized by Electron Backscattering Diffraction (EBSD) combined with serial sectioning. Mathematical equations on fission gas diffusion and heat conduction were studied and derived to simulate the fission gas transport under GB effect. Verification models showed that 2-D elements can be used to model GBs to reduce the number of elements. The effect of each variable, including fuel stoichiometry, temperature, GB diffusion, triple junction diffusion and GB thermal resistance, is verified, and they are coupled in multi-physics simulations to study the transport of fission gas at different radial location of a fuel pellet. It was demonstrated that the microstructural model can be used to incorporate the effect of different physics to study fission gas transport. The results suggested that the GB effect is the most significant at the edge of fuel pellet where the temperature is the lowest. In the high temperature region, the increase in bulk diffusivity due to excess oxygen diminished the effect of GB diffusion.
The modeling work began with 3-D microstructure reconstruction for three uranium oxide samples with different oxygen stoichiometry: UO2.00 UO2.06 and UO2.14. The 3-D models were created based on the real microstructure of depleted UO2 samples characterized by Electron Backscattering Diffraction (EBSD) combined with serial sectioning. Mathematical equations on fission gas diffusion and heat conduction were studied and derived to simulate the fission gas transport under GB effect. Verification models showed that 2-D elements can be used to model GBs to reduce the number of elements. The effect of each variable, including fuel stoichiometry, temperature, GB diffusion, triple junction diffusion and GB thermal resistance, is verified, and they are coupled in multi-physics simulations to study the transport of fission gas at different radial location of a fuel pellet. It was demonstrated that the microstructural model can be used to incorporate the effect of different physics to study fission gas transport. The results suggested that the GB effect is the most significant at the edge of fuel pellet where the temperature is the lowest. In the high temperature region, the increase in bulk diffusivity due to excess oxygen diminished the effect of GB diffusion.
ContributorsLim, Harn Chyi (Author) / Peralta, Pedro (Thesis advisor) / Jiang, Hanqing (Committee member) / Solanki, Kiran (Committee member) / Arizona State University (Publisher)
Created2014