Повышение быстродействия расчета напряженно-деформированного состояния

Размещено на http://www.allbest.ru/

Повышение быстродействия расчета напряженно-деформированного состояния методом дискретных элементов

Н.О. Шумейко

А.Е. Андреев

Компьютерное моделирование является в настоящее время одним из основных способов исследования динамики сложных механических систем. Оно является неотъемлемым этапом проектирования, оптимизации параметров конструкций и механизмов и широко применяется в различных областях науки и техники, таких как железнодорожный и автомобильный транспорт, авиастроение, робототехника. В связи с этим актуально создание и развитие программных комплексов, специализирующихся на компьютерном моделировании [1,2]. Одним из таких комплексов является программная система формирования решений уравнений нелинейной динамики (ФРУНД), предназначенная для моделирования динамических процессов в машинах и конструкциях. ФРУНД позволяет решать широкий круг задач анализа и синтеза машиностроительных конструкций, одна из которых -- моделирование системы тел.

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

Для оценки характеристик движения всей механической системы в целом используется метод моделирования динамики системы абсолютно твёрдых тел со связями [3]. Метод моделирования систем тел со связями применим к механической системе, состоящей из абсолютно твёрдых тел и связей между ними. Multi-bodysystem (MBS) модель рассматривается как основная, а модели, описывающие отдельные тела как вспомогательные [4]. Можно представить структуру MBS модели в виде графа, вершины которого соответствуют телам, а ребра соответствуют связям между телами, то для различных моделей могут получаться графы с циклами и графы в виде деревьев. программный нелинейный дискретный вычисление

Уравнения движения тела, освобожденного от связей с локальной системой координат в центре масс:

, (1)

где - масса i-го элемента, - ускорение, - вектор результирующей силы (сумма активных и реактивных сил), - матрица проекций угловых скоростей в системе координат i-го элемента, - вектор результирующего момента, - тензор инерции, - вектор углового ускорения, - вектор угловой скорости[5]. Верхний индекс обозначает систему координат, в которой задана величина. Нижний индекс показывает порядковый номер элемента модели.

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

,

где - вектор множителей Лагранжа, - матрица инерции, x- вектор абсолютных координат, - вектор правых частей, - вектор правых частей уравнения связей [6].

В матричном виде формула запишется как

.

Матрица системы сильно разрежённая [7].

Обзор программного комплекса

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

В архитектуре программного комплекса помимо ядра можно выделить две большие подсистемы: графическую оболочку FShell для работы с моделью и модули пакета ФРУНД.

Основные функции графической оболочки:

-- редактирование каркасной геометрии тел

-- включение в отдельные тела вспомогательных решателей

-- добавление CAD геометрии к телам

-- задание граничных условий в виде поверхностей CAD геометрии и

генерация регулярной ортогональной сетки

-- редактирование параметров модели

-- анимация результатов моделирования мультифизической модели

Ядро (компонент FCore) выполнено в виде консольного приложения, которое выполняет различные операции с моделью ключами командной строки. Интерфейс командной строки - это единственный интерфейс, который реализует ядро, он предоставляется графической оболочке или пользователю [8].

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

В результате решения СЛАУ получаются значения вторых производных от координат и значения множителей Лагранжа. После решения СЛАУ производится численное интегрирования системы обыкновенных дифференциальных уравнений. Для численного интегрирования применяются явные и неявные методы. Среди явных для рассматриваемой задачи наиболее распространённый метод Рунге-Кутты 4 порядка, среди неявных - метод Адамса-Башфорта-Мултона по схеме прогноза-коррекции. В данной работе используется метод Рунге-Кутты [9, 10].

После завершения работы программы данные расчетов записываются в файл, который можно открыть с помощью FShell, воспользоваться визуализацией модели и пронаблюдать процесс деформации по итерациям. Показана модель размером 10х10х10 (куб). В правой части формы визуализации отображена шкала напряженности, от 0 до определенного максимального значения. Цвета точек соответствуют значению напряженности из спектра (шкалы).

Оценка производительности отдельных частей метода

Для оценки производительности отдельных частей метода были проведены тесты на моделях размерностями 10х10х10, 20х20х20, 30х30х30, 40х40х40, 50х50х50, 100х100х100. По одному прогону на итерацию, в тесте одна итерация. Для тестирования использовался ПК с процессором AMD PhenomII x6.

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

Таблица 1. Результаты тестов производительности

В среднем вычисление коэффициентов занимает 4.7%, а вычисление вторых производных 63.5% общего времени.

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

Для распараллеливания была выбрана часть метода, отвечающая за пересчет вторых производных, использовалась библиотека OpenMP. Тестирование проводилось на ПК с процессором AMD PhenomII x6.

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

Таблица 2. Ускорение одной итерации

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

При расчетах малых моделей наблюдается меньшее ускорение. Это происходит из-за накладных расходов на обеспечение параллельного выполнения.

Как правило, расчеты модели не ограничиваются одной итерацией.

Таблица 3. Ускорение метода с различным количеством итераций

В таблице - общее время выполнения, - общее время выполнения с параллельной вставкой.

В среднем ускорение () составляет 2.43 раза. При этом вне зависимости от количества итераций значение ускорения не сильно отклоняется от среднего.

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

Литература

1. Кадомцев М.И., Ляпин А.А., Тимофеев С.И. К вопросам построения эффективных алгоритмов расчета системы «сооружение-грунт» // «Инженерный вестник Дона», 2012, №1. -ivdon.ru/magazine/archive/n1y2009/250.

2. М. И. Кадомцев, А. А. Ляпин, Шатилов Ю.Ю. Вибродиагностика строительных конструкций // «Инженерный вестник Дона», 2012, №3. -ivdon.ru/magazine/ archive/n3y2012/941

3. Виттенбург Й. Динамика систем твердых тел / Под ред. В. В. Румянцев.-- М.: Мир, 1980.-- С. 294.

4. Z. Qi, Y. Xu, X. Luo, S. Yao Recursive formulations for multibody systems with frictional joints based on the interaction between bodies // Multibody System Dynamics.-- 2010.-- Vol. 24.-- pp. 133-166.

5. Fleissner F., Gaugele T., Eberhard P. Applications of the discrete element method in mechanical engineering // Multibody System Dynamics.--2007.-- Vol. 18.-- pp. 81-94.

6. Б.Л. Глушак, В.Ф. Куропатенко, С.А. Новиков, В.М. Фомин Исследование прочности материалов при динамических нагрузках.-- Новосибирск: Наука. Сиб. отд-ние, 1992.-- С. 295

7. Писсанецки С. Технология разреженных матриц / Под ред. Х. Д. Икрамов.-- М. : Мир, 1988.-- С. 411.

8. Горобцов А. С., Карцов С. К., Кувшид Р. П. Комплекс ФРУНД - инструмент исследования динамики автомобиля // Автомобильная промышленность, 2005. № 4. С. 27-28.

9. Гетманский В. В., Горобцов А. С. Расчет динамики систем тел и физических процессов в отдельных телах на кластере // Математические методы в технике и технологиях, ММТТ-25: сб. тр. XXVмеждунар.науч. конф.-- Т. 3.-- ВолгГТУ: Саратов, 2012.-- С. 171-174.

10. Гетманский В. В., Сергеев Е. С., Горобцов А. С. Способы ускорения процедуры расчета при междисциплинарном моделировании // Матер.междунар. конф. СКТ-2010. -- Т. 1. -- Таганрог : Изд-во ТТИ ЮФУ, 2010.-- С. 199-222.

Размещено на Allbest.ru