Решение плоской задачи упругой среды, армированной тремя семействами волокон Текст научной статьи по специальности «Математика»

Вычислительные технологии

Том 10, Специальный выпуск, 2005

РЕШЕНИЕ ПЛОСКОЙ ЗАДАЧИ УПРУГОЙ СРЕДЫ, АРМИРОВАННОЙ ТРЕМЯ

___ и ___

СЕМЕЙСТВАМИ ВОЛОКОН

Н. А. Федорова Красноярский государственный технический университет, Россия

e-mail: ran@krsk.info

Systems of equations governing the plane problem for the elastic media reinforced with possible combinations of three nonextensible and equally stretchable fibres are derived. Some exact solutions of model [1, 2] are constructed using the invariant solutions of the derived partial differential equations.

Рассматриваются плоские задачи упругости для сред, армированных тремя семействами волокон. Пусть армирование выполнено волокнами постоянного поперечного сечения. Для описания композита используется структурная модель [3]. Модель содержит алгебраические и дифференциальные уравнения относительно интенсивностей армирования Uk(x,y), компонент тензора деформаций eij(x,y), деформаций в волокнах первого, второго и третьего семейств ek(x,y), напряжений в волокнах первого, второго и третьего семейств — Ok(x, y), осредненных напряжений оц(x,y), где x,y — декартовы координаты, ^k(x, y) — углы армирования, индексы i,j = 1, 2, k = 1, 2, 3. В рамках принятых обозначений при условии постоянства поля температур T исходная система имеет вид

(Uk 1ki),i + (Uk 1k2),2 = 0; (1)

e11l2l + e22lk2 + 2e12lk1lk2 = ek; (2) e11,22 + e22,11 = 2е12 ,12- (3)

Использованы обозначения: eT = Q^T, e°k = ek + eT, lk1 = cos , lk2 = sin , aa — коэффициенты линейного расширения материала k-го семейства волокон (k = 1, 2, 3). Символы ,1, ,2 означают частное дифференцирование по координатам x, y соответственно. Правая часть в (2) учитывает как случай равнодеформированных (ek = const, ek = const + e^), так и случай нерастяжимых (ek = 0, e0 = еT) семейств волокон и их возможные комбинации (некоторые из семейств волокон равнодеформируемы, другие — нерастяжимы).

Осредненные напряжения оц (x, y) имеют вид

3

Оц = ОоЦ + OkUkIkilkj, 0 = 1 - (U1(x, y) + U2(x, y) + U3 (x, y)). (4)

k=1

© Институт вычислительных технологий Сибирского отделения Российской академии наук, 2005.

В (4) напряжения в связующем определяются по формула

м

E E

aÍi = (i _ v2) (£ii + v£jj _ «c(l + v)T), aij = (1 + v) £ij, j =3 - г, г = 1, 2.

Здесь E, v, ac — соответственно модуль Юнга, коэффициент Пуассона и коэффициент температурного расширения связующего материала; Ek — модули Юнга k-го семейства волокон. Напряжения atj должны удовлетворять уравнениям равновесия

ai¿,i + a¿2,2 = bi, i = 1, 2. (5)

Правые части в (5)

bi = _ ^(1 _ 0)рс + £ Wfcр^ Ft

являются компонентами массовой распределенной нагрузки по направлениям прямоугольной декартовой системы координат. Здесь рс, рк — массовые плотности материалов связующего и волокон k-го семейства; Ft — компоненты удельной распределенной нагрузки, действующей на единицу массы.

К системе (1)-(5) присоединяются граничные условия на контуре. Уравнение контура Г задано в параметрическом виде: x = ^(s), y = 'ф(в), s — некоторый параметр. Причем Г = Гр U Ги. На контуре Гр заданы статические условия с нормальными и касательными усилиями pn(s), pT(s) соответственно:

ann2 + o-22n| + 2ai2niП2 = Pn(s), (022 _ an)niU2 + ^(П _ n|)= Pt(s). (6)

На другой части контура Ги заданы кинематические условия для перемещений ui, u2:

ui(r„) = u0(s), U2(ru)= u0(s). (7)

В (6), (7) u0(s), u0(s), pn(s), pT (s) — известные функции; ni = cos в, n2 = sin в, а в — угол, задающий направление внешней нормали к Г. С учетом (4) граничные условия (6) принимают вид

Uiai cos2(^i _ в) + ^2^2 cos2(^2 _ в) + ^303 cos2(^3 _ в) +

+ 0[шз(£и + V£22 _ LT) cos2 в + Ш3 (£22 + V£n _ LT) sin2 в +

+m4£i2 sin в cos в] = Pn(s), (8)

uiai sin2(^i _ в) + ^2^2 sin 2(^2 _ в) + ^303 sin 2(^3 _ в) +

+ 0m3 (£22 _ £ii + V(£n _ £22)) sin2в + 20 m.4£i2 cos2в = 2рт(s).

В (8) использованы обозначения для констант

11

LT = ac(1 + v)T, m3 = Emi, m4 = Em2, mi =--, m2 = .

1 _ v2 1 + v

Интенсивности о>2, задаются на той части контура Гш, где волокна входят в конст рукцию:

Общие ограничения для интенсивностей армирования имеют вид

) = w*(s), ) = (s), ) = Wg(s). (9)

0 < < 0.7, 0 = 1 _ _ о>2 _ о>3.

В дальнейшем при расчетах во всех случаях для удобства работы в системе производится ее обезразмеривание: декартовы координаты x, y относятся к некоторому характерному линейному размеру, напряжения — к характерному модулю Юнга используемых материалов. В настоящей работе сформулированы разрешающие системы для следующих случаев армирования композита тремя семействами волокон.

I. Пусть все три семейства волокон равнодеформируемы, т. е. £& = const, выполняется закон Гука = E£&, правая часть в (2) принимает вид = const + a^T. Разрешающая система уравнений строится следующим образом. Интенсивности o>i, удовлетворя-

ют уравнениям (1). Неизвестные компоненты деформаций входят в алгебраические (2) и дифференциальные (3) уравнения. Подстановка выражений для осредненных напряжений из (4) в уравнения равновесия (5) приводит к соотношениям

E

(1 + V)

E

(1 —V»)

0-1 а>1 cos <£idsl<£l + cos ^2+ 03W3 (00,5^3)^3^3 + EE (1+7) ( + (1—V2)

02^2(- Sin ^2)д52^2 + 01^1 (- Sin + (л ^ ^ (^£12),2 +

+ 03^з(- sin ^3)ds3^3 + --T^T(0(£ц + V£22 - LT))д = 61, (10)

+ , Л»£12) 1 + Тл-rT(0(£22 + V£11 - LT)) 2 = b2.

При дифференцировании напряжений используются соотношения (1), что позволяет упростить выражения. В уравнениях (10) введены обозначения для производных

dsfc^fc = cos ^fc+ sin ^fc^fe)2-

Здесь углы армирования — неизвестные функции.

Совокупность уравнений (1)-(3), (10) образует замкнутую систему девяти уравнений относительно девяти неизвестных: £¿¿, ^, . Удобно ввести три новых функции:

(Х y) = tg^(x ^(x y) = n/2.

Тогда (1) перепишется как

(wfc),1 + (wfc),2 + wfc ( (Zk),22 - (Zk)'12 ) = 0. (11)

V1 + zk 1 + zk /

Алгебраические уравнения (2) запишутся так:

£11 + zk£22 + 2zfc £12 = (1 + z2). (12)

В (10) соотношения для производных также преобразуются:

sindsfc^fc = sgnZfc . +k 2(^fc'1 + Zfc^fc'2),

1 + zk

cos ^fc dsfc ^fc = sgn zfc—■—2 ((zfc) 1 ^fc'1 + (^fco).

1 + z;

Функция sgn = ±1 в зависимости от знака тангенса угла армирования = tg , значения углов армирования содержатся в интервале 0 < < п. Причем

Ь = 1+1, Ь = .

С учетом приведенных формул соотношения (10) запишутся как

Z2 / \ Zi , Л

^ (z2,1 + z2 z2, 2) - , 2Л2 (z1,1 + Z1Z1,2) +

(1 + Z2 ) (1 + Z1 )

Е ¡3

+ (1+^) (°£12)'2 - ^3^ (1 + ¡2)2 (-3,1 + ¡3¡3,2) + Е

+ Тл-^ №п + ^£22 - ЬТ)) 1 = 61, (13)

(1 - У )

а1Ш1 (ТГ^ (^+ + а2"2 (¡2+ "2"2'2) +

Е л ¡3 / 1 \

+ тг;—^(°£12) 1 + —^ —¡3,1 + ¡3-3 ,2 +

Е

+ 7,-^ (0(£22 + - Ьт)) = 62.

В (13) компоненты деформаций находятся из (2) через

222222 Д — —¡2 ¡3 — ¡1^2 + ¡1^3 + ¡3 ¡2 — ¡3^1 + ¡2^1,

0 ^2 \ -0 ^/2 ^/2 .0 ^/2 ^/2 .0 .~2 ^2 I /—0 .~2 ^/2 .0 ^/2 ,~2

_ —£2z3z1 + £2^3¡2¡1 — £1^3¡1 — ¡2 + ¡2 — £3^3

£11 — д +

£0 ^/2 .~2 -0 ^/2 | .0 ,-^2 | .0 ¿2 ,—-0 ^2 \ -0 ,~2 + д .

Для остальных компонент тензора деформаций выражения аналогичны. После их подстановки в (13) и (3) получим три уравнения относительно трех неизвестных функций ¡1, ¡2, ¡3, причем из (13) — два нелинейных дифференциальных уравнения первого порядка, из (3) — нелинейное дифференциальное уравнение второго порядка.

Рассмотрим частные случаи. Если потребовать равенства деформаций в волокнах £1 — £2 — £0 — £0 (£0 — заданная функция), то решение системы (2) принимает вид

£12 — 0, £ц — £0, £22 — £0.

В этом случае уравнение совместности деформаций — уравнение Лапласа для задавае-0

мои £

£,11 + £,22 =

Из уравнении (10) и (11) должны определиться о>2, о>3, z1, z2, z3, что возможно при наложении дополнительных условии на один из углов армирования. Например, предполагается заданным, тогда интенсивность армирования вдоль первого семеиства волокон определяется из (11) при k =1. В частности, при = x интенсивность находится по

, F (y + ln(cos x)) F1(sin y/ex)

формуле = -, при = y — по формуле = -. Ироизволь-

cos x ex

ные функции F, F1 определяются из граничных условии. Для вычисления о>2, о>3, z2, z3

получим следующую систему нелинейных дифференциальных уравнении:

(z2) ,2 (z2),

М,1 + z2(^2),2 + ^2 ( - Z2 1 + ^

/Ni /Ai I (z3) ,2 (z3),1

Ы,1 + z3(^3),2 + ^3 ( - Z3Y+Z2

(Z2,1 = Z2 Z2,2) - (¿3,1 + ¿3*3,2) + (14)

(1 + z2 ) (1 + z3 )

+ EKi ( ) EKi +

+ ^-2 (-^2,1 - <^3д) = Z-2 ,1 + bl,

1 - V2 1 - V2

^ (Т^ (Z2 ^ =Z2Z2 2) + -3-3 (Т^ (z3'1 + Z3Z3'0 +

+ EK2 ( ) EK2 + b

+ ^-2 (-^2,2 - ^3,2) = i-2,2 + 02,

1 - V2 1 - V2

где

K = £ц + V£22 - LT, K2 = £22 + V£n - LT.

Правые части в (14) — известные функции.

II. Пусть равнодеформируемые семейства волокон расположены так, что углы армирования, например — известные функции, а третий угол армирования — произвольная функция координат. Математически задача при такой укладке становится переопределенной, и в случае построения решения необходимо вводить произвол, чтобы выполнить все уравнения модели (1)-(3) и уравнения равновесия (5). Из (12) выразим компоненты деформаций £11, £22 через £12:

£11 = вц£12 + б12,

£22 = e21£12 + e22, В дальнейшем используются обозначения

в11 = -02/01 tg2 - 2 tg ^1, б12 = £1 (1 + tg ^1) - 03/01 tg2 ^1,

e21 = 02/01, 622 = 03/01, (15)

где коэффициенты 01,02,03 вычисляются через заданные углы армирования и условия равнодеформируемости по формулам

01 = 1 - tg2 ^1/ tg2 ^2, 02 = 2(tg ^1/ tg2 ^2 - 1/ tg ^2),

03 = (£0(1 + tg ^2) - £0(1 + tg ^1))/ tg ^2. (16)

Из третьего уравнения (12) компонента деформации £12 определяется через неизвестную функцию Z3 = tg ^3:

_ b12 + b22z3 П7ч

£12 = -:-^—, (17)

en + e21z| + 2z3

где b12 = -e12£^, b22 = £0 - e22. С учетом (15) уравнение совместности деформаций примет вид

e11£12,22 + e21£12,11 = 2£12 ,12. (18)

Знак детерминанта 5 = e11e21 - 1 уравнения (18) определяет его тип при заданных значениях углов армирования

5 = -(z2 - z1)2 (19)

5 (Z1 + Z2)2. (19)

Пусть = arctg x, = arctg y. Решение уравнения (18)

x3

£12 = ^ + (1 - ln y)y + C1 (x + y) + C2 (20)

о

содержит две произвольные константы Ci, C2. Подстановка найденного решения (20) в граничные условия (8) позволяет решить две задачи: 1) установить ограничения на нагрузку; 2) для заданной нагрузки найти соответствующее уравнение граничного контура. Решаем задачу 2. Не ограничивая общности, предположим, что pT = 0, pn — заданная функция. Пусть угол ß, определяющий контур, неизвестен. После подстановки компонент деформаций £11, £12, £22 в (8) получим систему двух уравнений относительно неизвестных значений на граничном контуре и ß. Из первого уравнения (8) выразим

COS2(^3 - ß) = [pn - WiO"! cos2(^1 - ß) - ^2^2 COS2(^2 - ß) --0(шз (en + V£22 - LT) COS2 ß - m.3 (£22 + V£n - LT) sin2 ß -

-m4£12 sin ß cos ß)]/a>30"3. (21)

Правую часть обозначим Cp. Во второе уравнение входит sin2(^3 - ß), с учетом (21) и известных тригонометрических формул запишем

sin 2(^3 - ß) = ±2CHI -1 - 1.

После подстановки найденных значений второе уравнение (8) зависит от неизвестной ß:

P1 cos2 2ß + P2 sin2 2ß + P3 cos 2ß sin 2ß + p4 cos 2ß + p5 cos 2ß + p6 = 0. (22)

Коэффициенты p1, p2, p3, p4, p5, p6 выражаются через известные значения ш1, ш2, о>3, ^1, на граничном контуре. Применим универсальную тригонометрическую подстановку — выразим sin2ß, cos2ß через tg ß и получим полное уравнение четвертой степени относительно tg ß c новыми коэффициентами p*, i = 1, 6, определяемыми через p¿:

p* tg4 ß + p2 tg3 ß + p3 tg2 ß + p* tg ß + p* = 0. (23)

Уравнение —X = tg ß определяет уравнение контура y = f (x) в декартовой системе координат, отвечающее решению (20).

III. Пусть все армирующие семейства волокон нерастяжимы. Тогда в (2) правая часть имеет вид £°к = . Однако в силу идеализации модели волокон (абсолютно твердое тело) теперь а1, а2, а3 — неизвестные функции координат Поэтому в (10) в случае нерастяжимых волокон появятся производные от этих функций:

0"1;1 (о>1 cos2 ^1) + 01)2(о>1 cos sin ^1) + СТ2)1(^2 cos2 ^2) + +^2,2(^2 cos ^2 sin ^2) + ^3,1(^3 cos2 ^3) + ^3,2^3 cos ^3 sin ^3) --^1^1 sin - ^2^2 sin - ^3^3 sin +

E E

+ (T+V) ("£12 ),2 + (1 - v 2) (°(£11 + V22 - LT )),1 = b1,

a1;1 (ш1 cos sin + а1,2(ш1 sin2 + a2;1(a>2 cos sin ^2) +

+^2,2(^2 sin2 ^2) + ^3,1(^3 cos ^3 sin ^3) + ^3,2(^3 sin2 ^3) +

+^1^1 cos + ^2^2 cos + ^3^3 cos +

EE + (1 + ^) ("£12)'1 + (1 - V2) ("(£22 + V£11 - LT)),2 = 62.

Совокупность уравнений (1)—(3), (24) не является замкнутой разрешающей системой относительно названных переменных ец, ^, , ок, и для решения задачи следует вводить дополнительные предположения, например о способах укладки семейств армирующих волокон.

Приведем частные решения в случае нерастяжимых семейств волокон. Рассмотрим возможные "особые случаи", возникающие при вычислении компонент деформаций через тангенсы углов армирования из системы алгебраических уравнений (12), и укажем, какие при этом получаются решения исходной задачи.

А. Правая часть (12) е^ = акТ, поэтому при нулевой температуре система становится однородной. Условие ее разрешимости Д = 0, а именно

1 zf 2z1 1 z12 z1

A— 1 zf 2z2 —2 0 (zf - zf) (z2 - z1) — 0,

1 zf 2z3 0 (zf - zf) (z3 - z1)

приводит к требованию, чтобы траектории двух семейств совпадали: ^ = ¿2 или ^ = ¿3.

Б. Пусть Т = 0. Систему (12) можно рассматривать как систему трех линейных неоднородных уравнений относительно трех компонент деформаций ец. При выборе способа укладки семейств волокон арматуры необходимо рассмотреть возможность Д = 0 и одновременно Дк = 0. Например, для к = 1

Ai

£0(1 + zf) zf zi £0(1 + zf) zf Zf £0(1 + zf) zf za

-2

za zf £0(1 + zf)

zi zf £0(1 + zf)

0 0 £0 (1 + zf) - £1(1 + zf)

Следовательно, Д = 0, Дк = 0 при условии, что деформации в первом и втором семействе волокон равны и = ¿2, либо при условии, что равны деформации в первом и третьем семействе волокон и = ¿1, либо равны деформации в третьем и втором семействе волокон, а ¿з = ¿2.

Пусть ¿1 = ¿2. Решение системы (12) неединственно, выражения для компонент деформаций через неизвестную функцию запишутся в виде

е1^1 — 2е12^1 — к11

£22 —

zi z4

£ii — kii - zf£22j £i2 — const. Коэффициент kii вычисляется по формуле

k11 — £3zf — 2z3£12 + £0.

Уравнение совместности (3) в этом случае примет вид

— zf £22,22 + £22,11 — 0.

После вычисления производных

dii z3 + d12zl + d13z1 + d14

(25)

£22,11 —

£22,22

(zf - zf)4

+ d 13zi + di 4

(zl - z3)4

( )2 + P 11+ P 12 z (zi, 1) + Tzi z2)lz1, 11,

( z1 - z3 )

(z )2 + P11 + P12 z

(z1,2) + (zi ziTl ,22 ( z1 - z3 )

из (25) получим уравнение для неизвестной

¿11^ + ¿ц¿1 + ¿13^1 + ¿14 ( )2 + Р11 + Р12 _

/ 2 2\4 \г1,1) + / 2 2\2 ¿1,11

(¿I2 _ ¿32)4 (¿I2 _ ¿32)2

21 ¿11^1 + ¿12^1 + ¿13¿1 + ¿14/ х2 р11 + р12 I _п

-7 2-2Т4-(г1,2) + Т~2-¿1,22 ) — (26)

\ ^ _ ¿32)4 (¿Г _ ¿32)2

в разрешенном относительно старших производных виде —

B(zf)(z1, 1 ) + B(zf)

Z1, 11 - Z32Z1,22 + (zi, 1 )2 + (zi,2)2 = 0. (27)

Знак детерминанта S = —z3 — 1 < 0 для уравнения (27) определяет его гиперболический тип для любого значения заданного угла армирования В (27) использованы следующие обозначения:

A(z3) = dnz3 + ¿^zj2 + d13Z1 + ¿14, B (zf) = C11Z3 + C12Z4 + C13Z3 + C14Z2 + C15Z1 + C16, P11 = 2kn — 2e0z3, P12 = 2s12z2,

¿11 = —4p12, ¿12 = —2P11Z3 + 4pnz2, ¿13 = 4p12z3, ¿14 = puzf,

C11 = P11, C12 = P12, C13 = — 2pnzf, C14 = —2p12Z2, C15 = P11 z|, C16 = P12Z4.

При дополнительных условиях равнонапряженности третьего волокна о3 = const и задании угла армирования получена разрешающая система четырех уравнений относительно четырех неизвестных z15 01 а1, а2. Она состоит из (27), (12) при k =1 и двух уравнений (24), которые при сформулированных условиях упрощаются к виду

(^1,1 + 0~2дМ ^ 2 + sgn Z1(^1,2 + ^2,2)0^-^^-2

1 + Z1 1 + Z1

— sgn Z1 (0 + 02)01 Z1 2.2 (Z1,1 + Z1Z1,2) = (1 + z!2 )2

E

= b1 — Тл-(°(£11 + V£22 — LT)) 1,

(1 — v )

sgn Z1(^1,1 + ^2,1)0^ ^ 2 + (^1,2 + ^2,2)0^ Z1 2 +

z1

Л_ f±

(1 + z?)4 Z1'

E

= b2 — --^(0(£22 + V£11 — LT)) 2.

(28)

1 + z2 ' 1 ->^1 + z2

+ sgn z1(^1 + 0-2)01 ^ z1 ( — z1,1 + z^) =

(1 _ V2)

Система распадается, вначале находится ¿1 из (27), затем с его использованием определяется из (11), затем осуществляется переход к решению (28).

При нахождении частных решений нелинейного уравнения (27) используется алгоритм построения инвариантных решений уравнений в частных производных на основе групп преобразований, допускаемых данными уравнениями, описанный в [4]. Этот алгоритм строит определяющие уравнения для изовектора с использованием форм Картана. Уравнение в частных производных преобразуется к эквивалентной системе дифференциальных

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

Например, при = п/4 и постоянных значениях уравнение (27) запишется как

4г3 _ 2г2 — 4г — 1

— ^22 — + Ц — 2^ — + г, + 1 ((^>2 + (г1,2>2) = (29)

После применения алгоритма получено частное решение

X

г(х, у) = + у — х) + к, (£) + к2(£ + у — х))^ + ^\(у — х).

о

Произвольные функции к,(х), к2(х), находятся из граничных условий. Граничная задача для г, ставится на части контура, где волокна входят в конструкцию. Используются граничные условия (8). При принятых в данном пункте предположениях и в = п/2 они запишутся так

z

2

(Ш101 + Ш202) 1 ^ 2 + (1 - Ш1 - Ш2 - Ш3)т3 1 + z1

2£lzl - 2£12z1 - k11

£lzl - 2£12z1 - k11 + zi - z32

+ V kn - z

zi z3

pn(s) - 0.5Ш303,

(Ш101 + Ш202) sin 2^1 = 2pr (s) + 2(1 - Ш1 - Ш2 - Ш3)т4 + 03Ш3.

После нахождения ш101 + ш202 из второго уравнения и подстановки в первое получается кубическое уравнение относительно z1 на рассматриваемом контуре, оно имеет вид

0.5Bz1(zj! - zf) + Om3(£1z2 - 2£12z1 - k11 + vk11(zj! - z|) --z32(£?z2 - 2£12z1 - kn)) = (p„(s) - 0.5)(z2 - zf), B = 2pT + 2Qm4£12 - ш303.

Выберем в качестве связующего материала алюминий (E = 71 ГПА, v = 0.31), пусть он армирован стальными волокнами (E = 200 ГПА, ak = 10.510-6 K-1) [5]. Температура T = 10, интенсивности на границе ш1 = 0.1, ш2 = 0.1, ш3 = 0.3. Если зададим нормальное давление = 6.203, уравнение относительно z1 примет вид

0.2694z3 - 2.3034z? - 0.8114z1 = 0.

его корни 0; 8.8861; -0.3388. Приведем частное решение для (29), удовлетворяющее, например, первому корню (z1 = 0) на границе прямоугольной пластины размера x е [п/2,п], y е [п/2, 3п/2]:

- cos(x + y) 3cos(x - y) z1(x, y) = -2--cos x + sin y +--2--+ 1 + sin(x - y).

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

Список литературы

[1] Nemirovsky Ju., Feodorova N.A. The mathematical analysis of permiting systems of the equations of a flat problem of the reinforsed environments // KORUS-2002. Proc. of the 6th Intern. Symp. on Science and Technology. Novosibirsk, 2002. P. 195-197.

[2] Nemirovsky Ju., Feodorova N.A. Flat problem of the elastic environment reinforced with three families of fibres // K0RUS-2005. Proc. of the 9th Intern. Symp. on Science and Technology. Novosibirsk, 2005. P. 506-511.

[3] Nemirovsky Yu.V. On the elastic-plastic behaviour of the reinforced layer // Intern. J. Mech. Sci. 1970. Vol. 12. P. 898-903.

[4] Олвер П. Приложения групп Ли к дифференциальным уравнениям. М.: Мир, 1989. 635 с.

[5] Композиционные материалы: Справочник / В.В. Васильев, В.Д. Протасов, В.В. Болотин и др. М.: Машиностроение, 1990. 512 с.

Поступила в редакцию 2 ноября 2005 г.