Influence of physicochemical factors on rheology of human neutrophils.

The effects of variations in temperature, pH, and osmolality on the rheological properties of human neutrophils were determined by studying the cell deformation in response to aspirational pressure applied via a micropipette. The time history of the deformation was analyzed by the use of a standard solid viscoelastic model consisting of an elastic element K1 in parallel with a Maxwell element (an elastic element K2 in series with a viscous element mu). With changes in temperature over a range of 9-40 degrees C, only mu varied inversely with temperature, while K1 and K2 did not show significant alterations. Variations in pH over the range of 5.4-7.8 did not significantly affect the viscoelastic coefficients, but K1 and mu rose at pH 8.4. An increase in osmolality caused all three coefficients to rise, but a decrease in osmolality had relatively little effect on the coefficients. These changes in response to physicochemical variations serve to provide insights into the viscoelastic properties of neutrophils and their possible roles in health and disease.     

Comparison of the viscoelastic properties of normal hepatocytes and hepatocellular carcinoma cells under cytoskeletal perturbation

The viscoelastic properties of both hepatocytes and hepatocellular carcinoma (HCC) cells were measured by means of a micropipette aspiration technique. Experimental results were analyzed with a three-element standard linear solid model, in which an elastic element, K1, is in parallel with a Maxwell element composed of another elastic element, K2, in series with a viscous element, mu. Further, we investigated the relevance of viscoelastic properties of these two types of cells to the cytoskeleton structures by treating cells with three cytoskeletal perturbing agents, namely cytochalasin D (CD), colchicine (Col) and vinblastine (VBL). The results showed that the elastic coefficients, but not viscous coefficient of HCC cells (K1 = 103.6 +/- 12.6 N m-2, K2 = 42.5 +/- 10.4 N m-2, mu = 4.5 +/- 1.9 Pa s, n = 30), were significantly higher than the corresponding values for hepatocytes (K1 = 87.5 +/- 12.1 N m-2, K2 = 33.3 +/- 10.3 N m-2, mu = 5.9 +/- 3.0 Pa s, n = 24). Upon treatment with CD, the viscoelastic coefficients of both hepatocytes and HCC cells decreased uniformly, with magnitudes for the decrease in elastic coefficients of HCC cells (K1: 68.7 to 81.7 N m-2, 66.3 to 78.9%; K2: 34.5 to 37.1 N m-2, 81.2 to 87.3%) larger than those for normal hepatocytes (K1: 42.6 to 49.8 N m-2, 48.7 to 56.9%; K2: 17.2 to 20.4 N m-2, 51.7 to 61.3%). There was a smaller decrease in the viscous coefficient of HCC cells (2.0 to 3.4 Pa s, 44.4 to 75.6%) than that for hepatocytes (3.0 to 3.9 Pa s, 50.8 to 66.1%). Upon treatment with Col and VBL, the elastic coefficients of hepatocytes generally increased or tended to increase while those of HCC cells decreased. The differences in either the pattern or the magnitude of the effect of cytoskeletal perturbing agent on the viscoelastic properties between HCC cells and hepatocytes might possibly reflect differences in the state of the cytoskeleton structure and function, or in the cells' sensitivity to perturbing agent treatment between these two types of cells. Changes in the viscoelastic properties of cancer cells might well affect tumor cell invasion and metastasis as well as interactions between tumor cells and their micro-mechanical environments.     

An orthotropic viscoelastic model for the passive myocardium: continuum basis and numerical treatment

This study deals with the viscoelastic constitutive modeling and the respective computational analysis of the human passive myocardium. We start by recapitulating the locally orthotropic inner structure of the human myocardial tissue and model the mechanical response through invariants and structure tensors associated with three orthonormal basis vectors. In accordance with recent experimental findings the ventricular myocardial tissue is assumed to be incompressible, thick-walled, orthotropic and viscoelastic. In particular, one spring element coupled with Maxwell elements in parallel endows the model with viscoelastic features such that four dashpots describe the viscous response due to matrix, fiber, sheet and fiber-sheet fragments. In order to alleviate the numerical obstacles, the strictly incompressible model is altered by decomposing the free-energy function into volumetric-isochoric elastic and isochoric-viscoelastic parts along with the multiplicative split of the deformation gradient which enables the three-field mixed finite element method. The crucial aspect of the viscoelastic formulation is linked to the rate equations of the viscous overstresses resulting from a 3-D analogy of a generalized 1-D Maxwell model. We provide algorithmic updates for second Piola–Kirchhoff stress and elasticity tensors. In the sequel, we address some numerical aspects of the constitutive model by applying it to elastic, cyclic and relaxation test data obtained from biaxial extension and triaxial shear tests whereby we assess the fitting capacity of the model. With the tissue parameters identified, we conduct (elastic and viscoelastic) finite element simulations for an ellipsoidal geometry retrieved from a human specimen.     

On the measurement of shear elastic moduli and viscosities of erythrocyte plasma membranes by transient deformation in high frequency electric fields.

We present a new method to measure the shear elastic moduli and viscosities of erythrocyte membranes which is based on the fixation and transient deformation of cells in a high-frequency electric field. A frequency domain of constant force (arising by Maxwell Wagner polarization) is selected to minimize dissipative effects. The electric force is thus calculated by electrostatic principles by considering the cell as a conducting body in a dielectric fluid and neglecting membrane polarization effects. The elongation A of the cells perpendicular to their rotational axis exhibits a linear regime (A proportional to Maxwell tension or to square of the electric field E2) at small, and a nonlinear regime (A proportional to square root of Maxwell tension or to the electric field E) at large extensions with a cross-over at A approximately 0.5 micron. The nonlinearity leads to amplitude-dependent response times and to differences of the viscoelastic response and relaxation functions. The cells exhibit pronounced yet completely reversible tip formations at large extensions. Absolute values of the shear elastic modulus, mu, and membrane viscosity, eta, are determined by assuming that field-induced stretching of the biconcave cell may be approximately described in terms of a sphere to ellipsoid deformation. The (nonlinear) elongation-vs.-force relationship calculated by the elastic theory of shells agress well with the experimentally observed curves and the values of mu = 6.1 x 10(-6) N/m and eta = 3.4 x 10(-7) Ns/m are in good agreement with the micropipette results of Evans and co-workers. The effect of physical, biochemical, and disease-induced structural changes on the viscoelastic parameters is studied. The variability of mu and eta of a cell population of a healthy donor is +/- 45%, which is mainly due to differences in the cell age. The average mu value of cells of different healthy donors scatters by +/- 18%. Osmotic deflation of the cells leads to a fivefold increase of mu and 10-fold increase of eta at 500 mosm. The shear modulus mu increases with temperature showing that the cytoskeleton does not behave as a network of entropy elastic springs. Elliptic cells of patients suffering from elliptocytosis of the Leach phenotype exhibit a threefold larger value of mu than normal discocytes of control donors. Cross-linking of the spectrin by the divalent S-H agents diamide (1 mM, 15 min incubation) leads to an eightfold increase of mu whereas eta is essentially constant. The effect of diamide is reversed after treatment with S-S bond splitting agents.     

An elastic network model based on the structure of the red blood cell membrane skeleton.

A finite element network model has been developed to predict the macroscopic elastic shear modulus and the area expansion modulus of the red blood cell (RBC) membrane skeleton on the basis of its microstructure. The topological organization of connections between spectrin molecules is represented by the edges of a random Delaunay triangulation, and the elasticity of an individual spectrin molecule is represented by the spring constant, K, for a linear spring element. The model network is subjected to deformations by prescribing nodal displacements on the boundary. The positions of internal nodes are computed by the finite element program. The average response of the network is used to compute the shear modulus (mu) and area expansion modulus (kappa) for the corresponding effective continuum. For networks with a moderate degree of randomness, this model predicts mu/K = 0.45 and kappa/K = 0.90 in small deformations. These results are consistent with previous computational models and experimental estimates of the ratio mu/kappa. This model also predicts that the elastic moduli vary by 20% or more in networks with varying degrees of randomness. In large deformations, mu increases as a cubic function of the extension ratio lambda 1, with mu/K = 0.62 when lambda 1 = 1.5.     

Role of physical properties of resistance vessel wall in regulating peripheral blood flow--II. Structural and mechanical behavior

In order to examine the structural and mechanical properties of the vessel wall resistance when subjected to autoregulatory flow control, a mechanical model for the vascular wall was derived from a mathematical model. The mechanical model was an analogue model which connected in series the Maxwell model (elastic modulus: K3) with the parallel elements of Hill's model (elastic modules: K2) and Hooke's elastic model (elastic modulus: K1); it was also mathematically equivalent to the Spring model (see part I). The structural and mechanical properties of the resistance vessel wall were characterized by the three elastic moduli (K1, alpha K2 and K3) [mmHg]. The parameter alpha was a modification factor of the elastic modulus K2 given by the myogenic mechanism. After a numerical analysis of the experimental data given by the mechanical model, we confirmed that the arterial pressure range for autoregulatory flow controls shifted to the upper region with an increase of the elastic modulus K1 and the flow regulation was reduced.     

Comparative study of viscoelastic arterial wall models in nonlinear one-dimensional finite element simulations of blood flow

It is well known that blood vessels exhibit viscoelastic properties, which are modeled in the literature with different mathematical forms and experimental bases. The wide range of existing viscoelastic wall models may produce significantly different blood flow, pressure, and vessel deformation solutions in cardiovascular simulations. In this paper, we present a novel comparative study of two different viscoelastic wall models in nonlinear one-dimensional (1D) simulations of blood flow. The viscoelastic models are from papers by Holenstein et al. in 1980 (model V1) and Valdez-Jasso et al. in 2009 (model V2). The static elastic or zero-frequency responses of both models are chosen to be identical. The nonlinear 1D blood flow equations incorporating wall viscoelasticity are solved using a space-time finite element method and the implementation is verified with the Method of Manufactured Solutions. Simulation results using models V1, V2 and the common static elastic model are compared in three application examples: (i) wave propagation study in an idealized vessel with reflection-free outflow boundary condition; (ii) carotid artery model with nonperiodic boundary conditions; and (iii) subject-specific abdominal aorta model under rest and simulated lower limb exercise conditions. In the wave propagation study the damping and wave speed were largest for model V2 and lowest for the elastic model. In the carotid and abdominal aorta studies the most significant differences between wall models were observed in the hysteresis (pressure-area) loops, which were larger for V2 than V1, indicating that V2 is a more dissipative model. The cross-sectional area oscillations over the cardiac cycle were smaller for the viscoelastic models compared to the elastic model. In the abdominal aorta study, differences between constitutive models were more pronounced under exercise conditions than at rest. Inlet pressure pulse for model V1 was larger than the pulse for V2 and the elastic model in the exercise case. In this paper, we have successfully implemented and verified two viscoelastic wall models in a nonlinear 1D finite element blood flow solver and analyzed differences between these models in various idealized and physiological simulations, including exercise. The computational model of blood flow presented here can be utilized in further studies of the cardiovascular system incorporating viscoelastic wall properties.     

Leukocyte deformability: finite element modeling of large viscoelastic deformation.

An axisymmetric deformation of a viscoelastic sphere bounded by a prestressed elastic thin shell in response to external pressure is studied by a finite element method. The research is motivated by the need for understanding the passive behavior of human leukocytes (white blood cells) and interpreting extensive experimental data in terms of the mechanical properties. The cell at rest is modeled as a sphere consisting of a cortical prestressed shell with incompressible Maxwell fluid interior. A large-strain deformation theory is developed based on the proposed model. General non-linear, large strain constitutive relations for the cortical shell are derived by neglecting the bending stiffness. A representation of the constitutive equations in the form of an integral of strain history for the incompressible Maxwell interior is used in the formulation of numerical scheme. A finite element program is developed, in which a sliding boundary condition is imposed on all contact surfaces. The mathematical model developed is applied to evaluate experimental data of pipette tests and observations of blood flow.     

Passive transverse mechanical properties of skeletal muscle under in vivo compression

The objective of the present study is to determine the passive transverse mechanical properties of skeletal muscle. Compression experiments were performed on four rat tibialis anterior muscles. To assess the stress- and strain-distributions in the muscle during the experiment, a plane stress model of the cross section was developed for each muscle. The incompressible viscoelastic Ogden model was used to describe the passive muscle behaviour. The four material parameters were determined by fitting calculated indentation forces on measured indentation forces. The elastic parameters, mu and alpha, were 15.6+/-5.4 kPa and 21.4+/-5.7, respectively. The viscoelastic parameters, delta and tau, were 0.549+/-0.056 and 6.01+/-0.42 s. When applying the estimated material parameters in a three-dimensional finite element model, the measured behaviour can be accurately simulated.     

Local measurements of viscoelastic parameters of adherent cell surfaces by magnetic bead microrheometry.

A magnetic bead microrheometer has been designed which allows the generation of forces up to 10(4) pN on 4.5 micron paramagnetic beads. It is applied to measure local viscoelastic properties of the surface of adhering fibroblasts. Creep response and relaxation curves evoked by tangential force pulses of 500-2500 pN (and approximately 1 s duration) on the magnetic beads fixed to the integrin receptors of the cell membrane are recorded by particle tracking. Linear three-phasic creep responses consisting of an elastic deflection, a stress relaxation, and a viscous flow are established. The viscoelastic response curves are analyzed in terms of a series arrangement of a dashpot and a Voigt body, which allows characterization of the viscoelastic behavior of the adhering cell surface in terms of three parameters: an effective elastic constant, a viscosity, and a relaxation time. The displacement field generated by the local tangential forces on the cell surface is visualized by observing the induced motion of assemblies of nonmagnetic colloidal probes fixed to the membrane. It is found that the displacement field decays rapidly with the distance from the magnetic bead. A cutoff radius of Rc approximately 7 micron of the screened elastic field is established. Partial penetration of the shear field into the cytoplasm is established by observing the induced deflection of intracellular compartments. The cell membrane was modeled as a thin elastic plate of shear modulus mu * coupled to a viscoelastic layer, which is fixed to a solid support on the opposite side; the former accounts for the membrane/actin cortex, and the latter for the contribution of the cytoskeleton to the deformation of the cell envelope. It is characterized by the coupling constant chi characterizing the elasticity of the cytoskeleton. The coupling constant chi and the surface shear modulus mu * are obtained from the measured displacements of the magnetic and nonmagnetic beads. By analyzing the experimental data in terms of this model a surface shear modulus of mu * approximately 2 . 10(-3) Pa m to 4 . 10(-3) Pa m is found. By assuming an approximate plate thickness of 0.1 micron one estimates an average bulk shear modulus of mu approximately (2 / 4) . 10(-4) Pa, which is in reasonable agreement with data obtained by atomic force microscopy. The viscosity of the dashpot is related to the apparent viscosity of the cytoplasm, which is obtained by assuming that the top membrane is coupled to the bottom (fixed) membrane by a viscous medium. By application of the theory of diffusion of membrane proteins in supported membranes we find a coefficient of friction of bc approximately 2 . 10(9) Pa s/m corresponding to a cytoplasmic viscosity of 2 . 10(3) Pa s.     

Local rheology of human neutrophils investigated using atomic force microscopy

During the immune response, neutrophils display localized mechanical events by interacting with their environment through the micro-vascular transit, trans-endothelial, and trans-epithelial migration. Nano-mechanical studies of human neutrophils on localized nano-domains could provide the essential information for understanding their immune responsive functions. Using the Atomic Force Microscopy (AFM)-based micro-rheology, we have investigated rheological properties of the adherent human neutrophils on local nano-domains. We have applied the modified Hertz model to obtain the viscoelastic moduli from the relatively thick body regions of the neutrophils. In addition, by using more advanced models to account for the substrate effects, we have successfully characterized the rheological properties of the thin leading and tail regions as well. We found a regional difference in the mechanical compliances of the adherent neutrophils. The central regions of neutrophils were significantly stiffer (1,548 ± 871 Pa) than the regions closer to the leading edge (686 ± 801 Pa), while the leading edge and the tail (494 ± 537 Pa) regions were mechanically indistinguishable. The frequency-dependent elastic and viscous moduli also display a similar regional difference. Over the studied frequency range (100 to 300 Hz), the complex viscoelastic moduli display the partial rubber plateau behavior where the elastic moduli are greater than the viscous moduli for a given frequency. The non-disparaging viscous modulus indicates that the neutrophils display a viscoelastic dynamic behavior rather than a perfect elastic behavior like polymer gels. In addition, we found no regional difference in the structural damping coefficient between the leading edge and the cell body. Thus, we conclude that despite the lower loss and storage moduli, the leading edges of the human neutrophils display partially elastic properties similar to the cell body. These results suggest that the lower elastic moduli in the leading edges are more favorable for the elastic fluctuation of actin filaments, which supports the polymerization of the actin filaments leading to the active protrusion during the immune response.     

Models of skeletal muscle to explain the increase in passive stiffness in desmin knockout muscle

Absence of desmin in skeletal muscle was found to induce an increase in passive stiffness. The present study aimed at developing rheological models of passive muscle to explain this stiffening. Models were elaborated by using experimental data depicting muscle viscoelastic behaviour. The experimental protocol included stepwise extension tests applied on control and desmin knockout soleus muscles from mice. Linear and non-linear models were composed of elastic and viscous elements. They were constructed with the aim at taking the presence or absence of desmin into account by simulating desmin as an elastic element. Furthermore, associated adaptation of connective tissues in absence of desmin was modelled as an additional elastic element. Differences in passive behaviour induced by absence of desmin were predicted by using a linear model and a non-linear one. The non-linear model was selected because: (1) it is able to predict experimental viscoelastic kinetics accounting for the increase in passive stiffness in muscles lacking desmin, (2) its design is consistent with morphological data, and (3) stiffness characteristics of its elements are in accordance with the literature. Finally, this modelling approach demonstrates that both absence of desmin and adaptation of connective tissue are required to explain the increase in passive stiffness in desmin knockout muscles.     

Shear viscoelasticity of suspensions of biological cells with viscoelastic membrane II

To discuss the relaxation phenomena of biological cell suspensions, we calculate the complex intrinsic viscosity of dispersions of spherical cells with viscoelastic membrane as a function of the frequency taking account of interfacial tension at both the interfaces of the membrane. The Maxwell model and two kinds of the three-parameter models are used to describe the viscoelasticity of the cell membrane. The results are computed mainly for the Maxwell model similarly in case of the Voigt Model (Abe, K., Takano, Y. and Sakanishi, A. Biorheology 21 405-414, 1984). The computed results of the four models, the Voigt, the Maxwell and the two kinds of the three-parameter viscoelastic models, are compared with the experimental data.     

Red cell extensional recovery and the determination of membrane viscosity.

A theory of membrane viscoelasticity developed by Evans and Hochmuth in 1976 is used to analyze the time-dependent recovery of an elongated cell. Before release, the elongated cell is the static equilibrium where external forces are balanced by membrane elastic force resultants. Upon release, the cell recovers its initial shape with a time-dependent exponential behavior characteristic of the viscoelastic solid model. It is shown that the model describes the time-dependent recovery process very well for a time constant in the range of 0.1-0.13 s. The time constant is the ratio membrane surface viscosity eta:membrane surface elasticity mu. Measurements for the shear modulus mu of 0.006 dyne/cm give a value for the surface viscosity of red cell membrane as a viscoelastic solid material of eta = mu tc = (6-8) X 10(-4) poise . cm.     

An orthotropic viscoelastic material model for passive myocardium: theory and algorithmic treatment

This contribution presents a novel constitutive model in order to simulate an orthotropic rate-dependent behaviour of the passive myocardium at finite strains. The motivation for the consideration of orthotropic viscous effects in a constitutive level lies in the disagreement between theoretical predictions and experimentally observed results. In view of experimental observations, the material is deemed as nearly incompressible, hyperelastic, orthotropic and viscous. The viscoelastic response is formulated by means of a rheological model consisting of a spring coupled with a Maxwell element in parallel. In this context, the isochoric free energy function is decomposed into elastic equilibrium and viscous non-equilibrium parts. The baseline elastic response is modelled by the orthotropic model of Holzapfel and Ogden [Holzapfel GA, Ogden RW. 2009. Constitutive modelling of passive myocardium: a structurally based framework for material characterization. Philos Trans Roy Soc A Math Phys Eng Sci. 367:3445–3475]. The essential aspect of the proposed model is the account of distinct relaxation mechanisms for each orientation direction. To this end, the non-equilibrium response of the free energy function is constructed in the logarithmic strain space and additively decomposed into three anisotropic parts, denoting fibre, sheet and normal directions each accompanied by a distinct dissipation potential governing the evolution of viscous strains associated with each orientation direction. The evolution equations governing the viscous flow have an energy-activated nonlinear form. The energy storage in the Maxwell branches has a quadratic form leading to a linear stress–strain response in the logarithmic strain space. On the numerical side, the algorithmic aspects suitable for the implicit finite element method are discussed in a Lagrangian setting. The model shows excellent agreement compared to experimental data obtained from the literature. Furthermore, the finite element simulations of a heart cycle carried out with the proposed model show significant deviations in the strain field relative to the elastic solution.     

Cell motility and local viscoelasticity of fibroblasts

Viscoelastic changes of the lamellipodial actin cytoskeleton are a fundamental element of cell motility. Thus, the correlation between the local viscoelastic properties of the lamellipodium (including the transitional region to the cell body) and the speed of lamellipodial extension is studied for normal and malignantly transformed fibroblasts. Using our atomic force microscopy-based microrheology technique, we found different mechanical properties between the lamellipodia of malignantly transformed fibroblasts (H-ras transformed and SV-T2 fibroblasts) and normal fibroblasts (BALB 3T3 fibroblasts). The average elastic constants, K, in the leading edge of SV-T2 fibroblasts (0.48 +/- 0.51 kPa) and of H-ras transformed fibroblasts (0.42 +/- 0.35 kPa) are significantly lower than that of BALB 3T3 fibroblasts (1.01 +/- 0.40 kPa). The analysis of time-lapse phase contrast images shows that the decrease in the elastic constant, K, for malignantly transformed fibroblasts is correlated with the enhanced motility of the lamellipodium. The measured mean speeds are 6.1 +/- 4.5 microm/h for BALB 3T3 fibroblasts, 13.1 +/- 5.2 microm/h for SV-T2 fibroblasts, and 26.2 +/- 11.5 microm/h for H-ras fibroblasts. Furthermore, the elastic constant, K, increases toward the cell body in many instances which coincide with an increase in actin filament density toward the cell body. The correlation between the enhanced motility and the decrease in viscoelastic moduli supports the Elastic Brownian Ratchet model for driving lamellipodia extension.     

Passive deformation analysis of human leukocytes

The following analysis presents an experimental and theoretical study of the passive viscoelastic behavior of human leukocytes. Individual neutrophils in EDTA were observed both during their partial aspiration into a small micropipette and after expulsion from a large micropipette where the cell had been totally aspirated and deformed into a sausage shape. To analyze the data, a passive model of leukocyte rheology has been developed consisting of a cortical shell containing a Maxwell fluid which describes the average properties of the cell cytoplasm. The cortical shell represents a crosslinked actin layer near the surface of the cell and is assumed to be under pre-stressed tension. This model can reproduce the results of experiments using micropipette for both short-time small deformation and slow recovery data after large deformation. In addition, a finite element scheme has been established for the same model which shows close agreement with the analytical solution.     

Neutrophil transit times through pulmonary capillaries: the effects of capillary geometry and fMLP-stimulation

The deformations of neutrophils as they pass through the pulmonary microcirculation affect their transit time, their tendency to contact and interact with the endothelial surface, and potentially their degree of activation. Here we model the cell as a viscoelastic Maxwell material bounded by constant surface tension and simulate indentation experiments to quantify the effects of (N-formyl-L-methionyl-L-leucyl-L-phenylalanine (fMLP)-stimulation on its mechanical properties (elastic shear modulus and viscosity). We then simulate neutrophil transit through individual pulmonary capillary segments to determine the relative effects of capillary geometry and fMLP-stimulation on transit time. Indentation results indicate that neutrophil viscosity and shear modulus increase by factors of 3.4, for 10(-9) M fMLP, and 7.3, for 10(-6) M fMLP, over nonstimulated cell values, determined to be 30.8 Pa.s and 185 Pa, respectively. Capillary flow results indicate that capillary entrance radius of curvature has a significant effect on cell transit time, in addition to minimum capillary radius and neutrophil stimulation level. The relative effects of capillary geometry and fMLP on neutrophil transit time are presented as a simple dimensionless expression and their physiological significance is discussed.     

Mechanical behavior of human embryonic stem cell pellet under unconfined compression

As a prelude to the understanding of mechanotransduction in human embryonic stem cell (hESC) differentiation, the mechanical behavior of hESCs in the form of cell pellet is studied. The pellets were tested after 3 or 5 weeks of cell culture in order to demonstrate the effect of the duration of cell culture on the mechanical properties of the pellets. A micromechanical tester was used to conduct unconfined compression on hESC pellet, and experimental, numerical, and analytical methods were combined to determine the mechanical properties of hESC pellet. It is assumed that the mechanical behavior of hESC pellets can be described by an isotropic, linear viscoelastic model consisting of a spring and two Maxwell units in parallel, and the Poisson’s ratio of the hESC pellet is constant based on pellet deformation in the direction perpendicular to the compression direction. Finite element method (FEM) simulation was adopted to determine the values of Poisson’s ratio and the five parameters contained in the viscoelastic model. The variations of Poisson’s ratio and the initial elastic modulus are found to be larger compared with those of the four other parameters. Results show that longer duration of cell culture leads to higher modulus of hESC pellet. The effect of pellet size error on the values of mechanical parameters determined is studied using FEM simulation, and it is found that the effect of size error on Poisson’s ratio and initial elastic modulus is much larger than that on the other parameters.