
Анализ оптимальных структур композитных стержней (120
..pdf
Copyright ОАО «ЦКБ «БИБКОМ» & ООО «Aгентство Kнига-Cервис»
σ(0) |
= g(i)ε(0) |
+ g(i)ε(0) |
; |
|
1i |
11 1i |
12 2i |
|
|
σ2(0)i |
= g12(i)ε1(0)i |
+ g22(i)ε2(0)i |
; |
(28) |
τ(0) |
= g(i)γ(0), |
|
|
|
12i |
66 12i |
|
|
|
где коэффициенты матрицы жесткости каждого слоя в его естественной системе координат находятся по формулам (1);
5) используется критерий прочности для монослоя (24). Поскольку напряжения в слоях, как и средние напряжения в пакете, пропорциональны значению параметра нагрузки, предельное значение этого параметра показывает, во сколько раз необходимо увеличить напряжения (28) для того, чтобы в одном из слоев выполнилось в виде равенства одно из условий (24):
|
i |
|
|
|
(i) |
|
|
(i) |
|
|
|
|
|
||
|
|
|
σ1 |
|
|
|
σ2 |
|
(i) |
|
|
||||
|
|
|
|
|
σ1i |
|
|
σ2i |
|
|
|
|
|||
|
|
|
|
|
|
|
τ12i |
|
|||||||
Pпр = min min |
|
(0) |
, |
(0) |
, |
|
F12 |
|
, |
(29) |
|||||
(0) |
|
||||||||||||||
где |
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
= |
|
|
|
|
|
|
|
|
|
||||
|
σ1 |
(i)+ |
|
|
|
(0) |
|
|
|
|
|
||||
|
|
|
|
F(i)−, |
если σ(0) |
< 0; |
|
|
|
||||||
|
F1 |
|
, если σ1i |
> 0; |
|
|
|
||||||||
|
(i) |
|
|
1 |
|
|
|
|
1i |
|
|
|
|
|
|
|
(i) |
|
|
F2(i)−, если σ2(0)i |
< 0; |
|
|
(30) |
|||||||
|
σ2 |
= |
|
(i)+ |
, |
|
|
(0) |
> 0. |
|
|
|
|||
|
|
|
F2 |
|
если σ2i |
|
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Если в формулу (25) подставить расчетное значение силы N, то величина Pпр, определенная согласно (29), имеет смысл коэффициента запаса прочности. Требование сохранения прочности при заданной нагрузке N трансформируется в требование
Pпр 1. |
(31) |
Типичный график прочности перекрестно армированного стержня приведен на рис. 9, а. График состоит из четырех ясно различимых участков, соответствующих различным механизмам разрушения, согласно критерию (24):
21

Copyright ОАО «ЦКБ «БИБКОМ» & ООО «Aгентство Kнига-Cервис»
Рис. 9. Типичные зависимости прочности (а) и общей устойчивости (б) перекрестно армированных стержней из углепластика от угла армирования
•при малых углах армирования разрушение происходит вследствие достижения предела прочности при сжатии в направлении армирования;
•при увеличении угла армирования механизмом первого разрушения является растяжение поперек волокон (для некоторых материалов этот участок может отсутствовать);
•при дальнейшем увеличении угла армирования первое разрушение происходит от сдвига;
•заключительная часть графика соответствует первому разрушению от сжатия поперек волокон.
1.6.2. Устойчивость стержня
Устойчивость стержня, нагруженного осевой сжимающей силой (общая устойчивость), оценивается с использованием гипотезы Тимошенко, которая обеспечивает учет поперечного сдвига при закритическом изгибе стержня [3]. Формула для расчета критической нагрузки Nкр при граничных условиях, показанных на рис. 1,
имеет вид |
π2EI |
|
|
|
|||
Nкр = |
|
|
, |
(32) |
|||
L2 |
1 + |
π2EI |
|
||||
|
|
|
|||||
|
L2GF |
|
|
где EI и GF — жесткости стержня согласно (7) и (8).
22
Copyright ОАО «ЦКБ «БИБКОМ» & ООО «Aгентство Kнига-Cервис»
Требование сохранения несущей способности при заданной нагрузке N сводится к условию
Nкр/N 1. |
(33) |
Типичный график запаса общей устойчивости (см. (32), (33)), для перекрестно армированного стержня приведен на рис. 9, б.
1.6.3. Устойчивость оболочки
Тонкостенный трубчатый стержень способен терять устойчивость не только по стержневым формам (когда ось его искривляется, а поперечные сечения остаются недеформированными), но и как цилиндрическая оболочка (т. е. с сохранением прямолинейной оси и искажением геометрии сечений). Такие формы относятся к
местной устойчивости конструкции.
Для однородных изотропных оболочек оценка критической нагрузки ведется по формулам, выведенным из рассмотрения осесимметричных форм потери устойчивости. Для таких оболочек критические нагрузки осесимметричной и неосесимметричной потери устойчивости совпадают [3]. Однако для анизотропных конструкций осесимметричная и неосесимметричная потери устойчивости наступают при разных нагрузках, причем в зависимости от структуры армирования наиболее опасной может быть как одна, так и другая потеря. Нижеследующие формулы справедливы для оболочек, торцы которых подкреплены шпангоутами.
Для расчета на устойчивость композитной оболочки необходимо знать технические константы жесткости материала каждого
слоя: E1(i), E2(i), G(12i), ν(12i) E1(i)ν(21i) = E2(i)ν(12i) . Предваритель-
но должны быть рассчитаны по формулам (1), (2) коэффициенты матрицы жесткости каждого слоя и определены мембранные, смешанные и изгибные жесткости стенки оболочки:
n |
|
|
Bxx = gxx(i)hi; |
|
|
i=1 |
(34) |
|
n |
||
|
||
Bxy = gxy(i)hi; |
|
|
i=1 |
|
|
|
23 |
Copyright ОАО «ЦКБ «БИБКОМ» & ООО «Aгентство Kнига-Cервис»
|
|
|
|
|
|
|
|
|
|
|
|
|
|
n |
|
|
|
|
|
|
Byy = |
|
gyy(i)hi; |
|
|
||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
i=1 |
|
|
|
|
|
(34) |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
n |
|
|
|
|
|
|
Bss = |
|
gss(i)hi; |
|
|
||||||||||||||||
|
|
|
||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
i=1 |
|
|
|
|
|
|
Cxx = |
1 |
|
|
|
n |
(i) |
zi2 − zi2−1 |
; |
||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
i=1 gxx |
||||||||
|
|
2 |
||||||||||||||||||
|
1 |
|
|
|
|
(i) |
|
2 |
|
2 |
|
|
||||||||
|
1 |
|
|
|
|
n |
(i) |
|
2 |
|
2 |
|
|
|||||||
|
|
|
|
|
|
|
|
|||||||||||||
Cxy = |
|
|
2 |
|
|
|
i=1 gxy |
zi |
− zi−1 |
; |
||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
n |
|
|
|
|
|
|
(35) |
|
1 |
|
|
|
|
|
|
|
|
(i) |
|
2 |
|
2 |
|
|
||||
Cyy = |
|
|
2 |
|
|
|
i=1 gyy |
zi |
− zi−1 |
; |
||||||||||
|
1 |
n |
(i) |
|
|
|
|
|||||||||||||
|
|
|
|
|
||||||||||||||||
Css = |
|
2 |
|
|
|
|
i=1 gss zi − zi−1 ; |
|||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
n |
|
|
zi3 − zi3−1 ; |
||||
Dxx = |
|
|
|
|
i=1 gxx |
|||||||||||||||
3 |
||||||||||||||||||||
|
1 |
|
(i) |
3 |
|
3 |
|
|
||||||||||||
|
1 |
|
n |
(i) |
|
3 |
|
3 |
|
|
||||||||||
|
|
|
|
|
||||||||||||||||
Dxy = |
|
|
|
3 |
i=1 gxy |
zi |
− zi−1 ; |
|||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
n |
|
|
|
|
|
|
(36) |
|
1 |
|
|
|
|
|
|
(i) |
|
3 |
|
3 |
|
|
||||||
Dyy = |
|
|
3 |
i=1 gyy |
zi |
− zi−1 |
; |
|||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
n |
|
|
|
|
|
|
|
Dss = |
|
i=1 gss zi − zi−1 . |
||||||||||||||||||
3 |
||||||||||||||||||||
В выражениях (35) и |
(36) величины |
z |
|
|||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
i−1 и zi суть координа- |
ты внутренней и внешней поверхностей i-го слоя, отсчитываемые от координатной поверхности оболочки [2]. Положение координатной поверхности зависит от структуры многослойного пакета. Если эта структура симметрична относительно срединной поверхности (как, например, рассматриваемая в домашнем задании перекрестно армированная структура), то координатная поверхность выбирается совпадающей с серединной поверхностью оболочки. В этом случае все смешанные жесткости (35) оказываются равными нулю. Для перекрестно армированной структуры зависимости (34) и (36) упрощаются:
Bxx = gxxh; |
|
Bxy = gxyh; |
Byy = gyyh; Bss = gssh; |
(37) |
||||||
Dxx = |
gxxh3 |
; |
Dxy = |
gxyh3 |
; Dyy = |
gyyh3 |
; |
Dss = |
gssh3 |
. (38) |
|
|
|
|
|||||||
12 |
|
12 |
12 |
|
12 |
|
24

Copyright ОАО «ЦКБ «БИБКОМ» & ООО «Aгентство Kнига-Cервис»
Критическая нагрузка осесимметричной потери устойчивости при осевом сжатии описывается простой формулой [3]
N(ос) = 4π |
|
|
|
(39) |
|
D |
E h, |
||||
y |
|
|
xx y |
|
|
|
кр |
|
|
в которой величина E определяется согласно (4).
Критическое значение числа полуволн в осевом направлении, соответствующего осесимметричной потери устойчивости, может
быть оценено по формуле |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
m0 |
= π |
R2 |
(BxxDxx− |
C2 |
). |
(40) |
|||
|
L |
4 |
B |
B |
B2 |
|
|
|
|
|
|
|
xx yy |
xy |
|
|
|
− xx
Для симметричных относительно срединной поверхности (в частности, перекрестно армированных) структур выражение (40)
упрощается: |
|
|
|
|
|
|
m0 |
= π 4 |
|
|
|
(41) |
|
|
DxxR2 . |
|||||
|
|
L |
|
Eyh |
|
При вычислениях значениe величины m0 должно быть округлено до ближайшего целого числа.
Требование сохранения несущей способности при заданной нагрузке N при учете только осесимметричных форм потери устойчивости сводится к условию
kNкр(ос)/N 1, |
(42) |
где k — коэффициент устойчивости, значениe которого задается в домашнем задании.
При учете неосесимметричных форм потери устойчивости заранее указать наиболее опасную форму невозможно. В этом случае необходимо составить алгоритм, осуществляющий перебор целочисленных значений числа полуволн в осевом направлении (m = 1,2,3,...) и числа волн по окружности (n = 2, 3, 4, . . . ).
Для каждой пары параметров волнообразования m и n вычисляется собственное значение параметра нагрузки Pmn по формуле
|
a11 |
a12 |
a13 |
|
|
|
Det |
a21 |
a22 |
a23 |
= 0, |
(43) |
|
|
a31 |
a32 |
a33 −b33 Pmn |
|
|
25 |
|
|
|
|
|
|
Copyright ОАО «ЦКБ «БИБКОМ» & ООО «Aгентство Kнига-Cервис»
где |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
a |
11 |
= λ2 |
B |
+ λ2 |
B |
; |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
|
|
m |
|
|
xx |
|
|
n |
|
ss |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
a12 = a21 = λmλn (Bxy + Bss); |
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||
a |
22 |
= λ2 |
B |
+ λ2 |
|
B |
ss |
; |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
n |
|
|
yy |
|
|
m |
|
|
|
|
|
|
|
|
|
|
|
|
R |
|
; |
(44) |
|||
a13 |
= a31 |
|
= −λm |
λm2 Cxx + λn2 (Cxy + 2Css) + |
|
||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Bxy |
|
|
|
a23 |
= a32 |
|
= −λn |
λn2Cyy + λm2 (Cxy + 2Css) + |
R |
|
; |
|
|||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Byy |
|
|
|
|
a |
33 |
= |
Byy |
+ λ4 D |
xx |
+ 2λ2 |
λ2D + λ4D |
yy |
+ 4λ2 |
λ2D ; |
|||||||||||||||||
|
|
||||||||||||||||||||||||||
|
|
R2 |
|
m |
|
|
m |
n |
xy |
n |
|
|
m |
|
n |
ss |
|||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
b |
|
= λ2 |
|
|
N |
|
+ λ2pR; |
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||
33 |
|
m 2πR |
|
n |
|
|
mπ |
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
|
|
|
|
|
|
|
|
|
|
|
λm |
= |
; |
λn = |
|
n |
. |
|
|
|
|
|
(45) |
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
L |
|
|
R |
|
|
|
|
|
|
|
Величина p в последнем выражении (44) представляет собой значение давления, действующего на боковую поверхность оболочки (в домашнем задании p = 0). В формулы (44) подставляются положительные значения сжимающей силы N и внешнего давления p.
Формула (43) получена на основе теории пологих оболочек [3] и справедлива в тех случаях, когда число волн по окружности достаточно велико (это условие обычно выполняется для тонкостенных оболочек).
Перебор параметров волнообразования должен производиться
восевом направлении по крайней мере до значения m0 (см. (40) или (41)). Что касается числа n, то обычно достаточно 20—30 волн
вокружном направлении.
Предельное значение параметра нагрузки по критерию местной устойчивости (неосесимметричные формы)
P(неос) |
= |
k |
min |
{ |
P |
mn} |
, |
(46) |
пред |
|
· m,n |
|
|
|
где k — коэффициент устойчивости, значениe которого задается в домашнем задании.
По смыслу величина Pпред(неос) из (46) представляет собой запас устойчивости по отношению к действующим нагрузкам. Требова-
26

Copyright ОАО «ЦКБ «БИБКОМ» & ООО «Aгентство Kнига-Cервис»
ние сохранения несущей способности при учете неосесимметричных форм потери устойчивости сводится к условию
Pпред(неос) 1. |
(47) |
Типичные графики устойчивости перекрестно армированных стержней приведены на рис. 10. График рис. 10, а показывает только осесимметричные формы, что соответствует по формуле (39). На графике рис. 10, б показаны все формы: как осесимметричные, так и неосесимметричные (учет осесимметричных форм может быть проведен и непосредственно в рассмотренном выше алгоритме, если начинать перебор значений по параметру n не со значения 2, а с 0). Видно, что середина графика рис. 10, б соответствует осесимметричной потере устойчивости, тогда как при уменьшении или увеличении угла армирования конструкция теряет устойчивость по неосесимметричным формам.
Рис. 10. Типичные зависимости осесимметричной устойчивости оболочки (а) и устойчивости с учетом всех форм (б) для перекрестно армированных стержней из углепластика от угла армирования
Если при построении графика вместо плавной кривой, аналогичной кривой на рис. 10, б получаются резкие скачки, то это означает, что при переборе было учтено недостаточное число волн по оси или в окружном направлении.
27
Copyright ОАО «ЦКБ «БИБКОМ» & ООО «Aгентство Kнига-Cервис»
2. МЕТОДИКА ВЫПОЛНЕНИЯ ДОМАШНЕГО ЗАДАНИЯ
По итогам решения домашнего задания оформляется отчет, который включает в себя номер варианта и исходные данные, а также необходимые формулы, графики и результаты вычислений в соответствии с алгоритмом, приведенным ниже. Пример оформления приведен в разд. 3.
Оформленное домашнее задание представляется на защиту, во время которой студент должен быть готов ответить на вопросы по материалу, приведенному в разд. 1.
Исходными данными задачи являются:
•радиус тонкостенного стержня R;
•длина стержня L;
•толщина стенки стержня h;
•характеристики однонаправленного композита, из которого
состоит стержень:
◦модуль упругости однонаправленного материала в направлении армирования E1;
◦модуль упругости однонаправленного материала в поперечном направлении E2;
◦модуль сдвига однонаправленного материала G12;
◦коэффициент Пуассона однонаправленного материала ν12;
◦КЛТР однонаправленного материала в направлении армирования α1;
◦КЛТР однонаправленного материала в поперечном направлении α2;
◦коэффициент поглощения энергии однонаправленного материала в направлении армирования ψ1;
◦коэффициент поглощения энергии однонаправленного материала в поперечном направлении ψ2;
◦коэффициент поглощения энергии однонаправленного материала при сдвиге ψ6;
◦плотность однонаправленного материала ρ;
◦предел прочности однонаправленного материала при растяжении вдоль волокон F1+;
◦предел прочности однонаправленного материала при сжатии вдоль волокон F1−;
◦предел прочности однонаправленного материала при растяжении в поперечном направлении F2+;
28
Copyright ОАО «ЦКБ «БИБКОМ» & ООО «Aгентство Kнига-Cервис»
◦предел прочности однонаправленного материала при сжатии в поперечном направлении F2−;
◦предел прочности однонаправленного материала при сдвиге
F12;
•значение коэффициента местной (оболочечной) устойчивости k;
•наименьшее допустимое значение одной из трех жесткостных характеристик стержня (EF, GIp или EI);
•наибольшее допустимое значение одной из характеристик термического деформирования (продольный КЛТР αx или продольное перемещение торцов при нагреве на заданную величину
T− uL);
•наименьшее допустимое значение одной из трех частот соб-
ственных колебаний стержня f1(пр), f1(кр) или f1(изг) ; |
•наименьшее допустимое значение одной из четырех характеристик демпфирования стержня ψx, ψs, ψ(1изг) или q1(изг) ;
•значениe осевой сжимающей нагрузки, при которой должна быть обеспечена несущая способность конструкции.
Кроме того, задано, что стержень должен быть образован намоткой одного семейства перекрестно армированных слоев (т. е. имеет перекрестно армированную структуру) и является тонкостенным. При нагружении стержня необратимое деформирование не допускается. При расчете местной устойчивости используется теория пологих оболочек.
Необходимо определить диапазон допустимых значений угла армирования, исходя из требований, предъявляемых к конструкции.
При выполнении домашнего задания рекомендуется такая последовательность действий:
1) построить графики всех характеристик, к которым предъявляются требования. Для построения использовать зависимости
(1)—(46) и привести их в отчете. Для контроля правильности построения использовать рис. 2—10. Рекомендуется строить графики с использованием системы MathCAD или любой другой системы программирования. Точность графиков должна быть достаточной для корректного определения структурных параметров (в системе
29
Copyright ОАО «ЦКБ «БИБКОМ» & ООО «Aгентство Kнига-Cервис»
MathCAD рекомендуется использовать не менее 40—50 точек для построения каждой линии);
2)отметить на графиках требуемые значения каждой характеристики и соответствующие им структурные параметры (углы армирования). При необходимости уточнить значения структурных параметров численным расчетом (необходимая точность — не менее 0,5◦);
3)записать допустимые диапазоны значений угла армирования, исходя из каждого требования к конструкции;
4)определить допустимый диапазон значений угла армирования с учетом всей совокупности предъявляемых к конструкции требований.
3.ПРИМЕР ВЫПОЛНЕНИЯ ДОМАШНЕГО ЗАДАНИЯ
Исходные данные. Проектируется цилиндрический трубчатый стержень из углепластика, длиной L = 1,8 м и радиусом R = 75 мм. Стержень должен быть образован намоткой одного семейства перекрестно армированных слоев и иметь толщину h = 1,6 мм. Определить диапазон допустимых значений угла намотки, исходя из предъявляемых к конструкции требований:
•продольная жесткость стержня — не менее 20 МН;
•продольный КЛТР не более 5 · 10−6 K−1;
•низшая частота крутильных колебаний не менее 800 Гц;
•коэффициент диссипации при изгибных колебаниях на низшей собственной частоте — не менее 1,8 %;
•необходимость обеспечения несущей способности при нагрузке осевого сжатия — 70 кН.
При расчетах надо использовать следующие характеристики однонаправленного углепластика:
•модуль упругости однонаправленного материала в направлении армирования E1 = 140 ГПа;
•модуль упругости однонаправленного материала в поперечном направлении E2 = 9,6 ГПа;
•модуль сдвига однонаправленного материала G12 = 4,6 ГПа;
•коэффициент Пуассона однонаправленного материала ν12 =
=0,3;
30