UNIVERSIDAD NACIONAL DE TRUJILLO

FACULTAD DE INGENIERÍA

ESCUELA PROFESIONAL DE INGENIERÍA INDUSTRIAL

Optimización del valor de acero requerido para el concreto armado utilizando

el método de Newton Raphson y el método de la búsqueda de la sección

dorada.

Informe Final Grupo 2

MONOGRAFÍA

ALVAREZ QUIROZ, MARJORY BRISSET

CÁRDENAS ALCÁNTARA, LUIS ENRIQUE
LEÓN SANTILLAN, NATHALY NOHELIA
PAJARES BERNABÉ, LUIS CARLOS
QUISPE CASTILLO, LIZANDRO
VARGAS ZAPATA, ALEXIS (COORDINADOR)
MÉTODOS NUMÉRICOS

CICLO NIVELACIÓN

GELDRES MARCHENA TEODORO ALBERTO

TRUJILLO – 07/04/2022
 1

ÍNDICE
A. INTRODUCCIÓN 3

B. DESARROLLO TEMÁTICO 4

1. CONCRETO REFORZADO 4

2. VIGAS DE SECCIÓN RECTANGULAR DE CONCRETO REFORZADO 5

2.1. VIGAS 6

3. VIGAS RECTANGULARES SIMPLEMENTE ARMADAS 6

4. MÉTODO DE NEWTON-RAPHSON 9

4.1. ANÁLISIS DE CONVERGENCIA DEL MÉTODO DE NEWTON

RAPHSON 10

5. SITUACIÓN PROBLEMÁTICA 10

6. SOLUCIÓN DEL PROBLEMA 11

C. DISCUSIÓN 12

D. CONCLUSIONES Y RECOMENDACIONES 12

E. REFERENCIAS BIBLIOGRÁFICAS 13

APÉNDICE 16

ANEXO 1 17

ANEXO 2 18

ANEXO 3 20

ANEXO 4 21
 2

ÍNDICE DE TABLAS

Tabla de áreas del acero 21

ÍNDICE DE FIGURAS

Ilustración 1 Elementos estructurales de concreto armado (Colección Digital

UANL. Cap.1, p. 4) 5

Ilustración 2 Elementos estructurales de concreto armado (Colección Digital

UANL. Cap.1, p. 4)6

Ilustración 3 Dimensiones de una viga armada 10

Ilustración 4 Diagrama de flujo 17

Ilustración 5 Corrida del programa 20

 3

A. INTRODUCCIÓN

Este informe presenta el valor óptimo de acero mínimo en vigas de concreto

armado simple o doblemente armadas sometidas a flexión. Se utilizó la norma
ACI 318-04 para el diseño de la viga de concreto armado, ya que es un
elemento estructural importante en una edificación, por este motivo es
necesario usar normas y criterios al momento de diseñar la sección de viga.

El acero ha jugado un papel protagónico en el desarrollo de nuestra sociedad,

y un papel protagónico con una trayectoria de muchos años no ha perdido un
ápice de importancia con el paso del tiempo, sino todo lo contrario. Estamos
hablando de una industria dinámica que tiene mucho alcance global, es
fundamental en áreas como la construcción en la fabricación en general. En la
industria manufactura y de construcción, el acero juega un papel muy
importante siendo un material crítico porque ningún otro reúne sus
características que son: resistencia, plasticidad y versatilidad. Según
Worldsteel más de la mitad de acero producido en el mundo se destina a
edificios e infraestructuras.

Según L. R. Prieto-Serrano en su trabajo titulado “Método de la deformación

unitaria para diseño de vigas rectangulares doblemente reforzadas” Se
determina que, en el diseño de vigas armadas en compresión, la deformación
unitaria óptima de la armadura extrema en tracción bajo el estado último de
tracción es 0.005. En casos especiales, como aquellos en los que se permite
la redistribución de momentos, es necesario asegurarse de que. Con un
mínimo de 0,0075, la demanda de acero aumenta significativamente, casi un
30% más que la requerida para una deformación unitaria de 0,005.

En el artículo de Chacón R. (2009) en su trabajo “Resistencia de vigas

armadas de acero sometidas a cargas concentradas en secciones no
rigidizadas” se ha observado que, según diferentes normativas a nivel
mundial, las comprobaciones de la resistencia de vigas metálicas sometidas a
 4

cargas concentradas difieren ligeramente en los modelos teóricos en los que

se basan y en los resultados finales de las predicciones que produce cada uno
de ellos. Para determinar la carga última de fatiga para vigas sometidas a
cargas concentradas, diferentes propuestas normativas requieren la
verificación de dos casos conceptualmente distintos.

Además de esto, se utilizará el método de Newton-Raphson y la sección

dorada para obtener el área de acero requerida para la tensión, que es una
parte esencial del diseñ[Link] comparar cual de los dos métodos es
el que alcanza un menor número de iteraciones y por lo consiguiente es el
más eficaz para hallar nuestra incógnita. De esta manera reconocer la
importancia y aplicación de los métodos numéricos en la ingeniería industrial.

B. DESARROLLO TEMÁTICO

1. CONCRETO REFORZADO

Concreto al cual se le ha adicionado un refuerzo de acero en barras,

mallas electrosoldadas, pernos con cabeza y fibras de acero
deformadas dispersas para absorber los esfuerzos que el concreto por
su propia condición no lo puede hacer, pero entendiéndose que el
trabajo de los dos materiales es de conjunto, es decir, a partir de la
compatibilidad de deformaciones de los dos materiales. (Segura
Franco, 2011)
El concreto por sus características es el material ideal para elementos
que se encuentran expuestos a esfuerzos de compresión. Sin embargo,
Arqhys Arquitectura (2012) menciona que, en eventos externos como
sismos, pueden generar esfuerzos multidireccionales que producen
tanto tensión como compresión; debido a la capacidad a resistir tensión
del concreto cercana al 10% de su capacidad a compresión, resulta
necesario reforzarlo con un soporte o elemento que pueda compensar
esta carencia, y ese soporte lo brinda el acero estructural que posee
una resistencia a la tensión de más de 100 veces la del concreto.
 5

Un elemento de concreto reforzado debe tener una cantidad

balanceada de concreto y acero, ya que si existen elementos con
exceso de acero estos fallarán a tensión, y en este puede presentarse
un aplastamiento del concreto antes de que el acero llegue a fluir. Por lo
contrario, si un elemento no tiene suficiente acero, éste colapsará ante
la presencia de la primera grieta.

Las estructuras de concreto reforzado se usan en construcciones

modernas y económicas, así mismo, consta de miembros individuales
que se apoyan transmitiendo cargas a la estructura

En la Figura 1 observamos un sistema de elementos estructurales de

concreto reforzados como: losa, columnas, vigas, zapatas.

2. VIGAS DE SECCIÓN RECTANGULAR DE CONCRETO REFORZADO

Una viga de hormigón es rectangular cuando, en su sección a

compresión, una de las dos dimensiones es al menos el doble de la
dimensión menor analizada y su longitud supera a las otras dos
 6

dimensiones.

2.1. VIGAS

Según Rodes J., menciona la siguiente teoría sobre este tema:

Una viga es un miembro estructural cuyas cargas aplicadas son
principalmente perpendiculares al eje, por lo que el diseño es
principalmente para flexión y cortante; si la carga no es
perpendicular, se genera una fuerza axial, pero esta no es crítica
en el diseño porque es débil . La deformación se ve afectada por
una variedad de factores, que incluyen:

● Longitud libre del elemento entre apoyos

● Las condiciones de continuidad en los extremos de los
elementos
● La forma de la sección resistente
● El material a emplear
● La magnitud de la carga

Según Lefevre R. y Villar D., (2015) es su libro “Resistencia de

Materiales” nos plantean esta definición:

La deformación por flexión siempre resulta en la tendencia de

una sección longitudinal a doblarse bajo una carga aplicada.
Esta curvatura, que es esencialmente curvatura, varía según las
condiciones de cada situación particular.

Las vigas de hormigón armado son las más utilizadas en la

construcción porque son más baratas y el método de
 7

construcción es más común, aunque el proceso de construcción

es más lento y complicado.

VIGAS RECTANGULARES SIMPLEMENTE ARMADAS

La característica principal de una viga armada simple es que solo tiene

armadura en la zona portante, por lo que sabemos que la armadura
está a tracción y el hormigón a compresión. (Universidad Autónoma De
Coahuila, sf.)

La viga puede presentar fallas, y se pueden dar de dos maneras:

● La primera se puede dar cuando el acero de refuerzo llega a su

límite elástico aparente o límite de fluencia
● El segundo tipo de falla sucede cuando el concreto alcanza su
límite de ruptura sin que el acero haya logrado superar su fatiga.

Para evitar a las vigas reforzadas y a las vigas balanceadas, el

reglamento ACI 318 – 04 limita el porcentaje de refuerzo al 75%
del valor correspondiente a las secciones balanceadas. Por otra
parte, también las vigas con porcentajes muy pequeños, suelen
fallar súbitamente; para evitar ese riesgo el reglamento ACI
318-04 exige que el porcentaje mínimo en miembros sujetos a
flexión sea de: (Cedeño, 2017)
14.5
ρ= 𝐹𝑦

Por equilibrio de fuerzas

'
𝐶 = 𝑇𝐶 = 0. 85𝐹 𝑐β1𝑏𝐶𝑇 = 𝐴𝑠. 𝐹𝑦

Por lo tanto:

' 𝐴𝑠 ' 𝑐
0. 85𝐹 𝑐β1𝑏𝐶 = 𝐴𝑠. 𝐹𝑦 𝑏𝑑
𝐹𝑦 = 0. 85𝐹 𝑐β1 𝑑
 8

Del diagrama de deformaciones, aceptando las condiciones de viga

balanceada:

𝐹𝑦
ε𝑐 = 0. 003 ε𝑦 = 𝐸𝑠

𝑐 ε𝑐 0.003 6115
𝑑
= ε𝑐+ε𝑦
= 𝐹𝑦 = 6115+𝐹𝑦
⎡0.003+ ⎤
⎢ 6 ⎥
⎣ 2.039×10 ⎦

Por lo tanto

'
0.85𝐹 𝑐β1 6115
ρ𝑏 = 𝐹𝑦
× 6115+𝐹𝑦

La expresión anterior representa el valor del porcentaje de refuerzo

en la sección balanceada de una viga. El ACI 318-08 limita el
porcentaje máximo aplicable a miembros sujetos a flexión a 75.

'
0.85𝐹 𝑐β1 6115
ρ𝑏 = 𝐹𝑦
× 6115+𝐹𝑦

El asignar a fs el valor de Fy. Se está considerando que el acero

fluye y la viga es sobrereforzada:
𝐴𝑠.𝐹𝑦
𝑐= '
0.85𝐹 𝑐β1𝑏𝐶

Si llamamos:

𝐴𝑠 ρ𝐹𝑦
ρ= 𝑏𝑑
𝑐= ' ×𝑑
0.85𝐹 𝑐β1

Por lo tanto:

2
𝑀𝑢 = 𝑏𝑑 𝐹'𝑐ω. (1 − 0. 59ω)

Estableciendo el momento último en función del acero de refuerzo

se produce de la siguiente manera, refiriéndose a la figura 3 y
empleando la cuña rectangular de Whitney:
 9

(
𝑀𝑢 = 𝐴𝑠 𝐹𝑦 𝑑 −
𝑎
2 )
Ambas expresiones del momento último el reglamento las propone
afectadas de un coeficiente de seguridad, para las vigas vale 0,9,
por lo que quedarían finalmente:

2 '
( ⎣
)
𝑀𝑢 = ∅ 𝑏. 𝑑 𝐹 𝑐ω. (1 − 0. 59ω) 𝑀𝑢 = ∅⎡⎢𝐴𝑠 𝐹𝑦(𝑑 −
𝑎
2
⎤⎥
⎦

𝐴𝑠 𝐹𝑦
Donde:𝑎 = '
0.85𝐹 𝑐𝑏

Y Reemplazando la ecuación de “a” en “Mu” obtenemos:

2 2
∅𝑓*𝑓𝑦*𝐴𝑠
1.7*𝑓𝑐*𝑏
− ∅𝑓 * 𝑓𝑦 * 𝑑 * 𝐴𝑠 + 𝑀𝑢 = 0

Y de aquí mediante métodos numéricos calculamos el As requerido.

(Cedeño, 2017)

3. MÉTODO DE NEWTON-RAPHSON

Según Costa Vigil (2011) el método de Newton (conocido también como

el método de Newton-Raphson o el método de Newton-Fourier) es un
algoritmo eficiente para encontrar aproximaciones de los ceros o raíces
de una función real. También puede ser usado para encontrar el
máximo o mínimo de una función, encontrando los ceros de su primera
derivada.

Según Vargas (2018) menciona en su artículo titulado “Método de

Newton Raphson”:

Sabemos que por el teorema de Taylor que para un Xo ∈ (𝑎, 𝑏) tal que
F´(Xo) es diferente de cero y además si se cumple que tanto la función
como la derivada son continuas en el intervalo (a,b), entonces:

( ) (
𝑓(𝑥) = 𝑓 𝑥0 + 𝑥 − 𝑥0 𝑓 𝑥0 + … ) '( )
 10

Buscamos el punto donde f(x1) = 0, entonces:

( ) ( ) '( )
0 = 𝑓 𝑥0 + 𝑥1 − 𝑥0 𝑓 𝑥 0 + …

Por lo tanto, podemos concluir que:

𝑓(𝑥0)
𝑥1 ≈ 𝑥0 − 𝑓'(𝑥0)

De forma general:

𝑓(𝑥𝑛)
𝑥𝑛+1 = 𝑥𝑛 − 𝑓'(𝑥𝑛)

[Link]ÁLISIS DE CONVERGENCIA DEL MÉTODO DE NEWTON

RAPHSON

El método de Newton Raphson es un método abierto,es decir, no

hay un criterio general de convergencia La única manera de
alcanzar la convergencia es seleccionar un valor inicial lo
suficientemente cercano a la raíz buscada.

4. MÉTODO DE LA SECCIÓN DORADA:

Según Huiman (2008),menciona que la búsqueda de la sección dorada

es una técnica simple de búsqueda de una sola variable de propósito
general. Es similar en esencia al enfoque de la bisección para localizar
raíces.

La clave para hacer eficiente este procedimiento es la mejor elección de

los puntos intermedios. Esta meta se puede alcanzar al especificar que
las siguientes dos condiciones se cumplan.
 11

Sustituyendo la primera ecuación en la segunda tenemos:

Si tomamos el recíproco y R= l1/l2 , se llega a:

La cual se puede resolver para la raíz positiva

De forma general:

5. SITUACIÓN PROBLEMÁTICA
En una construcción es necesario una viga simplemente armada (solo
acero en tracción) de 30 cm x 60 cm que soporte un momento último de 15
tn-m, para esta estructura se usará un concreto con f’c = 280 kg/cm2 y
 12

acero con fy = 4200 kg/cm2. Calcularemos el acero requerido (tracción)

para el diseño.

6. SOLUCIÓN DEL PROBLEMA

Haremos uso de la siguiente ecuación:

2 2
∅𝑓*𝑓𝑦*𝐴𝑠
1.7*𝑓𝑐*𝑏
− ∅𝑓 * 𝑓𝑦 * 𝑑 * 𝐴𝑠 + 𝑀𝑢 = 0

Vamos a asumir que d=53 cm, siempre y cuando consideramos el

recubrimiento y posible diámetro de la varilla de acero. Además, el
momento Mu debe ser convertido a kg-cm ya que f’c y fy están dados en
kg/cm2, por ese motivo multiplicamos a Mu por 105:

Se obtiene:

2 2
0.9*4200 *𝑥 5
1.7*280*30
− 0. 9 * 4200 * 53 * 𝑥 + 15 * 10 = 0

Y calculando obtenemos:

2
1111. 76470588 * 𝑥 − 200340 * 𝑥 + 1500000 = 0 … (*)
 13

La ecuación (*) es del modelo de una cuadrática, por lo tanto tendrá 2

raíces y comenzamos iterando desde el X0=0 ya que el As debe ser un
valor mínimo.

C. DISCUSIÓN

❖ Para la resolución de esta situación problemática se aplicó usando

el método de Newton-Raphson y el de la sección dorada,
obteniendo los siguientes resultados:
● Aplicando el método de Newton-Raphson obtuvimos la respuesta a
nuestra incógnita en la cuarta iteración con un error menor al
0.0007%, siendo de 7.8273 cm2.A partir de ello, podemos decir que
podría usarse un acero Nº 10 porque cuenta con un área de 7.917
cm2. (Anexo 4)
● Aplicando el método de la sección dorada obtuvimos la respuesta a
nuestra incógnita en la doceava iteración con un error menor al
0.0002%, siendo de 7.7002 cm2.

❖ Analizando ambos resultados expuestos anteriormente podemos

deducir que el método de newton Raphson es eficiente, por el
motivo de obtener la respuesta en un menor número de iteraciones,
en comparación con el método de la búsqueda de la sección dorada.
De esta manera corroboramos lo que nos menciona Altamirano. Y.
Ochoa. D. (2015). acerca de dicho método: “El método de newton es
eficiente en la solución de sistemas de ecuaciones no lineales,
converge muy rápidamente y proporciona una muy buena precisión
en los resultados. El método se emplea en la solución de problemas
académicos y en problemas del mundo real.”
 14

Sin embargo, existen algunas desventajas al emplear este método,

al no tener un criterio general de convergencia como lo tienen otros
métodos, puede no converger dependiendo de la naturaleza de la
función y de la complejidad de la derivada.

D. CONCLUSIONES Y RECOMENDACIONES

● El acero en tracción requerido para el diseño de la viga sometida a

flexión debe tener un área mínima de 7.8273cm2, es decir, se necesitará
un acero #10 para que cubra con este requerimiento.

● El método de Newton-Raphson es más eficiente a comparación del

método de la búsqueda de la sección dorada al hallar solución a
ecuaciones no lineales.

● Los métodos numéricos en la ingeniería industrial juegan un papel muy

importante ya que su objetivo principal es encontrar soluciones
aproximadas a problemas complejos haciendo uso de operaciones
simples.

● Los métodos numéricos en la ingeniería industrial,son una importante

herramienta en varios ámbito,uno de ellos es la optimizació[Link] esta
manera nos permite llegar a realizar correctamente balances de tiempo,
costo y calidad; además actualmente la tecnología juega un papel
decisivo al momento de crear algoritmos que permitan acercarnos a
resultados reales.

● En un enfoque práctico y económico,es importante saber el óptimo valor

requerido que necesitamos usar en estructuras de acero, esto nos
permite tener una mayor flexibilidad en el uso del edificio,permitiendo
construir más rápido y nos permite tener un mejor diseño interior sin
tener que incrementar la inversión. De esta manera los intereses se
reducen durante la construcción,se puede hacer uso de la
 15

infraestructura en un menor tiempo y los ingresos productos llegan

antes en beneficio del flujo de capital.

● Tener tesis e informes de referencia para poder comparar los

resultados obtenidos.

E. REFERENCIAS BIBLIOGRÁFICAS

American Concrete Institute. Requisitos de reglamento para concreto

estructural y comentario (ACI 318-08). Estados Unidos: ACI, 2008. 518

p. ISBN: 978-0-87031-264-9 Católica de Colombia. Bogotá, D.C

[Link]

%20DE%20GRADO%20Christian%20Camilo%20Torres%[Link]

Asmat, C. (2020, 11 februari). INFLUENCIA DE REFORZAMIENTO CON

CFRP EN LA DUCTILIDAD DE VIGAS DE CONCRETO ARMADO.

AGER - Ingeniería Estructural. Geraadpleegd op 8 april 2022, van

[Link]

ilidad-de-vigas-de-concreto-armado/

Altamirano. Y. & Ochoa. D. (2015). Ventajas y Desventajas de Métodos

de Bisección, Secante y Newton Raphson. UTT Recuperado de:

[Link]

ajas-de-mtodos-de-biseccin-secante-y-newton-raphson

L. R. Prieto-Serrano, “Método de la deformación unitaria para diseño de

vigas rectangulares doblemente reforzadas”, Fac. Ing., vol. 24 (40), pp.

31–40, Sep.-Dic. 2015.

 16

[Link]

Cedeño, R. (2017). Análisis De Viga Simplemente Armada.

[Link]

e-Armada

Costa.V (2011). APPLICATION OF NUMERICAL METHODS TO

SOLVE NONLINEAR. Universidad Ricardo Palma.

[Link]

to-solve-

Chacón, Rolando, Mirambell, Enrique, & Real, Esther. (2009).

Resistencia de vigas armadas de acero sometidas a cargas

concentradas en secciones no rígidas. Boletín Técnico , 47(1), 1-16.

Recuperado en 08 de abril de 2022, de

[Link]

000100001&lng=es&tlng=es.

Equipo de Redacción Profesional. (2012). Concreto Armado. Portal de

Arquitectura. [Link]

Garcia, A. (2014). Modelación de capacidad esperada en vigas de

sección rectangular hasta tres metros de longitud de concreto reforzado

al ser ensayadas en laboratorio con ayuda de métodos finitos.

Universidad de San Carlos de Guatemala.

[Link]
 17

a/[Link]

Luévanos. A, López. S y Medina. M (2018). “Optimización de vigas de

concreto reforzado para secciones rectangulares con experimentos

numéricos”. Computación y Sistemas, 22(2), pp. 599–606.

[Link]

Mendieta, L. (2017). Aplicación de concreto reforzado con fibras de

acero en losas de contrapiso para viviendas de interés social. [Tesis de

pregrado, Universidad Militar Nueva Granada].

[Link]

[Link]?sequence=

Rosales L. y Ordoñez J. (2008). “Análisis del comportamiento de una

viga de concreto armado, ensayada a flexión, reforzada con barras de

fibra de carbono”. Facultad de Ingeniería, ingenieria Civil, Universidad

San carlos de Guatemala.

[Link]

Universidad Autónoma De Coahuila. (s.f.). Diseño de Vigas

Rectagulares. Coahuila

Vargas. J (2018). Método de Newton Raphson.

[Link]

[Link]

Viga Rectangular Simplemente Armada. UADY. Facultad de Ingeniería.

[Link]
 18

Vigas de concreto armado: consejos para hacer un buen trabajo. (2020,

3 september). Construyendo Seguro. Geraadpleegd op 8 april 2022,

van

[Link]

os-para-hacer-un-buen-trabajo/

APÉNDICE

1. DIAGRAMA DE FLUJO
1.1. G2-DIAGRAMA DE FLUJO- MÉTODO DE NEWTON RAPHSON
(pdf)
1.2. G2-DIAGRAMA DE FLUJO-MÉTODO DE LA BÚSQUEDA DE LA
SECCIÓN DORADA (pdf)

2. CODIFICACIÓN
2.1. NEWTON-RAPHSON
3. MONOGRAFÍA
3.1. InformeFinalG2 (Word)
3.2. InformeFinalG2 (pdf)
4. CAPTURAS DE REUNIONES GRUPALES
4.1. IMAGEN 1
4.2. IMAGEN 2
5. PRUEBA DE ESCRITORIO

Link de documentación complementaria:

[Link]

Kq6w3Eo56
19
 20

ANEXO 1

Ilustración 4 Diagrama de flujo del método de Newton Raphson

‌
 21

ANEXO 2
Codificación del Método de Newton – Raphson

% Definimos "x" como simbolico

syms x
% TITULO DEL METODO DEL PROYECTO
disp('===============================================================')
disp(' METODO: NEWTON RAPHSON')
disp('===============================================================')
% Menu de opciones para elegir metodo, punto a buscar e ingresar datos.
disp(' Menu de opciones')
disp('---------------------------------------------------------------')
disp(' 1) Punto Minimo')
disp(' 2) Punto Maximo')
disp(' ')
opc2 = input('Elija la opcion del punto que desea encontrar: ');
disp('---------------------------------------------------------------')
disp(' Informacion general')
disp('---------------------------------------------------------------')
f = input(' Introduzca la funcion f(x): ');
err = input(' Introduzca el porcentaje de error (%): ');
Maxiter = input(' Introduzca el numero maximo de iteraciones: ');
x0 = input(' Ingrese el valor inicial x0: ');
disp(' ')
disp('===============================================================')
% Obtenemos las derivadas de la funcion
df = diff(f);

%==========================================================================
% DESARROLLO DEL METODO =======================================
%==========================================================================
% Utilizamos la funcion "inline" para poder reemplazar valores en las
% funciones
f = inline(f);
df = inline(df);

% Definimos el error actual como 100%

ea = 100;
% Inicializacion del numero de iteraciones
j=0;
% Escribir informacion inicial de los valores
%fprintf(' i xi f(xi) df(xi) error \n')
%fprintf('%2d %10.5f %10.5f %10.5f\n',j,x0,f(x0),df(x0));
% Bucle de ejecucion
while (ea > err) && (j < Maxiter)
% Aproximamos la raiz con la formula correspondiente
x = x0 - (f(x0)/df(x0));
% Calculamos el % de error
ea = (abs(x-x0)/x)*100;
% Escribir informacion de los valores
%fprintf('%2d %10.5f %10.5f %10.5f %10.5f\n',j+1,x,f(x),df(x),ea)
% Almacenamos el nuevo valor de x para la siguiente iteracion
x0=x;
% Aumentamos en 1 el numero de iteraciones
j=j+1;
end
disp(' RESULTADOS')
 22

disp('---------------------------------------------------------------')
disp(' ')
if df(x0) > 0
if opc2 == 1
disp(['Valor Minimo : ',num2str(x0)])
else
disp('No se encontro un valor maximo dentro del rango inicial;')
disp(['sin embargo, el valor minimo hallado fue: ',num2str(x0)])
end
else
if opc2 == 1
disp('No se encontro un valor minimo dentro del rango inicial;')
disp(['sin embargo, el valor maximo hallado fue: ',num2str(x0)])
else
disp(['Valor Maximo : ',num2str(x0)])
end
end
disp(['Punto para el óptimo : ',num2str(f(x0))])
disp(['Numero de iteraciones: ',num2str(j)])
disp(['Error porcentual : ',num2str(ea),' %'])
 23

ANEXO 3

Ilustración 5 Corrida del programa

ANEXO 4
 24

Codificación del Método de Sección Dorada

opcion = menu('ELEGIR EL MÉTODO','seccion dorada');

if opcion==1

disp('Metodo de sección dorada')

clear

syms x

f(x)=input('Ingrese la función de la cual desea hallar su máximo

f(x)=');

xl=input('Ingrese el extremo izquierdo del intervalo xl:');

xu=input('Ingrese el extremo derecho del intervalo xu:');

MaxIter=input('Ingrese el número maximo de iteraciones Numit:');

tol=input('Ingrese el criterio de pare tol:');

d=((sqrt(5)-1)/2)*(xu-xl);

x1=xl+d;

x2=xu-d;

fprintf( '\n\n i xl f(xl) x2 f(x2)

x1 f(x1) xu f(xu) d\n')

fprintf('%3d
%12.6f%12.6f%12.6f%12.6f%12.6f%12.6f%12.6f%12.6f%12.6f\n',...

1,xl,f(xl),x2,f(x2), x1,f(x1), xu, f(xu),d)

k=1;

while k<=MaxIter

if f(x2)>f(x1)

xl=xl;

xu=x1;

else
 25

xl=x2;

xu=xu;

end

d=((sqrt(5)-1)/2)*(xu-xl);

x1=xl+d;

x2=xu-d;

fprintf('%3d
%12.6f%12.6f%12.6f%12.6f%12.6f%12.6f%12.6f%12.6f%12.6f\n',...

k+1,xl,f(xl),x2,f(x2), x1,f(x1), xu,

f(xu),d);

if abs(x2-x1)<tol

fprintf(' El maximo se encuentra en el punto

x=%8.4f \n',x2)

fprintf(' El maximo es f(x)= %8.4f \n',f(x2))

fprintf('El numero maximo de iteraciones necesarias

es %4d\n',k)

disp('Procedimiento terminado satisfactoriamente');

break;

else

k=k+1;

if k>MaxIter

fprintf('El metodo fracaso³ despues de %d

iteraciones \n',MaxIter)

disp('Procedimiento terminado sin éxito');

end

end
 26

end

end

ANEXO 5
 27

ANEXO 6

Tabla 2. Áreas del acero

Barra Diámetro Peso Área Perímetro

Num Pulg Cm Kg/m 𝐶𝑚

2 Cm

#2 1/4 0.635 0.250 0.317 2

#3 3/8 0.953 0.559 0.713 3

#4 1/2 1.270 0.993 1.267 4

#5 5/8 1.588 1.552 1.981 5

#6 3/4 1.905 2.235 2.850 6

#7 7/8 2.223 3.042 3.881 7

#8 1 2.540 3.973 5.067 8

#9 1 1/8 2.858 5.028 6.415 9

#10 1 1/4 3.175 6.207 7.917 10

#11 1 3/8 3.493 7.511 9.583 11

#12 1 1/2 3.810 8.938 11.401 12