Cap 1 Eloy
Cap 1 Eloy
Cap 1 Eloy
Año sabático 2016-2017. AUTOR: DR. JUAN ELOY BENITO HERNÁNDEZ Página i
AGRADECIMIENTOS.
Este libro de texto para Señales y Sistemas de control clásico con soporte de
Espacio de Estado fue preparado para los estudiantes del sexto semestre de la
Academia de Control de la carrera de Ingeniería en comunicaciones y Electrónica
dependiente de la Escuela Superior de Ingeniería Mecánica y Eléctrica- Zacatenco,
gracias al apoyo del Año Sabático que me concedió el Instituto Politécnico Nacional
durante el periodo 2016-2017.
Año sabático 2016-2017. AUTOR: DR. JUAN ELOY BENITO HERNÁNDEZ Página ii
PREFACIO.
Año sabático 2016-2017. AUTOR: DR. JUAN ELOY BENITO HERNÁNDEZ Página iii
INTRODUCCIÓN
En el mundo actual existen una infinidad de sistemas que requieren ser controlados:
una torre de destilación, un reactor químico, intercambiadores de calor, hornos de
fundición, un radar, el timón de una embarcación, un elevador, una impresora y un
número tan grande de sistemas que sería imposible enumerarlos aquí.
𝑞(𝑡) = 𝑘𝑣 𝑥(𝑡)
Donde 𝑞(𝑡) gasto
𝑥(𝑡) posición de apertura de la válvula
𝑘𝑣 ganancia de la válvula
Año sabático 2016-2017. AUTOR: DR. JUAN ELOY BENITO HERNÁNDEZ Página iv
φ
i excitación
Año sabático 2016-2017. AUTOR: DR. JUAN ELOY BENITO HERNÁNDEZ Página v
Los sistemas no lineales son representados por ecuaciones diferenciales no
lineales, cuando esto es posible, en otras ocasiones se representan con modelos
físicos o gráficos.
Un ejemplo clásico de una ecuación deferencial no lineal es la ecuación de Van der
Pool
𝑦̈ + 𝜖 (𝑦 2 − 1)𝑦̇ + 𝑦 = 0
Donde 𝜖 es constante y negativa.
Los sistemas pueden ser sumamente complejos por estar constituidos de
subsistemas relacionados entre sí. En el capítulo uno de este este libro
estudiaremos diferentes tipos de sistemas de acuerdo a su naturaleza física:
modelos mecánicos, eléctricos, electromecánicos, hidráulicos, neumáticos,
térmicos y otros especiales que consideramos de interés.
Año sabático 2016-2017. AUTOR: DR. JUAN ELOY BENITO HERNÁNDEZ Página vi
INDICE
AGRADECIMIENTOS ii
PREFACIO iii
INTRODUCCION iv
CAPÍTULO UNO. 1
MODELADO DE SISTEMAS.
CAPITULO DOS.
Año sabático 2016-2017. AUTOR: DR. JUAN ELOY BENITO HERNÁNDEZ Página vii
Orden a una entrada escalón. 207
Estados. 235
CAPITULO TRES.
CAPITULO CUATRO.
Año sabático 2016-2017. AUTOR: DR. JUAN ELOY BENITO HERNÁNDEZ Página viii
4.3 Margen de ganancia, margen de fase mediante el análisis grafico
de Bode y Nyquist. 304
CAPITULO CINCO.
BIBLIOGRAFÍA: 356
Año sabático 2016-2017. AUTOR: DR. JUAN ELOY BENITO HERNÁNDEZ Página ix
CAPITULO UNO
MODELADO DE SISTEMAS.
A mediados del siglo XIX y en pleno apogeo de la Revolución Industrial James Watt
desarrolla un regulador de velocidad de una máquina de vapor. Las máquinas de
vapor se emplearon como amplificadores de potencia y se les encontraron múltiples
usos. Maxwell en el año de 1868 desarrolla el primer tratamiento matemático de un
mecanismo de control con retroalimentación. En 1926 en Alemania, Stein escribe
sobre sistemas retroalimentados en la generación de vapor.
Durante la Segunda Guerra Mundial se invento uno de los sistemas clásicos del
control, el radar. Durante las dos décadas precedentes, a eta época se dio origen a
las principales ideas del control clásico reforzadas por el desarrollo de la aviación,
la electrónica, la teoría de circuitos, etc. Nyquist desarrolla un trabajo sobre
estabilidad de sistemas retroalimentados lineales, motivado primordialmente por el
análisis de redes de comunicación. Hazen desarrolla el primer trabajo enfocado
directamente a la teoría de los servomecanismos. La palabra servo fue usada por
primera vez (servus del latín esclavo o sirviente).
Año sabático 2016-2017. AUTOR: DR. JUAN ELOY BENITO HERNÁNDEZ Página 1
de regulación con la sensitividad, estabilidad y velocidad de acción. Mason realiza
un trabajo pionero en el establecimiento de los modelos y analogías de los procesos.
Después de 1945 el control queda constituido como una ciencia. Impulsado por la
presión de la guerra el MIT Radiation Laboratory participó significativamente en este
desarrollo. Después de la guerra Taplin y Harris establecen el concepto de función
de transferencia en estado estacionario (frecuencia). Garner y Barnes popularizan
en 1942 el uso del cálculo con transformadas de Fourier y Laplace. Muchas de las
ideas se toman de las aportaciones de algunos ingenieros en comunicaciones como
Bode.
Se abre la era aeroespacial (1956). La guerra fría y la conquista del espacio hacen
que Rusia y los Estados Unidos de Norteamérica inicien lo que se ha dado en llamar
el control moderno (control óptimo, control adaptivo, sistemas de autoaprendizaje,
control numérico). El genio matemático Norberto Winer introduce la cibernética
(kybernet del griego timonel o gobernar).
Año sabático 2016-2017. AUTOR: DR. JUAN ELOY BENITO HERNÁNDEZ Página 2
Sistema de Alrededores. El concepto de sistema está relacionado con el campo
en el cual concentramos nuestro interés de estudio y análisis. Dentro de él hay un
conjunto de entes relacionados por vínculos más o menos complejos y cuyo
comportamiento deseamos, en el mejor de los casos, predecir y/o controlar. Cuando
esto es posible, nosotros introducimos excitaciones al sistema para poder observar
su respuesta, el problema esta, en que la mayoría de las veces existen otras
entradas o perturbaciones al sistema que complican nuestro análisis. Nuestro poder
de predicción esta en la posibilidad de crear un modelo del sistema. Entre más se
ajuste el comportamiento de nuestro modelo al comportamiento del sistema,
diremos que nuestro conocimiento de él es mayor. El desarrollo de la ciencia está
basado en el mejoramiento de estos modelos de los sistemas que se presentan en
la naturaleza, basta con mencionar el modelo atómico que evoluciono de un objeto
sólido e indivisible a una compleja estructura cuyos últimos constituyentes son los
quarks.
En ingeniería nuestros modelos pueden ser de diferentes tipos: físicos los cuales no
están restringidos a la mera modificación de las dimensiones espaciales, su
comportamiento debe ser lo más apegado al sistema bajo las condiciones en las
cuales este opera o va a funcionar, como ejemplo tenemos los modelos
aerodinámicos de las alas de un avión en un túnel de viento o el modelo de una
presa en un sistema hidráulico; gráficos los cuales nos permiten obtener relaciones
de transferencia entre variables del sistema, por ejemplo las características de
saturación de materiales ferromagnéticos que relacionan la corriente de
magnetización y el flujo o las características de un diodo que relacionan el voltaje
de ánodo-cátodo con la corriente de cátodo. El modelo más empleado en el campo
de la ingeniería son los modelos matemáticos, los cuales consisten en escribir una
ecuación diferencial lineal o no lineal y cuyo orden n esta relacionado con la
complejidad del sistema.
Los sistemas establecen una relación de entradas y salidas las cuales podemos
representar en forma de bloques como ilustra la figura
Año sabático 2016-2017. AUTOR: DR. JUAN ELOY BENITO HERNÁNDEZ Página 3
y = S [U].
Y es la señal de salida.
1.2.1 clasificación.
Año sabático 2016-2017. AUTOR: DR. JUAN ELOY BENITO HERNÁNDEZ Página 4
Esta clasificación corresponde a la manera en que podemos considerar los
parámetros del sistema en relación a sus propiedades espaciales, por ejemplo un
sistema como el circuito RC, para su análisis consideramos las propiedades de
resistividad y capacitancia concentradas en un punto; en tanto que en una línea de
transmisión estas propiedades se deben considerar distribuidas a lo largo de toda
la línea, lo que modifica totalmente el análisis.
U(t) Y(t)
S
U(n) Y(n)
S
U[n]
U(t)
T(seg) n
Año sabático 2016-2017. AUTOR: DR. JUAN ELOY BENITO HERNÁNDEZ Página 5
Un sistema es invariante en el tiempo, si un corrimiento en t en la señal de entrada
causa un corrimiento en t en la señal de salida. Específicamente si y (t) ó y [n] es la
salida de un sistema invariante (continuo o discreto), cuando u (t) ó u [n] es la
entrada; entonces y (t-to) ó y [n-no] es la salida cuando u (t-to) ó u [n-no] se aplica.
Como ejemplos de sistemas invariantes en el tiempo tenemos el circuito RC,
considerando que sus parámetros de resistencia y capacitancia no varíen con
respecto al tiempo y las mismas entradas defasadas respecto al tiempo producirán
iguales salidas defasadas en ese mismo tiempo. El mismo sistema sirve de ejemplo
para ilustrar un sistema variante en el tiempo, ti tomamos en cuenta que las
condiciones de temperatura afectan a estos parámetros. Los cohetes espaciales
consumen la mayor parte de su combustible al ser lanzados, por lo que su
parámetro masa varía significativamente con el tiempo, lo que lo convierte en un
sistema variante en el tiempo.
u2 (t) = u1 (t - to)
y1 (n) = n u1 [n]
Si u2 [n] = u1 [n - no]
y2 [n] = n u1 [n - no]
Año sabático 2016-2017. AUTOR: DR. JUAN ELOY BENITO HERNÁNDEZ Página 6
Los sistemas estáticos o sin memoria son llamados así, ya que su salida en el
tiempo t depende de la entrada en ese mismo tiempo t. Una simple resistencia es
un ejemplo de un sistema estático, la corriente depende únicamente del valor
presente del voltaje y no de valores pasados o futuros. Los sistemas dinámicos o
con memoria son aquellos cuya salida en el tiempo t depende de la entrada en el
tiempo t y de otros valores de tiempo. Un ejemplo de este tipo de sistema es la
resistencia acompañada de un capacitor, el voltaje de salida en el capacitor,
dependerá de valores presentes del voltaje de entrada y de otros valores anteriores
de carga del capacitor. Los sistemas dinámicos se caracterizan por tener un
elemento de almacenamiento de energía, en sistemas eléctricos son bobinas y
capacitores, en sistemas mecánicos son inercias y resortes, etc. los que permiten
guardas la información.
1.2.1.7 Invertibles.
Año sabático 2016-2017. AUTOR: DR. JUAN ELOY BENITO HERNÁNDEZ Página 7
En un sistema anticipativo para un instante t2, el valor de la salida y (t2) depende de
algún valor de u en el intervalo [t2, oo].
0
0 4 8 12 t(seg)
T(seg) 0 2 4 6
y (t) =u (2 t) y (0) = u (0) = 1 y (2) = u (4) =1 ó .5 y (4) = u (8) =.5 y (6) = u (12) =1
y(t)
1
.5
0 2 4 6 t(seg)
t (seg) 0 8 16 24
y(t) =u (2 t) y (0) = u (0) = 1 y (2) = u (4) = 1 ó .5 y (4) = u (8) = .5 y (6) = u (12)= 1
y (t)
Año sabático 2016-2017. AUTOR: DR. JUAN ELOY BENITO HERNÁNDEZ Página 8
1
.5
0
8 16 24 t(seg)
T (seg) 0 8 16 24
Si u (t ) B para toda t
y(t ) = h ( )u(t − )d
y (t ) h( ) u (t − )d
y (t ) B h( )d
Año sabático 2016-2017. AUTOR: DR. JUAN ELOY BENITO HERNÁNDEZ Página 9
Un sistema discreto es estable, si es absolutamente sumable; un sistema continuo
es estable si es absolutamente integrable. Existen diferentes formas de visualizar la
estabilidad de un sistema, por ejemplo si se deja evolucionar el sistema en una
vecindad de un punto de equilibrio y el sistema regresa a ese punto, el sistema es
estable. Si por el contrario el sistema se aleja de él, el sistema es inestable. Existe
un caso crítico y es cuando el sistema se queda oscilando entorno al punto de
equilibrio. Un péndulo con fricción es un ejemplo de un sistema estable ya que al
apartar al péndulo de su posición de equilibrio, este regresa nuevamente a él; si
consideramos un péndulo ideal sin fricción, este oscilaría eternamente entorno al
punto de equilibrio.
Año sabático 2016-2017. AUTOR: DR. JUAN ELOY BENITO HERNÁNDEZ Página 10
U 1 (t)+ U 1 (t) y 1 (t)+ y 2 (t)
Ku(t) Ky(t)
S
y n (t) + a1 y n -1
(t) + … + a n-1 y 1 (t) + a n y (t) = b0 un (t) + b1 un – 1 (t) + … + b n-1 u1
(t) + b n u (t)
Donde los coeficientes a 1,…, a n, b0, b1, …, b n son constantes o cuando más función
del tiempo.
Año sabático 2016-2017. AUTOR: DR. JUAN ELOY BENITO HERNÁNDEZ Página 11
1.- Para los siguientes sistemas indicar si son invariantes o variantes en el tiempo,
continuos o discretos, causales o no causales y lineales o no lineales.
2.- Para los siguientes sistemas indicar si son continuos o discretos, estáticos o
dinámicos e invertibles o no invertibles.
a) y (t) = 2 u (t)
b) y (t) = u2 (t)
Su sistema inverso es z (t) = () = +u (t)
c) y [n] = () [k]
y [n - 1] = … + u [n – 2] + u [n – 1]
Año sabático 2016-2017. AUTOR: DR. JUAN ELOY BENITO HERNÁNDEZ Página 12
Restando y [n] – y [n – 1] = u [n]
Por lo que z [n] = {y [n] – y [n – 1]} y el sistema es dinámico, discreto e invertible.
U[n] +
- y [n]
Elevador al cuadrado
4.- Si y (t) = x (0.1 t + 1) y la entrada esta dada por la figura siguiente, indique si el
sistema es causal o no causal y grafique y (t) para esta entrada.
U(t)
20
10
0
2 4 t (seg)
t(seg) -10 10 30
Y(t) y(-10)=x(0)=20 ó 0 y (10)= x(2)=20y(30)=x(4)=10
Y(t)
Año sabático 2016-2017. AUTOR: DR. JUAN ELOY BENITO HERNÁNDEZ Página 13
20
10
0
10 0 10 30 t(seg)
V C Vc
V 1
Vc =
R + 1 / Cs Cs
V
Vc(s)=
RCs + 1
1
Si v (t)=(t) entonces V(s)= 1 y Vc (s)=
RCs + 1
Tomando transformada inversa vc (t)=1/RC e −1 / RC u(t)
Año sabático 2016-2017. AUTOR: DR. JUAN ELOY BENITO HERNÁNDEZ Página 14
Obtenga la respuesta del sistema si la entrada es un escalón unitario empleando la
integral de convolución.
Vc (t) = e − ( t − ) d
Vc (t)=1-e − t
Probar la estabilidad del sistema, dando una condición inicial y dejando que
evolucione a un punto de equilibrio.
dq(t )
Ri (t ) = v(t ) R = v(t )
dt
Año sabático 2016-2017. AUTOR: DR. JUAN ELOY BENITO HERNÁNDEZ Página 15
1 1
C
i (t )dt = v(t ) q(t ) = v(t )
C
di(t ) d 2 q(t )
L = v(t ) L =v(t)
dt dt 2
Los problemas que vamos a considerar se pueden resolver mediante las leyes de
Kirchhoff:
Ley de voltajes que establece que la suma algebraica de todos los voltajes de
las ramas alrededor de cualquier malla de un circuito es cero en todos los instantes
de tiempo.
Ley de corrientes que establece que la suma algebraica de todas las corrientes de
las ramas que concurren en un nodo es cero para todos los instantes de tiempo.
Año sabático 2016-2017. AUTOR: DR. JUAN ELOY BENITO HERNÁNDEZ Página 16
Los parámetros básicos de un sistema mecánico son: la inercia la cual es
determinada por la masa de un cuerpo, también llamada por esto masa inercial y
que mide la tendencia de un cuerpo a oponerse a cambios en su velocidad, el
coeficiente de fricción viscosa que mide la oposición a moverse en un medio fluido
como el aire, aceite, etcétera y la compliancia que mide la oposición a la
deformación de los cuerpos elásticos como resortes, ligas, etcétera. Estos
parámetros pueden ser descritos mediante ecuaciones diferenciales lineales. A
continuación se dan estas ecuaciones considerando como excitación a la fuerza,
como respuesta a la velocidad y como variable independiente el tiempo.
1. Masa o Inercia.
du (t )
Su ecuación es f (t ) = m
dt
d 2 x(t )
f (t ) = m
dt 2
Año sabático 2016-2017. AUTOR: DR. JUAN ELOY BENITO HERNÁNDEZ Página 17
Ns
Donde D coeficiente de fricción viscoso
m
Su ecuación es f (t ) = Du (t )
dx(t )
f (t ) = D
dt
Compliancia.
m
Donde K compliancia elástica
N
1
Su ecuación es f (t ) = u (t )dt
K
1
f (t ) = x(t )
K
Para el caso rotacional el par es la variable de excitación, la velocidad angular y la
posición angular son las respuestas y el tiempo es la variable independiente.
Año sabático 2016-2017. AUTOR: DR. JUAN ELOY BENITO HERNÁNDEZ Página 18
1.3.1.3. Sistemas Rotacionales.
d 2 (t )
(t ) = J
dt 2
Su ecuación es (t ) = D (t )
d (t )
(t ) = D
dt
Compliancia angular.
Año sabático 2016-2017. AUTOR: DR. JUAN ELOY BENITO HERNÁNDEZ Página 19
1
Donde K compliancia angular
Nm
1
K
Su ecuación es (t ) = (t )dt
1
(t ) = (t )
K
Aplicando el llamado principio de D’Alambert que establece que “la suma de las
fuerzas externas aplicadas a un cuerpo es igual a la suma de las fuerzas resistivas
Año sabático 2016-2017. AUTOR: DR. JUAN ELOY BENITO HERNÁNDEZ Página 20
al movimiento. Entendiendo por fuerzas resistivas las que presentan las inercias,
fricciones y compliancias.
du (t ) 1
f (t ) = m + Du (t ) + u (t )dt
dt K
1
f ( s) = msu( s ) + Du ( s) + u ( s)
Ks
1
s
u ( s) m
=
f ( s) D 1
s2 + s +
m Km
x( s )
Para obtener la función de transferencia , recordamos que el
f ( s)
desplazamiento es la integral de la velocidad e integrar en el dominio de Laplace es
dividir entre s .
1
x( s ) m
=
f ( s) D 1
s + s+
2
m Km
Año sabático 2016-2017. AUTOR: DR. JUAN ELOY BENITO HERNÁNDEZ Página 21
En este caso las dos masas representan dos puntos rígidos y por lo tanto
escribimos una ecuación por cada uno.
du1 (t ) 1
K1
f (t ) = m1 + D1 (u1 (t ) − u 2 (t )) + (u1 (t ) − u 2 (t ))dt
dt
du 2 (t ) 1 1
0 = m2
dt
+ D1 (u 2 (t ) − u1 (t )) +
K1 (u 2 (t ) − u1 (t ))dt + D2 u 2 (t ) +
K2
u 2 (t )dt
1 1
f ( s) = (m1 s + D1 + )u1 ( s) − ( D1 + )u 2 ( s)
K1 s K1 s
1 1 1
0 = −( D1 + )u1 ( s ) + (m2 s + D1 + D2 + + )u 2 ( s )
K1 s K1 s K 2 s
1
D1 +
K1 s
Despejando de la segunda ecuación u 2 ( s) = u (s)
1 1 1
m2 s + D1 + D2 + +
K1 s K 2 s
u1 ( s )
Sustituyendo en la primer ecuación y despejando a
f (s)
Año sabático 2016-2017. AUTOR: DR. JUAN ELOY BENITO HERNÁNDEZ Página 22
1 1
m2 s + D1 + D2 + +
u1 ( s) K1 s K 2 s
=
f ( s) 1 1 1 1 2
(m1 s + D1 + )( m2 s + D1 + D2 + + ) − ( D1 + )
K1 s K1 s K 2 s K1 s
Para
1
D1 +
u 2 ( s) K1 s
=
f (s) 1 1 1 1 2
(m1 s + D1 + )( m2 s + D1 + D2 + + ) − ( D1 + )
K1 s K1 s K 2 s K1 s
1 D + D2 1 1
s( s 2 + 1 s+( + ))
u1 ( s) m1 m2 K 1 m2 K 2 m2
=
f ( s)
Donde
D1 D2 D1 3 1 1 1 D + D2 2 1 D D 1
= s4 + ( + + )s + ( + + + 1 )s + ( 1 + 2 )s +
m2 m2 m1 m2 K1 m2 K 2 m1 K1 m1m2 m1m2 K 2 K1 K1 K 2 m1m2
Y
1 1
s( D1 s + )
u 2 ( s) m1m2 K1
=
f ( s)
Año sabático 2016-2017. AUTOR: DR. JUAN ELOY BENITO HERNÁNDEZ Página 23
En este caso tenemos nuevamente dos puntos rígidos, por lo que escribimos dos
ecuaciones
d1 (t ) 1
K 2
(t ) = J 1 + (1 (t ) − 2 (t ))dt
dt
d 2 (t ) 1
K 2
0 = J 2 + ( 2 (t ) − 1 (t ))dt + D 2 2 (t )
dt
1 1
( s) = ( J 1 s + ) 1 ( s ) − 2 ( s)
K 2 s K 2 s
1 1
0=− 1 ( s ) + ( J 2 s + D 2 + ) 2 ( s )
K 2 s K 2 s
1
K 2 s
2 ( s) = 1 ( s )
1
J 2 s + D 2 +
K 2 s
Año sabático 2016-2017. AUTOR: DR. JUAN ELOY BENITO HERNÁNDEZ Página 24
Sustituyendo en la primera de las ecuaciones
1
( )2
1 K 2 s
( s) = ( J 1 s + )1 ( s ) − 1 ( s )
K 2 s 1
J 2 s + D 2 +
K 2 s
1
J 2 s + D 2 +
1 ( s ) K 2 s
=
( s) ( J1s +
1
)( J 2 s + D 2 +
1
)−(
1 2
)
K 2 s K 2 s K 2 s
Simplificando
1 2 D 2 1
s + s+
1 ( s ) J1 J1 J 2 K 2 J 1 J 2
=
( s) D 1 1 D 2
s3 + 2 s2 + ( + )s +
J 2 J 2 K 2 J 1 K 2 K 2 J 1 J 2
La otra transferencia es
1
2 ( s) K 2 J 1 J 2
=
( s) D 2 1 1 D 2
s3 + s2 + ( + )s +
J 2 J 2 K 2 J 1 K 2 K 2 J 1 J 2
Año sabático 2016-2017. AUTOR: DR. JUAN ELOY BENITO HERNÁNDEZ Página 25
1.3.1.4. Sistemas de acoplamiento mecánico.
Sistema de balancín.
Considerando que los desplazamientos en los extremos del balancín son muy
pequeños, se debe cumplir
x1 x2
=
a b
Año sabático 2016-2017. AUTOR: DR. JUAN ELOY BENITO HERNÁNDEZ Página 26
x1 x 2
t = t ó
u1 (t ) u 2 (t )
=
a b a b
u1 (t ) a
=
u 2 (t ) b
f1 (t )u1 (t ) = f 2 (t )u2 (t )
Es claro que la ecuación anterior indica que cualquier variación de fuerza implica
una variación inversa de la velocidad. Generalmente suelen emplearse como
amplificadores de fuerza.
f1 (t ) b
=
f 2 (t ) a
1 (t )d1 = 2 (t )d 2
Año sabático 2016-2017. AUTOR: DR. JUAN ELOY BENITO HERNÁNDEZ Página 27
Donde 1 (t ) y 2 (t ) son las velocidades angulares de las ruedas y d1 y d 2 son
sus diámetros.
1 (t )1 (t ) = 2 (t )2 (t )
1 (t ) d1
=
2 (t ) d 2
1 (t ) n1
=
2 (t ) n2
1 (t ) n2
=
2 (t ) n1
Año sabático 2016-2017. AUTOR: DR. JUAN ELOY BENITO HERNÁNDEZ Página 28
D1
du1 (t ) 1
K1
f (t ) = m1 + D1u1 (t ) + u1 (t )dt + f1 (t )
dt
du 2 (t ) 1
K2
f 2 (t ) = m2 + D2 u 2 (t ) + u 2 (t )dt
dt
f1 (t ) b u1 (t ) a
= =
f 2 (t ) a u 2 (t ) b
1
f ( s ) = (m1 s + D1 + )u1 ( s ) + f1 ( s )
K1 s
1
f 2 ( s ) = ( m 2 s + D2 + )u 2 ( s )
K2s
Eliminando f1 (t ) , f 2 (t ) y u 2 (t ) mediante
Año sabático 2016-2017. AUTOR: DR. JUAN ELOY BENITO HERNÁNDEZ Página 29
f1 (t ) b u1 (t ) a
= =
f 2 (t ) a u 2 (t ) b
1 b 1
f ( s) = (m1 s + D1 + )u1 ( s) + ( ) 2 (m2 s + D2 + )u1 ( s)
K1 s a K2s
b b 1 b 1
f ( s) = [( m1 + ( ) 2 m2 ) s + ( D1 + ( ) 2 D2 ) + ( + ( )2 )]u1 ( s)
a a K1 s a K 2 s
1
f ( s) = [meq s + Deq + ]u1 ( s)
K eq s
1
s
u1 ( s ) meq
=
f ( s) Deq 1
s2 + s+
meq meq k eq
b
Donde meq = m1 + ( ) 2 m2
a
b
Deq = D1 + ( ) 2 D2
a
1 1 b 1
= + ( )2
K eq K1 a K2
b 1
( ) s
u 2 (s) a meq
=
f (s) Deq 1
s2 + s+
meq meq k eq
Año sabático 2016-2017. AUTOR: DR. JUAN ELOY BENITO HERNÁNDEZ Página 30
Escriba las ecuaciones diferenciales del sistema de engranes y obtenga sus
transferencias.
d1 (t )
(t ) = J 1 + 1 (t )
dt
d 2 (t ) d 2 (t ) 1
KC
2 (t ) = J 2 + JC + DC 2 (t ) + 2 (t )dt
dt dt
Eliminando 1 (t ) , 2 (t ) y 2 (t )
1 (t ) n1
=
2 (t ) n2
1 (t ) n2
=
2 (t ) n1
Año sabático 2016-2017. AUTOR: DR. JUAN ELOY BENITO HERNÁNDEZ Página 31
n1 2 1
(s) = J 1 s1 (s) + ( ) [ J 2 s + J C s + DC + ]1 ( s)
n2 KC s
n1 2 n n 1
(s) = [[ J 1 + ( ) ( J 2 + J C )]s + ( 1 ) 2 DC + ( 1 ) 2 ]1 (s)
n2 n2 n2 K C s
1
( s) = [ J eq s + Deq + ]1 ( s)
K eq s
1
s
1 ( s ) J eq
=
( s) Deq 1
s2 + s+
J eq J eq K eq
n1 2
Donde J eq = J 1 + ( ) (J 2 + J C )
n2
n1 2
Deq = ( ) DC
n2
1 n 1
= ( 1 )2
K eq n2 K C
n1 1
(
) s
2 ( s) n2 J eq
=
(s) Deq 1
s2 + s+
J eq J eq K eq
Año sabático 2016-2017. AUTOR: DR. JUAN ELOY BENITO HERNÁNDEZ Página 32
Consideremos los siguientes sistemas mecánico traslacional, eléctrico de mallas y
eléctrico de nodos.
Año sabático 2016-2017. AUTOR: DR. JUAN ELOY BENITO HERNÁNDEZ Página 33
Si se analizan los tres sistemas anteriores se observa que sus ecuaciones son
enteramente similares y por lo tanto su solución es de la misma forma. Podríamos
presentar aquí una serie de sistemas cuyas ecuaciones resultan enteramente
similares a las anteriores.
Esta serie de sistemas de muy diversa naturaleza cuyas ecuaciones diferenciales
que los definen son de la misma forma, les llamamos sistemas análogos. La ventaja
que presentan estas analogías es la de pensar en términos de algún tipo de sistema
que se nos facilite y traducir conclusiones que de ellos saquemos a sistemas que
sean análogos a los primeros. Por ejemplo, si sabemos que la energía de la masa
mu 2
m debido a su velocidad u (t ) está dada por y determinamos que la masa es
2
análoga a la inductancia y la velocidad a la corriente en un circuito de mallas, la
Li 2
energía en una inductancia estará dada por . En un circuito de nodos la masa
2
es análoga a la capacitancia C y la velocidad al voltaje, por lo que la energía del
Cv 2
capacitor será .
2
En esta forma manejamos conceptos como ancho de banda, resonancia,
amortiguamiento, frecuencia natural, constante de tiempo, etcétera en forma
indistinta para sistemas como los mostrados anteriormente.
Los sistemas anteriores nos permiten establecer dos tipos de analogías la llamada
fuerza-voltaje y la llamada fuerza-corriente. La primera corresponde a la analogía
de un sistema mecánico con un circuito de mallas y la segunda de un sistema
mecánico con un circuito de nodos. El siguiente cuadro resume estas analogías
tanto para las variables como para los parámetros.
Año sabático 2016-2017. AUTOR: DR. JUAN ELOY BENITO HERNÁNDEZ Página 34
Para el sistema mecánico de la figura dibuje sus diagramas análogos fuerza-
voltaje y fuerza-corriente. Escriba sus ecuaciones en el dominio del tiempo y de
Laplace.
Analogía fuerza-voltaje.
En este tipo de analogía cada punto rígido del sistema mecánico, representado
por las masas m1 y m2 las cuales tienen una velocidad diferente u1 (t ) y u 2 (t ) , se
traduce en una malla con una corriente diferente.
Se observa que los elementos mecánicos comunes a los dos puntos rígidos deben
quedar situados en las ramas comunes a las mallas de su análogo eléctrico. En este
tipo de analogía por lo tanto, nunca pude quedar una masa entre dos mallas.
Analogía fuerza-corriente.
Año sabático 2016-2017. AUTOR: DR. JUAN ELOY BENITO HERNÁNDEZ Página 35
Ahora los puntos rígidos dan origen al número de nodos, con potenciales diferentes.
Los elementos comunes a dos puntos rígidos quedan entre dos nodos, por lo que
las masas no pueden quedar entre dos nodos.
Las ecuaciones del sistema mecánico pueden ser escritas fácilmente a partir de
cualquiera de sus análogos eléctricos.
du1 (t )
f (t ) = m1 + D1 (u1 (t ) − u 2 (t ))
dt
du 2 (t ) 1
K2
0 = D1 (u 2 (t ) − u1 (t )) + m2 + D2 u 2 (t ) + u 2 (t )dt
dt
En el dominio de Laplace.
Año sabático 2016-2017. AUTOR: DR. JUAN ELOY BENITO HERNÁNDEZ Página 36
D1u1 ( s)
u 2 ( s) =
1
m2 s + D1 + D2 +
K2s
2
D1
f ( s ) = (m1 s + D1 − )u1 ( s)
1
m2 s + D1 + D2 +
K2s
1 2 D1 + D2 1
s + s+
u1 ( s) m1 m1 m2 m1 m2 K 2
=
f ( s) D D D 1 DD D1
s 3 + ( 1 + 2 + 1 )s 2 + ( + 1 2 )s +
m1 m2 m2 m2 K 2 m1 m2 m1 m2 K 2
La otra transferencia es
D1
s
u 2 ( s) m1 m2
=
f ( s) D D D 1 DD D1
s 3 + ( 1 + 2 + 1 )s 2 + ( + 1 2 )s +
m1 m2 m2 m2 K 2 m1 m2 m1 m2 K 2
Año sabático 2016-2017. AUTOR: DR. JUAN ELOY BENITO HERNÁNDEZ Página 37
1.3.2.2. Analogías electromecánicas de elementos de acoplamiento.
En ambos sistemas se considera el caso ideal, esto es, sistemas que transmiten
íntegramente la energía, en un caso potencia mecánica y en otro potencia eléctrica.
Recordando las relaciones que guardan las fuerzas y las velocidades en relación a
las distancias al punto de apoyo.
f1 (t ) b
=
f 2 (t ) a
u1 (t ) a
=
u 2 (t ) b
Año sabático 2016-2017. AUTOR: DR. JUAN ELOY BENITO HERNÁNDEZ Página 38
Donde las relaciones son:
v1 (t ) N1
=
v2 (t ) N 2
i1 (t ) N 2
=
i2 (t ) N1
Año sabático 2016-2017. AUTOR: DR. JUAN ELOY BENITO HERNÁNDEZ Página 39
1.3.3. Analogias Rotacionales.
1 (t ) n1
=
2 (t ) n2
1 (t ) n2
=
2 (t ) n1
Año sabático 2016-2017. AUTOR: DR. JUAN ELOY BENITO HERNÁNDEZ Página 40
Como ejemplo de sistemas de acoplamiento dibujemos las analogías fuerza-
voltaje y fuerza-corriente del siguiente sistema rotacional.
El sistema mecánico tiene cuatro puntos rígidos y por lo tanto sus sistemas
análogos deben de tener cuatro mallas o cuatro nodos según el tipo de analogía.
En las siguientes figuras se representan los circuitos análogos, indicando los
parámetros mecánicos a que corresponden.
Año sabático 2016-2017. AUTOR: DR. JUAN ELOY BENITO HERNÁNDEZ Página 41
Las ecuaciones del sistema mecánico rotacional son:
(t ) 1
(t ) = J + D(t ) + [(t ) − 1 (t )]dt
dt K
1
0 = [1 (t ) − (t )]dt + 1 (t )
K
1
KC
2 (t ) = [ 2 (t ) − C (t )]dt
1 d C (t )
KC
0= [ C (t ) − 2 (t )]dt +J C + DC C (t )
dt
1 (t ) n1 1 (t ) n2
= =
2 (t ) n2 2 (t ) n1
Año sabático 2016-2017. AUTOR: DR. JUAN ELOY BENITO HERNÁNDEZ Página 42
1.3.4 Problemas propuestos
Año sabático 2016-2017. AUTOR: DR. JUAN ELOY BENITO HERNÁNDEZ Página 43
inercia despreciable. Dibuje sus analogías fuerza-voltaje y fuerza-corriente y
escriba las ecuaciones mecánicas del sistema.
D1
s
U 2 ( s) m1 m2
1. =
F ( s) D D D 1 DD D1
s 3 + [ 1 + 2 + 1 ]s 2 + [ + 1 2 ]s +
m2 m2 m1 K 2 m2 m1 m2 K 2 m1 m2
D a+b
( )s
U 2 ( s) mC a
2. =
U ( s) a + b 2 D DC 1
s 2 + (( ) + )s +
a mC mC K C mC
Año sabático 2016-2017. AUTOR: DR. JUAN ELOY BENITO HERNÁNDEZ Página 44
1 1
3. ( s ) = ( D1 + D2 + )( s ) − ( D1 + ) 1 ( s ) − D2 2 ( s )
K1 s K1 s
1 1 1 1
0 = −( D1 + )( s) + ( J 1 s + D1 + + ) 1 ( s ) − 2 (s)
K1 s K1 s K 2 s K2s
1 1
0 = − D2 ( s ) − 1 ( s ) + ( J 2 s + + D2 ) 2 ( s )
K2s K2s
n1 1
( )
2 ( s) n2 K 0 J eq
4. =
( s ) Deq 1
s2 + s+
J eq K eq J eq
Año sabático 2016-2017. AUTOR: DR. JUAN ELOY BENITO HERNÁNDEZ Página 45
Año sabático 2016-2017. AUTOR: DR. JUAN ELOY BENITO HERNÁNDEZ Página 46