2
Armação de Búzios – RJ – Brasil October 21 st to 25 th , 2018 Método dos Elementos Finitos e Software Livre no cálculo do campo gerado por Estimulação Magnética Transcraniana R. L. Q. Ribeiro * , M. A. F. Feitosa Universidade de Pernambuco, Recife, Brasil *e-mail: [email protected] Antecedentes, Motivação e Objetivos. A Estimulação Magnética Transcraniana (EMT) é uma técnica não-invasiva que possibilita o estímulo de partes do sistema nervoso e é utilizada em diversos tratamentos. Consiste na passagem de um pulso de corrente por uma bobina, para produzir um campo magnético com taxa de variação temporal alta o suficiente para provocar o surgimento de um potencial de ação no nervo. Em trabalhos anteriores, desenvolvidos por um dos autores, verificou-se, através de algoritmo numérico, que pulsos de corrente de amplitude insuficiente para obtenção de estímulos neuronais, quando truncados a uma pequena fração do tempo total de descarga, podem gerar estímulos efetivos. Essa descoberta permite o desenvolvimento de equipamentos mais compactos, eficientes e com potencial de funcionamento em altas taxas de repetição. O objetivo desse artigo é a análise da efetividade do campo elétrico gerado pelos pulsos truncados, num modelo realista para a cabeça humana, através do uso do método dos elementos finitos e utilizando softwares gratuitos. Métodos. Para o cálculo do campo gerado pelo pulso de corrente truncado percorrendo uma bobina de EMT, é necessário um modelo da geometria da cabeça considerando todas as suas camadas e condutividades, bem como um da bobina a ser utilizada. Para modelar a cabeça é utilizado o pacote de softwares FreeSurfer e o FSL BET/BETsurf, que segmentam o cérebro a partir de imagens de ressonância magnética do tipo T1 e T2, que possuem contrastes diferentes para tecidos diferentes. Com o auxílio do programa SimNIBS, o método dos elementos finitos nos permite calcular a intensidade do campo e a densidade de corrente num esquema de malhas volumétricas complexas. O SimNIBS foi desenvolvido por pesquisadores alemães e dinamarqueses, para sistema operacional GNU. As diferentes bobinas a serem utilizadas são modeladas pelo GNU Octave, uma linguagem computacional, desenvolvida para computação matemática. Foram modeladas e simuladas quatro bobinas diferentes. Na simulação é possível modificar configurações de distância da bobina para a cabeça e intensidade do pulso. Após escolhido um ponto na cabeça para centralizar a projeção, informa-se a condutividade de cada elemento e define-se como será considerada a anisotropia da cabeça. Na Figura 1 tem-se uma prévia do campo elétrico na matéria cinzenta, na tela do SimNIBS, onde também é exibida a janela de determinação das condutividades e os valores de densidade de corrente no escalpo (a), no crânio (b), na matéria branca (c), no liquido cérebro espinal e na matéria cinzenta (e), para uma bobina em forma de 8, de 70mm. Resultados. Quatro tipos de bobinas foram simulados, duas circulares e duas em forma de 8, com diâmetros diferentes. O pulso de corrente simulado teve dl/dt de 3x10 6 A/s. É possível comprovar através das simulações a preferência pelas bobinas formato de 8 sobre as circulares quando se compara a focalização de uma com a outra. Os valores (mínimo - máximo) do campo elétrico induzido e da densidade de corrente para as 4 bobinas variaram respectivamente de acordo com a Tabela 1. Calculou-se os resultados a partir de valores de condutividade definidos, pelo programa,

Método dos Elementos Finitos e Software Livre no cálculo ...media.cbeb.org.br/media/uploads/s/Metodo_dos_Elementos_Finitos_e... · em forma de 8, de 70mm. Resultados. Quatro tipos

  • Upload
    buitram

  • View
    220

  • Download
    0

Embed Size (px)

Citation preview

Page 1: Método dos Elementos Finitos e Software Livre no cálculo ...media.cbeb.org.br/media/uploads/s/Metodo_dos_Elementos_Finitos_e... · em forma de 8, de 70mm. Resultados. Quatro tipos

ArmaçãodeBúzios–RJ–Brasil

October21stto25th,2018

MétododosElementosFinitoseSoftwareLivrenocálculodocampogeradoporEstimulaçãoMagnéticaTranscraniana

R.L.Q.Ribeiro*,M.A.F.FeitosaUniversidadedePernambuco,Recife,Brasil

*e-mail:[email protected]

Antecedentes, Motivação e Objetivos. A Estimulação Magnética Transcraniana (EMT) é umatécnica não-invasiva que possibilita o estímulo de partes do sistema nervoso e é utilizada emdiversos tratamentos. Consiste na passagem de um pulso de corrente por uma bobina, paraproduzirumcampomagnéticocomtaxadevariação temporalaltao suficienteparaprovocarosurgimentodeumpotencialdeaçãononervo.Emtrabalhosanteriores,desenvolvidosporumdosautores, verificou-se, através de algoritmo numérico, que pulsos de corrente de amplitudeinsuficienteparaobtençãodeestímulosneuronais,quandotruncadosaumapequenafraçãodotempo total de descarga, podem gerar estímulos efetivos. Essa descoberta permite odesenvolvimentodeequipamentosmaiscompactos,eficientesecompotencialdefuncionamentoemaltastaxasderepetição.Oobjetivodesseartigoéaanálisedaefetividadedocampoelétricogeradopelos pulsos truncados, nummodelo realista para a cabeça humana, através do uso dométododoselementosfinitoseutilizandosoftwaresgratuitos.Métodos. Para o cálculo do campo gerado pelo pulso de corrente truncado percorrendo umabobina de EMT, é necessário ummodelo da geometria da cabeça considerando todas as suascamadas e condutividades, bem como um da bobina a ser utilizada. Para modelar a cabeça éutilizadoopacotedesoftwaresFreeSurfereoFSLBET/BETsurf,quesegmentamocérebroapartirde imagensde ressonânciamagnéticado tipoT1eT2,quepossuemcontrastesdiferentesparatecidos diferentes. Com o auxílio do programa SimNIBS, o método dos elementos finitos nospermite calcular a intensidade do campo e a densidade de corrente num esquema demalhasvolumétricascomplexas.OSimNIBSfoidesenvolvidoporpesquisadoresalemãesedinamarqueses,parasistemaoperacionalGNU.AsdiferentesbobinasaseremutilizadassãomodeladaspeloGNUOctave, uma linguagem computacional, desenvolvida para computação matemática. Forammodeladas e simuladas quatro bobinas diferentes. Na simulação é possível modificarconfiguraçõesdedistânciadabobinapara a cabeçae intensidadedopulso.Apósescolhidoumponto na cabeça para centralizar a projeção, informa-se a condutividade de cada elemento edefine-se como será considerada a anisotropia da cabeça. Na Figura 1 tem-se uma prévia docampo elétrico na matéria cinzenta, na tela do SimNIBS, onde também é exibida a janela dedeterminaçãodascondutividadeseosvaloresdedensidadedecorrentenoescalpo(a),nocrânio(b),namatériabranca(c),noliquidocérebroespinalenamatériacinzenta(e),paraumabobinaemformade8,de70mm.Resultados.Quatrotiposdebobinasforamsimulados,duascirculareseduasemformade8,comdiâmetrosdiferentes.Opulsodecorrentesimuladotevedl/dtde3x106A/s.Épossívelcomprovaratravésdas simulações apreferênciapelasbobinas formatode8 sobre as circularesquando secompara a focalização de uma com a outra. Os valores (mínimo - máximo) do campo elétricoinduzidoedadensidadedecorrenteparaas4bobinasvariaramrespectivamentedeacordocomaTabela1.Calculou-seosresultadosapartirdevaloresdecondutividadedefinidos,peloprograma,

Page 2: Método dos Elementos Finitos e Software Livre no cálculo ...media.cbeb.org.br/media/uploads/s/Metodo_dos_Elementos_Finitos_e... · em forma de 8, de 70mm. Resultados. Quatro tipos

ArmaçãodeBúzios–RJ–Brasil

October21stto25th,2018

emS/m.Essesvalores,apesardenãotermosasconfiguraçõesexatasutilizadasporeles,sãobempróximosdosencontradosporYoungSunChoetal(DOI:10.1109/IEMBS.2010.5626096).

DiscussãoeConclusões.Assimulaçõesrealizadasatéomomento,ecomparaçãocomaliteratura,comprovamaeficiênciadométodoedosprogramasescolhidos.Avantagemdeseremsoftwarelivre é permitirem amodificação não apenas dos parâmetros disponibilizados,mas tambémdageometria das bobinas, e o principal para essa pesquisa, da configuração do pulso de correnteaplicadonassimulações.Equaçõesquedescrevemospulsos truncadosdecorrenteestãosendodesenvolvidaseintegradasaossoftwaresnuméricosparasimulações.Comessemétodonuméricoespera-se mais uma constatação da eficácia desses pulsos em modelos detalhados da cabeçahumana.

FiguraseTabelas.Figura1:TeladoSimNIBSeosvaloresdedensidadedecorrentenoescalpo(a),

nocrânio(b),namatériabranca(c),noliquidocérebroespinalenamatériacinzenta(e).

Tabela1:Valores(min.–máx.)docampoelétricoinduzidoedadensidadedecorrente.

Bobinas Campo(V/m) Densidade(A/m²)

Figurade8a 8.36x10-4–15.1 2.07x10-4–10.9

Figurade8b 7.09x10-4–14.9 1.09x10-4–12.3

Circular70mm 2.11x10-4–8.5 0.15x10-4–6.35

Circular90mm 4.16x10-4–9.73 0.14x10-4–8

Palavras Chave. Estimulaçãomagnética transcraniana,Método dos elementos finitos, SoftwareLivre,Cálculodecampoelétrico.