Obtención de caracteristica estática de sistemas de control automatico - simulink

Depósito
× = = − = = 1 × 3 × × (tan ) + × × × × tan + × × + +

Caracteristica estatica deposito 45

40

Altura teorica Qreq=0 m3/s Altura experimental Qreq=0 m3/s Altura teorica Qreq=2 m3/s35

Altura experimental Qreq= 2 m3/s Altura teorica Qreq= 4 m3/s Altura experimental Qreq= 4 m3/s

30

Altura de equilibrio (m3/s)

25

20

15

10

5

0

0

0.5

1

1.5

2Caudal de entrada equilibrio (m3/s)

2.5

3

3.5

4

Intercambiador
ρ × Cp × V ×
dT = (Cp × ρ × Q1 × T1 ) − (Cp × ρ × Q2 × T2 ) − U × A × (T2 − Tref ) dt

ρ × Cp × Q2 × T2 + ρ × Cp ×QB × T0 = ρ × Cp × Q3 × T3
Q0 = Q2 + Q3 = Q1 + QB
eq 1

Q1 = Q2
  + Q0 − Q1eq × T0eq  

T3eq

eq  ρ × Cp × Q1eq × T1eq + U × A × Tref Q ×  ρ × Cp × Q1eq + U × A  = eq Q0

(

)Caracteristica estatica del intercambiador con un incremento de caudal de entrada 110 Tª salida teorica incremento incremento 0 Tª salida experimental incremento 0 100 Tª salida teorica incremento +0.1m3/s Tª salida experimental incremento +0.1 m3/s 90 Tsalida equilibrio (ºC) Tª salida teorica incremento -0.1m3/s Tª salida teorica -0.1 m3/s 80

70

60

50

40

0

0.1

0.2

0.30.4
3

0.5

0.6

0.7

Fext equilibrio (m /s)

Tsalida equilibrio (ºC)

Caracteristica estatica del intercambiador con un incremento de temperatura a la entrada 130 Tª salida teoricaincremento 0ºC 120 Tª salida experimental incremento 0ºC Tª salida teorica incremento +20ºC 110 Tª salida experimental incremento +20ºC Tª salida teorica incremento -20ºC Tª salida experimental incremento-20ºC 100 90 80 70 60 50 40

0

0.05

0.1

0.15

0.2

0.25

0.3
3

0.35

0.4

0.45

0.5

Fext equilibrio(m /s)

Motor
ui (t ) = La × dii (t ) + Ra × ii (t ) + fcem(t ) dtfcem(t ) = K b × ω m (t ) Pm (t ) = K m × ii (t ) × K e × ie = K a × ii (t )

Pm (t ) = J m ×
Pleq (t ) =

dω m (t ) + Bm × ω m (t ) + Pleq (t ) dt

d ω L (t ) 1  + BL × ω L (t ) −...