Явные модели распространения изгибных краевых волн в тонких полубесконечных ортотропных пластинах Текст научной статьи по специальности «Математика»

УДК 534.121.1

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

Е. Л. Коссович

Саратовский государственный университет E-mail: [email protected]

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

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

Explicit Models for Flexural Edge Waves in Thin Orthotropic Plates

E. L. Kossovich

Analysis of flexural edge wave propagation in thin plates is presented. Several problems of semi-infinite plates vibrations are solved. These plates are assumed to be orthotropic. Some basic features of flexural edge wave propagation are found using the constructed explicit parabolic-elliptic models. They extract the localized wave contribution into the overall solution.

Key words: orthotropy, deflection, shear force, bending moment, explicit dual parabolic-elliptic model.

ВВЕДЕНИЕ

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

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

Рассмотрим колебания свободного торца полубесконечной ортотропной пластины, конфигурация которой представлена на рис. 1.

Будем считать, что основные направления ортотропных свойств пластины совпадают с осями 0х, 0у. Тогда уравнение, описывающее прогиб — такой пластины, имеет вид

д4 w

Dx дХ* +2(Dl +2Dxy)

д4 w d4w d2w

xy; dxw + Dy дУ4 +2phд2 =

(1)

у* ^

?h ' _____ Х1

/ , « —

Рис. 1. Полубесконечная пластина, декартова система координат

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

Бх

.ху > 0, Бх + Бу > 0, БхБу - Б1 > 0.

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

© Коссович Е. Л., 2013

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

Для свободного от нагрузок края имеем:

д2т ^ д2т ^ д3т „ ^ , д3т

Бу—г + =0, Бу—- + (Б1 +4 БХу )7Г^- =0,

ду2 дх2 ду3 дх2ду

для случая приложенного к торцу изгибающего момента краевые условия имеют вид ^ д2т ^ д2т ,, , ^ д3т „ ^ , д3т

БУ + Б 6X2 = -м"(х'*:), Б дуз + (Б +4Бху) дХ2дУу =0 (2)

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

^ д2т ^ д2т ^ д3т , ^ , д3т

Бу—г + Бг—г = 0, Бу— + (Бг + 4Бху) = -Щ,

ду2 дх2 ду3 дх2 ду

где М" — изгибающий момент, приложенный к краю пластины, а — обобщенная перерезывающая сила Кирхгофа.

Для ортотропной пластины дисперсионное уравнение, решением которого является скорость краевой изгибной волны, получается способом, аналогичным описанному в [1]. Оно имеет вид

(Б х Б у - Б 2) - Б у с4 + 4Б ху V Б х Б у - Б у с4 = 0, (3)

где с = V-гвХ"/р, р — параметр преобразования Фурье, примененного к пространственной координате у в уравнениях исследуемой задачи, 5 — частотный параметр из преобразования Лапласа по временной координате,

А" =

=

л" = ТЮ1

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

При получении уравнения (3) вводились следующие безразмерные параметры и переменные:

х = у = Нц, Ь = Тт, т = НШ *,

Бх = ББх, Бу = ББ у, Бг = ББ ь Бху = ББху.

Решения уравнения (3) имеют вид

1/4 , __х 1/4

ск = (Б хБ у - (^ Б 2 + 4 Б 2у - 2Б ху , С2 = (Я хБ у - Ц Б 2 + 4Б 2у + 2Б ху . (4)

Здесь ck соответствует коэффициенту скорости изгибной краевой волны для случая ортотропной пластины, а с2 — второй корень уравнения, не имеющий физического смысла. Из равенств (4) видно, что коэффициент скорости изгибной краевой волны зависит от четырех жесткостей Бх, Бу, Б1 и Бху.

Исследуем закономерности возникновения и распространения краевой изгибной волны в случае, когда свободный торец пластины загружен изгибающим моментом. В этом случае решается задача (1), (2). Решение для образа прогиба Ш(п) имеет следующий вид:

' -гзАпл/О

О Л/

ММ" с2 А1(с)в -г-+ Д2 (с)е-1-г-вV

с4 _ с4

-гвА"\/Бу с Ск

=--АЦ^-х (с) (5)

где М"(р, в) — образ безразмерного изгибающего момента, приложенного к краю пластины,

1

ас = —¡= \ Б1 + 2Б ху - V (Б 1 + 2Б ху )2 - Б хБ у + с4, (6)

зАпл / и

вс =

Б 1 + 2Бху + V (Я 1 + 2Бху)2 - БхБу + с4,

(7)

Д1 (с) = -

2£ ху -\/ (Б 1 + 2£ ху )2 - Я х Я у + с4

вс,

Д2 (С) =

2£ ху + V 1 + 2£ ху )2 - я х я у + С4

Д(с) = (вс - «с)(с4 - с|),

ас,

X(с) = (Дх^у - Б2) - с4 - )х)у - с4.

(8)

(9) (10) (11)

В случае приложенной к торцу тонкой ортотропной пластины перерезывающей силы решение задачи строится для угла поворота вертикального сечения пластины у(х,у,Ь) = дт/ду относительно оси 0х.

Общее решение задачи для образа угла поворота У"(п) записывается в виде

)(п) =

Мо

Д1(с)е"

сп + Д2(с)е-

-гвХс

с4 с4

Д(с)

-вс П

-X (с),

(12)

где

Дз (с) = -Д4(с) =

2Б ху ^(Б1 + 2Б ху )2 - £ х Б у + с4

ас,

2Б ху - V (Б1 + 2Б ху )2 - £ х Б у + с4

вс

Вклад краевой изгибной волны обусловлен вкладом полюсов знаменателя с = ±ск в общее решение задачи о прогибе. Путем асимптотического анализа решений (5)-(12) можно построить следующие модели, явно описывающие основные закономерности распространения краевой волны. Для каждого из случаев краевого нагружения модель будет состоять из двух уравнений: параболического уравнения, описывающего распространение волны вдоль края пластины, и эллиптического, характеризующего ее затухание вглубь.

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

ск4 д4 те 2рН д2 т

+

Яу дх4 Бу д12

=ое1Г

д2Мо Бу дх2 '

где те(х,1) — краевой прогиб, связанный с краевой изгибной волной Коненкова,

д« * = _Д1 (ск ) + Д2 (ск)

Д(ск)

ск = {Пх Бу - (А/Щ +4!2у - 2!ху)

1/4

Д1 (ск ) = -

2Бху -а/ (Б1 + 2Бху)2 - БхБу + с%4

= (БхОу - иБ2 +4Б2у + 2Бху)

вс(ск),

1/4

Д2 (ск ) =

ас(ск) =

2Бху + у (Б1 + 2Бху )2 - Бх Бу + с%4

Д(ск ) = (вс (ск) - ас (ск ))(ск 4 - с*4),

ас (ск),

у/Ш

Б1 + 2Бху - л/(Б1 + 2Бху )2 - Бх Бу + с%4,

вс(ск) =

—7= \/ Б1

л/Б^

Б1 + 2Бху ^/(Б1 + 2Бху )2 - Бх Бу + ск

1

Лпл / И>

Апл / И>

о л/

2

а

с

2

2

с

2

1

1

Для области вне торца соответствующая модель описывается следующим уравнением:

д4 wi

d4win

д4 win

Dx - тН ^ + 2(D 1 + 2Dxy)+ Dy--^ =Q,

Dy ) dx4

dy4

где win(x, y) — прогиб внутри пластины, вызванный краевой изгибной волной Коненкова, с краевыми условиями на торце

f г \ п д2Win п д2we

Win(x,0) = We(x), Dy ду2 = -Di-

дх2

В случае приложенной к краю перерезывающей силы уравнение, описывающее поведение на торце пластины, имеет вид

c%4 д4ve 2ph д2v,

+

Dy дx4 Dy д12

=qf

д2 No

Dy "ox2'

/ n(2) * Д3 (c*k ) + A4 (cl)

где ve (x, t) — угол поворота на торце пластины и Qe =- ^-—

A(ck)

Аз (cl) = -

А4 (cl) =

2Dxy + V(Di + 2Dxy )2 - Dx Dy + ci

2Dxy - (Di + 2Dxy )2 - Dx Dy + ck

«c (cl),

ec(ci),

Для области вне торца пластины соответствующее уравнение записывается следующим образом:

Dx

д\п Dy) дx4

+ 2(Di +2Dxy)

д4 Vi

дx2дyy

+ D

д4 Vi,

y ду4

= 0,

где — угол поворота, вызыванный вкладом краевой изгибной волны внутри пластины. Соот-

ветствующие краевые условия на торце имеют вид

Vin (x, 0) = Ve (x),

д 2?

Dy^ = -(Di +4Dxy^

д 2? 'xy) дx2

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

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

Рис. 2. Трехмерные профили краевой изгибной волны Коненкова в ортотропной пластине (связанная параболическая эллиптическая модель): а — Г 1 = 0.99

4

*

а

0.01 0.008 0.006 0.004 0.002 0

-0.002 -0.004 -0.006 -0.008 -0.0^

о .40 -30

б

Рис. 2. Окончание. б — ВХу = 5

О _40 -30

б

Рис. 3. Трехмерные профили полного прогиба ортотропной пластины (точное решение): а — В1 = 0.99; б — Вху = 5

ЗАКЛЮЧЕНИЕ

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

Работа выполнена при финансовой поддержке РФФИ (проект 11-01-00545-а).

а

Библиографический список

1. Коненков Ю. К. Об изгибной волне «рэлеевского» типа // Акустический журнал. 196g. Т. б, вып. 1. С. 124-12б.[Konenkov Yu. K. A Rayleigh-Type Flexural Wave // Soviet Phys. Acoustics. 1960. Vol. 6, iss. 1. P. 122-123.]

2. Norris A. N. Flexural edge waves // J. of Sound and Vibration. 1994. Vol. 171. P. 571-573.

3. Thompson I., Abrahams I. D., Norris A. N. On the existence of flexural edge waves on thin orthotropic plates // J. Acoust. Soc. America. 2002. Vol. 112. P. 1756-1765.

4. Zakharov D. D., Becker W. Rayleigh type bending waves in anisotropic media // J. Sound and Vibration. 2003. Vol. 261. P. 805-818.

УДК 629

АНАЛИТИЧЕСКОЕ РЕШЕНИЕ ЛИНЕЙНЫХ ДИФФЕРЕНЦИАЛЬНЫХ УРАВНЕНИЙ ОШИБОК БИНС, ФУНКЦИОНИРУЮЩЕЙ В НОРМАЛЬНОЙ ГЕОГРАФИЧЕСКОЙ СИСТЕМЕ КООРДИНАТ, ДЛЯ СЛУЧАЯ ДВИЖЕНИЯ ВДОЛЬ ЭКВАТОРА С ПОСТОЯННОЙ СКОРОСТЬЮ НА ПОСТОЯННОЙ ВЫСОТЕ

М. Ю. Логинов1, М. Г. Ткаченко2, Ю. Н. Челноков3

1 Институт проблем точной механики и управления РАН, Саратов E-mail: [email protected] 2Автономный университет Мехико E-mail: [email protected] 3Саратовский государственный университет E-mail: [email protected]

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

Ключевые слова: инерциальная навигация, бесплатформенные инерциальные навигационные системы, ошибки БИНС, аналитическое решение уравнений ошибок, экватор.

Analytical Solution of Linear Differential Error Equations of Strapdown Inertial Navigation System, Functioning in the Normal Geographic Reference Frame, for the Case of an Object, Following the Geographical Equator

M. Yu. Loginov, M. G. Tkachenko, Yu. N. Chelnokov

Analytical solution of linear differential error equations of the strapdown inertial navigation system, functioning in the normal geographic reference frame, for the object, following the Earth equator with constant speed and on the constant height, is derived. The solution is represented in the form, which is convenient for the analysis. The roots of the auxiliary equation are derived in the explicit form. Obtained results can be used, for example, for analysis of the accuracy of strapdown inertial navigation system.

Key words: inertial navigation, strapdown inertial navigation systems, strapdown INS errors, analytical solution of error equations, equator.

ВВЕДЕНИЕ

При построении алгоритмов функционирования БИНС используются так называемые уравнения идеальной работы БИНС, т.е. дифференциальные и функциональные соотношения, связывающие проекции векторов кажущегося ускорения и абсолютной угловой скорости объекта, измеряемые чувствительными элементами БИНС (при условии их идеального функционирования), с навигационными параметрами (координатами местонахождения и проекциями скорости) и параметрами ориентации. Возможны различные варианты таких уравнений [1-4]. В данной работе используются уравнения идеального функционирования БИНС в нормальной географической системе координат (НГСК), в которых в качестве промежуточных кинематических параметров ориентации используются параметры Эйлера (Родрига-Гамильтона) [4-10]. Эти уравнения используются в настоящее время для построения высокоточных алгоритмов функционирования современных отечественных БИНС, построенных на волоконно-оптических или лазерных гироскопах и кварцевых акселерометрах.

В работе [11] для этих уравнений выведены полные и линеаризованные дифференциальные урав-

© Логинов М. Ю., Ткаченко М. Г., Челноков Ю. H., 2G13

69