Теория и проектирование гидропневмоприводов практикум
.pdfПо аналогии с (6.12) запишем дифференциальное уравнение движения поршня исполнительного цилиндра в общем виде
m |
d2 z |
p |
F |
p |
F |
|
P |
sgn dz k |
в |
dz |
c |
c z . (6.15) |
||
п |
dt |
2 |
|
2 п |
|
шт шт |
|
|
тр |
dt |
|
0 |
1 |
|
|
|
|
|
|
|
|
|
|
|
dt |
|
|
||
Пренебрегая сопротивлением сливной магистрали (давление в полости слива в гидроцилиндре отсутствует) и силами трения в исполнительном гидроцилиндре, получим
m |
d2 z |
p |
F |
c |
c z . |
(6.16) |
|
п |
dt |
2 |
|
2 п |
0 |
1 |
|
|
|
|
|
|
|
|
|
Используя уравнение баланса расходов (6.7) для узла У2, установим связь между параметрами движения столба рабочей жидкости x в трубопроводе и параметрами движения поршня z исполнительного цилиндра
Qвх2 Qвых2 ,
где Qвх2 , Qвых2 – соответственно расходы на входе и выходе узла У2, причем
Qвх2 Qтр; |
|
Qвых2 Qцил, |
|||||||||||
где Qтр, Qцил – соответственно расходы в трубопроводе и в гидро- |
|||||||||||||
цилиндре. |
|
|
|
|
|
|
|
|
|
|
|
|
|
Q f |
|
|
dx |
; |
|
Q F |
dz |
. |
|||||
|
|
|
|
|
|||||||||
тр |
тр dt |
|
|
|
цил |
п dt |
|||||||
Уравнение баланса расходов для узла У2 приобретает вид |
|||||||||||||
|
f |
|
|
dx |
|
F |
dz |
, |
(6.17) |
||||
|
тр dt |
|
|||||||||||
откуда |
|
|
|
п dt |
|
|
|
||||||
|
dx |
|
Fп |
dz . |
|
|
|
||||||
|
|
|
(6.18) |
||||||||||
|
|
|
|
||||||||||
|
|
dt |
|
fтр dt |
|
|
|
||||||
71
Обозначим
a |
Fп |
(6.19) |
|
||
0 |
fтр |
|
|
|
и продифференцируем уравнение (6.18)
d2 x |
d2 z |
(6.20) |
dt2 |
a0 dt2 . |
Уравнение движение поршня главного цилиндра в общем виде
mпг d2 z Paг Pcг,
dt2
где mпг – масса подвижных элементов, приведенная к поршню главного цилиндра;
Раг и Рсг – соответственно сумма активных сил и сил сопротивления, действующих на поршень главного цилиндра.
Активной силой будет усилие на штоке Pпед(t), а силами сопротивления – силы сухого и вязкого трения и сила давления в поршневой полости. Сила давления
Pcг Pдг p1F1.
Если пренебречь трением и инерционностью приведенных к поршню главного цилиндра подвижных масс цилиндра (mпг = 0), из уравнения движения поршня главного цилиндра
p |
Pпед(t) |
. |
(6.21) |
|
|||
1 |
F1 |
|
|
|
|
||
Подставляем (6.21) в уравнение движения жидкости (6.14) и с уче-
том (6.18) и (6.20) получим
a a |
d2 x |
|
Pпед(t) |
p |
|
a2a |
dx 2 sign |
dx |
a a dx . |
(6.22) |
|||
|
|
F |
|
dt |
|
||||||||
0 1 |
dt |
2 |
|
|
2 |
0 |
2 dt |
|
0 3 dt |
|
|||
|
|
|
1 |
|
|
|
|
|
|
|
|
|
|
72
Разделим левую и правую часть уравнения (6.16) на Fп, просуммируем полученное уравнение с уравнением (6.22) и получим нелинейное дифференциальное уравнение второго порядка, описывающее динамику привода с главным цилиндром на входе (см. рис. 6.5),
mпFп
a a |
|
d |
2 |
z |
|
P |
|
(t) |
|
|
пед |
|
|||||
0 1 |
dt2 |
|
|
|
|
|||
|
|
F |
|
|||||
|
|
|
|
|
|
|
1 |
|
a0a3 ddzt
2 |
dz 2 |
dz |
|
||
a0 a2 |
|
|
sign |
|
|
|
dt |
dt |
(6.23) |
||
c0 c1z . Fп
6.2.2. Моделирование привода с распределителем на входе
Составим уравнения динамики для схемы по рис. 6.6.
Процесс моделирования движения жидкости в трубопроводе между узлами У1 и У2 и процесс моделирования движения поршня гидроцилиндра осуществляются так же, как и для расчетной схемы по рис. 6.5 (см. раздел 6.2.1), в результате чего получаются дифференциальные уравнения:
a |
d2 x |
p |
p |
a |
|
dx 2 sign |
dx |
a |
dx |
; |
||||||
1 |
dt |
2 |
1 |
|
|
2 |
|
2 |
|
|
|
|
3 |
|
||
|
|
|
|
|
|
|
|
dt |
|
dt |
|
dt |
|
|||
|
|
|
m |
|
d2 z |
p |
F |
c c z . |
|
|
|
|||||
|
|
|
п |
dt |
2 |
|
|
2 |
п |
0 |
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
(6.24)
(6.25)
Из уравнения (6.17) баланса расходов рабочей жидкости для узла У2 получаем выражения, устанавливающие связь между параметрами движения столба рабочей жидкости x в трубопроводе и параметрами движения поршня z исполнительного цилиндра:
|
dx |
a |
|
dz |
; |
(6.26) |
|
dt |
|
||||
|
|
0 dt |
|
|||
d2 x |
|
|
d2 z |
(6.27) |
||
dt2 |
a0 dt2 . |
|||||
73
Для формирования давления p1 на входе в трубопровод воспользуемся уравнением баланса расходов для узла У1:
Qвх1 = Qвых1,
где Qвых1 – расход жидкости в трубопроводе;
Qвх1 – расход жидкости через распределитель, который можно рассчитать как расход через дросселирующее отверстие.
Q |
f |
|
dx |
; |
Q |
μS |
|
2 pр |
, |
тр dt |
|
|
|||||||
вых1 |
|
|
вх1 |
|
р |
ρ |
|||
где μ – коэффициент расхода распределителя;
Sp – площадь рабочего окна распределителя, причем S р является функцией времени Sp Sp (t) и зависит от величины окна h(t),
которая моделируется как типовое входное воздействие;
pp – перепад давления на напорном канале распределителя.
pр pmax p1.
Для золотникового распределителя при отсутствии перекрытия можно записать
S(t) πDзh t .
где Dз – диаметр золотниковой пары. Таким образом,
Q |
μπD h t |
2(pmax p1) |
, |
|
|||
вх1 |
з |
ρ |
|
|
|
||
а уравнение баланса расходов жидкости для узла У1 принимает вид
μπD h t |
2(pmax p1) |
f |
dx |
, |
(6.28) |
|
ρ |
тр dt |
|||||
з |
|
|
|
74
откуда
p1 pmax a10
h2 (t)
где
dx 2 |
dx |
, |
||
|
|
sign |
|
|
dt |
dt |
|
||
0,5ρf 2 a10 (μπDзтр)2 .
(6.29)
(6.30)
Подставим (6.29) в уравнение (6.24)
|
a |
d2 x |
p |
|
p |
|
|
|
|
a10 |
|
a |
|
dx 2 sign |
dx |
a |
|
dx |
. (6.31) |
|||||||||||
|
|
2 |
max |
|
|
|
|
|
|
2 |
|
|
|
|
||||||||||||||||
|
|
1 |
dt |
|
2 |
|
|
|
|
2 |
(t) |
|
|
|
|
|
|
|
3 |
dt |
|
|||||||||
|
|
|
|
|
|
|
|
h |
|
|
|
|
|
dt |
|
|
dt |
|
|
|
||||||||||
С учетом (6.26) и (6.27) уравнение (6.31) приобретает вид |
|
|||||||||||||||||||||||||||||
a a |
d2z |
p |
|
p |
a2 |
|
|
a10 |
|
|
a |
|
|
dz 2 sign |
dz |
a a |
|
dz |
. (6.32) |
|||||||||||
|
2 |
max |
|
|
|
|
|
2 |
|
|
|
3 |
|
|||||||||||||||||
0 1 |
dt |
|
|
2 |
0 |
|
|
2 |
(t) |
|
|
|
|
|
|
0 |
dt |
|
||||||||||||
|
|
|
|
|
|
|
h |
|
|
|
|
|
dt |
|
dt |
|
|
|
|
|||||||||||
Разделим левую и правую часть уравнения (6.25) на Fп , просуммируем полученное уравнение с уравнением (6.32) и получим нелинейное дифференциальное уравнение второго порядка, описывающее динамику привода с распределителем на входе (см. рис. 6.6)
mпFп
|
d2 z |
|
|
a2 |
a |
|
|
|
|
|
dz |
|
2 |
|
|
dz |
|
|
|||
a a |
|
|
p |
|
|
|
10 |
a |
|
|
|
|
sign |
|
|
||||||
2 |
max |
|
|
2 |
|
|
|
|
|
|
|
||||||||||
0 1 |
|
0 |
|
2 |
(t) |
|
|
|
|
|
|||||||||||
|
dt |
|
|
|
|
h |
|
|
|
dt |
|
|
|
dt |
|
(6.33) |
|||||
|
|
|
|
|
dz |
|
c0 c1z |
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
|
|
|
a0a3 dt |
|
|
|
F |
|
. |
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
п |
|
|
|
|
|
|
|
|
|
|
|
|
6.3. Моделирование гидропривода с учетом сжимаемости рабочей жидкости
6.3.1. Общий подход к моделированию
Степень сложности математической модели гидропривода зависит от количества принимаемых во внимание сосредоточенных объемов жидкости и мест (узлов), где сконцентрированы эти объемы. Напри-
75
мер (рис. 6.7), можно допустить, что объем жидкости V1, заполняющей главный цилиндр, сосредоточен в узле У1, объем жидкости V3, заполняющей исполнительный цилиндр, сосредоточен в узле У3, а объем жидкости V2, заполняющей трубопровод, – в узле У2. Соот-
ветствующие коэффициенты податливости (pi) являются функциями давления в узлах. Обозначения на рис. 6.7 соответствуют обозначениям математической модели без учета податливости.
|
y |
хтр1 |
|
|
хтр2 |
|
|
z |
|
|
|
|
|
|
|||
|
|
|
R1 |
p2 |
|
|
|
Fп |
Pвх(t) |
p1 |
m1 |
m2 |
R2 |
p3 |
P |
||
|
|
|
|
|
|
z |
||
|
У1 |
|
fтр1, lтр1 |
У2 |
fтр2, lтр2 У3 |
|
||
|
F1 |
(p1) |
|
(p2) |
|
(p3) |
||
|
|
|
|
|||||
Рис. 6.7. Динамическая схема гидропривода с тремя податливостями
Для каждого i-го участка трубопровода, на которые разбивается привод, уравнение баланса давлений в общем виде запишется следующим образом (см. раздел 6.2.1):
p |
p |
a |
d2 xтрi |
a |
|
dxтрi |
2 sign |
dxтрi |
a |
|
dxтрi |
, |
(6.35) |
|||
|
2 |
|
|
|
|
3i |
|
|||||||||
i |
i 1 |
1i |
dt |
|
2i |
dt |
|
dt |
|
dt |
|
|||||
где |
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
a1i lтрi ; |
|
|
|
|
|
|
|
(6.36) |
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
a |
2i |
0,5 |
0,443k i lтрi |
; |
(6.37) |
|||
|
|
|||||||
|
|
i |
fтрi |
|
||||
|
|
|
|
|
|
|||
|
|
a |
|
27,5 lтрi |
. |
(6.38) |
||
|
|
|
||||||
|
|
1i |
|
|
fтрi |
|
||
|
|
|
|
|
|
|||
Для данной схемы составляются два уравнения такого типа, так как гидропривод разбит на два участка.
76
Необходимо также составить два уравнения движения поршней гидроцилиндров, принимая во внимание допущения при моделировании.
Уравнение баланса объемных расходов для каждого из узлов
Qвхi Qвыхi Qдi 0,
где Qвхi , Qвыхi – входящий и выходящий из узла расходы рабочей
жидкости;
Qдi – расход, затраченный на деформацию сосредоточенного в узле объема жидкости.
Эти расходы можно определить по выражениям:
Q |
F dy |
; |
Q |
|
F dz |
; |
|
Q |
|
Q |
f |
|
dxтр1 |
; |
||||
|
|
|
|
dt |
||||||||||||||
вх1 |
1 dt |
|
вых3 |
|
п dt |
|
|
вых1 |
|
вх2 |
|
|
|
тр1 |
|
|||
|
|
|
|
|
|
dxтр2 |
|
|
|
|
ψ p |
|
dp |
|
|
|||
|
Q |
Q |
|
f |
|
|
|
; |
Q |
V |
|
i . |
|
|||||
|
|
тр2 dt |
|
|
|
|||||||||||||
|
вых2 |
|
вх3 |
|
|
|
дi |
|
i |
|
i |
|
dt |
|
|
|||
Можно сделать допущения, что в узле У1 сосредоточен переменный объем жидкости, заполняющий главный цилиндр, в узле У2 – объем жидкости в трубопроводе, в узле У3 – переменный объем жидкости в исполнительном гидроцилиндре. Таким образом,
V1 F1( y0 y); |
V2 fтр1lтр1 fтр2lтр2; |
V3 Fп(z0 z), |
где y0 , z0 – соответственно расстояние между поршнем и крышкой
главного иисполнительного гидроцилиндровв начальномположении. Путем преобразования исходных уравнений, моделирующих внутренние процессы в гидроприводе, можно получить замкнутую систему обыкновенных дифференциальных уравнений, позволяющую
исследовать качество переходного процесса в приводе.
6.3.2. Модель привода с одним сосредоточенным объемом жидкости и главным цилиндром на входе
Рассмотрим динамическую схему гидроцепи с одним сосредоточенным объемом жидкости. Предположим, что весь объем жидкости сосредоточен возле исполнительного цилиндра (рис. 6.8).
77
y |
xтр |
|
|
z |
|
|
|
R |
|
Fп |
|
Pвх(t) |
m |
p2 |
Pz |
||
|
|||||
p1 |
|
|
|||
|
У1 |
fтр, lтр |
У2 |
|
|
F1 |
|
|
|
||
|
|
|
(p2) |
||
|
|
|
|
Рис. 6.8. Динамическая схема гидропривода с учетом сжимаемости рабочей жидкости
Запишем уравнение баланса расходов для узла У2
Qвх2 Qвых2 Qд2 |
0, |
(6.39) |
||||||||||
где |
|
|
|
|
dxтр |
|
|
|
|
|||
Q |
|
f |
|
; |
|
(6.40) |
||||||
|
|
|
dt |
|||||||||
вх2 |
|
|
тр |
|
|
|
|
|||||
Q |
|
|
F |
|
dz |
; |
|
|
(6.41) |
|||
|
|
|
|
|
|
|
||||||
|
вых2 |
|
п dt |
|
|
|
||||||
Q |
V |
ψ p |
|
dp2 |
. |
(6.42) |
||||||
|
||||||||||||
д2 |
|
2 |
|
|
2 |
|
dt |
|
||||
|
|
|
|
|
|
|
|
|||||
Предполагаем, что в узле У2 сосредоточен весь объем V2 жидкости, заполняющей напорную часть привода (трубопровод и переменные объемы в главном и исполнительном цилиндрах),
V2 fтрlтр F1( y0 y) Fп(z0 z),
Подставим (6.40)…(6.42) в (6.39) и с учетом (6.42) получим
|
dxтр |
|
dz |
fтрlтр F1( y0 y) Fп(z0 z) ( p2 ) |
dp |
|
|
fтр |
|
Fп |
|
2 |
0. |
||
dt |
dt |
||||||
|
|
|
dt |
|
(6.43)
(6.44)
78
Уравнение движения поршня исполнительного цилиндра по ана-
логии с (6.15):
m |
d2 z |
p |
F |
p |
F |
|
P |
sgn dz k |
в |
dz |
c |
c z . |
||
п |
dt |
2 |
|
2 п |
|
шт шт |
|
|
тр |
dt |
|
0 |
1 |
|
|
|
|
|
|
|
|
|
|
|
dt |
|
|
||
Если пренебречь сопротивлением сливной магистрали ( pшт 0), инерционностью подвижных масс поршня исполнительного гидро-
цилиндра |
(m 0) и силами трения в нем |
(P |
sgn dz |
k |
dz |
0), |
|
п |
тр |
dt |
|
в dt |
|
последнее уравнение приобретает вид |
|
|
|
|
|
|
откуда |
p2Fп c0 c1z, |
|
|
|
|
|
p2 c0 c1z . |
|
|
|
|
|
|
|
|
|
|
(6.45) |
||
|
Fп |
|
|
|
|
|
Из уравнениядвижения поршня главного цилиндра (см. раздел 6.2.1), пренебрегая трением и инерционностью приведенных к поршню главного цилиндра подвижных масс,
p |
Pпед(t) |
. |
(6.46) |
|
|||
1 |
F1 |
|
|
|
|
||
УравнениебалансадавленийдляучасткаУ1-У2 поаналогиис(6.15):
|
a |
|
d2 xтр |
|
|
p |
p |
|
a |
|
|
dxтр |
2 sign |
dxтр |
a |
|
dxтр |
. (6.47) |
||||||||||||||||||||||||||||||
|
|
|
|
2 |
|
|
2 |
2 |
|
|
|
|
|
|
|
|
3 |
|
||||||||||||||||||||||||||||||
|
1 |
|
dt |
|
|
|
|
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
dt |
|
|
|
|
|
dt |
|
|
|||||||||||||||||
где |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
dt |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
0,443k l |
|
|
|
|
|
|
|
|
|
27,5 l |
|
|
|
|
|
|||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
тр |
|
|
|
|
|
|
|
тр |
|
|
|
||||||||||||||||||
|
a |
l |
тр |
; |
|
a |
2 |
|
0,5 |
|
|
|
|
|
|
|
|
; |
a |
|
|
|
|
|
|
|
|
|
. |
|
|
|||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||||||
|
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
fтр |
|
|
|
|
|
|
3 |
|
|
|
|
|
|
fтр |
|
|
|
|
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
Подставляем (6.45) и (6.46) в (6.47). В результате получим |
|
|||||||||||||||||||||||||||||||||||||||||||||||
|
d2 xтр |
|
|
P (t) |
|
|
c |
c z |
|
|
|
|
|
|
dxтр |
2 |
|
|
|
|
dxтр |
|
|
|
|
dxтр |
|
|
||||||||||||||||||||
a |
|
|
|
пед |
|
|
|
|
0 |
|
|
1 |
a |
|
|
|
|
|
|
|
|
sign |
|
|
|
|
|
|
|
a |
|
|
|
|
|
. |
(6.48) |
|||||||||||
|
|
|
F |
|
|
|
F |
|
|
|
dt |
|
|
|
dt |
|
|
|
|
|
dt |
|||||||||||||||||||||||||||
1 dt2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
2 |
|
|
|
|
|
|
|
|
|
|
3 |
|
|
|
||||||||||||||||||
|
|
|
|
|
|
1 |
|
|
|
|
|
п |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
79
Дифференцируем (6.45) |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
dp2 |
|
c1 |
|
|
dz . |
|
|
|
|
|
|
|
|
|
(6.49) |
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
F |
|
|
|
|
|
|
|
|
|
|
|
||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
dt |
|
|
|
|
|
dt |
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
п |
|
|
|
|
|
|
|
|
|
|
|
|
|
||
Уравнение неразрывности потока рабочей жидкости для узла У1 |
||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
F1y fтрz, |
|
|
|
|
|
|
|
|
|
|
|
||||||||||
откуда |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
F1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
y |
|
|
z. |
|
|
|
|
|
|
|
|
|
|
(6.50) |
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
fтр |
|
|
|
|
|
|
|
|
|
|
||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
Из выражения (6.44) с учетом (6.9), (6.45), (6.49) и (6.50) |
|
|
|
|||||||||||||||||||||||||||||||||
|
|
|
|
|
dxтр |
|
|
dz |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
F |
|
|
|
|
|
|
|
|
|||
|
|
f |
тр |
|
|
|
F |
|
|
|
|
f |
l |
тр |
F ( y |
0 |
|
1 |
z) F (z |
z) |
|
|
|
|||||||||||||
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||||||||
|
|
|
|
|
|
|
п |
|
|
|
|
тр |
|
1 |
|
|
|
|
|
|
п 0 |
|
|
|
|
|
|
|||||||||
|
|
|
|
|
dt |
|
|
|
dt |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
fтр |
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
c |
|
|
|
|
|
|
|
|
|
dz |
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1 |
|
|
dt |
0, |
|
|
|
|
|
|
|
|||||||||||
|
|
|
|
|
|
|
|
|
Е |
|
F a |
|
p |
c c z |
|
|
|
|
|
|
|
|
||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
а0 |
|
п |
|
|
0 |
|
|
1 |
|
|
|
|
|
|
|
|
|
|
|
||||||
откуда |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
fтрlтр F1( y0 |
F1z |
fтр) Fп(z0 z) c1 |
|
|
|
|
|
||||||||||||||||||||||
dz |
f |
|
|
F |
|
|
dxтр |
. |
(6.51) |
|||||||||||||||||||||||||||
|
|
тр |
|
п |
|
|
|
|
|
|
|
|
FпЕа0 |
ap |
c0 c1z |
|
|
|
|
|
|
|||||||||||||||
dt |
|
|
|
|
|
|
|
|
|
|
|
|
|
dt |
|
|||||||||||||||||||||
Таким образом, в окончательную систему уравнений, моделирующую рассматриваемую гидравлическую цепь (см. рис. 6.8), войдут уравнения (6.48) и (6.51), а именно:
|
d2 xтр |
|
P (t) |
c |
c z |
|
|
dxтр 2 |
|
dxтр |
|
|
dxтр |
|
||||
a |
|
|
пед |
0 |
1 |
a |
|
|
|
|
sign |
|
|
|
a |
|
|
; |
|
F |
|
dt |
dt |
|
dt |
||||||||||||
1 dt2 |
|
|
F |
|
2 |
|
|
3 |
|
|||||||||
|
|
|
1 |
|
п |
|
|
|
|
|
|
|
|
|
|
|
|
|
dz |
|
|
|
|
|
fтрlтр F1( y0 F1z fтр) Fп(z0 z) c1 |
|
dxтр |
|
f |
|
F |
|
|
. |
||||
dt |
|
тр |
|
п |
|
FпЕа0 ap c0 c1z |
|
dt |
|
|
|
|
|
|
|
||||
80
