VIGA SIMPLEMENTE APOYADA CON CARGA UNIFORMEMENTE
DISTRIBUIDACOMPARACION DE LA SOLUCION NUMERICA DE LA
FLECHA MAXIMA FRENTE A LA SOLUCION ANALITICA
1. INTRODUCCIÓN
Dentro del campo de la ingeniería existe la famosa viga simplemente apoyada con carga
uniformemente distribuida.
La cual tiene la siguiente representación:
Esta viga simplemente apoyada representa distintas situaciones reales como por ejemplo:
Una viga de hormigón armado.
Viga de hormigón armado
La carga distribuida en la viga simplemente apoyada provoca la presencia de reacciones de
apoyo, momentos y la flecha máxima de la viga.
La flecha máxima es un fenómeno bastante estudiado dentro la ingeniería, puesto que por
las características del acero este tiende a romperse cuando sobre pasa a su límite de
fluencia. Por esta razón es muy importante el estudio de la flecha máxima en vigas.
En la resolución analítica de una ecuación diferencial muchas veces se hace un proceso
largo y engorroso, algunas de estas no admiten una solución analítica.
Los métodos numéricos nos dan una excelente aproximación a la solución de una ecuación
diferencial de problemas de contorno.
2. PLANTEAMIENTO DEL PROBLEMA
Se resolverá numéricamente la ecuación diferencial que representa a la línea elástica o
flecha de una viga simplemente apoyada con carga uniformemente distribuida.
Con los datos:
L=600 [cm ]
q=15 [kg/cm]
E=215 000 [kg/cm2] modulo de elasticidad del concreto
Sección de viga:
b= 20 cm base de la viga
h= 40 cm altura de la viga
I= 106666.6667 [cm4] momento de inercia de la viga
Para esta viga simplemente apoyada con carga uniformemente distribuida se desea saber el
valor de la flecha máxima la cual se encuentra en L/2=300 [cm]
La ecuación diferencial de la línea elástica o flecha está dada por:
La flecha máxima se presenta en L/2 es decir en L=300 cm
También es necesario realizarse una pregunta de investigación la cual es:
¿Qué aproximación tendrá la solución numérica de la flecha máxima de una viga
simplemente apoyada con carga uniformemente distribuida frente a la solución analítica?
3. OBJETIVOS
3.1. OBJETIVO GENERAL
Determinar la flecha máxima de una viga simplemente apoyada con carga
uniformemente distribuida con el método numérico Runge Kutta de segundo
orden.
3.2. OBJETIVOS ESPECIFICOS
Realizar un programa en Matlab GUIDE del método de Runge kutta de
segundo orden que permita la solución numérica de la flecha máxima en una
viga simplemente apoyada con carga uniformemente distribuida.
Validación del método numérico de Runge Kutta de segundo Orden
comparando con la solución analítica de la flecha máxima de una viga
simplemente apoyada con carga uniformemente distribuida.
4. MARCO TEORICO DEL PROBLEMA
DEFORMACIÓN EN VIGAS
Denominaremos línea elástica a la curva que forma la fibra neutra una vez cargada la viga,
considerando que ésta se encontraba inicialmente recta.
SUPUESTOS BASE.
Para establecer una serie de relaciones al interior de la sección, indicamos que se trata de
una viga, cuyo material se encuentra solicitado dentro del rango de proporcionalidad entre
tensiones y deformaciones, y en donde se admite la conservación de las caras planas. Dicho
en otra forma, donde se cumplen la ley de Hooke y la hipótesis de Bernouilli-Navier.
LEY DE HOOKE.
Establece que la relación entre la tensión y la deformación unitaria es una constante y se
denomina módulo de elasticidad.
La ecuación diferencial de la línea elástica o flecha de una viga viene dada por:
Donde:
M: momento a lo largo de la longitud en función de x
E: módulo de elasticidad concreto
I: momento de inercia de la sección de la viga
Se tiene la siguiente viga simplemente apoyada con carga uniformemente distribuda:
En la cual la ecuación de la línea elástica o flecha máxima esta presentada por la siguiente
ecuación diferencial:
Cuya solución analítica de la ecuación diferencial es:
Flecha máxima se produce cuando la cual es:
5. MARCO TEORICO DEL METODO NUMERICO
El método asociado con los nombres Runge (1885), Kutta (1901), consisten en obtener un
resultado al que se podría llegar al utilizar un número finito de términos de una serie de
Taylor de la forma:
con una aproximación en la cual se calcula de una fórmula del tipo:
A continuación se derivará sólo el caso más simple, cuando p = 1, para ilustrar el
procedimiento del caso general, ya que los lineamientos son los mismos.
Por otro lado, b puede manejarse más libremente y expresarse yi + bh , sin pérdida de
generalidad, como una ordenada arriba o abajo de la ordenada que da el método de Euler
simple
Con
Para obtener los parámetros desconocidos, se expande primero f (xi + μh, yi + λk0) en
serie de Taylor (obviamente mediante el desarrollo de Taylor de funciones de dos
variables).
Se sustituye en la ecuación anterior quedando:
Esta última ecuación se arregla en potencias de h, y queda
Para que los coeficientes correspondientes de coincidan se requiere:
Resolviendo estas ecuaciones se obtiene:
Reemplazamos estos valores obteniendo:
Conocida como algoritmo de Runge Kutta de segundo orden punto medio.
6. DESARROLLO
Se desea resolver la siguiente ecuación diferencial por el método de Runge Kutta segundo
Orden Punto Medio
La cual representa la ecuacion diferencial de la línea elástica de una viga
Para resolver el problema se hará uso del software Matlab en el cual se realizó un scrip:
syms x y z
clc
format long g
f=z
g=1/(215000*20*40^3/12)*(4500*x-7.5*x^2+0.000000001*y+0.000000001*z)
x0=0
x1=300
y0=0
z0=-15*600^3/(24*215000*20*40^3/12)
f=inline(f);
g=inline(g);
n=10
i=1;
h=(x1-x0)/n
while i<=n
disp('+++++++++++++++++++++++++++++++++++++++++++++++++++++++')
A=[x0 y0 z0]
k1=h*f(z0)
l1=h*g(x0,y0,z0)
B=[x0+h/2 y0+k1/2 z0+l1/2]
k2=h*f(z0+l1/2)
l2=h*g(x0+h/2,y0+k1/2,z0+l1/2)
C=[x0+h/2 y0+k2/2 z0+l2/2]
k3=h*f(z0+l2/2)
l3=h*g(x0+h/2, y0+k2/2, z0+l2/2)
D=[x0+h y0+k3 z0+l3]
k4=h*f(z0+l3)
l4=h*g(x0+h, y0+k3, z0+l3)
K=[k1;k2;k3;k4]
L=[l1;l2;l3;l4]
y1=y0+1/6*(k1+2*k2+2*k3+k4)
z1=z0+1/6*(l1+2*l2+2*l3+l4)
x1=x0+h
x(i)=x1;
y(i)=y1;
z(i)=z1;
x0=x1;
y0=y1;
z0=z1;
i=i+1;
disp('++++++++++++++++++++++++++++++++++++++++++++++++++++++++')
disp(' iteracion nueva')
end
disp('la solucion sera')
w=[x' y' z'];
w=vpa(w,9)
plot(x,y)
PRESENTACION DE RESULTADOS
Haciendo correr el script de Matlab nos da como resultado la flecha máxima numerica:
Resultado analítico de la flecha máxima:
Para un análisis del error Vs número de iteraciones se calculó distintos valores de la flecha
máxima numérica
N flecha
numero de Flecha maxima maxima error
iteraciones analitica calculada Error porcentual
2 -1.103742733 -1.2279138 0.12417106 -11.25
-
3 -1.103742733 -1.152798 0.04905524 4.44444485
4 -1.103742733 -1.1296117 0.02586897 -2.3437497
-
5 -1.103742733 -1.1196366 0.0158939 1.44000015
-
6 -1.103742733 -1.1144736 0.01073083 0.97222176
-
7 -1.103742733 -1.1114657 0.00772298 0.69970807
-
8 -1.103742733 -1.1095633 0.00582052 0.52734363
-
9 -1.103742733 -1.1082849 0.00454215 0.41152226
-
10 -1.103742733 -1.1073851 0.00364235 0.32999964
-
11 -1.103742733 -1.1067281 0.00298533 0.27047308
-
12 -1.103742733 -1.1062338 0.00249109 0.22569453
-
13 -1.103742733 -1.1058528 0.00211003 0.19117018
-
20 -1.103742733 -1.1046119 0.0008692 0.07874996
-
25 -1.103742733 -1.1042937 0.00055099 0.04991988
-
30 -1.103742733 -1.1041229 0.00038018 0.03444435
-
40 -1.103742733 -1.1039549 0.00021213 0.01921888
50 -1.103742733 -1.1038778 0.0001351 -0.0122399
-
60 -1.103742733 -1.1038362 9.3507E-05 0.00847181
-
80 -1.103742733 -1.1037951 5.2387E-05 0.00474631
-
100 -1.103742733 -1.1037762 3.3447E-05 0.00303033
-
150 -1.103742733 -1.1037576 1.4817E-05 0.00134243
-
200 -1.103742733 -1.1037511 8.317E-06 0.00075353
-
500 -1.103742733 -1.1037441 1.327E-06 0.00012023
Cuyo grafico es:
NUMERO DE ITERACIONES VS ERROR
0
0 100 200 300 400 500 600
-2
-4
-6
-8
-10
-12
En la iteración numero 60 presenta un error =0.00009
CONCLUSIONES
Se validó el resultado de la flecha máxima numérica con el resultado de la flecha
máxima analítica con n=60 número de iteraciones ya se alcanza un
error=9.3507E-05 la cual ya es aceptable.
Los métodos numéricos aplicados a problemas de la ingeniería brindan excelentes
aproximaciones hacia la solución real, por esta razón es muy importante su estudio.
El método de Runge Kutta punto medio presenta una convergencia cuadrática.
RECOMENDACIONES
Se recomienda realizar el cálculo numérico con el método de Runge Kutta de cuarto orden
y comparar la velocidad de convergencia con el método ya visto.
BIBLIOGRAFIA
Chapra (2005). Métodos Numéricos para [Link]. Mc GrawnHill.
