|
|
||||||||
Ann Thorac Surg 2002;73:1122-1129
© 2002 The Society of Thoracic Surgeons
a Department of Mechanical Engineering, Indian Institute of Technology, Madras, Chennai, India
b Department of Cardiothoracic Surgery, Ramachandra Medical College, Porur, Chennai, India
Accepted for publication October 17, 2001.
* Address reprint requests to Dr Balakrishnan, Professor and Head, Department of Cardiothoracic Surgery, Ramachandra Medical College, Porur, Chennai 600 116, India
e-mail: cvskrb{at}giasmd01.vsnl.net.in
| Abstract |
|---|
|
|
|---|
Methods. Dynamic finite element analysis was used along with a geometrically accurate model of the aortic valve and the sinuses. Shell elements along with proper contact conditions were also used in the model. Pressure patterns during the cardiac cycle were given as an input, and a linear elastic model was assumed for the material.
Results. We found that aortic root dilation starts before the opening of the leaflet and is substantial by the time leaflet opens. Dilation of the root alone helps in opening the leaflet to about 20%. The equivalent stress pattern shows an instantaneous increase in stress at the coaptation surface during closure. Stresses increase as the point of attachment is approached from the free surface.
Conclusions. The complex interplay of the geometry of the valve system can be effectively analyzed using a sophisticated dynamic finite element model. Results not previously brought out by the earlier static analysis shed new light on the root/valve interaction.
| Introduction |
|---|
|
|
|---|
Sophisticated analysis of the complex interplay between the aortic root and the leaflets, as well as the development of stresses and strains in the leaflet during opening and closing, is now possible using powerful computer techniques. This may aid in a better understanding of the dynamics of the aortic valve complex and hopefully will help in designing better substitutes for replacing the diseased aortic valve. More importantly, these might prove invaluable in predicting patterns of failure of replacement devices or repair techniques in the workbench setting rather than experiencing it in the clinical situation, with unfortunate consequences.
There have been several attempts to understand the mechanics and function of the aortic valve and its surrounding geometric structures through finite element analysis, a powerful mathematical tool [13]. These models have used static finite element analysis. Dynamic analysis has been used to study stresses in an artificial leaflet by Thornton [4]; but only a small part of the cardiac cycle was used for the study, which does not bring out the dynamic aspects of the leaflet/aortic root interplay. As pointed out by Thubrikar in the discussion of work by David and colleagues [5], there is a need to understand the geometry of the valve in a dynamic state. Also, no other author has studied the aortic valve using a geometrically accurate system of leaflet and the sinuses with a dynamic finite element model. The work undertaken in this study brings out certain important features that are not captured by the static finite element model.
| Material and methods |
|---|
|
|
|---|
|
|
One of the major issues in modeling an aortic valve (or any other biological organ) is the definition of the stress-free state. As pointed out by Xie and colleagues [7], an unloaded organ (in their study, a vein) is not in a stress-free state. Clearly, further work is required to study the stress-free state or to quantify the residual stress state of an aortic leaflet. The major effect of the residual stress state is the actual material behavior, which is highly nonlinear. In this work, the leaflets were assumed to be stress free in the open position.
Although the finite element model is complete and detailed, the material used in this study is linear elastic. The actual material model ideally to be used is hyperelastic, taking into account the anisotropy. Though ABAQUS supports hyperelastic material model, the published data [1], essentially uniaxial in nature is not accurate enough to fit an Ogden model or Mooney-Rivlin model [6]. A proper material characterization requires biaxial tension results as well. This aspect may be the main limitation of this work. However, in this study, the modulus has been chosen to mimic the actual behavior as closely as possible.
The elastic modulus in this work has been taken to be 1 Mpa for the leaflet [1]. The stress-strain relation is essentially bilinear in nature, with a transition region [1, 8]. The modulus varies by a factor of nearly 150 between the pre- and posttransition strains [1]. The range of the modulus is also very high even within a region. Grande and colleagues [8] have assumed the leaflets to be in the posttransition region and have used a corresponding modulus value. However, within a leaflet there are areas where the strains are in the transition region or even in the pretransition region. Hence it is difficult to choose a value of modulus for the analysis, and from the data given in [1], 1 MPa appears reasonable. The modulus of the sinus was taken to be 2 MPa. Experiments conducted on porcine aortic root wall [9] show that the behavior of the aortic root is highly anisotropic. Youngs modulus values measured are about 1 to 2 MPa. Hence, in this work 2 MPa has been assumed.
Such assumptions will not have a serious effect on stress, as the deformation is controlled by pressure acting on the aortic system. However, the calculated displacements will be conservative for a higher modulus. A Poissons ratio of 0.3 was chosen to overcome the numerical difficulty encountered at values near 0.5. The density of the leaflets was taken to be 1.1 g/mL and that of the sinus 2.0 g/mL.
Finite element analysis requires restraint conditions and pressure to be given as inputs for the problem. The following pressure-time relationship is assumed in this work [10].
The surfaces were classified as aortic surface, ventricular surface, and leaflet surface. The pressure on the aortic surface was modeled as two ramps, with the first ramp accounting for a pressure drop from 120 mm Hg to 80 mm Hg and the second ramp indicating an increase of pressure from 80 to 120 mm Hg. The complete cycle is shown in Figure 2. The pressure applied to the ventricular region is also shown. Again, the pressure variations are modeled as a number of linear ramp loadings. The pressure applied in the leaflet is, in principle, the pressure difference between the ventricle and the aorta. Figure 2 also brings out this pressure difference. The pressure cycle followed in this work is explained with respect to Figure 2. The start of the analysis, T0 indicates a particular stage in the diastolic cycle. T1 indicates the end of the diastolic cycle. The pressure in the left ventricle increases from this time. At T2, which equals approximately 0.4 seconds, the leaflet pressure attains a positive value and the leaflet opens. At T3 (approximately 0.47 seconds), the aortic pressure overcomes the ventricular pressure and the leaflets start closing. Although the entire cardiac cycle is for 0.83 seconds, in the actual execution of this program, the time of closure is reduced to compress the time required to run the problem.
|
t is given by the Courant conditions as:
t
Lmin /Cd, where Lmin is the smallest element dimension in the mesh and Cd is the dilatational wave speed in the material. Because extremely fine mesh has been used in this study, the stable step size is very small. Nevertheless, no artificial increase in density has been resorted to increase
t. The program was run in a twin processor SGI/ORIGIN 200 workstation (Silicon Graphics, Inc., Mountain View, CA, USA). Typical processor time is about 20 Cpu hours.
| Results |
|---|
|
|
|---|
Aortic root and leaflet deformation
The deformation of the leaflet and aortic wall present an interesting spectacle. Figure 3
provides a guide for understanding the pattern of deformation.
|
|
|
|
|
|
|
|
|
| Comment |
|---|
|
|
|---|
Grande-Allen and colleagues [2] used a finite element model to study the effect of sinuses on the stresses developing in leaflets in valve-sparing root replacements. There have been other studies to elucidate various aspects of the functioning of aortic valve [8, 12]. Cacciola and colleagues [3] used similar techniques to show that a stentless trileaflet valve has 75% lower stress levels in the leaflet as compared with stented valves. However, all of these studies have been carried out using static finite element model. We used a dynamic finite element model of the aortic valve, and this has brought out certain features not previously captured, especially with regard to the complex interplay between the aortic root and the leaflets. We studied two main aspects of aortic valve function in a typical cycle, namely, aortic root and leaflet deformation; and stresses in the various parts of the leaflet.
Aortic root and leaflet deformation
The most striking finding was that the leaflets begin to open even before any positive pressure is applied, due primarily to the effect of aortic root dilatation. In fact, substantial increase in diameter during systole occurs before the actual opening of the aortic valve. This finding is striking and was confirmed by a separate pseudo-analysis in which no pressure was applied to the leaflets but only to the aortic root. This type of parametric testing, which is easily possible, is indeed one of the strengths of this technology. Interestingly, this finding of the effect of aortic root dilatation in aiding of the leaflet opening, has been confirmed independently using animal models and sonomicrometry [5, 13]. Although the complex interplay between the aortic leaflets and the aortic root has long been obvious, this is the first time that this relationship has been strikingly demonstrated in a computer simulation model. In fact, when the compliance of the wall is increased during the parametric study, the leaflets begin to prolapsean event that is known to happen in patients with annuloaortic ectasia causing aortic incompetence. It is thus not altogether unexpected to encounter instances of leaflet thickening following the David technique of including the aortic valve in a rigid cylinder [2]. The other striking finding is the violent shudder in the leaflets during opening and closing. It is therefore not surprising to find leaflet substitutes made of disparate material failing over a period of time, inside the aortic root, when used as replacement devices.
von-Mises stress in the leaflets
The other aspect that was studied was the development of stress in various parts of the leaflets. The important observation is that there is an instantaneous increase in the stress at the coaptation surface during closure. Also, there is an increase in stress as the point of attachment of the leaflet to the root is approached from the free surface (Fig 11). Although it is not possible to make a direct comparison with the work of Grande-Allen and colleagues [2] and Grande and associates [12], the stress states are in the same range as reported in those works.
Conclusions
With increasing sophistication in the surgical treatment of aortic root diseases, including in the application of valve-sparing root replacements and aortic valve repairs, there is clearly a need for a precise understanding of the complex interplay of the various components of the aortic root. Historically, surgical treatment of valvular heart disease has been based on designing a strategy or a valve substitute, applying it in a clinical setting, and then waiting (sometimes for years) to see if it will work. Some of the earlier failures of replacement devices could have been predicted if sophisticated techniques to test them had been available. With the widespread availability of powerful computers and supportive software, precise mathematical modeling using techniques such as finite element analysis make it possible to model and study the aortic valve in a dynamic state. These techniques also allow one to observe the most minute changes taking place in the leaflets and root, as well as the stresses developing in the leaflets in a manner that is not possible using animal or other models. Use of these techniques will enable the testing of newer devices, development of new leaflet materials, as well as newer types of stents for stent-mounted valves and repair techniques at the computer workbench before clinical application, so that better long-term results can be expected.
| Acknowledgments |
|---|
|
|
|---|
| References |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
L. Weltert, M. D. de Tullio, L. Afferrante, A. Salica, R. Scaffa, D. Maselli, R. Verzicco, and R. De Paulis Annular dilatation and loss of sino-tubular junction in aneurysmatic aorta: implications on leaflet quality at the time of surgery. A finite element study Interact CardioVasc Thorac Surg, March 27, 2013; (2013) ivt116v1. [Abstract] [Full Text] [PDF] |
||||
![]() |
L. M. de Heer, R. P. J. Budde, J. van Prehn, W. P. T. M. Mali, L. W. Bartels, P. R. Stella, L. A. van Herwerden, J. Kluin, and K. L. Vincken Pulsatile Distention of the Nondiseased and Stenotic Aortic Valve Annulus: Analysis With Electrocardiogram-Gated Computed Tomography Ann. Thorac. Surg., February 1, 2012; 93(2): 516 - 522. [Abstract] [Full Text] [PDF] |
||||
![]() |
A. N. Azadani, S. Chitsaz, P. B. Matthews, N. Jaussaud, J. Leung, T. Tsinman, L. Ge, and E. E. Tseng Comparison of Mechanical Properties of Human Ascending Aorta and Aortic Sinuses Ann. Thorac. Surg., January 1, 2012; 93(1): 87 - 94. [Abstract] [Full Text] [PDF] |
||||
![]() |
A. Mangini, M. G. Lemma, M. Soncini, E. Votta, M. Contino, R. Vismara, A. Redaelli, and C. Antona The aortic interleaflet triangles annuloplasty: a multidisciplinary appraisal Eur J Cardiothorac Surg, October 1, 2011; 40(4): 851 - 857. [Abstract] [Full Text] [PDF] |
||||
![]() |
C. Martin, T. Pham, and W. Sun Significant differences in the material properties between aged human and porcine aortic tissues Eur J Cardiothorac Surg, July 1, 2011; 40(1): 28 - 34. [Abstract] [Full Text] [PDF] |
||||
![]() |
C. A. Conti, A. Della Corte, E. Votta, L. Del Viscovo, C. Bancone, L. S. De Santo, and A. Redaelli Biomechanical implications of the congenital bicuspid aortic valve: A finite element study of aortic root function from in vivo data J. Thorac. Cardiovasc. Surg., October 1, 2010; 140(4): 890 - 896.e2. [Abstract] [Full Text] [PDF] |
||||
![]() |
H.-J. Schafers, F. Langer, P. Glombitza, T. Kunihara, R. Fries, and D. Aicher Aortic valve reconstruction in myxomatous degeneration of aortic valves: Are fenestrations a risk factor for repair failure? J. Thorac. Cardiovasc. Surg., March 1, 2010; 139(3): 660 - 664. [Abstract] [Full Text] [PDF] |
||||
![]() |
A. Shiran, S. Adawi, M. Ganaeem, and E. Asmer Accuracy and reproducibility of left ventricular outflow tract diameter measurement using transthoracic when compared with transesophageal echocardiography in systole and diastole Eur Heart J Cardiovasc Imaging, March 1, 2009; 10(2): 319 - 324. [Abstract] [Full Text] [PDF] |
||||
![]() |
N. Gundiah, K. Kam, P. B. Matthews, J. Guccione, H. A. Dwyer, D. Saloner, T. A.M. Chuter, T. S. Guy, M. B. Ratcliffe, and E. E. Tseng Asymmetric Mechanical Properties of Porcine Aortic Sinuses Ann. Thorac. Surg., May 1, 2008; 85(5): 1631 - 1638. [Abstract] [Full Text] [PDF] |
||||
![]() |
T. Mihaljevic, M. R. Sayeed, S. C. Stamou, and S. Paul Pathophysiology of Aortic Valve Disease , January 1, 2008; 3(2008): 825 - 840. [Full Text] |
||||
![]() |
A. Cheng, P. Dagum, and D. C. Miller Aortic root dynamics and surgery: from craft to science Phil Trans R Soc B, August 29, 2007; 362(1484): 1407 - 1419. [Abstract] [Full Text] [PDF] |
||||
![]() |
R. Fries, T. Graeter, D. Aicher, H. Reul, C. Schmitz, M. Bohm, and H.-J. Schafers In vitro comparison of aortic valve movement after valve-preserving aortic replacement J. Thorac. Cardiovasc. Surg., July 1, 2006; 132(1): 32 - 37. [Abstract] [Full Text] [PDF] |
||||
![]() |
S. Hashim and D. Richens Finite element method in cardiac surgery Interact CardioVasc Thorac Surg, February 1, 2006; 5(1): 5 - 8. [Abstract] [Full Text] [PDF] |
||||
![]() |
V. C. Sripathi, R. K. Kumar, and K. R. Balakrishnan Further insights into normal aortic valve function: role of a compliant aortic root on leaflet opening and valve orifice area Ann. Thorac. Surg., March 1, 2004; 77(3): 844 - 851. [Abstract] [Full Text] [PDF] |
||||
![]() |
M. Handke, G. Heinrichs, F. Beyersdorf, M. Olschewski, C. Bode, and A. Geibel In vivo analysis of aortic valve dynamics by transesophageal 3-dimensional echocardiography with high temporal resolution J. Thorac. Cardiovasc. Surg., June 1, 2003; 125(6): 1412 - 1419. [Abstract] [Full Text] [PDF] |
||||
![]() |
T. Mihaljevic, S. Paul, L. H. Cohn, and A. Wechsler Pathophysiology of Aortic Valve Disease , January 1, 2003; 2(2003): 791 - 810. [Full Text] |
||||
![]() |
M. K. Pasque A different kind of "total artificial heart": the interactive, computer-based human heart model Ann. Thorac. Surg., April 1, 2002; 73(4): 1032 - 1034. [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
| ANN THORAC SURG | ASIAN CARDIOVASC THORAC ANN | EUR J CARDIOTHORAC SURG |
| J THORAC CARDIOVASC SURG | ICVTS | ALL CTSNet JOURNALS |