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

Прикладная механика композитов

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

Динамика композитов с трещинами

18Г

a UQ = 0. Продольное направление композита выбрано вдоль оси 2 . Соответствующее напряженное состояние можно опре­ делить по формулам

^ = 2 G [ i ( ^ - 4 f ) + T ^ V V ] ,

 

 

 

 

 

 

 

 

 

(7)

’- = о [ £ ( 2 £

-

&

)

+ К

£

+ * ) ] -

тг0 = T0Z = 0.

 

 

 

 

 

 

 

В цилиндрической системе

(г, 0, г ) :

 

 

 

 

dr2

'

г

dr

d z 2

 

 

 

Аналогично получаем основные волновые уравнения:

/ д 2ф ,

1 д ф

,

д 2ф \ ___

<Э2ф

д г 2

дг

 

dz2 )

d t 2

/ d2ty , J_ Jhj)__ г}з

,

d2\|) N

(8)

2 Ч d r 2 1 г dr

r2 1 d z 2 )

d t2

Величины ci и c2 определены уравнениями (4).

2.1.3. Сдвиг при кручении. Допустим, что цилиндрическое композитное тело динамически скручивается относительно оси 2 так, что имеется перемещение только в направлении оси 0:

UQ = w (г , 2, t) , ur = Uz = 0.

(9>

Существуют две ненулевые компоненты тензора напряжений:

( d w

w \

d w

(Ю >

T'S = G (.—

- — )■

Te* = G-^-

Это приводит к волновому уравнению

э

( d 2w

,. 1

ddw

ww . dd2w \

__ d 2w

( И >

C2 ^

' r

dr

r2 ' d z 2 )

d t 2

где c2 определено вторым из уравнений (4).

2.2. ПЛОСКАЯ ТРЕЩИНА, ПАРАЛЛЕЛЬНАЯ ВОЛОКНАМ

Рассмотрим однонаправленный волокнистый композит, показанный на рис. 1. Дефект в форме плоской трещины ши­ риной 2а находится в матрице со свойствами Gm, vm и pm, тогда как волокно обладает свойствами Gf, Vf и pf, а компа-

182

Дж. Си

зит свойствами GCt vc и рс!). Задачи о нормальном и сдвиго­ вом ударах в плоскости излагаются так, что суперпозиция со­ ответствующих им решений описывает нагружение в произ­ вольном направлении.

2.2.1. Нормальный удар. Дефект, имеющийся в матрице, испытывает мгновенный нормальный удар с амплитудой оо-

Волокно

Рис. 1. Удар по однонаправленному композиту, имеющему трещину в матрице, (а) Трещина, параллельная волокнам; (Ь) локальная зона у вер­ шины трещины.

Математически это можно задать с помощью единичной сту­ пенчатой функции Хевисайда H{t). Используя систему коор­ динат, показанную на рис. 1(a), получим следующие гранич­ ные условия при t > 0;

(<Уу)т = — <*0Н (0,

0 ^ | х | < а, У = 0.

( 12)

 

{х ху)т ==

Из симметрии относительно оси х вытекают следующие ус­ ловия:

(Uy)m = (xxy)m = 0, U I > а, у = 0.

(13)

Предполагается, что слой матрицы, содержащий трещину толщиной 2h, идеально связан с окружающим материалом, свойства которого описываются свойствами композита GCt vc и рс. Условия непрерывности перемещений и напряжений при у — ±Л для любых значений х и t имеют вид

(^jc)m (Цх)с>

(и„)т =

К )с 1 для всех х и у = dbh. (14)

{®у)т == (Оу)с,

(Хху)т =

(Хху)с)

о Упругие константы композита ствам составляющих его компонентов симостей, выведенных в работе [9].

Gc, vc можно рассчитать по свой­ Gm, vm и Gf, Vf с помощью зави­

Динамика композитов с трещинами

183

Эта задача была решена в работе [9] с помощью метода интегрального преобразования, разработанного в [3]. Уста­ новлено, что поле динамических напряжений у вершины тре­ щины обладает сингулярностью типа единицы, деленной на квадратный корень из радиального расстояния г, отсчитывае­ мого от вершины трещины, а угловое распределение выра­ жается через угол 0. Выражения для напряжений имеют вид

(ff*>m =

^

CO4

0( l - s i n { e s in f e ) +

 

K )m =

- ^ = - cos|

e ( l +

sin-j-e s i n | e^)+

(15)

( * x y ) m =

 

COS -

0 Sin -

0 COS - e +

 

где коэффициент интенсивности напряжений k\{t)

является

функцией времени:

 

 

 

 

 

ki(t) =

f\[(c2)mt/a]a0^/a.

(16)

Члены более высокого порядка по степеням г в уравнениях (15) незначительны, так как стремятся к нулю, когда г при­ ближается к вершине трещины (рис. 1(b)). Компонента на­

пряжения ( o z)m равна Vm[(Ox)m+ (Ог/)т].

Численные значения величины f\[{c2 )mt/a\ из уравнения (16) графически изображены при различных отношениях Gc/Gm и a/h на рис. 2 и 3. Неоднородности, являющиеся ре­ зультатом различия коэффициентов Пуассона и плотностей,

малы.

Следовательно, оправдано

предположение, что vc «

« vm= 0,29 и рс ~ PmГрафики

зависимости_нормирован-

ного

коэффициента

интенсивности

/г, (0/(о0 V a ) от

безраз­

мерного временного

параметра (C2 )mt/a при а/1г=

1,0 пока­

заны на рис. 2. Видно, что ординаты всех кривых сначала возрастают до максимума, а затем несколько уменьшаются.

Наиболее высокий

пик соответствует однородному случаю

Gc = Gm-

параметра — отношения длины трещины

Влияние другого

к толщине слоя матрицы — показано для случая Gc/Gm — 10 на рис. 3. Видно, что коэффициент интенсивности динамиче­ ских напряжений k\ (t) уменьшается с увеличением длины трещины. Существует тенденция к стабилизации трещины с увеличением ее длины. Излом на кривой, соответствующей, параметру a/h = 2,0 в области малых значений t, характе­ ризует прибытие волн, отраженных от поверхности раздела. Эта особенность не наблюдается на кривой, соответствующей значению a/h = 0,5, так как поверхность раздела более уда-

184

Дж. Си

Рис. 2. Зависимость коэффициента интенсивности напряжений от времени в случае нормального удара для различных отношений модулей сдвига.

лена от трещины и отраженные волны являются более слабыми.

2.2.2. Сдвиговый удар. Соответствующую задачу о сдвиго­ вом ударе с амплитудой то также можно решить методом, предложенным в работе [9]. Граничные условия для / > О определим в виде

(* х у )т = — Т0Я (/),

а, У = О,

(17)

К ) т = 0

О I х | <

,

 

 

а условия кососимметрии в виде

 

 

М т =

( о у )т = о , 1*1 > а ,

у = 0 .

(18)

Предположения о непрерывности перемещений и напряжений

(14)

остаются в силе.

С

помощью

локальной системы полярных координат г

и 0 (рис. 1 (Ь))

можно выразить сингулярные члены локаль­

ных напряжений через коэффициент интенсивности напряже­ ний £ г(0 :

W m =

- -^ = - sin j

0 (2 + cos j 0 cos | e) +

/ \

^2 (0 ■ 1 Q

1 Л

3 Q I

( 19)

(®Л> =

S i n - 6 c o s j 0 c o s - e +

CWm =

cos^ 0

( l - sin

sin I

e) + ...,

Рис. 3. Зависимость коэффициента интенсивности напряжений от времени в случае нормального удара для различных отношений длины трещины к.

шагу между волокнами.

(o2)mt/a

Рис. 4. Зависимость

коэффициента интенсивности напряжений от времени

в случае сдвигового

удара для различных отношений длины трещины

 

к шагу между волокнами.

186 Дж Си

где k2{t) определяется формулой

(0 = / 2 [ ( c 2 ) m t/a] Т 0 У а .

( 2 0 )

.Для вычисления числовых значений величины f2[{c2)mt/a\ использованы те же отношения a/h и Gc/G m. Сравнивая ре­ зультаты для k2{t) (рис. 4) с результатами для k\ (t) (рис. 3), видим, что характер кривых одинаков, за исключением того, что k2(t)> ki(t). Следовательно, сдвиговый удар более опа­ сен с точки зрения растрескивания матрицы.

В общем случае нормальный и сдвиговый удары могут происходить одновременно. Начало разрушения определяет­

ся комбинацией величин k\ (t)

и k2(t), достигающих

некото­

рого критического значения, скажем

 

 

 

 

 

F[k{{i),

k2(t)] =

Fc-

 

 

 

(21)

Конкретный вид уравнения

(21) обсуждается

ниже

вместе

•с критерием плотности энергии деформирования.

 

2.3. ПЛОСКАЯ ТРЕЩИНА, ПЕРПЕНДИКУЛЯРНАЯ ВОЛОКНАМ

Другой возможный

случай

ориентации

дефекта — трещи­

на, перпендикулярная

волокнам (рис. 5). Остальные условия

Матрица

 

такие же, как в разд. 2.2

Волокно

 

 

для

трещины,

параллельной

IT

 

 

волокнам.

 

В

рассматривае­

 

 

мом

случае

волокна

парал­

 

 

 

лельны оси у и, следователь­

 

 

 

но,

условия

непрерывности

 

 

 

перемещений

и напряжений

 

 

 

(14)

относятся к плоскостям

 

 

 

x = ± h

для

всех

значе­

 

 

 

ний у.

 

 

 

 

 

 

 

2.3.1.

 

 

 

 

 

 

 

жение. Пока трещина нахо­

 

 

 

дится в материале матрицы,

 

 

 

функциональная

зависи­

Композит

 

 

мость для

поля динамиче­

^с^с’Рс

 

 

ских напряжений у вершины

Рис. 5. Трещина, перпендикулярная

трещины определяется урав­

волокнам, в случае удара.

 

нениями (15). Поворот тре­

 

 

 

щины на

90°

влияет

на ди­

намические характеристики передачи нагрузки, поэтому из­ меняется интенсивность поля динамических напряжений:

*1 ( 0 = g 1 [(c 2)m tja] а 0 V а.

(2 2 )

Динамика композитов с трещинами

187

Рис. 6. Зависимость коэффициента интенсивности напряжений от времен» в случае нормального удара для трещины, перпендикулярной волокнам-

Граничные условия для трещины по-прежнему определяются уравнениями (12).

=

Численные значения

g\[(c2)mt/a\

Для

отношений a//i =

0,2; 0,4;

...; 0,8 при

условии, что

рс «

pm, vc « vm = 0,29

и

Gc/Gm =

10, приведены на рис. 6.

Для соответствующих

значений a/h

и Gc/G m кривые на рис. 6 лежат выше кривых

на

рис. 3.

В

отличие

от статического нагружения область,

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

2.3.2. Сдвиг. В работе [9] показано, что для трещины, изображенной на рис. 5, при выполнении условий (17) и (18) уравнение (19) для k2{t) принимает вид

(0 = ё2 [(c2)m t/a] То Vfl •

(23)

Рисунок 7 снова свидетельствует о том, что локальное дина­ мическое напряжение осциллирует со временем и имеет мак­ симум, соответствующий моменту прихода волны, отражен-

J______________I______________I_____________ L _

0,0

Z .0

4,0

6,0

8,0

W l *

Рис. 7. Зависимость коэффициента интенсивности напряжений от времени в случае сдвигового удара для трещины, перпендикулярной волокнам.

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

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

2.4. КРУГОВОЙ ДЕФЕКТ МЕЖДУ РАЗОРВАННЫМИ ВОЛОКНАМИ

Волновые уравнения (8) и (11) применяются в упруго­ динамических задачах, обладающих осевой симметрией или симметрией вращения. Это имеет место в случае волокон с круглым поперечным сечением, подверженных однородному растяжению или кручению. На рис. 8(a) показано место ве­ роятного начала разрушения, где концы двух разорванных круглых волокон образуют зазор толщиной 2h. В этом зазоре имеется круговой дефект радиуса а, связанный с системой прямоугольных декартовых координат х, у, z (рис. 8(a)).

2.4Л. Мгновенное растяжение. Круговой дефект мгновенно растягивается в направлении z так, что поле перемещений,

Динамика композитов с трещинами

189

Рис. 8. Мгновенное растяжение и кручение, прилагаемые к трещине, рас­ положенной между разорванными волокнами, (а) Круговой дефект между разорванными волокнами; (Ь) компоненты напряжений в цилиндрической системе координат.

описываемое уравнением (6), не

зависит

от координаты

0,

a UQ = 0. Требуемые граничные условия

в плоскости 2

=

0

имеют вид

 

 

 

 

 

 

(оrz)m =

<*0Н it),

 

а,

2= 0,

(24)

(Vrz)m =

0 ^ г <

0,

 

 

 

 

 

а из условий симметрии вытекает, что

 

 

 

 

{Дг)т ==(^rz)m == 9,

Г^

2 = 0.

(25)

Предполагается непрерывность перемещений и напряже­ ний на поверхностях раздела, где матрица соприкасается с концами волокон, т. е.

(^r)m — (^r)f>

(^z)m

(^z)fl для всех г и г = ± А . (26)

(<*г)т=

rz)m ==

rz)fI

С помощью метода интегральных преобразований в ра­ боте [9] определены следующие локальные напряжения в

190

 

 

 

 

 

Дж. Си

 

 

 

элементе, показанном на рис. 8(b):

 

 

 

(°>)m =

k\ — cos — 0 (1 — sin — 0 sin — 0^ +

• • •»

V r’ m

V2r

2

V

 

2

2

)

 

M m =

 

2 Vm COS у

0 +

 

 

 

 

(aJ"’ =

‘7

r

cosi

e( 1+

sini

0sinf

e) +

(27>

/

4

k\

(t)

 

1

Q

. 1

Л

3 Л .

 

 

v

T „ L =

— 17= r-

cos — 0

sin — 0COS — 0 +

 

 

rz m

л/2л

 

2

2

 

2

 

 

(^0z)m

(^r0)m

 

 

 

 

 

 

 

 

Отметим, что для кругового дефекта зависимости динамиче­ ских напряжений от г и 0 имеют такой же вид, как и для

плоского (15). Снова

уграницы

трещины

выполняется

условие плоской деформации, так как

соотношение (ов)т =

= vm[{ог)т + [ог)т] справедливо для уравнений

(27). Коэф­

фициент интенсивности напряжений k\(t) зависит от времени:

 

 

k\{t) = ^np\[{c2)mtla)o(i'yJa.

(28)

На рис.

9

показана зависимость величины

p i [ { c 2) m t / a ] от

времени

для a/h = 0,5; 1,0 и 2,0 в предположении, что рт =

= pf, vm «

v/ = 0,29 и Gf/Gm = 10,0. Видно что вначале зна­

чения k\{t)

очень быстро возрастают, достигая максимума»

а затем уменьшаются до величины ki = 2o0'\/a/n, соответ­ ствующей статическому решению. Когда размер дефекта ста­ новится больше величины зазора, пиковые значения k\{t) достигаются при меньших /, что объясняется более ранним приходом волн, отраженных от концов волокон.

2.4.2.Мгновенное скручивание. Допустим, что композит,

изображенный на рис. 8(a), внезапно скручивается так, чго и волокна, и матрица испытывают только угловое перемеще­ ние UQ (см. уравнение (9)). Волновое уравнение (11) следует решать при следующих граничных условиях:

(т:вг)т= — {г1/а)т1л/а,

0 < г < а , 2 = 0;

^

M m = 0,

г > а , 2 = 0,

 

где г\ — радиальное расстояние от плоскости х, у. Изменение сдвигового напряжения в уравнении (29) соответствует мгновенному приложению крутящего момента. Непрерывность сдвигового напряжения и углового перемещения задается ус-