Matching Items (24)
Filtering by
- All Subjects: Materials Science
- All Subjects: Finite element method
- 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
The focus of this investigation includes three aspects. First, the development of nonlinear reduced order modeling techniques for the prediction of the response of complex structures exhibiting "large" deformations, i.e. a geometrically nonlinear behavior, and modeled within a commercial finite element code. The present investigation builds on a general methodology, successfully validated in recent years on simpler panel structures, by developing a novel identification strategy of the reduced order model parameters, that enables the consideration of the large number of modes needed for complex structures, and by extending an automatic strategy for the selection of the basis functions used to represent accurately the displacement field. These novel developments are successfully validated on the nonlinear static and dynamic responses of a 9-bay panel structure modeled within Nastran. In addition, a multi-scale approach based on Component Mode Synthesis methods is explored. Second, an assessment of the predictive capabilities of nonlinear reduced order models for the prediction of the large displacement and stress fields of panels that have a geometric discontinuity; a flat panel with a notch was used for this assessment. It is demonstrated that the reduced order models of both virgin and notched panels provide a close match of the displacement field obtained from full finite element analyses of the notched panel for moderately large static and dynamic responses. In regards to stresses, it is found that the notched panel reduced order model leads to a close prediction of the stress distribution obtained on the notched panel as computed by the finite element model. Two enrichment techniques, based on superposition of the notch effects on the virgin panel stress field, are proposed to permit a close prediction of the stress distribution of the notched panel from the reduced order model of the virgin one. A very good prediction of the full finite element results is achieved with both enrichments for static and dynamic responses. Finally, computational challenges associated with the solution of the reduced order model equations are discussed. Two alternatives to reduce the computational time for the solution of these problems are explored.
ContributorsPerez, Ricardo Angel (Author) / Mignolet, Marc (Thesis advisor) / Oswald, Jay (Committee member) / Spottswood, Stephen (Committee member) / Peralta, Pedro (Committee member) / Jiang, Hanqing (Committee member) / Arizona State University (Publisher)
Created2012
Description
The goal of this research is to compare the mechanical properties of CP-Ti and Ti-O and to understand the relationship between a material's microstructure and its response to fatigue. Titanium has been selected due to its desirable properties and applicability in several engineering fields. Both samples are polished and etched in order to visualize and characterize the microstructure and its features. The samples then undergo strain-controlled fatigue tests for several thousand cycles. Throughout testing, images of the samples are taken at zero and maximum load for DIC analysis. The DIC results can be used to study the local strains of the samples. The DIC analysis performed on the CP-Ti sample and presented in this study will be used to understand how the addition of oxygen in the Ti-O impacts fatigue response. The outcome of this research can be used to develop long-lasting, high strength materials.
ContributorsRiley, Erin Ashland (Author) / Solanki, Kiran (Thesis director) / Oswald, Jay (Committee member) / School of Art (Contributor) / Mechanical and Aerospace Engineering Program (Contributor) / Barrett, The Honors College (Contributor)
Created2016-05
Description
This thesis presents a study of Boron Nitride (BN) and Copper (Cu)/BN multilayer thin films in terms of synthesis, chemical, structural, morphological, and mechanical properties characterization.
In this study, the influence of Ar/N₂ flow rate in synthesizing stoichiometric BN thin films via magnetron sputtering was investigated initially. Post magnetron sputtering, the crystalline nature and B:N stoichiometric ratio of deposited thin films were investigated by X-ray diffraction (XRD) and X-ray Photoelectron Spectroscopy (XPS) respectively. Thicknesses revealed by ellipsometry analysis for nearly stoichiometric B:N thin films and their corresponding deposition times were used for estimating BN interlayer deposition times during the deposition of Cu/BN multilayer thin films. To characterize the microstructure of the synthesized Cu/BN multilayer thin films, XRD and scanning electron microscopy (SEM) have been used. Finally, a comparison of nanoindentation measurements on pure Cu and Cu/BN multilayer thin films having different number of BN interlayers were used for studying the influence of BN interlayers on improving mechanical properties such as hardness and elastic modulus.
The results show that the stoichiometry of BN thin films is dependent on the Ar/N₂ flow rate during magnetron sputtering. An optimal Ar/N₂ flow rate of 13:5 during deposition was required to achieve an approximately 1:1 B:N stoichiometry. Grazing incidence and powder XRD analysis on these stoichiometric BN thin films deposited at room temperature did not reveal a phase match when compared to hexagonal boron nitride (h-BN) and cubic boron nitride (c-BN) reference XRD patterns. For a BN thin film deposition time of 5 hours, a thickness of approximately 40 nm was achieved, as revealed by ellipsometry. XRD and microstructure analysis using scanning electron microscopy (SEM) on pure Cu and Cu/BN thin films showed that the Cu grain size in Cu/BN thin films is much finer than pure Cu thin films. Interestingly, nanoindentation measurements on pure Cu and Cu/BN thin films having a similar overall thickness demonstrated that hardness and Young’s modulus of the films were improved significantly when BN interlayers are present.
ContributorsCaner, Sumeyye (Author) / Rajagopalan, Jagannathan (Thesis advisor) / Oswald, Jay (Committee member) / Solanki, Kiran (Committee member) / Arizona State University (Publisher)
Created2023
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
This research seeks to answer the question if there is a singular relationship between stishovite nucleation and the atomistic structure of the preshocked amorphous SiO$_2$. To do this a stishovite manufacturing method is developed in which 1,152 samples were produced. The majority of these samples did crystallize. The method was produced through two rounds of experiments and fine-tuning with the pressure damp, temperature damp, shock pressure using an NPHug fix, and sample origin. A new random atomic insertion method was used to generate a new and different SiO$_2$ amorphous structure not before seen within the research literature. The optimal values for shock were found to be 60~GPa for randomly atom insertion samples and 55~GPa for quartz origin samples. Temperature damp appeared to have a slight effect optimizing at 0.05~ps and the pressure damp had no visible effect, testing was done with temperature damp from 0.05 to 0.5~ps and pressure damp from 0.1 to 10.0~ps. There appeared to be significant randomness in crystallization behavior. The preshocked and postnucleated samples were transformed into Gaussian fields of crystal, mass, and charge. These fields were divided and classified using a cut-off method taking the number of crystals produced in portions of each simulation and classifying each potion as nucleated or non-nucleated. Data in which some nucleation but not a critical amount was present was removed constituting 2.6\% to 20.3\% of data in all tests. A max method was also used which takes only the maximum portions of each simulation to classify as nucleating. There are three other variables tested within this work, a sample size of 18,000 or 72,728~atoms, Gaussian variance of 1 or 4~\AA, and Convolutional neural network (CNN) architecture of a garden verity or all convolution along with the portioning classification method, sample origination, and Gaussian field type. In total 64 tests were performed to try every combination of variable. No significant classifications were made by the CNNs to nucleation or non-nucleation portions. The results clearly confirmed that the data was not abstracting to atomistic structure and was random by all classifications of the CNNs. The all convolution CNN testing did show smoother outcomes in training with less fluctuations. 59\% of all validation accuracy was held at 0.5 for a random state and 84\% was within $\pm0.02$ of 0.5. It is conclusive that prenucleation structure is not the sole predictor of nucleation behavior. It is not conclusive if prenucleation structure is a partial or non-factor within nucleation of stishovite from amorphous SiO$_2$.
ContributorsChristen, Jonathan Scorr (Author) / Oswald, Jay (Thesis advisor) / Muhich, Christopher (Committee member) / Zhuang, Houlong (Committee member) / Arizona State University (Publisher)
Created2021
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
Stress corrosion cracking (SCC) is a materials degradation phenomena resulting from a combination of stress and a corrosive environment. Among the alphabet soup of proposed mechanism of SCC the most important are film-rupture, film-induced cleavage and hydrogen embrittlement.
This work examines various aspects of film-induced cleavage in gold alloys for which the operation of hydrogen embrittlement processes can be strictly ruled out on thermodynamic grounds. This is so because in such alloys SCC occurs under electrochemical conditions within which water is stable to hydrogen gas evolution. The alloy system examined in this work is AgAu since the corrosion processes in this system occur by a dealloying mechanism that results in the formation of nanoporous gold. The physics behind the dealloying process as well as the resulting formation of nanoporous gold is today well understood.
Two important aspects of the film-induced cleavage mechanism are examined in this work: dynamic fracture in monolithic nanoporous gold and crack injection. In crack injection there is a finite thickness dealloyed layer formed on a AgAu alloy sample and the question of whether or not a crack that nucleates within this layer can travel for some finite distance into the un-corroded parent phase alloy is addressed. Dynamic fracture tests were performed on single edge-notched monolithic nanoporous gold samples as well as “infinite strip” sample configurations for which the stress intensity remains constant over a significant portion of the crack length. High-speed photography was used to measure the crack velocity. In the dynamic fracture experiments cracks were observed to travel at speeds as large as 270 m/s corresponding to about 68% of the Raleigh wave velocity. Crack injection experiments were performed on single crystal Ag77Au23, polycrystalline Ag72Au28 and pure gold, all of which had thin nanoporous gold layers on the surface of samples. Through-thickness fracture was seen in both the single crystal and polycrystalline samples and there was an indication of ~ 1 μm injected cracks into pure gold. These results have important implications for the operation of the film-induced cleavage mechanism and represent a first step in the development of a fundamental model of SCC.
This work examines various aspects of film-induced cleavage in gold alloys for which the operation of hydrogen embrittlement processes can be strictly ruled out on thermodynamic grounds. This is so because in such alloys SCC occurs under electrochemical conditions within which water is stable to hydrogen gas evolution. The alloy system examined in this work is AgAu since the corrosion processes in this system occur by a dealloying mechanism that results in the formation of nanoporous gold. The physics behind the dealloying process as well as the resulting formation of nanoporous gold is today well understood.
Two important aspects of the film-induced cleavage mechanism are examined in this work: dynamic fracture in monolithic nanoporous gold and crack injection. In crack injection there is a finite thickness dealloyed layer formed on a AgAu alloy sample and the question of whether or not a crack that nucleates within this layer can travel for some finite distance into the un-corroded parent phase alloy is addressed. Dynamic fracture tests were performed on single edge-notched monolithic nanoporous gold samples as well as “infinite strip” sample configurations for which the stress intensity remains constant over a significant portion of the crack length. High-speed photography was used to measure the crack velocity. In the dynamic fracture experiments cracks were observed to travel at speeds as large as 270 m/s corresponding to about 68% of the Raleigh wave velocity. Crack injection experiments were performed on single crystal Ag77Au23, polycrystalline Ag72Au28 and pure gold, all of which had thin nanoporous gold layers on the surface of samples. Through-thickness fracture was seen in both the single crystal and polycrystalline samples and there was an indication of ~ 1 μm injected cracks into pure gold. These results have important implications for the operation of the film-induced cleavage mechanism and represent a first step in the development of a fundamental model of SCC.
ContributorsChen, Xiying (Author) / Sieradzki, Karl (Thesis advisor) / Jiao, Yang (Committee member) / Oswald, Jay (Committee member) / Crozier, Peter (Committee member) / Peralta, Pedro (Committee member) / Arizona State University (Publisher)
Created2016
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