Москва, г. Москва и Московская область, Россия
Москва, г. Москва и Московская область, Россия
ГРНТИ 55.01 Общие вопросы машиностроения
ГРНТИ 55.13 Технология машиностроения
Предложены обобщенные алгоритмы моделирования технических поверхностей сложной формы с применением метода аксонометрии, а также моделирования технологических процессов намотки.
композиционные материалы, цифровые пространственные модели, аксонометрия, ортогональная система координат, геометрия поверхности, кривизна геодезической параллели
Введение
В настоящее время в самых разных областях промышленного производства широкое применение все чаще находят композиционные материалы, или композиты. Арматурой в подобных материалах могут быть волокна различной формы, нити, ленты и т.п. Сложность изготовления изделий из них требует создания математического аппарата для определения внутренней геометрии поверхности объекта, на которую производится намотка, и соответствующего математического аппарата для самого процесса намотки. В данной работе решение основывается на привлечении цифровой пространственной модели, создаваемой с помощью метода аксонометрии, что позволяет иметь системную логическую структуру решения поставленной задачи.
Схематический анализ решения задачи получения геометрических моделей процессов намотки изделий из композиционных материалов
Проектированию сложных поверхностей из композиционных материалов, армированных волокнами с применением процессов намотки, посвящены работы Нартовой Л.Г., Якунина В.И. [1], Аюшеева Н.Н. [2; 3], Битюкова Ю.И., Калинина В.А. [4-6], Павловой С.В. [2; 3]. В них были разработаны геометрические модели процесса намотки, учитывающие реальную структуру ленты, а в работах Аюшеева Т.В. [1; 2] рассмотрено также применение систем технического зрения для повышения точности реализации процесса намотки. В работе Битюкова Ю.И. представлен математический аппарат для геометрического моделирования технологических процессов намотки и выкладки с разложением кубического сплайна по базисным сплайнам на равномерной сетке. Данная задача коррелирует с задачами технологии процесса намотки исходя из условий обеспечения прочности конструкции и т.п. Все это отражено на схеме соответствующей компьютерной модели технологических процессов намотки (рис. 1).
Получение такой модели носит естественный итерационный характер, который определяется требованиями технического задания на создание изделия.
При решении данной задачи возникают различные варианты решения. Например, необходим выбор условия пересечения ленты при намотке (рис. 2) и схемы контакта лент (рис. 3).
В настоящей работе предлагается обобщенный метод геометрического обеспечения процесса намотки, основанный на привлечении метода аксонометрии. В известной литературе этому уделено мало внимания, хотя полученная методика и является основополагающей при решении таких задач.
Моделирование технических поверхностей сложной формы на основе метода аксонометрии (реконструированной аксонометрической цифровой модели)
В качестве начального условия задана цифровая интерпретация точечных преобразований пространственного континуума, основанная на линейных преобразованиях. Эти преобразования задаются в общем виде функциями
х¢=f(х,у,z), y¢=j(х,у,z), z¢=Y(х,у,z).
В зависимости от аналитической структуры этих функций получаются преобразования различных классов. Например, если преобразующие функции имеют вид
то получается параллельная аксонометрия, обозначенная как результат аффинных (линейных) преобразований элементов пространства. На ее основе предложен новый вычислительный метод для практической реализации алгоритма построения параллельной аксонометрии (рис. 4). Даны натуральная прямоугольная система координат Охуz и фигура Ф, жестко связанная с этой системой. Оси координат совпадают с направлениями основных измерений фигуры Ф, называемой оригиналом.
Заданы также натуральные и аксонометрические масштабы измерений вдоль осей Ох, Оу, Оz - ex, ey, ez, ex¢, ey¢, ez¢. Схема построения вытекает из рис. 4. Это параллельная аксонометрия. В ней выполняются точные и вполне определяемые метрические соотношения. Так, существует метрическое соотношение, характеризующее любую параллельную аксонометрическую систему:
u2+ v2+ w2=2+ctg2j, (1)
где u=
v=
w=
j=(sÙП¢) - угол, составленный направлением проецирования s и плоскостью аксонометрических проекций П¢.
Соотношение (1) включает в себя в качестве частного вида прямоугольную аксонометрию (j = 900), т.е.
u2+ v2+ w2=2. (2)
Здесь уместно отметить основные свойства аксонометрической проекции. Аксонометрические координаты точки А, измеренные аксонометрическими масштабными единицами, численно всегда равны натуральным координатам точки А.
Теоретические изыскания в области аксонометрических систем разных классов содержат ответ и на главный вопрос - о расположении аксонометрических осей. Решение обеспечивает наиболее общая формулировка теоремы Польке - Шварца: всякий невырождающийся полный четырехугольник можно рассматривать как параллельную проекцию тэтраэдра наперед заданной формы. В эти исследования наибольший вклад внесли советские геометры: Н.Ф. Четверухин, Е.А. Глазунов, Н.М. Бескин.
Теперь более упрощенная формулировка общей теоремы звучит так: в косоугольной аксонометрии аксонометрические оси на плоскости чертежа и единичные отрезки на них могут быть выбраны совершенно произвольно. Это означает, что, взяв на плоскости П¢ три проходящие через одну точку несовпадающие прямые Ox¢, Oy¢, Oz¢ и отложив на них три отрезка произвольной длины (конечной и отличной от нуля), следует утверждать, что полученная фигура может рассматриваться как некоторая параллельная проекция трех взаимно перпендикулярных осей координат Охуz с отложенными на них соответственно равными единичными отрезками ОЕх=ОЕу=ОЕz=е. Из этой теоремы вытекает чрезвычайно важное следствие: любая аксонометрическая система (Ох¢у¢z¢, еx¢, еy¢, еz¢) в общем случае определяется пятью независимыми параметрами - тремя аксонометрическими единичными масштабами и двумя углами между аксонометрическими осями.
Теперь, после представления обзора фундаментальных свойств параллельной аксонометрии, можно дать цифровую характеристику процесса построений (схемы, алгоритмы) с привлечением этого метода. Это, по существу, означает, что все графические операции в рассматриваемом случае должны иметь соответствующую аналитическую интерпретацию, представленную в численных символах линейных отображений точек пространства.
Итак, пусть в пространстве задана натуральная ортогональная система координат Охуz. Сопутствующие построения выполнены на рис. 5, где направление проецирования обозначено вектором s, а плоскость проецирования (аксонометрическая плоскость проецирования П¢) задана уравнением
Ах+Ву+Сz+D=0. (3)
Проецирующий луч для точки О определен уравнением
Исходя из (3, 4) определяются координаты точки О¢(хо, уо, zо).
Аналогично, решая последовательно систему (3, 4) и используя заданные точки А, В, С, можно определить координаты точек О¢(хо¢, уо¢, zо¢), А¢(хА¢, уА¢, zА¢), В¢(хВ¢, уВ¢, zВ¢), С¢(хС¢, уС¢, zС¢). Показатели искажения u, v, w находятся из трех соотношений вида
Этот алгоритм позволяет определять необходимые параметры аксонометрического чертежа на плоскости П¢. Теперь способ построения цифровой модели заданных оригиналов сводится автоматически к нахождению проекций их точек – вычислению их аксонометрических координат. Такая задача решается на основании свойств линейных преобразований пространства, задаваемых системой
х¢=a1х+b1y+c1z+d1,
y¢=a2х+b2y+c2z+d2, (5)
z¢=a3х+b3y+c3z+d3.
В системе (5) присутствует 12 существенных коэффициентов, однозначно ее определяющих. Для их вычисления следует задать четыре пары соответствующих точек (О, А, В, С, О¢, А¢, В¢, С¢), положение которых в данном отображении известно.
Порядок действий в предлагаемом алгоритме состоит из следующих операций. Вводится в рассмотрение и решается система
хо¢=a10+b10+c10+d1,
yо¢=a20+b20+c20+d2, (6)
zо¢=a30+b30+c30+d3.
Следовательно, d1= хо0¢, d2= у0¢, d3= z0¢.
Аналогично для точек В, В¢, C, C¢ получается
b1= хB¢- х0¢; b2= yB¢- y0¢; b3= zB¢- z0¢;
c1= хc¢- х0¢; c2= yc¢- y0¢; c3= zc¢- z0¢.
Теперь окончательно соотношения, определяющие необходимые линейные преобразования, запишутся в виде
х¢=(х¢А-х¢0)х+(х¢В-х¢0)у+(х¢C-х¢0)z+x¢0,
у¢=(у¢А-у¢0)х+(у¢В-у¢0)у+(у¢C-у¢0)z+у¢0,
z¢=(z¢А-z¢0)х+(z¢В-z¢0)у+(z¢C-z¢0)z+z¢0
и могут быть реализованы известными вычислительными методами.
Дальнейшие теоретические изыскания связаны с результатами классической дифференциальной геометрии, позволяющей ввести метрику на поверхности и осуществить переход к соответствующим вычислительным моделям геометрических образов применительно к решению технических задач разного уровня. Это прежде всего первая (основная) квадратичная форма на заданной поверхности r=r(u,v)=ix(u, v)+jy(u,v)+kz(u,v), записанная в виде
dS2=E(u,v)du2+2F(u,v)dudv+G(u,v)dv2, (7)
где коэффициенты E, F, G - функции параметров u, v:
tutu=E(u,v), tutv=F(u,v), tvtv=G(u,v).
Выражение (7) имеет в теории поверхностей огромное значение, так как позволяет ввести метрику на аналитически заданной поверхности. Целесообразно заметить, что, зная первую квадратичную форму на поверхности, можно измерять не только длины дуг кривой, но и углы между ними (рис. 6). Для этого используется соотношение
Следовательно, угол j между координатами u и v в какой-то текущей точке Мi определяется по формуле
так как du¹0, dv=0 (смещение вдоль линии u), du=0, dv¹0 (смещение вдоль линии v).
Следствие из всего изложенного: если заданы E, F, G как функции u, v, то, даже если не задано уравнение поверхности, можно определить:
- длины дуг кривых;
- угол между ними;
- площади замкнутых контуров D областей (отсеков) поверхности, т.е.
Совокупность этих задач определяет внутреннюю геометрию поверхности. Понятие кривизны линии на поверхности связано с так называемой второй квадратичной формой поверхности
Ldu2+2Mdudv+Ndv2,
Если рассматривать величину ко как кривизну нормального сечения поверхности, то теперь она определяется по формуле
Моделирование технологических процессов намотки технических поверхностей сложной формы
Новый вычислительный метод приемлем применительно к технологическому процессу намотки поверхности, заданной в любой аналитической форме. Начальные данные позволяют реализовать основополагающие теоретические результаты в области классической дифференциальной геометрии и получить нужный алгоритм. Создана цифровая модель технологического процесса укладки, основанная на метрических свойствах геодезических кривых, выбранных на исходной поверхности.
Классические методы определения кривизны линии положены в основу моделирования технологических процессов намотки в системе авиастроения. Модель укладки на основную поверхность имеет свои размеры - ширину d и толщину h. Следует вывести алгоритм построения однонаправленных нитей относительно заданной кривой армирования (исходной кривой).
Предполагается, что основная кривая определяется уравнениями в параметрической форме rk(t)=
(uk(t),vk(t)), где известные функции uk(t),vk(t) задают изменение параметров u, v на поверхности оправки вдоль этой кривой:
u=uk(t), v=vk(t), to£t£tk .
В этом случае модель укладки с различными геометрическими характеристиками в поверхностном сечении содержит текущие точки, имеющие вполне определенный геометрический смысл. В данном случае они принадлежат геодезической линии на исходной поверхности. Базовое экспериментальное свойство геодезических линий состоит в том, что в малых окрестностях поверхности геодезические параллели всегда можно рассматривать как линии v некоторой полугеодезической системы (рис. 7). При этом отрезок геодезической линии АВ между точками А и В на ней дает кратчайшее расстояние между ними по сравнению со всеми другими кривыми АВ на поверхности, соединяющими эти точки, однако при условии, что точка В берется в достаточно малой окрестности точки А.
Следует заметить, что множество геодезических линий на поверхности ведет себя в известном смысле так, как множество прямых на плоскости: через каждую точку по каждому направлению можно провести геодезическую линию. Это вполне естественно, так как прямые на плоскости представляют собой частный случай геодезических линий на поверхности.
Итак, цифровое моделирование процесса укладки состоит в применении следующего алгоритма заданных последовательных операций (рис. 8). Средняя нить укладывается точно по основной кривой т (кривой армирования), а нить, отстоящая от средней на параметр d, будет укладываться относительно геодезической параллельной кривой армирования ортогонально ей. Для этого выбирается текущая точка M(t) для дуги 
на геодезической линии. Т.е. кривая намотки есть геометрическое место точек 
= (t,d), положение которых определяется выбранным параметром d. Если d=0, то получается исходная (начальная) кривая армирования. Итак, начальные условия организации процесса укладки включают в себя следующие компоненты: исходная поверхность Ф, исходная кривая т, параметр d на геодезической параллели, соответствующие параметры полосы (ленты).
Зная уравнение геодезической параллели относительно кривой армирования, соответствующей выбранному параметру d, можно описать детально весь рабочий процесс моделирования, обусловленный приложенным геометрическим аппаратом.
Математическое обеспечение процесса моделирования заключается в следующем. Выбирается произвольная текущая точка Mi на исходной кривой т; u=u(s), v=v(s) - функции, задающие кривую l на исходной поверхности, определенную уравнением t(s)=t(u(s), v(s)). В зависимости от знака s уравнения линии т на поверхности имеют вид
Здесь uD(0)=uk(t), vD(0)=vk(t).
Найдя решения u(s), v(s) из системы (8), удовлетворяющей начальным условиям (9), можно получить уравнение искомой линии т¢. Целесообразно заметить, что предложенной модели построения соответствует инструментальный метод организации управляющих программ, используемых при решении поставленной задачи.
Одним из основных параметров отклонения кривой т¢ на поверхности Ф, отвечающим за равновесное состояние нити, служит угол j. Так, для равновесия нити необходимо и достаточно выполнение условия útgjú£m, где m - коэффициент скольжения материала ленты вдоль поверхности Ф. Этот параметр можно определить из общей формулы классической внутренней геометрии поверхностей:
где kn= kn(t,d)= - кривизна геодезической параллели. Тогда неравенство útgjú£m обеспечивает равновесность текущих нитей m¢I и возможность определения из этого неравенства размеров ленты (útgjn(t,d)ú£m, tÎ[t0,tk], dÎ[ что, в свою очередь, обеспечивает оптимальность всего технологического процесса.
В случае использования полугеодезической системы координат имеет место соотношение ds2=du2+Gdv2 для любой ее дуги s. Ранее приведен алгоритм для реализации однослойной намотки. Для многослойной намотки теоретическим обоснованием ее осуществления могут служить некоторые канонические определения и свойства геометрии поверхностей. С этой целью вводится понятие гауссовой кривизны, выраженной через коэффициенты первой и второй квадратичных форм:
В данном случае имеется в виду слабодеформируемая лента m0. Известно, что при изгибании поверхности гауссова кривизна k в каждой ее точке остается постоянной, т.е. k=const. Следовательно, для каждого нового слоя этот признак остается в силе. Новый слой – это лента из множества ортогональных геодезических кривых, геодезических окружностей, алгоритм построения которых известен.
Заключение
Даются обобщенные алгоритмы цифрового моделирования технических поверхностей сложной формы с применением метода аксонометрии. Наличие такого рабочего аппарата позволяет представить общую математическую модель проектирования сложных поверхностей из композиционных материалов, армированных волокнами. Это позволяет ввести метрику на поверхности и осуществить переход к соответствующим вычислительным моделям геометрических образов сложных поверхностей. При создании цифрового метода технологического процесса укладки лент на соответствующую поверхность сложной формы приведен алгоритм процесса многослойной намотки с учетом гауссовой кривизны поверхности. Предложенная модель может служить методической основой при разработке новых конкретных изделий с использованием процесса намотки, а также при обучении студентов старших курсов специальным разделам САПР.
На основании изложенного можно утверждать, что особенностью поверхности постоянной кривизны является столь же полная однородность ее внутренней геометрии, как и для обыкновенной плоскости (имеется в виду ее самоналожение с тремя степенями свободы).
1. Нартова, Л.Г. Начертательная геометрия / Л.Г. Нартова, В.И. Якунин. - 2005. - С. 243.
2. Аюшеев, Т.В. Постановка задачи геометрического моделирования текстильной оболочки из композиционных материалов / Т.В. Аюшеев, С.В. Павлова, А.Ю. Булгатов // Проблемы меха-ники современных машин. - 2015. - Т. 3. - С. 79-84.
3. Павлова, С.В. К вопросу геометрического проектирования изделий со сложной формой поверхности / С.В. Павлова // Омский научный вестник. Серия «Приборы, машины и техноло-гии». - 2010. - № 1 (87). - С. 17-19.
4. Битюков, Ю.И. Численный анализ схемы укладки ленты переменной ширины на технологическую оправку в процессе намотки конструкций из композиционных материалов / Ю.И. Битюков, В.А. Калинин // Механика композиционных материалов и конструкций. - 2010. - Т. 16. - № 2. - С. 276-290.
5. Битюков, Ю.И. Расчет характеристик схемы укладки ленты в процессе намотки оболочек из композиционных материалов при его геометри-ческом моделировании с помощью гладкого отображения прямоугольника в пространство / Ю.И. Битюков // Омский научный вестник. - 2010. - № 2 (90). - С. 128-131.
6. Битюков, Ю.И. Представление кубического сплайна в виде разложения по В-сплайнам на равномерной сетке с явно выписанными коэф-фициентами разложения / Ю.И. Битюков // Ес-тественные и технические науки. - 2010. - № 3. - С. 304-308.
