The Idealized Smeared Reinforcement (ISR) method for the optimization of concrete sections: a survey of the state-of-the-art and analysis of potential computational approaches



==“El Perfil idealizado de acero”, un método computacional para el diseño de acero de refuerzo en estructuras de concreto==

==“The Idealized steel profile”, a computational method for the design of reinforcing steel in concrete structures==

<span style="text-align: center; font-size: 75%;">Luis Fernando Verduzco Martínez[1], Alejandro Hernández Martínez, Humberto Esqueda Oliva [3]

</span>

<span style="text-align: center; font-size: 75%;">[1]lf.verduzcomartinez@ugto.mx

[2] alejandro.hernandez@ugtomx.onmicrosoft.com 

[3] esquedah@yahoo.com  <span id="fnc-1"></span>[[#fn-1|<sup>1</sup>]]   Resúmen

Se presenta el desarrollo y modelación matemática y numérica de un nuevo método computacional llamado “El perfil idealizado de acero” para el diseño de acero de refuerzo en estructuras de concreto reforzado. En el presente trabajo se mostrará su aplicación para columnas de concreto prismáticas sólidas de sección rectangular y circular sometidas a flexo-compresión biaxial, de acuerdo a los componentes y criterios de diseño y análisis para este tipo de estructuras, adaptando la modelación numérica para ambas geometrías, así como respetando los criterios que establecen las <span id='citeF-3'></span>[[#cite-3|[3]]].

Se probará al final dicho método para el diseño de diversos modelos experimentales de concreto reforzado bajo ciertas combinaciones de cargas a partir de un análisis estructural previo, haciendo comparaciones de los resultados de acuerdo a eficiencia estructural y geometría en la variación de varios paramétros de diseño y objetivos específicos.

Palabras clave: criterios de diseño, acero de refuerzo, concreto reforzado, formulación matemática, métodos computacionales    Abstract

It is presented the mathematical and numerical devolpement and modelling of a new computational method named “The idealized steel profile” for the design of reinforcing steel in concrete structures. Hereby it will be shown its application for solid prismatic concrete columns of both rectangular and circular cross sections subjected to flexo-compression stresses in two directions, according to the design components and analysis criteria stablished for this type of structures following the <span id='citeF-3'></span>[[#cite-3|[3]]].

At the end such method will be tested for various experimental structural models under certain load combinations from a structural analysis made previously, making comparisons between the results according to certain design criteria and requirements such as structural efficiency, costs and geometry following various specific design goals.

keywords: design criteria, reinforcing steel, reinforced concrete, mathematical formulation, computational methods

<span id="fn-1"></span>
<span style="text-align: center; font-size: 75%;">([[#fnc-1|<sup>1</sup>]]) Universidad de Guanajuato, Aula CIMNE-UG, Guanajuato, Guanajuato, México 2019. </span>

==1 Introducción==

En los proyectos de diseño y análisis estructural de elementos de concreto reforzado se requiere de análisis numéricos que evocan muchas iteraciones para llegar a un resultado que sea tanto estructural como económicamente eficiente; con el uso de métodos numéricos podemos llegar a estas soluciones muy rápidamente.

El presente trabajo está orientado precisamente a usar programación y métodos numéricos en la practica diaria en el diseño de análisis de infraestructura civil, en éste caso se determinará el acero de refuerzo necesario en columnas de concreto para que resistan ciertas condiciones de carga.

En la práctica común, la determinación del acero de refuerzo necesario en elementos de concreto se realiza mediante aproximaciones empíricas basadas en la experiencia; el proceso se vuelve bastante tardado y cansado.

Expondremos pues, el desarrollo de un nuevo método computacional, cuya principal función es la de determinar una configuración de varillas de acero de refuerzo para elementos de concreto sometidos a flexo-compresión, que resulte la más eficiente estructuralmente o la más económica. Para ello el diseñador únicamente tendrá que definir las dimensiones en función de los requerimientos estructurales y/o de espacio. Pueden ser estos elementos columnas, pilas, pilotes o incluso dados de cimentación, con secciones circulares y rectangulares sólidas. Usualmente, por cuestiones constructivas se requieren geometrías de secciones que no sean complejas y que sean uniformes

Como resultados de la programación de este método se obtienen tanto propuestas como recomendaciones de disposiciónes y arreglos del refuerzo, es decir, no solo se obtienen las opciones de número de varillas con su respectivo diámetro, sino también el acomodo que pudieran tener éstas e indicando además la eficiencia crítica estructural para cada arreglo. El programa selecciona ya sea el arreglo más eficiente etructuralmente o el más económico (de acuerdo a las preferencias del constructor) para cada columna. Diversos softwares comerciales de análisis y diseño estructural no hacen esto, de ellos se obitenen únicamente resultados de un análisis con las geometrías y el refuerzo que el diseñador propone inicialmente.

Las aplicaciones y alcances que puede tener éste enfoque de análisis son de gran potencial en el diseño de estructuras de concreto, ya que podría adaptarse a cualquier geometría, cambiando la forma en que se diseñan y analizan elementos de concreto reforzado de hoy en adelante, permitiéndo optimizar los volúmenes de materiales usados, tal que los costos de construcción sean mínimos.

==2 Objetivo==

Desarrollar un método computacional que sea capaz de determinar los posibles arreglos de acero de refuerzo para cualquier número de columnas de concreto de secciones circulares y rectangulares sólidas que se trate, y poder seleccionar la opción más eficiente o económica para cada una, optimizando los volumenes de aceo de refuerzo requeridos, mediante la programación de dicho método en un lenguaje de programación.

==3 Metodología general==

Se presenta la siguiente metodología de manera general para el desarrollo de las funciones del programa de optimización (consideraciones para análisis estructural, optimización, exportación de reportes y resultados).

* Desarrollar la formulación matemática y numérica para la programación en un lenguaje de programación, tanto para secicones de geometría ciruclar y rectangular, de acuerdo a los criterios de diseño y análisis que establecen las <span id='citeF-3'></span>[[#cite-3|[3]]]
* Hacer un estudio de precios unitarios y costos de proveedores de acerero, así como de equipo y mano de obra en general para armar el acero de refuerzo
* Insertar estos costos en el programa al momento del diseño
* Realizar las simulaciones y ejecución de experimentos requeridos para analizar el comportamiento de diferentes modelos estructurales de columnas ante diferentes condiciones de carga, haciendo comparaciones entre los diversos resultados referente a costos finales, eficiencias estructurales y geometrías

==4 Materiales y métodos==

===4.1 Métodos numéricos computacionales con aplicación a diseño de estructuras===

Los métodos numéricos constituyen técnicas mediante las cuales es posible formular problemas matemáticos con resolución mediante operaciones aritméticas. Con el desarrollo de las computadoras digitales eficientes y rápidas, el papel de los métodos numéricos en ingeniería aumento considerablemente.

Antes de las computadoras se gastaba bastante energía en la técnica misma de solución, en lugar de usarla en la definición del problema. Esto debido al trabajo monótono que se requería para obtener resultados numéricos aceptables. Ahora, al usar la potencia de una comutadora se obtienen soluciones directamente, sin tener que recurrir a simplifaciones o soluciones análiticas que pueden resultar complejas. Gracias a esto, es posible dar más importancia ahora a la interpretación de la solución y su incorporación al sistema total del problema.

El desarrollo de los métodos numéricos más famosos por aquellos grandes científicos, matemáticos e ingenieros han servido de inspiración para muchos investigadores en la ingeniería y ciencia más que para desarrollar la aplicación de los existentes a problemas, desarrollar los suyos propios con la potencia de las computadoras.

Se han desarrolaldo desde el uso de las computadoras innumerables métodos computacionales en el análisis y diseño de las estructuras. Tal vez el más famoso y usado de todos es el Método del Elementos Finito <span id='citeF-6'></span>[[#cite-6|[6]]] cuyas aplicaciones van desde el análisis elástico de elementos estructurales, hasta la simulación numérica de fluidos. Desde su primera aparición en 1956 el mismo método ha evulocionado y desarrollado para crear diferentes versiones del mismo dependiendo de su aplicación.

Sin embargo, de manera general, cada vez que se quiere hacer una mejora en un análisis y/o diseño estructural en cuánto a rápidez y aproximaciones de soluciones deseadas se crea un nuevo método numérico computacional. Ya que con solo hecho de implementar la potencia de una computadora para encontrar una solución de un problema se requiere de un ajuste de datos, desarrollo de modelos matemáticos, encontrar coeficientes o iterar operaciones aritméticas, y eso, en sí, lo vuelve un método numérico computacional.

==5 Antecedentes de construcción con concreto reforzado==

Tan solo desde el punto de vista ambiental, en la construcción, la tecnología del concreto armado se antepone a filosofías de sustentabilidad; la manipulación de aceros de construcción ASTM A615 (barras lisas y corrugadas) en forma desmedida desde la concepción de los diseños de refuerzos previos hasta su empleo en la construcción origina altos grados de desperdicios de acero, por efectos de cortes, doblados ineficinetes, o simplemente por una falta de diseño óptimo de este en las estructuras, tal que cumpla con todas las restricciones de eficiencia estructural al mismo tiempo que se use el mínimo posible.

===5.1 El acero de refuerzo en la industria de la construcción===

En obras de concreto armado, se generan entre el 7% al 27% de desperdicios de aceros ASTM A615 en forma de barras corrugadas. Estas barras se fabrican tanto de acero laminado en caliente como acero trabajado en frío. Los diámetros usuales de barras producidas en México varían de <math>\frac{1}{4}pulg.</math>, a <math>1\frac{1}{2} pulg.</math>, aunque algunos productores han fabricado barras corrugadas de <math>\frac{5}{16}pulg</math>, <math>\frac{5}{33}pulg</math> y <math>\frac{3}{16}pulg</math>.

El acero de refuerzo se produce además en otras más diversas formas para todo tipo de elementos de concreto, también las hay barras lisas, barras de acero helicoidales para la forticiación y reforzamiento de rocas, taludes y suelos a manera de perno de fijación. Incluso se fabrican como mallas cuadriculares (transpuestas en doble sentido) principalemtne para estructuras de losas y cimentaciones.

El tener un mejor control en la cantidad de acero que se usa para este tipo de estructuras en la construcción es primordial hoy en día para mitigar los impactos ambientales y económicos a que se enfrenta un país, pues la industria de la construcción siempre será de gran peso e influencia en cualquier región del mundo. Es imprescindible el desarrollo de tecnología que pueda favorecer el control en diseño y construcción de estas estructuras.

==6 Columnas de concreto reforzado==

Es el elemento estructural vertical empleado para sostener la carga de la edificación o estructura principal que se trate. Se utiliza principalemente por la libertad que proporciona para distribuir espacios al tiempo que cumple con la función de soportar el peso de la construcción.

Las formas, armados y las específicaciones de las columnas estarán en razón directa del tipo de esfeurzos y condiciones en general a que estas sean expuestas.

===6.1 Clasificación de columnas de concreto===

Las columnas de concreto reforzado pueden clasificarse como columnas cortas, intermedias y columnas largas.<span id="fnc-2"></span>[[#fn-2|<sup>1</sup>]] Es necesario tener esto en cuenta para entender cómo funciona el programa, ya que el análisis de las columnas preceden del mismo mecanismo de falla de estos tipos de columnas.

Columnas cortas: La carga que pueden soportar está regida por las dimensiones de su sección transversal y por la resistencia de los materiales de que está construida.

Columnas intermedias: La falla es por una combinación de aplastamiento y pandeo.

Columnas largas: La capacidad de carga axial en estás columnas se ve reducida a causa de los momentos secundarios resultates debidos a la deformacion por flexión de la columna.

El presente trabajo trata únicamente con columnas cortas.

<span id="fn-2"></span>
<span style="text-align: center; font-size: 75%;">([[#fnc-2|<sup>1</sup>]]) Jack C. McCormac, Rusell H. Brown, “Diseño de concreto reforzado”, 14va edición, Alfaomega(2015), p.257.</span>

===6.2 Tipos de columnas de concreto===

El tipo de columnas de concreto depende de la forma de la sección  transversal, así como de su tipo de refuerzo <span id="fnc-3"></span>[[#fn-3|<sup>1</sup>]]. Hay columnas de concreto con refuerzo en forma de varillas longitudinales de acero, con restricción lateral (helicoidal o con estribos cerrados).

Generalmente las columnas con estribos tienen una sección transversal rectangular, mientras que las columnas zunchadas suelen tener sección trasnversal circular, pero también pueden fabricarse con secciones rectangulares, octogonales, entre otras formas. Las espirales en comparación de los estribos incrementan en mayor medida la resistencia debido al aumento de la eficiencia del efecto de confinamiento del refuerzo transversal; aunque, en algunos casos, aumentando también los costos de construcción.

<span id="fn-3"></span>
<span style="text-align: center; font-size: 75%;">([[#fnc-3|<sup>1</sup>]]) Jack C. McCormac, Rusell H. Brown, “Diseño de concreto reforzado”, 14th edition, Alfaomega (2015), p. 258.</span>

===6.3 Conceptos generales de diseño===

Es necesario entender las hipótesis empleadas en el diseño de concreto reforzado, para comprender de donde provienen algunos factores, y números incluídos en el análisis.

====6.3.1 Carga axial====

Si se somete un espécimen de concreto simple con relación alto ancho de dos a una carga axial de compresión, la carga máxima se llegará a una deformación unitaria del orden de 0.002. Arbitrariamente se toma como 100 por ciento la resistencia de un espécimen con relación de esbeltez igual a dos. <span id='citeF-4'></span>[[#cite-4|[4]]]

Para relaciones de esbeltez mayores que dos, la resistencia baja, hasta llegar al 85 por ciento, aproximadamente. Por consiguiente, la resistencia de un elemento de concreto simple sujeto a compresión axial puede estimarse como el producto del 85 por ciento del esfuerzo medido en cilindro de control <math>{f'}_{c}</math>, ensayado en las mismas condiciones, multiplicado por el área de la sección transversal del elemento. Este factor de reducción, 0.85, es solo un promedio de resultados de ensayes en miembros colados verticalmente.

Cuando se le adiciona el refuerzo longitudinal a un espécimen de concreto simple y se le añade también el refuerzo transversal necesario para mantener las barras longitudinales en su posición durante el colado, la carga máxima se obtendrá ahora a una deformación unitaria del orden de 0.0021. La falla, en cambio, se produce a una deformación unitaria de 0.003 o 0.004, si el ensaye es de corta duración.

La resistencia adicional sobre la de un prisma de concreto simple debido a la adición del refuerzo longitudinal en compresión se puede estimar como el producto del área del acero por el esfuerzo de fluencia <math>f_{y}</math>

Por lo tanto, la resistencia o carga máxima a compresión que un prisma de concreto con refuerzo longitudinal puede alcanzar esta dada por la siguiente expresión:.

<span id="eq-1"></span>
{| class="formulaSCP" style="width: 100%; text-align: left;" 
|-
| 
{| style="text-align: left; margin:auto;width: 100%;" 
|-
| style="text-align: center;" | <math>  P_{oc} = (0.85)f{'}_{c}(A_{c}-A_{s})+f_{y}(A_{s}) </math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | (1)
|}

Cabe resaltar que se está considerando el área de acero en la contribución de la resistencia del concreto a compresión (primer término de la  Ecuación [[#eq-1|1]], que puede influir considerablemente en los cálculos.

Por otra parte, la carga máxima a tensión que el elemento puede soportar está dada por:

<span id="eq-2"></span>
{| class="formulaSCP" style="width: 100%; text-align: left;" 
|-
| 
{| style="text-align: left; margin:auto;width: 100%;" 
|-
| style="text-align: center;" | <math>  P_{ot} = (f_{y}(A_{s})) </math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | (2)
|}

En la que solo interviene el acero, ya que el concreto se agrieta y no contribuye a la resistencia.

====6.3.2 Flexo-compresión====

====Consideraciones generales de análisis====

<span id="fnc-4"></span>[[#fn-4|<sup>1</sup>]]

<br />a. Se aborda el análisis mediante propuestas de geometría de la sección de la columna que permanecerán fijas, haciendo variar únicamente el área de acero de refuerzo.
<br />b. Un elemento puede alcanzar su resistencia bajo innumerables combinaciones de carga axial y momento flexionante. Estas combinaciones varían desde una capa axial máxima (<math>P_{oc}</math>) sin momento de flexión, hasta un momento (<math>M_{o}</math>) sin fuerza axial.

<br />c. El lugar geométrico de las combinaciones de carga axial y momento flexionante con las que un elemento puede alcanzar su resistencia se representa gráficamente por medio de un ''diagrama de interacción" Figura [[#img-1|1]].
<br />d. Cuando al aumentar la carga externa, el momento y la carga axial crecen en la misma proporción, la historia de carga queda representada por una recta desde el origen, con una pendiente igual al cociente <math display="inline">P/M=1/e</math>. Figura [[#img-1|1]].
<br />e. Existen solo dos modos principales de falla de elementos sujetos a flexo-compresión: falla en compresión y falla en tensión:

<br />1. Falla a compresión: Se produce por aplastamiento del concreto. El acero del lado más comprimido fluye, en tanto que el lado opuesto no fluye en tensión
<br />2. Falla en tensión: Se produce cuando el acero de un lado fluye en tensión antes de que se produzca el aplastamiento del concreto en el lado opuesto más comprimido.<p> 

</p>

<br />f. El diagrama de interacción de un elemento puede obtenerse a partir de las hipótesis descritas para el cálculo de la resistencia de elementos sujetos a flexión pura, considerando que ahora la sumatoria de fuerzas debe ser igual a la carga P aplicada.<p>  

</p>       NOTA: Lo anterior aplica para cualquier geometría.

<div id='img-1'></div>
{| class="floating_imageSCP" style="text-align: center; border: 1px solid #BBB; margin: 1em auto; width: 100%;max-width: 100%;"
|-
|[[Image:Review_975456380192-fig4_1.png|396px|Diagrama de interacción típico. Dibujo propio.]]
|- style="text-align: center; font-size: 75%;"
| colspan="1" | '''Figura 1:''' Diagrama de interacción típico. Dibujo propio.
|}

====Fórmula de Bresler====

Bresler desarrolló una expresión muy simple para columnas rectangulares para calcular los valores máximos de la carga de compresión que actúa a excentricidades <math>e_{x}</math> y <math>e_{y}</math> en secciones rectangulares con refuerzo simétrico. Ecuación [[#eq-3|3]]. Dicha solución surgió ante la necesidad de evitar realizar los cálculos para determinar una superficie de interacción<span id="fnc-5"></span>[[#fn-5|<sup>2</sup>]], reduciendo el problemas a una combinación de soluciones más simples; dos de flexo-compresión en un plano de simetría y una de compresión axial.

<span id="eq-3"></span>
{| class="formulaSCP" style="width: 100%; text-align: left;" 
|-
| 
{| style="text-align: left; margin:auto;width: 100%;" 
|-
| style="text-align: center;" | <math>  \frac{1}{P_{R}}=\frac{1}{P_{rx}}+\frac{1}{P_{ry}}-\frac{1}{P_{oc}} </math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | (3)
|}

Donde: <math>P_{R}</math>=Carga normal resistente sobre la sección transversal del elemento actuando a excentricidades <math>e_{x}</math> y <math>e_{y}</math>.

<math>P_{rx}</math>=Carga normal resistente a una excentricidad de <math>e_{y}</math>.

<math>P_{ry}</math>=Carga normal resistente a una excentricidad de <math>e_{x}</math>.

La Ecuación [[#eq-3|3]] verifica los ensayes disponibles dentro del 20 por ciento de aproximación, y representa una familia de planos que aproximan los puntos de la superficie de interacción. Figura [[#img-2|2]].  <div id='img-2'></div>
{| class="floating_imageSCP" style="text-align: center; border: 1px solid #BBB; margin: 1em auto; width: 100%;max-width: 100%;"
|-
|[[Image:Review_975456380192-fig4.png|420px|Superficie de interacción típica. Tomada de:]]
|- style="text-align: center; font-size: 75%;"
| colspan="1" | '''Figura 2:''' Superficie de interacción típica. Tomada de:
|}
Otro enfoque para analizar columnas sujetas a carga axial y flexión en dos planos Ecuación [[#eq-4|4]] a es partir de la cual se desprende otra expresión simplificada para calcular la eficiencia de un elemento cuando la relación <math>P_{R}/P_{oc}<0.1</math>. Ecuación [[#eq-9|9]] (que es la se presenta en las NTC-2017 [ntc17]) a las cuáles se hará referencia a continuación.

<span id="eq-4"></span>
{| class="formulaSCP" style="width: 100%; text-align: left;" 
|-
| 
{| style="text-align: left; margin:auto;width: 100%;" 
|-
| style="text-align: center;" | <math>   (\frac{P_{u}-P_{nb}}{P_{ot}-P_{nb}})+(\frac{M_{ux}}{M_{nbx}}^{1.5})+(\frac{M_{uy}}{M_{nby}}^{1.5})=1.0 </math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | (4)
|}

Donde: <math>P_{u}</math>=Carga axial nominal aplicada

<math>P_{nb}</math>=Resistencia nominal a carga axial en la condición balanceada.

<math>M_{nbx}</math> y <math>M_{nby}</math>=Momentos nominales resistentes en la condición balanceada alrededor de los ejes  X y Y recpectivamente.

<math>M_{ux}</math> y <math>M_{uy}</math>=Momentos nominales aplicados alrededor de los ejes  X y Y respectivamente.

<span id="fn-4"></span>
<span style="text-align: center; font-size: 75%;">([[#fnc-4|<sup>1</sup>]]) González Cuevas, Fco. Robles Fernández,“Aspectos básicos del Concreto Reforzado”, 4ª edición, Limusa (2005). Capítulo 6, p. 127-155.</span>

<span id="fn-5"></span>
<span style="text-align: center; font-size: 75%;">([[#fnc-5|<sup>2</sup>]]) Es el espacio geométrico de los valores de carga axial que la sección de un elemento es capaz de resistir, en todos sus planos.</span>

===6.4 Normativa===

Anteriormente se presentaron hipótesis generales de diseño con concreto, porque son de tales hipótesis en que los reglamentos de construcción se basan. Para desarrollar éste proyecto se hizo referencia a las Normas Técnicas complementarias para el diseño y construcción de estructuras de concreto del Reglamento de Construcciones de la Ciudad de México (NTC-2017), que aunque se basan en las hipótesis anteriores no presentan tal cual los mismos criterios de diseño, y que son cuestionados mucho por la sociedad ingenieril. Veamos que proponen estas normas.

====6.4.1 Carga axial====

En las NTC-2017 se considera necesario hacer una modificación en el valor de <math>f{'}_{c}</math> mediante Factores de carga <math>F_{c}</math> y Factores de Resistencia <math>F_{R}</math><span id="fnc-6"></span>[[#fn-6|<sup>1</sup>]]. En el apartado de concreto se específica que para el cálculo de resistencias se utilice una resistencia reducida a la compresión del concreto denominada <math>f^{''}c</math>. Cuyo valor es:

<span id="eq-5"></span>
{| class="formulaSCP" style="width: 100%; text-align: left;" 
|-
| 
{| style="text-align: left; margin:auto;width: 100%;" 
|-
| style="text-align: center;" | <math>  f{''}c=0.85(f{'}c) </math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | (5)
|}

De modo que la ecuación Ecuación [ [[#eq-1|1]]] se transforma a:

<span id="eq-6"></span>
{| class="formulaSCP" style="width: 100%; text-align: left;" 
|-
| 
{| style="text-align: left; margin:auto;width: 100%;" 
|-
| style="text-align: center;" | <math>  P_{oc}=F_{R}((f{''}c)(A_{c}-A_{s})+f_{y}(A_{s})) </math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | (6)
|}

De modo que la ecuación Ecuación [ [[#eq-2|2]]] se transforma a:

<span id="eq-7"></span>
{| class="formulaSCP" style="width: 100%; text-align: left;" 
|-
| 
{| style="text-align: left; margin:auto;width: 100%;" 
|-
| style="text-align: center;" | <math>  P_{ot}=F_{R}(f_{y})(A_{s}) </math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | (7)
|}

Para una mejor interpretación de los factores mecionados se puede hacer referencia a la Figura [[#img-3|3]], que es el bloque de transformación equivalente de esfuerzos de las hipótesis de las NTC-2017.

<div id='img-3'></div>
{| class="floating_imageSCP" style="text-align: center; border: 1px solid #BBB; margin: 1em auto; width: 100%;max-width: 100%;"
|-
|[[Image:Review_975456380192-Fig5.png|576px|Hipótesis de las NTC-2017 sobre la distribución de deformaciones y esfuerzos en la zona de compresión. Adaptada de: <span id='citeF-3'></span>[[#cite-3|[3]]]]]
|- style="text-align: center; font-size: 75%;"
| colspan="1" | '''Figura 3:''' Hipótesis de las NTC-2017 sobre la distribución de deformaciones y esfuerzos en la zona de compresión. Adaptada de: <span id='citeF-3'></span>[[#cite-3|[3]]]
|}
Donde: <math> C=ab(f{''}_{c}) </math> <math> 0.65<=(\beta _{1}=1.05-\frac{f{'}_{c}}{1400})<=0.85 </math>

====6.4.2 Compresión y flexión en dos direcciones====

<span id="fnc-7"></span>[[#fn-7|<sup>2</sup>]]   Se usará la expresión siguiente (que es derivada de la fórmula de Bresler) para el cálculo de la carga máxima actuante en la sección a cierta excentricidad.

<span id="eq-8"></span>
{| class="formulaSCP" style="width: 100%; text-align: left;" 
|-
| 
{| style="text-align: left; margin:auto;width: 100%;" 
|-
| style="text-align: center;" | <math>   P_{R}=\frac{1}{\frac{1}{P_{rx}}+\frac{1}{P_{ry}}-\frac{1}{P_{oc}}} </math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | (8)
|}

Para valores de <math>\frac{P_{R}}{P_{oc}}<0.1</math>, se usará la expresión siguiente:

<span id="eq-9"></span>
{| class="formulaSCP" style="width: 100%; text-align: left;" 
|-
| 
{| style="text-align: left; margin:auto;width: 100%;" 
|-
| style="text-align: center;" | <math>  \frac{M_{ux}}{M_{rx}}+\frac{M_{uy}}{M_{ry}}<=1.0 </math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | (9)
|}

Donde: <math>M_{ux}</math> y <math>M_{uy}</math> son los momentos actuantes alrededor de los ejes X, Y, respectivamente.

<math>M_{rx}</math> y <math>M_{ry}</math> son los momentos resistentes de diseño alrededor de los mismos ejes.

<span id="fn-6"></span>
<span style="text-align: center; font-size: 75%;">([[#fnc-6|<sup>1</sup>]]) El factor de reducción <math>F_{R}</math>, que para el caso de columnas con carga axial es de 0.8 <span id='citeF-3'></span>[[#cite-3|[3]]]</span>

<span id="fn-7"></span>
<span style="text-align: center; font-size: 75%;">([[#fnc-7|<sup>2</sup>]]) González Cuevas, Fco. Robles Fernández, “Aspectos básicos del Concreto Reforzado”, 4ª edición, Limusa (2005). Capítulo 6, p. 148.</span>

==7 El método del perfil idealizado de acero==

Empleando las hipótesis de diseño anteriormente descritas se ha idealizado el acero de refuerzo de un elemento como se muestra. Figura [[#img-4|4]]   <div id='img-4'></div>
{| class="floating_imageSCP" style="text-align: center; border: 1px solid #BBB; margin: 1em auto; width: 100%;max-width: 100%;"
|-
|[[Image:Review_975456380192-Fig6.png|193px|Sección rectangular idealizada de concreto reforzado. Dibujo propio]]
|- style="text-align: center; font-size: 75%;"
| colspan="1" | '''Figura 4:''' Sección rectangular idealizada de concreto reforzado. Dibujo propio
|}
Para estas secciones se desarrolló un método puramente matemático con ecuaciones que definieran la resistencia de la columna de acuerdo a la ubicación del eje neutro de esfuerzos en la sección siguiendo paramétros de geometrías que resulta un tanto complejo para programar. Por ello también se desarrolló un método numérico que resulta más sencillo de programar y se meustra a continuación.

Tomando como referencia la siguiente figura Figura [[#img-5|5]]:  <div id='img-5'></div>
{| class="floating_imageSCP" style="text-align: center; border: 1px solid #BBB; margin: 1em auto; width: 100%;max-width: 100%;"
|-
|[[Image:Review_975456380192-perfilacero1.png|308px|Sección rectangular idealizada de concreto reforzado. Dibujo propio]]
|- style="text-align: center; font-size: 75%;"
| colspan="1" | '''Figura 5:''' Sección rectangular idealizada de concreto reforzado. Dibujo propio
|}
Donde: <math>t</math> = espesor de la sección transversal del acero de refuerzo idealizado como un refuerzo constante alrededor de la sección.

===7.1 Secciones rectangulares-modelo numérico===

Se discretizará el perfil idealizado de acero como se muestra a continuación Figura [[#img-6|6]]:  <div id='img-6'></div>
{| class="floating_imageSCP" style="text-align: center; border: 1px solid #BBB; margin: 1em auto; width: 100%;max-width: 100%;"
|-
|[[Image:Review_975456380192-numerico_rec_modelo.png|257px|Modelo geométrico para la metodología numérica de análisis:  Dibujo propio.]]
|- style="text-align: center; font-size: 75%;"
| colspan="1" | '''Figura 6:''' Modelo geométrico para la metodología numérica de análisis:  Dibujo propio.
|}
Donde: <math>dA_{ac}=dL(t)</math>

<math>L=h-2rec</math>

<math>dL=\frac{L}{nElementos}</math>

<math>nElementos=</math> Número de elementos

====7.1.1 Cálculo de la resistencia del acero====

Metodología      Para la determinación de la resistencia de carga axial y flexión se sigue la siguiente metodología:

* Determinar la distancia de cada segmento de perfil con respecto a la fibra superior más alejada de la sección transversal del elemento de concreto <math display="inline">d</math>.<p>   Para las partes superior e inferior del perfil se tiene:

<span id="eq-10"></span>
{| class="formulaSCP" style="width: 100%; text-align: left;" 
|-
| 
{| style="text-align: left; margin:auto;width: 100%;" 
|-
| style="text-align: center;" | <math>  F_{a}=\varepsilon _{a}E_{ac}t(b-2rec) </math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | (10)
|}

<span id="eq-11"></span>
{| class="formulaSCP" style="width: 100%; text-align: left;" 
|-
| 
{| style="text-align: left; margin:auto;width: 100%;" 
|-
| style="text-align: center;" | <math>  F_{b}=\varepsilon _{b}E_{ac}t(b-2rec) </math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | (11)
|}

{| class="formulaSCP" style="width: 100%; text-align: left;" 
|-
| 
{| style="text-align: left; margin:auto;width: 100%;" 
|-
| style="text-align: center;" | <math> M_{a,b}=F_{a}(d_{1}-\frac{1}{2}h)+F_{b}(d_{2}-\frac{1}{2}h) </math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | (12)
|}

Donde: <math>\varepsilon _{a}=(1-\frac{d_{1}}{c})(0.003)</math></p><p> <math>\varepsilon _{b}=(\frac{d_{2}}{c}-1)(0.003)</math></p><p> <math>d_{1}=rec+\frac{1}{2}(t)</math></p><p> <math>d_{2}=h-rec-\frac{1}{2}(t)</math></p><p> <math>E_{ac}=\frac{fy}{\varepsilon _{y}}=\frac{4200 \frac{kg}{cm^2}}{0.0021}</math></p><p> <math>E_{ac}=2e10 \frac{kg}{cm^{2}}</math>       Mientras que para los costados del perfil se tiene:

{| class="formulaSCP" style="width: 100%; text-align: left;" 
|-
| 
{| style="text-align: left; margin:auto;width: 100%;" 
|-
| style="text-align: center;" | <math> F_{R}=2\sum _{i=1}^{nElementos}E_{ac}\varepsilon (dA_{ac}) </math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | (13)
|}

{| class="formulaSCP" style="width: 100%; text-align: left;" 
|-
| 
{| style="text-align: left; margin:auto;width: 100%;" 
|-
| style="text-align: center;" | <math> M_{R}=\sum _{i=1}^{nElementos}-2E_{ac}\varepsilon (dA_{ac})(\frac{1}{2}h-d) </math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | (14)
|}

Donde: <math>t</math> = Espesor del perfil</p><p> <math>rec</math> = recubrimiento</p><p> <math>h</math> = peralte de la sección</p><p> <math>\varepsilon _{a}=(1-\frac{d}{c})(0.003)</math></p><p>       </p>
* A partir de ahí se cálcula <math display="inline">\varepsilon </math> en función de c, cuidando únicamente la siguiente restricción <math display="inline">[-0.0021<\varepsilon <0.0021]</math>.<p>  </p>
* Y finalmente calcular la fuerza y momento resultantes resistentes:<p>

{| class="formulaSCP" style="width: 100%; text-align: left;" 
|-
| 
{| style="text-align: left; margin:auto;width: 100%;" 
|-
| style="text-align: center;" | <math> F_{R}=F_{a}+F_{b}+2\sum _{i=1}^{nElementos}E_{ac}\varepsilon (dA_{ac}) </math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | (15)
|}

{| class="formulaSCP" style="width: 100%; text-align: left;" 
|-
| 
{| style="text-align: left; margin:auto;width: 100%;" 
|-
| style="text-align: center;" | <math> M_{R}=F_{a}(d_{1}-\frac{1}{2}h)+F_{b}(d_{2}-\frac{1}{2}h)+\sum _{i=1}^{nElementos}-2E_{ac}\varepsilon (dA_{ac})(\frac{1}{2}h-d) </math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | (16)
|}

</p>

====7.1.2 Cálculo de la resistencia para el concreto====

Lo que se pretende aquí es encontrar el área neta del concreto bajo esfuerzos para cualquier valor en la profundidad del eje neutro Figura [[#img-8|8]]. Al ser una sección rectangular, el análisis se simplifica a un modelo matemático puro, ya que el área variará uniformemente a lo largo de todo el peralte de la sección, ya sea respecto al eje X o Y. A partir de esa área neta se determinará el volumen del bloque de esfuerzos-deformaciones resultante de compresión en unidades de Fuerza, y posteriormente el momento de flexión resultante respecto al eje del centroide de la sección como sigue:

<div id='img-7'></div>
{| class="floating_imageSCP" style="text-align: center; border: 1px solid #BBB; margin: 1em auto; width: 100%;max-width: 100%;"
|-
|[[Image:Review_975456380192-fig21.png|465px|Diagrama de bloque equivalente de esuferzos de compresión en el concreto. Dibujo propio.]]
|- style="text-align: center; font-size: 75%;"
| colspan="1" | '''Figura 7:''' Diagrama de bloque equivalente de esuferzos de compresión en el concreto. Dibujo propio.
|}

<span id="eq-17"></span>
{| class="formulaSCP" style="width: 100%; text-align: left;" 
|-
| 
{| style="text-align: left; margin:auto;width: 100%;" 
|-
| style="text-align: center;" | <math>  C=-abf{''}c </math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | (17)
|}

<span id="eq-18"></span>
{| class="formulaSCP" style="width: 100%; text-align: left;" 
|-
| 
{| style="text-align: left; margin:auto;width: 100%;" 
|-
| style="text-align: center;" | <math>  M_{c}=-C(\frac{h}{2}-\frac{a}{2}) </math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | (18)
|}

<div id='img-8'></div>
{| class="floating_imageSCP" style="text-align: center; border: 1px solid #BBB; margin: 1em auto; width: 100%;max-width: 100%;"
|-
|[[Image:Review_975456380192-fig22.png|398px|Diagrama tridimensional del bloque equivalente de esfuerzos de compresión en el concreto. Dibujo propio.]]
|- style="text-align: center; font-size: 75%;"
| colspan="1" | '''Figura 8:''' Diagrama tridimensional del bloque equivalente de esfuerzos de compresión en el concreto. Dibujo propio.
|}

==8 Columnas circulares==

Análogo a la idealización del acero de refuerzo de columnas rectangulares, se ha idealizado el acero de refuerzo en éste tipo de columnas como se muestra en la siguiente figura Figura [[#img-9|9]].   <div id='img-9'></div>
{| class="floating_imageSCP" style="text-align: center; border: 1px solid #BBB; margin: 1em auto; width: 100%;max-width: 100%;"
|-
|[[Image:Review_975456380192-fig23.png|241px|Idealización de una sección circular de concreto reforzado. Dibujo propio.]]
|- style="text-align: center; font-size: 75%;"
| colspan="1" | '''Figura 9:''' Idealización de una sección circular de concreto reforzado. Dibujo propio.
|}
Donde: <math>t</math> = Espesor del perfil

<math>rec</math> = recubrimiento

<math>d_{ma}</math> = diámetro de la sección del elemento       Nota: El recubrimiento se mide desde el borde de la sección de la columna al centro del espesor del refuerzo ficticio.

A diferencia de la geometría rectangular, para la ciruclar resultó bastante complejo analizar la distribución de esfuerzos mediante el enfoque de análisis planteado originalmente (determinando la geometría de las configuraciones de estados de deformaciones, por lo que simplificó con el siguiente planteamiento:

===8.1 Cálculo de la resistencia del acero===

Lo que se hará es tratar el perfil de acero como una línea circular y subdividirla en pequeños segmentos con espesor <math>t</math>, Figura [[#img-10|10]], aunque en realidad cada segmento se tratará como un punto al que le corresponderá una cierta condición de esfuerzos dependiendo de su ubicación en el plano cartesiano de referencia y del caso del cual se trate.   <div id='img-10'></div>
{| class="floating_imageSCP" style="text-align: center; border: 1px solid #BBB; margin: 1em auto; width: 100%;max-width: 100%;"
|-
|[[Image:Review_975456380192-fig27.png|499px|Esquema de análisis para las condiciones de esfuerzos en el acero de refuerzo. Dibujo propio.]]
|- style="text-align: center; font-size: 75%;"
| colspan="1" | '''Figura 10:''' Esquema de análisis para las condiciones de esfuerzos en el acero de refuerzo. Dibujo propio.
|}
Donde: <math>dA_{ac}=dL(t)</math>

<math>L=\Pi (d_{ma}-2rec)</math>

<math>dL=\frac{L}{n}</math>

<math>n=</math> Número de elementos 

====8.1.1 Metodología====

Para la determinación de la resistencia de carga axial y flexión se sigue la siguiente metodología:

* Determinar la distancia de cada segmento de perfil con respecto a la fibra más alejada de la sección transversal del elemento de concreto <math display="inline">d</math>.<p>  </p>
* A partir de ahí se cálcula <math display="inline">\varepsilon </math> en función de c, cuidando únicamente la siguiente restricción <math display="inline">[-0.0021<\varepsilon <0.0021]</math>.<p>  </p>
* Y finalmente calcular la fuerza y momento resultantes resistentes:<p>

{| class="formulaSCP" style="width: 100%; text-align: left;" 
|-
| 
{| style="text-align: left; margin:auto;width: 100%;" 
|-
| style="text-align: center;" | <math> F_{R}=\sum _{i=1}^{nElementos}E_{ac}\varepsilon (dA_{ac}) </math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | (19)
|}

{| class="formulaSCP" style="width: 100%; text-align: left;" 
|-
| 
{| style="text-align: left; margin:auto;width: 100%;" 
|-
| style="text-align: center;" | <math> M_{R}=\sum _{i=1}^{nElementos}-E_{ac}\varepsilon (dA_{ac})(\frac{1}{2}d_{ma}-d) </math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | (20)
|}

</p>

===8.2 Casos de análisis para el concreto===

A diferencia de las columnas de sección rectangular, en las circulares el diferencial de área al cambiar el eje neutro no será uniforme.

También podremos transformar ésta área de compresión en un bloque equivalente de esfuerzos, y para poder hacer ésto será necesario cálcular el centroide del área real de compresión en cada posición del eje neutro y a partir de ahí calcular el momento de flexión resistente de la sección. Figura [[#img-11|11]].  <div id='img-11'></div>
{| class="floating_imageSCP" style="text-align: center; border: 1px solid #BBB; margin: 1em auto; width: 100%;max-width: 100%;"
|-
|[[Image:Review_975456380192-fig47.png|539px|Diagrama de distribución de esfuerzos para el concreto. Dibujo propio.]]
|- style="text-align: center; font-size: 75%;"
| colspan="1" | '''Figura 11:''' Diagrama de distribución de esfuerzos para el concreto. Dibujo propio.
|}
Donde: <math>a=\beta (c)</math>

<math>f{''}_{c}=F_{R}(0.85)f^{'}c</math>      Para calcular ésta fuerza de compresión tomaremos el siguiente sistema de referencia: Figura [[#img-12|12]], para definir integrales para calcular el área efectiva de compresión del concreto en cualquier posición del eje neutro.  <div id='img-12'></div>
{| class="floating_imageSCP" style="text-align: center; border: 1px solid #BBB; margin: 1em auto; width: 100%;max-width: 100%;"
|-
|[[Image:Review_975456380192-fig48.png|567px|Sistema de referencia cartesiano para el análisis de la distribución de esfuerzos de compresión en el concreto. Dibujo propio.]]
|- style="text-align: center; font-size: 75%;"
| colspan="1" | '''Figura 12:''' Sistema de referencia cartesiano para el análisis de la distribución de esfuerzos de compresión en el concreto. Dibujo propio.
|}
Subdividiremos el análisis en dos casos para simplificación:

Cuando <math>a<\frac{1}{2}d_{ma}</math>

{| class="formulaSCP" style="width: 100%; text-align: left;" 
|-
| 
{| style="text-align: left; margin:auto;width: 100%;" 
|-
| style="text-align: center;" | <math> F_{co}=f{''}c \int _{-z(y_{co}=a)}^{z(y_{co}=a)}(a-y_{co})\mathrm{d}z </math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | (21)
|}

{| class="formulaSCP" style="width: 100%; text-align: left;" 
|-
| 
{| style="text-align: left; margin:auto;width: 100%;" 
|-
| style="text-align: center;" | <math>\begin{array}{l} F_{co}\cdot Y_{co}=f{''}c \int _{-z(y_{co}=a)}^{z(y_{co}=a)}(a-y_{co})...\\ (\frac{1}{2}d_{ma}-\frac{1}{2}(a+y_{co})))\mathrm{d}z \end{array}</math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | (22)
|}

Cuando <math>a>=\frac{1}{2}d_{ma}</math>

{| class="formulaSCP" style="width: 100%; text-align: left;" 
|-
| 
{| style="text-align: left; margin:auto;width: 100%;" 
|-
| style="text-align: center;" | <math>\begin{array}{l} F_{co}=f{''}c (\Pi (\frac{1}{2}d_{ma})^{2}-...\\ ...-\int _{-z(y_{co}=d_{ma}-a)}^{z(y_{co}=d_{ma}-a)}((d_{ma}-a)-y_{co})\mathrm{d}z) \end{array}</math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | (23)
|}

{| class="formulaSCP" style="width: 100%; text-align: left;" 
|-
| 
{| style="text-align: left; margin:auto;width: 100%;" 
|-
| style="text-align: center;" | <math>\begin{array}{l} F_{co}\cdot Y_{co}=f{''}c (-\int _{-z(y_{co}=d_{ma}-a)}^{z(y_{co}=d_{ma}-a)}((d_{ma}-...\\ ...-a)-y_{co})(-(\frac{1}{2}d_{ma}-...\\ ...-\frac{1}{2}(d_{ma}-a-y_{co}))\mathrm{d}z) \end{array}</math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | (24)
|}

==9 Geometría y propiedades mecánicas de los materiales==

Recordando que el objetivo es determinar un arreglo de refuerzo tal que cumpla con un rango de eficiencia aceptable ante ciertos elementos mecánicos actuantes al que la estructura concreto estaría sometidos, manteniendo constante la geometría de las sección propuesta inicialmente. Por lo que nuestros datos de entrada serán el número de elementos a diseñar, para cada elemento se indicará la resistencia a compresión simple del concreto <math>f{'}_{c}</math>, la resistencia a la fluencia del acero <math>f_{y}</math>, y el recubrimiento <math>rec</math>.

El programa también necesitará las fuerzas actuantes en la columna para su diseño. Dichas fuerzas se determinan normalmente mediante un análisis estructural, el cual es común que se realice en algún software para tal fin. A partir de éste punto se determinará un valor inicial del espesor del perfil idealizado de acero <math>t</math>, acomodado en la secicón de concreto de tal manera que se respete el recubrimiento deseado, para así comenzar con las iteraciones de actualización de los diagramas de interacción de las columnas hasta que la condición mecánica actuante más crítica sea cubierta.

==10 Cálculo de los diagramas de interacción==

Se establece un cierto número de puntos para determinar los diagramas de interacción, determinado los valores <math>P_{ot}</math>, <math>P_{oc}</math> como los límites del rango de las fuerzas del diagrama de interacción, pudiendo extraer una fuerza respectiva a cada punto del mismo que tendrá el papel de raíz en el método numérico a emplear para el cálculo de su respectivo momento resistente.

El método numérico empleado para la aproximación de ráices es el llamado ''Método de la Falsa posición. <span id='citeF-2'></span>[[#cite-2|[2]]]

===10.1 Cálculo de la eficiencia mecánica-estructural===

Para cada condición de fuerza actuante en la columna se determinará la respectiva eficiencia estructural que hay entre ésta condición y la resistencia que el diagrama de interacción índica.

La diferencia en el cálculo de éstas eficiencias en columnas circulares y rectangulares es que en las rectangulares se deberán extraer de los diagramas de interacción calculados en el sentido de <math>x</math> y <math>y</math> la carga resistente y el momento resistente (<math>M_{rx}, M_{ry}, P_{rx}, P_{ry}</math>  [p. actual]), y para las circulares al ser su sección simétrica respecto a cualquier eje, se calculará la eficiencia respecto a una sola dirección, extrayendo la carga resistente y el momento resistente (<math>M_{r}, P_{r}</math>) del diagrama de interacción calculado.

Para lo anterior se empleará geometría analítica con el sistema cartesiano de los propios diagramas de interacción Figura [[#img-13|13]] calculando el punto de intersección entre dos rectas imaginarias A y B, tal que A partirá del origen pasando por alguna condición de fuerza actuante prolongándose hasta el borde del diagrama de interacción, y B será tal que una el punto anterior <math>(M_{i}, P_{i})</math> y posterior <math>(M_{i+1}, P_{i+1})</math> del punto en el cual la recta A intersecta al diagrama de interacción.   <div id='img-13'></div>
{| class="floating_imageSCP" style="text-align: center; border: 1px solid #BBB; margin: 1em auto; width: 100%;max-width: 100%;"
|-
|[[Image:Review_975456380192-diagintreferencia.png|375px|Diagrama de interacción en el plano cartesiano de referencia para el programa.]]
|- style="text-align: center; font-size: 75%;"
| colspan="1" | '''Figura 13:''' Diagrama de interacción en el plano cartesiano de referencia para el programa.
|}
De modo que el cálculo del momento resistente para una respectiva condición se cálcula como:

<span id="eq-25"></span>
{| class="formulaSCP" style="width: 100%; text-align: left;" 
|-
| 
{| style="text-align: left; margin:auto;width: 100%;" 
|-
| style="text-align: center;" | <math>  M_{r}=\frac{P_{i+1}+(\frac{P_{i}-P_{i+1}}{M_{i+1}-M_{i}})}{\frac{P_{u}}{M_{u}}-\frac{P_{i+1}-P_{i}}{M_{i+1}-M_{i}}} </math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | (25)
|}

Y su fuerza resistente respectiva como:

<span id="eq-26"></span>
{| class="formulaSCP" style="width: 100%; text-align: left;" 
|-
| 
{| style="text-align: left; margin:auto;width: 100%;" 
|-
| style="text-align: center;" | <math>  P_{r}=\frac{P_{u}}{M_{u}}M_{r} </math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | (26)
|}

Ver demostración Anexo 1. p. actual

Ahora ya se pueden aplicar las fórmulas que índican las NTC-2017 p. actual, para calcular las eficiencias de cada condición mecánica actuante.

===10.2 Cálculo de área requerida de acero===

Ya determinada la eficiencia mecánica estructural para cada condición de carga actuante se determinará a continuación la más crítica, y a partir de ésta se cálculará el espesor del perfil de acero simplemente iterando el valor de <math>t</math> con un cierto incremento <math>dt</math>, hasta llegar a un rango de eficiencia crítica requerida, establecida a priori, pues está será la eficiencia que nos servirá de comparación para seguir iterando o no. Figura [ [[#img-14|14]]]

<div id='img-14'></div>
{| class="floating_imageSCP" style="text-align: center; border: 1px solid #BBB; margin: 1em auto; width: 100%;max-width: 100%;"
|-
|[[Image:Review_975456380192-condcritica.png|359px|Condición de carga crítica de diseño.]]
|- style="text-align: center; font-size: 75%;"
| colspan="1" | '''Figura 14:''' Condición de carga crítica de diseño.
|}

==11 Determinación del acero de refuerzo==

Una vez que se ha determinado el área de acero a emplear para que se alcance la eficiencia deseada, se debe transformar esa área ficticia de acero en varillas corrugadas, distribuyéndolas uniformemente en la sección (obedeciendo las hipótesis de diseño de la Sección [[#6.3.2 Flexo-compresión|6.3.2]]), por lo que sólo se admitirán números pares de varillas.

El programa genera un arreglo matricial de propuestas de número de varillas, desde 4 (para columnas rectangulares) o 6 (para columnas circulares), calculando para cada diámetro el número de varillas requeridas, tomando como referencia los tipos de varillas disponibles comercialmente.

El usuario podrá elegir no sólo el diámetro del acero de refuerzo, sino también la distribución de éstas en la sección (solo para las columnas rectangulares, ya que para circulares la distribución es simétrica respecto a cualquier eje de simetría), de manera que siempre se respete la separación mínima o distancia libre entre varillas <math>S_{min}=\frac{3}{2}TMA</math>. Aunque por defecto lo que usualmente se busca es la opción que resulte la más eficiente o la más económica. De manera que el usuario también podrá elegir que criterio es el más conveniente.

Cabe mencionar que la distribución del acero de refuerzo tiene influencia en la eficiencia mecánica de los elementos, y para cada opción de número de varillas/diámetro le corresponderá un buen número de configuraciones de éstas en las sección, por lo que el programa hará también un análisis de eficiencia de cada opción disponible, pero ahora tratando al acero de refuerzo precisamente como varillas de acero, para verificar cuáles cumplen con el rango requerido de eficiencia que se establecio a priori, una vez hecho esto, el programa seleccionará la opción más eficiente.    </span>
{|  class="floating_tableSCP wikitable" style="text-align: center; margin: 1em auto;min-width:50%;"
|+ style="font-size: 75%;" |<span id='table-1'></span>Tabla. 1 Tabla de tipos de varillas comerciales en México
|-
|   Tipo (octavo) 
| <math>diam(pulg)</math> 
| <math>diam(cm)</math> 
|  área<math>(cm^{2}</math>
|-
|   <span style="font-size: 75%;"><math>\# 3</math> </span>
|  <span style="font-size: 75%;"><math>0.375</math> </span>
|  <span style="font-size: 75%;"><math>0.9525</math> </span>
|  <span style="font-size: 75%;"><math>0.7126</math></span>
|-
|  <span style="font-size: 75%;"><math>\# 4</math> </span>
|  <span style="font-size: 75%;"><math>0.5</math> </span>
|  <span style="font-size: 75%;"><math>1.27</math> </span>
|  <span style="font-size: 75%;"><math>1.2668</math></span>
|-
|  <span style="font-size: 75%;"><math>\# 5</math> </span>
|  <span style="font-size: 75%;"><math>0.625</math> </span>
|  <span style="font-size: 75%;"><math>1.5875</math> </span>
|  <span style="font-size: 75%;"><math>1.9793</math></span>
|-
|  <span style="font-size: 75%;"><math>\# 6</math> </span>
|  <span style="font-size: 75%;"><math>0.75</math> </span>
|  <span style="font-size: 75%;"><math>1.905</math> </span>
|  <span style="font-size: 75%;"><math>2.8502</math></span>
|-
|  <span style="font-size: 75%;"><math>\# 8</math> </span>
|  <span style="font-size: 75%;"><math>1.0</math> </span>
|  <span style="font-size: 75%;"><math>2.54</math> </span>
|  <span style="font-size: 75%;"><math>5.067</math></span>
|-
|  <span style="font-size: 75%;"><math>\# 9</math> </span>
|  <span style="font-size: 75%;"><math>1.125</math> </span>
|  <span style="font-size: 75%;"><math>2.8575</math> </span>
|  <span style="font-size: 75%;"><math>6.413</math></span>
|-
|  <span style="font-size: 75%;"><math>\# 10</math> </span>
|  <span style="font-size: 75%;"><math>1.25</math> </span>
|  <span style="font-size: 75%;"><math>3.175</math> </span>
|  <span style="font-size: 75%;"><math>7.917</math></span>
|-
|  <span style="font-size: 75%;"><math>\# 12</math> </span>
|  <span style="font-size: 75%;"><math>1.5</math> </span>
|  <span style="font-size: 75%;"><math>3.81</math> </span>
|  <span style="font-size: 75%;"><math>11.400</math></span>

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

==12 Diseño de experimento==

Para cada tipo de geometría se probó un modelo estructural de experimentación con diferentes valores de <math>f'{c}</math>.

===12.1 Columnas rectangulares===

En un análisis estructural hecho previamente se proponen las siguientes geometrías y <math>f'{c}</math>:

====12.1.1 Datos de entrada====

</span>
{|  class="floating_tableSCP wikitable" style="text-align: center; margin: 1em auto;min-width:50%;"
|+ style="font-size: 75%;" |<span id='table-2'></span>Tabla. 2 Datos de entrada para el diseño de acero de refuerzo en columnas rectangulares
|-
| <math>b(cm)</math> 
| <math>h(cm)</math> 
| <math>rec(cm)</math> 
| <math>f'{c}(cm^{2}</math>
|-
|   <span style="font-size: 75%;"><math>40</math> </span>
|  <span style="font-size: 75%;"><math>75</math> </span>
|  <span style="font-size: 75%;"><math>5</math> </span>
|  <span style="font-size: 75%;"><math>280</math></span>

|}
<span style="text-align: center; font-size: 75%;">     Con las siguientes combinaciones de carga para ambas estructuras:

</span>
{|  class="floating_tableSCP wikitable" style="text-align: center; margin: 1em auto;min-width:50%;"
|+ style="font-size: 75%;" |<span id='table-3'></span>Tabla. 3 Combinaciones de carga resultantes de un análisis estructural previo para ambos modelos estructurales.
|-
| <math>P(Ton)</math> 
| <math>M_{x}(Ton-m)</math> 
| <math>M_{y}(Ton-m)</math>
|-
|   <span style="font-size: 75%;"><math>-35</math> </span>
|  <span style="font-size: 75%;"><math>22</math> </span>
|  <span style="font-size: 75%;"><math>45</math></span>
|-
|  <span style="font-size: 75%;"><math>36</math> </span>
|  <span style="font-size: 75%;"><math>25</math> </span>
|  <span style="font-size: 75%;"><math>32</math></span>
|-
|  <span style="font-size: 75%;"><math>-304</math> </span>
|  <span style="font-size: 75%;"><math>19</math> </span>
|  <span style="font-size: 75%;"><math>65</math></span>
|-
|  <span style="font-size: 75%;"><math>-46</math> </span>
|  <span style="font-size: 75%;"><math>12</math> </span>
|  <span style="font-size: 75%;"><math>75</math></span>
|-
|  <span style="font-size: 75%;"><math>-187</math> </span>
|  <span style="font-size: 75%;"><math>10</math> </span>
|  <span style="font-size: 75%;"><math>45</math></span>

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

====12.1.2 Resultados====

Una vez hecho el análisis de resistencia con las cargas y datos dados se obtienen los siguientes diagramas de interacción, tabla de eficiencias estructurales para cada modelo estructural con los grosores resultantes del perfil idealizado de acero, distribución óptima en costo y eficiencia de varillas de acero en cada modelo estructural y sus respectivas tablas de eficiencia con dicha opción de varillado:

Eficiencia estructural con espesores de perfiles idealizados de acero resultantes:

<div id='img-15'></div>
{| class="floating_imageSCP" style="text-align: center; border: 1px solid #BBB; margin: 1em auto; width: 100%;max-width: 100%;"
|-
|[[Image:Review_975456380192-ef_tabla_t_colrec01.png|375px|Tabla de eficiencias con espesor de perfil (t) resultante-Modelo estructural 01.]]
|- style="text-align: center; font-size: 75%;"
| colspan="1" | '''Figura 15:''' Tabla de eficiencias con espesor de perfil (t) resultante-Modelo estructural 01.
|}
Diagramas de interacción con espesor de perfil (t) resultante para cada modelo:

Nota:Los puntos rojos representan las condiciones de carga en ese sentido.  <div id='img-16'></div>
{| class="floating_imageSCP" style="text-align: center; border: 1px solid #BBB; margin: 1em auto; width: 100%;max-width: 100%;"
|-
|[[Image:Review_975456380192-diagIntPerfilX_rec_col01.png|525px|Diagrama de interacción en X con espesor de perfil (t) resultante-Modelo estructural 01.]]
|- style="text-align: center; font-size: 75%;"
| colspan="1" | '''Figura 16:''' Diagrama de interacción en X con espesor de perfil (t) resultante-Modelo estructural 01.
|}
<div id='img-17'></div>
{| class="floating_imageSCP" style="text-align: center; border: 1px solid #BBB; margin: 1em auto; width: 100%;max-width: 100%;"
|-
|[[Image:Review_975456380192-diagIntPerfilY_rec_col01.png|525px|Diagrama de interacción en Y con espesor de perfil (t) resultante-Modelo estructural 01.]]
|- style="text-align: center; font-size: 75%;"
| colspan="1" | '''Figura 17:''' Diagrama de interacción en Y con espesor de perfil (t) resultante-Modelo estructural 01.
|}
Opciones de distribución de varillas disponibles resultantes:

<div id='img-18'></div>
{| class="floating_imageSCP" style="text-align: center; border: 1px solid #BBB; margin: 1em auto; width: 100%;max-width: 100%;"
|-
|[[Image:Review_975456380192-opciones_var_rec_01.png|288px|Tabla de opciones disponibles de varillado-Modelo 01.]]
|- style="text-align: center; font-size: 75%;"
| colspan="1" | '''Figura 18:''' Tabla de opciones disponibles de varillado-Modelo 01.
|}
Distribución de varillado óptima en costo y eficiencia

<div id='img-19'></div>
{| class="floating_imageSCP" style="text-align: center; border: 1px solid #BBB; margin: 1em auto; width: 100%;max-width: 100%;"
|-
|[[Image:Review_975456380192-seccion_col01_rec.png|525px|Distribución de opción óptima de varillado en costo y eficiencia estructural correspondiente a la opción 3 Tabla [[#img-18|18]]-Modelo 01]]
|- style="text-align: center; font-size: 75%;"
| colspan="1" | '''Figura 19:''' Distribución de opción óptima de varillado en costo y eficiencia estructural correspondiente a la opción 3 Tabla [[#img-18|18]]-Modelo 01
|}
Tomando como referencia las distribuciones de varillas en la sección del modelo con las tablas de opciones disponibles de varillado, se puede notar que el programa al tomar la opción más económica, hace la selección tomando en cuenta también la mejor eficiencia estructural, que aunque no es la más eficiente de manera absoluta, si lo es con respecto al costo más bajo, ya que se pueden presentar mismos costos óptimos con diferentes eficiencias estrucrales.

Diagramas de interacción para las opciones de varillado óptimas en costo para cada modelo estructural:

<div id='img-20'></div>
{| class="floating_imageSCP" style="text-align: center; border: 1px solid #BBB; margin: 1em auto; width: 100%;max-width: 100%;"
|-
|[[Image:Review_975456380192-diagIntVarX_rec_col01.png|525px|Diagrama de interacción en X con varillas resultantes-Modelo estructural 01.]]
|- style="text-align: center; font-size: 75%;"
| colspan="1" | '''Figura 20:''' Diagrama de interacción en X con varillas resultantes-Modelo estructural 01.
|}
<div id='img-21'></div>
{| class="floating_imageSCP" style="text-align: center; border: 1px solid #BBB; margin: 1em auto; width: 100%;max-width: 100%;"
|-
|[[Image:Review_975456380192-diagIntVarY_rec_col01.png|525px|Diagrama de interacción en Y con varillas resultantes-Modelo estructural 01.]]
|- style="text-align: center; font-size: 75%;"
| colspan="1" | '''Figura 21:''' Diagrama de interacción en Y con varillas resultantes-Modelo estructural 01.
|}

===12.2 Columnas circulares===

====12.2.1 Datos de entrada====

</span>
{|  class="floating_tableSCP wikitable" style="text-align: center; margin: 1em auto;min-width:50%;"
|+ style="font-size: 75%;" |<span id='table-4'></span>Tabla. 4 Datos de entrada para el diseño de acero de refuerzo en columnas circulares
|-
| <math>D(cm)</math> 
| <math>rec(cm)</math> 
| <math>f'{c}(cm^{2}</math>
|-
|   <span style="font-size: 75%;"><math>40</math> </span>
|  <span style="font-size: 75%;"><math>5</math> </span>
|  <span style="font-size: 75%;"><math>300</math></span>

|}
<span style="text-align: center; font-size: 75%;">     Con las siguientes combinaciones de carga para ambas estructuras:

</span>
{|  class="floating_tableSCP wikitable" style="text-align: center; margin: 1em auto;min-width:50%;"
|+ style="font-size: 75%;" |<span id='table-5'></span>Tabla. 5 Combinaciones de carga resultantes de un análisis estructural previo para ambos modelos estructurales.
|-
| <math>P(Ton)</math> 
| <math>M_{x}(Ton-m)</math>
|-
|   <span style="font-size: 75%;"><math>-35</math> </span>
|  <span style="font-size: 75%;"><math>22</math></span>
|-
|  <span style="font-size: 75%;"><math>36</math> </span>
|  <span style="font-size: 75%;"><math>25</math></span>
|-
|  <span style="font-size: 75%;"><math>-304</math> </span>
|  <span style="font-size: 75%;"><math>19</math></span>
|-
|  <span style="font-size: 75%;"><math>-46</math> </span>
|  <span style="font-size: 75%;"><math>12</math></span>
|-
|  <span style="font-size: 75%;"><math>-187</math> </span>
|  <span style="font-size: 75%;"><math>10</math></span>

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

====12.2.2 Resultados====

Una vez hecho el análisis de resistencia con las cargas y datos dados se obtienen los siguientes diagramas de interacción, tabla de eficiencias estructurales para cada modelo estructural con los grosores resultantes del perfil idealizado de acero, distribución óptima en costo y eficiencia de varillas de acero en cada modelo estructural y sus respectivas tablas de eficiencia con dicha opción de varillado:

Eficiencia estructural con espesores de perfiles idealizados de acero resultantes:

<div id='img-22'></div>
{| class="floating_imageSCP" style="text-align: center; border: 1px solid #BBB; margin: 1em auto; width: 100%;max-width: 100%;"
|-
|[[Image:Review_975456380192-ef_cir_01_perfil.png|237px|Tabla de eficiencias con espesor de perfil (t) resultante-Modelo estructural 01.]]
|- style="text-align: center; font-size: 75%;"
| colspan="1" | '''Figura 22:''' Tabla de eficiencias con espesor de perfil (t) resultante-Modelo estructural 01.
|}
Diagramas de interacción con espesor de perfil (t) resultante para cada modelo:

<div id='img-23'></div>
{| class="floating_imageSCP" style="text-align: center; border: 1px solid #BBB; margin: 1em auto; width: 100%;max-width: 100%;"
|-
|[[Image:Review_975456380192-diagIntPerfil_cir01.png|525px|Diagrama de interacción con espesor de perfil (t) resultante-Modelo estructural 01.]]
|- style="text-align: center; font-size: 75%;"
| colspan="1" | '''Figura 23:''' Diagrama de interacción con espesor de perfil (t) resultante-Modelo estructural 01.
|}
Opciones de distribución de varillas disponibles resultantes:

<div id='img-24'></div>
{| class="floating_imageSCP" style="text-align: center; border: 1px solid #BBB; margin: 1em auto; width: 100%;max-width: 100%;"
|-
|[[Image:Review_975456380192-ef_opc_cir_01.png|196px|Tabla de opciones disponibles de varillado-Modelo 01.]]
|- style="text-align: center; font-size: 75%;"
| colspan="1" | '''Figura 24:''' Tabla de opciones disponibles de varillado-Modelo 01.
|}
Distribución de varillado óptima en costo y eficiencia

<div id='img-25'></div>
{| class="floating_imageSCP" style="text-align: center; border: 1px solid #BBB; margin: 1em auto; width: 100%;max-width: 100%;"
|-
|[[Image:Review_975456380192-seccion_cir01.png|525px|Distribución de opción óptima de varillado en costo y eficiencia estructural correspondiente a la opción 2 Tabla [[#img-24|24]]-Modelo 01]]
|- style="text-align: center; font-size: 75%;"
| colspan="1" | '''Figura 25:''' Distribución de opción óptima de varillado en costo y eficiencia estructural correspondiente a la opción 2 Tabla [[#img-24|24]]-Modelo 01
|}
Diagramas de interacción para las opciones de varillado óptimas en costo para cada modelo estructural:

<div id='img-26'></div>
{| class="floating_imageSCP" style="text-align: center; border: 1px solid #BBB; margin: 1em auto; width: 100%;max-width: 100%;"
|-
|[[Image:Review_975456380192-diagIntVar_cir01.png|525px|Diagrama de interacción con varillas resultantes-Modelo estructural 01.]]
|- style="text-align: center; font-size: 75%;"
| colspan="1" | '''Figura 26:''' Diagrama de interacción con varillas resultantes-Modelo estructural 01.
|}

==13 Discusión y conclusiones==

El método aparte de ser muy fácil de programar en una computadora ofrece una nueva forma para analizar y diseñar estructuras de concreto reforzado. Sin duda que aún queda mucho por desarrollar y este puede es el comienzo de nuevas investigaciones en material de optimización de estrucruras de concreto, ya que como tal este método también se implementó en le presente trabajo una optimización, aunque con un método muy sencillo de búsqueda, más sin embargo los resultados obtenidos son los óptimos de acuerdo a los requerimientos de diseño.

==14 Agradecimientos==

Agradezco a la Universidad de Guanajuato por esta oportunidad de trabajo, así como especial agradecimiento al Aula CIMNE-UG dónde se desarrolló el programa computacional como tal en lenguaje C principalmente, con MatLab para la parte gráfica, y Dynamo-Revit para la programación visual paramétrica de los elementos estructurales y su acero de refuerzo. A mi director de tesis el Dr. Alejandro Hernández Martínez y a mi tutor proyectos en el Aula CIMNE el M.C. Humberto Esqueda Oliva.

==15 Anexo 1. Demostraciones del desarrollo del programa==

===15.1 Geometría analítica para el cálculo de eficiencias mecánicas===

Ecuación de la recta A:

{| class="formulaSCP" style="width: 100%; text-align: left;" 
|-
| 
{| style="text-align: left; margin:auto;width: 100%;" 
|-
| style="text-align: center;" | <math> y=\frac{y_{u}}{x_{u}}x </math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | (27)
|}

Ecuación de la recta B:

{| class="formulaSCP" style="width: 100%; text-align: left;" 
|-
| 
{| style="text-align: left; margin:auto;width: 100%;" 
|-
| style="text-align: center;" | <math> y=\frac{y_{i+1}-y_{i}}{x_{i+1}-x_{i}}x+(y_{i+1}+(\frac{y_{i}-y_{i+1}}{x_{i+1}-x_{i}})x_{i+1}) </math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | (28)
|}

Igualando las ecuaciones anteriores por compatibilidad:

<span id="eq-29"></span>
{| class="formulaSCP" style="width: 100%; text-align: left;" 
|-
| 
{| style="text-align: left; margin:auto;width: 100%;" 
|-
| style="text-align: center;" | <math>  \frac{y_{u}}{x_{u}}x=\frac{y_{i+1}-y_{i}}{x_{i+1}-x_{i}}x+(y_{i+1}+(\frac{y_{i}-y_{i+1}}{x_{i+1}-x_{i}})x_{i+1}) </math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | (29)
|}

Despejando x de la ecuación anterior:

<span id="eq-30"></span>
{| class="formulaSCP" style="width: 100%; text-align: left;" 
|-
| 
{| style="text-align: left; margin:auto;width: 100%;" 
|-
| style="text-align: center;" | <math>  x=x_{r}=\frac{y_{i+1}+(\frac{y_{i}-y_{i+1}}{x_{i+1}-x_{i}})}{\frac{y_{u}}{x_{u}}-\frac{y_{i+1}-y_{i}}{x_{i+1}-x_{i}}} </math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | (30)
|}

Sustituyendo la Ecuación [ [[#eq-25|25]]] en Ecuación [ [[#eq-29|29]]].

<span id="eq-31"></span>
{| class="formulaSCP" style="width: 100%; text-align: left;" 
|-
| 
{| style="text-align: left; margin:auto;width: 100%;" 
|-
| style="text-align: center;" | <math>  y=y_{r}=\frac{y_{u}}{x_{u}}x_{r} </math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | (31)
|}

Donde: <math>y_{u}=P_{u}</math>

<math>x_{u}=M_{u}</math>

<math>y_{r}=P_{r}</math>

<math>x_{r}=M_{r}</math>

<math>y_{i}=P_{i}</math>

<math>x_{i}=M_{i}</math>

<math>y_{i+1}=P_{i+1}</math>

<math>x_{i+1}=M_{i+1}</math>

==16 Referencias bibliográficas==

===BIBLIOGRAFÍA===

<div id="cite-1"></div>
'''[1]'''  Adolfo Iván Jiménez P. y Dra. Sonia Elda Ruiz Gómez, Factores óptimos de cargas para el diseño de columnas esbeltas, Tesis de Licenciatura, UNAM, México, 2013   <div id="cite-2"></div>
'''[[#citeF-2|[2]]]'''  Steven C. Chapra, Raymond P. Canale, Métodos numéricos para ingenieros, Quinta edición, McGraw Hill  <div id="cite-3"></div>
'''[[#citeF-3|[3]]]'''  Gaceta Oficial de la Ciudad de México, Normas técnicas complementarias de la Ciudad de México, 2017, Normas Técnicas Complementarias para Diseño y Construcción de Estructuras de Concreto 2017  <div id="cite-4"></div>
'''[[#citeF-4|[4]]]'''  González Cuevas, Fco. Robles Hernández, Aspectos básicos del concreto reforzado, 4ta edición, Limusa 2005  <div id="cite-5"></div>
'''[5]''' Jack McCormac y Rusell H. Brows, Diseño de concreto reforzado, 14va edición, Alfaomega, 2015   <div id="cite-6"></div>
'''[[#citeF-6|[6]]]'''  Turner, M. J., Clough, R. W., Martin H. C. and Topp, L. J. Stiffness and Deflection Analysis of Complex Structures . Journal of the Aeronautical Sciences, Vol. 23 No. 9, 1956 pp. 805-823.  <div id="cite-7"></div>
'''[7]'''  Yoeber Castro Atau, “Sistematización de detalles, habilitación y armado de acero ASTM A615 para construcciones de concreto armado: impacto técnico, económico y ambiental”, Universidad Nacional de San Cristobal de Huamanga, 2010  

</span>