Добавил:
Upload Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

MV_Mat_Modeli_2014

.pdf
Скачиваний:
22
Добавлен:
09.02.2016
Размер:
2.42 Mб
Скачать

31

До протоколу необхідно додати роздруківку графіка кривих рівного виходу, який побудовано згідно рівняння регресії.

Застосування пакета програми PLAN для складання багатофакторних планів експериментів. Програма багатофакторного планування експериментів PLAN оформлена у вигляді виконуваного файлу plan.exe, файлу вихідних даних isx.txt і файлу результатів rez.txt. Крім цього, на кафедрі ТЗЗ розроблена спеціальна програма побудови кривих рівного виходу за результатами роботи програми PLAN. Програма побудови кривих рівного виходу складається з файлу початкових даних dgr-rez.txt, виконуваного файлу graf3cga.exe і декількох службових файлів.

3 СКЛАДАННЯ МАТЕМАТИЧНИХ МОДЕЛЕЙ ТИПОВИХ ПРОЦЕСІВ ГАЛУЗІ АНАЛІТИЧНИМИ І КОМБІНОВАНИМИ МЕТОДАМИ

3.1 МАТЕМАТИЧНА МОДЕЛЬ МЕХАНІЧНИХ ПРОЦЕСІВ

Л а б о р а т о р н а р о б о т а № 5

МАТЕМАТИЧНЕ МОДЕЛЮВАННЯ ПРОЦЕСУ ДОЗУВАННЯ КОМПОНЕНТІВ КОМБІКОРМІВ

М е т а р о б о т и

Вивчити застосування аналітичних методів для опису механічних процесів. Скласти математичний опис вагового дозатора комбікормового заводу

для визначення його продуктивності і витрат енергії на дозування.

З а в д а н н я н а п і д г о т о в к у д о л а б о р а т о р н о ї р о б о т и

1.Вивчити навчальний матеріал [1, стор. 206–227; 2, стор. 105–107; 4, стор. 43–53; 4, стор. 56–60] і теоретичний матеріал лекцій. Звернути увагу на:

загальний порядок складання математичних моделей аналітичним методом; залежності, застосовувані для опису процесу дозування.

2.Відобразити в теоретичній частині протоколу:

мету роботи; конструкцію шнекового живильника;

формули для визначення продуктивності установки і енерговитрат на дозування.

3.Визначити етапи складання моделі процесу дозування.

4.Підготувати модель для визначення продуктивності шнекового дозатора і витрат енергії на дозування.

5.Використовувати відомі аналітичні залежності для складання математичного опису процесу дозування.

6.Виконати за вихідними даними індивідуального завдання розрахунки із визначення продуктивності шнекового дозатора і витрат енергії на дозування.

7.Зробити перевірку розрахунків на ЕОМ.

32

Т е о р е т и ч н а ч а с т и н а

У багатьох випадках можна одержати математичний опис технологічного процесу, не вдаючись до експериментів, а лише ґрунтуючись на відомих теоретичних залежностях, тобто аналітичними методами. Якщо ж поряд з відомими теоретичними залежностями використовуються і дані експериментів, то говорять про комбіновані методи складання математичного опису і математичних моделей.

Механічні процеси здрібнювання, сепарування, дозування, змішування, гранулюванню широко застосовуються в харчовій промисловості. У більшості випадків ці процеси добре вивчені, тобто теоретичні залежності, що описують їх, відомі.

Складання математичного опису аналітичними методами здійснюють у декілька етапів:

1.Постановка задачі. Тут необхідно скласти формулювання завдання, указавши вид технологічного процесу, об'єкт моделювання (використовуване обладнання), а також призначення математичної моделі.

2.Вибір показників якості математичної моделі. Ці показники прямо визначаються призначенням моделі.

3.Визначення параметрів, що впливають на об'єкт.

4.З'ясування залежностей, що зв'язують вхідні фактори і вихідні параметри.

Розглянемо процес складання математичного опису аналітичними методами на прикладі процесу дозування, характерного для комбікормового виробництва.

Етап 1. Постанова задачі. Скласти математичний опис вагового дозатора (вузла дозування) комбікормового заводу для визначення його продуктивності і витрат енергії на дозування.

Етап 2. Вибір показників якості. Вибираємо наступні показники якості: Gy – продуктивність вузла дозування;

Е – енерговитрати на дозування.

Етап 3. Визначення параметрів, що впливають на об'єкт: nк – кількість компонентів у комбікормі;

mi – маса і-го компонента в дозованій порції; Gi – продуктивність і-го живильника;

М – вантажопідйомність вагів.

Етап 4. З'ясування залежностей, що зв'язують вхідні фактори і вихідні параметри (показники якості).

У прикладі об'єктом буде ваговий дозатор сипучих компонентів. Кількість компонентів nк = 3. Схема установки наведена на рис. 5.1.

Дозатор складається з вагів (1), трьох бункерів (2), трьох шнекових живильників (3). Живильники включаються по черзі і працюють визначений час, що

33

залежить від процентного вмісту дозуючого компонента в комбікормі і від продуктивності самого живильника. Кожен живильник працює по черзі в двох режимах: 90 % дозованого компонента відмірюється в основному режимі і 10 % – у режимі досипання (з меншою продуктивністю для збільшення точності дозування).

Після того, як усі живильники відпрацюють, виміряна порція сировини надходить на наступну операцію – змішування.

2

2

3 3

1

а)

2

2 1

б)

3

2

Таким чином, загальний час дозування визначається тривалістю змішування продукту (щоб не допускати перебоїв у технологічному процесі).

Час змішування складає 6 хв., тому загальний час дозування всіх трьох компонентів не повинен перевищувати 5 хв.

Рисунок 5.1– Схема вагового дозатора

а) вигляд збоку; б) вигляд зверху 1– ваги; 2 – бункер; 3 – шнековий живильник

Р о з р а х у н к о в і ф о р м у л и :

 

1. Продуктивність установки

 

Gу = 24 nц Mη, т/добу,

(5.1)

де nц – кількість циклів дозування за годину (nц = 10, тому що кількість циклів дозування визначається кількістю циклів змішування);

η = 0,9 – К.К.Д. використання вагів.

2. Час дозування і-го компонента

τц = Στi, (5.2) де τi – час дозування і-го компонента (повний час дозування, що складається з часу роботи живильника в основному режимі τi1 та в режимі досипання

τi2);

τі = τi1 + τi2;

(5.3)

 

 

34

 

 

τI1

=

60 0,9 MI

;

(5.4)

 

1000 GI1

 

 

 

 

 

τI2

=

 

60 0,1 MI

,

(5.5)

 

1000 GI2

 

 

 

 

 

де 0,9 і 0,1 – коефіцієнти, що враховують процентне відношення маси до-

зованого компонента в основному режимі та у режимі досипання (90% і 10%, відповідно).

3. Маса і-го компонента

 

 

 

 

mI =

aI

M

η ,

(5.6)

 

 

 

 

 

100

 

 

де аі – процентний вміст і-го компонента в рецепті комбікорму.

 

4. Продуктивність шнека

 

 

 

 

G = Qvγ,

(5.7)

де Qv – об'ємна продуктивність шнека, м3/год;

 

γ – об'ємна маса сировини, т/м3. Приймають:

 

– для зернової сировини γ = 0,5 т/м3;

 

– для борошнистої сировини (висівки) γ = 0,3 т/м3;

 

– для мінеральної сировини γ = 1,2 т/м3.

 

Qv = 3600 F v, м3/год,

(5.8)

де F – площа перерізу шнека, який зайнятий сировиною, м2;

 

v – швидкість переміщення сировини шнеком, м/c.

 

F =

πD2

ψ ,

(5.9)

 

 

 

 

 

 

 

 

 

4

 

 

 

 

 

 

де D – діаметр шнека, м;

 

 

 

 

ψ = 0,4 — коефіцієнт заповнення шнека.

 

V =

S Nв

, м/c,

(5.10)

 

 

 

60

 

 

 

де S – крок гвинта шнека, м. Приймають S = D;

 

nв – частота обертання вала шнека.

 

Nв =

Nч

 

,

 

 

(5.11)

I

 

 

 

 

 

 

 

 

 

де nч – частота обертання вала електродвигуна;

і – передаточне число редуктора (приймають і = 30). Частоту nч приймають:

для основного режиму nч = 2800 хв.–1,

для режиму досипання nэ = 1420 хв.–1. Конструкція живильника показана на рис. 5.2.

Діаметр шнека D залежить від його марки. У дозаторах даного типу вста-

новлюють шнеки ПШ–200, ПШ-320 чи ПШ–400, де число означає діаметр шнека в міліметрах. Для кожного компонента необхідно вибрати свій живильник, з огляду на то, що чим діаметр більший, тим вища продуктивність дозатора і тем менший час і точність дозування. Рекомендують для мінеральної сировини ви-

35

користовувати дозатор ПШ-200, а для зернової і мучнистої – ПШ-320 або ПШ400.

Бункер для сировини

А – А

А

Сировина

S

Редуктор

Привідний

електродвигун

L

А

 

 

 

Ваговий бункер

 

 

 

 

 

 

 

 

 

 

Рисунок 5.2 – Конструкція живильника

6. Енерговитрати на дозування Е, кВт год.

3

 

E = ( NI1 τI1 +NI2 τI2 ) ,

(5.12)

I=1

 

де Ni – потужність, споживана електродвигуном у даному режимі, кВт; τі1 – час дозування в основному режимі, г; τі2 – час дозування в режимі досипання, г.

 

Gi

 

,

(5.13)

Nij =

 

G L K

м

 

 

3600

 

 

 

де g = 9,8 м2/с – прискорення вільного падіння;

L – довжина шнека (відстань від центра випускного отвору бункера до центра випускного отвору шнека у ваговий бункер), м;

Км = 1,2 – коефіцієнт надлишку потужності.

Величину Nij окремо розраховують для двох режимів дозування.

Зернова сировина є основною, тому для неї потрібно передбачити великий бункер (2×2 м чи 3×3 м – розміри брати з завдання).

Розміри бункерів для борошнистої і мінеральної сировини у всіх варіантах однакові і складають 1×1 м. Величина L повинна бути більше половини ширини бункера і кратної 0,5 м.

П и т а н н я д о с а м о к о н т р о л ю

1.Наведіть основні фізико-механічні закономірності, використовувані для встановлення математичних зв'язків між вхідними і вихідними параметрами об'єкта, який модулюється.

2.Наведіть приклади використання математичних моделей у харчовій технології.

3.Основні положення системного підходу до складання моделі механічних процесів.

36

4. Складання математичного опису процесу дозування на комбікормовому заводі аналітичним методом.

П о р я д о к в и к о н а н н я л а б о р а т о р н о ї р о б о т и

1.Відповісти викладачу на питання передлабораторного контролю.

2.У ході лабораторної роботи з вихідних даних відповідно до варіанта скласти модель процесу дозування компонентів комбікорму.

3.Провести аналітичний розрахунок продуктивності дозатора і витрат електроенергії на процес дозування.

4.Перевірити правильність розрахунків на ЕОМ.

П р а к т и ч н а ч а с т и н а

Умова задачі

Скласти математичний опис аналітичним методом процесу дозування для визначення добової продуктивності вагового дозатора і витрат електроенергії на процес дозування.

В а р і а н т и і н д и в і д у а л ь н и х з а в д а н ь

Вихідні дані вибирають з Додатка Д. За останньою цифрою номера залікової книжки визначають вантажопідйомність вагів, розміри бункера для зернових компонентів. На перетинанні рядків, обумовлених за останньою і передостанньою цифрами номера залікової книжки, визначають вміст компонентів (зернового, мінерального й висівок) у рецептурі комбікорму. Вміст висівок (В) визначають за формулою (5.14).

В = 100 – (З + М), % (5.14) де З – зернові компоненти; М – мінеральні компоненти

Порядок складання математичного опису аналітичним методом роботи вагового дозатора

1.Вибрати вихідні параметри.

2.Вибрати вхідні параметри.

3.Скласти параметричну схему.

4.На підставі відомих вихідних даних визначити Gy , формула (5.1).

5.Знайти масу кожного компонента mi у дозованій порції, формула (5.6).

6.Знайти площу перерізу шнека, зайнятого сировиною, формула (5.9).

7.Визначити частоту обертання вала шнека для основного режиму і режиму досипання, швидкість переміщення сировини шнеком, формули (5.11), (5.10).

8.Визначити об'ємну продуктивність шнека для кожного i-того компонента, формула (5.8).

9.Продуктивність шнека для кожного i-того компонента, формула (5.7). 10.Час дозування i-того компонента, (5.4) , (5.5), (5.3) і (5.2). 11.Потужність, споживану електродвигуном у даному режимі для i-того

компонента, формула (5.13).

12.Загальні витрати електроенергії на дозування, формула (5.12). 13.Перевірка правильності опису вагового дозатора на ЕОМ.

37

3 . 2 Т Е П Л О М А С О О Б М І Н Н І П Р О Ц Е С И

Л а б о р а т о р н а р о б о т а № 6

МАТЕМАТИЧНИЙ ОПИС МОДЕЛЕЙ ТЕПЛОМАСООБМІННИХ ПРОЦЕСІВ КОМБІНОВАНИМИ МЕТОДАМИ

М е т а р о б о т и

Вивчити застосування комбінованих методів для опису тепломасообмінних процесів. Відповідно до індивідуального завдання скласти математичний опис шахтної прямоточної зерносушарки аналітичним методом.

З а в д а н н я н а п і д г о т о в к у д о л а б о р а т о р н о ї р о б о т и

1.Повторити навчальний матеріал [6, стор.64–76].

2.Вивчити навчальний матеріал [1, стор.295–320; 3, стор.142–146; 5, стор. 84–103] і теоретичний матеріал лекцій. Звернути увагу на:

загальний порядок складання математичного опису аналітичним методом процесу сушіння в шахтній прямоточній зерносушарці;

залежності, застосовувані для опису тепломасообмінних процесів, що характеризують роботу шахтної зерносушарки і фізичний зміст величин, що входять у ці залежності.

2.Відобразити в теоретичній частині протоколу: мету роботи;

підготувати модель тепловологопереносу в шахтній прямоточній зерносушарці.

3.Визначити етапи складання моделі процесу сушіння.

4.Використовувати відомі аналітичні залежності для складання математичного опису процесу сушіння.

5.Виконати за вихідними даними індивідуального завдання розрахунки із визначення витрати сушильного агента і витрат умовного і натурального палива на сушіння зерна.

6.Зробити перевірку розрахунків на ЕОМ.

Т е о р е т и ч н а ч а с т и н а

Комбіновані методи математичного опису застосовують при недостатніх знаннях про об'єкт моделювання. У цьому випадку основну частину математичного опису складають аналітичним способом, а незначна частина - заснована на експериментальних даних.

У харчовій промисловості тепломасообмінні процеси мають величезне значення. До цих процесів відносять процеси сушіння, випарювання, ректифікації й інші. Загальним для всіх цих процесів є наявність двох станів чи фаз компонентів і перенесення компонентів з однієї фази в іншу. Перенесення речовин здійснюється спочатку в межах однієї фази (внутрішнє перенесення), потім речовина переборює поверхню розділу фаз (зовнішнє перенесення) і розподіляється в межах другої фази. Усередині кожної фази речовини переносяться, головним чином, молекулярною дифузією і конвективним перенесенням. Кожна фаза може складатися з одного чи декількох компонентів.

38

Розглянемо комбінований метод математичного опису стосовно до технологічного процесу сушіння зерна в шахтній прямоточній сушарці.

Сушіння є типовим тепломасообмінним процесом, при якому тепло переміщається від сушильного агента до матеріалу, що висушується, а волога — навпаки, від матеріалу, що висушується, до сушильного агента. Сушінню піддають різні харчові продукти, напівфабрикати і сировину (зерно, солод, кормові і харчові дріжджі, крохмаль і ін.). Процес сушіння обумовлений підведенням теплоти до продукту, що висушується, завдяки якому відбувається випаровування вологи. Як сушильний агент застосовують повітря, перегріту пару і топкові гази, що насичуються вологою, яка випаровується висушуваного матеріалу.

Складання математичного опису шахтної прямоточної зерносушарки здійснюють у декілька етапів.

Етап 1 – Постанова задачі.

Необхідно скласти формулювання завдання, указавши вид технологічного процесу, об'єкт моделювання (використовуване устаткування), а також призначення математичної моделі.

Потрібно скласти математичний опис шахтної прямоточної зерносушарки для визначення витрати сушильного агента і витрат палива на сушіння.

Етап 2 – Вибір показників якості.

Показниками якості будуть:

V – об'ємна витрата сушильного агента, м3/год.; ВТ – масова витрата палива, кг/год.

Етап 3 – Визначення вхідних факторів, що впливають на об'єкт.

Ці фактори можна розділити на декілька груп.

1. Фактори, що характеризують атмосферне повітря: t0 – температура, °С;

ϕ0 – відносна вологість, %.

2. Фактори, що характеризують об'єкт сушіння — сире зерно: G1 – витрата сирого зерна на вході в сушильну камеру, т/год.; ω1 – вологість зерна на вході в сушильну камеру, %.

3. Фактори, що характеризують сушильний агент: t1 – температура, °С.

Етап 4 – З'ясування залежностей, що зв'язують вхідні фактори і показники якості моделі.

Для опису процесу сушіння зерна в шахтній зерносушарці можна скористатися декількома рівняннями. У першу чергу це рівняння матеріального і теплового балансів.

1. Для визначення об'ємної витрати сушильного агента складемо рівняння матеріального балансу сушильної камери по волозі (рис. 6.1).

На рис. 6.1 позначені:

G1, G2 – витрата зерна на вході і виході сушильної камери відповідно, кг/год; ω1 і ω2 – вологість зерна на вході і виході сушильної камери відповідно, %; L – масова витрата сухого повітря в складі сушильного агента, кг/год.;

d1 і d2 – вологовміст сушильного агента, г/кгсух.в.

39

 

 

 

Волога сирого зерна

G1W1

 

 

 

100

 

 

 

 

 

 

 

 

 

 

 

Волога сушильного агента

 

 

 

 

Волога відпрацьованого

 

 

 

 

сушильного агента

 

Сушильна

 

 

 

D1

 

 

 

 

 

 

 

 

D2

L

 

камера

 

 

 

 

 

 

 

 

L

 

 

 

 

1000

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1000

 

Волога просушеного зерна G2ω2

100

Рисунок 6.1 – Матеріальний баланс сушильної камери за вологістю

Рівняння матеріального балансу сушильної камери по волозі має вигляд

G1

 

ω1

+ L

D1

 

= G2

 

ω2

+ L

 

D2

,

(6.1)

100

 

 

100

1000

 

1000

 

 

 

 

Перетворивши рівняння (6.1), одержимо

G1

 

ω1

− G2

 

ω2

= L

D2 D1

.

100

100

 

 

 

1000

 

Ліва частина рівняння представляє вологу, яка випарена із зерна в процесі сушіння. Позначимо цю вологу через W, кг/год

W = L

D2 D1

.

(6.2)

 

1000

 

 

З рівняння (6.2) можна визначити масові витрати сушильного агента L

L =

1000 W

.

(6.3)

 

d2 − d1

Убагатьох випадках зручно використовувати питомі величини, віднесені

до 1 кілограма випаруваної вологи. Позначимо питому масову витрату сушильного агента як l, кгсух. повітря /кгвип. вологи

l =

L

=

1000

.

(6.4)

 

 

Wd2 − d1

Таким чином, знаючи вологовміст сушильного агента на вході сушильної камери і на її виході, можна обчислити питомі масові витрати сушильного агента. Однак такий параметр, як масові витрати сушильного агента на практиці не застосовується, тому що його вимірювання викликає великі труднощі. Набагато легше піддається контролю інший параметр – об'ємні витрати сушильного агента. Тому слід перейти від масових витрат L до об'ємних витрат V. Для того,

щоб здійснити такий перехід, введемо поняття питомого об'єму vпит

 

 

V = L Vпит .

(6.5)

Значення питомого об'єму обчислюють за формулою

 

vпит =

 

287 (273 + t1 )

,

(6.6)

 

 

 

 

B −

Pп

 

 

 

 

 

 

 

100

 

 

 

де t1 – температура повітря чи сушильного агента, °С;

40

В – атмосферний тиск, Па; РП – парціальний тиск водяної пари, що міститься у повітрі су-

шильного агента, Па; ϕ – відносна вологість повітря чи сушильного агента, %.

Оскільки нас цікавлять витрати сушильного агента, то потрібно підставляти у формулу (6.6) параметри сушильного агента.

Значення В в лабораторній роботі приймають рівним 760 мм рт. ст.

Тиск, виражений у міліметрах ртутного стовпчика і тиск, виражений у Паскалях, зв'язані наступним співвідношенням

PП = 133,322 Pмм рт.ст.

(6.7)

Визначення парціального тиску водяної пари РП у сушильному агенті здійснюється по-різному, в залежності від температури сушильного агента t1:

при t1 > 100 °C РП = В, а при t1 > 100 °C

РП = РН,

 

де РН – парціальний тиск насиченої пари при тій же температурі.

Парціальний тиск насиченої водяної пари в сушильному

агенті РН,

мм рт.ст., можна визначити за формулою

 

 

lg(Pн ) = 0,622 +

7,5 T1

 

(6.8)

 

 

 

 

 

238 + T1

 

 

 

 

 

 

 

Відносну вологість сушильного агента ϕ1 , % визначають з рівняння, що

зв'язує відносну вологість повітря і його вологовміст

 

d1 = 6,22

ϕ1Pн

.

 

(6.9)

 

 

 

 

В −

ϕ1Pн

 

 

 

 

 

 

 

 

100

 

 

 

 

 

З рівняння (6.9) для сушильного агента

 

 

ϕ1 = 100

 

 

 

D1B

.

(6.10)

 

Pн (622 + D1 )

 

 

 

 

 

Таким чином, для визначення об'ємних витрат сушильного агента V залишилося знайти W, d1 та d2.

З рівняння матеріального балансу сушильної камери можна показати, що

 

W = G1 − G2 ,

(6.11)

де

G2

= G1

100 − ω1

.

(6.12)

 

 

 

 

100 − ω2

 

Вологість зерна на виході із сушильної камери ω2 у лабораторній роботі приймають рівною 14 %, а вологовміст сушильного агента d1 – рівним вологовмісту атмосферного повітря d0.

2. Для визначення вологовмісту відпрацьованого сушильного агента d2 і витрат палива на сушіння ВТ складемо рівняння теплового балансу сушильної камери (рис. 6.2).

На рис. 6.2 позначені: с1 і с2 – питома теплоємність зерна на вході і виході сушильної камери відповідно, кДж/кг К;Θ1 і Θ2 – температура зерна на вході і виході сушильної камери відповідно, °С; I1 і I2 – питома ентальпія сушильного агента, кДж/кг.

Соседние файлы в предмете [НЕСОРТИРОВАННОЕ]