Инженерный анализ методом конечных элементов (мкэ)

6. Ошибки метода конечных элементов Критерии устойчивости, сходимости и точности в основном определяются погрешностями различного рода операций, проводимых в МКЭ. Наряду с обычными ошибками округления и погрешностью приближенных методов линейной алгебры, применяемых в МКЭ, есть и ошибки, имеющие непосредственное отношение к методу конечных элементов. Разбиение области на КЭ не является единственным. Зависимость расчета от выполняемого пользователем выбора (построения) сетки КЭ и трудность оценки точности получаемых результатов является основными недостатками метода. Погрешности метода конечных элементов связаны с: – ошибки дискретизации, являющиеся результатом различий между действительной геометрией рассчитываемой области и ее аппроксимацией системой конечных элементов; – ошибки аппроксимации, обусловленные разностью между действительным распределением искомых функций в пределах КЭ и их представлением с помощью аппроксимирующих функций.   Ошибки дискретизации уменьшаются с увеличением числа конечных элементов и соответственно с уменьшением их размеров, причем они стремятся к нулю, когда размер элемента стремится к нулю.  Эти ошибки уменьшаются    и с применением криволинейных элементов на соответствующих границах    области. Ошибки аппроксимации не обязательно уменьшаются по мере уменьшения размеров элементов или повышения степени аппроксимации, поэтому могут ухудшать сходимость к точному решению или даже приводить  к расходимости.   Однако общий метод оценки(универсальный и теоретически обоснованный) погрешности МКЭ на сегодня отсутствует, а точное решение в реальных задачах обычно не известно. Поэтому наиболее часто для оценки погрешности используют следующий прием: выполняют несколько расчетов при различных разбиениях области КЭ, по результатам этих расчетов строится зависимость рассчитанных напряжений (перемещений, деформаций) от размера элемента, затем выполняется экстраполяция на случай размера элемента, стремящегося к нулю. Однако эти ошибки аппроксимации можно свести к минимуму,   если при построении аппроксимирующих функций обеспечить:   1) непрерывность искомой функции и ее производных при переходе через границу КЭ до степени  m–1 включительно (m – наибольший порядок производных искомой функции содержащихся в функционале); 2) выполнение условий полноты, т. е. при уменьшении размеров КЭ аппроксимирующие функции должны обеспечить стремление значений искомой функции, а также ее производных к постоянным значениям; 3) выполнение условий совместности искомой функции и частично ее производных на границе между смежными элементами; 4) приближенное удовлетворение условий совместности не основных переменных (например, напряжений, если основные неизвестные – перемещения) на границах КЭ, а также граничных условий в рассматриваемой области;   5) исключение концентрации напряжений в КЭ, если в рассматриваемой области такие концентрации заведомо отсутствуют;  6) при перемещениях КЭ как жесткого целого в нем не должны возникать     деформации.  Требование полноты аппроксимирующих функций необходимо для учета смещения КЭ как жесткого целого и обеспечения состояния постоянных деформаций в элементе. Механический смысл совместности заключается в непрерывности основных неизвестных на смежных границах соседних КЭ. В сложных эрмитовых элементах выполнение условий совместности достигается сложнее.  Между тем имеются случаи, когда несовместные элементы дают очень хорошие результаты при быстрой сходимости решения к точному. Отметим еще одну важную с точки зрения практики расчетов особенность метода. МКЭ (в рассмотренной постановке) подбирает поле перемещений так, чтобы минимизировать некоторый функционал, имеющий энергетический смысл. Поэтому точность определения упругой энергии, запасенной в конструкции при заданных нагрузках, оказывается выше, чем точность определения перемещений. Точность определения напряжений оказывается ниже, чем точность определения перемещений, поскольку напряжения определяются по деформациям, получаемым дифференцированием перемещений, и ошибки численного дифференцирования могут играть заметную роль. С учетом ошибок округления ситуация оказывается более сложной: при большом числе элементов N решение может расходится из-за накапливающихся ошибок округления, даже если условия сходимости выполняются. Эти ошибки наиболее существенны, если конечные элементы сильно вытянуты или имеют углы, величина которых близки к 0º или 180º. В этом случае расчет напряженно-деформированного состояния элемента становится плохо обусловленным(часть вблизи очень острого угла "не чувствует", что происходит в остальном элементе). С целью не допустить здесь больших ошибок разработчики пакетов КЭ обычно ограничивают отношение сторон элемента и величины углов; в пакеты вводятся специальные средства проверки элементов, рекомендующие пользователю - если необходимо - перестроить сетку или делающие это автоматически. Наилучшим в этом смысле являются КЭ в виде правильных многоугольников (квадрат, равносторонний треугольник, куб, правильный тетраэдр); приемлемыми являются элементы с отношением сторон до - примерно 1:4 и углами от 25º до 155º. В настоящее время, при отсутствии общей теории оценки погрешности МКЭ, в реальных задачах практически единственным методом оценки точности является построение зависимостей типа показанной на рис 14 (зависимость результатов - напряжений, перемещений - от разбиения) и их экстраполяция с учетом опыта. 7. Преимущества и недостатки В настоящее время область применения метода конечных элементов очень обширна и охватывает все физические задачи, которые могут быть описаны дифференциальными уравнениями. Наиболее важными преимуществами метода конечных элементов, благодаря которым он широко используется, являются следующие: Свойства материалов смежных элементов не должны быть обязательно одинаковыми. Это позволяет применять метод к телам, составленных из нескольких материалов. Криволинейная область аппроксимирована с помощью прямолинейных элементов или описана точно с помощью криволинейных элементов. Таким образом, метод можно использовать не только для областей с "хорошей" формой границы. Размеры элементов могут быть переменными. Это позволяет укрупнить или измельчить сеть разбиения области на элементы, если в этом есть необходимость. Указанные выше преимущества метода могут быть использованы при составлении достаточно общей программы для решения частных задач определенного класса. Главный недостаток метода конечных элементов заключается в необходимости составления вычислительных программ и применения вычислительной техники. Вычисления, которые требуется проводить при использовании метода конечных элементов, слишком громоздки для ручного счета даже в случае решения очень простых задач. МКЭ является очень ресурсоемким методом с точки зрения затрать вычислительных ресурсов. 8. Пример применения FEA систем 8.1. Строительство Изначально МКЭ разрабатывался как метод строительной механики, но получил широкое применение и в других областях. Система Autodesk Robot Structural Analysis является специализированной системой анализа строительных конструкций с применением МКЭ на этапах монтажа и эксплуатации. Многие современные проекты проходят подготовительные этапы анализа в данной системе, к примеру, производился анализ деловой центр "Москва - Сити". 8.2. Краш тест В последнее время математические методы прогнозирования поведения автомобиля в различных видах ДТП занимают всё больше времени в объёме, отведённом на разработку автомобиля. Мировые автогиганты имеют в своём арсенале мощнейшие компьютерные средства для дополнения натурных экспериментов виртуальными испытаниями. Использование такого рода программ позволяет снизить затраты на проведение огромного количества испытаний как отдельных узлов автомобиля, так и краш-тестов автомобиля в целом. Стоимость одного краш-теста автомобиля составляет от 150 до 200 тысяч долларов, экспериментальной модели автомобиля до 2 млн. долларов, в то же время виртуальный краш-тест стоит 5-7 тысяч долларов При разработке новой модели автомобиля производители проводят 150-200 виртуальных краш-тестов и 5-6 реальных. При проведение виртуальных краш-тестов производители используют в основном расчётные системы такие, как ANSYS, LS-DYNA, Abaqus, РАМ-Crash (ESI Group), основанные на применении методов конечных элементов (МКЭ). Эти программы представляют собой удобную среду для проведения численных расчётов, дополняя эксперименты, проводимые в краш-лабораториях. В системе LS-DYNA была создана модель прямого удара автомобиля Plymouth Neon 1996 года выпуска на скорости 51 км/ч в недеформируемый барьер. 8.3. Пример исследования в монолитном строительстве Теперь рассмотрим поведение крупнощитовой опалубки серии ГАММА (3м x 1.2м) отечественной фирмы Техноком-БМ при различных видах нагрузки, смоделированной в системе SolidWorks Simulation. Модель представляют собой сборку из двух деталей: Металлическая конструкция, состоящая из рамы и ребер жесткости (в дальнейшем – Каркас щита) Палуба из ламинированной фанеры С помощью модуля Simulation был произведен статический анализ созданной модели. Задание граничных условий: -В соответствии с номенклатурой, в приведенной таблице указаны материалы и номера их ГОСТов, используемые для изготовления опалубки данной серии. Палуба - Фанера общего назначения; Профиль - Сталь 10. В следующей таблице приведены параметры(необходимые для задания материальных уравнений связи напряжения и деформации), запрашиваемые системой для проведения данного вида анализа. -При исследовании потребовалось изучение физических условий эксплуатации опалубочных элементов в монолитном строительстве. В зависимости от различных условий нагрузка может видоизменяться. Нагрузка была задана трех видов: равномерное нагрузка. В качестве значения используется предельно допустимое давление согласно номенклатуре 90кН/м2; треугольная нагрузка (левый нижний угол слайда); Усеченная треугольная нагрузка (правый нижний угол слайда). по следующему закону (P=ρh). Для бетонной смеси используемой в отечественном монолитном строительстве ρ=2500 кг/м2. -Кроме того, важно знать как в реальных условиях устанавливается и закрепляется модель или определить граничные условия с точки зрения МКЭ. Для имитации реальных видов креплений в модели было установлено 2 типа креплений в 4 местах установки тяжей: Имитирующий стяжку и позволяющий поворот только вокруг своей оси. Имитирующий опору гайки и позволяющий перемещение только по плоскости прилегания. В системе SolidWorks Simulation сетка создается автоматическим построителем, поэтому пользователь может управлять только стартовыми параметрами для создания сетки и не может работать с каждым элементом, поэтому в данной системе есть высокая вероятность создания не качественных элементов (слишком вытянутых). В данной системе этот недостаток постарались компенсировать использованием элементов высокого порядка. В SolidWorks Simulation для оболочек используются элементы - плоские треугольники или параболический треугольный элемент (элемент второго порядка), а для твердых тел - тетраэдры чернового качества или параболические тетраэдры высокого качества. Использования такого рода элементов позволяет повысить точность, но усложняет дальнейший расчет. Средствами системы Simulation, для запуска анализа МКЭ, была создана трехмерная расчетная сетка, представленная на слайде и ее укрупненный вариант. Сетка в данном примере комбинированная. Так как каркас достаточно тонкий и не позволяет создать сетку достаточной точности в несколько рядов элементов, то каркас задается как оболочка которая дискретезируется плоскими треугольными элементами второго порядка, а палуба, как твердое тело, тетраэдрами высокого качества. На следующем слайде представлены эпюры перемещения опалубки при воздействия на нее разного вида нагрузок. Как видно в зависимости от формы давления видоизменялась и форма деформации опалубки. В результате выполнения анализа в системе Simulation было выявлено, что, вероятней всего, ни один из данных видов нагрузок опалубка с данными характеристиками материалов и конструкцией не выдерживает в местах установки тяжей. Заключение Метод конечных элементов является одним из наиболее распространенных методов проведения инженерных анализов разрабатываемых деталей. С помощью данного метода можно исследовать детали любой сложности и неоднородными материалами. В большинстве современных САПР есть возможность проведение инженерного анализа, с применением МКЭ. Популярность МКЭ обусловлена его результативностью и возможностью решения широкого спектра задач.

Похожие:

	Моделирование течения ледников методом конечных элементов ...		Рабочая программа составлена в соответствии с требованиями фгос во... Салтанова Т. В. Метод конечных элементов в расчётах прочности. Учебно-методический комплекс. Рабочая программа для студентов направления...
	Повышение точности технологических систем горизонтальных координатно-расточных... ...		Решение заданий из №128, 129 у доски по схеме ...
	Министерство образования Российской Федерации Санкт Петербургский... Задачи курса: Изучить основные математические результаты и методы, лежащие в основе метода конечных элементов и других вариационных...		Прогнозирование трещиностойкости бетона на основе метода конечных элементов Реальное строение материала и особенности его поведения под нагрузкой отражено в структурных теориях прочности. Однако практическое...
	Методические указания к выполнению расчетно-графического задания... В методических указаниях изложены краткие сведения о пк «лира», о методе конечных элементов, реализуемом в системе «лира», рассмотрен...		Образовательная программа «Компьютерные и информационные науки» Форма... Ук-1 Способность к критическому анализу и оценке современных научных достижений, генерированию новых идей при решении исследовательских...
	Физический факультет Сравнительное исследование парамагнитных свойств образцов, полученных методом пиролиза аэрозолей и золь-гель методом 19		Окислительно-восстановительные реакции Уметь определять степени окисления элементов в простых и сложных веществах, различать понятия: степень окисления, составлять уравнения...
	При некоторых хронических бактериальных заболеваниях с успехом применяется Сравнительный и количественный анализ белков мочи и перитонеальной жидкости, разделенных методом гель		Методические указания к лабораторным работам по дисциплине «Управление... Методические указания к лабораторным работам по дисциплине «Управление проектами» для студентов и слушателей факультета «Инженерный...
	Ее Величества Реакции Цель: продолжить формировать умения учащихся в определении степени окисления элементов, процессов окисления и восстановления в овр;...		Хорошавин Лев Борисович Докт техн наук реферат Периодической системе элементов, но и дополнительно по кластерам химических элементов, определять прогнозные свойства новых элементов...
	Курсовая работа по курсу «Теория случайных потоков» на тему «Анализ... «Анализ надёжности электроснабжения подстанции «Новая» методом случайных потоков»		Методические рекомендации к самоанализу урока с позиций системно деятельностного подхода Основой любого анализа урока может служить структурный анализ. Структурный (поэтапный) анализ – это выявление и оценка доминирующих...