Информационная система
«Ёшкин Кот»

XXXecatmenu

РЕСПУБЛИКАНСКИЕ СТРОИТЕЛЬНЫЕ НОРМЫ

 

ИНЖЕНЕРНЫЕ ИЗЫСКАНИЯ ДЛЯ СТРОИТЕЛЬСТВА.
СОСТАВЛЕНИЕ ПРОГНОЗА ИЗМЕРЕНИЙ
ТЕМПЕРАТУРНОГО РЕЖИМА ВЕЧНОМЕРЗЛЫХ
ГРУНТОВ ЧИСЛЕННЫМИ МЕТОДАМИ

 

РСН 67-87

Госстрой РСФСР

 

 

ГОСУДАРСТВЕННЫЙ КОМИТЕТ РСФСР ПО ДЕЛАМ СТРОИТЕЛЬСТВА

 

РАЗРАБОТАНЫ производственным объединением «Стройизыскания» Госстроя РСФСР совместно с институтом Фундаментпроект Минмонтажспецстроя СССР.

Исполнители: канд. геол.-мин. наук А.И. Левкович («Стройизыскания»), канд. геол.-мин. наук М.А. Минкин, инж. С.П. Дмитриева (Фундаментпроект), Ю.А. Попов («Стройизыскания»).

ВНЕСЕНЫ И ПОДГОТОВЛЕНЫ К УТВЕРЖДЕНИЮ производственным объединением по инженерно-строительным изысканиям («Стройизыскания») Госстроя РСФСР.

Вводятся впервые.

Государственный комитет РСФСР

по делам строительства (Госстрой РСФСР)

Республиканские строительные нормы

РСН 67-87

Госстрой РСФСР

Инженерные изыскания для строительства. Составление прогноза измерений температурного режима вечномерзлых грунтов численными методами

 

Настоящие Нормы устанавливают технические требования к производству работ по составлению прогноза измерений температурного режима грунтов для инженерно-геологического обоснования строительства новых, реконструкции и расширения действующих промышленных предприятий, зданий и сооружений, объектов сельскохозяйственного назначения, а также городов, поселков и сельских населенных пунктов на вечномерзлых грунтах.

Требования Норм не распространяются на составление прогноза для обоснования строительства гидротехнических, магистральных линейных и подземных сооружений, а также мостовых переходов.

Требованиями Норм следует руководствоваться при инженерно-геологическом обосновании строительства на участках талых грунтов, расположенных в пределах районов распространения вечномерзлых грунтов, а также при необходимости оценки динамики сезоннопромерзающего слоя в районах с сезонным промерзанием.

Составление прогноза изменения температурного режима грунтов необходимо производить в соответствии с требованиями к порядку составления прогноза мерзлотных инженерно-геологических условий, регламентированному РСН 31-83, а также с общими требованиями к прогнозу изменения инженерно-геологических условий, установленными действующими нормативными документами по изысканиям и проектированию оснований и фундаментов.

Внесены ПО «Стройизыскания» Госстроя РСФСР

Утверждены постановлением Государственного комитета РСФСР по делам строительства
от 20 августа 1987 г. № 152

Срок введения в действие
1 января 1988 г.

1. ОБЩИЕ ПОЛОЖЕНИЯ

1.1. Составление прогноза изменений температурного режима грунтов является необходимым элементом инженерно-геологического обоснования строительства (реконструкции, расширения) объектов народного хозяйства в районах распространения вечномерзлых грунтов.

Составление прогноза измерения температурного режима грунтов производится изыскательской организацией при участии проектной организации-заказчика.

Участие проектной организации заключается:

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

в представлении необходимых исходных данных для каждого варианта прогноза;

в совместном с изыскательской организацией обсуждении результатов прогноза.

1.2. Составление прогноза изменений температурного режима грунтов производится для ограниченного определенным образом грунтового массива путем последовательного расчета температурных полей в этом массиве, соответствующих любым заданным моментам времени от начала расчета. Температурный режим (совокупность последовательных температурных полей в этом массиве, соответствующих любым заданным моментам времени от начала расчета. Температурный режим (совокупность последовательных температурных полей) в грунтовом массиве рассчитывается как результат задаваемых на весь период расчета прогноза тепловых воздействий на верхней, боковых и нижней границах грунтового массива.

1.3. Расчет температурного режима грунтов производится для конкретных инженерно-геологических разрезов с учетом естественных изменений температурного режима воздуха и радиационного баланса дневной поверхности, естественных и техногенных условий теплообмена на поверхности грунтов (снежный и растительных покровы, насыпи, асфальтовые покрытия и т.п.) и техногенных источников и стоков тепла (здания и сооружения).

Расчет температурного режима грунтов производится на основе следующих данных:

материалов инженерно-геологических изысканий (инженерно-геологические разрезы с выделением классификационных разновидностей грунтов, физические и теплофизические свойства последних, естественные температуры грунтов и др.);

технического задания (пространственное размещение проектируемых объектов и проектный температурный режим в них, в том числе температурный режим в проветриваемых подпольях, вентилируемых насыпях и подвалах; при необходимости - вертикальная планировка территории; данные о заглублении проектируемых объектов, физические и теплофизические характеристики материала фундаментов, насыпей, полов и других теплоизоляционных покрытий и др.);

справочных материалов (температурный режим воздуха, радиационный баланс дневной поверхности, мощность и плотность снежного покрова и др.).

1.4. Расчетный срок (время) прогноза определяется расчетным сроком эксплуатации проектируемых объектов, для инженерно-геологического обоснования строительства которых составляется прогноз, и указывается в техническом задании.

Расчет прогноза может быть прекращен ранее в случаях стабилизации температурного режима грунтов в исследуемом грунтовом массиве. Под стабилизацией температурного режима в данном случае понимается наступление динамического равновесия температур в исследуемом массиве с учетом характера изменений тепловых воздействий на верхней границе массива.

1.5. Настоящие Нормы устанавливают правила расчета изменений температурного режима грунтов численным методом решения уравнения нестационарной теплопроводности с фазовыми переходами грунтовой влаги.

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

При составлении прогноза учитываются:

неоднородность состава, свойств и состояния грунтов в исследуемой грунтовой области;

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

изменение граничных тепловых условий во времени и пространство;

локальные источники и стоки тепла, расположенные внутри исследуемой грунтовой области.

1.6. Программа PROGNOZ имеет 5 модификаций, изложенных в п. 2.3., составлена на языке «ФОРТРАН-IV» и предназначена для использования на ЕС ЭВМ с операционной системой ОС ЕС MVT 6.1 и с памятью не менее 300 килобайт.

1.7. Пакет программ поставляется на магнитной ленте пользователя. Организация-держатель подлинника - институт «Фундаментпроект».

2. АЛГОРИТМ И СТРУКТУРА ПРОГРАММЫ

2.1. Математическая постановка задачи

Процесс распространения тепла в грунте в трехмерном пространстве описывается уравнением

                          (1)

где U(M,t) - температура грунта в точке М(х,y,z) в момент t;

H(M,U,t) - энтальпия (теплосодержание), отнесенная к единице объема грунта;

l(M,U)- коэффициент теплопроводности грунта;

F(M) - плотность тепловых внутренних источников и стоков в единице объема.

Энтальпия является функцией температуры, времени и координат. Так как рассмотренные ниже выводы верны для всех точек пространства, то будем рассматривать энтальпию только как функцию температуры. С учетом теплоты фазовых переходов в грунте, энтальпия равна:

где d(x) - дельта-функция Дирака

- объемная теплоемкость талого грунта;

- эффективная теплоемкость мерзлого грунта;

- удельная теплоемкость сухого грунта;

- плотность сухого грунта;

- удельная теплоемкость воды;

- удельная теплоемкость льда;

c- теплота фазовых переходов;

- суммарная влажность грунта в долях к весу абсолютно сухого грунта;

- незамерзшая вода при температуре x, принимается в виде

- коэффициенты, задающие кривую незамерзшей воды при

Задача рассматривается в параллелепипеде(т.е. ), на гранях которого задаются краевые условия:

на верхней границе                                                        (2)

а) температура окружающей среды:

или

б) теплопоток

или

в) теплообмен по закону Ньютона

на нижней границе                                                        (3)

температура окружающей среды

на боковых границах                                                     (4)

постоянный теплопоток

   

Начальная температура грунта известна по всем параллелепипеде D:

                   (5)

Требуется найти температуру грунта (функция U) непрерывную в D, удовлетворяющую уравнению (1), краевым условиям (2)-(4) и начальному условию (5).

2.2. Алгоритм решения задачи

Предлагаемым алгоритмом решения задача (1)-(5) решается энтальпийным конечно-разностным методом по явной

В прямоугольнике D вводится произвольная прямоугольная неравномерная разностная сетка с шагамии временная сетка с шагами:

        (6)

Применяя интегро-интерполяционный метод (метод баланса) построения однородных разностных схем, разностное уравнение, аппроксимирующее уравнение (1) по явной схеме на сетке (6) будет иметь вид

                           (7)

           

        (8)

где- теплосодержание элементав момент времени;

- теплопоток соответственно через верхнюю, нижнюю и четыре боковые грани каждого элемента;

- изменение энергии внутренних источников в объеме за время.

Для краевых условий I рода соответствующееграничных элементов равны:

где;

для краевых условий II рода

гдезаданная величина теплового потока;

для краевых условий III рода

,

где- температура внешней среды;

- термическое сопротивление;

- коэффициент конвективного теплообмена.

При решении задачи на шаге e по известной температуре и энтальпии He в каждом элементе определяются энтальпия на шаге e+1 по формуле (8). Так как существует взаимно однозначное соответствие между энтальпией и температурой, то находим температуру каждого элемента на слое e+1. Затем определяем энтальпию He+2 и т.д.

Расчетные формулы для определения энтальпии и температуры по известной энтальпии следующие:

По определению полная энтальпия

где- теплоемкость грунта (рис. 1)

 

 

Рис. 1. График изменения теплоемкости по температуре

— — — — —    без учета фазовых переходов в спектре отрицательных температур

____________    то же с учетом фазовых переходов в спектре отрицательных температур

Все дальнейшие расчеты проводятся с учетом фазовых переходов в спектре отрицательных температур. Незамерзшая вода учитывается в виде .

Критические значения энтальпии:

.

Здесь и далее- условная температура в градусах Цельсия, при которой теплосодержание (энтальпия) принимается равной нулю. (Абсолютный нуль К=-273°С).

.

Энтальпия определяется по следующим формулам (рис. 2):

если U0>U*+e, то

;

если U0<U*+e, то

.

 

 

Рис. 2. График изменения энтальпии по температуре

По известной энтальпии температура определяется:

если H(U0)>H2 , то

;

если, то

;

если, то U* определяется линейной интерполяцией по значениям табличной функции энтальпии, составляемой для каждого слоя грунта.

2.3. Структура программы

Алгоритм, описанный выше, реализован в программе PROGNOZ. Связь между частями программы представлена на рис. 3.

 

 

Рис. 3. Межмодульные связи в программе PROGNOZ

Программа PROGNOZ имеет 5 модификаций:

1. PROGNOZ-3S

Решается задача по трехмерной расчетной схеме. Главная (управляющая) программа MAIN составляется для каждого варианта прогноза. В программе MAIN вводятся исходные массивы начальной температуры и некоторые другие исходные данные. Размерности массива, индивидуальные для каждого варианта прогноза, задаются в операторах описания.

Из программы MAIN происходит обращение к подпрограмме WNS. В ней определяются коэффициенты А, В, и С для вычисления количества незамерзшей воды в каждом слое грунта, по формуле где U£0 - температура грунта. Для контроля выдается на печать таблица получаемых значений Ww от U при.

Из подпрограммы WNS получены массивы А(N), В(N), С(N) передаются в MAIN.

В подпрограмме GRUND определяются одно-, двух- или трехмерные температурные поля в заданном массиве грунтов. Здесь используются подпрограмма-функция HF и подпрограмма FIND. HF (N,U1) вычисляет энтальпию 1 м3 грунта в N-ом слое при температуре U1<URK (критическая температура фазовых переходов). Для расчета функция HF использует физико-механические и теплофизические свойства грунтов, введенные подпрограммой BLOCK DATA.

Подпрограмма FIND (HN1, N1,UL1) определяет температуру UL1 в слое грунта с номером N1 по энтальпии HN1, если HN1<H1(N1), т.е. еслисоответствует энтальпии мерзлого грунта. FIND использует общую область TABL, в которой содержатся таблицы значений энтальпии по температуре от 0 до -50°С по всем грунтовым слоям, вычисленным с помощью функции HF.

Подпрограмма ISTO (i, j, K, T, F) учитывает влияние источников тепла и охлаждающих установок, находящихся внутри массива грунта. Результат работы подпрограммы присваивается простой переменной F, где F - тепло, выделяемое внутренними источниками (стоками) в единицу времени (час) в фиксированном элементе (i,j,K) грунтового массива в момент времени T.

Алгоритм подпрограммы может быть составлен как по экспериментальным данным, так и по теоретическим расчетам применительно к конкретному случаю.

Для обращения к подпрограмме ISTO в подпрограмме GRUND (рекомендуемое приложение 2) следует заменить оператор А=0.0 с номером 0000399 на обращение к подпрограмме ISTO:

CALL ISTO (i,j,K,T,F).

В подпрограмме GRUND выводятся на печать значения рассчитанного температурного поля.

Пример текста подпрограммы ISTO приведен в рекомендуемом приложении 2.

2. PROGNOZ-3N

Задача решается по трехмерной расчетной схеме. Отличие от PROGNOZ-3S в том, что вся исходная информация вводится только с перфокарт, изменений текста не требуется. Максимальное количество элементов в расчетной области - 3000. Ввод данных осуществляется в управляющей программе MAIN3N и подпрограмме GLVN3N.

Подпрограмма BLOCK DATA отсутствует.

3. PROGNOZ-2S

Эта подпрограмма отличается от PROGNOZ-3S тем, что рассчитывает только двух- и одномерные температурные поля. Программы MAIN и BLOCK DATA составляются для каждого варианта прогноза. Расчет двухмерных задач по модификации PROGNOZ-2S сокращает время расчета по отношению к PROGNOZ-3S приблизительно 1,5 раза.

4. PROGNOZ-2N

Аналогично PROGNOZ-2S рассчитывается только двух- и одномерные температурные поля. Вся исходная информация вводится только с перфокарт. Ввод данных осуществляется в подпрограмме MAIN2N и подпрограмме GLVN2N. Подпрограмма BLOCK DATA отсутствует. Максимальное количество элементов в расчетной области - 3000.

5. PROGNOZ-L

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

Ввод данных аналогичен PROGNOZ-2N. Максимальное количество элементов (с учетом фиктивных) в расчетной области - 3000.

3. ПОДГОТОВКА ИСХОДНЫХ ДАННЫХ ДЛЯ РАСЧЕТА

3.1. Расчетная область

Расчетная область исследования может быть задана трехмерной, двухмерной или одномерной. Так как в грунтах теплообмен всегда пространствен, т.е. трехмерен, то в трехмерной расчетной области может решаться любая задача прогнозирования. В этом случае любая точка М грунта определяется координатами (x, y, z) в прямоугольной системе координат. В ряде случаев задача может быть сведена к двухмерной и даже одномерной, что более выгодно и удобно при расчете на ЭВМ.

Алгоритм составлен и может быть применен для расчетных областей, ограниченных прямоугольной верхней границей, в контурах которой фиксированы в плане тепловыделяющие объекты. В трехмерном случае область исследования параллелепипед, в двухмерном - прямоугольник и в одномерном - прямая. Для создания объема в двухмерной и одномерной областях задаются по направлениям, в которых размеры не определяются, единичные размеры, равные 1 м.

Область исследования фиксируется прямоугольной системой координат, с центром 0 в одной из вершин параллелепипеда и направлениями осей: z - в вертикальном направлении, x и y в горизонтальных. Направление горизонтальных осей x и y совпадает с осевыми линиями тепловыделяющих объектов. Величина ребер параллелепипеда определяет размеры области исследования. Его грани являются границами области исследования.

Различаются верхняя, нижняя и четыре боковые границы (соответственно верхняя, нижняя и четыре боковых грани параллелепипеда).

Теплопоток через боковые границы области исследования должен быть построен во времени и по всей плоскости грани. Исходя из этого, определяются размеры области исследования.

На практике чаще применяется условие отсутствия теплопотока через боковые границы. Тогда для тепловыделяющих сооружений расстояние от их контура в плане до боковой границы области исследования должно быть в 2,5-3 раза больше, чем размеры самого здания в том же направлении.

Положение нижней границы (глубина области исследования) выбирается ниже глубины распространения годовых колебаний температуры в грунтах и зависит от интенсивности тепловыделения на верхней границе, а также срока прогнозного расчета.

При естественных условиях теплообмена на верхней границе, или без наличия тепловыделяющих сооружений эта глубина должна быть не менее 15-20 м. В остальных случаях она должна быть на 10 м больше максимальной глубины чаши оттаивания, определенной расчетом по СНиП II-18-76.

При определении области исследования следует исходить также из экономии времени счета на ЭВМ. Для этого следует использовать возможную симметричность задачи, а также возможность сведения трехмерной области с двухмерной или одномерной.

Если задача симметрична, или может быть сведена к симметричной, то в трехмерной задаче плоскость (или плоскости) симметрии, а в двухмерной ось (или оси) симметрии должны быть границами области исследования с условием отсутствия теплообмена на них.

3.2. Разбивка области исследования

Область исследования разбивается на прямоугольные элементы произвольных размеров (рис. 4).

 

 

Рис. 4. Пример разбивки области исследования

Разбивка производится горизонтальными и вертикальными плоскостями (разбивочными плоскостями), параллельными соответствующим границам области в следующем порядке:

на верхней и нижней границах области исследования выделяются зоны с различным режимом теплообмена грунта с окружающей средой. На каждом из этих границ может быть выделено до 9 таких зон;

границы выделенных зон с различным режимом теплообмена должны быть параллельны боковым границам области исследования. Если граница выделенной зоны не отвечает этому требованию, то она заменяется ступенчатой линией со «ступеньками», аппроксимирующими ее контур и удовлетворяющими указанному требованию;

после построения границ зон через каждую прямую, образующую «ступеньки», как на верхней, так и на нижней границах области исследования проводятся разбивочные плоскости;

криволинейные границы литологических разностей заменяются ступенчатыми, с плоскостями образующих их «ступенек», параллельными границам области исследования;

проводятся разбивочные плоскости через плоскости полученных «ступенек»;

проводятся остальные разбивочные плоскости так, чтобы в зоне возможных фазовых переходов грунтовой влаги размеры ребер получаемых прямоугольных элементов в направлении теплообмена составляли от 0,5 до 2 м. Размер ребер тех же элементов в перпендикулярных направлениях может быть соответственно в 2-3 раза больше. Вне зоны возможных фазовых переходов и влияния тепловыделяющих объектов размер элементов в направлении теплообмена может быть постепенно увеличен до 6-8 м.

По каждому направлению число элементов не должно превышать 100.

В результате проведенного разбиения область исследования состоит из прямоугольных элементов разных размеров. Каждая выделенная литологическая разность представлена целым числом элементов. В контурах зон с различным режимом теплообмена на верхней и нижней границах содержится целое число элементов (в плане).

Полученное разбиение остается постоянным на все время проведения расчета.

За каждым элементом закрепляется трехзначный индекс (i, j, к), где i определяет номер элемента по направлению оси z, j - по направлению оси х и к - по направлению оси у.

Фиксируется число расчетных элементов по направлению осей.

Обозначим их соответственно NI, NJ, NK. Произведение P=NI´NJ´NK определяет число расчетных элементов во всей области исследования.

При решении двухмерной задачи NK =1. При решении одномерной NK =1 и NJ =1.

3.3. Начальные условия

В каждом полученном элементе области исследования задаются начальные условия (U*), т.е. температуры грунта, соответствующие времени начала расчета. Эти температуры назначаются в геометрических центрах элементов и считывается, что каждый элемент характеризуется температурой его центра. Последнее положение также относится к любым температурным полям, полученным в процессе расчета задачи.

Начальное распределение температур в области исследований назначается по данным термокаротажных работ. Время проведения этих работ считается моментом начала расчета.

Начальная температура (U0) ни в одном из элементов не может назначаться равной температуре начала фазовых переходов (U*). В этом случае задается температура (U*)±e, где e=0,001°С.

Начальная температура грунта (U0) записывается в табл. 1 по «разрезам», сделанным параллельно оси j в области исследования для каждого фиксированного значения K, где 1£K£NK.

Таблица 1

K=1

i

j

j =1

j =2

 

j =NJ

i=1

U (1, 1, 1)

U (1, 2, 1)

· · ·

U (1, NJ, 1)

i=2

U (2, 1, 1)

U (2, 2, 1)

· · ·

U (2, NJ, 1)

i=3

U (3, 1, 1)

U (3, 2, 1)

· · ·

U (3, NJ, 1)

·

·

·

·

·

·

·

·

·

 

·

·

·

i =Ni

U (Ni, 1, 1)

U (1, 2, 1)

· · ·

U (Ni, NJ, 1)

K=2

i

j

j =1

j =2

 

j =NJ

i=1

U (1, 1, 2)

U (1, 2, 2)

· · ·

U (1, NJ, 2)

i=2

U (2, 1, 2)

U (2, 2, 2)

· · ·

U (2, NJ, 2)

i=3

U (3, 1, 2)

U (3, 2, 2)

· · ·

U (3, NJ, 2)

·

·

·

·

·

·

·

·

·

 

·

·

·

i =Ni

U (Ni, 1, 2)

U (Ni, 2, 2)

· · ·

U (Ni, NJ, 2)

и т.д.

K=NK

i

j

j =1

j =2

 

j =NJ

i=1

U (1, 1, NK)

U (1, 2, NK)

· · ·

U (1, NJ, NK)

i=2

U (2, 1, NK)

U (2, 2, NK)

· · ·

U (2, NJ, NK)

i=3

U (3, 1, NK)

U (3, 2, NK)

· · ·

U (3, NJ, NK)

·

·

·

·

·

·

·

·

·

 

·

·

·

i =Ni

U (Ni, 1, NK)

U (Ni, 2, NK)

· · ·

U (Ni, NJ, NK)

3.4. Граничные условия

Граничные условия задаются на каждое границе области исследований на все время расчета температурного режима.

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

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

На каждой из боковых границ области исследования задается по одному граничному условию.

Различают граничные условия I, II и III рода.

Условие I рода - известна температура грунта на поверхности границы;

Условие II рода - известен теплопоток через поверхность границы.

Условие III рода - теплообмен через границу области исследования определяется температурой окружающей среды и коэффициентом теплообмена, т.е. происходит по закону Ньютона.

Q=s(Un-tB),

где Q - количество приходящего к поверхности грунта или уходящего от него тепла, кДж (ккал);

s- коэффициент теплоотдачи с поверхности грунта, Вт/(м2·°С), ккал/ /(м2´´°С·ч);

Un - температура поверхности грунта, °С;

tB - температура окружающей среды, °С.

Граничные условия верхней границы изменяются с периодом повторения Тпер, назначаемым в исходных данных расчета. На практике наиболее часто употребляется Тпер=1 год=8760 часов.

Период счета разбит на интервалы. Если период год, то интервалов 12 (12 месяцев) или 36 (36 декад). Количество интервалов задается в исходных данных переменной Li. На верхней границе области исследования, в каждой выделенной зоне с различным режимом теплообмена грунта с внешней средой, для каждого интервала периода счета задаются следующие средние значения за этот интервал величин, определяющих теплообмен в зависимости от рода краевых условий в зоне:

для условий I рода задается температура поверхности грунта, °С;

для условий II рода задается величина теплопотока; кДж (ккал);

для условий III рода задаются:

tB - температура окружающей среды, °С;

Rпов - термическое сопротивление поверхности грунта, м2·°С/Вт, м2·ч´ ´°С/ккал;

aпов- коэффициент теплообмена, Вт/(м2·°С), ккал/(м2·ч·°С).

На любой момент времени от начала счета значения этих величин определяются линейной интерполяцией между значениями их в соседних интервалах.

На нижней границе задаются условия теплообмена I рода. При этом в каждой выделенной зоне значения температуры могут быть только постоянными на все время проведения расчета.

На каждой их четырех боковых границ задаются условия II рода. Теплопоток на каждой границе постоянен на все время проведения расчета и по всей плоскости грани. Для граничных условий верхней границы заполняется табл. 2. При этом каждой выделенной зоне присваивается постоянный номер от 1 до М, где М £ 9.

Таблица 2

Номер зоны

Год краевых условий в зоне

Характеристики

Месяц (или декада)

I

II

¼

XII

1

III

Температура окружающей среды t, °С

 

 

 

 

Термическое сопротивление R, м2·°С/Вт (м2·ч·°С/ /ккал)

 

 

 

 

Коэффициент теплообмена a, Вт/(м2·°С), ккал/(м2´ ´ч·°С)

 

 

 

 

2

I

Температура на поверхности грунта t, °С

 

 

 

 

3

II

Теплопоток a, кДж (ккал)

 

 

 

 

4

III

Температура окружающей среды t, °С

 

 

 

 

Термическое сопротивление R, м2·°С/Вт (м2·ч·°С/ /ккал)

 

 

 

 

Коэффициент теплообмена a, Вт/(м2·°С), ккал/(м2´ ´ч·°С)

 

 

 

 

 

 

 

 

 

 

 

 

 

и т.д. для каждой выделенной зоны

После закрепления номеров за каждой зоной краевых условий верхней границы заполняется табл. 3. Каждое значениеэтой таблицы имеет вид:

,

где KOL=0, если во всех элементах (i, j, K), где 1£i£NI, а j, K фиксированы, внутренние источники отсутствуют;

KOL=1, если внутренние источники присутствуют хотя бы в одном из элементов (i, j, K), где 1£i£NI, а j, K фиксированы;

NP - номер типа источника, NP=1,2,...,9. Если источник отсутствует NP=1;

MD - номер зоны краевых условий верхней границы, в которой находится элементобласти исследования;

MDK - номер зоны краевых условий нижней границы для соответствующего элементаобласти исследования.

Таблица 3

j

K

K=1

K=2

· · ·

K =NK

j=1

KL (1, 1)

KL (1, 2)

· · ·

KL (1, NK)

j=2

KL (2, 1)

KL (2, 2)

· · ·

KL (2, NK)

·

·

·

·

·

·

·

·

·

 

·

·

·

j=NJ

KL (NJ, 1)

KL (NJ, 2)

· · ·

KL (NJ, NK )

Для двухмерной задачи табл. 3 состоит из одного столбца для K=1, для одномерной из одного числа KL(1,1).

3.5. Физические и теплофизические свойства

Расчетная область может содержать до 10 литологических разностей, характеризующихся следующими свойствами, учитывающимися алгоритмом:

- плотность сухого грунта кг/м3;

- весовая влажность грунта, доли единицы;

- удельная теплоемкость скелета, кДж/(кг·°С), ккал/(кг·°С);

- коэффициент теплопроводности грунта в талом состоянии Вт/(м´ ´°С), ккал/(м·ч·°С);

U0 - температура начала фазовых переходов, °С;

WW - количество незамерзшей воды в диапазоне отрицательных температур, доли, задаваемой кривой

,

где A, B, C - коэффициенты, определяющие количество незамерзшей воды в данной литологической разности.

Все вышеуказанные свойства должны быть определены для каждой литологической разности по лабораторным данным, справочным или нормативным документам. Они остаются неизменными для всего времени счета.

Составляется табл. 4. При этом каждой литологической разности присваивается постоянный номер от 1 до NSL, где NSL £ 10.


Таблица 4

Номер литологической разности

,

кг/м3

,

ккал/(кг·°С)

,

доли

, Вт/(м·°С),

ккал/(м·ч·°С)

, Вт/(м·°С),

ккал/(м·ч·°С)

, °С

WW

доли

1

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

·

·

·

 

 

 

 

 

 

 

 

 

 

 

 

 

 


В графе WW помещается информация о незамерзшей воде в данном слое. Если известна кривая, то есть уже известны коэффициенты A, B, C ее задающие, то записываются значения коэффициентов. Если известны лабораторные данные, то приводятся они. Например: и так далее не менее трех значений. Если в расчете WW=0, то коэффициенты, задающие кривую незамерзшей воды в данном слое будут равны:.

После закрепления номеров за каждой литологической разностью заполняется табл. 5, структура которой повторяет табл. 1. Таблицей 5 за каждым элементом (i, j, K) , полученным при разбивке области исследования, закрепляется номер той литологической разности, которой он принадлежит.

Таблица 5

K=1

i

j

j =1

j =2

 

j =NJ

i=1

N (1, 1, 1)

N (1, 2, 1)

· · ·

N (1, NJ, 1)

i=2

N (2, 1, 1)

N (2, 2, 1)

· · ·

N (2, NJ, 1)

·

·

·

·

·

·

·

·

·

 

·

·

·

i =Ni

N (Ni, 1, 1)

N (1, 2, 1)

· · ·

N (Ni, NJ, 1)

K=2

i

j

j =1

j =2

 

j =NJ

i=1

N (1, 1, 2)

N (1, 2, 2)

· · ·

N (1, NJ, 2)

i=2

N (2, 1, 2)

N (2, 2, 2)

· · ·

N (2, NJ, 2)

·

·

·

·

·

·

·

·

·

 

·

·

·

i =Ni

N (Ni, 1, 2)

N (Ni, 2, 2)

· · ·

N (Ni, NJ, 2)

и т.д.

K=NK

i

j

j =1

j =2

 

j =NJ

i=1

N (1, 1, K)

N (1, 2, K)

· · ·

N (1, NJ, NK)

i=2

N (2, 1, K)

N (2, 2, K)

· · ·

N (2, NJ, NK)

·

·

·

·

·

·

·

·

·

 

·

·

·

i =Ni

N (Ni, 1, K)

N (Ni, 2, K)

· · ·

N (Ni, NJ, NK)

Каждое значение N(i, j, K) в табл. 5 может быть только целым числом от 1 до NSL, где NSL£10 - число выделенных литологических разностей.

3.6. Дополнительные данные, необходимые
для проведения расчета

Для проведения расчета необходимо также определить значение следующих величин, задействованных алгоритмом.

Шаг по времени. Размеры элементов и шаг по времени связаны между собой следующим соотношением:

где Dt - шаг по времени, ч;

i - номер грунтовой разности;

NSL - число выделенных грунтовых разностей;

- теплоемкость мерзлого грунта i-го слоя при , кДж/(м3·°С), ккал/(м3·°С);

- количество незамерзшей воды, доли;

- удельная теплоемкость сухого грунта, кДж/(кг·°С), ккал/(кг·°С);

- удельная теплоемкость льда, кДж/(кг·°С), ккал/(кг·°С);

n- коэффициент объемного расширения, n=1,1;

- весовая влажность грунта, доли;

- коэффициент теплопроводности мерзлого грунта, Вт/(м·°С), ккал/ /(м·ч·°С);

n - мерность области исследования, 1£n£3;

Dh - размер ребра элемента.

Шаг по времени может быть изменен в процессе счета NT раз, где NT³1.

Время измерения расчетного шага от начала расчета задается одномерным массивомразмерностью.

где- время измерения шага от начала счета в часах.

Отсчет происходит от моментов начала счета, т.е. от 0.

Определяется Ткон (час) - время окончания решения задачи. Задаются Т1, Т2, Т3, ..., ТР - время выдачи результатов расчета на печать (в часах). Время последней выдачи результата TP обязательно меньше TКОН на один шаг по времени. В противном случае счет может прекратиться раньше, чем будет получен последний результат.

Определяются интервалы, через которые происходит очередная выдача результатов на печать.

и т.д.

Для последующего ввода в ЭВМ формируется одновременный массив HPRI размерностью INT. Размерность INT массива зависит от того, с постоянным шагом выдаются результаты расчета или нет.

INT=2 - выдача производится с постоянным шагом;

INT=P+1 - выдача производится с переменным шагом.

Если INT=2, то

Если INT=P+1, то

и т.д.

Задание печати распределения температур

Для удобства обработки информации, получаемой в результате расчета, фиксируются необходимые для анализа температурного режима «линии резервов», по которым выдается на печать распределение температур.

Положение «разрезов», параллельных плоскости ZOX, задается координатойэлементов, составляющих этот разрез.

Формируется массив целых чисел KP размерностью KW, где KW - число выдаваемых на печать «разрезов» по направлению. Если по этому направлению ни одного разреза выдавать на печать не надо, то KW=1 и KP(1)=0. В противном случаеи т.д. , где

Положение «разрезов», параллельных плоскости ZOY, задается координатой j элементов, составляющих этот разрез.

Формируется массив целых чисел jP размерностью jW, где jW - число выдаваемых на печать «разрезов» по направлению. Если по этому направлению ни одного разреза выдавать на печать не надо, то jW=1 и jP(1)=0. В противном случаеи т.д. где

Положение «разрезов», параллельных плоскости XOY, задается координатой i элементов, составляющих этот разрез.

Формируется массив целых чисел iP размерностью iW, где iW - число выдаваемых на печать «разрезов» по направлению.

Если по этому направлению ни одного разреза выдавать на печать не надо, тои. В противном случаеи т.д. , где.

4. ПРОВЕДЕНИЕ РАСЧЕТА

4.1. Вызов и загрузка программы в ЭВМ

Каждая из модификаций программы PROGNOZ состоит из нескольких модулей, объединенных главной программой MAIN. Подпрограммы-модули GRUND, WNS, HF, FIND содержатся в личной библиотеке исходных модулей (или библиотеке загрузочных модулей).

Ниже приводится пример пакета вызова программы PROGNOZ-3S из библиотеки исходных модулей. Задание требует не менее 300 килобайт памяти и использует оптимизирующий транслятор (ФОРТРАН-ОП) системы ОС ЕС.

Здесь:

имя - биб - имя личной библиотеки исходных модулей;

имя - тома - имя тома, на которой размещена библиотека.

// JOB

// EXEC FORTHCLG

// FORT. SYSIN DD DDNAME=PKAR

//                DD DSN=<ИМЯ - БИБ>(GRUND),DISP=SHR,

//                          VOL=SER=<ИМЯ-ТОМА>,UNIT=SYSDA

//                DD DSN=<ИМЯ - БИБ>(WNS),DISP=SHR,

//                          VOL=SER=<ИМЯ-ТОМА>,UNIT=SYSDA

////              DD DSN=<ИМЯ - БИБ>(HF),DISP=SHR,

//                          VOL=SER=<ИМЯ-ТОМА>,UNIT=SYSDA

//                DD DSN=<ИМЯ - БИБ>(FIND),DISP=SHR,

//                          VOL=SER=<ИМЯ-ТОМА>,UNIT=SYSDA

// PKAP DD*

<текст программы MAIN на перфокартах>

<текст подпрограммы BLOCK DATA на перфокартах>

// GO. SYSIN DD*

<исходные данные>

/*

//

Текст подпрограммы для модификации PROGNOZ-3S приведен в обязательном приложении 3.

4.2. Составление программы MAIN и подпрограммы BLOCK DATA

Программа MAIN в модификациях PROGNOZ-3S и PROGNOZ-2S составляется для каждого варианта расчета. Пример программы MAIN для расчета контрольного примера (приложение 1) приведен в обязательном приложении 3.

Программа MAIN содержит обращение к основной подпрограмме GRUND алгоритма. Перед обращением к подпрограмме GRUND в MAIN должны быть определены все формальные аргументы подпрограммы, приведенные в табл. 6.

Таблица 6

Идентификатор

Фактическое значение

NSL

Число выделенных литологических разностей (NSL£10)

NI

Число элементов по направлению i, (NI£100)

NJ

Число элементов по направлению j, (Nj£100)

NK

Число элементов по направлению K, (NK£100)

NT

Число изменений шага по времени (NT£10)

LI

Число интервалов времени в периоде, во время которых задаются средние значения характеристик внешней среды на верхней границе (Li£36)

M

Число зон с различными краевыми условиями на верхней границе (M£9)

M1

Число зон с различными краевыми условиями на нижней границе (M1£9)

TKON

Время окончания счета (в часах)

A1(NSL), B1(NSL)

Одномерные массивы размерностью NSL

C1(NSL)

Коэффициенты, задающие кривую незамерзшей воды по слоям

N(NI,NJ,NK)

Трехмерный массив размерностью NI´NJ´NK

Номера слоев грунта в соответствующих элементах (см. табл. 5)

KL(NJ,NK)

Двухмерный массив размерностью NJ´NK

Массив номеров зон краевых условий по элементам верхней и нижней границ (см. табл. 3)

U(NI,NJ,NK)

Трехмерный массив размерностью NI´NJ´NK

Начальная температура грунта по элементам (см. табл. 1)

IW

Число выдаваемых на печать «разрезов» по направлению i (см. п.3.6)

JW

Число выдаваемых на печать «разрезов» по направлению j (см. п.3.6)

KW

Число выдаваемых на печать «разрезов» по направлению K (см. п.3.6)

INT

Размер массива HPRI (см. п.3.6.)

HPRI(INT)

Одномерный массив размерности INT

Задает интервалы выдачи на печать температуры грунта (час) (см. п.3.6.)

H(NI,NG,NK)

Трехмерный массив размерностью NI´NJ´NK

Значение энтальпии по элементам на текущем временном слое

UN(NI,NG,NK)

Трехмерный массив размерностью NI´NJ´NK

Температура грунта по элементам на следующем временном слое

HNDV(NI,NG,NK)

Трехмерный массив размерностью NI´NJ´NK

Значение энтальпии по элементам на следующем временном слое

Nsi

Nsi =1, если единица измерения количества теплоты - кДж

 

Nsi ¹1, если единица измерения количества теплоты - ккал

Оператором DIMENSION определяется размерность массивов H, UN, HNDV, U, A1, B1, C1. Ввод массивов N, KL, U, HPRI осуществляется операторами описания типа (INTEGER, REAL), а также операторами присваивания.

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

Обращение к подпрограммам WNS (NSL, A1, B1, C1) присутствует в том случае, если коэффициенты A1, B1, C1, определяющие кривую незамерзшей воды по грунтовым разностям, неизвестны. В противном случае эти коэффициенты вводятся операторами присваивания или операторами описания типа.

Обращение к подпрограмме GRUND присутствует обязательно и имеет вид:

CALL GRUND (NSL, NI, NJ, NK, NT, LI, M, M1, TKON, A1, B1, C1, N, KL, U, IW, JW, KW, INT, HPRI, H, UN, HNOV, NSI).

Простые переменные NSL, NI, NJ, NK, NT, LI, M, M1, TKON, IW, JW, KW, INT, NSI при обращении к подпрограмме заменяются своими фактическими значениями (см. табл. 6).

Массивы H, UN, HNOV определяются только размерностью в операторе DIMENSION.

В подпрограмме BLOCK DATA (см. приложение 3) оператором DATA вводятся значения по слоям следующих физических и теплофизических характеристик грунта соответствующими идентификаторами:

ВСК (10) - плотность сухого грунта, кг/м3;

ССК (10) - удельная теплоемкость сухого грунта, кДж/(кг·°С), ккал/(кг´ ´°С);

WC (10) - суммарная влажность грунта в долях к весу абсолютно сухого грунта (доли);

UKR (10) - температура начала фазовых переходов.

Значения характеристик берутся из табл. 4.

Текст программы MAIN и подпрограммы BLOCK DATA записывается на стандартных бланках для записи текста на FORTRANe с последующей набивкой на перфокарты.

4.3. Подготовка к вводу группы данных с перфокарт

Начальные данные, вводимые с перфокарт описаны в порядке ввода в табл. 7 и 8.

По табл. 7 вводятся данные, необходимые для расчета коэффициентов кривой незамерзшей воды подпрограммой WNS, для каждой выделенной грунтовой разности. Если в программе MAIN обращения к подпрограмме WNS нет, то числовые данные для ввода готовятся только по табл. 8.

Таблица 7

Порядковый номер оператора ввода

Идентификатор

Формат

Единица измерения

Значение

1

T(3,NSL)

1256.2

°C

По три значения температуры, для которых известно количество незамерзшей воды, для каждой грунтовой разности в порядке их нумерации (см. табл. 4)

2

W(3,NSL)

1256.2

доли ед.

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

Таблица 8

Порядковый номер оператора ввода

Идентификатор

Формат

Единица измерения

Значение

3

TLM(NSL)

1356.2

Вт/(м·°С),

ккал/(м·ч·°С)

Одномерный массив размерностью NSL. Теплопроводность мерзлого грунта по слоям в порядке их нумерации (см. табл. 4)

 

TLT(NSL)

 

Вт/(м·°С),

ккал/(м·ч·°С)

Одномерный массив размерностью NSL. Теплопроводность талого грунта по слоям в порядке их нумерации (см. табл. 4)

 

WI(N)

 

м

Одномерный массив размерностью NI. Размер элементов по направлению I

 

NJ(NJ)

 

м

Одномерный массив размерностью NJ. Размер элементов по направлению J

 

HK(NK)

 

м

Одномерный массив размерностью NK. Размер элементов по направлению K

 

SHT(NT)

 

ч

Одномерный массив размерностью NT. Размер шага по времени (см. п.3.6)

4

TR(NT)

6F12.2

ч

Одномерный массив размерностью NT. Граница (время) измерения шага по времени (см. п.3.6)

5

TB(LI,M)

12F6.2

°С

Двухмерный массив размерностью LI´M.

Температура внешней среды на верхней границе области исследования при периодическом ее изменении по времени и зонам краевых условий в порядке их нумерации (см. табл. 2)

 

RB(LI,M)

 

м2·°С/Вт,

м2·ч·°С/ккал

Двухмерный массив размерностью LI´M.

Термическое сопротивление поверхности на верхней границе при периодическом изменении по времени и зонам краевых условий в порядке их нумерации (см. табл. 2)

 

AB(LI,M)

 

Вт/(м2·°С),

ккал/(м2·ч·°С)

Двухмерный массив размерностью LI´M.

Коэффициент конвективного теплообмена на верхней границе при периодическом изменении по времени и зонам краевых условий в порядке их нумерации (см. табл. 2)

 

TBIK(M1)

 

°С

Одномерный массив размерностью M1.

Температура внешней среды по зонам нижней границы в порядке их нумерации

 

TBJ

 

Вт/(м2·°С),

ккал/(м2·ч·°С)

Простая переменная. Теплопоток на левой боковой границе области исследования

 

TBJK

 

Вт/(м2·°С),

ккал/(м2·ч·°С)

Простая переменная. Теплопоток на правой боковой границе области исследования

 

TBK

 

Вт/(м2·°С),

ккал/(м2·ч·°С)

Простая переменная. Теплопоток на ближней боковой границе области исследования

 

TBKK

 

Вт/(м2·°С),

ккал/(м2·ч·°С)

Простая переменная. Теплопоток на дальней боковой границе области исследования

6

MR(M)

2713

б/р

Одномерный массив размерностью M.

Род краевых условий по зонам верхней границы в порядке их нумерации (см. табл. 2)

 

MK(5)

 

б/р

Одномерный массив размерностью 5.

Род краевых условий соответственно по нижней, ближней, дальней, левой и правой границам

 

JP(JW)

 

б/р

Одномерный массив размерностью JW.

Номера элементов по направлению J, задающие «разрезы», для которых выводится на печать распределение температур (см. п. 3.6). Если печать не нужна, то JP(1)=0

 

KP(KW)

 

б/р

Одномерный массив размерностью KW.

Номера элементов по направлению K, задающие «разрезы», для которых выводится на печать распределение температур (см. п. 3.6). Если печать не нужна, то KP(1)=0

 

IP(IW)

 

б/р

Одномерный массив размерностью IW.

Номера элементов по направлению I, задающие «разрезы», для которых выводится на печать распределение температур (см. п.3.6). Если печать не нужна, то IP(1)=0

Все переменные, соответствующие одному порядковому номеру оператора ввода, вводятся по единому формату.

4.4. Выходные данные

При работе программы за каждый шаг расчета по времени обновляется трехмерный массив температур UN(i,j,K), содержащий значения температур в центрах элементов исследуемого массива грунта. На заданный в исходных данных (п. 3.6) момент времени может быть выведено на печать двухмерное температурное поле, представляющее собой разрез массива U, определенный исходными параметрами (см. п. 3.6). Если в некотором элементе происходят фазовые переходы, то есть температура грунта равна температуре начала фазовых переходов (U*), то для этого элемента печатается число, показывающее, какая его часть находится в талом состоянии, по формату: 7ХХ.ХХ, где 7 - условная цифра, а ХХ.ХХ - процентное содержание талого грунта в элементе.

Кроме того, на печать выводятся номер элемента, через который проведен разрез, расстояние от его центра до начала координат и время на моменты печати. Для идентификации точек температурного поля двухмерная матрица окаймляется номерами элементов и значениями расстояния от их центров до начала координат.

Пример печати выходных данных приведен на рис. 5.

 

 

Рис. 5. Пример печати выходных данных

4.5. Подготовка данных для расчета по программе PROGNOZ-2S

Исходные данные готовятся согласно пп. 4.2-4.3 настоящего раздела, но следующие переменные имеют постоянные значения:

NK=1

iW=1

jW=1

KW=1

HK(1)=1

Ввод массивов jP, KP, iP отсутствует.

4.6. Подготовка данных для расчета по программе PROGNOZ-3N

Вся исходная информация вводится с перфокарт. Данные готовятся согласно табл. 9. Все переменные, соответствующие одному порядковому номеру оператора ввода, вводятся по единому формату. Максимальное число расчетных элементов 3000, то есть (Ni´Nj´NK £ 3000).

Таблица 9

Порядковый номер оператора ввода

Формат

Идентификатор

Единица измерения

Значение

1

3I5

Ni

б/р

Число элементов по направлению i

Nj

б/р

Число элементов по направлению j

NK

б/р

Число элементов по направлению K

2

16I5

NSL

б/р

Число выделенных литологических разностей (NSL£10)

 

NT

б/р

Число изменений шага по времени (NT£10)

 

Li

б/р

Число интервалов времени в периоде, во время которых задаются средние значения характеристик внешней среды на верхней границе (Li£36)

 

M

б/р

Число зон с различными краевыми условиями на верхней границе (М£9)

 

M1

б/р

Число зон с различными краевыми условиями на нижней границе (М£9)

 

iW

б/р

Число выдаваемых на печать «разрезов» по направлению i (см. п.3.6)

 

jW

б/р

Число выдаваемых на печать «разрезов» по направлению j (см. п.3.6)

 

KW

б/р

Число выдаваемых на печать «разрезов» по направлению К (см. п.3.6)

 

iNT

б/р

Размер массива (см. п.3.6)

 

NSt

б/р

=1, если количество теплоты измеряется в кДж

¹1, если количество теплоты измеряется в ккал

 

 

NGR

б/р

Число элементов, задающих «ступеньки» на верхней границе области исследования. Примечание. Присутствует только в алгоритме PROGNOZ-L

3

12F6.2

T(3,NSL)

°C

По три значения температуры, для которых задается количество незамерзшей воды, для каждой грунтовой разности в порядке их нумерации (см. табл. 4)

4

12F6.2

W(3,NSL)

доли ед.

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

5

13F6.2

Hi(Ni)

м

Одномерный массив размерностью Ni.

Размер элементов по направлению i

6

13F6.2

Hj(Nj)

м

Одномерный массив размерностью Nj.

Размер элементов по направлению j

7

13F6.2

HK(NK)

м

Одномерный массив размерностью NK.

Размер элементов по направлению K

8

13F6.2

U(Ni,Nj,NK)

°С

Трехмерный массив размерностью Ni´ ´Nj´NK.

Начальная температура грунта по элементам (см. табл. 1).

Каждый столбец таблицы набивать на новой перфокарте

9

16I5

N(Ni,Nj,NK)

б/р

Трехмерный массив размерностью Ni´ ´Nj´NK.

Номера слоев грунта в соответствующих элементах (см. табл. 5).

Каждый столбец таблицы набивать на новой перфокарте

10

12F6.2

TB(Li,M)

°С

Двухмерный массив размерностью Li´ ´M.

Температура внешней среды на верхней границе области исследования при периодическом ее изменении по времени и зонам краевых условий в порядке их нумерации (см. табл. 2)

 

 

RB(Li,M)

м2·°С/Вт,

м2·ч·°С/ккал

Двухмерный массив размерностью Li´ ´M.

Термическое сопротивление поверхности на верхней границе при периодическом ее изменении по времени и зонам краевых условий в порядке их нумерации (см. табл. 2)

 

 

AB(Li,M)

Вт/(м2·°С),

ккал/(м2´

´ч·°С)

Двухмерный массив размерностью Li´ ´M.

Коэффициент конвективного теплообмена на верхней границе при периодическом изменении по времени и зонам краевых условий в порядке их нумерации (см. табл. 2)

11

16I5

KL(Nj,NK)

б/р

Двухмерный массив размерностью Nj´ ´NK.

Массив номеров зон краевых условий по элементам верхней и нижней границ (см. табл. 3)

Примечание. В алгоритме PROGNOZ-L в массиве KL номер зоны краевых условий по элементам верхней границы задается равным 1.

11a

615

S(NGR,6)

б/р

Двухмерный массив размерностью NGR´6.

Задает номер зоны краевых условий по верхней границе по граням граничных элементов образующих «ступеньку» (см. п. 4.8).

Примечание. Массив S присутствует только в алгоритме PROGNOZ-L

12

F6.1,

5F6.2

G(NSL,6)

 

Значение теплофизических характеристик грунта по выделенным слоям. Двухмерный массив размерностью NSL´6 (см. табл. 4).

G(i,1)

кг/м3

Плотность сухого грунта i-го слоя

G(i,2)

кДж/(кг·°С)

Удельная теплоемкость сухого грунта i-го слоя

G(i,3)

доли

Суммарная влажность грунта i-го слоя в долях к весу абсолютно сухого грунта

G(i,4)

Вт/(м·°С)

Теплопроводность мерзлого грунта i-го слоя

G(i,5)

Вт/(м·°С),

ккал/(м2´

´ч·°С)

Теплопроводность талого грунта i-го слоя

G(i,6)

°С

Температура начала фазовых переходов i-го слоя

13

6F12.2

SHT(t)

ч

Размер шага по времени с границей времени его изменения 1£ t £NT

TR(i)

ч

TKON

ч

Время окончания счета

14

6F12.2

HPRi(iNT)

ч

Одномерный массив размерностью iNT £ 50.

Задает интервалы выдачи на печать температуры грунта (час) (см. п. 3.6)

15

13F6.2

TBiK(M1)

°С

Одномерный массив размерностью MI.

Температура внешней среды по зонам нижней границы в порядке их нумерации

TBj

Вт/(м2·°С),

ккал/(м2·ч)

Простая переменная. Величина теплопотока на левой боковой границе области исследования

TBjK

Вт/(м2·°С),

ккал/(м2·ч)

Простая переменная. Величина теплопотока на правой боковой границе области исследования

TBK

Вт/(м2·°С),

ккал/(м2·ч)

Простая переменная. Величина теплопотока на ближней боковой границе области исследования

TBKK

Вт/(м2·°С),

ккал/(м2·ч)

Простая переменная. Величина теплопотока на дальней боковой границе области исследования

16

16I5

MR(M)

б/р

Одномерный массив размерностью М.

Род краевых условий по зонам верхней границы в порядке их нумерации (см. табл. 2)

17

16I5

MK5

б/р

Одномерный массив размерностью 5.

Род краевых условий соответственно на нижней , ближней, дальней, левой и правой границах

18

16I5

jP(jW)

б/р

Одномерный массив размерностью jW.

Номера элементов по направлению j, задающие «разрезы», для которых выводится на печать температурное поле (см. п. 3.6).

Если печать не нужна то jP(1)=0

19

16I5

KP(KW)

б/р

Одномерный массив размерностью КW.

Номера элементов по направлению К, задающие «разрезы», для которых выводится на печать температурное поле (см. п. 3.6).

Если печать не нужна то КP(1)=0

20

16I5

iP(iW)

б/р

Одномерный массив размерностью iW.

Номера элементов по направлению i, задающие «разрезы», для которых выводится на печать температурное поле (см. п. 3.6).

Если печать не нужна то iP(1)=0

4.7. Подготовка данных для расчета по программе PROGNOZ-2N

Исходные данные готовятся по табл. 9. Следующие переменные имеют постоянные значения:

NK=1

iW=1

jW=1

KW=1

HK(1)=1

TBK=0

TBKK=0

Операторы ввода с порядковыми номерами 18, 19, 20 отсутствуют. Максимальное число расчетных элементов 3000, т.е. (Nj´NK £ 3000).

4.8. Подготовка данных для расчета по программе PROGNOZ-L

Данные готовятся согласно п. 4.7 настоящего раздела. Значение начальной температуры в фиктивных элементах (элементах, дополняющих область исследования до прямоугольной) задаются равными условному числу 99.99.

Краевые условия на верхней ступенчатой границе определяются номером зоны краевых условий на левой, верхней и правой гранях граничных элементов, формирующих верхнюю границу области исследования, и задаются массивом S размерностью (NGR´6) где NGR - число граничных элементов (рис. 6).

0

1

2

3

4

5

6

7

8

9

j

1

 

 

 

 

5

6

 

 

15

 

2

 

 

3

4

 

7

 

 

14

 

3

1

2

 

 

 

8

 

 

13

 

4

 

 

 

 

 

9

 

 

12

 

5

 

 

 

 

 

 

10

11

 

 

6

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

i

 

 

 

 

 

 

 

 

 

Рис. 6. Пример задания ступенчатой области исследования

 

- фиктивные элементы

 

4

- граничные элементы, 4 - порядковый номер элемента

S (d,1) - порядковый номер элемента, S (d,1)=d;

S (d,2) - координата элемента по направлению i 1£ S (d,2) £ Ni;

S (d,3) - координата элемента по направлению j 1£ S (d,3) £ Nj;

S (d,4) - номер зоны краевых условий на левой грани элемента, 0£ S (d,4)£ £ 9. S (d,4)=0, если эта грань не является поверхностью теплообмена с внешней средой;

S (d,5) - номер зоны краевых условий на верхней грани элемента, 0£ S (d,5) £ 9. S (d,5)=0, если эта грань не является поверхностью теплообмена с внешней средой;

S (d,6) - номер зоны краевых условий на правой грани элемента, 0£ S (d,6)£ £ 9. S (d,6)=0, если эта грань не является поверхностью теплообмена с внешней средой;

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

Ввод двухмерного массива S осуществляется по строкам оператором ввода 11-а (см. табл. 9).

Например, на верхней границе ступенчатой области (см. рис. 6) выделено 3 зоны с различными условиями теплообмена:

I зона - верхние грани 1, 2, 3, 10, 11 элементов;

II зона - верхние грани 5, 6, 15 элементов;

III зона - левые грани 3, 5, 12, 13, 14, 15 и правые грани 6, 7, 9 элементов.

Тогда массив S имеет вид:

S (d,1)

S (d,2)

S (d,3)

S (d,4)

S (d,5)

S (d,6)

1

3

1

0

1

0

2

3

2

0

1

0

3

2

3

3

1

0

4

2

4

0

1

0

5

1

5

3

2

0

6

1

6

0

2

3

7

2

6

0

0

3

8

3

6

0

0

3

9

4

6

0

0

3

10

5

7

0

1

0

11

5

8

0

1

0

12

4

9

3

0

0

13

3

9

3

0

0

14

2

9

3

0

0

15

1

9

3

2

0

Приложение 1

Рекомендуемое

Пример расчета

Требуется определить температурный режим, сформировавшийся в грунте в контуре отапливаемого производственного здания размером 20´ ´12 м без проветриваемого подполья через 10 лет после начала эксплуатации.

В процессе эксплуатации здания с одной из его сторон на части прилегающей территории в зимнее время будут образовываться снежные надувы (увеличение максимальной мощности снега на 0,4 м по сравнению с естественным).

Рассматривался режим теплообмена на прилегающей территории, наименее благоприятный для эксплуатации здания. Предполагалось, что снежные надувы образовывались симметрично со всех сторон здания. В результате определялся температурный режим под 1/4 частью здания (угол 6´ ´10 метров), образованной осями симметрии.

Область исследования представляет собой прямоугольный массив грунта, ограниченный глубиной 17 м и прямоугольной областью 21´15 м на верхней и нижней границах.

Размер элементов в метрах по направлению осей выбран для расчета следующий:

по оси z: 2, 2, 2, 2, 3, 3, 3;

по оси x: 3, 3, 2, 2, 2, 3, 3, 3;

по оси y: 3, 3, 3, 2, 2, 2.

Всего расчетных элементов: 7´8´6=336.

Шаг по времени 121,67 часа.

На верхней границе выделено три зоны с различными краевыми условиями. На нижней границе, на глубине 17 м температура принята постоянной, равной - 1,4 °С. На боковых границах - условие теплоизоляции, теплопоток равный 0.

Расположение зон краевых условий на верхней границе приведено на рис. 1 настоящего приложения.

 

 

Рис. 1. Разбивка верхней границы области исследования
на элементы
z зоны краевых условий

Условия теплообмена по месяцам и по зонам верхней границы приведены в табл. 1 приложения.


Таблица 1

Номер зоны

Род краевых условий

Характеристика

Месяц

I

II

III

IV

V

VI

VII

VIII

IX

X

XI

XII

1

III

Температура tn =t+Dt (°C)

-24,6

-23,2

-19,2

-9,5

3,7

13,1

18,0

12,4

4,6

-4,6

-16,4

-22,4

Термическое сопротивление R, м2·ч·°С/ккал

2,6

2,7

2,8

3,0

1,7

-

-

-

-

1,3

2,1

2,3

Коэффициент теплообмена a, ккал/(м2·ч·°С)

13,8

12,3

12,3

12,0

12,1

12,6

13,4

14,0

14,3

13,4

13,5

12,7

2

III

Температура t (°C)

t=20°C= Const

Термическое сопротивление R, м2·ч·°С/ккал

R=1,8= Const

Коэффициент теплообмена a, ккал/(м2·ч·°С)

a=20.0= Const

3

III

Температура tn =t+Dt (°C)

-24,6

-23,2

-19,2

-9,5

3,7

13,1

18,0

12,4

4,6

-4,6

-16,4

-22,4

Термическое сопротивление R, м2·ч·°С/ккал

4,4

4,4

4,7

5,0

2,8

-

-

-

-

2,1

3,5

3,9

Коэффициент теплообмена a, ккал/(м2·ч·°С)

13,8

12,3

12,3

12,0

12,1

12,6

13,4

14,0

14,3

13,4

13,5

12,7


Начальное распределение температуры для всей области исследования меняется только по глубине.

Глубина, м

1

3

5

7

9,5

12,5

15,5

°С

-0,1

-0,4

-1,0

-1,2

-1,4

-1,4

-1,4

Начальные условия приведены на январь месяц.

 

 

Рис. 2. Разрез грунта с выделенными литологическими разностями

На рис. 2 приложения изображен схематический разрез грунта, в котором выделены три различные грунтовые разности, физические и теплофизические свойства которых приведены в табл. 2. Каждая грунтовая разность по всем трем направлениям содержит целое число элементов.


Таблица 2

Номер слоя

Описание грунта

Характеристика

rd, кг/м3

Cd, ккал/(м·ч·°С)

Wtot, доли

lM, ккал/(м·ч·°С)

lT, ккал/(м·ч·°С)

WW, доли

1

Суглинок буро-коричневый, ожелезненный, с включением гравия и кальки

1390

0,22

0,25

1,3

1,15

WW (-0,3)=0,14

 

 

 

 

 

WW (-1,0)=0,12

 

 

 

 

 

WW (-10,0)=0,08

2

Супесь буровато-коричневая, сильно ожелезненная с прослоями суглинка и песка

1520

0,24

0,22

1,55

1,55

WW (-0,3)=0,12

 

 

 

 

 

WW (-1,0)=0,08

 

 

 

 

 

WW (-10,0)=0,05

3

Песок пылеватый, желтоко-

коричневый, буро-серый,

ожелезненный, оторфован

1500

0,23

0,27

2,35

2,15

WW (0)=0,0


Распределение температур требуется выдавать на печать на каждый год от начала эксплуатации по разрезам: АА¢, ВВ¢, СС¢, ДД¢, ЕЕ¢, FF¢, KK¢ (см. рис. 1).

Для поставленной задачи прогноза текст программы MAIN и BLOCK DATA для расчета по программе PROGNOZ-3S, числовая информация на бланках, заполняемая по табл. 7 и 8 раздела 4 даны в соответственно в обязательном приложении 3 и в настоящем.

Пример обработки полученных результатов приведен на рис. 3 (вклейка).


 

 

Рис. 3. Результаты расчета температурного режима
грунта по программе
PROGNOZ


 

Исходные данные для расчета контрольного примера
приложения 1 по программе
PROGNOZ-3S

 

Перфорировал _______ ___                 Проверил_____________

                                                                перфорацию

 

Перфорировал _______ ___                 Проверил_____________

                                                                перфорацию

 

Перфорировал _______ ___                 Проверил_____________

                                                                перфорацию

 


Выдача на печать для контроля ввода начальных данных,
подготовленных для расчета по программе
PROGNOZ-3S
контрольного примера приложения 1

 

-0.3

-1.0

-10.0

-0.3

-1.0

-10.0

0.0

0.0

0.0

 

 

 

 

0.14

0.12

0.08

0.12

0.08

0.05

0.0

0.0

0.0

 

 

 

 

1.30

1.55

2.35

1.15

1.45

2.15

2.0

2.0

2.0

2.0

3.0

3.0

3.0

3.0

3.0

2.0

2.0

2.0

3.0

3.0

3.0

3.0

3.0

3.0

2.0

2.0

2.0

121.67

 

 

 

 

 

 

 

 

 

 

 

876000.0

 

 

 

 

 

 

 

 

 

 

 

-24.6

-23.2

-19.2

-9.5

3.7

13.1

18.0

12.4

4.6

-4.6

-16.4

-22.4

 

20.0

20.0

20.0

20.0

20.0

20.0

20.0

20.0

20.0

20.0

20.0

20.0

 

-24.6

-23.2

-19.2

-9.5

3.7

13.1

18.0

12.4

4.6

-4.6

-16.4

-22.4

 

2.6

2.7

2.8

3.0

1.7

0.0

0.0

0.0

0.0

1.3

2.1

2.3

 

1.8

1.8

1.8

1.8

1.8

1.8

1.8

1.8

1.8

1.8

1.8

1.8

 

4.4

4.4

4.7

5.0

2.82

0.0

0.0

0.0

0.0

2.14

3.57

3.59

 

13.8

12.3

12.3

12.0

12.1

12.6

13.4

13.0

14.3

13.4

13.5

12.7

 

20.0

20.0

20.0

20.0

20.0

20.0

20.0

20.0

20.0

20.0

20.0

20.0

 

13.8

12.3

12.3

12.0

12.1

12.6

13.4

13.0

14.3

13.4

13.5

12.7

 

-1.4

0.0

0.0

0.0

0.0

 

 

 

 

 

 

 

 

3

3

3

1

2

2

2

2

1

4

6

8

1

4

6

0

 

 

 

 

 

 

/*

 

 

 

 

 

 

 

 

 

 

 

 

//

 

 

 

 

 

 

 

 

 

 

 

 


Приложение 2

Рекомендуемое

Пример текста подпрограммы ISTO

 

С                ПОДПРОГРАММА ISTO

С       ОПРЕДЕЛЯЕМ КОЛИЧЕСТВО ВЫДЕЛЕННОЙ ИЛИ

С       ПОГЛОЩЕННОЙ F В КАЛЛ (ИЛИ КДЖ) В ЕДИНИЦУ

С       ВРЕМЕНИ (ЧАС) В ФИКСИРОВАННОМ ЭЛЕМЕНТЕ С

С       КООРДИНАТАМИ (I,J,K)

С       Т - ТЕКУЩЕЕ ВРЕМЯ ОТ НАЧАЛА РАСЧЕТА

SUBROUTINE ISTO (I,J,K,T,F)

IF (T.GT.152.0).AND.(T.(T.154.0)) GO TO 10

GO TO 11

10 IF ((((K.EQ.3).AND.((J.EQ.3).OR.(J/EQ.5))).OR((J.EQ.4).AND.

# (K.EQ.5))).AND.(I.LE.7)) GO TO 12

GO TO 11

12 IF (I.LE.6) F=16666.6

IF (I.EQ.7) F=25000.0

11 OCNTINUE

RETURN

END

Приложение 3

Обязательное

Текст программы PROGNOZ-3S с результатами
расчета контрольного примера приложения 1

 

Раздел MAIN-3S содержит подпрограммы MAIN и BLOCK DATA, составленные для решения контрольного примера приложения 1.

Вызов программы на счет из библиотеки исходных модулей имеет вид:

// JOB

// CLASS=C, MSGLEVEL=1, TIME=15

// EXEC FORTHCLG, PARM, FORT=(¢OPT=02¢),

// REGION, FORT=300K, REGION, LKED=100 K, REGION,

// GO=100 K

// FORT. SYSIN DD DSN=BIBL 4,S (MAIN3S),

//                          VOL=SER=DTN033, DISP=SHR, UNIT=5061

//                DD DSN=BIBL 4,S (GRUN3S),

//                          VOL=SER=DTN033, DISP=SHR, UNIT=5061

//                DD DSN=BIBL4,S (WNS),

//                          VOL=SER=DTN033, DISP=SHR, UNIT=5061

//                DD DSN=BIBL4,S (FIND),

//                          VOL=SER=DTN033, DISP=SHR, UNIT=5061

//                DD DSN=BIBL4,S (HF),

//                          VOL=SER=DTN033, DISP=SHR, UNIT=5061

Здесь : BIBL4,S - имя библиотеки исходных модулей;

BIBL4,S - имя тома.

Результаты расчета контрольного примера представлены на 1-й и 5-й расчетный год.

*** MAIN3S ***

1

 

 

DIBENSION H(7,8,6),UR(7,8,6),HRDV(7,8,6),J(7,8,6),

2

 

#

A1(3),B1(3),C1(3)

3

 

 

IAIEGER NN(7,8,6)/336*3/XL(3,6)/14*111,2*131,2*111,12

4

 

#

*131,2*114*121/131,3*111,4*121,4*111/

5

 

 

REAL HPRI(2/-1.0,8760.0/

6

 

 

CALLWNS(3,A1,B1,C1)

7

 

 

TKUR=17601.0

8

 

 

DO 1 K=1,8

9

 

 

DO 1 J=1,8

10

 

1

N(1,J,K)=1

11

 

 

DO 2 K=1,3

12

 

 

DO 2 J=3,8

13

 

 

DO 2 I=2,5

14

 

2

N(I,J,K)=2

15

 

 

DO 3 K=1,5

16

 

 

DO 3 J=1,8

17

 

 

U(1,J,K)=-0,1

18

 

 

U(2,J,K)=-0,4

19

 

 

U(3,J,K)=-1,0

20

 

 

U(4,J,K)=-1,2

21

 

 

U(5,J,K)=-1,4

22

 

 

U(6,J,K)=-1,4

23

 

3

U(7,J,K)=-1,4

24

 

 

CALL GRAND(1,7,5,6,7,12,3,1,TKDH,A1,B1,C1,H,KL,U,

25

 

#

1,4,3,2,HPRI,DV,HRDV

26

 

 

STOP

27

 

 

END

 

C

 

 

 

C

 

 

 

C

 

 

 

C

 

 

 

C

 

 

28

 

 

BLOCK DATA

29

 

 

COMMON/RBIS/BCK(10),CCK(10),WC(10),UKR(10)

30

 

 

DATA BCK/1390,0.1520,0.1500.0,7*0/. CCK/.22,0.24,0.23

31

 

#

7*0/,WC/0.25,0.27,7*0.0/,UKR/10*0.0/

32

 

 

END

*** GRAN3S ***

1

 

 

SLBROLTINE GRUND(SL,NI,NJ,NK,NT,LI,M1,IK,AX,KL

2

 

#

,IW,JW,KW,HPRI,H,HNDV,SI)

3

 

 

DIMENSION H(NI,NJ,NK),Z(100),R(100),SK(100),

4

 

1

CIAL(10),CI(10), TLN(10),TLT(10)HI(100),HJ(100)HK(100)

5

 

2

TR(10),U(NI,NJ,NK)TB(36,10)RB(36,10)AB(36,10)

6

 

3

KL(NJ,NK),NJ(10),NK(10),JP(100),IP(100),KP(100)

7

 

4

AX(A,SL),BX(B,SL),CX(C,SL),UN(NI,NJ,NK)

8

 

 

COMMO\ /BCK(10),CCK(10),WC(10),HKP(10)

9

 

 

COMMO\ /A(10),B(10),C(10)

10

 

 

COMMO\ /RV,CL,CB,TFP

11

 

 

COMMO\ /FT(50,10),F510(10,10),F25(15,10),F02(20,15),

12

 

#

H1(10),H2(10)

13

 

 

RV=1

14

 

 

CL=2.0234

15

 

 

CB=4.0368

16

 

 

TFP=337.944

 

 

 

PP=3.6

17

 

 

IF (RSI.ET.1) GO TO 1

18

 

 

CL=0.5

19

 

 

CB=1.0

20

 

 

PP=1.0

21

1

 

TK=273.0

22

 

 

TPER=3760.0

23

 

 

DO 180 I=1,HSL

24

 

 

A(I)=AX(I)

25

 

 

B(I)=BX(I)

26

180

 

C(I)=CX(I)

27

 

 

DO 216 K=1,HSL

28

 

 

DO 217 IN=1,SJ

29

 

 

UL1=-1

30

217

 

FI(IN,K)=NF(K,UL1)

31

 

 

DO 218 IN=1,10

32

 

 

UL1=-5.0-0.5*IN

33

218

 

F510(I1,K)=IF(K,UL1)

34

 

 

DO 219 IN=1,15

35

 

 

UL1=-2.0-0.2*IN

36

219

 

F25(IN,K)=HI(K,UL1)

37

 

 

DO 220 IN=1,20

38

 

 

UL1=-0.1*IN

39

220

 

FCZ(IN,K)=IN(K,UL1)

40

216

 

CONTINUE

41

181

 

FORMAT (13F5.2)

42

182

 

FORMAT (12F6.2)

43

183

 

FORMAT (6F12.2)

44

184

 

FORMAT (27IN)

45

 

 

READ 131,(TLH(I),I=1,NSL),(TLT(I),I=1,NSL),(HI(I),I=1,NI),

46

 

#

(HJ(I),II=1,NJ),(HK(I),I=1,NK),(SLT(I),I,NI)

47

 

 

READ 133,(TR(I),I=1,NI)

48

 

 

READ 132, ((IB(I,J),I=1,LI),J=1,N),((RF(I,J),I=1,LI),J=1,N),

49

 

#

((AB(I,J),I=1,LI),J=1,N),(TBIK(I),I=1,N1),TRJ,TBJK,TBK,TB

50)

 

 

READ 134,(HR(I),I=1,N),(NK(I),I=1,5),(JP(I),J=1,JP),(KP(I),

51

 

#

I=1,KP),(IP(I),I=1,IP)

52

 

 

T=0.0

53

 

 

IT=1

54

 

 

KT=1

55

 

 

T2=TPER/LI

56

 

 

HI=SHI(IT)

57

133

 

FORMAT (5X,2IK=,I3,2X,F7.2,1H,2X,I=,F2.2,12X,1H=,

58

 

1

2X,I3,7YПЕРИОДАБ,I3,9H’ИНТЕГРАЛ’F7.2,4HЧАСА)

59

134

 

FORMAT (2X,12HБЛОКИ,15,16,17)

60

135

 

FORMAT (2X,12HМЕТРЫ,13,14,15)

61

136

 

FORMAT (2X,15,17F7.2)

62

160

 

FORMAT (5A,HJ=,13,2X,17F.2

63

 

 

IH=,2X,I3,7HПЕРИОДА,I3,9H’ИНТЕГРАЛА’F7.2,4НЧАСА)

64

161

1

FORMAT (2X,H2’,БЛОК’,15,16I7)

65

162

 

FORMAT (5X,2HПЕРИОДА,I3,9HИНТЕРВАЛА,

66

 

 

F7.2,4HЧАСА)

67

163

 

FORMAT (2X,H2МЕТРЫ,13F3.2)

68

 

 

DO 5 I=1,NSL

69

 

 

H1(I)=HF(2,UKR(I))-15.0

70

 

 

H2(I)=HF(1)+TFR*(CCK(I)*BCK(I)*(WC(I)-A(I)/(B(I)-UKP

 

 

 

(I))-C(I))+15.0

71

 

 

CIAL(I)=BCK(I)*(CCK(I)+WC(I)+CB)

72

5

 

CW(I)=CCK(I)+CL*HC(I)+C(I)*(CB-RV*CL)

73

 

 

DO 5 K=1,K

74

 

 

DO 6 J=1,J

75

 

 

DO 5 I=1,I

76

 

 

DV=NI(I)*NJ(J)*NK(K)

77

 

 

N1=N(I,J,K)

78

 

 

UL=(I,J,K)

79

 

 

IF (UL.LT.(NKR(N1)-0.0001)) GO TO 7

80

 

 

IF (UL.LT.(NKR(N1)+0.0001)) GO TO 9

81

 

 

GO TO 3

82

 

 

H(I,J,K)=HF(N1,UL)

83

 

 

GO TO 3

84

 

 

H(I,J,K)=(O.5*(H1(N1)+H2(N1))

85

 

 

GO TO 3

86

 

 

H(I,J,K)=(H(N1)+CTIL(N1)*(UL-UKZ(N1)))

87

 

 

CONTINUE

88

 

 

IPR=Z

89

 

 

IPR=HPRI(Z)

90

 

 

SK(1)=HK(J)

91

 

 

IF (NK.EO.1) GO TO 175

92

 

 

DO 132 K=2,NK

93

132

 

SN(K)=SK(K)

94

175

 

DO 138 K=1,NK

95

178

 

SN(K)=SK(K)-0.5*HK(K)

96

 

 

R(1)=KJ(1)

97

 

 

IF (NJ.EO.1,GO TO 176

98

 

 

DO 139 J=2,HJ

99

130

 

R(J)=R(J-1)*HJ(J)

100

176

 

DO 140 J=1,HJ

101

140

 

R(J)=R(J)-0.5*HJ(J)

102

 

 

Z(1)=HI(1)

103

 

 

IF (NI.EO.1) GO TO 177

104

 

 

DO 141 I=2,HI

105

141

 

Z(I)=Z(I-1)*HI(I)

106

172

 

DO 142 I=1,HI

107

142

 

Z(I)=Z(I)-0.5*HI(I)

108

11

 

FORMAT (20X,17HРЕШАЕТСЯ ЗАДАЧА 3,13,11Н

 

 

 

ИЗММЕРЕНИЯХ)

109

12

 

FORMAT (15X,16HНАПРАВЛЕНИЕ ОСЕЙ,,1НК)

110

13

 

FORMAT (140,30Х,6Н0----J)

111

14

 

FORMAT (31X,1H/)

112

22

 

FORMAT (30X,1HI)

113

23

 

FORMAT (10X,22HШАГ ПО ВРЕМЕНИ В ЧАСАХ,

 

 

 

(8F11.2))

114

24

 

FORMAT (12X,17HГРАНИЦА ИЗМЕРЕНИЯ,5Х, (8F11.2))

115

25

 

FORMAT (140,9X,17HЧИСЛО БЛОКОВ ПО I,I4)

116

26

 

FORMAT (2X,12HНОМЕР БЛОКА,I4,15I7)

117

27

 

FORMAT (2X,12HРАЗМЕР БЛОКА,16F7.2)

118

28

 

FORMAT (2X,12HКООРД. ЦЕНТРА,16F7.2)

119

29

 

FORMAT (140,9X,17HЧИСЛО БЛОКОВ ПО J,I4)

120

30

 

FORMAT (140,9X,17HЧИСЛО БЛОКОВ ПО K,I4)

121

33

 

FORMAT (140,14X,21HУСЛОВИЯ НА ГРАНИЦАХ)

122

34

 

FORMAT (3X,3HK/J,I4,23I5)

123

35

 

FORMAT (2X,I4,24I5)

124

36

 

FORMAT (5X,32HШИФР XYZW ЕСЛИ X=0 - НЕТ

 

 

 

ИСТОЧНИКА)

125

37

 

FORMAT (20X,25HY - НОМЕР ТИПА ИСТОЧНИКА)

126

38

 

FORMAT (20X,49HZ - НОМЕР ЗОНЫ КРАЕВЫХ

 

 

 

УСЛОВИЙ НА ВЕРХНЕЙ ГРАНИЦЕ)

127

39

 

FORMAT (20X,18H - ТОЖЕ НА НИЖНЕЙ)

128

40

 

FORMAT (140,21X,18H - НА ВЕРХНЕЙ ГРАНИЦЕ)

129

41

 

FORMAT (2X,19HНОМЕР ЗОНЫ/МЕСЯЦ,I4,12I7)

130

43

 

FORMAT (2X,16,4X,5HТЕМП ВОЗД.,15F7.2)

131

46

 

FORMAT (12X,9HСОПРОТИВЛ.,15F7.2)

132

47

 

FORMAT (12X,9HТЕПЛООТД.,15F7.2)

133

48

 

FORMAT (2X,24HПРОДОЛЖИТЕЛЬНОСТЬ МЕСЯЦА,

 

 

 

F9.2,5H ЧАСА)

134

49

 

FORMAT (2X,24HПЕРИОДИЧНОСТЬ ИЗМЕРЕНИЯ,

 

 

 

F9.2,5H ЧАСА)

135

121

 

FORMAT (140,14X,29HТЕМПЕРАТУРА НА ВЕРХНЕЙ

 

 

 

ГРАНИЦЕ)

136

122

 

FORMAT (5X,10HНОМЕР ЗОНЫ,14,15I7)

137

123

 

FORMAT (5X,10HТЕМПЕРАТУРА,15F7.2)

138

124

 

FORMAT (15X,10H)

139

125

 

FORMAT (5X,10HTBJ=,F7.2,2X,5HTBJ<=,F7.2)

140

 

1

1F7.2

141

126

 

FORMAT (10X,11HНОМЕР СЛОЯ,2X)

142

 

*

25HТЕПЛОПРОВ,3X)

143

127

 

FORMAT (5X,15,F12.2,5F7.2,3F7.3)

144

128

 

FORMAT (140,15X,28HНАЧАЛЬНАЯ ТЕМПЕРАТУРА

 

 

 

ГРУНТА)

145

130

 

FORMAT (15X,16I7)

146

131

 

FORMAT (5X,25HШАГ ГЕНЕРАТОРА РЕЗУЛЬТАТА

 

 

 

(10F9.2))

147

132

 

FORMAT (5X,2HK=,13,2X,F7.2)

148

164

 

FORMAT (5X,2HJ=,13,2X,F7.2)

149

110

 

FORMAT (5X,2HI=,13,2X,F7.2)

150

 

 

L=2

151

 

 

IF ((NK.ST.1,ND.(NJ.ST.1)) GO TO 150

152

 

 

IF ((NK.EO.1,ND.(NJ.EO.1)) L=1

153

 

 

DO TO 151

154

150

 

L=3

155

151

 

PRINT 11,L

156

 

 

PRINT 12

157

 

 

PRINT 13

158

 

 

DO 152 I=I,4

159

152

 

PRINT 14

160

 

 

PRINT 22

161

 

 

PRINT 23,(SHT(I),I=1,NI)

162

 

 

PRINT 24,(TR(I),I=1,NT)

163

 

 

PRINT 25,I

164

 

 

DO 153 L=1,NI,15

165

 

 

LK=L+14

166

 

 

IF (LK.ST.(I)LK=NI)

167

 

 

PRINT 26,(I,I=L,LK)

168

 

 

PRINT 27,(I,I=L,LK)

169

153

 

PRINT 28,(Z(I),I=L,LK)

170

 

 

PRINT 29,J

171

 

 

DO 154 L=1,HJ,15

172

 

 

LK=L+14

173

 

 

IF (LK.ST.J)LK=NJ)

174

 

 

PRINT 26,(J,J=L,LK)

175

 

 

PRINT 27,(J,J=L,LK)

176

154

 

PRINT 28,(R(J),J=L,LK)

177

 

 

PRINT 30,K

178

 

 

DO 155 L=1,NK,15

179

 

 

LK=L+14

180

 

 

IF (LK.ST.NK) LK=NK

181

 

 

PRINT 26,(K,K=L,LK)

182

 

 

PRINT 27,(NK(K),K=L,LK)

183

155

 

PRINT 28,(SK(K),K=L,LK)

184

 

 

PRINT 33

185

 

 

DO 157 L=1,NJ,24

186

 

 

LK=L+23

187

 

 

IF (LK.GT.NJ) LK=NJ

188

 

 

PRINT 34,(J,J=L,LK)

189

 

 

DO 157 K=1,NK

190

 

 

I=NK=K+1

191

157

 

PRINT 35,(I,(KL(J,I),J=L,LK)

192

 

 

PRINT 36

193

 

 

PRINT 37

194

 

 

PRINT 38

195

 

 

PRINT 39

196

 

 

PRINT 40

197

 

 

DO 154 K=1,NK

198

 

 

PRINT 41,(I,I=1,LI)

199

 

 

PRINT 43,(I,TB(I,K),I=1,LI)

200

 

 

PRINT 46,(I,RB(I,K),I=1,LI)

201

158

 

PRINT 47,(I,AB(I,K),I=1,LI)

202

 

 

PRINT 48,TZ

203

 

 

PRINT 49,TPER

204

 

 

PRINT 121

205

 

 

PRINT 122,(I,I=1,MI)

206

 

 

PRINT 123,(TBIK(I),I=1,MI)

207

 

 

PRINT 124

208

 

 

PRINT 125,TBJ,TBJK,TBK,TBKK

209

 

 

PRINT 126

210

 

 

DO 159 K=1,NSL

211

 

 

PRINT 127,K,BCK(K),CCK(K),WC(K),UKR(K),TLM(K),

212

 

2

TLT(K),l(K),B(K),C(K)

213

159

 

CONTINUE

214

 

 

PRINT 129

215

 

 

IF (KP(1).EO,J) GO TO 160

216

 

 

DO 170 LP=1,KP

217

 

 

K=KP(JP)

218

 

 

PRINT 132,K,SK(K)

219

 

 

DO 170 L=1,HJ,16

220

 

 

LK=L+15

221

 

 

IF (LK.GT.NJ) LK=NJ

222

 

 

PRINT 134,(J,J=L,LK)

223

 

 

PRINT 135,(R(J),J=L,LK)

224

 

 

DO 170 I=1,NI

225

 

 

PRINT 136,(I,Z(I),(U(I,J,K),J=L,LK))

226

170

 

PRINT 130,(H(I,J,K),J=L,LK)

227

160

 

IF (JP(1).EO.0) GO TO 171

228

 

 

DO 172 LP=1,JN

229

 

 

J=JN(LP)

230

 

 

PRINT 164,J,R(J)

231

 

 

DO 172 L=1,NK,16

232

 

 

LK=L+15

233

 

 

IF (LK.GT.NK) LK=NK

234

 

 

PRINT 161,(K,K=L,LK)

235

 

 

PRINT 135,(SK(I),I=L,LK)

236

 

 

D0 172 I=1,NI

237

 

 

PRINT 136,(I,Z(I),(U(I,J,K),K=L,LK))

238

172

 

PRINT 130,(H(I,J,K),K=L,LK)

239

171

 

IF (IP(1).EO.O) GO TO 173

240

 

 

D0 174 LP=1,IP

241

 

 

I=IP(LP)

242

 

 

PRINT 110,I,Z(I)

243

 

 

D0 174 L=1,NK,16

244

 

 

LK=L+15

245

 

 

IF (LK.GT.NK) LK=NK

246

 

 

PRINT 161,(K,K=L,LK)

247

 

 

PRINT 163,(SK(J),J=L,LK)

248

 

 

DO 174 J=1,NJ

249

 

 

PRINT 136,(J,R(J),(U(I,J,K),K=L,LK))

250

174

 

PRINT 130,(H(I,J,K),K=L,LK)

251

173

 

PRINT 131,(PRHI(I),I=2,INT)

252

10

 

T=T+NT

253

201

 

CONTINUE

254

 

 

IF (T.GT.TKON) GOTO 15

255

 

 

IF ((1-KT*TPC).GT.0.0) KT=KT+1

256

 

 

IF (T-TR(IT) 16,17,17

257

17

 

IT=IT+1

258

 

 

HI=SHI(IT)

259

16

 

T3=T-(KT-1).TPER

260

 

 

IBR=1

261

18

 

IF (T3-IBR*T2) 19,19,20

262

20

 

IBR=IBR+1

263

 

 

GO TO 13

264

19

 

T4=T3-(IBR-1)*T2

265

21

 

D0 31 K=1,NK

266

 

 

D0 32 J=1,NJ

267

 

 

IK=KL(J,K)

268

 

 

MC=IK/10-IK/100*10

269

 

 

MCK=IK/10*10

270

 

 

DO 50 I=1,NI

271

 

 

H1=,(I,J,K)

272

 

 

HN=H(I,J,K)

273

 

 

UL=U(I,J,K)

274

 

 

IF (HN-H1(NI) 51,52,52

275

51

 

TL=TLX(N1)

276

 

 

GO TO 53

277

52

 

IF (HN-H2(NI) 45,45,44

278

45

 

TL=0.5*(TLX(N1)+TLT(N1)

279

 

 

GO TO 53

280

44

 

TL=TLT(N1)

281

53

 

IF (I.GE.0) GO TO 54

282

55

 

GO TO 56

283

 

 

TE10=IB(IBR-1,MC)

284

 

 

GO TO 57

285

56

 

TE10=IB(LI,MC)

286

57

 

U1=TB10+(TB(IBR,MC)-TB10)/T2*T4

287

 

 

IJ=BR(10)

288

 

 

GO TO (58,59,60),IJ

289

58

 

R1=TL/TI(1)*(J1-CL)*2

290

 

 

GO TO 66

291

59

 

R1=U1

292

 

 

GO TO 66

293

60

 

GO TO (61),IBR

294

 

 

RB10=RB(IBR-1,MC)

295

 

 

AB10=AB(IBR-1,MC)

296

 

 

GO TO 62

297

61

 

RB10=RB(LI,MD)

298

 

 

AB10=AB(LI,MD)

299

62

 

RB11=RB10+(RB(IBR,MD)-RB10)/T2*T4

300

 

 

AB11=AB10+(AB(IBR,MD)-AB10)/T2*T4

301

 

 

R1=2/(TI(1)/TL+2*(1/AB11+RB11))*(H1-UL)

302

 

 

GO TO 66

303

54

 

U1=U(I-1,J,K)

304

 

 

N2=H(I-1,J,K)

305

 

 

HN2=H(I-1,J,K)

306

 

 

IF (HN2-H1(N2)) 63,64,64

307

63

 

TL1=TLT(N2)

308

 

 

GO TO 65

309

64

 

TL1=TLT(N2)

310

 

 

IF (HN2.LS.H2(N2)) TL1=0.5*(TLN(N2)+TL1)

311

65

 

UX=HI(I-1)+HI(I)

312

 

 

R1=2*(HI(I-1)*TL1+HI(I)*TL)/)UX+LX)*(H1-UL)

313

66

 

IF (I-I1) 67,68,67

314

67

 

U2=U(I+1,J,K)

315

 

 

H2=H(I+1,J,K)

316

 

 

HN2=H(I+1,J,K)

317

 

 

IF (HN3-H1(N2)) 69,70,70

318

69

 

TLZ=TLH(N2)

319

 

 

GO TO 71

320

70

 

TLZ=TLT(N2)

321

 

 

IF (HNR.LE.12(N2)) TLZ=0.5*(TLN(N2)+TL2

322

71

 

UX=HI(I+1)+HI(I)

323

 

 

R2=2*(HI(I+1)*TLZ+HI(I)*TL)/(UX*LX)*(U2-UL)

324

 

 

GO TO 72

325

68

 

U2=UB.K(M,X)

326

 

 

IF (MX(1)-1) 73,74,73

327

74

 

R2=TL/HI(I)*(U2-UL)*2

328

 

 

GO TO 72

329

73

 

R2=U2

330

72

 

IF (J.UE.1) GO TO 75

331

 

 

U1=TBJ

332

 

 

IF (MK(2)-1) 76,77,76

333

77

 

R3=TL/HJ(1)-(U3-UL)*2

334

 

 

GO TO 70

335

76

 

R3=TBJ

336

 

 

GO TO 70

337

75

 

UJ=U(I,J-1,K)

338

 

 

N2=N(I,J-1,K)

339

 

 

HN2=H(I,J-1,K)

340

 

 

IF (HN2-H1(N2)) 80,81,81

341

80

 

TL3=TLT(N2)

342

 

 

GO TO 82

343

81

 

TL3=TLT(N2)

344

 

 

IF (HN2.LE.12(N2)) TL3=0.5*(TLM(N2)+TL3)

345

82

 

UX=HJ(J-1)+HJ(J)

346

 

 

R3=2*(HJ(J-1)*TL3+HJ(J)*TL)/(UX-LX)*(U3-UL)

347

79

 

IF (J.LE.HJ) GO TO 83

348

 

 

U4=TBJK

349

 

 

IF (MX(3)-1) 84,85,84

350

85

 

R4=TL/HJ(NJ)*(U4-UL)*2

351

 

 

GO TO 90

352

84

 

R4=TBJK

353

 

 

GO TO 90

354

83

 

U4=U(I,J+1,K)

355

 

 

N2=N(I,J+1,K)

356

 

 

HN2=H(I,J+1,K)

357

 

 

IF (H/2-H1(N2)) 87,88,88

358

87

 

TL4=TLT(N2)

359

 

 

GO TO 89

360

88

 

TL4=TLT(N2)

361

 

 

IF (HN2.LE.H2(N2)) TL4=0.5*(TLM(N2)+TL4)

362

89

 

UX=HJ(J+1)+HJ(J)

363

 

 

R4=2*(HJ(J+1)*TL4+HJ(J)*TL)/(UX*UX)*(U4-UL)

364

90

 

IFF (K.NE.1) GO TO 91

365

 

 

U5=TBK

366

 

 

IF (MX(4)-1) 92,93,92

367

93

 

R5=TL/HK(1)*(U5-UL)*2

368

 

 

GO TO 98

369

92

 

R5=TBK

370

 

 

GO TO 98

371

91

 

U5=U(I,J,K-1)

372

 

 

N2=N(I,J,K-1)

373

 

 

HN2=H(I,J,K-1)

374

 

 

IF (HN2-H1(N2)) 95,96,96

375

95

 

TL5=TLT(N2)

376

 

 

GO TO 97

377

96

 

TL5=TLT(N2)

378

 

 

IF (HN2.LE.H2(N2)) TL5=0.5*(TLM(N2)+TL5)

379

97

 

UX=HK(K-1)+HK(K)

380

 

 

R5=2*(HK(K-1)*TL5+HK(K)*TL)/(UX*UX)*(U5-UL)

381

98

 

IF (K.NE.NK) GO TO 99

382

 

 

U6=TBKK

383

 

 

IF (MK(5)-1) 100,101,100

384

101

 

R6=TL/MK(NK)*(U6-UL)*2

385

 

 

GO TO 103

386

100

 

R6=TBKK

387

 

 

GO TO 103

388

99

 

U6=U(I,J,K+1)

389

 

 

N2=N(I,J,K-1)

390

 

 

HK2=H(I,J,K+1)

391

 

 

IF (HK2-H(N2)) 104,105,105

392

104

 

TL6=TLT(N2)

393

 

 

GO TO 106

394

105

 

TL6=TLT(N2)

395

 

 

IF (HN2.LE.H2(N2)) TL6=0.5*(TLM(N2)+TL6)

396

106

 

UX=HK(K+1)*HK(K)

397

 

 

R6=2*(HK(K+1)*TL6+HK(K)*TL)/(UX*UX)*(U6-UL)

398

103

 

CONTINUE

399

 

 

F=0.0

400Р

108

 

HNDV(I,J,K)=HJ+HT*((R1+R2)/HI(I)+(R3+R4)/HJ(J)+

 

 

 

(R5+R6)/HK(K)+1)*NP

402

 

 

HK1=HNDV(I,J,K)

403

 

 

IF (HK1-H1(H1)) 111,112,112

404

112

 

IF (HK1-H2(H2)) 113,113,114

405

113

 

UL1=UKR(N1)

406

 

 

GO TO 115

407

114

 

UL1=(HN1-H2(N1))/CTAL(N1)+UKR(N1)

408

 

 

GO TO 115

409

111

 

IF ((A(N1).EO.0.0).AND.(C(N1).EO.0.0)) GO TO 4

410

 

 

CALL FIND(HN1,N1,UL1)

411

 

 

GO TO 115

412

4

 

UL1=-HN1/H1(N1)*TK+TK

413

115

 

UN(I,J,K)=UL1

414

50

 

CONTINUE

415

32

 

CONTINUE

416

31

 

CONTINUE

417

 

 

DO 42 K=1,NK

418

 

 

DO 42 J=1,NJ

419

 

 

DO 42 I=1,NI

420

 

 

U(I,J,K)=UN(I,J,K)

421

 

 

H(I,J,K)=HNDV(I,J,K)

422

42

 

CONTINUE

423

 

 

IF (T.LT.(TPR-0.01)) GO TO 10

424

 

 

IF (HPRI(1).GT.0.0) IPR=IPR+1

425

 

 

TPR=TPR+HPRI(IPR)

426

3

 

KTP=KT-1

427

 

 

IBRP=IBR-1

428

 

 

IF (KP(1).EO.0) GO TO 143

429

 

 

DO 143 LP=1,KW

430

 

 

K=KP(LP)

431

 

 

DO 145 I=1,NI

432

 

 

DO 145 J-1,NJ

433

 

 

N1=N(I,J,K)

434

 

 

HN1=H(I,J,K)

435

 

 

IF ((H1(N1).LE.HN1).AND.(HN1.LE.H2(H1))

436

 

1

UK(I,J,K)=700.0+100*(HN1-H1(N1))/(H2(N1)-H1(N1))

437

145

 

CONTINUE

438

 

 

PRINT 133,K,SK(K),T,KTP,I8RP,T4

439

 

 

DO 143 L=1,NJ,16

440

 

 

LK=L*15

441

 

 

IF (LK.GT.NJ) LK=NJ

442

 

 

PRINT 134,(J,J=L,LK)

443

 

 

PRINT 135,(R(J),J=L,LK)

444

 

 

DO 143 I=1,NI

445

 

 

PRINT 136,(I,Z(I),(UN(I,J,K),J=L,LK))

446

143

 

CONTINUE

447

 

 

IF (JP(1).EO.0) GO TO 165

448

 

 

DO 165 LP=1,JN

449

 

 

J=JP(LP)

450

 

 

DO 148 K=1,NK

451

 

 

DO 148 I=1,NI

452

 

 

N1=N(I,J,K)

453

 

 

HN1=H(I,J,K)

454

 

 

IF ((H1(N1).LE.HN1).AND.(HN1.LE.H2(N1)))

455

 

1

UN(I,J,K)=700.0+100*(HN1-H1(N1))/(H2(N1)-H1(N1))

456

148

 

CONTINUE

457

 

 

PRINT 160,J,R(J),T,KTP,I8RP,T4

458

 

 

DO 165 L=1,NK,16

459

 

 

LK=L+15

460

 

 

IF (LK.GT.NK) LK=NK

461

 

 

PRINT 161,(K,K=L,LK)

462

 

 

PRINT 135,(SK(I),I=L,LK)

463

 

 

DO 166 I=1,N1

464

166

 

PRINT 136,(I,Z(I),(UN(I,J,K),K=L,LK))

465

165

 

CONTINUE

466

 

 

IF (IF(1).EO.0) GO TO 167

467

 

 

DO 167 LP=1,IW

468

 

 

I=IP(LP)

469

 

 

DO 149 K=1,NK

470

 

 

DO 149 J=1,NJ

471

 

 

N1=N(I,J,K)

472

 

 

HN1=H(I,J,K)

473

 

 

IF ((H1(N1).LE.HN1).AND.(HN1.LE.H2(N1)))

474

 

1

UN(I,J,K)=700.0+100*(HN1-H1(N1))/(H2(N1)-H1(N1))

475

149

 

CONTINUE

476

 

 

PRINT 162,I,Z(I),T,KTP,I8RP,T4

477

 

 

DO 167 L=1,NK,16

478

 

 

LK=L+15

479

 

 

IF (LK.GT.NK) LK=NK

480

 

 

PRINT 161,(K,K=L,LK)

481

 

 

PRINT 163,(SK(J),J=L,LK)

482

 

 

DO 168 J=1,NJ

483

168

 

PRINT 136,(J,Z(J),(UN(I,J,K) K=L,LK))

484

167

 

CONTINUE

485

 

 

GO TO 10

486

15

 

CONTINUE

487

 

 

RETURN

488

 

 

END

*** WNS ***

1

 

 

SUBROLTINE WNS(N,A,B,C)

2

 

 

DIMENSION A(N),B(N),C(N),T(3,10),V(3,10),Z(10,11),U(11)

3

11

 

FORMAT (12F6,2)

4

 

 

READ 11,((T(I,J),I=1,3),J=1,’)

5

 

 

READ 11,((T(I,J),I=1,3),J=1,’)

7

8

 

FORMAT (22l,7l¢,(11)=,F6.2,3X,ZHUR(T2)=,F5.2,3X,

 

 

 

ZHUR(T3)=,F6.2)

8

 

 

PRINT 16

13

 

 

U(1)=-30.0

14

 

 

DO 6 J=2,6

15

6

 

U(J)=U(J-1)+5.0

16

 

 

DO 5 J=7,11

17

5

 

U(J)=U(J-1)+1.0

18

 

 

DO 10 I=1,4

19

13

 

W2=W(2,I)-W(3,I)

20

 

 

U1=T(1,I)-T(3,I)

21

 

 

W1=W(1,I)-W(3,I)

22

 

 

U2=T(2,I)-T(3,I)

23

 

 

IF (((-0.0001.LT.W2).AND.(W2.LT.0.0001)).OR.((-0.0001.LT.

24

 

#

W1).AND(W1.LT.0.0001))) GO TO 2

25

 

 

GO TO 24

26

2

 

A(I)=0.0

27

 

 

B(I)=0.0001

28

 

 

C(I)=W(1,I)

29

 

 

GO TO 4

30

74

 

B(I)=(T(2,I*W2+U1-T(1,I)+W1*U2)/(W2*U1-W1*U2)

31

 

 

A(I)=-1*(B(I)-T(1,I))*(B(I)-T(((3,I))/U1

32

 

 

C(I)=-(3,I)-A(I)/B(I)-T(3,I)

41

4

 

DO 3 J=1,11

42

3

 

Z(I,J)=A(I)/B(I)-U(J))+C(I)

43

10

 

CONTINUE

44

16

 

FORMAT (T41,’НЕЗАМЕРЗШАЯ ВОДА ПО СЛОЯМ,)

45*

17

 

FORMAT (2H0,T25,’ВХОДНЫЕ ДАННЫЕ’,30X,

46

 

#

‘ВЫЧИСЛЕННЫЕ КОЭФФИЦИНТЫ <PHBO’)

47

18

 

FORMAT (1H),11X,4HСЛОЙ,I3,6X,4HT1=.F6.2,6X,4HT2=

48

 

#

,F6.2X,4HT3,2HA(,12,2H)=,F6.4,4X,2HB(,IZ,ZI)=,F8.4)

49

19

 

FORMAT (22X,7HW(T1)=,F6.2,3X,7HW(T2)=,F62,3X,7HW

50

 

#

(T3)=,F6.2,3X,2HA(,I2,2H)=,F8.4,4X,2HB(,I2,2H)=,F8.4)

51

 

 

PRINT 17

52

 

 

DO 20 I=1,4

53

 

 

PRINT 18,(I,T(1,I),T(2,I),T(3,I))

54

20

 

PRINT 19,(W(1,I),W(2,I),W(3,I),I,A(I),I,B(I),I,C(I))

55

21

 

FORMAT (2H0,T40,’КОНТРОЛЬНАЯ ТАБЛИЦА ЗНАЧЕНИЙ WN ОТ Т’)

56

22

 

FORMAT (1H0,14X,2HT=,7X,11F7,2)

57

23

 

FORMAT ((11X,4НСЛОЙ,I3,6X,11F7.2)

58

 

 

PRINT 21

59

 

 

PRINT 22,U

60

 

 

PRINT 23,((I,(Z(I,J),J=1,11)),I=1,N)

61

 

 

RETURN

62

 

 

END

*** FIND ***

1

 

 

SOHROLTINL FIND(HN1,N1,UL1)

11

 

 

COMMO\ /TABL/ FT(50,10),FS10(1),F25(15,10),F02(20,

21

 

#

10),H1(10),H2(10)

31

 

 

IF (HN1.GT/IT(30,N1) GO TO 201

41

 

 

IF (HN1.GT/IT(40,N1) GO TO 202

51

 

 

ML=41

61

 

 

ML=50

71

 

 

GO TO 203

81

202

 

ML=31

91

 

 

ML=40

101

203

 

DO 204 IN=ML,MU

111

 

 

IF (HN1.GT/IT(IN,N1) GO TO 205

121

204

 

CONTINUE

131

205

 

WA=FT(IN,N1)

141

 

 

WB=FT(IN-1,N1)

151

 

 

XA=IN

161

 

 

XB=IN-1

171

 

 

GO TO 200

181

201

 

IF (HN1.GT/IT(10,N1) GO TO 206

191

 

 

IF (HN1.GT/IT(20,N1) GO TO 207

201

 

 

ML=21

211

 

 

ML=30

221

 

 

GO TO 203

231

207

 

ML=11

241

 

 

ML=20

251

 

 

GO TO 203

261

208

 

IF (HN1.GT/IT(2,N1) GO TO 208

271

 

 

IF (HN1.GT/IT(5,N1) GO TO 209

281

 

 

DO 210 IN=2,10

291

 

 

IF (HN1.GE.F510(IN,N1)) GO TO 211

301

210

 

CONTINUE

311

211

 

WA=F510(IN,N1)

321

 

 

WB=F510(IN-1,N1)

331

 

 

XA=-5.0-0.5*IN

341

 

 

XB=XA+0.5

351

 

 

GO TO 200

361

209

 

DO 212 IN=2,15

371

 

 

IF (HN1.GE.F25(IN,N1)) GO TO 213

381

212

 

CONTINUE

391

213

 

WA=F25(IN,N1)

401

 

 

WB=F25(IN-1,N1)

411

 

 

XA=-2.0-0.2*IN

421

 

 

XB=XA+0,2

431

 

 

GO TO 200

441

208

 

IF (HN1.GE.F02(1,N1)) GO TO 216

451

 

 

GO TO 217

461

216

 

WA=F02(1,N1)

471

 

 

WB=N1(N1)

481

 

 

XA=-0.1

491

 

 

XB=0.0

501

 

 

GO TO 200

511

217

 

DO 214 IN=2,20

521

 

 

IF (HN1.GE.F02(IN,N1)) GO TO 213

531

214

 

CONTINUE

541

215

 

WA=F02(IN,N1)

551

 

 

WB=F02(IN-1,N1)

561

 

 

XA=-0.1+IN

571

 

 

XB=XA+0,1

581

200

 

UL1=(HN1-WA)/(WB-WA)+(XB-XA)+XA

591

 

 

RETURN

601

 

 

END

*** HF ***

1

 

FUNCTION HF(N,U1)

11

 

COMMO\ /RFIS/BCK(10),CCK(10),WC(10),UKR(10)

21

 

COMMO\ /CUEI/A(10),B(10),C(10)

22

 

COMMO\ /CONST/RV,CL,CB,TFR,TK

31

 

D=B(N)

41

 

D1=TK-J

51

 

D2=U1-TK

61

 

HF=BCK(N)*((CCK(N)+RV*CL+WC(A)+C(B)+(CB-R/*CL))*

71

3

D1+\((RV*CL-CB)*ALOG(ABC(U1-D)/C1))*D2/((J1-D)*D1)))

81

 

RETURN

91

 

END


Решается задача в 3 измерениях

направление осей К

0----J

/

/

/

/

I

Шаг по времени в часах                     121.67

Граница измерений             876000.00

Число блоков по I       7

Номер блока

1

2

3

4

5

6

7

Размер блока

2.00

2.00

2.00

2.00

3.00

3.00

3.00

Коорд. Центра

1.00

3.00

5.00

7.00

9.50

12.50

15.50

Число блоков по J      8

Номер блока

1

2

3

4

5

6

7

8

Размер блока

3.00

3.00

2.00

2.00

2.00

3.00

3.00

3.00

Коорд. центра

1.50

4.50

7.00

9.00

11.00

13.50

16.50

19.50

Число блоков по K     6

Номер блока

1

2

3

4

5

6

Размер блока

3.00

3.00

3.00

2.00

2.00

2.00

Коорд. центра

1.50

4.50

7.50

10.00

12.00

14.00

 

Условия на границах

K/J

1

2

3

4

5

6

7

8

6

121

121

121

121

111

111

111

111

5

121

121

121

121

131

111

111

111

4

131

131

131

131

131

131

111

111

3

111

111

131

131

131

131

131

131

2

111

111

111

111

111

111

131

131

1

111

111

111

111

111

111

111

111

Шифр XYZW если X=0 - нет источника

Y - номер типа источника

Z - номер зоны краевых условий на верхней границе

W - тоже на нижней

На верхней границе

Номер зоны / месяц

1

2

3

4

5

6

7

8

9

10

11

12

1

темп. возд.

-4.00

-23.20

-19.20

 -9.50

  3.70

13.10

10.00

12.40

  4.60

 -4.60

 16.40

-22.40

 

сопротивл.

  2.00

   2.70

   2.80

  3.00

  1.70

  0.0

  0.0

  0.0

  0.0

  1.30

   2.10

   2.30

 

теплоотд.

13.80

 12.30

 12.30

12.30

12.10

12.40

13.40

13.00

14.30

13.40

 13.50

 12.70

Номер зоны / месяц

1

2

3

4

5

6

7

8

9

10

11

12

2

темп. возд.

20.00

 20.00

 20.00

20.00

20.00

20.00

20.00

20.00

20.00

20.00

 20.00

 20.00

 

сопротивл.

  1.00

   1.00

   1.00

  1.00

  1.00

  1.00

  1.00

  1.00

  1.00

  1.00

   1.00

   1.00

 

теплоотд.

20.00

 20.00

 20.00

20.00

20.00

20.00

20.00

20.00

20.00

20.00

 20.00

 20.00

Номер зоны / месяц

1

2

3

4

5

6

7

8

9

10

11

12

2

темп. возд.

-24.60

-23.20

-19.20

 -9.50

  3.70

13.10

10.00

12.40

  4.60

 -4.00

-16.40

-22.40

 

сопротивл.

    4.40

   4.40

   4.70

  5.00

  2.07

  0.0

  0.0

  0.0

  0.0

  2.14

   3.57

   3.09

 

теплоотд.

  13.80

 12.30

 12.30

12.00

12.10

12.40

13.40

13.00

14.30

13.00

 13.50

 12.70

Продолжительность месяца           730.00

Периодичность измерения           8760.00

         Температура на нижней границе

Номер зоны                 1

Температура                -1.40

         На боковых

TBJ=0.0     TBJK=0.0  TBKH=0.0 TBKK=0.0

Номер слоя

ОБ   ВЕС   СК

С    УД

W   C

T   WJ

ТЕПЛ   М

ТЕПЛ   Т

A

B

C

1

1390.00

0.22

0.25

0.0

1.30

1.15

0.127

1.487

0.069

2

1520.00

0.24

0.22

0.0

1.55

1.45

0.045

0.301

0.046

3

1500.00

0.23

0.27

0.0

2.35

2.15

0.0

1.00

0.0

Начальная температура грунта

         К=1            1.50 м

Блоки по J

1

2

3

4

5

6

7

8

По I / метры

1.50

4.50

7.00

9.00

11.00

13.50

16.50

19.50

1

1.00

-0.10

-0.10

-0.10

-0.10

-0.10

-0.10

-0.10

-0.10

 

 

1

1

1

1

1

1

1

1

2

3.00

-0.40

-0.40

-0.40

-0.40

-0.40

-0.40

-0.40

-0.40

 

 

3

3

2

2

2

2

2

2

3

5.00

-1.00

-1.00

-1.00

-1.00

-1.00

-1.00

-1.00

-1.00

 

 

3

3

2

2

2

2

2

2

4

7.00

-1.20

-1.20

-1.20

-1.20

-1.20

-1.20

-1.20

-1.20

 

 

3

3

2

2

2

2

2

2

5

9.50

-1.40

-1.40

-1.40

-1.40

-1.40

-1.40

-1.40

-1.40

 

 

3

3

2

2

2

2

2

2

6

12.50

-1.40

-1.40

-1.40

-1.40

-1.40

-1.40

-1.40

-1.40

 

 

3

3

3

3

3

3

3

3

7

15.50

-1.40

-1.40

-1.40

-1.40

-1.40

-1.40

-1.40

-1.40

 

 

3

3

3

3

3

3

3

3

         К=4            10.00 м

Блоки по J

1

2

3

4

5

6

7

8

По I / метры

1.50

4.50

7.00

9.00

11.00

13.50

16.50

19.50

1

1.00

-0.10

-0.10

-0.10

-0.10

-0.10

-0.10

-0.10

-0.10

 

 

1

1

1

1

1

1

1

1

2

3.00

-0.40

-0.40

-0.40

-0.40

-0.40

-0.40

-0.40

-0.40

 

 

3

3

2

2

2

2

2

2

3

5.00

-1.00

-1.00

-1.00

-1.00

-1.00

-1.00

-1.00

-1.00

 

 

3

3

2

2

2

2

2

2

4

7.00

-1.20

-1.20

-1.20

-1.20

-1.20

-1.20

-1.20

-1.20

 

 

3

3

2

2

2

2

2

2

5

9.50

-1.40

-1.40

-1.40

-1.40

-1.40

-1.40

-1.40

-1.40

 

 

3

3

2

2

2

2

2

2

6

12.50

-1.40

-1.40

-1.40

-1.40

-1.40

-1.40

-1.40

-1.40

 

 

3

3

3

3

3

3

3

3

7

15.50

-1.40

-1.40

-1.40

-1.40

-1.40

-1.40

-1.40

-1.40

 

 

3

3

3

3

3

3

3

3

         К=6            14.00 м

Блоки по J

1

2

3

4

5

6

7

8

По I / метры

1.50

4.50

7.00

9.00

11.00

13.50

16.50

19.50

1

1.00

-0.10

-0.10

-0.10

-0.10

-0.10

-0.10

-0.10

-0.10

 

 

1

1

1

1

1

1

1

1

2

3.00

-0.40

-0.40

-0.40

-0.40

-0.40

-0.40

-0.40

-0.40

 

 

3

3

2

2

2

2

2

2

3

5.00

-1.00

-1.00

-1.00

-1.00

-1.00

-1.00

-1.00

-1.00

 

 

3

3

2

2

2

2

2

2

4

7.00

-1.20

-1.20

-1.20

-1.20

-1.20

-1.20

-1.20

-1.20

 

 

3

3

2

2

2

2

2

2

5

9.50

-1.40

-1.40

-1.40

-1.40

-1.40

-1.40

-1.40

-1.40

 

 

3

3

2

2

2

2

2

2

6

12.50

-1.40

-1.40

-1.40

-1.40

-1.40

-1.40

-1.40

-1.40

 

 

3

3

3

3

3

3

3

3

7

15.50

-1.40

-1.40

-1.40

-1.40

-1.40

-1.40

-1.40

-1.40

 

 

3

3

3

3

3

3

3

3

         J=1             1.50 м

Блоки по K

1

2

3

4

5

6

По I / метры

1.50

4.50

7.00

10.00

12.00

14.00

1

1.00

-0.10

-0.10

-0.10

-0.10

-0.10

-0.10

 

 

1

1

1

1

1

1

2

3.00

-0.40

-0.40

-0.40

-0.40

-0.40

-0.40

 

 

3

3

3

3

3

3

3

5.00

-1.00

-1.00

-1.00

-1.00

-1.00

-1.00

 

 

3

3

3

3

3

3

4

7.00

-1.20

-1.20

-1.20

-1.20

-1.20

-1.20

 

 

3

3

3

3

3

3

5

9.50

-1.40

-1.40

-1.40

-1.40

-1.40

-1.40

 

 

3

3

3

3

3

3

6

12.50

-1.40

-1.40

-1.40

-1.40

-1.40

-1.40

 

 

3

3

3

3

3

3

7

15.50

-1.40

-1.40

-1.40

-1.40

-1.40

-1.40

 

 

3

3

3

3

3

3

         J=4             9.00 м

Блоки по K

1

2

3

4

5

6

По I / метры

1.50

4.50

7.00

10.00

12.00

14.00

1

1.00

-0.10

-0.10

-0.10

-0.10

-0.10

-0.10

 

 

1

1

1

1

1

1

2

3.00

-0.40

-0.40

-0.40

-0.40

-0.40

-0.40

 

 

2

2

2

2

2

2

3

5.00

-1.00

-1.00

-1.00

-1.00

-1.00

-1.00

 

 

2

2

2

2

2

2

4

7.00

-1.20

-1.20

-1.20

-1.20

-1.20

-1.20

 

 

2

2

2

2

2

2

5

9.50

-1.40

-1.40

-1.40

-1.40

-1.40

-1.40

 

 

2

2

2

2

2

2

6

12.50

-1.40

-1.40

-1.40

-1.40

-1.40

-1.40

 

 

3

3

3

3

3

3

7

15.50

-1.40

-1.40

-1.40

-1.40

-1.40

-1.40

 

 

3

3

3

3

3

3

         J=6             13.50 м

Блоки по K

1

2

3

4

5

6

По I / метры

1.50

4.50

7.00

10.00

12.00

14.00

1

1.00

-0.10

-0.10

-0.10

-0.10

-0.10

-0.10

 

 

1

1

1

1

1

1

2

3.00

-0.40

-0.40

-0.40

-0.40

-0.40

-0.40

 

 

2

2

2

2

2

2

3

5.00

-1.00

-1.00

-1.00

-1.00

-1.00

-1.00

 

 

2

2

2

2

2

2

4

7.00

-1.20

-1.20

-1.20

-1.20

-1.20

-1.20

 

 

2

2

2

2

2

2

5

9.50

-1.40

-1.40

-1.40

-1.40

-1.40

-1.40

 

 

2

2

2

2

2

2

6

12.50

-1.40

-1.40

-1.40

-1.40

-1.40

-1.40

 

 

3

3

3

3

3

3

7

15.50

-1.40

-1.40

-1.40

-1.40

-1.40

-1.40

 

 

3

3

3

3

3

3

         J=8             19.50 м

Блоки по K

1

2

3

4

5

6

По I / метры

1.50

4.50

7.00

10.00

12.00

14.00

1

1.00

-0.10

-0.10

-0.10

-0.10

-0.10

-0.10

 

 

1

1

1

1

1

1

2

3.00

-0.40

-0.40

-0.40

-0.40

-0.40

-0.40

 

 

2

2

2

2

2

2

3

5.00

-1.00

-1.00

-1.00

-1.00

-1.00

-1.00

 

 

2

2

2

2

2

2

4

7.00

-1.20

-1.20

-1.20

-1.20

-1.20

-1.20

 

 

2

2

2

2

2

2

5

9.50

-1.40

-1.40

-1.40

-1.40

-1.40

-1.40

 

 

2

2

2

2

2

2

6

12.50

-1.40

-1.40

-1.40

-1.40

-1.40

-1.40

 

 

3

3

3

3

3

3

7

15.50

-1.40

-1.40

-1.40

-1.40

-1.40

-1.40

 

 

3

3

3

3

3

3

Шаг печати результата час      8760.00

К=1       1.50 м        Время Т=8760.16 часов       1 Периода 0 Интервала       0.16 часа

Блоки по J

1

2

3

4

5

6

7

8

По I / метры

1.50

4.50

7.00

9.00

11.0

13.50

16.50

19.50

1

1.00

761.25

761.55

761.91

762.71

763.01

733.20

763.53

763.66

2

3.00

-0.38

-0.44

-0.62

-0.59

-0.57

-0.57

-0.55

-0.55

3

5.00

-0.66

-0.73

-0.92

-0.91

-0.90

-0.59

-0.89

-0.89

4

7.00

-0.90

-0.96

-1.11

-1.12

-1.11

-1.11

-1.11

-1.11

5

9.50

-1,13

-1.17

-1.27

-1.28

-1.29

-1.29

-1.29

-1.29

6

12.50

-1,28

-1.30

-1.33

-1.34

-1.34

-1.35

-1.35

-1.35

7

15.50

-1,37

-1.37

-1.38

-1.38

-1.38

-1.38

-1.38

-1.38

К=4       10.00 м      Время Т=8760.16 часов       1 Периода 0 Интервала       0.16 часа

Блоки по J

1

2

3

4

5

6

7

8

По I / метры

1.50

4.50

7.00

9.00

11.0

13.50

16.50

19.50

1

1.00

796.28

796.39

795.93

793.60

783.28

776.98

764.07

763.82

2

3.00

702.42

701.25

-3.35

-0.38

-0.45

-0.50

-0.54

-0.54

3

5.00

-0.43

-0.50

-0.78

-0.83

-0.86

-0.87

-0.88

-0.89

4

7.00

-0.74

-0.82

-1.23

-1.09

-1.10

-1.11

-1.11

-1.11

5

9.50

-1,04

-1.09

-1.23

-1.27

-1.29

-1.29

-1.29

-1.29

6

12.50

-1,24

-1.25

-1.31

-1.33

-1.34

-1.35

-1.35

-1.35

7

15.50

-1,35

-1.36

-1.37

-1.38

-1.38

-1.38

-1.38

-1.38

К=6       14.00 м      Время Т=8760.16 часов       1 Периода 0 Интервала       0.16 часа

Блоки по J

1

2

3

4

5

6

7

8

По I / метры

1.50

4.50

7.00

9.00

11.0

13.50

16.50

19.50

1

1.00

5.62

5.66

5.85

4.78

786.13

765.26

763.34

763.06

2

3.00

727.45

726.35

-0.06

-0.11

-0.40

-0.54

-0.57

-0.57

3

5.00

-0.39

-0.45

-0.69

-0.77

-0.84

-0.88

-0.89

-0.90

4

7.00

-0.70

-0.78

-3.99

-1.07

-1.10

-1.11

-1.11

-1.11

5

9.50

-1,01

-1.07

-1.21

-1.26

-1.28

-1.29

-1.29

-1.29

6

12.50

-1,22

-1.28

-1.30

-1.32

-1.34

-1.35

-1.35

-1.35

7

15.50

-1,35

-1.36

-1.37

-1.38

-1.38

-1.38

-1.38

-1.38

J=1        1.50 м        Время Т=8760.16 часов       1 Периода 0 Интервала       0.16 часа

Блоки по K

1

2

3

4

5

6

По I / метры

1.50

4.50

7.50

10.00

12.0

14.00

1

1.00

761.25

762.36

766.28

796.28

4.54

5.62

2

3.00

-0.38

-0.35

-0.26

702.42

719.71

727.45

3

5.00

-0.66

-0.62

-0.54

-0.43

-0.40

-0.39

4

7.00

-0.90

-0.56

-0.88

-0.74

-0.72

-0.78

5

9.50

-1,13

-1.11

-1.67

-1.64

-1.62

-1.01

6

12.50

-1,28

-1.27

-1.23

-1.24

-1.23

-1.22

7

15.50

-1,37

-1.36

-1.36

-1.33

-1.35

-1.35

J=4        9.00 м        Время Т=8760.16 часов       1 Периода 0 Интервала       0.16 часа

Блоки по K

1

2

3

4

5

6

По I / метры

1.50

4.50

7.50

10.00

12.0

14.00

1

1.00

762.71

763.37

778.01

793.60

3.83

4.70

2

3.00

-0.59

-0.57

-0.49

-0.38

-0.15

-0.11

3

5.00

-0.91

-0.90

-0.87

-0.83

-0.79

-0.77

4

7.00

-1.12

-1.11

-1.10

-1.09

-1.07

-1.07

5

9.50

-1,28

-1.28

-1.28

-1.27

-1.27

-1.26

6

12.50

-1,34

-1.34

-1.33

-1.33

-1.32

-1.32

7

15.50

-1,38

-1.38

-1.38

-1.38

-1.38

-1.38

J=6        13.50 м      Время Т=8760.16 часов       1 Периода 0 Интервала       0.16 часа

Блоки по K

1

2

3

4

5

6

По I / метры

1.50

4.50

7.50

10.00

12.0

14.00

1

1.00

763.20

763.88

776.87

776.98

765.60

765.26

2

3.00

-0.57

-0.54

-0.49

-0.50

-0.52

-0.54

3

5.00

-0.89

-0.89

-0.87

-0.87

-0.88

-0.88

4

7.00

-1.11

-1.11

-1.11

-1.11

-1.11

-1.11

5

9.50

-1,29

-1.29

-1.29

-1.29

-1.29

-1.29

6

12.50

-1,35

-1.35

-1.35

-1.35

-1.35

-1.35

7

15.50

-1,38

-1.38

-1.38

-1.38

-1.38

-1.38

J=8        19.50 м      Время Т=8760.16 часов       1 Периода 0 Интервала       0.16 часа

Блоки по K

1

2

3

4

5

6

По I / метры

1.50

4.50

7.50

10.00

12.0

14.00

1

1.00

763.66

776.71

776.73

763.82

763.27

763.86

2

3.00

-0.55

-0.49

-0.49

-0.54

-0.57

-0.57

3

5.00

-0.89

-0.87

-0.87

-0.89

-0.89

-0.90

4

7.00

-1.11

-1.11

-1.11

-1.11

-1.11

-1.11

5

9.50

-1,29

-1.29

-1.29

-1.29

-1.29

-1.29

6

12.50

-1,35

-1.35

-1.35

-1.35

-1.35

-1.35

7

15.50

-1,38

-1.38

-1.38

-1.38

-1.38

-1.38

К=1       1.50 м        Время Т=43800.54 часов     5 Периода 0 Интервала       0.54 часа

Блоки по J

1

2

3

4

5

6

7

8

По I / метры

1.50

4.50

7.00

9.00

11.0

13.50

16.50

19.50

1

1.00

770.37

770.57

771.16

771.22

771.20

771.28

771.83

772.03

2

3.00

714.69

709.65

-0.08

-0.12

-0.13

-0.13

-0.12

-0.11

3

5.00

-0.25

-0.29

-0.42

-0.58

-0.53

-0.55

-0.54

-0.54

4

7.00

-0.48

-0.53

-0.66

-0.74

-0.78

-0.81

-0.81

-0.82

5

9.50

-0,75

-0.89

-0.90

-0.97

-1.01

-1.04

-1.06

-1.06

6

12.50

-1,03

-1.26

-1.11

-1.15

-1.17

-1.20

-1.21

-1.21

7

15.50

-1,28

-1.29

-1.30

-1.31

-1.32

-1.33

-1.34

-1.34

К=4       10.00 м      Время Т=43800.54 часов     5 Периода 0 Интервала       0.54 часа

Блоки по J

1

2

3

4

5

6

7

8

По I / метры

1.50

4.50

7.00

11.00

11.0

13.50

16.50

19.50

1

1.00

0.69

0.76

0.85

0,46

794.41

785.72

772.91

772.38

2

3.00

1.01

1.06

1.20

0.95

0.59

-0.01

-0.09

-0.10

3

5.00

710.10

708.61

-0.28

-0.16

-0.29

-0.48

-0.52

-0.54

4

7.00

-0.33

-0.37

-0.51

-0.60

-0.69

-0.77

-0.80

-0.81

5

9.50

-0,67

-0.71

-0.83

-0.92

-0.97

-1.02

-1.05

-1.06

6

12.50

-0,99

-1.02

-1.08

-1.12

-1.16

-1.19

-1.20

-1.21

7

15.50

-1,27

-1.28

-1.30

-1.31

-1.32

-1.33

-1.33

-1.34

К=6       14.00 м      Время Т=43800.54 часов     5 Периода 0 Интервала       0.54 часа

Блоки по J

1

2

3

4

5

6

7

8

По I / метры

1.50

4.50

7.00

9.00

11.0

13.50

16.50

19.50

1

1.00

7.92

8.02

8.25

6.67

799.62

775.78

771.71

771.12

2

3.00

3.17

3.26

3.63

2.92

1.00

701.03

-0.12

-0.14

3

5.00

793.32

796.56

0.87

0.63

-0.14

-0.45

-0.54

-0.56

4

7.00

-0.50

-0.32

-0.35

-0.40

-0.62

-0.75

-0.80

-0.82

5

9.50

-0,64

-0.63

-0.79

-0.88

-0.95

-1.01

-1.05

-1.06

6

12.50

-0,98

-1.01

-1.07

-1.11

-1.15

-1.18

-1.20

-1.21

7

15.50

-1,26

-1.27

-1.29

-1.30

-1.32

-1.33

-1.33

-1.34

J=1        1.50 м                  Время Т=43800.54 часов     5 Периода 0 Интервала       0.54 часа

Блоки по K

1

2

3

4

5

6

По I / метры

1.50

4.50

7.50

10.00

12.0

14.00

1

1.00

770.37

770.71

775.00

0.69

6.47

7.92

2

3.00

714.69

716.33

736.25

1.01

2.53

3.17

3

5.00

-0.25

-0.23

-0.18

710.18

760.45

793.32

4

7.00

-0.48

-0.46

-0.41

-0.33

-0.31

-0.30

5

9.50

-0,75

-0.73

-0.70

-0.67

-0.65

-0.64

6

12.50

-1,03

-1.02

-1.01

-0.99

-0.98

-0.98

7

15.50

-1,28

-1.28

-1.27

-1.27

-1.27

-1.26

J=4        9.00 м                  Время Т=43800.54 часов     5 Периода 0 Интервала                0.54 часа

Блоки по K

1

2

3

4

5

6

По I / метры

1.50

4.50

7.50

10.00

12.0

14.00

1

1.00

771.22

772.42

788.84

0.46

5.55

6.67

2

3.00

-0.12

-0.09

0.11

0.95

2.30

2.92

3

5.00

-0.50

-0.48

-0.40

-0.163

701.66

0.63

4

7.00

-0.74

-0.73

-0.68

-0.60

-0.54

-0.48

5

9.50

-0,97

-0.96

-0.94

-0.92

-0.90

-0.88

6

12.50

-1,15

-1.14

-1.13

-1.12

-1.12

-1.11

7

15.50

-1,31

-1.31

-1.31

-1.31

-1.31

-1.30

J=6        13.50 м      Время Т=438000.54 часов   5 Периода 0 Интервала       0.54 часа

Блоки по K

1

2

3

4

5

6

По I / метры

1.50

4.50

7.50

10.00

12.0

14.00

1

1.00

771.28

772.50

785.02

785.72

776.63

775.78

2

3.00

-0.13

-0.10

-0.03

-0.01

-0.00

701.03

3

5.00

-0.55

-0.53

-0.49

-0.48

-0.46

-0.45

4

7.00

-0.81

-0.80

-0.78

-0.77

-0.76

-0.75

5

9.50

-1,04

-1.04

-1.03

-1.02

-1.02

-1.01

6

12.50

-1,20

-1.19

-1.19

-1.19

-1.18

-1.18

7

15.50

-1,33

-1.33

-1.33

-1.33

-1.33

-1.33

J=8        19.50 м      Время Т=438000.54 часов   5 Периода 0 Интервала       0.54 часа

Блоки по K

1

2

3

4

5

6

По I / метры

1.50

4.50

7.50

10.00

12.0

14.00

1

1.00

772.03

784.57

784.59

772.38

771.40

771.12

2

3.00

-0.11

-0.04

-0.04

-0.10

-0.13

-0.14

3

5.00

-0.54

-0.51

-0.51

-0.54

-0.55

-0.56

4

7.00

-0.82

-0.80

-0.80

-0.81

-0.82

-0.82

5

9.50

-1,06

-1.06

-1.06

-1.06

-1.06

-1.06

6

12.50

-1,21

-1.21

-1.21

-1.21

-1.21

-1.21

7

15.50

-1,34

-1.34

-1.348

-1.34

-1.34

-1.34

 


 

СОДЕРЖАНИЕ

1. Общие положения

2. Алгоритм и структура программы

2.1. Математическая постановка задачи

2.2. Алгоритм решения задачи

2.3. Структура программы

3. Подготовка исходных данных для расчета

3.1. Расчетная область

3.2. Разбивка области исследования

3.3. Начальные условия

3.4. Граничные условия

3.5. Физические и теплофизические свойства

3.6. Дополнительные данные, необходимые для проведения расчета

4. Проведение расчета

4.1. Вызов и загрузка программы в ЭВМ

4.2. Составление программы MAIN и подпрограммы BLOCK DATA

4.3. Подготовка к вводу группы данных с перфокарт

4.4. Выходные данные

4.5. Подготовка данных для расчета по программе PROGNOZ-2S

4.6. Подготовка данных для расчета по программе PROGNOZ-3N

4.7. Подготовка данных для расчета по программе PROGNOZ-2N

4.8. Подготовка данных для расчета по программе PROGNOZ-L

Приложение 1

Пример расчета

Исходные данные для расчета контрольного примера приложения 1 по программе PROGNOZ-3S

Выдача на печать начальных данных подготовленных для расчета по программе PROGNOZ-3S контрольного примера приложения 1

Приложение 2

Пример текста программы ISTO

Приложение 3

Текст программы PROGNOZ-3S с результатами расчета контрольного примера приложения 1

 

 



© 2013 Ёшкин Кот :-)