Abstract
Extrusion-based bioprinting is an enabling biofabrication technique that is used to create heterogeneous tissue constructs according to patient-specific geometries and compositions. The optimization of bioinks as per requirements for specific tissue applications is an essential exercise in ensuring clinical translation of the bioprinting technologies. Most notably, optimum hydrogel polymer concentrations are required to ensure adequate mechanical properties of bioprinted constructs without causing significant shear stresses on cells. However, experimental iterations are often tedious for optimizing the bioink properties. In this work, a nonlinear finite element modeling approach has been undertaken to determine the effect of different bioink parameters such as composition, concentration on the range of stresses being experienced by the cells in the bioprinted construct. The stress distribution of the cells at different parts of the constructs has also been modeled. It is found that both bioink chemical compositions and concentrations can substantially alter the stress effects experienced by the cells. Concentrated regions of softer cells near pore regions were found to increase stress concentrations by almost three times compared with stress generated in cells away from the pores. The study provides a method for rapid optimization of bioinks, design of bioprinted constructs, as well as toolpath plans for fabricating constructs with homogenous properties.