Отчет по лабораторной работе Моделирование гидравлической системы в нестационарном режиме.

Российский химико-технологический университет

имени Д. И. Менделеева

Кафедра информатики и компьютерного моделирования

Отчет по лабораторной работе № 2

Моделирование гидравлической системы в нестационарном режиме .

Вариант № 14

Выполнил: студент гр. П-45

Золотарев Д. С.

Проверил: Тамбовцев И. И.

Москва 2010

Схема гидравлической системы:

Основные допущения, накладываемые на модель:

1) Система работает в нестационарном режиме;

2) Рассматривается только гидравлика процесса;

3) Сопротивление трения потоку пренебрежимо мало;

4) Жидкость несжимаема и имеет постоянную плотность;

5) Газ - идеальный pV=nRT;

6) Режим - изотермический T=const.

Система уравнений математического описания:

A) Уравнения расходов жидкости через клапаны:

1. v₁ = k₁ (p₁ - p₇ )^1/2

2. v₂ = k₂ (p₂ – p₈ ) ^1/2

3. v₃ = k₃ (p₇ - p₃ ) ^1/2

4. v₄ = k₄ (p₈ - p₄ ) ^1/2

5. v₅ = k₅ (p₈ - p₅ ) ^1/2

6. v₆ = k₆ (p₈ - p₆ ) ^1/2

7. v₇ = k₇ (p₈ – p₇ ) ^1/2

B) Определение давлений жидкости и газа в емкостях:

8. p₇ = p₉ + ρgh₁

9. p₉ = P_н H₁ / (H₁ - h₁ )

10. p₈ = p₁₀ + ρgh₂

11. p₁₀ = P_н H₂ / (H₂ - h₂ )

C) Дифференциальные уравнения заполнения емкостей:

12. dh₁ /dt = (v₁ - v₃ ­ – v₇ ) / s₁

13. dh₂ /dt = (v₂ – v₄ – v₅ – v_{6 +} v ₇ ) / s₂

D) Начальные условия:

14. h₁ (t=0) = h₁ ⁰

15. h₂ (t=0) = h₂ ⁰

Информационная матрица системы уравнений:

N п/п	v1	v2	v3	v4	v5	v6	v7	p7	p8	p9	p10	h1	h0 1	h2	h0 2
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15

Программа

Option Explicit

Option Base 1

Public culc As Integer

Const m As Byte = 2 ' Количество решаемых ОДУ

Const np% = 10, nk% = 7, nv% = 7, g! = 9.8

Public del! ' шаг

Private k!

Private y0!(m), x0! ' Начальное значение

Private n% ' Число шагов

Private kv% ' Кратность вывода

Private x1!, y1!(m), i%

Private ki!, ki1!, ki2!, vm!(nk), v!(nk), ak!(nk), p!(np)

Private hg!(m), h!(m), s!(m), pr!(m), ro!, pn!, x!, ipr%, ipr1%

Sub dydx(x As Single, y() As Single, pr!())

' расчет производных для текущих(pr) x,y

Dim j%

For j = 1 To m: h(j) = y(j): Next j

p(9) = pn * hg(1) / (hg(1) - h(1)): p(10) = pn * hg(2) / (hg(2) - h(2))

p(7) = p(9) + ro * g * h(1) * 0.000001: p(8) = p(10) + ro * g * h(2) * 0.000001

v(1) = ak(1) * Sgn(p(1) - p(7)) * Sqr(Abs(p(1) - p(7)))

v(2) = ak(2) * Sgn(p(2) - p(8)) * Sqr(Abs(p(2) - p(8)))

v(3) = ak(3) * Sgn(p(7) - p(3)) * Sqr(Abs(p(7) - p(3)))

v(4) = ak(4) * Sgn(p(8) - p(4)) * Sqr(Abs(p(8) - p(4)))

v(5) = ak(5) * Sgn(p(8) - p(5)) * Sqr(Abs(p(8) - p(5)))

v(6) = ak(6) * Sgn(p(8) - p(6)) * Sqr(Abs(p(8) - p(6)))

v(7) = ak(7) * Sgn(p(8) - p(7)) * Sqr(Abs(p(8) - p(7)))

pr(1) = (v(1) - v(3) + v(7)) / s(1): pr(2) = (-v(4) - v(5) - v(6) + v(2) - v(7)) / s(2)

For j = 1 To nk: vm(j) = v(j) * ro: Next j

With Worksheets("TMP")

' TMP out

'.Cells(ipr, 1) = i: .Cells(ipr, 2) = "p(5-7)": .Cells(ipr, 6) = "vm"

'.Cells(ipr, 3) = p(5): .Cells(ipr, 4) = p(6): .Cells(ipr, 5) = p(7)

'.Cells(ipr, 7) = vm(1): .Cells(ipr, 8) = vm(2): .Cells(ipr, 9) = vm(3)

'.Cells(ipr, 10) = vm(4): .Cells(ipr, 11) = vm(5): ipr = ipr + 1

'.Cells(ipr, 2) = "x": .Cells(ipr, 4) = "y(1-2)": .Cells(ipr, 7) = "pr(1-2)"

'.Cells(ipr, 3) = x: .Cells(ipr, 5) = y(1): .Cells(ipr, 6) = y(2)

'.Cells(ipr, 8) = pr(1):

Внимание, отключите Adblock

Вы посетили наш сайт со включенным блокировщиком рекламы!
Ссылка для скачивания станет доступной сразу после отключения Adblock!

Скачать