Comparação entre IDEA StatiCa Connection e ANSYS
Introdução
O objetivo desta comparação é demonstrar que os resultados fornecidos pelo IDEA StatiCa são semelhantes ou conservadores em relação a software de elementos finitos de uso geral de terceiros. Foi selecionado o software ANSYS [1], considerado comprovado e fiável. O modelo ANSYS foi criado e avaliado pelo prof. Ing. Jiří Kala, Ph.D. na Universidade de Tecnologia de Brno, República Checa.
Duas ligações selecionadas foram modeladas no ANSYS, versão 19.2, e no IDEA StatiCa Connection, versão 21.1.4.
Modelo ANSYS
Existem várias formas de modelar ligações de aço. O objetivo foi verificar que os elementos finitos utilizados no modelo IDEA StatiCa funcionam conforme previsto. Por isso, o modelo no ANSYS foi criado com a intenção de reproduzir o modelo gerado automaticamente no IDEA StatiCa. Embora isto seja relativamente simples para chapas, é muito complexo para soldaduras e parafusos, pois estes contêm elementos finitos especiais que simulam a resistência e as propriedades carga-deformação de parafusos e soldaduras de acordo com as normas, como o Code of Practice for the Structural Use of Steel 2011 [2]. Estão também ligados às chapas por uma montagem sofisticada de restrições multipontos e outros elementos complementares.
O modelo ANSYS foi criado utilizando elementos de casca SHELL 181 na sua linha de eixo. O SHELL 181 é um elemento de casca isoparamétrico de 4 nós com seis graus de liberdade em cada nó. Foram utilizadas cinco camadas de integração ao longo da espessura da casca. As chapas, as soldaduras, bem como as cabeças e porcas dos parafusos foram simuladas por este elemento com o critério de cedência de von Mises. As tensões de cedência das chapas foram de 275 MPa para chapas de aço com espessura inferior ou igual a 16 mm e de 265 MPa para chapas com espessura superior a 16 mm.
A simulação de ligações soldadas é uma tarefa complexa. Foi utilizado um modelo combinado de soldadura preparado por Turlier [3]. Este consiste num elemento de casca inclinado que simula a soldadura. A sua espessura e largura são iguais à espessura de garganta. Além disso, contém elementos de casca com modelo de material elástico que ligam o elemento de casca inclinado à malha dos elementos de casca que simulam as chapas ao longo da espessura da chapa. O critério de cedência de von Mises é tipicamente verificado para o elemento de casca inclinado da soldadura. Não é ideal para a comparação porque o modelo de cálculo da soldadura é simplificado e algumas tensões na soldadura não são consideradas.
Figura 1: Modelo combinado de soldadura
O contacto é descrito por um par de elementos CONTA 174 e TARGE 170 entre a superfície da placa de extremidade e o banzo da coluna, e também entre a cabeça do parafuso (porca) e as chapas. Foi utilizado o algoritmo de contacto pelo método de Lagrange aumentado e pesquisa de contacto nos pontos de Gauss. Foi utilizado um coeficiente de atrito de 0,3. A diferença para o contacto sem atrito é pequena. Como não é assumida pré-carga nos parafusos, a placa de extremidade move-se inicialmente para baixo, em pequenos incrementos de carga, sendo restringida pelos parafusos quando a tensão de contacto ainda é muito baixa. Foi utilizada a variante superfície-superfície do contacto flexível.
O parafuso foi modelado pelo elemento BEAM 188 com material elástico e área correspondente. O parafuso está fixo em ambas as extremidades nos elementos de casca que simulam a cabeça do parafuso e a porca. Elementos adicionais asseguraram o posicionamento do parafuso nos furos das chapas.
Foram criadas várias variantes com diferentes configurações de contacto. O contacto é, por natureza, uma característica geometricamente não linear. Foi encontrada uma solução para uma análise com grandes deformações, em que as equações de equilíbrio são atualizadas no modelo deformado; no entanto, foi também encontrada uma solução para análise com pequenas deformações. Foram testados contactos com coeficientes de atrito de 0,3 e 0,0. Estas variantes serviram para quantificar e minimizar o risco de imprecisões da análise numérica. As variantes acima mencionadas forneceram resultados coerentes e semelhantes. A avaliação detalhada foi realizada apenas com um modelo fiável que corresponde aos métodos do modelo IDEA StatiCa comparado.
Foi utilizado um solver direto de matriz esparsa para a análise. A análise não linear utiliza o método de Newton-Raphson completo. Foi utilizada a seleção automática de incrementos de carga. A carga inicial corresponde a 0,01 da carga aplicada; o incremento de carga mínimo e máximo é de 0,002 e 1, respetivamente. O número máximo de iterações em cada passo é 22.
Exemplo 1
O Exemplo 1 é uma ligação viga-coluna. A secção transversal da viga é UB 686 x 254 x 125. A secção transversal da coluna é UC 356 x 406 x 235 e está encastrada em ambas as extremidades. O aço S275 é utilizado para todos os elementos e chapas. A ligação é dimensionada como placa de extremidade com oito parafusos M45 de grau 10.9. A carga fatorada na viga na ligação é calculada:
- My = 920 kNm
- Vz = 460 kN
Modelo ANSYS
O modelo ANSYS tem uma coluna com um comprimento de 5 216 mm, o que corresponde ao modelo IDEA StatiCa. A coluna está encastrada em ambas as extremidades. A viga é modelada como uma consola com um comprimento de 2 000 mm (do nó à extremidade) e carregada pela força descendente de 460 kN, distribuída igualmente pelos nós que simulam a alma da viga. Ao contrário do modelo IDEA StatiCa, a coluna e a viga são criadas por elementos de casca em todo o comprimento. No IDEA StatiCa, apenas a ligação é criada por elementos de casca. Para o restante dos elementos, são utilizados os elementos condensados.
Foram utilizados 43 076 elementos SHELL 181 para criar o modelo. O modelo de análise tinha 259 326 equações com uma largura de banda de 144. Foram necessários 12 subpassos e 31 iterações para concluir a análise.
Figura 2: Modelo de casca ANSYS – vista geral
Figura 3: Pormenor da ligação
Figura 4: Pormenor das soldaduras e parafusos
Figura 5: Elementos de casca apresentados com a sua espessura – vista lateral
Figura 6: Elementos de casca apresentados com a sua espessura – axonometria
Comparação de resultados
Figura 7: Modelo ANSYS – tensão de von Mises – axonometria
Figura 8: Modelo IDEA StatiCa – tensão de von Mises – axonometria
Figura 9: Modelo ANSYS – tensão de von Mises – vista lateral
Figura 10: Modelo IDEA StatiCa – tensão de von Mises – vista lateral
Figura 11: Banzo traseiro da coluna – tensão de von Mises
Figura 12: Banzo frontal da coluna (na placa de extremidade) – tensão de von Mises
Figura 13: Alma da coluna – tensão de von Mises
Figura 14: Placa de extremidade – tensão de von Mises
Figura 15: Enrijecedores da coluna (chapas de continuidade)
A distribuição da tensão de von Mises em ambos os modelos é praticamente idêntica. As pequenas diferenças são atribuídas à malha mais refinada no modelo ANSYS e às diferenças na modelação de parafusos, soldaduras e contactos. Note-se que o IDEA StatiCa utiliza uma escala constante, enquanto a escala no ANSYS varia.
As tensões de pico são também muito semelhantes, como se pode observar nas tabelas seguintes. Uma diferença ligeiramente maior verifica-se na deformação plástica de pico da placa de extremidade. Tal deve-se novamente à malha mais refinada e às diferenças na modelação de parafusos e soldaduras.
Tabela 1: Tensões e deformações nas chapas – ANSYS
| Material | Espessura [mm] | \(\sigma\) [MPa] | \( \varepsilon_{pl} \) [-] | |
| C-bfl1 | S275 | 30.2 | 265 | 0.3 |
| C-tfl1 | S275 | 30.2 | 214 | |
| C-w1 | S275 | 18.4 | 265 | 0.1 |
| b-bfl1 | S275 | 16.2 | 265 | 0.07 |
| B-tfl1 | S275 | 16.2 | 265 | 0.05 |
| B-w1 | S275 | 11.7 | 275 | 0.01 |
| EP1 | S275 | 40 | 267 | 0.9 |
| STIFF1a | S275 | 18 | 201 | |
| STIFF1b | S275 | 18 | 201 | |
| STIFF1c | S275 | 18 | 118 | |
| STIFF1d | S275 | 18 | 118 |
Tabela 2: Tensões e deformações nas chapas – IDEA StatiCa
As diferenças são maiores no caso dos parafusos. No IDEA StatiCa, as forças nos parafusos nas duas filas superiores são sempre maiores. Tal deve-se à ocorrência de forças de alavanca. Isto é muito provavelmente causado pela maior rigidez dos parafusos no IDEA StatiCa à tração e por um contacto mais rígido. As forças de alavanca no IDEA StatiCa tendem a diminuir quando os parafusos estão mais carregados, cedem, deformam-se mais e as tensões nos contactos dissipam-se. O comportamento do perfil em T no IDEA StatiCa e a ocorrência de forças de alavanca são descritos, por exemplo, aqui. As diferenças entre as forças de corte podem ser atribuídas à diferença entre os contactos. O contacto no modelo ANSYS utiliza o coeficiente de atrito comummente usado de 0,3. Por outro lado, o IDEA StatiCa utiliza contacto sem atrito, que é a hipótese mais conservadora.
Tabela 3: Forças nos parafusos – ANSYS
| Força de tração [kN] | Força de corte [kN] | |
| B1 | 304 | 83 |
| B2 | 304 | 83 |
| B3 | 334 | 44 |
| B4 | 334 | 44 |
| B5 | 34.6 | 71 |
| B6 | 34.6 | 71 |
| B7 | 37.1 | 37 |
| B8 | 37.1 | 37 |
Tabela 4: Forças nos parafusos – IDEA StatiCa
As soldaduras são difíceis de avaliar no ANSYS devido à presença de tensões que são desprezadas no dimensionamento. No entanto, foi alcançada uma boa concordância entre o IDEA StatiCa e o ANSYS. Em geral, as tensões nas soldaduras mais significativas, como as soldaduras da viga à placa de extremidade, são ligeiramente superiores no IDEA StatiCa, o que significa que o seu dimensionamento é mais conservador. No caso de algumas soldaduras de enrijecedores, as tensões no ANSYS foram superiores.
Tabela 5: Tensões nas soldaduras – comparação entre ANSYS e IDEA StatiCa
| Item | Aresta | a [mm] | ANSYS fw [MPa] | IDEA StatiCa fw [MPa] |
| EP1 | B-bfl 1 | ◢10.0◣ | 202.1 | 217.6 |
| ◢10.0◣ | 207.5 | 218.4 | ||
| EP1 | B-tfl 1 | ◢10.0◣ | 214.1 | 217.5 |
| ◢10.0◣ | 196.4 | 216.6 | ||
| EP1 | B-w 1 | ◢6.0◣ | 215.1 | 218.2 |
| ◢6.0◣ | 215.1 | 218.2 | ||
| C-bfl 1 | STIFF1a | ◢8.0◣ | 106.3 | 144.6 |
| ◢8.0◣ | 206.2 | 190.6 | ||
| C-w 1 | STIFF1a | ◢8.0◣ | 201.1 | 68.6 |
| ◢8.0◣ | 61.0 | 65.9 | ||
| C-tfl 1 | STIFF1a | ◢8.0◣ | 90.4 | 76.3 |
| ◢8.0◣ | 65.1 | 60.8 | ||
| C-bfl 1 | STIFF1b | ◢8.0◣ | 195.1 | 191.8 |
| ◢8.0◣ | 129.2 | 145.5 | ||
| C-w 1 | STIFF1b | ◢8.0◣ | 207.1 | 65.9 |
| ◢8.0◣ | 63.6 | 68.7 | ||
| C-tfl 1 | STIFF1b | ◢8.0◣ | 110.0 | 60.8 |
| ◢8.0◣ | 86.5 | 76.3 | ||
| C-bfl 1 | STIFF1c | ◢8.0◣ | 157.5 | 162.2 |
| ◢8.0◣ | 135.2 | 158.1 | ||
| C-w 1 | STIFF1c | ◢8.0◣ | 29.4 | 67.6 |
| ◢8.0◣ | 28.2 | 65.8 | ||
| C-tfl 1 | STIFF1c | ◢8.0◣ | 54.4 | 51.8 |
| ◢8.0◣ | 74.4 | 66.5 | ||
| C-bfl 1 | STIFF1d | ◢8.0◣ | 137.6 | 159.8 |
| ◢8.0◣ | 161.1 | 163.7 | ||
| C-w 1 | STIFF1d | ◢8.0◣ | 87.9 | 65.8 |
| ◢8.0◣ | 92.4 | 67.6 | ||
| C-tfl 1 | STIFF1d | ◢8.0◣ | 65.4 | 66.5 |
| ◢8.0◣ | 54.2 | 51.8 |
Exemplo 2
O Exemplo 2 é uma ligação viga-coluna. A viga tem uma secção transversal UB 356 x 127 x 33. A coluna tem uma secção transversal UC 254 x 254 x 73 e está encastrada na base. Todo o aço utilizado é de grau S275. A ligação com placa de extremidade está equipada com seis parafusos M24 de grau 8.8. A carga fatorada na viga na ligação é calculada:
- My = 100 kNm
- Vz = 100 kN
Modelo ANSYS
O modelo ANSYS tem uma coluna com um comprimento de 1 606 mm, o que corresponde ao modelo IDEA StatiCa. A coluna está encastrada na base. A viga é modelada como uma consola com um comprimento de 1 000 mm (do nó à extremidade) e carregada por uma força descendente de 100 kN distribuída uniformemente pelos nós que simulam a alma da viga. Ao contrário do modelo IDEA StatiCa, a coluna e a viga são criadas por elementos de casca em todo o comprimento.
Foram utilizados 5 036 elementos SHELL 181 para criar o modelo. Tal resultou em 25 152 equações com uma largura de banda de matriz de 126. Para concluir a análise, foram necessários 11 subpassos e 22 iterações.
Figura 16: Modelo ANSYS – axonometria
Figura 17: Modelo ANSYS – pormenor na ligação
Figura 18: Modelo ANSYS – com espessura dos elementos de casca
Figura 19: Modelo ANSYS – vista lateral com espessura dos elementos de casca
Figura 20: Modelo ANSYS – vista lateral com espessura dos elementos de casca – pormenor da ligação
Comparação de resultados
Figura 21: ANSYS – Axonometria – tensão de von Mises
Figura 22: IDEA StatiCa – Axonometria – tensão de von Mises
Figura 23: Banzo da coluna na placa de extremidade – tensão de von Mises
Figura 24: Banzo da coluna na placa de extremidade – deformação plástica
Figura 25: Alma da coluna – tensão de von Mises
Figura 26: Enrijecedores da coluna – tensão de von Mises
Figura 27: Placa de extremidade – tensão de von Mises
Figura 28: Banzos da viga – tensão de von Mises
Figura 29: Alma da viga – tensão de von Mises
Figura 30: Alma da viga – deformações plásticas
A distribuição da tensão de von Mises em ambos os modelos é praticamente idêntica. As pequenas diferenças são atribuídas à malha mais refinada no modelo ANSYS e às diferenças na modelação de parafusos, soldaduras e contactos. Note-se que o IDEA StatiCa utiliza uma escala constante, enquanto a escala no ANSYS varia.
As tensões de pico são também muito semelhantes, como se pode observar nas tabelas seguintes. Uma diferença ligeiramente maior verifica-se na deformação plástica de pico da placa de extremidade. Tal deve-se novamente à malha mais refinada e às diferenças na modelação de parafusos e soldaduras.
Tabela 6: Tensões e deformações nas chapas – ANSYS
| Material | Espessura [mm] | \(\sigma\) [MPa] | \(\varepsilon_{pl}\) [-] | |
| C-bfl1 | S275 | 14.2 | 174 | |
| C-tfl1 | S275 | 14.2 | 275 | 0.386 |
| C-w 1 | S275 | 8.6 | 275 | 0.026 |
| B-bfl 1 | S275 | 8.5 | 246 | |
| B-tfl1 | S275 | 8.5 | 260 | |
| B-w1 | S275 | 6 | 275 | 0.077 |
| EP2 | S275 | 20 | 264 | |
| Stiff1a | S275 | 10 | 155 | |
| Stiff1b | S275 | 10 | 155 | |
| Stiff1c | S275 | 10 | 264 | |
| Stiff1d | S275 | 10 | 264 |
Tabela 7: Tensões e deformações nas chapas – IDEA StatiCa
Tabela 8: Forças nos parafusos – ANSYS
| Força de tração | Força de corte | |
| B1 | 104.2 | 14.7 |
| B2 | 104.2 | 14.7 |
| B3 | 47.1 | 14.3 |
| B4 | 47.1 | 14.3 |
| B5 | 12.1 | 21 |
| B6 | 12.1 | 21 |
Tabela 9: Forças e verificações dos parafusos – IDEA StatiCa
As diferenças são maiores no caso dos parafusos. No IDEA StatiCa, as forças nos parafusos são sempre superiores, exceto na fila inferior de parafusos. Tal deve-se à ocorrência de forças de alavanca. Isto é muito provavelmente causado pela maior rigidez dos parafusos no IDEA StatiCa à tração e por um contacto mais rígido. As forças de alavanca no IDEA StatiCa tendem a diminuir quando os parafusos estão mais carregados, cedem, deformam-se mais e as tensões nos contactos dissipam-se. As diferenças entre as forças de corte podem ser atribuídas à diferença entre os contactos. O contacto no modelo ANSYS utiliza o coeficiente de atrito comummente usado de 0,3. Por outro lado, o IDEA StatiCa utiliza contacto sem atrito, que é a hipótese mais conservadora.
Figura 31: Tensões de contacto entre a placa de extremidade e o banzo da coluna no ANSYS
Figura 32: Tensões de contacto entre a placa de extremidade e o banzo da coluna no IDEA StatiCa
As soldaduras são difíceis de avaliar no ANSYS devido à presença de tensões que são desprezadas no dimensionamento. No entanto, foi alcançada uma boa concordância entre o IDEA StatiCa e o ANSYS. Em geral, as tensões no IDEA StatiCa são ligeiramente superiores, o que significa que o seu dimensionamento é mais conservador.
Tabela 10: Tensões nas soldaduras
| Elemento | Soldadura | a [mm] | ANSYS fw [MPa] | IDEA StatiCa fw [MPa] |
| EP2 | B-bfl 1 | ◢6.0◣ | 218.0 | 215.7 |
| ◢6.0◣ | 166.5 | 215.7 | ||
| EP2 | B-tfl 1 | ◢6.0◣ | 129.2 | 120.7 |
| ◢6.0◣ | 88.3 | 135.9 | ||
| EP2 | B-w 1 | ◢5.0◣ | 219.1 | 215.6 |
| ◢5.0◣ | 219.1 | 215.6 | ||
| C-bfl 1 | STIFF1a | ◢4.0◣ | 40.8 | 41.5 |
| ◢4.0◣ | 60.8 | 57.3 | ||
| C-w 1 | STIFF1a | ◢4.0◣ | 47.5 | 61.2 |
| ◢4.0◣ | 37.9 | 57.5 | ||
| C-tfl 1 | STIFF1a | ◢4.0◣ | 167.1 | 137.2 |
| ◢4.0◣ | 111.0 | 105.7 | ||
| C-bfl 1 | STIFF1b | ◢4.0◣ | 62.7 | 57.2 |
| ◢4.0◣ | 41.8 | 41.4 | ||
| C-w 1 | STIFF1b | ◢4.0◣ | 47.5 | 57.6 |
| ◢4.0◣ | 66.4 | 61.2 | ||
| C-tfl 1 | STIFF1b | ◢4.0◣ | 120.2 | 105.4 |
| ◢4.0◣ | 167.4 | 136.9 | ||
| C-bfl 1 | STIFF1c | ◢4.0◣ | 58.8 | 32.2 |
| ◢4.0◣ | 30.8 | 30.8 | ||
| C-w 1 | STIFF1c | ◢4.0◣ | 83.2 | 80.9 |
| ◢4.0◣ | 65.4 | 82.4 | ||
| C-tfl 1 | STIFF1c | ◢4.0◣ | 174.0 | 215.8 |
| ◢4.0◣ | 164.3 | 214.3 | ||
| C-bfl 1 | STIFF1d | ◢4.0◣ | 19.6 | 30.8 |
| ◢4.0◣ | 20.9 | 32.2 | ||
| C-w 1 | STIFF1d | ◢4.0◣ | 73.9 | 82.4 |
| ◢4.0◣ | 96.6 | 80.9 | ||
| C-tfl 1 | STIFF1d | ◢4.0◣ | 163.3 | 214.0 |
| ◢4.0◣ | 173.6 | 215.8 |
Resumo
Duas ligações viga-coluna foram dimensionadas no IDEA StatiCa e comparadas com o ANSYS. As ligações de aço podem ser modeladas de muitas formas. O objetivo não era comparar diferentes técnicas de modelação, mas verificar o modelo de análise do IDEA StatiCa. Por isso, foi utilizada uma técnica de modelação semelhante no ANSYS – elementos de casca para chapas e soldaduras, e elementos de viga para parafusos. A malha era mais densa no modelo ANSYS e não continha quaisquer elementos especiais, como restrições multipontos ou elementos com critérios de rotura baseados em normas, neste caso, o Código de Hong Kong. As diferenças entre os modelos ANSYS e IDEA StatiCa são atribuídas a estas diferenças de modelação. No entanto, as diferenças são muito pequenas; os padrões de tensão e deformação plástica são praticamente idênticos. A principal diferença está nas forças nos parafusos, onde o IDEA StatiCa fornece forças de tração superiores, ou seja, resultados mais conservadores do que o ANSYS. As tensões nas soldaduras são difíceis de determinar, ao contrário do IDEA StatiCa, onde são utilizados elementos finitos especiais conformes com os requisitos normativos de dimensionamento. Em geral, verificou-se uma boa concordância entre as tensões nas soldaduras. As tensões nas soldaduras foram ligeiramente superiores no IDEA StatiCa, o que significa que o dimensionamento é conservador.
Referências
[1] Ansys® Mechanical Enterprise, Release 19.2
[2] Hong Kong Buildings Department, Code of Practice for Structural Use of Steel 2011 (2021 Edition), disponível em https://www.bd.gov.hk/doc/en/resources/codes-and-references/code-and-design-manuals/SUOS2011.pdf
[3] Turlier D., Klein P., Bérard F. ¨Seam Sim¨ method for seam weld structural assessment within a global structure FEA. Proc. Int. Conf. IIW2010 Istanbul (Turkey). AWST 651-658, 2010.