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

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

CC BY
62
15
i Надоели баннеры? Вы всегда можете отключить рекламу.
Область наук
i Надоели баннеры? Вы всегда можете отключить рекламу.
iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.
i Надоели баннеры? Вы всегда можете отключить рекламу.

Текст научной работы на тему «Нелинейные осесимметричные продольные волны в неоднородной цилиндрической оболочке»

Полученные результаты позволяют сделать следующие выводы.

Прогиб достигает наибольшего значения в точках внешнего контура

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

на контуре спая ; в случаях I, III, IV - изгибающий момент Мд3' на внутреннем контуре Ь\3). Величина же его зависит кроме этого и от ширины составляющих колец. Так, например, в круглой плите (а® = 5, 42) = 3, 43) = 1, а{3) = 0.25) для случая V - Мтах = 2.56 • т, если же в этой плите уменьшить а® до 4 , то Мтях =3.39 -т.

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

1. Лехницкий С. Г. Анизотропные пластинки. М.: Гостехиздат, 1957.

2. Копнина В. И. Чистый изгиб составной анизотропной плиты // Механика деформируемых сред. Саратов: Изд-во Сарат. ун-та, 1985. Вып. 9. С. 45 - 53.

3. Копнина В. И., Меглинский В. В. Квазирегулярность бесконечной системы в задаче об изгибе составной анизотропной плиты // Механика деформируемых сред Саратов: Изд-во Сарат. ун-та, 1974. Вып. 1. С. 57 - 64.

УДК 539.3:534:532.5

Л. И. Могилевич

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

Рассмотрим цилиндрическую оболочку для модели Кирхгофа-Лява, усиленную в двух направлениях взаимно перпендикулярными ребрами, расположенными симметрично по обе стороны срединной поверхности. Считая материал обшивки нелинейно упругим, примем для него кубическую зависимость интенсивности напряжений от интенсивности деформаций [1]

сг, = Ее1 - те] (¡ = х,у), (1)

где Е - модуль Юнга; т - константа материала, определяемая из опытов на растяжение и сжатие; сг - интенсивность напряжений; е,- интенсивность деформаций.

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

м.

=^[л (*, + л*,)];", = Т~2"[71 (*, - = = {х* + Я*,)] - = + " (2)

2(1 +у") -"„е,

При этом связи между деформациями и перемещениями имеют стандартный вид

ди дУ IV ди дУ

£х " дх '£у ~ ду

у 'ду ' дх'

дУ_

32Ж ду

д21¥

дУ

>ЛЛЛ' т

(3)

дх* 'лу я ' ду2 ,Лзу 2Я дхду' Здесь и, У, \У - перемещения срединной поверхности вдоль осей х, у и по нормали; Ех,Ег - модули Юнга материалов, из которых изготовлены ребра, параллельные у и х; - моменты инерции сечений этих ребер относительно линий, проходящих через их центры тяжести; с1^с12 - расстояния между ребрами в направлении у их; И - толщина обшивки оболочки; /л -коэффициент Пуассона.

Уравнения движения, записанные в усилиях и моментах, имеют следующий вид:

дЫх ^дТ у д2Ц. дЫу дТ _у нд2У. дх ду g дг2 ду дх g дг2 '

д2Мг

дх2 +2

д2Н

д2му 1 диу д

---У. +---У-+-

ду2 К ду дх

у . д2\У -Л-

дх ду

ду

дх ду

(4)

дхду g 6Г2

здесь у - удельный вес обшивки; g - ускорение свободного падения. Подставляя (1) - (3) в (4), получим уравнения движения в перемещениях Переход к безразмерным переменным (в этих уравнениях)

I

I

у(1-м2)Г

(5)

где А - амплитудный параметр возмущения, / - характерная длина волны, позволяет выявить четыре малых параметра подобия

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

Рассмотрим случай, когда

Примехгяя метод возмущений, представим решение в виде

= х* - с?,т = еУ ,т] = е2у ,и* =и0 + +..., IV' = + еЩ + ...У = е(У0 +еУ1 +...) . Учтем осевую симметрию задачи и потребуем выполнения условия

(7)

У,/? _

М,/4

В нулевом приближении получим

(8)

(9)

В следующем приближении будем иметь для компоненты продольной де-

ди0

формации Ф =

эф

дт

а?

1

Э3Ф 35Ф -+ со—г = 0.

д?

(10)

Полученное уравнение (10) ранее в литературе по теории солитонов и в приложениях не упоминалось.

Уравнение (10) имеет точное решение солитонного типа

л/10 асо 5

1 П}е 2 Гр*

-лГЧ +-л —Т т

2\<о 25 V со

\ у

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

^ = Т •

\Её

ЛИТЕРАТУРА

1. Землянухин А. И., Могилевич Л. И Нелинейные волны деформаций в цилиндрических оболочках // Изв. вузов. Сер. Прикладная нелинейная динамика. 1995. Т. 3. № 1. С. 52-58.

УДК 629

А. В. Молоденков, Я. Г. Сапунков

АНАЛИТИЧЕСКОЕ РЕШЕНИЕ ИМПУЛЬСНОЙ ЗАДАЧИ ОПТИМАЛЬНОГО РАЗВОРОТА КОСМИЧЕСКОГО АППАРАТА*

Рассматривается задача оптимального разворота сферически симметричного космического аппарата (КА) с ограниченной и импульсной тягой при произвольных граничных условиях по угловому положению и угловой скорости в кватернионной постановке. Функционал, определяющий качество переходного процесса, объединяет два критерия: время и интегральную величину импульсов тяги, затраченных на процесс управления. Постановка задачи имеет вид

2А = А_ою, (1)

ш = М, (2)

(3)

Л(0) = Ло, й(0 ) = ш0> (4)

Л(Т) = Лт, ш(Т) = аГх, (5)

_ з

где Л(Г) = А,0(0 + ' кватернион, описывающий угловое положе-

к=1

з

ние КА, ю(?)= ^G>k(t)ik - вектор угловой скорости КА - являются пере-ы\

менными состояния; M(t) - [Ml(f),M2(t),M3(?)]7 - вектор управляющего момента. Функции A(i), oa(i), Af(f) подчинены требованиям задачи Пон-трягинского типа; о означает кватернионное произведение. Требуется определить оптимальное управление Мот системой (1), (2) при ограничении на модуль управления (3), доставляющее минимум функционалу

г

/ = /(а1 + а2|АфйГ, о

* Работа выполнена при поддержке Российского фонда фундаментальных исследований, проект № 99-01-00192.

i Надоели баннеры? Вы всегда можете отключить рекламу.