Контактная задача о кручении круглым штампом трансверсально-изотропного упругого полупространства с неоднородным трансверсально‑изотропным покрытием Текст научной статьи по специальности «Физика»

ФИЗИКО-МАТЕМАТИЧЕСКИЕ НАУКИ

УДК 539.3

Контактная задача о кручении круглым штампом трансверсально-изотропного упругого полупространства с неоднородным трансверсально-изотропным покрытием1

А. С. Васильев, Е. В. Садырин

(Донской государственный технический университет),

И. А. Федотов

(Технологический университет Тшване, Претория, ЮАР)

Рассмотрена осесимметричная контактная задача теории упругости о кручении круговым штампом трансвер-сально-изотропного полупространства с неоднородным покрытием. С помощью техники интегральных преобразований решение задачи сведено к решению интегрального уравнения. Построено явное выражение трансформанты ядра интегрального уравнения для однородного трансверсально-изотропного покрытия. Предложена схема численного построения трансформанты ядра в случае произвольного закона изменения модуля сдвига по глубине покрытия. Используя аппроксимацию трансформанты ядра некоторыми специальными функциями, построено двухсторонне асимптотически точное приближённое аналитическое решение задачи. Проанализирован частный случай изменения модулей сдвига в трансверсально-изотропном покрытии, позволяющий построить явное замкнутое решение задачи. Численно исследуется модельная задача о кручении биоматериала. Считается, что покрытие представляет собой хрящевую ткань, а свойства подложки соответствуют кортикальной кости.

Ключевые слова: контактные задачи, теория упругости, кручение, неоднородность, анизотропия.

Введение. Первые работы по моделированию задачи о кручении упругого полупространства круглым штампом относятся к середине прошлого века. Задача о кручении однородного упругого полупространства круглым штампом впервые была сформулирована и решена в динамической постановке Рейснером и Сагочи [1]. В работе Снеддона [2] данная задача, с использованием техники интегральных преобразований, сведена к решению интегрального уравнения.

Грилицкий В. Д. [3] построил решение задачи о кручении круглым штампом изотропной двухслойной среды и ортотропного упругого слоя в виде ряда по степеням б 1, где б — отношение толщины первого слоя к радиусу штампа. Kassir М. К. [4] построил решение задачи о кручении неоднородного изотропного полупространства в случае степенного изменения модуля сдвига по глубине. Решение задачи для изотропного полупространства с неоднородным покрытием при произвольным законе изменения модуля сдвига по глубине получено Айзиковичем С. М. [5]. Erguven М. Е. рассмотрел задачу о кручении трансверсально-изотропного упругого полупространства с однородным покрытием [6].

Задача активно изучается по настоящее время, в частности, Wang Y.-S. использовал кусочно-линейную аппроксимацию трансформанты ядра для сведения задачи к решению интегрального уравнения [7]. В работе [8] решена задача о кручении трансверсально-изотропного пьезоэлектрического однородного полупространства.

1 Работа выполнена при поддержке Минобрнауки РФ (ГК № 11.519.11.3028, соглашения №14.В37.21.1632, 14.В37.21.1131, 14.132.21.1693) и РФФИ (12-07-00639_а).

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

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

Постановка задачи. Недеформируемый круглый штамп с плоским основанием жёстко сцеплен с верхней гранью Г упругого неоднородного трансверсально-изотропного полупространства £2. С полупространством связана цилиндрическая система координат г,ср,г. Штамп контактирует с полупространством по поверхности г = О, г < а .

К штампу приложен крутящий момент М, ось которого совпадает с осью г (рис. 1). Под действием этого момента штамп повернётся относительно оси г на угол е, вызвав деформацию кручения О.. Модули сдвига С„р, ^полупространства с глубиной изменяются по законам:

С(г) -Н <г <0 ^ _[ -Н<г< О

|32

(1)

где /др, /ф/

некоторые функции, определяющие закон изменения модулей сдвига в покрытии, причём — непрерывная, а /Ц>7 — непрерывно дифференцируемая; параметры З1, Зг задают скачок упругих свойств на границе покрытие-подложка.

Вне штампа поверхность Г не нагружена. При сделанных предположениях граничные условия задачи имеют вид:

г >а г <а.'

Учитывая, что слой подвержен лишь деформации скручивания, ненулевым является лишь смещение вдоль оси ср:

и2=иг= 0, £/ф =и(г,г).

Напряжения, в соответствии с законом Гука, примут вид:

г = 0: о =т =0,

К=о,

К=гг'

(2)

а=а=а=т =0,т = Є (г)

г ф г к 1 фг фг V /

5£Уф т

Зг '

дг

Ф

Г

(3)

Введём обозначения:

Тфг

ту1' ~Н<г<0

т(2), г<-Н,

!Г\7 > '

"ф =

и^,-Н<г<О

и

(2)

г<-Н

При г ^ оо и г ^ -оо напряжения исчезают. Считаем, что перемещения и напряжения сопрягаются на границе изменения закона неоднородности:

2 = -Н: т«=т«, и®=и®. (4)

Требуется определить закон распределения контактных касательных напряжений под штампом

Т^1=0=Та(Г)' Г^а-

Сведение задачи к решению интегрального уравнения. Уравнения равновесия при данной постановке принимают вид:

^ £Й? ди(1)

+ ,ф2 ф + С

Г дг

М-")

«ч1 1>2 Ф

дг2

\ ди® , *чг'

Г дГ V дг2

д2и^

4- = 0, - Н <г <0,

дг*

Ы2)

Л Ф

д2и®

+ ^фг(-/У)_|- = 0, г <-/7.

(5)

(6)

Для сведения задачи к решению интегрального уравнения используем интегральное преобразование Ханкеля:

Обозначим

СО

и^{г,г) 7=1,2.

(7)

(8)

Подставив (7), (8) в (5) и (6), и приравнивая к нулю в полученных соотношениях подынтегральные выражения, получим систему обыкновенных дифференциальных уравнений относительно функций 1/;(у,*),7=1,2:

^Г(У^) + у(^)^1'(У^)-У2к(г)Ц(у,г) = 0, -Н<г< 0, (9)

^'(у,г)-у2к(-Я)^2(у,г) = 0, -оо <г<-Н. (10)

Общее решение уравнения (10) имеет вид:

и2{ч,г) = В1{ч)е

У^к(-Н)

А(у)е

Из условия ограниченности перемещений при г ^ -оо имеем В2 (у) = 0, тогда

иг(у,*) = в,(

Условия сопряжения (4) можно записать в виде:

^ {у,-Н) = и2 (у,-Н), и[ (у,-Н) = р2^' (у,

Или, используя явный вид функции и2 (у,г) (11), получим:

(П)

Из (3) при г = 0 и (7) следует:

и^'~ни2,£Щ.

и,(ч,-Н) Р2"1 ' ’

и;(у,0)‘Т(у)/в1К(0).

(12)

(13)

Здесь Ду) связана с та (г) преобразованием Ханкеля:

СО

г(у) = / Тэ (р)Л (ру)р^р ■ (14)

О

Введём вспомогательную функцию и/ (у,г):

Щ (у,г) = у и, (у ,1)Т1 (у)^ф (0)^(0). (15)

Ввиду однородности, уравнение (9) относительно функции ^ (у,г) не изменится. Условия (12),

(13) относительно примут вид:

(16)

^Г(У,°) = У-До) ■ (17)

Далее ^ (у,0) будем обозначать /-*(у). Затем, используя условие (2), можно записать интегральное уравнение нашей задачи:

".и.ф.(4.0)3,(/■*)*</* = , * 1 г(¥)Г(У)ЛИ)<#У = ГЕ, г<а.

0 V

Используя (14), уравнение можно представить в виде:

а со

|та (p)pjr (у) Л (гу)Л (PY)^P = re^G,!p (0)G!pz (0),

г <а.

Сделаем замену переменных:

Wy = u, К = Н/а, г' = r/a, р' = p/a,Z.(t/) = Г (и//У),т(р') = тв (р'з).

Далее штрихи опускаем, тогда интегральное уравнение задачи примет вид:

1 со

j т (р) pj L (u) {игК1) Jj (урЛ 1) dudp = kre^G^ (0) Gw (0), r< 1. (18)

о о

Таким образом, для построения интегрального уравнения задачи необходимо построить функцию UI (у,0) из краевой двухточечной задачи (9), (16), (17).

Построение трансформанты ядра интегрального уравнения. Для некоторых частных случаев изменения модуля сдвига в покрытии удаётся получить явное выражение для трансформанты ядра L(u). Например, для покрытия, в котором модули сдвига постоянны:

г К- -Н<2<й (<?£„ -H<z<0

G"( -»<z<-Н, »'( )_|g*,

В этом случае решение уравнения (9) можно построить в явном виде, и трансформанта ядра принимает вид:

1 Gc б -Gs б \GC

L(u\ = 1+|je_________ = °фг°1 5 = 5 =

{ ) 1_[je-2uHbl ' V GCb +GSX ' 1 JGC ' 2

г*-' 1 (рг 2 \ Ф ^

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

Для произвольного закона изменения модуля сдвига по глубине значение трансформанты ядра может быть получено численно. Введём обозначения:

со, =и; (у,г), со2 =6/;* (у,г).

Уравнение (9) перепишем в виде линейной системы обыкновенных дифференциальных уравнений первого порядка:

/ ч1 2 2 / \ -Н <z< О , С02 = -v(z)со2 + у k(z)coi;

или в матричной форме:

«Г 0 О

= Аш, -Н <z < 0, А = dz

Y2k(z) -v(z)_

со =

со,

V 2 /

Краевые условия (16), (17) принимают вид:

ЩЩ11z-H = 3iVyp-H)' “21z=o = Y^M ■ (19)

Решение краевой двухточечной задачи при фиксированном у построим, используя метод моделирующих функций [9], предложенный для построения трансформант ядер динамических контактных задач. Ищем 63(y,z) в виде:

c6(y,z) = 61(y)a(y,z)evz.

Вектор a(y,z) определяется из задачи Коши при фиксированном у:

da . ■ „

— = Аа-уа, -H<z< О, dz

| z=-h = 1/Зг' ^21 z=-h = Y^/K (_^) ■

Затем, используя второе условие (19), находим ^ (у) и тем самым сох (у,0):

6i(y) = Y321 (у,0)Л/к(0),

Z.(y/W) = Г (у) = У; (у,0) = щ (у,0) = у^/к^О)^ (у,0)а2-1 (у,0).

Аналогично работе [10] можно показать, что трансформанта ядра L(u) обладает свойствами:

m=JA-(0!40) ' |im/-(£y)=i- (20)

Построение решения задачи. Будем говорить, что функция f[u) принадлежит классу V\N, где N— натуральное число, если она представима в виде:

f И = П^2 +%' а,<в, *С,А,* AJ,Bi ф В} при i *j.

/=1 U + Dj

Для трансформанты ядра Ln(u)<e Пд, можно построить [11] решение интегрального уравнения (18) в форме:

Т(Г)Л£^ЛЧМЧ(^^ + |С,7(Г,4Л)) ■ (21)

где постоянные С; определяются из системы линейных алгебраических уравнений:

У С D f— —1 + ^ + - 0 к - 1 2 • • • N

Здесь использованы обозначения:

лЧ sh Ar rs\\A .) ch Atdt

^М)=—+ ,,

y/l-r2 |l + Vl - Г2 ) Г л/f2 - Г2 +Vf2 - Г2 j p(A,B) = (AcUA + BsUA)/(B2 -A2).

29

Таким образом, для построения решения исходной контактной задачи достаточно аппроксимировать трансформанту ядра интегрального уравнения функцией класса П№ Для построения аппроксимации можно использовать алгоритм, описанный в работе [11], где показано, что решение (21) является асимптотически точным при Л—>0 и Л—>оо.

Решение задачи для частного случая изменения модулей сдвига. Пусть для модулей сдвига выполнены равенства:

G'V(z)G'fz(z) = c2= const, р1 = р2=1. (22)

Тогда первое равенство (8) можно переписать в виде:

^Щzj = GrЛz)/c ■ (23)

Учитывая (23), уравнение (9) принимает вид:

-|f77T1/;(v'z) + ^T^_^G''W1/l<Y'z) = 0' <24>

GrAz) GrAz) с

Общее решение уравнения (24) имеет вид:

U;(y,z) = C1(y)eF^/c +С2(у)е-™, (25)

где F(z) = jGrv(z)dz.

Продифференцируем U*(yrz) ПО Z.

^и; (v,z) = lG„(z)[ct (¥)е™‘ -с, (¥)<г™‘].

Тогда из краевого условия (16), учитывая (23), следует:

С2 (у) = 0. (26)

Из краевого условия (17) и равенств (23), (26) получаем следующее:

C1(y)evF(0)=l. (27)

Из (25), (26), (27) следует, что для любого у выполнено тождество

£у;(у,0) = 1.

Тогда

L(u).1.

Таким образом, если для модулей сдвига трансверсально-изотропного полупространства с неоднородным трансверсально-изотропным покрытием выполнены равенства (22), то трансформанта ядра интегрального уравнения тождественно равна единице. Следовательно, решение задачи совпадает с решением для однородного полупространства с модулем сдвига равным с и имеет вид:

/ \ 4 г

т(г) = -ЕС--------

Численные примеры. Рассмотрим трансверсально-изотропную упругую подложку с модулями сдвига(-/У) = 5, G (~Н) = 6 и пусть модули сдвига на поверхности покрытия равны

Gr4> (0) = (0) = 0,608. Такие значения модулей сдвига соответствуют свойствам хряща (покры-

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

1. frv (z) = 5,03 - 4,422e5z/w , fw (z) = 6,037 - 5,428e5z/w ;

2. frv (z) = 0,608 - 4,392z/W, ^ (z) = 0,608 -5,392z/H ;

3. ^(z) = ^(z) = 0,608;

4. frv{z) = f4Z{z) = 0,118 + 0,489e5z/w.

Эти законы описывают трансверсально-изотропные покрытия, модули сдвига которых возрастают экспоненциально (закон 1) и линейно (закон 2); однородное изотропное покрытие (закон 3) и изотропное покрытие (закон 4), упругие свойства которого экспоненциально убывают, а на границе покрытие-подложка наблюдается скачок модулей сдвига С„р, 6ф7 в 41 и 49 раз соответственно (31=41, 32=49). Графики изменения модулей сдвига по глубине для законов 1—4 изображены на рисунке 2.

0.1 0.6 (7,. 5 6 0.1 0.6 (7,. 5 6

<Р<,

Рис. 2. Изменение модулей сдвига по глубине в материалах 1—4

Трансформанты ядер, соответствующие покрытиям 1—4, изображены на рисунке 3. Построены аппроксимации трансформант ядер функциями класса П№ их погрешности равны 0,8 %,

0,82 %, 0,64 % и 0,49 % соответственно. Под погрешностью аппроксимации понимаем следующее выражение:

Д, = тах ^ (£/)//- (£/) -1| х 100 %.

10 10 10 10 10 10 и 10

Рис. 3. Трансформанты ядер для материалов 1—4

Введём величину:

тге, (Л,г) = т(Л,г)/т„от(г), где Л = Н/а — геометрический параметр задачи, характеризующий относительную толщину покрытия; т(Л,л) — контактные напряжения под штампом для рассматриваемой неоднородной среды; тИот (л) — контактные напряжения однородной изотропной среды с модулем сдвига, равным (-Н)С^(-Н).

Величина тге| характеризует распределение контактных касательных напряжений под штампом для неоднородной среды по сравнению с однородной средой.

О 0.25 0.5 0.75 г 0.99

в) Л=0.014

Рис. 4. Контактные напряжения для законов 1—4

Графики величины тге| для материалов 1—4 в случае относительно толстого покрытия

(Л=100), тонкого покрытия (Л« 1/100) и покрытия, сравнимого по размерам с размерами штампа (Л=1), изображены на рисунке 4. Из графиков видно, что неоднородность покрытия оказывает существенное влияние на решение контактной задачи, особенно для покрытий малой толщины. Наиболее существенное перераспределение контактных напряжений по сравнению со случаем однородного полупространства наблюдается на краю штампа (/"—>1).

0.13-т

0.12-

0 0.25 0.5 0.75 г 0.99

а) Л=1

0.105

0

0.25 0.5

б)Л=100

Заключение. Построено приближённое аналитическое решение контактной задачи о кручении трансверсально-изотропного упругого полупространства с неоднородным трансверсально-изотропным покрытием, эффективное для любых относительных толщин покрытий (значений параметра Л). Показано, что в случае, когда модули сдвига трансверсально-изотропного полупространства изменяются в соответствии с условием (22), распределение контактных напряжений под штампом совпадает со случаем однородного изотропного полупространства.

На основании высокой точности аппроксимации трансформант ядер интегральных уравнений (менее 0,8 %) можно утверждать о малой погрешности построенных решений контактной задачи даже для покрытий сложной структуры (материал 4).

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

Библиографический список

1. Reissner, Е. Forced torsional oscillations of an elastic half-space / E. Reissner, H. F. Sagoci // Journal of Applied Physics. — 1944. — Vol. 15. — № 9. — Pp. 652—654.

2. Sneddon, I. N. The Reissner—Sagoci problem / I. N. Sneddon // Proceedings of the Glasgow Mathematical Association. — 1966. — Vol. 7. — № 3. — Pp. 136—144.

3. Грилицкий, Д. В. Кручение двухслойной упругой среды / Д. В. Грилицкий // Прикладная механика. — 1961. — Т. 7. — № 1. — С. 89—94.

4. Kassir, М. К. The Reissner—Sagoci problem for a non-homogeneous solid / М. K. Kassir // International Journal of Engineering Science. — 1970. — Vol. 8. — № 10. — Pp. 875—885.

5. Айзикович, С. М. Кручение круглым штампом неоднородного полупространства / С. М. Айзикович // Расчёт оболочек и пластин : сб. трудов. — Ростов-на-Дону : РИСИ, 1978. — С. 156-169.

6. Erguven, М. Е. Torsion of a transversely isotropic elastic layer bonded to a transversely isotropic half-space / М. E. Erguven // International Journal of Engineering Science. — 1986. — Vol. 24. — № 9. - Pp. 1501-1509.

7. Liu, Tie-Jun. Reissner—Sagoci problem for functionally graded materials with arbitrary spatial variation of material properties / Tie-Jun Liu, Yue-Sheng Wang // Mechanics Research Communications. — 2009. — Vol. 36. — № 3. — Pp. 322—329.

8. Xiong, Su-ming. The Reissner—Sagoci problem for transversely isotropic piezoelectric halfspace / Su-ming Xiong, Guang-zheng Ni, Peng-fei Hou // Journal of Zhejiang University SCIENCE A. — 2005. - Vol. 6. - № 9. - Pp. 986-989.

9. Бабешко, В. А. Методы построения матриц Грина для стратифицированного упругого полупространства / В. А. Бабешко, Е. В. Глушков, Н. В. Глушкова // Журнал вычислительной математики и математической физики. — 1987. — Т. 27. — № 1. — С. 93—101.

10. Айзикович, С. М. О свойствах функций податливости, соответствующих слоистому и непрерывно-неоднородному полупространству / С. М. Айзикович, В. М. Александров // Доклады академии наук СССР. — 1982. — Т. 266. — № 1. — С. 40—43.

11. Айзикович, С. М. Асимптотические решения контактных задач теории упругости для неоднородных по глубине сред / С. М. Айзикович // Прикладная математика и механика. — 1982. — Т. 46. — № 1. — С. 148—158.

Материал поступил в редакцию 08.10.2012.

References

1. Reissner, Е., Sagoci, Н. F. Forced torsional oscillations of an elastic half-space. Journal of Applied Physics, 1944, vol. 15, no. 9, pp. 652—654.

2. Sneddon, I. N. The Reissner—Sagoci problem. Proceedings of the Glasgow Mathematical Association, 1966, vol. 7, no. 3, pp. 136—144.

3. Grilitskiy, D. V. Kruchenie dvuxslojnoj uprugoj sredy'. [Double-layer elastic medium torsion.]

Prikladnaya mexanika, 1961, vol. 7, no. 1, pp. 89—94 (in Russian).

4. Kassir, М. K. The Reissner—Sagoci problem for a non-homogeneous solid. International Jour-

nal of Engineering Science, 1970, vol. 8, no. 10, pp. 875—885.

5. Aizikovich, S. M. Kruchenie krugly'm shtampom neodnorodnogo poluprostranstva. [Round-die torsion of inhomogenous half-space.] Raschet obolochek i plastin : sb. trudov.[Analysis of shells and plates: coll.of sci.papers.] Rostov-on-Don : Rostov State University of Civil Engineering, 1978, pp. 156— 169 (in Russian).

6. Erguven, M. E. Torsion of a transversely isotropic elastic layer bonded to a transversely isotropic half-space. International Journal of Engineering Science, 1986, vol. 24, no. 9, pp. 1501—1509.

7. Liu, Tie-Jun, Wang, Yue-Sheng. Reissner—Sagoci problem for functionally graded materials with arbitrary spatial variation of material properties. Mechanics Research Communications, 2009, vol. 36, no. 3, pp. 322—329.

8. Xiong, Su-ming, Ni, Guang-zheng, Hou, Peng-fei The Reissner—Sagoci problem for transversely isotropic piezoelectric half-space. Journal of Zhejiang University SCIENCE A., 2005, vol. 6, no. 9, pp. 986—989.

9. Babeshko, V. A., Glushkov, E. V., Glushkova, N. V. Metody' postroeniya matricz Grina dlya stratificirovannogo uprugogo poluprostranstva. [Green matrices construction methods for stratified elastic half-space.] Zhurnal vy'chislitel'noj matematiki i matematicheskoj fiziki, 1987, vol. 27, no. 1, pp. 93—101 (in Russian).

10. Aizikovich, S. M., Alexandrov, V. M. O svojstvax funkcij podatlivosti, sootvetstvuyushhix sloistomu i neprery'vno-neodnorodnomu poluprostranstvu. [On receptance function properties corresponding to layered and continuous-inhomogeneous half-space.] Doklady' akademii nauk SSSR, 1982, vol. 266, no. 1, pp. 40—43 (in Russian).

11. Aizikovich, S. M. Asimptoticheskie resheniya kontaktny'x zadach teorii uprugosti dlya ne-odnorodny'x po glubine sred. [Asymptotic solutions to contact problems of elasticity theory for inho-mogeneous in depth media.] Prikladnaya matematika i mexanika, 1982, vol. 46, no. 1, pp. 148—158 (in Russian).

CONTACT PROBLEM ON TORSION OF TRANSVERSELY ISOTROPIC ELASTIC HALF-SPACE WITH INHOMOGENEOUS TRANSVERSELY ISOTROPIC COATING BY ROUND DIE1 A. S. Vasilyev, E. V. Sadyrin

(Don State Technical University),

I. A. Fedotov

(Tshwane University of Technology, Pretoria, South Africa)

An axisymmetric contact problem of the elasticity theory on the torsion of the transversely isotropic half-space with the inhomogeneous coating by the round die is considered. The problem solution is reduced to the integral equation solution using the integral transformations technique. The explicit expression of the kernel transform of the integral equation for the homogeneous transversely isotropic coating is constructed. The numerical construction scheme of the kernel transform in the case of the arbitrary law of variation of the coating depth shear modulus is suggested. The bilateral asymptotically exact approximate analytical problem solution is constructed using the kernel transform approximation by some special functions. A special case of variation of the shear moduli in the transversely isotropic coating that allows constructing an explicit dosed problem solution is analyzed. The model problem on the biomaterial torsion is studied numerically. The coating is considered a cartilaginous tissue, and the substrate properties correspond to a cortical bone.

Keywords: contact problems, theory of elasticity, torsion, inhomogeneity, anisotropy.

1 The research is done with the support from the Ministry of Education and Science of the Russian Federation (State Contracts No. 11.519.11.3028, Agreements Nos. 14.V37.21.1632, 14.V37.21.1131, 14.132.21.1693) and the RFFI (No. 12-07-00639_a).