Study of Non-Linear Deformation of Peripheral Stent Mechanics Using Computational Approach
AFFILIATIONS
AFFILIATIONS
Corresponding author (Address):
AFFILIATIONS
AFFILIATIONS
AFFILIATIONS
AFFILIATIONS
doi:10.17303/jber.2016.102
Citation: Hua Yang (2016) Study of Non-Linear Deformation of Peripheral Stent Mechanics Using Computational Approach. J Biomed Eng 1: 1-10.
Stent implants, specifically those implanted in Superficial Femoral Artery (SFA) in the lower limbs, are prone to high failure rates. The failures can occur due to abnormal behavior of the blood flow, mechanical forces from the arterial wall imposed on the stent, as well as leg movements such as bending, torsion, compression, and elongation that can create dynamic arterial environment imposed on the stent. In order to prevent high failure rates of peripheral stents there is a need to understand the unique biomechanical environment and the characteristics of the SFA arterial segment where stent is implanted. This study presents several steps that need to be taken into consideration before the complex and combined stent-artery-environment model can be created and analyzed. Results show finite element computational approach that demonstrates stent characteristics upon insertion inside the artery, creation of arterial segment with geometry and characteristics from real patient data, and theoretical approach to accurately impose blood flow onto the artery-stent model. This study recommends procedure for analyzing biomechanical environment of the stent-artery system using real artery geometry. Results show that stent design and artery specific geometry play a critical role in stent performance.
Keywords: Peripheral stent; Mechanical behavior; Non-linear deformation; Finite elements modeling; Stent-artery system; Restenosis
Cardiovascular disease is often related to atherosclerosis which can cause narrowing, rupture or erosion of the arterial wall, with eventual reduction (i.e. > 50% blockage) or occlusion (i.e. 100% blockage) of blood flow [1]. Among several different treatments is stent implant procedure, which is the insertion of a metal and/or plastic tubular structure that has the ability to expand into a cylindrical shape, either by use of balloon or a self-expansion mechanism [1,2]. It is reported that about 700,000 stent procedures (in coronary and peripheral arteries together) are done in the U.S. annually with about 50% stent failure rate in the peripheral arteries reported [3]. High failure rate of stent implants is a result of complex dynamic environment when delivered to a region of complex peripheral artery anatomy. Finite element analysis and numerical simulation tests are increasingly important for accurate design and manufacturing of stents in order to meet market demands and improve patient care. The possibilities of Finite Element Analysis (FEA) and Computational Fluid Dynamics (CFD), combined with contemporary imaging techniques, greatly facilitate stent research. Tremendous amount of research has been done on coronary artery stenting, on stent geometry, and on estimating artery-stent interaction however, very little is reported on superficial artery stenting and the process of using real artery geometry from patient data for the computational models [4-7].
Most studies estimate the geometry of the artery as cylinder with uniform diameter across its length where in reality arteries have curves and tortuosity specific to the patient anatomy [1,2,4]. Several aspects need to be considered in order for computational models to accurately estimate the artery-stent biomechanical environment and the issues causing stent failure. Stent geometry and loading conditions need to be accurately defined using principles of solid mechanics [8]. The real artery geometry is typically obtained by using Computed Tomography (CT) images, or intravascular ultrasound (IVUS) images, and reconstructing them using the Mimics Innovation Suite Software that will give 3D artery model with real geometry [9]. The constructed 3D model is the transferred to finite elements modeling software package coupled with stent model and various analyses are performed. To complete the numerical analysis especially the nonlinear analysis, the finite element method requires a lot of information, including the geometry and the material properties of the stents, the effect of the catheter and the artery, and the accurate application of boundary and loading conditions. MSC commercial software allows for proper analysis of these inherently difficult (e.g. non-homogeneous, non-linear) 3D models [10]. In our previous study [2] we have used MSC Software to validate the linear deformation of specific stent designs also analyzed theoretically and here we use similar computational approach to simulate nonlinear deformation behavior of the stent. In this study the appropriate finite element models are developed in Patran and solved with NASTRAN [11]. Prevalent stents have numerous design patterns, the widespread interest focused on the properties of stent designs in terms of multiple types geometry [1,5,6]. The elevated stress region largely depends on the strut shape of stents in terms of mechanical stress corresponding to the clinical observations [12,13]. The results from this study show that the strut length, separation angle and cell asymmetry are crucial to the performance of stent design. Finite element modeling approach is used to evaluate the deformation in balloon expanded stent under internal pressure first. Results from this study show the nonlinear behavior of the stent, and the plastic region with the maximum equivalent stress noted to occur at the 4.5 millimeter (mm) expansion diameter. Approach for generating artery 3D models from real patient data and imposing dynamic blood flow behavior is also presented. A fundamental understanding and step by step process for generating stent-artery-environment models have limited reports in the literature published to date.
Stent expansion is analyzed assuming struts within the stent expand at a constant rate via a plastic hinge mechanism. For each single hinge-strut unit, the beam theory (explained in prior study [2]) continuously contributes to the understanding of the complex and nonlinear inflation of expandable stent. The finite element model estimates and reproduces the characteristics during crimp and inflation mechanism. In the stent model, the stent geometry consists of several rows of identical parts, one intermediate row at each end (the same geometry, but different strut width), and one final row at each end that has the same basic geometry (and increased strut width) with the exception that there is no ‘connector’ on the outboard ends (as there are no additional rows to be attached to).
The repetitive nature of the geometry allows for a theoretical model or a finite element model of a portion of the entire stent to accurately characterize the stress state within the entire stent when inflated. This approach allows the finite element mesh to be finer in the region that is modeled, as opposed to using a similar number of elements for an entire stent. For example, the assumed strain formulation allows the element to accurately capture bending characteristics with as little as one element through the thickness. Similarly, in this study is assumed that due to symmetry unit structure analysis will reflect the deformation of the entire stent, where expansion, deflection and stress/strain will be also predicted. 3mm diameter stent is inflated to 4.5mm which produces large mechanical loads with nonlinear characteristics and these are considered in our model. In real application, initially the stent is an empty and collapsed balloon loaded on a guide wire, known as a balloon catheter, the combo is passed into the narrowed arteries and then inflated to a 4.5mm using pressures approximately 75 to 500 times normal blood pressure [14]. Additionally, most healthy arteries yield an average hoop stress of approximately 100 KPa which gives an idea what is the stent undergoing once placed inside the artery. By simulating the inflated stent structure, the maximum stress will be evaluated as a result of stent inflation to diameter of up to 4.5 mm. A Von Mises yield condition, commonly used with metallic alloys, along with a kinematic hardening law are used in the analysis. When the material is strained beyond the yield point, the yield surface must change for continued strain hardening. The kinematic hardening law is used for the analysis herein, to take into account the Bauschinger effect present during reversal of the stress field which occurs between the compression and inflation stages. Bauschinger effect describes that pre-straining in any direction, as defined by the principal axis of the strain tensor, will introduce an anisotropy for further deformation in any other direction [15]. The intensity of this pre-strain associated anisotropy is at maximum when the direction of further straining is opposite to that of the pre-strain. A finite element test model was also created in order to verify that the elastic-plastic material model used in the analysis accurately represents the supplied material data. Generally, inflation of the stent may be performed using direct pressure applied to inside diameter of the stent or manually expanding the stent through prescribed boundary conditions. Using MSC.NASTRAN nonlinear finite element analysis (implicit solution) (MSC. Software. Nastran, CA, USA) will analyze the stent with rhombus structure shown in Figure 1.
In the model plastic hinge will be imposed at each connecting joint. Stent expansion is analyzed assuming struts within the stent will expand at a constant rate via a plastic hinge mechanism. The strut unit is the analytical model considered and illustrated as presented in Figure 2.
The modeling approach of expansion considers a bending stent strut as pairs of cantilevers loaded by a point force at the connection part, applied perpendicular to the initial cantilever position. During the expansion procedure, the eight or sixteen struts (depending on the stent in question) in a circle grow along the sphere are viewed as bending struts [15]. The circumferentially connecting parts are deformed without bending. Then the stent expands δ_z in the axial direction and δ_c in the circumferential direction (Figure 2a). Furthermore, the connection region with plastic hinge on the struts represents radius of curvature (r) and the deformed length of strut ( l^’) in the plastic region (Figure 2b).
The analysis of the stent inflation requires a model capable of handling large deformations and large strains. In addition, nonlinear material properties must be considered to accurately capture the deformation process. Therefore, the Lagrangian and the Eulerian descriptions are considered for the large deformations of the stents [17]. The Lagrangian or material approach uses the initial coordinates of each particle to describe the deformation whereas the Eulerian or spatial description uses the coordinates of the particles in the deformed state to describe the deformation, when elasticity effects is prominent [17]. Here we are using the Lagrangian approach. The stress and strain calculations are based upon true (Cauchy) stress and true (logarithmic) strain measures [8]. The axial stress, circumstantial stress and radial stress are considered present under anatomical bending. The artery wall pressure depends on the blood properties. Arterial pressure of the systemic circulation in terms of the systolic pressure over diastolic pressure is represented as compressive stress [20]. During the nonlinear analysis, Lagrangian method is considered with respect to the un-deformed configuration. The inflation of the stent which causes large deformation is modeled using an elastic-plastic material model, and inelastic strains. Solution of nonlinear equations, t+ΔtR is externally applied loads, t+ΔtF is nodal point forces corresponding to internal element stresses presented by the relationship in equation 1 [17]:
t+ΔtR- t+ΔtF=0 (1)
Equation 1 is set to equal equation 2 through the relationship:
0tk∙Δu=t+ΔtR- t+ΔtF(i+1)=0 (2)
where the iteration, i=1,2,3,… 0tk means the stiffness matrix of the structure at initial state, and Δu is the displacement. At convergence, the compatibility is satisfied with stress-strain law and equilibrium of nodal point and local force resulting in the following equilibrium equation 3:
t+ΔtRR= t+ΔtF (3)
This relates to the updated Lagrangian method being used during inelastic deformations, and the stress and strain calculations are based upon true (Cauchy) stress and true (logarithmic) strain measures. The 3D model is meshed with quadratic 10-noded tetrahedral elements (TETRA10). Approximate mesh size = 0.001 mm with five or six elements across the strut width is defined. Total number of elements 35,777 and total number of nodes 100,113 for the 3D stent model analyzed. Strut angle 20°. The stent material is 316L stainless steel with properties listed in Table 1. . Boundary conditions are applied at each end. One end fixed which allows for axial displacement and the line along the total length of the stent, constrains displacement in transverse direction. The other end is free. 83KPa pressure is the load which is assumed upon stent being inserted in the artery and inflated to a maximum diameter of 4.5mm. This is done by pressurization of an elastomeric balloon material inserted inside the tube (i.e. stent) which mimics the inflation activity.
Despite that stent implants are efficient in restoring normal blood flow over time new blockage can form called restenosis [1]. Research studies based on clinical data have concluded that in addition to stent design, the geometry of the artery can contribute to the development of restenosis [4,21,22]. For example several conditions such as twisting, compression, elongation and flexion as shown in Figure 3 can be noted for artery. These conditions contribute to the stent deformation and performance over time and this is why mirroring real artery geometry for modeling accuracy is of great importance.
In this study IVUS images were obtained from UT Southwestern Medical Center in Dallas, Texas for a patient with specific conditions and using slicing function and custom developed code at UNT 3D artery model was constructed along with the type of blood flow in the specific segment of the artery as shown in Figure 4. We can see from Figure 4 that the real artery has special contours and curves along its walls that vary from typical cylindrical shapes with uniform geometry used in many modeling studies.
The rheological behavior of blood is quantified by non-Newtonian viscosity. The computational blood flow model is taken as three-dimensional and transient flow in deforming blood vessels. Newtonian fluid viscosity is always constant against the shear rate. By comparison, non-Newtonian fluid viscosity changes depending on the shear rate. Literature demonstrates the relationship between blood viscosity and the shear rate for Newtonian and four different non-Newtonian models [23]. For our approach popular Carreau non-Newtonian model [23] is selected with features zero-shear plateau at low rates and a power-law region at high rates as demonstrated by the equation 4:
μ(|γ|)=μ∞+(μ0-μ∞)[1+(λγ)2](n-1)/2) (4)
μ0 = 0.056 Pa∙s viscosity at zero shear rate,
λ=3.131 time constant, and
n=0.3568 power law index.
Reynolds number is set as 300 for incompressible, homogeneous and sinusoidal pulsatile flow and artery wall is viewed as rigid body. The Carreau model describes pseudoplastic flow with asymptotic viscosities at zero μ0 (Newtonian) and infinite μ∞(n< 1 non-Newtonian) shear rates and with no yield stress; λ is a constant with units of time , after λ=1/γ, non-Newtonian behavior is important; n=0.3568 is power law index as a measure of the deviation of the fluid form Newtonian (n=1, Newtonnian fluid; n>1, shear-thicking(dilatant fluid); n< 1, shearthinning (pseudo-plastics)) (Ewoldt). Following steps 1- 7 in supplementary data section the average velocity and the shear stress are derived in the control volume.
Figure 5 shows results of expanded stent as the balloon is inflated, and each of the struts open. The expandable stage is noted by the green configuration and when the stent is released and allowed to spring back to a ‘free state’ noted by the blue configuration. Cross-sectional view of expanded versus initial diameter is shown in the upper right corner of Figure 5. The equivalent stress was calculated to be around 498MPa shown on the contour of model at around 50% expansion in diameter (i.e. from 3.0mm to 4.5mm).
Figure 6illustrates close view of the expanded configuration with respect to the maximum equivalent stress value of the structure, of around 498MPa. The plastic regions (shown with yellow and green color in Figure 6) are concentrated along each hinge joint as expected and predicted by the analytical model shown in Figure 2b and again in Figure 6 to the right.
Using rheological equations of blood (presented as equations 4 in this paper and 1-7 from supplementary data section), the viscosity versus shear rate relationship is predicted by calculating the following parametric values summarized in Table 2.
Figure 7aillustrates the velocity distribution of the blood flow in artery cross-section and the maximum velocity at the centerline of the artery where the contour showing approximately 0.685 m/s. The blood viscosity distribution in the artery is also shown in Figure 7b.Figure 7c illustrates calculated pressure gradient value of approximately 1.5 kPa/m to correspond to the ideal mean velocity of approximately 0.35m/s also listed in Table 2. a) b)
During stent deformation analysis it is noticed that the dominant change of loading appears to cause an increase tensile load in the connection parts. The stress is primarily caused as a result of bending as the struts open and straighten out during stent inflation. The maximum equivalent stress of 498 MPa occurred at the 4.5mm expansion diameter and was found to be less than the stainless steel ultimate strength of 515 MPa indicating a safe stent design throughout the expansion range. The true-stress values are typically larger than the respective engineering stress measured values, which is commonly used with material test data. This must be considered when comparison of maximum stress results reported in this study are to be made with engineering test data. Generally, the experimental data from a uniaxial tension test is expressed in terms of true stress vs. true strain, not engineering stress or strain. The maximum stress values were found to occur in localized regions of the stent. These regions were generally found along the inner radii of each of the ‘corner’ links connecting each of the longitudinal struts. Stress values throughout the whole stent were typically much lower. As to be expected, the stress in each of the radius sections exceeded the elastic limit of the material (207MPa), resulting in plastic (or permanent) deformation of the stent upon removal of the loading. If the stress not exceeded the elastic limit, removal of the load would result in the stent recovering its initial configuration. It was also noticed that the apparent blood viscosity decreases as the flow rate increases. The greater the flow is, the greater the rate is where one lamina of blood shears against an adjacent lamina of blood. The greater tendency for blood accumulation in the axial lamina at higher flow rates is partly responsible for this non-Newtonian behavior. However, it is possible that at very slow flow rates, the suspended cells tend to form aggregates, which increases blood viscosity. As flow is increased, this aggregation decreases, and so does the viscosity of blood.
In summary, the study demonstrated that the geometric parameters and structure can play main factors in the performance of the stent implant. A fundamental understanding of the mechanics of the unit shape involving compression/inflation activity and their modeling analysis give important insights for optimize stent designs. Combination of clinical performance reports with FEA modeling approach can simulate future real physical situation, using refined numerical studies and more complex cases such as stent-artery or stent-flow interactions. Using the insights gained from this study will allow us to analyze stent mechanical behavior and life-cycle under artery and blood flow effect for future cases. The future work can lead into personalized patient care where computational models can match stent configuration and material suitable with the patient disease and health state.
This work was supported by internal funds from UNT. We thank the team at UTSW Medical Center for providing valuable feedback and stent samples for analyses as well as IVUS images.
The rate-of-deformation tensor D and the shear rate γ :
D=1/2 ((∂u_j)/(∂x_i )+(∂u_i)/(∂x_j )) (1)
γ =√(2D:D ) (2)
The shear stress is simply written in terms of a non-Newtonian viscosity μ(γ):
τ=μ(γ)D (3)
In Carreau model, considering u=0,v=0, velocity won’t change, ∂u/∂x=0,∂v/∂y=0,∂w/∂z=0.
Then, the momentum balance equation becomes: ρ (∂v)/∂t+ρ(v ∙)v =-P+∙τ (4)
Momentum in x-direction is derived: ρu ∂u/∂x+ρv ∂u/∂y =-∂P/∂x+∂/∂x (μ ∂u/∂x)+ ∂/∂y (μ ∂u/∂y)
Momentum in y-direction is derived: ρu ∂v/∂x+ρv ∂v/∂y=-∂P/∂y+∂/∂x (μ ∂v/∂x)+ ∂/∂y (μ ∂v/∂y)
Momentum in z-direction is derived: ρu ∂w/∂x+ρv ∂w/∂y=-∂P/∂z+∂/∂x (μ ∂w/∂x)+ ∂/∂y (μ ∂w/∂y)
∂P/∂z≠0 Then, equation is simplified as: ∂P/∂z=∂/∂x (μ ∂w/∂x)+ ∂/∂y (μ ∂w/∂y)
γ=√((∂w/∂x)^2+(∂w/∂y)^2 ) (5)
Then, to compute the blood flow through two ends constricted artery, the solver is implemented to solve sequentially the momentum balance equation. The coupling between velocity and pressure drop is achieved through algorithm. The continuity control volume (i.e. cross-section of the artery) is analyzed as discretized domain shown in Figure 8a. For 2D model, the viscosity at the control volume faces is evaluated precisely by setting the grid points and draw control volume faces between neighboring points as shown in Figure 1b.