| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
АКЦИОНЕРНОЕ ОБЩЕСТВО ЗАКРЫТОГО ТИПА НАУЧНО-ПРОИЗВОДСТВЕННАЯ ФИРМА «ЦКБА» РУКОВОДЯЩИЙ ДОКУМЕНТ РД 26-07-25-97 РУКОВОДЯЩИЙ ДОКУМЕНТ
Дата введения 4 мая 1998 года 1. ОБЛАСТЬ ПРИМЕНЕНИЯНастоящий руководящий документ распространяется на запорную, регулирующую и предохранительную арматуру, предназначенную для установки в системах атомных электростанций, опытных и исследовательских ядерных реакторах, в системах и установках химических, нефтехимических и других производств с максимальным рабочим давлением до 20 МПа (200 кгс/см2) и с температурой рабочей среды до 873 К (600 °С) и устанавливает методику проведения теплового расчёта арматуры DN 10 до 1200 мм включительно.
ЛИСТ УТВЕРЖДЕНИЯ Руководящего документа «Методика расчёта температурных полей трубопроводной арматуры». РД 26-07-25-97 Заместитель Генерального директора АОЗТ НПФ «ЦКБА» по научной работе Ю.И. Тарасьев Главный конструктор В.В. Ширяев Заместитель директора НТИЦ Р.И. Хасанов Начальник отдела 161 А.А. Косарев Начальник лаборатории 154 В.В. Никитин Начальник сектора Г.И. Сергевнина Ведущий инженер-исследователь Н.С. Косых
2. НОРМАТИВНЫЕ ССЫЛКИВ настоящем руководящем материале использованы ссылки на следующие нормативные документы: РД 302-07-122-89. «Методика расчёта температурных полей арматуры для сред с температурой до 873 К (600 °С)». СТП 07.81-634-96 «Методические указания по проведению теплового расчёта задвижек для АЭС». СТП 07.81-631-92. «Методика проведения теплового расчёта высокопараметрической арматуры клапанного типа». 3. ОБЩИЕ ПОЛОЖЕНИЯ3.1. Тепловой расчёт арматуры должен проводиться с целью определения температурных полей, необходимых для оценки работоспособности её теплонапряжённых узлов и деталей, а также выбора конструкционных материалов с ограниченной теплостойкостью. 3.2. Методика теплового расчёта должна решать линейные и нелинейные задачи стационарной и нестационарной теплопроводности при заданных граничных условиях для Кmax = 300 узловых точек. 3.3. Тепловой расчёт должен включать в себя: а) задачу расчёта; б) исходные данные; в) условия расчёта; г) тепловую модель; д) программу теплового расчёта на электронных вычислительных машинах (ЭВМ); е) расчет; ж) заключение. 3.4. Программа теплового расчёта составлена на алгоритмическом языке «фортран» и реализована на персональных ЭВМ типа IВМ РС/AТ. Требования к РС: - наличие оперативной памяти не менее 536 КБ; - наличие дискового пространства для хранения информации до 2 МБ. 3.5. Исходные данные и условия расчёта методики изложены в разделах 4, 5 настоящего стандарта. 4. ИСХОДНЫЕ ДАННЫЕ4.1. Для проведения теплового расчёта необходимы следующие данные: а) температурное поле на начальный момент времени t = 0, К; б) закон изменения температуры рабочей и окружающей сред (линейный, нелинейный) или скорость изменения температуры рабочей и окружающих сред, b, К/ч; в) Тн - температура рабочей среды на начальный момент времени t = 0, К; г) Тк - температура рабочей среды на конец нестационарного режима, К; д) Та - температура окружающей среды, К; е) геометрические размеры расчетного изделия, м; ж) марки материалов деталей; и) теплофизические характеристики: l - коэффициент теплопроводности материалов деталей, вт/м×ч×К; a0 - коэффициент теплоотдачи со стороны рабочей среды, вт/м2×К, aа - коэффициент теплоотдачи от наружных поверхностей изделия в окружающую среду, вт/м2×К, r - плотность материалов, кг/м3; с - удельная теплоёмкость материалов, вт/кг×К; л) e1 - абсолютная невязка, К; м) e2 - относительная невязка, %; н) w - коэффициент релаксации. 5. УСЛОВИЯ РАСЧЁТА5.1. При проведении теплового расчета за начальные условия следует принять температурное поле на начальный момент времени t = 0. 5.2. Граничные условия на наружных и внутренних поверхностях изделия (соприкасающихся с рабочей и окружающей средами) следует задавать в виде температуры рабочей и окружающей сред и соответствующих коэффициентов теплоотдачи (g = а×Dt). 5.3. Контактные термические сопротивления между соприкасающимися деталями не учитываются. 5.4. При проведении теплового расчёта следует рассматривать осесимметричную двухмерную нестационарную задачу в цилиндрической системе координат. 5.5. Тепловой расчёт следует проводить при переменных теплофизических характеристиках материалов (l, с, r) на каждом шаге по времени Dt. 6. МЕТОДИКА ТЕПЛОВОГО РАСЧЁТА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.
Примечание. При отсутствии дополнительной тепловой модели 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. Пример оформления теплового расчета приведен в приложении И.
ПРИЛОЖЕНИЯ Приложение АNG1 - >085< - количество внутр. узлов основной тепловой модели ND - >085< - номер первого узла граничной области доп. тепловой модели NG2 - >085< - номер последнего узла граничной области доп. тепловой модели N - >085< - номер последнего узла внутр. области допол. тепловой модели NH - >301< - номер первой точки среды L - >302< - номер последней точки среды KNN - >004< - количество среднемассовых объектов (KNN > 0)
Приложение БKVD => 22 < -- количество шагов по времени ММ => 100 < -- количество итераций DT1 => 0,100 < -- абсолютная невязка DT2 => 0,0010 < -- относительная невязка RK => 1,800 < -- коэффициент релаксации
Приложение В
(1 СТРОКА - ПУСТАЯ, 2 И 3 - ИНФОРМАЦИОННЫЕ) *>05<* - количество используемых материалов
Приложение Г
Приложение Д= > 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 Приложение Ж
КОЛИЧЕСТВО МАТЕРИАЛОВ = 5 Т1 = .0 Т2 = 120.0 Т3 = 240.0 Т4 = 360.0
геометрические характеристики узлов NG1 = 85 ND = 85 NG2 = 85 N = 85 NH = 301 L = 302 KNN = 4
KVD = 22 ММ = 100 RK = 1,800 DT1 = ,100 DT2 = ,0010
ВРЕМЯ = .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
ТЕМРЕ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Тепловой расчет клапана запорного 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.
3.2. Температура окружающей среды - 20 °С. 3.3. Режимы для расчёта приведены в таблице 2.
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 минут.
Среднемассовые температуры деталей фланцевого соединения клапана DN 100 в режимах 20 ↔ 190 °С. Скачок. Материал корпуса-стойки - ст. 20
5. ЗАКЛЮЧЕНИЕ В результате теплового расчёта узла фланцевого соединения клапана запорного DN 100 получены температуры в узловых точках тепловой модели фланцевого соединения, из которых определены среднемассовые температуры деталей соединения, а также осевые и радиальные перепады температур в нестационарном и стационарном режимах, необходимые при расчёте температурных напряжений. СОДЕРЖАНИЕ
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
© 2013 Ёшкин Кот :-) |