Research Articles

Stress Relaxation Near the Tip of a Stationary Mode I Crack in a Poroelastic Solid

[+] Author and Article Information
Chung-Yuen Hui

Field of Theoretical and Applied Mechanics,
Cornell University,
Ithaca, NY 14853

Rong Long

Department of Mechanical Engineering,
University of Colorado,
Boulder, CO 80303

Jing Ning

Department of Mechanical and Aerospace Engineering,
Cornell University,
Ithaca, NY 14853

Contributed by the Applied Mechanics Division of ASME for publication in the JOURNAL OF APPLIED MECHANICS. Manuscript received March 5, 2012; final manuscript received June 29, 2012; accepted manuscript posted July 25, 2012; published online January 22, 2013. Assoc. Editor: Chad Landis.

J. Appl. Mech 80(2), 021014 (Jan 22, 2013) (9 pages) Paper No: JAM-12-1093; doi: 10.1115/1.4007228 History: Received March 05, 2012; Revised June 29, 2012; Accepted July 25, 2012

We study the short time transient stress and pore pressure fields near the tip of a stationary crack when a sudden load is applied to a poroelastic solid. These fields are determined using a small scale “yielding” (SSY) analysis where the stress relaxation due to fluid flow is confined to a small region near the crack tip. They are found to exhibit the usual inverse square root singularity characteristic of cracks in linear elastic solids. Analysis shows that these fields are self-similar; the region of stress relaxation that propagates outward from the crack tip is proportional to Dct, where Dc is the cooperative diffusion coefficient and t is time. The pore pressure at the crack tip vanishes immediately after loading. The stress intensity factor at the crack tip is found to be reduced by a factor of 1/[2(1-v)], where v is the Poisson's ratio of the drained solid. Closed form approximations are found for the pore pressure and the trace of the effective stress. These approximate analytical solutions compare well with finite element results.

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


Riedel, H., and Rice, J. R., 1980, “Tensile Cracks in Creeping Solids”, ASTM STP 700, American Society for Testing and Materials, Philadelphia, pp. 112–130.
Tada, H., Paris, P. C., and Irwin, G. R., 2000, The Stress Analysis of Cracks Handbook, 3rd ed., ASME Press, New York.
Galli, M., Fornasiere, E., Cugnoni, J., and Oyen, M. L., 2011, “Poroviscoelastic Characterization of Particle-Reinforced Gelatin Gels Using Indentation and Homogenization,” J. Mech. Behav. Biomed. Mater., 4, pp. 610–617. [CrossRef] [PubMed]
Hu, Y., Zhao, X., Vlassak, J. J., and Suo, Z., 2010, “Using Indentation to Characterize the Poroelasticity of Gels,” Appl. Phys. Lett., 96, pp. 121904–121906. [CrossRef]
Hui, C. Y., Lin, Y. Y., Chuang, F. C., Shull, K. R., and Lin, W. C., 2006, “A Contact Mechanics Method for Characterizing the Elastic Properties and Permeability of Gels,” J. Polym. Sci., Part B: Polym. Phys., 44, pp. 359–370. [CrossRef]
Lin, Y. Y., and Hu, B. W., 2006, “Load Relaxation of a Flat Rigid Circular Indenter on a Gel Half Space,” J. Non-Cryst. Solids, 352, pp. 4034–4040. [CrossRef]
Scherer, G. W., 1989, “Drying Gels VIII. Revision and Review,” J. Non-Cryst. Solids, 109, pp. 171–182. [CrossRef]
Scherer, G. W., 1992, “Crack-Tip Stress in Gels,” J. Non-Cryst. Solids, 144, pp. 210–216. [CrossRef]
Ruina, A., 1978, “Influence of Coupled Deformation-Diffusion Effects On the Retardation of Hydraulic Fracture,” M.Sc. thesis, Brown University, Providence, RI.
Baumberger, T., Caroli, C., and Martina, D., 2006, “Solvent Control of Crack Dynamics in a Reversible Hydrogel,” Nature Mater., 5, pp. 552–555. [CrossRef]
Baumberger, T., Caroli, C., and Martina, D., 2006, “Fracture of a Biopolymer Gel as a Viscoplastic Disentanglement Process,” Eur. Phys. J. E, 21, pp. 81–89. [CrossRef]
Seitz, M. E., Martina, D., Baumberger, T., Krishnan, V. R., Hui, C. Y., and Shull, K. R., 2009, “Fracture and Large Strain Behavior of Self-Assembled Triblock Copolymer Gels,” Soft Matter, 5, pp. 447–456. [CrossRef]
Biot, M. A., 1941, “General Theory of Three-Dimensional Consolidation,” J. Appl. Phys., 12, pp. 155–164. [CrossRef]
Rice, J. R., and Cleary, M. P., 1976, “Some Basic Stress Diffusion Solutions for Fluid-Saturated Elastic Porous Media With Compressible Constituents,” Rev. Geophys. Space Phys., 14, pp. 227–241. [CrossRef]
Biot, M. A., 1955, “Theory of Elasticity and Consolidation for a Porous Anisotropic Solid,” J. Appl. Phys., 26, pp. 182–185. [CrossRef]
Knowles, J. K., and Sternberg, E., 1973, “An Asymptotic Finite-Deformation Analysis of the Elastostatic Field Near the Tip of a Crack,” J. Elast., 3, pp. 67–107. [CrossRef]
Stephenson, R., 1982, “The Equilibrium Field Near the Tip of a Crack for Finite Plane Strain of Incompressible Elastic Materials,” J. Elast., 12, pp. 65–99. [CrossRef]
Krishnan, V. R., Hui, C. Y., and Long, R., 2008, “Finite Strain Crack Tip Fields in Soft Incompressible Elastic Solids,” Langmuir, 24, pp. 14245–14253. [CrossRef] [PubMed]


Grahic Jump Location
Fig. 1

Schematic illustration of the SSY problem. The poro-elastic specimen on the left is loaded in Mode I tension. For short times, flow is confined in the region Ω which is very small compared with typical specimen dimensions. As a result, the crack can be modeled as semi-infinite with the elastic KI field prescribed at r=∞.

Grahic Jump Location
Fig. 2

(a) Contour plot of the pore pressure in a small region surrounding the crack tip at t=1.8 × 10-4. To avoid large numbers, the stress intensity factor is set to 1. Maximum pore pressure occurs at x=1.98Dct≈0.02 (Dc=1/1.8), y = 0. (b) Close up view of flow directions near maximum pressure. Arrows indicate flow direction.

Grahic Jump Location
Fig. 3

Finite element pore pressure distribution directly ahead of the crack tip (θ = 0) for v=0.4 at two different times

Grahic Jump Location
Fig. 4

Normalized pore pressure versus similarity variable r/(2Dct). Symbols are finite element results for v=0.4 at two different times. Solid line is obtained using Eq. (17).

Grahic Jump Location
Fig. 5

Normalized trace of effective stress directly ahead of the crack tip (θ = 0) versus similarity variable η/2 for two different Poisson ratios v=0.3,0.4 and at two different times t=10-4,10-5. Symbols are finite element results and the solid line is the analytical result given by Eq. (18).

Grahic Jump Location
Fig. 6

Crack opening displacement versus distance behind the crack tip r for two different Poisson's ratios v=0.3,0.2 and at two different times t=10-4,10-5. Symbols are finite element results and the solid line is Eq. (25).

Grahic Jump Location
Fig. 7

(a) Angular variations of u¯1 at r=δ for v=0.3. Dashed lines are finite element results, the solid line is obtained using Eq. (26a). (b) Angular variations of u¯2 at r=δ for v=0.3. The solid line is obtained using Eq. (26b) and the dashed line is the FEM result.

Grahic Jump Location
Fig. 8

Angular variation of p¯/erf(η/2) versus θ for v=0.3. Solid line is the analytical result Eq. (17).

Grahic Jump Location
Fig. 9

Angular variation of σ¯/erfc(η/2) versus θ for v=0.3. The dashed line is given by finite element results and the solid line is obtained using Eq. (18).

Grahic Jump Location
Fig. 10

Finite element meshes used in calculation



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