Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                
Scribd Logo
Cargar

¡Te damos la bienvenida a Scribd!

  • 238 vistas

    Marco Teorico

    Cargado por

    Rodrigonava
    Este documento describe el método de diferencias finitas para resolver problemas de placas. 1) El método fue propuesto en 1920 como una alternativa a la ecuación de Lagrange y plantea una serie de ecuaciones algebraicas lineales que definen los desplazamientos en puntos de la placa. 2) Se divide la placa en una malla y los desplazamientos en los vértices se consideran incógnitas de un sistema de ecuaciones. 3) El método permite aproximar derivadas mediante valores de la función en puntos vecinos y representar ecuaciones diferenciales

    Copyright:

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

    Marco Teorico

    Cargado por

    Rodrigonava
    238 vistas12 páginas
    Este documento describe el método de diferencias finitas para resolver problemas de placas. 1) El método fue propuesto en 1920 como una alternativa a la ecuación de Lagrange y plantea una serie de ecuaciones algebraicas lineales que definen los desplazamientos en puntos de la placa. 2) Se divide la placa en una malla y los desplazamientos en los vértices se consideran incógnitas de un sistema de ecuaciones. 3) El método permite aproximar derivadas mediante valores de la función en puntos vecinos y representar ecuaciones diferenciales

    Título original

    Marco teorico.docx

    Derechos de autor

    © © All Rights Reserved

    Formatos disponibles

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

    ¿Le pareció útil este documento?

    ¿Este contenido es inapropiado?

    Este documento describe el método de diferencias finitas para resolver problemas de placas. 1) El método fue propuesto en 1920 como una alternativa a la ecuación de Lagrange y plantea una serie de ecuaciones algebraicas lineales que definen los desplazamientos en puntos de la placa. 2) Se divide la placa en una malla y los desplazamientos en los vértices se consideran incógnitas de un sistema de ecuaciones. 3) El método permite aproximar derivadas mediante valores de la función en puntos vecinos y representar ecuaciones diferenciales

    Copyright:

    © All Rights Reserved
    Descargue como DOCX, PDF, TXT o lea en línea desde Scribd
    Descargar como docx, pdf o txt
    238 vistas12 páginas

    Marco Teorico

    Cargado por

    Rodrigonava
    Este documento describe el método de diferencias finitas para resolver problemas de placas. 1) El método fue propuesto en 1920 como una alternativa a la ecuación de Lagrange y plantea una serie de ecuaciones algebraicas lineales que definen los desplazamientos en puntos de la placa. 2) Se divide la placa en una malla y los desplazamientos en los vértices se consideran incógnitas de un sistema de ecuaciones. 3) El método permite aproximar derivadas mediante valores de la función en puntos vecinos y representar ecuaciones diferenciales

    Copyright:

    © All Rights Reserved
    Descargue como DOCX, PDF, TXT o lea en línea desde Scribd
    Descargar como docx, pdf o txt
    Está en la página 1de 12

    PROYECTO FINAL DE ELASTICIDAD CIV 308

    DIFERENCIAS FINITAS

    La solución de placas mediante diferencias finitas ha sido tradicionalmente un método muy


    empleado en la literatura.
    En las últimas décadas el auge del método de elementos finitos ha superado notablemente al
    desarrollo de soluciones utilizando diferencias finitas, pero ésta es aun una metodología
    importante conceptualmente y que puede programarse con mucha facilidad.

    Este método fue propuesto por Nielsen en 1920. Se trata de una solución alternativa a la
    ecuación de Lagrange. Se plantea mediante una serie de ecuaciones algebraicas lineales
    simultáneas, que definen la flecha de un número finito de puntos de la superficie de la placa.
    Una vez que se conocen las flechas de dichos puntos, se pueden obtener los momentos y los
    cortantes.

    Aplicación del cálculo por diferencias finitas. División de una placa mediante una malla
    cuadrada

    La utilización de medios informáticos permite la resolución de sistemas con un número


    significativamente grande de ecuaciones. Para su utilización en la práctica, se genera una
    malla que divide a la placa y se adapta a su contorno, considerando los desplazamientos en
    sus vértices como incógnitas de un sistema de ecuaciones lineales. En función de estos
    desplazamientos pueden expresarse las derivadas y, por tanto, los esfuerzos y la ecuación de
    equilibrio de la placa.

    Este método es de aplicación más general pero no permite abordar problemas geométricos
    complejos, como huecos, espesores variables o formas arbitrarias. Su uso, al igual que el de la
    ecuación diferencial, se limita generalmente a la creación de tablas de dimensionamiento
    rápido, para la obtención de esfuerzos en placas muy regulares apoyadas linealmente en sus
    contornos y bajo cargas uniformes o triangulares.

    1.4.1 Derivadas en función de diferencias finitas

    En una expansión en serie de Taylor se conoce el valor de la función f(x) en un punto


    determinado x = x0, e interesa conocer el valor de la función en un estado vecino
    caracterizado por x = x0 + Δx. Con referencia a la Figura 3.5.a, se puede escribir la expansión

    donde las derivadas se deben evaluar en el estado x = x0. Vale decir que si se conoce el valor
    de la función y de sus derivadas en un estado, se pueden obtener valores de la función para
    estados vecinos.

    Figura 3.5: (a) Función f(x); (b) Función f(x) en diferentes puntos x.

    Para valores de incrementos Δx pequeños, cada término no nulo en la expansión incide menos
    que el anterior. En situaciones prácticas, la ecuación (3.12) se trunca en determinado término,
    y el error que se comete en la aproximación es del orden del primer término siguiente que no
    fue incluido en la aproximación. Por ejemplo, se dice que el orden del error en la expresión
    (3.12) es O(Δx4).
    El método de diferencias finitas opera de manera inversa al desarrollo de Taylor: Se trata de
    averiguar cuánto valen las derivadas de una función a partir de valores de la función en
    determinados estados. Por ejemplo, supongamos que la expansión (3.12) se trunca luego del
    término Δx. Se tendrá así:
    Despejando la derivada primera se tiene:

    Vale decir que si se conoce el valor de la función f en los estados x = x0 y x = x0+Δx, entonces
    se puede evaluar la derivada de la función en x = x0. Una aproximación de este tipo, que
    utiliza valores de f en x0 y en x0 +Δx se denomina diferencias finitas hacia adelante.
    También puede obtenerse una aproximación de la derivada hacia atrás, y resulta en la forma

    Por último, pueden usarse expresiones de diferencias finitas centrales, en las que la derivada
    se aproxima como:

    En lugar de referirnos a las coordenadas x en cada caso, podemos identificar una secuencia de
    puntos en x y designarlos como n, n + 1, n + 2, etc., como se muestra en la Figura 3.5.b. Así se
    tendrán las derivadas primeras aproximadas como:
    Las derivadas de orden superior pueden computarse utilizado en forma recursiva las
    diferencias finitas hacia adelante y hacia atrás. Así, con diferencias hacia adelante y
    considerando que el espaciamiento Δx entre puntos que definen las derivadas es constante, se
    tiene:

    Sustituyendo por diferencias finitas hacia atrás


    se llega a

    Las derivadas cuartas pueden calcularse como

    Sustituyendo, se llega a:
    Derivadas parciales en función de diferencias finitas

    Cuando se trata de problemas bidimensionales, las derivadas se transforman en derivadas


    parciales en función de variables x1, x2. Supondremos un conjunto de puntos ubicados según
    una cuadrícula cuadrada o rectangular, para los cuales se encuentra definida una función f(x1,
    x2) (Figura 3.6.a). La grilla de puntos que se utiliza para evaluar derivadas en un punto es el
    conjunto de puntos vecinos que resultan afectados, y a este conjunto se lo denomina estrella.
    Una estrella en diferencias finitas es equivalente al concepto de elemento en el método de
    elementos finitos.

    Figura 3.6: (a) Función f(x1, x2); (b) Estrella utilizada para definir derivadas cuarta sobre una
    grilla ortogonal de puntos.

    Las derivadas direccionales no difieren de las derivadas en una dirección, pero las derivadas
    parciales cruzadas se deben calcular nuevamente. Utilizando la notación de la Figura 3.6.b, es
    fácil demostrar que las derivadas cuartas evaluadas en el punto 1 resultan:

    Las ecuaciones anteriores indican que para representar las derivadas cuartas se necesitan 13
    puntos, incluyendo al propio punto n = 1 en el que se desea realizar la evaluación. La
    numeración adoptada en la Figura 3.6.b es de tipo local, y sirve para identificar puntos con
    respecto a la posición del punto 1 al centro de la estrella.

    Representación de ecuaciones diferenciales en el dominio utilizando diferencias finitas

    Con la aproximación numérica de derivadas parciales es posible construir aproximaciones


    locales de ecuaciones diferenciales. Consideremos la ecuación diferencial en el dominio de
    una placa.

    donde p = p3 es la carga perpendicular a la superficie media de la placa y f = u3 es el


    desplazamiento transversal de la placa.
    Consideremos una malla de puntos igualmente espaciados en la direcciones x1 y x2, como se
    muestra en la Figura 3.6.b. Para cada punto de la malla es posible aproximar la ecuación
    diferencial mediante los valores de la función f en puntos vecinos. Se tiene así:

    La ecuación anterior puede simplificarse si se escribe uno de los espaciamientos de la grilla en


    función del otro:

    Y resulta:
    Es común escribir el operador de la ecuación en forma de molécula, cuyo centro corresponde
    al punto en el cual se obliga a cumplir la condición de diferencias finitas.

    En el caso particular que el espaciamiento en ambas direcciones sea el mismo

    Se tiene:

    Para el caso de malla cuadrada se tiene:

    Representación de condiciones de contorno utilizando diferencias finitas


    Las ecuaciones anteriores resultan adecuadas para puntos dentro del dominio de una placa
    alejados de los bordes. Pero a medida que nos acercamos a un borde los puntos necesarios no
    estarán disponibles porque caerían fuera de la placa. Hay dos maneras de tratar puntos
    cercanos al contorno: La primera es utilizando allí expresiones de diferencias finitas hacia un
    lado en lugar de centrales. La segunda es generar un conjunto de puntos auxiliares, ficticios,
    que quedarían ubicados fuera del dominio físico de la placa y servirían para dar apoyo a los
    operadores centrales
    desarrollados en esta sección. En esta presentación preferiremos la segunda opción por dos
    motivos: Sencillez de la formulación y evitar de tener errores diferentes por el empleo de
    operadores diferentes. Las condiciones de contorno de un problema de placas son también
    condiciones diferenciales, y se pueden representar de manera aproximada usando diferencias
    finitas. Cada tipo de contorno requiere de un tratamiento diferente:

     Contorno simplemente apoyado

    Las condiciones de este tipo de apoyo son desplazamiento nulo y momentos (o curvaturas)
    nulos en dirección perpendicular al contorno.

    (3.27)

    La primera condición es simple de aplicar, mientras que la segunda condición se puede


    representar utilizando diferencias finitas hacia un lado (hacia el interior de la placa). La
    condición de curvatura nula puede representarse también suponiendo un conjunto de puntos
    auxiliares que están fuera de la placa, tal como se indica en la Figura 3.7.a. Para ellos resulta

    (3.28)

    De manera que el punto auxiliar 3 tiene un valor de la función desplazamiento f igual y de


    signo contrario al de un punto 2 a igual distancia del borde pero hacia el interior.

     Contorno empotrado

    En un contorno empotrado se tiene


    (3.29)

    La condición de rotación nula puede representarse también suponiendo un conjunto


    de puntos auxiliares que están fuera de la placa, tal como se indica en la Figura
    3.7.b. Para ellos resulta:

    (3.30)

    Figura 3.7: Malla de puntos utilizados para (a) borde simplemente apoyado. (b)
    Borde empotrado. (c) Condicion de simetría

    De manera que el punto auxiliar 3 tiene un valor de la función desplazamiento f igual y del
    mismo signo que un punto a igual distancia del borde pero hacia el interior de la placa, tal
    como el punto 2.

     Contorno libre

    En un borde libre se tiene que el momento y el corte son nulos


    (3.31)

    Las expresiones que proporcionan los desplazamientos de los puntos auxiliares no son tan
    sencillas como las anteriores. Con referencia al arreglo de puntos de la Figura 3.6, si el punto 1
    se encuentra sobre el borde mismo, entonces los puntos ficticios tiene valores de la función
    f12, f3, f13 sobre la primera línea y f7 sobre la segunda. Se puede demostrar que las
    condiciones resultan en:

    Esto modifica la molécula de diferencias finitas. Para borde paralelo al eje x2 la molécula que
    corresponde a una distancia Δx1 desde el borde se representa como:

    La molécula que se aplica cuando el punto esta sobre el borde mismo es la siguiente:
    Condiciones de simetría

    Cuando existe una condición de simetría, entonces los desplazamientos a uno y otro lado de la
    línea de simetría son iguales. Con respecto a la Figura 3.7.c, se tiene

    (3.32)

    Solución de los desplazamientos en placas

    Las ecuaciones anteriores, que representan ecuaciones diferenciales de manera aproximada


    mediante una malla de puntos y diferencias finitas, pueden utilizarse para obtener la solución
    aproximada de los desplazamientos en puntos de una placa.

    Para una placa rectangular se deben realizar los siguientes pasos:

    1) Se considera una serie de puntos distribuidos sobre una malla rectangular, con
    espaciamiento constante en cada dirección, Δx1 y Δx2.

    2) Se aplica el operador que representa la ecuación diferencial en cada punto de la malla


    dentro del dominio de la placa. Se tendrá así un conjunto de condiciones en función de valores
    de la función desplazamiento en los puntos de la malla. Cerca del borde, los operadores
    requerirán el uso de puntos auxiliares que están fuera del dominio físico de la placa.

    3) Se aplican las condiciones de contorno, con las que se determinan los valores de los
    desplazamientos auxiliares fuera del dominio de la placa.
    4) Se soluciona un sistema de ecuaciones algebraicas con incógnitas que son desplazamientos
    en los nudos de la malla.

    5) Se evalúan los momentos en puntos de la placa a partir de los desplazamientos.

    Momentos en la placa a partir de los desplazamientos

    Una vez que se conocen los valores de los desplazamientos en los puntos seleccionados de la
    placas, es posible calcular cuanto valen los momentos utilizando las ecuaciones constitutivas y
    cinemáticas. Las relaciones entre momentos y desplazamientos son:

    Con referencia a la numeración interna de la posición de los puntos con respecto a un punto 1,
    adoptada en la Figura 3.6.b, se tienen las aproximaciones mediante diferencias finitas
    siguientes:

    También podría gustarte