Finite Element Formulation of Multiphasic Shell Elements for Cell Mechanics Analyses in FEBio

Abstract

With the recent implementation of multiphasic materials in the open-source finite element (FE) software FEBio, three-dimensional (3D) models of cells embedded within the tissue may now be analyzed, accounting for porous solid matrix deformation, transport of interstitial fluid and solutes, membrane potential, and reactions. The cell membrane is a critical component in cell models, which selectively regulates the transport of fluid and solutes in the presence of large concentration and electric potential gradients, while also facilitating the transport of various proteins. The cell membrane is much thinner than the cell; therefore, in an FE environment, shell elements formulated as two-dimensional (2D) surfaces in 3D space would be preferred for modeling the cell membrane, for the convenience of mesh generation from image-based data, especially for convoluted membranes. However, multiphasic shell elements are yet to be developed in the FE literature and commercial FE software. This study presents a novel formulation of multiphasic shell elements and its implementation in FEBio. The shell model includes front- and back-face nodal degrees-of-freedom for the solid displacement, effective fluid pressure and effective solute concentrations, and a linear interpolation of these variables across the shell thickness. This formulation was verified against classical models of cell physiology and validated against reported experimental measurements in chondrocytes. This implementation of passive transport of fluid and solutes across multiphasic membranes makes it possible to model the biomechanics of isolated cells or cells embedded in their extracellular matrix (ECM), accounting for solvent and solute transport.