Формат файлу SIR (версія SIC) - Інтегроване програмне забезпечення для моделювання каналів SIC ^ 2 та їх
Цей файл ASCII створюється EDISIR або за допомогою будь-якого редактора файлів ASCII. Є певні розділи, якими керує SIRENE, але не повністю EDISIR. EDISIR може читати та писати їх, але не вводити на інтерфейсах. Тому ці елементи повинен вводити редактор файлів ASCII. Він містить гідравлічні дані для програми перерахунку перехідних процесів SIRENE. Цей файл буде безпосередньо призначений для SIRENE, немає проміжного двійкового еквівалента на відміну від обчислення стаціонарного стану (DON двійковий еквівалент файлу ASCII FLU).
Цей файл призначений для старих версій SIC 4.33 або попередніх версій. Ці файли також можна імпортувати з версії SIC 5 та пізніших версій (у режимі інтерфейсу, але не в пакетному режимі). Крім того, якщо у вас є старі проекти SIC версії 4, ви можете захотіти зрозуміти цей формат для введення відповідних даних у проект SIC версії 5. Саме з цієї причини цей формат, який застарів, тим не менше описується. Попередні версії версії 5 одночасно використовували файли цього формату для супроводу переходу з версії 4 на версію 5, але це вже не так.
Спеціальний порядок написання цього файлу даних не дотримується. Різні типи можливих ліній описані нижче, кожен рядок починається з певної емблеми. З іншого боку, для імпортування даних до об’єкта цей об’єкт повинен бути вже визначений (наприклад: клапан повинен бути визначений перед введенням закону відкриття як функція часу).
Цей файл містить всю гідравлічну інформацію, необхідну для перехідного розрахунку програми SIRENE. Геометричні дані містяться в окремих файлах, таких самих, як і для стаціонарного стану. Коли одиниця не вказана, використовуються одиниці SI.
Для файлів версії 4 всі закони функцій часу розуміються у відносному вираженні, тобто якщо існує різниця між початковим значенням у початковій лінії води та значенням на початковий момент цього закону, тоді Закон буде зміщений цією різницею, щоб мати безперервність між початковим режимом і перехідним процесом, який потрібно імітувати.
Для файлів версії 5 такий спосіб обробки законів в обчисленні змінився, і ця зміна законів більше не виконується, тоді закони вважаються абсолютними. Це справедливо для всіх законів, за винятком законів побічних витрат як функції часу, які завжди компенсуються можливим початковим відхиленням, якщо таке є. Для версії 5 редактори даних управляють будь-якими автоматичними регулюваннями законів порівняно з початковими значеннями.
Назва_навчання
TITM, VERFIC, IFLAGX, IFLAGQ
ФОРМАТ ('!', A40, 1X, F5.3, 1X, I1, 1X, I1)
TITM назва дослідження,
Версія файлу VERFIC (3.0, 3.1, 4.0, 4.32, 4.33 тощо),
IFLAGX довгий режим координат x (0 якщо ні, 1 якщо так, тобто довгі координати x будуть помножені на 1000),
IFLAGQ режим високого потоку (0 якщо ні, 1 якщо так, тобто витрата буде помножена на 1000)
Досягає і Вузли
IB, N1, N2, TITB
ФОРМАТ ('N', 3I3, A40)
Номер досягнення IB,
Номер N1 висхідного вузла досяжності IB,
Номер вузла N2 досяжності IB за течією,
TITB назва охоплення
З версії 5 формат модифікується таким чином:
ФОРМАТ (‘N’, 3I3, A40, 2I5, 2F10.5)
з тими ж даними, що і вище, та 4 новими даними:
IFBPC (I, 1), IFBPC (I, 2), PCAM (I), PCAV (I)
IFBPC (I, 1) = 1, якщо вихідні навантаження рівні, = 0, якщо шанси рівні,
IFBPC (I, 2) = 1, якщо навантаження рівні нижче, = 0, якщо номінали рівні,
Перепад тиску PCAM (I) на вході в діапазон,
Перепад тиску PCAV (I) на виході з досяжності
Параметри розрахунку 1
IFLTI, IREP, TREP, IDHA, NEWT, NQUAS, IDISCR
ФОРМАТ ("B", 1X, I1, 1X, I1, 1X, F10.0, 1X, I1,1X, I1,1X, I3,1X, I1)
IFLTI flag_terme_inertie (1: з класичним терміном інерції, 2: локальний потік, 0 глобальний потік),
IREP flag_continuation_fichiers_res/rcs (0: ні),
TREP recovery_time (теплий старт) за хв,
Індекс автоматичного перепаду тиску IDHA для розширень (0: ні, 1: так), цей прапор передається від того, що вказаний у файлі .flu,
NEWT = 0, якщо фіксована точка, 1 якщо Newton, 2 якщо Fixed point + Newton,
NQUAS = 0, якщо чистий Ньютон, N, якщо Квазі-Ньютон: обчислення частоти похідної,
Індекс дискретизації IDISCR (0: класичний, 1: однорідний до постійного)
З версії 5 формат модифікується таким чином:
IFLTI, IREP, TREP, IDHA, INISV, NLSV, IDISCR
INISV = 0 класичне значення ініціалізації попереднього кроку,
INISV = 1 варіант ініціалізації попереднього кроку,
NLSV = 0 лінійний розрахунок,
Нелінійний розрахунок NLSV> 0 дає оновлення якобіана всі ітерації NLSV
Параметри розрахунку 2
IPS, IFR, IRCS, DT, TDEB, TFIN, XX, WW, XLAMB, ITER1, ERZ0, ITER2, EPSI
ФОРМАТ ('C', I3, 1X, 2 (I1, 1X), F10.2, 1X, 2 (F10.0, 1X), 2 (F4.2, 1X), F8.6, 1X, 2 (I4, 1X, F8.6))
IPS pas_de_stockage_resultat (кожні 1 або n не обчислюються, щоб обмежити розмір .RES-файлів),
IFR flag_de_regulation (0: без використання регулюючого файлу .reg, 1: with),
Прапор IRCS повного файлу .RCS (0: без створення цього файлу, 1: зі створенням цього файлу, необхідного для подальшого гарячого перезапуску, тобто перехідного обчислення з іншого перехідного обчислення),
DT no_time (у хвилинах),
TDEB час_ початку розрахунку (у хвилинах),
TFIN time_end розрахунку (у хвилинах),
XX верхній_ліміт ефективних_потоків_заборів,
WW нижчий_ліміт ефективного_потоку_до_розеток,
Лямбда-коефіцієнт контрактності XLAMB (0 = 4,17),
IDNS = зсув розділу для Z2 (корисно для клапанів AVIS, AVIO та змішаних клапанів) у кількості секцій розрахунку (Версії> = 4.17),
IFOL = рівняння (0: Гуссар, 1: Гуссар для положення та Cemagref для потоку, 2: W, заданий користувачем, та Cemagref для потоку),
IFOC = тип (0: Amil, 1: Avio Haute Chute, 2: Avio BC, 3: Avis HC, 4: Avis BC),
IFLO = модель,
ISF = режим розрахунку (0: нелінійний, 1: лінійний).
Назви вузлів, сума дебетів та розміри на вузлах
TC (NOMAX), ITN (NOMAX), QAE (NOMAX), ZNO (NOMAX), I
ФОРМАТ ('S', A6, I3, 2F10.5, 7X, I3)
TC_name_name (буквено-цифрове),
Тип вузла ITN = 1, якщо вища_модель, 2 якщо проміжна_вузол, -2 якщо нижча_модель,
QAE sum_debits = (Somme_Q_partant - Somme_Q_arrivant) на вузлі, що не стоїть за моделем, або швидкість витоку на вузлі aval_model,
ZNO сума_котів,
Я номер вузла
Розетки визначаються їх геометричними характеристиками
I, KP (I), KVP (I), QPO (I), QPR (I), X1P (I), X2P (I), X3P (I), X4P (I), X5P (I), X6P (I ), X7P (I), X8P (I), X9P (I), IFPL (I)
ФОРМАТ ('P', 3I3, 2F7.3, F6.2, 3F10.5, 3F6.2, F7.3, F6.2, I2)
I вузол (не вище_моделі та не_аваль_модель)
(KP) cas_prise =
1 Q (t) розрахунок відкриття клапана W
2 Розрахунок перемички Q (t) Z_seuil
3 Розрахунок перемички Q (t) width_threshold
4 розрахунок клапана Q
5 розрахування площини Z (t) Q
6 апм Q розрахунок
7 канальне числення Q
8 Q (z) розрахунок Q
9 плата L (t) обчислення Q
(KVP) cas_aval_prise =
1 Z_aval_ постійне утримання
2 deversoir_aval_prise
3 Q (z) _aval_prise
закони Q (z) на захопленнях вводяться в аналітичній формі:
Q = Q_base * ((Z-Z_seuil)/(Z_base-Z_seuil)) ** показник
(QPO) debit_prise_ob Об’єктив
(QPR) debit_prise_reel
(X1P) коефіцієнт_дебет (KP = 1, 2, 3, 4, 5, 6, 7, 9) або показник (KP = 8)
(X2P) Z_threshold
(X3P) поріг_ширини (KP = 1, 2, 3, 4, 5, 6, 9) або діаметр (KP = 1, 4, 7) або Q_base (KP = 8)
(X4P) відкриття (KP = 1, 4, 6) або Z_seuil_mini (KP = 2) або Z_base (KP = 8)
(X5P) максимальне_відкриття (KP = 1, 4 та X5P> 0)
або вказує на круговий клапан (KP = 1, 4 та X5P = 0) або Z_seuil_maxi (KP = 2)
або максимальна_ширина_порогового значення (KP = 3) або перемикач (KP = 6, 7, 8)
(X6P) коефіцієнт_дебіт_aval_prise (KVP = 2) або показник_aval_prise (KVP = 3)
(X7P) Z_aval_prise (KVP = 1) або Z_seuil_aval_prise (KVP = 2, 3)
(X8P) width_seuil_aval_prise (KVP = 2) або Q_base_aval_prise (KVP = 3)
(X9P) Z_base_aval_prise (KVP = 3)
IFPL = рівняння (0: Cemagref або 1: Cunge)
Закони Q (z) на модельному низхідному вузлі
I, Z, Q
ФОРМАТ ("A", I3, 6X, 2F10.3)
Я node_aval,
Z cote_aval (м),
Q debit_aval (м3/с)
Налаштування дисплея
КОДЕМ, БРАНМ, НАЗВИ, НБМ, НСМ, XLONG, ZREF
ФОРМАТ (‘Y’, A1, A8, A6, 2I3, F8.0, F8.2) у версії 4.32
ФОРМАТ (‘Y’, A1, A8, A6, I3, I4, F8.0, F8.2) у версії 4.33
КОДЕМ-код,
ІМЯ_ФІЛІЇ BRANM,
NODE ім'я_вузла,
Охоплення НБМ,
NSM number_relative_section,
XLONG абсциса,
ZREF cote_refer
Варіації відкриття клапана
IB, IS, KS, IO, T, X, K, L, ITYPOUV
ФОРМАТ ('U', 4I3, 2F10.3,3I3)
Діапазон IB,
Є відносний_секційний_ номер,
KS num._singular._relatif,
Тип конструкції IO: 1 (водостік), 2 (прямокутний клапан) та 3 (клапан Gec-Alsthom) з версії 4.07 включено (до версії 4.07 у клапанах керували лише клапанами)
T час (хв),
Х відкриття,
K номер закону,
L номер закону,
ITYPOUV (0: крок, 1: пандус)
Гідрограф або закон, прийнятий на вузлі моделі, що не стоїть за потоком
CH1, J, T, X, I, L, IOD, NOE, ITYPOFF
ФОРМАТ ("H", A1, 2X, I3, 3X, 2F10.3, I3, I5, 3I3)
CH1 char_loi_prise: 'Q' flow_ outlet, 'W' opening_ outlet, 'S' socket_switch, 'L' outlet_threshold, 'Z' Z_seuil_ outlet,
J noeud_non_aval,
T час (хв),
X закон_значення,
Я номер закону,
L номер пункту закону,
IOD (I) тип штекера,
Режим розрахунку NOE (I, 3) при зйомці,
ITYPOFF (0: крок, 1: пандус)
Цільові витрати
J, T, Q, I, L
ФОРМАТ ('O', I3, 6X, 2F10.3, 2I3)
J вузол,
T час (хв),
Q дебет_об'єктив,
Я номер закону,
L номер пункту закону
Закон Локера (S (z) у вузлі)
I, Z, S, X
ФОРМАТ (‘K’, I3, 6X, 3F10.3)
Я вузол,
Z розмірність,
S площа (га),
X швидкість випаровування (см/добу)
Бічний гідрограф (керований Сірене, але не Едісіром)
Якщо початковий загальний бічний потік, зазначений у цьому бічному гідрографі для цього досяжності, не дорівнює початковому загальному побічному потоку, що надходить із розрахунку стаціонарного стану, гідрограф буде переведений для забезпечення безперервності при запуску між стійким станом та перехідним процесом режим. Це справедливо для версій 4 SIC, але також і для версії 5.
З іншого боку, наразі неможливо ввести термін інфільтрації - випаровування як функцію протягу. Це може бути змодельовано приблизним чином шляхом послідовного відведення в режимі розрахунку потоку. Нам здається, це обмеження є дещо обмежувальним, оскільки оцінка терміну інфільтрації - випаровування дуже делікатне з дуже важливим фактором невизначеності, і ілюзорно розглядати закон q = f (z) дасть більш реалістичний.
Цей розділ "бічний гідрограф як функція часу" незабаром буде повністю керований Edisir (введення в графічні інтерфейси) у пізнішій версії. Edisir вже дозволяє читати та записувати його у файл, але не заповнювати інтерфейсами (починаючи з версії 3.8). Sirene це чудово справляється з точки зору розрахунків.
IB, K, T, Q
ФОРМАТ (‘L’, 2I3, 3X, 2F10.3)
Номер досягнення IB,
K прапор кроку (0) або рампи (1)
T час (хв),
Q дебетовий_бічний (м3/с)
Лімніграм на модельному низхідному вузлі (керований Sirene, але не Edisir)
Цей розділ дозволяє замінити закон Q (z), введений для обчислення в стаціонарному стані, законом z (t), що дозволяє, наприклад, представити вплив припливу та хвилі.
Незабаром ним керуватиме Edisir у пізнішій версії. Edisir вже дозволяє читати та записувати його у файл, але не заповнювати інтерфейсами.
Значення z обробляються відносно порівняно з першим значенням. Наприклад, якщо ви хочете збільшити розмірність z на 1м за мить 2h = 120mn, все, що вам потрібно зробити, це ввести:
Z 2 0 0,000 0,000
Z 2 0 120 000 1000
Попередження: якщо ви вводите закон z (t), вам слід пам’ятати про видалення закону Q (z), оскільки в іншому випадку він буде дотриманий, оскільки він має пріоритет над z (t). Закон Q (z) також має пріоритет над кошиком, який закривав би канал. Якщо цей бін ввести з інтерфейсів, то Q (z) буде правильно видалено з файлу .SIR. З іншого боку, якщо цей контейнер вводиться вручну безпосередньо у файл .SIR, не забудьте видалити Q (z), щоб він більше не використовувався.
Обчислення з z (t) набагато швидше, ніж з Q (z), оскільки не існує нелінійних ітерацій для обчислень у нижчому за рівнем вузлі. З іншого боку, цей нижній стан z (t) іноді вважається "жорстким" гідравлічно і може генерувати коливання потоків у цій зоні, наприклад у випадку, коли гідрограф, що стоїть вище за течією, надходить на вузол нижче за течією з az (t ) постійний. Цей нижній стан z (t) насправді не дуже фізичний, оскільки "поверхневого пінцета" не існує. Приплив фізичний, але це не просте встановлення рейтингу, а процес, який також пов’язує потоки та рейтинги через складні фізичні процеси. У випадку синусоїдального типу за течією z (t), що імітує приплив, апріорі паразитних коливань немає.
Навіть якщо це можливо на рівні файлу .SIR, не бажано мати як bin, так і z (t) для одного і того ж низхідного вузла. У цьому випадку можна було б накласти потік, а також розмірність, яка створює числові проблеми. У пізнішій версії цей режим буде явно заборонений, і в цьому випадку саме z (t) буде ефективним, а запис не працює.
ФОРМАТ (‘Z’, 2I3, 3X, 2F10.3)
I номер node_aval,
Прапор K кроку (0) або рампи (1) [опція з версії 3.7]
T час (хв),
Z cote_aval