Пример: распечатка множителя L (или его части) неупорядоченной и упорядоченных

алгоритмами СМ и RCM матрицы жесткости NRC=3 для Example 2.

Множитель L. (пример таблицы)

Далее анализ и выводы, при необходимости подкреплённые графиками.

П.6. Недостатки САЕ Sigma и подсистемы Preprocessor, предложения по совершенствованию комплекса Sigma и его подсистем.

Студент отмечает недостатки и достоинства препроцессора и системы Sigma и предлагает пути их совершенствования. Обязательными также являются указания на ошибки в методических указаниях к КР, предложениями иных формулировок и т.д.

Отчет по П.6 должен содержать:

Достоинства и недостатки препроцессора и предложения по их устранению.

Достоинства и недостатки системы Sigma и предложения по их устранению.

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

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

За правильное выполнение 1-ой обязательной части КР проставляется оценка "удовлетворительно" и студент вправе завершить её на этом этапе.

2. Вторая часть курсовой работы в 6-го семестра.

За выполнение одного из пунктов второй части КР окончательная оценка за КР 6-го семестра увеличивается на один балл, двух пунктов – два балла.

П.7. Качество сетки КЭ и влияние значений входных параметров на результаты решения задачи МКЭ.

В П.7 студент:

• определяет качество оптимизированной и неоптимизированной сеток проекта 2 КР6;

• оценивает влияние оптимизации сетки на результаты решения задачи;

• отбирает наилучшую сетку КЭ.

Введение

Оценка качества сетки конечных элементов

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

Поэтому приходится искать способы оценки качества сетки КЭ по косвенным признакам, при этом не только после её построения, но и, желательно, в процессе триангуляции.

Очевидным критерием является минимальная площадь КЭ. Понятно, чем меньше площадь конечных элементов, тем больше модель приближается к реальному объекту и тем более точно рассчитывается задача, если, конечно, накапливаемая при машинном округлении ошибка не становится критической. Оценка минимальной площади конечных элементов может оцениваться не по всему рассчитываемому объекту, а только по наиболее ответственной его части, например в области наибольших значений напряжений или в области, которая более всего интересует исследователя.

Дальше идёт интегральная вариация этого критерия – средняя площадь КЭ, рассчитываемая как по всему объекту, так и по его части.

Следующие критерии оценки качества сетки КЭ так или иначе анализируют форму конечных элементов.

Установлено, например, что ошибка метода КЭ при решении на треугольной сетке обратно пропорциональна величине синуса минимального угла в элементах сетки.

Отсюда в качестве естественного критерия ка

чества сетки можно принять синус минимального угла в КЭ сетки (или просто угол).

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

Оба эти простейших критерия слабо затрагивает форму большинства КЭ, изменение плотности КЭ в разных подобластях рассчитываемого объекта и поэтому могут использоваться только как вспомогательные.

Вообще, численных рекомендаций по форме КЭ немного.

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

Для формы КЭ, а с ним – интегральной оценки качества полученного сеточного разбиения строятся различные измерители (см. табл.).

Одним из таких измерителей является коэффициент формы, который вычисляется следующим образом. Для каждой стороны элемента Li определяется площадь идеального элемента такой величины (для равностороннего треугольника она равна 0,433(Li)^2, а для квадрата — (Li)^2), и затем эти площади осредняются.

Отношение этой осредненной «идеализированной» площади к реальной площади элемента принимается в качестве меры качества формы КЭ.

В литературе встречаются такие рекомендации:  

• углы КЭ должны быть не менее 30º и не более 150º;  

• не рекомендуется использовать треугольные элементы, в которых присутствует угол меньше 15º.

• отношение же сторон КЭ не более 1:10.

Другие измерители по оценке формы КЭ можно разработать на основе отношения площади КЭ к площади, вписанной в него окружности или на основе отношения радиуса вписанной в треугольник окружности к радиусу описанной.

Известно, например, что для любого треугольника , где R - радиус описанной окружности, а r - радиус вписанной окружности. Причем равенство имеет место для правильного треугольника.

Можно использовать некоторые известные соотношения в треугольнике. Например, отношение площади треугольника, вершинами которого являются точки касания вписанной окружности, к площади данного треугольника АВС меньше или равно . При этом, равно – для равностороннего треугольника.

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

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

Наиболее объективная интегральная оценка качества сеток осуществляется с использованием метода штрафных функций. В данном задании П.7 КР она не используется.

В П.7 для оценки качества сетки надо подобрать критерии, согласно простому большинству которых сетка будет обладать преимуществом или, наоборот, уступать сравниваемой сетке.

В качестве основных критериев в П.7 должны использоваться критерии минимального угла и среднего минимального угла, благо информация по эти критериям получена при выполнении П.3 К этим двум критериям необходимо привлечение ещё каких-либо критериев или критерия.

Содержание работы по П.7