CC BY
102
УДК 539.3
АНАЛИЗ АЛГОРИТМОВ ИССЛЕДОВАНИЯ УСТОЙЧИВОСТИ ТОНКОСТЕННЫХ ОБОЛОЧЕК
В. В. Карпов
Санкт-Петербургский государственный архитектурно-строительный университет, кафедра прикладной математики и информатики E-mail: karpov.v.v@lan.spbgasu.ru
Рассматривается три варианта алгоритмов исследования устойчивости тонкостенных оболочек: алгоритм, основанный на методе Ритца и итерационных процессах; алгоритм, основанный на методе наискорейшего спуска; алгоритм, основанный на методе продолжения решения по параметру. Дается анализ результатов исследования устойчивости оболочек, полученных с использованием этих алгоритмов.
Ключевые слова: алгоритмы устойчивости оболочек, алгоритм метода Ритца, алгоритм метода наискорейшего спуска, алгоритм метода продолжения по параметру.
Analysis of Algorithms Study of Stability of Thin-Walled Shells V. V. Karpov
Saint-Petersburg State University of Architecture and Civil Engineering,
Chair of Applied Mathematics and Informatics E-mail: vvkarpov@lan.spbgasu.ru
We consider three variants of algorithms for studying the stability of thin-shell: An algorithm based on the Ritz method and iterative processes, an algorithm based on the method of steepest descent, the algorithm based on a method of extending the solution to the parameter. Analyzes the results of the study of shells produced using these algorithms.
Key words: algorithms for the stability of shells, the algorithm of the Ritz method, the algorithm is the method of steepest descent algorithm of the continuation parameter.
При исследовании устойчивости тонкостенных оболочек решаются нелинейные задачи. При решении нелинейно-упругих задач и задач ползучести задача исследования устойчивости оболочек существенно усложняется. Дифференциальные уравнения равновесия весьма громоздки, поэтому целесообразнее использовать методы минимизации функционала полной энергии деформации.
Функционал полной энергии деформации подкрепленных оболочек вращения при решении линейно-упругих задач будет иметь вид [1]
а Ь
Э = 271--Г2) I I {^)(£Х + 2^£х£у + £Уу + №\1ху + к(Фх — $1)2 + к(Фу — @2)2)+
2(І - M2)
00
+2S(єxXi + MЄxX2 + єуX2 + МєуXi + 2MiYxyXi2) + ^^ (Xi + 2MXiX2 + x2 + 4MiXi2)-
2(І - M)(Pxи + PyV + qW)}ABdxdy; (1)
E
Здесь
І ди І дА тт, І л2 І дV І дВ тт, І
a nx + TBly - kxW +2 $v єу = b ly + AB ax и - kyW +2
І 1V І ди І дА І дВлт п л
YJ.. = —-------------------------1-- и--- V + в і 02 •
lxy А дх В ду АВ ду АВ 1x '
[І 1W 1 N л [І 1W 1 ТД
01 = \Tax + kxP)' 02 = "(В~дї + kyV)'
X =І д/x + дА. X =І д/у + дВ.
XI A 1x + АВ ду/у’ Х2 В ду + АВ 1x xi
І д/у І c)/x І ґ дА / дВ /
T~3x + В ~ду - TB \3y x +
Рх, Ру, д — компоненты внешней нагрузки вдоль осей х, у, г; А, В — параметры Ляме; Е, ^ модуль упругости и коэффициент Пуассона материала оболочки.
© Карпов В. В., 2012
63
Изв. Сарат. ун-та. Нов. сер. 2012. Т. 12. Сер. Математика. Механика. Информатика, вып. 1 Высота и расположение ребер задается функцией
т п пт
Н(х,у) = Е Ь? 6(х — х?) + Е Ь ${У — Уг) — ЕЕ Ьг? 6(х — X? )6(у — уг).
? = 1 г=1 г=1 ? = 1
где Ьг, Ь? — высота ребер параллельных осям у и х соответственно, Ьг? = шт{Ьг,Ь?}, 5(х — х?), д(у — уг) — единичные столбчатые функции, равные единице в местах присоединения ребер и равные нулю вне таких мест. Контакт ребра и обшивки происходит по полосе.
Таким образом, толщина всей конструкции равна Ь + Н. Если Н > 0, то оболочка подкреплена ребрами или накладками, а если Н < 0, то она ослаблена вырезами;
Н/2+И Н/2+И
Ё = ! Лг = Н, £ = ! гЛг = Н (Ь + Н),
Ь/2 Н/2
Н/2+И
1 = J г2 Лг = 0.25Ь2Н + 0.5ЬН2 + 1Н3.
Ъ/2
Наиболее простой в реализации алгоритм исследования устойчивости оболочки основан на методе Ритца. При этом получается система нелинейных алгебраических уравнений. Для решения таких систем применяется метод итераций, когда все нелинейные члены переносятся в правую часть.
Искомые функции перемещений и(х,у), V(х,у), Ш(х,у) и углов поворота нормали фх (х,у), 4у (х,у) представляем в виде
и (х,у) = £ и (I )Х ЦТ )¥ 1(1), V (х.у) = £ V (I )Х 2(1 )У 2(1), Ш (х.у) = £ Ш (I )Х 3(1 )У3(Т),
1=1 1=1 1=1
пп
Фх: (х,у) = ^2 р3 (I )Х 4(1 )У4(Ї), 4’у (х,у) = 52 ™ (Т)Х5(Т)У5(Т), (2)
1 = 1 1=1
где и(I), V(I), Ш(I), РБ(I), РМ(I) — неизвестные числовые параметры или функции переменной £ при решении задач ползучести; X1 — X5 — аппроксимирующие функции переменной х, удовлетворяющие при х = 0 и х = а заданным краевым условиям; У1 — У5 — аппроксимирующие функции переменной у, удовлетворяющие краевым условиям при у = 0, у = Ь.
В соответствии с методом Ритца, подставив (2) в (1), получим систему нелинейных алгебраических уравнений
дЭ =0. ™ =0, ^ = 0, ^ = 0, ^ = 0,
ди(1) дV(l) дШ(1) дРБ(1) дРМ(1)
которую кратко можно записать в виде
^л(Х ) — срР = ^и(Х).
Здесь X = (и(I).V(I),Ш(I),РБ(I),РМ(I))т, ^Л(X) — линейная часть, ^и(Х) — нелинейная часть уравнений; срР — нагрузочный член.
При последовательном увеличении нагрузки д сначала решается геометрически линейная задача
^л (X) — срР = 0,
затем методом итерации находится решение геометрически нелинейной задачи
*Л X) — срР = ^и (Хг-1). (3)
Процесс итераций заканчивается, когда прогибы на предыдущей и последующей итерациях практически совпадают. Нагрузка, при которой процесс итераций расходится, принимается за критическую нагрузку.
В качестве примера рассматриваемого алгоритма для исследования устойчивости подкрепленных пологих оболочек рассмотрим три варианта оболочек:
1) а = Ь = 600, к5 = кп = 953.64;
2) а = Ь = 600, к5 = кп = 476.82;
3) а = Ь = 600, к5 = кп = 238.41.
Здесь а = а/Ь, Ь = Ь/Ь, к^ = а2кх/Ь, кп = Ь2ку/Ь.
______________________ Ш — а4 д
Прогиб в безразмерном виде и безразмерная нагрузка будут иметь вид Ш = —, Р =
Ь
Оболочки могут быть подкреплены 18 или 36 ребрами высотой 3Ь и шириной 2Ь.
На рис. 1, а-в представлены графики «нагрузка Р - прогиб Ш» в центре оболочки (кривая с
индексом 1, кривая с индексом 2 - это прогиб Ш(0.25,0.25)). Кривые с номером 1 соответствуют оболочке без ребер, 2 - оболочке, подкрепленной 18 ребрами, 3 - оболочке, подкрепленной 36 ребрами. Расчеты проведены при N = 9.
Как видно из рис. 1, а-в, по рассмотренной методике находится только первая из возможных критических нагрузок. При этой нагрузке процесс итерации расходится, и закритическое поведение конструкции исследовать не удается. Следует отметить, что такая методика удобна и при решении нелинейно-упругих задач (физически нелинейных) и задач ползучести.
ЕЬ4'
б в Рис. 1. Кривые зависимости «Р — Ш» для вариантов оболочки 1 (а); 2 (б); 3 (в), подкрепленной различным
числом ребер
Тонкостенные оболочки могут иметь несколько форм потери устойчивости, когда прохлопывается только часть оболочки. Причем такое прохлопывание может происходить несколько раз перед тем, как оболочка прохлопнется полностью. Для исследования местных и общих форм потери устойчивости необходимо иметь возможность исследовать поведение конструкции и в закритической области. Такую возможность дает алгоритм, основанный на методе наискорейшего спуска [2, 3].
Если разложения искомых функций в ряды (2) подставить в функционал (1), то получим нелинейную функцию от неизвестных параметров и(I), V(I), Ш(I), РБ(I), РМ(I) или кратко Ф(Х,д). При каждом значении нагрузки д находится минимум этой функции по методу наискорейшего спуска,
Изв. Сарат. ун-та. Нов. сер. 2012. Т. 12. Сер. Математика. Механика. Информатика, вып. 1
когда последовательно уточняется значение Х, начиная с некоторого начального значения Х0
д Ф(Хг-1,д)
Хг — Хг-1 — Іг
дХ
Если алгоритм, основанный на методе Ритца и итерационных процессах, требует многократного решения систем линейных алгебраических уравнений с изменяющейся в процессе итераций правой частью, что занимает достаточно большое время расчета на ЭВМ, то алгоритм, основанный на методе наискорейшего спуска, свободен от этого недостатка. Некоторая сложность возникает с нахождением параметра іг на каждой итерации, но эта задача решается методом золотого сечения.
На рис. 2, а-в представлены графики «нагрузка д - прогиб Ш» оболочки в характерных точках (в центре, четверти и наибольший прогиб) для панелей цилиндрических оболочек. Рассматривается несколько вариантов стальных цилиндрических панелей с радиусом г — 5.4 м и толщиной Ь — 0.01 м, разными размерами а вдоль оси х и углами разворота ук (рад) вдоль оси у.
МПа
Итах -О- Щ0.5,0.5) Щ0.25,0.25)
Рис. 2. График «нагрузка д - прогиб Ш»: а = 10 м, ук = 1.57 (а); а = 10 м, ук = 3.14 (б); - а = 20 м, ук = 1.57 (в)
б
Как видно из рис. 2, а, б, для панели с параметром а = 10 м происходит один раз перескок в новое равновесное состояние и после потери устойчивости прогиб оболочки в различных точках начинает монотонно возрастать. Для панели с параметром а = 20 м при у^ = 1.57 происходит несколько перескоков в новое равновесное состояние (несколько потерь устойчивости) после чего прогибы в разных точках начинают монотонно возрастать.
Таким образом, этот алгоритм позволяет находить как местные, так и общие формы потери устойчивости оболочек, но не позволяет находить нижние критические нагрузки и точки бифуркации. Эта методика также позволяет решать как физически нелинейные задачи, так и задачи ползучести.
Решения нелинейных задач не единственно. Существует при одной и той же нагрузке несколько форм равновесных состояний. Существуют точки пересечения кривых «нагрузка - прогиб», которые называют точками бифуркации. Вблизи таких точек может наступить перескок с одной формы равновесия на другую. Поэтому необходимо при исследовании устойчивости оболочек находить не только критические нагрузки, но и точки бифуркации.
Алгоритм, основанный на методе продолжения решения по параметру, позволяет находить как верхние, так и нижние критические точки кривой «нагрузка - прогиб» и, кроме того, находить точки бифуркации [2,4].
Уравнения равновесия (3) кратко можно еще записать в виде
Г (Х,Р) = 0.
Продифференцировав это уравнение по параметру Р получим
дГ йХ дГ п
зХТр + дР = 0■ (4)
йХ
Уравнение (4)) линейное относительно -^р. Если при Р = Р0 известно X = Х0, то, используя для решения краевой задачи для уравнения (4)) метод Эйлера, получим расчетную схему
Хг+1 = Хг + АХг, Рг+1 = Рг + АРг-,
где АРг задается, а АХг находится из линейного уравнения
Гх (Хг, Рг )АХг + Гр (Хг, Рг)АРг = 0.
Иа рис. 3 представлены графики «нагрузка Р — прогиб Ш» в центре оболочки (сплошные линии) и в четверти (пунктирные линии) для пологой оболочки с параметрами а = Ь = 120Н, = кп = 32,
подкрепленной двумя ребрами, проходящими параллельно сторонам оболочки и пересекающимися в центре (кривая 2) и 8 ребрами (по 4 в каждом направлении, кривая 8) высотой 3к и шириной 2Н. Индекс Т означает, что ребра имеют высоту 15Н и ширину 3Л,. Для оболочки, подкрепленной 8 ребрами, происходит общая потеря устойчивости при Р = 1750. При этом прогиб при потере устойчивости составляет около двух толщин оболочки, т. е. половину от величины стрелы подъема оболочки (для рассматриваемой оболочки стрела подъема составляет АН). Для оболочки, подкрепленной 2 ребрами, вначале наступает местная потеря устойчивости при Р = 800 (прохлопываются четыре панели между ребрами и краем оболочки, точка М переходит в точку М1, а точка N — в точку М1), затем при Р = 1185 происходит общая потеря устойчивости (оболочка прохлопывается полностью). Для оболочки, подкрепленной более жесткими ребрами (кривая с индексом Т), местная потеря устойчивости происходит, а общей потери устойчивости может не произойти. Путем анализа знака определителя матрицы Якоби 1 можно определить точки бифуркации. В точках бифуркации определитель этой матрицы меняет знак.
Методика, основанная на методе продолжения решения по параметру, сложна для решения физически нелинейных задач и задач ползучести.
Изв. Сарат. ун-та. Нов. сер. 2012. Т. 12. Сер. Математика. Механика. Информатика, вып. 1
Рис. 3. Зависимость нагрузки от прогиба при различных формах потери устойчивости
Таким образом, анализируя рассмотренные алгоритмы решения задач устойчивости тонкостенных оболочек, можно сделать следующие выводы.
1. Алгоритм, основанный на методе Ритца и итерационных процессах, наиболее прост в реализации на ЭВМ. Он удобен при расчете напряженно-деформированного состояния (НДС) строительных конструкций, так как может быть применен с любого значения нагрузки. Кроме того, конструкции, применяемые в строительстве, могут работать при нагрузках, не только меньших первой критической нагрузки, но и до нагрузок, при которых наступают пластические деформации, или деформации ползучести. Применение такого алгоритма показано в работах автора и его учеников [1,2,5].
2. Алгоритм, основанный на методе наискорейшего спуска, позволяющий определить местные и общие формы потери устойчивости, впервые был применен в исследовании устойчивости подкрепленных оболочек вращения в работах [2, 3]. Процесс нахождения НДС оболочки может начинаться с любого значения нагрузки. Этот алгоритм удобен в тех задачах, где требуется знать весь процесс последовательной потери устойчивости оболочки.
3. В дополнение к возможностям перечисленных выше алгоритмов скажем об алгоритме, основанном на методе продолжения решения по параметру и позволяющем находить еще нижние критические нагрузки и точки бифуркации. Такой алгоритм применялся в работах автора и его учеников [1,2,4], а также в работах [6, 7].
4. Алгоритмы, основанные на методе Ритца и итерационных процессах и методе наискорейшего спуска при соответствующем подборе числа разложений искомых функций в ряды, дают решение с достаточно высокой точностью. Точность алгоритма, основанного на методе продолжения решения по параметру, зависит не только от числа функций разложений искомых функций в ряды, но и от величины шага нагружения Ад. Для повышения точности расчетов необходимо или брать достаточно малое значение Ад, если для решения линеаризированных уравнений применяется метод Эйлера, или использовать более точные методы для линеаризации уравнений.
Библиографический список
1. Карпов В. В. Математическое моделирование, алго- де наискорейшего спуска // Математическое модели-
ритмы исследования модели, вычислительный экспери мент в теории оболочек. СПб., 2006. 330 с.
2. Карпов В. В. Прочность и устойчивость подкреплен ных оболочек : в 2 ч. Ч. 1: Модели и алгоритмы ис- 4. Карпов В. В., Игнатьев О. В., Сальников А. Ю. следования прочности и устойчивости подкрепленных Нелинейные математические модели деформирования
рование и краевые задачи : тр. седьмой Всерос. науч. конф. с междунар. участием : в 2ч. Ч. 1. Самара, 2010. С. 47-51.
оболочек вращения. М., 2010. 286 с.
3. Баранова Д. А., Карпов В. В. Алгоритмы исследования устойчивости оболочек, основанные на мето-
оболочек переменной толщины и алгоритмы их исследования. М.; СПб., 2002. 420 с.
5. Карпов В. В., Кудрявцев В. К. Устойчивость ребри-
стых пологих оболочек при длительном нагружении // VII Всесоюз. конф. по теории пластин и оболочек. М., Вестн. ВолгГАСУ. Сер. Строительство и архитектура. 1970. С. 276-279.
2006. Вып. 6 (21). С. 53-57. 7. Карпов В. В., Петров В. В. Уточнение решений при
6. Карпов В. В., Филатов В. Н. Закритические дефор- исследовании шаговых методов в теории гибких пла-мации гибких пластин в температурном поле с учетом стин и оболочек // Изв. АН СССР. МТТ. 1975. № 5. изменения свойств материала от температуры // Тр. С. 189-191.
УДК 533.6
ПРИМЕНЕНИЕ ОБОБЩЕННЫХ ФУНКЦИЙ В ДИНАМИЧЕСКИХ КОНТАКТНЫХ ЗАДАЧАХ АЭРОУПРУГОСТИ КРЫЛА
Е. П. Лукащик
Кубанский государственный университет, Краснодар, кафедра информационных технологий E-mail: lep_9091 @mail.ru
В настоящей работе рассматривается круг вопросов, связанных с анализом деформированного состояния тонкого упругого крыла при колебаниях в ограниченной акустической среде. В качестве аппарата математических исследований выбрана теория обобщенных функций. Результаты проведенного численного эксперимента подтвердили наличие в акустической среде демпфирующих сил и резонансных явлений, обусловленных гибкими свойствами крыла.
Ключевые слова: ограниченная акустическая среда, тонкое крыло, обобщенные функции, демпфирующие силы, динамическая неустойчивость.
Application of Generalized Functions in Dynamical Contact Problems of Wing Aeroelasaticity
Е. P. Lukaschik
Kuban State University, Krasnodar,
Chair of Information Technology E-mail: lep_9091@mail.ru
The range of problems dealing with analysis of deformed state of thin elastic wing at the oscillations in limited acoustic medium is considered in this article. The theory of generalized functions was chosen as an instrument for the mathematical research. By results of performed numerical experiment the existence of damp forces in the acoustic medium and resonance effects caused by elastic properties of the wing was confirmed.
Keywords: limited acoustic media, thin wing, generalized functions, damp forces, dynamical instability.
В классической аэроупругости представлены исследования в основном по сравнительно жестким, слабо деформируемым крыльям. Однако в последнее время и в технике, и в научных разработках значительное внимание уделяется упругим несущим поверхностям. Увеличивается диапазон рабочих высот для летательных аппаратов, чаще используются малые высоты, где особенно существенно влияние турбулентности атмосферы. Задачи контроля и управления за реальными процессами аэроупругости становятся предметом обсуждения на международных научных школах, семинарах по гидродинамике больших скоростей, по моделям механики сплошной среды и др. Во многих работах ставится цель найти объективно существующие аналитические зависимости, которые могут представлять интерес с точки зрения управления аэродинамическими процессами [1]. Однако стремление учесть известные свойства материалов требует комплексных исследований явлений, зачастую с привлечением сведений из различных смежных областей.
Решение рассмотренной в данной работе контактной задачи аэроупругости поднимает две взаимно связанные проблемы. Первая из них касается определения внешних гидродинамических сил (при некотором заданном движении поверхности упругого тела), а вторая — оценки деформаций, возникающих в конструкции при уже «известной» нагрузке (выраженной через гидродинамические нагрузки и зависящей от параметров самой конструкции). Для связанной задачи аэродинамики рассмотрена математическая модель, основанная на волновом уравнении, описывающем малые колебания жидкости, и на линейном уравнении изгиба тонких пластин при наличии заданных постоянных усилий в срединной плоскости. Такое приближение позволяет несколько упростить анализ сложных механических явлений, обусловленных взаимодействием аэродинамических и упругих сил, в некоторых
© Лукащик Е. П., 2G12
69
