О вариационной задаче для деформации параметрической поверхности

Формулировка проблемы достижения условия непрерывности G и описание соответствующих уравнений для решения этой задачи. Функционалы "сдвиг кривой" и Квази-G1. Решение вариационных задач без ограничений в соответствии с теоремой Ферма, описание алгоритма.

Рубрика Математика
Вид статья
Язык русский
Дата добавления 21.06.2018
Размер файла 89,3 K

Отправить свою хорошую работу в базу знаний просто. Используйте форму, расположенную ниже

Студенты, аспиранты, молодые ученые, использующие базу знаний в своей учебе и работе, будут вам очень благодарны.

Размещено на http://www.allbest.ru/

О вариационной задаче для деформации параметрической поверхности

Берзин Д.В.

Кандидат физико-математических наук, доцент, Финансовый университет при Правительстве Российской Федерации, Москва

Аннотация

Допустим, мы решили трансформировать поверхность минимальным образом, чтобы соблюсти определенные условия непрерывности с ее соседями. В настоящей работе сформулирована проблема достижения условия непрерывности G и выписаны соответствующие уравнения для решения этой задачи.

Ключевые слова: информационные технологии, условие непрерывности G, NURBS, оптимальное решение

Abstract

Suppose we want to deform a base surface of a face to achieve certain continuity condition (for example, G continuity) with the given neighbors at common edges. We formulate a problem for preserving G continuity and give equations for a deformation that changes the surface geometry as little as possible.

Keywords: Information technology, G continuity, NURBS, optimal solution.

Предположим, что грань F0 окружена некоторым количеством соседних граней F1, F2, …. Мы хотим деформировать (начальную) базовую поверхность грани F0 для того, чтобы достичь некоторого условия непрерывности (например, условия G ) с данными соседними поверхностями у общих ребер. Такая деформация должна минимально изменять геометрию поверхности.

1. Функционал “сдвиг кривой”

Обозначим векторы начальных и измененных (деформированных) контрольных точек как соответственно. Рассмотрим кривую c(t), которая принадлежит (или располагается рядом с) начальной (не деформированной) поверхности S(P).пусть

w = w(c, t) = (u(t),(v(t)) (1)

uv-кривая, соответствующая c (t). Рассмотрим класс трехмерных кривых с фиксированным w и переменной P:

(2)

Рассмотрим функционал

(3)

который выражает расстояние (или максимальный зазор) между начальной и деформированной кривыми. Назовем этот функционал “сдвиг кривой”.

2. Другие функционалы

Рассмотрим еще два типа функционалов: H(P) и G(P). Функционал H(P), который мы назовем “сдвиг контрольной точки” выражает расстояние между множеством контрольных точек P и P. H(P) служит для контроля отклонения деформированной поверхности от начальной. Функционал “непрерывности” G(P) служит для того, чтобы выполнялось условие непрерывности, например, условие G№ с соседними гранями.

3. Квази-G№

Вместо условия G№, мы можем попытаться достичь несколько другого (и в некотором смысле более строгого условия, чем G№), которое, тем не менее, приведет к линейности в вариационной задаче. Пусть E - это произвольная, но фиксированная точка на некотором ребре, которое является общим для грани F0 и соседней грани F1.

Рассмотрим плоскость р , которая является касательной в точке E к базовой поверхности грани F1. Пусть - соответствующие касательные векторы (исходящие из точки E и направленные вдоль u и v соответственно) к начальной базовой поверхности S(P) грани F0. Спроектируем на плоскость р, получим пару векторов соответственно. Теперь мы можем выписать функционал “непрерывности” для этого условия в точке E:

(4)

где - соответствующие касательные векторы к деформированной поверхности S(P), а б и в являются действительными переменными. Соответственно, функционал непрерывности для множества точек Щ - это

.

4. Вариационная задача

Теперь мы можем выписать функционал “общей ошибки”

(5)

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

F(P)> min (6)

Эта вариационная задача без ограничений (см. [1]) может быть решена в соответствии с теоремой Ферма:

(7)

где - это решение данной задачи.

5. Замечания

При нашем подходе, количество узлов и количество контрольных точек остаются теми же самыми после деформации. Возможно, это ограничение не позволит достичь точного условия непрерывности и сохранения граничных кривых в пределах заданных допусков. Необходимо измерить непрерывность и сдвиги кривых, и, если необходимо, вставить дополнительные узлы в исходную поверхность, и после этого начать деформацию снова.

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

Условие “квази-G№” - это, конечно, не то же самое, что G№, но мы ожидаем, что в общем (6) “заставит” соответствующие касательные плоскости приблизиться к желаемым положениям.

6. Алгоритм

Приведем ниже описание алгоритма для достижения нужной деформации (см. также [2]). Чтобы вычислить новые расположения контрольных точек, мы должны выполнить следующие основные шаги:

1) Составить вектор из (N+1) контрольных точек для начальной поверхности. ферма функционал алгоритм теорема

2) Отобрать (d+1) точку выборки, чтобы сохранить положения граничных кривых . Будем называть эти точки “G0 точки выборки”.

3) Отобрать (g+1) точку выборки, чтобы сохранить непрерывность . Будем называть эти точки “G1 точки выборки”. В нашем исполнении, множество G1 точек выборки является подмножеством G0 точек выборки.

4) Вычислим “желаемую” касательную плоскость в каждой G1 точке выборки и получим 2(g+1) соответствующих проекций касательных векторов к начальной поверхности на касательную плоскость: . Другими словами, получаем пару трехмерных векторов для каждой G1 точки выборки.

5) Вычислим два вектора, каждый из которых состоит из (g+1) константы: .

6) Вычислим (N+1)-вектор для каждой G0 точки выборки используя базисные функции B-сплайна.

7) Вычислим пару (N+1)-векторов для каждой G1 точки выборки используя базисные функции B-сплайна и их производные.

8) Вычислим три (N+1)Ч(N+1) матрицы: .

9) Вычислим (N+1)Ч(N+1) матрицу .

10) Нам больше не нужны матрицы и мы можем освободить соответствующую память.

11) Вычислим три вектора, каждый из которых состоит из (N+1) константы: .

12) Нам больше не нужна матрица и вектор контрольных точек , и мы можем освободить соответствующую память.

13) Вычислим шесть (N+1)Ч(g+1) матриц .

14) Вычислим шесть (g+1)Ч(N+1) матриц .

15) Освободим память, выделенную для каждой контрольной точки.

16) Составим матрицу A. Освободим соответствующую память.

17) Составим вектор констант правой части B.

18) Решим систему уравнений, используя методы для систем с разреженными матрицами.

19) Получим множество новых контрольных точек.

20) Построим новую поверхность.

Литература

1. S. V. Fomin, I. M. Gelfand “Calculus of Variations” // Dover Publications, 2000

2. Berzin D.V. “Surface deformation with geometric constraints” // Международныйнаучно-исследовательскийжурнал = Research Journal of International Studies, №8 (15) 2013, part 1, p.6

Размещено на Allbest.ru


Подобные документы

  • Форма для ввода целевой функции и ограничений. Характеристика симплекс-метода. Процесс решения задачи линейного программирования. Математическое описание алгоритма симплекс-метода. Решение задачи ручным способом. Описание схемы алгоритма программы.

    контрольная работа [66,3 K], добавлен 06.04.2012

  • Описание метода потенциалов Математическая постановка задачи об оптимальных перевозках. Метод решения задачи об оптимальных перевозках средствами Ms Excel. Постановка параметрической транспортной задачи, ее математическое и компьютерное моделирование.

    курсовая работа [802,5 K], добавлен 21.10.2014

  • Особенности решения задач Диофантовой "Арифметики", которые решаются с помощью алгебраических уравнений или системы алгебраических уравнений с целыми коэффициентами. Характеристика великой теоремы Ферма, анализ и методы приминения алгоритма Евклида.

    реферат [36,8 K], добавлен 03.03.2010

  • Описание метода сведения краевой задачи к задаче Коши. Решение системы из двух уравнений с четырьмя неизвестными. Метод Рунге-Кутта. Расчет максимальной погрешности и выполнение проверки точности. Метод конечных разностей. Описание полученных результатов.

    курсовая работа [245,2 K], добавлен 10.07.2012

  • Решение эллиптических и параболических дифференциальных уравнений в частных производных. Суть метода Кранка-Николсона и теории разностных схем для теплопроводности. Построение численных методов с помощью вариационных принципов, описание Matlab и Mathcad.

    курсовая работа [1,4 M], добавлен 13.03.2011

  • Структура текстовой задачи. Условия и требования задач и отношения между ними. Методы и способы решения задач. Основные этапы решения задач. Поиск и составление плана решения. Осуществление плана решения. Моделирование в процессе решения задачи.

    презентация [247,7 K], добавлен 20.02.2015

  • Понятие функционала и оператора. Задачи, приводящие к экстремуму функционала, и необходимые условия его минимума. Связь между вариационной и краевой задачами. Функционалы, зависящие от нескольких функций. Вариационные задачи с подвижными границами.

    курсовая работа [313,3 K], добавлен 23.05.2010

  • Сущность методов сведения краевой задачи к задаче Коши и алгоритмы их реализации на ПЭВМ. Применение метода стрельбы (пристрелки) для линейной краевой задачи, определение погрешности вычислений. Решение уравнения сшивания для нелинейной краевой задачи.

    методичка [335,0 K], добавлен 02.03.2010

  • Квази-средние как обобщение классических средних величин. Квази-средние и функциональные уравнения. Решение некоторых функциональных уравнений. Характеристическое свойство квази-средних. Квази-средние и выпуклые функции.

    дипломная работа [412,7 K], добавлен 08.08.2007

  • Понятие Диофантовых уравнений, их сущность и особенности, методика и этапы решения. Великая теорема Ферма и порядок ее доказательства. Алгоритм решения иррациональных уравнений. Метод поиска Пифагоровых троек. особенности решения уравнения Каталана.

    учебное пособие [330,2 K], добавлен 23.04.2009

Работы в архивах красиво оформлены согласно требованиям ВУЗов и содержат рисунки, диаграммы, формулы и т.д.
PPT, PPTX и PDF-файлы представлены только в архивах.
Рекомендуем скачать работу.