CC BY-NC
9Построение тела слоистой структуры с граничными поверхностями Безье
Аюшеев Тумэн Владимирович
доктор технических наук
доцент, кафедра Инженерная и компьютерная графика, Восточно-Сибирский государственный
университет технологий и управления
670013, Россия, республика Бурятмя, г. Улан-Удэ, уп. Клюквенная, 35
И [email protected] Булычев Роман Николаевич
аспирант, кафедра Инженерная и компьютерная графика, Восточно-Сибирский государственный
университет технологий и управления
670013, Россия, республика Бурятмя, г. Улан-Удэ, уп. Ключевская, 40 в, оф строение 1
Мотоикин Петр Владимирович
кандидат технических наук
доцент, кафедра Инженерная и компьютерная графика, Восточно-Сибирский государственный
университет технологий и управления
670013, Россия, республика Бурятмя, г. Улан-Удэ, уп. 1-Я университетская, 11, оф. ДНТ ВСГУТУ
Бубеев 1/Ьннокектий Трофимович
кандидат технических наук
доцент, кафедра Инженерная и компьютерная графика, Восточно-Сибирский государственный
университет технологий и управления
670023, Россия, республика Бурятмя, г. Улан-Удэ, ул. 40 Лет Победы, 23, кв. 1
Статья из рубрики "Компьютерная графика, обработка изображений и распознавание образов"
Аннотация.
Объектом исследования является численное моделирование процесса инкрементальной формовки листового металла при изготовлении тонкостенных деталей на оборудовании с числовым программным управлением. Особое внимание уделяется задаче построения поверхности деформирования листовой заготовки, под воздействием рабочего инструмента. Предметом исследования является решение задачи построения гладкого трехпараметрического тела слоистой структуры с граничными поверхностями Безье по заданным узлам интерполяции. Тело строится по каркасу, состоящему из совокупности криволинейных восьмиугольных порций, у которых противоположные границы имеют различные параметрические длины. Методом исследования является способ построения трехпараметрического тела слоистой структуры с граничными поверхностями Безье с помощью применения технологии параметрического твердотельного моделирования. Новизна исследования заключается в разработке способа построения гладкого трехпараметрического тела слоистой структуры на произвольном каркасе в форме
Безье.Основные выводы. Использование кубической функции Безье с параметризацией по длине дуги для описания тела многослойной конструкции дает дополнительные преимущества в управлении формой и внутреннего пространства тела, чем конструкция, построенная на применении обобщенной интерполяции Эрмита.Параметрические трехмерные тела такого типа могут быть эффективно использованы в конечных элементах с нерегулярными границами, чтобы избежать трудностей в построении модели для подобных конструкций.
Ключевые слова: твердотельное моделирование, трехпараметрическое тело, порция тела, тело слоистой структуры, интерполяция, кривая Безье, поверхность Безье, послойное деформирование, инкрементальная формовка, многослойная конструкция
DOI:
10.7256/2454-0714.2018.2.26291
Дата направления в редакцию:
16-05-2018
Дата рецензирования:
24-05-2018
Дата публикации:
25-05-2018
Введение
При численном моделировании процесса инкрементальной формовки листового металла при изготовлении тонкостенных деталей на оборудовании с программным управлением возникает задача построения гладкой поверхности деформирования листовой заготовки,
проходящей через заданные точки - узлы интерполяции -Ш. При этом способе обработки материала давлением, в отличие от обычной штамповки, процесс формовки выполняется постепенно (послойно) в несколько этапов. Форма поверхности для каждого слоя задается определенным образом, чтобы обеспечить требуемую плавную деформацию листового материала без разрывов и заданной толщины. Эта поверхность образуется при движении рабочего инструмента по замкнутому контуру. В результате постепенного осаживания листовой заготовки под воздействием инструмента образуется поверхность слоя. При этом получаемые промежуточные поверхности слоев плавно меняются от одной граничной поверхности (плоская листовая заготовка) до другой (криволинейная поверхность изделия). Количество слоев зависит от требуемой точности получения внутреннего, или внешнего теоретического контура изделия. При численном
моделировании такой многослойной конструкции используют поверхностные модели [2-4]. Существенным недостатком такого подхода является необходимость преобразования каркаса исходной поверхности при описании поверхности каждого слоя. Кроме того, с увеличением количества слоев растет объем перерабатываемой информации.
Для математического описания таких конструкций наиболее эффективным
представляется использовать метод параметрического твердотельного моделирования. В этом методе тело представляется совокупностью ограничивающих его объем криволинейных восьмиугольных порций, граничные поверхности которых заданы параметрически. Многослойную конструкцию можно рассматривать как трехпараметрическое тело, где два параметра будут использованы для описания формы
поверхности, а третий параметр - для описания толщины тела. Используя метод можно построить уравнение порции трехпараметрического тела, обеспечивая необходимую непрерывность по производным на граничных поверхностях тела. Основной недостаток метода, основанного на применении порций с граничными бикубическими поверхностями Кунса в системе автоматизированного проектирования, связан с необходимостью задания касательных векторов и векторов смешанных производных, что может оказаться трудной задачей для проектировщика -6, 7].
Разработан ряд методов, позволяющих преодолеть эту трудность. Наиболее удобным является метод, основанный на применении кривой и поверхности в форме Безье. Наиболее широко используют бикубические поверхности Безье. Это позволяет получать составные поверхности требуемой гладкости при минимуме затрат на вычисление координат и производных в них, а также хранение структуры определяющих форму параметров.
В статье описывается метод построения уравнения порции тела для любой непрямоугольной области. Тело слоистой структуры можно строить по заданному каркасу, не преобразуя его и не размножая информации о нем.
Постановка задачи
Пусть нам задан массив узловых точек поверхности тела моделируемого объекта. Построена сетка кривых, определяющая поверхность тела. Сетка кривых делит поверхность на четырехугольные порции, которые ограничены и , V , м кривыми. Пусть параметрическая длина этих кривых изменяется в пределах от 0 до 1. Тогда г(ц ,у ), 0<и ,у <1, представляет внутреннее пространство порции тела, а г(/ ,у ), г(ц j ), ^и ,у ,к ), / ,к =0;1, определяют ее известные граничные поверхности (рис. 1).
г(1.0.0)
Рисунок 1. Порция трехпараметрического тела
Необходимо определить функцию r(u ,у ), которая при и =/ , V , м =к представляет нужную граничную поверхность порции тела, у которых противоположные границы имеют различные параметрические длины. При описании граничных кривых и поверхностей
тела использовать параметрическую кубическую кривую Безье. Определить положение контрольных точек из условия непрерывности второй производной на границах порции тела.
Описание порции трехпараметрического тела
В общем случае трехмерный радиус-вектор порции трехпараметрического тела степени (3,3,3) можно записать в виде
где Ь/ (£ )
многогранника.
базисные функции, Р ijk
вершины характеристического
Порция тела общего вида представляет собой криволинейный каркас, состоящий из 8 угловых точек, 12 граничных кривых и 6 граничных поверхностей (рис. 1). При описании граничных кривых и поверхностей порции тела многослойной конструкции в качестве базисных функций использовались функции Безье 3-степени. Рассмотрим их более подробно.
Описание граничной кривой порции тела
Формула кубической кривой Безье имеет вид
где Р/ - вершины характеристической ломаной. Точки Ро и Р3 являются концами кубической кривой Безье, точки Р1 и Р2 называют контрольными. Кривая Безье не проходит через контрольные точки, но они влияют на ее форму. В частности, прямые Р оР 1 и Р 2Р 3 являются касательными к кривой в ее концах (рис. 2).
Рисунок 2. Сегмент кубической кривой Безье
Формула (2) определяет кубическую интерполяционную функцию Безье на отрезке [0,1]. Говорят, что параметрическая длина этой кривой равна единице.
Для дальнейшего потребуется формула для кубической параметрической кривой Безье, у которой параметр и меняется от 0 до h , где h - расстояние между соседними узлами интерполяции. Введем преобразование s = uh . Тогда сегмент такой кривой можно записать так
к V я
и2 \ Щ - #
(3)
Теперь решим задачу проведения гладкой кривой Безье через заданные узлы интерполяции, путем их склеивания по второй производной.
Вычислим первую и вторую производные кривой Безье в ее концевых точках. Дважды дифференцируя по s равенство (3), получаем
Подставив значения s =0 и s ^ в выражение (4), получим:
Введем следующие обозначения: узлы интерполяции обозначим через д0, д1,..., дл , кривую Безье, соединяющую узлы д/ , и д/ +1, обозначим через г/ , а ее контрольные точки - через а/ , в/ , / =0,1,...,л -1.
Рассмотрим теперь г/ - кривую Безье, соединяющую узлы интерполяции д/ и д/ +1,
параметризованную значениями из отрезка -0,|11, где hi =|д/ +1-д/ |. Обозначим в равенствах (5) и (6) контрольные точки через а/ и в/ , получим
Потребуем, чтобы
Используя выражения (7) и (8), получим
Сгруппируем члены с а и в в левую часть, а д в правую. Тогда получим следующую систему уравнений для определения а' , в'-^:
Чтобы определить неизвестные величины а' , в' , надо задать еще два дополнительных
условия, так называемые краевые условия. Потребуем, чтобы вторая производная в начальной и конечной точках равнялась нулю. Получаем уравнения
которые через точки записываются как
Система уравнений (9), (10) вместе с (11) и (12) дает возможность определить контрольные точки кубического параметрического сплайна Безье. На каждом участке тогда сплайн будет определен и его можно рассчитывать по формуле (3).
Используя формулу (3), запишем уравнения граничных кривых порции тела:
Описание граничной поверхности порции тела
Уравнение порции бикубической поверхности Безье имеет вид [9]:
где Р у - вершины характеристического многогранника (рис. 3). Уравнения граничных поверхностей порции тела имеют следующий вид:
Рисунок 3. Порция бикубической поверхности Безье
г(0: р) =
3 А м 2 Зг1 м
—
V ¿00, 4» 1 ¿и. ¿00
Ро» Рди Р«п Рооа
Рдш Р;п Рлз
р р р р Г0Щ
Роза Р(Ш о? Рч
Описание порции тела для произвольного каркаса
Уравнение порции трехпараметрического тела (1) для произвольного каркаса в матричном виде будет иметь следующий вид:
Результаты экспериментов
Вычислительный эксперимент показал, что предложенная методика моделирования трехпараметрического для произвольного каркаса может быть использована при решении поставленной прикладной задачи. В рассмотренной модели определяются положения внешних и внутренних поверхностей тела, как при начальном условии, так и при деформации формы с помощью изменения положения вершин характеристического многогранника. Наличие двух промежуточных (контрольных) поверхностей между двумя противоположными граничными поверхностями тела дает возможность дополнительно управлять формой и внутренней плотностью тела, не изменяя при этом положения угловых точек каркаса тела. При описании тел составной формы кубическими функциями обеспечивается требуемый порядок гладкости при математическом описании многослойных конструкций, получаемых методом инкрементальной формовки из листового материала.
Для проведения вычислительного эксперимента были написаны подпрограммы с использованием встроенного языка программирования математического пакета программ MathCAD. Пример построения порции тела слоистой структуры представлен на рис. 4.
Рисунок 4. Пример построения порции тела слоистой структуры Заключение
Разработаны подпрограммы в среде MathCAD, использующие параметрические тела в форме Безье для моделирования внутренней многослойной структуры и внешней формы трехпараметрического тела. Использование кубической функции Безье с параметризацией по длине дуги для описания тела дает дополнительные преимущества в управлении формой и внутреннего пространства тела, чем конструкция, построенная на применении обобщенной интерполяции Эрмита. Проведенные примеры вычислительных экспериментов показывают, как легко определить и изменить форму и положение промежуточных поверхностей внутри тела.
Авторы рассматривают преимущество этого метода как способность определять граничные и внутренние начальные условия перед численным моделированием многослойных конструкций, получаемых различными технологическими способами, в частности, при инкрементальной формовке листового металла, под воздействием рабочего инструмента. Параметрические трехмерные тела такого типа могут быть эффективно использованы в конечных элементах с нерегулярными границами, чтобы избежать трудностей в построении модели для подобных конструкций.
Библиография
1. Micari F., Ambrogio G., Filice L. Shape and dimensional accuracy in single point incremental Forming: state of the art and future trends // Journal of Materials Processing Technology. 2007. V. 191(1-3). P. 390-395.
2. Guangcheng Z., Jinbo X., Xiaofan S., Xun Z., Chuankai L. Forming process of automotive body panel based on incremental forming // Metallurgical and Mining Industry. 2015. V. 1. I. 12. P. 350-357.
3. Kim T.J., Yang D.Y. Improvement of formability for the incremental sheet metal forming process // International Journal of Mechanical Sciences.2000.V. 42. I. 7. P.1271-1286.
4. Lanouar B.A., CamilleR., Arnaud D., Mohammed N. Simplified numerical approach for incremental sheet metal forming process // Engineering Structures. 2014. V. 62-63. P. 75-86.
5. Аюшеев Т.В. Метод построения сплошных тел с применением обобщенной интерполяции Эрмита // Информационные технологии. 2005. №
6. С. 27-32. б.Фокс А., Пратт М. Вычислительная геометрия, применение в проектировании и на производстве. - М.: Мир, 1982. - 304 с.
7. Ли К. Основы САПР (CAD/CAM/CAE). - СПБ.: Питер, 2004. - 560 с.
8. Durikovic R., Czanner S. Modelling with three types of Coons Bodies // International Journal of Modelling & Simulation. 2004. Vol. 24, № 2, P. 97-101.
9. Голованов Н.Н. Геометрическое моделирование. М.: Изд-во физико-математической литературы, 2002. - 472 с.
10. Борисенко В.В. Построение оптимального сплайна Безье // Фундаментальная и прикладная математика. 2016. Том 21. № 3. С. 57-72.
