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

книги / Математическое моделирование процессов механической обработки

..pdf
Скачиваний:
5
Добавлен:
12.11.2023
Размер:
1.97 Mб
Скачать

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

Выражение для изображения тока I1(р) определим следующим образом:

I

p

E p

 

,

(44)

Zвх p

 

1

 

 

 

 

 

 

 

где Zвх p – входное операторное сопротивление цепи.

 

 

 

 

 

 

 

R

 

1

 

p2

L R C p L R

 

 

Zвх p p L

 

R

 

 

 

 

 

 

 

 

,

 

p C

 

 

 

 

 

p R C 1

 

 

 

 

 

 

 

 

 

p C

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

E

 

p R C 1

 

 

 

p

E

 

p

2

L R C p L R

 

 

 

 

 

 

 

 

 

 

I

 

 

 

 

 

 

 

 

.

 

 

 

 

 

 

 

 

 

 

 

1

 

p

 

 

p R C 1

 

 

 

p3 L

R C p2 L p R

 

 

 

 

 

 

 

 

 

Получено уравнение в операторном виде, описывающее изменение тока в цепи (тягового усилия на ходовом винте привода подач исполнительного органа) при нулевых начальных условиях.

Полученные уравнения позволяют подробно проанализировать работу механизма (как, впрочем, и любого другого механизма) станка в следующих режимах:

1.Режим пуска двигателя, когда суппорт станка из состояния покоя приводится в движение, т.е. пуск привода на холостом ходу, начальные условия в этом случае равны нулю, так как массы не двигаются, валы не растянуты.

2.Режим переключения скорости движения суппорта станка

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

3.Режим торможения, когда тормозным устройством в приводе скорость движения суппорта изменяется от какого-то определенного значения до нуля. Начальные условия в этом случае не равны нулю.

51

4.Режим, когда на суппорт начинает действовать сила резания (скажем, PZ) – при врезании резца в заготовку, закрепленную

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

5.Режим стационарной работы привода, когда на суппорт действует сила резания PZ, изменяющаяся, например, по косинусоидальному гармоническому закону. Этому режиму будет соответствовать электрическая схема, в которой к выходу цепи приложен источник тока с гармоническим воздействием. Начальные условия в этом случае не равны нулю, внутренние источники будут оказывать существенное воздействие на протекание стационарного процесса в приводе станка.

Проверка адекватности метода электроаналогии

Рассмотрим вынужденные колебания линейной системы с одной степенью свободы с учетом трения. Принципиальная схема такой системы приведена на рис. 3. В этой схеме тело массой m удерживается упругой связью с жесткостью С с учетом трения (демпфером).

1. Классическое решение методом теории колебаний. Рассмотрим воздействие на систему с вязким трением внезапной нагрузки:

0 (t 0),

P t P0 t 0 .

В этом случае уравнение движения груза массой m имеет вид

m x C x R P0 ,

где m – масса груза (тела); с – жесткость пружины; R – сила трения. Зависимость силы трения R от смещения или скорости дви-

жения определяется физической природой трения. Наиболее про-

52

стым случаем является так называемое вязкое трение, когда сила трения пропорциональна скорости движения:

R a x.

(45)

Тогда уравнение будет иметь вид

m x C x a x P0 ,

при этом можно получить уравнение вида

 

 

 

 

 

 

x 2 n x p2 x P ,

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

где n

a

;

p2

C

;

p

C

 

– угловая частота собственных

2m

m

m

 

 

 

 

 

 

 

 

колебаний соответствующей консервативной системы.

Решение полученного уравнения будет определяться следующей формулой (промежуточных выкладок не приводим; графическая интерпретация решения приведена на рис. 24):

 

 

 

 

 

 

 

 

P

 

 

 

 

n

 

 

 

 

 

 

 

 

 

x

0

 

1 e n t cos p1 t

 

 

 

sin p1

t ,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

C

 

 

 

 

p1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

где p

 

p2 n2

угловая

частота

 

затухающих колебаний;

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

– период этих колебаний, т.е. время между двумя последо-

 

1

 

p1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

вательными максимальными отклонениями AK и AK+1 составляет

 

 

 

 

 

 

 

 

 

 

 

 

AK

en 1 ;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

AK 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

A

 

 

 

 

 

 

 

 

 

 

 

 

 

ln

 

 

K

 

n 1

– логарифмический декремент колебания.

 

 

 

 

 

 

 

 

AK 1

 

 

 

 

 

 

 

 

 

 

 

53

В данном случае

n

2 n

 

 

 

 

2 n

 

,

 

 

 

 

 

 

 

 

1

p1

 

 

p2 n2

 

 

 

p1

 

 

 

p

 

 

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1 2

 

 

 

 

 

 

 

 

 

 

 

2 2

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

Рис. 24. Затухание колебаний линейной системы

Решим пример при следующих численных значениях величин параметров системы с одной степенью свободы:

m = 20 кг; P0 =100 Н;

α = 100 (Н с)/м; С = 500 Н/м.

При этом

 

 

 

 

 

 

 

p

500

 

5 рад/с;

20

 

 

 

 

54

n

100

2,51/с;

2 20

 

 

p1 25 6,25 4,331/с;

n 0,577. p1

Получим уравнение

x 100500 1 e 2,5 t cos 4,33t 0,577sin 4,33t .

Окончательное уравнение

x 0,2 1 e 2,5 t cos 4,33t 0,577sin 4,33t .

Для получения выражения для скорости движения груза продифференцируем данное уравнение по времени:

x V dxdt d 0,2 0,2e 2,5 t cos 4,33t 0,577sin 4,33t .

Не приводя промежуточных выкладок, запишем окончательный результат:

x V 1,154e 2,5 tsin 4,33t .

Решение операторным способом по методу электроаналогии

Механическая цепь линейной системы с одной степенью свободы с учетом трения приведена на рис. 4. По механической цепи составлена эквивалентная электрическая схема, которая приведена на рис. 5. Операторная схема замещения с нулевыми начальными условиями приведена на рис. 25.

55

Рис. 25. Операторная схема замещения системы с одной степенью свободы

Определим входное операторное сопротивление электриче-

ской цепи Zвх(р):

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Zвх p

p L R 1

 

p L R

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

p L R p C

p L

R

 

 

 

 

 

 

 

 

p C

 

p L R

 

 

p2 L R C p L R

 

 

 

 

p L R

 

 

 

 

.

p C p L R

 

p C p L R

p2 L R C p L R

Найдем изображение напряженияU(p):

U p I p Zвх

p

I

 

 

 

p L R

 

p p2 L R C p L R

 

 

 

 

 

 

I L R

 

.

 

 

 

 

 

 

 

 

 

p2

L R C

 

 

 

 

p L R

 

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

Переход осуществим по формуле разложения

N p

m

N p

 

 

 

 

k

 

e pk t .

M p

M pk

k 1

 

56

Исходные данные для расчета принимаем те же, что и в первом случае:

С = 20 кг;

I = 100 Н;

R = 0,01 м/(Н с);

L = 2 10–9м/Н.

В соответствии с исходными данными выражение для изображения U(p) будет иметь вид

U p

 

2 10 3

 

 

N p

.

 

10 3 p2 2 10 3 p 0,01

M p

0, 4

 

 

Найдем первую производную характеристического уравнения М(р):

M p 0, 2 10 3 p 2 10 3.

Решаем характеристическое уравнение М(р) = 0:

M p 0, 4 10 3 p2 2 10 3 p 0,01 0.

Найдем

корни характеристического уравнения M(p) = 0

(квадратное уравнение):

 

 

 

 

 

 

 

 

 

 

 

 

 

p

 

2 10

3

 

 

2 10 3 2 0, 4 10 3 4 0,01

 

 

 

 

 

 

 

 

1 2

 

0,8 10 3

 

 

0,8 10 3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2,5 j 4,33.

 

 

 

Подставим значения корней в выражение для первой произ-

водной M p1 2 :

 

 

 

 

 

 

 

M p1 2 0,8 10 3 2,5 j 4,33 2 10 3 j 3, 47 10 3.

Преобразуем значения M p1 2 j 3, 47 10 3

в показатель-

ную форму (рис. 26).

 

 

 

 

 

 

 

57

Рис. 26. Комплексная плоскость

Любое комплексное число +jb в нашем конкретном случае можно представить в виде

j b p e j ; arctg b ,

так как j b 3, 47 10 3 , следовательно, = 90 , тогда

3,47 10 3 3,47 10 3 e j 90 .

Подставим полученное выражение в формулу разложения:

U t

2

10 3

e 2,5 j 4,33 t 0,577e j 90 e 2,5 j 4,33 t .

3, 47

 

10 3 e j 90

 

 

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

e j cos j sin ;

(46)

cos

e j e j

(47)

 

 

;

2

 

 

 

 

 

 

sin

 

e j e j

(48)

 

 

;

 

2

 

 

 

 

 

e j cos j sin .

(49)

58

Тогда

Ut 2 0,577e 2,5 t cos 90 4,33t 1,154e 2,5 t sin4,33t.

Внатуральном выражении уравнение будет иметь вид

x t V t 1,154e 2,5 tsin4,33t.

Графическая интерпретация полученного уравнения приведена на рис. 27.

Рис. 27. Переходный процесс в системе с одной степенью свободы

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

Кроме того, получим решение в операторном виде при ненулевых начальных условиях, если таковые имеются на индуктивности L (например, пружина предварительно сжата или растянута) и емкости С; начальные условия на емкости равны нулю, так как скорость движения груза в начальный период времени t = 0

59

равна нулю (V = 0). Операторная схема замещения при ненулевых начальных условиях будет иметь вид, приведенный на рис. 28.

Рис. 28. Операторная схема замещения системы с одной степенью свободы при ненулевых начальных условиях

Пользуясь законами Ома и Кирхгофа, найдем изображение напряжения U(p) в цепи.

Индексы сопротивлений Z(p) и напряжений Е(р) по токам цепи:

 

 

 

 

 

 

 

 

 

Z1 p p L,

 

 

 

 

 

 

 

E p

I

 

 

p L i 0 L I L i 0 L,

 

 

 

 

 

1

 

 

 

 

 

p

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Z2 p

 

p L R

 

,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

p L R

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

E2 p

 

E1 p

Z2 p

 

I L i 0 L

 

p L R

 

 

 

Z

p

 

 

 

p L

 

 

 

p L R

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

I L R i 0 L R

,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

p L R

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

1

 

 

 

 

 

 

p L R

 

Z3 p Z2 p

 

 

 

 

 

 

Z2 p

 

 

 

 

 

 

 

 

 

,

 

 

 

 

 

 

 

 

p2 L R C

p L R

 

 

p C

 

 

p C

 

 

60

Соседние файлы в папке книги