Simulating Stent Expansion Using Finite Element CAD

Simulating cardiovascular stent expansion requires integrating Computer-Aided Design (CAD) and Finite Element Analysis (FEA) to predict how a metallic or polymeric scaffold deforms inside an artery. This article explains the end-to-end computational process: constructing micro-scale CAD geometries, assigning non-linear material models, configuring contact mechanics between the balloon, stent, and arterial wall, and utilizing explicit FEA solvers to analyze deployment risks such as strut fracture, tissue trauma, and radial recoil.

1. Geometric CAD Modeling

The simulation begins by generating an accurate, parametric 3D CAD model of the stent in its unexpanded or crimped state. Given the repeating cellular architecture of stents, designers typically model a single unit cell and pattern it circumferentially and longitudinally. Stents feature sub-millimeter strut thicknesses, often ranging from 60 to 120 microns, requiring precise spline definition and rounded edges to avoid artificial stress singularities. The CAD assembly also incorporates auxiliary bodies: a cylindrical or folded angioplasty balloon and a multilayered, curved arterial wall segment with optional atherosclerotic plaque geometry.

2. Defining Non-Linear Material Models

Stent expansion involves severe plastic deformation or phase transformation, necessitating advanced material constitutive laws:

  • Balloon-Expandable Stents (Stainless Steel 316L, Cobalt-Chromium): Modeled using non-linear elastoplasticity with isotropic or kinematic hardening. The material yields under internal balloon pressure, holding its shape via plastic deformation after balloon deflation.
  • Self-Expanding Stents (Nitinol): Modeled using superelastic constitutive formulations (such as the Auricchio model). This accounts for the stress-induced martensitic phase transformation during crimping and the subsequent phase reversal to austenite during unconstrained expansion at body temperature.
  • Arterial Tissue and Plaque: Modeled as anisotropic, hyperelastic materials (e.g., using Holzapfel-Gasser-Ogden or Mooney-Rivlin models) to accurately mimic the stiffness, collagen fiber orientation, and compliance of the intima, media, and adventitia layers.

3. Meshing Strategy

Because the struts undergo high bending strains, mesh selection is critical to avoid numerical issues like shear locking. High-quality hexahedral (brick) elements with reduced integration are preferred for structured strut regions, typically requiring at least three to four elements through the strut thickness to accurately resolve bending moments. The arterial tissue and plaque are meshed with refined tetrahedral or hexahedral elements, with local refinement applied to inner lumen surfaces where contact occurs.

4. Boundary Conditions and Contact Mechanics

The expansion phase introduces complex multi-body interactions:

  • Contact Formulations: Penalty-based or augmented Lagrange contact algorithms are implemented to govern interactions across three distinct interfaces: balloon-to-stent, stent-to-vessel, and strut-to-strut (self-contact during crimping). Friction coefficients are applied to prevent artificial slipping.
  • Loading Steps: The simulation is divided into chronological phases: crimping (radial displacement boundary applied to external rigid plates), tracking/insertion, deployment (applying radial pressure to the balloon interior or retracting a delivery sheath for Nitinol), and deflation/recoil (removing internal pressure to evaluate vessel pushback).

5. Solving via Non-Linear FEA

Due to extreme geometric non-linearities, large strains, and abrupt contact changes, simulations typically utilize non-linear FEA solvers such as Abaqus/Explicit or ANSYS LS-DYNA. While the actual physical deployment occurs over several seconds, explicit dynamic simulations are often mass-scaled or executed over compressed timeframes to reduce computational cost. Engineers verify quasi-static equilibrium by ensuring the kinetic energy of the system remains negligible (typically below 5% to 10%) compared to the total internal strain energy.

6. Post-Processing and Performance Evaluation

The solved model yields quantitative biomechanical metrics used to optimize stent design and ensure regulatory compliance:

  • Dogboning: Measuring whether the proximal and distal stent ends expand before the central section, which can cause vessel dissection.
  • Radial Recoil: Calculating the percentage of diameter loss after balloon deflation due to arterial elasticity.
  • Foreshortening: Tracking the reduction in stent axial length as diameter increases.
  • Stress and Strain Concentrations: Identifying regions exceeding the ultimate tensile strength of the metal or the threshold for arterial rupture.
  • Fatigue Resistance: Exporting peak stress-strain states across cardiac pulsation cycles into Goodman or Soderberg fatigue frameworks to predict structural durability over millions of cycles.