CC BY
32
References
1. Uglov A. L., Erofeev V. I., Smirnov A. N. Akus-ticheskij kontrol' oborudovaniya pri izgotovlenii i ehkspluatacii [Acoustic control equipment in manufacturing and operation]. Moscow, Nauka, 2009. 279 p. (in Russian).
2. Glagolevskij B. A., Moskalenko I. B. Nizkochastot-nye akusticheskie metody kontrolya v mashinos-troenii [Low-frequency acoustic methods of inspection in mechanical engineering]. Leningrad, Mashinostroenie, 1977. 203 p. (in Russian).
3. Akhtyamov A. M. Teoriya identifikacii kraevyh uslovij i eyo prilozheniya [The theory of identification of boundary conditions and applications]. Moscow, Fizmatlit, 2009. 272 p. (in Russian).
4. Ahtyamov A. M., Muftahov A. V., Ahtyamo-va A. A. Ob opredelenii zakrepleniya i nagruzhen-nosti odnogo iz koncov sterzhnya po sobstvennym chastotam ego kolebanij [Determination of loading and fixation of one end of the rod based on the natural frequencies]. Vestn. Udmurtsk. un-ta. Matem. Mekh. Komp'yut. nauki, 2013, vol. 3, pp. 114-129 (in Russian).
5. Vatulyan A. O., Vasilev L. V. Ob opredelenii parametrov uprugogo zakrepleniya neodnorodnoj
balki [Determination of the parameters of an inho-mogeneous beam elastic fixation]. Ekologicheskij vestnik nauchnyh centrov CHEHS, 2015, no. 3, pp. 14-20 (in Russian).
6. Philippov A. P. Kolebaniya deformiruemyx system [Oscillations of deformable systems]. Moscow, Mashinostroenie, 1970. 736 p. (in Russian).
7. Kristensen R. Vvedenie v teoriyu vyazkoupru-gosti [Introduction to the of viscoelasticity theory]. Moscow, Mir, 1974. 340 p. (in Russian).
8. Bogachev I. V., Vatulyan A. O. Obratnye koe-hfficientnye zadachi dlya dissipativnyh operatorov i identifikaciya svojstv vyazkouprugih materialov [The inverse coefficient problem for dissipative operators and identification of the properties of vis-coelastic materials]. Vladikavkazskij matematich-eskij zhurnal, 2012, vol. 14, iss. 3, pp. 31-44 (in Russian).
9. Kalitkin N. N. Chislennye metody : ucheb. poso-bie [Numerical methods. Tutorial]. St. Petersburg, BHV-Petersburg, 2011. 586 p. (in Russian).
Please cite this article in press as:
Vatulyan A. O., Vasilev L. V. Determination of Attaching Parameters of Inhomogeneous Beams in the Presence of Damping. Izv. Saratov Univ. (N.S.), Ser. Math. Mech. Inform., 2016, vol. 16, iss. 4, pp. 449-456 (in Russian). DOI: 10.18500/1816-9791-2016-16-4-449-456.
УДК 539.3
ИЗГИБ МНОГОСВЯЗНЫХ АНИЗОТРОПНЫХ ПЛИТ С КРИВОЛИНЕЙНЫМИ ОТВЕРСТИЯМИ
С. А. Калоеров1, А. И. Занько2
1 Калоеров Стефан Алексеевич, доктор физико-математических наук, профессор кафедры теории упругости и вычислительной математики, Донецкий национальный университет, [email protected]
2Занько Алена Игоревна, аспирант кафедры теории упругости и вычислительной математики, Донецкий национальный университет, [email protected]
Предложен приближенный метод определения напряженного состояния тонких плит с криволинейными отверстиями, заключающийся в использовании комплексных потенциалов теории изгиба анизотропных плит, аппроксимации контуров отверстий дугами эллипсов и берегами прямолинейных разрезов, конформных отображений, в представлении комплексных потенциалов рядами Лорана и определении неизвестных коэффициентов рядов обобщенным методом наименьших квадратов. Изотропные плиты рассматриваются как частный случай анизотропных плит. Численные исследования проведены для плиты с квадратными и треугольными отверстиями. Исследованиями установлена высокая степень точности получаемых результатов. Отмечены значительные отличия известных в литературе результатов от полученных данным подходом.
Ключевые слова: анизотропная пластинка, изотропная пластинка, криволинейное отверстие, комплексные потенциалы, обобщенный метод наименьших квадратов.
DOI: 10.18500/1816-9791 -2016-16-4-456-464
ВВЕДЕНИЕ
Для плит с круговыми или эллиптическими контурами в свое время были разработаны методы решения задач определения напряженного состояния, основывающиеся на представлении исходных
комплексных потенциалов рядами Лорана и по полиномам Фабера и определении их коэффициентов методом рядов [1-3]. Эти методы с некоторой достоверностью давали значения искомых величин, правда, для достаточно неблизких расстояний между контурами. Для случаев же плит с криволинейными отверстиями имеются лишь отдельные исследования, результаты которых далеки от истинных.
Нашими исследованиями установлено, что в случае плоской задачи теории упругости для изотропной пластинки с отверстием в форме многоугольника из существующих методов наиболее точные результаты дает метод, использующий функции Колосова - Мусхелишвили с построением конформных отображений
.«>=4+Е $ ,
"(0 = Я(С +> I, (1)
где Я, ск — постоянные, определяемые из геометрических характетик контура отверстия. При этом для достижения достаточной точности в отображающей функции нужно оставлять довольно большое количество слагаемых, а комплексные потенциалы определять из граничных условий методом интегралов типа Коши [4]. Но для анизотропной пластинки с криволинейным отверстием точки при соответствующих аффинных преобразованиях и конформных отображениях переходят в различные точки, и метод интегралов типа Коши не применим. Эта тенденция сохраняется и в случае многосвязных пластин. Существующие для указанных случаев методы решения задач, связанные с применением метода малого параметра для преобразования уравнений контуров и использованием метода рядов при удовлетворении граничным условиям, приводят к значительным погрешностям получаемых результатов и на сегодняшний день мало эффективны при решении современных задач инженерной практики. В наших работах [5,6] предложено аппроксимировать криволинейные контуры дугами эллипсов и берегами прямолинейных разрезов, а для удовлетворения граничным условиям использовать дискретный метод наименьших квадратов [5] или обобщенный метод наименьших квадратов (ОМНК) [6]. Наиболее простым при численной реализации оказывается ОМНК. В данной статье последний подход использован для решения задач об изгибе анизотропных плит с произвольными криволинейными контурами. При этом изотропные плиты рассматривались как частный случай анизотропных плит. Численные исследования проведены для плит с квадратным и треугольным отверстиями. Исследовано влияние упругих характеристик материалов плит и геометрических характеристик отверстий на значения изгибающих моментов.
1. ПОСТАНОВКА И МЕТОД РЕШЕНИЯ ЗАДАЧИ
Рассмотрим занимающую многосвязную область, тонкую анизотропную плиту с эллиптическими отверстиями (I = 1,£) с полуосями , центрами в точках (х01 , у01 > и углами ^ между положительным направлением оси Ох и полуосями щ (рис. 1).
Плита находится под действием распределенных поперечных усилий р1 и изгибающих моментов на контурах и моментов М^0, Му°, Ну, действующих на бесконечности. В общем случае эллипсы могут переходить в прямолинейные разрезы, располагаться произвольно относительно друг друга, в том числе касаться, пересекаться, образуя контуры сложной конфигурации, аппроксимируемые дугами эллипсов и берегами прямолинейных разрезов. Рис. 1. Многосвязная плита
Определение напряженно-деформированного состояния рассматриваемой плиты при использовании комплексных потенциалов теории изгиба [1,3,7] сводится к нахождению комплексных потенциалов (гк> (к = 1, 2) из соответствующих граничных условий. Эти функции определены в многосвязных областях 5к, получаемых из заданной области 5 аффинными преобразованиями:
гк = х + Дк у, (2)
где дк — корни характеристического уравнения
Б22Д4 + 4^26Д3 + 2 (£>12 + 2£ее> Д2 + 4^16Д + Би = 0, (3)
Бу = ЬуБ0 — жесткости материала плиты;
Ьц = (а22абб - а2б) /А, &12 = (а1ба2б - а^абб> /А,
bie = («i2«26 - «i6«22) /А, ^22 = («ii«66 - «i6) /А, &26 = («12«16 - «26«ii) /А, bee = («il«22 - «У /А,
А =
«11 «12 «16 «12 «22 «26 «16 «26 «66
(4)
а^ — коэффициенты деформаций материала, = 2к3/3, к — полутолщина плиты.
Многосвязные области 5к ограничены контурами Ьк1, соответствующими контурам Х области 5 при указанных аффинных преобразованиях. Используя методы конформных отображений и разложений функций в ряды Лорана для производных этих функций, получаем выражения вида [5,7,8]
L о
Wk (zk) = ГкZk + (AkiZk + Bki)ln (zk - Zki) + n
i=1 —1 Z
в которых rk — постоянные, определяемые из системы уравнений
«kin
_ _ Z n
l=1 n— 1 Zkl
(5)
k=1
2 Re ^ ^ = C11 MX° + C21 My00 + C31 HXy, 2 Re ^ Mkrk = C12M£° + C22MO + C32HXO,
k=1
2 2 1 2 Re J] ^krk = C13MO + C23My00 + C33H0У, 2 Re ^ — rk = 0, (6)
k—1
k=1
Mk
C11 = (2D22 D66 - /А1, C21 = (2D16D26 - 2D12D66) /А1,
C31 = (2D12D26 - 2D12D22) /А1, C12 = (D12D26 - D16D22) /А1, C22 = (2D12D16 - 2D11D26) /А1, C32 = (DuD22 - D22) /А1, C13 = (2D16D26 - 2D12D66) /А1, C23 = (2D12D66 - 2D262 /А1,
D11 2D12 2D16
А1 = D12 2D26 2D22
D16 2D66 D26
Akl, Bkl — постоянные, определяемые из решения систем уравнений
2
k—1
2 Re ^ 1, Mk, Mk, — iAkl = 0, 0, 0,
Mk
P
2nD
11
2Re ^ ( 1, Mk ,Mk, -— ) iBkl k—1 Mk
= 0, 0, -
Mxl 2nD
M
yl
22
2nD
11
(7)
(8)
Р и Мх1, Ыу1 — главный вектор и компоненты главного момента относительно начала координат усилий и моментов, приложенных к контуру Х, (к1 — переменные, определяемые из конформных отображений
ШкЛ
Zk = Zkl + Rkl (kl +
(kl J
(9)
внешности единичного круга |£ki1 ^ 1 на внешности эллипсов Lki,
Zkl = Xoi + Дк yol,
Rki = [«i (cos + Дк sin ) + ibi (sin - Дк cos /2,
mki = [«i (cos ^i + Дк sin ^i) - ibi (sin ^i - Дк cos /2Rki,
(10)
ак1п — неизвестные постоянные, определяемые из граничных условий. В случае прямолинейного разреза в последних соотношениях нужно принять Ь = 0.
Конформные отображения можно построить и по координатам концов осей эллипсов Х (хАг ,уАг), (хв1 , Ув), (жСг ,УСг), (хВг ). В этом случае в формулах (10) нужно принять
а = 1V(хА - )2 + (УАг - Ус )2, &г = ^\/(хвг - )2 + (Увг - Увг )2,
L
2
2
xoi =
XA, + XC,
2
Vol =
Уа, + VC,
2
pi = arctg
Уа, - VC, xA, — xC,
(11)
При решении задач возникает необходимость аппроксимации части криволинейного контура в виде сторон угла в окрестности его вершины Аг (хАг ,уА1) дугой окружности Ц, вписанной в этот угол с центром на биссектрисе угла (рис. 2).
Обозначим углы между положительным направлением оси Ох и сторонами угла Аг через ^ и . Тогда угол Аг будет равен — ^. Поэтому угол наклона указанной биссектрисы к оси Ох = (^ + /2. Выбрав на сторонах угла отрезки 8г от вершины угла для радиуса вписанной окружности и координат ее центра, будем иметь
a = Sitg ((pi+1 - pi) /2), |O'All = y!6f + a2,
x0i = xA, - |O'Ai1 sin((Pi + Pi + 1) /2) ,
Voi = Уа, - |O'Ai1 cos ((pi + pi + 1) /2) .
Рис. 2. Аппроксимация части криволинейного контура дугой окружности
Зная координаты x'0i, y0i центра O[ вписанной окружности Ц, его радиус a[ и угол между осью Ox и биссектрисой угла Ai и принимая, что a[ = Ц, по формулам (9) и (10) найдем функцию, отображающую внешность единичного круга на внешность окружности Ц, вписанной в угол Ai.
Неизвестные постоянные akin, входящие в функции (5), будем определять из граничных условий на контурах Li [1,3,7]
2 Re Y^ 9kii Wk (zk) = fii
i = 1, 2,
(12)
k=1
в которых
gki1 = Pk / №k, gki2 = qk,
fi1 = - fi dx - mi dy + cidx + оц, fi2 = - mi dx - fi dy + cidy + C2i,
рк = Бц +2^16дк + Б12ц\, ак = Б12 + 2Б2б^к + Б22^1, ¡1 (5) = —у Рг (з)
о
сг — вещественные, с^ — комплексные постоянные.
Граничным условиям (12) будем удовлетворять обобщенным методом наименьших квадратов [9]. Для многосвязных областей этим условиям удобнее удовлетворять в дифференциальной форме. Продифференцировав условия (12) по дуге контура, выбрав на контурах плиты набор точек ит (х1т,у1т) (ш = 1,Ы1) и удовлетворив в них этим условиям, получим
2 L те
2 Re S S S gkiiSk Pkin (tkim)
dfii (t im)
k=1 =1 n=1
L
-2 Re SJ gkiiTk + ^
k=1
=1
Aki In (zk - zki) +
ds
(Aki zk + Bki)
zk - zk
l = 1, L, m = 1, Mi, i = 1, 2,
(13)
S = dzk ' =__n_
где Sk = ds ■ - С1 Лы K2 - m,,)
tki x m + №k Vi 7
Систему (13) дополним уравнениями однозначности прогиба [7] для каждого отверстия
2
2 Re iaki 1Rki = 0, l = 1,L.
k=1
(14)
s
s
s
s
s
2
После нахождения псевдорешений дополненной системы уравнений (13) с использованием сингулярных разложений [6] постоянные akln, а следовательно, и функции Wk (zk), будут известными и по ним можно вычислять изгибающие моменты и перерезывающие силы в любой точке плиты [1,3,7]
2 2 (Mx ,My ,HXy) = -2 Re 'rk) W'k (zk), (Nx,Ny) = -2 Re ^(^k sk, -Sk) W''' (zk). (15)
k=1 k=1
Здесь rk = Die + 2DeeMk + ^26^1, Sk = -Die - (D12 + 2Dee) Mk - 3D26Mk - D22Mk• 2. РЕШЕНИЕ ЧАСТНЫХ ЗАДАЧ
Приведенный подход был использован для решения задач о напряженном состоянии плит с отверстиями различной конфигурации. Далее описаны некоторые из полученных результатов для плит с квадратными и треугольными отверстиями. При проведении расчетов количество членов в рядах (5) и точек на контурах Ml, для которых составлялись уравнения (13), увеличивалось до тех пор, пока граничные условия на контурах не удовлетворялись с достаточно высокой степенью точности (пока значения моментов на площадках, касательных к контурам, не удовлетворяли условию Mn < 10-3). Для такого удовлетворения граничным условиям, как показали численные исследования, необходимо было в указанных рядах оставлять от 40 до 50 членов, на каждом из контуров брать от 200 до 500 «коллокационных точек». Далее описаны некоторые из полученных результатов для плиты из КАСТ-В изотропный (материал М1) [10], стеклопластика косоугольной намотки (материал М2) [10], сосны (материал М3) [1]. Коэффициенты деформаций aj для этих материалов приведены в табл. 1.
Таблица 1
Постоянные материалов
Материал ац • 10-4, МПа-1 а22 • 10-4, МПа-1 ai2 • 10-4, МПа-1 aee • 10 4, МПа-1 an/a22
М1 72.100 72.100 -8.600 161.500 1.000
М2 10.000 2.800 -0.770 27.000 3.571
М3 2.381 0.100 -0.024 1.333 23.810
При этом для того, чтобы использовать приведенное выше решение для изотропной плиты коэффициент деформации материала 012, равный -8,600 ■ 10-4, в расчетах заменялся на -8,605 ■ 10-4, т. е. на значение, незначительно отличающееся от табличного. В противном случае при решении задачи возникнет деление на ноль, так как корни характеристического уравнения (3) для изотропного материала будут двукратными и равными % и —г.
Пусть плита ослаблена квадратным отверстием со сторонами длиной 2а, параллельными осям координат (рис. 3). Контур отверстия свободен от усилий. На бесконечности изгибающие моменты соответственно равны Му° = ту, М= ИХу = 0. Стороны квадрата будем представлять внешними берегами разрезов Ь\, Х4 вдоль сторон, с кон-
цами в точках А, В, С, Б. В этом случае в приведенном решении нужно принять
С
D
L = 4, A (a, a), B (-a, a), C (-a, -a), D (a, -a),
Рис. 3. Квадратное отверстие в плите
ai — a2 — a3 — a4 — a.
Расчеты проводились без аппроксимации сторон вблизи вершины углов и с их аппроксимацией. При проведении расчетов вместо равных нулю полуосей Ь считалось, что Ь/а; = 10-4. Это уменьшает вычислительные погрешности. В рассматриваемом случае если обозначать через 6 отрезки на сторонах от вершин углов до точек их закругления, будем иметь
а; = 6, 1в'1А11 = //26, х01 = а — 6, у01 = а — 6, х02 = —а + 6, у02 = а — 6, х0з = —а + 6, у03 = —а + 6, х04 = а — 6, у04 = —а + 6.
В табл. 2 приведены значения изгибающих моментов М3 по контуру отверстия на площадках, перпендикулярных ему (Му — для вертикальных сторон, Мх — для горизонтальных сторон) в точках с указанными координатами с точностью до интенсивности приложенных моментов ту как множителя.
Таблица 2
Значения моментов по контуру квадратного отверстия в плите из различных материалов
Г Координаты точек Без закругления в углах С закруглением при 6 Без закругления в углах С закруглением при 6
0.01 0.005 0.01 0.005
Материал М1 Материал М2
0 1.0001, 0.0000 1.3173 1.3159 1.3169 1.5115 1.5052 1.5159
15 1.0001, 0.2588 1.3301 1.3421 1.3394 1.5250 1.5598 1.5429
30 1.0001, 0.5736 1.4058 1.4117 1.4141 1.6238 1.6423 1.6459
44.9956 0.9998, 1.0000 13.4327 - - 16.8064 - -
45 0.9971, 0.9971 - 3.5483 - - 2.2968 —
45 0.9985, 0.9985 - - 4.3247 - - 2.7517
45.0044 1.0000, 0.9998 -2.6345 - - -2.5382 - -
60 0.5736, 1.0001 0.2633 0.2535 0.2531 0.2357 0.2286 0.2316
75 0.2588, 1.0001 0.2687 0.2639 0.2598 0.2452 0.2335 0.2331
90 0.0000, 1.0001 0.2681 0.2622 0.2673 0.2479 0.2488 0.2524
Окончание табл. 2
Без С закруглением
Г Координаты закругления при 6
точек в углах 0.01 0.005
Материал М3
0 1.0001, 0.0000 1.6034 1.6002 1.6085
15 1.0001, 0.2588 1.6379 1.6522 1.6475
30 1.0001, 0.5736 1.7621 1.7475 1.7498
44.9956 0.9998, 1.0000 21.3145 - -
45 0.9971, 0.9971 -- 2.7712 -
45 0.9985, 0.9985 - - 3.0246
45.0044 1.0000, 0.9998 -0.7640 - -
60 0.5736, 1.0001 0.0733 0.0709 0.0712
75 0.2588, 1.0001 0.0745 0.0702 0.0709
90 0.0000, 1.0001 0.0739 0.0759 0.0743
Рассматривались случаи квадрата с прямыми углами (без закруглений), а также квадрат с закругленными углами вблизи вершин с малыми значениями радиусов (равных в данном случае отрезкам от вершин углов) закруглений 6. В табл. 2 координаты точек и значения 6 отнесены к полудлине а стороны квадрата.
На рис. 4 изображены графики распределения этих моментов для точек отрезков ЕА и на сторонах. Для остальных частей сторон квадрата значения моментов легко восстанавливаются по имеющейся симметрии. Кривые 1, 2, 3 соответствуют плите из материалов М1, М2, М3.
Из данных табл. 2 и рис. 4 видно, что при приближении к вершинам квадрата значения моментов стремятся к бесконечности, вдоль горизонтальных сторон к —го, вдоль вертикальных сторон — к +го, причем рост моментов вдоль вертикальных сторон больше, чем убывание вдоль горизонтальных сторон. Аппроксимация сторон углов вблизи вершин снимает особенность моментов в этих точках, причем чем меньше радиус закругления, тем больше значения моментов в вершинах (закругленных). Чем выше «степень анизотропии» (отличие ац/а22 от 1), тем больше значения моментов вдоль вертикальной стороны и меньше вдоль горизонтальной стороны.
ту 2.5 2.0 1.5 1.0 0.5 0.0
-0.5 Е
У
— ——- 4 I-1
, 3 2 1
А
Е
Рис. 4. Графики распределения изгибающих моментов в плите с одним квадратным отверстием
Б-,
С-
У
Заметим, что с целью сравнения результатов с известными из литературы для плиты из березовой фанеры с квадратным отверстием при всестороннем изгибе [11] исследования проводились и для этого случая. Значения изгибающих моментов, полученные в работе [11], оказались значительно отличными от наших вблизи вершин. Они достаточно близки к полученным нами лишь в точках, близких к серединам сторон квадрата.
Пусть бесконечная плита имеет два квадратных отверстия со сторонами длины 2а, параллельными осям координат (рис. 5). Стороны квадратов будем представлять внешними берегами разрезов Ь\, Ь2, Ьз, Ь5, Ьг, Ь8 вдоль сторон, т.е. с концами Л^ А\, В\, С\, £>1 и А2, В2, С2, -£>2- Расстояние между отверстия-
ми обозначим через с. В этом случае Ь = 8, (а, а), В\ (—а, а), х С\ (—а, —а), £>1 (а, —а), А2 (За + с, а), В2 (а + с, а), С2 (а + с, —а),
В2 (За + с, —а), «1 = а2 = аз = а4 = а$ = а,б = а? = «8 = а.
Расчеты проводились без аппроксимации сторон вблизи вер-рис. 5. Два квадратных отверс- шин углов и с их аппроксимацией. При проведении расчетов вме-тия в плите сто равных нулю полуосей Ь; считалось, что Ь^/а^ = 10-4.
В рассматриваемом случае если обозначать через 6 отрезки на сторонах от вершин углов до точек их закругления, будем иметь
а; = 6, |0'А; | = //26,
А с.
'01
= а - 6 У01 = а - 6
'02
'03
'05
= —а + 6, = 3а + с — 6,
У03 = —а + 6 У05 = а — 6
= —а + 6, = а — 6,
04
У02 = а — 6
У04 = —а + 6
'06
= —а + с + 6, у06 = а — 6,
'07
= а + с + 6, у07 = —а + 6, х08 = 3а + с — 6, у08 = —а + 6.
Численные исследования для плиты с двумя квадратными отверстиями были проведены для различных значений с/а, отношения расстояния между отверстиями с к полудлине их сторон а. Как показали расчеты, с уменьшением расстояния между отверстиями (отношения с/а) значения изгибающих моментов вблизи контуров и в точках перемычки растут.
Если расстояние между отверстиями более трех длин сторон (с/а > 6), то влияние одного отверстия на значения моментов около другого незначительно и им можно пренебречь. Для близких расстояний между контурами особенно большие значения моментов возникают в точках перемычки
вблизи контуров отверстий.
На рис. 6 изображены графики изменения изгибающих моментов Му/ту вблизи точки перемычки Е1 (а, 0) на контуре левого отверстия в зависимости от с/а. Кривые 1, 2, 3 соответствуют плите из материалов М1, М2, М3. Видно, что с уменьшением с/а значения моментов резко растут, особенно в плите из материала с сильной «степенью анизотропии» М3.
Пусть теперь плита с треугольным отверстием в виде равнобедренного треугольника с высотой 2а по оси абсцисс и углом при вершине ^ (рис. 7) изгибается моментами Му° = ту. В этом случае вершины треугольника располагаются в точках А (2а, 0), В (0, 2аtg (р/2)), С (0, —2а tg (р/2)).
Стороны треугольника будем представлять внешними берегами разрезов Ьь Ь2, Ь3 вдоль сторон, т. е. с концами А, В, С. В этом случае Ь = 3, а1 = а3 = 2а/ соб (^/2), а2 = 4а tg (^/2).
Численные исследования были проведены для различных значений угла ^ при вершине треугольника. Как показали расчеты, с уменьшением угла ^ значения изгибающих моментов вблизи боковых сторон треугольника уменьшаются, вблизи
Иу
3.0
2.5
2.0
1.5
\ 3
\ V2
0.0
0.4
0.8
1.2
1.6 с/а
Рис. 6. Графики изменения моментов Му в точке перемычки Е\ плиты с двумя квадратными отверстиями
У В
о
с
Рис. 7. Треугольное отверстие в плите
Му
15
10
основания растут. Если угол ^ равен нулю, то моменты в окрестности точки А получаются такими же, как в плите с одним эллиптическим отверстием-трещиной. Это видно из данных рис. 8, где приведены графики изменения Му/ту в точке оси Ох вблизи вершины А зависимости от угла
Кривые 1, 2, 3 соответствуют плите из материалов М1, М2, М3. Видно, что с уменьшением угла ^ значения моментов растут и при ^ = 0 они получаются такими же, как в плите с одним эллиптическим отверстием-трещиной, причем это значение не зависит от анизотропии материала. Последнее согласуется с известным фактом для пластинки с трещиной в плоской задаче, когда ау вблизи вершины на продолжении трещины (и коэффициент
интенсивности напряжений ^), чему при изгибе плиты соответствует Му, не зависит от анизотропии материала [5]. Как показали расчеты, в точках плиты на продолжении высоты, удаленных от вершины А на расстояние более высоты 2а, значения моментов такие же, как в плите без отверстия, т.е. влиянием отверстия на напряженное состояние можно пренебречь.
V Л 2
0
те/6
те/3
те/2
2те/3 Ф, рад
Рис. 8. График изменения изгибающих моментов вблизи вершины А треугольника
Библиографический список
1. Лехницкий С. Г. Анизотропные пластинки. М. : Гостехиздат, 1957. 463 с.
2. Космодамианский А. С. Напряженное состояние анизотропных сред с отверстиями или полостями. Киев ; Донецк : Вища шк., 1976. 200 с.
3. Меглинский В. В. Некоторые задачи изгиба тонких многосвязных анизотропных плит // Некоторые задачи теории упругости о концентрации напряжений и деформации упругих тел : сб. науч. тр. Саратов : Изд-во Сарат. ун-та, 1967. Вып. 3. С. 97-127.
4. Савин Г. Н Распределение напряжений около отверстий. Киев : Наук. думка, 1968. 887 с.
5. Калоеров С. А., Горянская Е. С. Двумерное напряженно-деформированное состояние многосвязного анизотропного тела // Концентрация напряжений. Механика композитов : в 12 т. Киев : А. С. К., 1998. Т. 7. С. 10-26.
6. Калоеров С. А., Добряк Д. А. Термоупругое состояние кусочно-однородной анизотропной пластин-
9.
10.
11.
ки // В1сн. Донец. ун-ту. Сер. А, Природ. науки. 2006. Вип. 2. С. 77-88.
Калоеров С. А. Комплексные потенциалы теории изгиба многосвязных анизотропных плит // Тео-рет. и прикладная механика. 2012. № 4 (50). С. 115-136.
Калоеров С. А., Авдюшина Е. В., Мироненко А. Б. Концентрация напряжений в многосвязных изотропных пластинках. Донецк : Изд-во ДонНУ, 2013. 440 с.
Форсайт Дж., Малькольм М., Моулер К. Машинные методы математических вычислений. М. : Мир, 1969. 280 с.
Космодамианский А. С., Калоеров С. А. Температурные напряжения в многосвязных пластинках. Киев ; Донецк : Вища шк., 1983. 160 с. Космодамианский А. С., Митраков В. А. Изгиб конечной анизотропной пластинки с криволинейным отверстием // Прикладная механика. 1976. Т. 12, № 12. С. 96-99.
7
Образец для цитирования:
Калоеров С. А., Занько А. И. Изгиб многосвязных анизотропных плит с криволинейными отверстиями // Изв. Сарат. ун-та. Нов. сер. Сер. Математика. Механика. Информатика. 2016. Т. 16, вып. 4. С. 456-464. 001: 10.18500/1816-9791-2016-16-4-456-464.
Bending of Multiconnected Anisotropic Plates with the Curvilinear Holes
S. A. Kaloerov1, A. I. Zanko2
1 Stephen A. Kaloerov, Donetsk National University, 24, Universitetskaya str., 83055, Donetsk, Ukraine, [email protected] 2Alena I. Zanko, Donetsk National University, 24, Universitetskaya str., 83055, Donetsk, Ukraine, [email protected]
An approximate method for determination of the stress state of thin plates with curvilinear holes, consisting in the use of the complex potential theory of bending of anisotropic plates, approximating the contours of holes by ellipse arcs and straight sections, the use of conformal mapping, presentation of complex potentials by Laurent series and determining the unknown series coefficients of the generalized least squares method. Isotropic plates are considered as a special case of anisotropic plates. Numerical studies carried out for plates with square or triangular holes. Studies have established a high degree of accuracy of the results. Significant differences were noted in the literature known from the real results obtained by this approach.
Key words: anisotropic plate, isotropic plate, curved hole complex potentials, generalized least squares method.
References
1. Lekhnitskii S. G. Anisotropic Plates. New York, Springer-Verlag, 1984, 546 p.
2. Kosmodamianskii A. S. Napriazhennoe sostoian-ie anizotropnykh sred s otverstiiami ili polostiami [Stressed state of anisotropic media with holes or cavities]. Donetsk, Vishcha shk., 1976, 200 p. (in Russian).
3. Meglinskii V. V. Nekotorye zadachi izgiba tonkikh mnogosviaznykh anizotropnykh plit [Some problems of bending of thin multi-connected anisotrop-ic plates]. Nekotorye zadachi teorii uprugosti o kontsentratsii napriazhenii i deformatsii uprugikh tel : sb. nauch. tr. [Some problems of elasticity theory about stress concentration and deformation of elastic bodies. Collection of Papers], Saratov, Saratov Univ. Press, 1967, iss. 3, pp. 97-127 (in Russian).
4. Savin G. N. Raspredelenie napriazhenii okolo otverstii [Stress distribution near the holes]. Kiev, Nauk. dumka, 1968, 887 p. (in Russian).
5. Kaloerov S. A., Gorianskaia E. S. Dvumer-noe napriazhenno-deformirovannoe sostoianie mnogosviaznogo anizotropnogo tela [The two-dimensional stress-strain state of a multi-connected anisotropic body]. Kontsentratsiia napriazhenii. Mekhanika kompozitov : v 12 t. [The mechanics of composites. Stress concentration]. Kiev, A. S. K., 1998, vol. 7, pp. 10-26 (in Russian).
6. Kaloerov S. A., Dobriak D. A. Termouprugoe sos-toianie kusochno-odnorodnoi anizotropnoi plastin-
ki [Thermoelastic state of piecewise-homogeneous anisotropic plate]. Visn. Donets. un-tu, Ser. A, Prirod. nauki [Bulletin of Donetsk University, Series A, Natural Sciences], 2006, iss. 2, pp. 77-88 (in Russian).
7. Kaloerov S. A. Kompleksnye potentsialy teorii izgiba mnogosviaznykh anizotropnykh plit [Complex potentials the theory of bending of multi-connected anisotropic plates].Teoret. i prikladnaia mekhani-ka [Theoretical and Applied Mechanics], 2012, no. 4 (50), pp. 115-136 (in Russian).
8. Kaloerov S. A., Avdiushina E. V., Mironenko A. B. Kontsentratsiia napriazhenii v mnogosviaznykh izotropnykh plastinkakh [Stress concentration in multi-connected anisotropic plates]. Donetsk, Don-NU Press, 2013, 440 p. (in Russian).
9. Forsaythe G. E., Malcolm M. A., Moler C. B. Computer methods for mathematical computations. New Jersey, Prentice-Hall, 1977, 280 p.
10. Kosmodamianskii A. S., Kaloerov S. A. Temper-aturnye napriazheniia v mnogosvi-aznykh plas-tinkakh [Thermal stresses in multi-connected plates]. Kiev, Donetsk, Vishcha shk., 1983, 160 p. (in Russian).
11. Kosmodamianskii A. S., Mitrakov V. A. Izgib konechnoi anizotropnoi plastinki s krivolineinym otverstiem [Bend of finite anisotropic plate with a curvilinear hole]. Prikladnaia mekhanika [Applied Mechanics], 1976, vol. 12, no. 12, pp. 96-99 (in Russian).
Please cite this article in press as:
Kaloerov S. A., Zanko A. I. Bending of Multiconnected Anisotropic Plates with the Curvilinear Holes. Izv. Saratov Univ. (N.S.), Ser. Math. Mech. Inform., 2016, vol. 16, iss. 4, pp. 456-464 (in Russian). DOI: 10.18500/1816-97912016-16-4-456-464.
