基于细观结构有限元模型的环氧沥青混合料间接拉伸试验研究

运用有限元方法建立环氧沥青混合料细观结构模型,对其间接拉伸试验(IDT)进行数值模拟.首先借助图像处理技术得到由集料和沥青砂浆组成的环氧沥青混合料二相细观结构,并通过蠕变试验获取沥青砂浆常温下的黏弹性材料参数,最后结合有限元手段建立包含集料、砂浆等在内的混合料细观结构有限元模型.数值模拟结果表明,有限元计算的混合料劲度模量与实际IDT试验结果吻合较好,通过改变加载方向、加载速率等参数,发现对混合料细观结构的劲度模量以及局部点位应力均造成一定影响,分析主要原因可能是由沥青混合料的内部结构分布不均匀性以及沥青砂浆的黏弹性特点所造成.研究成果可为微观有限元方法进一步推广应用于不同条件下沥青混合料微观力学响应仿真提供理论依据.     

Viscoelastic Model for Discrete Element Simulation of Asphalt Mixtures

This paper presents a viscoelastic model of asphalt mixtures with the discrete element method, where the viscoelastic behaviors of asphalt mastics (fine aggregates, fines, and asphalt binder) are represented by a Burger’s model. Aggregates are simulated with irregular shape particles consisting of balls bonded together by elastic contact models, and the interplaces between aggregates are filled with balls bonded with viscoelastic Burger’s model to represent asphalt mastic. Digital samples were prepared with the image analysis technique. The micromechanical model was developed with four constitutive laws to represent the interactions at contacts of discrete elements (balls) within an aggregate, within mastic, between an aggregate and mastic, and between two adjacent aggregates. Each of these constitutive laws consists of three parts: a stiffness model, a slip model, and a bonding model in order to provide a relationship between the contact force and relative displacement and also in order to describe slipping and tensile strength at a particular contact. The relationship between the microscale model input and macroscale material properties was derived, and an iterative procedure was developed to fit the dynamic modulus test data of asphalt mastic with Burger’s model. The favorable agreement between the discrete element prediction and the lab results on dynamic moduli and phase angles indicates that the viscoelastic discrete element model developed in this study is very capable of simulating constitutive behavior of asphalt mixtures.     

Microstructural Finite-Element Analysis of Influence of Localized Strain Distribution on Asphalt Mix Properties

Hot mix asphalt (HMA) is a composite material that consists of mineral aggregates, asphalt binders, and air voids. A finite element model of the HMA microstructure is developed to study the influence of localized strain distribution on the HMA mechanical response. Image analysis techniques are used to capture the HMA microstructure. Due to limitations on the image resolution, the microstructure is divided into two phases: aggregates larger than 0.3 mm and mastic (binder and aggregates smaller than 0.3 mm). A viscoelastic constitutive relationship is used to represent the mastic phase of the HMA microstructure. The mastic viscoelastic properties are obtained from the results of testing asphalt binders in a dynamic shear rheometer and microstructure analysis of idealized mastic. A step-wise finite element procedure is employed in order to account for the influence of the localized high strains on the mastic viscoelastic properties, and HMA mechanical response. The mastic and binder elements of the microstructure are shown to exhibit high strain values within the nonlinear viscoelastic range. The HMA viscoelastic properties are calculated at different strain levels, and the results are compared with experimental data obtained from the frequency sweep shear test.     

Three-Dimensional Discrete Element Models for Asphalt Mixtures

The main objective of this paper is to develop three-dimensional (3D) microstructure-based discrete element models of asphalt mixtures to study the dynamic modulus from the stress-strain response under compressive loads. The 3D microstructure of the asphalt mixture was obtained from a number of two-dimensional (2D) images. In the 2D discrete element model, the aggregate and mastic were simulated with the captured aggregate and mastic images. The 3D models were reconstructed with a number of 2D models. This stress-strain response of the 3D model was computed under the loading cycles. The stress-strain response was used to predict the asphalt mixture’s stiffness (modulus) by using the aggregate and mastic stiffness. The moduli of the 3D models were compared with the experimental measurements. It was found that the 3D discrete element models were able to predict the mixture moduli across a range of temperatures and loading frequencies. The 3D model prediction was found to be better than that of the 2D model. In addition, the effects of different air void percentages and aggregate moduli to the mixture moduli were investigated and discussed.     

Visualization and Simulation of Asphalt Concrete with Randomly Generated Three-Dimensional Models

The objective of this study is to visualize and simulate microscale properties of asphalt concrete with three-dimensional discrete element models under mechanical loading. The microstructure of the asphalt concrete sample was composed of three ingredients: coarse aggregates, sand mastic (a combination of fines, fine aggregates, and asphalt binder), and air voids. Coarse aggregates were represented with the irregular polyhedron particles which were randomly created with an algorithm developed for this study. The gaps among the irregular particles were filled with air voids and discrete elements of sand mastic. The mechanical behaviors of the three ingredients were simulated with specific constitutive models at different contacts of discrete elements. Based on the geometric and mechanical models, visualization and simulation of asphalt mixtures were conducted in this study.     

Prediction of Damage Behaviors in Asphalt Materials Using a Micromechanical Finite-Element Model and Image Analysis

A study of the micromechanical damage behavior of asphalt concrete is presented. Asphalt concrete is composed of aggregates, mastic cement, and air voids, and its load carrying behavior is strongly related to the local microstructural load transfer between aggregate particles. Numerical simulation of this micromechanical behavior was accomplished by using a finite-element model that incorporated the mechanical load-carrying response between aggregates. The finite-element scheme used a network of special frame elements each with a stiffness matrix developed from an approximate elasticity solution of the stress and displacement field in a cementation layer between particle pairs. Continuum damage mechanics was then incorporated within this solution, leading to the construction of a microdamage model capable of predicting typical global inelastic behavior found in asphalt materials. Using image processing and aggregate fitting techniques, simulation models of indirect tension, and compression samples were generated from surface photographic data of actual laboratory specimens. Model simulation results of the overall sample behavior and evolving microfailure/fracture patterns compared favorably with experimental data collected on these samples.     

Continuum Microviscoelasticity Model for Aging Basic Creep of Early-Age Concrete

We propose a micromechanics model for aging basic creep of early-age concrete. Therefore, we formulate viscoelastic boundary value problems on two representative volume elements, one related to cement paste (composed of cement, water, hydrates, and air), and one related to concrete (composed of cement paste and aggregates). Homogenization of the “nonaging” elastic and viscoelastic properties of the material’s contituents involves the transformation of the aforementioned viscoelastic boundary value problems to the Laplace-Carson (LC) domain. There, formally elastic, classical self-consistent and Mori-Tanaka solutions are employed, leading to pointwisely defined LC-transformed tensorial creep and relaxation functions. Subsequently, the latter are back-transformed, by means of the Gaver-Wynn-Rho algorithm, into the time domain. Temporal derivatives of corresponding homogenized creep and relaxation tensors, evaluated for the current maturation state of the material (in terms of current volume fractions of cement, water, air, hydrates, and aggregates; being dependent on the hydration degree, as well as on the water-cement and aggregate-cement ratios) and for the current time period since loading of the hydrating composite material, allow for micromechanical prediction of the aging basic creep properties of early-age concrete.     

Parameter Identification Procedure for Heterogeneous Viscoelastic Composites Using Iterative Functions

This paper presents a new numerical procedure for the determination of the viscoelastic compliance properties of a matrix phase from a simple three-point bending test on a composite beam. The composite is modeled as elastic inclusions randomly dispersed throughout a viscoelastic matrix. It is also assumed that the spatial distribution of the inclusions in the composite is known or can be determined. Zevin’s method of iterative functions is proposed for the determination of the matrix properties. Following a detailed explanation of the proposed scheme, a numerical verification is performed using three-dimensional finite-element (FE) analysis simulations. The proposed scheme was applied to the experimentally obtained creep compliance of the asphalt concrete beam. The obtained viscoelastic properties of the asphalt binder matrix phase were used as input into the FE model to simulate the behavior of the composite beam. An excellent comparison between the experimental data and the predicted beam deflections was observed. This shows that the proposed method is robust and it can be implemented to solve identification problems for viscoelastic composite materials.     

Determining Hardness and Elastic Modulus of Asphalt by Nanoindentation

Nanoindentation is a relatively new technique which has been used to measure nanomechanical properties of surface layers of bulk materials and of thin films. In this study, micromechanical properties such as hardness and Young’s modulus of asphalt binders and asphalt concrete are determined by nanoindentation experiments. Indentation tests are conducted on a base binder and two polymer-modified performance grade (PG) binders such as PG-70-22 and PG76-28. In addition, two Superpave asphalt mixes such as SP-B and SP-III are designed using these PG binders, and the corresponding mixes are compacted to prepare asphalt concrete. Aggregate, matrix (Materials Passing No. 4 sieve) and mastic (Materials Passing No. 200 sieve) phases of each asphalt concrete sample are indented using both Berkovich and Spherical indenters. In nanoindentation, an indenter penetrates into asphalt material and the load (milli-Newton) and the depth (nanometers) of indentation are recorded continuously. Indentation load versus displacement data are analyzed using Oliver and Pharr method to measure hardness and Young’s modulus. The unloading data of base binder is a straight line and therefore could not be analyzed using Oliver and Pharr’s method. However, the indentation data of the PG grade binders are successfully analyzed. Young’s modulus value is less than 3 GPa for mastic, 3 to 12 GPa for matrix, and greater than 12 GPa for aggregate studied herein. Based on the hardness data, mastic is 2 to 15 times softer than matrix materials, and matrix is 10 times softer than aggregate materials. The fact that the properties of the mastic can be measured while in the mixture, this study has great potential for realistic characterization of asphalt mixture components. In this study, spherical indenter is found to be suitable for asphalt binders based on the fact that the spherical indenter produces higher indentation depths than the Berkovich indenter. The study contributes significantly to the use of nanoindentation for transportation material characterization.     

Discrete-Element Modeling: Impacts of Aggregate Sphericity, Orientation, and Angularity on Creep Stiffness of Idealized Asphalt Mixtures

Hot-mix asphalt (HMA) contains a significant amount of mineral aggregate, approximately 95% by weight and 85% by volume. The aggregate sphericity, orientation, and angularity are very important in determining HMA mechanical behaviors. The objective of this study is to investigate the isolated effects of the aggregate sphericity index, fractured faces, and orientation angles on the creep stiffness of HMA mixtures. The discrete-element method was employed to simulate creep compliance tests on idealized HMA mixtures. Two user-defined models were used to build 102 idealized asphalt-mix digital specimens. They were the R-model and the A-model, short for a user-defined rounded aggregate model and a user-defined angular aggregate model, respectively. Of the 102 digital specimens, 84 were prepared with the R-model to investigate the effects of aggregate sphericity and orientation, whereas the remaining 18 were built with the A-model to address the effect of aggregate angularity. A viscoelastic model was used to capture the interactions within the mix specimens. It was observed that (1) as the sphericity increased, the creep stiffness of the HMA mixture increased or decreased, depending on the angles of aggregate orientation; (2) as the angle of aggregate orientation increased, the creep stiffness of the HMA mixture increased, with the rate depending on the sphericity index values; and (3) compared with the sphericity index and orientation angles, the influence of aggregate fractured faces was insignificant.     

Creep Behavior of Hot-Mix Asphalt due to Heavy Vehicular Tire Loading

The objective of this study is to characterize the creep behavior of hot-mix asphalt (HMA) at intermediate (20°C) and at high temperatures (40°C). To accomplish this objective, a nonlinear time-hardening creep model, characterized through laboratory testing, was incorporated into a three-dimensional (3D) finite-element (FE) model, which was used to calculate permanent creep strains after applying repetitive vehicular loading cycles at the pavement surface. Two different tire configurations were simulated representing a typical dual-tire assembly and a newly introduced wide-base tire (dual-tire: 275/80R22.5 and wide-base tire: 455/55R22.5). Results of the 3D FE model were successfully verified against pavement response measurements in the field at the Virginia Smart Road. While the elastic or linear viscoelastic FE model may not simulate permanent deformation or shear creep strains after repetitions of vehicular loading, a nonlinear time-hardening creep model could predict primary rutting damage in HMA and shear creep strains at the edge of the tire imprint caused by different tire configurations.     

Viscoelastic Modeling and Field Validation of Flexible Pavements

The objective of this study was to characterize hot-mix asphalt (HMA) viscoelastic properties at intermediate and high temperatures and to incorporate laboratory-determined parameters into a three-dimensional finite element (FE) model to accurately simulate pavement responses to vehicular loading at different temperatures and speeds. Results of the developed FE model were compared against field-measured pavement responses from the Virginia Smart Road. Results of this analysis indicated that the elastic theory grossly underpredicts pavement responses to vehicular loading at intermediate and high temperatures. In addition, the elastic FE model could not simulate permanent deformation or delayed recovery, a known characteristic of HMA materials. In contrast, results of the FE viscoelastic model were in better agreement with field measurements. In this case, the average error in the prediction was less than 15%. The FE model successfully simulated retardation of the response in the transverse direction and rapid relaxation of HMA in the longitudinal direction. Moreover, the developed model allowed predicting primary rutting damage at the surface and its partial recovery after load application.     

Determination of Shear and Bulk Moduli of Viscoelastic Solids from the Indirect Tension Creep Test

Because of its efficiency in analyzing complex viscoelastic problems, the finite-element (FE) analysis has been widely used to identify the time- and rate-dependent effects of viscoelastic materials on various structural conditions. When performing the FE analysis on a viscoelastic structure, most FE programs require fundamental material properties, shear and bulk moduli, of the given viscoelastic material as their input. However, the shear and bulk modulus tests are difficult to perform, so they have been commonly estimated from a single material test on the basis of the assumption that the Poisson’s ratio of viscoelastic materials is a time-independent constant. Such an assumption, however, might not be suitable because the Poisson’s ratio of the viscoelastic materials is also a function of time. Therefore, this study developed computation algorithms for determining the time-dependent Poisson’s ratio and shear and bulk moduli of asphalt mixtures, which have been well recognized as a viscoelastic material, by employing the indirect tension testing system. The shear and bulk moduli determined by the developed approach appear to be reasonable and realistic. Their applicability and reliability were also verified by comparing experimental data to the results of the FE analysis performed on the same circular specimen as that used in the indirect tension creep test.     

Micromechanical Modeling of the Viscoelastic Behavior of Asphalt Mixtures Using the Discrete-Element Method

This paper presents a methodology for analyzing the viscoelastic response of asphalt mixtures using the discrete-element method (DEM). Two unmodified (neat) and seven modified binders were mixed with the same aggregate blend in order to prepare the nine hot mix asphalt (HMA) mixtures used in this study. The HMA microstructure was captured using images of vertically cut sections of specimens. The captured grayscale images were processed into black and white images representing the mastic and the aggregate phases, respectively. These microstructure images were used to represent the DEM model geometry. Rheological data for the nine binders were obtained using the dynamic shear rheometer. These data were used to estimate the parameters of the viscoelastic contact models that define the interaction among the mix constituents. The DEM models were subjected to sinusoidal loads similar to those applied in the simple performance test (SPT). The DEM model predictions compared favorably with the SPT measurements. However, the simulation results tended to overpredict the dynamic modulus, E*, for mixtures made with neat binders and underpredict E* for those that consisted of modified binders. The DEM models gave mix phase angles, ?mix, higher than the experimental measurements.     

Development of an Elastoviscoplastic Microstructural-Based Continuum Model to Predict Permanent Deformation in Hot Mix Asphalt

Permanent deformation in hot mix asphalt is caused by a combination of densification (decrease in volume and hence increase in density) and shear deformation. The primary objective of this paper is to develop an elastoviscoplastic model that accounts for the influence of important microstructure properties such as anisotropy and damage on permanent deformation. The model incorporates a yield surface based on the Drucker-Prager function that is modified to capture the influence of stress state on the material response. Also, parameters that reflect the directional distribution of aggregates and damage density in the microstructure are included in this yield surface model. The elastoviscoplastic model is converted into a numerical formulation and is implemented in finite element (FE). The FE model is used in this study to simulate experimental measurements under different confining pressures and strain rates.     

Accelerated Discrete-Element Modeling of Asphalt-Based Materials with the Frequency-Temperature Superposition Principle

This paper presents a methodology to reduce the computation time for discrete-element (DE) modeling of asphalt-based materials, based on the frequency-temperature superposition principle. Laboratory tests on the dynamic modulus of asphalt sand mastics and asphalt mixtures were conducted at temperatures of -5, 4, 13, and 21°C and frequencies of 1, 5, 10, and 25?Hz, respectively. The test results of the asphalt sand mastics were fitted with the Burger’s model to obtain the microparameters for DE models. To reduce the computation time of the DE modeling, the regular loading frequencies were amplified to virtual frequencies. Simultaneously, the Burger’s model parameters (microparameters in DE models) of asphalt sand mastic at regular frequencies were modified to those at virtual frequencies on the basis of the frequency-temperature superposition principle. Because the virtual frequencies were much larger than the regular frequencies, the computation time was significantly reduced by conducting the DE modeling with the virtual frequencies and the corresponding modified Burger’s model parameters. The modeling work, which typically takes several months or years with the traditional methods, only took a few hours or less in this study.     

Thermomechanical Framework for the Constitutive Modeling of Asphalt Concrete

This study is concerned with the constitutive modeling of asphalt concrete. Unlike most constitutive models for asphalt concrete that do not take into account the evolution of the microstructure of the material, this study incorporates the evolution of the microstructure by using a framework that recognizes that a body’s natural configurations can evolve as the microstructure changes. The general framework, on which this study is based, is cast within a full thermomechanical setting. In this paper, we develop models within the context of a mechanical framework that stems from the general framework for models based on the full thermodynamic framework and the resulting equations represent a nonlinear rate type viscoelastic model. The creep and stress relaxation experiments of Monismith and Secor are used for validating the efficacy of the model, and it is found that the predictions of the theory agree very well with the available experimental results. The advantages of using such a framework are many, especially when one wants to model the diverse mechanical and thermodynamic response characteristics of asphalt and asphalt concrete.     

Computational Constitutive Model for Predicting Nonlinear Viscoelastic Damage and Fracture Failure of Asphalt Concrete Mixtures

A computational constitutive model was developed to predict damage and fracture failure of asphalt concrete mixtures. Complex heterogeneity and inelastic mechanical behavior are addressed by the model by using finite-element methods and elastic–viscoelastic constitutive relations. Damage evolution due to progressive cracking is represented by randomly oriented interface fracture, which is governed by a newly developed nonlinear viscoelastic cohesive zone model. Computational simulations demonstrate that damage evolution and failure of asphalt concrete mixtures is dependent on the mechanical properties of the mixture. This approach is suitable for the relative evaluation of asphalt concrete mixtures by simply employing material properties and fracture properties of mixture components rather than by performing expensive laboratory tests recursively, which are typically required for continuum damage mechanics modeling.     

Linear Elastic and Viscoelastic Deformation Behavior of Ice

The dependence of overall elastic and viscoelastic behavior of polycrystals on properties of hexagonal single crystals of ice is predicted using established methods of micromechanical analysis. The self-consistent and bounding methods are used first to relate elastic behavior of S2 columnar and granular ice to single crystal properties. The results for S2 ice are then extended to viscoelastic behavior in order to relate overall creep behavior, as well as other time and rate-dependent behavior, to crystal parameters. Good agreement with experimental results is shown. Emphasis in this paper is on linear viscoelastic behavior of S2 ice due only to intragranular creep. The effect of this creep and grain boundary sliding on microcracking is considered briefly in the concluding remarks.     

Finite-Element Analysis of Hot Mix Asphalt Microstructure Using Effective Local Material Properties and Strain Gradient Elasticity

This paper presents the development of a methodology for microstructure analysis and modeling of hot mix asphalt (HMA). This methodology relies on using effective local material properties and strain gradient theory in the finite element analysis of HMA microstructure. The effective local properties are calculated using an analytical micromechanical model that captures the influence of percent of particles on the microscopic response of HMA. Strain gradient elasticity is used in order to account for the effect of particle size in the finite element analysis. The autocorrelation function and the moving window technique are used to determine the microstructure characteristic length scales that are used in strain gradient elasticity. A number of asphalt mixes with different aggregate types and size distributions are analyzed in this paper.