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 ...

efn.uncor.edu
from efn.uncor.edu More from this publisher
15.07.2014 Views

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

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

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!