Full Text

Введение Вариационные методы играют важную роль во многих разделах механики и применимы как при решении линейных, так и нелинейных задач [1-7]. Эти методы особенно эффективны при решении вариационных задач с различного рода ограничениями. Однако при решении вариационных задач со сложными областями определения оптимизируемых функционалов возникают известные затруднения, связанные с отысканием базисной системы функций. Частично данные трудности преодолеваются на основе вариационно-разностных методов [3, 8-11]. Заметим в этом контексте, что эффективное решение широкого класса краевых задач со сложными областями достигается на основе применения метода конечных элементов [12-15]. Однако метод конечных элементов, основанный на построении глобальной матрицы жесткости и решении систем линейных уравнений большой размерности, оказывается не столь эффективным в задачах с локальными ограничениями (в частности, для задач механики контактного взаимодействия, включая задачи с препятствиями, а также задач для сложных конструкций с предельными условиями на перемещения и напряжения), которые требуют разработки нетрадиционных подходов [16-23]. В связи с этим представляется целесообразным разработка методов, основанных на сочетании конечно-элементных аппроксимаций и алгоритмов локального варьирования, использующих преимущества метода конечных элементов и метода локальных вариаций. Заметим, что классический метод локальных вариаций [24-27] оказался во многих случаях очень эффективным как метод математического программирования. Этот метод требовал минимальной информации для выполнения улучшающих вариаций решения и оказывался особенно эффективным с ограничениями на функции состояния (фазовые переменные). В настоящей работе развивается альтернативный подход к решению вариационных задач с ограничениями типа неравенств, не требующий использования методов математического программирования. Данный подход основан на методе конечных элементов с аппроксимациями Ритца-Галеркина для искомых функций. В предлагаемом подходе выполняются локальные вариации искомых функций, и проводится сравнение вариаций минимизируемого функционала. 1. Алгоритм метода конечных элементов с локальными вариациями Предлагаемый алгоритм метода конечных элементов с локальными вариациями рассматривается в контексте следующей задачи. Требуется найти функцию , определенную в области в плоскости , удовлетворяющую ограничению , и минимизировать функционал где a является заданной функцией ее аргументов, множество задано в трехмерном пространстве , причем нижними индексами обозначены частные производные. В частности, в качестве условий (1) будут рассмотрены двусторонние ограничения Здесь , причем и являются заданными функциями в области . Для упрощения дальнейших рассмотрений декомпозируем исходную область на равные прямоугольные области (ячейки) прямыми линиями. Имеем где - число ячеек. Неизвестные величины функции в узлах сетки обозначены как , т.е. . Решение вариационной задачи представляется в виде где - число узлов искомых величин , а основные функции определяются в малых подобластях с ячейками в вершинах. Базисные функции принадлежат Рис. 1. Базисная функция Fig. 1. The basis function пространству Соболева () и удовлетворяют существенным граничным условиям. Семейство базисных функций удовлетворяет следующим требованиям (рис. 1): 1) 2) только для тех элементов, которые соединены с узлами с узлами 3) 4) являются сужением на , где Таким образом, искомое решение удовлетворяет следующему равенству: Заметим, что билинейный полином применяется в качестве ограничивающей функции для прямоугольных ячеек (элементов), коэффициенты которой должны удовлетворять третьему условию соотношения (7). Множество узлов и элементов, принадлежащих области , составляют конечно-элементную сеть. При этом мы получаем базисную систему функций, которые являются полиномиальными и имеют малую область определения. Ограничения (1) и (4) трансформируются в рамках конечно-элементных аппроксимаций как где . Представленный интеграл (3) приближенно заменяется суммой интегралов по ячейкам, полностью принадлежащим области (число ячеек обозначается через ). Имеем , . В итоге преобразований рассматриваемая задача (1), (2), (3) и (4) может быть сформулирована как задача отыскания , удовлетворяющих ограничениям (10) (или (11)) и доставляющих минимум функционалу (12), (13). Решение задачи минимизации функционалов (12), (13) будем строить последовательными приближениями. В качестве начального приближения для возьмем любой набор чисел, удовлетворяющих (10) (или (11)). При выборе начального приближения желательно учесть имеющуюся априорную информацию о предполагаемом решении, различные качественные, физические или иные соображения. Удачный выбор начального приближения может ускорить сходимость метода. Приведем предлагаемый метод последовательных приближений. Зададимся некоторым малым положительным числом , которое будем дальше называть шагом варьирования. На каждом шаге итерационного алгоритма операции последовательно выполняются для всего набора узлов, при этом для каждого узла операция выполняется только один раз. Для нахождения нового значения для каждого узла рассмотрим в качестве новых возможных значений старое значение , а также величины и . Отметим, что при вычислении функционала потребуется пересчет в (12), (13) компоненты максимум только для четырех ячеек, связанных с , т.е. это слагаемое связано с ячейками, лежащими полностью в области и содержащими точку () в качестве узла. Определим эту величину как . В результате будем производить замену старого значения на величины или , так, если для данных величин выполняются ограничения (10) (или (11)) и происходит уменьшение функционала , т.е. или , и выполняется неравенство то величины заменяются значением или величиной , т.е. . Затем алгоритм переходит к следующему узлу. В противном случае в качестве решения принимается старое значение . Число вычислительных операций может быть сокращено, если в качестве приоритетного направления изменения узлового значения будет приниматься направление, на котором происходит максимальное уменьшение функционала для предыдущего узла. После завершения итераций по всем узлам мы повторяем процедуру, пока не будет достигнута полная сходимость для выбранных значений ,, другими словами, пока значения не перестанут меняться. После полной сходимости итераций с данными , , уменьшим шаг варьирования вдвое и снова будем проводить итерации описанным выше образом. Будем считать, что приближенное решение на указанной сетке получено, когда будут обеспечены условия полной сходимости итерационного процесса для некоторого малого , удовлетворяющего условию . Для получения уточненного решения можно использовать более мелкую сетку с большим числом узлов за счет уменьшения величин , . Начальным приближением решения на новой сетке служит уже полученное ранее решение. Итерационный процесс повторяется для новой сетки с уменьшающимися значениями , как было описано выше. Процесс решения заканчивается при полной сходимости итераций для некоторых достаточно малых ,,, удовлетворяющих (18). 2. Упругопластическая задача Моделирование упругопластического кручения прямых стержней будем проводить на основе деформационной теории пластичности неупрочняющегося материала (для односвязной области). В качестве примера применения описанного выше алгоритма приведем решение задачи кручения упругопластического стержня, записываемой в виде где - модуль сдвига; - предел текучести при чистом сдвиге; - угол закрутки, приходящийся на единицу длины стержня; - односвязная область поперечного сечения стержня; - граница области , - функция сдвиговых напряжений, такая, что а - сдвиговые компоненты тензора напряжений соответственно в плоскостях (соответствует направлению по оси, перпендикулярной плоскости ). Минимум функционала (19) разыскивается на классе функций , удовлетворяющих условиям (20), (21). Если решение задачи (19)-(21) найдено, то скручивающий момент подсчитывается по формуле . Далее используются безразмерные переменные и обозначения где - характерный размер области , a - параметр задачи. В последующем тильда у безразмерных переменных опускается. В безразмерных переменных основные соотношения рассматриваемой задачи записываются в виде , , , . Решение задачи упругопластического кручения значительно упрощается, если воспользоваться следующим преобразованием. Вместо дифференциального неравенства (25) можно воспользоваться непосредственным неравенством наложенным на функцию состояния . Здесь решение более простой геометрической задачи чисто пластического кручения, описываемого следующими соотношениями: . Это так называемая аналогия Надаи [28,29], приводящая к построению «пластической крыши», описываемой функцией . Результаты вычислений представлены в таблице для прямоугольных областей поперечных сечений при различных углах закрутки стержня и для различных отношений длин сторон прямоугольников . Скручивающий момент μ для различных α и k Torsion moment μ for various α and k k α0 μ0 α = 20/3 α = 10 α = 20 α = 40 α = 80 1 2,9629 0,1025 0,1560 0,1624 0,1657 0,1664 0,1666 1,5 2,3585 0,1701 0,2780 0,2859 0,2903 0,2913 0,2916 2 2,1505 0,2418 0,3999 0,4093 0,4149 0,4162 0,4165 3 2,0305 0,3988 0,6418 0,6559 0,6640 0,6660 0,6665 4 2,0021 0,7285 1,1269 1,1490 1,1625 1,1654 1,1662 Приведем сначала качественное описание решения. При малых значениях ограничения не достигаются в области . В этом случае решение соответствует чисто упругому кручению. Пластические зоны появляются, когда Рис. 2. Зависимость момента кручения от угла закрутки Fig. 2. The torsion moment as a function of twist angle Этот угол закрутки находится из решения задачи об упругом кручении стержня. В интервале как следует из (28), имеем . При пластические области увеличиваются с увеличением . Момент кручения является монотонно возрастающей функцией (рис. 2), а соответствует чисто пластической закрутке (). На рис. 3 и 4 показано развитие пластических областей и 3D-графики функции для прямоугольных областей c и . При проведении расчетов использовалась сетка размером 100x100, а шаг варьирования изменялся от до . В качестве начального приближения выбиралась Численное решение находилось для задачи с ограничением, заданным выражением (30). Полученные пластические зоны обозначены зеленым цветом, упругие - белым. Приведенные ниже рисунки для функции напряжений показывают, как при выполнении расчетов искомая функция напряжений выходит на наложенные на нее ограничения, и как следствие зоны выхода функции на ограничения являются зонами возникновения пластичности (обозначены темно-зеленым цветом) Зоны, где ограничения не достигаются, зоны упругого поведения материала обозначены светло-зеленым цветом. x1 x1 x2 x1 x2 Рис. 3. Пластические зоны и 3D-графики функции в случае и Fig. 3. Plastic zones and 3D-graphs of function for and x2 x1 x2 x1 x1 x2 Рис. 4. Пластические зоны и 3D-графики функции в случае и Fig. 4. Plastic zones and 3D-graphs of function for and Рассмотрим задачу упругопластического кручения стержня круглого поперечного сечения , имеющего радиальный разрез . Здесь используются цилиндрические (полярные координаты) . Выражение, определяющее «пластическую крышу», т.е. функцию , имеет вид [30] и функция (рис. 5) для . Рис. 5. Функция пластических напряжений Fig. 5. Function of plastic stresses Решение неклассической вариационной задачи минимизации функционала производилось на классе функций, удовлетворяющих ограничениям где функция определялась формулой (37), включающей зависимость от угла закрутки . (Функционал (38) приведен в безразмерной форме по аналогии с (29)). Численное решение с применением описанного в работе алгоритма разыскивалось на конечно-элементной сетке: Все вычисления производились при а шаг варьирования изменялся от до . В качестве начального приближения выбиралась . Полученные при этом распределения пластических (зеленый цвет) и упругих областей (белый цвет) показаны на рис.6 для значений . При расчете взаимодействие берегов разреза не учитывалось, и контактная задача с трением не решалась. Рис. 6. Распределения пластических и упругих областей при Fig. 6. Distributions of plastic and elastic regions at Полученные численные решения соответствуют результатам, приведенным в [30]. Сравнение по времени решения упругопластических задач кручения, рассмотренных в данной работе, с методами, приведенными в [30], показало существенное сокращение необходимого времени (3-5 раз) с учетом многих факторов реализации алгоритма. 3. Некоторые замечания и обобщения Выше была описана простейшая версия сочетания конечно-элементных аппроксимаций и метода локальных вариаций. Рассмотрим кратко некоторые возможные модификации предложенного алгоритма. Так, например, в рамках предложенного алгоритма возможны как другие способы дискретизации (триангуляции) области определения функционала, так и другие полиномиальные базисные функции (функции формы). Шаг варьирования на рассматриваемой текущей итерации может меняться при переходе от одного узла к другому. При этом могут рассматриваться алгоритмы с локально оптимальными шагами варьирования , выбираемыми из условия максимального убывания функционала по четырем ячейкам, имеющим общий узел в точке . Отметим теперь некоторые обобщения формулировки задачи (1)-(4), для которых описанный алгоритм, основанный на конечно-элементных аппроксимациях и использовании локального варьирования, может эффективно применяться. 1. Искомая функция , входящая в соотношения (1)-(4), может быть вектором-функцией с компонентами В этом случае улучшающие вариации выполняются последовательно для каждой компоненты при фиксированных значениях остальных компонент. 2. Описанный подход применим к вариационным задачам, в которых интегральный функционал (2), (3) зависит от производных более высокого порядка от искомой функции. Для вычисления этих производных и выполнения локальных вариаций функции в одном узле сетки требуется вычисление большего числа слагаемых в сумме (12) или использование полиномиальных базисных функций более высокой степени. Последнее означает применение конечно-элементных аппроксимаций. 3. Данный алгоритм применим к многомерным вариационным задачам, в которых искомая функция зависит от трех и большего числа переменных (), а интеграл (3) берется по области соответствующей размерности. В этом случае алгоритм решения вариационной задачи модифицируется естественным образом. В качестве элементарных ячеек принимаются параллелепипеды соответствующей размерности, а при локальной вариации в одном узле рассматриваются изменения интегралов по большему числу ячеек (), лежащих в окрестности данного узла. 4. Предложенный алгоритм легко адаптируется для нерегулярных сеток разбиения области, в том числе для сеток с переменным шагом разбиения.

About the authors

N V Banichuk

Ishlinsky Institute for Problems in Mechanics of the Russian Academy of Sciences; Moscow Institute of Physics and Technology (State University)

E V Makeev

Ishlinsky Institute for Problems in Mechanics of the Russian Academy of Sciences

References

Statistics

Views

Abstract - 187

PDF (Russian) - 55