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

La <strong>de</strong>rivada segunda <strong>de</strong> u respecto a x dos veces (curvatura <strong>de</strong>l eje medio) es<br />

⎡ ⎤<br />

u 1<br />

u ′′ = d2 u d 2 ξ<br />

dξ 2 dx = 4 [ ]<br />

φ<br />

1 2 L 2 ′ ξξ, ϕ 1′ ξξ, φ 2′ ξξ, ϕ 2′ ⎢ β 1<br />

⎥<br />

ξξ ⎣<br />

} {{ }<br />

u 2 ⎦ = B (ξ) u<br />

B(ξ)<br />

β 2<br />

} {{ }<br />

u<br />

similarmente para la función <strong>de</strong> peso<br />

La integral<br />

∫<br />

v ′′ = d2 v<br />

dξ 2 d 2 ξ<br />

L<br />

dx = 4 [<br />

φ<br />

1 2 L 2 ′ ξξ, ϕ 1′ ξξ, φ 2′ ξξ, ϕ 2′ ξξ<br />

} {{ }<br />

B(ξ)<br />

]<br />

⎡ ⎤<br />

v 1<br />

⎢ θ 1<br />

⎥<br />

⎣ v 2 ⎦<br />

θ 2<br />

} {{ }<br />

v<br />

= B (ξ) v<br />

d 2 v<br />

dx EI d2 u<br />

2 dx<br />

∫L<br />

dx = 2 vT B T (ξ) (EI) B (ξ) dx u = v T K u<br />

} {{ }<br />

K<br />

haciendo el cambio <strong>de</strong> variable en la integral (dx = L dξ e integrando entre -1 y 1) obtenemos K<br />

2<br />

⎡<br />

⎤<br />

12 6L −12 6L<br />

K = EI<br />

4L 2 −6L 2L 2<br />

L 3 ⎢ sim. 12 −6L<br />

⎥<br />

⎣<br />

⎦<br />

4L 2<br />

4.3.1. Condiciones <strong>de</strong> Dirichlet no-homogéneas<br />

Resulta importante notar como se tratan las condiciones <strong>de</strong> contorno esenciales no homogéneas<br />

(<strong>de</strong>splazamiento prescriptos), sea por ej. (β 1 = ¯β 1 ). Recor<strong>de</strong>mos que en el método <strong>de</strong> Galerkin las<br />

condiciones <strong>de</strong> contorno no homogéneas se satisfacían mediante una solución particular, en este<br />

caso la aproximación en el elemento resulta<br />

u (ξ) = φ 1 (ξ) u 1 + ϕ 1 (ξ) ¯β 1 + φ 2 (ξ) u 2 + ϕ 2 (ξ) β 2<br />

don<strong>de</strong> ϕ 1 (ξ) ¯β 1 es ahora nuestra solución particular y el elemento queda entonces con sólo tres<br />

parámetros in<strong>de</strong>pendientes. La matriz <strong>de</strong> rigi<strong>de</strong>z <strong>de</strong>l elemento queda ahora reducida a 3 × 3 y<br />

resulta <strong>de</strong><br />

[ v 1 , v 2 , θ 2] ∫ ( ) ⎡ ⎤<br />

⎡ ⎤<br />

φ<br />

16<br />

1′ ξξ<br />

⎣ φ 2′ ⎦<br />

L 4 ξξ (EJ) [ ]<br />

u 1<br />

φ 1′ ξξ, φ 2′ ξξ, ϕ 2′ ξξ dx ⎣ u 2 ⎦<br />

ϕ 2′ ξξ<br />

β 2<br />

que es equivalente a eliminar la 2 fila y la segunda columna <strong>de</strong> la matriz completa<br />

⎡<br />

⎤<br />

12 −12 6L<br />

K = EI<br />

L 3 ⎢ sim. 12 −6L<br />

⎥<br />

⎣<br />

⎦<br />

4L 2<br />

notar que para ello la función <strong>de</strong> peso se ha escrito ahora<br />

62<br />

v (ξ) = φ 1 (ξ) v 1 + φ 2 (ξ) v 2 + ϕ 2 (ξ) θ 2

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

Saved successfully!

Ooh no, something went wrong!