Gzoh

Descargar como docx, pdf o txt
Descargar como docx, pdf o txt
Está en la página 1de 4

 Modelo estático (Teoría de Bernoulli)

1
𝑚 ∗ 𝑣 2 = 𝑚 ∗ 𝑔 ∗ ℎ → 𝑣√2𝑔 ∗ ℎ
2
𝑄 = 𝑠 ∗ 𝑣 = 𝑠√2𝑔 ∗ ℎ

 Flujo por una válvula

𝑄𝑣 = 𝑘𝑣 ∗ 𝐴𝑠 √∆𝑃
Donde:

Qv = Flujo a través de la válvula

Kv= constante

∆𝑃=Presión diferencial a través de una salida P2-P1

 La salida es proporcional

𝑄𝑆 = 𝐾2 ∗ 𝑎2 ∗ √2 ∗ 𝑔 ∗ ℎ

 Entrada

𝑄𝑒 = 𝐾1 ∗ 𝑎1

 Se controlará la válvula de entrada, considerando 𝑎2 como perturbación


 Se tomarán en cuenta los siguientes parámetros.
0.05 𝑚2
𝐴𝑟𝑒𝑎 𝑑𝑒𝑙 𝑡𝑎𝑛𝑞𝑢𝑒 = 0.5𝑚2 𝑘1 = 𝑎1 = 0.6
𝑠

𝑚2
𝐻𝑚𝑎𝑥 𝑑𝑒𝑙 𝑡𝑎𝑛𝑞𝑢𝑒 = 1𝑚 𝑘2 = 0.015 𝑠
𝑎2 = 0.5

 Como ecuación principal tenemos.


𝑑ℎ
𝐴∗ = 𝑄𝑒 − 𝑄𝑠 (1)
𝑑𝑡
Si 𝑄𝑒 = 𝑘1 ∗ 𝑎1 ∧ 𝑄𝑠 = 𝑘2 ∗ 𝑎2 (2) y (3)

Reemplazamos las ecuaciones (2) y (3) en (1)


𝑑ℎ
𝐴∗ = 𝑘1 ∗ 𝑎1 − 𝑘2 ∗ 𝑎2 √2 ∗ 𝑔 ∗ ℎ
𝑑𝑡
Esta ecuación no es lineal puesto que la raíz le impide su linealidad. Se necesitará un punto de
equilibrio

𝑄𝑒 = 𝑄𝑠

𝑘1 ∗ 𝑎1 = 𝑘2 ∗ 𝑎2 ∗ √2 ∗ 𝐻𝑜
2
𝑘1 ∗ 𝑎1
𝐻𝑜 = ( )
𝑘2 ∗ 𝑎2 ∗ √2 ∗ 𝑔

Reemplazamos los parámetros

𝐻𝑜 = 0.8 m

Se linealizará por Taylor manteniendo 𝑎2 𝑐𝑡𝑒


𝑑ℎ
𝑓(𝑎1 , ℎ) = 𝐴 ∗
𝑑𝑡
𝜕𝑓 𝜕𝑓
𝑓(𝑎, ℎ) = 𝑓(𝑎
̅̅̅, ̅
1 ℎ) + ∆𝑎1 + ∆ℎ
𝜕𝑎1 𝜕𝑎1

Si 𝑎1 , ∆ℎ = ℎ − ℎ̅
∆𝑎1 = 𝑎1 − ̅̅̅

𝑘2 ∗ 𝑎2 ∗ √2 ∗ 𝑔
𝑎1 − 𝑘2 ∗ 𝑎2 √2 ∗ 𝑔 ∗ ℎ̅ + 𝑘1 ∗ ∆𝑎1 −
𝑓(𝑎1 , ℎ) = 𝑘1 ∗ ̅̅̅ ∗ ∆ℎ
2√ℎ
𝑑∆ℎ 𝑘2 ∗ 𝑎2 √2 ∗ 𝑔
𝐴∗ = 𝑘1 ∗ ∆𝑎1 − ∗ ∆ℎ
𝑑𝑡 2 ∗ √ℎ
Aplicamos Laplace

𝑘2 𝑎2 √2 ∗ 𝑔
𝐴 ∗ 𝑠 ∗ ℎ(𝑠) = 𝑘1 ∗ 𝑎1 (𝑠) − ℎ(𝑠)
2√ℎ
Salida sobre entrada del sistema
ℎ(𝑠) 𝑘1
=
𝑎1 (𝑠) 𝑘2 𝑎2 √2 ∗ 𝑔
𝐴∗𝑠+
2√ℎ
Reemplazamos los parámetros
0.1
𝐺(𝑠) =
𝑠 + 0.0375
Aplicamos un retenedor de orden cero (zoh)

𝑧−1 1 0.1
𝐺𝑧𝑜ℎ (𝑧) = ( )𝒵( ( ))
𝑧 𝑠 𝑠 + 0.0375

𝑧−1 8 1 1
𝐺𝑧𝑜ℎ (𝑧) = ( )( ( − ))
𝑧 3 𝑠 𝑠 + 0.0375
8 𝑧−1
𝐺𝑧𝑜ℎ (𝑧) = (1 − )
3 𝑧 − 𝑒 −0.0375𝑇
Si T=0.1
8 𝑧−1
𝐺𝑧𝑜ℎ (𝑧) = (1 − )
3 𝑧 − 0.9962
0.0304
𝐺𝑧𝑜ℎ (𝑧) =
3𝑧 − 2.9886
0.01013
𝐺𝑧𝑜ℎ (𝑧) =
𝑧 − 0.9962

Obtención de valores para cada n tiempo de muestreo con el método de la división directa

A partir de la fórmula de zoh se obtiene la función de transferencia que responde a un escalón.


Como:

1 𝑧
𝒵( ) = ( )
𝑠 𝑧−1
Entonces la función en tiempo discreto que se analiza es:

𝑧−1 1 0.1 𝑧
𝐹(𝑧) = ( )𝒵( ( )) ( )
𝑧 𝑠 𝑠 + 0.0375 𝑧−1

1 0.1
𝐹(𝑧) = 𝒵 ( ( ))
𝑠 𝑠 + 0.0375

8 1 1
𝐹(𝑧) = ( ( − ))
3 𝑠 𝑠 + 0.0375

8 𝑧 𝑧
𝐹(𝑧) = ( ( − ))
3 𝑧 − 1 𝑧 − 𝑒 −0.00375

0.098148 𝑧
𝐹(𝑧) = (( ))
(𝑧 − 1)(𝑧 − 0.96319)

0.098148 𝑧 −1
𝐹(𝑧) = (( ))
(1 − 1,96319𝑧 −1 + 0.96319𝑧 −2 )

0.098148𝑧 −1⁄1 − 1,96319𝑧 −1 + 0.96319𝑧 −2

0.098148𝑧 −1
0 + 0.19268𝑧 −2 − 0.0945351𝑧 −3⁄
1 − 1,96319𝑧 −1 + 0.96319𝑧 −2
0.098148𝑧 −1 + 0.19268𝑧 −2
0.28374𝑧 −3 ⁄1 − 1,96319𝑧 −1 + 0.96319𝑧 −2

0.098148𝑧 −1 + 0.19268𝑧 −2 + 0.28374𝑧 −3

Entonces al hacer las 3 primeras iteraciones se obtiene:

T F(nT)
0.1 0.098148
0.2 0.19268
0.3 0.28374

También podría gustarte