Matching Items (40)
Filtering by

Clear all filters

152030-Thumbnail Image.png
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

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
151523-Thumbnail Image.png
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

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
153438-Thumbnail Image.png
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

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.
ContributorsPanwar, Ajay (Author) / Oswald, Jay (Thesis advisor) / Dooley, Kevin (Committee member) / Chen, Kangping (Committee member) / Arizona State University (Publisher)
Created2015
153495-Thumbnail Image.png
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

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.
ContributorsSong, Pengchao (Author) / Mignolet, Marc P (Thesis advisor) / Chattopadhyay, Aditi (Committee member) / Oswald, Jay (Committee member) / Arizona State University (Publisher)
Created2015
153325-Thumbnail Image.png
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

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.
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
153327-Thumbnail Image.png
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

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.
ContributorsSahoo, Saurav (Author) / Mignolet, Marc Paul (Thesis advisor) / Chattopadhyay, Aditi (Committee member) / Oswald, Jay (Committee member) / Arizona State University (Publisher)
Created2014
153180-Thumbnail Image.png
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 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.
ContributorsBadwe, Nilesh (Author) / Sieradzki, Karl (Thesis advisor) / Peralta, Pedro (Committee member) / Oswald, Jay (Committee member) / Mahajan, Ravi (Committee member) / Arizona State University (Publisher)
Created2014
154124-Thumbnail Image.png
Description
The objective of this research is to develop robust, accurate, and adaptive algorithms in the framework of the extended finite element method (XFEM) for fracture analysis of highly heterogeneous materials with complex internal geometries. A key contribution of this work is the creation of novel methods designed to automate the

The objective of this research is to develop robust, accurate, and adaptive algorithms in the framework of the extended finite element method (XFEM) for fracture analysis of highly heterogeneous materials with complex internal geometries. A key contribution of this work is the creation of novel methods designed to automate the incorporation of high-resolution data, e.g. from X-ray tomography, that can be used to better interpret the enormous volume of data generated in modern in-situ experimental testing. Thus new algorithms were developed for automating analysis of complex microstructures characterized by segmented tomographic images.

A centrality-based geometry segmentation algorithm was developed to accurately identify discrete inclusions and particles in composite materials where limitations in imaging resolution leads to spurious connections between particles in close contact.To allow for this algorithm to successfully segment geometry independently of particle size and shape, a relative centrality metric was defined to allow for a threshold centrality criterion for removal of voxels that spuriously connect distinct geometries.

To automate incorporation of microstructural information from high-resolution images, two methods were developed that initialize signed distance fields on adaptively-refined finite element meshes. The first method utilizes a level set evolution equation that is directly solved on the finite element mesh through Galerkins method. The evolution equation is formulated to produce a signed distance field that matches geometry defined by a set of voxels segmented from tomographic images. The method achieves optimal convergence for the order of elements used. In a second approach, the fast marching method is employed to initialize a distance field on a uniform grid which is then projected by least squares onto a finite element mesh. This latter approach is shown to be superior in speed and accuracy.

Lastly, extended finite element method simulations are performed for the analysis of particle fracture in metal matrix composites with realistic particle geometries initialized from X-ray tomographic data. In the simulations, particles fracture probabilistically through a Weibull strength distribution. The model is verified through comparisons with the experimentally-measured stress-strain response of the material as well as analysis of the fracture. Further, simulations are then performed to analyze the effect of mesh sensitivity, the effect of fracture of particles on their neighbors, and the role of a particles shape on its fracture probability.
ContributorsYuan, Rui (Author) / Oswald, Jay (Thesis advisor) / Chawla, Nikhilesh (Committee member) / Liu, Yongming (Committee member) / Solanki, Kiran (Committee member) / Chen, Kangping (Committee member) / Arizona State University (Publisher)
Created2015
155972-Thumbnail Image.png
Description
The exceptional mechanical properties of polymers with heterogeneous structure, such as the high toughness of polyethylene and the excellent blast-protection capability of polyurea, are strongly related to their morphology and nanoscale structure. Different polymer microstructures, such as semicrystalline morphology and segregated nanophases, lead to coordinated molecular motions during deformation

The exceptional mechanical properties of polymers with heterogeneous structure, such as the high toughness of polyethylene and the excellent blast-protection capability of polyurea, are strongly related to their morphology and nanoscale structure. Different polymer microstructures, such as semicrystalline morphology and segregated nanophases, lead to coordinated molecular motions during deformation in order to preserve compatibility between the different material phases. To study molecular relaxation in polyethylene, a coarse-grained model of polyethylene was calibrated to match the local structural variable distributions sampled from supercooled atomistic melts. The coarse-grained model accurately reproduces structural properties, e.g., the local structure of both the amorphous and crystalline phases, and thermal properties, e.g., glass transition and melt temperatures, and dynamic properties: including the vastly different relaxation time scales of the amorphous and crystalline phases. A hybrid Monte Carlo routine was developed to generate realistic semicrystalline configurations of polyethylene. The generated systems accurately predict the activation energy of the alpha relaxation process within the crystalline phase. Furthermore, the models show that connectivity to long chain segments in the amorphous phase increases the energy barrier for chain slip within crystalline phase. This prediction can guide the development of tougher semicrystalline polymers by providing a fundamental understanding of how nanoscale morphology contributes to chain mobility. In a different study, the macroscopic shock response of polyurea, a phase segregated copolymer, was analyzed using density functional theory (DFT) molecular dynamics (MD) simulations and classical MD simulations. The two models predict the shock response consistently up to shock pressures of 15 GPa, beyond which the DFT-based simulations predict a softer response. From the DFT simulations, an analysis of bond scission was performed as a first step in developing a more fundamental understanding of how shock induced material transformations effect the shock response and pressure dependent strength of polyurea subjected to extreme shocks.
ContributorsLi, Yiyang (Author) / Oswald, Jay (Thesis advisor) / Rajan, Subramaniam D. (Committee member) / Solanki, Kiran (Committee member) / Chamberlin, Ralph (Committee member) / Rajagopalan, Jagannathan (Committee member) / Arizona State University (Publisher)
Created2017
153841-Thumbnail Image.png
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

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.
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