CC BY
41
маемых сред // Проблемы механики сплошной среды Автомодельные решения уравнений динамики идеаль-
(к 60-летию акад. В.П. Мясникова): Сб. науч. тр. Вла- ного упругопластического тела при условии пластич-
дивосток, 1996. C. 116-127. ности Треска // Журн. прикл. мех. и техн. физ. 1984.
11. Быковцев Г.И., Колокольчиков А.В., Сыгуров П.Н. № 6. С. 148-156.
УДК 539.3
ЧИСЛЕННОЕ ИССЛЕДОВАНИЕ НАПРЯЖЕННО-ДЕФОРМИРОВАННОГО СОСТОЯНИЯ В ЗАДАЧАХ ИЗГИБА ТОНКОЙ АНИЗОТРОПНОЙ ПРЯМОУГОЛЬНОЙ ПЛАСТИНКИ
П.Ф. Недорезов
Саратовский государственный университет,
кафедра математической теории упругости и биомеханики
E-mail: [email protected]
В рамках гипотез Кирхгофа рассматривается задача статического поперечного изгиба тонкой прямоугольной пластинки из анизотропного материала, у которого в каждой точке имеется одна плоскость упругой симметрии, параллельная срединной плоскости пластинки. Предполагается, что вид граничных условий вдоль каждой из сторон контура не меняется. Двумерная краевая задача для определения прогиба модифицированным методом сплайн-коллокации сводится к краевой задаче для системы обыкновенных дифференциальных уравнений, которая решается численно.
Приведены результаты числовых расчетов для двух вариантов закрепления граничного контура.
Ключевые слова: анизотропная пластинка, изгиб пластинки, напряженно-деформированное состояние, модифицированный метод сплайн-коллокации.
Numerical Study of Stress-Strain State of a Thin Anisotropic Rectangular Plate
P.F. Nedorezov
Saratov State University,
Chair of Mathematical Theory of Elasticity and Biomechanics E-mail: [email protected]
Static bending of a thin rectangular anisotropic plate is considered in the framework of Kirchhoff hypotheses. At each point of the plate there is one plane of elastic symmetry parallel to the middle plane of the plate. It is assumed that the type of boundary conditions does not change along each of the straight sides. By applying of a modified method of spline collocation the two-dimensional boundary value problem for the determination of deflection is reduced to a boundary value problem for the system of ordinary differential equations, which is solved numerically. The results of numerical calculations for two variants of the boundary conditions on the contour of the plate are presented.
Key words: anisotropic plate, bending of plate, stress-strain state, modified method of spline collocation.
1. ПОСТАНОВКА ЗАДАЧИ. ОСНОВНЫЕ УРАВНЕНИЯ И СООТНОШЕНИЯ
Рассматривается статический изгиб прямоугольной пластинки с размерами в плане а х Ь и малой толщиной Н (рис. 1) под действием поперечной нагрузки интенсивности д(х,у). Пластинка изготовлена из анизотропного материала, у которого в каждой точке имеется одна плоскость упругой симметрии, параллельная срединой плоскости. Деформации считаются малыми и подчиняются закону Гука.
Уравнение для определения прогиба — точек срединной плоскости такой пластинки в рамках классической теории Кирхгофа приводится в [1]. Это уравнение после перехода к безразмерным переменным £ = х/а, п = у/Ь и ряда преобразований для безразмерного прогиба Ш(£, п) = -(х, у)/Н может быть записано в виде
Рис. 1. Анизотропная пластинка
д4 W п д4 W п.
—- + 2Y2 + 2(Yi + Y4)c
дГ
д^3 дп
д4 W
д^2 дп2
+ 27з с3
д4 W д^дп3
+ Y5 с4
д4 W дп4
= qd^/Dh.
(1)
Здесь с = а/6, (г, ^ = 1, 2, 6) — соответствующие жесткости относительно осей х и у, = = ОцЛ,а-4, 71 = О12/Оц, 72 = 2О1в/Оп, 7э = 2026/Оц, 74 = 2Овв/Оц, 7б = ^22/Оп.
Внутренние усилия и моменты в пластинке через производные от функции Ш(ё, п) выражаются по формулам [1]
* /д2Ш 2д2Ш д2
( д2 Ш 2 д2 Ш д2
Г Ж + 75+ 73сЖц)
Му = —а^Цъ^ + 75 с^ + 73 с
п _ 2_* ( д2Ш 2д2Ш 0 д2
Нху = —^и^^ + 73с^ + 274дёд^; '
(д3Ш 1Г д3Ш . д3Ш 3 д3Ш\
М = -аОЧЖ + 1 2^ + (71 + 74)с2^ + 0-5^3с V) ,
^ Лг д3Ш . д3Ш 1г 2 д3Ш 3д3
Му = + (71 + 74)сдё^ + 1-573с + дёдП2 + ^с )
ДГ* /д3Ш 0 д3Ш . 0 , 2 д3Ш 3 д3Ш
М = -а°*Ч + 2"2 сдё2дп + (71 + 274)с2 дёдП2 + *с V
ДГ* ъ* ( д3Ш . 0 , д3Ш 0 2 д3Ш 3 д3Ш\
му = —+ Ь + 274)сдё^ + 273с2+ V) ■ (2)
Решение уравнения (1) должно быть подчинено граничным условиям, вид которых определяется условиями нагружения или деформирования контура пластинки.
Для численного определения напряженно-деформированного состояния (НДС) такой пластинки применяется модифицированный метод сплайн-коллокации (ММСК) [2], который в отличие от классического метода (КМСК) [3] позволяет решать задачи при любых внешних воздействиях на контуре. Единственное ограничение состоит в том, что вид условий вдоль каждой стороны граничного контура предполагается неизменным.
2. АППРОКСИМАЦИЯ ФУНКЦИИ Ш(ё, п) ПРИ РАЗЛИЧНЫХ ВАРИАНТАХ УСЛОВИЙ НА СТОРОНАХ ё = 0 И ё = 1
Следуя [2], будем искать функцию Ш(ё,п) в виде разложения
N+2
ш(ё,п)= £ в*,,-(ё)ш(п). (3)
, = -2
Здесь В5,, (ё) — нормализованные В — сплайны пятой степени [4], построенные на равномерной сетке ёг = (г = —5,М + 5), Нх = 1/М; нумерация сплайнов ведется по среднему узлу интервала-носителя ё,-3 < ё < ё,+3, на котором сплайн В5,,(ё) отличен от нуля. Функции Ш-к(п) и WN+к(п) должны быть подобраны так, чтобы удовлетворялись граничные условия на сторонах ё = 0 и ё = 1. Уравнения для остальных функций (п) получаются из требования, чтобы уравнение (1) при подстановке выражения (3) выполнялось в точках коллокации ё = ёг* (г = 0, N).
Рассмотрим ход решения задачи для нескольких вариантов граничных условий на стороне ё = 0: a) сторона ё = 0 деформирована заданным образом
Ш (0.0; п) = (п), дШ (д0.0; п) =©*(п). (4а)
Подстановка выражения (3) в условия (4a) приводит к системе двух линейных алгебраических уравнений, из которой следует
2
Ш-р(п) = то,р(п) + ^ аР'кШк(п) (Р = 1, 2), (5а)
к=0
где ар,к — известные коэффициенты, а то,р выражаются через -О (п) и @о(п).
b) сторона £ = 0 нагружена моментами и усилиями заданной интенсивности
Мх(0.0; п) = -а2^т0(п), N0(0.0; п) = -аО1°1р0(п)- (4Ь)
Тогда из (4Ь) и (3) получаются два дифференциальных соотношения, которые преобразуются к виду = шо,р(п) + ^ «Р,к^ + с-1 5- + с-2 5- Лр,гШ(п) (р = 1, 2) (5Ь)
с известными коэффициентами ар,к, вР,к, Лр,ки функциями то,р(п).
c) если при £ = 0 заданы смешанные граничные условия
Ш(0.0; п) = -0 (п), Мх (0.0; п) = -а2Я^тОМ, (4с)
то при подстановке разложения (3) в условия (4с) получатся алгебраическое уравнение, из которого определяется Ш-2(п), и дифференциальное соотношение для функции Ш-1(п)
Ш-2 (п) = -О (п) -
о(п) - £ В5,к(£о)Шк(п)) /В5,-2(£о),
(5с)
к=-1 ) 22
- = то(п) + к= «М^ + с =- ^Ш^(п).
Как и в предыдущих случаях, коэффициенты а1;к, ^ и функция то(п) известны. Совершенно аналогично из условий на стороне £ = 1 получаются формулы и соотношения, определяющие функции
Шм+к(п) (к = 1, 2):
а) если на стороне £ = 1 задаются прогиб и угол поворота
Ш (1.0; п) = (п), ^^ п) =©М (п), (6а)
то
Шм+р(п) = тм,Р(п) + У^ -кШм-к(п) (р = 1, 2); (7а)
к=о
Ь) при нагруженной стороне £ = 1
Мх(1.0; п) = -а2Я01 т^(п), №(1.0; п) = -аЯ^(п), (6Ь)
¿2 Шм+р / \ , ¿2 Шм-к , _1 Л й ¿Шм-г ,
¿п2 = тм,РаР,м-к ¿п2 + с вР,м-г , +
р / \ \ ^ гг 1\ - к -1 \ ^ Д и,МУ-г
----(п) , ^^рм-к + с ^ в --
к=о ¿п г=-2
2
+с-2 ^ Лр,м-гШм-г(п) (р =1, 2); (7Ь)
-2
г=-2
с) в случае смешанных граничных условий на стороне £ = 1
Ш(1.0; п) = (п), Мх(1.0; п) = -а2^ т^(п), (6с)
для функций +к(п) (к = 1, 2)имеем
Шм+2(п)= (-0(п) - £ В5,г(£м)Шг(пН /В5,м+2(£м), V г=м-2 /
+1 , , ^ -1 м+1 ТТГ , , —г— = тм(п)+ +с (п).
¿п г=м-2 ¿п 5=м-2
(7с)
2
Различные комбинации из приведенных выше выражений для функций Ш-к (п) и Щм+к (п) (к = = 1, 2) позволяют по формуле (3) получить разложения функции Щ(ё, п), удовлетворяющие разным вариантам граничных условий при ё = 0 и ё = 1.Эти разложения могут быть представлены в виде
ш (ё, п) = £ ъ (ё )Ш ( п) + м (ё, п), (8)
где функция М(ё, п) известным образом выражается через характеристики внешнего воздействия на сторонах ё = 0 и ё = 1 и В — сплайны В5,г(ё) (г = —2, — 1, N + 1, N + 2).
Отдельный интерес представляют однородные граничные условия (4а), (46), (4с) и (6а), (66), (6с). В этом случае функция Щ(ё, п) в виде (8) при значениях а'о = 0, jN = N и
2
ък (ё) = В5,к (ё) + ^ ар,к В5,-Р(ё), ъ (ё) = В5,, (ё) а = з, N — з) р=1
2
^-к(ё) = B5,N-к(ё) + Е -кB5,N+к(ё) (к = 0, 2) р=1
удовлетворяет условиям жесткого закрепления сторон ё = 0 и ё = 1.
Если ао = —2, jN = N+2, ъ (ё) = В5,, (ё) (а = ^^Ж+2), а функции Щ-Р(п) и (п) (р = 1, 2) удовлетворяют однородным уравнениям (56) и (66), то из (8) получается представление Щ(ё, п) для свободных от нагрузки сторонах ё = 0 и ё = 1. При ао = —1, ам = N + 1,
ъ-(ё) = В 1(ё ) [В5,-2 (ёо )В5>в (ё) — В5,-2(ё)В5>в(ёо)],
1
^-(ё) = -Б-ТТТ^+2(ёN)B5)N-(ё) — B5)N+2(ё)B51N)] (« = —1, 2),
В5^+2 (sN )
ъ (ё) = В5,, (ё) (а = з, N — з)
и функциях Ш-1 ( п) и +1(п), удовлетворяющих однородным уравнениям (5с) и (6с), формула (8) определяет Щ(ё, п) при свободно опертых сторонах ё = 0 и ё = 1.
В перечисленных случаях однородных условий в формуле (8) функцию М(ё, п) следует полагать равной нулю.
3. СИСТЕМЫ РАЗРЕШАЮЩИХ УРАВНЕНИЙ. ФОРМУЛИРОВКА КРАЕВЫХ ЗАДАЧ
Количество слагаемых в разложении (8) функции Щ(ё, п) и, следовательно, число неизвестных функций Щ ( п), подлежащих определению, меняется в зависимости от вида граничных условий при ё = 0 и ё = 1. Минимальное число неизвестных Щ (п), равное N + 1, получается в случае граничных условий (4а) и (6а). Максимальное число неизвестных Щ (п) в разложении (8) равно N + 5 и получается при заданных на сторонах ё = 0 и ё = 1 внешних усилиях, в случае смешанных граничных условий при ё = 0 и ё = 1 число неизвестных Щ (п) будет N + 3.
В предлагаемой методике переход от двумерной краевой задачи для функции Щ(ё, п) к одномерной задаче для неизвестных Щ (п) выполняется методом коллокации. С этой целью вводится в рассмотрение система точек коллокации ё = ёг (г = 0, N),ёN < 1 и требуется, чтобы результат подстановки (8) в (1) выполнялся вдоль прямых ё = ёг (г = 0, N). Полученную таким образом систему N + 1 обыкновенных дифференциальных уравнений четвертого порядка относительно функций Щ ( п) условно будем называть системой «А».
В случае, когда на сторонах £ = 0 и £ = 1 заданы смещения — условия (4а) и (6а), система «А» является разрешающей системой уравнений для определения всех неизвестных Ш} (п) (? = 0, N). Эта система уравнений стандартным приемом преобразуется в систему обыкновенных дифференциальных уравнений первого порядка, которая в векторной форме записывается в виде
где компонентами неизвестного вектора У(п) = (уг (п)} (я = 0,4Ж + 3) являются функции Ш} (п) (^ = 0, N) и их производные до третьего порядка включительно, а матрица А и вектор ^ известны.
Требуя, чтобы граничные условия на сторонах п = 0 и п = 1 выполнялись в точках коллокации £ = £г (я = 0, N) для значений компонентов вектор-функции У(п), получаем систему 4Ж + 4 соотношений, которую будем условно называть условиями «Б». Условия «Б» являются в рассматриваемом случае граничными условиями для вектор-функции У(п) и в векторной форме записываются в виде
где компоненты матриц Нк и векторов вк (к = 1, 2) соответствующей размерности известны.
Если при £ = 0 и £ = 1 задается распределение внешних усилий — условия (4Ь) и (6Ь), то к уравнениям системы «А» следует добавить уравнения (5Ь) и (7Ь) и из уравнений системы «А» исключить третьи и четвертые производные функций Ш-р(п) и Шм+р(п) (р = 1, 2). Выражения для этих производных получаются из (5Ь) и (7Ь) дифференцированием. Полученная система также преобразуется к виду (9), где компонентами вектора неизвестных У(п) = (уг(п)} (я = 0,4N + 11) являются функции Ш}(п) и их производные до третьего порядка включительно, а также функции Ш-р(п), Шм+Р(п) (р = 1, 2) и их первые производные. Выполнение граничных условий при п = 0 и п = 1 требуется не только в точках коллокации £ = £г (я = 0, N), но и в концевых точках отрезков £ = 0 и £ = 1. Таким образом, к условиям «Б» добавляются еще восемь условий (по два условия в каждой угловой точке). Соответствующие граничные условия для вектора У(п) по-прежнему представляются в виде (10).
В случае граничных условий смешанного типа — условия (4с) и (6с) — система уравнений «А» с добавленными к ней уравнениями (5с) и (7с), вообще говоря, составляет полную систему уравнений для функций (п) (з = -+ 1). Этой системе эквивалентна системе 4N + 6 обыкновенных дифференциальных уравнений первого порядка, для однозначного решения которой необходимо иметь 4N + 6 граничных условий. Как и ранее, 4N + 4 условия составляют условия «Б». Недостающие условия естественно формулировать в угловых точках пластинки. Однако в каждой из этих точек можно записать по два условия, при этом условия для Ш будут выполнены автоматически. Поэтому вопрос о том, какие два из четырех оставшихся условий выбрать, имеет неоднозначное решение.
Чтобы избежать подобной ситуации, за систему разрешающих уравнений в этом случае примем систему «А», дополненную продифференцированными по п уравнениями (5с) и (7с). Порядок такой системы равен 4N + 8, и в качестве граничных условий для нее можно взять 4N + 4 условия «Б», к которым добавлены оставшиеся четыре условия в угловых точках. Такой подход делает постановку краевой задачи для указанной системы разрешающих уравнений математически корректной. Компонентами вектора У(п) в этом случае будут функции Ш} (п) (^ = 0, N), их производные до третьего порядка включительно, а также функции Ш-1 (п), Шм+1 (п) и их первые производные. Полученная таким способам краевая задача снова записывается в виде (9), (10).
Аналогичным образом строятся краевые задачи, если на сторонах £ = 0 и £ = 1 задаются граничные условия разного типа.
Численное решение краевых задач (9), (10) выполняется методом дискретной ортогонализации, который обеспечивает высокую точность результатов.
¿У — —
= АУ (п)+ ^ (п),
(9)
Н1У (0)= в1, Н2 У (1) = в2,
(10)
Рис. 2. Ортотропная пластинка
4. ОБ ИЗГИБЕ ПЛАСТИНКИ ИЗ ОРТОТРОПНОГО МАТЕРИАЛА
Рассмотрим частный случай прямоугольной пластинки, изготовленной из ортотропного материала. Предполагается, что в каждой точке одно главное направление упругости перпендикулярно срединной плоскости, а два других — направления x' и у' (рис. 2) — составляют со сторонами пластинки угол —. Пластинка изгибается поперечной нагрузки интенсивности q(x,y).
Если угол — равен нулю, то дифференциальное уравнение для прогиба точек срединной поверхности не содержит слагаемых с нечетными производными по переменным £ и п и 72 = 7з = 0 [1]. Применение ММСК для численного определения НДС такой пластинки детально изложено в работе [5].
При — = 0 для ортотропной пластинки остается в силе [1] уравнение (1). Однако теперь коэффициенты Yj (i = 1, 5) не являются независимыми, так как жесткости Dj (i, j = 1, 2,6) относительно осей x и у выражаются через главные жесткости Di, D2, Dk, D3 = Di v2 + 2Dk, вычисленные относительно осей x' и у' по известным формулам [1].
В этом случае НДС пластинки определяется в соответствии с методикой, изложенной в разд. 2 и 3.
5. ПРИМЕРЫ ЧИСЛОВЫХ РАСЧЕТОВ
Изложенная методика была применена при численном исследовании НДС квадратных пластинок (a = b = 1м) под действием распределенной статической нагрузки интенсивности q(x, у) = q0 = const. Предполагалось, что края £ = 0 и п = 0 жестко закреплены, а края £ = 1 и п = 1 свободны; угловая точка £ = п = 1 либо свободна — задача 1, либо подкреплена шарниром — задача 2. Вычисления выполнены для пластинок из ортотропного материала АГ-4с, у которого [6] E1 = 2.1 ■ 104 МПа, E2 = 1.6 ■ 104 МПа, G = 0.42 ■ 104 МПа, v2 = 0.07.
Значения функции W в долях величины q0/D* (D* = E1 ho/12(1 — v1v2)) в некоторых характерных точках пластинки при — = 0, п/12, п/6, п/3 приведены в табл. 1 (D* = 1.761 * 103h0). Для задачи 2 в столбце значений maxW в скобках указаны координаты точек, где достигается максимум W.
Таблица 1
Значения функции W (материал АГ-4с)
Задача Ф max W W(0.5; 0.5) W (0.5; 1.0) W(1.0; 0.5)
1 0.0 5.953E-2 1.081E-2 2.545E-2 2.670E-2
п/12 5.356E-2 1.073E-2 2.396E-2 2.938E-2
п/6 5.398E-2 1.126E-2 2.484E-2 2.864E-2
п/4 5.158E-2 1.101E-2 2.488E-2 2.437E-2
2 0.0 8.138E-3 (1.00; 0.59) 5.047E-3 6.197E-3 7.726E-3
п/12 8.191E-3 (1.0; 0.53) 5.231E-3 7.230E-3 8.155E-3
п/6 8.181E-3 (0.74; 0.67) 5.626E-3 7.427E-3 7.632E-3
п/4 7.879E-3 (0.70; 0.70) 5.628E-3 6.670E-3 6.632E-3
Аналогичные результаты для пластинок из дельта-древесины, которая обладает более сильной анизотропией по сравнению с АГ-4с, с постоянными Е1 = 3.05 ■ 104 МПа, Е2 = 0.467 ■ 104 МПа, С = 0.22 ■ 104 МПа, = 0.02 [1], представлены в табл.2 (Я0 = 2.626 ■ 103^).
Приведенные данные позволяют оценить влияние ориентации сторон пластинки по отношению к главным направлениям ортотропии и степень анизотропии материала на прогибы пластинки.
Таблица 2
Значения функции Ш (дельта-древесина)
Задача V max W W(0.5; 0.5) W (0.5; 1.0) W(1.0; 0.5)
1 0.0 1.0832E-1 2.165E-2 4.060E-2 5.789E-2
п/12 1.377E-1 2.261E-2 5.207E-2 7.374E-2
п/6 1.5151E-1 2.114E-2 5.5689-2 7.385E-2
п/4 2.735E-2 3.604E-2 1.492E-1 1.034E-1
2 0.0 3.129E-2 (1.00; 0.58) 1.275E-2 6.545E-3 3.010E-2
п/12 3.313E-2 (1.00; 0.55) 1.326E-2 1.016E-2 3.247E-2
п/6 2.506E-2 (1.00; 0.59) 1.190E-2 1.149E-2 2.387E-2
п/4 4.443E-2 (0.61; 1.00) 1.501E-2 4.032E-2 1.428E-2
Библиографический список
1. Лехницкий С.Г. Анизотропные пластинки. М.: ГИТТЛ, 1957. 463 с.
2. Недорезов П.Ф., Шевцова Ю.В., Ромакина О.М. Модифицированный метод сплайн-коллокации в задачах изгиба прямоугольных пластинок // Математическое моделирование и краевые задачи: Тр. Второй Всерос. науч. конф. Самара: СамГТУ, 2005. Ч. 1. С. 203-209.
3. Григоренко Я.М., Крюков Н.Н. Решение задач теории пластин и оболочек с применением сплайн-функций (Обзор) // Прикл. механика. 1995. Т. 31, № 6. С. 3-26.
4. Завьялов Ю.С., Квасов Ю.И., Мирошниченко В.Л. Методы сплайн-функций. М.: Наука, 1980. 352 с.
5. Шевцова Ю.В. Применение модифицированного метода сплайн-коллокации в задачах изгиба ортотропных прямоугольных пластинок // Проблемы прочности элементов конструкций под действием нагрузок и рабочих сред: Сб. науч. тр. Саратов: Изд-во Сарат. гос. техн. ун-та, 2005. С. 146-149.
6. Амбарцумян С.А. Теория анизотропных пластин. М.: Наука, 1967. 268 с.
