The purpose of the present study was to develop a computationally efficient finite element model that could be useful for parametric analysis of the biphasic poroviscoelastic (BPVE) behavior of articular cartilage under various loading conditions. The articular cartilage was modeled as the BPVE mixture of a porous, linear viscoelastic, and incompressible solid and an inviscid and incompressible fluid. A finite element (FE) formulation of the BPVE model was developed using two different algorithms, the continuous and discrete spectrum relaxation functions for the viscoelasticity of the solid matrix. These algorithms were applied to the creep and stress relaxation responses to the confined compression of articular cartilage, and a comparison of their performances was made. It was found that the discrete spectrum algorithm significantly saved CPU time and memory, as compared to the continuous spectrum algorithm. The consistency analysis for the present FE formulation was performed in comparison with the IMSL, a commercially available numerical software package. It was found that the present FE formulation yielded consistent results in predicting model behavior, whereas the IMSL subroutine produced inconsistent results in the velocity field, and thereby in the strain calculation.

1.
Buckwalter, J. A., Rosenberg, L. C., and Hunziker, E. B., 1990, “Articular cartilage: Composition, structure, response to injury and methods of facilitating repair,” in: Articular Cartilage and Knee Joint Function, Ewing, J. W., ed., Raven Press, New York, pp. 19–56.
2.
Cost
T. L.
, and
Becker
E. B.
,
1970
, “
A multidata method of approximate Laplace transform inversion
,”
Int, J. Num. Meth. Eng.
, Vol.
2
, pp.
207
219
.
3.
Crump
K. S.
,
1976
, “
Numerical inversion of Laplace transforms using a Fourier series approximation
,”
J. Assoc. Comp. Machinery
, Vol.
23
, pp.
89
96
.
4.
Eberliardt
A. W.
,
Keer
L. M.
,
Lewis
J. L.
, and
Vithoontien
V.
,
1990
, “
An analytical model of joint contact
,”
ASME JOURNAL OF BIOMECHANICAL ENGINEERING
, Vol.
112
, pp.
407
413
.
5.
Fung, Y. C., 1965, Foundations of Solid Mechanics, Prentice-Hall, Englewood Cliffs, NJ.
6.
Fung, Y. C., 1972, “Stress-strain-history relations of soft tissues in simple elongation,” in: Biomechanics: Its Foundations and Objectives, Fung, Y. C. Perrone, N., and Anliker, M., eds., Prentice-Hall, Englewood Cliffs, NJ, pp. 181–208.
7.
Fung, Y. C., 1994, Biomechanics: Mechanical Properties of Living Tissues, Springer-Verlag, New York.
8.
Hayes
W. C.
, and
Bodine
A. J.
,
1978
, “
Flow-independent viscoelastic properties of articular cartilage matrix
,”
J. Biomech.
, Vol.
11
, pp.
407
419
.
9.
Hodge
W. A.
,
Fijan
R. S.
,
Carlson
K. L.
,
Burgess
R. G.
,
Harris
W. H.
, and
Mann
R. W.
,
1986
, “
Contact pressures in the human hip joint measured in vivo
,”
Proc. Natl. Acad. Sci.
, Vol.
83
, pp.
2879
2883
.
10.
Holmes
M. H.
,
1986
, “
Finite deformation of soft tissue: Analysis of a mixture model in uni-axial compression
,”
ASME JOURNAL OF BIOMECHANICAL ENGINEERING
, Vol.
108
, pp.
372
381
.
11.
Kwan
M. K.
,
Lai
W. M.
, and
Mow
V. C.
,
1990
, “
A finite deformation theory for cartilage and other soft hydrated connective tissues—I. Equilibrium results
,”
J. Biomeck
, Vol.
23
(
2
), pp.
145
155
.
12.
Kwan
M. K.
,
Lin
T. H.-C.
, and
Woo
S. L.-Y.
,
1993
, “
On the viscoelastic properties of the anteromedial bundle of the anterior cruciate ligament
,”
J. Biomech.
, Vol.
26
, pp.
447
452
.
13.
Lai
W. M.
, and
Mow
V. C.
,
1980
, “
Drag-induced compression of articular cartilage during a permeation experiment
,”
Biorheology
, Vol.
17
, pp.
111
123
.
14.
Lai
W. M.
,
Mow
V. C.
, and
Roth
V.
,
1981
, “
Effects of nonlinear strain dependent permeability and rate of compression on the stress behavior of articular cartilage
,”
ASME JOURNAL OF BIOMECHANICAL ENGINEERING
, Vol.
103
, pp.
61
66
.
15.
Mak
A. F.
,
1986
a, “
The apparent viscoelastic behavior of articular cartilage—The contributions from the intrinsic matrix viscoelasticity and interstitial fluid flows
,”
ASME JOURNAL OF BIOMECHANICAL ENGINEERING
, Vol.
108
, pp.
123
130
.
16.
Mak
A. F.
,
1986
b, “
Unconfined compression of hydrated viscoelastic tissues: A biphasic poroviscoelastic analysis
,”
Biorheology
, Vol.
23
, pp.
371
383
.
17.
Math/Library, 1992, “FORTRAN Subroutines for Mathematical Applications,” IMSL, Sugar Land, TX.
18.
Mow
V. C.
,
Kuei
S. C.
,
Lai
W. M.
, and
Armstrong
C. G.
,
1980
, “
Biphasic creep and stress relaxation of articular cartilage in compression: theory and experiments
,”
ASME JOURNAL OF BIOMECHANICAL ENGINEERING
, Vol.
102
, pp.
73
83
.
19.
Mow
V. C.
, and
Lai
W. M.
,
1980
, “
Recent development in synovial joint biomechanics
,”
SIAM Review
, Vol.
22
(
3
), pp.
275
317
.
20.
Mow
V. C.
,
Mak
A. F.
,
Lai
W. M.
,
Rosenberg
L. C.
, and
Tang
L. H.
,
1984
a, “
Viscoelastic properties of proteoglycan subunits and aggregates in varying solution concentrations
,”
J. Biomech
, Vol.
17
, pp.
325
338
.
21.
Mow
V. C.
,
Holmes
M. H.
, and
Lai
W. M.
,
1984
b, “
Fluid transport and mechanical properties of articular cartilage: a review
,”
J. Biomech.
, Vol.
17
, pp.
377
394
.
22.
Mow
V. C.
,
Ateshian
G. A.
, and
Spilker
R. L.
,
1993
, “
Biomechanics of diarthrodial joints: A review of twenty years of progress
,”
ASME JOURNAL OF BIOMECHANICAL ENGINEERING
, Vol.
115
, pp.
460
467
.
23.
Myers
B. S.
,
McElhaney
J. H.
, and
Doherty
B. J.
,
1991
, “
The viscoelastic responses of the human cervical spine in torsion: Experimental limitations of quasilinear theory, and a method for reducing these effects
,”
J. Biomeck
, Vol.
24
, pp.
811
817
.
24.
Puso, M. A., Weiss, J. A., and Schauer, D. A., 1995, “Finite element implementation for anisotropic quasilinear viscoelasticity,” Advances in Bioengineering, ASME, pp. 57–58.
25.
Setton, L. A., and Mow, V. C., 1992, “A generalized biphasic poroviscoelastic model for articular cartilage: Theory and experiment,” Advances in Bioengineering, ASME, pp. 589–592.
26.
Setton
L. A.
,
Zhu
W.
, and
Mow
V. C.
,
1993
, “
The biphasic poroviscoelastic behavior of articular cartilage: Role of the surface zone in governing the compressive behavior
,”
J. Biomeck
, Vol.
26
, pp.
581
592
.
27.
Spilker
R. L.
, and
Suh
J.-K.
,
1990
, “
Formulation and evaluation of a finite element model for the biphasic model of hydrated soft tissue
,”
Comp. Struc.
, Vol.
35
(
4
), pp.
425
439
.
28.
Spirt
A. A.
,
Mak
A. F.
, and
Wassell
R. P.
,
1989
, “
Nonlinear viscoelastic properties of articular cartilage in shear
,”
J. Orthop. Res.
, Vol.
7
, pp.
43
49
.
29.
Suh
J. K.
,
Spilker
R. L.
, and
Holmes
M. H.
,
1991
, “
A penalty finite element analysis for nonlinear mechanics of biphasic hydrated soft tissue under large deformation
,”
Int. J. Num. Metk. Eng.
, Vol.
32
, pp.
1411
1439
.
30.
Tschoegl, N. W., 1989, The Phenomenological Theory of Linear Viscoelastic Behavior: An Introduction, Springer-Verlag, Berlin.
31.
Viidik
A.
,
1968
, “
A rheological model for uncalcified parallel-fibered collagenous tissues
,”
J. Biomech
, Vol.
1
, pp.
3
11
.
32.
Zhu, W., and Mow, V. C., 1990, “Viscometric properfies of proteoglycan solutions at physiological concentration,” in: Biomechanics of Dialhrodial Joints, Mow, V. C., Ratcliffe, A., and Woo, S. L.-Y., eds., Springer-Verlag, New York, pp. 313–344.
33.
Zhu
W.
,
Mow
V. C.
,
Koob
T. J.
, and
Eyre
D. R.
,
1993
, “
Viscoelastic shear properties of articular cartilage and the effects of glycosidase treatments
,”
J. Orthop. Res.
, Vol.
11
, pp.
771
781
.
This content is only available via PDF.
You do not currently have access to this content.