CC BY
22
ТРЕЩИНА СО СВЯЗЯМИ МЕЖДУ БЕРЕГАМИ НА ГРАНИЦЕ РАЗДЕЛА ДОРОЖНОГО ПОКРЫТИЯ И УПРУГОГО ОСНОВАНИЯ
Ш.Г. ГАСАНОВ, доктор технических наук, доцент Азербайджанский технический университет, г. Баку, E-mail: irakon63@hotmail. com
Исследуется плоская задача о трещине-расслоении, возникающей на границе раздела дорожного покрытия, сцепленного с упругим основанием из другого упругого материала, когда к поверхности покрытия приложена нормальная нагрузка. Задача о равновесии трещины на границе соединения материалов при действии внешних нормальных нагрузок и усилий в связях, препятствующих ее раскрытию, сводится к системе нелинейных сингулярных интегро-дифференциальных уравнений с ядром типа Коши. Условие предельного равновесия трещины с концевой зоной формулируется с учетом критерия предельной вытяжки связей.
КЛЮЧЕВЫЕ СЛОВА: дорожное покрытие, трещина, упругое основание, кусочно-квадратичная аппроксимация, механика разрушения.
Анализ [1, 2] состояния дорожных покрытий на упругом основании выявил, что материалам таких соединений свойственны трещиновидные несплошности. Между дорожным покрытием и упругим основанием могут появляться переходные зоны, в которых физико-механические свойства материала отличаются от свойств соединяемых материалов. Эти повреждения на границе раздела сред покрытия и упругого основания могут иметь как естественное происхождение (расслоения, включения, поры), так и вызываться технологическими процессами. Несмотря на важное значение отмеченных факторов на прочность многослойных структур до настоящего времени эти вопросы не нашли должного в методах расчета и проектирования дорожной одежды. Следовательно, разработка расчетных моделей исследования повреждений дорожного покрытия представляет актуальную проблему.
Примем следующие упрощающие предположения относительно работы дорожного покрытия:
1) дорожное покрытие является неразрезной полосой бесконечной длины неизменного поперечного сечения, лежащей на сплошном упругом основании;
2) вертикальные силы приложены в плоскости симметрии дорожного покрытия, а боковые и продольные силы не влияют на изгибающий момент и на напряженно-деформированное состояние, вызванное процессом контактирования колеса с дорожным покрытием.
Для расчета напряженно-деформированного состояния пары дорожное покрытие - упругое основание приходим к следующей задаче теории упругости.
Пусть в декартовых координатах х, y имеем двухслойное тело, состоящее из покрытия толщины h с упругими характеристиками G1 (модуль сдвига) и ¡и1 (коэффициент Пуассона), сцепленной с упругой полуплоскостью с характеристиками G2 и ц2 (рис. 1).
Рассмотрим задачу механики разрушения для двухслойного тела, когда к наружной поверхности приложена нормальная нагрузка p(x). На границе поверхности дорожного покрытия y = h , |x| <<х> касательное напряжение равно нулю. Считаем, что под действием этой нагрузки на линии раздела упругих сред y = 0 в силу недостаточной адгезионной прочности между покрытием и полуплоскостью (основанием) образуется трещина. Принимается, что при действии нормальной нагрузки на поверхности покрытия взаимодействие между
поверхностями в концевых зонах трещины характеризуется нормальными и касательными напряжениями сцепления.
Рис. 1. Расчетная схема задачи Выделим части трещины длиной d (концевые зоны), примыкающие к ее
вершинам {I - d < |х| < I, у = о), в которых берега трещины взаимодействуют,
так что это взаимодействие сдерживает раскрытие трещины. Взаимодействие берегов концевой зоны моделируется путем введения между берегами трещины связей, имеющих заданную диаграмму деформирования. Физическая природа таких связей и размеры концевой зоны зависят от вида соединяемых материалов. Считается, что закон деформирования связей задан. В общем случае он представляет нелинейный закон деформирования [3 - 5].
При действии внешней нормальной нагрузки в связях, соединяющих берега трещины на границе раздела сред, возникают нормальные qy(x) и касательные qxy(x) усилия. Величины этих напряжений заранее неизвестны и подлежат определению в процессе решения задачи.
Используем принцип суперпозиции. Рассматриваемое состояние двухслойного тела можно представить в виде суммы двух напряженных состояний:
1) адгезионное соединение материалов без трещины под действием нормальной нагрузки р(х) на наружной поверхности покрытия;
2) адгезионное соединение материалов с трещиной на границе раздела сред, при этом на берега трещины сносятся напряжения, имеющие место на этой линии, где расположена трещина, в первом напряженном состоянии.
Граничные условия задачи для первого напряженного состояния запишутся в виде (верхний индекс 1 соответствует покрытию, верхний индекс 2 - полуплоскости):
при у = h при у = 0
^ =- Р( х);
г(1) = 0
ху
и(1) + ш(1) = и(2) + ш(2);
^(1) + г-т(1) = ^(2) + iт(2) у ху у ху
(1) (2)
где i2 = -1, при у ^ -да перемещения и напряжения исчезают.
Для решения краевой задачи (1), (2) используем четыре функции Папкови-ча-Нейбера Fm {х, у) {п, т = 1, 2). Из них две для покрытия (верхний индекс 1) и две для полуплоскости (верхний индекс 2).
Напряжения и перемещения выражаются через функции Папковича-Ней-бера по известным формулам [6]
<о , 4CF(m) C2 F(m) C2 F(m) y = 2(1 - ¡Um - y d F
2Gm dy Cy Cy
2
T(m) C xy _
2G„ Cx
m
(3)
яр (») яр (т) (1 - 2цт )р2(т) - у ^
Су су
^ = _®Т-уС^; ^ .(3 - 4^,)^- У^ .
Сх Сх сУ сУ
В виду симметрии задачи по х, используем cos-преобразование Фурье. Примем, что
ад
Р1(1) = ^[Ashay + 5chay]cosaxda ,
о
ад
Р2(1) = j[Cshay + Dchay]аcosaxda , (4)
F1(2) = JEea cosaxda , F2(2) = JFeaya cosaxda .
0 0 Удовлетворяя функциями (3), (4) граничным условиям (1), (2), получим систему шести линейных алгебраических уравнений относительно шести неизвестных функций A(a), 5(a), C (a), D(a), E (a), F (a):
P
2(1 - u )(Cchah + Dshah)- Ashah - Bchah - ah(Cshah + Dchah) =---, (5)
2Ga
(1 - 2u )(Cshah + Dchah) - Achah - Bshah - ah(Cchah + Dshah) = 0, B = E , (3 - 4u )D - A = (3 - 4u2)F - E,
^[2(1 - u )C - B] = G2 [2(1 - ¡2 )F-E], ^[(1 - 2u )D - A] = G2 [(1 - 2u )F - E ],
2 1
где трансформанта преобразования Фурье p = — J p(x) cos axdx .
* 0
Решая алгебраическую систему уравнений (5) методом последовательного исключения неизвестных, находим коэффициенты A (a), B(a), C (a), D(a), E (a), F (a).
С помощью соотношений (3), (4) находим компоненты напряжений <rly и т' на IX < I, y = 0
(x,0) = -J -2ah ja(1 -ax)ea + ¿(1 + ax)e~a(4h-x) +
n 0 a + be + te [ +1 [r e-a(2h-x) + R2 e -a(2h+x) ]j cosahda, (6)
1> (x,0) = - 1 , ^\ -2ah ja(1 - 2u + ax)e-- + n 0 a + be + te h [
+1 [дзех) + R4е-а(2Ъ+х)]^sm ahda .
2
Здесь a = 3 - + 10к + Зк2 - 12^1к -12ц2к +к - 4^2к2, Ь = 3 - 4^1 - 6к + 4^2к + 4^1к + Зк2 - 4^2к2,
г = 10-24^ -4к +16^ + 8ка2/2 + 4h2а2 - 16ка2/2 - 12к 1а1/2 + + 16^2к2а2/2 + 8^2к2 - 6к2 + - 16^2к + 8^2к, R1 = а1 (а/г +1)+ с(а/г + 2) - г, R2 = а1 (а/ -1) + с(2-а/)-г, R3 = а1ах + с(1 + ах)-г, R4 = ах(с -а1) + г,
с = 2 + 4к - 6к2 - 8^2к + 8ц2к2
= С2
а1 = 4а/ + 8ка/- 12к2а/- 16^2ка/ + 16к2ц2а/ . Граничные условия на берегах трещины для второго напряженного состоя-
ния будут
^ - =- —у -
-у - = -(-1 -
(-1 - ) ,
(—у - ^ )+ Чу (х) - Чху (X),
х < I - d I - d < |х| < ^.
(7)
Условия (1) - (7) должны быть дополнены соотношением, связывающим раскрытие берегов трещины и усилия в связях. Это уравнение, без потери общности, в поставленной задаче можно представить в виде
и+ (х, 0) - и" (х, 0) -1 [и + (х, 0) - и- (х, 0)] = с(х, -)[ду (х) - 1дху (х)] , (8)
где — - модуль вектора усилий в связях; С(х,ст) - эффективная податливость
связей, зависящая от натяжения связей, (и+-и) - нормальная составляющая
раскрытия берегов в концевой области трещины, (и+ - и ) - касательная составляющие раскрытия берегов в концевой области трещины.
Эффективную податливость связей в соотношении (8) из соображений теории размерности можно представить так
С (х,-) = ^ (х,-),
Е,
(9)
где Е(х,ст) - безразмерная функция; Еэ - эффективный модуль упругости связей; Н - линейный размер, пропорциональный толщине зоны неоднородности на участке соединения материалов.
Из-за того, что длина трещины мала по сравнению с размерами двухслойного тела, то компоненты раскрытия берегов трещины можно представить в форме
и 1х
(х, 0)-и (х, 0) = — С0Е(х,—)чу(х), и+(х, 0)-и (х, 0) = — С0Е(х,—)д (х).
Е
Е
с = Н
С I,
где параметр С0 характеризует относительную податливость связей в концевой зоне трещины.
Приближенные представления для раскрытия берегов трещины при произвольной нагрузке на ее берегах могут быть получены на основе решения для раскрытия трещины при действии сосредоточенной силы, приложенной на поверхности трещины [7]. В случае двухслойного тела в дальнейшем будет удобнее использовать выражение для производных раскрытия трещины.
Выражение для производных перемещений берегов трещины при действии произвольных нагрузок на ее поверхности, приведенное в [8], после преобразований и симметрии искомых функций, примет вид
д
[и+(х,0)-и-(х,0)-i (и+(х, 0) - и (х, 0))] = А £+4 (^ + — )
1Р
+ i-„ -
2 А
+
4 - х
у1е2 - х21 4+х
\(х-, (*)-(#4
0 4 - х
+
(10)
е^ 4 -х2 Здесь
(х^у (У)-Уху (фу
Чу (х) - Чху (х) = (Ру (х) - iPxy (х))(^ ,
— (х) - ^х) )(ТГх^ ,
А = -4
—у - ^ =
1 ( к + 1 к. + 1
(11)
+
а
к1,2 = 3 - 4 ^1,2 , Р =
2 /
1пЯ
2ж :
4 =
а2 к1 + а1к 2 + а2
Отметим, что чу (х), Ру (х) - четные функции, а чху (х), Рху (х) - нечетные.
Функции ру (х) и рху (х) подлежат определению из решения задачи.
Для коэффициентов интенсивности напряжений от действия внешних нагрузок и связей в концевой зоне, найдем
К1 + КД =
с/%Р
м
,1/2+р
;«
я п>»+^
(12)
с/%Р
М
. ч1/2+р
12+РД12/[4+4 1 [чу у)+^(ук.
4
Для определения напряжений в связях в концевой зоне трещины получим интегро-дифференциальное уравнение. Перейдем к безразмерным переменным т = х/4 , г = 4/4 . Подставляя формулы (8) - (9) в (10) после некоторых преобразований, получаем систему нелинейных интегро-дифференциальных уравнений относительно неизвестных функций ру (т) и рху (т):
-[е(т, -)(чу (т) - iЧy (т))]- а( )(чу (т) + iЧy (т))+
С0 д Еэ дт1
2 А
+ -
(13)
ж
л/Т-т211
1 -т
+т
-Р 1
/К(т, г)тру (х) - ирху (x)]dt
1-^ 4
=41^41т)--
1 + 4
2 А
1 -т
У1 + т
Р1
/ К (т, г (
^у (?) - Итх
1(х)^.
Здесь ядро К (т, г) интегрального оператора в (13) имеет вид
К (т, г ) =
1
г2-т2
(14)
Решение системы уравнений позволяет определить неизвестные функции (т) и рх (т) и, тем самым, напряжения в связях чу (х), Чху (х).
е
е
2
Методика численного решения.
Полученные уравнения (13) представляют собой систему нелинейных сингулярных интегро-дифференциальных уравнений с ядрами типа Коши. Для их решения используем коллокационную схему кусочно-квадратичной аппроксимацией неизвестных функций. Перейдем к безразмерным переменным
Чу(0= Чу(*)/Ро; Чху(т) = Чху(х)/Ро; Ру(0 = Ру(х)/Ро; Р*у(т) = Р*у(х)/Ро> где
р0 - силовой фактор. Преобразуем систему уравнений (13) к форме удобной
для численного решения. Неизвестные функции Ру (т) и Ру (т) ищем в виде
Ру(т)=лТУ^Г, РХу(т)=лТУ^Г. (15)
Проведем дифференцирование по переменной т в левой части уравнения (13). После соответствующих преобразований, умножения обеих частей уравнения на (л/Т-7)/сов^т) и разделения действительных и мнимых частей, найдем
& (т)
А, {т,а)^Вг] {т,а)/, (т)+а {К,(т,г)/}.(,)Л, = Dl (т) , (16)
ат \-djt
(I, 1 = 1, 2),
А.(т,*) = (1 -тУ{т,а)[8г] + Яу^(т)] , Я, = (1 -8,)sgn(I-,),
В (т,ст)= г(т,°) -(1 -т)
¿у| М- Ыт)-1 V яу Г
1 1 + т 2) 11 2 1 + т
«*>"я)-+я ^
К,(т,,) = Нтй[т(2-,)+,(, -+ Я,^(т)), Н(т,;
л/1 + т , -т
Dl (т) представляют собой действительную и мнимую, соответственно, правые части уравнения (13).
Первый шаг итерационного решения уравнения (16) заключается в решении этого уравнения для линейно-упругих связей. На каждой последующей итерации уравнения (16) решается для квазиупругих связей вида (8) с эффективной податливостью вдоль концевой зоны трещины и зависящей от величины модуля вектора усилий в связях, найденного на предыдущем шаге решения. Вычисление эффективной податливости осуществляется подобно нахождению секущего модуля в методе переменных параметров упругости [9]. Процесс итерации заканчивается, как только усилия в связях, найденные на двух последовательных итерациях, мало отличаются друг от друга.
Для анализа предельного равновесия трещины-расслоения с концевой зоной необходимо рассмотреть два условия разрушения: 1) условие разрыва связей на краю концевой зоны и 2) условие роста вершины трещины.
Для условия продвижения трещины использовали силовое [10] условие
К, + гкц = к1с + гк11с. (17)
Из этого соотношения с учетом формулы (12) найдем соотношение, связующее внешнюю нормальную нагрузку и напряжения в связях.
В качестве первого условия разрушения использовали критерий критического раскрытия трещины. Считали, что разрыв связей на краю концевой зоны (х0 = I - Л) происходит при выполнении следующего условия
|и+(хо)-и (хо)-¡{и+(хо)-и (хо= 5С, (18)
где 8С - предельная вытяжка связи.
При использовании силового условия роста трещины, характеристики состояния предельного равновесия трещины-расслоения в двухслойном теле определяются из совместного решения системы уравнений (16), (17), (18).
На рис. 2 приведены графики распределения нормальных напряжений в связях чу/р0 (р0 - силовой фактор) вдоль концевой зоны трещины при Н/1 = 0,1 для различных относительных размеров концевой зоны а/1 .
9у/Р<> . .
4,0
3,0
2,0
0 0,25 0,50 0.75 ^
Рис. 2
Усилия в связях всегда максимальны на краю концевой зоны при постоянной вдоль концевой области податливости связей. Характер изменения касательных напряжений вдоль концевой зоны подобен изменению нормальных напряжений, но абсолютные значения касательных напряжений значительно ниже. На основании расчетов в таблице приведены значения абсолютных максимумов модуля вектора усилий параметров для различных значений параметров к и р при Н/1 = 0,1; Еэ = 1,5 ГПа.
к Р Р0
1 0 6,888
2 0,053351 6,814
5 0,110318 6,603
10 0,138312 6,407
50 0,166633 6,214
100 0,170675 6,172
Анализ численных результатов показывает, что с увеличением параметра к (1 < к < 1000) происходит перераспределение компонент усилий. Касательные усилия в связях значительно возрастают, а нормальные снижаются. При к ^ да (р ^ 0,175) нормальные и касательные компоненты напряжений сравнимы по величине.
На рис. 3 приведены зависимости отношения К11 / К1 от размера концевой зоны трещины а/7.
0,5 , 0,3\ \ \ 0,15
а/1 = 0,7 \
а/1 = 0,9
0,9
0,6
0,3
0 0,25 0,50 0.75 dj t
Рис. 3
1. Бабков В.Ф., Андреев О.В. Проектирование автомобильных дорог. - Изд. 2-е. -М.: Транспорт. 1987. - 368 с.
2. Щясянов Ш.Щ., Пириев Й.М., Гараисаейв Н.М. Автомобил йол-ларынын няглиййат-истисмар эюстяриъиляринин йцксялдилмяси. -Бакы: Азярбайъан няшриййаты, 2009. - 360 с.
3. Budianscky B., Evans A.G., Hutchinson J.W. Fiber-matrix de bonding effects on cracking in aligned fiber ceramic composites// Int. J. Solid Structures, 1995, v.32, No 3-4, p.315-328.
4. Ji H., de Gennes P.G. Adhesion via Connector Molecules: The Many-stitch Problem. Macromolecules, 1993, vol. 26. P. 520-525.
5. Goldstein R. V., Bakirov V.F., Perelmuter M.N. Modeling of the adhesion strength and fracture kinetics of the microelectronic package polymer-polymer joints // Int. Proc. Inst. Phys. Technol. Russian Ac. A Sci. v. 13. Modeling and Simulation of submicron Technology and Devices, 1997, p. 115-125.
6. Уфлянд Я.С. Интегральные преобразования в задачах теории упругости. - М.; Л.: Изд-во АН СССР, 1963. - 367 с.
7. Rice J.R., Sih G.C. Plane problems of cracks in dissimilar media // Trans. ASME, J. Appl. Mech., 1965, v. 32. P. 218-224.
8. Слепян Л.И. Механика трещин. - Л.: Судостроение, 1981, 295 с.
9. Биргер И.А. Расчет конструкций с учетом пластичности и ползучести // Изв. АН СССР. Механика, 1965, №2. С. 113 - 119.
10. Черепанов Г.П. Механика хрупкого разрушения. - М.: Наука, 1974. 640с.
A CRACK WITH INTERFACIAL BONDS BETWEEN FACE ON BORDER OF THE UNIT OF THE ROAD COVERING AND THE ELASTIC BASIS
Hasanov Sh.H.
The plane problem about the crack - stratification, arising on border of the unit of the road covering linked to the elastic basis from other elastic material when normal loading is enclosed to a surface of a covering is investigated. The problem on equilibrium of a crack on border of connection of materials at action of external normal loadings and tractions in the bonds interfering its disclosing, is reduced to system nonlinear singular the integro-differential equations with a nucleus such as Cauchy. The condition of limit equilibrium of a crack with a tip zones is formulated, taking account in view of criterion of limit traction of the bonds.
KEY WORDS: road covering, crack, elastic base, criterion of limit traction of the bonds, integro-differential equations.
