книги / Прикладные задачи устойчивости многослойных композитных оболочек
..pdfП р и л о ж е н и е А
РАСЧЁТ ЖЕСТКОСТЕЙ НЕОДНОРОДНЫХ ПО
ТОЛЩИНЕ МНОГОСЛОЙНЫХ И ТРЁХСЛОЙНЫХ
КОМПОЗИТНЫХ КОНСТРУКЦИЙ
А.1. Неоднородные по толщине многослойные оболочки
А. 1.1. Структура и жёсткости многослойного ортотропного пакета. Рассматривается в общем случае неоднородный, в частно сти, многослойный пакет общей толщиной h, состоящий из N слоёв (рис. А.1). Нумерация слоёв производится, начиная с внутренней по верхности оболочки. Каждый к-й слой толщиной h^ (к = 1,2,... , N) имеет однородную структуру, его нейтралвная (срединная) поверхности отстоит от внутренней поверхности оболочки на расстояние а,(к\
Полагается, что все слои являются ортотропными, их оси ортотропии совпадают с линиями главнв1х кривизн оболочки. В этом случае характеристики упругости каждого к-то слоя в тангенциаль ной плоскости ху определяются четырьмя модулями упругости Е ^
(г = 1,2, 3,4): |
Е |
, Е ^ |
— модули упругости |
в направлениях осей |
||
х (i = 1) и у |
(i = 2); Е ^ = G |
— модуль сдвига в тангенциальной |
||||
плоскости (i = 3); |
Е ^ |
= i / ^ E ^ |
= v ^ E ^ |
— |
модули, связанные |
|
с влиянием коэффициентов Пуассона щ и U2 (г = |
4). |
Сопротивление слоёв поперечным сдвигам характеризуется величи
нами G ( i = 1,2), которые представляют собой модули сдвига А;-го слоя в плоскостях xz (i = 1) и yz (i = 2).
Здесь и в дальнейшем нижний индекс г относится к направлению анизотропии, а верхний индекс к — к номеру слоя.
Сопротивление ортотропного к-ro слоя в рамках модели прямой
линии определяется жесткостями в [ к^ |
, |
на |
растяжение-сжатие |
в двух направлениях, жёсткостью В ^ |
= |
В ^ на |
сдвиги в тангенци |
альной плоскости, изгибными жесткостями D\k\ |
в двух направ |
||
лениях, крутильной жёсткостью D^ |
|
, а также влиянием коэф |
фициентов Пуассона. Кроме того, сопротивление поперечным сдвигам в плоскостях xz, yz характеризуется сдвиговыми жесткостями к \ к\
К ^ . Перечисленные жёсткости находятся с помощью известных со отношений [31]:
А.1. Неоднородные no толщине многослойные оболочки |
223 |
оболочки (рис. АЛ). Тогда для расчёта жесткостей получаются следу ющие зависимости:
Д |
Ei (z) dz; |
С\ = |
Ei (z)z dz; |
|
z0- h |
|
z0- h |
Di |
Ei (z) z2dz; |
Ei (z) |
Ei |
= |
|||
z0- h |
|
1 - Ui (z) v2 (z) ’ |
|
|
|
||
E 3 (z) |
= G (z); E 4(z) = vi(z)E2(z) = v2(z)E\(z); |
||
|
i = |
1, 2,3,4. |
(A.2) |
Здесь Bi — жёсткости многослойного пакета на растяжение-сжатие в двух направлениях (i = 1, 2), на сдвиги в тангенциальной плоскости (г = 3) и жёсткости на растяжение-сжатие, связанные с коэффициен тами Пуассона (* = 4); Di — некоторые изгибные жёсткости в двух направлениях (i = 1, 2), крутильная жёсткость (г = 3), изгибные жёст кости, связанные с коэффициентами Пуассона (г = 4), вычисленные относительно координатной поверхности ZQ. Величины Q являются статическими моментами соответствующих жесткостей Д относитель но оси, проходящей на уровне ZQ\ э т и величины характеризуют взаим ное влияние изгибных и мембранных компонент НДС в соотношениях закона Гука [32].
Как следует из (А.2), жёсткости Д многослойного пакета на растяжение-сжатие и сдвиги в тангенциальной плоскости не зависят от величины ZQ, т. е. от выбора координатной поверхности. В то же время статические моменты Ci жесткостей и изгибные и крутильные жёсткости Di существенно зависят от величины zo.
Можно показать, что для многослойного пакета с ортотропными слоями имеется четыре положения координатной поверхности Zi (г = 1, 2,3,4), для которых соответствующие статические моменты Ci обращаются в нуль, а изгибные и крутильные жёсткости прини мают минимальные значения. Именно минимальными жесткостями
определяется |
сопротивление многослойной конструкции изгибам |
|
и кручениям [30]. |
|
|
Принимая |
это |
во внимание, из (А.2) получаем следующие за |
висимости для определения минимальных жесткостей многослойного ортотропного пакета:
h Z i
Bi |
Ei (z) dz; |
Di |
Ei (z) z2dz; |
Ci = 0; |
|
о |
{ h - Z i ) |
|
|
|
1 |
h |
|
|
|
Ei (z) z dz; |
1, 2,3,4. |
(A.3) |
|
|
Z i |
|||
|
Д |
о |
|
|
|
|
|
|
224 Прил. А. Неоднородные по толщине композитные конструкции
Жёсткости Ki многослойного пакета на поперечные сдвиги в плос костях xz (г = 1) и yz (i = 2) определяются зависимостями [31]:
_ 1_ |
Cf(z) |
|
Zi |
dz; Ci{z) |
Ei (z) z dz\ |
||
Ж |
Gi(z) |
|
Z |
|
|
|
|
|
i |
= 1 , 2 . |
(A.4) |
Учитывая, что в многослойных конструкциях модули упругости изменяются по толщине ступенчато, целесообразно (при N > 3) в при ведённых соотношениях интегрирование заменить суммированием:
Bi |
= |
£ |
# |
|
|
|
|
|
k= 1 |
|
|
|
|
Ci |
= |
N |
В гк){ (z0 - |
а ^ ) |
= BiZo - Z B ^ a W - , |
|
£ |
||||||
|
|
k= 1 |
|
|
k= |
1 |
D i |
= |
D CA |
+ D *-, |
|
|
(A.5) |
|
|
|
||||
Df |
|
|
D* = £ |
B\k) (z0 - |
a(fc))2; |
|
|
|
i = 1,2, 3,4. |
|
|
|
|
Величина |
|
D f l |
представляет собой сумму собственных изгибно- |
|||
крутильных жесткостей |
слоёв, составляющих пакет. Величина D * |
|||||
— сумма моментов жесткостей В ^ |
слоёв относительно выбранной оси, |
|||||
отстоящей от внутренней поверхности на расстояние ZQ. |
Из (А.5) видно, что жёсткости Bi не зависят от выбора величины ZQ, т. е. от выбора координатной поверхности. В то же время стати ческие моменты Ci и изгибно-крутильные жёсткости Di существенно содержат величину ZQ. Для каждого номера i (i = 1,2,3,4) можно найти такое положение zo = z* координатной поверхности, при кото ром соответствующие статические моменты равны нулю. Приравнивая в (А.5) величину C i нулю, получим:
1 N |
N |
п ( к ) |
|
zi = f r ^ 2 |
a{k)B*k) = JE |
a{k)- B - ’ i = 1 >2’3’4- |
(А-6) |
Из (А.5) следует, что величина С* равна нулю, если ZQ = z». Легко показать, что если изгибно-крутильные жёсткости D i (i = 1, 2,3,4) вы числять относительно осей, проходящих на уровне zit то эти жёсткости принимают минимальные значения. Минимальные жёсткости не зави сят от выбора координатной поверхности, являются собственными чис лами многослойного пакета и определяют сопротивление многослойной конструкции изгибам и кручению.
А.1. Неоднородные no толщине многослойные оболочки |
225 |
Подставляя в (А.5) вместо ZQ соответствующие величины z* из (А.6), можно найти минимальные изгибные и крутильную жёсткости:
А |
= |
А ?л + D*; |
|
|
|
|
А СЛ= |
Е Д |
М : |
D * = |
j r B p } ( Z i - |
a (fc) ) 2; |
|
|
|
к=1 |
|
|
k=1 |
(A.7) |
A |
= |
N |
и л |
A = 0; |
i = 1,2,3,4. |
|
£ A |
( > |
|
k= 1
Заменяя в (A.4) интегрирование послойным суммированием, мож но получить зависимости для определения жесткостей многослойного пакета на поперечные сдвиги:
1 |
1 |
|
( c j k)hW) |
■ |
к ^ |
б * |
hW ■ |
||
|
А? |
£ |
|
к[к) |
|||||
K i |
|
’ |
-*4 |
v |
’ |
||||
Cf> = £ |
В\з) (Zi - aW) = |
|
+ (zt - a(fc)) |
|
(A.8) |
||||
|
B\(*). |
||||||||
|
J=1 |
|
|
|
|
|
|
|
|
c (0) = |
C |
( N) = |
Q. |
1, |
2. |
|
|
|
|
|
|
|
|
|
|
|
|
Вслучае однородных оболочек все четыре значения Z* сливаются
водно: z\ = Z2 = Z3 = Z4 = h/2. При этом срединная поверхность одновременно является и нейтральной, т.е. вычисленные относительно неё статические моменты равны нулю (С\ = С2 = С3 = Сц = О).
Точки в поперечном сечении многослойного пакета, отстоящие от внутренней поверхности на расстояниях Z* (г = 1,2,3,4) и вычислен ные по формулам (А.З) или (А.6), принято называть центрами жестко стей. Четыре центра жесткостей соответствуют четырём независимым характеристикам упругости ортотропных пластин и оболочек: Е\, А .
G, E \ V2 = E 2V \.
Таким образом, в отличие от однородных конструкций, многослой ные оболочки, составленные из ортотропных слоёв, в силу незави симости четырёх постоянных упругости ортотропного слоя и произ вольности расположения слоёв в пакете в общем случае не имеют единой нейтральной поверхности. Вместо этого имеется четыре опре делённых поверхности z* (г = 1, 2,3,4), относительно которых соот ветствующие статические моменты А равны нулю, а соответствующие изгибно-крутильные жёсткости А принимают минимальные значения. По аналогии с однородными оболочками четыре поверхности, соответ ствующие уровням Zj, можно условно назвать «нейтральными» [32].
Единая нейтральная поверхность существует в многослойной обо
лочке только в том случае, когда правая |
часть (А.6) |
не |
зависит |
|
от индекса г. В частности, |
это условие |
выполняется |
для |
оболо |
чек, составленных из слоёв, |
одинаковых по механическим |
свойствам, |
8 С.Н. Сухинин
226 Прил. А. Неоднородные по толщине композитные конструкции
но разной толщины. Оно выполняется также, если многослойный па кет симметричен по толщине относительно срединной поверхности. В обоих этих случаях нейтральная поверхность совпадает со средин ной. Из (А.6) следует, что многослойная ортотропная оболочка имеет единую нейтральную поверхность, если выполняется условие:
D (&) |
D W |
n (*0 |
r>(k) |
|
|
|
n l |
_ d 2 |
_ д 3 |
_ n 4 |
k = 1,2,3,... |
,N. |
|
B \ |
B 2 |
B 3 |
B 4 |
|||
|
|
Можно попытаться найти и другие частные случаи, когда много слойный пакет, составленный из ортотропных слоёв, имеет единую нейтральную поверхность.
Как известно, во многих случаях коэффициенты Пуассона оказыва ют незначительное влияние на результаты расчёта НДС, устойчивости и колебаний тонкостенных конструкций. Поэтому при решении инже нерных задач обычно вводятся приведённые значения коэффициентов Пуассона, одинаковые для всего пакета. Приведенные коэффициенты Пуассона получаются из зависимостей (А.5) при г = 4:
|
1 |
SL |
„ ш( к ) ((к\) |
1 |
N |
г > ( к ) ( к ) |
; |
v\ |
В ^ |
" = -б - |
Е |
в \ щ. |
|||
|
|
fc=l |
|
В 2 |
К— 1 |
|
|
^2 = |
4 - |
Е |
в \ к)4 к) = |
- 5 - jt |
B (2k)v\k) |
(А.9) |
|
|
|||||||
|
|
к= 1 |
|
|
к — 1 |
|
|
В 11/2 = В21/\.
Если за коэффициенты Пуассона принять величины (А.9), то вы полняются следующие зависимости, получаемые из (А.6):
|
, |
N |
N |
|
N |
z 1= |
24 ж |
£ |
BS4«W = % £ в ^ |
(к) = ж |
£ в ^ т - |
|
|
к= 1 |
fc=l |
ВА |
к= 1 |
|
|
|
Z3 = T r E B3fc)«(fc) |
|
(АЛО) |
|
|
|
Вз ы |
|
|
Отсюда следует, что если коэффициенты Пуассона одинаковы для всего пакета, то число центров жесткостей z* (число «нейтральных» поверхностей) снижается до двух. Один из них соответствует изгибу, другой — кручению многослойного пакета.
А Л .2. Расчёт характеристик упругости многослойных орто тропных пакетов с произвольным армированием. Структура мно гослойного композитного пакета представляет собой в общем случае произвольно расположенные слои следующих четырёх видов:
—слой, образованный армированием в некотором 1-м направлении
суглом укладки </? = 0° (например, в осевом направлении цилиндриче ской оболочки) — рис. А.2;
228 Прил, А. Неоднородные по толщине композитные конструкции
Т а б л и ц а A.l
О
о
о
|
н |
|
|
|
н |
|
|
|
К |
|
|
|
S |
|
|
9- |
X |
|
|
с* |
|
|
|
t=3 |
m |
|
|
о |
|
|
|
О |
Он |
|
|
sH |
s |
|
|
^ |
S |
|
|
|
Он |
|
|
си |
JS |
2 |
К |
g |
§ |
О |
|
о |
& * |
ч |
|
^ |
с |
а |
|
S |
° |
<и о |
|
я |
н |
» |
4 |
а |
о |
° |
|
СЗ * |
|
||
о |
к |
к |
|
s |
о |
|
|
=SО |
С? |
|
|
d |
о |
о |
|
*3 |
S |
|
|
о |
|
\о |
|
н |
|
^ |
|
Номер |
слоя |
|
|
- с |
о |
о |
л? |
-sF |
? |
о |
о |
о |
о |
||
- |
|
|
-iF |
-sT |
? |
о |
о |
о |
о |
о |
о |
о |
|
-iT |
F? |
-sT |
? |
о |
о |
ъ |
о |
О |
|
- |
|
|
£ |
-SF |
? |
|
|
|
|||
- |
|
-х |
|
-sT |
? |
- с |
-iT |
н |
л? |
-SF |
? |
Н |
н |
н |
н |
н |
|
- |
н |
н |
-iF |
-sT |
? |
н |
н |
н |
н |
||
|
0 |
0 |
|
1 |
? |
ч |
О |
"э- |
"э- |
||
О |
о |
||||
-н |
|
|
-н |
|
-н |
'~е |
-iF |
-iF |
JF |
-SF |
? |
^СЗ |
^сз |
^СЗ |
'~СЗ |
'~СЗ |
|
|
-iF |
-iF |
JF |
-SF |
? |
|
|
|
' i |
' i |
' i |
~ : |
-£ : |
-£ : |
со :: |
-ie : |
й; |
АЛ. Неоднородные по толщине многослойные оболочки |
229 |
Аналогично, номера &2, £3,^4 в табл. А.1 соответствуют слою, арми рованному в перпендикулярном направлении (р = 90°), двойному слою («бислою») с углом армирования ± р и ортотропному слою. Заметим, что за толщину h(k^ «бислоя» следует принимать суммарную толщину двух монослоёв, составляющих «бислой».
В соответствии с исходными данными, помещенными в табл. А.1, производится расчёт многослойного, неоднородного по толщине компо зитного пакета. Однако в табл. А.1 заданы лишь свойства ортотропного слоя под номером к^. Свойства упругости перекрестно армированных (±</?) и однонаправленных (р = 0°, р = 90°) слоёв остаются неизвест ными.
Расчет неизвестных характеристик упругости для перекрестно ар мированных слоёв, а также для слоёв, армированных в направлениях 1 {р = 0°) и 2 {р = 90°), приводится ниже.
Характеристики упругости двойного слоя с перекрестным ар мированием под углами ± р . Перекрестное армирование с углами укладки ± р характеризуется тем, что на каждый слой, уложенный под углом р, накладывается такой же ответный слой, уложенный под углом (—р). Такая пара слоёв составляет ортотропный «бислой», оси упругости которого совпадают с направлениями 1 и 2 (рис. А.З).
Для каждого монослоя рассматриваемого «бислоя» задаются тол щина h, модули упругости Е\р, Е 2 Р в направлениях вдоль и поперёк армирования, коэффициенты v\p («сужение» при растяжении вдоль армирования) и щ р («укорочение» при растяжении поперек армирова ния), модуль сдвига Gp в плоскости монослоя, а также модули Gip, G2P поперечных сдвигов вдоль и поперек армирования соответственно.
Эти заданные или определённые экспериментально величины яв ляются исходными данными для получения характеристик упругости ортотропного «бислоя». Характеристики упругости «бислоя» определя ются следующими величинами:
Е\у Е2 — модули упругости в первом и втором направлениях соот ветственно;
щ, 1/2 — коэффициенты поперечной деформации (Пуассона);
Е1и2 — E2V1;
G — модуль сдвига в поверхности (12);
G\, G2 — модули поперечных сдвигов в первом и втором направле ниях соответственно.
В указанных величинах в данном случае опущены верхние индек сы & .
Необходимые характеристики упругости «бислоя» находятся с ис пользованием результатов [131:
тр |
а ( \ |
\ |
^ |
12 |
\ |
^ 1 2 |
; |
|
|
Ei |
= Ац(1 |
- щ1/2у, |
v\ = — |
|
V2 = -r- |
|
|||
|
|
|
Л22 |
|
Ац |
|
(А.11) |
||
Е2 |
= А22( 1 |
—щг'г); |
Е 11/2 = -Егщ = Ai2(l |
- щг'г); |
|||||
|
G\2 = G = А 33.
230 Прил. А. Неоднородные по толщине композитные конструкции
А33 = |
Gv cos22p |
+ | |
[E\v + E 2ip - 2E\ipv2!p] sin2 2p\ |
||
L~'l<p |
1 |
(p^2cp |
> |
±yZcp -i |
(А. 12) |
|
1 |
|
1 |
^\ср^2ср |
Модули поперечных сдвигов определяются соотношениями:
G\ = G\tpCos2ip + G2ipsin2p\
(A.13)
G2 = G\v sin2p + G2ipcos2p.
Для полной характеристики «бислоя» необходимо задать его об щую толщину hv=2 h = hM'3') (h — толщина монослоя) и расстояние av = а середины «бислоя» от опорной поверхности — в рассматри ваемом случае от внутренней поверхности оболочки.
Характеристики упругости монослоя с осевой и кольцевой армировками под углами р = 0° и <р = 9 0 °. Армировку с углом ар мирования <р= 0° составляют слои, уложенные вдоль оси 1 (рис. А.2). Из (АЛ 1), (А.12) и (А.13) при р = 0 получаем:
Е\ = E\v \ |
Е 2 = E2v; |
V\ = v\9\ |
v2 = v2lf\ |
G = Gtp] |
v |
G2 = G2lp\ |
(A.14) |
G\ = G\v \ |
(p = 0°). |
||
Таким образом, |
характеристики упругости монослоя совпадают |
с соответствующими характеристиками армирующих элементов. В этих формулах опущены верхние индексы к\.
Армировку с углом армирования р = 90° составляют слои, уложен ные в направлении 2 (рис. А.2). Из формул (А.11), (А.12) и (А.13) при р = 90° следуют значения характеристик упругости монослоя второго
направления: |
|
|
|
Е 1= E2ip’, |
Е2 = E \v \ |
v\ = v2^\ |
v2 — v\y\ |
G = Gip\ |
G \= G 2ip; |
G2 = G^-, |
(A.15) |
(p = 0°). |
Здесь опущены верхние индексы к2 .
Полученные по формулам (А.11), (А.13), (А.14) и (А.15) харак теристики слоёв, составляющих многослойный, неоднородный по тол щине композитный пакет, вносят в табл. А.1, которая содержит все необходимые исходные данные для расчёта жесткостей многослойного