Using Hob Offset to Balance Dynamic Strength in Spur Gears*

An analytical study of the effect of hob offset on dynamic tooth strength of spur gears is presented. The study was limited to equal and opposite offset values applied to the pinion and gear to maintain the standard operating center distance. The analysis presented is performed using a new version of the NASA gear dynamics code DANST.

The operating speed of a transmission has a significant influence on the amount of hob offset required to equalize dynamic stresses in the pinion and gear. In the transmission studied, at low speeds, the optimum hob offset was found to fluctuate within a range. At higher speeds, the optimum value is constrained by the minimum allowed thickness at the tip of the pinion tooth. For gears that must operate over a range of speeds, an average offset value can be used. Spur gears designed with the procedure presented here can have significant improvements in load capacity.     

Mixed convection heat transfer from a horizontal channel with protruding heat sources

A numerical investigation is carried out to study fluid flow and heat transfer characteristics of conjugate mixed convection from a two dimensional horizontal channel with four protruding heat sources mounted on one of the finite thick channel walls. The flow is assumed as laminar, hydrodynamically and thermally developing. Water and FC70 are the fluids under consideration. The geometric parameters such as spacing between the channel walls (S), size of protruding heat sources (L_h×t_h), thickness of substrate (t) and spacing between heat sources (b) are fixed. Results are presented to show the effect of parameters such as Re_S, Gr_S^*, Pr, k_p/k_f and k_s/k_f on fluid flow and heat transfer characteristics. Using the method of asymptotic expansions, correlations are also presented for the maximum temperature of heat source.     

齿向修形对滤波减速器润滑性能的影响分析

综合考虑了滤波减速器齿向修形参数、真实齿面粗糙度和瞬态效应等因素,建立了轮齿混合润滑数学模型,数值计算了不同修形参数值对应不同啮合点的最大压力和中心膜厚,分析了齿面粗糙度和转速对润滑性能的影响.结果表明:修形参数r和Ry均存在一个优化范围,使得轮齿表面最大油膜压力显著降低,边缘效应弱化,而中心膜厚则随着r和Ry的增大而逐渐增大;未修形轮齿边缘油膜压力受粗糙度的影响而急剧增大,边缘效应更加显著,修形后轮齿的边缘效应得到了明显改善,因此,轮齿修形也因粗糙表面的存在而显得更加重要;随着转速逐渐降低,轮齿表面的平均油膜厚度逐渐变小,接触比逐渐增大,轮齿表面由弹流润滑逐渐转为混合润滑,最后演变为边界润滑.     

渐开线斜齿轮非稳态弹流润滑数值模拟研究

建立了渐开线斜齿轮啮合的弹流润滑计算模型,将斜齿圆柱齿轮啮合的齿面接触等效为有限长线接触的弹流润滑问题.考虑斜齿轮啮合的实际因素,将斜齿轮啮合过程中的等效曲率半径和齿面载荷的变化反映到弹流润滑计算模型中,应用统一Reynolds方程方法求得轮齿在1个完整啮合周期内的瞬时弹流润滑数值解.结果表明：斜齿轮啮合线上各点处的膜厚、压力均有较大不同,各接触点处的油膜厚度受综合曲率半径的影响较大;斜齿轮传动非稳态效应相对较弱;小齿轮齿根附近和节点位置处润滑状态较差;适当增大压力角可以改善齿轮的润滑.     

Computational approach to design face-milled spiral bevel gear drives with favorable conditions of meshing and contact

The micro-geometry of the tooth surfaces of spiral bevel and hypoid pinions has to be fine-adjusted to obtain enhanced meshing and contact characteristics during the meshing process with their corresponding mating gears. In this paper, a new methodology is proposed to design face-milled spiral bevel gear drives to, firstly, derive favorable orientation and dimensions of the contact pattern between the mating surfaces of the gears and, secondly, obtain a predesigned parabolic function of negative transmission errors with limited magnitude of maximum transmission errors. The proposed approach is based on the definition of the desired topography for the active surfaces of the pinion followed by a numerical derivation of their finishing machine-tool settings through a bound-constrained optimization algorithm. Increasing mechanical strength and reducing the levels of noise and vibration of face-milled spiral bevel gear drives constitute the main objectives of the proposed design process. A numerical example is provided to illustrate the applicability of the developed theory .     

Sensitivity Study of Simulation Parameters Controlling CO<Subscript>2</Subscript> Trapping Mechanisms in Saline Formations

The primary purpose of this study is to understand quantitative characteristics of mobile, residual, and dissolved CO₂ trapping mechanisms within ranges of systematic variations in different geologic and hydrologic parameters. For this purpose, we conducted an extensive suite of numerical simulations to evaluate the sensitivities included in these parameters. We generated two-dimensional numerical models representing subsurface porous media with various permutations of vertical and horizontal permeability (k _v and k _h), porosity (f{\phi}), maximum residual CO₂ saturation (S_gr^max{S_{\rm gr}^{\max}}), and brine density (ρ _br). Simulation results indicate that residual CO₂ trapping increases proportionally to k_v, k_h, S_gr^max{k_{\rm v}, k_{\rm h}, S_{\rm gr}^{\max}} and ρ _br but is inversely proportional to f.{\phi.} In addition, the amount of dissolution-trapped CO₂ increases with k _v and k _h, but does not vary with f{\phi } , and decreases with S_gr^max{S_{\rm gr}^{\max}} and ρ _br. Additionally, the distance of buoyancy-driven CO₂ migration increases proportionally to k _v and ρ _br only and is inversely proportional to k_h, f{k_{\rm h}, \phi } , and S_gr^max{S_{\rm gr}^{\max}} . These complex behaviors occur because the chosen sensitivity parameters perturb the distances of vertical and horizontal CO₂ plume migration, pore volume size, and fraction of trapped CO₂ in both pores and formation fluids. Finally, in an effort to characterize complex relationships among residual CO₂ trapping and buoyancy-driven CO₂ migration, we quantified three characteristic zones. Zone I, expressing the variations of S_gr^max{S_{\rm gr}^{\max}} and k _h, represents the optimized conditions for geologic CO₂ sequestration. Zone II, showing the variation of f{\phi} , would be preferred for secure CO₂ sequestration since CO₂ has less potential to escape from the target formation. In zone III, both residual CO₂ trapping and buoyancy-driven migration distance increase with k _v and ρ _br.     

Enhanced model of load distribution along the line of contact for non-standard involute external gears

The presence of undercut at the tooth root, non-equal addendum on pinion and wheel, non-standard tooth height or non-standard center distance may have decisive influence on the load distribution along the line of contact of spur and helical gear teeth. The curve of variation of the meshing stiffness along the path of contact, quite symmetric respect the midpoint of the interval of contact, loses its symmetry for non-standard geometries and operating conditions. As a consequence, the critical contact points for bending and wear calculations may be shifted from their locations for standard gears. In this paper, a non-uniform model of load distribution along the line of contact of standard spur and helical gears, obtained from the minimum elastic potential criterion, has been enhanced to fit with the meshing conditions of the above mentioned non-standard cylindrical gear pairs. The same analytical formulation of the initial model may be used for the non-standard gears by considering appropriate values of a virtual contact ratio, which are also presented in the paper.     

Compensation of errors of alignment and contact pattern repositioning in hypoid gears with low crossing shaft angle

The investigation of the influence of different errors of alignment on the mesh behavior of hypoid gears with low crossing shaft angle has been carried out. As a result, a simplified methodology of compensation of errors of alignment and correction of the contact pattern by controlled axial displacements to the pinion and wheel of the hypoid gear set has been proposed. The shaft angle error and the offset error mainly influence the position of the contact pattern in profile direction. The pinion and wheel axial position errors influence the position of contact pattern along the tooth trace direction and change the position slightly along the tooth profile direction. All the assembly errors increase the peak-to-peak value of transmission errors, except positive values of the pinion axial position error that allows decreasing the peak-to-peak transmission error slightly in an appropriate range. Finite element analysis has been performed to verify the compensating ability of the developed methodology.

     

Three-dimensional photoelastic study of stresses in rack gears

A three-dimensional photoelastic analysis using the stress-freezing and slicing techniques was employed in studying the stress distribution and concentration across the thickness of a thick rack-gear tooth. Full-size photoelastic models of about 25-mm thickness, having various teeth parameters were machined from PLM-4B plates. A special mounting fixture was designed to hold the models in the loading frame and a knife-edged line load was used in loading the models. In the analysis, the method of oblique incidence was adopted for the separation of the principal stresses. The intension was to determine the magnitude and location of the maximum stress at the tensile fillet and to establish the stress-concentration factors for various geometric tooth configurations.From this study, it was concluded that the fillet stresses depend largely upon the pressure angle, fillet radius and the position of the load. Maximum values of the tensile fillet stresses occur at the middle plane of the tooth thickness. Furthermore, the value of the stress-concentration factor increases with the increase of the tooth thickness. In general, this investigation gave values of stresses much higher than those values calculated by the simple-flexure theory.     

Comparison of energy release rate and minimum strain energy density for potential failure of composite with bi-material interface

An analytical method is developed to describe the fields of stress and displacement in a bi-material strip specimen with an edge interfacial crack. All of the basic governing equations, boundary conditions on crack surfaces and conditions of continuity along the interface are satisfied by the eigenfunction expansion method. The other boundary conditions are satisfied by the generalized variational principle. The stress intensity factors are calculated for determining the energy release rate and minimum strain energy density factor S_min that is used the strain energy density criterion. Problems with oscillatory singularity and contact zone are discussed. Not only the effects of bi-material modulus ratio, thickness ratio, Poisson's ratio and crack length to S_min, but also the influences of bi-material modulus ratio, thickness ratio to phase angle are presented. Among these parameters, particular situations where S_min become jeopardously high and lead to failure are discussed.     

Estimation of three-phase relative permeabilities

Starting from the statistical structural model of Alemánet al. (1988), we have developed an alternative to Stone's (1970, 1973; Aziz and Settari, 1979) methods for estimating steady-state, three-phase relative permeabilities from two sets of steady-state, two-phase relative permeabilities. Our result reduces to Stone's (1970; Aziz and Settari, 1979) first method, when the steady-state, two-phase relative permeability of the intermediate-wetting phase with respect to either the wetting phase or the nonwetting phase is a linear function of the saturation of the intermediate-wetting phase. As the curvature of either of these relative permeability functions increases, the deviation of our result from Stone's (1970; Aziz and Settari, 1979) first method increases. Currently, there are no data available that are sufficiently complete to form the basis of a comparison between our result and either of the methods of Stone (1970, 1973; Aziz and Settari, 1979).Notation a free parameter in Equation (19) - B(m, n) Beta function defined by Equation (17) - F ^(w), F^(nw) defined by Equations (31) and (27), respectively - G ⁽ⁱ⁾ defined by Equations (37) and (39) - H ⁽ⁱ⁾ defined by Equations (38) and (40) - k ⁽ⁱ⁾ three-phase relative permeability fo phasei - k ^(i)* defined by Equations (34) through (36) - k ^(i,j) relative permeability to phasei during a two-phase flow with phasej, possibly in the presence of an immobile phase - k ^(i,j)* defined by analogy with Equations (41) and (42) - k ^(i,j)** defined by Equations (49), (50), (53), and (54) - k _max ⁽ⁱ⁾ defined by Equation (11) - k ₁₉₇₀ ^(iw) defined by Equation (10) - k ₁₉₇₃ ^(iw)* defined by Equation (58) - k ₁₉₇₃ ^(iw) defined by Equation (13) - L length and diameter of cylindrical averaging surfaceS - L _t length of an individual capillary tube enclosed byS - L _t ^* defined by Equation (19) - L _t,min length of pore whose radius isR _max - N total number of pores contained within the averaging surfaceS - p ₁ ⁽ⁱ⁾ ,p ₂ ⁽ⁱ⁾ pressure of phasei at entrance and exit of averaging surfaceS, respectively - p defined by Equation (21) - p _c ^(i,j) capillary pressure function - p _c ^(i,j)* defined by Equations (23), (29), and (32) - p ⁽ⁱ⁾ intrinsic average of pressure within phasei defined by Alemánet al. (1988) - R pore radius - R ^* defined by Equation (18) - R _max maximum pore radius that occurs withinS - s ⁽ⁱ⁾ local saturation of phasei - s ^(i)* defined by Equation (7) - s _min ⁽ⁱ⁾ minimum or immobile saturation of phasei - S averaging surface introduced in local volume averaging - V ⁽ⁱ⁾ volume of phasei occupying the pore space enclosed byS Greek Letters , parameters in the Beta distribution defined by Equation (16) - ^(w), ^(nw) functions of only the wetting phase saturation and the non-wetting phase saturation, respectively. Introduced in Equation (6) - ^(i,j) interfacial tension between phasesi andj - (x) Gamma function - defined by Equation (57) - , spherical coordinates in system centered upon the axis of the averaging surfaceS - _max maximum value of , 45 °, in view of assumption (9) - ^(i,j) contact angle between phasesi andj measured through the displacing phase - ^(w),^(nw) functions of only the wetting phase saturation and the non-wetting phase saturation, respectively. Introduced in Equation (12) Other gradient operator Amoco Production Company, PO Box 591 Tulsa, OK 74102, U.S.A.     

Stabilization using fractional-order PI and PID controllers

This paper presents a solution to the problem of stabilizing a given fractional dynamic system using fractional-order PIλ and PI^λD^μ controllers. It is based on plotting the global stability region in the (k _p, k _i)-plane for the PI^λ controller and in the (k _p , k _i , k _d)-space for the PIλDμ controller. Analytical expressions are derived for the purpose of describing the stability domain boundaries which are described by real root boundary, infinite root boundary and complex root boundary. Thus, the complete set of stabilizing parameters of the fractional-order controller is obtained. The algorithm has a simple and reliable result which is illustrated by several examples, and hence is practically useful in the analysis and design of fractional-order control systems.     

Errata: Analytical Treatment of Some Extended Problems in Structural Optimization

ABSTRACT

The recursive Newton-Euler dynamic formulation is used to derive the equations of motion of a manipulator with harmonic drives. The derivation is general, in that the harmonic drive is viewed as a separate body that forms a kinematically closed chain with two contiguous links, leading to a comprehensive dynamic model. The harmonic drive is modeled as a flexible and rigid gear with a high gear reduction ratio. Under different modeling assumptions, the effects of gear flexibility and dynamic coupling are examined using a seven-degree-of-freedom manipulator.     

Variation of crack-opening stresses in three-dimensions: finite thickness plate

Crack growth and closure behavior of a center cracked finite thickness plate subjected to constant amplitude cyclic load is investigated by means of a three-dimensional elastic-plastic finite-element analysis. Results are obtained for initial half crack length c_i to half plate thickness t ratios of c_i/t = 3.891 and 1.465 which shall be referred to, respectively, as thin and thick plate. A constant amplitude load with R = S_min/S_max = 0.1 and S_max/σ₀ = 0.25 is applied, where S stands for the stress amplitude and σ₀ the effective yield stress. Crack closure for the thinner plate is found to be largest at and near the free plate surface and to decrease toward the interior during the unloading portion of cyclic loading. The closure pattern stabilizes at the interior and exterior regions, respectively, for c_i/t = 3.981 at 0.34S_max and 0.56S_max and for c_i/t = 1.465 at 0.26S_max and 0.46S_max.A load-reduced displacement technique was used to determine crack-opening stresses at specified locations in the plate from the displacements calculated after 7th cycle (using unloading and reloading portions of cyclic loading). All locations were on the plate exterior surface and were located behind the crack tip and at the centerline of the crack. The opening stresses at the specified points as certain percentage of the maximum stress amplitude were obtained.     

Reflection, transmission and excitation of SH-surface waves by a discontinuity in mass-loading on a semi-infinite elastic medium

The scattering of an SH-wave by a discontinuity in mass-loading on a semi-infinite elastic medium is investigated theoretically. The incident wave is either a plane body wave or a plane SH-surface wave. The problem is reduced to a Wiener-Hopf problem for the scattered wave. In this problem the amplitude spectral density of the particle displacement occurs as unknown function. Special attention is given to the numerical values of the surface wave contributions to the scattered field.Nomenclature x ₁, x ₂, x ₃ Cartesian coordinates - , polar coordinates in x ₁, x ₃-plane - volume mass density - surface mass density of mass-loading - , Lamé constants - U scalar wave function, defined by (2.1) - c _S speed of propagation of uniform shear waves in bulk medium (c _S=(/)^1/2) - angular frequency - t time - k _S wave number of uniform shear waves (k _S=/c _S) - reduced specific acoustic impedance of mass-loading (=k _S /) - k _m wave number of SH-surface wave (k _m=k _S(1+ ²)^1/2) - _1,2,3 partial differentiation with respect to x _1,2,3 - ⁱ angle between x ₃-axis and direction of propagation of incident body wave - ⁱ wave number in horizontal direction of incident body wave ( ⁱ=k _S sin( ⁱ)) - ⁱ wave number in vertical direction of incident body wave ( ⁱ=k _S cos( ⁱ)) - C _1,2 complex amplitude of surface wave excited by a body wave - R reflection factor of surface wave, when surface wave is incident - T transmission factor of surface wave, when surface wave is incident - S particle displacement vector The research presented in this paper has been carried out with partial financial support from the Delfts Hogeschoolfonds.     

The Effect of Wettability on Three-Phase Relative Permeability

We study three-phase flow in water-wet, oil-wet, and fractionally-wet sandpacks. We use CT scanning to measure directly the oil and water relative permeabilites for three-phase gravity drainage. In an analogue experiment, we measure pressure gradients in the gas phase to determine the gas relative permeability. Thus we find all three relative permeabilities as a function of saturation. We find that the gas relative permeability is approximately half as much in a oil-wet medium than in an water-wet medium at the same gas saturation. The water relative permeability in the water-wet medium and the oil relative permeability in the oil-wet medium are similar. In the water-wet medium the oil relative permeability scales as k _roS _o ⁴ for S _o>S _or, where S _or is the waterflood residual oil saturation. With octane as the oil phase, k _roS _o ² for S _o<S _or, while with decane as the oil phase, k _ro falls sharply for S _o<S _or. The water relative permeability in the oil-wet medium resembles the oil relative permeability in the water-wet medium for a non-spreading oil such as decane. These observations can be explained in terms of wetting, spreading, and the pore scale configurations of fluid.     

Soil parameters to predict the performance of off-road vehicles

Soil parameters c, , k_c, k and n were determined by use of a Penetro-Shear apparatus. This device utilizes the rotating and penetrating motion of a circular plate. The performance of a full scale M113 track is predicted using soil parameters calculated from the Penetro-Shear Apparatus data and compared to experimental determination of the drawbar pull-weight ratio vs.per cent of track slip curves.     

Outer layer similarity in fully rough turbulent boundary layers

Turbulent boundary layer measurements were made on a flat plate covered with uniform spheres and also on the same surface with the addition of a finer-scale grit roughness. The measurements were carried out in a closed return water tunnel, over a momentum thickness Reynolds number (Re) range of 3,000–15,000, using a two-component, laser Doppler velocimeter (LDV). The results show that the mean profiles for all the surfaces collapse well in velocity defect form. Using the maximum peak to trough height (R_t) as the roughness length scale (k), the roughness functions (U⁺) for both surfaces collapse, indicating that roughness texture has no effect on U⁺. The Reynolds stresses for the two rough surfaces also show good agreement throughout the entire boundary layer and collapse with smooth wall results outside of the roughness sublayer. Quadrant analysis and the velocity triple products show changes in the rough wall boundary layers that are confined to y<8k_s, where k_s is the equivalent sand roughness height. The present results provide support for Townsends wall similarity hypothesis for uniform three-dimensional roughness. However, departures from wall similarity may be observed for rough surfaces where 5k_s is large compared to the thickness of the inner layer.     

Conjugate natural convection in a square porous cavity

 Steady-state conjugate natural convection in a square cavity filled with a porous medium is studied numerically in this paper. The enclosure consists of two horizontal conductive walls of finite thickness and two vertical walls at different uniform temperatures. The focus is on the role of solid-fluid conductivity ratio, k, on the flow and heat transfer characteristics and the average Nusselt number, , over the vertical hot and cold walls of the cavity for a limited set of particular parameters. It was shown that the interface temperature, θ_w, along the top of the solid wall decreases with the increase in the wall conductivity k. Also, the values of decreases with the increase of the values of the parameter k. Comparison with known results from the open literature when the wall thickness of the horizontal solid walls is neglected (non-conjugate problem) is excellent. Received on 4 April 2000     

Soret effect and slow mass diffusion as a catalyst for overstability in Marangoni-Bénard flows

We study the onset of time dependent Marangoni-Bénard convection in binary mixtures subject to Soret effect by numerical computation of linear instability thresholds in infinite fluid layers and two-dimensional boxes. The calculations are done for positive Marangoni numbers (Ma > 0) and negative Marangoni Soret parameters S _M = –(D _S γ _c )/(Dγ _T ) where D _S and D are the Soret and mass diffusion coefficients, respectively, and γ _T , γ _c are the first derivatives of the surface tension with respect to temperature and concentration. Our purpose is to understand why for particular choices of Prandtl and Schmidt numbers, the increase of the stabilizing solutal contribution leads to a decrease of the critical temperature difference, a phenomenon already reported by Chen & Chen [5] and Skarda et al. [12] For various choices of Prandtl and Schmidt numbers we analyze the evolution of the critical Marangoni number Ma _c , critical wavenumber k _c and angular frequency ω _c with S _M and compute the corresponding eigenvectors. We next propose a physical mechanism which explains how the stabilizing solutal contribution acts as a catalyst for overstability. Finally, we extend our results to two dimensional boxes of small aspect ratio.