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

Метод конечных элементов (МКЭ) - это математический метод решения дифференциальных уравнений, позволяющий с высокой точностью моделировать реальность, в том числе, в таких задачах как прочность и деформация конструкций и их элементов.

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

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

 Последовательность процедур алгоритма МКЭ может быть представлена в следующем виде:

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

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

2. Выбор вариационного принципа.

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

3. Выбор аппроксимирующих функций.

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

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

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

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

4. Реализация вариационного принципа.

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

 - методом непосредственного сложения жесткостей;

 - методом конгруэнтного преобразования;

 - при помощи конечно-разностных операторов.

5. Учет граничных условий.

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

6. Решение системы алгебраических уравнений.

7. Определение деформаций и напряжений.

После определения узловых перемещений в соответствии с известными соотношениями теории упругости могут быть определены деформации и напряжения.

На рисунке 1  приведен алгоритм проведения прочностных расчетов элементов пластины  полотна конвейера с использованием широко известного пакета моделирования и конечно-элементного анализа MSC.visualNastran.

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

Основу MSC.visualNastran  составляет отработанная технология элементов и надежные численные методы.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 


Рисунок 1 - Алгоритм проведения прочностных расчетов элементов

пластины  полотна конвейера

 

 

На рисунке 2 представлена рабочая среда моделирования MSC.visualNastran 4D.

 

 

Рисунок 2 - Рабочая среда моделирования MSC.visualNastran 4D

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

На рисунках 3 – 7  представлено распределение напряжений, деформаций  и перемещений пластины при падении груза массой 160кг.

Распределение напряжений 160 кг

 

Рисунок 3  - Распределение напряжений

 

 

Распределение напряжений с сеткой 160 кг

 

Рисунок 4 - Распределение напряжений с сеткой

 

Распределение деформаций с деформированной формой 160 кг

 

Рисунок 5 - Распределение деформаций  с деформированной формой 

 

Распределение перемещений с деформированной формой 160 кг

 

 

Рисунок 6 - Распределение перемещений с деформированной формой

Распределение запаса прочности 160 кг

 

Рисунок 7  - Распределение запаса прочности 160кг