Descargar

Algunas aplicaciones del MEF a problemas elípticos

Enviado por Ana Guerrero


Partes: 1, 2

    edu.red

    Introducción al Método de los Elementos Finitos 2 Algunas aplicaciones del MEF a problemas elípticos

    Problema de Elasticidad • Consideremos un cuerpo elástico, isótropo y homogéneo B, que ocupa el dominio acotado O?R3, con frontera • G=Gt?Gu tal que GtnGu=Ø y área Gu>0. El cuerpo B está sometido a una fuerza volumétrica f, •

    • y a una fuerza superficial t aplicada sobre Gt. Se supone B fijo a lo largo de Gu.

    Se busca determinar

    – desplazamientos u.

    – deformaciones e, dependientes de los desplazamientos de acuerdo a la cinemática de la deformación.

    – tensiones s, dependientes de las deformaciones de acuerdo a la ley constitutiva del material. O O Gu Gt

    edu.red

    ?u =0 sobre Gu 1? ?ui ?u j ? ??x ?xi ? ? = ??e d +µe ?,µ? , ctes. de Lamé. , µ = ?: coeficiente de Poisson 1+? (1+?)(1-2?) n =?s n Introducción al Método de los Elementos Finitos 3 Problema de Elasticidad • • Nota: en adelante, usaremos las siguientes convenciones de notación: Ecuación de equilibrio

    CB Dirichlet (despl. impuesto)

    CB Neumann (tracción impuesta) ? ?divs + f =0 en O ? ?sn = t sobre Gt O O Gu Gt n Ecuaciones de clausura

    – Cinemáticas: asumiendo pequeñas deformaciones:

    eij = ? + ? 2? j – Constitutivas: asumiendo comportamiento elástico lineal (ley de Hooke): , – derivada parcial: ui, j = ?ui ?xj 3 – sumatoria: sij j ij j j=1 +

    E E? E :módulo de Elasticidad 3 sij kk ij ij, k=1 con ? =

    edu.red

    { con V= v:v?? ? ? y v = 0 en Gu ?H (O) 1 ? ? dx+? ? -? 2 +v vi, j j,i ? – sij de pequeñas 2 -? sij ij ? ivids+? fividx = 0 Introducción al Método de los Elementos Finitos 4 Forma variacional del problema de Elasticidad • Dado se llega a haciendo } 3 (V) Hallar u?V/a(u,v) = L(v), ?v?V (D) sij, j + fi = 0 en O, i, j =1,2,3 O O G O O Gt O sij, jvidx+ fividx = 0 sijnjvids+ fividx = 0 -?Osijvi, j e (v)dx+ t d +µe t v ds+ O Gt i i O T. de Green Ley de Hooke CB dx+?Gt tivids+?O fividx = 0

    dx+?Gt tivids+?O fividx = 0 sijvi, j +s jivj,i O Simetría de s Cinemática

    O deformaciones L(v) a(u,v)

    edu.red

    i.e.,L(v) = ? v V ,?v?V,?? +. a(u,v) =? ??ui,i ij ij ? ij(v)dx =? ??divudivv +µeij ij(v)?dx continua, i.e., a(u,v) =? u v ,?u,v?V,? ? H (O) H (O) V-elíptica, i.e., a(v,v) =a v ,?v?V,a ? , con v + H (O) Demo.: a(v,v) = ?? (divv) dx+µ? eij(v)eij(v) dx H (O) =? vi H (O) = µ? eij ij(v) dx = µc v H (O) Introducción al Método de los Elementos Finitos 5 •

    • – –

    – . 1 1 2 2 . 3

    i-1 1 1 1

    2 es simétrica, i.e., a(u,v) = a(v,u),?u,v?V. Forma variacional del problema de Elasticidad (cont.)

    Se puede demostrar que la forma lineal L(v) =?O fividx+?G tivids es continua, t

    Se puede demostrar que la forma bilineal d +µe (u)?e (u)e O O 1 2 2 +

    . + O O

    O , c? (v)e Desigualdad de Korn

    edu.red

    { } Vh 1 (K)? ? ,?K?Th = v:v?V y v K ?? ?P H (O) Introducción al Método de los Elementos Finitos 6 •

    • 1 u-uh = Ch u H2(O) MEF aplicado al problema de Elasticidad

    Consideremos el problema de Elasticidad en O?R3. Sea Th={K} una malla de tetraedros de O. Definimos el espacio de EF 3

    El MEF aplicado al problema de Elasticidad consiste en Hallar uh ?Vh /a(uh,v) = L(v), ?v?Vh La solución uh?Vh satisface

    edu.red

    { } V = v:v?? ?H0 1(O)? ? y divv = 0 en O Introducción al Método de los Elementos Finitos 7 • Problema de Stokes

    Consideremos las ecuaciones de Stokes para el flujo estacionario de un fluido Newtoniano incompresible encerrado en un dominio O?R3, sometido a una fuerza volumétrica f : •

    • sij, j + fi = 0 en O, Balance de cant. de movto. en O, Ley const. de fluido Newtoniano ui,i = 0 en O, Condición de incompresibilidad ui = 0 sobre G, CB Dirichlet u: velocidad sij = 2µeij(u)- pdij s: tensión p: presión µ: viscosidad -µ?ui + p,i = fi en O, Balance de cant. de movto. p/fluido Newtoniano Definimos el espacio de funciones de prueba 3

    Luego, podemos llevar el problema de Stokes a la forma variacional (V) Hallar u?V/a(u,v) = L(v), ?v?V

    edu.red

    fi i v dx = µ? ?ui i v dx-µ? ?ui Introducción al Método de los Elementos Finitos 8 • •

    • Dado que µ >0, se demuestra (ídem problema de Poiss

    Partes: 1, 2
    Página siguiente