АКЦИОНЕРНОЕ ОБЩЕСТВО ЗАКРЫТОГО ТИПА
НАУЧНО-ПРОИЗВОДСТВЕННАЯ ФИРМА «ЦКБА»
РУКОВОДЯЩИЙ
ДОКУМЕНТ
МЕТОДИКА РАСЧЁТА ТЕМПЕРАТУРНЫХ ПОЛЕЙ
ТРУБОПРОВОДНОЙ АРМАТУРЫ
РД 26-07-25-97
РУКОВОДЯЩИЙ ДОКУМЕНТ
Методика
расчёта температурных полей |
РД 26-07-25-97 Вводится впервые |
Дата введения 4 мая 1998 года
Настоящий руководящий документ распространяется на запорную, регулирующую и предохранительную арматуру, предназначенную для установки в системах атомных электростанций, опытных и исследовательских ядерных реакторах, в системах и установках химических, нефтехимических и других производств с максимальным рабочим давлением до 20 МПа (200 кгс/см2) и с температурой рабочей среды до 873 К (600 °С) и устанавливает методику проведения теплового расчёта арматуры DN 10 до 1200 мм включительно.
Председатель Технического Комитета М.И. Власов «05» ____02______ 1998 г. |
УТВЕРЖДАЮ Генеральный директор В.А. Айриев «___» ___________ 199 г. |
ЛИСТ УТВЕРЖДЕНИЯ
Руководящего документа «Методика расчёта температурных полей трубопроводной арматуры».
РД 26-07-25-97
Заместитель Генерального
директора АОЗТ НПФ «ЦКБА»
по научной работе Ю.И. Тарасьев
Главный конструктор В.В. Ширяев
Заместитель директора НТИЦ Р.И. Хасанов
Начальник отдела 161 А.А. Косарев
Начальник лаборатории 154 В.В. Никитин
Начальник сектора Г.И. Сергевнина
Ведущий инженер-исследователь Н.С. Косых
СОГЛАСОВАНОРуководитель инспекции __________ В.А. Лавров «17» _января__ 1997 г. |
СОГЛАСОВАНО
Подпись А.А. Шаталов (Письмо № 11-11/307 от 9.09.98 г.) |
В настоящем руководящем материале использованы ссылки на следующие нормативные документы:
РД 302-07-122-89. «Методика расчёта температурных полей арматуры для сред с температурой до 873 К (600 °С)».
СТП 07.81-634-96 «Методические указания по проведению теплового расчёта задвижек для АЭС».
СТП 07.81-631-92. «Методика проведения теплового расчёта высокопараметрической арматуры клапанного типа».
3.1. Тепловой расчёт арматуры должен проводиться с целью определения температурных полей, необходимых для оценки работоспособности её теплонапряжённых узлов и деталей, а также выбора конструкционных материалов с ограниченной теплостойкостью.
3.2. Методика теплового расчёта должна решать линейные и нелинейные задачи стационарной и нестационарной теплопроводности при заданных граничных условиях для Кmax = 300 узловых точек.
3.3. Тепловой расчёт должен включать в себя:
а) задачу расчёта;
б) исходные данные;
в) условия расчёта;
г) тепловую модель;
д) программу теплового расчёта на электронных вычислительных машинах (ЭВМ);
е) расчет;
ж) заключение.
3.4. Программа теплового расчёта составлена на алгоритмическом языке «фортран» и реализована на персональных ЭВМ типа IВМ РС/AТ.
Требования к РС:
- наличие оперативной памяти не менее 536 КБ;
- наличие дискового пространства для хранения информации до 2 МБ.
3.5. Исходные данные и условия расчёта методики изложены в разделах 4, 5 настоящего стандарта.
4.1. Для проведения теплового расчёта необходимы следующие данные:
а) температурное поле на начальный момент времени t = 0, К;
б) закон изменения температуры рабочей и окружающей сред (линейный, нелинейный) или скорость изменения температуры рабочей и окружающих сред, b, К/ч;
в) Тн - температура рабочей среды на начальный момент времени t = 0, К;
г) Тк - температура рабочей среды на конец нестационарного режима, К;
д) Та - температура окружающей среды, К;
е) геометрические размеры расчетного изделия, м;
ж) марки материалов деталей;
и) теплофизические характеристики:
l - коэффициент теплопроводности материалов деталей, вт/м×ч×К;
a0 - коэффициент теплоотдачи со стороны рабочей среды, вт/м2×К,
aа - коэффициент теплоотдачи от наружных поверхностей изделия в окружающую среду, вт/м2×К,
r - плотность материалов, кг/м3;
с - удельная теплоёмкость материалов, вт/кг×К;
л) e1 - абсолютная невязка, К;
м) e2 - относительная невязка, %;
н) w - коэффициент релаксации.
5.1. При проведении теплового расчета за начальные условия следует принять температурное поле на начальный момент времени t = 0.
5.2. Граничные условия на наружных и внутренних поверхностях изделия (соприкасающихся с рабочей и окружающей средами) следует задавать в виде температуры рабочей и окружающей сред и соответствующих коэффициентов теплоотдачи (g = а×Dt).
5.3. Контактные термические сопротивления между соприкасающимися деталями не учитываются.
5.4. При проведении теплового расчёта следует рассматривать осесимметричную двухмерную нестационарную задачу в цилиндрической системе координат.
5.5. Тепловой расчёт следует проводить при переменных теплофизических характеристиках материалов (l, с, r) на каждом шаге по времени Dt.
6.1. В основу методики теплового расчёта положено уравнение теплового баланса. Решение данного уравнения осуществлено методом конечных разностей (МКР).
6.2. Тепловая модель
6.2.1. Тепловая модель рассчитываемого изделия основана на дискретном представлении элементарных участков изделия с помощью схем замещения из термических сопротивлений.
6.2.2. В основе тепловой модели (см. рис. 1) принята осесимметричная половина сечения по оси штока (оси патрубков - в зависимости от того, какую часть изделия необходимо рассчитать) рассчитываемого изделия.
6.2.3. Разбивку на элементарные участки следует производить таким образом, чтобы избежать значительных искажений геометрии изделия.
6.2.4. Центральным точкам элементарных участков тепловой модели присваивается порядковый номер.
6.2.5. Верхняя и нижняя границы тепловой модели определяются исходя из конструктивных особенностей рассчитываемого изделия, а также в зависимости от того, температурное поле какого узла изделия интересует конструктора (заказчика). Проточная часть корпуса в тепловую модель выступающей части изделия не включается. При необходимости определения температурного поля нижней части корпуса (зоны проточной части) составляется отдельная тепловая модель. Принцип её составления аналогичен рассмотренному в настоящем разделе.
6.2.6. На нижней границе тепловой модели тепловой поток в осевом направлении задаётся равный нулю (Q = 0).
Тепловая модель клапана запорного DN 100
1 - фланец корпуса; 2 - шпилька; 3 - фланец крышки; 4 - втулка сильфонной сборки.
6.2.7. При составлении тепловой модели фланцевого соединения изделия шпильки и гайки заменяются эквивалентными цилиндрами, имеющими общую с фланцевым соединением ось. Воздушный зазор между шпилькой и верхним фланцем не учитывается.
6.2.8. Эквивалентный цилиндр гайки выбирается таким образом, чтобы площадь его поперечного сечения и площадь боковой поверхности были равны соответственно суммарной площади поперечных сечений и суммарной площади боковых поверхностей гаек фланцевого соединения.
6.2.9. Шпильки заменяются тремя эквивалентными цилиндрами, расположенными:
а) в зоне верхнего фланца;
б) между нижним и верхним фланцами;
в) в зоне нижнего фланца.
6.2.10. Эквивалентные цилиндры шпилек выбираются таким образом, чтобы площади поперечных сечений и площади их боковых поверхностей были равны соответственно суммарной площади поперечных сечений и суммарной площади боковых поверхностей шпилек фланцевого соединения.
6.2.11. Тепловая модель должна включать внешние узловые точки, предназначенные для задания граничных условий. Внешним узловым точкам присваиваются номера 301, 302 и т.д.
6.3. Описание программы теплового расчёта на ЭВМ
6.3.1. Расчётная тепловая схема
6.3.1.1. В основу расчетной тепловой схемы для определения температурного поля элементарного участка положено конечно-разностное пятиточечное уравнение теплового баланса, решённое относительно температуры в узловой точке для момента времени t:
(6.1)
где Qm-i, Qn-i, Qp-i, Qq-i - тепловые потоки между рассчитываемым участком и соседними участками m, n, p, q, соответствующие моменту времени t;
ti - температура в узловой точке i в момент времени t;
- температура в узловой точке i в момент времени t - Dt;
Dt - шаг по времени.
6.3.1.2. Расчёт температуры в узловой точке каждого рассчитываемого элементарного участка в момент времени следует производить по формуле:
(6.2)
где i = 1, 2, 3 ... k - номера узловых точек;
ci - удельная теплоёмкость материала участка, Вт/кг×К;
- суммарная термическая проводимость соседних элементов;
Rm-i, Rn-i, Rp-i, Rq-i - термические сопротивления между узловыми точками, К/Вт.
6.3.1.3. Число уравнений должно равняться числу внутренних узлов тепловой модели.
6.3.1.4. Для каждого рассчитываемого участка следует определять термические сопротивления в радиальном и осевом направлениях, а при наличии соседних наружных узловых точек следует определять термические сопротивления теплообмена.
6.3.1.5. Расчёт термических сопротивлений в радиальном направлении для каждой узловой точки тепловой модели следует производить по формуле:
(6.3)
где l - коэффициент теплопроводности материала участка, вт/м;
hz - размер участка в осевом направлении, м;
ri - внутренний радиус участка, м;
ri+1 - внешний радиус участка, м.
6.3.1.6. Расчёт термических сопротивлений в осевом направлении для каждой узловой точки тепловой модели следует производить по формуле:
(6.4)
где hz, l, ri+1, ri - величины, определённые в п. 6.3.1.5.
6.3.1.7. Расчёт термических сопротивлений теплообмена для наружных узловых точек следует производить по формуле:
(6.5)
где a - коэффициент теплоотдачи со стороны рабочей и окружающей сред, вт/м2×К,
S - площадь элементарного участка, омываемого рабочей или окружающей средами, м2.
6.3.1.8. Расчёт коэффициента теплоотдачи от торцевых и боковых поверхностей изделия в окружающую среду следует выполнять по методике, изложенной в РД 302-07-122-89.
6.3.1.9. Расчёт коэффициента теплоотдачи от рабочей среды к внутренним поверхностям проточной части и горловины изделия следует выполнять по методикам, изложенным в:
а) СТП 07.81-631-92 - для арматуры клапанного типа;
б) СТП 07.81-634-96 - для задвижек.
6.3.1.10. Геометрические характеристики эквивалентных цилиндров, заменяющих гайки и шпильки, определяются следующим образом:
а) расчёт расстояния от узловой точки эквивалентного цилиндра гаек до оси симметрии определяется по формуле:
(6.6)
где z - число гаек;
r2 - радиус гайки, м.
б) расчёт половины толщины эквивалентного цилиндра гаек производится по формуле:
(6.7)
в) расстояние от узловой точки эквивалентного цилиндра части шпилек, расположенной в зонах верхнего и нижнего фланцев до оси симметрии, принимается равным радиусу окружности центров отверстий под шпильки до оси симметрии;
г) расчёт половины толщины эквивалентного цилиндра части шпилек, расположенной в зоне верхнего и нижнего фланцев, производится по формуле:
(6.8)
где z - число шпилек;
rш1,2 - соответственно радиусы шпилек в зоне верхнего и нижнего фланцев, м.
д) расчёт половины толщины эквивалентного цилиндра части шпилек, расположенной между верхним и нижним фланцами, производится по формуле:
(6.9)
где rш2 - радиус шпильки в зоне между верхним и нижним фланцами, м.
е) расчёт расстояния от узловой точки эквивалентного цилиндра части шпилек, расположенной между верхним и нижним фланцами, производится по формуле:
(6.10)
6.3.2. Программа теплового расчёта трубопроводной арматуры, предназначенная для реализации на ЭВМ, включает в себя:
а) алгоритм теплового расчёта изделия;
б) общие сведения о программе;
в) логическую структуру;
г) подготовку и описание входных данных;
д) рабочие массивы и вспомогательные переменные;
е) подготовку и описание выходных данных.
6.3.3. Алгоритм теплового расчёта изделия
6.3.3.1. Основное расчётное уравнение (6.2) для определения на ЭВМ температуры в узловой точке имеет вид:
ti = Amtm + Antn + Aptp + Aqtq + ... + Aiti-Dt (6.11)
где Аm, Аn, Аp, Aq, Ai - коэффициенты уравнения:
(6.12)
6.3.3.2. Решение системы уравнений (6.12) следует производить релаксационным методом Зейделя, т.е. решение на каждом шаге по времени считается законченным, если температурное поле в следующем приближении не отличается от предыдущего на величину более заданной.
6.3.3.3. При описании программы теплового расчёта сетка узловых точек имеет вид, представленный на рис. 2, где i = 1, 2, 3 ... k - номер узловой точки;
N1, N2, N3, N4 - номера соседних узловых точек;
- расстояние от узловой точки до границы участка в радиальном направлении;
RT - расстояние от узловой точки до оси симметрии;
DH = hz - расстояние от узловой точки до границы участка в осевом направлении;
ri+1, ri - внешний и внутренний радиусы участка.
6.3.4. Общие сведения о программе
Программа содержит следующую информацию:
а) используемые операторы DIMENSION, READ, FORMAT, PRINT, DO, IF, GOTO, CONTINUE, STOP, END;
б) используемые стандартные программы FORT, LKED, ALOG.
6.3.5. Логическая структура программы
6.3.5.1. Логическая структура программы теплового расчёта содержит:
а) единую блок-схему теплового расчёта;
б) порядок расположения блоков в программе;
в) описание блоков, составляющих программу.
6.3.5.2. Единая блок-схема теплового расчёта представлена на рис. 3.
6.3.5.3. Порядок расположения блоков в программе представлен на рис. 4.
6.3.5.4. Каждый блок, составляющий программу теплового расчёта, содержит следующую информацию:
а) блок 1, представленный на рис. 5, содержит операторы, приведенные в п. 6.3.4. В блоке 1 необходимо описать все массивы, которые встречаются в программе расчёта. Исходное температурное поле образуется путём засылки начального значения температуры рабочей среды во все ячейки массива температуры. Помимо ввода исходных данных и формирования температурного поля в блоке 1 производится печать исходных данных;
б) блок 2, представленный на рис. 6, служит для определения граничных условий, необходимых для расчёта температурного поля в данный момент времени. При определении граничных условий используется линейная интерполяция температуры среды и коэффициентов теплоотдачи между заданными значениями для различных моментов времени;
в) блок 3, представленный на рис. 7, предназначен для расчёта коэффициентов основного уравнения Ат, Ап, Ар, Аq, Аi на каждом шаге времени. При расчёте коэффициентов Ат, Ап, Ар, Аq, Аi используется температурное поле, полученное на предыдущем шаге времени. В начале линейным интерполированием определяются значения li, Сi, затем значения Rm-i, Rn-i, Rp-i, Rq-i и выполняется расчёт величин , Ат, Ап, Ар, Аq, Аi. Все указанные в этом блоке расчёты производятся в цикле (для каждой расчётной точки);
г) блок 4, представленный на рис. 8, предназначен для решения системы уравнений релаксационным методом Зейделя. Решение считается законченным, если температурное поле в следующем приближении не отличается от предыдущего на величину более заданной (см. п. 6.3.3.2).
6.3.6. Подготовка и описание входных данных
6.3.6.1. Исходными данными для расчёта являются:
а) геометрия элементарных участков тепловой модели рассчитываемого изделия;
б) граничные условия расчёта;
в) теплофизические характеристики материалов деталей и элементов, составляющих изделие;
г) начальные условия расчёта.
6.3.6.2. Ввод информации осуществляется в текстовые файлы средствами текстового редактора, загруженного в ПЭВМ (LEXICON):
а) в файле «geom.dat» (приложение А) задаются:
- общее количество элементарных участков тепловой модели;
- номера сред, окружающих изделие (рабочая и внешняя);
- количество среднемассовых точек температурного поля;
- геометрические характеристики элементарных участков тепловой модели;
- номера материалов элементарных участков тепловой модели (в соответствии с номерами из файла «tx01.dat»);
- номера элементарных участков тепловой модели, среднемассовые температуры которых необходимо определить;
б) в файле «wr01.dat» (приложение Б) задаются:
- параметры временного интервала проведения теплового расчёта (шаги времени);
- количество итераций расчёта;
- скорость изменения температуры рабочей среды;
- граничные условия на внутренних и наружных поверхностях изделия. Задаются в виде температуры рабочей и окружающей сред, коэффициентов теплоотдачи на поверхностях, омываемых рабочей и окружающей средами;
в) в файле «tx01.dat» (приложение В) задаются:
- данные по теплофизическим характеристикам материалов, применяемым в рассчитываемой конструкции трубопроводной арматуры. Значения теплофизических характеристик приводятся при ряде температур. Каждому материалу, введённому в данный файл, соответствует присвоенный номер;
- количество применяемых материалов;
г) в файле «tn01.dat» (приложение Г) содержится таблица массива начальных температур для узлов тепловой модели;
д) в файл «osn1.dat» (приложение Д) вводится набор текстовых файлов с входной информацией: «geom01.dat», «wr01.dat», «tx01.dat», «tn01.dat» для конкретного теплового расчёта методом конечных разностей. Ввод осуществляется средствами редактора «Norton».
Блок-схема теплового расчёта
Файл «osn1.dat» с входной информацией зафиксирован в программе расчёта - файл «filt.for» (приложение Е).
6.3.6.3. Обозначения исходных данных, принятые в тексте программы, представлены в таблице 1.
Обозначение в программе |
Наименование величины |
|
- |
KN |
Характер расчёта (начало или продолжение) |
- |
KVD |
Количество различных шагов по времени |
- |
KM |
Количество материалов в изделии |
- |
NG1 |
Количество внутренних узлов в тепловой модели |
- |
ND |
Номер первичного узла граничной области доп. тепловой модели |
- |
NG2 |
Номер последнего узла граничной области доп. тепловой модели |
- |
N |
Номер последнего узла внутренней области доп. тепловой модели |
- |
NH |
Номер первой точки среды |
- |
L |
Номер последней точки среды |
- |
MM |
Максимальное число итераций |
e1 |
DT1, C |
Абсолютная невязка |
e2 |
DT2 |
Относительная невязка |
w |
RK |
Коэффициент релаксации |
m, n, p, q |
N1 [I:NG1], N2 [I:NG1], N3 [I:NG1], N4 [I:NG1], NI [ND:N] |
Массивы номеров соседних участков (узлов) |
ri + (ri+1 - ri)/2 |
RT [I:NG1], RT [ND:N], м DR [I:NG1], DR [ND:N], м |
Массив радиусов узлов Массив расстояний от узловых точек до границ участков в радиальном направлении |
Dt |
D [I:KVD], c |
Массив шагов по времени |
KD [I:KVD] |
Массив количества шагов |
|
- |
MM [I:NG1], MM [ND:N] |
Массив номеров участков материалов |
- |
TG1, TG2, TG3, TG4, K |
Значения температуры, при которых задаются l, r, с |
l |
TR1 [I:KM], TR2 [I:KM], TR3 [I:KM], TR4 [I:KM], Вт/м×К |
Массивы теплопроводности материалов при различных температурах |
r × c |
РR1 [I:KM], РR2 [I:KM], РR3 [I:KM], РTR4 [I:KM], Вт/м3×К |
Массивы произведения плотности материалов на их удельную теплоёмкость при различных температурах |
- |
МКТ [I:KTS] |
Массив количества линейных участков изменения температуры среды [КТS = L - NH + 1] |
- |
MKK [I:KTS] |
Массив количества линейных участков изменения коэфф. теплоотдачи |
- |
VTS [I:KTS, I:K], c |
Двухмерный массив начальных моментов времени, соответствующих линейным участкам изменения температуры среды [К = МК] |
t0, ta |
T [I:KTS, I:KT], K |
Двухмерный массив температуры среды |
- |
VKT [I:KT, I:KK], c |
Двухмерный массив начальных моментов времени, соответствующих линейным участкам изменения коэффициентов теплоотдачи [КК = МКК] |
a0, aa |
TK [I:KTS, I:KK], Вт/м2×К |
Двухмерный массив коэффициентов теплоотдачи |
- |
T [I:NG1], T [ND:N], K |
Начальное температурное поле |
NT1 [I:KTG], NT2 [I:KTG] |
Массивы номеров узлов основной тепловой модели для определения температуры на границе дополнительной тепловой модели [КТG = NG2 - ND + 1] |
Примечание. При отсутствии дополнительной тепловой модели ND = NG2 = N массивы N1, N2, N3, N4, RT, DR, DН, ММ, Т должны содержать информацию только об основной тепловой модели. Массивы Т1, Т2 не вводятся.
6.3.6.4. Рабочие массивы и вспомогательные переменные, принятые в тексте программы, приведены в таблице 2.
6.3.6.5. Трансляция программы в счёт осуществляется запуском командного файла «с.bat».
6.3.6.6. Выходные данные выдаются ПЭВМ в виде машинограммы, включающей в себя:
а) исходные данные;
б) результаты вычислений. Результаты расчёта записываются в файл «niki.txt» (приложение Ж).
6.3.6.7. Выходные величины приведены в таблице 3.
6.3.6.8. Пример оформления теплового расчета приведен в приложении И.
Обозначение в программе |
Наименование величины |
|
Am, An, Ap, Aq, Ai |
A1 [1:NG1], A2 [1:NG1], A3 [1:NG1], A4 [1:NG1], A5 [1:NG1], A1 [ND:N] |
Массивы коэффициентов системы уравнений |
t, ti |
TH [1:NG1], TH [ND:N], K |
Температурное поле предыдущего шага по времени |
li |
TP [1:NG1], TP [ND:N], Вт/м×K |
Массив теплопроводности материалов участков при данном температурном поле |
a0, aa |
TP Вт/м2×К |
Массив коэффициентов теплоотдачи в момент времени t |
ts |
SAM (I), K |
Массив средних температур в узлах |
tm, tn, tp, tq |
T1, T2, T3, T4, K |
Температуры соседних узлов |
Rm-i, Rn-i, Rp-i, Rq-i |
R1, R2, R3, R4, K/ккал |
Термические сопротивления между соседними узлами |
ti |
TT, K |
Температура узла в момент времени t |
ri |
M |
Номер итерации |
Pt |
MN |
Количество несовпадений |
- |
1I, 1K, K, KK, TD, JD, 1I, 1KK, DT, DTG, DTP, DPK, DVR, KDV, DTP, DPK, DT, DTG, A, B, C, Ka, P1, P2, P3, S, SMMAX1, SMMAX2 |
Вспомогательные переменные |
ci |
CT [1:NG1], T [ND:N], Вт/кг×К |
Массив теплоёмкостей участков при данном температурном поле |
Обозначение в программе |
Наименование величин |
|
t |
VRM, VR, мин |
Время |
t |
T [I:NG], T [ND:N], K |
Температурное поле |
t0, ta |
T [NH:L], K |
Температура среды |
ПРИЛОЖЕНИЯ
NG1 - >085< - количество внутр. узлов основной тепловой модели
ND - >085< - номер первого узла граничной области доп. тепловой модели
NG2 - >085< - номер последнего узла граничной области доп. тепловой модели
N - >085< - номер последнего узла внутр. области допол. тепловой модели
NH - >301< - номер первой точки среды
L - >302< - номер последней точки среды
KNN - >004< - количество среднемассовых объектов (KNN > 0)
N1 |
N2 |
N3 |
N4 |
RT |
DR |
DH |
NM |
|
1 |
301 |
2 |
6 |
0 |
0,060500 |
0,003500 |
0,075000 |
03 |
2 |
1 |
3 |
7 |
0 |
0,069000 |
0,005000 |
0,075000 |
03 |
3 |
2 |
4 |
8 |
302 |
0,078500 |
0,004500 |
0,075000 |
03 |
4 |
3 |
5 |
9 |
302 |
0,085000 |
0,002000 |
0,075000 |
03 |
5 |
4 |
302 |
10 |
302 |
0,096000 |
0,009000 |
0,075000 |
03 |
6 |
301 |
7 |
22 |
1 |
0,060500 |
0,003500 |
0,075000 |
03 |
7 |
6 |
8 |
23 |
2 |
0,069000 |
0,005000 |
0,075000 |
03 |
8 |
7 |
9 |
24 |
3 |
0,078500 |
0,004500 |
0,075000 |
03 |
9 |
8 |
10 |
25 |
4 |
0,085000 |
0,002000 |
0,075000 |
03 |
10 |
9 |
302 |
26 |
5 |
0,096000 |
0,009000 |
0,075000 |
03 |
11 |
0 |
12 |
14 |
0 |
0,011000 |
0,011000 |
0,015000 |
01 |
12 |
11 |
13 |
15 |
301 |
0,025500 |
0,003500 |
0,015000 |
05 |
13 |
12 |
301 |
16 |
301 |
0,032500 |
0,003500 |
0,015000 |
05 |
14 |
0 |
15 |
17 |
11 |
0,011000 |
0,011000 |
0,015000 |
01 |
15 |
14 |
16 |
18 |
12 |
0,025500 |
0,003500 |
0,015000 |
05 |
16 |
15 |
301 |
19 |
13 |
0,032500 |
0,003500 |
0,015000 |
05 |
17 |
0 |
18 |
29 |
14 |
0,011000 |
0,011000 |
0,002500 |
01 |
18 |
17 |
19 |
30 |
15 |
0,025500 |
0,003500 |
0,002500 |
05 |
19 |
18 |
20 |
31 |
16 |
0,032500 |
0,003500 |
0,002500 |
05 |
20 |
19 |
21 |
32 |
301 |
0,041000 |
0,005000 |
0,002500 |
01 |
21 |
20 |
22 |
33 |
301 |
0,051500 |
0,005500 |
0,002500 |
01 |
22 |
21 |
23 |
34 |
6 |
0,060500 |
0,003500 |
0,002500 |
03 |
23 |
22 |
24 |
302 |
7 |
0,069000 |
0,005000 |
0,002500 |
03 |
24 |
23 |
25 |
302 |
8 |
0,078500 |
0,004500 |
0,002500 |
03 |
25 |
24 |
26 |
27 |
9 |
0,085000 |
0,002000 |
0,002500 |
03 |
26 |
25 |
302 |
302 |
10 |
0,096000 |
0,009000 |
0,002500 |
03 |
27 |
302 |
302 |
28 |
25 |
0,085000 |
0,002000 |
0,008500 |
02 |
28 |
302 |
302 |
49 |
27 |
0,085000 |
0,002000 |
0,008500 |
02 |
29 |
0 |
30 |
35 |
17 |
0,011000 |
0,011000 |
0,008500 |
01 |
30 |
29 |
31 |
36 |
18 |
0,025500 |
0,003500 |
0,008500 |
05 |
31 |
30 |
32 |
37 |
19 |
0,032500 |
0,003500 |
0,008500 |
05 |
32 |
31 |
33 |
38 |
20 |
0,041000 |
0,005000 |
0,008500 |
03 |
33 |
32 |
34 |
39 |
21 |
0,051500 |
0,005500 |
0,008500 |
03 |
34 |
33 |
302 |
40 |
22 |
0,060500 |
0,003500 |
0,008500 |
03 |
35 |
0 |
36 |
41 |
29 |
0,011000 |
0,011000 |
0,008500 |
01 |
36 |
35 |
37 |
42 |
30 |
0,025500 |
0,003500 |
0,008500 |
05 |
37 |
36 |
38 |
43 |
31 |
0,032500 |
0,003500 |
0,008500 |
05 |
38 |
37 |
39 |
44 |
32 |
0,041000 |
0,005000 |
0,008500 |
05 |
39 |
38 |
40 |
45 |
33 |
0,051500 |
0,005500 |
0,008500 |
05 |
40 |
39 |
302 |
46 |
34 |
0,060500 |
0,003500 |
0,008500 |
05 |
41 |
0 |
42 |
51 |
35 |
0,011000 |
0,011000 |
0,010000 |
01 |
42 |
41 |
43 |
52 |
36 |
0,025500 |
0,003500 |
0,010000 |
03 |
43 |
42 |
44 |
53 |
37 |
0,032500 |
0,003500 |
0,010000 |
03 |
44 |
43 |
45 |
54 |
38 |
0,041000 |
0,005000 |
0,010000 |
03 |
45 |
44 |
46 |
55 |
39 |
0,051500 |
0,005500 |
0,010000 |
03 |
46 |
45 |
47 |
56 |
40 |
0,060500 |
0,003500 |
0,010000 |
03 |
47 |
46 |
48 |
57 |
302 |
0,069000 |
0,005000 |
0,010000 |
03 |
48 |
47 |
49 |
58 |
302 |
0,078500 |
0,004500 |
0,010000 |
03 |
49 |
48 |
50 |
59 |
28 |
0,085000 |
0,002000 |
0,010000 |
02 |
50 |
49 |
302 |
60 |
302 |
0,096000 |
0,009000 |
0,010000 |
03 |
51 |
0 |
52 |
61 |
41 |
0,011000 |
0,011000 |
0,004500 |
01 |
52 |
51 |
53 |
62 |
42 |
0,025500 |
0,003500 |
0,004500 |
04 |
53 |
52 |
54 |
63 |
43 |
0,032500 |
0,003500 |
0,004500 |
04 |
54 |
53 |
55 |
64 |
44 |
0,041000 |
0,005000 |
0,004500 |
04 |
55 |
54 |
56 |
65 |
45 |
0,051500 |
0,005500 |
0,004500 |
03 |
56 |
55 |
57 |
66 |
46 |
0,060500 |
0,003500 |
0,004500 |
03 |
57 |
56 |
58 |
302 |
47 |
0,069000 |
0,005000 |
0,004500 |
03 |
58 |
57 |
59 |
302 |
48 |
0,078500 |
0,004500 |
0,004500 |
03 |
59 |
58 |
60 |
302 |
49 |
0,085000 |
0,002000 |
0,004500 |
02 |
60 |
59 |
302 |
302 |
50 |
0,096000 |
0,009000 |
0,004500 |
03 |
61 |
0 |
62 |
67 |
51 |
0,011000 |
0,011000 |
0,007000 |
01 |
62 |
61 |
63 |
68 |
52 |
0,025500 |
0,003500 |
0,007000 |
04 |
63 |
62 |
64 |
69 |
53 |
0,032500 |
0,003500 |
0,007000 |
04 |
64 |
63 |
65 |
70 |
54 |
0,041000 |
0,005000 |
0,007000 |
04 |
65 |
64 |
66 |
71 |
55 |
0,051500 |
0,005500 |
0,007000 |
03 |
66 |
65 |
302 |
72 |
56 |
0,060500 |
0,003500 |
0,007000 |
03 |
67 |
0 |
68 |
73 |
61 |
0,011000 |
0,011000 |
0,007000 |
01 |
68 |
67 |
69 |
74 |
62 |
0,025500 |
0,003500 |
0,007000 |
04 |
69 |
68 |
70 |
75 |
63 |
0,032500 |
0,003500 |
0,007000 |
04 |
70 |
69 |
71 |
76 |
64 |
0,041000 |
0,005000 |
0,007000 |
04 |
71 |
70 |
72 |
77 |
65 |
0,051500 |
0,005500 |
0,007000 |
03 |
72 |
71 |
302 |
78 |
66 |
0,060500 |
0,003500 |
0,007000 |
03 |
73 |
0 |
74 |
79 |
67 |
0,011000 |
0,011000 |
0,007500 |
01 |
74 |
73 |
75 |
80 |
68 |
0,025500 |
0,003500 |
0,007500 |
04 |
75 |
74 |
76 |
81 |
69 |
0,032500 |
0,003500 |
0,007500 |
04 |
76 |
75 |
77 |
82 |
70 |
0,041000 |
0,005000 |
0,007500 |
04 |
77 |
76 |
78 |
83 |
71 |
0,051500 |
0,005500 |
0,007500 |
03 |
78 |
77 |
302 |
84 |
72 |
0,060500 |
0,003500 |
0,007500 |
03 |
79 |
0 |
80 |
0 |
73 |
0,011000 |
0,011000 |
0,010000 |
01 |
80 |
79 |
81 |
0 |
74 |
0,025500 |
0,003500 |
0,010000 |
04 |
81 |
80 |
82 |
0 |
75 |
0,032500 |
0,003500 |
0,010000 |
04 |
82 |
81 |
83 |
0 |
76 |
0,041000 |
0,005000 |
0,010000 |
04 |
83 |
82 |
84 |
0 |
77 |
0,051500 |
0,005500 |
0,010000 |
03 |
84 |
83 |
85 |
0 |
78 |
0,060500 |
0,003500 |
0,010000 |
03 |
85 |
84 |
302 |
0 |
302 |
0,074500 |
0,010500 |
0,010000 |
03 |
NI |
NL |
КI |
К6 |
массив для определения среднемассовой температуры |
||||
1 |
01 |
22 |
26 |
05 |
||||
2 |
25 |
59 |
59 |
25 |
||||
3 |
46 |
56 |
60 |
50 |
||||
4 |
32 |
38 |
40 |
34 |
||||
NТ1 |
NТ2 |
массив узлов основной тепловой модели для определения температуры на границе дополнительной тепловой модели |
||||||
1 |
... |
... |
KVD => 22 < -- количество шагов по времени
ММ => 100 < -- количество итераций
DT1 => 0,100 < -- абсолютная невязка
DT2 => 0,0010 < -- относительная невязка
RK => 1,800 < -- коэффициент релаксации
D |
KD |
|
1 |
0,5 |
10 |
2 |
5,0 |
5 |
3 |
10,0 |
3 |
4 |
20,0 |
3 |
5 |
30,0 |
6 |
6 |
60,0 |
5 |
7 |
120,0 |
5 |
8 |
300,0 |
8 |
9 |
600,0 |
6 |
10 |
1800,0 |
6 |
11 |
3600,0 |
5 |
12 |
0,5 |
10 |
13 |
5,0 |
5 |
14 |
10,0 |
3 |
15 |
20,0 |
3 |
16 |
30,0 |
6 |
17 |
60,0 |
5 |
18 |
120,0 |
5 |
19 |
300,0 |
8 |
20 |
600,0 |
6 |
21 |
1800,0 |
6 |
22 |
3600,0 |
5 |
MKT |
MKK |
5 |
3 |
1 |
1 |
ВРЕМЯ |
Тсред |
0,0 |
20,0 |
10,0 |
190,0 |
36000,0 |
190,0 |
36010,0 |
20,0 |
72000,0 |
20,0 |
ВРЕМЯ |
АLFA |
0,0 |
1900,0000 |
36000,5 |
1100,0000 |
36010,0 |
1100,0000 |
ВРЕМЯ |
Тсред |
0,0 |
20,0 |
ВРЕМЯ |
АLFA |
0,0 |
7,0000 |
КЛАПАНА ЗАПОРНОГО DN 100 16.07.96 г. (крышка, корпус из ст. 20) исп. Сергевнина Г.И. |
(1 СТРОКА - ПУСТАЯ, 2 И 3 - ИНФОРМАЦИОННЫЕ)
*>05<* - количество используемых материалов
0000,0 |
0120,0 |
0240,0 |
0360,0 |
0000,00 |
0120,00 |
0240,00 |
0360,00 |
|
материал |
теплопроводность |
произведение р · с |
||||||
1 о8х18н10т |
11,340 |
13,150 |
14,940 |
10,640 |
862,00 |
942,00 |
984,00 |
1010,00 |
2 хн35вт |
10,600 |
12,800 |
14,900 |
16,600 |
1000,00 |
1005,00 |
1020,00 |
1040,00 |
3 сталь 20 |
43,970 |
41,490 |
39,010 |
36,510 |
742,00 |
885,00 |
990,00 |
1080,00 |
4 сталь 45 |
51,900 |
50,300 |
47,100 |
43,700 |
894,00 |
931,00 |
998,00 |
1076,00 |
5 бронза |
64,700 |
76,400 |
86,300 |
93,000 |
0600,00 |
0600,00 |
0600,00 |
0600,00 |
КN = >0< -- характер расчетов МАССИВ НАЧАЛЬНЫХ ТЕМПЕРАТУР |
||||||||||
..1 |
..2 |
..3 |
..4 |
..5 |
..6 |
..7 |
..8 |
..9 |
..0 |
|
Т(001 ... 010) |
020,0 |
020,0 |
020,0 |
020,0 |
020,0 |
020,0 |
020,0 |
020,0 |
020,0 |
020,0 |
Т(011 ... 020) |
020,0 |
020,0 |
020,0 |
020,0 |
020,0 |
020,0 |
020,0 |
020,0 |
020,0 |
020,0 |
Т(021 ... 030) |
020,0 |
020,0 |
020,0 |
020,0 |
020,0 |
020,0 |
020,0 |
020,0 |
020,0 |
020,0 |
Т(031 ... 040) |
020,0 |
020,0 |
020,0 |
020,0 |
020,0 |
020,0 |
020,0 |
020,0 |
020,0 |
020,0 |
Т(041 ... 050) |
020,0 |
020,0 |
020,0 |
020,0 |
020,0 |
020,0 |
020,0 |
020,0 |
020,0 |
020,0 |
Т(051 ... 060) |
020,0 |
020,0 |
020,0 |
020,0 |
020,0 |
020,0 |
020,0 |
020,0 |
020,0 |
020,0 |
Т(061 ... 070) |
020,0 |
020,0 |
020,0 |
020,0 |
020,0 |
020,0 |
020,0 |
020,0 |
020,0 |
020,0 |
Т(071 ... 080) |
020,0 |
020,0 |
020,0 |
020,0 |
020,0 |
020,0 |
020,0 |
020,0 |
020,0 |
020,0 |
Т(081 ... 090) |
020,0 |
020,0 |
020,0 |
020,0 |
020,0 |
020,0 |
020,0 |
020,0 |
020,0 |
020,0 |
= > biblal\tx26.dat < = - файл с теплофизичес. характеристиками материалов
= > biblal\ge26.dat < = - файл с номерами и размерами всех узлов
= > biblal\wr26.dat < = - файл с массивом шагов по времени и темпер. среды
= > biblal\tn16.dat < = - файл с начальной температурой
= > biblal\niki.txt < = - файл для результатов вычислений
DIMENSION N1(300), N2(300), N3(300), N4(300), NM(300), RT(300),
*DR(300), DH(300), A1(300), A2(300), A3(300), A4(300), A5(300), T(315)
*, TH(300), TP(315), CT(300), D(30), KD(30), TP1(25), TP2(25), TP3(25),
*TP4(25), PR1(25), PR2(25), PR3(25), PR4(25), MKT(15), MKK(15),
*VTS(15, 20), TS(15, 20), VKT(15, 15), TK(15, 15), NT1(20), NT2(20)
DIMENSION SUM(99), SAM(99), OB(99), NI(50), KI(50), NL(50), KP(50)
character*75 fdi1, fdi2, fdi3, fdi4, fdi5
character*15 idi1, idi2, idi3, idi4, idi5
character*3 ioi1, ioi2, ioi3
open (5, file = 'biblal\osn1.dat', status = 'old')
read (5, '(4x, a, /, 3(4x, a, /), 4x, a)') idi1, idi2, idi3, idi4, idi5
write (*, '(5x, a, 2x, a)')' файл с материалами ->', idi1
write (*, '(5x, a, 2x, a)')' файл с геометрией ->', idi2
write (*, '(5x, a, 2x, a)')' файл с шагами времени ->', idi3
write (*, '(5x, a, 2x, a)')' файл с нач. температурой ->', idi4
write (*, '(5x, a, 2x, a)')' выходной файл ->', idi5
wase = 0.9E13
smmax1 = 0.
smmax2 = 0.
close (5)
open (5, file = idi1, status = 'old')
read (5, '(5(3x, a, /))') fdi1, fdi2, fdi3, fdi4, fdi5
write (*, '(7(3x, a, /))') fdi1, fdi2, fdi3, fdi2, fdi4, fdi2, fdi5
read (5, '(4x, i2)') KM
write (*, '(25x, a, 1x, i2/)') 'КОЛИЧЕСТВО МАТЕРИАЛОВ = ', KM
read (5, '(5x, a, 10x, a)') ioi1, ioi2
read (5, '(14x, 4(lx, f6.1), /, 1x)') TG1, TG2, TG3, TG4
write (*, '(5x, 4(8x, a, f6.1))'i) 'T1 = ', TG1, 'T2 = ', TG2, 'T3 = ', TG3,
*'T4 =', TG4
open (6, file = idi5, status = 'old')
write (6, '(7(3x, a, /))') fdi1, fdi2, fdi3, fdi2, fdi4, fdi2, fdi5
write (6, '(25x, a, lx, i2, /)') 'КОЛИЧЕСТВО МАТЕРИАЛОВ = ', KM
write (6, '(lx, 4(8x, a, f6.1), /)') 'T1 = ', TG1, 'T2 = ', TG2, 'T3 = ', TG3,
*'T4 = ', TG4
write (6, '(8x, a, 12x, a, /)') 'теплопроводность',
*'плотность * теплоемкость'
do 900 I = 1, KM
read (5, '(14x, 4(1x, f6.3), 2x, 4(1x, f7.2))') TP1(I), TP2(I), TP3(I),
*TP4(I), PR1(I), PR2(I), PR3(I), PR4(I)
write (6, '(2x, i2, 4(2x, f6.3), 4x, 4(1x, f7.2))') I, TP1(I), TP2(I),
*TP3(I), TP4(I), PR1(I), PR2(I), PR3(I), PR4(I)
900 write (*, '(2x, i2, 4(2x, f6.3), 4x, 4(1x, f7.2)) ') I, TP1(I), TP2(I),
*TP3(I), TP4(I), PR1(I), PR2(I), PR3(I), PR4(I)
close (5)
open (5, file = idi2, status = 'old')
read (5, '(7(5x, i3, /), /, /)') NG1, ND, NG2, N, NH, L, KNN
write (6, '(2x, /, 12x, A)') 'геометрические характеристики узлов'
write (6, '(4x, /, 7(3x, a, i3), /)') 'NG1 = ', NG1, 'ND = ', ND, 'NG2 = ', NG2,
*'N = ', N, 'NH = ', NH, 'L = ', L, 'KNN = ', KNN
write (*, '(7(3x, а, i3))') 'NG1 = ', NG1, 'ND = ', ND, 'NG2 = ', NG2,
*'N = ', N, 'NH = ', NH, 'L = ', L, 'KNN = ', KNN
do 901 I = 1, NG1
read (5, '(6x, 4(3x, i3), 3(3x, f8.6), 3x, i2)') N1(I), N2(I), N3(I), N4(I),
*RT(I), DR(I), DH(I), NM(I)
901 write (6, '(3x, i3, 4(3x, i3), 3(3x, f8.6), 3x, i2)') I, N1(I), N2(I), N3(I),
*N4(I), RT(I), DR(I), DH(I), NM(I)
IF (ND.EQ.N) GO TO 903
902 write (6, '(3x, i3, 4(3x, i3), 3(3x, f8.6), 3x, i2)') I, N1(I), N2(1), N3(1),
*N4(I), RT(I), DR(I), DH(I), NM(I)
903 continue
read (5, '(9x, a, /, /)') ioi3
if (ioi3.EQ.ioi2) go to 904
write (*, '(5x, a)') '-> лишние узлы в файле ge01.dat - N1, N2 ... '
go to 999
904 do 905 I 1, KNN
read (5, '(6x, 4(3x, i3))') NI(I), NL(I), KI(I), KP(I)
write (6, '(3x, i3, 4(3x, a, i3))') I, 'NI = ', NI(I), 'NL = ', NL(I), 'KI = ', *KI(I), 'KP', KP(I)
905 write (*, '(3x, i3, 4(3x, a, i3))') I, 'NI = ', NI(I), 'NL = ', NL(I), 'KI = ', *KI(I), 'KP = ', KP(I)
read (5, '(9x, a, /, /)') ioi3
if( ioi3.EQ.ioi2) go to 906
write (*, '(5x, a) ') '-> лишние записи в файле ge01.dat - NI, NL... '
go to 999
906 IF(ND.EQ.N) GO TO 908
KTG = NG2 - ND + 1
do 907 I = 1, KTG
read (5, '(6x, 2(3x, i3))') NT1(I), NT2(I)
write (6, '(3x, i2, 2(3x, a, i3))') I, 'NT1 = ', NT1(I), 'NT2 = ', NT2(I)
907 write (*, '(3x, i3, 2(3x, a, i3))') I, 'NT1 = ', NT1(I), 'NT2 = ', NT2(I)
908 close (5)
open (7, file = 'prn', status = 'old')
open (5, file = idi3, status = 'old')
read (5, '(7x, i6, /, 7x, i6, /, 7x, f6.3, /, 7x, f6.4, /, 7x, f6.3, /, /, /)')
*KVD, MM, DT1, DT2, RK
WRITE (6, '(2x, /, 2X, 2(2X, A, I4), 2(2X, A, F6.3), 2X, A, F6.4, /)') 'KVD = '
*, KVD, 'MM = ', MM, 'RK = ', RK, 'DT1 =', DT1, 'DT2 = ', DT2
WRITE (*, '(2X, 2(2X, A, I4), 2(2X, A, F6.3), 2X, A, F6.4)') 'KVD = ', KVD,
*'MM = ', MM, 'RK = ', RK, 'DT1 = ', DT1, 'DT2 = ', DT2
do 912 I = 1, KVD
READ (5, '(9X, F7.1, 3X, I2)') D(I), KD(I)
WRITE (6, '(9X, I3, 3X, F7.1, 3X, I2)') I, D(I), KD(I)
912 WRITE (*, '(9X, I3, 3X, F7.1, 3X, I2)') I, D(I), KD(I)
read (5, '(3x, a, /, /)') ioi3
if (ioi3.EQ.ioi2) go to 913
write (*, '(5x, a)') '-> лишние записи в файле wr01.dat - D, KD'
go to 999
913 KTS = L - NH + 1
DO 914 I = 1, KTS
READ (5, '(3X, I3, 3X, I3)') MKT(I), MKK(I)
WRITE (6, '(9X, I3, 3X, A, I3, 3X, A, I3)') I, 'MKT = ', MKT(I), 'MKK = ', MKK(I)
914 WRITE (*, '(9X, I3, 3X, A, I3, 3X, A, I3)') I, 'MKT = ', MKT(I), 'MKK = ', MKK(I)
read (5, '(3x, a, /, /)') ioi3
if (ioi3.EQ.ioi2) go to 915
write (*, '(5x, a)') '-> лишние записи в файле wr01.dat - MКТ, МКК'
go to 999
915 DO 920 I = 1, KTS
DO 916 K = 1, MKT(I)
read (5, '(3x, f7.1, 3x, f5.1)') VTS(I, K), TS(I, K)
WRITE (6, '(3X, A, F7.1, 3X, A, F5.1)') 'VTS = ', VTS(I, K), 'TS = ', TS(I, K)
916 WRITE (*, '(3X, A, F7.1, 3X, A, F5.1)') 'VTS = ', VTS(I, K), 'TS = ', TS(I, K)
read (5, '(3x, a, /, /)') ioi3
if (ioi3.EQ.ioi2) go to 917
write (*, '(5x, 12, a)') I, 'лишние записи в файле wr01.dat - VTS, TS'
918 WRITE (*, '(3X, A, F7.1, 2X, A, F10.4)') 'VKT = ', VKT(I, K), 'TK = ', TK(I, K)
read (5, '(3x, a, /, /}') ioi3
if (ioi3.EQ.ioi2) go to 919
write {*, '(5x, 12, a)') I, 'лишние записи в файле wr01.dat - VКТ, ТК'
919 CONTINUE
920 CONTINUE
close (5)
IXOD = 1
open (5, file = idi4, status = 'old')
do 909 I = 1, 300
909 T(I) = TS(1, 1)
read (5, '(7x, i1, /, /, /)') KN
write (*, '(10x, a, i2) ') 'KN = ', KN
M = 0
KPR = 0
777 KPR = KPR + 1
KPR1 = KPR*10
IF (KPR1.GT.N) GO TO 778
GO TO 777
778 IF (KN.EQ.0) GO TO 911
write (6, '(25x, a, /)') 'НАЧАЛЬНОЕ РАСПРЕДЕЛЕНИЕ ТЕМПЕРАТУРЫ'
do 910 I = 1, 30
K = 10*I
READ (5, '(16X, 10(1X, F5.1))') T(K-9), T(K-8), T(K-7), T(K-6),
*T(K-5), T(K-4), T(K-3), T(K-2), T(K-1), T(K)
910 write (6, '(16X, 10(1X, F5.1))') T(X-9), T(K-8), T(K-7), T(K-6),
*T(K-5), T(K-4), T(K-3), T(K-2), T(K-l), T(K)
911 continue
107 close (5)
VR = VTS(1,1)
DO 189 JD = 1, KVD
DVR = D(JD)
KVD = KD(JD)
DO 189 ID = 1, KDV
С ИНТЕРПОЛЯЦИЯ
VR = VR + DVR
DO 130 I = 1, KTS
IT = I + NH - 1
K = MKT(I)
KK = MKK(I)
IF (K.EQ.1.OR.VR.GT.VTS(I, K)) GO TO 97
IK = K - 1
DO 98 J = 1, IK
A = VTS(I, J)
B = VTS(I, J + 1)
IF (VR.LE.A.OR.B.LT.VR) GO TO 98
T(IT) = TS(I, J) + (TS(I, J + l) - TS(I, J))*(VR - A)/(B - A)
98 CONTINUE
GO TO 99
97 T(IT) = TS(I, K)
99 IF (KK.EQ.1.OR.VR.GT.VKT(I, KK)) GO TO 124
IKK = KK - 1
DO 122 JK = 1, IKK
A = VKT(I, JK)
B = VKT(I, JK + 1)
IF (VR.LE.A.OR.B.LT.VR) GO TO 122
TP(IT) = TR(I, JK) + (TK(I, JK + 1) - TK(I, JK))*(VR-A)/(B-A)*)
124 ТP(IT) TК(I, КК)
130 CONTINUE
K1 = 1
K2 = NG1
77 DO 96 I = K1, K2
K = NM(I)
TT = T(I)
IF (TT.LT.TG2) GO TO 93
IF (TG2.LE.TT.AND.TT.LE.TG3) GO TO 94
DTG = TG4 - TG3
DT = TT - TG3
A = TP3(K)
B = PR3(K)
DTP = TP4(K) - A
DPR = PR4(K) - B
GO TO 95
93 DTG = TG2 - TG1
DТ = TT - TG1
A = TP1(K)
B = PR1(K)
DTP = TP2(K) - A
DPR = PR2(K) - B
GO TO 95
94 DTG = TG3 - TG2
DT = TT - TG2
A = TP2(K)
B = PR2(K)
DTP = TR3(K) - A
DPR = PR3(K) - B
95 C = DT/DTG
TP(I) = DTP*C + A
P = DPR*C + B
96 CT(I) = P*3.1416*8.*RT(I)*DR(I)*DH(I)
С РАСЧЕТ КОЭФФИЦИЕНТОВ
DO 37 I = K1, K2
TH(I) = T(I)
A = RT(I)
B = DR(I)
C = DH(I)
P = 4*3,1416*TP(I)
P1 = P*C
P2 = P*A*B
K = N1(I)
if (K.EQ.0) go to 9
GO TO 11
9 R1 = 0.9E13
GO TO 14
11 R48 = A/(A - B)
R1 = LOG(R48)/P1
P3 = 4.*3.1416*C*TP(K)
if (K.GT.N) go to 13
12 R1 = RH + LOG(1.+DR(K)/RT(K))/P3
GO TO 14
13 R1 = R1 + 1./P3/(A - B)
14 K = N2(I)
if (K.EQ.0) go to 15
GO TO 16
Р3 = 4.*3.1416*С*ТР(К)
if (K.GT.N) go to 19
18 R2 = R2 + LOG(RT(K)/(RT(K) - DR(K)))/P3
GO TO 22
19 R2 = R2 + 1./P3/(A + B)
22 K = N3(I)
if (K.EQ.0) go to 24
GO TO 25
24 R3 = 0.9E13
GO TO 29
25 R3 = C/P2
Р3 = 4.*3.1416*ТР(К)
if (K.GT.N) go to 28
27 R3 = R3 + DH(K)/P3/RT(K)/DR(K)
GO TO 29
28 R3 = R3 + 1./P3/A/B
29 K = N4(I)
if (K.EQ.0) go to 31
GO TO 32
31 R4 = 0.9E13
GO TO 36
32 R4 = C/P2
P3 = 4.*3.1416*TP(K)
if (K.GT.K) go to 34
33 R4 = R4 + DH(K)/P3/RT(K)/DR(K)
GO TO 36
34 R4 = R4 + 1./P3/A/B
36 S = 1./R1 + 1./R2 + 1./R3 + 1./R4 + 3600.*CT(I)/DVR
A5(I) = 3600.*CT(I)/DVR/S
A1(I) = 1./R1/S
A2(I) = 1./R2/S
A3(I) = 1./R3/S
A4(I) = 1./R4/S
37 CONTINUE
C РАСЧЕТ ТЕМПЕРАТУР
IF (K1.NE.1) K1 = NG2 + 1
M = 1
39 NH = 0
DO 66 I = K1, K2
K = N1(I)
if (K.EQ.0) go to 41
GO TO 42
41 T1 = 0
GO TO 44
42 T1 = T(K)
44 K = N2(I)
if (K.EQ.0) go to 45
GO TO 46
45 T2 = 0.
GO TO 48
46 T2 = T(K)
48 K = N3(I)
if (K.EQ.0) go to 49
GO TO 51
49 T3 = 0.
GO TO 52*)
85 continue
DO 701 I = NH, L
701 WRITE (6, '(10Х, А, I3, A, F6.1)') 'TEMPERATURA', I, ' SREDY =', T(I)
do 700 I = 1, KNN
write (6, '(10x, A, a, f6.1, A, I3, A, I3)') 'среднемассовая',
*'температура = ', SAM(I), ' для NI = ', NI(I), ' и KI = ', KI(I)
700 continue
С
IF (knn.ne.4) GO TO 189
IF (knn.ne.4) GO TO 189
SAM13 = SAM(1) - SAM(3)
SAM42 = SAM(4) - SAM(2)
IF (SMMAX1.GT.abs(SAM13)) GO TO 888
SMMAX1 = SAM13
vra = vr
888 CONTINUE
IF (SMMAX2.GT.abs(SAM42)) GO TO 889
SMMAX2 = SAM42
vri = vr
889 CONTINUE
write (6, '(10X, A, a, f6.1)') 'перепад температур',
*'t1 - t3 = ', SAM13
write (6, '(10X, A, a, f6.1)') 'перепад температур',
*'t4 - t2 = ', SAM42
IF (ND.EQ.N) GO TO 189
189 CONTINUE
IF (ND.EQ.N) GO TO 193
101 continue
111 continue
193 CONTINUE
999 CONTINUE
IF (knn.ne.4) GO TO 1899
С
write (6, '(10X, A, a, f6.1)') 'МАКСИМАЛЬНЫЙ ПЕРЕПАД',
*'t1 - t3 = ', SMMAX1
write (6, '(10x, A, a, f6.1)') 'В МОМЕНТ ВРЕМЕНИ',
*' = ', vra
write (6, '(10X, A, a, f6.1)') 'МАКСИМАЛЬНЫЙ ПЕРЕПАД',
*'t4 - t2 = ', SMMAX2
write (6, '(10X, A, a, f6.1)') 'В МОМЕНТ ВРЕМЕНИ',
*' = ', vri
1899 STOP
END
С SIM(KK) = (RT(KI(KK)) + DR(KI(KK)))**2. - (RT(NI(KK)) - DR(NI(KK)))**2
SAM(KK) = SUM(KK)/OB(KK)
310 CONTINUE
IF (ND.EQ.N) GO TO 81
IF (K1.NE.1) GO TO 81
DO 80 K = 1, KTG
IT = ND+ K - 1
I = NT1(K)
J = NT2(K)
80 T(IT) = (T(I)+T(J))/2
K1 = ND
K2 = N
GO TO 77
КЛАПАНА ЗАПОРНОГО DN 100 16.07.96 г. (крышка, корпус из ст. 20) исп. Сергевнина Г.И. |
КОЛИЧЕСТВО МАТЕРИАЛОВ = 5
Т1 = .0 Т2 = 120.0 Т3 = 240.0 Т4 = 360.0
плотность * теплоемкость |
||||||||
1 |
11,340 |
13,150 |
14,940 |
16,640 |
862,00 |
942,00 |
984,00 |
1010,00 |
2 |
10,600 |
12,800 |
14,900 |
16,600 |
1000,00 |
1005,00 |
1020,00 |
1040,00 |
3 |
43,970 |
41,490 |
39,010 |
36,510 |
742,00 |
885,00 |
990,00 |
1080,00 |
4 |
51,900 |
50,300 |
47,100 |
43,700 |
894,00 |
931,00 |
998,00 |
1076,00 |
5 |
64,700 |
76,400 |
86,300 |
93,000 |
600,00 |
600,00 |
600,00 |
600,00 |
геометрические характеристики узлов
NG1 = 85 ND = 85 NG2 = 85 N = 85 NH = 301 L = 302 KNN = 4
301 |
2 |
6 |
0 |
,060500 |
,003500 |
,075000 |
3 |
|||||||
2 |
1 |
3 |
7 |
0 |
,069000 |
,005000 |
,075000 |
3 |
||||||
3 |
2 |
4 |
8 |
302 |
,078500 |
,004500 |
,075000 |
3 |
||||||
4 |
3 |
5 |
9 |
302 |
,085000 |
,002000 |
,075000 |
3 |
||||||
5 |
4 |
302 |
10 |
302 |
,096000 |
,009000 |
,075000 |
3 |
||||||
6 |
301 |
7 |
22 |
1 |
,060500 |
,003500 |
,075000 |
3 |
||||||
7 |
6 |
8 |
23 |
2 |
,069000 |
,005000 |
,075000 |
3 |
||||||
8 |
7 |
9 |
24 |
3 |
,078500 |
,004500 |
,075000 |
3 |
||||||
9 |
8 |
10 |
25 |
4 |
,085000 |
,002000 |
,075000 |
3 |
||||||
10 |
9 |
302 |
26 |
5 |
,096000 |
,009000 |
,075000 |
3 |
||||||
11 |
0 |
12 |
14 |
0 |
,011000 |
,011000 |
,015000 |
1 |
||||||
12 |
11 |
13 |
15 |
301 |
,025500 |
,003500 |
,015000 |
5 |
||||||
13 |
12 |
301 |
16 |
301 |
,032500 |
,003500 |
,015000 |
5 |
||||||
14 |
0 |
15 |
17 |
11 |
,011000 |
,011000 |
,015000 |
1 |
||||||
15 |
14 |
16 |
18 |
12 |
,025500 |
,003500 |
,015000 |
5 |
||||||
16 |
15 |
301 |
19 |
13 |
,032500 |
,003500 |
,015000 |
5 |
||||||
17 |
0 |
18 |
29 |
14 |
,011000 |
,011000 |
,002500 |
1 |
||||||
18 |
17 |
19 |
30 |
15 |
,025500 |
,003500 |
,002500 |
5 |
||||||
19 |
18 |
20 |
31 |
16 |
,032500 |
,003500 |
,002500 |
5 |
||||||
20 |
19 |
21 |
32 |
301 |
,041000 |
,005000 |
,002500 |
1 |
||||||
21 |
20 |
22 |
33 |
301 |
,051500 |
,005500 |
,002500 |
1 |
||||||
22 |
21 |
23 |
34 |
6 |
,060500 |
,003500 |
,002500 |
3 |
||||||
23 |
22 |
24 |
302 |
7 |
,069000 |
,005000 |
,002500 |
3 |
||||||
24 |
23 |
25 |
302 |
8 |
,078500 |
,004500 |
,002500 |
3 |
||||||
25 |
24 |
26 |
27 |
9 |
,085000 |
,002000 |
,002500 |
3 |
||||||
26 |
25 |
302 |
302 |
10 |
,096000 |
,009000 |
,002500 |
3 |
||||||
27 |
302 |
302 |
28 |
25 |
,085000 |
,002000 |
,008500 |
2 |
||||||
28 |
302 |
302 |
49 |
27 |
,085000 |
,002000 |
,008500 |
2 |
||||||
29 |
0 |
30 |
35 |
17 |
,011000 |
,011000 |
,008500 |
1 |
||||||
30 |
29 |
31 |
36 |
18 |
,025500 |
,003500 |
,008500 |
5 |
||||||
31 |
30 |
32 |
37 |
19 |
,032500 |
,003500 |
,008500 |
5 |
||||||
32 |
31 |
33 |
38 |
20 |
,041000 |
,005000 |
,008500 |
3 |
||||||
33 |
32 |
34 |
39 |
21 |
,051500 |
,005500 |
,008500 |
3 |
||||||
34 |
33 |
302 |
40 |
22 |
,060500 |
,003500 |
,008500 |
3 |
||||||
35 |
0 |
36 |
41 |
29 |
,011000 |
,011000 |
,008500 |
1 |
||||||
36 |
35 |
37 |
42 |
30 |
,025500 |
,003500 |
,008500 |
5 |
||||||
37 |
36 |
38 |
43 |
31 |
,032500 |
,003500 |
,008500 |
5 |
||||||
38 |
37 |
39 |
44 |
32 |
,041000 |
,005000 |
,008500 |
5 |
||||||
39 |
38 |
40 |
45 |
33 |
,051500 |
,005500 |
,008500 |
5 |
||||||
40 |
39 |
302 |
46 |
34 |
,060500 |
,003500 |
,008500 |
5 |
||||||
41 |
0 |
42 |
51 |
35 |
,011000 |
,011000 |
,010000 |
1 |
||||||
42 |
41 |
43 |
52 |
36 |
,025500 |
,003500 |
,010000 |
3 |
||||||
43 |
42 |
44 |
53 |
37 |
,032500 |
,003500 |
,010000 |
3 |
||||||
44 |
43 |
45 |
54 |
38 |
,041000 |
,005000 |
,010000 |
3 |
||||||
45 |
44 |
46 |
55 |
39 |
,051500 |
,005500 |
,010000 |
3 |
||||||
46 |
45 |
47 |
56 |
40 |
,060500 |
,003500 |
,010000 |
3 |
||||||
47 |
46 |
48 |
57 |
302 |
,069000 |
,005000 |
,010000 |
3 |
||||||
48 |
47 |
49 |
58 |
302 |
,078500 |
,004500 |
,010000 |
3 |
||||||
49 |
48 |
50 |
59 |
28 |
,085000 |
,002000 |
,010000 |
2 |
||||||
50 |
49 |
302 |
60 |
302 |
,096000 |
,009000 |
,010000 |
3 |
||||||
51 |
0 |
52 |
61 |
41 |
,011000 |
,011000 |
,004500 |
1 |
||||||
52 |
51 |
53 |
62 |
42 |
,025500 |
,003500 |
,004500 |
4 |
||||||
53 |
52 |
54 |
63 |
43 |
,032500 |
,003500 |
,004500 |
4 |
||||||
54 |
53 |
55 |
64 |
44 |
,041000 |
,005000 |
,004500 |
4 |
||||||
55 |
54 |
56 |
65 |
45 |
,051500 |
,005500 |
,004500 |
3 |
||||||
56 |
55 |
57 |
66 |
46 |
,060500 |
,003500 |
,004500 |
3 |
||||||
57 |
56 |
58 |
302 |
47 |
,069000 |
,005000 |
,004500 |
3 |
||||||
58 |
57 |
59 |
302 |
48 |
,078500 |
,004500 |
,004500 |
3 |
||||||
59 |
58 |
60 |
302 |
49 |
,085000 |
,002000 |
,004500 |
2 |
||||||
60 |
59 |
302 |
302 |
50 |
,096000 |
,009000 |
,004500 |
3 |
||||||
61 |
0 |
62 |
67 |
51 |
,011000 |
,011000 |
,007000 |
1 |
||||||
62 |
61 |
63 |
68 |
52 |
,025500 |
,003500 |
,007000 |
4 |
||||||
63 |
62 |
64 |
69 |
53 |
,032500 |
,003500 |
,007000 |
4 |
||||||
64 |
63 |
65 |
70 |
54 |
,041000 |
,005000 |
,007000 |
4 |
||||||
65 |
64 |
66 |
71 |
55 |
,051500 |
,005500 |
,007000 |
3 |
||||||
66 |
65 |
302 |
72 |
56 |
,060500 |
,003500 |
,007000 |
3 |
||||||
67 |
0 |
68 |
73 |
61 |
,011000 |
,011000 |
,007000 |
1 |
||||||
68 |
67 |
69 |
74 |
62 |
,025500 |
,003500 |
,007000 |
4 |
||||||
69 |
68 |
70 |
75 |
63 |
,032500 |
,003500 |
,007000 |
4 |
||||||
70 |
69 |
71 |
76 |
64 |
,041000 |
,005000 |
,007000 |
4 |
||||||
71 |
70 |
72 |
77 |
65 |
,051500 |
,005500 |
,007000 |
3 |
||||||
72 |
71 |
302 |
78 |
66 |
,060500 |
,003500 |
,007000 |
3 |
||||||
73 |
0 |
74 |
79 |
67 |
,011000 |
,011000 |
,007500 |
1 |
||||||
74 |
73 |
75 |
80 |
68 |
,025500 |
,003500 |
,007500 |
4 |
||||||
75 |
74 |
76 |
81 |
69 |
,032500 |
,003500 |
,007500 |
4 |
||||||
76 |
75 |
77 |
82 |
70 |
,041000 |
,005000 |
,007500 |
4 |
||||||
77 |
76 |
78 |
83 |
71 |
,051500 |
,005500 |
,007500 |
3 |
||||||
78 |
77 |
302 |
84 |
72 |
,060500 |
,003500 |
,007500 |
3 |
||||||
79 |
0 |
80 |
0 |
73 |
,011000 |
,011000 |
,010000 |
1 |
||||||
80 |
79 |
81 |
0 |
74 |
,025500 |
,003500 |
,010000 |
4 |
||||||
81 |
80 |
82 |
0 |
75 |
,032500 |
,003500 |
,010000 |
4 |
||||||
82 |
81 |
83 |
0 |
76 |
,041000 |
,005000 |
,010000 |
4 |
||||||
83 |
82 |
84 |
0 |
77 |
,051500 |
,005500 |
,010000 |
3 |
||||||
84 |
83 |
85 |
0 |
78 |
,060500 |
,003500 |
,010000 |
3 |
||||||
85 |
84 |
302 |
0 |
302 |
,074500 |
,010500 |
,010000 |
3 |
||||||
1 |
N1 = 1 |
NL = 22 |
KI = 26 |
KP = 5 |
|
|||||||||
2 |
N1 = 25 |
NL = 59 |
KI = 59 |
КР = 25 |
|
|||||||||
3 |
N1 = 46 |
NL = 56 |
KI = 60 |
KР = 50 |
|
|||||||||
4 |
N1 = 32 |
NL = 38 |
KI = 40 |
КР = 34 |
|
|||||||||
KVD = 22 ММ = 100 RK = 1,800 DT1 = ,100 DT2 = ,0010
,5 |
10 |
|
|
2 |
5,0 |
5 |
|
3 |
10,0 |
3 |
|
4 |
20,0 |
3 |
|
5 |
30,0 |
6 |
|
6 |
60,0 |
5 |
|
7 |
120,0 |
5 |
|
8 |
300,0 |
8 |
|
9 |
600,0 |
6 |
|
10 |
1800,0 |
6 |
|
11 |
3600,0 |
5 |
|
12 |
,5 |
10 |
|
13 |
5,0 |
5 |
|
14 |
10,0 |
3 |
|
15 |
20,0 |
3 |
|
16 |
30,0 |
6 |
|
17 |
60,0 |
5 |
|
18 |
120,0 |
5 |
|
19 |
300,0 |
8 |
|
20 |
600,0 |
6 |
|
21 |
1800,0 |
6 |
|
22 |
3600,0 |
5 |
|
1 |
МКТ = 5 |
МКК = 3 |
|
2 |
МКТ = 1 |
МКК = 1 |
|
VTS = ,0 |
ТS = 20,0 |
||
VTS = 10,0 |
ТS = 190,0 |
||
VTS = 36000,0 |
ТS = 190,0 |
||
VTS = 36010,0 |
ТS = 20,0 |
||
VTS = 72000,0 |
ТS = 20,0 |
||
VTK = ,0 |
ТК = 1900,0000 |
||
VTK = 36000,5 |
ТК = 1100,0000 |
||
VTK = 36010,0 |
ТК = 1100,0000 |
||
VTS = ,0 |
ТS = 20,0 |
||
VTS = ,0 |
ТК = 7,0000 |
ВРЕМЯ = .5 секунд ПОЛЕ ТЕМПЕРАТУР кол. итераций = 8
ТЕМРЕRATURA 301 SREDY = 28.5
ТЕМРЕRATURA 302 SREDY = 20.0
среднемассовая температура = 20,0 для NI = 1 и КI = 26
среднемассовая температура = 20,0 для NI = 25 и КI = 59
среднемассовая температура = 20,0 для NI = 46 и КI = 60
среднемассовая температура = 20,0 для NI = 32 и КI = 40
перепад температур t1 - t3 = ,0
перепад температур t4 - t2 = ,0
ВРЕМЯ 1,0 секунд ПОЛЕ ТЕМПЕРАТУР кол. итераций = 10
ТЕМРЕRATURA 301 SREDY = 37,0
ТЕМРЕRATURA 302 SREDY = 20,0
среднемассовая температура = 20,1 для NI = 1 и КI = 26
среднемассовая температура = 20,0 для NI = 25 и КI = 59
среднемассовая температура = 20,0 для NI = 46 и КI = 60
среднемассовая температура = 20,0 для NI = 32 и КI = 40
перепад температур t1 - t3 = ,1
перепад температур t4 - t2 = ,0
ВРЕМЯ = 1,5 секунд ПОЛЕ ТЕМПЕРАТУР кол. итераций = 13
ТЕМРЕRATURA 301 SREDY = 45,5
ТЕМРЕRATURA 302 SREDY = 20,0
среднемассовая температура = 20,2 для NI = 1 и КI = 26
среднемассовая температура = 20,0 для NI = 25 и КI = 59
среднемассовая температура = 20,0 для NI = 46 и КI = 60
среднемассовая температура = 20,0 для NI = 32 и КI = 40
перепад температур t1 - t3 = ,2
перепад температур t4 - t2 = ,0
ВРЕМЯ = 2,0 секунд ПОЛЕ ТЕМПЕРАТУР кол. итераций = 13
ТЕМРЕRATURA 301 SREDY = 54,0
ТЕМРЕRATURA 302 SREDY = 20,0
среднемассовая температура = 20,3 для NI = 1 и КI = 26
среднемассовая температура = 20,0 для NI = 25 и КI = 59
среднемассовая температура = 20,0 для NI = 46 и КI = 60
среднемассовая температура = 20,1 для NI = 32 и КI = 40
перепад температур t1 - t3 = ,3
перепад температур t4 - t2 = ,1
ВРЕМЯ = 2,5 секунд ПОЛЕ ТЕМПЕРАТУР кол. итераций = 14
ТЕМРЕRATURA 301 SREDY = 62,5
ТЕМРЕRATURA 302 SREDY = 20,0
среднемассовая температура = 20,5 для NI = 1 и КI = 26
среднемассовая температура = 20,0 для NI = 25 и КI = 59
среднемассовая температура = 20,0 для NI = 46 и КI = 60
среднемассовая температура = 20,1 для NI = 32 и КI = 40
перепад температур t1 - t3 = .5
перепад температур t4 - t2 = .1
ВРЕМЯ = 3,0 секунд ПОЛЕ ТЕМПЕРАТУР кол. итераций = 14
ТЕМРЕRATURA 301 SREDY = 71,0
ТЕМРЕRATURA 302 SREDY = 20,0
среднемассовая температура = 20,7 для NI = 1 и КI = 26
среднемассовая температура = 20,0 для NI = 25 и КI = 59
среднемассовая температура = 20,0 для NI = 46 и КI = 60
среднемассовая температура = 20,1 для NI = 32 и КI = 40
перепад температур t1 - t3 = ,7
перепад температур t4 - t2 = ,1
ВРЕМЯ = 3,5 секунд ПОЛЕ ТЕМПЕРАТУР кол. итераций = 14
ТЕМРЕRATURA 301 SREDY = 79,5
ТЕМРЕRATURA 302 SREDY = 20,0
среднемассовая температура = 20,9 для NI = 1 и КI = 26
среднемассовая температура = 20,0 для NI = 25 и КI = 59
среднемассовая температура = 20,0 для NI = 46 и КI = 60
среднемассовая температура = 20,2 для NI = 32 и КI = 40
перепад температур t1 - t3 = ,9
перепад температур t4 - t2 = ,2
ВРЕМЯ = 4,0 секунд ПОЛЕ ТЕМПЕРАТУР кол. итераций = 15
ТЕМРЕRATURA 301 SREDY = 88,0
ТЕМРЕRATURA 302 SREDY = 20,0
среднемассовая температура = 21,2 для NI = 1 и КI = 26
среднемассовая температура = 20,0 для NI = 25 и КI = 59
среднемассовая температура = 20,0 для NI = 46 и КI = 60
среднемассовая температура = 20,3 для NI = 32 и КI = 40
перепад температур t1 - t3 = 1,2
перепад температур t4 - t2 = ,3
ВРЕМЯ = 4,5 секунд ПОЛЕ ТЕМПЕРАТУР кол. итераций = 15
ТЕМРЕRATURA 301 SREDY = 96,5
ТЕМРЕRATURA 302 SREDY = 20,0
среднемассовая температура = 21,5 для NI = 1 и КI = 26
среднемассовая температура = 20,0 для NI = 25 и КI = 59
среднемассовая температура = 20,0 для NI = 46 и КI = 60
среднемассовая температура = 20,4 для NI = 32 и КI = 40
перепад температур t1 - t3 = 1,5
перепад температур t4 - t2 = ,4
ВРЕМЯ = 60,0 секунд ПОЛЕ ТЕМПЕРАТУР кол. итераций = 23
ТЕМРЕRATURA 301 SREDY = 190,0
ТЕМРЕRATURA 302 SREDY = 20,0
среднемассовая температура = 67,5 для NI = 1 и КI = 26
среднемассовая температура = 25,7 для NI = 25 и КI = 59
среднемассовая температура = 23,2 для NI = 46 и КI = 60
среднемассовая температура = 56,6 для NI = 32 и КI = 40
перепад температур t1 - t3 = 44,3
перепад температур t4 - t2 = 30,9
ВРЕМЯ = 80,0 секунд ПОЛЕ ТЕМПЕРАТУР кол. итераций = 24
ТЕМРЕRATURA 301 SREDY = 190,0
ТЕМРЕRATURA 302 SREDY = 20,0
среднемассовая температура = 79,2 для NI = 1 и КI = 26
среднемассовая температура = 29,6 для NI = 25 и КI = 59
среднемассовая температура = 25,7 для NI = 46 и КI = 60
среднемассовая температура = 64,3 для NI = 32 и КI = 40
перепад температур t1 - t3 = 53,6
перепад температур t4 - t2 = 34,7
ВРЕМЯ = 210,0 секунд ПОЛЕ ТЕМПЕРАТУР кол. итераций = 22
ТЕМРЕRATURA 301 SREDY = 190,0
ТЕМРЕRATURA 302 SREDY = 20,0
среднемассовая температура = 129,0 для NI = 1 и КI = 26
среднемассовая температура = 57,2 для NI = 25 и КI = 59
среднемассовая температура = 44,3 для NI = 46 и КI = 60
среднемассовая температура = 92,5 для NI = 32 и КI = 40
перепад температур t1 - t3 = 84,7
перепад температур t4 - t2 = 35,2
ВРЕМЯ = 240,0 секунд ПОЛЕ ТЕМПЕРАТУР кол. итераций = 22
ТЕМРЕRATURA 301 SREDY = 190,0
ТЕМРЕRATURA 302 SREDY = 20,0
среднемассовая температура = 136,3 для NI = 1 и КI = 26
среднемассовая температура = 63,1 для NI = 25 и КI = 59
среднемассовая температура = 48,6 для NI = 46 и КI = 60
среднемассовая температура = 96,6 для NI = 32 и КI = 40
перепад температур t1 - t3 = 87,7
перепад температур t4 - t2 = 33,6
ВРЕМЯ = 270,0 секунд ПОЛЕ ТЕМПЕРАТУР кол. итераций = 22
ТЕМРЕRATURA 301 SREDY = 190,0
ТЕМРЕRATURA 302 SREDY = 20,0
среднемассовая температура = 142,6 для NI = 1 и КI = 26
среднемассовая температура = 68,6 для NI = 25 и КI = 59
среднемассовая температура = 52,8 для NI = 46 и КI = 60
среднемассовая температура = 100,4 для NI = 32 и КI = 40
перепад температур t1 - t3 = 89,9
перепад температур t4 - t2 = 31,8
ВРЕМЯ = 300,0 секунд ПОЛЕ ТЕМПЕРАТУР кол. итераций = 22
ТЕМРЕRATURA 301 SREDY = 190,0
ТЕМРЕRATURA 302 SREDY = 20,0
среднемассовая температура = 148,1 для NI = 1 и КI = 26
среднемассовая температура = 73,7 для NI = 25 и КI = 59
среднемассовая температура = 56,8 для NI = 46 и КI = 60
среднемассовая температура = 103,8 для NI = 32 и КI = 40
перепад температур t1 - t3 = 91,2
перепад температур t4 - t2 = 30,1
ВРЕМЯ = 360,0 секунд ПОЛЕ ТЕМПЕРАТУР кол. итераций = 22
ТЕМРЕRATURA 301 SREDY = 190,0
ТЕМРЕRATURA 302 SREDY = 20,0
среднемассовая температура = 156,5 для NI = 1 и КI = 26
среднемассовая температура = 82,6 для NI = 25 и КI = 59
среднемассовая температура = 64,4 для NI = 46 и КI = 60
среднемассовая температура = 109,8 для NI = 32 и КI = 40
перепад температур t1 - t3 = 92,1
перепад температур t4 - t2 = 27,1
ВРЕМЯ = 600,0 секунд ПОЛЕ ТЕМПЕРАТУР кол. итераций = 21
ТЕМРЕRATURA 301 SREDY = 190,0
ТЕМРЕRATURA 302 SREDY = 20,0
среднемассовая температура = 174,9 для NI = 1 и КI = 26
среднемассовая температура = 108,5 для NI = 25 и КI = 59
среднемассовая температура = 89,6 для NI = 46 и КI = 60
среднемассовая температура = 127,5 для NI = 32 и КI = 40
перепад температур t1 - t3 = 85,2
перепад температур t4 - t2 = 19,0
ВРЕМЯ = 1200,0 секунд ПОЛЕ ТЕМПЕРАТУР кол. итераций = 19
ТЕМРЕRATURA 301 SPEDY = 190,0
ТЕМРЕRATURA 302 SPEDY = 20,0
среднемассовая температура = 184,7 для NI = 1 и КI = 26
среднемассовая температура = 138,6 для NI = 25 и КI = 59
среднемассовая температура = 126,6 для NI = 46 и КI = 60
среднемассовая температура = 151,2 для NI = 32 и КI = 40
перепад температур t1 - t3 = 58,1
перепад температур t4 - t2 = 12,6
ВРЕМЯ = 1500,0 секунд ПОЛЕ ТЕМПЕРАТУР кол. итераций = 19
ТЕМРЕRATURA 301 SREDY = 190,0
ТЕМРЕRATURA 302 SREDY = 20,0
среднемассовая температура = 185,6 для NI = 1 и КI = 26
среднемассовая температура = 146,2 для NI = 25 и КI = 59
среднемассовая температура = 137,1 для NI = 46 и КI = 60
среднемассовая температура = 158,1 для NI = 32 и КI = 40
перепад температур t1 - t3 = 48,5
перепад температур t4 - t2 = 11,8
ВРЕМЯ = 36025,0 секунд ПОЛЕ ТЕМПЕРАТУР кол. итераций = 19
ТЕМРЕRATURA 301 SREDY = 20,0
ТЕМРЕRATURA 302 SREDY = 20,0
среднемассовая температура = 173,3 для NI = 1 и КI = 26
среднемассовая температура = 167,0 для NI = 25 и КI = 59
среднемассовая температура = 167,4 для NI = 46 и КI = 60
среднемассовая температура = 165,7 для NI = 32 и КI = 40
перепад температур t1 - t3 = 5,9
перепад температур t4 - t2 = -1,4
ВРЕМЯ = 36030,0 секунд ПОЛЕ ТЕМПЕРАТУР кол. итераций = 17
ТЕМРЕRATURA 301 SREDY = 20,0
ТЕМРЕRATURA 302 SREDY = 20,0
среднемассовая температура = 170,5 для NI = 1 и КI = 26
среднемассовая температура = 166,8 для NI = 25 и КI = 59
среднемассовая температура = 167,2 для NI = 46 и КI = 60
среднемассовая температура = 163,0 для NI = 32 и КI = 40
перепад температур t1 - t3 = 3,3
перепад температур t4 - t2 = -3,8
ВРЕМЯ = 36050,0 секунд ПОЛЕ ТЕМПЕРАТУР кол. итераций = 20
ТЕМРЕRATURA 301 SREDY = 20,0
ТЕМРЕRATURA 302 SREDY = 20,0
среднемассовая температура = 160,6 для NI = 1 и КI = 26
среднемассовая температура = 165,2 для NI = 25 и КI = 59
среднемассовая температура = 166,3 для NI = 46 и КI = 60
среднемассовая температура = 153,5 для NI = 32 и КI = 40
перепад температур t1 - t3 = -5,7
перепад температур t4 - t2 = -11,7
ВРЕМЯ = 36060,0 секунд ПОЛЕ ТЕМПЕРАТУР кол. итераций = 20
ТЕМРЕRATURA 301 SREDY = 20,0
ТЕМРЕRATURA 302 SREDY = 20,0
среднемассовая температура = 155,9 для NI = 1 и КI = 26
среднемассовая температура = 164,1 для NI = 25 и КI = 59
среднемассовая температура = 165,6 для NI = 46 и КI = 60
среднемассовая температура = 149,5 для NI = 32 и КI = 40
перепад температур t1 - t3 = -9,8
перепад температур t4 - t2 = -14,6
ВРЕМЯ = 36080,0 секунд ПОЛЕ ТЕМПЕРАТУР кол. итераций = 21
ТЕМРЕRATURA 301 SREDY = 20,0
ТЕМРЕRATURA 302 SREDY = 20,0
среднемассовая температура = 147,1 для NI = 1 и КI = 26
среднемассовая температура = 161,6 для NI = 25 и КI = 59
среднемассовая температура = 164,1 для NI = 46 и КI = 60
среднемассовая температура = 143,1 для NI = 32 и КI = 40
перепад температур t1 - t3 = -16,9
перепад температур t4 - t2 = -18,5
ВРЕМЯ = 36100,0 секунд ПОЛЕ ТЕМПЕРАТУР кол. итераций = 21
ТЕМРЕRATURA 301 SREDY = 20,0
ТЕМРЕRATURA 302 SREDY = 20,0
среднемассовая температура = 138,9 для NI = 1 и КI = 26
среднемассовая температура = 158,6 для NI = 25 и КI = 59
среднемассовая температура = 162,3 для NI = 46 и КI = 60
среднемассовая температура = 137,6 для NI = 32 и КI = 40
перепад температур t1 - t3 = -23,3
перепад температур t4 - t2 = -21,0
ВРЕМЯ = 36120,0 секунд ПОЛЕ ТЕМПЕРАТУР кол. итераций = 21
ТЕМРЕRATURA 301 SREDY = 20,0
ТЕМРЕRATURA 302 SREDY = 20,0
среднемассовая температура = 131,2 для NI = 1 и КI = 26
среднемассовая температура = 155,4 для NI = 25 и КI = 59
среднемассовая температура = 160,3 для NI = 46 и КI = 60
среднемассовая температура = 132,9 для NI = 32 и КI = 40
перепад температур t1 - t3 = -29,1
перепад температур t4 - t2 = -22,5
ВРЕМЯ = 36180,0 секунд ПОЛЕ ТЕМПЕРАТУР кол. итераций = 23
ТЕМРЕRATURA 301 SREDY = 20,0
ТЕМРЕRATURA 302 SREDY = 20,0
среднемассовая температура = 110,9 для NI = 1 и КI = 26
среднемассовая температура = 145,3 для NI = 25 и КI = 59
среднемассовая температура = 153,8 для NI = 46 и КI = 60
среднемассовая температура = 122,0 для NI = 32 и КI = 40
перепад температур t1 - t3 = -42,9
перепад температур t4 - t2 = -23,4
ВРЕМЯ = 57600,0 секунд ПОЛЕ ТЕМПЕРАТУР кол. итераций = 1
ТЕМРЕRATURA 301 SREDY = 20.0
ТЕМРЕRATURA 302 SREDY = 20.0
среднемассовая температура = 20,0 для NI = 1 и КI = 26
среднемассовая температура = 20,8 для NI = 25 и КI = 59
среднемассовая температура = 21.1 для NI = 46 и КI = 60
среднемассовая температура = 20,9 для NI = 32 и КI = 40
перепад температур t1 - t3 = -1,1
перепад температур t4 - t2 = ,1
21,4 |
21,4 |
21,4 |
21,3 |
21,3 |
21,2 |
21,5 |
21,4 |
21,4 |
21,3 |
|
Т 7. |
21,3 |
21,3 |
21,5 |
21,4 |
21,4 |
21,4 |
21,3 |
21,3 |
21,5 |
21,4 |
Т 8. |
21,4 |
21,4 |
21,3 |
21,3 |
21,2 |
20,0 |
20,0 |
20,0 |
20,0 |
20,0 |
ТЕМРЕRATURA 301 SREDY = 20,0
ТЕМРЕRATURA 302 SREDY = 20,0
среднемассовая температура = 20,0 для NI = 1 и КI = 26
среднемассовая температура = 20,7 для NI = 25 и КI = 59
среднемассовая температура = 21,1 для NI = 46 и КI = 60
среднемассовая температура = 20,9 для NI = 32 и КI = 40
перепад температур t1 - t3 = -1,0
перепад температур t4 - t2 = ,1
ВРЕМЯ = 72000,0 секунд ПОЛЕ ТЕМПЕРАТУР кол. итераций = 1
ТЕМРЕRATURA 301 SREDY = 20,0
ТЕМРЕRATURA 302 SREDY = 20,0
среднемассовая температура = 20,0 для NI = 1 и КI = 26
среднемассовая температура = 20,6 для NI = 25 и КI = 59
среднемассовая температура = 20,8 для NI = 46 и КI = 60
среднемассовая температура = 20,7 для NI = 32 и КI = 40
перепад температур t1 - t3 = -,8
перепад температур t4 - t2 = ,1
МАКСИМАЛЬНЫЙ ПЕРЕПАД t1 - t3 = 92,1
МАКСИМАЛЬНЫЙ ПЕРЕПАД t4 - t2 = 37,7
Тепловой расчет клапана запорного DN 100 (в качестве примера принят клапан сильфонный НГ 26526-100), выполнен методом конечных разностей с использованием СТП 07.81-633-96 «Методика теплового расчета трубопроводной арматуры для АЭС», а также СТП 07.81-633-92 «Методика проведения теплового расчета высокопараметрической арматуры клапанного типа». Расчеты проведены на персональной ЭВМ.
Тепловые расчёты проведены для клапана, корпус и крышка которого выполнены из ст. 20.
1. ЗАДАЧА РАСЧЕТА
Тепловой расчёт клапана запорного DN 100 выполнен с целью получения среднемассовых температур деталей фланцевого соединения.
2. УСЛОВИЯ РАСЧЕТА
2.1. Контактные термические сопротивления между соприкасающимися деталями не учитывались.
2.2. При тепловом расчете клапана рассматривалась осесимметричная двухмерная нестационарная задача в цилиндрической системе координат.
2.3. Теплофизические характеристики приняты постоянными для каждого элементарного участка.
2.4. Клапан без теплоизоляции.
2.5. Для расчётов выбраны из ТУ 26-07-1407-86 наиболее жесткие режимы работы клапана.
3. ИСХОДНЫЕ ДАННЫЕ
3.1. Основные технические данные и характеристики клапана приведены в таблице 1.
Температура рабочей среды, °С |
Рабочее давление, МПа |
Вид рабочей среды |
Материал основных деталей |
Основание (документ) |
|
НГ26526-100 |
до 200 |
4,0 |
вода |
корпус-стойка - ст. 20 |
ТУ 26-07-1407-86 |
3.2. Температура окружающей среды - 20 °С.
3.3. Режимы для расчёта приведены в таблице 2.
Основание (Документ) |
||||
Начальная температура, °С |
Конечная температура, °С |
Скорость изменения температуры, °С/c |
Скорость движения среды, м/с |
|
20 |
190 |
скачок |
5 |
ТУ 26-07-1407-86 |
4. РАСЧЕТ
4.1. Исходные данные для расчёта коэффициента теплоотдачи а1 от рабочей среды к внутренней поверхности проточной части клапана приведены в таблице 3, расчет коэффициента а1 представлен в таблице 4.
4.2. Коэффициент теплоотдачи от рабочей среды к внутренней поверхности горловины а2 принят равным 0,1 а1 исходя из конструктивных особенностей корпусов клапанов и на основании СТП 07.81-631-92 «Методика проведения теплового расчёта высокопараметрической арматуры клапанного типа».
4.3. Коэффициент теплоотдачи с наружной поверхности клапанов а3 = 7,8 ккал/м2×ч×С рассчитан по РД 302-07-122-89 «Методика расчета температурных полей арматуры для сред с температурой до 373 К (600 °С)».
4.4. Тепловая модель клапана запорного представлена на черт. 1 настоящего документа. На левой части чертежа показаны детали узла фланцевого соединения под номерами от 1 до 4, среднемассовые температуры которых необходимо определить для расчёта температурных напряжений.
4.5. Результаты расчётов среднемассовых температур деталей узла фланцевого соединения клапана приведены в таблице 5.
Температурное поле в стационарном режиме соответствует времени, равном t = 300 минут.
Нагрев 20 → 190 °С |
Охлаждение 190 → 20 °С |
|
Рабочая среда |
вода |
вода |
Температура рабочей среды Т, °C |
190 |
20 |
Давление рабочей среды Рр, МПа |
4,0 |
4,0 |
Скорость движения рабочей среды w, м/с |
5,0 |
5,0 |
Кинематическая вязкость рабочей среды v, м2/с |
0,15´10-6 |
0,55´10-6 |
Коэффициент теплопроводности рабочей среды, ккал/м×ч×С |
0,577 |
0,555 |
Критерий Прандтля Рr |
0,94 |
3,54 |
Диаметр условного прохода d, м |
100 |
100 |
Расчётные значения |
||
Нагрев 20 → 190 °С |
Охлаждение 190 → 20 °С |
|
Критерий Рейнольдса Re = w×d/v |
3,33´106 |
0,9´106 |
Критерий Нуссельта Nu = 0,021(Re)6,8 ´ Рr6,43 |
3354 |
2094 |
Коэффициент теплоотдачи a1 = Nu × l/d |
19352 |
11621 |
Среднемассовые температуры деталей фланцевого соединения клапана DN 100 в режимах 20 ↔ 190 °С. Скачок. Материал корпуса-стойки - ст. 20
РЕЖИМ НАГРЕВА от 20 до 190 °С |
РЕЖИМ ОХЛАЖДЕНИЯ от 191 до 20 °С |
|||||||||||
Номер деталей фланцевого соединения |
Температурные перепады Dt, °C |
Номер деталей фланцевого соединения |
Температурные перепады Dt, °C |
|||||||||
1 |
2 |
3 |
4 |
1 |
2 |
3 |
4 |
|||||
Номера точек тепловой модели |
Номера точек тепловой модели |
|||||||||||
1...10 22...26 |
25, 27, 28, 49, 59 |
46...50 56...60 |
32...34 38...40 |
Dt1-3 |
Dt4-2 |
1...10 22...26 |
25, 27, 28, 49, 59 |
46...50 56...60 |
32...34 38...40 |
Dt1-3 |
Dt4-2 |
|
0 |
20 |
20 |
20 |
20 |
0 |
0 |
186,5 |
167,5 |
167,5 |
176,5 |
19 |
3,9 |
0,5 |
46,1 |
21,4 |
20,6 |
39,6 |
25,4 |
18,2 |
170,5 |
166,8 |
167,2 |
163 |
3,3 |
-3,6 |
1 |
67,5 |
25,7 |
23,2 |
56,6 |
44,3 |
30,9 |
155,9 |
164,1 |
165,6 |
119,5 |
-9,8 |
-11,6 |
2 |
98,6 |
38,1 |
31,1 |
75,8 |
67,7 |
37,7 |
131,2 |
155,4 |
160,3 |
132,9 |
-29,1 |
-22,5 |
3 |
120,5 |
51 |
39,9 |
87,7 |
80,6 |
36,7 |
110,9 |
145,3 |
153,8 |
122 |
-42,9 |
-23,4 |
6 |
156,5 |
82,6 |
64,4 |
109,3 |
92,1 |
27,1 |
69,3 |
117,4 |
133,5 |
10,4 |
-61,2 |
-16,9 |
10 |
174,9 |
108,5 |
89,6 |
127,5 |
85,2 |
19 |
42 |
91,1 |
110,1 |
61,7 |
-68 |
-9,3 |
30 |
186,1 |
151,9 |
145,1 |
163,1 |
40,9 |
11,2 |
21,5 |
41,1 |
49,7 |
39,7 |
-28,2 |
-1,1 |
60 |
186,7 |
165,7 |
164,6 |
175 |
22,1 |
9,3 |
20,2 |
23,9 |
25,5 |
23,8 |
-5,3 |
-0,1 |
300 |
186,5 |
167,5 |
167,5 |
176,5 |
19 |
8,9 |
20 |
20,8 |
21,1 |
20,9 |
-1,1 |
-0,1 |
5. ЗАКЛЮЧЕНИЕ
В результате теплового расчёта узла фланцевого соединения клапана запорного DN 100 получены температуры в узловых точках тепловой модели фланцевого соединения, из которых определены среднемассовые температуры деталей соединения, а также осевые и радиальные перепады температур в нестационарном и стационарном режимах, необходимые при расчёте температурных напряжений.
СОДЕРЖАНИЕ
6. Методика теплового расчёта. 3 Приложение И. Пример оформления расчёта клапана запорного dn 100. 38 |