Mechanistic modeling of fracture in asphalt mixtures under compressive loading

When an asphalt mixture is subjected to a destructive compressive load, it experiences a sequence of three deformation stages, as follows: the (1) primary, (2) secondary, and (3) tertiary stages. Most literature research focuses on plastic deformation in the primary and secondary stages, such as prediction of the flow number, which is in fact the initiation of the tertiary stage. However, little research effort has been reported on the mechanistic modeling of the damage that occurs in the tertiary stage. The main objective of this paper is to provide a mechanistic characterizing method for the damage modeling of asphalt mixtures in the tertiary stage. The preliminary study conducted by the writers illustrates that deformation during the tertiary flow of the asphalt mixtures is principally caused by the formation and propagation of cracks, which was signaled by the increase of the phase angle in the tertiary phase. The strain caused by the growth of cracks is the viscofracture strain, which can be obtained by conducting the strain decomposition of the measured total strain in the destructive compressive test. The viscofracture strain is employed in the research reported in this paper to mechanistically characterize the time-dependent fracture (viscofracture) of asphalt mixtures in compression. By using the dissipated pseudostrain energy-balance principle, the damage density and true stress are determined and both are demonstrated to increase with load cycles in the tertiary stage. The increased true stress yields extra viscoplastic strain, which is the reason why the permanent deformation is accelerated by the occurrence of cracks. To characterize the evolution of the viscofracture in the asphalt mixtures in compression, a pseudo J-integral Paris' law in terms of damage density is proposed and the material constants in the Paris' law are determined, which can be employed to predict the fracture of asphalt mixtures in compression. © 2013 American Society of Civil Engineers.

Characterization of viscoplastic yielding of asphalt concrete

A temperature and strain rate dependent yield surface model was proposed to characterize the viscoplastic yielding of asphalt concrete. Laboratory tests were conducted on specimens that have two binders, two air void contents, and three aging periods. Strain decomposition was performed to obtain viscoplastic strain and stress-pseudostrain curves were constructed to determine the model parameters accurately and efficiently. Results indicate that a stiffer asphalt concrete has greater cohesion and strain hardening amplitude, both of which decline as temperature increases or strain rate decreases. The temperature and strain rate factors of the yield surface can be accurately determined solely by the peak stress of the strength tests. © 2013 Elsevier Ltd. All rights reserved.

Characterizing permanent deformation and fracture of asphalt mixtures by using compressive dynamic modulus tests

Permanent deformation and fracture may develop simultaneously when an asphalt mixture is subjected to a compressive load. The objective of this research is to separate viscoplasticity and viscofracture from viscoelasticity so that the permanent deformation and fracture of the asphalt mixtures can be individually and accurately characterized without the influence of viscoelasticity. The undamaged properties of 16 asphalt mixtures that have two binder types, two air void contents, and two aging conditions are first obtained by conducting nondestructive creep tests and nondestructive dynamic modulus tests. Testing results are analyzed by using the linear viscoelastic theory in which the creep compliance and the relaxation modulus are modeled by the Prony model. The dynamic modulus and phase angle of the undamaged asphalt mixtures remained constant with the load cycles. The undamaged asphalt mixtures are then used to perform the destructive dynamic modulus tests in which the dynamic modulus and phase angle of the damaged asphalt mixtures vary with load cycles. This indicates plastic evolution and crack propagation. The growth of cracks is signaled principally by the increase of the phase angle, which occurs only in the tertiary stage. The measured total strain is successfully decomposed into elastic strain, viscoelastic strain, plastic strain, viscoplastic strain, and viscofracture strain by employing the pseudostrain concept and the extended elastic-viscoelastic correspondence principle. The separated viscoplastic strain uses a predictive model to characterize the permanent deformation. The separated viscofracture strain uses a fracture strain model to characterize the fracture of the asphalt mixtures in which the flow number is determined and a crack speed index is proposed. Comparisons of the 16 samples show that aged asphalt mixtures with a low air void content have a better performance, resisting permanent deformation and fracture. © 2012 American Society of Civil Engineers.

Microstructure-based inherent anisotropy of asphalt mixtures

Asphalt mixtures have been demonstrated to be anisotropic materials in both laboratory and field tests. The anisotropy of asphalt mixtures consists of inherent anisotropy and stress-induced anisotropy. In previous work, the inherent anisotropy of asphalt mixtures was quantified by using only the inclination angles of the coarse aggregate particles in the asphalt mixtures. However, the inclination of fine aggregates also has a contribution to the inherent anisotropy. Moreover, the contribution to the inherent anisotropy of each aggregate may not be the same as in the previous work but will depend on the size, orientation, and sphericity of the aggregate particle. This paper quantifies the internal microstructure of the aggregates in asphalt mixtures by using an aggregate-related geometric parameter, the vector magnitude. The original formulation of the vector magnitude, which addresses only the orientation of coarse aggregates, is modified to account for not only the coarse aggregate orientation, but also the size, orientation, and sphericity of coarse and fine aggregates. This formulation is applied to cylindrical lab-mixed lab-compacted asphalt mixture specimens varying in asphalt binder type, air void content, and aging period. The vertical modulus and the horizontal modulus are also measured by using nondestructive tests. A relationship between the modified vector magnitude and the modulus ratio of the vertical modulus to the horizontal modulus is developed to quantify the influence of the inherent microstructure of the aggregates on the anisotropy of the mixtures. The modulus ratio is found to depend solely on the aggregate characteristics including the inclination angle, size, and sphericity, and it is independent of the asphalt binder type, air void content, and aging period. The inclination angle, itself, proves to be insufficient to quantify the inherent anisotropy of the asphalt mixtures. © 2011 American Society of Civil Engineers.

Anisotropic viscoelastic properties of undamaged asphalt mixtures

A test protocol and a data analysis method are developed in this paper on the basis of linear viscoelastic theory to characterize the anisotropic viscoelastic properties of undamaged asphalt mixtures. The test protocol includes three nondestructive tests: (1) uniaxial compressive creep test, (2) indirect tensile creep test, and (3) the uniaxial tensile creep test. All three tests are conducted on asphalt mixture specimens at three temperatures (10, 20, and 30°C) to determine the tensile and compressive properties at each temperature and then to construct the master curve of each property. The determined properties include magnitude and phase angle of the compressive complex modulus in the vertical direction, magnitude and phase angle of the tensile complex modulus, and the magnitude and phase angle of the compressive complex modulus in the horizontal plane. The test results indicate that all tested asphalt mixtures have significantly different tensile properties from compressive properties. The peak value of the master curve of the tensile complex modulus phase angle is within a range from 65 to 85°, whereas the peak value of the compressive moduli phase angle in both directions ranges from 35 to 55°. In addition, the undamaged asphalt mixtures exhibit distinctively anisotropic properties in compression. The magnitude of the compressive modulus in the vertical direction is approximately 1.2 to ̃2 times of the magnitude of the compressive modulus in the horizontal plane. Dynamic modulus tests are performed to verify the results of the proposed test protocol. The test results from the proposed test protocol match well with those from the dynamic tests. © 2012 American Society of Civil Engineers.

A new creep test method for asphalt mixtures

A new creep test, Partial Triaxial Test (PTT), was developed to study the permanent deformation properties of asphalt mixtures. The PTT used two duplicate platens whose diameters were smaller than the diameter of the cylindrical asphalt mixtures specimen. One base platen was centrally placed under the specimen and another loading platen was centrally placed on the top surface of the specimen. Then the compressive repeated load was applied on the loading platen and the vertical deformation of the asphalt mixture was recorded in the PTTs. Triaxial repeated load permanent deformation tests (TRT) and PTTs were respectively conducted on AC20 and SMA13 asphalt mixtures at 40°C and 60°C so as to provide the parameters of the creep constitutive relations in the ABAQUS finite element models (FEMs) which were built to simulate the laboratory wheel tracking tests. The real laboratory wheel tracking tests were also conducted on AC20 and SMA13 asphalt mixtures at 40°C and 60°C. Then the calculated rutting depth from the FEMs were compared with the measured rutting depth of the laboratory wheeling tracking tests. Results indicated that PTT was able to characterize the permanent deformation of the asphalt mixtures in laboratory. The rutting depth calculated using the parameters estimated from PTTs' results was closer to and showed better matches with the measured rutting than the rutting depth calculated using the parameters estimated from TRTs' results. Main reason was that PTT could better simulate the changing confinement conditions of asphalt mixtures in the laboratory wheeling tracking tests than the TRT.

A generalized Drucker–Prager viscoplastic yield surface model for asphalt concrete

A generalized Drucker–Prager (GD–P) viscoplastic yield surface model was developed and validated for asphalt concrete. The GD–P model was formulated based on fabric tensor modified stresses to consider the material inherent anisotropy. A smooth and convex octahedral yield surface function was developed in the GD–P model to characterize the full range of the internal friction angles from 0° to 90°. In contrast, the existing Extended Drucker–Prager (ED–P) was demonstrated to be applicable only for a material that has an internal friction angle less than 22°. Laboratory tests were performed to evaluate the anisotropic effect and to validate the GD–P model. Results indicated that (1) the yield stresses of an isotropic yield surface model are greater in compression and less in extension than that of an anisotropic model, which can result in an under-prediction of the viscoplastic deformation; and (2) the yield stresses predicted by the GD–P model matched well with the experimental results of the octahedral shear strength tests at different normal and confining stresses. By contrast, the ED–P model over-predicted the octahedral yield stresses, which can lead to an under-prediction of the permanent deformation. In summary, the rutting depth of an asphalt pavement would be underestimated without considering anisotropy and convexity of the yield surface for asphalt concrete. The proposed GD–P model was demonstrated to be capable of overcoming these limitations of the existing yield surface models for the asphalt concrete.

High performance liquid level monitoring system based on polymer fiber Bragg gratings embedded in silicone rubber diaphragm

Liquid-level sensing technologies have attracted great prominence, because such measurements are essential to industrial applications, such as fuel storage, flood warning and in the biochemical industry. Traditional liquid level sensors are based on electromechanical techniques; however they suffer from intrinsic safety concerns in explosive environments. In recent years, given that optical fiber sensors have lots of well-established advantages such as high accuracy, costeffectiveness, compact size, and ease of multiplexing, several optical fiber liquid level sensors have been investigated which are based on different operating principles such as side-polishing the cladding and a portion of core, using a spiral side-emitting optical fiber or using silica fiber gratings. The present work proposes a novel and highly sensitive liquid level sensor making use of polymer optical fiber Bragg gratings (POFBGs). The key elements of the system are a set of POFBGs embedded in silicone rubber diaphragms. This is a new development building on the idea of determining liquid level by measuring the pressure at the bottom of a liquid container, however it has a number of critical advantages. The system features several FBG-based pressure sensors as described above placed at different depths. Any sensor above the surface of the liquid will read the same ambient pressure. Sensors below the surface of the liquid will read pressures that increase linearly with depth. The position of the liquid surface can therefore be approximately identified as lying between the first sensor to read an above-ambient pressure and the next higher sensor. This level of precision would not in general be sufficient for most liquid level monitoring applications; however a much more precise determination of liquid level can be made by linear regression to the pressure readings from the sub-surface sensors. There are numerous advantages to this multi-sensor approach. First, the use of linear regression using multiple sensors is inherently more accurate than using a single pressure reading to estimate depth. Second, common mode temperature induced wavelength shifts in the individual sensors are automatically compensated. Thirdly, temperature induced changes in the sensor pressure sensitivity are also compensated. Fourthly, the approach provides the possibility to detect and compensate for malfunctioning sensors. Finally, the system is immune to changes in the density of the monitored fluid and even to changes in the effective force of gravity, as might be obtained in an aerospace application. The performance of an individual sensor was characterized and displays a sensitivity (54 pm/cm), enhanced by more than a factor of 2 when compared to a sensor head configuration based on a silica FBG published in the literature, resulting from the much lower elastic modulus of POF. Furthermore, the temperature/humidity behavior and measurement resolution were also studied in detail. The proposed configuration also displays a highly linear response, high resolution and good repeatability. The results suggest the new configuration can be a useful tool in many different applications, such as aircraft fuel monitoring, and biochemical and environmental sensing, where accuracy and stability are fundamental. © (2015) COPYRIGHT Society of Photo-Optical Instrumentation Engineers (SPIE). Downloading of the abstract is permitted for personal use only.

Weak form equation–based finite-element modeling of viscoelastic asphalt mixtures

The objective of this study is to demonstrate using weak form partial differential equation (PDE) method for a finite-element (FE) modeling of a new constitutive relation without the need of user subroutine programming. The viscoelastic asphalt mixtures were modeled by the weak form PDE-based FE method as the examples in the paper. A solid-like generalized Maxwell model was used to represent the deforming mechanism of a viscoelastic material, the constitutive relations of which were derived and implemented in the weak form PDE module of Comsol Multiphysics, a commercial FE program. The weak form PDE modeling of viscoelasticity was verified by comparing Comsol and Abaqus simulations, which employed the same loading configurations and material property inputs in virtual laboratory test simulations. Both produced identical results in terms of axial and radial strain responses. The weak form PDE modeling of viscoelasticity was further validated by comparing the weak form PDE predictions with real laboratory test results of six types of asphalt mixtures with two air void contents and three aging periods. The viscoelastic material properties such as the coefficients of a Prony series model for the relaxation modulus were obtained by converting from the master curves of dynamic modulus and phase angle. Strain responses of compressive creep tests at three temperatures and cyclic load tests were predicted using the weak form PDE modeling and found to be comparable with the measurements of the real laboratory tests. It was demonstrated that the weak form PDE-based FE modeling can serve as an efficient method to implement new constitutive models and can free engineers from user subroutine programming.

Energy-based damage and fracture framework for viscoelastic asphalt concrete

A framework based on the continuum damage mechanics and thermodynamics of irreversible processes using internal state variables is used to characterize the distributed damage in viscoelastic asphalt materials in the form of micro-crack initiation and accumulation. At low temperatures and high deformation rates, micro-cracking is considered as the source of nonlinearity and thus the cause of deviation from linear viscoelastic response. Using a non-associated damage evolution law, the proposed model shows the ability to describe the temperature-dependent processes of micro-crack initiation, evolution and macro-crack formation with good comparison to the material response in the Superpave indirect tensile (IDT) strength test.

Towards asphalt mixture morphology evaluation with the virtual specimen approach

The morphology of asphalt mixture can be defined as a set of parameters describing the geometrical characteristics of its constituent materials, their relative proportions as well as spatial arrangement in the mixture. The present study is carried out to investigate the effect of the morphology on its meso- and macro-mechanical response. An analysis approach is used for the meso-structural characterisation based on the X-ray computed tomography (CT) data. Image processing techniques are used to systematically vary the internal structure to obtain different morphology structures. A morphology framework is used to characterise the average mastic coating thickness around the main load carrying structure in the structures. The uniaxial tension simulation shows that the mixtures with the lowest coating thickness exhibit better inter-particle interaction with more continuous load distribution chains between adjacent aggregate particles, less stress concentrations and less strain localisation in the mastic phase.

Effect of mixing sequence on the workability and performance of asphalt mixtures

In Sweden, during recent years, a new type of mixing protocol has been applied, in which the order of mixing is changed from the conventional method. Improved workability and diminished mixing and compaction energy needs have been important drivers for this. Considering that it is the mastic phase, which is modified by changing the mixing order, it provides an interesting case study for explaining the mechanisms of workability in connection with the mastic phase. To do so, an analytical viscosity framework was combined with a mixture morphology framework to upscale to the mixing level and tribology principles to explain the interaction between the mastic and the aggregates. From the mastic viscosity protocol, it was found that the mixing order significantly affects the resulting mastic viscosity. To analyse the effect of this on the workability and resulting mixture performance, X-ray computed tomography was used to analyse mixtures produced by the two different mixing sequences. Mechanical testing was utilised to determine the long-term mechanical performance. In this part of the study, mastic viscosity as a function of particle concentration and distribution was directly coupled to improved mixture workability and enhanced long-term performance.

Using overlay test to evaluate fracture properties of field-aged asphalt concrete

Field material testing provides firsthand information on pavement conditions which are most helpful in evaluating performance and identifying preventive maintenance or overlay strategies. High variability of field asphalt concrete due to construction raises the demand for accuracy of the test. Accordingly, the objective of this study is to propose a reliable and repeatable methodology to evaluate the fracture properties of field-aged asphalt concrete using the overlay test (OT). The OT is selected because of its efficiency and feasibility for asphalt field cores with diverse dimensions. The fracture properties refer to the Paris’ law parameters based on the pseudo J-integral (A and n) because of the sound physical significance of the pseudo J-integral with respect to characterizing the cracking process. In order to determine A and n, a two-step OT protocol is designed to characterize the undamaged and damaged behaviors of asphalt field cores. To ensure the accuracy of determined undamaged and fracture properties, a new analysis method is then developed for data processing, which combines the finite element simulations and mechanical analysis of viscoelastic force equilibrium and evolution of pseudo displacement work in the OT specimen. Finally, theoretical equations are derived to calculate A and n directly from the OT test data. The accuracy of the determined fracture properties is verified. The proposed methodology is applied to a total of 27 asphalt field cores obtained from a field project in Texas, including the control Hot Mix Asphalt (HMA) and two types of warm mix asphalt (WMA). The results demonstrate a high linear correlation between n and −log A for all the tested field cores. Investigations of the effect of field aging on the fracture properties confirm that n is a good indicator to quantify the cracking resistance of asphalt concrete. It is also indicated that summer climatic condition clearly accelerates the rate of aging. The impact of the WMA technologies on fracture properties of asphalt concrete is visualized by comparing the n-values. It shows that the Evotherm WMA technology slightly improves the cracking resistance, while the foaming WMA technology provides the comparable fracture properties with the HMA. After 15 months aging in the field, the cracking resistance does not exhibit significant difference between HMA and WMAs, which is confirmed by the observations of field distresses.

Implementation of pseudo J-integral based Paris’ law for fatigue cracking in asphalt mixtures and pavements

Pavement analysis and design for fatigue cracking involves a number of practical problems like material assessment/screening and performance prediction. A mechanics-aided method can answer these questions with satisfactory accuracy in a convenient way when it is appropriately implemented. This paper presents two techniques to implement the pseudo J-integral based Paris’ law to evaluate and predict fatigue cracking in asphalt mixtures and pavements. The first technique, quasi-elastic simulation, provides a rational and appropriate reference modulus for the pseudo analysis (i.e., viscoelastic to elastic conversion) by making use of the widely used material property: dynamic modulus. The physical significance of the quasi-elastic simulation is clarified. Introduction of this technique facilitates the implementation of the fracture mechanics models as well as continuum damage mechanics models to characterize fatigue cracking in asphalt pavements. The second technique about modeling fracture coefficients of the pseudo J-integral based Paris’ law simplifies the prediction of fatigue cracking without performing fatigue tests. The developed prediction models for the fracture coefficients rely on readily available mixture design properties that directly affect the fatigue performance, including the relaxation modulus, air void content, asphalt binder content, and aggregate gradation. Sufficient data are collected to develop such prediction models and the R2 values are around 0.9. The presented case studies serve as examples to illustrate how the pseudo J-integral based Paris’ law predicts fatigue resistance of asphalt mixtures and assesses fatigue performance of asphalt pavements. Future applications include the estimation of fatigue life of asphalt mixtures/pavements through a distinct criterion that defines fatigue failure by its physical significance.