Use of high pressures for solving problems of hypersonic aerodynamics

The use of high pressures in a hypersonic aerodynamic experiment is founded on physical grounds. Calculation results of Mach and Reynolds numbers reachable at the line of gas condensation are given as functions of the temperature and pressure in the plenum chamber. Approaches to solving problems of designing ultrahigh-pressure facilities that ensure outflow with pressures up to20,000 atm are described. These problems include the stop of the first-stage piston at the point of maximum pressure, suppression of the reaction force, provision of normal operation of seals of the moving piston, and reduction of friction forces in the seals. The principles considered are used in an actually operating facility. Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 41, No. 5, pp. 103–114, September–October, 2000.     

具有封严蓖齿转子系统的动力稳定性分析

封严蓖齿的气弹力对转子系统的稳定性起重要作用，本文建立了具有封严蓖齿的转子系统动力稳定性的分析方法，该方法利用传递矩阵法建立系统的主导方程，将封严蓖齿的气弹力以等效刚度及等效阻尼系数的方式引入转子系统的主导方程，用数值积分法解主导方程，求得系统的动力特性，根据扰动解判断系统的稳定性。本文所建立的分析软件，也可计算转子系统的不平衡响应、临界转速、突加不平衡响应等。文中对一模型转子进行了算例分析。     

Fast evaluation of friction forces in traction-drive contacts including thermal effects

Summary  In automotive traction drives, power is transmitted by friction forces. The friction forces result from the shear stresses developed in lubricated and highly loaded contacts between rolling bodies. Due to the kinematics of a traction drive, shear velocities occur in both the rolling direction and perpendicular to it. Due to these shear velocities and by normal pressure, the lubricant is forced to build up shear stresses. The increase of the shear stresses may be modelled by a nonlinear viscous element. The describing differential equations are coupled by the equivalent shear stress, which defines the nonlinear behaviour of the element. A fast method is described to evaluate the coupled differential equations. By using a known analytical approximation for the equivalent shear stress, the differential equations are decoupled and can be solved analytically. In an iterative procedure the equivalent shear stress is updated, and the complete solution is found. The iterative method is extended to account for thermal effects in the contact. Received 17 June 1999; accepted for publication 26 October 1999     

Non-steady-state elastic–plastic behavior of metal forming (rolling) with a generalized friction law

An analytical and generalized friction law is formulated and a three-dimensional large-strain non-steady-state elastic–plastic finite element analysis has been performed for rolling process. The contact/friction problem at the interface between the workpiece and the rolls is treated rigorously by using this new friction law. The numerical results, including the evolutions of roll torque and roll forces, and the interfacial normal and shear stresses, are presented and discussed.     

Prediction of draft forces in cohesionless soil with the Discrete Element Method

The Discrete Element Method (DEM) is applied to predict draft forces of a simple implement in cohesionless granular material. Results are compared with small-scale laboratory tests in which the horizontal force is measured at a straight blade. This study is focused on the case of cohesionless material under quasi-static conditions.The DEM requires the calibration of the local contact parameters between particles to adjust the bulk material properties. The most important bulk property is the angle of internal friction ?. In the DEM, the shear resistance is limited in the case of spherical particles due to excessive particle rotations. This is cured by retaining rotations of the particles. Although this is known to prevent the material from developing shear bands, the model still turns out to be capable of predicting the reaction force on the blade.In contrast to empirical formulas for this kind of application, the DEM model can easily be extended to more complex tool geometries and trajectories. This study helps to find a simple and numerically efficient setup for the numerical model, capable of predicting draft forces correctly and so allowing for large-scale industrial simulations.     

Mechanics of advancing pin-loaded contacts with friction

This paper considers finite friction contact problems involving an elastic pin and an infinite elastic plate with a circular hole. Using a suitable class of Green's functions, the singular integral equations governing a very general class of conforming contact problems are formulated. In particular, remote plate stresses, pin loads, moments and distributed loading of the pin by conservative body forces are considered. Numerical solutions are presented for different partial slip load cases. In monotonic loading, the dependence of the tractions on the coefficient of friction is strongest when the contact is highly conforming. For less conforming contacts, the tractions are insensitive to an increase in the value of the friction coefficient above a certain threshold. The contact size and peak pressure in monotonic loading are only weakly dependent on the pin load distribution, with center loads leading to slightly higher peak pressure and lower peak shear than distributed loads. In contrast to half-plane cylinder fretting contacts, fretting behavior is quite different depending on whether or not the pin is allowed to rotate freely. If pin rotation is disallowed, the fretting tractions resemble half-plane fretting tractions in the weakly conforming regime but the contact resists sliding in the strongly conforming regime. If pin rotation is allowed, the shear traction behavior resembles planar rolling contacts in that one slip zone is dominant and the peak shear occurs at its edge. In this case, the effects of material dissimilarity in the strongly conforming regime are only secondary and the contact never goes into sliding. Fretting tractions in the forward and reversed load states show shape asymmetry, which persists with continued load cycling. Finally, the governing integro-differential equation for full sliding is derived; in the limiting case of no friction, the same equation governs contacts with center loading and uniform body force loading, resulting in identical pressures when their resultants are equal.     

Strain and work-hardening in a coulomb-type granular medium

Equations are given for the elastoplastic strain of a granular medium together with experimental relationships for sands. Equations for the characteristics are drawn up in general form for the two-dimensional case. It is shown that experimental data on the dilatancy rate as a function of the angle of internal friction reflect the condition for orthogonality of the characteristics of the velocity distribution to the direction of the dry-friction forces in the sliding areas. The flow in a shear tester is discussed. The calculation from the universal relationships agrees with experiment.     

Response of an elastic half-space with power-law nonhomogeneity to static loads

In this paper a series of problems for an isotropic elastic half-space with power-law nonhomogeneity are considered. The action of surface vertical and horizontal forces applied to the half-space is studied. A part of the paper deals with the case of zero-valued surface shear modulus (for positive values of the power determining the nonhomogeneity). This condition leads to simple solutions for two-dimensional (2D) case when radial distribution of stresses exists for surface loads concentrated along an infinite line. Corresponding results for the three-dimensional (3D) case are constructed on the basis of the relationships between 2D and 3D solutions developed in the paper. A more complicated case, in which the shear modulus at the surface of the half-space differs from zero, is treated using fundamental solutions of the differential equations for Fourier–Bessel transformations of displacements. In the paper the fundamental solutions are built in the following two forms: (a) a combination of functions expressing displacements of the half-space under the action of vertical and horizontal forces in the case of zero surface shear modulus, and (b) a representation of the fundamental solutions using confluent hypergeometric functions. The results of numerical calculation given in the paper relate to Green functions for the surface vertical and horizontal point forces.     

三元乙丙密封材料不同工况下的摩擦性能

采用SRV摩擦磨损试验机对用于核主泵副密封的三元乙丙橡胶材料与具有DLC涂层的钢配偶件进行摩擦性能试验,获得该橡胶密封材料的摩擦系数随润滑状态、载荷、温度和相对滑动速度等的变化规律,并通过摩擦表面的形貌分析揭示了其摩擦机理,可为核主泵机械密封性能分析及副密封选材提供依据.     

Simulating shear behavior of a sandy soil under different soil conditions

Understanding of soil shear behavior is very important in the field of agricultural machinery and soil dynamics. In this study, a discrete element model was developed using a simulation tool, Particle Flow Code in Three Dimensions (PFC^3D). The model simulates direct shear tests of soil and predicts soil shear behavior, in terms of shear forces and displacements. To determine and calibrate model parameters (stiffness of particles, strength and stiffness of bond between particles), laboratory direct shear tests were conducted to examine effects of soil moisture content and bulk density on shear behaviors of a sandy soil. Three soil moisture levels (0.02%, 13.0%, and 21.5%) and four bulk density levels (0.99, 1.28, 1.36, and 1.50 Mg/m³) were used in the tests. The test results showed that in general drier and denser soil conditions produced higher shear forces. Based on the test results, the bond strengths of the model particles were determined from soil cohesion and internal friction angle. The model particle stiffness was calibrated based on the yield forces from the tests. The calibrated particle stiffness varied from 1.0 × 10³ to 8.2 × 10³ N/m, depending on soil moisture and density levels. The bond stiffness calibrated was 1.0 × 10⁷ Pa/m for all soil conditions.     

Calibration of granular material parameters for DEM modelling and numerical verification by blade-granular material interaction

The discrete element method (DEM) is a promising approach to model blade-granular material interactions. The accuracy of DEM models depends on the model parameters. In this study, a calibration process was developed to determine the parameter values. The particle size was the same as the real material and the particle shape was modelled using two spherical particles rigidly clumped together to form a single grain. Laboratory shear tests and compressions tests were used to determine the material internal friction angle and stiffness, respectively. These tests were replicated numerically using DEM models with different sets of particle friction coefficients and particle stiffness values. The shear test results are found to be dependent on both the particle friction coefficient and the particle stiffness. The compression test results show that it is only dependent on the particle stiffness. The combination of shear test and compression test results can be used to determine a unique set of particle friction and particle stiffness values. The calibration process was validated experimentally and numerically by modelling a blade moving through granular material. Results show that the forces acting on the blade can be accurately modelled with DEM and the maximum error is found to be 26%. The relative particle-blade displacements were used to predict the position and shape of the shear lines in front of the blade. A good qualitative correlation was achieved between the experiments and the DEM simulations.     

颤振环境MoS2/Ag薄膜碰撞滑动接触摩擦性能研究

基于空间机构的运动特性,考虑空间颤振环境的影响,采用粗粒化分子动力学研究MoS₂/Ag薄膜的碰撞滑动接触摩擦性能,建立颤振环境碰撞滑动接触摩擦的粗粒化分子动力学模型,对比了纯Ag和MoS₂/Ag薄膜的摩擦性能,研究了初始碰撞速度、滑动速度以及空间温度对碰撞滑动接触摩擦过程的影响.结果表明：与纯Ag相比,MoS₂/Ag薄膜表现出更优异的摩擦性能;压头碰撞速度对动能有一定的贡献,初始碰撞速度的增加会增大压头压入基体的深度,使得平均摩擦力增大;滑动速度的增加会加剧原子间的相互剪切摩擦,使平均摩擦力增加;MoS₂/Ag薄膜在100～500 K温度范围内表现出良好的摩擦性能,当空间温度为600 K时,其摩擦性能降低,并伴随着MoS₂膜的破裂.     

干气密封滑动摩擦界面切向接触刚度分形模型

为揭示干气密封滑动摩擦界面高频微幅自激摩擦振动规律,用分形参数表征摩擦界面形貌特性,根据重新建立的微凸体接触变形方式,以及对非协调弹性体在切向力作用下初始滑动问题的研究,建立了干气密封滑动摩擦界面切向接触刚度分形模型. 通过数值对切向接触刚度的影响因素进行了分析,研究结果表明:切向接触刚度随分形维数、真实接触面积和材料特性系数的增大而增大;切向接触刚度随特征尺度、摩擦系数的增大逐渐减小. 相比于分形维数、特征尺度和材料特性系数对切向接触刚度的影响,摩擦系数的影响相对较小. 这些研究结果为进一步研究干气密封高频微幅自激摩擦振动奠定了基础.      

Contact problem for a plane crack under a normally incident antiplane shear wave

The contact-interaction problem for a stationary plane crack with friction between its edges under the action of a normal (to the crack plane) harmonic shear wave is addressed. Antiplane deformation conditions are considered. The distribution of contact forces and displacement discontinuity of crack edges are studied Published in Prikladnaya Mekhanika, Vol. 43, No. 5, pp. 138–142, May 2007.     

子午线轮胎接触变形的结构有限元分析

针对195／60R14子午线轮胎建立了三维非线性有限元模型，着重研究了额定充气压力及静载荷作用下帘线承受拉应力和剪应力的基本特征。计算结果表明，接地区域摩擦力呈斜对称分布，反映了轮胎中帘线-橡胶复合材料存在变形耦合效应；冠带层、带束层、胎体帘线应力分布较为复杂，载荷变化对其应力水平和分布影响较大；在胎肩部位应力较高，且随载荷变化局部帘线应力变化剧烈，在承受交变载荷时，易形成层间剥落。分析结果有助于预测轮胎的使用性能，可以针对性地应用于因轮胎结构设计引起的质量损坏，某些对轮胎使用性能不利的受力状态可通过结构的优化设计来克服。     

Brownian dynamics simulation of a bead-spring chain model with configuration-dependent anisotropic mobility

The paper presents numerical studies of a bead-spring chain model for concentrated polymer solutions and melts. Anisotropic configuration-dependent mobility tensors determine the interaction forces between the molecular chains. The analysis consists of the solution of the Langevin equation in which the Brownian forces are numerically simulated. For a chain of 10 beads connected by Hookean springs the rheological properties are studied for the case of a sudden onset of a steady shear. The results flow illustrate the effects of the anisotropy of the forces on the beads and of different shear rates. The experimentally observed overshoot in the shear growth function and the first normal stress function is simulated by the model, albeit not quantitavely.     

Two-dimensional sliding frictional contact of functionally graded materials

A multi-layered model for sliding frictional contact analysis of functionally graded materials (FGMs) with arbitrarily varying shear modulus under plane strain-state deformation has been developed. Based on the fact that an arbitrary curve can be approached by a series of continuous but piecewise linear curves, the FGM is divided into several sub-layers and in each sub-layers the shear modulus is assumed to be a linear function while the Poisson's ratio is assumed to be a constant. In the contact area, it is assumed that the friction is one of Coulomb type. With this model the fundamental solutions for concentrated forces acting perpendicular and parallel to the FGMs layer surface are obtained. Then the sliding frictional contact problem of a functionally graded coated half-space is investigated. The transfer matrix method and Fourier integral transform technique are employed to cast the problem to a Cauchy singular integral equation. The contact stresses and contact area are calculated for various moving stamps by solving the equations numerically. The results show that appropriate gradual variation of the shear modulus can significantly alter the stresses in the contact zone.     

结构面抗剪强度参数对岩质边坡稳定的影响

边坡工程中,结构面抗剪强度是非常重要的力学参数,合理地选择确定抗剪强度的方法,常常要分析c、φ值在边坡稳定中的作用。根据极限平衡原理,对10~30m高的平面滑动型边坡不同抗剪强度的控稳结构面c、φ值在边坡稳定中的作用进行了定量分析。研究表明c值对稳定系数的影响随坡高增大逐渐减弱,φ值对稳定系数的影响随坡高增大逐渐增强。得到了不同坡高下c、φ值边坡稳定作用等值曲线,该曲线很容易确定抗剪强度参数c、φ值在边坡稳定中的作用。     

Non-linear dynamic response of a rotating radial Timoshenko beam with periodic pulse loading at the free-end

Consideration is given to the dynamic response of a Timoshenko beam under repeated pulse loading. Starting with the basic dynamical equations for a rotating radial cantilever Timoshenko beam clamped at the hub in a centrifugal force field, a system of equations are derived for coupled axial and lateral motions which includes the transverse shear and rotary inertia effects, as well. The hyperbolic wave equation governing the axial motion is coupled with the flexural wave equation governing the lateral motion of the beam through the velocity-dependent skew-symmetric Coriolis force terms. In the analytical formulation, Rayleigh-Ritz method with a set of sinusoidal displacement shape functions is used to determine stiffness, mass and gyroscopic matrices of the system. The tip of the rotating beam is subjected to a periodic pulse load due to local rubbing against the outer case introducing Coulomb friction in the system. Transient response of the beam with the tip deforming due to rub is discussed in terms of the frequency shift and non-linear dynamic response of the rotating beam. Numerical results are presented for this vibro-impact problem of hard rub with varying coefficients of friction and the contact-load time. The effects of beam tip rub forces transmitted through the system are considered to analyze the conditions for dynamic stability of a rotating blade with intermittent rub.     

激光加工微凹坑轴表面对唇形密封泵吸及摩擦特性的影响

为了研究轴表面形貌对径向唇形密封泵吸作用及摩擦特性的影响,利用Nd:YAG固体脉冲激光对旋转轴表面进行微凹坑织构化处理,在油封密封试验机上,针对不同转速、不同表面形貌等工况,完成了一系列对比试验,获取了试验过程中的泵吸量、摩擦力以及密封圈温升等参数.并通过扫描电镜对试验后各密封唇表面进行了分析.结果表明:三角形微凹坑顶点朝向油液侧的表面能够增强泵吸作用,且转速对泵吸率影响很大,而其他形状和参数的微凹坑表面不利于泵吸作用的产生;与未织构面相比,圆形凹坑面积占有率为7%和21%的表面具有良好的减摩效果,三角形凹坑顶点沿周向的两表面间减摩效果较好.