Method Article

Modellazione a parametro raggruppato e elemento finito dell'insufficienza cardiaca con frazione di espulsione preservata

DOI:

10.3791/62167

⸱

February 13th, 2021

 ,  , 
1Institute for Medical Engineering and Science, Massachusetts Institute of Technology, 2Health Science and Technology Program, Harvard/Massachusetts Institute of Technology, 3Department of Mechanical Engineering, Massachusetts Institute of Technology

In This Article

Summary

$$\rightleftharpoonup{xx}$$ $$\longleftharp{xx}$$, $$\longrightharp{xx}$$,

Questo lavoro introduce due modelli computazionali di insufficienza cardiaca con frazione di espulsione preservata basata su un approccio a parametro forfait e analisi degli elementi finiti. Questi modelli sono utilizzati per valutare i cambiamenti nell'emodinamica del ventricolo sinistro e la relativa vascolarizzazione indotta dal sovraccarico di pressione e dalla ridotta conformità ventricolare.

Abstract

$$\rightleftharpoonup{xx}$$ $$\longleftharp{xx}$$, $$\longrightharp{xx}$$,

Gli sforzi scientifici nel campo della modellazione computazionale delle malattie cardiovascolari si sono in gran parte concentrati sull'insufficienza cardiaca con ridotta frazione di espulsione (HFrEF), trascurando ampiamente l'insufficienza cardiaca con frazione di espulsione preservata (HFpEF), che è diventata più recentemente una forma dominante di insufficienza cardiaca in tutto il mondo. Motivati dalla scarsità di HFpEF nelle rappresentazioni silico, in questo documento vengono presentati due distinti modelli computazionali per simulare l'emodinamica dell'HFpEF derivante dal sovraccarico di pressione ventricolare sinistro. In primo luogo, un modello a parametro raggruppato orientato agli oggetti è stato sviluppato utilizzando un risolutore numerico. Questo modello si basa su una rete a forma di Windkessel a dimensione zero (0D), che dipende dalle proprietà geometriche e meccaniche degli elementi costitutivi e offre il vantaggio di bassi costi computazionali. In secondo luogo, è stato utilizzato un pacchetto software di analisi degli elementi finiti (FEA) per l'implementazione di una simulazione multidimensionale. Il modello FEA combina modelli multifisica tridimensionali (3D) della risposta cardiaca elettromeccanica, deformazioni strutturali ed emodinamica basata sulla cavità fluida e utilizza un modello semplificato di parametro raggruppato per definire i profili di scambio del flusso tra diverse cavità del fluido. Attraverso ogni approccio, sono stati simulati con successo sia i cambiamenti emodinamici acuti che cronici nel ventricolo sinistro e nella vascuola prossimale risultanti dal sovraccarico di pressione. In particolare, il sovraccarico di pressione è stato modellato riducendo l'area dell'orifizio della valvola aortica, mentre il rimodellamento cronico è stato simulato riducendo la conformità della parete ventricolare sinistra. Coerentemente con la letteratura scientifica e clinica dell'HFpEF, i risultati di entrambi i modelli mostrano (i) un'elevazione acuta del gradiente di pressione transaortica tra il ventricolo sinistro e l'aorta e una riduzione del volume dell'ictus e (ii) una diminuzione cronica del volume ventricolare sinistro diastolico finale, indicativo di disfunzione diastolica. Infine, il modello FEA dimostra che lo stress nel miocardio HFpEF è notevolmente più elevato che nel tessuto cardiaco sano durante tutto il ciclo cardiaco.

Introduction

$$\rightleftharpoonup{xx}$$ $$\longleftharp{xx}$$, $$\longrightharp{xx}$$,

L'insufficienza cardiaca è una delle principali cause di morte in tutto il mondo, che si verifica quando il cuore non è in grado di pompare o riempire adeguatamente per stare al passo con le esigenze metaboliche del corpo. La frazione di espulsione, cioè la quantità relativa di sangue immagazzinata nel ventricolo sinistro che viene espulso ad ogni contrazione viene utilizzata clinicamente per classificare l'insufficienza cardiaca in (i) insufficienza cardiaca con frazione di espulsione ridotta (HFrEF) e (ii) insufficienza cardiaca con frazione di espulsione preservata (HFpEF), per frazioni di espulsione inferiori o superiori al 45%,rispettivamente 1,2,3. I sintomi dell'HFpEF si sviluppano spesso in risposta al sovraccarico di pressione ventricolare sinistra, che può essere causato da diverse condizioni tra cui stenosi aortica, ipertensione e ostruzione del tratto di deflusso ventricolaresinistro 3,4,5,6,7. Il sovraccarico di pressione guida una cascata di aberrazioni molecolari e cellulari, portando all'ispessimento della parete ventricolare sinistra (rimodellamento concentrico) e, in definitiva, all'irrigidimento o alla perdita diconformità delle pareti 8,9,10. Questi cambiamenti biomeccanici influenzano profondamente l'emodinamica cardiovascolare in quanto si traducono in un'elevata relazione pressione-volume diastolico finale e in una riduzione del volume end-diastolico11.

La modellazione computazionale del sistema cardiovascolare ha avanzato la comprensione delle pressioni e dei flussi sanguigni sia in fisiologia che in malattia e ha favorito lo sviluppo di strategie diagnostiche e terapeutiche12. Nei modelli silico sono classificati in modelli a bassa o alta dimensione, con i primi che utilizzano metodi analitici per valutare proprietà emodinamiche globali a bassa domanda computazionale e il secondo fornisce una descrizione multiscala e multifisica più ampia della meccanica cardiovascolare e dell'emodinamica nel dominio 2D o 3D13. La rappresentazione di Windkessel a parametro forfait è la più comune tra le descrizioni a bassa dimensione. Basato sull'analogia del circuito elettrico (legge di Ohm), questo imita il comportamento emodinamico complessivo del sistema cardiovascolare attraverso una combinazione di elementi resistivi, capacitivi e induttivi14. Un recente studio di questo gruppo ha proposto un modello windkessel alternativo nel settore idraulico che consente la modellazione di cambiamenti nella geometria e nella meccanica di grandi vasi-camere cardiache e valvole-in modo più intuitivo rispetto ai tradizionali modelli analogici elettrici. Questa simulazione è sviluppata su un risolutore numerico orientato agli oggetti (vedi tabella dei materiali)e può catturare la normale emodinamica, gli effetti fisiologici dell'accoppiamento cardiorespiratorio, il flusso sanguigno guidato dalle vie respiratorie nella fisiologia a cuore singolo e i cambiamenti emodinamici dovuti alla costrizione aortica. Questa descrizione espande le capacità dei modelli a parametro raggruppato offrendo un approccio fisicamente intuitivo per modellare uno spettro di condizioni patologiche tra cui l'insufficienzacardiaca 15.

I modelli ad alta dimensione si basano su FEA per calcolare l'emodinamica spatiotemporale e le interazioni fluido-struttura. Queste rappresentazioni possono fornire descrizioni dettagliate e accurate del campo sanguigno locale; tuttavia, a causa della loro bassa efficienza computazionale, non sono adatti per studi sull'intero alberocardiovascolare 16,17. Un pacchetto software (vedi tabella dei materiali)è stato utilizzato come piattaforma FEA anatomicamente accurata del cuore umano adulto a 4 camere, che integra la risposta elettromeccanica, le deformazioni strutturali e l'emodinamica basata sulla cavità fluida. Il modello di cuore umano adattato comprende anche un semplice modello a parametro grumoso che definisce lo scambio di flusso tra le diverse cavità del fluido, nonché una caratterizzazione meccanica completa del tessutocardiaco 18,19.

Sono stati formulati diversi modelli di insufficienza cardiaca a parametro forfait e FEA per catturare anomalie emodinamiche e valutare strategie terapeutiche, in particolare nel contesto di dispositivi meccanici di assistenza circolatoria per HFrEF20,21,22,23,24. Un'ampia gamma di modelli di parametri raggruppati 0D di varie complessità ha quindi catturato con successo l'emodinamica del cuore umano in condizioni fisiologiche e HFrEF attraverso l'ottimizzazione di sistemi elettrici analogici Windkessel a due otre elementi 20,21,23,24. La maggior parte di queste rappresentazioni sono modelli uni- o biventricolari basati sulla formulazione di elastanza variabile nel tempo per riprodurre l'azione contrattile del cuore e utilizzare una relazione pressione-volume end-diastolica non lineare per descrivere il riempimento ventricolare25,26,27. Modelli completi, che catturano la complessa rete cardiovascolare e imitano sia l'azione di pompaggio atriale che ventricolare, sono stati utilizzati come piattaforme per il test dei dispositivi. Tuttavia, sebbene esista un significativo corpus letterario intorno al campo dell'HFrEF, pochissimi nei modelli silico di HFpEF sono statiproposti 20,22,28,29,30,31.

Un modello a bassa dimensionalità dell'emodinamica HFpEF, recentemente sviluppato da Burkhoffetal. Inoltre, utilizzano la loro piattaforma in silico per valutare la fattibilità di un dispositivo circolatorio meccanico per HFpEF, pionieristica ricerca computazionale di HFpEF per studi fisiologici e sviluppo di dispositivi. Tuttavia, questi modelli rimangono incapaci di catturare i cambiamenti dinamici nei flussi sanguigni e nelle pressioni osservate durante la progressione della malattia. Un recente studio di Kadry et al. Il loro lavoro fornisce un'analisi emodinamica completa della disfunzione diastolica basata sulle proprietà attive e passive del miocardio. Allo stesso modo, la letteratura dei modelli ad alta dimensione si è concentrata principalmente su HFrEF19,33,34,35,36,37. 33 ha proposto un modello FEA di fluido cardiaco-elettromeccanica completamente accoppiato per prevedere il profilo emodinamico HFrEF e l'efficacia di un dispositivo di assistenza ventricolare sinistro (LVAD). Questo modello biventricolare (o a due camere) utilizzava un circuito accoppiato di Windkessel per simulare l'emodinamica del cuore sano, HFrEF e HFrEF con supporto LVAD33,37.

Allo stesso modo, Sack etal. La loro geometria biventricolare è stata ottenuta dai dati di risonanza magnetica (MRI) di un paziente, e la mesh a elementi finiti del modello è stata costruita usando la segmentazione dell'immagine per analizzare l'emodinamica di un ventricolo destro in fallimento supportato da VAD35. Sono stati sviluppati approcci cardiaci FEA a quattro camere per migliorare l'accuratezza dei modelli del comportamento elettromeccanicodel cuore 19,34. In contrasto con le descrizioni biventricolari, i modelli a quattro camere derivati dalla risonanza testrizzazione del cuore umano forniscono una migliore rappresentazione dell'anatomiacardiovascolare 18. Il modello cardiaco utilizzato in questo lavoro è un esempio consolidato di modello FEA a quattro camere. A differenza dei modelli FEA a parametro forfait e biventricolare, questa rappresentazione cattura i cambiamenti emodinamici man mano che si verificano durante la progressionedella malattia 34,37. 34 , ad esempio, utilizzavala stessa piattaforma per implementare un modello numerico di crescita del rimodellamento osservato in HFrEF e HFpEF. Tuttavia, questi modelli valutano gli effetti dell'ipertrofia cardiaca solo sulla meccanica strutturale e non forniscono una descrizione completa dell'emodinamica associata.

Per far fronte alla mancanza di HFpEF nei modelli silico in questo lavoro, il modello a parametro forfait precedentemente sviluppatoda questo gruppo 15 e il modello FEA sono stati riadattati per simulare il profilo emodinamico di HFpEF. A tal fine, verrà dimostrata per la prima volta la capacità di ogni modello di simulare l'emodinamica cardiovascolare al basale. Verranno quindi valutati gli effetti del sovraccarico di pressione ventricolare sinistro indotto dalla stenosi e della ridotta conformità ventricolare sinistra dovuta al rimodellamento cardiaco- un tipico segno distintivo dell'HFpEF.

Access restricted. Please log in or start a trial to view this content.

Protocol

$$\rightleftharpoonup{xx}$$ $$\longleftharp{xx}$$, $$\longrightharp{xx}$$,

1. Modello a parametro forfait 0D

  1. Configurazione della simulazione
    NOTA: Nell'ambiente del risolutore numerico (vedere la tabella dei materiali), costruire il dominio come illustrato nella figura 1. Questo è composto dal cuore a 4 camere, dalla parte superiore del corpo, dall'addome, dalla parte inferiore del corpo e dai compartimenti toracici, così come dalla vascucolatura prossimale, tra cui l'aorta, l'arteria polmonare e la vena cavae superiore e inferiore. Gli elementi standard utilizzati in questa simulazione fanno parte della libreria idraulica predefinita. I dettagli sono disponibili nei file supplementari.
    1. Naviga nella libreria idraulica per trovare gli elementi richiesti: tubazione idraulica, camera idraulica a volume costante, resistenza lineare, pompa centrifuga, valvola di ritegno, orifizio ad area variabile e fluido idraulico personalizzato.
      1. Far cadere gli elementi della tubazione idraulica nel workspace.
        NOTA: Questi rappresentano perdite di attrito, nonché conformità alle pareti e comprimibilità dei fluidi nei vasi sanguigni e nelle camere cardiache. Attraverso questo blocco, la perdita di pressione viene calcolata usando la legge di Darcy-Weisbach, mentre il cambiamento di diametro dovuto alla conformità alle pareti dipende dalla costante di proporzionalità di conformità, dalla pressione luminare e dalla costante di tempo. Infine, la comprimibilità del fluido è definita dal modulo di massa del mezzo.
      2. Inserire gli elementi della camera idraulica a volume costante per definire la conformità delle pareti e la comprimibilità del fluido.
        NOTA: Questo blocco non tiene conto delle perdite di pressione dovute all'attrito.
      3. Aggiungere gli elementi di resistenza lineare per definire la resistenza al flusso.
        NOTA: Questo è indipendente dalle proprietà geometriche della vascucolatura, analogamente all'elemento resistivo utilizzato nei modelli elettrici analogici Windkessel. Altri blocchi, come la pompa centrifuga, la valvola di ritegno, l'orifizio ad area variabile e gli elementi del fluido idraulico personalizzati devono essere inseriti per generare l'input di pressione desiderato nel sistema, modellare gli effetti delle valvole cardiache sul flusso sanguigno e definire le proprietà meccaniche del sangue. Attraverso questi elementi, il comportamento del sistema cardiovascolare sia in fisiologia che in malattia può essere pienamente catturato. Il segnale di ingresso per la pompa centrifuga è disponibile nella figura S1A.
      4. Modellare la contrattilità di ogni camera cardiaca attraverso l'elemento della camera di conformità variabile personalizzato.
        NOTA: accetta la conformità come segnale di ingresso definito dall'utente variabile nel tempo e si basa sul modello di elastanza variabile nel tempo (Figura S1B-D).
    2. Fornire i parametri relativi ad ogni elemento, come mostrato nella tabella S1, che si trova anche in Rosalia etal.
    3. Inserire un elemento Sequenza ripetuta segnale fisico (PS) per ciascuno dei blocchi che richiedono un segnale di ingresso definito dall'utente variabile nel tempo: la pompa LV, gli elementi di conformità a conformità variabile e i blocchi orifizi ad area variabile.
      NOTA: I segnali di ingresso utilizzati per questa simulazione sono disponibili nella figura S1.
    4. Selezionare il risolutore implicito ODE 23t predefinito ed eseguire la simulazione per 100 s per raggiungere uno stato stazionario.

2. Il modello FEA

  1. Configurazione della simulazione
    NOTA: il modello FEA utilizza un'analisi elettro-meccanica accoppiata in sequenza. In questo modello, l'analisi elettrica viene condotta per prima; quindi i potenziali elettrici risultanti vengono utilizzati come fonte di eccitazione nella seguente analisi meccanica. Pertanto, l'impostazione di simulazione contiene due domini di lavoro: i domini elettrico (ELEC) emeccanico( MECH ), predefiniti nel software di simulazione FEA (Table of Materials)18. Pertanto, la sezione seguente descrive solo il flusso di lavoro di analisi. Il modello FEA utilizza le seguenti subroutine utente HETVAL, VUANISOHYPERe UAMP per la modellazione di materiali elettrici emeccanici 18.
    1. Spostarsi nel dominio ELEC per eseguire l'analisi elettrica utilizzando la procedura di temperatura predefinita nel modulo Standard.
      1. Utilizzare un passaggio di analisi singola denominato BEAT. Impostare la durata del ciclo cardiaco su 500 ms e applicare un impulso potenziale elettrico a un set di nodi che rappresenta il nodo sinoatriale (SA)(set di nodi: R_Atrium-1.SA_NODE).
      2. Esaminare la forma d'onda elettrica predefinita, che varia da -80 mV a 20 mV su 200 ms con la definizione di ampiezza del passo uniforme, come descritto nella guidadel modello 18. Utilizzate i valori di default delle costanti del materiale nell'analisi elettrica per regolare il ritardo AV.
      3. Avviare il modulo Processo e creare un processo denominato heart-elec.
    2. Una volta completata la configurazione dell'analisi elettrica, navigare nel dominio MECH per eseguire l'analisi meccanica basata sulla cavità del fluido.
      NOTA: La simulazione meccanica viene eseguita dopo l'analisi elettrica e i potenziali elettrici risultanti vengono utilizzati come fonte di eccitazione per l'analisi meccanica. L'analisi meccanica contiene più passaggi.
      1. Utilizzare i tre passaggi principali denominati PRE-LOAD, BEAT1e RECOVERY1. Nella fase PRE-LOAD, esaminare le condizioni limite dello stato pre-stressato del cuore. Utilizzare 0,3 s come tempo di passaggio per aumentare linearmente la pressione nelle camere del fluido.
        NOTA: I valori predefiniti della pressione della cavità del fluido sono riportati nella tabella S3. Lo stato pre-stressato del cuore era già definito nella normale configurazione di simulazione cardiaca e le condizioni iniziali del nodo sono fornite nei file di simulazione esterni, come elencato nella tabella S5. Il ricalcolo dello stato di sollecitazione zero mediante la simulazione meccanica inversa è richiesto ogni volta che la condizione limite viene modificata, come spiegato nei passaggi 3.2.2-3.2.4.
      2. Nel passaggio BEAT1, utilizzate 0,5 s come tempo di passaggio per simulare la contrazione.
      3. Nella fase RECOVERY1, selezionare 0,5 s per il rilassamento cardiaco e il riempimento ventricolare per una frequenza cardiaca di 60 bpm.
      4. Abilitare i passaggi successivi, BEATX e RECOVERYX, per simulare più di un ciclo cardiaco per raggiungere uno stato stazionario.
        NOTA: Tre cicli cardiaci saranno sufficienti per raggiungere lo stato stazionario. Un ciclo della simulazione viene completato in ~ 8 h su un processore a 24 core (3,2 GHz x 24).
      5. Avviare il modulo Processo e creare un processo denominato heart-mech, abilitando l'opzione di doppia precisione.
  2. Esaminare il modello Windkessel a parametro forfait semplificato
    NOTA: Il dominio meccanico del modello FEA ha un modello di flusso sanguigno, che si basa su un circuito semplificato a parametro forfait e viene creato come combinazione di cavità di fluido a base di superficie e scambi di fluidi come si vede nella figura 218.
    1. Utilizzate la rappresentazione di Windkessel menzionata nella nota precedente per eseguire la simulazione. Esaminare la rappresentazione del modello del flusso sanguigno per regolare i valori degli elementi resistivi e capacitivi rispettivamente per le resistenze al flusso e le conformità strutturali.
    2. Esaminare la rappresentazione degli elementi finiti 3D di quattro camere cardiache e assicurarsi che le loro posizioni geometriche siano accurate.
    3. Controllare l'assemblaggio del cuore e passare al modulo Interazione per regolare i valori di conformità e contrattilità di ciascuna delle quattro camere cardiache.
      NOTA: i valori predefiniti nel modulo Interazione sono configurati per simulare un ciclo di battito cardiaco umano sano idealizzato18.
    4. Esaminare le seguenti cavità del fluido idrostatico nel modulo interazione, CAV-AORTA, CAV-LA, CAV-LV, CAV-PULMONARY_TRUNK, CAV-RA, CAV-RV, CAV-SVC, CAV-ARTERIAL-COMP, CAV-PULMONARY-COMP e CAV-VENOUS-COMP(Tabella S3).
    5. Utilizzare le camere di conformità (CAV-ARTERIAL-COMP, CAV-PULMONARY-COMP e CAV-VENOUS-COMP) come volumi cubici in quanto rappresentano la conformità delle circolazioni arteriose, venose e polmonari.
    6. Collegare tre volumi cubici di conformità a una molla messa a terra e rivedere il valore di rigidità per modellare la risposta pressione-volume nelle circolazioni arteriose, venose e polmonari.
    7. Controllare le seguenti definizioni di scambio di fluidi tra le cavità del fluido idrostatico: arterioso-venoso, atrio venoso-destro, ventricolo destro-destro, sistema ventricolo-polmonare destro, atrio polmonare-sinistro del sistema, ventricolo sinistro-atrio-sinistro e ventricolo-aorta sinistro (Tabella S4).
    8. Regolare il coefficiente di resistenza viscosa per modificare il modello del flusso sanguigno in ogni collegamento di scambio di fluidi (per ulteriori informazioni sull'effetto di resistenza viscosa), vedere i file supplementari.
  3. Simulazione multifisica
    1. Individuare il file di database CAE nella directory di lavoro.
      NOTA: il modello FEA in questo protocollo viene recapitato nel database ed è denominato LH-Human-Model-Beta-V2_1.cae.
    2. Inserire i file di input, oggetto e libreria nella directory di lavoro per eseguire la simulazione. Vedere la tabella S5 per l'elenco completo dei file di input e raccolta.
    3. Avviare il software di simulazione del modello FEA (vedere table of materials).
      NOTA: consultare il provider di software per verificarne la compatibilità con le versionisuccessive 18.
    4. Esaminare le parti, l'assemblaggio e le condizioni limite nei domini ELEC e MECH, come descritto nelle sezioni 2.2 e 2.3.
    5. In primo luogo, eseguire il processo di simulazione elettrica denominato heart-elec, come descritto nella sezione 2.1.1.3. Ispezionare visivamente i risultati del potenziale elettrico per verificare che la simulazione dell'elec cardiaco sia stata eseguito come previsto. Assicurarsi quindi che il file di risultato heart-elec.odb si trova nella directory di lavoro.
    6. Passare alla seconda fase di simulazione passando al dominio MECH. Esaminare i valori delle costanti del materiale utilizzate nella simulazione meccanica per modellare la risposta cardiaca passiva e attiva desiderata.
    7. Assicurate che i file della libreria dei materiali per l'analisi meccanica utilizzino il nome della stringa HYBRID. Per modificare la risposta del materiale delle camere cardiache, regolate il file di materiale ibrido appropriato o sostituite l'intera risposta del materiale definendo un nuovo comportamento del materiale nella sezione Materiali del modulo CAE.
      NOTA: Informazioni dettagliate sulle leggi costitutive integrate sono disponibili nella guida per l'utente18.
    8. Nella fase PRE-LOAD, impostare le pressioni delle cavità idrostatiche per ottenere il comportamento fisiologico desiderato. Utilizzare l'opzione di ampiezza liscia integrata per salire da zero al livello di pressione desiderato, come descritto al passaggio 2.1.2.1.
    9. Disabilitare le condizioni limite di pressione definite al punto 2.1.2.1 per eseguire il modello del flusso sanguigno con un volume sanguigno complessivo costante all'interno del sistema di circolazione. Eseguire il processo di simulazione denominato heart-mech, come descritto nella sezione 2.1.2.5.

3. Stenosi della valvola aortica

NOTA: La stenosi aortica è spesso un driver di HFpEF in quanto porta a sovraccarico di pressione e, in definitiva, a rimodellamento concentrico e perdita di conformità della parete ventricolare sinistra. L'emodinamica osservata nella stenosi aortica spesso progredisce verso quelle osservate in HFpEF.

  1. Il modello a parametro raggruppato
    1. Modificare il segnale di ingresso nell'elemento della sequenza ripetuta PS rispetto alla valvola aortica, situata nel compartimento ventricolare sinistro. Simulare una riduzione dell'area dell'orifizio pari al 70% rispetto alla linea di base (tabella S6).
      NOTA: I valori di ingresso rappresenteranno l'area orifizio della valvola stenotica durante ogni battito cardiaco. Il valore dell'area dell'orifizio può essere facilmente regolato moltiplicando il vettore dei valori di uscita iniziale dell'elemento PS della valvola aortica per un valore decimale corrispondente all'area dell'orifizio finale rispetto al suo valore originale. In questo lavoro, un fattore di 0,3 è stato utilizzato per raggiungere il 70% di costrizione.
  2. Il modello FEA
    1. Modificate la definizione di scambio fluido del parametro LINK-LV-ARTERIAL.
      NOTA: Questo parametro possiede un coefficiente di resistenza viscoso sintonizzato sul flusso sanguigno tra il ventricolo sinistro e l'aorta. L'area di scambio efficace può essere modificata per regolare il flusso sanguigno e creare l'appropriato modello di stenosi aortica (Tabella S7).
    2. Individuare la cartella della casella degli strumenti e copiare i file all'interno di tale cartella nella directory di lavoro principale.
    3. Eseguire una simulazione meccanica inversa eseguendo i file della casella deglistrumenti 18. A tal fine, modificare le pressioni di aspirazione del ventricolo sinistro e dell'atrio sinistro a 6 mmHg nella cavità del fluido per regolare il loro stato volumetrico iniziale per il modello di stenosi aortica. Eseguire la inversePreliminary.py di lavoro.
      NOTA: Il ricalcolo dello stato di sollecitazione zero utilizzando la simulazione meccanica inversa è necessario ogni volta che viene modificata la condizione limite.
    4. Una volta completata la simulazione meccanica inversa, eseguire le funzioni di post-elaborazione: calcNodeCoords.py e straight_mv_chordae.py. Utilizzate i valori di default per gli altri parametri di flusso ed eseguite una nuova simulazione meccanica come descritto nella sezione 2.1.2.5.

4. Emodinamica HFpEF

NOTA: Per simulare gli effetti del rimodellamento cronico, sono state modificate le proprietà meccaniche del cuore sinistro.

  1. Il modello a parametro raggruppato
    1. Modificare la conformità diastolica ventricolare sinistra dell'elemento di conformità LV per imitare l'irrigidimento della parete a causa del sovraccarico di pressione, utilizzando il valore della conformità end-diastolica nella tabella S8.
      NOTA: si supponga che la conformità sia drop linearmente da end-systole a end-diastole.
    2. Aumentare la resistenza alle perdite della pompa LV a 18 × 106 Pa s m-3 (Tabella S8)per catturare le elevate pressioni ventricolari sinistra osservate in HFpEF.
  2. Il modello FEA
    1. Modificate le proprietà del materiale attivo della geometria del ventricolo sinistro. Aumentare il componente di rigidità per ottimizzare la risposta attiva del tessuto che influisce sui componenti di sollecitazione nelle direzioni della fibra e del foglio nel modello costitutivo.
      1. Modificate la risposta materiale del ventricolo sinistro nel file mech-mat-LV_ACTIVE file.
        NOTA: L'entità della rigidità per la camera ventricolare sinistra può essere sintonizzata per fornire gli effetti di conformità diastolici appropriati.
      2. Aumentare i parametri di rigidità a e b nella formulazione iperelastica anisotropica per catturare la maggiore risposta di rigidità per la fisiologia HFpEF.
      3. Nella fase PRE-LOAD, impostare le pressioni della cavità fluida del ventricolo sinistro e dell'atrio sinistro a 20 mmHg.
      4. Eseguire una simulazione meccanica inversa per ottenere lo stato volumetrico del ventricolo sinistro e dell'atrio. Esportare le coordinate nodali dal file heart-mech-inverse.odb 18.
      5. Eseguire le funzioni di post-elaborazione: calcNodeCoords.py e straight_mv_chordae.py, come descritto nel passaggio 3.2.4. Individuare i nuovi file di input nodali nella directory di lavoro ed eseguire una nuova simulazione meccanica, come descritto al punto 2.1.2.5.

Access restricted. Please log in or start a trial to view this content.

Results

$$\rightleftharpoonup{xx}$$ $$\longleftharp{xx}$$, $$\longrightharp{xx}$$,

I risultati delle simulazioni di base sono illustrati nella figura 3. Ciò descrive le forme d'onda di pressione e volume del ventricolo sinistro e dell'aorta (Figura 3A) nonché il circuito fotovoltaico ventricolare sinistro (Figura 3B). I due modelli in silico mostrano un'emodinamica ventricolare aortica e sinistra simile, che si trovano all'interno dell'intervallo fisiologico. Piccole differenze nella risposta prevista dalle due pi...

Access restricted. Please log in or start a trial to view this content.

Discussion

$$\rightleftharpoonup{xx}$$ $$\longleftharp{xx}$$, $$\longrightharp{xx}$$,

Le piattaforme lumped-parameter e FEA proposte in questo lavoro hanno ricapitolazione l'emodinamica cardiovascolare in condizioni fisiologiche, sia nella fase acuta del sovraccarico di pressione indotto dalla stenosi che nell'HFpEF cronico. Catturando il ruolo che il sovraccarico di pressione svolge nelle fasi acute e croniche dello sviluppo HFpEF, i risultati di questi modelli sono in accordo con la letteratura clinica di HFpEF, inclusa l'insorgenza di un gradiente di pressione transaortica a causa della stenosi aortica...

Access restricted. Please log in or start a trial to view this content.

Disclosures

$$\rightleftharpoonup{xx}$$ $$\longleftharp{xx}$$, $$\longrightharp{xx}$$,

Non ci sono conflitti di interesse associati a questo lavoro.

Acknowledgements

$$\rightleftharpoonup{xx}$$ $$\longleftharp{xx}$$, $$\longrightharp{xx}$$,

Riconosciamo i finanziamenti del programma Harvard-Massachusetts Institute of Technology Health Sciences and Technology e del SITA Foundation Award dell'Institute for Medical Engineering and Science.

Access restricted. Please log in or start a trial to view this content.

Materials

List of materials used in this article
NameCompanyCatalog NumberComments
Abaqus SoftwareDassault Systè mes Simulia Corp.Versione utilizzata: 2018; software di simulazione FEA
HETVALDassault Systè mes Simulia Corp.Versione utilizzata: 2018
Libreria idraulica (isotermica)MathWorksVersione utilizzata: 2020a
Living Heart Human ModelDassault Systè mes Simulia Corp.Versione utilizzata: V2_1, piattaforma FEA anatomicamente accurata di cuore umano adulto a 4 camere
MATLABMathWorksVersione utilizzata: 2020a, solutore numerico orientato agli oggetti
SIMSCAPE FLUIDSMathWorks
UAMPDassault Systè mes Simulia Corp.Versione utilizzata: 2018
VUANISOHYPERDassault Systè mes Simulia Corp.Versione utilizzata: 2018

References

$$\rightleftharpoonup{xx}$$ $$\longleftharp{xx}$$, $$\longrightharp{xx}$$,
  1. Borlaug, B. A., Paulus, W. J. Heart failure with preserved ejection fraction: Pathophysiology, diagnosis, and treatment. European Heart Journal. 32 (6), 670-679 (2011).
  2. Borlaug, B. A., Kane, G. C., Melenovsky, V., Olson, T. P. Abnormal right ventricular-pulmonary artery coupling with exercise in heart failure with preserved ejection fraction. European Heart Journal. 37 (43), 3293-3302 (2016).
  3. Borlaug, B. A. Evaluation and management of heart failure with preserved ejection fraction. Nature Reviews Cardiology. 17 (9), 1-15 (2020).
  4. Carabello, B. A., Paulus, W. J. Aortic stenosis. The Lancet. 373 (9667), 956-966 (2009).
  5. Lam, C. S. P., Donal, E., Kraigher-Krainer, E., Vasan, R. S. Epidemiology and clinical course of heart failure with preserved ejection fraction. European Journal of Heart Failure. 13 (1), 18-28 (2011).
  6. Omote, K., et al. Left ventricular outflow tract velocity time integral in hospitalized heart failure with preserved ejection fraction. ESC Heart Failure. 7 (1), 167-175 (2020).
  7. Samson, R., Jaiswal, A., Ennezat, P. V., Cassidy, M., Jemtel, T. H. L. Clinical phenotypes in heart failure with preserved ejection fraction. Journal of the American Heart Association. 5 (1), (2016).
  8. Weber, K. T., Brilla, C. G., Janicki, J. S. Myocardial fibrosis: Functional significance and regulatory factors. Cardiovascular Research. 27 (3), 341-348 (1993).
  9. Borbély, A., et al. Cardiomyocyte stiffness in diastolic heart failure. Circulation. 111 (6), 774-781 (2005).
  10. Borlaug, B. A., Lam, C. S. P., Roger, V. L., Rodeheffer, R. J., Redfield, M. M. Contractility and Ventricular Systolic Stiffening in Hypertensive Heart Disease. Insights Into the Pathogenesis of Heart Failure With Preserved Ejection Fraction. Journal of the American College of Cardiology. 54 (5), 410-418 (2009).
  11. Penicka, M., et al. Heart Failure With Preserved Ejection Fraction in Outpatients With Unexplained Dyspnea. A Pressure-Volume Loop Analysis. Journal of the American College of Cardiology. 55 (16), 1701-1710 (2010).
  12. Owen, B., Bojdo, N., Jivkov, A., Keavney, B., Revell, A. Structural modelling of the cardiovascular system. Biomechanics and Modeling in Mechanobiology. 17 (5), 1217-1242 (2018).
  13. Zhou, S., et al. A review on low-dimensional physics-based models of systemic arteries: Application to estimation of central aortic pressure. BioMedical Engineering Online. 18 (1), 41(2019).
  14. Sagawa, K., Lie, R. K., Schaefer, J. Translation of Otto frank's paper "Die Grundform des arteriellen Pulses" zeitschrift für biologie 37. Journal of Molecular and Cellular Cardiology. 22 (1899), 253-254 (1990).
  15. Rosalia, L., Ozturk, C., Van Story, D., Horvath, M., Roche, E. T. Object-oriented lumped-parameter modeling of the cardiovascular system for physiological and pathophysiological conditions. Advanced theory and simulations. , (2021).
  16. Lopez-Perez, A., Sebastian, R., Ferrero, J. M. Three-dimensional cardiac computational modelling: METHODS, features and applications. BioMedical Engineering Online. 14, 35(2015).
  17. Xie, X., Zheng, M., Wen, D., Li, Y., Xie, S. A new CFD based non-invasive method for functional diagnosis of coronary stenosis. BioMedical Engineering Online. 17 (1), 36(2018).
  18. Abaqus Dassault, S. SIMULIA living heart human model user documentation. , (2017).
  19. Baillargeon, B., Rebelo, N., Fox, D. D., Taylor, R. L., Kuhl, E. The living heart project: A robust and integrative simulator for human heart function. European Journal of Mechanics, A/Solids. 48, 38-47 (2014).
  20. Moscato, F., et al. Use of continuous flow ventricular assist devices in patients with heart failure and a normal ejection fraction: a computer-simulation study. The Journal of Thoracic and Cardiovascular Surgery. 145 (5), 1352-1358 (2013).
  21. Fresiello, L., Meyns, B., Di Molfetta, A., Ferrari, G. A Model of the Cardiorespiratory Response to Aerobic Exercise in Healthy and Heart Failure Conditions. Frontiers in Physiology. 7 (189), (2016).
  22. Moscato, F., et al. Left ventricle afterload impedance control by an axial flow ventricular assist device: a potential tool for ventricular recovery. Artificial Organs. 34 (9), 736-744 (2010).
  23. Colacino, F. M., Moscato, F., Piedimonte, F., Arabia, M., Danieli, G. A. Left ventricle load impedance control by apical VAD can help heart recovery and patient perfusion: a numerical study. Asaio Journal. 53 (3), 263-277 (2007).
  24. Gu, K., et al. Lumped parameter model for heart failure with novel regulating mechanisms of peripheral resistance and vascular compliance. Asaio Journal. 58 (3), 223-231 (2012).
  25. Suga, H., Sagawa, K., Kostiuk, D. P. Controls of ventricular contractility assessed by pressure-volume ratio, Emax. Cardiovascular Research. 10 (5), 582-592 (1976).
  26. Fernandez de Canete, J., Saz-Orozco, P. d, Moreno-Boza, D., Duran-Venegas, E. Object-oriented modeling and simulation of the closed loop cardiovascular system by using SIMSCAPE. Computers in Biology and Medicine. 43 (4), 323-333 (2013).
  27. Heldt, T., Shim, E. B., Kamm, R. D., Mark, R. G., et al. Computational modeling of cardiovascular response to orthostatic stress. Journal of Applied Physiology. 92 (3), 1239-1254 (2002).
  28. Granegger, M., et al. A Valveless Pulsatile Pump for the Treatment of Heart Failure with Preserved Ejection Fraction: A Simulation Study. Cardiovascular Engineering and Technology. 10 (1), 69-79 (2019).
  29. Hay, I., Rich, J., Ferber, P., Burkhoff, D., Maurer, M. S. Role of impaired myocardial relaxation in the production of elevated left ventricular filling pressure. American Journal of Physiology-Heart and Circulatory Physiology. 288 (3), 1203-1208 (2005).
  30. Kadry, K., et al. Biomechanics of diastolic dysfunction: a one-dimensional computational modeling approach. American Journal of Physiology-Heart and Circulatory Physiology. 319 (4), 882-892 (2020).
  31. Luo, C., Ramachandran, D., Ware, D. L., Ma, T. S., Clark, J. W. Modeling left ventricular diastolic dysfunction: classification and key indicators. Theoretical Biology & Medical Modelling. 8, 14(2011).
  32. Burkhoff, D., et al. Left atrial decompression pump for severe heart failure with preserved ejection fraction: theoretical and clinical considerations. JACC: Heart Failure. 3 (4), 275-282 (2015).
  33. Ahmad Bakir, A., Al Abed, A., Stevens, M. C., Lovell, N. H., Dokos, S. A Multiphysics Biventricular Cardiac Model: Simulations With a Left-Ventricular Assist Device. Frontiers in Physiology. 9 (1259), (2018).
  34. Genet, M., Lee, L. C., Baillargeon, B., Guccione, J. M., Kuhl, E. Modeling pathologies of diastolic and systolic heart failure. Annals of Biomedical Engineering. 44 (1), 112-127 (2016).
  35. Sack, K. L., et al. Investigating the Role of Interventricular Interdependence in Development of Right Heart Dysfunction During LVAD Support: A Patient-Specific Methods-Based Approach. Frontiers in Physiology. 9 (520), (2018).
  36. Baillargeon, B., et al. Human cardiac function simulator for the optimal design of a novel annuloplasty ring with a sub-valvular element for correction of ischemic mitral regurgitation. Cardiovascular Engineering and Technology. 6 (2), 105-116 (2015).
  37. Sack, K. L., et al. Partial LVAD Restores Ventricular Outputs and Normalizes LV but not RV Stress Distributions in the Acutely Failing Heart in Silico. The International Journal of Artificial Organs. 39 (8), 421-430 (2016).
  38. Baumgartner, H., et al. Echocardiographic assessment of valve stenosis: EAE/ASE recommendations for clinical practice. Journal of the American Society of Echocardiography. 22 (1), 1-23 (2009).
  39. Rajani, R., Hancock, J., Chambers, J. The art of assessing aortic stenosis. Heart. 98, 14(2012).
  40. Vahanian, A., et al. Guidelines on the management of valvular heart disease: The Task Force on the Management of Valvular Heart Disease of the European Society of Cardiology. European Heart Journal. 28 (2), 230-268 (2007).
  41. Matiwala, S., Margulies, K. B. Mechanical approaches to alter remodeling. Current Heart Failure Reports. 1 (1), 14-18 (2004).
  42. NIH Clinical Trials Registry. ImCardia for DHF to Treat Diastolic Heart Failure (DHF) Patient a Pilot Study (ImCardia). , (2011).

Access restricted. Please log in or start a trial to view this content.

Reprints and Permissions

$$\rightleftharpoonup{xx}$$ $$\longleftharp{xx}$$, $$\longrightharp{xx}$$,

Request permission to reuse the text or figures of this JoVE article

Request Permission

Tags

Heart FailurePreserved Ejection FractionLumped Parameter ModelFinite Element ModelingCardiovascular HemodynamicsComputational MethodsMedical DevicesHydraulic ElementsWall ComplianceFluid CompressibilityContractilitySimulation SetupElectrical AnalysisSinoatrial NodeBoundary ConditionsCardiac Cycle