Research Papers

Determination of the Reduced Creep Function of Viscoelastic Compliant Materials Using Pipette Aspiration Method

[+] Author and Article Information
Yan-Ping Cao

AML, Institute of Biomechanics and Medical
Department of Engineering Mechanics,
Tsinghua University,
Beijing 100084, China
e-mail: caoyanping@tsinghua.edu.cn

Guo-Yang Li, Man-Gong Zhang, Xi-Qiao Feng

AML, Institute of Biomechanics and Medical
Department of Engineering Mechanics,
Tsinghua University,
Beijing 100084, China

1Corresponding author.

Manuscript received February 14, 2014; final manuscript received March 7, 2014; accepted manuscript posted March 13, 2014; published online April 1, 2014. Editor: Yonggang Huang.

J. Appl. Mech 81(7), 071006 (Apr 01, 2014) (7 pages) Paper No: JAM-14-1068; doi: 10.1115/1.4027159 History: Received February 14, 2014; Revised March 07, 2014; Accepted March 13, 2014

Determining the mechanical properties of soft matter across different length scales is of great importance in understanding the deformation behavior of compliant materials under various stimuli. A pipette aspiration test is a promising tool for such a purpose. A key challenge in the use of this method is to develop explicit expressions of the relationship between experimental responses and material properties particularly when the tested sample has irregular geometry. A simple scaling relation between the reduced creep function and the aspiration length is revealed in this paper by performing a theoretical analysis on the aspiration creep tests of viscoelastic soft solids with arbitrary surface profile. Numerical experiments have been performed on the tested materials with different geometries to validate the theoretical solution. In order to incorporate the effects of the rise time of the creep pressure, an analytical solution is further derived based on the generalized Maxwell model, which relates the parameters in reduced creep function to the aspiration length. Its usefulness is demonstrated through a numerical example and the analysis of the experimental data from literature. The analytical solutions reported here proved to be independent of the geometric parameters of the system under described conditions. Therefore, they may not only provide insight into the deformation behavior of soft materials in aspiration creep tests but also facilitate the use of this testing method to deduce the intrinsic creep/relaxation properties of viscoelastic compliant materials.

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


Ahn, B., and Kim, J., 2010, “Measurement and Characterization of Soft Tissue Behavior With Surface Deformation and Force Response Under Large Deformations,” Med. Image Anal., 14(2), pp. 138–148. [CrossRef]
Wang, M., 2003, “Developing Bioactive Composite Materials for Tissue Replacement,” Biomaterials, 24(13), pp. 2133–2151. [CrossRef]
Desprat, N., Richert, A., Simeon, J., and Asnacios, A., 2005, “Creep Function of a Single Living Cell,” Biophys. J., 88(3), pp. 2224–2233. [CrossRef]
Discher, D. E., Janmey, P., and Wang, Y. L., 2005, “Tissue Cells Feel and Respond to the Stiffness of Their Substrate,” Science, 310(5751), pp. 1139–1143. [CrossRef]
Butcher, J. T., McQuinn, T. C., Sedmera, D., Turner, D., and Markwald, R. R., 2007, “Transitions in Early Embryonic Atrioventricular Valvular Function Correspond With Changes in Cushion Biomechanics That are Predictable by Tissue Composition,” Circ. Res., 100(10), pp. 1503–1511. [CrossRef]
Vaziri, A., and Gopinath, A., 2008, “Cell and Biomolecular Mechanics in Silico,” Nat. Mater., 7(1), pp. 15–23. [CrossRef]
Lv, S., Dudek, D. M., Cao, Y., Balamurali, M. M., Gosline, J., and Li, H., 2010, “Designed Biomaterials to Mimic the Mechanical Properties of Muscles,” Nature, 465(7294), pp. 69–73. [CrossRef]
Li, B., Cao, Y. P., Feng, X. Q., and Gao, H., 2012, “Mechanics of Morphological Instabilities and Surface Wrinkling in Soft Materials: A Review,” Soft Matter, 8(21), pp. 5728–5745. [CrossRef]
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]
Kranenburg, J. M., Tweedie, C. A., van Vliet, K. J., and Schubert, U. S., 2009, “Challenges and Progress in High—Throughput Screening of Polymer Mechanical Properties by Indentation,” Adv. Mater., 21(35), pp. 3551–3561. [CrossRef]
Oyen, M. L., 2007, “Sensitivity of Polymer Nanoindentation Creep Measurements to Experimental Variables,” Acta Mater., 55(11), pp. 3633–3639. [CrossRef]
Smith, S. B., Finzi, L., and Bustamante, C., 1992, “Direct Mechanical Measurements of the Elasticity of Single DNA Molecules by Using Magnetic Beads,” Science, 258(5085), pp. 1122–1126. [CrossRef]
Ziemann, F., Rädler, J., and Sackmann, E., 1994, “Local Measurements of Viscoelastic Moduli of Entangled Actin Networks Using an Oscillating Magnetic Bead Micro-Rheometer,” Biophys. J., 66(6), pp. 2210–2216. [CrossRef]
Sato, M., Theret, D. P., Ohshima, N., Nerem, R. M., and Wheeler, L. T., 1990, “Application of the Micropipette Technique to the Measurement of Cultured Porcine Aortic Endothelial Cell Viscoelastic Properties,” ASME J. Biomech. Eng., 112(3), pp. 263–268. [CrossRef]
Discher, D. E., Mohandas, N., and Evans, E. A., 1994, “Molecular Maps of Red Cell Deformation: Hidden Elasticity and In Situ Connectivity,” Science, 266(5187), pp. 1032–1035. [CrossRef]
Aoki, T., Ohashi, T., Matsumoto, T., and Sato, M., 1997, “The Pipette Aspiration Applied to the Local Stiffness Measurement of Soft Tissues,” Ann. Biomed. Eng., 25(3), pp. 581–587. [CrossRef]
Hochmuth, R. M., 2000, “Micropipette Aspiration of Living Cells,” J. Biomech., 33(1), pp. 15–22. [CrossRef]
Trickey, W. R., Lee, G. M., and Guilak, F., 2000, “Viscoelastic Properties of Chondrocytes From Normal and Osteoarthritic Human Cartilage,” J. Orthop. Res., 18(6), pp. 891–898. [CrossRef]
Zhou, E. H., Lim, C. T., and Quek, S. T., 2005, “Finite Element Simulation of the Micropipette Aspiration of a Living Cell Undergoing Large Viscoelastic Deformation,” Mech. Adv. Mater. Struct., 12(6), pp. 501–512. [CrossRef]
Merryman, W. D., Guilak, F., Sacks, M. S., and Bieniek, P. D., 2009, “Viscoelastic Properties of the Aortic Valve Interstitial Cell,” ASME J. Biomech. Eng., 131(4), p. 041005. [CrossRef]
Guevorkian, K., Colbert, M. J., Durth, M., Dufour, S., and Brochard-Wyart, F.2010, “Aspiration of Biological Viscoelastic Drops,” Phys. Rev. Lett., 104(21), p. 218101. [CrossRef]
Zhao, R., Sider, K. L., and Simmons, C. A., 2011, “Measurement of Layer-Specific Mechanical Properties in Multilayered Biomaterials by Micropipette Aspiration,” Acta Biomater., 7(3), pp. 1220–1227. [CrossRef]
Zhou, E. H., Xu, F., Quek, S. T., and Lim, C. T., 2012, “A Power-Law Rheology-Based Finite Element Model for Single Cell Deformation,” Biomech. Model Mech., 11(7), pp. 1075–1084. [CrossRef]
Zhao, R., Wyss, K., and Simmons, C. A., 2009, “Comparison of Analytical and Inverse Finite Element Approaches to Estimate Cell Viscoelastic Properties by Micropipette Aspiration,” J. Biomech., 42(16), pp. 2768–2773. [CrossRef]
Cao, Y. P., Ji, X. Y., and Feng, X. Q., 2010, “Geometry Independence of the Normalized Relaxation Functions of Viscoelastic Materials in Indentation,” Philos. Mag., 90(12), pp. 1639–1655. [CrossRef]
Cao, Y. P., Zhang, M. G., and Feng, X. Q., 2013, “Indentation Method for Measuring the Viscoelastic Kernel Function of Nonlinear Viscoelastic Soft Materials,” J. Mater. Res., 28(6), pp. 806–816. [CrossRef]
Baaijens, F. P., Trickey, W. R., Laursen, T. A., and Guilak, F., 2005, “Large Deformation Finite Element Analysis of Micropipette Aspiration to Determine the Mechanical Properties of the Chondrocyte,” Ann. Biomed. Eng., 33(4), pp. 494–501. [CrossRef]
Lakes, R. S., 2009, Viscoelastic Materials, Cambridge University Press, New York.
Barenblatt, G. I., 1996, Scaling, Self-Similarity, and Intermediate Asymptotics: Dimensional Analysis and Intermediate Asymptotics, Cambridge University Press, New York.
Pesic, P., 2005, Sky in a Bottle, MIT, Cambridge, MA.
Christensen, R. M., 1982, Theory of Viscoelasticity: An Introduction, Academic, New York.
DS Simulia Corp., 2010, “ABAQUS User's Manual,” Version 6.10, Dassault Systemes, Providence, RI.
Simo, J. C., 1987, “On a Fully Three-Dimensional Finite-Strain Viscoelastic Damage Model: Formulation and Computational Aspects,” Comput. Meth. Appl. Mech. Eng., 60(2), pp. 153–173. [CrossRef]
Mooney, M., 1940, “A Theory of Large Elastic Deformation,” J. Appl. Phys., 11(9), pp. 582–592. [CrossRef]
Rivlin, R. S., 1948, “Large Elastic Deformations of Isotropic Materials. IV. Further Developments of the General Theory,” Philos. Trans. R. Soc., A, 241(835), pp. 379–397. [CrossRef]
Fung, Y. C., Fronek, K., and Patitucci, P., 1979, “Pseudoelasticity of Arteries and the Choice of Its Mathematical Expression,” Am. J. Physiol., 237(5), pp. H620–H631.
Arruda, E. M., and Boyce, M. C., 1993, “A Three-Dimensional Constitutive Model for the Large Stretch Behavior of Rubber Elastic Materials,” J. Mech. Phys. Solids, 41(2), pp. 389–412. [CrossRef]
Zhang, M. G., Cao, Y. P., Li, G. Y., and Feng, X. Q., “Pipette Aspiration of Hyperelastic Compliant Materials: Theoretical Analysis, Simulations, and Experiments” (unpublished).
Chien, W. Z., 1947, “Large Deflection of a Circular Clamped Plate Under Uniform Pressure,” Acta Phys. Sin., 7(2), pp. 102–107.
Wang, Q. M., Mohan, A. C., Oyen, M. L., and Zhao, X. H., 2014, “Separating Viscoelasticity and Poroelasticity of Gels With Different Length and Time Scales,” Acta Mech. Sin., 30(1), pp. 20–27. [CrossRef]


Grahic Jump Location
Fig. 1

(a) Schematic plot of the creep function; (b) illustration of the aspiration test, negative pressure is imposed to the surface of a tested material through a cylindrical tube with the inner radius R0 and outer radius Re

Grahic Jump Location
Fig. 4

Numerical experiments used to illustrate the usefulness of the analytical solution given by Eq. (10). (a), (c), and (e) represent the finite element models for examples A, B, and C, respectively. A1 and A2 in (a) and (c) are fixed to eliminate the rigid-body motion of the tested materials, other points can be selected which has on effects on the identified results. (b), (d), and (f) are the deformed configurations. The displacements at Bi(i=1,2,3) are recorded as the aspiration lengths.

Grahic Jump Location
Fig. 2

Finite element models used to examine the effects of the variation in the contact area on the identified reduced creep function using Eq. (10). (a) and (b) represent the computational model and the deformation of the tested material for case I; (c) and (d) represent the computational model and the deformation of the tested material for case II; (e) and (f) are the computational model and the deformed configuration of the tested material for case III.

Grahic Jump Location
Fig. 3

A comparison of the identified reduced creep functions using Eq. (10) with the actual solution input in the finite element simulations:  J˜ (t) = 1 + (1 − exp (− t /10)) for different cases

Grahic Jump Location
Fig. 7

(a) Illustration of the loading procedure in an aspiration creep test, where tr is the time at which the negative pressure reach the maximum value and then is kept as constant; (b) schematic plot of the generalized Maxwell model, where Ke,K1,K2,…,Kj are spring constants and η1,η2,…,ηj are viscosity coefficients; (c) normalized aspiration lengths given by finite element simulations and those predicted using Eqs. (10) and (13)

Grahic Jump Location
Fig. 8

Experimental data used to illustrate the difference between the result identified using Eq. (10) and that given by Eq. (13). (a) Variation of the aspiration length with time given by the experiments of Ref. [20]; (b) normalized aspiration length l˜(t) at the creep stage.

Grahic Jump Location
Fig. 5

Pipette aspiration test of soft layer bonded to rigid substrate. (a) Finite element model; (b) deformed configuration of the tested material, the displacement at the center point of the upper surface is recorded as the aspiration length.

Grahic Jump Location
Fig. 6

A comparison of the identified reduced creep functions using Eq. (10) and the aspiration lengths given by finite element simulations with the actual solution. The actual solution is taken here as J˜(t) = 1 + 0.1(1 − exp (−t)) + 0.9(1 − exp (−t/18)). Here, the actual solution represents the reduced creep function input in the finite element simulations.



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