Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

Теория и проектирование гидропневмоприводов практикум

.pdf
Скачиваний:
1
Добавлен:
29.11.2025
Размер:
836.37 Кб
Скачать

По аналогии с (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 PP,

dt2

где mпг – масса подвижных элементов, приведенная к поршню главного цилиндра;

Раг и Рсг – соответственно сумма активных сил и сил сопротивления, действующих на поршень главного цилиндра.

Активной силой будет усилие на штоке Pпед(t), а силами сопротивления – силы сухого и вязкого трения и сила давления в поршневой полости. Сила давления

PPдг 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

Соседние файлы в предмете [НЕСОРТИРОВАННОЕ]