Расчет периода первого тона собственных колебаний растения злаковой культуры Текст научной статьи по специальности «Электротехника, электронная техника, информационные технологии»

УДК 631.354

РАСЧЕТ ПЕРИОДА ПЕРВОГО ТОНА СОБСТВЕННЫХ КОЛЕБАНИЙ РАСТЕНИЯ ЗЛАКОВОЙ КУЛЬТУРЫ

© 2010 М.А. Гутров

Челябинская государственная агроинженерная академия

Поступила в редакцию 27.03.2010

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

Ключевые слова: злаковая культура, свободные колебания

При анализе технологического процесса зерноуборочных жаток часто возникает вопрос о характере взаимодействия рабочих органов жатки с убираемыми растениями. С технологической точки зрения важно, чтобы рабочие органы жатки входили в контакт с растениями без удара. В частности, одним из критериев выбора режима работы мотовила (рис. 1) является отсутствие потерь свободным зерном, которое является результатом удара планки мотовила по растениям. Для исключения нежелательных технологических явлений в работе мотовила на практике выбирают такую траекторию его движения (рис. 2), при которой величина абсолютных скоростей точек планки не превышает допустимых значений [1], а векторы скорости направлены вдоль стеблей (вертикально вниз) в момент входа планки в стеблестой. Движение любой точки планки мотовила может быть задано системой двух алгебраических уравнений

y(t)= V • t + R

(

cos

- arcsin

V

R -a

\

+ a -1

— R -

1 —

V2

<t ) = l

+ -

V

(

— R - sin

a

— arcsin

V

R -a

R - a

+ a-1

где yOz - декартова система координат, начало которой совмещено с основанием растения; Ve -переносная скорость мотовила; R - радиус мотовила; шу - относительная угловая скорость мотовила; l - высота растения; t - время.

Величина и направление абсолютной скорости планки мотовила является главным критерием, на основе которого в механике сельскохозяйственных машин контакт планки с растением

Гутров Михаил Александрович, кандидат технических наук, доцент кафедры «Сопротивление материалов». E-mail: gma74@mail.ru

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

ось вращения,

планка мотовила

*

Рис. 1. Внешний вид мотовила зерноуборочной жатки

Рис. 2. Траектория движения планки мотовила и

схема взаимодействия с растениями: 1 - точка начала контакта планки и растений; 2 -точка завершения контакта

2

С позиций механики деформируемого твердого тела ударное взаимодействие обладает двумя признаками: кратковременность процесса и возникновение весьма больших ударных сил. Последний признак, по всей видимости, является следствием быстротечности процесса удара и возникающих при этом значительных ускорений в соударяющихся телах. Кратковременность действия значительных по величине сил инерции приводит к возникновению в телах различных волновых процессов. Кильчевский Н.А. в монографиях, посвященных теории соударения твердых тел [2], определяет термин «удар» как быстропротекающий механический процесс. Под продолжительностью быстропротекающего процесса понимается такой интервал времени, который значительно меньше периода собственных колебаний одного или нескольких взаимодействующих тел. Таким образом, величина периода первого тона колебаний разделяет ударное и неударное взаимодействие и служит основным критерием возникновения ударных явлений при контакте двух и более тел.

Применительно к механике взаимодействия планки мотовила и убираемых растений можно утверждать, что их контакт сопровождается ударом, если период собственных колебаний отдельного растения либо планки мотовила будет значительно больше времени удара. Следовательно, чтобы выбрать безударный режим работы мотовила необходимо знать 3 параметра: продолжительность времени контактного взаимодействия планки с растением, а также периоды низших круговых колебаний планки и отдельного растения. Период собственных колебаний планки или растения можно найти путем решения задачи о свободных колебаниях соответствующего тела. Из 2 рассматриваемых тел (планка и растение) необходимо выбрать то, которое имеет наибольший период собственных колебаний на первой частоте. Из теории колебаний следует, что более жесткому телу соответствует более высокие частоты собственных колебаний и, соответственно, меньшие значения периодов колебаний. Говоря о поперечных колебаниях тел, следует отметить, что конструкционная из-гибная жесткость поперечного сечения мотовила, а также жесткость его отдельных планок, безусловно, намного больше жесткости поперечного стебля на изгиб. Следовательно, для оценки технологического процесса контактного взаимодействия планки с растением на предмет возникновения ударных явлений необходимо в первую очередь сравнивать период собственных колебаний растения на первом тоне с продолжительностью времени его соударения планкой.

Для расчета собственных колебаний растения необходимо принять модель, которая бы адекватно описывала механику деформирования реального растения. В курсе лекций учебной

дисциплины «Уборочные машины», которые читаются для будущих инженеров агропромышленного комплекса, а также в специальной литературе посвященной исследованию технологического процесса, растение не рассматривается как деформируемое тело. В работе [5] показано, что на основе численного решения уравнений равновесия механики геометрически нелинейно деформируемых стержней можно адекватно решить задачу контактного взаимодействия растений с планкой мотовила зерноуборочной жатки. Отдельное растение моделировалось в виде стержня с распределенными массовыми и геометрическими параметрами, жестко заделанного одним торцом. Следовательно, для определения периода и частот собственных колебаний стержня необходимо решить соответствующие уравнения динамики гибких стержней. В первом приближении рассмотрим задачу свободных колебаний консольного стержня в линейной постановке, соответствующей малым перемещениям стержня. С этой целью используем модель растения в виде жестко закрепленного одним концом стержня с равномерно распределенными параметрами, а колос имитируем сосредоточенной массой на свободном торце стержня (рис. 3). Расчет собственных частот сводится к общеизвестному решению задачи о свободных изгиб-ных колебаний стержня - расчету дифференциального уравнения четвертого порядка в частных производных

д4 у т д2 у

—4 +--Т = 0

дz4 EJx д12

(1)

где Ш - жесткость поперечного сечения стебля; т - интенсивность распределенной массы стебля.

Рис. 3. Расчетная схема растения с учетом наличия колоса, т0 - масса колоса

Для рассматриваемой расчетной схемы уравнение (1) имеет решение в замкнутой форме [3]. Полагая, что уравнение упругой линии стержня можно представить в виде у( г, г ) = У (г)Т (г ), получают систему двух обыкновенных дифференциальных уравнений

Т (г) + а2Т (г )= 0

У (г У -а4У (г ) = 0

где а4=а>2т/Е/х; ю - частота собственных колебаний стержня.

Поиск решения второго уравнения системы ищем в виде суммы произведений постоянных интегрирования С, и функций Крылова Кг.

4

У (г) = ^ С Кг. С учетом краевых условий ин-

тегрирования

У(0)=0; У'(0)=0; У"(/)=0;

У'"(/)=-т0ю2У(Е/)-1 получаем систему четырех алгебраических уравнений относительно неизвестных С

У (0) = С1К1 + С2 К2 + С3К3 + С4 К4 = 0

^4 4

У '(0) = С1К4 + С2 К1 + С3 К2 + С4 К3 = 0 У "(I) = С1К3 + С2 К4 + С3 К1 + С4 К2 = 0

У"'(/) = С1К2 + С2 К3 + С3 К4 + С4 К1 = -

т0а У

а3 Е/ (2)

Первые два уравнения системы дают нулевые значения констант С1=С2=0. Поиск определителя системы, составленной из третьего и четвертого уравнения системы (2), приводит к частотному трансцендентному уравнению (3), расчет которого относительно неизвестного а возможен на основе графических методов, либо численных методов решения.

К,--- К4

1 1 К, 4

а К К 2 К

К 4 - К"К

= -к

Или

ек(а/ )ео8 (а1) +1

а 8Ъ(а/) еов(а/) - ек(а/) бш (а1)

= -к.

(3)

где кт=т0/т - коэффициент, выражающий отношение массы колоса т0 и распределенной массы стебля т.

Для установления факта наличия ударных явлений более всего нас интересует частота первого тона Ю\, которая влияет на величину периода свободных колебаний растения Т. Сравнивая величину Т1 и время взаимодействия стебля с

рабочим органом можно классифицировать их механическое взаимодействие. Приведенную на рис. 3 модель растения можно усовершенствовать. Введем на свободном торце стержня «толстый» стержень, который будет моделировать колос (рис. 4). Данное изменение в расчетной схеме позволяет учесть не только массу колоса, но и его геометрические параметры. /0 - длина колоса, /0=т0/20/12 - момент инерции колоса относительно оси, перпендикулярной плоскости рисунка и проходящей через центр масс колоса. Такая модель имеет больше физического сходства с реальным растением, чем модель, приведенная на рис. 3.

Рис. 4. Расчетная схема растения, учитывающая массовые и геометрические параметры колоса

Введенные в расчетную схему изменения требуют интегрирования уравнения (1) при других краевых условиях на свободном конце стержня. Краевые условия интегрирования уравнения (1) для расчетной схемы (рис. 4) для жесткой заделки будут равны. У(0)=0; У'(0)=0. На свободном конце стержня

У'(/) = /0а2У '(Е/ )-1, У'"(/) = - т0р 2У (Е/ )-1.

С учетом того, что К (0) ^ 0 ^ С1 = С2 = 0 система алгебраических уравнений для расчета постоянных интегрирования С, будет иметь вид.

С 3 К1 + С 4 К 2 =

/ 0а' а2 Е/

-У'

т0®

С 3 К4 + С 4 К] =---—

3 4 4 1 а 3 Е/

У

(4)

Расчет определителя системы (4), с учетом соотношений ю2=а4Е//т, /0=т0/20/12, У=С3К3+СК4, У=С3К2+С4К3, приводит к поиску корней частотного трансцендентного уравнения в виде.

1=1

1

к - к + ктоК4 к

2 К 4 + ктаКъ 1

к К1 + ктаК4 К

К,--К,

3 К4 + ктаКъ 2

= кк2ос2 — т 1 12

(5)

где kl=lJl - безразмерный коэффициент, выражает отношение длины колоса к длине стебля растения.

Для поиска корней уравнения (3) и (5) необходимо знать коэффициенты кг и кт, которые являются геометрическими параметрами растений злаковых колосовых культур. С учетом экспериментальных сведений [4, 5] примем отношение массы колоса и распределенной массы

стебля равным кт=0,861-1,205. Для середины диапазона имеем кт=1,03. Соотношение средних значений длины колоса и длины стебля для 5 сортов пшеницы составляет ^=0,089-0,117 [4]. Примем в качестве исходного параметра среднее значение кг=0,10. Для указанных средних значений коэффициентов кг и кт численно решим уравнения (3, 5) в системе МаШСАБ. Корни частотных уравнений будем искать для диапазона высот растений злаковых колосовых культур (пшеница, рожь, ячмень)0,5-1,5 м [16]. Первые 8 корней уравнения (3) для фиксированных значений высоты растений I приведены в таблице.

Таблица. Значения корней частотного уравнения (3)

1 ,м а а а3 а4 а5 а6 а7 а8

0,75 1,557 5,34 9,490 13,659 17,838 22,020 26,204 30,390

1,00 1,240 4,028 7,132 10,255 13,387 19,659 22,797 25,936

1,25 1,037 3,24 5,717 8,213 10,716 13,223 15,732 18,242

По расчетным данным, приведенным в таблице видно, что величина корней частотного уравнения а{ существенно зависит от длины стебля. Так как нас интересует только первая частота колебаний растения, то построим графики а\=а\(Г) значений корней частотных уравнений (3, 5) (рис. 5).

Рис. 5. Зависимость значения первого корня а1 частотного уравнения (3, 5) от высоты растения I

Варьирование геометрическими коэффициентами кг и кт в пределах вышеуказанного диапазона приводит к изменению значений функции а1=а1(1) в пределах 4%. Из графиков следует, что учет в расчете инерции поворота колоса несущественно приводит к уменьшению значения а! не более чем на 0,11%. Математически уравнение (5) значительно проще, чем уравнение (3). Следовательно в дальнейших расчетах можно пользоваться расчетной схемой растения, изображенной на рис. 3. По расчетным значениям корней а, уравнений (3, 5), вычислим ряд значений круговых частот стержня:

0 = а

Ю_ т

(6)

Период свободных колебаний стержня, вычислим по зависимости:

Т =

2п

0

2п

(7).

Для расчета частоты свободных колебаний необходимо знать жесткость поперечного сечения стебля. Изгибная жесткость поперечного сечения Б1 реального стебля является функцией его длины I, однако изменение жесткости в пределах длины одного растения незначительно. По экспериментальным сведениям [4] межвидовое и межсортовое различие величины модуля упругости более существенно.

Примем величину изгибной жесткости поперечного сечения стебля пшеницы Б1=0,022-0,045 Нм2. [5]. Для указанного диапазона жесткости по функции а!=а!(1) и формулам (6, 7) построим график Т1=Т1(1) (рис. 6). Приведенные графики являются предельными кривыми. Они охватывают весь существующий диапазон жесткости поперечного сечения стебля и высот растений злаковый колосовых культур. На практике для конкретного сорта и региона возделывания злаковых культур диапазон высот может быть более узким [15]. Для практического использования аппроксимируем результаты численного расчета (рис. 5) аналитической функцией в виде полинома третьей степени

2

а(1 ) = Х ЬгГ

(8)

Рис. 6. Зависимость периода собственных колебаний на первой частоте от высоты растения I и жесткости EJ поперечного сечения стебля

Коэффициенты Ьг полинома (8) для предельных значений диапазона изгибных жесткостей поперечного сечения стебля равны

Ь= (4,39; -6,42; 4,35; -1,08)

(9)

Максимальная относительная погрешность аппроксимации не превышает 1,01%. С учетом (8) период свободных колебаний (7) запишем в виде:

7 =■

2п

X Ь1г л —

(10)

Взаимодействие рабочего органа жатки и растения будет носить ударный характер, если время контакта будет значительно меньше, чем период собственных колебаний растения на первой частоте. Следовательно, если растение злаковой колосовой культуры взаимодействует с рабочим органом менее чем Т1=0,27-2,26 с, то такое взаимодействие можно классифицировать как удар. Однако для этого необходимо знать реальную продолжительность времени взаимодействия растения с рабочим органом жатки. В первом приближении продолжительность времени механического контакта планки с растением рассчитаем из уравнений движения точек планки мотовила. Согласно технологическим требованиям к работе мотовила, абсолютная скорость ее планки в момент времени ^ начала контакта с растениями не должна иметь горизонтальную составляющую, т.е. Уy=dy/dt=0. Время г2 контакта растения с планкой обусловлено моментом среза растения режущим аппаратом.

Положим для момента времени г2 вертикальная составляющая скорости планки отсутствует Vz=dz/dt=0. В этом случае продолжительность времени контакта, будет равна

1

(

Аг = г 2 - = —

®г V

п

--агсБт

2

Л

(11)

где Х=Яюу/Уе - безразмерный коэффициент, Х=1,1-1,8 [1] По графику (рис. 7) можно установить, что продолжительность времени контакта планки с растениями существенно зависит от рабочей скорости движения жатки и лежит в диапазоне Дг=0,06-0,28 с.

Период собственных колебаний растения на первом тоне 7^=0,27-2.26 с на порядок превышает диапазон времени контакта планки и растений, что позволяет утверждать об их ударном взаимодействии. Запишем равенство Дг=Т\. С учетом формул (10, 11) и Х=Яюу/Уе имеем:

Я

(

е V

п

-- агс81п

2

Л

2п

X Ь1

.(12)

Рис. 7. Продолжительность времени контакта планки с растениями при различных значениях коэффициента X

Из формулы (12) выразим скорость Уе движения жатки по полю, получим

У =

■ - агсзт

X Ь1г

(13)

Практическая значимость полученной функции (13) заключается в том, что она устанавливает взаимозависимость между геометрическими параметрами, физико-механическими свойствами убираемых растений и скоростью движения зерноуборочной жатки (рис. 8). Из рисунка следует, что график функции Уе=Уе(Х) содержит экстремум. Координату точки Хэ

г=0

г=0

г=0

г=0

экстремума функции (13) с учетом того, что Я, тс, I, Е/, Аэ не равны нулю, вычислим из уравнения:

dV„

= arcsin

_1_

X

+ -

1

X 1 -

= -* = 0 1 2 X?

(14)

Уравнение (14) является трансцендентным. В диапазоне значений технологического коэффициента А=1,1-1,8 уравнение имеет единственный корень Аэ=1,533. На практике это значение соответствует предельной рабочей скорости Утах жатки и ее максимальной производительности при условии отсутствия ударных явлений в контактном взаимодействии планок мотовила и убираемых растений.

Рис. 8. Графики предельных значений рабочих скоростей зерноуборочной жатки при различных значениях коэффициента А

Согласно приведенным графикам (рис. 8), максимальная рабочая скорость движения жатки по полю не должна превышать Vmax=0,3-3 км/ч. На практике технологическая скорость движения современных жаток может достигать 8-12 км/ч [1], что приводит к потерям зерна.

Выводы: полученные математические выражения (8-10, 13) позволяют классифицировать механическое взаимодействие планки и растения, а также разработать практические рекомендации к выбору рабочей скорости движения жатки с целью минимизации потерь зерна при взаимодействии растения и мотовила.

СПИСОК ЛИТЕРАТУРЫ:

1. Иванцов, В.И. Валковые жатки / В.И. Иванцов, О.И. Солошенко. - М.: Машиностроение, 1984. -200 с., ил.

2. Кильчевский, Н.А. Динамическое контактное сжатие твердых тел. Удар. - К.: Наукова думка, 1976. - 311 с.

3. Бидерман, В.Л. Теория механических колебаний. - М: Высшая школа, 1980. - 408 с.

4. Лукьянова, И.В. Анализ видовых и сортовых особенностей устойчивости стеблей злаковых культур к полеганию с учетом их физико-механических свойств и архитектоники для использования в селекции. Автореф. на соис. уч. ст. докт. биол. наук. ФГОУ ВПО «Кубанский государственный аграрный университет»., Краснодар, 2008. - 51 с.

5. Гутров, М.А. Определение параметров зерноуборочных жаток на основе математического моделирования технологического процесса. Дисс. на соис. уч. ст. канд. техн. наук. ФГОУ ВПО «Челябинский государственный агроинженерный университет»., Челябинск, 1999. - 211 с.

CALCULATION THE FIRST TONE PERIOD OF GRAIN CROP

OWN OSCILLATIONS

© 2010 M.A.Gutrov Chelyabinsk State Agroengineering Academy

Calculation of free oscillations period of grain crop free oscillations of eared grass on the lowest frequency is lead. Are drawn outputs on character of technological interaction of a reel and removed plants.

Key words: grain crop, free oscillations

Mikhail Gutrov, Candidate of Technical Sciences, Associate Professor at the Department "Resistance of Materials". E-mail: gma74@mail.ru