Research on pavement longitudinal crack propagation under non-uniform vehicle loading

The top-down crack (TDC) has become the major cracking mechanism in thick pavement structure, especially for the perpetual pavement. The heavy-vehicle load condition plays a key role in the fracture characteristics of pavement cracks. A three-dimensional finite element (FE) model of the tire tread rubber-block and the pavement is established to describe the stress–strain field of the pavement with TDC. The three-directional friction condition and non-uniform distribution between tire and pavement are especially considered. Then the orthogonal Design of Experiment method is applied to discuss the effect law of fracture characteristics for the longitudinal crack affected by the multiple loading parameters simultaneously. Based on the statistics and data analysis of the diverse test results, it is found there is a nonlinear relationship between the equivalent stress intensity factor of the pavement and the load parameters.The longitudinal distance (from the load location to the center of the crack port) has a great influence on the equivalent stress intensity factor. The equivalent stress intensity factor of the crack port is more than double the crack tip buried in pavement, which extends in II/III composite type mostly. The simulation results can be used as reference for the further study on the fracture mechanism of pavement cracks and their control technique.     

车、路的相互作用下沥青路面动力学特性分析

文中采用二自由度四分之一汽车悬架模型模拟汽车系统,依据弹性层状体系理论,建立路面结构的三维有限元分析模型,考虑车路相互作用,采用ANSYS有限元分析软件对移动车辆荷载作用下路面各结构层中的位移、应力、应变进行了模拟。计算分析了行车速度、悬架刚度、悬架阻尼、轮胎刚度和轮胎阻尼五个参数对路面动力响应的影响。结果表明：沥青面层处于三向受压状态,层内切应力是引起其破坏的主要原因;最大水平拉应力和最大横向拉应力均发生在路面结构的基层和底基层结合处;车速对路面动态响应的影响规律很复杂,应考虑车辆模型和路面不平整度,并划分速度区间加以探讨;路面动态响应随轮胎刚度、悬架刚度和悬架阻尼的增大而减小。上述结论对于深入分析路面结构动力响应与疲劳损坏以及研究车辆与路面相互作用的机理有重要价值。     

Backcalculation of pavement layer elastic modulus and thickness with measurement errors

This paper presents a backcalculation method for pavement layer elastic modulus and thickness. The effect of deflection measurement errors on the backcalculated results is also considered. The falling weight deflectometer (FWD) data are generated by applying a load to the pavement while calculating deflection at various fixed distances from the load centre. The measurement errors in FWD data are simulated by perturbing the theoretical deflections. Using these data, a backcalculation technique based on the improved genetic algorithm is proposed. In order to deal with the measurement errors, besides the common root mean square, a new objective function called area value with correction factor is introduced to the backcalculation algorithm. Numerical examples for two- and four-layer pavement structures are presented, which show the capability of the proposed method in backcalculation of pavement layer modulus and thickness.     

Effect of full-size and down-scaled accelerated traffic loading on pavement behavior

Accelerated pavement testing (APT) is an effective testing procedure to evaluate asphalt pavements. With APT it is possible to determine and measure the structural response and pavement performance under a controlled, accelerated damage accumulation in a compressed period of time. However, different types of APT technologies can lead to different results. Full-size loading devices simulate road traffic accurately, but are expensive, while down-scaled size simulators are cost effective, nevertheless further away from reality. In this work, two types of APT mobile load simulators with different loading characteristics are compared with respect to pavement response in the field and in the laboratory. The MLS10 is a full-size simulator, whereas the MMLS3 is a one-third scale device. The relationship between the devices was studied in terms of the measured strains induced by both machines in the same pavement. Therefore, a testing field was instrumented with strain gauges and first trafficked with MLS10. Later, a slab of the instrumented pavement was cut off the road and tested in the laboratory with the smaller MMLS3. Furthermore, the structure of the pavement was modelled with a viscoelastic finite element method model and the moving loads of both machines were simulated considering size, speed and approximate footprints of their tires. As for the pavement materials, the properties of the different asphalt layers were determined in the laboratory. Experimentally acquired strain data were used to validate the models. Stress fields under different loading and environmental conditions were analysed and compared. The evaluation shows that the models can predict the pavement response under different loading conditions. However, they still need to be improved to increase the accuracy under different conditions. Further, the analysis of the strains show that both load simulators induce a different stress–strain situation and scaling of the pavement should be considered.     

无网格法模拟复合型疲劳裂纹的扩展

本文提出了用无网格Galerkin法模拟构件在复合变形作用下疲劳裂纹扩展路径并预估其疲劳寿命的方法。该法能够自然模拟疲劳裂纹的扩展,不需要网格重构,避免了裂纹扩展过程中的精度受损。应用无网格数值结果计算了J积分和应力强度因子IK和IIK;按照最大周向应力理论获得了裂纹扩展偏斜角。基于最小应变能密度因子理论,确定了裂纹扩展量aD,并能获得疲劳载荷的循环周数ND。文末对数值模拟结果和实验拟合结果进行了对照。     

基于遗传神经网络的履带行驶系统载荷识别方法

针对煤矿掘进机器人履带行驶系统工作环境恶劣,载荷无法直接有效获取这一工程实际问题,提出了基于遗传神经网络的振动信号载荷识别方法。构建了遗传算法(GA)优化BP(back propagation)神经网络载荷识别模型,采用路试法试验采集了履带小车的5组振动加速度数据和单组应力载荷数据,探讨路面不平度频率和驱动轮啮频等对履带车振动和应力载荷的影响规律;借助快速傅里叶变换(FFT)对原始应力载荷数据进行去噪处理,依据履带小车行驶平顺性指标,利用sym8小波函数对振动加速度信号进行5层特征提取以提高载荷识别的精度,然后将5组小波变换分解的加速度数据和滤波后的应力载荷数据分别作为GA-BP神经网络的输入和输出进行训练及验证,揭示了履带行驶系统运动过程中振动与应力载荷之间的关系。研究结果表明,路面不平度频率、驱动轮啮频及转频为小车振动的主要频率成分,路面不平度引起的振动频率为13.765 Hz,驱动轮啮频为68.25 Hz,转频为3.25 Hz。多组试验得到的BP神经网络最佳隐含层神经元数为63,GA-BP神经网络识别的应力载荷与期望应力载荷具有较高吻合度,相对误差为4.5%,验证了该方法的有效性...     

Evaluation of the structure-induced rolling resistance (SRR) for pavements including viscoelastic material layers

This article addresses the question of whether the viscoelastic behavior of asphalt pavements, as part of the factors that impact the rolling resistance, might lead to substantially higher energy consumption of road traffic, as compared to non-dissipative pavements. In the context of sustainable development (fuel consumption, gas emission), this is a current issue for the pavement design community. This problem is tackled by a theoretical approach which consists in computing the structure-induced rolling resistance (SRR) for a vehicle by using the mechanical response of a layered viscoelastic medium excited by moving loads (i.e. tires of the vehicle). The structure-induced power dissipation is shown to be proportional to the integral over contours of the applied loads of the pavement deflection, times the outward normal to these contours. Following the developed approach, the determination of the SRR force engendered by the structure-induced dissipation is obtained from the computation of the deflection of the viscoelastic pavement to an assumed pressure distribution. Such an approach applied to a thick asphalt pavement allows the SRR for a typical road structure to be evaluated as a function of temperature and speed. A non-dimensional analysis is also presented to extend the results to some other contexts. Under the assumptions made for the derivation of the SRR, the main conclusion of the case study is that the extra fuel consumption for a vehicle resulting from the viscoelastic behavior of asphalt pavements is small.     

水泥混凝土路面角隅应力分析

鉴于我国刚性、半刚性基层上水泥混凝土路面角隅断裂多发的状况,采用弹性地基上不等平面尺寸双层结构模型分析其产生的原因。面层角隅在单个矩形荷载作用下,讨论了半空间地基和文克勒(Winkler)地基路面结构应力、弯沉的异同,给出了按应力和弯沉等效时地基参数之间的换算关系;分析了基层超宽对面层和基层自身应力的影响,得到了基层超宽系数的近似回归式;探讨了面层角隅受荷和温度梯度共同作用下面层应力采用叠加方法和耦合方法计算的差异;初步给出了板底脱空(层间接触刚度弱化)、邻板接缝传荷对路面结构应力、弯沉影响的规律。     

基于非线性模型的沥青路面动力响应研究

根据路基变形的非线性及沥青路面具有明显粘弹性的特点,将沥青路面简化为非线性粘弹性地基上的粘弹性无限长梁,建立了移动载荷作用下非线性粘弹性梁系统动力响应数学模型。利用Adomian分解法和小波变换法得到求解稳态响应解析解的新方法。通过实际道路参数对沥青路面动力响应进行了数值仿真,研究了车速、车辆轴载、路面材料及温度对沥青路面动力响应的影响规律。结果表明该文提出的计算方法简便、快捷,是求解该类非线性动力响应问题的一种有效方法;非线性模型更能准确地反映重载及超载时沥青路面实际结构受力状态。     

Dynamic response of top-down cracked asphalt concrete pavement under a half-sinusoidal impact load

A 3D finite element analysis model of cracked asphalt pavement is established by the FEM software ABAQUS. Based on dynamics mechanics, fracture mechanics and finite element theory, this paper studies the influence of various vehicle speeds, crack location, crack depth, damping ratio etc. on the dynamic response. The results show that the surface deflection, the maximum tensile strain at the bottom of the asphalt layer, and the maximum shear stress of the asphalt layer decreased with the increase in vehicle velocity when there is no crack in the pavement. No matter where the transverse position of the crack is the stress intensity factors increase with the increase in crack depth and decrease exponentially with the increase in longitudinal distance between the vehicle center and the crack. In the case of the crack locating in the center of wheel clearance, the surface deflection increases with the crack depth increasing. But if the crack is at the edge of the wheel track, there will be a critical vehicle velocity where the surface deflection is smaller than the asphalt pavement without crack if the vehicle velocity is above it. The maximum tensile strain at the bottom of the asphalt layer and the maximum shear stress of the asphalt layer are also smaller than the asphalt pavement without crack. The maximum tensile strain and the maximum shear stress decrease with the damping ratio increasing. So the increase in damping ratio can help to alleviate the propagation of cracks.     

动荷载作用下含裂缝公路结构体的应力强度因子

以沈阳-大连高速公路为工程背景，基于弹性动力学理论，采用平面应变有限单元法，分析了车辆荷载对含裂缝路面体的动态作用，分析过程中，车辆荷载简化为正弦分布柔性荷载；路面结构体计算模型抽象为平面应变模型；路面结构体为弹性的连续介质，为了反映裂尖应力，位移场的奇异性和减少模型网格数，在裂尖环向设置了奇异单元。通过计算得到裂尖的位移场，由位移外插得到I-型应力强度因子随加载时间的变化规律。同时探讨了初始裂缝长度和公路结构材料阻尼比的变化对I-型应力强度因子分布规律的影响，为路面体的动态破坏研究提供了一定的理论参考。     

Application of the weight function method to two-dimensional elastodynamic fracture mechanics

In this paper, the weight function method is used for two-dimensional mixed-mode crack analyses of clastostatic and elastodynamic problems. By the use of the Laplace transformation method and an indirect boundary element method, the dynamic stress intensity factors for a finite sheet containing a central or an edge crack are evaluated. A Green's function method is introduced which depends on the weight function for an impulsive applied load. The Green's function can be used to determine stress intensity factors for arbitrary time dependence of the boundary conditions. The stress intensity factors obtained by the weight function method are compared where possible, with existing solutions.     

Statistical modelling of predicted non‐stationary vehicle vibrations

This paper describes a method to predict the vertical vibrations of road vehicles from measured pavement profiles. It discusses the limitations of current methods used for analysing and simulating vehicle vibrations and shows that more accurate characterization and simulation of the transport environment must take into account the non‐stationary nature of road vehicle vibrations. Vertical vibrations for typical transport vehicles under various operating conditions and pavement profiles are predicted using a computer model of the vehicle characteristics and analysed to produce the spectral and statistical characteristics. The paper also presents an improved method to compute the vibration intensity by using a dynamic segmentation data reduction technique. The effectiveness of the procedure to characterize the non‐stationarity of random vehicle vibrations is demonstrated. Finally, the paper deals with the statistical distribution of the vibration intensity and demonstrates how it can be adapted to a technique for the simulation the non‐stationary nature of random vehicle vibrations. Copyright © 2002 John Wiley & Sons, Ltd.     

Stress intensity factor at the surface and at the deepest point of a semi-elliptical surface crack in plates under stress gradients

The weight function method is used to calculate stress intensity factors for a semi-elliptical surface crack in a plate exposed to stress gradients. Starting from a reference load and stress intensity factor an approximate reference displacement field is calculated analytically. The present method allows to calculate stress intensity factors with minimal numerical effort at the deepest point and at the surface. Comparisons with FEM-results from the literature are presented to show satisfying agreement.     

Constant stress intensity factors through closed-loop control

A new way of obtaining a constant stress intensity factor is achieved for any test specimen geometry subjected to closed-loop control loading. In contrast to using only load or displacement control the method draws on combined feedback from both displacement and load sensing, reducing the variation in the stress intensity factor by two decades compared to that if tested under constant load or displacement. A change in the signal mix ratio for a rectangular compact tension specimen is equivalent to changing the angle of a tapered compact tension specimen. This method can eliminate the need for the use of ‘complex’ geometries or for geometries in which measured crack lengths are used in conjunction with a computer or some other means to adjust the loading continuously for achieving a constant stress intensity factor.     

基于路面等级匹配的重卡包装物运输激励研究

目的利用传感器搭配数据采集卡收集路面不平频谱,将路面信息匹配标准路面等级以及车辆行驶速度。方法通过Matlab将标准PSD与随机激励匹配路面等级,同时建立1/2四轴重卡动力学模型,并用能量法建立动力学方程。将收集的路面不平频谱对应到相应的路面等级,再结合车辆速度的设定,最后求解得到车辆受激励后轮胎的动载位移频谱,分析得出被运输包装物的半挂车平板动载位移。结果重卡运输前轴轮胎在A级路面以60 km/h的车速经过该路面的动载位移量在0.8和9.8 s时达到峰值,且路面响应位移不超过6 cm。结论求得被运输包装物的所受激励频谱,为被包装物的运输振动安全性研究提供支撑,可结合具体被运输包装物的脆值理论,提供被运输物品发生运输损坏的数值仿真。     

Evaluation of fracture parameters in continuously nonhomogeneous piezoelectric solids

A contour integral method is developed for computation of stress intensity and electric intensity factors for cracks in continuously nonhomogeneous piezoelectric body under a transient dynamic load. It is shown that the asymptotic fields in the crack-tip vicinity in a continuously nonhomogeneos medium is the same as in a homogeneous one. A meshless method based on the local Petrov-Galerkin approach is applied for computation of physical fields occurring in the contour integral expressions of intensity factors. A unit step function is used as the test functions in the local weak-form. This leads to local integral equations (LBIEs) involving only contour-integrals on the surfaces of subdomains. The moving least-squares (MLS) method is adopted for approximating the physical quantities in the LBIEs. The accuracy of the present method for computing the stress intensity factors (SIF) and electrical displacement intensity factors (EDIF) are discussed by comparison with available analytical or numerical solutions.     

On design methods for bolted joints in composite aircraft structures

The problems related to the determination of the load distribution in a multirow fastener joint using the finite element method are discussed. Both simple and more advanced design methods used at Saab Military Aircraft are presented. The stress distributions obtained with an analytically based method and an FE-based method are compared. Results from failure predictions with a simple analytically based method and the more advanced FE-based method of multi-fastener tension and shear loaded test specimens are compared with experiments. Finally, complicating factors such as three-dimensional effects caused by secondary bending and fastener bending are discussed and suggestions for future research are given.     

An Interaction Integral Method for Computing Fracture Parameters in Functionally Graded Magnetoelectroelastic Composites

A contour integral method is developed for the computation of stress intensity, electric and magnetic intensity factors for cracks in continuously nonhomogeneous magnetoelectroelastic solids under a transient dynamic load. It is shown that the asymptotic fields in the crack-tip vicinity in a continuously nonhomogeneos medium are the same as in a homogeneous one. A meshless method based on the local Petrov-Galerkin approach is applied for the computation of the physical fields occurring in the contour integral expressions of intensity factors. A unit step function is used as the test functions in the local weak-form. This leads to local integral equations (LIEs) involving only contour-integrals on the surfaces of subdomains. The moving least-squares (MLS) method is adopted for approximating the physical quantities in the LIEs. The accuracy of the present method for computing the stress intensity factors (SIF), electrical displacement intensity factors (EDIF) and magnetic induction intensity factors (MIIF) are discussed by comparison with numerical solutions for homogeneous materials.     

Interface fracture properties of a bimaterial ceramic composite

In this investigation, the interface fracture toughness is measured for a pair of ceramic clays which are joined together. The Brazilian disk specimen, which provides a wide range of mode mixity, is employed to measure these properties. Calibration equations relating the stress intensity factors to the applied load and geometry are determined by means of the finite element method and the M-integral. The effect of residual stresses is accounted for by employing a weight function to obtain the contribution to the stress intensity factors. Total stress intensity factors are obtained by superposition. These are employed to determine the critical interface energy release rate as a function of mode mixity from critical data obtained from tests carried out on the Brazilian disk specimens. An energy release rate fracture criterion is compared to the experimental results for .