Um simples modelo empírico para quantificar a deformação residual de uma liga Cu-Al-Be-Nb com memória de forma

A simple empirical model to quantify the residual strain of a Cu-Al-Be-Nb shape memory alloy

RESUMO

O conhecimento amplo do comportamento termomecânico das ligas com memória de forma (LMF) faz-se necessário para obtenção de ligas com elevada histerese térmica, o que torna possível eliminar a necessidade de meios criogênicos em muitas aplicações. Desta forma, os modelos matemáticos surgem como uma alternativa eficaz para descrição deste comportamento termomecânico. No presente trabalho foi desenvolvido um modelo matemático capaz de mensurar a deformação residual para uma liga Cu-Al-Be-Nb em função das variáveis: deformação, temperatura e taxa de deformação. Observa-se que os resultados obtidos por meio do modelo apresentam boa concordância quando comparado com os dados experimentais. Além disso, verificou-se que o modelo é capaz de prever as condições necessárias para se obter a máxima deformação residual.

Palavras-chave Deformação residual; Taxa de deformação; Modelo empírico

ABSTRACT

The wide knowledge of the thermomechanical behavior of shape memory alloys (SMAs) is necessary to obtain alloys with high thermal hysteresis, which makes it possible to eliminate the need for cryogenic means in many applications. In this way, the mathematical models appear as an effective alternative to describe this thermomechanical behavior. In the present work, a mathematical model was developed able to measure the residual strain for a Cu-Al-Be-Nb alloy as a function of the variables: strain, temperature and strain rate. It is observed that the results obtained through the model are presented with good agreement when compared with the experimental data. Furthermore, it was found that the model is capable of predicting the conditions necessary to obtain the maximum residual strain.

Keywords Residual strain; Strain rate; Empirical model

1. INTRODUÇÃO

Os materiais inteligentes, cada vez mais, vêm estimulando uma série de pesquisas voltadas às diversas áreas do conhecimento. Em particular, tem as ligas com memória que apresentam a capacidade de retornarem a uma forma predefinida mediante um processo termomecânico adequado, após aparente deformação plástica.

As ligas com memória de forma (LMF) têm sido amplamente utilizadas em campos como engenharia civil, medicina, engenharia mecânica, dentre outras áreas. Também podem ser utilizadas em estruturas de aeronaves, dispositivos de tratamentos cirúrgicos, como músculos robóticos, bem como em outras áreas [¹].

Algumas aplicações das ligas com memória de forma exigem meios criogênicos para serem concretizadas [²]. Esta necessidade onera consideravelmente os custos de operação para consolidação da aplicação. Visando eliminar a necessidade de meios criogênicos, nas últimas décadas, tem se buscado LMF que apresentem elevada histerese térmica.

O mecanismo da pré-deformação tem sido empregado como uma alternativa eficaz para se conseguir elevar a histerese térmica das LMF [³, ⁴, ⁵]. Este procedimento baseia-se na estabilização da martensita induzida por tensão após o processo de pré-deformação. Esta martensita estabilizada (associada à uma parcela da deformação residual) é convertida termicamente em austenita com temperaturas superiores ao A_f, aumentando desta forma, a histerese térmica da liga.

O efeito da taxa de deformação e da temperatura, durante a pré-deformação, na deformação residual para a liga com memória de forma Cu-Al-Be-Nb foi investigado por Oliveira et al. [⁶]. Verificou-se que é possível conseguir elevação da deformação residual em temperaturas um pouco acima do M_s e com taxa de deformação, da ordem de 10^-5 s^-1.

Em qualquer que seja a aplicação é necessário que se tenha compreensão do comportamento da liga, para que seja possível otimizar o uso de seu potencial [⁷]. Consequentemente, se faz clara a necessidade da implementação de modelos numéricos que sejam capazes de simular, de modo preciso, o comportamento termomecânico dessas ligas.

Com o passar dos anos, muitos cientistas e engenheiros vêm contribuindo com a vasta literatura sobre o comportamento experimental das ligas com memória de forma [⁷]. Visando uma análise mais precisa do comportamento destes materiais, tem se tornado cada vez maior o interesse no desenvolvimento de modelos matemáticos capazes de descrevê-lo de maneira adequada, permitindo explorar todo o seu potencial.

Um esforço notável foi dedicado nas últimas décadas para o desenvolvimento de modelos matemáticos capazes de entender os mecanismos físicos subjacentes à transformação martensítica em LMF e facilitar o projeto e a fabricação de dispositivos baseados nestas ligas [⁸].

O comportamento termomecânico das ligas com memória de forma pode ser descrito por modelos constitutivos por meio de abordagens microscópicas ou macroscópicas [⁹]. Alternativamente, modelos empíricos ou fenomenológicos são relativamente muito simples e servem como desenvolvimento de dispositivos e sistemas inteligentes [¹⁰]. Os modelos empíricos são obtidos pela regressão de um número limitado de dados experimentais [¹¹].

A grande maioria dos modelos encontrados na literatura foi desenvolvida para condições de carregamento quase estático ignorando a geração de calor interno durante o carregamento [¹²]. Recentemente, tem havido um crescente interesse pelo desenvolvimento de modelos que incorporam os efeitos da taxa de carregamento no comportamento termomecânico das LMF [¹², ¹³, ¹⁴, ¹⁵].

Lei et al [¹⁶] propuseram uma equação exponencial para calcular um fator residual da memória de forma (associado a deformação residual) em função do número de ciclos, da temperatura e da taxa de deformação. A equação proposta possibilita prever a variação do fator residual para diferentes condições termomecânicas. Além disso, um modelo constitutivo macromecânico foi proposto para prever o comportamento mecânico cíclico a temperatura constante.

Um modelo micromecânico (baseado em um algoritmo de autômatos celulares e na teoria da plasticidade de Ramberg-Osgood) para prever acúmulo de deformação residual em função do número de ciclos foi proposto por Paradis et al., [¹⁷]. O modelo apresentou capacidade para reproduzir efeitos relacionados as LMF tais como: superelasticidade, ciclagem térmica e efeito memória de forma.

Observando a importância da avaliação de mecanismos que permitam elevação da histerese térmica das LMF, visando aplicações com baixo custo de operação, faz-se necessário o conhecimento amplo do comportamento termomecânico destas ligas, principalmente avaliar a influência da taxa de deformação. O desenvolvimento de um modelo matemático capaz de prever a influência da taxa de deformação na histerese térmica das LMF torna-se de fundamental importância para o projeto e fabricação de dispositivos LMF de baixo custo. Além disso, cabe destacar que, um modelo matemático desta natureza reduz os custos relativos a realização de ensaios termomecânicos que exigem meios criogênicos.

Neste sentido, o principal objetivo deste trabalho foi desenvolver um modelo empírico capaz de mensurar a deformação residual para uma liga Cu-Al-Be-Nb em função das variáveis: deformação, temperatura e taxa de deformação. A validação do modelo foi obtida comparando os resultados numéricos obtidos através do modelo com resultados experimentais encontrados em Oliveira et al [⁶].

2. CARACTERÍSTICAS DO MODELO

O modelo proposto é empírico e baseado nos dados experimentais obtidos por Oliveira et al [⁶] para uma liga Cu-Al-Be-Nb. Verificam-se primeiramente as relações de dependência individuais da deformação residual com a deformação aplicada, temperatura e taxa de deformação. Em seguida, foi proposto uma equação para quantificar a deformação residual que engloba os efeitos simultâneos da deformação, temperatura e taxa de deformação.

2.1 Relações de dependência:

Baseado nos dados experimentais observou-se que cada uma das grandezas: deformação (ε), temperatura (T) e taxa de deformação ($\dot{\epsilon }$) exercem influência sobre a deformação residual. Assim, a partir de algumas considerações matemáticas, é possível expressar a deformação residual em função dessas variáveis, ou seja:

$${\epsilon }_r=f\left(\epsilon ,\dot{\epsilon },T\right)$$ (1)

Para obter a expressão desejada, inicialmente, é necessário considerar a influência isolada de cada variável no processo. Consequentemente, é possível considerar que:

$${\epsilon }_r=f\left(\epsilon \right)$$ (2)

A Figura 1 mostra um diagrama esquemático que ilustra a influência da deformação aplicada na deformação residual observada nos dados experimentais.

Figura 1
Representação esquemática do comportamento da deformação residual com a deformação aplicada.

Considerando inicialmente a influência da deformação no comportamento da deformação residual, é possível reescrever a equação (2) da seguinte forma:

$${\epsilon }_r=A\epsilon +B$$ (3)

Por meio da análise dos dados experimentais obtidos por Oliveira et al [⁶] verificou-se que os coeficientes angulares (A) e linear (B) da equação (3) mostraram-se dependentes da taxa de deformação. Assim, será feito as seguintes considerações:

Se $\dot{\mathrm{\epsilon }}=\dot{\mathrm{\epsilon }}1={\dot{\mathrm{\epsilon }}}_{\mathrm{qs}}$, então ${\mathrm{\epsilon }}_{\mathrm{r}}={\mathrm{A}}_1\mathrm{\epsilon }+{\mathrm{B}}_1$. Mas, se $\dot{\mathrm{\epsilon }}={\dot{\mathrm{\epsilon }}}_2>{\dot{\mathrm{\epsilon }}}_{\mathrm{qs}}$, então ${\mathrm{\epsilon }}_{\mathrm{r}}={\mathrm{A}}_2\mathrm{\epsilon }+{\mathrm{B}}_2$. Desta forma, considerando que os parâmetros A e B variam com a taxa de deformação, é possível escrever que:

$\mathrm{A}=\mathrm{A}\left(\dot{\mathrm{\epsilon }}\right)$ e $\mathrm{B}=\mathrm{B}\left(\dot{\mathrm{\epsilon }}\right)$. Logo, a equação 3 pode ser reescrita da seguinte forma:

$${\epsilon }_r=A\left(\dot{\epsilon }\right)\epsilon +B\left(\dot{\epsilon }\right)$$ (4)

As expressões que estabelecem $\mathrm{A}=\mathrm{A}\left(\dot{\mathrm{\epsilon }}\right)$ e $\mathrm{B}=\mathrm{B}\left(\dot{\mathrm{\epsilon }}\right)$ são:

$$A\left(\dot{\epsilon }\right)\epsilon =A_1+\left(\frac{A_2-A_1}{{\dot{\epsilon }}_2-{\dot{\epsilon }}_{qs}}\right)\left(\dot{\epsilon }-{\dot{\epsilon }}_{qs}\right)$$ (5)

$$B\left(\dot{\epsilon }\right)\epsilon =B_1+\left(\frac{B_2-B_1}{{\dot{\epsilon }}_2-{\dot{\epsilon }}_{qs}}\right)\left(\dot{\epsilon }-{\dot{\epsilon }}_{qs}\right)$$ (6)

A Figura 2 mostra um diagrama esquemático que ilustra a influência da temperatura na deformação residual observada nos dados experimentais [⁶]. Os resultados indicam que a máxima deformação residual é obtida na condição de taxa de deformação quase estática e na temperatura A_s. Observa-se ainda que fixando a deformação imposta e a taxa de deformação, a deformação residual diminui de forma linear com o aumento da temperatura, conforme observado na Figura 2.

Figura 2
Representação esquemática do comportamento da deformação residual com a temperatura.

Acrescentando a influência da temperatura na equação (4)se tornará:

$${\epsilon }_r={\epsilon }_{r_1}-\frac{\partial {\epsilon }_r}{\partial T}{\left(T-A_S\right)}_{.}$$ (7)

onde A_s é temperatura inicial da transformação da martensita em austenita.

A inclinação $\frac{\partial {\epsilon }_r}{\partial T}$pode ser obtida da seguinte forma:

$$\frac{\partial {\epsilon }_r}{\partial T}=\frac{{\epsilon }_{r_1}-{\epsilon }_r}{A_S-T}$$ (8)

Porém $\frac{\partial {\epsilon }_r}{\partial T}$varia com ε, assim para ε = ε₁, considerando$\left(\frac{\partial {\epsilon }_r}{\partial T}\right)1=c_1$, e

para $$\left(\frac{\partial {\epsilon }_r}{\partial T}\right)2=c_2$$ . Portanto,

$${\epsilon }_r={\epsilon }_{r_1}-C\left(\epsilon \right){\left(T-A_S\right)}_{.}$$ (9)

Onde:

$$C=C\left(\epsilon \right)=c_1-\left(\frac{c2-c1}{\epsilon 1-\epsilon 2}\right)\left(\epsilon -{\epsilon }_1\right)$$ (10)

Assim, o modelo desejado é representado através da equação:

$${\epsilon }_r=A\left(\dot{\epsilon }\right)\epsilon +B\left(\dot{\epsilon }\right)-C\left(\epsilon \right)\left(T-A_S\right)$$ (11)

em que A e B são constantes dependentes da taxa de deformação e C é uma constante que depende da deformação aplicada.

De acordo com a equação (11) é possível prever a variação da deformação residual com as condições termomecânicas impostas ao material (temperatura, deformação aplicada e taxa de carregamento).

2.2 Comparação do modelo com os dados experimentais

Na Figura 3 apresenta a comparação entre os resultados obtidos pelo modelo proposto e os experimentos realizados por Oliveira et al [⁶], que utilizaram uma liga superelástica Cu-11,8Al-0,6Be-0,5Nb (% em peso) para avaliar a influência das taxas de deformação, temperatura e deformação aplicada nos valores de deformação residual. Tais deformações residuais foram obtidas por meio de ensaios de tração com ciclo de carga/descarga. As condições termomecânicas estudadas no trabalho de Oliveira et al [⁶] foram: deformações aplicadas de 4%, 5%, 6% e 7% em diferentes temperaturas, com taxas de deformação de 0,0005s^-1 e 0,00001s^-1 aplicadas durante a etapa de carregamento. Os resultados obtidos pelo modelo são bastante próximos ao resultado experimental, conforme observado na Figura 3. As maiores diferenças nos resultados, entre o modelo proposto e os dados experimentais, foram encontradas para as temperaturas de 268 e 323 K.

Figura 3
Comparação entre modelo proposto e os resultados experimentais: deformação residual x temperatura para deformação de 6% e taxa de deformação de 0,0005 s^-1

Comparações entre o modelo proposto e os dados experimentais, referentes à influência da deformação imposta, estão representadas nas Figuras 4, 5 e 6. Na Figura 4 a temperatura foi fixada em 298 K, a taxa de deformação em 0,0005s^-1 e a deformação imposta variou entre 4 e 6%. Para a Figura 5 a temperatura foi fixada em 243 K, a taxa de deformação em 0,0005s^-1 e a deformação imposta variou entre 4 e 7%. Na Figura 6 a temperatura foi fixada em 243 K, a taxa de deformação em 0,00001 s^-1 e a deformação imposta variou entre 4 e 7%. Por meio da análise comparativa entre as curvas dos dados experimentais e os dados fornecidos pelo modelo, é possível observar que há uma boa concordância.

Figura 4
Comparação entre modelo proposto e os resultados experimentais: deformação residual x deformação imposta para temperatura de 298 K e taxa de deformação de 0,0005 s^-1.

Figura 5
Comparação entre o modelo proposto e os resultados experimentais obtidos por Oliveira et al [⁶]: deformação residual x deformação imposta para temperatura de 243 K e taxa de deformação de 0,0005 s^-1.

Figura 6
Comparação entre modelo proposto e os resultados experimentais obtidos por Oliveira et al [⁶]: deformação residual x deformação imposta para temperatura de 243 K e taxa de deformação de 0,00001 s^-1.

Vale ressaltar, que o modelo foi capaz de prever com boa aproximação a deformação residual em função das condições termomecânicas imposta ao material. Apesar das condições termomecânicas para obtenção dos resultados das Figuras 4, 5 e 6 serem diferentes, é possível constatar que os incrementos de deformação residual (provocado pelo aumento da deformação aplicada) observados nos dados experimentais não apresentam diferenças consideráveis com os obtidos pelo modelo proposto. Destaca-se que a influência da deformação aplicada na deformação residual deve ser verificada fixando-se a taxa de deformação aplicada e a temperatura do ensaio.

Analisando as Figuras 5 e 6 é possível verificar o efeito da taxa de deformação na deformação residual. É possível observar nas Figuras 5 e 6 que para uma mesma deformação aplicada, o modelo é capaz de reproduzir o comportamento da deformação residual com a variação da taxa de deformação. Para uma mesma condição termomecânica aplicada, por exemplo para uma deformação 7% e uma temperatura de 243 K, as deformações residuais obtidas experimentalmente para as taxas de deformação de 0,0005 s^-1 e 0,00001 s^-1 foram de 3,5% e 4,1%, respectivamente. Já os valores destas deformações residuais previstas pelo modelo para a mesma condição termomecânica foram de 3,57% e 4,08% para as taxas de deformação de 0,0005 s^-1 e 0,00001 s^-1, respectivamente. Pode-se notar que o aumento da deformação residual, provocado pelo aumento da taxa de deformação, observado nos resultados experimentais é similar ao encontrado pelo modelo proposto.

3. CONCLUSÕES

Sobre o modelo empírico desenvolvido neste trabalho, é possível concluir que:

1. O modelo empírico proposto para quantificar a deformação residual é uma equação linear em função das variáveis: taxa de deformação, temperatura e deformação aplicada.

2. O modelo apresenta boa concordância com os dados experimentais.

3. O modelo permite prever as condições termomecânicas que maximizam à deformação residual.

AGRADECIMENTOS

Os autores agradecem a CAPES e ao CNPq pelo apoio financeiro.

BIBLIOGRAFIA

• ¹ JANI, J. M., LEARY, M., SUBIC A., et al, “A review of shape memory alloy research, applications and opportunities”, Materials & Design, v. 56, pp. 1078-1113, 2014.
• ² LI, Z., CAO, L. F., et al, “Deformation behavior of β1 phase in polycrystalline CuAlMnZnZr shape memory alloy with wide transformation hysteresis. Materials Science and Engineering: A”, v.403, n.1–2, pp.227-233, 2005.
• ³ PIAO, M., OTSUKA, K., MIYAZAKI, S., et al, “Mechanism of the A_S temperature increase by pre-deformation in thermoelastic alloys”, Materials Transactions JIM, v. 34, pp. 919-929, 1993.
• ⁴ HE, X. M., RONG, L. J., YAN, D. S., et al, “TiNiNb wide hysteresis shape memory alloy with low niobium content”, Materials Science and Engineering A, v. 371, pp. 193–197, 2004.
• ⁵ ZHAO, X., YAN, X., YANG, Y., et al, “Wide hysteresis NiTi(Nb) shape memory alloys with low Nb content (4.5 at.%)”, Materials Science and Engineering A, v. 438-440, pp. 575–578, 2006.
• ⁶ OLIVEIRA, D. F. D., et al, “Assessment of Pipe Coupling by Using the Recovery of Stress-Induced Martensites in Superelastic Cu-11.8Al-0.6Be-0.5Nb Alloy”, Journal of Materials Engineering and Performance, abril 2017.
• ⁷ BRINSON, L. C., “One-Dimensional Constitutive Behavior of Shape Memory Alloys: Thermomechanical Derivation with Non-Constant Material Functions and Redefined Martensite Internal Variable”, Journal of Intelligent Material Systems and Structures, Gottingen, pp. 229-242, abril 1993.
• ⁸ CISSE, C., ZAKI, W., ZINEB, T. B. “A review of constitutive models and modeling techniques for shape memory alloys”, International Journal of Plasticity, pp. 244-284, 2016.
• ⁹ PAIVA, A., SAVI, M.A. “An Overview of Constitutive Models for Shape Memory Alloys”, Mathematical Problems in Engineering, 2006.
• ¹⁰ KHANDELWAL, A., BURAVALLA, V. “Models for shape memory alloy behavior: An overview of modeling approaches”, Int J Struct Changes, v.1, pp. 111–148, 2009.
• ¹¹ OTSUKA, K. AND WAYMAN, C. M. (1999). Shape memory materials, Cambridge university press
• ¹² KADKHODAEI, M., RAJAPAKSE, R. K. N. D., MAHZOON, M., et al, “Modeling of the cyclic thermomechanical response of SMA wires at different strain rates”, Smart Materials and Structures, v. 16, n. 6, pp. 2091–2101, 2007.
• ¹³ ROH, J.-H. “Thermomechanical Modeling of Shape Memory Alloys with Rate Dependency on the Pseudoelastic Behavior”, Mathematical Problems in Engineering, p. 9, 2014.
• ¹⁴ AURICCHIO, F., FUGAZZA, D., DESROCHES, R. “Rate-dependent thermo-mechanical modelling of superelastic shape-memory alloys for seismic applications”, Journal of Intelligent Material Systems and Structures, v. 19, n. 1, pp. 47–61, 2008.
• ¹⁵ VITIELLO, A., GIORLEO, G., MORACE, R. E. “Analysis of thermomechanical behaviour of Nitinol wires with high strain rates”, Smart Materials and Structures, pp. 215-221, 2005.
• ¹⁶ LEI, H., ZHOU, B., WANG, Z., et al, “A macromechanical constitutive model of shape memory alloys under uniaxial cyclic loading”, Journal of Mechanics , v. 28, pp. 469–477, 2012.
• ¹⁷ PARADIS, A., TERRIAULT, P., BRAILOVSKI, V. “Modeling of residual strain accumulation of NiTi shape memory alloys under uniaxial cyclic loading”, Comput Mater Sci, v. 47, n.2, pp. 373–383, 2009.

Datas de Publicação

• Publicação nesta coleção
  24 Jul 2020
• Data do Fascículo
  2020

Histórico

• Recebido
  09 Jan 2019
• Aceito
  12 Jul 2019