universidad nacional de c´ordoba - Facultad de Ciencias Exactas ...
universidad nacional de c´ordoba - Facultad de Ciencias Exactas ...
universidad nacional de c´ordoba - Facultad de Ciencias Exactas ...
Create successful ePaper yourself
Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.
aparentes <strong>de</strong> resonancia que no son reales. Otra dificultad <strong>de</strong> este tipo <strong>de</strong> métodos es que<br />
requieren un gran esfuerzo <strong>de</strong> mo<strong>de</strong>lación, ya que requieren la discretización <strong>de</strong>l volumen<br />
<strong>de</strong> continuo que ro<strong>de</strong>a a la estructura hasta el bor<strong>de</strong> ficticio. A pesar <strong>de</strong> las mencionadas<br />
dificulta<strong>de</strong>s, los métodos numéricos <strong>de</strong> dominio permiten la solución <strong>de</strong> una gran variedad<br />
<strong>de</strong> problemas con geometrías arbitrarias.<br />
Otra forma <strong>de</strong> solucionar problemas <strong>de</strong> interacción dinámica continuo-estructura es<br />
mediante métodos analíticos (por ejemplo [37]). Sin embargo, éstos son solamente viables<br />
para problemas <strong>de</strong> valores en el contorno <strong>de</strong> geometría sencilla. Los métodos <strong>de</strong> elementos<br />
<strong>de</strong> contorno combinan las ventajas <strong>de</strong> los métodos analíticos con las <strong>de</strong> los métodos<br />
numéricos <strong>de</strong> volumen, es <strong>de</strong>cir que son métodos numéricos aplicables a geometrías arbitrarias<br />
y que no requieren <strong>de</strong> la mo<strong>de</strong>lación <strong>de</strong>l continuo, sino <strong>de</strong> las superficies en don<strong>de</strong><br />
existen fuerzas <strong>de</strong> interacción.<br />
3.2 Formulación general <strong>de</strong>l IBEM<br />
El comportamiento elastodinámico en régimen estacionario <strong>de</strong> un medio elástico continuo<br />
(sin consi<strong>de</strong>rar fuerzas <strong>de</strong> volumen) pue<strong>de</strong> ser <strong>de</strong>scripto mediante la ecuación <strong>de</strong> Helmholtz:<br />
(λ + µ) ∂ ∇ · u + µ∇ 2 u + ω 2 ρ u = 0 (3.1)<br />
∂x i<br />
[ ] T<br />
don<strong>de</strong> λ y µ son las constantes elásticas <strong>de</strong> Lamé y u = u x u y u z<br />
es el vector <strong>de</strong>splazamiento.<br />
El planteo <strong>de</strong>l problema se completa con la introducción <strong>de</strong> las condiciones <strong>de</strong><br />
bor<strong>de</strong> naturales y esenciales. La solución a la ecuación 3.1 pue<strong>de</strong> encontrarse en términos<br />
<strong>de</strong> fuentes distribuidas, φ (X), <strong>de</strong> la siguiente manera (por ejemplo, [1], [5], [32]):<br />
∫<br />
u (x,ω) = G (x, X, ω) · φ (X,ω) dΩ (X) (3.2)<br />
don<strong>de</strong> G (x, X, ω) es el tensor <strong>de</strong> Green cuyas componentes están <strong>de</strong>finidas como los <strong>de</strong>splazamientos<br />
resultantes a la acción <strong>de</strong> una fuente cuya resultante total es unitaria (esto<br />
56