Análise numérica de um provete para a caracterização do comportamento à fractura da madeira sob solicitações de modo misto I/III

 

A.M.J.P. Barreto^a, M.F.S.F. de Moura^b*, R.D.S.G. Campilho^b

 

^a Instituto de Engenharia Mecânica e Gestão Industrial, Campus da FEUP,

Rua Dr.Roberto Frias,400, 4200-465 Porto, Portugal

^b Departamento de Engenharia Mecânica e Gestão Industrial, Faculdade de Engenharia da Universidade do Porto, Rua Dr. Roberto Frias, s/n, 4200-465 Porto, Portugal

*mfmoura@fe.up.pt

 

RESUMO: Neste trabalho realizou-se o estudo numérico de um provete adequado para a caracterização à fractura da madeira sob solicitações de modo misto I/III. O provete tem uma geometria similar à utilizada no ensaio Edge Crack Torsion e permite a obtenção de diferentes combinações de modo misto I/III. A distribuição das componentes de libertação de energia na frente de fenda foi estimada recorrendo ao método do fecho virtual de fenda. A variação da distância entre os actuadores permite de uma forma versátil alterar a combinação de modos pretendida. Recorreu-se ainda a um modelo coesivo de modo misto com o objectivo de simular a iniciação e a propagação considerando três combinações diferentes de modo misto. A aplicação dos métodos do fecho virtual de fenda e da calibração da flexibilidade propiciou resultados consistentes, validando assim o ensaio proposto.

Palavras chave: Madeira, Fractura, Modo Misto I+III

 

ABSTRACT: In this work a numerical analysis of a specimen adequate for wood fracture characterization under mixed-mode I+III loading was performed. The specimen has a geometry similar to the one used in the Edge Crack Torsion test and provides different combinations of mixed-mode I+III. The distribution of strain energy release rates at the crack front was estimated through the virtual crack close technique. The variation of the distance between loading actuators allows changing the mixed-mode ratio easily. A cohesive mixed-mode damage model was used to simulate damage initiation and growth considering three different mode ratios. The application of the virtual crack close technique and compliance calibration methods originate consistent results which validate the proposed test.

Keywords: Wood, Fracture, Mixed-mode I+III

 

1. Introdução

A aplicação estrutural da madeira tem vindo a aumentar, sendo por isso fundamental o estudo aprofundado do comportamento mecânico do material garantindo uma maior segurança para aumentar a confiança dos projectistas. Neste contexto adquire especial relevância o estudo do comportamento da madeira à fractura. A madeira é um material ortotrópico cuja estrutura anatómica evidencia três planos de simetria: a direcção longitudinal (L) das fibras, a direcção radial (R) dos anéis de crescimento e a direcção tangencial (T) aos anéis de crescimento. Devido às aplicações típicas da madeira os sistemas de propagação RL e TL são os mais propensos à propagação de fendas. A sua caracterização à fractura requer a determinação das taxas críticas de libertação de energia em modo I (G_I), modo II (G_II) e modo III (G_III). Existem inúmeros estudos na bibliografia sobre o comportamento da madeira sob solicitações em modo I [1-4] e em modo II [5-7]. De um modo geral pode-se afirmar que o ensaio Double Cantilever Beam (DCB) para o modo I e os ensaios End-Notched Flexure (ENF) e End-Loaded Split (ELS) para o modo II são os mais utilizados e propiciam medições rigorosas das respectivas taxas críticas de libertação de energia. Todavia o estudo de materiais não homogéneos sob solicitações em modo III é escasso, devido às enormes dificuldades inerentes à obtenção de uma solicitação de modo III puro. Donaldson (1988) [8], Ling (1998) [9] e Martin (1989) [10] estudaram o provete Split Cantilever Beam (SCB) em materiais compósitos, com o objectivo de prever o seu comportamento à fractura em puro modo III. Todos concluíram que havia uma grande componente de modo II nos bordos do provete, o que impossibilitava a propagação da fenda em puro modo III. O provete SCB não era, portanto, o indicado para este tipo de análise. Alternativamente, o provete Edge Crack Torsion (ECT) idealizado por Lee (1993) [11] tem sido utilizado para o estudo da fractura interlaminar de compósitos sob solicitações de modo III. Li et al. (1997) [12] estudaram este ensaio usando vários provetes com diferentes comprimentos de fenda (a) para calibrar a flexibilidade, tendo verificado que no centro do provete existe apenas modo III, como se pretende. Ratcliffe (2004) [13] estudou provetes de compósitos pré-impregnados unidirecionais para determinar a viabilidade do teste ECT para o cálculo da taxa de libertação de energia em modo III. Os resultados indicaram que o crescimento da delaminagem em modo III se inicia no centro do provete e os valores medidos da taxa crítica de libertação de energia em modo III dependem em grande parte do comprimento da delaminagem. De Morais et al. (2009) [14] usaram o ensaio ECT para a caracterização à fractura de compósitos de carbono/epóxido. Os autores procuraram identificar a geometria ideal do provete e o respectivo método de tratamento de resultados. Os resultados demonstraram uma não-linearidade considerável antes do ponto de carga máxima o que impede uma definição rigorosa do ponto de iniciação. No entanto, as análises feitas às zonas de não-lineariedade mostram dois prováveis pontos de iniciação correspondendo a valores de G_IIIc entre 850 and 1100 J/m² para todas as geometrias testadas. de Moura et al. (2009) [15] realizaram uma análise numérica do ensaio ECT para caracterizar a tenacidade de compósitos laminados. Para colmatar algumas incoerências inerentes à aplicação da Teoria das Placas, foi desenvolvido um novo método de tratamento de resultados. Concluiu-se que o teste ECT pode ser útil na caracterização à fractura dos compósitos laminados especialmente para comprimentos de fenda de 18 mm que corresponde a 41% do comprimento de fenda normalizado, estando estes valores sempre dependentes da geometria e do material do provete.

Embora a caracterização à fractura em modos puros seja útil, as aplicações estruturais da madeira estão, de um modo geral, sob solicitações de modo misto. De facto, tendo em consideração a ortotropia do material e a diferença de rigidez e resistência nas diferentes direcções, a fractura raramente ocorre por rotura das fibras da madeira, uma vez que a sua resistência é muito superior à existente nas direcções radial e tangencial. Consequentemente, a fenda tende a propagar-se paralelamente à direcção das fibras, independentemente da sua orientação inicial induzindo assim condições de fractura em modo misto. É assim fundamental a caracterização do material sob solicitações de modo misto. As solicitações de modo misto (I/II) foram objecto de estudo por alguns autores. Jernkvist (2001) [16] recorreu aos ensaios DCB com braços assimétricos e Single Edge-Notched Tensile, cujos provetes foram entalhados com uma fenda inclinada relativamente à aplicação da carga. Tchegg et al. (1994) [17] utilizaram o ensaio Wedge Splitting (WS) assimétrico que permite a alteração do rácio de modo misto I/II por modificação do ângulo das cunhas assimétricas. Oliveira et al. (2007) [18] realizaram um estudo numérico sobre o ensaio Mixed-Mode Bending (MMB) aplicado à fractura da madeira em modo misto I/II. Os mesmos autores concluíram que este ensaio é aplicável e apresenta a vantagem de propiciar uma fácil variação do rácio de modo misto. Estes mesmos autores realizaram ensaios experimentais e análises numéricas sobre os provetes Single Leg Bending (SLB) e End Loaded Split (ELS) e concluíram que estes ensaios, embora só providenciem um valor constante do rácio de modo misto I/II, apresentam uma vantagem importante relacionada com a simplicidade da sua execução.

Apenas existe um estudo de Ehart et al. (1998) [19] sobre solicitações de modo misto (I/III). Os seus autores realizaram testes de solicitações de modo misto (I/III) em materiais derivados da madeira, utilizando uma máquina bi-axial. Foi usado um provete designado por CT (Compact Tension), o qual se baseia numa geometria com dois entalhes laterais que diminuem a presença do modo II espúrio. Contudo, este ensaio, além de requerer uma máquina bi-axial, é de difícil execução e não elimina completamente a presença de modo II espúrio.

Neste trabalho foi idealizado e estudado numericamente um novo ensaio designado por Modified Edge Crack Torsion (MECT) que, como o próprio nome indica, se baseia no ensaio ECT. Recorrendo ao método Virtual Crack Closure Technique (VCCT) obtiveram-se os perfis de distribuição das três componentes da taxa de libertação de energia (G_I, G_II, G_III) na extremidade da fenda. Verificou-se que a solução proposta permite a variação do rácio de modo misto com facilidade, bastando para tal alterar a posição dos actuadores de carregamento. Procedeu-se ainda à simulação numérica por elementos finitos da iniciação e propagação de dano usando modelos coesivos.

 

2. Ensaio MECT

A concepção inicial do provete Edge Crack Torsion (ECT) proposta por Lee (1993) [11] consiste numa placa, com uma fenda no seu plano médio, submetida a torção que origina uma solicitação predominante de modo III. Para isso, este provete está sujeito a uma solicitação combinada com três pinos fixos e um actuador, originando assim um carregamento de corte em modo III. Para além do modo III na zona central do provete existe modo II espúrio desde os pinos até aos bordos exteriores.

 

Figura 1 - Geometria do provete ECT usado por Silva et al. (2006).

 

Baseado neste ensaio (ECT), idealizou-se um novo ensaio MECT (Modified Edge Crack Torsion) em que os dois pinos posteriores se mantêm fixos tal como no ECT, substituindo os dois pinos anteriores por dois actuadores paralelepipédicos que induzem uma solicitação de abertura da pré-fenda (ver Figura 2). A alteração da combinação de modo misto I/III consegue-se pela alteração da distância entre os elementos actuadores (L´) e das suas dimensões.

 

Figura 2 - Provete MECT.

 

3. Combinação de modos I/III

O objectivo deste estudo consistia na definição de geometrias que permitissem a obtenção de diferentes combinações de modos, particularmente em três situações distintas: predominância de modo I, predominância de modo III e um rácio unitário, ou seja presença equitativa de modo I e modo III na frente de fenda. Neste contexto, foram feitas inúmeras simulações variando a distância entre os elementos actuadores (L´) com o objectivo de calcular as componentes da taxa crítica de libertação de energia em modo I e em modo III que propiciassem as combinações pretendidas. As restantes dimensões (Tabela 1) foram mantidas constantes. As simulações foram realizadas recorrendo ao software Abaqus^® usando elementos sólidos tridimensionais de oito nós e a elementos finitos de interface de oito nós previamente desenvolvidos por de Moura et al. (1997) e Gonçalves et al. (2000). A malha de elementos finitos usada neste trabalho tem 39236 elementos sólidos contínuos e 5580 elementos de interface, sendo particularmente refinada nas regiões da extremidade da fenda, dos punções actuadores e dos suportes. Na Figura 3 apresenta-se uma imagem da deformada obtida durante a solicitação.

 

Tabela 1- Dimensões em milímetros do provete MECT.

 

Figura 3 – Deformada do provete MECT com L´=12 mm.

 

Os perfis das taxas de libertação de energia na frente de fenda do provete MECT foram obtidos recorrendo ao método do fecho virtual de fenda (Virtual Crack Closure Technique (VCCT)). Nesta técnica, as tensões e os deslocamentos relativos dos nós dos elementos finitos de interface são utilizados para obter as componentes das taxas de libertação de energia (G_i, i=I, II e III)

 

 

Nestas equações s _j3, t _j31 e t _j32 representam as tensões na frente de fenda (Figura 4). Por sua vez, u_kt, v_kt and w_kt são os deslocamentos dos nós localizados na face superior imediatamente antes da extremidade da fenda (nós kt na Figura 4), enquanto que u_kb, v_kb e w_kb são os deslocamentos dos nós homólogos da face inferior (nós kb, na Figura 4). A principal vantagem desta técnica reside no facto das distribuições de G_i poderem ser obtidas numa única análise por elementos finitos.

 

Figura 4 - Esquema de localização dos nós locais, utilizado para aplicação do método VCCT.

 

Os valores das propriedades mecânicas da madeira de Pinus pinaster Ait. introduzidos nas simulações numéricas encontram-se nas Tabelas 2 e 3 (Silva et al. (2006)).

 

Tabela 2 - Propriedades elásticas da madeira de Pinus pinaster Ait. (Silva et al. (2006)).

 

Tabela 3 - Propriedades de resistência e fractura da madeira de Pinus pinaster Ait. (Silva et al. (2006)).

 

Calcularam-se então as áreas circunscritas pelas curvas respeitantes a G_I e a G_III de modo a avaliar-se a relação de modo misto (ver Figura 5). Os valores dessas áreas divididas pela energia total (G_T=G_I+G_II+G_III) correspondem à percentagem de G_I e G_III face a G_T. A Tabela 4 contém as diferentes percentagens em função das distâncias entre os elementos actuadores na fase inicial do carregamento, ou seja sem qualquer dano.

 

Figura 5 – Perfis de distribuição das taxas de libertação de energia na frente de fenda para L´=28 mm.

 

Tabela 4 - Diferentes percentagens de G_I, G_II e G_III para diferentes valores de L´.

 

Conclui-se que, qualquer que seja a distância entre os elementos actuadores, o efeito do modo II é praticamente nulo e portanto a simulação ocorre em modo misto (I/III), como se pretende. Note-se que os valores de L´ escolhidos propiciam combinações de modos que contemplam os requisitos iniciais, isto é, G_I>G_III, G_I≈G_III e G_III>G_I.

 

4. Modelo coesivo

Com o objectivo de simular a iniciação e a propagação da fenda no provete MECT recorreu-se a um modelo coesivo de modo misto previamente desenvolvido [20]. O modelo de modo puro (Figura 6) baseia-se numa relação constitutiva linear entre as tensões (s) e os deslocamentos relativos (d_r). O modelo de modo misto baseia-se no mesmo princípio mas contempla a existência dos três modos de solicitação. A iniciação do dano é simulada recorrendo a um critério quadrático de tensões

 

 

onde s_i  e s_u,i (i = I, II, III) representam, respectivamente, as tensões e as resistências locais nos três modos. A propagação do dano é simulada pelo critério energético linear

 

 

sendo G_i e G_ic (i = I, II, III) as taxas de libertação de energia e os respectivos valores críticos, respectivamente. A área do triângulo mais pequeno da Figura 6 representa a energia libertada em cada modo, enquanto que a área do triângulo maior corresponde à respectiva energia crítica de fractura. A descrição detalhada do método pode ser encontrada na referência [21].

 

Figura 6 – Modelo coesivo em modo puro e em modo misto.

 

5. Resultados

Numa primeira análise recorreu-se ao método do VCCT como método de tratamento de resultados para obtenção da taxa critica de libertação de energia em modo misto, usando o critério energético linear (equação (3)). Assim, utilizando a equação (3) e os resultados da Tabela 4 calcula-se uma estimativa de G_T em função da distância entre os elementos actuadores, considerando os respectivos valores das taxas críticas de libertação de energia (Tabela 3). Considerando, por exemplo uma distância entre os elementos actuadores de 12 mm e recorrendo à equação (3), desprezando a componente de modo II, e à Tabela 4 tem-se que

 

 

de onde se obtém G_T = 0,321 N/mm. Aplicando o mesmo procedimentos aos restantes casos obtêm-se os valores apresentados na Tabela 5.

 

Tabela 5 - Valores de G_T usando o método VCCT.

 

Este procedimento apresenta como principal vantagem a sua simplicidade de aplicação. Todavia, apresenta também uma desvantagem que pode adulterar os resultados. De facto, nada garante que o perfil de distribuição de energias na frente de fenda se mantenha constante para valores de carregamento próximos da propagação. O desenvolvimento de uma zona de processo de fractura pode eventualmente alterar o rácio entre modos, ou seja os valores da Tabela 4 podem não ser exactamente os mesmos aquando da iniciação. Neste contexto, usou-se um método alternativo baseado na equação de Irwin-Kies

 

 

Este método requer a calibração da flexibilidade em função de diferentes comprimentos de fenda iniciais (a₀), ou seja, é necessário o estabelecimento da relação C₀=f(a₀), o que experimentalmente requer o fabrico de vários provetes com diferentes a₀. Numericamente procedeu-se à simulação de provetes com cinco comprimentos de fenda inicial. Verificou-se que a relação C₀=f(a₀) é bem reproduzida através de uma aproximação linear (Figura 7).

 

Figura 7 - Aproximação linear à curva representativa da calibração da flexibilidade.

 

De seguida procedeu-se a uma análise de propagação de dano usando o modelo coesivo de modo misto já descrito. A figura 8 apresenta a curva força-deslocamento típica destas simulações. Verificou-se que a propagação ocorre no início da não linearidade da curva força-deslocamento com um valor de força próximo da força máxima. Um dos problemas intrínsecos a este provete relaciona-se com o facto de não ser evidente a definição da dimensão da frente de fenda onde ocorre dissipação de energia (parâmetro B na equação (5)). Este parâmetro foi estimado através dos perfis de distribuição de energia na frente de fenda (Figura 5), onde se verifica que a dissipação de energia na frente de fenda está confinada a uma zona restrita.

 

Figura 8 – Curva P-d obtida para L´ = 22 mm.

 

A Tabela 6 contém os valores de G_T para as três distâncias entre os elementos actuadores (L’) analisados pelo VCCT. Comparando com os resultados da Tabela 5 verifica-se que existem diferenças não desprezáveis, o que se atribui ao menor rigor do VCCT como já foi discutido.

 

Tabela 6 - Valores calculados pela equação de Irwin-Kies.

 

Conclusões

Neste trabalho realizou-se um estudo numérico sobre o provete MECT aplicado ao estudo do comportamento da madeira sob solicitações de modo misto I/III. Este provete permite uma fácil variação das combinações de modos, bastando para tal alterar a distância entre os actuadores. A presença de modo II espúrio é desprezável, o que constitui uma importante vantagem tendo em consideração o fim em vista. Aplicaram-se dois métodos diferentes de tratamento de resultados tendo-se constatado que o VCCT, embora de simples aplicação não propicia resultados satisfatórios. É assim necessário recorrer ao método baseado na calibração da flexibilidade que implica a realização de vários provetes com diferentes comprimentos de fenda inicial. É também necessária uma análise prévia por elementos finitos para a determinação da extensão da frente de fenda onde ocorre dissipação de energia (parâmetro B da equação de Irwin-Kies), recorrendo ao método VCCT. De qualquer modo, conclui-se que o ensaio proposto satisfaz os requisitos necessários para a caracterização da madeira à fractura sob solicitações de modo misto I+III.

 

Agradecimentos

Os autores agradecem à Fundação para a Ciência e Tecnologia por suportar o trabalho aqui apresentado, através do projecto de investigação PDTC/EME-PME/64839/2006.

 

Referências

[1] Reiterer, A., Sinn, G., Stanzl-Tschegg, S. E., Fracture characteristics of different wood species under mode I loading perpendicular to the grain. Materials Science and Engineering A332 29–36, 2002.        [ Links ]

[2] Vasic, S., Smith, I., Bridging crack model for fracture of spruce. Engineering Fracture Mechanics,69:745-760, 2002.

[3] Dourado, N., Morel, S., de Moura, M.F.S.F., Valentin, G., Morais, J., Comparison of fracture properties of two wood species through cohesive crack simulations. Composites Part A: 39, 415–27, 2007.

[4] de Moura, M.F.S.F., Morais, J.J.L., Dourado, N., A new data reduction scheme for mode I wood fracture characterization using the DCB test. Engineering Fracture Mechanics 75:3852–3865, 2008.

[5] Silva, M.A.L., de Moura, M.F.S.F., Morais, J.J.L., Numerical analysis of the ENF test for mode II wood fracture, Composites Part A: 37,1334–1344, 2006.

[6] de Moura, M.F.S.F., Silva, M.A.L., de Morais, A.B., Morais J.J.L.. Equivalent crack based mode II fracture characterization of wood. Engineering Fracture Mechanics,73, 978-993, 2006.

[7] Yoshihara, H., Mode II R-curve of wood measured by 4-ENF test. Engineering Fracture Mechanics 71:2065–77, 2004.

[8] Donaldson, S. L., Mode III interlaminar fracture characterization of composite materials. Composites Science & Technology, 32: 225-49, 1988.

[9] Ling, C. L., Characterization of delamination in advanced composite materials under Mode III loading conditions. MSc Thesis, The Faculty of the School of Engineering of the Air Force Institute of Technology, Air University, USAF, December 1998.

[10] Martin, R. H., Evaluation of the split cantilever beam for mode III delamination testing. NASA Technical Memorandum 101562, March 1989.

[11] Lee, S. M., An Edge Crack Torsion Method for Mode III Delamination Fracture Testing. Journal of Composites Technology & Research, 15, 193-201, 1993.

[12] Li, J., Lee, S. M., Lee, E. W., O´Brien, T. K., Evaluation of the edge crack torsion (ECT) test for mode III interlaminar fracture toughness of laminated composites. Journal of Composites Technology & Research, 9, 174-183, 1997.

[13] Ratcliffe, J.G., Characterization of the edge crack torsion (ECT) test for mode III fracture toughness measurement of laminated composites. NASA/TM-2004-213269.

[14] de Morais, A.B., Pereira, A.B., de Moura, M.F.S.F., Magalhães, A.G., Mode III interlaminar fracture of carbon/epoxy laminates using the edge crack torsion (ECT) test, Composites Science and Technology, 69:670–676, 2009.

[15] de Moura, M. F. S. F., Fernandez, M.V.C., de Morais, A.B., Campilho, R.D.S.G., Numerical analysis of the edge crack torsion test for mode III interlaminar fracture of composite laminates, Engineering Fracture Mechanics 76:469-478, 2009

[16] Jernkvist, L., Fracture of wood under mixed mode loading II. Experimental investigation of Picea abies. Engineering Fracture Mechanics 68:565-576, 2001.

[17]  Tschegg, E. K., Elser, M., Stanzl-Tschegg, S.E., Biaxial Fracture Tests on Concrete – Development and Experience. Cement & Concrete Composites 17:57-75, 1995.

[18] Oliveira, J.M.Q., de Moura, M.F.S.F., Silva, M.A.L., Morais, J.J.L., Numerical analysis of the MMB test for mixed-mode I/II wood fracture. Composites Science and Technology, 67:1764-1771, 2007.

[19] Ehart R.J.A., Stanzl-Tschegg S.E., Tschegg E.K., Crack free interaction and mixed mode fracture of wood composites during mode III loading. Engineering Fracture Mechanics 61:253-278, 1998.

[20] de Morais, A.B., de Moura, M.F.S.F., Assessment of initiation criteria used in interlaminar fracture tests of composites. Engineering Fracture Mechanics, 72:2615-2627, 2005.

[21] Campilho, R.D.S.G., de Moura, M.F.S.F., Domingues, J.J.M.S., Modelling single and double-lap repairs on composite materials. Composites Science and Technology, 65:1948-1958, 2005.