UNIVERSIDAD DE HUANUCO E.A.

P INGENIERIA CIVIL

METODOS NUMERICOS Ing: Carlos Alberto Jara Trujillo

UNIVERSIDAD DE HUANUCO E.A.P INGENIERIA CIVIL

INTRODUCCION
Los métodos numéricos, en contraposición a los analíticos, se emplean para hallar
soluciones aproximadas. Esa es la diferencia clave. Se usan en los casos que no
resultan prácticos o posibles hallar una solución analítica, como por ejemplo para
resolver una integral de la que no existe o no se conoce la primitiva. Suelen ser
rápidos y sencillos de programar, siendo este su punto fuerte. A pesar de ser
aproximaciones, dado que pueden calcularse hasta obtener la precisión deseada,
proporcionan resultados perfectamente válidos para la mayoría de situaciones.
Existen, además, métodos muy generales aplicables a familias completas de casos.
Por ejemplo, el método de NewtonRhapson se puede utilizar para hallar raíces de la
mayoría de funciones, sin necesidad de programar sistemas que sepan distintos
métodos de resolución dependientes del caso. Desde tiempos ancestrales el papel del
ingeniero ha sido básicamente el mismo, tratar de conocer e interpretar los
mecanismos de la naturaleza para así poder modificarla al servicio del hombre. Para
ello ha utilizado sus conocimientos, intuición, experiencia y los medios naturales a los
que en cada momento ha tenido disponibles. Con el gran poder de cómputo que se
tiene en estos días, el ingeniero dispone de grandes ventajas para poder llevar a cabo
su misión y abordar cada día retos más ambiciosos en la solución de nuevos
problemas, cuyos aspectos políticos, económicos, científicos o tecnológicos pueden
tener un mayor impacto en la mejora de la calidad de vida del hombre. Encontramos
así aplicaciones de los métodos numéricos en los ámbitos más diversos desde
sectores tecnológicos tan clásicos como la ingeniería estructural o la aerodinámica de
aviones, hasta aplicaciones más sofisticadas como ingeniería de alimentos, ingeniería
médica, diseño de fármacos, biología, construcción o diseño de estructuras. En la
actualidad, gracias a la gran evolución que han tenido los métodos numéricos y su
implementación en potentes computadoras, es posible, por ejemplo, modelar el cálculo
de vigas o diseñas grandes urbanizaciones con cargas diferentes, permite calcular el
porcentaje de error.

METODOS NUMERICOS Ing: Carlos Alberto Jara Trujillo

UNIVERSIDAD DE HUANUCO E.A.P INGENIERIA CIVIL

OBJETIVOS
OBJETIVO GENERAL

 Realizar un estudio de los métodos numéricos y su aplicación en la ing. civil:

estructuras

OBJETIVOS ESPECIFICOS

 Ampliar los conocimientos generales de los métodos numéricos conocidos.

 Determinar la aplicación de los métodos numéricos a la Ingeniera Civil.

 Determinar la importancia de la utilización de los métodos numéricos.

 Identificar los modelos matemáticos que facilitan el diseño y cálculos de

estructuras.

METODOS NUMERICOS Ing: Carlos Alberto Jara Trujillo

UNIVERSIDAD DE HUANUCO E.A.P INGENIERIA CIVIL

MARCO TEORICO METODOS NUMERICOS:

Los métodos numéricos son técnicas mediante las cuales es posible formular
problemas matemáticos de tal forma que puedan resolverse usando operaciones
aritméticas. Hay muchos tipos de métodos numéricos, y comparten una característica
común: invariablemente se deben realizar un buen número de tediosos cálculos
aritméticos. Los métodos numéricos son herramientas muy poderosas para a solución
de problemas. Pueden manejar sistemas de ecuaciones grandes, no linealidades y
geometrías complicadas, comunes en la ingeniería. También es posible que se utilice
software disponible comercialmente que contenga métodos numéricos. El uso
inteligente de estos programas depende del conocimiento de la teoría básica de estos
métodos; además hay muchos problemas que no pueden plantearse al emplear
programas hechos, conociendo bien los métodos numéricos se puede diseñar
programas propios y así no comprar software costoso. Al mismo tiempo se aprende a
conocer y controlar los errores de aproximación que son inseparables de los cálculos
numéricos a gran escala. Los métodos numéricos son adecuados para la solución de
problemas comunes de ingeniería, ciencias y administración, utilizando computadoras
electrónicas. En el proceso de solución de problemas por medio de computadoras se
requieren los pasos siguientes:

 Especificación del problema: Con esto se indica que se debe identificar

perfectamente el problema y sus limitaciones, las variables que intervienen y
los resultados deseados.

 Análisis: es la formulación de la solución del problema denominada también

algoritmo, de manera que se tenga una serie de pasos que resuelvan el
problema y que sean susceptibles de ejecutarse en la computadora.

 Programación: Este paso consiste en traducir el método de análisis o algoritmo

de solución expresándole como una serie detallada de operaciones.

 Verificación: Es la prueba exhaustiva del programa para eliminar todos los

errores que tenga de manera que efectúe lo que desea los resultados de
prueba se comparan con soluciones conocidas de problemas ya resueltos.
 Documentación: Consiste en preparar un instructivo del programa de manera
que cualquier persona pueda conocer y utilizar el programa. Producción: Es la

METODOS NUMERICOS Ing: Carlos Alberto Jara Trujillo

 UNIVERSIDAD DE HUANUCO E.A.P INGENIERIA CIVIL

última etapa en la que solo se proporcionan datos de entrada del programa

obteniéndose las soluciones correspondientes.

1. La ecuacióni=9 e t sin 2 πt−3.5 , describe una corriente oscilatoria en un circuito

eléctrico, donde t se expresa en segundos. Determine todos los valores de t de modo
que i = 3. Con valores iniciales de [0,1] con un error de 10−4

Solucion

f ( x )=9 et sin 2 πt−3.5

Método bisección

i a b xr f(a) f(b) f(xr) f(a)*f(xr) error

1 0 1 0.5 -3.5 -3.5 -3.5 12.25  
2 0.5 1 0.75 -3.5 -3.5 -3.5 12.25 33.33333333
3 0.75 1 0.875 -3.5 -3.5 -3.5 12.25 14.28571429
4 0.875 1 0.9375 -3.5 -3.5 -3.5 12.25 6.666666667
5 0.9375 1 0.96875 -3.5 -3.5 -3.5 12.25 3.225806452
6 0.96875 1 0.984375 -3.5 -3.5 -3.5 12.25 1.587301587
7 0.984375 1 0.9921875 -3.5 -3.5 -3.5 12.25 0.787401575
8 0.9921875 1 0.99609375 -3.5 -3.5 -3.5 12.25 0.392156863
9 0.99609375 1 0.99804688 -3.5 -3.5 -3.5 12.25 0.195694716
10 0.99804688 1 0.99902344 -3.5 -3.5 -3.5 12.25 0.097751711
11 0.99902344 1 0.99951172 -3.5 -3.5 -3.5 12.25 0.048851979
12 0.99951172 1 0.99975586 -3.5 -3.5 -3.5 12.25 0.024420024
13 0.99975586 1 0.99987793 -3.5 -3.5 -3.5 12.25 0.012208522
14 0.99987793 1 0.99993896 -3.5 -3.5 -3.5 12.25 0.006103888
15 0.99993896 1 0.99996948 -3.5 -3.5 -3.5 12.25 0.003051851
16 0.99996948 1 0.99998474 -3.5 -3.5 -3.5 12.25 0.001525902
17 0.99998474 1 0.99999237 -3.5 -3.5 -3.5 12.25 0.000762945
18 0.99999237 1 0.99999619 -3.5 -3.5 -3.5 12.25 0.000381471
19 0.99999619 1 0.99999809 -3.5 -3.5 -3.5 12.25 0.000190735

Xr = 0.99999809 con un Et = 0.0001907

METODOS NUMERICOS Ing: Carlos Alberto Jara Trujillo

UNIVERSIDAD DE HUANUCO E.A.P INGENIERIA CIVIL

2. La ecuación de Ergun, que se da abajo, sirve para describir el flujo de un líquido a

través de un lecho empacado. ∆P es la caída de presión, r es la densidad del fluido, Gₒ
es la velocidad másica (el cociente del flujo de masa dividido entre el área de la sección
transversal), Dₑ es el diámetro de las partículas dentro del lecho, µ es la viscosidad del
fluido, L es la longitud del lecho y ϵ es la fracción vacía del lecho.

Dₑ 3
∗ε
∆ pρ L (1−ε )
=150
2
Gₒ (1−ε ) DₑGₒ
( )
μ

Dado los siguientes valores:

∆ pρ Dₑ
=¿ 10
Gₒ 2 L

DₑGₒ
=1000
μ

Mediante el método de la biseccion con un error en un intervalo de [1,2] con un error

de10−4.

Solucion:

Sustituyendo en los parámetros verdaderos tenemos

ε3 1−ε
10 =150 +1.75
1−ε 1000

Esto puede organizarse y expresarse como un problema de raíz.

ε3
f ( ε )=0.15 ( 1−ε ) +1.75−10 =0
1−ε

Grafica: Nos indica que existe raíz.

METODOS NUMERICOS Ing: Carlos Alberto Jara Trujillo

 UNIVERSIDAD DE HUANUCO E.A.P INGENIERIA CIVIL

ε3
Usando el método de la bisección. f ( ε )=0.15 ( 1−ε ) +1.75−10 =0
1−ε

Con rangos [1,2].

i a b xr f(a) f(b) f(xr) f(a)*f(xr) error

1 0 1 0.5 -3.5 -3.5 -3.5 12.25  
2 0.5 1 0.75 -3.5 -3.5 -3.5 12.25 33.33333333
3 0.75 1 0.875 -3.5 -3.5 -3.5 12.25 14.28571429
4 0.875 1 0.9375 -3.5 -3.5 -3.5 12.25 6.666666667
5 0.9375 1 0.96875 -3.5 -3.5 -3.5 12.25 3.225806452
6 0.96875 1 0.984375 -3.5 -3.5 -3.5 12.25 1.587301587
7 0.984375 1 0.9921875 -3.5 -3.5 -3.5 12.25 0.787401575
8 0.9921875 1 0.99609375 -3.5 -3.5 -3.5 12.25 0.392156863
9 0.99609375 1 0.99804688 -3.5 -3.5 -3.5 12.25 0.195694716
10 0.99804688 1 0.99902344 -3.5 -3.5 -3.5 12.25 0.097751711
11 0.99902344 1 0.99951172 -3.5 -3.5 -3.5 12.25 0.048851979
12 0.99951172 1 0.99975586 -3.5 -3.5 -3.5 12.25 0.024420024
13 0.99975586 1 0.99987793 -3.5 -3.5 -3.5 12.25 0.012208522
14 0.99987793 1 0.99993896 -3.5 -3.5 -3.5 12.25 0.006103888
15 0.99993896 1 0.99996948 -3.5 -3.5 -3.5 12.25 0.003051851
16 0.99996948 1 0.99998474 -3.5 -3.5 -3.5 12.25 0.001525902
17 0.99998474 1 0.99999237 -3.5 -3.5 -3.5 12.25 0.000762945
18 0.99999237 1 0.99999619 -3.5 -3.5 -3.5 12.25 0.000381471
19 0.99999619 1 0.99999809 -3.5 -3.5 -3.5 12.25 0.000190735

METODOS NUMERICOS Ing: Carlos Alberto Jara Trujillo

UNIVERSIDAD DE HUANUCO E.A.P INGENIERIA CIVIL

3. El volumen V del líquido contenido en un tanque esférico de radio r está relacionado

con la profundidad h del líquido por:

π k 2 (3 r−h)
V=
3

Determine h para r = 1 m y V = 0.75 m³, mediante el método de la bisección con un

error de 10−4. En un intervalo de [-1 , 2]

Solucion:

Entonces despejando tenemos.

π h2 (3−h)
f ( h )=0.5−
3

METODOS NUMERICOS Ing: Carlos Alberto Jara Trujillo

 UNIVERSIDAD DE HUANUCO E.A.P INGENIERIA CIVIL

Mediante el método de la bisección obtenemos la raíz.

i a b xr f(a) f(b) f(xr) f(a)*f(xr) error

-
1 -1 2 0.5 -3.6887902 -3.6887902 0.15449847 0.56991244  
-
2 0.5 2 1.25 0.15449847 -3.6887902 -2.3634308 0.36514644 93.4629578
-
3 1.25 2 1.625 -2.3634308 -3.6887902 3.30222705 7.80458512 28.4291852
-
4 1.625 2 1.8125 3.30222705 -3.6887902 -3.5852465 11.839298 7.89400269
-
5 1.8125 2 1.90625 -3.5852465 -3.6887902 3.66204141 13.1293212 2.09705201
- -
6 1.90625 2 1.953125 3.66204141 -3.6887902 3.68199515 13.4836187 0.54192723
- -
7 1.953125 2 1.9765625 3.68199515 -3.6887902 3.68707796 13.5758032 0.13785468
- -
8 1.9765625 2 1.98828125 3.68707796 -3.6887902 3.68836046 13.5992725 0.03477153
- -
9 1.98828125 2 1.99414063 3.68836046 -3.6887902 3.68868256 13.6051909 0.0087321
- -
10 1.99414063 2 1.99707031 3.68868256 -3.6887902 3.68876327 13.6066767 0.00218797
- -
11 1.99707031 2 1.99853516 3.68876327 -3.6887902 3.68878347 13.607049 0.00054762
- -
12 1.99853516 2 1.99926758 3.68878347 -3.6887902 3.68878852 13.6071421 0.00013698

Entonces la raíz de la función es : Xr = 1.99926758

con un error total de 0.00013%

METODOS NUMERICOS Ing: Carlos Alberto Jara Trujillo

UNIVERSIDAD DE HUANUCO E.A.P INGENIERIA CIVIL

4. En ingeniería estructural, la fórmula de la secante define la fuerza por unidad de área,

P/A, que ocasiona la tensión máxima σ m en una columna que tiene una razón de
esbeltez L/k dada es:

P σm
=
A ec
1+
( )
k2
sec ⁡¿ ¿

donde
ec/k2 = razón de excentricidad,
E = módulo de elasticidad.
Si para una viga de acero, E = 200 000 MPa, ec/k2 = 0.4 y σ m = 250 MPa, calcule P/A
para L/k = 50. Recuerde que sec x = 1/cos x, resolver con el método de la bisección con
un error de 10−4 con intervalos de [15, 25].

Solucion:

Despejamos “X” tenemos

P
X=
A

Tenemos la función:

40
f ( x )= −x
1 x
1+0.25 sec ⁡[
√
2 30000
∗100]

METODOS NUMERICOS Ing: Carlos Alberto Jara Trujillo

 UNIVERSIDAD DE HUANUCO E.A.P INGENIERIA CIVIL

Aplicamos bisección con la función : con los intervalos [15,25]

40
f ( x )= −x
1 x
1+0.25 sec ⁡[
√
2 30000
∗100]

i a b xr f(a) f(b) f(xr) f(a)*f(xr) error

1 15 25 20 25.5714923 16.9673246 20.9052551 534.578571  
2 20 25 22.5 20.9052551 16.9673246 18.7491921 391.956645 11.4994982
3 22.5 25 23.75 18.7491921 16.9673246 17.7815436 333.389578 5.4418701
4 23.75 25 24.375 17.7815436 16.9673246 17.3483844 308.481055 2.49682714
5 24.375 25 24.6875 17.3483844 16.9673246 17.1501015 297.526554 1.15616215
6 24.6875 25 24.84375 17.1501015 16.9673246 17.0565841 292.522148 0.54827731
7 24.84375 25 24.921875 17.0565841 16.9673246 17.0113955 290.156298 0.26563724
8 24.921875 25 24.9609375 17.0113955 16.9673246 16.9892168 289.010286 0.13054567
9 24.9609375 25 24.9804688 16.9892168 16.9673246 16.9782344 288.446906 0.064685
10 24.9804688 25 24.9902344 16.9782344 16.9673246 16.9727704 288.167674 0.03219301
11 24.9902344 25 24.9951172 16.9727704 16.9673246 16.9700452 288.02868 0.01605881
12 24.9951172 25 24.9975586 16.9700452 16.9673246 16.9686843 287.95934 0.00801994
13 24.9975586 25 24.9987793 16.9686843 16.9673246 16.9680043 287.924709 0.0040076
14 24.9987793 25 24.9993896 16.9680043 16.9673246 16.9676644 287.907403 0.00200321
15 24.9993896 25 24.9996948 16.9676644 16.9673246 16.9674945 287.898752 0.00100146
16 24.9996948 25 24.9998474 16.9674945 16.9673246 16.9674095 287.894428 0.00050069
17 24.9998474 25 24.9999237 16.9674095 16.9673246 16.9673671 287.892266 0.00025034
18 24.9999237 25 24.9999619 16.9673671 16.9673246 16.9673458 287.891185 0.00012517

METODOS NUMERICOS Ing: Carlos Alberto Jara Trujillo