Pode selecionar diferentes opções no programa para apresentar os resultados não suavizados ou suavizados. Estas encontram-se descritas no manual online do RFEM (em inglês). As opções de suavização que permitem avaliar as forças internas e tensões da superfície também são explicadas neste artigo técnico. O exemplo do artigo seguinte demonstra como é que as forças internas dos elementos finitos são utilizadas para determinar os valores durante a suavização. Além da posição de um elemento finito, a teoria das placas – Mindlin ou Kirchhoff – também tem um papel importante neste processo.
Exemplo: placa de acordo com Mindlin e Kirchhoff
É criada uma placa como modelo bidimensional com as dimensões 4 m ⋅ 3 m. A placa está apoiada de forma articulada nas linhas de borda mais curtas de ambos os lados. A rotação das linhas de borda mais compridas é restringida por um encastramento. Para evitar os efeitos da deformação transversal, o coeficiente de Poisson do material é definido como zero. Para comparar a determinação das forças internas de acordo com as teorias de flexão de placas de Kirchhoff e Mindlin, é selecionada uma espessura de placa de 80 cm. Assim sendo, é admitida uma relação d/L = 0,2 na zona limite entre as duas teorias de flexão (ver https://www.dlubal.com/pt/apoio-tecnico-e-formacao/apoio-tecnico/faq/003158).
Para simplificar, o tamanho da malha de EF é definido para 1 m. Assim, existem 4 ⋅ 3 elementos finitos.
Sobre a placa atua uma carga de p = 5 kN/m². O peso próprio não é aplicado no caso de carga.
Forças internas não suavizadas
Os resultados do núcleo computacional para as forças internas vx e mx são considerados para uma secção longitudinal no centro da placa. O diagrama das forças internas não suavizadas elemento a elemento pode ser representado pela opção de visualização "Não contínuo". De acordo com as teorias de flexão de Mindlin e Kirchhoff, estão disponíveis os diagramas apresentados na Figura 02.
Suavização padrão para nós interiores
Para nós de EF que se encontram dentro de uma superfície, a média aritmética é, inicialmente, calculada a partir dos valores nodais dos elementos finitos adjacentes. Para esta abordagem, os elementos devem estar numa superfície e no mesmo lado de uma possível linha interna. O nó não pode ser um nó definido pelo utilizador na superfície.
Para a posição de descontinuidade em x = 1,00 m, os valores são os seguintes internamente (o resultado deste passo não está disponível para a saída).
Mindlin:
Kirchhoff:
Suavização para nós de borda
Nas bordas de uma superfície (opção de suavização "Contínuo nas superfícies") ou de um modelo ("Totalmente contínuo"), não existem elementos finitos adjacentes que possam ser utilizados para a média dos valores nodais. Portanto, é utilizada uma abordagem diferente, realizada em dois passos.
No primeiro passo, os valores médios são calculados para os nós que não estão localizados na borda de uma superfície ou de um modelo. Os valores dos nós nas bordas são calculados de forma a manter os valores originais no centro dos elementos finitos. No segundo passo, são determinados os valores médios para os nós das bordas.
Para o nó de borda em x = 0,00 m, resultam os seguintes valores.
Mindlin:
Kirchhoff:
Forças de corte
Segundo Mindlin, as forças de corte são calculadas como primeira derivação da deformação.
Estes valores são determinados pelo núcleo computacional e utilizados diretamente. A suavização das forças de corte é então realizada conforme descrito acima, dependendo da posição do nó de EF na superfície.
Na teoria de flexão de acordo com Kirchhoff, as forças de corte são calculadas como terceira derivação da deformação.
Esta determinação do valor utilizando a terceira derivação diminui significativamente a precisão. Por essa razão, as forças de corte do núcleo computacional não são utilizadas, sendo determinadas a partir das derivações dos momentos de uma abordagem melhorada.
Os momentos mx, my e mxy das equações representam os valores suavizados que são determinados de acordo com os métodos descritos acima. Comparativamente ao núcleo computacional, estes permitem obter valores mais precisos para as forças de corte.
Momentos
Enquanto os valores de momento nos pontos de integração correspondem aos valores teóricos, a extrapolação resulta numa perda de precisão quando se utiliza a função paraboloide hiperbólica para a suavização. Um paraboloide hiperbólico apenas aproxima a distribuição dos momentos no elemento. Por isso, é utilizado um algoritmo melhorado que substitui a extrapolação por um método integral avançado de força de corte. Uma vez que a distribuição da força de corte no elemento de superfície é assumida como a forma de um paraboloide hiperbólico (superfície de segunda ordem), a integral desta superfície representa uma superfície de terceira ordem que apresenta o diagrama de momentos com maior precisão.
Esta abordagem corresponde às equações descritas acima para determinar as forças de corte através das derivações dos momentos. Os momentos são então determinados com as seguintes equações.
As constantes de integração mx,0 e my,0 são calculadas a partir das condições do valor no centro do elemento.
S representa a superfície do elemento.
Os valores das forças internas são determinados de acordo com os métodos descritos acima, as integrais por uma integração numérica. Os momentos assim determinados são depois suavizados de acordo com o método para nós interiores ou nós de borda.
Foram obtidos os seguintes valores neste exemplo segundo a teoria de flexão de Mindlin.
Primeiro elemento:
Segundo elemento:
Foram obtidos os seguintes com a teoria de flexão segundo Kirchhoff.
Primeiro elemento:
Segundo elemento:
No processo de suavização, os passos descritos acima são realizados gradualmente no programa. Os valores suavizados podem ser representados graficamente com a opção de visualização "Contínuo nas superfícies".
.png?mw=350&hash=b023d6c658e181cb7d69028c0f3994dedab96fc5)
.jpg?mw=350&hash=bd355674e865e37d5548c329b8f5f0153cefe276)
.png?mw=600&hash=49b6a289915d28aa461360f7308b092631b1446e)