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

Вестн. Сам. гос. техн. ун-та. Сер. Физ.-мат. науки. 2014. № 1(34). С.9—18

Дифференциальные уравнения

УДК 517.988.67

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

Т. Е. Бадокина

Национальный исследовательский

Мордовский государственный университет имени Н. П. Огарёва,

Россия, 430005, Саранск, ул. Большевистская, 68.

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

Ключевые слова: сверхзвуковой поток газа, прогиб пластины, аэроупругость, бифуркация, уравнение разветвления.

1. Постановка задачи. Рассматривается класс двухточечных краевых задач [1,2] для обыкновенных дифференциальных уравнений (ОДУ) четвертого порядка, описывающих прогиб удлиненной упруго опертой пластины в сверхзвуковом потоке газа, сжимаемой / растягиваемой внешними усилиями.

В безразмерных переменных задача о прогибе пластины описывается уравнением [1,2]

/ w// \ И Г1

Х2(——-Tw"+f3ow = kK(w',M,K)+9w" [(1+w/2)1/2-1 ]dx. (1)

V(1+ w/2)3/2/ J о

ISSN: 2310-7081 (online), 1991-8615 (print); doi: http://dx.doi.org/10.14498/vsgtu1258 © 2014 Самарский государственный технический университет.

Образец цитирования: Т. Е. Бадокина, “Модели многопараметрических бифуркационных задач для обыкновенных дифференциальных уравнений четвертого порядка” // Вестн. Сам. гос. техн. ун-та. Сер. Физ.-мат. науки, 2014. № 1 (34). С. 9-18. doi: 10.14498/vsgtu1258.

Сведения об авторе: Татьяна Евгеньевна Бадокина, ассистент, каф. прикладной математики, дифференциальных уравнений и теоретической механики.

E-mail address: [email protected]

9

Т. Е. Бадокина

Общий подход к решению такого типа задач рассматривается при кинематических граничных условиях

w(0) = w' (0) = 0, w''(1) = w(3)(1)=0. (2)

Здесь w = w(x) —прогиб пластины; 0 ^ x1 в d, —ж < y1 < ж, x = xi/d, 0 ^ x ^ 1 — прямоугольные координаты;

K (w', M, к)

1

1 + -—1 Mw' 2 .

2к/(к— 1)

— при одностороннем обтекании;

K (w', M, к)

1 — Mw' 2 .

2к/(к—1)

1 + -—1 Mw' 2 .

2к/(к— 1)

— при двустороннем обтекании;

х2 = h2 Т = qd

Х 12(1 — p2)d2, Eh

в

1 , = Pod

1 — p2, k = Eh

где d — ширина пластины, h — её толщина, E — модуль Юнга, p — коэффициент Пуассона, M = v/c^ — число Маха (v — скорость потока газа, —скорость звука в невозмущенном газе), к — показатель политропы, po—давление, во — коэффициент жёсткости основания. Интегральное слагаемое учитывает дополнительное усилие в срединной плоскости при прогибе.

Для вычисления малых изгибных форм в окрестностях критических значений бифуркационных параметров (число Маха M = Mo + £i, сжимающее (растягивающее) усилие T = То + £2) применяются методы теории бифуркаций [3]. Зависимость дифференциального уравнения от бифуркационных параметров, выраженная через корни соответствующих характеристических уравнений линеаризации, которые можно считать заданными точно, позволяет определить критические бифуркационные поверхности и кривые, в окрестностях точек которых строится асимптотика разветвляющихся решений в виде сходящихся по малым параметрам рядов, и тем самым решать бифуркационную задачу о прогибе пластины в точной постановке. Наибольшие трудности возникают при исследовании линеаризованной спектральной задачи, фредгольмовость которой, подобно [4], доказывается построением соответствующей функции Грина, выполненным по стандартной схеме, изложенной в монографии [5].

В известных нам работах к бифуркационным задачам аэроупругости, как правило, применялся метод Галёркина. Методами теории бифуркаций задача о прогибе прямоугольной пластины исследована в работе [6]. Прогиб удлиненной пластины при учете только одного бифуркационного параметра — числа Маха — в работах П. А. Вельмисова и Б. В. Логинова [7], П. А. Вельмисова и С. В. Киреева [8], в которых применялся метод групповых преобразований Ц. На, позволяющий сводить двухточечные граничные задачи для ОДУ четвёртого порядка к задаче Коши. В работе [9] исследована наиболее простая модельная бифуркационная задача для удлиненной пластины с двумя

10

Модели многопараметрических бифуркационных задач ...

бифуркационными параметрами (сжатие/растяжение и число Маха) при использовании методов теории ветвления решений нелинейных уравнений [3]. Эта же техника применялась в работе [10].

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

L(w) = x2w(4) — Tw" + aw' + e0w = 0, a = 1(2)ккМ (3)

с граничными условиями (2), которой отвечает характеристическое уравнение

Х2Л4 — ТА2 + аА + во = 0. (4)

Здесь множитель 1(2) в параметре а соответствует одно(дву-)стороннему обтеканию пластины сверхзвуковым потоком газа.

Соответствующая сопряжённая линеаризованная задача строится стандартными методами [5,11] —интегрированием по частям на промежутке [0,1] квадратичной формы L(w) ■ ш — и имеет вид

х2ш(4) — Тш'' — аш' + в0ш = 0 (5)

с граничными условиями

ш(0) = ш'(0) = 0, х2ш''(1) — Тш(1) = 0, х2ш(3)( 1) — Тш'(1) — аш(1) = 0.

При исследовании алгебраического уравнения (5) с тремя параметрами Т/х2, а/х2 и в0/х2 использован метод Штурма разделения корней алгебраического уравнения. Установлено, что уравнение (4) может иметь корни одного из следующих видов:

A) —y ± гё\, y ± i$2 (y, 81, 82 > 0, 82 > 81);

B) —ai, —«2, y ± i8 (ai, «2, y, 8 > 0, ai > «2);

C) —a ± в1, a ± в2 (a, въ в2 > 0, в2 < в1 < a);

A-B) —a, —a, a ± i8 (a, 8 > 0);

B-C) —a, —2y + a, y, Y (a, Y > 0, y < a < 2y).

Вид корней характеристического уравнения выбран с учётом равенства нулю суммы корней таким образом, что в невырожденных случая все искомые коэффициенты и функции являются функциями трёх аргументов, в вырожденных случаях — двух аргументов.

Асимптотика разветвляющихся решений по малым параметрам £1, £2 в точке бифуркации (Т0,М0) строится для случаев существования бифуркационных кривых.

Линеаризованная в точке бифуркации задача (1), (2) определяет фредгольмов оператор

B : C4+а[0,1] ^ Ca[0,1]

с одномерным подпространством нулей N(B) = span{^(x)} и дефектным подпространством N*(B) = span{^(x)}. Применение к уравнению (1) в разложении аналитической нелинейности

11

Т. Е. Бадокина

Bw = x2w(4) — Tow'' + aw' + e0w = R(w, e) = '3

= X2(3w'2w(4) + 3w''3 + 9w'w"w"'} + eiw" — 1(2)knw'e2 + ^w'' / w'2dx—

\2 ) 2 Jо

kK(K + 3M2w'2 + kK(K + 1) Moe2w'2 + k'K£ + 1) M,3,/» + ...

4

kK(K +1) M03w'3

6

(верхняя (нижняя) строка последнего слагаемого отвечает одностороннему (двустороннему) обтеканию пластины потоком газа) регуляризатора Шмидта [3] представляет его в виде системы

Bw = R(w, e) + £z, £ — (w, 7) = 0.

В аналитическом случае решение первого уравнения представляется в виде ряда

w = wioo £ + woioei + wooie2 + ^ wai;a£ai ea.

ai + |a|>1

Подстановка его во второе уравнение системы даёт уравнение разветвления (УР) Шмидта

L(£,e) = £ — (w(e),Y) = 0,

эквивалентное уравнению (1) в том смысле, что оно имеет то же количество малых решений, представимых в виде сходящихся рядов по тем же дробным степеням малых параметров.

При одностороннем обтекании главная часть УР, определяемая методом диаграммы Ньютона [3], имеет вид

L(£, e) = L2oo£2 + Liio£ei + Lioi£e2 + • • • = 0,

где

L2oo = kK(K4+ ^ Mo2(^'2,^), Liio = — (<£'', ф), Lioi = k«(^'y0);

Г = B-i, B = B + ( ■ ,y)z; у и z — биортогональные системы к <p £ N(B) и ф £ N*(B) соответственно. Здесь и далее

(a,b)

j a(x)b(x)dx. о

При L2oo = 0 решение задачи представляется рядом, сходящимся в малой окрестности ei =0, e2 = 0:

w(x) = — LiioeL+ Lioig2 ^(x) + o(|ei1, |e2|).

L2oo

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

12

Модели многопараметрических бифуркационных задач ...

Главная часть УР для двустороннего обтекания пластины сверхзвуковым потоком газа имеет вид

L(£, £) = £эоо{3 + Lno^ei + L101 2 + ... = 0,

где

L

3оо

кк(к + /г3 / ,j3

6

M„V3, й) - X2 (3 Л<2) + .V'3 + WV", «) -

---\<р" w'2dx,^

о

L200 = 0 Liio = —(v",'Ф), Lioi = vо(v','Ф), °o = 2ккМо. При L300 = 0 асимптотика сходящегося решения выражается формулой

Lii0£i + Li0i£2l1/2

w(x) = ±

L

300

<^(x) + O(|ei|, Ы),

где знаки £i и £2 определяются условием неотрицательности подкоренного выражения. Следовательно, при L300 = 0 имеет место бифуркация типа «вилки».

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

А. Характеристическое уравнение имеет две пары комплексно-сопряжённых корней (—7 ± i5i, у ± i62), и решение линеаризованной задачи даётся соотношением

w(x) = e Yx(ci cos(5ix) + c2 sin(^ix)) + eYX (c3 cos(62x) + c4 sin(62x)).

Множество точек бифуркации определяется равенством нулю определителя А матрицы граничных условий:

А = М2 [(72 + 62 )e-2Y + (Y2 + ^2)e27 ] +

+ sin 6i j ((y2 — ^2)^4 + (3y4 — + 672£M + 3y4^2 + Y2^2 — 4y6) sin 62+

+ 2y^2(^4 + (Y2 + 3^i)^2 — Y2^2 — 4Y2) cos 62} +

+ 26i cos 6^y((y2 — 3^2)^2 — — Y262 + 4Y2) sin 62+

+ 62 ((3y2 — б^ + 3y262 + 7y4) cos 62}.

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

^ = {(Y, 6i, 62) I 62 > 6i, y > 0, 6i > 0, 62 > 0}

13

Т. Е. Бадокина

определитель А > 0, то есть прогиб отсутствует.

B. Корнями характеристического уравнения являются два отрицательных и пара комплексно-сопряжённых чисел. С учётом равенства нулю суммы корней вводится замена a1 = 7 + a, a2 = 7 — a, a > 0. Тогда решению

w(x) = c1e-(Y+a>x + C2e-(Y-a)x + eYX (c3 cos(Sx) + C4 sin(Sx)) отвечает определитель матрицы граничных условий

А = 2aS [(a2 — y2)2e 27 + (y2 + S2)2e2Y| +

+ (a + y)2e-a sinS(4y4 — Y2(a2 + 3S2) + 4YaS2 + a2S2 — S4)+ + 2S cos S(4y3 — aY2 — Y(a2 — S2) + aS2)

— (a — y)2ea sinS(4y4 — Y2(a2 + 3S2) — 4YaS2 + a2S2 — S4) + + 2S cos S(4y3 — aY2 — Y(a2 — S2) — aS2)

> 0.

Лемма. На множестве Q = {(a, y, S) | a G (0, y),Y > 0, S > 0} прогиб отсутствует.

Доказательство. Рассмотрим определитель

А = 2aS [(a2 — Y2)2e-2Y + (y2 + S2)2e2Y ] +

+ (a + Y)2e-a [A sin S + B cos S] — (a — Y)2ea [C sin S + D cos S],

причём

A2 + B2 = C2 + D2 =

= (y2 + S2)2(S4 + 2S2(4y2 + a2) + (4y2 — a2)2) <

< (y2 + S2)2(S2 + 4y2 + a2)2,

тогда

А > 2aS [(a2 — Y2)2e-2Y + (y2 + S2)2e2Y] —

— (a + Y)2e-a(Y2 + S2)(S2 + 4y2 + a2) —

— (a — y )2ea(Y2 + S2)2(S2 + 4y 2 + a2).

Слагаемое 2aS(a2 — y2)2e-2Y > 0 является малым и его можно не учитывать. Так как y > a, для А выполняется неравенство

А > (y2 + S2) [2aS(Y2 + S2)e2a —

— (a — Y)2ea (S2 + 4y2 + a2) — (a + Y)2e-a(S2 + 4y2 + a2)].

Разложение экспонент в ряды приводит определитель к виду

14

Модели многопараметрических бифуркационных задач ...

ОО ОО

А = Е = Е П [2”+1а«(12 + i2)-

n=0 n=0

- (y - a)2(i2 + 4y2 + a2) - (-1)n(a + y)2(i2 + 4y2 + a2)]. Для нечётного n очевидно, что

an > ai(y2 + i2) + 4aY(i2 + 4y2 + a2) > 0.

Если n — чётное, то

an = 2ai(Y2 + i2)2n - (y2 + a2)(i2 + 4y2 + a2).

Начиная с некоторого k все члены ряда будут положительны. В силу неограниченного возрастания функции 2n с ростом числа n будем иметь, что afc+i > afc. Таким образом, все первые k - 1 отрицательных члена компенсируются последующими, начиная с k-того. Соответственно А > 0, и, следовательно, прогиб отсутствует.

C. Характеристическое уравнение имеет два отрицательных и два положительных корня (-a ± въ a ± в2), его коэффициенты имеют вид

T

X2

в2 + в! + 2a > 0, —2

X2

2a(e2 - в|), ^2

X2

(a2 - в2)(a2 - вг).

Таким образом, в этом случае возможно только растяжение пластины. Из теоремы Виета следует, что a > в1, в1 > в2. Определитель матрицы граничных условий для

w(x) = e-ax(ci е-в1Х + С2бв1Х) + eax(cae-e2x + С4вв2 х)

следующий:

А = 4в^2((a2 - в2)2e-2a + (a2 - в|)2е2“) -- (4a2 + (в1 - в2)2)х

х ((a + вl)2(a - в2)2e-(вl+в2) + (a - в^2^ + в2)2eвl+в2) +

+ (4a2 - (в1 + в2)2) х

х ((a + в1)2(a + в2)2e-(вl-в2) + (a - в^2^ - в2)2eвl-в2).

При a > в1 > в2 он положительный (А > 0).

В силу ограниченного объёма статьи доказательство отсутствия прогиба здесь не приводится.

A-B. Для вырожденного случая, когда уравнение (4) имеет отрицательный корень кратности два и пару комплексно-сопряжённых корней (-a, -a, a ± ii), линеаризованное уравнение (3) имеет решение

w(x) = e-ax (c1 + c2x) + eax (c3 cos(ix) + c4 sin(ix)).

15

Т. Е. Бадокина

Определитель матрицы граничных условий является функцией двух переменных:

А = a4Se-2a + 5(а2 + 52)2e2a+

+ а(8а4 — 4а5 + 3а352 — 2а252 + а54 + 254) sin5+

+ 2а25 (7а2 — 4а3 — а52 — 352) cos 5.

Бифуркационная кривая р = р(а, 5) состоит из точек (а, 5), в которых А = 0. Численный анализ показывает, что в рассматриваемом случае таких точек нет.

Замечание. Контроль правильности вычислений определителей матриц граничных условий осуществляется предельным переходом от случая А при 5i ^ 0 либо от случая В при а ^ у с последующим исследованием этих определителей в предельных случаях.

B-C. Характеристическое уравнение имеет два отрицательных и положительный корень кратности два (—а, — 2у + а, у, у), и соответствующее решение (3) записывается так:

w(x) = cie-ax + С2 e-(2Y-a)x + C3eYX + C4 xeYX.

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

А = у[а2(3у3 + 2у 2(а + 4) — ау(а — 2) — 2а2)е-а + 2у 3(а — у )eY ]e7+

+ (2у — а)2 [2а2(а — y)e-a — у(3у3 + 2у2(2 + а) —

— ау (а — 6) + 2а2^7 ]e-27+a.

На множестве Q = {(а,у) | а £ (у, 2у)} прогиб отсутствует, так как А > 0.

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

Работа выполнена при поддержке ФЦП «Научные и научно-педагогические кадры инновационной России» (соглашение 14.В37.21.0373 от 14.11.2012).

This work was supported by Federal Target Program “Research and scientific-pedagogical personnel of innovative Russia in 2009-2013”, the agreement 14.В37.21.0373 on 14/11/2012.

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

1. В. В. Болотин, Неконсервативные задачи теории упругой устойчивости, М.: Физматлит, 1961. 339 с.; V. V. Bolotin, Non-conservative problems of the theory of elastic stability, Oxford, Pergamon Press, 1963, 320 pp.

2. А. С. Вольмир, Устойчивость деформируемых систем, М.: Наука, 1967. 984 с. [A. S. Vol’mir, Ustoychivost’ deformiruyemykh sistem [Stability of deformable systems], Moscow, Nauka, 1967, 984 pp. (In Russian)]

3. М. М. Вайнберг, В. А. Треногин, Теория ветвления решений нелинейных уравнений, М.: Наука, 1969. 524 с.; M. M. Vainberg, V. A. Trenogin, The Theory of Branching of Solutions of Nonlinear Equations, Groningen, Wolters-Noordhoff, 1974.

16

Модели многопараметрических бифуркационных задач ...

4. B. V. Loginov, T. E. Badokina, O. V. Makeeva, “Green functions construction for divergence problems in aero-elasticity”, ROMAI Jornal, 2008, vol. 4, no. 2, pp. 33-44.

5. М. А. Наймарк, Линейные дифференциальные операторы, М.: Наука, 1969. 528 с.; M. A. Naimark, Linear differential operators, V. 1, Elementary theory of linear differential operators, London, Toronto, Harrap, 1968; V. 2, Linear differential operators in Hilbert space, London, Toronto, Harrap, 1968.

6. Б. В. Логинов, О. В. Кожевникова, “Вычисление собственных изгибных форм и асимптотики разветвляющихся решений бифуркационной задачи о дивергенции пластины”// Известия РАЕН, 1998. Т. 2, №3. С. 112-120. [B. V. Loginov, O. V. Kozhevnikova, “Computation of eigen-bending forms and branching solutions asymptotics for bifurcation problem on rectangular plate divergence”, Izvestiya RAEN,, 1998, vol. 2, no. 3, pp. 112-120. (In Russian)].

7. П. А. Вельмисов, Б. В. Логинов, “Метод групповых преобразователей и ветвление решений в двухточечных граничных задачах аэроупругости” / Материалы Между-нар. конф. “Дифференциальные уравнения и их приложения” (Саранск, 20-22 декабря 1994 г.), Саранск, 1995. С. 120-125. [P. A. Vel’misov, B. V. Loginov, “Group transformtion method and solutions branching in two-point boundary value problems of aeroelasticity”, Materialy Mezhdunar. konf. “Differentsial’nyye uravneniya i ikh prilozheniya” [Materials of the International Conference “Differential Equations and their Applications”], Saransk, Mordovia Univ., 1995, pp. 120-125. (In Russian)].

8. П. А. Вельмисов, С. В. Киреев, А. О. Кузнецов, “Устойчивость пластины в сверхзвуковом потоке газа” // Вестн. Ульяновск. гос. техн. ун-та. Сер. Естеств. науки, 1999. №1. С. 44-51. [P. A. Vel’misov, S. V. Kireyev, A. O. Kuznetsov, “Plate stability in supersonic gas flow”, Vestn. Ul’yanovsk. Gos. Tekh. Univ. Ser. Estestv. Nauki, 1999, no. 1, pp. 44-51. (In Russian)].

9. B. V. Loginov, A. V. Tsyganov, O. V. Kozhevnikova, “Strip-plate divergence as bifurcational problem with two spectral parameters”, Proceedings of International Symposium on Trends in Applications of Mathematics to Mechanics (Seeheim, Germany, August 22-28, 2004), Aachen, Shaker Verlag, 2005, pp. 235-246.

10. Т. Е. Бадокина, Б. В. Логинов, Ю. Б. Русак, “Построение асимптотики решений нелинейной краевой задачи для дифференциального уравнения четвертого порядка с двумя бифуркационными параметрами” // Известия Иркутского государственного университета. Сер. Математика, 2012. Т. 5, №1. С. 2-12. [T. E. Badokina, B. V. Loginov, Yu. B. Rusak, “Construction of the asymptotics of solutions of nonlinear boundary value problems for fourth order differential equation with two bifurcation parameters”, Izvestiya Irkutskogo gosudarstvennogo universiteta. Ser. Matematika, 2012, vol. 5, no. 1, pp. 2-12. (In Russian)].

11. Э. Камке, Справочник по обыкновенным дифференциальным уравнениям, М.: Наука, 1971. 576 с. [E. Kamke, Spravochnik po obyknovennym differentsial’nym uravneniyam [Handbook on Ordinary Differential Equations], Moscow, Nauka, 1971, 576 pp. (In Russian)]

Поступила в редакцию 03/IX/2013; в окончательном варианте — 24/I/2014; принята в печать — 21/II/2014.

17

T. E. Badokina

MSC: 34B08, 58E07

MODELS OF MULTIPARAMETER BIFURCATION PROBLEMS FOR THE FOURTH ORDER ORDINARY DIFFERENTIAL EQUATIONS

T. E. Badokina

Ogarev Mordovia State University,

68, Bol’shevistskaya st., Saransk, 430005, Russian Federation.

We consider the problem of computing the bifurcating solutions of nonlinear eigenvalue problem for an ordinary differential equation of the fourth order, describing the divergence of the elongated plate in a supersonic gas flow, compressing (extending) by external boundary stresses on the example of the boundary conditions (the left edge is rigidly fixed, the right one is free). Calculations are based on the representation of the bifurcation parameter using the roots of the characteristic equation of the corresponding linearized operator. This representation allows one to investigate the problem in a precise statement and to find the critical bifurcation surfaces and curves in the neighborhood of which the asymptotics of branching solutions is being constructed in the form of convergent series in the small parameters. The greatest difficulties arise in the study of the linearized spectral problem. Its Fredholmness is proved by constructing the corresponding Green's function and for this type of problems it is performed for the first time.

Keywords: supersonic gas flow, buckling, aeroelasticity, bifurcation, branching equation.

Received 03/IX/2013;

received in revised form 24/I/2014;

accepted 21/II/2014.

ISSN: 2310-7081 (online), 1991-8615 (print); doi: http://dx.doi.org/10.14498/vsgtu1258 © 2014 Samara State Technical University.

Citation: T. E. Badokina, “Models of Multiparameter Bifurcation Problems for the Fourth Order Ordinary Differential Equations”, Vestn. Samar. Gos. Tekhn. Univ., Ser. Fiz.-Mat. Nauki [J. Samara State Tech. Univ., Ser. Phys. & Math. Sci.], 2014, no. 1(34), pp. 9-18. doi: 10.14498/vsgtu1258. (In Russian)

Author Details: Tat’yana E. Badokina, Assistant, Dept. of Applied Mathematics, Differential Equations and Theoretical Mechanics.

E-mail address: [email protected]