Matching Items (17)
Filtering by
- All Subjects: Mechanical Engineering
- Genre: Masters Thesis
- Creators: Oswald, Jay
Description
Cavitation erosion is a significant cause of wear in marine components, such as impellers, propellers or rudders. While the erosion process has been widely studied on metals, the effect of cavitation on polymers is not well-understood. The stress response in metals differs greatly from that of polymers, e.g. rate and temperature effects are far more important, thus damage and wear mechanisms of polymers under cavitating flows are significantly different. In this work, heat-driven failure caused by viscous dissipation and void nucleation resulting from tensile stresses arising from stress wave reflections are investigated as two possible material failure mechanisms.
As a first step in developing a fundamental understanding of the cavitation erosion process on polymer surfaces, simulations are performed of the collapse of individual bubbles against a compliant surface e.g. metallic substrates with polyurea coatings. The surface response of collapse-driven impact loads is represented by a idealized, time-dependent, Gaussian pressure distribution on the surface. A two-dimensional distribution of load radii and durations is considered corresponding to characteristic of cavitating flows accelerated erosion experiments. Finite element simulations are performed to fit a response curve that relates the loading parameters to the energy dissipated in the coating and integrated with collapse statistics to generate an expected heat input into the coating.
The impulsive pressure, which is generated due to bubble collapse, impacts the material and generates intense shock waves. The stress waves within the material reflects by interaction with the substrate. A transient region of high tensile stress is produced by the interaction of these waves. Simulations suggests that maximum hydrostatic tension which cause failure of polyurea layer is observed in thick coating. Also, the dissipated viscous energy and corresponding temperature rise in a polyurea is calculated, and it is concluded that temperature has influence on deformation.
As a first step in developing a fundamental understanding of the cavitation erosion process on polymer surfaces, simulations are performed of the collapse of individual bubbles against a compliant surface e.g. metallic substrates with polyurea coatings. The surface response of collapse-driven impact loads is represented by a idealized, time-dependent, Gaussian pressure distribution on the surface. A two-dimensional distribution of load radii and durations is considered corresponding to characteristic of cavitating flows accelerated erosion experiments. Finite element simulations are performed to fit a response curve that relates the loading parameters to the energy dissipated in the coating and integrated with collapse statistics to generate an expected heat input into the coating.
The impulsive pressure, which is generated due to bubble collapse, impacts the material and generates intense shock waves. The stress waves within the material reflects by interaction with the substrate. A transient region of high tensile stress is produced by the interaction of these waves. Simulations suggests that maximum hydrostatic tension which cause failure of polyurea layer is observed in thick coating. Also, the dissipated viscous energy and corresponding temperature rise in a polyurea is calculated, and it is concluded that temperature has influence on deformation.
ContributorsPanwar, Ajay (Author) / Oswald, Jay (Thesis advisor) / Dooley, Kevin (Committee member) / Chen, Kangping (Committee member) / Arizona State University (Publisher)
Created2015
Description
Recent studies of the occurrence of post-flutter limit cycle oscillations (LCO) of the F-16 have provided good support to the long-standing hypothesis that this phenomenon involves a nonlinear structural damping. A potential mechanism for the appearance of nonlinearity in the damping are the nonlinear geometric effects that arise when the deformations become large enough to exceed the linear regime. In this light, the focus of this investigation is first on extending nonlinear reduced order modeling (ROM) methods to include viscoelasticity which is introduced here through a linear Kelvin-Voigt model in the undeformed configuration. Proceeding with a Galerkin approach, the ROM governing equations of motion are obtained and are found to be of a generalized van der Pol-Duffing form with parameters depending on the structure and the chosen basis functions. An identification approach of the nonlinear damping parameters is next proposed which is applicable to structures modeled within commercial finite element software.
The effects of this nonlinear damping mechanism on the post-flutter response is next analyzed on the Goland wing through time-marching of the aeroelastic equations comprising a rational fraction approximation of the linear aerodynamic forces. It is indeed found that the nonlinearity in the damping can stabilize the unstable aerodynamics and lead to finite amplitude limit cycle oscillations even when the stiffness related nonlinear geometric effects are neglected. The incorporation of these latter effects in the model is found to further decrease the amplitude of LCO even though the dominant bending motions do not seem to stiffen as the level of displacements is increased in static analyses.
The effects of this nonlinear damping mechanism on the post-flutter response is next analyzed on the Goland wing through time-marching of the aeroelastic equations comprising a rational fraction approximation of the linear aerodynamic forces. It is indeed found that the nonlinearity in the damping can stabilize the unstable aerodynamics and lead to finite amplitude limit cycle oscillations even when the stiffness related nonlinear geometric effects are neglected. The incorporation of these latter effects in the model is found to further decrease the amplitude of LCO even though the dominant bending motions do not seem to stiffen as the level of displacements is increased in static analyses.
ContributorsSong, Pengchao (Author) / Mignolet, Marc P (Thesis advisor) / Chattopadhyay, Aditi (Committee member) / Oswald, Jay (Committee member) / Arizona State University (Publisher)
Created2015
Description
The football helmet is a device used to help mitigate the occurrence of impact-related traumatic (TBI) and minor traumatic brain injuries (mTBI) in the game of American football. The current design methodology of using a hard shell with an energy absorbing liner may be adequate for minimizing TBI, however it has had less effect in minimizing mTBI. The latest research in brain injury mechanisms has established that the current design methodology has produced a helmet to reduce linear acceleration of the head. However, angular accelerations also have an adverse effect on the brain response, and must be investigated as a contributor of brain injury.
To help better understand how the football helmet design features effect the brain response during impact, this research develops a validated football helmet model and couples it with a full LS-DYNA human body model developed by the Global Human Body Modeling Consortium (v4.1.1). The human body model is a conglomeration of several validated models of different sections of the body. Of particular interest for this research is the Wayne State University Head Injury Model for modeling the brain. These human body models were validated using a combination of cadaveric and animal studies. In this study, the football helmet was validated by laboratory testing using drop tests on the crown of the helmet. By coupling the two models into one finite element model, the brain response to impact loads caused by helmet design features can be investigated. In the present research, LS-DYNA is used to study a helmet crown impact with a rigid steel plate so as to obtain the strain-rate, strain, and stress experienced in the corpus callosum, midbrain, and brain stem as these anatomical regions are areas of concern with respect to mTBI.
To help better understand how the football helmet design features effect the brain response during impact, this research develops a validated football helmet model and couples it with a full LS-DYNA human body model developed by the Global Human Body Modeling Consortium (v4.1.1). The human body model is a conglomeration of several validated models of different sections of the body. Of particular interest for this research is the Wayne State University Head Injury Model for modeling the brain. These human body models were validated using a combination of cadaveric and animal studies. In this study, the football helmet was validated by laboratory testing using drop tests on the crown of the helmet. By coupling the two models into one finite element model, the brain response to impact loads caused by helmet design features can be investigated. In the present research, LS-DYNA is used to study a helmet crown impact with a rigid steel plate so as to obtain the strain-rate, strain, and stress experienced in the corpus callosum, midbrain, and brain stem as these anatomical regions are areas of concern with respect to mTBI.
ContributorsDarling, Timothy (Author) / Rajan, Subramaniam D. (Thesis advisor) / Muthuswamy, Jitendran (Thesis advisor) / Oswald, Jay (Committee member) / Mignolet, Marc (Committee member) / Arizona State University (Publisher)
Created2014
Description
Monte Carlo simulations are traditionally carried out for the determination of the amplification of forced vibration response of turbomachine/jet engine blades to mistuning. However, this effort can be computationally time consuming even when using the various reduced order modeling techniques. Accordingly, some investigations in the past have focused on obtaining simple approximate estimates for this amplification. In particular, two of these have proposed the use of harmonic patterns of the blade properties around the disk as an approximate alternative to the many random patterns of Monte Carlo analyses. These investigations, while quite encouraging, have relied solely on single degree of freedom per sector models of the rotor.
In this light, the overall focus of the present effort is a revisit of harmonic
mistuning of rotors focusing first the confirmation of the previously obtained findings with a more detailed model of the blisk in both conditions of an isolated blade-dominated resonance and of a veering between blade and disk dominated modes. The latter condition cannot be simulated by a single degree of freedom per sector model. Further, the analysis will consider the distinct cases of mistuning due to variations of material properties (Young's modulus) and geometric properties (geometric mistuning). In the single degree of freedom model, both mistuning types are equivalent but they are not, as demonstrated here, in more realistic models. The difference arises because changes in geometry induce not only changes in natural frequencies of the blades alone but of their modes and the importance of these two sources of variability is discussed with both Monte Carlo simulation and harmonic mistuning results.
The present investigation focuses also on the possible extension of the harmonic mistuning concept and of its quantitative information that can be derived from such analyses. From it, a novel measure of blade-disk coupling is introduced and assessed in comparison with the coupling index introduced in the past. In conclusions, the low cost of harmonic mistuning computations in comparison with full Monte Carlo simulations is
demonstrated to be worthwhile to elucidate the basic behavior of the mistuned rotor in a random setting.
In this light, the overall focus of the present effort is a revisit of harmonic
mistuning of rotors focusing first the confirmation of the previously obtained findings with a more detailed model of the blisk in both conditions of an isolated blade-dominated resonance and of a veering between blade and disk dominated modes. The latter condition cannot be simulated by a single degree of freedom per sector model. Further, the analysis will consider the distinct cases of mistuning due to variations of material properties (Young's modulus) and geometric properties (geometric mistuning). In the single degree of freedom model, both mistuning types are equivalent but they are not, as demonstrated here, in more realistic models. The difference arises because changes in geometry induce not only changes in natural frequencies of the blades alone but of their modes and the importance of these two sources of variability is discussed with both Monte Carlo simulation and harmonic mistuning results.
The present investigation focuses also on the possible extension of the harmonic mistuning concept and of its quantitative information that can be derived from such analyses. From it, a novel measure of blade-disk coupling is introduced and assessed in comparison with the coupling index introduced in the past. In conclusions, the low cost of harmonic mistuning computations in comparison with full Monte Carlo simulations is
demonstrated to be worthwhile to elucidate the basic behavior of the mistuned rotor in a random setting.
ContributorsSahoo, Saurav (Author) / Mignolet, Marc Paul (Thesis advisor) / Chattopadhyay, Aditi (Committee member) / Oswald, Jay (Committee member) / Arizona State University (Publisher)
Created2014
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
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
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
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