Интервальный конечный элемент - Interval finite element

Максимальное напряжение по Мизесу в задаче о плоских напряжениях с интервальными параметрами (рассчитывается градиентным методом).

В численный анализ, то интервальный метод конечных элементов (интервал МКЭ) это метод конечных элементов который использует параметры интервала. Интервальный МКЭ может применяться в ситуациях, когда невозможно получить достоверные вероятностные характеристики конструкции. Это важно в бетонных конструкциях, деревянных конструкциях, геомеханике, композитных конструкциях, биомеханике и во многих других областях.[1] Цель интервального конечного элемента - найти верхнюю и нижнюю границы различных характеристик модели (например, стресс, смещения, поверхность текучести и т. д.) и использовать эти результаты в процессе проектирования. Это так называемая схема наихудшего случая, которая тесно связана с расчет предельного состояния.

Дизайн наихудшего случая требует меньше информации, чем вероятностный дизайн однако результаты более консервативны [Köylüoglu and Елисаков 1998].[нужна цитата ]

Применение интервальных параметров к моделированию неопределенности

Рассмотрим следующее уравнение:

где а и б находятся действительные числа, и .

Очень часто точные значения параметров а и б неизвестны.

Предположим, что и . В этом случае необходимо решить следующее уравнение

Существует несколько определений множества решений этого уравнения с интервальными параметрами.

Единый набор решений

При таком подходе решением является следующий набор

Это наиболее популярный набор решений интервального уравнения, и этот набор решений будет применен в этой статье.

В многомерном случае множество единых решений намного сложнее. Множество решений следующей системы линейные интервальные уравнения

показано на следующем рисунке

Решение set.png

Набор точных решений очень сложен, поэтому необходимо найти наименьший интервал, который содержит набор точных решений.

Набор решений-3.png

или просто

где

Смотрите также [1]

Набор параметрических решений интервальной линейной системы

Метод интервальных конечных элементов требует решения системы уравнений, зависящих от параметров (обычно с симметричной положительно определенной матрицей). Пример набора решений общей системы уравнений, зависящих от параметров.

показан на картинке ниже.[2]

Множество решений зависимой от параметров системы уравнений

Алгебраическое решение

В этом подходе x является номер интервала для которого уравнение

доволен. Другими словами, левая часть уравнения равна правой части уравнения. В этом частном случае решение имеет вид потому что

Если неопределенность больше, т.е. , тогда потому что

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

Метод

Рассмотрим PDE с интервальными параметрами

где - вектор параметров, принадлежащих заданным интервалам

Например, уравнение теплопередачи

где параметры интервала (т.е. ).

Решение уравнения (1) можно определить следующим образом

Например, в случае уравнения теплопередачи

Решение очень сложно из-за того, что на практике более интересно найти наименьший возможный интервал, который содержит точное множество решений .

Например, в случае уравнения теплопередачи

Метод конечных элементов приводит к следующей системе алгебраических уравнений, зависящей от параметров

где K это матрица жесткости и Q это правая часть.

Интервальное решение можно определить как многозначную функцию

В простейшем случае описанную выше систему можно рассматривать как систему линейные интервальные уравнения.

Также возможно определить интервальное решение как решение следующей задачи оптимизации

В многомерном случае внутреннее решение можно записать как

Интервальное решение против вероятностного решения

Важно знать, что параметры интервала дают разные результаты, чем равномерно распределенные случайные величины.

Параметр интервала учитывать все возможные распределения вероятностей (для ).

Для определения параметра интервала необходимо знать только верхний и нижняя граница .

Расчет вероятностных характеристик требует знания большого количества экспериментальных результатов.

Можно показать, что сумма n интервальных чисел равна раз больше, чем сумма соответствующих нормально распределенных случайных величин.

Сумма п номер интервала равно

Ширина этого интервала равна

Рассматривать нормально распределенная случайная величина Икс такой, что

Сумма п нормально распределенная случайная величина - это нормально распределенная случайная величина со следующими характеристиками (см. Шесть Сигм )

Можно считать, что ширина вероятностного результата равна 6 сигм (ср. Шесть Сигм ).

Теперь мы можем сравнить ширину интервала результата и вероятностного результата.

Из-за этого результаты интервального конечного элемента (или в целом анализа наихудшего случая) могут быть завышены по сравнению со стохастическим анализом бедренной кости (см. Также распространение неопределенности Однако в случае не вероятностной неопределенности невозможно применить чисто вероятностные методы. Поскольку вероятностные характеристики в этом случае точно не известны [ Елисаков 2000].

Можно рассматривать случайные (и нечеткие случайные величины) с интервальными параметрами (например, со средним интервалом, дисперсией и т. Д.). Некоторые исследователи используют интервальные (нечеткие) измерения в статистических расчетах (например, [2] ). В результате таких расчетов мы получим так называемые неточная вероятность.

Неточная вероятность понимается в очень широком смысле. Он используется как общий термин для обозначения всех математических моделей, которые измеряют вероятность или неопределенность без точных числовых вероятностей. Он включает как качественные (сравнительная вероятность, частичное упорядочение предпочтений,…), так и количественные режимы (интервальные вероятности, функции убеждений, верхнее и нижнее предвидение,…). Неточные вероятностные модели необходимы в задачах вывода, когда релевантная информация является скудной, расплывчатой ​​или противоречивой, а также в задачах принятия решений, где предпочтения также могут быть неполными. [3].

Простой пример: моделирование растяжения, сжатия, деформации и напряжения)

TensionCompression.JPG

1-мерный пример

в напряжение -сжатие проблема, следующие уравнение показывает взаимосвязь между смещение ты и сила п:

где L длина, А - площадь поперечного сечения, а E является Модуль для младших.

Если модуль Юнга и сила неопределенны, то

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

Рассчитайте экстремальные значения смещения следующим образом:

Рассчитать напряжение используя следующую формулу:

Вычислите производную деформации, используя производную от перемещений:

Рассчитайте экстремальные значения смещения следующим образом:

Также можно рассчитать предельные значения деформации, используя перемещения

тогда

Та же методология может быть применена к стресс

тогда

и

Если рассматривать стресс как функцию напряжения, тогда

тогда

Конструкция безопасна при стрессе меньше заданного значения т.е.

это условие верно, если

После расчета мы знаем, что это соотношение выполняется, если

Пример очень простой, но он показывает применение интервальных параметров в механике. Интервальные МКЭ используют очень похожую методологию в многомерных случаях [Pownuk 2004].

Однако в многомерных случаях связь между неопределенными параметрами и решением не всегда монотонна. В этом случае должны применяться более сложные методы оптимизации.[1]

Многомерный пример

В случае напряжения -сжатие уравнение равновесия имеет следующий вид:

где ты это смещение, E является Модуль для младших, А площадь поперечного сечения, а п является распределенной нагрузкой.Чтобы получить уникальное решение, необходимо добавить соответствующие граничные условия, например

Если Модуль для младших E и п неопределенны, то интервальное решение можно определить следующим образом

Для каждого элемента МКЭ уравнение можно умножить на тестовую функцию v

где

После интеграция по частям получим уравнение в слабой форме

где

Введем набор точек сетки , где - количество элементов, а функции линейной формы для каждого элемента FEM

где

левая конечная точка элемента, левая конечная точка элемента с номером «е». Приближенное решение в "е" -м элементе представляет собой линейную комбинацию функций формы

После подстановки в слабую форму уравнения получим следующую систему уравнений

или в матричной форме

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

где

- глобальная матрица жесткости,

- вектор решения,

это правая часть.

В случае проблемы растяжения-сжатия

Если пренебречь распределенной нагрузкой п

После учета граничных условий матрица жесткости имеет следующий вид

Правая часть имеет следующий вид

Предположим, что модуль Юнга E, площадь поперечного сечения А и нагрузка п неопределенны и принадлежат некоторым интервалам

Интервальное решение можно определить вычислением следующим образом

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

Результатом расчетов являются интервальные смещения

Предположим, что смещения в столбце должны быть меньше некоторого заданного значения (из соображений безопасности).

Неопределенная система безопасна, если интервальное решение удовлетворяет всем условиям безопасности.

В данном конкретном случае

или просто

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

Метод интервальных конечных элементов может применяться для решения задач, в которых недостаточно информации для создания надежной вероятностной характеристики конструкций [ Елисаков 2000]. Метод интервальных конечных элементов может быть также применен в теории неточная вероятность.

Метод комбинирования конечных точек

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

Метод комбинирования конечных точек дает решение, которое обычно бывает точным; к сожалению, этот метод имеет экспоненциальную вычислительную сложность и не может быть применен к задачам со многими интервальными параметрами [Neumaier 1990].

Метод расширения Тейлора

Функция можно расширить с помощью Серия Тейлор В простейшем случае ряды Тейлора используют только линейную аппроксимацию.

Верхнюю и нижнюю границы решения можно рассчитать по следующей формуле

Метод очень эффективен, но не очень точен.
Для повышения точности можно применить разложение Тейлора более высокого порядка [Pownuk 2004].
Этот подход также может быть применен в интервальный метод конечных разностей и метод интервальных граничных элементов.

Градиентный метод

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

если тогда
если тогда

Экстремальные значения решения можно рассчитать следующим образом

Во многих конструкторских приложениях этот метод дает точное решение.
Если решение не монотонное, оно обычно разумно. Для повышения точности метода можно применять тесты на монотонность и анализ чувствительности более высокого порядка. Метод может быть применен для решения линейных и нелинейных задач вычислительной механики [Pownuk 2004]. Применение метода анализа чувствительности к решению задач гражданского строительства можно найти в следующей статье [М.В. Рама Рао, А. Павунук и И. Скална 2008].
Этот подход также может быть применен в интервальный метод конечных разностей и метод интервальных граничных элементов.

Элемент за элементом метод

Муханна и Маллен применили поэлементную формулировку к решению уравнения конечных элементов с интервальными параметрами [Muhanna, Mullen 2001]. Используя этот метод, можно получить решение с гарантированной точностью в случае стропильных и каркасных конструкций.

Методы возмущений

Решение матрица жесткости и вектор нагрузки можно расширить с помощью теория возмущений. Теория возмущений приводит к приближенному значению интервального решения [Qiu, Елисаков 1998]. Метод очень эффективен и может быть применен к большим задачам вычислительной механики.

Метод поверхности отклика

Можно приблизить решение используя поверхность отклика. Затем можно использовать поверхность отклика для получения интервального решения [Akpan 2000]. Используя метод поверхности отклика, можно решить очень сложную задачу вычислительной механики [Beer 2008].

Чистые интервальные методы

Некоторые авторы пытались применить чистые интервальные методы к решению задач конечных элементов с интервальными параметрами. В некоторых случаях можно получить очень интересные результаты, например [Попова, Янков, Бонев 2008]. Однако в целом метод дает очень завышенные результаты [Kulpa, Pownuk, Skalna 1998].

Параметрические интервальные системы

[Popova 2001] и [Skalna 2006] ввели методы решения системы линейных уравнений, в которой коэффициенты являются линейными комбинациями интервальных параметров. В этом случае можно получить очень точное решение интервальных уравнений с гарантированной точностью.

Смотрите также

использованная литература

  1. ^ а б «Архивная копия». Архивировано из оригинал на 2011-10-05. Получено 2008-10-12.CS1 maint: заархивированная копия как заголовок (ссылка на сайт)
  2. ^ Попова Е. А. Набор параметрических решений интервальной линейной системы. В архиве 2010-01-27 на Wayback Machine
  • U.O. Акпан, Т. Коко, И. Орисамолу, Б.К. Галлант, Практический нечеткий анализ конструкций с помощью конечных элементов, Конечные элементы в анализе и проектировании, 38, стр. 93–111, 2000.
  • М. Бир, Оценка несовместимых инженерных данных, Третий семинар по надежным инженерным вычислениям (REC08) Технологический институт Джорджии, 20–22 февраля 2008 г., Саванна, Джорджия, США.
  • Демпстер, А. П. (1967). «Верхняя и нижняя вероятности, индуцированные многозначным отображением». Анналы математической статистики 38 (2): 325-339. [4]. Проверено 23 сентября 2009 г.
  • Анализ неопределенности в гражданском строительстве, В. Феллин, Х. Лессманн, М. Обергуггенбергер и Р. Вайдер (ред.), Springer-Verlag, Берлин, 2005 г.
  • И. Елисаков, Возможные ограничения вероятностных методов в технике. Обзоры прикладной механики, Том 53, № 2, стр. 19–25, 2000 г.
  • Главачек, И., Хлебоун, Дж., Бабушка, И.: Проблемы с неопределенными исходными данными и метод наихудшего сценария. Эльзевир, Амстердам (2004)
  • Кёйлюоглу, У., Исаак Елисаков; Сравнение стохастических и интервальных конечных элементов, применяемых к сдвиговым каркасам с неопределенными свойствами жесткости, Computers & Structures Volume: 67, Issue: 1-3, April 1, 1998, pp. 91–98
  • Кульпа З., Павнюк А., Скална И., Анализ линейных механических конструкций с неопределенностями с помощью интервальных методов. Компьютерная механика и инженерные науки, т. 5. 1998. С. 443–477.
  • Д. Моенс и Д. Вандепитт, Теория интервальной чувствительности и ее применение к анализу огибающей частотной характеристики неопределенных структур. Компьютерные методы в прикладной механике и технике Vol. 196, № 21-24, 1 апреля 2007 г., стр. 2486–2496.
  • Мёллер Б., Бир М. Нечеткая случайность - неопределенность в гражданском строительстве и вычислительной механике, Springer, Berlin, 2004.
  • Р.Л. Муханна, Р.Л. Маллен, Неопределенность в задачах механики - интервальный подход. Журнал инженерной механики, Том 127, № 6, 2001 г., 557-556
  • А. Ноймайер, Интервальные методы для систем уравнений, Cambridge University Press, Нью-Йорк, 1990.
  • Попова Е. О решении параметризованных линейных систем. W. Kraemer, J. Wolff von Gudenberg (Eds.): Научные вычисления, Проверенные числа, интервальные методы. Kluwer Acad. Издательство, 2001, с. 127–138.
  • Е. Попова, Р. Янков, З. Бонев: Ограничение отклика механических конструкций с неопределенностями по всем параметрам. В R.L.Muhannah, R.L.Mullen (Eds): Proceedings of the NSF Workshop on Reliable Engineering Computing (REC), Svannah, Georgia USA, 22-24 февраля 2006 г., стр. 245-265
  • A. Pownuk, Численное решение нечетких дифференциальных уравнений в частных производных и их применение в вычислительной механике, Нечеткие дифференциальные уравнения с частными производными и реляционные уравнения: характеристика и моделирование коллектора (М. Никравеш, Л. Заде и В. Коротких, ред.), Исследования в области нечеткости и мягкие вычисления, Physica-Verlag, 2004, стр. 308–347.
  • А. Павунук, Эффективный метод решения крупномасштабных инженерных задач с интервальными параметрами на основе анализа чувствительности, Материалы семинара NSF по надежным инженерным вычислениям, 15–17 сентября 2004 г., Саванна, Джорджия, США, стр. 305–316
  • М.В. Рама Рао, А. Павунук и И. Скална, Расчет напряжений монолитной железобетонной балки с неопределенными структурными параметрами, Семинар NSF по надежным инженерным вычислениям, 20–22 февраля 2008 г., Саванна, Джорджия, США, стр. 459–478
  • И. Скална, Метод внешнего интервального решения систем линейных уравнений, линейно зависящих от интервальных параметров, Надежные вычисления, Том 12, номер 2, апрель, 2006 г., стр. 107–120
  • З. Цю и И. Елисаков, Антиоптимизация конструкций с большими неопределенными, но неслучайными параметрами с помощью интервального анализа Компьютерные методы в прикладной механике и технике, Том 152, выпуски 3-4, 24 января 1998 г., страницы 361-372
  • Бернардини, Альберто, Тонон, Фульвио, Граничная неопределенность в гражданском строительстве, Springer 2010
  • Бен-Хаим Ю., Елисаков И., 1990, Выпуклые модели неопределенности в прикладной механике. Издательство Elsevier Science, Нью-Йорк
  • Валлиаппан С., Фам Т.Д., 1993, Нечеткий конечно-элементный анализ основания на упругой грунтовой среде. Международный журнал численных и аналитических методов в геомеханике, том 17, стр. 771–789
  • Елисаков И., Ли Ю. В., Старнес Дж. Х., 1994, Детерминированный метод прогнозирования влияния неизвестных, но ограниченных модулей упругости на изгиб композитных конструкций. Компьютерные методы в прикладной механике и технике, Том 111, с. 155–167.
  • Валлиаппан С. Фам Т.Д., 1995, Анализ упругопластических конечных элементов с нечеткими параметрами. Международный журнал численных методов в инженерии, 38, стр. 531–548.
  • Рао С.С., Сойер Дж. П., 1995, Нечеткий метод конечных элементов для анализа неточно определенных систем. Журнал AIAA, Том 33, № 12, стр. 2364–2370
  • Кёйлюоглу Х.У., Чакмак А., Нильсен С.Р.К., 1995, Отображение интервалов в строительной механике. В кн .: Спанос, под ред. Вычислительная стохастическая механика. 125-133. Балкема, Роттердам
  • Муханна, Р. Л. и Р. Л. Маллен (1995). "Разработка интервальных методов для нечеткости в механике сплошной среды" в материалах 3-го Международного симпозиума по моделированию и анализу неопределенности и ежегодной конференции Североамериканского общества обработки нечеткой информации (ISUMA – NAFIPS'95), IEEE, 705–710

внешние ссылки