Прикладная механика композитов
..pdfДинамика композитов с трещинами |
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ф |
1Ч д г 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) |
|
{х ху)т == 0»
Из симметрии относительно оси х вытекают следующие ус ловия:
(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) соответствует мгновенному приложению крутящего момента. Непрерывность сдвигового напряжения и углового перемещения задается ус-