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

УДК 678.7

РАЗРАБОТКА ОБОБЩЕННЫХ МАТЕМАТИЧЕСКИХ МОДЕЛЕЙ ПОЛИМЕРНЫХ КОМПОЗИЦИОННЫХ МАТЕРИАЛОВ, ПРИМЕНЯЕМЫХ В АВИАСТРОЕНИИ

© 2019 Г.В. Дмитриенко, Е.Н. Згуральская, Г.Л. Ривин, А.А. Федоров

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

г. Ульяновск

Статья поступила в редакцию 22.08.2019

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

Ключевые слова: математическая модель, полимерные композиционные материалы, авиастроение.

Исследование выполнено при финансовой поддержке РФФИ и Правительства Ульяновской области в рамках научного проекта № 18-48-732005\18

ВВЕДЕНИЕ

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

Дмитриенко Герман Вячеславович, доктор технических наук, профессор кафедры «Самолетостроение». E-mail: dmitrienko.german@yandex.ru Згуральская Екатерина Николаевна, старший преподаватель кафедры «Самолетостроение». Георгий Леонидович Ривин, кандидат технических наук, доцент кафедры «Самолетостроение». E-mail: avia@ulstu.ru

Федоров Александр Александрович кандидат технических наук, доцент кафедры «Самолетостроение».

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

ОПИСАНИЕ ЗАДАЧИ

Проведя всесторонний системный анализ, установлено:

- ПКМ характеризуются следующими параметрами: e - диэлектрическая проницаемость материла, tgd - тангенс угла диэлектрических потерь, d - проводимость (поверхностная проводимость);

- ПКМ имеют различные температурные характеристики.

Исходя из перечисленного ставиться следующая задача: составление более точных математических моделей: электродинамических и тепловых моделей для полного охвата ПКМ используемого на производстве, для проведение диагностики их характеристик в свободном пространстве.

ПУТИ РЕШЕНИЯ ЗАДАЧИ

Для составления математических моделей электродинамики применяется более строгий с электродинамической точки зрения математический метод, метода интегральных уравнений (МИУ). Данный метод активно развивался и применяется для радиоволновых измерений в свободном пространстве [9-21], также МИУ применим для различных нетрадиционных задач. Так, Васильевым Е.Н. МИУ - развивался в антенной технике, для диэлектрических и проводящих тел имеющих осевую симметрию [17]. Кислюком М.Ж. МИУ - использовался в методе вторичных волн для решения граничных задач в свободном пространстве, в волноводах и резонаторах [18]. С помощью МИУ были получены точные аналитические решения многих практических задач: дифракции, электродинамики и антенной техники, которые раньше решались приближенными методами и не имели точных аналитических решений.

Составление ИУ производится из уравнений Максвелла, используя специальные функции: вектора Герца или векторные потенциалы. Эти функции должны удовлетворять условию излучения на бесконечность, которое математиче-

Е (р (2 АЕ + graddivAE)- гоАм;

(1)

ски записывается как: Ит г

ду дг

+ к и = о. миу

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

ЭЛЕКТРОДИНАМИЧЕСКИЕ МОДЕЛИ ПКМ

Процедура составления ИУ производиться через поверхностные токи, наводимые на поверхности образца ПКМ зондирующей электромагнитной волной. В ИУ поверхностные токи задаются через векторные электрический или магнитный потенциалы. Диэлектрические параметры ПКМ и характеристики отражающей поверхности задают через поверхностный импеданс. Информационным параметром является отраженная волна, [21-25]. Результаты измерений обрабатываются численными методами в строгой электродинамической постановке.

Для сокращения размерности задачи применяется методика составления ИУ для тел вращения. Напряженности электрического и магнитного полей выражается через возбуждающие токи.

Н (р) = (к2 Ам + graddivAM) + гогАЕ,

где к - волновое число; £а - абсолютная диэлектрическая проницаемость; М'а - абсолютна магнитная проницаемость; м> -круговая частота.

АЕ =|3Е (д)}{р, д)1оч,

Ам =|3м (д)о(р, д, (2)

х

где АЕ, Ам - векторные электрические и магнитные потенциалы; ¡Е(ц), ]м(ц) - электрический и магнитный токи d точке истока; 0(р,ц)

^ Ехр(- №г)

- функция Грина О =----, г - рассто-

4яг

яние между точками наблюдения р и истока ц; р - точка наблюдения; ц - точка истока; 5 - граница поверхности Б материала или поверхности датчика.

Эти соотношения удовлетворяют волновому уравнению Гельмгольца и граничным условиям излучения на бесконечность [17,21]. Сделав подстановку (2) в (1) получим

Е(р) = И-]Щ1аЗЕ (др(р, д) + — (др(р, д))-

1П1Р.

- тсг(.1и Ш(р, ч№оя ;

(3)

Н (р) = г! - (д)в(р, д) + — (д р(р, д)) +

I I }—■«

реБ; О--

+ го( (д Р(р, д))) . (4)

Ехр(- ); К = . 4жг ' ^ а

Полученные уравнения эквивалентны ИУ для диэлектрических тел, в которых электрический ток протекает по поверхности или в скин-слое проводящих тел. Образец ПКМ по своим электродинамическим свойствам относится к импедансным материалам, имеющим большое затухание. При этом можно считать, что токи наводимые внутри материала ПКМ имеют малую амплитуду и создают малую составляющую поля во внешнем пространстве, поэтому ими можно пренебречь. На основании изложенного выше интегрирование производится по поверхности, а не по объему [17]. В связи с громоздкостью преобразований промежуточные выкладки опускаем.

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

1. Считаем, что электромагнитная волна, проникающая в глубь ПКМ перпендикулярна внешней поверхности, и направлена по одной из координатных осей (ось 2).

2. В силу правила перестановочной двойственности уравнений Максвелла, следующих из их симметрии, тензорные функции хЕМ, НЕМ связаны между собой следующими выражениями.

ЬЕ (р, д,£а, ра )=-^М (р, д,еа, ра)

( Ч,8а , Ра )=£Е ( Ч,8а , Ра ) .

Тогда используя утверждение, что вектора Е и Н удовлетворяют уравнению Гельмгольца, можно сделать переход от трехкратного интеграла по поверхности Б к двукратному интегралу по поверхности Ь.

Е (р ) = ¡£Е (р, д,еа \1Е (д) +

ь

+ Г (р, д,еа, Ма)м (д), (8)

Н (р ) = | ЪЕ (р, д,еа, \1Е (д) +

ь

+ Нм ( д,еа, Ма )м (д)1. (9)

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

В зависимости от параметров е, tgd и d все типы ПКМ можно условно разделить на группы материалов, соответственно три основных модели для проведения диагностики:

1. полубесконечного тела;

2.тонкого тела;

3. многослойного тела.

Модель полубесконечного тела

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

представим описание для волноводного датчика. Обозначим Ь1 - внешнюю границу поверхности волноводной конструкции, Ь2 - внешнюю границу поверхности измеряемого образца ПКМ. Тогда электромагнитное поле внутри области 1 измерительного датчика в точке наблюдения р согласно (5) и (6) представляется в виде

Е(р)= | (я >11 +

" (10) + \ [[ (я ) + (я)].

Н(р )= \ Г Iм (я )1 +

г " т (11)

+ \ [ (я)+ъм7м (я)].

где Ь1 - граница волноводной конструкции, Ь2 -граница измеряемого образца ПКМ.

Используя истокообразное представление (10) и применим граничное условие на поверхности образца Ь2 (р1Ь2) в виде, эквивалентном виде

~м =-3Е ], I = ^VРо /, (12)

получим интегральное уравнение относительно тока ]Е первой гармоники. Полное поле внутри образца ПКМ на большом удалении от границы раздела ПКМ и свободное пространство равно нулю.

Рис.1. Расположение вспомогательной поверхности

Это условие будет выполняется, если на вспомогательной поверхности (рис. 1), расположенной внутри тела образца ПКМ или волно-водного датчика, обращается в ноль тангенциальная составляющая полного магнитного поля в выражении (11). Для получения интегрального уравнения точка наблюдения р помещается на вспомогательную поверхность Бф, расположенную внутри образца ПКМ, умножаем истокообразное уравнение на нормаль в точке наблюдения р, направление нормали в сторону источника излучения, а правую часть приравнивается к нулю, получаем:

О = п х \ ^ (я )1 +

и (13)

+ Й | Ь1ЕТ1Е (Я)-(я ЦЕ (я )}11,

Ь2

где 2 =

^ О О 2,

- матрица поверхностного им-

педанса.

ИУ на границе Ь1 относительно магнитного тока получается из истокообразного представления (8). Опуская точку наблюдения р на границу Ь1 датчика и умножая на нормаль в точке р, получается следующее интегральное уравнение

(р)+1 (я >11

+

и

+1 ет (я )ЧГ2[п (я ЦЕ (я)] = о.

(14)

Ь2

В выражении (14) наличие внеинтегрально-го члена представляет поле падающей волны на границе Ь1 в задаче на отражение. Полученная система из двух интегральных уравнений (13) и (14) описывает электромагнитное поле в измерительном волноводном датчике с образцом ПКМ. Система является полной, т.к. неизвестным в ней являются только ток 3^ на границе Ь2. Полученные уравнения записаны в векторной форме, делая переход к скалярной форме, уравнения (13) и (14) решаются численно, вычисляется неизвестный ток 3^ на границе Ь2. Составление интегральных уравнений для конусного рупора, такое же отличие заключается в поверхности интегрирования, поверхность волновода заменяется на поверхность рупора. Составленные интегральные уравнения на поверхности Е в результате позволяют получить распределение реального существующего тока. Составленные интегральные уравнения на поверхности Е' используются как вспомогательные, и применяется когда возникают трудности с неустойчивостью решения.

Модель электрически тонкого тела

Структура измерительного радиоволнового датчика на базе круглого волновода, МИУ позволяет численно рассчитать комплексный коэффициент отражения для рупора произвольной

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

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

Составляется система интегральных уравнений для внешней (освещенной) и тыльной (теневой) границ тонкого тела, используя уравнения (8) и (9). На внешней поверхности (11) поле зондирующей волны создается током 3м, на тыльной поверхности (12), прошедшая волна будет описываться наведенными токами и 32. Процесс составления интегральных уравнений такой же, как для полубесконечной модели, учитывая на границах Ь1 и Ь2 граничные условия для импедансного тела. Запишем систему уравнений:

-(Р1)+ | (* >11 +

1 ы

+ | И (Ч1 )-^[й (Я1ЦЕ (Я1)] = 0;

Ь2

+1 № (2 )1+

ы

+1 (() (*2)-(ЧГ (<ь Ь^ (<ь )] = о;

(15)

(16)

У

Ни , 1 Ъ

1 _ 1 2

а)

б)

Рис. 2. Конструкция волноводного датчика для измерений на прохождение: а) - на круглом волноводе; б) - на круглом волноводе с коническим рупором

^(р2) +1Ь^ (Я2 >11 + + I( -Ь^)(я2)-( ->^¡5(я2>3*(я2)}11 = о.Ч

(17)

где р1, точки на границе Ь1, а р2 ц2 на границе 12.

Уравнение (15) для внешней поверхности Ь1 аналогично уравнениям для однородного полубесконечного тела. Последний интегральный член в левой части его учитывает отражение от тыльной поверхности. Уравнения (16), (17) описывают поле прошедшей волны, при этом первые интегральные члены учитывают влияние внешней поверхности. На обоих поверхностях выполняются граничные условия (12). Полученные уравнения представляют собой уравнения Фредгольма второго рода, которые решаются численными методами. Описываемый ток, в интегральных уравнениях на освещенной поверхности представляется как = JМад + JМтр, на

затененной поверхности 3 м = Jм .

По предложенной математической модели был проведен контроль образца ПКМ. По результатам измерений было определено значение относительной диэлектрической проницаемости материала У-ПУ е=940±30% и тангенса угла диэлектрических потерь tg5=1,3±40%, имеющих модуль коэффициента отражения 0,94.

Многослойные модели ПКМ

В предыдущих моделях считали, что поле, проникающее вглубь измеряемого образца ПКМ мало, и его воздействием можно пренебречь в случае, когда закон распределения токов внутри ПКМ не влияет на закон распределения токов на поверхности ПКМ. Эти модели (бесконечной пластины) пригодны для материалов, имеющих большое затухание. В случае, когда материалы имеют высокое значение коэффициента отражения и малую геометрическую толщину, а токи имеющиеся внутри ПКМ не малые, составляется многослойная модель для ПКМ как в [17,21], в которой протекающие внутри токи влияют на токи, протекающие на поверхности ПКМ. Такая модель применима для материалов имеющих сложную внутреннюю структуру и большое значение tgd близкое к 1 и более. В случае на отражение уточненная многослойная математическая модель имеет вид рис.3, на прохождение рис. 4.

Модель образца ПКМ представляется виде многослойной структуры состоящей из однородных слоев, количество слоев определяется дискретом разбиения по е . Интегральные уравнения для многослойной структуры составляется на основе уравнений полученных ранее (8) и (9) для каждого слоя, используя обе его границы. Нумерация слоев производится в соответствии с

рис. 3 и рис. 4. Методика составления интегральных уравнений для многослойной структуры следующая. Составление интегральных уравнений для многослойной среды производиться как для одного слоя, учитывая распределение токов на границах раздела. Сшивка производится по граничным условиям Щукина-Леонтовича т.к. ПКМ обладает е , ц - электрической и магнитной проницаемостями и ст - проводимостью. Образец разбит на N - слоев, их границы не имеют общих точек. В каждом слое могут находиться сторонние источники, нумерация слоев производится по направлению от источника наблюдения в глубь образца ПКМ. Составление уравнений производится, для слоев К и К+1 имеющих общую границу БК используя (8), (9) к К-му слою применяют теорему эквивалентности. Нормаль на границах поверхностей раздела слоев направлена в сторону первоисточников.

Для К слоя уравнения имеют вид:

+

+

~ЕК (р) = | [ )+$ м!км )]}

^ У

Ук

\ [[ (я ) + (я)} + Е* (р)

УК-1

реБк

±нк (р )= \ [ (я)+ьм!км (я)}

+

+

\ [ (я) + (я)]} + Н* (р)

реБк (18)

Интегральные уравнения составлены в матричной форме для слоя К, аналогично записываются для слоя К+1, интегралы берутся по поверхностям БК и БК+1 как (18).

^Ек+1(р)= \ [[ Н^З+ЛяК

^ У

УК+1

+ \ [ (Я )+№ (я)} + Е^+1 (р)

Ук

РеБк+1

:^Нк+1(р)= \ [(я)+(я)}

УК+1

\ [ (я)+(я К+н^Ы

+

+

+

(19)

Для получения возбуждаемых токов на поверхности уравнение (18) домножаются вектор-но на нормаль и на ек+1 , а (19) - на нормаль и ек , затем уравнения (18) и (19) складываются. Получаются уравнения токов, протекающих на границе раздела слоев.

У

К-1

П1

п2

П1

п1+1

У Структура НКМ

П1+2

плоскость прилегания вспомогательные поверхности в НКМ

- образец НКМ

Рис. 3. Многослойные модели ПКМ на отражение: а) - математическая; б) - физическая модель

П1

п2

П;

П;+1

П;+2

/Структура НКМ

П|+1

П)+2

"У

ОООООСа

плоскость прилегания вспомогательные поверхности в НКМ

> образец НКМ

Рис. 4. Многослойные модели ПКМ на прохождение: а) - математическая; б) - физическая модель

2е

[.Е? (р)]+"%^ [п (р).^ (р)]+

?кМ (Р ) =

е К + е К+1 е К + еК+1

+ п(р)х | {к%Е - ек+,%Е(я) + [ек- ек+]^м ]ГКМ (я)] +

ек

+ п(р)х |ек%^ (я)+ек%м!кМ-1 (яК - п(Р)х

Ек-,

х |ек+1%Е^к+1 (Я)+ек+1%м!км+1 (я)<4

Ек+1

Р^Ек

7кЕ (р)=-^— [[ (р [н^ (р)]+

- К +- К+1 - К +- К+1

+ п (р)х | { киЕ -- к+у ] (я)+[- [ -- к+1им ] (я )К + + п(р)х |-кьЕ^-1 (я)+-кьм7км-1 (яМа, - п^)х

х | -к+.ьЕ (я) + -к^"7КМ+: (яК

Ек+1

Р^Ек

В описании процесса на отражение можно заменить внутреннюю структуру одним слоем, и в описании на прохождение также можно заменить на один слой конечной толщины. Составление ИУ производится для каждого слоя, для границ и Е№1 интегральные уравнения

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

^ (р)= 2[п(р),Н^т (р)]+

+ 2п(р)х | ЪЕД1Е (я)ёач +2п(р)х

е0

х | [ (я) + ЬМГ1М (я)К , реЕ0

(р)=

2е. е, + е,

[)ДСТ (р)]+

2е,

[р),Ё2ст (р)]н

^х|{[ -е^(я)+ [е2%м -е1%м]Г1М(я)}

+ е 2 Е,

п(р)х | %Е:Е (я)аач

е + е J ч

Т С2 Е0

РеЕ1

+

+

+ J{hE -_,hE] (q)+[ -_,hM](q)]dc

_ + _ i

+ _?__п(р)х JhEJE(q)d, _ + _2 i

o„

Р

МИУ применим к для измеряемых образцов ПКМ с гладкой рабочей поверхностью (по отношению к размеру зондирующей волны), с математической точки зрения точен.

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

Алгоритм решения интегральных уравнений

Решение интегральных уравнений производится численными методами. Наиболее распространенным для интегральных уравнений является проекционный численный метод (метод Галеркина) [18-20]. Решение находится в виде

N

суммы: и(у) = ^(у), которая подставля-

п=1

ется в интегральное уравнение. Функция %п(у) является начальной функцией некоторой полной системы. Далее выбирается другая система функций ц.(х) - (1=1, 2, ..., N , на каждую из них домножаем интегральное уравнение и интегрируем в пределах (с, d). В результате получается система линейных алгебраических уравнений (СЛАУ):

N

Xап0гп = Ь, г = 1, 2, ..., N, (20)

где

d b

Ягп = (ХУЬ (^ (у\lxdy ,

с а

а

ь=| / (*к- (х ^ .

с

Метод хорош если функции %п(у) при небольшом N хорошо аппроксимируют решение, а интегралы (22) считаются аналитически. Функции и -л. могут быть одинаковыми.

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

составляющие Бг , Е электромагнитного поля в волноводе. В круглом волноводе поле основной волны И11 представляется как

.2 Г .. \

Er — J—-J1

_11

_

V r0

Sin(()Exp(- Jhz) .

E( = J~-J1

_

_

(29)

V ro

Cos((p)Exp(- jhz).

где w - угловая частота колебаний; pia - абсолютная магнитная проницаемость; цп - корень функции Бесселя; J , J[ - функция Бесселя первого порядка и ее производная; r0 - радиус волновода; h - волновое число; r, z, j - составляющие цилиндрической системы координат. Операторы rot, grad, div интегральных уравнений вида:

H(p)— | j- jmeaJM (q)G(p, q) + ■ graddiv (jM (q )G(p, q))+

+

rot(JE (q )G(p, q ))

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

Е7д = £ Е7 ехр(тр),

^ 2л

где ЕГ = — | ЕГ ехр(- тф^ф.

2я 0

Количество членов обычно берется 5-7, для получения устойчивого решения. В измерительном датчике распространяется только волна И11 , поэтому делается преобразование для первой гармоники, получается

E —

1 ^Л

2 r_

Л

f \ _ r

V r0

E( — j V_aro J/

2 _

_

\

V ro

(30)

В элементах для тензорной функции Грина все расстояния нормируются, умножаются на к - волновое число свободного пространства, для приведения к одной размерности в (2.32).

El =■

El —

К

Jl

J'

V r0

_1

V Г

Делая подстановку, переходим к итоговым

r

r

n—1

r

2

r

0

r

r

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

Составленные системы уравнений позволяют производить расчеты, когда в качестве информационного параметра использования комплексный коэффициент отражения, для этой цели задаваемый магнитный ток представляется как сумма токов падающего и отраженного 3м = 3Мад + 3Мтр. В этом случае в систему добавляется неизвестный отраженный ток, избыточность составленной системы позволяет его вычислить. Решение задачи представлено для симметричного возбуждения, возбуждение производится волной основного типа H11. В итоговых выражениях две составляющие диагональной матрицы поверхностного импеданса с помощью (2.1) выражается через составляющие полей: 11 = Ек/Иу , 12 =- Еу/Ик . При нулевом импедансе получаем выражения для идеально проводящей поверхности.

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

ЗАКЛЮЧЕНИЕ

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

По приведенной теории был произведен контроль и определено значение относительной диэлектрической проницаемости образца ПКМ, материала У-ПУ е=950±30% и тангенса угла диэлектрических потерь tgd=1,3±40%, имеющих модуль коэффициента отражения 0,94.

В результате проведенного рассмотрения процесса диагностики e ПКМ МИУ, производит усложнение способа обработки результатов контроля, выигрыш за счет усложнения математического аппарата в случае с ПКМ. Усложнение открывает более широкие возможности по автоматизации и обработке результатов измерений, использую высокопроизводительную компьютерную технику.

СПИСОК ЛИТЕРАТУРЫ

1. Буланов, И.М. Технология ракетных и аэрокосмических конструкций из композиционных материалов. / И.М. Буланов, В.В. Воробей. - М.: Изд-во МГТУ им. Баумана, 1998. - 514 с.

2. Тарнопльский, Ю.М. Пространственно - армированные композиционные материалы. Справочник. / Ю.М. Тарнопльский, И.Г. Жигун, В.А. Поляков. - М.: Машиностроение, 1987. - 224 с.

3. Волокнистые и дисперсноупрочные композиционные материалы. Справочник. / В.В. Васильев. -М.: Изд-во Наука, 1976. - 214 с.

4. Углеродные волокна и углекомпозиты. // Под ред. А.А. Берлина. - М.: Мир, 1988. - 270 с.

5. Углеродные волокна и углекомпозиты. / Под ред. Э. Фатцер. - М.: Мир, 1988. - 336 с.

6. Емец Ю.П. Электрические характеристики композиционных материалов с регулярной структурой./ Ю.П. Емец. - Киев: Наук. думка, 1986. -191 с.

7. Макалистер, Л. Многонаправленные углерод-углеродные композиты. / Л. Макалистер, У. Лак-ман // Прикладная механика композитов. - М.: Изд-во Мир, 1989. - С. 226-294.

8. Заявка 2772520 Франция, МПК6Н 01 О 17/00 Material composite structural absorbant les ondes radar et utilisation d>un tel materiau: / Escarmant Jean Francois; Giat Ind. SA.- /9715681; Заявл. 11.12.97; Опубл. 18.6.99.

9. Шуньлянь, Ч. Рассеяние волн поверхностями с периодической структурой./ Ч. Шуньлянь, Г. Жэньао // ТИИЭР. - 1981. - Т.69. - №9. - С. 43-56.

10. Whitman G., Schwering F. Scatterring by periodic metal surfaces approach. "IEEE MTT-s Int. Microwave Symp., Cherry Hill, N.S., 1976, Dig.Techn.Pap." New York, N/Y., 1976, pp.362-364(англ)

11. Wu Te-Kao, TsaiLeonardL. Scatteringfrom arbitrarily-shaped lossy dielectric bodies of revolution. "Radio Sci" 1977, t.12, №5, pp.7098-718(англ)

12. Matsuhara Masanori "Дэнси дзёхо цусин гаккай ромбунси, Trans. Inst. Electron., Int and Commun. Eng. Jap.", 1987, -C70,№4, pp.487-492 (яп)

13. Кулижский, А.В. Распространение сигнала в двухслойной среде со случайными крупномасштабными неоднородностями./ А.В. Кулижский, М.В. Тинин // Изв. ВУЗов. Сер. Радиофизика. - 1993. -Т.36. - №6. - С.475-483.

14. Комаров, С.А. Рассеяние радиоволн на плоскослоистом полупространстве с шероховатой границей / С.А. Комаров, А.И. Якушев // Радиотехника и электроника. - 1998. - Т.43. - №6. - С. 650-656.

15. Дандаров, В.А. Численное решение 2-х мерного интегрального уравнения Фредгольма в задаче о распространении радиоволн над нерегулярной земной поверхностью / В.А. Дандаров //Радиотехника и электроника. - 2002. - Т.47. - №11. -С.1329-1334

16. Benali A., Chandezon J., Fontaine J. A new theory for scattering of electromagnetic waves from conducting or dielectric rough surfaces //IEEE transactions on antennas and propagation, 1992, vol.40, no.2. pp. 141-148

17. Васильев, Е.Н. Возбуждение тел вращения / Е.Н. Васильев- М.: Радио и связь, 1988. - 271 с.

18. Кислюк, М.Ж Метод вторичных волн в задачах электродинамики / М.Ж. Кислюк. - Л.: 1970. - 151 с.

19. Малушков, Г.Д. Методы решения задач электромагнитного возбуждения тел вращения./ Г.Д. Малушков // Изв. ВУЗов. Сер. Радиофизика. - 1975. - Т.18. - №11. - С. 1563-1587.

20. Никольский, В.В. Проекционные методы в элек-

тродинамике // Прикладная электродинамика./ В.В. Никольский - М.: Высшая школа. - 1977. -вып.1. - С.4-50.

21. Хижняк, Н.А. Интегральные уравнения макроскопической электродинамики / Н.А. Хижняк. -Киев: Наук. думка, 1986. - 256 с.

22. Неганов, В.А. Линейная макроскопическая электродинамика / В.А. Неганов, С.Б. Раевский, Г.П. Яровой. - М.: Радио и связь, т. 1, 2000. - 512 с.

23. Неганов, В.А. Электродинамика и распространение радиоволн / В.А. Неганов, О.В. Осипов, С.Б. Раевский, Г.П. Яровой. - М.: Радио и связь, 2005. - 648 с.

24. Петров, Б.М. Электродинамика и распространение радиоволн / Б.М. Петров. - М.: Горячая линия-Телеком, 2003. - 558 с.

25. Никольский, В.В. Электродинамика и распространение радиоволн / В.В. Никольский, Т.И. Никольская. - М.: Наука, 1989. - 452 с.

DEVELOPMENT OF GENERALIZED MATHEMATICAL MODELS OF POLYMERIC COMPOSITE MATERIALS USED IN THE AIRCRAFT

© 2019 G.V. Dmitrienko, E.N. Zguralskaya, G.L. Rivin, A.A. Fedorov

Ulyanovsk State Technical University, Separate Structural Unit "Institute of Aviation Technologies and Management",

Ulyanovsk

The article discusses mathematical models for diagnosis. polymer composite materials by radio wave methods.

Keywords: mathematical model, polymer composite materials, aircraft construction.

German Dmitrienko, Doctor of Engineering Sciences, Associate Professor, Professor at the Department of Aircraft Engineering E-mail: dmitrienko.german@yandex.ru Ekaterina Zguralskaya, Senior Lecturer at the Department of Aircraft Engineering.

Georgy Rivin, Candidate of Engineering Sciences, Associate Professor at the Department of Aircraft Engineering. E-mail: avia@ulstu.ru

Alexander Fedorov, Candidate ofTechnical Sciences, Associate Professor at the Department of Aircraft Engineering.