CC BY
124
Заключение
В данной работе показан подход к разработке процедуры автоматизированной многопараметрической оптимизации управления двигателями внутреннего сгорания (ДВС), работающих на квазистационарных режимах в составе транспортных средств, с целью снижения токсичности и улучшения топливной экономичности. Оптимальные режимы работы ДВС характеризуются парами значений крутящего момента и частоты вращения, определенные для каждого характерного значения мощности в результате оптимизации.
В качестве критерия оптимальности был выбран интегральный показатель, объединяющий в себе критерии оптимальности по токсическим и экономическим показателям работы ДВС.
Были представлены результаты оптимизации параметров управления ДВС для серии экспериментальных данных, которые показали высокую скорость и надежность выбранного пути и алгоритма.
Следующие работы в рамках данного проекта будут посвящены дальнейшей разработке программно-аппаратного комплекса оптимизации параметров управления ДВС.
Литература
1. "GMMC-Combustion Engine Optimization Project," 2011. [Online]. Available: http://www.fcc.chalmers.se/comp/projects/combustion-engine-optimisation.
2. FRONTIER Co, "GT-SUITE User's Conference," in Optimization of Automotive Control Parameters with FRONTIER, 2001.
3. Гладков Л.А. Генетические алгоритмы: Учебное пособие, 2-е издание ed., Москва: Физматлит, 2006.
4. Мэтьюз Д. Численные методы. Использование MATLAB., Вильямс., 2001.
5. Buhmann M.D. Radial Basis Functions: Theory and Implementations, Cambridge University, 2003.
6. Смирнов А.Б., Петриченко Д.А., Шустров Ф.А., Татарников А.П. "Оптимизация параметров управления ДВС на статических режимах работы на основе многопараметрической нейросетевой модели" Современные проблемы науки и образования - 2011.-№6.(приложение "Технические науки"). (Электронный журнал) URL: http://online.rae.ru/880 (дата обращения: 19.01.2012 г.).
7. Хрипач Н.А., Лежнев Л.Ю., Папкин Б.А., Тингаев Н.В., Шустров Ф.А., Татарников А.П. " Механотронные системы управления газообменом и сгоранием " Современные проблемы науки и образования - 2011. -№6. (приложение "Технические науки"). (Электронный журнал) URL: http://online.rae.ru/889 (дата обращения: 19.01.2012 г.).
8. Ипатов А.А., Кутенёв В.Ф., Каменев В.Ф., Фомин В.М., Хрипач Н.А., Лежнев Л.Ю., Артёмов А. А. " Разработка автомобилей с гибридной силовой установкой, работающей на водородных видах топлива " Труды НАМИ / ГНЦ РФ ФГУП «НАМИ». - Вып. №242: Комбинированные энергоустановки автотранспортных средств: сб. научн. ст. - М., 2009. -с. 26 - 66.
9. Ипатов А.А., Хрипач Н.А., Лежнев Л.Ю., Кириллов В.А., Папкин Б.А. "Автономные системы выработки тепловой и электрической энергии на биотопливе" Энергия: экономика, техника, экология. - 2010. - № 3. - с. 6-12.
Математическое моделирование гидропривода вентилятора для системы
охлаждения автомобильного двигателя
Труханов К. А.
МГТУ им. Н.Э.Баумана, trukhanov@mail.ru
Аннотация. В статье приведена математическая модель гидропривода вентилятора для охлаждения автомобильного двигателя, уравнения, блок схемы соответствующих элементов гидросистемы. В качестве модели потока при описании про-
цесса теплообмена принята модель идеального перемешивания. Получена математическая модель гидропривода и соответствующие ей переходные процессы.
Ключевые слова: гидропривод вентилятора, охлаждение двигателя, математическое моделирование, математическая модель гидропривода вентилятора для системы охлаждения автомобильного двигателя.
Современные методы проектирования любой гидросистемы основаны на подробных знаниях о процессах, протекающих в рассматриваемых системах. Такими знаниями служат математическая модель и описание системы. Моделирование системы особенно важно - это быстрее и дешевле, чем проводить физический эксперимент. Полученные результаты в ходе математического моделирования могут быть использованы для оптимизации системы, определение режимов эксплуатации, причин возможных отказов.
Принятой при анализе стала схема гидропривода вентилятора для системы охлаждения автомобильного двигателя с дискретным управлением гидрораспределителя [1]. Описание гидросистемы проводилось в виде уравнений и блок-схем отдельных элементов, входящих в ее математическую модель.
Уравнение, определяющее подачу насоса:
=ар ■Кчр -цур, (1)
где: ар - угол наклона диска, либо блока насоса, Кдр - коэффициент подачи насоса, Пур - объемный КПД насоса.
Предполагается, что ар и Пур - входные сигналы, а 0Р - выходной сигнал, соответственно.
Блок-схема уравнения (1) приведена на рисунке 1.
ар Ор
-
Рисунок 1 - Блок-схема уравнения 1
Уравнение переливного клапана, для упрощения математической модели которого, принята статическая характеристика:
Ор = кzb ■(р - Рь X при р > р^ (2)
Qz = 0, прир <рь, (3)
в уравнении (2) и (3): ^ - расход через клапан, Кь - коэффициент крутизны статической характеристики клапана, р - рабочее давление системы, рь - давление срабатывания клапана. Блок-схема уравнения (2), (3) приведена на рисунке 2.
Р
Рь
Рисунок 2 - Блок-схема переливного клапана
Уравнение, описывающее сжимаемость жидкости, можно определить, используя приближенное соотношение:
а=^ • Ар ,
Е, А
(4)
где: Он - расход жидкости, вызванный сжимаемостью жидкости на участке насос-
распределитель, Уц - объем жидкости в трубопроводах и элементах гидросистемы на участке насос-распределитель, Ец - модуль объемной жидкости, р - давление жидкости на участке насос-распределитель (приняли, что трубопроводы короткие, поэтому пренебрегаем инерцией жидкости и потерями давления в них из-за сопротивления трения [2]), ШШ - оператор дифференцирования, применяемый ниже во всех уравнениях.
Соотношение (4) можно использовать для определения величины давления:
Е г
Р = ~-Т-I Qs ■ Ж, (5)
Ks
Блок-схема уравнения (5) приведена на рисунке 3.
Рисунок 3 - Блок-схема уравнения 5
Уравнение гидромотора:
Qm
К цт ' ® т
П ут
(6)
в уравнении (6): Qm - расход жидкости, идущей к гидромотору, Кдт - рабочий объем гидромотора (далее ГМ), приведенный к одному радиану угла поворота вала (характерный объем), юш - угловая скорость вращения вала ГМ, Пуш - объемный КПД насоса.
Крутящий момент, развиваемый гидромотором, может быть описан уравнением:
-Мт К цт ' ^рт ' П тт, (7)
где: Мт - крутящий момент на валу ГМ, Кдт - характерный объем ГМ, Арт - перепад давления в ГМ, Арт =ра-ръ, пшш - механический КПД ГМ. Блок-схемы уравнений (6) и (7) приведены на рисунках 4 и 5, соответственно.
К> 111 \ „
*
л™ —
(У1
Рисунок 4 - Блок-схема уравнения 6 Рисунок 5 - Блок-схема уравнения 7
Уравнение нагрузки на валу ГМ:
Крутящий момент Мт, создаваемый на валу ГМ, равен сумме моментов от действующих сил.
Уравнение моментов имеет вид:
Мт = М1 + Мв + Мь,, (8)
где: М] - момент инерционных нагрузок, Мв - момент трения, Мь - момент, создаваемый вентилятором.
После подстановки значений, указанных моментов, уравнение (8) принимает вид:
г2
Мь = Ра ■А/ -Яы -Га',,
¥а = (0,06....0,08) - 3 Оа 2 ,,
т - (Я2 + г2)
]г
2
(9)
(10)
(11)
Мт =
Л ш,
Л
+ Вт ■шт + МТ
(12)
в уравнениях (9) - (12): ра - плотность воздуха при средней его температуре в радиаторе (см. формулу 40), Лг - фронтальная поверхность вентилятора, Яы - внешний радиус лопастей вентилятора, г - внутренний радиус лопастей вентилятора, т - масса вентиляторного колеса, Уа - осевая скорость на выходе вентилятора, 1т - момент инерции двигателя и элементов конструкции нагрузки (вентилятора), Вт - коэффициент сопротивления, ют - угловая скорость вращения вала ГМ (принимаем, что скорость вращения вала ГМ равна скорости вращения вентилятора), Юг - угловая скорость вращения вала вентилятора.
Соотношение (12) используется для определения угловой скорости вала ГМ.
Блок-схема уравнения (12) дана на рисунке 6.
Потоком рабочей жидкости, подводимой к ГМ вентилятора, управляет дискретный распределитель при наличии внешнего воздействия (сигнала). Таким сигналом в рассматриваемом случае является электрический сигнал (24 В), подаваемый на катушку устройства (электромагнит, далее ЭМ). ЭМ состоит из (рисунок 8а) цилиндрического корпуса-магнитопровода 1, в котором соосно расположена обмотка 3. Внутри расположен цилиндрический якорь 2. Магнитный поток проходит через два зазора: рабочий плоский 4 и паразитный радиальный 5 [3].
Математическая модель дискретного распределения представлена следующими уравнениями [4]:
ХР = М ■ а = М ■
ЛУ
--М ■
Л 2 х
Лх
Ре - К5 ■ (х - хо) - В ~ = М
ли2 '
Л 2 х
Л 2 х Л2
М
Ре - Кь
V
лх Л
(х - х0) - В —
1
Ре = - ■ г 2 ■ ф(х),
(13)
(14)
(15)
(16)
в уравнениях (13) - (16): х - линейное перемещение золотника, х0 - предварительное смещение золотника, создаваемое пружиной, М - масса золотника, а - ускорение золотника, У -скорость золотника, К - коэффициент жесткости пружины, В - коэффициент трения золотника электромагнита, Ре - электромагнитная сила, г - ток контура управления, Ф(х) - функция, определяемая уравнением (23).
Блок-схема уравнения (15) приведена на рисунке 7.
М
О"
в
Рисунок 6 - Блок-схема уравнения 12
Рисунок 7 - Блок-схема уравнения 15
Рисунок 8 - Электромагнит с плоским стопом [3]: а) конструктивная схема, б) зависимость развиваемой силы от перемещения якоря. 1 - корпус-магнитопровод; 2 - якорь; 3 - обмотка; 4 - осевой зазор; 5 - радиальный зазор.
Электрическая часть подсистемы ЭМ описана следующими уравнениями:
и 0 = Я ■ г +
Л X
Л
X - Ь( х) ■ г, и0 = Я ■ I + Ь(х)— +1 ■ Ф(х) ■—,
Ь( х)
N:
ц 0 ■ п ■ а ■ Л ■ N
2 ' х ^
Я
м
V а + х J
х
V а + х J
£ =
Ям -■
ц0а■Л■N
&
ё
га + хЛ
Ф( х) -1
ц0 ■ п ■а■Л
Л (Ь( х))
Лх
V х (а + хУ
Лх
и0 - Я ■ г -г ■ Ф(х) ■ — 0 Лг
(17)
(18)
(19)
(20) (21) (22)
(23)
(24)
Лг ь(х)
в уравнениях (17) - (24): г - ток контура управления, х - линейное перемещение золотника, Я - сопротивление катушки электромагнита, Ям - сопротивление потерь контура электромагнитной катушки (магнитное сопротивление), и0 - напряжение питания катушки электромагнита, цо - магнитная постоянная, а - величина ферромагнитного шунта (осевое расстояние), ё - величина радиального зазора, Л - диаметр штока, Ь(х) - индуктивность катушки электромагнита, N - число витков электромагнитной катушки, X - магнитный поток.
Рисунок 9 - Блок-схема уравнения 24
Ку = • бШа-
(26)
Раздел 1. Наземные транспортные средства, энергетические установки и двигатели. Блок-схема уравнения (24) дана на рисунке 9.
Уравнение расхода жидкости через распределитель в линии распределитель-гидромотор
= КУ-Х-^р - ра , (25)
12 ь,
в уравнениях (25), (26), К - удельная проводимость золотника, х - линейное перемещение золотника, - коэффициент расхода, р - плотность рабочей жидкости, р - давление в трубопроводе перед распределителем, ра - давление в трубопроводе перед ГМ, С - диаметр золотника, а - угол конусности запорного элемента.
Уравнения баланса расхода жидкости, поступающей к гидромотору:
Qv1 - Qm = 0б1, (27)
<28)
си
где: - расход жидкости, связанный с сжимаемостью жидкости на участке распределитель-гидромотор, У§1 - объем жидкости в трубопроводах и элементах гидросистемы на участке распределитель-гидромотор, Б8 - модули объемной жидкости, ра - давление жидкости на участке распределитель-гидромотор.
Е г
Ра =~ГГ~ - | Qsl -с, (29)
Ksl *
Уравнение расхода жидкости, пропускаемой распределителем в линию гидромотор-распределитель
QV2 = Ку-х-^рь - р1, (30)
где: Ку - удельная проводимость золотника, см. уравнение (26), х - линейное перемещение золотника, рь - давление в линии за гидромотором, рг - давление на сливе (подпор в баке).
Блок-схемы уравнений (25) и (30) представлены на рисунках 10, 11, соответственно.
Qsl =
к
р , Рь ^ + г--.
Ра ( —► 11Л1 - (г * Р1 „ \и\ (Г
Рисунок 10 - Блок-схема уравнения 25 Рисунок 11 - Блок-схема уравнения 30
Уравнение баланса расхода жидкости, идущей от гидромотора:
Qm - Qv2 = QS2, (31)
Qs 2 = ^ х сь., Ес с
(32)
в уравнениях (31), (32), Qs2 - расход жидкости, связанный с сжимаемостью жидкости на участке гидромотор-распределитель, У§2 - объем жидкости в трубопроводах и элементах гидросистемы на участке гидромотор-распределитель, Б8 - модули объемной жидкости, рь - давление жидкости на участке гидромотор-распределитель.
Е г
рь =77^ \ Qs2 -С, (33)
К
5 2
Блок-схемы уравнений (29) и (33) представлены на рисунках 12 и 13, соответственно.
Рисунок 12 - Блок-схема уравнения 29 Рисунок 13 - Блок-схема уравнения 33
Гидромотор крепится на специальном фланце с боковой стороны радиатора, вращает вентилятор, прогоняющий воздух через радиатор. Расход воздуха, проходящий через радиатор, зависит от скорости вращения вентилятора, ю^
Математическое описание процесса в теплообменных аппаратах удобно записать в виде выражения, характеризующее изменение температуры в потоке теплоносителя во времени, обусловленное, во-первых, движением потока и, во-вторых, теплопередачей [5, 6].
Для рассматриваемой системы принято, что двигатель автомобиля работает с максимальной мощностью, на максимальных оборотах. Кроме того, предполагаем, что радиатор будет установлен за кабиной, поэтому потоком набегающего воздуха можно пренебречь.
Принято также, что структура потока теплоносителя соответствует модели идеального перемешивания.
Схема контура системы охлаждения автомобильного двигателя показана на рисунке 14:
Рисунок 14 - Схема контура системы охлаждения автомобильного двигателя: ДВС - двигатель внутреннего сгорания (далее ДВС);
Трубопровод 1 и Трубопровод 2, соединяющие радиатор и ДВС.
Схемы к уравнениям теплового баланса для радиатора и ДВС представлены на рисунках 15 и 16, соответственно.
Уравнение баланса энергии радиатора, исходя из баланса количество тепла, отводимого от двигателя и передаваемого от воды к охлаждающему воздуху, можно применить уравнения:
Qa = Qc, (34)
ma ' Ca ' (Ta,e - Ta,i ) = Mc ' Cc ' (Tc,i - Tc,e ), (35)
Mc ■ cc
dTc,, dt
Mc ' Cc ' (Tc,i Tc,e ) - £ ' Ma ' ca ' (Tc,i Taj ),
(36)
—Ты,тс
!
„ 1 Радиатор, „ 1
* 1 тг * 1
гг
Рисунок 15 - Схема к уравнению теплового баланса для радиатора.
Рисунок 16 - Схема к уравнению теплового баланса для ДВС.
Количество тепла, отводимое от радиатора охлаждающим воздухом:
йа = Оа -р а • С с
&Та
О а = %• (К^ -г 2) • • Ьы • цг • фт(аы )• соБ(а ы) • Юf,
Оа=кг•ю,
(37)
(38)
(39)
'а / '
в уравнениях (34) - (39): Оа - производительность вентилятора [7], КЬ1 - внешний радиус лопастей вентилятора, г - внутренний радиус лопастей вентилятора, - число лопастей вентилятора, ЬЬ1 - ширина лопастей по хорде, п - КПД вентилятора, аы - угол наклона лопасти к направлению воздушного потока, юг = юш - угловая скорость вращения вала вентилятора, kf - коэффициент пропорциональности производительности вентилятора, ра - плотность воздуха при средней его температуре в радиаторе, найденная по уравнению
р0 0,1-106 1 Л кг Ра = „ = _ _ = 1,07 — , (40)
К • т
* * ж
287•325
м
здесь: йС - количество отводимого тепла теплоносителем, ЕИ_ех - поверхность охлаждения радиатора, ТС - температура теплоносителя в радиаторе, Та - температура воздуха, проходящая через радиатор.
КИ_ех - коэффициент теплопередачи, определяемый с помощью формулы К =■ 1
Л И ех
1
+
8ь
+
1
(41)
а с /А,1 /а а
в уравнении (41): ас - коэффициент теплоотдачи от теплоносителя к стенке радиатора, аа -коэффициент теплоотдачи от стенок радиатора воздуху, 51 - толщина стенки, - коэффициент теплопроводности металла трубок радиатора.
Если структура потока теплоносителя соответствует модели идеального перемешивания, то для математического описания этого потока можно использовать уравнение (с учетом источника тепла, возникающего в потоке за счет теплопередачи):
Шс'Сс
= Шс-Сс • (Тс4 -Тсе) + Шс-дт,
(42)
F ■ Kh ex -кТ
qT
mc
или Qc — Fh_ex ■ Kh_ex ■ (Tc - Ta )
(43)
здесь: Qc - количество теплоты отводимое от двигателя охлаждающей жидкостью, Qa - количество теплоты переносимое воздухом, m a - массовый расход воздуха через радиатор, m c - массовый расход теплоносителя в радиаторе, mc - масса охлаждающей жидкости, содержащейся в радиаторе, ca - теплоемкость воздуха, cc - теплоемкость теплоносителя (этиленгликоля), Сс - теплопроводность, Ta,e - температура воздушного потока на выходе радиатора, TGii - температура воздушного потока на входе радиатора, Tc,e - температура охлаждающей жидкости на выходе радиатора, Tci - температура охлаждающей жидкости на выходе радиатора, 8 - эффективность радиатора, Ga - производительность вентилятора, определяется ранее, исходя из конструктивных особенностей, dTa/dt - изменение температуры воздуха в радиаторе, mc *qT -интенсивность теплообмена, AT = (Т1-Т2) - разность температур первичного и вторичного теплоносителей.
Для теплообменного аппарата, в котором структура потоков соответствует моделям идеального перемешивания, математическая модель представляет собой систему двух уравнений, записанных для первичного и вторичного теплоносителей:
mc ■ cc
ma ■ ca
dTc,e dt
'dTa,. dt
— m c
— m a
Cc ■ (Tc,i Tc,e ) Fh_ex ' Kh_ex ' (Tc,e Ta,e ),
Ca ■ (Ta,i Ta,e ) + Fh_ex ' Kh_ex ' (Tc,e Ta,e
(44)
Величина ¥}1_ех*К}1_ех*(Т1-Т2) имеет знак «минус» в уравнении описания потока теплоносителя, который отдает тепло, и знак «плюс», если тепло воспринимается теплоносителем. Математическая модель исследуемого теплообменника должна включать:
• уравнения теплового баланса;
• уравнения теплоотдачи;
• уравнение переходного режима (динамическая характеристика), показывающее, что при теплообмене в неустановившемся потоке, когда dQ1 и dQ2 , часть тепла аккумулируется стенкой в любой момент времени переходного режима [8].
Соответствующие уравнения записаны в виде
2 dQc
T —
T —
^ a ,e
a i ■ Fi 2
dt
+ 2 ■ Tst Tci,
dQa
a 2 ■ F2 dt
+ 2 Tst-T,
-1 st
I dTst — Tst — Tst 0 —
dQc dQa
1 st 0
dQc dt dQa dt
Gst ■ cst 0
dt dt
dt,
m c
— m a
C ■(T ■ -T ) ■ C ■(T -T )
(45)
где: Тщ, Тзю - текущая и начальная температура стенки радиатора, G^ - вес стенки, сл - теплоемкость стенки радиатора, а^ а2 - коэффициенты теплоотдачи теплоносителей. Схема к уравнению теплового баланса для двигателя представлена на рисунке 16.
Уравнение, описывающее тепловой баланс в двигателе
А ЛТ Л ■
ее --Шс ■ Сс ■(■ -) + ()е
Ше ■ С с
у Ж J
(46)
здесь: ше - масса теплоносителя в двигателе, Те^ и Те,е - температура теплоносителя на входе
в двигатель и выходе из него, (Ое - мощность тепловыделения двигателя.
Уравнения, описывающие тепловой баланс в трубопроводах представлены двумя системами уравнений (47) и (48), соответственно, для "Трубопровод 1" и "Трубопровод2":
Г ^ \
ШР1 ■ Сс
(СТ,
у Л J
= Ш с ■ Сс ■(Те,е - Тс; ) + 0,
рЪ
С0р1 Лг
Шс ■ Сс ■ (Те,е Тс■ ),
Т ■ = ■
1 П 1
С(р1
а р1 ■ ¥р1 Л
+ 2 ■ Т — Т
~ ^ 1 рЫ 1 е,е
(47)
ЫТ, ^
Шр2 ■ сс
у Сг J
Шс ■ Сс ■ (Тс,е - Те,г ) + 0р2,
С0р 2 Л
Ш
Т
^ о
^с ■ Сс ■ (Тс,е Те,/ ),
2 С0р2
а р2 ■ Рр2 Сг
+ 2 ■ Т — Т
' ^ 1 р21 с,е
(48)
в системах уравнений (47) - (48): шр1, шр2 - массы охлаждающей жидкости, содержащиеся в 1-ом трубопроводе, полагаем, что шр1 ~ шр2, (Ор1, (Ор2 - мощность тепловыделения 1-ым трубопроводом.
Размерности всех величин, приведенных в статье и указанных в таблице 1, даны в системе СИ.
Таблица 1.
Значения используемых величин.
Кчг % К2ь Рь V, Е, г V г "■тт Чтт 1. вг 1*Ы
18 6,679* 10"5 0,97 4,667*10'10 24,0*10® 33.8*10"6 !,65*10'6 4,13К*10"6 0,91 4,138*10'6 0,892 0,39 1,95 0,325
г ш Ц Ьы Оы Чг Аг Кьс М и» Й)
0,146 8,024 9 0,02 0,648 0,35 0,132 54,427 75 0 0,3 912 24
В К V N а (1 $ ¿V Р НУ V. V* Р| а
22 24 0,0237 200 0,01 0,015 0,001 8*10'? 890 0,82 4,712*10'4 26*10"* 0,27*106 32
са I ■ «й «1 а2 О'е V,
3840 1005 60 915 75 55 157000 7*10'3 1.45 1,45
Из полученных результатов видно, что для охлаждения двигателя мощностью 272 кВт (370 л.с.) необходим гидропривод с потребляемой мощностью 23 кВт. Данный привод обеспечивает поддержание допустимой температуры автомобильного двигателя в установленных границах +85-95°С. Промоделированные условиях работы - максимальной нагрузка на а/м двигатель и максимальная температура окружающей среды - +60°С (Установившаяся температура а/м двигателя составила +90,41°С).
Рисунок 17 - Переходные процессы Рисунок 1S - Изменение температуры: А -
системы: А - перемещение температура охлаждающей жидкости на выходе
золотника, Б - обороты вентилятора. из двигателя.
Полученная математическая модель может быть использована для анализа систем охлаждения автомобильных двигателей любой мощности с применением гидростатического привода.
Mатематическая модель пригодна для анализа системы охлаждения а/м двигателя при разных режимах работы а/м двигателя и разной окружающей температуры. Схема может быть использована для оптимизации систем охлаждения ДВС с гидроприводом.
В полученной математической модели коэффициенты а1, а2 - коэффициенты теплоотдачи теплоносителей, были взяты из опыта проектирования теплообменных аппаратов. Для более точного представления их необходимо вычислять, используя теорию подобия, определяя коэффициент Nu.
Литература
1. Труханов К. A. Гидропривод вентилятора для системы охлаждения автомобильного двигателя. // Журнал «Известия MГТУ «MAM№> Выпуск №2 (12) 2011.
2. Попов Д. Н. Динамика и регулирование гидро- и пневмосистем. -M.: Mашиностроение, 1977 Учеб. для машиностроительных вузов 424 с. с ил.
3. Семенов С.Е. Электромеханические преобразователи электрогидравлических следящих приводов. -M.: MГТУ им. Н.Э.Баумана, 1997. - 27 с.: ил.
4. Janusz Zieba. Simulation of a Solenoid Actuator for a Device for Investigating Dynamic Air Permeability Through Flat Textile Products. // Technical University of Lodz. April/June 2003, Vol.11, No. 2(41).
5. Mошенцев Ю.Л. Теплообменные аппараты ДВС: Учебное пособие. - Николаев: Николаевская областная типография, 2006. - 431с, ил.
6. Кригер A.M. Жидкостное охлаждение автомобильных двигателей Текст. M.E. Дискин, A.JI. Новенников, В.И. Пикус M. Mашиностроение, 19В5. 176 с.
7. M. H. Salah, T. H. Mitchell, J. R. Wagner, Ph.D., PE, and D. M. Dawson, Ph.D. Nonlinear Control Strategy for Advanced Vehicle Thermal Management Systems. // Automotive Research Laboratory Departments of Mechanical and Electrical Engineering Clemson University, Clem-son, SC. 2007.
8. Турупалов В.В., Чернышев Н.Н., Прядко А.А. Математическое описание процесса теплообмена в противоточных теплообменных аппаратах. // Научные труды ДонНТУ. Выпуск
Рисунок 19 - Структурная схема гидропривода системы охлаждения автомобильного
двигателя с гидроприводом вентилятора: А - схема моделирования гидропривода; Б - схема моделирования дискретного гидрораспределителя; В - схема моделирования
теплового радиатора системы
