Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:
Численные методы решения задач строительства на ЭВМ..pdf
Скачиваний:
45
Добавлен:
15.11.2022
Размер:
16.27 Mб
Скачать

(6.83)

dx j

\d y

6.5.3. Метод Ритца

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

\F {x ,y ,y 'У ,..)dx

(6.84)

А'о

Идея метода Ритца заключается в том, что искомую функцию у(х) доставляющую экстремум функционалу,

представляют в виде линейной комбинации функций щ(х):

п

=

=

+ а2Ф2(*) + - + яиФнМ -

(6.85)

 

/=1

 

 

Число п зависит от требуемой точности.

Здесь срi(x) - координатные функции, выбираемые вполне определенным образом, а именно: так, чтобы функция у(х) удовлетворяла граничным условиям задачи (или хотя бы части из них). Точность решения в большой степени зависит от удачного подбора этих функций и, вообще говоря, возрастает с увеличением их числа. Подбор этих функций требует определенного навыка. Для линейных задач чаще всего применяют полиномы или

тригонометрические функции.

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

Подставляем решение (6.85) в функционал (6.84). После интегрирования и подстановки пределов он становится функцией, зависящей от параметров aL:

а< ) =

>а2 >•••> ) •

(6.86)

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

¥_ = 0; f = 0;

= 0.

(6.87)

да\

да„

 

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

Рассмотрим применение метода Ритца на примере изгиба стержня постоянного сечения под действием равномерной нагрузки q (рис.6.18).

Пример 6.13. Запишем выражение полной потенциальной

 

энергии

изгибаемой

балки

q

(функционал):

 

* * * * * М

П = j[± £/x( /

)2 - qxy]dx

(6.88)

-И-------------------[

О

 

 

 

В соответствии с равенством (6.72) принимаем

Рис.6.18. К примеру 6.13

прогиб балки

 

 

 

П

 

 

 

у{х) =

.

 

 

/=1

 

(6.89) Координатные функции %(х) в рассматриваемом примере должны удовлетворять всем или части граничных условий, например, ^до) = 0 и

^/(0) = 0(прогиб и угол поворота в жесткой заделке равны 0). Можно принять:

 

V/+I

 

^ w = fy ] ; р2М = (уЛ

= | у

(6.90)

Решим сначала задачу с одним параметром, положив

 

У ( х ) = я,

X\ 2

(6.91)

 

Тогда,

у" (л) =

.

Внося эти выражения в функционал (6.88),

получаем

 

 

 

 

 

 

 

 

 

 

 

 

\ 2 .

2EJX 2

qr

 

л-^ 'я М Ч т

dx =г— а , ---- г-а,.

 

 

 

13

 

З/2

 

Из условия (6.87)

дП

4EJX

ql

Л

 

 

да,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

находим

 

«1 =

*/4

 

 

 

 

(6.92)

И уравнение прогиба балки

\2Е1Г

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

= 12EJX

 

 

 

(6.93)

 

 

 

 

 

 

 

Максимальный прогиб балки на свободном конце при х=1

 

 

 

 

ф 4

 

 

 

 

(6.94)

 

 

Л1) = 12EJr

 

 

 

 

 

 

 

 

 

что отличается от точного значения

 

 

 

 

 

 

 

yQ) =

¥ 4

 

 

 

(6.95)

примерно на 17%.

 

 

8EJS

 

 

 

 

 

 

 

 

 

 

 

 

Для уточнения решения примем

 

 

 

 

 

у(х) = а

у !

+а.

х ^

 

(6.96)

 

 

 

 

 

 

 

\ h

 

 

Подставляем это

выражение

 

и

выражение 2-й

производной

."'-л _

о. ^xaL

в функционал (6.84):

 

 

у ( х ) = 1 >

+>

о 1_+ 6 ^ 1\2

 

 

 

 

 

я = ^ Н /2-

 

 

 

+ a2(y j

 

2 Д /2

/3

^-<7

 

 

 

 

 

 

 

 

Дифференцируя П по aj и а2получаем систему уравнений

дП

дах

дП = EJ

 

 

 

 

dx - О

да2

 

 

 

 

 

или

 

3

 

q f

 

ci\ + — а7

= -------

 

 

1 2

2

\2EJ

 

 

ax+ 2a2 = Я14

 

 

 

 

 

2AEJX J

Отсюда находим:

«I

5

gl4

 

24

EJX

 

 

 

 

 

 

 

 

/4

 

 

 

 

qr

 

 

«2 =• 12EJr

 

Подставляя а, и a2в уравнение (6.96), получаем

 

 

 

5

gl2x2

glx*

 

У(х) = 24

EJr

12EJX

(6.97)

(6.98)

Наибольшее

значение прогиба у{1) =

qi4

что совпадает

с точным

 

 

8£/„

 

 

значением

(6.95). Таким образом,

добавление одного

параметра

существенно повысило точность решения.

 

 

6.6. Решение дифференциальных уравнений с

использованием электронных таблиц Microsoft Excel

6.6.1. Решение задачи Коши методом Эйлера

Решить задачу Коши у =2ху, у(0) = 1

(6.99)

методом Эйлера для двух разностных сеток с шагом й=0,2 и h=0,1.

Сравнить полученные результаты с известным точным решением

у =е Проанализировать полученные приближенные решения,

сравнить их с имеющимся точным решением и сделать вывод о продолжении или прекращении итерационного процесса.

Расчетная схема для h=0,1 строится аналогично. Полученное приближенное решение имеет вид

 

 

 

 

 

 

 

 

 

Таблица 6.2.

Xi

0

0,1

0,2

0,3

0,4

0,5

0,6

0,7

0,8

0,9

1

Yi(li/2)

1,000 1,008 1,037 1,089 1,166 1,273 1,419 1,613 1,871 2,214 2,673

Сравнение полученных приближенных решений приведено на рис.6.23 и разница между ними не превосходит 5%.

Метод FyHre-Кутта ( h=0.2)

Рис.6.23.Сравнение результатов расчета

Таким образом, в качестве решения исходной задачи Коши принимаем сеточную функцию таблицы 6.2.

6.6.2. Решение краевой задачи методом конечных разностей

Методом конечных разностей найти решение краевой задачи (6.62) (пример 6.11) на отрезке хе[1, 2] с шагом //=0,25 и с шагом А=0,125. Сравнить полученные результаты и сделать вывод о необходимости поиска следующего приближения или о прекращении счета.