Методы оперативной идентификации тепловых режимов в конструкциях летательных аппаратов тема диссертации и автореферата по ВАК РФ 05.07.11, кандидат технических наук Геджадзе, Игорь Юрьевич

• Геджадзе, Игорь Юрьевич
• кандидат технических наук
• 1998

Оглавление диссертации кандидат технических наук Геджадзе, Игорь Юрьевич

ОГЛАВЛЕНИЕ

Введение

Глава 1. Алгоритм наблюдения (линейная модель)

1.1. Постановка задачи

1.2. Метод решения линейной ОЗТ

1.3. Анализ точности решения и правило выбора параметра регуляризации

1.4. Описание алгоритма и моделирование

1.5. Основные результаты главы

Глава 2. Задача наблюдения: особенности измерительной

схемы

2.1. О необходимости планирования

2.2. О практических условиях единственности

2.3. Оптимизация измерительной схемы

2.4. Границы адекватности одномерной постановки

2.5. Основные результаты главы

Глава 3. Алгоритм наблюдения (квазилинейная модель)

3.1. Постановка задачи

3.2. Метод последовательных приближений

3.3. О сходимости метода последовательных приближений

3.4. Особенности практической реализации

3.5. Алгоритм и результаты моделирования

3.6. Обработка данных теплового эксперимента

3.7. Основные результаты главы

Глава 4. Безитерационный метод решения нелинейной ОЗТ

4.1. Постановка задачи

4.2. Решение расширенной линейной ОЗТ

4.3. Численный анализ

4.4. О выборе параметров

4.5. Основные результаты главы

Заключение

Приложение А. Экономичный алгоритм решения двухпараметри-

ческой задачи наименьших квадратов

А.1. Постановка задачи

А.2. Декомпозиция

А.З. Задача наименьших квадратов

А.4. Алгоритм

Список литературы

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

Введение диссертации (часть автореферата) на тему «Методы оперативной идентификации тепловых режимов в конструкциях летательных аппаратов»

- 3-Введение

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

Перспективным направлением в исследовании и отработке тепло-нагруженных конструкций, а также технологий производства, связанных с реализацией требуемых тепловых процессов, является методология, основанная на решении обратных задач теплообмена (ОЗТО) [1]. В дальнейшем под ОЗТО будем понимать в основном задачи диагностики и идентификации параметров тепловых процессов. Соответственно трем основным способам теплообмена рассматриваются обратные задачи: теплопроводности (ОЗТ), конвективного теплообмена, радиационного теплообмена. ОЗТ, в свою очередь, делятся на граничные, ретроспективные, геометрические и коэффициентные. Граничные задачи заключаются в нахождении функций и параметров, входящих в граничные условия, ретроспективные - в нахождении начальных условий, геометрические - в реконструировании геометрических характеристик границы области, коэффициентные - в нахождении коэффициентов уравнений в частных производных, описывающих процесс тепло-

обмена. Данная методология успешно развивается на протяжении двух последних десятилетий и нашла практическое применение в различных областях техники. Основная причина столь широкого применения связана, в частности, с возможностью учета эффектов нестационарности, нелинейности, и пространственного характера теплообменных процессов. Кроме того, данный подход позволяет проводить экспериментальные исследования в условиях, близких к натурным или непосредственно на этапе эксплуатации, является более информативным по сравнению с классическими методами, что позволяет ускорить проведение экспериментальных работ и тем самым уменьшить цикл отработки нового изделия. В последние годы проводится также работа по созданию автоматизированных систем, реализующих методологию ОЗТО, объединенных под управлением единого монитора в рамках персонального автоматизированного рабочего места (ПАРМ) инженера-исследователя тепловых процессов [30].

За последнее десятилетие произошел беспрецендентный скачок в развитии вычислительных средств. Так, на сегодняшний день, большинство алгоритмов ОЗТО, требующих (по понятиям десятилетней давности) немалых вычислительных ресурсов, могут быть легко реализованы на персональном компьютере. Весьма впечатляют и возможности современных спецвычислителей, в частности, сигнальных процессоров 3-го поколения. Поэтому, не вызывает удивления тот факт, что одна из наиболее значимых тенденций в развитии современной техники заключается в насыщении технических объектов электронными системами, осуществляющими управление на всех уровнях и но все более сложным алгоритмам, т.е. в общем повышении уровня "интеллекта" технических систем. С учетом этой тенденции можно говорить также о системах автоматического управления (САУ) тепловыми процессами, в том числе о замкнутых, т.е. реализующих принцип обратной связи. Уточним, что речь идет о существенно нестационарных, быстротекущих тепловых процессах, несводимых к регулярным режи-

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

В качестве примера систем, где целесообразно использование замкнутых САУ можно привести различные варианты активных систем тепловой защиты (ТПЗ), таких как конвективное, пленочное и пористое охлаждение теплонагруженных конструкций. Последнюю рассмотрим более подробно [22]. Пористое охлаждение предполагается использовать, в частности, для летательных аппаратов, входящих в атмосферу с гиперзвуковой скоростью, в качестве альтернативы различным видам разрушаемой тепловой защиты, при использовании которых может изменяться геометрия носовой части объекта и, следовательно, условия обтекания, со всеми вытекающими отсюда проблемами для управления и обеспечения требуемой точности движения. Пористая ТПЗ может представлять собой набор пористых вставок, с заданным шагом расположенных на охлаждаемой поверхности ЛА, через которые продувает-

ся холодный газ. За счет большой поверхности теплообмена в пористом теле происходит интенсивный теплосъем с пористого каркаса, а охлаждение несущей поверхности происходит в результате оттока тепла в направлении пористых вставок. Кроме того, происходящий при этом вдув газа в набегающий поток приводит к увеличению толщины пограпслоя и уменьшению трения, что служит еще одним важным механизмом защиты. Однако, данная картина является лишь первым приближением, и может быть сильно искажена при протекании в пристеночном слое химических реакций, а также присутствием радиационной составляющей теплового потока, которая в ряде случаев вообще может оказаться доминирующей. Таким образом, имеет место комплексный физический процесс, настолько сложный, что существующие математические модели дают не более чем качественное представление о нем. Конкретная реализация данного процесса обусловлена влиянием множества внешних случайных факторов, от состояния атмосферы до флуктуаций угла атаки. В рассматриваемом примере, по-видимому, целесообразно применение замкнутой САУ, задача которой заключается в удержании температуры поверхности на уровне, гарантирующем ее неразрушение. Для работы такой системы необходимо знать в текущий момент времени и для заданного набора точек поверхности: температуру, тепловой поток по направлению нормали к поверхности, тепловой поток вдоль поверхности в направлении пористых вставок. Если перенос энергии излучением играет значительную роль необходимо также разделять конвективную и радиационную составляющие теплового потока. Все перечисленные величины могут быть определены на основании косвенных температурных измерений с использованием методов решения ОЗТ, однако для этого необходима такая модификация этих методов, которая позволила бы решать задачу наблюдения, т.е. идентифицировать заданные параметры в реальном масштабе времени. В качестве управляющих параметров могут служить расход охладителя и его химический состав. Использование подобной системы позво-

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

Другая область, где может оказаться необходимым управлять тепловыми нагрузками по принципу следящей системы или осуществлять непрерывную во времени температурную диагностику объекта - это тепловые эксперименты и испытания. К примеру, требуется воспроизводить на стендах заданные временные зависимости тепловых потоков или температур в образцах и моделях, получая информацию от контрольных датчиков. Для моделирования нестационарных тепловых потоков большой мощности образец материала может быть помещен в высокотемпературную дозвуковую газовую струю, вытекающую, например, из сопла плазмотрона, а интенсивность теплового воздействия может изменяться за счет регулирования по заданной программе расстояния от образца до сопла плазмотрона. Однако, первоначально необходимо провести тарировку, т.е. получить зависимости температуры струи и коэффициента теплоотдачи от координаты вдоль струи. Впоследствии, любое изменение условий эксперимента, например, смена сопла, изменение газового состава смеси, формы или размера образца ведет к необходимости заново проводить тарировку. В данном примере целесообразно построение замкнутой САУ, задача которой заключается в отработке заданной зависимости теплового потока к поверхности образца (в этом случае САУ входит в состав автоматизированной системы). Для этого в текущий момент времени необходимо знать температуру газовой струи и коэффициент теплоотдачи, либо непосредственно тепловой поток в образец. Как и прежде, перечисленные параметры могут быть получены (идентифицированы) с использованием соответствующих модификаций методов решения ОЗТ. Управляющим параметром служит координата расположения образца в газовой струе. Если же моделируемые процессы столь скоротечны, что быстродейст-

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

При проведении экспериментальных исследований целесообразно производить экспресс-обработку экспериментальных данных, даже при условии, что более полные и точные результаты будут получены по завершении эксперимента. Это позволяет следить за ходом испытаний и контролировать проходящие процессы, например, с целью выявления внештатной ситуации и своевременного их прекращения. Для этого, к примеру, требуется знать в текущий момент времени температуру поверхности объекта в наиболее опасной точке или сечении. В качестве иллюстрации рассмотрим проблему измерения высоких температур и тепловых потоков большой мощности. В настоящее время для измерений в основном используются следующие термопары: вольфрам-вольф-рамренивые до 2500К, платино-платинородиевые до 1800-1900К, хро-мель-алюмелевые до 1600-1700К, и хромель-копелевые до 1100К. Предположим, что задача состоит в измерении температуры металли-

кг

ческого расплава и для этой цели используется температурный зонд, который на короткое время вводится в расплав (при этом темпеатура поверхности зонда приблизительно становится равной температуре расплава), а затем убирается и охлаждается до начальной температуры (рис. 1, кр. 1). Допустим, что чувствительный элемент зонда можно рассматривать как пластину заданной толщины, на теплоизолированной поверхности которой располагается термопара. Как показано на рис. 1, при соответствующем выборе длительности зондирования температура в точке заделки термопары составляет порядка 25% от температуры поверхности зонда (кр. 2), однако этого вполне достаточно чтобы получить оценку температуры поверхности (кр. 3) из решения граничной ОЗТ.

1.2 г*т/Г'ТТ ' ' '""*'""

/ тих

10 —ГТ^Т---

08—£П---

о.б-4---

0.4-—-Т--

02—Г/^-^;—

0.0 ШШШ1\У--Ч......»■■■■.............^

-0.2'-----

0.0 0.4 0.8 1.2 1.6 2.0

Рис. 1. Оценка температуры поверхности зонда из решения ОЗТ.

(При расчетах использовались теплофизические свойства материала пластины характерные для молибдена. Для измерения температуры расплава стали в мартеновской печи можно использовать чувствительный элемент из молибдена, покрытый керамической пленкой на основе дисилицида молибдена или двуокиси циркония. Верхняя граница температуры поверхности такого элемента в агрессивной среде составляет около 2300К). Таким образом, использование методов решения ОЗТ позволяет определять значения высоких температур, в действительности измеряя температуры существенно более низкие. В рассматриваемом примере это позволяет применить вместо вольфрам-вольфрамре-ниевых термопар высокотехнологичные хромель-алюмелевые. При этом результирующая точность определения искомой температуры может оказаться существенно более высокой, чем при непосредственном измерении в высокотемпературной зоне. Следует также отметить, что если интервал зондирования мал настолько, что измеряемую температуру среды и коэффициент теплообмена на поверхности зонда можно считать стационарными на интервале, то нет необходимости ждать, пока температура поверхности зонда выровняется с температурой среды. В рассматриваемом случае задача САУ может состоять в управлении движением зонда, исходя из условия наилучшей идентификации искомых параметров, при ограничениях на температуру поверхности зонда и температуру в точках заделки термопар.

Т/Тта аг

|Д1

Заключение диссертации по теме «Тепловые режимы летательных аппаратов», Геджадзе, Игорь Юрьевич

4.5. Основные результаты главы

В данной главе рассматривается метод решения задачи наблюдения температурного поля, которое задается квазилинейным уравнением теплопроводности. Метод состоит в линеаризации исходной нелинейной постановки посредством замены нелинейных членов модели некоторыми неопределенными возмущениями, которые аппроксимируются линейными комбинациями базисных функций. Коэффициенты этих линейных комбинаций также считаются неизвестными величинами и подлежат идентификации наряду с искомыми краевыми условиями. Специфика задачи состоит в некорректности исходной постановки и необходимости регуляризации в виде (4.5). Выведено выражение для ошибки оценивания (4.33), и исследовано поведение компонент этой ошибки в зависимости от значений параметров а и 0 для ряда примеров. Показано, что в некоторой области изменения параметров ошибка оценивания мала, т.е. среди множества квазирешений, минимизирующих (4.5) существуют достаточно близкие к точному. Проблема, таким образом, состоит в нахождении оптимальных значений параметров. Один из возможных вариантов приводится в гл. 4.5, однако следует признать, что вопрос выбора параметров, а скорее эффективного алгоритма выбора параметров требует дальнейших исследований.

Замечание 1. Задача (4.2)-(4.4) может заменять задачу (3.8),

3.9) на первом этапе метода последовательных приближений (3.8)-(3.12), что обеспечивает скорейшую сходимость метода, однако данная замена эффективна лишь на первых приближениях, и в ряде случаев обеспечивает получение решения на первой же итерации.

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

3.10), функция рассматривается просто как дополнительный шумовой вход и может интерпретироваться как погрешность задания модели. В этом случае для выбора параметра регуляризации предлагается использовать обобщенный принцип невязки [29]. Еще один способ построения алгоритма устойчивого к ошибкам модели рассматривается в работе [31]. Способ строится на более тонком учете свойств невязки как функции, однако по сути не отличается от обобщенного принципа невязки. Не вдаваясь в подробности просто отметим, что использование этого принципа означает попытку подавить влияние (7(тд) только за счет увеличения а. На рис. 27-30 задаче (3.10) соответствуют значения коэффициентов при максимальном значении 0 = 104. Нетрудно видеть, что эффективность такого подхода значительно ниже, поскольку он приводит к чрезмерному росту смещения.

Замечание 3. Задача (4.2)-(4.4) может возникнуть не только при линеаризации нелинейной ОЗТ. Рассмотрим, к примеру, линейную задачу, однако предположим, что наблюдается температурное поле не бесконечной пластины, а цилиндра заданного радиуса. На боковой поверхности цилиндра происходит теплообмен такой интенсивности, что одномерная модель является неадекватной. В этом случае мы можем оставаться в рамках одномерной модели если введем в правую часть уравнения теплопроводности функцию 0{т,х), которая имеет смысл дополнительного источника или стока тепла. Далее необходимо использовать предлагаемый алгоритм, позволяющий наилучшим образом подавить влияние 0{%,х) на результат. При этом, для выбора оптимальных значений параметров необходима некоторая дополнительная информация относительно функции <2(т,х). Эта информация, однако, может носить достаточно общий и приблизительный характер, и может быть получена на основе каких-либо априорных соображений, либо в результате дополнительных измерений и т.д.

Замечание 4. Для наблюдения нестационарного температурного поля в случае линейной модели процесса мы используем, по сути дела, регуляризированный линейный цифровой фильтр, передаточные функции которого зависят от параметра регуляризации. Изложенный выше способ линеаризации позволяет использовать линейный цифровой фильтр также в случае, когда модель процесса нелинейна, однако передаточные функции этого фильтра зависят уже от двух параметров: параметра регуляризации и весового параметра, отвечающего за подавление влияния неопределенной возмущающей функции. Очевидно, что этот метод является приближенным.

-116 -Заключение

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

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

Выполненные в диссертационной работе исследования имеют целью специальную адаптацию методов решения ОЗТ для приложения к задачам наблюдения и экспресс-обработки результатов измерений. Наиболее существенными из полученных в работе результатов являются следующие:

1. Предложен и обоснован метод построения последовательной процедуры решения граничной ОЗТ, состоящий в переходе к гранично-ретроспективной постановке задачи на малом временном интервале, предшествующем текущему моменту времени. Данный метод наиболее эффективен именно для существенно нестационарных (несводимых к регулярным режимам) тепловых процессов, характерных для конструкций ЛА. Способ имеет ряд "технологических" преимуществ и наиболее подходит для создания системы наблюдения.

2. Впервые в практике решения ОЗТ для оценки точности решения использован подход, основанный на анализе передаточных функций регуляризированного фильтра. Тем самым, получен объективный критерий качества алгоритмов, применяемый по ходу работы для оптимизации и выбора параметров.

3. Предложен метод выбора параметра регуляризации, актуальный в случае, когда оценка вычисляется на базе выборки малого объема. Метод основан на статистическом критерии [11], однако область изменения параметра ограничивается исходя из необходимости достижения некоторой гарантированнной точности решения.

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

5. Предложен эффективный алгоритм решения нелинейной гранично-ретроспективной ОЗТ на базе метода последовательных приближений и разработаны основные детали его использования. Необходимость решения нелинейной задачи связана с тем, что теплофизические характеристики большинства материалов, используемых в конструкциях JIA, зависят от температуры.

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

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

8. С учетом рекомендаций и результатов работы в МАИ разработан и изготовлен охлаждаемый датчик теплового потока, рассчитанный на длительное наблюдение нестационарных потоков мощностью до

1.5 х 104 Квт/м2 и содержащих частоты до 5Гц. Комбинация двух таких датчиков, использующих различные режимы охлаждения, позволяет наблюдать параметры конвективного теплообмена.

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

Теперь еще раз подчеркнем некоторые специфические проблемы, возникающие при решении задачи наблюдения, которые могут составить предмет дальнейших исследований. Прежде всего отметим необходимость введения запаздывания при наблюдении граничного условия. Если расчетное время можно сократить упрощая алгоритм, оптимизируя вычислительный процесс, или просто применяя более быстродействующее устройство, то величина вынужденного запаздывания выбирается только исходя из требований к точности. Поэтому далее встает вопрос о построении эффективного управления по ретроспективной оценке. Еще одна проблема заключается в необходимости иметь достаточно хорошее представление о наблюдаемых процессах и измерительных шумах. К сожалению, построить систему наблюдения, которая позволила бы одинаково хорошо восстанавливать любые граничные условия при произвольных шумах, невозможно. Сказанное безусловно относится и к алгоритмам, работающим по полной выборке, которые рассматриваются нами как эталон, однако для последовательного алгоритма это гораздо более актуально. В то же время можно предполагать, что объектом наблюдения вряд ли станут совершенно неизученные тепловые процессы, поэтому задача скорее будет состоять в наблюдении отклонений от ожидаемого хода событий. То же относится и к вопросу изучения погрешности измерений. Можно видеть (см. рис. 6(6)), что существует ярко выраженная область наиболее опасных частот в шуме, которая с увеличением интервала оценивания "размазывается" по частотной оси. Наличие такой области может иметь как положительные, так и отрицательные последствия в зависимости от факта реального присутствия этих частот в шуме. Поэтому возникает проблема анализа шума с целью выявления и возможного устранения источников, генерирующих наиболее опасные компоненты. С другой стороны, при уменьшении интервала оценивания снижается влияние на точность оценивания низких частот, в том числе систематических составляющих и вековых членов (уход нулей). Все эти особенности должны быть учтены при проектировании как ЧЭ датчиков, так и всей последующей обрабатывающей цепочки.

Список литературы диссертационного исследования кандидат технических наук Геджадзе, Игорь Юрьевич, 1998 год

Список литературы

1. Алифанов О.М. Обратные задачи теплообмена. М., 1988.

2. Алифанов О.М., Артюхин Е.А. Определение граничных условий в процессе тепловых газодинамических испытаний. //ТВТ. 1978. Т. 16. №4. С.819-825.

3. Алифанов О.М., Артюхин Е.А., Румянцев C.B. Экстремальные методы решения некорректных задач. Москва: Наука, 1988.

4. Алифанов О.М., Артюхин Е.А., Логинов С.Н., Малоземов В.В. К вопросу решения обратной задачи теплопроводности методом динамической фильтрации. //ИФЖ, T.XLI, №5, 1981.

5. Андерсон Т. Статистический анализ временных рядов. Москва: Мир. 1976.

6. Бек Дж., Блакуэлл Б., Сент-Клер Ч. Некорректные обратные задачи теплопроводности. М.: Мир. 1989.

7. Богуславский И.А. Робастный алгоритм рекуррентной фильтрации. //ДАН СССР. 1988. Т.300. №6.

8. Богуславский И.А. Аппроксимация задачи нелинейной фильтрации последовательностью линейных задач; мера робастности фильтрации. //Изв. РАН. Техническая кибернетика. 1994. №1. С. 85-99.

9. Будник С.А., Гусева Л.И., Шибин А.Г. Анализ схемы измерений температуры для определения комплекса характеристик теплозащитного покрытия. //ИФЖ. 1989. Т.56. №3. С432-440.

10. Воеводин В.В. О методе регуляризации. //ЖВМ и МФ. 1969. Т.9., №3. С. 673-675

И. Воскобойников Ю.Е., Томсонс Я.И. Построение регуляризованно-го решения одной обратной задачи теплопроводности при случайных ошибках в исходных данных. // ИФЖ. 1977. Т.ЗЗ. №6. 12. Геджадзе И.Ю., Шутяев В.П. Обоснование метода возмущений для квазилинейной задачи теплопроводности. //ЖВМ и МФ. 1998. Т.38 №6. С.948-955.

13. Жуковский E.JI. Статистическая регуляризация алгебраических систем уравнений. //ЖВМ и МФ. 1972. Т.12. №1. С188-191

14. Карслоу У., Егер Д. Теплопроводность твердых тел. М.: Наука, 1964.

15. Коздоба Л.А. Методы решения нелинейных задач теплопроводности. М.: Наука, 1975.

16. Ландис Е.М. Некоторые вопросы качественной теории эллиптических и параболических уравнений. //Успехи математических наук. 1959. Т.14. вып. 1(85). С.67-85

17. Лоусон Ч., Хенсон Р. Численное решение задач метода наименьших квадратов. М., Наука. 1986.

18. Марпл С.Л. Цифровой спектральный анализ и его приложения. М.: Мир. 1990.

19. Мацевитый Ю.М., Мултановский A.B. Идентификация в задачах теплопроводности. Киев: Наукова думка. 1982.

20. Мишин В.П., Алифанов О.М. Повышение качества отработки теп-лонагруженных конструкций и обратные задачи теплообмена. 4.II. Практические приложения. //Машиноведение. 1986. №6. С. 11-21.

21. Морозов В.А. О принципе невязки при решении операторных уравнений методом регуляризации. //ЖВМ и МФ. 1968. Т.8. №8. С. 295-309.

22. Полежаев Ю.В., Юревич Ф.Б. Тепловая защита. М.: Энергия, 1976.

23. Симбирский Д.Ф. Температурная диагностика двигателей. Киев: Техника. 1976.

24. Тарасов Р.П. Численное обращение причинных операторов Винера-Хопфа и метод цифровых фильтров в обратных задачах теплопроводности. //ЖВМ и МФ. 1993. Т.ЗЗ. №11. С.1602-1625.

25. Тихонов А.Н., Арсенин В.Я. Методы решения некорректных задач. М.: Наука, 1986.

26. Тихонов А.Н., Гончарский A.B., Степанов В.В., Ягола А.Г. Регу-ляризирующие алгоритмы и априорная информация. М.: Наука, 1983.

27. Хуанг Т.С., Эклунд Дж.-О. Быстрые алгоритмы в цифровой обработке изображений. М.: Радио и связь, 1984.

28. Шафиев Р.А., Кугель М.Я. О двупараметрическом методе регуляризации L-псевдообращения и выборе параметров. //Известия АН АзССР. Серия физико-технических и математических наук. №6. 1984. С. 24-29.

29. Ягола А. Г. О выборе параметра регуляризации по обобщенному принципу невязки. //ДАН СССР. 1979. Т.245. №1. С.71-73.

30. Alifanov О.М., Budnik S.A., Doroshin V.A., Mikhailov V.V., Nenaro-komov A.V. An Experimental-computational System Based on Inverse Methods for Investigating Thermally Loaded Structure Testing: Conception, Realization and Results. Preliminary Proceedings - Vol. 2, 2nd International Conference on Inverse Problems in Engineering: Theory and Practice. 9-14 June, Le Croisic, France, 1996.

31. Fadale T.D., Nenarokomov A.V., Emery A.F. Uncertainties in Parameter Estimation: the Inverse Problem. //Int. J. Heat Mass Transfer. Vol.38. №3. pp. 511-518. 1995.

32. Hensel E.C., Hills R.G. A Space Marching Finite Difference Algorithm for the One Dimensional Inverse Conduction Heat Transfer Problem, ASME Paper No. 84-HT-48, NY, August 1984.

33. Kurzhanski А.В., Khapalov A.Yu. An Observation Theory For Distributed-Parameter Systems. // II AS A, 1990.

34. H.J. Reinhardt, Dinh Nho Hao. A Sequential Conjugate Gradient Method for the Stable Numerical Solution to Inverse Heat Conduction Problems. Материалы конференции "Обратные задачи и идентификация динамических систем". Август 1994. Санкт-Петербург.

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

1. Алифанов О.М., Артюхин Е.А., Геджадзе И.Ю. Последовательное регуляризированное решение одной обратной задачи теплопроводности. //Доклады РАН. Уравнения математической физики. Принята к публикации 20.05.97.

2. Alifanov О.М, Gejadze I.Yu. Thermal Loads Identification Technique for Materials and Structures in Real Time. //Acta Astronáutica. Vol. 41, Nos 4-10, pp. 255-265, 1997.

3. Алифанов O.M.,Геджадзе И.Ю. Об одном методе оперативной идентификации тепловых нагрузок. // ИФЖ. 1998. Т. 71. №1. С.30-40.

4. Алифанов О.М., Геджадзе И.Ю. Об одном перспективном направлении в технике тепловых измерений. //"Проблемы машиностроения", Харьков. 1998. Т.1. №1.

5. Артюхин Е.А., Геджадзе И.Ю. Об одном методе решения задачи наблюдения для нестационарного температурного поля (линейный случай). //Изв. РАН. Теория и системы управления. 1998. №4.

6. Геджадзе И.Ю., Шутяев В.П. Об одном методе решения задачи наблюдения для нестационарного температурного поля (нелинейный случай). // Изв. РАН, Теория и системы управления. 1998. №6.

7. Геджадзе И.Ю., Шутяев В.П. Обоснование метода возмущений для квазилинейной задачи теплопроводности. // ЖВМ и МФ. 1998. Т.38. №6. С.948-955.

8. Геджадзе И.Ю. Быстрый алгоритм решения двухпараметрической задачи наименьших квадратов. // ЖВМ и МФ. 1998. Т.39 №12