15.07.2014 Views

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

SHOW MORE
SHOW LESS

You also want an ePaper? Increase the reach of your titles

YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.

El momento torsor se obtiene aplicando la ec. (1.39)<br />

T = 2<br />

∫ 3 ∫ 2<br />

−3<br />

−2<br />

ˆϕ dy dx = 74,26<br />

que pue<strong>de</strong> compararse con el valor exacto T = 76,4. El valor <strong>de</strong> la tensión <strong>de</strong> corte máxima resulta<br />

en |τ| = 3,02, mientras que el valor exacto es |τ| = 2,96.<br />

De los ejemplos presentados <strong>de</strong>staquemos nuevamente que el método <strong>de</strong> Galerkin produce<br />

(generalmente) sistemas <strong>de</strong> ecuaciones simétricos, a diferencia <strong>de</strong>l método <strong>de</strong> colocación que da<br />

matrices no simétricas. Cuando hay que resolver gran<strong>de</strong>s sistemas este aspecto resulta <strong>de</strong> singular<br />

importancia.<br />

1.4.1.2. Ejercicios<br />

Ejercicio N ◦ 4: a partir <strong>de</strong>l Ejemplo N ◦ 1<br />

a) Obtener la solución exacta <strong>de</strong> la ecuación diferencial planteada.<br />

b) Completar los <strong>de</strong>talles en la obtención <strong>de</strong> las expresiones generales <strong>de</strong> k lm y f l . Verificar los<br />

valores para el método <strong>de</strong> colocación y Galerkin. Para este último caso resulta útil recordar que<br />

I lm =<br />

∫ L<br />

0<br />

sin<br />

( ) lπx<br />

L<br />

sin<br />

( mπx<br />

)<br />

L<br />

{ L<br />

dx = 2<br />

, l = m<br />

0 l ≠ m<br />

c) Obtener una aproximación utilizando M = 1.<br />

d) Graficar la solución exacta y las aproximaciones obtenidas para M = 1 y 2 y los dos métodos<br />

<strong>de</strong> <strong>residuos</strong> pon<strong>de</strong>rados utilizados.<br />

En una tabla comparar los valores <strong>de</strong> la solución exacta y las aproximaciones obtenidas (con<br />

M = 2) en los puntos x = 1/3 y x = 2/3.<br />

Ejercicio N ◦ 5: a partir <strong>de</strong>l Ejemplo N ◦ 2<br />

a) Verificar las expresiones <strong>de</strong> k lm y f l en las ec. (1.40).<br />

b) Las curvas ϕ = cte tienen un significado importante en el problema que se trata ya que en<br />

ellas el módulo <strong>de</strong> la tensión <strong>de</strong> corte<br />

τ =<br />

√<br />

σ 2 zx + σ 2 zy =<br />

√ (∂ϕ<br />

∂y<br />

) 2 (<br />

+ − ∂ϕ ) 2<br />

∂x<br />

se mantiene constante y la dirección <strong>de</strong>l vector resultante coinci<strong>de</strong> con la tangente a dicha curva.<br />

La tangente <strong>de</strong> la función ϕ = cte se obtiene <strong>de</strong> la <strong>de</strong>rivada total <strong>de</strong> ϕ respecto a x<br />

dϕ<br />

dx = ∂ϕ<br />

∂x + ∂ϕ<br />

∂y<br />

dy<br />

dx = 0<br />

Entonces<br />

1) graficar las curvas resultante ϕ = cte que pasan por los puntos <strong>de</strong> abcisas: x = 0,5, 1, 1,5, 2, 2,5, 3<br />

y or<strong>de</strong>nada y = 0.<br />

2) Verificar que el valor <strong>de</strong> la tangente a la curva ϕ =cte. obtenida a partir <strong>de</strong> la anterior,<br />

coinci<strong>de</strong> con la pendiente que se obtiene como relación <strong>de</strong> las componentes <strong>de</strong> la tensión <strong>de</strong> corte<br />

τ.<br />

3) Obtenga el valor <strong>de</strong> la tensión <strong>de</strong> corte máxima haciendo pasar una parábola y encontrando<br />

su máximo.<br />

Ejercicio ( N ◦ 6: La distribución <strong>de</strong>l momento flector M (x) en una viga cargada con w (x) =<br />

πx<br />

)<br />

sin por unidad <strong>de</strong> longitud, satisface la ecuación<br />

L<br />

10<br />

d 2 M<br />

dx 2<br />

= w(x)

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

Saved successfully!

Ooh no, something went wrong!