Добавил:
Upload Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

МДТТ дз1

.pdf
Скачиваний:
11
Добавлен:
23.03.2016
Размер:
608.99 Кб
Скачать

Государственное образовательное учреждение высшего профессионального образования

«Московский государственный технический университет имени Н.Э. Баумана»

(МГТУ им. Н.Э. Баумана)

ФАКУЛЬТЕТ

Специальное машиностроение

КАФЕДРА

Космические аппараты и ракеты-носители

Д О М А Ш Н Е Е З А Д А Н И Е № 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

 

‰qq

 

 

Мы можем определить модуль упругости в направлении оси ):

= ‰ − 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