Scribd Logo
Cargar

¡Te damos la bienvenida a Scribd!

  • 58 vistas

    Trabajo Diferencias Finitas

    Cargado por

    Sergio Ballesteros
    El documento describe un problema de oscilaciones armónicas de un resorte con una masa unida. Se resuelve el problema analíticamente usando las ecuaciones diferenciales y condiciones de frontera, obteniendo la solución x(t)=10cos(4t). También se resuelve numéricamente usando el método de diferencias finitas, generando una matriz de ecuaciones que es resuelta en MATLAB. Los resultados numéricos muestran un error menor al 1% en comparación con la solución analítica.

    Copyright:

    © All Rights Reserved
    Descargue como ODT, PDF, TXT o lea en línea desde Scribd

    Trabajo Diferencias Finitas

    Cargado por

    Sergio Ballesteros
    58 vistas8 páginas
    El documento describe un problema de oscilaciones armónicas de un resorte con una masa unida. Se resuelve el problema analíticamente usando las ecuaciones diferenciales y condiciones de frontera, obteniendo la solución x(t)=10cos(4t). También se resuelve numéricamente usando el método de diferencias finitas, generando una matriz de ecuaciones que es resuelta en MATLAB. Los resultados numéricos muestran un error menor al 1% en comparación con la solución analítica.

    Descripción original:

    Diferencias finitas

    Derechos de autor

    © © All Rights Reserved

    Formatos disponibles

    ODT, PDF, TXT o lea en línea desde Scribd

    ¿Le pareció útil este documento?

    ¿Este contenido es inapropiado?

    El documento describe un problema de oscilaciones armónicas de un resorte con una masa unida. Se resuelve el problema analíticamente usando las ecuaciones diferenciales y condiciones de frontera, obteniendo la solución x(t)=10cos(4t). También se resuelve numéricamente usando el método de diferencias finitas, generando una matriz de ecuaciones que es resuelta en MATLAB. Los resultados numéricos muestran un error menor al 1% en comparación con la solución analítica.

    Copyright:

    © All Rights Reserved
    Descargue como ODT, PDF, TXT o lea en línea desde Scribd
    Descargar como odt, pdf o txt
    58 vistas8 páginas

    Trabajo Diferencias Finitas

    Cargado por

    Sergio Ballesteros
    El documento describe un problema de oscilaciones armónicas de un resorte con una masa unida. Se resuelve el problema analíticamente usando las ecuaciones diferenciales y condiciones de frontera, obteniendo la solución x(t)=10cos(4t). También se resuelve numéricamente usando el método de diferencias finitas, generando una matriz de ecuaciones que es resuelta en MATLAB. Los resultados numéricos muestran un error menor al 1% en comparación con la solución analítica.

    Copyright:

    © All Rights Reserved
    Descargue como ODT, PDF, TXT o lea en línea desde Scribd
    Descargar como odt, pdf o txt
    de 8

    ESCUELA DE INGENIERA MECNICA, UNIVERSIDAD INDUSTRIAL DE SANTANDER

    METODO DE LOS ELEMENTOS FINITOS


    28 DE OCTUBRE DE 2016, I SEMESTRE ACADMICO DE 2016

    SERGIO FERNANDO BALLESTEROS OSORIO


    201307
    TRABAJO DIFERENCIAS FINITAS

    PROBLEMA:
    El problema consiste en tirar hacia abajo una masa unida a un resorte de 10[cm] de longitud
    respecto a la posicin de equilibrio, sujetarla hasta que t=0 y soltarla desde el reposo en ese
    instante.
    Figura1. Esquema del problema

    Fuente: https://wilfridomtz.files.wordpress.com/2014/09/capitulo-5.pdf
    Ecuacin diferencial:
    d2 x
    +w 2 x=0
    d t2
    d2 x
    +16 x =0
    2
    dt

    Donde:

    Condiciones de contorno:
    x ( 1)=10
    x ' (1)=0

    1) ANALITICAMENTE:
    Al aplicar las condiciones de frontera:
    x ( t ) =C 1 C o s ( 4 t ) +C 2 S e n ( 4 t )

    Se obtiene:
    x ( 1 ) =10=C 1 * 1+ C 2 * 0

    Entonces:
    C1 =10

    Resultando:
    x ( t )=10 C o s ( 4 t ) +C 2 S e n ( 4 t )

    Como:
    x ( t )=40 S e n ( 4 t ) + 4 C 2 C o s ( 4 t )

    Se tiene:
    '
    x ( 1 )= 0= 4 C 2

    As que:
    C 2=0

    Obteniendo la ecuacin para el movimiento del sistema resulta:


    x ( t )=10 C o s ( 4 t )

    2) DIFERENCIAS FINITAS

    d2 x

    d t2

    x k +1 x k
    x k x k 1

    h
    h
    h

    d x x k+1 2 x k + x k1

    d t2
    h2
    2

    Reemplazando en la ecuacin diferencial inicial:

    xk +1 2 x k + x k 1
    d x
    2
    2
    +w x=
    + w x k =0
    2
    2
    dt
    h
    2

    Si se ordenase:
    1
    2
    1
    x k 1 +
    +w 2 x k + 2 x k +1= 0
    h2
    h2
    h

    Para compactar la expresin anterior:


    M x k1 + P x k + M x k+1= 0

    Donde:
    M=

    1
    h2

    P=

    2
    2
    +w
    h2

    Para las condiciones de contorno:


    y

    t 1=0

    t N =b

    Obtendramos de la siguiente manera el valor de h que corresponde a la medida de las


    divisiones:

    h=

    b
    N 1

    donde N corresponde el nmero de divisiones a realizar

    Para realizar la matriz de ecuaciones, se toma un punto anterior a

    t1 ,

    entendiendo el concepto

    que la derivada de una ecuacin en un punto, es la recta tangente en dicho punto:

    '

    x ( 1 )=

    x 2 x0
    2h

    '
    x 0= x 2 2 h x ( 1 )

    De esta manera, podemos encontrar las ecuaciones que rigen el sistema:


    Si k=1:
    Si k=2:
    Si k=3:

    M x 0 + P x 1+ M x 2 =0
    M x 1 + P x 2 + M x 3 =0
    M x 2 + P x 3 + M x 4=0

    Si k=k

    M x k1 + P x k + M x k+1=0

    Sustituyendo:
    x ( 1)=10
    '
    x 0= x 2 2 h x ( 1 )

    ( x 2 2 h x' ( 1 ) ) M + 10 P+ M x 2=0
    10 M + P x 2 + M x 3 =0
    M x 2 + P x 3 + M x 4=0
    M x k1 + P x k + M x k+1=0

    Despejando las constantes:

    2 M x 2 =( 2 h x ' ( 1 ) ) M 10 P
    P x 2+ M x 3 =10 M
    M x 2 + P x 3 + M x 4=0
    M x 3 + P x 4 + M x 5=0
    M x k1 + P x k + M x k+1= 0

    Obteniendo finalmente la matriz de ecuaciones de manera general:

    2M

    x2

    x3

    x4

    M
    0
    0

    P
    M
    0

    M
    P
    M

    SOLUCIONADO CON MATLAB:

    x5

    xk

    (2hx'(1))M10P
    -10P
    =

    0
    0
    0

    clc
    close all
    clear all
    b=6;
    N=800;
    h=b/(N-1);
    Omega2=16;
    alfa=10;
    alfaprima=0;
    M=1/h^2;
    P=Omega2-2/(h^2);
    A=zeros(N-1,N-1);
    for contador=1:N-1
    A(contador,contador)=M;
    if(contador<N-1)
    A(contador+1,contador)=P;
    end
    if(contador<N-2)
    A(contador+2,contador)=M;
    end
    end
    A(1,1)=2*M;
    ladoderecho=zeros(N-1,1);
    ladoderecho(1)=2*alfaprima*h*M-alfa*P;
    ladoderecho(2)=-alfa*M;
    solucion=inv(A)*ladoderecho;
    x=[alfa;solucion];
    t=(0:h:b)';
    plot(t,x)

    RESULTADOS OBTENIDOS:

    a) Analtico:
    Mediante la ecuacin:
    x ( t )=10 C o s ( 4 t )

    t
    0
    0,0075
    0,015
    0,0225
    0,03
    0,0375
    0,045
    0,0525
    0,06
    0,0675

    x
    10
    9,995500
    337
    9,982005
    399
    9,959527
    33
    9,928086
    359
    9,887710
    779
    9,838436
    928
    9,780309
    147
    9,713379
    749
    9,637708
    964

    0,075
    0,0825

    b) Diferencias finitas

    9,553364
    891
    9,460423
    435

    Porcentaje de error entre el mtodo exacto y el mtodo mediante diferencias finitas

    t
    0
    0,0075
    0,015
    0,0225
    0,03
    0,0375
    0,045
    0,0525
    0,06

    x Exacta
    10
    9,995500
    337
    9,982005
    399
    9,959527
    33
    9,928086
    359
    9,887710
    779
    9,838436
    928
    9,780309
    147
    9,713379
    749

    x
    Diferenci
    as
    10
    9,995488
    729
    9,981958
    986
    9,959422
    978
    9,927901
    039
    9,887421
    609
    9,838021
    211
    9,779744
    417
    9,712643
    807

    % Error
    0
    0,000116
    138
    0,000464
    972
    0,001047
    761
    0,001866
    621
    0,002924
    544
    0,004225
    435
    0,005774
    154
    0,007576
    572

    0,0675
    0,075
    0,0825

    ENLACE ARCHIVO .m
    https://mega.nz/#!eBlSTIKD

    9,637708
    964
    9,553364
    891
    9,460423
    435

    9,636779
    924
    9,552221
    215
    9,459043
    974

    0,009639
    635
    0,011971
    454
    0,014581
    395

    También podría gustarte