CC BY
52
УДК 539.3
МЕТОД СПЛАЙН-КОЛЛОКАЦИИ И ЕГО МОДИФИКАЦИЯ В ЗАДАЧАХ СТАТИЧЕСКОГО ИЗГИБА ТОНКОЙ ОРТОТРОПНОЙ ПРЯМОУГОЛЬНОЙ ПЛАСТИНКИ
О.М. Ромакина, Ю.В. Шевцова*
Саратовский государственный университет, кафедра компьютерной алгебры и теории чисел, *кафедра геометрии
E-mail: [email protected], [email protected]
Приводится методика численного определения напряженно-деформированного состояния (НДС) изгибаемой тонкой орто-тропной пластинки при нетрадиционных способах закрепления ее краев. Результаты числовых расчетов для трех различных материалов позволяют оценить влияние анизотропии материала и способа закрепления пластинки на ее НДС.
Ключевые слова: метод сплайн-коллокации, пластинка, изгиб пластинки.
Spline-Collocation Method and its Modification in the Problems of Static Bending of Thin Orthotropic Rectangular Plate
O.M. Romakina, Yu.V. Shevtsova*
Saratov State University,
Chair of Computer's Algebra and Theory of Numbers, * Chair of Geometry
E-mail: [email protected], [email protected]
A numerical method for determining the stress-strain state (SSS) of a bended thin rectangular plate with non-classical boundary conditions is presented. Numerical results for three different materials can be used to estimate the influence of the material anisotropy and boundary conditions on its SSS.
Keywords: spline-collocation method, plate, bending of plate.
1. ПОСТАНОВКА ЗАДАЧИ. ОСНОВНЫЕ УРАВНЕНИЯ И СООТНОШЕНИЯ
Задача статического изгиба тонкой орторопной пластинки в системе координат Охуг, начало которой помещено в срединной плоскости, а координатные оси параллельны главным направлениям анизотропии, сводится к нахождению прогиба т(х,у), удовлетворяющего уравнению [1]
» д4 т +2Е д4 т + ид4 т ( ) (1)
дх4 дх2 ду2 ду4
Изгибающие моменты Мх, Му, крутящий момент Нху и перерезывающие силы Их, Иу выражаются через функцию т(х,у) по формулам [1]
д2 т д2 т \ _ / д2 т д2 т \ „ _ д2 т
Mx = -DAlT4 + , My = -D2\ + , Hxy = -2Dk
О 9 1 2 о 9 M irJ-y 2 I " 1 о О 1 О О I 1 "iy k о о ?
дх2 ду2 J \ дх2 ду2 J dxdy
N J d д3 w D д3 w \ N ( D д3 w D д3 w
x1
œ3 f) œf) о i2 i У v r) œ2 Я о t
дх3 дхду2 / у \ дх2ду ду
* = + (^3 + 2»^) , * = - (сдз + + - (2)
Здесь »г = ЕгЬ3 /12(1 — и2) (г = 1, 2) — жесткости изгиба; = ОН3/12 — жесткость кручения; »3 = »\и2+2»к; Ег, V (г = 1, 2) — модули Юнга и коэффициенты Пуассона для главных направлений анизотропии.
Граничные условия для функции т(х, у) записываются в соответствии со способом закрепления контура пластинки.
В дальнейшем рассматриваются прямоугольные пластинки с размерами в плане а х Ь, у которых края параллельны главным направлениям ортотропии.
В безразмерных переменных £ = х/а, п = у/Ь для безразмерного прогиба Ш(£,п) = т(х,у)/Н из уравнения (1) получаем
д4 Ш од4 Ш лд4 Ш
» о? дёэп?+»с4-опт = «ы (3)
© О.М. Ромакина, Ю.В. Шевцова, 2010
О.М. Ромакина, Ю.В. Шевцова. Метод сплайн-коллокации и его модификация
а формулы (2) принимают вид
* 2
Mx = -D* а
д2 W
+ V2 С2
,д2 W дп2
My = -D*а2 ( v1
д2 W д2 W
дС2
+ с2
1П* д3 W 2 д3 W Nx = -а ( D*^^ + D*c2
дС3
д3 W
NX = -MD*^ + (D* + 2D* )c2
дСдп2 д3 W дС2 дп
д п2
д3 W
Hxy = 2D* а2 с
д2 W д Сд п
N" = -a^D3 + D*c2
д3 W д п3
* * * д3 W * 2 д3 W
N* = -ас (D* + 2D* + D* с2
дС2дп
д п3
(4)
В формулах (3) и (4) обозначено Б* = Б^/а4 (г = 1, 2,3, к), с = а/Ь.
2. РЕШЕНИЕ МЕТОДОМ СПЛАЙН-КОЛЛОКАЦИИ
Для прямоугольной ортотропной пластинки точные аналитические решения краевых задач для уравнения (3) известны только в случаях опертого контура пластинки и когда две противоположные стороны подкреплены шарнирами при произвольных условиях на остальной части контура [1]. При других, более сложных способах закрепления применяются различные приближенные методы. В частности, для пластинки, у которой стороны £ = 0 и £ = 1 как-то закреплены, численное решение можно получить эффективным методом сплайн-коллокации [2]. Согласно этому методу решение уравнения (3) ищется в виде
N
w (С,п) = £ < (C)Wj (п), j=0
(5)
где через р^- (£) (^ = 0, N) обозначены линейные комбинации нормализованных В-сплайнов пятой степени В5,г(£) (г = —2,Ж + 2) [3], построенных на соответствующем расширении равномерной сетки, узлы которой определены по правилу £ = г/Ж (г = 0, Ж). Эти комбинации подбираются так, чтобы автоматически были выполнены условия при £ = 0 и £ = 1, и имеют вид: в случае свободно опертых сторон
ро(£) = 12В5,-2(£) — 3В5,- 1 (£) + В5,о(£), р 1 (£) = —В5,- 1 (£) + Ва, 1 (£),
Р2(£) = —Ва,-2(£) + Ва,2(£), р^-2(£) = Ва,^-2 (£) — Ва,*+2 (£), р^-1 (£) = — Ва^-1(£) + Ва,^+1(£), р^ (£) = Ва,^ (£) — 3Ва,*+1(£) + 12Ва,*+2(£);
при жесткой заделке
<о(С) = 41, 25B5,-2(С) - 0.12505,-1 (С) + £5,0 (С), <1 (С) = 32, 5£,,-2(С) - 2.25Д;,-1(С) + £5,1 (С),
<2(С) = 2, 25В5,-2(С) - 0.12505,-1 (С) + £5,2(С),
<N — 2(С) = B5,N-2(С) - 0.125B5,N+1(С) + 2.25£,,N+2(С),
<n-1(С) = B5,N-1(С) - 2.25b5,N+1(С) +32.5b5,N+2(С),
<N(С) = B5,N(С) - 4.125b5,N+1 (С) + 41.25b5,N+2(С);
при любом способе закрепления <(С) = B5,j(С) (j = 3, N - 3).
Система N +1 обыкновенных дифференциальных уравнений четвертого порядка для неизвестных функций Wj (п) (j = 0, N) получается из требования, чтобы результат подстановки выражения (5) в уравнение (3) выполнялся в точках коллокации С = С* (i = 0, N). Точки коллокации выбираются по правилу: С* = 6 + t/N (i = 0,N/2 - 1), С*/2 = С^2, С* = 6 - t/N (i = N/2 + 1,N), 0 < t < 1. Эта система стандартным приемом преобразуется в записанную в нормальной форме Коши систему уравнений первого порядка, которая в векторной форме может быть представлена в виде
dY — —
— = AY (п) + F (п), (6)
где У = Шг,
а Ш а2 Ш Ш
с!г! ' а г/2 ' а-^3
, (г = 0, N) — вектор неизвестных, а матрица А = (аг,8} и вектор
^ = {/г(п)}, г, 5 = 0, N имеют известные компоненты, которые выражаются через значения функций (£) и их производных и функции д(£, п) в точках коллокации.
Граничные условия для вектора У получаются из граничных условий на сторонах п = 0 и п = 1 и всегда могут быть записаны в виде
ЯхУ (0)= 61, Н2 У (1) = 62,
(7)
с прямоугольными матрицами Н1 и Н2 размерности (2N + 2) х (4N + 4) и (4N + 4) — мерными векторами е1 и е2, компоненты которых известны.
Краевая задача (6), (7) решается численно методом дискретной ортогонализации, который, как показывают решения многочисленных тестовых задач, обеспечивает получение практически точных результатов.
Указанный подход был применен при изучении напряженно-деформированного состояния (НДС) двух пластинок, у одной из которых контур жестко защемлен (задача 1), а у другой стороны £ = 0 и п = 0 закреплены жестко, а при £ = 1 и п = 1 выполняются условия шарнирного опирания (задача 2). При постоянной нагрузке интенсивности д(х, у) = до = 1кН/м2 вычисления проводились для квадратных пластинок а = Ь = 1м, Н = 0.01м из ортотропных материалов АГ-4с (Е1 =2.1 ■ 104МПа, Е2 = 1.6 ■ 104МПа, С = 0.42 ■ 104МПа, ^ = 0.07 [4]) и дельта-древесины (Е1 = 3.05 ■ 104МПа, Е2 = 0.467 ■ 104МПа, С = 0.22 ■ 104МПа, ^ = 0.02[1]) и для стальной пластинки (Е = 2 ■ 105МПа, V = 0.03 ). Графики изменения величины прогиба Ш вдоль диагонали £ = п представлены для задачи 1 на рис. 1, а, для задачи 2 — на рис. 1, б. На этих рисунках кривые 1 соответствуют материалу АГ-4с, кривые 2 — дельта-древесине, кривые 3 — стали.
-- 1 - 2 — 3
//
Ч\
//
//
ч
0 0.1 0.2 0.3 0.4 0.5 0.Б 0.7 0.8 0.9 1
а
Рис. 1. Изменение величины прогиба Ш вдоль диагонали: а - задача 1, б - задача 2
0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 б
Максимальные значения прогиба Ш и изгибающих моментов Мх и Му, а также значения моментов в серединах сторон контура пластинки приведены для задачи 1 в табл. 1, а для задачи 2 — в табл. 2. В табл. 2 в скобках для каждого значения указаны координаты точки, в которой эти значения получаются; значения в первой строке соответствуют материалу АГ-4с, во второй — дельта-древесине, в третьей — стали.
Таблица 1
Максимальные значения прогиба и изгибающих моментов для ортотропной пластинки.
Задача 1. Метод КМСК
тах Ш = Ш(0.5, 0.5) тах Мх = Мх(0.0, 0.5) тах Му = Му (0.5, 0.0) Мх(0.5, 0.5) Му (0.5, 0.5)
9.165Е-1 5.927Е-4 3.691Е-4 2.541Е-4 1.446Е-4
9.482Е-1 8.207Е-4 3.423Е-5 4.089Е-4 9.915Е-6
6.905Е-2 5.131Е-4 5.131Е-4 2.289Е-4 2.289Е-4
О.М. Ромакина, Ю.В. Шевцова. Метод сппайн-коппокации и его модификация
Таблица 2
Максимальные значения прогиба и изгибающих моментов для ортотропной пластинки.
Задача 2. Метод КМСК
max W max Mx max My Mx(0.5, 0.5) My (0.5, 0.5)
1.703E+0 (0.575;0.575) 8.307E-4 (0.0:0.575) 5.197E-4 (0.575;0.0) 3.484E-4 2.015E-4
1.838E+0 (0.575;0.550) 1.176E-3 (0.0;0.55) 3.220E-4 (0.575;0.0) 4.930E-5 1.552E-5
1.201E-1 (0.575;0.575) 6.940E-4 (0.0;0.575) 6.940E-4 (0.575;0.0) 3.042E-4 3.042E-4
3. МОДИФИЦИРОВАННЫЙ МЕТОД СПЛАЙН-КОЛЛОКАЦИИ
Если у пластинки закреплены одна или две смежных стороны, то метод сплайн-коллокации в классическом виде [2] неприменим, так как не удается подобрать такие комбинации B-сплайнов, чтобы тождественно удовлетворялись граничные условия на двух противоположных сторонах. В этом случае численное решение можно получить модифицированным методом сплайн-коллокации [4]. Детальное изложение этого метода применительно к ортотропным пластинкам приведено в [5].
Для определенности будем считать, что закреплена сторона £ = 0. Тогда согласно [4, 5] решение краевой задачи для уравнения (3) ищется в виде
N+2
W(£,n) = £ ъ (£)W (п). j=0
(8)
Здесь функции ^(£) = 0, N + 2), как и в классическом варианте метода сплайн-коллокации, должны обеспечить выполнение граничных условий при £ = 0, а функции ЭД^+1 (п) и ЭД^+2(п) определяются из алгебраических при задании на стороне £ = 1 прогиба и угла поворота или дифференциальных уравнений, если при £ = 1 заданы изгибающий момент и перерезывающая сила.
В обоих указанных случаях по-прежнему требуется выполнение в точках коллокации уравнения (3) и граничных условий на сторонах п = 0 и п = 1, а полученная таким образом система обыкновенных уравнений и граничные условия преобразуются к виду (6), (7).
Отметим, что при задании на стороне £ = 1 нагрузки к уравнениям четвертого порядка, полученным после подстановки (8) в (3), добавляются два уравнения второго порядка, которые получаются из граничных условий при £ = 1. Выполнение граничных условий на сторонах п = 0 и п = 1 в этом случае требуется не только в точках коллокации, но и в концевых точках стороны £ = 1. Тогда вектор неизвестных имеет компоненты
Y (п) = Wi,
dWi d2 Wi d3 Wi dn ' dn2 ' dn3
, Wn+1, WN+2,
d Wn+i dn
d Wn+2 \
dn )
(i = 0, N),
соответственно увеличивается размерность матриц A, H1, H2 и векторов F, e1, e2.
Совершенно аналогично решается задача, если сторона £ = 0 не закреплена, а например, свободна или загружена усилиями заданной интенсивности. В этом случае безразмерный прогиб W(£, п) ищется в виде
N+2
w(£,n)= £ B5,j(£)Wj(n), j = -2
а функции W-i (n) и W-2(n) определяются из дифференциальных уравнений второго порядка, которые получаются из граничных условий при £ = 0. Модифицированный метод сплайн-коллокации был применен для вычисления прогибов и изгибающих моментов в пластинках из стеклопластика АГ-4с и дельта-древесины и стали от действия распределенных усилий q(x, y) = q0 = const. Предполагалось, что угловые точки подкреплены шарнирами, а остальная часть граничного контура пластинок свободна от нагрузки и закрепления. Графики изменения прогиба вдоль диагонали £ = n для этого случая показаны на рис. 2, нумерация кривых принята такой же, как и на рис. 1, в силу симметрии графики построены для значений £ < 0.5. Значения прогиба W в долях q0/D* в характерных точках приведены в табл. 3, а изгибающих моментов Mx и Mx — в табл. 4.
О 0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 0.45 0.5
Рис. 2. Задача 3. Изменение величины прогиба Ш вдоль диагонали
Таблица 3
Значения прогиба ортотропной пластинки в характерных точках. Метод ММСК
max W W(0.5, 0.5) W(0.0, 0.5) W(0.0, 0.5)
1.901E+1(0.5;0.5) 1.901E+1 1.479E+1 1.247E+1
4.119E+1(0.0;0.5) 4.043E+1 4.119E+1 1.283E+1
1.393E+0(0.5;0.5) 1.393E+0 9.697E-1 9.697E-1
Таблица 4
Значения изгибающих моментов ортотропной пластинки в характерных точках.
Метод ММСК
max Mx = Mx (0.5; 0.0) max My = My (0.0; 0.5) Mx(0.5,0.5) My (0.5, 0.5)
1.947E-3 1.759E-3 8.573E-4 9.610E-4
2.909E-3 1.430E-3 4.753E-4 1.149E-3
1.504E-3 1.504E-3 1.117E-3 1.117E-3
Приведенные графики и таблицы позволяют оценить влияние анизотропии материала и способа закрепления пластинки на ее НДС.
Библиографический список
1. Лехницкий С.Г. Анизотропные пластинки. М.: ГИТТЛ, 1957. 463 с.
2. Григоренко Я.М., Крюков Н.Н. Решение задач теории пластин и оболочек с применением сплайн-функций (Обзор) // Прикл. механика. 1995. Т. 31, № 6. С. 3-26.
3. Завьялов Ю.С., Квасов Ю.И., Мирошниченко В.Л. Методы сплайн-функций. М.: Наука, 1980. 352 с.
4. Недорезов П.Ф., Шевцова Ю.В., Ромакина О.М. Модифицированный метод сплайн-коллокации в зада-
чах изгиба прямоугольных пластинок // Математическое моделирование и краевые задачи: Тр. Второй Все-рос. науч. конф. Самара: Изд-во Самар. гос. техн. ун-та, 2005. Ч.1. С. 203-209.
5. Шевцова Ю.В. Применение модифицированного метода сплайн-коллокации в задачах изгиба ортотропных прямоугольных пластинок // Проблемы прочности элементов конструкций под действием нагрузок и рабочих сред: Сб. науч. тр. Саратов: Изд-во Сарат. гос. техн. ун-та, 2005. С. 146-149.
