УДК 630*377.4
Ю. А. Ким, М. Т. Насковец, Н. И. Жарков, В. И. Гиль
Белорусский государственный технологический университет
ВЗАИМОДЕЙСТВИЕ КОЛЕСНОГО ДВИЖИТЕЛЯ ХОДОВЫХ СИСТЕМ С ГРУНТОВЫМ ОСНОВАНИЕМ
В связи с ростом массы транспортных средств, вызванных требованиями повышения их производительности, возрастают нагрузки на дорожное покрытие со стороны ходовых систем. Это, в свою очередь, приводит к разрушению опорной поверхности, увеличению сил сопротивления передвижению и увеличению энергозатрат. Поэтому вопросы, связанные с процессом взаимодействия колеса с опорной поверхностью по-прежнему являются актуальными.
В статье рассмотрен вопрос взаимодействия движителя ходовых систем (колеса) с грунтовым основанием. При взаимодействии пневматического колеса с почвогрунтом происходит взаимная деформация контактирующих тел. Если известны форма и размеры жесткого колеса, то при заданной максимальной глубине его погружения можно определить как форму поверхности контакта, так и закон распределения вертикальных перемещений частиц почвы. Разработана математическая модель, описывающая этот процесс, отличающаяся от существующих тем, что форма поверхности контакта не задается, а рассчитывается исходя из условий взаимодействия.
В результате проведенных исследований получено решение, позволяющее определять величину и характер распределения напряжений, а также момент сопротивления качению колеса и его зависимость от свойств взаимодействующих тел.
Ключевые слова: ходовая система, колесный движитель, грунтовое основание, объемная деформация, поверхность контакта, гармоническая функция, ряды Фурье.
Yu. A. Kim, M. T. Naskovets, N. I. Zharkov, V. I. Gil
Belarusian State Technological University
INTERACTION OF THE WHEEL DRIVING SYSTEM WITH GROUND BASIS
In connection with the growth of energy saturation of agricultural machines, row tractors and special purpose tractors, their mass increases, and consequently the loads of their running systems increase on the supporting surface. As a consequence, there are such negative phenomena as overgrazing of soil above permissible norms, there are energy costs for movement.
The article considers the issue of the interaction of the propulsion system of the undercarriage systems (wheels) with the soil base. When the pneumatic wheel interacts with the soil, mutual deformation of the contacting bodies occurs. If the shape and dimensions of the hard wheel are known, then for a given maximum depth of its immersion, you can determine both the shape of the contact surface and the law of distribution of vertical movements of soil particles. A mathematical model has been developed that describes this process, which differs from the existing ones in that the shape of the contact surface is not specified, but calculated based on the interaction conditions.
As a result of the research, a solution was obtained that allows one to determine the magnitude and nature of the stress distribution, as well as the moment of rolling resistance of the wheel and its dependence on the properties of interacting bodies.
Key words: undercarriage, wheel mover, soil base, volumetric deformation, contact surface, harmonic function, Fourier series.
Введение. Исследователь, указывает профессор Орнатский Н. В. [1], посвящающий свои силы изучению грунтов, неизбежно сталкивается с большой сложностью процессов их формирования и с чрезвычайным разнообразием их свойств. Поэтому к изучению грунтов необходимо подходить со строгим учетом их индивидуальных особенностей. Однако было бы неправильным, ориентируясь на эту необходимость, не использовать возможности рассмотрения некоторых общих зависимостей, вытекающих из применения
к грунтам основных физических законов, имеющих всеобъемлющий характер и не исключающих из сферы своего действия ни одного материального тела, в том числе и грунта. Реальные грунты всегда лишь в принципе следуют механической схеме, так как многообразие факторов, влияющих на особенности поведения грунта в тех или иных условиях, не может быть учтено даже самой совершенной схемой, при этом, являясь приемлемой и безусловно полезной в своем, сравнительно узком, кругу частных применений.
В зависимости от скорости деформирования грунтовые основания могут вести себя по-разному, а именно при небольшом интервале времени действия нагрузки наиболее ярко проявляются упругие свойства, и наоборот, при длительном воздействии больших по величине статических нагрузок грунт проявляет реологические свойства.
Однако реальные почвогрунты не являются чисто упругими и чисто пластичными. Один и тот же почвогрунт при различных условиях деформации проявляет свойства как упругости, так и пластичности [2, 3]. Деформация почвогрунтов сопровождается разрушением частиц и структурных агрегатов, выжиманием воды и воздуха из пор. Относительная доля того или иного из приведенных выше явлений зависит от состояния почвогрунта. Сжатие воды и защемленние объемов воздуха определяют процессы упругой деформации, фильтрацию воздуха и воды, а также относительные сдвиги частиц и разрушение структурных элементов - необратимые [4, 5, 6].
При рассмотрении вопроса воздействия колесного движителя на грунтовое дорожное покрытие последнее, ввиду кратковременности действия нагрузки, можно с достаточной степенью точности рассматривать как упругое основание.
Основная часть. Несмотря на известные недостатки, модель упругого полупространства заняла ведущее место в задачах статики конструкций на упругом основании и широко применяется для исследования напряженно-деформированного состояния грунтового массива методами теории упругости. Поскольку в задачах динамики имеют место малые деформации, быстро изменяющиеся во времени, применение в этом случае к грунту схемы упругого полупространства с моделью линейно деформируемой среды является еще более оправданным, чем в задачах статики [7]. При малых нагрузках деформация прямо пропорциональна им, то есть грунт в этом случае проявляет себя как линейно деформируемое тело.
Взаимодействие колесного движителя с почвой сопровождается большими скоростями деформации грунта [8]. Значение динамического модуля определяется по формуле
Ед = Е (1+ ),
(1)
где О - коэффициент, зависящий от вида грунта, равен для связного 0,09-0,11, для несвязного 0,06-0,08; ук - скорость поступательного движения колеса. Справедливо при ук < 5,5 м/с.
В случае, если верхний слой опорной поверхности достаточно мощный, последнюю можно рассматривать как бесконечное полупространство.
В процессе взаимодействия каждая частица почвы перемещается по сложной траектории. При этом чем больше размеры колеса, тем больше преобладают вертикальные перемещения в направлении oz. Обозначим перемещения частиц почвы по направлениям осей координат ох, оу, oz соответственно п\, И2, из, а также для краткости Х1 = х, Х2 = у, хз = z. Перемещения и для упругого полупространства должны удовлетворять уравнениям равновесия Ляме [9].
и
Р1(у)
--V
В
С
Рис. 1. Схема расположения координат для однородного полупространства
80
Уиг + (к +1)— = 0, г = 1, 2, 3, (2)
8хг
где к = X / ц, а X и ц - коэффициенты Ляме:
X = -
Е о
Е
(1 + о )(1 - 2о); ^ 2 (1 + о)
здесь У - оператор Лапласа; о - коэффициент поперечной деформации; Е - модуль деформации почвы; 0 - объемная деформация.
0 = V ^ .
1Й 8Хг
Общее решение уравнений (2) имеет вид [9]
И =Фг
8хг
(3)
где фг - произвольные гармонические функции.
Функция Y определяется из уравнения 8z
к +1 ^ 5фг-
к + 3 5x
Почвогрунт по физическим свойствам не является однородной средой. Модуль деформации Е может считаться постоянной величиной только в пределах отдельных слоев. Представление грунтового основания в виде многослойной анизотропной среды с числом слоев больше двух делает математическую модель более адекватной, но приводит к значительным трудностям вычислительного характера, а также громоздким уравнениям для вычисления коэффициентов ряда.
Определим форму поверхности контакта с почвой при заданной наибольшей глубине погружения Н, жесткого колеса, имеющего форму и размеры рассматриваемого пневматического при давлении воздуха в шине Р, равном нулю. Перемещения частиц почвы на расстояние у от средней плоскости вращения колеса определятся из выражения
г (х, у ) = Н (у)-Я (у ) + л] [ Я (у )]2 - X2 , (4)
где Н (у) - наибольшая глубина погружения колеса в почву в том же сечении; Я (у) - радиус сечения, расположенного на расстоянии от средней плоскости вращения колеса.
На рис. 2 изображена схема к определению функции деформации почвы под жестким колесом.
X
X
Рис. 2. Схема к определению функции деформации
Искомые функции фг-, входящие в уравнение (3), являются гармоническими. Поэтому решение удобно представить в виде двойного тригонометрического ряда Фурье [10]. Функции, разложенные в ряд Фурье, являются периодическими, поэтому математически представляется, что вне области интегрирования (в областях, смежных с областью пятна контакта колеса с почвой) действуют такие же нагрузки (рис. 3).
Рис. 3. Схема действия нагрузки при области интегрирования, равной длине контакта
В действительности это не так. Поэтому область интегрирования необходимо расширить. Чем шире область интегрирования, тем меньше влияние соседних нагрузок А и С на напряжения и деформации от нагрузки В (рис. 4). По принципу Сен Венана напряжения и деформации в массиве убывают по мере удаления от места приложения нагрузки [11]. Вследствие этого область изменения х и y представим в виде прямоугольника со сторонами 2z и 2H: z = 4Rk, H = 4B, Rk и B - соответственно радиус и половина ширины колеса.
Разложим функцию деформации почвы под жестким колесом f(x, y) в двойной тригонометрический ряд Фурье по х и y.
f (х, У) = ~40 + 1É a0iCi (У) + 2É ak0Ск (Х) +
4 2 г =1 2 k=1
+2Éb0iSl(y) + ÉÉ[аыск(х)Сг(y) +
2 11l
к=1 1=1 1
+ bkiSi (У)ск (X) + 2X акоSk ( X) +
2 к=1
+X X КС ( y)Sk ( X) + Pk-S- ( У ) Sk ( X)] ,
к=1 i=1
где для краткости, обозначено
, , . ( knx ) , , ( knx ) Sk (X ) = sin I-—I, ck (x ) = COS I-—I ;
(5)
S, (y )=sin I f I, c- (y )=COS If
iny
F1(y)
F2(y)
Рис. 4. Схема действия нагрузки при области интегрирования большей, чем пятно контакта
z
z
z
z
Коэффициенты ряда вычисляются по квадратурным формулам:
где
1 z H
aki =—í í f (xyК (x)dxcг (y)dy;
- z - H
1 z H
ъи = -¡h í íf (x' y) ck (x )dxS< (y )dy;
- z - H
1 z H
акг =—í íf (x.y)c (x)dxSk (x)dx;
- z-H
1 z H
акг = — í í f (x. у )ci (x )dxSk (x )dx• (6)
- z - H
Для нахождения коэффициентов ряда необходимо значения f (x, y), найденные из (4), подставить в формулы (6).
Контур поверхности контакта ограничен двумя кривыми (рис. 1), уравнения которых x¡ = F¡ (y), x2 = F2 (y) получены из выражения (4) при условии, что f (x, y) = 0. Уравнения кривых ADC и ABC:
X] (у) = -/ 2( Як-ау 4)( Н 0 -ау4)8-(Н - у4)2 82; х2 (у) = у!(Як-ау4)2-(Як -Н0)2 ,
где Як - ау4 - уравнение профиля шины; 5 -коэффициент упругого восстановления почвы. Деформация почвы на поверхности происходит лишь в области контакта, пределы интегрирования в формулах (6) по х от ^ до z следует заменить на пределы от до ^2(у).
Значения функции / (х, у) в интервалах изменения х от ^ до Е\(у) и от ^г(у) до z равны 0. Аналогично по у: вместо от —Н до Н подставить от -В до В.
Представим гармоническую функцию фз(х, у, z) в виде ряда Фурье:
Фз ( y, z ) = Yy АгСг (y ) + BoiSi (У ) ( ) +
г=1
+ ЕЁ[ АкгСг (x )Сг (y ) + B^ (x ) S (y )]e (z ) +
k=1 г=1
+ Y[CkoSk (x) + DkoCk (x)]ek (z) +
k=1
+ Y Y (y )Sk (x) + DklSt (y )Sk (x)] ekz. (7)
k=1 г =1
e. (z) = exp| -Hn 1 z
;(z)=exp
- + -
H2
ek (z) = exP^- Y)z
(8)
Коэффициенты ряда функции фз(х, у) находятся с учетом условия
Фз (x, У,0) = f (x, y) + c ( + Pw )
(9)
Функция / (х, у) определяется из уравнений (5) и (6). Остаётся найти фз(х, у, 0) и Р .
Для этого Р надо выразить через фз(х, у, 0),
используя уравнение объемной деформации 00 при z = 0.
Давление контакта определим по формуле
Pzzo =^00 + 2ц1
5z
где 00 - объемная деформация.
Из трёх функций фг положим ф2, поскольку перемещениями частиц грунта в направлении оси у на поверхности, ввиду их малости, пренебрегаем. Оставшихся двух функций достаточно, чтобы удовлетворить остальные краевые условия. Перемещения частиц почвы в контакте с пневматическим колесом по оси ох связаны с перемещением по оси oz уравнением ф^х, у, z) = = Рфз(х, у, z), где в - коэффициент пропорциональности, который определяется из условия (4). Касательные усилия на поверхности выражаются в интегральной форме, а не в каждой точке, поэтому такое представление функции ф1 является обоснованным.
После всех преобразований получим:
00 =(1 + Л)Ш( + (^)0]-
Pz. = (%)o + (^0 П,
6x
5z
(10)
где т = X(1 + А) + 2цА; п = (1 + A)(X + 2ц); А =- k + ¡ ; X, ц - коэффициенты Ляме.
k + 3
гЭфз.
Фз (x, У, 0) = f (x, y) + (—■3)0 pmc -dx
+ (nc + cPw. Значения подставим в условие (4).
m
Mf = & К ) (АР )к •
к=1
Полученное решение позволяет оценивать влияние Pw на величину момента сопротивления качению Ы/, глубину колеи h, величину и характер распределения давления в контакте Ргг0, а также на форму поверхности контакта пневматического колеса с почвой.
Заключение. В общем случае грунт является анизотропным, многослойным, вязко-упру-гопластичным телом. Учет названных свойств делает математическую модель более адекватной, но при этом решение задачи значительно усложняется. В различных условиях свойства грунта проявляются по-разному. При расчете напряженного состояния грунтового массива, являющегося основанием фундамента строительного сооружения при длительном воздействии больших величин статических нагрузок, грунт проявляет такие свойства, как пластичность, ползучесть. При кратковременном действии небольших по величине нагрузок грунт ярче проявляет упругие свойства. Именно такой режим нагружения характерен для процесса взаимодействия ходовых систем транспортных средств с грунтовым основанием, когда глубина распространений напряжений невелика. В этом случае грунтовое основание можно считать упругим изотропным полупространством. В результате проведенных исследований получено решение, позволяющее определять величину и характер распределения напряжений, а также момент сопротивления качению колеса и его зависимость от свойств взаимодействующих тел.
Литература
1. Орнатский Н. В. Механика грунтов. М.: Изд-во МГУ, 1962. 435 с.
2. Гуськов В. В. Тракторы. Часть II. Теория. Минск: Вышэйшая школа, 1977. 384 с.
3. Кацыгин В. В. Основы теории выбора оптимальных параметров мобильных сельскохозяйственных машин и орудий // Вопросы земледельческой механики. 1975. Т. 8. С. 3-142.
4. Бабков В. Ф., Безрук В. М. Основы грунтоведения и механика грунтов. М.: Высшая школа, 1976. 325 с.
5. Золотаревская Д. И. Зависимость между сжимающими напряжениями и осадкой почв // Механизация и электрификация сельского хозяйства. 1980. № 2. С. 30-32.
6. Взаимодействие ходовых систем тракторов с почвой. В. В. Кацыгин [и др.] // Механизация и электрификация сельского хозяйства. 1983. № 1. С. 18-19.
7. Ксеневич И. П., Скотников В. А., Ляско М. И. Ходовая система - почва - урожай. М.: Агропром-издат, 1985. 304 с.
8. Скойбеда А. Т. Автоматизация ходовых систем колесных машин. Минск: Наука и техника, 1979. 277 с.
9. Филоненко-Бородич М. М. Теория упругости. М.; Л.: ОГИЗ, 1947. 295 с.
10. Романовский П. И. Ряды Фурье. М.: Наука, 1973. 334 с.
11. Седов Л. И. Механика сплошной среды. М.: МГУ, 1968. Ч. 1. 436 с.
12. Рабочев И. С., Камнев В. I., Маслов А. С. Уменьшение отрицательного воздействия мобильных агрегатов на почву // Вестник сельскохозяйственной науки. 1979. № 4. С. 90-94.
13. Гапоненко В. С. О путях снижения уплотняющего воздействия машинно-тракторных агрегатов на почву // Влияние сельскохозяйственной техники на почву. 1981. С. 56-61.
Введем обозначение: Ы = -40 + еР№, где а00 -
первый член ряда функции f (х, у). Сделано это с целью учета влияния Рк, а также для изменения величины Ы по глубине. Разложим постоянную в ряд.
Ы =+ 2 £ [^ (у) + ¿8* 0% (х) +
4 2 1 =1 г=1
1 ¿ЪЛ (у) + XXXX[8^ (х)с (у) +
2 1 =1 к=1 1 =1
+ У кА (у )Ск (х)] + 2 ¿е к 0 Лк (х) +
2 1=1
+ ¿¿екС (у)8к (х) + ЦкД (у)Лк (х).
к =1 1 =1
После всех преобразований найдем давление в контакте колеса с опорной поверхностью Р . Вертикальная нагрузка О определяется из
условия О = . Площадь пятна контак-
та разбивается на прямоугольники с площадями АРк, давление в центре каждого определяется как
О = £ (Р^) к (АР) к,
к=1
где т - число прямоугольников.
Момент сопротивления качению колеса определяется по формуле
14. Омельянов А. Е. О применении пневматических колес на сельхозмашинах // Сельхозмашина. 1948. № 5. С. 15-18.
15. Алексейчик Н. А., Будько Ю. В., Терехов Б. А. Повышение проходимости сельскохозяйственных машин. Минск: Урожай, 1979. 139 с.
References
1. Ornatskiy N. V. Mekhanika gruntov [Soil mechanics]. Moscow, MGU Publ., 1962. 435 p.
2. Guskov V. V. Traktory. Chast' II. Teoriya [Tractors. Part II. Theory]. Minsk, Vysheyshaya shkola Publ., 1977. 384 p.
3. Katsygin V. V. Fundamentals of the theory of choosing the optimal parameters of mobile agricultural machinery and implements. Voprosy zemledel'cheskoy mekhaniki [Agricultural mechanics issues], 1975, vol. 8, pp. 3-142 (In Russian).
4. Babkov V. F., Bezruk V. M. Osnovy gruntovedeniya i mekhanika gruntov [Fundamentals of soil science and soil mechanics]. Moscow, Vysshaya shkola Publ., 1976. 325 p.
5. Zolotarevskaya D. I. The relationship between compressive stresses and soil sediment. Mekha-nizatsia i elektrifikatsia sel'skogo khozyaystva [Agricultural mechanization and electrification], 1980, no. 2, pp. 30-32 (In Russian).
6. ^tsygin V. V., Orda А. I., Аfanas'yev N. I. The interaction of tractor driving systems with soil. Mekhanizatsiya i elektrifikatsia sel'skogo khozyaystva [Agricultural mechanization and electrification], 1983, no. 1, pp. 18-19 (In Russian).
7. Ksenevich I. P., Skotnikov V. A., Lyasko M. I. Khodovya sistema - pochva - urozhay [Running system - soil - harvest]. Moscow, Agropromizdat Publ., 1985. 304 p.
8. Skoybeda А. Т. Avtomatizatsyia khodovykh system kolesnykh mashin [Automation of undercarriage systems for wheeled vehicles]. Minsk, Nauka i tekhnika Publ., 1979. 277 p.
9. Filonenko-Borodich M. M. Teoriyauprugosti [Theory of elasticity]. Moscow; Leningrad, OGIZ Publ., 1947. 295 p.
10. Romanovskiy P. I. Ryady Fur'ye [Fourier Series]. Moscow, Nauka Publ., 1973. 334 p.
11. Sedov L. I. Mekhanika sploshnoy sredy. Chast' I [Continuum mechanics. Part 1]. Moscow, MGU Publ., 1968. 436 p.
12. Rabochev I. S., Kamnev V. I., Maslov A. S. Reduction of the negative impact of mobile aggregates on soil. Vestnik sel'skokhozyaystvennoy nauki [Bulletin of Agricultural Science], 1979, no. 4, pp. 90-94 (In Russian).
13. Gaponenko V. S. On ways to reduce the sealing effect of machine-tractor aggregates on the soil. Vliyaniye sel'skokhozyaystvennoy tekhniki na pochvu [Influence of agricultural machinery on soil], 1981, pp. 56-61 (In Russian).
14. Omel'yanov A. Ye. On the use of pneumatic wheels on agricultural machines. Sel'khozmashina [Agricultural machinery], 1948, no. 5, pp. 15-18 (In Russian).
15. Alekseychik N. A., Budko Yu. V., Terekhov B. A. Povysheniye prokhodimosti sel'sko-khozyaystvennykh mashin [Raising the passability of agricultural machines]. Minsk, Urozhay Publ., 1979. 139 p.
Информация об авторах
Ким Юрий Алексеевич - кандидат технических наук, доцент кафедры инженерной графики. Белорусский государственный технологический университет (220006, г. Mинск, ул. Свердлова, 13а, Республика Беларусь). E-mail: y.kim@belstu.by
Насковец Михаил Трофимович - кандидат технических наук, доцент, профессор кафедры лесных машин, дорог и технологий лесопромышленного производства. Белорусский государственный технологический университет (220006, г. Mинск, ул. Свердлова, 13а, Республика Беларусь). E-mail: naskovets@belstu.by
Жарков Николай Иванович - кандидат технических наук, ведущий инженер кафедры лесных машин, дорог и технологий лесопромышленного производства. Белорусский государственный технологический университет (220006, г. Mинск, ул. Свердлова, 13а, Республика Беларусь). E-mail: zharkov@belstu.by
Гиль Виталий Иванович — кандидат технических наук, старший преподаватель кафедры инженерной графики. Белорусский государственный технологический университет (220006, г. Mинск, ул. Свердлова, 13а, Республика Беларусь). E-mail: v.gil@belstu.by
Information about the authors
Kim Yuriy Alekseуevich - PhD (Engineering), Assistant Professor, the Department of Engineering Drawing. Belarusian State Technological University (13a, Sverdlova str., 220006, Minsk, Republic of Belarus). E-mail: y.kim@belstu.by
Naskovets Mikhail Trofimovich - PhD (Engineering), Associate Professor, Professor, the Department of Logging Machinery, Forest Roads and Timber Production Technology. Belarusian State Technological University (13a, Sverdlova str., 220006, Minsk, Republic of Belarus). E-mail: naskovets@belstu.by
Zharkov Nikolay Ivanovich - PhD (Engineering), Leading Engineer, the Department of Logging Machinery, Forest Roads and Timber Production Technology. Belarusian State Technological University (13a, Sverdlova str., 220006, Minsk, Republic of Belarus). E-mail: zharkov@belstu.by
Gil Vitaliy Ivanovich - PhD (Engineering), Senior Lecture, the Department of Engineering Drawing. Belarusian State Technological University (13a, Sverdlova str., 220006, Minsk, Republic of Belarus). E-mail: v.gil@belstu.by
Поступила 15.10.2019