2012 Механика № 2
УДК 539.4: 678.067
А.Н. Полилов, Н.А. Татусь
Институт машиноведения им. А.А. Благонравова Российской академии наук,
Москва, Россия
ЭКСПЕРИМЕНТАЛЬНОЕ ОБОСНОВАНИЕ КРИТЕРИЕВ ПРОЧНОСТИ ВОЛОКНИСТЫХ КОМПОЗИТОВ, ПРОЯВЛЯЮЩИХ НАПРАВЛЕННЫЙ ХАРАКТЕР РАЗРУШЕНИЯ
Волокнистые композиты со слабой полимерной матрицей или непрочными поверхностями раздела «волокно-матрица» обнаруживают некоторые специфические направленные механизмы разрушения (расслоения, расщепления, скол вдоль волокон, полосу сдвига - кинк и др.), и для численного описания этих видов разрушения необходимы нетрадиционные подходы, новые специальные модели разрушения и критерии прочности, учитывающие взаимодействие нормальных и касательных напряжений на плоскости разрушения. Показано, что использование простых линейных критериев прочности позволяет лучше интерпретировать экспериментальные данные и оценивать в опытах прочностные параметры целостности композитных структур, которые можно использовать при проектных расчетах композитных элементов конструкций.
Ключевые слова: разрушение, критерии прочности, укладка волокон, расслоение, циклическое нагружение, волокнистые композиты, полоса сдвига.
A.N. Polilov, N.A. Tatus
Institute of Machines Science of Russian Academy of Sciences,
Moscow, Russian Federation
EXPERIMENTAL SUBSTANTIATION OF STRENGTH CRITERIA FOR FRP SHOWING DIRECTIONAL TYPE OF FRACTURE
FRP (fiber-reinforced-plastics) with weak polymeric matrix and/or fiber-matrix interface demonstrate some specific directional modes of failure behavior (delamination, splitting, failure along the fibers, kink, etc.) and to describe numerically these failure modes it is necessary to use some non-traditional approaches, new specific models of fracture processes and strength criteria taking into account the normal-shear stresses interaction on the fracture planes. It's shown that using simple linear fracture criteria it's possible to understand the experimental data more correctly. The application of new proposed FRP fracture criteria makes it possible to estimate experimentally some structure integrity parameters being used in more correct design of composite elements of constructions.
Keywords: fracture, strength criteria, laying fiber, delamination, fatigue, fiber reinforced plastics, kink.
Введение
Многие авторы [1, 2, 3, 5, 8, 21 и др.] формулировали критерии прочности анизотропных композитов в виде некоторого полинома от инвариантов, представляющих собой свертки тензора напряжений ъ у с
так называемыми тензорами прочности четных рангов: 11 = ВГ!,оГ!,, 12 = = Луи ° у °1с1 и т.д. Число констант Вге для плоского случая равно трем, АуМ - шести, и нахождение такого большого числа констант из ограниченного числа опытов, как правило, оказывается неустойчивым: экспериментальные данные описываются одинаково хорошо различным набором констант, и некоторые из них приходится назначать произвольно (например, равными 0). Если ограничиться вторым инвариантом, критерий прочности приобретает вид
11+¡в=1, (1)
где в = 2 - Гольденблат-Копнов [3]; в = 1 - Цай-Ву [21]), и предельная
кривая в плоскости -а2 в случае растяжения в двух направлениях напряжениями и а2 будет представлять собой эллипс.
Приближение предельной кривой эллипсом для ортогонально армированного пластика со слабой матрицей вряд ли может считаться удовлетворительным. Действительно, пренебрегая прочностью полимерного связующего, мы получаем ортогональную решетку, которая может разрушиться в результате разрыва волокон в том или ином направлении. Решетка может выдерживать напряжения только вдоль волокон, но не под углом к ним, когда возникают касательные напряжения. В соответствии с этим условие разрушения в плоскости - а2
будет представляться прямоугольником. Конечно, вводя в уравнение (1) инварианты более высоких порядков, можно этот прямоугольник аппроксимировать с любой степенью точности, но число констант и соответствующие трудности их экспериментального определения при этом неизмеримо возрастут.
С другой стороны, тензорный полином (1) определяет гладкую поверхность в пространстве напряжений, и мы ничего не можем ска-
зать о смене механизмов разрушения. В то же время разрушение, например, ортогонально армированных пластиков при растяжении под углом к волокнам (см. п. 1) происходит по плоскости, содержащей одно или другое семейство армирования. Учитывая направленный характер разрушения волокнистого композита, естественно принять в качестве условия прочности некоторое предельное соотношение между напряжениями, действующими в плоскости возможного разрушения, включающей в себя одно из направлений армирования. Обозначая через ап и тп нормальное и касательное напряжение на этой плоскости, будем формулировать условие прочности в виде [7, 9]
фК , тп) =1 (2)
1. Растяжение ортогонально армированных пластиков
Для ортогонально армированных пластиков естественно применить условие (2), поставленное на одной из двух перпендикулярных друг другу площадок, содержащих волокна одного из семейств армирования. Простейшее предположение о виде функции ф в (2) состоит в том, что она линейна, и тогда можно выписать два критерия, поставленных на одной ( п1 ) или на другой ( п2) плоскости возможного разрушения,
CTn1 + m1Tn1 =СТ(°°) = с1, (3)
°n2 + m2 Тп 2 = а( 90°) = С2. (4)
Формулы преобразования тензора напряжений при повороте координатных осей на угол а
аn1 = о1 sin2 а + а2 cos2 а + 2т sin а cos а, ап2 = а1 cos2 а + а2 sin2 а - 2т sin а cos а, (5)
тп1 = тп2 = x(cos2 а - sin2 а) + (а2 -a1)sin а cos а
используем для случая растяжения под углом а к одному из семейств волокон о1 = т = 0, о 2 = о (а). Тогда из критериев (3), (4) получаем две
зависимости прочности о (а) от угла вырезки образца а :
c *
а(а) =---2------- --------- при а <а (6)
cos а + m1sin а cos а
а(а) = — ----------------------------------------------—- при а >а . (7)
sin а + m2sin а cos а
При нагружении под углом а* происходит переход от одной плоскости разрушения к другой, и а* находится из равенства (6) и (7), приводящего к квадратному уравнению относительно 1§а*
а* = arctg2 (m1s - m2 +*J(m2 - m,s)2 + 4s, (8)
где 8 = с2 / с1.
1.1. «Равнопрочный» пластик (0/90)
Как видно из рис. 1 (сплошные линии), зависимость (6) хорошо описывает изменение прочности ортогонально армированных пластиков при растяжении под углами от 0° до 45°.
Рис. 1. Зависимости прочности ортогонально армированных пластиков от направления
вырезки образцов: х - углепластик;-----линейный критерий; Д - боропластик;
------квадратичный; ° - стеклопластик; - • - • - инвариантный
При одинаковом числе слоев N и Ы2 в двух ортогональных направлениях можно считать с1 = с2, щ = т2 и а*=45° . Формулы (6) и (7) описывают симметричную относительно 45° зависимость прочности от угла: о (а) = а (90°- а) (рис. 2, в). При Ы2 = 0 имеем случай однонаправленного материала. А при неравных N и Ы2 можно считать прочности вдоль волокон пропорциональными числу слоев с1 / с2 = N / Ы2. Если допустить также, что сдвиговая прочность вдоль
волокон слабо зависит от числа слоев (это - по смыслу - прочность полимерной матрицы), то m2 / тг = c1/ с2 = N1/ N2, и окажется, что для известной структуры ортогонального армирования (известное число слоев N и N2) всю немонотонную зависимость прочности от направления растяжения (рис. 2, г) можно описать на основании всего лишь двух экспериментов, например для углов вырезки 0° и 45°. Угол смены плоскости разрушения нужно определять из условия (8).
На рис. 2, а данные с рис. 1 для углепластика нанесены в координатах ап1 = а(а) cos2 а, тп1 =т(а) cos а sin а . Прямая линия (3) в этих
координатах, приводящая к зависимости (6), одинаково хорошо описывает эксперимент, как и участок эллипса:
(9)
Критерий (9) при подстановке в (5) приводит к выражению
а (а) = (cos4а/а*2 + sin2а-cos2а/т*2) . (10)
Расчеты по формуле (10) приведены на рис. 1 штриховыми линиями, а штрих-пунктирными для стеклопластика показана обработка по формуле типа (1) из [2], хуже согласующаяся с экспериментом и приводящая к более громоздким расчетным выражениям. По своей простоте и удобству обработки результатов в координатах а„-т„ предпочтение следует отдать линейным зависимостям (3), (4). Интересно отметить, что кривые (6) и (7) немонотонны и для равнопрочного (0/90) композита образуют ю-образную нижнюю огибающую, т.е. угол
u /" \ t 1 * минимума кривой (6) а1 = arctg mj лежит левее а , а угол минимума
кривой (7) а2 = 2arctg(-m2) - правее а*.
1.2. Однонаправленный материал (0)
Наиболее естественно применение условия типа (2) к однонаправленным композитам [9, 15, 17, 22]. Предположение о линейном виде функции ф приводит к критериям (3) и (4), где, по сути,
C = о(0°) - прочность вдоль волокон; с2 = о(90°) - прочность при рас-
тяжении поперек волокон, хотя более правильно считать с1, с2, как и т1, т2, феноменологическими параметрами, определяемыми при обработке экспериментальных данных по прочности образцов, вырезанных в различных направлениях. Наиболее просто определять эти параметры, проводя по методу наименьших квадратов прямую линию через экспериментальные точки, нанесенные в координатах аи1 - ти1 или
°п2 -тп2. При малых углах а < а* происходит разрыв волокон, и зависимость прочности имеет вид (6). Часто в литературе по механике композитов полагают т1 = 0, что приводит к противоречащему опыту
результату: прочность растет с ростом угла разориентации. Условие (4)
до1
лишено этого логического порока: — < 0 при 0 < а<а2, где а'2 - кода
ордината минимума зависимости (7).
Критический угол (8), при котором происходит переход от разрыва волокон к разрушению матрицы, определяется равенством (6) и (7), из которого а* « агС^ (в/т2) при с2 / с1 =в<<1. Для однонаправленных пластиков угол а * мал (например, для углепластиков а(0°) «1000 МПа, а(90°) « 25 МПа, т2 « 0,5, а* = 3° ), и точное экспериментальное определение т1 в крайне узком диапазоне углов 0-3° вряд ли возможно. Поэтому удобно упростить зависимость (6), полагая
о(а) = о(0°) при а<а*. (11)
Максимальная погрешность критерия (11) при вычислении а(а)
по сравнению с (6) при а = а* составляет всего лишь 2 %, что вполне допустимо с учетом разброса экспериментальных данных. Естественно считать а* «оптимальным» углом малой разориентации волокон (аналог свилеватой и потому нерасщепляющейся древесины); при этом прочность вдоль волокон практически не снижается, а сопротивление росту продольных трещин расщепления возрастает многократно, исключая расщепления материала при циклических нагрузках и вблизи концентраторов напряжений. Зависимость (7) для однонаправленного пластика немонотонна, и наименьшая прочность не соответствует растяжению под 90° к волокнам. Угол а2 = -^агС^ (-т2), соответствую-
щий минимальной прочности, лежит в диапазоне 45-90°, поскольку параметр т2 положителен.
120
100
80
60
40
20
а(а)
;2 1
1 1
' і
ГЧ. А
я/4
л/2
а
Рис. 2. Предельная кривая для углепластика (0/90) в координатах ая-хя (а) и зависимости прочности ортогонально армированных композитов от направления вырезки образцов по критериям прочности (6) - кривая 1 и (7) - кривая 2; б - однонаправленный КМ (0), в - «равнопрочный» (0/90), г - (04/90і)
Ранее [9, 17] предлагался также квадратичный вариант критерия разрушения поверхности раздела:
V
= 1.
ча2 у
(12)
Различие в критериях (12) и (7) при обработке экспериментов незначительно. Удобство линейного критерия (4) состоит в простоте расчетных формул и надежности метода определения констант т2, с2 путем проведения прямой линии через опытные точки в координатах сти2 -ти2 (см. рис. 2, а), рассчитанные через критическое напряжение
а(а) по формулам (5). Отметим, что если через те же точки провести эллипс (12), то равенства одинаковых по смыслу величин о(90°) = о2*
б
а
в
г
и особенно с2 / т2 = т* не будут выполняться. Таким образом, значение
прочности внутрислойного сдвига, определенное при растяжении, зависит от используемого критерия, что необходимо учитывать в расчетах.
1.3. Критерии циклической прочности монослоя (0)
Простейший, чисто феноменологический, подход к построению критерия разрушения в условиях длительного или циклического нагружения, а также при воздействии различных факторов: температуры, влажности, агрессивной среды, облучения состоит в замене параметров критерия прочности (4) - с2 и т2 - на функционалы от истории нагружения, в частности на функции от числа циклов N, поскольку вид разрушения при циклическом растяжении однонаправленных композитов такой же, как при квазистатическом. В [15] по результатам испытаний на циклическое растяжение прямоугольных образцов из однонаправленного углепластика, вырезанных под различными углами к направлению армирования, строили для разного числа циклов: 104, 105, 106 прямые линии в координатах gn2 - тn1 (типа изображенной на рис. 2, а),
по которым определяли значения параметров для каждого уровня долговечности. На рис. 3, а приведены определенные таким образом значения m2(N) и с2(N), которые достаточно точно аппроксимируются следующими функциями:
с 2 (N) = с 2 (1 - 0,1 lg N); m2 (N) = m2 (1 - 0,19^lg N). (13)
Подставляя (13) в (7), получаем зависимость прочности однонаправленных углепластиков от направления растяжения и от числа циклов. В многоцикловой области (N > 103) значение т2 = 0,3 становится практически постоянным, что позволяет получить весьма простое выражение
а(а,N)= . С(1-0-1lgN---------, (14)
sin а + 0,3sin а cos а
описывающее, как видно из рис. 3, б, целую серию экспериментальных зависимостей критического числа циклов от приложенных циклических напряжений для разных углов разориентации.
1 10 104 10й N
а
о(а, АО, МПа
150
100
50
X X 1 X
X X X
— д сг. о /2 — < оо >оо
3 . . <
V 7-5 ■ V ^ У ,
10
104
10
10 N
б
Рис. 3. Зависимости параметров (1 - с2,2- т2) линейного критерия (4) (а) и прочности однонаправленного углепластика (14) (б) от направления растяжения а и от логарифма числа циклов N 1 -а = 10°,2 -а = 22,5°,3 -а = 45°, 4 -а = 90°
Основной эффект - немонотонные и негладкие зависимости прочности от направления нагружения, описываемые с помощью критериев, поставленных на площадках возможного разрушения, связанных с направлениями укладки волокон. Традиционные тензорнополиномиальные критерии прочности описывают гладкую предельную поверхность в пространстве напряжений (многомерный эллипсоид), и при этом нет возможности говорить об условиях того или иного механизма разрушения, что является известным недостатком при последовательной схеме послойного анализа, когда в зависимости от вида разрушения монослоя надо либо полностью исключать его из рассмотрения (после разрыва волокон), либо полагать нулем характеристики матрицы (после выполнения условия скола вдоль волокон).
2. Двухосное растяжение и кручение косоугольно намотанных труб ( ± а )
Вид разрушения намоточных композитных труб отличается от плоских образцов из-за отсутствия «эффекта перерезанных нитей» и из-за переплетения жгутов при намотке (нет четкой слоистой структуры). Поэтому в трубах происходит не скол вдоль волокон, а искажение решетки арматуры (п. 2.2). Меняется и вид критерия прочности, хотя подход формально остается таким же.
2.1. Феноменологические критерии прочности труб
На рис. 5 приведены данные из [13] для прочности при растяжении (рис. 4, а) и при кручении (рис. 4, б) стеклопластиковых труб со схемой намотки волокон (±а ). Формальные выражения для проекции напряжений на плоскость, содержащую направление волокон, для осевого растяжения напряжением о(а) можно представить как
Рис. 4. Зависимости прочности при растяжении (а) и при кручении (б) стеклопластиковых труб от угла намотки волокон ±а
О 50 100 150 200 250 300
Рис. 5. Данные по прочности стеклопластиковых труб в координатах си - тп при растяжении (•) и кручении (°)
an = a(a)sin2 a; тп = a(a) sin a cos a, (15)
для кручения со средним по толщине стенки касательным напряжением т(а):
Формальность подхода состоит в том, что при смене угла армирования а мы, по сути, имеем дело с другим материалом, а не просто с изменением направления нагружения, как рассматривалось выше.
Точки, соответствующие разрушению при растяжении, на рис. 5 можно аппроксимировать какой-либо кривой второго порядка, для простоты - полуокружностью:
откуда с использованием (15) получается чрезвычайно простая зависимость
трубы с намоткой (+45°).
Данные для кручения на рис. 5 также можно аппроксимировать участком эллипса, а для простоты - полуокружностью:
В (20) 2ск соответствует прочности на кручение трубы с намоткой (+45°).
Зависимости (18), (20) (кривые 1 на рис. 4, а, б) вследствие своей очевидной простоты очень удобны для инженерных расчетов и для выбора рациональных углов армирования. Учет ненулевой прочности связующего, существенный для углов намотки, близких к 90°, можно осуществить линейным преобразованием этих зависимостей (кривые 2 на рис. 4, а, б)):
ап = x(a)sin2a; тп =т(а)cos2a.
(16)
(17)
а(а) = 2срС£а = а(45°) С^а. (18)
По смыслу константа 2ср равна прочности а (45°) на растяжение
(18)
(19)
и с подстановкой (16) в (19) получаем
т (а) = 2ск sin 2а = т( 45°) sin 2а.
(20)
(а) = [а (45°) - а (90°)] х ctga + а (90°),
(21)
т(а) = [т(45°)-т( 90°)] х sin 2а + т( 90°). (22)
Зависимость (21) справедлива при а>а* =arctg(а(0°)/а(45°)); при а < а* считаем, что о (а) = а (0°).
2.2. Модель ромба из нерастяжимых волокон
В [11] для косоугольно армированной трубы была обоснована модель ячейки материала в виде ромба из нерастяжимых стержней (рис. 6), повороту которых препятствует распределенный момент р, действующий со стороны матрицы. Сила осевого растяжения пропорциональна длине малой диагонали Pz = 2аzl sin а, смещение вершины
Az = (2lcosа) с1а = -2Ыпаёа. Величина суммарного момента, препятствующего повороту стержней, считается пропорциональной площади ромба M = 2^,12sinacosa. Приравнивая работу внешней силы работе момента при повороте стержней на угол са, получим 2а J2 sin2 аёа = 2ц12 sin а cos аёа.
Рис. 6. Схема деформирования элемента косоугольно армированного композита при двухосном растяжении
Считая, что разрушение соответствует критическому значению ц* распределенного момента, характеризующему способность сопротивляться искажению решетки, находим выражение <зг (а) = ,
в точности совпадающее с формулой (18), которая независимо была получена из обработки экспериментальных данных.
При двухосном растяжении продольные аг (осевое нагружение)
и окружные а0 (внутреннее давление) напряжения, как видно из рис. 6,
стремятся исказить решетку в разных направлениях, что обусловливает (рис. 7) упрочняющий эффект при двухосном растяжении [6] и приводит к зависимости
(45°)
,(45°)
= 1.
(23)
Прямые скобки - обозначение абсолютной величины - говорят о смене знака выражения в зависимости от соотношения между напряжениями, и это приводит к двум параллельным предельным прямым линиям 1, 2 (см. рис. 7) на плоскости ае-а2, проходящим через точки на осях аг (а), а0 (а), соответствующим действию только осевых или окружных напряжений.
а б
Рис. 7. Экспериментальные данные по прочности стеклопластиковых труб с намоткой (±45°) (а) [6] и (±55°); (б) [16, 20] при комбинации осевого растяжения и внутреннего
давления
Для симметричной намотки (±45°) растягивающие и окружные напряжения, казалось бы, равнозначны, но окружные напряжения, например, в баллонах создаются внутренним давлением, и потеря несущей способности определяется не разрывом волокон, а потерей герметичности. Поэтому в реальности экспериментальные значения аг (45°)
и а0 (45°) оказываются различными. Если же ради методической простоты пренебречь этим различием, то уравнение (23) чрезвычайно упрощается:
|а z tga - Gectga| = а г (45°) (24)
и описывает две прямые линии 1, 2 в плоскости oe-oz (см. рис. 7). Прямая 3 соответствует разрыву волокон, когда перекашивание ромба невозможно. Для заданного отношения приложенных напряжений легко указать оптимальный угол намотки, обращающий в ноль разность в прямых скобках. Это означает невозможность перекашивания ромба и перенесение всех усилий только на волокна. Наибольшая прочность соответствует оптимальному углу намотки для заданного отношения приложенных напряжений tgP= oe/oz .
a* = arrtg^gP^ (45°)ae (45°). (25)
Например, в цилиндрической части баллона для сжатого газа («труба с заглушками») окружные напряжения вдвое превышают осевые oe = 2о z, и оптимальный угол намотки имеет известное из «нитяной аналогии» значение a* = arctgV2 = 54°44'.
Прочность при углах, близких к оптимальному, определяется разрывом волокон, а не перекашиванием ромба, и участок 3 на рис. 7 с достаточной точностью может быть заменен прямой, проходящей через точки на осях gz (0°), ое (90°), соответствующие нагружению вдоль волокон. Уравнение прямой 3 можно принять в виде
а0 cos2 a = az, a0sin2 a=ae, gz +Ge =g0 =gz (0°). (26)
Вообще говоря, для приближенного описания зависимостей прочности композитных труб от угла намотки волокон в области gz • Ge > 0 достаточно всего двух экспериментов для определения
прочности gz (45°) при растяжении трубы с намоткой (±45°) и прочности gz (0°) на растяжение вдоль волокон. Прочность при окружных напряжениях можно в первом приближении считать равной прочности на растяжение при соответственных углах намотки gz (a) = Ge (90° - a).
Далее через точки g z (a) и Ge (a), оцененные по формулам типа
(18), для каждого угла намотки a можно провести прямые линии под углом р = arctg (tg2a) (см. рис. 7), а используя gz (0°), построить прямую 3,
ограничивающую прочность по условию разрыва волокон. В области осевого сжатия добавляется один экспериментальный параметр а С2 (а), и уравнение предельной прямой 4 (рис. 7, б) имеет вид
~1 (27)
ог(а) ое(а)
Как видно из рис. 7, простая модель ячейки волокон, приводящая к кусочно-линейной предельной поверхности, хорошо описывает эксперименты на трубах с различными углами намотки.
При двухосном растяжении труб напряжениями ог и ое экспериментальные данные также хорошо аппроксимируются в координатах ог - ое участком повернутого эллипса, проходящим через точки на
осях ог (а) и ое (а), соответствующие прочностям при одноосном нагружении, и через точку на линии оптимального нагружения = 1§2а,
связанную из (26) с прочностью вдоль волокон о (0°):
о2 t + о0 оzае
2
о
2
2t2 -
о"(45°)(і +12)
о2(0°)
= t 2о2 (45 °), (28)
где о (0°), о (45°) - прочности на растяжение труб с намотками (0°) и (±45°); і = 1§а . Действительно, уравнение повернутого эллипса в осях ог - ое имеет вид
2 2 0z . 0Є 0z0Є = і
A2 Б2 C2 '
При одноосном растяжении (ае = 0) ^ а г (а) = а( 45°)/1 = A.
При чистом внутреннем давлении (az = 0)^ае(а) =
= а( 45°)-t = B.
При нагружении в «оптимальном» направлении ае = а( 0°) sin2 а =
а(0°) , ч 2 а(0)t2 1 2 (l +12)2
= V '• а, = а(0 )cos2 а= V ’ • — =-------------*-----
(1 +12) 2 1 + *2 ’С2 а2 (45°) а2 (0°) I2 '
Основной эффект: упрочнение косоугольно армированных труб или сосудов давления при двухосном растяжении, наличие оптималь-
ных углов намотки для данного напряженного состояния можно описать на основе модели ромбической элементарной ячейки. Для симметричной намотки ( ±а ) при двухосном напряженном состоянии критерий (28) позволяет проводить расчет напряжений и анализ прочности для каждой пары слоев. Это дает большие преимущества по сравнению с традиционным послойным методом, когда напряжения пересчитываются в каждом монослое, уложенном под некоторым углом к выбранной оси детали. Наклонный монослой не обладает упругой симметрией по отношению к осям детали, а симметричная пара слоев - орто-тропна, поэтому расчет по парам слоев значительно надежнее, так как более достоверны данные по прочности намоточных труб, чем данные, полученные на плоских образцах, вырезанных под различными углами к направлению укладки волокон. Изложенные критерии хорошо согласуются с экспериментами и дают возможность простых инженерных оценок критических напряжений и оптимальных схем армирования.
3. Кручение и сжатие с образованием кинка -наклонной полосы сдвига
3.1. Критерий наклонного разрушения при кручении
Некоторые композиты с разориентированным армированием разрушаются при кручении по плоскости, составляющей некоторый угол а с продольной осью 2. Ось х направлена в плоскости пластины, ось у -в направлении толщины.
За критерий разрушения можно по аналогии с (9) принять достижение зависящего от направления а критического значения квадратичной комбинации нормальных и касательных напряжений на этой плоскости
Г Л2 Г Л2
= С (а). (29)
Касательные напряжения т2уг и т^ и выражаемые через них а2п и тП достигают своих наибольших значений на контуре поперечного сечения. Рассматривая межслойные касательные напряжения т, в основном влияющие на разрушение, можно выразить по аналогии с (16) нормальное <5п и касательное тп напряжения:
с п Т п
п + п
V V Т* V
ап =Т^ 81П2а; Тп = Т^ с°з2а. (30)
После подстановки выражений (30) в критерий (29) можно найти угол а, соответствующий минимуму касательного напряжения т .
Для иллюстрации примем простую зависимость прочности от угла наклона площадки разрушения: С2 (а) = С (1 + б1п а), т.е. прочность
минимальна при межслойном разрушении (а = 0° ) и максимальна -при разрыве поперек волокон (а = 90°). Если принять для примера
Г т 2
I — I = 2, нахождение наклона площадки разрушения а сведется к
V о
, - 1 + б1п а „ д/ _
минимизации функции: / = т С-------------. Условие -^—= 0 приводит
1 + б1п 2а да
к уравнению 12s4 +16^3 - 4^2 - 8^ +1 = 0, корень которого ^ = б1п а0 = 0,12 и а0 = 7°, что хорошо согласуется с наклоном площадки разрушения при кручении образцов из стеклопластика.
3.2. Схема образования кинка при сжатии
На рис. 8 представлена схема распространенного вида разрушения при сжатии композитов или древесины путем образования полосы сдвига, «кинка» [19] (термин, применяемый уже и в отечественной литературе без кавычек). Часто считают, что сдвиговое разрушение происходит под углом 45°, т.е. вдоль направления наибольших касательных напряжений, однако учет влияния нормальных и касательных напряжений на плоскости разрушения и зависимости прочности от угла может привести к любым углам наклона кинка.
Ниже рассмотрена простая энергетическая модель возникновения кинка под углом а к продольной оси. При повороте волокон в полосе сдвига шириной 8 считается, что их несущая способность сохраняется до тех пор, пока они не станут перпендикулярны границе кинка: предельное состояние показано на рис. 8, б. Работа постоянной приложенной силы Р = а-И - Ь (где Ь - толщина образца) на перемещении Д = 8 (1 -б1п а) затрачивается на работу разрушения матрицы, пропорциональную объему кинка и удельной работе разрушения С. Приравнивая эти работы - силы и разрушения, - получим
оЪИЬ (1 - sin а) = ОЪЬ
h
sin а
(31)
до
Из условия минимума напряжения — = 0 находится угол, под
5а
которым «энергетически выгодно» образование кинка: либо а = 90°, что соответствует смятию по торцу, либо sin а* = 2, а* = 30° (а не 45°, как традиционно считают по аналогии с пластичными металлами).
Рис. 8. Схема образования полосы сдвига -кинка - при сжатии
Основной эффект: кроме удовлетворения чисто познавательного интереса, определение угла кинка а* помогает выбрать допустимую длину рабочей части образца на кручение и на сжатие: > h ctg а*.
4. Изгиб слоистых балок: два вида разрушения
4.1. Нелинейный критерий межслойного разрушения
Испытания на изгиб коротких композитных балок стандартизованы и широко используются для определения условной межслойной сдвиговой прочности, которую принято считать равной наибольшему касательному напряжению в момент разрушения прямоугольных балок:
3P
Т0 =
4Ь • h
(32)
Однако из опытов эта прочность зависит от отношения пролета I к высоте балки И и не может надежно применяться для расчетов реальных изгибаемых композитных элементов. Для объяснения такого
эффекта и для получения расчетных зависимостей в [12] был предложен линейный критерий межслойного разрушения:
(33)
■ тт ху = с.
где т, с - экспериментально определяемые параметры материала.
Считаем, что межслойное разрушение (рис. 9, а), происходит на расстоянии у0 от нейтральной оси, где достигает наибольшего значения линейная комбинация касательных т и нормальных изгибных а х
напряжений, которые вследствие неоднородности композита также дают вклад в межслойное разрушение.
Наиболее опасным оказывается сечение под приложенной нагрузкой, где в балочном приближении распределения напряжений при трехточечном изгибе имеют вид
3Р
т ху г. 1.3
Ък
Л
3Р • I • у
ЪИ
(34)
дР
Координата у0 найдется из условия — = 0 при подстановке (34)
ду
в (33):
Уо =
2т
(35)
Поскольку разрушение должно начинаться внутри образца, координата у0 < И/2, а следовательно, критерий (33) справедлив, пока
¡/к < т.
После подстановки (35) в (34) и в (33) с заменой критической силы Р на условную сдвиговую прочность (32) получаем согласующуюся с экспериментальными данными (рис. 9) зависимость
с • т
2 ¡2 т +—г
к2
или для удобства обработки эксперимента - линейную зависимость
1 т 1 (I'2
= -+-------1т1 . (37)
т0 с т • с V И )
На рис. 9, б приведена обработка экспериментальных данных по
изгибу стеклопластиковых балок в осях — -1 — | , и видно, что в пре-
т0 V И )
— Л2
1 I 2
делах 1 — 1 < т , в которых справедлив критерий (33), данные хорошо
V П )
ложатся на прямую линию (37).
Константу с можно считать прочностью по нормальным напряжениям при изгибе длинных балок; отношение 1/И = т разделяет области разрушения от расслоения и от растяжения. Как видно из рис. 9, б (светлые кружки), при 1/И > т прямая линия в координатах — -^-И
с • И
т0 VИ
21
(38)
хорошо согласуется с экспериментом. Зависимости (36) и (38) стыкуются при 1/И = т. Более того, в этой точке совпадают их производные, эти две кривые касаются друг друга, происходит плавный переход одной зависимости в другую. Гипербола (38) остается всегда выше, чем кривая (36). Это удобно при чисто компьютерной обработке экспериментов. Задав параметрически кривые (38) - для разрушения от нормальных напряжений и (36) - для межслойного разрушения, можно, меняя параметр с, «спускать» кривую (38) вплоть до касания с кривой (36). Точка касания определит значения параметров т и с. Предложенный подход дает более точное описание зависимости межслойной прочности от размеров балки, чем традиционный критерий т0 = т0 (штриховая линия на рис. 9, а), и позволяет рассчитывать несущую способность композитных элементов, подверженных действию комбинации изгибающих моментов и перерезывающих сил.
т 0 =
1/Тд, МПа-0,03
,-і
0,01
/ V
/ .
у V
юс т2 2 J0 ^00
0 5 10« 15
(Щ
т
а б
Рис. 9. Схема испытаний и зависимости условной сдвиговой прочности т0 от отношения пролета балки / к его высоте к для однонаправленного стеклопластика (а) и те же
экспериментальные данные в координатах — -| — I (•) и — -( — | (о) (б)
То ^ к) То ^ к)
4.2. Межслойная прочность при циклическом изгибе
По аналогии с подходом, примененным в п. 1.3, для описания циклической межслойной прочности можно заменить параметры критерия (33) т и с на функции от числа циклов т(Ы), с(Ы). По результатам испытаний на циклический изгиб коротких балок из углепластика в [14] была показана параллельность предельных прямых в коорди-
1 (/Л2
натах ——г-1 — I (рис. 10, б), что говорит о независимости
То N I к )
произведения т(Ы )• с(Ы ) = т(і)- с(і) от числа циклов. В качестве кривой Велера для изгибной прочности с (Ы) можно принять полученное при растяжении тех же углепластиков уравнение
с (Ы ) = с (і) (1,04- 0,035^ N). (39)
В результате изменение условной межслойной прочности т0 (32)
в зависимости от числа N циклов и от длины образца описывается на основании лишь результатов статических испытаний и независимых экспериментов на циклическое растяжение:
, . с (1) т (1)
Т0 (,//к) =-------------------------------------------------— (40)
т2 (1)(1,04 - 0,035^ N)-2 +1 —)
Обработка данных по формуле (40) приведена на рис. 10, а, расхождение с экспериментом не превышает 7 %. Возможна и обратная процедура - оценка кривой усталости для растяжения по данным циклических испытаний коротких образцов, что приводит к экономии материала и снижению необходимой мощности испытательного оборудования, а также упрощению крепления образцов, в том числе при температурных испытаниях.
т0, МПа
а(0), МПа
1 -1 / И = 4; 2 -1 / И = 5; 3 -1 / И = 7,5;
4 -1 /И = 9,5; 5-1 /И = 12,5 о - кривая усталости для однонаправленного углепластика при растяжении
1 - N = 100,2 -104,3 -105,4 -106 циклов
Рис. 10. Зависимости условной межслойной прочности от числа циклов: — (а)
Т0
от ^И] при различном числе циклов N (б)
Основной эффект: зависимость условной сдвиговой прочности от пролета балки для образцов данной толщины ставит под сомнение корректность ее определения при стандартных испытаниях на изгиб коротких балок с использованием формулы (32). Критерий, учитывающий совместное влияние нормальных и касательных напряжений, позволяет описать этот эффект и получить гладкую зависимость критической нагрузки от длины балки. Для изгиба криволинейных образцов или сегментов колец подобный подход описан в [10].
Заключение
Критерии прочности, учитывающие направленный характер разрушения волокнистых композитов, позволяют более правильно интерпретировать экспериментальные данные и рассчитывать на прочность композитные элементы конструкций.
Работа выполнена при финансовой поддержке РФФИ (грант № 12-08-00259-а).
Библиографический список
1. Алфутов Н.А., Зиновьев П. А., Попов В.Г. Расчет многослойных пластин и оболочек из композиционных материалов // Библиотека расчетчика. - М.: Машиностроение, 1984. - 264 с.
2. Ашкенази Е.К. Анизотропия машиностроительных материалов. -Л.: Машиностроение, 1969. - 111 с.
3. Сопротивление стеклопластиков / В.Л. Бажанов, И.И. Голь-денблат, В. А. Копнов, А. Д. Поспелов, А.М. Синюков. - М.: Химия, 1968. - 304 с.
4. Лехницкий С.Г. Кручение анизотропных и неоднородных стержней. - М.: Наука, 1971. - 240 с.
5. Малмейстер А.К., Тамуж В.П., Тетерс Г.А. Сопротивление жестких полимерных материалов. - Рига: Зинатне, 1967. - 398 с.
6. Мартиросян М.М. Упрочнение ориентированного стеклопластика при двухосном растяжении // Механика полимеров. - 1976. -№ 6. - С. 1025-1029.
7. Махмутов И.М., Полилов А.Н. Критерий прочности для ортогонально армированных композитов с полимерной матрицей // Изв. АН СССР. Механика твердого тела. - 1976. - № 5. - С. 160-169.
8. Полилов А.Н. Неупругие свойства композиционных материалов // Новое в зарубежной науке. Механика. Вып. 16. - М.: Мир, 1978. - 296 с.
9. Полилов А.Н. Критерий разрушения поверхности раздела в однонаправленных композитах // Изв. АН СССР. Механика твердого тела. - 1978. - № 2. - С. 116-119.
10. Полилов А.Н. Определение прочности при изгибе криволинейных образцов // Машиноведение. - 1984. - № 1. - С. 54-60.
11. Полилов А.Н. Схема расчета прочности косоугольно армированных композитов при плоском напряженном состоянии // Механика композитных материалов. - 1980. - № 2. - С. 221-226.
12. Полилов А.Н., Хохлов В.К. Расчетный критерий прочности композитных балок при изгибе // Машиноведение. - 1979. - № 2. - С. 53-57.
13. Исследование прочности намоточных эпоксидных угле- и стеклопластиков при кручении, растяжении и поперечном изгибе / Ю.Н. Работнов, И.Н. Данилова, А.Н. Полилов, Т.В. Соколова,
И.С. Карпейкин, М.В. Вайнберг // Механика полимеров. - 1978. -№ 2. - С. 219-225.
14. Критерий межслойной прочности углепластиков при циклических нагрузках / Ю.Н. Работнов, В.П. Когаев, А.Н. Полилов,
B.Б. Стрекалов // Механика композитных материалов. - 1982. - № 6. -
C. 983-986.
15. Усталостная прочность однонаправленных углепластиков при растяжении под углом к направлению армирования / Ю.Н. Работнов, В.П. Когаев, А.Н. Полилов, В.Б. Стрекалов, А.М. Думанский // Механика композитных материалов. - 1985. - № 1. - С. 44-47.
16. Lamination theory in the prediction of failure envelopes for filament wound materials subjected to biaxial loading / G.C. Eckold, D. Lead-better, P.D. Soden, P.R. Griggs // Composites. - 1978. - Vol. 9, no. 4. -P. 243-246.
17. Hashin Z., Rotem A., A fatigue failure criterion for fiber reinforced materials // Journal of Composite Materials. - 1973. - Vol. 7, no. 5. -P. 448-464.
18. Narayanaswami R. Evaluation of the tensor polynomial and Hoffman strength theories for composite materials // Journal of Composite Materials. - 1977. - Vol. 11, no. 5. - P. 366-377.
19. Parry T.V., Wroski A.S. Kinking and tensile, compressive and interlaminar shear failure in carbon-fiber-reinforced plastic beams tested in flexure // J.Mater.Science. - 1981. - Vol. 16, no. 5. - P. 439-450.
20. The strength of a filament wound composites under biaxial loading / P.D. Soden, D. Leadbetter, P.R. Griggs, G.C. Eckold // Composites, 1978. - Vol. 9, no. 4. - P. 247-250.
21. Tsai S.W., Wu E.M. A general theory of strength for anisotropic materials // J. Composite Materials. - 1971. - Vol. 5, no. 1. - P. 58-67.
22. Полилов А.Н., Татусь Н.А. Критерии прочности полимерных волокнистых композитов, описывающие некоторые экспериментально наблюдаемые эффекты // Проблемы машиностроения и автоматизации. -2008. - № 3. - С. 103-110.
References
1. Alfutov N.A., Zinoviev P.A., and Popov V.G. Raschet mnogo-sloynykh plastin i obolochek iz kompozitsionnykh materialov [Designof multilayered composite plates and shells]. Biblioteka raschetchika. Moscow: Mashinostroenie, 1984, 264 s.
2. Ashkenazi E.K. Anizotropiya mashinostroitelnykh materialov [Anisotropy of engineering materials]. Leningrad: Mashinostroenie, 1969, 111 p.
3. Bazhanov V.L., Goldenblat I.I., Kopnov V.A. Pospelov, A.D., Sinyukov A.M. Soprotivlenie stekloplastikov [Strengthofglass-fiber-rein-forced plastics]. Moscow: Khimiya, 1968, 304 p.
4. Lekhnitskii S.G. Kruchenie anizotropnykh i neodnorodnykh sterzhney [Torsionofanisotropic and inhomogeneous rods]. Moscow: Nauka, 1971, 240 p.
5. Malmeyster A.K., Tamuzh V.P., Teters G.A. Soprotivlenie zhest-kikh polimernykh materialov [Strength of rigid polymers]. Riga: Zinatne, 1967, 398 p.
6. Martirosyan M.M. Uprochnenie orientirovannogo stekloplastika pri dvukhosnom rastyazhenii [Strengthening of oriented glass-fiber-reinforced plasticsunder biaxial tension]. Mekhanikapolimerov, 1976, no. 6, p. 10251029.
7. Mahmutov I.M., Polilov A.N. Kriteriy prochnosti dlya ortogonalno armirovannykh kompozitov s polimernoy matritsey [Strength criterion for cross-ply reinforced composites with polymeric matrix]. Izv. AN SSSR, Mekhanika tverdogo tela, 1976, no. 5, pp. 160-169.
8. Neuprugie svoystva kompozitsionnykh materialov [Nonelastic properties of composite materials]. Novoe v zarubezhnoy nauke. Seriya: Mekhanika. Vyp. 16. Moscow: Mir, 1978, 296 p.
9. Polilov A.N. Kriteriy razrusheniya poverkhnosti razdela v od-nonapravlennykh kompozitakh [Fracture criterion for interface in unidirectional composites]. Izv. AN SSSR, Mekhanika tverdogo tela, 1978, no. 2, p. 116-119.
10. Polilov A.N. Opredelenie prochnosti pri izgibe krivolineynykh obraztsov [Strength estimation using curvilinear beam bending]. Mashinov-edenie, 1984, no. 1, pp. 54-60.
11. Polilov A.N. Skhema rascheta prochnosti kosougolno armirovannykh kompozitov pri ploskom napryazhennom sostoyanii [Model for strength estimationfor angle-ply composites under plane stress state]. Mek-hanika kompozitnykh materialov, 1980, no. 2, pp. 221-226.
12. Polilov A.N. Khokhlov V.K. Raschetnyy kriteriy prochnosti kompozitnykh balok pri izgibe [Design strength criterion for composite beams under bending]. Mashinovedenie, 1979, no. 2, pp. 53-57.
13. Polilov A.N. Khokhlov V.K. Issledovanie prochnosti namotoch-nykh epoksidnykh ugle- i stekloplastikov pri kruchenii, rastyazhenii i pop-erechnom izgibe [Design strength criterion for composite beams under bending]. Mekhanika polimerov, 1978, no. 2, pp. 219-225.
14. RabotnovYu.N., Kogaev V.P., Polilov A.N., Strekalov V.B. Kriteriy mezhsloynoy prochnosti ugleplastikov pri tsiklicheskikh nagruz-kakh [Interlaminar strength criterionfor CFRP under cyclic loading]. Mekhanika kompozitnykh materialov, 1982, no. 6, pp. 983-986.
15. Rabotnov Y.N., Kogaev V.P., Polilov A.N., Strekalov V.B., Du-mansky A.M. Ustalostnaya prochnost odnonapravlennykh ugleplastikov pri rastyazhenii pod uglom k napravleniyu armirovaniya [Fatigue strength of unidirectional carbon fiber reinforced plastics in tension at an angle to the direction of reinforcement]. Mekhanika kompozitnykh materialov, 1985, no. 1, pp. 44-47.
16. Eckold G.C., Leadbetter D., Soden P.D., Griggs P.R. Lamination theory in the prediction of failure envelopes for filament wound materials subjected to biaxial loading. Composites, 1978, vol. 9, no. 4, pp. 243-246.
17. Hashin Z., Rotem A., A fatigue failure criterion for fiber reinforced materials. Journal of Composite Materials, 1973, vol. 7, no. 5, pp. 448-464.
18. Narayanaswami R. Evaluation of the tensor polynomial and Hoffman strength theories for composite materials. Journal of Composite Materials, 1977, vol. 11, no. 5, pp. 366-377.
19. Parry T.V., Wroski A.S. Kinking and tensile, compressive and interlaminar shear failure in carbon-fiber-reinforced plastic beams tested in flexure. J.Mater.Science, 1981, vol. 16, no.5, pp. 439-450.
20. Soden P.D., Leadbetter D., Griggs P.R., Eckold G.C. The strength of a filament wound composites under biaxial loading. Composites, 1978, vol. 9, no. 4, pp. 247-250.
21. Tsai S.W., Wu E.M. A general theory of strength for anisotropic materials. J. Composite Materials, 1971, vol. 5, no. 1, pp. 58-67.
22. Polilov A.N., Tatus N.A. Kriterii prochnosti polimernykh vo-loknistykh kompozitov, opisyvayushchie nekotorye eksperimentalno na-blyudaemye effekty [Strength criteria forfiber reinforced plastics described some of the experimentally observed effects]. Problemy mashinostroeniya i avtomatizatsii, 2008, no. 3, pp. 103-110.
Об авторах
Полилов Александр Николаевич (Москва, Россия) - доктор технических наук, профессор, заведующий лабораторией безопасности и прочности композитных конструкций Института машиноведения им. А.А. Благонравова Российской академии наук (101990, Москва, М. Харитоньевский пер. 4., e-mail: polilov@imash.ru).
Татусь Николай Алексеевич (Москва, Россия) - кандидат технических наук, научный сотрудник лаборатории безопасности и прочности композитных конструкций Института машиноведения им. А.А. Благонравова Российской академии наук (101990, Москва, М.Харито-ньевский пер. 4., e-mail: nikalet@mail.ru).
About the authors
Polilov Alexander Nikolaevich (Moscow, Russian Federation) -Doctor of Technical Sciences, Professor. The Head of Laboratory of Safety and Strength of Composite Structures Institute of Machines Science of Russian Academy of Sciences (101990, Moscow, Russian Federation, Malyi Kharitonyevsky per. 4, e-mail: polilov@imash.ru).
Tatus Nikolay Alexeevich (Moscow, Russian Federation) - Ph. D. in Technical Sciences, Researcher of Laboratory of Safety and Strength of Composite Structures, Institute of Machines Science of Russian Academy of Sciences, Moscow (101990, Moscow, Russian Federation, Malyi Kharitonyevsky per. 4, e-mail: nikalet@mail.ru).
Получено 15.05.2012