Research Papers

Numerical Modeling of the Nonlinear Elastic Response of Filled Elastomers via Composite-Sphere Assemblages

[+] Author and Article Information
Taha Goudarzi

e-mail: goudarz2@illinois.edu

Oscar Lopez-Pamies

e-mail: pamies@illinois.edu
Department of Civil
and Environmental Engineering,
University of Illinois Urbana–Champaign,
IL 61801-2352

An extension of this approach to nonlinear, though convex, problems appears to have been first carried out by Barrett and Talbot [7] in the context of two-phase dielectrics.

The analogous kinematically admissible approximation is discussed in the Appendix.

In other words, the composite sphere is a composite structure and not a composite material.

The initial shear moduli of standard reinforcing fillers (e.g., silica) are typically four orders of magnitude larger than those of standard elastomers (e.g., silicone).

Manuscript received July 21, 2012; final manuscript received January 18, 2013; accepted manuscript posted January 28, 2013; published online July 19, 2013. Assoc. Editor: Martin Ostoja-Starzewski.

J. Appl. Mech 80(5), 050906 (Jul 19, 2013) (10 pages) Paper No: JAM-12-1339; doi: 10.1115/1.4023497 History: Received July 21, 2012; Revised January 18, 2013; Accepted January 28, 2013

This paper proposes an effective numerical method to generate approximate solutions for the overall nonlinear elastic response of isotropic filled elastomers subjected to arbitrarily large deformations. The basic idea is first to idealize the random microstructure of isotropic filled elastomers as an assemblage of composite spheres and then to generate statically admissible numerical solutions, via finite elements, for these material systems directly in terms of the response of a single composite sphere subjected to affine stress boundary conditions. The key theoretical strengths of the method are discussed, and its accuracy and numerical efficiency assessed by comparisons with corresponding 3D full-field simulations. The paper concludes with a discussion of straightforward extensions of the proposed method to account for general classes of anisotropic microstructures and filler-elastomer interphasial phenomena, features of key importance in emerging advanced applications.

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


Leblanc, J. L., 2010, Filled Polymers: Science and Industrial Applications, CRC Press, Boca Raton, FL.
Smallwood, H. M., 1944, “Limiting Law of the Reinforcement of Rubber,” J. Appl. Phys., 15, pp. 758–766. [CrossRef]
Guth, E., 1945, “Theory of Filler Reinforcement,” J. Appl. Phys., 16, pp. 20–25. [CrossRef]
Hill, R., 1972, “On Constitutive Macrovariables for Heterogeneous Solids at Finite Strain,” Proc. R. Soc. Lond. A, 326, pp. 131–147. [CrossRef]
Lopez-Pamies, O., Goudarzi, T., and Nakamura, T., 2013, “The Nonlinear Elastic Response of Suspensions of Rigid Inclusions in Rubber: I—An Exact Result for Dilute Suspensions,” J. Mech. Phys. Solids, 61, pp. 1–18. [CrossRef]
Hashin, Z., 1962, “The Elastic Moduli of Heterogeneous Materials,” ASME J. Appl. Mech., 29, pp. 143–150. [CrossRef]
Barrett, K. E., and Talbot, D. R. S., 1996, “Bounds for the Effective Properties of a Nonlinear Two-Phase Composite Dielectric,” Proc. 8th Int. Symp. Continuum Models and Discrete Systems, Varna, Bulgaria, June 11–16, 1995, World Scientific, Singapore.
Carpi, F., De Rossi, D., Kornbluh, R., Pelrine, R., and Sommer-Larsen, P., 2008, Dielectric Elastomers as Electromechanical Transducers, Elsevier, New York.
Danas, K., Kankanala, S. V., and Triantafyllidis, N., 2012, “Experiments and Modeling of Iron-Particle-Filled Magnetorheological Elastomers,” J. Mech. Phys. Solids, 60, pp. 120–138. [CrossRef]
Ogden, R. W., 1978, “Extremum Principles in Non-Linear Elasticity and Their Application to Composites—I Theory,” Int. J. Solids. Struct., 14, pp. 265–282. [CrossRef]
Lee, S. J., and Shield, R. T., 1980, “Variational Principles in Finite Elastostatics, J. Appl, Math. Phys., 31, pp. 437–453. [CrossRef]
Khisaeva, Z. F., and Ostoja-Starzewski, M., 2006, “Mesoscale Bounds in Finite Elasticity and Thermoelasticity of Random Composites,” Proc. R. Soc. Lond. A, 462, pp. 1167–1180. [CrossRef]
Ogden, R. W., 1997, Non-Linear Elastic Deformations, Dover Publications, Mineola, NY.
Ramier, J., 2004, “Comportement Mécanique d’élastomérs Chargés, Influence de l'adhésion Charge–Polymére, Influence de la Morphologie,” Ph. D. dissertation, Institut National des Science Appliquées de Lyon, France.
Milton, G. W., 2002, The Theory of Composites (Cambridge Monographs on Applied and Computational Mathematics, Vol. 6), Cambridge University Press, Cambridge, UK.
Willis, J. R., 1989, “The Structure of Overall Constitutive Relations for a Class of Nonlinear Composites, IMA J. Appl. Math., 43, pp. 231–242. [CrossRef]
Dassault Systemes Simulia Corp., 2011, ABAQUS Version 6.11 Documentation, Providence, RI.
Lopez-Pamies, O., 2010,“A New I1-Based Hyperelastic Model for Rubber Elastic Materials,” Comptes Rendus Mecanique, 338, pp. 3–11. [CrossRef]
Lopez-Pamies, O., Goudarzi, T., and Danas, K., 2013, “The Nonlinear Elastic Response of Suspensions of Rigid Inclusions in Rubber: II—A Simple Explicit Approximation for Finite-Concentration Suspensions,” J. Mech. Phys. Solids, 61, pp. 19–37. [CrossRef]
Bornert, M., Stolz, C., and Zaoui, A., 1996, “Morphologically Representative Pattern-Based Bounding in Elasticity,” J. Mech. Phys. Solids, 44, pp. 307–331. [CrossRef]
Bertoldi, K., and Lopez-Pamies, O., 2012, “Some Remarks on the Effect of Interphases on the Mechanical Response and Stability of Fiber-Reinforced Elastomers,” ASME J. Appl. Mech., 79, p. 031023. [CrossRef]
Qu, M., Deng, F., Kalkhoran, S. M., Gouldstone, A., Robisson, A., and Van Vliet, K. J., 2011, “Nanoscale Visualization and Multiscale Mechanical Implications of Bound Rubber Interphases in Rubber-Carbon Black Nanocomposites,” Soft Matter, 7, pp. 1066–1077. [CrossRef]


Grahic Jump Location
Fig. 1

(a) Electron micrograph of a styrene-butadiene rubber filled with an isotropic distribution of silica particles [14] and (b) its idealization as a composite-sphere assemblage (CSA) in the undeformed configuration. All the composite spheres in the assemblage are homothetic in that they have the same ratio of inner-to-outer radius Ri/Ro = c1/3.

Grahic Jump Location
Fig. 4

Contour plots of the maximum principal logarithmic strain for a composite sphere with c = 0.15, Neo-Hookean matrix, and 104–times stiffer Neo-Hookean particle subjected to affine uniaxial stress (24) with S¯ = diag(s¯1 > 0,0,0); the undeformed configuration is also depicted for comparison purposes. The overall stretch in the direction of applied stress is λ¯1 = 3.5.

Grahic Jump Location
Fig. 2

Schematic illustrating that the overall response of a CSA subjected to affine stress boundary conditions can be variationally approximated by the overall response of a corresponding single composite sphere subjected to the same affine stress boundary conditions. Specifically, the approximation is such that the total elastic energy W¯ of the CSA is bounded from below by the total elastic energy W¯S of the single composite sphere.

Grahic Jump Location
Fig. 3

Three representative meshes in the undeformed configuration for a composite sphere with particle concentration c = 0.15: (a) coarse mesh with 28,400 elements, (b) fine mesh with 102,600 elements, and (c) very fine mesh with 260,800 elements

Grahic Jump Location
Fig. 5

The normalized initial effective shear modulus μ¯/μ of isotropic incompressible elastomers filled with random isotropic distributions of rigid particles. Plots are shown for: (i) the CSA approximation μ¯S, (ii) full-field FE simulations, (iii) the Hashin–Shtrikman lower bound μ¯HS/μ = (2+3c)/(2-2c), and (iv) the Reuss lower bound μ¯R/μ = 1/(1-c), as functions of the concentration of particles c.

Grahic Jump Location
Fig. 6

Macroscopic response of Neo-Hookean rubber filled with an isotropic distribution of rigid particles of concentration c = 0.15. Part (a) displays the entire effective stored-energy function in terms of the macroscopic principal stretches λ¯1 and λ¯2, whereas part (b) shows the energy along axisymmetric loading conditions with λ¯1 = λ¯, λ¯2 = λ¯-1/2. Results are shown for the CSA approximation W¯S in both parts, and for corresponding full-field FE simulations for isotropic distributions of spherical particles in part (b).

Grahic Jump Location
Fig. 9

(a) Electron micrograph of a magnetorheological elastomer with iron particles distributed anisotropically in chainlike structures [9] and (b) its idealization as an ellipsoidal assemblage of possibly nonspherical particles (CEA). All the composite ellipsoids in the assemblage are homothetic in that they are scaled-up or scaled-down versions of each other. Part (b) also illustrates schematically the straightforward incorporation of bound rubber into the CEA idealization.

Grahic Jump Location
Fig. 7

Macroscopic response of filled Neo-Hookean rubber with various values of concentration of particles c under: (a) uniaxial compressive, (b) uniaxial tensile, (c) pure shear, and (d) simple shear loading conditions. Plots are shown for the CSA approximation and corresponding full-field FE simulations for isotropic distributions of spherical particles.

Grahic Jump Location
Fig. 10

Comparisons between the kinematically admissible approximation (A1) and (A2), denoted as CSA-F¯, and the statically admissible approximation Eqs. (17)–(20), denoted as CSA-S¯, for the overall response of filled Neo-Hookean rubber. Part (a) shows results for the normalized initial shear modulus μ¯/μ as a function of particle concentration c, while part (b) shows stress-deformation results for c = 0.15 under uniaxial tension.

Grahic Jump Location
Fig. 8

Macroscopic response of filled silicone rubber with various values of concentration of particles c under: (a) uniaxial compressive, (b) uniaxial tensile, (c) pure shear, and (d) simple shear loading conditions. Plots are shown for the CSA approximation and corresponding full-field simulations for isotropic distributions of spherical particles.



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