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

Система автоматизированного моделирования стрелового крана Монография Омск

.pdf
Скачиваний:
59
Добавлен:
19.04.2019
Размер:
3.16 Mб
Скачать

где W 1;

W

 

γж з

s

ζз γж

Q ;

 

 

 

 

 

 

 

11

12

 

 

g Sз

 

g

Sз2

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

γж з

 

 

 

ζз γж з

 

2

 

 

 

 

 

W

g Sз2

 

Q

 

 

g Sз3

Q ; W 0; W 1; W 0.

13

 

1

 

 

 

1

21

 

22

23

Разветвлениегидролинии

Математическая модель гидравлического разветвления:

P2 P1 2ζ2g γSж22 Q22 ;

P3 P1 2ζ3g γSж32 Q32;

Q2 = Q1 Q3;

Q3 = Q1 Q4;

(3.2.13)

где P1, Q1 – соответственно давление и расход жидкости на входе в тройник; P2, Q2, P3, Q3 – давление и расход жидкости на выходах тройника; S2, S3 – площади проходных сечений выходных гидролиний тройника; 2, 3 – коэффициенты гидравлических сопротивлений.

Математическая модель разветвления гидролиний в виде ГМП:

 

 

 

 

 

 

P2

W11

W12

 

 

 

 

 

 

 

 

 

 

 

 

 

P1 ;

 

 

 

 

 

 

 

Q2

 

W21

W22

 

(3.2.14)

 

 

 

 

 

 

P3

W31

W32

Q1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Q3

 

W41

W42

 

 

где W 1; W

 

 

ζ2 γж

Q

; W 0; W

1; W

1;

 

 

 

11

12

 

g

S22

2

21

 

22

31

 

 

 

ζ3 γж

 

 

 

 

 

 

 

 

W32

 

Q3; W41 0; W42 1; P1 dP1; Q1 dQ1;P2 dP2 ;

g S32

 

 

 

 

 

 

 

 

 

 

 

 

 

Q2 dQ2 ; P3 dP3 ; Q3 dQ3 .

70

Соединениегидролиний

Математическая модель гидравлического соединительного тройника:

 

 

1

 

 

 

 

ζ

1

γ

ж

Q2

 

ζ

2

γ

ж

Q2

 

 

P

 

 

P

P

 

 

 

 

 

 

 

;

2

2 g S2

2 g S2

3

 

 

1

2

 

1

 

2

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

Q3 = Q1 + Q2;

 

 

 

 

 

 

 

(3.2.15)

где P1, P2 – давление жидкости входных гидролиний; Q1, Q2 – расходы жидкости входных гидролиний; P3, Q3 – соответственно давление и расход в выходной гидролинии; S1, S2 – площади проходных сечений входных гидролиний; 1, 2 – коэффициенты гидравлических сопротивлений гидролиний.

Математическая модель соединения гидролиний в виде ГМП:

 

 

 

 

 

 

 

 

 

 

 

 

P1

 

 

 

 

P

 

W

W W W

 

Q

 

 

(3.2.16)

 

 

3

 

11

 

12

13

14

 

1 ;

 

 

 

Q3

 

 

W21

W22 W23

W24

P2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Q2

 

 

 

где W

0,5; W

 

 

ζ1 γс

 

Q ; W

0,5; W

 

 

 

ζ2 γс

 

Q ;

2 g S2

 

2 g S

 

11

12

 

1

13

 

14

 

2

2

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

2

 

W21 0,5; W22 0; W23 0,5; W24 0; P1 dP1; Q1 dQ1; P2 dP2 ; Q2 dQ2 ; P3 dP3 ; Q3 dQ3 .

Дросселирующий элемент

Математическая модель дросселирующего элемента:

P P

 

γж др

Q

 

ζдр γж

Q2

;

(3.2.17)

g Sдр

2 g Sдр2

2 1

 

1

 

1

 

 

 

 

Q2 = Q1;

 

 

 

(3.2.18)

где P1, P2 – давление соответственно на входе и выходе дросселя; Q1, Q2 – расходы жидкости соответственно на входе и выходе дросселя;

71

Sдр – площадь проходного сечения дросселя; др – коэффициент гидравлического сопротивления дросселя.

Математическая модель дросселирующего элементав виде ГМП:

 

P2

 

W11

W12

W13

 

 

 

 

 

 

 

 

 

 

 

 

 

 

P1

 

,

 

(3.2.19)

 

 

 

 

 

 

 

 

Q1

 

 

 

Q2

 

 

W21

W22 W23

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

SЉр

 

 

 

 

 

 

 

γж

др

 

ζдр γж

 

 

 

 

γж др

 

ζдр γж др

2

 

где W 1; W

 

 

 

 

s

 

Q; W

 

 

 

 

Q

 

Q

;

 

g S

 

g S2

g S2

g S3

11

12

 

 

 

1

 

13

 

1

1

 

 

 

 

 

др

 

др

 

 

 

 

 

 

др

 

др

 

 

W21 0; W22 1; W23 0.

Рассмотрим пример расчетной схемы гидропривода самоходного грузоподъемного крана с телескопической стрелой. Расчетная схема в целом состоит из функциональных блоков: дизельного двигателя, насосной установки, разветвления гидролинии, гидролинии, двух гидромоторов (поворота платформы и лебедки груза), двух гидроцилиндров (телескопирования стрелы и подъема стрелы) (рис. 3.4).

ТЕЛЕСКОПИРОВАНИЕ СТРЕЛЫ

ПОВОРОТ

ПЛАТФОРМЫ

 

ЛЕБЕДКА

ДВС

ГРУЗОВАЯ

ПОДЪЕМ, ОПУСКАНИЕ СТРЕЛЫ

Рис. 3.4. Расчетная схема гидропривода грузоподъемного крана

72

73

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

X4

U5

 

 

 

F5

X5

U6

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

О5

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

F2

 

 

 

 

 

 

X4

U8

 

 

 

F8

X8

U9

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

U1

 

 

 

 

 

 

 

X2

U3

 

 

 

 

 

 

 

О8

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

X

1 U

 

 

 

 

 

X

 

U4

 

 

 

 

 

 

 

 

 

 

О1

2

 

О2

 

О3

3

О4

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

F11

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

X4

U11

 

 

 

 

X11 U12

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

О11

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

F14

X4 U14

 

X14 U15

О14

 

 

 

 

 

 

 

 

 

 

 

 

 

F7

 

 

 

 

 

 

 

 

О6

 

X

6

U

7

 

 

 

X7

 

 

 

 

 

О7

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

F10

 

 

 

X

 

U

 

 

 

 

 

 

 

X10

 

О9

9

10

 

 

 

 

 

 

О10

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

X

 

U

 

 

 

 

F13

X13

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

О12

12

13

 

О13

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

F16

 

X

 

U

 

 

X16

О15

15

16

 

О16

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Рис. 3.5. Расчетная схема гидропривода грузоподъемного крана в виде ГМП

Представим расчетную схему, описанную выше, в виде многомерного гидравлического многополюсника (рис. 3.5).

Составим математическую модель гидропривода грузоподъемного крана, представленного на расчетной схеме на рис.3.5.

Уравнения ГМП для каждого элемента расчетной схемы имеют

вид:

О1 : X1 W1XU U1;

О2 : X2 W2XU U2 W2XF F2; О3 : X3 W3XU U3;

О4 : X4 W4XU U4; X4 W4XU U4; X4 W4XU U4; X4 W4XU U4;

О5 : X5 W5XU U5 W5XF F5; О6 : X6 W6XU U6;

О7 : X7 W7XU U7 W7XF F7; О8 : X8 W8XU U8 W8XF F8; О9 : X9 W9XU U9;

О10 : X

10 W10XU

U10 W10XF F10;

 

О11 : X

11 W11XU

U11 W11XF F11;

 

О12 : X

12 W12XU

U12;

 

 

О13

: X

13

W13XU

U13 W13XF F13;

 

О14

: X

14

W14XU

U14 W14XF F14

;

 

О15

: X

15

W15XU

U15;

 

 

О16

: X

16

W16XU

U16 W16XF F16

.

(3.2.20)

Установим соответствие между векторами входных и выходных величин, связывающих ГМП между собой:

X1 U2; X2 U3; X3 U4; X4 U5; X5 U6; X6 U7; X4 U8;

X8 U9; X9 U10; X4 U11; X11 U12; X12 U13; X4 U14;

X14 U15; X15 U16.

(3.2.21)

74

Получим математическую модель гидропривода грузоподъемного крана:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

U

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

F

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

F

 

 

 

 

 

 

 

 

 

 

 

 

5

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

W1(s) W2(s)

W3(s) W4(s)

0

0

0

0

0

0

F7

 

X7

 

 

 

 

W(s)W(s)

0

0

W(s)W(s)

0

0

0

0

F8

 

X

10

5

6

0

0

7

8

W (s) W (s)

0

0

 

 

 

X

W(s)W (s)

0

0

 

 

13

9

10

0

0

0

0

11

12

W (s)W (s)

F10

 

X

W (s)W (s)

0

0

 

16

13

14

 

 

 

 

 

 

15

16

 

 

 

 

 

 

 

 

 

 

 

 

 

 

F

 

 

 

 

 

 

 

 

 

 

 

 

11

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

F13

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

F14

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

F

 

 

 

 

 

 

 

 

 

 

 

 

16

 

 

 

 

 

 

 

 

 

 

 

(3.2.22)

где W1(s) W7XU W6XU W5XU W4XU W3XU W2XU W1XU ; W2(s) W7XU W6XU W5XU W4XU W3XU W2XF;

W3(s) W7XU W6XU W5XF; W4(s) W7XF;

W5(s) W10XU W9XU W8XU W4XU W3XU W2XU W1XU; W6(s) W10XU W9XU W8XU W4XU W3XU W2XF;

W7(s) W10XU W9XU W8XF; W8(s) W10XF;

W9(s) W13XU W12XU W11XU W4XU W3XU W2XU W1XU; W10(s) W13XU W12XU W11XU W4XU W3XU W2XF;

W11(s) W13XU W12XU W11XF; W12(s) W13XF;

W13(s) W16XU W15XU W14XU W4XU W3XU W2XU W1XU ; W14(s) W16XU W15XU W14XU W4XU W3XU W2XF;

75

W15(s) W16XU W15XU W14XF;

W16(s) W16XF .

Данные уравнения гидропривода позволяют проводить исследование всех векторов гидросистемы.

3.3. Моделирование управляющих и возмущающих воздействий на динамическую систему стрелового грузоподъемного крана

Стреловой грузоподъемный кран испытывает воздействия различной природы: ветровые нагрузки, действие внешних сил, управляющие воздействия.

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

Кран должен работать в любую погоду, в том числе и ветреную. Предельное давление ветра, при котором еще возможна и безопасна работа крана, определяет ветровую нагрузку рабочего состояния, которая учитывается при проверке грузовой устойчивости крана против опрокидывания. Кран в нерабочем состоянии рассчитывается на давление ветра, имеющего место при шторме, буре или урагане, которое учитывается при проверке собственной устойчивости крана против опрокидывания [9].

Ветровые нагрузки определяются в соответствии с ГОСТ 1451-77. Ветровая нагрузка на кран определяется как сумма статической и динамической составляющих. Статическая составляющая, соответствующая установившейся скорости ветра, должна быть учтена во всех случаях. Динамическая составляющая, вызываемая пульсацией скорости ветра, должна быть учтена при проверке устойчивости кранов

против опрокидывания [9].

Скорость ветра характеризуется пульсацией, что предопределяет динамичность ветрового воздействия.

Динамическое давление ветра q связано с плотностью воздуха=1,225 кг/м3 и его скоростью v, направленной параллельно поверхности земли, следующей формулой:

 

v2

1,225 v2

2

 

 

q

 

 

 

0,613 v

 

,

(3.3.1)

2

2

 

 

 

 

 

 

 

76

где скорость ветра выбирается в зависимости от района бывшего

СССР (территория СССР была разбита на 7 районов). Распределенная ветровая нагрузка pВ (Па) на единицу расчетной

площади элемента конструкции или груза определяется по формуле:

PВ = qkc(1+ ),

(3.3.2)

где q – динамическое давление ветра; k – коэффициент, учитывающий изменение динамического давления по высоте; с – коэффициент аэродинамической силы; 1+ = 1+mП 1,2 ... 1,3 – динамический коэффициент, учитывающий пульсацию ветрового воздействия mП и коэффициент динамичности , являющийся функцией периода собственных колебаний T(с).

В нормальных условиях эксплуатации для рабочего состояния крана принимают

kc(1+ ) =1,6 ... 2,5. (3.3.3)

Ветровая нагрузка на конструкцию крана, ее отдельные элементы и груз определяются по формуле:

PВ pВ Fi ,

(3.3.4)

где Fi – площадь, воспринимающая давление ветра, под которой понимают теневую площадь. Площадь груза, если она не известна, определяется в зависимости от грузоподъемности крана Q (т):

Fi 3,2 Q .

(3.3.5)

Для моделирования влияния внешних сил при исследовании динамических систем часто используются простейшие детерминированные функции:

- гармонический сигнал F(t) A j( t ψ) или F(t ) A sin( t );

-импульс f(t ) dl0(t ); dt

-линейная функция f(t)=kt+a;

 

0,

при t 0

 

- ступенчатое воздействие

 

;

(3.3.6)

f(t)

 

 

при t 0

 

 

1,

 

77

0, при t 0

- трапециевидная функция f(t) kt, при 0 t t1 ,

 

 

 

при t t1

h,

где k – коэффициент, характеризующий крутизну склона трапеции и зависящий от скорости.

Внешние воздействия могут задаваться в виде массивов значений t и f(t) с интерполяцией значений в промежуточных точках; либо значениями в проекциях на оси систем координат.

Управление моделью (решением системы уравнений) гидропривода осуществляется генерированием управляющих воздействий: текущих значений подачи топлива ДВС, параметров регулирования гидронасосов и гидромоторов, положения (площадей проходных сечений) золотников и дросселей, давлений гидроцилиндров.

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

4. РАЗРАБОТКА СИСТЕМЫ АВТОМАТИЗИРОВАННОГО МОДЕЛИРОВАНИЯ СТРЕЛОВОГО ГРУЗОПОДЪЕМНОГО КРАНА

4.1. Структура системы автоматизированного моделирования стрелового грузоподъемного крана

Система автоматизированного моделирования СГК предназначена для моделирования стреловых грузоподъемных кранов с целью исследования рабочих процессов, статики и динамики СГК, анализа аварийных ситуаций, анализа устойчивости СГК. Областью применения системы могут быть научные исследования и разработки стрелового гидравлического крана.

Она позволяет без знания языков программирования и численных методов, в понятных и простых для восприятия терминах, с использованием библиотек типовых элементов производить построение модели динамической системы СГК с учетом влияния внешней среды, технологических условий работы.

78

Структурная схема системы автоматизации моделирования СГК приведена на рис. 4.1.

Интерфейс САМ позволяет организовать в наглядной форме ввод основных параметров и характеристик динамической системы СГК, построение структурной схемы гидропривода, задание параметров моделирования, возмущающих и управляющих воздействий, параметров вывода результата.

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

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

Интерфейс САМ СГК

Коррекция

параметров

 

 

 

 

 

 

Программа

Программа

 

Программа

 

Программа про-

 

формирования

 

формирования

 

ведения дина-

 

оформления ре-

модели механи-

 

модели подсис-

 

мического ис-

 

зультатов мо-

ческой подсис-

 

темы гидропри-

 

следования

 

делирования

темы

 

вода

 

(решение систем

 

 

 

 

 

 

уравнений)

 

 

 

 

 

 

 

 

 

Библиотеки типовых элементов механической подсистемы

Библиотеки типовых элементов гидропривода и схем гидропривода

Рис. 4.1. Структурная схема САМ СГК

Библиотеки типовых элементов созданы для хранения и выборки элементов механической подсистемы и подсистемы гидропривода.

79