ASU Electronic Theses and Dissertations
This collection includes most of the ASU Theses and Dissertations from 2011 to present. ASU Theses and Dissertations are available in downloadable PDF format; however, a small percentage of items are under embargo. Information about the dissertations/theses includes degree information, committee members, an abstract, supporting data or media.
In addition to the electronic theses found in the ASU Digital Repository, ASU Theses and Dissertations can be found in the ASU Library Catalog.
Dissertations and Theses granted by Arizona State University are archived and made available through a joint effort of the ASU Graduate College and the ASU Libraries. For more information or questions about this collection contact or visit the Digital Repository ETD Library Guide or contact the ASU Graduate College at gradformat@asu.edu.
Displaying 1 - 10 of 39
Filtering by
- Language: English
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
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
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
In this dissertation, the results of our comprehensive computational studies of disordered jammed (i.e., mechanically stable) packings of hard particles are presented, including the family of superdisks in 2D and ellipsoids in 3D Euclidean space. Following a very brief introduction to the hard-particle systems, the event driven molecular dynamics (EDMD) employed to generate the packing ensembles will be discussed. A large number of 2D packing configurations of superdisks are subsequently analyzed, through which a relatively accurate theoretical scheme for packing-fraction prediction based on local particle contact configurations is proposed and validated via additional numerical simulations. Moreover, the studies on binary ellipsoid packing in 3D are briefly discussed and the effects of different geometrical parameters on the final packing fraction are analyzed.
ContributorsXu, Yaopengxiao (Author) / Jiao, Yang (Thesis advisor) / Oswald, Jay (Committee member) / Liu, Yongming (Committee member) / Arizona State University (Publisher)
Created2014
Description
This research examines several critical aspects of the so-called "film induced cleavage" model of stress corrosion cracking using silver-gold alloys as the parent-phase material. The model hypothesizes that the corrosion generates a brittle nanoporous film, which subsequently fractures forming a high-speed crack that is injected into the uncorroded parent-phase alloy. This high speed crack owing to its kinetic energy can penetrate beyond the corroded layer into the parent phase and thus effectively reducing strength of the parent phase. Silver-gold alloys provide an ideal system to study this effect, as hydrogen effect can be ruled out on thermodynamic basis. During corrosion of the silver-gold alloy, the less noble metal i.e. silver is removed from the system leaving behind a nanoporous gold (NPG) layer. In the case of polycrystalline material, this corrosion process proceeds deeper along the grain boundary than the matrix grain. All of the cracks with apparent penetration beyond the corroded (dealloyed) layer are intergranular. Our aim was to study the crack penetration depth along the grain boundary to ascertain whether the penetration occurs past the grain-boundary dealloyed depth. EDS and imaging in high-resolution aberration corrected scanning transmission electron microscope (STEM) and atom probe tomography (APT) have been used to evaluate the grain boundary corrosion depth.
The mechanical properties of monolithic NPG are also studied. The motivation behind this is two-fold. The crack injection depth depends on the speed of the crack formed in the nanoporous layer, which in turn depends on the mechanical properties of the NPG. Also NPG has potential applications in actuation, sensing and catalysis. The measured value of the Young's modulus of NPG with 40 nm ligament size and 28% density was ~ 2.5 GPa and the Poisson's ratio was ~ 0.20. The fracture stress was observed to be ~ 11-13 MPa. There was no significant change observed between these mechanical properties on oxidation of NPG at 1.4 V. The fracture toughness value for the NPG was ~ 10 J/m2. Also dynamic fracture tests showed that the NPG is capable of supporting crack velocities ~ 100 - 180 m/s.
The mechanical properties of monolithic NPG are also studied. The motivation behind this is two-fold. The crack injection depth depends on the speed of the crack formed in the nanoporous layer, which in turn depends on the mechanical properties of the NPG. Also NPG has potential applications in actuation, sensing and catalysis. The measured value of the Young's modulus of NPG with 40 nm ligament size and 28% density was ~ 2.5 GPa and the Poisson's ratio was ~ 0.20. The fracture stress was observed to be ~ 11-13 MPa. There was no significant change observed between these mechanical properties on oxidation of NPG at 1.4 V. The fracture toughness value for the NPG was ~ 10 J/m2. Also dynamic fracture tests showed that the NPG is capable of supporting crack velocities ~ 100 - 180 m/s.
ContributorsBadwe, Nilesh (Author) / Sieradzki, Karl (Thesis advisor) / Peralta, Pedro (Committee member) / Oswald, Jay (Committee member) / Mahajan, Ravi (Committee member) / Arizona State University (Publisher)
Created2014
Description
Al 7075 alloys are used in a variety of structural applications, such as aircraft wings, automotive components, fuselage, spacecraft, missiles, etc. The mechanical and corrosion behavior of these alloys are dependent on their microstructure and the environment. Therefore, a comprehensive study on microstructural characterization and stress-environment interaction is necessary. Traditionally, 2D techniques have been used to characterize microstructure, which are inaccurate and inadequate since the research has shown that the results obtained in the bulk are different from those obtained on the surface. There now exist several techniques in 3D, which can be used to characterize the microstructure. Al 7075 alloys contain second phase particles which can be classified as Fe-bearing inclusions, Si-bearing inclusions and precipitates. The variation in mechanical and corrosion properties of aluminum alloys has been attributed to the size, shape, distribution, corrosion properties and mechanical behavior of these precipitates and constituent particles. Therefore, in order to understand the performance of Al 7075 alloys, it is critical to investigate the size and distribution of inclusions and precipitates in the alloys along with their mechanical properties, such as Young's modulus, hardness and stress-strain behavior. X-ray tomography and FIB tomography were used to visualize and quantify the microstructure of constituent particles (inclusions) and precipitates, respectively. Microscale mechanical characterization techniques, such as nanoindentation and micropillar compression, were used to obtain mechanical properties of inclusions. Over the years, studies have used surface measurements to understand corrosion behavior of materials. More recently, in situ mechanical testing has become more attractive and advantageous, as it enables visualization and quantification of microstructural changes as a function of time (4D). In this study, in situ X-ray synchrotron tomography was used to study the SCC behavior of Al 7075 alloys in moisture and deionized water. Furthermore, experiments were performed in EXCO solution to study the effect of applied stress on exfoliation behavior in 3D. Contrary to 2D measurements made at the surface which suggest non-uniform crack growth rates, three dimensional measurements of the crack length led to a much more accurate measurement of crack growth rates.
ContributorsSingh, Sudhanshu Shekhar (Author) / Chawla, Nikhilesh (Thesis advisor) / Alford, Terry (Committee member) / Solanki, Kiran (Committee member) / Oswald, Jay (Committee member) / Arizona State University (Publisher)
Created2015
Description
Fracture phenomena have been extensively studied in the last several decades. Continuum mechanics-based approaches, such as finite element methods and extended finite element methods, are widely used for fracture simulation. One well-known issue of these approaches is the stress singularity resulted from the spatial discontinuity at the crack tip/front. The requirement of guiding criteria for various cracking behaviors, such as initiation, propagation, and branching, also poses some challenges. Comparing to the continuum based formulation, the discrete approaches, such as lattice spring method, discrete element method, and peridynamics, have certain advantages when modeling various fracture problems due to their intrinsic characteristics in modeling discontinuities.
A novel, alternative, and systematic framework based on a nonlocal lattice particle model is proposed in this study. The uniqueness of the proposed model is the inclusion of both pair-wise local and multi-body nonlocal potentials in the formulation. First, the basic ideas of the proposed framework for 2D isotropic solid are presented. Derivations for triangular and square lattice structure are discussed in detail. Both mechanical deformation and fracture process are simulated and model verification and validation are performed with existing analytical solutions and experimental observations. Following this, the extension to general 3D isotropic solids based on the proposed local and nonlocal potentials is given. Three cubic lattice structures are discussed in detail. Failure predictions using the 3D simulation are compared with experimental testing results and very good agreement is observed. Next, a lattice rotation scheme is proposed to account for the material orientation in modeling anisotropic solids. The consistency and difference compared to the classical material tangent stiffness transformation method are discussed in detail. The implicit and explicit solution methods for the proposed lattice particle model are also discussed. Finally, some conclusions and discussions based on the current study are drawn at the end.
A novel, alternative, and systematic framework based on a nonlocal lattice particle model is proposed in this study. The uniqueness of the proposed model is the inclusion of both pair-wise local and multi-body nonlocal potentials in the formulation. First, the basic ideas of the proposed framework for 2D isotropic solid are presented. Derivations for triangular and square lattice structure are discussed in detail. Both mechanical deformation and fracture process are simulated and model verification and validation are performed with existing analytical solutions and experimental observations. Following this, the extension to general 3D isotropic solids based on the proposed local and nonlocal potentials is given. Three cubic lattice structures are discussed in detail. Failure predictions using the 3D simulation are compared with experimental testing results and very good agreement is observed. Next, a lattice rotation scheme is proposed to account for the material orientation in modeling anisotropic solids. The consistency and difference compared to the classical material tangent stiffness transformation method are discussed in detail. The implicit and explicit solution methods for the proposed lattice particle model are also discussed. Finally, some conclusions and discussions based on the current study are drawn at the end.
ContributorsChen, Hailong (Author) / Liu, Yongming (Thesis advisor) / Jiao, Yang (Committee member) / Mignolet, Marc (Committee member) / Oswald, Jay (Committee member) / Solanki, Kiran (Committee member) / Arizona State University (Publisher)
Created2015
Description
In this thesis, a FORTRAN code is rewritten in C++ with an object oriented ap-
proach. There are several reasons for this purpose. The first reason is to establish
the basis of a GPU programming. To write programs that utilize GPU hardware,
CUDA or OpenCL is used which only support C and C++. FORTRAN has a feature
that lets its programs to call C/C++ functions. FORTRAN sends relevant data to
C/C++, which in turn sends that data to OpenCL. Although this approach works,
it makes the code messy and bulky and in the end more difficult to deal with. More-
over, there is a slight performance decrease from the additional data copy. This is
the motivation to have the code entirely written in C++ to make it more uniform,
efficient and clean. The second reason is the object oriented feature of the C++. The
“abstraction”, “inheritance” and “run-time polymorphism” features of C++ provide
some form of classes and objects, the ability to build new abstractions, and some
form of run-time binding, respectively. In recent years, some of popular codes has
been rewritten in C++ which were initially in FORTRAN. One of these softwares is
LAMMPS.
In this code the level set equation is solved by RLSG method to track the interface in
two phase flow. In gas/fluid flows, the surface tension is important and only exists at
the interface. Therefore, the location and some geometric features of interface need
to be evaluated which can be achieved by solving the level set equation.
proach. There are several reasons for this purpose. The first reason is to establish
the basis of a GPU programming. To write programs that utilize GPU hardware,
CUDA or OpenCL is used which only support C and C++. FORTRAN has a feature
that lets its programs to call C/C++ functions. FORTRAN sends relevant data to
C/C++, which in turn sends that data to OpenCL. Although this approach works,
it makes the code messy and bulky and in the end more difficult to deal with. More-
over, there is a slight performance decrease from the additional data copy. This is
the motivation to have the code entirely written in C++ to make it more uniform,
efficient and clean. The second reason is the object oriented feature of the C++. The
“abstraction”, “inheritance” and “run-time polymorphism” features of C++ provide
some form of classes and objects, the ability to build new abstractions, and some
form of run-time binding, respectively. In recent years, some of popular codes has
been rewritten in C++ which were initially in FORTRAN. One of these softwares is
LAMMPS.
In this code the level set equation is solved by RLSG method to track the interface in
two phase flow. In gas/fluid flows, the surface tension is important and only exists at
the interface. Therefore, the location and some geometric features of interface need
to be evaluated which can be achieved by solving the level set equation.
ContributorsSafarkhani, Salar (Author) / Herrmann, Mrcus (Thesis advisor) / Oswald, Jay (Committee member) / Rykczewski, Konrad (Committee member) / Arizona State University (Publisher)
Created2015