книги / Одноковшовые погрузчики
..pdfшинство элементов основаны на плоском эскизе, по которому создается базовый трехмерный объект. Последовательное наращивание 3D-объектов и позволяет в итоге получить желаемый результат (рис. 14).
Основы метода конечных элементов
Метод конечных элементов (МКЭ) в настоящее время является стандартом при решении задач механики твердого тела посредством численных алгоритмов. МКЭ занял лидирующее положение благодаря возможности моделировать широкий круг объектов и явлений. Абсолютное большинство конструктивных элементов, узлов и конструкций, изготовленных из самых разнообразных материалов, имеющих различную природу, могут быть рассчитаны посредством МКЭ. При этом нужно учитывать неизбежные при любой численной аппроксимации условности и погрешности. Поэтому вопрос соответствия между расчетной моделью и реальностью является основным при использовании программ анализа. Несмотря на то, что такие программы имеют более или менее подробную документацию, они все равно остаются в определенной степени черными ящиками. Это означает определенную непредсказуемость результатов, а также некоторый произвол в их интерпретации. Следовательно, качество заключений, принимаемых на основе результатов, всецело зависит от квалификации, а также, применительно к расчету на прочность, принципиального знакомства с основами МКЭ.
В основе метода лежит дискретизация (разбивка на более мелкие, обособленные, но взаимосвязанные между собой элементы) объекта с целью решения уравнений механики сплошной среды в предположении, что эти соотношения выполняются в пределах каждой из элементарных областей. Эти области называются конечными элементами. Они могут соответствовать реальной части пространства, как, например, пространст-
111
венные элементы (рис. 15, а, б), или же быть математической абстракцией, как элементы стержней, балок, пластин или оболочек (рис. 15, в). В пределах конечного элемента назначаются свойства ограничиваемого им участка объекта (это могут быть, например, характеристики жесткости и прочности материала, плотность и т.д.) и описываются поля интересующих величин (применительно к механике твердого тела это перемещения, деформации, напряжения и т.д.). Параметры из второй группы
аб
в
Рис. 15. Объемный конечный элемент
112
назначаются в узлах элемента, а затем вводятся интерполирующие функции (разность между «соседними» значениями), посредством которых соответствующие значения можно вычислить в любой точке внутри элемента или на его границе. Задача математического описания элемента сводится к тому, чтобы связать действующие в узлах факторы. В механике сплошной среды это, как правило, перемещения и усилия.
Построение программы МКЭ
Имея математический аппарат для получения матриц жесткости конечных элементов, приведения нагрузок, приложенных к поверхности или в объеме элемента к усилиям в узлах, а также решения обратных задач: вычисления полей деформаций и напряжений в объеме элемента на базе перемещений в узлах – можно построить алгоритм МКЭ. Приведем один из его вариантов для решения задач в линейной постановке [9].
1.Производится дискретизация объема, занимаемого деталью или сборкой на элементы, или, иначе, строится сетка конечных элементов. Для объемного тела область делят на тетраэдры с гранями, аппроксимируемыми линейными (линейная зависимость от координат) или параболическими функциями координат. Для поверхностных моделей – на плоские (линейная) или криволинейные (параболическая зависимость) треугольники.
2.Для пространственных конечных элементов степенями свободы являются перемещения в направлении осей локальной системы координат элемента. Для конечных элементов оболочек к трем перемещениям в каждом узле добавляются по три угла поворота нормали к срединной поверхности области, аппроксимируемой элементом, относительно тех же осей.
3.Определяются зависимости для преобразования перемещений и углов поворота в узлах к глобальной системе координат.
113
4.Вычисляются матрицы жесткости конечных элементов.
Вформулы для расчета компонентов матриц жесткости конечных элементов помимо координат узлов входят модули упругости и коэффициенты Пуассона материалов. Другими словами, если анализируется сборка, то в зависимости от принадлежности элемента детали при расчете матриц жесткости элементов используются соответствующие характеристики жесткости материала.
5.Полученные матрицы жесткости с использованием зависимостей для перехода от локальных систем координат элемента в глобальные преобразуются в глобальную систему координат.
6.Матрицы жесткости, представленные в глобальных координатах, объединяются в глобальную матрицу жесткости
[K].
7.Назначенные пользователем граничные условия, статические и кинематические, приводятся к нагрузкам и перемещениям в узлах, выраженным в глобальной системе координат, и включаются в столбец усилий [F].
8.Полученная линейная система уравнений вида [K]·[∆] = = [F] решается относительно столбца перемещений. Это наиболее трудоемкий этап расчета. Для решения используются итерационные или прямые методы. Матрица жесткости, как правило, хранится в компактной форме, структура которой определяется до этапа ее заполнения матрицами жесткости элементов.
9.Для каждого конечного элемента, с учетом перемещений (углы поворота) в узлах и аппроксимирующих функций, рассчитываются деформации. Если элементы линейные – деформации в пределах элементов постоянные, если элементы параболические – деформации изменяются линейно. На основе деформаций вычисляются напряжения в элементах. При необходимости (функция программы) напряжения в узлах смежных элементов осредняются (это также весьма ответственный этап,
114
по-разному решаемый в различных программах) с последующим пересчетом напряжений в пределах каждого элемента.
10. На основе компонентов напряженно-деформированно- го состояния и параметров прочности материала (материалов) производится вычисление эквивалентных напряжений по ка- кому-либо критерию прочности.
В среде SolidWorks эти расчеты ведутся при помощи встроенной программы COSMOSWorks.
COSMOSWorks – приложение к SolidWorks, предназна-
ченное для решения задач механики деформируемого твердого тела методом конечных элементов. Программа использует геометрическую модель детали или сборки SolidWorks для формирования расчетной модели. Интеграция с SolidWorks дает возможность минимизировать операции, связанные со специфическими особенностями конечно-элементной аппроксимации. Назначение граничных условий производится в привязке к геометрической модели. Такими же особенностями обладают и процедуры представления результатов. Развитие программы происходит как в направлении увеличения вычислительных возможностей, так и, в значительной степени, за счет расширения функциональности, связанной с имитацией расчетных ситуаций, возникающих в повседневной инженерной практике.
COSMOSWorks позволяет выполнять следующие виды моделирования:
–статический анализ в упругой постановке с расчетом отдельных деталей по пространственной или оболочечной модели, а также сборок в трехмерной постановке с учетом взаимодействия деталей;
–расчет собственных частот и соответствующих им форм для деталей в твердотельном или оболочечном представлении,
атакже сборок с неподвижными деталями;
–расчет величин критических нагрузок потери устойчивости и соответствующих им форм для деталей в твердотель-
115
ном или оболочечном представлении, а также сборок с неподвижными деталями;
–тепловой расчет с учетом явлений теплопроводности, конвекции, излучения, но без учета движения сред;
–термоупругий анализ на базе результатов теплового расчета;
–параметрическая оптимизация по критерию минимизации/максимизации массы, объема, собственных частот и критической силы;
–имитация деформирования конструкции с учетом физической и геометрической нелинейности, а также ввиду изменения нагрузок и температуры во времени;
–моделирование эффекта падения конструкции на жесткую поверхность;
–усталостный расчет с учетом кривых усталости, формы кривой нагрузки, а также линейной гипотезы суммирования повреждений.
Все эти типы анализа могут быть связаны с одним и тем же объектом SolidWorks.
Последовательность расчета
COSMOSWorks требует соблюдения базовой канвы алгоритма метода конечных элементов, предоставляя внутри каждого этапа определенную свободу в последовательности шагов подготовки модели и рассмотрения результатов.
Предполагаемая цепочка событий для расчета в упругой постановке:
1.Создание анализа определенного типа и определение его настроек. Последние могут быть изменены в любой момент перед выполнением расчета.
2.Заполнение, если необходимо, таблицы параметров, определяющей набор величин, которые могут изменяться (кон-
116
кретно – для которых могут назначаться списки значений)
входе расчета.
3.Подготовка исходных данных внутри заданного анализа:
• назначение материала детали или деталям;
• назначение кинематических граничных условий;
• назначение статических граничных условий;
• создание сетки.
4.Связывание, в случае необходимости, параметров из таблицы параметров с соответствующими анализами.
5.Выполнение расчета.
6.Обработка результатов:
•создание необходимых диаграмм;
•анализ диаграмм;
•экспорт результатов.
Расчет ковша фронтального погрузчика (рис. 16–21)
Рис. 16. Ковш
117
Информация о модели:
Имя модели: ковш Активная конфигурация: По умолчанию
Твердые тела
|
<L_MdInf_SldBd_Nm/> |
|
Рассматри- |
|
Объемные |
|
Путь документа/ |
|
|
|
вается как |
|
свойства |
|
Дата изменения |
|
|
|
|
|
||||||
|
|
|
Твердое |
|
Масса: 984,907 kg |
|
C:\Users\Alex\ |
|
|
|
|
|
|
|
|||
|
|
|
тело |
|
Объем: 0,12791 m^3 |
|
Desktop\Стрела |
|
|
|
|
|
|
Плотность: |
|
погрузчика |
|
|
|
|
|
|
7700 kg/m^3 |
|
3D\ковш.SLDPRT |
|
|
|
|
|
|
Масса: 9652,09 N |
|
Dec 20 2013 |
|
|
|
|
|
|
|
|
|
|
Свойства исследования:
Имя исследования |
Статический анализ 1 |
Тип анализа |
Статический анализ |
Тип сетки |
Сетка на твердом теле |
Тепловой эффект: |
Вкл |
Термический параметр |
Включить тепловые нагрузки |
Температура при нулевом напряжении |
298 Kelvin |
Включают эффекты давления жидкости из |
Выкл |
SolidWorks Flow Simulation |
|
Тип решающей программы |
FFEPlus |
Влияние нагрузок на собственные частоты: |
Выкл |
Мягкая пружина: |
Выкл |
Инерционная разгрузка: |
Выкл |
118 |
|
Несовместимые параметры связи |
|
|
Автоматические |
|
||
Большие перемещения |
|
|
Выкл |
|
||
Вычислить силы свободных тел |
|
|
Вкл |
|
||
Трение |
|
|
|
Выкл |
|
|
Использовать адаптивный метод: |
|
|
Выкл |
|
||
Папка результатов |
|
|
|
Документ SolidWorks |
||
Единицы измерения: |
|
|
|
|
||
|
|
|
|
|
||
Система единиц измерения: |
|
СИ (MKS) |
|
|||
Длина/Перемещение |
|
|
mm |
|
||
Температура |
|
|
Kelvin |
|
||
Угловая скорость |
|
|
Рад/сек |
|
||
Давление/Напряжение |
|
N/m^2 |
|
|||
Свойства материала: |
|
|
|
|
||
|
|
|
|
|||
Ссылка на модель |
|
Свойства |
Компоненты |
|||
|
|
Имя: |
Легированная сталь |
Твердое те- |
||
|
|
Тип модели: |
Линейный |
ло (ковш) |
||
|
|
|
Упругий |
|
||
|
|
|
Изотропный |
|
||
|
|
Критерий прочно- |
Неизвестно |
|
||
|
|
сти по умолчанию: |
|
|
|
|
|
|
Предел текучести: |
6.20422e+008 N/m^2 |
|
||
|
|
Предел прочности |
7.23826e+008 N/m^2 |
|
||
|
|
при растяжении: |
|
|
|
|
|
|
Модуль упругости: |
2.1e+011 N/m^2 |
|
||
|
|
Коэффициент |
0.28 |
|
||
|
|
Пуассона: |
|
|
|
|
|
|
Массовая плот- |
7700 kg/m^3 |
|
||
|
|
ность: |
|
|
|
|
|
|
Модуль сдвига: |
7.9e+010 N/m^2 |
|
||
|
|
Коэффициент теп- |
1.3e-005 /Kelvin |
|
||
|
|
лового расшире- |
|
|
|
|
|
|
ния: |
|
|
|
|
|
|
|
|
|
|
|
Данные кривой:N/A |
|
|
|
|
|
|
|
|
|
|
|
|
119 |
Нагрузки и крепления:
|
|
Имя крепления |
Изображение |
|
|
|
|
Данные крепления |
|||||||||
|
|
крепления |
|
|
|
|
|||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
Объекты: |
1 грани |
|||||
Зафиксированный-1 |
|
|
|
|
|
|
|
Тип: |
Зафиксированная |
||||||||
|
|
|
|
|
|
|
|
|
|
|
|
геометрия |
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Результирующие силы |
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
Компоненты |
|
X |
|
|
Y |
|
|
|
Z |
|
|
|
Результирующая |
|
|
|
|
Сила реакции(N) |
|
229353 |
|
|
-200423 |
|
|
-2643.39 |
|
|
|
304597 |
|
||
|
|
Реактивный |
|
0 |
|
|
0 |
|
|
0 |
|
|
|
|
0 |
|
|
|
|
момент(N.m) |
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Объекты: 4 грани |
||||||
|
|
|
|
|
|
|
|
|
|
|
|
Тип: Зафиксированный |
|||||
|
Зафиксированный |
|
|
|
|
|
|
|
|
|
|
|
шарнир |
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
шарнир-1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
Результирующие силы |
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
Компоненты |
|
X |
|
|
Y |
|
|
|
Z |
|
|
|
Результирующая |
|
|
|
|
Сила реак- |
|
-304349 |
|
|
330343 |
|
|
2640.61 |
|
|
|
449179 |
|
||
|
|
ции(N) |
|
|
|
|
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Реактивный |
|
0 |
|
|
0 |
|
|
0 |
|
|
|
|
0 |
|
|
|
|
момент(N.m) |
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Имя нагрузки |
Загрузить |
|
|
|
|
Загрузить данные |
|||||||||
|
|
изображение |
|
|
|
|
|||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
Объекты: |
1 грани |
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Тип: |
Приложить |
||
|
|
Сила-1(Подается |
|
|
|
|
|
|
|
|
|
|
|
|
нормальную |
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
силу |
|||
на режущую кромку |
|
|
|
|
|
|
|
|
Значение: |
||||||||
|
|
|
|
|
|
|
|
150000 N |
|||||||||
|
|
ковша) |
|
|
|
|
|
|
|
|
|||||||
|
|
|
|
|
|
|
|
|
Угол фазы: |
0 |
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
Единицы deg |
|||||
|
|
|
|
|
|
|
|
|
|
|
измерения: |
|
|
120