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.

6.7.3.2. Fuerzas <strong>de</strong> contorno<br />

El trabajo virtual <strong>de</strong> las fuerzas <strong>de</strong> contorno, sólo se consi<strong>de</strong>ra sobre aquellos elementos que<br />

tengan un lado o cara coinci<strong>de</strong>nte con el contorno don<strong>de</strong> se conocen las fuerzas exteriores. Sea<br />

entonces un elemento cualquiera sobre el contorno <strong>de</strong>l cuerpo que tiene cargas actuantes sobre una<br />

<strong>de</strong> sus caras. Dicha cara estará <strong>de</strong>finida por un subconjunto <strong>de</strong> los nudos <strong>de</strong>l elemento NC. Notar<br />

que <strong>de</strong>bido a las exigencias impuestas sobre las funciones <strong>de</strong> interpolación, los <strong>de</strong>splazamientos<br />

sobre una cara <strong>de</strong>l elemento <strong>de</strong>pen<strong>de</strong>n sólo <strong>de</strong> los nudos sobre la cara es <strong>de</strong>cir sobre el subconjunto<br />

NC. Los <strong>de</strong>splazamientos virtuales pue<strong>de</strong>n entonces escribirse:<br />

∑NC<br />

δu =<br />

I=1<br />

φ I δu I<br />

De la misma forma es posible <strong>de</strong>scribir la carga externa<br />

∑NC<br />

f = φ I f I<br />

⎡<br />

⎤<br />

φ 1 φ 2 φ NC<br />

f = ⎣ φ 1 φ 2 ... ... ... φ NC ⎦<br />

φ 1 φ 2 φ NC } {{ }<br />

⎢<br />

⎣<br />

I=1<br />

¯Φ(ξ)<br />

⎡<br />

f 1 1<br />

f 1 2<br />

f 1 3<br />

...<br />

...<br />

f NC<br />

1<br />

f2<br />

NC<br />

f NC<br />

3<br />

⎤<br />

⎥<br />

⎦<br />

} {{ }<br />

f e<br />

Reemplazando en la expresión <strong>de</strong>l trabajo virtual <strong>de</strong> las fuerzas <strong>de</strong> contacto, tenemos que:<br />

∫<br />

∫<br />

S e<br />

f δu dS e = δu T e<br />

S e<br />

˜Φ T ˜Φ dS f e = δu T e<br />

˜M f e = δu T e<br />

g e<br />

Similarmente al caso <strong>de</strong> la matriz <strong>de</strong> rigi<strong>de</strong>z, se obtiene un vector global <strong>de</strong> cargas nodales<br />

equivalentes r ensamblando las contribuciones <strong>de</strong> las fuerzas másicas elementales y <strong>de</strong> las fuerzas<br />

<strong>de</strong> contacto. El trabajo virtual externo pue<strong>de</strong> escribirse en forma compacta como:<br />

∫<br />

∫<br />

− F δu dv − f δu dS σ = −δu T G r<br />

v<br />

S σ<br />

Sumando entonces trabajo virtual interno mas externo e igualando a 0.<br />

∫<br />

∫<br />

∫<br />

σ ij δε ij dv − F δu dv − f δu dS σ<br />

∼ = δu<br />

T<br />

G K u G − δu T G r = 0<br />

v<br />

v<br />

S σ<br />

Finalmente la condición impuesta por el P.T.V. requiere que los δu G puedan tomar cualquier<br />

valor en forma in<strong>de</strong>pendiente lo que conduce al siguiente sistema <strong>de</strong> ecuaciones lineales algebraicas<br />

simultáneas<br />

K u G = r<br />

6.8. Elemento cuadrilátero <strong>de</strong> cuatro nodos<br />

Como ejemplo <strong>de</strong> un elemento sencillo obtendremos la matriz <strong>de</strong> rigi<strong>de</strong>z y el vector <strong>de</strong> cargas<br />

nodales equivalentes <strong>de</strong> un elemento cuadrilátero <strong>de</strong> 4 nodos para estados planos (problemas<br />

bidimensionales)<br />

120

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

Saved successfully!

Ooh no, something went wrong!