Matching Items (22)
Filtering by
- All Subjects: Mechanical Engineering
- Genre: Doctoral Dissertation
- Creators: Oswald, Jay
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
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
An approach for modeling resistance spot welding of thin-gauge, dissimilar metal sheets with high electrical conductivity is presented in this work. In this scenario, the electrical and thermal contact resistances play a dominant role in heat generation and temperature evolution within the workpieces; these interactions ultimately control the weld geometry. Existing models are limited in modeling these interactions, especially for dissimilar and thin-gauge metal sheets, and at higher temperatures when the multiphysics becomes increasingly interdependent. The approach presented here uses resistivity measurements, combined with thermal modeling and known bulk resistance relationships to infer the relationship between electrical contact resistance and temperature for each of the different material interfaces in the welding process. Corresponding thermal contact resistance models are developed using the Wiedemann-Franz law combined with a scaling factor to account for nonmetallic behavior. Experimental and simulation voltage histories and final weld diameter were used to validate this model for a Cu/Al/Cu and a Cu/Al/Cu/Al/Cu stack-ups. This model was then used to study the effect of Ni-P coating on resistance spot welding of Cu and Al sheets in terms of weld formation, mechanical deformation, and contact resistance. Contact resistance and current density distribution are highly dependent on contact pressure and temperature distribution at the Cu/Al interface in the presence of alumina. The Ni-P coating helps evolve a partially-bonded donut shaped weld into a fully-bonded hourglass-shaped weld by decreasing the dependence of contact resistance and current density distribution on contact pressure and temperature distribution at the Cu/Al interface. This work also provides an approach to minimize distortion due to offset-rolling in thin aluminum sheets by optimizing the stiffening feature geometry. The distortion is minimized using particle swarm optimization. The objective function is a function of distortion and smallest radius of curvature in the geometry. Doubling the minimum allowable radius of curvature nearly doubles the reduction in distortion from the stadium shape for a quarter model. Reduction in distortion in the quarter model extends to the full-scale model with the best design performing 5.3% and 27% better than the corresponding nominal design for a quarter and full-scale model, respectively.
ContributorsVeeresh, Pawan (Author) / Oswald, Jay (Thesis advisor) / Carlson, Blair (Committee member) / Hoover, Christian (Committee member) / Rajagopalan, Jagannathan (Committee member) / Solanki, Kiran (Committee member) / Arizona State University (Publisher)
Created2022
Description
Aromatic polymers, with benzene-like rings in their main chains, include materials such as polyurea, an amorphous elastomer capable of dissipating large amounts of energy under dynamic loading, which makes it a promising coating for defensive systems. Although computational research exists that investigates the atomic-level response of polyurea and other amorphous aromatic polymers to extreme conditions, there is little experimental work to validate these models 1) at the atomic-scale and 2) under high pressures characteristic of extreme dynamic loading. Understanding structure-property relationships at the atomic-level is important for polymers, considering many of them undergo pressure and temperature-induced structural transformations, which must be understood to formulate accurate predictive models. This work aims to gain a deeper understanding of the high-pressure structural response of aromatic polymers at the atomic-level, with emphasis into the mechanisms associated with high-pressure transformations. Hence, atomic-level structural data at high pressures was obtained in situ via multiangle energy dispersive X-ray diffraction (EDXD) experiments at the Advanced Photon Source (APS) for polyurea and another amorphous aromatic polymer, polysulfone, chosen as a reference due to its relatively simple structure. Pressures up to 6 GPa were applied using a Paris Edinburgh (PE) hydraulic press at room temperature. Select polyurea samples were also heated to 277 °C at 6 GPa. The resulting structure factors and pair distribution functions, along with molecular dynamics simulations of polyurea provided by collaborators, suggest that the structures of both polymers are stable up to 6 GPa, aside from reductions in free-volume between polymer backbones. As higher pressures (≲ 32 GPa) were applied using diamond anvils in combination with the PE press, indications of structural transformations were observed in both polymers that appear similar in nature to the sp2-sp3 hybridization in compressed carbon. The transformation occurs gradually up to at least ~ 26 GPa in PSF, while it does not progress past ~ 15 GPa in polyurea. The changes are largely reversible, especially in polysulfone, consistent with pressure-driven, reversible graphite-diamond transformations in the absence of applied temperature. These results constitute some of the first in situ observations of the mechanisms that drive pressure-induced structural transformations in aromatic polymers.
ContributorsEastmond, Tyler (Author) / Peralta, Pedro (Thesis advisor) / Hoover, Christian (Committee member) / Hrubiak, Rostislav (Committee member) / Mignolet, Marc (Committee member) / Oswald, Jay (Committee member) / Arizona State University (Publisher)
Created2022
Description
Intelligent engineering designs require an accurate understanding of material behavior, since any uncertainties or gaps in knowledge must be counterbalanced with heightened factors of safety, leading to overdesign. Therefore, building better structures and pushing the performance of new components requires an improved understanding of the thermomechanical response of advanced materials under service conditions. This dissertation provides fundamental investigations of several advanced materials: thermoset polymers, a common matrix material for fiber-reinforced composites and nanocomposites; aluminum alloy 7075-T6 (AA7075-T6), a high-performance aerospace material; and ceramic matrix composites (CMCs), an advanced composite for extreme-temperature applications. To understand matrix interactions with various interfaces and nanoinclusions at their fundamental scale, the properties of thermoset polymers are studied at the atomistic scale. An improved proximity-based molecular dynamics (MD) technique for modeling the crosslinking of thermoset polymers is carefully established, enabling realistic curing simulations through its ability to dynamically and probabilistically perform complex topology transformations. The proximity-based MD curing methodology is then used to explore damage initiation and the local anisotropic evolution of mechanical properties in thermoset polymers under uniaxial tension with an emphasis on changes in stiffness through a series of tensile loading, unloading, and reloading experiments. Aluminum alloys in aerospace applications often require a fatigue life of over 109 cycles, which is well over the number of cycles that can be practically tested using conventional fatigue testing equipment. In order to study these high-life regimes, a detailed ultrasonic cycle fatigue study is presented for AA7075-T6 under fully reversed tension-compression loading. The geometric sensitivity, frequency effects, size effects, surface roughness effects, and the corresponding failure mechanisms for ultrasonic fatigue across different fatigue regimes are investigated. Finally, because CMCs are utilized in extreme environments, oxidation plays an important role in their degradation. A multiphysics modeling methodology is thus developed to address the complex coupling between oxidation, mechanical stress, and oxygen diffusion in heterogeneous carbon fiber-reinforced CMC microstructures.
ContributorsSchichtel, Jacob (Author) / Chattopadhyay, Aditi (Thesis advisor) / Dai, Lenore (Committee member) / Ghoshal, Anindya (Committee member) / Huang, Huei-Ping (Committee member) / Jiao, Yang (Committee member) / Oswald, Jay (Committee member) / Arizona State University (Publisher)
Created2022
Description
Over the past few decades there has been significant interest in the design and construction of hypersonic vehicles. Such vehicles exhibit strongly coupled aerodynamics, acoustics, heat transfer, and structural deformations, which can take significant computational efforts to simulate using standard finite element and computational fluid dynamics techniques. This situation has lead to development of various reduced order modelling (ROM) methods which reduce the parameter space of these simulations so they can be run more quickly. The planned hypersonic vehicles will be constructed by assembling a series of sub-structures, such as panels and stiffeners, that will be welded together creating built-up structures.In this light, the focus of the present investigation is on the formulation and validation of nonlinear reduced order models (NLROMs) of built-up structures that include nonlinear geometric effects induced by the large loads/large response. Moreover, it is recognized that gaps between sub-structures could result from the these intense loadings can thus the inclusion of the nonlinearity introduced by contact separation will also be addressed. These efforts, application to built-up structures and inclusion of contact nonlinearity, represent novel developments of existing NLROM strategies. A hat stiffened panel is selected as a representative example of built-up structure and a compact NRLOM is successfully constructed for this structure which exhibited a potential internal resonance. For the investigation of contact nonlinearity, two structural models were used: a cantilevered beam which can contact several stops and an overlapping plate model which can exhibit the opening/closing of a gap. Successful NLROMs were constructed for these structures with the basis for the plate model determined as a two-step process, i.e., considering the plate without gap first and then enriching the corresponding basis to account for opening of the gap. Adaptions were then successfully made to a Newton-Raphson solver to properly account for contact and the associated forces in static predictions by NLROMs.
ContributorsWainwright, Bret Aaron (Author) / Mignolet, Marc P (Thesis advisor) / Oswald, Jay (Committee member) / Peralta, Pedro (Committee member) / Spottswood, Stephen (Committee member) / Rajan, Subramaniam (Committee member) / Arizona State University (Publisher)
Created2021
Description
Owing to the surge in development of endovascular devices such as coils and flow diverter stents, doctors are inclined to approach surgical cases non-invasively more often than before. Treating brain aneurysms as a bulging of a weakened area of a blood vessel is no exception. Therefore, promoting techniques that can help surgeons have a better idea of treatment outcomes are of invaluable importance.
In order to investigate the effects of these devices on intra-aneurysmal hemodynamics, the conventional computational fluid dynamics (CFD) approach uses the explicit geometry of the device within an aneurysm and discretizes the fluid domain to solve the Navier-Stokes equations. However, since the devices are made of small struts, the number of mesh elements in the boundary layer region would be considerable. This cumbersome task led to the implementation of the porous medium assumption. In this approach, the explicit geometry of the device is eliminated, and relevant porous medium assumptions are applied. Unfortunately, as it will be shown in this research, some of the porous medium approaches used in the literature are over-simplified. For example, considering the porous domain to be homogeneous is one major drawback which leads to significant errors in capturing the intra-aneurysmal flow features. Specifically, since the devices must comply with the complex geometry of an aneurysm, the homogeneity assumption is not valid.
In this research, a novel heterogeneous porous medium approach is introduced. This results in a substantial reduction in the total number of mesh elements required to discretize the flow domain while not sacrificing the accuracy of the method by over-simplifying the utilized assumptions.
In order to investigate the effects of these devices on intra-aneurysmal hemodynamics, the conventional computational fluid dynamics (CFD) approach uses the explicit geometry of the device within an aneurysm and discretizes the fluid domain to solve the Navier-Stokes equations. However, since the devices are made of small struts, the number of mesh elements in the boundary layer region would be considerable. This cumbersome task led to the implementation of the porous medium assumption. In this approach, the explicit geometry of the device is eliminated, and relevant porous medium assumptions are applied. Unfortunately, as it will be shown in this research, some of the porous medium approaches used in the literature are over-simplified. For example, considering the porous domain to be homogeneous is one major drawback which leads to significant errors in capturing the intra-aneurysmal flow features. Specifically, since the devices must comply with the complex geometry of an aneurysm, the homogeneity assumption is not valid.
In this research, a novel heterogeneous porous medium approach is introduced. This results in a substantial reduction in the total number of mesh elements required to discretize the flow domain while not sacrificing the accuracy of the method by over-simplifying the utilized assumptions.
ContributorsYadollahi Farsani, Hooman (Author) / Herrmann, Marcus (Thesis advisor) / Frakes, David (Thesis advisor) / Chong, Brian (Committee member) / Peet, Yulia (Committee member) / Oswald, Jay (Committee member) / Arizona State University (Publisher)
Created2018
Description
Improved knowledge connecting the chemistry, structure, and properties of polymers is necessary to develop advanced materials in a materials-by-design approach. Molecular dynamics (MD) simulations can provide tremendous insight into how the fine details of chemistry, molecular architecture, and microstructure affect many physical properties; however, they face well-known restrictions in their applicable temporal and spatial scales. These limitations have motivated the development of computationally-efficient, coarse-grained methods to investigate how microstructural details affect thermophysical properties. In this dissertation, I summarize my research work in structure-based coarse-graining methods to establish the link between molecular-scale structure and macroscopic properties of two different polymers. Systematically coarse-grained models were developed to study the viscoelastic stress response of polyurea, a copolymer that segregates into rigid and viscous phases, at time scales characteristic of blast and impact loading. With the application of appropriate scaling parameters, the coarse-grained models can predict viscoelastic properties with a speed up of 5-6 orders of magnitude relative to the atomistic MD models. Coarse-grained models of polyethylene were also created to investigate the thermomechanical material response under shock loading. As structure-based coarse-grained methods are generally not transferable to states different from which they were calibrated at, their applicability for modeling non-equilibrium processes such as shock and impact is highly limited. To address this problem, a new model is developed that incorporates many-body interactions and is calibrated across a range of different thermodynamic states using a least square minimization scheme. The new model is validated by comparing shock Hugoniot properties with atomistic and experimental data for polyethylene. Lastly, a high fidelity coarse-grained model of polyethylene was constructed that reproduces the joint-probability distributions of structural variables such as the distributions of bond lengths and bond angles between sequential coarse-grained sites along polymer chains. This new model accurately represents the structure of both the amorphous and crystal phases of polyethylene and enabling investigation of how polymer processing such as cold-drawing and bulk crystallization affect material structure at significantly larger time and length scales than traditional molecular simulations.
ContributorsAgrawal, Vipin (Author) / Oswald, Jay (Thesis advisor) / Peralta, Pedro (Committee member) / Chamberlin, Ralph (Committee member) / Solanki, Kiran (Committee member) / Jiao, Yang (Committee member) / Arizona State University (Publisher)
Created2016
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
The objective of this research is to develop robust, accurate, and adaptive algorithms in the framework of the extended finite element method (XFEM) for fracture analysis of highly heterogeneous materials with complex internal geometries. A key contribution of this work is the creation of novel methods designed to automate the incorporation of high-resolution data, e.g. from X-ray tomography, that can be used to better interpret the enormous volume of data generated in modern in-situ experimental testing. Thus new algorithms were developed for automating analysis of complex microstructures characterized by segmented tomographic images.
A centrality-based geometry segmentation algorithm was developed to accurately identify discrete inclusions and particles in composite materials where limitations in imaging resolution leads to spurious connections between particles in close contact.To allow for this algorithm to successfully segment geometry independently of particle size and shape, a relative centrality metric was defined to allow for a threshold centrality criterion for removal of voxels that spuriously connect distinct geometries.
To automate incorporation of microstructural information from high-resolution images, two methods were developed that initialize signed distance fields on adaptively-refined finite element meshes. The first method utilizes a level set evolution equation that is directly solved on the finite element mesh through Galerkins method. The evolution equation is formulated to produce a signed distance field that matches geometry defined by a set of voxels segmented from tomographic images. The method achieves optimal convergence for the order of elements used. In a second approach, the fast marching method is employed to initialize a distance field on a uniform grid which is then projected by least squares onto a finite element mesh. This latter approach is shown to be superior in speed and accuracy.
Lastly, extended finite element method simulations are performed for the analysis of particle fracture in metal matrix composites with realistic particle geometries initialized from X-ray tomographic data. In the simulations, particles fracture probabilistically through a Weibull strength distribution. The model is verified through comparisons with the experimentally-measured stress-strain response of the material as well as analysis of the fracture. Further, simulations are then performed to analyze the effect of mesh sensitivity, the effect of fracture of particles on their neighbors, and the role of a particles shape on its fracture probability.
A centrality-based geometry segmentation algorithm was developed to accurately identify discrete inclusions and particles in composite materials where limitations in imaging resolution leads to spurious connections between particles in close contact.To allow for this algorithm to successfully segment geometry independently of particle size and shape, a relative centrality metric was defined to allow for a threshold centrality criterion for removal of voxels that spuriously connect distinct geometries.
To automate incorporation of microstructural information from high-resolution images, two methods were developed that initialize signed distance fields on adaptively-refined finite element meshes. The first method utilizes a level set evolution equation that is directly solved on the finite element mesh through Galerkins method. The evolution equation is formulated to produce a signed distance field that matches geometry defined by a set of voxels segmented from tomographic images. The method achieves optimal convergence for the order of elements used. In a second approach, the fast marching method is employed to initialize a distance field on a uniform grid which is then projected by least squares onto a finite element mesh. This latter approach is shown to be superior in speed and accuracy.
Lastly, extended finite element method simulations are performed for the analysis of particle fracture in metal matrix composites with realistic particle geometries initialized from X-ray tomographic data. In the simulations, particles fracture probabilistically through a Weibull strength distribution. The model is verified through comparisons with the experimentally-measured stress-strain response of the material as well as analysis of the fracture. Further, simulations are then performed to analyze the effect of mesh sensitivity, the effect of fracture of particles on their neighbors, and the role of a particles shape on its fracture probability.
ContributorsYuan, Rui (Author) / Oswald, Jay (Thesis advisor) / Chawla, Nikhilesh (Committee member) / Liu, Yongming (Committee member) / Solanki, Kiran (Committee member) / Chen, Kangping (Committee member) / Arizona State University (Publisher)
Created2015