CC BY
36
Механика
УДК 539.3
ЗАДАЧА О ЦЕНТРАЛЬНОЙ ПРОДОЛЬНОЙ ТРЕЩИНЕ НОРМАЛЬНОГО ОТРЫВА С НАПОЛНИТЕЛЕМ В ПОЛОСЕ
Н.Н. Антоненко1
Предложен способ решения задачи о центральной продольной трещине нормального отрыва с наполнителем в полосе. Для решения задачи использовано интегральное преобразование Фурье. Задача сведена к интегро-дифференциальному уравнению относительно функции, связанной со скачком вертикальных перемещений на берегах трещины. Приведены результаты численных расчетов, которые иллюстрируют влияние наполнителя трещины, толщины и упругих характеристик полосы на коэффициенты интенсивности напряжений.
Ключевые слова: полоса; трещина; наполнитель; коэффициенты интенсивности напряжений; интегральное преобразование Фурье; интегро-дифференциальное уравнение.
Введение
Рассматривается задача о напряженно-деформированном состоянии полосы, содержащей прямолинейную трещину, параллельную границам полосы. В литературе рассматриваются два типа задач о продольных трещинах в полосе. К первому типу можно отнести задачи о трещине в полосе, границы которой свободны от напряжений, а к берегам трещины приложены нагрузки [1-7]. Ко второму - задачи о трещинах со свободными берегами с различными видами условий на границах полосы [1]. В вышеуказанных работах рассматриваются трещины типа разрезов, однако на практике трещины могут возникать в результате непроклейки слоистых элементов конструкций. Дефекты такого вида можно считать заполненными неким клейким веществом [8]. В данной статье рассматривается задача о центральной продольной трещине с наполнителем в полосе. В рамках этой модели предполагается, что скачки вертикальных перемещений на берегах трещины пропорциональны нормальным напряжениям на ее берегах [8].
Постановка задачи
Рассмотрим упругую полосу толщины 2Н. Материал полосы будем характеризовать двумя упругими характеристиками: модулем сдвига [I и коэффициентом Пуассона V . Нижняя граница полосы жестко заделана, а к верхней границе приложена нормальная сосредо-ченная нагрузка. Введем две локальные декартовые системы координат (г = 1,2) так, как показано на рис. 1.
На серединной плоскости полосы при г1 = Н (z2 = 0), имеется продольная трещина с наполнителем, занимающая область х е (—а; а). Напряжения и перемещения, которые
относятся к верхней половине полосы (х е (—¥ +»), 0 < z1 < Н ), будем обозначать нижним индексом 1, а к нижней (хе (—¥ +¥), 0 < z2 < Н ) - 2.
Граничные условия:
^ (х,0) = х) , (х,0) = 0, П2 (х,Н) = 0, ^ (х,Н) = 0. (1)
Условия на серединной плоскости полосы:
I
3
Рис. 1
1 Антоненко Нина Николаевна - кандидат физико-математических наук, старший преподаватель, кафедра общей математики, Запорожский национальный технический университет. E-mail: antonenkonina@i.ua
О-,2 (x,0) = szl (^h ) , -х22 ( X,0)=txz1 (^ h) ,
ГА(х), IX < а,
и2 (х,0) = и1 (х,h), я2 (х,0)-я1 (х,h) = < (2)
[0, > а,
—л (х,h) = е/ (х), |х| < а , (3)
где А(х) = \/а2 - х2/(х), /(х)е , /(±а)ф 0, е - коэффициент, который характеризует
наполнитель.
Требуется определить скачок вертикальных перемещений на берегах трещины и коэффициенты интенсивности напряжений (КИН).
_ _ хк _ _ а ~ О ~ h т е
Введем безразмерные величины хк = — , гк = —, а = —, О =—, h = —, т = —, е = —
и = и , 1~ = Я, -хвк = , = , где I, Е - характерные величины ([/] = м, [Е] = Па). В
дальнейшем тильды над функциями и переменными будем опускать, считая, что все преобразования выполняются над безразмерными величинами.
Метод решения
Для решения задачи воспользуемся интегральным преобразованием Фурье:
1 ¥ ^
/(Х)= | /(х)в'Ххёх, /(х) = | /(X)е~'ХхёХ.
_¥ 2- _¥
В пространстве трансформант компоненты напряженно-деформированного состояния полосы могут быть представлены в виде линейных комбинаций вспомогательных функций, которые связаны с напряжениями и перемещениями точек ее верхней границы такими формулами [9]:
а = —г (Х,0), Р = (Х,0), 7 = т(Х,0), 3 = -Х-хг (Х,0), (4)
р
где 5 = -Хи (X,,), Ж = рЯ(X,2), р =1X1.
Трансформанты напряжений и перемещений полосы имеют вид [9]:
2т Ж(X,г) = ((2-а)р,-ар,еЬр,)а + 2(-ар,р, + еЬр,)0 +
+2 ((1 -а) р,-юр, еЬ р, )у-ар, р1§, (5)
2т 5(X,z) = ар,р,а + 2((1 -а)р, + ар,еЬр,)Ь +
+2 (ар, р, + еЬ р, )у+(( 2-а) р, + ар, еЬ р, )3, (6)
(X,,) = (еЬрI-ар,рI)а+2а( р,-р,еЬр,)Ь-2ар,р?у-
-((1 -а) р, + аpz еЬ р, )3, (7)
IX
—--Х2 (X,,) = (-(1 -а)р, + ар,еЬр,)а + 2ар,р,р +
+2а($Ь р, + р, еЬ р, )у+( еЬ р, + ар, р, )3. (8)
Запишем условия (1)-(2), с учетом формул (4), в пространстве трансформант Фурье:
«1 = 0, 31 = 0, «2 = —!1 (X,h), ¿2 = -^^ (X,h),
р
Г2 =т^1 (X, h), ь-тжх (X, h )=трм (X), ^ (X, h ) = 0, ж2 (X, h ) = 0,
а
где М(X) = | А(I)е^сИ .
а
Антоненко Н.Н. Задача о центральной продольной трещине
нормального отрыва с наполнителем в полосе
Используя формулы (5)-(8), из последних соотношений найдем ß, g, a2, ß2, g2, ^ • Трансформанты нормальных напряжений в точках срединной плоскости полосы sz1 (X,h) примут такой вид:
s (X, h )= pAi(p) QA (p >, (9)
D (p >
где A1 (p)=-mw(a1 + a2 e~2pi + a3 e"4p + a4 e"6p -a1 e"8p),
A2 (p ) = b1e~p1 + b2e~3 p1 + b3 e"5 p1 + b4 e"7p, D (p ) = a1 + ke~4p1 + a1 e"8p , a1 = w2 - 2w, a2 = 2w(-w+2wp1 + 2wp2), a3 = 8p1 (w + w2 + 2w2p2), a4 = 2w(w+2wp1 - 2wp2),
b1 = -a1 (1 + p1), b2 = 2w + w2 + 3w2p1 + 4w2p2, b3 = 2w + w2 - 3w2p1 + 4w2pj2, b4 = -a1 (1 - p1),
22
wi = 0,5/((1 -n)), p1 = pÄ , W = 1 -w, k = -2(2 + a1 + 8w2p2).
Используя связь между оригиналами и трансформантами нормальных напряжений (9) в точках срединной плоскости полосы и последнее из условий (2), получим интегральное уравнение задачи:
cf (X )= I A (t )K (t - x) dt - QL (x),
pcf ч 7
—a
¥ pA1 (p) ¥ a2 (p)
где K(z) = I-cos pzdp, L(z) = I —cos pzdp .
1 ' J0 D (p) 1 'J D (p)
Выделим в ядре интегрального уравнения регулярную Ki (z) и сингулярную K2 (z) части:
aa
pcf (х )= I A (t )K1 (t — x) dt — j«J A (t) K2 (t — x) dt — gL (x), (10)
—a —a
¥ f A (p) 1 ¥ A1 (p)
где K1 (z) = I p I —j—f- + jw cos pzdp , K2 (z) = I p cos pzdp, lim —= — jw. J 1 D (p) I 0 p®¥ D (p)
о V ' /
Интеграл K2 (z) расходится, поэтому его будем понимать в следующем смысле:
+¥ +¥ 2 2 1 г г q — z 1 K2 (z )= I p cos pzdp = lim I p cos pze~qpdp = lim -- = —-.
0 q®0+0 0 q®0+0 (q2+z2 )2 z2
ö A( t)
Поскольку I --dt при x = t имеет неинтегрируемую особенность, то воспользуемся фор-
- a (t - x)
мальным равенством:
=—¡лмЩ=—А!^)а + гЖШ = . (1,)
1 (г — х)2 J 4 ; I г — х ) г — х -1 г — х -1 г — х
—а\ 1 —а 4 у —а —а —а
Математическое обоснование (11) приведено в [10]. В (11) использован тот факт, что берега трещины должны смыкаться, т.е. Л(±а) = 0 .
Окончательно интегральное уравнение задачи (10) принимает вид:
/
а __а (\/а2 — 12/(I))
ко/ (х )= N а2 — г2 / (*) Я (г — х) Ж + [®Р-— Ж — QL (х). (12)
^ г — х
—а —а
Неизвестную функцию /( х) будем искать в виде:
/(х) = Еа,и21 (х), (13)
г=0
где и21 (х) - полиномы Чебышева второго рода.
Для определения неизвестных « раскладываем левую и правую части уравнения (12) в линейные комбинации полиномов Чебышева и приравниваем коэффициенты при полиномах одинакового порядка. Количество членов в линейных комбинациях выбираем из условия, чтобы решения, полученные на п -м и п + 1-м шагах, отличались на некоторую наперед заданную величину.
[40^/ (г))
t - х
Используя асимптотические оценки, полученные в [8], для интеграла |
-а
получаем
та/(а )-у/2а
—z1 (х, h ) =-Ц=--+ О (1), г = х - а, при х ® а + 0.
2\1 г
Формула для вычисления коэффициентов интенсивности напряжений принимает вид:
К1 = та/(а)4яа , где функция / ( х ) вычислена по формуле (13).
Численные результаты
Численные расчеты проведены для трещины длины 2а = 2 . Полоса находится под действием нормальной сосредоточенной силы 0 = 1. Ниже в табл. 1-4 приведены результаты, которые иллюстрируют зависимость КИНов от упругих характеристик полосы, ее полуширины и наполнителя трещины.
Таблица 1
Зависимость коэффициентов интенсивности напряжения от толщины полосы (V = 0,3, т = 1)
//а 2 3 4 5 10 100
К1 е = 1 0,11205 0,08416 0,06611 0,05408 0,02788 0,00282
К' , е=0 0>1 -а 0,30996 0,21577 0.16476 0,13299 0,06733 0,00676
Таблица 2
Зависимость коэффициентов интенсивности напряжения от модуля сдвига полосы (V = 0,3, !<1 а = 10)
т 1 2 3 4 5 10
КI е = 1 0,02788 0,03945 0,04577 0,04976 0,05250 0.05890
К' , е=0 0>1 -а 0,06733 0,06733 0,06733 0,06733 0,06733 0,06733
Таблица 3
Зависимость коэффициентов интенсивности напряжения от коэффициента Пуассона полосы (т = 1, /'/а = 10)
V 0,1 0,2 0,3
К, . г— , е = 1 0^1 -а 0,02375 0,02563 0,02788
К< , е=0 0^1 -а 0,06697 0,06709 0,06733
Антоненко Н.Н.
Задача о центральной продольной трещине нормального отрыва с наполнителем в полосе
Таблица 4
Зависимость коэффициентов интенсивности напряжения от коэффициента, который характеризует наполнитель с (п = 0,3, т = 1, Ча = 10)
c 1 2 3 4 5 10
Kj q4PÎ 0,027SS 0,01757 0,012S2 0,01009 0,00S32 0,00443
Из табл. 1-4 можно сделать выводы:
1) увеличение полуширины полосы и коэффициента c приводит к уменьшению КИНов;
2) к увеличению КИНов приводят увеличение модуля сдвига и коэффициента Пуассона полосы;
3) для случая трещины, берега которой свободны от напряжений ( c = 0 ), упругие характеристики полосы практически не влияют на КИНы.
Литература
1. Murakami, Y. Stress intensity factors handbook: in 2 Vol. / Y. Murakami. - Pergamon Press, 19S7. - Vol. 1. - 1566 p.
2. Сметанин, Б.И. Некоторые задачи о щелях в упругом клине и слое / Б.И. Сметанин // Инж. ж. МТТ. - 196S. - № 2. - С. 115-122.
3. Саврук, М. П. Двумерные задачи упругости для тел с трещинами / М.П. Саврук. - Киев: Наук. думка, 19S1. - 324 с.
4. Александров, В.М. Продольная трещина в преднапряженном тонком упругом слое со свободными границами / В.М. Александров, Б.И. Сметанин // ПММ. - 2005. - Т. 69, Вып. 1. -С.150-159.
5. Fichter, W.B. Stresses at the tip of a longitudinal crack in a plate strip / W.B. Fichter. - Washington: National Aeronautics and Space Administration, 1967. - 55 p.
6. Александров, В.М. Продольная трещина в ортотропной упругой полосе со свободными гранями / В.М. Александров // Изв. РАН. МТТ. - 2006. - № 1. - С. 115-124.
7. Пожарский, Д.А. Асимптотические решения смешанных задач для упругой полосы и клина / Д.А. Пожарский, А.А. Молчанов. - Вестник Донского государственного технического университета. - 2010. - Т. 10. - С. 447-454.
S. Антоненко, Н.М. Моделювання трщини нормального вщриву з наповнювачем на межi багатошарового пакета та твплощини / Н.М. Антоненко, 1.Г. Величко // Прикл. проблеми мех. i мат. - 2011. - Вип. 9. - С. 141-149.
9. Ткаченко 1.Г. Двомiрна мшана задача термопружност для багатошаровоï основи / 1.Г. Ткаченко // Прикладш проблеми мехашки i математики. - 2005. - Вип. 3. - С. 70-7S.
10. Лифанов, И.К. Метод сингулярных интегральных уравнений и численный эксперимент (в математической физике, аэродинамике, теории упругости и дифракции волн) / И.К. Лифанов. -М.: Янус, 1995. - 520 с.
Поступила в редакцию 22 января 2015 г.
Bulletin of the South Ural State University Series "Mathematics. Mechanics. Physics" _2015, vol. 7, no. 3, pp. 65-70
THE PROBLEM OF THE CENTRAL LONGITUDINAL OPENING MODE CRACK WITH THE FILLER IN THE STRIP
N.N. Antonenko
The solution to the problem of the central longitudinal opening mode crack with the filler in the strip is given. Within this model we suppose that disturbances of vertical displacement at crack edges are proportional to the normal stress at its edges. Fourier integral transformation is used to solve the problem. The problem is reduced to integral differential equation as function which is connected with disturbances of vertical displacement at crack edges. On the basis of numerical results the conclusion is made that the increase of the half-width of the strip and coefficient which is characterized by the filler leads to the reduction stress intensity coefficients (SIC), the increase of modulus of rigidity and Pois-son's strain ratio will lead to the increase of stress intensity coefficients, elastic behaviour of the strip doesn't affect stress intensity coefficients for the cracks with stress-free edges.
Keywords: strip; crack; filler; stress intensity coefficients; Fourier integral transformation; integral differential equation.
References
1. Murakami Y. Stress intensity factors handbook: in 2 Vol. Pergamon Press, 1987. Vol. 1. 1566 p.
2. Smetanin B.I. Nekotorye zadachi o shchelyakh v uprugom kline i sloe (Some of the problem of cracks in an elastic wedge and layer). Inzh. zh. MTT. 1968. no. 2. pp. 115-122. (in Russ.).
3. Savruk M.P. Dvumernye zadachi uprugosti dlya tel s treshchinami (Two-dimensional elasticity problems for the bodies with cracks). Kiev, Nauk. dumka Publ., 1981. 324 p. (in Russ.).
4. Aleksandrov V.M., Smetanin B.I. Prikladnaya matematika i mekhanika (Journal of Applied Mathematics and Mechanics). 2005. Vol. 69. Issue 1. pp. 150-159. (in Russ.).
5. Fichter W.B. Stresses at the tip of a longitudinal crack in a plate strip. Washington, National Aeronautics and Space Administration, 1967. 55 p.
6. Aleksandrov V.M. Izv. RAN. MTT. 2006. no. 1. pp. 115-124. (in Russ.).
7. Pozharskiy D.A., Molchanov A.A. Vestnik Donskogo gosudarstvennogo tekhnicheskogo uni-versiteta. 2010. Vol. 10. pp. 447-454. (in Russ.).
8. Antonenko N.M., Velichko I.G. Prikladnye problemy matematiki i mekhaniki (Applied Problems of Mathematics and Mechanics). 2011. Issue 9. pp. 141-149. (in Ukrainian).
9. TKaneHKO I.r. Prikladnye problemy matematiki i mekhaniki (Applied Problems of Mathematics and Mechanics). 2005. Issue 3. pp. 70-78. (in Ukrainian).
10. Lifanov I.K. Metod singulyarnykh integral'nykh uravneniy i chislennyy eksperiment (v mate-maticheskoy fizike, aerodinamike, teorii uprugosti i difraktsii voln) (The method of singular integral equations and numerical experiment (in mathematical physics, physics of aerodynamics, elasticity and wave diffraction)). Moscow, Yanus Publ., 1995. 520 p. (in Russ.).
Received 22 January 2015
1 Antonenko Nina Nikolayevna is Cand. Sc. (Physics and Mathematics), Senior Lecturer, Department of General Mathematics, Zaporizhzhya National Technical University. E-mail: antonenkonina@i.ua
