Research Papers

Multifunctional Optimization of Viscoelastic Materials Subjected to Spherical Impact

[+] Author and Article Information
Martin Herrenbrück

Fachgebiet Computational Mechanics,
Technische Universität München,
München 80333, Germany
e-mail: herrenbruck@gmail.com

Fabian Duddeck

Fachgebiet Computational Mechanics,
Technische Universität München,
München 80333, Germany
e-mail: duddeck@tum.de

Roman Lackner

Arbeitsbereich für Materialtechnologie,
Leopold-Franzens-Universität Innsbruck,
Innsbruck 6020, Austria
e-mail: Roman.Lackner@uibk.ac.at

1Corresponding author.

Contributed by the Applied Mechanics Division of ASME for publication in the JOURNAL OF APPLIED MECHANICS. Manuscript received June 13, 2015; final manuscript received September 3, 2015; published online September 28, 2015. Assoc. Editor: Weinong Chen.

J. Appl. Mech 82(12), 121009 (Sep 28, 2015) (8 pages) Paper No: JAM-15-1314; doi: 10.1115/1.4031554 History: Received June 13, 2015; Revised September 03, 2015

The event of rigid solids impacting a viscoelastic body is encountered in many engineering disciplines. However, for this problem no analytical solution is available, making the proper design of, e.g., protective bodies a difficult task. As a remedy, generally valid solutions in form of nondimensional response curves are presented in this paper, serving as reference for validation purposes or as design charts for the performance-oriented development of impact absorbers. Hereby, the problem of a frictionless rigid sphere impinging a viscoelastic half-space is investigated by finite-element analyses. The general applicability of the results is assured by a transformation to dimensionless problem parameters. For this purpose, the analytical solution by Hertz (1881, “Über die Berührung fester elastischer Körper,” J. die Reine Angew. Math., 1882(92), pp. 156–171) for the purely elastic impact is taken into account. The chosen nondimensional format leads to a reduced number of system parameters, allowing for a compact representation by so-called master curves. From these, optimal material characteristics are found for three different design objectives.

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


Argatov, I. I. , 2013, “ Mathematical Modeling of Linear Viscoelastic Impact: Application to Drop Impact Testing of Articular Cartilage,” Tribol. Int., 63, pp. 213–225. [CrossRef]
Hertz, H. , 1881, “ Über die Berührung fester elastischer Körper,” J. die Reine Angew. Math., 1882(92), pp. 156–171.
Pao, Y.-H. , 1955, “ Extension of the Hertz Theory of Impact to the Viscoelastic Case,” J. Appl. Phys., 26(9), pp. 1083–1088. [CrossRef]
Radok, J. , 1957, “ Visco-Elastic Stress Analysis,” Q. Appl. Math., 15(2), pp. 198–202.
Lee, E. H. , and Radok, J. R. M. , 1960, “ The Contact Problem for Viscoelastic Bodies,” ASME J. Appl. Mech., 27(3), pp. 438–444. [CrossRef]
Hunter, S. C. , 1960, “ The Hertz Problem for a Rigid Spherical Indenter and a Viscoelastic Half-Space,” J. Mech. Phys. Solids, 8(4), pp. 219–234. [CrossRef]
Graham, G. A. C. , 1965, “ The Contact Problem in the Linear Theory of Viscoelasticity,” Int. J. Eng. Sci., 3(1), pp. 27–46. [CrossRef]
Yang, W. H. , 1966, “ The Contact Problem for Viscoelastic Bodies,” ASME J. Appl. Mech., 33(2), pp. 395–401. [CrossRef]
Ting, T. C. T. , 1966, “ The Contact Stresses Between a Rigid Indenter and a Viscoelastic Half-Space,” ASME J. Appl. Mech., 33(4), pp. 845–854. [CrossRef]
Graham, G. A. C. , 1967, “ The Contact Problem in the Linear Theory of Viscoelasticity When the Time Dependent Contact Area Has Any Number of Maxima and Minima,” Int. J. Eng. Sci., 5(6), pp. 495–514. [CrossRef]
Ting, T. C. T. , 1968, “ Contact Problems in Linear Theory of Viscoelasticity,” ASME J. Appl. Mech., 35(2), pp. 248–254. [CrossRef]
Hunter, S. C. , 1967, “ The Solution of Boundary Value Problems in Linear Viscoelasticity,” Mechanics and Chemistry of Solid Propellants: Proceedings of the Fourth Symposium on Naval Structural Mechanics, A. C. Eringen , H. Liebowitz , S. L. Koh, and J. M. Crowley , eds., Pergamon Press, New York, pp. 257–295.
Ogibalov, P. , Koltunov, M. , and Soldatov, M. , 1969, “ Formulation and Solution of Boundary Value Problems of Polymer Mechanics,” Mech. Compos. Mater., 5(1), pp. 43–51.
Kalker, J. , 1977, “ A Survey of the Mechanics of Contact Between Solid Bodies,” ZAMM, 57(5), pp. T3–T17. [CrossRef]
Calvit, H. H. , 1967, “ Numerical Solution of the Problem of Impact of a Rigid Sphere Onto a Linear Viscoelastic Half-Space and Comparison With Experiment,” Int. J. Solids Struct., 3(6), pp. 951–960. [CrossRef]
Aboudi, J. , 1979, “ The Dynamic Indentation and Impact of a Viscoelastic Half-Space by an Axisymmetric Rigid Body,” Comput. Methods Appl. Mech. Eng., 20(2), pp. 135–150. [CrossRef]
Sabin, G. C. W. , 1987, “ The Impact of a Rigid Axisymmetric Indentor on a Viscoelastic Half-Space,” Int. J. Eng. Sci., 25(2), pp. 235–251. [CrossRef]
Aksel, N. , 1986, “ On the Impact of a Rigid Sphere on a Viscoelastic Half-Space,” Arch. Appl. Mech., 56(1), pp. 38–54.
Chen, C. P. , and Lakes, R. S. , 1990, “ Design of Viscoelastic Impact Absorbers: Optimal Material Properties,” Int. J. Solids Struct., 26(12), pp. 1313–1328. [CrossRef]
Argatov, I. , 2012, “ An Analytical Solution of the Rebound Indentation Problem for an Isotropic Linear Viscoelastic Layer Loaded With a Spherical Punch,” Acta Mech., 223(7), pp. 1441–1453. [CrossRef]
Findley, W. N. , Lai, J. S. , and Onaran, K. , 1976, Creep and Relaxation of Nonlinear Viscoelastic Materials—With an Introduction to Linear Viscoelasticity, North Holland Publishing Company, Amsterdam.
Christensen, R. , 1982, Theory of Viscoelasticity: An Introduction, Academic Press, New York.
Lakes, R. S. , 1998, Viscoelastic Solids, CRC Press, Boca Raton, FL.
Deresiewicz, H. , 1968, “ A Note on Hertz's Theory of Impact,” Acta Mech., 6(1), pp. 110–112. [CrossRef]
Johnson, K. L. , 1985, Contact Mechanics, Cambridge University Press, Cambridge, UK.
ansys 12.1., 2009, Theory Reference for the Mechanical APDL and Mechanical Applications, ANSYS, Canonsburg, PA.
Cheng, Y.-T. , and Cheng, C.-M. , 2004, “ Scaling, Dimensional Analysis, and Indentation Measurements,” Mater. Sci. Eng.: R: Rep., 44(4–5), pp. 91–149. [CrossRef]
Herrenbrück, M. , 2013, “ Finite-Element Based Determination of Response Spectra of Viscoelastic Materials Subjected to Low-Velocity Impact Loading,” Ph.D. thesis, Technische Universität München, Munich, Germany.
Ashby, M. F. , 2005, Materials Selection in Mechanical Design, Butterworth-Heinemann, Boston, MA.


Grahic Jump Location
Fig. 1

Rheological models for linear viscoelasticity: (a) the Maxwell element and (b) the standard linear solid

Grahic Jump Location
Fig. 2

Finite-element model where use of rotational symmetry was made

Grahic Jump Location
Fig. 3

Example results (normalized to unit range) for the penetration, velocity, kinetic energy, and acceleration over time

Grahic Jump Location
Fig. 4

Maxwell model—each symbol represents the result from one numerical simulation: (a), (c), and (e) dimensional results and (b), (d), and (f) master curves obtained by the use of dimensionless abscissa and ordinate. (a) Dependence of uve on τ, (b) Πu over Πτθ, (c) dependence of vfinal on τ, (d) Πrest over Πτθ, (e) dependence of ave on τ, and (f) Πa over Πτθ.

Grahic Jump Location
Fig. 5

Standard linear solid, where five different values for α were used: (a), (c), and (e) dimensional results and (b), (d), and (f) master curves obtained by the same scaling as in Fig. 4. The Hertzian solution refers to the initial stiffness E0. The horizontal lines given in (b) and (f) for small Πτθ indicate the limit values (1−α)−0.4 and (1−α)+0.4, respectively. (a) Dependence of uve on τ, (b) Πu over Πτθ, (c) dependence of vfinal on τ, (d) Πrest over Πτθ, (e) dependence of ave on τ, and (f) Πa over Πτθ.

Grahic Jump Location
Fig. 6

Standard linear solid (α=0.9): dimensional results from a numerical parameter study for varying characteristic times τ in the range from 0.001 s to 4.0 s. For each of the seven simulations, the acceleration is plotted versus the absolute value of the penetration, where the physical time is a curve parameter. From this, the maximum values for acceleration and penetration are marked by diamonds. For comparison, the dotted line shows the relationship between uel and ael according to the Hertz solution.

Grahic Jump Location
Fig. 7

Standard linear solid: plot of Πa over Πu (i.e., nondimensional representation), where the solid line indicates the Hertzian solution using E=(1−α)E0 for a varying parameter α

Grahic Jump Location
Fig. 8

Maxwell model: multi-objective function D(Πτθ) for different weights wu=1−wa, where the symbols represent result points from numerical simulations corresponding to D=du (i.e., wu=1) and D=da (i.e., wu=0), respectively

Grahic Jump Location
Fig. 9

Optima according to the L2 and the L∞ norm: (a) plot of the normalized distances du and da over the design variable Πτθ. The sum of the squared distances and its square root (L2 norm) are minimum at Πτθ=0.301. The intersection of du and da is located at Πτθ=0.261. (b) L2 and L∞ optimum in the objective function space. All designs are Pareto optimal.

Grahic Jump Location
Fig. 10

Maxwell model: Identification of Πτθ=τ/θ for three different constraints on the maximum penetration, i.e., three different values for du, and the respective result for da

Grahic Jump Location
Fig. 11

Maxwell model: relation between du, uve, da, and ave for the problem parameters listed in Table 3

Grahic Jump Location
Fig. 12

Maxwell model: Sensitivity of uve and ave with respect to τ for the cases ① (uve100%=5 mm, full lines), ② (uve100%=7 mm, dashed lines), and ③ (uve100%=12 mm, dotted lines), see Table 5



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