Properties of dynamic rupture and energy partition in a solid with a frictional interface

We study properties of dynamic ruptures and the partition of energy between radiation and dissipative mechanisms using two-dimensional in-plane calculations with the finite element method. The model consists of two identical isotropic elastic media separated by an interface governed by rate- and state-dependent friction. Rupture is initiated by gradually overstressing a localized nucleation zone. Different values of parameters controlling the velocity dependence of friction, the strength excess parameter and the length of the nucleation zone, lead to the following four rupture modes: supershear crack-like rupture, subshear crack-like rupture, subshear single pulse and supershear train of pulses. High initial shear stress and weak velocity dependence of friction favor crack-like ruptures, while the opposite conditions favor the pulse mode. The rupture mode can switch from a subshear single pulse to a supershear train of pulses when the width of the nucleation zone increases. The elastic strain energy released over the same propagation distance by the different rupture modes has the following order: supershear crack, subshear crack, supershear train of pulses and subshear single pulse. The same order applies also to the ratio of kinetic energy (radiation) to total change of elastic energy for the different rupture modes. Decreasing the dynamic coefficient of friction increases the fraction of stored energy that is converted to kinetic energy. General considerations and observations suggest that the subshear pulse and supershear crack are, respectively, the most and least common modes of earthquake ruptures.     

Impulsive motion of a suspension: Effect of antisymmetric stresses and particle rotation

A two-phase continuum theory (two-fluid model) for a suspension of rigid spherical particles in a Newtonian fluid is applied to investigate theoretically the flow induced by impulsive motion of an infinite flat plate. Consideration of rotational intertia of the particles gives rise to an antisymmetric part of the volume averaged stress tensor of the continuous phase. The influence of particle rotation and of antisymmetric stresses of the continuous phase, which depend on the relative rotational motion between the particles and the ambient fluid, on the motion of each phase and on the skin friction is examined.Approximate solutions to the equations, corresponding to the physical situation of large and small particle slip, are obtained by power series expansions for small and large times.     

Nonlinear dynamics of an elastic rod with frictional impact

A model is presented for the impact with friction of a flexible body in translation and rotation. This model consists of a system of nonlinear differential equations which considers the multiple collisions as well as frictional effects at the contacting end, and allows one to predict the rigid and elastic body motion after the impact. The kinetic energy is derived by utilizing a generalized velocity field theory for elastic solids. The model uses a dry coefficient of friction and a nonlinear contact force. We introduce a finite number of vibrational modes to take into account the vibrational behavior of the body during impact. The vibrations, the multiple collisions, and the angle of incidence angle, are found to be important factors for the kinematics of frictional impact. Analytical and experimental results were compared to establish the accuracy of the model.     

Nonlinear dynamics of oriented elastic solids

Based on a continuum model for oriented elastic solids the set of nonlinear dispersive equations derived in Part I of this work allows one to investigate the nonlinear wave propagation of the soliton type. The equations govern the coupled rotation-displacement motions in connection with the linear elastic behavior and large-amplitude rotations of the director field. In the one-dimensional version of the equations and for two simple configurations an exhaustive study of solitons is presented. We show that the transverse and/or longitudinal elastic displacements are coupled to the rotational motion so that solitons, jointly in the rotation of the director and the elastic deformations, are exhibited. These solitons are solutions of a system of linear wave equations for the elastic displacements which are nonlinearly coupled to a sine-Gordon equation for the rotational motion. For each configuration, the solutions are numerically illustrated and the energy of the solitions is calculated. Finally, some applications of the continuum model and the related nonlinear dynamics to several physical situations are given and additional more complex problems are also evoked by way of conclusion.     

Wave propagation in a non-ideal dusty gas

In this paper we have studied the behavior of wave motion as propagating wavelets and their culmination into shock waves in a non-ideal gas with dust particles. In the absence of non-ideal effect the gas satisfies an equation of state of Mie–Gruneisen type. An expansion wave resulting from the action of receding piston is considered and the solutions to this problem showing effects of dust particles and non-idealness are obtained. The propagation of weak waves is considered and the flow variables in the region bounded by the piston and the characteristic wave front are found out. The expansive action of a receding piston undergoing an abrupt change in velocity is discussed. Cases of central expansion fan and shock fronts are studied and the solutions up to first order in the physical plane are obtained. The effects of non-idealness and dust particles are discussed in each case.     

基于离散单元法和人工神经网络的近壁颗粒动力学特征研究

颗粒与壁面的相互作用往往对颗粒流动具有显著影响. 为研究颗粒与壁面作用机理, 对滚筒内颗粒流动过程进行离散单元法(DEM)数值模拟. 基于模拟结果统计分析靠近壁面处颗粒的运动特征, 结果表明, 小摩擦系数时颗粒平动和旋转速度均近似满足正态分布, 但由于壁面影响, 摩擦系数增大时颗粒沿滚筒轴向的旋转速度偏离正态分布, 颗粒动力学理论推导壁面边界条件时应考虑速度正态分布的修正及速度脉动的各向异性. 采用人工神经网络(ANN)构建了颗粒无因次旋转温度、滑移速度和平动温度之间的函数模型, 进而可以在常规双流模型壁面边界条件中考虑颗粒旋转的影响. 基于DEM模拟及结果分析可以为壁面边界条件的理论构造和半经验修正提供基础数据和封闭模型.      

Wave-modulated orbits in rate-and-state friction

A frictional spring-block system has been widely used historically as a model to display some of the features of two slabs in sliding frictional contact. Putelat et al. (2008) [7] demonstrated that equations governing the sliding of two slabs could be approximated by spring-block equations, and studied relaxation oscillations for two slabs driven by uniform relative motion at their outer surfaces, employing this approximation. The present work revisits this problem. The equations of motion are first formulated exactly, with full allowance for wave reflections. Since the sliding is restricted to be independent of position on the interface, this leads to a set of differential-difference equations in the time domain. Formal but systematic asymptotic expansions reduce the equations to differential equations. Truncation of the differential system at the lowest non-trivial order reproduces a classical spring-block system, but with a slightly different “equivalent mass” than was obtained in the earlier work. Retention of the next term gives a new system, of higher order, that contains also some explicit effects of wave reflections. The smooth periodic orbits that result from the spring-block system in the regime of instability of steady sliding are “decorated” by an oscillation whose period is related to the travel time of the waves across the slabs. The approximating differential system reproduces this effect with reasonable accuracy when the mean sliding velocity is not too far from the critical velocity for the steady state. The differential system also displays a period-doubling bifurcation as the mean sliding velocity is increased, corresponding to similar behaviour of the exact differential-difference system.     

Heat transfer from rotating finned heat exchangers with different orientation angles

The local and average heat transfer characteristics of spoke like fins that extend outward from a rotating shaft have been determined experimentally. The experiments encompassed a number of geometrical parameters, including the length and chord of the fins, the number of fins deployed around the circumference of the shaft and the orientation angles of the fin. The experiments cover a wider range of rotational speeds, which varies from 25 up to 2,000 rpm. Three wire heat flux sensors have been used in conjunction with a slip ring apparatus to evaluate the local and average heat transfer coefficients. The output results indicated that, the heat transfer transition on rotating fins occurs at Reynolds number lower than encountered on the stationary rectangular fins in crossflow. In general, with non zero incidence angle, the rotating system acts as a fan and creates axial air motion, which enhance the heat transfer rate. However, the effect of orientation angle reduces with increasing the rotational speed. The Nusselt number data are independent of the number of fins in the circumferential array at high rotational speed and are weakly dependent at low Reynolds numbers. To facilitate the use of the results for design, correlations were developed which represent the fin heat transfer coefficient as a continuous function of the investigated independent parameters.     

DYNAMIC ANALYSIS OF ECCENTRIC WHEELS ALONG SLOPE FROM PURE ROLLING TO CONTINUOUS BOUNCING

研究了第十一届全国周培源大学生力学竞赛团体赛中轮子跳高试题的偏心圆轮方案,建立了无初速度释放的偏心圆轮沿斜面纯滚动到连续跳动的动力学理论模型,分析了弹跳现象的产生条件,研究了纯滚动到连续弹跳的动力学过程,给出了斜面倾角、初始释放角和摩擦系数等因素对跳跃行为的影响。     

GAS-KINETIC UNIFIED ALGORITHM FOR BOLTZMANN MODEL EQUATION IN ROTATIONAL NONEQUILIBRIUM AND ITS APPLICATION TO THE WHOLE RANGE FLOW REGIMES

基于过去开展稀薄自由分子流到连续流气体运动论统一算法框架,采用转动惯量描述气体分子自旋运动,确立含转动非平衡效应各流域统一玻尔兹曼模型方程.基于转动能量对分布函数守恒积分,得到计及转动非平衡效应气体分子速度分布函数方程组,使用离散速度坐标法对分布函数方程所依赖速度空间离散降维;应用拓展计算流体力学有限差分方法,构造直接求解分子速度分布函数的气体动理论数值格式;基于物面质量流量通量守恒与能量平衡关系,发展计及转动非平衡气体动理论边界条件数学模型及数值处理方法,提出模拟各流域转动非平衡效应玻尔兹曼模型方程统一算法.通过高、低不同马赫数1：5～25氮气激波结构与自由分子流到连续流全飞行流域不同克努森数（9×10^-4～10）Ramp制动器、圆球、尖双锥飞行器、飞船返回舱外形体再入跨流域绕流模拟研究,将计算结果与有关实验数据、稀薄流DSMC模拟值等结果对比分析,验证统一算法模拟自由分子流到连续流再入过程高超声速绕流问题的可靠性与精度.     

Modeling and theoretical analysis of a novel ratcheting-type cam-based infinitely variable transmission system

An infinitely variable transmission (IVT) is a system that allows for a continuous (non-discrete) variation (including zero) in transmission ratio between two rotating elements. In this paper, a novel ratcheting-type IVT mechanism is presented and its geometrical design and kinematic analysis are studied in details. The proposed system contains two identical units. Each unit includes a cam with a follower, oscillatory slotted links pivoted at a shaft that can be moved vertically by a hydraulic ram (alterable transmission ratio), and a grooved wheel with an actuating rod. The input rotational motion is converted through each unit to an oscillatory angular motion of controlled amplitude. This resulting motion is rectified using a ratchet to get a unidirectional output rotational motion. Therefore, the system output motion will have a different velocity and acceleration than those of the system input. The kinematic analysis revealed that the transmission ratio can be varied continuously in a range from zero to infinity. The analysis also showed that, for particular transmission ratios, the system gives uniform output (angular velocity and acceleration) for a corresponding uniform input.     

Laser-scanning particle image velocimetry applied to a delta wing in transient maneuver

A laser scanning technique, which utilizes a galvanometer scanner to produce particle-image photographs, is employed to investigate the flow over a delta wing undergoing pitching maneuvers at a high angle of attack. Use of a unique forcing system and a large-scale prism arrangement allow characterization of the instantaneous velocity field over the entire crossflow plane at a desired angle of attack.Contours of constant streamwise vorticity are calculated from the crossflow velocity field at various pitching rates. The vorticity distribution occurring during the pitch-up motion differs substantially from that on the stationary wing at the same angle of attack. During the pitch-up motion, the leading-edge vortex is remarkably coherent, in contrast to the disordered structure on the stationary wing. During the corresponding pitch-down motion, the vorticity distribution is quite similar to that on the stationary wing at the same angle of attack. This behavior is evident for a range of pitching rates.     

雨刮器非线性摩擦振动的稳定性与分岔特性

针对雨刮器建立2自由度非线性摩擦振动动力学模型,基于复模态理论计算复特征值并进行稳定性及其对刮刷速度的依赖性分析;通过数值计算分析摩擦振动对刮刷速度的分岔特性,并利用相轨迹、庞加莱映射、频谱特性分析不同刮刷速度下的非线性振动现象.研究发现:摩擦-速度特性的负斜率是导致系统不稳定的根本原因,增大刮刷速度有利于提高系统的稳定性;在高、低刮速区,随着刮刷速度的下降,系统振动形态遵循周期→准周期→混沌的演化规律,并会伴随显著的粘滑振动;仅高速区的周期振动和非振动条件下,刮刷时无附加的粘滑振动.     

Resonant dynamics in a rotordynamic system with nonlinear inertial coupling and shaft anisotropy

The response of a nonlinear, damped Jeffcott rotor with anisotropic stiffness is considered in the presence of an imbalance. For sufficiently small external torque or large imbalance, resonance capture or rotordynamic stall can occur, whereby the rotational velocity of the shaft is unable to increase beyond the fundamental resonance between the rotational and translational motion. This phenomena provides a mechanism for energy transfer from the rotational to the translational mode. Using the method of averaging a reduced-order model is developed, valid near the resonance, that describes this resonant behavior. The equilibrium points of these averaged equations, which correspond to stationary solutions of the original equations and rotordynamic stall, are described as the applied torque, damping, and anisotropy vary. As the anisotropy increases, assumed to arise from increasing shaft cracks, the torque required to eliminate the possibility of stall increases. However, when the system is started with zero initial conditions, the minimum torque required to pass through the resonance is approximately constant as the anisotropy increases. The predictions from the reduced-order model are verified against numerical simulations of the original equations of motion.     

Formation of separation zones in an asymmetric motion of a body in an elastic medium

We obtain an analytical solution of the problem on the motion of a body with wedge-shaped nose in an elastic medium for the case in which a medium separation zone may occur near the nose owing to asymmetry. The character of the dependence of the separation region length on the body velocity, the nose opening angle, the motion asymmetry degree, and the friction coefficient is found. It is shown that if the body moves at a velocity greater than the transverse wave velocity, then there is a limit velocity at which the separation region near the nose of the body disappears.     

活塞二阶运动及活塞裙部摩擦的瞬态分析

为分析不同工况下内燃机活塞-缸套系统的摩擦学和动力学行为,本文从内燃机系统动力学分析入手,结合流体润滑方程,建立了考虑系统变惯量影响的活塞二阶运动与摩擦学行为的瞬态分析模型.以某型号单缸机为例,计算和比较了整个启动过程中曲轴转速、二阶运动及摩擦功耗的变化;对比分析了启动时第一周期与稳定运行时某一周期内活塞二阶运动和活塞裙部摩擦学行为的差别.结果表明:该模型可应用于内燃机动力学和摩擦学行为的稳态及瞬态分析.     

An elasto-plastic theory of dislocation and disclination fields

A linear theory of the elasto-plasticity of crystalline solids based on a continuous representation of crystal defects – dislocations and disclinations – is presented. The model accounts for the translational and rotational aspects of lattice incompatibility, respectively associated with the presence of dislocations and disclinations. The defects content relates to the incompatible plastic strain and curvature tensors. The stress state is described by using the conjugate variables to strain and curvature, i.e., the stress and couple-stress tensors. Defect motion is described by two transport equations. A dynamic interplay between dislocations and disclinations results from a disclination-induced source term in the transport of dislocations. Thermodynamic guidance provides the driving forces conjugate to dislocation and disclination velocity in a continuous context, as well as admissible constitutive relations for the latter. When dislocation and disclination velocity vanish, the model reduces to deWit’s elasto-static theory of crystal defects. It also reduces to Acharya’s linear elasto-plastic theory for dislocation fields when the disclination density is ignored. The theory is intended for use in instances where rotational defects matter, such as grain boundaries. To illustrate its applicability, a finite high-angle tilt boundary is modeled using a disclination dipole and its behavior under tensile loading normal to the boundary is shown.     

Pulsatory Turbulent Compressible Fluid Flow and Propagation of Pressure Waves in a Channel

The properties of the damping coefficient and phase velocity of propagation of small-amplitude pressure waves as functions of the oscillation frequency are investigated for the turbulent flow of a weakly compressible fluid in a circular pipe. The wall friction is found by solving numerically the equation of motion and the relaxation equations for the turbulent shear stress and viscosity which provide the basis for a turbulent transfer model developed for unsteady conditions. The properties are explained in terms of an analysis of the calculated data on turbulent transfer. The results obtained are compared with experiments.     

La dynamique des tremblements de terre vue à travers le séisme de Landers du 28 juin 1992

We study the influence of initial conditions and of friction laws on the propagation of dynamic rupture during the earthquake of 28 June 1992 in Landers, California. We model this earthquake solving the elastodynamic wave equation by a finite difference method and we model observed accelerograms in order to get a better knowledge of the dynamic rupture process of this earthquake. In our model rupture propagates spontaneously under the simultaneous control of the initial stress field and friction. We model friction by a simple slip-weakening law. Finally, we inverted the initial stress field and the friction law from the radiation produced by 1992 Landers earthquake using a trial-and-error method. The dynamic model obtained by trial-and-error inversion provides a very satisfactory fit between synthetics and strong motion data. Rupture history and duration of the Landers earthquake are in good agreement with previous kinematic inversion results, without introducing major changes in final slip distribution on the fault. The solution of the dynamic inverse problem is non-unique because this problem is intrinsically ill-posed. Two complementary mechanical models were inverted in order to model the Landers earthquake, and to reproduce the seismic data. The first model corresponds to the asperity model in which only initial stress distribution is heterogeneous. The second model is a barrier model in which the initial stress was perfectly uniform while rupture resistance was heterogeneous. To cite this article: S. Peyrat et al., C. R. Mecanique 330 (2002) 235–248.