Capítulo 1 Métodos de residuos ponderados Funciones de prueba ...
Capítulo 1 Métodos de residuos ponderados Funciones de prueba ... Capítulo 1 Métodos de residuos ponderados Funciones de prueba ...
La derivada segunda de u respecto a x dos veces (curvatura del eje medio) es ⎡ ⎤ u 1 u ′′ = d2 u d 2 ξ dξ 2 dx = 4 [ ] φ 1 2 L 2 ′ ξξ, ϕ 1′ ξξ, φ 2′ ξξ, ϕ 2′ ⎢ β 1 ⎥ ξξ ⎣ } {{ } u 2 ⎦ = B (ξ) u B(ξ) β 2 } {{ } u similarmente para la función de peso La integral ∫ v ′′ = d2 v dξ 2 d 2 ξ L dx = 4 [ φ 1 2 L 2 ′ ξξ, ϕ 1′ ξξ, φ 2′ ξξ, ϕ 2′ ξξ } {{ } B(ξ) ] ⎡ ⎤ v 1 ⎢ θ 1 ⎥ ⎣ v 2 ⎦ θ 2 } {{ } v = B (ξ) v d 2 v dx EI d2 u 2 dx ∫L dx = 2 vT B T (ξ) (EI) B (ξ) dx u = v T K u } {{ } K haciendo el cambio de variable en la integral (dx = L dξ e integrando entre -1 y 1) obtenemos K 2 ⎡ ⎤ 12 6L −12 6L K = EI 4L 2 −6L 2L 2 L 3 ⎢ sim. 12 −6L ⎥ ⎣ ⎦ 4L 2 4.3.1. Condiciones de Dirichlet no-homogéneas Resulta importante notar como se tratan las condiciones de contorno esenciales no homogéneas (desplazamiento prescriptos), sea por ej. (β 1 = ¯β 1 ). Recordemos que en el método de Galerkin las condiciones de contorno no homogéneas se satisfacían mediante una solución particular, en este caso la aproximación en el elemento resulta u (ξ) = φ 1 (ξ) u 1 + ϕ 1 (ξ) ¯β 1 + φ 2 (ξ) u 2 + ϕ 2 (ξ) β 2 donde ϕ 1 (ξ) ¯β 1 es ahora nuestra solución particular y el elemento queda entonces con sólo tres parámetros independientes. La matriz de rigidez del elemento queda ahora reducida a 3 × 3 y resulta de [ v 1 , v 2 , θ 2] ∫ ( ) ⎡ ⎤ ⎡ ⎤ φ 16 1′ ξξ ⎣ φ 2′ ⎦ L 4 ξξ (EJ) [ ] u 1 φ 1′ ξξ, φ 2′ ξξ, ϕ 2′ ξξ dx ⎣ u 2 ⎦ ϕ 2′ ξξ β 2 que es equivalente a eliminar la 2 fila y la segunda columna de la matriz completa ⎡ ⎤ 12 −12 6L K = EI L 3 ⎢ sim. 12 −6L ⎥ ⎣ ⎦ 4L 2 notar que para ello la función de peso se ha escrito ahora 62 v (ξ) = φ 1 (ξ) v 1 + φ 2 (ξ) v 2 + ϕ 2 (ξ) θ 2
pues debe satisfacer las condiciones homogéneas de contorno (θ = 0 donde β = ¯β ⇒ θ 1 = 0). La solución particular contribuye al término independiente de la forma [ v 1 , v 2 , θ 2] ∫ ( ) ⎡ ⎤ φ 1′ 16 ξξ ⎣ φ 2′ ⎦ L L 4 ξξ (EI) ϕ 1′ ξξ dx ¯β 1 ϕ 2′ ξξ que no es otra cosa que la segunda columna de la matriz original multiplicada por el valor conocido ¯β 1 . Ejercicio 1: Calcular los autovalores y autovectores de la matriz K Ejercicio 2: Calcular el término independiente debido a una carga uniforme 4.3.2. Formulación débil a partir del Principio de Trabajos Virtuales Si bien el método de residuos ponderados se presenta como una técnica numéricas consistente para resolver una ecuación diferencial con valores en el contorno, resulta ilustrativo y muy conveniente asociarla con principios físicos conocidos, de tal forma de poder dar una interpretación conceptual a distintos aspectos del método. Recordemos rápidamente el principio de trabajos virtuales, definamos primero un sistema virtual de desplazamientos δu. Conceptualmente un desplazamiento virtual es un incremento posible de desplazamientos a partir de una posición dada. Posible de ocurrir significa en este caso dos cosas 1. que satisfaga las ecuaciones de compatibilidad interna del problema, para el problema en estudio esto significa que la derivada segunda d2 u ds 2 exista (y sea finita). 2. que satisfaga las condiciones esenciales homogéneas del problema, es decir que tomado como desplazamiento incremental a partir de un campo de desplazamientos que satisface las condiciones esenciales no-homogéneas, la suma también satisfaga estas últimas. δu = 0 donde u = ū El principio de trabajos virtuales (desplazamientos virtuales) dice que un sistema de fuerzas internas (momentos y esfuerzos de corte en este caso) está en equilibrio con un conjunto de acciones externas (cargas y momentos) si para cualquier (es decir para todo) desplazamiento virtual se satisface que el trabajo virtual de las fuerzas internas es igual al de las fuerzas externas. Escribamos el principio de trabajos virtuales para el problema en estudio: ∫ T.V.I. = L ( ) d 2 ∫ δu M dx = δu q(x) dx + dx 2 L ( ) dδu dx ] L ¯M 0 reemplazando la ecuación constitutiva en el primer miembro, el T.V.I. resulta ∫ ( ) ( ) d 2 δu d 2 u T.V.I. = (EI) dx dx 2 dx 2 L + δu ¯Q ] L 0 = T.V.E. donde podemos reconocer la misma forma que la formulación débil obtenida por residuos ponderados con la diferencia que aquí v = δu. Si ahora comparamos las condiciones que impusimos a la función de peso v en residuos ponderados con las condiciones impuestas a los desplazamientos virtuales δu vemos que no hay ninguna diferencia y ambos formulaciones conducen a las mismas ecuaciones. Notar que estamos hablando de las ecuaciones de trabajos virtuales en forma discreta, es decir que buscamos la solución u dentro de una familia discreta de funciones (descripta por las funciones de forma utilizadas y sus parámetros asociados) y exigimos que se satisfaga la igualdad T.V.I. = T.V.E. para esa misma familia de funciones δu y no para cualquier función admisible. 63
- Page 9 and 10: La solución de los sistemas conduc
- Page 11 and 12: Una viga de longitud 1 está simple
- Page 13 and 14: Notar que K es no simétrica aun cu
- Page 15 and 16: integrando por partes el primer té
- Page 17 and 18: La sustitución de la aproximación
- Page 19 and 20: Figura 2 Ejercicio 21. a) Problema
- Page 21 and 22: Capítulo 2 El método de elementos
- Page 23 and 24: Figura 1 Aproximación de una funci
- Page 25 and 26: Figura 2 Comportamiento de tres fun
- Page 27 and 28: 2.4.1. Propiedades de la matriz K y
- Page 29 and 30: Figura 5 Descripción local y globa
- Page 31 and 32: 2.4.3.2. Cálculo del vector de car
- Page 33 and 34: de funciones. Una medida universalm
- Page 35: Figura 6 Curvas del error utilizand
- Page 38 and 39: eferiremos como ” ley de conserva
- Page 40 and 41: ambos miembros haciendo tender el p
- Page 42 and 43: 1. Condiciones de borde esenciales,
- Page 44 and 45: 2. Las condiciones de borde no pued
- Page 46 and 47: de tal modo que todo punto x en el
- Page 48 and 49: N 1 (ξ) = (ξ − ξ 2) (ξ 1 −
- Page 50 and 51: modo de ubicar un nodo en todo punt
- Page 52 and 53: Figura 5 Malla de N nudos y N-1 ele
- Page 54 and 55: 3. Condiciones naturales de Neumann
- Page 56 and 57: (c) Considere las condiciones (ii)
- Page 58 and 59: los polinomios de interpolación in
- Page 62 and 63: 4.3.3. Formulación a partir del Pr
- Page 64 and 65: ¯K = ∫ L 0 BT D B d¯x 1 = ¯K =
- Page 66 and 67: Las cargas normales al eje de la vi
- Page 68 and 69: 2. θ 1 (ξ) = ∑ 2 I=1 N I (ξ)
- Page 70 and 71: 4.5.7. Cambio de base La expresión
- Page 72 and 73: Figura 5 Viga con deformación de c
- Page 74 and 75: 4.6.1. Matriz de rigidez de una vig
- Page 76 and 77: Donde las funciones de forma son: N
- Page 78 and 79: (una para cada W J ) en función de
- Page 80 and 81: Ejercicio 1-Sea el problema de conv
- Page 82 and 83: o directamente las coordenadas noda
- Page 84 and 85: donde r es el residuo que se quiere
- Page 86: ⎡ K 3−4 = ⎢ ⎣ K 4−5 = ⎡
- Page 89 and 90: Figura 1 Conducción del calor en 2
- Page 91 and 92: 5.3. Forma variacional del problema
- Page 93 and 94: 5.5. Flujo en un medio poroso El fl
- Page 95 and 96: ortotropía coinciden con las direc
- Page 97 and 98: 5.6.5. Forma débil de la ecuación
- Page 99 and 100: Como φ es conocida sobre el contor
- Page 101 and 102: La formulación débil que resulta
- Page 103 and 104: 5.9.2. Formulación débil usando r
- Page 105 and 106: ⎡ σ = ⎢ ⎣ ⎤ σ 11 σ 22 σ
- Page 107 and 108: Figura 7 Teoría de placas clásica
- Page 109 and 110: Los esfuerzos de corte transversal
pues <strong>de</strong>be satisfacer las condiciones homogéneas <strong>de</strong> contorno (θ = 0 don<strong>de</strong> β = ¯β ⇒ θ 1 = 0). La<br />
solución particular contribuye al término in<strong>de</strong>pendiente <strong>de</strong> la forma<br />
[<br />
v 1 , v 2 , θ 2] ∫ ( ) ⎡ ⎤<br />
φ 1′<br />
16<br />
ξξ<br />
⎣ φ 2′ ⎦<br />
L L 4 ξξ (EI) ϕ 1′ ξξ dx ¯β 1<br />
ϕ 2′ ξξ<br />
que no es otra cosa que la segunda columna <strong>de</strong> la matriz original multiplicada por el valor conocido<br />
¯β 1 .<br />
Ejercicio 1: Calcular los autovalores y autovectores <strong>de</strong> la matriz K<br />
Ejercicio 2: Calcular el término in<strong>de</strong>pendiente <strong>de</strong>bido a una carga uniforme<br />
4.3.2. Formulación débil a partir <strong>de</strong>l Principio <strong>de</strong> Trabajos Virtuales<br />
Si bien el método <strong>de</strong> <strong>residuos</strong> pon<strong>de</strong>rados se presenta como una técnica numéricas consistente<br />
para resolver una ecuación diferencial con valores en el contorno, resulta ilustrativo y muy conveniente<br />
asociarla con principios físicos conocidos, <strong>de</strong> tal forma <strong>de</strong> po<strong>de</strong>r dar una interpretación<br />
conceptual a distintos aspectos <strong>de</strong>l método.<br />
Recor<strong>de</strong>mos rápidamente el principio <strong>de</strong> trabajos virtuales, <strong>de</strong>finamos primero un sistema virtual<br />
<strong>de</strong> <strong>de</strong>splazamientos δu. Conceptualmente un <strong>de</strong>splazamiento virtual es un incremento posible<br />
<strong>de</strong> <strong>de</strong>splazamientos a partir <strong>de</strong> una posición dada. Posible <strong>de</strong> ocurrir significa en este caso dos<br />
cosas<br />
1. que satisfaga las ecuaciones <strong>de</strong> compatibilidad interna <strong>de</strong>l problema, para el problema en<br />
estudio esto significa que la <strong>de</strong>rivada segunda d2 u<br />
ds 2 exista (y sea finita).<br />
2. que satisfaga las condiciones esenciales homogéneas <strong>de</strong>l problema, es <strong>de</strong>cir que tomado como<br />
<strong>de</strong>splazamiento incremental a partir <strong>de</strong> un campo <strong>de</strong> <strong>de</strong>splazamientos que satisface las<br />
condiciones esenciales no-homogéneas, la suma también satisfaga estas últimas.<br />
δu = 0 don<strong>de</strong> u = ū<br />
El principio <strong>de</strong> trabajos virtuales (<strong>de</strong>splazamientos virtuales) dice que un sistema <strong>de</strong> fuerzas<br />
internas (momentos y esfuerzos <strong>de</strong> corte en este caso) está en equilibrio con un conjunto <strong>de</strong> acciones<br />
externas (cargas y momentos) si para cualquier (es <strong>de</strong>cir para todo) <strong>de</strong>splazamiento virtual se<br />
satisface que el trabajo virtual <strong>de</strong> las fuerzas internas es igual al <strong>de</strong> las fuerzas externas.<br />
Escribamos el principio <strong>de</strong> trabajos virtuales para el problema en estudio:<br />
∫<br />
T.V.I. =<br />
L<br />
( ) d 2 ∫<br />
δu<br />
M dx = δu q(x) dx +<br />
dx 2 L<br />
( ) dδu<br />
dx<br />
] L<br />
¯M<br />
0<br />
reemplazando la ecuación constitutiva en el primer miembro, el T.V.I. resulta<br />
∫ ( ) ( )<br />
d 2 δu d 2 u<br />
T.V.I. =<br />
(EI) dx<br />
dx 2 dx 2<br />
L<br />
+ δu ¯Q ] L<br />
0 = T.V.E.<br />
don<strong>de</strong> po<strong>de</strong>mos reconocer la misma forma que la formulación débil obtenida por <strong>residuos</strong> pon<strong>de</strong>rados<br />
con la diferencia que aquí v = δu. Si ahora comparamos las condiciones que impusimos a<br />
la función <strong>de</strong> peso v en <strong>residuos</strong> pon<strong>de</strong>rados con las condiciones impuestas a los <strong>de</strong>splazamientos<br />
virtuales δu vemos que no hay ninguna diferencia y ambos formulaciones conducen a las mismas<br />
ecuaciones.<br />
Notar que estamos hablando <strong>de</strong> las ecuaciones <strong>de</strong> trabajos virtuales en forma discreta, es <strong>de</strong>cir<br />
que buscamos la solución u <strong>de</strong>ntro <strong>de</strong> una familia discreta <strong>de</strong> funciones (<strong>de</strong>scripta por las funciones<br />
<strong>de</strong> forma utilizadas y sus parámetros asociados) y exigimos que se satisfaga la igualdad T.V.I. =<br />
T.V.E. para esa misma familia <strong>de</strong> funciones δu y no para cualquier función admisible.<br />
63