Development of the Views of Yu.N. Rabotnov on Strength Criteria of Composites

Cover Page

Cite item

Full Text

Abstract

Besides the well-known fundamental works of Yu.N. Rabotnov in the field of hereditary elasticity and creep theory, one of the aspects of his scientific activity was the mechanics of composite materials and, in particular, a new class of strength criteria for structurally anisotropic composites proposed by him.

The main feature of Yuri Rabotnov’s approach was not an attempt to construct a uniform smooth limiting surface in the stress-space, but taking into account real fracture mechanisms, which, as a rule, are directional in nature.

Now such approaches become crucial in calculation algorithms modeling the fracture process with taking into account the degradation of elastic and strength properties, but in the period of the first Yu.N. Rabotnov’s papers they were pioneering and caused certain discussions. Development and application of some of Rabotnov’s proposed types of strength criteria for fiber-reinforced composites in tension, compression and complex stress state are discussed in this anniversary paper.

Full Text

  1. Введение. Академик Юрий Николаевич Работнов, безусловно, относится к плеяде выдающихся ученых-механиков ХХ столетия [1]. К наиболее известным областям, в которых Ю.Н. Работновым были получены классические фундаментальные результаты, относятся техническая теория оболочек (с анализом краевого эффекта), теория наследственной упругости (1948 г.) [2–4] в приложении к проблеме ползучести и механика повреждаемости, введение параметра поврежденности [5]. Полученные им результаты мирового уровня привели к тому, что Такео Ёкобори, создавая Международный конгресс по разрушению (ICF), пригласил Работнова на должность одного из директоров Конгресса в ранге вице-президента. В конце жизни Юрий Николаевич получил выкованное на медном листе свидетельство признания его Почетным президентом этого Конгресса.

После одной из конференций Конгресса в Мюнхене (1972 г.) Юрий Николаевич окончательно увлекся механикой разрушения [6], и на Международной конференции в МГУ заявил, что не теории упругости, пластичности, ползучести, а именно – механика разрушения и механика композитов являются наиболее интересными и перспективными областями механики твердого тела. Тогда же по инициативе Работнова начался перевод и издание фундаментальных коллективных монографий [7, 8], сыгравших огромную роль в развитии этих областей науки в нашей стране.

Нельзя не отметить педагогическую ценность его учебников [9, 10], первый из которых Валентин Валентинович Новожилов в приветственном письме назвал “жемчужиной среди многочисленных учебников по сопротивлению материалов”, а второй по праву можно считать лучшей (вероятно, в мире) энциклопедической монографией по механике твердого тела.

Возможно, это счастливая случайность, но привезенная С.Т. Милейко из Англии книга [11] чрезвычайно понравилась Юрию Николаевичу и была мгновенно переведена и издана. Его заинтересовала задача Гордона об остановке трещины границей раздела [12], и он включил полученный результат в монографию [10]. А далее – развивался интерес к особым механизмам расслоения и расщепления композитов [13, 14], и последняя его прижизненная статья [15] была посвящена новой модели множественного расщепления с выпучиванием, приводящим сжимаемое звено бамбука (или композитную трубу) к форме, напоминающей “китайский фонарик”.

Оригинальный взгляд имел Юрий Николаевич и на методы построения предельных поверхностей в пространстве напряжений, приводя иногда строчку Павла Когана “Я с детства не любил овал, я с детства – угол рисовал”. Работнов обосновал новый тип критериев прочности [16] (разд. 1–3), отличающийся от популярных в то время тензорно-полиномиальных критериев, и с обоснованием новых критериев он выступал и на Съезде по механике в Киеве [17], и на Конгрессе в Канаде [18], и на советско-японском симпозиуме в МГУ [19]. Это была еще одна область, которая занимала его в последние годы – механика разрушения композитов, обоснование особых видов критериев прочности конструктивно анизотропных материалов [20–23].

Последняя книга Юрия Николаевича по механике разрушения [24] была издана (и позднее – переиздана) на основе перевода его лекций с французского, которым он блестяще владел еще с детства. В 1979 г. Юрий Николаевич читал во Франции по приглашению университета Пуатье цикл лекций по механике разрушения металлов и композитов. Лекции были столь увлекательными и блестящими, что, по воспоминаниям очевидцев, слушатели (большей частью профессора, а не только студенты) часами не отпускали лектора, засыпая его бесконечными вопросами и просьбами рассказать еще, на которые Юрий Николаевич откликался, не считаясь со временем. Возможно, именно эти перегрузки и привели в скором времени к тяжелому удару, от которого он так и не смог полностью оправиться.

Книгу научного наследия [25] мы составляли только из статей, написанных его рукой (без соавторов), и поэтому все мы можем вновь и вновь насладиться не только глубокими идеями, но и его блестящим русским языком.

Цель данной юбилейной статьи – изложить некоторые идеи Ю.Н. Работнова в области механики разрушения композитов и показать их развитие в современной научной литературе.

  1. Предельные поверхности для совместно работающих упруго-пластических структур. Многие известные авторы (Мизес–Хилл, Малмейстер, Цай-Ву, Гольденблат–Копнов и др.) формулировали критерии прочности анизотропных материалов в виде полинома от инвариантов, представляющих собой свертки тензора напряжений с “тензорами прочности” четных рангов: I1=Arsσrs,I2=Bijklσijσkl,I=Cijklmnσijσklσmn и т. д.

Если ограничиться вторым инвариантом, критерий приобретает вид:

I1α+I2β+I3γ+...=1 (2.1)

Число констант Ars даже для плоского случая равно трем, Bijkl – шести (в трехмерном случае 6 и 14, соответственно), и нахождение такого большого числа констант прочности из ограниченного числа опытов оказывается неустойчивым: экспериментальные данные описываются одинаково хорошо совершенно различным набором констант, и некоторые из них по произволу приходится назначать равными нулю.

Параметры полиномиальных критериев типа (2.1) можно было бы считать просто коэффициентами регрессионного анализа, подбираемыми по критерию наилучшего соответствия доступному набору экспериментальных данных. Тогда, казалось бы, снимаются вопросы: какие параметры полагать равными нулю, взаимодействием между какими напряжениями можно пренебречь, и просто нужно отбрасывать некоторые прочностные параметры за их малостью. Но основное противоречие при этом не исчезает. Определяемые параметры остаются, по сути, не свойствами материала, а функционалами от имеющихся экспериментальных данных и метода их обработки, и определение большого числа параметров остается неустойчивым.

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

Альтернативный подход к построению критериев разрушения был предложен Ю.Н. Работновым [10, 16, 24] на основе модели совместно работающих упруго-пластических структур, в данном случае – волокон и матрицы. Приближение предельной кривой в плоскости σ1σ2 эллипсом для однонаправленного (рис. 1, а) или ортогонально армированного (рис. 1, б) пластика со слабой матрицей вряд ли может считаться удовлетворительным. Действительно, пренебрегая прочностью матрицы, мы получаем ортогональную решетку, которая может выдерживать только напряжения вдоль волокон, но не под углом к ним, и предельная поверхность в координатах σ1σ2τ12 представляет собой прямоугольник в плоскости σ1σ2. Учет ненулевого предела текучести матрицы превращает плоский прямоугольник в “диванную подушку” (рис. 1, б). Разумеется, вводя в уравнение (2.1) инварианты более высоких порядков, можно эту поверхность с плоскими гранями аппроксимировать с любой степенью точности, но число констант и соответствующие трудности их экспериментального определения при этом неизмеримо возрастут. Отсылая читателя за математическими формулировками к статьям [10, 16, 24, 25] Ю.Н. Работнова, можно для волокнистых композитов пояснить геометрию построения сложной предельной поверхности, которая получается движением центра предельной поверхности для одного компонента по предельной поверхности для другого компонента.

 

Рис. 1. Предельные поверхности в виде: а – закругленного цилиндра для однонаправленного композита, б – “диванной подушки” (ottoman cushion) для ортогонально армированного пластика; в — трехмерного косоугольного параллелепипеда для армирования в трех и более направлениях.

 

Например, предельная поверхность для ряда параллельных волокон имеет вид отрезка вдоль оси σ1 (рис. 1, а). Добавляя предельную поверхность для изотропной матрицы в виде эллипсоида, получаем предельную поверхность для однонаправленного пластика в виде цилиндра с закругленными концами. Добавление ортогонального семейства волокон, для которого предельная “поверхность” представляет собой отрезок вдоль оси σ2, эквивалентно движению цилиндра вдоль этой оси и превращению поверхности в “диванную подушку” (рис. 1, а) с плоскими гранями и закругленными краями. Добавление третьего семейства волокон, то есть движение этой “подушки” вдоль третьего направления приведет к поверхности в виде существенно трехмерного косоугольного параллелепипеда с закругленными ребрами (рис. 1, в). Такая структура армирования (с тремя и более семействами волокон) лишена главного недостатка ортогонально армированных и однонаправленных пластиков – слабого сопротивления сдвигу и растяжению под углом к волокнам.

Применительно к представленным на рис. 1 предельным поверхностям легко вывести, например, зависимость критического напряжения σ(φ) от угла φ между направлением растяжения и направлением укладки волокон (ось 1). Проекции “вектора” нагружения в координатах σ1σ2τ12 выражаются в виде:

σ1=σ(φ)cos2φ,σ2=σ(φ)sin2φ,τ12=σ(φ)sinφcosφ (2.2)

Термин “вектор” используется для наглядности и взят в кавычки, так как компоненты тензора второго ранга преобразуются при повороте системы координат не по векторному правилу. Этот “вектор” (2.2) может “протыкать” боковую поверхность цилиндра на рис. 1, а, уравнение которой имеет вид:

σ2σ(90)2+τ12τ(0)2=1, (2.3)

где σ(90),τ(0) – прочности при растяжении ортогонально волокнам и на сдвиг вдоль волокон. Условие такого “протыкания” при подстановке (2.2) в (2.3) определяет выражение критерия прочности по условию разрушения матрицы при растяжении под углом к волокнам:

σφ=sin4φ/σ(90)2+sin2φcos2φ/τ(0)21/2 (2.4)

Если для наглядности принять прочности матрицы на сдвиг и на растяжение одинаковыми: σ(90)=τ(0)=τ, то из (2.4) получится совершенно простая зависимость, которая подтверждалась в экспериментах:

σ(φ)=τsinφ (2.5)

“Вектор” нагружения может протыкать также сферическую (при σ(90)=τ(0)=τ) концевую часть цилиндра, уравнение которой имеет вид:

σ1σ(0)2+σ22+τ122=τ, (2.6)

где σ(0)(τ) – прочность при растяжении вдоль волокон. Подстановкой (2.2) получим из (2.6) квадратное уравнение для определения σ(φ):

Aσ2(φ)2Bσ(φ)+C=0, (2.7)

где A=c2+s4, B=c2σ(0), C=σ2(0)τ2, c=cosφ, s=sinφ.

Решение уравнения (1.7) получается довольно громоздким, но поскольку разрыв волокон происходит лишь при весьма малых углах φ, легко показать, что σ(φ)σ(0) при φ0, и для малых углов φ можно с достаточной точностью считать:

σ(φ)σ(0) (2.8)

Равенство напряжений (2.5) и (2.8) определяет угол φ* перехода от разрыва волокон к разрушению матрицы, и можно получить удобную оценку рационального угла малой разориентации: φarcsinσ(90)/σ(0), при котором продольная прочность практически не снижается, а сопротивление расщеплению резко растет. Подобная разориентация волокон создается в стволе дерева, и поэтому свилеватая древесина, оставаясь прочной в продольном направлении, имеет гораздо более высокую, чем прямослойная древесина, прочность в поперечном направлении, и поэтому ее очень трудно расколоть.

Для ортогонально армированного пластика (рис. 1, б) “вектор” нагружения может “протыкать” предельную поверхность либо по плоскости τ12=τ(0), и тогда зависимость прочности имеет вид

σφ=τ0/sinφcosφ, (2.9)

либо по скругленному краю, и тогда приближенно можно считать (в зависимости от стороны “подушки”):

σ(φ)σ(0) или σ(φ)σ(90) (2.10)

Угол φ* смены площадки разрушения определяется равенством выражений (2.9) и (2.10):

φ12arcsin2τ0/σ0 (2.11)

Полученные на основе подхода Ю.Н. Работнова выражения критериев прочности однонаправленных и ортогонально армированных пластиков оказываются чрезвычайно простыми по сравнению с тензорно-полиномиальными критериями и достаточно хорошо согласуются с экспериментами по определению прочности в различных направлениях.

  1. Критерии прочности типа Кулона – Мора. Ю.Н. Работнов предложил [17–19] формулировать критерии прочности, учитывающие направленный характер разрушения волокнистых композитов, в виде предельного соотношения между нормальным σn и касательным τn напряжениями на площадке возможного разрушения, содержащей направление одного из семейств волокон:

Φσn,τn=1 (3.1)

Простейшее предположение [24, 26] о виде функции Φ в (3.1) состоит в том, что она линейна, и тогда, например, для ортогонально армированного пластика можно выписать два критерия на одной (n1) или на другой (n1) плоскости возможного разрушения:

σni+miτni=ci;i=1,2,c1=σ0,c2=σ90 (3.2)

Используя формулы (2.2) преобразования тензора напряжений при повороте координатных осей на угол φ, получаем из (3.2) для случая одноосного растяжения две зависимости прочности σ(φ) от угла вырезки образца φ:

σφ=c1cos2φ+m1sinφcosφ при φφ (3.3)

σφ=c2sin2φ+m2sinφcosφ при φ>φ (3.4)

При угле φ* происходит переход от одной плоскости разрушения к другой, и этот угол находится из равенства (3.3) и (3.4), приводящего к квадратному уравнению относительно tgφ*, из которого находится:

φ=arctg12m1εm2+m2m1ε2+4εarctgεm2ε1;ε=c2c1 (3.5)

Зависимости (3.3) и (3.4) представлены на рис. 2 для разных соотношений между числами волокон в ортогональных направлениях и, в частности, для однонаправленного композита, для которого сохраняется и оказывается наиболее естественным тот же подход в предположении, что число волокон в ортогональном направлении равно нулю. Из (3.5) при малых углах φ<φ происходит разрыв волокон, и зависимость прочности имеет вид (3.3). При φφ происходит разрушение матрицы по критерию (3.4).

 

Рис. 2. Схема нагружения (а) и зависимости прочности ортогонально армированных композитов от направления вырезки образцов: б – однонаправленный (0); в – “равнопрочный” (0/90); г – (04/901). Линия 1 – (2.3); 2 – (2.4).

 

Для однонаправленных композитов угол φ очень мал, например, для углепластиков: σ(0)1000 МПа, σ(90)25÷40 МПа, m20,5 φ=3÷4. Надежное экспериментальное определение m1 в таком крайне узком диапазоне углов вряд ли возможно, поэтому для инженерных приложений удобно упростить зависимость (3.3), полагая σφ=σ0 при φφ. Погрешность при таком вычислении σ(φ) по сравнению с (3.3) в диапазоне φ<φ не превышает 4%, что вполне допустимо с учетом точности эксперимента и разброса свойств композитов.

На рис. 3 проиллюстрирован метод определения параметров линейных критериев прочности: экспериментальные данные для ортогонально армированного углепластика нанесены в координатах σn=σ*(φ)cos2φ, τn=σ*(φ)sinφcosφ.

 

Рис. 3. Предельные кривые для углепластика (0/90) в координатах σn-τn.

 

Прямая линия (3.2) в этих координатах одинаково хорошо согласуется с экспериментом, как и участок эллипса:

σnσ(0)2+τnτ2=1, (3.6)

соответствующий следующей зависимости прочности:

σφ=cos4φ/σ20+sin2φcos2φ/τ201/2 (3.7)

Но проведение прямой линии через экспериментальные точки в координатах σnτn всегда нагляднее и надежнее, поэтому предпочтение следует отдать линейным критериям (3.2), приводящим к расчетным зависимостям (3.3), (3.4). Напомним, что критическое соотношение между нормальным и касательным напряжением в теориях Кулона и Мора тоже, как правило, предполагалось линейным, но в теории Мора угол площадки разрушения требует определения, а для композитов эта площадка должна содержать направление одного из семейств волокон, и она определяется простым перебором.

Естественно считать φ (3.5) оптимальным углом малой разориентации волокон (аналог свилеватой, и потому нерасщепляющейся древесины). При таком угле разориентации прочность вдоль волокон практически не снижается, а сопротивление росту продольных трещин возрастает многократно, исключая расщепления композитного материала при циклических нагрузках и вблизи концентраторов напряжений. Подобным образом поступает и природа: древесные волокна растут не строго вдоль ствола, а совершают полный оборот вокруг ствола на длине примерно 6–10 м. При этом угол разориентации (3°÷5°) оказывается весьма близким к полученным выше оценкам.

Простота критерия (3.2) позволяет легко распространить его на область усталостного разрушения [20, 21], полагая параметры критерия функционалами от истории нагружения или – при постоянной амплитуде – функциями от числа циклов.

С помощью критериев, учитывающих направленный характер разрушения, можно объяснить и описать немонотонные и негладкие поверхности прочности: в простейшем случае – зависимости прочности от направления нагружения. Традиционные тензорно-полиномиальные критерии описывают гладкую предельную поверхность в пространстве напряжений (многомерный эллипсоид), поэтому сложно установить условия возникновения того или иного вида разрушения, что является известным недостатком для послойного анализа, когда в зависимости от вида разрушения монослоя надо либо полностью исключать его из рассмотрения (после разрыва волокон), либо полагать практически нулем характеристики матрицы (после выполнения условия разрушения поверхности раздела). Сейчас подобные критерии Хашина с установлением коэффициентов деградации свойств в зависимости от вида разрушения считаются общепризнанными, но Ю.Н. Работнов обратил на это внимание еще 50 лет назад.

  1. Модель ромба для косоугольно армированных труб. Ю.Н. Работнов предложил [23] модель разрушения композитной ячейки в виде ромба из нерастяжимых нитей (рис. 4), повороту которых препятствует момент M=2μl2sinαcosα, возникающий из-за сдвига матрицы и пропорциональный площади ромба.

 

Рис. 4. Схема деформирования элемента косоугольно армированного композита.

 

Сила осевого растяжения, приложенная к вершине ромба, пропорциональна длине малой диагонали Pz=2σzlsinα. Толщину ромба можно принять за единицу: она сокращается в конечном уравнении энергетического баланса. Смещение вершины при изменении угла: Δ=2lcosα'dα=2lsinαdα. Знак “минус” означает, что с ростом осевых напряжений σz угол α убывает. Приравнивая работу осевой силы 12PzΔ работе момента Mdα, препятствующего изменению угла dα, получим

σzα=μctgα=σz45ctgα=σz45/t, (4.1)

где критическое значение распределенного момента μ соответствует прочности на растяжение трубы с намоткой ±45°σz(45)=μ. Удивительно, но модель Ю.Н. Работнова привела в точности к такому же соотношению (4.1), полученному ранее в [22] чисто феноменологически при обработке экспериментов по растяжению стеклопластиковых трубок.

Наиболее интересно применение данной модели для двухосного растяжения. Окружная сила Pθ=2σθlcosα от действия окружных напряжений σθ стремится, напротив, увеличить α, и приведенный выше анализ остается в силе с точностью до замены α на 90°α. Как видно на рис. 4, осевые σz и окружные σθ напряжения пытаются исказить ромб в разных направлениях, что дает (рис. 5) значительный упрочняющий эффект и приводит к следующему условию разрушения:

σztgασθctgα=σz45 (4.2)

Прямые скобки – обозначение абсолютной величины – говорят о смене знака выражения в скобках в зависимости от того, какой член в этой разности больше, что приводит на плоскости σθσz к двум параллельным предельным прямым линиям 1, 2 (рис. 5), проходящим через точки σzασθα на осях.

Для заданного отношения приложенных напряжений, как показано в [26], легко указать оптимальный угол намотки, обращающий в ноль разность в прямых скобках в (3.2). Это означает невозможность перекашивания ромба и перенесение всех усилий только на волокна. Например, в цилиндрической части баллона для сжатого газа (“труба с заглушками”) окружные напряжения примерно вдвое превышают осевые σθ=2σz, и оптимальный угол намотки α* имеет известное из “нитяной аналогии” значение: tg2α=2α=arctg2=5444'.

Отрезок прямой 3 замыкает прямоугольную предельную поверхность, и при углах, близких к оптимальному, когда перекашивание ромба невозможно, происходит разрыв волокон, и реализуется наибольшая прочность. Критерий Работнова, приводящий к кусочно-линейной предельной поверхности, несмотря на свою удивительную простоту, хорошо согласуется с экспериментами (рис. 5) и позволяет легко проектировать оптимальную структуру армирования сосудов давления, намотанных симметричными парами слоев под углами ±αi.

 

Рис. 5. Предельные поверхности для двухосного нагружения стеклопластиковых труб: с намоткой: ±55° (а), ±72.5° (б), ±35° (в).

 

  1. Этюд об образовании кинка при сжатии. Обычно считают, что полосы сдвига (линии Чернова–Людерса) располагаются в плоском образце при сжатии (или растяжении) под углом 45° – по направлению наибольших касательных напряжений. Но Юрий Николаевич увидел, что в экспериментах на сжатие стеклопластиков возникает совершенно другой угол направленного разрушения. Он сформулировал изящную модель, результат которой до изумления точно совпал с экспериментом. Согласно примененному Работновым энергетическому критерию, кинк должен возникать (и возникает) под углом α=30° к продольной оси [26]. Правда, он решил эту задачу не публиковать, посчитав ее неким частным случаем.

Термин кинк (kink – полоса сдвига) уже и в отечественной литературе применяется без кавычек. Образование кинка представляет одну из новых и очень, как выяснилось, важных проблем расчета композитных конструкций. Например, в металлах при действии сжимающих напряжений рост трещин около отверстий не происходит, и надрезы, перпендикулярные направлению сжатия, считаются неопасными. Но в волокнистых композитах, именно при сжатии, возможно появление кинка около отверстий. Кинк – опасный вид растущей (линейной) зоны разрушения, но методы оценки несущей способности волокнистых композитов при сжатии обоснованы все еще недостаточно. И связано это в первую очередь со сложностью и многообразием механизмов разрушения.

Простая модель образования кинка при сжатии композитов (или древесины), проиллюстрированная на рис. 6, основана на следующих допущениях.

 

Рис. 6. Схема образования кинка при сжатии.

 

Считается, что при повороте волокон в кинке высотой δ, сопровождаемом разрушением матрицы, несущая способность сохраняется до тех пор, пока потерявшие устойчивость волокна не станут перпендикулярны прямолинейной границе кинка. Это предельное состояние показано на рис. 6, б. При образовании кинка происходит уменьшение высоты образца на величину Δ=δ(1sinα), и работа внешней приложенной силы сжатия P=σhw (где h и w – толщина и ширина образца) на этом перемещении затрачивается на работу разрушения матрицы, пропорциональную произведению объема кинка δhw/sinα на некоторую (неизвестную) удельную работу разрушения G. Приравнивая работу силы и работу разрушения, получим

σδhw1sinα=Gδwhsinα (5.1)

Из условия минимума напряжения в (5.1): σ/α=0[(1sinα)sinα]'=0 находится угол α*, под которым энергетически выгодно образование кинка: это либо α*=90°, что соответствует смятию по торцу, либо sinα=12α=30° – и это значение с высокой точностью подтверждается в опытах на сжатие тканых стеклопластиков.

Кроме удовлетворения чисто познавательного интереса, определение угла кинка позволяет выбрать допустимую длину L рабочей части образца на сжатие, чтобы разрушение происходило в пределах рабочей зоны: L>hctgα*1.7h. Модель, предложенная Работновым, лишь открывает путь к созданию расчетных методик, может быть, поэтому он решил, что ее еще рано публиковать.

  1. Замедленное разрушение композитов при сжатии в результате потери устойчивости волокон. Один из характерных видов разрушения волокнистых композитов при сжатии – это внутренняя потеря устойчивости волокон.

В нулевом приближении из уравнения продольного изгиба шарнирно опертой пластины толщиной h и длиной L (ширину полагаем равной единице, так как она во все соотношения входит линейно и сокращается в окончательных выражениях)

B2vx2+Pv=0,B=112Efh3,P=σh (6.1)

можно найти эйлерово напряжение при потере устойчивости, принимая прогиб в виде vx=v0sin(πx/L),

σe=π2Efh212L2 (6.2)

В композите без поддержки матрицы длинные и тонкие волокна теряли бы устойчивость при очень низких сжимающих напряжениях, поэтому для оценки прочности волокнистых композитов на сжатие Работновым принята модель композитной среды в виде чередующихся слоев из материалов волокон и матрицы, и слои волокон толщиной hf предполагаются идеально упругими с модулем Юнга Ef.

В первом приближении слои матрицы толщиной hm предполагаются подчиняющимися при сдвиге линейному закону Гука: τ=Gmγ. Считается, что вся сжимающая нагрузка P воспринимается упругими волокнами и напряжение в них: σf=σc1+hm/hf=σc/ψf, где σc – среднее напряжение в композите,  ψf – объемная доля волокон. Направив ось x по средней линии одного из слоев, рассмотрим возможность симметричного выпучивания слоев в одну сторону (рис. 7).

 

Рис. 7. Схема внутренней потери устойчивости волокон (слоев) при сжатии.

 

Находящиеся между ними слои более слабой матрицы получают сдвиговую деформацию γ=v(x)/x, где vx=v0sinλx – прогиб упругих слоев. В результате сдвига возникают касательные напряжения τ с двух сторон изогнутого слоя, действие которых препятствует выпучиванию и сводится к распределенному по длине моменту m=τhf=Gmγhf на единицу ширины. Распределенный момент действует с обратным знаком относительно момента от приложенной продольной силы Р, и уравнение изгиба (6.1) принимает вид:

B3vx3+Pvxm=0 (6.3)

После подстановки прогиба и сокращения на λcosλx получим Pe=Bλ2+Gmhf, и по аналогии с (6.2)

σce=112λ2ψfEfhf2+Gmψf (6.4)

Из (6.4) видно, что при линейном законе упругости для матрицы минимум критического напряжения обеспечивает наибольшая длина волны выпучивания L=π/λ, и основной вклад в эйлерово напряжение для композита дает модуль сдвига матрицы, так как он намного больше эйлерова напряжения (6.2) для волокон.

Во втором приближении принимается нелинейное соотношение между сдвиговыми деформациями и касательными напряжениями в матрице: τ=τ0γ1/m. При этом уравнение (6.3) принимает вид:

d3vdx3+k2dvdxβ1mdvdx1m=0, (6.5)

где k2=P/B, β1/m=τ0hf/B. Степень  в выражении для  выбрана для удобства интегрирования уравнения (6.5).

Чем больше прогиб слоев, тем выше сдвиговая деформация матрицы и возникающие в ней касательные напряжения. Поэтому можно амплитуду прогиба принять пропорциональной касательному напряжению: v=τsinλx.

Тогда уравнение (6.5) перепишется в виде:

τλcosλxk2λ2=βτλcosλx1/m (6.6)

Домножая обе части уравнения (6.6) на cosλxdλx и интегрируя по λx от π/2 до +π/2, получим, что интеграл от левой части (6.6) пропорционален π/2, а от правой – числовому коэффициенту – интегралу от дробной степени косинуса:

μ1m=2ππ2+π2cos1m+1λxd(λx). Тогда уравнение (6.6) примет более простой вид

k2=λ2+Kλ1+1/m;K=βμτ1m1/m (6.7)

Напомним, что  больше единицы, k2 пропорционален эйлеровой нагрузке Pe, и из (6.7) следует, что существует критическое значение длины волны выпучивания, реализующее минимум эйлеровой нагрузки:

λ*=K(m1)2mm3m1

Подставляя λ* в (6.7), можно найти минимальное эйлерово напряжение для критической длины волны при внутреннем выпучивании слоистой структуры.

В третьем приближении Работновым была рассмотрена модель длительной прочности (запаздывающего разрушения – delayed fracture). Предполагалось, что волокна нечувствительны ко времени нагружения, а деформирование матрицы при сдвиге подчиняется закону установившейся ползучести

τ=τ0γ˙γ01/n (6.8)

Здесь τ0 и γ0 – константы ползучести, одна из которых может быть выбрана по произволу. При симметричном выпучивании слои матрицы получают сдвиговую деформацию γ=v(x,t)/x, где vx,t – прогиб упругих слоев. Принимая за безразмерное время t истинное время, умноженное на γ0, получаем из (6.8)

m=τhf=τ0hf2vxt1/n, и уравнение (6.5) преобразуется к виду

3vx3+k2vxβ1/n2vxt1/n=0 (6.9)

После разделения переменных:

τnμβdτ=λn1k2λ2ndt, (6.10)

где обозначено μ1n=2ππ/2π/2(cosz)1n+1dz;  z=λx. Если в начальном состоянии ось сжимаемого слоя описывалась уравнением v0=τ0sinλx, то из (6.10):

(τ1nτ01n)(1n)1μβ=λn1k2λ2nt, и можно найти критическое время, при котором прогиб v и напряжение τ обращаются в бесконечность:

tc=μβλτ0n1k2λ2nn1 (6.11)

В реальной структуре композита всегда имеются начальные искривления волокон с разными длинами волн, поэтому для оценки долговечности следует выбрать такое значение λ*, которому соответствует минимальное критическое время из (6.11). Условие tc/λ=0 дает:

λ2=n13n1k2 (6.12)

Если обозначить L0 критическую длину шарнирно опертых упругих слоев (без матрицы), для которых данная сила P является эйлеровой, то длина L* полуволны потери устойчивости композита с матрицей, находящейся в условиях ползучести, дается формулой

L*=L03n1n1 (6.13)

Минимальное значение критического времени получается при постановке λ* из (6.12) в (6.11):

tmin=3n1BP3n12μβτ02nτ0nn1n+12 (6.14)

Предложенная Работновым довольно простая модель длительного разрушения композита при сжатии позволяет оценить влияние ползучести матрицы на критическое время (6.14) и предсказывает наличие характерной длины (6.13) волны внутреннего выпучивания волокон.

  1. Современное развитие идей Ю.Н. Работнова. Фундаментальные результаты Работнова в области механики композитов, в частности, изложенные в [10, 24] и в посвященной 100-летию Учителя книге [26] известны как в нашей стране, так и за рубежом. Его непосредственные ученики: академики И.Г. Горячева, Б.Д. Аннин, чл.- корр. РАН С.А. Шестериков, Е.В. Ломакин, доктора наук С.Т. Милейко, Ю.В. Немировский, Г.И. Брызгалин, А.В. Березин, Ю.В. Суворова, А.А. Мовчан, В.И. Астафьев, А.М. Думанский, А.Н. Полилов – передавали и передают идеи Юрия Николаевича ученикам второго поколения: С.И. Алексеева, И.В. Викторова, Н.А. Татусь, К.А. Хвостунков – всех не перечислишь.

Анализ литературы последних лет показывает огромное число ссылок на основополагающие монографии Работнова [2, 4, 10]. В списке литературы приведены некоторые статьи последних лет: по нелинейному деформированию и ползучести [27–30], по длительной прочности [31, 32], по применению и развитию моделей наследственной упругости [33–37]. В свое время Работнов интересовался задачами о концентрации напряжений [5, 13], и его модели получали развитие [12, 38]. В частности, он предложил модель “сдвигового анализа” для оценки концентрации напряжений около захватов при растяжении плоских образцов из композитов со слабой полимерной матрицей, и со временем эта модель была доведена до конкретных расчетных оценок [39, 40].

Неоднократно высказывал Юрий Николаевич соображения о сходстве композитов и биологических материалов типа древесины, но лишь с освоением мощных компьютеров удалось решить интересные задачи моделирования системы криволинейного армирования в зоне сучка [41–43], что позволяет с применением аддитивных технологий [44] проектировать принципиально новые, биоподобные (bio-inspired) виды соединений композитных деталей. Статьи [13, 15, 21] породили многие работы, посвященные механизмам разрушения композитов в виде образования трещин расслоения или расщепления [45, 46].

Вслед за Ю.Н. Работновым его ученики обратили свой интерес к биомеханике [47, 48], и в частности, к возможности создания ветвящихся упругих элементов, напоминающих крону дерева с постоянной суммарной площадью ветвей (правило Леонардо) [49, 50].

И наконец, маленькая заметка с идеей формулировки простого критерия прочности для модели ромба из нерастяжимых стержней [23] оказалась чрезвычайно плодотворной и нашла свое развитие в серии работ по рациональному армированию композитных сосудов давления [51–54].

Не только прямые ученики Юрия Николаевича (и их аспиранты) развивали его идеи. Многие известные представители Сибирской и Якутской научных школ были подготовлены в Новосибирском Академгородке и в аспирантуре Института машиноведения РАН и по-прежнему считают академика Ю.Н. Работнова одним из своих Учителей, что неоднократно отмечается в их публикациях [55].

Отметим, что только по базе Scopus за последние 5 лет (2018–2023 гг.) число ссылок на труды Ю.Н. Работнова в мире превышает 800, а за все прошедшие годы – сотни тысяч. Малая часть из них отмечена в приведенном кратком обзоре.

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

Вместо эпилога. Незаменимых людей нет, но бывают люди, потеря которых кажется невосполнимой.

×

About the authors

A. N. Polilov

Blagonravov Mechanical Engineering Research Institute

Author for correspondence.
Email: polilovan@mail.ru
Russian Federation, Moscow

References

  1. Lomakin E.V., Polilov A.N. Hereditary in science and in living (to 100-years jubilee of Yu.N. Rabotnov // Vestn. RAS, 2014, vol. 84, no. 4, pp. 60–64 (in Russian)
  2. Rabotnov Yu.N. Creep of Structural Elements. Moscow: Nauka, 2014. 752 p. (in Russian)
  3. Rabotnov Yu.N., Mileiko S.T. Low-Time Creep. Moscow: Nauka, 1970. 223 p. (in Russian)
  4. Rabotnov Yu.N. Elements of Hereditary Mechanics of Solids. Moscow: Nauka, 1977. 384 p. (in Russian)
  5. Rabotnov Yu.N. Effect of stress concentration on long-time strength // Mech. of Solids, 1967, no. 3, pp. 36–41. (in Russian)
  6. Rabotnov Yu.N., Polilov A.N. Problems of fracture of materials // Mehanika, 1974, no. 5, pp. 17–32. (in Russian)
  7. Fracture. Vol. 7. / Transl. ed. Rabotnov Yu.N. Moscow: Mir, 1976. 636 p. (in Russian)
  8. Composite Materials. Vol. 5. Fracture and Fatigue. Moscow: Mir, 1978. 486 p. (in Russian)
  9. Rabotnov Yu.N. Materials Resistance. Moscow: Fizmatlit, 1962. 455 p.
  10. Rabotnov Yu.N. Mechanics of Deformable Solids. Moscow: Nauka, 1988. 712 p. (in Russian)
  11. Gordon J. New Strength Theory. Moscow: Mir, 1971. 272 p. (in Russian)
  12. Polilov A.N. Crack arrest by week interface // Mech. of Solids, 1974, no. 1, pp. 68–72. (in Russian)
  13. Polilov A.N., Rabotnov Yu.N. Fracture of low shear strength composites near edge notches // Mech. of Solids, 1976, no. 6, pp. 112–119. (in Russian)
  14. Polilov A.N., Rabotnov Yu.N. The development of delaminations under compression of composites // Mech. of Solids, 1983, no. 4, pp. 166–171. (in Russian)
  15. Rabotnov Yu.N., Polilov A.N. On the fracture of composite pipes by shape of Chinese lantern // Mech. of Compos. Mater., 1983, no. 3, pp. 548–550. (in Russian)
  16. Rabotnov Yu.N. Elasto-plastic state of composite structure // in: Problems of Hydro-Dynamics and Mechanics of Continuous Media. Moscow: Nauka, 1969. pp. 411–415. (in Russian)
  17. Polilov A.N., Rabotnov Yu.N. Strength criterion for FRP // in: Proc. IV All-Union Conf. on Theor. &Appl. Mech. Kiev: 1976. p. 104. (in Russian)
  18. Rabotnov Yu.N., Polilov A.N. Strength criteria for fiber reinforced plastics // Advances in Research on the Strength and Fracture of Materials. 4th Int. Conf. Fract. (ICF 4), Waterloo, Canada, June 1977. Vol. 3B. Applications and Non-Metals. Sess. 3. Composites / Ed. by Taplin D.M.R. N.Y.: Pergamon, 1978. pp. 1059–1067.
  19. Rabotnov Yu.N., Polilov A.N. Strength criteria for reinforced plastics // in: Composite Materials. First USSR–Japanese Symp., 1977, Moscow, MSU, 1979, pp. 375–384.
  20. Rabotnov Yu.N., Polilov A.N., Kogaev V.P., et al. Fatigue strength of unidirectional CFRP under tension at the angle to direction of reinforcement // Mech. Compos. Mater., 1985, no. 1, pp. 44–47. (in Russian)
  21. Rabotnov Yu.N., Polilov A.N., Kogaev V.P. et al. Interlaminar strength criterion for CFRP under cyclic loading // Mech. Compos. Mater., 1982, no. 6, pp. 983–986. (in Russian)
  22. Rabotnov Yu.N., Danilova I.N., Polilov A.N. et al. Strength investigation of wounded epoxy CFRP and GFRP under torsion, tensile and 3-point bending load // Mech. Polym., 1978, no. 2, pp. 219–225. (in Russian)
  23. Rabotnov Yu.N. On strength of two-directions reinforced composites // Mech. Polym., 1978, no. 5, pp. 832–834 (in Russian)
  24. Rabotnov Yu.N. Introduction into Fracture Mechanics. Moscow: Nauka, 1987, 80 p. (in Russian)
  25. Rabotnov Yu.N. Problems of Deformable Solids Mechanics. Selected Papers. Moscow: Nauka, 1991. 196 p. (in Russian)
  26. Polilov A.N. Etudes on Mechanics of Composites. Moscow: Nauka, 2015. 320 p. (in Russian)
  27. Dumanski A. M., Alimov M.A., Terekhin A.V. Experiment- and computation-based identification of mechanical properties of fiber reinforced polymer composites // J. Phys. Conf. Ser., 2019, iss. 1, art. no. 1158.
  28. Ashpiz E.S., Suvorova Yu.V., Alexeeva S.I. et al. Modelling of creep process of geo-grids at long-time loading // Zavod. Lab. Diagn. Mater., 2006, no. 1, pp. 49–54. (in Russian)
  29. Khokhlov A.V. The analysis of own properties of creep curves under step loading, being born from non-linear Rabotnov’s equations for visco-elastic materials // Vestn. Bauman MSTU. Natural Sci., 2017, no. 3, pp. 93–123. (in Russian)
  30. Maslov B.P. Nonlinear hereditary creep of isotropic composites of random structure // Int. J. Appl. Mech., 2022, vol. 58, no. 1, pp. 75–90.
  31. Lokoshenko A.V., Fomin L.V., Larin N.S. Creep rupture strength of rod stretched in an aggressive medium with various two-connected forms of their cross-section // Mech. of Solids, 2021, vol. 56, no. 7, pp. 1277–1294.
  32. Radchenko V.P., Afanas’eva E.A., Saushkin M.N. Prediction high temperature and long-term strength of visco-plastic material using a leader sample // Vestn. Samar. Gos. Univ., Ser. Fiz.-Mat. Nauki, 2023, vol. 27, no. 2, pp. 292–308. (in Russian)
  33. Suvorova Yu.V. On Rabotnov’s non-linear hereditary-type equation and its applications // Mech. of Solids, 2004, no. 1, pp. 174–181. (in Russian)
  34. Khokhlov A.V. The analysis of properties of relaxation curves with initial stage of ramp-deformation, being born from Rabotnov’s non-linear hereditary theory // Mech. Compos. Mater., 2018, no. 4, pp. 687–708.
  35. Shitikova M.V. Impact of rigid sphere on an infinite viscoelastic Kirchhoff-Love plate considering volume and shear relaxations // Vestn. St. Pb. Gos. Univ. Ser. Math., Mech., Astron., 2023, vol. 10 (68), no. 1, pp. 139–154. (in Russian)
  36. Vlasov D.D., Polilov A.N., Tatus N.A. On the problems of creep curves prediction for polymeric composites using vibration tests // Mashinostr. i Inzh. Obraz., 2021, no. 3–4, pp. 37–48. (in Russian)
  37. Vlasov D.D., Polilov A.N. The possibility of creep prediction of viscoelastic polymer composites using frequency dependences of complex modulus components // Mech. Compos. Mater., 2022, vol. 58, no. 1, pp. 274–289.
  38. Tatus N.A., Polilov A.N., Vlasov D.D. The effect of holes on the strength reduction of composite specimens with different structure of fiber reinforcement // Zavod. Lab. Diagn. Mater., 2022, vol. 88, no. 4, pp. 58–65. (in Russian)
  39. Polilov A.N., Arutunova A.S., Tatus’ N.A. Effect of stress concentration near grips on composite strength realization under tension // Zavod. Lab. Diagn. Mater., 2020, vol. 86, no. 11, pp. 48–59.
  40. Polilov A.N., Tatus’ N.A. Estimating stress concentration near grips under tension of polymer composite samples // J. Machin. Manuf.&Reliab., 2020, vol. 49, no. 5, pp. 446–456.
  41. Malakhov A.V., Polilov A.N. Design of composite structures reinforced curvilinear fibres // Compos. Pt. A. Appl. Sci.&Manuf., 2016, vol. 87, pp. 23–31.
  42. Malakhov A.V., Polilov A.N. Design algorithm of rational fiber trajectories in arbitrarily loaded composite plate // J. Machin. Manuf.&Reliab., 2017, vol. 46, no. 5, pp. 479–487.
  43. Malakhov A.V., Polilov A.N., Li D, Tian X. Increasing the bearing capacity of composite plates in the zone of bolted joints by using curvilinear trajectories and a variable fiber volume fraction // Mech. Compos. Mater., 2021, vol. 57, no. 3, pp. 287–300.
  44. Polilov A.N., Tatus’ N.A. Technological mechanics of composites, based on learning of Nature experience in design of bio-materials&structures // Probl. Mashinostr. i Avtomat., 2021, no. 1, pp. 59–85. (in Russian)
  45. Plitov I.S., Polilov A.N. Rational dimensions of segments of bamboo stems and composite tubes subjected to compression, flexure, and torsion // J. Machin. Manuf.&Reliab., 2015, vol. 44, no. 3, pp. 239–248.
  46. Polilov A.N., Vlasov D.D., Tatus’ N.A. Precise criterion of delamination of composite beam under bending // Zavod. Lab. Diagn. Mater., 2023, vol. 89, no. 10, pp. 63–73.
  47. Polilov A.N., Tatus’ N.A. Biomechanical principles of optimum design of composite constructions // Mashinostr. i Inzh. Obraz., 2018, no. 3 (56), pp. 19–26. (in Russian)
  48. Polilov A.N., Tatus’ N.A. Biomechanics of Strength of Fibrous Composites. Moscow: Fizmatlit, 2018. 328 p. (in Russian)
  49. Polilov A.N., Tatus N.А. Designing branching or shaped composite elements by analogy with the structure of treetops // J. Machin. Manuf.&Reliab., 2017, vol. 46, no. 4, pp. 385–393.
  50. Polilov A.N., Tatus N.А., Tian X., Arutjunova A.S. Equistrong branchy composite beams with a constant total area of variable elliptic cross sections // Mech. Compos. Mater., 2019, vol. 55, no. 3, pp. 465–482.
  51. Polilov A.N., Sklemina O.Yu., Tatus’ N.А. Ply-by-ply calculation method and failure criteria of composite pressure vessels wounded with symmetric pair of layers. Pt. 1. Features of ply-by-ply design method for composite structure made of symmetric pair of layers // Mashinostr. i Inzh. Obraz., 2020, no. 3 (64), pp. 21–30. Pt. 2. Simple strength criteria for angle-ply reinforced tubes under two-axis tension // Mashinostr. i Inzh. Obraz., 2020, no. 4 (65), pp. 14–25. (in Russian)
  52. Polilov A.N., Vlasov D.D., Sklemina O.Yu., Tatus’ N.A. Strength criteria of obliquely wound composite tubes under biaxial tension // Strength of Mater., 2021, vol. 53, no. 5, pp. 765–774. https://doi.org/10.1007/s11223-021-00342-7
  53. Polilov A.N., Sklemina O.Yu., Tatus’ N.A. Design method of reinforcement structure with symmetric pairs of layers by the example of composite gas tank // Mech. Compos. Mater., 2022, vol. 57, no. 6, pp. 769–784. https://doi.org/10.1007/s11029-022-09998-x
  54. Polilov A.N., Vlasov D.D., Sklemina O.Yu. Tatus’ N.A. Evaluation of the strength of a composite cylinder for compressed gas // J. Machin. Manuf.&Reliab., 2022, vol. 51, no. 1, pp. 46–54. https://doi.org/10.3103/S1052618822010083
  55. Load Capacity and Safety of Metal-Composite Gas Tanks for Space Apparatus / Ed. by Moskvichev V.V., Testoedov N.A. Novosibirsk. Nauka, 2021. 440 p. (in Russian)

Supplementary files

Supplementary Files
Action
1. JATS XML
2. Fig. 1. Limit surfaces in the form of: a - rounded cylinder for unidirectional composite, b - “sofa cushion” (ottoman cushion) for orthogonally reinforced plastic; c - three-dimensional oblique parallelepiped for reinforcement in three and more directions.

Download (102KB)
3. Fig. 2. Loading scheme (a) and strength dependences of orthogonally reinforced composites on the direction of specimen cutting: b - unidirectional (0); c - “equal-strength” (0/90); d - (04/901). Line 1 - (2.3); line 2 - (2.4).

Download (331KB)
4. Fig. 3. Limit curves for CFRP (0/90) in coordinates .

Download (107KB)
5. Fig. 4. Deformation scheme of the obliquely reinforced composite element.

Download (112KB)
6. Fig. 5. Limit surfaces for biaxial loading of FRP pipes: with winding: ±55° (a), ±72.5° (b), ±35° (c).

Download (281KB)
7. Fig. 6. Schematic of kink formation under compression.

Download (157KB)
8. Fig. 7. Schematic of internal loss of stability of fibers (layers) under compression.

Download (188KB)

Copyright (c) 2024 Russian Academy of Sciences

Согласие на обработку персональных данных с помощью сервиса «Яндекс.Метрика»

1. Я (далее – «Пользователь» или «Субъект персональных данных»), осуществляя использование сайта https://journals.rcsi.science/ (далее – «Сайт»), подтверждая свою полную дееспособность даю согласие на обработку персональных данных с использованием средств автоматизации Оператору - федеральному государственному бюджетному учреждению «Российский центр научной информации» (РЦНИ), далее – «Оператор», расположенному по адресу: 119991, г. Москва, Ленинский просп., д.32А, со следующими условиями.

2. Категории обрабатываемых данных: файлы «cookies» (куки-файлы). Файлы «cookie» – это небольшой текстовый файл, который веб-сервер может хранить в браузере Пользователя. Данные файлы веб-сервер загружает на устройство Пользователя при посещении им Сайта. При каждом следующем посещении Пользователем Сайта «cookie» файлы отправляются на Сайт Оператора. Данные файлы позволяют Сайту распознавать устройство Пользователя. Содержимое такого файла может как относиться, так и не относиться к персональным данным, в зависимости от того, содержит ли такой файл персональные данные или содержит обезличенные технические данные.

3. Цель обработки персональных данных: анализ пользовательской активности с помощью сервиса «Яндекс.Метрика».

4. Категории субъектов персональных данных: все Пользователи Сайта, которые дали согласие на обработку файлов «cookie».

5. Способы обработки: сбор, запись, систематизация, накопление, хранение, уточнение (обновление, изменение), извлечение, использование, передача (доступ, предоставление), блокирование, удаление, уничтожение персональных данных.

6. Срок обработки и хранения: до получения от Субъекта персональных данных требования о прекращении обработки/отзыва согласия.

7. Способ отзыва: заявление об отзыве в письменном виде путём его направления на адрес электронной почты Оператора: info@rcsi.science или путем письменного обращения по юридическому адресу: 119991, г. Москва, Ленинский просп., д.32А

8. Субъект персональных данных вправе запретить своему оборудованию прием этих данных или ограничить прием этих данных. При отказе от получения таких данных или при ограничении приема данных некоторые функции Сайта могут работать некорректно. Субъект персональных данных обязуется сам настроить свое оборудование таким способом, чтобы оно обеспечивало адекватный его желаниям режим работы и уровень защиты данных файлов «cookie», Оператор не предоставляет технологических и правовых консультаций на темы подобного характера.

9. Порядок уничтожения персональных данных при достижении цели их обработки или при наступлении иных законных оснований определяется Оператором в соответствии с законодательством Российской Федерации.

10. Я согласен/согласна квалифицировать в качестве своей простой электронной подписи под настоящим Согласием и под Политикой обработки персональных данных выполнение мною следующего действия на сайте: https://journals.rcsi.science/ нажатие мною на интерфейсе с текстом: «Сайт использует сервис «Яндекс.Метрика» (который использует файлы «cookie») на элемент с текстом «Принять и продолжить».