CC BY
41
УДК 539.376
ПРИМЕНЕНИЕ МЕТОДОВ ТЕОРИИ ВОЗМУЩЕНИЙ В ПЛОСКОЙ ЗАДАЧЕ РАВНОНАПРЯЖЁННОГО АРМИРОВАНИЯ МЕТАЛЛОКОМПОЗИТНЫХ ПЛАСТИН ПРИ УСТАНОВИВШЕЙСЯ ПОЛЗУЧЕСТИ
А. П. Янковский
Институт теоретической и прикладной механики им. С. А. Христиановича СО РАН,
630090, г. Новосибирск, ул. Институтская, 4/1.
E-mails: nemirov@itam.nsc.ru, shulgin@itam.nsc.ru
Методами теории возмущений решена задача равнонапряжённого армирования пластин, нагруженных в своей плоскости и 'работающих в условиях установившейся ползучести. Показана возможность существования нескольких альтернативных решений рассматриваемой задачи, которые можно надежно выделять, используя разработанный алгоритм. Построены конкретные проекты равнонапряжённого армирования двусвязной пластины при разных типах её нагружения.
Ключевые слова: пластины, установившаяся ползучесть, рациональное проектирование, равнонапряжённое армирование, малый параметр, теория возмущений.
Введение. При проектировании металлокомпозитных пластин целесообразно осуществлять их оптимизацию. Одним из прочностных критериев рационального армирования тонкостенных конструкций является требование равнонапряженности силовых элементов — волокон (проволок) — вдоль их траекторий, так как при этом несущая способность арматуры используется наиболее полно.
При длительной эксплуатации изделия в случае статического термосилового нагружения подавляющую часть времени металлокомпозитная конструкция работает в условиях установившейся ползучести [1], поэтому актуальной является проблема равнонапряженного армирования (РА) тонкостенных конструкций, работающих в условиях установившейся ползучести всех фаз композиции.
Постановка задачи РА металлокомпозитных пластин в условиях установившейся ползучести, нагруженных в своей плоскости и усиленных волокнами постоянного поперечного сечения, и математический анализ соответствующей системы разрешающих уравнений приведены в [2], где показано, что получающаяся система квазилинейных уравнений смешанно-составного типа и соответствующие ей нелинейные граничные условия не позволяют надеяться на построение широкого круга аналитических решений и требуют разработки обоснованных численных процедур решения. В то же время отмеченные в [2] особенности разрешающей системы уравнений не позволяют произвести простой перенос известных методов численного интегрирования квазилинейных систем уравнений [3] на решение рассматриваемой задачи. Наличие в задаче естественного малого параметра (см. ниже) позволяет надеяться на
Андрей Петрович Янковский (д.ф.-м.н.), ведущий научный сотрудник, лаборатория физики быстропротекающих процессов.
эффективность применения хорошо зарекомендовавших себя методов теории возмущений [4 и др.].
Настоящая работа посвящается разработке и исследованию специфики применения методов теории возмущений для построения и анализа особенностей решений плоской задачи РА металлокомпозитных пластин в условиях установившейся ползучести.
1. Система разрешающих уравнений и граничные условия. Для удобства дальнейшего изложения систему разрешающих уравнений задачи РА и соответствующие ей граничные условия [2] будем записывать в безразмерном виде, предполагая, что в пластину постоянной толщины h внедрено N семейств армирующих проволок. В дальнейшем все размерные функции и величины будем помечать сверху чертой, а соответствующие им безразмерные функции и величины — обозначать теми же символами, но без черты. Тогда замкнутая система разрешающих уравнений, содержащая в качестве неизвестных функций безразмерные скорости vi, V2 перемещений при установившейся ползучести, безразмерную температуру в и параметры РА Шк, фк, состоит из [2]: уравнений равновесия
(-1)*^ ОкШкlkjds (фк,фк) + eAi (v, ш, в) = -Xi (x, w), j = 3 - i, i = 1, 2; (1)
к
условий постоянства поперечных сечений волокон
ds (фк ,Шк)+ Шкдп (фк ,фк) = 0, k = 1, 2,..., N; (2)
уравнения стационарной плоской теплопроводности
(Л11в,1 + Л12в,2) д + (Л21в, 1 + Л22в,2),2 + № (в^> — в) = — Q (x, ш) ; (3)
условий равнонапряженности арматуры
ds (фк, V1)cos фк + ds (фк,V2) sin фк = Ск(в) = /к (о к,в), (4)
ак = const, k = 1, 2,... , N,
где
A* (v, w, в) = 2 [ago (H, в) (2vi,i + Vj,j)] ,* + [ago (H, в) (v*,j + Vj,*)] ,j ,
j = 3 — i, Xi, (x, w) = Xo (x) a + ^ Х(к) (x) шк, i = 1,2, a = 1 — Q,
к
Q = ^ Шк, Q (x, w) = aQo (x) + £ Шк Qk (x), 1к1 = cos фк, 1к2 = sin фк,
кк
ds (фк, •) = (-)д cos фк + 0,2 sin фк, dn (фк, •) = ds (фк + п/2, •) =
= - Од sin фк + 0,2 cos фк, k = 1,2,..., N, Лij = Q-1^ Шк x (5)
к
X |[Q (Ак - Ao) + Ao] ^i^j + (-1)i+jlкslкrAкAo [Q (Ao - Ак) + Ак] 1 j ,
s = 3 - i, r = 3 - j, i,j = 1,2, v = {v1,v2} , x = {ж1,ж2} ,
Ш = {ui, Ш2, ■ ■ ■ , COn} , я = 2^/^д + ^1,1^2,2 + ^2,2 + °>25 0l,2 + ^2,l)2;
суммирование производится по указанному индексу от 1 до N; нижний индекс после запятой означает частное дифференцирование по соответствующей переменной Хі.
Пусть область С, занимаемая конструкцией в плане, ограничена контуром Г, тогда на одной части (обозначим её Гр) этого контура могут быть заданы статические [2]
У, CkWk cos2 ^k - в) + eDn (v, и, в) = pn
Л (б)
> усkWk sin 2 (фk - в) + £Dt (v, и, в) = 2рт, (Xl,X2) Є Гр, k
где
Dn (v, и, в) = 2ag0 (H, в) [Vl,l (1 + n^ + V2,2 (1 + n2) + (Vl,2 + V2,l) П1П2] ,
Dt (v, и, в) = 2ag0 (H, в) [2 (V2,2 - Vl,l) П1П2 + (Vl,2 + V2,l) (nl - n2)] , (7)
nl = cos в, n2 = sin в,
на другой части (обозначим её Г*) контура Г — кинематические
Vi (Г*)= Vе (Г), i = 1, 2 (В)
и на всем контуре Г (= Гр U Г*) —тепловые
X [(Ліі в,1 + Л12 в,2) П1 + (Л21 в,1 + Л22 в,2) П2 + q] + 7 (в - ве) = 0 (9)
граничные условия. Кроме того, на части контура Г (обозначим её Г^), на которой волокна k-того семейства «входят» в пластину, необходимо задать краевые условия для интенсивностей армирования [2]:
Wk (Г^)= Wok, k = 1, 2,... , N. (10)
Решение краевой задачи РА (1)-(4), (б), (В)-(10) должно удовлетворять физическим ограничениям [2]
0 < Wk (k = 1, 2,..., N), 0 < Q < w* < 1. (11)
В уравнениях и соотношениях (1)-(11) использованы следующие формулы обезразмеривания и обозначения:
Ck = Ck/СЪ Ck = Ck/СЪ ^ ^/Л^
fk (Ck, в) = fk (Л*Ck, Л*в) /Л* = fk fa, в) /Л*,
go (H, в) = go (№ Л*в) /B* = go (H, л) /B*,
H = H/С*, Vi = Vi/(DC*), Xi = Xi/D ,
Xi = DXi/л *, x° = DXC/Л *, Xf) = DX^/л *, (12)
Ле = Ло/Л *, в = л/л *, ве = ле/л*, вте = лте/л*,
Q = D2Q/^ *л*), Qk = D2Qk/(л*в*), Qe = D2Qe/fa л*),
q = Dg/^ *л*), jU = 2D^/(hA *), Pn = pn/C *, Pt = Pt / C *,
Vo = v0/(DC*), Л, = Лі,-/л*, є = B*/B*, i, j = 1, 2, k = 1, 2,... , N;
Хг, у — компоненты векторов приведённых объёмных (массовых) нагрузок и скорости установившейся ползучести точек пластины по направлениям Щ (г = 1, 2) прямоугольной декартовой системы координат; В — характерный
размер пластины в плане; X0, Х(к) —компоненты удельных объёмных (массовых) нагрузок, действующих на связующее и арматуру к-того семейства; Щ, Ск — напряжение и скорость деформации ползучести в к-том семействе арматуры, связанные известной зависимостью Ск = /к (Ск, Щ; Щ (Й,0) — заданная функция, являющаяся коэффициентом пропорциональности между интенсивностью касательных напряжений Т и интенсивностью скоростей деформаций Й в связующем: Т = $о (Й, 0) Й; Шк, фк - интенсивность и угол (отсчитываемый от направления Ж1) армирования волокнами к-того семейства; ш* — предельно допустимая удельная суммарная плотность армирования (на практике ш* ^ 0,7); 0 , 0^, — температура конструкции и температура окружающей среды на лицевых плоскостях пластины; Д — коэффициент конвективного теплообмена между связующим и окружающей средой на лицевых поверхностях пластины; Яо, Як — плотности мощностей внутренних источников тепла в связующем и волокнах к-того семейства; Я — приведённая плотность мощности внутренних источников тепла в композиции; Що, 0к — коэффициенты теплопроводности связующего и арматуры к-того семейства (в случае учёта термочувствительности фазовых материалов зависят от 0 ); рп, рт — нормальные и касательные контурные напряжения соответственно; Щ — заданные компоненты скорости перемещения точек контура Г; 0о — заданная на контуре Г температура; 0 — тепловой поток через торцевую (боковую) поверхность конструкции; %, 7 — функции переключения, позволяющие задавать на Г разные тепловые граничные условия; в — угол, задающий направление внешней нормали к Г; Шок — заданные на контуре Г^к краевые значения функций Шк; 0* — характерное значение температуры конструкции; А* — характерное значение коэффициента теплопроводности связующего; а *, С * — характерные значения напряжения и скорости деформации ползучести в первом семействе арматуры, связанные линейной зависимостью а* = В *С *; В * — характерное значение коэффициента пропорциональности $о (Й, 0) для связующего.
В случае степенного закона С = В0т, связывающего скорость деформации С ползучести и напряжение а в материале каждой фазы композиции, имеем зависимости [1]
4о (Й ,0) = Ьо (0) ЙК<0-1,
Ъо = (Уз^+^ВоГ^0, д0 (0) = 1 /т0 (0) , (13)
/к (Ск, Щ) = Вк (0) акт"^, к = 1, 2,... , Ж,
где Во (0, то (0), Вк (0), Шк (0) — известные характеристики установившейся ползучести материалов связующего и арматуры к-того семейства, зависящие заданным образом от температуры 0 .
При выборе определенных материалов фаз композиции параметр е (см. (12)) может быть малым. Действительно, пусть металлокомпозитная пластина получена путём внедрения в медное (Си) связующее стальной проволоки
У8А, тогда при температуре 9 ~ 200 °С характеристики материалов фаз композиции (см. (13)) имеют значения [5, 6]:
т-1
Си: то = 2,16, Во = 5,63 ■ 10-9 (МПа)-т0 ■ ч
У8А: ті = 24,98, В1 = 1,054 ■ 10-84 (МПа)-т1 ■ ч
-1
(14)
Аппроксимируем степенные законы установившейся ползучести рассматриваемых материалов с параметрами (14) двухзвенными кусочно-линейными функциями
ак =
*(1)
а
Вкч&, ?(2)/
+ Вк>(£к - С^), Ск >СГ > 0, к = 0,1.
0 < Ск < Ск (і)
(і)
(15)
На рис. 1 кривые 1, 2 характеризуют степенные зависимости а к ~ Ск (к = 0,1) с параметрами (14) для меди и проволоки У8А соответственно. Двухзвенные ломаные 1;, 2; (штриховые линии) аппроксимируют кривые 1, 2 и определяются параметрами (см. (15)):
Си:
В01) = 2,647 ■ 1012 Па ■ ч,
а
(1)
С01 = 1,852 г(1)
10-5 ч-1,
Во(2) = 5,369
У8А: В( ) = 1,611 ■ 1015 Па ■ ч, а
(1)
С 1(1) = 9,097
10-7 ч-1,
В 1(2) = 1,257
= 49,005 МПа,
1011 Па = 1,466 ГПа,
1012 Па
ч;
ч.
(16)
(Параметры (16) определены методом наименьших квадратов на интервале
о < С < 10 4 ч 1, при этом варьировались величины сг , £(1). При принятом на рис. 1 масштабе по оси ординат линии 1, 1 визуально почти не отличаются.)
Если в (12) с учётом (16) принять а * = а (1), С і = С Г;, В* = В0
(1)
?(1)
(1)
В* =
= , то в рассматриваемой задаче е = В*/В£ = 0,00164 — весьма малый
параметр. (Подобные оценки малости параметра е можно получить и в случаях использования в качестве арматуры вольфрамовых проволок [7], а в качестве связующего — алюминиевых сплавов [8] и др.)
а, МПа 1200 Н
800
400
0
1,1'
0 2 4 6 8 £ • 105, 1/ч
Рис. 1. Кривые установившейся ползучести фазовых материалов
В [2] показано, что система разрешающих уравнений (1)—(4) является системой квазилинейных уравнений смешанно-составного типа [9], причём траектории армирования, определяемые углами фк, являются её действительными характеристиками. Структура уравнений (1)—(4) и граничных условий (6) (с учётом соотношений (5), (7)) показывает, что задачи определения параметров РА фк, шк, температурного поля 9 и поля скоростей VI, г2 установившейся ползучести связаны, и решать их необходимо совместно. Последнее обстоятельство существенно усложняет разработку методов решения нелинейной краевой задачи РА (1)-(4), (6), (8)-(10).
2. Применение методов теории возмущений. В настоящее время теория систем квазилинейных уравнений смешанно-составного типа разработана недостаточно полно [9], чтобы можно было аналитически исследовать свойства решений системы разрешающих уравнений (1)-(4) в общем случае, поэтому целесообразно воспользоваться какими-либо упрощениями этой системы и соответствующих ей граничных условий.
Важная особенность системы разрешающих уравнений задачи РА состоит в наличии в уравнениях равновесия (1) и в статических граничных условиях (6) малого параметра е. Это позволяет разыскивать решение плоской задачи РА пластин при помощи методов теории возмущений [4 и др.], например, метода разложения решения в виде асимптотических рядов по степеням малого параметра е [4] или следующего итерационного процесса (который качественно схож с итерационным методом определения основного момент-
ного напряженного состояния в теории оболочек [10]). Пусть ф^, ^к^, г(в), 9(в) (г = 1, 2, 1 ^ к ^ N) — з-тые приближения неизвестных функций, тогда (з + 1)-е приближения будем получать путём интегрирования уравнений
j = 3 — i, i = 1, 2;
(17)
— 0(r)) = -Q (*, ^(r)) ; (18)
ak = const, k = 1, 2,...,N, (ж^ж2) € G
при граничных и краевых условиях
E^kW^ cos2 (ф^ - ft) = Pn - (v(s), <^(s),0(s)) ,
k
(20)
^>kwkr) sin2 (ф^ - = 2Pt - eDr (v(s), <^(s),0(s)) , (Ж1,Ж2) € Гр;
k
(Г)= (Г), i = 1, 2;
(21)
Ц(лЦ^ + AfcW?) ni + (а21)+ Л^М) П2 + q] +
+ y (Vr) - 0o) = 0; (22)
^kr) (r^k )= w0k, k = 1, 2, ...,N, r = s + 1. (23)
Для начала итерационного процесса нулевое приближение выберем в виде
vj0) = 0 (i = 1, 2), 0(о) = const w 1 (например, 0(о) = 0*/0* = 1). (24)
Здесь v(s), w(s) —векторные функции, подобные v, w (см. (5)); ikl =
= cos ф^\ 1fc2) = sin ^kr); выражения для коэффициентов Л(Г) получаются
из Лу (см. (5)) заменой ^, wk их r-тыми приближениями ф^, w(r).
Уравнения итерационного процесса (17)—(23) получаются из соответствующих уравнений разрешающей системы следующим образом. В уравнениях равновесия (1) и в статических граничных условиях (6) в операторы Aj, Dn, Dr, имеющие в качестве сомножителя малый параметр е, подставляются приближения неизвестных функций, определенные на предыдущей итерации. Поэтому слагаемые, имеющие сомножитель е, на каждой итерации известны и в (17), (20) перенесены в правую часть. Задание начального приближения в виде (24) аналогично отбрасыванию в (17), (20) на первой итерации (г = 1) слагаемых, содержащих е. Остальные уравнения итерационного процесса на каждой итерации формально совпадают с соответствующими соотношениями разрешающей системы уравнений (см. п. 1).
Проанализируем систему (17)—(19) при различных количествах семейств арматуры (N). Пусть в конструкцию внедрено одно семейство волокон (N=1). Тогда три квазилинейных уравнения (17) содержат две неизвестные функции
ф(г), w(r), следовательно, подсистема (17) переопределена и для существования её решения необходимо выполнение некоторых условий совместности, которые свяжут между собой внешние нагрузки Xi, X2 [11]. С другой стороны, если система (17) совместна и её решение известно, то уравнение (18) замкнуто относительно приближения температуры 0(r), но при ф(г), 0(r), известных из (17), (18), получаем одно линейное уравнение (19) относительно двух неиз-
(r) (r)
вестных функций v1 , v2 . Таким образом, подсистема (19) при N = 1 оказывается незамкнутой относительно приближений скоростей установившейся ползучести, а значит, при использовании одного семейства волокон решение задачи РА методом (17)—(24) получить нельзя.
Если в конструкцию внедрены два семейства арматуры (N = 2), то для каждого r ^ 1 система (17)—(19) распадается на три подсистемы: первая подсистема (17) состоит из четырех квазилинейных уравнений и замкнута относительно приближений параметров РА ф^, ^k^ вторая подсистема (18) состоит из одного уравнения и при ф^, wkr), известных из (17), замкнута относительно приближения температуры 0(r); третья подсистема (19) состоит из двух уравнений и при ф^, 0(r), известных из (17), (18), замкнута относительно приближений скоростей vkr) (к = 1, 2). Аналогично распадаются граничные и краевые условия (20)—(23): условия (20), (23) состоят из четырёх
равенств и замкнуты относительно контурных значений функций фк , шк ;
(г) (г)
тепловое условие (22) при фк , шк , известных из решения краевой задачи (17), (20), (23), задаёт на Г граничные значения функции 0(г) или условия по тепловому потоку, или граничные условия III рода, а кинематические условия
(21) определяют на Г контурные значения двух функций (к = 1, 2).
Таким образом, при N = 2 итерационный процесс (17)-(24) позволяет расщепить ранее связанные задачи определения параметров РА, температурного поля и поля скоростей установившейся ползучести на отдельные подзадачи, которые можно последовательно интегрировать на каждой итерации.
При внедрении в пластину более двух семейств арматуры ^ ^ 3) система (17)-(19) в целом замкнута относительно неизвестных функций, но её
подсистема (17), содержащая только приближения параметров РА ф^,
(к = 1, 2,..., N), не замкнута, так как состоит из N + 2 уравнений и содержит 2Nнеизвестных функций. Следовательно, при N ^ 3 система (17)—(19) не распадается так же просто, как это было при N = 2, на отдельные подсистемы, определяющие приближения параметров РА, температуры и скоростей ползучести. При исследовании типа системы (17)—(19) выясняется, что её характеристический определитель тождественно равен нулю. Это вызвано тем, что N (^ 3) уравнений (19) содержат частные производные только
от двух функций г1г), г2г), а значит, должны выполняться некоторые условия совместности этих уравнений. Покажем, что следствиями подсистемы (19) являются некоторые конечные и дифференциальные зависимости между приближениями углов РА ф(г) (1 ^ к ^ N), которые определяют условия совместности этой системы. Для этого перепишем условия РА (19) в форме (см. (1.4) в [2])
(г),(г),(г) Лг)Лг)1{г)
С11 /к! /к!
+ С:
(25)
где
с(г) = I э 2
( )
(г)
г
«кг)=Ск («("0.
(26)
Рассмотрим уравнения (25) как систему линейных алгебраических уравнений (СЛАУ) относительно «У"). Выберем из (25) три уравнения, определитель которых отличен от нуля в некоторой окрестности точки (Ж1,Ж2). Пусть для простоты это будут три первых уравнения. Тогда, используя правило Крамера, из них можно однозначно определить функции
С(г) =
А,
у
ф1г) ,ф('),фЗг),
( )
А (фГ),ф(г),ф(г)
и = 1, 2,
(27)
где
А11 = 2ёе!
С( ) Ск
( ) ( ) к2 к2
( ) ( ) к1 к2
А12 = А21 = ёе!
( ) ( ) к1 к1
( ) ( ) С( ) к2 к2 Ск
А22 = 2 ёе!
( ) ( ) С( ) ( ) ( ) к1 к1 Ск к1 к2
А = 2 ёе!
( ) ( ) ( ) ( ) ( ) ( )
1 к1 1к1 1к2 1 к2 1 к! 1к2
к = 1, 2, 3
(здесь использована построчная запись определителей по к). Следовательно,
(г)
приближения скоростей «У деформаций ползучести являются нелинейными
функциями от фкг), 0(г) (к = 1, 2, 3).
Если количество семейств арматуры больше трех, то, подставив выражения (27) в уравнения (25) при 4 ^ к ^ N, получим N—3 конечных нелинейных
зависимости между приближениями углов РА фПг) (1 ^ п ^ N) и температуры 0(г). Кроме того, функции С" должны удовлетворять безразмерному условию совместности скоростей деформаций ползучести [1]
С(г) + С(г) = 2С(г) . (29)
«11,22 + «22,11 = 2«12,12. (29)
Подстановка выражений (27) в уравнение совместности (29) даёт ещё одну, но уже дифференциальную, зависимость между фкг), $(г) (к = 1,2,3).
Получающиеся при этом конечные и дифференциальные связи совместно с (17), (18) образуют замкнутую систему уравнений относительно приближений параметров РА фкг), (к = 1,2, ...^) и температуры 0(г), а для
определения приближений скоростей г(г), г(г) ползучести можно использовать два любых уравнения (19) или первые соотношения (26) (в зависимости от вида кинематических граничных условий).
Следовательно, преобразованная таким образом система (17)—(19) позволяет вновь на каждой итерации расщепить ранее связанные подзадачи определения параметров РА, температурного поля и поля скоростей установившейся ползучести в пластине и интегрировать их последовательно. Однако основная трудность, возникающая при таком подходе, заключается в том, что при подстановке выражений (27) в уравнение (29) и при дальнейшем раскрытии определителей (28) получим квазилинейное уравнение, содержащее несколько тысяч слагаемых. Очевидно, что анализ такого уравнения потребует неоправданно больших затрат.
Таким образом, проведенное рассуждение показывает, что при внедрении в конструкцию более двух семейств арматуры итерационный процесс (17)— (24) не приводит к принципиальному упрощению исходной системы разрешающих уравнений задачи РА.
Обратимся к более подробному анализу случая укладки двух семейств волокон, так как при этом, как было показано выше, итерационный процесс позволяет существенно упростить процедуру решения задачи РА. Далее в настоящем исследовании будем считать N = 2. (С практической точки зрения условие N = 2 не является обременительным, так как в композитных пластинах несущие элементы — арматура — должны воспринимать соизмеримые напряжения на площадках, различно ориентированных по двум направлениям.)
Правые части в двух первых уравнениях (17) можно трактовать как возмущенные объёмные (массовые) нагрузки, при этом система (17) формально совпадает с системой (3.4), (3.5) в [2], которая там подробно исследована. Здесь же отметим лишь основные свойства системы (17). Согласно [2], эта система квазилинейных уравнений с частными производными первого порядка имеет две двукратные действительные характеристики, определяемые
углами фДг) (k = 1, 2), но не подпадает ни под один тип существующей на сегодняшний день классификации систем уравнений с частными производными [12]. Система (17) обладает свойствами, присущими как гиперболическим (ее нельзя считать гиперболической, так как она не приводится к характеристической форме [2, 13]), так и параболическим (ее нельзя считать параболической, так как в общем случае в силу ф(г) = ф(г) характеристики разных семейств различны) системам уравнений. Однако при отсутствии объёмных нагрузок (Xj = 0, i = 1, 2) на первой итерации (г = 1) в силу выбора начального приближения в виде (24) система (17) распадается на две замкнутые (при каждом k) подсистемы параболического типа (см. (3.14) в [2]). С другой стороны, характеристики системы (17) являются линиями слабых разрывов [2, 13], поэтому, например, решение задачи Коши для этой системы может быть получено лишь в области определенности решения (об этом свидетельствуют и многочисленные расчёты, проведенные автором), что характерно для решений задачи Коши для систем гиперболических уравнений [13].
(r)
Так как углы фд задают траектории РА, а фд — г-тые приближения к фк, то и характеристики системы (17) можно трактовать как r-тые приближения к траекториям армирования. Если итерационный процесс сходится, то характеристики системы (17) совпадут в пределе с траекториями РА, которые, в свою очередь, являются характеристиками системы разрешающих уравнений (1)—(4), как было показано в [2].
Для квазилинейных систем уравнений с действительными характеристиками наиболее естественно и просто ставится задача Коши, поэтому обратимся к вопросу постановки задачи Коши для системы (17). Эта система удовлетворяет условиям теоремы Коши—Ковалевской [12], а значит, при аналитических входных данных решение задачи Коши для нее существует и единственно. Для более широких классов входных данных теоремы существования для системы (17) автору не известны.
Покажем, что статические граничные условия (20) совместно с краевыми
условиями (23) задают краевые условия для функций фДг). Для этого будем считать, что арматура обоих семейств «входит» в пластину на той части контура, где заданы статические граничные условия (Гр = Г^д, k = 1, 2). Тогда
систему трансцендентных уравнений (20), замкнутую относительно фДг) (Гр) (k = 1, 2), путём тригонометрических преобразований можно привести к виду
b + c _,r2,.,2 I ^2, ,2
COS (2<р2 - Ф) =
ь° + c° — + CT°W02
2(T2W02\/b2 + С2 b — CT2W02 sin 2^2
tg 2(^1 =
c — CT2W02 cos 2^2 '
Jk (Гр)
ф((г) (ГР)= ^k + в (k = 12).
cos ф = . , этф = . ,
л/Ъ2 + с2 л/Ъ2 + с2
b = 2pr — єDr (v(s), ^(s), 0(s)) , c = 2P„ — cti^0( — ^2^02,
Pn = Pn — єDra (v(s), <^(s), 0(s)) ,
r = s + 1.
Функция в в (30) задаёт направление внешней нормали к контуру Гр, поэтому, чтобы волокна входили на нем в область С, необходимо решение системы (30) разыскивать на открытом интервале £ (§) ^)- На этом интервале первое уравнение (30) в зависимости от значений ф и правой части может иметь до двух различных корней, второе же уравнение при известном ^>2 имеет только один корень. Следовательно, система (30) может иметь до
двух различных наборов решений относительно (Гр).
Таким образом, статические граничные условия и краевые условия для интенсивностей армирования задают два набора альтернативных краевых условий для системы (17), поэтому в классе аналитических входных данных задача Коши для этой системы может быть корректно поставлена на контуре Гр = Г^, к = 1, 2.
При отсутствии распределённых объёмных нагрузок (Х = 0, г = 1,2) задача Коши (17), (23), (30) на первой итерации (г = 1) может быть проинтегрирована аналитически (см. (3.14)-(3.24) в [2]). При этом из (17) с учётом
(24) в общем случае 8Ш ^ф1г) — ф2г)) = 0 следует
д* (ф^ф^) =0, к = 1, 2, (31)
а значит, характеристики системы (17) (а также и траектории РА в первом приближении) являются прямыми линиями, определяемыми контурными значениями углов ф^1 (Гр), известными из (30).
Замечание 1. Как уже отмечалось, в общем случае нагружения пластины (X = 0) теоремы существования решения задачи Коши (17), (23), (30) в более широких, чем аналитический, классах входных данных автору не известны. Однако при отсутствии массовых нагрузок (Х = 0, г = 1, 2) решение этой задачи Коши на первой итерации существует, если краевые значения для функций фк принадлежат классу гладких функций, а краевые значения для (т. е. ^о^) принадлежат классу непрерывных функций, заданных на части контура Гр. (Об этом свидетельствует структура решения (3.19), (3.21)-(3.24)
в [2].)
Если после интегрирования задачи Коши (17), (23), (30) её область определенности покрывает всю область С, занимаемую конструкцией в плане, то в каждой точке (Ж1,Ж2) € С известны функции Фк , иук . При выполнении
для функций ограничений, аналогичных (11), уравнение (18) является эллиптическим уравнением второго порядка [2] относительно приближения температуры 0(г), которому соответствует тепловое граничное условие (22). Такие граничные задачи и методы их интегрирования хорошо изучены [3, 1416], поэтому не будем останавливаться на обсуждении этих вопросов более подробно.
Если функции ф(г),0(г) (к = 1,2) известны из решения краевых задач (17), (23), (30) и (18), (22), то подсистема (19) является линейной гиперболической системой уравнений с частными производными первого порядка относительно приближений скоростей ползучести, причём её характеристики определяются углами ф(г) и совпадают с характеристиками системы (17). В качестве краевых условий для подсистемы (19) выступают кинематические
граничные условия (21), заданные на части контура rv, где, по предположению, волокна «выходят» из пластины.
Так как область определенности решения задачи Коши (17), (23), (30) покрывает всю область G и характеристики подсистем (17), (19) совпадают, то после интегрирования задачи Коши (19), (21) в каждой точке (ж1, Ж2) € G
(r)
будут известны функции vi , i = 1, 2. (Методы интегрирования задачи Коши для гиперболических систем уравнений достаточно хорошо разработаны [3, 13], поэтому не будем останавливаться на этом вопросе более подробно.)
При отсутствии распределённых объёмных нагрузок (Xj = 0) и постоянстве температуры (0(1) = 0О = const) на первой итерации (г = 1) можно построить аналитическое решение задачи Коши (19), (21). Действительно, выше уже отмечалось, что при Xj = 0 (i = 1,2) характеристики системы (17), а значит, и системы (19) на первой итерации прямолинейны, т. е. выполняются равенства (31). Поэтому систему (19) при 0(1) = 0О = const (в силу постоянства правых частей: = const) можно записать в форме инвариан-
тов Римана [13] с нулевыми правыми частями, после чего легко построить решение задачи Коши (19), (21):
v(1) cos ф^ + v21} sin ф^1 - (®1 cos ф^1 + X2 sin ф^) Cfc =
= v0 (q) cos фк1} (q) + v0 (q) sin фЕ) (q) -- (n1(q) cos ^(q) + П2(q) sinф(1)(q^ Cfc, k = 1,2, (32)
где Cfc = / (afc,0o) = const, = const; ф^ (q) —значения функций ф<^1) на контуре rv, известные из решения задачи Коши (17), (23), (30); (?) —
функции, задающие точки контура rv: ж& = (q), k = 1, 2; q — естественный
параметр вдоль контура. (Так как характеристики системы (19) при Xj = = 0 на первой итерации прямолинейны, то в уравнениях (32) имеют место
равенства ф^ (ж1, Ж2) = ф^ (?).)
Таким образом, на каждой итерации решение задачи РА сводится к последовательному интегрированию задачи Коши (17), (23), (30), граничной задачи (18), (22) и задачи Коши (19), (21). Такое неклассическое сведение краевой задачи РА к последовательности решений нескольких задач Коши возникает из-за того, что система разрешающих уравнений (1)—(4) совместно с граничными и краевыми условиями (6), (8)—(10) образуют краевую задачу с сингулярным возмущением [4], а итерационный процесс (17)—(24) относится к разряду методов теории возмущений.
Для проверки сходимости итерационного процесса целесообразно r-тые приближения неизвестных функций подставить в систему разрешающих уравнений и соответствующие ей граничные и краевые условия (п. 1), а затем проконтролировать нормы возникающих при этом невязок. Следует отметить, что невязки возникают лишь в уравнениях равновесия (1) и статических граничных условиях (6) и имеют вид:
r(r)
о; = е
A (v(r),^<r),0(r)) - A (v(s),^<s),0(s)) (ж1, ж2) € G, i = 1, 2;
— е[^п (^(г),^(г),0(г)) - £п («(в),и(в),0(в))
4г) — Абт (Уг),^(г),0(гЛ - яТ (^(5),
(33)
(жі,ж2) є Гр; г = 8 + 1.
Эти равенства показывают, что определение невязок при использовании итерационного процесса не требует дополнительных вычислительных затрат, так как все выражения, входящие в (33), необходимо вычисляются для построения решения на (г + 1)-й итерации. Автор считает, что на настоящий момент времени контроль норм невязок (33) является наиболее удобным практическим критерием проверки сходимости процесса (17)—(24).
Замечание 2. Вместо кинематических граничных условий (8), (21) могут быть заданы некоторые смешанные условия, позволяющие определить скорости ползучести на контуре. Например, при задании на части Г контура Г, кривизна к которой нигде не равна нулю, нормального напряжения рп и касательной скорости V0 установившейся ползучести можно определить г-тое приближение для нормальной составляющей скорости гп, если задача Коши (17), (23), (30) уже проинтегрирована на г-той итерации. Действительно, используя на Г условия РА (25) (к = 1, 2) и первое статическое граничное условие (6), которое с учётом (7), (26) принимает вид
2ш(г)#0 (я(г), 0(г)) ^1і) (1 + пі) + {22 (1 + «2) + 2£(2)ПіП2
= рп —
со«2 - в), (Жі,Ж2) є Г8, (34)
(г)
можем, как из СЛАУ, определить на части контура Г8. Тогда приближение скорости тангенциальной деформации ползучести на Г будет иметь вид
{(г) — г «іп2 в + г со«2 в — г «іп 2в (35)
{ТТ = Си йіи2 в + {2? СО82 в - {£? йіп 2в.
Эта скорость деформации связана со скоростями точек контура [1, 17]:
{ТТ — (в,г0) +
кг,
(г)
где гПг) — г-тое приближение нормальной составляющей вектора скорости смещений точек контура Г5. Если на Г всюду к — 0, то из последнего равенства можно определить
{ТТ - дп (в, гТ0) , (Ж1,Ж2) Є Г8, г — 8 + 1.
(36)
Следовательно, задание на части контура Г нормального напряжения рп и касательной скорости г0 ползучести позволяет на каждой итерации определить приближение нормальной составляющей гПг) скорости ползучести на
п
к
этом контуре, а значит, и краевые условия для подсистемы (19), если кривизна части контура rs, на которой волокна выходят из пластины, нигде не равна нулю. Подчеркнем, что в этом случае в окрестности контура rs может возникнуть краевой эффект (вследствие сингулярного возмущения исходной
краевой задачи РА). Действительно, в уравнении (34) коэффициенты при содержат в качестве сомножителя малый параметр е, поэтому после разрешения СЛАУ (25), (34) может оказаться, что и достигают по модулю больших значений на части контура rs. Это обстоятельство может породить краевой эффект в связующем и отрицательно сказаться на сходимости итерационного процесса (17)-(24), (36). (Краевой эффект в арматуре не возникает, так как она всюду равнонапряжена: стд = const, k = 1, 2.)
Так как равенства (23), (30) могут определять на контуре Гр до двух альтернативных наборов краевых условий для системы (17) (выбор того или иного набора этих условий определяется только волей проектировщика или из каких-либо дополнительных соображений), удовлетворяющей условиям теоремы Коши—Ковалевской [12], а также в силу единственности решения граничной задачи (18), (22) [14] и задачи Коши (19), (21) [13] при помощи итерационного метода (17)-(24) можно построить до двух различных решений исходной задачи РА, если, конечно, они существуют. Этими решениями можно управлять за счет варьирования контурных функций wofc в краевых условиях (23). (Так как выбор того или иного решения всецело зависит от воли проектировщика, то в процессе решения задачи РА при использовании разработанного итерационного метода не может произойти «несанкционированного» перехода от одного решения к другому, например за счет погрешности численных схем и т. п.) Такое многообразие альтернативных проектов РА выгодно с практической точки зрения, так как позволяет выбрать из их совокупности проекты, обладающие теми или иными эксплуатационными свойствами, а также наиболее удобные с точки зрения технологической реализации.
Замечание 3. Если в пластине с РА-структурой оба семейства волокон изготовлены из одного материала (fi (ст, 0) = /2 (ст, 0)) и напряжения в них одинаковы (cti = СТ2), а на контуре Гр = Г^д в краевых условиях (23) выполняются равенства woi = W02, то РА-структура в обоих возможных проектах совпадает. (С математической точки зрения два решения отличаться лишь тем, что траектории РА разных семейств меняются в них местами; с технологической точки зрения это один и тот же проект.)
Замечание 4. Как уже отмечалось в начале настоящего пункта, решение рассматриваемой задачи РА можно построить и с помощью разложения всех неизвестных функций в асимптотические ряды по малому параметру [4]
(x) = Е (x) em-i, (x) = Е wkm) (x) em-i,
m=1 m=1
(m) („\ „.m-1 (37)
0 (x)= E 0(m) (x) 1, V (x) = E V (x) 1,
m=1 m=
i = 1, 2, k = 1, 2,... , N.
Асимптотический анализ системы разрешающих уравнений (1)—(4) и соответствующих ей граничных (6), (8), (9) и краевых (10) условий с учётом
разложений (37) показывает:
1) первые слагаемые (ф(1), Ш(1), 0(1), ) в рядах (37) полностью совпа-
дают с решением системы (17)—(24) на первой (г = 1) итерации;
2) если на второй (г = 2) итерации метода (17)—(24) левые части уравнений (17)—(20), (22) линеаризовать в окрестности решения для первой итерации, то получающаяся система и краевые условия будут совпадать с уравнениями и краевыми условиями, определяющими коэффициенты при вторых слагаемых (фк2), шД2, 0(2), -и(2)) в разложениях (37).
На основе столь существенного сходства уравнений итерационного процесса (17)—(24) и метода асимптотических рядов (37) все результаты исследований, полученные ранее в настоящем пункте при анализе уравнений (17)—(24), можно перенести и на соответствующие цепочки уравнений, получающиеся после подстановки (37) в исходную систему разрешающих уравнений и краевые условия задачи РА (п. 1). Однако итерационный процесс (17)-(24) обладает рядом преимуществ перед методом асимптотических рядов:
1) правые части в уравнениях (17)-(23) на каждой итерации имеют одну и ту же достаточно простую структуру (в цепочках же уравнений метода асимптотических рядов (37) структура правых частей с увеличением номера уточнения т существенно усложняется);
2) для построения г-того приближения по методу итераций используется информация лишь о предыдущем (г — 1)-м приближении (в правых же частях уравнений метода асимптотических рядов используется информация о всех предыдущих п-тых уточнениях, п = 1, 2,... , г — 1);
3) характеристики подсистем (17), (19) изменяются от итерации к итерации и при сходящемся итерационном процессе совпадают в пределе с траекториями армирования (в уравнениях же метода асимптотических рядов (37) характеристики в соответствующих подсистемах уравнений фиксированы и совпадают с характеристиками подсистем (17), (19) на первой (г = 1) итерации);
4) для проверки сходимости обоих обсуждаемых методов целесообразно г-тое приближения по методу итераций и г-е частичные суммы типа (37) подставить в систему разрешающих уравнений и соответствующие ей краевые условия (п. 1) и проконтролировать нормы возникающих при этом невязок, причём в случае метода асимптотических рядов (37) эти невязки возникают во всех уравнениях и краевых условиях исходной системы (п. 1) и с увеличением номера уточнения (с увеличением г) структура их существенно усложняется, в случае же итерационного процесса невязки возникают только в уравнениях (1), (6) и имеют достаточно простую структуру (33).
Указанные преимущества позволяют отдать предпочтение итерационному процессу (17)-(24) при практическом решении задач РА пластин в условиях установившейся ползучести.
Замечание 5. Если задача РА уже решена (известны функции фд, шд, 0, ■ид, к = 1, 2), то по формулам (1.2)—(1.4) из [2] можно определить напряженное состояние в связующем и в композиции (сту) пластины в условиях установившейся ползучести. Кроме того, зная фд, шд, 0 (к = 1, 2), можно определить и начальное напряженное состояние сту в пластине (его можно вычислить, например, методами, изложенными в [18, 19]). Решение же о неустановившейся
ползучести на первой её стадии можно получить приближённо, используя методику, предложенную в [1]. Для этого представим текущие осреднённые напряжения в пластине в виде
(j = (1 - т (t)) (ij + т (t) (ij, i, j = 1, 2,
где т (t) — монотонная функция времени t, возрастающая от значения т = О в начальный момент времени t = О до т = 1 при t ^ то. Функцию т (t) можно определить исходя из экстремальных принципов теории ползучести
(см. § 34 в [1]). В случае изотропных пластин реализация такого подхода
не вызывает особых трудностей (так как можно почти дословно повторить все рассуждения § 34 в [1]), а в случае сложно армированных пластин (как это имеет место в задачах РА) требует проведения дополнительного более тщательного исследования, выходящего за рамки настоящей работы.
З. Обсуждение результатов расчётов. В качестве примера рассмотрим расчёт РА-структур при установившейся ползучести металлокомпозитной двусвязной пластины, ограниченной контурами, уравнения которых заданы в безразмерной полярной (r, ф) системе координат:
Гр = Гл = Г^2 : Гі (<£>) = О,Б - О,ОБ cos (2^ - п/2),
Г, : r2 (^>) = 1 + 0,08cos (2^>), О ^ ^ < 2п.
Внешний контур Г, пластины жестко закреплен (vQ (^>) = vQ (^>) = 0), а на внутреннем контуре Гр заданы равномерно распределённые контурные нагрузки или рт. Распределенные объёмные нагрузки отсутствуют (Xi = 0,
i = 1, 2), температурное поле предполагается однородным (О = 0Q = const). Пластина изготовлена из меди и армирована двумя семействами стальной проволоки У8А (механические характеристики установившейся ползучести фаз композиции имеют значения (13), (14)). Плотности армирования проволками обоих семейств на внутреннем контуре имеют значения
W1 (Гр) = w2 (Гр) = WQ1 (^) = WQ2 (^>) = 0,3Б. Гр = Г^1 = ТШ2.
На рис. 2 изображена РА-структура, полученная в случае нагружения внутреннего контура давлением = —Б00 МПа (рт = 0) при задании уровня напряжений в арматуре <ті = (2 = —1 ГПа. Согласно замечанию 3, в этом случае структура, соответствующая второму решению задачи РА, полностью совпадает с изображенной на рис. 2. (Об этом свидетельствуют и результаты проведенного расчёта.)
На рис. 3 изображены две альтернативные РА-структуры, полученные в случае нагружения внутреннего контура касательной нагрузкой рт = = —2Б0 МПа (pn = 0) при задании уровня напряжений в арматуре Ті = = —(2 = 1 ГПа. При этом существуют оба решения задачи РА, и, как видно из рис. 3, структуры армирования в них существенно отличаются. Однако, кроме структуры, отличие заключается и в напряженном состоянии, возникающем в связующем. Так, в проекте, изображенном на рис. 3, б, максимум интенсивности напряжений в связующем на 43,75 % больше той же величины в проекте, представленном на рис. 3, а. Кроме того, в указанных проектах отличаются поля скоростей установившейся ползучести и общий расход арматуры, который в случае, изображенном на рис. 3, б, на 44,28 % меньше, чем в РА-структуре, приведённой на рис. 3, а.
Рис. 2. Структура равнонапряженного армирования двусвязной пластины, нагруженной
на внутреннем контуре давлением
Рис. 3. Два альтернативных проекта равнонапряженного армирования двусвязной пластины, нагруженной на внутреннем контуре касательными нагрузками
Подчеркнём, что возможность построения двух альтернативных РА-про-ектов определяется тем, что задачи РА относятся к разряду обратных задач механики композитов (а точнее, к специфическим коэффициентным задачам), которые, как известно, могут иметь неединственное решение.
Возможность построения двух решений задачи РА можно продемонстрировать на следующем простейшем примере. Рассмотрим прямоугольную удлиненную пластину (в идеале бесконечной длины) шириной l, в которой закрепление, нагружение, температурное поле и структура не изменяются в продольном направлении ж2. (Пренебрегая торцевыми эффектами, в этом случае получим решение задачи РА, зависящее только от переменной жь) Тогда при однородном температурном поле (в = во = const), отсутствии объёмных
нагрузок (Xi = 0, i = 1, 2) и задании на продольном контуре x1 = l > 0 постоянной касательной нагрузки рт = const = 0 (pn = 0) из (1), (6), (7) получим
(в = 0)
Е OkWk cos2 фк + 4еадо (H, в) = Pn = 0,
k
Е OkWk sin 2фк + 2еадо (H, в) V2,i = 2рт, (38)
к
в = в0 = const, 0 ^ x1 ^ l.
Условие РА (4) принимает вид
v1,1 cos2 фк + 0,5v2,1 sin2^k = £k = const, k = 1, 2. (39)
Рассмотрим РА-структуры, в которых W1 = W2, ф1 = —ф2, O1 = — O2, £1 = = — £2, тогда из (39) получим ^1,1 = 0 и первое уравнение (38) выполняется тождественно, а второе преобразуется в
o1w1 sin 2^ + еад0 (H, в0) v2,1 = pT, 0 ^ x1 ^ l. (40)
В силу малости параметра е отбросим в первом приближении второе слагаемое в левой части равенства (40), тогда получим
o1w1 sin 2^ = pT, 0 ^ x1 ^ l. (41)
Отсюда, задав изначально рт, о1 и w1 (так как w1 = w01 = const, см. (10)), можем определить угол армирования
ф\ = — arcsin (———— J , 0 ^ x\ ^ I. (42)
2 \01W1 /
Однако помимо (42) решением уравнения (41) является и ф\ = \ — ф 1, причём эти решения совпадают лишь в частном случае ф\ = ф\ = j, т. е., согласно (42), при рт = o1w1. В общем же случае (|рт| < |o1 W1|) уравнение (41) имеет два альтернативных решения (^ = ^1), которым соответствуют две различных РА-структуры (как это имеет место на рис. 3).
Работа выполнена при поддержке РФФИ (код проекта 07-08-00152-а).
БИБЛИОГРАФИЧЕСКИЙ СПИСОК
1. Качанов Л. М. Теория ползучести. — М.: Физматгиз, 1960. — 456 с.
2. Немировский Ю. В., Янковский А. П. Равнонапряженное армирование металлокомпозитных пластин волокнами постоянного поперечного сечения в условиях установившейся ползучести // Механика композитных материалов, 2008. — Т. 44, №1. — C. 1134.
3. Самарский А. А. Теория разностных схем. — М.: Наука, 1989. — 616 с.
4. Найфе А. Введение в методы возмущений. — М.: Мир, 1984. — 536 с.
5. Писаренко Г. С., Можаровский Н. С. Уравнения и краевые задачи теории пластичности и ползучести: Справочное пособие. — Киев: Наукова думка, 1981. — 496 с.
6. Композиционные материалы: Справочник / ред. Д. М. Карпинос. — Киев: Наукова думка, 1985. — 592 с.
7. Карпинос Д. М., Невгод В. А., Тучинский Л. И. и др. Ползучесть и длительная прочность вольфрамовых проволок// Пробл. прочности, 1972. — №1. — C. 70-73.
8. Соснин О. В. О ползучести слабо упрочняющихся материалов при нестационарных температурно-силовых режимах// Пробл. прочности, 1972. — №1. — C. 74-77.
9. Джураев Т. Д. Краевые задачи для уравнений смешанного и смешанно-составного типов. — Ташкент: Фан, 1979. — 238 с.
10. Новожилов В. В., Черных К. Ф., Михайловский Е. И. Линейная теория тонких оболочек. — Л.: Политехника, 1991. — 656 с.
11. Камке Э. Справочник по дифференциальным уравнениям в частных производных первого порядка. — М.: Наука, 1966. — 260 с.
12. Петровский И. Г. Лекции об уравнениях с частными производными. — М.: Физматгиз, 1961. — 400 с.
13. Рождественский Б. Л., Яненко Н. Н. Системы квазилинейных уравнений. — М.: Наука, 1969. — 592 с.
14. Бицадзе А. В. Краевые задачи для эллиптических уравнений второго порядка. — М.: Наука, 1966. — 204 с.
15. Беляев Н. М., Рядно А. А. Методы теории теплопроводности: Учеб. пособие для вузов. В 2-х частях. Ч. 1. — М.: Высш. школа, 1982. — 327 с.
16. Беляев Н. М., Рядно А. А. Методы теории теплопроводности: Учеб. пособие для вузов. В 2-х частях. Ч. 2. — М.: Высш. школа, 1982. — 304 с.
17. Демидов С. П. Теория упругости. — М.: Высш. школа, 1979. — 432 с.
18. Немировский Ю. В., Резников Б. С. Прочность элементов конструкций из композитных материалов. — Новосибирск: Наука, 1986. — 168 с.
19. Немировский Ю. В., Янковский А. П. Термоупругопластический изгиб сложно армированных пластин // Механика композитных материалов, 2005. — Т. 41, №6. — С. 715742.
Поступила в редакцию 27/11/2009; в окончательном варианте — 16/УШ/2009.
MSC: 74K20, 74R20
APPLICATION OF METHODS OF THE THEORY OF PERTURBATIONS IN FLAT PROBLEM OF EQUALLY-STRESSED REINFORCING OF METAL-COMPOSITE PLATES AT STEADY CREEP
A. P. Yankovskii
Khristianovich Institute of Theoretical and Applied Mechanics,
Siberian Branch of the Russian Academy of Sciences,
630090, Novosibirsk, Institutskaya st., 4/1.
E-mails: nemirov@itain.nsc.ru, shulgin@itam.nsc.ru
The problem of equally-stressed reinforcings of the plates loaded in their plane and working in condition of steady creep is solved by methods of the theory of perturbations. The possibility of existence of several alternative solutions of the problem under consideration is demonstrated which can be selected reliably, using the algorithm developed. Explicit projects of equally-stressed reinforcings of a doubly-connected plate are constructed at different types of its load.
Key words: plates, steady creep, rational projection, equal-stress reinforcing, small parameter, perturbations theory.
Original article submitted 27/II/2009; revision submitted 16/VIII/2009.
Andrey P. Yankovskii (Dr. Sci. (Phys. & Math.)), Leading Research Scientist, Physics of Fast Process Laboratory.
