Three-dimensional finite element analysis of residual stress in arteries

Calculation of residual stress in arteries, using the analytical approach has been quite valuable in our understanding of its critical role in vascular mechanics. Stresses are calculated at the central section of an infinitely long tube by imposing a constant axial stretch while deforming the artery from the stress-free state to its unloaded state. However, segments used to perform opening-angle measurements have finite lengths. Further, the stress-free artery configuration is assumed to be circular. Experiments show that they are slightly noncircular. The numerical approach to residual stress calculation can allow us to study both these issues. Using 3D cylindrical geometries and an isotropic material model, we investigated how segment length can affect residual stress calculations and identified the appropriate segment length for experiments. Further, we recorded and used the true noncircular stress-free state of an artery segment, computed the residual stress distribution, and compared it to that from a similar, but circular segment. Our findings suggest that segment length must be ten times the wall thickness for it to be long enough. We also found that the circularity assumption may be a reasonable approximation for typical arteries.