DOI 10.22337/2077-9038-2018-2-137-142
Расчет конструкций переменной толщины методом наискорейшего спуска
В.В.Петров, СГТУ им. Гагарина Ю.А., Саратов Д.А.Пименов, СГТУ им. Гагарина Ю.А., Саратов Р.В.Мищенко, СГТУ им. Гагарина Ю.А., Саратов
В статье обсуждается методика расчёта элементов конструкций переменной толщины. Такие конструктивные элементы описываются дифференциальными уравнениями с переменными коэффициентами, для реализации которых необходимо иметь надёжный метод расчёта, позволяющий получить достаточно точное решение. Наиболее эффективным методом для расчёта таких конструкций является метод наискорейшего спуска, разработанный Л.В. Канторовичем. В рамках этой статьи идея метода изложена на примере решения задач изгиба балки и пластинки переменной толщины.
Приведена последовательность расчёта конструкции переменной толщины МНС на примере статически неопределимой балки, где в качестве начального приближения использовалось уравнение изгиба балки постоянного поперечного сечения. Затем этот метод был обобщён на более сложную двумерную конструкцию - пластинку переменной толщины. Была решена проблема построения начального приближения при решении дифференциального уравнения в частных производных. В качестве примера была рассмотрена квадратная в плане пластинка, шарнирно опёртая по контуру. Результаты расчёта сравнивались с результатами, полученными методом конечных разностей.
При решении конкретных задач методом наискорейшего спуска было выявлено его принципиальное отличие от прямых методов, таких как, например, Ритца-Тимошенко, Бубнова-Галёркина, которое заключается в том, что последовательные приближения при решении задач получаются не в априорно выбранной форме, а в форме, определяемой самой задачей. В МНС решение корректируется качественно в процессе реализации метода, а при решении задачи вариационными методами мы выбираем аппроксимирующую функцию и тем самым задаём конфигурацию решения.
Использование МНС позволяет получать конечные формулы для определения напряжённо деформируемого состояния конструкций переменной толщины, что позволит оперативно осуществлять их вариантное проектирование.
Ключевые слова: пластинка, метод наискорейшего спуска, метод конечных разностей, ограниченный оператор, невязка решения.
Calculation of Constructions of Variable Thickness by
Method of Steepest Descent
V.V.Petrov, SSTU, Saratov
D.A.Pimenov, SSTU, Saratov
R.V.Mishenko, SSTU, Saratov
The technique of calculation of structural elements of variable thickness is discussed in the article. Such constructive
elements are described by differential equations with variable coefficients, for the implementation of which it is necessary to have a reliable calculation method that allows obtaining a fairly accurate solution. The most effective method for calculating such structures is the steepest descent method, developed by L.V. Kantorovich. Within the framework of this article, the idea of the method is set forth in the example of solving problems of bending a beam and a plate of variable thickness.
A sequence is given for calculating the design of a variable MSD thickness using the example of a statically indeterminate beam, where the bending equation for a beam of constant cross section was used as the initial approximation. Then this method was generalized to a more complex two-dimensional construction - a plate of variable thickness. The problem of constructing the initial approximation for solving a partial differential equation was solved. As an example, we considered a square plate in plan, hinged on the contour. The results of the calculation were compared with the results obtained by the finite difference method.
In solving specific problems by the method of steepest descent, it was revealed that it differs from direct methods, such as, for example, Ritz-Timoshenko, Bubnov-Galerkin, which consists in the fact that successive approximations in solving problems are not obtained in a priori chosen form, but in the form , determined by the problem itself. In the MSD, the solution is corrected qualitatively in the course of implementing the method, and when solving the problem by variational methods, we choose the approximating function and thereby set the solution configuration.
The use of the MSD makes it possible to obtain finite formulas for determining the stress-strain state of structures of variable thickness, which will allow them to quickly implement their variant design.
Keywords: plate, method of steepest descent, finite-difference method, bounded operator, misclosure of decision.
В строительной практике широкое применение находят конструкции, которые по соображениям архитектурной выразительности или по экономическим критериям имеют переменную толщину. Такие конструктивные элементы описываются дифференциальными уравнениями с переменными коэффициентами, для реализации которых необходимо иметь надёжный метод расчёта, позволяющий получить достаточно точное решение, что полезно при вариантном про-
ектировании. Весьма эффективным представляется метод, предложенный в 1948 году Л.В. Канторовичем и названный им градиентным методом или методом наискорейшего спуска (МНС) [1]. В последующие годы был опубликован ряд работ, где строгими математическими методами была доказана состоятельность метода и выявлены его потенциальные возможности в теоретических исследованиях задач, касающихся нахождения минимума квадратичных функционалов и связанных с ним линейных проблем. Целью данной статьи является изучение возможностей этого метода при расчёте конструкций переменной толщины.
Метод был разработан Л.В. Канторовичем применительно к квадратичным функционалам, то есть к ограниченным операторам. Известно, что поиск экстремума функционала эквивалентен решению соответствующего дифференциального уравнения, то есть исследованию неограниченных операторов.
Л.В. Канторовичу принадлежит идея В-ограниченных операторов, которая состоит в том, что сложный линейный неограниченный оператор А можно ограничить простым линейным неограниченным оператором В. Л.В. Канторовичем доказана теорема: если оператор А положительно определён и В - ограничен, то есть выполняются условия:
т(Ви,и)<(Аи,и)<М(Ви,и), 0<т<М <+оо, (1)
то имеет место В-сходимость процесса наискорейшего спуска к решению исходного уравнения с быстротой геометрической прогрессии. Эта теорема доказана для случая, когда рассматриваемая задача описывается линейными дифференциальными уравнениями. Здесь и далее используется принятое обозначение скалярного произведения двух функцийД иД2 в виде (Д , Д).
В соответствии с этой теоремой возникает проблема построения симметричного, положительно полуограниченного оператора В, родственного оператору А. Область определения оператора В должна совпадать с областью определения оператора А, кроме того должно выполняться неравенство (1). Рядом авторов показано, что при решении задач строительной механики В-ограниченность оператора А обеспечивается, если упругая система, описываемая оператором В, является более жёсткой по сравнению с исходной, описываемой оператором А [2]. Поэтому оператор В для конкретной задачи можно построить исходя из инженерных соображений, а быстрота сходимости метода будет зависеть от удачного выбора оператора В.
Рассмотрим последовательность расчёта конструкции переменной толщины МНС на примере статически неопределимой балки. На первом этапе находим начальное приближение в виде решения уравнения изгиба балки постоянного поперечного сечения, которое в безразмерном виде запишется следующим образом
Лп / „ч
Г =/>(£)' (2)
здесь дифференциальный оператор, где,/0 посто-
янная по длине жёсткость балки; ы0 - начальное приближение прогиба; р(£) - заданная поперечная нагрузка; £ = х / Ь - безразмерная координата, где Ь - длина балки. Подставляя начальное приближение в уравнение изгиба балки переменной толщины, определяем невязку решения Д£):
Р(£) = Аи0-р{{) =
с?_
-р(*)> (3)
где Л= ,„2
- исходный неограниченный опера-
тор рассматриваемой задачи. В предположении, что ось координат £ находится в верхней плоскости балки, переменную жёсткость определяем по формуле:
г о г
№ = к |
->,(4)1
Ьёг = Е—кНг), 12 у '
Ви0 -J0
где Е - модуль упругости материала балки.
На втором этапе находим корректирующую функцию Z1(£) из дифференциального уравнения
На третьем этапе вычисляем величину градиента спуска £1 по формуле (вг г)
где скалярные произведения (В1г,1- энергия оператора В, и (А1^ Z1) - энергия оператора А, определяются по формулам:
(б) (7)
Прогиб в первом приближении МНС с учётом поправки решения запишется следующим образом:
и\ (£) ~ ио (£) ~ (£) • (8)
Далее процесс итераций повторяется до достижения необходимой точности решения по формуле:
Принципиальное отличие метода наискорейшего спуска от прямых методов, например, Ритца-Тимошенко, Бубнова-Галёркина заключается в том, что последовательные приближения получаются не в априорно выбранной форме, а в форме, определяемой самой задачей. При решении задачи вариационными методами мы выбираем аппроксимирующую функцию и тем самым задаём конфигурацию решения. В МНС решение качественно корректируется при реализации процедуры метода.
В качестве примера рассмотрим задачу изгиба балки под действием равномерно распределённой нагрузки р0, представленную на рисунке 1, где нижнее очертание балки представляет собой квадратную параболу.
Рис. 1. Балка переменной толщины
Уравнение безразмерной толщины балки имеет вид:
где Я =
К - К
■ безразмерный параметр относительной тол-
Решаем дифференциальное уравнение (4) и находим корректирующую функцию 1г с точностью до параметра относительной толщины А:
Рассмотрим изложенный алгоритм расчёта применительно к двухмерной задаче - изгибу пластинки переменной толщины [3]. Так как эта задача описывается уравнением в частных производных, то трудности возникают уже на начальном этапе - построении начального приближения. В качестве примера рассмотрим квадратную в плане пластинку, шарнирно опёртую по контуру. Нагрузка равномерно распределённая. Уравнение изгиба пластинки постоянной толщины в безразмерном виде запишем следующим образом
где приняты следующие безразмерные величины:
<Э4
ь а ' Ъ д
э4
д4
щины, с помощью которого можно управлять напряжённо деформированным состоянием конструкции. В качестве примера рассмотрим балку, один конец которой шарнирно опёрт, а другой жёстко защемлён. В этом случае решение уравнения (2) имеет следующий вид:
По формуле (3) определяем невязку решения:
где а = Ь - размеры пластинки в плане; р10
, ¿>0 = 12(1
_ ,.2\Я0а ' £Уг4
Вычисляем энергию оператора B по формуле (6):
Вычисляем энергию оператора А по формуле (7)
(Аг„г1) = р20/,(л) ю-3, /<(Я) = -0,6Л' + 5,7Я8 -24,4Л7 + 60,6/1® -93,5Л3 + 89,6Я4 -49,5Я3 +12,4Я2.
Далее по формуле (5) вычисляем величину градиента спуска и записываем прогиб балки переменной толщины по формуле (8).
На рисунке 2 представлены эпюры изгибающих моментов в рассмотренных балках переменной толщины. Результаты расчёта, полученные в третьем приближении МНС, сравнивались с результатами, полученными МКР при разбиении балки на 512 частей. Результаты практически совпали.
Результаты расчётов показали, что с увеличением параметра А происходит существенное перераспределение величины изгибающих моментов из центральной зоны к защемлённому краю. С увеличением стрелы подъёма в центральной части балки возникает практически безмоментная зона. Изгибающий момент в защемлении балки увеличился почти в два раза.
толщина пластинки на контуре; E - модуль упругости; P0 - безразмерная нагрузка, ^ - коэффициент Пуассона. Следовательно, оператор B будет иметь вид бигармонического оператора.
С учётом принятых обозначений безразмерных величин уравнение изгиба пластинки переменной толщины в безразмерном виде представим следующим образом:
Следовательно, оператор А, будет иметь вид:
А = У2[Л3(<Г,?7)У2], (12)
гдф переменная толщина пластинки. Вид переменной цилиндрической жёсткости зависит от способа задания переменной толщины. Если, например, переменная толщина пластинки описывается с помощью синусоидального веларо-ида [4], то функция h(£, щ) имеет вид:
г/) = 1 - Я эт(71%)бш(лт]).
Поперечное сечение пластины идентично изображённой на рисунке 1 балке переменной толщины. Для этого случая граничных условий начальное приближение прогиба можно найти методом двойных тригонометрических рядов. Для этого разложим нагрузку в ряд по синусам
(13)
где т = 1, 2, 3,..., п = 1, 2, 3,... - целые числа, ат = тп, Рп = пп, P - коэффициенты Фурье, вычисляемые по формуле:
Рис. 2. Эпюры изгибающих моментов при различных А
Ртп (14)
о о
Начальное приближение прогиба и0(£ ф ищем в виде ряда: И0(£7 ) = (15)
где Стп - амплитуда прогиба, вычисляемая по формуле:
Р....,
(16)
(21)
тсх(т2 +гг)~ Невязка решения имеет вид:
Корректирующая функция Z1 есть решение дифференциального уравнения
которое по виду совпадает с уравнением (10), но правая часть этого уравнения представляет собой «фиктивную» нагрузку, роль которой выполняет невязка решения. Это уравнение решаем методом двойных тригонометрических рядов. Для этого разложим невязку решения (17) в двойной ряд по синусам:
где Fmn коэффициенты Фурье, вычисляемые по формуле: Величина градиента спуска е1 определяется по формуле:
1 (лад)
где скалярные произведения {В1г, ZJ - энергия оператора В, и (AZ^ Z1) - энергия оператора A, определяются по формулам:
(5Z1,Z1) = |{[V4Z1]Z1^, (22)
о о
(^zi,Z,)=JJ{V2[/!3(Í^)V2Z1]Z1}^. (23)
о о
В первом приближении прогиб и0(£ ф определяется по формуле:
р 4 п
Невязка решения/7^, ц) примет вид: F (£ r¡) = Р0 (2, U3 - 5,5Л2 + 3,51Л - 3,25) sin (я£) sin (щ). Найдём корректирующую функцию Zfé, ф:
Вычисляем энергии операторов B и A:
Используя формулы (18) и (21), найдём поправку к прогибу:
8(32,8Я3 6,08л-2Л2 t 38.4Я 3,6л-2)
В первом приближении МНС прогиб пластины примет вид:
8(32,8Я3 -б,075я-2Я2 +38,4Я-3,6л-2
"i (!.'?)=
с„-
г103
sin(^)sin(^). (25)
"" л-б(16,4Я!-в,015л-А1 + 76,82-3,6л2
На рисунке 3 представлены эпюры прогибов, нормированные к единице в центре пластинки, а на рисунке 4 эпюры изгибающих моментов вдоль линии ц = 0,5 для пластинки переменной толщины во втором приближении, где результаты совпали с результатами, полученными МКР с сеткой 32*32. Сравнение кривых показывает качественные изменения, происходящие в эпюрах прогибов с ростом параметра X.
С ростом параметра X наблюдается существенное перераспределение изгибающих моментов по длине балки (см. рис. 4). Видно, что с уменьшением толщины пластинки в её центральной части формируется безмоментная зона, а в зонах, примыкающих к контуру, происходит значительное увеличение изгибающих моментов.
На рисунке 5 показана функция ф(Х), полученная в трёх приближениях МНС. Решение, полученное методом конечных разностей, совпало со вторым приближением МНС. Прогиб пластинки переменной толщины имеет вид:
и{%,Т1) = ф(Х)$ т(я£)8т(л77) (26)
0.1
0.2
0.3
0.4
I ч
0 0.2 0.4 0.6 0.8
1
Рис. 3. Эпюры прогибов, нормированных к единице при различных X
—«- -1-1-»
1 1 1 1 1 1
VV 0.9 i i
ч 1 ч\ i
0.1N sV Ч.
1 1 1
ж'' 16,4Д' - 6,ОНх'Л2 + 76,8Л - 3,6л-2
-/'010'sin^)sin(;r//). (24)
М(^0.5)-10"2
Рис. 4. Эпюры изгибающих моментов при различных X
// fi-
МНС ш fi ■ ft; tf
\ ft //.•'
л
мн СП V
м ЯС1
0.09 0.08 0.07 0.06 0.05 0.04 0.03 0.02 0.01
О
0.05 0.15 0.25 0.35 0.45 0.55 0.65 0.75 0.85 А
Рис. 5. Изменение амплитуды прогиба ф при различных величинах X
Для удобства расчётов можно аппроксимировать амплитуду прогиба пластинки ф(Х) с помощью подходящих математических формул.
Например, для получения конечной формулы можно аппроксимировать амплитуду ф(Х) полиномиальной зависимостью
Если кривую, полученную в третьем приближении МНС, аппроксимировать показательной, экспоненциальной или степенной зависимостью, то получим следующие формулы:
Использование МНС позволило получить конечную формулу (26) для определения прогиба пластинки, шарнирно опёртой по периметру, что позволит оперативно осуществлять вариантное проектирования при различных параметрах относительной толщины пластинки.
В решении балки и пластинки переменной толщины содержится параметр X, с помощью которого можно управлять
напряжённо деформированным состоянием конструкции. Принятое уравнение переменной толщины h(^), или h(^), конструкции будет влиять на величину концентраторов напряжений на опоре. Например, если использовать для описания переменной толщины половину эллипса, то можно создать плавный переход на границе конструкции, что обеспечит снижение концентраторов напряжений.
В заключении отметим следующие важные преимущества данного метода: во-первых, отпадает необходимость решать уравнения с переменными коэффициентами, а во-вторых, при применении МНС решение качественно улучшается в процессе решения, и нет необходимости создавать аппроксимирующие функции, к виду которых весьма чувствительны популярные вариационные методы.
Литература
1. КанторовичЛ.В. Функциональный анализ и прикладная математика / Л.В. Канторович // УМН. - 1948. - Т. 3. - № 6 (28). - С. 89-185.
2. ДеркачёвА.А. Общая теория метода мажорантной упругой системы / А.А. Деркачёв. - Душанбе, 1963.
3. Петров В.В. Теория расчета пластинок и оболочек / В.В. Петров. - Саратов: СГТУ им. Гагарина Ю.А., 2014. - С. 163.
4. Кривошапко С.Н. Аналитические поверхности / С.Н. Кривошапко, В.Н. Иванов, С.Н. Халаби. - М.: Наука, 2006. - С. 544.
Literatura
1. Kantorovich L.V. Funktsional'nyj analiz i prikladnaya matematika / L.V. Kantorovich // UMN. - 1948. - T. 3. - № 6 (28). - S. 89-185.
2. Derkachev A.A. Obshhaya teoriya metoda mazhorantnoj uprugoj sistemy / A.A. Derkachev. - Dushanbe, 1963.
3. Petrov V.V. Teoriya rascheta plastinok i obolochek / V.V. Petrov. - Saratov: SGTU im. Gagarina Yu.A., 2014. - S. 163.
4. Krivoshapko S.N. Analiticheskie poverhnosti / S.N. Krivoshapko, V.N. Ivanov, S.N. Halabi. -M.: Nauka, 2006. - S. 544.
Петров Владилен Васильевич, 1935 г.р. (Саратов). Доктор технических наук, профессор, академик РААСН. Заведующий кафедрой «Теория сооружений и строительные конструкции» Саратовского государственного технического университета им. Гагарина Ю.А. (410054, Саратов, ул.Политехническая, 77. СГТУ). Сфера научных интересов: разработка методов расчёта строительных конструкций с учётом нелинейных факторов. Тел.: +7 (905) 382-21-23. E-mail: vvp@sstu.ru.
Пименов Дмитрий Алексеевич, 1991 г.р. (Саратов). Аспирант кафедры «Теория сооружений и строительные конструкции» Саратовского государственного технического университета им. Гагарина Ю.А. (410054, Саратов, ул.Политехническая, 77. СГТУ). Сфера научных интересов: разработка методов расчёта строительных конструкций с учётом нелинейных факторов. Тел.: +7 (937) 632-19-12. E-mail: scorpions91@inbox.ru.
2 2018 141
Мищенко Роман Викторович, 1991 г.р. (Саратов). Аспирант кафедры «Теория сооружений и строительные конструкции» Саратовского государственного технического университета им. Гагарина Ю.А. (410054, Саратов, ул.Политехническая, 77. СГТУ). Научные интересы: разработка методов расчёта строительных конструкций с учётом нелинейных факторов. Тел.: +7 (987) 338-40-22. E-mail: roman_radon4@mail.ru.
Petrov Vladilen Vasilyevich, born in 1935. (Saratov). Academician of Russian Academy of Architecture and Construction Sciences, Doctor of Technical Sciences, Professor, Head of the Department of "Theory of Structures and Building Constructions" Saratov State Technical University by Gaga rin Yu.A. (410054, Saratov, Politechnicheskaya st., 77. SSTU). Scientific interests: development of methods for calculating building structures taking into account nonlinear factors. The address: 410031, Saratov, Naberezhnaya Kosmonavtov, 8, apt. 11. Phone: +7 (905) 382-21-23. E-mail: vvp@sstu.ru.
Pimenov Dmitry Alekseyevich, born in 1991. (Saratov). Post-graduate student of the chair "Theory of Structures and Building Constructions" of the Saratov State Technical University by Gagarin Yu.A. (410054, Saratov, Politechnicheskaya st., 77. SSTU). Scientific interests: development of methods for calculating building structures taking into account nonlinear factors. Address: 413102, Engels, Komarova, 29. Telephone: +7 (937) 632-19-12. E-mail: scorpions91@inbox.ru.
Mischenko Roman Viktorovich, born in 1991. (Saratov). Post-graduate student of the chair "Theory of Structures and Building Constructions" of the Saratov State Technical University by Gagarin Yu.A. (410054, Saratov, Politechnicheskaya st., 77. SSTU). Scientific interests: development of methods for calculating building structures taking into account nonlinear factors. Address: Saratov, Rostovskaya, 38, Apt. 111. Phone: +7 (987) 338-40-22. E-mail: roman_radon4@mail.ru.