close

Вход

Забыли?

вход по аккаунту

?

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

код для вставкиСкачать
168
НАУЧНЫЕ ВЕДОМОСТИ
Серия Математика. Физика. 2016 № 20 (241). Выпуск 44
_________________________________________________________________
УДК 517.958
ВЫЧИСЛИТЕЛЬНАЯ СХЕМА РЕШЕНИЯ НЕКОТОРЫХ ЗАДАЧ
ИЗ ТЕОРИИ СОПРОТИВЛЕНИЯ МАТЕРИАЛОВ
THE COMPUTATIONAL SCHEME FOR SOLUTION OF SOME PROBLEMS
FROM THE THEORY OF STRENGTH OF MATERIALS
Н.А. Чеканов 1, И.Н. Беляева 1, Н.Н. Чеканова 2
N.A. Chekanov 1, I.N. Belyaeva 1, N.N. Chekanova 2
1)Белгородский
национальный исследовательский университет,
Россия, 308015, г.Белгород, ул. Победы, 85
Belgorod National Research University,
85 Pobedy St, Belgorod, 308015, Russia
2)
Харьковский институт банковского дела Университета банковского дела НБУ,
Украина, 61174, г. Харьков, прт. Победы, 55
Kharkov Institute of Banking of National University of Banking,
55 av. Pobedy, Kharkov, 61174, Ukraine
E-mail: сhekanov@bsu.edu.ru;ibelyaeva@bsu.edu.ru; chekanova76@list.ru
Аннотация. Разработан алгоритм и составлена программа в среде MAPLE для решения обыкновенных
дифференциальных уравнений IV порядка, в общем, в виде обобщенных степенных рядов. Дифференциальные
уравнения могут содержать регулярные особые точки. Приведены примеры решений дифференциальных
уравнений IV порядка из теории сопротивления материалов, которые показывают эффективность разработанной программы.
Resume. The developed algorithm and program in MAPLE of the solution ordinary differential equations IV
order , in general, in the form of generalized power series are developed. The differential equation could consist the
singular regular points. Some examples of the solution differential equations IV order in theory of resistance of materials are presented, that show the efficiency of the developed program.
Ключевые слова: дифференциальные уравнения IV порядка, алгоритм, математическое моделирование, обобщенно степенные ряды, регулярные особые точки
Key words: differential equations IV order, algorithm, mathematic modeling, generalized power series, singular regular points
Введение
Эффективным методом интегрирования дифференциальных уравнений является поиск решений в виде степенных рядов и последующим нахождением коэффициентов этих рядов [1]. Но на
практике при конкретных расчетах приходится сталкиваться с большим объемом вычислений при
нахождении неизвестных коэффициентов степенных рядов, причем сложность нахождения решений увеличивается в тех точках, в которых имеются особенности. Однако использование современных компьютеров вместе с пакетами программ для аналитических вычислений, таких как Maple,
Mathematica, Reduce и другие, позволяют достаточно быстро выполнить необходимые расчеты по
построению решений линейных дифференциальных уравнений в виде рядов, причем до очень
больших порядков.
В задачах математической физики наиболее часто возникают линейные дифференциальные
уравнения второго и четвертого порядков. Следует также отметить, что решения некоторых нели-
НАУЧНЫЕ ВЕДОМОСТИ
Серия Математика. Физика. 2016. № 20 (241). Выпуск 44
169
________________________________________________________________
нейных дифференциальных уравнений могут быть выражены через линейно независимые решения
соответствующих линейных дифференциальных уравнений. Например, решение нелинейного уравнения Ермакова-Милна-Пинни строится из двух независимых решений линейного обыкновенного
дифференциального уравнения второго порядка.
В работе представлена вычислительная схема решения обыкновенных линейных дифференциальных уравнений четвертого порядка в виде обобщенных степенных рядов с использованием
системы компьютерной алгебры Maple. Также с помощью составленной программы были найдены
линейно-независимые решения для ряда конкретных дифференциальных уравнений.
Вычислительная схема
Рассмотрим дифференциальное уравнение четвертого порядка
y ( IV ) ( x)  P3 ( x) y ( x)  P2 ( x) y ( x)  P1 ( x) y ( x)  P0 ( x) y( x)  0 .
(1)
В случае, если коэффициенты-функции P3 ( x) , P2 ( x) , P1 ( x) , P0 ( x) не содержат регулярных
особых точек и являются голоморфными функциями в окрестности точки x  x0
P3 ( x) 


k 0
P1 ( x) 


k 0
p k(3) ( x  x0 ) k , P2 ( x) 
p k(1) ( x  x0 ) k , P0 ( x) 

 pk(2) ( x  x0 ) k ,
k 0

 pk(0) ( x  x0 ) k ,
(2)
k 0
то линейно независимые решения y1 , y 2 , y 3 и y 4 могут быть представленными виде степенных
рядов:
y1 ( x)  1 

 ck(1) ( x  x0 ) k ,
(3а)
k 2
y 2 ( x)  x  x 0 

 ck(2) ( x  x0 ) k ,
(3б)
k 2
y 3 ( x)  ( x  x 0 ) 2 / 2 

 ck(3) ( x  x0 ) k ,
(3в)
k 2
y 4 ( x)  ( x  x 0 ) 3 / 3 

 ck(4) ( x  x0 ) k ,
(3г)
k 2
(1)
( 2)
( 3)
( 4)
Коэффициенты c k , c k , c k , c k
определяются единственным образом посредством под-
становки рядов (3) в уравнение (1) и приравниванием к нулю коэффициентов при различных степенях независимой переменной в левой части полученного равенства.
При наличии полюсов не выше четвертого порядка в точке x  x0 вид решений (3) будет
иным в зависимости от корней определяющего уравнения (см., например, [1-4]). Из теории обыкновенных дифференциальных уравнений [2] известно, что для того чтобы уравнение, в частности, вида (1) имело в окрестности особой точки x  x0 хотя бы одно частное решение в виде обобщенного
степенного ряда
170
НАУЧНЫЕ ВЕДОМОСТИ
Серия Математика. Физика. 2016 № 20 (241). Выпуск 44
_________________________________________________________________
y ( x)  ( x  x 0 ) 

 c k ( x  x0 ) k ,
(c0  0) ,
(4)
k 0
где показатель

есть некоторое постоянное число, достаточно, чтобы это уравнение имело вид

y ( IV ) ( x) 

k 0

p k(3)
x  x0
y ( x) 

k 0


p k( 2)
( x  x0 )
2
y ( x) 
k 0

 pk(0)
p k(1)
( x  x0 )
3
y ( x) 
k 0
( x  x0 ) 4
y ( x)  0 .
(5)
Показатель  находится из определяющего уравнения:
 (   1) (   2) (   3)   (   1) (   2) p03   (   1) p02   p10  p00  0 .
(6)
Пусть 1 ,  2 ,  3 и  4 есть корни уравнения (6). Тогда, если корни определяющего уравнения 1 ,  2 ,  3 и  4 независимы, и никакие два из них не различаются на целое число, то каждому
числу  соответствует определенная последовательность коэффициентов c k , а всего получается
четыре независимых решения, образующих фундаментальную систему
y1 ( x)  ( x  x0 ) 1

 ck(1) ( x  x0 ) k ,
k 0
y 2 ( x)  ( x  x 0 )  2

 ck(2) ( x  x0 ) k ,
k 0
y 3 ( x)  ( x  x 0 ) 3

 ck(3) ( x  x0 ) k ,
k 0
y 4 ( x)  ( x  x 0 )  4

 ck(4) ( x  x0 ) k ,
k 0
1
2
3
(c0(1)  0) ,
(7а)
(c0( 2)  0) ,
(7б)
(c0(3)  0) ,
(7в)
(c0( 4)  0) .
(7г)
4
Коэффициенты c k , c k , c k и c k определяются подстановкой рядов (7) в уравнение (5), при
1
2
3
4
этом коэффициенты c 0 , c 0 , c 0 и c 0 остаются произвольными (далее положим их равными единице).
Если найденные четыре значения  таковы, что два или несколько различаются на целое
число, то они могут быть расположены в виде независимых последовательностей
1, 2 , ,  1 ,
 ,  1, ,  1
,

(8)
так, чтобы величины в каждой последовательности различались только на целые числа, и их вещественные части образовали бы невозрастающую последовательность. Только первый член каждой
последовательности дает решение вида (4), поскольку, например, любой член    k последовательности     1 равен   или меньше его на положительное целое число.
Рассмотрим одну из последовательностей показателей дифференциального уравнения,
например, последовательность
НАУЧНЫЕ ВЕДОМОСТИ
Серия Математика. Физика. 2016. № 20 (241). Выпуск 44
171
________________________________________________________________
1 ,  2 , ,  1 ,
которая так расположена, что если    , то      положительное целое число или нуль. Поскольку эти показатели не обязательно равны, они могут быть разделены на подпоследовательности
так, чтобы члены каждой подпоследовательности были равны между собой. Так, предположим, что
1   2     i 1 соответствуют корню кратности i ;  i   i 1     j 1 соответствуют корню
кратности j  i ;  j   j 1     k 1 соответствуют корню кратности k  j и т.д. до тех пор пока
ряд не будет исчерпан.
Рассмотрим показатель 1 в первой подпоследовательности. В данном случае возникает
подпоследовательность
i решений
y1 ( x)  w1 ( x, 1 ) ,
y2 ( x)  w1 ( x, 1 ) ln x  w2 ( x, 1 )
………………………………………………….
yi 1 ( x)  w1 ( x, 1 )(ln x) i 1  (i  1)w2 ( x, 1 )(ln x) i 2    wi 1 ( x, 1 ) ,
где ws ( x,  s )  ( x  x0 )
s
(9)

 cks ( x  x0 ) k , (c0s  0, s  1,2,3,...) . Наличие w1 ( x, 1 )(ln x) r 1 в Wr показы-
k 0
вает, что эти i решений линейно-независимы.
Рассмотрим показатель  i второй подпоследовательности, ему соответствуют j  i решений
yi ( x)  w1 ( x,  i )(ln x) i  i w2 ( x,  i )(ln x) i 1    wi ( x,  i )
………………………………………………………...
y j 1 ( x)  w1 ( x,  i )(ln x) j 1  ( j  1)w2 ( x,  i )(ln x) j 2    w j 1 ( x,  i ) . (10)
Аналогично подпоследовательность индекса j дает k  j решений и т.д. до тех пор пока ряд
не будет исчерпан.
Так как функции y1 ( x) , y 2 ( x) , y3 ( x) и y 4 ( x) являются линейно независимыми решениями
уравнения (5), то с их помощью находим общее решение
y( x)  C1  y1 ( x)  C2  y2 ( x)  C3  y3 ( x)  C4  y4 ( x) .
(11)
Разработанная программа [5] позволяет находить решения дифференциальных уравнений
четвертого порядка в виде степенных рядов, в общем, произвольной степени n , но ограниченной
возможностями конкретной вычислительной машины.
Решение задач из теории сопротивления материалов
При помощи разработанной программы найдены решения двух задач из теории сопротивления материалов.
Пример 1. Основным уравнением устойчивости плоской формы изгиба двутавровой балки,
к которой приложен изгибающий момент и внешняя сила, является линейное дифференциальное
уравнение четвертого порядка [6]
172
НАУЧНЫЕ ВЕДОМОСТИ
Серия Математика. Физика. 2016 № 20 (241). Выпуск 44
_________________________________________________________________
1
 ( IV ) ( x) 
a
2
 ( x) 
1
d4
 ( x)  0 ,
где a , d - геометрические и прочностные характеристики балки,  (x) - узел закручивания балки
вдоль ее длины. В данном уравнении коэффициенты-функции равны P3 ( x)  0 , P2 ( x)  
P1 ( x)  0 , P0 ( x)  
1
1
a2
,
, не имеют особых точек.
d4
С помощью разработанной программы [5] получаем четыре линейно независимых решения первые
члены которых приведены ниже
1 ( x)  1 
1
24 d
4
1
x4 
2
720 a d

4
(d 4  a 4 )
x6 
40320a d
(d 8  4 a 4 d 4  3 a 8 )
87178291200 a
 2 ( x)  x 
4
1
120 d

4
d
12
1
x5 
2
5040a d
(d  3 a d
8
10
4
4
4
 3 ( x) 
6
8
x7 
12
(d 8  3 a 4 d 4  a 8 )
x10 
8
479001600a d
(d 12  5 a 4 d 8  6 d 4 a 8  a 12 )
20922789888000a
a )
6227020800a d
(d 4  2 a 4 )
3628800a d
x14 
8
8
8
x8 
(d 4  a 4 )
4
362880a d
8
x9 
12
d
16
(d 4  2 a 4 )
6
39916800a d
(d  4 a d 4  3 a 8 )
8
x13 
8
12
x12 
x16 ,
x11 
4
1307674368000 a
10
d
12
x15 ,
(d 4  a 4 ) 6 (d 4  2 a 4 ) 8 (d 8  3 a 4 d 4  a 8 ) 10 (d 8  4 a 4 d 4  3 a 8 ) 12
x2
1
4

x

x 
x 
x 
x 
2 24 a 2
720 a 4 d 4
40320a 6 d 4
3628800a 8 d 8
479001600a10 d 8

(d 12  5 a 4 d 8  6 d 4 a 8  a12 )
87178291200 a
 4 ( x) 
12
d
12
x14 
(d 12  6 a 4 d 8  10 d 4 a 8  4 a12 )
20922789888000a
14
d
12
x16 ,
(d 4  2 a 4 ) 9 (d 8  3 a 4 d 4  a 8 ) 11
(d 4  a 4 ) 7
x3
1
5

x

x

x 
x 
3 60 a 2
2520a 4 d 4
181440a 6 d 4
19958400a 8 d 8

(d 8  4 a 4 d 4  3 a 8 )
3113510400a 10 d 8
x 13 
(d 12  5 a 4 d 8  6 d 4 a 8  a 12 )
653837184000 a 12 d 12
x 15 .
НАУЧНЫЕ ВЕДОМОСТИ
Серия Математика. Физика. 2016. № 20 (241). Выпуск 44
173
________________________________________________________________
Пример 2. Уравнение
y ( IV ) ( x) 
2F

y ( x) 
M2
2
y ( x) 
F2
2
y ( x)  0
описывает вал, на который действует крутящий момент M и сила сжатия F , приложенная вдоль
оси вала, вал имеет постоянное поперечное сечение и  - жесткость вала на изгиб [7].
Данное уравнение не имеет особых точек. С помощью разработанной программы [5] получены четыре линейно независимых решения первые члены которых приведены ниже
y1 ( x)  1 

y 3 ( x) 
40320
F2
120 2
F3
60
3
x5 
6
x5 
F3
360 3
362880
6
F 2 (4 F 2  M 2 )
720
x8 
x6 
F 2 (16 F 4  12 M 2 F 2  M 4  F 2 2 )
4
x6 
1260
5
x7 
F 3 (16 F 4  16 M 2 F 2  3 M 4  2 F 2 2 )
181440
F 2 (4 F 2  M 2 )
x9 
F 3 (2 F 2  M 2 )
5040 4
x7 
7
F 3 (2 F 2  M 2 )
10080 5
x8 
F 3 (16 F 4  16 M 2 F 2  3 M 4  2 F 2 2 )
1814400
7
x9 ,
x10 ,
x2
M 2 4 F M 2 5 (4 M 2 F 2  M 4  F 2 2 ) 6 F (4 M 2 F 2  2 M 4  F 2 2 ) 7

x 
x 
x 
x 
2 24 2
60 3
720 4
2520 5


y 4 ( x) 
24
2
x4 
F 2 (16 F 4  12 M 2 F 2  M 4  F 2 2 )
y 2 ( x)  x 

F2
(16 F 4 M 2  12 M 4 F 2  4 F 4 2  M 6  2M 2 F 2 2 )
40320 6
x8 
F (16 F 4 M 2  16 M 4 F 2  4 F 4 2  3 M 6  4 M 2 F 2 2 )
181440
7
x9 ,
x3
F 4 (4 F 2  M 2 ) 5 F (2 F 2  M 2 ) 6 (16 F 4  12 M 2 F 2  M 4  F 2 2 ) 7

x 
x 
x 
x 
3 6
60 2
90 3
2520 4


F (16 F 4  16 M 2 F 2  3 M 4  2 F 2 2 )
10080
5
x8 
(64 F 6  80 F 4 M 2  24 M 4 F 2  12 F 4 2  M 6  2 M 2 F 2 2 )
181440
6
x9 .
174
НАУЧНЫЕ ВЕДОМОСТИ
Серия Математика. Физика. 2016 № 20 (241). Выпуск 44
_________________________________________________________________
Заключение
В данной работе представлена вычислительная схема решения обыкновенных дифференциальных уравнений четвертого порядка в виде обобщенно степенных рядов. На основе разработанного алгоритма составлена программа в Maple и решены две задачи из теории сопротивления материалов. В дальнейшем планируется построение функции Грина для дифференциальных уравнений IV
порядка.
Список литературы
1. Матвеев Н.М. 1963. Методы интегрирования обыкновенных дифференциальных уравнений. М.:
Высшая школа: 546.
Matveev N.M. 1963. Methods integrating of ordinary differential equations. M., Vysschhaya schkola: 546.
2. Трикоми Ф. 1962. Дифференциальные уравнения. М. Издательство иностранной литературы: 351.
Tricomi F. 1961. Differential equations. Turin: Blackie & son limited : 348.
3. Сансоне Дж. 1953. Обыкновенные дифференциальные уравнения. Т.1, М.: Изд-во иностранной литературы : 346.
Sansone J. 1948. Ordinary differential equations. V.1 : 346.
4. Айнс Э. Л. 1939. Обыкновенные дифференциальные уравнения, ОНТИ: Государственное научнотехническое издательство Украины : 719.
Ince E.L. 1939. Differential equations. London: University Press: 700.
5. Беляева И.Н., Чеканов Н.А., Чеканова Н.А. 2016. Программа символьно-численного интегрирования
линейного дифференциального уравнения четвертого порядка в виде обобщенных рядов. Свидетельство о государственной регистрации программы для ЭВМ №2016611952.
Belyaeva I.N., Chekanov N.A., Chekanova N.N. 2016. Program of symbol-numeric integration of linear differential equation of four order. Patent of RU, Program for ECM, №2016611952.
6. Феппль А., Феппль Л. 1936. Сила и деформация. Прикладная теория упругости, том 2, М-Л.: НКТЛ:
409 с.
Föppl A., Föppl L. 1928. Drang und Zwang, v.II., München - Berlin: Verlag von R. Oldenbourg: 385.
7. Коллатц Л. 1968. Задачи на собственные значения с техническими приложениями, М.: Наука, гл.
ред. физ.-мат. Литературы: 504.
Collatz L. 1963. Eigenwertaufgaben mit technischen anwendungen, Leipzig: Akademishe verlagsgesellshaft,
Geest & Portig K. : 495.
Документ
Категория
Без категории
Просмотров
4
Размер файла
1 613 Кб
Теги
решение, сопротивления, материалы, вычислительной, некоторые, задачи, схема, теория
1/--страниц
Пожаловаться на содержимое документа