Применение методов теории возмущений в задачах равнонапряжённого армирования изгибаемых металлокомпозитных пластин, работающих в условиях установившейся ползучести Текст научной статьи по специальности «Математика»

Механика деформируемого твёрдого тела

УДК 519.376

ПРИМЕНЕНИЕ МЕТОДОВ ТЕОРИИ ВОЗМУЩЕНИЙ В ЗАДАЧАХ РАВНОНАПРЯЖЁННОГО АРМИРОВАНИЯ ИЗГИБАЕМЫХ МЕТАЛЛОКОМПОЗИТНЫХ ПЛАСТИН, РАБОТАЮЩИХ В УСЛОВИЯХ УСТАНОВИВШЕЙСЯ ПОЛЗУЧЕСТИ

А. П. Янковский

Институт теоретической и прикладной механики им. С. А. Христиановича СО РАН,

Россия, 630090, г. Новосибирск, ул. Институтская, 4/1.

E-mails: lab4nemir@rambler.ru

Методами теории возмущений решена задача равнонапряжённого армирования изгибаемых металлокомпозитных пластин, 'работающих в условиях установившейся ползучести. Показана возможность существования нескольких альтернативных решений 'рассматриваемой задачи, которые можно надёжно выделять, используя разработанный алгоритм. Построены конкретные проекты равнонапряжённого армирования двусвязной пластины при разных плотностях армирования на внутреннем контуре.

Ключевые слова: изгибаемые пластины, установившаяся ползучесть, рациональное проектирование, равнонапряжённое армирование, малый параметр, теория возмущений.

Введение. Настоящая работа продолжает исследования, опубликованные в [1,2], где была сформулирована задача равнонапряжённого армирования (РА) изгибаемых металлокомпозитных пластин, работающих в условиях установившейся ползучести всех компонентов композиции. В [1] показано, что система разрешающих уравнений такой задачи РА является квазилинейной и относится к смешанно-составному типу [3]. Особенности этой системы уравнений и соответствующих ей нелинейных граничных условий, выявленные в [1,2], не позволяют надеяться на построение широкого круга аналитических решений и требуют разработки обоснованных численных процедур решения, причём указанные особенности не позволяют произвести и простой перенос известных методов численного интегрирования квазилинейных систем уравнений [4,5] на решение рассматриваемой задачи. Наличие в задаче естественного малого параметра (см. ниже) позволяет надеяться на эффективность применения хорошо зарекомендовавших себя методов теории возмущений [6 и др.].

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

Андрей Петрович Янковский (д.ф.-м.н., проф.), ведущий научный сотрудник, лаб. физики быстропротекающих процессов.

1. Система разрешающих уравнений и граничные условия. Для удобства дальнейшего изложения систему разрешающих уравнений задачи РА и соответствующие ей граничные условия [1] запишем в безразмерном виде, предполагая, что в равномерно нагретую пластину постоянной толщины 2Л, внедрено N семейств армирующих проволок постоянного поперечного сечения. (Все размерные функции и величины будем помечать сверху волной, например /г., а соответствующие им безразмерные функции и величины — обозначать теми же символами, но без волны.) Тогда замкнутая система разрешающих уравнений, содержащая в качестве неизвестных функций безразмерную скорость прогиба установившейся ползучести и и параметры РА сок, фк, Ак, состоит [1] из уравнения равновесия

[дп(грк,Ак) - А2к/шк] — еВ(у, со) = -р(х1,х2); (1)

к

условий постоянства поперечных сечений волокон

д.3(фк,шк)+шкдп(фк,фк) = 0, к = 1,2,... , АГ; (2)

условий равнонапряжённости арматуры

дз(фк,д3(фк,у)) ~ Акдп(фк, ,у)/шк = -£% = сог^ ф 0, к = 1,2,..., Ж, (3)

где

Ак = шкд3(фк,фк), к = 1,2,... ,N]

(4)

B{v,w) =

i=l,2

+

, гг

+ 2

l-^WfcU№ lV,n

,12

U = 3 - г),

дз(Фк, • ) = (• XicosV’fc + ( • ),2Sin фк,

дп(Фк, ■) = д8(фк + п/2, •) = -(•) ,1 sill + (•),2 cos фк (1 ^ к ^ N),

(5)

К = 2У у2,П + ^11^,22 + ^22 + v2,12> Ш = {^1, W2, • • • , 10N} 5

суммирование по к производится от 1 до N] нижний индекс после запятой означает частное дифференцирование по соответствующей переменной Хі, г = = 1,2.

Пусть область G, занимаемая пластиной в плане, ограничена контуром Г, тогда на одной части этого контура (обозначим её Гр) могут быть заданы статические [1]

^2 а°кшк cos2(фк -/?) — eDM(v, со) = Мп, к

-^2a°kAksm(tpk ~ /3) + 0,5^2 aldT(u)ksm2(tpk -/?))- (6)

- eCF(v,co) - eDF(v,co) = Fnz, (жі,ж2) Є Г.

pi

где

Ом{у, ш) = (і - ^ шк) 1 [(2у,п +у,22) сое2 /3+

к

+ (2-У,22 + У,1і) єіп2 /З + У,12 віп 2/3] ,

Ср(у,ш) = ^ Пі| (і - '^2шк^^°~1(2у,іі + ь,м) + і=1,2 к ,г

+ (і - ^^)кМ0“Чі2 .}, (7)

к ,3

Ор(у,ш) = ^ (-1)г+1Пі|(^1 0,5(г?,22 - у,п) віп 2/3+

г=1,2

+ V,і2 сов 2/3] > (І = 3 — г), пі = сов/3, пг = єіп/З,

із

дт(-) = 9П(Р, •);

на другой части (обозначим её Г^) контура Г — кинематические граничные условия:

у(Т^=у0, д3((3,у) = вп, [хі, ж2) Є Г„. (8)

Кроме того, на части контура Г (обозначим её Гд.), на которой волокна к-того семейства «входят» в пластину, необходимо задать краевые условия для интенсивностей армирования [1]:

^кі^к) ^Ок, к 1,2,..., N.

(9)

Решение краевой задачи РА (1)—(4), (6), (8), (9) должно удовлетворять физическим ограничениям [1]

(10)

В уравнениях и соотношениях (1)—(10) использованы следующие формулы обезразмеривания и обозначения:

р 0 /? 7?

О _ к лО _ <*к с0 — Д (л°\тк —

к /ъ і „ лхО’ ^>к Гп’ ^к к к) > ^

(2 + /^)5-0’

і

Хі = р

, Мп =

м„

Ь2р

к2а°г кЧ°

р —

> 1 ПХ —

.ок.

кщ кв,

У0 = вп =

о’

є =

4 Ме?

2 + /Хо <7?

МО

г = 1,2, Л = 1,2,...,ЛГ;

4 ^о(С?

МО

4 Б0

(П)

(12)

Во = (\/Зто+1Во) = const, В к = В к^к = const,

1 1 ^13)

НО = —, Цк = — (1 ^ к ^ N);

т0 тк

Во, то, Вк, mu — известные характеристики (константы) материалов связующего и арматуры fc-того семейства, определяющие степенной закон установившейся ползучести = Бга(Т™п ({„ — скорость деформации, ап — напряжение, 0 ^ п ^ N) фазовых материалов композиции; v — скорость прогиба установившейся ползучести; р— заданная распределённая поперечная нагрузка; Х\, Х2 — прямоугольные декартовы координаты точек срединной плоскости пластины; D — характерный размер пластины в плане; <т^ = const — заданное значение напряжения в равнонапряжённой арматуре fc-того семейства на верхней лицевой поверхности пластины (жз = h)] ^ = const — соответствующее значение скорости деформаций установившейся ползучести арматуры fc-того семейства; шк, фк — плотность и угол (отсчитываемый от направления х\) армирования проволокой fc-того семейства; w* = const — предельно допустимая суммарная плотность армирования (на практике w* ~ 0,7); Мп, Fnz — изгибающий момент и приведённая поперечная сила Кирхгофа, заданные на Гр; /3 — угол, задающий направление внешней нормали к Г и откладываемый от направления х\] Vo, дп — скорость прогиба на и производная от скорости прогиба по направлению внешней нормали к контуру Г^; Wqк — заданные на Г& плотности армирования.

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

Пусть металлокомпозитные пластины получены путем внедрения в медное (Си) или алюминиевое (АМгЗ) связующее стальной проволоки У8А, тогда при температуре Т ~ 200 °С характеристики фазовых материалов таких композиций (см. (13)) имеют следующие значения [7-9]:

Си

т0 = 2,16, В0 = 5,63 • 10“у (МПа)

-9 /'МГГсЛ-"г0 . т1-1

АМгЗ : т0 = 6,045, В0 = 1,025 • 10“14 (МПа)“т° • ч"1, (14)

У8А : т\ = 24,98, Вх = 1,054 • 10“84 (МПа)”™1 • ч-1.

На рис. 1, а изображены обратные степенные зависимости

Ск{1к) = (В~1£к)**к, к = 0,1

при параметрах ползучести (13), (14). Кривая 1 соответствует меди, линия

2 —сплаву АМгЗ, кривая 3 —стальной проволоке У8А. На рис. 1, б изображены зависимости в(^) (см. (12)—(14)), где кривая 1 соответствует Си-У8А-композиции, а линия 2 - АМгЗ-У8А-композиции.

Как видно из рис. 1, б, в достаточно большом интервале изменения € (0,10-4] (ч-1) значения параметра £ для рассматриваемых металлокомпозиций имеют порядок 10-2, т. е. являются малыми. При этом величины

<т, МПа

є

1200

400

800

3

0,02

0,04

0,03

0,01

0 0,2 0,4 0,6 0,8 | • 104,ч-1

О 0,2 0,4 0,6 0,8 и • 104, ч-1

б

а

Рис. 1. Степенные аппроксимации установившейся ползучести некоторых фазовых материалов (а) и зависимости малого параметра от скорости деформаций позучести в равнонапряжённой арматуре (б)

(входящие в те же уравнения (1), (6), что и е), согласно (11), (13), (14), для арматуры из проволоки У8А равны 0,98, т.е. сомножители |сг£| |ь|. Следова-

тельно, параметр е в уравнениях (1), (6) действительно можно рассматривать как малый. (Подобные оценки малости параметра е можно получить и в случае использования в качестве арматуры вольфрамовых проволок [10 и др.]) В [1] показано, что система разрешающих уравнений (1)—(4) является системой квазилинейных уравнений смешанно-составного типа [3], причём траектории армирования, определяемые углами ф^, являются её двукратными действительными характеристиками. Структура уравнений (1)—(4) и граничных условий (6) (с учётом соотношений (5), (7)) показывает, что задачи определения параметров РА ф^, шк и скорости прогиба установившейся ползучести V связаны и решать их необходимо совместно. Последнее обстоятельство существенно усложняет разработку методов решения нелинейной краевой за-

2. Применение методов теории возмущений. В настоящее время теория систем квазилинейных уравнений смешанно-составного типа разработана недостаточно полно [3], чтобы можно было аналитически исследовать свойства решений системы разрешающих уравнений (1)-(4) в общем случае, поэтому целесообразно воспользоваться какими-либо упрощениями этой системы и соответствующих ей граничных условий.

Важная особенность системы разрешающих уравнений задачи РА состоит в наличии в уравнении равновесия (1) и в статических граничных условиях (6) малого параметра е. Это позволяет разыскивать решение задачи РА изгибаемых пластин при помощи методов теории возмущений [6 и др.], например метода разложения решения в виде асимптотических рядов по степеням малого параметра е [6] или нижеследующего итерационного процесса (который качественно схож с итерационным методом определения основного момент-ного напряжённого состояния в теории оболочек [11]).

Пусть ф[*\ А^\ (1 ^ к ^ Ы)— в-тое приближения неизвест-

ных функций, тогда (в + 1)-е приближения получим путем интегрирования уравнений

дачи РА (1)—(4), (6), (8), (9).

Лг)

дв(.‘Фк)^'Фк)) = ~Тг)’ к = 1,2,... ,ЛГ; (16)

шк

д3{ф{^) = 0, к = 1,2,... ,ЛГ; (17)

^(4Г) > ^(Г)) - ^ 9га(^г) ^ = СОП^ + °> П

“к (18)

д3{ф([\у^) = у1г\ к = 1,2,... ,ЛГ, (Ж1,Ж2)еС

при граничных и краевых условиях

£*МГ) сое2(^г) -/?) = Мп + еОмЬ^,и('з]), к

-'52°кАк) 8[п-Р) = РП2+еСР(у(8\шМ)+еВР(у(8\шМ)-

к

-0,5^о-^г(4г)8т2(^г) -/3)), (жьж2) € Гр; к

У{'г){х1,х2) = ^о(Г^),

у1г\хг,х2) = 8т(^г) - /3)дт(у0) +сов(ф^ — /3)0П(Г„), (жьж2) € Г*,;

(19)

(20)

Ш^г)(ж1,ж2) = шок(Тк), (Х1,Х2) етк, к = 1,2,..., АГ, г = в + 1, (21)

где —введённые для удобства функции; о/8) —векторная функция, аналогичная ш (см. (5)).

Для начала итерационного процесса нулевое приближение выберем в виде

у{0) = 0. (22)

Уравнения итерационного процесса (15)—(21) получаются из соответствующих уравнений разрешающей системы следующим образом. В уравнении равновесия (1) и в статических граничных условиях (6) в операторы В, Им, Ср, В>р, имеющие в качестве сомножителя малый параметр е, подставляются приближения неизвестных функций, определённые на предыдущей итерации. Поэтому слагаемые, имеющие сомножитель е, на каждой итерации известны и в (15), (19) перенесены в правую часть. Задание начального приближения в виде (22) аналогично отбрасыванию в (15), (19) на первой итерации (г = 1) слагаемых, содержащих е. Остальные уравнения итерационного процесса на каждой итерации формально совпадают с соответствующими соотношениями разрешающей системы уравнений (см. п. 1).

Проанализируем систему (15)—(18) при различных количествах семейств арматуры (-/V). Пусть в пластину внедрено одно семейство волокон (Ж = 1). Тогда три квазилинейных уравнения (15)—(17) замкнуты относительно трех

неизвестных функций ф^\ \ определяющих г-тые приближения па-

раметров РА. Если эти функции известны из каких-то соображений, то два линейных уравнения (18) замкнуты относительно двух функций у^г\ У^г\

определяющих г-тое приближение скорости прогиба установившейся ползучести. Аналогично распадаются граничные и краевые условия (19)—(21): условия (19), (21) состоят из трех равенств и замкнуты относительно контурных

значений функций ф^, , с, кинематические же условия (20) задают на

кромке Г„ контурные значения функций У^г\ Следовательно, при N = 1 итерационный процесс (15)—(22) позволяет расщепить ранее связанные задачи определения параметров РА и скорости прогиба в поперечно изгибаемой пластине на отдельные подзадачи, которые можно последовательно интегрировать на каждой итерации.

Если в пластину внедрены два семейства арматуры (Ж = 2), то система (15)—(18) в целом замкнута относительно неизвестных функций фк\ А^\ \ у(г\ Ук\ но её подсистема (15)—(17), содержащая приближения только

параметров РА ф^\ А^\ (к = 1,2), не замкнута, так как состоит из пяти уравнений и содержит шесть неизвестных функций. Следовательно, при N = 2 система (15)—(18) не распадается столь же просто, как это было при N = 1, на несвязанные подзадачи определения приближений параметров РА и скорости прогиба. При исследовании типа системы (15)—(18) выясняется, что её характеристический определитель тождественно равен нулю. Это обстоятельство вызвано тем, что два последних уравнения (18) содержат частные производные только от одной функции у(г\ а значит, должно выполняться некоторое условие совместности этих уравнений, которое можно получить за счёт перекрестного дифференцирования двух последних равенств (18) по

направлениям фт (т = 3 — к) [12]. Однако такое условие совместности со-(V) (г)

держит функции фк , Уки их частные производные, т.е. преобразованная соответствующим образом система (18) совместно с (15)—(17), по-прежнему не расщепляет задачу определения приближений скорости прогиба и параметров РА на отдельные подзадачи. Следовательно, при N = 2 применение итерационного процесса для решения задачи РА изгибаемых пластин приводит лишь к упрощению структуры разрешающих уравнений и граничных условий, но качественные особенности системы уравнений итерационного процесса остаются прежними: для каждого г ) 1 в (15)—(22) задачи определения приближений параметров РА и скорости прогиба связаны между собой; соответствующие системы уравнений имеют составной тип [13]; система (15)—(18) является квазилинейной. Все это существенно усложняет разработку методов интегрирования уравнений итерационного процесса по сравнению со случаем N = 1.

Если в пластину внедрены более двух семейств волокон (Ж ^ 3), то (по тем же причинам, что и при N = 2) характеристический определитель системы (15)—(18) тождественно равен нулю. Однако в этом случае можно легко получить условия совместности подсистемы (18), содержащие только приближения параметров РА, что позволяет замкнуть подсистему (15)—(18) относительно функций фк \ А^\ (1 ^ к ^ И), а значит, и расщепить исходную

задачу РА на её составляющие, как при N = 1. Получим указанные зависимости, переписав условия РА (18) в исходной форме (см. (12) в [1]):

^111к11к1 + У%1к2 1к2 + 2у%1к11к2 = = соп^ ^ °> ю < ЛГ, (23)

23

i(r) Лг) Лг) • i (r) где lKk{ =cos^;, lli =S1I1

Рассмотрим (23) как систему линейных алгебраических уравнений относительно v^j. Выберем из (23) три уравнения, определитель которых отличен от нуля в некоторой окрестности точки (х\,х2)- Пусть для конкретности это будут три первых уравнения, тогда, используя правило Крамера, из них можно однозначно определить

i,j = 1,2,

(24)

где

Ап = —2det А12 = А21 = - det А22 = -2 det А = 2 det

lk2 lk2

i(r)i(r) lk 1 lk2

i(r)Ar) hkl hkl

i(r)Ar) k2 k2

i(r)Ar) Lkl Lkl

i(r)Ar) hkl hkl

i(r)i(r) k2 k2

/W/W lk 1 lk2

i(r)i(r) hkl k2

k = 1,2,3.

(25)

производные являются функциями от ф^

(Здесь использована построчная запись определителей по к.) Следовательно,

к = 1,2,3.

Если количество семейств волокон больше трех, то, подставив выражения (24) в уравнения (23) при 4 ^ к ^ -/V, получим N — 3 нелинейные конечные

зависимости между функциями фп ^ (1 ^ п ^ -/V). Кроме того, производные

у^ должны удовлетворять условиям совместности деформаций:

(v[n),2 =

(26)

Подстановка выражений (24) в (26) даёт ещё две, но уже дифференци-

(г)

альные зависимости между фк , к = 1, 2, 3. Получающиеся в итоге конечные и дифференциальные связи совместно с (15)-(17) образуют замкнутую систему уравнений относительно приближений параметров РА ф^\ А^к \ оо^к\ 1 ^ к ^ N. Следовательно, преобразованная таким образом система (15)—

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

Однако такой подход имеет следующую особенность. Пусть уже известны приближения фк^ (1 ^ к ^ -/V), тогда из (24) с учётом (25) будут известны

(г) чгз1

во всей области О, занимаемой пластиной в плане, будут известны полные дифференциалы (1у^ = у^йхг + у^2(1х2- Интегрируя эти выражения, задав для них начальные условия всего в одной точке [14], получим производные у^ , г = 1,2, которые отвечают условиям = ^21 > поэтому выражение

с1у^ = у^^Хх + у('т^Х2 будет полным дифференциалом функции Интегрирование последнего выражения (при задании начального условия всего

производные у\£, удовлетворяющие условиям совместности (26). А значит.

в одной точке) позволит определить во всей области О г-тое приближение скорости прогиба Очевидно, что случаи закрепления пластины, когда кинематические граничные условия для скорости прогиба и её частных производных первого порядка заданы всего в одной точке (или каждое в своей точке) кромки, практически не встречаются.

Таким образом, проведенное рассуждение показывает, что при внедрении в пластину более одного семейства арматуры итерационный процесс (15)—(22) не приводит к принципиальному упрощению исходной системы разрешающих уравнений задачи РА.

Проанализируем более подробно уравнения итерационного процесса (15)-(22) в случае укладки одного семейства равнонапряжённой арматуры (Ж = 1). Характеристическое уравнение подсистемы (15)—(17), замкнутой относительно функций ф^\ \ имеет вид

— х'2 сов •01^)2(ж2 + сое ф^) = 0, Г = в + 1, (27)

где производная х'2 = (1х2)(1х\ задаёт направление характеристики.

Равенство (27) означает, что система (15)—(17) имеет две ортогональные действительные характеристики, причём одна из них, задаваемая углами

(г)

ф\ , является двукратной. Наличие только действительных характеристик у этой системы указывает на её сходство с системами уравнений гиперболического типа, рассмотренными в [15]. Однако если уравнения (15), (16) записаны в инвариантах Римана, то уравнение (17) не удаётся ни записать в инвариантах Римана, ни привести к характеристической форме [15], что не позволяет отнести систему (15)—(17) к гиперболическому типу (в смысле определения, введенного в [15]). Кратность же одной из характеристик указывает на близость этой системы к системам параболического типа [16].

Действительно, если задать некоторые значения функций ф^\

на двукратной характеристике системы (15)—(17), определяемой углами ф^\ и попытаться определить частные производные первого порядка от этих функ-

/ (г)

ции, то окажется, что от щ эти производные можно определить непосредственно из системы (15)—(17), а для определения производных от А^\ с по направлениям, не касательным к этой характеристике (например, по нормальному направлению д^ф^, •)), необходимо дважды продифференцировать спо направлению двукратной характеристики, т.е. необходимо использовать вторую производную д2(ф^\ Сс1^)- Известно, что подобными особенностями обладают системы уравнений параболического типа [16].

Если же задать значения функций ф^\ на характеристике, Опре-

ем)

деляемой углами ф\ +7г/2, и вновь попытаться определить частные производные первого порядка от этих функций, то окажется, что от ф^\ сэти

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

Следовательно, система (15)—(17) не подпадает ни под один из типов су-

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

Указанная особенность системы (15)—(17) приводит к тому, что характеристики системы могут быть линиями слабых разрывов. Кроме того, решение задачи Коши для нее может быть построено только в области определённости решения (об этом свидетельствуют многочисленные расчёты, проведенные автором), что характерно для решений задачи Коши для систем уравнений гиперболического типа [15].

Если функции ф^\ известны из системы (15)—(17), то система

уравнений (18) является линейной системой параболического типа, замкнутой относительно функций У^г\ характеристическое уравнение которой имеет вид

(sinф^ — х'2 cos ф^)2 = 0, г = s + 1. (28)

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

Уравнения (27), (28) показывают, что на каждой итерации подсистемы (15)—(17) и (18) имеют совпадающие двукратные характеристики, которые в силу квазилинейности подсистемы (15)—(17) изменяются от итерации к ите-

(г)

рации и определяются углами ф\ . Так как функция ф\ задаёт траектории

(г)

РА, а ф\ —r-тое приближение к ф\, то двукратные характеристики систем (15)—(17) и (18) можно трактовать как r-тые приближения к траекториям РА. Если итерационный процесс (15)—(22) сходится, то двукратные характеристики совпадут в пределе с траекториями РА, которые, в свою очередь, являются двукратными характеристиками системы разрешающих уравнений задачи РА (1)-(4) (см. (21) в [1]).

Для квазилинейных систем уравнений с действительными характеристиками наиболее естественно и просто ставится задача Коши, поэтому обратимся к вопросу постановки задачи Коши для систем (15)—(17) и (18). Эти подсистемы удовлетворяют условиям теоремы Коши—Ковалевской [18], поэтому при аналитических входных данных решение задачи Коши для них существует и единственно. Для более широких классов входных данных теоремы существования для этих подсистем автору неизвестны.

Покажем, что статические граничные условия (19) совместно с краевым

условием (21) задают контурные значения функций ф^\ Для этого будем считать, что арматура «входит» в пластину на той части контура Г, где заданы статические граничные условия (Гр = Ti). Тогда на контуре Гр из

(21) известны контурные значения функции входящей в (19). Исключив эту функцию из первого равенства (19), получим замкнутое относительно

/ (г)

щ трансцендентное уравнение, решение которого имеет вид

ф^ = ±aT<^{-yJ^+eDJ$^‘),uls)l} +(3, (ж1,ж2) €ГР. (29)

Здесь угол (3 задаёт направление внешней нормали к контуру Гр, поэтому, чтобы ортогональные характеристики системы (15)—(17) входили в область О, занимаемую пластиной в плане, необходимо решение уравнения (29) разыскивать на открытом интервале (ф^ — [3) € (тг/2,7г)и(7г, 37г/2). Это же требование приводит к выбору знака (—) перед радикалом в (29). Если краевые значения функции -ф^ (Гр) уже известны из уравнения (29), то второе граничное условие (19) определяет контурные значения функции

+ eDpiv^, ш^) - 0,5а?9г [wqi sin2(V’ir') - f3)\ ), (хі,х2) ЄГР. (ЗО)

4Г)(ГР) = “——Ай---------~{Fnz+eCF(v^,uj^) +

<7° Sin(-0i — (3)

Равенства (29), (30) имеют смысл при условиях sin^^ — (3) ф 0,

cos(,0^ — f3) ф 0, Woi Ф 0, причём два первых из этих неравенств возможны лишь при условии

0 < [Mn + eDm(v^ , w*^)] /(cr^Woi) < 1, (xi,x2) еГр. (31)

При выполнении ограничений (31) равенство (29) за счёт выбора знаков

(±) определяет два корня г/^(Гр), которые попадают в указанный выше интервал и задают два направления, симметричные относительно нормали к контуру Гр. При известной контурной функции г/^(Гр) равенство (30) однозначно определяет А^(Гр). Следовательно, статические граничные условия

(19) и краевое условие (21) задают на контуре Гр два набора альтернативных краевых условий для подсистемы (15)—(17). А значит, в классе аналитических входных данных задача Коши для этой подсистемы может быть корректно сформулирована.

Косвенным доказательством этого может служить следующий частный случай — случай чистого изгиба пластины, когда р = 0, Fnz = 0. При Мп = = const ф 0, Woi = const ф 0 из (15)—(17), (29), (30) в силу теоремы Коши-

Ковалевской следует, что на первой итерации = 0, а подсистема (16), (17) формально совпадаёт с (17), (31) в [19], т. е. траектории РА в первом при-

(V) (г)

ближении являются прямыми линиями, а функции ф\ , ш\ определяются аналитическими выражениями (3.14)—(3.24) в [20].

Замечание 1. Как уже отмечалось, в общем случае нагружения пластины (р ф 0, Fnz ф 0, Мп ф const) теоремы существования решения задачи Коши (15)—(17), (21), (29), (30) в более широких, чем аналитический, классах входных данных автору не известны. Однако в рассмотренном выше случае чистого изгиба решение этой задачи Коши на первой итерации существует, если краевые значения для функций ф^ принадлежат классу гладких

функций, а краевые значения для (т. е. ш^к) принадлежат классу непрерывных функций, заданных на части контура Гр. (Об этом свидетельствует структура решения (3.19), (3.21)-(3.24) в [20].)

Если после интегрирования задачи Коши (15)—(17), (21), (29), (30) её область определённости покрывает всю область G, занимаемую пластиной

в плане, то в каждой точке {х\,х2) € G известны функции \ А^\ ш[ \ В этом случае подсистема (18), согласно (28), является линейной системой уравнений параболического типа [16]. Для этой подсистемы также естественно ставится задача Коши, причём в качестве краевых условий для функций v(r\ Vвыступают кинематические граничные условия (20), в которых функции г^Гг,), 9n(Tv) предполагаются аналитическими. Напомним, что по высказанному ранее условию волокна «выходят» из пластины на части контура Г„, где заданы кинематические граничные условия (20). В связи с этим для однозначной формулировки задачи Коши для системы (18) необходимо, чтобы характеристики системы (15)—(17) (см. (27)) нигде не касались контура IV

Таким образом, на каждой итерации решение задачи РА сводится к последовательному интегрированию задачи Коши (15)—(17), (21), (29), (30) с учётом (31) и задачи Коши (18), (20). Такое неклассическое сведение краевой задачи РА к последовательности решений двух задач Коши возникает из-за того, что система разрешающих уравнений (1)—(4) совместно с граничными и краевыми условиями (6), (8), (9) образует краевую задачу с сингулярным возмущением [6], а итерационный процесс (15)—(22) относится к разряду методов теории возмущений.

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

§(г) = _е[Б(^(г),и;(г)) -£(«(*>,о/*))], (хъх2) € G,

5м = -г[-°м(г>(г),и>(г)) - DM

5Р = -e[CF{v^\^r)) - CF{v^s\^s)) + DF{v^r\^r))-

- Df(v^s\u}^)], (xi,x2)eTp, r = s +1,

где нужно учесть соотношения (5), (7).

Равенства (32) показывают, что определение невязок при использовании итерационного процесса не требует дополнительных вычислительных затрат, так как все выражения, входящие в (32), необходимо вычисляются для построения решения на (г + 1)-й итерации. Автор считает, что на настоящий момент времени контроль норм невязок (32) является наиболее удобным практическим критерием проверки сходимости процесса (15)—(22).

Замечание 2. Так как равенства (21), (29), (30) с учётом (31) могут определять на контуре Гр до двух альтернативных наборов краевых условий для системы (15)—(17) (выбор того или иного набора этих условий определяется только волей проектировщика или из каких-либо дополнительных соображений), удовлетворяющей условиям теоремы Коши—Ковалевской [18], а также в силу единственности решения задачи Коши (18), (20) (в случае аналитических входных данных) при помощи итерационного метода (15)—(22) можно построить до двух различных решений исходной задачи РА, если, конечно, они существуют. Этими решениями можно управлять за счёт варьирования

контурных функций шок в краевых условиях (21). (Так как выбор того или иного решения всецело зависит от воли проектировщика, то в процессе решения задачи РА при использовании разработанного итерационного метода не может произойти «несанкционированного» перехода от одного решения к другому, например за счёт погрешности численных схем и т. п.) Такое многообразие альтернативных проектов РА выгодно с практической точки зрения, так как позволяет выбрать из их совокупности проекты, обладающие теми или иными эксплуатационными свойствами, а также наиболее удобные с точки зрения технологической реализации.

Замечание 3. Два корня ф^[\ определяемые равенством (29), задают два направления, симметричных относительно нормали к контуру Гр. Поэтому анализ уравнений итерационного процесса показывает, что если форма пластины, её нагружение и закрепление обладают симметрией относительно некоторой оси, лежащей в плоскости (ж1,ж2), то два получающихся решения задачи РА будут симметричными относительно этой оси. С технологической точки зрения эти два решения представляют собой один и тот же проект РА. Действительно, если уже изготовлен один из указанных проектов, то второй можно получить поворотом вокруг оси симметрии на угол 7Г. (Эти результаты подтверждаются конкретными расчётами, например в осесимметричном случае [2].)

Замечание 4. Как уже отмечалось в начале настоящего пункта, решение рассматриваемой задачи РА можно построить и с помощью разложения всех неизвестных функций в асимптотические ряды по малому параметру [6]:

СЮ СЮ

Хт){„ „ \ т—1

Фк(хі,ж2) = ^Фк (хі,Х2)єт шк{жі,ж2) = ^2^к (жі,ж2)є”

т=1 т=1

ос ос

Ак{жі,ж2) = ^ А('™\х1,х2)єт~1, у(х1,х2) = ^2 У(-т')(х1,х2)єт~1, (33)

т=1 т=1

ос ос

\('!!П\х1,х2)ет~1, ,(Ж1,Ж2)=

т= 1 т=1

ос

14(жьж2) = ^2 У^т\х1,х2)ет~1, Ук = д3(фк,ь), к = 1,2,...,Ж.

т= 1

Асимптотический анализ системы разрешающих уравнений (1)—(4) и соответствующих ей граничных (6), (8) и краевых (9) условий с учётом разложений (33) показывает следующее:

1) первые слагаемые {ф^, А^\ У^) в рядах (33) полностью

совпадают с решением системы (15)—(22) на первой (г = 1) итерации;

2) если на второй (г = 2) итерации метода (15)—(22) левые части уравнений

(15)—(19), (21) и правую часть второго уравнения (20) линеаризовать в окрестности решения для первой итерации, то получающаяся система и краевые условия будут совпадать с уравнениями и краевыми условиями, определяющими коэффициенты при вторых слагаемых (ф^, ,

у<у2\ У^) в разложениях (33).

На основе столь существенного сходства уравнений итерационного процесса (15)—(22) и метода асимптотических рядов (33) все результаты исследований, полученные ранее в настоящем пункте при анализе уравнений (15)—(22),

можно перенести и на соответствующие цепочки уравнений, получающие-ся после подстановки (33) в исходную систему разрешающих уравнений и краевые условия задачи РА (п. 1). Однако итерационный процесс (15)—(22) обладает рядом преимуществ перед методом асимптотических рядов:

1) правые части в уравнениях (15)—(20) на каждой итерации имеют одну и ту же достаточно простую структуру (в цепочках же уравнений метода асимптотических рядов (33) структура правых частей с увеличением номера уточнения т существенно усложняется);

2) для построения г-того приближения по методу итераций используется информация лишь о предыдущем (г — 1)-м приближении (в правых же частях уравнений метода асимптотических рядов используется информация о всех предыдущих п-х уточнениях: п = 1, 2,... , г — 1);

3) характеристики подсистем (15)—(17), (18) изменяются от итерации к итерации и при сходящемся итерационном процессе двукратные характеристики совпадают в пределе с траекториями армирования (в уравнениях же метода асимптотических рядов (33) характеристики в соответствующих подсистемах уравнений фиксированы и совпадают с характеристиками подсистем (15)—(17), (18) на первой итерации; см. (27), (28) при г = 1);

4) для проверки сходимости обоих обсуждаемых методов целесообразно г-тые приближения по методу итераций и г-тые частичные суммы типа (33) подставить в систему разрешающих уравнений и соответствующие ей краевые условия (п. 1) и проконтролировать нормы возникающих при этом невязок, причём в случае метода асимптотических рядов (33) эти невязки возникают во всех уравнениях и краевых условиях исходной системы (п. 1) и с увеличением номера уточнения (с увеличением г) структура их существенно усложняется, в случае же итерационного процесса невязки возникают только в уравнениях (1), (6) и имеют достаточно простую структуру (32).

Указанные преимущества позволяют отдать предпочтение итерационному процессу (15)—(22) при практическом решении задач РА изгибаемых металлокомпозитных пластин, работающих в условиях установившейся ползучести.

Замечание 5. Если задача РА уже решена (известны функции у),

то по формулам (3)—(5) из [21] можно определить напряжённое состояние в связующем и в композиции (сгу) пластины, изгибаемой в условиях установившейся ползучести. Кроме того, зная ф\, ш\, можно определить и начальное напряжённое состояние в пластине (его можно вычислить, например, методами, изложенными в [22,23]). Решение же о неустановившейся ползучести на первой её стадии можно получить приближенно, используя методику, подробно описанную в замечании 2 в [21].

Особенности уравнений итерационного процесса (15)—(22), которые обсуждались в настоящем пункте, требуют привлечения высокоточных устойчивых методов численного интегрирования задач Коши (15)—(17), (21), (29), (30) и (18), (20), поэтому для интегрирования этих задач использовались обобщённые методы Рунге—Кутты четвёртого порядка точности [24]. (Более подробное обсуждение этого вопроса выходит за рамки настоящего исследования.)

Решение задачи РА должно удовлетворять физическим ограничениям (10). Так как на каждой итерации уравнение (17) формально совпадает с (2), то, как показано в [25], при <^01^1) > 0 (см. (9), (21)) всюду в области определённости решения задачи Коши для подсистемы (15)—(17) будет выполняться неравенство и)^\х1,х2) > 0. Следовательно, при 6^01 (Г 1) > 0 (см. (29)) на каждой итерации, а значит, и в пределе, если итерационный процесс сходится, первые N неравенств (10) выполняются автоматически, поэтому в процессе решения задачи РА необходимо контролировать выполнение лишь последнего неравенства (10), нарушение которого означает, что при технологической укладке волокна выпучиваются из элементарных армированных слоев из-за недостатка места.

3. Обсуждение результатов расчётов. В качестве примера рассмотрим расчёт РА-структур при установившейся ползучести изгибаемой металлокомпозитной двусвязной пластины, ограниченной контурами, уравнения которых заданы в безразмерной полярной (г, ф) системе координат:

Гр = Г1 : Г\{(р) = 0,5 + 0,052 сов(3^ — 37г/2),

Г^: Г2 (<р) = 0,9 — 0,15 сое (р, 0^(р<2тт.

По внешнему контуру Г^ пластина жёстко защемлена (г?о(^)=0, вп(<р)=0, см. (8)), а на внутреннем контуре Гр задан равномерно распределённый безразмерный изгибающий момент Мп = 0,55, Рпх = 0 (см. (6)). Безразмерная равномерно распределённая поперечная нагрузка характеризуется значением р = 0,08 (см. (1)). Пластина изготовлена из меди и армирована одним семейством стальной проволоки У8А (механические характеристики установившейся ползучести фаз композиции имеют значения (13), (14)). Напряжения в арматуре на лицевых поверхностях по модулю имеют значения | а® | = = 1,2 ГПа. Согласно поведению кривой 3 на рис. 1, а и линии 1 на рис. 1, б, малый параметр е = 2,28 • 10_3.

Плотность армирования проволокой У8А на внутреннем контуре имеет значение

ал(Гр)=а;о1Ы = 0,7, Гр = Гь (34)

На рис. 2 изображены две альтернативные структуры, соответствующие двум разным решениям рассматриваемой задачи РА. Проект, изображенный на рис. 2, а, определяется выбором знака (—) в (29), а на рис. 2, б — знаком (+). В этих двух проектах существенно различаются не только РА-структуры, но и податливости пластин и напряжённо-деформированное состояние в связующем материале. Общий расход арматуры в проектах, изображенных на рис. 2, примерно одинаков, однако при схеме укладки волокон, приведённой на рис. 2, а, максимальная по модулю скорость прогиба установившейся ползучести на 9 % больше, чем при укладке, приведённой на рис. 2, б. Более того, максимальное значение интенсивности скоростей деформаций ползучести в связующем пластины, изображенной на рис. 2, а, на порядок больше, чем в пластине на рис. 2, б.

Согласно замечанию 2, РА-проектами можно управлять за счёт варьирования контурных значений плотностей армирования Сс>01 (см. (9)). Продемонстрируем это, задав на внутреннем контуре, в отличие от (34), переменную

плотность армирования

Ш1(ГР) = и)01(<р) = 0,68 + 0,1 сов((р + 0,85), Гр = Г1, 0 ^ <р < 2тт. (35)

На рис. 3 изображены два альтернативных РА-проекта при задании краевого условия (9) в виде (35). Общий расход арматуры в этих проектах также примерно одинаков и на 3,2 % меньше, чем в проектах, изображенных на рис. 2. Однако в случае (35) скорость прогиба пластины, изображенной на рис. 3, а, на 19,7 % больше, чем пластины на рис. 3, б, а интенсивность скоростей деформаций в связующем пластины, приведённой на рис. 3, а, в 7,9 раза больше, чем в пластине на рис. 3, б.

На рис. 2, 3 траектории РА изображены таким образом, что между двумя соседними траекториями на внутреннем контуре Гр = Г1 входит фиксиро-

Рис. 2. Два альтернативных проекта равнонапряжённого армирования двусвязной пластины, армированной на внутреннем контуре с постоянной плотностью (34)

а б

Рис. 3. Два альтернативных проекта равнонапряжённого армирования двусвязной пластины, армированной на внутреннем контуре с переменной плотностью (35)

ванное количество арматуры. Это позволяет получить некоторое наглядное представление о действительном распределении плотностей армирования в РА-проектах, приведённых на этих рисунках. При совмещении рис. 2, а и 3, а, а также рис. 2, б и 3, б можно увидеть, что точки «входа» изображенных траекторий армирования на рис. 3 сместились против часовой стрелки относительно тех же точек и траекторий РА на рис. 2. Кроме того, изменилась и кривизна траекторий как по знаку, так и по модулю. В целом, продемонстрированные проекты качественно отражают все многообразие решений задачи РА изгибаемых пластин, работающих в условиях установившейся ползучести.

Возможность построения двух альтернативных РА-проектов при заданном краевом условии (9) (в частности, (34), (35)) определяется тем, что задачи РА относятся к разряду обратных задач механики композитов (а точнее, к специфическим коэффициентным задачам), которые, как известно, могут иметь неединственное решение.

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

Работа выполнена при поддержке РФФИ (код проекта 11-08-00186-а).

БИБЛИОГРАФИЧЕСКИЙ СПИСОК

1. А. 77. Янковский, “О некоторых свойствах решения задачи равнонапряженного армирования изгибаемых металлокомпозитных пластин, работающих в условиях установившейся ползучести”// Вестн. Сам. гос. техн. ун-та. Сер. Физ.-мат. науки, 2011. №2(23). С. 62-73. [А. P. Yankovskii, “On some properties of equal-stress problem solution reinforcement bending the metal-composite plates working in steady creep conditions” // Vestn. Samar. Cos. Tekhn. Univ. Ser. Fiz.-Mat. Nauki, 2011. no. 2(23). Pp. 62-73].

2. A. 77. Янковский, “Равнонапряженное армирование кольцевых изгибаемых металлокомпозитных пластин, работающих в условиях установившейся ползучести”// Вестн. Сам. гос. техн. ун-та. Сер. Физ.-мат. науки, 2010. №5(21). С. 42-54. [А. P. Yankovskii, “Equal-Stress Reinforcement the Ring Bending Metal-Composites Plates Working in Conditions of Steady Creep” // Vestn. Samar. Cos. Tekhn. Univ. Ser. Fiz.-Mat. Nauki, 2010. no. 5(21). Pp. 42-54].

3. Т. Д. Джураев, Краевые задачи для уравнений смешанного и смешанно-составного типов. Ташкент: Фан, 1979. 239 с. [Т. D. Dzhuraev, Boundary value problems for equations of mixed and mixed-composite types. Tashkent: Fan, 1979. 239 pp.]

4. А. А. Самарский, Теория разностных схем. М.: Наука, 1989. 616 с. [A. A. Samarskii, Theory of Difference Schemes. Moscow: Nauka, 1989. 616 pp.]

5. О. C. Zeinkiewicz, R. L. Taylor, The finite element method. Vol. 1: The Basis. Oxford: Butterworth-Heinemann, 2000. xviii+689 pp.

6. A. Nayfeh, Introduction to Perturbation Techniques. New York: Wiley-Interscience, 1981. xiv+519 pp.; русск. пер.: А. Найфе, Введение в методы возмущений. М.: Мир, 1984. 536 с.

7. Г. С. Писаренко, 77. С. Можаровский, Уравнения и краевые задачи теории пластичности и ползучести: Справочное пособие. Киев: Наукова думка, 1981. 496 с. [С. S. Pisarenko, N. S. Mozharovskii, Equations and Boundary-Value Problems of the Theory of Plasticity and Creep: Handbook. Kiev: Naukova Dumka, 1981. 496 pp.]

8. А. Ф. Никитенко, Ползучесть и длительная прочность металлических материалов. Новосибирск: Новосиб. госуд. архит.-строит. унив., 1997. 278 с. [A. F. Nikitenko, Creep

and Long-Term Strength of Metal Materials. Novosibirsk: Novosib. Gosud. Arhkit. Stroit. Univ., 1997. 278 pp.]

9. Композиционные материалы: Справочник / ред. Д. М. Карпинос. Киев: Наукова думка, 1985. 592 с. [ Composite Materials: Handbook/ ed. D. M. Karpinos. Kiev: Naukova Dumka, 1985. 592 pp.]

10. Д. М. Карпинос, В. А. Невгод, JI. И. Тучинский, Л. H. Клименко, Л. P. Вишняков, В. А. Беспятый, В. X. Кадыров, “Ползучесть и длительная прочность вольфрамовых проволок”// Пробл. прочности, 1972. №1. С. 70-73; англ. пер.: D. М. Karpinos, V. A. Nevgod, L. I. Tuchinskii, L. N. Klimenko, L. R. Vishnyakov, V A. Bespyatyi, V. Kh. Kadyrov, “Creep and durability of tungsten wire” // Strength of Materials, 1972. Vol. 4, no. 1. Pp. 72-75.

11. В. В. Новожилов, К. Ф. Черных, Е. И. Михайловский, Линейная теория тонких оболочек. Л.: Политехника, 1991. 656 с. [V. V. Novozhilov, К. F. Chernykh, Е. I. Mikhaylovskiy, Linear Theory of Thin Shells. Leningrad: Politekhnika, 1991. 656 pp.]

12. E. Kamke, Differentialgleichungen. Losungsmethoden und Losungen. Vol. II: Partielle Differentialgleichungen erster Ordnung fur eine gesuchte Funktion. Leipzig, 1959; русск. пер.: Э. Камке, Справочник по дифференциальным уравнениям в частных производных первого порядка. М.: Наука, 1966. 260 с.

13. А. Джураев, Системы уравнений составного типа. М.: Наука, 1972. 228 с. [A. Dzhuraev, Systems of Equations of Composite Type. Moscow: Nauka, 1972. 228 pp.]

14. Г. М. Фихтенголъц, Курс дифференциального и интегрального исчисления. Т. 1. М.: Наука, 1970. 607 с. [С. М. Fikhtengol’ts, A Course of Differential and Integral Calculus. Vol. 1. Moscow: Nauka, 1970. 607 pp.]

15. Б. Л. Рождественский, H. H. Яненко, Системы квазилинейных уравнений. М.: Наука, 1969. 592 с. [В. L. Rozhdestvenskii, N. N. Yanenko, Systems of Quasi-Linear Equations. Moscow: Nauka, 1969. 592 pp.]

16. О. А. Ладыженская, Линейные и квазилинейные уравнения параболического типа. М.: Наука, 1967. 736 с. [О. A. Ladyzhenskaya, Linear and Quasi-Linear Equations of Parabolic Type. Moscow: Nauka, 1967. 736 pp.]

17. Л. М. Качанов, Основы теории пластичности. М.: Наука, 1969. 420 с. [L. М. Kachanov, Fundamentals of the Theory of Plasticity. Moscow: Nauka, 1969. 420 pp.]

18. И. Г. Петровский, Лекции об уравнениях с частными производными. М.: Физматгиз, 1961. 400 с. [I. С. Petrovskii, Lectures on Partial Differential Equations. Moscow: Fizmatgiz, 1961. 400 pp.]

19. А. П. Янковский, “Применение методов теории возмущений в плоской задаче равнонапряжённого армирования металлокомпозитных пластин при установившейся ползучести”// Вестн. Сам. гос. техн. ун-та. Сер. Физ.-мат. науки, 2009. №2(19). С. 53-71. [A. P. Yankovskii, “Application of methods of the theory of perturbations in flat problem of equally-stressed reinforcing of metal-composite plates at steady creep” // Vestn. Samar. Cos. Tekhn. Univ. Ser. Fiz.-Mat. Nauki, 2009. no. 2(19). Pp. 53-71].

20. Ю. В. Немировский, А. П. Янковский, “Равнонапряжённое армирование металлокомпозитных пластин волокнами постоянного поперечного сечения в условиях установившейся ползучести”// Мех. композ. матер., 2008. Т. 44, №1. С. 11-34; англ. пер.: Yu. V. Nemirovskii, А. P. Yankovskii, “Equal-stressed reinforcement of metal-composite plates with fibers of constant cross section in steady-state creep” // Mech. Compos. Mater., 2008. Vol. 44, no. 1. Pp. 9-24.

21. Ю. В. Немировский, А. П. Янковский, “Установившаяся ползучесть слоисто-волокни-стых изгибаемых металлокомпозитных пластин” // Вестн. Сам. гос. техн. ун-та. Сер. Физ.-мат. науки, 2008. №2(17). С. 66-76. [Yu. V. Nemirovskii, А. P. Yankovskii, “The steady creep layer-fibrous bendings metal-composites plates” // Vestn. Samar. Cos. Tekhn. Univ. Ser. Fiz.-Mat. Nauki, 2008. no. 2(17). Pp. 66-76].

22. Ю. В. Немировский, Б. С. Резников, Прочность элементов конструкций из композитных материалов. Новосибирск: Наука, 1986. 168 с. [Yu. V. Nemirovskii, В. S. Reznikov,

Strength of Structural Elements Made of Composite Materials. Novosibirsk: Nauka, 1986. 168 pp.]

23. Ю. В. Немировский, А. П. Янковский, “Термоупругопластический изгиб сложно армированных пластин”// Мех. композ. матер., 2005. Т. 41, №6. С. 715-742; англ. пер.: Yu. V. Nemirovskii, А. P. Yankovskii, “Thermoelastoplastic Bending of Complexly Reinforced Plates” // Mech. Compos. Mater., 2005. Vol. 41, no. 6. Pp. 477-496.

24. Ю. В. Немировский, А. П. Янковский, “Обобщение методов Рунге-Кутты и их применение к интегрированию начально-краевых задач математической физики” // Сиб. журн. вычисл. матем., 2005. Т. 8, №1. С. 57-76. [Yu. V. Nemirovskii, А. P. Yankovskii, “Generalization of the Runge-Kutta methods and their application tointegration of initial-boundary value problems of mathematical physics” // Sib. Zh. Vychisl. Mat., 2005. Vol. 8, no. 1. Pp. 57-76].

25. Ю. В. Немировский, А. П. Янковский, “О некоторых особенностях уравнений оболочек, армированных волокнами постоянного поперечного сечения” // Механика композиционных материалов и конструкций, 1997. Т. 3, №2. С. 15-39. [Yu. V. Nemirovskii, А. P. Yankovskii, “On some features of the equations of shells reinforced with fibres of constant crosssection” // Mekhanika Kompozitsionnykh Materialov і Konstrukciy, 1997. Vol. 3, no. 2. Pp. 15-39].

Поступила в редакцию 31/1/2013; в окончательном варианте — 27/11/2013.

MSC: 74K20, 74R20

APPLICATION OF METHODS OF THE PERTURBATION THEORY TO PROBLEM OF EQUALLY-STRESSED REINFOCING OF BENDING METAL-COMPOSITE PLATES IN CONDITIONS OF STEADY-STATE CREEP

A. P. Yankovskii

Khristianovich Institute of theoretical and applied mechanics the Siberian Branch of the Russian Academy of Science,

4/1, Institutskaya St., Novosibirsk, 630090, Russia.

E-mails: lab4nemir@rambler.ru

The problem of equally-stressed reinforcing of bending metal-composite plates in conditions of steady-state creep, is solved by the help of methods of the perturbation theory. The opportunity of existence of several alternative solutions of a considered problem which can be reliably determined, using the developed, algorithm is shown. Concrete projects are constructed for equally-stressed reinforcing doubly connected, plates with different density of reinforcing on an internal contour.

Key words: bending plates, steady-state creep, rational designing, equally-stressed reinforcing, small parameter, perturbation theory.

Original article submitted 31/1/2013; revision submitted 27/11/2013.

Andrey P. Yankovskii (Dr. Sci. (Phys. & Math.)), Leading Research Scientist, Lab. of Fast Processes Physic.