Ejercicio de Fijacion Sistema Transitorio #2

Descargar como pdf o txt
Descargar como pdf o txt
Está en la página 1de 41

Respuesta de Régimen Transitorio

EJERCICIO 4.1.

Calcular los valores de K y Kh para que el sistema tenga una respuesta con un
sobreimpulso del 20% y un tiempo de 1sg.
R(s) G(s) Y(s)
+-
H(s)
K
G (s) ; H(s) 1 K hs
s(s 1)

K
s(s 1) K
M (s)
K (1 K h s) s2 (1 KK h )s K
1
s(s 1)

2
1
Mp e 0.2 = 0.45
2
tp 1sg wd rad / sg wd wn 1 w n = 3.52 rad/sg
wd
3.52 2 K
M(s) 2
M (s)
s 3.168s 3.52 2 s 2
(1 KK h )s K
K 3.52 2 12.39

1 KK h 3.16

Kh 2.16 / 12.39 0.178

1
Problemas de Ingeniería de Sistemas: Sistemas Continuos. Conceptos básicos.

EJERCICIO 4.2.

Dada la respuesta del siguiente sistema a una entrada escalón unidad, calcular k, f y M.
x
k
M F
f

3.285

1 w 2n
G (s) 2
G (s) C
Ms fs k s2 2 w ns w 2n

Comparando ambas expresiones se tiene:

wn k M; f (2 Mk ) ; C 1 / k

c( t p ) c( ) 3.285 3
Mp 0.095
c( ) 3

2
1
Mp e 0.095 = 0.6

tp 2sg wn 1.96 rad / sg


2
wn 1

x( ) lim sX(s) 3
s 0

2
Respuesta de Régimen Transitorio

1 1
lim sG (s) lim 2
3 C 1/ k 3
s 0 s s 0 Ms fs k

k 1 / 3N / m

k 1
wn k M 1.96 M 2
1.96 3 1.96 2
M 0.087 Kg

1
f (2 Mk ) 0.6 0.6 2 Mk 0.6 2 0.087
3
f 0.204 N /( m / sg)

EJERCICIO 4.3.

Para el sistema de la figura siguiente, calcular el valor de K y p para que cumpla:

Mp 5% ts 4sg
R(s) Y(s)
G(s)
+-
K
G (s)
s(s p)

K w2n
M (s ) M (s )
s2 ps K s2 2 w ns w2n

jw
2
1
Mp e 0.05 0.69

Re
ts 4sg w n 1.13 rad / sg
wn -0.78
ArcCos( ) 46.36º

wn 0.69 * 1.13 0.78

K
M(s) 2
s ps k

K w 2n 1.27

p 2 n 2 0.69 1.13 1.56

3
Problemas de Ingeniería de Sistemas: Sistemas Continuos. Conceptos básicos.

EJERCICIO 4.4.

25
Para el sistema de la figura siguiente, donde G (s)
s(s 6)
Calcular: - Tiempo de subida
- Tiempo de pico.
- Sobreimpulso.
- Tiempo de asentamiento.

R(s) Y(s)
G(s)
+-

25
G (s) s(s 6) 25
M (s)
1 G (s) 25 s 2
6s 25
1
s(s 6)

wn 5rad / sg; 0.6

- Tiempo de Subida:
tr
wd

ArcCos( ) ArcCos(0.6) 0.927 rad / sg

2
wd wn 1 5 1 0.6 2 4rad / sg

0.927
tr 0.55sg
4

- Tiempo de Pico:
tp
wd

tp 0.785sg
4

- Sobreimpulso:
2
1
Mp e

0.6
1 0.6 2
Mp e 0.095 9.5%

4
Respuesta de Régimen Transitorio

Tiempo de asentamiento:
ts
wn

ts 1.05sg
0.6 * 5

EJERCICIO 4.5.

Para el sistema del ejercicio 1.1. suponiendo:

M 1kg; f 1N /(m / sg);


k 1N / m;
F( t ) 1N;

Obtener la respuesta temporal del sistema ante un escalón unitario

Del ejercicio 1.1 se tiene:


X(s) 1
2
F(s) s M fs k

Aplicando una fuerza escalón unitario la salida del sistema será:


1 1
X(s) 2
(s s 1) s

1 s 0.5 0.5 0.866


X(s)
s (s 0.5) 2 0.866 2 0.866 (s 0.5) 2 0.8662

5
Problemas de Ingeniería de Sistemas: Sistemas Continuos. Conceptos básicos.

Aplicando tablas de transformadas:

0.5 t 0.5 t
x(t) 1 e Cos(0.866t ) 0.576e Sen (0.866 t )

EJERCICIO 4.6.

Para el sistema del ejercicio 1.3. suponiendo:

M 1 Kg; M 0.5 Kg; f = f = 1 N m/sg ; k =1N m


1 2 1 2

Calcular la respuesta temporal del sistema ante una entrada escalón unitario.

En el ejercicio 2.3 se obtuvo las funciones de transferencia que relacionaban las variables
del sistema:
V1 (s) M 2s f 2 k s
F(s) sM 1 f 1 f 2 M 2 s f 2 k s f 22

V2 (s) f2
F( s ) sM 1 f1 f 2 M 2s f 2 ks f 22

Sustituyendo por los valores:

V1 (s) 0.5s 1 1 s 0.5s 2 s 1


3 2
F(s) s 2 0.5s 1 1 s 1 0.5s 2s 2s 2

6
Respuesta de Régimen Transitorio

0.5s 2 s 1 1
V1 (s) 3 2
(0.5s 2s 2s 2) s

0.5 0.212 0.288s 0.086


V1 (s) 2
s s 3.13 s 0.435 1.043 2

Aplicando las tablas de transformadas de Laplace:

3.13t 0.435 t
v 1 ( t ) 0.5 0.212e 0.353 cos(1.043t 0.613)e

V2 (s) 1 s
F(s) s 2 0.5s 1 1 s 1 0.5s 3 2s 2 2s 2

s 1
V2 (s)
(0.5s 3 2s 2 2s 2) s

0.239 0.239s 0.541


V2 (s) 2
s 3.13 s 0.435 1.043 2

Aplicando las tablas de transformadas de Laplace:

3.13t 0.435 t
v 2 ( t ) 0.239e 0.663 sen(1.043t 0.369)e

0.7

0.6 v1(t)

0.5

0.4

0.3
v2(t)
0.2

0.1

-0.1
0 2 4 6 8 10
Tiempo (s)

7
Problemas de Ingeniería de Sistemas: Sistemas Continuos. Conceptos básicos.

EJERCICIO 4.7.

Para el sistema del ejercicio 1.10. si se aplica una fuerza f(t) de 1 Newton, calcular:

- Posición del deslizador en estado estacionario.


- Posición máxima a la que llega.

En el ejercicio 1.10. se obtuvo la función de transferencia del sistema linealizada en torno del
punto de reposo donde x = 1:
F(s)
X(s) 2
3 s 10 s 10

Valor en régimen estacionario:


1
F(s) s
lim s X(s) lim s lim s 0.1
s 0 s 0 3 s2 10 s 10 s 0 3 s 2 10 s 10

Luego
x est 1 0.1 1.1

Valor máximo:
1
X(s) 3
F(s) 10 10
s2 s
3 3

Comparando esta expresión con la de un sistema de segundo orden se tiene:

10
n
3

10
2 n
3

10 3 30
0.9128
2 3 10 6

0.9128
2
1 0.9128 2 3
Mp e 1
e 0.88 10 0.088%

3
X max 11
. 11
. 0.88 10 11001
.

8
Respuesta de Régimen Transitorio

EJERCICIO 4.8.

1) Obtener las funciones de transferencia G1, G2, G3, G4.y total del sistema.
2) Introducir un controlador proporcional que haga al sistema comportarse como si fuese de
2º orden, con sobreimpulso, ante entrada escalón, del 4.321%

G1(s)
+
R(s) + G3(s) G4(s)
- +
G2(s)

Para obtener la función de transferencia del sistema de la figura se han realizado una serie
de ensayos a los bloques G1, G2, G3 y G4 que lo constituyen.
Al bloque G1 se le ha sometido a una señal de entrada rampa unidad, obteniéndose la
respuesta de la gráfica:
4

0
0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1

Al bloque G2 se le ha introducido una señal escalón unidad, respondiendo como se puede


ver en la gráfica:
1.5

0.5

0
0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1

9
Problemas de Ingeniería de Sistemas: Sistemas Continuos. Conceptos básicos.

Al bloque G3 se le ha introducido una señal impulso unidad, obteniéndose:

0.14

0.12
0.11
0.1

0.08

0.06

0.04

0.02

0
0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1

Al bloque G4 se le ha introducido una señal escalón unidad, obteniéndose


1

0.9

0.8

0.7

0.6

0.5

0.4

0.3

0.2

0.1

0
0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1

1) El bloque G1(s) corresponde con un derivador con ganancia 3: G1 (s) 3s

El bloque G2(s) corresponde con una ganancia unitaria: G 2 (s) 1

0.11 1
El bloque G3(s) corresponde con un integrador de ganancia 0.11: G 3 (s)
s 9s

1
El bloque G4(s) corresponde con un integrador de ganancia unidad: G 4 (s)
s

Luego la cadena directa del sistema quedará:

1
G T (s) (G 1 (s) G 2 (s)) G 3 (s) G 4 (s) (3s 1)
9s 2

10
Respuesta de Régimen Transitorio

3s 1
G T (s) 9s 2 K (3s 1) K (3s 1)
M (s)
1 G T (s) 3s 1 2
9s 3Ks K 9 s 2 3 Ks K
1 K 2
9s 9 9

2) Suponiendo que el cero no afecta a la respuesta:

2 2
1 2 2 · 2 2
Mp e 0.04321 1
1
1

2
2 1 0.7071

K
n
3

3
K 2 n
9

3 K
K 2 K 2
9 3

K 2

EJERCICIO 4.9.

Obtener la respuesta del sistema del ejercicio 1.9. ante una entrada escalón unitario. Analizar
el posible uso de un sistema equivalente reducido.

En el ejercicio 1.9. se obtuvo la función de transferencia del sistema que relaciona el


desplazamiento angular del eje de salida con el desplazamiento angular del eje de referencia:
Ea(s) (s)
(s) E(s) m
c
(s)
r Amplific. Motor Engrane
+-

c 50000
3 2
r s 50.5s 1725s 50000

La transformada de la salida c, para esta entrada escalón unitario es:


50000 1
c (s ) 3 2
s 50.5s 1725s 50000 s

11
Problemas de Ingeniería de Sistemas: Sistemas Continuos. Conceptos básicos.

Expansión en fracciones parciales:

1 0.54 s 35.3 22.34 35.3


c (s) 0.458
s s 39.1 (s 5.7) 2 35.32 35.3 (s 5.7) 2 35.32

y la respuesta del sistema ante una entrada escalón queda:


39.1t 5.7 t 22.34 5.7 t
c (t) 1 0.54e 0.458 Cos(35.3t )e Sin (35.3t )e
35.3
Que de forma gráfica sería:

La respuesta es muy similar al de un sistema de segundo orden subamortiguado, con las


siguientes características:

Mp = 41%; tr = 0.07sg.; tp = 0.11sg.

El motivo de que la respuesta sea similar a la de un sistema de segundo orden está en que,
para este valor de ganancia del amplificador, la función de transferencia de lazo cerrado tiene
dos polos complejos conjugados en -5.7 j35.3 (dominantes) y uno real en -39.1. De esta
forma podemos ver que el polo del eje real tiene una constante de tiempo pequeña comparada
con la de los otros dos polos y por eso la forma de la respuesta es dominada por los
complejos.

Se podía pensar en aproximar el sistema por una simplificación con solo los polos
complejos, como la siguiente:
1278.6
M (s) 2
s 11.4s 1278.6

El problema está en que será una aproximación muy burda ya que la relación entre las
distancias al eje imaginario, de los polos complejos con el polo real, es menor de 10 (6.8)

12
Respuesta de Régimen Transitorio

luego no sería despreciable.

Para comprobar que no se puede realizar la aproximación (en principio) comparamos la


respuesta, ante un escalón unidad, del sistema real con la del sistema reducido:

Se puede ver como en la parte transitoria se produce un error apreciable.

EJERCICIO 4.10.

Para el sistema del ejercicio 1.12., suponiendo que el amplificador K toma el valor 0.6,
dibujar de forma aproximada la respuesta del sistema x(t) ante una entrada escalón unitario en
la referencia de posición r(t). Indicar los valores numéricos de los principales parámetros de
dicha respuesta.

En el ejercicio 1.12. se obtuvo el siguiente diagrama de bloques del sistema:


Bomba Embolo
Amplificador
r(t) e(t) v(t) p(t) x(t)
K
1.2 0.25
2
+_ 2s 1 5s 2.5s 1
c(t)
Transductor
10

1
R (s)
s

1.2 0.25 0.3


G (s) 2 3
2s 1 5s 2.5s 1 10s 10s 2 4.5s 1

13
Problemas de Ingeniería de Sistemas: Sistemas Continuos. Conceptos básicos.

0.3
0.6
K G (s) 10s 3
10s 2 4.5s 1 0.18
M (s)
1 K G (s)H(s) 0.3 10s 3
10s 2 4.5s 2.8
1 0.6 3
10
10s 10s 2 4.5s 1

Los polos de la función de transferencia de lazo cerrado se encuentran en:

p1, 2 0.0718 j0.567


p3 0.856

Como p3 está lo suficientemente alejado para poder despreciarlo frente a p1 y p2 se


considerará el sistema equivalente de segundo orden formado por los dos polos dominantes p1
y p2.
Para estos polos se obtiene los principales parámetros del sistema:

2
p1, 2 j d n j n 1 0.0718 j0.567

n 0.0718

2
n 1 0.567

2
0.1266 1.01604 0.01604 0.1256
2
1

0.0718
n 0.5717
0.1256

cos ar cos 1.445rad

1.445
tr 2.99s
d 0.567

tp 5.54s
d 0.567

ts 43.75s
0.0718

0.1256
2
1 1 0.1256
Mp e e 0.6718 67.18%

Valor final:
0.18 1
lim sM (s)R (s) s 3 2
0.064
s 0 10s 10s 4.5s 2.8 s

14
Respuesta de Régimen Transitorio

La figura muestra la respuesta la respuesta del sistema completo.


Respuesta al escalón unitario
0.1
0.09
0.08
0.07
0.06
Amplitud

0.05
0.04
0.03
0.02
0.01
0
0 10 20 30 40 50 60 70 80
Tiempo(s)

El sistema equivalente considerado para el cálculo de los parámetros ha sido:


k red k red
M red (s)
(s 0.0718)2 0.567 2
s 2
0.1436s 0.3266

k 1 k red
lim sM (s)R (s) s 2
0.064 0.064 k red 0.0209
s 0 s 0.1436s 0.3266 s 0.3266

0.0209
M red (s) 2
s 0.1436s 0.3266
La respuesta del sistema reducido es:
Respuesta al escalón del equivalente reducido
0.12

0.1

0.08
Amplitud

0.06

0.04

0.02

0
0 10 20 30 40 50 60 70 80
Tiempo (s)

15
Problemas de Ingeniería de Sistemas: Sistemas Continuos. Conceptos básicos.

Comparando ambas respuestas en la misma gráfica pueden observarse las diferencias


existentes al despreciar el polo más alejado del origen.

Comparativa de las respuestas de ambos sistemas


0.12

0.1

0.08
Amplitud

0.06

0.04

0.02

0
0 10 20 30 40 50 60 70 80
Tiempo(s)
La respuesta en línea continua corresponde con el sistema total y la de línea discontinúa
con el equivalente reducido.

EJERCICIO 4.11.

Para el sistema de control cuyo diagrama de bloques se muestra en la figura:

X(s) E(s) U(s) Y(s)


K G(s)
+-

H(s)

3(s 0.66) 10
G (s) H(s)
s 2 4s 5 (2s 8)(s 6)

1. Obtener la transformada inversa de la respuesta del sistema para una entrada escalón
unitario cuando K=1.
2. Calcular un sistema reducido equivalente al de lazo cerrado, para K= 0.1
3. Calcular los parámetros de la respuesta del sistema para una entrada escalón unitario:
- Tiempo de subida (tr)
- Tiempo de pico (tp)
- Máximo sobreimpulso (Mp)
- Tiempo de asentamiento (ts)

16
Respuesta de Régimen Transitorio

1. Transformada inversa.
3s 2
K· G (s) s 4s 5 2 (3s 2)(s 4)(s 6)
M (s)
1 K· G (s)H(s) 3s 2 5 (s 2
4s 5)(s 4)(s 6) 15s 10
1 2
s 4s 5 (s 4)(s 6)

1 (3s 2)(s 4)(s 6) 3(s 0.667)(s 4)(s 6)


Y (s )
s (s 2 4s 5)(s 4)(s 6) 15s 10 s(s 1.546)(s 7.132)((s 2.661) 2 2.170 2 )

Cuya descomposición en fracciones simples es:

0.369 0.561 0.0699 0.86s 1.209


Y(s)
s s 1.546 s 7.132 (s 2.661) 2 2.1702

0.86s 1.209 s 1.406 s 1.406 1.255 1.255


0.86 0.86
(s 2.661)2 2.1702 (s 2.661) 2 2.1702 (s 2.661) 2 2.1702

s 1.406 1.255 1.255 s 2.661 1.255


0.86 0.86
(s 2.661)2 2.170 2 (s 2.661)2 2.1702 (s 2.661)2 2.170 2

s 2.661 1.255 2.170


0.86
(s 2.661)2 2.170 2 2.170 (s 2.661) 2 2.1702

s 2.661 2.170
0.86 0.5783
(s 2.661)2 2.170 2 (s 2.661) 2 2.170 2

Aplicando las tablas de transformada de Laplace:

1 1 at
1 e
s s a

s a at at
e cos t e sen t
(s a ) 2 2
(s a ) 2 2

0.369 0.561 0.0699 s 2.661 2.170


Y(s) 0.86 0.5783
s s 1.546 s 7.132 (s 2.661) 2 2.170 2 (s 2.661) 2 2.1702

1.546 t 7.138 t 2.661t


y( t ) 0.369 0.561e 0.0699e 0.86e cos(2.170 t ) 0.5783 sen(2.170t )

1.546 t 7.138 t 2.661t


y( t ) 0.369 0.561e 0.0699e 0.9934e cos(2.170 t 0.5243)

17
Problemas de Ingeniería de Sistemas: Sistemas Continuos. Conceptos básicos.

2. Sistema reducido equivalente. Factorizando M(s) para K=0.1 se tiene:

3s 2
0.1
K· G (s) s 2
4s 5 (3s 2)(s 4)(s 6)
M (s)
1 K· G (s)H(s) 3s 2 5 (s 2
4s 5)(s 4)(s 6) 15s 10
1 0.1 2
s 4s 5 (s 4)(s 6)

0.1(3s 2)(s 4)(s 6) 0.3(s 0.66)(s 4)(s 6)


M (s) 4
s 14s 3 69s 2 147.5s 121 ((s 2.16) 2 0.97 2 )(s 3.47)(s 6.2)

Un sistema aproximado, aunque de forma grosera, sería:

0.2184 0.2184
M ' (s) 2 2 2
(s 2.16) 0.97 s 4.32s 5.6

donde se ha ajustado la ganancia del equivalente reducido para tener la misma ganancia
estática que el sistema original (0.039).

Un modelo reducido más ajustado sería:

0.321(s 0.66) 0.321(s 0.66)


M ' ' (s)
(s 2.16) 2 0.97 2 s 2 4.32s 5.6

La diferencia está en el cero en –0.66, porque suma a la respuesta del primer equivalente la
respuesta del sistema de segundo orden a un impulso:

0.321s 0.321 0.66 1


Y' ' (s) M (s) X(s) 2 2
s 4.32s 5.6 s 4.32s 5.6 s

0.321 0.321 0.66 1


Y' ' (s) 2 2
s 4.32s 5.6 s 4.32s 5.6 s

Respuesta al impulso
0.06

0.05

0.04

0.03

0.02

0.01

0
0 0.5 1 1.5 2 2.5 3
Tiempo (s.)

18
Respuesta de Régimen Transitorio

Respuesta al escalón del sistema original y equivalentes reducidos


0.07

0.06
M(s)
0.05
M’(s)
0.04

0.03 M’’(s)

0.02

0.01

0
0 0.5 1 1.5 2 2.5 3
Tiempo (s.)

3. Parámetros obtenidos sobre el equivalente reducido:

2
n 5.6 n 5.6 2.36 rad / s

4.32
2 n 4.32 0.91 arccos arccos 0.91 24.5º 0.427 rad
2 2.36

2
d n 1 2.36 1 0.912 0.978 rad/s

tp 3.21s
d 0.978

0.427
tr 2.77s
d 0.978

ts 1.46s
n 0.91 2.36

2
1
Mp e 0.00101 0.1%

19
Problemas de Ingeniería de Sistemas: Sistemas Continuos. Conceptos básicos.

EJERCICIO 4.12.

El siguiente sistema se utiliza para eliminar durante su transporte parte de la humedad de


una sustancia que se almacena en un silo:

Aire
Vm(t) + Vref(t)
A
(t) _
Vh(t)

h(t)

La señal de error e(t) compara la referencia de humedad Vref(t) con la obtenida del sensor
Vh(t), se amplifica por un valor A, y se utiliza Vm(t) como entrada a un motor de CC que
actúa sobre una válvula que gobierna el paso de una corriente de aire caliente. La ganancia del
amplificador es de A Voltios/Voltio.

La relación entre la posición del eje del motor (t) en grados y la tensión aplicada al mismo
Vm(t) viene dada en la siguiente gráfica:
Vm(t) (t)
1 1

t t
10
La relación entre la humedad del material y(t) (en %) y la posición del eje del motor (t)
(en grados) viene dada por la respuesta impulsional:
(t)

Area=1
t

20
Respuesta de Régimen Transitorio

La humedad es detectada a la salida de la cinta Vh (voltios) Tiempo


transportadora mediante un instrumento empleado (segundos)
como sensor que genera tensión por unidad de 0 0
variación del grado de humedad h. Para calcular la 7.89 0.5
función de transferencia del sensor de humedad, cuya 12.64 1
ecuación diferencial que rige su funcionamiento es 15.56 1.5
17.31 2
d Vh(t) 18.36 2.5
k h(t) a Vh(t)
dt 19 3
19.4 3.5
se realizó el ensayo de someter al sensor a un 19.64 4
cambio brusco de 10 unidades en la humedad 19.78 4.5
obteniéndose los siguientes datos que se muestran en 19.87 5
la tabla adjunta. 19.92 5.5
19.95 6
19.97 6.5
20 8

Dibujar el diagrama de bloques del sistema completo que tiene como entrada la tensión de
referencia Vref(s) y como salida la humedad del material H(s), indicando todas las variables y
funciones de transferencia.

Analizando el sistema representado se van obteniendo las relaciones entre las diferentes
variables del sistema:
E(s) Vref (s) Vh(s)

Vm(s)
Vm(s) A E (s) A
E (s)
La relación entre la posición del eje del motor (t) en grados y la tensión aplicada al mismo
Vm(t) se obtiene de la gráfica correspondiente. En ella se puede ver que cuando la entrada
Vm(t) es un escalón unitario, la salida (t) es una rampa de pendiente 0.1. Luego dicho
elemento puede representarse como un integrador cuya función de transferencia corresponde
con:
0.1
(s) Vm(s)
s

(s) 0.1
Vm(s) s

La relación entre la humedad del material y(t) (en %) y la posición del eje del motor (t)
(en grados) viene dada por la respuesta impulsional representada en la gráfica. Esta respuesta
corresponde a un sistema de primer orden de la forma:

H(s) k
(s) s

21
Problemas de Ingeniería de Sistemas: Sistemas Continuos. Conceptos básicos.

La respuesta temporal de un sistema de primer orden ante una entrada impulso es de la


forma: h (t ) k e t

El valor inicial de la respuesta es 0.1:

h (0) k e0 k 0.1 k 0 .1

El área debe ser igual a 1:


t
t e 1 k 0.1
k e k k 0 1 0.1
0 0

Luego la función de transferencia es:


H(s) 0.1
(s) s 0.1

Para calcular la función de transferencia del sensor de humedad, se tiene la ecuación


diferencial que rige su funcionamiento:
d Vh(t)
k h(t) a Vh(t)
dt

Aplicando transformada de Laplace se tendrá:

k H(s) sVh(s) a Vh(s)

k H(s) (s a) Vh(s)

Vh(s) k k'
H(s) s a s 1

La tabla proporciona información sobre la respuesta de este sistema ante una entrada
escalón de amplitud 10. Como el sistema corresponde con uno de primer orden, de los datos
de la tabla obtenemos el valor de la constante de tiempo y de la ganancia k’. El valor al que
se estabiliza la respuesta para una entrada escalón de 10 unidades es 20.

Para el sistema de primer orden el valor al que se estabiliza la respuesta es:

k' 10
Vh lim s 10k'
s 0 s 1 s

20
Luego: 10k ' 20 k' 2
10
La constante de tiempo corresponde con el tiempo transcurrido hasta que la respuesta
alcanza el 63.2% del valor final. El 63.2 % de 20 es 12.64 que mirando en la tabla
corresponde con t=1 segundo.

22
Respuesta de Régimen Transitorio

Vh(s) 2
Por tanto la función de transferencia es:
H(s) s 1
El diagrama de bloques correspondiente al sistema completo será:
Vref(s) E(s) Vm(s) 0.1 (s) 0.1 H(s)
A
+_ s s 0.1
Vh(s)

2
s 1

EJERCICIO 4.13.

La figura muestra el esquema del control de velocidad del motor de combustión interna
mediante la utilización de una servoválvula hidráulica.
Suministro
Servoválvula
Hidráulico Entrada
Hidráulica
Out In Out Combustible
Amp. dif.
Er i
+ Válvula lineal
Eo -
Ya
Ya qr
wo Motor wo
Carga
Tacómetro Combustion T

Se sabe que el amplificador diferencial tiene una ganancia regulable, de valor K y que la
función que realiza es: i K (E r E o ) . La válvula lineal presenta una relación constante
entre el flujo de mezcla de combustible (qr)a su salida y su desplazamiento (Ya), igual a
10000. Además, se somete al conjunto Motor-Carga a una entrada escalón de amplitud
unidad, obteniéndose la respuesta de la gráfica siguiente:

23
Problemas de Ingeniería de Sistemas: Sistemas Continuos. Conceptos básicos.

Se dispone así mismo de una gráfica dada por el fabricante de la servoválvula hidráulica,
que muestra la respuesta de la misma ante una entrada escalón unidad de intensidad:

- Se ha sometido a la tacodinamo a una entrada rampa de pendiente 500 r.p.m. (velocidad en


el eje creciente linealmente), y se ha obtenido que su respuesta viene dada por la expresión:
E o (rampa 500) 5 t 0.5 0.5 e 10 t
Se pide:
1. Diagrama de Bloques del sistema indicando la posición de la planta, actuador, regulador,
sensor, realimentación, señal de error, de actuación, de entrada, de salida,...etc.
2. Función de transferencia de cada bloque.
3. Función de transferencia de lazo cerrado.
4. Indicar mediante observación de la FTLC como será su respuesta (1er orden, 2º orden,
sobreamortiguado, etc.)

1. Diagrama de bloques.

Señal de Señal de control Señal de actuación


referencia Amplificador
diferencial Actuador Planta
Er i Ya qr Motor / 0
r
Regulador Servo- Válvula
Pot. +_ válvula lineal. Carga

Señal de error Señal de salida


E0
Tacómetro
Señal de Sensor
realimentación

2. Funciones de transferencia.

Motor-carga:

La gráfica corresponde con la respuesta de un sistema de primer orden donde la constante

24
Respuesta de Régimen Transitorio

de tiempo leída de la gráfica es aproximadamente 4 segundos y la ganancia K del sistema


es 0.1.
0 ( s) K 01
. 0.025
Q r (s) s 1 4s 1 s 0.25
Servoválvula:

La gráfica corresponde con la respuesta de un sistema de segundo orden subamortiguado,


donde leyendo de la gráfica se puede ver que el tiempo de subida es de 0.4 segundos y el
tiempo de pico de 0.55 segundos.

tp 0.55 d 5.71rad / s
d tp

tr 0.4 tr d 0.857 rad


d

cos cos 0.857 0.654

2 d
d n 1 n 7.54rad / s
2
1

Luego la función de transferencia será:

Ya (s) K1 2n
I(s) s2 2 ns 2
n

donde K1 corresponde con la ganancia de régimen estacionario.

Ya (s) 0.01 56.85 0.5685


2 2
I(s) s 9.86s 56.85 s 9.86s 56.85

Ya (s) 1
2
I(s) s 16s 100

Válvula lineal:
Q r (s)
10000
Ya (s)

Tacómetro:
10 t
E0 (t) 5t 0.5 0.5e

Tomando transformadas:

5 0.5 0.5 5s 50 0.5s 2 5s 0.5s 2 0


E 0 (s)
s2 s s 10 s 2 (s 10) 2
s (s 10)

25
Problemas de Ingeniería de Sistemas: Sistemas Continuos. Conceptos básicos.

Como esta señal es la respuesta del tacómetro ante entrada rampa de pendiente 500:

500
0 (s)
s2

E 0 (s) 50 s2 0.1
0 (s ) s 2 (s 10) 500 s 10

Otra posible forma de obtener la función de transferencia sería derivando la respuesta del
sistema ante la rampa para obtener la respuesta del sistema ante un escalón:

dE 0 ( t ) 10 t 10 t
5 5e 5(1 e )
dt

Esta respuesta corresponde a un sistema de primer orden con constante de tiempo 0.1 y
ganancia 0.01 (al ser la entrada de amplitud 500).

E 0 (s) K2 0.01 0.1


0 (s ) s 1 0.1s 1 s 10

Detector de error:
I(s) K (E r (s) E 0 (s))

3. Función de transferencia de lazo cerrado.

1 0.025 250K
G T (s) K 2
1000 2
s 16s 100 s 0.25 (s 16s 100)(s 0.25)

0.1
H(s)
s 10

250K
2
G T (s) (s 16s 100)(s 0.25)
M (s)
1 G T (s) H(s) 250K 0.1
1
(s 2 16s 100)(s 0.25) s 10

250K (s 10) 250K (s 10)


M (s) 2 4 3
(s 16s 100)(s 0.25)(s 10) 25K s 26.25s 266.5s 2 1065s 250 25K

4. Como se puede ver en la expresión anterior, la posición de los polos de la función de


transferencia de lazo cerrado depende del valor que tome K. Con lo cual es muy dificil decir
por mera observación cual sería el comportamiento. Habría que dar valores a K o dibujar el
lugar de las raíces.

26
Respuesta de Régimen Transitorio

EJERCICIO 4.14.

Se desea diseñar un controlador de altura para un globo aerostático.


Aire caliente
en interior, a
temperatura T

Quemador

Elevación por
Viento (W) Z

Para ello se dispone de un altímetro, cuya función de transferencia es la unidad, y de la


dinámica del sistema que, debido a la dificultad de encontrar las ecuaciones de la dinámica
del sistema con precisión, se ha obtenido en base a varios experimentos.
El primer experimento tenía la misión de mostrar la relación entre el Flujo de Calor (Q)
aplicado al interior del globo mediante un quemador y la Temperatura que se alcanza en su
interior. Para ello se ha realizado un ensayo de forma que se ha aplicado una señal de entrada
tipo impulsional de amplitud 10(Kcal/sg), obteniéndose la Gráfica Nº 1.
Posteriormente se ha realizado otro ensayo para conocer la relación entre la Temperatura
interior del globo y la altura (Z) que alcanza el mismo. Para esto se ha introducido una señal
impulsional de amplitud 5, obteniéndose la Gráfica Nº 2.

Area=2

Gráfica 1 Gráfica Nº 2

Por último, se ha comprobado que, para caracterizar completamente la dinámica del globo,
es necesario tener en cuenta la modificación de la altura producida por la componente
ascensional del viento (W).
Calcular:
- Diagrama de bloques de la planta.
- Función de transferencia (Altura/Calor de entrada) de la planta.
- Diagrama de bloques en lazo cerrado, indicando la señales de entrada al sistema y
si estas son controlables o no.

27
Problemas de Ingeniería de Sistemas: Sistemas Continuos. Conceptos básicos.

- Diagrama de bloques de la planta.


W(s)

Q(s) T(s) Z1(s) Z(s)


G1(s) G2(s)

Como se observa en la figura anterior, la planta estará formada, según los datos aportados,
por dos bloques que representan la dinámica del sistema y que permiten modelar las
variaciones de la temperatura interior y de la altura del globo. (Comentar que las ecuaciones
que en este ejercicio se han considerado como absolutas son, en realidad “variaciones”).
Además, se incluye la entrada de perturbación producida por la componente vertical del
aire.

- Función de transferencia (altura/calor de entrada) de la planta.


Para calcular la función de transferencia de la planta es necesario obtener la de cada uno de
los bloques que la componen:
a) Función G1(s):
La respuesta temporal de la gráfica Nº 1 es típica de los sistemas de primer orden ante una
entrada impulso, luego será de primer orden. Entonces:

k 1 k at
G 1 (s ) g1 ( t ) y( t ) (t) y( t ) L k e
s a s a

Para t=0, se comprueba que y( t ) k 10 .


Y calculando el área bajo la curva:
at 10 10
area 10 e dt 0 1
0
a a
Y como el área es dato, queda que a=5.
Y teniendo en cuenta que la amplitud del impulso de entrada ha sido 10, la función de
transferencia G1(s) queda:
T(s) 1
G 1 (s)
Q(s) s 5
b) Función G2(s):
La expresión matemática de la curva de la gráfica Nº 2 puede ser aproximada por:
t
y( t ) 0.03 u ( t ) 0.03 e 0.02

Donde el valor de la constante de tiempo del termino exponencial se ha obtenido a partir de


la gráfica, comprobando el tiempo transcurrido para llegar al 63,2% del valor final.
Tomando transformadas de Laplace con C.I. nulas:

0.03 0.03 15
.
Y(s)
s s 50 s(s 50)

Esta es la transformada de la respuesta del sistema G2 ante una señal de entrada impulso de
amplitud 5, luego si dividimos entre la amplitud de la señal de entrada se obtiene la función
de transferencia:

28
Respuesta de Régimen Transitorio

Z 1 (s) 0.3
G 2 (s)
T(s) s(s 50)

Entonces, la función de transferencia de la planta queda:

Z 1 (s) 0.3
G (s)
Q(s) s(s 50)(s 5)

- Diagrama de bloques en lazo cerrado:


W(s)
Controlador
Zref E(s) Q(s) T(s) Z1(s) Z(s)
1 G1(s) G2(s)

En este sistema se observan dos entradas Zref y W. La primera de ellas es la entrada


mediante la cual el operador del globo modifica la posición, es totalmente controlable.
En cuanto a la señal W, es no controlable (además de aleatoria y observable) ya que depende
de la velocidad que alcance el viento en cada instante.

EJERCICIO 4.15.

Dibujar, indicando un número de puntos que permita observarlo correctamente, la respuesta


2(s 1)
del sistema cuya función de transferencia de lazo cerrado es M (s) , ante una entrada
( 2s 1)
tal como la que se muestra en la figura:

+2

0 6 12 18 2 Tiempo

-2

Descomponemos la señal de entrada en intervalos de tiempo y analizamos el tipo de señal.


Así observamos que desde t = 0 a t =6-, la señal de entrada es un escalón de amplitud 2:

Entrada R(s) 2(s 1) Respuesta C(s)


M (s)
2s 1

29
Problemas de Ingeniería de Sistemas: Sistemas Continuos. Conceptos básicos.

luego:
2(s 1) 2 4(s 1)
C1 (s) M (s) R 1 (s)
2s 1 s s(2s 1)

que descomponiendo en fracciones simples, quedará:

4(s 1) A B
s(2s 1) s 2s 1

4(s 1)
A s 4
s(2s 1) s 0

4(s 1)
B (2s 1) 4
s(2s 1) s 0.5

4(s 1) 4 4 4 2
s(2s 1) s 2s 1 s s 0.5

Haciendo ahora la transformada inversa, tendremos:

0 .5 t
c1 ( t ) 4 2e
dando valores a t:

t 0 c1 ( t ) 4.000

t 1 c1 ( t ) 2.7869

t 2 c1 ( t ) 3.2642

t 3 c1 ( t ) 3.5537

t 6 c1 ( t ) 3.900

A partir de este instante la entrada r2(t) es un escalón de valor -4. Luego:

2(s 1) 4 8(s 1) 4(s 1)


C 2 (s) M (s) R 2 (s)
2s 1 s s(2s 1) s(s 0.5)

Expandiendo en fracciones simples:

4(s 1) A B 8 4
s(s 0.5) s (s 0.5) s s 0.5

Hallamos la transformada inversa, resultando:

30
Respuesta de Régimen Transitorio

0.5 t
c2 (t ) 8 4e

Ahora el origen está en el punto (6,0) y la gráfica parte del valor 3.900, luego con estas
consideraciones, construimos la tabla:

t 0 punto(6 ,0) c2 (t ) 8 4 4 ordenada : 4 3.900 0.1

t 1 punto(7,0) c2 (t ) 8 2.4261 5.5739 ordenada : 5.5739 3.9 1.6739

t 2 punto(8,0) c2 (t ) 8 1.4715 6.5285 ordenada : 6.5285 3.9 2.6285

t 3 punto(9,0) c2 (t ) 8 0.8925 7.1075 ordenada : 7.1075 3.9 3.2075

t 6 punto(12 ,0) c2 (t ) 8 0.1991 7.8 ordenada : 7.8 3.900 3.9

Y volvería a repetirse el proceso, quedando la gráfica de la respuesta:

+4
+3.9

+2

0.1
-0.1 6 12 18 24 Tiempo

-2

-3.9
-4

31
Problemas de Ingeniería de Sistemas: Sistemas Continuos. Conceptos básicos.

EJERCICIO 4.16.

Para el sistema que se muestra en la figura siguiente, se han realizado una serie de
experiencias con el fin de obtener las funciones de transferencia de cada uno de los bloques
que lo componen:
W1(s)
+
R(s) C(s)
Amplificador Regulador Actuador Planta
+ +
-
+
sensor
+
W2(s)

Se ha observado que la integral de la respuesta del regulador, cuando se introduce una


señal tipo escalón unitario, viene dada por:
t
3 0.2
3t e
5

Para conocer la respuesta del sensor, se ha introducido también un escalón unidad y se ha


graficado su salida, siendo ésta la que se muestra a continuación:

2.296

1.5
Respuesta

0.5

0
0 0.5 1 1.5 2 2.5 3 3.5 4
1.226
Tiempo (seg.)

Atendiendo a los datos del fabricante, se sabe que la función de transferencia del actuador
es: (s + 3), y se ha obtenido el flujograma de la planta, que es el que se muestra:

1
R'(s)
4
1 31
2 /(s+1) /(s+1)

6 7
5 1
C'(s)

32
Respuesta de Régimen Transitorio

Obtener las funciones de transferencia de cada uno de los bloques que componen el
sistema.

Regulador

Si se deriva la integral de la respuesta, se tendrá la respuesta, luego:

3 5t 5t 5t
c( t ) 3 e ( 5) 3 3e 3 (1 e )
5

Si se hace la transformada de Laplace de esta respuesta, se tiene:

1 0.2 1 0.2s 1 0.2s 1


C(s) 3 R (s) G (s) G (s) 3 3
s 0.2s 1 s s (0.2s 1) s (0.2s 1)

Despejando G(s):
3 15
G (s)
0.2s 1 s 5
Sensor

La respuesta es la de un sistema de segundo orden, que se puede expresar en su forma


canónica como:
K 2n
G (s)
s2 2 ns 2
n

De la gráfica de la respuesta, se tiene que:

- El tiempo de pico (tp) es 1.2 segundos.

2.296 2
- El sobreimpulso (Mp) es 0.148 14.8% .
2

A partir de estos valores, se pueden obtener los de n y , de acuerdo con las siguientes
relaciones:
2
1
Mp e 0.148 0.52

3.678
tp 1.226 n 3 rad seg
2 2
d n 1 n 1 0.52 n

luego sustituyendo en la ecuación de G(s), quedará:

K 32 18
G (s) 2
;K ganancia del sistema 2
s 2 0.52 3 s 32 s 2
3.12 s 9

33
Problemas de Ingeniería de Sistemas: Sistemas Continuos. Conceptos básicos.

Planta

La relación entre la salida C'(s) y la entrada R'(s), viene dada por el flujograma. La
obtención de la función de transferencia puede verse en el ejercicio 3.14.:

C' (s) 56 s 2 145 s 103


M ' (s)
R ' (s) (s 1) 2

Sustituyendo en el diagrama de bloques, tendremos:


W1(s)
E(s)
+ C(s)
R(s)
2
15 56 s 145 s 103
K s 3
+ s 5 + (s 1) 2
-
+ 18
2
s 3.12 s 9
+
W2(s)

EJERCICIO 4.17.

0.1
Se sabe que G0(s) ante una entrada R (s) tiene una respuesta:
s
respuesta
28

0
tiempo

y que G1(s) ante una entrada escalón de amplitud A = 2 tiene una respuesta:

respuesta

0 tiempo

34
Respuesta de Régimen Transitorio

de G2(s) se conoce su respuesta impulso unitario:

respuesta
3

Área=1

tiempo

de G3(s) se conoce también su respuesta impulso unitario:

Respuesta
0.3333

0.1326
tiempo

0.6

y que G4(s) ante una entrada rampa unitaria tiene como respuesta:

Calcular dichas funciones de transferencia.

Comenzando por G0(s):

La entrada es un escalón de valor 0.1, y según la gráfica de la respuesta, se ve que es un


escalón de amplitud 28, luego hemos amplificado la entrada un valor de:
28
G 0 (s ) G 0 (s) 280
0.1

35
Problemas de Ingeniería de Sistemas: Sistemas Continuos. Conceptos básicos.

Para G1(s), puede verse en la gráfica que la respuesta es un escalón de valor 2 y como esta
respuesta se produce ante una entrada escalón de amplitud 2, la función de transferencia G1(s)
es la unidad:
G 1 (s ) 1

De G2(s), se sabe por el tipo de respuesta que se trata de un sistema de primer orden, que se
define por la constante de tiempo y ganancia como:

k
G 2 (s )
Ts 1

k t
La ecuación de la respuesta es: c( t ) e T
, para t 0
T

Integrando entre 0 e , se obtiene el valor del área encerrada, y de ahí el valor de la


ganancia.
t
k t k t k e T k kT
Área 1 e T
dt e T
dt ( T) 0 1 k
0
T T0 T 1T T T
0
Luego la ganancia:
k=1
Del valor de la ordenada en el origen se obtiene la constante de tiempo, ya que:

Para t = 0:
k 0 1 1
c(0) 3 3 e T
T
T T 3
Luego será:
1
G 2 (s )
s
1
3

Para G3(s) la respuesta tiene por ecuación:


t
K
g 3 ( t) e T
T

y se dan dos puntos de dicha gráfica, por lo que se puede hallar las dos incógnitas que se
tienen en la ecuación:
K
Para t = 0 g 3 ( 0) 0.3333
T
K 0.6 0.6
Para t = 0.6 g 3 (0.6) 0.1326 e T 0.3333 e T
T

0.6
0.1326 0.6
ln e T
ln 0.9217 T 0.651
0.3333 T
Luego:
0.217
G 3 (s )
1 0.651s

36
Respuesta de Régimen Transitorio

G4(s):

En la gráfica vemos que se trata de la respuesta de un sistema de primer orden a una


entrada rampa unitaria, siendo el sistema de la forma:
1
G 4 (s )
Ts 1

siendo T el error en el infinito, luego de la gráfica:

0.08
T = 0.08 segundos
0.08

Luego:

1
G 4 (s )
0.08s 1

T 2TT 3T

EJERCICIO 4.18.

Obtener la señal de salida del sistema de control de la figura,


R2 (s)
+
R1 (s) 2 + 1 Y(s)
+_ s 3 s

cuando se le aplican las entradas r1(t) y r2(t).


r1 (t) r2 (t)
3

t
t -2

Aplicando superposición de las dos entradas se calculará la salida.


La salida cuando r2(t)=0 es:
R1 (s) 2 Y(s)
1
+_ s 3 s

37
Problemas de Ingeniería de Sistemas: Sistemas Continuos. Conceptos básicos.

R1 (s) 2 Y(s)

+_ s(s 3)

2
Y(s) s(s 3) 2
R 1 (s) 2 s 2
3s 2
1
s(s 3)

2 2 3 6
Y(s) R 2 R 1 (s) 2 2
2 ( s) 0
s 3s 2 s 3s 2 s s(s 3s 2)

y cuando r1(t)=0:
R2 (s) Y(s)
1
+ s
+

2
1
s 3

1
Y(s) s s 3
R 2 (s) 2 s 2
3s 2
1
s(s 3)
s 3 s 3 2 2s 6
Y(s) R 2 R 1 (s) 2 2
1 ( s) 0
s 3s 2 s 3s 2 s s(s 3s 2)

Aplicando superposición:

6 2s 6 2s 2
Y(s) Y(s) R Y(s) R 2 2 2 2
2 ( s) 0 1 ( s) 0 s(s 3s 2) s(s 3s 2 ) s(s 3s 2) s 3s 2

Se obtiene ahora la transformada inversa de la respuesta:

2 2 A B
Y(s) 2
s 3s 2 (s 1)(s 2 ) s 1 s 2
2 2
A 2
(s 2 ) s 1 1

2 2
B 2
(s 1) s 2
1

38
Respuesta de Régimen Transitorio

2 2
Y(s)
s 1 s 2
t 2t 2t
y (t ) 2e 2e 2 (e e t)

EJERCICIO 4.19.

Para el sistema del ejercicio 3.21. cuyo diagrama de bloques se muestra a continuación:
n(t) N(s)

r(t) +
c(t)
2
+ s 3
R(s) + C(s)
_

3
s 2

donde la entrada r(t) es un escalón de amplitud 2, y la perturbación n(t) es de tipo exponencial


decreciente, de amplitud unidad y constante de tiempo 1/3 segundos, calcular la respuesta
temporal en ausencia y presencia de ruido.

En el ejercicio 3.21. se tiene el cálculo de las funciones de transferencia.

La función de transferencia salida/referencia es:

C(s) 2( s 2 )
R (s ) s 2 5s 12

Y la función de transferencia salida/ruido:


C(s) s 2 5s 6
R (s) s 2 5s 12

Se calculará la transformada inversa de Laplace de la salida en cada una de las situaciones.

Sin ruido, (N(s) = 0):


1
L C( s ) c ( t)

4(s 2) A Bs C
C(s) 2 2
s(s 5s 12) s s 5s 12

En donde, A, B y C se determinan a partir de:

4s 8 A s 2 5 s 12 Bs C s

39
Problemas de Ingeniería de Sistemas: Sistemas Continuos. Conceptos básicos.

Si se comparan los coeficientes de s2, s y s0, términos de los dos miembros de la ecuación
anterior, se obtiene:

0 A B, 4 5A C, 8 12A

2 2
Resolviendo el sistema, se tiene: A C B
3 3
Por tanto,
21 2 s 2 1
C(s) 2 2
3s 3 5
2
23 3 5
2
23
s s
2 2 2 2

5 5
s 1
2 1 2 2
2 2
3 s 5
2
23 5
2
23
s s
2 2 2 2

5 23
s
2 1 2 7 2 2
2 2
3 s 5
2
23 2 23 5
2
23
s s
2 2 2 2

Luego,
5t 5t
2 23 7 23 23
c( t ) 1 e 2 cos t e 2 sen t
3 2 23 2

1
Con ruido, N (s)
s 3

1 1 1
L C R ( s ) C N (s ) L C R (s ) L C N (s ) c (t)

1 1 4 (s 2)
L C R (s ) L 2
ya resuelta anteriormente
s (s 5 s 12)

5 5
s 2
1 1 s 2 1 2 2
L C N (s ) L L
s2 5 s 12 5
2
23
2

s
2 2

40
Respuesta de Régimen Transitorio

5 23
s
1 2 1 5 2 2
L 2
L 2 2
5
2
23 2 23 5
2
23
s s
2 2 2 2

23
5t
23 1 23 2
e 2 cos t L 2
2 23 5
2
23
s
2 2

5t 5t
23 23 23
e 2 cos t e 2 sen t
2 23 2

Luego:
5t 5t 5t
2 2 23 14 23 23 23
c (t) e 2 cos t e 2 sen t e 2 cos t
3 3 2 3 23 2 2

5t
23 23
e 2 sen t
23 2

5t 5t
2 1 23 11 23 23
c( t ) e 2 cos t e 2 sen t
3 3 2 3 23 2

41

También podría gustarte