Parâmetro-caroço e Modelagem de elemento finito de insuficiência cardíaca com fração de ejeção preservada

Resumo

Este trabalho introduz dois modelos computacionais de insuficiência cardíaca com fração de ejeção preservada baseada em uma abordagem de parâmetros e análise de elementos finitos. Esses modelos são utilizados para avaliar as alterações na hemodinâmica do ventrículo esquerdo e vasculatura relacionada induzida pela sobrecarga de pressão e diminuição da conformidade ventricular.

Resumo

Os esforços científicos no campo da modelagem computacional de doenças cardiovasculares têm se concentrado em grande parte na insuficiência cardíaca com redução da fração de ejeção (HFrEF), amplamente negligenciando a insuficiência cardíaca com a fração de ejeção preservada (HFpEF), que mais recentemente se tornou uma forma dominante de insuficiência cardíaca em todo o mundo. Motivados pela escassez de HFpEF em representações de silico, dois modelos computacionais distintos são apresentados neste artigo para simular a hemodinâmica do HFpEF resultante da sobrecarga de pressão ventricular esquerda. Primeiro, um modelo de parâmetros amontoados orientados a objeto foi desenvolvido usando um solucionador numérico. Este modelo é baseado em uma rede derivada de Windkessel em dimensões zero (0D), que depende das propriedades geométricas e mecânicas dos elementos constitutivos e oferece a vantagem de baixos custos computacionais. Em segundo lugar, foi utilizado um pacote de software de análise de elementos finitos (FEA) para a implementação de uma simulação multidimensional. O modelo FEA combina modelos multifísicos tridimensionais (3D) da resposta cardíaca eletromecânica, deformações estruturais e hemodinâmica à base de cavidade fluida e utiliza um modelo simplificado de parâmetros fixos para definir os perfis de troca de fluxo entre diferentes cavidades fluidas. Através de cada abordagem, tanto as alterações hemodinâmicas agudas e crônicas no ventrículo esquerdo quanto na vasculatura proximal resultantes da sobrecarga de pressão foram simuladas com sucesso. Especificamente, a sobrecarga de pressão foi modelada pela redução da área orifício da válvula aórtica, enquanto a remodelagem crônica foi simulada pela redução da conformidade da parede ventricular esquerda. Consistente com a literatura científica e clínica da HFpEF, os resultados de ambos os modelos mostram (i) uma elevação aguda do gradiente de pressão transaréltica entre o ventrículo esquerdo e a aorta e uma redução no volume do AVC e (ii) uma diminuição crônica no volume ventricular esquerdo terminante, indicativo de disfunção diastólica. Finalmente, o modelo FEA demonstra que o estresse no miocárdio HFpEF é notavelmente maior do que no tecido cardíaco saudável durante todo o ciclo cardíaco.

Introdução

A insuficiência cardíaca é uma das principais causas de morte em todo o mundo, que ocorre quando o coração é incapaz de bombear ou preencher adequadamente para acompanhar as demandas metabólicas do corpo. A fração de ejeção, ou seja, a quantidade relativa de sangue armazenada no ventrículo esquerdo que é ejetada a cada contração é usada clinicamente para classificar a insuficiência cardíaca em (i) insuficiência cardíaca com fração de ejeção reduzida (HFrEF) e (ii) insuficiência cardíaca com fração de ejeção preservada (HFpEF), para frações de ejeção inferiores ou superiores a 45%,^{respectivamente 1,2,3}. Os sintomas do HFpEF frequentemente se desenvolvem em resposta à sobrecarga de pressão ventricular esquerda, que pode ser causada por várias condições, incluindo estenose aórtica, hipertensão e obstrução do trato de saída ventricular esquerda^3,4,5,6,7. A sobrecarga de pressão impulsiona uma cascata de aberrações moleculares e celulares, levando ao espessamento da parede ventricular esquerda (remodelagem concêntrica) e, finalmente, ao enrijecimento da parede ou perda de conformidade^8,9,10. Essas alterações biomecânicas afetam profundamente a hemodinâmica cardiovascular, pois resultam em uma relação elevada de aumento do volume de pressão diastólica e na redução do volume diastólico final¹¹.

A modelagem computacional do sistema cardiovascular avançou na compreensão das pressões e fluxos sanguíneos tanto na fisiologia quanto na doença e tem fomentado o desenvolvimento de estratégias diagnósticas e terapêuticas¹². Em modelos silico são classificados em modelos de baixa ou alta dimensão, com o primeiro utilizando métodos analíticos para avaliar propriedades hemodinâmicas globais com baixa demanda computacional e este último fornecendo uma descrição multiescala e multifísica mais extensa da mecânica cardiovascular e hemodinâmica no domínio 2D ou 3D¹³. A representação de Windkessel com parâmetros é a mais comum entre as descrições de baixa dimensão. Com base na analogia do circuito elétrico (lei de Ohm), isso imita o comportamento hemodinâmico geral do sistema cardiovascular através de uma combinação de elementos resistivos, capacitivos e indutivos¹⁴. Um estudo recente deste grupo propôs um modelo alternativo windkessel no domínio hidráulico que permite a modelagem de mudanças na geometria e mecânica de grandes vasos-câmaras cardíacas e válvulas - de forma mais intuitiva do que os modelos analógicos elétricos tradicionais. Esta simulação é desenvolvida em um solucionador numérico orientado a objetos (ver a Tabela de Materiais) e pode capturar a hemodinâmica normal, efeitos fisiológicos do acoplamento cardiorrespiratório, fluxo sanguíneo respiratório na fisiologia uni-cardíaca e alterações hemodinâmicas devido à constrição aórtica. Esta descrição expande-se sobre as capacidades de modelos de parâmetros agrupados, oferecendo uma abordagem fisicamente intuitiva para modelar um espectro de condições patológicas, incluindo insuficiência cardíaca¹⁵.

Modelos de alta dimensão são baseados na FEA para calcular hemodinâmica esposteral e interações de estrutura fluida. Essas representações podem fornecer descrições detalhadas e precisas do campo de fluxo sanguíneo local; no entanto, devido à sua baixa eficiência computacional, eles não são adequados para estudos de toda a árvore cardiovascular^16,17. Um pacote de software (ver a Tabela de Materiais) foi empregado como uma plataforma FEA anatomicamente precisa do coração humano adulto de 4 câmaras, que integra a resposta eletromecânica, deformações estruturais e hemodinâmica baseada em cavidade fluida. O modelo de coração humano adaptado também compreende um modelo simples de parâmetros amontoados que define a troca de fluxo entre as diferentes cavidades fluidas, bem como uma caracterização mecânica completa do tecido cardíaco^18,19.

Vários modelos de metametro e FEA de insuficiência cardíaca foram formulados para capturar anormalidades hemodinâmicas e avaliar estratégias terapêuticas, particularmente no contexto de dispositivos de assistência circulatória mecânica para HFrEF^{20,21,22,23,24}. Uma ampla gama de modelos de parâmetros 0D de várias complexidades capturou com sucesso a hemodinâmica do coração humano em condições fisiológicas e HFrEF através da otimização de dois ou três elementos sistemas analógicos elétricos Windkessel^20,21,23,24. A maioria dessas representações são modelos uni ou biventriculares baseados na formulação de tempo variando-elastância para reproduzir a ação contratil do coração e usar uma relação de pressão-volume não linear para descrever o preenchimento ventricular^25,26,27. Modelos abrangentes, que capturam a complexa rede cardiovascular e imitam a ação de bombeamento atrial e ventricular, têm sido usados como plataformas para testes de dispositivos. No entanto, embora exista um corpo significativo de literatura no campo da HFrEF, muito poucos em modelos silico de HFpEF foram propostos^{20,22,28,29,30,31}.

Um modelo de baixa dimensão da hemodinâmica HFpEF, recentemente desenvolvido por Burkhoff et al.³² e Granegger et al.^28,pode capturar os laços de volume de pressão (PV) do coração de 4 câmaras, recapitulando totalmente a hemodinâmica de vários fenótipos de HFpEF. Além disso, utilizam sua plataforma em silico para avaliar a viabilidade de um dispositivo circulatório mecânico para hfpef, pesquisa computacional pioneira do HFpEF para estudos de fisiologia, bem como desenvolvimento de dispositivos. No entanto, esses modelos permanecem incapazes de captar as mudanças dinâmicas nos fluxos sanguíneos e pressões observadas durante a progressão da doença. Um estudo recente de Kadry et al.³⁰ captura os vários fenótipos de disfunção diastólica ajustando o relaxamento ativo do miocárdio e a rigidez passiva do ventrículo esquerdo em um modelo de baixa dimensão. Seu trabalho fornece uma análise hemodinâmica abrangente da disfunção diastólica com base nas propriedades ativas e passivas do miocárdio. Da mesma forma, a literatura de modelos de alta dimensão tem focado principalmente em HFrEF^{19,33,34,35,36,37}. Bakir et al.³³ propuseram um modelo FEA de fluido cardíaco totalmente acoplado para prever o perfil hemodinâmico HFrEF e a eficácia de um dispositivo de assistência à esquerda-ventricular (LVAD). Este modelo biventricular (ou duas câmaras) utilizou um circuito windkessel acoplado para simular a hemodinâmica do coração saudável, HFrEF e HFrEF com suporte LVAD^33,37.

Da mesma forma, Sack et al.³⁵ desenvolveram um modelo biventricular para investigar disfunção ventricular direita. Sua geometria biventricular foi obtida a partir dos dados de ressonância magnética (RM) do paciente, e a malha de elemento finito do modelo foi construída usando segmentação de imagem para analisar a hemodinâmica de um ventrículo direito com falha suportada pelo VAD³⁵. Abordagens cardíacas FEA de quatro câmaras foram desenvolvidas para melhorar a precisão dos modelos do comportamento eletromecânico do coração^19,34. Em contraste com as descrições biventriculares, modelos de quatro câmaras derivadas da ressonância magnética do coração humano fornecem uma melhor representação da anatomia cardiovascular¹⁸. O modelo cardíaco empregado neste trabalho é um exemplo estabelecido de um modelo FEA de quatro câmaras. Ao contrário dos modelos FEA de parâmetros e biventriculares, essa representação captura alterações hemodinâmicas à medida que ocorrem durante a progressão da doença^34,37. Genet et al.³⁴, por exemplo, utilizaram a mesma plataforma para implementar um modelo de crescimento numérico da remodelagem observada no HFrEF e HFpEF. No entanto, esses modelos avaliam os efeitos da hipertrofia cardíaca apenas na mecânica estrutural e não fornecem uma descrição abrangente da hemodinâmica associada.

Para suprida a falta de HFpEF em modelos de silico neste trabalho, o modelo de parâmetros agrupados anteriormente desenvolvido por este grupo¹⁵ e o modelo FEA foram readaptados para simular o perfil hemodinâmico da HFpEF. Para isso, a capacidade de cada modelo de simular hemodinâmica cardiovascular na linha de base será demonstrada pela primeira vez. Os efeitos da sobrecarga de pressão ventricular esquerda induzida por estenose e da diminuição da conformidade ventricular esquerda devido à remodelagem cardíaca - uma marca típica do HFpEF - serão então avaliados.

Protocolo

1. Modelo de parâmetro amontoado de 0D

1. Configuração de simulação
   NOTA: No ambiente solucionador numérico (ver a Tabela de Materiais),construa o domínio conforme mostrado na Figura 1. Este é composto pelo coração de 4 câmaras, na parte superior do corpo, no corpo abdominal, na parte inferior do corpo e nos compartimentos torácicos, bem como da vasculatura proximal, incluindo a aorta, a artéria pulmonar e a cavae venae superior e inferior. Os elementos padrão utilizados nesta simulação fazem parte da biblioteca hidráulica padrão. Os detalhes podem ser encontrados nos Arquivos Suplementares.
   1. Navegue pela biblioteca hidráulica para encontrar os elementos necessários: tubulação hidráulica, câmara hidráulica de volume constante, resistência linear, bomba centrífuga, válvula de verificação, orifício de área variável e fluido hidráulico personalizado.
      1. Jogue os elementos do gasoduto hidráulico no espaço de trabalho.
         NOTA: Estes são responsáveis por perdas de atrito, bem como conformidade com a parede e compressão de fluidos em vasos sanguíneos e câmaras cardíacas. Através deste bloco, a perda de pressão é calculada utilizando-se a lei Darcy-Weisbach, enquanto a mudança de diâmetro devido à conformidade da parede depende da constante de proporcionalidade de conformidade, da pressão luminal e da constante de tempo. Finalmente, a compressão do fluido é definida pelo módulo a granel do meio.
      2. Insira os elementos de câmara hidráulica de volume constante para definir a conformidade da parede e a compressão dos fluidos.
         NOTA: Este bloco não leva em conta as perdas de pressão devido ao atrito.
      3. Adicione os elementos de resistência linear para definir a resistência ao fluxo.
         NOTA: Isto é independente das propriedades geométricas da vasculatura, análogamente ao elemento resistivo utilizado nos modelos elétricos analógicos windkessel. Outros blocos, como a bomba centrífuga, a válvula de verificação, o orifício de área variável e os elementos hidráulicos personalizados devem ser inseridos para gerar a entrada de pressão desejada ao sistema, modelar os efeitos das válvulas cardíacas no fluxo sanguíneo e definir as propriedades mecânicas do sangue. Através desses elementos, o comportamento do sistema cardiovascular tanto na fisiologia quanto na doença pode ser totalmente capturado. O sinal de entrada para a bomba centrífuga pode ser encontrado na Figura S1A.
      4. Modele a contratação de cada câmara cardíaca através do elemento de câmara de conformidade variáveis personalizada.
         NOTA: Isso aceita a conformidade como um sinal de entrada definido pelo usuário e é baseado no modelo de elastância de variação de tempo(Figura S1B-D).
   2. Fornecer os parâmetros relativos a cada elemento, como mostrado na Tabela S1,também encontrado em Rosalia et al.¹⁵
   3. Insira um elemento de sequência de repetição de sinal físico (PS) para cada um dos blocos que requerem um sinal de entrada definido pelo usuário: a bomba LV, os elementos de conformidade de conformidade variável e os blocos orifícios de área variável.
      NOTA: Os sinais de entrada utilizados para esta simulação podem ser encontrados na Figura S1.
   4. Selecione o solucionador implícito ODE 23t padrão e execute a simulação para 100 s para atingir um estado estável.

2. O modelo FEA

1. Configuração de simulação
   NOTA: O modelo FEA utiliza uma análise elétrica-mecânica acoplado em sequência. Neste modelo, a análise elétrica é realizada primeiro; em seguida, os potenciais elétricos resultantes são usados como fonte de excitação na análise mecânica a seguir. Portanto, a configuração de simulação contém dois domínios de trabalho: os domínios elétricos (ELEC) e mecânico(MECH),que são predefinidos no software de simulação FEA(Tabela de Materiais)¹⁸. Assim, a seção a seguir descreve apenas o fluxo de trabalho de análise. O modelo FEA utiliza as seguintes subrotinas do usuário HETVAL, VUANISOHYPERe UAMP para a modelagem de material elétrico e mecânico¹⁸.
   1. Navegue pelo domínio ELEC para realizar análises elétricas usando o procedimento de temperatura predefinido no módulo Padrão.
      1. Use uma etapa de análise única chamada BEAT. Defina a duração do ciclo cardíaco para 500 ms e aplique um pulso potencial elétrico a um conjunto de nó representando o nó sinoatrial (SA)(conjunto de nó: R_Atrium-1.SA_NODE).
      2. Revise a forma de onda elétrica padrão, que varia de -80 mV a 20 mV acima de 200 ms com a definição de amplitude de passo liso, conforme descrito no guia do modelo¹⁸. Use os valores padrão das constantes materiais na análise elétrica para ajustar o atraso av.
      3. Inicie o módulo Job e crie um trabalho chamado heart-elec.
   2. Uma vez concluída a configuração de análise elétrica, navegue pelo domínio MECH para realizar a análise mecânica baseada na cavidade do fluido.
      NOTA: A simulação mecânica é realizada após a análise elétrica, e os potenciais elétricos resultantes são usados como fonte de excitação para a análise mecânica. A análise mecânica contém múltiplas etapas.
      1. Use os três passos principais chamados PRE-LOAD, BEAT1e RECOVERY1. Na etapa PRÉ-LOAD, revise as condições de limite do estado pré-estressado do coração. Use 0,3 s como o tempo de passo para aumentar linearmente a pressão nas câmaras de fluidos.
         NOTA: Os valores predefinidos de pressão da cavidade do fluido são mostrados na tabela S3. O estado pré-estressado do coração já foi definido na configuração normal de simulação cardíaca, e as condições iniciais do nó são fornecidas nos arquivos de simulação externa, conforme listado na Tabela S5. O recálculo do estado de estresse zero utilizando a simulação mecânica inversa é necessário sempre que a condição de limite é modificada, conforme explicado nas etapas 3.2.2-3.2.4.
      2. Na etapa BEAT1, use 0,5 s como o tempo de passo para simular contração.
      3. Na etapa RECOVERY1, selecione 0,5 s para relaxamento cardíaco e preenchimento ventricular para uma frequência cardíaca de 60 bpm.
      4. Habilite as etapas subsequentes, BEATX e RECOVERYX,simular mais de um ciclo cardíaco para atingir um estado estável.
         NOTA: Três ciclos cardíacos serão suficientes para atingir o estado estável. Um ciclo da simulação é concluído em ~8 h em um processador de 24 núcleos (3,2 GHz x 24).
      5. Inicie o módulo Job e crie um trabalho chamado heart-mech, permitindo a opção de dupla precisão.
2. Revise o modelo windkessel de parâmetros fixos simplificados
   NOTA: O domínio mecânico do modelo FEA possui um modelo de fluxo sanguíneo, que é baseado em um circuito simplificado de parâmetros e é criado como uma combinação de cavidades de fluidos à base de superfície e trocas de fluidos, como visto na Figura 2¹⁸.
   1. Use a representação windkessel mencionada na nota acima para executar a simulação. Revise a representação do modelo de fluxo sanguíneo para ajustar os valores dos elementos resistivos e capacitivos para resistências de fluxo e conformidades estruturais, respectivamente.
   2. Revise a representação do elemento finito 3D de quatro câmaras cardíacas e garanta que suas posições geométricas sejam precisas.
   3. Verifique a montagem do coração e mude para o módulo Interação para ajustar os valores de conformidade e contrailidade de cada uma das quatro câmaras cardíacas.
      NOTA: Os valores padrão no módulo Interação são configurados para simular um ciclo de batimentos cardíacos humanos idealizado¹⁸.
   4. Revise as seguintes cavidades de fluido hidrostático no módulo de Interação, 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(Tabela S3).
   5. Utilize as câmaras de conformidade (CAV-ARTERIAL-COMP, CAV-PULMONARY-COMP e CAV-VENOUS-COMP) como volumes cúbicos, pois representam a conformidade das circulações arterial, venosa e pulmonar.
   6. Conecte três volumes cúbicos de conformidade a uma mola aterrada e revise o valor da rigidez para modelar a resposta ao volume de pressão nas circulações arterial, venosa e pulmonar.
   7. Verifique as seguintes definições de troca de fluidos entre as cavidades de fluidos hidrostáticos: átrio arterial-venoso, venoso-direito, ventrículo direito-direita, sistema ventrículo direito-pulmonar, átrio do sistema pulmonar esquerdo, ventrículo esquerdo-esquerdo e ventrículo esquerdo-aorta(Tabela S4).
   8. Ajuste o coeficiente de resistência viscoso para modificar o modelo de fluxo sanguíneo em cada link de troca de fluidos (consulte arquivos suplementares para obter mais informações sobre o efeito de resistência viscosa).
3. Simulação multifísica
   1. Localize o arquivo do banco de dados cae no diretório de trabalho.
      NOTA: O modelo FEA neste protocolo é entregue no banco de dados e é nomeado como LH-Human-Model-Beta-V2_1.cae.
   2. Insira os arquivos de entrada, objeto e biblioteca no diretório de trabalho para executar a simulação. Consulte tabela S5 para obter a lista completa de arquivos de entrada e biblioteca.
   3. Lançar o software de simulação de modelo FEA (ver a Tabela de Materiais).
      NOTA: Consulte o provedor de software para compatibilidade com versões^{posteriores 18}.
   4. Revise as condições de partes, montagem e limite nos domínios ELEC e MECH, conforme descrito nas seções 2.2 e 2.3.
   5. Primeiro, execute o trabalho de simulação elétrica chamado heart-elec,conforme descrito na seção 2.1.1.3. Inspecione visualmente os resultados potenciais elétricos para verificar se a simulação coração-elec funcionou como esperado. Em seguida, certifique-se de que o arquivo de resultado heart-elec.odb está no diretório de trabalho.
   6. Mova-se para a segunda fase de simulação mudando para o domínio MECH. Reveja os valores das constantes materiais utilizadas na simulação mecânica para modelar a resposta cardíaca passiva e ativa desejada.
   7. Certifique-se de que os arquivos da biblioteca de materiais para a análise mecânica usem o nome de sequência HÍBRIDA. Para modificar a resposta material das câmaras cardíacas, ajuste o arquivo de material híbrido apropriado ou substitua toda a resposta do material, definindo um novo comportamento material na seção Materiais no módulo CAE.
      NOTA: Informações detalhadas sobre as leis constitutivas incorporadas podem ser encontradas no guia do usuário¹⁸.
   8. Na etapa PRÉ-LOAD, defina as pressões das cavidades hidrostáticas para obter o comportamento fisiológico desejado. Use a opção de amplitude suave incorporada para aumentar de zero para o nível de pressão desejado, conforme descrito na etapa 2.1.2.1.
   9. Desativar as condições de limite de pressão definidas em 2.1.2.1 para executar o modelo de fluxo sanguíneo com um volume sanguíneo global constante dentro do sistema de circulação. Execute o trabalho de simulação chamado heart-mech, como descrito na seção 2.1.2.5.

3. Estenose da válvula aórtica

NOTA: A estenose aórtica é frequentemente um driver de HFpEF, pois leva à sobrecarga de pressão e, finalmente, à remodelação concêntrica e perda de conformidade da parede ventricular esquerda. A hemodinâmica observada na estenose aórtica frequentemente progride para aqueles vistos no HFpEF.

1. O modelo de parâmetros amontoados
   1. Modifique o sinal de entrada no elemento sequência repetindo PS em relação à válvula aórtica, localizada no compartimento ventricular esquerdo. Simule uma redução da área do orifício igual a 70% em relação à linha de base(Tabela S6).
      NOTA: Os valores de entrada representarão a área do orifício da válvula estenótica durante cada batimento cardíaco. O valor da área do orifício pode ser facilmente ajustado multiplicando o vetor de valores de saída inicial do elemento PS da válvula aórtica por um valor decimal correspondente à área de orifício final em relação ao seu valor original. Neste trabalho, utilizou-se um fator de 0,3 para atingir 70% de constrição.
2. O modelo FEA
   1. Modifique a definição de troca de fluidos do parâmetro LINK-LV-ARTERIAL.
      NOTA: Este parâmetro possui um coeficiente de resistência viscoso sintonizado no fluxo sanguíneo entre o ventrículo esquerdo e a aorta. A área de troca eficaz pode ser modificada para ajustar o fluxo sanguíneo e criar o modelo de estenose aórtica apropriada(Tabela S7).
   2. Localize a pasta caixa de ferramentas e copie os arquivos dentro dessa pasta para o diretório principal de funcionamento.
   3. Realize uma simulação mecânica inversa executando os arquivos da caixa de ferramentas¹⁸. Para isso, altere as pressões de sucção do ventrículo esquerdo e o átrio esquerdo para 6 mmHg na cavidade do fluido para ajustar seu estado volumoso inicial para o modelo de estenose aórtica. Execute a função inversePreliminary.py.
      NOTA: O recálculo do estado de estresse zero usando a simulação mecânica inversa é necessário sempre que a condição de limite é modificada.
   4. Uma vez concluída a simulação mecânica inversa, execute as funções de pós-processamento: calcNodeCoords.py e straight_mv_chordae.py. Use os valores padrão para os outros parâmetros de fluxo e execute uma nova simulação mecânica conforme descrito na seção 2.1.2.5.

4. Hemodinâmica HFpEF

NOTA: Para simular os efeitos da remodelagem crônica, as propriedades mecânicas do coração esquerdo foram modificadas.

1. O modelo de parâmetros amontoados
   1. Modifique a conformidade diastólica ventricular esquerda do elemento de conformidade LV para imitar o endurecimento da parede devido à sobrecarga de pressão, utilizando o valor da conformidade diastólica final na Tabela S8.
      NOTA: Assuma a conformidade para cair linearmente do sístole final para a estatamtole final.
   2. Aumente a resistência ao vazamento da bomba LV para 18 × 10⁶ Pa s m^-3 (Tabela S8) para capturar as pressões ventriculares esquerdas elevadas observadas no HFpEF.
2. O modelo FEA
   1. Edite as propriedades do material ativo da geometria do ventrículo esquerdo. Aumente o componente de rigidez para ajustar a resposta do tecido ativo que afeta os componentes de estresse nas instruções de fibras e folhas no modelo constitutivo.
      1. Modifique a resposta material do ventrículo esquerdo no arquivo mech-mat-LV_ACTIVE.
         NOTA: A magnitude da rigidez da câmara ventricular esquerda pode ser ajustada para fornecer os efeitos de conformidade diastólicas adequados.
      2. Aumente os parâmetros de rigidez a b na formulação hiperelástica anisotropica para capturar a resposta de rigidez aumentada para a fisiologia do HFpEF.
      3. Na etapa PRÉ-LOAD, defina as pressões de cavidade fluida do ventrículo esquerdo e o átrio esquerdo para 20 mmHg.
      4. Realize uma simulação mecânica inversa para obter o estado volumoso do ventrículo esquerdo e do átrio. Exporte as coordenadas nodal do arquivo heart-mech-inverse.odb ¹⁸.
      5. Executar as funções pós-processamento: calcNodeCoords.py e straight_mv_chordae.py, conforme descrito na etapa 3.2.4. Localize os novos arquivos de entrada nodal no diretório de trabalho e realize uma nova simulação mecânica, conforme descrito na seção 2.1.2.5.

Resultados

Os resultados das simulações da linha de base são ilustrados na Figura 3. Isso retrata as formas de onda de pressão e volume do ventrículo esquerdo e da aorta (Figura 3A),bem como o laço fotovoltaico esquerdo(Figura 3B). Os dois modelos de silico apresentam hemodinâmica aórtica e ventricular esquerda semelhante, que estão dentro da faixa fisiológica. Pequenas diferenças na resposta prevista pelas duas plataformas podem ser notadas d...

Discussão

As plataformas de parâmetros e FEA propostas neste trabalho recapitularam a hemodinâmica cardiovascular em condições fisiológicas, tanto na fase aguda da sobrecarga de pressão induzida por estenose quanto na HFpEF crônica. Ao captar o papel que a sobrecarga de pressão desempenha nas fases agudas e crônicas do desenvolvimento da HFpEF, os resultados desses modelos estão de acordo com a literatura clínica da HFpEF, incluindo o surgimento de um gradiente de pressão transatróptica devido à estenose aórtica, aumento da pressã...

Materiais

Lista de materiais utilizados neste artigo

Nome	Empresa	Número de catálogo	Comentários
Abaqus Software	Dassault Systè mes Simulia Corp.		Versão utilizada: 2018; Software de simulação de FEA
HETVAL	DASSAULT Systè mes Simulia Corp.		Versão usada: 2018
Biblioteca Hidráulica (Isotérmica)	MathWorks		Versão usada: 2020a
Living Heart Human Model	Dassault Systè mes Simulia Corp.		Versão usada: V2_1, plataforma FEA anatomicamente precisa de coração humano adulto de 4 câmaras
MATLAB	MathWorks		Versão usada: 2020a, solucionador numérico orientado a objetos
SIMSCAPE FLUIDS	MathWorks
UAMP	Dassault Systè mes Simulia Corp.		Versão usada: 2018
VUANISOHYPER	Dassault Systè mes Simulia Corp.		Versão usada: 2018