π
<-
Chat plein-écran
[^]

cccccccc2


Hierarchy of files

 Downloads
 Files created online(24914)
 TI-Nspire
(18491)

 mViewer GX Creator Lua(13257)

DownloadTélécharger


LicenceLicense : Non spécifiée / IncluseUnspecified / Included

 TéléchargerDownload

Actions



Vote :

ScreenshotAperçu


Tester en ligne !

Informations

Catégorie :Category: mViewer GX Creator Lua TI-Nspire
Auteur Author: CristianLazoQ
Type : Classeur 3.6
Page(s) : 33
Taille Size: 1.77 Mo MB
Mis en ligne Uploaded: 12/07/2019 - 19:20:23
Uploadeur Uploader: CristianLazoQ (Profil)
Téléchargements Downloads: 1
Visibilité Visibility: Archive publique
Shortlink : http://ti-pla.net/a2270076

Description 

MODELADO, SIMULACION Y CONTROL
DE SISTEMAS DINAMICOS

P.F.PULESTON y F.VALENCIAGA




Nota: Este apunte tiene por objetivo principal introducir al modelado,
simulación y control de sistemas dinámicos empleando Matlab. Gran
parte del material aquí presentado esta basado en los Control Tutorial
for Matlab de la Universidad de Michigan, Demos y Manuales de
Matlab. Abundante información adicional y soluciones a muchos otros
problemas vinculados al control de sistemas, procesamiento de señal y
calculo matricial en general, puede encontrarse en dicha bibliografía.



Parte II. Ejemplo de Aplicación Al Control De
Sistemas

Sección 1: Modelado de un Motor de CC
Ecuaciones físicas del sistema
Requerimientos de Diseño
Representación en MATLAB y respuesta de lazo cerrado


Ecuaciones físicas del sistema
Un actuador mecánico muy difundido es el motor de CC. Provee directamente
movimiento rotacional y, adecuadamente acondicionado, movimiento traslacional.
El circuito eléctrico de armadura y el diagrama mecánico rotacional, se muestran
en la figura:
Para el ejemplo se consideraron los siguientes parámetros:

* momento de inercia del sistema (J) = 0.01 kg.m^2/s^2
* coeficiente de roce (b) = 0.1 Nms
* constante de fuerza electromotriz (K=Ke=Kt) = 0.01 Nm/Amp
* resistencia de armadura (R) = 1 ohm
* inductancia de armadura (L) = 0.5 H
* entrada (V): Fuente de Tensión
* posición del eje: Θ
* Se supone rotor y eje rígidos.
La cupla (T) está relacionada con la corriente de armadura y la fem (e) con la
velocidad de rotación, según las ecuaciones:

T = Kt ⋅ i


e = Ke ⋅ θ

siendo ambas constantes iguales (Kt=Ke=K)

En base a la ley de Newton y la ley de Kirchoff, resultan las siguientes ecuaciones
diferenciales que describen la dinámica del sistema:

•• •
J ⋅ θ+ b ⋅ θ = K ⋅ i

di •
L⋅ + R ⋅i = V − K ⋅θ
dt
1. Función de Transferencia
Aplicando la Transformada de Laplace y haciendo cero las condiciones iniciales,
las ecuaciones del sistema quedan expresadas en el dominio de s:

s ⋅ ( J ⋅ s + b ) ⋅ Θ(s ) = K ⋅ I (s )

( L ⋅ s + R ) ⋅ I ( s ) = V − K ⋅ s ⋅ Θ( s )
Eliminando I(s) se obtiene la transferencia entre la entrada de tensión de armadura

V y la velocidad de rotación Θ como salida:


Θ K
=
V ( J ⋅ s + b ) ⋅ (L ⋅ s + R ) + K 2

2. Espacio de Estados
La descripción del sistema de estados en el dominio temporal puede obtenerse

definiendo las variables físicas velocidad de rotación θ(t ) y corriente de armadura
i(t), como variables de estado, la tensión de armadura v(t) como entrada y la
velocidad de rotación como salida:

 b K  •
− J
 •  
d θ
 = J  ⋅ θ  +  0  ⋅ v
dt  i  − K R  
  −   i  1 L 
 L L

Requerimientos de Diseño
El motor sin compensar puede rotar solamente a 0,1rad/s con una entrada de 1 Volt
(ver simulación de la planta a lazo abierto). Uno de los requerimientos es que en
estado estacionario presente un error respecto de la velocidad deseada menor que
el 1%. Dinámicamente se espera un tiempo de establecimiento de 2 seg y un
sobrepaso menor que el 5% para evitar daños en la máquina. Es decir:

• tiempo de establecimiento de 2 seg
• sobrepaso menor que el 5%

• Error de estado estacionario 1%

Representación en Matlab y respuesta de lazo cerrado
1. Función de Transferencia
Para representar la función de transferencia es necesario considerar los polinomios
numerador y denominador:

num = K
den = ( J ⋅ s + b ) ⋅ ( L ⋅ s + R ) + K 2
Creando un archivo .m:

J=0.01;
b=0.1;
K=0.01;
R=1;
L=0.5;
num=K;
den=[(J*L) ((J*R)+(L*b)) ((b*R)+K^2)];

Para observar la respuesta al escalón del sistema a lazo abierto, basta con agregar
al archivo los siguientes comandos y ejecutarlo:

Step(num,den,0:0.1:3)
title('Step Response for the Open Loop System')

Resultando:
De la figura se observa que a lazo abierto se obtiene una salida 10 veces más chica
que la deseada (0,1 rad/s) y 3 seg de establecimiento, sin cumplir las
especificaciones.

2. Espacio de Estados
Análogamente se puede crear el siguiente archivo .m y ejecutarlo, resultando la
misma figura.

J=0.01;
b=0.1;
K=0.01;
R=1;
L=0.5;
A=[-b/J K/J
-K/L -R/L];
B=[0
1/L];
C=[1 0];
D=0;

step(A, B, C, D)

Nota: dependiendo la versión de Matlab la función step puede requerir de la
creación de un sistema SYS:
SYS=SS(A,B,C,D);

step(SYS)
Sección 2: Método de Diseño de un PID para
Control de Velocidad de un Motor CC
Control Proporcional
Control PID
Sintonía


Recordando las ecuaciones dinámicas y transferencia del motor de CC:

s ⋅ ( J ⋅ s + b ) ⋅ Θ(s ) = K ⋅ I (s )

( L ⋅ s + R ) ⋅ I ( s ) = V − K ⋅ s ⋅ Θ( s )

Θ K
=
V ( J ⋅ s + b ) ⋅ (L ⋅ s + R ) + K 2

resulta el siguiente diagrama esquemático del sistema compensado:




El sistema a lazo cerrado debe cumplir las especificaciones establecidas
anteriormente:

• tiempo de establecimiento de 2 seg

• sobrepaso menor que el 5%

• Error de estado estacionario 1%

El objetivo de esta sección será diseñar un controlador PID que permita verificar
las especificaciones. Como primer paso crearemos el siguiente archivo .m:
J=0.01;
b=0.1;
K=0.01;
R=1;
L=0.5;
num=K;
den=[(J*L) ((J*R)+(L*b)) ((b*R)+K^2)];

Correspondiente al sistema a lazo abierto y recordando que la transferencia del PID
es de la forma:

KI K D ⋅ s2 + K P ⋅ s + K I
KP + + KD ⋅ s =
s s



Control Proporcional
Probemos primero un control proporcional con ganancia 100. Para ello hay que
agregar los siguientes comandos al final del archivo .m:

Kp=100;
numa=Kp*num;
dena=den;

Para determinar la transferencia de lazo cerrado se puede utilizar el comando
cloop:

[numac,denac]=cloop(numa,dena);

siendo los vectores numac y denac el numerador y denominador de la transferencia
a lazo cerrado

Nota: dependiendo la versión de Matlab la función cloop puede no existir y ser
necesario recurrir al comando feedback:

SYS1=tf(numa,dema)
SYS2=tf(1,1)
SYSlc=feedback(SYS1,SYS2);
Para ver la respuesta de lazo cerrado al escalón de entrada basta con agregar al
archivo .m los siguientes comandos:

t=0:0.01:5;
step(numac,denac,t)
title('Step response with Proportion Control')

Resultando:




Control PID
Puede observarse que el error de estado estacionario y el sobrepaso son excesivos.
Es sabido que agregando un término integral se elimina el error de estado
estacionario al escalón, mientras que un término derivativo, adecuadamente
sintonizado, puede reducir el sobrepaso. Probemos entonces con un controlador
PID con bajas ganancias KI y KD. Para ello hay que modificar el archivo .m de la
siguiente manera:

J=0.01;
b=0.1;
K=0.01;
R=1;
L=0.5;
num=K;
den=[(J*L) ((J*R)+(L*b)) ((b*R)+K^2)];

Kp=100;
Ki=1;
Kd=1;

numc=[Kd, Kp, Ki];
denc=[1 0];
numa=conv(num,numc);
dena=conv(den,denc);
[numac,denac]=cloop(numa,dena);
step(numac,denac)
title('PID Control with small Ki and Kd')

resultando




Sintonía
Ahora el tiempo de establecimiento es muy largo, por lo que incrementaremos KI a
200. Modificando este valor en nuestro archivo y ejecutándolo se obtiene:
siendo una respuesta mucho más rápida que la anterior, pero a costa de aumentar el
sobrepaso. En consecuencia, para reducirlo, incrementaremos el valor de KD a 10
en el archivo .m:




Resultando:

Kp=100,
Ki=200,
Kd=10,
una sintonía adecuada para cumplir los requerimientos de diseño.
Sección 3: Diseño del Compensador Empleando
el Método de Lugar de Raíces
Graficación del Diagrama de Lugar de Raíces
Obtención de la ganancia usando el comando rlocfind
Compensador por atraso de fase
Graficación de la respuesta de lazo cerrado


En esta sección trataremos el diseño de un compensador para el motor de CC,
basándonos en el método del lugar de raíces. Consideremos las mismas
especificaciones de diseño. Como en la sección anterior el primer paso es crear un
archivo .m que contenga el modelo de la planta a lazo abierto.




J=0.01;
b=0.1;
K=0.01;
R=1;
L=0.5;
num=K;den=[(J*L) ((J*R)+(L*b)) ((b*R)+K^2)];

Graficación del Diagrama de Lugar de Raíces
La idea principal del método es compensar agregando polos y/o ceros
adecuadamente de modo tal que la respuesta a lazo ...

Archive contentsContenu de l'archive

Action(s) SizeTaille FileFichier
1.64 Ko KB readme.txt
607.03 Ko KB cccccccc2/21-30.tns
600.70 Ko KB cccccccc2/11-20.tns
526.93 Ko KB cccccccc2/01-10.tns
160.83 Ko KB cccccccc2/31-33.tns

Pub / Ads

-
Search
-
Featured topics
Offre TI-Planet/Jarrety pour avoir la TI-83 Premium CE avec son chargeur pour 79,79€ port inclus !
Offre TI-Planet/Jarrety pour avoir la TI-Nspire CX CAS à seulement 130€ TTC port inclus!
Jailbreake ta TI-Nspire avec Ndless et profite des meilleurs jeux et applications !
123
-
Donations / Premium
For more contests, prizes, reviews, helping us pay the server and domains...

Discover the the advantages of a donor account !
JoinRejoignez the donors and/or premium!les donateurs et/ou premium !


Partner and ad
Notre partenaire Jarrety 
-
Stats.
329 utilisateurs:
>311 invités
>13 membres
>5 robots
Record simultané (sur 6 mois):
6892 utilisateurs (le 07/06/2017)
-
Other interesting websites
Texas Instruments Education
Global | France
 (English / Français)
Banque de programmes TI
ticalc.org
 (English)
La communauté TI-82
tout82.free.fr
 (Français)