Scientific journal
International Journal of Applied and fundamental research
ISSN 1996-3955
ИФ РИНЦ = 0,593

STUDY ON PROCESS OF CREEP OF SOIL AFTER SHORT-TERM INCREASE OF DYNAMIC LOAD

Zolotarevskaya D.I. 1
1 Russian State Agrarian University – Moscow Timiryazev Agricultural Academy
2191 KB
Creep processes in cohesive soils occurring after the initial impact on these soils loads for sinusoidal as well as linear laws. Viscoelastic properties of the soil describes first-order differential equation that relates the compressive stress, the rate of change of compressive stresses and relative deformation of the compression. Considered viscoelastic properties of soddy-podzolic light loamy soil identified in field tests. The study took into account that the density of soil is a function of the depth of deformable layer. Based on statistical processing experimental data, dependence of initial soil density as a function of depth adopted quadratic. Performed computer experiments, in which found indicators that can estimate quantify the change of the deformation, stress, soil density when its creep after short-term increase for dynamic loads. Identified impact on these figures the initial soil density, humidity, increasing load speed and other factors.
viscoelastic properties of soils
deformation
compressive stress
soil compaction
creep
mathematical modeling

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

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

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

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

При влажностях w почвы, меньших ее полной влагоемкости, и при сжимающих напряжениях σ, меньших предела ее прочности σпр, почва под действием нагрузки уплотняется и упрочняется. Неуплотненные почвы при таких w являются вязкоупругими [3, 9]. Настоящая работа посвящена исследованию и математическому моделированию ползучести и уплотнения почв, находящихся в таком состоянии.

Время t взаимодействия с почвой движителей мобильных машин и рабочих органов сельскохозяйственных машин и орудий при выполнении полевых работ весьма мало. Поэтому важно исследовать изменение вязкоупругих свойств и плотности почв при таких малых t.

Ранее нами в работе [11] был предложен метод расчета, позволяющий находить показатели напряженно-деформированного состояния почвы в различных условиях нагружения и приращение плотности почвы на различной глубине, возникающее в результате действия внешней нагрузки; исследован процесс ползучести почвы.

Цели данной работы состояли в следующем: развить полученные в работе [11] результаты исследования и математического моделирования деформирования и уплотнения почвы с учетом ее вязкоупругих свойств; уточнить и применить предложенный метод расчета показателей напряженно-деформированного состояния и уплотнения почвы, изменяющихся при кратковременном динамическом возрастании нагрузки, а также расчета показателей последующего процесса ползучести почвы.

Материалы и методы исследования

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

Результаты исследований и их обсуждение

Результаты экспериментального исследования вязкоупругих свойств почвы

На полях опытного хозяйства РГАУ – МСХА им. К.А. Тимирязева были проведены экспериментальные исследования реологических свойств, закономерностей деформирования и уплотнения почв в результате действия динамических нагрузок при работе колесных тракторов МТЗ-82 [4–9]. В данной работе в расчетах использованы экспериментальные данные, полученные в работе [4]. Содержание этих опытов отражено также в работе [10] . Исследуемая почва – дерново-подзолистая легкосуглинистая. На экспериментальном поле верхний деформируемый слой почвы, распространенный на глубину H = 1,0 м, расположен на практически недеформируемом основании. Поверхность почвы горизонтальна. Проведен ряд серий опытов, каждая из которых состояла из шести последовательных проходов по почве трактора.

В работе [11] выполнено исследование для почвы, верхний деформируемый слой которой до воздействия на нее штамповой нагрузки распространен на глубину H = 0,51 м.

Введем ось Oy, направленную вертикально вниз, начало которой O расположим на поверхности почвы. Плотность верхнего деформируемого слоя почвы перед проходами по ней машин (начальная) переменна по глубине y.

В опытах, описанных в работах [4, 5] , на ряде участков поля, отобранных случайным образом, находили плотность ρ и влажность w почвы перед каждой серией опытов в слоях 0 – 0,1, ...., 0,9 – 1,0 м. Для получения аналитической зависимости ρ(y) в качестве представителя каждого слоя почвы принимали его середину. В результате статистической обработки экспериментальных данных из работ [4, 5] нами получены квадратичные зависимости плотности ρ почвы от глубины y

zol07.wmf,

zol08.wmf, (1)

где k1 и k2 ≠ 0 – коэффициенты, ρ0 – свободный член. Для почвы опытного поля (при w = 15,5 % , ρ0 = 1,1412 г/см3, k1 = 1,8762 г/(см3 м), k2 = 1,1374 г/(см3 м2); корреляционное отношение в зависимости (1) равно 0,9861. Высокое корреляционное отношение свидетельствует о весьма большой точности нахождения плотности почвы по формуле (1). При zol09.wmf почву можно уплотнить, так как ρ(y) < ρпр, ρ(H) = ρпр, где ρпр – плотность предельно уплотненной почвы при σ < σпр. В данной работе зависимость начальной плотности почвы от глубины y описана уравнением (1). Зависимость (1) применена также в работах [6–8].

На основании результатов исследований [4 – 9] закономерность сжатия почвы в направлении оси Oy при каждом фиксированном y будем моделировать дифференциальным уравнением

zol10.wmf, (2)

где ε – относительная деформация сжатия; t – время, с; zol12.wmf и zol13.wmf, с-1 и МПа – характеристики реологических (вязкоупругих) свойств почвы. Границы применимости этого уравнения для конкретных почв выявляют по экспериментальным данным. Адекватность уравнения (2) при моделировании закономерности сжатия легкосуглинистой почвы известного гранулометрического состава при w = 15–26 % и при определенных значениях t, σ, ε, а также супесчаной почвы подтверждена экспериментально [4–9].

Характеристика p, с -1 зависит от круговой частоты ω, с-1 гармонического процесса деформирования, возникающего, в частности, в результате деформирования почвы при качении цилиндрического штампа или колеса. В расчетных формулах, полученных на основе определяющего уравнения (2), ее целесообразно представить в виде zol15.wmf, где g – преобразованная безразмерная характеристика вязкоупругих свойств почвы.

Путем статистической обработки большого числа экспериментальных значений g и q, полученных на опытном поле при различных ρ, ω и w, найдены уравнения регрессии для определения при каждом фиксированном y характеристик g и q исследованной почвы:

zol16.wmf (3)

zol17.wmf (4)

с коэффициентами множественной корреляции, равными соответственно 0,7931 и 0,7528. Теснота связей в этих уравнениях достаточно большая, поэтому они могут быть применены в расчетах для исследованной почвы и для других почв, имеющих близкие к ней физические свойства. Уравнения (3) и (4) соответствуют следующим интервалам значений влияющих факторов: ρ = 1,14–1,88 г/см3, w = 15–26 %, ω = 0,93–5,01 с-1. В уравнениях (3) и (4) при y = 0 плотность почвы ρ = ρ0.

Результаты теоретического исследования ползучести почвы после кратковременного возрастания динамической нагрузки

Пусть на поверхность почвы посредством штампа передается внешняя нагрузка. Ширина штампа принята достаточно большой, в этом случае деформацию почвы можно считать плоской. На основе математического моделирования процесса распространения в почве при ее нагружении вязкоупругих затухающих волн деформации сжатия почвы дифференциальным уравнением с частными производными четвертого порядка и с переменными коэффициентами нами были получены формулы, позволяющие найти глубину Hp распространения деформации сжатия почвы, приращение плотности и плотность почвы на различной глубине. Выполненные нами исследования показали, что при H ≤ 1 м глубина Hp = H, поэтому в данной работе приняли Hp = H = 1 м. Относительная деформация сжатия почвы в каждый момент времени t равна

zol19.wmf, (5)

где zol20.wmf – переменная абсолютная деформация сжатия (осадка) почвы.

В теории вязкоупругости чаще рассматривают ползучесть в деформируемых средах при постоянных уровнях напряжений, создаваемых в результате ступенчатого статического изменения нагрузок, то есть принимают, что каждый определенный уровень напряжений достигается мгновенно [12]. Нами учтено, что внешняя нагрузка прикладывается к почве не мгновенно, а возрастает в течение некоторого времени.

В механике грунтов выявляют в основном их строительные свойства, характеризующие деформирование грунтов как оснований фундаментов сооружений, поэтому ползучесть грунтов исследуют в длительные промежутки времени, измеряемые часами и днями [2, 14]. Нами исследована ползучесть почв в малые промежутки времени, возникающая после возрастания сжимающих напряжений σ в почве в секунды и их доли.

Ползучесть представляет собой процесс изменения во времени деформаций в материале или в деформируемой среде при постоянном напряжении. С целью исследования процесса ползучести нами рассмотрены две стадии изменения напряженно-деформированного состояния почвы. Первая стадия соответствует времени zol21.wmf. На этой начальной стадии изменяются напряжения в почве и ее деформации: zol22.wmf, zol23.wmf. При zol24.wmf σ и ε возрастают: zol25.wmf, zol26.wmf.

Вторая стадия изменения напряженно-деформированного состояния почвы (стадия ползучести) соответствует zol27.wmf. На этой стадии zol28.wmf, а деформации с течением времени изменяются: zol29.wmf.

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

1. Начальное возрастание сжимающих напряжений по синусоидальному закону и последующая ползучесть почвы.

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

zol30.wmf, zol31.wmf, (6)

где σm – амплитуда напряжения. Подставив zol32.wmf в уравнение (2) и учитывая начальное условие: при zol33.wmf относительная деформация сжатия почвы ε = 0, найдем

zol34.wmf,

zol35.wmf. (7)

Получили формулу для определения приращения плотности почвы на различной глубине. При Hp = H и квадратичной зависимости (1) начальной плотности почвы от глубины приращение плотности почвы на глубине zol36.wmf приближенно равно

zol37a.wmf

zol37b.wmf, (8)

где μ – коэффициент поперечного расширения почвы, Приняли, что на глубине H почва до действия на нее внешней нагрузки имела максимально возможную плотность, поэтому zol38.wmf.

В результате действия внешней нагрузки с течением времени t изменяются деформации почвы, глубина распространения деформации сжатия почвы, характеристики g и q ее вязкоупругих свойств. Осадка почвы и новая глубина распространения ее деформируемого слоя перед началом стадии ползучести соответственно равны: zol39.wmf и zol40.wmf. Характеристики вязкоупругих свойств почвы перед началом ее ползучести вследствие изменения плотности почвы равны g1 и q1.

б) Стадия ползучести почвы соответствует zol41.wmf. На этой стадии zol42.wmf; zol43.wmf. Подставив zol44.wmf в уравнение (2) и учитывая начальное условие: при zol45.wmf относительная деформация zol46.wmf, получим формулу, характеризующую процесс ползучести почвы:

zol47.wmf, zol48.wmf. (9)

Величину ε0, входящую в формулу (9), находим по формуле (7) при значении zol49.wmf.

Уравнение (9) при постоянных значениях его параметров описывает линейную ползучесть. Оно показывает, что при zol50.wmf деформации сжатия почвы с течением времени неограниченно возрастают по линейному закону. Такого вида закономерность изменения деформации при ползучести характерна для идеальной вязкоупругой среды Максвелла. Однако почва не является идеальной вязкоупругой средой. Из опытов известно, что при постоянных σ < σпр почва с течением времени уплотняется и упрочняется до постоянных значений плотности и прочности, ее деформации стабилизируются [3]. Эта закономерность уплотнения почвы при ползучести учтена в разработанном нами методе расчета.

На стадии ползучести за новое начало отсчета глубины – величины yн (то есть за новую поверхность почвы) принимаем координату zol51.wmf. Плотность почвы при yн = 0 в начале процесса ползучести равна

zol52.wmf. (10)

Зная приращения плотности почвы на различной глубине, найдем новые значения ее плотности ρн(yн) для всех yн∈[0, Hн]. В соответствии с полученными ранее результатами [8, 9] зависимость плотности почвы от глубины yн аппроксимирована квадратичной функцией вида (1), но с измененными входящими в нее параметрами.

При изменении плотности почвы изменяются ее вязкоупругие свойства. Эти свойства определяются новыми значениями характеристик g и q ее вязкоупругих свойств. Характеристики g и q вязкоупругих свойств исследовавшейся почвы находим по формулам (3) и (4) при zol53.wmf, то есть при zol54.wmf. С течением времени в связи с уплотнением почвы характеристики ее вязкоупругих свойств изменяются.

2. Начальное возрастание сжимающих напряжений по линейному закону и последующая ползучесть почвы.

а) Пусть сжимающие напряжения, передаваемые на поверхность почвы посредством плоского штампа и действующие в направлении оси Oy, со временем zol55.wmf возрастают по линейному закону

zol56.wmf, zol57.wmf, (11)

где v – линейная скорость изменения сжимающих напряжений, м/с; M – коэффициент пропорциональности, кПа/м. В расчетах получено, что численно М = 1. Если zol58.wmf, то zol59.wmf; если zol60.wmf, то zol61.wmf; zol62.wmf. Рассмотрим деформирование почвы вдоль оси zol63.wmf.

Подставив выражение скорости изменения напряжений zol64.wmf в уравнение (2) и учитывая начальное условие: при zol65.wmf относительная деформация zol66.wmf, найдем,

zol67.wmf,

zol68.wmf. (12)

б) Стадия ползучести почвы соответствует zol69.wmf. На этой стадии zol70.wmf; zol71.wmf. Процесс ползучести почвы характеризуется формулой (9) с переменными значениями характеристик вязкоупругих свойств почвы. Величину ε0, входящую в эту формулу, находим в данном случае по формуле (12) при значении zol72.wmf.

Характеристики g и q, а также g1 и q1 вязкоупругих свойств почвы и показатели ее уплотнения при zol73.wmf находим аналогично тому, как это сделано для случая начального возрастания сжимающих напряжений по закону (6).

Метод расчета изменения показателей напряженно-деформированного состояния почвы в результате кратковременного начального возрастания динамической нагрузки и последующей ползучести

На основе выполненного математического моделирования нами разработан метод расчета показателей изменения напряженно-деформированного состояния почвы и показателей ее уплотнения в случае кратковременного возрастания при zol74.wmf начальных сжимающих напряжений σ(t) по законам (6) и (11) и последующей ползучести почвы при зависимости (1) начальной плотности почвы от глубины. В расчетах учитывается, что деформация стабилизируется в течение времени tст стабилизации деформации в результате того, что принимается во внимание изменение с течением t < tст плотности ρ(y) почвы и, следовательно, характеристик g(y) и q(y) ее вязкоупругих свойств.

Исходными для выполнения каждого расчета являются следующие данные: Н, глубина zol75.wmf, zol76.wmf, zol77.wmf, zol78.wmf, w, σm, t0, коэффициенты уравнений регрессии zol79.wmf, zol80.wmf (для исследовавшейся легкосуглинистой почвы это коэффициенты уравнений (3) и (4)), ω (в случае возрастания σ(t) при zol81.wmf по синусоидальному закону (6)) или v (в случае возрастания σ(t) при zol84.wmf по линейному закону (11)).

Характерные особенности предложенного метода расчета состоят в следующем. Величины деформации сжатия почвы, характеристик g и q ее вязкоупругих свойств, показатели напряженно-деформированного состояния и уплотнения почвы при действии сжимающих напряжений в почве являются непрерывными функциями t. Приближенно считаем, что в малые промежутки времени zol85.wmf (zol86.wmf эти величины постоянны и их изменение происходит дискретно. В расчетах приняли zol87.wmf 0,2 с.

Рассмотрим малые значения времени ti. На первой стадии изменения напряженно-деформированного состояния почвы при zol88.wmf находим в случае возрастания сжимающих напряжений по синусоидальному закону (6): напряжение zol89.wmf, zol90.wmf (по формуле (7) при zol91.wmf), zol92.wmf, zol93.wmf по формуле (8), zol94.wmf. zol95.wmf по формуле (10), zol96.wmf и zol97.wmf – старое и новое значения плотности почвы в слое 0 – 0,1 м, характеристики g1 и q1 вязкоупругих свойств почвы при zol98.wmf – новом значении ρ0. Для zol99.wmf – нового значения глубины, соответствующей глубине Н1, новая плотность почвы равна zol100.wmf.

На второй стадии изменения напряженно-деформированного состояния почвы (стадии ползучести) приближенно считаем, что сжимающие напряжения постоянны в малые промежутки времени zol101.wmf. Рассматриваем моменты времени zol102.wmf, zol103.wmf, zol104.wmf, zol105.wmf, ..., zol106.wmf. При каждом zol107.wmf, соответствующем ползучести почвы, вычисляем, воспользовавшись формулой (9), приращение zol108.wmf относительной деформации сжатия почвы, zol109.wmf – приращение абсолютной деформации сжатия почвы, приращение суммарной относительной zol110.wmf и абсолютной zol111.wmf деформации сжатия почвы, новую глубину распространения деформируемого слоя почвы, находим новую зависимость вида (1) плотности почвы от глубины. При zol112.wmf и zol113.wmf = 0 имеем

zol114.wmf, (13)

где zol115.wmf – относительная деформация сжатия почвы (суммарная).

С течением t характеристика q возрастает, а характеристика g убывает. При этом zol116.wmf, zol117.wmf, zol118.wmf. Деформацию почвы считаем стабилизированной, если при zol119.wmf получим zol120.wmf, в других слоях плотность почвы также достигает своего предельного значения ρпр. Если в m-ом расчете получим zol121.wmf, то это означает, что далее с течением времени плотность почвы практически не увеличивается, наступает условная стабилизация ее свойств.

Аналогично выполняются расчеты в случае линейного закона (11) возрастания сжимающих напряжений при zol122.wmf и последующей ползучести почвы.

Результаты компьютерных экспериментов. Анализ полученных результатов

Нами разработаны компьютерные программы, позволяющие найти по полученным формулам и алгоритмам показатели, характеризующие напряженно-деформированное состояние и уплотнение почвы в различные фиксированные моменты времени при кратковременном начальном динамическом возрастании сжимающих напряжений в течении времени zol123.wmf и последующей ползучести почвы. По этим программам выполнены компьютерные эксперименты, в которых найдены указанные показатели для дерново-подзолистой легкосуглинистой почвы, вязкоупругие свойства которой исследованы в опытах, описанных в работах [4, 5]. Эти свойства, которые отражены также в данной работе, описывает уравнение (2). Характеристики вязкоупругих свойств почвы перед началом ее деформирования находим по уравнениям регрессии (3) и (4). Опыты по исследованию ползучести почвы выполнили для двух режимов начального деформирования почвы при zol124.wmf: 1) по синусоидальному закону (6) (опыты группы А); 2) по линейному закону (11) (опыты группы Б).

В случае начального возрастания сжимающих напряжений по синусоидальному закону (6) исследовали влияние на ползучесть и уплотнение почвы: 1) ее начальной плотности, 2) круговой частоты ω изменения сжимающих напряжений при zol125.wmf, 3) влажности w почвы, 4) времени t0 начального возрастания сжимающих напряжений, 5) амплитуды напряжений σm. Для выявления характера и количественной оценки влияния этих факторов выполнили базовый расчет и пять серий однофакторных компьютерных экспериментов (серии 1, А – 5, А).

Базовый расчет провели по следующим исходным данным: H1 = 0,3 м; H = 1 м; ρ0 = 1,14 г/см3; ρ(H1) = 1, 602 г/см3; ρ(H) = 1,88 г/см3; w = 15 %; t0 = 0,2 с; σm = 100 кПа; ω = 2 с-1. В каждой из пяти серий компьютерных экспериментов варьировали по одному из основных влияющих факторов, а остальные факторы принимали такими же, как в базовом расчете.

В случае начального возрастания сжимающих напряжений по линейному закону (11) исследовали влияние на ползучесть и уплотнение почвы: 1) ее начальной плотности, 2) скорости v возрастания сжимающих напряжений при zol126.wmf, 3) влажности w почвы, 4) времени t0. Для выявления характера и количественной оценки влияния этих факторов выполнили базовый расчет и четыре серии однофакторных компьютерных экспериментов (серии 1, Б – 4, Б). Базовый расчет провели по следующим исходным данным: H1 = 0,3 м; H = 1 м; ρ0 = 1,14 г/см3; ρ(H1) = 1, 602 г/см3; ρ(H) = 1,88 г/см3; w = 15 %; t0 = 7 с; v = 3 м/с. В каждой из четырех серий компьютерных экспериментов варьировали по одному из основных влияющих факторов, а остальные факторы принимали такими же, как в базовом расчете.

В опытах серий 1, А и 1, Б выявили влияние начальной плотности почвы на исследуемые показатели. Каждая из этих серий состояла из пяти опытов. Начальную плотность почвы в каждом из опытов этих серий принимали по данным табл. 1.

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

На рис. 1 и рис. 2 представлены соответственно зависимости относительной деформации сжатия почвы и приращения относительной деформации сжатия почвы от времени t в процессе ее ползучести в случае начального возрастания сжимающих напряжений при zol127.wmf по синусоидальному закону (6) при различных значениях начальной плотности почвы.

zolot1.wmf

Рис. 1. Зависимости при ползучести относительной деформации сжатия почвы от времени при различных значениях начальной плотности почвы (w = 15 %, начальное возрастание сжимающих напряжений по синусоидальному закону (6), σm = 100 кПа, ω = 2 с-1, t0 = 0,2 с). 1, 2, 3, 4 и 5 – при ρ(y) по данным опытов 1–5, приведенным в табл. 1

zolot2.wmf

Рис. 2. Зависимости при ползучести приращения относительной деформации сжатия почвы от времени при различных значениях начальной плотности почвы (w = 15 %, начальное возрастание сжимающих напряжений по синусоидальному закону (6), σm = 100 кПа, ω = 2 с-1, t0 = 0,2 с), 1, 2, 3, 4 и 5 – при ρ(y) по данным опытов 1–5, приведенным в табл. 1

Таблица 1

Начальная плотность почвы в ряде опытов

№ Опыта

Плотность почвы при различных y, г/см3

ρ0

ρ(0,5)

ρ (1,0)

1

2

3

4

5

1,14

1,2

1,3

1,4

1,5

1,602

1,62

1,65

1,70

1,75

1,88

1,88

1,88

1,88

1,88

Соответствующие кривые получены нами также для случая изменения σ(t) при zol128.wmf по линейному закону (11). Представленные на рис. 1 зависимости ε(t), а также другие, полученные в результате расчетов данные, показывают, что при одном и том же значении σ0 большим значениям начальной плотности почвы соответствуют меньшие значения ε0. С возрастанием начальной плотности почвы снижаются величины ε(t), соответствующие одному и тому же значению t, снижаются величины предельного значения относительной деформации сжатия εпред почвы.

В табл. 2 приведены данные, отражающие результаты опытов серии 1, Б. Опыты характеризуют деформации и уплотнение почвы при различных значениях ее начальной плотности в случае начального возрастания сжимающих напряжений при zol129.wmf по линейному закону (11) и последующей ползучести до условной стабилизации деформации. В этой таблице даны начальные значения плотности почвы ρ(0,05) в слое 0–0,1 м, численные значения показателей, полученных при возрастании напряжений по линейному закону: относительных деформаций сжатия почвы ε0, плотности ρн(0,05) почвы в слое 0–0,1 м при t0 = 7 с. Приведены также значения показателей, характеризующих последующую ползучесть почвы: новое значение плотности почвы ρн(0,05) в слое 0–0,1 м при t = 20 с, время условной стабилизации деформаций tст, ρст(0,05), а также приращения ∆ρ(0,05) при tст, ∆ε(tст).

Таблица 2

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

Показатели

При ρ0, г/см3

1,14

1,2

1,3

1,4

1,5

В результате возрастания сжимающих напряжений в течение времени t0 = 7 c

σ0, кПа

21

21

21

21

21

ε0

0,2679

0,2256

0,1849

0,1476

0,1258

ρ(0,05), г/см3

1,2534

1,2856

1,4439

1,4533

1,5458

h0, см

0,2113

0,1841

0,1514

0,1286

0,1117

ρн(0,05), г/см3 при t0 = 7 c

1,2636

1,2897

1,3731

1,4620

1,5313

В результате последующей ползучести

ρн(0,05), г/см3

при t = 20 с

1,4715

1,4805

1,4923

1,5434

1,6009

tст, с

108

107

105

91

69

ρст(0,05), г/см3

1,6762

1,6730

1,6699

1,6698

1,6698

∆ ρ(0,05), г/см3

При tст

4,04∙10–4

4,03∙10–4

4,03∙10–4

4,03∙10–4

4,02∙10–4

∆ε(tст)

1,14∙10–4

1,14∙10–4

1,18∙10–4

1,17∙10–4

1,21∙10–4

Таблица 3

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

Показатели

При скорости v, м/с

1

2

4

В результате возрастания напряжений в течение времени t0 = 7 с

σ0, кПа

7

14

28

ε0

0,6045

0,6147

0,9675

ρ(0,05), г/см3

1,2648

1,2648

1,2648

h0, см

3,767

3,793

4,917

ρн(0,05), г/см3 при t0 = 7 с

1,4039

1,4577

1,4682

В результате последующей ползучести

ρн(0,05), г/см3

при t = 20 с

1,4921

1,5745

1,5834

tст, с

274

144

35,5

ρст(0,05), г/см3

1,78

1,75

1,72

∆ρ(0,05), г/см3

при tст

4,61∙10–4

6,86∙10–4

7,73∙10–4

∆ε(tст)

1,03∙10–4

1,54∙10–4

2,36∙10–4

Примечание. Начальная плотность почвы – по данным опыта 1 из табл. 1.

В опытах серии 2, А выявили влияние на исследуемые показатели круговой частоты ω в случае начального возрастания сжимающих напряжений при zol132.wmf по синусоидальному закону (6) . В опытах серии 2, Б выявили влияние на исследуемые показатели линейной скорости v в случае начального возрастания сжимающих напряжений при zol133.wmf по линейному закону (11). Каждая из этих серий состояла из пяти опытов.

В табл. 3 приведены данные, отражающие результаты опытов серии 2, Б. Опыты характеризуют деформации и уплотнение почвы при одинаковых значениях ее начальной плотности в случае начального возрастания сжимающих напряжений при zol134.wmf по линейному закону (11) и последующей ползучести до условной стабилизации деформации. В этой таблице даны начальные значения плотности почвы ρ(0,05) в слое 0–0,1 м, численные значения показателей, полученных в случае возрастания сжимающих напряжений по закону (11): относительных деформаций сжатия почвы ε0, плотности ρ(0,05) почвы в слое 0–0,1 м при t0 = 7 с. Приведены также значения показателей, характеризующих последующую ползучесть почвы: новое значение плотности почвы ρн(0,05) в слое 0–0,1 м при t = 20 с, время условной стабилизации деформаций tст, zol135.wmf, а также приращения ∆ρ(0,05) при tст, zol136.wmf. Аналогичные данные получены в опытах серии 2, А.

Результаты расчетов показывают, что с ростом v, а также с ростом ω относительные деформации почвы возрастают до своего предельного значения εпред, а величины εпред и время tст снижаются. На рис. 3 и рис. 4 приведены кривые, характеризующие соответственно зависимости приращения ∆ρ(0,05) плотности почвы и плотности ρ(0,05) почвы в слое 0–0,1 м от времени t при различных значениях круговой частоты ω в случае начального возрастания сжимающих напряжений по синусоидальному закону (6). Из рисунков видно, что увеличение ω приводит к возрастанию ∆ρ(0,05) и ρ(0,05).

zolot3.wmf

Рис. 3. Зависимости при ползучести приращения плотности почвы в слое 0–0,1 м от времени при различных значениях круговой частоты ω (w = 15 %, начальное возрастание сжимающих напряжений по синусоидальному закону (6), σm = 100 кПа, ρ(y) – по данным опыта 1 из табл. 1, t0 = 0,2 с), 1, 2 и 3 – при ω = 1; 2 и 4 с-1

zolot4.wmf

Рис. 4. Зависимости при ползучести плотности почвы в слое 0–0,1 м от времени при различных значениях круговой частоты ω (w = 15 %, начальное возрастание сжимающих напряжений по синусоидальному закону (6), σm = 100 кПа, ρ(y) – по данным опыта 1 из табл, 1, t0 = 0,2 с), 1, 2, 3 и 4 – при ω = 1; 3, 4 и 5 с-1

В опытах серий 3, А и 3, Б нами получены данные, характеризующие изменение исследуемых показателей во времени при различных значениях влажности w почвы. Результаты расчетов показывают, что с ростом w деформации почвы возрастают до εпред, а tст снижается.

На рис. 5 представлены кривые, характеризующие зависимости плотности ρ(0,05) почвы в слое 0–0,1 м от времени t при различных значениях влажности w почвы. Как видно из этого рисунка, при одних и тех же значениях t большим значениям w соответствуют большие значения плотности почвы; увеличение w приводит к возрастанию ρ(0,05) и снижению tст.

zolot5.wmf

Рис. 5. Зависимости при ползучести плотности почвы в слое 0–0,1 м от времени при различных значениях влажности w почвы (начальное возрастание сжимающих напряжений по синусоидальному закону (6), σm = 100 кПа, ω = 2 с-1, t0 = 0,2 с, ρ(y) – по данным опыта 1 из табл. 1), 1, 2, 3 и 4 – при w = 15, 17, 19 и 22 %

Анализ полученных результатов показал, что основное влияние на все рассмотренные показатели оказывает начальная плотность почвы. Существенное влияние оказывают также скорость возрастания сжимающих напряжений при zol138.wmf и влажность почвы.

Отмеченный характер влияния начальной плотности почвы на ее ползучесть и уплотнение обусловлен тем, что при увеличении плотности почвы характеристика вязкоупругих свойств почвы q возрастает, причем zol139.wmf, где E – модуль упругости почвы, а характеристика g убывает, причем zol140.wmf. Свойства почвы приближаются к упругим.

Отмеченный характер влияния влажности почвы на ее ползучесть и уплотнение обусловлен тем, что при увеличении w характеристика вязкоупругих свойств почвы g возрастает, а характеристика q убывает, причем zol141.wmf. Свойства почвы приближаются к текучим.

Заключение

На основе математического моделирования вязкоупругих свойств почвы дифференциальным уравнением (2) получены аналитические зависимости и алгоритмы, позволяющие расчетным путем находить показатели напряженно-деформированного состояния почвы при различных законах возрастания с течением времени zol142.wmf сжимающих напряжений, а также рассчитывать приращение плотности почвы на различной глубине, возникающее в результате возрастания сжимающих напряжений при zol143.wmf и при ползучести почвы. Исследован процесс ползучести почвы после кратковременного начального возрастания сжимающих напряжений. Расчеты выполнены с использованием полученных нами в полевых испытаниях экспериментальных данных о вязкоупругих свойствах дерново-подзолистой легкосуглинистой почвы известного гранулометрического состава в определенном интервале изменения ее влажности и скорости возрастания сжимающих напряжений при zol144.wmf.

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

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

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