Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

книги / Методы оценки трещиностойкости конструкционных материалов

..pdf
Скачиваний:
14
Добавлен:
12.11.2023
Размер:
12.05 Mб
Скачать

Рис. 1. Схематическое изображение сечений тела с трещиной плоскостью, перпендикулярной к плоскости трещины (а) и совпадающей с плоскостью трещины (б).

размер области тела в окрестности вершины трещины, где материал деформирован за предел упругости, т. е. характерный линейный размер области предразрушения. Для некоторых классов матери­ алов эта область представляет собой вытянутую вдоль плоскости расположения трещины клинообразную область (см. рис. 1, а).

Разрушение, которое реализуется распространением трещины при выполнении условий автомодельности ее зоны предразрушения (характерный линейный размер зоны предразрушения мал по сравнению с размерами трещины и тела, х0 а), называется квазихрупким [145]. Если размеры зоны предразрушения х0 (пласти­ ческой зоны) соизмеримы с размерами трещины и тела, а пласти­ ческие деформации пе заполняют всего кинематического сечения разрушения, то разрушение в таком теле будем называть квазивязким.

В рамках концепций теории трещин [9, 82, 118, 145] разруше­ ние твердого тела рассматривается как процесс развития трещины. Такая теория более точно (по сравнению с классическими [138]) описывает явления разрушения. Она предусматривает следующий механизм разрушения. В процессе деформации тела в окрестности контура трещины возникает зона предразрушения. Развитие и фор­ мирование этой зоны происходит до тех пор, пока напряжения или деформации в ней не достигнут определенного (предельного) значения, т. е. пока не наступит предельно-равновесное состояние. В этот момент и происходит распространение (старт) трещины в зону предразрушения.

Поскольку предельно-равновесное состояние соответствует стартовому (начальному) моменту распространения трещины, оно характеризует только локальное разрушение тела. Математиче­ ские условия, которые описывают наступление предельно-равно­ весного состояния, называются критериями локального разру­ шения.

Формулировка критериев локального разрушения зависит от модельного представления зоны предразрушения. Остановимся подробнее на некоторых моделях локального разрушения твердых тел. С этой целью рассмотрим трехмерное тело, ослабленное плос­ кой трещиной с контуром L (рис. 1, б) и введем следующие обозна­ чения: а — характерный линейный размер трещины; х0 — харак­ терный линейный размер области предразрушения по нормали п к контуру трещины; Oraz — цилиндрическая система координат, выбранная так, что плоскость z = 0 совпадает с плоскостью трещи­ ны (случай сечения такого тела плоскостью, проходящей через ось Oz, показан на рис. 1, а); Н0 (а) — радиус-вектор контура трещи­ ны; R (а) — радиус-вектор линии пересечения поверхности зоны предразрушения с плоскостью z = 0 (см. рис. 1, б)\ Р — параметр внешней нагрузки, которая приложена симметрично относительно плоскости трещины. Имея в виду изложенное, рассмотрим неко­ торые основные модели механики хрупкого разрушения.

1. Модель Гриффитса — Ирвина. Эта модель [193] представля­ ет собой упругий континиум с трещиной (трещинами), для которого

максимальный линейный размер

области предразрушения х0

мал по сравнению с характерным

линейным размером трещины а

(или каким-то другим линейным размером тела), т. е. выполняются

условия автомодельности зоны предразрушения: х0

a; R0 (а) «

« R (а). В рамках названной модели считают, что

зона предраз­

рушения локализируется в малой окрестности у вершины трещины. Компоненты сгрр, стгр, агг тензора напряжений в локальной системе координат Ог0$ (см. рис. 1, а) в окрестности рассматривае­ мой макротрещины нормального разрыва можно записать так

[82,

145]:

 

а у = у ш ; /«(Р) + °(1),

где

Кх = К1 (Р1 а) — коэффициент интенсивности напряжений,

который зависит от конфигурации тела, размеров__^^фи£урации трещины, величины действующей нагрузки Р, но не зависит от координат г и Р; fa (|5) — известные функции угла Р; 0 (1) — огра­ ниченная величина при г0 ->■ 0.

Условием локального разрушения тела (страгивания трещины) является равенство (критерий Ирвина)

K 1 ( P „ a ) = K Ul

(1.1)

где Р* — предельное значение внешней

нагрузки для тела

с трещиной, т. е. такое минимальное значение нагрузки Р, при достижении которого наступает распространение трещины; К\с — значение трещиностойкости материала (максимальное значение коэффициента в уловиях плоской деформации.

Напряженное состояние в теле считается симметричным отно­ сительно плоскости трещины и описывается методами линейной теории упругости. При этом можно показать [82, 145], что коэффи­

циент интенсивности на­ пряжений К 1определяется следующей формулой:

К г = lim [У 2nr аг (г, а, 0)],

(1 . 2)

где г — малое расстояние по нормали п между неко­ торой точкой тела М и контуром трещины L (см. рис. 1, б)\ oz (г, а, 0) —

нормальные растягивающие напряжения oz(г, ос, z),

действующие в плоскости

*о&

рис 2. Схематическое представление тупиковой части трещины,

Сопротивление материала распространению трещины (локаль­ ное свойство материала) оценивается в этом случае только од­ ной характеристикой К\с (или Кс в случае обобщенного плоского напряженного состояния). Эта характеристика определяется непо­ средственно из эксперимента (см. гл. VI).

2. бк-Модель. Эта модель представляет собой

[82]

линейно­

упругое тело со следующими

дополнительными

свойствами:

1) максимальные растягивающие

напряжения в теле не

превосхо­

дят величины о0 — предела хрупкой прочности материала; 2) связь между напряжениями и деформациями определяется законом Гука, если упругие растягивающие напряжения не достигают величины а0; 3) в теле возникает зона предразрушения, если максимальные растягивающие напряжения, вычисленные на основе линейной теории упругости, превосходят величину сг0; 4) области предразру­ шения рассматриваются как трещины-разрезы (рис. 2) размера х0 = R (а) R0 (а), противоположные берега которых притяги­ ваются друг к другу с напряжением а0, если расстояние между берегами такой трещины не превосходит некоторой величины

6К, и не взаимодействуют между собой,

если расстояние больше,

чем бк. Величины а0 и бк связаны равенством

ст06к = 2у,

(1.3)

где у — плотность энергии разрушения материала, т. е. энергия, расходуемая на образование единицы новой поверхности в данном материале.

Из условий непрерывности напряжений в деформируемом теле с трещинами в рамках 6к-модели следует [82], что

(1.4)

Здесь vn (г) = ип [г, о, Р, R0 (a), R (а), а0] — нормальная состав­ ляющая вектора смещения точек берегов трещины-разреза, вы­ численная методами линейной теории упругости в рамках сформу­

лированной модели. Равенство (1.4) определяет значение функции R (ос) для области предразрушения.

Сформулированную модель характеризуют, помимо упругих коэффициентов (закона Гука), еще две величины: о,0 и бк. Рассмат­ риваемая модель описывает развитие области предразрушения в окрестности контура трещины (увеличение размера х0), а усло­ вия распространения трещины (локального разрушения) сводятся к следующему равенству:

2v„ [О, О, Р*, R0(a), R (а), ст0] = 6К.

(1.5)

Если характерный линейный размер х0 области предразруше­ ния мал по сравнению с размером трещины, т. е. х0 ^ а, как это показано в работе [82], то из условия (1.5) и равенства (1.3) получа­ ем условие (Г.1).

Условие (1.5) является базовым для разработки локального критерия хрупкого разрушения в виде критического раскрытия трещины (KPT-критерий). Этот критерий сводится к следующему [82]: трещина начнет увеличиваться в длину (распространяться), как только ее раскрытие в вершине достигнет критического (для данного материала при заданных условиях) значения 6К.

3. Аналог бк-модели для упруго-пластического тела. Рассмот­ рим трехмерное тело, ослабленное плоской макротрещиной. Под макротрещиной здесь подразумевается такая трещина, для которой при заданных внешних нагружениях, размерах тела и его механи­ ческих характеристиках выполняются условия автомодельности зоны предразрушения в условиях предельного равновесия, т. е. имеет место понятие о тонкой структуре зоны предразрушения [145], когда а:0 « а и й 0 (а) « й (а).

Предельно-равновесное состояние такого тела можно опреде­ лить на основании критерия Гриффитса — Ирвина (1.1). Однако указанный критерий только интегрально учитывает свойства зоны предразрушения и не дает никакой информации о ее размерах, что необходимо для математического описания условий автомодель­ ности этих зон. Решение этой задачи целесообразно осуществить в рамках следующей расчетной схемы [1, 2].

Пусть внешняя нагрузка задана таким образом, что напряжен­ ное состояние тела симметрично относительно плоскости располо­ жения трещины, а материал тела считается упруго-пластическим, подчиняющимся условию пластичности Треска — Сен-Венана. В силу условий автомодельности зон предразрушения и симметрич­ ности напряженного состояния относительно плоскости располо­ жения трещины в достаточно малой окрестности ее контура будет осуществляться условие плоской деформации, которое описывает­ ся коэффициентом интенсивности напряжений К г. Как показыва­ ют экспериментальные данные [55, 163, 187], в случае плоской де­ формации пластические зоны локализуются главным образом вдоль некоторого слоя, направленного примерно под углом 45°—•

— 72° к плоскости расположения трещины. Поэтому зону предраз-

рушения в окрестности контура трещины представим следующим образом. Под уг­ лом 45— 72° к плоскости располо­ жения трещин размещены (рис. 3) достаточно узкие зоны пластического скольжения длиной Z. При решении задачи эти зоны рассматриваются как трещины-разрезы, к поверхности ко­

торых приложены касательные на-

1 пряжения т = -у огт, а скачок нормаль.

ных перемещений равен нулю.

 

 

Предельно-равновесное состояние те­

 

ла определяется KPT-критерием,

т. е.

Рис. 3. Схематическое изо­

на основании соотношений

 

 

бражение полос пластич­

2у[Д„(а*),

а*, 0] = б к;

 

ности в зоне предразруше-

^а[7?о((Х)>

O&t 0]а=ос* = 0,

(1.6)

ния.

 

где у Ш0 (а), а, 0] — нормальная составляющая и (г, a, z) вектора смещения точек тупиковой части трещины; а* — координатный угол точек контура трещины, в окрестностях которых наступает предельно-равновесное состояние.

Упруго-пластическая задача для неограниченного тела, ослаб­ ленного полубесконечным разрезом, в окрестности контура ко­ торого осуществляется условие плоской деформации, рассмотрена в работах [66,143]. Используя результаты этой работы, находим

2v [Д0(а), а, 0] =

0,444

(* V ) - ,

(1.7)

 

 

UT

 

где v — коэффициент Пуассона;

Е — модуль Юнга;

ат— пре­

дел текучести материала при растяжении.

Протяженность пластической области I в этом случае (см. рис.

3) вычисляют по формуле

 

г?-2

(1.8)

I =0,578 — 4 - .

Подставляя выражение (1.7) в соотношения (1.6), для определе­ ния предельного значения внешних нагрузок получаем уравнения

 

0,444 (1 — v2) К.\, =

8котЕ;

dKf (a)

=

о,

(1.9)

 

да

 

 

 

 

a =a *

 

 

где

Ки — значение

коэффициента интенсивности

напряжений

Кг при Р = Р% и a = a*.

 

случая

макротрещин

при

Поскольку для рассматриваемого

Р =

Р* должно быть

К\¥ =

Kic, то

на основании соотношений

(1.1), (1.8) и (1.9) между характеристиками трещиностойкости для упруго-пластических тел (К{С, бк) и критической длиной

2 7-281

17

(Z = Z*) пластической зоны можно установить такие соотношения!

K ic = 2,33loTV I*.

(1. 11)

При экспериментальном определении характеристики трещиностойкости материала К\с с использованием критериального урав­ нения (1.1) необходимо выбирать такие геометрические параметры образцов, для которых выполняются условия автомодельности воны предразрушения. Если длина Z* пластической зоны может быть установлена экспериментально, то условие (1.11) можно применить для нахождения К\с или проверки достоверности оп­ ределения этой характеристики другим путем, т. е. выполнения условий автомодельности.

Упомянутые три типа моделей лежат в основе формирования

и развития современной механики разрушения твердых тел. Мо­

дели второй и третьей групп еще не получили необходимого завер­

шения, однако они представляют значительный интерес, так как

охватывают самый широкий класс конструкционных материалов.

Модели такого типа составляют теоретическую перспективу раз­

вития современной механики разрушения. В плане практических

приложений указанных теоретических моделей важна разработка

методов определения характеристик стойкости материала против

хрупкого разрушения (К\С1 бк, у), а также

решение конкретных

задач для построения Ки и 6к-тарировок с

целью эффективного

определения этих характеристик.

3. Некоторые соотношения теории упругости

При рассмотрении задач механики хрупкого разрушения важным и достаточно трудным этапом является определение и анализ напряженно-деформированного состояния в упругом трехмерном теле, ослабленном дефектами типа трещин. Это объясняется тем, что в случае трехмерных задач отсутствует такой единый и эффективный аналитический аппарат, как метод Колосова — Мусхелишвили [72] в плоской теории упругости.

Исходные уравнения пространственных задач теории упругос­ ти и основные методы их решения сформулированы в ряде учебни­ ков и монографий по теории упругости (см., например, [59, 63, 78, 130]). Ниже выводятся лишь некоторые соотношения статики и динамики упругого тела, необходимые в дальнейшем для ис­ следования предельного равновесия квазихрупкого цилиндра, ослабленного внешней кольцевой трещиной.

1. Исходные уравнения статики трехмерного упругого тела. Пусть упругое тело занимает некоторый объем и на его поверх-

ах cos (п, х) +

Хуу cos (re, у) +

ххг cos (re, z) -- FПХ>

 

T

COS (n, x) +

 

COS (n, у) +

Tyz cos (n, z) = F^y\

(1 . 12)

Xy

 

Oy

 

 

 

Txzcos (n, я) +

Ty2cos (n, y) +

Gzcos (n, z) =

 

 

или

 

 

 

 

 

 

 

U =

U(S> V = V,( S )

w = w(S)

 

(1.13)

где ax, ay,

a2, тху, тХ2, ту2 — компоненты тензора

напряжении

в прямоугольной системе декартовых координат Oxyz; F^L F^y,

F $ — компоненты заданного вектора напряжения, действующего на поверхности тела S; п — внешняя нормаль к поверхности S; и, у, и; — компоненты (проекции) вектора упругих перемещений

и соответственно на оси Ox, Oy, Oz; u(S), y(S), iy(S) — значения функ* ций и, v, w в точках поверхности S.

Решение задачи теории упругости для такого тела сводится

к отысканию вектор-функции и упругих перемещений, которая при отсутствии массовых сил должна удовлетворять заданным на по­ верхности тела граничным условиям (1.12), (1.13) и векторному уравнению равновесия

 

grad div и + (1 — 2v) Аи = О,

(1.14)

где А =

— трехмерный оператор

Лапласа.

Зависимость между функциями и, v, w и компонентами сгх, (Ту, а2, тХу, тХ2, Ту2 определяется на основании закона Гука, со­ отношений Коши и выражается такими формулами:

О

ъ

ии . UU , UIV

— относительное

Здесь ц — модуль сдвига;

0 =

 

объемное расширение.

Если на поверхности тела S заданы условия (1.12), то опреде­ ление упругого равновесия тела составляет первую основную гра­ ничную задачу. Кроме первой основной задачи в теории упругос­ ти значительный интерес представляет и вторая основная гранич­ ная задача, т. е. определение упругого равновесия тела, когда на его поверхности S заданы условия (1.13). Наконец, во многих случаях (контактные задачи, задачи теории трещин и т, д.) большое

значение имеет основная смешанная граничная задача, т. е. когда на одной части поверхности тела заданы смещения, а на осталь­ ной — напряжения.

Вдальнейшем, при рассмотрении предельного равновесия тел

стрещинами, будут необходимы решения только основной сме­ шанной граничной задачи теории упругости.

Общее решение уравнения (1.14) было получено П. Ф. Папковичем [78], а позднее Г. Нейбером [74] в следующей форме:

2|хи = grad (ср 4- гЧ?) — 4 (1 — v) Y,

(1.16)

где 4я — гармоническая вектор-функция; ф — гармоническая ска­

лярная функция; г — радиус-вектор точки тела.

Решение Папковича — Нейбера (1.16) включает четыре скаляр­ ные функции, а именно: ф и три скалярные компоненты Т 1?

Ч'з вектор-функции Y. Из решения (1.16) легко вывести более прос­ тое решение для случая, когда касательные напряжения равны нулю в плоскости, которую можно выбрать в качестве координат­ ной плоскости z = 0. При этом искомое решение будет содержать только одну гармоническую функцию. Действительно, учитывая выражения (1.15) и (1.16), касательные напряжения тхг и туг можно представить в следующем виде:

Ъ*

+ 2 ( 1 — V)

 

(1.17)

 

+ 2 (1 — v)

 

 

где

 

dVj

 

дУ2

дЧз

F = ( l - 2 v ) V 3- ^ _ x

- у -

dz

dz

dz

 

 

 

Чтобы напряжения (1.17) были равны нулю в плоскости 2 = 0,

необходимо в (1.16)

положить

 

*1 =

^ = 0; ( l - 2 v ) Y 3 = - g - .

(1.18)

Следовательно, при выполнении равенств (1.18) функции ф и

Т* можно представить в виде

 

 

 

 

 

Ф =

(1 — 2v) X;

 

f =

(o,0,

- g - ) ,

 

где X — гармоническая функция.

 

 

 

Тогда решение уравнения (1.14) будет иметь вид

 

2|Ш= z

d 4

+

(1 — 2v) dX

;

 

 

 

dxdy

 

 

 

 

 

2\iv = 2

d 4

+

d - * ) - £

;

(1.19)

dydz

о

. 6

4d 4

0 /,

. d%

 

 

 

= z ~%dz2----- 2 (! ~

v) “dT •

 

Рис. 4. Изображение компонентов вектора пе­ ремещения в цилиндрической системе коорди­ нат.

Подставляя выражения (1.19) в со­ отношении (1.15), для определения на­ пряжений получаем формулы

ох = z

dx2dz

,

и л д дгх

ач

 

 

 

 

 

 

o v = z

dPx

+

(1 — 2v)

— 2v

 

ay*dz

^

 

04

dz2

Or = z- d 4

~ " dz3

d3X

Tyz Z dydz2 *

84

Txy — Z - dxdydz

дЧ dz2 *

03X TX2 = 2 ■dxdz2

84

+ ( l - 2v) дхду

(1 .20)

В цилиндрической системе координат г, a, z (рис. 4) при осевой симметрии напряженно-деформированного состояния решение (1.19) будет иметь такую форму:

2[iur = z

вч

д%

drdz

■ + ( I - 2V) TarT ;

2)xuz = z

дЧ

( 1. 21)

•2(1 - v ) - g - .

az2

Эти соотношения приводят в свою очередь к следующим выраже­ ниям для компонентов тензора напряжений:

.ь - b

II

 

Q N

 

II

Trz =

Z

d 4

+

( 1

- 2

дЧ

■2v

дЧ

dr2dz

 

 

 

дг2

 

dz2

дЧ

 

дЧ

 

 

 

(1 .22)

dz3

 

dz2

 

 

 

 

 

 

d 4

»

lar ^az — 0*

 

 

drdz2

 

 

 

 

T

 

 

 

 

 

 

 

 

 

d 4

Из формул (1.22) следует, что если производная - ^ 2:- остает­

ся конечной при z = 0, то тгг =

0 при z =

0, а нормальные напря­

жения и перемещения в плоскости z = 0

определяются так:

= ~ (“Й5") 1=0

= — 2 (1

V) ( - £ - ) г=0 •

( L 2 3 )

При определении упругого равновесия цилиндра с внешней кольцевой трещиной в дальнейшем будем использовать соотноше­ ния (1.21) — (1.23),

2. Распространение упругих волн. Теперь коротко рассмотрим случай, когда приложенные нагрузки изменяются со временем так, что в теле возникают динамические напряжения. При этом урав­ нение движения упругого тела можно записать в векторной форме так;

Аи + ■

grad div и = ± - g

- ,

(1.24)

где и — вектор перемещения; Л — плотность

упругой

среды;

t — время.

 

 

 

Вдальнейшем при решении динамической задачи для цилиндра

свнешней кольцевой трещиной необходимыми будут решения урав­ нения (1.24) в случаях плоской деформации и осевой симметрии тела.

Если в теле компонента смещения w равна нулю, а компоненты

и, v зависят только от х , у , но не от z (плоская деформация, парал­ лельная плоскости Оху), то уравнение (1.24) для этого случая в скалярной форме можно записать так:

д2и_

д2и

 

 

1

/

д2и

 

д2у

\ _

_Х_ д2и в

 

дх2

ду2

 

1— 2v

^ дх2

* дхду J

ц

dt2 *

(1.25)

д2и

д2У

,

 

1

/

д2и .

д2и

\ _

X

д2и

 

~д&

ду2

 

\— 2v

\

ду2

 

дхду ]

ц

dt2

 

Решение системы

дифференциальных

уравнений (1.25)

ищем

в форме

 

дфх ,

 

 

 

 

 

 

 

 

 

 

и =

 

ду

V =

 

 

ду

 

дх

(1.26)

 

 

дх

"*

 

 

 

 

 

где фх, 4j5x — неизвестные функции.

 

 

находим

соотношения

Подставляя

(1.26)

в уравнения (1.25),

д

2 —2у

/ а2фх

, _£^фх_\____Х_ а2фх 1

 

 

дх

1— 2v

^

дх2

 

ду2

J

 

р

а<2 J

 

 

а

[

ах2

 

____ *

*4* 1

п.

 

 

~

ау

[

а^2

р

 

л 2 J

 

(1.27)

а

г 2 —2v /

а%

а2фХ \

 

\

а2фХ i

 

ду

[ 1— 2v (а х 2

ар2

j

 

р

аг2

J

 

 

____а

г а2фх

 

а2фх

х

а2фх- i Л

 

 

 

cte

 

дх2

 

ду2

 

ус

dt2

J

 

 

 

Эти выражения показывают, что соотношения (1.26) представляют решение уравнении движения (1.25) в том случае, если <px и выбраны таким образом, что они являются решениями уравнений

л /«5PL .

д2<РА

d2(Pi

(1.28)

с1{ дх2 ^

ду2 )

&2

 

о /

\

 

 

2^ дх2

J