(39 intermediate revisions by 3 users not shown)
Line 1: Line 1:
<!-- metadata commented in wiki content
+
<!--
 
+
 
+
 
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
 
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
 
<big>'''NUMERICAL IMPLEMENTATION FOR CONCEPTION OF STRUT AND TIE MODELS IN REINFORCED CONCRETE STRUCTURES'''</big></div>
 
<big>'''NUMERICAL IMPLEMENTATION FOR CONCEPTION OF STRUT AND TIE MODELS IN REINFORCED CONCRETE STRUCTURES'''</big></div>
Line 20: Line 18:
 
e-mail: [mailto:arturladeira@gmail.com arturladeira@gmail.com], [mailto:amilton@ufop.edu.br amilton@ufop.edu.br]
 
e-mail: [mailto:arturladeira@gmail.com arturladeira@gmail.com], [mailto:amilton@ufop.edu.br amilton@ufop.edu.br]
 
-->
 
-->
 +
==Resumo==
 +
 +
O Modelo de Bielas e Tirantes pode ser uma excelente alternativa para o dimensionamento de elementos estruturais em concreto armado submetidos a estado plano de tensão e para regiões que apresentem descontinuidade de ordem geométrica ou estática, substituindo procedimentos empíricos por uma metodologia racional de projeto. Para tornar a concepção do modelo menos dependente da experiência do projetista, o presente artigo tem como objetivo aliar a técnica de otimização topológica ESO (''Evolutionary Structural Optimization'') ao Método dos Elementos Finitos para geração automática dos modelos de bielas e tirantes. O critério de evolução do método de otimização topológica adotado considera a eliminação de elementos menos solicitados em termos de tensão, a partir de uma análise elástico-linear. Nesse contexto, é possível obter soluções otimizadas de problemas complexos envolvendo o concreto estrutural. São apresentados três exemplos numéricos para comprovação e validação das formulações e técnicas implementadas,  cujos modelos de bielas e tirantes obtidos apresentam boa concordância em relação às respostas encontradas em trabalhos científicos precursores sobre o tema.
 +
 +
'''Palavras-chave''': Modelo de bielas e tirantes, método dos elementos finitos, otimização topológica, concreto armado
  
 
==Abstract ==
 
==Abstract ==
Line 25: Line 28:
 
The Strut-and-Tie Model can be na excellent alternative for the design of reinforced concrete structural elements submitted to plane stress state and for regions with geometric or static discontinuity, replacing empirical procedures with a rational design methodology. To make the design of the model less dependent on the designer's experience, this article aims to combine the topological optimization technique ESO (Evolutionary Structural Optimization) with the Finite Element Method for the automatic generation of strut-and-tie models. The evolution criterion of the adopted topological optimization method considers the elimination of less stressed elements in terms of stress, from an elastic-linear analysis. In this context, it is possible to obtain optimized solutions of complex problems involving structural concrete. Three numerical examples are presented to prove and validate the formulations and techniques implemented whose strut-and-tie models present a good agreement in relation to the answers found in precursor scientific works on the subject.
 
The Strut-and-Tie Model can be na excellent alternative for the design of reinforced concrete structural elements submitted to plane stress state and for regions with geometric or static discontinuity, replacing empirical procedures with a rational design methodology. To make the design of the model less dependent on the designer's experience, this article aims to combine the topological optimization technique ESO (Evolutionary Structural Optimization) with the Finite Element Method for the automatic generation of strut-and-tie models. The evolution criterion of the adopted topological optimization method considers the elimination of less stressed elements in terms of stress, from an elastic-linear analysis. In this context, it is possible to obtain optimized solutions of complex problems involving structural concrete. Three numerical examples are presented to prove and validate the formulations and techniques implemented whose strut-and-tie models present a good agreement in relation to the answers found in precursor scientific works on the subject.
  
'''Keywords''' '': Strut-and-Tie Model, Finite Element Method, Topological Optimization, Reinforced Concrete.''
+
'''Keywords''': Strut-and-tie model, finite element method, topological optimization, reinforced concrete
 
+
==Resumo==
+
 
+
O Modelo de Bielas e Tirantes pode ser uma excelente alternativa para o dimensionamento de elementos estruturais em concreto armado submetidos a estado plano de tensão e para regiões que apresentem descontinuidade de ordem geométrica ou estática, substituindo procedimentos empíricos por uma metodologia racional de projeto. Para tornar a concepção do modelo menos dependente da experiência do projetista, o presente artigo tem como objetivo aliar a técnica de otimização topológica ESO (''Evolutionary Structural Optimization'') ao Método dos Elementos Finitos para geração automática dos modelos de bielas e tirantes. O critério de evolução do método de otimização topológica adotado considera a eliminação de elementos menos solicitados em termos de tensão, a partir de uma análise elástico-linear. Nesse contexto, é possível obter soluções otimizadas de problemas complexos envolvendo o concreto estrutural. São apresentados três exemplos numéricos para comprovação e validação das formulações e técnicas implementadas,  cujos modelos de bielas e tirantes obtidos apresentam boa concordância em relação às respostas encontradas em trabalhos científicos precursores sobre o tema.
+
 
+
'''Palavras-chave''' '': Modelo de Bielas e Tirantes, Método dos Elementos Finitos, Otimização Topológica, Concreto Armado.''
+
  
==1 Introdução==
+
==1. Introdução==
  
 
Para regiões ou elementos estruturais nos quais a Hipótese de Bernoulli não descreva adequadamente o comportamento estrutural ou distribuição de tensões, pode-se recorrer a outras alternativas de dimensionamento, tais como o Método dos Elementos Finitos (MEF) e o Modelo de Bielas e Tirantes (MBT). Nessas regiões, denominadas na literatura de ‘Regiões D’ (''Descontinuity''), as tensões de cisalhamento são significativas e a distribuição de deformações é não linear. Como exemplo, citam-se elementos estruturais como vigas-parede, consolos, sapatas, nós de pórticos, blocos rígidos de fundação sobre estacas, furos em vigas e dentes Gerber.
 
Para regiões ou elementos estruturais nos quais a Hipótese de Bernoulli não descreva adequadamente o comportamento estrutural ou distribuição de tensões, pode-se recorrer a outras alternativas de dimensionamento, tais como o Método dos Elementos Finitos (MEF) e o Modelo de Bielas e Tirantes (MBT). Nessas regiões, denominadas na literatura de ‘Regiões D’ (''Descontinuity''), as tensões de cisalhamento são significativas e a distribuição de deformações é não linear. Como exemplo, citam-se elementos estruturais como vigas-parede, consolos, sapatas, nós de pórticos, blocos rígidos de fundação sobre estacas, furos em vigas e dentes Gerber.
Line 39: Line 36:
 
Estudos pioneiros envolvendo o MBT tiveram origem no início do século XX. Ritter e Morsch propuseram, a partir de resultados experimentais, a analogia clássica do modelo de treliças para o dimensionamento a cisalhamento de vigas fletidas de concreto armado. Segundo esse modelo, admite-se a substituição da viga original por uma treliça equivalente definida a partir da distribuição de tensões. As barras tracionadas representam campos de tensão de tração (tirantes), enquanto as barras comprimidas representam campos de tensão de compressão (bielas).
 
Estudos pioneiros envolvendo o MBT tiveram origem no início do século XX. Ritter e Morsch propuseram, a partir de resultados experimentais, a analogia clássica do modelo de treliças para o dimensionamento a cisalhamento de vigas fletidas de concreto armado. Segundo esse modelo, admite-se a substituição da viga original por uma treliça equivalente definida a partir da distribuição de tensões. As barras tracionadas representam campos de tensão de tração (tirantes), enquanto as barras comprimidas representam campos de tensão de compressão (bielas).
  
Foi a partir dos trabalhos desenvolvidos por Schlaich ''et al''. [1], entretanto, que o tema ganhou forte impulso. Além das vigas inicialmente analisadas, esses pesquisadores estenderam a aplicação do modelo de bielas e tirantes a outros tipos de elementos estruturais. Abordaram temas como procedimentos para definição das regiões com e sem descontinuidade, geração dos modelos de treliças no interior do contínuo de concreto, cálculo dos esforços internos, diretrizes para verificação das tensões nas bielas e regiões nodais e cálculo e detalhamento da armadura necessária.
+
Foi a partir dos trabalhos desenvolvidos por Schlaich et al. [1], entretanto, que o tema ganhou forte impulso. Além das vigas inicialmente analisadas, esses pesquisadores estenderam a aplicação do modelo de bielas e tirantes a outros tipos de elementos estruturais. Abordaram temas como procedimentos para definição das regiões com e sem descontinuidade, geração dos modelos de treliças no interior do contínuo de concreto, cálculo dos esforços internos, diretrizes para verificação das tensões nas bielas e regiões nodais e cálculo e detalhamento da armadura necessária.
  
 
Entretanto, a não unicidade do modelo topológico, torna a concepção dependente da experiência e da sensibilidade estrutural do projetista para representar o fluxo interno de tensões. Sendo a armadura calculada e distribuída conforme o modelo topológico definido para o elemento estrutural de concreto armado, a correta definição desse não somente gera economia como também está relacionada à segurança.
 
Entretanto, a não unicidade do modelo topológico, torna a concepção dependente da experiência e da sensibilidade estrutural do projetista para representar o fluxo interno de tensões. Sendo a armadura calculada e distribuída conforme o modelo topológico definido para o elemento estrutural de concreto armado, a correta definição desse não somente gera economia como também está relacionada à segurança.
Line 49: Line 46:
 
Para a realização do objetivo principal pode-se destacar os seguintes objetivos específicos: implementação de um elemento finito triangular de três nós para análise em estado plano de tensão; implementação de uma rotina de otimização topológica dentro de um Programa de Elementos Finitos. A escolha do elemento triangular de três nós para a análise numérica de estruturas em estado plano de tensões, deve-se ao fato de que esse elemento requer uma discretização do contínuo bastante detalhada, permitindo assim definir as regiões de compressão e tração do modelo bielas e tirantes com mais refinamento. Uma vez que, em algumas etapas das análises a técnica de evolução utilizada no processo de otimização topológica na estrutura plana de concreto consiste em eliminar o elemento o da malha de elementos finitos, por isso a exigência de uma malha bem refinada, em consequência um elemento mais pobre em termos das funções de forma.
 
Para a realização do objetivo principal pode-se destacar os seguintes objetivos específicos: implementação de um elemento finito triangular de três nós para análise em estado plano de tensão; implementação de uma rotina de otimização topológica dentro de um Programa de Elementos Finitos. A escolha do elemento triangular de três nós para a análise numérica de estruturas em estado plano de tensões, deve-se ao fato de que esse elemento requer uma discretização do contínuo bastante detalhada, permitindo assim definir as regiões de compressão e tração do modelo bielas e tirantes com mais refinamento. Uma vez que, em algumas etapas das análises a técnica de evolução utilizada no processo de otimização topológica na estrutura plana de concreto consiste em eliminar o elemento o da malha de elementos finitos, por isso a exigência de uma malha bem refinada, em consequência um elemento mais pobre em termos das funções de forma.
  
==2 Formulação do Elemento Finito Triangular Linear==
+
==2. Formulação do elemento finito triangular linear==
  
 
Para geração do modelo de bielas e tirantes considerando análise linear, foi implementado a técnica ESO que é uma alternativa ao rigor matemático de métodos clássicos de otimização. O algoritmo baseia-se em uma análise via MEF e, como o presente estudo trata de elementos estruturais de concreto armado submetidos a estado plano de tensões, foi necessária, portanto, a implementação de um elemento finito bidimensional, cuja formulações é descrita neste item.
 
Para geração do modelo de bielas e tirantes considerando análise linear, foi implementado a técnica ESO que é uma alternativa ao rigor matemático de métodos clássicos de otimização. O algoritmo baseia-se em uma análise via MEF e, como o presente estudo trata de elementos estruturais de concreto armado submetidos a estado plano de tensões, foi necessária, portanto, a implementação de um elemento finito bidimensional, cuja formulações é descrita neste item.
Line 55: Line 52:
 
O elemento finito implementado neste trabalho é o elemento triangular de três nós e dois graus de liberdade por nó, denominado  ''Constant Strain Triangle'' (CST), usado na simulação numérica pelo MEF, baseado em deslocamentos. No caso da análise linear, o material é considerado homogêneo, isotrópico e linear.
 
O elemento finito implementado neste trabalho é o elemento triangular de três nós e dois graus de liberdade por nó, denominado  ''Constant Strain Triangle'' (CST), usado na simulação numérica pelo MEF, baseado em deslocamentos. No caso da análise linear, o material é considerado homogêneo, isotrópico e linear.
  
Segundo o Princípio dos Trabalhos Virtuais (PTV), aplicando-se um campo de deformação virtual compatível ao elemento triangular tem-se: δWint = δWext, isto é, o trabalho virtual externo é igual ao trabalho virtual interno. O trabalho virtual interno, δWint, pode ser escrito como:
+
Segundo o Princípio dos Trabalhos Virtuais (PTV), aplicando-se um campo de deformação virtual compatível ao elemento triangular tem-se: <math display="inline">\delta W_{int}= \delta W_{ext}</math>, isto é, o trabalho virtual externo é igual ao trabalho virtual interno. O trabalho virtual interno, <math display="inline">\delta W_{int}</math>, pode ser escrito como:
  
 
{| class="formulaSCP" style="width: 100%;width: 100%;text-align: center;"  
 
{| class="formulaSCP" style="width: 100%;width: 100%;text-align: center;"  
Line 62: Line 59:
 
{| style="text-align: center; margin:auto;width: 100%;"
 
{| style="text-align: center; margin:auto;width: 100%;"
 
|-
 
|-
| <math display="inline">\delta W_{int}=\underset{V}{\int\!\int\!\int }{\sigma }_{ij}\delta {\epsilon }_{ij}dV</math>  
+
| <math display="inline">\delta W_{int}=\underset{V}{\displaystyle \int\!\int\!\int }{\sigma }_{ij}\delta {\epsilon }_{ij}dV</math>  
 
|}
 
|}
|  style="text-align: right;vertical-align: top;width: 5px;text-align: right;white-space: nowrap;"|(1)
+
|  style="text-align: right;vertical-align: center;width: 5px;text-align: right;white-space: nowrap;"|(1)
 
|}
 
|}
  
 
+
onde  <math display="inline">\delta </math> é o operador variacional, <math>{\sigma }_{ij}</math> é o estado tensional real em um ponto qualquer no elemento, e  <math display="inline">\delta {\epsilon }_{ij}</math> é o estado de deformação virtual em um ponto qualquer no elemento. A partir do tensor de deformação de Green-Lagrange e desprezando-se as tensões no plano de normal na direção <math>z</math>, chega-se à Eq. (2) para trabalho virtual interno, onde <math>t</math> é a espessura do elemento triangular
onde  <math display="inline">\delta </math> é o operador variacional, σij é o estado tensional real em um ponto qualquer no elemento, e  <math display="inline">\delta {\epsilon }_{ij}</math> é o estado de deformação virtual em um ponto qualquer no elemento. A partir do tensor de deformação de Green-Lagrange e desprezando-se as tensões no plano de normal na direção ''z'',chega-se à Eq. (2) para trabalho virtual interno, onde ''t'' é a espessura do elemento triangular.
+
  
 
{| class="formulaSCP" style="width: 100%;width: 100%;text-align: center;"  
 
{| class="formulaSCP" style="width: 100%;width: 100%;text-align: center;"  
Line 75: Line 71:
 
{| style="text-align: center; margin:auto;width: 100%;"
 
{| style="text-align: center; margin:auto;width: 100%;"
 
|-
 
|-
| <math>\delta W_{int}=t\underset{A}{\int\!\int }[\delta u_{,x}{\sigma }_x+</math><math>\delta v_{,y}{\sigma }_y+(\delta u_{,y}+\delta v_{,x}){\tau }_{xy}]dA.</math>  
+
| <math>\delta W_{int}=t\underset{A}{\int\!\int }\left[\delta u_{,x}{\sigma }_x+\delta v_{,y}{\sigma }_y+(\delta u_{,y}+\delta v_{,x}){\tau }_{xy}\right]dA.</math>  
 
|}
 
|}
|  style="vertical-align: top;width: 5px;text-align: right;white-space: nowrap;"|(2)
+
|  style="vertical-align: center;width: 5px;text-align: right;white-space: nowrap;"|(2)
 
|}
 
|}
 
  
 
Para o elemento triangular de três nós é usual adotar as funções de interpolação em coordenadas naturais, as quais podem ser escritas como:
 
Para o elemento triangular de três nós é usual adotar as funções de interpolação em coordenadas naturais, as quais podem ser escritas como:
  
{| style="width: 100%;"  
+
{| class="formulaSCP" style="width: 100%;width: 100%;text-align: center;"  
 
|-
 
|-
| style="text-align: center;vertical-align: top;width: 92%;"|<math>N_1={\xi }_1=\xi =\frac{\left(y_2-y_3\right)}{2A}x+</math><math>\frac{\left(x_3-x_2\right)}{2A}y+\frac{(y_3-y_2)x_2-(x_3-x_2)y_1}{2A}</math>
+
|  
|  style="vertical-align: top;"|(3a)
+
{| style="text-align: center; margin:auto;width: 100%;"
 
|-
 
|-
|  style="text-align: center;vertical-align: top;width: 92%;"|<math>N_2={\xi }_2=\eta =\frac{\left(y_3-y_1\right)}{2A}x+</math><math>\frac{\left(x_1-x_3\right)}{2A}y+\frac{(x_3-x_1)y_1-(y_3-y_1)x_1}{2A}</math>  
+
|  style="text-align: center;vertical-align: top;"|<math>N_1={\xi }_1=\xi =\frac{\left(y_2-y_3\right)}{2A}x+</math><math>\frac{\left(x_3-x_2\right)}{2A}y+\frac{(y_3-y_2)x_2-(x_3-x_2)y_1}{2A}</math>  
|  style="vertical-align: top;"|(3b)
+
|}
 +
|  style="vertical-align: center;width: 5px;text-align: right;white-space: nowrap;"|(3a)
 
|-
 
|-
|  style="text-align: center;vertical-align: top;width: 92%;"|<math>N_3={\xi }_3=1-\xi -\eta </math>  
+
|
|  style="vertical-align: top;"|(3c)
+
{| style="text-align: center; margin:auto;width: 100%;"
 +
|-
 +
|  style="text-align: center;vertical-align: top;"|<math>N_2={\xi }_2=\eta =\frac{\left(y_3-y_1\right)}{2A}x+</math><math>\frac{\left(x_1-x_3\right)}{2A}y+\frac{(x_3-x_1)y_1-(y_3-y_1)x_1}{2A}</math>
 +
|}
 +
|  style="vertical-align: center;width: 5px;text-align: right;white-space: nowrap;"|(3b)
 +
|-
 +
|
 +
{| style="text-align: center; margin:auto;width: 100%;"
 +
|-
 +
|  style="text-align: center;vertical-align: top;"|<math>N_3={\xi }_3=1-\xi -\eta </math>  
 +
|}
 +
|  style="vertical-align: center;width: 5px;text-align: right;white-space: nowrap;"|(3c)
 
|}
 
|}
  
 
+
onde <math>x_1</math>, <math>x_2</math>, <math>x_3</math>, <math>y_1</math>, <math>y_2</math>, <math>y_3</math> são as coordenadas cartesianas do elemento e <math>A</math> sua área.
onde ''x''<sub>1</sub>, ''x''<sub>2</sub>, ''x''<sub>3</sub>, ''y''<sub>1</sub>, ''y''<sub>2</sub>, ''y''<sub>3 </sub>são as coordenadas cartesianas do elemento e ''A'' sua área.
+
  
 
Definindo o vetor de deslocamentos nodais por  <math>q={\left\{\begin{array}{cccccc}
 
Definindo o vetor de deslocamentos nodais por  <math>q={\left\{\begin{array}{cccccc}
Line 102: Line 108:
 
\end{array}\right\}}^T</math> e representando as funções de interpolação em coordenadas naturais pelo vetor coluna  <math>{\Phi }^T=\left(\begin{array}{ccc}
 
\end{array}\right\}}^T</math> e representando as funções de interpolação em coordenadas naturais pelo vetor coluna  <math>{\Phi }^T=\left(\begin{array}{ccc}
 
N_1 & N_2 & N_3
 
N_1 & N_2 & N_3
\end{array}\right)</math> , definem-se as equações aproximadas dos deslocamentos associados aos deslocamentos nodais q que podem ser escritas na forma matricial conforme a Eq. (4).
+
\end{array}\right)</math>, definem-se as equações aproximadas dos deslocamentos associados aos deslocamentos nodais q que podem ser escritas na forma matricial conforme a Eq. (4)
  
 
{| class="formulaSCP" style="width: 100%;width: 100%;text-align: center;"  
 
{| class="formulaSCP" style="width: 100%;width: 100%;text-align: center;"  
Line 117: Line 123:
 
\end{array}\right]q\mbox{.}</math>  
 
\end{array}\right]q\mbox{.}</math>  
 
|}
 
|}
|  style="text-align: right;vertical-align: top;width: 5px;text-align: right;white-space: nowrap;"|(4)
+
|  style="text-align: right;vertical-align: center;width: 5px;text-align: right;white-space: nowrap;"|(4)
 
|}
 
|}
  
 
+
Na Eq. (4), <math>O</math> é um vetor coluna nulo com três termos. Sendo os deslocamentos <math>u</math> e <math>v</math> funções dos deslocamentos nodais, os seus variacionais podem ser escritos a partir da seguinte expressão: <math display="inline">\delta a=\delta q^T\left(\frac{\partial a}{\partial q}\right)</math> , sendo  <math display="inline">\textbf{q}</math> o vetor dos deslocamentos nodais,  <math display="inline">\partial </math> o operador diferencial e <math>a</math> substituído por <math>u</math> e <math>v</math> Substituindo esses variacionais na Eq. (2) chega-se ao trabalho virtual de um elemento triangular dado pela equação a seguir
Na Eq. (4), O é um vetor coluna nulo com três termos. Sendo os deslocamentos ''u'' e ''v'' funções dos deslocamentos nodais, os seus variacionais podem ser escritos a partir da seguinte expressão: <math display="inline">\delta a=\delta q^T\left(\frac{\partial a}{\partial q}\right)</math> , sendo  <math display="inline">\boldsymbol{q}</math> o vetor dos deslocamentos nodais,  <math display="inline">\partial </math> o operador diferencial e ''a'' substituído por ''u'' e ''v''. Substituindo esses variacionais na Eq. (2) chega-se ao trabalho virtual de um elemento triangular dado pela equação a seguir.
+
  
 
{| class="formulaSCP" style="width: 100%;width: 100%;text-align: center;"  
 
{| class="formulaSCP" style="width: 100%;width: 100%;text-align: center;"  
Line 128: Line 133:
 
{| style="text-align: center; margin:auto;width: 100%;"
 
{| style="text-align: center; margin:auto;width: 100%;"
 
|-
 
|-
| <math>\delta W_{int}=\delta \boldsymbol{q}^Tt\underset{A}{\int\!\int }[\frac{\partial u_{,x}}{\partial \boldsymbol{q}}{\sigma }_x+</math><math>\frac{\partial v_{,y}}{\partial \boldsymbol{q}}{\sigma }_y+</math><math>(\frac{\partial u_{,y}}{\partial \boldsymbol{q}}+\frac{\partial v_{,x}}{\partial \boldsymbol{q}}){\tau }_{xy}]dA.</math>  
+
| <math>\delta W_{int}=\delta \textbf{q}^Tt\underset{A}{\int\!\int }\left[\frac{\partial u_{,x}}{\partial \textbf{q}}{\sigma }_x+\frac{\partial v_{,y}}{\partial \textbf{q}}{\sigma }_y+\left(\frac{\partial u_{,y}}{\partial \textbf{q}}+\frac{\partial v_{,x}}{\partial \textbf{q}}\right){\tau }_{xy}\right]dA.</math>  
 
|}
 
|}
|  style="vertical-align: top;width: 5px;text-align: right;white-space: nowrap;"|(5)
+
|  style="vertical-align: center;width: 5px;text-align: right;white-space: nowrap;"|(5)
 
|}
 
|}
  
 
+
sendo <math display="inline">\textbf{f}_{ext}</math> o vetor de forças externas nas direções dos graus de liberdade do elemento triangular, o trabalho virtual externo é dado por  <math display="inline">\delta W_{ext}=\delta \textbf{q}^T\textbf{f}_{ext}</math>,  onde  <math display="inline">\textbf{f}_{ext}</math> é o vetor de forças externas aplicadas diretamente na direção dos graus de liberdade do elemento, e as forças nodais equivalentes obtidas a partir do carregamento externo atuando no contorno do elemento. Da condição <math display="inline">\delta W_{ext}=\delta W_{int}</math> tem-se:
Sendo <math display="inline">\boldsymbol{f}_{ext}</math> o vetor de forças externas nas direções dos graus de liberdade do elemento triangular, o trabalho virtual externo é dado por  <math display="inline">\delta W_{ext}=\delta \boldsymbol{q}^T\boldsymbol{f}_{ext}</math> ,  onde  <math display="inline">\boldsymbol{f}_{ext}</math> é o vetor de forças externas aplicadas diretamente na direção dos graus de liberdade do elemento, e as forças nodais equivalentes obtidas a partir do carregamento externo atuando no contorno do elemento. Da condição <math display="inline">\delta W_{ext}=\delta W_{int}</math> tem-se:
+
  
 
{| class="formulaSCP" style="width: 100%;width: 100%;text-align: center;"  
 
{| class="formulaSCP" style="width: 100%;width: 100%;text-align: center;"  
Line 141: Line 145:
 
{| style="text-align: center; margin:auto;width: 100%;"
 
{| style="text-align: center; margin:auto;width: 100%;"
 
|-
 
|-
| <math>\delta \boldsymbol{q}^Tt\underset{A}{\int\!\int }[\frac{\partial u_{,x}}{\partial \boldsymbol{q}}{\sigma }_x+</math><math>\frac{\partial v_{,y}}{\partial \boldsymbol{q}}{\sigma }_y+</math><math>(\frac{\partial u_{,y}}{\partial \boldsymbol{q}}+\frac{\partial v_{,x}}{\partial \boldsymbol{q}}){\tau }_{xy}]dA=</math><math>\delta q^Tf_{ext}</math>  
+
| <math>\delta \textbf{q}^Tt\underset{A}{\int\!\int }\left[\frac{\partial u_{,x}}{\partial \textbf{q}}{\sigma }_x+\frac{\partial v_{,y}}{\partial \textbf{q}}{\sigma }_y+\left(\frac{\partial u_{,y}}{\partial \textbf{q}}+\frac{\partial v_{,x}}{\partial \textbf{q}}\right){\tau }_{xy}\right]dA=</math><math>\delta q^Tf_{ext}</math>  
 
|}
 
|}
|  style="text-align: right;vertical-align: top;width: 5px;text-align: right;white-space: nowrap;"|(6)
+
|  style="text-align: right;vertical-align: center;width: 5px;text-align: right;white-space: nowrap;"|(6)
 
|}
 
|}
  
 
+
Como a expressão acima deve ser válida para qualquer campo de deslocamento virtual  <math display="inline">\delta {\bf q}</math>, segue que <math display="inline">\textbf{f}_{int}-\textbf{f}_{ext}=\mbox{0}</math>, onde  <math display="inline">\textbf{f}_{int}</math> é o vetor de forças internas definido por:
Como a expressão acima deve ser válida para qualquer campo de deslocamento virtual  <math display="inline">\delta \boldsymbol{q}</math> , segue que <math display="inline">\boldsymbol{f}_{int}-\boldsymbol{f}_{ext}=\mbox{0}</math> , onde  <math display="inline">\boldsymbol{f}_{int}</math> é o vetor de forças internas definido por:
+
  
 
{| class="formulaSCP" style="width: 100%;width: 100%;text-align: center;"  
 
{| class="formulaSCP" style="width: 100%;width: 100%;text-align: center;"  
Line 154: Line 157:
 
{| style="text-align: center; margin:auto;width: 100%;"
 
{| style="text-align: center; margin:auto;width: 100%;"
 
|-
 
|-
| <math>\mbox{f}_{int}=t\underset{A}{\int\!\int }[\frac{\partial u_{,x}}{\partial \boldsymbol{q}}{\sigma }_x+</math><math>\frac{\partial v_{,y}}{\partial \boldsymbol{q}}{\sigma }_y+</math><math>(\frac{\partial u_{,y}}{\partial \boldsymbol{q}}+\frac{\partial v_{,x}}{\partial \boldsymbol{q}}){\tau }_{xy}]dA</math>  
+
| <math>f_{int}=t\underset{A}{\int\!\int }\left[\frac{\partial u_{,x}}{\partial \textbf{q}}{\sigma }_x+\frac{\partial v_{,y}}{\partial \textbf{q}}{\sigma }_y+\left(\frac{\partial u_{,y}}{\partial \textbf{q}}+\frac{\partial v_{,x}}{\partial \textbf{q}}\right){\tau }_{xy}\right]dA</math>  
 
|}
 
|}
|  style="text-align: right;vertical-align: top;width: 5px;text-align: right;white-space: nowrap;"|(7)
+
|  style="text-align: right;vertical-align: center;width: 5px;text-align: right;white-space: nowrap;"|(7)
 
|}
 
|}
 
  
 
A Eq. (7) pode ser escrita na forma matricial como:
 
A Eq. (7) pode ser escrita na forma matricial como:
Line 167: Line 169:
 
{| style="text-align: center; margin:auto;width: 100%;"
 
{| style="text-align: center; margin:auto;width: 100%;"
 
|-
 
|-
| <math>\boldsymbol{f}_{int}=t\underset{A}{\int\!\int }\left[\begin{array}{c}
+
| <math>\textbf{f}_{int}=t\underset{A}{\int\!\int }\left[\begin{array}{c}
 
{\sigma }_x{\Phi }_{,x}+{\tau }_{xy}{\Phi }_{,y}\\
 
{\sigma }_x{\Phi }_{,x}+{\tau }_{xy}{\Phi }_{,y}\\
 
{\sigma }_y{\Phi }_{,y}+{\tau }_{xy}{\Phi }_{,x}
 
{\sigma }_y{\Phi }_{,y}+{\tau }_{xy}{\Phi }_{,x}
 
\end{array}\right]\mbox{ }dA</math>  
 
\end{array}\right]\mbox{ }dA</math>  
 
|}
 
|}
|  style="text-align: right;vertical-align: top;width: 5px;text-align: right;white-space: nowrap;"|(8)
+
|  style="text-align: right;vertical-align: center;width: 5px;text-align: right;white-space: nowrap;"|(8)
 
|}
 
|}
  
 
+
Utilizando o método de Newton-Rapshon na solução do problema  <math display="inline">\textbf{f}_{int}-\textbf{f}_{ext}=\mbox{0}</math> é necessária a determinação da derivada dessa expressão em relação aos deslocamentos nodais, obtendo assim a matriz de rigidez tangente. Sendo  <math display="inline">\textbf{f}_{ext}</math> constante em relação aos deslocamentos nodais, a matriz de rigidez tangente é dada por  <math display="inline">\textbf{K}=\partial \textbf{f}_{int}/\partial \textbf{q}</math>, que, após manipulações algébricas, pode ser escrita de acordo com a Eq. (9)
Utilizando o método de Newton-Rapshon na solução do problema  <math display="inline">\boldsymbol{f}_{int}-\boldsymbol{f}_{ext}=\mbox{0}</math> é necessária a determinação da derivada dessa expressão em relação aos deslocamentos nodais, obtendo assim a matriz de rigidez tangente. Sendo  <math display="inline">\boldsymbol{f}_{ext}</math> constante em relação aos deslocamentos nodais, a matriz de rigidez tangente é dada por  <math display="inline">\boldsymbol{K}=\partial \boldsymbol{f}_{int}/\partial \boldsymbol{q}</math> , que, após manipulações algébricas, pode ser escrita de acordo com a Eq. (9).
+
  
 
{| class="formulaSCP" style="width: 100%;width: 100%;text-align: center;"  
 
{| class="formulaSCP" style="width: 100%;width: 100%;text-align: center;"  
Line 183: Line 184:
 
{| style="text-align: center; margin:auto;width: 100%;"
 
{| style="text-align: center; margin:auto;width: 100%;"
 
|-
 
|-
| <math>\boldsymbol{K}=t\underset{A}{\int\!\int }\left[\begin{array}{c}
+
| <math>\textbf{K}=t\underset{A}{\int\!\int }\left[\begin{array}{c}
{\Phi }_{,x}{\left(\frac{\partial {\sigma }_x}{\partial \boldsymbol{q}}\right)}^T+{\Phi }_{,y}{\left(\frac{\partial {\tau }_{xy}}{\partial \boldsymbol{q}}\right)}^T\\
+
{\Phi }_{,x}{\left(\frac{\partial {\sigma }_x}{\partial \textbf{q}}\right)}^T+{\Phi }_{,y}{\left(\frac{\partial {\tau }_{xy}}{\partial \textbf{q}}\right)}^T\\
{\Phi }_{,y}{\left(\frac{\partial {\sigma }_y}{\partial \boldsymbol{q}}\right)}^T+{\Phi }_{,x}{\left(\frac{\partial {\tau }_{xy}}{\partial \boldsymbol{q}}\right)}^T
+
{\Phi }_{,y}{\left(\frac{\partial {\sigma }_y}{\partial \textbf{q}}\right)}^T+{\Phi }_{,x}{\left(\frac{\partial {\tau }_{xy}}{\partial \textbf{q}}\right)}^T
\end{array}\right]\mbox{ }dA</math>  
+
\end{array}\right]dA</math>  
 
|}
 
|}
|  style="text-align: right;vertical-align: top;width: 5px;text-align: right;white-space: nowrap;"|(9)
+
|  style="text-align: right;vertical-align: center;width: 5px;text-align: right;white-space: nowrap;"|(9)
 
|}
 
|}
  
 +
onde, <math>\Phi_{,x}</math> e <math>\Phi_{,y}</math> são as derivadas das funções de forma com relação a <math>x</math> e <math>y</math>, respectivamente,  <math>{\sigma }_x</math> e  <math>{\sigma }_y</math> são as tensões normais nas direções x e y, respectivamente,  <math>{\tau }_{xy}</math> é a tensão cisalhante, <math>\bf q</math> é o vetor de deslocamento nodal, <math>t</math> é a espessura e <math>A</math> a área do elemento.
  
onde, Φ,x e Φ,y são as derivadas das funções de forma com relação a ''x ''e ''y'', respectivamente,  <math>{\sigma }_x</math> e  <math>{\sigma }_y</math> são as tensões normais nas direções x e y, respectivamente,  <math>{\tau }_{xy}</math> é a tensão cisalhante, '''q''' é o vetor de deslocamento nodal, ''t'' é a espessura e ''A'' a área do elemento.
+
Na Eq. (9), a derivada da tensão normal na direção <math>x</math> em relação aos deslocamentos nodais é dada por:
 
+
Na Eq. (9), a derivada da tensão normal na direção ''x ''em relação aos deslocamentos nodais é dada por:
+
  
 
{| class="formulaSCP" style="width: 100%;width: 100%;text-align: center;"  
 
{| class="formulaSCP" style="width: 100%;width: 100%;text-align: center;"  
Line 203: Line 203:
 
| <math>\frac{\partial {\sigma }_x}{\partial q}=\frac{\partial {\sigma }_x}{\partial {\epsilon }_x}\frac{\partial {\epsilon }_x}{\partial q}+</math><math>\frac{\partial {\sigma }_x}{\partial {\epsilon }_y}\frac{\partial {\epsilon }_y}{\partial q}</math>  
 
| <math>\frac{\partial {\sigma }_x}{\partial q}=\frac{\partial {\sigma }_x}{\partial {\epsilon }_x}\frac{\partial {\epsilon }_x}{\partial q}+</math><math>\frac{\partial {\sigma }_x}{\partial {\epsilon }_y}\frac{\partial {\epsilon }_y}{\partial q}</math>  
 
|}
 
|}
|  style="width: 5px;text-align: right;white-space: nowrap;"|(10)
+
|  style="width: 5px;text-align: right;vertical-align: center;white-space: nowrap;"|(10)
 
|}
 
|}
  
 
+
As derivadas das deformações lineares em relação aos deslocamentos nodais são dadas pelas Eqs. (11) e (12)
As derivadas das deformações lineares em relação aos deslocamentos nodais são dadas pelas Eq. (11) e (12).
+
  
 
{| class="formulaSCP" style="width: 100%;width: 100%;text-align: center;"  
 
{| class="formulaSCP" style="width: 100%;width: 100%;text-align: center;"  
Line 219: Line 218:
 
\end{array}\right]</math>  
 
\end{array}\right]</math>  
 
|}
 
|}
|  style="text-align: right;vertical-align: top;width: 5px;text-align: right;white-space: nowrap;"|(11)
+
|  style="text-align: right;vertical-align: center;width: 5px;text-align: right;white-space: nowrap;"|(11)
 
|-
 
|-
 
|  
 
|  
Line 229: Line 228:
 
\end{array}\right]</math>  
 
\end{array}\right]</math>  
 
|}
 
|}
|  style="text-align: right;vertical-align: top;width: 5px;text-align: right;white-space: nowrap;"|(12)
+
|  style="text-align: right;vertical-align: center;width: 5px;text-align: right;white-space: nowrap;"|(12)
 
|}
 
|}
  
  
Substituindo as Eq. (11) e (12) na derivada da tensão normal na direção ''x'' em relação aos deslocamentos nodais (Eq. (10)), tem-se:
+
Substituindo as Eqs. (11) e (12) na derivada da tensão normal na direção ''x'' em relação aos deslocamentos nodais (Eq. (10)), tem-se:
  
 
{| class="formulaSCP" style="width: 100%;width: 100%;text-align: center;"  
 
{| class="formulaSCP" style="width: 100%;width: 100%;text-align: center;"  
Line 245: Line 244:
 
\end{array}\right]</math>  
 
\end{array}\right]</math>  
 
|}
 
|}
|  style="text-align: right;vertical-align: top;width: 5px;text-align: right;white-space: nowrap;"|(13)
+
|  style="text-align: right;vertical-align: center;width: 5px;text-align: right;white-space: nowrap;"|(13)
 
|}
 
|}
  
 
+
De forma análoga à descrita para tensão normal na direção <math>x</math>, pode-se chegar às Eqs. (14) e (15) para as derivadas em relação aos deslocamentos nodais das outras tensões atuantes no elemento
De forma análoga à descrita para tensão normal na direção ''x'', pode-se chegar às Eq. (14) e (15) para as derivadas em relação aos deslocamentos nodais das outras tensões atuantes no elemento.
+
  
 
{| class="formulaSCP" style="width: 100%;width: 100%;text-align: center;"  
 
{| class="formulaSCP" style="width: 100%;width: 100%;text-align: center;"  
Line 260: Line 258:
 
{\Phi }_{,y}
 
{\Phi }_{,y}
 
\end{array}\right]</math>
 
\end{array}\right]</math>
 
 
 
|}
 
|}
|  style="text-align: right;vertical-align: top;width: 5px;text-align: right;white-space: nowrap;"|(14)
+
|  style="text-align: right;vertical-align: center;width: 5px;text-align: right;white-space: nowrap;"|(14)
 
|-
 
|-
 
|  
 
|  
Line 273: Line 269:
 
\end{array}\right]</math>  
 
\end{array}\right]</math>  
 
|}
 
|}
|  style="text-align: right;vertical-align: top;width: 5px;text-align: right;white-space: nowrap;"|(15)
+
|  style="text-align: right;vertical-align: center;width: 5px;text-align: right;white-space: nowrap;"|(15)
 
|}
 
|}
  
 
+
A seguir são determinadas as derivadas das funções de forma em relação aos eixos <math>x</math> e <math>y</math>.
A seguir são determinadas as derivadas das funções de forma em relação aos eixos ''x ''e'' y''.
+
  
 
{| class="formulaSCP" style="width: 100%;width: 100%;text-align: center;"  
 
{| class="formulaSCP" style="width: 100%;width: 100%;text-align: center;"  
Line 298: Line 293:
 
\end{array}\right\}</math>  
 
\end{array}\right\}</math>  
 
|}
 
|}
|  style="text-align: right;vertical-align: top;width: 5px;text-align: right;white-space: nowrap;"|(16)
+
|  style="text-align: right;vertical-align: center;width: 5px;text-align: right;white-space: nowrap;"|(16)
 
|-
 
|-
 
|  
 
|  
Line 317: Line 312:
 
\end{array}\right\}</math>  
 
\end{array}\right\}</math>  
 
|}
 
|}
|  style="text-align: right;vertical-align: top;width: 5px;text-align: right;white-space: nowrap;"|(17)
+
|  style="text-align: right;vertical-align: center;width: 5px;text-align: right;white-space: nowrap;"|(17)
 
|}
 
|}
  
 
+
Na representação isoparamétrica, as coordenadas cartesianas <math>x</math> e <math>y</math> são relacionadas com as coordenadas paramétricas <math>\xi</math> e <math>\eta</math>. Então, para mudança do domínio de integração <math>dA= dxdy</math> para  <math>d\xi d\eta </math> utiliza-se a relação  <math>dA=det\mbox{J}d\xi d\eta </math> onde <math>\bf J</math> é a matriz jacobiano da transformação das coordenadas <math>x</math> e <math>y</math> para as coordenadas paramétricas  <math>\xi </math> e  <math>\eta </math> dada pela equação seguir. Dessa forma, <math>det{\bf J} = 2A</math>
Na representação isoparamétrica, as coordenadas cartesianas ''x ''e'' y'' são relacionadas com as coordenadas paramétricas ξ e η. Então, para mudança do domínio de integração d''A''= d''x''d''y ''para  <math>d\xi d\eta </math> utiliza-se a relação  <math>dA=det\mbox{J}d\xi d\eta </math> onde '''J''' é a matriz jacobiano da transformação das coordenadas ''x ''e'' y'' para as coordenadas paramétricas  <math>\xi </math> e  <math>\eta </math> dada pela equação seguir. Dessa forma, det'''J '''= 2''A''.
+
  
 
{| class="formulaSCP" style="width: 100%;width: 100%;text-align: center;"  
 
{| class="formulaSCP" style="width: 100%;width: 100%;text-align: center;"  
Line 333: Line 327:
 
\end{array}\right]</math>  
 
\end{array}\right]</math>  
 
|}
 
|}
|  style="text-align: right;vertical-align: top;width: 5px;text-align: right;white-space: nowrap;"|(18)
+
|  style="text-align: right;vertical-align: center;width: 5px;text-align: right;white-space: nowrap;"|(18)
 
|}
 
|}
  
 
+
==3. Otimização estrutural evolucionária==
==3 Otimização Estrutural Evolucionária==
+
  
 
A técnica ESO surge como uma alternativa ao rigor matemático de métodos clássicos de otimização. Apresenta uma base teórica simples cujo fundamento consiste na inserção de vazios na estrutura através da eliminação gradual dos elementos menos solicitados do domínio durante o processo de evolução. Portanto, para se obter a configuração ótima deve-se agregar ao estudo um nível de análise estrutural dependente de um domínio discreto, o que torna o emprego do MEF uma etapa do algoritmo de otimização.
 
A técnica ESO surge como uma alternativa ao rigor matemático de métodos clássicos de otimização. Apresenta uma base teórica simples cujo fundamento consiste na inserção de vazios na estrutura através da eliminação gradual dos elementos menos solicitados do domínio durante o processo de evolução. Portanto, para se obter a configuração ótima deve-se agregar ao estudo um nível de análise estrutural dependente de um domínio discreto, o que torna o emprego do MEF uma etapa do algoritmo de otimização.
Line 343: Line 336:
 
Neste trabalho a representação matemática técnica ESO baseia-se no conceito de tensão, isto é, o nível máximo de tensão na estrutura, obtido por análises via MEF, é tomado como um indicador do nível de eficiência de cada elemento. Elementos com baixo nível de tensão são, portanto, sistematicamente removidos da estrutura. A cada iteração novos elementos ineficientes são eliminados da malha e o procedimento se repete até que o campo de tensão atuante em todo o domínio seja praticamente constante e muito próximo da tensão admissível do material ou que seja atingida a restrição de volume mínimo.
 
Neste trabalho a representação matemática técnica ESO baseia-se no conceito de tensão, isto é, o nível máximo de tensão na estrutura, obtido por análises via MEF, é tomado como um indicador do nível de eficiência de cada elemento. Elementos com baixo nível de tensão são, portanto, sistematicamente removidos da estrutura. A cada iteração novos elementos ineficientes são eliminados da malha e o procedimento se repete até que o campo de tensão atuante em todo o domínio seja praticamente constante e muito próximo da tensão admissível do material ou que seja atingida a restrição de volume mínimo.
  
O critério de remoção é feito comparando-se a tensão de von Mises de cada elemento com a tensão de von Mises máxima existente em toda a estrutura. Portanto, no fim de cada iteração todos os elementos que atendem à Eq. (19) serão eliminados. A forma de retirada do elemento ocorre atribuindo-se baixos valores para seu módulo de elasticidade longitudinal (''E''=10<sup>-12</sup>). Desse modo, evita-se o remalhamento da estrutura.
+
O critério de remoção é feito comparando-se a tensão de von Mises de cada elemento com a tensão de von Mises máxima existente em toda a estrutura. Portanto, no fim de cada iteração todos os elementos que atendem à Eq. (19) serão eliminados. A forma de retirada do elemento ocorre atribuindo-se baixos valores para seu módulo de elasticidade longitudinal (<math>E=10^{-12}</math>). Desse modo, evita-se o remalhamento da estrutura
  
 
{| class="formulaSCP" style="width: 100%;width: 100%;text-align: center;"  
 
{| class="formulaSCP" style="width: 100%;width: 100%;text-align: center;"  
Line 352: Line 345:
 
| <math>{\sigma }_e^{\upsilon M}<RR_i.\mbox{ }{\sigma }_{Max}^{\upsilon M}</math>  
 
| <math>{\sigma }_e^{\upsilon M}<RR_i.\mbox{ }{\sigma }_{Max}^{\upsilon M}</math>  
 
|}
 
|}
|  style="text-align: right;vertical-align: top;width: 5px;text-align: right;white-space: nowrap;"|(19)
+
|  style="text-align: right;vertical-align: center;width: 5px;text-align: right;white-space: nowrap;"|(19)
 
|}
 
|}
 
  
 
onde:
 
onde:
Line 360: Line 352:
 
<math display="inline">{\sigma }_{e}^{vM}</math>= tensão de von Mises no elemento analisado;
 
<math display="inline">{\sigma }_{e}^{vM}</math>= tensão de von Mises no elemento analisado;
  
''RRi''= razão de rejeição na i-ésima iteração (0< ''RRi ''<1,0);
+
<math>RRi</math>= razão de rejeição na i-ésima iteração (<math>0< RRi <1,0</math>);
  
 
<math display="inline">{\sigma }_{Max}^{vM}</math>= máxima tensão de von Mises da iteração.
 
<math display="inline">{\sigma }_{Max}^{vM}</math>= máxima tensão de von Mises da iteração.
  
A razão de rejeição é usada para retardar o processo de remoção do elemento. O ciclo de remoção ocorre até que não possam mais ser removidos elementos para um dado valor de ''RRi.'' Quando isto ocorre, um estado de equilíbrio é alcançado. O processo evolucionário é redefinido adicionando-se à ''RRi'' uma razão de evolução, ''ER''. A razão de rejeição é atualizada de acordo com a Eq. (20).
+
A razão de rejeição é usada para retardar o processo de remoção do elemento. O ciclo de remoção ocorre até que não possam mais ser removidos elementos para um dado valor de <math>RRi</math>. Quando isto ocorre, um estado de equilíbrio é alcançado. O processo evolucionário é redefinido adicionando-se à <math>RRi</math> uma razão de evolução, <math>ER</math>. A razão de rejeição é atualizada de acordo com a Eq. (20)
  
 
{| class="formulaSCP" style="width: 100%;width: 100%;text-align: center;"  
 
{| class="formulaSCP" style="width: 100%;width: 100%;text-align: center;"  
Line 373: Line 365:
 
| <math>RR_{i+1}=RR_i+ER</math>  
 
| <math>RR_{i+1}=RR_i+ER</math>  
 
|}
 
|}
|  style="text-align: right;vertical-align: top;width: 5px;text-align: right;white-space: nowrap;"|(20)
+
|  style="text-align: right;vertical-align: center;width: 5px;text-align: right;white-space: nowrap;"|(20)
 
|}
 
|}
  
 
+
O valor inicial da razão de rejeição (<math>RR_0</math>) é definido de forma empírica pelo usuário. Entretanto, segundo Querin [3], para garantir melhor convergência, os valores de <math>RR_0</math> e <math>ER</math> devem ser de aproximadamente 1%. O processo se repete enquanto a estrutura não atingir o volume final, <math>VF</math>, definido pelo usuário, ou seja:
O valor inicial da razão de rejeição (''RR<sub>0</sub>'') é definido de forma empírica pelo usuário. Entretanto, segundo Querin [3], para garantir melhor convergência, os valores de ''RR<sub>0</sub> ''e ''ER'' devem ser de aproximadamente 1%. O processo se repete enquanto a estrutura não atingir o volume final, ''VF'', definido pelo usuário, ou seja:
+
  
 
{| class="formulaSCP" style="width: 100%;width: 100%;text-align: center;"  
 
{| class="formulaSCP" style="width: 100%;width: 100%;text-align: center;"  
Line 386: Line 377:
 
| <math>VR<(1-VF)\mbox{ }.\mbox{ }VT</math>  
 
| <math>VR<(1-VF)\mbox{ }.\mbox{ }VT</math>  
 
|}
 
|}
|  style="text-align: right;vertical-align: top;width: 5px;text-align: right;white-space: nowrap;"|(21)
+
|  style="text-align: right;vertical-align: center;width: 5px;text-align: right;white-space: nowrap;"|(21)
 
|}
 
|}
  
 
+
onde <math>VR</math> é o volume  retirado acumulado até aquela iteração, <math>VF</math> é o volume final expresso em percentual (por exemplo, <math>VF= 0,4</math>  implica que a retirada de elementos cessará quando o volume da estrutura atingir 40% do seu volume total inicial) e <math>VT</math> o volume inicial total da estrutura.
onde ''VR'' é o volume  retirado acumulado até aquela iteração, ''VF'' é o volume final expresso em percentual (por exemplo, ''VF'' = 0,4 implica que a retirada de elementos cessará quando o volume da estrutura atingir 40% do seu volume total inicial) e ''VT'' o volume inicial total da estrutura.
+
  
 
Matematicamente, o ESO pode ser escrito como:
 
Matematicamente, o ESO pode ser escrito como:
Line 399: Line 389:
 
{| style="text-align: center; margin:auto;width: 100%;"
 
{| style="text-align: center; margin:auto;width: 100%;"
 
|-
 
|-
| <math>\mbox{D }(j)=\lbrace \begin{array}{c}
+
| <math>\mbox{D }(j)=\left\{ \begin{array}{c}
 
\mbox{D}_0,\\
 
\mbox{D}_0,\\
 
\mbox{0,}
 
\mbox{0,}
Line 414: Line 404:
 
\Gamma \\
 
\Gamma \\
 
\Gamma \, {}'
 
\Gamma \, {}'
\end{array}</math>  
+
\end{array}\right.</math>  
 
|}
 
|}
|  style="text-align: right;vertical-align: top;width: 5px;text-align: right;white-space: nowrap;"|(22)
+
|  style="text-align: right;vertical-align: center;width: 5px;text-align: right;white-space: nowrap;"|(22)
 
|}
 
|}
 
  
 
onde:
 
onde:
  
'''D'''(j) = matriz constitutiva do ponto j <math display="inline">\in \,</math> Ω;
+
<math>{\bf D}(j)</math> = matriz constitutiva do ponto <math display="inline"> j\in \, \Omega</math>;
  
'''D'''<sub>0</sub> = matriz constitutiva inicial;
+
<math>{\bf D}_0</math> = matriz constitutiva inicial;
  
Ω = <math display="inline">\Gamma +\Gamma '</math> = domínio da estrutura;
+
<math>\Omega = \Gamma +\Gamma '</math> = domínio da estrutura;
  
<math>\Gamma =\left\{\Omega \mbox{/}\left(\frac{{\sigma }_e^{vM}}{{\sigma }_{m\acute{a}x}^{vM}}\right)\geq {\mbox{RR}}_i\right\}</math> conjunto dos elementos que não serão removidos.
+
<math>\Gamma =\left\{\Omega /\left(\frac{{\sigma }_e^{vM}}{{\sigma }_{m\acute{a}x}^{vM}}\right)\geq {\mbox{RR}}_i\right\}</math> conjunto dos elementos que não serão removidos.
  
<math>\Gamma \, {}'=\left\{\Omega \mbox{/}\left(\frac{{\sigma }_e^{vM}}{{\sigma }_{m\acute{a}x}^{vM}}\right)<{\mbox{RR}}_i\right\}</math> , conjunto dos elementos que serão removidos.
+
<math>\Gamma \, {}'=\left\{\Omega /\left(\frac{{\sigma }_e^{vM}}{{\sigma }_{m\acute{a}x}^{vM}}\right)<{\mbox{RR}}_i\right\}</math>, conjunto dos elementos que serão removidos.
 +
 
 +
Portanto, o algoritmo ESO apresenta a seguinte marcha, representada no fluxograma da [[#img-1|Figura 1]]:
  
Portanto, o algoritmo ESO apresenta a seguinte marcha, representada no fluxograma da Figura 1:
 
  
 
1º Passo: discretização do domínio e aplicação das condições de contorno e ações prescritas;
 
1º Passo: discretização do domínio e aplicação das condições de contorno e ações prescritas;
Line 438: Line 428:
 
2º Passo: análise da estrutura via MEF e cálculo das tensões principais e tensões de von Mises em cada elemento;
 
2º Passo: análise da estrutura via MEF e cálculo das tensões principais e tensões de von Mises em cada elemento;
  
3º Passo: retirar os elementos que satisfaçam a Eq. 19, dentro de um limite pré-definido de volume (p%);
+
3º Passo: retirar os elementos que satisfaçam a Eq. (19), dentro de um limite pré-definido de volume (p%);
  
 
4º Passo: repetir os passos 2 e 3 até que seja atingido o equilíbrio;
 
4º Passo: repetir os passos 2 e 3 até que seja atingido o equilíbrio;
  
5º Passo: acréscimo da razão de rejeição conforme a Eq. 20 e iniciar nova retirada de elementos repetindo os passos 2, 3 e 4.
+
5º Passo: acréscimo da razão de rejeição conforme a Eq. (20) e iniciar nova retirada de elementos repetindo os passos 2, 3 e 4.
 
+
<div id='img-1'></div>
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
+
{| style="text-align: center; border: 1px solid #BBB; margin: 1em auto; width: auto%;max-width: auto;"
[[Image:Draft_Ladeira_717005243-image51.png|600px]] </div>
+
|-
 
+
| style="padding:10px;" | [[Image:Draft_Ladeira_717005243-image51.png|600px]]
<div id="_Toc21935919" class="center" style="width: auto; margin-left: auto; margin-right: auto;">
+
|- style="text-align: center; font-size: 75%;"
Figura 1 Fluxograma do algoritmo ESO em nível de tensão.</div>
+
| colspan="1" style="padding-bottom:10px;" | '''Figura 1'''. Fluxograma do algoritmo ESO em nível de tensão
 +
|}
  
==4 Critério de Escoamento do von Mises==
+
==4. Critério de escoamento do von Mises==
  
 
O critério de remoção é baseado na tensão de von Mises de cada elemento na i-ésima iteração, que em termos de tensões principais pode ser calculada conforme a Eq. (23):
 
O critério de remoção é baseado na tensão de von Mises de cada elemento na i-ésima iteração, que em termos de tensões principais pode ser calculada conforme a Eq. (23):
Line 478: Line 469:
  
  
==5 Exemplos Numéricos==
+
==5. Exemplos numéricos==
  
 
Neste item são apresentadas algumas aplicações considerando o método ESO implementado. São feitas análises elástico-lineares em estruturas submetidas a estado plano de tensão. O material é considerado homogêneo e isotrópico.
 
Neste item são apresentadas algumas aplicações considerando o método ESO implementado. São feitas análises elástico-lineares em estruturas submetidas a estado plano de tensão. O material é considerado homogêneo e isotrópico.
  
Para cada exemplo são definidas as propriedades mecânicas do material e o domínio inicial de projeto. Além disso, são especificados os seguintes parâmetros: ''RR'' (razão de rejeição), ''ER'' (razão de evolução), ''VF'' (volume final desejado), ''VI ''(volume máximo retirado por iteração), ''VR'' (volume total retirado até uma dada iteração), número da iteração e malha de elementos finitos adotada. Nesse item podem ser feitos os agradecimentos às instituições de fomento à pesquisa, às empresas privadas patrocinadoras do trabalho, aos pesquisadores e profissionais que auxiliaram na execução dos ensaios e fornecimento de dados e materiais.
+
Para cada exemplo são definidas as propriedades mecânicas do material e o domínio inicial de projeto. Além disso, são especificados os seguintes parâmetros: <math>RR</math> (razão de rejeição), <math>ER</math> (razão de evolução), <math>VF</math> (volume final desejado), <math>VI</math> (volume máximo retirado por iteração), <math>VR</math> (volume total retirado até uma dada iteração), número da iteração e malha de elementos finitos adotada. Nesse item podem ser feitos os agradecimentos às instituições de fomento à pesquisa, às empresas privadas patrocinadoras do trabalho, aos pesquisadores e profissionais que auxiliaram na execução dos ensaios e fornecimento de dados e materiais.
  
 
===5.1 Estrutura de Michell===
 
===5.1 Estrutura de Michell===
  
O primeiro exemplo a ser apresentado trata-se de uma viga simplesmente apoiada submetida a uma carga concentrada conforme indicado na Figura 2a, comumente chamada na literatura de estrutura de Michell. A solução analítica é mostrada na Figura 2b. O material adotado foi o aço, cujo módulo de elasticidade ''E=''200GPa, coeficiente de Poisson ''ʋ=''0.3 e espessura igual a 1mm.
+
O primeiro exemplo a ser apresentado trata-se de uma viga simplesmente apoiada submetida a uma carga concentrada conforme indicado na [[#img-2|Figura 2]]a, comumente chamada na literatura de estrutura de Michell. A solução analítica é mostrada na [[#img-2|Figura 2]]b. O material adotado foi o aço, cujo módulo de elasticidade <math>E=200</math>GPa, coeficiente de Poisson <math> \nu =0.3</math> e espessura igual a 1mm.
 +
<div id='img-2'></div>
 +
{| style="text-align: center; border: 1px solid #BBB; margin: 1em auto; width: auto%;max-width: auto;"
 +
|-
 +
| style="padding:10px;" | [[Image:Draft_Ladeira_717005243-image55.png|600px]]
 +
|- style="text-align: center; font-size: 75%;"
 +
| colspan="1" style="padding-bottom:10px;" | '''Figura 2'''. (a) Domínio inicial para estrutura de Michell. (b) Solução analítica
 +
|}
  
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
 
[[Image:Draft_Ladeira_717005243-image55.png|600px]] </div>
 
  
<div id="_Toc21935924" class="center" style="width: auto; margin-left: auto; margin-right: auto;">
+
O domínio foi discretizado numa malha de elementos triangulares de <math>96\times 40</math>. O processo evolucionário teve início com uma razão de rejeição (<math>RR</math>) de 1% e uma razão de evolução (<math>ER</math>) de 0,75%. O volume retirado (<math>VR</math>) de 60% do volume inicial e o volume máximo retirado por iteração (<math>VI</math>) de 1,75%. As [[#img-3|Figuras 3]]a, [[#img-3|3]]b e [[#img-3|3]]c apresentam a evolução da estrutura, enquanto que a [[#img-3|Figura 3]]d ilustra o modelo de bielas e tirantes obtido, onde as bielas estão representadas em azul e os tirantes em vermelho.
<span style="text-align: center; font-size: 75%;">Figura 2 (a) Domínio inicial para estrutura de Michell (b) Solução analítica</span></div>
+
<div id='img-3'></div>
 +
{| style="text-align: center; border: 1px solid #BBB; margin: 1em auto; width: 60%;"
 +
|-
 +
| style="padding:10px;" | [[Image:Draft_Ladeira_717005243-image56.png|468px]]
 +
|- style="text-align: center; font-size: 75%;"
 +
| colspan="1" style="padding:10px;" | '''Figura 3'''. (a) Iteração 22, <math>RR=2,5%</math>, <math>VR=5,6%</math>. (b) Iteração 69, <math>RR=4,75%</math>, <math>VR=18,6%</math>. (c) Iteração 175, <math>RR=10,75%</math>, <math>VR=60,0%</math>. (d) Modelo de bielas e tirantes
 +
|}
  
O domínio foi discretizado numa malha de elementos triangulares de 96x40. O processo evolucionário teve início com uma razão de rejeição (''RR'') de 1% e uma razão de evolução (''ER'') de 0,75%. O volume retirado (''VR'') de 60% do volume inicial e o volume máximo retirado por iteração (''VI'') de 1,75%. As Figura 3a, 3b e 3c apresentam a evolução da estrutura, enquanto que a Figura 3d ilustra o modelo de bielas e tirantes obtido, onde as bielas estão representadas em azul e os tirantes em vermelho.
+
===5.2 Consolo curto===
  
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
+
Este exemplo apresenta um consolo curto projetado para suportar uma carga pontual de 500kN. A [[#img-4|Figura 4]] traz as dimensões da estrutura em milímetros. O módulo de elasticidade do concreto foi tomado igual a <math> E=28567</math> MPa, coeficiente de Poisson <math> \nu =0,15</math>  e a espessura assumida como sendo igual a 300mm.
[[Image:Draft_Ladeira_717005243-image56.png|468px]] </div>
+
  
<div id="_Toc21935925" class="center" style="width: auto; margin-left: auto; margin-right: auto;">
+
<div id='img-4'></div>
<span style="text-align: center; font-size: 75%;">Figura 3 (a) Iteração 22, ''RR=''2,5%'', VR=''5,6% (b) Iteração 69, ''RR=''4,75%'', VR=''18,6% (c) Iteração 175, ''RR=''10,75%'', VR=''60,0% (d) Modelo de bielas e tirantes</span></div>
+
{| style="text-align: center; border: 1px solid #BBB; margin: 1em auto; width: 30%;;"
 +
|-
 +
| style="padding:10px;" | [[Image:Draft_Ladeira_717005243-image57.png|244px]]
 +
|- style="text-align: center; font-size: 75%;"
 +
| colspan="1" style="padding-bottom:10px;" | '''Figura 4'''. Consolo Curto (Liang et al. [4])
 +
|}
  
===5.2 Consolo Curto===
 
  
Este exemplo apresenta um consolo curto projetado para suportar uma carga pontual de 500kN. A Figura 4 traz as dimensões da estrutura em milímetros. O módulo de elasticidade do concreto foi tomado igual a ''E=''28567MPa, coeficiente de Poisson ''ʋ ''= 0,15 e a espessura assumida como sendo igual a 300mm.
+
Este exemplo foi estudado por Liang et al. [4] e Almeida et al. [5]. A [[#img-5|Figura 5]] apresenta a evolução da estrutura e o correspondente modelo de bielas e tirantes sugerido pelos autores, no qual as bielas estão representadas por linha pontilhada e os tirantes por linha cheia.
  
[[Image:Draft_Ladeira_717005243-image57.png|144px]]
+
<div id='img-5'></div>
 +
{| style="text-align: center; border: 1px solid #BBB; margin: 1em auto; width: auto%;max-width: auto;"
 +
|-
 +
| style="padding:10px;" | [[Image:Draft_Ladeira_717005243-image58.png|600px]]
 +
|- style="text-align: center; font-size: 75%;"
 +
| colspan="1" style="padding-bottom:10px;" | '''Figura 5'''. (a), (b), (c) Processo evolucionário.  (d) Modelo de bielas e tirantes (Liang et al. [4])
 +
|}
  
<div id="_Toc21935926" class="center" style="width: auto; margin-left: auto; margin-right: auto;">
 
<span style="text-align: center; font-size: 75%;">Figura 4 Consolo Curto (Liang ''et al.''[4])</span></div>
 
  
Este exemplo foi estudado por Liang ''et al''. [4] e Almeida ''et al''. [5]. A Figura 5 apresenta a evolução da estrutura e o correspondente modelo de bielas e tirantes sugerido pelos autores, no qual as bielas estão representadas por linha pontilhada e os tirantes por linha cheia.
+
No presente trabalho a estrutura foi modelada usando uma malha com 5664 elementos triangulares de três nós com 25 mm de lado. Para obter a topologia ótima apresentada na [[#img-6|Figura 6]], os parâmetros adotados foram: <math>RR = 0,01</math>, <math>ER=0,05</math>, <math>VF=45%</math>, <math>VI=0,05</math>. Os campos de compressão (bielas) são representados em azul e os campos de tração (tirantes) em vermelho.
  
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
+
<div id='img-6'></div>
[[Image:Draft_Ladeira_717005243-image58.png|600px]] </div>
+
{| style="text-align: center; border: 1px solid #BBB; margin: 1em auto; width: auto%;max-width: auto;"
 +
|-
 +
| style="padding:10px;" | [[Image:Draft_Ladeira_717005243-image59.png|500px]]
 +
|- style="text-align: center; font-size: 75%;"
 +
| colspan="1" style="padding:10px;" | '''Figura 6'''. (a) Iteração 25, <math>RR=11,0%</math>, <math>VR=15,0%</math>. (b) Iteração 31,  <math>RR=11,0%</math>, <math>VR=35,0%</math>.  (c) Iteração 47, <math>RR=16,0%</math>, <math>VR=55,0%</math>
 +
|}
  
<div id="_Toc21935927" class="center" style="width: auto; margin-left: auto; margin-right: auto;">
+
===5.3 Pilar de ponte===
<span style="text-align: center; font-size: 75%;">Figura 5 (a), (b), (c) Processo evolucionário  (d) Modelo de bielas e tirantes (Liang ''et al.'' [4])</span></div>
+
  
No presente trabalho a estrutura foi modelada usando uma malha com 5664 elementos triangulares de três nós com 25 mm de lado. Para obter a topologia ótima apresentada na Figura 6, os parâmetros adotados foram: ''RR = ''0,01,'' ER=''0,05, ''VF=''45%'', VI=''0,05. Os campos de compressão (bielas) são representados em azul e os campos de tração (tirantes) em vermelho.
+
Este exemplo foi proposto por Liang et al. [6] e estudado também por Almeida et al. [7]. Trata-se de um pilar de ponte projetado para suportar quatro cargas concentradas de 2750kN transferidas por quatro vigas de aço-concreto. O pilar tem espessura de 1,5m e é admitido como sendo engastado na fundação. A geometria, condições de contorno e ações do problema estão indicados na [[#img-7|Figura 7]], com dimensões em milímetros e kN.
  
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
+
<div id='img-7'></div>
[[Image:Draft_Ladeira_717005243-image59.png|600px]] </div>
+
{| style="text-align: center; border: 1px solid #BBB; margin: 1em auto; width: auto%;max-width: auto;"
 +
|-
 +
| style="padding:10px;" | [[Image:Draft_Ladeira_717005243-image60.png|354px]]
 +
|- style="text-align: center; font-size: 75%;"
 +
| colspan="1" style="padding-bottom:10px;" | '''Figura 7'''. Domínio de projeto do pilar de ponte (Liang et al. [6])
 +
|}
  
<div id="_Toc21935928" class="center" style="width: auto; margin-left: auto; margin-right: auto;">
+
O domínio foi discretizado numa malha refinada com 12260 elementos finitos triangulares. As propriedades do material isotrópico utilizado são módulo de elasticidade <math>E = 28,6</math>GPa e coeficiente de Poisson <math>\nu = 0,15</math>.
<span style="text-align: center; font-size: 75%;">Figura 6 (a) Iteração 25, ''RR = ''11,0%, ''VR = ''15,0%  (b) Iteração 31,  ''RR = ''11,0%, ''VR = ''35,0%  (c) Iteração 47, ''RR= ''16,0''%'','' VR=''55,0%</span></div>
+
  
===5.3 Pilar de Ponte===
+
O processo evolucionário teve início com uma razão de rejeição (<math>RR</math>) igual a 4% e uma razão de evolução (<math>ER</math>) de 2%. O volume final desejado igual a 48% do volume inicial e taxa de retirada de material por iteração (<math>VI</math>) igual a 1,75%. A [[#img-8|Figura 8]] apresenta o MBT alcançado. As regiões em vermelho e azul indicam, respectivamente, regiões de tração, tirantes, e regiões de compressão, bielas.
  
Este exemplo foi proposto por Liang ''et al''. [6] e estudado também por Almeida ''et al''. [7]. Trata-se de um pilar de ponte projetado para suportar quatro cargas concentradas de 2750kN transferidas por quatro vigas de aço-concreto. O pilar tem espessura de 1,5m e é admitido como sendo engastado na fundação. A geometria, condições de contorno e ações do problema estão indicados na Figura 7, com dimensões em milímetros e kN.
+
<div id='img-8'></div>
 
+
{| style="text-align: center; border: 1px solid #BBB; margin: 1em auto; width: 60%;"
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
+
|-
[[Image:Draft_Ladeira_717005243-image60.png|354px]] </div>
+
| style="padding:10px;" | [[Image:Draft_Ladeira_717005243-image61.png|528px]]  
 
+
|- style="text-align: center; font-size: 75%;"
<div id="_Toc21935959" class="center" style="width: auto; margin-left: auto; margin-right: auto;">
+
| colspan="1" style="padding:10px;" | '''Figura 8'''. (a) Modelo de bielas e tirantes obtido no presente trabalho - Iteração 164, <math>RR=18%</math>, <math>VR=52%</math>. (b) Topologia ótima proposta por Liang et al. [6]
<span style="text-align: center; font-size: 75%;">Figura 7 Domínio de projeto do pilar de ponte (Liang ''et al''., [6])</span></div>
+
|}
  
O domínio foi discretizado numa malha refinada com 12260 elementos finitos triangulares. As propriedades do material isotrópico utilizado são módulo de elasticidade ''E'' = 28,6GPa e coeficiente de Poisson ʋ = 0,15.
+
==6. Conclusões==
 
+
O processo evolucionário teve início com uma razão de rejeição (''RR'') igual a 4% e uma razão de evolução (''ER'') de 2%. O volume final desejado igual a 48% do volume inicial e taxa de retirada de material por iteração (''VI'') igual a 1,75%. A Figura 8 apresenta o MBT alcançado. As regiões em vermelho e azul indicam, respectivamente, regiões de tração, tirantes, e regiões de compressão, bielas.
+
 
+
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
+
[[Image:Draft_Ladeira_717005243-image61.png|528px]] </div>
+
 
+
<div id="_Toc21935960" class="center" style="width: auto; margin-left: auto; margin-right: auto;">
+
<span style="text-align: center; font-size: 75%;">Figura 8 (a) Modelo de bielas e tirantes obtido no presente trabalho - Iteração 164, ''RR''=18%, ''VR''=52% (b) Topologia ótima proposta por Liang ''et al.'' [6]</span></div>
+
 
+
==6 Conclusões==
+
  
 
O objetivo do artigo é a apresentação de uma formulação numérica para verificar o fluxo de tensões em estruturas bidimensionais para geração automática do modelo de bielas e tirantes que auxiliem o projetista na concepção e posterior dimensionamento de estruturas de concreto armado. O algoritmo de otimização topológica ESO foi, portanto, utilizado com esse propósito. O uso do Método dos Elementos Finitos, considerando análise linear elástica, se torna conveniente devido ao critério de remoção adotado para o ESO.
 
O objetivo do artigo é a apresentação de uma formulação numérica para verificar o fluxo de tensões em estruturas bidimensionais para geração automática do modelo de bielas e tirantes que auxiliem o projetista na concepção e posterior dimensionamento de estruturas de concreto armado. O algoritmo de otimização topológica ESO foi, portanto, utilizado com esse propósito. O uso do Método dos Elementos Finitos, considerando análise linear elástica, se torna conveniente devido ao critério de remoção adotado para o ESO.
Line 557: Line 564:
 
Para o dimensionamento, alerta-se para o fato de que o modelo de bielas e tirantes deve atender aos procedimentos impostos por códigos normativos, principalmente no que se refere à utilização de armadura mínima, que podem ser maiores que os tirantes concebidos.
 
Para o dimensionamento, alerta-se para o fato de que o modelo de bielas e tirantes deve atender aos procedimentos impostos por códigos normativos, principalmente no que se refere à utilização de armadura mínima, que podem ser maiores que os tirantes concebidos.
  
==7 Agradecimentos==
+
==Agradecimentos==
  
 
Os autores agradecem à CAPES, CNPq, FAPEMIG, Fundação Gorceix e PROPEC/UFOP pelo apoio financeiro e suporte necessários para realização desta pesquisa.
 
Os autores agradecem à CAPES, CNPq, FAPEMIG, Fundação Gorceix e PROPEC/UFOP pelo apoio financeiro e suporte necessários para realização desta pesquisa.
  
==8 Referências bibliográficas==
+
==Referências==
 +
 
 +
<div class="auto" style="text-align: left;width: auto; margin-left: auto; margin-right: auto;font-size: 85%;">
  
[1] Schlaich, J., Schäfer, K., Jennewein, M. Toward a Consistent Design of Structural Concrete. PCI Journal, vol. 32, n. 3, pp.75-150, 1987.
+
[1] Schlaich J., Schäfer K., Jennewein M. Toward a consistent design of structural concrete. PCI Journal, 32(3):75-150, 1987.
  
[2] Xie, Y. M.; Steven, G. P. A simple evolutionary procedure for structural optimization. Computers&Structures, v. 49, Issue5, p. 885-896, 1993.
+
[2] Xie Y.M., Steven G.P. A simple evolutionary procedure for structural optimization. Computers & Structures, 49(5):885-896, 1993.
  
[3] Querin, O.M., Evolutionary Structural Optimization stress based formulation and implementation. PhD dissertation, University of Sydney,1997.
+
[3] Querin O.M., Evolutionary structural optimization stress based formulation and implementation. PhD dissertation, University of Sydney, 1997.
  
[4] Michell, A. G. M''. ''The limits of economy of material in frame-structures. Philosophy Magazine, 8, p.589–597, 1904.
+
[4] Michell A.G.M. The limits of economy of material in frame-structures. Philosophy Magazine, 8:589–597, 1904.
  
[5] Liang, Q. Q.; Xie Y. M.; Steven, G.P. Topology optimization of strut-and-tie models in reinforced concrete structures using an evolutionary procedure. ACI Structural  Journal. 97(2):322–330, 2000.
+
[5] Liang Q.Q., Xie Y.M., Steven G.P. Topology optimization of strut-and-tie models in reinforced concrete structures using an evolutionary procedure. ACI Structural  Journal, 97(2):322–330, 2000.
  
[6] Almeida, V. S.; Simonetti, H. L. and Neto, L. O. Comparative analysis of strut-and-tie models using Smooth Evolutionary Structural Optimization. Engineering Structures, 56, 1665–1675, 2013a.
+
[6] Almeida V.S., Simonetti H.L., Neto L.O. Comparative analysis of strut-and-tie models using Smooth Evolutionary Structural Optimization. Engineering Structures, 56:1665–1675, 2013a.
  
[7] Liang, Q. Q., Uy, B. e Steven, G.P. Performance-Based Optimization for Strut-and-Tie modeling of Structural Concrete. Journal of Structural Engineering. v.128, n. 6, pp. 815-823, 2002.
+
[7] Liang Q.Q., Uy B., Steven G.P. Performance-based optimization for strut-and-tie modeling of structural concrete. Journal of Structural Engineering, 128(6):815-823, 2002.
 +
</div>

Latest revision as of 12:06, 24 February 2021

Resumo

O Modelo de Bielas e Tirantes pode ser uma excelente alternativa para o dimensionamento de elementos estruturais em concreto armado submetidos a estado plano de tensão e para regiões que apresentem descontinuidade de ordem geométrica ou estática, substituindo procedimentos empíricos por uma metodologia racional de projeto. Para tornar a concepção do modelo menos dependente da experiência do projetista, o presente artigo tem como objetivo aliar a técnica de otimização topológica ESO (Evolutionary Structural Optimization) ao Método dos Elementos Finitos para geração automática dos modelos de bielas e tirantes. O critério de evolução do método de otimização topológica adotado considera a eliminação de elementos menos solicitados em termos de tensão, a partir de uma análise elástico-linear. Nesse contexto, é possível obter soluções otimizadas de problemas complexos envolvendo o concreto estrutural. São apresentados três exemplos numéricos para comprovação e validação das formulações e técnicas implementadas, cujos modelos de bielas e tirantes obtidos apresentam boa concordância em relação às respostas encontradas em trabalhos científicos precursores sobre o tema.

Palavras-chave: Modelo de bielas e tirantes, método dos elementos finitos, otimização topológica, concreto armado

Abstract

The Strut-and-Tie Model can be na excellent alternative for the design of reinforced concrete structural elements submitted to plane stress state and for regions with geometric or static discontinuity, replacing empirical procedures with a rational design methodology. To make the design of the model less dependent on the designer's experience, this article aims to combine the topological optimization technique ESO (Evolutionary Structural Optimization) with the Finite Element Method for the automatic generation of strut-and-tie models. The evolution criterion of the adopted topological optimization method considers the elimination of less stressed elements in terms of stress, from an elastic-linear analysis. In this context, it is possible to obtain optimized solutions of complex problems involving structural concrete. Three numerical examples are presented to prove and validate the formulations and techniques implemented whose strut-and-tie models present a good agreement in relation to the answers found in precursor scientific works on the subject.

Keywords: Strut-and-tie model, finite element method, topological optimization, reinforced concrete

1. Introdução

Para regiões ou elementos estruturais nos quais a Hipótese de Bernoulli não descreva adequadamente o comportamento estrutural ou distribuição de tensões, pode-se recorrer a outras alternativas de dimensionamento, tais como o Método dos Elementos Finitos (MEF) e o Modelo de Bielas e Tirantes (MBT). Nessas regiões, denominadas na literatura de ‘Regiões D’ (Descontinuity), as tensões de cisalhamento são significativas e a distribuição de deformações é não linear. Como exemplo, citam-se elementos estruturais como vigas-parede, consolos, sapatas, nós de pórticos, blocos rígidos de fundação sobre estacas, furos em vigas e dentes Gerber.

Estudos pioneiros envolvendo o MBT tiveram origem no início do século XX. Ritter e Morsch propuseram, a partir de resultados experimentais, a analogia clássica do modelo de treliças para o dimensionamento a cisalhamento de vigas fletidas de concreto armado. Segundo esse modelo, admite-se a substituição da viga original por uma treliça equivalente definida a partir da distribuição de tensões. As barras tracionadas representam campos de tensão de tração (tirantes), enquanto as barras comprimidas representam campos de tensão de compressão (bielas).

Foi a partir dos trabalhos desenvolvidos por Schlaich et al. [1], entretanto, que o tema ganhou forte impulso. Além das vigas inicialmente analisadas, esses pesquisadores estenderam a aplicação do modelo de bielas e tirantes a outros tipos de elementos estruturais. Abordaram temas como procedimentos para definição das regiões com e sem descontinuidade, geração dos modelos de treliças no interior do contínuo de concreto, cálculo dos esforços internos, diretrizes para verificação das tensões nas bielas e regiões nodais e cálculo e detalhamento da armadura necessária.

Entretanto, a não unicidade do modelo topológico, torna a concepção dependente da experiência e da sensibilidade estrutural do projetista para representar o fluxo interno de tensões. Sendo a armadura calculada e distribuída conforme o modelo topológico definido para o elemento estrutural de concreto armado, a correta definição desse não somente gera economia como também está relacionada à segurança.

A otimização topológica (OT) é comumente vista como um método computacional para lançar estruturas a partir da distribuição ótima de material em uma determinada região do espaço. Neste artigo, isso é feito através de uma combinação do Método dos Elementos Finitos, de um modelo para o comportamento do material e de técinas de otimização. Ou seja, o domínio de projeto é discretizado em uma malha refinada de elementos finitos, de modo que se possa analisar seu comportamento e, então, é distribuído material de forma racionalizada através de algoritmos de otimização.

É nesse contexto que se insere o objetivo deste trabalho, isto é, fornecer uma ferramenta eficaz e confiável para a geração automática do modelo de bielas e tirantes via OT, definindo a melhor configuração a ser adotada para a análise. Será adotada uma técnica de otimização de layout, proposta inicialmente por Xie e Steven [2] denominado Otimização Estrutural Evolucionária (Evolutionary Structural Optimization – ESO). A essência do método consiste na remoção gradual de regiões menos solicitadas, com base num critério de penalidade baseado em tensões equivalentes de von Mises. Isto é, elementos com tensões abaixo de um determinado limite são removidos da malha a cada iteração num processo denominado “hard-kill”. Dessa forma, é possível obter uma estrutura ótima para um dado volume remanescente.

Para a realização do objetivo principal pode-se destacar os seguintes objetivos específicos: implementação de um elemento finito triangular de três nós para análise em estado plano de tensão; implementação de uma rotina de otimização topológica dentro de um Programa de Elementos Finitos. A escolha do elemento triangular de três nós para a análise numérica de estruturas em estado plano de tensões, deve-se ao fato de que esse elemento requer uma discretização do contínuo bastante detalhada, permitindo assim definir as regiões de compressão e tração do modelo bielas e tirantes com mais refinamento. Uma vez que, em algumas etapas das análises a técnica de evolução utilizada no processo de otimização topológica na estrutura plana de concreto consiste em eliminar o elemento o da malha de elementos finitos, por isso a exigência de uma malha bem refinada, em consequência um elemento mais pobre em termos das funções de forma.

2. Formulação do elemento finito triangular linear

Para geração do modelo de bielas e tirantes considerando análise linear, foi implementado a técnica ESO que é uma alternativa ao rigor matemático de métodos clássicos de otimização. O algoritmo baseia-se em uma análise via MEF e, como o presente estudo trata de elementos estruturais de concreto armado submetidos a estado plano de tensões, foi necessária, portanto, a implementação de um elemento finito bidimensional, cuja formulações é descrita neste item.

O elemento finito implementado neste trabalho é o elemento triangular de três nós e dois graus de liberdade por nó, denominado Constant Strain Triangle (CST), usado na simulação numérica pelo MEF, baseado em deslocamentos. No caso da análise linear, o material é considerado homogêneo, isotrópico e linear.

Segundo o Princípio dos Trabalhos Virtuais (PTV), aplicando-se um campo de deformação virtual compatível ao elemento triangular tem-se: , isto é, o trabalho virtual externo é igual ao trabalho virtual interno. O trabalho virtual interno, , pode ser escrito como:

(1)

onde é o operador variacional, é o estado tensional real em um ponto qualquer no elemento, e é o estado de deformação virtual em um ponto qualquer no elemento. A partir do tensor de deformação de Green-Lagrange e desprezando-se as tensões no plano de normal na direção , chega-se à Eq. (2) para trabalho virtual interno, onde é a espessura do elemento triangular

(2)

Para o elemento triangular de três nós é usual adotar as funções de interpolação em coordenadas naturais, as quais podem ser escritas como:

(3a)
(3b)
(3c)

onde , , , , , são as coordenadas cartesianas do elemento e sua área.

Definindo o vetor de deslocamentos nodais por e representando as funções de interpolação em coordenadas naturais pelo vetor coluna , definem-se as equações aproximadas dos deslocamentos associados aos deslocamentos nodais q que podem ser escritas na forma matricial conforme a Eq. (4)

(4)

Na Eq. (4), é um vetor coluna nulo com três termos. Sendo os deslocamentos e funções dos deslocamentos nodais, os seus variacionais podem ser escritos a partir da seguinte expressão: , sendo o vetor dos deslocamentos nodais, o operador diferencial e substituído por e Substituindo esses variacionais na Eq. (2) chega-se ao trabalho virtual de um elemento triangular dado pela equação a seguir

(5)

sendo o vetor de forças externas nas direções dos graus de liberdade do elemento triangular, o trabalho virtual externo é dado por , onde é o vetor de forças externas aplicadas diretamente na direção dos graus de liberdade do elemento, e as forças nodais equivalentes obtidas a partir do carregamento externo atuando no contorno do elemento. Da condição tem-se:

(6)

Como a expressão acima deve ser válida para qualquer campo de deslocamento virtual , segue que , onde é o vetor de forças internas definido por:

(7)

A Eq. (7) pode ser escrita na forma matricial como:

(8)

Utilizando o método de Newton-Rapshon na solução do problema é necessária a determinação da derivada dessa expressão em relação aos deslocamentos nodais, obtendo assim a matriz de rigidez tangente. Sendo constante em relação aos deslocamentos nodais, a matriz de rigidez tangente é dada por , que, após manipulações algébricas, pode ser escrita de acordo com a Eq. (9)

(9)

onde, e são as derivadas das funções de forma com relação a e , respectivamente, e são as tensões normais nas direções x e y, respectivamente, é a tensão cisalhante, é o vetor de deslocamento nodal, é a espessura e a área do elemento.

Na Eq. (9), a derivada da tensão normal na direção em relação aos deslocamentos nodais é dada por:

(10)

As derivadas das deformações lineares em relação aos deslocamentos nodais são dadas pelas Eqs. (11) e (12)

(11)
(12)


Substituindo as Eqs. (11) e (12) na derivada da tensão normal na direção x em relação aos deslocamentos nodais (Eq. (10)), tem-se:

(13)

De forma análoga à descrita para tensão normal na direção , pode-se chegar às Eqs. (14) e (15) para as derivadas em relação aos deslocamentos nodais das outras tensões atuantes no elemento

(14)
(15)

A seguir são determinadas as derivadas das funções de forma em relação aos eixos e .

(16)
(17)

Na representação isoparamétrica, as coordenadas cartesianas e são relacionadas com as coordenadas paramétricas e . Então, para mudança do domínio de integração para utiliza-se a relação onde é a matriz jacobiano da transformação das coordenadas e para as coordenadas paramétricas e dada pela equação seguir. Dessa forma,

(18)

3. Otimização estrutural evolucionária

A técnica ESO surge como uma alternativa ao rigor matemático de métodos clássicos de otimização. Apresenta uma base teórica simples cujo fundamento consiste na inserção de vazios na estrutura através da eliminação gradual dos elementos menos solicitados do domínio durante o processo de evolução. Portanto, para se obter a configuração ótima deve-se agregar ao estudo um nível de análise estrutural dependente de um domínio discreto, o que torna o emprego do MEF uma etapa do algoritmo de otimização.

Neste trabalho a representação matemática técnica ESO baseia-se no conceito de tensão, isto é, o nível máximo de tensão na estrutura, obtido por análises via MEF, é tomado como um indicador do nível de eficiência de cada elemento. Elementos com baixo nível de tensão são, portanto, sistematicamente removidos da estrutura. A cada iteração novos elementos ineficientes são eliminados da malha e o procedimento se repete até que o campo de tensão atuante em todo o domínio seja praticamente constante e muito próximo da tensão admissível do material ou que seja atingida a restrição de volume mínimo.

O critério de remoção é feito comparando-se a tensão de von Mises de cada elemento com a tensão de von Mises máxima existente em toda a estrutura. Portanto, no fim de cada iteração todos os elementos que atendem à Eq. (19) serão eliminados. A forma de retirada do elemento ocorre atribuindo-se baixos valores para seu módulo de elasticidade longitudinal (). Desse modo, evita-se o remalhamento da estrutura

(19)

onde:

= tensão de von Mises no elemento analisado;

= razão de rejeição na i-ésima iteração ();

= máxima tensão de von Mises da iteração.

A razão de rejeição é usada para retardar o processo de remoção do elemento. O ciclo de remoção ocorre até que não possam mais ser removidos elementos para um dado valor de . Quando isto ocorre, um estado de equilíbrio é alcançado. O processo evolucionário é redefinido adicionando-se à uma razão de evolução, . A razão de rejeição é atualizada de acordo com a Eq. (20)

(20)

O valor inicial da razão de rejeição () é definido de forma empírica pelo usuário. Entretanto, segundo Querin [3], para garantir melhor convergência, os valores de e devem ser de aproximadamente 1%. O processo se repete enquanto a estrutura não atingir o volume final, , definido pelo usuário, ou seja:

(21)

onde é o volume retirado acumulado até aquela iteração, é o volume final expresso em percentual (por exemplo, implica que a retirada de elementos cessará quando o volume da estrutura atingir 40% do seu volume total inicial) e o volume inicial total da estrutura.

Matematicamente, o ESO pode ser escrito como:

(22)

onde:

= matriz constitutiva do ponto ;

= matriz constitutiva inicial;

= domínio da estrutura;

conjunto dos elementos que não serão removidos.

, conjunto dos elementos que serão removidos.

Portanto, o algoritmo ESO apresenta a seguinte marcha, representada no fluxograma da Figura 1:


1º Passo: discretização do domínio e aplicação das condições de contorno e ações prescritas;

2º Passo: análise da estrutura via MEF e cálculo das tensões principais e tensões de von Mises em cada elemento;

3º Passo: retirar os elementos que satisfaçam a Eq. (19), dentro de um limite pré-definido de volume (p%);

4º Passo: repetir os passos 2 e 3 até que seja atingido o equilíbrio;

5º Passo: acréscimo da razão de rejeição conforme a Eq. (20) e iniciar nova retirada de elementos repetindo os passos 2, 3 e 4.

Draft Ladeira 717005243-image51.png
Figura 1. Fluxograma do algoritmo ESO em nível de tensão

4. Critério de escoamento do von Mises

O critério de remoção é baseado na tensão de von Mises de cada elemento na i-ésima iteração, que em termos de tensões principais pode ser calculada conforme a Eq. (23):

(23)


Para o caso de estado plano de tensão, com

(24)


5. Exemplos numéricos

Neste item são apresentadas algumas aplicações considerando o método ESO implementado. São feitas análises elástico-lineares em estruturas submetidas a estado plano de tensão. O material é considerado homogêneo e isotrópico.

Para cada exemplo são definidas as propriedades mecânicas do material e o domínio inicial de projeto. Além disso, são especificados os seguintes parâmetros: (razão de rejeição), (razão de evolução), (volume final desejado), (volume máximo retirado por iteração), (volume total retirado até uma dada iteração), número da iteração e malha de elementos finitos adotada. Nesse item podem ser feitos os agradecimentos às instituições de fomento à pesquisa, às empresas privadas patrocinadoras do trabalho, aos pesquisadores e profissionais que auxiliaram na execução dos ensaios e fornecimento de dados e materiais.

5.1 Estrutura de Michell

O primeiro exemplo a ser apresentado trata-se de uma viga simplesmente apoiada submetida a uma carga concentrada conforme indicado na Figura 2a, comumente chamada na literatura de estrutura de Michell. A solução analítica é mostrada na Figura 2b. O material adotado foi o aço, cujo módulo de elasticidade GPa, coeficiente de Poisson e espessura igual a 1mm.

Draft Ladeira 717005243-image55.png
Figura 2. (a) Domínio inicial para estrutura de Michell. (b) Solução analítica


O domínio foi discretizado numa malha de elementos triangulares de . O processo evolucionário teve início com uma razão de rejeição () de 1% e uma razão de evolução () de 0,75%. O volume retirado () de 60% do volume inicial e o volume máximo retirado por iteração () de 1,75%. As Figuras 3a, 3b e 3c apresentam a evolução da estrutura, enquanto que a Figura 3d ilustra o modelo de bielas e tirantes obtido, onde as bielas estão representadas em azul e os tirantes em vermelho.

Draft Ladeira 717005243-image56.png
Figura 3. (a) Iteração 22, , . (b) Iteração 69, , . (c) Iteração 175, , . (d) Modelo de bielas e tirantes

5.2 Consolo curto

Este exemplo apresenta um consolo curto projetado para suportar uma carga pontual de 500kN. A Figura 4 traz as dimensões da estrutura em milímetros. O módulo de elasticidade do concreto foi tomado igual a MPa, coeficiente de Poisson e a espessura assumida como sendo igual a 300mm.

Draft Ladeira 717005243-image57.png
Figura 4. Consolo Curto (Liang et al. [4])


Este exemplo foi estudado por Liang et al. [4] e Almeida et al. [5]. A Figura 5 apresenta a evolução da estrutura e o correspondente modelo de bielas e tirantes sugerido pelos autores, no qual as bielas estão representadas por linha pontilhada e os tirantes por linha cheia.

Draft Ladeira 717005243-image58.png
Figura 5. (a), (b), (c) Processo evolucionário. (d) Modelo de bielas e tirantes (Liang et al. [4])


No presente trabalho a estrutura foi modelada usando uma malha com 5664 elementos triangulares de três nós com 25 mm de lado. Para obter a topologia ótima apresentada na Figura 6, os parâmetros adotados foram: , , , . Os campos de compressão (bielas) são representados em azul e os campos de tração (tirantes) em vermelho.

Draft Ladeira 717005243-image59.png
Figura 6. (a) Iteração 25, , . (b) Iteração 31, , . (c) Iteração 47, ,

5.3 Pilar de ponte

Este exemplo foi proposto por Liang et al. [6] e estudado também por Almeida et al. [7]. Trata-se de um pilar de ponte projetado para suportar quatro cargas concentradas de 2750kN transferidas por quatro vigas de aço-concreto. O pilar tem espessura de 1,5m e é admitido como sendo engastado na fundação. A geometria, condições de contorno e ações do problema estão indicados na Figura 7, com dimensões em milímetros e kN.

Draft Ladeira 717005243-image60.png
Figura 7. Domínio de projeto do pilar de ponte (Liang et al. [6])

O domínio foi discretizado numa malha refinada com 12260 elementos finitos triangulares. As propriedades do material isotrópico utilizado são módulo de elasticidade GPa e coeficiente de Poisson .

O processo evolucionário teve início com uma razão de rejeição () igual a 4% e uma razão de evolução () de 2%. O volume final desejado igual a 48% do volume inicial e taxa de retirada de material por iteração () igual a 1,75%. A Figura 8 apresenta o MBT alcançado. As regiões em vermelho e azul indicam, respectivamente, regiões de tração, tirantes, e regiões de compressão, bielas.

Draft Ladeira 717005243-image61.png
Figura 8. (a) Modelo de bielas e tirantes obtido no presente trabalho - Iteração 164, , . (b) Topologia ótima proposta por Liang et al. [6]

6. Conclusões

O objetivo do artigo é a apresentação de uma formulação numérica para verificar o fluxo de tensões em estruturas bidimensionais para geração automática do modelo de bielas e tirantes que auxiliem o projetista na concepção e posterior dimensionamento de estruturas de concreto armado. O algoritmo de otimização topológica ESO foi, portanto, utilizado com esse propósito. O uso do Método dos Elementos Finitos, considerando análise linear elástica, se torna conveniente devido ao critério de remoção adotado para o ESO.

Assim, um domínio estendido é inicialmente discretizado numa malha refinada de elementos finitos e, iterativamente, o algoritmo conduz à configuração ótima que representa o modelo de bielas e tirantes procurado.

Finalmente, pode-se afirmar que os exemplos numéricos apresentados estão de acordo com os resultados encontrados na literatura científica, o que validam a efetividade da formulação implementada.

Para o dimensionamento, alerta-se para o fato de que o modelo de bielas e tirantes deve atender aos procedimentos impostos por códigos normativos, principalmente no que se refere à utilização de armadura mínima, que podem ser maiores que os tirantes concebidos.

Agradecimentos

Os autores agradecem à CAPES, CNPq, FAPEMIG, Fundação Gorceix e PROPEC/UFOP pelo apoio financeiro e suporte necessários para realização desta pesquisa.

Referências

[1] Schlaich J., Schäfer K., Jennewein M. Toward a consistent design of structural concrete. PCI Journal, 32(3):75-150, 1987.

[2] Xie Y.M., Steven G.P. A simple evolutionary procedure for structural optimization. Computers & Structures, 49(5):885-896, 1993.

[3] Querin O.M., Evolutionary structural optimization stress based formulation and implementation. PhD dissertation, University of Sydney, 1997.

[4] Michell A.G.M. The limits of economy of material in frame-structures. Philosophy Magazine, 8:589–597, 1904.

[5] Liang Q.Q., Xie Y.M., Steven G.P. Topology optimization of strut-and-tie models in reinforced concrete structures using an evolutionary procedure. ACI Structural Journal, 97(2):322–330, 2000.

[6] Almeida V.S., Simonetti H.L., Neto L.O. Comparative analysis of strut-and-tie models using Smooth Evolutionary Structural Optimization. Engineering Structures, 56:1665–1675, 2013a.

[7] Liang Q.Q., Uy B., Steven G.P. Performance-based optimization for strut-and-tie modeling of structural concrete. Journal of Structural Engineering, 128(6):815-823, 2002.

Back to Top

Document information

Published on 06/04/20
Accepted on 29/03/20
Submitted on 08/12/19

Volume 36, Issue 2, 2020
DOI: 10.23967/j.rimni.2020.03.008
Licence: CC BY-NC-SA license

Document Score

5

Views 242
Recommendations 2

Share this document

claim authorship

Are you one of the authors of this document?