Theses and Dissertations
This collection includes both ASU Theses and Dissertations, submitted by graduate students, and the Barrett, Honors College theses submitted by undergraduate students.
Displaying 1 - 6 of 6
Filtering by
- Creators: Oswald, Jay
Description
In this dissertation, the results of our comprehensive computational studies of disordered jammed (i.e., mechanically stable) packings of hard particles are presented, including the family of superdisks in 2D and ellipsoids in 3D Euclidean space. Following a very brief introduction to the hard-particle systems, the event driven molecular dynamics (EDMD) employed to generate the packing ensembles will be discussed. A large number of 2D packing configurations of superdisks are subsequently analyzed, through which a relatively accurate theoretical scheme for packing-fraction prediction based on local particle contact configurations is proposed and validated via additional numerical simulations. Moreover, the studies on binary ellipsoid packing in 3D are briefly discussed and the effects of different geometrical parameters on the final packing fraction are analyzed.
ContributorsXu, Yaopengxiao (Author) / Jiao, Yang (Thesis advisor) / Oswald, Jay (Committee member) / Liu, Yongming (Committee member) / Arizona State University (Publisher)
Created2014
Description
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
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 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.
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
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 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
Description
This investigation is focused on the consideration of structural uncertainties in nearly-straight pipes conveying fluid and on the effects of these uncertainties on the dynamic response and stability of those pipes. Of interest more specifically are the structural uncertainties which affect directly the fluid flow and its feedback on the structural response, e.g., uncertainties on/variations of the inner cross-section and curvature of the pipe. Owing to the complexity of introducing such uncertainties directly in finite element models, it is desired to proceed directly at the level of modal models by randomizing simultaneously the appropriate mass, stiffness, and damping matrices. The maximum entropy framework is adopted to carry out the stochastic modeling of these matrices with appropriate symmetry constraints guaranteeing that the nature, e.g., divergence or flutter, of the bifurcation is preserved when introducing uncertainty.
To support the formulation of this stochastic ROM, a series of finite element computations are first carried out for pipes with straight centerline but inner radius varying randomly along the pipe. The results of this numerical discovery effort demonstrate that the dominant effects originate from the variations of the exit flow speed, induced by the change in inner cross-section at the pipe end, with the uncertainty on the cross-section at other locations playing a secondary role. Relying on these observations, the stochastic reduced order model is constructed to model separately the uncertainty in inner cross-section at the pipe end and at other locations. Then, the fluid related mass, damping, and stiffness matrices of this stochastic reduced order model (ROM) are all determined from a single random matrix and a random variable. The predictions from this stochastic ROM are found to closely match the corresponding results obtained with the randomized finite element model. It is finally demonstrated that this stochastic ROM can easily be extended to account for the small effects due to uncertainty in pipe curvature.
To support the formulation of this stochastic ROM, a series of finite element computations are first carried out for pipes with straight centerline but inner radius varying randomly along the pipe. The results of this numerical discovery effort demonstrate that the dominant effects originate from the variations of the exit flow speed, induced by the change in inner cross-section at the pipe end, with the uncertainty on the cross-section at other locations playing a secondary role. Relying on these observations, the stochastic reduced order model is constructed to model separately the uncertainty in inner cross-section at the pipe end and at other locations. Then, the fluid related mass, damping, and stiffness matrices of this stochastic reduced order model (ROM) are all determined from a single random matrix and a random variable. The predictions from this stochastic ROM are found to closely match the corresponding results obtained with the randomized finite element model. It is finally demonstrated that this stochastic ROM can easily be extended to account for the small effects due to uncertainty in pipe curvature.
ContributorsShah, Shrinil (Author) / Mignolet, Marc P (Thesis advisor) / Liu, Yongming (Committee member) / Oswald, Jay (Committee member) / Arizona State University (Publisher)
Created2017
Description
A method for modelling the interactions of dislocations with inclusions has been developed to analyse toughening mechanisms in alloys. This method is different from the superposition method in that infinite domain solutions and image stress fields are not superimposed. The method is based on the extended finite element method (XFEM) in which the dislocations are modelled according to the Volterra dislocation model. Interior discontinuities are introduced across dislocation glide planes using enrichment functions and the resulting boundary value problem is solved through the standard finite element variational approach. The level set method is used to describe the geometry of the dislocation glide planes without any explicit treatment of the interface geometry which provides a convenient and an appealing means for describing the dislocation. A method for estimating the Peach-Koehler force by the domain form of J-integral is considered. The convergence and accuracy of the method are studied for an edge dislocation interacting with a free surface where analytical solutions are available. The force converges to the exact solution at an optimal rate for linear finite elements. The applicability of the method to dislocation interactions with inclusions is illustrated with a system of Aluminium matrix containing Aluminium-copper precipitates. The effect of size, shape and orientation of the inclusions on an edge dislocation for a difference in stiffness and coefficient of thermal expansion of the inclusions and matrix is considered. The force on the dislocation due to a hard inclusion increased by 8% in approaching the sharp corners of a square inclusion than a circular inclusion of equal area. The dislocation experienced 24% more force in moving towards the edges of a square shaped inclusion than towards its centre. When the areas of the inclusions were halved, 30% less force was exerted on the dislocation. This method was used to analyse interfaces with mismatch strains. Introducing eigenstrains equal to 0.004 to the elastic mismatch increased the force by 15 times for a circular inclusion. The energy needed to move an edge dislocation through a domain filled with circular inclusions is 4% more than that needed for a domain with square shaped inclusions.
ContributorsVeeresh, Pawan (Author) / Oswald, Jay (Thesis advisor) / Jiang, Hanqing (Committee member) / Liu, Yongming (Committee member) / Arizona State University (Publisher)
Created2016