CC BY
22
УДК 517.9
DOI: 10.14529/mmp140212
НАЧАЛЬНО-КОНЕЧНАЯ ЗАДАЧА ДЛЯ ЛИНЕЙНОЙ СТОХАСТИЧЕСКОЙ МОДЕЛИ ХОФФА
Е.А. Солдатова
Линейная модель Хоффа, исследующая в линейном приближении динамику выпучивания двутавровых балок в конструкции, представляет собой множество линейных одномерных уравнений Хоффа, заданных на ребрах геометрического графа с условиями непрерывности и баланса потоков в его вершинах. Ранее детерминированная модель изучалась в разных аспектах многими специалистами. Стохастическая модель изучается впервые. В качестве метода исследования используется классический подход Ито - Стратоновича - Скорохода, распространенный на гильбертовы пространства и уравнения соболевского типа. Основной результат - теорема об однозначной разрешимости начально-конечной задачи с аддитивным белым шумом, под которым понимается обобщенная производная K-винеровского процесса. Решение представлено в виде формул, допускающих постановку вычислительных экспериментов.
Ключевые слова: начально-конечное условие, линейные уравнения Хоффа; стохастические уравнения соболевского типа; геометрический граф; виннеровский процесс; аддитивный белый шум.
В [1] впервые было рассмотрено линейное стохастическое уравнение соболевского типа
Ldu = Mudt + NdW (1)
с начальным условием Коши. Здесь операторы L, M, N € L(U; F), где U и F - веществен-
ные сепарабельные гильбертовы пространства, а оператор M (L^-ограничен, p € {0} U N. Между тем было отмечено, что для уравнений соболевского типа более «естественным> (по сравнению с условием Коши) являются другие начальные условия, например условие Шо-уолтера - Сидорова [2]. В [3] впервые была рассмотрена начально-конечная задача, ошибочно названная «задачей Веригина>. Впоследствии ошибка, дезориентирующая многих исследователей, была исправлена. Начально-конечное условие
Pin(u(0) - uo) = Pfin(u(r) - ur) = 0, т € R+, (2)
где Pin, Pfin € L(U) - относительно спектральные проекторы, было использовано во множестве случаев, имеющих прикладное значение (см. обзор [4]), и даже было применено в теории оптимального управления [5]. Поскольку условие (2) приобретает все большую популярность [6], то предлагается впредь именовать его «условием Свиридюка - Загребиной> ради купирования возможной путаницы.
Пусть выполняется условие на L-спектр aL(M) оператора M, гарантирующее существование проекторов Pin и Pfin
aL(M) = fin(M) U &Ln(M), причем fin(M) содержится в ограниченной области ..
D С C с кусочно гладкой границей 7, причем 7 П aL(M) = 0.
Заметим, что Pin + Pfin = P - единица разрешающей группы однородного уравнения (1). Кроме того, (SZ) гарантирует существование проекторов Qin, Qfin € L(F), которые в сумме равны единице Q группы, порожденной однородным уравнением (1) на F
В последнем слагаемом правой части (1) символом W = W (t) обозначен F-значный K-винеровский процесс
го
W(t) = Vxkвк(t)<Pk,
k=0
t Є R+, где {вк(t)} - последовательность независимых броуновских движений, [Xk} - последовательность собственных значений самосопряженного и положительно определенного ядерного оператора K Є L(U), занумерованных по невозрастанию с учетом их кратности, а [^к} - последовательность собственных векторов.
В дальнейшем считаем оператор M (L, 0)-ограниченным, а также выполненным условие
QN = N. (*)
Следуя идеологии [1], [4], выпишем «мягкое решение> (mild solution) задачи (1), (2)
u(t) = Utnuo + Sin i Ui~sL~1 QinNW(s)ds + L~lQNW(t) +
Jo
+и%ПТ + Sfin £ UfnSL- inQfinNW(s)ds - L-1nQfinNW(t).
Теорема 1. Пусть оператор M (L, 0)-ограничен и выполнены условия (SZ) и (*). Тогда для любых независимых друг от друга и от K-винеровского процесса W U-значных гауссовых случайных величин uo и иТ существует единственное решение и Є CL2 задачи (1), (2).
Пусть теперь G = G(V, E) - конечный связный ориентированный граф, где V = {Vi} -множество вершин, а E = {Ei} - множество ребер, причем каждое ребро Ej имеет длину lj Є R+ и площадь поперечного сечения dj Є R+. В вершинах V графа G зададим условия «непрерывности> и «баланса потока> соответственно
Uj(0,t) = Uk(0,t) = um(lm,t) = Un(ln,t),Ej,Ek Є Ea(Vi),Em,En Є Еш(Vi), (3)
^ djujx(0,t) - ^ dkukx(lk,t) = 0, (4)
j:Ej&Ea(Vi) k: Ek&E- (Vi)
где через Ea(w')(Vi) обозначено множество ребер с началом (концом) в вершине Vi, t Є R. Если граф состоит из одного нециклического ребра (т.е. вершин у графа две), то условие (3) отсутствует, а условие (4) превращается в условие Неймана. Если же ребро циклическое (т.е. вершина у графа одна), то условия (3), (4) превращаются в условия согласования. Заметим еще, что в контексте условий (3), (4) «отсутствовать> не значит «быть равным нулю>. Например, если в вершину Vi все ребра «входят>, то первые два равенства в (3) и уменьшаемое в (4) именно «отсутствуют>, а не равны нулю.
Теперь на графе G с условиями (3), (4) рассмотрим линейную стохастическую модель Хоффа
Xduj + dujxx = auj dt + Nj dWj, (5)
описывающую динамику выпучивания двутавровых балок, находящихся под постоянной нагрузкой в конструкции со случайным внешним воздействием. Здесь случайный процесс uj = uj(x,t), (x,t) Є (a,b) x R, характеризует отклонение j-той балки от положения равновесия; параметры Х Є R+, а Є R характеризуют нагрузку и свойства материала балок, dWj = dWj(t) отвечает случайной нагрузке на j-ю балку (аддитивный белый шум).
Чтобы редуцировать (3) - (5) к (1), введем в рассмотрение гильбертово пространство Я = {и = (иі,и2,... ,из,...) : из Є №2>(0,1з), и выполнено (3)} со скалярным произведением {и, у) = ¿з I (иЗхУЗх + из Уз )йх. Затем построим еще одно гильбертово про-
Е ее 0
странство 1^2(О) = {д = (д\, д2,..., дг,...) : дг Є ^(0,^)}. со скалярным произведением і
(д, К) = Е ¿3 дзКз¿х, причем заметим плотность и непрерывность (даже компактность) Е еЕ 0
вложения Я ^ Ь2(С). Отождествим Ь2(С) со своим сопряженным, и через $ обозначим сопряженное относительно двойственности {■, ■) пространство к Я. Очевидно, $ - гильбертово пространство, причем вложение Я ^ $ компактно.
і
Фиксируем а Є М+ и формулой {Аи, у) = ¿з / (изхУзх—аизУз)4х определим оператор
Е ее 0
А : Я ^ $. Поскольку СіЦиЦ^ < {Аи,и) < С2І|и||и при всех и Є Я и некоторых Ск Є М+, к =
1, 2, то линейный оператор А : Я ^ $ биективен и непрерывен. Отсюда по теореме Банаха следует существование оператора А-1 : $ ^ Я. Поскольку вложение Я ^ $ компактно, то оператор А-1 Є £($) является компактным. Значит, спектр оператора А вещественен, дискретен, конечнократен и сгущается только к +ж. Более того, {Аи,и) > соп8^|и||д, где 1Ы1я - норма пространства Я. Отсюда следует, что спектр &(А) С М+. Наконец, фиксируем а, X Є М и построим операторы Ь = (X — а)І + А, М = а(а1 + А). Справедлива
Лемма 1. (і) При любых а, X Є М операторы Ь, М Є £(Я; $), причем спектр <г(Ь)
оператора Ь вещественен, дискретен, конечнократен и сгущается только к +ж.
(іі) При любых а, X Є М \ {0} оператор М (Ь, 0)-ограничен.
Пусть {щ} - собственные значения оператора Ь, занумерованные по невозрастанию с учетом их кратности; а {^} - соответствующие им ортонормированные в смысле Ь2(О) функции. Построим проекторы и разрешающую группу задачи (3) - (5) соответственно
{I, если 0 Є а(Ь); ( I, если 0 Є а(Ь);
I — Е (■ ,¥к)<Рк, если 0 Є <г(Ь); Я =4 I — ^ (■ ,рк)рк, если 0 Є <г(Ь);
Vк=^-а \ Хк=Х—а
го
и * = £ 'в^( ■ ,фк )?к, (6)
к=1
где штрих у знака суммы означает отсутствие членов ряда с номерами к такими, что Рк = X + а. (Заметим, что, несмотря на «:похожесть>, проекторы Р и Я определены на разных пространствах).
Теперь у нас все готово для определения свободного члена в правой части уравнения (1) по правым частям уравнений (5). Вернемся к оператору А, и в качестве оператора К возьмем оператор Грина А-1. Его собственные значения Хк = (Vк — а)-1, занумерованные по невозрастанию, сходятся к точке нуль. Асимптотическое поведение собственных чисел {Рк — а} оператора А изучено плохо, поэтому мы нуждаемся в гипотезе
Рк ~ к2 при к ^ ж, (7)
которая заведомо выполняется, если граф О представляет собой цепочку последовательно соединенных ребер. Считая, что (7) выполнено, заключаем, что К - ядерный оператор.
Наконец, построим K-винеровский процесс
ГО
W (t) = ^2' y/\kßk(t)<Pk, (8)
k=l
где штрих означает то же самое, что и в (6). В силу леммы 1 существование этого штриха гарантирует выполнение (*).
Стало быть, все символы абстрактного решения определены и значит, справедлива
Теорема 2. Пусть а, X € R \ {0}. Тогда для любой U-значной случайной величины щ, независимой от K-винеровского процесса (8) существует единственное мягкое решение задачи (2) - (5).
В заключение объясним свободные члены уравнения (5). По построению пространство L2(G) является конечной прямой суммой попарно ортогональных подпространств. Обозначим через nj ортопроектор на j-ое подпространство. Так вот, Wj = njW, где W = W(t) (8).
Литература
1. Загребина С.А. Линейные уравнения соболевского типа с относительно р-ограниченными операторами и аддитивным белым шумом / С.А. Загребина, Е.А. Солдатова // Известия Иркутского государственного университетата. Серия: Математика. - 2013. - Т.6, №1. - С. 20-34.
2. Свиридюк Г.А., Загребина С.А. Задача Шоуолтера - Сидорова как феномен уравнений соболевского типа / Г.А. Свиридюк, С.А. Загребина // Известия Иркутского государственного университетата. Серия: Математика. - 2010. - Т.3, №1. - С. 104-125.
3. Свиридюк Г.А. Задача Веригина для линейных уравнений соболевского типа с относительно р-секториальными операторами / Г.А. Свиридюк, С.А. Загребина // Дифференциальные уравнения. - 2002. - Т.38, №12. - С. 1646-1652.
4. Загребина С.А. Начально-конечные задачи для неклассических моделей математической физики / С.А. Загребина // Вестник Южно-Уральского государственного университетата. Серия: Математическое моделирование и программирование. - Челябинск, 2013. - Т.6, №2. - С. 5-24.
5. Манакова Н.А. Оптимальное управление решениями начально-конечной задачи для линейной модели Хоффа / Н.А. Манакова, А.Г. Дыльков // Математические заметки. -
2013. - Т.94, №2. - С. 225-Ц236.
6. Замышляева А.А. Оптимальное управление решениями начально-конечной задачи для уравнения Буссинеска - Лява / А.А. Замышляева, О.Н. Цыпленкова // Вестник ЮжноУральского государственного университетата. Серия: Математическое моделирование и программирование. - Челябинск, 2012. - №5 (264), вып. 11. - C. 13-24.
Екатерина Александровна Солдатова, кафедра «Уравнения математической физики>, Южно-Уральский государственный университет (г. Челябинск, Российская Федерация), [email protected].
Поступила в редакцию 17 марта 2014 г.
Bulletin of the South Ural State University. Series "Mathematical Modelling, Programming & Computer Software",
2014, vol. 7, no. 2, pp. 124-128.
MSC 35K70, 60H30 DOI: 10.14529/mmp140212
The Initial-Final Value Problem for the Linear Stochastic Hoff Model
E.A. Soldatova, South Ural State University, Chelyabinsk, Russian Federation, [email protected].
The stochastic linear Hoff model of buckling of I-beam constructions amounts to a set of linear one-dimensional Hoff equations defined on the edges of a geometric graph with continuity and balance-of-flows conditions at its vertices. The deterministic model has been studied in various aspects by many mathematicians. We study the stochastic model for the first time. Our tool is the classical Ito - Stratonovich - Skorokhod approach extended to Hilbert spaces and Sobolev-type equations. The main result is an existence and uniqueness theorem for solutions to the initial-final value problem with additive white noise, understood as the generalized derivative of the K-Wiener process. The formulas expressing the solution are suitable for computer simulations.
Keywords: initial-final value problem; linear Hoff equations; stochastic Sobolev-type equations; geometric graph; Wiener process; additive white noise.
References
1. Zagrebina S.A., Soldatova E.A. The Linear Sobolev-Type Equations with Relatively p-Bounded Operators and Additive White Noise. The Bulletin of Irkutsk State University. Series "Mathematics", 2013, vol. 6, no. 1, pp. 20-34. (in Russian).
2. Sviridyuk G.A., Zagrebina S.A. The Showalter - Sidorov Problem as Phenomena of the Sobolev-Type Equations. The Bulletin of Irkutsk State University. Series "Mathematics", 2010, vol. 3, no. 1, pp. 51-72. (in Russian).
3. Sviridyuk G.A., Zagrebina S.A. Verigin’s Problem for Linear Equations of the Sobolev Type with Relatively p-Sectorial Operators. Differential Equations, 2002, vol. 38, no. 12, pp. 1745-1752. DOI: 10.1023/A:1023812213901
4. Zagrebina S. A. The Initial-Finite Problems for Nonclassical Models of Mathematical Physics. Bulletin of the South Ural State University. Series "Mathematical Modelling, Programming & Computer Software, 2013, vol. 6, issue 2, pp. 5-24. (in Russian)
5. Manakova N.A., Dylkov A. G. Optimal Control of the Solutions of the Initial-Finish Problem for the Linear Hoff Model. Mathematical Notes, 2013, vol. 94, issue 2, pp. 220-230. DOI: 10.1134/S0001434613070225
6. Zamyshlyaeva A.A., Tsyplenkova O.N. The Optimal Control over Solutions of the Initialfinish Value Problem for the Boussinesque-Love Equation Bulletin of the South Ural State University. Series "Mathematical Modelling, Programming & Computer Software", 2012, vol. 5 (264), issue 11, pp. 13-24. (in Russian)
Received March 17, 2014
