Поведение конструкций из композитных материалов
..pdfАналогично шесть основных уравнений, включающих Мх, Му , Мху, Qx, Qy и w, получаются из (3.12), (3.14), (3.15), (3.65) - (3.67). Во всех случаях в связи с симметрией отсутствует связь между нагружением в плоскости пластины и в поперечном направлении. Хорошо известно и часто наблюдается на практике, что нагрузки, действующие в плоскости пластины, вызывают потерю устойчивости, т.е. поперечные прогибы, которые обычно приводят к разрушению. Для разъяснения этого пара докса отметим, что в данной главе рассматривается линейная теория упругости, а задача устойчивости, по-существу, является нелинейной. Для краткости изложения положения нелинейной теории упругости здесь не приводятся; с ними можно ознакомиться по работам [1 , 2 ].
После преобразования (3.83) с учетом нагрузок, действующих в плос кости пластины, получим
1>,— 7 |
Э4>v |
Э 4и> |
32w |
|
|
+ 2Z>3 |
+ А > Т Т |
-р(х, у) +Мх |
|
|
|
Эх4 |
Эх2Э^ 2 |
э У |
" з ? |
|
|
|
|
|
Э2W |
W |
(3.91) |
|
|
|
+ 2N,*■>’ дхду + К |
ду2 |
|
|
|
|
|
где видна связь между нагрузками в плоскости пластины и поперечным прогибом.
Критические нагрузки, подобно собственным частотам, независимы от поперечных нагрузок, которые в данном случае не учитываются. Однако на практике возможны случаи, когда совместное действие поперечных нагрузок и нагрузок, действующих в плоскости пластины, может привести к разрушению прежде, чем оно наступит в результате потери устойчивости. Но критическая нагрузка не зависит от поперечной нагрузки так же, как и частота собственных колебаний. Отметим, что при конструировании элементов, подверженных действию сжимающих нагрузок, которые могут вызвать потерю устойчивости, следует отдавать предпочтение пластинам, обладающим срединной плоскостью симметрии (Ву = 0). В противном случае возникает эффект продольного изгиба.
Обратившись теперь к уравнению (3.81), запишем дифференциаль ное уравнение устойчивости для композитной пластины под действием только осевой нагрузки (IVх). Отбрасьюая р(х, у ) , получим
_ |
Э4w |
Э4w |
|
(3.92) |
|
D}— т + 2D3~ |
|
||||
1 |
Эх |
4 |
J • |
а / |
дх2 |
|
|
дх2ду |
Опять принимая, что форма потери устойчивости для слоистой пласти ны из композитного материала соответствует форме потери устойчивости изотропной пластины с такими же граничными условиями, для шарнир но опертой по всем четырем краям слоистой пластины зададим
Е Атп sin |
. |
Л7ГУ |
sin—— |
||
а |
о |
|
т= 1 П=1 |
|
|
Подставляя (3.93) в (3.92), запишем выражение для критической |
||
нагрузки 7VX с/.: |
|
|
Л |
(3-94) |
Отметим, что вследствие однородности уравнения (3.92) |
величина |
Ат п , как уже отмечалось, не может быть определена. При определении нагрузки Nx сг существенным является только наименьшее ее значение. Однако мы не знаем, какая величина т или п соответствует наименьшему значению критической нагрузки. Все величины п стоят в числителе, поэто му наименьшему Nx сг соответствует п = 1. Величина т входит в различ ные члены, включающие изгибные жесткости (Dlt D2, D3) и отношение длины пластины к ее ширине а/Ь. Поэтому заранее не очевидно, какая величина т соответствует нижнему пределу Nx^ сг. Вместе с тем, для конкретной пластины, эта величина легко устанавливается.
Таким образом, в отличие от задач, рассмотренных в предыдущих разделах, в задачах на собственные значения определяются частоты соб ственных колебаний конструкции и критические нагрузки. Поскольку частоты собственных колебаний и величины критических нагрузок яв ляются собственными значениями, соответствующие формы колебаний и потери устойчивости являются собственными функциями. Следова тельно, при анализе любой конструкции необходимо определение четырех параметров: максимального прогиба, максимальных напряжений, частот собственных колебаний (если динамическое нагружение действует на изделие или передается ему) и критических нагрузок (если они являются сжимающими).
Что можно сказать о собственных колебаниях и критических нагруз ках для пластин из композитных материалов, граничные условия для которых иные, нежели шарнирное опирание? Для пластин, рассмотренных в данном разделе, двойные тригонометрические ряды были использованы для представления форм колебаний и потери устойчивости при условии, что края пластины шарнирно оперты.
Все комбинации форм колебаний для балок могут быть применены к пластинам с другими граничными условиями. Эти вопросы были рас смотрены в работе [7]. Все производные и интегралы для балочных функций сведены в таблицы для облегчения их использования и приведе ны в работах [8, 9]. Аналогичные выражения для форм колебаний и потери устойчивости можно найти во многих источниках, они также могут использоваться вместо (3.86) и (3.93).
Необходимо еще раз отметить, что собственные частоты и критические нагрузки определяются в данном разделе без учета поперечных сдвиго вых деформаций, и являются поэтому приближенными величинами. Однако ввиду относительной простоты этих расчетов, они могут быть рекомендованы как предварительные. Если учитывать поперечные сдви говые деформации, как это будет сделано в следующем разделе, то вели чины частот собственных колебаний и критической нагрузки будут иметь меньшие значения, чем определены здесь. Таким образом, определение критических нагрузок без учета поперечной сдвиговой деформации не идет в запас прочности.
3.11. СТАТИЧЕСКИЙ И ДИНАМИЧЕСКИЙ АНАЛИЗ ПЛАСТИН ИЗ КОМПОЗИТНЫХ МАТЕРИАЛОВ С УЧЕТОМ ВЛИЯНИЯ ПОПЕРЕЧНЫХ СДВИГОВЫХ ДЕФОРМАЦИЙ
Все результаты, приведенные в предыдущих разделах данной главы, получены без учета влияния поперечной сдвиговой деформации. Поэтому их следует рассматривать как приближенные и в виду их относительной простоты, рекомендовать к использованию в предварительных расчетах для первоначального определения размеров конструкции.
В данном разделе при рассмотрении пластин из композитного материа ла, подверженных действию статических и динамических нагрузок, учи тываются деформации поперечного сдвига.
Пренебрежение поперечной сдвиговой деформацией означает, как это было установлено раньше, что
(3.95)
Поскольку перемещения в плоскости пластины и(х, у , z) и v(x, у, z) образуются в результате смещения и поворота нормали к срединной плоскости, они имеют вид
и(х, у, z) = uQ(x, y) + za(x, у) |
^ щ |
v(x, у , z) = v0(x, y) + zfi(x, у)
с учетом (3.95) получили
(3.97)
И так как
«V- |
Эа |
ЭД |
к - |
Ч ай + |
Э^ |
9 х . К, |
9у |
||||
|
|
'' |
2 \ ду |
Эх / |
|
то следует, что |
|
|
|
|
|
|
d2w |
Э 2w |
|
|
32w |
Кд ~ |
Эх2 ’ |
ку = ду2 |
’ |
|
(3.99) |
|
ЭхЭу |
как это было показано раньше в (3.23).
Обратившись теперь к уравнениям равновесия и определяющим урав
нениям, получим для случая статического нагружения |
|
|||||
dN |
|
3Nxv |
|
|
|
|
Эх* + |
ду |
+ Т'* |
Т2^ = 0 |
(ЗЛ00)- |
||
3Nxr |
|
дК. |
|
|
|
|
дх |
+ |
by |
+ V |
- V - » |
<ЗЮ1 ) |
|
до |
|
3Q.. |
|
|
(3.102) |
|
~ t |
+ |
э; + р ( х . у ) - 0 |
||||
дМ |
|
3M xl. |
|
h г |
(3.103) |
|
Эх |
+ |
ду |
в * + 2 [т" + Т2' ] - ° |
|||
|
||||||
3М хг |
ЭМ, |
|
и |
(3.104) |
||
Эх |
+ Эу |
|
2 [т,, + т2 „ ] - 0 |
Для пластины, имеющей срединную симметрию (Ву =0) и при отсутст-
вии смешанных жесткостей [( |
) 16 = ( )г 4 = ( )45 = 0], определяю- |
щие уравнения имеют вид1 |
|
Nx = A ue°x + Ane°v |
(3.105) |
|
(3.106) |
Nxy = 2Л66с°ху |
(3.107) |
Мх = DU KX + Dn Ky |
(3.108) |
Му = DU KX4- Г>22ку |
(3.109) |
1 Здесь предполагается, что exz Ф J0 и Су2 Ф 0, т.е. углы поворота нормали к срединной поверхности пластины а и р уже не связаны с ее прогибом равенствами (3.97). Прим. ред.
Q x “ 2 J4 55€ X J — 2 Л 55 |
_ |
Эи> \ |
(3.111) |
|
“ + |
Эх) |
|||
|
|
|||
Qy ~ 2^44€y, = 2^44 |
|
|
(3.112) |
Поскольку пластина имеет срединную плоскость симметрии, здесь отсутствуют смешанные эффекты и усилия, действующие в плоскости пластины, не связаны с поперечными нагрузками, прогибами и углами поворота. Поэтому для анализа изгиба могут быть использованы урав нения (3.102), (3.104) и (3.108), (3.112), составляющие систему из восьми уравнений с восемью неизвестными.
Подставляя выражения (3.108), (3.112) в (3.102), (3.104) и исполь зуя равенства (3.98), получим следующую систему дифференциальных уравнений для слоистой композитной пластины, находящейся под дейст
вием |
поперечной нагрузки, для |
которой |
(Вц = 0) ( |
) 1 6 = ( |
) 2 6 = |
|||||
= ( |
) 45 = 0 и нет поверхностных сдвиговых напряжений. |
|
|
|||||||
|
Э2а |
Э2а |
, |
ч |
Э2В |
|
|
|
(3.113) |
|
D u — 2 +D66— |
+ ( D u + D6b) ^ - - 2 A — |
^ |
1 = 0 |
|||||||
|
дх |
|
ду |
66' |
3JC3>> |
|
|
|
|
|
|
|
|
|
+ °66J + |
|
Э2/? |
|
|
|
(З.П4) |
(*>12 + Аб) ъх ъу |
|
эу 2 |
2А“ |
|
|
|||||
|
|
|
|
|||||||
2А 55 |
0а |
0 2и> |
0 / 3 |
|
»t" vr |
+Р(х, +) = 0 |
|
(3.115) |
||
дх + |
дх2 |
+ 2/144 |
+ зу 2 |
|
||||||
|
|
ду |
|
|
|
|
||||
При решении |
этих уравнений |
в |
работе |
[1 0 ] |
был использован метод |
Навье. Для шарнирно опертой по всем четырем краям пластины, находя щейся под действием поперечной нагрузки, были использованы функции
|
> )= Е |
у с |
|
. |
ттгх |
. |
птгу |
(3.116) |
Н х , |
пт |
sin------ |
sin |
ь |
||||
|
|
а |
|
|||||
|
1Иа 1 п=\ |
|
|
|
|
|
|
«(*- у ) =
Р(х, у )=
ОО |
ОО |
ттгх |
. птгу |
(3.117) |
Е |
У А |
cos |
sin— |
|
1яы /*тп |
а |
о |
|
|
Ш=1 |
п=\ |
|
|
|
ОО |
ОО |
. ттгх |
птгу |
(3.118) |
Е Е в тп |
sin |
cos . |
||
а |
О |
|
/>(•*, j')^ E I |
f t l T T X |
sin |
П7Ту |
(3.119) |
Q» sin- |
|
m =1/1=1
которые удовлетворяют граничным условиям по краям: w = 0, bajbx = =0,н ах= 0,я и Э)3/Э^ = 0 на^ = 0, Ь.
Подставляя эти функции в исходные дифференциальные уравнения, получим
^ |
1 1 |
^ 1 2 |
^ 1 3 |
|
|
|
|
|
/ в " ’ ) |
|
|
^ 2 2 |
^ 2 3 |
|
^ |
1 3 |
L j 3 |
^ 3 3 |
. [ c „ , J |
|
О |
II |
О |
\ Q m n /
( 3 . 1 2 0 )
где <7m/J —коэффициент разложения поперечной нагрузки.
L \ \ i X ^ + + 2 4 5 5 ; L i 2 = (Z)12 + D 66)Xm\ n ;
Ь \ Ъ= 2Ass^m \ |
Ьгг = В еь \гт + D22 ^ + 2444; |
|
L 2 3 = |
2444X^1 |
Z/33 = 24 5 5X ^ + 2Л 44Х^ |
Решая систему (3.120), найдем |
||
J |
= (L \2L2I |
L 22L ]3)qmn |
m n |
|
det |
|
|
|
||
B„,„ = |
( ^ 12^13 “ |
L u L 23)qmn |
|
det |
|||
|
|||
Cmn = |
( ^ 11^22 “ |
^ 12)^/? |
|
det |
|||
|
(3.121)
(3.122)
(3.123)
где det - |
определитель матрицы в уравнении (3.120). |
|
Решив |
задачу |
нахождения а, 0 и w, можно_определить кривизны: |
кх = (Эа/Э*), |
= (ЪР1ду),кху = 1/2 [(Эй/Э7 + др/ду)]. |
Затем для слоистой пластины напряжения в каждом слое могут быть
найдены следующим образом: |
|
|
||||
|
|
e .i |
e.2 |
0 |
|
(3.124) |
a > |
= |
e.2 |
022 |
0 |
z |
|
CT* > . |
A |
0 |
0 |
2ёбб. A |
|
_ K x y |
Если поперечная нагрузка р(х, у) распределена по поверхности, то:
|
4 |
еа гь |
|
ттгх |
. птгу |
, , |
(3.125) |
|
|
Ть)0 )0р{х' у) s m T ~ |
Sin~ |
dxdy |
|||||
Когда эта нагрузка является постоянной, имеем |
|
|||||||
Я т п |
4р(х, |
У) (1 |
—cos ттг)(1 —cos птг) |
(3.126) |
||||
|
тптг |
|
|
|
|
|
|
|
Для сосредоточенной силы Р, приложенной в точке х =£ и у =rj: |
|
|||||||
|
4Р . гтт£ |
. |
ттт) |
|
|
|
||
<lm»~ |
ab |
Sm |
a |
Sm |
b |
|
|
(3.127) |
Для нагрузки, действующей на прямоугольный участок со сторонами
и и г |
и центром |
с координатами £ и т? (рис. 3.4), величина qmn опреде- |
||||||
|
|
|
|
1 |
m |
Ф |
1 |
|
ляется следующим образом: для п/Ъ Ф — |
и — |
— |
|
|||||
|
|
|
|
v |
а |
|
и |
|
|
л , |
ч |
. ( т т \ \ |
( n7rv\ |
I ттги\ |
|
||
|
4р(ху) 8Ш( — ) c o s ( ^ - ) |
c o s ( - ^ - ) |
(3.128) |
|||||
Я т п |
|
|
|
|
|
|
|
|
щяп/Ь= 1 /v |
и т/а= 1 /и |
|
|
|
|
|
||
Ч т п |
= О |
|
|
|
|
|
|
(3.129) |
|
|
|
|
|
|
|
||
|
Э2а |
Э2а |
, |
. Э20 |
„ , |
/_ |
Эи> \ |
,8 2а |
£>" Эл2 + ° ЬЬ9 v2 + ( ° ' 2 + 1>м |
дхду |
2А>5\ а+ дх) |
1 bt2 |
|||||
|
|
by |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
(3.129) |
Для любого иного поперечного нагружения qmn определяется из урав нения (3.125).
Если необходимо найти собственные частоты слоистой композитной пластины со срединной плоскостью симметрии (Вц = 0) и при ( ) \ 6 =
= ( )2 6 = ( |
) 4 s = 0 , но с учетом поперечных сдвиговых деформаций, |
|
то в уравнении |
(3.115) надо положитьр(х, у) =J) и ввести - ph(d*w/dt2) |
|
в правую часть. Кроме того, поскольку а. и j3 являются |
независимыми |
|
переменными, угловые колебательные движения элемента |
нормали к сре- |
динной поверхности пластины приведут к появлению в правых частях уравнений (3.113) и (3.114) членов /(Э2а/Э^2) и /(Э20/Эг2). Таким образом,
|
д2« |
п |
|
э2^ |
? j |
п |
9VV |
д2/3 |
(Z),2 + Dbb) Ъ хЪ у+ 0« Ъ х г |
|
и Ъу* |
4 |
J2+ |
ду |
= I |
||
|
dt2 |
|||||||
|
|
|
|
|
|
|
|
(3.130) |
. З а |
3 2w |
|
30 |
32w\ |
, d2w |
|
(3.131) |
|
5s| a ^ + |
э ! 7 |
+ 2/1d4 —— (- |
Р |
3/J |
|
|||
|
ду |
+ э / Г |
|
|
Рис. 3.4. Функции распределения нагрузки, действующей на элементарный участок пластины
Рис. 3.5. Различные функции распределения нагрузки:
а - синусоидальная, F (Г) |
= F 0 sin |
(яг/г х ) |
б - ступенчатая, F (Г) = F 0; в - триан- |
|
гулярная, F(f) |
= F 0 (1 - |
t/t0); г - |
экспоненциальная (нагрузка взрывом), F(t) = |
|
= F ne ~at\ д - |
ступенчатая триангулярная |
(нагрузка от ядерного взрыва), F (Г) = |
||
= /7о (1 |
F(f) = F, (1 -f/f,) |
|
где р определяется равенством (3.85) и
/ = рЛ3/12 |
(3.132) |
Аналогично процедуре, проделанной в разд. 3.5, для шарнирно опертой пластины примем:
|
|
|
|
|
|
тттх |
. |
Л7ГУ |
iw, |
(3.133) |
»(х, у, |
0 |
= |
£ |
£ |
С£„ sin------ |
Sin— -г— |
~ |
|||
s |
Ъ |
е |
|
|||||||
|
|
|
m= 1 я=1 |
|
а |
|
|
|
||
|
|
|
|
|
|
|
|
|
||
a ( x , y , t ) |
= |
£ |
E ^ . c o s |
W 7TX |
. |
П7Гу |
iw/ |
(3.134) |
||
------ |
sin—г 1- |
е |
||||||||
|
|
|
т =1п =1 |
|
а |
|
Ь |
|
|
|
|
|
|
|
|
|
|
|
|
||
—/ |
ч |
|
0 0 |
оо |
|
ттгх |
|
птгу |
■ |
(3.135) |
|
|
Е |
|
|
||||||
Р(х, у, о = £ |
Д7 « п — — COS~T ~ е |
т =1 /7 =1
Подставляя эти выражения в основные уравнения для динамического нагружения, получим систему однородных уравнений для определения собственных частот колебаний:
(3.136)
где величины L без штриха были определены ранее, и
, |
ънъ |
3 |
, |
р Л 3 2 |
1,1 “ 1,1 |
^7 |
“ m»’ ^ |
2 “ 122 - |
“ мл- |
Лз'з = Z-зз -
Как было показано в работе [2], три значения собственных частот можно получить, решая (3.136) для каждой пары величин т и п . Однако есть две частоты, которые являются значительно большими по сравнению с третьей, поскольку они связаны с инерцией поворота нормального эле мента и определяются правыми частями выражений (3.129) и (3.130). Эти частоты редко являются важными при анализе поведения конструк ции. Если ими пренебречь, то = L п и Ь}2 = L 22 и квадрат оставшейся собственной частоты может быть найден из уравнения
= \ Q L 33 + 2 ^ 1 2 ^ 2 3 ^ 1 3 “ ^ 2 2 ^ 1 3 ~ ^ П ^ 2 з ] / p h Q |
(3.137) |
где Q - L I \L2 2 - LI 2 , а также
A' |
L \2L 23 |
L22 L Jз |
-- |
c |
|
S*HI It |
||
|
_ ^12^13 |
(3.138) |
д , |
^11^23^/ |
|
^ m - l n - P |
7 |
|
Если влиянием поперечной сдвиговой деформации пренебречь, то выра |
|||
жение (3.137) |
будет идентично |
возведенному |
в квадрат выражению |
(3.87). |
|
|
|
В работе [2] |
также получены |
решения для |
слоистой композитной |
шарнирно опертой пластины, рассмотренной в данном разделе, при дина
мическом поперечном нагружении р(х, |
у, |
t). Инерция поворота при этом |
|||||||||
не учитывается и используется интеграл |
P(t) |
типа |
свертки, |
который |
|||||||
также известен |
как интеграл |
суперпозиции и интеграл Дюхамеля, |
т.е. |
||||||||
« ( * > |
у, О |
|
|
|
|
|
|
|
|
|
|
, |
i f |
У |
( i - . ) * * - „ L u - L |
11L n ) |
ш |
» |
si„ 2 g f . ( , ) |
|
|||
|
|
|
|
|
|
|
|
|
|
(ЗЛЩ |
|
|
|
|
|
L n L n — L U L22 |
. |
питх |
гиту |
p ( t ) |
|||
|
|
|
i ( |
t ) b |
|
—— — sin----------- C O S — г - |
|||||
|
|
|
|
е |
|
|
а |
Ь |
|
||
|
|
|
т = \ п = 1 V |
7 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
(3.140) |
|
где |
|
|
|
|
|
|
|
|
|
|
|
P(r) = jfV (r) |
sin ымя(/ - |
r)d r |
|
|
|
|
|
(3.141) |
где qmn — коэффициенты разложения поперечной нагрузки в ряд вида (3.119). Таким образом, если для заданной поперечной распределенной
нагрузки р(х, у, t) |
решение дается в виде (3.138), (3.140), то компонен |
|
ты кривизны кх, |
ку и кху для пластины могут быть найдены по выра |
|
жениям (3.99), а напряжения в каждом слое согласно (3.124). |
||
Функция |
P(t) |
была аналитически найдена для нескольких видов |
нагружения |
в соответствии с рис. 3.5. |