Efficient Simulation of Gear Contacts Including Transient Elastohydrodynamic Effects

2013 ◽  
Vol 135 (3) ◽  
Author(s):  
Peter Fietkau ◽  
Bernd Bertsche

This paper describes an efficient transient elastohydrodynamic simulation method for gear contacts. The model uses oil films and elastic deformations directly in the multibody simulation, and is based on the Reynolds equation including squeeze and wedge terms as well as an elastic half-space. Two transient solutions to this problem, an analytical and a numerical one, were developed. The analytical solution is accomplished using assumptions for the gap shape and the pressure in the middle of the gap. The numerical problem is solved using multilevel multi-integration algorithms. With this approach, tooth impacts during gear rattling as well as highly loaded power-transmitting gear contacts can be investigated and lubrication conditions like gap heights or type of friction may be determined. The method was implemented in the multibody simulation environment SIMPACK. Therefore it is easy to transfer the developed element to other models and use it for a multitude of different engineering problems. A detailed three-dimensional elastic multibody model of an experimental transmission is used to validate the developed method. Important values of the gear contact like normal and tangential forces, proportion of dry friction, and minimum gap heights are calculated and studied for different conditions. In addition, pressure distributions on tooth flanks as well as gap forms are determined based on the numerical solution method. Finally, the simulation approach is validated with measurements and shows good consistency. The simulation model is therefore capable of predicting transient gear contact under different operating conditions such as load vibrations or gear rattling. Simulations of complete transmissions are possible and therefore a direct determination of transmission vibration behavior and structure-borne noise as well as of forces and lubrication conditions can be done.

Author(s):  
Yin Gao ◽  
Mike McHenry ◽  
Brad Kerchof

Cut spike fasteners, used with conventional AREMA rolled tie plates and solid sawn timber ties, are the most common tie and fastener system used on North American freight railroads. Cut spikes are also used to restrain tie plates that incorporate an elastic rail fastener — that is, an elastic clip that fastens the rail to the tie plate. Elastic fasteners have been shown to reduce gage widening and decrease the potential for rail roll compared to cut spike-only systems. For this reason, elastic fastener systems have been installed in high degree curves on many railroads. Recent observations on one Class I railroad have noted broken cut spikes when used with these types of tie plates in mountainous, high degree curve territory. Broken screw spikes and drive spikes on similar style plates have also been observed. In this paper, a simulation method that integrates a vehicle-track system dynamics model, NUCARS®, with a finite element analysis model is used to investigate the root causes of the broken spikes. The NUCARS model consists of a detailed multibody train, wheel-rail contact parameters, and track model that can estimate the dynamic loading environment of the fastening system. For operating conditions in tangent and curve track, this loading environment is then replicated in a finite element model of the track structure — ties, tie plates, and cut spikes. The stress contours of the cut spikes generated in these simulations are compared to how cut spikes have failed in revenue service. The tuning and characterization of both the vehicle dynamics multibody model and the finite element models are presented. Additionally, the application of this approach to other types of fastening systems and spike types is discussed. Preliminary results have identified a mechanism involving the dynamic unloading of the tie plate-to-tie interface due to rail uplift ahead of the wheel and the resulting transfer of net longitudinal and lateral forces into the cut spikes. Continued analysis will attempt to confirm this mechanism and will focus on the severity of these stresses, the effect of increased grade, longitudinal train dynamics, braking forces, and curvature.


2019 ◽  
Vol 9 (2) ◽  
pp. 192-197
Author(s):  
Somrita Ghosh ◽  
Aritra Acharyya

Background: The time and frequency responses of Multiple Quantum Barrier (MQB) nano-scale Avalanche Photodiodes (APDs) based on Si~3C-SiC material system have been investigated in this final part. Methods: A very narrow rectangular pulse of pulse-width of 0.4 ps has been used as the input optical pulse having 850 nm wavelength incidents on the p+-side of the MQB APD structures and corresponding current responses have been calculated by using a simulation method developed by the authors. Results: Finally the frequency responses of the devices are obtained via the Fourier transform of the corresponding pulse current responses in time domain. Conclusion: Simulation results show that MQB nano-APDs possess significantly faster time response and wider frequency response as compared to the flat Si nano-APDs under similar operating conditions.


2020 ◽  
Vol 12 (1) ◽  
pp. 703-717
Author(s):  
Yin Wei ◽  
Wang Jiaqi ◽  
Bai Xiaomin ◽  
Sun Wenjie ◽  
Zhou Zheyuan

AbstractThis article analyzes the technical difficulties in full-section backfill mining and briefly introduces the technical principle and advantages of backfilling combined with caving fully mechanized mining (BCCFM). To reveal the strata behavior law of the BCCFM workface, this work establishes a three-dimensional numerical model and designs a simulation method by dynamically updating the modulus parameter of the filling body. By the analysis of numerical simulation, the following conclusions about strata behavior of the BCCFM workface were drawn. (1) The strata behavior of the BCCFM workface shows significant nonsymmetrical characteristics, and the pressure in the caving section is higher than that in the backfilling section. φ has the greatest influence on the backfilling section and the least influence on the caving section. C has a significant influence on the range of abutment pressure in the backfilling section. (2) There exits the transition area with strong mine pressure of the BCCFM workface. φ and C have significant effect on the degree of pressure concentration but little effect on the influence range of strong mine pressure in the transition area. (3) Under different conditions, the influence range of strong mine pressure is all less than 6 m. This article puts forward a control strategy of mine pressure in the transition area, which is appropriately improving the strength of the transition hydraulic support within the influence range (6 m) in the transition area according to the pressure concentration coefficient. The field measurement value of Ji15-31010 workface was consistent with numerical simulation, which verifies the reliability of control strategy of the BCCFM workface.


1998 ◽  
Vol 11 (1) ◽  
pp. 570-570
Author(s):  
Johan Holmberg ◽  
Lennart Lindegren ◽  
Chris Flynn

We use the Hipparcos survey to derive an improved model of the local galactic structure. The availability of parallaxes for all the stars permits direct determination of stellar distributions, eliminating the basic indeterminacy of classical methods based on star counts. Hipparcos gives for the first time a truly three-dimensional view of the solar vicinity, and a complete, homogeneous and highly accurate set of magnitudes and colours. This means that new techniques can be applied in the treatment of the data which place strong constraints on a model that tries to describe the local Galactic structure. Here we investigate how well a static model of low complexitycan describe the Hipparcos observations. The interpretation of the Hipparcos data is complicated by various observational errors and selection effects that are hard to treat correctly. We do not try to correct the data, but instead use a model and subject this model to the same observational errors and selection effects. A model catalogue is created that can be compared with the observed catalogue directly in the observational domain, thereby eliminating the effects from various biases. Many features in the HR diagram are for the first time seen in field stars thanks to Hipparcos, such as the slanted red giant clump, previously seen in rich old open clusters such as Berkeley 18. This and other features ofthe observed HR diagram are well reproduced by the model thanks to the rather detailed modelling of the joint Mv/B — V distribution. Actually, separate distributions were derived for the three different components, disk, thick disk and halo, using the kinematic characteristics of the components to discriminate between them.


1971 ◽  
Vol 38 (4) ◽  
pp. 899-905 ◽  
Author(s):  
L. B. Freund

Three-dimensional wave propagation in an elastic half space is considered. The half space is traction free on half its boundary, while the remaining part of the boundary is free of shear traction and is constrained against normal displacement by a smooth, rigid barrier. A time-harmonic surface wave, traveling on the traction free part of the surface, is obliquely incident on the edge of the barrier. The amplitude and the phase of the resulting reflected surface wave are determined by means of Laplace transform methods and the Wiener-Hopf technique. Wave propagation in an elastic half space in contact with two rigid, smooth barriers is then considered. The barriers are arranged so that a strip on the surface of uniform width is traction free, which forms a wave guide for surface waves. Results of the surface wave reflection problem are then used to geometrically construct dispersion relations for the propagation of unattenuated guided surface waves in the guiding structure. The rate of decay of body wave disturbances, localized near the edges of the guide, is discussed.


2017 ◽  
Vol 2017 ◽  
pp. 1-9 ◽  
Author(s):  
Zhengzhi Wang ◽  
Chunling Zhu

In view of the rotor icing problems, the influence of centrifugal force on rotor blade icing is investigated. A numerical simulation method of three-dimensional rotor blade icing is presented. Body-fitted grids around the rotor blade are generated using overlapping grid technology and rotor flow field characteristics are obtained by solving N-S equations. According to Eulerian two-phase flow, the droplet trajectories are calculated and droplet impingement characteristics are obtained. The mass and energy conservation equations of ice accretion model are established and a new calculation method of runback water mass based on shear stress and centrifugal force is proposed to simulate water flow and ice shape. The calculation results are compared with available experimental results in order to verify the correctness of the numerical simulation method. The influence of centrifugal force on rotor icing is calculated. The results show that the flow direction and distribution of liquid water on rotor surfaces change under the action of centrifugal force, which lead to the increasing of icing at the stagnation point and the decreasing of icing on both frozen limitations.


2017 ◽  
Vol 34 (5) ◽  
pp. 1551-1571 ◽  
Author(s):  
Ming Xia

Purpose The main purpose of this paper is to present a comprehensive upscale theory of the thermo-mechanical coupling particle simulation for three-dimensional (3D) large-scale non-isothermal problems, so that a small 3D length-scale particle model can exactly reproduce the same mechanical and thermal results with that of a large 3D length-scale one. Design/methodology/approach The objective is achieved by following the scaling methodology proposed by Feng and Owen (2014). Findings After four basic physical quantities and their similarity-ratios are chosen, the derived quantities and its similarity-ratios can be derived from its dimensions. As the proposed comprehensive 3D upscale theory contains five similarity criteria, it reveals the intrinsic relationship between the particle-simulation solution obtained from a small 3D length-scale (e.g. a laboratory length-scale) model and that obtained from a large 3D length-scale (e.g. a geological length-scale) one. The scale invariance of the 3D interaction law in the thermo-mechanical coupled particle model is examined. The proposed 3D upscale theory is tested through two typical examples. Finally, a practical application example of 3D transient heat flow in a solid with constant heat flux is given to illustrate the performance of the proposed 3D upscale theory in the thermo-mechanical coupling particle simulation of 3D large-scale non-isothermal problems. Both the benchmark tests and application example are provided to demonstrate the correctness and usefulness of the proposed 3D upscale theory for simulating 3D non-isothermal problems using the particle simulation method. Originality/value The paper provides some important theoretical guidance to modeling 3D large-scale non-isothermal problems at both the engineering length-scale (i.e. the meter-scale) and the geological length-scale (i.e. the kilometer-scale) using the particle simulation method directly.


2014 ◽  
Vol 1008-1009 ◽  
pp. 850-860 ◽  
Author(s):  
Zhou Wei Zhang ◽  
Jia Xing Xue ◽  
Ya Hong Wang

A calculation method for counter-current type coil-wound heat exchanger is presented for heat exchange process. The numerical simulation method is applied to determine the basic physical parameters of wound bundles. By controlling the inlet fluid velocity varying in coil-wound heat exchanger to program and calculate the iterative process. The calculation data is analyzed by comparison of numerical result and the unit three dimensional pipe bundle model was built. Studies show that the introduction of numerical simulation can simplify the pipe winding process and accelerate the calculation and design of overall configuration in coil-wound heat exchanger. This method can be applied to the physical modeling and heat transfer calculation of pipe bundles in coil wound heat exchanger, program to calculate the complex heat transfer changing with velocity and other parameters, and optimize the overall design and calculation of spiral bundles.


Author(s):  
H. X. Liang ◽  
Q. W. Wang ◽  
L. Q. Luo ◽  
Z. P. Feng

Three-dimensional numerical simulation was conducted to investigate the flow field and heat transfer performance of the Cross-Wavy Primary Surface (CWPS) recuperators for microturbines. Using high-effective compact recuperators to achieve high thermal efficiency is one of the key techniques in the development of microturbine in recent years. Recuperators need to have minimum volume and weight, high reliability and durability. Most important of all, they need to have high thermal-effectiveness and low pressure-losses so that the gas turbine system can achieve high thermal performances. These requirements have attracted some research efforts in designing and implementing low-cost and compact recuperators for gas turbine engines recently. One of the promising techniques to achieve this goal is the so-called primary surface channels with small hydraulic dimensions. In this paper, we conducted a three-dimensional numerical study of flow and heat transfer for the Cross-Wavy Primary Surface (CWPS) channels with two different geometries. In the CWPS configurations the secondary flow is created by means of curved and interrupted surfaces, which may disturb the thermal boundary layers and thus improve the thermal performances of the channels. To facilitate comparison, we chose the identical hydraulic diameters for the above four CWPS channels. Since our experiments on real recuperators showed that the Reynolds number ranges from 150 to 500 under the operating conditions, we implemented all the simulations under laminar flow situations. By analyzing the correlations of Nusselt numbers and friction factors vs. Reynolds numbers of the four CWPS channels, we found that the CWPS channels have superior and comprehensive thermal performance with high compactness, i.e., high heat transfer area to volume ratio, indicating excellent commercialized application in the compact recuperators.


Author(s):  
S. S. N. Ahmad ◽  
C. Cole ◽  
M. Spiryagin ◽  
Y. Q. Sun

Implementation of a new bogie concept is an integrated part of the vehicle design which must follow a rigorous testing and validation procedure. Use of multibody simulation helps to reduce the amount of time and effort required in selecting a new concept design by analysing results of simulated dynamic behaviour of the proposed design. However, the multibody simulation software mainly looks at the dynamics of a single vehicle; hence, forces from the train configuration operational dynamics are often absent in such simulations. Effects of longitudinal-lateral and longitudinal-vertical interactions between rail vehicles have been found to affect the stability of long trains [1,2]. The effect of wedge design on the vertical dynamics of a bogie has also been discussed in [3,4]. It is important to apply the lateral and vertical forces from a train simulation into a single multibody model of a wagon to check its behaviour when operating in train configuration. In this paper, a novel methodology for the investigation of new bogie designs has been proposed based on integrating dynamic train simulation and the multibody vehicle modelling concept that will help to efficiently achieve the most suitable design of the bogie. The proposed methodology suggests that simulation of any configuration of bogie needs to be carried out in three stages. As the first stage, the bogie designs along with the wagon configurations need to be presented as a multibody model in multibody simulation software to test the suitability of the concept. The model checking needs to be carried out in accordance with the wagon model acceptance procedure established in [5]. As the second stage, the wagon designs need to be tested in train configurations using a longitudinal train dynamics simulation software such as ‘CRE-LTS’ [2], where a train set consisting of the locomotives and wagons will be simulated to give operational wagon parameters such as lateral and vertical coupler force components. As the third stage, the detailed dynamic analysis of bogies and wagons needs to be performed with a multibody software such as ‘Gensys’ where lateral and vertical coupler force components from the train simulation (second stage) will be applied on the multibody model to replicate the worst case scenario. The proposed methodology enhances the selection procedure of any alternate bogie concept by the application of simulated train and vehicle dynamics. The simulated case studies show that simulation of wagon dynamic behaviour in multibody software combined with data obtained from longitudinal train simulation is not only possible, but it can identify issues with a bogie design that can otherwise be overlooked.


Sign in / Sign up

Export Citation Format

Share Document