APENDICE

Partiendo de la ecuación de Navier-Stokes

𝜕𝑣⃗
𝜌( + (𝑣⃗ + 𝛻⃗⃗)𝑣⃗) = −𝛻⃗⃗𝑃 + 𝜌𝐾
⃗⃗ + (𝜆 + 𝜂)𝛻⃗⃗(𝛻⃗⃗ ⋅ 𝑣⃗) + 𝜌𝛾𝛻 2 𝑣⃗
𝜕𝑡

Flujo incompresible → ⃗∇⃗ ⋅ v

⃗⃗ = 0, ⃗⃗
k = gk̂ ya que es en 2D ρk
⃗⃗ se desprecia.

𝜕𝑣⃗
Flujo estacionario → =0
𝜕𝑡
𝜌(𝑣⃗ ⋅ 𝛻⃗⃗)𝑣⃗ = −𝛻⃗⃗𝑃 + 𝜌𝛾𝛻 2 𝑣⃗ ; donde 𝑣⃗ = 𝑢𝑖 + 𝑣𝑗

𝜕𝑢 𝜕𝑣 𝜕𝑢 𝜕𝑣
𝜕𝑥 𝜕𝑥 𝜕𝑃 𝜕𝑃 𝜕 𝜕 𝜕𝑥 𝜕𝑥
𝜌(𝑢 𝑣) ⋅ (𝜕𝑢 𝜕𝑣 )=(− 𝜕𝑥 − )+ 𝜌𝛾 (𝜕𝑥 𝜕𝑦) (𝜕𝑢 𝜕𝑣 )
𝜕𝑦
𝜕𝑦 𝜕𝑦 𝜕𝑦 𝜕𝑦

Así se obtienen las siguientes ecuaciones:

𝜕𝑢 𝜕𝑢 𝜕 𝜕𝑢 𝜕 𝜕𝑢 𝜕𝑃
𝜌𝑢 + 𝜌𝑣 − 𝜌𝛾 ( ) − 𝜌𝛾 ( ) + =0
𝜕𝑥 𝜕𝑦 𝜕𝑥 𝜕𝑥 𝜕𝑦 𝜕𝑦 𝜕𝑥
𝜕𝑣 𝜕𝑣 𝜕 𝜕𝑣 𝜕 𝜕𝑣 𝜕𝑃
𝜌𝑢 + 𝜌𝑣 − 𝜌𝛾 ( ) − 𝜌𝛾 ( ) + =0
𝜕𝑥 𝜕𝑦 𝜕𝑥 𝜕𝑥 𝜕𝑦 𝜕𝑦 𝜕𝑦
𝜕𝑢 𝜕𝑣
+ = 0 "𝑃𝑜𝑟 𝑠𝑒𝑟 𝑖𝑛𝑐𝑜𝑚𝑝𝑟𝑒𝑠𝑖𝑏𝑙𝑒"
𝜕𝑥 𝜕𝑦

Ahora multiplicando cada ecuación por una función de peso 𝑢̅, 𝑣̅ 𝑦 𝑝̅ respectivamente, e
integrando en el dominio.
𝜕𝑢 𝜕𝑢 𝜕 𝜕𝑢 𝜕 𝜕𝑢 𝜕𝑃
∬{𝜌𝑢 + 𝜌𝑣 − 𝜌𝛾 ( ) − 𝜌𝛾 ( ) + }𝑢̅𝑑𝑥𝑑𝑦 = 0
𝜕𝑥 𝜕𝑦 𝜕𝑥 𝜕𝑥 𝜕𝑦 𝜕𝑦 𝜕𝑥
𝑠

𝜕𝑣 𝜕𝑣 𝜕 𝜕𝑣 𝜕 𝜕𝑣 𝜕𝑃
∬{𝜌𝑢 + 𝜌𝑣 − 𝜌𝛾 ( ) − 𝜌𝛾 ( ) + }𝑣̅ 𝑑𝑥𝑑𝑦 = 0
𝜕𝑥 𝜕𝑦 𝜕𝑥 𝜕𝑥 𝜕𝑦 𝜕𝑦 𝜕𝑦
𝑠

𝜕𝑢 𝜕𝑣
∬{ + }𝑝̅ 𝑑𝑥𝑑𝑦 = 0
𝜕𝑥 𝜕𝑦
𝑠

Por otro lado, se tiene para u que:

𝜕 𝜕𝑢 𝜕𝑢 𝜕𝑢̅ 𝜕 𝜕𝑢 𝜕 𝜕𝑢 𝜕𝑢 𝜕𝑢̅ 𝜕 𝜕𝑢
(𝑢̅ ) = + 𝑢̅ ( ) → −𝑢̅ ( )= − (𝑢̅ )
𝜕𝑥 𝜕𝑥 𝜕𝑥 𝜕𝑥 𝜕𝑥 𝜕𝑥 𝜕𝑥 𝜕𝑥 𝜕𝑥 𝜕𝑥 𝜕𝑥 𝜕𝑥
𝜕 𝜕𝑢 𝜕𝑢 𝜕𝑢̅ 𝜕 𝜕𝑢 𝜕 𝜕𝑢 𝜕𝑢 𝜕𝑢̅ 𝜕 𝜕𝑢
(𝑢̅ ) = + 𝑢̅ ( ) → −𝑢̅ ( )= − (𝑢
⃗⃗ )
𝜕𝑦 𝜕𝑦 𝜕𝑦 𝜕𝑦 𝜕𝑦 𝜕𝑦 𝜕𝑦 𝜕𝑦 𝜕𝑦 𝜕𝑦 𝜕𝑦 𝜕𝑦
Lo mismo para v:
𝜕 𝜕𝑣 𝜕𝑣 𝜕𝑣̅ 𝜕 𝜕𝑣 𝜕 𝜕𝑣 𝜕𝑣 𝜕𝑣̅ 𝜕 𝜕𝑣
(𝑣̅ ) = + 𝑣̅ ( ) → −𝑣̅ ( )= − (𝑣̅ )
𝜕𝑥 𝜕𝑥 𝜕𝑥 𝜕𝑥 𝜕𝑥 𝜕𝑥 𝜕𝑥 𝜕𝑥 𝜕𝑥 𝜕𝑥 𝜕𝑥 𝜕𝑥
𝜕 𝜕𝑣 𝜕𝑣 𝜕𝑣̅ 𝜕 𝜕𝑣 𝜕 𝜕𝑣 𝜕𝑣 𝜕𝑣̅ 𝜕 𝜕𝑣
(𝑣̅ ) = + 𝑣̅ ( ) → −𝑣̅ ( )= − (𝑣̅ )
𝜕𝑦 𝜕𝑦 𝜕𝑦 𝜕𝑦 𝜕𝑦 𝜕𝑦 𝜕𝑦 𝜕𝑦 𝜕𝑦 𝜕𝑦 𝜕𝑦 𝜕𝑦
Entonces se tiene para la primera integral:
𝜕𝑢 𝜕𝑢 𝜕𝑃 𝜕𝑢 𝜕𝑢̅ 𝜕𝑢 𝜕𝑢̅ 𝜕 𝜕𝑢 𝜕 𝜕𝑢
∬( 𝜌𝑢 + 𝜌𝑣 + )𝑢̅𝑑𝑥𝑑𝑦 + ∬ 𝜌 𝛾( + )𝑑𝑥𝑑𝑦 − ∬ 𝜌 𝛾{ (𝑢̅ ) + (𝑢̅ )}𝑑𝑥𝑑𝑦 = 0
𝜕𝑥 𝜕𝑦 𝜕𝑥 𝜕𝑥 𝜕𝑥 𝜕𝑦 𝜕𝑦 𝜕𝑥 𝜕𝑥 𝜕𝑦 𝜕𝑦
𝑠 𝑆 𝑆

Por el teorema de la divergencia:

𝜕 𝜕𝑢 𝜕 𝜕𝑢
∬{ (𝑢̅ ) + ⃗⃗) ⋅ 𝑑𝑙⃗
(𝑢̅ )}𝑑𝑥𝑑𝑦 = ∫(𝑢̅𝛻⃗⃗𝑢
𝜕𝑥 𝜕𝑥 𝜕𝑦 𝜕𝑦
𝑠 𝐶

Puesto que se realizará sobre una trayectoria cerrada la última integral de este campo vectorial
sobre la curva es cero y lo mismo aplica para la segunda integral, entonces se tiene:
̅
𝜕𝑢 𝜕𝑢 ̅
𝜕𝑢 𝜕𝑢 𝜕𝑢 𝜕𝑢 𝜕𝑃
∬𝑠 𝜌𝛾{𝜕𝑥 𝜕𝑥 + 𝜕𝑦 𝜕𝑦 }𝑑𝑥𝑑𝑦 + ∬𝑠{𝜌𝑢 𝜕𝑥 + 𝜌𝑣 𝜕𝑦 + 𝜕𝑥 }𝑢̅𝑑𝑥𝑑𝑦 = 0…………………….1

𝜕𝑣 𝜕𝑣̅ 𝜕𝑣 𝜕𝑣̅ 𝜕𝑣 𝜕𝑣 𝜕𝑃
∬𝑠 𝜌𝛾{𝜕𝑥 𝜕𝑥 + 𝜕𝑦 𝜕𝑦}𝑑𝑥𝑑𝑦 + ∬𝑠{ 𝜌𝑢 𝜕𝑥 + 𝜌𝑣 𝜕𝑦 + 𝜕𝑥 }𝑣̅ 𝑑𝑥𝑑𝑦 = 0……………………….2

𝜕𝑢 𝜕𝑣
∬𝑠 (𝜕𝑥 + 𝜕𝑦) 𝑝̅ 𝑑𝑥 𝑑𝑦 = 0…………………………………………………………………………………….3

Ahora analizando el elemento en un sistema (𝜂, 𝜉) ∶

(𝑣𝑥 , 𝑣𝑦 ) → (𝑣𝜂 , 𝑣𝜉 )
Se propone: Ni= ai + bi𝜂 + 𝑐𝑖 𝜉
Función de forma que debe cumplir: N1(0,0)= 0 N1(1,0)= 0 N1(0,1)= 1
N2(0,0)= 0 N2(1,0)= 1 N2(0,1)= 1
N3(0,0)= 1 N3(1,0)= 0 N3(0,1)= 1
Se tienen las siguientes funciones de forma:
N1= 𝜉 , N2 = 𝜂 , N3= 1-𝜂 -𝜉

⃗⃗ = 𝑁𝑗 , ⃗⃗⃗⃗⃗
Ahora se tiene que : 𝑢 𝑉 = 𝑁𝑗 𝑦 𝑃⃗⃗ = 𝑁𝑗 también 𝑢 = ∑𝑖 𝑁𝑖 𝑢𝑖 , 𝑣 = ∑𝑖 𝑁𝑖 𝑣𝑖 y 𝑃 = ∑𝑖 𝑁𝑖 𝑃𝑖

Es decir: u=𝜉u1 + 𝜂u2 + (1- 𝜉- 𝜂)u3 , v= 𝜉v1 + 𝜂v2 + (1- 𝜉- 𝜂)v3 y p= 𝜉P1 + 𝜂P2 + (1- 𝜉- 𝜂)P3

Entonces las parciales con respecto a 𝜉 y 𝜂 son:

𝑢1 𝑣1 𝑝1
𝜕𝑢 𝜕𝑣 𝜕𝑝
𝜕𝜉
= (1,0, −1) (𝑢2 ) ,
𝜕𝜉
= (1,0, −1) (𝑣2 ) ,
𝜕𝜉
= (1,0, −1) (𝑝2 )
𝑢3 𝑣3 𝑝3
𝑢1 𝑣1 𝑝1
𝜕𝑢 𝜕𝑣 𝜕𝑝
𝜕𝜂
= (0,1, −1) (𝑢2 ) ,
𝜕𝜂
= (0,1, −1) (𝑣2 ) ,
𝜕𝜂
= (0,1, −1) (𝑝2 )
𝑢3 𝑣3 𝑝3
Por otro lado, se tiene que:

𝑥 = 𝜉(𝑥1 − 𝑥3 ) + 𝜂(𝑥2 − 𝑥3 ) + 𝑥3
𝑦 = 𝜉(𝑦1 − 𝑦3 ) + 𝜂(𝑦2 − 𝑦3 ) + 𝑦3
𝑥 − 𝑥3 𝑥13 𝑥23 𝜉
(𝑦 − 𝑦3 ) = ( 𝑦13 𝑦23 ) (𝜂 )
Donde 𝑥𝑖𝑗 = 𝑥𝑖 − 𝑥𝑗 así también 𝑦𝑖𝑗 = 𝑦𝑖 − 𝑦𝑗

𝑥13 𝑥23 −1 𝑥 − 𝑥3 𝜉
(𝑦 𝑦23 ) (𝑦 − 𝑦3 ) =(𝜂)
13

𝑦23 (𝑥−𝑥3 ) 𝑥23 (𝑦−𝑦3 )

Así : 𝜉 = (𝑥 -
13 𝑦23 −𝑥23 𝑦13 ) (𝑥13 𝑦23 −𝑥23 𝑦13 )

𝑥23 (𝑦−𝑦3 ) 𝑦13 (𝑥−𝑥3 )

y 𝜂 = (𝑥 -
13 𝑦23 −𝑥23 𝑦13 ) (𝑥13 𝑦23 −𝑥23 𝑦13 )

Si hacemos a=(𝑥13 𝑦23 − 𝑥23 𝑦13 )

𝜕𝜉 𝑦23 𝜕𝜂 𝑦13 𝜕𝜉 𝑥23 𝜕𝜂 𝑥13
Entonces las parciales son: = , =− , =− y =
𝜕𝑥 𝑎 𝜕𝑥 𝑎 𝜕𝑦 𝑎 𝜕𝑥 𝑎

Ahora para la ecuación de las velocidades en x es decir para u y 𝑢 ⃗⃗, usando las funciones de
forma 𝑁1 , 𝑁2 𝑦 𝑁3 para 𝑢
⃗⃗ , se tienen las siguientes integrales de la ec. 1:
 𝜕𝑢 𝜕𝑁1 𝜕𝑢 𝜕𝑁1 𝜕𝑢 𝜕𝑢 𝜕𝑃
∬𝑠 𝜌𝛾 (𝜕𝑥 𝜕𝑥
+ 𝜕𝑦 𝜕𝑦
) 𝑑𝑥 𝑑𝑦 + ∬𝑠 (𝜌𝑢 𝜕𝑥 + 𝜌𝑣 𝜕𝑦 + 𝜕𝑥 ) 𝑁1 𝑑𝑥 𝑑𝑦 = 0………….1.1

𝜕𝑢 𝜕𝑁2 𝜕𝑢 𝜕𝑁2 𝜕𝑢 𝜕𝑢 𝜕𝑃
∬𝑠 𝜌𝛾 (𝜕𝑥 𝜕𝑥
+ 𝜕𝑦 𝜕𝑦
) 𝑑𝑥 𝑑𝑦 + ∬𝑠 (𝜌𝑢 𝜕𝑥 + 𝜌𝑣 𝜕𝑦 + 𝜕𝑥 ) 𝑁2 𝑑𝑥 𝑑𝑦 = 0………….1.2

𝜕𝑢 𝜕𝑁3 𝜕𝑢 𝜕𝑁3 𝜕𝑢 𝜕𝑢 𝜕𝑃
∬𝑠 𝜌𝛾 (𝜕𝑥 𝜕𝑥
+ 𝜕𝑦 𝜕𝑦
) 𝑑𝑥 𝑑𝑦 + ∬𝑠 (𝜌𝑢 𝜕𝑥 + 𝜌𝑣 𝜕𝑦 + 𝜕𝑥 ) 𝑁3 𝑑𝑥 𝑑𝑦 = 0………….1.3

Ahora calculando el producto de parciales:

𝑢1 𝑢1 𝑢1
𝜕𝑢 𝜕𝑁1
=(𝜕𝑢 𝜕𝜉
+
𝜕𝑢 𝜕𝜂 𝜕𝜉 𝑦 𝑦 𝑦 𝑦
) ( ) = [(0,1, −1) (𝑢2 ) 23 + (0,1, −1) (𝑢2 ) (− 13)] 23 = 232 [(𝑦23 , −𝑦13 , 𝑦12 ) (𝑢2 )]
𝜕𝑥 𝜕𝑥 𝜕𝜉 𝜕𝑥 𝜕𝜂 𝜕𝑥 𝜕𝑥
𝑢3
𝑎
𝑢3
𝑎 𝑎 𝑎
𝑢3

Procediendo de forma semejante para los demás productos de parciales se obtiene:

𝑢1
𝜕𝑢 𝜕𝑁1 𝑥23
= - 𝑎2 𝑢
[(−𝑥23 , 𝑥13 , 𝑥21 ) ( 2 )]
𝜕𝑦 𝜕𝑦
𝑢3
𝑢1
𝜕𝑢 𝜕𝑁2 𝑦13
= - 𝑎2 𝑢
[(𝑦23 , −𝑦13 , 𝑦12 ) ( 2 )]
𝜕𝑥 𝜕𝑥
𝑢3
𝑢1
𝜕𝑢 𝜕𝑁2 𝑥13
= [(−𝑥23 , 𝑥13 , 𝑥21 ) (𝑢2 )]
𝜕𝑦 𝜕𝑦 𝑎2
𝑢3
𝑢1
𝜕𝑢 𝜕𝑁3 𝑦12
= [(𝑦23 , −𝑦13 , 𝑦12 ) (𝑢2 )]
𝜕𝑥 𝜕𝑥 𝑎2
𝑢3
𝑢1
𝜕𝑢 𝜕𝑁3 𝑥21
= [(−𝑥23 , 𝑥13 , 𝑥21 ) (𝑢2 )]
𝜕𝑦 𝜕𝑦 𝑎2
𝑢3
Ahora para la segunda integral, se calculan las parciales por la velocidad:
𝑢1
𝜕𝑢 𝜕𝑢 𝜕𝜉 𝜕𝑢 𝜕𝜂
𝑢 = (𝜉, 𝜂, 1 − 𝜉 − 𝜂) (𝑢2 ) ( + )
𝜕𝑥 𝑢3 𝜕𝜉 𝜕𝑥 𝜕𝜂 𝜕𝑥
𝑢1 𝑢1
1
= (𝜉, 𝜂, 1 − 𝜉 − 𝜂) (𝑢2 ) [𝑎 (𝑦23 , −𝑦13 , 𝑦12 ) (𝑢2 )]
𝑢3 𝑢3
𝑦23 𝑢12 − 𝑦13 𝑢1 𝑢2 + 𝑦12 𝑢1 𝑢3
1
= 𝑎 (𝜉, 𝜂, 1 − 𝜉 − 𝜂) (𝑦23 𝑢1 𝑢2 − 𝑦13 𝑢22 + 𝑦12 𝑢2 𝑢3 )
𝑦23 𝑢1 𝑢3 − 𝑦13 𝑢2 𝑢3 − 𝑦12 𝑢32
 −𝑥23 𝑢1 𝑣1 + 𝑥13 𝑢2 𝑣1 + 𝑥21 𝑢3 𝑣1
𝜕𝑢 1
𝑣 𝜕𝑦 = 𝑎 (𝜉, 𝜂, 1 − 𝜉 − 𝜂) (−𝑥23 𝑢1 𝑣2 + 𝑥13 𝑢2 𝑣2 + 𝑥21 𝑢3 𝑣2 )
−𝑥23 𝑢1 𝑣3 + 𝑥13 𝑢2 𝑣3 + 𝑥21 𝑢3 𝑣3
𝑝1 𝑝1
𝜕𝑝 𝜕𝑝 𝜕𝜉 𝜕𝑝 𝜕𝜂 𝑦23 𝑦
= 𝜕𝜉 𝜕𝑥 + 𝜕𝜂 𝜕𝑥 =(1,0, −1) (𝑝2 ) 𝑎 + (0,1, −1) (𝑝2 ) (− 𝑎13)
𝜕𝑥
𝑝3 𝑝3
𝑝1
1
= 𝑎 (𝑦23 , −𝑦13 , 𝑦12 ) (𝑝2 )
𝑝3
𝑝1
𝜕𝑝 1
𝑝
= 𝑎 (−𝑥23 , 𝑥13 , 𝑥21 ) ( 2 )
𝜕𝑦
𝑝3
Ahora sustituyendo todos los valores calculados en las ecuaciones, para las tres
ecuaciones de u con sus respectivas funciones de forma (1.1,1.2,1.3):
𝑢1 𝑢1
𝑦23 𝑥23
∬𝑠{𝜌𝛾 ( 𝑎2 [(𝑦23 , −𝑦13 , 𝑦12 ) (𝑢2 )] − 𝑎2 [(−𝑥23 , 𝑥13 , 𝑥21 ) (𝑢2 )])
𝑢3 𝑢3
𝑦23 𝑢12 − 𝑦13 𝑢1 𝑢2 + 𝑦12 𝑢1 𝑢3
1
+𝜌 (𝜉, 𝜂, 1 − 𝜉 − 𝜂) (𝑦23 𝑢1 𝑢2 − 𝑦13 𝑢22 + 𝑦12 𝑢2 𝑢3 ) 𝜉
𝑎
𝑦23 𝑢1 𝑢3 − 𝑦13 𝑢2 𝑢3 − 𝑦12 𝑢32

1 −𝑥23 𝑢1 𝑣1 + 𝑥13 𝑢2 𝑣1 + 𝑥21 𝑢3 𝑣1

+𝜌 (𝜉, 𝜂, 1 − 𝜉 − 𝜂) (−𝑥23 𝑢1 𝑣2 + 𝑥13 𝑢2 𝑣2 + 𝑥21 𝑢3 𝑣2 ) 𝜉
𝑎 −𝑥 𝑢 𝑣 + 𝑥 𝑢 𝑣 + 𝑥 𝑢 𝑣
23 1 3 13 2 3 21 3 3
𝑝1
1
+ (𝑦23 , −𝑦13 , 𝑦12 ) (𝑝2 ) 𝜉}𝑑𝑥 𝑑𝑦 = 0 ……………………………………………..…1.4
𝑎
𝑝3
𝑢1 𝑢1
𝑦13 𝑥13
∬𝑠{𝜌𝛾 (- 𝑎2 [(𝑦23 , −𝑦13 , 𝑦12 ) (𝑢2 )] + 𝑎2 [(−𝑥23 , 𝑥13 , 𝑥21 ) (𝑢2 )])
𝑢3 𝑢3
𝑦23 𝑢12 − 𝑦13 𝑢1 𝑢2 + 𝑦12 𝑢1 𝑢3
1
+𝜌 (𝜉, 𝜂, 1 − 𝜉 − 𝜂) (𝑦23 𝑢1 𝑢2 − 𝑦13 𝑢22 + 𝑦12 𝑢2 𝑢3 ) 𝜂
𝑎
𝑦23 𝑢1 𝑢3 − 𝑦13 𝑢2 𝑢3 − 𝑦12 𝑢32

1 −𝑥23 𝑢1 𝑣1 + 𝑥13 𝑢2 𝑣1 + 𝑥21 𝑢3 𝑣1

+𝜌 (𝜉, 𝜂, 1 − 𝜉 − 𝜂) (−𝑥23 𝑢1 𝑣2 + 𝑥13 𝑢2 𝑣2 + 𝑥21 𝑢3 𝑣2 ) 𝜂
𝑎 −𝑥 𝑢 𝑣 + 𝑥 𝑢 𝑣 + 𝑥 𝑢 𝑣
23 1 3 13 2 3 21 3 3
 𝑝1
1
+ (𝑦23 , −𝑦13 , 𝑦12 ) (𝑝2 ) 𝜂}𝑑𝑥 𝑑𝑦 = 0 ……………………………………………..…1.5
𝑎
𝑝3
𝑢1 𝑢1
𝑦12 𝑥21
∬𝑠{𝜌𝛾 ( 𝑎2 [(𝑦23 , −𝑦13 , 𝑦12 ) (𝑢2 )] + 𝑎2 [(−𝑥23 , 𝑥13 , 𝑥21 ) (𝑢2 )])
𝑢3 𝑢3
𝑦23 𝑢12 − 𝑦13 𝑢1 𝑢2 + 𝑦12 𝑢1 𝑢3
1
+𝜌 (𝜉, 𝜂, 1 − 𝜉 − 𝜂) (𝑦23 𝑢1 𝑢2 − 𝑦13 𝑢22 + 𝑦12 𝑢2 𝑢3 ) (1 + 𝜉)
𝑎
𝑦23 𝑢1 𝑢3 − 𝑦13 𝑢2 𝑢3 − 𝑦12 𝑢32

1 −𝑥23 𝑢1 𝑣1 + 𝑥13 𝑢2 𝑣1 + 𝑥21 𝑢3 𝑣1

+𝜌 (𝜉, 𝜂, 1 − 𝜉 − 𝜂) (−𝑥23 𝑢1 𝑣2 + 𝑥13 𝑢2 𝑣2 + 𝑥21 𝑢3 𝑣2 ) (1 + 𝜉)
𝑎 −𝑥 𝑢 𝑣 + 𝑥 𝑢 𝑣 + 𝑥 𝑢 𝑣 23 1 3 13 2 3 21 3 3
𝑝1
1
+ (𝑦23 , −𝑦13 , 𝑦12 ) (𝑝2 ) (1 + 𝜉)}𝑑𝑥 𝑑𝑦 = 0 …………………………………..…1.6
𝑎
𝑝3
𝜕𝑥 𝜕𝑥
𝜕𝜉 𝜕𝜂 𝑥13 𝑥23
Calculando el Jacobiano: J= |𝜕y 𝜕y
|=| y13 y23 |=𝑥13 y23 -𝑥23 y13 =a
𝜕𝜉 𝜕𝜂

Ahora sustituyendo en las integrales anteriores e integrando sobre la segunda se tiene:

𝑢1 1
𝜌𝛾 1−𝜉
2
(𝑦23 + 𝑥23 , -y23 y13 -𝑥23 𝑥13 , y23 y12 -𝑥21 𝑥23 ) (𝑢2 )
2
∫ ∫0 𝑑𝜂 𝑑𝜉
a
𝑢3 0

1
1-𝜉
𝑦23 𝑢21 − 𝑦13 𝑢1 𝑢2 + 𝑦12 𝑢1 𝑢3
+𝜌 [∫ ∫ (𝜉 2 , 𝜉𝜂, 𝜉 − 𝜉 2 − 𝜉𝜂)𝑑𝜂 𝑑𝜉 ] 𝑦23 𝑢1 𝑢2 − 𝑦13 𝑢22 + 𝑦12 𝑢2 𝑢3
0 2
0
(𝑦23 𝑢1 𝑢3 − 𝑦13 𝑢2 𝑢3 − 𝑦12 𝑢3 )
(( ))
1
1-𝜉 −𝑥23 𝑢1 𝑣1 + 𝑥13 𝑢2 𝑣1 + 𝑥21 𝑢3 𝑣1
+𝜌 [∫ ∫ (𝜉 2 , 𝜉𝜂, 𝜉 − 𝜉 2 − 𝜉𝜂)𝑑𝜂 𝑑𝜉 ] (−𝑥23 𝑢1 𝑣2 + 𝑥13 𝑢2 𝑣2 + 𝑥21 𝑢3 𝑣2 )
0
0 −𝑥23 𝑢1 𝑣3 + 𝑥13 𝑢2 𝑣3 + 𝑥21 𝑢3 𝑣3
𝑝1
1 1-𝜉
(𝑦23 , −𝑦13 , 𝑦12 ) (𝑝2 ) ∫ ∫0 𝜉𝑑𝜂 𝑑𝜉 = 0………………………………………………………………1.7
𝑝3 0
 𝑢1 1
𝜌𝛾 1−𝜉
(-y23 y13 -𝑥23 𝑥13 , 𝑦13 + 𝑥13 , 𝑥21 𝑥23 -y13 y12 ) (𝑢2 ) ∫ ∫0 𝑑𝜂 𝑑𝜉
2 2
a
𝑢3 0
1
1-𝜉
𝑦23 𝑢21 − 𝑦13 𝑢1 𝑢2 + 𝑦12 𝑢1 𝑢3
+𝜌 [∫ ∫ (𝜉𝜂, 𝜂2 , 𝜂 − 𝜉𝜂 − 𝜂2 )𝑑𝜂 𝑑𝜉 ] 𝑦23 𝑢1 𝑢2 − 𝑦13 𝑢22 + 𝑦12 𝑢2 𝑢3
0 2
0
(𝑦23 𝑢1 𝑢3 − 𝑦13 𝑢2 𝑢3 − 𝑦12 𝑢3 )
1
1-𝜉 −𝑥23 𝑢1 𝑣1 + 𝑥13 𝑢2 𝑣1 + 𝑥21 𝑢3 𝑣1
+𝜌 [∫ ∫ 2
(𝜉𝜂, 𝜂 , 𝜂 − 𝜉𝜂 − 𝜂 2 )𝑑𝜂
𝑑𝜉 ] (−𝑥23 𝑢1 𝑣2 + 𝑥13 𝑢2 𝑣2 + 𝑥21 𝑢3 𝑣2 )
0
0 −𝑥23 𝑢1 𝑣3 + 𝑥13 𝑢2 𝑣3 + 𝑥21 𝑢3 𝑣3
𝑝1
1 1-𝜉
(𝑦23 , −𝑦13 , 𝑦12 ) (𝑝2 ) ∫ ∫0 𝜉𝑑𝜂 𝑑𝜉 = 0………………………………………………………………1.8
𝑝3 0
𝑢1 1
𝜌𝛾 1−𝜉
2 2 𝑢
(y12 y23 -𝑥21 𝑥23 , 𝑥21 𝑥13 -y13 y12 , 𝑦12 + 𝑥12 ) ( 2 ) ∫ ∫0 𝑑𝜂 𝑑𝜉
a
𝑢3 0
1
1-𝜉 𝑦23 𝑢12 − 𝑦13 𝑢1 𝑢2 + 𝑦12 𝑢1 𝑢3
2 )𝑑𝜂
+𝜌 [∫ ∫ (𝜉(1 − 𝜂 − 𝜉), 𝜂(1 − 𝜂 − 𝜉), (1 − 𝜉 − 𝜂) 𝑑𝜉 ] (𝑦23 𝑢1 𝑢2 − 𝑦13 𝑢22 + 𝑦12 𝑢2 𝑢3 )
0
0 𝑦23 𝑢1 𝑢3 − 𝑦13 𝑢2 𝑢3 − 𝑦12 𝑢32
1
1-𝜉 −𝑥23 𝑢1 𝑣1 + 𝑥13 𝑢2 𝑣1 + 𝑥21 𝑢3 𝑣1
+𝜌 [∫ ∫ (𝜉(1 − 𝜂 − 𝜉), 𝜂(1 − 𝜂 − 𝜉), (1 − 𝜉 − 𝜂)2 )𝑑𝜂 𝑑𝜉 ] (−𝑥23 𝑢1 𝑣2 + 𝑥13 𝑢2 𝑣2 + 𝑥21 𝑢3 𝑣2 )
0
0 −𝑥23 𝑢1 𝑣3 + 𝑥13 𝑢2 𝑣3 + 𝑥21 𝑢3 𝑣3

𝑝1
1 1-𝜉
(𝑦23 , −𝑦13 , 𝑦12 ) (𝑝2 ) ∫ ∫0 (1 − 𝜂 − 𝜉)𝑑𝜂 𝑑𝜉 = 0……………………………………………………1.9
𝑝3 0
Resolviendo las integrales:
1
1−𝜉
1
∫ ∫ 𝑑𝜂 𝑑𝜉 =
0 2
0
1
1-𝜉
2 1 1
∫ ∫ (𝜉 2 , 𝜉𝜂, 𝜉 − 𝜉 2 − 𝜉𝜂)𝑑𝜂 𝑑𝜉 = ( , , )
0 24 24 24
0
1
1-𝜉
1 2 1
∫ ∫ (𝜉𝜂, 𝜂2 , 𝜂 − 𝜉𝜂 − 𝜂2 )𝑑𝜂 𝑑𝜉 = ( , , )
0 24 24 24
0
 1
1-𝜉
1 1 2
∫ ∫ (𝜉(1 − 𝜂 − 𝜉), 𝜂(1 − 𝜂 − 𝜉), (1 − 𝜉 − 𝜂)2 )𝑑𝜂 𝑑𝜉 = ( , , )
0 24 24 24
0

1
1-𝜉
1
∫ ∫ 𝜉𝑑𝜂 𝑑𝜉 =
0 6
0
1
1-𝜉
1
∫ ∫ 𝜂𝑑𝜂 𝑑𝜉 =
0 6
0
1
1-𝜉
1
∫ ∫ (1 − 𝜂 − 𝜉)𝑑𝜂 𝑑𝜉 =
0 6
0

Sustituyendo estos valores en 1.7,1.8 ,1.9 y uniendo estas ecuaciones resultantes se tiene:

𝑦223 + 𝑥223 -y13 y23 -𝑥13 𝑥23 y12 y23 -𝑥12 𝑥23 𝑢
1
𝜌𝛾 2 + 𝑥2
-y y -𝑥 𝑥 𝑦 -𝑥 𝑥 -y y [ 𝑢2]
2𝑎 13 23 13 23 13 13 21 23 12 13
𝑢
[ y12 y23 -𝑥12 𝑥23 -𝑥12 𝑥13 -y12 y13 𝑦212 + 𝑥212 ] 3
𝑦 𝑢2 − 𝑦 𝑢 𝑢 + 𝑦 𝑢 𝑢 −𝑥 𝑢 𝑣 + 𝑥 𝑢 𝑣 + 𝑥 𝑢 𝑣
23 1 13 1 2 12 1 3 23 1 1 13 2 1 21 3 1
2 1 1
𝜌 2
+ 24 [1 2 1] 𝑦23 𝑢1 𝑢2 − 𝑦13 𝑢2 + 𝑦12 𝑢2 𝑢3 −𝑥23 𝑢1 𝑣2 + 𝑥13 𝑢2 𝑣2 + 𝑥21 𝑢3 𝑣2
1 1 2 2
(𝑦23 𝑢1 𝑢3 − 𝑦13 𝑢2 𝑢3 − 𝑦12 𝑢3 −𝑥23 𝑢1 𝑣3 + 𝑥13 𝑢2 𝑣3 + 𝑥21 𝑢3 𝑣3 )
𝑦23 -𝑦13 𝑦12 𝑃1
1
+ 6 [𝑦23 -𝑦13 𝑦12 ] [𝑃2 ]=0
𝑦23 -𝑦13 𝑦12 𝑃3

Factorizando la matriz del segundo término :

𝑦23 𝑢21 − 𝑦13 𝑢1 𝑢2 + 𝑦12 𝑢1 𝑢3 −𝑥23 𝑢1 𝑣1 + 𝑥13 𝑢2 𝑣1 + 𝑥21 𝑢3 𝑣1

𝑦23 𝑢1 𝑢2 − 𝑦13 𝑢22 + 𝑦12 𝑢2 𝑢3 −𝑥23 𝑢1 𝑣2 + 𝑥13 𝑢2 𝑣2 + 𝑥21 𝑢3 𝑣2
2
(𝑦23 𝑢1 𝑢3 − 𝑦13 𝑢2 𝑢3 − 𝑦12 𝑢3 −𝑥23 𝑢1 𝑣3 + 𝑥13 𝑢2 𝑣3 + 𝑥21 𝑢3 𝑣3 )
𝑢1 (𝑦23 − 𝑦13 + 𝑦12 ) v1 (-x23 + x13 -x21 ) 𝑢1
=[𝑢2 (𝑦23 − 𝑦13 + 𝑦12 ) + v2 (-x23 + x13 -x21 )] [𝑢2 ]
𝑢3 (𝑦23 − 𝑦13 + 𝑦12 ) v3 (-x23 + x13 -x21 ) 𝑢3
Pero 𝑦23 − 𝑦13 + 𝑦12 =𝑦2 − 𝑦3 − 𝑦1 + 𝑦3 + 𝑦1 − 𝑦2 = 0 y -x23 + x13 -x21 =0
Por lo tanto, se obtiene la siguiente ecuación para u:

𝑦223 + 𝑥223 -y13 y23 -𝑥13 𝑥23 y12 y23 -𝑥12 𝑥23 𝑢
1
𝜌𝛾 2 2
-y13 y23 -𝑥13 𝑥23 𝑦13 + 𝑥13 𝑢
-𝑥21 𝑥23 -y12 y13 [ 2 ]
2𝑎
𝑢
[ y12 y23 -𝑥12 𝑥23 -𝑥12 𝑥13 -y12 y13 𝑦212 + 𝑥212 ] 3
𝑦23 -𝑦13 𝑦12 𝑃1
1
+ 6 [𝑦23 -𝑦13 𝑦12 ] [𝑃2 ]=0………………………………………………………………………………2.1
𝑦23 -𝑦13 𝑦12 𝑃3

Similarmente realizando el procedimiento para la ec. 1 se realiza lo mismo para la ec. 2 y

ya que son similares salvo que en lugar de u es v se obtiene lo siguiente:

𝑦223 + 𝑥223 -y13 y23 -𝑥13 𝑥23 y12 y23 -𝑥12 𝑥23 v
1
𝜌𝛾 2 2
-y13 y23 -𝑥13 𝑥23 𝑦13 + 𝑥13 v
-𝑥21 𝑥23 -y12 y13 [ 2 ]
2𝑎
v
[ y12 y23 -𝑥12 𝑥23 -𝑥12 𝑥13 -y12 y13 𝑦212 + 𝑥212 ] 3
-x23 x13 x12 𝑃1
1
+ 6 [-x23 x13 x12 ] [𝑃2 ]=0………………………………………………………………………………2.2
-x23 x13 x12 𝑃3

Para la ecuación 3 se realizan los siguientes cálculos:

𝜕𝑢 𝜕𝑣 𝜕𝑢 𝜕𝜉 𝜕𝑢 𝜕𝜂 𝜕v 𝜕𝜉 𝜕v 𝜕𝜂
∬ ( + ) 𝑝̅ 𝑑𝑥 𝑑𝑦 = ∬ ( + + + ) 𝑝̅ 𝑑𝑥 𝑑𝑦
𝜕𝑥 𝜕𝑦 𝜕𝜉 𝜕𝑥 𝜕𝜂 𝜕x 𝜕𝜉 𝜕y 𝜕𝜂 𝜕y
𝑠 𝑠

u1 v1
1 1
=∬𝑆 (𝑎 (𝑦23 , −𝑦13 , 𝑦12 ) ( 2 ) + 𝑎 (-x23 , x12 , x21 ) (v2 )) Nj 𝑎 𝑑𝜂 𝑑𝜉
u
u3 v3
Para cada función de forma Nj con j=1,2,3 se tiene:
u1 v1 1
1−𝜉
[(𝑦23 , −𝑦13 , 𝑦12 ) (u2 ) + (-x23 , x12 , x21 ) (v2 )] ∫ ∫ Nj 𝑑𝜉 𝑑𝜂
0
u3 v3 0

Pero como se calculó anteriormente:

1 1 1
1-𝜉 1-𝜉 1-𝜉
1
∫ ∫ 𝜉 𝑑𝜉𝑑𝜂 = ∫ ∫ 𝜂 𝑑𝜉𝑑𝜂 = ∫ ∫ (1 − 𝜂 − 𝜉)𝑑𝜂 𝑑𝜉 =
0 0 0 6
0 0 0
Así se tienen las siguientes ecuaciones:
u1 v1
1
6
[(𝑦23 , −𝑦13 , 𝑦12 ) (u2 ) + (-x23 , x13 , x21 ) (v2 )] = 0 para 𝑁1 = 𝜉
u3 v3
u1 v1
1
6
[(𝑦23 , −𝑦13 , 𝑦12 ) (u2 ) + (-x23 , x13 , x21 ) (v2 )] = 0 para 𝑁2 = 𝜂
u3 v3
u1 v1
1
[(𝑦23 , −𝑦13 , 𝑦12 ) (u2 ) + (-x23 , x13 , x21 ) (v2 )] = 0 para 𝑁3 = (1 − 𝜂 − 𝜉)
6 u3 v3
En forma matricial:
𝑦23 -𝑦13 𝑦12 u1 -x23 x13 x21 v1
1 1
[𝑦23 -𝑦13 𝑦12 ] [u2 ] + [-x23 x13 x21 ] [v2 ]=0………………………………………………2.3
6 6
𝑦23 -𝑦13 𝑦12 u3 -x23 x13 x21 v3
Y finalmente sumando las ecuaciones 2.1,2.2 y 2.3 se obtiene el siguiente ELEMENTO
MAESTRO:
𝑦23 2 + 𝑥23 2 −𝑦23 𝑦13 − 𝑥23 𝑥13 𝑦23 𝑦12 − 𝑥23 𝑥12 𝑦 −𝑦13 𝑦12 𝑢1
𝜌𝛾 1 23
{ [−𝑦23 𝑦13 − 𝑥23 𝑥13 𝑦13 2 + 𝑥13 2 −𝑥12 𝑥13 − 𝑦13 𝑦12 ] + [𝑦23 −𝑦13 𝑦12 ]} [𝑢2 ]
2𝑎 6 𝑦 −𝑦13 𝑦12 𝑢3
𝑦23 𝑦12 − 𝑥23 𝑥12 −𝑥12 𝑥13 − 𝑦13 𝑦12 𝑦12 2 + 𝑥12 2 23

𝑦23 2 + 𝑥23 2 −𝑦23 𝑦13 − 𝑥23 𝑥13 𝑦23 𝑦12 − 𝑥23 𝑥12 −𝑥23 𝑥13 𝑥21 𝑣1
𝜌𝛾 1
+ { [−𝑦23 𝑦13 − 𝑥23 𝑥13 𝑦13 2 + 𝑥13 2 −𝑥12 𝑥13 − 𝑦13 𝑦12 ] + [−𝑥23 𝑥13 𝑥21 ]} [𝑣2 ]
2𝑎 6 −𝑥 𝑥13 𝑥21 𝑣3
𝑦23 𝑦12 − 𝑥23 𝑥12 −𝑥12 𝑥13 − 𝑦13 𝑦12 𝑦12 2 + 𝑥12 2 23

𝑦 − 𝑥23 𝑥13 − 𝑦13 𝑦12 − 𝑥12 𝑝1 0

1 23
+ [𝑦23 − 𝑥23 𝑥13 − 𝑦13 𝑦12 − 𝑥12 ] [𝑝2 ] = [0]
6 𝑦 −𝑥 𝑥13 − 𝑦13 𝑦12 − 𝑥12 𝑝3
23 23 0