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

Анализ оптимальных структур композитных стержней (120

..pdf
Скачиваний:
28
Добавлен:
15.11.2022
Размер:
405.68 Кб
Скачать

Copyright ОАО «ЦКБ «БИБКОМ» & ООО «Aгентство Kнига-Cервис»

1.3.Характеристики термического деформирования

Вдомашнем задании могут быть установлены требования к коэффициенту линейного термического расширения (КЛТР) в направлении оси стержня или к величине смещения торца стержня при заданном изменении температуры.

Характеристики термического деформирования многослойного

ортотропного материала определяются коэффициентами линейно-

го термического расширения αx и αy или коэффициентами термических напряжений βx и βy [1].

Для расчета термоупругих характеристик пакета необходимо знать КЛТР и технические константы жесткости каждого слоя:

α(i), α(i), E(i), E(i), G(i), ν(i)

E(i)ν(i)

= E(i)ν(i)

. Следует отме-

1

2

1

2

12

12

1

21

2

12

 

тить, что величины α1 для большинства современных углепластиков и органопластиков отрицательны (при нагреве они не расширяются, а сжимаются в продольном направлении). Величины α2 положительны для всех волокнистых композитов.

Сначала по известным значениям КЛТР слоя вычисляются его коэффициенты термических напряжений в естественной системе координат:

β(i)

= g(i)α(i)

+ g(i)α(i)

;

 

1

11

1

12

2

 

(9)

β(i)

= g(i)α(i)

+ g(i)α(i)

 

,

 

2

12

1

22

2

 

 

где коэффициенты матрицы жесткости слоя определяются согласно

(1).

Затем находятся коэффициенты термических напряжений многослойного пакета:

βx =

βy =

n

 

 

 

 

 

 

 

(10)

n

(i)

 

2

(i)

2

 

β1

cos

 

ϕi + β2

sin

ϕi

hi;

 

 

 

i=1

 

β1(i) sin2

ϕi + β2(i) cos2

ϕi h˜i.

i=1

 

 

 

 

 

 

 

 

Последний шаг — вычисление КЛТР многослойного материала:

αx = (βx − νxy βy)/Ex;

(11)

αy = (βy − νyxβx)/Ey,

 

 

11

Copyright ОАО «ЦКБ «БИБКОМ» & ООО «Aгентство Kнига-Cервис»

где технические константы жесткости пакета Ex, Ey, νxy, νyx определяются по формулам (1)—(4).

Продольное перемещение одного торца стержня относительно

другого uL при заданном изменении температуры

T пропорцио-

нально величине αx:

 

uL = αxL T.

(12)

На рис. 4 показан типичный график продольного КЛТР перекрестно армированного углепластика.

Рис. 4. Зависимость продольного КЛТР и смещения торца перекрестно армированного стержня от угла армирования

Следует помнить, что отрицательное значение КЛТР однонаправленного материала в направлении армирования α1 для большинства современных углепластиков приводит к тому, что начальный участок графика до диапазона 42...44соответствует отрицательным КЛТР: длина такого стержня при нагреве уменьшается,

апри охлаждении увеличивается.

1.4.Низшие частоты собственных колебаний

Вдомашнем задании предусмотрены три вида требований к динамическим характеристикам проектируемого стержня. Ограничения могут накладываться на следующие характеристики:

• низшую частоту продольных колебаний стержня f1(пр);

12

Copyright ОАО «ЦКБ «БИБКОМ» & ООО «Aгентство Kнига-Cервис»

низшую частоту крутильных колебаний стержня f1(кр);

низшую частоту изгибных колебаний стержня f1(изг).

Для расчета указанных величин необходимо знать технические константы жесткости материала каждого слоя: E1(i), E2(i), G(12i), ν(12i)

E1(i)ν(21i) = E2(i)ν(12i)

, а также плотности слоев ρi.

Сначала по формулам (1)—(4) рассчитываются технические константы упругости, затем определяются жесткостные характеристики стержня (5)—(8).

Масса тонкостенного трубчатого стержня вычисляется по фор-

муле

 

 

 

 

 

 

 

 

 

 

M = 2πRhLρ,

(13)

где

n

 

 

 

 

 

 

 

 

 

 

 

 

 

˜

 

 

 

 

 

 

 

ρ =

ρihi

(14)

 

 

 

 

 

i=1

 

 

 

 

 

 

 

— средняя плотность материала.

 

 

 

 

 

 

 

 

Низшая частота продольных колебаний стержня с закреплен-

ными торцами

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

f1(пр) = 2

 

 

ML

.

 

 

(15)

 

 

1

 

 

 

EF

 

Низшая частота крутильных колебаний стержня с закреплен-

ными торцами

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

.

 

f1(кр) = 2

 

 

2

πML xy

(16)

 

1

 

 

 

RhG

 

 

На рис. 5 приведены типичные графики зависимостей низших частот продольных и крутильных колебаний перекрестно армированных стержней из углепластика.

Низшая частота изгибных колебаний композитного стержня определяется на основе гипотезы Тимошенко, которая обеспечивает учет поперечного сдвига при изгибе стержня [3]. Расчетная формула для граничных условий, показанных на рис. 1, имеет вид

f1(изг) =

π

 

 

EI

 

 

.

(17)

2L

 

π2EI

 

 

 

 

ML 1 +

 

 

 

 

 

 

 

L2GF

 

 

 

 

 

 

 

 

 

 

 

 

13

Copyright ОАО «ЦКБ «БИБКОМ» & ООО «Aгентство Kнига-Cервис»

Рис. 5. Типичные зависимости низших частот продольных (а) и крутильных (б) колебаний перекрестно армированных стержней из углепластика от угла армирования

Рис. 6. Типичная зависимость низшей частоты изгибных колебаний перекрестно армированного стержня от угла армирования

Формула (17) справедлива для ненагруженного стержня. Типичная зависимость низшей частоты изгибных колебаний

перекрестно армированного стержня из углепластика показана на рис. 6. В соответствии с формулой (17), основная часть графика

14

Copyright ОАО «ЦКБ «БИБКОМ» & ООО «Aгентство Kнига-Cервис»

подобна зависимости квадратного корня из продольного модуля упругости, как и на графике рис. 5, а. Однако при малых углах ϕ оказывается существенным влияние низкой сдвиговой жесткости GF, что и приводит к уменьшению собственных частот и появлению максимума (на рис. 6 пунктиром показаны частоты, рассчитанные без учета сдвига по гипотезе Бернулли).

1.5.Характеристики демпфирования

Вразличных вариантах домашнего задания предусмотрен учет следующих характеристик демпфирования композитных стержней:

• коэффициенты диссипации при продольных колебаниях;

• коэффициенты диссипации при крутильных колебаниях;

• коэффициенты диссипации при изгибных колебаниях на низшей частоте;

• мощности диссипации при изгибных колебаниях на низшей частоте.

Воснове расчетных моделей для определения характеристик демпфирования композитных конструкций лежит энергетический подход [4].

Коэффициентом диссипации, или коэффициентом поглощения, любой механической системы называется отношение энергии, рассеянной за один период гармонических затухающих колебаний, к максимальной упругой энергии в этом периоде [5]. Если коэффициент диссипации не зависит от амплитуды, то логарифмический декремент колебаний постоянен и последовательные амплитуды составляют геометрическую прогрессию. Именно такая картина характерна обычно при малых колебаниях композитных конструкций. Для большинства современных композитов коэффициенты диссипации практически не зависят от частоты.

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

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

15

Copyright ОАО «ЦКБ «БИБКОМ» & ООО «Aгентство Kнига-Cервис»

соответствующих типах колебаний. Мощность диссипации равна произведению коэффициента диссипации и числа периодов колебаний за единицу времени, т. е. частоты собственных колебаний. Как следует из определения, коэффициент диссипации — безразмерная величина. Обычно его выражают в процентах, а мощность диссипации — в процентах в секунду.

Для расчета коэффициентов диссипации многослойных композитов и состоящих из них конструкций необходимо знать технические константы жесткости и коэффициенты диссипации каждого

слоя: E1(i), E2(i), G(12i), ν(12i) E1(i)ν(21i) = E2(i)ν(12i) , ψ(1i), ψ(2i), ψ(6i). Для

расчета мощностей диссипации к указанным величинам добавляются также плотности материалов слоев ρi.

Для расчета характеристик демпфирования сначала следует определить коэффициенты матриц жесткости отдельных слоев (1), а также коэффициенты матриц упругодиссипативных характеристик (УДХ) слоев в своих естественных системах координат [1]:

p(i)

= ψ(i)g(i)

;

 

11

1

11

 

 

p(i)

= ψ(i)g(i)

;

(18)

22

2

22

 

 

p(12i) = 1/2 ψ(1i) + ψ(2i) g12(i); p(66i) = ψ(6i)g66(i).

Затем вычисляются коэффициенты матриц УДХ слоев в системе координат конструкции:

p(i)

= p(i) cos4 ϕ

i

+ p(i) sin4

ϕ

i

+

 

 

xx

11

 

 

 

22

 

 

 

 

 

 

 

 

(i)

(i)

(i)

 

 

 

 

 

 

 

 

+

2p(i) + 4p(i)

sin2 ϕ cos2

ϕ ;

 

 

 

12

66

 

 

 

 

 

i

i

pxy(i) = p11 +(pi)22

44p66

sin24

ϕi cos2 ϕi +

p(i) = p(i) sin4 ϕ

+ p(i) cos4

ϕ +

 

(19)

 

+

p12

 

 

sin

ϕi + cos

 

ϕi ;

 

yy

11

 

i

 

22

 

 

i

 

 

 

 

 

 

(i)

(i)

(i)

 

 

 

 

 

 

 

 

+

p(i)

+ 4p(i)

sin2 ϕ

i

cos2 ϕ ;

 

 

12

 

66

 

 

 

 

 

 

i

pss(i) = p11 +(pi)22 − 2p12

sin2 ϕi cos2 2 ϕi +

 

+ p66

sin2 ϕi − cos2 ϕi .

 

16

Copyright ОАО «ЦКБ «БИБКОМ» & ООО «Aгентство Kнига-Cервис»

После этого могут быть определены коэффициенты диссипации пакета:

ψx =

1

 

n

pxx(i)

− 2pxy(i)νxy + pyy(i)νxy2

h˜i;

 

 

 

Ex

 

i=1

 

 

1

 

 

˜

(20)

ψs =

 

 

n

 

(i)

 

 

 

 

 

pss hi,

 

G

 

 

 

 

 

 

xy

i=1

 

 

 

 

 

 

 

 

 

 

где Ex, Gxy, νxy — технические константы жесткости многослойного ортотропного пакета, которые вычисляются согласно (1)—(4).

Коэффициент диссипации при продольных колебаниях композитного стержня вычисляется по первой формуле (20). Для расчета коэффициента диссипации при крутильных колебаниях следует использовать вторую формулу (20).

При изгибных колебаниях стержней, для которых существенно влияние поперечного сдвига, коэффициенты диссипации различаются для разных форм колебаний. Для собственных колебаний формы справедлива формула

 

ψxEI

L ϑ,x2 dx + ψsGF

L

(ϑ − w,x)2 dx

ψ(изг) =

L

L

 

 

, (21)

 

 

0

 

 

0

 

 

 

EI

 

ϑ,x2 dx + GF

(ϑ − w,x)2 dx

 

 

 

 

 

 

 

 

 

0

0

 

 

 

 

где индекс после запятой обозначает дифференцирование по данной координате; w(x) — прогиб стержня; ϑ(x) — угол поворота сечения в соответствии с гипотезой Тимошенко; EI и GF — жесткости стержня согласно (7) и (8).

Для граничных условий, показанных на рис. 1, зависимости ω(x) и ν(x) выражаются простыми тригонометрическими функциями, и формула (21) для коэффициента диссипации на низшей собственной частоте приобретает вид

 

 

ψ

+

ψ

π2EI

 

 

 

 

 

 

 

 

ψ(изг)

=

x

 

 

s L2GF

.

(22)

 

 

 

 

 

1

 

1 +

 

π2EI

 

 

 

 

 

L2GF

 

Мощность диссипации при изгибных колебаниях на низшей частоте определяется произведением величин ψ(1изг) из (22)

17

Copyright ОАО «ЦКБ «БИБКОМ» & ООО «Aгентство Kнига-Cервис»

и f

(изг) из (17):

 

 

 

 

1

= ψ(изг)f

(изг).

 

 

q(изг)

(23)

 

1

1

1

 

На рис. 7 показаны типичные зависимости коэффициентов диссипации перекрестно армированных стержней из углепластика при продольных и крутильных колебаниях, а на рис. 8 — соответствующие коэффициенты диссипации и мощности диссипации при изгибных колебаниях.

Рис. 7. Типичные зависимости коэффициентов диссипации при продольных (а) и крутильных (б) колебаниях перекрестно армированных стержней из углепластика от угла армирования

Рис. 8. Типичные зависимости коэффициентов диссипации (а) и мощностей диссипации (б) при изгибных колебаниях перекрестно армированных стержней из углепластика от угла армирования

18

Copyright ОАО «ЦКБ «БИБКОМ» & ООО «Aгентство Kнига-Cервис»

Хорошо заметное различие между графиками рис. 7, а и 8, а объясняется влиянием сдвига, особенно существенным при малых углах армирования. Вид графика рис. 8, б становится понятен исходя из графика рис. 8, а и графика рис. 6.

1.6. Несущая способность композитных стержней

1.6.1. Прочность

Как уже отмечалось, при выполнении домашнего задания прочность композита рассчитывается по первому разрушению. Таким образом, предельные напряжения по данному критерию соответствуют началу необратимого деформирования (разрушения) слоев композитного материала.

Алгоритмы расчета прочности при пропорциональном нагружении рассмотрены в [1]. Они сводятся к расчету предельного значения параметра нагрузки при заданных значениях начальных напряжений. Для расчета необходимо знать технические константы

(i) (i)

 

 

1

2

12

ν12

1 ν21 =

жесткости материала каждого слоя

 

E(i),E

(i)

,G(i),

(i)

E(i) (i)

= E2 ν12

и пять пределов прочности: при одноосном растяже-

нии F1+(i) и сжатии F1−(i) в направлении армирования, при одно-

осном растяжении F2+(i) и сжатии F2−(i) в поперечном направлении, а также при сдвиге в осях естественной системы координат

слоя F12(i).

Прочность слоев оценивается в соответствии с критерием максимальных напряжений для монослоя [1]:

F1−(i) σ1(i) F1+(i);

 

F2−(i) σ2(i) F2+(i);

(24)

12(i)| F12(i)

 

(величины F1−(i) и F2−(i) считаются отрицательными).

Вектор начальных напряжений при нагружении стержня осевой сжимающей силой N имеет вид

19

Copyright ОАО «ЦКБ «БИБКОМ» & ООО «Aгентство Kнига-Cервис»

 

(0)

 

 

 

−N/2πRh

 

 

 

σx(0)

 

=

 

0

.

(25)

σy

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(0)

 

 

 

0

 

 

 

τxy

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Если необратимое деформирование материала недопустимо, расчет прочности по первому разрушению сводится к такой последовательности операций:

1)рассчитываются по формулам (1)—(4) технические константы жесткости многослойного ортотропного пакета;

2)определяются начальные деформации пакета, соответствующие начальным напряжениям (25), т. е. единичному значению параметра нагрузки:

ε(0)

= σx(0)

− νxy σy(0)

x

 

 

Ex

 

 

 

ε(0)

= σy(0)

− νyxσx(0)

y

 

 

Ey

 

 

 

γ(0)

=

τxy(0)

;

xy

 

Gxy

 

 

 

 

;

;

(26)

3) рассчитываются деформации отдельных слоев при единичном значении параметра нагрузки:

ε(0)

= ε(0) cos2

ϕ + ε(0) sin2

ϕ + γ(0) cos ϕ sin ϕ ;

 

1i

x

i

y

i

xy

i

i

 

ε2(0)i

= εx(0) sin2 ϕi + εy(0) cos2

ϕi − γxy(0) cos ϕi sin ϕi;

(27)

γ12(0)i = − εx(0) − εy(0)

sin2ϕi + γxy(0) cos2ϕi.

 

 

Поскольку в векторе начальных напряжений τ(0)xy = 0, то напряжения в частях слоя, армированных под углами +ϕi и −ϕi, одинаковы;

4) определяются напряжения в слоях пакета, соответствующие единичному значению параметра нагрузки (при заданной осевой силе):

20

Соседние файлы в предмете [НЕСОРТИРОВАННОЕ]