Study of the influence of vehicle parameters on tractive performance in deep snow

This paper describes an experimental study of tractive performance in deep snow, carried out with a new special skid steered tracked vehicle, developed by Bodin [1]. The vehicle design parameters studied include the influence of the ground clearance of the vehicle belly and the longitudinal location of the centre of gravity on tractive performance in deep snow, as well as the effect of initial track tension. The most important results from the test show that an increase in the ground clearance has a positive effect on the drawbar pull, originating from a greater increase in the thrust than in the track motion resistance and a slight decrease in the belly drag. Tests of the longitudinal location of the centre of gravity show that a location ahead of the midpoint of the track contact length is to be preferred. The drawbar pull increases with the centre of gravity moving forward. This is due to a reduced track motion resistance, a slight decrease in the belly drag and an almost constant vehicle thrust. The reason for the decreased track motion resistance and belly drag with the centre of gravity located ahead of the midpoint of the track contact length is a decreased vehicle trim angle.     

Tractive and braking performance of a flexible tracked tractor moving up and down weak sloped terrain

The objective of this study is to analyse the tractive and braking performance of a tractor travelling up and down a weak silty loam sloped terrain. The effects of track belt size on terrain-track system parameters were investigated experimentally, and the force and energy balances were clarified for the actual flexible tracked vehicle. The flexibility of the track belt has been precisely analysed as a function of track tension, loading and reloading properties of terrain, and contact pressure distributions. The results show that the optimum effective driving (or braking) force decreases with the increase of slope angle due to the decreasing vehicle weight component, while the thrust (or drag) decreases and the compaction resistance increases. The contact pressure distribution under the flexible track belt shows a triangular wavy pattern having peak values under each track roller. The shear resistance distribution has positive and negative peak values for the driving and braking states, respectively.     

Tractive performance of a bulldozer running on weak ground

To determine the tractive performance of a bulldozer running on weak ground in the driven state, the relations between driving force, drawbar pull, sinkage, eccentricity and slip ratio have been analysed together with each energy balance; effective input energy, sinkage deformation energy, slippage energy and drawbar pull energy. It is considered that the thrust is developed not only on the main straight part of the bottom track belt but also on parts of the front idler and rear sprocket, and the compaction resistance is calculated from the amount of slip sinkage. For a given vehicle and soil properties, it is determined that the drawbar pull increases directly with the slip ratio and reaches about 70% of the maximum driving force. The compaction resistance reaches about 13% of the maximum driving force. The sinkage of the rear sprocket, the eccentricity, and the trim angle increase with the increment of slip ratio due to the slip sinkage. These analytical results have been verified experimentally. After determining the optimum slip ratio to obtain a maximum effective tractive power, it is found that a larger optimum drawbar pull at optimum contact pressure could be obtained for a smaller eccentricity of vehicle center of gravity and a larger track length-width ratio under the same contact area.     

考虑摩擦因素的结合面切向接触刚度分形预估模型及其仿真分析

为建立更完善和精确的结合面接触刚度模型,本文根据分形理论和摩擦学原理,从微观角度建立了考虑摩擦因素的结合面切向接触刚度分形预估模型.通过数值仿真分析研究了接触载荷、分形维数、摩擦系数和接触面积等因素对结合面切向接触刚度的影响.分析结果表明:结合面切向接触刚度随法向载荷和分形维数的增加而增大,而随分形尺度参数的增大而减小;摩擦系数对结合面切向接触刚度的影响较大,不同实际接触面积下的切向刚度相差较大;当分形维数较小时,摩擦系数对结合面切向刚度的影响将降低.这些研究对于进一步开展结合面的动力学特性研究具有重要意义.     

滑轨导向式静/动态双轴拉伸实验技术

基于液压伺服高速加载系统,发展了一种材料双轴拉伸力学性能测试技术。利用锥面接触导向驱动方法,把加载锤竖直方向的驱动力转化为水平方向的双轴驱动力,从而实现对十字形试样平面双轴加载。借助有限元数值模拟手段优化了锥面接触角和十字形试样尺寸。当接触锥角为45°时,既有较好的水平驱动转化效率,同时又保持较小的接触力,确保水平驱动加载各组件在弹性变形范围内,可多次重复使用。确定了加载臂狭缝个数、狭缝与减薄区边缘长度和标距段厚度等试样设计关键参数,在十字形试样测试标距段内实现了均匀变形。设计了测力夹持一体化导杆和非接触光学全场应变测试系统,准确获得了试样的应力和应变。利用此平面双轴拉伸加载装置,开展2024-T351铝合金板单轴拉伸实验和激光探测同步性验证实验,验证装置设计的可行性;开展铝合金板材在不同加载速率下的双轴拉伸实验,得到在双轴加载下铝合金板材应力应变曲线,并与单轴加载下实验结果进行了对比分析。     

粗糙表面接触分形模型的提出与发展

自表面接触塑性变形模型问世以来，经过近４０年的发展，已经形成以分形几何理论为基础的Ｍ－Ｂ粗糙表面接触分形模型．Ｍ－Ｂ模型以分形参数代替统计学参数表征粗糙表面，推导出了实际接触面积与载荷的关系，以及实际弹性接触面积和实际塑性接触面积的计算公式，指出了影响接触面变形性质的因素与规律．由于分形参数的尺度独立性，可望利用Ｍ－Ｂ模型对接触面积的预测不受测量仪器分辨率和取样长度等因素的影响，故其比基于统计分析的Ｇ－Ｗ接触模型更为合理．尽管如此，Ｍ－Ｂ模型还有待完善，多方面的问题尚待进行深入研究     

载重子午线轮胎帘线受力的有限元分析

以10.00R20轮胎为原型,针对特定的断面轮廓形状和帘线结构形式建立了轮胎结构计算的有限元模型,着重研究了额定充气压力状态以及受静负荷作用时各部分帘线受力的基本特征。结果表明:(1)充气状态下胎冠部和胎圈部的帘线结构参数之间的相互影响很小,因而在分析设计中可以相对独立地考虑;(2)静负荷作用下帘线受力在接地面附近的分布情况较复杂,尤其是负荷较大时可能出现胎冠中部帘线张力不足甚至受压以及0°带束层边缘受拉力过大的情况,这些对轮胎使用性能不利的受力状况可以通过帘线结构的合理设计来克服。     

横摆角下汽车尾部湍流特征量的PIV试验分析

汽车尾部湍流场是汽车压差阻力的主要来源,在HD-2汽车模型风洞中,首先使用测力天平和测压系统,对横摆角工况下汽车模型的气动六分力和纵对称截面48个测点的表面压力进行了测量,然后利用PIV测量技术对模型在横摆角分别为0°、15°的尾部湍流场进行了测量,获得该模型尾流场的速度场、涡量场和雷诺应力流场信息,通过计算得出尾流场区域空间相关系数和湍流积分尺度。结果表明:在横摆角工况下,汽车模型尾部涡流的结构呈现向上发展的趋势;尾流场拖拽涡的范围和强度的增大导致了模型气动力出现较大的增加;湍流积分尺度的变化表明,尾部涡流区的分离噪声与涡流分离位置有关,在汽车尾部造型设计中,要尽量推迟尾部涡流的分离。     

三峡库区滑坡的时空分布特征与成因探讨

三峡库区地质灾害频繁发生,地质历史上大型滑坡广泛发育。近年库区滑坡灾害有逐年增加趋势。这些滑坡形成的堆积体成因复杂,不能单纯从传统的外动力或内动力作用来解释,给库区灾害的防治工程带来了很大困难。三峡库区大型古滑坡主要沿长江干流和支流的深切河谷两岸分布,与河流发育演化史密切相关。15万年以来,长江三峡地区经历了新构造期以来抬升速率最快的一次构造幕,山原期夷平面上升了100m,河流强烈下切,形成高陡岸坡,卸荷效应显著;同时,该时期三峡地区的气候受到青藏高原环境变化的影响,在间冰期可能存在相似的暖湿古气候环境,降水极为丰富。三峡库区的大型古滑坡在新构造快速抬升时期广泛发育,并且与暖湿多雨期相对应,表明库区滑坡的发育演化是新构造运动和气候变化内外动力耦合作用的结果。     

Some further studies on the mechanical properties of muskeg in relation to vehicle mobility

An experimental study of the mechanical properties of two types of muskeg was carried out. The results have further verified the basic features of a mathematical model for the pressure-sinkage characteristics of muskeg proposed previously. The procedure for deriving the value of the stiffness of the underlying peat from load-sinkage tests has been validated. The response of the muskeg to repetitive loading was measured. It has been found that these muskeg types exhibit noticeable elastic rebound and hysteresis during the unloading-reloading cycle. The average slope of the unloading and reloading line appears to increase with the load level at which the unloading-reloading cycle begins. The experimental results further confirm that from thw vehicle mobility as well as environmental protection standpoints, the tearing-off of the muskeg surface mat by vehicle running gear is to be avoided. When the mat is torn, the vehicle running gear will be in contact with the underlying peat which has lower bearing capacity and shear strength. This results in higher sinkage and lower traction.     

Development of a two-dimensional model of a compliant non-pneumatic tire

An analytical model for a compliant non-pneumatic tire on frictionless, rigid ground is presented. The tire model consists of a thin flexible annular band and spokes that connect the band to a rigid hub. The annular band is modeled using curved beam theory that takes into account deformations due to bending, shearing and circumferential extension. The effect of the spokes, which are distributed continuously in the model and act as linear springs, is accounted for only in tension, which introduces a nonlinear response. The quasi-static, two-dimensional analysis focuses on how the contact patch, vertical tire stiffness and rolling resistance are affected by the stiffness properties of the band and the spokes. A Fourier series representation of the shear strain in the annular band and the complex modulus of the material were used to predict rolling resistance due to steady state rolling. From the analysis point of view, when the wheel is loaded at its hub, the following three distinct regions develop: (1) a support region where the hub hangs by the spokes from the upper part of the flexible band, (2) a free surface region where the spokes buckle and have no effect, and (3) a contact region where the flexible band is supported by the ground without the effect of the spokes. The angular bounds of these three regions are determined by the spoke angle and the contact angle, which are respectively the angle at which the spokes start to engage in tension and the angle that defines the edge of contact. Closed-form expressions of contact stress, stress-resultants and displacements at the centroids of the cross-sections of the flexible band are expressed in terms of these angles, which must be determined numerically. A thorough parametric analysis of quantities of interest for the tire is presented, which can be used to help support the optimal and rational design of compliant non-pneumatic tires. The model was validated by comparison with two computational models using the commercial finite element software ABAQUS and by experimental rolling resistance data.     

基于碰撞模型的斜坡滚石颗粒速度预测

滑坡滚石灾害是西部山区常见的地质灾害类型,具有突发性和随机性强的特点,是山区地质灾害预测和防治工作的重点和难点. 本文基于颗粒接触理论,考虑影响斜坡滚石碰撞过程中的随机因素,建立了用于预测斜坡滚石颗粒碰撞后速度的理论模型. 根据冲量及冲量矩定理建立滚石颗粒碰撞基本方程,得到斜坡滚石颗粒碰撞后反弹速度的解析解. 结果表明:斜坡滚石碰撞后反弹速度的解析解包含了坡角、坡体上被碰颗粒速度以及角度、入射速度和角度以及撞击角度等随机因素;当考虑入射滚石颗粒与坡体上被碰颗粒的撞击角度变化时,模型预测结果与试验结果吻合较好;本文进一步预测了滚石颗粒碰撞后颗粒反弹线速度和角度以及反弹旋转角速度的概率分布情况. 结果显示,反弹颗粒速度和角度以及反弹旋转角速度的概率分布均服从高斯分布;当坡体上被碰颗粒速度和坡角发生变化时,其对反弹颗粒速度和角度以及反弹旋转角速度概率分布定性上没有影响,但是对概率分布的中心参数有显著影响.      

砂土地基锚板基础抗拔承载力PFC数值分析

在模型试验的基础上,采用PFC离散单元法对条形锚板基础在中密砂土地基中的抗拔性能进行了数值分析。数值模拟采用簇颗粒单元来模拟砂粒的不规则形状,颗粒级配根据模型试验福建标准砂的级配按照相似级配法生成,细观参数根据数值双轴试验确定。水平锚板数值模拟结果与模型试验结果基本一致。与模型试验结果相比,颗粒流数值模拟能得到颗粒间接触力链的分布及其演化规律,能从细观角度来探明宏观抗拔承载力特性的演化机理。在此基础上对倾斜锚板上拔过程进行了模拟,分析了锚板前后砂粒的运动趋势以及接触力链的演化规律,并与已有承载力结果进行了对比分析。     

CFD–DEM simulation of particle deposition characteristics of pleated air filter media based on porous media model

In this study, the three-dimensional physical model of pleated air filtration media was simplified to porous media model, and the calculation parameters of porous media were obtained based on experimental data. The model of V-shaped pleated air filter media is constructed, the height of the media pleat is 50 mm and the pleat thickness is 4 mm, the pleat angle is 3.7°. The Hertz-Mindlin contact model was modified by Johnson Kendall Roberts (JKR) adhesion contact model. The deposition process of particles in media was simulated based on computational fluid dynamics (CFD) theory and discrete element method (DEM). Results show that the CFD–DEM coupling method can be effectively applied to the macro research of pleated air filter media. The particles will form dust layer and dendrite structure on the fiber surface, and the dust layer will affect the subsequent air flow organization, and the dendrite structure will eventually form a “particle wall”. The formation of the “particle wall” will prevent the particles from moving further in the fluid domain, which makes area of pleated angle become the “low efficiency” part about the particle deposition. Compared with area of pleated angle, the particles are concentrated in the opening area and the middle area of the pleated to agglomerate and deposit.     

Crack kinking from an interface crack with initial contact between the crack lips

The authors recently theoretically studied crack kinking and opening from an initially closed crack (without friction) in some homogeneous medium. The same problem, but for an interface crack, is considered here. Comninou has shown that the asymptotic stress field prior to kinking is governed by a single, mode II stress intensity factor (SIF). Using this result, plus a homogeneity property of the problems of elastic fracture mechanics with unilateral contact envisaged, a change of scale, and two reasonable hypotheses, we establish the expression of the SIF at the tip of the small, open crack extension. It is shown that whatever the geometry of the external boundary and the crack and whatever the loading, these SIF depend solely upon the initial (mode II) SIF (in a linear way), the kink angle and Dundurs' parameters α and β. Using this result and Goldstein and Salganik's “principle of local symmetry” to predict the kink angle, one finds that it is independent of the loading but does depend upon Dundurs' parameters α and β. This contrasts with the case of an ordinary (initially closed) crack in some homogeneous medium, for which the kink angle was not only independent of the loading but an absolute constant. However, it is numerically found that the influence of the mismatch of elastic properties upon the kink angle is rather weak.     

凹边坡参数对边坡稳定性的影响

抽象出凹边坡的几何模型,采用量纲分析的方法,给出凹边坡几何参数和物理力学参数对边坡稳定系数的影响通式．在FLAC定义的边坡稳定系数意义下,经过数万次计算,采用二次回归的方法,确定了通式中的影响指数,形成经验公式．结果表明,受Mohr-Coulomb和拉应力两个强度准则控制的材料,形成的凹边坡的稳定性受无量纲粘聚力、边坡角、内摩擦角的控制．凹形边坡有利于边坡的稳定,边坡稳定性与岩土的抗拉强度关系不大．这一稳定系数公式对于形状规则均匀性好且符合Mohr-Coulomb破坏准则的凹边坡可以给出精度高的估计．对于其它边坡的稳定性评价还需要做有针对性的工作．     

基于基底长度的粗糙表面分形接触模型的构建与分析

为建立更为精确的粗糙表面接触模型,根据微凸体变形特征、分形理论以及摩擦的作用,从微观角度基于基底长度建立了粗糙表面分形接触模型.通过与其他粗糙表面接触模型和实验数据的比较,验证了本文模型的正确与合理,并由数值仿真分析了分形维数、接触载荷与真实接触面积之间的相互关系.分析结果表明:特征尺度一定时,要维持一定的真实接触面积,分形维数越大,所需要的力也越大;分形维数与特征尺度一定时,随着载荷的增加,接触面积也在增加;特征尺度与接触力一定时,随着分形维数的增大,真实接触面积在减小.该模型的建立为进一步研究粗糙表面的摩擦、磨损与润滑提供了理论依据.     

Infiltration of liquid droplets into porous media: Effects of dynamic contact angle and contact angle hysteresis

We present a detailed theory for infiltration, which accounts for a general model for the dynamic contact angle between the droplet and the porous medium as well as contact angle hysteresis, and analyze the resulting equations of motion. The theory shows that infiltration of droplets into dry porous media involves three phases due to contact angle hysteresis: (1) An increasing drawing area (IDA) phase during which the interface between the droplet and the porous medium increases, (2) a constant drawing area (CDA) phase during which the contact line of the droplet remains pinned, and (3) a decreasing drawing area (DDA) phase. The theory is based on the following assumptions: (1) The droplet has the shape of a spherical cap, (2) the porous medium consists of a bundle of vertical tubes of same size, and (3) the pressure within the droplet is uniform. We find that infiltration always consists of a cascade process formed by the IDA, CDA, and DDA phases, where the entire process may begin or end in any of the three phases. The entire process is formulated with four nondimensional parameters: Three contact angles (initial, advancing, and receding) and a porous permeability parameter. A comparison of our theory to experimental data suggests that one should use different parameterizations for the dynamic contact angle models of the IDA and DDA phases. In general, the IDA and DDA phases are described by integro-differential equations. A numerical-solution approach is presented for solving the dynamic equations for infiltration. The total time of infiltration and the time dependence of drawing area are critically affected by the occurrence of the IDA, CDA, and DDA phases as well as by the permeability. With ordinary differential equations (ODEs), we are able to approximate the IDA phase and to describe exactly infiltration processes that start out with the CDA or DDA phase.     

非线性Mohr－Coulomb体边坡的极限荷载

本文根据拟线性偏微分方程组的特征线理论，求得无重双曲、抛物Ｍｏｈｒ－Ｃｏｕｌｏｍｂ体边坡极限荷载的解析解．用两参数的双曲线和抛物线可以较好拟合由实验提供的破坏曲线．由其解发现屈服准则和坡角对极限荷载影响很大，理想塑性体的极限荷载也非最小．     

Traffic performance of a rubber-tracked vehicle travelling up and down a sloped pavement

The tractive and braking performances of a 40 kN rubber-tracked vehicle travelling up and down a sloped pavement depend on the grouser shape. The purpose of this paper is to find the most suitable grouser shape to obtain the maximum optimum effective tractive effort and the maximum optimum effective braking force and to clarify the several traffic performances of the vehicle travelling up and down sloped concrete and asphalt paved roads. As results, it is verified that the most suitable shape of rubber grouser is an equilateral trapezoid type of contact length 3 cm for concrete pavement and another of contact length 5 cm for asphalt pavement, respectively, and that the effective tractive effort and the effective braking force decrease with the increment of slope angle.