Научная статья на тему 'An Algorithm for Automatic Generation and Evaluation of Leaflet Apparatus Models for Heart Valve Prostheses'

An Algorithm for Automatic Generation and Evaluation of Leaflet Apparatus Models for Heart Valve Prostheses Текст научной статьи по специальности «Медицинские технологии»

CC BY
85
29
i Надоели баннеры? Вы всегда можете отключить рекламу.
Ключевые слова
heart valve prosthesis / simulation / algorithm for generation and evaluation of leaflet apparatuses

Аннотация научной статьи по медицинским технологиям, автор научной работы — P.S. Onischenko, K.Yu. Klyshnikov, Е.А. Ovcharenko, L.S. Barbarash

The aim of the study is to develop and verify an algorithm for automatic generation of leaflet apparatus models for prosthetic heart valves, to optimize the basic parameters of the models in order to minimize the stress-strain state and maximize the geometric area of the orifice. Materials and Methods. The suggested algorithm consists of three blocks: “Generator”, “Modeling”, “Analysis”. The first block creates a three-dimensional model of the leaflet apparatus using the specified parameters (height, radius, thickness, degree of “sagging”, angle of the free edge deviation). Numerical simulation of the apparatus functioning is further performed using the finite element method. Then, the statistical analysis of the von Mises stresses is done and the opening area of the design in question is calculated. Verification was performed by comparing quantitatively the lumen areas of the leaflet apparatus in the open state, obtained from the literature data for the Trifecta bioprosthesis (19, 21, and 23 mm in diameter), with the results of the described algorithm operation. Results. The verification of the algorithm has demonstrated the following deviations in the lumen area in the open state: 2.85% for 19 mm, 14.81% for 21 mm, and 23.17% for 23 mm models. This difference is due to the choice of the model material (no data could be found on the physical and mechanical properties of the pericardium used for the fabrication of the Trifecta bioprostheses). The generation of a large number of designs (n=1517) without fixation of certain geometry parameters has shown that thickness of the leaflet apparatus makes the greatest contribution to the degree of opening; its dependence on the thickness and arising peak von Mises stresses has been demonstrated. Of the valvular models obtained, 278 showed the opening degree greater than 80% and maximum peak von Mises stresses below 4 MPa for the proposed model of the pericardium, which is 65% below the ultimate strength of the material. Out of 278 leaflet models, 3 “optimal” designs were selected meeting the diameter criteria of 19, 21, and 23 mm. The loss index for them was 0.24, 0.19, 0.20 with the opening degrees of 88.28, 84.48, 88.12%, and maximum peak von Mises stresses of 3.62, 1.21, 1.87 MPa, respectively. Conclusion. The developed algorithm makes it possible to automatically generate three-dimensional models of the leaflet apparatus, numerically simulate the opening process using the finite element method, statistically analyze the results obtained, and calculate the lumen area. The algorithm was verified based on the data for the Trifecta bioprosthesis of three standard sizes. The presented algorithm can be used both for the research and development of various designs and for obtaining “optimal” models of sash devices.

i Надоели баннеры? Вы всегда можете отключить рекламу.
iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.
i Надоели баннеры? Вы всегда можете отключить рекламу.

Текст научной работы на тему «An Algorithm for Automatic Generation and Evaluation of Leaflet Apparatus Models for Heart Valve Prostheses»

Generation and Evaluation lodels for Heart Valve Prostheses

henko, Junior Researcher, Laboratory of New Biomaterials, Department of Experimental Medicine; u. Klyshnikov, PhD, Researcher, Laboratory of New Biomaterials, Department of Experimental Medicine; Е.А. Ovcharenko, PhD, Head of the Laboratory of New Biomaterials, Department of Experimental Medicine;

L.S. Barbarash, MD, DSc, Professor, Academician of the Russian Academy of Sciences, Chief Researcher

Research Institute for Complex Issues of Cardiovascular Diseases, 6 Sosnovy Blvd, Kemerovo, 650002, Russia

The aim of the study is to develop and verify an algorithm for automatic generation of leaflet apparatus models for prosthetic heart valves, to optimize the basic parameters of the models in order to minimize the stress-strain state and maximize the geometric area of the orifice.

Materials and Methods. The suggested algorithm consists of three blocks: "Generator", "Modeling", "Analysis". The first block creates a three-dimensional model of the leaflet apparatus using the specified parameters (height, radius, thickness, degree of "sagging", angle of the free edge deviation). Numerical simulation of the apparatus functioning is further performed using the finite element method. Then, the statistical analysis of the von Mises stresses is done and the opening area of the design in question is calculated.

Verification was performed by comparing quantitatively the lumen areas of the leaflet apparatus in the open state, obtained from the literature data for the Trifecta bioprosthesis (19, 21, and 23 mm in diameter), with the results of the described algorithm operation.

Results. The verification of the algorithm has demonstrated the following deviations in the lumen area in the open state: 2.85% for 19 mm, 14.81% for 21 mm, and 23.17% for 23 mm models. This difference is due to the choice of the model material (no data could be found on the physical and mechanical properties of the pericardium used for the fabrication of the Trifecta bioprostheses).

The generation of a large number of designs (n=1517) without fixation of certain geometry parameters has shown that thickness of the leaflet apparatus makes the greatest contribution to the degree of opening; its dependence on the thickness and arising peak von Mises stresses has been demonstrated. Of the valvular models obtained, 278 showed the opening degree greater than 80% and maximum peak von Mises stresses below 4 MPa for the proposed model of the pericardium, which is 65% below the ultimate strength of the material.

Out of 278 leaflet models, 3 "optimal" designs were selected meeting the diameter criteria of 19, 21, and 23 mm. The loss index for them was 0.24, 0.19, 0.20 with the opening degrees of 88.28, 84.48, 88.12%, and maximum peak von Mises stresses of 3.62, 1.21, 1.87 MPa, respectively.

Conclusion. The developed algorithm makes it possible to automatically generate three-dimensional models of the leaflet apparatus, numerically simulate the opening process using the finite element method, statistically analyze the results obtained, and calculate the lumen area. The algorithm was verified based on the data for the Trifecta bioprosthesis of three standard sizes. The presented algorithm can be used both for the research and development of various designs and for obtaining "optimal" models of sash devices.

Key words: heart valve prosthesis; simulation; algorithm for generation and evaluation of leaflet apparatuses.

How to cite: Onischenko P.S., Klyshnikov K.Yu., Ovcharenko Е.А., Barbarash L.S. An algorithm for automatic generation and evaluation of leaflet apparatus models for heart valve prostheses. Sovremennye tehnologii v medicine 2022; 14(4): 6, https://doi.org/10.17691/ stm2022.14.4.01

This is an open access article under the CC BY 4.0 license (https://creativecommons.org/licenses/by/4.0/).

Corresponding author: Pavel S. Onischenko, e-mail: onis.pavel@gmail.com

//////////////////////^^^^

6 CTM I 2022 I vol. 14 j No.4 P.S. Onischenko, K.Yu. Klyshnikov, E.A. Ovcharenko, L.S. Barbarash

Introduction

High prevalence of cardiovascular diseases, acquired valvular defects in particular, results in the increase of operations for their replacement both in Russia (about 10 thousand operations annually [1]) and in the world (over 250 thousand [1, 2]). In contrast to mechanical prostheses, the biological ones do not require anticoagulant therapy, provide hemodynamics similar in character to the native, and can be implanted by a transcatheter technique. However, despite these benefits [3, 4], 50% of the implanted prosthetic heart valves develop structural changes in 10-15 years (leaflet rupture, calcinosis, pannus formation, leaflet avulsion in the commissural zone), which impair their functions [5, 6]. These alterations are caused by a prolonged effect of the immune system (immune rejection) on the implant [7] and accumulated mechanical fatigue damage in the leaflet apparatus [6, 8].

To solve the problem of the immune response, new materials are actively being developed [9, 10], methods of biomaterial conservation and treatment are being studied [11-13]. Engineers and researchers also pay attention to the second aspect of prosthesis dysfunction, fatigue damage, in order to create new types of heart valve prostheses without these drawbacks. However, their development is a great and time-consuming task at every step beginning from the design and material choice up to the grounding the safety and durability of the finished product [14]. It is unknown which of the characteristics of the implant makes the greatest contribution to the total strength of the medical device. The authors of some works demonstrate a significant effect of stress field distribution on the risks of damage to the prosthetic biomaterial, exposed to a long-term functioning, both in numerical [15, 16] and natural investigations [17] recommending optimization of a leaflet design in order to reduce the amplitude of stresses. Frequently, the current approaches to such optimization of the leaflet apparatus are of a subjective or exploratory character. Thus, a "classic" method of designing the leaflet apparatus for heart valve prostheses consists in the iterative selection of the optimal geometric parameters in the cycle: designing the "initial" geometry — hydrodynamic testing in vitro; correction of the leaflet geometry — retesting, and so on. The current tools of design and engineering analysis make it possible to simplify and reduce this succession by computer simulation [15, 18, 19] and exclusion of ineffective leaflet models. However, the first stage of forming the "initial" geometry remains unchanged and bears the burden of subjective conception of designers/ technologists, which may be originally incorrect.

A proposed option to reduce the probability of dysfunctions consists in optimization of the geometry of the most mobile and important element in the valve, its leaflet apparatus, which experiences the greatest loads and deformations during the cardiac cycle. In

recent time, the researchers in different countries have attempted to develop the leaflet apparatus [20-23] using a more systemic approach to design based on the mathematical methods and algorithms of geometry optimization excluding the above-mentioned subjectivity. This approach was motivated by a wider choice of the "initial" forms of the leaflet apparatus, introduction of the qualitative characteristics of their selection, search for design-effectiveness interconnections for a given component in the perspective of developing more perfect models of prosthetic heart valves. However, such works have some disadvantages: parametrized lists of linear dimensions create difficulties in some cases for understanding differences between the models when only one parameter is changed; researchers use maximum stress as the only parameter characterizing the functionality of a specific prosthesis design; the technologies of simulating the liquid-solid interaction are rather time-consuming (as they need computation of liquid movement apart from the deformation of a solid body) and cannot be operatively used in clinical practice.

The proposed algorithm for designing and evaluating the models of the leaflet apparatus for heart valve prostheses allows one to create final geometries with clear and comprehensible characteristics, automatically perform numerical simulation of the functioning process with the possibility of using certain material properties and setting pressures acting on the leaflet from the aorta and left ventricle. Linear dimensions may be obtained by processing patient's CT images and a pressure gradient from the results of Doppler echocardiography [24].

The aim of the study is to describe the developed algorithm for automatic generation of the leaflet apparatus models for heart valves, their numerical analysis and processing of the results obtained; to generate a large number of designs and choose the optimal ones for the prostheses with diameters of 19, 21, and 23 mm.

Materials and Methods

Algorithm for automatic generation of the leaflet apparatus

The algorithm for generation and investigation of the leaflet apparatus models proposed in the present work represents a succession of numerical methods including three blocks:

1. Generator. The task of this block is to create a three-dimensional facet leaflet model in STL format on the basis of the geometric input data. Owing to the automatic realization, it is possible to generate a large file of 3D models varying input data within the specified ranges.

2. Simulation. The work of every 3D leaflet model created by the generator is numerically simulated. This stage uses the final element method to analyze the stress-strain state in the engineering analysis environment.

3. Analysis. The quantitative characteristics of the numerical simulation are automatically analyzed: the stress-strain state and effectiveness of the work of every simulated leaflet.

A detailed description of each stage is presented below.

"Generator" block. The developed algorithm receives geometric parameters of the designed leaflet apparatus from the user in a simple and understandable form: height, size, degree of leaflet dome sagging, its thickness, and angle of the free edge deviation. On the basis of the point cloud generated by these parameters, building of the leaflet node links in a single closed surface is performed. Then, a 3D mesh is built and files are created, which will be used for initiating a numerical experiment (see "Simulation" block). The "Generator" is realized in the form of its own algorithm using the package of applied

programs Matlab R2021a (The MathWorks, USA) for solving the tasks of technical computations.

The following set of geometric parameters is used as input data (Figure 1 (a)):

H-i (mm) — height of the proposed leaflet from its lower part up to the top point of the commissural strut;

^ (°) — angle of the free leaflet edge deviation;

R (%) — degree of leaflet dome "sagging";

Rin (mm) — radius of the proposed prosthesis;

T (mm) — leaflet thickness;

Sec (°) — the number of degrees occupied by one leaflet.

"Simulation" block. The opening phase of the proposed prosthesis has been considered within the frames of this work, which allowed us to neglect the possible contacts and to simplify the model up to the investigation of one leaflet. The files prepared on

Figure 1. Materials and methods of investigation:

(a) designations of parameters for building the geometry and surface to which boundary conditions were applied; (b) pressure profiles where aortic pressure is applied to the F2 surface and ventricular to F1 the phase of the leaflet apparatus opening is marked by a rectangle; (c) results of testing the pericardium samples on the universal testing machine; the red line shows material properties selected for numerical simulation; (d) a graphic presentation of the deformed leaflet apparatus (left) and the top view (right); grey line shows the leaflet contour in the initial position

//////////////////////^^^^

8 CTM j 2022 j vol. 14 j No.4 P.S. Onischenko, K.Yu. Klyshnikov, E.A. Ovcharenko, L.S. Barbarash

the generator were transmitted to the engineering analysis environment Abaqus/CAE (Dassault Systemes, France) in which the model leaflet opening process was automatically numerically analyzed.

Boundary conditions for numerical simulation were set for the surfaces (see Figure 1 (a)) within the time segment of 0-0.2 s (Figure 1 (b)), which corresponds to the phase of leaflet apparatus opening:

F-i — pressure from the left ventricle (Figure 1 (a), ventricular pressure);

F2 — pressure from the aorta (Figure 1 (a), aortic pressure);

F3 — part of the leaflet apparatus fixed from the motion in 3 degrees of freedom in cylindrical coordinates.

The model of the material for numerical simulation was obtained by testing the samples (n=10) of xenopericardium, from which a leaflet apparatus of clinical biological prostheses was fabricated, on the universal Z-series testing machine (Zwick/Roel, Germany) using the sensor with a nominal load of 50 N at the ambient temperature of 37°C maintained by a thermal chamber. Samples were cut on the cutting press using a knife of a specific shape (ISO 37, type 4). Sample thickness was measured with a TRP thickness gauge with ±0.01-mm margin of permissible error (pressing force below 1.5 N). The speed of the cross-arm movement during testing was 20 mm/min. The obtained data were exported as the ratio of the relative elongation (mm) to stress (MPa) endured by the sample.

Taking into consideration the variability of the biomaterial mechanical properties determined by its natural non-uniformity, the whole set of separate curves was averaged using a third-order polynomial function built by the method of normal equations. In this way, a stress-strain conditional average curve was obtained for this material, which was employed at the stage of numerical simulation. Goodness of fit degree between the approximating curve and initial data of the uniaxial tensile was estimated by the coefficient of determination R2, which for the third-degree polynomial was equal to 0.86, which may be considered satisfactory for the tasks of the present study. The resulted properties of the material being used further for numerical simulation are shown with a red line in Figure 1 (c).

"Analysis" block. Using the results of numerical simulation, the mean value, maximum von Mises stress, and standard deviation were calculated. The degree of opening (DO) of the leaflet apparatus was computed as a ratio of a leaflet area projection on the XY surface to the initial non-deformed state (Figure 1 (d)). Since the considered model contains one leaflet, the area of its projection in the non-deformed state will occupy 1/3 of the circumference area with the Rin radius. Consequently, the desired DO of the leaflet apparatus may be determined by the formula:

where Scirc is the area of the circumference with the Rin radius; S is the area of opening during simulation.

To evaluate integrally the designs under consideration, a loss index (LI) was introduced, the formalized form of which may be written down in the following way:

LI =V (100%-DO)2 + (SNm-Smax)2 ,

where Smax — maximum stress arising in the process of simulation; S|im — ultimate strength of the pericardium [25]. A simple interpretation may be presented by "the less, the better" expression. At D0=100% and zero maximum stresses, the LI will be equal to 0, in other cases it will tend to 1. At the increased values of Slim, the proposed method of assessment will show the value exceeding 1, which also demonstrates a failure of the considered leaflet apparatus design.

Algorithm verification

Despite the application of the validated engineering tools, i.e. building faceted bodies and numerical simulation on the Abaqus/CAE platform, the generation-simulation-analysis integral link requires quality and consistency testing. In this connection, generation of three leaflet geometries of the Trifecta commercial clinical prosthesis, 19, 21, and 23 mm in size (St. Jude Medical, Inc., USA) on the basis of changing the geometric characteristics of highly precise 3D models using the proposed algorithm became a separate block for investigation (Figure 2). These models were obtained by microtomographic scanning using the OREL-MT

Prosthesis model

Parameter 19 mm 21 mm 23 mm

9 (°) 0 0 0

H (mm) 11.0 11.5 12.0

R (%) 75 75 75

Rin (mm) 9.5 10.5 11.5

Sec (°) 119 119 119

T (mm) 0.5 0.5 0.5

DO =

S

W 3

■100%,

Figure 2. Three-dimensional reconstructions of clinical heart valve prostheses (three standard sizes of Trifecta prosthesis) obtained on the basis of microtomographic analysis, and quantitative measurements of their geometric parameters

experimental system (Tomsk, Russia). Using the obtained two-dimensional DICOM slices, facet bodies were reconstructed and their parameters measured (see the table in Figure 2).

These measurements of the 3D models were used as input data for the generation of three leaflets by means of the developed algorithm for each of the presented Trifecta bioprostheses at the stage of validation. The literature data were used as a material model for this case [26], since Trifecta bioprostheses are fabricated from the bovine xenopericardium stabilized with glutaraldehyde. The area of leaflet opening obtained during simulation by the proprietary algorithm was compared with the data from the experimental studies reported in the literature [27-29] and manufacturer's documentation for this bioprosthesis [30].

Investigation of a large set of leaflet geometries

The second step of our work was an attempt to assess the leaflet apparatus varying the parameters, the range of which was chosen to cover the majority of the existing tricuspid valve prostheses in aortic and tricuspid position. 1517 unique geometries were generated by our algorithm. They became the basis of numerical simulation. The following boundaries for changing the leaflet apparatus parameters were set:

H1=[10; 25] mm;

q>=[-30; 30]°;

R=[0; 100]%;

Rin=[7.5; 20.0] mm;

T=[0.1; 1.0] mm.

Results

Verification. In the process of verification, DO of the built algorithms for leaflet apparatus models were evaluated and the results obtained were compared with those of the clinical studies [27-29] and with the manufacturer's documentation [30] (Figure 3).

It has been established that the algorithm is capable of reproducing the process of opening the 19-mm Trifecta prosthesis with the error of 2.85%, 21-mm prosthesis —

260 -240 -

I 220 -

03

£ 200 -ro

. E 180-

140 -120 -

19 21 23

Standard size

Figure 3. Comparison of the values of opening areas for the leaflet apparatus models obtained in the process of numerical simulation using the proposed algorithm (black line) with the data derived from the literature sources [27-30]

The results are presented as median, 25 and 75 percentiles, minimum and maximum for each standard size

250

700

57

Degree of opening (%)

93

30

^ ,

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

ofoS

Maximum stress (MPa)

120

fa

Degree of opening (%)

^ ^ ^ ^ ^ Maximum stress (MPa)

Figure 4. Quantitative distribution of models by the "degree of opening" and "maximum stress" parameters:

Q3 and Q4 — from above, Q4 — from below

////////////////////^^^^

10 CTM J 2022 J vol. 14 J No.4 P.S. Onischenko, K.Yu. Klyshnikov, E.A. Ovcharenko, L.S. Barbarash

14.81%, and 23-mm — 23.17% according to the "area of leaflet apparatus opening" parameter. These differences are caused by the model of the material, which is "softer" than the one used in the fabrication of the reference prosthesis [26].

Investigation of a large set of leaflet geometries. Analyzing the results of 1517 models built, a significant scattering of the DO parameter (71.97±15.91%) has been noted. So, for 759 models (Q3 and Q4, DO>75.19%) the mean DO value was 82.68±4.8% and maximum arising stress — 6.25±7.28 MPa. In the "best" 380 leaflet apparatuses (Q4, DO>82.14%), this value was 86.76±2.96% with maximum stress of 6.71 ±8.08 MPa. The histograms of distribution by the degree of opening and maximum von Mises stress are shown in Figure 4.

Despite the fact that an essential part of the geometries shows the maximum von Mises stress less than 11.6 MPa (this value is a destruction threshold for the pericardium used in bioprostheses [25]), some models exceed this threshold. It speaks of the initially non-optimal parameters used in designing the leaflet apparatus, which will most likely result in dysfunction of the heart valve prosthesis.

A reasonable question is: which of the obtained models gives considerable DO (more than 80% in our case) and at the same time have the peak values below the destruction threshold. To find the solution, the results obtained have been sorted using the suggested conditions. As a result, 403 models were derived with the mean DO values (85.63±3.43%) and maximum stresses (3.81±2.72 MPa). The graphs of dependences of maximum stress or DO on geometric parameters (apart from a leaflet thickness, Figure 5) do not show any evident regularity and the trend line represents a horizontal line, that is why they are not shown in Figure 5. The leaflet thickness has influenced the DO considerably. This effect may be associated with physical and mechanical properties of the pericardium used for the material models.

It should be noted that 69% of all 403 selected leaflet apparatus models have peak values below 4 MPa, the average LI=0.25 (minimal — 0.16, maximal — 0.36), which gives space for the selection of "optimal" geometries. The data presented in Figure 5 allow us to suppose that:

1) the higher the stress, the greater the degree of opening;

2) the thicker the leaflet, the lower the degree of opening and, consequently, the smaller maximum arising stresses in the leaflet apparatus.

"Optimal" geometry. The main applied motivation for the creation of this algorithm is an automatic selection of such geometric parameters of the leaflet apparatus which provide the greatest performance — the area of opening at minimally possible stress amplitudes. Of the large set of the leaflet geometry designs selected in the process of investigation (see markers in Figure 5), 3 models (see the Table) for the 19-, 21-, and 23-mm prostheses are designated as conventionally "optimal".

<u Q

96 94 92 90 88 86 84 82 80 78

96 94 92 90 88

86 °-0t

84 o " " ,00 „„ o "o o 82 oo

f. /%.• 0

o8°

o o

o O o o

o o o o o o 8 °

- °cp o° ...fi______""—bo9""'"

% '

o <b

.....

..9

o

o« 8

o

° <9 o

8 10 12 Maximum stress(MPa)

80 78

0.1

0.2

0.3

0.4

0.5

0.6

0.7

0.8 0.9 1 Thickness(mm)

1.0 0.9 0.8 0.7 0.6 0.5 0.4 0.3 0.2 0.1 0

10 12 Maximumstress(MPa)

Figure 5. Distribution of geometries after the selection by the criteria "degree of opening" more than 80% and "maximum stress" less than 11.6 MPa

Dash line — exponential trend line. Rhomb, square, and triangle mark the location of leaflet apparatus designs with 19-, 21-, and 23-mm standard sizes, respectively (see "Optimal" geometry")

2

4

6

2

4

6

8

Characteristics of "optimal" designs of leaflet apparatuses

Parameter Prosthesis model

19 mm 21 mm 23 mm

9 (°) 5.17 -26.16 -7.94

H-] (mm) 15.22 18.81 22.52

Rin (mm) 9.38 10.49 11.56

Sec (°) 119 119 119

T (mm) 0.14 0.27 0.19

Degree of opening (%) 88.28 84.48 88.12

Smax (MPa) 3.62 1.21 1.87

Loss index 0.24 0.19 0.2

19

21

23

UUU

Figure 6. "Optimal" designs of leaflet apparatuses with 19-, 21-, and 23-mm standard sizes in the open state:

(a) top view, (b) side view, (c) isometric view

Figure 6 shows "optimal" leaflet apparatuses in the closed and open state. Stresses are localized in the commissural zones, but numerical values do not exceed critical for the biomaterial model used. Qualitatively, it is seen that these models provide complete and symmetric opening with maximization of the area.

Discussion

The algorithm of assessing the effect of geometric parameters of the heart valve prostheses on the peak stress and DO values, the most critical in selecting the

optimal leaflet shapes, is presented in this work. Verification of the algorithm operation was exemplified by the opening area of the leaflet apparatus of 19-, 21-, and 23-mm Trifecta prosthesis. This step has shown the convergence of the proposed algorithm with clinical data [27-29] and manufacturer's documentation [30].

A large number of leaflet apparatus models have been generated using various parameters. The results obtained show that in the models with greater lumen area in the open state, the arising stress increases and, in some combinations of the geometric parameters, the final design may be initially "non-viable" due to critical stresses in the valve leaflet material. However, using this algorithm at the stage of prosthesis design, it is possible to find such combinations of geometric parameters that the final item would be optimal by the two main criteria of durability in terms of biomechanics and functioning: maximal DO and minimal value of von Mises stresses. An additional benefit of this method is gaining a large bulk of information about the biomechanics of a specific design depending on the material used.

The results of this work agree with the data of other studies: a small change in one of the parameters can be critically reflected on the arising peak stresses and DO. The material of the leaflet apparatus and thickness of the leaflet itself play an important role in the quality of the prosthetic valve operation [20, 31].

Application prospects. In the recent decade, a lot of works appeared on the investigation of polymer prosthetic heart valves [9, 32-37]. The undisputable advantage of this material is the possibility of fabricating devices of a complicated shape, which is almost impossible to achieve for their biological analogs due to cutting of the latter from the sheet biomaterial. The presented algorithm of automatic generation and investigation of the obtained leaflet apparatus models may be applied for creation of optimal designs of polymer leaflets from the standpoint of biomechanics. Since the polymer-based heart valve prostheses are not limited by cutting from the plain sheet, this approach opens new possibilities for engineers: manufacturing the products by molding, immersion (dipping) [32, 33], or 3D printing [35]. These advantages over bioprostheses make it possible to use more complex and effective designs of the leaflet apparatus, which may be selected using the presented algorithm.

a

b

c

////////////////////^^^^

12 CTM 1 2022 1 vol. 14 J No.4 P.S. Onischenko, K.Yu. Klyshnikov, E.A. Ovcharenko, L.S. Barbarash

Limitations. Limitations of the described algorithm are worth mentioning. As the algorithm is a primary approximation to the complex development of the tools for engineer support in the field of designing heart valve prostheses, we used the model of one leaflet, which does not take into consideration the contacts with two adjacent leaflets in the diastolic phase. Further improvement of the algorithm will solve the task adding to the obtained results the area of cooptation and contact stress. The base of geometric and mechanical characteristics of leaflet designs may be employed as a starting point in creating the most suitable leaflet apparatuses for specific prosthetic models using a machine learning technique — the next step towards the development of the optimal geometry.

Conclusion

The algorithm for automatic generation and evaluation of leaflet apparatuses for heart valve prostheses has been developed and verified. Its verification was performed on the basis of Trifecta bioprosthesis of three standard sizes. Of the 1517 models obtained, 278 showed the opening degree greater than 80% and maximum peak von Mises stresses below 4 MPa for the proposed model of the pericardium, which is 65% below the ultimate strength of the material. The presented designs show the capability of the algorithm to generate "optimal" leaflet apparatuses if certain known geometric parameters are set.

Study funding. The work was supported by the Russian Science Foundation, grant No.21-75-10128, https://rscf.ru/project/21-75-10128/.

Conflicts of interest. The authors have no conflicts of interest to declare.

References

1. Bokeriya L.A., Milievskaya E.B., Kudzoeva Z.F., Pryanishnikov V.V., Skopin A.I., Yurlov I.A. Serdechno-sosudistaya khirurgiya — 2018. Bolezni i vrozhdennye anomalii sistemy krovoobrashcheniya [Cardiovascular surgery — 2018. Diseases and congenital anomalies of the circulatory system]. Moscow: Nauchnyy tsentr serdechno-sosudistoy khirurgii RAMN; 2019; 270 p.

2. Siddiqui R.F., Abraham J.R., Butany J. Bioprosthetic heart valves: modes of failure. Histopathology 2009; 55(2): 135-144, https://doi.org/10.1111/j.1365-2559.2008.03190.x.

3. Siregar S., de Heer F., Groenwold R.H.H., Versteegh M.I.M., Bekkers J.A., Brinkman E.S., Bots M.L., van der Graaf Y., van Herwerden L.A. Trends and outcomes of valve surgery: 16-year results of Netherlands Cardiac Surgery National Database. Eur J Cardiothorac Surg 2014; 46(3): 386397, https://doi.org/10.1093/ejcts/ezu017.

4. Isaacs A.J., Shuhaiber J., Salemi A., Isom O.W., Sedrakyan A. National trends in utilization and in-hospital outcomes of mechanical versus bioprosthetic aortic valve replacements. J Thorac Cardiovasc Surg 2015; 149(5): 1262-1269.e3, https://doi.org/10.1016/j.jtcvs.2015.01.052.

5. Capodanno D., Petronio A.S., Prendergast B., Eltchaninoff H., Vahanian A., Modine T., Lancellotti P., Sondergaard L., Ludman P.F., Tamburino C., Piazza N., Hancock J., Mehilli J., Byrne R.A., Baumbach A., Kappetein A.P., Windecker S., Bax J., Haude M. Standardized definitions of structural deterioration and valve failure in assessing long-term durability of transcatheter and surgical aortic bioprosthetic valves: a consensus statement from the European Association of Percutaneous Cardiovascular Interventions (EAPCI) endorsed by the European Society of Cardiology (ESC) and the European Association for Cardio-Thoracic Surgery (EACTS). Eur J Cardiothorac Surg 2017; 52(3): 408-417, https://doi.org/10.1093/ejcts/ezx244.

6. Rodriguez-Gabella T., Voisine P., Puri R., Pibarot P., Rodes-Cabau J. Aortic bioprosthetic valve durability: incidence, mechanisms, predictors, and management of surgical and transcatheter valve degeneration. J Am Coll Cardiol 2017; 70(8): 1013-1028, https://doi.org/10.1016/j. jacc.2017.07.715.

7. Kostyunin A.E., Yuzhalin A.E., Rezvova M.A., Ovcharenko E.A., Glushkova T.V., Kutikhin A.G. Degeneration of bioprosthetic heart valves: update 2020. J Am Heart Assoc 2020; 9(19): e018506, https://doi.org/10.1161/jaha.120.018506.

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

8. Abbasi M., Barakat M., Dvir D., Azadani A. Detailed stress analysis of Edwards-SAPIEN and Medtronic CoreValve devices. Is leaflet stress comparable to surgical Carpentier-Edwards PERIMOUNT Magna bioprosthesis? Struct Heart 2019; 3(Suppl 1): 192, https://doi.org/10.1080/24748706.2019. 1591103.

9. Oveissi F., Naficy S., Lee A., Winlaw D.S., Dehghani F. Materials and manufacturing perspectives in engineering heart valves: a review. Mater Today Bio 2019; 5: 100038, https://doi. org/10.1016/j.mtbio.2019.100038.

10. Bezuidenhout D., Williams D.F., Zilla P. Polymeric heart valves for surgical implantation, catheter-based technologies and heart assist devices. Biomaterials 2015; 36: 6-25, https:// doi.org/10.1016/j.biomaterials.2014.09.013.

11. Kudryavtseva Yu.A. Bioprosthetic heart valves. From idea to clinical use. Kompleksnye problemy serdecno-sosudistyh zabolevanij 2015; 4: 6-16.

12. Zhuravleva I.Y., Karpova E.V., Oparina L.A., Cabos N., Ksenofontov A.L., Zhuravleva A.S., Nichay N.R., Bogachev-Prokophiev A.V., Trofimov B.A., Karaskov A.M. Bioprosthetic xenopericardium preserved with di- and penta-epoxy compounds: molecular cross-linking mechanisms, surface features and mechanical properties. Patologiya krovoobrashcheniya i kardiokhirurgiya 2018; 22(3): 49-55, https://doi.org/10.21688/1681-3472-2018-3-56-68.

13. Kostyunin A.E., Glushkova T.V. Expression of matrix metalloproteinases 1, 2, 9, 12 in xenogenic tissues of epoxy-crosslinked bioprosthetic heart valves explanted due to dysfunction. Rossijskij kardiologiceskij zurnal 2020; 25(10): 3978, https://doi.org/10.15829/1560-4071-2020-3978.

14. Rotman O.M., Bianchi M., Ghosh R.P., Kovarovic B., Bluestein D. Principles of TAVR valve design, modelling, and testing. Expert Rev Med Devices 2018; 15(11): 771-791, https://doi.org/10.1080/17434440.2018.1536427.

15. Martin C., Sun W. Simulation of long-term fatigue damage in bioprosthetic heart valves: effects of leaflet and stent elastic properties. Biomech Model Mechanobiol 2014; 13(4): 759-770, https://doi.org/10.1007/s10237-013-0532-x.

16. Martin C., Sun W. Comparison of transcatheter aortic valve and surgical bioprosthetic valve durability: a fatigue

simulation study. J Biomech 2015; 48(12): 3026-3034, https:// doi.org/10.1016/j.jbiomech.2015.07.031.

17. Zhou H., Wu L., Wu Q. Structural stability of novel composite heart valve prostheses — fatigue and wear performance. Biomed Pharmacother 2021; 136: 111288, https://doi.org/10.1016/j.biopha.2021.111288.

18. Mohammadi H., Mequanint K. Prosthetic aortic heart valves: modeling and design. Med Eng Phys 2011; 33(2): 131147, https://doi.org/10.1016/j.medengphy.2010.09.017.

19. Praveen Kumar G., Mathew L. Three-dimensional computer-aided design-based geometric modeling of a new trileaflet aortic valve. Artif Organs 2010; 34(12): 1121-1124, https://doi.org/10.1111/j.1525-1594.2009.00973.x.

20. Li K., Sun W. Simulated transcatheter aortic valve deformation: a parametric study on the impact of leaflet geometry on valve peak stress. Int J Numer Method Biomed Eng 2017; 33(3): e02814, https://doi.org/10.1002/cnm.2814.

21. Hsu M.C., Kamensky D., Xu F., Kiendl J., Wang C., Wu M.C.H., Mineroff J., Reali A., Bazilevs Y., Sacks M.S. Dynamic and fluid-structure interaction simulations of bioprosthetic heart valves using parametric design with T-splines and Fung-type material models. Comput Mech 2015; 55(6): 1211-1225, https://doi.org/10.1007/s00466-015-1166-x.

22. Xu F., Morganti S., Zakerzadeh R., Kamensky D., Auricchio F., Reali A., Hughes T.J.R., Sacks M.S., Hsu M.C. A framework for designing patient-specific bioprosthetic heart valves using immersogeometric fluid-structure interaction analysis. Int J Numer Method Biomed Eng 2018; 34(4): e2938, https://doi.org/10.1002/cnm.2938.

23. Travaglino S., Murdock K., Tran A., Martin C., Liang L., Wang Y., Sun W. Computational optimization study of transcatheter aortic valve leaflet design using porcine and bovine leaflets. J Biomech Eng 2020; 142(1): 011007, https:// doi.org/10.1115/1.4044244.

24. Bokeriya L.A., Belal G. Methods for echocardiographic evaluation of aortic valve hemodynamics after prosthetics: methods and precautions. Kreativnaa kardiologia 2012; 6(1): 73-79.

25. Ovcharenko E.A., Klyshnikov K.U., Yuzhalin A.E., Savrasov G.V., Glushkova T.V., Vasukov G.U., Nushtaev D.V., Kudryavtseva Y.A., Barbarash L.S. Comparison of xenopericardial patches of different origin and type of fixation implemented for TAVI. Int J Biomed Eng Technol 2017; 25(1):

59, https://doi.org/10.1504/ijbet.2017.10007484.

ndo-Espitia N., Sanchez-Arevalo F., Hernandez-hanical assessment of bovine pericardium imaging, enhanced backscattering and analysis. Biomed Opt Express 2015; org/10.1364/boe.6.002953.

Bichi S., Biondi A., Blasio A., B., Gherli R., Leva C., Alfieri O., Antona C., A., Martinelli L.,

Menicanti L., Paolini G., Beghi C. St. Jude Medical Trifecta aortic valve: results from a prospective regional multicentre registry. J Cardiothorac Surg 2015; 10: 169, https://doi. org/10.1186/s13019-015-0379-6.

28. Permanyer E., Estigarribia A.J., Ysasi A., Herrero E., Semper O., Llorens R.St. Jude Medical Trifecta™ aortic valve perioperative performance in 200 patients. Interact Cardiovasc Thorac Surg 2013; 17(4): 669-672, https://doi.org/10.1093/ icvts/ivt270.

29. Mannam G., Mishra Y., Modi R., Gokhale A.G.K., Sethuratnam R., Pandey K., Malhotra R., Anand S., Borah A., Mukhopadhyay S., Shah D., Mahant T.S. Early hemodynamic performance of the Trifecta™ surgical bioprosthesis aortic valve in Indian patient population: 12 month outcomes of the EVEREST post-market study. J Cardiothorac Surg 2018; 13(1): 96, https://doi.org/10.1186/s13019-018-0783-9.

30. Trifecta™ Valve with Glide™ Technology (GT) instructions for use. St. Jude Medical; 2016; URL: https:// manuals.sjm.com/~/media/manuals/product-manual-pdfs/e/1/ e1157bdf-2fb8-4567-81ef-10f620e79af9.pdf.

31. Abbasi M., Barakat M.S., Dvir D., Azadani A.N. A noninvasive material characterization framework for bioprosthetic heart valves. Ann Biomed Eng 2019; 47(1): 97-112, https://doi. org/10.1007/s10439-018-02129-5.

32. Rahmani B., Ghanbari H., Tzamtzis S., Burriesci G., Seifalian A.M. Polymeric heart valves. In: Encyclopedia of biophysics. Roberts G.C.K. (editor). Springer-Verlag Berlin Heidelberg; 2013; p. 1893-1900, https://doi.org/10.1007/978-3-642-16712-6_702.

33. Powell S.K., Cruz R.L.J., Ross M.T., Woodruff M.A. Past, present, and future of soft-tissue prosthetics: advanced polymers and advanced manufacturing. Adv Mater 2020; 32(42): e2001122, https://doi.org/10.1002/adma.202001122.

34. Motta S.E., Falk V., Hoerstrup S.P., Emmert M.Y. Polymeric valves appearing on the transcatheter horizon. Eur J Cardiothorac Surg 2021; 59(5): 1057-1058, https://doi. org/10.1093/ejcts/ezab089.

35. Coulter F.B., Schaffner M., Faber J.A., Rafsanjani A., Smith R., Appa H., Zilla P., Bezuidenhout D., Studart A.R. Bioinspired heart valve prosthesis made by silicone additive manufacturing. Matter 2019; 1(1): 266-279, https://doi. org/10.1016/j.matt.2019.05.013.

36. Ghanbari H., Viatge H., Kidane A.G., Burriesci G., Tavakoli M., Seifalian A.M. Polymeric heart valves: new materials, emerging hopes. Trends Biotechnol 2009; 27(6): 359-367, https://doi.org/10.1016/jtibtech.2009.03.002.

37. Stasiak J.R., Serrani M., Biral E., Taylor J.V., Zaman A.G., Jones S., Ness T., De Gaetano F., Costantino M.L., Bruno V.D., Suleiman S., Ascione R., Moggridge G.D. Design, development, testing at ISO standards and: in vivo feasibility study of a novel polymeric heart valve prosthesis. Biomater Sci 2020; 8(16): 4467-4480, https://doi. org/10.1039/d0bm00412j.

/my/mmmy/mmmmmmy/mmmy/mmmy/mmmmmmi

i Надоели баннеры? Вы всегда можете отключить рекламу.