A 2-D finite element model for the intervertebral disc in which quadriphasic theory is coupled to the transport of solutes involved in cellular nutrition was developed for investigating the main factors contributing to disc degeneration. Degeneration is generally considered to result from chronic disc cell nutrition insufficiency, which prevents the cells from renewing the extracellular matrix and thus leads to the loss of proteoglycans. Hence, the osmotic power of the disc is decreased, causing osmomechanical impairments. Cellular metabolism depends strongly on the oxygen, lactate and glucose concentrations and on pH in the disc. To study the diffusion of these solutes in a mechanically or osmotically loaded disc, the osmomechanical and diffusive effects have to be coupled. The intervertebral disc is modeled here using a plane strain formulation at the equilibrium state under physiological conditions after a long rest period (called unloaded state). The correlations between solute distribution and various properties of healthy and degenerated discs are investigated. The numerical simulation shows that solute distribution in the disc depends very little on the elastic modulus or the proteoglycan concentration but greatly on the porosity, diffusion coefficient and endplate diffusion area. This coupled model therefore opens new perspectives for investigating intervertebral disc degeneration mechanisms.

In an attempt to understand the role of structural rearrangement onto the cell response during imposed cyclic stresses, we simulated numerically the frequency-dependent behavior of a viscoelastic tensegrity structure (VTS model) made of 24 elastic cables and 6 rigid bars. The VTS computational model was based on the nonsmooth contact dynamics (NSCD) method in which the constitutive elements of the tensegrity structure are considered as a set of material points that mutually interact. Low amplitude oscillatory loading conditions were applied and the frequency response of the overall structure was studied in terms of frequency dependence of mechanical properties. The latter were normalized by the homogeneous properties of constitutive elements in order to capture the essential feature of spatial rearrangement. The results reveal a specific frequency-dependent contribution of elastic and viscous effects which is responsible for significant changes in the VTS model dynamical properties. The mechanism behind is related to the variable contribution of spatial rearrangement of VTS elements which is decreased from low to high frequency as dominant effects are transferred from mainly elastic to mainly viscous. More precisely, the elasticity modulus increases with frequency while the viscosity modulus decreases, each evolution corresponding to a specific power-law dependency. The satisfactorily agreement found between present numerical results and the literature data issued from in vitro cell experiments suggests that the frequency-dependent mechanism of spatial rearrangement presently described could play a significant and predictable role during oscillatory cell dynamics.

Cell deformability and mechanical responses of living cells depend closely on the dynamic changes in the structural architecture of the cytoskeleton (CSK). To describe the dynamic reorganization and the heterogeneity of the prestressed multi-modular CSK, we developed a two-dimensional model for the CSK which was taken to be a system of tension and compression interactions between the nodes in a divided medium. The model gives the dynamic reorganization of the CSK consisting of fast changes in connectivity between nodes during medium deformation and the resulting mechanical behavior is consistent with the strain-hardening and prestress-induced stiffening observed in cells in vitro. In addition, the interaction force networks which occur and balance to each other in the model can serve to identify the main CSK substructures: cortex, stress fibers, intermediate filaments, microfilaments, microtubules and focal adhesions. Removing any of these substructures results in a loss of integrity in the model and a decrease in the prestress and stiffness, and suggests that the CSK substructures are highly interdependent. The present model may therefore provide a useful tool for understanding the cellular processes involving CSK reorganization, such as mechanotransduction, migration and adhesion processes.

A novel topology optimization model based on homogenization methods was developed for predicting bone density distribution and anisotropy, assuming the bone structure to be a self-optimizing biological material which maximizes its own structural stiffness. The feasibility and efficiency of this method were tested on a 2D model for a proximal femur under single and multiple loading conditions. The main aim was to compute homogenized optimal designs using an optimal laminated microstructure. The computational results showed that high bone density levels are distributed along the diaphysis and form arching struts within the femoral head. The pattern of bone density distribution and the anisotropic bone behavior predicted by the model in the multiple load case were both in good agreement with the structural architecture and bone density distribution occurring in natural femora. This approach provides a novel means of understanding the remodeling processes involved in fracture repair and the treatment of bone diseases.