17.04.2013 Views

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 ...

SHOW MORE
SHOW LESS

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

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!