Әдеби деректерге және нысанға орындалған алғашқы зерттеулерге сүйеніп және теңдеулер құрылымы температура үшін - сызықтық және қысымтүсуі үшін бейсызықты квадратты құрылым ретінде қабылдаймыз. Жүйе жүктелгенде қойылатын коэффициенттердің бастапқы мәндері бұл жобада регрессиялық талдау әдісі арқылы бағаланады. Бастапқы коэффициенттердің мәндері 6-шы кестеде келтірілген.
Басқару нысаны (2.1 бөлім) ретіндегі контактілік технологиялық процесі квазиорнықты нысандар жіктелігіне жатқызылады. Олардың орнықтылығы белгілі бір уақыт интервалымен (шамамен 10 минут) шектелген. Осыған контактілікті оптималдау мәселесін шешуді әрбір сатысында Кагмаж алгоритміне сәйкес үлгілерін бейімдеу керек. Кагмаж алгоритмі:
a j [nСH 1]СQ a j [n]СP
( y Э [n]СP y Р [n]) х Эj [n]
m
2
iСQ1
,
Э Р
есептелген мәндері;
xjэ - нысанның кіріс тәжірибелік айнымалылары;
Z - т ұрақты шама.
m
iСQ1
Бейімделу алгоритмі жұмысқа орнықтылық периоды соңында, жүйенің қызмет етуінің жалпы алгоритмін оптималдау блогын жүктер алдында қосылады. Төмендегі шарт орындалғанда ғана бейімделу алгоритмінің жұмысы аяқталады.
a j [nСH 1]СP a j [n]СOС .
Жоғарыда келтірілген (2.3.2 бөлімінде) контактілік аппараты бойымен газдарды оптималдау, тарату мәселесі теңдіктер мен теңсіздіктер типтегі шектеулі шартты экстремум есебі болып табылады. Осыған контактілікті бұл есепті шешу үшін процедура ұсынылған. Ол арқылы есеп оптималдау критерий өрнегіне штраф функциялары енгізу шартты емес экстремум есебіне әкелінеді. Мұнда ТВСО-1 оперативті жадыны үнемдеу мақсатында шешімнің эвристикалық процедурасы ұсынылған. Оның блок-схемасы төменде келтірілген (2.3 сурет). Блок 1. Таратушы ішкешен үшін ағымдағы ақпараттарды жинау жүзеге асырылады. Блок 2. Төмендегі формула бойынша металлургиялық газдардың жалпы саны анықталады.
nQ1
мұндағы Qn - электрпештерден (nСQ 1,2) , конвертордан (nСQ 3,6) , күкірт жағушы бөлімнен (n=7) газ шығындары. Блок 3. Төмендегі формула бойынша SO2 арқылы газдардың жалпы саны анықталады.
so2 P nQ1
мұндағы C n - электрпештерден (nСQ 1,2) , конверторлардан (nСQ 3,6) , күкірт жағатын бөлімдерден (n=7) шығатын газдардағы SO 2 концентрациясы.
Блок 4. j-ші контактілік аппаратындағы газ шығыны Qj жүзеге асырылады.
Блок 5. j-ші контактілік аппаратының і-ші қабатының шығысындағы температура есептеледі.
TkiСQ a0iСH a1iСH THiСH a2iСH CiСH a3iСH Qi
мұндағы a0 , a1 , a2 , a3 - математикалық үлгінің параметрлері.
Блок 6. j-ші контактілік аппаратының тотығу дәрежесі есептеледі.
iСQ1
мұндағыСLTi - і-ші қабаттың кірісі мен шығысындағы температулар айырмасы.
Блок 7. SO2 арқылы j-ші контактілік аппаратының өнімділігі анықталады. jСQ Q jСI C jСIСM j .
Блок 8. Төмендегі формула бойынша келесі агрегаттың нөмірі анықталады.
j = j+1.
Блок 9. Барлық агрегаттардың сұралғандығы талданады.
Блок 10. SO2 бойынша j-ші контактілік аппаратының өнімділіктерінің қосындысын есептеледі.
jСQ1
Блок 11. Газдардағы SO2 жалпы санынан SO2 өнімділіктері қосындысының ауытқуы есептеледі.
Блок 12. Есептелген ауытқу мәні қандай да бір шамамен (η) салыстырылады. Егер ауытқу мәні η шамадан аз болса, онда 13-ші блокқа өтеді, әйтпесе 14-ші блокқа өтеді.
Блок 13. Шектеулерді тексеру орындалады.
Q jСO Q jСO Q j
TKСO TKСO TK
Блок 14. Оптималдау критерий өрнегіне штраф функцияларын енгізу жүзеге асырылады.
мұндағыС1 ,С 2 ,С 3 ,С 4 - коэффициенттер.
Блок 15. Контактілік аппараттары бойындағы шығындардың оптималдау Q j мәндері шығарылады.
Басы Ақпарат жинау Шығарылатын газдың жалпы мөлшерін анықтау Qz SO2 бойынша шығарылатын газдың жалпы мөлшерін анықтау
j-ші БА бойындағы газ шығынын анықтау
Qj j-ші БА j-ші қабатының шығысындағы температураны есептеу
j-ші БА тотығу дәрежесін анықтау M j-ші БА өнімділігін анықтау
