Теоретические основы работы газостатических опор. Монография
.pdf
|
12αμasπdH |
|
2 |
|
κ +1 |
||
γ = |
( |
) |
2 (κ −1) |
E. |
|||
|
|
||||||
|
π H 3 ps |
κ + 1 |
|
|
|
||
Параметры γ, γ è γ называют конструктивными параметрами питателей.
Åñëè x = xd – решение уравнения (3.36) или уравнений (3.38) (в зависимости от типа питателя), то давление на кромке питателя
определяется по формуле pd |
= ps |
xd . При этом распределение дав- |
||||||
лений в смазочном слое находится на основе формулы (3.31): |
||||||||
|
|
p 2 |
− p 2 |
|
R |
|
||
p = |
|
d |
|
a |
ln |
o |
+ ps 2 . |
(3.39) |
|
|
E |
|
|
||||
|
|
|
|
|
r |
|
||
Несущая способность рассматриваемой дисковой опоры вычис- |
||||||||
ляется с помощью ЭВМ по формуле |
|
|||||||
|
|
|
R |
|
|
|
|
|
|
F = 2π ∫O |
( p − pa )rdr. |
|
|||||
|
|
|
0 |
|
|
|
|
|
Коэффициент при lnr в формуле (3.39) определяет коэффициент A в формуле (3.31):
A = − pd 2 − pa 2 .
E
Поэтому расход газа для дискового подпятника можно определить по формуле
Q = |
κπ H 3 |
|
p 2 |
− p |
2 |
. |
|
|
|
|
d |
a |
|||
|
2 |
|
|
|
|||
|
12μas |
|
|
E |
|
|
|
Расчет полосовой опоры с щелевым дросселем
Рассмотрим (рис. 3.3) прямоугольную полосовую опору протяженности T и ширины L, причем L << T. Пусть газ в смазочный слой подается через продольную капиллярную щель, имеющую ширину b и такую же протяженность, что и опора. Давление на входе в щелевой дроссель принимаем равным ps, а давление окружающей среды – pa. Принимаем толщины смазочного слоя H и дросселирующей щели c
61
постоянными. Для упрощения исследования пренебрегаем концевыми эффектами в опоре и на входе в щелевой дроссель.
|
|
Х |
У |
а |
б |
|
Ps 
Ðèñ. 3.3. Схема щелевой опоры
Åñëè îñü Ox направить поперек питающей щели, то комплексные потенциалы потоков в зазоре определятся, очевидно, по формуле (3.28), из которой следует, что для смазочного слоя в опоре справедливы равенства
|
|
p 2 |
|
− p 2 |
|
X + pd 2 , −l1 < X < 0, |
|||||||
|
|
|
|
|
d |
|
|
a |
|
||||
|
|
|
|
|
|
l1 |
|
|
|||||
p |
2 |
|
|
|
|
|
|
|
|
|
|
(3.40) |
|
|
= |
|
|
pd |
2 |
− pa |
2 |
|
|
||||
|
|
|
|
|
|
2 |
|
||||||
|
|
|
− |
|
|
|
|
|
|
X + pd |
, l2 > X > 0, |
||
|
|
|
|
|
l2 |
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
− p |
2)Y , −l < X < 0, |
||||
|
|
|
|
|
|
( p 2 |
|||||||
|
|
ψ = |
|
d |
|
|
a |
l |
1 |
||||
|
|
|
|
|
|
|
1 |
|
|||||
|
|
|
|
|
|
−( p 2− p 2) Y , l > X > 0. |
|||||||
|
|
|
|
|
|
|
d |
|
|
|
a l |
2 |
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
2 |
|
|
Масса газа, протекающего в единицу времени из участка питающей щели шириной dY в область смазочного слоя X < 0, согласно формулам (3.24) и (3.38), равна
dQ = |
κ H 3 |
|
( p 2 |
− p 2 ) |
dY |
, |
(3.41) |
24μa |
|
|
|||||
1 |
2 |
d |
a |
l |
|
||
|
s |
|
|
1 |
|
|
|
62
а в область X > 0
dQ |
|
= |
κ H 3 |
|
( p |
2 − p 2 ) |
dY |
. |
(3.42) |
|
24μas |
|
|
||||||
|
2 |
|
2 |
d |
a |
l2 |
|
||
Так как поток газа в дросселирующей щели полностью аналогичен потоку газа в смазочном слое, но происходит под действием давлений ps è pd, то нетрудно записать, что через участок щели протяженности dX в единицу времени протекает масса газа, равная
dM = |
κ H 3 |
|
( p 2 |
− p 2 ) |
dY |
. |
(3.43) |
24μas |
|
|
|||||
|
2 |
s |
d |
b |
|
||
На каждом элементе границы смазочный слой – питающая щель должно выполняться условие баланса расходов газа
dQ1 + dQ2 = dM,
из которого, учитывая формулы (3.39)–(3.41), следует уравнение для определения давления на выходе из щели
( |
1 |
+ |
1 |
)H 3 ( pd |
2 |
− pa |
2 ) = |
c3 |
( ps |
2 − pd |
2 ). |
(3.44) |
|
|
|
b |
|||||||||
|
l1 l2 |
|
|
|
|
|
|
|
||||
Легко находим
|
|
|
|
|
|
pd = ps |
γ + p |
2 |
, |
|
|
|
|
|
|
1 |
|
||
|
|
|
|
|
|
γ + 1 |
|
||
|
|
|
|
|
|
|
|
|
|
ãäå γ = |
|
l l |
c3 |
, p1 |
= pa / ps . |
|
|
|
|
|
1 |
2 |
|
|
|
|
|||
(l |
+ l |
)bH 3 |
|
|
|
||||
1 |
|
2 |
|
|
|
|
|
|
|
Несущая способность части полосовой опоры, лежащей в области X < 0, вычисляется по формуле
0 |
0 |
2 |
2 |
|
|
F1 = T ∫ pdX = T ∫ |
pd |
− pa |
X + pd |
2 dX . |
|
|
|
||||
− l1 |
− l1 |
|
l1 |
|
|
Интеграл легко вычисляeтся. Получаем
F1 = |
2 |
Tl1 |
( pd |
+ |
pa |
2 |
). |
3 |
|
|
|||||
|
|
|
|
pd + pa |
|||
63
Аналогично запишется формула для несущей способности части полосовой опоры, лежащей в области X > 0:
F2 = |
2 |
Tl2 |
( pd |
+ |
pa |
2 |
). |
3 |
|
|
|||||
|
|
|
|
pd + pa |
|||
Очевидно, что полная несущая способность полосовой опоры с щелевым поддувом равна
F = |
2 |
TL( pd |
+ |
pa |
2 |
). |
(3.45) |
|
|
|
|||||
|
3 |
|
|
pd + pa |
|
||
Массу газа, протекающего в единицу времени через смазочный слой полосовой опоры, можно определить путем суммирования по всей длине питающей щели элементарных расходов dM, вычисляемых по формуле (3.41), что дает
Q = M = |
κ c3 |
( ps 2 |
− pd |
2 ) |
T. |
(3.46) |
|
24μas |
2b |
|
|||||
|
|
|
|
||||
3.4. Метод источников и стоков
Если поддув газа в смазочный слой плоской опоры производится через питатели диафрагменного типа, то в ряде практически важных случаев удается построить комплексный потенциал потока в смазочном слое, используя метод источников и стоков. Для этого рассматривают задачу о течении газа в слое той же самой толщины, что и в рассчитываемой опоре, но больших, иногда неограниченных, размеров, чем смазочная пленка. Распределяют по всему этому слою системы источников и стоков таким образом, чтобы порожденный ими поток газа в области, соответствующей смазочному слою, полностью совпадал с потоком в этом слое. Комплексный потенциал системы источников и стоков определяется в виде суммы комплексных потенциалов источников, стоков и поступательного потока. Достаточно просто строится комплексный потенциал для бесконечно длинной полосовой опоры с одной дорожкой питателей малого диаметра, параллельной границам опоры. Он лежит в основе теории смазки газостатических подшипников различных типов, находящих применение в технике. Поэтому остановимся более подробно на выводе соответствующих формул.
64
Комплексный потенциал дорожки источников в плоскости
Используем декартовы координаты ξ, η. Îñü Oξ направляем вдоль границы полосы (рис. 3.4). Пусть ζ = ξ + iη. Рассматриваем случай однорядной опоры ширины L, когда все питатели одинаковые, имеют диаметры d и расположены равномерно с шагом B на линии η = l, так что их центры имеют координаты
ζk = ζ0 + Bk, k = 0, ±1, ±2,.. ,
ãäå ζ0 – координата центра питателя, выбираемого в качестве базового ζ0 = ξ0 + il.
η
2 /n
λ 2
L
|
1 |
|
λ |
O |
ξ |
Ðèñ. 3.4. Схема полосовой однорядной опоры
Построим комплексный потенциал, порожденный рассматриваемой дорожкой питателей в безграничном слое, имитируя каждый питатель точечным источником мощности Q′. Для источника, помещенного в точку ς = ςk, имеем (см. формулы (3.32) и (3.48))
w = − |
12μa |
2 |
Q′ ln(ζ − ζ |
|
) + c |
|
= |
|
||||||||||||
|
|
|
s |
|
|
|
|
|||||||||||||
|
κπ H 3 |
|
|
|
||||||||||||||||
|
k |
|
|
|
|
|
|
|
k |
|
|
k |
|
|
|
|||||
= − |
12μa 2 |
Q′ ln(ζ − ζ |
|
− Bk) + c |
|
= |
||||||||||||||
|
|
|
s |
|
|
|||||||||||||||
κπ H 3 |
0 |
k |
||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
= − |
12μa |
2 |
Q′ ln(1 − |
ζ − ζ |
k |
) + c |
, (k ≠ 0). |
|||||||||||||
|
|
s |
|
|
|
|||||||||||||||
κπ H |
|
|
Bk |
|
||||||||||||||||
|
3 |
|
|
|
|
|
|
|
k |
|
|
|
|
|
||||||
65
Найдем комплексный потенциал всей системы источников как сумму комплексных потенциалов каждого из них:
|
|
12μas |
′ |
|
|
|
|
|
|
|
|
|
ζ |
− ζ0 |
|
ζ − ζ0 |
|
||
|
|
2 |
|
|
|
|
|
|
∞ |
|
|
|
|
|
|
|
|
|
|
w = − |
|
|
Q {[ln(ζ − ζ |
0 |
) + |
∑ln[(1 − |
|
|
|
)(1 + |
|
)]} +c = |
|||||||
κπ H 3 |
|
|
Bk |
Bk |
|||||||||||||||
|
|
|
|
|
|
k =1 |
|
|
|
|
|
|
|
||||||
|
|
2 |
|
|
|
|
∞ |
|
|
|
ζ − ζ0 |
|
|
|
|
|
|
||
= − |
12μas |
Q′ ln{(ζ − ζ |
|
)∏(1 − ( |
)2 ]} + c. |
|
|
||||||||||||
|
0 |
|
|
|
|||||||||||||||
|
κπ H 3 |
|
k =1 |
|
|
|
|
Bk |
|
|
|
|
|
|
|||||
Так как имеет место равенство |
|
|
|
|
|
|
|||||||||||||
|
|
|
|
|
|
∞ |
|
|
|
β |
2 |
|
|
|
|
|
|
|
|
|
|
|
|
β∏(1 − |
|
|
|
|
) = sin β, |
|
|
|
|||||||
|
|
|
|
|
k 2π 2 |
|
|
|
|||||||||||
|
|
|
|
|
k =1 |
|
|
|
|
|
|
|
|
|
|
||||
то можно записать, что комплексный потенциал, порождаемый дорожкой источников в безграничном слое, определяется по формуле
w = − |
12μa |
2 |
Q′ ln sin |
π (ζ − ζ |
) |
+ D, |
s |
|
0 |
|
|||
κπ H 3 |
B |
|
||||
|
|
|
|
|||
ãäå D = D0 + iD* – комплексная постоянная.
Комплексный потенциал дорожки источников в полосе
Поток газа, порождаемый дорожкой питателей в полосе при условии, что на границах полосы давления постоянные и равны pa, может быть получен от распределенных по всей безграничной плоскости систем дорожек источников и стоков обильности Q, параллельных границам полосы, базовые источники которых
имеют координаты |
ζ m = ζ0 + 2mLi, (m = ± 1, ± 2, ...), а базовые сто- |
||||
êè – ζ m′ = |
|
|
(m′ = ± 1, ± 2, ...; |
|
|
ζ0 + 2m′Li, |
ζ0 = ξ0 − iη0 ) , и некоторого |
||||
поступательного потока, нормального к границам полосы. Здесь при построении искомого комплексного потенциала нельзя непосредственно суммировать все бесконечное множество комплексных потенциалов, так как такая сумма будет расходящейся. Поэтому рассмотрим сначала случай, когда в плоскости расположено ограниченное число (2K + 1) дорожек источников и столько же дорожек стоков.
Суммарный потенциал, порождаемый дорожками источников, равен
66
|
12μas |
′ |
π (ζ − ζ0 ) |
|
π (ζ − ζ − k ) |
|
π (ζ − ζ k ) |
|
|
2 |
|
|
K |
|
|
|
|
u = − |
|
Q {[ln sin |
|
+ ∑ln[sin |
|
sin |
|
]}. |
κπ H 3 |
|
|
|
|||||
1 |
|
B |
k =1 |
B |
|
B |
|
Используя формулу для произведения синусов, можно записать с точностью до постоянной
|
|
12μas |
2 |
|
π (ζ − ζ0 ) |
K |
cos |
2π (ζ − ζ0 ) |
|
||
|
|
′ |
|
B |
|
||||||
u |
= − |
|
|
Q {[ln sin |
|
+ ∑ln[1 − |
|
|
|
|
]}. |
κπ H 3 |
B |
|
|
4kL |
|||||||
1 |
|
|
k =1 |
|
ch |
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
B
Аналогично, суммарный потенциал, поставляемый дорожками стоков, равен
|
|
|
12μas |
2 |
|
π (ζ − |
|
|
|
cos |
2π (ζ − ζ 0 ) |
|
|||
|
|
|
′ |
ζ0 ) |
K |
|
B |
|
|||||||
u |
|
= |
|
|
|
Q {[ln sin |
|
|
|
+ ∑ln[1 − |
|
|
|
|
]}. |
2 |
|
κπ H 3 |
B |
|
ch |
4kL |
|||||||||
|
|
|
|
k =1 |
|
|
|||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
B
Введем еще в рассмотрение однородный поступательный поток, направленный поперек полосы, комплексный потенциал которого для удобства представляем в виде
u |
|
= |
12μa 2 |
Q′b i |
ζ |
+ e , |
(3.47) |
|
|
|
s |
|
|||||
|
|
κπ H 3 |
L |
|||||
|
3 |
|
|
k |
k |
|
||
ãäå bk è ek – постоянные, подлежащие определению. Рассмотрим теперь поток с суммарным потенциалом
|
|
|
|
|
|
|
|
|
|
uK = u1 + u2 + u3 = |
|
|
|
|
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
π (ζ − |
|
|
|
|
|
4kL |
|
2π (ζ − |
|
|
|
|
|
||
|
12μas |
2 |
|
|
|
ς |
|
sin |
ζ0 ) |
|
K |
ch |
− cos |
ζ0 ) |
|
|
|
||||||||
|
′ |
|
|
|
B |
|
|
|
B |
|
B |
|
|
|
|
|
|||||||||
= |
|
|
|
Q {b i |
|
+ ln |
|
|
|
|
|
+ ∑ln[ |
|
|
|
|
|
|
|
|
]} + e |
|
. |
||
|
κπ H 3 |
|
|
π (ζ − ζ0 ) |
|
|
4kL |
|
2π (ζ − ζ |
0 ) |
|
||||||||||||||
|
|
|
K |
|
L |
|
sin |
|
k =1 |
ch |
− cos |
|
K |
|
|||||||||||
|
|
|
|
|
|
|
|
|
B |
|
|
|
|
B |
|
B |
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
Определим постоянные bK è eK, выполняя следующие условия: |
|||||||||||||||||||||||
|
|
|
|
|
|
Re[uK (ξ0 )] = pa |
2 , |
|
Re[uK (ξ0 |
+ iL)] = pa |
2 . |
|
|
|
|
|
|
||||||||
67
Используя формулы sin (ia) = isha è cos (ia) = cha, находим
|
sh |
π |
[(2k + 1)L + l]sh |
π |
[(2k − 1)L − l] |
|
sh |
π |
[L + l] |
|||||
K |
B |
B |
|
B |
||||||||||
uk = ln ∏ |
|
|
|
|
|
|
|
+ ln |
|
|
= |
|||
|
π |
|
|
|
|
π |
|
|
|
π |
|
|||
k =1 |
sh |
[(2k + 1)L − l]sh |
[(2k − 1)L |
+ l] |
|
sh |
[L − l] |
|||||||
|
|
|
|
|||||||||||
B |
B |
|
B |
|||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
sh |
π |
[(2K + 1)L + l] |
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
= ln |
|
B |
|
|
. |
|
|
|
|
|
|
|
|
|
sh |
π |
[(2K + 1)L − l] |
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|||||
B
Искомый потенциал, отвечающий потоку газа, порожденному дорожкой питателей в полосе, получается предельным переходом при K → ∞. В этом случае
|
|
|
|
|
lim bK |
= |
2πl |
, |
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
B |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
w = lim u |
|
= |
12μa |
2 |
Q′f (ζ ), |
|
|
|
|
|
|
||||||||||||||
|
|
|
|
|
|
|
|
s |
|
|
|
|
|
|
|
|||||||||||||
|
K |
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||
|
|
|
|
|
|
|
|
|
|
κπ H 3 |
|
|
|
|
|
|
|
|
|
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
π (ζ − |
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
sin |
ζ |
0 ) |
|
|
2πl |
|
|
|
|
||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
B |
|
|
|
|
|
|
|
2 |
|
|||||||
ãäå |
f (ζ ) = v(ζ ,ζ0 ,ζ0 ) + ln |
|
|
|
|
|
|
|
|
|
+ |
iζ + pa |
, |
|||||||||||||||
sin |
|
|
π (ζ − ζ |
0 ) |
B |
|
||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
B |
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
1 − cos |
2π (ζ − |
|
|
|
4kL |
|
|
|
||||||||||||||||
|
|
|
∞ |
ζ0 ) |
/ ch |
|
|
|
||||||||||||||||||||
|
|
|
|
|
|
|
|
|
||||||||||||||||||||
|
v(ζ ,ζ0 , |
ζ0 ) = ∑ln[ |
|
|
|
|
|
|
|
|
|
B |
|
|
|
|
|
|
|
B |
]}. |
|
|
|||||
|
|
|
|
|
|
|
|
2π (ζ − ζ0 ) |
|
|
4kL |
|
|
|||||||||||||||
|
|
|
k =1 |
1 − cos |
/ ch |
|
|
|
||||||||||||||||||||
|
|
|
|
|
|
|
|
|||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
B |
|
|
|
|
|
|
|
B |
|
|
|
||
Отметим два свойства аналитической функции важных для дальнейших исследований.
(3.48)
v(ζ ,ζ0 ,ζ0 ) ,
68
1. На границе полосы η = 0 реальная часть функции v равна нулю, что можно записать в виде
Re[v(ζ ,ζ0 , |
|
|
|
ζ0 )] = 0. |
(3.49) |
||
В справедливости этого свойства легко убедиться простой проверкой, учитывая равенство
1
Re[v(ξ,ζ0 ,ζ0 )] = [v(ξ,ζ0 ,ζ0 ) + v(ξ,ζ0 ,ζ0 )]. 2
Из вида функции v (3.48) следует, что
v(ξ,ζ0 ,ζ0 ) = −v(ξ,ζ0 ,ζ0 ),
что и доказывает справедливость утверждения (3.49). 2. Имеет место равенство
B / 2+ξo |
|
Re ∫ v(ζ ,ζ0 ,ζ0 ) η= const dζ = 0. |
(3.50) |
− B / 2+ξo
Для доказательства этого свойства рассмотрим интеграл
Коши от функции v(ξ, ς 0 , ς0 ) по прямоугольному контуру ABCD, ограниченному прямыми η = 0, (AB); ξ = B/2 + ξ0, (BC); η = const, (CD); ξ = – B / 2 + ξ0, (DA). Реальная часть этого интеграла равна нулю в силу аналитичности неоднозначной логарифмической подинтегральной функции
Re ∫v(ζ ,ζ0 ,ζ0 )dζ = 0.
С
Можем записать
Re[ ∫v(ς , ς 0 , ς0 )dς + ∫v(ς , ς 0 , ς0 )dς + ∫v(ς , ς 0 , ς0 )dς + ∫v(ς , ς 0 , ς0 )dς ] = 0.
AB BC CD DA
В силу периодичности функции v ïî ξ имеем
∫v(ς , ς 0 , ς0 )dς + ∫v(ς , ς 0 , ς0 )dς = 0.
BC DA
69
Âсилу первого свойства подынтегральная функция в интеграле
∫v(ς , ς 0 , ς0 )dς равна нулю, поэтому и сам этот интеграл равен нулю.
AB
Следовательно, имеет место равенство
Re ∫ v(ς , ς 0 , ς0 )dς = 0.
CD
η = const на отрезке CD, а концы этого отрезка совпадают с пределами интегрирования в формуле (3.50), справедливость которой таким образом доказана.
3.5.Расчет полосовых однорядных газостатических опор
ñдискретным поддувом
Рассмотрим сперва однорядную полосовую опору протяженности T >> L, вдув газа в смазочный слой которой производится через n одинаковых питателей диафрагменного типа, расположенных равномерно по всей длине опоры в один ряд, отстоящий от боковой границы полосы на расстоянии L (ñì. ðèñ. 3.4).
При расчете такой опоры будем пренебрегать концевыми эффектами и рассматривать ее как часть бесконечно длинной полосовой опоры с одним рядом питателей, имеющем шаг B = T/n. В этом случае при H = const для комплексного потенциала потока в смазочном слое применима формула (3.48), а среднеквадратичое давление на кромке питателя может быть вычислено по формуле
p d 2 = |
1 |
Re |
|
∫ w(ζ ) |
|
dζ |
|
. |
(3.51) |
|||
|
|
|||||||||||
|
||||||||||||
|
πd |
|
ζ −ζO |
|
= d / 2 |
|
||||||
|
|
|
|
|||||||||
|
|
|
|
|
|
|||||||
Здесь контурные интегралы от аналитических в круге фун- |
||||||||||||
кций, входящих в w(ζ), вычисляем по теореме о среднем, а от функции, имеющей логарифмическую особенность в центре круга, – непосредственно. При вычислениях учитываем, что на кромке питателя
ζ − ζ0 |
= |
1 |
deiθ , |
|
dζ |
|
= |
1 |
|
|
dθ |
|
, |
|
|
|
|
|
|
||||||||||
|
|
|
||||||||||||
|
2 |
|
2 |
|
|
|
|
|
, |
|||||
ln sin π (ζ – ζ0 ) ≈ ln πd + iθ.
B 2B
70
