Применение функций чувствительности в задачах математического моделирования систем с распределенными параметрами. Часть 1

Московский государственный технический университет 
имени Н.Э. Баумана 
 
 
А.Ю. Бушуев 
 
 
Применение функций чувствительности  
в задачах математического моделирования  
систем с распределенными параметрами 
 
Часть 1 
 
Методические указания к курсовому 
и дипломному проектированию 
 
Москва 
Издательство МГТУ им. Н.Э. Баумана 
2011 

УДК 536.24 
ББК 16.4.1 
        Б90 
Рецензент С.В. Аринчев 
Бушуев А.Ю.  
Б90 
Применение функций чувствительности в задачах матема-  
тического моделирования систем с распределенными параметрами: метод. указания к курсовому и дипломному проектированию. – Ч. 1. – М. : Изд-во МГТУ им. Н.Э. Баумана, 
2011. – 44, [4] с. : ил.  
 
Рассмотрены вопросы применения функций чувствительности 
к различным задачам, возникающим в инженерной практике при 
проектировании технических систем, описываемых уравнениями в 
частных производных. 
Для студентов старших курсов, обучающихся по специальности «Прикладная математика». Методические указания могут быть 
использованы при выполнении курсовых и дипломных работ. 
Рекомендовано Учебно-методической комиссией НУК ФН 
МГТУ им. Н.Э. Баумана. 
УДК 536.24 
                                                                                                           ББК 16.4.1 
 
 МГТУ им. Н.Э. Баумана, 2011 

 
ВВЕДЕНИЕ 
Аппарат функции чувствительности применяют для анализа и 
синтеза технических систем. Анализ чувствительности основан на 
изучении влияния изменения параметров технических систем на 
их поведение. Например, при проектировании конструкции необходима информация о влиянии параметров проектирования на 
целевую функцию задачи и функции ограничений. Частные производные этих функций по проектным параметрам и называют 
функциями чувствительности (ФЧ). Плодотворность и эффективность данного подхода подтверждена разнообразным применением в технике. Так, ФЧ широко используют при оптимальном 
проектировании механических конструкций [1, 2], идентификации 
математических моделей сложных динамических систем [3], оптимальном планировании эксперимента [4]. Перечисленные направления исследований представлены в монографиях или научных статьях и сложны для первоначального знакомства 
студентов.  
В предлагаемой первой части методических указаний рассмотрено применение ФЧ в задачах моделирования систем с использованием уравнения теплопроводности. 
1. ПРИМЕНЕНИЕ ФУНКЦИЙ ЧУВСТВИТЕЛЬНОСТИ  
В ЗАДАЧЕ СИНТЕЗА МНОГОСЛОЙНОЙ 
КОНСТРУКЦИИ 
Многие энергетические установки, паровые и газовые турбины, парогенераторы, теплообменники элементы конструкций 
ракет и двигателей летательных аппаратов содержат дорогостоящие теплоизоляционные и теплозащитные материалы. Возможности повышения экономичности, надежности и ресурса 
 
3

работы таких установок в значительной мере зависят от совершенства применяемой теплозащиты, ее оптимизации. Поэтому 
актуальна проблема проектирования многослойной тепловой 
защиты конструкций по различным критериям. 
Неразрушающиеся теплозащитные покрытия (ТЗП) широко 
применяют на летательных аппаратах в том случае, когда их геометрические размеры в процессе гиперзвукового полета в плотных 
слоях атмосферы должны сохраняться. Типовые примеры летательных аппаратов такого типа – «Спейс шатл» и «Буран». 
Рассмотрим следующую проектную задачу. Требуется определить толщины слоев 


1
2
,
, ...,
m
h
h h
h

 конструкционного пакета 
(КП), подверженного воздействию высокотемпературной среды, из 
условия обеспечения равенства температур на границе слоев заданным значениям: 






1
2
1
2
ˆ
,
,...,
,
,...,
0,
i
m
i
m
I i
h h
h
T
h h
h
T
1,...,
,
0,
,
i
m
I i
n


 
(1.1) 
 







где 
i

 – функционал-рассогласование между максимальной температурой в i-м контролируемом узле 
i
T  и максимально допус




;  – 
тимым ее значением 

ˆ
,
I i
T
 К; 





ˆ
0,
,
max
,
i
i
i
T
T x
T x


время, с; 

 – момент времени достижения максимальной температуры в контролируемом узле; ˆ
 – правое значение временного 
интервала; ( )
I i  – номер границы КП, на которой действует i-е ограничение; m  – число контролируемых температур (число варьируемых толщин); n – число слоев в КП. 
Температурный режим конструкции находится из решения 
следующей краевой задачи: 


,
k k
k
T
T
c
T
T
x
x















 
1
,
k
k
x
x x

 
ˆ
0 , 
1,...,
k
n

, (1.2) 
 


0
0
,0
,
,
n
T x
T
x
x
x



 
(1.3) 
 




0,
0,
,
1,...,
1,
0
k
k
T x
T x
k
n







, 
(1.4) 
 
4 

0,
0,
,
k
k
k
k
T x
T x
T
T
x
x











 






 
(1.5) 
k
n
1,...,
1,
0,



где x – координата, м; 
k
с  – удельная теплоемкость k-го слоя КП, 
Дж/(кг·К); 
k
 – коэффициент теплопроводности k-го слоя КП, 
Вт/(м·К); 
k
 – плотность материала k-го слоя КП, кг/м3; 
0
T  – начальная температура конструкции, К. 
Граничные условия на поверхностях 
0
w  и 
n
w  КП в общем 
случае имеют вид 
 




0
0
1
,
w
T x
T
q
T
x





, 
(1.6) 
 




,
,
n
n
n
w
T x
T
q
T
x





 
(1.7) 
где 
0
w
q
, 
n
w
q
 – плотность тепловых потоков, подводимых к границам 
0
w  и 
n
w  КП, Вт/м2. 
Решение поставленной задачи находят в рамках двухконтурного итерационного алгоритма. При этом задача подбора толщин 
решается во внутреннем контуре, где применяется приближенная 
математическая модель, отличающаяся от (1.2) – (1.7) использованием фиксированных значений теплофизических характеристик 
для каждого слоя пакета. Кроме того, решение задачи прогрева 
здесь требуется найти для существенно меньшего числа узловых 
значений искомых функций. Корректировка этой математической 
модели осуществляется во внешнем контуре на базе решения исходной постановки задачи (1.2) – (1.7) при толщине слоев, найденных во внутреннем контуре, и должна обеспечивать идентичность температур в узлах приближенной математической модели, 
найденных при одинаковой толщине слоев пакета в исходной и 
упрощенной постановках. 
Указанная проблема может решаться разными способами. Так, 
в статье [5] она проводилась путем использования во внутреннем 
 
5