层间接触对沥青路面粘弹性动力响应影响研究

为分析沥青面层材料粘弹特性及路面各层间接触条件对沥青路面动力响应的影响,基于解析的方法,开展了层间非完全连续沥青路面粘弹性动力响应的求解工作．采用修正的Burgers模型定义沥青面层材料的粘弹性本构关系,考虑沥青路面层间接触条件,在车辆荷载作用下,建立沥青路面的理论计算模型;通过Laplace-Hankel积分变换将偏微分方程组转化为常微分方程组并对问题进行求解;采用转换矩阵表征层间接触条件,求得层间非完全连续沥青路面粘弹性动力响应的解析表达式．从沥青路面实例计算结果发现：修正的Burgers模型中的瞬时弹性模量参数是对弯沉计算结果影响最大的因素之一,路表弯沉随的增大呈下降趋势,特别是当较小时,这种趋势尤为明显;沥青混合料本构模型中的粘弹性修正系数B和黏性参数是影响路面路表弯沉计算结果的另两个重要因素,并且B和对弯沉峰值出现时间的影响具有相反的趋势;层间的非完全连续条件对沥青路面动态弯沉计算结果影响较大,并以面层与基层间的非完全连续对弯沉计算结果影响最为显著．     

Fractional Derivative Viscoelasticity at Large Deformations

A time domain viscoelastic model for large three-dimensional responses underisothermal conditions is presented. Internal variables with fractional orderevolution equations are used to model the time dependent part of the response. By using fractional order rate laws, the characteristics of the timedependency of many polymeric materials can be described using relatively fewparameters. Moreover, here we take into account that polymeric materials are often used in applications where the small deformations approximation does nothold (e.g., suspensions, vibration isolators and rubber bushings). A numerical algorithm for the constitutive response is developed and implemented into a finite element code forstructural dynamics. The algorithm calculates the fractional derivatives by means of the Grünwald–Lubich approach.Analytical and numerical calculations of the constitutive response in the nonlinearregime are presented and compared. The dynamicstructural response of a viscoelastic bar as well as the quasi-static response of athick walled tube are computed, including both geometrically and materiallynonlinear effects. Moreover, it isshown that by applying relatively small load magnitudes, the responses ofthe linear viscoelastic model are recovered.     

Discrete fracture modeling of asphalt concrete

A heterogeneous fracture approach is presented for modeling asphalt concrete that is composed of solid inclusions and a viscous matrix, and is subjected to mode-I loading in the fracture test configuration. A heterogeneous fracture model, based on the discrete element method (DEM), is developed to investigate various fracture toughening mechanisms of asphalt materials using a high-resolution image processing technique. An energy-based bilinear cohesive zone model is used to model the crack initiation and propagation of materials, and is implemented as a user-defined model within the discrete element method. Experimental fracture tests are performed to investigate various fracture behavior of asphalt concrete and obtain material input parameters for numerical models. Also, bulk material properties are necessary for each material phase for heterogeneous numerical models; these properties are determined by uniaxial complex modulus tests and indirect tensile strength tests. The main objective of this study is to integrate the experimental tests and numerical models in order to better understand the fracture mechanisms of asphaltic heterogeneous materials. Experimental results and numerical simulations are compared at different test conditions with excellent agreement. The heterogeneous DEM fracture modeling approach has the potential capability to understand various crack mechanisms of quasi-brittle materials.     

黏弹性体界面裂纹的冲击响应

研究两半无限大黏弹性体界面Griffith裂纹在反平面剪切突出载荷下,裂纹尖端动应力强度因子的时间响应,首先,运用积分变换方法将黏弹性混合黑社会问题化成变换域上的对偶积分方程,通过引入裂纹位错密度函数进一步化成Cauchy型奇异积分方程,运用分片连续函数法数值求解奇异积分方程,得到变换域内的动应力强度因子,再用Laplace积分变换数值反演方法,将变换域的解反演到时间域内,最终求得动应力强度因子的时间响应,并对黏弹性参数的影响进行分析。     

沥青混合料增量型热粘弹性本构关系研究

基于热粘弹性力学理论,就不同的温度条件下沥青混合料的应力松弛特征开展了试验研究,用广义Maxwell模型模拟沥青混合料的粘弹特性,应用热流变简单材料的时温等效原理对实验结果进行了分析和参数拟合,然后通过理论推导得到沥青混合料非定常和非均匀变温条件下增量型热粘弹性本构关系;在此基础上,给出利用此增量型应力应变本构关系进行沥青路面热粘弹性力学分析的数值实现方法,并给出一个计算实例。     

Nonlinearly Viscoelastic Nanoindentation of PMMA Under a Spherical Tip

The Oliver-Pharr method has been well established to measure Young’s modulus and hardness of materials without time-dependent behavior in nanoindentation. The method, however, is not appropriate for measuring the viscoelastic properties of materials with pronounced viscoelastic effects. One well-known phenomenon is the formation of unloading “nose” or negative stiffness during unloading that often occurs during slow loading-unloading in nanoindentation on a viscoelastic material. Most methods in literature have only considered the loading curve for analysis of viscoelastic nanoindentation data while the unloading portion is not analyzed adequately to determine the nonlinearly viscoelastic properties. In this paper, nonlinearly viscoelastic effects are considered and modeled using the nonlinear Burgers model. Nanoindentation was conducted on poly-methylmethacrylate (PMMA) using a spherical indenter tip. An inverse problem solving approach is used to allow the finite element simulation results to agree with the nanoindentation load–displacement curve during the entire loading and unloading stage. This approach has allowed the determination of the nonlinearly viscoelastic behavior of PMMA at submicron scale. In addition, the nanoindentation unloading “nose” has been captured by simulation, indicating that the negative stiffness in the viscoelastic material is the result of memory effect in time-dependent materials.     

Second-order two-scale computational method for ageing linear viscoelastic problem in composite materials with periodic structure

The correspondence principle is an important mathematical technique to compute the non-ageing linear viscoelastic problem as it allows to take advantage of the computational methods originally developed for the elastic case. However, the correspondence principle becomes invalid when the materials exhibit ageing. To deal with this problem, a second-order two-scale (SOTS) computational method in the time domain is presented to predict the ageing linear viscoelastic performance of composite materials with a periodic structure. First, in the time domain, the SOTS formulation for calculating the effective relaxation modulus and displacement approximate solutions of the ageing viscoelastic problem is formally derived. Error estimates of the displacement approximate solutions for SOTS method are then given. Numerical results obtained by the SOTS method are shown and compared with those by the finite element method in a very fine mesh. Both the analytical and numerical results show that the SOTS computational method is feasible and efficient to predict the ageing linear viscoelastic performance of composite materials with a periodic structure.     

黏弹性传动带1:3内共振时的周期和混沌运动

研究了参数激励作用下黏弹性传动带在1:3内共振时的周期解分岔和混沌动力学. 同时考虑传动带的线性外阻尼因素和材料内阻尼因素. 首先建立了具有线性外阻尼情况下的黏弹性传动带平面运动时的非线性动力学方程, 黏弹性材料的本构关系用Kelvin模型描述. 然后考虑黏弹性传动带的横向振动问题, 利用多尺度法和Galerkin离散法得到黏弹性传动带系统在1:3内共振时的平均方程. 最后利用数值模拟方法研究了黏弹性传动带系统的周期振动和混沌动力学, 得到了系统在不同参数下的混沌运动. 数值模拟结果说明黏弹性传动带系统存在周期分岔, 概周期运动及混沌运动.     

Predictions of nonlinear viscoelastic behavior using a hybrid approach

An approach for nonlinear viscoelastic characterization is presented which uses the combined measurements from creep and dynamic mechanical tests. Although the methodology should extend to several materials and geometries, this research concentrates on thin film polymers used in the manufacture of high altitude scientific balloons. Typically, the constitutive behavior of these materials is characterized through the use of linear viscoelastic techniques. Although this linear approach provides an accurate model for small strains or loads, these materials have been shown to be highly stress dependent and, consequently, it is necessary to identify this nonlinear behavior. Traditional creep measurements require extensive laboratory test times, yet the results obtained from dynamic mechanical analysis provide the capability to predict long term material performance without a lengthy experimentation program. However, dynamic mechanical methods are currently limited to linear response; thus, an approach is presented in which the stress-dependent behavior is derived from short-term creep measurements in a manner analogous to time-temperature superposition. Predictions of material response using linear and nonlinear approaches are compared with experimental results obtained from traditional long-term creep tests. Although linear pre-dictions deteriorate for large stresses, excellent agreement is shown for the nonlinear model.     

Homogenization of the viscoelastic heterogeneous materials with multi-coated reinforcements: an internal variables formulation

In this work, a new homogenization method to estimate the effective behavior of viscoelastic heterogeneous materials with multi-coated reinforcements is presented. Unlike classical methods that are based on the Laplace transform, the present internal variables formulation operates directly in the time domain. Using the Green’s function techniques, the micromechanical approach is based on establishing a new integral equation adapted to scale transition methods. Using this integral equation, we apply a generalized self-consistent scheme to determine the local stress concentration equations and the effective behavior of multi-coated inclusion-reinforced materials. To assess the reliability of our model, some applications to the isotropic viscoelastic heterogeneous materials with homothetic spherical inclusions are given. The model is applied to the case of two-phase and three-phase materials, and the results are compared to exact solutions. Results for three-phase materials are presented regarding the influence of soft and stiff viscoelastic interphase on the effective behavior of heterogeneous materials.     

Determination of Asphalt Binder Creep Compliance Using Depth-Sensing Indentation

Asphalt binders are common construction materials, however due to time- and temperature- dependence, their mechanical properties are often difficult to characterize. Several standard tests methods exist to describe their complex behavior. This paper presents an exploratory feasibility study of a flat-tip indentation testing to analyze the linear viscoelastic properties of asphalt binders. Depth-sensing indentation testing has been extensively used to characterize the properties of many engineering materials, however the applications to asphalt binders are very limited. This paper presents a simple solution for the creep compliance in tension derived for flat-tipped indenter. This solution was verified with the Finite Element Analysis and then applied to the experimental results from the indentation testing performed on one typical unmodified asphalt binder. The testing was conducted at three different low temperatures and under three different creep load levels to verity the linearity of the response, and to evaluate the robustness and applicability of the indentation method. Furthermore, the creep compliance functions determined from the indentation testing were compared with a more traditional 3-point bending experiments. The results show that there is a non-uniform discrepancy between the two testing methods, most likely due to nonlinear behavior of the asphalt binder at higher temperatures and micro-damage of the binder samples at lower temperatures. Other possible sources of error between indentation and 3-point bending are problems determining the initial tip-specimen contact surface and possible tip-specimen adhesion. It is concluded that flat-tipped indentation at low temperatures should be performed at lower load levels to avoid excessive stress concentrations that leads to micro-damage and nonlinear response of asphalt binders. Alternatively, asphalt binders at low temperatures could be evaluated using different indenter geometries, such as spherical or pyramidal, using corresponding parameter interpretation procedures.     

GENERAL SECOND ORDER FLUID FLOW IN A PIPE

ＧＥＮＥＲＡＬＳＥＣＯＮＤＯＲＤＥＲＦＬＵＩＤＦＬＯＷＩＮＡＰＩＰＥＨｅＧｕａｎｇｙｕ（何光渝）（ＤｅｐａｒｔｍｅｎｔｏｆＰｅｔｒｏｌｅｕｍＥｎｇｉｎｅｅｒｉｎｇ，Ｘｉ＇ａｎＰｅｔｒｏｌｅｕｍＩｎｓｔｉｔｕｔｅ，Ｘｉ＇ａｎ７１００６１，Ｐ．Ｒ．Ｃｈｉ...     

分数阶导数粘弹性模型的有限元法

在分析分数阶导数三元件模型理论的基础上,把分数阶导数三元件模型引入有限元模型中,推导出具有分数阶导数三元件本构关系的粘弹性结构动力学有限元格式。同时,应用分数阶导数型粘弹性结构动力学方程的数值算法求解了该有限元格式的数值解。并以二维沥青路面结构为例进行了路面动态粘弹性响应分析。算例分析表明,该方法能够正确有效地进行路面动态粘弹性分析。     

Nonlinear flutter of viscoelastic rectangular plates and cylindrical panels of a composite with a concentrated masses

The problem of flutter of viscoelastic rectangular plates and cylindrical panels with concentrated masses is studied in a geometrically nonlinear formulation. In the equation of motion of the plate and panel, the effect of concentrated masses is accounted for using the δ-Dirac function. The problem is reduced to a system of nonlinear ordinary integrodifferential equations by using the Bubnov-Galerkin method. The resulting system with a weakly singular Koltunov-Rzhanitsyn kernel is solved by employing a numerical method based on quadrature formulas. The behavior of viscoelastic rectangular plates and cylindrical panels is studied and the critical flow velocities are determined for real composite materials over wide ranges of physicomechanical and geometrical parameters.     

General dynamic equation and dynamical characteristics of viscoelastic Timoshenko beams

In this paper, a governing differential equation of viscoelastic Timoshenko beam including both extension and shear viscosity is developed in the time domain by direct method. To measure the complex moduli and three parameters of standard linear solid, the forced vibration technique of beam is successfully used for PCL and PMMA specimens. The dynamical characteristics of viscoelastic Timoshenko beams, especially the damping properties, are derived from a considerable number of numerical computations. The analyses show that the viscosity of materials has great influence on dynamical characteristics of structures, especially on damping, and the standard linear solid model is the better one for describing the dynamic behavior of high viscous materials.     

Wave propagation in a fractional viscoelastic Andrade medium: Diffusive approximation and numerical modeling

This study focuses on the numerical modeling of wave propagation in fractionally-dissipative media. These viscoelastic models are such that the attenuation is frequency-dependent and follows a power law with non-integer exponent within certain frequency regimes. As a prototypical example, the Andrade model is chosen for its simplicity and its satisfactory fits of experimental flow laws in rocks and metals. The corresponding constitutive equation features a fractional derivative in time, a non-local-in-time term that can be expressed as a convolution product whose direct implementation bears substantial memory cost. To circumvent this limitation, a diffusive representation approach is deployed, replacing the convolution product by an integral of a function satisfying a local time-domain ordinary differential equation. An associated quadrature formula yields a local-in-time system of partial differential equations, which is then proven to be well-posed. The properties of the resulting model are also compared to those of the Andrade model. The quadrature scheme associated with the diffusive approximation, and constructed either from a classical polynomial approach or from a constrained optimization method, is investigated. Finally, the benefits of using the latter approach are highlighted as it allows to minimize the discrepancy with the original model. Wave propagation simulations in homogeneous domains are performed within a split formulation framework that yields an optimal stability condition and which features a joint fourth-order time-marching scheme coupled with an exact integration step. A set of numerical experiments is presented to assess the overall approach. Therefore, in this study, the diffusive approximation is demonstrated to provide an efficient framework for the theoretical and numerical investigations of the wave propagation problem associated with the fractional viscoelastic medium considered.     

Dynamical behaviors of nonlinear viscoelastic thick plates with damage

Based on the deformation hypothesis of Timoshenko's plates and the Boltzmann's superposition principles for linear viscoelastic materials, the nonlinear equations governing the dynamical behavior of Timoshenko's viscoelastic thick plates with damage are presented. The Galerkin method is applied to simplify the set of equations. The numerical methods in nonlinear dynamics are used to solve the simplified systems. It could be seen that there are plenty of dynamical properties for dynamical systems formed by this kind of viscoelastic thick plate with damage under a transverse harmonic load. The influences of load, geometry and material parameters on the dynamical behavior of the nonlinear system are investigated in detail. At the same time, the effect of damage on the dynamical behavior of plate is also discussed.     

黏弹性功能梯度材料裂纹问题的有限元方法

针对组分材料体积含量任意分布的黏弹性功能梯度材料裂纹问题建立有限元分析途径. 通过Laplace变换,将黏弹性问题转化到象空间中求解,基于反映材料非均匀的梯度单元和裂纹尖端奇异特性的奇异单元计算象空间中的位移、应力和应变场,应用虚拟裂纹闭合方法得到应变能释放率,分别由应力和应变能释放率确定应力强度因子. 给出这些断裂参量在物理空间和象空间之间的对应关系,由数值逆变换求出其在物理空间的相应值. 文中分析两端均匀受拉的黏弹性边裂纹板条,首先针对松弛模量表示为空间函数和时间函数乘积的特殊梯度材料进行计算,结合对应原理验证方法的有效性. 然后分析组分材料体积含量具有任意梯度分布的情形,由Mori-Tanaka方法预测象空间中的等效松弛模量. 计算结果表明,蠕变加载条件下,应变能释放率随时间增加,其增大程度与黏弹性组分材料体积含量相关. 由于梯度材料的非均匀黏弹性性质,产生应力重新分布,导致应力强度因子随时间变化,其变化范围与组分材料的体积含量分布方式有关.     

Nonlinear vibrations of viscoelastic cylindrical shells taking into account shear deformation and rotatory inertia

The vibration problem of a viscoelastic cylindrical shell is studied in a geometrically nonlinear formulation using the refined Timoshenko theory. The problem is solved by the Bubnov–Galerkin procedure combined with a numerical method based on quadrature formulas. The choice of relaxation kernels is substantiated for solving dynamic problems of viscoelastic systems. The numerical convergence of the Bubnov–Galerkin procedure is examined. The effect of viscoelastic properties of the material on the response of the cylindrical shell is discussed. The results obtained by various theories are compared.     

Approximate and numerical analysis of nonlinear forced vibration of axially moving viscoelastic beams

Steady-state periodical response is investigated for an axially moving viscoelastic beam with hybrid supports via approximate analysis with numerical confirmation. It is assumed that the excitation is spatially uniform and temporally harmonic. The transverse motion of axially moving beams is governed by a nonlinear partial-differential equation and a nonlinear integro-partial-differential equation. The material time derivative is used in the viscoelastic constitutive relation. The method of multiple scales is applied to the governing equations to investigate primary resonances under general boundary conditions. It is demonstrated that the mode uninvolved in the resonance has no effect on the steady-state response. Numerical examples are presented to demonstrate the effects of the boundary constraint stiffness on the amplitude and the stability of the steady-state response. The results derived for two governing equations are qualitatively the same,but quantitatively different. The differential quadrature schemes are developed to verify those results via the method of multiple scales.