МДТТ дз1
.pdfГосударственное образовательное учреждение высшего профессионального образования
«Московский государственный технический университет имени Н.Э. Баумана»
(МГТУ им. Н.Э. Баумана)
ФАКУЛЬТЕТ |
Специальное машиностроение |
КАФЕДРА |
Космические аппараты и ракеты-носители |
Д О М А Ш Н Е Е З А Д А Н И Е № 1
По курсу:
Механика деформированного твердого тела
Вариант № 10
Студент |
СМ1-81 |
|
Г.С. Нахапетян |
|
(Группа) |
|
(И.О. Фамилия) |
Преподаватель |
СМ1 |
|
Б.С. Сарбаев |
|
(Кафедра) |
|
(И.О. Фамилия) |
Москва 2015
Условие задания
Прямой композитный стержень (см. рис. 1) изготовлен из многослойного композиционного |
|||
материала со схемой армирования ± / 90 / 0 . Упругие характеристики таковы: |
|||
= 150 ГПа |
= 8 ГПа |
= 5 ГПа |
= 0.25 |
Геометрические размеры поперечного сечения показаны на рис. 2. Длина стержня и величины внешней нагрузки приведены в исходных данных.
Для указанного стержня необходимо:
1)Определить перемещение и осевую силу
2)Воспользоваться точным аналитическим решением
3)Построить эпюры перемещения и осевой силы
4)Определить сечение, в котором действует наибольшая осевая сила
5)Для данного сечения рассчитать напряжения и деформации в монослоях
Исходные данные
Рис. 1 – Расчетная схема
|
Рис. 2 – Профиль поперечного сечения |
|
|
||
= 0.003 |
! |
= 0.008 = 0.7 м |
= 30 |
||
|
|
|
Решение
Пункт 1
Составление аналитического решения для определения перемещения и осевой силы
Для решения поставленной задачи воспользуемся уравнением статики, записанного для плоского деформированного твердого тела:
|
'( |
+ |
'( |
+ + = 0 |
|
||
|
& ') |
') |
|
||||
|
%'( |
+ |
'( |
+ + = 0 |
|
||
|
$ ') |
') |
|
||||
+ , + − объемные силы, действующие в направлении ) , ) . |
= 0), а также учитывая, |
||||||
C учетом того, что задача линейная, в направлении ) (т.е. ( = ( = + |
|||||||
что ( = ( ) , получим: |
.( |
+ + = 0 |
1 |
|
|||
Осевая сила есть: |
.) |
|
|
||||
|
|
|
/ = ( |
|
|||
− площадь поперечного сечения |
|
|
|
||||
|
( = |
∙ 1 |
|
||||
|
|
|
|
||||
|
|
|
|
|
.2 |
|
|
Объединив выражения, получим: |
|
|
|
1 = .) |
|
||
|
|
|
|
.2 |
|
|
|
|
|
|
|
|
2 |
|
|
|
( = .) |
|
|||||
|
|
|
|
|
.2 |
3 |
|
|
/ = .) |
|
|||||
Так как рассматривается лишь одно направление ) , то для упрощения записей примем: |
|||||||
) = ); ( = ( |
; = |
; 1 = 1 ; 2 = 2 = 2; |
+ = + ; |
||||
Получим: |
.( |
|
|
|
|
|
|
|
+ + |
= 0 |
4 |
|
|||
|
.) |
|
|||||
|
|
|
|
.2 |
5 |
|
|
|
( = .) |
|
|||||
|
5 = 678 |
9: |
; |
|
|||
|
97 |
|
Исходя из граничных и начальных условий удобно использовать различные виды уравнения (4). Выведем их:
|
.( |
. ( |
|
|
./ |
|
Первое уравнение получим путем умножения площади поперечного сечения на (4): |
||||||
Или: |
< .) + + = = |
.) |
+ + |
= |
.) + = 0 |
|
95 |
|
|
|
|
|
|
|
+ >7 = ?; |
|
@ |
|||
− |
97 |
|
||||
Линейная распределенная сила в направлении оси ) |
|
|||||
Второе уравнение получим путем подстановки (6) в (7): |
|
|
||||
|
678 |
9A: |
|
|
B |
|
|
97A + >7 = ? |
|
Применим уравнения (6)-(8) непосредственно к нашей задаче. Для этого, разобьем стержень на 2 участка:
Запишем граничные и начальные условия и условие стыковки участков, необходимые для определения констант интегрирования:
) = 0; 2 = 0 |
|
9 |
|
|
) = 3 ; / |
= −! |
= / |
10 |
11 |
) = 2 ; 2 = 2 |
; / |
|
Участок 1
=
..)/ + = 0 / = ..2)
Интегрируем:
/ = − ) + C
2 = − 2 ) + C ) + C
Используя 9 :
0 = − 2 0 + C ∙ 0 + C => EA = ?
Участок 2
= 0
..)/ = 0 / = ..2)
Интегрируем: |
|
|
|
|
|
|
|
|
/ = C |
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
Используя 10 : |
|
|
|
|
|
|
2 = C ) + CF |
|
|
|
||||||||
|
|
|
|
|
|
/ |
= EG = −H |
|
|
|
||||||||
|
|
|
|
|
|
|
|
|
|
|||||||||
|
|
|
|
|
|
|
|
Стыковка |
|
|
|
|
|
|
||||
Для начала перепишем полученные значения 2 , 2 , / , / для наглядности: |
||||||||||||||||||
|
|
|
|
|
|
|
/ |
= − ) + C |
|
|
|
|||||||
|
|
|
|
|
|
|
1 |
|
|
|
+ C )J |
|
|
|
||||
|
|
|
2 = I− 2 ) |
|
|
|
|
|||||||||||
|
|
|
|
|
|
|
|
/ = −! |
|
|
|
|
|
|
||||
Используя 11 : |
|
|
|
|
|
2 = − ! ) + CF |
|
|
|
|||||||||
1 |
|
2 |
|
|
|
|
|
|
|
! |
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
||||||||
|
K− 2 |
|
+ C ∙ 2 L = − 2 + CF |
|||||||||||||||
|
|
− ∙ 2 + C = −! |
|
|
|
|
|
|||||||||||
|
|
|
|
|
|
|
|
|
: |
|
|
|
|
|
||||
|
|
|
|
EO = A>P − H |
|
|
|
|
||||||||||
|
1 |
|
|
|
|
|
→ : |
! |
|
|
||||||||
|
|
K− 2 2 |
|
|
+ 2 − ! ∙ 2 L = − |
|
2 + CF |
|||||||||||
|
|
|
|
|
|
|
|
|
! |
|
|
! |
||||||
|
|
CF = −2 |
|
+ 4 |
|
− 2 |
|
+ 2 |
|
|
|
|
ER = A |
>PA |
|
|
|
|
678 |
||
Подставив численные значения из «исходные данные», получим: |
|||||
/ ) |
|
) |
|||
|
|
= − − 0.002 |
|||
2 ) = − |
) |
||||
2 − 0.002) |
|||||
/ ) |
= −0.008 = STUVW |
||||
|
|||||
2 ) |
= −0.008 ∙ ) + 0.006 ∙ |
||||
Упростим полученные выражения, перейдя к)новой переменной: |
|||||
Y = ; |
Y 0,3 |
||||
5O [ |
= −?. ??G ∙ [ − ?. ??A |
||||
678 |
|||||
:O [ = −?. ??O\ ∙ P ∙ [A − ?. ??A ∙ P ∙ [ = −O. ?\ ∙ O?]G ∙ [A − O. R ∙ O?]G ∙ [ |
|||||
5A [ |
= −?. ??B = ^_`ab |
||||
678 |
:A [ = −\. ; ∙ O?]G ∙ [ + R. A ∙ O?]G
Пункт 2
Построение эпюр перемещения и осевой силы
Обратим внимание на то, что участок ) 2 , 3 является опасным (т.е. в нем наблюдается max |
|||
осевой силы /cd). Рассмотрим сечение ) = 3 Y = 3 . Значение осевой силы: |
|||
/cd = −8 ∙ 10 |
] |
8 |
∙ 10] |
|
∙ ∙ = − 3 |
∙ 10] ≈ −2.667 ∙ |
Пункт 3
Расчет напряжений и деформаций в монослоях в опасном сечении
Распишем поэтапно решение поставленной задачи. Рассмотрим исходные данные для решения |
|||||||||||
задачи. Нам задана следующая схема армирования: |
|
|
|
|
|
||||||
|
|
±30 / 90 / 0 |
|
|
|||||||
Иначе говоря, нам заданы углы намотки f каждого монослоя в многослойном КМ, а также их |
|||||||||||
объемные содержания gf и количество U. Также известны характеристики материалов каждого |
|||||||||||
монослоя f, f, f, f, f = hhkjij |
f. Также нам дано усредненное значение вектора |
||||||||||
напряжения l(сn = o( , (p, (qrт, которое выглядит следующим образом: |
|||||||||||
( = |
/cd |
|
|
|
|
|
|
(p |
= (q = 0 |
||
|
|
= −2.667 ∙ |
|
|
|||||||
= 7 ∙ 80 + 60 − 7 ∙ 10]t = 9.31 ∙ 10]F |
м |
||||||||||
|
|
|
( ≈ −1.738 ∙ 10] ∙ |
|
|
||||||
|
|
|
|
|
Н |
( = Па |
|
|
|||
|
|
|
= м |
0 т |
|||||||
l(сn = ∙ −1.738 ∙ 10] , |
0, |
||||||||||
Теперь перейдем непосредственно к этапам решения: |
|
|
|
||||||||
Этап 1. Определяем матрицу жесткости. |
|
y |
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
с = v w f xf w f тgf |
|
||||||||
Здесь: |
|
|
|
|
fz |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
sin 2 f |
|
||
|
|
} cos |
f |
sin |
f |
… |
|||||
w f = |
|
|
|
2 |
|||||||
|
sin |
f |
cos f |
− sin22 f |
ƒ |
||||||
|
|
{− sin 2 f |
sin 2 f |
cos 2 f |
|||||||
w f − матрица поворота, связывающая напряжения в i-м монослое |
|||||||||||
|
xf |
|
|
‡f |
‡f ∙ f |
0 |
|
||||
|
= † ‡f |
∙ f |
|
|
‡f |
|
0 ˆ |
|
|||
|
|
|
|
|
0 |
|
|
0 |
|
|
|
xf − матрица, характеризующая упругость i-го монослоя |
|
|
|||||||||
|
|
|
f, f = 1 − f ∙ f |
|
|
||||||
|
|
|
‡ |
|
|
f, f |
|
|
|
||
В нашем случае матрица жесткости выглядит следующим образом: |
|
||||||||||
|
|
|
|
|
‰ |
‰ p |
0 |
|
|
||
|
|
с = †‰ p |
‰pp |
0 ˆ |
|
|
|||||
|
|
|
|
|
0 |
0 |
|
|
|
Мы можем определить модуль упругости в направлении оси ):
= ‰ − ‰ p
‰pp
И, соответственно, мы можем определить неизвестную величину :
= 0.003 ∙ ∙ − по условию задания
Этап 2. Определяем матрицу упругих податливостей.
•с = с ]
Этап 3. Определяем усредненный вектор деформаций. l1сn = •с l(сn
Этап 4. Определяем вектор деформаций в i-м монослое. l1fn = w f ] l1сn
w f − матрица поворота, связывающая деформации в i-м монослое
w f ] = w f т
Этап 5. Определяем вектор напряжений в i-м монослое. l(fn = xf l1fn
Поэтапное решение задачи представлено в системе MathCad prime 3.0.
Ввод исходных данных
схема армирования
|
|
|
|
n 12 |
i 1 n |
|
ϕ 30 |
ϕ −30 |
ϕ 90 |
ϕ |
0 |
|
|
1 |
4 |
7 |
|
10 |
|
|
ϕ 30 |
ϕ 90 |
ϕ 0 ϕ 0 |
δ 1 |
∑δ =1 |
||
2 |
5 |
8 |
|
11 |
i 12 |
i i |
ϕ −30 |
ϕ 90 |
ϕ 0 |
ϕ |
0 |
|
|
3 |
6 |
9 |
|
12 |
|
|
|
|
|
|
0.524 |
|
|
|
|
|
0.524 |
|
|
|
|
−0.524 |
||
|
|
|
|
−0.524 |
|
|
|
|
|
1.571 |
|
|
|
π |
|
1.571 |
|
ϕ |
ϕ |
= |
|
||
radi |
i |
180 |
|
1.571 |
|
|
|
|
|
0 |
|
|
|
|
|
0 |
|
|
|
|
|
0 |
|
|
|
|
|
0 |
|
|
|
|
|
0 |
|
E10 |
|
E1 |
=150.502 GPa |
|
1−ν ν |
||||
|
|
21 |
||
|
|
12 |
||
E20 |
|
E2 |
=8.027 GPa |
|
1−ν ν |
||||
|
|
21 |
||
|
|
12 |
l 0.7 m
A 9.31 10−4 m2
E1 150 GPa
E2 8 GPa
G12 5 GPa
ν12 0.25
ν21 ν12 E2 =0.013
E1