![]() Главная страница Случайная лекция ![]() Мы поможем в написании ваших работ! Порталы: БиологияВойнаГеографияИнформатикаИскусствоИсторияКультураЛингвистикаМатематикаМедицинаОхрана трудаПолитикаПравоПсихологияРелигияТехникаФизикаФилософияЭкономика ![]() Мы поможем в написании ваших работ! |
Входная информация для самопроверки. Для изучения данной темы Вам необходимо восстановить в памятиДля изучения данной темы Вам необходимо восстановить в памяти: - из курса прикладной математики - понятия и операции: линейной алгебры, векторной алгебры; понятия: интеграла, кратного интеграла, минимума функции нескольких переменных; - из курса сопротивления материалов – понятия: напряжений, деформаций, обобщённых усилий, работы обобщёных внутренних сил, работы нагрузки; деформативные свойства материалов; -из настоящего спецкурса – понятия: функционала, координатной функции, координатной системы, линейной аппроксимации, функционала энергии, невязки приближённого решения, ортогональности функций, последовательности решений, идея метода Ритца, идея метода Бубнова-Галёркина, идеи методов линеаризации. . 7.2. Содержание темы
Метод конечных элементов.Методы, изложенные в лекции ММ-6, которые мы будем называть традиционными прямыми методами, в своё время сыграли значительную роль в численном решении ряда важных практических задач. Однако, при использовании этих методов расчётчики столкнулись с рядом трудноразрешимых проблем. Перечислим некоторые из них. 1. Подбор координатных систем для неканонических областей. Расчётная схема конструкции может представлять собой область одно-, двух- и трёхмерного пространства. Так, расчётная схеме стержня - отрезок прямой линии - является областью прямой, то-есть, одномерного пространства; расчётная схема изгибаемой плиты или нагруженной в своей плоскости пластинки является областью плоскости, то-есть, двумерного пространства; равно, как и областью криволинейной поверхности (то-есть, также двумерного пространства) является оболочка. Расчётная схема массива (например, столбчатого фундамента или основания) - это область трёхмерного пространства. Область двух- или трёхмерного пространства будем относить к какой-либо фиксированной системе координат и называть канонической, если её граница состоит из участков, каждый из которых проходит: - для двумерной области - по координатной линии -для трёхмерной области - по координатной поверхности Для прямоугольной декартовой системы координат на плоскости это будет прямоугольник, в пространстве - паралеллепипед; для полярной системы координат - часть сектора, выделенная сторонами 2. Скачки или резкие изменения переменных коэффициентов дифференциального уравнения задачи.При расчёте конструкций это обстоятельство возникает, например, в случаях скаков или резких изменений свойств материалов или геометрии конструкций. При этом следует ожидать скачков или резкого изменения производных функции, описывающей перемещения, а придать координатным функциям способность отражать наличие таких скачков или резких изменений невозможно, так как эти функции должны иметь несколько непрерывных производных порядка не менее заданного. 3. Скачки или резкие изменения правой части неоднородного дифференциального уравнения задачи.Это обстоятельство возникает, например, когда интенсивность распределённой нагрузки описывается функцией, имеющей такой характер, вплоть до такой абстракции, как сосредоточенная сила или сосредоточенный момент. Сосредоточенной силе соответствует случай бесконечно большой интенсивности нагрузки на бесконечно малой области её локализации при конечном значении равнодействующей. Сосредоточенный момент моделируется парой состредоточенных сил с бесконечно малым плечом (величина момента приэтом конечна). . Для строгого описания таких нагрузок требуется аппарат обобщённых функций, но это выходит за рамки нашего спецкурса. У координатных функций такие особенности должны иметь их производные хотя бы того же порядка, что и порядок уравнения задачи, то0есть, мы опять сталкиваемся с проблемой несоответствия ожидаемых свойств решения и свойств координатных функций. 3. Сложность задачи конструирования координатной системы в общем случае.Для решения этой задачи в общем случае требуется достаточно высокая квалификация, она для большого количества классов задач не поддаётся стандартизации и требует эвристического подхода. 4. Трудность учёта различных краевых условий.При наличии таких краевых условий, как податливые опоры, задании на одном участке граници нескольких типов опор, применении локальных опор, изменении нагрузок на край в пределах участка границы и в ряде других случаев конструирование функции
Эти проблемы любой степени сложности достаточно легко преодолеваются в методе конечных элементов (МКЭ) – едва ли не самом мощном и универсальном численном методе для решения в первую очередь краевых задач для дифференциальных уравнений (обыкновенных и в частных производных) и систем таких уравнений. Основная идея метода конечных элементов. Метод конечных элементов основан на использовании нового принципа конструирования координатных систем. В традиционных вариантах прямых методов координатные функции задавались на всей области определения оператора задачи, границы которой задавались краевыми условиями и были на этой области достаточно гладкими. Это как раз и порождало практически все перечисленные выше противоречия. В отличие от этого в МКЭ используются принципиально новые координатные функции. Их основной отличительный признак состоит в том, что они удовлетворяют заданным условиям не на всей области определения оператора задачи, а на её достаточно малой (но всё же конечной) подобласти, за пределами которой они тождественно равны нулю. Иными словами, можно сказать, что они заданы на этой малой подобласти.
Таким образом, носителем каждой координатной функции метода конечных элементов является замыкание подобласти, на которой она определена, другими словами, объединение подобласти и её границы.
![]() Таким образом, в качестве координатных функций в МКЭ выбираются финитные функции. Воосоздание искомой функции на её носителе в МКЭ выполняется по а) дифференциальному уравнению (системе дифференциальных уравнений) или вариационному принципу, как в традиционных прямых методах; б) краевым условиям в отдельных точках границы носителя, а не на всей границе, как это было в традиционных прямых методах. Точки границы носителя, в которых задаются краевые условия, называются узлами. Как правило, в узлах задаются главные краевые условия. Их достаточно для определения искомой функции (в случае системы уравнений – искомой вектор-функции). Примечание. Иногда рассматриваются случаи расположения узлов внутри носителя; мы такими случаями не занимаемся. Отсюда естественным образом можно перейти к ключевому в МКЭ понятию – конечный элемент (КЭ).
КЭ могут иметь различную форму, и это позволяет более точно (хоть в ряде случаев всё равно приближённо) отразить форму области на расчётной схеме, а также некоторые априори известные особенности искомого решения. Примеры конечных элементов для расчётных схем конструкций различной размерности приведены на рис. 7.1. Здесь здесь на рис. 7.1а показан КЭ одномерной области (стержня), на рис. 7.1б – прямоугольный КЭ двумерной области (изгибаемой плиты или плосконапряжённой пластины), на рис. 7.1в – их же треугольный КЭ, на рис. 7.1г, д – КЭ трёхмерной области (массива) в виде паралеллепипеда или тетраэдра. На этих рисунках 1 – узлы, 2 – границы, 3 – рёбра границ. В криволинейных системах координат границы КЭ также удобно принимать криволинейными. Граница одномерного КЭ – всегда точка; она же является узлом. Граница двумерного КЭ - линия, которая состоит из гладких, обычно, прямых, участков; разделяющая их точка называется узлом. Граница трёхмерного КЭ – поверхность, которая состоит из гладких, часто - плоских, участков. Линии, разделяющие эти участки, называются ребрами; точка пересечения рёбер называется узлом. Узел трёхмерного элемента, таким образом, это, чаще всего, вершина многогранника (например, паралеллепипеда, тетраэдра), двумерного – вершина многоугольника (например, прямоугольника, треугольника). Таким образом, конечный элемент – ограниченный носитель координатной функции, узел конечного элемента является узлом носителя координатной функции. Можно сказать, что у одномерного КЭ понятия узла и границы совпадают; у двумерного КЭ границей является замкнутая линия, имеющая изломы (скачки касательной), и эти изломы являются узлами; у трёхмерного КЭ границей является замкнутая поверхность, эта поверхность имеет изломы, называемые рёбрами, а изломы рёбер называются узлами. Повторим, что для воссоздания координатной функции в пределах КЭ необходимо указать уравнение или вариационный принцип задачи и значения функции в узлах.
Для выполнения расчёта методом конечных элементов вся область, на которой задана краевая задача (вся область, занятая конструкцией) разбивается на отдельные КЭ таким образом, что - вся область занята КЭ без разрывов; - любые два КЭ либо не имеют общих точек, либо имеют общий участок границы (включая узлы), который может выродиться в одну точку (узел); никакие внутренние точки двух и более КЭ не могут быть общими (пересечение двух КЭ либо пусто, либо состоит только (!) из точек их общего участка границы). Можно сказать, что конечные элементы либо не соприкасаются, либо соприкасаются на их общем участке границы без нахлёста. Иногда используются конечноэлементные разбивки с зазорами между элементами, однако требуется, чтобы: 1) зазоры в узлах отсутствовали; 2) зазоры между узлами стремились к нулю при неограниченном сгущении узлов. В частности, такое разбиение на КЭ выполняется иногда при расчёте оболочек, очерченных по криволинейным поверхностям, с использованием плоских КЭ. Часто бывает удобным, чтобы участки границы двумерной области были направлены вдоль координатных линий (особенно, если таких участков четыре), участки границы трёхмерной области совпадали с координатными поверхностями, а рёбра – с координатными линиями. Обычно наиболее сложным вопросом является разбивка двумерной или трёхмерной области задания задачи в пределах её приопорной зоны, так как границу области часто невозможно составить из участков границы КЭ. В таких случаях приходится приближённо аппроксимировать границу области участками границы КЭ, но так, чтобы зазор между ними стремился к нулю вместе с размером самого большого КЭ. При этом внешние узлы КЭ надо располагать на границе области. Таким образом мы получаем конечноэлементную аппроксимацию границы области. У КЭ не может быть узлов, которые бы не совмещались с узлами соседних элементов или узлами конечноэлементной аппроксимации границы области. При этом в процессе разбивки заданной области на КЭ на указанной области выделяются точки, на которых совмещаются либо узлы соприкасающихся КЭ, либо узлы КЭ и границы; эти узлы называются узлами рассчитываемой системы (точнее, её расчётной схемы). Другие правила разбивки сетки узлов будут изложены после описания МКЭ.
Рис.7.2. Пример разбивки плоской области на конечные элементы
![]()
Если вычертить схему разбивки области на КЭ, отметить на ней узлы, занумеровать и те, и другие, а затем указать на КЭ информацию о свойствах системы (отражённых в коэффициентах этого уравнения), на узлах - данные о взаимодействии с внешней средой (отражённые в свободном члене дифференциального уравнения и в краевых условиях), то мы получим конечноэлементную расчётную схему системы. Методы определения этих данных будут сообщены ниже. Нумерация КЭ и узлов позволяет организовать ввод остальной информации в компьютер. Обычно широко используемые программы, реализующие МКЭ, ориентируясь на конкретную область приложений (прочностные расчёты, расчёты теплопередачи, расчёты электрических или гидравлических сетей и т.п.) содержат библитеки КЭ, позволяющие вычислить определённым образом организованные данные о естественных краевых условиях КЭ при заданных его главных краевых усовиях для различных типичных часто встречающихся случаев. При выполнении конкретного расчёта необходимо выбрать тип КЭ, проставить его размеры и указать величины параметров, характеризующих физические свойства материала. Информация о главных и естественных краевых условиях КЭ по определённым правилам приведена к его узлам – в каждом узле указывается некоторая функция искомого решения и его производных, обычно – сама искомая функция и её производные.
Одна из центральных идей МКЭ состоит в том, что приближённое решение должно полностью определяться конечным количеством параметров, указываемых в каждом из конечного числа узлов системы. Таким образом, мы состояние континуальной системы описываем дискретным множеством параметров; в этом и состоит дискретизация задачи в МКЭ. Количество параметров, указываемых в узлах, определяется порядком дифференциального уравнения и применяемым методом поиска коэффициентов линейной аппроксимации решения (в нашем случае – метод Ритца или метод Галёркина).. Каждому параметру каждого узла ставится в соответствие одна координатная функция. Она равна единице в этом узле и нулю - в остальных узлах ЕЭ. Благодаря этому а) линейная аппроксимация искомой функции включает слагаемое в виде произведения параметра на значение координатной функции, то-есть, узловой параметр является коэффициентом этой линейной аппроксимации и одновременно значением аппроксимации в узле; б) в других узлах параметр умножается на нулевое значение координатной функции, то-есть, не влияет на значение аппроксимации в других узлах. Действительно, выделим одно слагаемое линейной аппроксимации (см. лекцию ММ-6, формула (6.4))
где Отсюда легко заметить, что, если линейная аппроксимация уже найдена, то узловой параметр данного узла равен значению аппроксимации в этом узле и не влияет на значения аппроксимации в других узлах. В случае одномерной области в качестве параметра могут выступать искомая функция, её обыкновенные производные и их линейные комбинации, а в случае двумерных и трёхмерных областей – функции и её частные производные, в том числе – смешанные, а также их линейные комбинации. Наибольший порядок указанных производных на единицу меньше порядка дифференциального уравнения в методе Галёркина или порядка дифференциального оператора в методе Ритца (если он сформулирован в энергетически сопряжённых переменных - момент - кривизна, нормальная сила - осевое удлинение и т.п.). В случае применения интегрирования по частям при использовании метода Галёркина этот порядок может быть уменьшен вдвое. При наличии общего участка границы у двух смежных КЭ координатная система должна обеспечить возможность приравнивания хотя бы главных краевых условий данной задачи. В задачах строительной механики обычно это перемещения точек границы и (при наличии изгибного состояния) повороты нормали в этой точке вокруг границы. Таким образом, в условиях сопряжения КЭ участвуют искомая функция и (при наличии изгиба) её производные вдоль нормали к границе. Для того, чтобы обеспечить возможность приравнивания производных, взятых по нормали к границе в сопрягаемых КЭ, на всём их общем участке границы, приходится использовать в качестве узловых параметров смешанные производные. Однако, это требование иногда нарушается и возникают так называемые несогласованные координатные функции. Ещё одно условие, накладываемое на координатные функции в МКЭ: при стремлении размеров КЭ к нулю они должны обеспечивать возможность соответствующим подбором узловых параметров получать любые постоянные (в том числе – нулевые) значения искомой функции и её производных (порядка менее порядка дифференциального оператора) на КЭ и его границе. Часто наиболее удобными для использования в качестве координатных функций являются многочлены. В случае одномерной области полные многочлены степени Каждая координатная функция имеет вид многочлена и его коэффициенты должны обеспечивать выполнение заданного набора граничных условий в узлах (когда один из узловых параметров равен единице, а осталбные – нулю, а данная координатная функция равна единице на одном узле и нлую – на остальных). Для определения его коэффициентов в простых случаях удобно использовать метод неопределённых коэффициентов – записывается общий вид такого многочлена и определяются дифференцированием все его производные, участвующие в формулировании граничных условий КЭ (если производные участвуют в формулировании граничных условий). Это позволяет записать граничные условия в узлах КЭ, выраженные через общий вид многочлена и его производных, и получить систему Следует отметить, что в различных источниках при изложении МКЭ координатные функции называются также функциями формы или функциями распределения. Пример 1. Дано: одномерная область, порядок дифференциального уравнения Общий вид многочлена первой степени
Для первой функции:
откуда Для второй функции
откуда
Пример 2. Дано: одномерная область, порядок дифференциального уравнения Общий вид кубического многочлена Первый из них удовлетворяет условиям
второй –
третий –
четвёртый –
Составляем поочерёдно для каждой из этих координатных функций систему уравнений, отражающих граничные условия на узлах:
Вектора решений этих систем представим в виде таблицы (табл. 7.1).
Таблица 7.1.
Введём удобное для анализа свойств координатных функций относительное независимое переменное
Пример 3. Дано: двумерная область, порядок дифференциального уравнения по каждой координате
Первая из них удовлетворяет условиям вторая - третья – четвёртая – Составляем поочерёдно для каждой из этих координатных функций систему уравнений, отражающих граничные условия на узлах:
Вектора решений этих систем представим в виде таблицы (табл. 7.1). Таблица 7.1.
Введём удобные для анализа свойств координатных функций относительные независимые переменные
Работа с несогласованными координатными функциями требует, как правило, более глубокого анализа, иногда с физическим осмыслением пренебрегаемых условий сопряжения, поэтому примеры таких функций мы не приводим.
Существенно, что в наших задачах количество главных и естественных краевых условий одинаково и это энергетически сопряжённые величины – интеграл от произведения естественного условия на вариацию главного условия равен вариации работы обобщённого усилия (естественного условия) на соответствующем возможном обобщённом перемещении (главном условии). В процессе расчёта подбираются такие значения главных краевых условий в узлах (одинаковые для всех узлов КЭ, которые совмещены с данным узлом системы), чтобы в этих же узлах выполнялись (в приближённом смысле) естественные краевые условия. В задачах строительной механики это означает, что узлы должны находиться в равновесии, при этом условия равновесия могут записываться в трёх эквивалентных вариантах: либо в виде равенства нулю обощённой силы, действующей на узел, либо в виде равенства нулю работы обобщённой силы на соответствующем обобщённом перемещении, либо как условие минимума функционала энергии. Если среди краевых условий есть главные и естественные, то, как мы выяснили, главные условия должны выполняться точно, а естественные – приближённо. В прямых методах это означает, что главные условия должны выполняться уже во время конструирования координатной системы, а приближённое соблюдение естественных уловий должно обеспечиваться в результате расчёта..Исходя из этого, условия сопряжения подобластей в узлах обеспечивается: - строгим равенством величин, задающих главные условия (в строительной механике – обобщённых перемещений нормали и её поворотов для балок и плит и перемещений для массивов); - возможностью выполнения любых заданных естественных условий в форме условий сопряжения соответствующим подбором величин, при помощи которых формулируются главные условия; связь между главными и естественными условиями может быть приближённой (из-за того, что используется только конечное число координатных функций), но любому сочетанию естественных условий должно соответствовать создающее его сочетание главных условий. Строгое равенство главных условий в налагаемых друг на друга узлах смежных элементов (можно сказать – в совмещаемых узлах смежных элементов) обеспечивается в алгоритме (и программе) тем, что неизвестными являются приближённые значения искомой функции в узлах системы, и этим значениям приравниваются узловые параметры КЭ для их узлов, попадающих в эти узлы системы. Рассмотрим, какой вид приобретают коэффициенты системы уравнений методов Ритца или Галёркина, когда эти методы сочетаются с МКЭ. Согласно традиции, установившейся в литературе по МКЭ, будем обозначать эти коэффициенты как Вначале выберем КЭ, занимающий всю рассматриваемую область. Введём обозначения:
Этот интеграл в случае неодномерной области – кратный, Количество элементов множества
![]() В ряде случаев различные КЭ одной конечноэлементной сетки могут быть сориентированы по-разному. В этих случаях целесообразно вводить на каждом КЭ свою систему координат. Такая система координат называетсяместной(локальной) системой координат данного КЭ. Если ввести общую(глобальную) систему координат, то взаимное положение этой системы и каждой из местных систем может быть охарактеризовано углом взаимного поворота
где В аналитической геометрии доказывается, что
Тогда из (7.4) с учётом (7.6) можно записать
Подставим (7.7) в (7.3), после чего - преобразованное выражение (7.3) в (7.5)
где
В МКЭ принято определять естественные граничные условия КЭ, приведенные к его узлам (узловые обобщённые силы КЭ в случае задачи строительной механики)
Теперь перейдём к случаю сетки КЭ с произвольным количеством элементов. Узлы, принадлежащие какому-либо КЭ либо совпадающие с его узлами, в МКЭ называют инцидентнымиему. Соответствие между номерами КЭ и номерами инцидентных им узлов системы задаётся расчётчиком в исходных данных в виде специальной таблицы инциденций. Ненулевому значению любого узлового параметра узла системы, инцидентного нескольким КЭ, соответствует в этих КЭ аппроксимация решения
где Пусть теперь в некотором узле системы с глобальным номером
где Учитывая, что за пределами носителя координатные функции тождественно равны нулю, (7.11) можно записать в виде
По аналогии с формулами (6.30), (6.38) и сучётом (7.8) десь следует иметь иметь ввиду,запишем блоки матрицы коэффициентов МКЭ в форме методов Ритца или Бубнова-Галёркина (опять согласно традиции МКЭ заменим обозначение коэффициентов на
Здесь Для столбцов свободных членов аналогично (7.9) получаем выражение для
Для матрицы коэффициентов всей системы уравнений для определения узловых параметров рассчитываемого объекта и её столбца свободных членов естественно получаем, используя (7.13) и (7.14)
Упомянутая система уравнений в матричной записи имеет вид
Это разрешающая система уравнений МКЭ.
![]()
Обычно при рациональной нумерации узлов и конечных элементов разность Решая систему (7.17), получаем приближённые значения искомой функции в узлах системы
где Её производные порядка до
где Обычно необходимы не только значения искомой функции в отдельных узлах системы, получаемые в результате решения системы (7.17), но также и значения функции и её производных (или их линейных комбинаций) между узлами, то-есть, внутри КЭ. В этих случаях для каждого КЭ строятся аппроксимации искомой вектор-функции (7.11). В случае решения граничных задач для систем уравнений с несколькими неизвестными функциями идея МКЭ сохраняется, но требуется соответствующее обобщение расчётных соотношений. В теории МКЭ доказывается, что при правильной разбивке сетки КЭ погрешность метода падает при уменьшении диаметра (наибольшего размера) самого большого КЭ, или сгущении узлов системы, и приближённое решение стремится к точному. При рациональной разбивке сетки узлов адекватность расчётной схемы МКЭ растёт при сгущении узлов, в первую очередь - в окрестности источников возмущения решения (входящие углы, проёмы, местные воздействия, резкие изменения геометрии системы или свойств материала ит.д.). Следует подчеркнуть, что предварительная (априорная) оценка погрешности практически невозможна, и единственная возможность получить приближённое решение с погрешностью,не превосходящей заданную – организавать последовательнось процессов решения задачи при сгущении узлов и измельчении КЭ от расчёта к расчёту. При грамотной подготовке исходных данных последовательность полученных приближённых решени должна сходится к точному решению, а разность двух решений может послужить для оценки погрешности (апостериорной), если, например, сгущение узлов осуществлять только путём добавления новых узлов и сравнивать вектора значений неизвестнойо функции в тех узлах, которые есть в нескольких последовательных сетках узлов. Методы решения нелинейных задач в строительной механикеприобретают специфические формы. Рассмотрим их на примере нелинейных задач, в которых дискретизация выпоняется методом Бубнова-Галёркина в форме МКЭ. Общим для этих задач является зависимость матрицы жёсткости системы от решения Уравнение (7.17) принимает вид
Конкретизируем пути применения изученных нами ранее методов линеаризации (лекции ММ-4 и ММ-5). 1. Метод простых итераций.Запишем (7.20) в виде
Это попрежнему уравнение, так как оно не разрешено относительно
В механике деформируемого твёрдого тела, в строительной механике такой итерационный процесс называется методом переменных параметров упругости, так как он предусматривает пересчёт коэффициентов уравнения (7.20) на каждой итерации в соответствии с заданными зависимостями напряжений от деформаций. Естественно, приведение к виду (7.21) нам понадобилось только для пояснения связи между понятиями метод простых итерацийи метод переменных параметров упругости. При выполнении реальных расчётов записывать решение в таком виде нецелесообразно, нужно просто решать систему уравнений, ставших линейными после подстановки предыдущей итерации в выражение для
2. Модифицированный метод Ньютона.Использование обычного метода Ньютона при выполнении реальных расчётов практически невозможно из-за отсутствия аналитических выражений для отмеченных выше нелинейных зависимостей и необходимости их определения в большом количестве точек – узлов системы. Его замена аналогом – многомерным вариантом метода секущих в отмеченных случаях также чрезмерно трудоёмка. В то же время модифицированный метод Ньютона в ряде случаев может быть применён (и действительно применяется). Запишем (7.20) в виде
где Найдём матрицу Якоби левой части (7.24)
Её значение в точке естественного состояния (
(7.26) Теперь мы можем записать итерационный процесс модифицированного метода Ньютона
Введём обозначения
Таким образом, итерационный процесс модифицированного метода Ньютона принял форму метода фиктивных нагрузок Ильюшина
где Иногда этот метод называют методом упругих решений, но в настоящее время чаще это название употребляют для обоих описанных выше методов. Применяемые иногда методы начальных напряжений и начальных деформаций также приводятся к виду метода фиктивных нагрузок. 3. Шагово-итерационные методы. Итерационный метод может применяться в случае, если порождаемый им итерационный процесс в условиях решаемой задачи сходится. В общем случае одним из необходимых условий, обеспчивающих сходимость, является принадлежность начального приближения области сходимости неподвижной точки оператора задачи (то-есть, решения данной нелинейной задачи). Кроме того, при выполнении практических расчётов важно, чтобы итерационный процесс сходился достаточно быстро. В то же время, чем выше нагрузка и, следовательно, чем выше уровень напряжений и деформаций, тем сильнее отклоняются упомянутые выше нелинейные зависимости от линейных и тем меньшей и труднее оцениваемой оказывается область сходимости и тем медленнее сходится вычислительный процесс. Для ускорения этого процесса нужно начальное приближение располагать как можно ближе к искомому (неизвестному заранее!) решению. Этого можно достигнуть, если в качестве начального приближения при данном уровне нагрузки использовать пусть приближённое, но достаточно точное решение, которое можно получить при несколько меньшей нагрузке. Можно построить процесс постепенного пошаговогоувеличения нагрузки, начиная со столь малой нагрузки, что начальным приближением для неё могло бы служить решение, полученное в предположении линейноупругой работы материала. При этом окончательное (на последней итерации) решение, полученное на предыдущем шаге нагрузки, принимается в качестве начального приближения для последующего
Дата добавления: 2014-04-15; просмотров: 472; Нарушение авторских прав ![]() Мы поможем в написании ваших работ! |