Determination of structural composite elements strength through the use of the refined plastic‐hinge method

==Resumo==

O presente trabalho tem como finalidade implementar e avaliar uma estratégia generalizada baseada no método da compatibilidade de deformações (MCD) para a obtenção das curvas de resistência e de início de plastificação de estruturas mistas de aço e concreto. Devido à anisotropia do concreto, a definição da resistência de seções transversais que possuem tal material em sua composição não é tão simples como nas estruturas metálicas. Recorre‐se ao método da rótula plástica refinado (MRPR) para realizar, via método dos elementos finitos (MEF), análises avançadas, ou seja, considerando conjuntamente os efeitos não lineares geométricos e físicos. Desse modo, a perda gradual da rigidez no MRPR é determinada através da combinação de esforço normal e momento fletor na seção transversal. Para a solução das equações não lineares global e local (a nível da seção transversal) utiliza‐se o método iterativo de ''Newton‐Rahpson''  acoplado a estratégias para passar por pontos críticos. Apresentam‐se análises numéricas avançadas de sistemas estruturais mistos de aço e concreto através da formulação implementada. Os resultados obtidos são comparados com respostas numéricas e experimentais presentes na literatura.

==Abstract==

This work presents the computational implementation and evaluation of a generalized strategy based on Strain Compatibility Method to obtain the initial and full yield curves of steel‐concrete composite structures. Due to anisotropy of concrete the evaluation of the full yield capacity of cross sections that have this material in its compostion is not so simple as for steel structures. Use is made of the Refined Plastic Hinge Method (RPHM) to perform, via Finite Element Method (FEM), the advanced analysis considering simultaneously second order and inelasticity effects. Thus, the gradual decrease of the stiffness in RPHM is determined by a section combination of axial force and bending moment. The Newton‐Raphson method with path‐following strategies is adopted for solving the nonlinear global and local (in cross section level) equations. Advanced numerical analysis of composite steel‐concrete structural systems are performed through the presented numerical formulation and the results obtained here are compared with experimental and numerical data given in the literature.

==Palavras‐chave==

Análise não linear ; Estruturas mistas ; Método dos elementos finitos ; Método da compatibilidade de deformações ; Método da rótula plástica refinado ; Curvas de resistência

==Keywords==

Nonlinear Analysis ; Composite Structures ; FEM ; Strain Compatibility Method ; Refined Plastic‐Hinge Method ; Yield Curves

==1. Introdução==

Na concepção dos projetos estruturais na engenharia civil, os projetistas iniciam os cálculos baseados em 3 premissas: segurança, tempo (elaboração dos projetos e execução da obra) e economia. Para se atingir o equilíbrio entre essas variáveis há uma gama de materiais e métodos disponíveis. De forma isolada, o aço e o concreto são os materiais mais utilizados na construção civil.

A associação do aço e do concreto visa, fundamentalmente, o melhor aproveitamento físico e mecânico desses materiais [[#bib0135|[1]]] . Dessa forma, as estruturas mistas de aço e concreto apresentam vantagens no que diz respeito à resistência, rigidez, proteção dos elementos metálicos (incêndio e corrosão), custo‐benefício, etc. As estruturas mistas também possuem vantagens durante a execução. Durante a cura do concreto, esse pode ser suportado por perfis metálicos e fôrmas, reduzindo gastos com escoramentos e aumentando o espaço livre de circulação na obra.      

Apesar do uso generalizado de sistemas mistos de aço e concreto, procedimentos numéricos precisos para a sua análise e concepção não evoluíram da mesma forma como para estruturas de aço. Como consequência, as referências normativas possuem falhas quando aplicadas ao projeto seguro e econômico de estruturas mistas [[#bib0140|[2]]] . Landesmann [[#bib0145|[3]]]  cita que, nas análises de estruturas com alta deslocabilidade, normas de projeto [[#bib0150|[4]]] ; [[#bib0155|[5]]]  ;  [[#bib0160|[6]]]  não são capazes de descrever o comportamento real da estrutura. Além disso, Weng e Yen [[#bib0165|[7]]]  demonstram experimentalmente que as diferentes formas como as normas de projeto abordam as estruturas mistas também podem aproximar ou distanciar‐se do comportamento estrutural real.      

Fong [[#bib0170|[8]]]  destaca que muitas normas tendem a recomendar o uso de uma análise de segunda ordem e um método de projeto como forma de se obter resultados mais precisos e de maneira eficiente. Alguns trabalhos recentes tiveram como objetivo o propósito de obter formulações numéricas para análise avançada de estruturas mistas de aço e concreto, com base no método da rótula plástica refinado (MRPR) [[#bib0145|[3]]]  ;  [[#bib0170|[8]]] . Por sua vez, o MRPR [[#bib0175|[9]]]  é descrito como método de análise de plasticidade concentrada, já que a consideração da plastificação do material é avaliada nos pontos nodais de cada elemento finito.      

Para a avaliação da resistência de elementos mistos de aço e concreto através do método da compatibilidade de deformações (MCD), 2 estratégias de solução se destacam. Abordado no presente trabalho, em um raro acoplamento ao MRPR, o método de Newton‐Raphson [[#bib0180|[10]]]  ;  [[#bib0185|[11]]]  usa como variáveis de ajuste a deformação axial no centroide plástico (CP) da seção e a curvatura em um dos eixos principais (problema 2 D). Determina‐se então a deformação axial em cada ponto discreto da seção e utilizam‐se as respectivas relações constitutivas para o somatório de tensões. O procedimento consiste na construção da relação momento‐curvatura (problema não linear), onde, para um dado esforço axial, o ponto limite de momento fletor encontrado reflete a combinação Normal‐Momento responsável pela plastificação da seção. Tal metodologia é abordada na Seção 2 deste trabalho. Adeptos dos métodos quasi‐Newton, alguns autores [[#bib0190|[12]]]  ;  [[#bib0195|[13]]]  utilizam esse procedimento de solução para problemas não lineares acoplado ao método da falsa posição. O processo consiste na determinação da posição da linha neutra e sua orientação, no caso de problemas tridimensionais, que são os parâmetros ajustados no ciclo iterativo. Com isso, podem‐se obter os esforços internos resistentes.      

Tem‐se como objetivo neste trabalho acoplar o MCD ao MRPR e, juntamente com o conceito de rigidez homogeneizada [[#bib0200|[14]]] , proporcionar à plataforma ''Computational System of Advanced Structural Analysis''[[#bib0205|[15]]]  – CS‐ASA – a possibilidade de realizar análises avançadas de estruturas mistas. Além disso, para ampliar a gama de possibilidades de análise, prescrições normativas  [[#bib0155|[5]]]  ;  [[#bib0210|[16]]]  para obtenção da capacidade resistente dos elementos foram também implementadas. Torna‐se assim possível avaliar o comportamento estrutural de elementos mistos em relação às considerações normativas, comparando‐o com resultados obtidos através da formulação generalizada, proposta através do acoplamento MRPR/MCD.      

O problema estrutural não linear global, no contexto do método dos elementos finitos (MEF) e MRPR/MCD, é resolvido por meio da seguinte equação:

{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;" 
|-
| <math>\boldsymbol{\mbox{F}}_i(\boldsymbol{\mbox{U}}\mbox{,}\boldsymbol{\mbox{P}}\mbox{,}\boldsymbol{\mbox{S}}_h)\cong \lambda \boldsymbol{\mbox{F}}_r</math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | ( 1)
|}

sendo '''F'''<sub>''i''</sub>  o vetor de forças internas, que é escrito em função do vetor de deslocamentos nodais '''U''' , esforços internos nos membros, '''P''' , e do parâmetro de avaliação da degradação da rigidez à flexão da seção, '''S'''<sub>''h''</sub> . O vetor de forças externas é definido pelo produto entre o parâmetro de carga, λ, e '''F'''<sub>''r''</sub> , que é o vetor de forças externas de referência.      

Para cada passo de carga, a primeira etapa na solução de (1) envolve a avaliação da chamada solução incremental predita (Δ'''λ'''<sup>0</sup>  e Δ'''U'''<sup>0</sup> ). Para isso, o incremento inicial do parâmetro de carga, Δ'''λ'''<sup>0</sup> , é determinado através de alguma estratégia, que neste trabalho é feita através da técnica do deslocamento generalizado [[#bib0215|[17]]] . Em seguida, pode‐se determinar o incremento inicial dos deslocamentos nodais, Δ'''U'''<sup>0</sup> , através de uma analogia ao sistema '''KU'''  = '''F''' . Porém, em um problema estrutural não linear, dificilmente a solução predita define uma configuração de equilíbrio do sistema. Dessa forma, numa segunda etapa, são necessárias correções nessa solução predita, que neste trabalho é feita através do método de Newton‐Raphson e técnicas de continuidade para obtenção do traçado completo da trajetória de equilíbrio. A correção do parâmetro de carga ao longo do processo iterativo é feita aqui através da norma mínima dos deslocamentos residuais [[#bib0220|[18]]]  acoplada ao método de Newton‐Raphson.      

O presente artigo está organizado em mais 4 seções. Na Seção 2 é apresentado o MCD, juntamente com as relações constitutivas utilizadas e o modelo de discretização da seção transversal. Mais adiante, na Seção 3, a formulação numérica baseada no MRPR [[#bib0175|[9]]]  é detalhada, incluindo o conceito de rigidez homogeneizada, na montagem da matriz de rigidez do elemento. São simulados alguns problemas na Seção 4 para validação da metodologia proposta, que vão desde a obtenção da capacidade resistente de seções transversais isoladas até ao comportamento de sistemas estruturais bidimensionais. A formulação utilizada permite fazer comparações com prescrições normativas tais como o AISC [[#bib0210|[16]]]  e o Eurocode 4 [[#bib0155|[5]]] . Por fim, alguns aspectos importantes observados nos exemplos estudados são levantados na Seção 5.      

==2. Método da compatibilidade de deformações==

Esse método busca acoplar a configuração deformada da seção às relações constitutivas dos materiais que a compõem. Deve‐se definir então uma forma de discretização para que as deformações da seção sejam capturadas da forma mais eficiente possível. A divisão da área em fatias (problema bidimensional), como ilustrado na [[#fig0005|figura 1]] , fornece bons resultados [[#bib0195|[13]]]  ;  [[#bib0200|[14]]] , mas é sensível ao grau de refinamento imposto pelo analista.

<span id='fig0005'></span>

{| style="text-align: center; border: 1px solid #BBB; margin: 1em auto; max-width: 100%;" 
|-
|


[[Image:draft_Content_419355684-1-s2.0-S0213131515000656-gr1.jpg|center|493px|Discretização da seção transversal em fatias: (a) seção transversal, (b) ...]]


|-
| <span style="text-align: center; font-size: 75%;">

Figura 1.

Discretização da seção transversal em fatias: (a) seção transversal, (b) discretização em fatias.

</span>
|}

Neste trabalho, optou‐se por utilizar o método iterativo de Newton‐Raphson com uma estratégia de incrementos constantes para a obtenção da relação momento‐curvatura. Sendo assim, as relações tensão‐deformação dos materiais foram utilizadas. Para o aço, foi adotada uma relação constitutiva elastoplástica perfeita, ilustrada na [[#fig0010|figura 2]] , em que:

{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;" 
|-
| <math>\sigma =\begin{array}{c}
-f_y\mbox{ }\mbox{se}-{\epsilon }_u\leq \epsilon \leq -{\epsilon }_y\\
E_a\epsilon \mbox{ }\mbox{se}-{\epsilon }_y<\epsilon <{\epsilon }_y\\
f_y\mbox{ }\mbox{se}\mbox{ }{\epsilon }_y\leq \epsilon \leq {\epsilon }_u
\end{array}</math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | ( 2)
|}

sendo ''f''<sub>''y''</sub>  a resistência de escoamento, ''E''<sub>''a''</sub>  o módulo de elasticidade do aço, ''ɛ''<sub>''y''</sub>  e ''ɛ''<sub>''u''</sub>  as deformações de escoamento e última, respectivamente.

<span id='fig0010'></span>

{| style="text-align: center; border: 1px solid #BBB; margin: 1em auto; max-width: 100%;" 
|-
|


[[Image:draft_Content_419355684-1-s2.0-S0213131515000656-gr2.jpg|center|228px|Relação constitutiva do aço.]]


|-
| <span style="text-align: center; font-size: 75%;">

Figura 2.

Relação constitutiva do aço.

</span>
|}

Para o concreto, a relação constitutiva é definida através de 2 trechos. Quando comprimido, o diagrama parábola retângulo é adotado ([[#fig0015|fig. 3]] a), limitado pela deformação última do concreto à compressão, ''ɛ''<sub>''cu''</sub> . Há ainda a possibilidade da consideração, ou não, da resistência à tração dada por trecho linear até que a deformação de início de fissuração, ''ɛ''<sub>''tr''</sub> , seja atingida, juntamente com a máxima resistência à tração, ''f''<sub>''cr''</sub>  ([[#fig0015|fig. 3]] b). Dessa forma, tem‐se:

{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;" 
|-
| <math>\sigma =\begin{array}{c}
-f_c\mbox{ }\mbox{se}\mbox{ }{\epsilon }_{cu}\leq \epsilon \leq {\epsilon }_{ci}\\
\left[\frac{2\epsilon }{{\epsilon }_{ci}}-{\left(\frac{\epsilon }{{\epsilon }_{ci}}\right)}^2\right]f_c\mbox{ }\mbox{se}\mbox{ }{\epsilon }_{ci}<\epsilon \leq 0\\
E_{tr}\epsilon \mbox{ }\mbox{se}\mbox{ }0\leq \epsilon \leq {\epsilon }_{cr}\\
f_{cr}\frac{{\alpha }_1{\alpha }_2^2}{1+\sqrt{500\epsilon }}\mbox{ }\mbox{se}\mbox{ }\epsilon >{\epsilon }_{cr}
\end{array}</math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | ( 3)
|}

em que ''f''<sub>''c''</sub>  e ''f''<sub>''cr''</sub>  são as resistência máxima à compressão e à tração do concreto, respectivamente. Vale ressaltar que ''f''<sub>''cr''</sub>  é dada como uma função de ''f''<sub>''c''</sub> , sendo obtida por 1,4(''f''<sub>''c''</sub> /10)<sup>2/3</sup> . ''E''<sub>''c''</sub> , ''E''<sub>''cs''</sub>  e ''E''<sub>''tr''</sub>  são, respectivamente, os módulos de elasticidade tangente inicial, secante e à tração do concreto, ''ɛ''<sub>''ci''</sub>  é a deformação limite do regime parabólico do concreto à compressão, ''α''<sub>1</sub>  é o fator que considera a interação das barras longitudinais com o concreto, e ''α''<sub>2</sub>  é o fator que considera os efeitos de duração da carga e cargas cíclicas. Os parâmetros ''α''<sub>1</sub>  e ''α''<sub>2</sub>  são obtidos de forma experimental, porém Chiorean [[#bib0185|[11]]]  utiliza valores fixos em suas análises (''α''<sub>1</sub>  = 1 e ''α''<sub>2</sub>  = 0,75). Dessa forma esses mesmos valores serão adotados neste trabalho.

<span id='fig0015'></span>

{| style="text-align: center; border: 1px solid #BBB; margin: 1em auto; max-width: 100%;" 
|-
|


[[Image:draft_Content_419355684-1-s2.0-S0213131515000656-gr3.jpg|center|503px|Relação constitutiva do concreto: (a) compressão, (b) tração.]]


|-
| <span style="text-align: center; font-size: 75%;">

Figura 3.

Relação constitutiva do concreto: (a) compressão, (b) tração.

</span>
|}

===2.1. Relação momento‐curvatura===

Uma vez discretizada a seção transversal, 2 variáveis são relevantes para se empregar o MCD acoplado ao método de Newton‐Raphson: a área das fatias e sua respectiva posição. Essa segunda é referenciada ao CP, já que alguns pesquisadores [[#bib0180|[10]]] ; [[#bib0190|[12]]]  ;  [[#bib0225|[19]]]  apontam que dessa forma os problemas de convergência são minimizados.      

A posição do CP da seção é obtida através da seguinte expressão [[#bib0230|[20]]] :

{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;" 
|-
| <math>y_{CP}=\frac{y_cA_c\mbox{ }\gamma \mbox{ }f_{cd}+y_aA_af_{yd}+y_bA_bf_{ybd}}{A_c\mbox{ }\gamma \mbox{ }f_{cd}+A_af_{yd}+A_bf_{ybd}}</math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | ( 4)
|}

com ''A''<sub>''c''</sub> , ''A''<sub>''a''</sub>  e ''A''<sub>''b''</sub>  sendo as áreas de concreto, do perfil metálico e da armadura longitudinal; ''γ''  é um coeficiente que indica se o concreto está confinado (''γ''  = 1,0) ou não (''γ''  = 0,85); ''y''<sub>''c''</sub> , ''y''<sub>''a''</sub>  e ''y''<sub>''b''</sub>  são as coordenadas do centroide, das respectivas áreas no eixo das ordenadas, em relação a um sistema de referência; e ''f''<sub>''cd''</sub> , ''f''<sub>''yd''</sub>  e ''f''<sub>''ybd''</sub>  são a resistência de cálculo do concreto, as resistências de escoamento de cálculo do perfil e das armaduras, respectivamente.      

A [[#fig0020|figura 4]]  ilustra a distribuição de deformações na seção mista para uma combinação de esforços axial e de momento fletor. Note que a deformação axial na ''i'' ‐ésima fatia, ''ɛ''<sub>''i''</sub> , é dada por uma função linear já que foi adotado que os elementos apresentarão interação total. Escreve‐se assim:

{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;" 
|-
| <math>{\epsilon }_i={\epsilon }_0+\phi y_i</math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | ( 5)
|}

em que ''y''<sub>''i''</sub>  é a distância entre os CP da fatia analisada e da seção transversal, ''ɛ''<sub>0</sub>  é a deformação axial no CP da seção e ''ϕ'' , a respectiva curvatura.

<span id='fig0020'></span>

{| style="text-align: center; border: 1px solid #BBB; margin: 1em auto; max-width: 100%;" 
|-
|


[[Image:draft_Content_419355684-1-s2.0-S0213131515000656-gr4.jpg|center|365px|Campo de deformações bidimensional.]]


|-
| <span style="text-align: center; font-size: 75%;">

Figura 4.

Campo de deformações bidimensional.

</span>
|}

Para a notação matricial adotada a seguir, as variáveis ''ɛ''<sub>0</sub>  e ''ϕ''  serão posições do vetor de deformações '''X'''  = [''ɛ''<sub>0</sub>''ϕ'' ]<sup>''T''</sup> . Chiorean [[#bib0185|[11]]]  ressalta que adotando '''X'''  = '''0'''  na primeira iteração, a convergência é atingida mais rapidamente. Numericamente, pode‐se dizer que o equilíbrio da seção é obtido quando a seguinte equação é satisfeita:

{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;" 
|-
| <math>\boldsymbol{\mbox{F}}\left(\boldsymbol{\mbox{X}}\right)=</math><math>\boldsymbol{\mbox{f}}_{\mbox{ext}}-\boldsymbol{\mbox{f}}_{\mbox{int}}\cong \boldsymbol{0}</math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | ( 6)
|}

em que o vetor de forças externas '''f'''<sub>''ext''</sub>  é dado pelo esforço axial, ''N'' , e pelo momento fletor, ''M'' , ou seja:

{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;" 
|-
| <math>\boldsymbol{\mbox{f}}_{ext}=\left[\begin{array}{c}
N\\
M
\end{array}\right]</math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | ( 7)
|}

Já o vetor de forças internas é dado por expressões integrais clássicas para o esforço axial, ''N''<sub>''int''</sub> , e momento fletor, ''M''<sub>''int''</sub> . Uma vez conhecidas as áreas, ''A''<sub>''i''</sub> , e posições, ''y''<sub>''i''</sub> , de cada fatia, a integral passa a ser o somatório descrito como:

{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;" 
|-
| <math>\boldsymbol{\mbox{f}}_{\mbox{int}}=\left[\begin{array}{c}
N_{\mbox{int}}\cong \sum_{i=1}^{n_{fat}}\sigma \left[{\epsilon }_i\left({\epsilon }_0,\phi \right)\right]A_i\\
M_{\mbox{int}}\cong \sum_{i=1}^{n_{fat}}\sigma \left[{\epsilon }_i\left({\epsilon }_0,\phi \right)\right]y_iA_i
\end{array}\right]</math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | ( 8)
|}

Embora seja eficiente iniciar o processo com '''X  = 0                            ''' , a convergência só é atingida na primeira iteração se as forças externas forem nulas. Assim, para a iteração seguinte, ''k''  + 1, o vetor de deformações é dado por:

{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;" 
|-
| <math>\boldsymbol{\mbox{X}}^{k+1}=\boldsymbol{\mbox{X}}^k+</math><math>\boldsymbol{\mbox{F}}^{{'}}{\left(\boldsymbol{\mbox{X}}^k\right)}^{-1}\boldsymbol{\mbox{F}}\left(\boldsymbol{\mbox{X}}^k\right)</math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | ( 9)
|}

na qual '''F’'''  é a matriz de rigidez tangente da seção transversal ou matriz Jacobiana do problema não linear expresso na Equação (6), isto é:

{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;" 
|-
| <math>\boldsymbol{\mbox{F}}^{{'}}=\left(\frac{\partial \boldsymbol{\mbox{F}}}{\partial \boldsymbol{\mbox{X}}}\right)\cong \left[\begin{array}{cc}
\sum_{i=1}^{n_{fat}}E_{T,i}A_i & \sum_{i=1}^{n_{fat}}E_{T,i}A_iy_i\\
\sum_{i=1}^{n_{fat}}E_{T,i}A_iy_i & \sum_{i=1}^{n_{fat}}E_{T,i}A_iy_i^2
\end{array}\right]</math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | ( 10)
|}

em que ''E''<sub>''T,i''</sub>  são os módulos de elasticidade tangentes dos materiais presentes em cada fatia.      

O critério de convergência adotado neste trabalho [[#bib0185|[11]]]  é baseado na razão das normas Euclidianas do vetor de forças desequilibradas, '''F''' , e do vetor de forças externas, '''f'''<sub>''ext''</sub> . Assim, escreve‐se:

{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;" 
|-
| <math>\frac{\left\|\boldsymbol{\mbox{F}}\right\|}{\left\|\boldsymbol{\mbox{f}}_{ext}\right\|}\leq Tol</math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | ( 11)
|}

com ''Tol''  sendo a tolerância, assumida igual a 10<sup>−5</sup> , como proposto por Chiorean [[#bib0185|[11]]] .      

Na [[#tbl0005|tabela 1]]  é detalhado, na forma de algoritmo, o processo de obtenção da relação momento‐curvatura.

<span id='tbl0005'></span>

{| class="wikitable" style="min-width: 60%;margin-left: auto; margin-right: auto;"
|+

Tabela 1.

Estratégia numérica para obtenção da relação momento‐curvatura

|-

! '''1.'''
| Leitura dos dados da seção transversal e dos materiais
|-

! '''2.'''
| Discretização da seção em fatias
|-

! '''3.'''
| Obtenção do centroide plástico CP (Equação 4)
|-

! '''4.'''
| Translação do sistema de referência para o CP
|-

! '''5.'''
| Inicialização: '''X'''  = '''0'''
|-

! '''6.'''
| para cada incremento de momento fletor faça
|-

! '''7.'''
|  Monta‐se '''f'''<sub>''ext''</sub>
|-

! '''8.'''
|  para k = 1, ''nmax''  faça                                                    
|-

! '''9.'''
|   Determina‐se ''ɛ''  (Equação 5)                                                    
|-

! '''10.'''
|   Monta‐se '''f'''<sub>''int''</sub>  (Equação 8)                                                    
|-

! '''11.'''
|   Calcula‐se '''F''' ('''X''' ) (Equação 6)                                                    
|-

! '''12.'''
|   se <nowiki>|</nowiki><nowiki>|</nowiki>'''F''' <nowiki>|</nowiki><nowiki>|</nowiki>/<nowiki>|</nowiki><nowiki>|</nowiki>'''f'''<sub>''ext''</sub> <nowiki>|</nowiki><nowiki>|</nowiki> = ''Tol''  → Pare o processo iterativo e siga para linha 16                                                    
|-

! '''13.'''
|   Monta‐se a matriz de rigidez tangente da seção '''F’'''  (Equação 10)                                                    
|-

! '''14.'''
|   Verificação da singularidade de '''F’'''
|-

! '''15.'''
|   se '''F’'''  é singular → Momento último resistente encontrado–parar processo                                                    
|-

! '''16.'''
|   Correção do vetor de deformações '''X'''  (Equação 9)                                                    
|-

! '''17.'''
|  fim para
|-

! '''18.'''
| fim para
|}

===2.2. Curvas de resistência===

Quando, para um dado esforço axial, é atingido o momento máximo da relação momento‐curvatura, tem‐se a plastificação total da seção. Define‐se então, que esse par de esforços é um ponto da curva de resistência.

A curva de início de plastificação também é obtida da relação momento‐curvatura. Quando a primeira fatia da seção apresenta deformação axial, ''ɛ'' , maior que a deformação de início de escoamento do aço (do perfil e das barras longitudinais) e/ou a deformação de início de plastificação do concreto, a fatia inicia o processo de degradação e consequentemente a seção perde rigidez. O momento responsável por esse fato é tido como momento de início de plastificação.      

===2.3. Curva de início de plastificação simplificada===

Para viabilizar o uso do MRPR com algumas metodologias de normas de projeto foi necessário introduzir um limite para o regime elástico, combinando esforço normal e momento fletor. Fong e Chan [[#bib0235|[21]]]  propuseram que uma deformação de início de degradação para o concreto poderia ser determinada relacionando a lei de Hooke com a relação constitutiva do concreto sob compressão, como ilustrada na [[#fig0025|figura 5]] .

<span id='fig0025'></span>

{| style="text-align: center; border: 1px solid #BBB; margin: 1em auto; max-width: 100%;" 
|-
|


[[Image:draft_Content_419355684-1-s2.0-S0213131515000656-gr5.jpg|center|228px|Relação constitutiva do concreto à compressão com deformação de início de ...]]


|-
| <span style="text-align: center; font-size: 75%;">

Figura 5.

Relação constitutiva do concreto à compressão com deformação de início de plastificação alterada.

</span>
|}

Com a análise do trecho parabólico definido na Equação (3), define‐se que a deformação ''ɛ''<sub>''cin''</sub>  é igual à metade de ''ɛ''<sub>''ci''</sub> , e consequentemente ''f''<sub>''cin''</sub>  é dado por 0,75''f''<sub>''c''</sub> . Dessa forma, os esforços isolados, axiais e momentos fletores, que correspondem ao início de degradação da rigidez à flexão da seção são determinados através da soma de tensões elásticas limites desprezando a resistência à tração do concreto.      

==3. Modelagem numérica via método da rótula plástica refinado==

A definição da capacidade resistente dos elementos estruturais descrita na seção anterior foi introduzida à plataforma computacional CS‐ASA [[#bib0205|[15]]] . Esse programa, que vem sendo desenvolvido ao longo dos anos, já é consolidado para análises não lineares de estruturas de aço, utilizando metodologias para análise da plasticidade concentrada, nos pontos nodais, tal como o MRPR [[#bib0175|[9]]] . As considerações feitas na Seção 2, juntamente com o conceito de rigidez homogeneizada, viabilizam a análise não linear de estruturas mistas de aço e concreto.      

As hipóteses descritas a seguir são consideradas na modelagem das estruturas a serem analisadas:
* interação total entre os materiais presentes na seção transversal;
* todos os elementos são inicialmente retos e prismáticos, e o empenamento é desconsiderado;
* considera‐se que o sistema estrutural é perfeitamente travado no eixo ortogonal ao plano da estrutura, caracterizando um problema bidimensional;
* efeitos de instabilidade localizados são desprezados;
* grandes deslocamentos e rotações de corpo rígido são permitidos;
* os efeitos de deformações causadas pelo cisalhamento serão desprezados.

Na modelagem do sistema estrutural via MEF utiliza‐se o elemento finito de pórtico plano delimitado pelos pontos nodais ''i''  e ''j'' , como mostrado na  [[#fig0030|figura 6]] . Também são ilustrados os esforços internos, ''M''<sub>''i''</sub> , ''M''<sub>''j''</sub>  e ''P'' , e os respectivos graus de liberdade ''θ''<sub>''i''</sub> , ''θ''<sub>''j''</sub>  e ''δ'' , referentes ao sistema corrotacional. São inseridas molas fictícias, de comprimento nulo, nas extremidades do elemento para simular a degradação da rigidez à flexão através do parâmetro ''S''<sub>''h''</sub> .

<span id='fig0030'></span>

{| style="text-align: center; border: 1px solid #BBB; margin: 1em auto; max-width: 100%;" 
|-
|


[[Image:draft_Content_419355684-1-s2.0-S0213131515000656-gr6.jpg|center|304px|Elemento finito de viga‐coluna com molas fictícias referenciado ao sistema ...]]


|-
| <span style="text-align: center; font-size: 75%;">

Figura 6.

Elemento finito de viga‐coluna com molas fictícias referenciado ao sistema corrotacional.

</span>
|}

O parâmetro ''S''<sub>''h''</sub> , que é dado por unidade de rigidez à flexão por radianos, é definido dentro de 3 domínios, como ilustrado na [[#fig0035|figura 7]] . Quando em regime elástico, ''S''<sub>''h''</sub>  é tido como infinito, em termos numéricos adota‐se 10<sup>16</sup> . Tal regime é limitado pela curva de início de plastificação. Uma vez atingido o regime plástico (situação onde as forças internas atingem a curva de resistência), tem‐se a degradação completa da rigidez à flexão. Nesse caso, ''S''<sub>''h''</sub>  é tomado, numericamente, como 10<sup>−10</sup> .

<span id='fig0035'></span>

{| style="text-align: center; border: 1px solid #BBB; margin: 1em auto; max-width: 100%;" 
|-
|


[[Image:draft_Content_419355684-1-s2.0-S0213131515000656-gr7.jpg|center|313px|Diagrama de interação esforço normal momento fletor: definição da rigidez à ...]]


|-
| <span style="text-align: center; font-size: 75%;">

Figura 7.

Diagrama de interação esforço normal momento fletor: definição da rigidez à flexão das molas fictícias.

</span>
|}

Entre as curvas, de início de plastificação e de resistência, considera‐se que a perda de rigidez ocorre de forma gradativa através da seguinte equação:

{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;" 
|-
| <math>S_h=\frac{{\left(EI\right)}_{\mbox{comp}}}{L}\left(\frac{M_{pr}-M}{M-M_{er}}\right)</math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | ( 12)
|}

em que ''L''  é o comprimento do elemento finito; ''M''<sub>''pr''</sub>  e ''M''<sub>''er''</sub>  são, respectivamente, os momentos resistente último e de início de plastificação, definidos na Seção 2; e (''EI'' )<sub>''comp''</sub>  é a rigidez à flexão da seção mista homogeneizada, que para seções duplamente simétricas é dada por:

{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;" 
|-
| <math>{\left(EI\right)}_{\mbox{comp}}=\eta \left[{\left(EI\right)}_\mbox{a}+\right. </math><math>\left. {\left(EI\right)}_\mbox{b}+{\left(\mu EI\right)}_\mbox{c}\right]</math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | ( 13)
|}

em que os subscritos ''a'' , ''b''  e ''c''  referem‐se, respectivamente, ao perfil de aço, às armaduras e ao concreto; ''μ''  o fator de redução da rigidez flexional de concreto devido à fissuração; e ''η''  um coeficiente de redução da rigidez global. Esses últimos são tidos como dados de entrada, já que há uma gama de possibilidades de valores disponíveis na literatura  [[#bib0155|[5]]] ; [[#bib0170|[8]]]  ;  [[#bib0210|[16]]] . No caso de lajes apoiadas sobre perfis I , usa‐se a seguinte expressão [[#bib0240|[22]]] :

{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;" 
|-
| <math>{\left(EI\right)}_{\mbox{comp}}=E_\mbox{a}\left(0,4\mbox{I}_{comp}+\right. </math><math>\left. 0,6{\mbox{I'}}_{comp}\right)</math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | ( 14)
|}

sendo ''I''<sub>''comp''</sub>  e ''I’''<sub>''comp''</sub>  os momentos de inércia da seção mista em regiões de momento positivo e de momento negativo, respectivamente, e ''E''<sub>''a''</sub>  é módulo de elasticidade do aço do perfil.      

Para o elemento ilustrado na [[#fig0030|figura 6]] , a relação força‐deslocamento é expressa por [[#bib0170|[8]]] ; [[#bib0175|[9]]] ; [[#bib0195|[13]]] ; [[#bib0200|[14]]]  ;  [[#bib0205|[15]]] :

{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;" 
|-
| <math>\left\{\begin{array}{c}
\Delta P\\
\Delta M_i\\
\Delta M_j
\end{array}\right\}=\left[\begin{array}{ccc}
{\left(EA\right)}_{comp}/L & 0 & 0\\
0 & S_{hi}-S_{hi}^2(k_{jj}+S_{hj})/\beta  & (S_{hi}S_{hj}k_{ij})/\beta \\
0 & (S_{hi}S_{hj}k_{ji})/\beta  & S_{hj}-S_{hj}^2(k_{ii}+S_{hi})/\beta 
\end{array}\right]\left\{\begin{array}{c}
\Delta \delta \\
\Delta {\theta }_i\\
\Delta {\theta }_j
\end{array}\right\}\mbox{ }\mbox{ou}\mbox{ }\Delta \boldsymbol{\mbox{f}}_c=</math><math>\boldsymbol{\mbox{K}}_c\Delta \boldsymbol{\mbox{u}}_c</math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | ( 15)
|}

em que ''β''  = (''S''<sub>''hi''</sub>  + ''k''<sub>''ii''</sub> )(''S''<sub>''hj''</sub>  + ''k''<sub>''jj''</sub> ) − ''k''<sub>''ji''</sub>''k''<sub>''ij''</sub> . Os subscritos ''i''  e ''j''  indicam as extremidades do elemento e o subscrito ''c''  o sistema de coordenadas utilizado (corrotacional); (''EA'' )<sub>''comp''</sub>  é a rigidez axial da seção transversal mista homogeneizada; Δ''P''  e Δ''M''  são, respectivamente, o esforço normal e o momento fletor incrementais; Δδ e Δθ são os incrementos de deformação axial e rotação nodais.      

Os termos ''k''<sub>''ii''</sub> , ''k''<sub>''ij''</sub> , ''k''<sub>''ji''</sub>  e ''k''<sub>''jj''</sub>  são responsáveis por simular os efeitos de segunda ordem que, neste trabalho, são definidos como [[#bib0215|[17]]] :

{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;" 
|-
| <math>\mbox{k}_{\mbox{ii}}=\mbox{k}_{\mbox{jj}}=\frac{4{\left(\mbox{EI}\right)}_{\mbox{comp}}}{\mbox{L}}\mbox{+}\frac{2PL}{15}\mbox{+}\frac{4\mbox{P}\left(\mbox{I}_{\mbox{comp}}\right)}{\mbox{L}\left(\mbox{A}_{\mbox{comp}}\right)}\mbox{ }e{\mbox{ k}}_{ij}=</math><math>k_{ji}=\frac{2{\left(EI\right)}_{comp}}{L}-\frac{PL}{30}+</math><math>\frac{2P\left(I_{comp}\right)}{L\left(A_{comp}\right)}</math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | ( 16)
|}

sendo ''A''<sub>''comp''</sub>  e ''I''<sub>''comp''</sub>  propriedades geométricas da seção mista homogeneizada.      

Uma vez formadas rótulas plásticas na estrutura, um incremento de carga faz com que as forças internas tendem a violar o domínio resistente naquele ponto. Faz‐se necessária então uma estratégia para que isso não ocorra. Assim, o esforço axial é mantido constante e o momento fletor é reduzido de forma a não violar a curva de resistência. Essa estratégia é denominada ''return mapping''[[#bib0205|[15]]] . Uma vez definida a relação força‐deslocamento no sistema corrotacional, deve‐se levá‐la ao sistema global. Maiores detalhes sobre a transformação de sistemas e o ''return mapping''  podem ser vistos em Silva  [[#bib0205|[15]]] .      

==4. Exemplos numéricos==

Nesta seção, as estratégias numéricas apresentadas são utilizadas visando o estudo de sistemas estruturais mistos de aço e concreto. Os exemplos avaliados são divididos em 4 grupos: relação momento‐curvatura, curvas de resistência, coluna isolada e pórtico simples. Ressalta‐se que os resultados encontrados são comparados com os apresentados na literatura, obtidos de forma numérica e/ou experimental. Além disso, é feito um estudo comparativo entre a metodologia generalizada de obtenção da resistência (MCD) e as descritas em normas de projetos.

===4.1. Relação momento‐curvatura===

O objetivo deste primeiro exemplo é a construção da relação momento‐curvatura de uma seção mista parcialmente revestida, ilustrada na [[#fig0040|figura 8]] a, que foi inicialmente analisada por Caldas [[#bib0180|[10]]] . O aço do perfil I  possui módulo de elasticidade de 205.000 MPa e resistência de escoamento de 250 MPa. O concreto utilizado foi confeccionado para atingir uma resistência de compressão característica, ''f''<sub>''ck''</sub> , de 20 MPa quando o mesmo apresentar deformação entre ‐ 0,0035 e ‐ 0,002. Os coeficientes de minoração da resistência do aço e do concreto são tomados, respectivamente, como 1,10 e 1,40.

<span id='fig0040'></span>

{| style="text-align: center; border: 1px solid #BBB; margin: 1em auto; max-width: 100%;" 
|-
|


[[Image:draft_Content_419355684-1-s2.0-S0213131515000656-gr8.jpg|center|497px|Relações momento‐curvatura da seção mista parcialmente revestida: (a) maior ...]]


|-
| <span style="text-align: center; font-size: 75%;">

Figura 8.

Relações momento‐curvatura da seção mista parcialmente revestida: (a) maior inércia, (b) menor inércia.

</span>
|}

As cargas axiais fixas aplicadas na seção para a construção de cada uma das relações momento‐curvatura são frações da carga axial última resistente denominada ''P''<sub>''u''</sub> , tais como 20, 40 e 60%. Considerando um campo de deformações uniforme igual a ‐ 0,002 em todas as 15 fatias da seção discretizada, tem‐se ''P''<sub>''u''</sub>  = −6.382,177 kN.      

Ressalta‐se que a tração no concreto foi desprezada e foi adotada uma relação constitutiva elástica‐perfeitamente plástica, com ''ɛ''<sub>''u''</sub>  tomado igual a 0,01, para descrever o comportamento do aço.      

Os resultados obtidos são ilustrados nas [[#fig0040|figura 8]] a e b para os eixos de maior e menor inércia, respectivamente.      

Pode‐se observar na [[#fig0040|figura 8]]  uma boa convergência entre os resultados obtidos pelo presente trabalho e por Caldas [[#bib0180|[10]]] . Esse autor utiliza para essa mesma análise a curvatura como variável incremental e, além disso, o critério de parada do processo incremental e iterativo é definido pelos domínios de deformação do concreto. Destaca‐se que a estratégia de verificação da singularidade da matriz de rigidez da seção para finalizar o processo realizado nesta pesquisa é eficiente.      

===4.2. Capacidade resistente de uma seção de concreto armado===

Liu et al. [[#bib0195|[13]]]  estudaram, por meio de métodos ''quasi'' ‐Newton, a capacidade resistente da seção de concreto armado ilustrada na  [[#fig0045|figura 9]] . As propriedades dos materiais que compõem a seção são apresentadas na [[#tbl0010|tabela 2]] . A tração no concreto é desprezada na análise denominada concreto sem tração (CST) e considerada na análise concreto com tração (CCT), e o aço tem seu comportamento descrito por meio de uma relação tensão‐deformação elástica perfeitamente plástica, com ''ɛ''<sub>''u''</sub>  = 0,01.

<span id='fig0045'></span>

{| style="text-align: center; border: 1px solid #BBB; margin: 1em auto; max-width: 100%;" 
|-
|


[[Image:draft_Content_419355684-1-s2.0-S0213131515000656-gr9.jpg|center|188px|Seção de concreto armado (dados em mm).]]


|-
| <span style="text-align: center; font-size: 75%;">

Figura 9.

Seção de concreto armado (dados em mm).

</span>
|}

<span id='tbl0010'></span>

{| class="wikitable" style="min-width: 60%;margin-left: auto; margin-right: auto;"
|+

Tabela 2.

Propriedades dos materiais

|-

! Concreto
| ''f''<sub>''ck''</sub>
| 40 MPa
|-

| 
| ''E''<sub>''c''</sub>
| 35 GPa
|-

| 
| ''γ''<sub>''c''</sub>
| 1,5
|-

! Barras de aço
| ''f''<sub>''ys''</sub>
| 355 MPa
|-

| 
| ''E''<sub>''b''</sub>
| 205 GPa
|-

| 
| ''γ''<sub>''s''</sub>
| 1,15
|-

| 
| ''ϕ''
| 20 mm
|}

O objetivo deste exemplo é avaliar a capacidade resistente da seção de concreto armado no estado limite último (ELU), para a combinação de esforços normal e momento fletor, tanto no eixo de maior quanto no de menor inércia. Verifica‐se na [[#fig0050|figura 10]] a uma boa convergência nos resultados encontrados no presente trabalho e os apresentados por Liu et al. [[#bib0195|[13]]]  por meio da metodologia denominada bloco de tensão equivalente, sendo a resistência à tração do concreto desprezada. Nota‐se que a consideração da tração na relação constitutiva do concreto amplifica a resistência da seção em aproximadamente 2% na região acima do esforço normal de compressão igual a 5.000 kN. Isso se deve ao fato de que nesse trecho da curva de resistência algumas fatias da seção passam a apresentar deformações de tração.

<span id='fig0050'></span>

{| style="text-align: center; border: 1px solid #BBB; margin: 1em auto; max-width: 100%;" 
|-
|


[[Image:draft_Content_419355684-1-s2.0-S0213131515000656-gr10.jpg|center|px|Capacidade resistente da seção de concreto armado: (a) comparação com a ...]]


|-
| <span style="text-align: center; font-size: 75%;">

Figura 10.

Capacidade resistente da seção de concreto armado: (a) comparação com a literatura, (b) Influência da discretização.

</span>
|}

Além da comparação com o resultado presente na literatura, fez‐se um estudo da sensibilidade à discretização da seção transversal. O objetivo desse estudo é avaliar o número mínimo de fatias para que a resposta seja aceitável e, assim, reduzir o número de ''loops''  necessários para a montagem da matriz de rigidez e do vetor de forças internas da seção. Essa avaliação se torna importante na simulação de sistemas estruturais com grande número de graus de liberdade, uma vez que é necessária a obtenção das curvas de resistência em cada ponto nodal da estrutura.      

Foram testadas 3 discretizações para esse estudo: 8, 16 e 25 fatias. Como ilustrado na [[#fig0050|figura 10]] b, pode‐se concluir que a capacidade resistente obtida com o menor índice de refinamento já é suficiente. Outra conclusão que deve ser ressaltada é que índices de refinamento menores que 8 fatias não satisfazem o critério de parada proposto. Isso ocorre pois as fatias intermediárias não atingem valores de deformação suficientes para a singularização da matriz de rigidez. Pode‐se afirmar que para um número de fatias inferior a 8 a curva de resistência seria determinada quando uma fatia apresentar deformação superior à deformação última de um dos materiais componentes da seção. Em outras palavras, parte da seção ainda estaria em regime elástico.      

===4.3. Capacidade resistente de uma seção mista de aço e concreto===

Neste exemplo, uma seção mista totalmente revestida com concreto é estudada ([[#fig0055|fig. 11]] ). As propriedades dos materiais componentes da seção são dadas na [[#tbl0015|tabela 3]]  e os resultados obtidos aqui são confrontados com os fornecidos em Liu et al. [[#bib0195|[13]]] . Além disso, são feitas comparações com as metodologias apresentadas em normas de projeto [[#bib0155|[5]]]  ;  [[#bib0210|[16]]]  e a metodologia simplificada para a curva de início de plastificação apresentada no item 2.3.

<span id='fig0055'></span>

{| style="text-align: center; border: 1px solid #BBB; margin: 1em auto; max-width: 100%;" 
|-
|


[[Image:draft_Content_419355684-1-s2.0-S0213131515000656-gr11.jpg|center|224px|Seção mista (dados em mm).]]


|-
| <span style="text-align: center; font-size: 75%;">

Figura 11.

Seção mista (dados em mm).

</span>
|}

<span id='tbl0015'></span>

{| class="wikitable" style="min-width: 60%;margin-left: auto; margin-right: auto;"
|+

Tabela 3.

Propriedades dos materiais

|-

! Concreto
| ''f''<sub>''ck''</sub>
| 35 MPa
|-

| 
| ''E''<sub>''c''</sub>
| 34 GPa
|-

| 
| ''γ''<sub>''c''</sub>
| 1,5
|-

! Barras de aço
| ''f''<sub>''ys''</sub>
| 460 MPa
|-

| 
| ''E''<sub>''b''</sub>
| 205 GPa
|-

| 
| ''γ''<sub>''s''</sub>
| 1,15
|-

| 
| ''ϕ''
| 32 mm
|-

! Perfil
| ''f''<sub>''y''</sub>
| 355 MPa
|-

| 
| ''E''<sub>''a''</sub>
| 200 GPa
|-

| 
| ''γ''<sub>''s''</sub>
| 1,05
|}

Assim como no exemplo anterior, o MCD acoplado ao método de Newton‐Raphson feito no presente trabalho é satisfatório se comparado com o resultado de Liu et al. [[#bib0195|[13]]] , como visto na [[#fig0060|figura 12]] . Na literatura, a curva de resistência é apresentada apenas no trecho referente à combinação do esforço axial de compressão e de momentos positivos. É importante salientar que a consideração da resistência à tração no concreto não é significativa, pois a análise dos pontos das curvas indica um acréscimo médio de 0,5% na resistência nos trechos em que parte da seção encontra‐se tracionada.

<span id='fig0060'></span>

{| style="text-align: center; border: 1px solid #BBB; margin: 1em auto; max-width: 100%;" 
|-
|


[[Image:draft_Content_419355684-1-s2.0-S0213131515000656-gr12.jpg|center|px|Capacidade resistente da seção mista de aço e concreto: (a) comparação com a ...]]


|-
| <span style="text-align: center; font-size: 75%;">

Figura 12.

Capacidade resistente da seção mista de aço e concreto: (a) comparação com a literatura, (b) influência da discretização.

</span>
|}

Foram realizadas novamente simulações para diferentes discretizações da seção, como visto na [[#fig0060|figura 12]] b. Diferentemente da seção em concreto armado, uma pequena divergência ocorre nos trechos próximos à carga axial nula. Quanto menor o refinamento, nota‐se que a capacidade resistente é levemente superestimada, como mostrado na [[#tbl0020|tabela 4]] .

<span id='tbl0020'></span>

{| class="wikitable" style="min-width: 60%;margin-left: auto; margin-right: auto;"
|+

Tabela 4.

Momentos resistentes para carga axial P = ‐170 kN para cada discretização

|-

! Eixo
! ''M''<sub>''Rd,8''</sub> /''M''<sub>''Rd,16''</sub>
! ''M''<sub>''Rd,8''</sub> /''M''<sub>''Rd,25''</sub>
! ''M''<sub>''Rd,16''</sub> /''M''<sub>''Rd,25''</sub>
|-

! Maior inércia
| 1,003
| 1,005
| 1,002
|-

! Menor inércia
| 1,017
| 1,017
| 1,0
|}

Na [[#fig0065|figura 13]]  são ilustradas as curvas de resistência obtidas com o MCD, pelo EUROCODE 4 [[#bib0155|[5]]]  e pelo AISC [[#bib0210|[16]]] , além da curva de início de plastificação aproximada. Vale ressaltar que ao mesmo tempo que o AISC [[#bib0210|[16]]]  é conservador, o EUROCODE 4 [[#bib0155|[5]]]  apresenta pontos com resistência superior ao MCD. Na mesma figura é mostrada a boa convergência entre a curva de início de plastificação por aproximação linear e a obtida pelo MCD. Assim, pode‐se acoplá‐la aos métodos normativos viabilizando as análises via MRPR.

<span id='fig0065'></span>

{| style="text-align: center; border: 1px solid #BBB; margin: 1em auto; max-width: 100%;" 
|-
|


[[Image:draft_Content_419355684-1-s2.0-S0213131515000656-gr13.jpg|center|373px|Curvas de resistência: MCD x prescrições normativas.]]


|-
| <span style="text-align: center; font-size: 75%;">

Figura 13.

Curvas de resistência: MCD x prescrições normativas.

</span>
|}

===4.4. Colunas mistas isoladas===

Bridge [[#bib0245|[23]]]  realizou testes experimentais em pilares tubulares retangulares e seus resultados têm sido estudados por vários pesquisadores [[#bib0235|[21]]] ; [[#bib0250|[24]]]  ;  [[#bib0255|[25]]] . No presente trabalho, 3 pilares com excentricidade em apenas um dos eixos são simulados, como ilustrado na [[#fig0070|figura 14]] . Os protótipos se diferenciam basicamente pela sua geometria, excentricidade da carga e propriedades dos materiais constituintes. A coluna denominada C1 apresenta menor comprimento, maior seção transversal e materiais mais resistentes do que as demais. Já as colunas C2 e C3 são semelhantes, apresentando diferença apenas na excentricidade da carga vertical aplicada. Esses dados são apresentados na [[#tbl0025|tabela 5]] . Foi adotado o módulo de elasticidade do aço igual a 205.000 MPa. Para o concreto, o módulo de elasticidade secante é facilmente obtido pela relação entre a resistência máxima de compressão e a deformação ''ɛ''<sub>''ci''</sub> .

<span id='fig0070'></span>

{| style="text-align: center; border: 1px solid #BBB; margin: 1em auto; max-width: 100%;" 
|-
|


[[Image:draft_Content_419355684-1-s2.0-S0213131515000656-gr14.jpg|center|271px|Modelo idealizado para análise de pilares isolados com excentricidade.]]


|-
| <span style="text-align: center; font-size: 75%;">

Figura 14.

Modelo idealizado para análise de pilares isolados com excentricidade.

</span>
|}

<span id='tbl0025'></span>

{| class="wikitable" style="min-width: 60%;margin-left: auto; margin-right: auto;"
|+

Tabela 5.

Propriedades geométricas e dos materiais

|-

! Espécime
! ''L'' (mm)                                                    
! ''e'' (mm)                                                    
! ''b'' (mm)                                                    
! ''h'' (mm)                                                    
! ''t'' (mm)                                                    
! ''f''<sub>''y''</sub> (kN/mm<sup>2</sup> )                                                    
! ''f''<sub>''c''</sub> (kN/mm<sup>2</sup> )                                                    
|-

! C1
| 2.130
| 38
| 203,7
| 203,9
| 9,96
| 0,291
| 0,0302
|-

! C2
| 3.050
| 38
| 152,5
| 152,3
| 6,48
| 0,254
| 0,035
|-

! C3
| 3.050
| 64
| 152,5
| 152,3
| 6,48
| 0,254
| 0,035
|}

As trajetórias de equilíbrio obtidas com o CS‐ASA são ilustradas, juntamente com os resultados da literatura, na [[#fig0075|figura 15]] a. As curvas de resistência são determinadas por meio do MCD e de normas [[#bib0155|[5]]]  ;  [[#bib0210|[16]]] . Também são apresentadas as cargas últimas atingidas pelas colunas C1, C2 e C3, experimentalmente [[#bib0245|[23]]]  (''P''<sub>''Test''</sub> ), por Fong e Chan [[#bib0235|[21]]]  (''P''<sub>''FC''</sub> ) e pelo CS‐ASA (''P''<sub>''MCD''</sub> ), na [[#tbl0030|tabela 6]] . Pode‐se verificar a boa concordância entre os resultados obtidos no presente trabalho e os apresentados na literatura.

<span id='fig0075'></span>

{| style="text-align: center; border: 1px solid #BBB; margin: 1em auto; max-width: 100%;" 
|-
|


[[Image:draft_Content_419355684-1-s2.0-S0213131515000656-gr15.jpg|center|536px|Trajetórias de equilíbrio: pilares C1, C2 e C3 (a) comparação com a literatura, ...]]


|-
| <span style="text-align: center; font-size: 75%;">

Figura 15.

Trajetórias de equilíbrio: pilares C1, C2 e C3 (a) comparação com a literatura, (b) comparação MCD x normas de projeto.

</span>
|}

<span id='tbl0030'></span>

{| class="wikitable" style="min-width: 60%;margin-left: auto; margin-right: auto;"
|+

Tabela 6.

Cargas últimas obtidas (em kN) e comparações com a literatura

|-

! Espécime
! ''P''<sub>''Test''</sub>
! ''P''<sub>''FC''</sub>
! ''P''<sub>''MCD''</sub>
! ''P''<sub>''MCD''</sub> /''P''<sub>''Test''</sub>
! ''P''<sub>''MCD''</sub> /''P''<sub>''FC''</sub>
|-

! C1
| 1.956
| 1.938,7
| 1.956,32
| 1,000
| 1,009
|-

! C2
| 680
| 723,5
| 699,18
| 1,028
| 0,966
|-

! C3
| 513
| 548,6
| 544,16
| 1,060
| 0,992
|-

| colspan="4" | Média
| 1,029
| 0,989
|-

| colspan="4" | Desvio padrão
| 0,030
| 0,022
|}

As trajetórias de equilíbrio ilustram que o comportamento do elemento estrutural ao longo do histórico de carregamento, obtido de forma numérica, retrata o comportamento real obtido em laboratório por Bridge [[#bib0245|[23]]] . Em termos estatísticos, verifica‐se uma diferença média entre os resultados numéricos (CS‐ASA) e experimentais de apenas 3%. Destaca‐se também a pequena dispersão dos resultados através do desvio padrão, [[#tbl0030|tabela 6]] .      

Em relação aos resultados de Fong e Chan [[#bib0235|[21]]] , que também utilizaram o MRPR, vale destacar que esses autores empregaram diferentes: formulação geométrica não linear e avaliação da capacidade resistente da seção transversal.      

As análises utilizando as recomendações do AISC [[#bib0210|[16]]]  são conservadoras para estruturas mistas. Como consequência, as cargas máximas encontradas por simulações usando a norma americana, ''P''<sub>''AISC''</sub> , são menores do que as das análises considerando o EUROCODE 4 [[#bib0155|[5]]]  e o MCD. Já a norma europeia, em média, superestima a carga última do elemento estrutural, ''P''<sub>''EC4''</sub> . Porém, em termos estatísticos, apresenta resultados mais próximos ao MCD, como mostrado nas  [[#tbl0035|Tabela 7]]  ;  [[#tbl0040|Tabela 8]] .

<span id='tbl0035'></span>

{| class="wikitable" style="min-width: 60%;margin-left: auto; margin-right: auto;"
|+

Tabela 7.

Cargas últimas obtidas com o EUROCODE 4 [[#bib0155|[5]]]

|-

! Espécime
! ''P''<sub>''EC4''</sub>
! ''P''<sub>''EC4''</sub> /''P''<sub>''Test''</sub>
! ''P''<sub>''EC4''</sub> /''P''<sub>''FC''</sub>
! ''P''<sub>''EC4''</sub> /''P''<sub>''MCD''</sub>
|-

! C1
| 1.912,56
| 0,978
| 0,987
| 0,978
|-

! C2
| 725,18
| 1,067
| 1,002
| 1,037
|-

! C3
| 563,33
| 1,098
| 1,027
| 1,035
|-

| colspan="2" | Média
| 1,048
| 1,005
| 1,017
|-

| colspan="2" | Desvio padrão
| 0,062
| 0,020
| 0,034
|}

<span id='tbl0040'></span>

{| class="wikitable" style="min-width: 60%;margin-left: auto; margin-right: auto;"
|+

Tabela 8.

Cargas últimas obtidas com o AISC [[#bib0210|[16]]]

|-

! Espécime
! ''P''<sub>''AISC''</sub>
! ''P''<sub>''AISC''</sub> /''P''<sub>''Test''</sub>
! ''P''<sub>''AISC''</sub> /''P''<sub>''FC''</sub>
! ''P''<sub>''AISC''</sub> /''P''<sub>''MCD''</sub>
|-

! C1
| 1.685,70
| 0,862
| 0,870
| 0,862
|-

! C2
| 618,39
| 0,909
| 0,855
| 0,884
|-

! C3
| 476,12
| 0,928
| 0,868
| 0,875
|-

| colspan="2" | Média
| 0,900
| 0,864
| 0,874
|-

| colspan="2" | Desvio padrão
| 0,034
| 0,008
| 0,011
|}

Na [[#fig0075|figura 15]] b são apresentadas as trajetórias de equilíbrio obtidas com as curvas de resistência de normas [[#bib0155|[5]]]  ;  [[#bib0210|[16]]] .      

===4.5. Pórtico com pilares metálicos e viga mista===

Huu e Kim [[#bib0130|[26]]]  analisaram o pórtico simples biengastado apresentado na [[#fig0080|figura 16]] . O vão livre entre os pilares metálicos com seção W12 × 50 é de 8 m e a viga possui uma laje apoiada, com espessura de 0,102 m e largura 1,219 m, sobre um perfil W12 × 27. Eles ainda consideraram a interação total entre a laje e o perfil. O sistema estrutural é modelado com concreto de resistência máxima de compressão ''f''<sub>''c''</sub>  igual a 16 MPa e o aço com resistência de escoamento ''f''<sub>''y''</sub>  e módulo de elasticidade iguais a 252,4 MPa e 200.000 MPa, respectivamente. Para esse exemplo, as curvas de resistência são definidas através do MCD. A tração no concreto é desprezada e sua deformação de início de plastificação na compressão é tomada igual a ‐ 0,002. Cargas concentradas de igual intensidade, 150 kN, são aplicadas de forma proporcional e simultaneamente no centro do vão livre da viga e no topo do pilar posicionado à esquerda. Duas discretizações foram testadas: a malha um apresenta uma divisão regular dos membros; já a malha 2 apresenta um refinamento maior, com pontos nodais concentrados em regiões de momentos fletores máximos.

<span id='fig0080'></span>

{| style="text-align: center; border: 1px solid #BBB; margin: 1em auto; max-width: 100%;" 
|-
|


[[Image:draft_Content_419355684-1-s2.0-S0213131515000656-gr16.jpg|center|481px|Pórtico com viga mista: geometria, carregamento e discretização.]]


|-
| <span style="text-align: center; font-size: 75%;">

Figura 16.

Pórtico com viga mista: geometria, carregamento e discretização.

</span>
|}

Na [[#fig0085|figura 17]] , o fator de carga, λ, e o deslocamento lateral no topo do pilar direito são comparados com dados apresentados por Chiorean [[#bib0185|[11]]] . Conclui‐se que as simulações pelo MRPR, utilizando as curvas de resistência obtidas através do MCD, estão condizentes, mais uma vez, com a resposta presente da literatura [[#bib0185|[11]]] . Observa‐se uma perda de rigidez brusca próximo ao fator de carga igual a 0,6. O motivo para tal observação se deve às relações constitutivas adotadas, que implicam uma proximidade das curvas de início de plastificação e de resistência. Nesse ponto, a avaliação indica a formação de uma rótula plástica na extremidade direita da viga (região de momento negativo), local onde o perfil metálico trabalha praticamente isolado. Vale ressaltar que o aumento do refinamento da malha, concentrando pontos nodais em locais estratégicos, não influenciou significativamente na resposta obtida.

<span id='fig0085'></span>

{| style="text-align: center; border: 1px solid #BBB; margin: 1em auto; max-width: 100%;" 
|-
|


[[Image:draft_Content_419355684-1-s2.0-S0213131515000656-gr17.jpg|center|367px|Trajetória de equilíbrio do pórtico analisado.]]


|-
| <span style="text-align: center; font-size: 75%;">

Figura 17.

Trajetória de equilíbrio do pórtico analisado.

</span>
|}

As [[#fig0090|figuras 18]] a e 18b mostram, respectivamente, os diagramas de esforços internos solicitantes e a deformada do pórtico, referentes à carga máxima atingida na análise, λ = 1,085. Pode‐se observar o momento fletor de 166,5 kNm atuando no nó 9, onde ocorre a formação da rótula plástica. Dessa forma, um incremento das forças externas faz com que os esforços solicitantes superem a capacidade resistente do elemento. Para que essa violação não comprometa a análise, foi necessário, portanto, retornar a combinação de esforços (M,N) para o limite resistente da seção através do ''return mapping''[[#bib0205|[15]]] .

<span id='fig0090'></span>

{| style="text-align: center; border: 1px solid #BBB; margin: 1em auto; max-width: 100%;" 
|-
|


[[Image:draft_Content_419355684-1-s2.0-S0213131515000656-gr18.jpg|center|510px|Diagrama de esforços e configuração deformada do pórtico simples quando se ...]]


|-
| <span style="text-align: center; font-size: 75%;">

Figura 18.

Diagrama de esforços e configuração deformada do pórtico simples quando se atinge a carga máxima. (a) Diagramas de esforços internos (λ = 1,085), (b) configuração deformada (λ = 1,085).

</span>
|}

==5. Considerações finais==

No presente trabalho foi apresentada uma metodologia generalizada para a obtenção da capacidade resistente de elementos mistos de aço e concreto. A formulação é baseada no MCD e utiliza o método de Newton‐Raphson para a determinação da relação momento‐curvatura. Sua generalidade permite, que através das relações constitutivas apresentadas, seções puramente de aço ou de concreto armado submetidas à combinação de esforço axial e de momento fletor sejam também modeladas.

Na análise de estruturas mistas, os métodos normativos podem tanto superestimar quanto subestimar a capacidade resistente da seção, uma vez que o MCD é tido como mais próximo à realidade. Observa‐se claramente o conservadorismo do AISC [[#bib0210|[16]]]  ante uma avaliação mais realista proporcionada pelo MCD. Essa característica é oposta à observada na curva de resistência do EUROCODE 4 [[#bib0155|[5]]] .      

O uso do MRPR acoplado a curvas de resistência propostas por normas de projeto foi testado. Fez‐se necessária a introdução de uma curva simplificada para indicar o início da degradação da rigidez à flexão da seção. Além disso, as rigidezes axial e flexional foram determinadas através do conceito de homogeneização da seção transversal. A comparação com resultados experimentais [[#bib0245|[23]]]  e numéricos [[#bib0235|[21]]]  presentes na literatura validam a metodologia proposta, uma vez que os resultados via MCD foram mais próximos aos obtidos em laboratório, quando comparados com curvas de normas de projeto [[#bib0155|[5]]]  ;  [[#bib0210|[16]]] .      

Como etapas futuras desta pesquisa, pretende‐se estudar o comportamento de elementos mistos com interação parcial e pórticos mistos com ligações semirrígidas. Além disso, pretende‐se estender as formulações e simulações ao contexto da análise dinâmica não linear de estruturas mistas.

==Agradecimentos==

Os autores agradecem ao CNPq, CAPES, FAPEMIG e PROPEC/PROPP/UFOP o apoio recebido para a realização deste trabalho.

==Referências==

<ol style='list-style-type: none;margin-left: 0px;'><li><span id='bib0135'></span>
[[#bib0135|[1]]] M.D. Denavit, Characterization of behavior of steel‐concrete composite members and frames with applications for design. Doctoral thesis, Graduate College of the University of Illinois Urbana‐Champaign, Illinois, United States of America, 2012.</li>
<li><span id='bib0140'></span>
[[#bib0140|[2]]] C.K. Iu, M.A. Bradford, W.F. Chen; Second‐order inelastic analysis of composite framed structures based on the refined plastic hinge method; Eng Struct, 31 (2009), pp. 799–813</li>
<li><span id='bib0145'></span>
[[#bib0145|[3]]] A. Landesmann; Plastic‐hinge approach for inelastic analysis of steel‐concrete framed structures; J Constr Steel Res, 66 (2010), pp. 323–334</li>
<li><span id='bib0150'></span>
[[#bib0150|[4]]] AISC; Specification for Structural Steel Buildings; Chicago (2005)</li>
<li><span id='bib0155'></span>
[[#bib0155|[5]]] Eurocode 4, Design of composite steel and concrete structures Part 1: General rules an rules for buildings. Comite Europeen de Normalisation (CEN), Brussels, Belgium, 2004.</li>
<li><span id='bib0160'></span>
[[#bib0160|[6]]] NBR 8800, Projeto de Estruturas de Aço e de Estruturas Mistas de Aço e Concreto de Edifícios. Associação Brasileira de Normas Técnicas (ABNT), Rio de Janeiro, Brasil, 2008.</li>
<li><span id='bib0165'></span>
[[#bib0165|[7]]] C.C. Weng, S.L. Yen; Comparisons of concrete‐encased composite column strength provision of ACI code and AISC specifications; Eng Struct, 24 (2002), pp. 59–72 Elsevier</li>
<li><span id='bib0170'></span>
[[#bib0170|[8]]] M. Fong, Second‐order analysis of imperfect light weight and composite structures. Doctoral thesis, Department of Civil and Structural Engineering, Hong Kong Polytechnic University, Hong Kong, China, 2012.</li>
<li><span id='bib0175'></span>
[[#bib0175|[9]]] S.L. Chan, P. Chui; Non‐linear static and cyclic analysis of steel frames with semi‐rigid connections; Elsevier, Oxford (2000)</li>
<li><span id='bib0180'></span>
[[#bib0180|[10]]] R.B. Caldas, Análise numérica de pilares mistos aço‐concreto, Dissertação de Mestrado, Programa de Pós‐Graduação em Engenharia Civil, Deciv/EM/UFOP, Ouro Preto, MG, Brasil, 2004.</li>
<li><span id='bib0185'></span>
[[#bib0185|[11]]] C.G. Chiorean; A computer method for nonlinear inelastic analysis of 3 D composite steel‐concrete frame structures; Eng Struct, 57 (2013), pp. 125–152</li>
<li><span id='bib0190'></span>
[[#bib0190|[12]]] S. Chen, J.G. Teng, S.L. Chan; Design of biaxially loaded short composite columns of arbitrary section; J Struct Eng, 127 (6) (2001), pp. 678–685</li>
<li><span id='bib0195'></span>
[[#bib0195|[13]]] S.W. Liu, Y.P. Liu, S.L. Chan; Advanced analysis of hybrid steel and concrete frames part 1: Cross‐section analysis technique and second order analysis; J Constr Steel Res, 70 (2012), pp. 326–336</li>
<li><span id='bib0200'></span>
[[#bib0200|[14]]] S.L. Chan, S.W., Liu, Y.P. Liu, Advanced analysis of hybrid frames structures by refined plastic hinge approach, In: Uy, B.; Tao,Z.; Mashiri, F.; Zhu, X.; Mirza, O.; Tan, E. L. (Ed.). Steel and Composite Structures–Proceedings of 4 th International Conference. 2010.</li>
<li><span id='bib0205'></span>
[[#bib0205|[15]]] A.R.D. Silva, Sistema Computacional para a Análise Avançada Estática e Dinâmica de Estruturas Metálicas, Tese de Doutorado, Programa de Pós‐Graduação em Engenharia Civil, Deciv/EM/UFOP, Ouro Preto, MG, Brasil, 2009.</li>
<li><span id='bib0210'></span>
[[#bib0210|[16]]] A.I.S.C.; Specification for Structural Steel Buildings; Chicago (2010)</li>
<li><span id='bib0215'></span>
[[#bib0215|[17]]] Y.B. Yang, S.B. Kuo; Theory & Analysis of Nonlinear Framed Structures; Prentice Hall (1994)</li>
<li><span id='bib0220'></span>
[[#bib0220|[18]]] S.L. Chan; Geometric and material nonlinear analysis of beam‐columns and frames using the minimum residual displacement method; Int J Numer Meth Eng, 26 (1988), pp. 2657–2669</li>
<li><span id='bib0225'></span>
[[#bib0225|[19]]] M.G. Sfakianakis; Biaxial bending with axial force or reinforced, composite and repaired concrete section of arbitrary shape by fiber model and computer graphics; Adv Eng Softw, 33 (2002), pp. 227–242</li>
<li><span id='bib0230'></span>
[[#bib0230|[20]]] K. Roik, R. Bergmann; Design method for composite columns with unsymmetrical croos‐sections; J Contr Steel Res, 15 (1990), pp. 153–168</li>
<li><span id='bib0235'></span>
[[#bib0235|[21]]] M. Fong, S.L. Chan, Advanced analysis of steel concrete composite beam‐columns by refined plastic hinge method, International Journal of Structural Stability and Dynamics, 12 (2012).</li>
<li><span id='bib0240'></span>
[[#bib0240|[22]]] J. Liew, H. Chen, N. Shanmugam; Inelastic analysis of steel frames with composite beams; J Struct Eng, 127 (2001), pp. 194–202</li>
<li><span id='bib0245'></span>
[[#bib0245|[23]]] R. Bridge; Concrete filled steel tubular columns; Civil Eng Transact, 18 (1976), pp. 127–133</li>
<li><span id='bib0250'></span>
[[#bib0250|[24]]] S.W. Liu, Y.P. Liu, S.L. Chan; Advanced analysis of hybrid steel and concrete frames part 2: Refined plastic hinge and advanced analysis; J Constr Steel Res, 70 (2012), pp. 337–349</li>
<li><span id='bib0255'></span>
[[#bib0255|[25]]] R. Gonçalves, J. Carvalho; An efficient geometrically exact beam element for composite columns and its application to concrete encased steel I sections; Eng Struct, 75 (2014), pp. 213–224</li>
<li><span id='bib0130'></span>
[[#bib0130|[26]]] C.N. Huu, S.E. Kim; Pratical nonlinear analysis of steel‐concrete composite frames using fiber‐hinge method; J Constr Steel Res, 74 (2012), pp. 90–97</li>
</ol>