Механика
УДК 531.01+531.552
МНОГОПАРАМЕТРИЧЕСКОЕ СЕМЕЙСТВО ФАЗОВЫХ ПОРТРЕТОВ В ДИНАМИКЕ ТВЕРДОГО ТЕЛА, ВЗАИМОДЕЙСТВУЮЩЕГО
СО СРЕДОЙ
М. В. Шамолин1
Исследуется задача плоскопараллельного движения однородного симметричного твердого тела, взаимодействующего со средой лишь через плоский участок его внешней поверхности. При построении силового поля используется информация о свойствах струйного обтекания в условиях квазистационарности, при этом движение среды не изучается, а рассматривается такая задача динамики твердого тела, в которой его характерное время движения относительно своего центра масс соизмеримо с характерным временем движения самого центра масс.
Ключевые слова: твердое тело, сопротивляющаяся среда, система с переменной диссипацией, фазовый портрет.
A problem of plane-parallel motion of a homogeneous symmetric rigid body interacting with a medium only through a flat region of its outer surface is studied. The force field is constructed on the basis of data on the properties of jet flow under quasistationarity conditions. The motion of the medium is not studied. The problem of rigid body dynamics is considered for the case when the characteristic time of motion of the body relative to its center of mass is comparable with the characteristic time of motion of this mass center.
Key words: rigid body, resisting medium, system with variable dissipation, phase portrait.
Работа представляет собой исследование задачи движения симметричного твердого тела, взаимодействующего со средой лишь через плоский участок его поверхности. При построении силового поля используется информация о свойствах струйного обтекания в условиях квазистационарности, при этом движение среды не изучается [1—3]. В отличие от [3, 4], где зависимость момента силы воздействия среды от угловой скорости пренебрежимо мала, в данной работе в соответствии с экспериментом [5] учитываются эффекты от влияния вращательных производных момента гидродинамических сил по угловой скорости [6, 7].
С практической точки зрения представляет существенный интерес изучение устойчивости прямолинейного поступательного (невозмущенного) движения, при котором скорости точек тела перпендикулярны плоской пластине. Полное нелинейное исследование продиктовано необходимостью отыскания таких условий, при которых существуют колебания ограниченной амплитуды возле невозмущенного движения [5].
В предлагаемой работе применяется разработанная автором методика полного нелинейного исследования диссипативных динамических систем определенного вида, благодаря чему получено новое многопараметрическое семейство фазовых портретов с предельными циклами на фазовом цилиндре, состоящее из бесконечного множества топологически неэквивалентных фазовых портретов.
1. Постановка задачи. Пусть однородное тело массы т совершает плоскопараллельное движение в среде, при этом передняя часть поверхности — плоская пластина AB (рис. 1) — находится в условиях струйного обтекания. В отсутствие касательных сил воздействие среды на пластину сводится к силе S (приложенной в точке N) и ей ортогональной. Остальная часть поверхности может быть размещена внутри объема, ограниченного струйной поверхностью, срывающейся с края пластины. Похожее явление возникает после входа тела в воду [2, 5]. Сила тяжести считается пренебрежимо малой по сравнению с сопротивлением. координат Dxyz (ось z перпендикулярна плоскости рисунка, Dzx — плоскость симметрии). Тогда существует режим невозмущенного движения, перпендикулярного пластине AB. При этом перпендикуляр Dx, опущенный из центра тяжести C тела на пластину, принадлежит
D
ЧА
У
Рис. 1
Свяжем с пластиной систему
1 Шамолин Максим Владимирович — доктор физ.-мат. наук, проф., вед. науч. сотр. НИИ механики МГУ; e-mail: shamolin@imec.msu.ru.
линии действия силы S, а при возмущении данного режима скорость v точки D отклоняется от оси Dx на угол а (атаки).
Введем фазовые координаты: v = |v|, а, Q — угловая скорость (AB = А). Примем далее квадратичную зависимость величины силы S от скорости: S = s\v2. При этом обычно считают, что si = pPcx/2, где р — плотность среды, P — площадь пластины. Безразмерный коэффициент лобового сопротивления cx зависит от а, числа Струхаля и других величин, при этом в статических моделях их считают параметрами. Мы же вводим безразмерную фазовую переменную "типа Струхаля" и = QA/v, а также функцию s = si sgncos а. Ограничимся случаем зависимости cx только от а, т.е. будем считать s функцией а, а yN = DN — функцией пары а, и.
При невозмущенном движении а(Ь) = 0, u(t) = 0, поэтому yN при малых а, и примем в виде
yN = А(ка — Ни), (1)
где к и h — некоторые постоянные. Зависимостью s от а в силу симметрии тела пренебрегаем. Теперь модель содержит параметры s = si, к, h, определяемые формой пластины, где лишь первый параметр размерный [2, 5].
Положение тела на плоскости зададим координатами (xo,yo) точки D и углом ф. Полярные координаты (v, а) конца вектора скорости точки D и величина Q связаны с переменными (x0, y*,ф*, Ф) следующим образом:
ф = Q, x0 = v cos(а + ф), y0 = v sin(а + ф). (2)
Состояние системы определяется функциями (и,а, Q,xo,Уо,Ф), первые три из которых рассматриваются в качестве квазискоростей. Ясно, что тройка координат (xo,Уо,Ф) циклическая, что приводит к понижению порядка системы. Динамическая же часть образует замкнутую систему (а = DC, I — центральный момент инерции)
v* cos а — а^v sin а — Qv sin а + aQ2 = —s^v2/^,, v* sin а + а^v cos а + Qv cos а — aQ* = 0,
(3)
IQ* = yN (а,u)s(а)v2, и = QA/v.
2. Классы функций, определяющих воздействие среды. В систему (3) входят функции yN (а, и), s, определяющие воздействие среды. Если, в частности, yN = у(а), то вид y устанавливается через эксперимент [2-5]. Цель работы — учет влияния вращательных производных силы воздействия среды по угловой скорости, что требует введения в функции воздействия дополнительного аргумента. Это — нетривиальная задача моделирования. Учитывая соотношение (1), функцию yN берем в виде yN(а, и) = yN(а, Q/v) = у(а) — HQ/v (H > 0 [2, 5]), при этом
IQ* = F (аф2 — Hs^)Qv, (4)
где F(а) = у(а^(а).
Система (3) и уравнение (4) содержат функции F(а), s^), явный вид которых даже для пластин простой формы аналитически описать затруднительно, поэтому используется прием "погружения" задачи в более широкий класс задач, учитывающий лишь качественные свойства функций F(а), s^). Построить же классы {у}, {s} помогает результат Чаплыгина, получившего функции у(а), s^) для плоскопараллельного обтекания пластины бесконечной длины [8]:
у(а) = у0(а) = A sin а £ {у}, A> 0; s^) = s0^) = B cos а £ {s}, B > 0. (5)
Сочетая (5) с экспериментом [2, 4, 7, 8], формально опишем данные классы, состоящие из гладких 2п-периодических функций (у нечетная, а s четная), таких, что у(а) > 0 при а £ (0,п), у'(0) > 0, у'(п) < 0 ({у} = Y); s^) > 0 при а £ (0,п/2), s^) < 0 при а £ (п/2,п), s(0) > 0, s'(n/2) < 0 ({s} = Я). Функции у и s меняют знак при замене а на а + п. Далее, введенная в (4) функция F гладкая нечетная п-периодическая, такая, что F(а) > 0 при а £ (0,п/2), F'(0) > 0, F'(п/2) < 0 ({F} = Ф). Типичным представителем класса Ф является функция [8] F = Fo (а) = AB sin а cos а £ Ф.
3. Ключевой вопрос для приложений. В связи с отмеченной ранее неустойчивостью можно поставить вопрос: существуют ли угловые колебания оси симметрии тела конечной (ограниченной) амплитуды? Сформулируем более общо: найдется ли такая пара функций у и s, чтобы для некоторого решения уравнений движения выполнялось ограничение 0 < o:(t) < а* < п/2 начиная с некоторого момента t = ti? Ранее показано [3, 4], что, когда величины yN и s зависят лишь от угла атаки (т.е. H = 0),
13 ВМУ, математика, механика, №3
для любой допустимой пары функций у(а) и s(a) во всем диапазоне (0 < а < п/2) конечных углов атаки в системе отсутствуют нетривиальные ограниченные решения. Таким образом, положительный ответ на поставленный выше вопрос возможен в случае учета зависимости момента силы воздействия среды от угловой скорости. Как будет показано, при некоторых предположениях можно ожидать положительного результата.
4. Приведенные уравнения движения. Вводя новую независимую переменную по формуле dq = v dt, первые два уравнения системы (3) и уравнение (4) перепишем в виде
v' = Ф(а, u)v, (6)
а = —ш + аш2 sin а + aF(а) cos а/I + s(a) sin а/m — аИшв(а) cos а/I,
ш' = F(a)/I + аш3 cos а — awF(а) sin а/I + шв(а) cos а/m — Ишз(а)/1 + аИш2з(а) sin а/I, (7) Ф(а, ш) = —аш2 cos а + aF (а) sin а/I — s^) cos а/m — аИшз(а) sin а/I, (...)' = d/dq, ш = d^/dq.
Система (7) образует независимую подсистему второго порядка. Исследуем устойчивость ее тривиального решения, введя обозначения A = y'N(0), B = s(0), n0 = F(0)/I и безразмерные параметры
Hi = 2B/mn0 > 0, ц2 = an0 > 0, /л3 = BH/In0 > 0. (8)
Предложение 1. Пусть выполнено неравенство
H — H2\ < 2. (9)
Тогда при H3 > Hi + H2 (соответственно при H3 < Hi + H2) тривиальное решение системы (7) асимптотически устойчиво (соответственно является отталкивающим).
Для выяснения возможности рождения предельного цикла около начала координат исследуем устойчивость тривиального решения системы (7) при критическом сочетании параметров Нз = Hi + Н2. Для этого вычислим вспомогательный индекс [3, 9, 10] (s2 = s''(0), ¡3 = F'''(0)):
Ind = -5/X2 - Hi + Ф- + 2 — - ^ + 8/X2 --1 + /X'/4 Л2 (5в - 3s2)U + 12/x2f/2
In0 mn0 In0 (l + H2(Hi + H2))2 In0
H2 + Hi/2
U =
l + H2(Hi + H2)'
Предложение 2. Если 1пё < 0 (соответственно 1пё > 0) и выполнено неравенство (9), то начало координат фазовой плоскости системы (7) при критическом сочетании параметров Нз = Н1 + Н2 является слабоустойчивым (соответственно неустойчивым) фокусом,.
Теорема 1. Пусть для системы (7) выполнено неравенство (9). Тогда справедливы следующие утверждения.
1) Если 1пё < 0, то для любого фиксированного Н2 найдутся такие 5\, > 0, что при Нз € (Н2, Н2 + ¿1) начало координат является сильноустойчивым фокусом,; при Нз € (Н2 — $2,Н2) начало координат является силънонеустойчивым фокусом,, окруженным устойчивым предельным циклом, размер которого растет с уменьшением, Нз от Н2 до Н2 ~ $2, как \/\Н2 ~ Нз|-
2) Если 1пё > 0, то для любого фиксированного Н2 найдутся такие 51, 52 > 0, что при Нз € (Н2 — 52,Н2) начало координат является сильнонеустойчивым фокусом; при Нз € (Н2,Н2 + ¿1) начало координат является сильноустойчивым фокусом,, окруженным неустойчивым предельным циклом, размер которого 'растет с увеличением Нз от Н2 до Н2 + $2, как л/\Н2 — /хз|•
Проверить выполнение условия Нз > Н1 + Н2 несложно. Проверка же условия 1пё < 0 затруднительна, поскольку для каждого конкретного тела не только явный вид, но и старшие производные даже в отдельных точках пары функций (уN,в) неизвестны.
5. Точки покоя систем. У системы третьего порядка (6), (7) существуют положения равновесия, которые заполняют одномерные многообразия. Поэтому задача о точках покоя разбивается на две — о точках покоя системы третьего порядка в фазовом пространстве W = Е+ {г} х Я2{а, и} и о точках покоя системы (7) на плоскости Я2{а,и} или на цилиндре Б1 {а шоё2п} х Я1^}.
5.1. Точки покоя системы третьего порядка. Искомые точки покоя зададим системами
а = п/2+п1, I € Z; и = 0, г = г1; (10)
а = —п/2 + 2п1, l £ Z; и = —1/а, а = 0, v = v2; (11)
а = п/2 + 2п1, l £ Z; и = 1/а, а = 0, v = va, (12)
которые соответствуют в фазовом пространстве W прямым, сплошь заполненным точками покоя системы (ср. с [3]).
5.2. Тривиальные и нетривиальные точки покоя системы второго порядка. У системы (7) существуют точки покоя на плоскости В2{а,и}, которые задаются следующими совокупностями:
а = 2пк, к £ Z; и = 0; (13)
а = (2к + 1)п, к £ Z; и = 0; (14)
а = п/2 + п1, l £ Z; и = 0; (15)
а = —п/2 + 2п1, l £ Z; и = —1/а, а = 0; (16)
а = п/2 + 2nl, l £ Z; и = 1/а, а = 0. (17)
Системы (15)-(17), и только они, задают точки покоя, в которые ортогонально проектируются многообразия положений равновесия (10)-(12) соответственно. Системы (13), (14), и только они, задают точки, в которые ортогонально проектируются частные решения системы (6), (7) вида
v(q) = v° ехр^—— q*J , v° = v(0), a(q)=7rk, к G Z, u(q)=0,
при к = 0 соответствующие невозмущенному движению.
Точки покоя (13)-(17) назовем тривиальными положениями равновесия (ТПР). Наряду с ТПР могут существовать точки покоя, не лежащие на прямых {(а, и) £ R2 : и sin а cos а = 0} и соответствующие нетривиальным частным решениям системы (2), (3); их будем называть нетривиальными (НПР).
5.3. Нетривиальные частные решения. Рассмотрим систему (7) при условиях (5). Аналогично может быть рассмотрена система общего вида (7). Если НПР существуют, то они удовлетворяют следующим равенствам:
и = а sin а, [аи2 — аП2 + cos2 а — B cos а/m + а(1 — аа) = 0; (18)
а = а± = n0
¡12 +Из ± ^(№+/хз)2-4 /(2(1 + М2Мз))- (19)
Следствие. При выполнении условий (9), (5) система (7) не имеет НПР. Замечание 1. Если у системы (7) и существуют НПР, то они удовлетворяют уравнению
и2
i + ^я s(a)
I cos а
и
а F(а) 1Г s(а)
+ - Н - sin а
1 F (а) , ч
+ -—^sina = 0. (20)
I cos а
I cos а I cos а
Дискриминант уравнения (20) отрицателен, если почти для любого а выполнено неравенство
^F(а) + Hs^) sin а)2 < 4IF^) sin а cos а. (21)
Если же (21) не имеет места, то и у системы (7) могут существовать положения равновесия (ао, ио), удовлетворяющие системе (18), (19). Они индуцируют наличие следующих нетривиальных решений:
v(q) = v0eKq, v0 = v(0), k< 0, а^) = а0, а0 = а(0) = пк, к £ Z, и(д) = и0, и0 = и(0).
Замечание 2. У системы (2), (3) существуют частные решения вида
v(t) = vo(1 — voKt)-i, v0 = v(0), аф) = ао, ао = а(0), Q(t) = d^(t)/dt = Qo(1 — voKt)-i, Qo = Q(0);
^(t) = фо — Q ln(1 — voKt)/vo К, фо = ф(0); x(t) = xo + vo sin(аo + фо )/Qo — vo sin(аo + фо + Qo ln(1 — vont)/vo/k)/Qo , xo = x(0); y (t) = yo — vo cos(аo + Ф0 )/Qo + vo cos^o + фо + Qo ln(1 — voKt)/vo/k)/Qo , Уо = y(0).
14 ВМУ, математика, механика, №3
Замечание 3. Траектория точки О, соответствующая последним частным решениям, — окружность с центром в (хо + го 8ш(ао + ро)/0о,уо — го со8(ао + )/0о) и радиусом го/0о.
5.4. Топологическая классификация ТПР. 1) Точки (13) являются: седлами, если 2Н1(Нз — Н2) > 4 + н2; отталкивающими, если 2н1 (Нз — Н2) < 4 + Н2, Нз < Н1 + Н2; притягивающими, если 2н1(нз — Н2) < 4 + Н1, Нз > Н1 + Н2. 2) Точки (14) являются: седлами, если 2Н1(Н2 — Нз) > 4 + н2; отталкивающими, если 2Н1(Н2 — Нз) < 4 + н2, Нз > —Н1 + Н2; притягивающими, если 2н1(Н2 — Нз) < 4 + н2, Нз < —Н1 + Н2. 3) Точки (15) — седла. 4) Точки (16), (17) являются притягивающими: устойчивыми фокусами, если Н1Н2 < 4; устойчивыми узлами, если Н1Н2 > 4; устойчивыми вырожденными узлами, если Н1Н2 = 4.
6. Расслоения фазового пространства, его симметрии. Фазовые траектории системы (6), (7) лежат на поверхностях, представляющих собой двумерные цилиндры. В частности, если существует во всем фазовом пространстве первый интеграл системы (7), то он является функцией (а,и), а поэтому задает семейство цилиндров в пространстве W = Я+ {г} х Я2{а, и}.
Замечание 4. Поскольку йф/М = гй<р/йд, то в пространстве Я+ {г} х Я2{а, 0} системы (3) данному первому интегралу будет соответствовать первый интеграл, который зависит от всех трех фазовых переменных. Поэтому удобнее строить фазовый портрет системы (6), (7) в пространстве W с помощью портрета системы (7) на плоскости Я2{а, и}, который будет являться ортогональной проекцией на плоскость {(а, и, г) € Яз : г = еопв^. Тогда возможно поднятие траекторий с плоскости в пространство W.
Замечание 5. Движение возможно лишь в области В' = {(а,и,г) € Яз : г > 0} С Яз{а,и,г}. После замены р = 1пг полученное векторное поле в фазовом пространстве Яз{а,и,р} не зависит от р и ортогонально проектируется на семейство плоскостей {(а,и,р) € Яз : р = еопв^.
Замечание 6. Чтобы поднять траектории в трехмерное пространство с плоскости, необходимо выяснить знак проекции поля в Яз{а,и,р} на ось р либо на ось г в области В' и рассмотреть цилиндры М = М°иЖ, М° = {(а,ш,у) е в' : 4>(а,ш) сов'1 а = 0}, М = Мг и М2 и ..М3 = {(а,ш,у) е В' : а = п/2 + пз}, разрезающие фазовое пространство на области, в каждой из которых проекция поля на ось р имеет свой знак. На М проекция на ось г нулевая.
Замечание 7. Для любых ^ € Ф, в € £ векторное поле системы (7) обладает свойством центральной симметрии относительно точек (пк, 0), к € Z: в координатах (а, и) поле меняет направление при замене (а + пк,и) (—а + пк, —и), к € поэтому поле системы (6), (7) обладает симметрией относительно лучей {(а, и, г) € В' : а = пк, и = 0}.
7. Классификация фазовых портретов системы на двумерном цилиндре для некоторой области параметров. 7.1. Введение в классификацию фазовых портретов. Поведение поля системы возле точек (0, 0) и (п, 0) определяет основные особенности глобального портрета. Изучим те системы, которые допускают выполнение неравенства (21) лишь на ТПР.
Бесконечномерное пространство параметров системы обозначим через ,11. В работе [11] показано, что существуют и единственны траектории, уходящие на бесконечность. Их а- и и-предельными множествами являются бесконечно удаленные точки (пк, к € Z .В [3, 9] установлено, что отсутствуют замкнутые
траектории или замкнутые кривые из траекторий, охватывающих цилиндр, при этом у системы (7) также отсутствуют замкнутые фазовые характеристики, охватывающие цилиндр.
В [3] доказано существование замкнутых кривых из траекторий, стягиваемых по цилиндру в точку, например существование неизолированных периодических траекторий или предельных циклов. Эта задача сводится к отысканию замкнутых кривых из траекторий лишь вокруг изолированных точек покоя индекса 1.
Основной проблемой классификации портретов является изучение поведения устойчивых и неустойчивых сепаратрис имеющихся гиперболических седел.
7.2. Некоторые свойства ключевых фазовых траекторий системы в области ,1. Введем семейство полос на цилиндре П(а1 ,а2) = {(а,и) € Я2 : < а < а2}, при этом П(-п/2п12) = П, П(п12>зп12) = П'.
Следующее утверждение относится ко всему бесконечномерному пространству параметров и базируется на наличии топографической системы Пуанкаре (ТСП) в полосе П [12, 13].
Предложение 3. Устойчивая сепаратриса в полосе П' для точки (п/2,0) имеет в качестве а-предельного множества точку (п, 0) [3] (см. рис. 2, а—в).
Следствие. У системы (7) в полосе П' не существует никаких замкнутых характеристик.
8. Многопараметрическое семейство фазовых портретов системы на двумерном цилиндре. Динамическая система (7) имеет на двумерном фазовом цилиндре следующие положения равновесия:
(0, 0) и (п, 0), (п/2, 0) и (3п/2, 0), (п/2,1/а) и (—п/2,— 1/а).
(22)
Теорема 2. В пространстве параметров .]\ существует область , (которая далее и рассматривается) положительной меры, соответствующая следующему поведению траекторий системы: 1) других положений равновесия, кроме (22), система не имеет; 2) в полосе П' система не имеет замкнутых фазовых характеристик; 3) в полосе П около положения равновесия (0, 0) при изменении параметров системы (8) может происходить бифуркация рождения единственного устойчивого цикла из слабого фокуса.
Замечание 8. Области параметров ,2 соответствует следующее поведение траекторий возле положений равновесия (22): 1) положение равновесия (п, 0) является отталкивающим; 2) положение равновесия (0,0) является отталкивающим, если ¡1 > ¡э — ¡2, и притягивающим, если ¡1 ^ ¡э — ¡12; если ц 1 = Цэ — ¡2, то оно является слабым притягивающим фокусом. Замкнутые характеристики системы для области параметров ,2 могут существовать лишь в полосе П [3, 9].
Рассмотрим ключевой вопрос о глобальном поведении следующих сепаратрис:
1) выходящей из точки (п/2, 0) в полосу П'; 2) входящей в точку (—п/2, 0) из полосы П; 3) выходящей из точки (п/2, 0) в полосу П.
Теорема 3. Глобальное поведение любых двух сепаратрис 1-3 независимо [3,10]. рис ^
Выберем в качестве пары ключевых сепаратрисы 1 и 2.
Определение 1. Индексом к\ сепаратрисы 1 назовем рациональное число, выбираемое из множества {г € Q : г = 1/4 + т,г = 1/2 + т,т € N0}. Скажем, что к1 = г, если сепаратриса 1 имеет в качестве ^-предельного множества точку (2пг, 1/а) в случае г = 1/4 + т и точку (2пг — 0, в случае г = 1/2 + т.
Определение 2. Индексом к2 сепаратрисы 2 назовем натуральное число ], выбираемое из множества {] € N : ] = 1, 2, 3, 4, 5}. Скажем, что к2 = ], если сепаратриса 2 имеет в качестве а-предельного множества точку (0, 0) или устойчивый предельный цикл, ] = 1; точку (п/2, 0), ] = 2; точку (п, 0), ] = 3; точку (—0, —то), ] = 4; точку (— п, 0), ] = 5.
Теорема 4. Для любого к из (возможно, усеченной) области определения допустимо соответствующее поведение сепаратрис 1 и 2 (определения корректны, см. также [4]).
Теорема 4 позволяет сделать вывод: любое достаточно малое возмущение системы, описывающей физический маятник, которое дает искомую систему, бесконечно много раз перестраивает глобальный тип невозмущенного фазового портрета. Портреты, индекс к которых равен (1/4, 2)*, (1/4, 4)*, (1/4, 5), показаны соответственно на рис. 2, а, б, в. Звездочкой помечены портреты, имеющие предельные циклы. Для системы (7) при выполнении условий (5) имеем трехпараметрическое семейство портретов, зависящее от параметров (8). Двухпараметрическое семейство, построенное в [4], не содержит предельных циклов. Но эти два семейства объединяет тот факт, что каждым значениям безразмерных параметров задачи соответствует пара независимых индексов (в данном случае к1, к2), "кодирующих" тип портрета.
15 ВМУ, математика, механика, №3
Работа выполнена при финансовой поддержке РФФИ, грант № 08-01-00231-а.
СПИСОК ЛИТЕРАТУРЫ
1. Ламб Г. Гидродинамика. М.: Физматгиз, 1947.
2. Локшин Б.Я., Привалов В.А., Самсонов В.А. Введение в задачу о движении тела в сопротивляющейся среде. М.: МГУ, 1986.
3. Шамолин М.В. Методы анализа динамических систем с переменной диссипацией в динамике твердого тела. М.: Экзамен, 2006.
4. Шамолин М.В. Новое двухпараметрическое семейство фазовых портретов в задаче о движении тела в среде // Докл. РАН. 1994. 337, № 5. 611-614.
5. Ерошин В.А., Самсонов В.А., Шамолин М.В. Модельная задача о торможении тела в сопротивляющейся среде при струйном обтекании // Изв. РАН. Механ. жидкости и газа. 1995. № 3. 23-27.
6. Бюшгенс Г.С., Студнев Р.В. Динамика продольного и бокового движения. М.: Машиностроение, 1969.
7. Бюшгенс Г.С., Студнев Р.В. Динамика самолета. Пространственное движение. М.: Машиностроение, 1988.
8. Чаплыгин С.А. О движении тяжелых тел в несжимаемой жидкости // Полн. собр. соч.: В 3 т. Т. 1. Л.: Изд-во АН СССР, 1933. 133-135.
9. Шамолин М.В. Замкнутые траектории различного топологического типа в задаче о движении тела в среде с сопротивлением // Вестн. Моск. ун-та. Матем. Механ. 1992. № 2. 52-56.
10. Шамолин М.В. Некоторые задачи дифференциальной и топологической диагностики. 2-е изд., перераб. и доп. М.: Экзамен, 2007.
11. Шамолин М.В. Существование и единственность траекторий, имеющих в качестве предельных множеств бесконечно удаленные точки, для динамических систем на плоскости // Вестн. Моск. ун-та. Матем. Механ. 1993. № 1. 68-71.
12. Шамолин М.В. Применение методов топографических систем Пуанкаре и систем сравнения в некоторых конкретных системах дифференциальных уравнений // Вестн. Моск. ун-та. Матем. Механ. 1993. № 2. 66-70.
13. Шамолин М.В. Пространственные топографические системы Пуанкаре и системы сравнения // Успехи матем. наук. 1997. 52, вып. 3. 177-178.
Поступила в редакцию 12.11.2008
УДК 539.3
К ЗАДАЧЕ ТЕОРИИ УПРУГОСТИ В ПЕРЕМЕЩЕНИЯХ
ДЛЯ ЦИЛИНДРИЧЕСКОГО СЛОЯ С СИЛЬНО РАЗЛИЧАЮЩИМИСЯ ХАРАКТЕРНЫМИ РАЗМЕРАМИ
Т. И. Гаряева1, Д. В. Георгиевский2
Проведен анализ главных членов общих асимптотических разложений решений первой краевой задачи трехмерной теории упругости в перемещениях (квазистатика, сжимаемость) для цилиндрического слоя. Естественным малым асимптотическим параметром является отношение толщины слоя к длине образующей. При этом радиус основания может быть любого "промежуточного", включая концы, порядка. Такой геометрией обладает, например, цилиндрическое тело, имеющее характерные макро-, микро- и наноразмеры по различным направлениям.
Ключевые слова: теория упругости, задача в перемещениях, тонкое тело, цилиндрический слой, асимптотическое решение, система главного приближения.
An analysis of the principal terms of the general asymptotic expansions for the solutions to the 3D elasticity boundary value problem in terms of displacements (quasistatic case, compressibility) for a cylindrical layer is performed. A ratio of the layer thickness to the height of the
1 Гаряева Татьяна Ивановна — асп. каф. механики композитов мех.-мат. ф-та МГУ; e-mail: georgiev@mech.math.msu.su.
2 Георгиевский Дмитрий Владимирович — доктор физ.-мат. наук, проф. каф. механики композитов мех.-мат. ф-та МГУ; e-mail: georgiev@mech. math.msu.su.