 Research
 Open Access
 Published:
Multiobjective optimization of coronary stent using Kriging surrogate model
BioMedical Engineering OnLine volume 15, Article number: 148 (2016)
Abstract
Background
In stent design optimization, the functional relationship between design parameters and design goals is nonlinear, complex, and implicit and the multiobjective design of stents involves a number of potentially conflicting performance criteria. Therefore it is hard and timeconsuming to find the optimal design of stent either by experiment or clinic test. Fortunately, computational methods have been developed to the point whereby optimization and simulation tools can be used to systematically design devices in a realistic timescale. The aim of the present study is to propose an adaptive optimization method of stent design to improve its expansion performance.
Methods
Multiobjective optimization method based on Kriging surrogate model was proposed to decrease the dogboning effect and the radial elastic recoil of stents to improve stent expansion properties and thus reduce the risk of vascular instent restenosis injury. Integrating design of experiment methods and Kriging surrogate model were employed to construct the relationship between measures of stent dilation performance and geometric design parameters. Expected improvement, an infilling sampling criterion, was employed to balance local and global search with the aim of finding the global optimal design. A typical diamondshaped coronary stentballoon system was taken as an example to test the effectiveness of the optimization method. Finite element method was used to analyze the stent expansion of each design.
Results
27 iterations were needed to obtain the optimal solution. The absolute values of the dogboning ratio at 32 and 42 ms were reduced by 94.21 and 89.43%, respectively. The dogboning effect was almost eliminated after optimization. The average of elastic recoil was reduced by 15.17%.
Conclusion
This article presents FEM based multiobjective optimization method combining with the Kriging surrogate model to decrease both the dogboning effect and radial elastic recoil of stents. The numerical results prove that the proposed optimization method effectively decreased both the dogboning effect and radial elastic recoil of stent. Further investigations containing more design goals and more effective multidisciplinary design optimization method are warranted.
Background
As the leading cause of mortality, cardiovascular disease is often related to atherosclerosis which caused by the progressive formation of plaque and eventually results in an obstruction (stenosis) for blood flow through the artery [1–5]. Compared to traditional treatments such as drugs and surgery for coronary artery diseases (narrowing or blockage of the coronary arteries), percutaneous transluminal coronary stenting with the aid of coronary balloon angioplasty is more widely adopted in clinical practice thanks to its high initial success rate, minimal invasive nature, and improved longterm effectiveness. A stent is a wire metal meshed tube placed in the vessel during coronary balloon angioplasty to offer radial strength and to overcome the acute elastic recoil. The stent is put over a balloon catheter and moved into stenosis segment. Then, it expands as the balloon is inflated to open the blocked vessel. After the balloon and catheter are removed, the stent remains in the vessel to act as a scaffold to help prevent arteries from becoming narrowed or blocked again. Nowadays, intravascular stents are routinely and successfully used in medical treatment, but it still needs to be improved. For example, instent restenosis remains the main obstacle for the development of stent. It is known that instent restenosis is caused by artery injury due to stent expansion and vascular inflammation to the stent struts. Therefore, scholarly efforts to improve stent expansion performance and reduce the injury of blood vessel caused by stent implantations in stent design optimization are of great importance.
A desirable stent should possess a number of excellent mechanical properties, including smaller dogboning ratio and smaller radial elastic recoil. The dogboning phenomenon caused by nonuniform balloonstent expansion has a significant impact on the development of thrombus and intimal hyperplasia [6, 7]. A larger dogboning ratio indicates a more serious warpage at the ending struts, which will cause mechanical damage to the vessel wall and results in instent restenosis [8–10]. Additionally, the radial elastic recoil due to elastic deformation of stent has a significant impact on the mechanical support of stent. It is believed that the stent design may affect stent expansion performance such as the dogboning ratio and radial elastic recoil. Thus, it is important to predict the dogboning effect and radial elastic recoil and to optimize the design before manufacturing the stent.
Computational simulation (e.g., finite element analysis (FEA)) can be a very useful tool to study the stent expansion [11–14]. Dumoulin and Cochelin [11] evaluated and characterized the mechanical properties and behaviors of a balloon expandable stent. Etave et al. [15] compared the mechanical performance of two types of stents. In terms of stent design, Migliavacca et al. [16, 17] and Beule et al. [1] assessed the mechanical properties and behavior of balloon expandable stents to determine how the FEA method could be used to optimize stent designs. It is easy to study the mechanical properties and analyze the effective factors, but it is difficult to find the globally optimal solution since the functional relationship between the geometrical parameters and dilation performance of stent is complex, nonlinear and implicit. For the traditional stent design method, a finite number of different designs are assessed and compared with each other to find an optimal design. This method is often used in industry when commercial demands restrict the time spent on developing a better product. However a limited number of discrete points fail to represent all the information in the design space and thus it is very difficult to find the optimal design through the traditional methods. For traditional methods, the designer has to formulate how many different designs can be tested in the available time and their work is also affected by the design methodology, problem fidelity, parameterization and measures of performance. The focus of the study is on the computational design of stent. However, since the stent is minute and boundary conditions for the expanding of stent in the vessel are complicated, it is relatively hard to apply finite element modelling in the optimization of stent. Actually not only the dependences mentioned above but also the availability and power of computers and software shall be considered when the traditional design methods are adopted.
Consequently, some approximation models are widely used in engineering to construct simplified approximations for analysis codes, especially when the analysis is hard and timeconsuming. Incorporating a shape optimization algorithm based on a proven convergence theory into the design process allows engineers to systematically identify the most favourable designs. An adaptive optimization method based on Kriging surrogate model was already proposed to eliminate the dogboning phenomenon by Li et al. [18]. A derivativefree optimization algorithm coupling computational fluid dynamics (CFD) was used for stent design by Gundert et al. [19]. In this paper, Kriging models were used as alternatives to the method of traditional secondorder polynomial response surfaces for constructing global approximations in stent optimization. As a semiparametric approach, the Kriging model [20, 21] is much more flexible than approaches based on parametric behavioral models.
Taking the consideration above in mind, we adopted the Kriging model to create an approximate functional relationship between the design objective and design parameters to replace the expensive reanalysis of the stent dogboning ratio and radial elastic recoil. The optimization iterations are based on the approximate relationship between the design objective and design parameters to reduce the high computational cost. An adaptive optimization method based on the Kriging surrogate model combing with modified rectangular grid (MRG) approach was proposed to minimize the radial elastic recoil and the dogboning effect of stent during the expansion process. Expected improvement (EI) function is employed in the adaptive process [18], which can balance local and global searches and then find the global optimal design even with a small sample size. The FEA solver of ANSYS was used to analyze the measurements of stent expansion performance.
Methods
Finite element model
A typical PalmazSchatz stent (shown in Fig. 1) was investigated in this study. The geometries and loading method of the stent was supported by Ref. [22]. A balloon that is 11.4 mm in length and 0.12 mm in thickness was placed inside the stent. The outer wall of the balloon is close to the inner wall of the stent. Geometric dimensions of the stent are shown in Fig. 1. A timerelated pressure (shown in Fig. 2) was loaded on the inner surface of balloon to stimulate the expanding process of the balloon stent system.
Bilinear elastic–plastic and hyperelastic (MooneyRivlin) materials were assumed for slotted tube stents and balloon. Data of the material properties used in this study was from previous studies [23, 22].
The balloonstent dilatation involves nonlinearities namely contact, elastoplasticity and large deformation. Thence, “solid 185” in ANSYS was used for threedimensional modeling of the stent. The element possesses plasticity, large deflection, large strain capabilities and mixed formulation capability for stimulating deformation of nearly incompressible elastoplastic materials. Shell 181 was used to model the balloon since it is wellsuited to be applied in cases that involve large strain nonlinear. Due to the symmetry of the entire structure and the loading pressure only, 1/16 of the model (1/8 in circumferential direction and 1/2 in longitudinal direction) was preformed to analyze the dogboning ratio and radial recoil of stent, as shown in Fig. 3. The balloon was modelled as a hyperelastic shell. The nodes at the distal end of balloon were constrained without rigid body displacement, while the nodes at the distal end of stent were free. Symmetry boundary conditions were applied to the symmetry parts of the stent and balloon. Pressure involving three load phases (shown in Fig. 2) was applied to the inner surface of the balloon. A facetoface contact between balloon and stent inner surface was considered and the friction between them was ignored. The stent expanded as the balloon was inflated, as shown in Fig. 4. The stent was expanded to a large permanent deformation due to the expansion of balloon and the plastic strain occurs in most parts of struts. Thence, the stent can stay in the stenotic artery permanently and holds the artery wall open to prevent restenosis or narrowing of coronary arteries after the balloon is deflated and withdrawn. In this study, the stent outer diameter of each design was dilated to a same diameter of 4.54 mm. Obviously, the pressures were varied with different stent geometries. The binarysearch method was adopted to find the pressure of each design of stent to dilate the proximal ends of it (marked in Fig. 3) to the nominal diameter after unloading of the balloon. This is to allow the stenotic segment to be opened in agreement with a health artery (diameter 4.54 mm in this study) after stent dilation.
Optimization problem
Generally, the dogboning effect exists throughout the expanding process and usually reaches its maximum in the beginning of loading [24], but the struts have no direct contact with the vessel wall. While, during the stent deformed from 25 to 32 ms, the dogboning ratio was relatively large [18] in this study and it could cause injury to the vascular wall. Therefore, the dogboning ratio of the stent can be defined as:
where \(d_{\text{radial}}^{\text{distal}}\) and \(d_{\text{radial}}^{\text{proximal}}\) denote the distal and proximal radial displacements respectively when it is at 32 ms. It is easy to find that when the distal of stent is overexpanded, DR is larger than 0; while when the proximal of stent is overexpanded, DR is less than 0. Therefore, in order to expand the stent uniformity along the longitudinal direction, we need to minimize the absolute value of DR.
For the typical diamond stent, bridge struts provide the main radial support of stenosis. Furthermore, the radial displacement and radial recoil of stent in these bridge struts are different only in the case that there is a uniform expansion of stent along its longitudinal direction. Therefore, four points on the bridge struts were chosen as the points of interest (shown in Fig. 1), and the average of radial recoil at these four points was taken as stent radial elastic recoil. Therefore, the stent radial elastic recoil can be defined as:
where \(R_{{{\text{POI(No}}.i )}}^{\text{loading}}\) and \(R_{{{\text{POI(No}}.i )}}^{\text{unloading}}\) are the radius of the stent at the point of interest (POI) (shown in Fig. 1) when it is at 32 and 42 ms respectively. Obviously, RER is the average of the radial elastic recoil at the POIs.
Our optimization objective is to find a set of design variables to reduce both DR and RER. A common approach in multiobjective optimization is to compute the weighted sum of all the objective functions. Therefore, the multiobjective problem can be transformed into a singleobjective problem and the objective function can be considered as:
Because the two objectives are mutually incompatible, it is not easy to choose appropriate weights ω _{1} and ω _{2}. Moreover, the two objectives have different scales. If we scale both DR and RER to the same range namely [0, 1], then we might be able to assign some reasonable weights. In order to do this, we need to improve the objective function. And the optimization of the coronary stent can be defined as follows:
where ω ∈ (0, 1), DR_{min} and DR_{max} are the minimum and maximum of DR in the samples respectively, RER _{min} and RER _{max} are the minimum and maximum of RER in the samples respectively. x is the design variables consisting of the geometrical parameters W _{ i }, i= 1,…, 5 and T of stent and L (the length of balloon), as shown in Fig. 1. In this study, ω = 0.5.
Optimization method
A finite element based multiobjective optimization method combining with Kriging surrogate model [18] was constructed for the stent optimization to improve stent expansion performance. Kriging was used to build the approximate functional relationship between the design objective and design variables. A modified rectangular grid (MRG) approach was adopted to generate the initial sample for Kriging. EI function was employed to balance the local and global search to find the global optimal design.
Kriging approximate method
The Kriging model is regarded as a method of functional modeling of a stochastic process. Therefore it is named as the “stochastic process model”, which is written as:
where, x ^{i} = {x ^{i}_{1} , x ^{i}_{2} ,…, x ^{i}_{ m } } denotes the ith sample point with a number of m variables; \(\hat{y}({\mathbf{x}}^{i})\) is an approximate function fitted to a number of n sample points; f(x ^{i}) is a linear or nonlinear function of x ^{i}; β is the regression coefficient to be estimated; and z(x ^{i}) is the stochastic function, with a mean of zero and a variance σ ^{2}. The spatial correlation function between stochastic functions is given by
where R(θ, x ^{i}, x ^{j}) is the Gaussian correlation function with θ, which characterizes the spatial correlation between two samples. Parameters can be estimated by maximizing the likelihood of sample points.
where f = [f _{1}, f _{2}, …, f _{ n }]. The estimates \(\hat{\beta }\) and \(\hat{\sigma }^{2}\) can be obtained by Eq. (7).
Predictor
A linear combination of the response values of sample Y could be used to estimate \(\hat{y}({\mathbf{x}}^{*} )\) of a new point x*
The mean squared error (MSE) of this predictor can be minimized by unbiased estimation, which gives
where
Therefore, the function value \(\hat{y}\left( {{\mathbf{x}}^{ * } } \right)\) at every new point x* can be predicted by using Eq. (9).
Sampling strategy
MRG approach was adopted to generate the sample points for constructing the Kriging model. The range of m design variables was defined as \(l_{j} \le x_{j} \le u_{j} ,\,\,j = 1,\, \ldots ,\,m\). The number of levels in the jth dimension is q _{ j }. Then the approach is performed as follows:

1.
Narrow the range of variables as
$$l_{j} \le x_{j} \le \hat{u}_{j} ,\,\,\hat{u}_{j} = u_{j}  \frac{1}{2}\frac{{u_{j}  l_{j} }}{{q_{j}  1}},\quad j = 1,\, \ldots ,\,m$$(11) 
2.
Perform rectangular grid (RG) sampling [25] in the narrowed space as
$$x_{j}^{i} = l_{j} + k_{j}^{(i)} \frac{{\hat{u}_{j}  l_{j} }}{{q_{j}  1}},\quad k_{j} = 0,\,\,1,\, \ldots ,\,q_{j}  1 \quad i = 1,2,\, \ldots, \, \prod\limits_{j = 1}^{m} {q_{j} }$$(12) 
3.
Add a stochastic movement of each sample point in each dimension as
$$\frac{{\alpha_{ij} }}{2}\frac{{u_{j}  l_{j} }}{{q_{j}  1}} \quad j = 1,\,2,\, \ldots ,\,m$$(13)
where \({\kern 1pt} {\kern 1pt} {\kern 1pt} \alpha_{ij} \in [ 0, 1]\), which is assumed to be normally distributed.
Compared to RG, MRG has several advantages such as preventing sample points from lying in boundary, which can provide more useful information for constructing the Kriging model. What’s more, MRG can ensure that the points have lower pairwise correlation between the factors and avoid the case of sample points spaced too close to each other. The distance between two arbitrary points must satisfy
Expected improvement (EI)
Generally, the response surface based optimization is to find the minimum of the response surface. But this method often results in a local minimum, even if iterations are performed in the search. Fortunately, an “expected improvement (EI)” function can be used to balance local and global search, which regarded as an effective global optimization (EGO) [21]. The EI method computes the extent of improvement of response value. For any point x which is not one of the already known sample points in the design space, the value of Y(x) is unknown. Thus, Y(x) can be considered as a random variable and assumed it is normally distributed with a mean \(\hat {\mathbf{y}} ({\mathbf{x}})\) and variance σ^{2} got from the Kriging predictor. If the current best optimization function value is Y _{min}, and Y(x) = Y _{min} − I, then an improved I will be obtained. The likelihood to achieve such an improvement is given by the normal density function
The expected improvement is the expected value of the improvement obtained by integrating over the following density:
With integration by parts, it can be obtained
where \(\Phi\) and \(\phi\) denote the normal cumulative distribution and density functions respectively. And
where Eq. (18) can be meaningless if σ(x) equals zero. Hence, it can be written as
The first term of Eq. (19) refers to the difference between the current minimum response value Y _{min} and the prediction \(\hat{y}\left( {\mathbf{x}} \right)\) at x. Hence, it is large when \(\hat{y}\left( {\mathbf{x}} \right)\) is small. The second term is the product of the root mean squared error (RMSE) σ(x) and the normal density function \(\phi (u)\). \(\phi (u)\) is large when σ(x) is large and \(\hat{y}\left( {\mathbf{x}} \right)\) is closed to Y _{min}. Thus, the expected improvement will be larger when the predicted value is smaller than Y _{min} and/or there is a lot of uncertainty associated with the prediction.
The convergence criterion
The optimization iteration stops when
where Δ_{1}, Δ_{2} and Δ_{3} are the convergence tolerances. Y _{max} and Y _{min} are the maximal and minimal function values of sample points respectively. f _{ k } and f _{ k−1} are the objective function values at the kth and k − 1th iteration, respectively. This stopping criterion enjoys the advantage that the user can set the “relative” tolerance Δ_{1} without considering the magnitude of the problem response.
Implementation of optimization procedure
The optimization algorithm for a coronary stent based on Kriging model combining DOE methods and EI function can be described as follows:
 Step 1:

Get a set of n _{ s } samples using MRG.
 Step 2:

Run ANSYS program with Binarysearch method to dilate stent at sample point i, i = 1,…, n _{ s } to nominal diameter and obtain \(d_{\text{radial}}^{\text{distal}}\) and \(d_{\text{radial}}^{\text{proximal}}\), \(R_{{{\text{POI(No}}.i )}}^{\text{loading}}\) and \(R_{{{\text{POI(No}}.i )}}^{\text{unloading}}\), i = 1, …, 4. Calculate f(x _{ i }) in problem (3) at the sample point i, i = 1, …, n _{ s. }
 Step 3:

Find the sample point with the minimum f(x _{ i }) as the initial point for the optimization.
 Step 4:

Get an approximate functional relationship between the design objective f(x) and design variables using Kriging surrogate model based on the trial samples. Calculate EI(x) based on f(x).
 Step 5:

Select optimization algorithm to implement the optimization design based on max EI and obtain the modified design x _{ k. }
 Setp 6:

Get the predictive value \(\hat{y}_{k}\) of x _{ k } based on Kriging and compute f(x _{ k }) by ANSYS program.
 Step 7:

The optimization iteration was stopped when a suitable level of convergence is reached and/or the available time for the optimization process is exhausted. The process of constructing and maximizing EI does not stop until the Euclidean norm between real value f(x _{ k }) and predictive value \(\hat{y}_{k}\) fall below a given tolerance, the Euclidean norm between current and previous iterates falls below a given tolerance, and the criterion stipulated in “The convergence criterion” section is reached. If not, then add the modified design into the set of samples and go to step 3.
Results
The absolute value of dogboning ratio and radial elastic recoil of PalmazSchatz stent were minimized by the multiobjective optimization method proposed in this paper. The initial trial samples which included the initial design and 30 samples generated by MRG (listed in Table 1) were selected to construct the Kriging surrogate model. The dogboning ratio and radial elastic recoil of stents for all trial samples are simulated by finite element method. EI function was adopted to balance local and global search in the design space. The optimization process started from the initial point which has the minimum value of design objective among all the sample points. 27 iterations were needed to obtain the optimal solution as shown in Fig. 5.
Optimization results in details
The optimization result was compared to the original design as shown in Table 2. The absolute value of the dogboning ratio at 32 ms was reduced by 94.21%, which indicates that the dogboning effect was almost eliminated. Moreover, although the absolute value of dogboning ratio after unloading at 42 ms was not considered in the optimization objective, it was reduced by 89.43%. The average of radial elastic recoil was reduced by 15.17%. Because the radial elastic recoil is not only related to the stent’s structure, but also connected with the materials and expansion process of the stent, it is hard to eliminate the radial elastic recoil completely. Both the dogboning and radial elastic recoil are two important features to evaluate the stent expansion performance, but some factors that influence them are contradictory to each other. From the optimal result, we can see that the optimal stent with the greater W _{5} results in higher radial force, which means that the optimal stent can better support the artery wall. This explains the decrease of radial elastic recoil of optimal stent. But the stent with a higher radial force is hard to be expanded. Generally, the stent with a higher radial force is hard to be dilated and the ends of it will open first during expansion. This phenomenon is referred to as the dogboning effect. However, the smaller W _{2}, W _{3}, W _{4} and T together with appropriate W _{1} and L result in lower dogboning, although they lead to higher radial elastic recoil. It shows our proposed multiobjective optimization method can effectively find a set of design variables that minimizes both dogboning and elastic recoil.
The radial displacement distributions of the original and optimal stents at 32 ms, which symmetrical displayed in the circumferential direction, are shown in Fig. 6. The diameters of the original and optimal stents at the proximal ends were dilated to a same diameter of 4.54 mm after the deflation of balloon. The radial displacement of the original stent at the distal ends was much larger than that at the proximal ends. While, the proximal and distal radial expansions of the optimal stent were similar and the dogboning ratio was almost 0, which indicates uniform expansion along the length of optimal stent.
Improvement of stent expansion process
The proximal and distal radius of both the original and optimal stents during the dilation process are shown in Fig. 7. The proximal radius of all the stents located at the sample points were expanded to the same nominal radius (2.27 mm) after the deflation of balloon. Figure 7 shows that the difference of the radius of optimal stent at the proximal and distal ends was smaller than that of the original one, particularly in the period from 25 to 32 ms. This indicates a uniform dilation of the optimal stent along its length. The reduction of the stent’s radius during the period from 25 to 42 ms demonstrates the radial elastic recoil of stent. Because the elastic recoil is not only related to the stent structure, but also related to the materials and expansion process, it wasn’t decreased completely in this study.
Figure 8 shows the dogboning ratio for original and optimal stents along with time of stent dilation. The dogboning effect reached its maximum at the prophase of loading stage, and was reduced and remained in an almost constant value after stent expansion from 25 to 32 ms of loading and from 32 to 42 ms of unloading. The radial of stent during these periods reached its maximum, thereby the contact between stent and artery wall expanded and intensified. The dogboning effect during these periods would cause serious instantaneous mechanical damage to the blood vessel. The maximum DR of the optimal stent was 8.78, while the maximum DR of the original stent was 13.9. It means that the maximum DR has been reduced by 36.83% after optimization. Moreover, the time of maximum DR of optimal stent is earlier than that of original stent, which is helpful to reduce the risk of mechanical damage caused by stent on vessel wall. During the period from 25 to 42 ms, DR of the optimal stent was almost 0, while DR of the original stent was about 0.17. DR of the optimal stent in this period is decreased, and it denotes smaller instantaneous mechanical damage caused by stent on vessel wall.
Discussions
A FEM based multiobjective optimization method combining with the Kriging surrogate model is proposed to reduce the dogboning effect and the radial elastic recoil of stent to improve the stent expansion performance. Our results show that the proposed optimization method could be used for stent design optimization effectively and conveniently. This provides a new method of stent design and represents a new direction of research. This optimization method combined with experimental verification can serve as a useful tool for stent design before manufacture.
In contrast to the expansive computational simulations employed in the comparison test studies [11, 12, 13, 15, 16, 17], the surrogate modeling approach in which response surface models (RSMs) were used to represent the relationship between design objectives and design variables [27]. Whilst most studies of stent design relate to multiple objectives, some articles only dealt with a single objective function. Harewood et al. [26] focused on radial stiffness of a single ring. Li et al. [28, 29] optimized stent dogboning and drug release, respectively. Grogan et al. [27] performed a single objective optimization for maximum radial strength. When considering multiple objectives, Pant et al. [30] and Bressloff [31] conducted FEA simulation to generate a range of multidisciplinary objectives. Pant et al. [30] constructed the Pareto fronts generated by treating each objective separately. Bressloff [31] recast the optimization as a constrained problem, wherein design improvement is sought in one objective while other objectives were considered as constraints. Multiobjective optimization of stent design involves a large number of design goals. It is difficult to find the optimal solution to improve all of them just by one of the methods to solve multiobjective problem, such as combining the design objectives in a single weighted objectives function, searching the Pareto fronts, and taking same design objectives as constraints. In future work, these methods can be used in combination under the premise of rational planning of design objectives and design variables of stent optimization systems, including stent auxiliary expansion, instent blood flow, drug release, and biomechanical response of vascular tissue, to improve the performance of stenting.
Some limitations of this study include: (a) FEA model of stent dilation does not contain blood vessels and thrombosis, (b) Balloon folding is not considered during its expansion process, (c) The results of optimal design have made some improvement of stenting performance, but it’s yet short of enough validation through experiment.
Conclusions
This article presents a FEM based multiobjective optimization method combining with the Kriging surrogate model to decrease both the dogboning effect and radial elastic recoil of stents. The Kriging surrogate model coupled with DOE methods was adopted to construct an approximate functional relationship between the objective function and geometries. The EI function was employed to balance local and global searches with the aim of finding out the global optimal design. The proposed optimization method effectively decreased both the dogboning effect and radial elastic recoil of stent. More issues of stent design should be considered and more effective multidisciplinary design optimization method should be investigated to continue our study.
References
 1.
Beule MD, Impe RV, Verhegghe B, Segers P, Verdonck P. Finite element analysis and stent design: reduction of dogboning. Technol Health Care. 2006;14:233–41.
 2.
Tang D, Yang C, Kobayashi S, Zheng J, Vito RP. Effect of stenosis asymmetry on blood flow and artery compression: a threedimensional fluidstructure interaction model. Ann Biomed Eng. 2003;31:1182–93.
 3.
Tang D, Yang C, Kobayashi S, Ku DN. Effect of a lipid pool on stress/strain distributions in stenotic arteries: 3D fluidstructure interactions (FSI) models. J Biomech Eng. 2004;126:363–70.
 4.
Qiao A, Liu Y, Guo Z. Wall shear stresses in small and large twoway bypass grafts. Med Eng Phys. 2006;28:251–8.
 5.
Qiao A, Liu Y. Numerical study of hemodynamics comparison between small and large femoral bypass grafts. Commun Numer Methods Eng. 2008;24:1067–78.
 6.
Lim D, Cho SK, Park WP, Kristensson A, Ko JY, Hassani STS, Kim HS, AlHassani STS, Kim HS. Suggestion of potential stent design parameters to reduce restenosis risk driven by foreshortening or dogboning due to nonuniform balloonstent expansion. Ann Biomed Eng. 2008;36:1118–29.
 7.
Mortier P, Beule MD, Carlier SG, Impe RV, Verhegghe B, Verdonck P. Numerical study of the uniformity of balloonexpandable stent deployment. J Biomech Eng. 2008;130:210–8.
 8.
McClean DR, Eigler NL. Stent design: implications for restenosis. Rev Cardiovasc Med. 2002;3:16–22.
 9.
Rogers C, Tseng DY, Squire JC, Edelman ER. Balloonartery interactions during stent placement: a finite element analysis approach to pressure, compliance, and stent design as contributors to vascular injury. Circulation. 1999;84:378–83.
 10.
Schulz C, Herrmann RA, Beihzrz C, Pasquantonio J, Ait E. Coronary stent symmetry and vascular injury determine experimental restenosis. Heart. 2000;83:462–7.
 11.
Dumoulin C, Cochelin B. Mechanical behaviour modelling of balloonexpandable stents. J Biomech. 2000;33:1461–70.
 12.
Tan LB, Webb DC, Kormi K, AlHassani STS. A method for investigating the mechanical properties of intracoronary stents using finite element numerical simulation. J Cardiol. 2001;78:51–67.
 13.
Lally C, Dolan F, Prendergast PJ. Cardiovascular stent design and vessel stresses: a finite element analysis. J Biomech. 2005;38:1574–81.
 14.
Holzapfel GA, Stadler M, Gasser TC. Mechanics of biological tissue. Heidelberg: Springer; 2005. p. 207–20.
 15.
Etave F, Finet G, Boivin M, Boyer JC, Rioufol G, Thollet G. Mechanical properties of coronary stents determined by using finite element analysis. J Biomech. 2001;34:1065–75.
 16.
Migliavacca F, Petrini L, Colombo M, Auricchio F, Pictrabissa R. Mechanical behavior of coronary stents investigated through the finite element method. J Biomech. 2002;35:803–11.
 17.
Migliavacca F, Petrini L, Massarotti P, Schievano S, Auricchio F, Dubini G. Stainless and shape memory alloy coronary stents: a computational study on the interaction with the vascular wall. Biomech Model Mechanobiol. 2004;2:205–17.
 18.
Li H, Wang X. Design optimization of balloonexpandable coronary stent. Struct Multidiscip Optim. 2013;48(4):837–47.
 19.
Gundert TJ, Marsden AL, Yang W, LaDisa JF. Optimization of cardiovascular stent design using computational fluid. J Biomech Eng. 2012;134(1):011002.
 20.
Lophaven SN, Nielsen HB, Sondergaard J. DACE—a Matlab Kriging toolbox; version 2; informatics and mathematical modelling. Kongens Lyngby: Technical University of Denmark; 2002.
 21.
Jones DR, Schonlau M, Welch WJ. Efficient global optimization of expensive blackbox functions. J Global Optim. 1998;13:445–92.
 22.
Chua SND, Donald BJM, Hashmi MSJ. Finiteelement simulation of stent expansion. J Mater Process Tech. 2002;120:335–40.
 23.
Chua SND, Donald BJM, Hashmi MSJ. Effects of varying slotted tube (stent) geometry on its expansion behaviour using finite element method. J Mater Process Tech. 2004;155–156:1764–71.
 24.
Kiousis DE, Wulff AR, Holzapfel GA. Experimental studies and numerical analysis of the inflation and interaction of vascular balloon catheterstent systems. Ann Biomed Eng. 2009;37:315–30.
 25.
Bressloff N, Ragkousis G, Curzen N. Design optimisation of coronary artery stent systems. Ann Biomed Eng. 2016;44:357–67.
 26.
Harewood F, Thornton R, Sharp P. Step change in design: exploring sixty stent design variations overnight. 2011. http://www.altairproductdesign.com.
 27.
Grogan J, Leen S, McHugh P. Optimizing the design of a bioabsorbable metal stent using computer simulation methods. Biomaterials. 2013;34(33):8049–60.
 28.
Li H, Qiu T, Zhu B, Wu J, Wang X. Design optimization of coronary stent based in finite element models. Sci World J. 2013;2013:630243.
 29.
Li H, Zhang Y, Zhu B, Wu J, Wang X. Drug release analysis and optimization for drugeluting stents. Sci World J. 2013;2013:827839.
 30.
Pant S, Limbert N, Curzen N, Bressloff N. Multiobjectives design optimisation of coronary stents. Biomaterials. 2011;32:7755–73.
 31.
Bressloff N. Multiobjective design of a biodegradable coronary artery stent. Stud Mechanobiol Tissue Eng Biomater. 2014;15:1–28.
Declarations
Authors’ contributions
HL, BZ, MW, DZ and AQ were responsible for the design, data collection and overall investigation, and established the optimization model. JG and ZL were responsible for the optimization method. All authors contributed to the writing of the paper. All authors read and approved the final manuscript.
Authors’ information
Wang’s group (Wang, Zhao and Li) has been doing research in stent design, mechanical properties analysis, structure optimization and micromolding, see website: http://gs1.dlut.edu.cn/Supervisor/Front/dsxx/new/Default.aspx?WebPageName=WangMJ. http://gs1.dlut.edu.cn/Supervisor/Front/dsxx/new/Default.aspx?WebPageName=zhaodanyang. Zhu has been doing research in computational mechanics, finite element method and materials science: http://gs1.dlut.edu.cn/Supervisor/Front/dsxx/new/Default.aspx?WebPageName=zhubao. Gu and Li have been doing research in optimization algorithms, highperformance computing and largescale scientific: http://denm.dlut.edu.cn/info/1589/2878.htm. Qiao’s group has been playing a leading role in graft, stent and blood flow dynamics of aortic dissection: http://life.bjut.edu.cn/szdw/jcrc/20151231/14515449774055993_1.html.
Acknowledgements
This work was supported by the National Natural Science Foundation of China (Grant no. 11502044, 81171107) and Postdoctoral Science Foundation of China (Grant no. 2014M561222).
Competing interests
The authors declare that they have no competing interests.
About this supplement
This article has been published as part of BioMedical Engineering OnLine Volume 15 Supplement 2, 2016. Computational and experimental methods for biological research: cardiovascular diseases and beyond. The full contents of the supplement are available online http://biomedicalengineeringonline.biomedcentral.com/articles/supplements/volume15supplement2 .
Funding
Publication charges for this article have been funded by Postdoctoral Science Foundation of China (Grant no. 2014M561222).
Author information
Affiliations
Corresponding author
Rights and permissions
Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
About this article
Cite this article
Li, H., Gu, J., Wang, M. et al. Multiobjective optimization of coronary stent using Kriging surrogate model. BioMed Eng OnLine 15, 148 (2016). https://doi.org/10.1186/s1293801602689
Published:
Keywords
 Stent
 Dogboning
 Radial elastic recoil
 Blackbox techniques
 Kriging surrogate model
 Design optimization