Parametric studies on dynamic stiffness of ball bearings supporting a flexible rotor

2019 ◽  
Vol 25 (15) ◽  
pp. 2175-2188 ◽  
Author(s):  
Tikam Chand Gupta
Keyword(s):  
Numerical Integration ◽  
Dynamic Stiffness ◽  
Mathematical Formulation ◽  
Rolling Element Bearing ◽  
Radial Clearance ◽  
Integration Scheme ◽  
Time Step ◽  
Average Value ◽  
Bearing Stiffness ◽  
Time Invariant

Most of the researchers in the field of dynamics of the rolling element bearing have considered bearing stiffness as time invariant and/or not related to dynamics of the bearing. In the present paper, the bearing stiffness has been taken as function of dynamic response at every time step of numerical simulation and a detailed parametric study is performed to investigate the effect of flexibility of the rotor shaft, rotational speed, and internal radial clearance on the instantaneous and average value of dynamic stiffness of the ball bearing. The mathematical formulation is based on the Timoshenko beam finite element theory. Gravity and bearing forces are considered as external forces acting on a free-free flexible shaft. A stable Newmark- β numerical integration scheme coupled with Newton–Raphson method is used for numerical integration and for convergence to an accurate value of bearing stiffness. The results showing the variation of different components of bearing stiffness as a function of time-invariant parameters has improved the understanding of the dynamic behavior of the bearing during motion. The variation pattern of bearing stiffness coefficients is observed to be sensitive to direction of rotation. The amplitude of periodic change of these coefficients increases with the increase of the stiffness ratio of shaft and the decrease of radial clearance.

Numerical Integration of Some Stiff Constitutive Models of Inelastic Deformation

1980 ◽  
Vol 102 (1) ◽  
pp. 92-96 ◽  
Author(s):  
Virendra Kumar ◽  
Mahesh Morjaria ◽  
Subrata Mukherjee
Keyword(s):  
Numerical Integration ◽  
Inelastic Deformation ◽  
Time Integration ◽  
Constitutive Models ◽  
Integration Scheme ◽  
Uniaxial Deformation ◽  
Time Step ◽  
One Step ◽  
Numerical Time Integration ◽  
Step Control

Several strategies for numerical time-integration of some stiff constitutive models of inelastic deformation are presented in this paper. Numerical results and comparisons are presented for the integration of one such model for the case of uniaxial deformation under various prescribed histories of stress or strain. A simple one step Euler type integration scheme with automatic time-step control, which can be easily adapted to the solution of multiaxial boundary value problems, appears promising.

Simulations of Gear-Rolling Element Bearing Interactions in Geared Transmissions

2003 ◽  
Author(s):  
F. Lahmar ◽  
P. Velex
Keyword(s):  
Equations Of Motion ◽  
Rolling Element Bearing ◽  
Integration Scheme ◽  
Time Step ◽  
Normal Contact ◽  
Contact Algorithm ◽  
Rolling Element ◽  
Non Linear ◽  
External Forces ◽  
Gear Rolling

The modular model of geared systems presented in this paper makes it possible to simultaneously account for the contact conditions in gears and rolling element bearings. Gears are modeled as two rigid cylinders connected by distributed mesh stiffnesses while ball and roller bearings contribute to the equations of motion as time-varying, non-linear external forces. Solutions are obtained by combining a Newmark time-step integration scheme, a Newton-Raphson method for ball bearing non-linearity and a normal contact algorithm that deals with the contact problem between the teeth. It is found that the static gear-bearing couplings are generally more important than the dynamic couplings with a significant influence of the gear on the bearing response. Finally, it is shown that, in certain conditions, bearings can generate non-linear parametric excitations of the same orders of magnitude as those associated with the meshing of helical gears.

scholarly journals Regional and seasonal truncation errors of trajectory calculations using ECMWF high-resolution operational analyses and forecasts

2017 ◽  
Author(s):  
Thomas Rößler ◽  
Olaf Stein ◽  
Yi Heng ◽  
Lars Hoffmann
Keyword(s):  
High Resolution ◽  
Numerical Integration ◽  
Transport Processes ◽  
Integration Scheme ◽  
Free Troposphere ◽  
Time Step ◽  
Truncation Errors ◽  
Trajectory Calculations ◽  
Integration Schemes ◽  
Runge Kutta

Abstract. Lagrangian particle dispersion models (LPDMs) are indispensable tools to study atmospheric transport processes. The accuracy of trajectory calculations, which form an essential part of LPDM simulations, depends on various factors. Here we focus on truncation errors that originate from the use of numerical integration schemes to solve the kinematic equation of motion. The optimization of numerical integration schemes to minimize truncation errors and to maximize computational speed is of great interest regarding the computational efficiency of large-scale LPDM simulations. In this study we analyzed truncation errors of six explicit integration schemes of the Runge Kutta family, which we implemented in the Massive-Parallel Trajectory Calculations (MPTRAC) model. The simulations were driven by wind fields of the latest operational analysis and forecasts of the European Centre for Medium-range Weather Forecasts (ECMWF) at T1279L137 spatial resolution and 3 h temporal sampling. We defined separate test cases for 15 distinct domains of the atmosphere, covering the polar regions, the mid-latitudes, and the tropics in the free troposphere, in the upper troposphere and lower stratosphere (UT/LS) region, and in the lower and mid stratosphere. For each domain we performed simulations for the months of January, April, July, and October for the years of 2014 and 2015. In total more than 5000 different transport simulations were performed. We quantified the accuracy of the trajectories by calculating transport deviations with respect to reference simulations using a 4th-order Runge-Kutta integration scheme with a sufficiently fine time step. We assessed the transport deviations with respect to error limits based on turbulent diffusion. Independent of the numerical scheme, the truncation errors vary significantly between the different domains and seasons. Especially the differences in altitude stand out. Horizontal transport deviations in the stratosphere are typically an order of magnitude smaller compared with the free troposphere. We found that the truncation errors of the six numerical schemes fall into three distinct groups, which mostly depend on the numerical order of the scheme. Schemes of the same order differ little in accuracy, but some methods need less computational time, which gives them an advantage in efficiency. The selection of the integration scheme and the appropriate time step should possibly take into account the typical altitude ranges as well as the total length of the simulations to achieve the most efficient simulations. However, trying to generalize, we recommend the 3rd-order Runge Kutta method with a time step of 170 s or the midpoint scheme with a time step of 100 s for efficient simulations of up to 10 days time based on ECMWF's high-resolution meteorological data.

Application of the Modified Rosenbrock Algorithm for Transient Rotordynamic Analysis

2012 ◽  
Author(s):  
Anand Srinivasan ◽  
Dhruv Kumar
Keyword(s):  
Numerical Integration ◽  
Equations Of Motion ◽  
Mathematical Formulation ◽  
Single Step ◽  
Degree Of Freedom ◽  
Step Size ◽  
Time Step ◽  
Time Step Size ◽  
Forcing Function ◽  
Transient Simulations

It is well known that transient rotordynamic analyses involve numerical integration of the equations of motion in order to study the response of the system under an applied forcing function. A common problem that arises in such simulations is the choice of step-size that needs to be used to obtain numerically stable results. Traditional numerical integration techniques such as the Runge-Kutta algorithms not only require splitting up second order differential equations as two first order equations, but also necessitate multiple integrations at each time-step, thus increasing the solution time. The Newmark-beta and Wilson-theta algorithms are some of the prevalent methods that have been used for transient simulations in rotordynamics. However, those single-step methods are only conditionally stable, and require iterations to converge to a solution at each time step, thus making it pseudosingle-step. In the more recent years, a modified form of the Rosenbrock algorithm has been proposed as a numerically stable and true single-step mathematical formulation for the integration of structural dynamics problems. In this paper, the modified Rosenbrock algorithm has been applied to a transient start-up multi-degree-of-freedom rotordynamics problem. A constant time step-size algorithm has been used for the simulations, and results of the transient analysis have been presented. The fact that a multi-degree-of-freedom system can be solved without condensation of the higher order modes makes the superior numerical damping characteristics of the algorithm become evident.

scholarly journals Numerical Integration Scheme Using Singular Perturbation Method

2013 ◽  
Author(s):  
Gibin Gil ◽  
Ricardo G. Sanfelice ◽  
Parviz E. Nikravesh
Keyword(s):  
Singular Perturbation ◽  
Perturbation Method ◽  
Numerical Integration ◽  
Equations Of Motion ◽  
Computation Time ◽  
Integration Scheme ◽  
Singular Perturbation Method ◽  
Time Step ◽  
Numerical Integration Scheme ◽  
Deformable Bodies

Some multi degree-of-freedom dynamical systems exhibit a response that contain fast and slow variables. An example of such systems is a multibody system with rigid and deformable bodies. Standard numerical integration of the resultant equations of motion must adjust the time step according to the frequency of the fastest variable. As a result, the computation time is sacrificed. The singular perturbation method is an analysis technique to deal with the interaction of slow and fast variables. In this study, a numerical integration scheme using the singular perturbation method is discussed, its absolute stability condition is derived, and its order of accuracy is investigated.

scholarly journals Three-Dimensional Elastodynamic Analysis Employing Partially Discontinuous Boundary Elements

Algorithms ◽  
2021 ◽  
Vol 14 (5) ◽  
pp. 129
Author(s):  
Yuan Li ◽  
Ni Zhang ◽  
Yuejiao Gong ◽  
Wentao Mao ◽  
Shiguang Zhang
Keyword(s):  
Side Effect ◽  
Domain Boundary ◽  
Three Dimensional ◽  
Boundary Integral ◽  
Integration Scheme ◽  
Computational Accuracy ◽  
Time Step ◽  
Corner Points ◽  
Discontinuous Elements ◽  
The Time Domain

Compared with continuous elements, discontinuous elements advance in processing the discontinuity of physical variables at corner points and discretized models with complex boundaries. However, the computational accuracy of discontinuous elements is sensitive to the positions of element nodes. To reduce the side effect of the node position on the results, this paper proposes employing partially discontinuous elements to compute the time-domain boundary integral equation of 3D elastodynamics. Using the partially discontinuous element, the nodes located at the corner points will be shrunk into the element, whereas the nodes at the non-corner points remain unchanged. As such, a discrete model that is continuous on surfaces and discontinuous between adjacent surfaces can be generated. First, we present a numerical integration scheme of the partially discontinuous element. For the singular integral, an improved element subdivision method is proposed to reduce the side effect of the time step on the integral accuracy. Then, the effectiveness of the proposed method is verified by two numerical examples. Meanwhile, we study the influence of the positions of the nodes on the stability and accuracy of the computation results by cases. Finally, the recommended value range of the inward shrink ratio of the element nodes is provided.

Influence of Clearances and Thermal Effects on the Dynamic Behavior of Gear-Hydrodynamic Journal Bearing Systems

2013 ◽  
Vol 135 (6) ◽  
Author(s):  
R. Fargère ◽  
P. Velex
Keyword(s):  
Degrees Of Freedom ◽  
Integration Scheme ◽  
Time Step ◽  
Specific Element ◽  
Normal Contact ◽  
Contact Algorithm ◽  
Gear Teeth ◽  
Proposed Model ◽  
Elastic Couplings ◽  
Definition Of

A global model of mechanical transmissions is introduced which deals with most of the possible interactions between gears, shafts, and hydrodynamic journal bearings. A specific element for wide-faced gears with nonlinear time-varying mesh stiffness and tooth shape deviations is combined with shaft finite elements, whereas the bearing contributions are introduced based on the direct solution of Reynolds' equation. Because of the large bearing clearances, particular attention has been paid to the definition of the degrees-of-freedom and their datum. Solutions are derived by combining a time step integration scheme, a Newton–Raphson method, and a normal contact algorithm in such a way that the contact conditions in the bearings and on the gear teeth are simultaneously dealt with. A series of comparisons with the experimental results obtained on a test rig are given which prove that the proposed model is sound. Finally, a number of results are presented which show that parameters often discarded in global models such as the location of the oil inlet area, the oil temperature in the bearings, the clearance/elastic couplings interactions, etc. can be influential on static and dynamic tooth loading.

scholarly journals Unbalance Response Prediction for Rotors on Ball Bearings Using Speed- and Load-Dependent Nonlinear Bearing Stiffness

2005 ◽  
Vol 2005 (1) ◽  
pp. 53-59 ◽  
Author(s):  
David P. Fleming ◽  
J. V. Poplawski
Keyword(s):  
Response Prediction ◽  
Rolling Element Bearing ◽  
Response Analysis ◽  
Ball Bearings ◽  
Moderate Amount ◽  
Constant Stiffness ◽  
Unbalance Response ◽  
Rolling Element ◽  
Bearing Stiffness ◽  
Element Bearing

Rolling-element bearing forces vary nonlinearly with bearing deflection. Thus, an accurate rotordynamic analysis requires that bearing forces corresponding to the actual bearing deflection be utilized. For this work, bearing forces were calculated by COBRA-AHS, a recently developed rolling-element bearing analysis code. Bearing stiffness was found to be a strong function of bearing deflection, with higher deflection producing markedly higher stiffness. Curves fitted to the bearing data for a range of speeds and loads were supplied to a flexible rotor unbalance response analysis. The rotordynamic analysis showed that vibration response varied nonlinearly with the amount of rotor imbalance. Moreover, the increase in stiffness as critical speeds were approached caused a large increase in rotor and bearing vibration amplitude over part of the speed range compared to the case of constant-stiffness bearings. Regions of bistable operation were possible, in which the amplitude at a given speed was much larger during rotor acceleration than during deceleration. A moderate amount of damping will eliminate the bistable region, but this damping is not inherent in ball bearings.

Vibration analysis of multiple degrees of freedom mechanical system with dry friction

2012 ◽  
Vol 227 (7) ◽  
pp. 1505-1514 ◽  
Author(s):  
SD Yu ◽  
BC Wen
Keyword(s):  
Mechanical System ◽  
Dry Friction ◽  
Degrees Of Freedom ◽  
Simple Procedure ◽  
Mathematical Problem ◽  
Equations Of Motion ◽  
Inequality Constraints ◽  
Integration Scheme ◽  
Time Step ◽  
Multiple Degrees Of Freedom

This article presents a simple procedure for predicting time-domain vibrational behaviors of a multiple degrees of freedom mechanical system with dry friction. The system equations of motion are discretized by means of the implicit Bozzak–Newmark integration scheme. At each time step, the discontinuous frictional force problem involving both the equality and inequality constraints is successfully reduced to a quadratic mathematical problem or the linear complementary problem with the introduction of non-negative and complementary variable pairs (supremum velocities and slack forces). The so-obtained complementary equations in the complementary pairs can be solved efficiently using the Lemke algorithm. Results for several single degree of freedom and multiple degrees of freedom problems with one-dimensional frictional constraints and the classical Coulomb frictional model are obtained using the proposed procedure and compared with those obtained using other approaches. The proposed procedure is found to be accurate, efficient, and robust in solving non-smooth vibration problems of multiple degrees of freedom systems with dry friction. The proposed procedure can also be applied to systems with two-dimensional frictional constraints and more sophisticated frictional models.

scholarly journals Investigation of Air Pocket Behavior in Pipelines Using Rigid Column Model and Contributions of Time Integration Schemes

Water ◽  
2021 ◽  
Vol 13 (6) ◽  
pp. 785
Author(s):  
Arman Rokhzadi ◽  
Musandji Fuamba
Keyword(s):  
Transient Flow ◽  
Time Integration ◽  
Driving Pressure ◽  
Implicit Time ◽  
Numerical Dissipation ◽  
Integration Scheme ◽  
Time Step ◽  
Column Model ◽  
Integration Schemes ◽  
Backward Euler

This paper studies the air pressurization problem caused by a partially pressurized transient flow in a reservoir-pipe system. The purpose of this study is to analyze the performance of the rigid column model in predicting the attenuation of the air pressure distribution. In this regard, an analytic formula for the amplitude and frequency will be derived, in which the influential parameters, particularly, the driving pressure and the air and water lengths, on the damping can be seen. The direct effect of the driving pressure and inverse effect of the product of the air and water lengths on the damping will be numerically examined. In addition, these numerical observations will be examined by solving different test cases and by comparing to available experimental data to show that the rigid column model is able to predict the damping. However, due to simplified assumptions associated with the rigid column model, the energy dissipation, as well as the damping, is underestimated. In this regard, using the backward Euler implicit time integration scheme, instead of the classical fourth order explicit Runge–Kutta scheme, will be proposed so that the numerical dissipation of the backward Euler implicit scheme represents the physical dissipation. In addition, a formula will be derived to calculate the appropriate time step size, by which the dissipation of the heat transfer can be compensated.

Sign in / Sign up

Export Citation Format

Share Document