ВведениеАктуальность проблемы. Комплексный анализ современных подходов к моделированию процесса уплотнения грунтов представляет стратегическую важность для реализации масштабных инфраструктурных проектов в Российской Федерации. В условиях активного развития транспортной системы и реализации национальных проектов "Безопасные качественные дороги" и "Комплексное развитие территорий" проблема эффективного уплотнения грунтов приобретает особую значимость [1]. Современные требования к качеству строительства дорожных одежд и земляных сооружений диктуют необходимость разработки новых, более точных методов прогнозирования поведения грунтов в процессе уплотнения [2,3].Особую актуальность данное исследование получает в свете задач по импортозамещению и развитию отечественных технологий в дорожном строительстве. Отсутствие надежных расчетных моделей, адаптированных к российским грунтовым условиям и современной уплотняющей технике, приводит к значительным экономическим потерям и снижению долговечности строительных объектов [4,5].Уплотнение грунтов представляет собой сложный физико-механический процесс, характеризующийся нелинейным поведением трехфазной дисперсной системы. Фундаментальная сложность задачи заключается в необходимости одновременного учета реологических свойств грунтовой массы, кинетики изменения структурных связей и динамики перераспределения фазовых компонентов [6,7].Исторический анализ позволяет выделить три основных этапа развития теории уплотнения [8-10]. Эмпирико-статистический этап характеризовался накоплением экспериментальных данных и разработкой первых нормативных зависимостей, где ключевым достижением стало введение метода стандартного уплотнения Проктора [9]. Аналитический этап связан с развитием теорий контактного взаимодействия на основе решений Герца-Беляева [11], однако эти модели оказались ограниченно применимы для описания поведения упруго-пластичных грунтовых сред. Современный комплексный этап ознаменовался появлением зонных моделей Калужского [12,13] и развитием компьютерного моделирования, однако даже эти подходы не позволяют адекватно описать процесс уплотнения на всех его стадиях.Существующие теоретические модели демонстрируют принципиальные расхождения с экспериментальными данными по достижению предельного состояния уплотнения [14-16]. Наблюдается противоречие между дискретным характером моделей и непрерывностью реального процесса уплотнения, а также между статическими параметрами моделей и динамической природой изменения свойств грунта.Исследование предлагает системный подход к анализу ограничений существующих моделей и обоснованию необходимости разработки адаптивной модели, учитывающей эволюцию деформационных характеристик грунта и накопление пластических деформаций. Практическая значимость работы определяется возможностью создания на ее основе интеллектуальных систем контроля уплотнения и оптимизированных технологических карт производства работ.Цель исследования - проведение критического анализа существующих теоретических моделей контактного взаимодействия вальца с грунтом и разработка концептуальных основ новой адаптивной модели, позволяющей корректно прогнозировать поведение грунта в процессе многократного нагружения.Материалы и методы исследованияМетодологической основой исследования выступил системный анализ современных теоретических представлений о механике процесса уплотнения грунтов. В работе применен комплексный подход, интегрирующий методы сравнительного моделирования, математического анализа и верификации теоретических положений экспериментальными данными.Научная новизна методологии заключается в разработке критериев оценки адекватности существующих расчетных моделей, основанных на сопоставлении теоретических предсказаний с результатами натурных испытаний. Для установления корреляционных зависимостей между параметрами уплотняющего воздействия и деформационными характеристиками грунтов использованы аппарат нелинейной регрессии и методы вариационной статистики.Экспериментальная составляющая исследования базируется на метаанализе репрезентативной выборки научных публикаций, охватывающей широкий спектр грунтовых условий и технологических режимов уплотнения. Верификация теоретических моделей осуществлялась посредством сравнения с данными инструментальных измерений контактных напряжений и деформационных параметров.Особенностью примененной методологии является учет реологических свойств грунтовой среды и их трансформации в процессе многократного нагружения. Для анализа кинетики изменения структурных связей использованы современные представления механики деформируемого твердого тела с учетом специфики трехфазных дисперсных систем.Математический аппарат исследования включает методы решения контактных задач теории упругости и пластичности, адаптированные для описания поведения грунтовых сред. Для оценки достоверности полученных результатов применены критерии статистической значимости и методы проверки репрезентативности экспериментальных данных.РезультатыПроведенное комплексное исследование механизма уплотнения грунтов позволило установить фундаментальные закономерности, определяющие эффективность процесса уплотнения в зависимости от физико-механических свойств грунтовых сред и параметров технологического воздействия. Системный анализ существующих теоретических моделей в сопоставлении с экспериментальными данными выявил необходимость разработки уточненных расчетных зависимостей, адекватно описывающих поведение грунта при многократном нагружении.В результате критического осмысления классических теорий контактного взаимодействия и современных подходов к моделированию процесса уплотнения сформирована целостная система уравнений, учитывающая реологические свойства грунтов, кинематические параметры уплотняющей техники и особенности распределения контактных напряжений в трехфазной дисперсной среде. Особое внимание уделено верификации теоретических положений экспериментальными данными, что обеспечило достоверность полученных результатов и возможность их практического применения.Разработанные модели и зависимости позволяют перейти от эмпирических оценок к научно обоснованному прогнозированию поведения грунтов в процессе уплотнения, что создает теоретическую базу для оптимизации технологических параметров и разработки энергоэффективных методов уплотнения в дорожном и промышленном строительстве.Экспериментальное определение оптимальных параметров уплотнения грунтовПроведенные исследования стандартного уплотнения по Проктору позволили установить четкие корреляционные зависимости между типом грунта и его уплотняемостью. Полученные данные (Таблица 1) демонстрируют существенные различия в поведении различных грунтовых сред под нагрузкой.Таблица 1 Результаты определения оптимальных параметров уплотнения грунтовTable 1Results of determination of optimal soil compaction parameters / Результаты определения оптимальных параметров уплотнения грунтовТип грунта / Soil typeОптимальная влажность, % / Optimal moisture content, %Максимальная плотность, г/см³ / Maximum density, g/cm³Диапазон изменения плотности / Density variation rangeПесчаный / Sandy7-111,60-1,95 / 1.60-1.95±0,18 г/см³ / ±0.18 g/cm³Супесчаный / Sandy loam9-141,65-1,95 / 1.65-1.95±0,15 г/см³ / ±0.15 g/cm³Суглинистый / Loam13-191,60-1,85 / 1.60-1.85±0,13 г/см³ / ±0.13 g/cm³Тяжелый суглинистый и глинистый / Heavy loam and clay18-241,55-1,70 / 1.55-1.70±0,08 г/см³ / ±0.08 g/cm³  Анализ распределения контактных напряженийЭкспериментальные исследования распределения контактных напряжений под рабочими органами уплотняющих машин выявили сложный характер взаимодействия в системе "валec-грунт". На рисунке 1 представлена схема распределения напряжений под колесом, демонстрирующая наличие трех характерных зон деформации.Установлено, что максимальные контактные напряжения достигаются в зоне пластических деформаций, где наблюдается наиболее интенсивное уплотнение грунта. При этом распределение напряжений по площади контакта имеет куполообразный характер с максимумом, смещенным в направлении движения катка.Верификация теоретических моделей контактного взаимодействияСравнительный анализ шести теоретических моделей показал различную степень их соответствия экспериментальным данным. Наибольшую точность продемонстрировала модифицированная теория Герца-Беляева, описываемая следующими зависимостями:Максимальное контактное напряжение определяется как:                              (1)где коэффициенты упругости рассчитываются по формуле: ,                                             (2)   Рисунок 1. Зависимость скорости изменения напряженного состояния грунта и размеры в плане поверхности контакта рабочих органов машин с уплотняемым грунтом [19]Figure 1. Dependence of the rate of change of the stressed state of the soil and the dimensions in terms of the contact surface of the working bodies of the machines with the compacted soil   Для практического применения предложена уточненная степенная зависимость:                            (3)  - коэффициенты пропорциональности, определяющие степень податливости материалов цилиндра и полупространства.Коэффициенты , зависят только от рода материалов. В общем случае коэффициент С входит в уравнение (10) где:  - напряжение;  - величина сжатия в данной точке.Трехзонная модель процесса уплотненияДетальное исследование распределения контактных давлений подтвердило адекватность трехзонной модели Калужского. На рисунке 2 представлена схема распределения нормальных контактных давлений, иллюстрирующая выделение трех характерных зон. Для каждой зоны установлены соответствующие расчетные зависимости:Для зоны I (сдвига):       (4)Для зоны II (пластических деформаций):                                          (5) где:К - безразмерный коэффициент.Для зоны III (разгрузки):                                              (6)где:К - коэффициент пропорциональности;  - показатель степени‚ зависящий от рода материала сдавливаемых тел. Рисунок 2. Схема распределения напряжений под колесом [19]Figure 2. Voltage distribution scheme under the wheelРазработка практических расчетных методикДля инженерных расчетов были предложены упрощенные эмпирические зависимости. Установлено, что среднее и максимальное контактные напряжения могут быть определены с достаточной для практики точностью по формулам:,                            (7)где:  - коэффициент, подобный коэффициент К в уравнении (6) В [18] при исследовании зависимости между погружением металлического шара в грунт и нагрузок‚ установлено‚ что эта зависимость подчинена следующей закономерности:,                                               (8)где:  - величина погружения шарового штампа;  - коэффициент пропорциональности.Исследование влияния скоростных факторовОсобое внимание было уделено исследованию влияния скорости на процесс уплотнения. Установлено, что изменение давлений во времени может быть описано уравнением:,                                            (9)где: Д - диаметр шара.Для определения эффективной скорости уплотнения предложены зависимости:                  (10)                           (11)где  - коэффициент перехода от среднего напряжения к максимальному‚ который зависит от характера распределения напряжений под штампом) аналогичны формуле, полученной на основе теории Герца-Беляева.Дополнительные аналитические зависимостиВ рамках исследования также были получены другие важные зависимости. Установлена связь между нагрузкой и погружением штампа:.                   (12)                  (13)здесь знак + относится к первой зоне, а знак - ко 2 и 3 зонам;  толщина уплотняемого слоя грунта под осью катка; толщина уплотняемого слоя грунта в месте ;  - толщина слоя грунта в исследуемой точке контакта. В первой зоне  может принимать значения от Н1 до Нк.. Во второй зоне Нк   может принимать значения от  до Нк. Нк лежит на границе 1-й и 2-й зон.  - коэффициент трения частиц грунта о каток; - коэффициент трения частиц грунта об основание; среднее напряжение сжатия; q - нормальное удельное давление катка на поверхность слоя;  - коэффициент, характеризующий распределяющую способность слоя, он зависит от глубины и расстояния от оси катка.На основе анализа уравнения (13) распределение нормальных контактных давлений изобразится эпюрой, изображенной на рисунке 3.  Рисунок 3. Схема распределения нормальных контактных давлений [26]Figure 3. Normal contact pressure distribution scheme according to Ya.A. Kaluzhsky Для среднего контактного напряжения при действии шарового штампа:             (14)