Full Text

Введение Большое количество искусственных и естественных материалов обладают, помимо композиционного состава, многоуровневой поровой структурой. В зависимости от условий функционирования материала его поры могут быть пустыми или содержать поровую жидкость (например, воду). Широко распространенными и хорошо известными представителями класса пористых водонасыщенных композитов являются горные породы, бетоны, костные ткани. Известно, что присутствующая в порах этих материалов вода или водный раствор способны оказывать значительное влияние на их механические свойства. В частности, поровая жидкость полагается ответственной за многочисленные явления преждевременного разрушения вмещающего материала при нагрузках, существенно меньших, чем предельно допустимые для этих же материалов, не насыщенных жидкостью [1-5]. Влияние поровой воды на механические характеристики (в частности, снижение прочности) пористых водонасыщенных композитов проявляется по-разному на различных временных масштабах. Так, на длительных временных промежутках вода приводит к развитию различных дефектов в твердофазном каркасе и его физико-химической деградации, известной как адсорбционное понижение прочности. Данный эффект определяет, в частности, снижение статической прочности цементного камня и бетона [3, 6, 7]. Другим важным фактором, определяющим механические свойства бетона в условиях динамического изменения напряженного состояния (т.е. на малых временных интервалах, когда скорость роста нагрузки превосходит скорость перераспределения флюида в твердофазном каркасе), является механическое влияние воды. Оно заключается в распирающем действии порового давления на твердофазный каркас и проявляется в снижении макроскопической прочности и повышении эффективного модуля всестороннего сжатия материала (а также других упругих модулей, связанных с ним) [8-10]. При этом, ввиду высокой подвижности воды, ее механическое влияние на напряженное состояние и механические характеристики пористого материала проявляется практически во всем интервале скоростей деформирования, в том числе в интервале скоростей, при которых материал каркаса является нечувствительным к динамике нагружения [9, 10]. Действительно, при механическом нагружении свободная (химически не связанная) вода перераспределяется в системе связанных несплошностей материала, что является важным каналом перераспределения и релаксации локальных напряжений. Данное влияние может быть значимым уже при малом объемном содержании воды в материале, составляющем несколько процентов. Механическое влияние поровой жидкости на эффективные динамические механические свойства наиболее критично для пористых материалов с хрупким каркасом. Величина вклада поровой жидкости в напряженное состояние и прочность хрупкого каркаса в значительной степени определяется, помимо объемного содержания жидкости и динамики нагружения, и параметрами поровой структуры (распределением пор по размерам, характерными масштабами пор, их связанностью и т.д.). Эта проблема актуальна не только для горных пород [11], но и для искусственных композиционных конструкционных материалов, включая тяжелые высокопрочные бетоны. Тяжелые бетоны традиционно применяются в качестве строительного материала для наиболее ответственных элементов дамб, плотин, опор различных мостов, вышек, портовых сооружений. Тяжелые бетоны характеризуются наличием двухмасштабной поровой структуры. Поскольку влияние поровой жидкости на величину напряжений и прочность каркаса во многом определяется возможностью ее перераспределения в поровом пространстве, каждая из поровых подсистем в бетоне вносит свой вклад в интегральный механический эффект жидкости. Отметим, что при этом важную роль играет степень связности поровых систем различных масштабов. Бетоны характеризуются кратно более высокими значениями прочности на сжатие в сравнении с прочностью на растяжение и традиционно предназначаются для функционирования в условиях действия сжимающих нагрузок. «Распирающее» влияние порового флюида на напряженное состояние каркаса также наиболее значительно в условиях сжатия. В связи со сказанным актуальной задачей является исследование и обобщение влияния порового флюида на величину динамической прочности на сжатие высокопрочных бетонов, в том числе выявление вклада каждой из поровых подсистем и их взаимной связанности. Результаты такого исследования имеют не только фундаментальную значимость, но и широкий спектр практических приложений. В частности, они могут быть востребованы при решении задач определения «оптимальных» параметров внутренней структуры бетонов, функционирующих в условиях влагонасыщения или в водной среде. Другим перспективным приложением, которое заслуживает отдельного упоминания, является оценка механического влияния воды на деградацию поверхностных слоев и износ поверхностных слоев бетонных элементов конструкций в крупных водоемах северных широт, в том числе в условиях арктического климата [12]. Такие опорные бетонные элементы находятся в условиях длительного контактного нагружения протяженными (десятки километров) ледовыми плитами толщиной до 1-2 м. При динамических подвижках ледовых массивов со скоростями до нескольких метров в секунду (например, при «льдотрясениях») локальные давления в пятнах контакта могут кратковременно достигать десятков мегапаскалей и даже превышать величину прочности бетона на одноосное сжатие. В условиях динамического нагружения поровое давление в поверхностных областях бетона в окрестностях пятен контакта многократно возрастает вследствие относительно низкой скорости фильтрации в поровых каналах и может провоцировать локальное разрушение. Определение условий локального разрушения при динамическом контактном (сжимающем) нагружении и оценка соответствующих значений динамической прочности актуальны для разработки и параметризации статистических моделей абразии бетонных сооружений массивными плавучими объектами (ледовыми массивами или массивными судами) [13-15]. Традиционными способами изучения динамического механического поведения водонасыщенных материалов в условиях сжатия являются лабораторные эксперименты [16-19] и мониторинг крупномасштабных объектов, функционирующих в водонасыщенном состоянии [15, 20-22]. Перспективной альтернативой является компьютерное моделирование динамического нагружения представительных образцов бетонов с использованием связанных механических моделей, учитывающих деформацию каркаса, перераспределение (фильтрацию) порового флюида и взаимовлияние этих процессов. Очевидным преимуществом такого подхода является возможность не только качественного, но количественного определения вкладов различных факторов поровой структуры в механическое поведение (включая разрушение) и динамическую прочность водонасыщенного бетона. В настоящей работе решение поставленной задачи осуществлялось путем компьютерного моделирования одноосного сжатия флюидонасыщенных представительных образцов бетона с двухмасштабной поровой структурой гибридным методом однородно деформируемых проницаемых дискретных элементов [23, 24]. Применялась связанная механическая модель водонасыщенных хрупких материалов, основанная на принципах модели пороупругости Био. Отметим, что применяемая модель учитывает механическое взаимовлияние твердофазного каркаса и порового флюида (через деформацию порового пространства и соответствующее изменение порового давления), а также фильтрацию флюида в поровом пространстве [25]. Решение задачи осуществлялось в два этапа. На первом этапе исследовалась зависимость величины динамической прочности водонасыщенного бетона на сжатие от параметров, определяющих мобильность флюида в поровом пространстве, и анализировалась возможность построения обобщенной зависимости в терминах комбинированного безразмерного параметра. Далее исследовались вклады микропоровой и капиллярной составляющих порового пространства (а также их интегрированности в единую проницаемую сеть) в величину динамической прочности бетонов. 1. Объект исследования Объектом изучения являлись мезомасштабные образцы тяжелого высокопрочного бетона, поровое пространство которых полагалось насыщенным жидкостью (в качестве «базовой» жидкости рассматривалась вода). Исходное значение порового давления жидкости полагалось равным атмосферному давлению. Под мезомасштабом здесь подразумевается характерный масштаб ключевых элементов внутренней структуры: агрегатов миллиметровых линейных размеров и «микропор» с характерными размерами сотни микрометров. На мезоскопическом уровне описания структура бетона представляет собой матрицу цементного камня (например, на основе высокопрочного бездобавочного портландцемента ПЦ 600 ДО, ПЦ 500 ДО) с упрочняющими крупными агрегатами заполнителя (щебень базальта и диабаза). Цементный камень является основным компонентом высокопрочных бетонов, его объемное содержание варьируется в интервале 30-60 %. В качестве заполнителя в работе рассматривался щебень базальта, который полагался однородным изотропным беспористым (непроницаемым) материалом. Типичная форма включений заполнителя (частиц базальтового щебня) - равноосная, средний размер включений - 10-15 мм. Одним из ключевых структурных параметров цементного камня, определяющих его эффективные прочностные и упругие свойства, является пористость. В зависимости от технологических режимов приготовления цементного камня и свойств компонентов его пористость может быть различной как в качественном, так и в количественном отношении. В качественном отношении различие может наблюдаться, в частности, по типу поровой структуры: например, одно- и двухмасштабная, сообщающаяся и изолированная пористость. В количественном отношении различие может касаться, в частности, величины объема порового пространства (приходящегося на пористость определенного масштабного уровня) и проницаемости (диаметра поровых каналов). В соответствии с литературными данными тяжелый бетон характеризуется двухуровневой поровой структурой, капиллярная пористость достигает ~6 % объема цементного камня, содержание микропор - до 4 % объема [26-28]. Характерный размер микропор составляет сотни микрометров, их форма, как правило, является округлой. Поперечные размеры капиллярных каналов-пор могут варьироваться в широких пределах (от 0,05 до 15 мкм), что определяется различными факторами, в частности, условиями приготовления и твердения бетона, а также соотношением воды и цемента. Изменение размеров этих пор является следствием отложения продуктов гидратации цемента в каналах-порах при затвердевании бетона, которое может приводить к многократному сужению их характерного просвета и соответствующему снижению проницаемости цемента (до 1-2 порядков величины) [26, 29, 30]. Просветность капилляров определяет мощность перераспределения поровой жидкости, обусловленного градиентом порового давления. Действительно, перераспределение жидкости в поровом пространстве материала способствует уменьшению максимальных локальных значений «распирающего» порового давления на каркас. В случае снижения характерного диаметра каналов уменьшается плотность мощности потока жидкости, что определяет более высокие значения порового давления. Другими словами, при одной и той же величине скорости внешнего механического нагружения в материале с меньшим характерным диаметром каналов поровое давление возрастает быстрее, и, соответственно, разрушение начинается при меньшем значении приложенной нагрузки. Сказанное определяет актуальность анализа влияния поперечного размера капилляров на динамическое поведение и величину динамической прочности водонасыщенного бетона. Отметим, что другими существенными факторами являются объемное содержание микропор в цементном камне и их связанность с капиллярной поровой системой. В частности, при интенсивном отложении продуктов гидратации цемента на стенках микропор может происходить «закупоривание» выходов капиллярных каналов, т.е. изоляция микропорового пространства. В соответствии со сказанным выше в работе исследовалось динамическое поведение модельных образцов бетона в условиях одноосного сжатия при различных значениях капиллярной проницаемости, содержания микропор и степени их «закупоренности». Рассмотрены несколько характерных случаев поровой структуры цементного камня. В первом случае его поровая сеть содержит только капиллярные поры. Во втором - поровое пространство содержит, помимо капиллярной сети, и микропоры, связанные с капиллярной пористостью в единую фильтрационную систему. В третьем случае микропоры полагаются «закупоренными», т.е. не сообщающимися с сетью капилляров (данный случай соответствует интенсивному отложению продуктов гидратации цемента на стенках микропор и «закупориванию» выходов капиллярных каналов). Во всех случаях поровое пространство цементного камня было заполнено водой (в том числе «закупоренные» микропоры). 2. Связанная механическая модель высокопрочного бетона с двухмасштабной поровой структурой На основе описанных в предшествующем разделе положений построена двумерная компьютерная модель представительного мезомасштабного объема высокопрочного бетона. Модель включает в рассмотрение два пространственных масштаба: мезоскопический и микроскопический. Мезомасштабные компоненты внутренней структуры бетона (цементный камень, агрегаты, микропоры) задаются явным образом (рис. 1). Пространственное размещение и определение геометрии агрегатов в модельных образцах осуществлялось с использованием реализованного исполнителями алгоритма диаграмм Воронова. Микропоры случайным образом распределялись в объеме матрицы цементного камня. Задание микропор цементного камня осуществляется удалением определенных областей материала из начальной беспористой структуры образца. При создании модели учитывались такие параметры мезомасштабной структуры, как объемное содержание агрегатов и микропор, их средний размер, форма и пространственное распределение в образце. В проведенных исследованиях объемное содержание базальтовых агрегатов принималось равным 10 %. Агрегаты полагались приблизительно равноосными, угловатыми. Объемное содержание микропор в цементном камне составляло 0 % (только капиллярная пористость) или 4 %. Микропористость цементного камня считалась моноразмерной, форма микропор - круглой (диаметр микропоры равен 500 мкм). Особенности внутренней структуры бетона на более низких масштабах (включая капиллярную пористость) учитываются эффективно (интегрально) через физико-механические свойства цементного камня, агрегатов и порового флюида. Эти свойства включают в себя применяемые определяющие соотношения и уравнения динамики для компонентов системы, а также конкретные значения коэффициентов и параметров этих соотношений. В качестве вычислительной основы для построения компьютерной модели использовался развиваемый авторами метод однородно деформируемых дискретных элементов [21]. Метод дискретных элементов базируется на представлении материала как ансамбля химически связанных между собой частиц, имеющих заданную массу, форму, объем и поверхность. Три последних могут изменяться в процессе взаимодействия элементов. Каждый элемент моделирует некоторую область (фрагмент) цементного камня или базальтового агрегата. Динамика моделируемого тела определяется численным решением системы уравнений движения Ньютона - Эйлера, каждое из которых описывает движение и поворот одной из частиц. Авторами используется наиболее распространенная реализация метода, в рамках которой уравнение движения частицы формулируется в упрощенной форме, отвечающей диску (в двумерной постановке задачи) или шару (в 3D случае) [31, 32]. Такое приближение имеет название приближения эквивалентных дисков/шаров. Основным преимуществом метода однородно деформируемых элементов в сравнении с классическим формализмом жестких (недеформируемых) элементов является учет изменения (как упругого, так и необратимого) формы и объема дискретного элемента в ходе нагружения. В рамках метода однородно деформируемых дискретных элементов напряжения и деформации в объеме элемента полагаются однородно распределенными и характеризуются тензорами усредненных напряжений и деформаций. Компоненты этих тензоров вычисляются как суперпозиция сил взаимодействия элемента с соседями и деформаций элемента в этих парах. Соотношения для сил взаимодействия элементов записываются в многочастичной форме и, в свою очередь, определяются значениями локальных напряжений. Такая взаимосвязь напряжений и сил взаимодействия позволяет естественным образом реализовывать различные модели упругости, пластичности и разрушения (в том числе с многопараметрическими критериями) в рамках формализма дискретных элементов. Микропоры Агрегаты Рис. 1. Пример структуры представительного мезомасштабного объема модельного бетона (содержание базальтовых агрегатов - 10 %, величина микропористости - 4 %) Fig. 1. The example of the structure of the representative mesoscale volume of the model concrete (the content of basalt aggregates - is 10 %, microporosity - is 4 %) Отметим, что развитая оригинальная реализация метода дискретных элементов (МДЭ) получила специальное название - метод подвижных клеточных автоматов (MCA). Данное преимущество определяет, в частности, возможность корректного моделирования механического отклика проницаемых флюидонасыщенных материалов с применением связанных моделей, эффективно учитывающих механическое взаимовлияние твердофазного каркаса и поровой жидкости, а также перераспределение (фильтрацию) жидкости в поровом пространстве материала. Для моделирования флюидонасыщенных материалов и сред методом однородно деформируемых элементов развит гибридный формализм, в рамках которого дискретный элемент рассматривается как пористый и проницаемый [24, 33]. Поровое пространство элемента описывается неявно и характеризуется интегральными параметрами, такими как объем порового пространства, пористость и проницаемость. Жидкость в поровом пространстве элемента учитывается неявно и также характеризуется интегральными параметрами: массой, плотностью, давлением и динамической вязкостью (здесь и далее полагаем поровую жидкость ньютоновской). В работе материал агрегатов (базальт) полагался изотропным упруго-хрупким, цементный камень рассматривался как изотропный упругопластический (величина неупругой деформации мала и сравнима с величиной упругой деформации). В качестве определяющего соотношения для данных компонентов в области упругой деформации использовался обобщенный закон Гука для изотропных материалов. Отметим при этом, что неупругое поведение мезомасштабных объемов цементного камня при нагружении обусловлено, главным образом, возникновением, накоплением и ростом повреждений микро- и наноскопического масштабов. Поэтому для описания релаксации напряжений и накопления неупругой деформации в этом компоненте бетона использовался неассоциированный закон пластического течения с критерием Мизеса - Шлейхера (модель Николаевского) [33], широко применяемый для хрупких природных и искусственных материалов [35, 36]. В рамках применяемой связанной модели описание механического взаимодействия твердофазного каркаса проницаемого дискретного элемента, моделирующего фрагмент цементного камня, и жидкости в его капиллярном поровом пространстве осуществлялось в рамках модели пороупругости Био. Поровая жидкость полагалась линейно сжимаемой. Основные определяющие соотношения, связывающие усредненные напряжения и деформации твердофазного каркаса элемента, а также давление жидкости в поровом (капиллярном) пространстве элемента имеют вид: (1) (2) (3) Здесь символ D обозначает приращение соответствующей переменной за шаг по времени численной схемы интегрирования уравнений движения элемента; sab и eab - компоненты тензоров усредненных напряжений и деформаций в объеме дискретного элемента; σmean - среднее напряжение; G и K - модули сдвига и всестороннего сжатия материала элемента в отсутствие поровой жидкости; dab - символы Кронекера, а = (1-K/Ks) - коэффициент пороупругости в модели Био (Ks - модуль всестороннего сжатия материала стенок каркаса); Vp и - текущее и исходное (в недеформированном элементе) значения порового объема; j и j0 - соответствующие значения пористости; V - объем элемента; и - равновесные значения плотности и давления жидкости при атмосферных условиях (в отсутствие механического ограничения объема жидкости); rfl - текущее значение плотности жидкости в поровом пространстве элемента; Kfl - модуль объемного сжатия жидкости. Влияние порового давления в объеме дискретного элемента, моделирующего цементный камень, на его неупругое механическое поведение учитывалось в рамках использования модифицированной формулировки критерия пластичности Мизеса - Шлейхера с использованием эффективного среднего напряжения Терцаги: (4) где α - коэффициент внутреннего трения; Y - сцепление (предел упругости материала при чистом сдвиге); - интенсивность напряжений. В рамках метода дискретных элементов локальное разрушение (на рассматриваемом масштабе) моделируется разрывом химической связи между элементами, на поверхности взаимодействия которых выполняется заданное критическое условие. В настоящем исследовании в качестве критерия разрушения цементного камня и базальтовых агрегатов применялся критерий Друккера-Прагера [37], также формулируемый в терминах эффективных напряжений Терцаги: (5) где β = σc/σt, σc и σt - значения прочности «сухого» материала на сжатие и растяжение. Для беспористого базальта применялась традиционная формулировка критерия Друккера - Прагера, в которой вместо используется σmean. Отметим, что применяемые критерии пластичности (4) и прочности (5) являются обобщениями классического критерия пластичности/разрушения Кулона - Мора. Их выбор обусловлен тем, что они позволяют адекватно описывать отклик хрупких материалов (в том числе флюидонасыщенных) в сложном напряженном состоянии. При описании перераспределения жидкостей в связанном поровом пространстве бетонов традиционно учитывают следующие механизмы: 1) фильтрация, обусловленная градиентом (перепадом) давления поровой жидкости; 2) просачивание вследствие капиллярных сил; 3) диффузия жидкости через наноразмерные несплошности, которые в некотором количестве всегда присутствуют в цементном камне и могут быть интегрированы в капиллярную сеть. В данной работе исследовался механический отклик бетонов при динамическом нагружении, поэтому вкладом диффузии жидкости через несплошности нанометрового диапазона пренебрегали. Мощность капиллярного механизма перераспределения также достаточно мала, но в то же время капиллярные силы могут обусловливать «запирание» фильтрации (фильтрация начинается, когда величина перепада давления превышает пороговое значение, которое для бетонов может варьироваться от атмосферного давления до 1 МПа и выше). Величина порогового давления обратно пропорциональна характерному поперечному размеру капиллярных каналов. В то же время коэффициенты этой связи крайне чувствительны как к составу цементного камня, так и к составу поровой жидкости. Кроме того, характерные значения порового давления при сжимающих нагрузках, сопоставимых с прочностью бетона, многократно превышают 1 МПа. Поэтому без существенного снижения общности результатов в настоящей работе в качестве движущей силы перераспределения жидкости в капиллярном поровом пространстве элементов рассматривался только градиент порового давления (вклады капиллярного течения и капиллярных сил не учитывались). Фильтрация описывалась на основе уравнения транспорта жидкости в форме, учитывающей конечную скорость распространения возмущений [38]: , где . (6) Здесь - плотность мощности потока жидкости; k - капиллярная проницаемость; h - вязкость флюида; tr - время релаксации; - скорость звука в жидкости, Уравнение (6) решалось методом конечных объемов на ансамбле дискретных элементов. Детальное описание численной реализации модели флюидонасыщенных материалов в рамках гибридного МДЭ приведено в работах [24, 33]. В проводимом исследовании использовались следующие значения механических параметров компонентов модельного бетона. Значения упругих констант базальта и цементного камня и параметры единой кривой упрочнения (для цементного камня) соответствовали данным, приведенным в [39]. Для базальта они составляют G = 60 ГПа и K = 166 ГПa, для цементного камня G = 17 ГПа и K = 62,5 ГПa. Как отмечалось выше, цементный камень моделировался как упругопластический материал. Величина сцепления Y принималась равной 18,5 МПа, величина коэффициента деформационного упрочнения - 14,1 ГПа. Для базальта использовались следующие значения параметров прочности в критерии разрушения (5): 260 МПа и 52 МПа, для цементного камня: 120 МПа и 30 МПа. Величина капиллярной проницаемости цементного камня k полагалась непосредственно связанной с капиллярной пористостью j (в настоящем исследовании j = 6 %) и характерным поперечным размером (эффективным диаметром) капиллярных поровых каналов dch известным соотношением k = jdch2 [40]. В проводимом исследовании величина эффективного диаметра капиллярных каналов варьировалась в различных расчетах в пределах от 0,1 до 1 мкм для анализа влияния степени отложения продуктов гидратации цемента в капиллярном поровом пространстве на динамическую прочность бетона. Жидкость в микропорах цементного камня описывалась эффективно, аналогично учету жидкости в капиллярном поровом пространстве элемента. Каждая микропора характеризовалась массой, плотностью и давлением находящейся в ней жидкости (использовалось приближение равнораспределения плотности жидкости в микропоре), текущий объем микропоры определялся на основе специального алгоритма, описанного в работе [24]. Начальное значение порового давления воды в капиллярных и микропорах (до начала нагружения) было равно величине атмосферного давления. Интегрированность микропор в капиллярную поровую сеть (т.е. возможность истечения флюида из микропор в капиллярную сеть или притока флюида из капиллярной сети) моделировалась параметрическим заданием величины эффективного диаметра капиллярных поровых каналов dch на стенках микропор. Рассматривались два предельных случая: величина dch равна значению того же параметра в объеме цементного камня (полная связанность порового пространства) и dch = 0 («закупоренные» микропоры). Межфазные границы описывались с использованием модели «бесконечно тонкой» межфазной границы [41] с использованием непотенциальной модели Жубеля - Бейлора [42, 43]. Основной характеристикой интерфейса в рамках данной модели является критерий прочности. Для границ раздела включений базальта и цементного камня использовался критерий прочности Друккера - Прагера, коэффициенты которого отвечали коэффициентам критерия разрушения для цементного камня. 3. Результаты моделирования и обсуждение 3.1. Условия нагружения и варьируемые параметры модели В работе проводились детальные параметрические исследования влияния свободной поровой воды на величину динамической прочности и разрушение мезомасштабных образцов высокопрочного бетона с одно- и двухмасштабной поровой структурой при одноосном сжатии. Двумерные бетонные образцы моделировались в приближении плосконапряженного состояния. В некотором приближении (в частности, в отношении жесткости образца) это соответствует условиям реального эксперимента на бетонных кубических образцах. Наличие естественной смазки на торцевых поверхностях образцов (прослойки жидкой воды наноразмерной толщины) учитывалось эффективно [26-28]. В модели это обеспечивалось путем разрешения поперечных смещений элементам пуансонов, контактирующих с образцом, и использования вязкой тангенциальной силы между контактирующими элементами пуансонов и поверхности образца (коэффициент сухого трения полагался равным нулю). Боковые поверхности образца являлись открытыми и обеспечивали возможность свободного истекания флюида из образца (окружающее образец пространство полагалось заполненным воздухом при атмосферном давлении). Моделировалось одноосное сжатие образцов с постоянной скоростью. При нагружении нижний (опорный) пуансон полагался неподвижным, верхний пуансон смещался в вертикальном направлении. Для плавного выхода образца на установившийся режим нагружения скорость верхнего (сжимающего) пуансона сначала возрастала по линейному закону до достижения заданного значения, а затем оставалась постоянной до разрушения образца. Интервал рассматриваемых скоростей нагружения определялся на основе оценки локальных скоростей деформации поверхностей бетонных сооружений при контакте с морскими льдами (или судами). Эти значения могут различаться на 3-6 порядков величины (от квазистатического до 100-102 с-1). Интервал значений скорости смещения верхнего пуансона V определялся на основе данных оценок скоростей деформации. При проведении исследования, помимо скорости деформации , варьировались также пропускная способность капиллярной поровой сети (ее проницаемость k), поперечный размер (ширина) образца и динамическая вязкость поровой жидкости h (в качестве которой потенциально могут выступать пресная или соленая вода, а также различные растворы). В расчетах скорость деформации (V - скорость сжатия образца, H - высота образца) варьировалась в пределах трех порядков величины (от 5×10-3 до 5 с-1), вязкость жидкости - в пределах порядка величины относительно вязкости соленой воды (2×10-4 до 3×10-3 Па×с), проницаемость - в пределах двух порядков величины (от 6×10-16 до 6×10-14 м2), ширина образца - в 4 раза (от 5 до 20 см). Анализировались величина прочности на сжатие и характер разрушения образцов. 3.2. Динамическая прочность бетонов с одномасштабной (капиллярной) поровой структурой Из результатов многочисленных теоретических и экспериментальных исследований хорошо известен нелинейный характер зависимости прочности на сжатие образцов хрупких флюидонасыщенных материалов от скорости деформирования [5, 11, 24, 33]. В квазистатическом режиме нагружения величина прочности стремится к верхнему пределу, отвечающему прочности ненасыщенного образца. По мере роста скорости сжатия величина прочности стремится к нижнему пределу, отвечающему прочности водонасыщенных гидроизолированных образцов. При этом размах изменения прочности может достигать нескольких десятков процентов. Качественно такой тренд описывается на основе рассмотрения конкуренции двух процессов. Первым из них является сжатие порового пространства материала, вызывающее рост порового давления. Данный эффект обусловливает более раннее (в сравнении с «сухими образцами») начало зарождения повреждений и развитие трещин в образцах и тем самым определяет пониженное значение прочности. Вторым, противоположно направленным процессом является фильтрация порового флюида в образце при сжатии и его истечение через свободные поверхности вследствие возникновения и роста градиента давления от центра образца к периферии. Соотношение этих конкурирующих процессов определяет скорость роста порового давления в материале, а следовательно, и конкретное значение прочности при заданной скорости нагружения. В исследованиях авторов, проведенных ранее на макроскопических образцах различных флюидонасыщенных хрупких материалов с одномасштабной пористостью, было показано [24, 33], что величина прочности является функцией не просто от скорости деформирования, но отношения скорости деформирования к характерной скорости перераспределения порового флюида в образце. Последняя определяется комбинацией материальных и геометрических параметров, включающей проницаемость твердофазного каркаса k, вязкость жидкости h, характерное расстояние фильтрации W (ширина образца) и характерную величину разности поровых давлений DP на расстоянии W (т.е. величину перепада давления на полуширине образца). Рассматриваемые в настоящей работе бетоны обладают двухмасштабной пористостью, при этом характерная скорость фильтрации определяется, главным образом, проницаемостью сети капиллярных пор k. Для изучения влияния проницаемости цементного камня, обусловленной связанной сетью капиллярных пор-каналов, на прочность бетона проведено моделирование одноосного сжатия модельных образцов с цементным камнем, содержащим только капиллярную пористость (микропоры отсутствуют). Типичная начальная структура одного из таких образцов соответствует представленной на рис. 1 только при нулевом значении микропористости. При проведении исследования варьировались скорость деформации , вязкость жидкости h, пропускная способность капиллярной поровой сети (ее проницаемость k), ширина образца W. Результаты моделирования показали, что прочность мезомасштабных образцов водонасыщенного бетона является однозначной функцией отношения характерных скоростей деформации и фильтрации. Установлено, что в качестве безразмерного параметра, характеризующего это отношение, может быть использован комбинированный параметр, формулируемый аналогично числу Дарси [11]: (7) Параметр (7) имеет тот же вид, что и число Дарси, однако между ними есть важное различие. В определении числа Дарси фигурирует актуальное (текущее) значение перепада порового давления DP (в рассматриваемом случае это разность давлений в центре образца и на свободной боковой поверхности в поперечном к оси нагружения направлении). В процессе деформирования образцов величина числа Дарси непрерывно возрастает (поровое давление в центре образца постоянно увеличивается, а на боковой поверхности остается постоянным и равным атмосферному). Поэтому в соотношении (7) в качестве DP используется константа, определяемая как максимальное значение перепада поровых давлений по Терцаги. Она равна разности оценки Терцаги максимально достижимого порового давления в цементном камне ( где j0 - исходная капиллярная пористость цементного камня) и давления на боковых границах образца (в данном случае 105 Па): Для моделируемого цементного камня DP ~ 7 МПа. Полученная зависимость прочности бетонных образцов от аналога числа Дарси хорошо аппроксимируется сигмоидальной кривой: (8) где - динамическая прочность бетонного образца на сжатие (в МПа); » 55,6 МПа - минимальное значение прочности бетонного образца, достигаемое в условиях пренебрежимо малой фильтрации (при высоких скоростях нагружения или малой проницаемости), » 64 МПа - максимальное значение прочности, отвечающее прочности сухого образца или водонасыщенного образца, деформируемого квазистатически (когда вытекание жидкости через боковые стенки поддерживает величину порового давления, равную атмосферному). Характерный размах изменения величины прочности составляет ( ) » 8 МПа или 13 % от величины прочности ненасыщенного бетона. Соответствующая зависимость приведена на рис. 2 (кривая 1). Можно видеть, что она характеризуется наличием трех характерных участков. Рис. 2. Зависимость прочности на сжатие водонасыщенных бетонных образцов с различной поровой структурой от аналога безразмерного числа Дарси: кривая 1 - одноуровневая поровая структура бетона (только капиллярные поры); кривая 2 - двухуровневая поровая структура бетона (микропоры и капиллярные поры). I-III - различные стадии, выделяемые на кривых (пояснения в тексте) Fig. 2. The dependence of the compressive strength of the water-saturated concrete samples with different pore structures on the analogue of the dimensionless Darcy number: curve (1) - is a single-level pore structure of the concrete (only capillary pores); curve (2) - is the two-level pore structure of concrete (micropores and capillary pores). I-III - various stages highlighted on the curves (explanations in the text) Первый участок (при Da < 8) соответствует такой комбинации значений параметров образца, флюида и скорости нагружения, при которой жидкость успевает перераспределяться в объеме каналов и вытекать из образца, поддерживая поровое давление на низком уровне. В этом случае фильтрация (истекание) жидкости является дополнительным механизмом релаксации напряжений в каркасе, прочность образца в этой области значений параметров близка к максимальной и изменяется слабо. Максимальное значение прочности достигается при Da ® 0 и составляет ~64 МПа. Третий участок (при Da > 143) соответствует таким значениям комбинации параметров , h, k и W, при которых жидкость не успевает перераспределяться в капиллярной сети в процессе нагружения бетона. В этом случае пористость фактически является закрытой, и жидкость создает избыточное поровое давление, максимально понижая прочность бетона. Поэтому в интервале Da > 143 прочность бетона почти не изменяется и близка к минимальному значению ~56 МПа. Второй участок зависимости , при 8 < Da < < 143, является переходным, на нем соотношение скоростей процессов деформации каркаса и фильтрации поровой жидкости изменяется в широких пределах. Анализ разрушения образцов, характеризующихся различными значениями числа Дарси, свидетельствует об их однотипном характере (рис. 3). а б Рис. 3. Разрушение модельных водонасыщенных образцов бетона с капиллярной пористостью и отсутствующей микропористостью при различных значениях безразмерного параметра, аналогичного числу Дарси: а - Da = 0,6; б - Da = 143 Fig. 3. Fracture of the model water-saturated concrete samples with the capillary porosity (without microporosity) at different values of a dimensionless parameter similar to the Darcy number (a) - Da = 0.6, (b) - Da = 143 Отметим при этом, что с ростом числа Дарси, соответствующим уменьшению роли процессов фильтрации порового флюида и увеличению скорости роста порового давления в деформируемом образце, возрастает «множественность» разрушения. Действительно, по мере повышения Da увеличивается количество магистральных и соединяющих их вторичных (локальных) трещин, в отдельных областях разрушение приобретает характер дробления материала цементного камня. 3.3. Вклад микропористости в величину динамической прочности бетона с двухмасштабной поровой структурой Для изучения сравнительного вклада микропор в величину прочности образцов бетонов с двухмасштабной пористостью были рассмотрены образцы, построенные на основе описанных выше образцов с капиллярной пористостью путем явного добавления в их структуру микропористости (см. рис. 1). Как отмечалось выше, диаметр микропор составлял 500 мкм, величина микропористости - 4 % общего объема цементного камня. Микропоры и капиллярные поры полагались связанными (объединенными в единую фильтрационную систему). Анализ результатов моделирования показал, что зависимость прочности водонасыщенных бетонов с двухмасштабной пористостью от аналога числа Дарси, как и в случае одномасштабной пористости, аппроксимируется сигмоидальной кривой: (9) где » 43,3 МПа, » 49,9 МПа (кривая 2 на рис. 2). Можно видеть, что влияние микропористости на процессы фильтрации флюида в поровой сети и изменение прочностных характеристик бетона проявляются во всем рассмотренном диапазоне изменения значений параметра Da (что отличается от бетонов, обладающих только капиллярной пористостью). В первую очередь, важно отметить существенное общее снижение величины динамической прочности образцов на сжатие. Так, максимальное и минимальное значения прочности значения ( и ) флюидонасыщенного бетона с двухмасштабной пористостью уменьшились примерно на одинаковую величину (13-14 МПа или около 25 %) в сравнении с аналогичными значениями прочности бетона, цементный камень которого содержит только капиллярную сеть. Это является следствием того, что микропоры в структуре бетона являются сильными локальными концентраторами напряжений, обусловливающими более ранее зарождение повреждений и трещин в цементном камне при нагружении и формирование магистральной трещины при меньших значениях приложенной нагрузки. Данный эффект микропор не связан с поровым флюидом и поэтому в приблизительно равной степени проявляется как при высоких, так и при низких значениях скорости деформирования (т.е. при Da ® 0 и Da > 143). В то же время в области промежуточных значений Da динамика изменения величины прочности образцов от к отличается от той же для бетонов, имеющих только капиллярную пористость (кривые 1 и 2 на рис. 2). В частности, протяженность участка I зависимости сокращается с 0 < Da < 8 (бетон с одноуровневой пористостью) до 0 < Da < 0,4 (бетон с двухуровневой пористостью). В то же время стадия III (минимальные значения прочности) начинается при тех же характерных значениях числа Дарси, что и в случае бетона с одноуровневой поровой системой (при Da > 143). Таким образом, появление пор крупного масштаба (микропор) обусловливает увеличение протяженности переходной стадии II более чем на порядок величины (с 8 < Da < 143 до 0,4 < Da < 143). Описанное изменение зависимости является следствием перераспределения жидкости между поровыми пространствами различных масштабов. Действительно, крупные поры (микропоры) являются сильными локальными концентраторами напряжений и деформируются сильнее капиллярных пор. Соответственно, в процессе деформирования часть поровой жидкости переходит в объем капиллярных пор, повышая в них поровое давление. Ввиду того что характерное расстояние между микропорами много меньше полуширины образца, транспорт порового флюида между микро- и капиллярным поровыми пространствами начинает оказывать влияние на напряженное состояние и прочность бетона уже в области значений Da, соответствующей малым скоростям деформации. Следствием сказанного является более «ранний» переход от стадии I к стадии II (снижение прочности), т.е. динамическая прочность бетона с микропорами начинает существенно снижаться при меньших скоростях деформирования. Последнее можно интерпретировать как меньшую стойкость бетонов, обладающих двухмасштабной поровой структурой, к динамическим нагрузкам в сравнении со статическими нагрузками. Анализ картин разрушения образцов бетона с двухуровневой пористостью (рис. 4) показал, что оно носит значительно более локализованный характер в сравнении с бетоном с одноуровневой пористостью (см. рис. 3). а б Рис. 4. Разрушение модельных водонасыщенных образцов бетона с двухуровневой поровой структурой при различных значениях безразмерного параметра, аналогичного числу Дарси: а - Da = 0,4; б - Da = 143. Микропоры интегрированы в капиллярную поровую сеть Fig. 4. The fracture of the model water-saturated concrete samples with a two-level pore structure at different values of a dimensionless parameter similar to the Darcy number: a - Da = 0.4, b - Da = 143. Micropores are integrated into the capillary pore network Данная закономерность (склонность к локализации разрушения) является типичной для хрупких материалов, содержащих в структуре относительно крупные включения с пониженными механическими характеристиками (в частности, поры). Исследование влияния процессов перераспределения поровой жидкости на особенности разрушения показало, что при малых значениях Da (соответствующих малым скоростям деформирования) система трещин, как правило, более разветвлена, чем при больших Da (соответствующих большим скоростям деформирования). Отметим также, что при малых Da разрушение агрегатов носит преимущественно интерфейсный характер («отслоение» включений от матрицы на отдельных участках межфазной границы), а сами агрегаты, как правило, остаются целыми. При больших значениях Da многократно возрастает количество поврежденных или разрушенных агрегатов. Таким образом, в водонасыщенных бетонах с двухуровневой пористостью величина динамической прочности и особенности разрушения определяются как объемным содержанием крупных пор, так и перераспределением поровой воды в капиллярной поровой системе (последнее характеризуется числом Дарси). В частности, зависимость прочности образцов от параметра Da аналогична той же в образцах с капиллярной пористостью, однако характеризуется другими численными значениями прочности и интервалов Da, отвечающих основным стадиям. Так, наличие в бетоне микропористости, интегрированной в капиллярную поровую сеть, увеличивает вклад порового флюида в снижение прочностных свойств в области Da, соответствующей квазистатическому нагружению. Влияние интегрированности микропор в капиллярную поровую сеть высокопрочных водонасыщенных бетонов исследовалось на основе сопоставления результатов моделирования для образцов с микропорами, сообщающимися с капиллярной поровой сетью, и образцов с «закупоренными» микропорами. Во всех случаях микропоры, как и капиллярные поры, были заполнены жидкостью. Анализ механического отклика бетонов с закупоренными микропорами показал, что динамическая прочность образцов при сжатии также хорошо аппроксимируется сигмоидальной функцией, как и в случае бетонов с микропорами, интегрированными в капиллярную поровую сеть: (10) где » 38,3 МПа, » 50,3 МПа. Сопоставление зависимостей для этих бетонов (рис. 5) показало, что при малых значениях Da, когда роль процессов фильтрации порового флюида пренебрежимо мала, связность микро- и капиллярного поровых пространств не оказывает влияние на величину прочности. Вклад данного фактора становится значимым на стадии II (при Da > 1,43-2,86) - прочность на сжатие бетона с изолированными микропорами снижается примерно вдвое быстрее, чем в случае бетона со связанными между собой поровыми пространствами. На стадии III (при Da > 143), соответствующей «большим» скоростям нагружения бетона и минимальному вкладу фильтрации флюида через капиллярную поровую сеть, разница значений прочности достигает 13-15 %. Таким образом, «изоляция» микропор от капиллярной поровой сети снижает динамическую прочность бетона, но при этом почти не изменяет интервалов значений комбинации параметров, входящих в Da, отвечающих различным стадиям кривой Рис. 5. Зависимость прочности на сжатие водонасыщенных бетонных образцов с двухмасштабной поровой структурой от величины аналога числа Дарси. Кривая 1 - микропоры интегрированы в капиллярную поровую сеть; кривая 2 - микропоры изолированы от капиллярной поровой сети («закупорены») Fig. 5. The dependence of the compressive strength of the water-saturated concrete samples with a two-scale pore structure on the value of the analog of the Darcy number. Curve (1) the micropores are integrated into the capillary pore network, curve (2) the micropores are isolated from the capillary pore network (“clogged”) Анализ картин разрушения бетонов подтвердил обнаруженные закономерности. Так, при малых значениях Da (Da < 2,86) характер разрушения бетонов с «изолированной» и «неизолированной» микропористостью почти одинаков (см. рис. 4, 6). При Da = 143 наблюдается различие, касающееся пути распространения разрушающих трещин и фрагментации образца. Так, образцы бетона с «закупоренными» микропорами разделяются на большое число приблизительно равноосных фрагментов разного размера (рис. 6, б). Разрушение бетона с микропорами, интегрированными в поровую сеть, носит существенно более локализованный характер - разделение на крупные фрагменты, вытянутые в направлении оси нагружения (рис. 4, б). а б Рис. 6. Разрушение модельных водонасыщенных образцов бетона с двухуровневой поровой структурой при различных значениях безразмерного параметра, аналогичного числу Дарси: а - Da = 0,4; б - Da = 143. В обоих образцах микропоры «изолированы» от капиллярной сети Fig. 6. The fracture of the model water-saturated concrete samples with a two-level pore structure at different values of a dimensionless parameter similar to the Darcy number: (a) - Da = 0.4, (b) - Da = 143. In both samples micropores are isolated from the capillary pore network Таким образом, изолированность микропор оказывает значимое влияние и на особенности разрушения в области значений Da, которая соответствует достаточно большим скоростям нагружения (Da > 2,86). Выводы На основе компьютерного моделирования механического поведения водонасыщенного бетона при одноосном сжатии выявлены два фактора, определяющие величину динамической прочности тяжелых водонасыщенных бетонов с двухмасштабной поровой структурой. Первый фактор - мобильность поровой жидкости в сети капиллярных пор, которая определяет скорость снижения и выравнивания напряжений в пористом каркасе за счет фильтрации. Он определяет нелинейный характер изменения прочности со скоростью нагружения и может быть эффективно охарактеризован логистической зависимостью динамической прочности от числа Дарси (Da). Второй фактор - связанность крупных микропор с капиллярной сетью. Он обусловливает снижение концентрации напряжений в микропорах (наиболее сильных концентраторах напряжений) за счет фильтрации «избыточной» поровой жидкости в капиллярную поровую сеть и одновременно способствует возрастанию концентрации напряжений в окружающем микропоры цементном камне (за счет роста порового давления в капиллярах). Действие данного фактора определяет расширение скоростного интервала, в котором величина прочности существенно зависит от скорости нагружения, более чем на порядок величины. Вклады перечисленных факторов в изменение величины динамической прочности водонасыщенного бетона являются аддитивными и приблизительно равнозначными, а их суммарный вклад при «высоких» скоростях деформирования достигает 25 %. Таким образом, прочность водонасыщенных бетонов определяется не только содержанием микропор как наиболее крупных концентраторов напряжений, но в не меньшей степени - их интегрированностью в капиллярную пористость. Рассмотренное в настоящей работе условие одноосного сжатия (его можно рассматривать как эквивалент трехосного испытания образца в сосуде с жидкостью, находящейся при атмосферном давлении) является намного более «мягким» в сравнении с типичными условиями функционирования представительных объемов бетона. Данное исследование является только первым шагом в направлении изучения прочности и характера разрушения представительных объемов водонасыщенных бетонов в реальных (стесненных) условиях нагружения. Стесненные условия нагружения реализуются как во внутренних областях нагружаемых массивных бетонных элементов, так и в приповерхностных областях бетона в окрестности пятен контактного взаимодействия с плавучими массивными объектами. Анализ поведения представительных объемов бетона в условиях стесненного сжатия и при различных значениях порового давления на границах (что отвечает различной глубине погружения в водную среду) является предметом дальнейших исследований. Результаты проведенного исследования обладают большой степенью общности и вносят важный вклад в развитие современных представлений о сложных многофакторных процессах, определяющих наступление предельного состояния флюидонасыщенных хрупких гетерогенных материалов с многомасштабной поровой структурой. В частности, эти результаты могут применяться при прогнозировании интенсивности процессов деградации и абразивного изнашивания поверхностных слоев морских бетонных сооружений, функционирующих в арктических условиях, при динамическом контактном нагружении морскими льдами и/или плавучими объектами.

About the authors

Ig. S Konovalenko

Institute of Strength Physics and Materials Science; National Research Tomsk Polytechnic University

E. V Shilko

Institute of Strength Physics and Materials Science; National Research Tomsk State University

Iv. S Konovalenko

National Research Tomsk Polytechnic University

References

Statistics

Views

Abstract - 621

PDF (Russian) - 434