Dinamica de Sistemas S1o y S2o
Dinamica de Sistemas S1o y S2o
Dinamica de Sistemas S1o y S2o
C (s) 1
= ( EC 1)
R (s ) Ts+1
En lo sucesivo, se analizan las respuestas del sistema a entradas como la función escalón
unitario, rampa unitaria e impulso unitario. Se supone que las condiciones iniciales son cero.
Obsérvese que todos los sistemas que tienen la misma función de transferencia presentarán
la misma salida en respuesta a la misma entrada. Para cualquier sistema físico dado, la
respuesta matemática recibe una interpretación física.
1 1
C ( s )=
Ts+1 s
1 T 1 1
C ( s )= − = − ( EC 2)
s Ts+ 1 s s + ( 1/ T )
La Ecuación 3 plantea que la salida c(t) es inicialmente cero y al final se vuelve unitaria. Una
característica importante de tal curva de respuesta exponencial c(t) es que, para t =T, el
valor de c(t) es 0.632, o que la respuesta c(t) alcanzó 63.2% de su cambio total. Esto se
aprecia con facilidad sustituyendo t =T en c(t). Es decir,
| |
❑
dc 1 1
¿ e−t /T = ( EC 4)
dt t−0 T t−0 T
1 1
C ( s )=
Ts+1 s ²
1 T T²
C ( s )= − + (EC 5)
s ² s Ts+ 1
e ( t )=r ( t )−c ( t )
¿ T (1−e−t /T )
Conforme t tiende a infinito, e−t /T se aproxima a cero y, por tanto, la señal de error e(t) se
aproxima a T o
e(∞)=T
Para la entrada impulso unitario, R(s) =1 y la salida del sistema de la Figura 1(a) pueden
obtenerse como
1
c ( s )= (EC 7)
Ts+1
En el análisis anterior, se demostró que, para la entrada rampa unitaria, la salida c(t) es
Para la entrada escalón unitario, que es la derivada de la entrada rampa unitaria, la salida
c(t) es
−t /T
c (t )=1−e , para t ≥ 0 [ vease EC 3 ]
Por último, para la entrada impulso unitario, que es la derivada de la entrada escalón
unitario, la salida c(t) es
1 −t /T
c ( t )= e , para t ≥ 0 [ vease EC 8 ]
T
Una comparación de las respuestas del sistema para estas tres entradas indica con claridad
que la respuesta a la derivada de una señal de entrada se obtiene diferenciando la
respuesta del sistema para la señal original. También se observa que la respuesta para la
integral de la señal original se obtiene integrando la respuesta del sistema para la señal
original y determinando las constantes de integración a partir de la condición inicial de salida
cero. Esta es una propiedad de los sistemas lineales e invariantes con el tiempo. Los
sistemas lineales y variables con el tiempo y los sistemas no lineales no poseen esta
propiedad.
Jċ+Bċ=T
Js²C(s)+BsC(s) =T(s)
C (s) 1
=
T ( s ) s( Js +B)
Utilizando esta función transformada, la Figura 5-5(a) se puede redibujar como se muestra
en la Figura 5(b), que se puede modificar como se muestra en la Figura 5(c). La función de
transferencia en lazo cerrado se obtiene entonces como
C (s) K K /J
= =
R (s ) Js ²+ Bs+ K B K
( )
s ²+
J
s+( )
J
Tal sistema en el que la función de transferencia en lazo cerrado posee dos polos se
denomina sistema de segundo orden. (Algunos sistemas de segundo orden pueden
contener uno o dos ceros.)
Ilustración 5:(a) Servosistema; (b) diagrama de bloques; (c) diagrama de bloques
simplificado.
C (s) K
= ( EC 9)
R (s ) Js ²+ Bs+ K
C (s) K/J
=
R (s )
[ s+
B
2J
+ (√ 2BJ )²− KJ ][ s+ 2BJ −√( 2BJ ) ²− KJ ]
Los polos en lazo cerrado son complejos si B²-4JK <0, y son reales si B²-4JK≥ 0. En el
análisis de la respuesta transitoria, es conveniente escribir
K B
=ω ²n =¿ 2ζ ω n =2� donde � se denomina atenuación; ω n ,
J J
frecuencia natural no amortiguada, y ζ, factor de amortiguamiento relativo del sistema. El
factor de amortiguamiento relativo ζ es el cociente entre el amortiguamiento real B y el
amortiguamiento crítico Bc=2 √ JK o bien
B B
ζ= =
Bc 2 √ JK
Ilustración 6: Sistema de segundo orden.
C (s) ω ²n
= (EC 10)
R (s ) s ²+ 2ζ ωn s+ ω ² n
Esta forma se denomina forma estándar del sistema de segundo orden. El comportamiento
dinámico del sistema de segundo orden se describe a continuación en términos de dos
parámetros ζ y �n. Si 0< ζ <1, los polos en lazo cerrado son complejos conjugados y
se encuentran en el semiplano izquierdo del plano s. El sistema, entonces, se denomina
subamortiguado y la respuesta transitoria es oscilatoria. Si ζ =0, la respuesta transitoria
no se amortigua. Si ζ =1, el sistema se denomina críticamente amortiguado. Los
sistemas sobreamortiguados corresponden a ζ >1.
Ahora se obtendrá la respuesta del sistema que aparece en la Figura 6 para una entrada
escalón unitario. Se considerarán tres casos diferentes: el subamortiguado (0< ζ <1), el
críticamente amortiguado ( ζ =1) y el sobreamortiguado ( ζ > 1) .
C ( s) ω ²n
=
R (s ) (s+ζ ω n +J ωd )( s+ ζ ωn−J ω d )
ω2n
C ( s )= (EC 11)
( s 2+2 ζ ω n s+ ω2n ) s
1 s+ 2ω ²n
C ( s )= −
s s ²+2 ζ ωn s +ω ² n
1 s+ ω ² n ζ ωn
− −
s ( s+ζ ω n ) +ω ² d ( s +ζ ω n )2+ ω ² d
2
L−1=
[ s +ζ ω n
( s +ζ ω n ) +ω d
2 2
] =e−ζ ω t cos ω d t
n
−1
L =
[ ωd
( s +ζ ω n ) +ω d
2 2
] =e
−ζ ωn t
sen ω d t
ζ
¿ 1−e
−ζ ω n t
( cos ω d t+
√1−ζ 2
sen ω d t
) ,
e−ζ ω t −1 √ 1−ζ
2
( )
n
ζ
¿e
−ζ ωn t
( cos ωd t +
√1−ζ 2
sen ωd t
) , para t≥0
Esta señal de error presenta una oscilación sinusoidal amortiguada. En estado estacionario,
o en t=∞, no existe un error entre la entrada y la salida. Si el factor de amortiguamiento
relativo f es igual a cero, la respuesta se vuelve no amortiguada y las oscilaciones continúan
indefinidamente. La respuesta c(t) para el caso del amortiguamiento cero se obtiene
sustituyendo ζ =0 en la Ecuación 12, lo cual produce
2) Caso críticamente amortiguado ( ζ =1): si los dos polos de C(s)/R(s) son casi iguales,
el sistema se aproxima mediante uno críticamente amortiguado. Para una entrada escalón
unitario, R(s) =1/s y C(s) se escribe como
ω ²n
C ( s )= EC 14
(s +ω n)² s
3) Caso sobreamortiguado ( ζ =1 ): en este caso, los dos polos de C(s)/R(s) son reales
negativos y diferentes. Para una entrada escalón unitario, R(s)=1/s y C(s) se escriben como
ω ²n
C ( s )= EC 16
( s +ζ ω n+ ωn √ ζ 2−1 )( s +ζ ωn−ωn √ ζ 2 −1) s
La transformada inversa de Laplace de la Ecuación 16 es:
ζ + √ ζ 2−1
¿
2
ζ −1 ωn t
ζ + √¿
¿
−¿
2 √ζ 2−1¿
1
c (t )=1+ ¿
ζ −√ ζ 2−1
¿
ζ 2−1ω n t
ζ− √¿
¿
−¿
2 √ ζ 2−1 ¿
−1
¿
ωn es t es t
( )
1 2
¿ 1+ − , para t ≥0 EC 17
2 √ ζ 2−1 s 1 s2
Cuando ζ es apreciablemente mayor que la unidad, uno de los dos exponenciales que
decaen disminuye mucho más rápido que el otro, por lo que el término exponencial que
decae más rápido puede pasarse por alto (corresponde a una constante de tiempo más
pequeña). Es decir, si −s 2 se localiza mucho más cerca del eje J ω que −s 1 (lo cual
significa que |s2|≪|s 1|¿ para una solución aproximada se puede no considerar −s 1 .
Esto se permite debido a que el efecto de −s 1 e n la respuesta es mucho más pequeño
que el de −s 2 , ya que el término que incluye s 1 en la Ecuación 17 se descompone
mucho más rápido que el término que tiene a s 2 Una vez desaparecido el término
exponencial que decae más rápido, la respuesta es similar a la de un sistema de primer
orden, y C(s)/R(s) se aproxima mediante:
Esta forma aproximada es una consecuencia directa de que los valores iniciales y los
valores finales tanto del C(s)/R(s) original como del aproximado coincidan. Con la función de
transferencia aproximada C(s)/R(s), la respuesta escalón unitario se obtiene como
ζ ω n−ω n √ ζ ²−1
C ( s )=
( s+ ζ ωn −ωn √ ζ ²−1 ) s
La respuesta del tiempo c(t) es, entonces
2
ζ −1
ζ− √¿
¿
ωn t
−¿
c (t )=1−e ¿
Esto proporciona una respuesta escalón unitario aproximada cuando uno de los polos de
C(s)/ R(s) puede pasarse por alto. La Figura 7 contiene una familia de curvas c(t) con
diversos valores de ζ , donde la abscisa es la variable adimensional ω n t . Las curvas
sólo son funciones de ζ y se obtienen a partir de
1. Tiempo de retardo, td
2. Tiempo de subida, tr
3. Tiempo pico, tp
4. Sobreelongación, m p
5. Tiempo de asentamiento, ts
3. Tiempo pico, t p : el tiempo pico es el tiempo requerido para que la respuesta alcance el
primer pico de sobreelongación.
C ( t p )−c ( ∞)
Porcentajede sobreelongación máxima= × 100
c (∞)
deben presentar respuestas de tiempo aceptables. (Esto significa que el sistema de control
debe modificarse hasta que la respuesta transitoria sea satisfactoria.) Obsérvese que todas
estas especificaciones no se aplican necesariamente a cualquier caso determinado. Por
ejemplo, para un sistema sobreamortiguado no se aplican los términos tiempo pico y
sobreelongación máxima. (En los sistemas que producen errores en estado estacionario
para entradas escalón, este error debe conservarse dentro de un nivel de porcentaje
especificado. En la Sección 9 se incluyen análisis detallados de los errores en estado
estacionario.)
c( t r ) =1, o que
t
ζ
(¿¿ r )=1=1−e−ζ ω t cos ω d t r + n r
( √ 1−ζ ²
sen ω d t r ) EC 18
c¿
ζ
cos ω d t r + sen ω d t r =0
√1−ζ ²
Como ωn = √ 1−ζ ² = ω d y ζ ω n=σ , se tiene
1−ζ ² ω d
tan ω d t r= √ =
ζ σ
ω
t r=
1
ωd
tan−1 d =
−σ ωd( )
π −β
( EC 19)
Tiempo pico t p : si se remite a la Ecuación 12, se obtiene el tiempo pico diferenciando c(t)
con respecto al tiempo y suponiendo que esta derivada es igual a cero. Como
dc ζ
dt
=ζ ωn e
−ζ ω t
cos ωd t +n
(
√1−ζ ²
sen ω d t
)
ζ ωd
+e
−ζ ω n t
( ω d sen ωd −
√ 1−ζ ²
cosω d t )
y los términos de coseno de esta última ecuación se cancelan uno al otro, dc/dt, evaluada en
t= t p , se simplifica a
dc ωd −ζ ω t
=( sen ωd t p ) e =0 n p
dt t =t p √1−ζ 2
sen ωd t p=0
o bien
ω d t p =0, π , 2 π , 3 π , …
π
t p= (EC 20)
ωd
M p=c ( t p )−1
ζ
¿−e−ζ ω (π / ω ) cos π +
n d
( √ 1−ζ ²
sen π )
¿ e−(σ / ω ) π=e−(ζ / √ 1−ζ )π (EC 21)
2
d
c ( t p ) −c (∞ )
M p=
c (∞)
e−ζ ω t −1 √ 1−ζ
2
( )
n
Las curvas 1 ± ( e−ζ ω t / √ 1−ζ 2 ) son las curvas envolventes de la respuesta transitoria para
n
una entrada escalón unitario. La curva de respuesta c(t) siempre permanece dentro de un
par de curvas envolventes, como se aprecia en la Figura 10. La constante de tiempo de
estas curvas envolventes es 1/ ζ ω n .
La velocidad de decaimiento de la respuesta transitoria depende del valor de la constante de
tiempo 1/ ζ ω n . Para un ω n determinado, el tiempo de asentamiento t s es una función
del factor de
Los resultados se muestran en la Figura 11. Para 0 ¿ ζ <¿ 0.9, si se utiliza el criterio del
2%, t s es aproximadamente cuatro veces la constante de tiempo del sistema. Si se
emplea el criterio del 5%, t s es aproximadamente tres veces la constante de tiempo.
Obsérvese que el tiempo de asentamiento alcanza un valor mínimo alrededor de ζ =0.76
(para el criterio del 2%) o de ζ =0.68 (para el criterio del 5%), y después aumenta casi
linealmente para valores grandes de ζ . Las discontinuidades en las curvas de la Figura
11 surgen debido a que un cambio infinitesimal en el valor de ζ puede provocar un
cambio finito en el tiempo de asentamiento.
Por conveniencia, cuando se comparan las respuestas de los sistemas, por lo general se
define el tiempo de asentamiento t s como
4 4
t s=4 T = = ( criterio del2 ) ( EC 22)
σ ζ ωn
3 3
t s=3 T = = ( criterio del 5 ) (EC 23)
σ ζ ωn
A partir del análisis anterior, es evidente que, para una respuesta rápida, ω n debe ser
grande. Para limitar la sobreelongación máxima M p y para reducir el tiempo de
asentamiento, el factor de amortiguamiento relativo ζ no debe ser demasiado pequeño.
La relación entre la sobreelongación en porcentaje M p y el factor de amortiguamiento
relativo ζ se presenta en la Figura 12. Obsérvese que, si el factor de amortiguamiento
relativo está entre 0.4 y 0.7, el porcentaje de sobreelongación máxima para la respuesta
escalón está entre 25 y 4%.
Ilustración 12: M p frente a la curva ζ.
Es importante darse cuenta de que las ecuaciones que se han obtenido para el tiempo de
subida, tiempo de pico, sobreelongación máxima y tiempo de asentamiento son válidas
únicamente para el sistema de segundo orden estándar definido por la Ecuación 10. Si el
sistema de segundo orden contiene uno o dos ceros, la forma de la curva de respuesta a un
escalón unitario será bastante diferente de las que se muestran en la Figura 7.
C (s) k
= (EC 24)
R (s ) Js ²+ ( B+ k k h ) s+k
B+k k h
ζ= ( EC 25)
2 √ kJ
Respuesta impulso de sistemas de segundo orden. Para una entrada impulso unitario
r(t), la transformada de Laplace correspondiente es la unidad, o R(s) =1. La respuesta
impulso unitario C(s) del sistema de segundo orden de la Figura 6 es
ω ²n
c ( s )=
s ²+2 ζ ωn s +ω ² n
Para 0≤ ζ <1,
ωn
c ( t )= e−ζ ω t sen ω n √ 1−ζ 2 t , para t ≥ 0(EC 26)
n
√ 1−ζ 2
Para ζ =1 ,
−ω n t
c (t )=ω ²n t e , para t ≥0 ( EC 27)
Para ζ >1 ,
ωn −(ζ − √ζ −1)ωn t
2 ωn −(ζ + √ ζ −1) ωn t
2
Obsérvese que, sin tomar la transformada inversa de Laplace de C(s), también se obtiene el
tiempo de respuesta c(t) diferenciando la respuesta escalón unitario correspondiente, ya que
la función impulso unitario es la derivada con respecto al tiempo de la función de escalón
unitario. En la Figura 14 aparece una familia de curvas de respuesta impulso unitario
obtenida mediante las Ecuaciones 26 y 27 con diversos valores de ζ . Las curvas c(t)/
ω n se dibujan frente a la variable adimensional ω n t y, por tanto, sólo son funciones de
ζ . Para los casos críticamente amortiguado y sobreamortiguado, la respuesta impulso
unitario siempre es positiva o cero; es decir, c(t) ≥0. Esto se aprecia en las Ecuaciones 27 y
28. Para el caso subamortiguado, la respuesta impulso unitario c(t) oscila alrededor de cero
y toma valores tanto positivos como negativos.
tan −1 √
1−ζ 2
ζ
t= , donde 0<ζ <1( EC 29)
ω n √ 1−ζ 2
[La Ecuación 29 se puede obtener igualando dc/dt a cero y resolviendo para t.] La
sobreelongación máxima es
−1 √ 1−ζ
2
c (t)max =ω n exp
( −ζ
√1−ζ 2
tan
ζ ) , donde 0< ζ <1(EC 30)