Research Papers: Design Automation

Setting Material Function Design Targets for Linear Viscoelastic Materials and Structures

[+] Author and Article Information
R. E. Corman

Department of Mechanical
Science and Engineering,
University of Illinois at Urbana-Champaign,
Urbana, IL 61801
e-mail: corman1@illinois.edu

Lakshmi Rao

Department of Industrial and Enterprise
Systems Engineering,
University of Illinois at Urbana-Champaign,
Urbana, IL 61801

N. Ashwin Bharadwaj

Department of Mechanical
Science and Engineering,
University of Illinois at Urbana-Champaign,
Urbana, IL 61801

James T. Allison

Assistant Professor
Department of Industrial and Enterprise
Systems Engineering,
University of Illinois at Urbana-Champaign,
Urbana, IL 61801

Randy H. Ewoldt

Assistant Professor
Department of Mechanical
Science and Engineering,
University of Illinois at Urbana-Champaign,
Urbana, IL 61801

Contributed by the Design Automation Committee of ASME for publication in the JOURNAL OF MECHANICAL DESIGN. Manuscript received March 9, 2015; final manuscript received January 25, 2016; published online March 21, 2016. Assoc. Editor: Nam H. Kim.

J. Mech. Des 138(5), 051402 (Mar 21, 2016) (12 pages) Paper No: MD-15-1206; doi: 10.1115/1.4032698 History: Received March 09, 2015; Revised January 25, 2016

Rheologically complex materials are described by function-valued properties with dependence on a timescale (linear viscoelasticity), input amplitude (nonlinear material behavior), or more generally both (nonlinear viscoelasticity). This complexity presents a difficulty when trying to utilize these material systems in engineering designs. Here, we focus on linear viscoelasticity and a methodology to identify the desired viscoelastic behavior. This is an early-stage design step to optimize target (function-valued) properties before choosing or synthesizing a real material. In linear viscoelasticity, it is not obvious which properties can be treated as independent design variables. Thus, it is nontrivial to select the most design-appropriate constitutive model, to be as general as possible, but not violate fundamental restrictions. We use the Kramers–Kronig constraint to show that frequency-dependent moduli (e.g., shear moduli G(ω) and G(ω)) cannot be treated as two independent design variables. Rather, a single function such as the relaxation modulus (e.g., K(t) for force-relaxation or G(t) for stress relaxation) is an appropriate function-valued design variable. A simple case study is used to demonstrate the framework in which we identify target properties for a vibration isolation system. Viscoelasticity improves performance. Different parameterizations of the kernel function are optimized and compared for performance. While parameterization may limit the generality of the kernel function, we do include a nonobvious representation (power law) that is found in real viscoelastic material systems and in the spring-dashpot paradigm would require an infinite number of components. Our methodology provides a means to answer the question, “What viscoelastic properties are desirable?” This ability to identify targeted behavior will be useful for subsequent stages of the design process including the selection or synthesis of real materials.

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


Autumn, K. , Liang, Y. A. , Hsieh, S. T. , Zesch, W. , Chan, W. P. , Kenny, T. W. , Fearing, R. , and Full, R. J. , 2000, “ Adhesive Force of a Single Gecko Foot-Hair,” Nature, 405(6787), pp. 681–685. [CrossRef] [PubMed]
Shadwick, R. E. , 1999, “ Mechanical Design in Arteries,” J. Exp. Biol., 202, pp. 3305–3313. [PubMed]
Denny, M. , 1980, “ The Role of Gastropod Pedal Mucus in Locomotion,” Nature, 285(5761), pp. 160–161. [CrossRef]
Caldwell , 1972, “ Vibration Damping in World Trade Center Using Viscoelastic Material,” 3M.
Mahmoodi, P. , Robertson, L. E. , Yontar, M. , Moy, C. , and Field, I. , 1987, “ Performance of Viscoelastic Dampers in World Trade Center Towers,” 3M.
de Gennes, P. G. , 1996, “ Soft Adhesives,” Langmuir, 12(19), pp. 4497–4500. [CrossRef]
Ewoldt, R. H. , 2014, “ Extremely Soft: Design With Rheologically Complex Fluids,” Soft Rob., 1(1), pp. 12–20. [CrossRef]
Yang, H. , and Chang, E. , 1997, “ The Role of Viscoelastic Properties in the Design of Pressure-Sensitive Adhesives,” Trends Polym. Sci., 5(11), pp. 380–384.
Nachbar, W. , and Schipmolder, J. B. , 1969, “ Optimization of a Viscoelastic Structure: The Seat-Belt Problem,” ASME J. Appl. Mech., 36(3), p. 565. [CrossRef]
Fischer, C. , Braun, S. A. , Bourban, P.-E. , Michaud, V. , Plummer, C. J. G. , and Månson, J.-A. E. , 2006, “ Dynamic Properties of Sandwich Structures With Integrated Shear-Thickening Fluids,” Smart Mater. Struct., 15(5), pp. 1467–1475. [CrossRef]
Lee, Y. S. , Wetzel, E. D. , and Wagner, N. J. , 2003, “ The Ballistic Impact Characteristics of Kevlar Woven Fabrics Impregnated With a Colloidal Shear Thickening Fluid,” J. Mater. Sci., 38(13), pp. 2825–2833. [CrossRef]
Nachbar, W. , and Schipmolder, J. B. , 1971, “ Optimization of a Standard Linear Viscoelastic Material for Use in a Seat Belt,” ASME J. Appl. Mech., 38(2), p. 532. [CrossRef]
Ferry, J. D. , 1980, Viscoelastic Properties of Polymers, 3rd ed., Wiley, New York.
Xu, H. , Li, Y. , Brinson, C. , and Chen, W. , 2014, “ A Descriptor-Based Design Methodology for Developing Heterogeneous Microstructural Materials System,” ASME J. Mech. Des., 136(5), p. 051007. [CrossRef]
Asbeck, A. , De Rossi, S. , Galiana, I. , Ding, Y. , and Walsh, C. , 2014, “ Stronger, Smarter, Softer: Next-Generation Wearable Robots,” Rob. Autom. Mag., 21(4), pp. 22–33. [CrossRef]
Shepherd, R. F. , Ilievski, F. , Choi, W. , Morin, S. A. , Stokes, A. A. , Mazzeo, A. D. , Chen, X. , Wang, M. , and Whitesides, G. M. , 2011, “ Multigait Soft Robot,” Proc. Natl. Acad. Sci. U. S. A., 108(51), pp. 20400–20403. [CrossRef] [PubMed]
Mahdavi, A. , Ferreira, L. , Sundback, C. , Nichol, J. W. , Chan, E. P. , Carter, D. J. D. , Bettinger, C. J. , Patanavanich, S. , Chignozha, L. , Ben-Joseph, E. , Galakatos, A. , Pryor, H. , Pomerantseva, I. , Masiakos, P. T. , Faquin, W. , Zumbuehl, A. , Hong, S. , Borenstein, J. , Vacanti, J. , Langer, R. , and Karp, J. M. , 2008, “ A Biodegradable and Biocompatible Gecko-Inspired Tissue Adhesive,” Proc. Natl. Acad. Sci. U. S. A., 105(7), pp. 2307–2312. [CrossRef] [PubMed]
Laulicht, B. , Langer, R. , and Karp, J. M. , 2012, “ Quick-Release Medical Tape,” Proc. Natl. Acad. Sci. U. S. A., 109(46), pp. 18803–18808. [CrossRef] [PubMed]
Kim, D.-H. , and Rogers, J. A. , 2008, “ Stretchable Electronics: Materials Strategies and Devices,” Adv. Mater., 20(24), pp. 4887–4892. [CrossRef]
Kim, D.-H. , Xiao, J. , Song, J. , Huang, Y. , and Rogers, J. A. , 2010, “ Stretchable, Curvilinear Electronics Based on Inorganic Materials,” Adv. Mater. (Deerfield Beach, Fla.), 22(19), pp. 2108–2124. [CrossRef]
Ashby, M. F. , 1999, Material Selection in Mechanical Design, Butterworth-Heinemann, Boston, MA.
Dieter, G. , 2009, Engineering Design: A Materials and Processing Approach, 4th ed., McGraw-Hill Higher Education, New York.
Olson, G. B. , 2000, “ Designing a New Material World,” Science, 288(5468), pp. 993–998. [CrossRef]
Cussler, E. L. , and Moggridge, G. D. , 2012, Chemical Product Design, 2nd ed., Cambridge University Press, Cambridge, UK.
Nelson, A. , and Ewoldt, R. H. , 2015, “ Design of Yield-Stress Fluids: A Rheology to Structure Inverse Problem,” (in press).
Nelson, A. , 2015, “ Expanding Yield-Stress Fluid Paradigms,” M.S. thesis, University of Illinois, Champaign, IL.
Grindy, S. C. , Learsch, R. , Mozhdehi, D. , Cheng, J. , Barrett, D. G. , Guan, Z. , Messersmith, P. B. , and Holten-Andersen, N. , 2015, “ Control of Hierarchical Polymer Mechanics With Bioinspired Metal-Coordination Dynamics,” Nat. Mater., 14(12), pp. 1210–1216. [CrossRef] [PubMed]
Chen, P. , Li, Q. , Grindy, S. , and Holten-Andersen, N. , 2015, “ White-Light-Emitting Lanthanide Metallogels With Tunable Luminescence and Reversible Stimuli-Responsive Properties,” J. Am. Chem. Soc., 137(36), pp. 11590–11593. [CrossRef] [PubMed]
Kim, H. M. , Michelena, N. F. , Papalambros, P. Y. , and Jiang, T. , 2003, “ Target Cascading in Optimal System Design,” ASME J. Mech. Des., 125(3), p. 474. [CrossRef]
Kim, H. M. , Rideout, D. G. , Papalambros, P. Y. , and Stein, J. L. , 2003, “ Analytical Target Cascading in Automotive Vehicle Design,” ASME J. Mech. Des., 125(3), p. 481. [CrossRef]
Choi, H. , McDowell, D. L. , Allen, J. K. , Rosen, D. , and Mistree, F. , 2008, “ An Inductive Design Exploration Method for Robust Multiscale Materials Design,” ASME J. Mech. Des., 130(3), p. 031402. [CrossRef]
Linsey, J. S. , Tseng, I. , Fu, K. , Cagan, J. , Wood, K. L. , and Schunn, C. , 2010, “ A Study of Design Fixation, Its Mitigation and Perception in Engineering Design Faculty,” ASME J. Mech. Des., 132(4), p. 041003. [CrossRef]
Hilton, H. H. , and Yi, S. , 1992, “ Analytical Formulation of Optimum Material Properties for Viscoelastic Damping,” Smart Mater. Struct., 1(2), p. 113. [CrossRef]
Hilton, H. H. , 2005, “ Optimum Linear and Nonlinear Viscoelastic Designer Functionally Graded Materials-Characterizations and Analysis,” Composites, Part A, 36(10), pp. 1329–1334. [CrossRef]
Prasad, J. , and Diaz, A. R. , 2008, “ A Concept for a Material That Softens With Frequency,” ASME J. Mech. Des., 130(9), p. 091703. [CrossRef]
Bharadwaj, N. A. , Allison, J. T. , and Ewoldt, R. H. , 2013, “ Early-Stage Design of Rheologically Complex Materials Via Material Function Design Targets,” ASME Paper No. DETC2013-13462.
Dealy, J. M. , 1995, “ Official Nomenclature for Material Functions Describing the Response of a Viscoelastic Fluid to Various Shearing and Extensional Deformations,” J. Rheol., 39(1), pp. 253–265. [CrossRef]
Bird, R. B. , Armstrong, R. , and Hassger, O. , 1987, Dynamics of Polymeric Liquids: Volume 1 Fluid Mechanics, 2nd ed., Wiley, New York.
Johnson, C. D. , 1995, “ Design of Passive Damping Systems,” ASME J. Mech. Des., 117(B), p. 171. [CrossRef]
Hilton, H. H. , Lee, D. H. , and El Fouly, A. R. A. , 2008, “ Generalized Viscoelastic Designer Functionally Graded Auxetic Materials Engineered/Tailored for Specific Task Performances,” Mech. Time-Depend. Mater., 12(2), pp. 151–178. [CrossRef]
Breuer, S. , and Onat, E. T. , 1962, “ On Uniqueness in Linear Elasticity,” Q. Appl. Math., 19(4), pp. 355–359.
Rabotnov, Y. N. , 1977, Elements of Hereditary Solid Mechanics (English Translation), Mir, Moscow, Russia.
Chambon, F. , and Winter, H. H. , 1987, “ Linear Viscoelasticity at the Gel Point of a Cross-Linking PDMS With Imbalanced Stoichiometry,” J. Rheol., 31(8), pp. 683–697. [CrossRef]
Winter, H. H. , and Chambon, F. , 1986, “ Analysis of Linear Viscoelasticity of a Cross-Linking Polymer at the Gel Point,” J. Rheol., 30(2), pp. 367–382. [CrossRef]
Koike, A. , Nemoto, N. , Inoue, T. , and Osaki, K. , 1995, “ Dynamic Light Scattering and Dynamic Viscoelasticity of Poly(Vinyl Alcohol) in Aqueous Borax Solutions. 1. Concentration Effect,” Macromolecules, 28(7), pp. 2339–2344. [CrossRef]
Corman, R. , 2015, “ Enabling Design With Rheological Complexity: Intuition and Optimization of Viscoelastic Materials,” M.S. thesis, University of Illinois at Urbana-Champaign, Champaign, IL.
Tsiklauri, D. , and Beresnev, I. , 2001, “ Non-Newtonian Effects in the Peristaltic Flow of a Maxwell Fluid,” Phys. Rev. E, 64(3), p. 036303. [CrossRef]
Cordobes, F. , Partal, P. , and Guerrero, A. , 2004, “ Rheology and Microstructure of Heat-Induced Egg Yolk Gels,” Rheol. Acta, 43(2), pp. 184–195. [CrossRef]
Trepat, X. , Deng, L. , An, S. S. , Navajas, D. , Tschumperlin, D. J. , Gerthoffer, W. T. , Butler, J. P. , and Fredberg, J. J. , 2007, “ Universal Physical Responses to Stretch in the Living Cell,” Nature, 447(7144), pp. 592–595. [CrossRef] [PubMed]
Winter, H. H. , and Mours, M. , 1997, “ Rheology of Polymers Near Liquid-Solid Transitions,” Neutron Spin Echo Spectroscopy, Viscoelasticity, Rheology (Advances in Polymer Science), Vol. 134, Springer, Berlin, pp. 165–234.
Kollmannsberger, P. , and Fabry, B. , 2011, “ Linear and Nonlinear Rheology of Living Cells,” Annu. Rev. Mater. Res., 41(1), pp. 75–97. [CrossRef]
Athans, M. , and Falb, P. L. , 2006, Optimal Control: An Introduction to the Theory and Its Applications, Dover Publications/McGraw-Hill, New York.
Betts, J. T. , 2010, Practical Methods for Optimal Control and Estimation Using Nonlinear Programming, SIAM, Philadelphia, PA.
Allison, J. T. , Guo, T. , and Han, Z. , 2014, “ Co-Design of an Active Suspension Using Simultaneous Dynamic Optimization,” ASME J. Mech. Des., 136(8), p. 081003. [CrossRef]
Karnopp, D. , 1995, “ Active and Semi-Active Vibration Isolation,” ASME J. Mech. Des., 117(B), p. 177. [CrossRef]
Inman, D. J. , 2013, Engineering Vibration, 4th ed., Prentice-Hall, New York.
Corman, R. E. , and Rao, L. G. , 2015, Linear Viscoelastic Design Optimization-File Exchange-MATLAB Central.
Donelan, M. A. , Hamilton, J. , and Hui, W. H. , 1985, “ Directional Spectra of Wind-Generated Waves,” Philos. Trans. R. Soc., A, 315(1534), pp. 509–562. [CrossRef]
Latheef, M. , and Swan, C. , 2013, “ A Laboratory Study of Wave Crest Statistics and the Role of Directional Spreading,” Proc. R. Soc. A, 469(2152), p. 20120696. [CrossRef]
Larson, R. G. , 1999, The Structure and Rheology of Complex Fluids, Oxford University Press, Oxford, UK.
Rubinstein, M. , and Colby, R. H. , 2003, Polymer Physics, Oxford University Press, Oxford, UK. [PubMed] [PubMed]
Mewis, J. , and Wagner, N. J. , 2012, Colloidal Suspension Rheology, Cambridge University Press, Cambridge, UK.
Rao, L. G. , and Allison, J. T. , 2015, “ Generalized Viscoelastic Material Design With Integro-Differential Equations and Direct Optimal Control,” ASME Paper No. DETC2015-46768.


Grahic Jump Location
Fig. 2

The introduction of a viscoelastic connection can change the design space from a discrete arrangement of linear springs and dashpots (left) to a single viscoelastic element with a relaxation kernel, K(t)

Grahic Jump Location
Fig. 1

Function-valued viscoelastic properties may enable novel or improved performance (performance example schematics depict vibration damping of structures, vibration isolation of a suspended mass, and pressure sensitive adhesives). An early-stage design question is “what properties are optimal?” (represented here by the left, curved, downward arrow). This paper focuses on a methodology for answering that question with linear viscoelastic materials/systems including the nontrivial problem of selecting design-appropriate modeling using material properties (represented by the right, curved, upward arrow). Successful identification of targeted properties then leads to microstructure and formulation design (material selection or material synthesis) [2528], which is beyond the scope of the work here.

Grahic Jump Location
Fig. 3

Parameterizations of the relaxation kernel, K(t), for a general linear viscoelastic fluid model (Fig. 4(c)) into numerous linear viscoelastic models to reduce the complexity of the design space. The viscoelastic models used to parameterize the relaxation modulus of the added viscoelastic component are: (i) a linear dashpot (solid-line), (ii) a Maxwell element (linearspring and dashpot connected in series) single mode (dashed–dotted line), (iii) multimode Maxwell element (short dashed–dotted line), and (iv) a critical gel model, mechanical spring-dashpot analog is not applicable (short-dashed line). Inset plot is double-log plot of the same curves with the Maxwell model characteristic timescale λ (as in Eq. (13)) and the critical gel exponent n (as in Eq. (16)) labeled for reference.

Grahic Jump Location
Fig. 4

Design of optimal viscoelastic vibration isolation for a one-dimensional spring-mass system (a); (b) the typical approach of an arrangement of springs and dashpots, and (c) a generalized viscoelastic element with relaxation kernel, K(t). The latter approach increases design freedom and identifies more optimal targets.

Grahic Jump Location
Fig. 7

The design space for (a) the single mode Maxwell element and (b) the critical gel to isolate vibrations in a single mass-spring system (given by Eqs. (13) and (16)). Each model increases the design space to two dimensions, compared to the single dimension of a simple linear dashpot. Higher saturation correlates with lower peak normalized accelerations, as given in the scale bar (right), under simple sinusoidal displacement forcing. The white circle identifies a global minimum of peak acceleration.

Grahic Jump Location
Fig. 8

(a) Optimized performance (acceleration amplitude) and (b) the corresponding optimal viscoelastic (VE) design (kernel function). The viscoelastic models used to parameterize the relaxation kernel of the added viscoelastic component (shown in Fig. 4(c)) are: (i) a linear dashpot, solid, (ii) a single mode Maxwell element, dashed–dotted lines, (iii) three-mode Maxwell element, short dashed–dotted lines, and (iv) a critical gel model, short-dashed lines. A multimode Maxwell model uses its additional degrees of freedom to achieve approximately the same relaxation kernel and performance as a single mode, which therefore overlap as the most optimal of these parameterizations.

Grahic Jump Location
Fig. 5

Design tradeoff from added dashpot in Fig. 4(b) vibration isolator, displacement amplitude (|x|), and acceleration amplitude (|x¨|) responses. Decreasing peak amplitude at a resonant frequency lends to worse signal attenuation at higherfrequencies. Damping coefficients (ζ=c/2mωn) of ζ=0.1,0.5,0.9 are shown. Response of the mass is plotted with respect to ω̃, a frequency normalized by the resonant frequency.

Grahic Jump Location
Fig. 6

Design involving a generalized viscoelastic element with relaxation kernel K(t) can eliminate high-frequency acceleration. Here shown with Maxwell element (dashed–dotted line), and critical gel (short-dashed line) compared to dashpot (solid line), or no additional component (dashed line).

Grahic Jump Location
Fig. 9

The simple vibration isolator problem is extended to optimize the performance for a range of natural frequencies of the system but a fixed K(t). Increasing the natural frequency of the system (ω1=k/m) can be conceptualized by increasing the spring constant k and/or decreasing the system mass, m, shown in (a)–(c).

Grahic Jump Location
Fig. 10

The performance (left) and viscoelastic (VE) design (right) for a Maxwell element and Critical gel element in the extended vibration isolation problem. Multimode and single mode system response was identical, thus only a single mode is shown. The system performance is measured by the acceleration response while the viscoelastic design corresponds to the kernel function, K(t), of the viscoelastic element. Dashed lines represent optimization for a reference nondimensional natural frequency ω̃*=1, where ω̃* is as described in Eq. (34); solid lines represent the results optimized for the range of ω̃*. In the left panels (performance), the values of the objective function (fourth order norm of the maximum accelerations of each curve), optimized for a reference ω̃*=1 (dashed lines) and a range of ω̃* values are shown as horizontal lines. The optimizer works to minimize this objective function.



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