Математическое моделирование в задачах. Учебное пособие
.pdfДеформация сплошной среды |
71 |
5.12.У ê à ç à í è å. Тензоры скоростей деформаций (а)—(г) совпадают с полученными в задаче 5.11.
|
|
|
|
0 |
|
0 |
|
|
0 |
|
kx1 |
1 |
|
|||
5.13. а) |
D x1 |
; x2; x3 |
= |
|
0 |
|
|
0 |
|
kx2 |
; |
|||||
|
|
|
B |
|
|
1 |
|
kx |
2 |
|
|
C |
|
|||
|
|
|
|
B |
|
kx |
|
|
0 |
|
C |
|
||||
|
|
|
|
@ |
|
|
|
|
|
|
|
A |
|
|||
|
|
|
|
0 |
|
0 |
|
|
kx2 |
0 |
|
|
|
|||
б) |
D x1 |
; x2; x3 |
= |
|
kx2 |
0 |
|
01 |
; |
|
|
|||||
|
|
|
B |
|
|
|
|
|
C |
|
|
|
||||
|
|
|
|
B |
|
|
|
|
|
|
C |
|
|
|
||
|
|
|
|
@ |
|
|
|
|
|
|
A |
|
|
|
||
|
|
|
|
|
|
|
|
0 |
|
|
0 |
|
0 |
|
|
|
|
D x1; x2; x3 = 1 |
0k |
|
0 |
01 |
: |
|
|
|
|||||||
|
|
|
|
|
|
|
|
0 |
|
k |
0 |
|
|
|
|
|
в) |
|
|
|
|
|
|
B |
|
|
|
C |
|
|
|
|
|
|
2 |
|
0 |
|
|
|
|
|
||||||||
|
|
|
|
|
|
|
B0 |
|
0C |
|
|
|
|
|||
|
|
|
|
|
|
|
@ |
|
|
|
A |
|
|
|
|
|
|
|
|
1 |
|
|
|
|
0 |
|
0 |
|
|
|
|
|
|
5.14. а) |
D (r; '; z) = |
00 |
|
|
|
|
r2 |
|
01 |
; |
|
|
|
|
|
|
|
|
B |
|
|
|
|
|
|
C |
|
|
|
|
|
|
|
|
|
|
B |
|
|
|
|
|
|
C |
|
|
|
|
|
|
|
|
|
@ |
|
|
|
|
|
|
A |
|
|
|
|
|
|
|
|
|
0 |
|
|
|
|
0 |
|
0 |
|
|
|
|
|
|
б) D (r; '; z) = 0; |
|
|
|
|
|
|
r21 |
|
|
|
|
|||||
|
|
|
0 0 |
|
|
|
0 |
|
|
|
|
|||||
|
|
|
B |
|
|
|
|
|
|
|
C |
|
|
|
|
|
в) D(r; '; z) = B 0 |
|
|
0 |
0 C; |
|
|
|
|
||||||||
|
|
|
@ |
|
|
|
|
|
|
|
A |
|
|
|
|
|
|
|
|
r2 |
|
|
0 |
0 |
|
|
|
|
|
||||
5.15. У ê à ç à í è å. См. ответы к задаче 5.2.
72 |
Ответы, указания, решения |
Модели жидкостей
6.2.@@t + u@x@ + v @y@ + w @@z + @u@x + @y@v + @w@z = 0:
6.3.а) да; б) нет; в) да; г) нет.
6.4.а) нет; б) да; в) да; г) нет.
6.5.нет.
6.7. |
v2 = |
Ax1 |
x1=0 = 0: |
r2 |
6.9.а) да; б) да.
Модели жидкостей |
73 |
6.10.
@v1 |
@v1 |
|
@v1 |
|
|
@p |
|
|
||||||
|
+v1 |
|
|
|
+ v2 |
|
|
|
= |
|
|
+ |
|
|
@t |
@x1 |
@x2 |
@x1 |
|
||||||||||
|
+ 2 |
@2v1 |
+ |
|
@2v1 |
|
+ |
@2v2 |
+ f1 |
|||||
|
@(x1)2 |
@(x2)2 |
@x1@x1 |
|||||||||||
@v2 |
@v2 |
|
@v2 |
|
|
@p |
|
|
||||||
|
+v1 |
|
|
+ v2 |
|
|
= |
|
+ |
|
||||
@t |
@x1 |
@x2 |
@x2 |
|
||||||||||
|
+ |
@2v2 |
|
|
|
@2v2 |
|
|
@2v1 |
+ f2 |
||||
|
|
+ 2 |
|
+ |
|
|||||||||
|
@(x1)2 |
@(x2)2 |
@x1@x1 |
|||||||||||
6.11.(а)—(г) div P = (0; 0; 0):
6.12.а) (0; 0; 4k ):
б) Р å ø å í è å. В задаче 5.13, б найден тензор скоростей деформации
|
|
|
|
0 |
0 |
|
kx2 |
0 |
|
D x1 |
; x2 |
; x3 |
= |
kx2 |
0 |
01 |
; |
||
|
|
|
|
B |
0 |
|
0 |
C |
|
|
|
|
|
B |
|
0C |
|
||
|
|
|
|
@ |
|
|
|
A |
|
который найден в п. (б) задачи 5.13.
Тензор напряжений для вязкой несжимаемой жидкости определяется как функция тензора скоростей деформаций и давления P = pI + 2 D: Для случая нулевого давления и известного тензора скоростей деформации тензор напряжений будет равен
|
|
|
0 |
0 |
2 |
kx2 |
0 |
|
P x1 |
; x2; x3 |
= |
2 kx2 |
0 |
01 |
: |
||
|
|
|
B |
0 |
|
0 |
C |
|
|
|
|
B |
|
0C |
|
||
|
|
|
@ |
|
|
|
A |
|
74 |
|
|
|
|
Ответы, указания, решения |
|||
Дивергенция тензора в |
декартовой системе координат |
|||||||
может быть найдена как |
|
|
|
|
|
|
|
|
0 |
@P11 |
+ |
@P12 |
+ |
@P13 |
1 |
||
@x1 |
|
@x2 |
|
@x3 |
||||
B |
|
|
|
|
|
|
|
C |
B |
|
|
|
|
|
|
|
C |
div P = B |
@P21 |
+ |
@P22 |
+ |
@P23 |
C; |
||
@x1 |
|
@x2 |
|
@x3 |
||||
B |
1 |
+ |
2 |
+ |
3 |
C |
||
B |
@P31 |
|
@P32 |
|
@P33 C |
|||
@ |
|
|
|
|
|
|
|
A |
B |
@x |
|
|
@x |
|
|
@x |
C |
что дает значение div P = ( 2k ; 0; 0).
в) (0; 0; 0):
6.13.( k ; 0; 0):
6.14.(0; k ; 0):
6.15.( 2k ; 0; 0):
6.16.( k ; 0; 0):
Модель упругого тела |
75 |
Модель упругого тела
7.1.
@v1 |
@v1 |
|
@v1 |
|
|
|
|
|
|
|
|
|
|||
|
+ v1 |
|
+ v2 |
|
|
= |
|
|
|
|
|
|
|
|
|
|
1 |
|
2 |
|
|
|
|
|
|
|
|
|
|||
@t |
@x |
|
@x |
|
|
|
|
|
|
|
|
|
|||
|
|
|
= ( + ) |
|
@2u1 |
|
|
@2u2 |
|
||||||
|
|
|
|
|
|
|
+ |
|
|
|
|||||
|
|
|
@(x1)2 |
|
@x1x2 |
||||||||||
|
|
|
+ |
@2u1 |
|
+ |
@2u1 |
|
|
|
|
||||
|
|
|
@(x1)2 |
@(x2)2 |
|
||||||||||
@v2 |
@v2 |
|
@v2 |
|
|
|
|
|
|
|
|
|
|||
|
+ v1 |
|
+ v2 |
|
|
= |
|
|
|
|
|
|
|
|
|
|
1 |
|
2 |
|
|
|
|
|
|
|
|
|
|||
@t |
@x |
|
@x |
|
|
|
|
|
|
|
|
|
|||
|
|
|
= ( + ) |
|
@2u1 |
|
|
@2u2 |
|
||||||
|
|
|
|
+ |
|
||||||||||
|
|
|
@x1x2 |
@(x2)2 |
|||||||||||
|
|
|
+ |
@2u2 |
|
+ |
@2u2 |
|
|
|
|
||||
|
|
|
@(x1)2 |
@(x2)2 |
|
||||||||||
+
+
7.2.
01
( + 2 )a |
0 |
0 |
C |
P = B 0 |
a |
0 |
@A
00 a
7.3.а) ( kt 1; 0; 0);
б) ( kt 1; kt 2; kt 3);
в) (0; kt 1; 0);
8u1 = 1(cos(k 3t) 1) + 2 sin(k 3t);
>
>
>
<
г) u2 = 1 sin(k 3t) + 2(cos(k 3t) 1);
>
>
>
:u3 = 0:
76 |
|
|
|
|
Ответы, указания, решения |
7.4. а) " = |
0 0 |
0 |
0 |
1 |
; |
|
kt |
0 |
0 |
C |
|
|
B |
|
|
|
|
|
@ |
|
|
A |
|
|
0 |
0 |
0 |
|
|
01
kt |
0 |
0 |
C; |
б) " = B 0 |
kt |
0 |
@A
0 0 kt
|
|
0 |
0 |
kt |
0 |
1 |
|
в) |
" = 1 |
B |
kt |
0 |
0 |
C |
; |
2 |
0 |
0 |
0 |
|
|||
|
|
@ |
|
|
|
A |
|
г) Р å ø å í è å. Вычислим компоненты ui вектора смещений u u1 = x1 1 = 1 cos k 3t 1 + 2 sin k 3t ;
u2 = x2 2 = 1 sin k 3t + 2 cos k 3t 1 ; u3 = x3 3 = 0:
Компоненты тензора малых деформаций в приближении x
рассчитываются по формуле
|
|
|
|
1 |
|
@ui |
@uj |
|||||
|
|
|
"ij = |
|
|
|
|
+ |
|
|
: |
|
|
|
|
2 |
@ j |
@ i |
|||||||
Получаем |
|
|
|
|
|
|
|
|
|
|||
|
|
@u1 |
|
|
|
|
|
|
@u2 |
|||
"11 |
= |
|
= cos k 3t 1; |
|
|
|
"22 = |
|
|
= cos k 3t 1; |
||
@ 1 |
|
|
|
@ 2 |
||||||||
"33 |
= |
@u3 |
= 0: |
|
|
|
|
|
|
|
|
|
@ 3 |
|
|
|
|
|
|
|
|
|
|||
Разложим cos( ) в ряд Тейлора и отбросим слагаемые выше первого порядка ( 2; 4 и т. д.), тогда получим
"11 |
=1 + |
1 |
k 3t |
2 |
+ |
1 |
k 3t |
4 |
+ ::: 1 |
0; |
|
|
|
|
|
|
|||||||
2! |
|
4! |
|
||||||||
"22 |
0; |
|
"33 = 0: |
|
|
|
|
|
|
||
Модель упругого тела |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
77 |
|||||
Найдем остальные элементы тензора. |
|
|
|
|
|||||||||||||||||
|
1 |
|
|
|
|
@u1 |
|
@u2 |
1 |
|
|
|
|
|
|||||||
"12 = "21 = |
|
|
|
|
|
|
|
|
+ |
|
|
|
= |
|
|
sin(k 3t) |
sin(k 3t) = 0; |
||||
2 |
|
|
@ 2 |
@ 1 |
2 |
||||||||||||||||
|
1 |
|
|
|
|
@u1 |
@u3 |
|
|
|
|
|
|
|
|||||||
"13 = "31 = |
|
|
|
|
+ |
|
= |
|
|
|
|
|
|
|
|||||||
2 |
|
@ 3 |
@ 1 |
|
|
|
|
|
|
|
|||||||||||
= |
kt |
1 sin(k 3t) + 2 cos(k 3t) |
= |
1 |
ktx2; |
||||||||||||||||
|
|
|
|
|
|||||||||||||||||
2 |
|
|
2 |
||||||||||||||||||
|
1 |
|
|
|
|
@u2 |
@u3 |
|
|
|
|
|
|
|
|||||||
"23 = "32 = |
|
|
|
|
+ |
|
= |
|
|
|
|
|
|
|
|||||||
2 |
|
@ 3 |
@ 2 |
|
|
|
|
|
|
|
|||||||||||
=kt2 1 cos(k 3t) 2 sin(k 3t) = 12ktx1:
Иокончательно
|
|
|
0 |
|
|
0 |
kx2t |
|
|
|
" = |
1 |
0 |
0 |
|
|
0 |
|
kx1t1 |
: |
|
|
|
|
|
|||||||
|
2 |
B |
2 |
t |
|
1 |
|
0 |
C |
|
|
|
Bkx |
|
kx t |
|
C |
|
|||
|
|
@ |
|
|
|
|
|
A |
|
|
7.5.
а) P = 0 ( |
0 |
|
|
kt |
0 |
1; |
||||
|
|
B |
|
+ 2 )kt |
|
0 |
0 |
C |
||
|
|
|
0 |
|
|
0 |
kt |
|||
|
|
@ |
|
|
|
|
|
|
|
A |
|
|
|
|
0 |
|
|
|
|
|
1 |
|
|
|
|
kt |
|
0 |
0 |
C; |
||
б) P = ( + 2 ) B 0 |
|
kt |
0 |
|||||||
|
|
|
|
@ |
|
|
|
|
kt |
A |
|
|
0 |
|
0 |
|
|
0 |
|
||
|
P = |
0 |
kt |
0 |
1 |
; |
|
|
||
в) |
B |
kt |
0 |
0 |
C |
|
|
|||
|
0 |
0 |
0 |
|
|
|
||||
|
|
@ |
|
|
|
A |
|
|
|
|
78 |
Ответы, указания, решения |
|
|
0 |
0 |
|
0 |
kx2t |
|
|
|
г) |
P = |
0 |
|
0 |
|
kx1t1 |
: |
||
|
B |
2 |
|
1 |
0 |
C |
|
||
|
|
Bkx t |
|
kx t |
|
C |
|
||
|
|
@ |
|
|
|
|
|
A |
|
7.6. а) P hni = ( ( + 2 )kt; 0; 0); P hmi = (0; kt; 0);
P hki = (0; 0; kt);
б) P hni = ( ( + 2 )kt; 0; 0); P hmi = (0; ( + 2 )kt; 0); P hki = (0; 0; ( + 2 )kt);
в) P hni = (0; kt; 0); P hmi = ( kt; 0; 0); P hki = (0; 0; 0);
г) P hni = (0; 0; 2kt); P hmi = (0; 0; kt); P hki = (2kt; kt; 0):
Библиографический список
1.Жермен П. Курс механики сплошных сред. Общая теория. М.: Высшая школа, 1983.
2.Ильюшин А. А. Механика сплошной среды. М.: МГУ, 1990.
3.Мeйз Дж. Теория и задачи механики сплошных сред. М.: Мир, 1974.
4.Самарский А. А., Михайлов А. П. Математическое моделирование. М.: Физматлит, 2002.
5.Седов Л. И. Механика сплошной среды: В 2 т. М.: Наука, 1973.
6.Хакимзянов Г. С., Чубаров Л. Б., Воронина П. В. Математическое моделирование. Новосибирск: НГУ, 2014.
7.Основания устройства и конструкция орудий и боеприпасов наземной артиллерии. Воениздат, 1976.
Учебное издание
Воронина Полина Владимировна,
Лапин Василий Николаевич
МАТЕМАТИЧЕСКОЕ МОДЕЛИРОВАНИЕ В ЗАДАЧАХ
Учебное пособие
Редактор С. В. Исакова
Подписано в печать 05.04.2023 г. Формат 60 x 84 1/16. Уч.-изд. л. 5. Усл. печ. л. 4,6.
Тираж 160 экз. Заказ • 46
Издательско-полиграфический центр НГУ. 630090, Новосибирск-90, ул. Пирогова, 2.
