Matching Items (8)
Filtering by
- All Subjects: Mechanical Engineering
- Genre: Doctoral Dissertation
Description
The focus of this investigation is on the optimum placement of a limited number of dampers, fewer than the number of blades, on a bladed disk to induce the smallest amplitude of blade response. The optimization process considers the presence of random mistuning, i.e. small involuntary variations in blade stiffness properties resulting, say, from manufacturing variability. Designed variations of these properties, known as intentional mistuning, is considered as an option to reduce blade response and the pattern of two blade types (A and B blades) is then part of the optimization in addition to the location of dampers on the disk. First, this study focuses on the formulation and validation of dedicated algorithms for the selection of the damper locations and the intentional mistuning pattern. Failure of one or several of the dampers could lead to a sharp rise in blade response and this issue is addressed by including, in the optimization, the possibility of damper failure to yield a fail-safe solution. The high efficiency and accuracy of the optimization algorithms is assessed in comparison with computationally very demanding exhaustive search results. Second, the developed optimization algorithms are applied to nonlinear dampers (underplatform friction dampers), as well as to blade-blade dampers, both linear and nonlinear. Further, the optimization of blade-only and blade-blade linear dampers is extended to include uncertainty or variability in the damper properties induced by manufacturing or wear. It is found that the optimum achieved without considering such uncertainty is robust with respect to it. Finally, the potential benefits of using two different types of friction dampers differing in their masses (A and B types), on a bladed disk are considered. Both A/B pattern and the damper masses are optimized to obtain the largest benefit compared to using identical dampers of optimized masses on every blade. Four situations are considered: tuned disks, disks with random mistuning of blade stiffness, and, disks with random mistuning of both blade stiffness and damper normal forces with and without damper variability induced by manufacturing and wear. In all cases, the benefit of intentional mistuning of friction dampers is small, of the order of a few percent.
ContributorsMurthy, Raghavendra Narasimha (Author) / Mignolet, Marc P (Thesis advisor) / Rajan, Subramaniam D. (Committee member) / Lentz, Jeff (Committee member) / Chattopadhyay, Aditi (Committee member) / Jiang, Hanqing (Committee member) / Arizona State University (Publisher)
Created2012
Description
With the increasing focus on developing environmentally benign electronic packages, lead-free solder alloys have received a great deal of attention. Mishandling of packages, during manufacture, assembly, or by the user may cause failure of solder joint. A fundamental understanding of the behavior of lead-free solders under mechanical shock conditions is lacking. Reliable experimental and numerical analysis of lead-free solder joints in the intermediate strain rate regime need to be investigated. This dissertation mainly focuses on exploring the mechanical shock behavior of lead-free tin-rich solder alloys via multiscale modeling and numerical simulations. First, the macroscopic stress/strain behaviors of three bulk lead-free tin-rich solders were tested over a range of strain rates from 0.001/s to 30/s. Finite element analysis was conducted to determine appropriate specimen geometry that could reach a homogeneous stress/strain field and a relatively high strain rate. A novel self-consistent true stress correction method is developed to compensate the inaccuracy caused by the triaxial stress state at the post-necking stage. Then the material property of micron-scale intermetallic was examined by micro-compression test. The accuracy of this measure is systematically validated by finite element analysis, and empirical adjustments are provided. Moreover, the interfacial property of the solder/intermetallic interface is investigated, and a continuum traction-separation law of this interface is developed from an atomistic-based cohesive element method. The macroscopic stress/strain relation and microstructural properties are combined together to form a multiscale material behavior via a stochastic approach for both solder and intermetallic. As a result, solder is modeled by porous plasticity with random voids, and intermetallic is characterized as brittle material with random vulnerable region. Thereafter, the porous plasticity fracture of the solders and the brittle fracture of the intermetallics are coupled together in one finite element model. Finally, this study yields a multiscale model to understand and predict the mechanical shock behavior of lead-free tin-rich solder joints. Different fracture patterns are observed for various strain rates and/or intermetallic thicknesses. The predictions have a good agreement with the theory and experiments.
ContributorsFei, Huiyang (Author) / Jiang, Hanqing (Thesis advisor) / Chawla, Nikhilesh (Thesis advisor) / Tasooji, Amaneh (Committee member) / Mobasher, Barzin (Committee member) / Rajan, Subramaniam D. (Committee member) / Arizona State University (Publisher)
Created2011
Description
Gels are three-dimensional polymer networks with entrapped solvent (water etc.). They bear amazing features such as stimuli-responsive (temperature, PH, electric field etc.), high water content and biocompatibility and thus find a lot of applications. To understand the complex physics behind gel's swelling phenomenon, it is important to build up fundamental mechanical model and extend to complicated cases. In this dissertation, a coupled large deformation and diffusion model regarding gel's swelling behavior is presented. In this model, free-energy of the total gel is constituted by polymer stretching energy and polymer-solvent mixing energy. In-house nonlinear finite element code is implemented with fast computational capability. Complex phenomenon such as buckling and healing of cracked gel by swelling are studied. Due to the wide coverage of polymeric materials and solvents, solvent diffusion in gels not only follows Fickian diffusion law where concentration map is continuous but also follows non-Fickian diffusion law where concentration map shows high gradient. Phenomenological model with viscoelastic polymer constitutive and concentration dependent diffusivity is created. The model well captures this special diffusion phenomenon such as sharp diffusion front and distinctive swollen and unswollen region.
ContributorsZhang, Jiaping (Author) / Jiang, Hanqing (Thesis advisor) / Peralta, Pedro (Committee member) / Dai, Lenore (Committee member) / Rajan, Subramaniam D. (Committee member) / Chawla, Nikhilesh (Committee member) / Arizona State University (Publisher)
Created2012
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 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
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
Hydrodynamic phenomena such as the Rayleigh-Taylor (RT) and Richtmyer-Meshkov (RM) instabilities can be described by exponential/linear growth of surface perturbations at a bimaterial interface when subjected to constant/impulsive acceleration. A challenge in designing systems to mitigate or exploit these effects is the lack of accurate material models at large dynamic strain rates and pressures. In particular, little stress-strain constitutive information at large strain rates and pressures is available for transient material phases formed at high pressures, and the continuum effect the phase transformation process has on the instability evolution. In this work, a phase-aware isotropic strength model is developed and partially validated with a novel RM-based instability experiment in addition to existing data from the literature. With the validated material model additional simulations are performed to provide insight into to the role that robust material constitutive behavior (e.g., pressure, temperature, rate dependence) has on RM instability and how RM instability experiments can be used to characterize and validated expected material behavior. For phase aware materials, particularly iron in this work, the simulations predict a strong dependence on the Atwood number that single phase materials do not have. At Atwood numbers close to unity, and pressures in the high pressure stability region, the high pressure phase dominates the RM evolution. However, at Atwood numbers close to negative one, the RM evolution is only weakly affected by the high-pressure phase even for shocks well above the phase transformation threshold. In addition to RM evolution this work looks at the closely related shock front perturbation evolution. Existing analytical models for isentropic processes in gases and liquids are modified for metal equation of states and plastic behavior for the first time. It is found that the presence of a volume collapsing phase transformation with increased pressure causes shock front perturbations to decay sooner, while plastic strength has the opposite effect which is significantly different from the effect viscosity has. These results suggest additional experimental setups to validate material models, or relevant material parameters that can be optimized for system design objectives, e.g., minimize feed through perturbations in inertial confinement fusion capsules.
ContributorsOpie, Saul (Author) / Peralta, Pedro (Thesis advisor) / Loomis, Eric (Committee member) / Oswald, Jay (Committee member) / Rajan, Subramaniam D. (Committee member) / Solanki, Kiran (Committee member) / Arizona State University (Publisher)
Created2017
Description
The focus of this dissertation is first on understanding the difficulties involved in constructing reduced order models of structures that exhibit a strong nonlinearity/strongly nonlinear events such as snap-through, buckling (local or global), mode switching, symmetry breaking. Next, based on this understanding, it is desired to modify/extend the current Nonlinear Reduced Order Modeling (NLROM) methodology, basis selection and/or identification methodology, to obtain reliable reduced order models of these structures. Focusing on these goals, the work carried out addressed more specifically the following issues:
i) optimization of the basis to capture at best the response in the smallest number of modes,
ii) improved identification of the reduced order model stiffness coefficients,
iii) detection of strongly nonlinear events using NLROM.
For the first issue, an approach was proposed to rotate a limited number of linear modes to become more dominant in the response of the structure. This step was achieved through a proper orthogonal decomposition of the projection on these linear modes of a series of representative nonlinear displacements. This rotation does not expand the modal space but renders that part of the basis more efficient, the identification of stiffness coefficients more reliable, and the selection of dual modes more compact. In fact, a separate approach was also proposed for an independent optimization of the duals. Regarding the second issue, two tuning approaches of the stiffness coefficients were proposed to improve the identification of a limited set of critical coefficients based on independent response data of the structure. Both approaches led to a significant improvement of the static prediction for the clamped-clamped curved beam model. Extensive validations of the NLROMs based on the above novel approaches was carried out by comparisons with full finite element response data. The third issue, the detection of nonlinear events, was finally addressed by building connections between the eigenvalues of the finite element software (Nastran here) and NLROM tangent stiffness matrices and the occurrence of the ‘events’ which is further extended to the assessment of the accuracy with which the NLROM captures the full finite element behavior after the event has occurred.
i) optimization of the basis to capture at best the response in the smallest number of modes,
ii) improved identification of the reduced order model stiffness coefficients,
iii) detection of strongly nonlinear events using NLROM.
For the first issue, an approach was proposed to rotate a limited number of linear modes to become more dominant in the response of the structure. This step was achieved through a proper orthogonal decomposition of the projection on these linear modes of a series of representative nonlinear displacements. This rotation does not expand the modal space but renders that part of the basis more efficient, the identification of stiffness coefficients more reliable, and the selection of dual modes more compact. In fact, a separate approach was also proposed for an independent optimization of the duals. Regarding the second issue, two tuning approaches of the stiffness coefficients were proposed to improve the identification of a limited set of critical coefficients based on independent response data of the structure. Both approaches led to a significant improvement of the static prediction for the clamped-clamped curved beam model. Extensive validations of the NLROMs based on the above novel approaches was carried out by comparisons with full finite element response data. The third issue, the detection of nonlinear events, was finally addressed by building connections between the eigenvalues of the finite element software (Nastran here) and NLROM tangent stiffness matrices and the occurrence of the ‘events’ which is further extended to the assessment of the accuracy with which the NLROM captures the full finite element behavior after the event has occurred.
ContributorsLin, Jinshan (Author) / Mignolet, Marc (Thesis advisor) / Jiang, Hanqing (Committee member) / Oswald, Jay (Committee member) / Spottswood, Stephen (Committee member) / Rajan, Subramaniam D. (Committee member) / Arizona State University (Publisher)
Created2020