CC BY
55
УДК 539.3:517.95:534.1
НЕСТАЦИОНАРНЫЕ ЗАДАЧИ ДИНАМИКИ ПЛАСТИН И ЦИЛИНДРИЧЕСКИХ ОБОЛОЧЕК ВРАЩЕНИЯ СТУПЕНЧАТОГО СЕЧЕНИЯ
Ю. П. Дьяченко1, Э. Я. Еленицкий2, Д. В. Петров1
1 Самарский государственный университет,
443011, Самара, ул. Академика Павлова, 1.
2 Самарский государственный архитектурно-строительный университет,
443001, Самара, ул. Молодогвардейская, 194.
E-mails: dyachenkoyuri@yandex.ru, elenit@list.ru, den51413@gmail.com
Предлагается методика точного (в рамках сформулированных допущений) нестационарного динамического 'расчёта конструируемых систем ступенчатопеременной толщины. При этом расчётная схема учитывает смещения срединных поверхностей сопрягаемых элементов. Аналитическое решение иллюстрируется результатами расчётов фундаментной части водонапорной плотины ГЭС при действии гидравлического удара и динамической реакцией составной круговой оболочки вращения при её импульсном нагружении.
Ключевые слова: нестационарные задачи, метод начальных параметров, конечное интегральное преобразование, плита переменного сечения, цилиндрическая оболочка.
Введение. Одним из эффективных аналитических методов решения нестационарных задач динамики для составных конструкций в замкнутой форме является метод непосредственного интегрирования и, в частности, его реализация в форме метода начальных параметров. Важное теоретическое и практическое значение имеет случай, когда рассматривая конструкция имеет переменную толщину, изменяющуюся непрерывно или ступенчато. Однако возможность применения такого аналитического подхода для решения нестационарных задач связана, в первую очередь, с необходимостью обеспечения условия ортогональности разложений, получаемых в результате разделения переменных. Другая особенность применения метода начальных параметров заключается в том, что в процессе интегрирования дифференциальных уравнений движения отдельных элементов системы получаются общие решения, содержащие как быстро возрастающие, так и быстро убывающие функции. В результате разрешающая система алгебраических уравнений для определения собственных значений оказывается плохо обусловленной. Отмеченные трудности удается устранить, если использовать метод ортогонализа-ции С. К. Годунова и различные варианты метода прогонки [1-3]. При этом применение этих методов ограничивалось исследованием только свободных колебаний симметричных систем.
Авторами предлагается подход, позволяющий получать в замкнутой форме решения нестационарных задач динамики пластин и оболочек вращения ступенчато-переменной толщины с конечной сдвиговой жёсткостью (теория
Дьяченко Юрий Петрович (к.ф.-м.н., доцент), доцент, каф. математического моделирования в механике. Еленицкий Эдуард Яковлевич (к.т.н., доцент), доцент, каф. сопротивления материалов и строительной механики. Петров Денис Владиславович, студент, каф. математического моделирования в механике.
типа С. П. Тимошенко). В отличие от традиционной схемы применения метода начальных параметров (когда начальные параметры имеют физический смысл перемещений и усилий), здесь в качестве последних используются произвольные константы одного из участков плиты (или оболочки) ступенчатопеременной толщины. Так же как и в [4], срединные плоскости соединяемых участков (пластин или оболочек) могут быть смещены относительно друг друга. При этом в математической модели учитывается механизм взаимного влияния компонент напряжённо-деформированного изгибного и мембранного состояний элементов. В предлагаемом алгоритме используются точные решения динамических задач для гладкой пластины и гладкой цилиндрической оболочки с конечной сдвиговой жёсткостью [5,6].
1. Математическая формулировка задачи. Математическая формулировка задачи для составной конструкции аналогична задаче для отдельного её элемента. Эта постановка должна содержать дифференциальные уравнения движения каждого участка конструкции, начальные и граничные условия. При этом последние должны содержать как краевые условия на внешних кромках составной конструкции, так и условия сопряжения её участков.
Для иллюстрации предлагаемой методики расчёта рассмотрим два примера конструируемых систем.
Сначала рассмотрим плиту ступенчатого сечения, состоящую из п элементов и имеющую шарнирное опирание на краях х = 0, х = I и произвольное— на гранях у = 0, у = Ь. Считаем, что конструкция лежит на упругом основании с коэффициентом постели 7 и _у-тый элемент плиты загружен произвольными динамическими нагрузками Р^{х,у,1) (рис. 1).
Напряжённо-деформирование состояние участков плиты определяются векторами
Е](х,у,1) = [Uj,Vj,Wj,OíXj,Oíyj\ , Р^{х,у,1) = [NXj, Sj, Qxj, Mxj, ,
&Лх,У,$ = 0з, ^}Т, (1)
где Fj — соответственно вектор-функции перемещений и усилий, отнесённых к срединной поверхности ^’-того участка, ] = 1,2,..., п.
Математическая формулировка задачи для составной конструкции при нулевых начальных условиях в стандартной форме имеет вид
где Л.,- — матрица дифференциальных операторов в частных производных, Gj—диагональная матрица инерционных коэффициентов. Элементы этих матриц приведены в работе [5]. Равенства (4) представляют кинематические и статические условия сопряжения смежных элементов, взаимное смещение которых учитывается матрицей В^.
В качестве второго примера рассмотрим составную цилиндрическую круговую оболочку вращения ступенчатого сечения (рис. 2), состоящую из п элементов. Считаем, что е-тый участок оболочки загружен произвольными осесимметричными динамическими нагрузками Ре(х,1). Напряжённо-деформированное состояние цилиндра определяется следующими векторами:
(2)
(3)
= BjAj(x,0,t),
(4)
Уе(х, І) = ре, \Ге, Фе]Т, Ре(х, І) = [Же.т, Яех, Мех]1
Ае(х,і) = [ЙД]Т,
где ие, 1¥е, Фе — компоненты вектора перемещений Уе; Мех, С,}ех, Мех ионенты вектора усилий Ре.
ком-
(5)
х
М -
Рис. 2. Схема цилиндрической оболочки ступенчато-переменной толщины
Математическая формулировка задачи для составной конструкции при нулевых начальных условиях имеет вид
где Ье, Се—соответственно матрица линейных дифференциальных операторов в частных производных по ж и диагональная матрица инерционных коэффициентов, которые определяются формулами [6]:
Здесь Ие1, Д=2, .]ег, Зе2 — соответственно жёсткость растяжения, цилиндрическая жёсткость, площадь сечения и момент инерции сечения оболочки; ге — радиус цилиндра; и, к — коэффициенты Пуассона и сдвига.
Равенства (8) представляют кинематические и статические условия сопряжения смежных элементов с номерами (е — 1) и е, взаимное смещение £е-1 которых учитывается матрицей Ве_\ [4,5].
Для замкнутой формулировки задачи к соотношениям (7)—(9) необходимо добавить условия, соответствующие конкретным способам закрепления торцевых концов цилиндрической оболочки (х = О, X = Ь).
2. Построение общего решения рассматриваемых задач для произвольного динамического воздействия. Сравнивая математические формулировки этих двух задач, видим, что визуально формулы (1)—(4) и (5)—(8) идентичны. Покажем возможность интегрирования начально-краевой задачи (1)—(4) путём применения метода начальных параметров в сочетании с методом конечных интегральных преобразований.
Принятые условия шарнирного опирания на краях х = 0, х = I позволяют использовать синус- и косинус-преобразования Фурье по переменной х:
где Фк{х) —диагональная матрица синусов и косинусов аргумента кттх/1 = т.
В результате применения преобразования (9) к исходной начально-крае-вой задаче, соотношения (1)—(4) принимают вид
(6)
(7)
к(1 — и) д 2 дх к(1 — г/)
О 5
к=1
ЦкШ)= ЦъЩк
(Ьзк + = Р3{(у,і), (Ю)
Л>,0)=0, = 0, (11
і=О
Ц_1к{а].1,і) = В]Цк{0,і), (12)
[і
^¿(у^)=Уо фк{х)Рз{х,у,і)(Іх.
Аналогичная процедура должна быть применена также для соответствующих граничных условий на рёбрах у = ао для первого участка, у = ап — для последнего участка конструкции.
С целью дальнейшего разделения переменных задачи (10)—(12) введём на сегментах [0, а^\ конечное интегральное преобразование по переменной у:
<р(Ьы,$ = '*Г, в]ы(хы,у)ОзЩк(у,г)<1у, (13)
і=і С©
^к(у,і) = ^^(А^,£)Л^(Аь,у)||^|Г2. (14)
г=1
Здесь выражение (13)—прямое преобразование, (14)—формула обращения, справедливая при выполнении обобщённого условия ортогональности:
№'тт п ч, (\\ZikW~2, при г = т;
.Jo ^ 0, при г ф т.
Применяя преобразование (13) к уравнениям (10) и начальным условиям (11), приходим к зависимостям
/ В]ы[Ьзк + О^—^б^у = р(\кг,г), (15)
3=1 -]0
<^(Аь,0) = 0, ф{ Аь,^)|4=0 = о, (16)
где р(Хкг^) — трансформанта нагрузки, определяемая формулой
га ^ г а, з -> ^
р(Хкг^) = ^2 В1ыЩк(У^)(1У-
з=1
Использование двух условий структурного алгоритма метода конечных интегральных преобразований [5] приводит зависимость (15) к следующему виду:
ф{Хы^) + \2ы1р{\ы,{) = р{\ы^), г = 1,2,.... (17)
Уравнение (17) с учётом начальных условий (16) представляет счётную задачу Коши для каждого тона колебаний, решение которой записывается в форме определённого интеграла.
Использование особого приёма в преобразовании выражений для внеин-тегральных членов, появляющихся при применении метода конечных интегральных преобразований, позволяет получить условия, при которых их сумма обращается в нуль. Эти условия в матричной форме имеют вид
— — = ] = 2, 3,..., п. (18)
Полученные равенства (18) являются инвариантными к исходным краевым условиям (4) и обеспечивают обобщённую ортогональность разложения (14).
В процессе преобразования выражения (15) также получены дифференциальные уравнения для форм колебаний системы, которые совместно с граничными условиями и представляют ядровую задачу по определению собственных функций разложения (14). Интегралы полученных дифференциальных уравнений могут быть представлены в матричной форме:
A-jki AjkiGjki■, (19)
где — матрица общих решений однородных уравнений, — вектор-столбец произвольных постоянных. Подставляя (19) в (18) и выполняя стандартные процедуры метода начальных параметров, выражаем искомые векторы произвольного ^’-того элемента плиты через аналогичный вектор первого участка, то есть
з
П [-Вту1т^(0)] Ат_ 1 ¿.¿(ат_ 1) (20)
т= 2
Заметим, что при наличии соотношения (20), достаточно сформулировать краевые условия на границах плиты (у = 0 при j = 1, у = ап при ] = п), что соответствует десяти граничным условиям задачи и обеспечивает замкнутую форму её решения.
При применении метода разделения переменных для интегрирования на-чально-краевой задачи (5)—(8) укажем другой путь, отличный от изложенного выше. С этой целью предварительно рассмотрим свободные колебания системы, полагая Ре(х,£) = 0. Используя известные процедуры получения решения в форме стоячих волн, приходим к системе однородных дифференциальных уравнений
(Ьег - Ш2Ое)Уег{х) = 0, ¿ = 1,2,.... (21)
Соотношения (21) получены в результате представления вектора (5) в виде
00 ->
Ае(х,г) = ^Лег(ж)Т*(£), (22)
г=1
где Тг—скалярные функции времени, общие для всех участков цилиндра; Лег (ж)—функции форм перемещений и усилий свободных колебаний е-ТОГО
элемента, определяемые как линейные комбинации известных частных решений однородных уравнений (21). Подставляя разложение (22) в (6), получаем с учётом зависимостей (21) следующие соотношения:
¿(Ті + ш2гТг)СеУег = Ре(х,і). (23)
г=1
Умножим каждую строку формул (23) соответственно на [/ег> И^ег, Фег) проинтегрируем полученные равенства по длине каждого элемента и просуммируем отдельно левые и правые части. В результате получим
СЮ
= £>г(£), (24)
г=1
где
га „1е га „1е
= У\ / У^СеУегОх, р:1 У I У^Р^Х.
е=1 70 е=1 70
Применяя аналогичную процедуру к условиям (7), получим
£ад(0) = о, ¿д*/ВД|4_0 = о. (25)
г=1 г=1
Соотношение обобщённой ортогональности собственных форм колебаний конструкции в данном случае принимает вид
к2 - = 0, (26)
Левая часть равенства (26) с учётом (21) может быть представлена сле-
дующим образом:
(Ч2 - ^)Пгз = ^2 (У^еУез - У11еУег)<1х. (27)
е=1
/О
Выполняя интегрирование по частям первого слагаемого в правой части (27), получим
П _» _» -> ->
к2 -
е=1
В результате подстановки пределов интегрирования и перегруппировки слагаемых приходим к зависимостям
Лег(0) = Ве— іЛеі (¿е_ і), (28)
которые являются инвариантными по отношению к исходным краевым условиям (8) и обеспечивают ортогональность разложения (22). Подстановка равенства (26) в (24), (25) приводит к задаче Коши для каждого ¿-того тона колебаний. В результате имеем
Ті + со2Ті = Рііі,)/Яц
с однородными условиями при £ = 0:
тг(0) = 0, Тг(1) |*=0 = о (г = 1,2,...).
Таким образом, все компоненты разложения (22) определяются на основе интегрирования уравнений для линейного осциллятора и решения краевой задачи Штурма—Лиувилля (21), (28).
Частные интегралы дифференциальных уравнений (21) известны [6]. Представим эти результаты в матричной форме
Авг(х) = АеЛ(х)Сег, (29)
где Аег — частные решения (21), записанные в элементарных функциях; Cei — вектор произвольных постоянных, имеющий вид
Сег — \Ceili Cei2, Сei3, С*ег4) Сег5) Сегб]
Т
Подставим соотношения (29) в условия сопряжения (28) и выполним стандартные процедуры метода начальных параметров, связанные с представлением вектора Сei произвольной ступени цилиндра через аналогичный вектор первого участка
е
Сег= [П В^1Атг(Ъ)\с1г.
т= 2
Заметим, что при наличии соотношений (29) достаточно сформулировать шесть краевых условий на торцах оболочки, что обеспечивает замкнутую форму решения задачи.
Предлагаемая интерпретация метода начальных параметров при решении нестационарных динамических задач конструируемых систем ступенчато-пе-ременного сечения не приводит к появлению в решении функций, имеющих быстро возрастающие и быстро убывающие составляющие. В результате отпадает необходимость в процедуре ортогонализации решения, как это делается в рамках метода С. К. Годунова [2], A.A. Абрамова [3] и других авторов только для симметричных систем при их свободных колебаниях.
Из-за ограниченности объема статьи приведём результаты расчёта по предлагаемой методике для цилиндрической оболочки, представляющей собой матрицу взрывной штамповки (см. рис. 2), имеющей жёсткое закрепление по концам. Нагрузка моделировалась следующим выражением [7]:
Pe(t) = P0exp(-9t), е = 1,2,3. (30)
Расчёты выполнены в двух вариантах для следующих исходных данных: ¿1 = ¿2 = h = 1 м, h\ = hs = 0,2 м, Л-2 = 0,1 м, г\ = Гз = 0,9 м, Г2 = 0,95 м, £12 = £23 = 0,05 м, Е = 2,1 • 105 МПа, р = 7800 кг/м3, v = 0,16, к = 0,86, Р0 = 103 кН, 9 = 1200 1/с.
Вариант 1 соответствует конструкции с гладкой наружной поверхностью, вариант 2 — симметричному изменению поперечного сечения составного цилиндра. В таблице приведена низкочастотная часть спектра свободных колебаний сооружения для обоих случаев.
Спектр частот свободных колебаний конструкции
Номер тона Частота, кГц (вариант 1) Частота, кГц (вариант 2)
1 0,924 0,957
2 0,938 0,959
3 1,068 1,100
4 1,264 1,294
5 1,475 1,529
6 1,893 1,974
7 2,346 2,425
8 2,702 2,775
9 3,188 3,318
10 3,785 3,894
Сравнение результатов показывает, что учёт в расчётной схеме смещения срединных поверхностей цилиндрических оболочек (вариант 1) приводит к снижению всех собственных значений системы, которое для основного тона колебаний составляет 3,6 %. В вычислениях напряжённо-деформированного состояния конструкции при действии нестационарной нагрузки (30) удерживалось двадцать тонов колебаний, что обеспечивало точность расчётов в 1 4- 3 %.
На рис. 3 приведены формы основного тона и первого обертона колебаний системы Ше\ и \¥е2, эпюра прогибов 1¥е и перемещений вдоль образующей ие срединных поверхностей элементов, а также эпюры усилий (¡¿еХ) Мех, N(.0, построенные вдоль образующей составного цилиндра в момент времени £ = ¿1 = 4,3 • 10-4 с. Сплошным и штриховым линиям соответствуют ре-
зультаты первого и второго вариантов расчёта. Как видно из первых двух графиков рис. 3, изгибные формы первых двух тонов колебаний системы для сравнительных случаев меняются местами. Поскольку все формы нормальных перемещений не имеют скачков на границах участков, на эпюре We разрывы также отсутствуют. Однако эпюра перемещений Ue существенным образом зависит от наличия смещений срединных поверхностей элементов. Эта особенность приводит к повышению уровня напряжённого состояния составной оболочки. Так, например, увеличение максимальных значений кольцевых усилий для сравниваемых вариантов составляет 13,25%. Отмеченное обстоятельство следует учитывать при проектировании конструкции.
Анализ осциллограмм усилий и перемещений показал, что закон изменения этих функций во времени соответствует закону изменения давления на оболочку. После снятия нагрузки её уровень напряжённо-деформированного состояния значительно падает.
В заключение следует отметить, что описанные в статье подходы имеют достаточно универсальный характер и могут быть использованы для нестационарного динамического расчёта различных стержневых, пластинчатых и оболочечных систем ступенчатого сечения.
Авторы выражают благодарность за поддержку на конференции «Математическая физика и её приложения - 2010» лабораторией математической физики СамГУ, грантами АВЦП 3341 и 10854 и контрактом ФЦП 2173.
БИБЛИОГРАФИЧЕСКИЙ СПИСОК
1. Бидерман В. Л. Прикладная теория механических колебаний. М.: Высш. шк., 1972. 416 с. [Biderman V. L. Applied Theory of mechanical vibrations. Moscow: Vyssh. shk., 1972. 416 pp.]
2. Годунов С. К. О численном решении краевых задач для систем линейных обыкновенных дифференциальных уравнений// УМН, 1961. Т. 16, №3(99). С. 171-174. [Godunov S. К. Numerical solution of boundary-value problems for systems of linear ordinary differential equations// Uspekhi Mat. Nauk, 1961. Vol. 16, no. 3(99). Pp. 171-174].
3. Абрамов А. А. О переносе граничных условий для систем линейных дифференциальных уравнений (вариант метода прогонки) // Ж. вычисл. мат. и мат. физики, 1961. Т. 1, №3. С. 542-545; англ. пер.: Abramov A. A. On the transfer of boundary conditions for systems of ordinary linear differential equations (a variant of the dispersive method) // U.S.S.R. Comput. Math. Math. Phys., 1962. Vol. 1, no. 3. Pp. 617-622.
4. Еленицкий Э. Я., Дьяченко Ю.П. Свободные колебания прямоугольной пластины ступенчатого сечения с конечной сдвиговой жёсткостью/ В сб.: Задачи со свободными границами и нелокальные задачи для нелинейных параболических уравнений; ред. Ю. А. Митропольский, А. А. Березовский. Киев: Ин-т математики НАН Украины, 1996. С. 17-20. [Elenitskiy Е. Ya., D’yachenko Yu. P. Free vibrations of a rectangular plate stepped section with finite shear stiffness / In: Problems with, free boundaries and nonlocal problems for nonlinear parabolic equations; eds. Yu. A. Mitropol’skiy, A. A. Berezovskiy. Kiev: In-t Matematiki NAN Ukrainy, 1996. Pp. 17-20].
5. Дьяченко Ю. П. Метод расчёта нестационарного воздействия на прямоугольные пластины ступенчато-переменной толщины // Вестн. Сам. гос. ун-та. Естественнонаучн. сер., 2008. №2(61). С. 136-159. [D’yachenko Yu. P. A method for computing the unsteady action on rectangular plates of stepwise varying thickness // Vestn. Samar. Cos. Univ. Estestvennonauchn. Ser., 2008. no. 2(61). Pp. 136-159].
6. Григолюк Э.И., Селезов И. Т. Неклассические теории колебаний стержней, пластин и оболочек/ Итоги науки и техники. Механика твёрдых деформируемых тел. Т. 5. М.: ВИНИТИ, 1973. 272 с. [Grigolyuk Е. Selezov I. Т. Nonclassical theories of rod, plate, and
shell vibrations / Results in Science and Technology. Mechanics of Deformable Solids. Vol. 5. Moscow: VINITI, 1973. 272 pp.]
7. Динамический расчёт сооружений на специальные воздействия: Справочник проектировщика/ ред. Б. Г. Коренев, И. М. Рабинович. М.: Стройиздат, 1981. 215 с. [ Dynamic calculation of structures for special actions: Designer’s Handbook / eds. B. G. Korenev, I.M. Rabinovich. Moscow: Stroyizdat, 1981. 215 pp.]
Поступила в редакцию 21/XII/2010; в окончательном варианте — 22/111/2011.
MSC: 35Q74; 74К20, 74К25
NON-STATIONARY PROBLEMS OF THE DYNAMICS OF STEPPED SECTION PLATES AND ROTATION CYLINDRICAL SHELLS
Y.P. D’yachenko1, E.J. Elenitskiy2, D. V. Petrov1
1 Samara State University,
1, Academician Pavlov St., Samara, 443011, Russia.
2 Samara State University of Architecture and Civil Engineering,
194, Molodogvardeyskaya St., Samara, 443001, Russia.
E-mails: dyachenkoyuri@yandex.ru, elenit@list.ru, den51413@gmail.com
The technique of exact non-stationary dynamic calculation of compound designed systems of in steps-variable thickness a method initial parameters is offered. The settlement scheme considers displacement of median surfaces of interfaced elements. As an example calculation of base plate of a dam GES and matrices of explosive punching is resulted at pulse influence.
Keywords: non-stationary problems, method of initial parameters, plate of step section, cylindrical shells.
Original article submitted 21/XII/2010; revision submitted 22/111/2011.
Yuriy P. D’yachenko (Ph.D. (Phys. & Math.)), Associate Professor, Dept, of Mathematical Methods in Mechanics. Eduard J. Elenitskiy (Ph.D. (Techn.)), Associate Professor, Dept, of Resistance of Materials and Constructional Mechanics. Denis V. Petrov, Student, Dept, of Mathematical Methods in Mechanics.
