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

Поведение конструкций из композитных материалов

..pdf
Скачиваний:
34
Добавлен:
15.11.2022
Размер:
10.12 Mб
Скачать

Аналогично шесть основных уравнений, включающих Мх, Му , Мху, 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 и нет поверхностных сдвиговых напряжений.

 

 

 

Э

Э

,

ч

Э

 

 

 

(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)

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.