Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

книги / Теория упругости. Ч. 1 Основные уравнения механики сплошных сред. Вариационные методы

.pdf
Скачиваний:
2
Добавлен:
12.11.2023
Размер:
5.05 Mб
Скачать

 

 

 

 

 

 

- IOI -

 

Здесь использованы

табличные

интегралы

 

/'cos.аха!х= ^

si*a x ;

Захс(х = ^ r v ^

л > г <?х,

 

 

 

 

(cosА х)п = A ” os (t x *

, (б . 16)

После вычисления выражение потенциальной энергии примет вид

 

П~

f

с'т

 

£

< « • " )

 

 

' е7' а * Ш

§ '

Для определения параметра

CL

вычисляем

 

§ " £J*a{0 !

 

откуда

 

 

 

а

 

££_/£_),

(6Л9)

 

 

 

 

 

 

 

 

 

 

 

л * ( л /

 

С учетом этого уравнение упругой линии балки принимает вид

 

 

 

 

■32’ /

^

п/>4

 

 

 

 

=л *

 

('-*»§)•

(6-го)'

л *

//■

 

 

 

 

Iл=е

I

- £ )ё(':созтУ°-”957щ-- (6-г1)

Сравнивая (6.21) с (6*13) устанавливаем! что оаибка составляет

4,5% от точного решения.

 

 

 

 

Определим величину m a x

напряжения:

 

e

” - _ I ^ \ I Z a U

t o s * ± ;

Г 6 .2 2 )

 

К

Ч

К

( 2 £ /

2 е

 

 

 

 

 

 

 

 

.25)

 

 

=-^т / / - — ) - 2 4 - = 0 .2 9 4 5 4 - ? 4 • (6.;

 

 

я * y

' x J

W.

К

 

- 102 -

Точное решение сопромата дает

^тах

. * £

- 0 5 - 3 ?

(6.24)

2W,

и' ° W,

 

Величина погрешности kl%. Следовательно, метод Ритца при опреде­ лении перемещения дает в рассматриваемом случае лучшие результа­ ты, чем при определении напряжения. Бели функцию V взять в ви­ де ряда

гг--

 

(6.25)

где Л =1,3,5,?,....,

то при к ~5 погрешность в

дос­

тигнет 0,03$, а в

8,1$6.

 

§6.3. Метод Бубнова - Гадеокина

Вряде случаев оказывается более предпочтительным не составлять

выражения потенциальной анергии системы, как в методе Ритца, а решать приближенно дифференциальные уравнения, например, урав­ нения равновесия.

Одним на наиболее распространенных методов приближенного ранения дифференциальных уравнений прикладной теория упругости является метод Бубнова-Галеркина. Положим, что решение задачи сво­

дится к дифференциальному (интегральному) уравнению

 

 

L - (х, иг, их', иг;...)-0.

(6.26)

Примем, что

 

 

u r ( a 7,а 2,а3 f. . х);

(6.27)

где

а п - произвольные параметры, жак и в методе Ритца. При

■вбых значениях этих параметров функция tt/fx) должна удовлетво­ рять геометрическим ,(накладываются на линейные и* угловые переме-

«ения, т.е. Щ ы Л

) н силовым (определяются моментами и перере­

зывающими силами,

т.е. иг*-иг")

граничным условиям. При реие-

яяи задач методом

Ритца бывает

доотаточно удовлетворить толь­

ко геометрическим

условии. В других методах, как это полазаю

- ю з - во нногпз: работах, необходимо. удовлетворить воем граничным уело-

в Ъ ш, как геометрическим, гак в силовым (метод Галеркина, Графф-

ца

и др). При выборе функции (6,27) следует руководствоваться -

в

первую

очередь видок ожидаемого решения* Как и в методе 'Рвтца,

фупздяя

иг должна быть как можно ближе к, точному решению при

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

По методу Бубнова - Галеркина функция (6.27) берется в виде

ряда

 

“ '•=‘2, Ъ (х) * аг Ъ ( х) '1, ■■'■+a » f n ( x) ,

(6.26)

удовлетворяющего всем граничнЙи условиям задачи. Далее и г по

(6„28) подставляем в уравнение (6*26).

 

2 ,(x ,a „ a e , . . ' . r a J * / C x ) .

(6.29)

Затек фупкцию-ошибяуj ? ( x ) уннояаеи поочередно ъъ.‘/ п %,■■■,

берем от этих произведений интегралы по всей области изменения

Л. о приравниваем каждый из нн£ нулю и получаем т а кш

'образом

езегецу уравнений по .числу, неизвестных а п :

 

6

 

/ у ( х ) у п ф - 0 .

(6.30)

а

 

Определяя на этой системы постоянные. (2П , подставляем их в прмнятуа функцпа (6.28) в получаем прнбляжённоё решение уравнения

(6о26)*

Так выглядит чисто внеяняя сторона метода Б.Г.Галаркнна.

Этот

метод был применен им в 1915 г.*,Б.Г .Галеротным было пока­

зано,

что путем только что отложенного приема можно подучить

весьма

точное решение многих задач 'строительной механик. В

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

решении многих задач прикладной теории упругости н механики. Еще раньше Галеркина в отзыве на работу Тямояенко К*Г«Буб­

нов. указал, что ранение расожатр)зваеиых задач можно получить,

не прибегая к понято»

потенциальной энергии системы* Он рекомен­

довал подставлять ряд

(6*28) в дифференциальное уравнение-, а

затем поступать опноавннм выше обрааом,-

 

Напоиннр, что семейство функций

(л ) , обладающее свой­

ством

 

- 104 -

6

 

 

< > * / )

(6.31)

a s

>

 

 

£ - * i

)

 

a

 

 

называется семейством ортогональных в

промежутке

<2 - 6 функций,

а -свойство (6.31) - свойством ортогональности. Этим свойством обладают, в частности, тригонометрические функции, некоторые по­ линома и др.

С учетом (6.31) условия Бубнова - Галеркина (6.30) модно рассматривать как.требование, чтобы функция - ошибка была орто­

гональна к функциям

. Если решение

найдено точно, то/*(х )Щ

следовательно, условия

(6.30) выполняются, так как функция,

тож­

дественно равная на интервале а - 8 нулю,

ортогональна вообще

но

всем функциям.» Однако

решая задачу приближенно, мы не можем

 

требовать, чтобы J> (x )

была ортогональна

ко всем функциями до­

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

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

Методы Ритца и Галеркина обладают различной чувствитель­ ностью к апроксимацин степени производных координатной функции.

В методе Ритца, например, необязательно

выполнять граничные усло­

вия в высших производных,-в то время как в

методе Галеркина за

этим необходимо внимательно

одедить; иначе

мы

рискуем подучить

большие погреинооти. Иногда

выбор функции

 

v f

в наиболее естест­

венном виде приводит к совершенно неверным результатам.

для иллюстрации метода Бубнова - Галеркина рассмотрим тот же пример с консольной балкой, решенный методом Ритца. Диффереи-

 

 

- 105 -

 

циальнсо уравнение упругой

линии балки будет

т е т ь вид

 

 

 

 

E 7 x l r "

- U ( e - s ) ‘ : 0 .

(6.32)

ЕСДД; HSE и прежде, принять 8а функцию

 

 

гг.а (г-СИg )

(6.33)

а подставить (6,33) в (6.32), то в результате преобразований

<6,28)-(6.30) полушш

 

 

 

/ [EJ*a

]('-со* й У ^ 0,

откуда

 

(6.34)

 

1__<S ^ / 6

4

 

 

 

(6.35)

 

4

 

 

что дает,более чей двухкратную ошибку.

 

 

Получанаоо раехоздзнве объясняется тем, что выбранная функ­

ция

V хорошо ограаает уравнение упругой линия балки и плохо -

его

производные (первую и особенно вторую). При репенни задачи

кетодон Рнтца вследствие зтого получаются расхождения лишь при определении производных (напряжений)* При ревенин же задачи ис­ тодом Бубнова - Галеркина большая ошибка ' имеет место не толь­ ко в производных, но а з самой функции. Это видно ив рис. 6.2, на которой рядом с истинными эпюрами изгибающих моментов (тг")

и перерезывающих cis ( v

w) пунктиром показаны

эпюры,

соответствую­

щие принятой функции (6.33). Разница, как видно, разительная*

Кроне того, функция V

не удовлетворяет всей граничный условиям

коисольвой балки (а именно: силовому условию

^

О ) , что

явхиетоя основной причиной расхождения. Следовательно, функция.

v должна быть выбрана

иначе*

' При решении задачи неходок Бубнова - Галеркина подбор ап-

рокснмирующей функции

практически удобнее начинать ео старшей

производной, входящей

в

дифференциальное уравнение (6.32). При-

- 106 -

ir

(6.-36)

Эта функция хорошо отрапазт

ВТОРУЮ ПРОИЗВОДНУЮ. П00Д9 1Щ-

тегрирования подучим

Из условий v - 0

j

, / '’ 0

1 - 0

инеем

 

 

 

/1=-—

 

;

 

 

Подставляя

V в уравнения

(6.32)

и производя

интегрирование

(6.30) по

£

от 0 до

f

, получим

 

 

 

 

a _ j$ £ _

 

60 *я * -rs 6ж '

 

 

 

~ 2 fJ r

'

6 * Я 3

ж г я

 

 

 

 

 

 

 

 

Откуда

 

 

 

 

 

 

 

 

 

 

-4

‘1§ Р / Н ' “ ‘ё

46.39).

 

 

 

 

 

 

WK

тах

'

Wr

 

 

Точность,

как видим, в

обоих случаях вполне

удовлетворительная

Дальнейшего уточнения можно было бы достичь путем принятия внес то V ' ряда

(6Л1)

- 107 -

§ 6.Л-» Салланвый метод

Этот иетод был применен в 1933'году Л.В.Канторовичем к кручению стержня прямоугольного профиля, В.Дунканом - к кру­

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

Идея смешанного метода заключается В'том, что искомые функ­ ции представляют в виде произведения двух функций. Однаиз функ­ ций - известная, подобрана так, чтобы частично удовлетворить "раничныи условиям. Другая не - неизвестная, должна зависеть от йеиьшего числа переменных и определяться при пойощи вариационно­ го качала или из решения уравнения Эйлера-Лагранжа. Таким обра­ зом, в случае решения задачи, зависящей от двух переменных! при приведении к обыкновенным дифференциальным уравнениям решение разыскивают в такой форме, когда в состав его входят неопреде­

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

Для выяснения сущности метода рассмотрнм пример: найти приб­ лиженное решение уравнения Пуассона (основное уравнение кручения

стержней)

= - /

для прямоугольника/-^; - S ; 6 ]

при условии

■и * О на

контуре. Будем решение

искать в форме

 

 

 

й = ( б 2~ у 2) / ( х

) \

(6.Д2)

оно удовлетворяет

граничному условию на прямых

. И в

матаналиэа известно, что задача Дирихле для уравнения Пуассона свяаана с задачей о минимуме интеграла Дирихле, соответствующего данной задаче

- 108 -

где интеграл

распространен

на плоскою, область

S , ограничен -

ную контуром

Г • Поставим теперь задачу найти

функцию й (х , у ) у

непрерывную в области S

вмеоте с частными производпши первого

и второго порядка, принимающую на контуре Г .заданные значсшш

а - ср( 5)

н дающую интегралу

(6.43) минимальное

значение.

В этом случае, составив' уравнение Эйлера,

убедимся, что

оно яв­

ляется уравнением Пуаооона.

 

 

 

 

 

 

 

Подставляя значение и

в

интеграл

(6.43),

имеем

 

J ( u ) ~ / d K ^ [

 

 

 

 

^ y ~ 2 f ( 6 2- y s) ] ^ .

(6.«)

Выполняя интегрирование па у

 

, приводим задачу

к минимуму прос­

того интеграла

 

 

 

 

 

 

 

 

 

 

Ф И

Ш

?

Ф

У

- т

6У У *

>

<«•«>

т.е. получили интеграл вида

 

 

 

 

 

 

 

 

 

 

X,

 

 

 

 

 

 

 

 

j -

f

f ( x ,y

, у ')Ы

х ,

.

х -О .

 

 

хо

Для такого интеграла уравнение Эйлера извеотно и записывается в виде

d

(6.4ё)

a r t - f я*

 

Соответствующее уравнение Эйлера для интеграла (6.45) будет следующим:

£ Ш 6У')-т6> +1 6^ 0’

f

? Ё г J

.

(6.47)

 

 

4 6 *

 

 

 

- 109

-

 

 

Роная $то уравнение с постоянный! коэффициентами,

находим

 

S ( * ) * c , . i A $ £ + b & J § f + £ -

(6 . « ,

 

Постоянные сг и

сг

определяются из условийf ( a . J = ^ (-(2)~

о,

обеспечивающих выполнение условия и ~ 0

на сторонахX ~ + Q # д0_

скольку решение,

очевидно, будет четной функцией х

, 10 с

= 0 -

для определения

cf

инеем

 

 

/

С,С/гУ §

Т * 2 ~ ° '

 

 

где

 

 

 

 

s A x -

е*-е~ х

 

 

2

сА х = -

 

 

Теперь подставляя найденное значение

с г я (6.48),

подучаем

/ ( * ) =

 

 

(6.49)

 

 

 

 

(6.50)

Отсюда, например, для крутящего моневта инеем

 

 

 

а /

c /iiM * - )

 

М = ш [ / Q ( x , y ) d x d y = 2 Ё в Г (/---- Ц ^ - Ы х Р ( б 3-

s

 

- м

ci H /

-*

 

 

 

 

(6.51)

Сравним это значение М (таблица) с точный (подучено при реше­ нии дифференциального уравнения Пуассона методой Фурье) и приб­ лиженным значением по методу Ритца (первое приближение)!

 

 

- н о

-

Методы

#•'

-

п

 

6

Z

Точный

2,244

3,659

Смешанный

2,234

3,653

Метод Ритца

2,222

3,555

и

4,213

4,209

4,000

а

в * * 0

5,335

5,333

4,444

Из таблицы видно, что наибольшая ошибка - около 0,5$ полу­

чается

для квадратного сечония; если

, совпадение с

точным

решением

црчти полное. Метод Ритца в

случае ^ = 1 да­

ет ошибку вдвое

большую, чем смешанный.

 

§ 6.5. Метод минимума квадратичного уклонения

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

и = у-fz , а

,, а 2 , . . Q „ ) f.

(6.52)

рассматриваем выражение функции

- ошибки

 

а

an ) = f f e ) ± o .

(G.53)

Постоянные <з,, ^ ^ / " ^ п о д б и р а ю т так, чтобы интеграл от квад­

рата /*(& ) имел бы минимальное значение. Это позволяет сделать функцию - ошибку наименее отклоняющейся от нуля в среднем. Если бы мы не возвели y f e ) в квадрат, то операция отыскания миниму­

ма от интеграла \ f\ ffa ) не имела бы смысла. Этот интеграл

мог бы быть равен нулю при функции - ошибке, сколь угодно отли­ чающейся от нуля. Квадрат леf f e ) имеет постоянно положитель­ ный знак. С таким же успехом можно было бы взять четвертую или любую другую четную степень J ? (g ) , но это привело бы к более сложный выкладкам.