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

2012 Механика №1

УДК 539.3

Ю.С. Суходолова, Н.А. Труфанов

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

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

Рассмотрен вариант метода конечных элементов для плоских задач теории упругости, основанного на вариационном принципе Кастильяно. Приведены количественные характеристики сходимости метода на примере решения задачи о растяжении пластины переменной нагрузкой. Даны рекомендации по реализации статических граничных условий в плоских задачах теории упругости.

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

Y.S. Sukhodolova, N.A. Trufanov

State National Research Polytechnic University of Perm, Perm, Russian Federation

ABOUT A FINITE ELEMENT BASED ON THE CASTIGLIANO VARIATIONAL PRINCIPLE FOR PLANE ELASTICITY PROBLEMS

We investigate the variant of the finite element method for plane elasticity problems, based on Castigliano variational principle. The quantitative characteristics of the convergence of the method on the example of solving the problem of tensile plate of variable load are resulted. We make recommendations for the implementation of the static boundary conditions in plane elasticity problems .

Keyword: finite element method, Castigliano variational principle, static boundary conditions.

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

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

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

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

Функционал дополнительной энергии упругой системы имеет

вид

У

(о,ь>

Рис. 1. Расчетная схема поставленной задачи

1 а Ь

Пс (0) = 2Е 11 + 0 - 2уо х0 у + 2 (1 + у) т2у ] йхёУ.

2Е

° °

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

их (О, у )= иу (х,°) = О, Го хх (а у ) = р (у),

■ (a, у) = Тху (x, Ь) = тху (°, у) = Тху (х,°) = 0уу ( Ь) = °.

(1)

(2)

(3)

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

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

д 2ф ду

2 ’ 0 уу

д2 ф дх7’

Т ху =-

д2 ф дхду

(4)

Тогда (1) примет вид

1 а Ь

П(0 )= 2Е

д2ф^ +^д 2ф^ ,ду2J +1дх' у

- 2у

(д2ф^( д2ф^

ду2

+ 2 (1 + у)

22

д ф дхду

(5)

бхёу.

Для дискретизации функционала (5) выберем прямоугольный конечный элемент с четырьмя узлами в углах (рис. 2). Поскольку по форме задача для функции напряжений (бигармоническое уравнение) соответствует задаче об отыскании прогиба прямоугольной пластины, а возможные типы граничных условий в обеих задачах одинаковы (условия на искомые функции и их первые и вторые производные по координатам), разумно исследовать возможности конечных элементов, апробированных для задач изгиба пластин. Известны рекомендации [1, 3] по применению прямоугольного конечного элемента с тремя неизвест-

дф^ (дф^

ными

=*'-(5 )•

ду

■в каждом узле. Однако в таком случае в

угловых точках элемента возникает противоречие, так как в них нельзя

д2ф = д2ф

обеспечить выполнение равенства

дхду дхду

чениях узловых неизвестных элемента [8].

при произвольных зна-

Стхх (0,0)Па

(0,1) (1,1)

і к

<! І І с —

(0,0)

(1,0) Тху (0,0), Па

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

Для преодоления этого затруднения примем смешанную производную в качестве одного из узловых параметров:

дф^ (дф^ ( д2ф ^

{«' =М2).•

І (“У У

йу I ’(дхду

Используя локальную систему координат, функцию напряжения Эри в элементе можно аппроксимировать выражением [8]

4 (

- дфІ — дфІ — д2ф

- + N

+ N

{ (С,п)}Г {5'}, (6)

дх ду дхду у

где ( п)} - столбец функций форм элемента;

{бе} = {5', 5 ^, 5 к, 5г} - столбец узловых неизвестных элемента;

N„ = N°. (С) N°. (п), N2, = N°. (С) N (п);

Мз, = N (СШп), N4, = N (С) N., (п).

Функции формы Nij в уравнениях (6) определяются как

N„1 (?) = N°. (?) = 1 -3С- + 2?3, N„1 (п) = N„4 (п) = 1 -3п2 + 2п3, Nll (?) = N,4 (0 = а (-2?2 +С3), Nll (п) = N„(4) = Ь (п-2п2 + п3),

N,2 (?) = N°з (?) = 3?2 - 2?3, N°. (п) = N,3 (п) = 3п2 - 2п3,

N,2 (?) = N,3 (О = - а (2 - С3), N,4 (п) = м„ (п) = -Ь (п2 - П ) .

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

Используя связь функции напряжения Эри (4), введем следующее обозначение:

[ * ]

_д!

ду2

д2

{'}'

дх2

д2

дудх

{¿V'}'

_± _д_

Ь2 дп2

4 д2

а2 дс

4 д2

аЬ д^дп

К }

{/У'}' {Я'}'

где Ье ] - матрица градиентов конечного элемента для функции напряжений в случае плоской задачи. Таким образом, функционал (5) для конечного элемента примет вид

П ' =

211 {5'}' [ *' ]' [ °" ][ *' ]{5' п = 0,

где Э 1 ^ - матрица констант упругой податливости:

1 - V2 -V (1 + V) °

-V (1 + V) 1 - V2 °

° ° 1 + V

в случае плоско-деформирован-

ного состояния;

1 -V

0

-V 1 0

0 0 1 + V

в случае плоско-напряженного состояния.

Согласно вариационному принципу Кастильяно требование ста-

ционарности функционала дополнительной энергии

дПс

д{5}

= 0, где

П с = IП/ , приведет нас к системе линейных алгебраических уравне-

е

ний [С]{8} = {Ж}, где {8} - глобальный столбец узловых неизвестных.

Основные трудности, возникающие при реализации данного варианта метода конечных элементов, заключаются в задании граничных условий в напряжениях. Граничные условия в перемещениях для данной вариационной формулировки являются естественными, и их учет труда не составляет. Выполнение главных статических граничных условий (3) может представлять собой достаточно сложную процедуру. Рассмотрим учет граничных условий на примере одной стороны элемента с узлами ' и -, параллельной оси Оу. В любой точке на этой д'ф Д

стороне а^ = —-. Дважды интегрируя это выражение и определяя

ду

константы интегрирования через ф и фу в узловых точках, получим

уу

-фу +ф я = ^

у'

уУ

-фг + ф- - фУ'Ь = Ца„.

у' у'

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

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

рировать по внешним границам элементов заданные условия в напря-

У У

жениях (вычислить интегралы вида: / (у) = |а^йу, g (у) = |/ (у) йу.).

0 0

Для сравнительного анализа результатов поставленная задача решена тремя способами:

1) МКЭ, основанным на вариационном принципе Кастильяно;

2) МКЭ, основанным на вариационном принципе Лагранжа (четырехугольный билинейный элемент);

3) путем прямой минимизации функционала дополнительной энергии упругого тела методом Ритца [9].

Для получения приближенного аналитического решения методом Ритца была выбрана аппроксимация решения на основе базисной функции вида

Ф = 2 £у 2 1 - 6 Ъ^] + ( - а 2 )2 (у 2 - ъ2 )2 (а1 + а2х 2 + аз у 2) ,

где а1, а2, а3 - неизвестные постоянные, £ = 20.

Подставляя все в функционал (5) и дифференцируя по неизвестным постоянным коэффициентам, придем к линейным уравнениям относительно а1, а2, а3, получим следующее их значение:

V V V

а, = 0,07983-—, а2 = 0,1250-—, а3 = 0,0182^-——.

1 а Ъ2 2 а Ъ 3 а Ъ

Приближенное решение задачи представлено на рис. 3. На рис. 4 представлены результаты, полученные методом конечных элементов в напряжениях на сетке 10 х10 узлов. Для сравнительного анализа метода конечных элементов в формулировках Лагранжа и Кастильяно представлено изменение напряжений в конкретных точках пластины в зависимости от числа узловых неизвестных в конечно-элементных сетках (рис. 5). МКЭ в формулировке Кастильяно при определении напряжений имеет значительно более высокую скорость сходимости.

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

а б в

Рис. 3. Приближенное аналитическое решение задачи. Изолинии равных значений напряжений: а - а^, Па , б - ауу, Па , в - тху, Па

а б в

Рис. 4. Решение методом конечных элементов в напряжениях. Изолинии равных значений напряжений: а - ахх, Па , б - а , Па , в - т , Па

ауу (°’°)’Па

аXX ( 0,0) > Па

Тху ( °,0), Па

а N б N 6

Рис. 5. Изменение напряжений с увеличением числа степеней свободы в конечноэлементном аналоге: 1 - решение на основе вариационного принципа Кастильяно,

2 - решение на основе вариационного принципа Лагранжа; а - ахх (0,0), Па ,

б - ауу ( ^ ,Па , в - Т*у ( 0^) , Па

Для формулировки в напряжениях эти ограничения отсутствуют. Рассмотрим результаты применения МКЭ в напряжениях к решению задачи о поперечном сжатии переменной распределенной нагрузкой изотропного бруса с коэффициентом Пуассона, равным 0,5: левая грань бруса закреплена жестко, на правой грани приложена параболи-

N

ческая сжимающая нагрузка, верхняя и нижняя грани свободны от нагрузок. На рис.6 представлены графики сходимости решения по нормальным напряжениям: изображены напряжения в центральной точке пластины в зависимости от числа узловых неизвестных. Решение практически сходится при 1600 степеней свободы, что соответствует конечно-элементной сетке 20 х 20 узлов.

ст»т (0,0),Па

ал,(0,0),Па

а N б N

Рис. 6. Изменение напряжений с увеличением числа степеней свободы. Решение методом конечных элементов в напряжениях: а - ст^ (0,0), Па , б - ст^ (0,0), Па

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

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

1. Зенкевич О.К. Метод конечных элементов в технике. - М.: Мир, 1976. - 541 с.

2. Сегерлинд Л. Дж. Применение метода конечных элементов. -М.: Мир, 1979. - 392 с.

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

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

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

6. Демидов С.П. Теория упругости: учебник для вузов. - М.: Высш. школа, 1979. - 432 с.

7. Победря Б.Е. Численные методы в механике деформируемого твердого тела. - М: Изд-во МГУ, 1995. - 368 с.

8. Д. Норри, Ж. де Фриз Введение в метод конечных элементов. -М.: Мир, 1981. - С. 205-207.

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

References

1. Zenkevich O. The finite element method in engineering science [Metod konechnykh elementov v tekhnike]: Transl. from eng. Moscow: Mir, 1976, 541 p.

2. Segerlind L.Dzh. Applied finite element analysis [Primenenie me-toda konechnykh elementov], Moscow: Mir, 1979, 392 p.

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

4. Shardakov I.N., Trufanov N.A., Matveenko V.P. Method of geometrical immersion in the theory of elasticity [Metod geometricheskogo pogruzheniya v teorii uprugosti], Ekaterinburg: UrO RAN, 1999, 298 p.

5. Kamenskikh A.A., Trufanov N.A., Matveenko V.P. Numerical realization of the geometrical immersion based on Castigliano variational principle [Chislennaya realizatsiya metoda geometricheskogo pogruzheniya na osnove variatsionnogo printsipa Kastilyano]. PSTUMechanics Bulletin -VestnikPGTU. Mekhanika, Perm, 2010, No. 3, pp. 5-18.

6. Demidov S.P. Theory of elasticity [Teoriya uprugosti], Moscow: Vyssh. shkola, 1979, 432 c.

7. Pobedrya B.E. Numerical Methods in the Theory of Elasticity and Plasticity [Chislennye metody v teorii uprugosti i plastichnosti]. Moscow: MGU Publ., 1995, 368 p.

8. Norri D., Zh. de Friz An Introduction to Finite Element Analysis [Vve-denie v metodkonechnykh elementov]. Moscow: Mir, 1981, pp. 205-207.

9. Timoshenko S.P., Guder Dzh. Theory of Elasticity [Teoriya uprugosti]: Transl. from eng. Under. ed. Shapiro G.S. Moscow: Nauka, 1979, pp. 269-272.

Об авторах

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

Труфанов Николай Александрович (Пермь, Россия) - доктор технических наук, профессор, заведующий кафедрой вычислительной математики и механики Пермского национального исследовательского политехнического университета (614990, г. Пермь, Комсомольский пр., 29, е-mail: [email protected]).

About the authors

Sukhodolova Yuliya Sergeevna (Perm, Russian Federation) - Student of Department of Computational Mathematics and Mechanics, State National Research Polytechnic University of Perm (614990, 29, Komsomol-sky prospect, Perm, Russian Federation, е-mail: [email protected]).

Trufanov Nikolay Aleksandrovich (Perm, Russian Federation) -Doctor of Technical Sciences, Professor, Head of Department of Computational Mathematics and Mechanics, State National Research Polytechnic University of Perm (614990, 29, Komsomolsky prospect, Perm, Russian Federation, e-mail: [email protected]).

Получено 10.02.2012