Рус Eng Cn Перевести страницу на:  
Please select your language to translate the article


You can just close the window to don't translate
Библиотека
ваш профиль

Вернуться к содержанию

Программные системы и вычислительные методы
Правильная ссылка на статью:

Построение тела слоистой структуры с граничными поверхностями Безье

Аюшеев Тумэн Владимирович

доктор технических наук

доцент, кафедра Инженерная и компьютерная графика, Восточно-Сибирский государственный университет технологий и управления

670013, Россия, республика Бурятия, г. Улан-Удэ, ул. Клюквенная, 35

Ayusheev Tumen Vladimirovich

Doctor of Technical Science

Associate Professor, Department of Engineering and Computer Graphics, East-Siberian State University of Technology and Management

670013, Russia, respublika Buryatiya, g. Ulan-Ude, ul. Klyukvennaya, 35

atv62@bk.ru
Булычев Роман Николаевич

аспирант, кафедра Инженерная и компьютерная графика, Восточно-Сибирский государственный университет технологий и управления

670013, Россия, республика Бурятия, г. Улан-Удэ, ул. Ключевская, 40 в, оф. строение 1

Bulychev Roman Nikolaevich

graduate student, Department of Engineering and Computer Graphics, East-Siberian State University of Technology and Management

670013, Russia, respublika Buryatiya, g. Ulan-Ude, ul. Klyuchevskaya, 40 v, of. stroenie 1

buly4ev-roman@ya.ru
Мотошкин Петр Владимирович

кандидат технических наук

доцент, кафедра Инженерная и компьютерная графика, Восточно-Сибирский государственный университет технологий и управления

670013, Россия, республика Бурятия, г. Улан-Удэ, ул. 1-Я университетская, 11, оф. ДНТ ВСГУТУ

Motoshkin Petr Vladimirovich

PhD in Technical Science

Associate Professor, Department of Engineering and Computer Graphics, East-Siberian State University of Technology and Management

670013, Russia, respublika Buryatiya, g. Ulan-Ude, ul. 1-Ya universitetskaya, 11, of. DNT VSGUTU

mpv_mpv@mail.ru
Бубеев Иннокентий Трофимович

кандидат технических наук

доцент, кафедра Инженерная и компьютерная графика, Восточно-Сибирский государственный университет технологий и управления

670023, Россия, республика Бурятия, г. Улан-Удэ, ул. 40 Лет Победы, 23, кв. 1

Bubeev Innokentii Trofimovich

PhD in Technical Science

Associate Professor, Department of Engineering and Computer Graphics, East-Siberian State University of Technology and Management

670023, Russia, respublika Buryatiya, g. Ulan-Ude, ul. 40 Let Pobedy, 23, kv. 1

it_bubeev@mail.ru

DOI:

10.7256/2454-0714.2018.2.26291

Дата направления статьи в редакцию:

14-05-2018


Дата публикации:

25-05-2018


Аннотация: Объектом исследования является численное моделирование процесса инкрементальной формовки листового металла при изготовлении тонкостенных деталей на оборудовании с числовым программным управлением. Особое внимание уделяется задаче построения поверхности деформирования листовой заготовки, под воздействием рабочего инструмента. Предметом исследования является решение задачи построения гладкого трехпараметрического тела слоистой структуры с граничными поверхностями Безье по заданным узлам интерполяции. Тело строится по каркасу, состоящему из совокупности криволинейных восьмиугольных порций, у которых противоположные границы имеют различные параметрические длины. Методом исследования является способ построения трехпараметрического тела слоистой структуры с граничными поверхностями Безье с помощью применения технологии параметрического твердотельного моделирования. Новизна исследования заключается в разработке способа построения гладкого трехпараметрического тела слоистой структуры на произвольном каркасе в форме Безье. Основные выводы. Использование кубической функции Безье с параметризацией по длине дуги для описания тела многослойной конструкции дает дополнительные преимущества в управлении формой и внутреннего пространства тела, чем конструкция, построенная на применении обобщенной интерполяции Эрмита. Параметрические трехмерные тела такого типа могут быть эффективно использованы в конечных элементах с нерегулярными границами, чтобы избежать трудностей в построении модели для подобных конструкций.


Ключевые слова:

твердотельное моделирование, трехпараметрическое тело, порция тела, тело слоистой структуры, интерполяция, кривая Безье, поверхность Безье, послойное деформирование, инкрементальная формовка, многослойная конструкция

Abstract: The object of the study is a numerical simulation of the process of incremental forming of sheet metal in the manufacture of thin-walled parts on equipment with numerical program control. Particular attention is paid to the problem of constructing the surface of deformation of a sheet blank, under the influence of a working tool.The subject of the study is the solution of the problem of constructing a smooth three-parameter body with a layered structure using Bézier boundary surfaces from given interpolation nodes. The body is built on a frame, consisting of a set of curvilinear octagonal portions, in which the opposite boundaries have different parametric lengths.The method of investigation is a method for constructing a three-parameter body with a layered structure using Bezier boundary surfaces using the technology of parametric solid modeling. The novelty of the study is to develop a method for constructing a smooth three-parameter body of a layered structure on an arbitrary frame in the form of Bézier.Main conclusions. The use of the cubic Bézier function with parameterization along the length of the arc for describing the body of a multilayer structure provides additional advantages in controlling the shape and internal space of the body than the construction constructed using the generalized Hermite interpolation.Parametric three-dimensional bodies of this type can be effectively used in finite elements with irregular boundaries to avoid difficulties in constructing a model for such structures.


Keywords:

solid modeling, three-parameter body, body portion, lamellar body, interpolation, Bézier curve, Bézier surface, lamellar deformation, incremental forming, multilayer structure

Введение

При численном моделировании процесса инкрементальной формовки листового металла при изготовлении тонкостенных деталей на оборудовании с программным управлением возникает задача построения гладкой поверхности деформирования листовой заготовки, проходящей через заданные точки – узлы интерполяции [1]. При этом способе обработки материала давлением, в отличие от обычной штамповки, процесс формовки выполняется постепенно (послойно) в несколько этапов. Форма поверхности для каждого слоя задается определенным образом, чтобы обеспечить требуемую плавную деформацию листового материала без разрывов и заданной толщины. Эта поверхность образуется при движении рабочего инструмента по замкнутому контуру. В результате постепенного осаживания листовой заготовки под воздействием инструмента образуется поверхность слоя. При этом получаемые промежуточные поверхности слоев плавно меняются от одной граничной поверхности (плоская листовая заготовка) до другой (криволинейная поверхность изделия). Количество слоев зависит от требуемой точности получения внутреннего, или внешнего теоретического контура изделия. При численном моделировании такой многослойной конструкции используют поверхностные модели [2-4]. Существенным недостатком такого подхода является необходимость преобразования каркаса исходной поверхности при описании поверхности каждого слоя. Кроме того, с увеличением количества слоев растет объем перерабатываемой информации.

Для математического описания таких конструкций наиболее эффективным представляется использовать метод параметрического твердотельного моделирования. В этом методе тело представляется совокупностью ограничивающих его объем криволинейных восьмиугольных порций, граничные поверхности которых заданы параметрически. Многослойную конструкцию можно рассматривать как трехпараметрическое тело, где два параметра будут использованы для описания формы поверхности, а третий параметр – для описания толщины тела. Используя метод [5], можно построить уравнение порции трехпараметрического тела, обеспечивая необходимую непрерывность по производным на граничных поверхностях тела. Основной недостаток метода, основанного на применении порций с граничными бикубическими поверхностями Кунса в системе автоматизированного проектирования, связан с необходимостью задания касательных векторов и векторов смешанных производных, что может оказаться трудной задачей для проектировщика [6, 7].

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

В статье описывается метод построения уравнения порции тела для любой непрямоугольной области. Тело слоистой структуры можно строить по заданному каркасу, не преобразуя его и не размножая информации о нем.

Постановка задачи

Пусть нам задан массив узловых точек поверхности тела моделируемого объекта. Построена сетка кривых, определяющая поверхность тела. Сетка кривых делит поверхность на четырехугольные порции, которые ограничены u, v, w кривыми. Пусть параметрическая длина этих кривых изменяется в пределах от 0 до 1. Тогда r(u,v,w), 0<u,v,w<1, представляет внутреннее пространство порции тела, а r(i,v,w), r(u,j,w), r(u,v,k), i,j,k=0;1, определяют ее известные граничные поверхности (рис. 1).

Рисунок 1. Порция трехпараметрического тела

Необходимо определить функцию r(u,v,w), которая при u=i, v=j, w=k представляет нужную граничную поверхность порции тела, у которых противоположные границы имеют различные параметрические длины. При описании граничных кривых и поверхностей тела использовать параметрическую кубическую кривую Безье. Определить положение контрольных точек из условия непрерывности второй производной на границах порции тела.

Описание порции трехпараметрического тела

В общем случае трехмерный радиус-вектор порции трехпараметрического тела степени (3,3,3) можно записать в виде [8]

гдеbi(t) - базисные функции, Pijk- вершины характеристического многогранника.

Порция тела общего вида представляет собой криволинейный каркас, состоящий из 8 угловых точек, 12 граничных кривых и 6 граничных поверхностей (рис. 1). При описании граничных кривых и поверхностей порции тела многослойной конструкции в качестве базисных функций использовались функции Безье 3-степени. Рассмотрим их более подробно.

Описание граничной кривой порции тела

Формула кубической кривой Безье имеет вид [9]:

где Pi - вершины характеристической ломаной. Точки P0 и P3 являются концами кубической кривой Безье, точки P1 и P2 называют контрольными. Кривая Безье не проходит через контрольные точки, но они влияют на ее форму. В частности, прямые P0P1 и P2P3 являются касательными к кривой в ее концах (рис. 2).

Рисунок 2. Сегмент кубической кривой Безье

Формула (2) определяет кубическую интерполяционную функцию Безье на отрезке [0,1]. Говорят, что параметрическая длина этой кривой равна единице.

Для дальнейшего потребуется формула для кубической параметрической кривой Безье, у которой параметр u меняется от 0 до h, где h - расстояние между соседними узлами интерполяции. Введем преобразование s=uh. Тогда сегмент такой кривой можно записать так

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

Вычислим первую и вторую производные кривой Безье в ее концевых точках. Дважды дифференцируя по s равенство (3), получаем

Подставив значения s=0 и s=h в выражение (4), получим:

Введем следующие обозначения: узлы интерполяции обозначим через g0, g1,…, gn, кривую Безье, соединяющую узлы gi, и gi+1, обозначим через ri, а ее контрольные точки – через αi, βi, i=0,1,…,n-1.

Рассмотрим теперь ri – кривую Безье, соединяющую узлы интерполяции gi и gi+1, параметризованную значениями из отрезка [0,h], где hi =|gi+1-gi|. Обозначим в равенствах (5) и (6) контрольные точки через αi и βi, получим

Потребуем, чтобы

Используя выражения (7) и (8), получим

Сгруппируем члены с α и β в левую часть, а g в правую. Тогда получим следующую систему уравнений для определения αj, βj [10]:

Чтобы определить неизвестные величины αj, βj, надо задать еще два дополнительных условия, так называемые краевые условия. Потребуем, чтобы вторая производная в начальной и конечной точках равнялась нулю. Получаем уравнения

которые через точки записываются как

Система уравнений (9), (10) вместе с (11) и (12) дает возможность определить контрольные точки кубического параметрического сплайна Безье. На каждом участке тогда сплайн будет определен и его можно рассчитывать по формуле (3).

Используя формулу (3), запишем уравнения граничных кривых порции тела:

Описание граничной поверхности порции тела

Уравнение порции бикубической поверхности Безье имеет вид [9]:

где Pij – вершины характеристического многогранника (рис. 3).

Уравнения граничных поверхностей порции тела имеют следующий вид:

Рисунок 3. Порция бикубической поверхности Безье

Описание порции тела для произвольного каркаса

Уравнение порции трехпараметрического тела (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. Кim 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. 6.Фокс А., Пратт М. Вычислительная геометрия, применение в проектировании и на производстве. – М.: Мир, 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.
References
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. Ayusheev T.V. Metod postroeniya sploshnykh tel s primeneniem obobshchennoi interpolyatsii Ermita // Informatsionnye tekhnologii. 2005. №
6. S. 27-32. 6.Foks A., Pratt M. Vychislitel'naya geometriya, primenenie v proektirovanii i na proizvodstve. – M.: Mir, 1982. – 304 s.
7. Li K. Osnovy SAPR (CAD/CAM/CAE). – SPB.: Piter, 2004. – 560 s.
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. Golovanov N.N. Geometricheskoe modelirovanie. M.: Izd-vo fiziko-matematicheskoi literatury, 2002. – 472 s.
10. Borisenko V.V. Postroenie optimal'nogo splaina Bez'e // Fundamental'naya i prikladnaya matematika. 2016. Tom 21. № 3. S. 57-72.