книги / Проектирование и расчёт крепи капитальных выработок
..pdfПри этом вместо множителя а* используется величина [431
a * _ ( j _a *^ 0~1 |
1~Ь5/] aibt |
(6 . 12) |
|
|
|
О + Vo) (1 — 06,) + |
1.56, (г-и > )1_“' ’ |
|
где t l — время, проходящее между проведением выработки и воз ведением крепи, с; ос* определяется по формуле (6.2).
Для выполнения всех указанных выше расчетов используются решения плоской контактной задачи, изображенной на рис. 6.1, при граничных условиях, отражающих равенство векторов напря жений и смещений на линии контакта L, и отсутствие внешних сил на внутреннем контуре L x. После представления полных напря жений в виде сумм начальных напряжений, обусловленных собст венным весом или действием сил тектонического происхождения, выражаемых формулами:
а(х0>= — Nxa |
1 +2к* |
1 - X * cos 2a J , |
|
|
|
|
2 |
|
|
|
|
- X * |
cos 2aj , |
|
|
|
|
|
|
TxJ — — jVja* - 1 |
— sin 2a, |
(6.13) |
и дополнительных напряжений ax, au, xxy, вызываемых наличием выработки, и введения комплексных потенциалов срt (z), i|'f (г) (t = 0,1), регулярных в соответствующих областях S{ и связан ных с напряжениями известными формулами Колосова—Мусхе- лишвили, поставленная контактная задача для определения до полнительных напряжений сводится к краевой задаче теории ана литических функций комплексного переменного, решение которой получено в работе [41 ]. Полные напряжения определяются сумми рованием найденных дополнительных напряжений с начальными.
При расчете крепи на действие давления подземных вод началь ные напряжения в массиве определяются формулой
о40) = 4 0) = — ув {Нг— х), |
(6.14) |
вследствие чего граничные условия в контактной задаче для до полнительных напряжений отражают наличие на внешнем кон туре L сечения крепи скачка нормальных напряжений сгр, равного давлению воды
вр = ffp — ув (Hi— x), |
(6.15) |
где индексами к и м отмечены соответственно нормальные напряже ния в кольце (крепи) и в среде (массиве); ув — объемный вес воды; Н 1 — уровень воды в массиве, отсчитываемый от начала коорди нат (если крепь водонепроницаемая).
93
Рис. 6.2. Расчетные схемы для определения напряженного состояния крепи при действии длинной произвольно направленной волны:
а— продольной, б — поперечной
Вслучае водопроницаемой крепи вместо Нг в расчете исполь зуется уровень Н\, соответствующий остаточному напору, опреде
ляемый по формуле [27 ]
|
Н\ = |
(6.16) |
|
In |
|
где ku |
k 0— коэффициенты фильтрации |
соответственно материала |
крепи и породы; R и R 0 — наружный и внутренний средние радиусы |
||
крепи; |
Rt — условный радиус области |
питания [34]. Указанная |
контактная задача решена в работе [41 ].
Расчет крепи на сейсмические воздействия землетрясений со стоит в определении наиболее неблагоприятного напряженного со стояния в каждом сечении конструкции из возможных при различ ных сочетаниях совместного действия длинных продольных и попе речных волн любого направления в плоскости поперечного сечения выработки. С этой целью рассматриваются две плоские квазистатические контактные задачи теории упругости, расчетные схемы которых даны на рис. 6.2, а, б.
Напряженное состояние крепи от действия падающей под про извольным углом а к оси симметрии выработки длинной продоль ной волны а(р) определяется на основании решения контактной задачи для кольца произвольной формы (с одной осью симметрии), подкрепляющего отверстие в среде из другого материала, испыты-
вающей на бесконечности двухосное сжатие (см. рис. 6.2, а) напря жениями Р и IP, зависящими от параметров волн:
P = - j - l b y c 1T0, |
l = |
(6.17) |
где kc — коэффициент сейсмичности, соответствующий баллу земле трясения и принимаемый равным 0,1 — для интенсивности 9 бал лов, 0,05 — для 8 баллов, 0,025 — для 7 баллов (181; с , — ско рость распространения продольных волн, определяемая измере ниями или вычисляемая по формуле
с1= * ! — |
------- ^ |
------- , |
(6.18) |
|
V |
У |
(l + v „ ) ( t - 2 v 0) |
|
Т0 — преобладающий период колебаний частиц породы. Напряженное состояние крепи от действия длинной, падающей
под произвольным углом а поперечной волны a(S) определяется на основании решения плоской квазистатической задачи для кольца
в среде, |
подверженной на |
бесконечности чистому сдвигу |
(см. |
||||||
рис. |
6.2, |
б). |
|
|
|
на бесконечности |
|
||
Касательные напряжения |
|
||||||||
|
|
|
|
|
Q = |
—— |
|
|
(6.19) |
|
|
|
|
|
|
1я |
|
|
|
где |
с2 — скорость |
распространения |
упругих поперечных |
волн, |
|||||
определяемая из |
соотношения |
|
|
|
|||||
|
|
с2= |
л |
/ |
----- ^ |
---- = С х д / - |
1 --2Vfl |
(6.20) |
|
|
|
|
V |
2V (1 I- v0) |
V |
2 ( 1 - v o ) |
|
Для определения наиболее неблагоприятного напряженного состояния от сейсмических воздействий сумма и разность получаю щихся из решения указанных задач выражений для нормальных тангенциальных напряжений, отражающие напряженное состояние крепи от совместного действия одновременно приходящих (худ ший случай) продольных и поперечных волн
СГ01 = |
-f- СТо5^» |
002 — |
—OQS) |
(6.21) |
в каждом сечении конструкции исследуются на экстремум по углу падения волн а путем решения уравнений
- ^ “ - = 0 , |
-^ ® L = o, |
(6.22) |
д а |
д а |
|
и для каждого сечения находятся то сочетание действия волн раз ного характера и тот угол их падения, при которых нормальные тангенциальные напряжения в данном сечении максимальны по абсолютной величине, т. е. расчет производится на основании ана литического построения огибающих эпюр нормальных тангенци альных напряжений. В этом случае за расчетные на сжатие и на растяжение принимаются эпюры усилий, соответствующих напря
жениям а0, максимальным по абсолютной величине, взятые со зна
ками «+» и «—».
Если крепь проектируется с допущением образования трещин, то фаза растяжения в продольной волне из рассмотрения исклю чается, и расчет производится на основании двух огибающих, по строенных по максимальным значениям сжимающих и растяги вающих нормальных тангенциальных напряжений от совместного действия продольных волн в фазе сжатия и поперечных волн. В этом случае для проверки прочности сечений на сжатие и на растяжение используются по две различные эпюры расчетных уси лий — изгибающих моментов М и продольных сил ЛЛ Описанный подход, предложенный и подробно обоснованный в работе [43], принят в ВСН-193—81 Минтрансстроя [18]. Алгоритм расчета по лучен также в работе [43]. Ниже приводится полный алгоритм, составленный в результате обобщения и модификации алгоритмов, опубликованных в [41 и 43], и порядок расчета замкнутой мо нолитной бетонной и железобетонной крепи на статические на грузки — действие собственного веса пород или тектонических сил в массиве и давление подземных вод, а также на сейсмические воздействия землетрясений. Алгоритм реализован в виде программы FOK1 для ЕС ЭВМ (автор программы — А. Н. Козлов).
6.2. Алгоритм |
расчета |
|
Перечень входных и выходных параметров расчета приведен |
||
соответственно в табл. 6.1 и 6.2. |
|
|
Расчет состоит из трех этапов: |
1 — построение конформного |
|
отображения; 2 — определение корней основных систем |
уравне |
|
ний; 3 — определение расчетных напряжений и усилий |
в крепи. |
6.2.1.первый этап — построение конформного отображения —
заключается в определении коэффициентов а0, а ъ , а5, харак теризующих форму и размеры внутреннего контура поперечного
сечения |
крепи. Конформное отображение |
строится численным об |
|
разом |
на основе автоматизации |
графоаналитического способа |
|
П. В. |
Мелентьева. При этом |
форма |
и размеры внутреннего |
контура поперечного сечения крепи задаются согласно позиции 1.1 табл. 6.1. Для получения конформного отображения в программу FOK1 включена подпрограмма CONF, разработанная В. Л. Кипеневым и А. А. Лянда *
После |
нахождения |
коэффициентов 7i0, |
аи |
аъ решается |
|||
уравнение |
пятой степени |
вида |
|
|
|
||
|
а0хь-|- (ах—1± — А) х4+ а2х3+ а3х2+ а4х -|- аъ= 0, |
(6.23) |
|||||
где / х — длина |
первого луча; А — толщина крепи |
в своде. |
|
||||
Далее отыскиваются |
величины |
|
|
|
|||
K i = ~ - . |
R = a0R\, |
Ъ = а 1+11г\+'/а0 |
( i =l , |
,4), |
(6-24) |
№ п/п |
Параметр |
1Общие данные
1.1.Форма и размеры внутреннего контура по
перечного |
сечения крепи |
(задаются длинами |
13 лучей, |
проведенных из |
центра описанной |
окружности под углом 15° друг к другу). Номера лучей отсчитываются от вертикальной оси симметрии, м
1.2. |
Толщина крепи в своде, м |
|
1.3. |
Модуль деформации материала крепи, МПа |
|
1.4. |
Коэффициент Пуассона материала крепи |
|
1.5. |
Модуль деформации породы, МПа |
|
1.6. |
Коэффициент Пуассона породы |
|
1.7. |
Модуль деформации упрочненной или ослаблен |
|
|
ной зоны, МПа |
|
1.8. |
Коэффициент Пуассона упрочненной или ослаб |
|
|
ленной зоны |
|
1.9. |
Толщина упрочненной или ослабленной зоны, м |
|
1.10. |
Параметры линейной наследственной ползуче |
|
|
сти, с"41-0^ |
|
1.11. |
Время |
стабилизации деформаций, с |
1.12. |
Время, |
проходящее между проходкой выра |
|
ботки и возведением крепи, с |
|
2 |
Дополнительные данные в зависимости от вида |
|
|
нагрузки |
|
2.1. |
Для расчета на действие собственного веса по |
|
2.1.1. |
род |
|
Коэффициент бокового давления пород в нена |
||
|
рушенном массиве |
|
2.1.2. |
Объемный вес породы, кН/м3 |
|
2.1.3. |
Коэффициент пористости пород |
|
2.1.4. |
Глубина заложения выработки, м |
|
2.1.5. |
Расстояние от забоя выработки до места возве |
|
|
дения |
крепи, м |
2.1.6. |
Угол наклона пластов к вертикали, градус |
|
2.2. |
Для расчета крепи в массиве, подверженном |
|
|
действию тектонических сил |
|
2.2.1. |
Величина большего главного напряжения, МПа |
|
2.2.2. |
Отношение главных начальных напряжений |
|
2.2.3. |
Угол наклона главных осей начальных напря |
|
|
жений к вертикали и горизонтали, градус |
|
2.2.4. |
Расстояние от забоя выработки до места возве |
|
|
дения |
крепи, м |
Условные
обозначения
^1» ^21•••^13
д
Vl
Vo
Er
V'
h
6/, a t
t h
X
У
e
H
l
Ф
N i
X*
a
l
№ п/п |
Параметр |
Условные |
обозначения |
2.3.Для расчета на действие давления подземных
вод
2.3.1.Уровень подземных вод от центра выработки, м
2.3.2.Коэффициенты фильтрации материала крепи и породы, м/с
2.4.Для расчета на сейсмические воздействия
2.4.1.Коэффициент сейсмичности
2.4.2.Преобладающий период колебаний частиц по роды, с
2.4.3.Объемный вес, породы, кН/м3
2.4.4.Коэффициент пористости пород
и 1 ki, k0
кс
То
V е-
П р и м е ч а н и я : 1. Параметры 1,7— 1,9 задаются при наличии упрочненной или ослабленной зоны вокруг выработки. 2. Параметры 1.10— 1.12 задаются в случае, если породы подвержены ползучести. 3. Параметры 2.1.3, 2.4.4 задаются для расчета крепи в водоносных породах. 4. Параметры 2.3.2 задаются для учета водопроницаемости
крепи. 5. В качестве входного параметра |
для случая 2.4 может задаваться скорость |
|
распространения продольных волн |
м/с. 6. Если известна величина средней нагрузки |
|
на крепь Рср, то она может задаваться |
вместо параметров 2.1.2 — 2.1.5, 2.2.1, 2.2,4. |
7. Если неизвестен параметр 2.4.2, то может быть приближенно принято значение 7*0=0,5 с.
Таблица 6.2
№ п/п |
Параметр |
Условные |
обозначения |
1 Координаты точек внутреннего и наружного контуров поперечного сечения крепи и средней линии, м
2Нормальные контактные напряжения, МПа
3Касательные контактные напряжения, МПа
4Нормальные тангенциальные напряжения на внутреннем контуре поперечного сечения кре пи, МПа
5Нормальные тангенциальные напряжения на внешнем контуре поперечного сечения крепи, МПа
6Изгибающие моменты, МН-м
7Продольные силы, МН
*ъ Уъ
*2. У2*
*о. Уо
Тр0
ст0внутр
ОДвнеш
М
N
где в качестве R\ используется единственный действительный ко рень уравнения (6.23), больший единицы.
Описанный первый этап является общим для расчета крепи на все перечисленные виды нагрузок.
6.2.2. Второй этап расчета — определение корней основных систем уравнений — имеет общую часть для всех видов нагрузок, состоящую в вычислении коэффициентов матриц основных систем уравнений. Свободные же члены уравнений, составление допол нительной системы и формулы для выполнения третьего этапа рас-
чета (определение напряжений и усилий в крепи) различны для каждого вида нагрузки.
Второй этап расчета состоит из следующих операций.
6.2.2.1. Определяются |
величины |
|
|
Р = |
. |
*o = 3 - 4 v 0, |
*i = 3—4vx (6.25) |
Ео ( 1 + v j )
инаходятся входящие в расчетные соотношения комбинации де формационных характеристик материала крепи и окружающей
выработку горной породы по формулам:
t = - L+ У?Р -, |
|
d= _ L n i_ ( |
|
1= |
1— d, |
s = 1—1. |
||||||
1 + * 1 |
|
|
|
l+X l |
|
|
|
|
(6.26) |
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
6 .2.2.2. Вычисляются вспомогательные величины по рекуррент |
||||||||||||
ным формулам: |
|
|
|
|
|
|
|
|
|
|
||
|
|
3 -1 |
|
|
|
|
|
|
|
|
|
|
h i = q i + Si+i 2 |
|
k q k h k + i + i , |
|
|
|
|
|
|
||||
|
|
|
|
^ |
|
|
|
(i = |
4, |
3, . |
,0 ; |
?o= 0) |
hi = 4iRi l + 6i+i X |
hqkR\ ^kJrX)hk+i+1» |
|
|
|
|
|
||||||
|
|
|
|
k = \ |
|
|
|
|
|
|
|
|
4 |
|
|
|
|
|
|
|
|
|
|
|
|
h-t = b i + E |
|
fykh-t+ u |
|
|
|
|
|
|
(б-27) |
|||
k=\ |
|
|
|
|
|
|
|
|
( i = |
1,2,. |
8) |
|
л1,=х,^+ E « r (A+1)/iL1+b |
|
|
||||||||||
|
|
|
|
|
||||||||
|
ft=l |
|
|
|
|
|
|
|
|
|
|
|
где |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
J |
1 |
при |
i = l |
|
f l |
при |
i< 4 |
|
||
* |
|
j |
0 |
при |
i Ф 1 |
* |
( 0 |
при |
i > 4 |
|
||
и величины |
|
|
|
|
|
|
|
|
|
|
|
|
wt = Rih’- h [ |
(/------ 8 , |
—7, . |
|
|
|
0 ,1 ,. |
4).(6.29) |
|||||
6 .2.2.3. Составляется |
матрица системы |
I десяти линейных ал |
||||||||||
гебраических |
|
уравнений вида |
|
|
|
|
|
|
||||
4 |
|
|
|
|
6 |
= |
dk |
(/2 = |
1 , • |
|
4), |
|
V = 1 |
|
|
|
|
v = 1 |
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
||
4 |
|
|
|
|
6 |
|
|
|
|
|
|
|
Е ^.vCv+ |
Е a'k.vav = d'k |
(ft = |
1, • |
|
6). |
(6.30) |
||||||
v = l |
|
|
|
|
V = 1 |
|
|
|
|
|
|
|
* В Тульском политехническом институте имеется также программа построения конформного отображения по координатам произвольного числа точек контура, разработанная В. Г. Гараичуком.
Коэффициенты при |
неизвестных |
cv |
(v = 1, |
4 ) |
в |
первых |
|
четырех уравнениях определяются по формуле |
|
|
|
||||
Ck, v = ^k, v {dR2\k+ 0 |
+ |
^dw—v—А+i» |
|
|
(6.31) |
||
где |
1 |
при |
k = V, |
|
|
|
|
|
|
|
(6.32) |
||||
кь.ч,= |
ПрИ |
к ф \ . |
|
|
|||
4 ,v — | Q |
|
|
|
||||
Коэффициенты при |
неизвестных |
av |
(v = 1, |
6 ) |
в |
первых |
|
четырех уравнениях определяются |
из |
соотношения |
|
|
|
||
|
|
|
2 |
|
|
|
|
&k, v— ^ ^6V_ |
k+i |
d 5^ 6v-\.ppWp—k+ihv+p+i^ • |
(6.33) |
Коэффициенты при неизвестных cy (v = 1, шести уравнениях вычисляются по формуле
Ck, v = dv^k-v^k-v+l-
Коэффициенты при неизвестных av (v = 1, шести уравнениях определяются по формуле
, 4) в остальных
(6.34)
, 6) в последних
= |
v ( ^ l 2^ + |
s) — v j^v+fc(cf/?i 2k + |
l) hv+k+l + |
t$v+kw v+k+i — |
||
|
|
d&k+i |
^v+pPWp+k+ihv+p+iJ • |
|
(6.35) |
|
6 .2.2.4. |
Вычисляются |
свободные |
члены системы уравнений 1 |
|||
в следующем порядке: |
на действие |
собственного |
веса пород |
|||
6.2.2.4.1.— для |
расчета |
|||||
1 Л 1 -J- A, T-J-2 |
1 -- ^ |
( * = 1, |
, |
4), |
||
dk = h — т— t f i ---------— % |
||||||
^ = 6*-!— |
|
|
( k = l , |
|
, 6 ). (6.36) |
|
6 .2.2.4.2.— для |
расчета |
на действие |
внешнего |
гидростатиче |
ского давления подземных вод вычисляются вспомогательные ве личины по формулам
Я2 = |
4 |
|
(6.37) |
Нъ = /?2 + Z |
V> |
||
|
v=1 |
|
|
где /?2 — единственный больший единицы действительный корень уравнения пятой степени (6.23), решенного при подстановке вместо li + Д в коэффициент при х* значения Н г (или Н\, вычисленного по формуле (6.16), если крепь водопроницаема).
Определяются коэффициенты:
cc1 = 2H 2 + q1q2 + 2q2q3+ 3q8qi ; |
а 2= ------ |
(1 + q i ~ q 1q 3— 2<72<7,); |
аз = |
----- £ - (?« — ftft); |
|
|
а 4 = ----- --- q3; |
а 5 = ----- ^ -q 4; |
|||||||||
Ро= — Яй |
|
$l = |
2Hiq1— q2— 2q1q2— 3q2q3— 3qiq3‘lq3qi ; (6.38) |
|||||||||||
Рг = |
2НгЦг — <7з------ (<7 i + |
Ц\ лг 3q4qz + |
^ЯгЯЯз + ^ЯгЯ*)' |
|||||||||||
|
Рз = |
|
|
|
|
|
|
|
2 |
|
4 |
|
|
4 |
|
2Н 2<7з — 94— Я1Я2 -------- Яг------— ЯхЯЯг-----“ ЯхЯь\ |
|||||||||||||
|
|
|
Р4 = |
2Н 2^ 4 ----- — Я2 — Я1Я3 ---------- Яг\ |
|
|||||||||
|
|
|
|
|
|
|
|
|
2 |
|
|
4 |
|
|
|
|
|
|
Р б — ---- ? 1 ? 4 -----? 2 ? 3 ----------~ |
^ |
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
О |
*— И |
||
Рв = |
— ЯгЯ4 ------ —Х |
& |
|
|
р7 = — w 4; |
|
||||||||
|
|
|
|
|
||||||||||
к свободные члены системы уравнений 1 по формулам: |
||||||||||||||
dx= ------ - — |
d Y |
v/iv+1oyv-------- 1— |
Y vPva\'H— |
d*1 — a i’ |
||||||||||
|
1 + H„ |
v=o |
v |
|
|
l + x x |
v=o |
|
|
</(l+x,) |
||||
do = |
|
F |
j |
v ' |
, |
, |
„ , |
|
|
1 |
4 |
(v+l)X |
||
|
|
d |
E |
(v+l)/jv+2oyv- |
. , „ |
E |
||||||||
|
|
1 -h XQ |
v=0 |
|
|
|
|
|
1-j- Xi v=o |
|
||||
|
|
|
|
|
X |
Pv+l^v H---- -- ------ - |
a 2, |
|
(6.39) |
|||||
|
|
|
|
|
|
|
|
|
|
d{ t + |
xx) |
|
|
|
|
|
do — |
|
|
d [ E |
(v + 2)fiy+swv+ M - i j - |
||||||||
|
|
|
|
l+xo |
|
|
|
|
|
|
|
|
||
|
|
1 |
[ |
4 |
|
|
|
|
|
|
|
. |
dxi — l |
|
|
1 + Xi |
E o (v + |
2) Pv+sO'v + P i^ - i] H-----------------«з; |
|||||||||||
|
|
|
|
|
|
|
|
|
|
|
d( 1+xO |
|||
|
|
d4 = |
1 + |
|
■d [З/14КУ0 “f* W—2^2 "Ь 2ш_!/г3] — |
|||||||||
|
|
|
|
x0 |
|
|
|
|
|
|
|
|||
|
|
1 |
|
|
|
|
|
|
|
|
|
|
d*! — l |
|
|
l+X ! |
^ E Q(V + 3) Pv+3^V+ Pl^-2 j |
d{ 1+ x,) ■«4; |
|||||||||||
|
|
|
F |
|
|
2 |
|
|
|
|
. |
2 |
vpvnyv+2 + |
|
|
|
|
|
d E |
vftv+1ayv+2----г т — E |
|||||||||
|
|
1 + x 0 |
v= l |
|
v_r* |
1 + Xi v~ l |
|
|
||||||
|
|
|
1 |
о |
/n - 2 |
, |
\ |
, |
/=■ |
л2 И г 2 + 0 = |
||||
|
|
|
|
Pi (#1 2 + |
xi) |
+ |
|
|||||||
|
|
1 -I- X i |
|
|
|
|
|
1 + X0 |
||||||
do = |
F |
|
dh2w4i- |
|
|
1 |
Pi^V |
|
--- |
p2 (^?l 4 + Xl) 4" |
||||
1 + x0 |
|
1 + X i |
|
|||||||||||
|
|
|
|
|
|
|
1 + X i |
|
|
|||||
|
|
|
|
|
|
|
|
F |
Л |
/ J n - 4 |
/), |
|
|
|
|
|
|
|
|
|
|
|
|
з ( « Г 4 + |
|
|
1 + X0
ds = |
1 |
|
|
|
|
|
К (dRT6 + l). |
||
1 +XX Рз (R I |
6 |
+ Xl) + |
|||||||
|
1 + Щ |
|
|||||||
йл = |
1 |
(-^1 8 H- Xi), |
^5 — |
Ps ( R T 10+ Xi)> |
|||||
P4 |
|||||||||
1 + |
Til |
|
|
|
|
|
1 + |
*1 |
|
|
|
da = |
1 |
P.(/?r,2 + >*i). |
|||||
|
|
|
|
||||||
где |
|
1 + |
|
*1 |
|
|
|
||
|
|
|
|
|
|
|
|
||
|
|
f = |
|
l - Z |
v |
9- |
(6.40) |
||
|
|
|
|
|
V = 1 |
|
|
||
6 .2 .2 .4.2 .— для |
расчета |
крепи |
в массиве |
пород, подверженном |
действию тектонических сил, или в массиве с мелкослоистым на клонным напластованием (главные оси начальных напряжений наклонены к вертикали и горизонтали под углом а) система урав нений I решается дважды с группами свободных членов (1а) и (16), которые вычисляются по формулам
dk (la) = XkRi |
4) |
|
( * = 1, . |
|
|
dk (16) = — <7* |
|
|
d'k (la) = |
6) |
(6.41) |
( k = \ , |
d'k{Щ = - К
В результате решения системы уравнений I получаются две группы корней:
cv(I, a) — v = |
1 , . |
4 , |
av(I, a)—v = |
1, . |
6 , |
cv(I, 6)— v = |
1 , . |
4 , |
av(I, 6)—v = |
1 , |
, 6 . |
Дополнительно составляется система уравнений II, матрица ко торой получается из коэффициентов матрицы системы I по фор мулам:
ст, v — с т, V — |
2 ( / - f d/?im) |
при т = v (т= 1, . |
4), |
||
I I |
|
|
|
|
|
яы = —01,1 — 2tw11 |
0 i , 2 = |
— 0 i , 2 — |
4 t o 2, |
|
|
02.1 = |
-- 02,1---2tWo, |
02,2 = |
02,2 |
4Й01, |
|
0 зл = —03.1— 2 to_b |
0 з ! г = |
— 03.2 — |
4 to o , |
|
|
04.1 == |
--- 0 4 ,1 -----2 t W —2> |
0 4 , 2 = |
04,2 |
4 / Ш _ 1 , |
|
01.1 = |
—0 ы —2 £to3 Н- Лз 0 “h |
**)]» |
|
|||
01.2 = |
—0 i,2—4 [to4 |
Л4 (/ -f- dR\ |
)], |
(6.42) |
||
02,i = |
—02,i —2 [to4 + /i4 (/ + d/?i |
)]> |
|
|||
'11 |
----02,2, |
|
= |
— 0 3 .3, |
|
|
02,2 = |
03,3 |
|
|
|||
'II |
--- 0 4 ,4, |
'II |
= |
— 05 ,5, |
|
|
04,4 = |
05,5 |
|
|
'II
06,6 = --- 00,6,