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

2. θ 1 (ξ) = ∑ 2 I=1 N I (ξ) θ I 1 3. u 2 (ξ) = φ 1 (ξ) u 1 2 + ϕ 1 (ξ) θ 1 3 + φ 2 (ξ) u 2 2 + ϕ 2 (ξ) θ 2 3 4. u 3 (ξ) = φ 1 (ξ) u 1 3 − ϕ1 (ξ) θ 1 2 + φ2 (ξ) u 2 3 − ϕ2 (ξ) θ 2 2 Donde puede verse que hemos usado para el comportamiento axial los polinomios de Lagrange de grado 1 y para el comportamiento flexional los polinomios de Hermite de orden 1. 4.5.5. Desarrollo Recordando que d = d dξ = d 2 dx 1 dξ dx 1 dξ L ( ) −1 dξ dx1 = = J −1 dx 1 dξ tenemos que las deformaciones generalizadas resultan ε = du 1 = N I 2 dx ′ ξ 1 L uI 1 = ( ) u 2 1 − u 1 1 1 L χ 1 = dθ 1 = N I 2 dx ′ ξ 1 L θI 1 = ( θ 2 1 − ) 1 θ1 1 L χ 3 = d2 u 2 dx 2 1 = [ 6ξ L ( ) = d2 u 2 dξ 2 dξ 2 = 4 [ ] φ 1 1 dx 1 L 2 ′ ξξu 1 2 + ϕ 1′ ξξθ 1 3 + φ 2′ ξξu 2 2 + ϕ 2′ ξξθ 2 3 ⎡ ⎤ −1 + 3ξ , , − 6ξ 2 L L , 1 + 3ξ ] u 1 2 ⎢ θ 1 3 ⎥ 2 L ⎣ u 2 ⎦ 2 θ 2 3 χ 2 = − d2 u 3 dx 2 1 [ 6ξ = − L ( ) 2 = − d2 u 3 dξ dξ 2 = − 4 [ φ 1 1 dx 1 L 2 ′ ξξ u1 3 − ϕ1′ ξξ θ1 2 + φ2′ ξξ u2 3 − ] ϕ2′ ξξ θ2 2 ⎡ ⎤ + 3ξ , −−1 , − 6ξ ] u 1 3 + 3ξ , −1 ⎢ θ 1 2 ⎥ 2 L L2 L ⎣ u 2 ⎦ 3 θ 2 2 ⎡ ⎢ ⎣ Que puede resumirse en la siguiente expresión matricial ⎤ ε χ 1 χ 3 χ 2 ⎥ ⎦ = 1 L ⎡ ⎢ ⎣ −1 +1 −1 +1 6ξ −1 + 3ξ − 6ξ L L −1 + 3ξ − 6ξ L 6ξ L +1 + 3ξ +1 + 3ξ ⎡ ⎤ ⎥ ⎦ ⎢ ⎣ u 1 1 u 1 2 u 1 3 θ 1 1 θ 1 2 θ 1 3 u 2 1 u 2 2 u 2 3 θ 2 1 θ 2 2 θ 2 3 ⎤ ⎥ ⎦ 70

}{{} ε =B (ξ) } {{ } }{{} a e 4×1 4×12 12×1 En tanto que la relaciones constitutivas pueden escribirse matricialmente como ⎡ ⎤ ⎡ ⎤ ⎡ ⎤ N EA ε ⎢ M t ⎥ ⎣ M 3 ⎦ = ⎢ αGI p ⎥ ⎢ χ 1 ⎥ ⎣ EI 3 ⎦ ⎣ χ 3 ⎦ M 2 EI 2 χ 2 σ }{{} 4×1 = D }{{} 4×4 ε }{{} 4×1 Finalmente la matriz de rigidez (en el sistema local) resulta de integrar a lo largo de la viga K }{{} L = 12×12 ∫ L 0 B T }{{} 12×4 D }{{} 4×4 }{{} B dx 1 = 4×12 ∫ 1 −1 B (ξ) T D B (ξ) L 2 dξ Notar que como máximo en B(ξ) hay polinomios de orden 1, luego en el producto B T D B hay como máximo polinomios de orden 2 en ξ, por lo cual bastan dos puntos de integración si se va a realizar una integración numérica. La integral puede hacerse en este caso en forma analítica sin ningún problema. La matriz de rigidez (cuya obtención se deja como ejercicio) es idéntica a la obtenida en los cursos de cálculo matricial de estructuras. Esto es así porque las funciones de interpolación utilizadas son capaces de reproducir la solución exacta de las ecuaciones diferenciales (homogéneas). 4.5.6. Término independiente (vector de cargas) Supongamos que las cargas están representadas localmente (es decir respecto al sistema local) y varían en forma lineal dentro del elemento ⎡ q (ξ) = ⎣ q ⎤ 1 2∑ q 2 ⎦ = N I (ξ) q I q 3 I=1 El trabajo virtual externo se puede escribir como ⎡ ∫ L ∫ L T.V.E. = δu T q (ξ) dx 1 = [δu 1 , δu 2 , δu 3 ] ⎣ q ⎤ 1 (ξ) q 2 (ξ) ⎦ dx 1 0 0 q 3 (ξ) ⎡ ⎤ ∫ L N 1 N 2 T = [δa e ] T ⎣ φ 1 ϕ 1 φ 2 ϕ 2 ⎦ 0 φ 1 0 −ϕ 1 φ 2 0 −ϕ 2 Donde F L × ⎡ = [δa e ] T } {{ } 1×12 ⎣ N 1 N 2 N 1 N 2 N 1 ∫ L 0 Φ T (ξ) } {{ } 12×3 N (ξ) } {{ } 3×6 dx 1 [ q 1 q 2 ] } {{ } 6×1 ⎡ ⎤ ⎦ dx 1 N 2 ⎢ ⎣ q 1 1 q 1 2 q 1 3 q 2 1 q 2 2 q 2 3 ⎤ ⎥ ⎦ = [δa e ] T F } {{ } }{{} L 1×12 12×1 es el término independiente referido al sistema de coordenadas locales. (4.3) 71

2. θ 1 (ξ) = ∑ 2<br />

I=1 N I (ξ) θ I 1<br />

3. u 2 (ξ) = φ 1 (ξ) u 1 2 + ϕ 1 (ξ) θ 1 3 + φ 2 (ξ) u 2 2 + ϕ 2 (ξ) θ 2 3<br />

4. u 3 (ξ) = φ 1 (ξ) u 1 3 − ϕ1 (ξ) θ 1 2 + φ2 (ξ) u 2 3 − ϕ2 (ξ) θ 2 2<br />

Don<strong>de</strong> pue<strong>de</strong> verse que hemos usado para el comportamiento axial los polinomios <strong>de</strong> Lagrange<br />

<strong>de</strong> grado 1 y para el comportamiento flexional los polinomios <strong>de</strong> Hermite <strong>de</strong> or<strong>de</strong>n 1.<br />

4.5.5. Desarrollo<br />

Recordando que<br />

d<br />

= d dξ<br />

= d 2<br />

dx 1 dξ dx 1 dξ L<br />

( ) −1<br />

dξ dx1<br />

=<br />

= J −1<br />

dx 1 dξ<br />

tenemos que las <strong>de</strong>formaciones generalizadas resultan<br />

ε = du 1<br />

= N I 2<br />

dx ′ ξ<br />

1 L uI 1 = ( )<br />

u 2 1 − u 1 1<br />

1<br />

L<br />

χ 1 = dθ 1<br />

= N I 2<br />

dx ′ ξ<br />

1 L θI 1 = ( θ 2 1 − ) 1<br />

θ1 1<br />

L<br />

χ 3 = d2 u 2<br />

dx 2 1<br />

=<br />

[ 6ξ<br />

L<br />

( )<br />

= d2 u 2 dξ 2<br />

dξ 2 = 4 [ ]<br />

φ<br />

1<br />

1<br />

dx 1 L 2 ′ ξξu 1 2 + ϕ 1′ ξξθ 1 3 + φ 2′ ξξu 2 2 + ϕ 2′ ξξθ 2 3<br />

⎡ ⎤<br />

−1 + 3ξ<br />

, , − 6ξ<br />

2<br />

L L , 1 + 3ξ ]<br />

u 1 2<br />

⎢ θ 1 3 ⎥<br />

2 L ⎣ u 2 ⎦<br />

2<br />

θ 2 3<br />

χ 2 = − d2 u 3<br />

dx 2 1<br />

[ 6ξ<br />

= −<br />

L<br />

( ) 2<br />

= − d2 u 3 dξ<br />

dξ 2 = − 4 [<br />

φ<br />

1<br />

1<br />

dx 1 L 2 ′ ξξ u1 3 − ϕ1′ ξξ θ1 2 + φ2′ ξξ u2 3 − ]<br />

ϕ2′ ξξ θ2 2<br />

⎡ ⎤<br />

+ 3ξ<br />

, −−1 , − 6ξ ]<br />

u 1 3<br />

+ 3ξ<br />

, −1 ⎢ θ 1 2 ⎥<br />

2<br />

L L2 L ⎣ u 2 ⎦<br />

3<br />

θ 2 2<br />

⎡<br />

⎢<br />

⎣<br />

Que pue<strong>de</strong> resumirse en la siguiente expresión matricial<br />

⎤<br />

ε<br />

χ 1<br />

χ 3<br />

χ 2<br />

⎥<br />

⎦ = 1 L<br />

⎡<br />

⎢<br />

⎣<br />

−1 +1<br />

−1 +1<br />

6ξ<br />

−1 + 3ξ − 6ξ<br />

L L<br />

−1 + 3ξ<br />

− 6ξ<br />

L<br />

6ξ<br />

L<br />

+1 + 3ξ<br />

+1 + 3ξ<br />

⎡<br />

⎤<br />

⎥<br />

⎦<br />

⎢<br />

⎣<br />

u 1 1<br />

u 1 2<br />

u 1 3<br />

θ 1 1<br />

θ 1 2<br />

θ 1 3<br />

u 2 1<br />

u 2 2<br />

u 2 3<br />

θ 2 1<br />

θ 2 2<br />

θ 2 3<br />

⎤<br />

⎥<br />

⎦<br />

70

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

Saved successfully!

Ooh no, something went wrong!