C:\Documents and Settings\Sever Serban\My Documents\carti ...
C:\Documents and Settings\Sever Serban\My Documents\carti ...
C:\Documents and Settings\Sever Serban\My Documents\carti ...
Create successful ePaper yourself
Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.
C<br />
A<br />
P<br />
I<br />
T<br />
O<br />
L<br />
U<br />
L<br />
1<br />
modelarea ªi simularea<br />
proceselor ªi sistemelor<br />
1.1. noÞiuni fundamentale<br />
Disciplina Tehnici de Modelare şi Simulare a Proceselor şi Sistemelor (TMSPS) se<br />
ocupă cu studiul metodelor şi tehnicilor prin care, procese din lumea reală, se<br />
reprezintă sub formă de modele matematice abstracte, pe baza cărora, prin simulare,<br />
se pot obţine informaţii pertinente asupra procesului analizat; utilitatea modelelor rezidă<br />
atât în analiza comportamentului dinamic din trecutul evoluţiei sale temporale, cât, mai<br />
ales, în scopul acţiunilor de prognoză a comportamentului său. Utilitatea cea mai<br />
elevată a modelelor de procese este reprezentată de accepţiunea de sistem -un model<br />
matematic cu o formă specifică de reprezentare- ce permite sinteza structurilor de<br />
conducere automată a proceselor implicate.<br />
Prin noţiunea de PROCES se înţelege un ansamblu de elemente concrete -uzual, de<br />
factură artificială- inclusiv fenomenele fizico- chimice implicate, cu o destinaţie<br />
economică concretă cum ar fi: producere de energie, conversie între diverse forme de<br />
energie, producerea de bunuri materiale, transport, etc. Este important faptul că, un<br />
proces, poate fi privit ca o structură de sine stătătoare, separabilă de mediul<br />
înconjurător, cu care dialoghează doar prin fluxurile de intrare în proces, respectiv de<br />
ieşire din acesta. Multitudinea proceselor -cu care ne putem confrunta într-un astfel de<br />
studiu- este foarte diversă, cele mai importante fiind:<br />
Procesele industriale sunt destinate producţiei de forme finite de energie sau<br />
bunuri materiale, reprezentând pseudototalitatea obiectelor ce fac obiectul<br />
modelării matematice, cu scopul declarat al realizării conducerii automate;<br />
PDF created with pdfFactory Pro trial version www.pdffactory.com
2 S. ERBAN MODELARE i SIMULARE<br />
Procesele economice sunt specifice dezvoltării societăţilor productive,<br />
economiilor naţionale şi chiar structurilor mondiale; caracterizate -uzualvaloric<br />
-monetar- sunt procese de o complexitate deosebită, caracterizarea lor<br />
presupunând cunoştinţe de specialitate dar şi sistemice deosebite. Este<br />
domeniul în care prognoza dezvoltării, obţinută pe baza modelării matematice,<br />
este cea mai intuitivă şi elocventă;<br />
Procesele biologice sunt specifice ştiinţelor vieţii care, privite din perspectiva<br />
modelării matematice, oferă perspective nebănuite pentru diagnoza<br />
deficienţelor de funcţionare precum şi a metodelor de corijare.<br />
Dacă este necesară o exemplificare a noţiunii de proces, cea mai sugestivă ilustrare este<br />
dată de cea a proceselor industriale. Astfel<br />
Exemplul 1.1. Grupul turbină -generator reprezentat în figura 1.1 este evident un<br />
proces industrial destinat producerii de energie electrică prin consumul unei alte<br />
forme de energie cum ar fi: energia termică a aburului -grup turbogenerator, energia<br />
hidraulica a apei -grup hidrogenerator, etc. Procesul poate fi privit ca izolat de mediul<br />
înconjurător, cu care dialoghează doar prin fluxul de energie de intrare (mărime cauză)<br />
în acest proces care este (se face referire la grupul turbogenerator) energia termică,<br />
caracterizată de debitul de abur, temperatura şi presiunea acestuia, respectiv fluxul de<br />
energie de ieşire (mărimea efect), caracterizat de tensiune, curent şi frecvenţă, mărimi<br />
ce determină puterea electrică<br />
produsă de generator.<br />
Exemplul 1.2. Procesul industrial<br />
din figura 1.2 este reprezentat de<br />
un amestecător în care, pe baza a<br />
două debite de intrare Q1, Q2, cu<br />
concentraţii în aceeaşi substanţă<br />
activă c1, c2 (şi care reprezintă<br />
fluxul de intrare w Fig. 1.1. Grup turbin-generator<br />
i ) se obţine un<br />
debit de ieşire Q de concentraţie c<br />
ca proces industrial<br />
(fluxul de ieşire we ), singurele cu<br />
care se produce interacţiunea sa cu mediul înconjurător.<br />
PDF created with pdfFactory Pro trial version www.pdffactory.com
Noiuni fundamentale 3<br />
Fig. 1.2. Amestectorul ca proces industrial<br />
Exemplul 1.3. Un algoritm de calcul implementat pe un calculator numeric, ce<br />
prelucrează “sincronizat” (toate semnalele sunt definite la aceleaşi momente de timp,<br />
precizate de ceasul de timp real) semnalele numerice este absolut similar unui proces,<br />
cu specificităţi datorate suportului discret al funcţiilor.<br />
Exemplul 1.4. Un proces economic este reprezentat, de exemplu, de balanţa<br />
legăturilor între ramurile unei economii, producţiile acestora fiind p i : 1=1, 2, ..., n.<br />
Un MODEL al unui proces este o reprezentare, într-o formă utilizabilă, a aspectelor<br />
esenţiale legate de funcţionalitatea acestuia, folosind mărimile specifice ale procesului.<br />
Este de la sine înţeles că orice model dă o reprezentare simplificată a procesului<br />
implicat, că întotdeauna obţinerea acestuia este orientată spre un grad prestabilit de<br />
complexitate. De exemplu, deşi funcţionalităţile generale sunt, practic, întotdeauna<br />
neliniare, se caută modele liniare care să exprime suficient de exact funcţionalitatea<br />
procesului. Modul de caracterizare a diverselor procese, generează diverse tipuri de<br />
modele:<br />
modele empirice, folosite în general pentru fizica şi chimia microscopică;<br />
modele fenomenologice la care fundamentale sunt legile fizice şi chimice<br />
specifice domeniului de apartenenţă a procesului vizat;<br />
modele funcţionale pentru care reprezentările se fac prin relaţii sau scheme<br />
funcţionale;<br />
modele matematice sunt acelea pentru care reprezentările funcţionalităţilor<br />
sunt ecuaţionistice (analitice).<br />
Prin noţiunea de SISTEM se înţelege o formă rafinată de model al unui proces, în<br />
sensul unificării denumirilor mărimilor, a notaţiilor precum şi a modalităţilor de<br />
PDF created with pdfFactory Pro trial version www.pdffactory.com
4 S. ERBAN MODELARE i SIMULARE<br />
reprezentare. Pentru sisteme se pot construi teorii generale de analiză sau sinteză, cu<br />
formulări ce nu mai sunt tarate de natura procesului concret (şi deci, a modelului<br />
ataşat). Sistemele se pot descompune în SUBSISTEME sub restricţia ca legăturile de<br />
interconectare să fie numai de factură informaţională (fără transfer de masă sau<br />
energie); în egală măsură, din sisteme, se pot realiza structuri complexe, prin conexiuni<br />
informaţionale, menite să confere ansamblului proprietăţi superioare faţă de sistemul<br />
originar.<br />
1.2. modelare ªi simulare<br />
Prin MODELARE înţelegem operaţia de determinare a unui model ataşat unui proces<br />
concret, într-o formă adecvată de exprimare şi precizie, ce permite exploatarea sa<br />
pentru scopurile propuse. Obişnuit, modelul se rescrie imediat sub formă de sistem, pur<br />
şi simplu prin renotarea mărimilor implicate precum şi prin adaptarea exprimării<br />
acestuia la un şablon specific.<br />
Cât priveşte SIMULAREA, aceasta este o metodă experimental- aplicativă prin care<br />
se testează un model (sistem) pentru a stabili gradul de adecvare a acestuia în raport<br />
cu procesul iniţial. Simularea se realizează prin implementarea sistemului (modelului)<br />
pe un dispozitiv de simulare care poate să fie specializat (ataşat unui sistem anume) sau<br />
universal (calculator analogic sau numeric). Este cazul să menţionăm dezvoltarea fără<br />
precedent a tehnicilor numerice de simulare, care pun la dispoziţia cercetătorilor<br />
programe dedicate acestei operaţii, cu o interfaţare prietenoasă accesibilă chiar şi<br />
începătorilor în domeniul informatic.<br />
Operaţiile de modelare şi simulare se subordonează schemei din fiura 1.3 care<br />
sugerează inclusiv necesitatea perfecţionării modelelor pentru a satisface o serie de<br />
criterii de calitate impuse.<br />
Fig.1.3. Structura cibernetică asociată activităţilor de modelare şi simulare<br />
Cât priveşte modalităţile de simulare pentru un proces fizic existent, acestea sunt<br />
evidenţiate în figura 1.4.<br />
PDF created with pdfFactory Pro trial version www.pdffactory.com
Noiuni fundamentale 5<br />
Fig.1.4. Schemă funcţională ce evidenţiază relaţia proces -simulare<br />
Se menţionează că, trecerea de la proces la sistem, se poate face pe două căi principale<br />
(vezi figura 1.4):<br />
# Prin experimente pe procesul (sau modelul) fizic existent, urmate de prelucrări<br />
adecvate ale datelor pentru a obţine exprimări analitice ale sistemelor asociate;<br />
# Prin modelare şi simulare, într-un proces iterativ, până la atingerea unui grad<br />
de acurateţe impus prin criterii adecvate.<br />
De asemenea, se atrage atenţia încă o dată că:<br />
• Sistemele sunt exprimări coerente şi aproximative a funcţionalităţilor<br />
proceselor;<br />
• Sunt preferate exprimările prin sisteme liniare, chiar dacă procesele modelate<br />
conţin neliniarităţi, datorită avantajelor pe care le conferă acestea în analiza şi<br />
sinteza structurilor automate de conducere;<br />
• Acceptarea, de facto, a exprimărilor sistemice aproximative va pretinde, în<br />
plus faţă de alte deziderate, ca structura de conducere să atenueze suficient<br />
incertitudinele modelului sistemic -parametrice sau de structură. De aceea<br />
structurile de conducere automată au ca element definitoriu reacţia (feedback),<br />
producerea erorii instantanee faţă de programul de conducere şi elaborarea<br />
comenzii printr-o prelucrare adecvată a acesteia.<br />
1.3. clasificarea sistemelor (modelelor)<br />
PDF created with pdfFactory Pro trial version www.pdffactory.com
6 S. ERBAN MODELARE i SIMULARE<br />
În tot ce urmează în continuare se practică următoarele notaţii specifice teoriei<br />
sistemelor:<br />
Mărimea de stare (stare): x n ;<br />
Mărimea de intrare (de com<strong>and</strong>ă): u m ;<br />
Mărimea de perturbaţie: v r ;<br />
Mărimea de ieşire (măsurată): y p ;<br />
Variabila independentă timp: t (mai exact t + ) sau, după caz, t (în<br />
fapt real t ).<br />
Există o serie întreagă de criterii de clasificare a sistemelor (modelelor). Într-o listă, ce<br />
nu are pretenţia că este exhaustivă, acestea sunt:<br />
1. După natura modelului (sistemului): fizice, fenomenologice (conceptuale),<br />
simbolice (formale) sau matematice (analitice).<br />
2. După dependenţa de timp a mărimilor:<br />
- sisteme statice (mărimile sunt constante -nu depind de timp)<br />
Yf X , U (1.1)<br />
- sisteme dinamice (toate mărimile depind de timp -vezi în continuare).<br />
3. După varianţa parametrilor sistemului:<br />
- sisteme variante (parametrii depind de timp)<br />
- sisteme invariante (parametrii săi sunt constanţi).<br />
4. După liniaritatea dependenţelor:<br />
- sisteme neliniare<br />
- sisteme liniare.<br />
5. După caracterul variabilei independente timp:<br />
- sisteme netede sau continue (t sau, mai exact, t + )<br />
- sisteme (cu timp) discrete (t sau, mai precis, t ).<br />
Cu aceste criterii se evidenţiază următoarele clase de sisteme:<br />
G Sisteme neliniare, netede şi variante<br />
x t f t,x t,ut,vt<br />
yt gt,x t<br />
cu x t <br />
G Sisteme neliniare, discrete şi variante<br />
d x t<br />
dt<br />
PDF created with pdfFactory Pro trial version www.pdffactory.com<br />
(1.2)
Noiuni fundamentale 7<br />
x t1 f t,x t,ut,vt<br />
yt gt,x t<br />
G Sisteme liniare, netede şi variante<br />
x t A t x t B t u t E t v t<br />
y t C t x t<br />
(1.3)<br />
cu A t nxn , B t nxm C t pxn , E t nxr (1.4)<br />
G Sisteme liniare, discrete şi variante (relaţia de mai jos este (1.5))<br />
x t1 A t x t B t u t E t v t<br />
y t C t x t<br />
G Sisteme neliniare, netede şi invariante<br />
x t f x t,ut,vt<br />
yt gxt<br />
cu x t <br />
d x t<br />
dt<br />
G Sisteme neliniare, discrete şi invariante<br />
x t1 f x t,ut,vt<br />
yt gxt<br />
G Sisteme netede, liniare şi invariante<br />
x t A xt Bu t Ev t<br />
y t Cxt<br />
cu<br />
G Sisteme discrete, liniare şi invariante<br />
x t1 A xt Bu t Ev t<br />
y t Cxt<br />
cu<br />
cu A t nxn , B t nxm<br />
C t pxn , E t nxr<br />
A nxn , B nxm<br />
C pxn , E nxr<br />
A nxn , B nxm<br />
C pxn , E nxr<br />
(1.6)<br />
(1.7)<br />
(1.8)<br />
(1.9)<br />
Se poate adăuga şi următoarele criterii:<br />
6. După dimensiunea mărimii de intrare şi/sau a celei de ieşire:<br />
- sisteme cu mai multe intrări şi mai multe ieşiri (multi I/E sau, încă, după<br />
denumirea din engleză MIMO -multiple inputs- multiple outputs) sunt toate<br />
sistemele menţionate anterior;<br />
- sisteme cu o intrare şi o ieşire (mono I/E sau, după literatura de limbă<br />
engleză, SISO -single input- single output) ce au ca particularitate m=1, u,<br />
r=1, v şi p=1 adică y ; se practică notaţiile particulare<br />
PDF created with pdfFactory Pro trial version www.pdffactory.com
8 S. ERBAN MODELARE i SIMULARE<br />
G Sisteme neliniare, netede sau discrete şi variante<br />
x / t f t,x t,u t,v t<br />
y t g t,x t<br />
cu x / t <br />
d x t<br />
dt<br />
, neted<br />
x t1 , discret<br />
G Sisteme neliniare, netede sau discrete şi invariante<br />
x / t f x t,u t,v t<br />
y t g x t<br />
cu x / t <br />
d x t<br />
dt<br />
, neted<br />
x t1 , discret<br />
G Sisteme liniare, netede sau discrete şi variante<br />
x / t A t x t b tut etvt<br />
yt c T t xt<br />
G Sisteme liniare, netede sau discrete şi invariante<br />
x / t A xt but evt<br />
yt c T xt<br />
cu<br />
cu<br />
(1.10)<br />
(1.11)<br />
A t nxn , b t n cT t 1xn , e t n (1.12)<br />
A nxn , b n<br />
c T 1xn , e n<br />
(1.13)<br />
7. După dependenţa variabilelor de o coordonată geometrică generalizată:<br />
- sisteme cu parametri concentraţi (finit dimensionale), sunt toate modelele<br />
prezentate anterior;<br />
- sistemele cu parametri distribuiţi se asociază proceselor la care mărimile<br />
caracteristice, la un moment de timp precizat, depind şi de o coordonată<br />
geometrică.<br />
Exemplul 1.1. Fie procesul fizic de încălzire a unei bare matalice pe la capetele sale,<br />
folosind temperaturile T 1 şi T 2 , aşa cum se prezintă în figura 1.5.<br />
PDF created with pdfFactory Pro trial version www.pdffactory.com
Noiuni fundamentale 9<br />
Fig.1.5. Distribuia temperaturii într-o bar<br />
înclzit pe la capete<br />
Intuitiv chiar, se vede că temperatura din bară este funcţie de timp dar şi de coordonata<br />
geometrică, aici cu o singură dimensiune . Modelarea procesului de încălzire<br />
presupune aplicarea ecuaţiei difuziei căldurii la cazul monodimensional, rezultând<br />
t,<br />
t<br />
a 2 t,<br />
2<br />
;t , L,L ;a <br />
c >0 (1.14)<br />
care se completează cu condiţia iniţială (specificată )<br />
0, f , L,L (1.15)<br />
şi condiţia la limită (la frontieră) scrisă, în acest caz, pe baza ipotezei Fourier<br />
t,<br />
<br />
<br />
L<br />
t,<br />
<br />
<br />
L<br />
1 T 1 t t,L<br />
2 T 2 t t, L<br />
(1.16)<br />
În aceste ecuaţii este conductibilitatea termică, este coeficient de transmitere a<br />
căldurii printr-o suprafaţă, este greutatea specifică, iar c este căldura specifică.<br />
Desigur că starea este x(t, ) = (t, ) iar comenzile sunt u 1 = T 1 , u 2 = T 2 dar, modelul<br />
matematic (şi deci sistemul dinamic asociat) este de o structură cu totul aparte. Având<br />
în vedere faptul că [ -L, L ] (o mulţime nenumărabilă) astfel de sisteme se<br />
numesc şi infinit dimensionale. <br />
8. După modul de caracterizare a modelului ataşat procesului analizat:<br />
- sistemele structurale sunt cele descrise pe stare, de forma celor exprimate<br />
anterior, deoarece starea x(t) este acceptată ca fiind asociată direct mărimilor<br />
interne procesului;<br />
PDF created with pdfFactory Pro trial version www.pdffactory.com
10 S. ERBAN MODELARE i SIMULARE<br />
- sisteme funcţionale sunt cele la care exprimarea se face intrare -ieşire (I/E<br />
sau I/O -în engleză), numai pe baza mărimilor de intrare şi de ieşire, fără a<br />
implica cu nimic structura internă a procesului; în acest sens se pot evidenţia<br />
G Sisteme liniare SISO, netede, caracterizate în domeniul timp<br />
n<br />
<br />
i0<br />
a i<br />
m<br />
d i y t<br />
dti <br />
j0<br />
b j<br />
d j u t<br />
dt j<br />
, n>m (1.17)<br />
G Sisteme liniare SISO, discrete, caracterizate în domeniul timp<br />
n<br />
<br />
i0<br />
m<br />
ai y t i <br />
j0<br />
b j u t j , n>m (1.18)<br />
G Sisteme liniare SISO, netede, caracterizate cu transformata Laplace (prin<br />
funcţie de transfer)<br />
H s b m s m b m1 s m1 b 1 sb 0<br />
a n s n a m1 s n1 a 1 sa 0<br />
<br />
m<br />
<br />
j0<br />
n<br />
<br />
i0<br />
b j s j<br />
a i s i<br />
, n>m (1.19)<br />
G Sisteme liniare SISO, discrete, caracterizate cu transformată Z (prin funcţie de<br />
transfer)<br />
H z b m z m b m1 z m1 b 1 zb 0<br />
a n z n a m1 z n1 a 1 za 0<br />
<br />
m<br />
<br />
j0<br />
n<br />
<br />
i0<br />
b j z j<br />
a i z i<br />
, n>m (1.20)<br />
ultimele reprezentând, în fapt, aplicarea transformatei Laplace respectiv Z modelelor<br />
din timp menţionate anterior.<br />
9. După modul de explicitare a modelului se obţin:<br />
- modele parametrice sunt toate cele prezentate până la acest moment,<br />
comunicarea modelului făcându-se prin ecuaţii formulate pe baza parametrilor<br />
acestuia;<br />
- modele neparametrice se precizează prin elemente globale, ce nu ţin de<br />
explitarea bazată pe valorile parametrilor modelului; astfel de modele sunt<br />
răspunsul cauzal la impulsul unitar are strucura din figura 1.6<br />
PDF created with pdfFactory Pro trial version www.pdffactory.com
Noiuni fundamentale 11<br />
Fig.1.6. Structura unui răspuns la impuls<br />
răspunsul indicial este ieşirea sistemului când intrarea este de formă treaptă<br />
unitară, forma tipică fiind cea din figura 1.7.<br />
Fig.1.7. Trei forme posibile pentru<br />
răspunsul indicial<br />
răspunsul în frecvenţă al modelului reprezentat, uzual, sub formă de hodograf<br />
(figura 1.8) sau caracteristici logaritmice (figura 1.9).<br />
PDF created with pdfFactory Pro trial version www.pdffactory.com
12 S. ERBAN MODELARE i SIMULARE<br />
1.4. limbaje numerice de simulare<br />
Fig.1.8. Hodograf, ca model neparametric<br />
Fig.1.9. Caracteristici logaritmice,<br />
ca model neparametri<br />
Simularea numerică a funcţionalităţilor sistemelor a cunoscut până în prezent o<br />
dezvoltare fără precedent. O serie întreagă de limbaje specializate sunt destinate<br />
PDF created with pdfFactory Pro trial version www.pdffactory.com
Noiuni fundamentale 13<br />
simulării celor mai diverse structuri de sisteme liniare sau neliniare, dintre care<br />
se menţionează:<br />
MATLAB -un program compact ce dispune de o bibliotecă matematică<br />
performantă, elaborată pe principii de precizie şi acurateţe a calculelor<br />
numerice matriciale, are ataşate toolbox-uri specializate pe domenii<br />
specifice cum ar fi: electrotehnică, chimie, analiza semnalelor, analiză<br />
spectrală, teoria sistemelor (răspunsul în timp al sistemelor liniare şi<br />
neliniare de o manieră generală sau particulară (răspuns indicial, răspuns<br />
la impuls, răspuns liber, etc.), răspuns în frecvenţă, sinteza conducerii<br />
automate clasice, sinteza compensatoarelor dinamice pentru sistemele<br />
reprezentate pe stare (legi de reacţie, estimatoare de stare), conducerea<br />
robustă a sistemelor, etc.; grafica oferită de acest program este de<br />
asemenea sugestivă, reprezentativă şi usor de apelat;<br />
SIMULINK-ul din MATLAB este un toolbox specializat pentru<br />
modelarea sistemelor dinamice, de factură liniară sau neliniară, cu<br />
inovaţia esenţială a reprezentării sistemelor pe blocuri, similare<br />
funcţionalităţilor analogice (sumatoare, integratoare, conexiuni, aparate<br />
de vizualizare a semnalelor -aparate de măsură, osciloscoape cu mai<br />
multe canale, oscilografe, etc.), pentru o adaptare uşoară la tehnica<br />
numerică a mulţilor specialişti care au utilizat -timp îndelungat- tehnica<br />
analogică de modelare;<br />
LABVIEW este un program numeric care permite simularea structurilor<br />
de conducere automată, respectând structura fizică a buclelor de<br />
automatizare, relevând toată problematica aferentă acestei activităţi:<br />
sisteme de achiziţie a datelor (traductoare şi senzori), filtrarea<br />
semnalelor achiziţionate din procesul condus, modelarea sistemelor de<br />
execuţie, a compensatoarelor dinamice precum şi vizualizarea<br />
răspunsurilor dinamice ale acestora;<br />
MATHCAD excelează prin grafica tridimensională fiind recom<strong>and</strong>at<br />
proiectărilor tehnologice;<br />
MATEMATICA la care calculul matematic este cel vizat cu preponderenţă,<br />
apărând elemente de certă particularitate cum ar fi cel de calcul formal al<br />
primitivelor unor funcţii st<strong>and</strong>ard;<br />
SIMNON;<br />
MAPLE;<br />
SIMAN<br />
ultimele fiind programe cu o răspândire zonală redusă.<br />
PDF created with pdfFactory Pro trial version www.pdffactory.com
14 S. ERBAN MODELARE i SIMULARE<br />
În cadrul acestei discipline, pentru modelarea sistemelor, se va utiliza numai programul<br />
MATLAB, inclusiv modulul denumit SIMULINK.<br />
PDF created with pdfFactory Pro trial version www.pdffactory.com