Finite element simulation of location- and time-dependent mechanical behavior of chondrocytes in unconfined compression tests.

Experimental evidence suggests that cells are extremely sensitive to their mechanical environment and react directly to mechanical stimuli. At present, it is technically difficult to measure fluid pressure, stress, and strain in cells, and to determine the time-dependent deformation of chondrocytes. For this reason, there are no data in the published literature that show the dynamic behavior of chondrocytes in articular cartilage. Similarly, the dynamic chondrocyte mechanics have not been calculated using theoretical models that account for the influence of cell volumetric fraction on cartilage mechanical properties. In the present investigation, the location- and time-dependent stress-strain state and fluid pressure distribution in chondrocytes in unconfined compression tests were simulated numerically using a finite element method. The technique involved two basic steps: first, cartilage was approximated as a macroscopically homogenized material and the mechanical behavior of cartilage was obtained using the homogenized model; second, the solution of the time-dependent displacements and fluid pressure fields of the homogenized model was used as the time-dependent boundary conditions for a microscopic submodel to obtain average location- and time-dependent mechanical behavior of cells. Cells and extracellular matrix were assumed to be biphasic materials composed of a fluid phase and a hyperelastic solid phase. The hydraulic permeability was assumed to be deformation dependent and the analysis was performed using a finite deformation approach. Numerical tests were made using configurations similar to those of experiments described in the literature. Our simulations show that the mechanical response of chondrocytes to cartilage loading depends on time, fluid boundary conditions, and the locations of the cells within the specimen. The present results are the first to suggest that chondrocyte deformation in a stress-relaxation type test may exceed the imposed system deformation by a factor of 3-4, that chondrocyte deformations are highly dynamic and do not reach a steady state within about 20 min of steady compression (in an unconfined test), and that cell deformations are very much location dependent.