Research Papers

Multiscale Modeling of Damage Evolution in Viscoelastic Bituminous Mixtures Subjected to Cyclic Loading

[+] Author and Article Information
Yong-Rak Kim

Associate Professor
Department of Civil Engineering,
224 Engineering Building,
Kyung Hee University,
Gyeonggi-do, Korea446-701
e-mail: ykim3@khu.ac.kr

Flavio V. Souza

Research and Development Director
Multimech Research and Development,
LLC Omaha, NE 68022
e-mail: fsouza@multimechrd.com

Taehyo Park

Department of Civil and
Environmental Engineering,
Hanyang University,
Seoul, Korea133-791
e-mail: cepark@hanyang.ac.kr

1Corresponding author.

Contributed by the Materials Division of ASME for publication in the JOURNAL OF ENGINEERING MATERIALS AND TECHNOLOGY. Manuscript received May 30, 2012; final manuscript received February 14, 2013; published online March 25, 2013. Assoc. Editor: Xi Chen.

J. Eng. Mater. Technol 135(2), 021005 (Mar 25, 2013) (9 pages) Paper No: MATS-12-1116; doi: 10.1115/1.4023773 History: Received May 30, 2012; Revised February 14, 2013

This study presents a multiscale computational model and its application to predict damage dependent mechanical behavior of bituminous mixtures subjected to cyclic loading. Two length scales (global and local) are two-way coupled in the model framework by linking a homogenized global scale to a heterogeneous local scale representative volume element. Based on the unique two-way coupled multiscaling and the use of the finite element technique incorporated with the material viscoelasticity and cohesive zone fracture, the model approach can successfully account for the effect of mixture heterogeneity, material viscoelasticity, and damage accumulation due to cracks in the small scale on the overall performance of larger scale mixtures or structures. This step requires only the properties of individual constituents. To demonstrate the model and its features, bending beam fatigue testing of a bituminous mixture, which is composed of elastic aggregates and viscoelastic bitumen, is simulated by altering the mixture's constituent properties. The model clearly presents progressive damage characteristics with repetitive loading cycles and the analysis clearly demonstrates the sensitivity of the approach to constituent material properties. The multiscale model presented herein is expected to drastically reduce time-consuming and expensive fatigue tests, which, when performed in the traditional manner, require many replicates and do not define the cause of microstructural fatigue, damage, and failure.

Copyright © 2013 by ASME
Your Session has timed out. Please sign back in to continue.


Masad, E., Niranjanan, S., Bahia, H., and Kose, S., 2001, “Modeling and Experimental Measurements of Localized Strain Distribution in Asphalt Mixes,” J. Transp. Eng., 127(6), pp. 477–485. [CrossRef]
Soares, J. B., Freitas, F. A., and Allen, D. H., 2003, “Crack Modeling of Asphaltic Mixtures Considering Heterogeneity of the Material,” Transp. Res. Rec., 1832, pp. 113–120. [CrossRef]
Abbas, A., Masad, E., Papagiannakis, T., and Shenoy, A., 2005, “Modeling of Asphalt Mastic Stiffness Using Discrete Elements and Micromechanics Analysis,” Int. J. Pavement Eng., 6(2), pp. 137–146. [CrossRef]
Dai, Q., Sadd, M. H., and You, Z., 2006, “A Micromechanical Finite Element Model for Linear and Damage-Coupled Viscoelastic Behaviour of Asphalt Mixture,” Int. J. Numer. Analyt. Meth. Geomech., 30(11), pp. 1135–1158. [CrossRef]
Kim, Y., Allen, D. H., and Little, D. N., 2007, “Computational Constitutive Model for Predicting Nonlinear Viscoelastic Damage and Fracture Failure of Asphalt Concrete Mixtures,” Int. J. Geomech., 7(2), pp. 102–110. [CrossRef]
Aragão, F. T. S., Kim, Y., Lee, J., and Allen, D. H., 2011, “Micromechanical Model for Heterogeneous Asphalt Concrete Mixtures Subjected to Fracture Failure,” J. Mater. Civ. Eng., 23(1), pp. 30–38. [CrossRef]
Feng, Z., Zhang, P., Guddati, M. N., Kim, Y. R., and Little, D. N., 2008, “Virtual Testing Procedure for Cracking Performance Prediction of HMA,” Proceedings of the 6th RILEM International Conference on Cracking in Pavements, Chicago, June 16–18, pp. 575–585.
Soares, R. F., Kim, Y., and Allen, D. H., 2008, “Multiscale Computational Modeling for Predicting Evolution of Damage in Asphaltic Pavements,” Proceedings of the 6th RILEM International Conference on Cracking in Pavements, Chicago, June 16–18, pp. 599–608.
Lutif, J., Souza, F. V., Kim, Y., Soares, J. B., and Allen, D. H., 2010, “Multiscale Modeling to Predict the Mechanical Behavior of Asphalt Mixtures,” Transp. Res. Rec., 2181, pp. 28–35. [CrossRef]
Souza, F. V., 2009, “Multiscale Modeling of Impact on Heterogeneous Viscoelastic Solids With Evolving Microcracks,” Ph.D. dissertation, University of Nebraska–Lincoln, Lincoln, NE.
Souza, F. V. and Allen, D. H., 2010, “Multiscale Modeling of Impact on Heterogeneous Viscoelastic Solids Containing Evolving Microcracks,” Int. J. Numer. Methods Eng., 82(4), pp. 464–504. [CrossRef]
Allen, D. H., 2001, “Homogenization Principles and Their Application to Continuum Damage Mechanics,” Compos. Sci. Technol., 61, pp. 2223–2230. [CrossRef]
Dugdale, D., 1960, “Yielding of Steel Sheets Containing Slits,” J. Mech. Phys. Solids, 8, pp. 100–104. [CrossRef]
Barenblatt, G. I., 1962, “Mathematical Theory of Equilibrium Cracks in Brittle Fracture,” Adv. Appl. Mech., 7, pp. 55–129. [CrossRef]
Needleman, A., 1987, “A Continuum Model for Void Nucleation by Inclusion Debonding,” Trans ASME, J. Appl. Mech., 54(3), pp. 525–531. [CrossRef]
Tvergaard, V., 1990, “Effect of Fiber Debonding in a Whisker-Reinforced Metal,” Mater. Sci. Eng., 125(2), pp. 203–213. [CrossRef]
Geubelle, P. H. and Baylor, J. S., 1998, “Impact-Induced Delamination of Composites: A 2D Simulation,” Composites, Part B, 29(5), pp. 589–602. [CrossRef]
Ortiz, M. and Pandolfi, A., 1999, “Finite-Deformation Irreversible Cohesive Elements for Three-Dimensional Crack Propagation Analysis,” Int. J. Numer. Methods Eng., 44(9), pp. 1267–1282. [CrossRef]
Yoon, C. and Allen, D. H., 1999, “Damage Dependent Constitutive Behavior and Energy Release Rate for a Cohesive Zone in a Thermoviscoelastic Solid,” Int. J. Fract., 96, pp. 55–74. [CrossRef]
Park, K., Paulino, G. H., and Roesler, J. R., 2009, “A Unified Potential-Based Cohesive Model of Mixed-Mode Fracture,” J. Mech. Phys. Solids, 57(6), pp. 891–908. [CrossRef]
Kim, Y., Allen, D. H., and Little, D. N., 2005, “Damage-Induced Modeling of Asphalt Mixtures Through Computational Micromechanics and Cohesive Zone Fracture,” J. Mater. Civ. Eng., 17(5), pp. 477–484. [CrossRef]
Nemat-Nasser, S. and Hori, M., 1993, Micromechanics: Overall Properties of Heterogeneous Materials, North Holland, New York.
Kouznetsova, V. G., 2002, “Computational Homogenization for the Multi-Scale Analysis of Multi-Phase Materials,” Ph.D. dissertation, Technische Universiteit Eindhoven, Eindhoven, The Netherlands.
Kim, Y., Souza, F. V., and Teixeira, J. L., 2013, “A Two-Way Coupled Multiscale Model for Predicting Damage-Associated Performance of Asphaltic Roadways,” Comput. Mech., 51, pp. 187–201. [CrossRef]
Kim, Y., Little, D. N., and Lytton, R. L., 2003, “Fatigue and Healing Characterization of Asphalt Mixtures,” J. Mater. Civ. Eng., 15(1), pp. 75–83. [CrossRef]
Cheng, D., Little, D. N., Lytton, R. L., and Holste, J. C., 2003, “Moisture Damage Evaluation of Asphalt Mixtures by Considering Both Moisture Diffusion and Repeated-Load Conditions,” Transp. Res. Rec., 1832, pp. 42–58. [CrossRef]


Grahic Jump Location
Fig. 1

Multiscale modeling concept with two length scales (local and global)

Grahic Jump Location
Fig. 2

A simple flow chart of the multiscale computational algorithm

Grahic Jump Location
Fig. 3

Stress relaxation modulus curves of the two sand mastics

Grahic Jump Location
Fig. 5

Vertical displacements observed at the bottom center of the beam

Grahic Jump Location
Fig. 6

The σxx observed from a global scale element located at the bottom center of the beam

Grahic Jump Location
Fig. 4

Overall modeling process and local-global finite element meshes

Grahic Jump Location
Fig. 7

Evolution of the Cxxxx and Cyyyy components of the tangent constitutive tensor

Grahic Jump Location
Fig. 8

Snapshots of contour plots resulting from the simulation of the B2 mixture



Some tools below are only available to our subscribers or users with an online account.

Related Content

Customize your page view by dragging and repositioning the boxes below.

Related Journal Articles
Related eBook Content
Topic Collections

Sorry! You do not have access to this content. For assistance or to subscribe, please contact us:

  • TELEPHONE: 1-800-843-2763 (Toll-free in the USA)
  • EMAIL: asmedigitalcollection@asme.org
Sign In