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 ...
Para resolver el problema se propone una función φ (función de tensión) que satisfaga en forma implícita las ecuaciones de equilibrio del problema en el dominio. Habíamos visto que en este caso la única ecuación de equilibrio no trivial es Definiendo φ tal que ∂σ xz ∂x + ∂σ yz ∂y = 0 τ xz = ∂φ ∂y τ yz = − ∂φ ∂x reemplazado en la ecuación de equilibrio del dominio la satisface en forma idéntica. A su vez si la reemplazamos en la ecuación de equilibrio en el contorno τ ν = 0 (5.4) resulta τ ν = ∂φ dy ∂y ds + ∂φ ∂x dx ds = dφ ds = 0 La expresión anterior indica que φ (s) =cte. en el contorno. Para dominios simplemente conexos basta con fijar un valor constante arbitrario para la función incógnita sobre todo el contorno, valor que se elige igual 0 por comodidad. Si escribimos ahora las deformaciones en términos de la función de tensión, tenemos γ xz = τ xz G xz = 1 G xz ∂φ ∂y = −βy + w′ x γ yz = τ yz G yz = − 1 G yz ∂φ ∂x = βx + w′ y Derivando la primera respecto a y , la segunda respecto a x ( ) ∂ 1 ∂φ = −β + ∂y G xz ∂y w′ xy − ∂ ( ) 1 ∂φ = β + ∂x G yz ∂x w′ yx y observando que debe cumplirse que para que la función de alabeo w sea compatible w′ xy = w′ yx, restando la segunda de la primera resulta ( ) ∂ 1 ∂φ + ∂ ( ) 1 ∂φ = −2β (5.6) ∂x G yz ∂x ∂y G xz ∂y Que es una ecuación de compatibilidad (w′ xy = w′ yx) en función de φ. Si el material es homogéneo e isótropo resulta la ecuación de Laplace con condiciones de contorno homogéneas (dominios simplemente conexos) ∇ · ∇φ = −2Gβ φ (s) = 0 5.6.7. Forma débil de la ecuación de compatibilidad Para obtener la forma débil multiplicamos la ecuación 5.6 por una función de peso ψ e integramos por partes el primer miembro ∫ S ∫ A ψ [ ∂ ∂x ( 1 ∂φ G yz ∂x ) + ∂ ( )] ∫ 1 ∂φ dA = −ψ2β dA (5.7) ∂y G xz ∂y A [( ) ( ) ] ∫ [( ) ( )] ∫ 1 ∂φ 1 ∂φ 1 ∂φ 1 ∂φ ψ ν x + ν y ds − ∇ψ · , dA = −2β ψ dA G yz ∂x G xz ∂y A G xz ∂y G yz ∂x A (5.8) 93
Como φ es conocida sobre el contorno (0) la función de peso ψ vale 0 sobre el contorno y la integral sobre el contorno se anula. En consecuencia la forma débil resulta ∫ A [ ∂ψ ∂x , ∂ψ ] ⎡ 1 ⎢ G 0 ⎣ xz ∂y 1 0 G yz ⎤ ⎡ ⎥ ⎢ ⎦ ⎣ ∂φ ∂x ∂φ ∂y ⎤ ⎥ ⎦ dA = −2β ∫ A ψ dA (5.9) En el caso de dominios multiplemente conexos debe cumplirse que la función de tensión sea constante en cada contorno. Fijando el valor φ = 0, en el contorno exterior, en cada contorno interno S i la función valdrá ¯φ i . Estos valores de ¯φ i son ‘a priori’ desconocidos y se obtienen de la solución numérica. Lo que debe imponerse es que en todos los puntos j de cada contorno interno i el valor de φ sea el mismo φ j = ¯φ i (∀j ∈ S i ) 5.7. Flujo potencial En el caso de flujos potenciales las condiciones que debe cumplir el campo de velocidades son dos 1. continuidad, asociado con que la divergencia del campo de velocidades u sea nula ∇ · u = 0 2. irrotacionalidad, asociado a que el rotor del campo de velocidades sea nulo ∇ × u = 0 Hay dos formas equivalentes de abordar el problema, con diferentes ventajas de acuerdo al problema que se intenta analizar. La variable fundamental del problema no es el campo de velocidades, sino que éste (como todos los flujos tratados hasta ahora) se derivan de una variable (escalar). 5.7.1. Función potencial Una primera posibilidad es utilizar como variable independiente al potencial ϕ, de esta manera, el campo de velocidades resulta u = −∇ϕ Naturalmente si el campo u se deriva de un potencial, entonces satisface en forma explícita la condición irrotacionalidad y lo único que resta imponer es que satisfaga continuidad, es decir ∇ · u = −∇ · ∇ϕ = 0 con lo cual se obtiene la ecuación de Laplace. Las condiciones de contorno que pueden imponerse en este caso son: 1. esenciales, es posible fijar el valor de ϕ (potencial hidráulico) 2. naturales, es posible fijar el valor de la velocidad normal al contorno u ν = ν · ∇ϕ 5.7.2. Función líneas de corriente La segunda posibilidad es utilizar como variable independiente la función línea de corriente ψ que es conjugada de la función potencial ϕ. El campo de velocidades u queda ahora definido por ⎡ ⎤ 94 u = [ u1 ] = u 2 ⎢ ⎣ ∂ψ ∂x 2 − ∂ψ ∂x 1 ⎥ ⎦
- Page 48 and 49: N 1 (ξ) = (ξ − ξ 2) (ξ 1 −
- Page 50 and 51: modo de ubicar un nodo en todo punt
- Page 52 and 53: Figura 5 Malla de N nudos y N-1 ele
- Page 54 and 55: 3. Condiciones naturales de Neumann
- Page 56 and 57: (c) Considere las condiciones (ii)
- Page 58 and 59: los polinomios de interpolación in
- Page 60 and 61: La derivada segunda de u respecto a
- Page 62 and 63: 4.3.3. Formulación a partir del Pr
- Page 64 and 65: ¯K = ∫ L 0 BT D B d¯x 1 = ¯K =
- Page 66 and 67: Las cargas normales al eje de la vi
- Page 68 and 69: 2. θ 1 (ξ) = ∑ 2 I=1 N I (ξ)
- Page 70 and 71: 4.5.7. Cambio de base La expresión
- Page 72 and 73: Figura 5 Viga con deformación de c
- Page 74 and 75: 4.6.1. Matriz de rigidez de una vig
- Page 76 and 77: Donde las funciones de forma son: N
- Page 78 and 79: (una para cada W J ) en función de
- Page 80 and 81: Ejercicio 1-Sea el problema de conv
- Page 82 and 83: o directamente las coordenadas noda
- Page 84 and 85: donde r es el residuo que se quiere
- Page 86: ⎡ K 3−4 = ⎢ ⎣ K 4−5 = ⎡
- Page 89 and 90: Figura 1 Conducción del calor en 2
- Page 91 and 92: 5.3. Forma variacional del problema
- Page 93 and 94: 5.5. Flujo en un medio poroso El fl
- Page 95 and 96: ortotropía coinciden con las direc
- Page 97: 5.6.5. Forma débil de la ecuación
- Page 101 and 102: La formulación débil que resulta
- Page 103 and 104: 5.9.2. Formulación débil usando r
- Page 105 and 106: ⎡ σ = ⎢ ⎣ ⎤ σ 11 σ 22 σ
- Page 107 and 108: Figura 7 Teoría de placas clásica
- Page 109 and 110: Los esfuerzos de corte transversal
- Page 111 and 112: Capítulo 6 Elementos finitos en do
- Page 113 and 114: el punto p (ξ, η), y calculamos s
- Page 115 and 116: Para el caso cuadrático (ξ, η) =
- Page 117 and 118: ∂u ∂x 2 = NN∑ I=1 [ ∂N I
- Page 119 and 120: Figura 5 Mapeamientos para funcione
- Page 121 and 122: Elemento (L ( 1 , L 2 , L 3 ) w l L
- Page 123 and 124: ahora dependiente sólo del valor d
- Page 125 and 126: ∫ v ∫ σ ij δε ij dv = v δε
- Page 127 and 128: 6.8.1. Funciones de interpolación,
- Page 129 and 130: Siendo todas las variables constant
- Page 131 and 132: (el I por ejemplo) con dicha fuente
- Page 133 and 134: Capítulo 7 Aspectos generales asoc
- Page 135 and 136: 7.2. Imposición de restricciones n
- Page 137 and 138: entonces ¯k ij = k ij + 1 ǫ a ia
- Page 139 and 140: Figura 2 viga flexible entre column
- Page 141 and 142: D Matriz diagonal (definida positiv
- Page 143 and 144: 7.5. Suavizado de Variables En el m
- Page 145 and 146: 138
- Page 147 and 148: Difícilmente se pueda encontrar un
Para resolver el problema se propone una función φ (función <strong>de</strong> tensión) que satisfaga en forma<br />
implícita las ecuaciones <strong>de</strong> equilibrio <strong>de</strong>l problema en el dominio. Habíamos visto que en este caso<br />
la única ecuación <strong>de</strong> equilibrio no trivial es<br />
Definiendo φ tal que<br />
∂σ xz<br />
∂x + ∂σ yz<br />
∂y = 0<br />
τ xz = ∂φ<br />
∂y<br />
τ yz = − ∂φ<br />
∂x<br />
reemplazado en la ecuación <strong>de</strong> equilibrio <strong>de</strong>l dominio la satisface en forma idéntica. A su vez si la<br />
reemplazamos en la ecuación <strong>de</strong> equilibrio en el contorno τ ν = 0 (5.4) resulta<br />
τ ν = ∂φ dy<br />
∂y<br />
ds + ∂φ<br />
∂x<br />
dx<br />
ds = dφ<br />
ds = 0<br />
La expresión anterior indica que φ (s) =cte. en el contorno. Para dominios simplemente conexos<br />
basta con fijar un valor constante arbitrario para la función incógnita sobre todo el contorno, valor<br />
que se elige igual 0 por comodidad.<br />
Si escribimos ahora las <strong>de</strong>formaciones en términos <strong>de</strong> la función <strong>de</strong> tensión, tenemos<br />
γ xz = τ xz<br />
G xz<br />
= 1<br />
G xz<br />
∂φ<br />
∂y = −βy + w′ x<br />
γ yz = τ yz<br />
G yz<br />
= − 1<br />
G yz<br />
∂φ<br />
∂x = βx + w′ y<br />
Derivando la primera respecto a y , la segunda respecto a x<br />
( )<br />
∂ 1 ∂φ<br />
= −β +<br />
∂y G xz ∂y<br />
w′ xy<br />
− ∂ ( ) 1 ∂φ<br />
= β +<br />
∂x G yz ∂x<br />
w′ yx<br />
y observando que <strong>de</strong>be cumplirse que para que la función <strong>de</strong> alabeo w sea compatible w′ xy = w′ yx,<br />
restando la segunda <strong>de</strong> la primera resulta<br />
( )<br />
∂ 1 ∂φ<br />
+ ∂ ( ) 1 ∂φ<br />
= −2β (5.6)<br />
∂x G yz ∂x ∂y G xz ∂y<br />
Que es una ecuación <strong>de</strong> compatibilidad (w′ xy = w′ yx) en función <strong>de</strong> φ. Si el material es homogéneo<br />
e isótropo resulta la ecuación <strong>de</strong> Laplace con condiciones <strong>de</strong> contorno homogéneas (dominios<br />
simplemente conexos)<br />
∇ · ∇φ = −2Gβ<br />
φ (s) = 0<br />
5.6.7. Forma débil <strong>de</strong> la ecuación <strong>de</strong> compatibilidad<br />
Para obtener la forma débil multiplicamos la ecuación 5.6 por una función <strong>de</strong> peso ψ e integramos<br />
por partes el primer miembro<br />
∫<br />
S<br />
∫<br />
A<br />
ψ<br />
[ ∂<br />
∂x<br />
( 1 ∂φ<br />
G yz ∂x<br />
)<br />
+ ∂ ( )] ∫<br />
1 ∂φ<br />
dA = −ψ2β dA (5.7)<br />
∂y G xz ∂y<br />
A<br />
[( ) ( ) ] ∫ [( ) ( )] ∫<br />
1 ∂φ 1 ∂φ<br />
1 ∂φ 1 ∂φ<br />
ψ<br />
ν x + ν y ds − ∇ψ ·<br />
,<br />
dA = −2β ψ dA<br />
G yz ∂x G xz ∂y<br />
A G xz ∂y G yz ∂x<br />
A<br />
(5.8)<br />
93