Vascular Stent Deployment FSI
Overview
先生、血管ステントの展開シミュレーションでFSIが必要な理由は何ですか？
Theory and Physics
A stent is a medical device that expands a metal mesh structure inside a blood vessel to widen a narrowed vessel. The blood flow pattern after stent deployment strongly depends on the stent shape and vessel wall deformation, and is directly linked to the risk of restenosis (neointimal hyperplasia). It is known that restenosis is more likely to occur in areas with low WSS (Wall Shear Stress), so it's important to accurately predict the blood flow field and wall stress using FSI.
What types of stents are there?
There are three main types. Balloon-expandable stents (BES: coronary artery stents), self-expanding stents (SES: Nitinol alloy, for carotid/peripheral use), and drug-eluting stents (DES: coated with restenosis-inhibiting drugs). Their mechanical behaviors differ, so the structural models also change.
Governing Equations
What system of equations is solved for stent deployment FSI?
It's a coupled problem of three structures (stent, balloon, vessel wall) and blood flow.
The stent is modeled with SUS316L or Nitinol (Hyperelasticity). The Auricchio-Taylor model is standard for Nitinol's hyperelasticity.
The transformation strain $\boldsymbol{\varepsilon}^{tr}$ occurs during austenite-martensite phase transformation.
The vessel wall uses the Holzapfel-Gasser-Ogden model, and blood flow uses the incompressible Navier-Stokes equations. If plaque is present, a linear elastic or plastic model is used to represent the hard layer.
How is the balloon modeled?
The balloon is a polymer thin film (e.g., polyamide) modeled with membrane elements. Internal pressure is increased stepwise to expand the stent from its crimped diameter to the target diameter. Sometimes the folded shape of the balloon is also reproduced, but it's a trade-off with computational cost.
Nitinol—How "Shape Memory" Sparked a Stent Revolution
Nitinol (Ni-Ti alloy), widely used in modern vascular stents, possesses two characteristics: shape memory effect and superelasticity. It transforms to the austenite phase near body temperature (37°C), allowing a stent inserted in a contracted state to expand naturally inside the vessel. To handle Nitinol in FSI theory, a hyperelastic constitutive law (such as the Brinson rule) is required, not the usual linear elastic model. Implementing this material model is one of the biggest factors raising the difficulty level of stent FSI analysis.
Physical Meaning of Each Term
- Structural-Thermal Coupling Term: Thermal expansion due to temperature change induces structural deformation, and deformation affects the temperature field. $\sigma = D(\varepsilon - \alpha \Delta T)$. 【Everyday Example】Railroad tracks expanding in summer, narrowing the gaps—a classic case of temperature rise → Thermal Expansion → stress generation. Warping of electronic circuit boards after soldering is also due to differences in thermal expansion coefficients between materials. Thermal stress occurs in engine cylinder blocks due to temperature differences between hot and cold parts, potentially leading to cracks.
- Fluid-Structure Interaction (FSI) Term: Bidirectional interaction where fluid pressure/shear force deforms the structure, and structural deformation changes the fluid domain. 【Everyday Example】Suspension bridge cables vibrating in strong wind (Vortex-Induced Vibration)—wind force shakes the structure, the shaken structure alters the wind flow, further amplifying vibration. Blood flow in the heart and elastic deformation of vessel walls, and aircraft wing flutter (aeroelastic instability) are also typical FSI problems. One-way coupling may suffice in some cases, but bidirectional coupling is essential for large deformations.
- Electromagnetic-Thermal Coupling Term: A feedback loop where Joule heating $Q = J^2/\sigma$ causes temperature rise, and temperature change alters electrical resistance. 【Everyday Example】Nichrome wire in an electric stove heats up (Joule heat) and glows red when current flows—temperature rise changes resistance, altering current distribution. Eddy current heating in IH cooking heaters and increased sag of power lines due to temperature rise are also examples of this coupling.
- Data Transfer Term: Interpolation resolves mesh mismatch between different physical fields. 【Everyday Example】When calculating "feels-like temperature" by combining "temperature data" and "wind data" in weather forecasting, interpolation is needed if observation points differ—similarly in CAE coupled analysis, structural and CFD meshes generally don't match, so the accuracy of data transfer (Interpolation) at the interface directly affects result reliability.
Assumptions and Applicability Limits
- Weak coupling assumption (one-way coupling): Effective when one physical field affects the other but the reverse is negligible.
- Cases requiring strong coupling: Large deformations in FSI, cases with strong temperature dependence in electromagnetic-thermal coupling.
- Separation of time scales: When characteristic times of each physical field differ significantly, subcycling can improve efficiency.
- Interface condition consistency: Ensure energy/momentum conservation at the coupling interface is satisfied numerically.
- Non-applicable cases: When three or more physical fields are strongly coupled simultaneously, monolithic methods may be necessary.
Dimensional Analysis and Unit Systems
| Variable | SI Unit | Notes / Conversion Memo |
|---|---|---|
| Thermal expansion coefficient $\alpha$ | 1/K | Steel: ~12×10⁻⁶, Aluminum: ~23×10⁻⁶ |
| Coupled interface force | N/m² (pressure) or N (concentrated force) | Verify force balance on fluid and structure sides. |
| Data transfer error | Dimensionless (%) | Interpolation accuracy depends on mesh density ratio. Below 5% is a guideline. |
Numerical Methods and Implementation
It's usually divided into two stages.
Stage 1: Stent Deployment (Structure Only)
1. Stent crimping (diameter reduction)
2. Expansion via balloon pressurization
3. Recoil (Springback) after balloon deflation
Stage 2: Post-Deployment FSI
1. Set the deformed shape after deployment as the initial condition
2. Import residual stresses
3. Execute FSI analysis for pulsatile blood flow
Importing residual stresses is important, I see.
Significant residual stresses exist in the stent after deployment; ignoring them drastically reduces the prediction accuracy of wall stress. Use Abaqus's *IMPORT function or Ansys's ICTRL/RESUME command to carry over the stress field from the previous analysis.
Contact Algorithm
How is contact between the stent and vessel wall handled?
There are three contact pairs: stent-balloon, stent-vessel wall, balloon-vessel wall.
How is contact between the stent and vessel wall handled?
There are three contact pairs: stent-balloon, stent-vessel wall, balloon-vessel wall.
| Contact Pair | Method | Friction Coefficient |
|---|---|---|
| Stent-Balloon | Surface-to-surface contact | 0.05–0.1 |
| Stent-Vessel Wall | Surface-to-surface contact | 0.1–0.2 |
| Stent-Plaque | Surface-to-surface contact | 0.2–0.3 |
In Abaqus, *CONTACT PAIR or General Contact is used; in LS-DYNA, *CONTACT_AUTOMATIC_SURFACE_TO_SURFACE is used. The Augmented Lagrangian method can suppress penetration better than the penalty method.
Mesh Strategy
Meshing the thin stent struts seems challenging.
Divide the strut cross-section (typically 80–120μm × 80–120μm) with hexahedral elements, at least 3×3. The entire stent becomes 500k–2M elements.
Locally refine the vessel wall near stent struts, using element sizes below 1/3 of the strut width. The fluid mesh needs elements placed even in the gaps between struts, often resulting in 1M–5M elements total.
With that scale, how long is the computation time?
Deployment analysis (quasi-static) typically takes 8–24 hours; FSI (several heartbeats) takes 24–72 hours. This is an estimate assuming implicit method parallel computation with 16–32 cores.
Two-Step Analysis of Stent Deployment—Solving Folding and Deployment Separately
FSI analysis of vascular stent deployment is typically performed in two steps: "structural analysis of the folding process" → "FSI analysis of deployment inside the vessel." The key point is carrying over the residual stresses generated during folding as initial conditions for the deployment analysis. If this is omitted and analysis starts from the deployed stent shape, the actual deployment force and load on the vessel can be underestimated by 20–30%. This two-step method is also increasingly required in medical device regulatory submission documents.
Monolithic Method
Solves all physical fields simultaneously as one system of equations. Stable for strong coupling but complex to implement and memory-intensive.
Partitioned Method (Partitioned Iterative Method)
Solves each physical field independently, exchanging data at the interface. Easy to implement and can utilize existing solvers. Suitable for weak coupling.
Interface Data Transfer
Nearest neighbor (simplest but low accuracy), projection (conservative), RBF interpolation (robust for mesh mismatch). Balance between conservation and accuracy is important.
Sub-iteration
Perform sufficient iterations within each coupling step to ensure interface condition consistency. Residual criteria are scaled based on typical values of each physical field.
Aitken Relaxation
Automatically adjusts the relaxation factor for coupling iterations. An adaptive method that prevents divergence from over-relaxation and accelerates convergence.
Stability Condition
Beware of the added mass effect (in fluid-structure coupling when structural density ≈ fluid density). If unstable, apply Robin-type interface conditions or the IQN-ILS method.
Analogy for Aitken Relaxation
Aitken relaxation is like "balancing a seesaw." If one side pushes too hard, the other side flies up, and the recoil causes it to push too hard again—Aitken relaxation automatically adjusts the pushing force to suppress this oscillation. When coupling iterations oscillate and fail to converge, it's an adaptive method that automatically adjusts the next correction amount based on the previous correction.
Practical Guide
1. Stent Shape Creation: Create the strut pattern from CAD in the deployed shape, then wrap it into a cylinder.
2. Vessel Shape Acquisition: Obtain patient-specific shape from CT/MRI, or use an IDEAL shape (straight tube + stenosis).
3. Material Definition: Stent (SUS316L: Elasto-Plastic / Nitinol: Auricchio hyperelasticity), Vessel wall (Holzapfel-Gasser-Ogden), Plaque (Elasto-Plastic)
4. Crimp ~ Deployment Analysis: Abaqus/Explicit or Abaqus/Standard
5. FSI Analysis of Deployed Shape: Ansys System Coupling or Abaqus co-simulation
6. Postprocessing: Stent von Mises stress (fatigue evaluation), WSS, OSI (restenosis risk evaluation)
Do you also do fatigue evaluation?
The FDA (U.S. Food and Drug Administration) mandates fatigue life evaluation for stents. Fatigue durability for 400 million cycles (equivalent to 10 years) is required. Plot mean stress and alternating stress on a Goodman diagram and confirm they are below the fatigue limit.
A safety factor SF = 2 or more is a common standard.
Common Pitfalls
Where are the tricky points in practice?
| Problem | Cause | Countermeasure |
|---|---|---|
| Stent not adhering to vessel wall | Overestimation of recoil | Set balloon overexpansion ratio (10–20%) |
| Strut fracture judgment | Mesh-dependent stress concentration | Use submodel method for local refinement |
| Unrealistic dog-boning | Balloon pressurization sequence | Stepwise division with center-first pressurization |
| WSS distribution asymmetric left-right | Influence of vessel curvature, bifurcation | Sufficient inlet development length |
What is dog-boning?
A phenomenon in balloon-expandable stents where the ends expand first and the center expands later. It causes excessive stress on the vessel wall at the stent ends, leading to restenosis. It is controlled by balloon compliance and stent cell design.
Related Topics
なった
詳しく
報告