|
|
||||||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Departments of 1Cardiac Sciences, 2Physiology and Biophysics, and 3Civil Engineering; and 4Libin Cardiovascular Institute of Alberta, University of Calgary, Calgary, Alberta, Canada
Submitted 18 January 2006 ; accepted in final form 1 June 2006
| ABSTRACT |
|---|
|
|
|---|
ventricular function, right ventricular overload; finite element analysis
Using finite element (FE) analysis, our group studied the behavior of the interventricular septum during diastole with acute pulmonary artery constriction (PAC) (17, 18). Under control conditions, the septum behaved as part of the LV and was under slight circumferential tension; the dominant stress was a low radial compression over the width of the septum. Previous FE analyses estimating normal wall stress (2, 3, 10) showed similar results.
The stress pattern changed dramatically with PAC, where RV pressure was higher than LV pressure (17, 18). Two-dimensional (2D) echocardiography showed that the septum flattened and often inverted toward the LV. Under these conditions, large quasi-circumferential compressive stresses developed in the septum, forming an arch-like pattern from the LV side of one RV insertion point to the LV side of the other, rising to the RV side of the septum at midspan. In this septal model, bending moments had to be applied at the insertion points (i.e., the ends of the model) with PAC to match the deformed shape of the model with that observed experimentally. If the entire cross section of the heart had been considered, these bending moments would have had to be balanced by equal and opposite bending moments in the ventricular free walls. Because the thickness of the LV is greater than that of the RV, it would be expected that most of these bending moments would be carried into the LV free wall, which would then become flatter. The extent of septal flattening and inversion with PAC may be dependent on the stiffness of the pericardium and the constraint it exerts. The concomitant transfer of moments into the LV free wall may also be affected by the pericardium. We used a 2D model to match the geometric information obtained from echocardiography. Results from this model, as described here, led us to blood flow experiments (6) that showed that the hydrostatic stresses induced by PAC reduced blood flow in specific areas of the myocardium; those blood flow experiments may help to explain signs of reduced blood flow in patients with normal coronary arteries.
Therefore, with the use of FE analysis, the goal of this study was to describe quantitatively the deformation of the septum and LV free wall with acute PAC and to define the role of the pericardium in this response.
| METHODS |
|---|
|
|
|---|
LV, RV, and aortic pressures were measured with 8-Fr micromanometer-tipped catheters (Millar Instruments, Houston, TX), which were inserted via peripheral vessels. Pericardial pressures were recorded with flat, fluid-filled balloon catheters loosely attached to the epicardium (12). The heart was repositioned in the pericardium, and the edges of the pericardium were reapproximated with loose interrupted sutures (21).
A pneumatic constrictor (1416 mm, In-Vivo Metrics, Healdsburg, CA) was placed on the pulmonary artery to increase RV pressure transiently and change the end-diastolic transseptal pressure gradient [LV end-diastolic pressure (LVEDP) RV end-diastolic pressure (RVEDP)]. Pressures were recorded simultaneously with 2D echocardiography (Sonos 1000, Hewlett-Packard, Palo Alto, CA); an LV minor-axis view at the level of the papillary muscles was recorded by using a 2.5-MHz transesophageal probe held on the surface of the RV free wall. The shapes and dimensions used for the FE analysis (i.e., LV and RV diameters and wall thicknesses) were derived from analysis of the echocardiographic images (Echo Analysis version 1.05, Advanced Measurements; Calgary, AB, Canada). Examples of the echo data used are shown in Fig. 1. End-diastolic echocardiographic images were digitized by manually tracing the endocardium and epicardium of the LV and the RV endocardial surface of the septum and its insertion points. The curvatures (i.e., 1/R, where R is the radius of curvature measured from the echo images) of the septum and the LV free wall were calculated. An LV area index was also calculated from the echo images as the product of the anterior-posterior and septum-to-free wall epicardial diameters.
|
FE model. A 2D FE model of an equatorial plane of the heart was developed by using PATRAN (version 5.4; MSC Software, Santa Ana, CA) and is shown in Fig. 2. The thicknesses of the LV and RV walls were typical of values at end diastole. The pericardium was modeled as congruent with the LV and RV free walls, separated from them tangentially near the RV insertion points. Slide-line elements were used to separate the free walls from the pericardium and to allow free sliding movement within the pericardium. All solid elements were eight-node quadrilaterals. Over 1,600 elements were used in the model. One node was fixed in space to allow the model to converge. Hyperelasticity was assumed for all materials, and geometric nonlinearity was utilized, because large displacements were expected.
|
The 2D model was chosen for its simplicity and convenience. As illustrated by our several previous experimental studies of ventricular interaction (1, 9, 11, 14, 23, 24), minor-axis dimensions can be used to analyze the phenomena of ventricular interaction in ways that are both useful and consistent.
The model was solved using Abaqus (version 6.27, Abaqus; Pawtucket, RI). As shown in the previous study (18), both linear and nonlinear stress-strain relationships produced similar results after analysis. Linear analysis was used because it required much less computer time.
Data recorded under control conditions and during PAC, with the pericardium closed and with the pericardium opened, were analyzed and modeled.
| RESULTS |
|---|
|
|
|---|
|
Results of the FE model for control conditions for one dog at end diastole are shown in Fig. 3A and Fig. 4, A and C. Under control conditions (pericardium closed), the displaced mesh is very similar to the initial unstressed mesh (Fig. 3A); the LV (including the septum) was quasi-circular.
|
|
As seen in the previous septum-only model (18), the dominant stress was radial compression (Fig. 4, A and C), with the inner endocardial layer showing greater compression than the outer epicardial layer. The pattern was uniform throughout the wall. The maximal compression seen was 10.5 mmHg (or 1.4 kPa; negative values indicate compression and positive values indicate tension), whereas mean pixel-average radial compression ranged between 3.2 ± 0.7 and 4.3 ± 0.7 mmHg. To obtain a solution for the model, we used a pericardial stiffness (0.060.090 MPa, or 452677 mmHg) that was an order of magnitude higher than the myocardial (muscle) stiffness (2246 mmHg). The pericardium was subject to circumferential tension of >150 mmHg (Fig. 4A). Tensions at the RV insertion points were low. The stress distribution was similar in all hearts analyzed.
With PAC, RVEDP exceeded LVEDP: RVEDP was 8.9 ± 2.8 mmHg, whereas LVEDP was 5.5 ± 2.5 mmHg, making the transseptal pressure gradient negative (3.4 ± 1.7 mmHg; Table 1). These values are similar to those observed by other groups in previous studies (5, 22). To reach a solution with FE analysis for PAC, the stiffness of the myocardium was increased slightly compared with control (i.e., it ranged from 0.006 to 0.009 MPa). End-diastolic pericardial pressure increased from 6.9 ± 0.8 to 8.5 ± 0.8 mmHg with PAC (P < 0.01). Pericardial model stiffness was also increased (ranging from 0.09 to 0.12 MPa).
Results of the FE model during PAC (pericardium closed) for one dog are shown in Fig. 3B and Fig. 4, B and D. The displaced mesh (Fig. 3B) shows the enlarged RV and the flattened and/or inverted septum that occurs with PAC. The LV free wall also shows some flattening when compared with control, a result of bending moments induced by the flattened septum.
Figure 4B shows a contour plot of the maximal principal stresses (exhibiting the greatest tensions), and Fig. 4D shows a contour plot of minimal principal stresses (exhibiting the greatest compressions). The circumferential tension seen in the pericardium with PAC was significantly higher than seen with control (Fig. 4, comparing A and B). Within the constraint of the pericardium, the septum flattened and, except for two dogs, inverted, with compressive circumferential stresses developing out to the insertion points (Fig. 4D). Again, this is similar to the results from the septum-only model (18). Rather than seeing a constant curvature (as with control conditions), there was a sharp bend near the RV insertions and flattening and inversion of the septum. The end-diastolic curvature of the septum changed from 0.276 ± 0.015 cm1 (concave to the LV) to 0.125 ± 0.011 cm1 with PAC (inverted into the LV, concave to RV; P < 0.01). The inverted septum and the larger RV size (Fig. 3B) caused greater pull and increased tension at the RV insertion points (Fig. 4, B and D). Not including the extreme tensions at the RV insertion that may be model dependent, the average stress over the region of the epicardial RV insertion points was seen to change from 4.1 ± 0.6 at control to +22.2 ± 1.4 mmHg with PAC. Maximal tensions in this region were as great as 75.8 ± 5.9 mmHg. The inverted septum also caused a zone of tension (to develop) on the endocardial side of the midseptum (with PAC) (regional average pixel stress changed from 5.4 ± 0.6 at control to +19.7 ± 2.0 mmHg; P < 0.01). The maximum individual-pixel tensile stress in this region was 36.8 ± 3.2 mmHg.
The maximum individual-pixel compressive stress (98.3 mmHg or 13.0 kPa) occurred on the LV side of the insertion points (Fig. 4D; minimal principal stresses). This maximal compression is almost 10 times that seen in the control model (Fig. 4, C vs. D), and this compression was circumferentially oriented rather than radially. The maximum individual-pixel compressive stress in the endocardial insertion points increased significantly from 5.3 ± 0.5 at control to 75.0 ± 7.7 mmHg during PAC.
Another zone of significant compression developed on the RV side of the inverted septum, although these stresses were not as great as those at the insertion points (Fig. 4D). The maximum individual-pixel compression on the RV side of the septum increased from 4.1 ± 0.7 to 46.6 ± 5.0 mmHg. The circumferential stresses extended beyond the insertions, into the LV free wall. The LV free wall flattened significantly (P < 0.01) but never inverted, with the curvature at midspan changing from 0.349 ± 0.023 (control) to 0.274 ± 0.015 cm1 (PAC). The average circumferential compression at the midpoint of the LV free wall epicardium, where flattening was most obvious, changed from 3.7 ± 0.7 to 28.1 ± 9.8 mmHg (P < 0.01). This stress distribution was similar for all hearts analyzed.
The hydrostatic component of stress is the mean of the principal stresses and tends to cause changes in tissue volume, as opposed to changes in shape, which are caused by changes in the deviatoric components of stress. Plots of the hydrostatic component of stress were shown in our previous study (6). Consistent with the results shown in the present paper, the hydrostatic stress plots showed regions of compression at the endocardial RV insertion points as well as the RV side of the midseptum.
The bending moments around the circumference of the LV were calculated from data recorded under normal conditions and during PAC; the FE model results are illustrated in Fig. 5 for an individual pericardium-closed case. Under control conditions (Fig. 5A), the bending moments around the LV circumference were small. Slight tension at the RV insertion points was consistent with the inflections in the LV endocardial circumference. Axial forces were uniform. With PAC, the patterns were magnified. Large bending moments developed in the septum, and tensions at the RV insertion points increased because of larger RV circumferential forces. As a result, the midpoint of the septum was displaced toward the LV, and, in this case, the septum inverted. The LV free wall flattened (Fig. 3, A vs. B) as a result of the induced bending.
|
|
With the pericardium open, PAC still produced a negative transseptal gradient (4.2 ± 1.1 mmHg), although LVEDP and RVEDP were significantly lower than during PAC with the pericardium closed. LVEDP decreased from 5.2 ± 0.9 (control, open pericardium) to 3.0 ± 0.6 mmHg during PAC while RVEDP increased from 4.0 ± 1.1 (control, open pericardium) to 7.1 ± 1.2 mmHg (see Table 1). During PAC with the pericardium open, the RV showed a large increase in size. In all cases, the septum flattened and in eight dogs it inverted, but there was no change in the LV free wall curvature associated with this flattening and/or inversion (Fig. 3D).
FE results of a PAC for one dog with the pericardium opened are shown in Fig. 6, B and D. When the pericardium was opened, the septum inverted significantly less with PAC than when it was closed. In five of the 14 dogs, the septum only flattened, whereas in the others, it inverted. Mean septal curvature during PAC was 0.007 ± 0.059 cm1, compared with 0.125 ± 0.011 cm1 during PAC with the pericardium closed.
Arch-like circumferential stresses still developed in the septum but were smaller than when the pericardium was closed (compare Fig. 6D with Fig. 4D). The maximal compressive stresses were again seen on the LV endocardial side of the insertion points but were smaller than with the pericardium-closed model. For this individual case (Fig. 6), a maximum compression of 48.6 mmHg occurred again at the insertion points, compared with 98.3 mmHg when the pericardium was closed. The maximum individual-pixel compression with PAC (22.7 ± 1.4) was also reduced with pericardium open (compared with 75.0 ± 7.7 for the pericardium-closed case). The compression within the septum (midpoint, RV side) was also reduced in magnitude over the pericardium-closed model. The overall size of the compression zones at both the insertion points and septum appeared to be reduced (Fig. 6D compared with the pericardium-closed model in Fig. 4D).
When the pericardium was open, PAC did not flatten the LV free wall significantly (0.326 ± 0.019 cm1 compared with 0.318 ± 0.020 cm1 for control, pericardium opened; not significantly different) nor did stresses change, in spite of substantial septal deformation.
| DISCUSSION |
|---|
|
|
|---|
Although our data were obtained during acute RV pressure loading, our findings may also apply to chronic RV pressure loading that is associated with similar septal flattening and inversion (16, 19). Substantial pericardial remodeling occurs with experimental intracardiac shunting (18) and in patients with chronic RV pressure overload, such as in chronic embolic thromboembolic pulmonary hypertension (19). However, those patients had enlarged RVs, their septums flattened (i.e., moved leftward during diastole), and, after embolectomy, septal motion returned to normal. This would suggest that stress patterns in those patients were not substantially different from those we observed in this study. Stress patterns with the pericardium open also support this: PAC still caused septal flattening, albeit to a lesser degree. Thus pericardial remodeling may attenuate the effects; however, it is likely that decreasing the transseptal pressure gradient still deforms the heart and causes similar stresses.
During PAC with the pericardium closed, the LV free wall flattened in a similar manner as the septum, although not to the same degree. This appears to be the result of the bending moments induced by the septal inversion. When the pericardium was opened, however, LV free wall flattening did not occur, probably because the bending moments were smaller than when the pericardium was closed and the angle of pull of the RV was increased. The effect of the changes caused by the shift in RV pull can be seen by comparing Figs. 4 and 7. With PAC at end diastole when the pericardium is taut, pericardial constraint causes a greater negative transseptal pressure gradient, greater septal inversion, and thus greater bending moments to be transferred to the LV free wall, compared with when the pericardium is open. Therefore, it is clear that data from studies performed with an open pericardium and FE or other analytic models may not accurately reflect structural and mechanical responses of the organ and/or individual tissues when the pericardium is normally closed. The potential role of the pericardium should be considered in any experimental study or analytic model of the heart.
To reach a solution under PAC conditions, the stiffnesses of the muscle and pericardium had to be increased, compared with the values needed under control conditions. Such increases had also been required with analysis of the septum-only model (18). We therefore hypothesized that the compressive stresses that developed with PAC might reduce coronary blood flow, perhaps by impeding the drainage of venous blood. Congestion of blood within the wall has been shown by others (19, 25) to increase tissue stiffness. If compression reduced blood flow, it might help to explain the apparent perfusion defects in patients with normal coronary arteries who have pulmonary hypertension, left bundle branch block, or other conditions that elevate RV pressure (13, 16, 20).
A coronary perfusion study using colored microspheres was performed to test this hypothesis . Under control conditions with the pericardium closed, there was similar flow in all regions of the LV, and these flow values compared well with previously published flow data (4, 8). With PAC, however, the flow pattern changed significantly. The areas showing the greatest decrease in flow occurred at the endocardial side of the RV insertion points as well as on the RV side of the septum. These zones of high flow reduction compared well with the zones of high circumferential compressive stress shown in Fig. 4D. The data show a strong correspondence between the circumferential compressive stresses and blood flow. Compressive stresses are greater on the endocardial side than the epicardial side, and this was also where there was reduced flow (except in zones opposite a flattening/inversion zone). This transmural difference in flow has also been demonstrated by other research groups (8, 22) but never explained in terms of compression changes. Thus our results provide a possible mechanism for septal hypoperfusion with RV pressure overload.
The following points summarize our conclusions.
1) With the pericardium closed, bending moments induced in the septum during acute PAC are transmitted into the LV free wall, flattening that structure. The decrease in constraint due to opening the pericardium allows the RV to enlarge considerably, changes the "pull" angle of the RV on the LV wall, and reduces the moments induced in the flattened (inverted) septum and in the LV free wall. With the pericardium open, the LV free wall does not flatten like it does during PAC with the pericardium closed.
2) The pericardium is a major modulator of the behavior of the structural elements of the heart during diastole. Its presence substantially affects the deformation of the heart during diastole and, consequentially, strains and stresses. Our results suggest that data from experiments and models in which the pericardium is open should be reconsidered in this light.
| GRANTS |
|---|
|
|
|---|
| ACKNOWLEDGMENTS |
|---|
C. Gibbons Kroeker held a postdoctoral research fellowship from the Natural Science and Engineering Research Council (Ottawa), N. G. Shrive is a Killam Memorial Professor, and J. V. Tyberg was a Scientist of the Alberta Heritage Foundation for Medical Research (Edmonton).
| FOOTNOTES |
|---|
The costs of publication of this article were defrayed in part by the payment of page charges. The article must therefore be hereby marked "advertisement" in accordance with 18 U.S.C. Section 1734 solely to indicate this fact.
| REFERENCES |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
P. S. Halvorsen, L. A. Fleischer, A. Espinoza, O. J. Elle, L. Hoff, H. Skulstad, T. Edvardsen, and E. Fosse Detection of myocardial ischaemia by epicardial accelerometers in the pig Br. J. Anaesth., January 1, 2009; 102(1): 29 - 37. [Abstract] [Full Text] [PDF] |
||||
![]() |
P. S. Halvorsen, A. Espinoza, L. A. Fleischer, O. J. Elle, L. Hoff, R. Lundblad, H. Skulstad, T. Edvardsen, H. Ihlen, and E. Fosse Feasibility of a three-axis epicardial accelerometer in detecting myocardial ischemia in cardiac surgical patients. J. Thorac. Cardiovasc. Surg., December 1, 2008; 136(6): 1496 - 1502. [Abstract] [Full Text] [PDF] |
||||
![]() |
P. P. Sengupta, V. K. Krishnamoorthy, W. P. Abhayaratna, J. Korinek, M. Belohlavek, T. M. Sundt III, K. Chandrasekaran, F. Mookadam, J. B. Seward, A. J. Tajik, et al. Disparate patterns of left ventricular mechanics differentiate constrictive pericarditis from restrictive cardiomyopathy. J. Am. Coll. Cardiol. Img., January 1, 2008; 1(1): 29 - 38. [Abstract] [Full Text] [PDF] |
||||
![]() |
Y. Liu, W. Zhang, and G. S. Kassab Effects of myocardial constraint on the passive mechanical behaviors of the coronary vessel wall Am J Physiol Heart Circ Physiol, January 1, 2008; 294(1): H514 - H523. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
| Visit Other APS Journals Online |