Дифференциально-геометрический метод образования развертывающихся поверхностей Текст научной статьи по специальности «Физика»

ТЕХНОЛОГИЯ МАШИНОСТРОЕНИЯ

УДК 514.182

К.Л. Панчук, А.С. Нитейский

ДИФФЕРЕНЦИАЛЬНО-ГЕОМЕТРИЧЕСКИИ МЕТОД ОБРАЗОВАНИЯ РАЗВЕРТЫВАЮЩИХСЯ ПОВЕРХНОСТЕЙ

Во многих областях практической деятельности человека применяются линейчатые поверхности: в судостроении при выполнении обшивки корпуса судна, в самолетостроении при построении теоретических моделей элементов горизонтального оперения, в архитектурно-строительном проектировании, при проектировании пространственно-шарнирных механизмов, используемых в конструкциях роботов и манипуляторов, при разработке орудий почвообработки и др. Множественным применениям линейчатых поверхностей соответствуют теоретические исследования в области образования и конструирования этих поверхностей. Одно из недостаточно изученных направлений исследований связано с образованием развертывающихся поверхностей.

В дифференциальной геометрии известно построение всевозможных развертывающихся поверхностей по их пространственной ортогональной траектории [1]. Для общего случая, когда кривая на развертывающейся поверхности не является ее ортогональной траекторией, образование этой поверхности не рассматривалось.

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

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

Рис.1. Положение образующей прямой в подвижном трехграннике кривой линии

Как следует из схемы на рис.1, положение ор-

та определяется векторным уравнением

l = (t cos3 + n sin3)cosa + b sin a ( 1)

При перемещении трехгранника вдоль кривой прямая линия описывает линейчатую поверхность. Будем считать, что угловые параметры а и в прямой зависят от лонгального параметра s кривой. Введение двух параметров а= a(s) и в = в (s) выделяет из комплекса прямых, пересекающих данную кривую, определенную линейчатую поверхность. Выделение развертывающейся поверхности из множества линейчатых поверхностей комплекса прямых происходит по известному в теории линейчатых поверхностей [1] условию

dl . ^ (t,l,—) = 0 . ds

(2)

Из уравнения (1) можно получить выражение для векторной производной единичного орта

dl -

— = -t ds

. d3 da

sin 3 cos a(-------------+ k) + cos 3 sin a---------------

yds J ds

+

+ n

cos 3 cos a(— + k) - sina(sin 3 —— + Z) ds ds

+ b cos a(z sin 3 + ——).

ds

(3)

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

Уравнение (2), с учетом (1) и (3), после соответствующих преобразований, принимает вид

т dL . nda . _

(t, l,—) = sin 3------sin a cos a cos 3X

ds ds (4)

(k + —3) + z(1 - cos2 a cos 32) = 0.

ds

Из последнего уравнения, при в=п/2 получаем выражение da= - x(s)ds, представляющее собой необходимое и достаточное условие существования развертывающейся поверхности, образуемой по ее ортогональной траектории [1].

Для случая плоской кривой (х=0, 0^0,

аф±п/2, вф0) получаем формулу

da

■ - ctg3(kds + d3) = 0 , (5)

sin a cos a

которая может быть положена в основу образования развертывающейся поверхности по принадлежащей ей плоской кривой. Этот случай рассмотрен в [2].

В уравнении (4) присутствуют переменные a(s), в (s), k(s) и X(s), что говорит о множественности образуемых развертывающихся поверхностей, определяемых как геометрией исходной кривой линии, так и положением образующей прямой линии относительно сопровождающего трехгранника этой кривой.

Примем для определенности следующие условия: в=сonst и X=Xk, где X - коэффициент с определенным числовым значением. В таком случае, принимая si^=A, cosв=B, получаем уравнение

kds = A -

B sin a cos a - A(1 - B2 cos2 a)

.. (6)

Введем для a(s) замену: t(s) = tg(a(s)), получая преобразованное выражение (6)

, 7 dt

kds = p—-------------, (7)

где

9

t + bt + m

A , B 2

p =--; b =-; m = 1 - B

A A .

(8)

Интегрирование (7) приводит к уравнению

fkds = fdp = ф = p arctg(t + b/2) + C, (9) h 1

h

где h = ( ^ И 4m - b

2

a = arctg

i ,h . b h ■ tg(-■ф)--P 2

(10)

2

(11)

2а 6а2

где а =сот1. Определим первые три производные радиус'-вектора (11) по переменной Хо.

2

-I - -хо г хо • -» -1 г хо ■

г = г + ]-° + к—; г = у- + к-

2а

_2

(12)

На основании уравнений (12) определим орты сопровождающего трехгранника кривой:

- г' - г" -п

t = = Ш1Г ; п = т^-г = Ш-2 г ;

\г' \г”

b = 7^4 = m3 V>r"].

Г ,r ]

(13)

В формулах ортов переменные параметры определяются следующим образом:

mi =

2a

2

9 9 ’ m2

Xq + 2a

2 , 2 a + xq

m3 =

2a3

2 о 2

Xq + 2a

(14)

Определим кривизну и кручение исходной кривой линии:

k =

_IV v ]_

4a3

V" I3 (2a2 + Xq2)2

Z =

(f'f'f"') 4a3

[f ", f "]2 “ (2a2 + x02)2

(15)

Очевидно Х=^к, Х=1, что указывает на принадлежность кривой линии к классу линий откоса, характеризующихся постоянным значением отношения х(з).к(з) [1]. Таким образом, параметры (8) в случае (Х=1) принимают значения

р = —А; Ь = —Б;

h = 1

(^)V4m - b2 = (У2)45

A2 -1.

Ф - текущий параметр -

угол наклона касательной к исходной кривой. Полагая С=0 и учитывая t(5) = tg («(5)), из (9) получаем уравнение

Для постоянного параметра h очевидно ограничение А = 5гп [ >т]0,2 . В итоге для исходной кривой формула (10) принимает вид:

a = arctg

h

htg (- ф) - B A

(16)

которое при текущем изменении углового параметра ф, на основании уравнения (1) позволяет описать развертывающуюся поверхность. Следует отметить, что значения И в уравнении (10) ограничены, поскольку 4т - Ь2>0 и зависят от значений углового параметра в и коэффициента X.

Рассмотрим пример. Пространственная кривая, на основе которой будет образована развертывающаяся поверхность, описывается уравнением

Если отсчитывать лонгальный параметр 5 кривой в сторону возрастания переменной Хо, то для его элемента можно записать уравнение

^ =71 + (у0(хо))2 + (хо))2 • ^о, для рассматриваемого примера принимающее вид

о 2 , 2

с1в = 2а-+Х2 • ох0 . (17)

2а 2

В этом случае формула для элемента угла наклона касательной будет следующей:

ёр = кйя = —2а—- • йхо . (18)

2а + х0

Интегрируя выражение dф, получаем

V Ох0 и \\х

Р = I -----— + Ро = 2иаг^ (- х0) |

•> х, ™ 10

a + -

+ Фо

2a

2

a

х

0

a

a

Поскольку Р0 = 2иаг^ (— • 0) = 0, где

а

ХоТ(0)=0, то окончательно формула для углового параметра ф принимает вид

ф = 2и ■ аг(— хог ) . а

(19)

где Хот - текущее значение абсциссы вершины сопровождающего трёхгранника кривой; и = 0, 7071. Формулы (19), (16), (13) и (1) позволяют получить развертывающуюся поверхность для рассматриваемого примера. В работе уравнения этой поверхности получены с помощью Мар1е-программы и имеют следующий вид.

е(/, Т) = ґ + Т

4.3(/2 + 8.0) 1

+ 3.4-

ґ

-1.0

/(V 1.0 + (м)2 лj-(м)ґ 2

д/1.0 + (м)2 (ґ2 + 8.0)

Гг

у(/, Т) = -0.25/2 + Т

- 2.2

- + 0.84

8.0 -1.0/2

- 4.0-

/2 + 8.0 /2 + 8.0 /(V 1.0+(м)2 у

(м)

д/1.0 + (м)2 (/2 + 8.0)

у(/, Т) = -0.25/2 + Т

(20)

- 2.2—^-------+ 0.8480 - Ш I /

- 4.0-

/2 + 8.0 >/1.0 + (м)2

(м)/

^1.0 + (м)2 (t2 + 8.0)

где М = —0.80tg (l.3arc tg (0.35t))—0.27, а = -2;

в=1 рад.; Л =1; параметр Хо и текущий параметр положения точки на прямолинейной образующей обозначены соответственно ^ и Т.

Для определения стрикции полученной поверхности развертывающейся поверхности необходимо воспользоваться ее уравнением [1]

р = г

■ I.

(21)

в котором векторная производная есть производная сложной функции

-, й1 й1 йа йф ^22)

йа йф йх0

После ввода обозначений

Ни В

— = а; — = и; — = а

А а 2

формулы (16) и (19) принимают вид

р = 2— • аг^(и • х0) . (23)

а = аг^ [ю — й • tg (др)\. (24)

Определяя по формулам по формулам (23) и ёр ёа

(24) производные ----- и-----, а по формулам (13)

ёх0 ёр

векторные производные

ё1 ё1 т, ёЬ

У =-----; п' =------; Ь' =-----,

ёх0 ёх0 ёх0

после подстановки их с учетом формул (11) и первой из формул (13) в итоговую формулу (21), получим в координатной форме параметрические уравнения стрикции развертывающейся поверхности, образованной на основе линии (11). Ввиду громоздкости этих уравнений, которые также получены при помощи Мар1е-программы, они не приводятся. На представленной на рис. 2 визуализации отсека полученной поверхности (20) исходная кривая (11) и стрикция (кривая, касательная ко всем образующим прямым линиям поверхности) изображены как кривые, выходящие за пределы отсека поверхности.

Рис. 2. Визуализация развертывающейся поверхности

В заключение отметим возможность обобщения рассмотренной геометрической модели для образования, как косых, так и собственно развертывающихся линейчатых поверхностей. Как известно, прямая линия в пространстве определяется четырьмя независимыми параметрами, например, а , в, У и 3 Введем в пространстве сопровождающего трехгранника кривой некоторую группу линейных преобразований, например, шести параметрическую группу движений. Для одного из движений этой группы установим зависимость его параметров от лонгального параметра 5 кривой. В этом случае получаем возможность непрерывно изменять параметры положения а($), в(я), У(я) и 3(?) образующей прямой линии относительно сопровождающего трехгранника кривой линии, то есть в каждом фиксированном положении трехгранника будем получать фиксированное в нем положение прямой линии и разным положениям

/

/

/

трехгранника будут взаимно однозначно соответ- сти позволяет получать развертывающиеся по-ствовать разные в нем положения прямой линии. верхности в рассматриваемой обобщенной геомет-Введение дополнительно условия развертываемо- рической модели.

СПИСОК ЛИТЕРАТУРЫ

1. Рашевский, П. К. Курс дифференциальной геометрии. - М.: Гостехиздат, 1956. - 420 с.

2. Нитейский, А. С. Конструирование торсовой поверхности методом подвижного трехгранника Френе // Омский научный вестник. - 2013. - № 2(120). - С. 151-153.

3. Двайт, Г. Б. Таблицы интегралов и другие математические формулы. М.: Наука. - 1978. - 224с.

Авторы статьи

Панчук

Константин Леонидович; докт.техн.наук, доц., зав. каф. проф. каф. «Инженерная геометрия и САПР» (Омский государственный технический университет), e-mail: Panchuk [email protected]

Нитейский Антон Сергеевич; аспирант каф.«Инженерная геометрия и САПР» (Омский государственный технический университет), e-mail: [email protected].