Matching Items (22)
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
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
154549-Thumbnail Image.png
Description
Polymer matrix composites (PMCs) are attractive structural materials due to their high stiffness to low weight ratio. However, unidirectional PMCs have low shear strength and failure can occur along kink bands that develop on compression due to plastic microbuckling that carry strains large enough to induce nonlinear matrix deformation. Reviewing

Polymer matrix composites (PMCs) are attractive structural materials due to their high stiffness to low weight ratio. However, unidirectional PMCs have low shear strength and failure can occur along kink bands that develop on compression due to plastic microbuckling that carry strains large enough to induce nonlinear matrix deformation. Reviewing the literature, a large fraction of the existing work is for uniaxial compression, and the effects of stress gradients, such as those present during bending, have not been as well explored, and these effects are bound to make difference in terms of kink band nucleation and growth. Furthermore, reports on experimental measurements of strain fields leading to and developing inside these bands in the presence of stress gradients are also scarce and need to be addressed to gain a full understanding of their behavior when UDCs are used under bending and other spatially complex stress states.

In a light to bridge the aforementioned gaps, the primary focus of this work is to understand mechanisms for kink band evolution under an influence of stress-gradients induced during bending. Digital image correlation (DIC) is used to measure strains inside and around the kink bands during 3-point bending of samples with 0°/90° stacking made of Ultra-High Molecular Weight Polyethylene Fibers. Measurements indicate bands nucleate at the compression side and propagate into the sample carrying a mixture of large shear and normal strains (~33%), while also decreasing its bending stiffness. Failure was produced by a combination of plastic microbuckling and axial splitting. The microstructure of the kink bands was studied and used in a microstructurally explicit finite element model (FEM) to analyze stresses and strains at ply level in the samples during kink band evolution, using cohesive zone elements to represent the interfaces between plies. Cohesive element properties were deduced by a combination of delamination, fracture and three-point bending tests used to calibrate the FEMs. Modeling results show that the band morphology is sensitive to the shear and opening properties of the interfaces between the plies.
ContributorsPatel, Jay K (Author) / Peralta, Pedro D (Thesis advisor) / Oswald, Jay (Committee member) / Jiang, Hanqing (Committee member) / Solanki, Kiran (Committee member) / Ayyar, Adarsh (Committee member) / Arizona State University (Publisher)
Created2016
154828-Thumbnail Image.png
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

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
155121-Thumbnail Image.png
Description
Finite element simulations modeling the hydrodynamic impact loads subjected to an elastomeric coating were performed to develop an understanding of the performance and failure mechanisms of protective coatings for cavitating environments.

In this work, two major accomplishments were achieved: 1) scaling laws were developed from hydrodynamic principles and numerical

Finite element simulations modeling the hydrodynamic impact loads subjected to an elastomeric coating were performed to develop an understanding of the performance and failure mechanisms of protective coatings for cavitating environments.

In this work, two major accomplishments were achieved: 1) scaling laws were developed from hydrodynamic principles and numerical simulations to allow conversion of measured distributions of pressure peaks in a cavitating flow to distributions of microscopic impact loadings modeling individual bubble collapse events, and 2) a finite strain, thermo-mechanical material model for polyurea-based elastomers was developed using a logarithmic rate formulation and implemented into an explicit finite element code.

Combining the distribution of microscopic impact loads and finite element modeling, a semi-quantitative predictive framework is created to calculate the energy dissipation within the coating which can further the understanding of temperature induced coating failures.

The influence of coating thickness and elastomer rheology on the dissipation of impact energies experienced in cavitating flows has also been explored.

The logarithmic formulation has many desired features for the polyurea constitutive model, such as objectivity, integrability, and additive decomposition compatibility.

A review and discussion on the kinematics in large deformation, including a comparison between Lagrangian and Eulerian descriptions, are presented to explain the issues in building rate-dependent constitutive models in finite strains.

When comparing the logarithmic rate with other conventional rates in test examples, the logarithmic rate shows a better conservation of objectivity and integrability.

The modeling framework was validated by comparing predictions against temperatures measured within coatings subjected to a cavitating jet.

Both the experiments and models show that the temperatures generated, even under mild flow conditions, raise the coating temperature by a significant amount, suggesting that the failure of these coatings under more aggressive flows is thermally induced.

The models show that thin polyurea coatings synthesized with shorter molecular weight soft segments dissipate significantly less energy per impact and conduct heat more efficiently.

This work represents an important step toward understanding thermally induced failure in elastomers subjected to cavitating flows, which provides a foundation for design and optimization of coatings with enhanced erosion resistance.
ContributorsLiao, Xiao (Author) / Oswald, Jay (Thesis advisor) / Liu, Yongming (Committee member) / Jiang, Hanqing (Committee member) / Rajan, Subramaniam D. (Committee member) / Platte, Rodrigo (Committee member) / Arizona State University (Publisher)
Created2016
155251-Thumbnail Image.png
Description
This investigation develops small-size reduced order models (ROMs) that provide an accurate prediction of the response of only part of a structure, referred to as component-centric ROMs. Four strategies to construct such ROMs are presented, the first two of which are based on the Craig-Bampton Method and start with a

This investigation develops small-size reduced order models (ROMs) that provide an accurate prediction of the response of only part of a structure, referred to as component-centric ROMs. Four strategies to construct such ROMs are presented, the first two of which are based on the Craig-Bampton Method and start with a set of modes for the component of interest (the β component). The response in the rest of the structure (the α component) induced by these modes is then determined and optimally represented by applying a Proper Orthogonal Decomposition strategy using Singular Value Decomposition. These first two methods are effectively basis reductions techniques of the CB basis. An approach based on the “Global - Local” Method generates the “global” modes by “averaging” the mass property over α and β comp., respectively (to extract a “coarse” model of α and β) and the “local” modes orthogonal to the “global” modes to add back necessary “information” for β. The last approach adopts as basis for the entire structure its linear modes which are dominant in the β component response. Then, the contributions of other modes in this part of the structure are approximated in terms of those of the dominant modes with close natural frequencies and similar mode shapes in the β component. In this manner, the non-dominant modal contributions are “lumped” onto the dominant ones, to reduce the number of modes for a prescribed accuracy. The four approaches are critically assessed on the structural finite element model of a 9-bay panel with the modal lumping-based method leading to the smallest sized ROMs. Therefore, it is extended to the nonlinear geometric situation and first recast as a rotation of the modal basis to achieve unobservable modes. In the linear case, these modes completely disappear from the formulation owing to orthogonality. In the nonlinear case, however, the generalized coordinates of these modes are still present in the nonlinear terms of the observable modes. A closure-type algorithm is then proposed to eliminate the unobserved generalized coordinates. This approach, its accuracy and computational savings, was demonstrated on a simple beam model and the 9-bay panel model.
ContributorsWang, Yuting (Author) / Mignolet, Marc P (Thesis advisor) / Jiang, Hanqing (Committee member) / Liu, Yongming (Committee member) / Oswald, Jay (Committee member) / Rajan, Subramaniam D. (Committee member) / Spottswood, Stephen M (Committee member) / Arizona State University (Publisher)
Created2017