Matching Items (41)
Filtering by
- All Subjects: Mechanical Engineering
- Creators: Oswald, Jay
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
In this thesis, a FORTRAN code is rewritten in C++ with an object oriented ap-
proach. There are several reasons for this purpose. The first reason is to establish
the basis of a GPU programming. To write programs that utilize GPU hardware,
CUDA or OpenCL is used which only support C and C++. FORTRAN has a feature
that lets its programs to call C/C++ functions. FORTRAN sends relevant data to
C/C++, which in turn sends that data to OpenCL. Although this approach works,
it makes the code messy and bulky and in the end more difficult to deal with. More-
over, there is a slight performance decrease from the additional data copy. This is
the motivation to have the code entirely written in C++ to make it more uniform,
efficient and clean. The second reason is the object oriented feature of the C++. The
“abstraction”, “inheritance” and “run-time polymorphism” features of C++ provide
some form of classes and objects, the ability to build new abstractions, and some
form of run-time binding, respectively. In recent years, some of popular codes has
been rewritten in C++ which were initially in FORTRAN. One of these softwares is
LAMMPS.
In this code the level set equation is solved by RLSG method to track the interface in
two phase flow. In gas/fluid flows, the surface tension is important and only exists at
the interface. Therefore, the location and some geometric features of interface need
to be evaluated which can be achieved by solving the level set equation.
proach. There are several reasons for this purpose. The first reason is to establish
the basis of a GPU programming. To write programs that utilize GPU hardware,
CUDA or OpenCL is used which only support C and C++. FORTRAN has a feature
that lets its programs to call C/C++ functions. FORTRAN sends relevant data to
C/C++, which in turn sends that data to OpenCL. Although this approach works,
it makes the code messy and bulky and in the end more difficult to deal with. More-
over, there is a slight performance decrease from the additional data copy. This is
the motivation to have the code entirely written in C++ to make it more uniform,
efficient and clean. The second reason is the object oriented feature of the C++. The
“abstraction”, “inheritance” and “run-time polymorphism” features of C++ provide
some form of classes and objects, the ability to build new abstractions, and some
form of run-time binding, respectively. In recent years, some of popular codes has
been rewritten in C++ which were initially in FORTRAN. One of these softwares is
LAMMPS.
In this code the level set equation is solved by RLSG method to track the interface in
two phase flow. In gas/fluid flows, the surface tension is important and only exists at
the interface. Therefore, the location and some geometric features of interface need
to be evaluated which can be achieved by solving the level set equation.
ContributorsSafarkhani, Salar (Author) / Herrmann, Mrcus (Thesis advisor) / Oswald, Jay (Committee member) / Rykczewski, Konrad (Committee member) / Arizona State University (Publisher)
Created2015
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 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.
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
Description
A new critical plane-energy model is proposed in this thesis for multiaxial fatigue life prediction of homogeneous and heterogeneous materials. Brief review of existing methods, especially on the critical plane-based and energy-based methods, are given first. Special focus is on one critical plane approach which has been shown to work for both brittle and ductile metals. The key idea is to automatically change the critical plane orientation with respect to different materials and stress states. One potential drawback of the developed model is that it needs an empirical calibration parameter for non-proportional multiaxial loadings since only the strain terms are used and the out-of-phase hardening cannot be considered. The energy-based model using the critical plane concept is proposed with help of the Mroz-Garud hardening rule to explicitly include the effect of non-proportional hardening under fatigue cyclic loadings. Thus, the empirical calibration for non-proportional loading is not needed since the out-of-phase hardening is naturally included in the stress calculation. The model predictions are compared with experimental data from open literature and it is shown the proposed model can work for both proportional and non-proportional loadings without the empirical calibration. Next, the model is extended for the fatigue analysis of heterogeneous materials integrating with finite element method. Fatigue crack initiation of representative volume of heterogeneous materials is analyzed using the developed critical plane-energy model and special focus is on the microstructure effect on the multiaxial fatigue life predictions. Several conclusions and future work is drawn based on the proposed study.
ContributorsWei, Haoyang (Author) / Liu, Yongming (Thesis advisor) / Jiang, Hanqing (Committee member) / Oswald, Jay (Committee member) / Arizona State University (Publisher)
Created2016
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
In this paper, at first, analytical formulation of J-integral for a non-local particle model (VCPM) using atomic scale finite element method is proposed for fracture analysis of 2D solids. A brief review of classical continuum-based J-integral and anon-local lattice particle method is given first. Following this, detailed derivation for the J-integral in discrete particle system is given using the energy equivalence and stress-tensor mapping between the continuum mechanics and lattice-particle system.With the help of atomistic finite element method, the J-integral is expressed as a summation of the corresponding terms in the particle system.
Secondly, a coupling algorithm between a non-local particle method (VCPM) and the classical finite element method (FEM) is discussed to gain the advantages of both methods for fracture analysis in large structures. In this algorithm, the discrete VCPM particle and the continuum FEM domains are solved within a unified theoretical framework. A transitional element technology is developed to smoothly link the 10-particles element with the traditional FEM elements to guaranty the continuity and consistency at the coupling interface. An explicit algorithm for static simulation is developed.
Finally, numerical examples are illustrated for the accuracy, convergence, and path-independence of the derived J-integral formulation. Discussions on the comparison with alternative estimation methods and potential application for fracture simulation are given. The accuracy and efficiency of the coupling algorithm are tested by several benchmark problems such as static crack simulation.
Secondly, a coupling algorithm between a non-local particle method (VCPM) and the classical finite element method (FEM) is discussed to gain the advantages of both methods for fracture analysis in large structures. In this algorithm, the discrete VCPM particle and the continuum FEM domains are solved within a unified theoretical framework. A transitional element technology is developed to smoothly link the 10-particles element with the traditional FEM elements to guaranty the continuity and consistency at the coupling interface. An explicit algorithm for static simulation is developed.
Finally, numerical examples are illustrated for the accuracy, convergence, and path-independence of the derived J-integral formulation. Discussions on the comparison with alternative estimation methods and potential application for fracture simulation are given. The accuracy and efficiency of the coupling algorithm are tested by several benchmark problems such as static crack simulation.
ContributorsZope, Jayesh (Author) / Liu, Yongming (Thesis advisor) / Oswald, Jay (Committee member) / Jiang, Hanqing (Committee member) / Arizona State University (Publisher)
Created2016
Description
Characterization and modeling of deformation and failure in metallic materials under extreme conditions, such as the high loads and strain rates found under shock loading due to explosive detonation and high velocity-impacts, are extremely important for a wide variety of military and industrial applications. When a shock wave causes stress in a material that exceeds the elastic limit, plasticity and eventually spallation occur in the material. The process of spall fracture, which in ductile materials stems from strain localization, void nucleation, growth and coalescence, can be caused by microstructural heterogeneity. The analysis of void nucleation performed from a microstructurally explicit simulation of a spall damage evolution in a multicrystalline copper indicated triple junctions as the preferred sites for incipient damage nucleation revealing 75% of them with at least two grain boundaries with misorientation angle between 20-55°. The analysis suggested the nature of the boundaries connecting at a triple junction is an indicator of their tendency to localize spall damage. The results also showed that damage propagated preferentially into one of the high angle boundaries after voids nucleate at triple junctions. Recently the Rayleigh-Taylor Instability (RTI) and the Richtmyer-Meshkov Instability (RMI) have been used to deduce dynamic material strength at very high pressures and strain rates. The RMI is used in this work since it allows using precise diagnostics such as Transient Imaging Displacement Interferometry (TIDI) due to its slower linear growth rate. The Preston-Tonks-Wallace (PTW) model is used to study the effects of dynamic strength on the behavior of samples with a fed-thru RMI, induced via direct laser drive on a perturbed surface, on stability of the shock front and the dynamic evolution of the amplitudes and velocities of the perturbation imprinted on the back (flat) surface by the perturbed shock front. Simulation results clearly showed that the amplitude of the hydrodynamic instability increases with a decrease in strength and vice versa and that the amplitude of the perturbed shock front produced by the fed-thru RMI is also affected by strength in the same way, which provides an alternative to amplitude measurements to study strength effects under dynamic conditions. Simulation results also indicate the presence of second harmonics in the surface perturbation after a certain time, which were also affected by the material strength.
ContributorsGautam, Sudrishti (Author) / Peralta, Pedro (Thesis advisor) / Oswald, Jay (Committee member) / Solanki, Kiran (Committee member) / Arizona State University (Publisher)
Created2016
Description
ABSTRACT
A large fraction of the total energy consumption in the world comes from heating and cooling of buildings. Improving the energy efficiency of buildings to reduce the needs of seasonal heating and cooling is one of the major challenges in sustainable development. In general, the energy efficiency depends on the geometry and material of the buildings. To explore a framework for accurately assessing this dependence, detailed 3-D thermofluid simulations are performed by systematically sweeping the parameter space spanned by four parameters: the size of building, thickness and material of wall, and fractional size of window. The simulations incorporate realistic boundary conditions of diurnally-varying temperatures from observation, and the effect of fluid flow with explicit thermal convection inside the building. The outcome of the numerical simulations is synthesized into a simple map of an index of energy efficiency in the parameter space which can be used by stakeholders to quick look-up the energy efficiency of a proposed design of a building before its construction. Although this study only considers a special prototype of buildings, the framework developed in this work can potentially be used for a wide range of buildings and applications.
A large fraction of the total energy consumption in the world comes from heating and cooling of buildings. Improving the energy efficiency of buildings to reduce the needs of seasonal heating and cooling is one of the major challenges in sustainable development. In general, the energy efficiency depends on the geometry and material of the buildings. To explore a framework for accurately assessing this dependence, detailed 3-D thermofluid simulations are performed by systematically sweeping the parameter space spanned by four parameters: the size of building, thickness and material of wall, and fractional size of window. The simulations incorporate realistic boundary conditions of diurnally-varying temperatures from observation, and the effect of fluid flow with explicit thermal convection inside the building. The outcome of the numerical simulations is synthesized into a simple map of an index of energy efficiency in the parameter space which can be used by stakeholders to quick look-up the energy efficiency of a proposed design of a building before its construction. Although this study only considers a special prototype of buildings, the framework developed in this work can potentially be used for a wide range of buildings and applications.
ContributorsJain, Gaurav (Author) / Huang, Huei-Ping (Thesis advisor) / Ren, Yi (Committee member) / Oswald, Jay (Committee member) / Arizona State University (Publisher)
Created2016
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
A previously developed small time scale fatigue crack growth model is improved, modified and extended with an emphasis on creating the simplest models that maintain the desired level of accuracy for a variety of materials. The model provides a means of estimating load sequence effects by continuously updating the crack opening stress every cycle, in a simplified manner. One of the significant phenomena of the crack opening stress under negative stress ratio is the residual tensile stress induced by the applied compressive stress. A modified coefficient is introduced to determine the extent to which residual stress impact the crack closure and is observed to vary for different materials. Several other literature models for crack closure under constant loading are also reviewed and compared with the proposed model. The modified model is then shown to predict several sets of published test results under constant loading for a variety of materials.
The crack opening stress is formalized as a function of the plastic zone sizes at the crack tip and the current crack length, which provided a means of approximation, accounting for both acceleration and retardation effects in a simplified manner. A sensitivity parameter is introduced to modify the enlarged plastic zone due to overload, to better fit the delay cycles with the test data and is observed to vary for different materials. Furthermore, the interaction effect induced by the combination of overload and underload sequence is modeled by depleting the compressive plastic zone due to an overload with the tensile plastic zone due to an underload. A qualitative analysis showed the simulation capacity of the small time scale model under different load types. A good agreement between prediction and test data for several irregular load types proved the applicability of the small time scale model under variable amplitude loading.
The crack opening stress is formalized as a function of the plastic zone sizes at the crack tip and the current crack length, which provided a means of approximation, accounting for both acceleration and retardation effects in a simplified manner. A sensitivity parameter is introduced to modify the enlarged plastic zone due to overload, to better fit the delay cycles with the test data and is observed to vary for different materials. Furthermore, the interaction effect induced by the combination of overload and underload sequence is modeled by depleting the compressive plastic zone due to an overload with the tensile plastic zone due to an underload. A qualitative analysis showed the simulation capacity of the small time scale model under different load types. A good agreement between prediction and test data for several irregular load types proved the applicability of the small time scale model under variable amplitude loading.
ContributorsVenkatesan, Karthik Rajan (Author) / Liu, Yongming (Thesis advisor) / Oswald, Jay (Committee member) / Jiang, Hanqing (Committee member) / Arizona State University (Publisher)
Created2016