Конечно-элементная реализация метода геометрического погружения на основе вариационного принципа Кастильяно для плоской задачи теории упругости Текст научной статьи по специальности «Физика»

2013 Механика № 1

УДК 539.3

H.A. Труфанов, Ю.С. Кузнецова

Пермский национальный исследовательский политехнический университет, Пермь, Россия

КОНЕЧНО-ЭЛЕМЕНТНАЯ РЕАЛИЗАЦИЯ МЕТОДА ГЕОМЕТРИЧЕСКОГО ПОГРУЖЕНИЯ НА ОСНОВЕ ВАРИАЦИОННОГО ПРИНЦИПА КАСТИЛЬЯНО ДЛЯ ПЛОСКОЙ ЗАДАЧИ ТЕОРИИ УПРУГОСТИ

Рассмотрен вариант метода геометрического погружения для плоских задач теории упругости, основанного на методе конечных элементов в напряжениях в рамках принципа минимума дополнительной работы упругой системы. Суть метода геометрического погружения заключается в сведении исходной задачи для линейно-упругого тела произвольной формы к итерационной последовательности задач теории упругости на некоторой канонической области. Сформулирована итерационная процедура для решения вариационного уравнения метода геометрического погружения, а также процедура построения его дискретного аналога с помощью метода конечных элементов в напряжениях для плоской задачи теории упругости в декартовой системе координат. Использован вариант конечного элемента в терминах функции напряжений для удовлетворения аппроксимирующих выражений уравнениям равновесия. Продемонстрировано практическое применение метода на примере решения плоской задачи для упругой пластины с прямоугольным вырезом. Получено достаточно хорошее соответствие результатов определения полей напряжений в сравнении с традиционным методом конечных элементов в перемещениях. Установлена практическая сходимость итерационной процедуры метода геометрического погружения. Уделено внимание способам задания статических граничных условий, являющихся главными для данной вариационной формулировки. Использован способ модификации матрицы податливости системы конечных элементов и метод множителей Лагранжа.

Ключевые слова: метод конечных элементов, вариационный принцип Кастильяно, метод геометрического погружения, итерационная процедура, граничные условия, метод множителей Лагранжа.

N.A. Trufanov, Y.S. Kuznetsova

Perm National Research Polytechnic University, Perm, Russian Federation

FINITE ELEMENT REALIZATION OF GEOMETRICAL IMMERSION METHOD ON THE BASIS OF KASTILYANO'S VARIATION PRINCIPLE FOR TWO-DIMENSIONAL PROBLEM THEORY ELASTICITY

The alternative geometrical immersion method for plane problems theory elasticity, based on a finite elements method in term of stresses within the principle minimal additional work of elastic system is considered. The geometrical immersion method consists in reducing an initial problem for linearly elastic body of any form to iterative sequence of problems theory elasticity on some initial area. Iterative procedure for the solution of the variation equation of the geometrical immersion method and a procedure of formulation of its discrete analog by means of the finite elements method in term of stresses for plane problems theory elasticity in the Cartesian system of coordinates are formulated. The alternative finite element method in terms of stresses function for satisfaction of approximating expressions to the balance equations is used. Practical application of the method on the example of the plane problems solution for the elastic plate with rectangular hole is shown. Rather good compliance of results of stresses fields definition in comparison with a traditional finite elements method in movements is obtained. Practical convergence of iterative procedure of the geometrical immersion method is established. We focus our attention to the problem of the static boundary conditions that are the main for this variation formulation. The way of modification of a flexibility matrix of the system of finite elements and Lagrange’s multipliers method is used.

Keywords: method of finite elements, Kastilyano's variation principle, method of geometrical immersion, iterative procedure, boundary conditions, method of multipliers Lagrange.

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

вательности задач для тела канонической конфигурации, при этом возможно применение конечных элементов в напряжениях простой формы [3].

1. Рассмотрим вариационную постановку задачи теории упругости в рамках вариационного принципа Кастильяно. Функционал дополнительной энергии упругой системы имеет вид

Пс (°)=1 ]Ч-% (<*)^|Ь ^ , (1)

2 V Би

где V - область, в которой разыскивается решение; ст, е - тензоры напряжений и деформаций с компонентами ст. {г), г. {г) соответственно;

г = (х,у,2)eV; и - заданные перемещения на поверхности Би; I■ - усилия на поверхности, где заданы кинематические граничные условия.

Функционал (1) определен на множестве статически допустимых полей напряжений ст. (г), удовлетворяющих в области V уравнениям

равновесия и статическим граничным условиям на части поверхности тела Я ( Я иБи = Б).

а V а и /

ст.. а. = Т., г е Б ,

V ] ‘ ’ о ’

где Т - заданные усилия на поверхности Ба; а. - компоненты вектора

внешней единичной нормали к поверхности тела.

Условие минимума функционала (1)

5ПС = 0

эквивалентно выполнению в области V уравнений совместности в напряжениях (уравнений Бельтрами-Митчела) и граничных условий в перемещениях на границе Би

Ц =и , г еБи.

Введем в рассмотрение некоторую каноническую область У0, включающую исходную область V произвольной формы, так что

V = V \JVA, где V - дополнение области V до области У0.

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

пс (6)=\& % (<*)+1 & % (а) й ^ - 1 ¡<3ц % (<?) й ^ -1 ь (<у)и, ,

2 V 2 уА 2 уА яш

пс (6)=1 Ч (д)- 2 (<^)й^д- / I . (2)

2 V 2 V Яи

Можно показать, что поля напряжений, минимизирующие функционал (2), в области V доставляют минимальное значение функционалу (1).

Условие минимума функционала (2) примет вид

шс (°)= |§ал- (°) - 15стл- (а)I (5дР> ^=0 (3)

V VA Би

Вариационное уравнение (3) будем решать, используя следующую итерационную процедуру:

{8а/*%(о^)dV-{8а/*%(о«)dVA-{I,(ба<м>)и(¡5,, = 0, (4)

V VA Би

о!0)= 0, к=1,2,3...,п,

где к - номер итерации.

Таким образом, на каждой итерации в (4) необходимо решать вариационное уравнение в канонической области У0, что можно с успехом сделать, используя существующие для этих целей формулировки МКЭ в напряжениях, например прямоугольные элементы для прямоугольных канонических областей [4]. Приведенный далее пример показывает, что итерационная процедура (4) сходится.

2. Рассмотрим процедуру построения дискретного аналога вариационного уравнения (4) с помощью метода конечных элементов в напряжениях для плоской задачи теории упругости в декартовой системе координат. Для дискретизации выберем прямоугольный конечный элемент с четырьмя узлами в углах (рис. 1).

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

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

{5у={ф./^| [VI

г г \ йх йу )\дхду )

мента описано в статье [4].

. Подробно использование данного эле-

Л

(0,1)

(0,0)

(1,1) ■о

(1,0) £

Рис. 1. Прямоугольный конечный элемент в локальной системе координат

Уделим особое внимание статическим граничным условиям, являющимся главными для данной вариационной формулировки. Для выполнения главных граничных условий можно использовать традиционный прием, заключающийся в модификации глобальной матрицы податливости, задавая значения функции напряжений и ее первых производных на границах. Такой подход вызывает затруднения, связанные с определением узловых неизвестных на границах, поэтому для удовлетворения статических граничных условий широко используется метод множителей Лагранжа [5, 6]. Применение данного метода иллюстрирует следующий пример: пусть на стороне элемента, параллельной оси Ох , действуют поверхностные силы Ту , как изображено на рис. 2.

Связь функции напряжений Эри и компонент тензора напряжений имеет вид

Э2ф

¿У’

уу

Э2ф

'~дх2.

ху

Э2ф

дхду

Ту ( * )

Рис. 2. Элемент, на сторону которого действуют поверхностные силы

Дважды проинтегрировав по переменной у на стороне элемента

• • д2ф „ / ч

узлами I и у выражение —2=Т (х) и определив константы интегриро-

с)зс

вания через ф и

V dх )

в узловых точках

Ф,-

Эф

V dх у

> Ф

р

,получим

V ¿х ) у

Эф | +Г 5ф

V йх )1

-ф.+ф .■

V йх )

= ¡Ту (х) ¿х

х

хУ х

а = \ТУ (х) йхйх.

(5)

(6)

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

[г М={4

(7)

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

хх

Согласно процедуре метода конечных элементов в напряжениях выражения, входящие в уравнение (4), можно представить для отдельного конечного элемента в матричном виде:

|5ст/+1)є.. (о*к+1)) с1У = 5{фе(к+1)}[ке

Vе

|5а/+% (й>к >) +1|}[к

Ф

г(к+1)

Ф

<к )

где

матрица податливости элемента, принадлежащего канони-

ческой области; к

матрица податливости элемента, принадлежа-

щего области дополнения; - вектор обобщенных узловых пере-

мещений элемента.

Проведя теперь типовую процедуру ансамблирования конечных элементов и учет статических граничных условий методом множителей Лагранжа, получим итерационную последовательность СЛАУ:

К №' И [о]

{<р< к *■>}'

{л<к

= М{ф(к >}

+І

{5}

где [К]0 =^[ке - глобальная матрица податливости канонической

е

области; [К]д=^[ке д - глобальная матрица податливости дополне-

е

ния; {ри}=^{р^ - глобальный вектор обобщенных узловых переме-

е

щений; - глобальный вектор узловых неизвестных; {А,} - вектор

множителей Лагранжа.

Условие остановки итерационной процедуры выбрано в виде

Ф

,(к)

< 0,001,

(8)

где в качестве нормы вектора |ф| выбрана норма Чебышева.

о

и

е

0

3. Рассмотрим задачу о плосконапряженном состоянии прямоугольной упругой пластины с прямоугольным вырезом при продольном ее растяжении переменной распределенной нагрузкой и заданными на контуре выреза компонентами вектора перемещений (рис. 3). Так как задача симметрична относительно осей Ох и Оу, будем рассматривать четверть пластины.

V

уа

(0,Ь)

А ~> —г

/ ^ —г

/ ^ —г

/ 1 —г

/

/ J —с.

(0, Ь/2)

(0,0)

Г( у)= 20

1-/

ь2

х

Рис. 3. Расчетная схема задачи

Граничные условия задачи:

их ( °> у ) = °> у е[Ь/2,Ь ],

иу (х,0) = 0, х е[а/2,а],

их (х, Ь/2) = 0,1Е -10 • х, иу (х, Ь/2) = 0,1Е -10, х е[0, а/2], их (а/2, у) = 0,4Е -10, иу (а/2, у) = 0,05Е -10 • у, у е[0, Ь/2],

°хх (a, у ) = Т (у ) , ^ху (a, у ) = 0, у е[0, Ь] ,

тху (х Ь) = ауу (х Ь) = 0, х е[0, а],

тху(0,у) = 0 уе[Ь/2,Ь], тху ( х,0 ) = 0, х е[ а/2, а ].

(9)

Для конкретизации зададим а = 8 м, Ь = 4 м, Е = 2Е + 11 Па - модуль Юнга, V = 0,35 - коэффициент Пуассона.

Согласно введенному выше условию V0 = V и ^ каноническая область будет представлять собой прямоугольник с координатами

вершин К0,0); (а,0); (а, Ь); (0, Ь)}. Следует отметить, что для канонической области граничные условия при х = 0 и у = 0 будут справедливы для всей стороны.

Построим итерационную процедуру для поставленной задачи в форме (4):

I)(8а „< ‘+,)в хх (а< ‘+1>) + бауу(‘+,)в уу (о< ‘+1)) + 5т „< ‘+1>у „ (а< ‘+1>)) <Му =

0 0

Ь/2 а/2

= I х}к+1)8хх(а(к}) + 5стуу(к+1)8уу(ст(к)) + 5тук+1)Уху(к}))¿хёу +

00 а/2 ч (10)

\ ( К (5°(к+1))^х ( х, Ь/ 2 ) + *у (5<^(к+1))иу ( х, V 2 )) ^ +

0

/2

\ (К (5°( к+1)) их (а12, у)+*у (5°(к+1)) иу (а12, у)) 4у-

0

0

Ь/2

+

Реализация соотношения (10) произведена методом конечных элементов, описанным в п. 2. Статические граничные условия (9) удовлетворяются с помощью метода множителей Лагранжа, граничные значения касательных напряжений задаются путем модификации глобальной матрицы податливости.

Все вычисления произведены в среде МЛТЬЛВ 7.0.1. Для выявления оптимальной сетки поставленная задача была решена методом конечных элементов в напряжениях. Исследование метода геометрического погружения, основанного на вариационном принципе Кастильяно, производилось при числе степеней свободы, равном 1284. Итерационная процедура сходится за 20 итераций при остановке по норме (8).

Поставленная задача также решена в программном комплексе ЛК8У8 с использованием формулировки в перемещениях при аналогичном числе узловых неизвестных (рис. 4).

Решение, полученное методом геометрического погружения в напряжениях, достаточно близко к решению поставленной задачи методом конечных элементов в перемещениях (рис. 5).

-6.408 -4.725 -3.041 -1.358 325346 2.009 3.692 5.376 7.059

%

.

-25.325 21.652 -17.978 14.305 -10.631 6.958 -3.285 388714 14.062

, Па

I___X

9.702 -7.972 6.242 -4.512 2.782 -1 .052 678025 2 .408 4 .138

б

Рис. 4. Распределение компонент тензора напряжений, полученных методом конечных элементов (ЛКБУБ) на основе вариационного принципа Лагранжа: а - касательное напряжение тху ; б - нормальное напряжение а;

в - нормальное напряжение ахх

а

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

Рис. 6 иллюстрирует характер сходимости итерационной процедуры метода геометрического погружения при выборе нулевого начального приближения для компонент тензора напряжений. Для сравнения приведены эпюры напряжений, найденные методом конечных элементов в перемещениях (кривые с крестиками). Результаты, представленные на рис. 6, соответствуют теоретическим представлениям, согласно которым решения, полученные на основе вариационных

принципов Лагранжа и Кастильяно, образуют вилку, внутри которой лежит точное решение задачи теории упругости.

°ху > Па

Рис. 5. Графики распределения компонент тензора напряжений: а, б, в - при у = Ы2 ; г, д, е — при х = а/ 2;-•— метод геометрического погружения

в напряжениях, * метод конечных элементов в перемещениях

Рис. 6. Сходимость итерационной процедуры с увеличением числа итераций: а - на примере ст^, Па на границе у = 0, x е[4,8]; б - на примере <зуу, Па

на границе у = 0, х е [4,8]; —*— решение, полученное методом конечных

элементов в перемещениях; *— решение, полученное методом

геометрического погружения на каждой итерации

Таким образом, продемонстрировано практическое применение метода геометрического погружения на основе вариационного принципа Кастильяно и его конечно-элементной реализации в напряжениях. Получено достаточно хорошее соответствие результатов определения полей напряжений в сравнении с традиционным МКЭ в перемещениях. Установлена практическая сходимость итерационной процедуры. Возможно развитие предлагаемого подхода на области с более сложной, в том числе криволинейной, геометрией.

Библиографический список

1. Галлагер Р. Метод конечных элементов. Основы. - М.: Мир, 1984. - 428 с.

2. Шардаков И.Н., Труфанов H.A., Матвеенко В.П. Метод геометрического погружения в теории упругости / УрО РАН. - Екатеринбург, 1999. - 298 с.

3. Каменских A.A., Труфанов H.A., Матвеенко В.П. Численная реализация метода геометрического погружения на основе вариационного принципа Кастильяно // Вестник ПГТУ. Механика. - Пермь: Изд-во Перм. гос. техн. ун-та, 2010. - № 3. - С. 5-18.

4. Суходолова Ю.С., Труфанов Н.А. О конечном элементе на основе вариационного принципа Кастильяно для плоских задач теории упругости // Вестник ПНИПУ. Механика. - Пермь: Изд-во Перм. гос. техн. ун-та, 2012. - № 1. - С. 168-178.

5. Nesrin Sarigul, Richard H. Gallagher Assumed stress function finite element method: two-dimensional elasticity // International journal for numerical methods in engineering - 1989. - Vol. 18. - P. 1577-1598.

6. Girija C.V. Vallabhan, Azene Muluneh A finite element model for plane elasticity problems using the complementary energy theorem // International journal for numerical methods in engineering. - 1982. - Vol. 18. -P. 291-309.

7. Тимошенко С.П., Гудьер Дж. Теория упругости. - М.: Наука, 1975. - 576 с.

References

1. Gallager R. Metod konechnykh elementov. Osnovy [Finite element analysis fundamentals]. Moscow: Mir, 1984. 428 p.

2. Shardakov I.N., Trufanov N.A., Matveenko V.P. Metod geomet-richeskogo pogruzheniya v teorii uprugosti [Method of geometrical immersion in the theory of elasticity]. Ekaterinburg: Ural’skoe otdelenie Rossiy-skoy Akademii Nauk, 1999. 298 p.

3. Kamenskikh A.A., Trufanov N.A., Matveenko V.P. Chislennaya realizatsiya metoda geometricheskogo pogruzheniya na osnove variatsion-nogo printsipa Kastilyano [Numerical realization of the geometrical immersion based on Castigliano variational principle]. Vestnik permskogo gosu-darstvennogo tekhnicheskogo universiteta. Mekhanika, 2010, no. 3, pp. 5-18.

4. Sukhodolova Y.S., Trufanov N.A. O konechnom elemente na os-nove variatsionnogo printsipa Kastilyano dlya ploskix zadach teorii uprugosti [About a finite element based on the castigliano variational principle for plane elasticity problems]. Vestnik permskogo gosudarstvennogo tekhnicheskogo universiteta. Mekhanika, 2012, no. 1, pp. 168-178.

5. Nesrin Sarigul, Richard H. Gallagher Assumed stress function finite element method: two-dimensional elasticity. International journal for numerical methods in engineering, 1989, vol. 18, pp. 1577-1598.

6. Girija Vallabhan C.V., Azene Muluneh. A finite element model for plane elasticity problems using the complementary energy theorem. International journal for numerical methods in engineering, 1982, vol. 18, pp. 291309.

7. Timoshenko S.P., Goodier J.N. Teoriya uprugosti [Theory of elasticity] Moscow: Nauka, 1975. 576 p.

Об авторах

Труфанов Николай Александрович (Пермь, Россия) - доктор технических наук, профессор, заведующий кафедрой вычислительной математики и механики Пермского национального исследовательского политехнического университета (614990, г. Пермь, Комсомольский пр., 29, e-mail: vmm@cpl.pstu.ac.ru).

Кузнецова Юлия Сергеевна (Пермь, Россия) - аспирант кафедры вычислительной математики и механики Пермского национального исследовательского политехнического университета (614990, г. Пермь, Комсомольский пр., 29, e-mail: Suhodolchik@mail.ru).

About the authors

Trufanov Nikolay Aleksandrovich (Perm, Russian Federation) -Doctor of Technical Sciences, Professor, Head of Department of Computational Mathematics and Mechanics, Perm National Research Polytechnic University (29, Komsomolsky av., 614990, Perm, Russian Federation, e-mail: vmm@cpl.pstu.ac.ru).

Kuznetsova Yuliya Sergeevna (Perm, Russian Federation) - Postgraduate student of Department of Computational Mathematics and Mechanics, Perm National Research Polytechnic University (29, Komsomolsky av., 614990, Perm, Russian Federation, e-mail: Suhodolchik@mail.ru).

Получено 7.02.2013