Search NASA⌕ Search

SEARCH · Search NASA

Results for “numerical modal analysis”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 199 records · Page 11

Perturbation solutions of combustion instability problems

A method involving approximate modal analysis using the Galerkin method followed by an approximate solution of the resulting modal-amplitude equations by the two-variable perturbation method (method of multiple scales) is applied to two problems of pressure-sensitive nonlinear combustion instability in liquid-fuel rocket motors. One problem exhibits self-coupled instability while the other exhibits mode-coupled instability. In both cases it is possible to carry out the entire linear stability analysis and significant portions of the nonlinear stability analysis in closed form. In the problem of self-coupled instability the nonlinear stability boundary and approximate forms of the limit-cycle amplitudes and growth and decay rates are determined in closed form while the exact limit-cycle amplitudes and growth and decay rates are found numerically. In the problem of mode-coupled instability the limit-cycle amplitudes are found in closed form while the growth and decay rates are found numerically. The behavior of the solutions found by the perturbation method are in agreement with solutions obtained using complex numerical methods.

Googerdy, A.↗

Current measurements in the intermittent-contact mode of atomic force microscopy using the Fourier method: a feasibility analysis

Atomic force microscopy (AFM) is an important tool for measuring a variety of nanoscale surface properties, such as topography, viscoelasticity, electrical potential and conductivity. Some of these properties are measured using contact methods (static contact or intermittent contact), while others are measured using noncontact methods. Some properties can be measured using different approaches. Conductivity, in particular, is mapped using the contact-mode method. However, this modality can be destructive to delicate samples, since it involves continuously dragging the cantilever tip on the surface during the raster scan, while a constant tip–sample force is applied. In this paper we discuss a possible approach to develop an intermittent-contact conductive AFM mode based on Fourier analysis, whereby the measured current response consists of higher harmonics of the cantilever oscillation frequency. Such an approach may enable the characterization of soft samples with less damage than contact-mode imaging. To explore its feasibility, we derive the analytical form of the tip–sample current that would be obtained for attractive (noncontact) and repulsive (intermittent-contact) dynamic AFM characterization, and compare it with results obtained from numerical simulations. Although significant instrumentation challenges are anticipated, the modelling results are promising and suggest that Fourier-based higher-harmonics current measurement may enable the development of a reliable intermittent-contact conductive AFM method.

14 SOLAR ENERGY↗

Aeroelastic Analysis of SUGAR Truss-Braced Wing Wind-Tunnel Model Using FUN3D and a Nonlinear Structural Model

Considerable attention has been given in recent years to the design of highly flexible aircraft. The results of numerous studies demonstrate the significant performance benefits of strut-braced wing (SBW) and trussbraced wing (TBW) configurations. Critical aspects of the TBW configuration are its larger aspect ratio, wing span and thinner wings. These aspects increase the importance of considering fluid/structure and control system coupling. This paper presents high-fidelity Navier-Stokes simulations of the dynamic response of the flexible Boeing Subsonic Ultra Green Aircraft Research (SUGAR) truss-braced wing wind-tunnel model. The latest version of the SUGAR TBW finite element model (FEM), v.20, is used in the present simulations. Limit cycle oscillations (LCOs) of the TBW wing/strut/nacelle are simulated at angle-of-attack (AoA) values of -1, 0 and +1 degree. The modal data derived from nonlinear static aeroelastic MSC.Nastran solutions are used at AoAs of -1 and +1 degrees. The LCO amplitude is observed to be dependent on AoA. LCO amplitudes at -1 degree are larger than those at +1 degree. The LCO amplitude at zero degrees is larger than either -1 or +1 degrees. These results correlate well with both wind-tunnel data and the behavior observed in previous studies using linear aerodynamics. The LCO onset at zero degrees AoA has also been computed using unloaded v.20 FEM modes. While the v.20 model increases the dynamic pressure at which LCO onset is observed, it is found that the LCO onset at and above Mach 0.82 is much different than that produced by an earlier version of the FEM, v. 19.

Bartels, Robert E.↗

Exploring Responses of Contact Kelvin Probe Force Microscopy in Triple-Cation Double-Halide Perovskites

Hybrid organic–inorganic perovskite (HOIP) solar cells have undergone unprecedentedly rapid improvement in power conversion efficiencies during the last decade, but controversy remains regarding the role of local phenomena and heterogeneity in grain by grain. These challenges in turn necessitate understanding the fundamental mechanisms operating at individual grains and exploring how these mechanisms give rise to the collective behaviors of the films. Although numerous efforts have targeted local phenomena in these materials using classical scanning probe microscopy modalities, the measurements remain difficult to interpret due to intrinsic coupling between possible ferroelectric, ion migration, bulk and surface electrochemical phenomena, and the effects of light on these behaviors. Here, we explore the electromechanical responses in (Cs 0.05 FA 0.85 MA 0.15 )Pb(I 0.85 Br 0.15 ) 3 perovskite thin film in the dark and under illumination using band excitation-piezoresponse force microscopy (BE-PFM). To identify and separate ferroic, electrostatic, and electrochemical responses, we utilize contact-Kelvin probe force microscopy (cKPFM). The cKPFM responses show weak variability with morphological features in the dark, while such variability significantly increases under illumination. Deconvolution of the cKPFM responses via principal component analysis (PCA) shows heterogeneity in grain to grain possibly due to variation in migration and segregation of halide ions in mixed perovskites at specific grains. Further, we propose an electrostatic mechanism based on the physics of semiconductors, which can directly link the observed cKPFM behaviors with the ionic and electronic density of states (DOS) and their time dynamics. Here, these findings suggest that cKPFM is a powerful probe of local ionic behaviors on the nanoscale, both in HOIPs and other materials.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Search-based model identification of smart-structure damage

This paper describes the use of a combined model and parameter identification approach, based on modal analysis and artificial intelligence (AI) techniques, for identifying damage or flaws in a rotating truss structure incorporating embedded piezoceramic sensors. This smart structure example is representative of a class of structures commonly found in aerospace systems and next generation space structures. Artificial intelligence techniques of classification, heuristic search, and an object-oriented knowledge base are used in an AI-based model identification approach. A finite model space is classified into a search tree, over which a variant of best-first search is used to identify the model whose stored response most closely matches that of the input. Newly-encountered models can be incorporated into the model space. This adaptativeness demonstrates the potential for learning control. Following this output-error model identification, numerical parameter identification is used to further refine the identified model. Given the rotating truss example in this paper, noisy data corresponding to various damage configurations are input to both this approach and a conventional parameter identification method. The combination of the AI-based model identification with parameter identification is shown to lead to smaller parameter corrections than required by the use of parameter identification alone.

Glass, B. J.↗

Rotor Design Options for Improving V-22 Whirl-Mode Stability

With current tilt-rotor technology, very stiff, thick wings of limited aspect ratio are essential to meet aeroelastic stability requirements, which severely limits cruise efficiency and maximum speed. Larger and more efficient tiltrotors; will need longer and lighter wings, for which whirl-mode flutter is a serious design issue. Numerous approaches to improving the whirl-mode airspeed boundary have been investigated, typically focusing on wing design, active stability augmentation, and variable geometry rotors. The research to be reported applies the much simpler approaches of sweeping the outboard blade sections and extending balance weights at the tip. Improving proprotor whirl-mode stability margins is an ongoing research activity at NASA Ames Research Center. Previous efforts focussed on the XV-15 rotor. This is the first application of these design principles to a full V-22 rotor model in CAMRAD II (a previous publication presented results for a preliminary analytical model). The analysis to be discussed in the proposed paper includes a dual load-path blade model, a compressible-flow aerodynamic model, a representation of swashplate/rotor modal coupling, and other new features. The paper will include discussions of the mechanisms by which blade sweep and tip-mass offsets improve whirl-flutter margins. The V-22 was modeled with CAMRAD II; the rotor model is shown. It is also shows an example swept blade derived from the V-22 rotor system. For this rotor, the primary significance of sweep is the improved whirl-flutter boundary, not the reduced Mach-number effects. An offset tip mass is also shown.

Acree, C. W., Jr.↗

Machine tool cross beam design, fabrication, and testing using metal big area additive manufacturing

This paper describes the application of metal Big Area Additive Manufacturing (mBAAM) to the fabrication of a machine tool cross beam. The replacement of a traditional box design weldment with a new design printed by wire arc additive manufacturing using the MedUSA system at Oak Ridge National Laboratory (ORNL) is detailed. This requires a new design strategy based on the unique mBAAM capabilities. The intent of the new design is to reduce mass, while maintaining the dynamic stiffness. To compare the two designs, the natural frequencies and mode shapes are measured using impact testing and predicted using finite element analysis. It is confirmed that the printed structure dynamics agreed with the numerical model predictions, which demonstrates that it is feasible to model a large-scale mBAAM part and understand its behavior prior to printing. Another notable outcome of this study is that the significant residual stress and distortion in the print indicate that knowledge gaps remain for widespread implementation of mBAAM.

42 ENGINEERING↗

Modeling and Simulation for Exercise Vibration Isolation and Stabilization System Design

The microgravity environment that crew members experience on orbit presents a well-known health challenge, particularly when it comes to loss of muscle and bone mass. To counteract these negative effects, exercise countermeasures play a critical role in the daily routine of the crew on the International Space Station (ISS). To help inform requirements for upcoming exploration missions such as the Gateway Program, a new device, called the European Enhanced Exploration Exercise Device (E4D), is being built by the European Space Agency through their contractor, the Danish Aerospace Company. The E4D is being demonstrated on the ISS and is unique from other current exercise devices in that it provides four separate modalities in a single device: resistive, cycle ergometry, seated aerobic rowing, and rope pulling. To support the integration of E4D on ISS, a passive Vibration Isolation and Stabilization (VIS) system was required by NASA, and this responsibility was given to the Johnson Space Center. This paper describes the end-to-end process of modeling, simulation, and analysis used to inform the mechanical design of the VIS system. The process begins with the collection of representative exerciser motion capture (MoCap) through ground-based testing with the developmental E4D in both the Prototype Immersive Technology (PIT) Laboratory and the Active Response Gravity Offload System (ARGOS) facility at the Johnson Space Center. These collected MoCap data are processed through human biomechanics modeling to create forcing functions as input to a multibody dynamics simulation of the combined E4D/VIS system, with numerous resulting outputs. These outputs include microgravity accelerations, overall system displacements, internal and transmitted loads, as well as collisions. Microgravity accelerations are compared for compliance against ISS requirements while displacements are used for sway space determinations on the design and volumetric constraints within the targeted module. Internal loads are supplied to the supporting stress analysis teams and external loads for structural loads and dynamics teams, both at NASA and ESA. Finally, contact and clearance analysis is performed using the simulation to eliminate potential design issues. To ensure that the elements of the multibody simulation were verified and validated, correlation against multiple VIS related ground hardware testbeds was performed and characterized. In addition to the isolation part of the VIS problem, stabilization is also key to the integrated performance and evaluation. Due to the difficulty in defining ISS requirements in this area, the stability of the exerciser was inspected via analysis. Loss of balance was defined analytically as occurring when the resultant force vector acting on the exerciser lies outside the base of support of the feet. Both the VIS and E4D teams have recently gone through their Critical Design Reviews (CDRs) and the iterative model-based approach has been integral to inform mechanical design, particularly in the case of the VIS. This same end-to-end approach is now being applied for the Gateway Program, where an Exploration Exercise Device (EED) derived from the E4D and notional VIS for the device are under concept development.

Countermeasures↗

Wake Instability Behind Isolated Trip Near the Leading Edge of the BOLT-II Configuration

The BOLT-II (Holden Mission) configuration is an extended version of the BOLT flight article and depicts hypersonic boundary-layer transition in the presence of multiple and potentially interacting instability mechanisms. Several numerical studies of the boundary-layer instabilities over these configurations have been reported in the recent literature, including our previous studies of the modal instability characteristics of boundary-layer streaks adjacent to the minor-axis symmetry plane of the BOLT configuration and the wake instabilities behind a diamond planform (“pizza-box”) trip along the symmetry plane on the secondary side of the BOLT-II configuration. The present work extends the latter study to a scaled version of the same trip that is located in the region of nonzero crossflow in the vicinity of the leading edge at X/L = 0.5. Collectively, the trips along the symmetry plane and near the leading edge are the focus of the NASA roughness experiment on the secondary side of the BOLT-II configuration. The laminar basic state computation at the nominal flight design condition of Re ∞ = 5.44 x 10 6 /m and Re ∞ = 2.5 x 10 6 /m shows that the leading-edge trip with k/δ ≈ 0.70 and planform-side-length-to-height ratio of b/k = 3.0 induces multiple asymmetric, longitudinal streaks within the trip wake. The most prominent streak among these resembles a finite amplitude crossflow vortex, and it supports the amplification of multiple families of unstable modes. The application of multi-dimensional instability analysis to the wake flow reveals that the most amplified unstable mode can achieve a peak N-factor of up to 20 by X/L = 0.80, indicating that the onset of transition is more than likely to occur within eighty percent of the model length. To the best of our knowledge, the present study represents the first analysis including nonparallel and curvature effects on hypersonic tripwake instabilities in the presence of boundary-layer crossflow over a three-dimensional configuration.

Boundary layer transition↗

Monitoring of mechanical changes in a pipe assembly with complex geometry using multi-mode acoustic signals

Guided wave ultrasonic testing has been widely used for nondestructive and long-range evaluation of diverse mechanical structures, such as pipes and plates. In general, a large number of guided wave modes exist even in the simplest geometries, and they exhibit highly dispersive characteristics. Thus, sophisticated excitation methods and elaborate modal analyses need to be implemented to excite, acquire, and process well-defined pure modes as part of such a measurement method. If the interrogated structure has complex geometry, for example, with cross-sectional variations or multiple joints, the utilization of pure modes becomes significantly more difficult. To overcome this difficulty, this paper proposes the use of multi-mode, multi-frequency guided waves for long-range inspection of structures with arbitrarily complex geometries. Furthermore, using broadband chirp signals and multi-mode time–frequency analysis in a prototypical complex pipe assembly that includes multiple pipe diameters, elbows, and flange connections, we experimentally show sensitive detection of mechanical changes (both material addition and loss). In addition, we demonstrate that the complexity of the structure can be modeled in the digital domain with sufficient detail to allow accurate numerical simulation of the acoustic response using this method.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Neural network surrogate models for equations of state

Equation of state (EOS) data provide necessary information for accurate multiphysics modeling, which is necessary for fields such as inertial confinement fusion. Here, we suggest a neural network surrogate model of energy and entropy and use thermodynamic relationships to derive other necessary thermodynamic EOS quantities. We incorporate phase information into the model by training a phase classifier and using phase-specific regression models, which improves the modal prediction accuracy. Our model predicts energy values to 1% relative error and entropy to 3.5% relative error in a log-transformed space. Although sound speed predictions require further improvement, the derived pressure values are accurate within 10% relative error. Our results suggest that neural network models can effectively model EOS for inertial confinement fusion simulation applications.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Buckling and vibration analysis of laminated panels using VICONOPT

The analysis aspects of the 23,000-line FORTRAN program VIPASA with constraints and optimization (VICONOPT) are described. Overall stiffness matrices assembled from the earlier exact vibration and instability of plate assemblies including shear and anisotropy (VIPASA) flat plate stiffnesses are optionally coupled by Langrangian multipliers to find critical buckling loads, or natural frequencies of undamped vibration, of prismatic assemblies of anisotropic flat plates with arbitrarily located point supports or simple transverse supporting frames. The longitudunal continuity of typical wing and fuselage panels is closely approximated because the solutions are for the infinitely long structure obtained by repeating a bay and its supports longitudinally. Any longitudinally invariant in-plane plate stresses are permitted, and very rapid solutions are guaranteed by numerous refinements, including multilevel substructuring and a method for repetitive cross sections that is exact for regular polygons used to represent cylinders. Modal displacements and stresses in or between plies of laminated plates are calculated and plotted, with values being recovered at all nodes of substructures. Comparison with usual approximate finite-element methods confirms that, for comparably converged solutions, VICONOPT is typically between 100 and 10(exp 4) times faster.

Kennedy, David↗

Ground Vibration Test Planning and Pre-Test Analysis for the X-33 Vehicle

This paper describes the results of the modal test planning and the pre-test analysis for the X-33 vehicle. The pre-test analysis included the selection of the target modes, selection of the sensor and shaker locations and the development of an accurate Test Analysis Model (TAM). For target mode selection, four techniques were considered, one based on the Modal Cost technique, one based on Balanced Singular Value technique, a technique known as the Root Sum Squared (RSS) method, and a Modal Kinetic Energy (MKE) approach. For selecting sensor locations, four techniques were also considered; one based on the Weighted Average Kinetic Energy (WAKE), one based on Guyan Reduction (GR), one emphasizing engineering judgment, and one based on an optimum sensor selection technique using Genetic Algorithm (GA) search technique combined with a criteria based on Hankel Singular Values (HSV's). For selecting shaker locations, four techniques were also considered; one based on the Weighted Average Driving Point Residue (WADPR), one based on engineering judgment and accessibility considerations, a frequency response method, and an optimum shaker location selection based on a GA search technique combined with a criteria based on HSV's. To evaluate the effectiveness of the proposed sensor and shaker locations for exciting the target modes, extensive numerical simulations were performed. Multivariate Mode Indicator Function (MMIF) was used to evaluate the effectiveness of each sensor & shaker set with respect to modal parameter identification. Several TAM reduction techniques were considered including, Guyan, IRS, Modal, and Hybrid. Based on a pre-test cross-orthogonality checks using various reduction techniques, a Hybrid TAM reduction technique was selected and was used for all three vehicle fuel level configurations.

Bedrossian, Herand↗

Natural Frequency Testing and Model Correlation of Rocket Engine Structures in Liquid Hydrogen - Phase I, Cantilever Beam

Many structures in the launch vehicle industry operate in liquid hydrogen (LH2), from the hydrogen fuel tanks through the ducts and valves and into the pump sides of the turbopumps. Calculating the structural dynamic response of these structures is critical for successful qualification of this hardware, but accurate knowledge of the natural frequencies is based entirely on numerical or analytical predictions of frequency reduction due to the added-fluid-mass effect because testing in LH2 has always been considered too difficult and dangerous. This fluid effect is predicted to be approximately 4-5% using analytical formulations for simple cantilever beams. As part of a comprehensive test/analysis program to more accurately assess pump inducers operating in LH2, a series of frequency tests in LH2 were performed at NASA/Marshall Space Flight Center's unique cryogenic test facility. These frequency tests are coupled with modal tests in air and water to provide critical information not only on the mass effect of LH2, but also the cryogenic temperature effect on Young's Modulus for which the data is not extensive. The authors are unaware of any other reported natural frequency testing in this media. In addition to the inducer, a simple cantilever beam was also tested in the tank to provide a more easily modeled geometry as well as one that has an analytical solution for the mass effect. This data will prove critical for accurate structural dynamic analysis of these structures, which operate in a highly-dynamic environment.

Brown, Andrew M.↗

Predicted vibration responses of Apollo structure and effects of pressure correlation lengths on response

Progressive wove and reverberant acoustic fields ore often used in the laboratory for qualification testing of aerospace structures that are subjected to unsteady aerodynamic environments which occur during flight through the atmosphere. Acoustic simulation of these in-flight environments i s achieved i f the structural vibration response to acoustic excitation is equivalent, over the frequency range of interest, to the response caused by aerodynamic turbulence. To determine the degree of response simulation that could be expected for such qualification testing of the Apollo Spacecraft, responses of a segment of the Apollo structure were estimated theoretically for both acoustic and flight environments. The vibration analysis presented in this report was performed for the Spacecraft Lunar Module Adaptor; and the three environments treated are turbulence at Mach 1 and Mach 2, a reverberant acoustic field, and a modified acoustic progressive wave field. The progressive wave field is assumed to be generated within a specially designed shroud of sixteen axially oriented ducts which is constructed around and contoured to the external skin of the vehicle. Each duct is open on the side adjacent to the skin and is driven acoustically at one end by an independent noise source, thus allowing for different pressure correlation patterns around the circumference of the vehicle. The purpose of the analysis was to determine the effects of surface pressure correlation lengths on the response characteristics of representative Apollo structure, and to select on this basis, the optimum acoustic environment for structural vibration qualification tests of the Apollo in the Spacecraft Acoustic Laboratory at MSC. The method of analysis employed consists of determining, independently, the space average, mean-square acceleration spectral density for each of the significant modes of vibration of an equivalent cylindrical shell, and then expressing the total structural response as a linear summation of these modal responses. Approximately 570 classical modes of a pinned-end cylinder were used in order to find the spectral response characteristics of the structure throughout the frequency range of 10 - 1,000 Hz. The effective forcing functions or joint acceptances of the various modes of vibration were computed and are presented graphically for each of the pressure excitation environments. A brief development of the response equations used is presented in the report. Space average acceleration response spectra were computed for twelve different cases including aerodynamic turbulence at Mach 1 and Mach 2, reverberant acoustic field, five duct correlation patterns, two cases of an axially damped progressive wave field, and two cases in which the structural damping constant was varied from the expected value. The response spectra were computed and plotted by a high speed digital computer far a flat excitation pressure spectrum; and the response spectra were numerically integrated to give one-third octave bond average responses. A complete set of digital computer programs has been developed for both flat plate and cylindrical shell structures for the above environments, and each requires run times of three minutes or less. Thus, the work presented herein, along with the results, constitutes a practical advancement to the existing state-of-the-art of vibration prediction.

R W White↗

Model fusion with physics-guided machine learning: Projection-based reduced-order modeling

The unprecedented amount of data generated from experiments, field observations, and large-scale numerical simulations at a wide range of spatiotemporal scales has enabled the rapid advancement of data-driven and especially deep learning models in the field of fluid mechanics. Although these methods are proven successful for many applications, there is a grand challenge of improving their generalizability. This is particularly essential when data-driven models are employed within outer-loop applications like optimization. In this work, we put forth a physics-guided machine learning (PGML) framework that leverages the interpretable physics-based model with a deep learning model. Leveraging a concatenated neural network design from multi-modal data sources, the PGML framework is capable of enhancing the generalizability of data-driven models and effectively protects against or inform about the inaccurate predictions resulting from extrapolation. We apply the PGML framework as a novel model fusion approach combining the physics-based Galerkin projection model and long- to short-term memory (LSTM) network for parametric model order reduction of fluid flows. We demonstrate the improved generalizability of the PGML framework against a purely data-driven approach through the injection of physics features into intermediate LSTM layers. Our quantitative analysis shows that the overall model uncertainty can be reduced through the PGML approach, especially for test data coming from a distribution different than the training data. Moreover, we demonstrate that our approach can be used as an inverse diagnostic tool providing a confidence score associated with models and observations. The proposed framework also allows for multi-fidelity computing by making use of low-fidelity models in the online deployment of quantified data-driven models.

42 ENGINEERING↗

Dynamics retrieval from stochastically weighted incomplete data by low-pass spectral analysis

Time-resolved serial femtosecond crystallography (TR-SFX) provides access to protein dynamics on sub-picosecond timescales, and with atomic resolution. Due to the nature of the experiment, these datasets are often highly incomplete and the measured diffracted intensities are affected by partiality. To tackle these issues, one established procedure is that of splitting the data into time bins, and averaging the multiple measurements of equivalent reflections within each bin. This binning and averaging often involve a loss of information. Here, we propose an alternative approach, which we call low-pass spectral analysis (LPSA). In this method, the data are projected onto the subspace defined by a set of trigonometric functions, with frequencies up to a certain cutoff. This approach attenuates undesirable high-frequency features and facilitates retrieving the underlying dynamics. A time-lagged embedding step can be included prior to subspace projection to improve the stability of the results with respect to the parameters involved. Subsequent modal decomposition allows to produce a low-rank description of the system's evolution. Using a synthetic time-evolving model with incomplete and partial observations, we analyze the LPSA results in terms of quality of the retrieved signal, as a function of the parameters involved. We compare the performance of LPSA to that of a range of other sophisticated data analysis techniques. We show that LPSA allows to achieve excellent dynamics reconstruction at modest computational cost. Finally, we demonstrate the superiority of dynamics retrieval by LPSA compared to time binning and merging, which is, to date, the most commonly used method to extract dynamical information from TR-SFX data.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Modeling, Modal Properties, and Mesh Stiffness Variation Instabilities of Planetary Gears

Planetary gear noise and vibration are primary concerns in their applications in helicopters, automobiles, aircraft engines, heavy machinery and marine vehicles. Dynamic analysis is essential to the noise and vibration reduction. This work analytically investigates some critical issues and advances the understanding of planetary gear dynamics. A lumped-parameter model is built for the dynamic analysis of general planetary gears. The unique properties of the natural frequency spectra and vibration modes are rigorously characterized. These special structures apply for general planetary gears with cyclic symmetry and, in practically important case, systems with diametrically opposed planets. The special vibration properties are useful for subsequent research. Taking advantage of the derived modal properties, the natural frequency and vibration mode sensitivities to design parameters are investigated. The key parameters include mesh stiffnesses, support/bearing stiffnesses, component masses, moments of inertia, and operating speed. The eigen-sensitivities are expressed in simple, closed-form formulae associated with modal strain and kinetic energies. As disorders (e.g., mesh stiffness variation. manufacturing and assembling errors) disturb the cyclic symmetry of planetary gears, their effects on the free vibration properties are quantitatively examined. Well-defined veering rules are derived to identify dramatic changes of natural frequencies and vibration modes under parameter variations. The knowledge of free vibration properties, eigen-sensitivities, and veering rules provide important information to effectively tune the natural frequencies and optimize structural design to minimize noise and vibration. Parametric instabilities excited by mesh stiffness variations are analytically studied for multi-mesh gear systems. The discrepancies of previous studies on parametric instability of two-stage gear chains are clarified using perturbation and numerical methods. The operating conditions causing parametric instabilities are expressed in closed-form suitable for design guidance. Using the well-defined modal properties of planetary gears, the effects of mesh parameters on parametric instability are analytically identified. Simple formulae are obtained to suppress particular instabilities by adjusting contact ratios and mesh phasing.

Parker, Robert G.↗