Файл: Никитенко В.Д. Подготовка программ для станков с числовым программным управлением.pdf

ВУЗ: Не указан

Категория: Не указан

Дисциплина: Не указана

Добавлен: 20.06.2024

Просмотров: 141

Скачиваний: 0

ВНИМАНИЕ! Если данный файл нарушает Ваши авторские права, то обязательно сообщите нам.

вписанного в него треугольника высотой б связана с ним следующим неравенством:

Sc e r M S - j - .

(13)

Раскроем левую и правую части неравенства (13). Площадь сегмента 5С 0 Г М можно определить, вычитая из площади эллиптического сектора площадь треугольника, образованного двумя соседними радиусами и хордой /:

"^сегм = = *"*сект ' "^треуг* ( ^ )

Воспользовавшись формулой для определения площади

сектора,

получим

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ф + Д ф

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Ф

 

 

 

 

 

 

 

 

 

и вычислив интеграл, будем иметь

 

 

 

 

 

 

 

 

* с е к т

= 4

[™*g ( - ^ ( Ф + М ) -

 

arctg ( * ? _ ) ] . (15)

Площадь

треугольника

 

 

равна

 

 

 

 

 

 

 

 

 

 

 

 

5•Чреу

гг

= 4"

r

'

/+l

S l n Д <

Р =

 

 

 

 

 

 

 

 

т р

 

2

 

r

 

 

 

 

 

 

 

 

 

 

 

 

 

= —

 

 

 

s i n АФ

 

 

 

 

 

 

 

 

 

 

(16)

 

 

 

 

2

V [ l — e2cos2cp] [I — eacosa(q> + Дф)] *

 

Из

этого же треугольника,

 

 

использовав

теорему

коси­

нусов,

определяем

длину

 

хорды

 

 

 

 

 

 

 

 

 

 

/ = ~\[r\ + r\+\ 2nri+i

cos Аф

=

 

 

 

 

 

2 - е 2

[cos2

(ф + Дф) + cos2 <р] —

 

 

 

 

• 2 cos Дф

 

+ e

cos

 

ф] [1 — е

 

cos

 

(ф + Дф)]

 

 

 

 

 

 

 

 

2

 

 

 

 

2

 

 

 

 

 

2

 

2

 

Л « Л 1

(17)

 

 

 

 

[1 — е3 cos2 ф] [ 1 — е2

cos2

(ф + Дф)]

Для облегчения записи введем следующие обозначения:

Vll—

e2 cos2 9]H— е2

 

cos2

(ф + Аф)]

 

= А;

 

 

arctg [ * < Ф + *•>>] - a r c t g [ i ^ L ] = В;

(18)

/ 2 — е2 [cos2

(ф + Аф) + cos2

 

ф] — 2A cos Аф = С.

 

27


Используя формулу (14) и принятые обозначения (18), находим из неравенства (13) значение стрелки прогиба б:

 

 

6

 

£ f l i ! z ^ l .

 

(19)

Значения

Аф,

удовлетворяющие

неравенству

(19),

используются

при

аппроксимации.

 

 

 

Алгоритм

решения выглядит следующим образом:

1)

с шагом изменения Аф вычисляют

значение

 

 

 

г.

= a

АВ — /г sin Am

 

 

 

 

 

F

g

 

 

 

до тех пор, пока при Аф;- = Аф0 +

nh не

будет

иметь

место

выполнение

неравенства F s£ б;

 

 

2) определяют угол аппроксимации и координаты

точки

аппроксимации

 

 

 

 

i

4>i = 21 Аф,.;

xt = r,cosq>,; yt = гi sin ф,;

3) вычисляют координатные приращения для преды­ дущей точки

 

 

bxi-i=\xl

xi_1\;

 

 

=

\У1 У1-11

(»' = 1, 2, 3

n);

4)

повторяется

расчет до

достижения

углом ф зна­

чения

я/2.

 

 

 

На основании расчетов, выполненных этими двумя способами с помощью ЭВМ «Наири-С», составлены таб­ лицы, позволяющие при заданных a, b и б определить коор­ динаты точек аппроксимации от 0 до л/2. Таблицы содер­ жат значения угла аппроксимации ф через шаг Аф, вы­ численный одним из приведенных выше способов, а также значения координат х и у в соответствующей промежуточ­

ной точке и

абсолютные величины их приращений Ах

и А(/ (табл.

6).

Содержание таблиц иллюстрируется графиком, при­ веденным на рис. 11, построенным для фиксированных

значений б = 0,01 мм и а = 10 мм, но для разных

зна­

чений

Ь.

 

На

рис. 12 представлены семейства кривых х — f (Ф)

и у = f (ф). Приведенные графики и таблицы могут

быть

28


 

 

 

 

 

 

 

 

Таблица 6

Аппроксимация

дуг эллипсов а =

10 мм, Ъ =

6 мм,

б = 0,01 мм

точки

Ф, ра д

Дф, р а д

X,

мм

&х, мм

у,

мм

Ау, мм

0

0,00

0,05

10,00000

0,03460

0,00000

0,49869

1

0,05

0,05

9,96540

0,10235

0,49869

0,49092

2

0,10

0,05

9,86305

0,16595

0,98961

0,47597

3

0,15

0,06

9,69709

0,27400

1,46557

0,54289

4

0,21

0,06

9,42309

0,34255

2,00846

0,50465

5

0,27

0,06

9,08054

0,39604

2,51311

0,46154

и т. д.

 

 

 

 

 

 

 

 

полезны для программистов при определении

шага ап­

проксимации

и величин

приращений

по

координатам

в случае программирования дуги эллипса, заданного своими полуосями. Кривая линия при задании для линей­ ного интерполятора аппро­

ксимируется ломаной, для

мм

х=х(ц>)

а*10пп

У'У(<Р)

S=0,01мм

линейно-кругового—лома­

 

 

 

 

 

ной или дугами

окружно­

 

 

 

 

сти,

для

параболическо­

\ \

b

 

 

го — параболами.

 

y

 

Ар

 

 

 

 

 

 

 

рад

 

 

1

 

 

 

 

0.1?

 

 

S*0,01Mм

 

 

 

 

 

 

 

 

 

 

 

0,08

 

 

b =8

 

'

\ \

\

 

6

r

N

 

 

 

 

 

 

 

 

 

 

Aо,*

b=2

 

 

0,00

 

 

b -2

 

 

 

0,i

0J

1,2 Ц.рад

0,8

 

1,2 4>,рад

Рис.

11. Зависимость

шага ап­

Рис. 12. Значения

координат

проксимации эллипса от угла ф

промежуточных точек при ап­

 

 

 

 

проксимации дуг эллипса

Количество участков между смежными узловыми точ­ ками зависит от характера аппроксимируемой кривой, заданной точности на аппроксимации и аппроксимиру­ ющей функции.

Большая часть контуров деталей в машиностроении (ориентировочно до 94%) составляется простейшими гео­ метрическими элементами: прямыми и дугами окружности.

29



Рис. 13. Параболическая интерполяция

Остальные геометрические элементы не принадлежат к простейшим. Они могут быть заданы аналитическим, табличным, графическим способом [14].

При аналитическом способе задания геометрический элемент задается либо уравнением вида у = f (х), либо параметрическими уравнениями у = у (t), х = х (t). При табличном задании точки контура задаются числовыми координатами, то есть известно, что при значениях абс­ циссы х = Л'0, хъ . . ., хп ордината принимает значения У о, Уъ • • -I Уп- При графическом способе задания дается

чертеж геометрического элемента или всего кон­ тура. Геометрический способ задания весьма не точен и применяется редко.

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

ское выражение простое. При задании функции у = / (х) таблицей значений и в случае, когда аналитическое вы­ ражение для / (х) известно, но является слишком слож­ ным, применяется замена функции / (х) ее интерполя­ ционным многочленом [10].

Геометрически задача отыскания функции / (х) по заданным ее частным значениям означает, что требуется

построить кривую, проходящую

через

точки

плоскостей

с координатами (,v„, у0), (хъ уг),

. . .,

(х„, уп),

а так как

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

Такого рода задачи называются задачами параболи­ ческой интерполяции. Точки х0, х ь . . ., хп называются узлами интерполяции (рис. 13). Многочлен F (х), удо­ влетворяющий поставленным условиям, называется ин­ терполяционным многочленом, а формулы для его построе­ ния — интерполяционными формулами. Задачи парабо-

30