Accommodating energy price volatility in life cycle cost analysis of asphalt pavements

Energy price is related to more than half of the total life cycle cost of asphalt pavements. Furthermore, the fluctuation related to price of energy has been much higher than the general inflation and interest rate. This makes the energy price inflation an important variable that should be addressed when performing life cycle cost (LCC) studies re- garding asphalt pavements. The present value of future costs is highly sensitive to the selected discount rate. Therefore, the choice of the discount rate is the most critical element in LCC analysis during the life time of a project. The objective of the paper is to present a discount rate for asphalt pavement projects as a function of interest rate, general inflation and energy price inflation. The discount rate is defined based on the portion of the energy related costs during the life time of the pavement. Consequently, it can reflect the financial risks related to the energy price in asphalt pavement projects. It is suggested that a discount rate sensitivity analysis for asphalt pavements in Sweden should range between –20 and 30%.

Crack initiation in asphalt mixtures under external compressive loads

Crack initiation was studied for asphalt mixtures under external compressive loads. High tensile localized stresses e direction of the external loads. A quantitative crack initiation criterion the edges of compressed air voids lead to the growth of wing cracks in thon was derived using pseudostrain energy balance principle. Bond energy is determined and it increases with aging and loading rate while decreases with temperature. Cohesive and adhesive cracking occur simultaneously and a method was proposed to determine the individual percentage. The crack initiation criterion is simplified and validated through comparing the predicted and measured compressive strength of the asphalt mixtures.

Fracture and viscoelastic properties of asphalt binders during fatigue and rest periods

Viscoelastic asphalt binder plays an important role in bonding individual aggregate particles and contributes to the durability and stability of asphalt pavement. When asphalt binder is subjected to cyclic loading, deformation and fracture may develop simultaneously within it, leading to the deterioration of material properties and eventually fatigue failure. Research has found that some degree of recovery may develop if rest periods are applied after fatigue deterioration. However, it is not clear whether such recovery is caused by fracture healing, viscoelastic recovery, or both. This paper presents an analysis to differentiate the contributions of fracture healing and viscoelastic recovery to the asphalt binder during rest periods. It also evaluates the damage caused by deformation and fracture during a fatigue process. It is found that viscoelastic recovery plays an important role in the instant increase in the dynamic shear modulus at the beginning of the rest period. The effect of fracture healing on dynamic shear modulus recovery is more dominant in the long term. A healing index is developed based only on the fracture healing of asphalt binder, excluding the effect of viscoelastic recovery. It can be used to evaluate the true healing properties of different asphalt binders. Copyright © 2014 by ASTM International.

Anisotropic characterization of crack growth in the tertiary flow of asphalt mixtures in compression

Asphalt mixtures exhibit primary, secondary, and tertiary stages in sequence during a rutting deterioration. Many field asphalt pavements are still in service even when the asphalt layer is in the tertiary stage, and rehabilitation is not performed until a significant amount of rutting accompanied by numerous macrocracks is observed. The objective of this study was to provide a mechanistic method to model the anisotropic cracking of the asphalt mixtures in compression during the tertiary stage of rutting. Laboratory tests including nondestructive and destructive tests were performed to obtain the viscoelastic and viscofracture properties of the asphalt mixtures. Each of the measured axial and radial total strains in the destructive tests were decomposed into elastic, plastic, viscoelastic, viscoplastic, and viscofracture strains using the pseudostrain method in an extended elastic-viscoelastic correspondence principle. The viscofracture strains are caused by the crack growth, which is primarily signaled by the increase of phase angle in the tertiary flow. The viscofracture properties are characterized using the anisotropic damage densities (i.e., the ratio of the lost area caused by cracks to the original total area in orthogonal directions). Using the decomposed axial and radial viscofracture strains, the axial and radial damage densities were determined by using a dissipated pseudostrain energy balance principle and a geometric analysis of the cracks, respectively. Anisotropic pseudo J-integral Paris' laws in terms of damage densities were used to characterize the evolution of the cracks in compression. The material constants in the Paris' law are determined and found to be highly correlated. These tests, analysis, and modeling were performed on different asphalt mixtures with two binders, two air void contents, and three aging periods. Consistent results were obtained; for instance, a stiffer asphalt mixture is demonstrated to have a higher modulus, a lower phase angle, a greater flow number, and a larger n1 value (exponent of Paris' law). The calculation of the orientation of cracks demonstrates that the asphalt mixture with 4% air voids has a brittle fracture and a splitting crack mode, whereas the asphalt mixture with 7% air voids tends to have a ductile fracture and a diagonal sliding crack mode. Cracks of the asphalt mixtures in compression are inclined to propagate along the direction of the external compressive load. © 2014 American Society of Civil Engineers.

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.

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.