Для нахождения первого приближения по формуле (39.10) потребуются значения рІ^ФоОО1)]- Чтобы вычислить эту функ цию, достаточно определить из (39.12) Фо(т}) = arg юо(Ç) и под ставить найденный аргумент в правую часть соотношения (39.11) вместо Ф. В результате будем иметь
G (#) = |
Л'*0/?1 = , |
6 0 1 |
(39.14) |
Задавая |
теперь в интервале |
(0, я) |
п узлов |
|
|
q, |
= |
2v—1 |
, |
0 |
|
|
Щ |
~2я~ я, |
v = l , 2, . |
|
построим |
четный |
|
тригонометрический полином |
G„_i (f>) = |
л—1 |
|
|
|
|
|
|
= ^ Лѵ cos vù, интерполирующий |
G(f>) в данных |
узлах. По- |
Ѵ=0
лином этот согласно формуле (39.14) не будет содержать чле
нов с нечетными индексами. Обозначая, далее, через |
G„_i (т>) |
тригонометрический полином, сопряженный с Gn-ifö), |
найдем |
первое приближение искомой отображающей функции в окон
чательном |
виде: |
|
|
|
|
|
|
<ві (а) = |
(ùg (o)gn-i |
(о), |
g„-i |
(а) = G„_, (а) + iG*+l |
(Ü). (39.15) |
|
Численные |
расчеты были проведены для значений |
полуосей |
|
2 |
|
9 |
|
|
|
|
|
а — —т= , о = -^=г ; |
было |
взято |
д = 1 2 . Для первого |
прнближе- |
ния было найдено |
выражение |
|
|
|
|
|
|
|
о (£) = Çe**+e V сг 2 ѵ £ 2 ѵ , |
|
(39.16) |
где |
|
|
|
|
ЙО = — 0,003604, |
|
|
а = |
0,127706, |
а0 = 0,976450, |
а., = 0,022810, |
ß = |
0,016136, |
а6 =0,003460, а8 =0,000703, |
а1 0 =0,000134. |
|
Кривая, являющаяся изображением окружности по первому |
приближению |
(39.16), отклоняется от заданного эллипса не бо |
лее, чем на 0,006345, а максимальное отклонение для началь
ного приближения имеет |
порядок 0,048524. Первая |
итерация |
сближает, таким образом, |
изображение |
нулевого приближения |
к истинному контуру примерно в восемь |
раз. Отметим, |
что ото |
бражение, найденное по методу Л. В. Канторовича в рассмат
риваемом примере, отходит от |
точного примерно на 0,02 |
(Л. В. Канторович и В. И. Крылов |
[1], стр. 441). |
Глава пятая
О ПЛОСКИХ ЗАДАЧА Х НЕСИММЕТРИЧНОЙ УПРУГОСТИ
Согласно концепции |
братьев Коссера (Е. Cosserat, F. Cosse- |
r a t [1])» учитывающей |
вращательное |
взаимодействие частиц ма |
териала, при изучении |
напряженного |
состояния твердого дефор |
мируемого континуума необходимо наряду с обычными напря
жениями (сила на единицу площади) вводить |
в рассмотрение |
и моментные напряжения (моменты на единицу |
площади). |
Наличие моментных напряжений предполагает несимметрич ность тензора (обычного) напряжения и приводит к рассмотре нию градиентов деформации высоких порядков. В случае плос кой деформации изотропного тела 1 ) появляется одна дополни тельная упругая постоянная, имеющая размерность длины и свя зывающая моменты с вызванными ими искривлениями.
Определяющие уравнения и граничные условия, соответст вующие тому или иному состоянию равновесия, имеют более сложный вид по сравнению с классическим случаем. Число гра ничных условий (геометрических, статических и смешанных), необходимых для определения упругого равновесия тела, равно трем. Число возможных комбинаций задания этих условий зна чительно больше числа тех же комбинаций в классической уп ругости.
Интерес к теории моментных напряжений значительно возрос за последнее время и различные аспекты ее стали предметом изучения многих авторов2 ). Основная часть имеющихся к нас тоящему времени солидного числа научных публикаций посвя щена общим теоретическим вопросам обоснования теории.
') Что мы и будем иметь в виду-во всем дальнейшем.
|
|
|
|
|
2 ) Укажем здесь |
на работы Мпндлниа |
[1—3], Миндлина и Тирстеиа [1], |
Трусделла и Туппна |
(Truesdell and Toupin |
[I]) и Тупина (Totipin [I]). Из ра |
бот советских авторов назовем |
работы |
Э. Л. Аэро и Е. В. Кувшинского [I] |
и Г. Н. Савина [3]. Довольно |
полную |
библиографию по указанным вопро |
сам можно найти в обзорной статье Г. Н. Савина и Ю. Ы. Немиша [1].
В настоящей главе мы приведем основные соотношения не симметричной теории упругости и укажем применение аппарата теории функций комплексного переменного к решению гранич ных задач в простейших случаях.
§ 40. Основные соотношения плоской моментной теории. Постановка граничных задач
1°. К обычным уравнениям равновесия тела, связывающим четыре компоненты напряжения, присоединяется одно уравнение равновесия моментов Мх и Ми. Система уравнений равновесия в случае отсутствия объемных сил здесь имеет вид')
да |
дт,,,. |
0, |
от,.,, |
да,. |
п |
дМ |
дМ„ . |
XIJX = |
п |
0. |
, л п .. |
-л* + |
- Г * = |
+ |
— ^ = |
0, |
-1- |
« - I - гХу - |
|
40.1) |
дх |
ду |
' |
дх ' |
ду |
|
дх |
ду |
|
|
|
|
Ввиду несимметричности тензора напряжений касательные напряжения имеют симметричную и антисимметричную состав ляющие T s и га:
Симметричная составляющая т.„ так же как и в классическом случае, вызывает обычную деформацию сдвига ъ х у ,
х = 2 ^ х у , |
г х у = ± 2 [дх ' ду |
(40.3) |
где ц. — модуль сдвига, |
а антисимметричная |
составляющая ха, |
уравновешиваемая моментнымн напряжениями, поворачивает элемент на угол ы,
Моментные напряжения Мх, Ми искривляют элемент, вследст вие чего появляются пропорциональные, по предположению, этим напряжениям кривизны хх, у.„, определяемые формулами
да |
<Эш |
|
% х ~~ ~дх~ ' |
КУ ~ ~ду ' |
(40.5) |
где В — постоянная, называемая |
модулем изгиба. |
|
За новую постоянную материала, характеризующую |
степень |
влияния моментных напряжений, здесь принята имеющая раз
мерность длины величина /, определяемая |
равенством |
1°- = — . |
(40.6) |
') Содержание настоящего пункта заимствовано |
из работМипдлииа [1, 3]. |
Поле напряжений (обычных и моментных) выражается через две функции напряжений U и V по формулам
d2U |
д-Ѵ |
|
|
_ |
d4J_ _, |
д^Ѵ_ |
|
|
а * = ду* |
дх ду |
' |
аУ |
~~ |
дх'- г |
дх ду ' |
(40.7) |
d*-U - d-W |
' ^ _ |
д'-Ц ^ dW_ |
|
дх |
ду |
ду"-' |
V |
= |
дх |
ду 1 |
дх2 ') |
|
|
|
дѴ_ |
|
|
|
|
|
(40.8) |
|
'' дх ' |
|
У |
ду |
|
|
|
|
|
|
|
Функции напряжении U и V удовлетворяют |
системе |
урав |
нений |
|
|
|
|
|
|
|
|
д-{Ѵ-РЬѴ) |
= |
2 ( 1 - Ѵ ) Р А Д І / , |
|
дх |
|
|
|
|
|
|
|
|
j-(V-l*AV) |
= |
2 ( \ - v ) |
l ^ A U . |
|
|
Из этих уравнений легко заметить, что функция U — ангармо ническая, а V удовлетворяет уравнению
АѴ—12А2Ѵ—0.
Для упругих смещений и, ѵ и напряжений (40.7), (40.8) имеют место следующие формулы общего комплексного пред ставления, установленные Р. Д. Миидлиным [3] ' ) :
2р. (и + іѵ) = хер (г) — zw' |
(z) |
— яр (z) — |
|
|
|
|
|
- 8 ( l - v ) |
I4p"(z) |
+ |
2 i ^ , |
|
|
|
|
|
|
|
dz |
|
ox + |
(jy |
= |
2 [w' (г) + |
VW)}, |
|
• ax + i (rxv |
+ |
x,JX) |
= |
2 [zw" |
(г) - f я|/ (г) |
+ |
|
8(1 - v ) / у " ( г ) + 2/ â 2 - |
|
+ |
= 8 t ( l -v)l2w"(z) |
+ 2^L, |
x = 3 - 4 v . [ ( 4 0 Л 1 ) |
|
|
dz |
j |
Здесь ер и яр обозначают потенциалы Колосова — Мусхелишвили,
2U = z~w(z)+zVJ) |
+ |
X(z) |
+ |
W), |
%'{*) = №), |
(40.12) |
v — коэффициент |
Пуассона, |
Я — решение уравнения |
Гельм- |
гольца |
|
|
|
|
|
|
|
|
АН |
— К*Н = 0, |
|
% = |
~ ; |
(40.13) |
а — обозначает |
операцию |
|
|
|
|
|
dz |
|
|
|
|
|
|
|
|
д |
|
|
|
д |
|
(40.14) |
|
д'г |
= ~\1х~ |
1' |
1 ду |
|
|
|
|
') В названной работе имеются опечатки.
Наконец, функция V допускает представление (см. там же)
у = - 4t ( 1 - V) Г- [ср' (г) - |
+ H. |
(40.15) |
Выпишем еще известные зависимости между упругими эле ментами в полярной H декартовой системах координат:
ог |
+ сто = ах |
-!- Oy, |
|
|
Т г в — Т д г |
|
|
(40.16) |
о в |
— сгг - f і (тг* + тег) = К |
— ст* -I-1 (тхѵ - f т„Л )] e2,'», |
|
|
|
|
|
(40.17) |
|
ог |
+ iü» = (и + |
іѵ) е-''*, 2 = x + iij re1 |
(40.18) |
Ha основании этих формул и формул (40.9) — (40.11) по лучим следующие формулы комплексного представления реше ний уравнений теории моментных напряжений в полярных ко ординатах:
|
ar |
- f о* = |
2 [ср/ (г) + |
ср' (г)], |
|
с в - |
аг + і (хгь + |
т*г ) = |
2 [г<р" (г) + |
і|/ (г) |
+ |
(40.19) |
|
- 8 ( 1 - ѵ ) / Ѵ " ( 2 ) - Ь 2 ^ г , |
|
|
|
|
|
8(1 —v) / У ( г ) + 2i |
dz |
Ol* |
(40.20) |
2ц (о, + ш») = |
[xcp (2) - 2 Ф ' (г) - xp (2) - 8 (1 - |
v) /чр" (г) |
+ |
|
|
|
- f 2t. дНdz |
|
(40.21) |
2°. Подсчитаем компоненты вектора напряжения Jncfs, Fn ds, действующие на элемент дуги какого-либо профиля AB, прове денного в плоскости Оху. Имеем
Х„ = |
о* cos пх -•- % ѵ COS tt!/, j |
(40.22) |
|
|
К,, = |
Txy COS ПХ + OyCOS /If/, |
|
где n обозначает нормаль к кривой AB, проведенную вправо по отношению к положительному на AB направлению. Замечая, что на AB
^ |
dij |
^ |
dx |
COS ПХ = |
—T' |
cos n y — |
|
|
as |
|
|
из формул (40.22) и (40.7) выводим
|
A " ~ ds 1 ду |
|
V |
— |
_ |
dU |
dV |
|
дх |
dx |
ду |
|
1 |
" ~~ |
ds |
|
|
|
|
|
Складывая предыдущие равенства, после умножения второго из них на і, и интегрируя полученную сумму по s, получим
2 - І - {U — іѴ) = |
if (X„ + iYn) ds + const. |
dz |
•„ |
Внесем сюда вместо U и V их выражения (40.12) и (40.15). Тогда получим следующее равенство, справедливое вдоль лю бого профиля в физической области: •
<p(z) + z с,/ (г) + -ф (z) + 8 ( 1 - |
V) /я ср" (г) - |
,. |
дН |
= |
2і |
|
=г f (X„ 4- iYn) ds 4- const. (40.23} 'о
Тем же путем, исходя из формул (40.8) и используя, пред ставление (40.15) для моментиого напряжения М„, действующе го на тот же профиль А В,
находим |
|
Мп = Мх |
cos юс 4- My cos inj, |
|
|
(40.24) |
|
|
|
|
|
|
|
|
|
|
|
|
Мп |
= |
- Re f - ^ - |
f 8 (1 - |
V) /2 |
<j?4i) - 2i ^ |
- |
] | . |
(40.25)- |
|
Приведем, |
наконец, |
формулу |
комплексного |
представления |
для |
угла поворота |
со, даваемого |
равенством (40.4). На |
основа |
нии формулы |
(40.9) |
имеем |
|
|
|
|
|
|
со = Im | — |
(и + Ц |
= - ^ - Im {(и 4- 1)Ф ' (г) + - ^ г я]. |
(40.26) |
|
Рассмотрим |
теперь |
интеграл |
|
|
|
|
|
|
|
|
' |
= |
f (Х,« и 4- У„ V + Ма со) ds, |
|
|
(40.27) |
где L обозначает контур физической области 5, предполагаемой |
конечной. |
|
|
|
|
|
|
|
|
|
|
|
|
Заменив под интегралом Хп, У„, Мп |
их значениями из формул |
(40.22), (40.24), преобразуем интеграл |
по формуле |
Грина и вос |
пользуемся уравнениями |
равновесия (40.1). Будем |
иметь |
, |
[' С Г да . |
дѵ , |
|
дѵ . |
du . |
|
|
|
|
|
|
4 |
(%< - |
те д ) со + |
М, - ^ - 4- Мѵ Щах dij. |
(40.28) |
|
Учитывая |
(40.4) |
и (40.3), находим |
|
|
|
|
тх„ |
+ т,, ^ |
|
+ (%< - |
т ж у ) со = |
( - * • 4- - І ^ ) = |
Ч^%- (40.29> |
Сумму первых двух членов под интегралом (40.28) выразим согласно закону Гука через деформации sx, гѵ, примем во.
внимание равенство (40.29), а для двух последних членов в ин теграле воспользуемся формулами (40.5). Тогда упомянутый интеграл представится в виде
|
/= 2\\Wdxdij, |
(40.30) |
|
s |
|
|
где |
|
|
|
2W == а.(е* + еу )2 + |
2р{г1 + гІ + 2г%) + 4ß [4 |
+ xf,). (40.31) |
Выражение W — неособенная |
положительная |
квадратичная |
форма совокупности |
компонент |
деформаций |
(как обычных, |
так и вызванных моментными напряжениями),— представляет собой потенциальную энергию упругой деформации, заключен ную в единице объема.
Отсюда немедленно следует равенство нулю всех шести ком понент напряжения ах, ау, ххѵ, хух, Мх, М„ во всей области, если на
ее границе выполняется одна из двух следующих |
групп условий: |
X,~Yn=Mn |
= 0, |
(40.32) |
и = о = (й = |
0. |
(40.33) |
По условиям (40.32) решение будет определено |
с точностью до |
жесткого смещения тела как целого, |
|
и = — гу-\-а, |
w = e.v4-ß, |
|
где е, ос, ß — постоянные.
Задачи определения упругого равновесия тела при неодно родных условиях (40.32) и (40.33) будем по аналогии с класси ческим случаем называть соответственно первой и второй основ ными задачами моментной теории упругости.
Выше мы доказали единственность решений этих задач для конечной области. Приведенное рассуждение с известным огра ничением поля напряжений на бесконечности справедливо и для случая бесконечных областей.
Так же, как и в классическом случае, для существования ре шения первой основной задачи в конечных областях необходимо выполнение условий статики твердого тела. В нашем случае эти условия имеют вид
f (Хя + iYn) ds = 0, |
\ (xYn - УХ„ + М„) ds = 0. |
(40.34) |
L |
L |
|
Они непосредственно следуют из уравнений равновесия |
(40.1). |
Указанные выше комплексные представления полей напряже ний, смещений и деформации позволяют теперь сформулировать основные граничные задачи теории моментных напряжений в
терминах теории |
функций. |
|
|
|
|
|
П е р в а я о с н о в н а я |
з а д а ч а . Найти |
аналитические |
в |
рассматриваемой |
области |
функции |
cp(z), ^(z) |
от |
комплексного |
аргумента |
|
z и |
регулярное |
в |
той же |
области |
решение |
H |
уравне |
ния |
(40.13) |
|
по |
граничным |
|
|
условиям |
|
|
|
|
|
|
|
|
ф (t) |
+ |
t Ф' (0 + |
Ф (О -+8(1 — V ) Z* ф " (О - |
,. дН |
= |
|
|
|
21 |
|
|
|
|
|
|
|
|
|
|
|
|
|
= М О + іМ О + const, |
(40.35) |
|
|
|
|
|
|
! ( 1 - ѵ ) Р Ф " ( 0 - 2 і | ' |
= М 0 и а L > |
|
|
где |
fi, |
/г, |
fo — |
заданные |
на |
|
контуре |
L |
вещественные |
функции |
от |
t, определяемые |
|
через |
внешние |
усилия |
и |
моменты, |
действующие |
на |
границу |
|
тела, |
|
равенствами |
|
|
|
|
|
|
|
|
|
|
|
|
|
fi |
+ ih |
= i Г (Хп |
+ |
îYn) |
ds, |
f0 = |
- M„. |
|
(40.36) |
|
Второй |
|
основной |
|
задаче |
|
соответствуют |
граничные |
|
условия |
|
|
|
х Ф |
(0 - t |
ф' (0 — т|з |
(0 — 8 (1 — ѵ ) / 2 |
ф " (0 + |
|
|
|
|
|
|
|
|
|
|
|
|
|
2*-=- = 2 | i ( g x |
- i - igs ), |
(40.37) |
|
|
|
|
|
|
|
|
|
|
|
|
а/ |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Im |
|
(х + |
1 ) Ф ' ( 0 |
т |
Ш |
= |
2j.tg-0 |
на L , |
|
|
|
|
|
|
|
|
2/2 |
|
|
|
|
|
|
|
go — заданные |
|
смещения |
|
точек |
контура |
L , |
g0- |
•за- |
данные |
на |
L |
значения |
|
угла |
|
|
поворота: |
|
|
|
|
|
|
|
|
|
|
u+iv=gl-\-ig2, |
|
|
a — go |
|
на |
L . |
|
(40.38) |
Наряду с условиями указанных двух основных задач можно, разумеется, ставить и другие граничные условия, например, задавать на границе области и, ѵ, Мп и т. д. Среди этих других задач мы укажем лишь одну, соответствующую соприкасанию тела с жестким профилем заданной формы при отсутствии сил трения (см. § 1). В качестве третьего граничного условия при мем здесь наличие углов поворота точек границы и запишем условия задачи в виде
Г = 0 , vn=?g(t), |
a = h(t) |
на L , |
(40.39) |
где Т — касательная составляющая |
вектора напряжения, а о , — |
нормальное |
смещение; |
g(t) |
и h(t) |
обозначают |
заданные |
па L |
функции. Эту задачу мы будем условно называть третьей |
задачей |
момеитной |
теории упругости. |
|
|
|
|
Исходя |
из формул |
(40.35) |
и (40.37), можно |
без особого |
тру |
да построить эффективные решения основных задач для облас тей, ограниченным одной окружностью или двумя концентриче скими окружностями. Некоторые простейшие задачи такого рода будут рассмотрены в следующих параграфах.
§ 41. Решение первой основной задачи для бесконечной плоскости с круговым отверстием1 )
Центр отверстия радиуса R возь'мем в начале координат. Условия (40.35) без несущественной постоянной в правой части примут в нашем случае вид
Ф(0 |
+tVTt) + W) + m ц |
Щ |
= |
h •!- i-h, |
|
R e ^ |
imp" {t) - |
2і |
ІЩ- |
= - /„, |
t = Re'*, |
|
|
|
|
ді |
|
|
где принято обозначение |
|
|
|
(41.1) |
|
|
|
|
|
|
ш = 8(1— ѵ)12. |
|
|
(41.2) |
Мы |
будем считать, |
что поле |
па |
бесконечности |
однородно2 ) |
и, кроме того, главный вектор усилий, приложенных к обводу
отверстия, равен |
нулю. |
|
|
|
Положим |
в области вне контура отверстия |
|
|
|
|
со |
|
со , |
|
|
ф(г) = |
Г г - ; - v a * * " * , а|,(г) = |
Г'г - f V ак г~к |
(Г = Г), |
(41.3) |
|
|
i |
|
о |
|
|
H (x, |
у) = |
v |
Ай/с* (Xr) |
г - /хм1 ?2 , |
= hk, |
(41.4) |
|
|
— со |
|
|
|
|
где Кп(Кг) обозначает функцию Бесселя второго рода от мни мого аргумента (функция Макдопальда).
Выражая операцию —— в полярной системе координат, при меняя ее к правой части выражения (41.4) и используя правило дифференцирования функции Кп (см., например, Ватсон [1]), получим
Тг = ~ 2- S Mft-i Kk frr) е^. |
(41.5) |
—со
Ряды для производных от функции (p(z) мы получим, оче видно, дифференцированием соответствующее число раз по z ряда (41.3) для ф(г). Составим все производные нужных поряд ков и положим, кроме того, на окружности
U (0 + Ш |
= |
Ъ |
Ake*\ |
|
|
— СО |
(41.6) |
|
|
со |
— /о (0 |
= |
2 |
с-л e'"ft!>, cc_fe = a.k. i |
') Задача для спободного от усилий отверстия'в однородном поле была решена в цитированной выше работе Миндлина [1]. Решения некоторых ана логичных задач имеются в работах Вейтцмана [I] и Хартренфта и Си [1].
2 Ï JVloMeimibie напряжения на бесконечности везде принимаем равными нулю: М[т) = M œ ) = 0.
Все предыдущие ряды будем считать равномерно сходящи мися на окружности и в этом предположении внесем их в усло вия (41.1). Мы будем иметь после некоторых приведении
со |
|
|
|
со |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
y i a |
k R - k e - i k u |
- y |
i ( k - 2 ) R - k + 2 |
[1 — m (ft - |
1) |
tf~2] ö A _ 2 |
е і |
к Ь 4- |
|
i |
CO |
|
|
з |
СО |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
+ |
Jj.^LR~K ß'*9 + |
X |
Уі |
Cfe-i |
Kk |
(Kr) е«» +2ГRe'» 4- P R e - ' » |
= |
|
|
|
|
|
|
|
|
|
|
|
|
|
— |
v |
A, |
|
|
|
|
(41.7) |
|
|
|
|
|
|
|
|
|
|
l + |
|
— |
±j rt-k 1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
l S |
tu IKk+i |
(Kr) |
+ |
|
|
|
|
|
|
|
|
|
|
Kk-i |
(Kr)]e'*> |
= H a t e ' f c » , |
|
|
(41.8) |
причем |
|
|
|
|
Ск |
= |
t'/lft, |
c_/e = |
— cA. |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Сравнение |
в |
этих равенствах |
коэффициентов при е'п * |
даст: |
|
|
|
|
|
|
|
a'Q — Яс^о = Л0 , |
j |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
a/CiC0 = a0 f |
|
f |
|
|
|
|
(41.9) |
(2 - |
л) Я - » + 2 |
[1 - |
m |
(я - |
1) Я~2 ] Н/ г _2 4- |
и; |
Д-» |
4- |
Яля _, /С„ |
= |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
= Л„ |
п > 1 , |
(41.10) |
//irt(« - fl)/? - » - 8 «„4 - »,[/C n + 2 +/C«Jc„ + 1 |
= |
- 2 i a „ + 1 > |
n > |
|
1, (41.11) |
|
|
|
|
а„ R-" |
- Кс~^Кп |
= |
Л_„, |
я > 1 . |
|
|
|
(41.12) |
п=\ |
При такой записи систем уравнений необходимо |
при |
|
первый |
член |
в |
левой |
части |
соотношения |
(41.10), |
содержащий |
множителем |
(и—2), считать равным |
нулю, |
а |
правые |
части |
Ai |
и Л_і равенств |
(41.10) |
и |
(41.12) —заменить |
соответственно |
на |
Ai— |
2TR |
и Л_, — Г ' R . |
|
|
|
|
|
|
|
|
|
ап, |
|
|
|
|
Система |
(41.11), (41.12) определит неизвестные |
с п + ь |
/ г ^ І , |
а система |
(41.10)—все |
неизвестные |
|
, начиная |
с |
/ г = 1 . Не |
достающие неизвестные а'0, с0 , ci определяются из формулы (41.9).
Решая совместно системы (41.11) и (41.12), находим
а,, =
(41.13)
д - (И-і) [2(-а„+ 1 / ? 3 - тп (я 4 1) Л _ „ ]
сп+\ =