Дифференциальное уравнение двигателя

Функциональная схема собственно двигателя как элемента ком­бинированного двигателя представлена на рис. 19, а. Условие статического равновесия определяется уравнением (1).

При нарушении установившегося режима крутящий момент двигателя М и момент сопротивления Мс получают приращения ?М и ?МС, в общем случае не равные между собой, поэтому

М + ?М ? Мс + ?Мc.

Избыток (или недостаток) энергии, появляющийся в связи с этим в системе двигатель—потребитель, расходуется на измене­ние скорости движения подвижных деталей этой системы. Спо­собность двигателя аккумулировать в себе механическую энергию определяется массой подвижных деталей двигателя и потребителя, т. е. их суммарным моментом инерции J, приведенным к оси вра­щения коленчатого вала двигателя.

Приведенный момент инерции двигателя и агрегатов, связан­ных с ним в работе (потребителей), можно определить в виде суммы

J=Jд + Jп ,

где Jд — приведенный момент инерции масс подвижных деталей двигателя; Jп — приведенный момент инерции масс подвижных деталей потребителя.

Из составных частей двигателя наибольшая доля приведен­ного момента инерции принадлежит маховику. Практика показы­вает, что с достаточной степенью точности в расчетах можно при­нять Jл ? (1,1 ?1,4) Jм.

Схема определения приведенного момента инерции

В зависимости от конструкции движущихся деталей приведен­ный момент инерции потребителя

где Je — момент инерции одной из вращающихся масс; nе — число вращающихся масс; ие — передаточное число передачи, связыва­ющей двигатель с одной из вращающихся масс.

При работе двигателя на гребной винт (судовые или авиацион­ные установки) находят приведенные моменты инерции масс греб­ного вала, передач и гребного винта.

Для выполнения расчета необходимо определить суммарную площадь сечений гребного винта на каком-либо произвольном радиусе (рис. 31).

Если известна площадь fi, то можно определить ширину ци­линдрической поверхности того же радиуса ri по соотношению

Взяв несколько таких радиусов, можно построить профиль диска, момент инерции которого приблизительно равен искомому. Однако профиль диска, построенного таким способом, криволи­нейный, поэтому для удобства расчетов сложный диск можно заменить набором k соответствующих цилиндрических колец (не менее 15). Тогда приведенный момент инерции гребного винта

Влияние воды на приведенный момент инерции погруженного в нее гребного винта на основании опытных данных оценивается увеличением приведенного момента инерции сухого винта на 25—40%. Приведенный момент инерции гребного винта при по­гружении в воду увеличивается тем больше, чем больше угловая скорость вала машины.

Приведенный момент инерции Jп потребителя в виде транспортного агрегата (автомобиль, трактор) определяется из условий равенства кинетической энергии приведенной массы

где ? — угловая скорость коленчатого вала; Ма — масса агре­гата; ? — скорость его движения; Твр — кинетическая энергия различных вращающихся масс в механизме транспортного агре­гата (передачи, валы, фрикционы и т. п.).

Величину Твр обычно оценивают приближенно относительно кинетической энергии основной массы Ма, поэтому

Суммирование полученных таким образом приведенных моментов инерции элементов силовой установки дает возможность опреде­лить приведенный к оси вращения коленчатого вала момент инер­ции J силовой установки двигателя в целом.

Если воспользоваться принципом д'Аламбера, то уравнение динамического равновесия системы двигатель—потребитель можно представить в виде уравнения (8).

Так как при неустановившихся режимах работы угловая ско­рость коленчатого вала, крутящий момент двигателя и момент сопротивления получают приращения соответственно ??, ?М и ?МС, то

где ?0, М0 и Мс0 — соответственно угловая скорость коленчатого вала, крутящий момент двигателя и момент сопротивления потре­бителя на выбранном равновесном режиме.

Подстановка выражений (15) в уравнение (8) приводит последнее к виду

Если учесть в полученном выражении уравнение статического равновесия (1), то раскрытия скобок

Статические характеристики потребителей (см. рис. 25) пока­зывают, что момент сопротивления Мс потребителя в общем слу­чае определяется скоростным режимом силовой установки — угловой скоростью со коленчатого вала и настройкой самого по­требителя — выбором его определенной статической характери­стики. В зависимости от типа потребителя под его настройкой по­нимают выбор либо угла атаки винта (судовые и авиационные условия), либо определенной передачи, или определенного про­филя и качества дороги в транспортных условиях и т. д. Если параметр, характеризующий настройку потребителя, обозначить ./V, то момент сопротивления потребителя определится функцио­нальной зависимостью

Характеристики на рис. 25 показывают, что зависимость (17) в общем случае является нелинейной. Поэтому для определения приращения ?Мс функцию (17) необходимо разложить в ряд Тейлора. В этом случае

При малых значениях отклонений ?? и ?N разложение можно линеаризовать, в результате чего приращение момента сопроти­вления ?Мс определится простым выражением вида

Для определения производных, входящих в выражение (18), можно воспользоваться характеристиками, например, приведенными на рис. 25. Производная дМс / д? при N = const определится непосред­ственно в виде тангенса угла наклона касательной к точке харак­теристики, соответствующей выбранному равновесному режиму. Для получения производной дМс / дN при ? = const характери­стики на рис. 25 необходимо перестроить в координаты Мс = = f (N) При ? = const.

Анализ параметров, влияющих на крутящий момент двигателя, показывает, что число этих параметров, номенклатура и удельный вес каждого из них зависят от типа двигателя. Так, крутящий момент дизеля связан с его эффективной мощностью Nе соотно­шением

M = 103Nc /?.

Мощность может быть определена в виде отношении часового расхода топлива Gч к удельному эффективному расходу топлива ge:

здесь Ни — теплотворная способность топлива; iд — число ци­линдров в двигателе; gц — цикловая подача топлива; ?е — эффек­тивный КПД; ?д — тактность двигателя.

Подстановка полученного выражения для подсчета эффектив­ной мощности в исходное выражение крутящего момента приводит последнее к виду

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

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

Изложенное свидетельствует о необходимости индивидуаль­ного анализа условий, вызывающих изменение крутящего момента применительно к двигателю определенного типа.

В качестве примеров рассмотрим дизель х наддувом и дизель без наддува.

Дизель с наддувом. Крутящий момент дизеля соотноше­нием (19) связан с цикловой подачей топлива и эффективным КПД, поэтому

Разложение этой функциональной зависимости в ряд Тейлора и последующая линеаризация дают

Частные производные, входящие в полученное выражение, в соответствии с формулой (19) имеют вид

Эффективный коэффициент полезного действия определяется произведением индикаторного ?i и механического ?м коэффициен­тов полезного действия двигателя:

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

Эксперименты показывают, что изменение скоростного режима дизеля при постоянном положении органа управления (внешняя и частичные скоростные характеристики) сопровождается незначи­тельным изменением ?е. Например, для дизеля 1Д6Н это измене­ние не превышает 1,5—2,0% в широком диапазоне скоростных ре­жимов.

При работе дизеля по нагрузочной характеристике (? = const) механический КПД при наличии наддува изменяется менее суще­ственно в области нагрузок, близких к номинальной. При средних и малых нагрузках влияние ?м на величину ?е более заметно. Однако во всем диапазоне нагрузочных режимов определяющее влияние на ?e оказывает индикаторный КПД. Основным фактором, определяющим значение ?i при изменении нагрузки в дизелях с наддувом, является коэффициент избытка воздуха ?. Поэтому для всех возможных установившихся режимов работы эффективный КПД дизеля с наддувом определяется двумя основными пара, метрами: ? — коэффициентом избытка воздуха и ? — угловой скоростью коленчатого вала.

В со­ответствии с этим для эффектив­ного КПД определяющей яв­ляется зависимость

показанная для одного из скоростных режимов на рис. 32. Разложение функции (22) в ряд Тейлора и последующая линеари­зация дают

Производные, входящие в это выражение, должны определяться по точке статической характеристики, соответствующей выбран­ному равновесному режиму. Как известно,

где l0 — количество воздуха, теоретически необходимое для сго­рания 1 кг топлива; Gвц — цикловая подача воздуха, определяемая выражением

Здесь Vh — рабочий объем цилиндра двигателя; pк — плотность воздуха перед впускными клапанами; ?? — коэффициент напол­нения цилиндра. Следовательно,

Параметры Vh и l0 являются постоянными для всех режимов работы двигателей, поэтому

После разложения полученной зависимости в ряд и последующей линеаризации с учетом формулы (24) можно получить

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

pк / pnк = const

позволяет установить, что

Коэффициент наполнения в основном зависит от угловой ско­рости ? коленчатого вала и давления наддува рк, поэтому

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

Подстановка соотношений (25)-(27) в выражение (23) приводит последнее к виду

Функциональные зависимости (20) и (28) показывают, что кру­тящий момент дизеля с наддувом определяется тремя основными параметрами: цикловой подачей топлива gц, давлением наддува рк и угловой скоростью коленчатого вала ?, т. е.

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

Производные, входящие в это разложение, в соответствии с формулами (20) и (28) определяются выражениями

Следует при этом подчеркнуть, что производная дМ / д? в выраже­ниях (30) и (31) учитывает изменение крутящего момента при изме­нении скоростного режима только за счет изменения эффективного КПД и коэффициента наполнения при неизменной цикловой подаче топлива. Чтобы подчеркнуть это обстоятельство, производная дМ / д? далее дополняется индексами ?е, ??:

Таким образом, определены приращения момента сопротивле­ния (18) и крутящего момента двигателя (30). Подставляя их в ис­ходное уравнение (16), последнее можно представить в виде

Анализ переходных процессов и сопоставление их между собой, оценка динамических свойств двигателя оказываются наиболее удобными, если в уравнение (32) вместо абсолютных значений координат ввести их относительные безразмерные значения (коор­динаты). Для осуществления такой замены следует выбрать базо­вые значения соответствующих параметров. Выбор этот может быть произвольным, однако наиболее часто в качестве базовых значений используют значения параметров выбранного равновес­ного режима, от которого отсчитываются отклонения; значения параметров номинального режима и некоторые другие.

Если базовыми принять значения координат на выбранном рав­новесном режиме, то в соответствии с этим относительные коорди­наты, определяемые уравнением (32), получат вид

здесь ? — безразмерное изменение угловой скорости коленчатого вала; q — безразмерное изменение цикловой подачи топлива; ? — безразмерное изменение давления наддува; ?д — безразмерное изменение настройки потребителя (изменение нагрузки двигателя).

Если ввести безразмерные координаты (33) в уравнение (32), то последнее после деления всех членов уравнения на коэффициент при q получит вид

Коэффициент Т данного уравнения, имеющий размерность - вре­мени, называют временем собственно дизеля:

безразмерный коэффициент к называют коэффициентом самовы­равнивания двигателя:

их называют коэффициентами усиления по давлению наддува 01 и настройке потребителя 02.

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

При работе дизеля на холостом ходу равновесный режим возмо­жен при равенстве индикаторного крутящего момента двигателя Мi моменту сил внутреннего сопротивления Мвт, т. е. при выпол­нении условия

M1 — Мвт = 0.

При нарушении равновесного режима холостой работы двигателя уравнение динамического равновесия (16) следует представить в форме

где ?Мi — отклонение индикаторного крутящего момента двига­теля от его значения на выбранном равновесном режиме; ?Mвт — отклонение момента внутренних сил сопротивления от его значе­ния на выбранном равновесном режиме. Индикаторный крутящий момент двигателя, так же как и эффективный, зависит от цикловой подачи топлива gц и почти не зависит от давления наддува в связи с небольшими дозами впрыскиваемого топлива. Момент внутрен­них сопротивлений при данном тепловом состоянии двигателя определяется только угловой скоростью коленчатого вала. Следова­тельно,

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

Дизель без наддува. Система воздухоснабжения дизеля без наддува на всех возможных режимах его работы предусматривает такое поступление воздуха в цилиндры, при котором происходит полное сгорание топлива. Если весь тракт подачи воздуха к ци­линдрам не имеет сопротивлений, которые могут отрицательно сказаться на коэффициенте наполнения при нарушении установив­шихся режимов, то коэффициент избытка воздуха такого дизеля на всех установившихся и неустановившихся режимах поддержи­вается на уровне, исключающем заметное изменение эффективного КПД (например, правее значения ? в точке D на рис. 32). В этом случае с достаточной степенью точности можно принять условие

существенно упрощающее формулы производных (31).

Если дополнительно учесть, что эффективный КПД слабо зави­сит от скоростного режима, и принять д?e / д? ? 0, то коэффициенты (35)—(39) дифференциального уравнения (34) будут иметь вид

В связи с этим уравнение дизеля без наддува получит вид

Линейное дифференциальное уравнение собственно дизеля, как и всех других элементов системы, целесообразно представлять в операторной форме. В этом случае операция дифференцирования по времени d / dt обозначается символом р. В соответствии с этим приемом вторая производная по времени обозначается символом p2 и т.д.:

Операцию интегрирования также можно выразить через символ в виде отношений

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

называется собственным оператором рассматриваемого элемента.

С учетом обозначения (46) запись дифференциальных уравнений существенно упрощается: