Добавил:
Upload Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:
Математика_Семестр3_МетодПособие.pdf
Скачиваний:
157
Добавлен:
05.06.2015
Размер:
1.75 Mб
Скачать

80

В качестве примера рассмотрим решение вышеизложенными методами задачи отыскания прогибов шарнирно опёртой балки постоянного поперечного сечения, нагруженной равномерно распределённой поперечной нагрузкой интенсивности q (см. рисунок 2).

Рис. 2

 

Задача сводится к решению дифференциального уравнения

 

EIyIV + q = 0

(2.124)

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

 

′′

 

y(0) = y (0) = 0,

(2.125)

y(l) = y (l) = 0.

′′

 

Ограничиваясь двучленным приближением, будем искать решение уравнения в виде

y = C U

1

+C U

2

= C sinπx

+C

2

sin

3πx

.

(2.126)

 

1

2

1

l

 

 

l

 

 

 

 

 

 

 

 

 

 

Выбранные функции удовлетворяют всем перечисленным выше требованиям. Они линейно независимы, непрерывно дифференцируемы и удовлетворяют граничным условиям (2.125).

Подставляя (2.126) в уравнение (2.124), получим невязку

 

 

 

 

EIπ 4

 

 

 

πx

 

 

 

3πx

+ q .

(2.127)

r(x,C

,C

 

) =

 

C

sin

 

+81C

 

sin

 

 

 

l 4

l

 

l

1

 

2

 

 

1

 

 

2

 

 

 

 

Решение методом Бубнова

Из условия ортогональности невязки (2.127) к выбранным координатным функциям согласно (2.120) следует

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

81

 

 

 

 

 

 

 

 

 

 

 

 

 

 

l

EIπ 4

 

 

 

πx

 

 

 

 

 

 

 

 

3πx

 

πx

 

 

 

 

 

 

 

 

 

C1 sin

 

 

+

81C2 sin

 

 

 

 

 

 

 

+ q sin

 

 

dx = 0,

 

 

l

4

 

 

l

 

 

 

 

l

 

l

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

l

EIπ

4

 

 

 

πx

 

 

 

 

 

 

 

 

3πx

 

3πx

 

 

 

 

 

 

C1 sin

 

 

 

+

81C2 sin

 

 

 

 

 

 

 

 

+ q sin

 

 

 

dx = 0.

 

l

4

 

 

 

 

l

 

 

 

 

 

l

 

l

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

После интегрирования получим систему

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

C

 

+

8ql

5

 

 

 

 

= 0,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

1

π 5EI

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

81l C2

 

2ql

5

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

+

 

 

= 0,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

3π

 

5

EI

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

решая которую, находим

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

C

= −

16ql 4

,

 

 

C

2

 

= −

4ql 4

 

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

π

5EI

 

 

 

 

 

 

 

 

 

243π 5EI

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Подставляя C1

и

C2

в искомое решение (2.126), получим при x = l/2

ymax = y(l / 2) = −0,01301 ql 4 .

EI

Точное значение максимального прогиба (с точностью до четырех значащих цифр) равно

ymax = y(l / 2) = −0,01302

ql 4

.

(2.128)

EI

 

 

 

Как видим, ошибка, полученная при решении методом Бубнова в двучленном приближении, составляет сотые доли процента.

Решение методом наименьших квадратов

Решение по-прежнему отыскиваем в форме (2.126). Интегральная квадратичная ошибка согласно (2.121) будет равна

 

 

 

 

 

 

 

 

b

2

 

 

 

 

 

 

 

 

l

EIπ 4

 

 

πx

 

3πx

 

2

S(C1,...,Cn )

=

r

 

(x,C1

, ,Cn )dx

=

 

 

 

 

 

C1 sin

 

 

+ 81C2 sin

 

 

+ q

dx =

 

 

l

4

l

 

l

 

 

 

 

 

 

 

 

a

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

EIπ

4

 

2

 

1

2

 

 

2

2

 

ql

4

 

 

 

 

4

 

 

 

108

 

 

 

 

 

 

 

 

 

 

 

 

 

81

 

 

 

 

 

 

 

 

 

 

 

 

 

 

=

 

4

 

 

 

 

2C1

+

 

C2

+

 

 

1

+

 

 

 

C1 +

 

C2

.

 

 

 

 

l

 

 

2

π 4

 

π

π

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

EI

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

82

 

 

 

 

 

 

Вычисляя производные

dS / dC1, dS / dC2

и приравнивая их нулю, полу-

чим систему

 

 

 

 

 

 

 

 

 

 

 

π C

+

4ql 4

= 0,

 

 

 

 

 

 

 

 

 

1

π

4EI

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ql 4

 

 

 

2

π C2

 

 

 

 

81

+108

 

 

 

 

= 0,

π

4

EI

 

 

 

 

 

 

 

 

 

 

 

решение которой дает те же значения постоянных C1 и C2 , что и в методе Бубнова:

C

= −

16ql 4

,

C

2

= −

4ql 4

.

 

 

1

 

π 5EI

 

 

 

243π 5EI

 

 

 

 

 

 

 

 

Поэтому результаты решения краевой задачи (2.124), (2.125) методами Бубнова и наименьших квадратов полностью совпадают.

Решение методом коллокаций

За точки коллокации выбираем точки x1 = l / 3 и

 

x2 = l / 2. Приравни-

вая нулю невязку (2.127) в точках коллокации, получим систему

 

 

 

C

 

 

 

 

 

 

ql 4

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3

 

+

 

 

= 0,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

2

 

 

π 4EI

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ql 4

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

C1

 

81C2 +

 

 

 

 

 

 

= 0.

 

 

 

 

 

 

 

 

 

π

4

EI

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

В результате находим

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2ql 4

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

ql

4

 

C

= −

 

,

 

C

 

 

=

 

 

 

3

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

3π 4EI

 

 

 

 

 

2

 

81 3 π

4 EI

 

По формуле (2.126) при

x =l / 2 с учëтом найденных констант получим

прогиб в среднем сечении балки

 

y(l / 2) = −0,01183

ql 4

.

 

 

EI

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Сравнение с точным решением (2.128) показывает, что ошибка, полученная при использовании метода коллокаций, составляет примерно 9% .

83

3. ОПЕРАЦИОННОЕ ИСЧИСЛЕНИЕ

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

3.1. Преобразование Лапласа

Преобразованием Лапласа называется преобразование, которое ставит в соответствие функции f(t) действительной переменной t функцию F(p) комплексной переменной p по формуле

 

F(p) = f (t)ept dt.

(3.1)

0

 

Несобственный интеграл в правой части формулы (3.1), зависящий от комплексного параметра p, называется интегралом Лапласа.

Интеграл сходится и действительно определяет собой некоторую функцию F(p), если подынтегральная функция f(t) удовлетворяет следующим условиям:

1.f(t) - кусочно - непрерывная функция;

2.f(t) = 0 при t < 0;

3.f(t) по абсолютной величине возрастает не быстрее заранее выбранной показательной функции, то есть можно найти такие постоянные M и α , что

f (t) < M eα t .

Число α называется показателем роста функции f(t).

Функция f(t) действительного аргумента t, удовлетворяющая перечисленным выше трем условиям, называется оригиналом , а функция F(p) комплексного аргумента p, определяемая согласно ( 3.1), называется изобра-

жением по Лапласу функции f(t)

или просто изображением .

Символически соответствие между оригиналом и изображением запи-

сывается обычно так:

 

f (t)=F(p)

или F(p)=f (t).

84

Совокупность всех оригиналов f(t) называется пространством оригиналов, а совокупность всех изображений F(p) – пространством изображений.

Изображения простейших функций

1. Пусть оригинал имеет вид:

0

при

t < 0,

f (t) =

при

t 0.

1

Эта функция широко применяется в приложениях и называется единичной функцией или функцией Хевисайда.

Изображение этой функции найдем, вычисляя интеграл Лапласа (3.1):

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

F(p) = ept dt = −

1

ep t

 

=

.

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

p

 

 

0

 

p

 

 

 

Следовательно,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

f (t) =1=

 

1

 

.

 

 

 

 

 

 

 

(3.2)

 

p

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2. Изображение показательной функции.

Пусть

f (t) = ea t , тогда

 

 

 

 

 

1

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

F(p) = eat ep t dt = e(pa)t dt = −

 

 

 

e(pa)t

=

.

 

p a

p a

0

0

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Следовательно,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

eat =

 

1

 

 

.

 

 

 

 

 

 

(3.3)

 

p a

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Нахождение изображений непосредственно по формуле (3.1) обычно затруднительно. Во многих случаях изображения по Лапласу могут быть найдены значительно проще, если воспользоваться свойствами преобразования Лапласа.

Свойства преобразования Лапласа

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