Метод конечных элементов (1061787), страница 8
Текст из файла (страница 8)
Во-вторых, должно оставаться справедливым суммирование ~3.4) и, если не добавлены поверхностные интегралы по границе между элементами [71, мы должны быть уверены, что величины, подобные 1" и д, остаются на них ограниченными. Это возможно, если производные наивысшего порядка от Ф, входящие в выражения для функций ~ и д, конечны.
Таким образом, можно сформулировать второй критерий. Критерий 2. Функции формы элемента 1Х1 должны быть выбраны так, чтобы ® и ее производные, на порядок более низкие, чем- производные, входящие в выражения для 1 и д, были непрерывны на границе между элементами. Этот критерий несложно объяснить, если представить себе что элементы разделяются между собой очень тонким слоем, который должен быть учтен в определяющих величину ~ интегралах и в котором происходит плавный переход между значениями неизвестных функций в смежных элементах. Если в пределе при стремлении толщины этого переходного слоя к нулю его вклад в ~ исчезает, то равенство 134) справедливо.
На фиг. 3.1 показана такая воображаемая переходная зона между двумя элементами, Представим себе, что скалярная функция ф определена так, что на границе раздела между элементами ее значения, полученные для каждого из элементов, равны. На фиг. 3.1,а приведен график функции, угол наклона которой принимает в переходной зоне конечное значение, хотя и претерпевает разрыв (фиг. 3.1, б).
Вторая производная в этой зоне имеет очень большое значение (фиг. 3.1,в) и стремится к бесконечности при уменьшении ширины зоны„ Обоби~ение понятия конечных элечентое Следовательно, если в выражение для т входит только первая производная, то, для того чтобы гарантировать отсутствие вклада в т от переходной зоны, достаточно обеспечить непрерывность только функции ф.
Однако если в функционал входит вторая производная, то справедливость равенства (3.4) не может быть гарантирована, поскольку (из-за умножения бесконечной величины на нулевую площадь) величина вклада в энергию этой зоны становится неопределенной. 31 Стрешипкя к- с Фиг. 3.1. Межэяеиентная зона, в которой вторая производная непрерывной фукции Ф при Лх-~-О становится неограниченной. Эти два критерия являются обобщением более частных критериев разд. 2.5, Здесь под сходимостью мы понимаем следующее: при бесконечном уменьшении размеров элементов ошибки в определении величины т стремятся к нулю, Это иногда позволяет сказать, что решение, полученное при одном разбиении, заведомо лучше решения, полученного при другом. Очевидно, что в смысле определения величины т 1формула (3.2)1 это утверждение справедливо, если функция формы первого типа разбиения включает все Функции 4ормы второго типа разбиения. Именно .такой случай возникает, когда новое разбиение получается последующим делением более крупных элементов.
Сходимость по ~ при этом будет монотонной. Это обстоятельство впервые было установлено Мелошем в 1963 г. 181. Глава 8 3.3. Переменные, не связанные с узлами Важно напомнить читателю, что при определении функции формы [равенство (3.2)~ указывалось, что (Ф)' содержит либс узловые значения неизвестных функций, либо некоторые параметры, характеризующие этот элемент. Для обеспечения непрерывности функции между элементами в соответствии с требованием критерия 2 (разд. 3,2) и описания физической сущности задачи обычно рассматриваются значения величин в узловых точках. Однако всегда можно ввести дополнительные функции, принимающие нулевые значения на границах элемента и не вызывающие нарушения непрерывности, и умножить их на некоторые параметры, по которым будет минимизироваться функционал.
Введение таких переменных, не связанных с узлами, может увеличить точность и иногда оказывается полезным Щ Поскольку в общем случае эти параметры связаны только с одним элементом, минимизация по ним может быть выполнена перед составлением ансамбля, а сами параметры исключаются из матриц элемента. Частным случаем неузловой переменной является хорошо известный множитель Лагранжа 131. Он вводится в том случае, когда на функцию ф наложено дополнительное условие, не вытекающее ни из краевых условий, ни из ограничений на функции формы.
Пусть дополнительное условие имеет вид 6(®) =О. (3.1О) В этом случае для решения задачи требуется минимизировать') величину х х+~6, (3.11) где Х вЂ” типичный дополнительный параметр. Если наложено несколько таких условий, то минимизируется величина ж'=-х+ Х к~а- (3.12) 1=1 ') С доказательством этого положения можно ознакомиться в книгах по вариационному исчислени1о. Очевидно, что дХ"/дх = 6, и, следовательно, ограничение является условием минимизации.
В точке экстремума т* = т, так как 6 = 0 (физический смысл Х становится ясным из соотношения Х = дХ*/дб). Теперь дополнительными параметрами задачи являются вели- чины Х;, которые опять можно связать с элементами или гра- ницами между элементами. Примеры использования таких мно- жителей Лагранжа будут приведены ниже. Обобщение понятия конечных элементов Интересно рассмотреть случай, когда функционал (3.1) представляет собой квадратичную форму и для ряда линейных ограничений на функцию Ф используются множители Лагранжа.
В общем случае можно выразить у через величины узловых параметров в виде х=-,РГ%1(Ф)+ ®'ФЬ (3.13) где Я вЂ” симметричная матрица. Вместо (3.10) линейные ограничения можно записать в матричной форме: И(%=0, (3.14) где ~6~ — постоянная матрица. Записывая множители Лагранжа в виде вектора 1Ц размерности, равной числу столбцов матрицы ~61 (т.
е. числу ограничений), имеем 1 (цт щ щ + (р)т (ф) + ( р~ р))т р~ (3 15) Составляя систему уравнений, аналогичную (3.6), но уже для двух множеств неизвестных, получаем д (й) ! ~Щ Щ~~ ((Ф~~ ~(Р1~ Здесь нужно отметить два обстоятельства.
Во-первых, система уравнений остается симметричной — факт, облегчающий применение обычных методов решения систем линейных алгебраических уравнений (гл. 20). Во-вторых, на диагонали матрицы появляются нули, что иногда затрудняет получение решения. В настоящей книге множители Лагранжа используются редко, исключение составляют случаи, когда приходится вводить некоторые ограничения.
Типичным примером этого может служить ситуация, когда при применении метода конечных элементов нарушаются условия непрерывности. В результате введения дополнительных условий в виде ограничений на параметры получается корректное решение. Однако в связи с тем, что такие ограничения приводят к увеличению общего числа неизвестных, появляются дополнительные трудности. 3.4. Другие подходы к методу конечных элементов Несмотря на то что приближенная минимизация функционала — самый распространенный способ подхода к методу конечных элементов, это никоим образом не означает, что такой Глава 8 подход является единственна возможным. Например, в первых работах по строительной механике строились чисто физические модели, и, хотя приходилось делать некоторые математические оговорки, касающиеся обоснования и сходимости использованных методов, зачастую получались неплохие инженерные решения.
Существуют и другие возможности, позволяющие математически получить основные соотношения метода конечных элементов непосредственно из дифференииальнь~х уравнений задачи, Они будут здесь кратко описаны. Возможные преимущества таких методов состоят в том, что: а) исчезает необходимость искать функциональный эквивалент известным дифференциальным уравнениям; б) зти методы могут быть распространены на задачи, для которых функционал либо вообще не существует, либо пока еще не получен 110). Рассмотрим задачу приближенного решения системы дифференциальных уравнений, которым должна удовлетворять неизвестная функция ® в области .К Запишем основное уравнение в виде А(®) =О, (3.17) а граничное условие на границе Я как С((Я) =О. (3.18) В'ЙЖ" =0„ (3.21) где К вЂ” любая функция координат.
Если число неизвестных параметров (Ф) равно и, то, выбрав а линейно независимых функций В';, запишем соответствующую систему уравнений ~ П7фИГ ~ Й~';А( ~й) (Ф)) ШГ =О, (3.22) Если пробная функция, удовлетворяющая граничным условиям, записана в общей форме ®. = ~л'ИФ), (3.19). где, как и прежде, [У~ является функцией координат, а 1Ф)— система а параметров„то в общем случае А(®,)=Я ~ О.
(3.20) Наилучшим решением будет то, которое дает во всех точках области Р наименьшую невязку Л Щ, Очевидно, это решение, можно получить, использовав то обстоятельство, что если невязка Л тождественно равна нулю всюду в области, то Обоби1ение понятия конечных элементов из которой может быть найдена функция (Ф), Зтот процесс называется методом взвешенных невязок, а йт; — весовой функцией. Выбор различных весовых функций приводит к различным классическим методам.
Коллокация в точке. В этом случае полагается, что В';=1 в некоторой точке ~ и рйвна нулю во всех остальных. При этом фактически основное дифференциальное уравнение удовлетворяется в и отдельных точках, Коллокация в подобласти. В этом методе считается, что 'йт; = 1 в некоторой подобласти и Ит~ — — О в остальной части области. Зто эквивалентно тому, что интеграл обращается в нуль в некоторых подобластях, число которых достаточно для того, чтобы получить необходимое число уравнений. Метод Галеркина, В этом случае 1й';= Уи т.
е. в качестве весовой функции выбирается функция формы, с помощью которой аппроксимируется решение. Зтот метод обычно приводит к наилучшим результатам. При использовании в любом из упомянутых методов соотношения (3.19), определяющего принятую аппроксимацию, можно выявить основные особенности метода конечных элементов. Во-первых, результирующая система уравнений будет иметь ленточный вид, так как влияние каждого параметра распространяется только на элементы, примыкающие к рассматриваемой узловой точке. Во-вторых (в предположении, что, как и ранее, границы между элементами пе дают никакого вклада), интегралы вычисляются для каждого элемента независимо, а затем полученные результаты суммируются. Очевидно, что правила получения коэффициентов для ансамбля будут такими же, как и в задачах строительной механики, если оператор 4 линеен 1см.
![](https://s.studizba.com/z.php?f=/templates/si/i/noavatar.png&w=100&h=100&t=1"){kind=avatar}
