CC BY
60
УДК 691.33
ТЕОРЕТИЧЕСКИЕ АСПЕКТЫ СТРУКТУРНО-ИМИТАЦИОННОГО МОДЕЛИРОВАНИЯ МАКРОСТРУКТУРЫ КОМПОЗИЦИОННЫХ СТРОИТЕЛЬНЫХ МАТЕРИАЛОВ
THEORETICAL ASPECTS OF STRUCTURAL AND SIMULATION MODELING OF THE MACROSTRUCTURE OF COMPOSITE BUILDING MATERIALS
Адилходжаев А.И., Махаматалиев И.М., Шаумаров С.С. Adilhodzhaev A.I. , Mahamataliev I.M., Shaumarov S.S.
Ташкентский институт инженеров железнодорожного транспорта (Ташкент, Узбекистан) Tashkent institute of railway engineers (Tashkent, Uzbekistan)
Аннотация. В настоящее время описание совместной работы разнородных по свойствам фаз многокомпонентных высококачественных бетонов в рамках структурно-имитационного моделирования цементных композитов является весьма актуальной задачей строительного материаловедения. В контексте решения вышеуказанных задач в данной статье приводятся теоретические аспекты структурно-имитационного моделирования на примере разработанной программы «Бетон-технология», позволяющей прогнозировать прочностные свойства многокомпонентного высококачественного бетона на макроструктурном уровне.
Ключевые слова: структурно-имитационное моделирование, композиционные строительные материалы, высококачественный бетон, композит, макроструктура, прочность.
Дата принятия к публикации: 15.07.2018
Дата публикации: 25.09.2018
Сведения об авторах:
Адилходжаев Анвар Ишанович - доктор технических наук, профессор, проректор по научной работе Ташкентского института инженеров железнодорожного транспорта; [email protected]
Махаматалиев Иркин Муминович - кандидат технических наук, профессор кафедры "Строительство зданий и промышленных сооружений" Ташкентского института инженеров железнодорожного транспорта; [email protected]
Шаумаров Саид Санатович - кандидат технических наук, доцент кафедры "Строительство зданий и промышленных сооружений" Ташкентского института инженеров железнодорожного транспорта; shoumarovss@gmail. com
Введение.
Несмотря на определенные достижения в производстве бетонов различного назначения и имеющийся значительный опыт их производства в Узбекистане и за рубежом ® ® 1
J Abstract. At the present time, the description of the
J joint work of multi-component, high-quality concretes
J in the framework of the structural and simulation mod-
J eling of cement composites is a very topical task of
J building materials science. In the context of the solu-
J tion of the above problems, theoretical aspects of struc-
J tural and simulation modeling are presented on the
J example of the developed program "Concrete-
J technology", which allows to predict the strength prop-
J erties of multi-component high-quality concrete at the
J macro structural level. J
J
J Keywords: structural simulation simulation, composite
J building materials, high-quality concrete, composite,
J macrostructure, strength. J
J
J Date of acceptance for publication: 15.07.2018
J Date of publication: 25.09.2018 J
J Authors' information:
J Anvar I. Adilkhodjaev - Doctor of Technical
J Sciences, Professor, Vice-Rector for scientific work at
J Tashkent Institute of Railway Transport Engineers
J (Uzbekistan); [email protected]
J Irkin M. Mahamataliev - Candidate of Technical
J Sciences, Professor of Department "Construction of
J buildings and industrial structures" at Tashkent Institute
J of Railway Transport Engineers (Uzbekistan);
J Said S. Shaumarov - Ph.D. (Eng.), Associate Professor
J of Department "Construction buildings and structures"
J at Tashkent Institute of Railway Engineers
J (Uzbekistan); [email protected] J
внедрение инновационных технологий в производство сборного и монолитного бетона и железобетона идёт очень медленно. В особенности это касается высокомарочных и многокомпонентных бетонов, обеспечиваю-
щих ресурсо- и энергосбережение как в процессе их производства, так и в период эксплуатации. Решение таких сложных задач связано, с одной стороны, с изысканием и расширением сырьевой базы материалов, имеющей свои особенности, оптимизацией составов, разработкой механизмов регулирования вопросов смесеобразования в увязке с технологическими режимами производства, а также с всесторонними исследованиями характера протекания структурных изменений на уровне микро- и макроструктуры и установлении взаимосвязи с физико-механическими и эксплуатационными показателями разрабатываемых материалов, с другой стороны.
Установление корреляционных зависимостей между параметрами формирующейся структуры и свойствами многокомпонентных (5-6-ти, а иногда и более составляющих) бетонов, технологических переделов получения материала с наперёд заданным комплексом свойств является весьма сложной задачей, для решения которой необходимо привлечение нетрадиционных методик, позволяющих охватить весь спектр задач: проектирование состава, исследования закономерностей и характера структурообразования, технологии приготовления, прогнозирования физико-механических и эксплуатационных показателей создаваемого материала. Попытки поиска научно обоснованного решения вышеуказанных задач с использованием неординарных методик привело к разработке и созданию полиструктурной теории композиционных строительных материалов [1, 2].
Полиструктурная теория композиционных строительных материалов как единая система научных представлений о закономерностях структурообразования, технологии и свойствах композиционных материалов строительного назначения получила заслуженное признание в научном мире.
Главная отличительная особенность полиструктурной теории состоит в том, что принцип полиструктурности выступает не только как классификационный фактор или методологический прием для объяснения механизма формирования структуры и свойств материала, но и как ключ к направленному исследо-
ванию возникающих между компонентами химических реакций с целью формирования требуемых физико-технических показателей материала, разработке рациональной технологии изготовления. В рамках полиструктурной теории композиционных строительных материалов впервые четко определены основные структурообразующие факторы для каждого структурного уровня и получены количественные зависимости свойств композитов от этих факторов [2].
Вопросы разработки строительного композита согласно основных постулатов полиструктурной теории композиционных строительных материалов рассматриваются сквозь призму последовательных этапов: физико-химические взаимодействия и превращения протекающие в формирующейся среде с учетом индивидуальных особенностей компонентов, их влияний на ход течения процессов формирования структуры, оптимизации технологических операций изготовления и т.д.
Такая многогранная задача, решение которой направлено на получение «идеального» композита с наперёд заданным комплексом свойств путем обработки большего массива переменных факторов, не может быть реализована обычным тривиальным подходом.
Анализ материалов исследований [3-6], проведенных к настоящему времени, показал, что реализация такой масштабной работы может быть выполнена путем моделирования структуры разрабатываемого композита. Полученные имитационные и виртуальные модели послужат основой для оптимизации физико-механических и технологических параметров проектируемого материала с учетом индивидуальных особенностей исходных компонентов и их объёмной концентрации, а использование математического аппарата и современных компьютерных программных комплексов в полной мере обеспечит решение оптимизационных задач путем варьирования большого числа переменных факторов [7, 8].
Моделирование позволяет получить представление о внутренних процессах, протекающих при формировании структуры, установить механизм влияния отдельных элементов на ход течения возникающих взаи-
модействий и объяснить особенности формирования структуры в целом, получить виртуальную картину формирующейся структуры композита [9-14]. В результате появляется возможность производить ранжирование и описание влияния различных внешних и внутренних факторов для формирования требуемых физико-механических, физико-технических и специальных свойств проектируемого материала.
Структурно-имитационное моделирование макроструктуры многокомпонентного высококачественного бетона.
В контексте вышеуказанных теоретических положений была разработана программа «Бетон-технология», позволяющая прогнозировать прочностные свойства тяжёлого многокомпонентного высококачественного бетона на макроструктурном уровне. Общий вид моделируемого объекта исследования приведен на рис. 1.
' ► х
М ' 1
Рис. 1. Моделируемый объект и его структурные компоненты
Как видно из рис. 1, в структуре исследуемого композита различают следующие элементы:
- крупный заполнитель;
- мелкий заполнитель;
- цементный камень;
- воздушные поры;
- различные дефекты.
При моделировании макроструктура бетона условно принята состоящей из двух блоков: крупного заполнителя и цементно-При моделировании макроструктура бетона условно принята состоящей из двух блоков: крупного заполнителя и цементно-песчаного раствора, а начальные дефекты макроструктуры бетона и его компонентов -с коллинеарными трещинами (рис. 2).
Рис.2. Геометрические параметры начальных дефектов структуры объекта
При создании модели учитывались следующие геометрические параметры начальных дефектов макроструктуры:
- г^ = const - радиус поры с координатами
Xi, Yj ;
- loij - начальная длина трещины, численно равная 0,184r;
- aij - ориентация трещины относительно вектора нагрузки (случайная величина, подчиняющаяся закону произвольного распределения в интервале от 0.. ,2п);
- Xi, Yj - координаты начальных дефектов макроструктуры, являющиеся независимыми случайными величинами, распределенными по равномерному закону по площади пластины (моделируемого образца);
- N - количество начальных дефектов в поле образца (Nmin = 30).
Включения (зерна крупного заполнителя) моделировались выпуклыми многогранниками (рис. 3.) со следующими геометрическими параметрами:
- Rij = const - радиус описанной вокруг многоугольника окружности;
- nij - число вершин многоугольника (n <6);
- © m - ориентация k-ой вершины многоугольника M-го включения относительно вектора нагрузки (случайная величина, распределенная по равномерному закону на интервале 0...2п);
- Xм, YM - координаты центра описанной окружности M-го включения;
- ф - концентрация включений;
- Кф - коэффициент формы включения.
^х
Y
Рис.3. Общий вид включений в моделируемом объекте
Концентрация включений ф определяется по формуле
м
S Sm m=1
S.
(1)
об
Коэффициент формы Кф определяется по формуле
KФ -
Lmax _ Lmax cos^1
Ашп ^тт 008 %2 где £тах, Lmin - максимальная и минимальная стороны многоугольника включения.
В модели учитывались следующие физические параметры включений:
- Ев - модуль упругости;
- [Лц - коэффициент Пуассона;
- К® - коэффициент интенсивности напряжений при нормальном разрыве;
- К^- коэффициент интенсивности напряжений при плоском сдвиге.
Контактная зона включения (сторона многоугольника включения с заданной шириной 3) является либо постоянной, либо случайной величиной, распределенной по произвольному закону.
Физические параметры включений:
- Кк1с - критический коэффициент интенсивности напряжений при нормальном разрыве;
- Кк11с - критический коэффициент интенсивности напряжений при плоском сдвиге принимались прямо пропорциональными аналогичному параметру матрицы:
V- к - Л Vм ■ - М 1с ;
Кк — Л КМ •
^11с ~ ^М^Пс 5
где MKTK, MKTM - микротвердость контактной зоны и матрицы соответственно (являются постоянными величинами).
В разработанной модели определяющим ист о чником, вызывающим деформацию матрицы и компонентов бетона, является внешняя нагрузка g.
При этом объектом анализа деформационных свойств бетона в модели являются:
- слияние изолированных дефектов и их развитие на границе с включениями;
- образование и развитие зигзаг-трещин, возникающих при огибании трещинами включений.
Лм -
ГМ
MKT
К
MKT
м
(2) (3)
где Sm - площадь M-го включения; Sa6 -площадь образца.
y
g
y
x
При моделировании на начальном этапе развития трещин (рис.4.) величина внешней нагрузки g , вызывающей развитие двух трещин изолированного начального дефекта макроструктуры, определяется по формуле (индекс г принимает значения из области матрицы, либо значение из области контактной зоны, либо значение из области включения):
л[ж¥^К1^а1 + а
*
Ч =
2
(4)
Концентрация включений ф для вершины
а2 (рис.4) составляет:
Ф(а2)= I
' Г2
3г
4\
а V
2%"
а}( г2
I
Зг
2%
4 Л
а1-%
а2 + %
х
2% 2% Я
а1-%
а2 + %
ё%. (5)
Условие слияния начальных дефектов с различными диаметрами определяется выражением
Ч> =
аг + аг-1,г+1 2с - (гг + Г-1.
> 0,46.
(6)
г +1)
Величина нагрузки, вызывающей развитие N (/ >2) слившихся начальных дефектов макроструктуры бетона или его компонентов (рис.5) описывается уравнением: - для вершины трещины а/
, =^пГгК1сУ^а1 + а
Ча/ =
У(а„)
для вершины трещины а\
* 4ж¥1К1с¥1^а1 + а
Ч 9(а1)
N
(7)
(8)
д —►
► X
-а
а.
М "2
Рис. 4. Начальный дефект макроструктуры бетона и его компонентов
Рис. 5. Геометрические параметры начальных дефектов макроструктуры бетона и
его компонентов в процессе слияния
4
д
Нагрузка д* , способствующая развитию
зигзаг-трещины, вычисляется по формулам:
- для вершины трещины а^
+ а
(р(ам, Ьм ) для вершины трещины а]
^КПоУ^а1 + а
где
Чщ Р(аъ К ) (Р(ам, ЬМ = (Р(ам ) +
N
+ -
К
м
К
£ Ат /, *).
Но т=1
^ - Ьт
aN + Ьт
ср(аь Ьм = рЫ +
+
К
Л м
\ К11о )
£ Ат («, /, *)
т=1
а1 - Ьп
(9)
(10)
(11)
(12)
aN + Ьт
К
ц= 1,15 Кпо
К
м о
(13)
при выходе трещины из матрицы во включение
Л=1,15
К
К 11о
К
(14)
Прямолинейное развитие трещины через включение (матрицу) наблюдается, если К!> Ка, I = 1, 2, 3.
Трещина останавливается на границе с включением, если
К! <К1а , I = 1, 2, 3, где Ка - принимается в зависимости от местоположения вершины трещины (I ~ В - во
включении; I ~ М - в матрице; I ~ К - в контактной зоне).
Ккоэффициент интенсивности нормальных разрывающих напряжений у вершины трещины К вычисляется по следующим формулам:
- для вершины трещины а]
д Р(а1, Ьм .
К! =
г1Л] л(а+^)
для вершины трещины aN
К =_ др(%> Ьм)
(15)
(16)
Развитие трещины на границе с включением моделируется следующим образом: при выходе трещины на грань включения (контактную зону) с ростом нагрузки g трещина будет развиваться прямолинейно (через включение) или тормозиться им.
Условие развитие трещины на границе с включением определяется параметром ] на основании неравенств:
- трещина огибает включение: -]< I;
- прямолинейное развитие трещины (К
>К1Ы): -]> I;
- остановка трещины контактной зоной включения (К < К!о1):
при выходе трещины из включения в матрицу
ж(а1 + %)
В разработанной структурно-имитационной модели процесса разрушения бетона и его компонентов на уровне макроструктуры формирующиеся взаимодействия между структурными элементами в процессе твердения бетона производятся по вышеуказанным формулам. В качестве дополнительного также принято условие, что трещина развивается параллельно приложенной нагрузке и ее рост прекращается при выходе на другую трещину.
Разработанная модель процесса разрушения цементного бетона и программное обеспечение «Бетон-технология» позволяют проследить эволюцию формирования деструктивных процессов, прогнозировать прочностные свойства тяжёлого бетона на макро-структурном уровне и получать виртуальное изображение динамики распределения напряжений в бетонном образце в процессе увеличения внешней нагрузки g вплоть до завершающего момента разрушения образца включительно (рис. 6).
Заключение.
Графическая визуализация результатов численного моделирования прочности цементного бетона с использованием программного комплекса «Бетон-технология» соответствует существующим представлениям о процессе разрушения цементных композиций и позволяет заключить, что предложенный метод структурно-имитационного моделирования процессов разрушения бетона при одноосном сжатии пригоден для теоретической оценки прочности рассматриваемой системы.
ДИНАМИКА РАЗРУШЕНИЯ БЕТОНА С ВКЛЮЧЕНИЯМИ (данные образца: А=400; В=300 мм; кол-во дефектов = 50)
A q=0
A q=5
Aq=10 мПа
Aq=15
A q=2 0
Aq=25
Физические характеристики
матрицы: K|Cm =2 MH/M3/2 Kpcm =8 MH/M3/2 Ем = 15000 мПа U = 0,25
Рз = 1000 кг/м3
Физические характеристики включений
kICb =1,5 МН/м3/2 Крсв = 6 МН/м3/2 Eb = 20000 мПа
Aq=46 - разлом
50 100 150 200 250 300 350
50 100 150 200 250 3С0 35
50 100 150 200 250 300 350
50 100 150 200 250 300 350
Рис. 6. Графическая визуализация результатов численного моделирования прочности цементного бетона с использованием программного комплекса «Бетон-технология»
Список литературы
1. Соломатов, В.И. Полиструктурная теория композиционных строительных материалов / В.И. Соломатов. - Ташкент: Фан, 1991. - 345 с.
2. Соломатов, В.И. Развитие полиструктурной теории композиционных строительных материалов / В.И. Соломатов // Изв. вузов. Строительство. - 1985. - № 8. - С.58-64.
3. Баженов, Ю.М. Достижения, проблемы и направления развития теории и практики строительного материаловедения / Ю.М. Баженов // Материалы X Академических чтений РААСН. Пенза, Казань, 2006. -С. 3-7.
4. Баженов, Ю.М. Основные подходы к компьютерному моделированию строительных композитов / Ю.М. Баженов // Строительные материалы. - 2006. - №7. - С. 2-4.
$ References
$ 1. Solomatov V.I. Polystructural theory of
$ composite building materials. Tashkent, Fan,
$ 1991. 345 p. (In Russian)
$ 2. Solomatov V.I. Development of the
$ polystructural theory of composite building
$ materials. Izvestiya vuzov. Stroitelstvo, 1985,
$ No.8, pp. 58-64. (in Russian)
$ 3. Bazhenov Yu.M. Achievements, prob-
$ lems and directions of development of theory
$ and practice of building materials science.
$ "Materialy XAkademicheskikh chteniy RAASN
$ [Materials of X academic readings of
$ RAASN]. Penza, Kazan, 2006, pp. 3-7. (In
$ Russian)
$ 4. Bazhenov Yu.M. Main approaches to
$ computer modeling of building composites.
$ Stroitelnye materially, 2006, No.7, pp. 2-4. (In
$ Russian)
5. Чермашенцев, В.М. Теоретические $ аспекты компьютерного моделирования $ эффективных композиционных материалов $ / В.М. Чермашенцев // Изв. вузов. Строи- $ тельство. - 2002. - №3. - С. 33-40. $
6. Харитонов, А.М. Исследование $ свойств цементных систем методом струк- $ турно-имитационного моделирования / $ А.М. Харитонов // Строительные материа- $ лы. - 2008. - №9. - С. 81-83. $
7. Adilhodjaev, A.I. About nature of inter- $ phase interactions of basalt aggregate and poly- $ carboxylate super plasticizer with cement in $ multi-component concrete / A.I. Adilhodjaev, $ I.M. Mahamataliev, V.M. Tsoy // "IBAUSIL" - $ "19 Internationale Baustofftaguung". - Wei- $ mar: 2015. - Р. 1211-1219. $
8. Адилходжаев, А.И. Вопросы управле- $ ния качеством при проектировании соста- $ вов многокомпонентных высококачествен- $ ных бетонов с минеральными наполнителя- $ ми / А.И. Адилходжаев, И.М. Махаматали- $ ев, В.М. Цой // Вестник ТашГТУ. - 2016. - $ №1. - С. 252-260. $
9. Воробьев, В.А. Основные задачи ком- $ пьютерного моделирования строительных $ композитов / В.А. Воробьев // Строитель- $ ные материалы. - 2006. - №7. - С. 19-21. $
10. Воробьев, В.А. Применение физико- $ математических методов в исследованиях $ свойств бетона / В.А. Воробьев. - М.: Выс- $ шая школа, 1977. - 345 с. $
11. Shaumarov, S.S. On the issue of in- $ creasing energetic efficiency of buildings in $ railway transport / S.S. Shaumarov // VIII In- $ ternational Conference "Transport Problems - $ 2016". - Katowice, 2016. - P. 522-532. $
12. Shaumarov, S.S. On peculiarities of $ formation of the thermal mode in operating $ panel buildings / S.S. Shaumarov, A.I. Adil- $ hodzhayev, F.F. Karimova // International con- $ ference «Science and practice: a new level of $ integration in the modern world». - London, $ 2018. P. 365-368. $
13. Shaumarov, S.S. On the method of es- $
timation of thermal protective properties of ex- $
ternal walls of operating residential buildings / $
S.S. Shaumarov, A.I. Adilhodzhayev, S.I. Kan- $
dakharov // International conference «Student $
science: research works». Part II. - San Fran- $
5. Chermashentsev V.M. Theoretical aspects of computer simulation of effective composite materials. Izvestiya vuzov. Stroitelstvo, 2002, No.3, pp. 33-40. (In Russian)
6. Kharitonov A.M. Investigation of the properties of cement systems by the method of structural-simulation simulation. Stroitelnye materialy, 2008, No. 9, pp.81-83. (In Russian)
7. Adilkhodjaev A.I., Makhamataliev I.M., Tsoy V.M. About nature of interphase interactions of basalt aggregate and polycarboxylate super plasticizer with cement in multi-component concrete. "IBAUSIL" 19 Internationale Baustofftaguung, Weimar, Bundesre-pablik Deutschland, 2015, pp.1211-1219.
8. Adilkhodjaev A.I., Mahamataliev I.M., Tsoy V.M. Questions of quality management in the design of multi-component high-quality concretes with mineral fillers. Vestnik Tash-GTU, 2016, No.1, pp. 252-260. (In Russian)
9. Vorobyev V.A. The main tasks of computer modeling of building composites. Stroitelnye materially, 2006, No.7, pp. 19-21. (In Russian)
10. Vorobiev V.A. Application of physical and mathematical methods in studies of concrete properties. Moscow, Vysshaya Shkola, 1977. 345 p. (In Russian)
11. Shaumarov S.S. On the issue of increasing energetic efficiency of buildings in railway transport. VIII International Conference "Transport Problems - 2016". Katowice, Poland, pp. 522-532.
12. Shaumarov S.S., Adilhodzhayev A.I., Karimova F.F. On peculiarities of formation of the thermal mode in operating panel buildings. International conference «Science and practice: a new level of integration in the modern world». London, 2018, pp. 365-368.
13. Shaumarov S.S., Adilhodzhayev A.I., Kandakharov S.I. On the method of estimation of thermal protective properties of external walls of operating residential buildings. International conference «Student science: research works». Part II. San Francisco, California, USA, 2018, pp. 103-106.
14. Adylhodzayev A.I., Shaumarov S.S. The issue of thermal renovation of infrastructure of railway transport is evaluated. X International Scientific Conference "Transport
cisco, 2018. - P. 103-106. $ Problems - 2018". Katowice, Poland, 2018,
14. Adylhodzayev, A.I. The issue of thermal $ pp. 13-18.
renovation of infrastructure of railway trans- $
port is evaluated / A.I. Adylhodzayev, S.S. $
Shaumarov // X International Scientific Con- $
ference "Transport Problems - 2018". - Ka- $
towice, 2018. - P. 13-18. $
