Известия РАН. Механика твердого тела, 2019, № 5, стр. 70-85

О МАСШТАБНОЙ ЗАВИСИМОСТИ ТРЕЩИНОСТОЙКОСТИ СТРУКТУРИРОВАННЫХ ТЕЛ (НА ПРИМЕРЕ МОРСКОГО ЛЕДЯНОГО ПОКРОВА)

Н. М. Осипенко *

Институт проблем механики им. А.Ю. Ишлинского РАН
Москва, Россия

* E-mail: osipnm@mail.ru

Поступила в редакцию 08.02.2019
После доработки 10.02.2019
Принята к публикации 19.02.2019

Полный текст (PDF)

Аннотация

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

Ключевые слова: разрушение, трещина, структура, сжатие

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

Одним из базовых положений линейной механики разрушения является допущение о малости концевой области трещины (в частности, трещины нормального разрыва) по сравнению с размерами трещины, и ее автономности. Это предполагает независимость механизма локального деформирования и разрушения в этой области в состоянии предельного равновесия (во внутренней задаче теории трещин) от локальных условий нагружения и геометрии внешней задачи, масштаб которой задается размером трещины и границами тела. Такая постановка задачи дает возможность характеризовать напряженное состояние в вершине трещины нормального разрыва одним параметром – коэффициентом интенсивности напряжений, предельная величина которого (KIC – для трещин нормального разрыва) характеризует сопротивляемость среды развитию трещин и является фиксируемым свойством материала. Вместе с тем существуют ситуации, в которых такие условия соблюдаются не для всех сценариев разрушения. К ним относится, в частности, изменение масштаба разрушения. Заметим, что влияние масштабного фактора на параметры прочности для таких материалов, как горные породы, бетон и лед являются объектом внимания достаточно давно [1]. Интерес к этой теме в последнее время связан, в том числе, с противоположными тенденциями в масштабных эффектах для различных свойств материалов – уменьшением удельной прочности при росте размеров нагружаемого тела и увеличением эффективной трещиностойкости в той же ситуации. Основные направления этих исследований обозначены в [2]. Отметим также недавнюю работу [3], в которой показанo масштабное изменение трещиностойкости снега.

2. Масштабные вариации при определении трещиностойкости морского льда. Рассмотрим имеющиеся данные об изменении эффективной трещиностойкости, связанные с изменением масштаба, на примере механики разрушения морского льда. Морской лед со столбчатозернистой структурой и порами, заполненными рассолом (тип S2 [4]) имеет несколько характерных масштабных уровней, которые влияют на развитие трещин в морском ледяном покрове, начиная с размера зерна и расстояния между капиллярными порами и дренажными каналами до размеров ледовых образований. Напряженное состояние в структуре льда при сжатии в масштабе окрестности одиночного капилляра близко к таковому в пластине с отверстием в однородно нагруженном упругом теле (плоская деформация). Такая система может служить отдельным структурным элементом, способным преобразовать приложенные к телу внешние напряжения. Модель разрушения пористого тела была в основных чертах предложена ранее [57]. Для большинства однородных и слабопористых материалов, для которых такое представление справедливо, характерно развитие трещин нормального разрыва вдоль оси нагружения на контуре пор [8]. Очаг разрушения образуется в области локальных растягивающих напряжений в окрестности контура отверстия радиуса (R) в направлении вдоль оси сжатия, что соответствует положению максимума нормальных упругих напряжений на контуре отверстия. В результате при внешнем сжатии в этой области образуется устойчиво развивающаяся трещина нормального разрыва (l). Разрушение происходит при слиянии трещин, растущих от соседних пор, когда величина L = l + R близка к половине среднего расстояния между порами. При таком условии в капиллярной структуре средняя пористость n может быть связана с соотношением размеров структуры.

(2.1)
$R{\text{/}}L\sim {{n}^{{1/2}}}$

Для оценки эффективной трещиностойкости морского льда в масштабе капиллярной пористой структуры ранее была предложена модель разрушения [9]. По условиям предельного равновесия симметричных микротрещин выходящих на контур одиночной поры взаимосвязь между прочностью при одноосном сжатии поперек оси капилляров и капиллярной пористостью может быть представлена в виде

(2.2)
${{\sigma }_{1}} = - \frac{{K_{{IC}}^{*}}}{{\sqrt {\pi L} }}\frac{1}{{1.1{{n}^{{1.65}}}\sqrt {(1 - {{n}^{{0.5}}})} }}$
где $K_{{IC}}^{*}$ – эффективный критический коэффициент интенсивности напряжений – трещиностойкость льда в масштабе единичной поры, n – капиллярная пористость.

Сравнение полученного выражения с эмпирическим соотношением, приведенным в [10] и [11] для столбчато-зернистого морского льда, сжатого поперек оси капилляров, показывает, что уровень эффективных значений $K_{{IC}}^{*} \approx 5\;кРа \cdot \sqrt м $.

Обратимся далее к экспериментальным данным критического коэффициента интенсивности напряжений – эффективной трещиностойкости льда ${{K}_{{1C}}}$ – при увеличении масштаба разрушения:

– При использовании малых (лабораторных) образцов (характерный размер L ~ ~ $0.05{\kern 1pt} - {\kern 1pt} 0.1$ м [9, 12, 13]) при температуре ∼ –10°С величина ${{K}_{{1С}}} \approx 60{\kern 1pt} - {\kern 1pt} 80\;кПа\sqrt м $. Образцы в этих испытаниях имели конфигурации, рекомендованные для конструкционных материалов.

– По данным испытаний, проведенных в [13, 14, 16], на фрагментах морского ледяного покрова толщиной до ∼1.8 м, при развитии трещин в льдинах с характерным размером в плане от 0.5 до 3 м величина ${{K}_{{1C}}}$ возрастает до значения ${{К}_{{1С}}}_{{ef}} \approx 250\;кПа\sqrt м $.

На рис. 1 в логарифмических координатах показана совокупность данных о влиянии масштаба разрушения на изменение эффективной трещиностойкости морского льда (□ – данные авторов [9, 12], – сводные данные из работ [14, 15]). Также отражено ограничение величины эффективной трещиностойкости при наибольших размерах образцов (в ситуации, когда размер плоского образца значительно больше его толщины ($L \gg h$)). При L > 10 м и $h \leqslant 0.8$ м трещиностойкость практически не меняется [13]. Далее будем интересоваться диапазоном условий, при которых проявляется зависимость эффективной трещиностойкости ${{K}_{{1C}}}$ от размера образца. Линейная аппроксимация соотношения $\ln {{K}_{{{\text{IC}}}}}$ и $\ln L$ на участке $L < 10$ м (прямая на рис. 1) дает соотношение, связывающее трещиностойкость и характерный размер в виде

(2.3)
$\ln {{K}_{{1C}}} \approx 5.1 + 0.5\ln L$

или

(2.4)
${{K}_{{1C}}} \approx m{{L}^{{0.5}}}[МПа\sqrt м ]$, $m = 0.164$ [МПа]
Рис. 1

Базовые эксперименты в масштабе лабораторного образца и выше, на которых была найдена зависимость эффективной трещиностойкости от размера образца [13, 14] выполнялись с соблюдением рекомендуемого геометрического подобия образцов (в плане) при увеличении размеров. При этом был отмечен ряд важных деталей. Увеличение эффективной трещиностойкости K1С сопровождалось ростом предельного раскрытия вершины надреза (радиус исходного надреза оставался постоянным) и множественным нарушением связей в области концентрации растягивающих напряжений [1618]. Это указывает на формирование увеличивающейся области предразрушения – концевой зоны трещины – в ситуации возрастания измеряемого параметра. Детали механизма такого явления представляют предмет отдельного исследования (относятся к внутренней задаче механики разрушения). Обсуждение этого вопроса для льда можно найти в [15]. Соглашаясь с выводами [15] о неправомерности предположения об одинаковости механизмов предразрушения, определяющих энергозатраты в концевых зонах при распространении трещины при такой разнице масштабов разрушения, можно по приведенным экспериментальным данным заключить, что имеется аналогия реакции материала на некоторый определяющий параметр нагружения в разных масштабах. Здесь уместно отметить, что этот вопрос нуждается в уточнении, в том числе и в силу недостаточности данных в диапазоне промежуточных размеров объектов (1–100 мм).

Отметим также возможный вариант обобщенного подхода (со стороны внешней задачи) в виде связи уровня трещиностойкости с размером образца и размером концевой области трещины. В силу такой связи при изменении размера эффективной концевой зоны (зоны предразрушения) по мере роста трещины в предельно-равновесном режиме в рамках внешней задачи для характеристики напряженного состояния в этой области перестает строго действовать асимптотика решений задачи теории упругости для малых расстояний от вершины трещины. Поэтому изменяются условия разрушения. Рассмотрим для немалой концевой области некоторые из них, полагая, что в рамках внешней задачи сохраняется упруго-хрупкий характер разрушения.

Возрастание размеров концевой зоны при увеличении размеров образца подтверждается измерениями раскрытия в вершине надреза (CTOD) к моменту страгивания трещины. Так по данным [13] при увеличении размера компактного образца морского льда от 0.5 до 3 м величина CTOD изменялась от 9–12 мкм до 36–45 мкм.

Для оценки характерного размера концевой обрасти в ледяном покрове в рамках подхода внешней задачи теории трещин воспользуемся известным соотношением для трещины в предельно-равновесном состоянии [19].

(2.5)
$r \approx \frac{1}{{2\pi }}{{\left( {\frac{{{{K}_{{1C}}}}}{{{{\sigma }_{t}}}}} \right)}^{2}}$

Напомним, что характерный линейный размер концевой области может быть использован и в качестве характеристики трещиностойкости (параметр L1C [20]).

(2.6)
L1C$ = {{\left( {\frac{{{{K}_{{1C}}}}}{{{{\sigma }_{t}}}}} \right)}^{2}}$

Для оценки локальной прочности при изменении масштаба концевой зоны воспользуемся соотношением Базанта [21, 22]

(2.7)
${{\sigma }_{t}} = \frac{\Sigma }{{{{{(1 + L{\text{/}}B)}}^{{1/2}}}}}$
где $\Sigma $ и В – постоянные (по данным [13] для морского льда $\Sigma \approx 0.59$ [МПа], $B \approx 0.37$ [м]).

Связывая размер L с размером образца, и учитывая (2.4), получим из (2.5) и (2.7) оценку размера концевой обрасти

(2.8)
$r \approx 0.033L(0.37 + L)$

Вид этой функции показан на рис. 2 (размеры даны в метрах). Можно видеть, что при характерных размерах ледяных образований порядка 1 м и выше ожидаемый эффективный размер концевой зоны, оставаясь значительно меньше размеров образца в плане, значительно больше, чем размеры структурных элементов материала (зерен, пор) и порядка толщины образцов. Поэтому при развитии трещин нормального разрыва в концевых зонах в этих масштабах могут иметь значение промежуточные структуры разрушения малого масштаба, образующиеся при взаимодействии между элементами исходной структуры льда, а при $L > 10$ м размер концевой зоны сопоставим с толщиной образца и в этой области меняется напряженное состояние.

Рис. 2

Поясним это на примере развития трещины в балке под действием изгибающего момента. Согласно [24] для этого случая нагружения имеем

(2.9)
$\begin{gathered} {{K}_{1}} = \frac{{6M}}{{t{{H}^{2}}}}\sqrt {\pi \ell } f\left( {\ell {\text{/}}H} \right) \\ f\left( {\ell {\text{/}}H} \right) = 1.22 - 1.40(\ell {\text{/}}H) + 7.33{{(\ell {\text{/}}H)}^{2}} - 13.08{{(\ell {\text{/}}H)}^{3}} + 14.0{{(\ell {\text{/}}H)}^{4}} \\ \end{gathered} $
где $\ell $ длина трещины, Н – толщина балки, t – ее ширина, М – изгибающий момент.

Полагая, как обычно, условие предельного равновесия ${{K}_{I}} = {{K}_{{IC}}}$ учтем взаимосвязь KIC и масштабного эффекта льда в форме (2.4). Будем считать в качестве параметра, определяющего масштаб процесса разрушения, размер неповрежденной части изгибаемой балки L$\left( {H - \ell } \right)$ на продолжении надреза. Критическая величина изгибающего момента определится соотношением

(2.10)
$M = \frac{{{{K}_{{IC}}}t{{H}^{2}}}}{{6\sqrt {\pi \ell } f(\ell {\text{/}}H)}} \approx \frac{{mt{{H}^{2}}\sqrt {H - \ell } }}{{6\sqrt {\pi \ell } f(\ell {\text{/}}H)}}$

На рис. 3 показана оценка критической величины изгибающего момента по (2.10) для балки шириной 1 м и толщиной 1.5 м с трещиной, имеющей размер $\ell $ (кривая 1). Для сравнения приведена оценка для варианта с уровнем трещиностойкости (KIC = = 200 кПа м–0.5 – кривая 2) и трещиностойкостью, измеренной в масштабе лабораторных образцов (KIC ∼ 80 кПа м–0.5– кривая 3).

Рис. 3

Кривая 1 занимает промежуточное положение. Поэтому, например, оценка глубины необходимого начального надреза для обламывания ледяной балки по данным лабораторного определения KIC (кривая 3) дает заниженную величину эффективной трещиностойкости. Соотношение критического изгибающего момента в варианте учета масштабного изменения трещиностойкости по (2.10) и при постоянной ее величине, соответствующей максимальному уровню KIC для рассматриваемого примера показано на рис. 4 для балки толщиной 1.5 м.

Рис. 4

Эксперименты по определению трещиностойкости горных пород, льда и других хрупких материалов обычно проводятся по методике с использованием локально нагружаемых “компактных образцов”, в том числе и образцов большого размера [13], либо при изгибе образцов – балочек. При этом обычно соблюдаются рекомендации по определенному соотношению размеров исходного надреза (зародышевого нарушения) и размеров образца.

Поэтому при пропорциональном увеличении масштаба – размеров образцов и исходных надрезов – фактические условия нагружения в области концевой области исходного надреза (трещины), такие как градиент внешних напряжений на продолжении исходного надреза, изменяются. В этом случае можно ожидать влияния этого фактора на эффективную трещиностойкость по мере развития индивидуальной трещины в образце, когда эффективный масштаб будет задаваться расстоянием до свободной границы. Отражением такого влияния может быть и отмеченное в [23] влияние границ нагружаемой области на измеряемые параметры. В такой ситуации представляет интерес учет влияния градиента напряжений от внешних нагрузок на параметры и механизм предразрушения в концевой области трещины.

3. Влияние градиента напряжений. Качественное объяснение эффекта для трещин с немалой концевой областью представляется следующим. Размер концевой области, в которой происходит затрата энергии при подготовке локального разрыва испытывает влияние вариаций поля внешних напряжений, приходящихся на эту область. В частности, при резком изменении внешних напряжений (возрастании) по мере приближения к вершине исходного нарушения в концевой области трещины, локальные процессы предразрушения стягиваются к этой вершине. Происходит уменьшение размера концевой области и, соответственно, уменьшение потребления энергии в окрестности вершины надреза. По мере выполаживания внешних напряжений размер концевой области растет, приближаясь к асимптотике, задаваемой линейной механикой разрушения (вариант с возрастанием напряжений на продолжении разреза в квазистатической постановке практически не реализуется).

Роль локальной неоднородности напряженного состояния в определении критических условий разрушения конструкций и природных объектов давно привлекает внимание. Среди подходов к этому вопросу, обзор которых можно найти, например, в [24] выделим область интереса к влиянию градиентов напряжений. Так в [25] показано, что неравномерность напряженного состояния приводит к снижению разрушающей способности напряжения в области максимальных значений. Согласно градиентным критериям для определения разрушающей нагрузки нужно сравнивать не максимальное значение первого главного напряжения ${{\sigma }_{1}}$ с пределом прочности материала ${{\sigma }_{в}}$ (принятого в качестве эквивалентного), а максимальное значение эффективного напряжения ${{\sigma }_{{ef}}}$, которое при неоднородном напряженном состоянии меньше эквивалентного. Поэтому эффективное напряжение ${{\sigma }_{{ef}}}$ пропорционально ${{\sigma }_{1}}$ и зависит от степени неравномерности поля напряжений в окрестности рассматриваемой точки тела и неоднородности материала. К настоящему времени вопрос о влиянии такого эффекта не может считаться исчерпанным. Коротко остановимся на известных подходах. В [26] предложена форма зависимости

(3.1)
${{\sigma }_{{ef}}} = {\text{ }}{{\sigma }_{1}}{\text{/}}f\left( {{{g}_{1}},{\text{ }}{{L}_{1}},{\text{ }}\beta } \right){\text{ }} < {\text{ }}{{\sigma }_{1}}$, ${{\sigma }_{{ef}}} = {{{{\sigma }_{1}}} \mathord{\left/ {\vphantom {{{{\sigma }_{1}}} {(1 - \beta + \sqrt {{{\beta }^{2}} + {{L}_{1}}{{g}_{1}}} )}}} \right. \kern-0em} {(1 - \beta + \sqrt {{{\beta }^{2}} + {{L}_{1}}{{g}_{1}}} )}}$
где ${{g}_{1}}$ – градиент первого главного напряжения.

В отличие от классических однопараметрических критериев градиентный критерий разрушения содержит еще параметры L1 и $\beta $. Параметр L1 находится из условия соответствия градиентного критерия с линейной механикой разрушения и сводится к известной задаче об одноосном растяжении плоскости с эллиптическим отверстием (перпендикулярно большой оси эллипса). Связь параметра L1 с критическим коэффициентом интенсивности напряжений ${{K}_{{{\text{I}}\,c}}}$ и пределом прочности ${{\sigma }_{в}}$ дается в виде

(3.2)
${{L}_{1}} = \frac{2}{\pi }{{K_{{{\text{I}}C}}^{2}} \mathord{\left/ {\vphantom {{K_{{{\text{I}}C}}^{2}} {\sigma _{в}^{{\text{2}}}}}} \right. \kern-0em} {\sigma _{в}^{{\text{2}}}}}$

Таким образом, параметр L1 выражается через известные характеристики материала и приравнивается к известному в механике разрушения критическому размеру дефекта в виде трещины Гриффитса. Безразмерный параметр β, изменяющийся от нуля до единицы, учитывает квазихрупкий характер разрушения: β = σв/(E${{\varepsilon }_{*}}$), где ${{\varepsilon }_{*}}$ – деформация в момент разрушения при одноосном однородном растяжении; E – модуль упругости.

В другом варианте градиентного подхода [27, 28] локальная прочность материала предполагается зависящей от размера зоны концентрации напряжений Le. Если размер Le достаточно велик по сравнению с размерами структурных составляющих материала, то величина локальной прочности мало отличается от величины предельного напряжения σ0, определенной в условиях однородного распределения напряжений. Наоборот, если Le сопоставим с размерами структурных элементов, их влияние на локальную прочность становится заметным. Причем это влияние тем больше, чем меньше размер Le по отношению к характерному размеру структуры материала L0. Таким образом, локальная прочность материала зависит не просто от размера зоны концентрации напряжений Le, а от соотношения L0/Le, которое характеризует масштаб в рассматриваемой задаче.

Приведем некоторые оценки, предполагая основным влияющим фактором градиент напряжений на потенциальной траектории трещины. Вернемся к компактному образцу, служившему прототипом образца для натурных испытаний в работах со льдом [13, 14]. Свяжем ${{K}_{{{\text{IC}}}}}$ с уровнем градиента внешнего поля напряжений в области вершины трещины. В компактном образце, нагружаемом сосредоточенными силами на внешнем конце надреза градиент напряжений на продолжении исходного надреза создается в основном изгибающим моментом в плоскости образца. В балочном приближении компактный образец можно представить в виде двух полубалочек, передающих изгибающий момент в область прямолинейной траектории трещины на продолжении исходного надреза.

(3.3)
$M \approx PL{\text{*}}$
где Р – сосредоточенное усилие от нагружающего устройства, $L{\text{*}}$ – плечо (в экспериментах на квадратном компактном образце эта величина близка к сумме длины исходного надреза и полуразмера оставшейся части образца. $L* \approx (\ell + L){\text{/}}2$, где $L$ – полный размер образца, $\ell $ – длина надреза (трещины).

По решению задачи теории упругости [29] для такого образца

(3.4)
$P = \frac{{{{K}_{I}}t\sqrt {L{\text{*}}} }}{{f(\ell {\text{/}}L*)}}$

В случае $\ell {\text{/}}L$ ∼ 0.5 величина $f(\ell {\text{/}}L)$ ∼ 9.66.

Изгибные напряжения на траектории трещины в том же приближении

(3.5)
$\sigma = \frac{{12My}}{{t{{{(L - \ell )}}^{3}}}}$

Градиент внешних напряжений

(3.6)
$\sigma {\text{'}} = \frac{{d\sigma }}{{dy}} = \frac{{6P(L + \ell )}}{{t{{{(L - \ell )}}^{3}}}} = \frac{{6(L + \ell ){{K}_{I}}\sqrt L }}{{f(\ell {\text{/}}L){{{(L - \ell )}}^{3}}}}$

Для варианта $\ell {\text{/}}L$ ∼ 0.5

(3.7)
$\sigma {\text{'}} \approx \frac{{7.45{{K}_{I}}}}{{{{L}^{{3/2}}}}}$

Согласно (2.4) и (3.7), в предельно-равновесном состоянии для морского льда получим

(3.8)
$\sigma ' \approx 1.2{\text{/}}L\;[{\text{MПа/м}}]$

Отсюда, используя (2.5), в рамках внешней задачи можно, таким образом, представить масштабную зависимость трещиностойкости морского льда в виде функции градиента внешних напряжений на продолжении трещины в момент страгивания трещины

(3.9)
${{K}_{{1C}}} \approx m{{L}^{{0.5}}}\sim m{{\left( {1.2{\text{/}}\left| {\sigma {\text{'}}} \right|} \right)}^{{0.5}}} = n{{\left( {\left| {\sigma {\text{'}}} \right|} \right)}^{{ - 0.5}}}$, $n \approx 0.2$ [(МПа)1.5]

Можно отметить, что если в (3.2) при определении характеристического размера L1 в качестве ${{K}_{{{\text{IC}}}}}$ используется эффективная трещиностойкость материала по (3.9), влияние градиента напряжений в (3.1) нивелируется. Механизм влияния градиента внешних напряжений на условия предельного равновесия в концевой области можно проиллюстрировать, представляя его в виде вклада вариации внешних напряжений в окрестности вершины трещины в оценку размеров концевой области. На ее границе будем по-прежнему считать напряжение близким к предельному ${{\sigma }_{в}}$, так что условие, определяющее размер концевой зоны запишется в виде

(3.10)
${{\sigma }_{в}} \approx {{f}_{1}}({{K}_{{IC}}},r) + {{f}_{2}}(\operatorname{var} \sigma ,r{\text{),}}\quad {{f}_{1}}({{K}_{{IC}}},r) \approx \frac{{{{K}_{{{\text{IC}}}}}}}{{\sqrt {2\pi r} }},$

Для случая немалой концевой области при постоянном градиенте напряжений можно записать приближенное соотношение ${{f}_{2}}(\operatorname{var} \sigma ,r{\text{)}}\sim \sigma {\text{'}}r$.

Отметим, что такая форма условия предельного равновесия по смыслу близка к предложенной в [30] модификации соотношения (2.7), в которое добавлен свободный член, отвечающий за дополнительные затраты энергии при изменении масштаба. Соотношение (3.10) при высоких уровнях градиента напряжений и размеров концевой области, (то есть в ситуации преобладания процессов деформирования, связанных с искажением поля напряжений в концевой области) стремится к асимптотике

(3.11)
$r\sim {{\sigma }_{в}}{\text{/}}\sigma {\text{'}}$

Для такого варианта уровень эффективной трещиностойкости определится, как

(3.12)

Отмечая совпадение формы взаимосвязи эффективной трещиностойкости с градиентом напряжений в (3.9) и (3.12) можно получить для морского льда оценку $n \approx \sqrt {2\pi \sigma _{в}^{3}} \approx 0.2$ [(МПа)1.5], и ${{\sigma }_{в}} \approx 0.185$ мПа, что близко к уровню прочности морского льда на разрыв [31].

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

Эксперименты проводились на компактных образцах размером (49 × 49 мм в плане). Толщина гипсовых образцов – 15 мм, парафиновых – 10 мм. Изменялась длина исходного надреза в серии экспериментов при сохранении остальных размеров образцов, что вызывало различие градиента напряжений при страгивании трещин. Результаты измерений показаны на рис. 5 (обозначения: o – гипс, + – парафин). Можно видеть устойчивую тенденцию снижения эффективной трещиностойкости при увеличении градиента напряжений. На рис. 6 показаны относительные изменения параметра ${{K}_{{1C}}}{\text{/}}\sigma {\text{'}}$ при росте $\sigma {\text{'}}$ для тех же образцов модельных материалов в сравнении с оценкой для льда (3.9) по данным, представленным выше (пунктир). Здесь ${{({{K}_{{{\text{IC}}}}}{\text{/}}\sigma {\text{'}})}_{0}}$ отражает значения параметров для экспериментов при минимальном уровне градиента напряжений в эксперименте. Можно видеть, близкие тенденции понижения относительного уровня значений при росте градиента напряжений для всех испытанных материалов, что позволяет считать значимой принятую взаимосвязь характеристик трещиностойкости структурированных (зернистых, пористых) материалов с масштабом градиента внешних напряжений.

Рис. 5
Рис. 6

4. Оценка возможных эффектов. Рассмотренная модель расширяет количество возможных сценариев развития разрушения для структурированных материалов с масштабной чувствительностью трещиностойкости. В качестве примера рассмотрим плоскую модельную задачу о начале распространения трещин в окрестности кругового отверстия в упругой пластине при статическом одноосном растяжении (рис. 7). Эта задача, в частности, моделирует разрушение морского льда, содержащего капиллярные поры, в окрестности капилляров. Ограничимся при этом приближением таких размеров концевой зоны, при которых градиент внешних напряжений по ее длине можно считать постоянным.

Рис. 7

Напомним, что в обычном представлении сценария разрушения при росте трещин, выходящих на контур отверстия при растяжении, по мере удаления вершины трещины от края отверстия коэффициент интенсивности напряжений в ее концевых областях монотонно возрастает, вследствие чего при достижении им величины ${{K}_{I}} = {{K}_{{IC}}}$ = const наступает неустойчивая фаза и переход к динамическому разрушению.

Изменение величины ${{K}_{{IC}}}$ под воздействием градиента напряжений меняет сценарий разрушения. Приведем некоторые оценки, основываясь на полученных соотношениях.

Распределение упругих напряжений в области концентрации напряжений в окрестности отверстия в плоскости имеет вид [32]

(4.1)
${{\sigma }_{y}} = {{\sigma }_{\infty }}\left( {2 + \frac{{{{R}^{2}}}}{{{{x}^{2}}}} + 3\frac{{{{R}^{4}}}}{{{{x}^{4}}}}} \right),$$\frac{{d{{\sigma }_{y}}}}{{dx}} = \sigma {\text{'}} = - 2{{\sigma }_{\infty }}\left( {\frac{{{{R}^{2}}}}{{{{x}^{3}}}} + 6\frac{{{{R}^{4}}}}{{{{x}^{5}}}}} \right)$
где R – радиус отверстия.

При формировании концевой зоны зародышевого нарушения (малой трещины) на границе отверстия основную роль играет градиент напряжений, создаваемый отверстием. Отсюда и из (3.9) следует характер изменения эффективной трещиностойкости для трещины на начальном этапе удаления от отверстия в виде

(4.2)
${{K}_{{1C}}} \approx n{{\left( {\left| {\sigma {\text{'}}} \right|} \right)}^{{ - 0.5}}} = {n \mathord{\left/ {\vphantom {n {\sqrt {2{{\sigma }_{\infty }}\left( {\frac{{{{R}^{2}}}}{{{{x}^{3}}}} + 6\frac{{{{R}^{4}}}}{{{{x}^{5}}}}} \right)} }}} \right. \kern-0em} {\sqrt {2{{\sigma }_{\infty }}\left( {\frac{{{{R}^{2}}}}{{{{x}^{3}}}} + 6\frac{{{{R}^{4}}}}{{{{x}^{5}}}}} \right)} }}$

На рис. 8b показан характер предполагаемого изменения эффективной трещиностойкости материала (льда) под влиянием градиента напряжений в окрестности отверстия. На рис. 8a показан характер изменения KIC при удлинении трещины для такой схемы нагружения (по [29]). Эта зависимость для малых величин $\left( {\ell {\text{/}}R - 1} \right)$ может быть приближенно представлена в виде

(4.3)
${{K}_{I}} \approx {{\sigma }_{\infty }}\sqrt \ell \sqrt {\frac{{\ell {\text{/}}R - 1}}{{\ell {\text{/}}R - 0.71}}} $
Рис. 8

Учитывая (4.2) и (4.3) запишем в таком же приближении условие предельного равновесия

${{K}_{{IC}}} = {{K}_{1}}({{\sigma }_{\infty }};\ell {\text{/}}R)$
$\begin{gathered} \hfill n{{\left( {\left| {\sigma {\text{'}}} \right|} \right)}^{{ - 0.5}}} \approx {{\sigma }_{\infty }}\sqrt \ell \sqrt {\frac{{\ell {\text{/}}R - 1}}{{\ell {\text{/}}R - 0.71}}} \\ \hfill {n \mathord{\left/ {\vphantom {n {\sqrt {2{{\sigma }_{\infty }}\left( {\frac{{{{R}^{2}}}}{{{{\ell }^{3}}}} + 6\frac{{{{R}^{4}}}}{{{{\ell }^{5}}}}} \right)} }}} \right. \kern-0em} {\sqrt {2{{\sigma }_{\infty }}\left( {\frac{{{{R}^{2}}}}{{{{\ell }^{3}}}} + 6\frac{{{{R}^{4}}}}{{{{\ell }^{5}}}}} \right)} }} = {{\sigma }_{\infty }}\sqrt \ell \sqrt {\frac{{\ell {\text{/}}R - 1}}{{\ell {\text{/}}R - 0.71}}} \\ \end{gathered} $
(4.4)
${{\sigma }_{\infty }} = {{{{n}^{{2/3}}}} \mathord{\left/ {\vphantom {{{{n}^{{2/3}}}} {\sqrt[3]{{2\left( {\frac{{{{R}^{2}}}}{{{{\ell }^{2}}}} + 6\frac{{{{R}^{4}}}}{{{{\ell }^{4}}}}} \right)\left( {\frac{{(\ell {\text{/}}R - 1)}}{{(\ell {\text{/}}R - 0.71)}}} \right)}}}}} \right. \kern-0em} {\sqrt[3]{{2\left( {\frac{{{{R}^{2}}}}{{{{\ell }^{2}}}} + 6\frac{{{{R}^{4}}}}{{{{\ell }^{4}}}}} \right)\left( {\frac{{(\ell {\text{/}}R - 1)}}{{(\ell {\text{/}}R - 0.71)}}} \right)}}}}$

Это соотношение иллюстрирует характер развития трещины в непосредственной близости от контура отверстия. В описываемом сценарии разрушения по мере роста трещины происходит возрастание эффективной трещиностойкости на фоне роста коэффициента интенсивности напряжений по мере удаления от границы отверстия.

В результате характер соотношения предельного равновесия для малой трещины в окрестности отверстия (рис. 9а) становится немонотонным. Появление экстремума в области $\ell {\text{/}}L$ ∼ 1.15 означает возможность торможения и остановки трещины на начальном этапе после ее инициирования и начального неустойчивого развития, что отражается и на изменении предельных внешних напряжений ${{\sigma }_{\infty }}$ (рис. 9b) по уровню критического напряжения (4.2). Необходимо напомнить, что рассматриваемый эффект относится только к диапазону высоких градиентов внешних напряжений. По мере его уменьшения эффект ослабления локальной трещиностойкости падает, ${{K}_{{ICef}}} \to {{K}_{{IC}}} \approx $ const, и сценарий разрушения возвращается к схемам с постоянным уровнем трещиностойкости. В данном варианте задачи это означает выход на неустойчивый режим развития трещины при большом удалении ее вершины от отверстия. Оценим реальную границу возможного устойчивого роста микротрещин в образцах морского льда для рассматриваемой ситуации. По данным [13] при проведении крупномасштабных экспериментов локальная прочность льда при растяжении (изгибе) ${{\sigma }_{t}}\sim {{\sigma }_{\infty }} \approx 0.5$ МПа. Для этого уровня напряжений величина ${{\sigma }_{\infty }}{\text{/}}n \approx 1.47$. По (26) можно оценить в тех же допущениях, что он соответствует устойчивому подрастанию трещины до размера $\ell {\text{/}}R \approx 2.3$. Поскольку трещина такого размера значительно влияет на распределение напряжений у отверстия, приводимая оценка является качественной оценкой сверху. При обозначенной в [13] солености льда на уровне 4–7 промилле пористость льда, заполненная рассолом, по данным [31] составляет n ∼ 0.04–0.06, отношение полурасстояния между порами к радиусу пор для такой пористости $\ell {\text{*/}}R \approx 4.1{\kern 1pt} - {\kern 1pt} 5.0$. Следовательно, устойчивое подрастание трещин до размера $\ell {\text{/}}R \approx 2.33$ должно смениться неустойчивым ростом до слияния с соседними при $\ell {\text{*/}}R \approx $ 4.1–5.0.

Рис. 9

Приведенные оценки показывают и возможный сценарий разрушения пористой структуры морского льда при растяжении в направлении поперек оси параллельных капилляров с образованием зоны предразрушения. В связи с этим упомянем результаты крупномасштабных экспериментов в [13], в которых было отмечено частичное сохранение несущей способности льда в объеме немалой концевой области магистрального разрыва, что свидетельствует о развитии в ней структуры разрушения в устойчивом режиме. Большую роль в развитии такой несущей структуры может играть этап устойчивого подрастания трещин в окрестности элементов исходной структуры материала (пор, кристаллов) при возрастании нагрузки за счет увеличения ${{K}_{{1Cef}}}$ при снижении уровня градиента напряжений в локальных процессах разрушения. При этом в пористой среде в окрестностях растущих микротрещин возможно провоцирование роста трещин у соседних пор. Это приводит к образованию систем (эшелонов) трещин, образующих связанную систему нарушений. Механизм образования упорядоченных систем трещин частично рассмотрен в [33] для вариантов ситуаций, в которых увеличение сопротивления развитию отдельной трещины связано с падением уровня коэффициента интенсивности напряжений от внешних нагрузок при постоянстве уровня трещиностойкости. Тот же эффект может быть получен, как было показано выше, и при увеличении коэффициента интенсивности напряжений, если возрастает локальная трещиностойкость материала. Описанные явления увеличивают многообразие возможных сценариев разрушения. Дополнительные вариации также дает влияние исходной нарушенности контура отверстия (поры). В самом деле, если контур поры достаточно гладкий (для нашего примера $\ell {}_{0}/R < 1.2$рис. 9b), инициирование разрушения происходит на падающем участке кривой критического напряжения. Поэтому сначала, после страгивания зародышевого нарушения, ожидается неустойчивый рост трещин до размера, определяемого восходящей ветвью кривой на рис. 9. Затем рост трещин проходит в устойчивом режиме, описанном выше. Размер области устойчивого режима развития трещин и, тем самым, возможность образования регулярной структуры множественного разрушения, таким образом, зависит от механического состояния контура капилляров. В частности, при начально-гладком контуре отверстия (при $\ell {}_{0}/R \to 1$) возможно отсутствие этапа устойчивого роста. В этом случае можно прогнозировать развитие магистрального разрыва без образования области предразрушения.

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

Полагая в (4.4) $\ell {\text{/}}R \approx {{\nu }^{{ - 0.5}}}$ , где $\nu $ – объем рассола в капиллярных порах, и аппроксимируя возрастающую ветвь функции ${{\sigma }_{\infty }}(\ell {\text{/}}R)$ линейной функцией, получим

(4.5)
${{\sigma }_{\infty }} = {{\sigma }_{t}} \approx {{n}^{{2/3}}}(0.17 + 0.55{{\nu }^{{ - 0.5}}}) \approx 0.058(1 + 3.24{{\nu }^{{ - 0.5}}})$

На рис. 10 показан вид этой функции (при n ∼ 0.2 [(МПа)1.5]) в сравнении с функцией, предложенной в [34] (кривая 1), а также данными экспериментов (точки на рисунке) и обобщающей функцией по [35] (кривая 2). Заметим, что среди многочисленных работ по исследованию прочности морского льда при изгибе нет единого представления о виде определяющих соотношений, связывающих прочность и пористость. Найденные зависимости влияния объема рассола варьируются от линейных функций [36] до экспоненциальных [34]. Введение процессов, характеризующих множественное разрушение при подготовке разрыва может облегчить анализ сценария разрушения. Эффект немонотонного изменения эффективной трещиностойкости, когда при страгивании микротрещины в условиях растяжения происходит скачкообразное ее подрастание, сменяющееся на некотором участке устойчивым ростом, может также служить основой механизма развития регулярной структуры разрушения в условиях крупномасштабной концентрации растягивающих напряжений. Существование этапа возрастания эффективной трещиностойкости может приводить к устойчивому подрастанию многих очагов разрушения, образующих упорядоченную структуру разрушения при взаимодействии трещин, растущих в окрестности отдельных элементов исходной структуры. В частности, это относится к развитию структуры предразрушения в немалой концевой зоне в окрестностях вершины магистрального разрыва или сдвига. Такой механизм дополняет возможные сценарии микроразрушения в описанных условиях, обычно относимые к влиянию статистической неоднородности прочностных свойств материалов.

Рис. 10

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

Такие области в структурированных материалах содержат нарушения, предшествующие развитию магистральных разрывов. К ним относятся, в частности, концевые зоны магистральных разрывов, в которых материал содержит множественные нарушения (микротрещины, приуроченные к элементам микроструктуры), сохраняя при этом некоторую долю несущей способности. Аналогичные явления наблюдаются при развитии магистрального разрыва в крупномасштабных фрагментах морского ледяного покрова, имеющего столбчатозернистую структуру и капиллярную пористость [15]. Описанный процесс объясняет генезис структуры концевой зоны крупномасштабных разрывов в ледяном покрове.

Этот механизм также создает принципиальную возможность ранней диагностики повреждений в окрестностях концентраторов напряжений для материалов, обладающих аналогичной тенденцией изменения эффективной трещиностойкости, поскольку появление малых трещин в условиях высокого градиента напряжений (и пониженной локальной трещиностойкости), обозначая инициирование разрушения, не приводит к их распространению в неустойчивом режиме. Работа выполнена при финансовой поддержке программы Президиума РАН № I.2.55П: “Арктика – научные основы новых технологий освоения, сохранения и развития”.

Список литературы

  1. Лавров В.В. Деформация и прочность льда. Л.: Гидрометеоиздат, 1969. 206 с.

  2. Hoover C.G., Bažant Z.P. Cohesive crack, size effect, crack band and work – of-fracture models compared to comprehensive concrete fracture tests // International journal of fracture. 2014. V. 187(1). P. 133–143.

  3. Frigo B. Fantillia A.P., Chiaia B. Size effect on fracture toughness of snow. XXIV Italian Group of Fracture Conference, 1–3 March 2017, Urbino, Italy. Procedia Structural Integrity 3. 2017. P. 261–268.

  4. Тышко К.П., Черепанов Н.В., Федотов В.И. Кристаллическое строение морского ледяного покрова. СПб.: Гидрометеоиздат. 2000. 66 с.

  5. Гольдштейн Р.В., Ладыгин В.М., Осипенко Н.М. Модель разрушения слабо пористого материала при сжатии и растяжении // Физико-технические проблемы разработки полезных ископаемых. 1974. № 1. С. 3–16.

  6. Гольдштейн Р.В. Разрушение при сжатии // Успехи механики. 2003. Т. 2. № 2. С. 3–20.

  7. Goldstein R.V., Osipenko N.M. Some questions on ice and ice cover fracture in compression. In: Ice-structure interaction. IUTAM-IAHR Symp. St. John’s, New Foundland, Canada, Berlin, Heidelberg: Springer-Verlag, 1991. P. 251–266.

  8. Lajtai E.Z., Lajtai V.N. The collapse of cavities // Int. J. Rock Mech. Min. Sci. & Geomech. Abstr. 1975. V. 12. P. 81–86.

  9. Гольдштейн Р.В., Осипенко Н.М. Некоторые вопросы механики прочности морского льда // Физическая мезомеханика. 2014. Т. 17. № 6. С. 59–69.

  10. Timko G.W., Trederking R.M.W. Compressive strength of ice sheets // Cold Regions Science and Technology. 1990. V. 17. P. 227–240.

  11. Schwarz J., Weeks W.F. Engineering properties of sea ice // Journal of Glaciology. 1977. V. 19. № 8I. P. 499–531.

  12. Гольдштейн Р.В., Осипенко Н.М. Механика разрушения и некоторые вопросы разрушения льда. В сб.: Механика и физика льда. М.: Наука, 1983. С. 65–94.

  13. Dempsey J.P., Adamson R.M., Mulmule S.V. Scale effects on the in-situ tensile strength and fracture of ice. Part II: First-year sea ice at Resolute, NWT // International Journal of Fracture. 1999. V. 95(1). P. 347–366.

  14. Dempsey J.P., DeFranco S.J., Adamson R.M., Mulmule S.V. Scale effects on the in-situ tensile strength and fracture of ice. Part I: Large grained freshwater ice at Spray Lakes Reservoir, Alberta // International Journal of Fracture. 1999. V. 95(1) P. 325–345.

  15. Dempsey J., Cole D.M., Wang S. Tensile fracture of a single crack in first-year sea ice // Philosophical Transactions of the Royal Society A: Mathematical, Physical and Engineering Sciences. V. 376(2129):20170346, September 2018. https://doi.org/10.1098/rsta.2017.0346

  16. Dempsey J.P. Research trends in ice mechanics // International Journal of Solids and Structures. 2000. V. 37. P. 131–153.

  17. Totman C.A., Uzorka O.E., Dempsey J.P., Cole D.M. Sub-size fracture testing of FY sea ice. Proceedings of the 6th International Conference on Fracture Mechanics of Concrete and Concrete Structures. V. 3. P. 1683–1690.

  18. Dempsey J.P., Xie Y., Adamson R.M., Farmer D.M. Fracture of a ridged multi-year Arctic sea ice floe // Cold Regions Science and Technology. 2012. V. 76–77. P. 63–68.

  19. Черепанов Г.П. Механика хрупкого разрушения. М.: Наука, 1974. 640 с.

  20. Вайншельбаум В.М., Гольдштейн Р.В. О материальном масштабе длины как мере трещиностойкости пластичных материалов и его роль в механике разрушения. Препринт. М.: ИПМ АН СССР, 1976. 70 с.

  21. Bazant Z.P., Planas J. Fracture and size effect in concrete and other quasibrittle materials. Boca Raton, U.S.A.: CRC Press, 1998. 616 p.

  22. Bažant Z.P. Size effect in blunt fracture: Concrete, rock, metal // ASCE Journal of Engineering Mechanics. 1984. V. 110. P. 518–535.

  23. Duan K., Hu X.-Z., Wittmann F.H. Size effect on fracture resistance and fracture energy of Concrete // Materials and Structures/Materiaux et Constructions. V. 36. March 2003. P. 74–80.

  24. Сукнев С.В. Нелокальные и градиентные критерии разрушения квазихрупких материалов при сжатии // Физическая мезомеханика. 2018. Т. 2. № 4. С. 22–32.

  25. Lajtai E.Z. Effect of tensile stress gradient on brittle fracture initiation. // Int. J. Rock Mech. Min. Sci. 1972. V. 9. № 5. P. 569–578.

  26. Леган М.А. О взаимосвязи градиентных критериев локальной прочности в зоне концентрации напряжений с линейной механикой разрушения // ПМТФ. 1993 Т. 34. № 4. С. 146–154.

  27. Новопашин М.Д., Сукнев С.В. Градиентные критерии предельного состояния // Вестник СамГУ, Естественнонаучная серия. 2007. № 4(54). С. 316–334.

  28. Сукнев С.В. Нелокальные и градиентные критерии разрушения квазихрупких материалов при сжатии // Физическая мезомеханика. 2018. Т. 21. № 4. С. 22–32.

  29. Справочник по коэффициентам интенсивности напряжений / Ред. Ю. Мураками. М.: Мир, 1990. Т. 1. 448 с.

  30. Bažant Z.P., Kazemi M.T., Hasegawa T., Mazars J. Size effect in Brazilian split-cylinder tests: Measurments and fracture analysis // ACI Materials Journal. 1991. V. 88. P. 325–332.

  31. Богородский В.В., Гаврило В.П. Лед. Л.: Гидрометеоиздат, 1980. 384 с.

  32. Тимошенко С.П., Гудьер Дж. Теория упругости. М.: Наука, 1975. 576 с.

  33. Гольдштейн Р.В., Осипенко Н.М. О модели разрушения структурированной среды в условиях сжатия // Изв. РАН. МТТ. 2010. № 6. С. 86–97.

  34. Timco G.W., O`Brien S. Flexural strength equation for sea ice // Cold Regions Science and Technology. 1994. V. 22. P. 285–298.

  35. Пэн. Ч. Экспериментальные исследования и анализ характеристик предела прочности морского льда на изгиб вдоль побережья Бохайского залива // Молодой ученый. 2017. № 27. С. 39–46.

  36. Dykins I. Tensile and flexural properties of saline ice. Proc. Int. Symp. on Physics of Ice. Munich. 1968. P. 9–14.

Дополнительные материалы отсутствуют.