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

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

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

5.4. Для цилиндрической оболочки из задачи 5.3 со свойствами материала, данны­ ми в задаче 5.1, принять, что волокна имеют осевое расположение. Края оболочки

защемлены

с обоих концов, нагрузка р = 0,7 МПа и Nx = 0. Определить: а) ах и

ое при х = 0;

б) ох и (70 при х = 380 мм, т.е. при х = L/2.

5.5. Рассмотреть круговую цилиндрическую оболочку длиной L = 1,27 м, радиу­

сом 254 мм

и с толщиной стенки И =0,5 мм, полученную из материала задачи 5.1,

где Dx = Dxx, DQ =D22. Какова длина зоны краевого эффекта для этой оболочки? 5.6. Рассмотреть оболочку, образованную квазиизотропным композитным мате­ риалом, со свойствами Ех = EQ =126 ГПа и v = 0,3, радиус оболочки 254 мм, а общая толщина стенки 7,62 мм. Длина оболочки 1,02 м. Определить: а) длину зоны крае­ вого эффекта; б) мембранные напряжения оболочки ах, если она нагружена осе­

вым растягивающим усилием, равным 45,3 КН?

Глава 6. ЭНЕРГЕТИЧЕСКИЕ МЕТОДЫ АНАЛИЗА КОНСТРУКЦИЙ ИЗ КОМПОЗИТНЫХ МАТЕРИАЛОВ

6.1. ВВЕДЕНИЕ

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

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

ивариационный принцип Рейснера. Последний широко использовался в работе

[1].Теорема минимума дополнительной энергии используется сравнительно редко. В данной главе будет описана и эффективно использоваться теорема минимума потенциальной энергии, поскольку она позволяет решать все задачи, связанные с расчетом конструкций из композитных материалов.

Как было отмечено выше, альтернативой постановки задачи основанной на диф­ ференциальных уравнениях для упругого тела, такого как балка, пластина или обо­

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

6.2. ТЕОРЕМА МИНИМУМА ПОТЕНЦИАЛЬНОЙ ЭНЕРГИИ

Потенциальная энергия любого упругого тела в общем виде может быть записана уравнением

(6.1)

где W - удельная энергия деформации (см. разд. 6.4); R —объем упру­ гого тела; 7), Uj и - /-той компоненты поверхностной нагрузки, дефор­ мации и объемной силы соответственно; ST - часть поверхности тела, на которой заданы нагрузки. Первый член в правой части выражения (6.1) представляет собой энергию деформации упругого тела, второй и третий - работу поверхностных и объемных сил.

Теорема минимума потенциальной энергии формулируется следующим образом: ”Из всех перемещений, удовлетворяющих условиям совмест­ ности и заданным граничным условиям, те, которые удовлетворяют уравнениям равновесия, соответствуют минимуму потенциальной энер­

гии”.

 

Математически это выражается в виде (5 -

символ вариации)

8V=0.

(6.2)

В задачах механики используются следующие три операции:

= 8 ( (Lc)’ Ъ { у 2) = 2Уй У’f syd x = 8f y dx

(6 -3)

Содержащаяся в (6.1) удельная энергия деформации W в декартовой системе координат имеет вид

И/ = i Gifu =

+ W r + 2 ° / . -

+ О*./*, + 0 * /* . + ау*у:

При использовании теоремы минимума потенциальной энергии в соот­ ношении (6.4) напряжения заменяются через деформации, а деформации через перемещения.

6.3.РАСЧЁТ БАЛОК

Вкачестве примера использования теоремы минимума потенциальной

энергии рассмотрим изгибаемую балку (рис. 6.1). Балка имеет длину L в направлении х, ширину Ъ и высоту h. Она подвергнута действию по­ перечной распределенной нагрузки q(x) в направлении z. Модуль упру­ гости изотропного материала балки Е. Напряжения и деформации свя­ заны соотношением

ох = Еех

 

(6.5)

а деформации и перемещения равенством

 

с

d и,

d2w

(6.6)

dx

z

 

d x 2

 

z

£t

q(x)

Рис. 6.1. Балка под действием поперечной

 

 

распределенной нагрузки q {х)

 

t t / t .4- f

t

поскольку для изгибаемой балки их = -z(d w /b x ). Обратившись к выра­ жениям (6.4) —(6.6) и помня, что в элементарной теории изгиба балок

ау =аг = аху = exz = eyz =oyz =0, получим

W = \ ал

 

 

(6.7)

Поэтому энергия деформации U имеет вид

f L г + Ь / 2 r + h / 2 .

( d 2 W \ 2

2 1 j 1

и = (

Щ

- п \

^ dzdydx

j О J - Ь /2 J - h / 2

\ d x I

 

ш

dx

 

(6.8)

2

 

 

 

 

где / = bh3/12. Аналогично, для работы поверхностных сил в выражении (6.1) имеем

;т(и ^ =

L

 

 

 

 

f

q(x)w(x)dx.

 

S j*

о

 

 

 

 

 

Тогда уравнение (6.1) примет вид

 

v ~ f i

L(

i ?

) i x

~

i L‘,( x ) w u ) '**

(б-5>

Решая уравнение (6.2) с учетом (6 3 ), найдем

 

SV= 0 = El

•'о

\ dx

I

d x - f Lq (x)S w (x)d x

(6.10)

 

 

Jo

 

Символ вариации (6) может быть введен под интеграл, поскольку операции варьирования и интегрирования взаимоперестановочны. Вариа­ ции Е, I и q(x) отсутствуют, так как они являются заданными величинами.

Интегрируя первый член правой части выражения (6.10) по частям, получим

V2

rL&2W d2(6w)-dx

 

 

= E I fro

dx 2

d x 2

 

 

 

 

d2w / dw

L

 

*Ld3w d(8w)

 

 

 

 

 

0

 

l

dx

 

 

 

Е1\отз*гЛ

 

d2w d w

L

El

d3w с

L

 

 

 

E I ^ - S

(£) 0

, 5w

0

 

d x 2

 

d x 3

- n

i L—

. 6vt>dx

 

 

( 6 .11)

 

•'o d x 4

 

 

 

 

Подставив выражение (6.11) в (6.10) и проведя перестановку, будем иметь

8V = 0 =

rd2w Л/ dw '

L

__d3w _

 

E l^ - 8

 

El

fiw

J

 

d x 2

 

0

d x 3

 

г

F .d 4w

,

. Swdx = 0

 

(6.12)

ax

 

 

 

 

 

Для того, чтобы интеграл обратился в нуль должно выполняться усло­ вие

Г ,d W ( \

(6.13)

E I — ^ = q{x)

ах

 

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

Отметим также, что для выполнения равенства (6.12) необходимо, чтобы и первые два члена были равны нулю. Отсюда при х = 0 и х = L, т.е. на каждом конце балки, или EI(d2w/dx2) = - Мх = 0, или производ­ ная dw/dx должна быть заданной, т.е. ее вариация равна нулю. Кроме того, или EI(d3w/dx3) = - Vx =0, или должен быть задан прогиб w.

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

6.4. ПРЯМОУГОЛЬНАЯ ПЛАСТИНА ИЗ КОМПОЗИТНОГО МАТЕРИАЛА НАХОДЯЩАЯСЯ В УСЛОВИЯХ ПОПЕРЕЧНОГО НАГРУЖЕНИЯ

И ГИДРОТЕРМИЧЕСКИХ ВОЗДЕЙСТВИЙ

Формулировка теоремы о потенциальной энергии

В работе [3] детально показано, как осуществляется анализ ком­ позитных панелей с учетом анизотропии, поперечной деформации сдвига, термических и гидротермических воздействий. Уравнения связи между напряжениями и деформациями были записаны в форме (236) , а входя­ щие в них величины определены ранее. Геометрические соотношения записываются в форме (2.41), а из них следуют равенства (2.42). Если пренебречь е2 и о2 , то определяющие уравнения для слоистой пластины упрощаются и записываются в форме (2.62). В работе [3] показано, что

идля рассматриваемой задачи этими величинами можно пренебречь. Используя теорему минимума потенциальной энергии, уравнение,

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

+ ау [ ( у - ауЬ Т - 0 УАт] + ахг[2ех :]

+ V- [2«>•-•] +'а*у [2е*>- «яу& Т - &,Am] }dzdA

(6.14)

Здесь А представляет собой площадь пластины, размеры ^которой: 0< х< а, 0 < у < 6 и - h/2<z < й /2.

Теперь, подставив определяющие уравнения (2.62) и соотношения между деформациями и перемещениями (2.40) в (6.14), получим

■ - Л 1Ш

Эм0 дР + В12 дх ду

'

п

Эц0

да

Р\\

/ Эос\2

Эм0

ди0

 

11

Эх

Эх

2

\ Эх )

12 Эх

Эу

Эу0 Эа

+ в , М

¥ ч - л .

 

а«0 а«0

Эи0 ди0

Эд>

Эх

16

Эх Э^

Эх Эх

 

12

Эх

 

+ D,

да да

да Э/3

+ В 16

Эи0 Эа

Эм0 Э/?

 

Эм0 Эа

Эи0 Эа

 

Эх ду + Эх Эх

 

Эх

ду

Эх

Эх

ду

Эх

Эх

Эх

, -^22 / dv0 \ 2

 

dv0 dfi

D22tdf} \ 2

 

dvn Э

 

dv0 dv0

 

 

и иго .

+ ~ \ J V J + В ^ Ь ~ \ Ь 1

 

26 ду ду

ду Эх

+ Я

Э^о Эа

Ч

М

Эи*, Э/?

Зу0

Э/?

+ D

3

a f + f

f

7 6

Эх

Эд>

ду

Эх *

ду

ду + Эх

ду

 

’2 6

Эд>

Эд>

Эх Эд>

 

- п

, -3w

,

и 9w

9w 9w

а2

Эн»

1 / Эи» \ 2

4 5

^ + “ 87 + ^ 3 7 + 8 j 9 7

+ а 55 , Т + “ э 7 + г ( э г ) .

+ Л4 4

£ +'SH(S)1 +4

3 (£),+£& +M&)’

+ Д'66

Эм0 Эа

Эu0

dfi dv0 Эа

dv0 Э/J

Эд> Эд>

Эд>

Эх + Эх Эу + Эх Эх

 

 

 

 

Эм0

-17 ("5+ N ~J ) -(■v +V )1(■<+ ) -If(■«5+«С)

-If(«/;+л/д)-( +

 

1|)(

™)

+ Г* + М* + Л/Г* - р (х, y)w (x, у) |сЫ

(6.15)

где все величины, исключая перемещения, определяются равенствами (2.50)—(2.53), (255)—(257) и

Т*=

- 1

 

7

Щ к [«,■]*[«/]* A T 2 (z, t)dz ,

 

2*=

1 hk _

j

М*=

-

Z

}*

/ ) * ,

 

2 *=

!* * _ !

 

__

ЛГ

 

А*

-

МГ*=

S

 

/

[G,y]*K] [&j]AT{z, t)Am(z, i)dz .

 

k=

1

hk A

 

Равенство

(6.15) является более общим. Без члена, соответствующего

поверхностной нагрузке, оно содержит тридцать членов, определяющих энергию деформации для панели из композитного материала. Продолжая обсуждение, обратимся к равенству (2.62). Очевидно, что если в слоистой пластине деформации растяжения и сдвига не связаны (А 16 = Л2в = 0), то два члена исчезают. Если не связано растяжение с кручением (В 16 = = В 2 6 = 0), то исчезают еще два члена. Аналогично, не будет еще двух членов; если нет связи изгиба и кручения ф 16 =£>26 = 0). Когда слоис­ тая конструкция симметрична относительно срединной плоскости, что является очень распространенным случаем, то исчезнет пдть членов, вклю­ чающих смешанные жесткости Вц.

Распределение влажности и температуры

Рассмотрев ранее качественно гидротермические воздействия, необхо­ димо перейти к количественному определению AT(z) и Aw(z) при заданных условиях на поверхностях. Очевидно, что они в общем случае являются функциями х и у. Однако в тонкостенных конструкциях, таких как пластины и оболочки, наиболее существенным является расп­ ределение температуры и влажности по толщине.

В работе [4] анализ распределения концентраций влаги по толщине тонкой балки, пластины и оболочки осуществляется на основе класси­ ческого уравнения диффузии:

л И Д 2 ) _ 9 ( Д т ) _ 0

~ Э Z 2 9 /

где Dz - постоянная диффузии влаги по толщине композита; t время. В случае, если концентрация влаги мгновенно возрастает На величину Аш0 при t = 0 на верхней и нижней поверхностях конструкции, то реше­ ние (6.16) дает

 

ОО

Aw(z, t) = Аm 0 +

£ m n cos anz

 

n-0

где

 

(2л + 1) TT

(6.18)

=

" тг \ 2и + 1

(6.19)

/

Если концентрация влаги возрастает на Ат0 при t = 0 только на верх­ ней поверхности тонкостенной конструкции, то решение имеет вид

Д/и(г, r) = Am0| l - £ Р„cos ?„(2 + Л /2) |

(6.20)

где

 

 

(2л + 1)тг

(6.21)

<7" “

 

И

 

 

4 f

( - 1 ) " с -„>0,

(6.22)

7Г^ (2л + 1)

 

Решения (6.17) и (6.20) справедливы также для мгновенного тепло­ вого воздействия, если Ат заменить на АТ и учесть, что Dz —коэффи­ циент температуропроводности.

Для установившегося процесса переходный режим затухает и решения вырождаются в существенно более простые линейные функции распреде­ ления по направлению z :

(6.23)

Очевидно, что такое же выражение будет иметь место и для AT(z) . По этим решениям могут быть найдены усилия и моменты, соответст­ вующие напряжениям, вызванным воздействием температуры и влаж­

ности.

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

Граничные условия

Найдены решения для трех видов граничных условий: все края шарнир­ но опертые, все края защемленные и два противоположных —шарнирно опертые, а два других - защемленные. Если необходимо учитывать несимметрию укладки слоев, то вследствие появления смешанных эффек­ тов эти граничные условия необходимо уточнить. Следуя работе [5], можно записать по два варианта граничных условий для шарнирно опер­ тых и защемленных краев:

51: wn = М„ = UQ= Nnl = О

 

52: w„ = М„ = N„ = u0l = О

(6.24)

Cl : w„ = a„ = u0„ = Nn, = 0

 

C2: w =a„ = Nn = u0l = 0

 

Здесь w - поперечный прогиб; u0 - перемещение срединной плос­ кости; М - момент; N - усилие; а - угол поворота; п - направление нормали к контуру; t —направление касательной к контуру.

Представление решения

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

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

В работах [6—8] рассмотрена задача изгиба пластины с учетом попереч­ ной деформации сдвига, при этом углы поворота а и 0 задавались как функции поперечного прогиба на основе анализа соответствующего реше­ ния для балки. Однако здесь будет использована форма представления перемещений и углов поворота, приведенная в работе [5].

Для граничных условий S 1:

w(х, / ) = I

I

sin \ „ , Х sin А„.у

(6.25)

w «>

« ( * > . У ) = Е

X, Гтп sin \ тх COS \„ у

 

'"odd

^'Hld

 

 

00

 

ЪО

 

(6.27)

P ( x , y ) = X

 

X

Л т„cos \ mx sin X„j>

^*odd ^<*dd

 

OO

oo

 

Иo(*. J O = E E Ц,,„ sin Xmjc cos X„.y

(6.28)

^ o d * )

^ i x J d

 

00

 

OO

 

(6.29)

v0(x, y )= £

 

£

Km„cosXmxsin X„>-

^ i x j i l

**(4 Jd

 

Для граничных Условий 52 первые три разложения сохраняются, а два последние момента имеют вид

 

0 0

0 0

 

“o(*> y) =

£

E ^m„cos Xmx sin \„ y

(6.30)

 

w,xJd **«чЫ

 

0o(*. J') =

00

oo

 

E E K,„ sin Xmjc cos X„.y

(6.31)

 

modd **t4ld

 

Для защемленных краев граничные условия С1 используются балоч­ ные функции, представленные в работе [9] и усовершенствованные в [5,10, 11]. Они имеют следующий вид:

ОО00

w(*. JO" X Е К,пФи„(х)ФЫ'(у)

т —odd nl4ld

00 00

“(*. У ) = Х X ГтпФат(х)Фап(х, у) m«4ld n«4ld

OO00

Р( х , у ) = ' Е Е Лфрт(х)фрт(х)фр„(у) m«Hldn.Hld

мо(*. .у) = Е Е

sin Xmjc cos Х„.у

W«4jd W,4ld

 

(6.32)

(6.33)

(6.34)

(6.35)