Development of a transient large strain contact method for biological heart valve simulations

A new 2D method to implement transient contact using Comsol Multiphysics (finite element analysis software that enables multiphysics simulations) is described, which is based on Hertzian contact. This method is compared to the existing (default) contact method that does not enable real transient simulations, but instead performs steady-state solutions where time is a constant. The two types of contact modelling have been applied to simple 2D biological heart valve models, undergoing strains in the region of 10% under 20 kPa pressure (applied over 0.3 s). Both the methods predicted comparable stress patterns, locations of peak stresses, deformations and directions of principal stress. The default contact method predicted slightly greater contact stresses, but spreads over a shorter surface length than the new contact method. The default contact method is useful for contact systems with little transient dependency, due to ease of use. However, where transient conditions are important the new contact method is preferred.