РАСЧЕТ СТАТИЧЕСКИХ ХАРАКТЕРИСТИК ПОДШИПНИКОВ ТУРБОАЕРЕГАТОВ БОЛЬШОЙ МОЩНОСТИ

РТМ 108.020.101-76

Издание официальное

РАЗРАБОТАН И ВНЕСЕН Центральным научно-исследовательским и проектно-конструкторским котлотурбинным институтом им. И. И. Ползунова

Директор

Заведующий базовым отраслевым отделом стандартизации Заведующий отделом паровых турбин Руководители темы:

Исполнители:

ПОДГОТОВЛЕН К УТВЕРЖДЕНИЮ Техническим управлением Министерства энергетического машиностроения

Начальник Технического управления    В.    П.    ПЛАСТОВ

Начальник отдела опытно-конструкторских и научно-исследовательских работ

по турбостроению    В.    К.    ЛЕОНТЬЕВ

УТВЕРЖДЕН И ВВЕДЕН в действие указанием Министерства энергетического машиностроения от 1 апреля 1976 г. № ПС-002/2421

Заместитель министра

PTM 108.020.101—76 Стр. 9

Н4;*| J+*K*60M-¥ti»i*|||(|)tJX*etlO ;

nv:«mf*i i

hi:=m i+xmicotfiNvi'tfVKaiNSiNvi)t3 ;

FOR K;*1 *TEp 1 UNTIL N do В|в|ц R(KI:«Hi i C(K);*H2 ;

F(Nlt«H9 I F(«+M):»H« :

FIN+NEllaMF |

Btkku-n+k1;жн    i

IF N»l TH*N $I Kxtt-Na+K];S-H4 ;

If N«N THEN BIKxn+K):*-H4 i END |

b(bn>h)i»a*BiNN-Hi : boo*iu,*.*) ;

FOR K:«1 ITIP 1 UMTlt NN £2 BtK)|kB{N|.K|R| {

FOOB7(R,B1,Ra) $

FOR 3 i «I STEP I U»T£L MM 22 *!«:«• i

FOR 3;«l «Tip 1 UNT U M 00 XOxn>n*J);*C{ J) ;

PlOBTlN.M,PI) |

FOB 3:*| 3TCP I UNTIL М3 22 Ft3);sFlJl+PllJ) ;

FIOCTIn.M.S.F.O,#) J

Р0034(o(я» *) ;

1_F HI <H2 THEN 0O_T£ ЦЕ1

IF P0J0K3.2) THEN во ТО MEj ; FOR KJsI^eTEP 1 UNTIL N 00 RESIN NlHJtsIF 03 THEN

ХХиу(К*М)-УУхУ1К1 ELSE XKFN2) I IF WfKJ>0 THEN 60 TO ME 1 ; ENO ;

RIkMlEE j JJ:r»| ; rem:=iv j

во то mbs ;

P1oe7(N,N,S,J»,XS.PU j

РОвОО(I I,|V,X,И)    ;

1V;bIV+NN+N3 i

EO| K:«l 411? 1 UNTIL N3 00

NSlNli=IUlK|-»PHKl ;

Ройзв(Е,х»,в) ;

FOR K|sl STEP 1 UNTIL ^MN DO xstKl:*b(k)

ENO ,

P0500(II» ЙЕМ»Х E » В S) t IF l»0 TNEN REN:-REM+NN+N3 }

IF EP_Si»i THEN P1041 (*,X$,V,XS>О !

END ;

52-12 *^E3 '

f5r K:*t sTEp i un^til n oo jnIKI :=t+Mu-l for k:*i step i unJil нз £o PiIkj-.=o ; P0S00(n,REM,Pt ) ;

PTM 108.020.101—76 Стр. 11

FOR i:=lNlM»-| ВТЕР -1 UNTIL 1 DO

BEGl_N r£m:=re*-nn-n* i POSOO( 1, REM, «,»»*> i IF EPS|>1 THEN Pl04t(i>X,P2) { p|OS7(M,N,9,P1.K,V) J FOR KU1 STEP I UHT)L N3 DO pi»n»!*iMkI4P«(iu j IF P050I(3(2) THEN gO_TO ME7 *

FOR Kj*1 STjP I UNTIL N qO BENIN Л|К):=£Г «3 t«IN

k*«piIk+hi.i(u«piu] else

pllK+NE) |

pHk+n*i:*RIk1 j £F ЛIHI<0 THEN

BENIN PI IК 11 spKK+N) :жр ЦК+Ы2) : *

® I

INI N I Sr I

END

IMP I PO»00(11,R*«(P1) i If EPSt> I THEH P»04l(l,P1) )

in; i

P0717 |

IF R THEN 33;bH-I j END j

begin длину PS( i:3),psx,psyi г.Зхм) ;

PO?17(*IM,POSOO,PO»01> i for i:*i step i uhta 3«_M do PSMt 11 ;>«*»)! I I J!*0 I FOR IS*1 STEP 1 UHTlt JlF 0 THEM M“1 ELSE 33 DO

BEeiu pe5O0|l,(l-l»X(HH+H3)+1,P1) i

*:«l s

pog з:*»,а<з go

•■•IN FS*J»I*SIH(tj>1)XN+1,JX *,TZ<R>.M> I psx I (d - i    ) :«-pe I Jlxcoss l*x

UHNi) I

pSlH <J-»)|«F|+I ):«PSIJ)*S!NS12X

ими I I

END

Ш !

HRl*J:»«x8lH(l»K.T8lRJ,P8x) ; HRtR+»|:««<»l*<i,M,T*lRJ,PS»)    ;

HR|lRl:*«xS|M(H+l,2*«,TSlR),P$X) ;

МН1«а+*1:а***|ИИ+1,**М_#Т$1й),Р8а) !

psulRl: = .s*siM(s*M«-b3*M,TS(Rl_>psx) »

PSUtR+al:я.5х*1М(ахИ+Ч ,3*M,T$JrJ,Р8Ю i IF R«10RLV>.S THIR eO TO RU8 i HR[2):=JKHRlii ; hr(4) :*2xHr(4) ; hrU2):*2xhri(2j i nR1(4):*2xhru4) ; PS012)_;cBxPSU(2j . P8W{4]:s2xPS0C4J i

PTM 108.020.tot—76 Стр. 13

flus:

end end i REAt

PROCEDURE CO lx) J VALUE X ; REAL X I

bEcin xx:»x;»hm*s«nIx) j w;«нх«со8(x) ; for r:m,z do BE6IN ННСХ.И) S PR6 ;

it ^T<H 1Й1Й

>161W P«oia):*p*y(4);=o ;

60 To HOC J

CNO

ENO I

mcc:    eat    i):«psw19)+p*^uH*    I

SR(>):«t)0M(6Rl n> i

«•it):*Ptulu^ptuia) ) SRI a):«if absurI il)»»-io then sam/eRin else If 60S(ORItI)«*”«0 THEM Zx ELSE If 6R(2)>»-1° ТНЕн 1•J70796 6L*E

-i.sfOf»* ;

U AM(Sft(2 ) )>*-<OANDAB8(eR{ 1J )>«- 10 THEM SR la I! sPOOf HSR 12) i8Rl⁹)) i If SRI2)>3<14159 THEM

0Rl2);*eRiai-o.assies* ; ss:«eR(B)-zx ;

IP POS01I1,I) THEN

рЮ41(Х,*8,вЯ12),вй{ U,1T) if titUs)<»-io the_n mas:=i else ip Aes(x/ss)<i0 then mas:=ioo ; cO:«s»/MAi emo i

PROCEDURE MH(X,»)    ;    VALUE    X,y ; REAL Х,У i

вес I n xx:*x : vv:=si-Tx(r-i) !

ENO ;

PROCEourE ^cc<*.“.Z*,ZVJ    ;    VA_LuE X,y ;

real x»v»zx»zv I BEDIM xx:»x I theta:*if abb(u)<»-iо then 1.570a e_lsE arctan(x/y) ; for n:«i,a do

B*e]M HH(X.V) I PRC ;

IF T«l* THEN

BECIH PSU{8J:=PSU_l4j:-o ;

QO TO MET

end

Ш >

zx:=psui1l+psuiaI ; zy:=psu(3)+psu141 END ;

PROCEDURE CF(X,6) 5 АЯЙАУ x,6 ;

ВЕС IN IF X(2J>10Rxl2J<T-l THEN х(2):*т-1+.05_х1т/2 ;

oo n пТ1нп t iili i = :> Sod i fH'XJd

f (cJ9-(Cnc=:tn)x { *eix*t* STSaS eo N lliHn I dlls i«:c Soli I on! die OX 0*

( tdnmxcniM ilSii М1йГ <n-*>(li !•>*« -Г So и STISn i Silt i*h Soj

{ ХЭ OX 09 N1HX 0»<M*OKMl>t Zi

! Uli is+s = :s 00 n lUiin » SSI» I*n Soj

< o*:« f n'»)J

l (s+l»:i lee'i e‘CN:t'на 1*Мн:« Ih'* iSSSS

( c'i SieilST ! 9 i»i5 «Togo < a iiniiioid f x ISiSS ! m'n SaiSIST t idi lSja Sda'H SniSZ ! (M'K'd'sdj'MMifxa iifloiaoid

* oSl

П    no    oe    8*i*:    r    io5

f m'HHMi

lo'leix'ttix* <T(xcae+^,,ov-Ui'i»)09tOd

Й5Й1 (t *1)lOfOd dT I *4(В 10+8 4 t1 19a:XHi I <C*|9‘|I}Э'(8|Х*П»X>30 ! xixso**: (t )x нзнГ 2Kle)x-i)-t<84ii jx dl

91 diD 9Л—ЮГОСО‘801 Wld

*1 |,0):=(НШ-С( I ))/сс ЛЗ t xl Л :*s

вне :

РОоЗт(Я> В 1 > В2) ;

РОоЛ<Я»в»н) »

8:>о } К;=к+1 ;

FOR I:*l STEP 1 UNTIL м 00

eeelN s:=s+Ht 1 Ц2 ; Х( 1) _:=xi i ]+H[ i) i END ;

SO TO REP ;

Ex:

eno ;

PROCIOURE Р» J

белы    hxi 11: 13xmhx),hzi 11: i3«NHzl ;

P0500(t,i,hzi ) f POR I I :«| STEP 1 UNTIL NHX 00 0E61N РОТЗГ (f/ILSE, I I'M,SR) }

H*11 I I I:ȣF SR(15J80 TH|N 0 ELSE Ip IT* 1 then CRI IJ ELSE

S«RT(OR(IS)) J FOR j;s3 STEP I UNTIL 12 DO

HxiСI!+Cj-sJ*nhx):=CRIJ) :

pOR j:=25,2S oo

HX I ( n-MJ-l4lxNHXJt*eu( j) ;

ENO j

Pt07l(NHx>NHZ,Hxi ,Нп,ЯЯ, 12,3,FALSE) i

PTM 108.020.101-76 Стр. 17

Р0160СЯ1,‘3L*)

for n:»i step 1 инти ннг do poibo(ai,'зхго.огх¹,hziш » ; Р01б0(Д|,ч*) ;

FOR 3;ri STEP 1 UHTU 12 Do

BE6IH FOB |i:=l STEP I UNTIL NHZ 00

POl*0(Ab‘X-»Z0.30‘ ,HH 11 1+JXMHZ n i

poieom,'»-’> {

END ;

Р01б0(Л1,^аЭц'>    ;

poo40(hzi) ;

FOR I I 1*1 STEP 1 UNTIL HHZ 00

BEGIN FOR J!*0 STEP 1 UNTIL 12 00

P01«0(At,‘X-2Z0.30',HZtfI1+JxHHZl) f P0160(R1,‘L’)

END ;

Р0160(Л1,*3L’) i END }

PROCEDURE stpr ;

BEGIN ARRAN P1,P2.P3( UNI.4I ,4E,H10;N] ,ae, WSlOlHl ; REAL TBR J PO7l7(POl6O,POS00iSlM) !

TSR:=TSfRJ ;

IF R=1 THEN

p0160(A1,'SL20X'НИЖНЯЯПОЛОВИНА * 2L¹ ) ELSE Р01«0(Д1,'StaOX'BEPxHRRnoAOBHHfl’aL’) I

Hi:*Hxx,.cot*(0)-y»MSiNS(oi _;

IF I THt_N J₃:*M i

1£ * Iff®

H4:«1+XX¥COS8<2xMxUM)-WVxs|NSI2xM*U*l ; FOJ I*.«О *T|P \ UHTU И 00

eanliawsii);«o i

рой I:*0 «tip i untu оз oo

IlliN И*;*t+*K*COSS(ta|xU_M) -yyXSlNSl»* |Я||М1 ;

Wt 0];a I/ИЯ ;

pouoiAt,‘wax’,И i peuoim.'Jtth’i

ir,iM(Fiia*x-iI'M «тех)) ;

IE »»• 12!!!

POlOOtl,(l-t)XfNH+W3)+M*PI,p2) J

Li '>¹ 1215

P0500(b(l*l)«(MN+M»l+MI+1,pl) ;

iE ^uja I2J2

POiOO |li] x (NN+N3 )4||1+ I, PS ) i

poi x!*■ *»eP » u*rn w и oo

mm if и** ?****»*•;

i! I|‘ ШЙ    t

if ^,e,a I2!j    |

*mi«»M®>*H«»tp3(K)-iM|K))/l4x TM) 1

It ⁱⁿⁱ I2!fi

РТМ 108.620Л01—76 Сто. 19

ЧПК] :*Hi/a-HH3xUxP*lKj-,*x Р»1ЯJ>/(|2*Т5Я) ;

JF 1*М-)ДН0» THEN HE(Hl:*H4/2-H4t3*(,fxPUK)-ax Paim)/na»^,TSR) {

Р0|*0(АЬ ’K-4Z0.40' .Pf(KJ) i END }

P0160(4), ЧЧ } ws(n:*siM(o.N,Tt(R),W) i ЯВ1) );*H2»3*UxP0l »>-.3xpam >/414“ TZlRl) !

end ;

FOR K;*T STEP I UNTU N DO 1FN0TR THEN

aetiN iт:=if inik)<jo then »*»nik)xuh else Ek(33NUH*InINI-JJ) }

Н4:*1+ХХ*<*Ов»С I т l-VXxS|N*( IT I j «eiKj:=h4/8 end ;

I euoi:»Hi/8 j

6RIR+20):*SIM(O_tH,T2IR•,i gRlR+22):s6lM(0|M,TSR>EB)^DLf2 I jFNOTR THEN FOR I :*J0+> ЦЦ I UmX£L N 00 " bE61N W10>:*H4/(>+X**«C0S5(2*l*0Mi-SNx S1NS12KIXUM1)ti }

0 — угол наклона линии центров О_нО_ц к вертикали; Re^cp — среднее число Рейнольдса; q* — расход смазки, см³/с; t — средняя безразмерная температура смазки; с — теплоемкость смазки, кг • см/(кг • °С);

Т —удельный вес смазки, кгс/см³.

Схема расточки вкладыша

2. СТАТИЧЕСКИЕ ХАРАКТЕРИСТИКИ СМАЗОЧНОГО СЛОЯ

2Л. К статическим характеристикам смазочного слоя относится коэффициент нагрузки, потери мощности на трение, расход смазки и средняя температура баббита.

2.2.    Метод расчета статических характеристик основан на интегрировании уравнения Рейнольдса для ламинарного смазочного слоя. Граничными условиями являются равенство нулю давления на контуре подшипника, а при обрыве смазочного слоя — равенство нулю давления и его градиента на линии обрыва. Предполагается, что в пределах половины подшипника имеется не более одной зоны обрыва.

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

2.4.    Потери мощности на трение при числе Рейнольдса более 1000 и расход смазки через нерабочие зоны вычисляются вручную.

3. ПОРЯДОК РАСЧЕТА

3.1. Производится расчет безразмерных характеристик смазочного слоя на ЭВМ по программе, приведенной в настоящем РТМ.

PTM 108.020.101—76 Стр. 3

3.2.    Вычисляется коэффициент нагрузки по формуле

_Pf_

*    LD\>-v>

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

q* = LD²w$q и потери мощности на трение

™*=    ₄ф".........

где q* — расход смазки (входной, выходной или торцевой), см³/с; q — безразмерный расход (из результатов расчета по программе приложения 1); w*—потери мощности на трение в смазочном слое, кгс-см; w — безразмерные потери мощности.

3.4.    Вычисляется среднее значение числа Рейнольдса

где б — радиальный зазор;

D — диаметр цапфы;

v — кинематическая вязкость смазки, см²/с.

3.5.    Если число Рейнольдса больше 1000, то потери мощности на трение вычисляются не так, как указано в п. 3.3, а по формуле

_w* _ ^я,^ц<д£>³Ке^сР/_т 4Ф ’

где /_т = —--коэффициент турбулентного трения.

1\С

3.6.    Расход смазки через нерабочие зоны

q* —    j [ 1 + X cos(0 + ®)] rf®,

<p₂

он складывается с расходом смазки через торцы подшипника, вычисленным в п. 3.3, где р_а — давление на входе в подшипник, кгс/см².

3.7. Для определения средней температуры смазки необходимо решить относительно Т уравнение

Т=Т₀^

_где f — безразмерная температура;

[а(Г) —вязкость смазки, взятая из температурно-вязкостной зависимости;

1*

Стр. 4 РТМ 108.020.101—76

Т —средняя температура смазочного слоя, соответствующая вязкости ц, °С;

Т₀ —температура на входе в подшипник, °С; с —теплоемкость смазки, кг * см/(кг • °С);

7 —удельный вес смазки, кгс/см³.

4. ПРОГРАММА ДЛЯ ВЫЧИСЛЕНИЙ НА ЭВМ

4.1.    Программа расчета статических характеристик смазочного слоя составлена на языке АЛГОЛ-60 с использованием стандартных программ ИС-2 и предназначена для трансляции при помощи транслятора ТА-1М на любую ЭВМ типа М-20, имеющую не менее трех кубов памяти. Текст программы помещен в приложении 1, перечень использованных стандартных программ ИС-2 — в приложении 2.

4.2.    Статические характеристики смазочного слоя вычисляются в виде функции эксцентриситета цапфы %. Для каждого эксцентриситета, начиная с наименьшего, находят угол 0 наклона линий центров О_нО_ц, обеспечивающий заданный угол действия нагрузки Р, после чего подсчитываются коэффициент нагрузки, потери мощности на трение и расход смазки. Начальное приближение для угла 9 задается в исходных данных. Если начальные приближения угла 0 неизвестны, то значение начального эксцентриситета задается лишь для наименьшего эксцентриситета.

4.3.    При построении разностной аппроксимации ввиду ограниченного объема памяти ЭВМ в окружном направлении подшипника (ось ф — см. чертеж) вводятся крупная и мелкая сетки. Количество делений сетки в осевом направлении подшипника постоянно. Мелкая сетка, полученная при кратной разбивке крупной сетки, обеспечивает точность решения. Разностные уравнения решаются методом матричной прогонки.

4.4.    Одновременно с вычислением статических характеристик с помощью программы определяют динамические коэффициенты смазочного слоя по РТМ 24.020.21—74.

4.5.    После вычисления статических и динамических характеристик в зависимости от эксцентриситета, в тех случаях, когда число эксцентриситетов не менее четырех, могут быть получены интерполяцией значения этих характеристик при заданных произвольных значениях коэффициентов нагрузки.

5. ИНСТРУКЦИЯ ПО ЗАПОЛНЕНИЮ БЛАНКА ИСХОДНЫХ ДАННЫХ

5.1. Для вычисления статических характеристик смазочного слоя в ЭВМ вводятся исходные величины в следующем порядке: М — число крупных разбивок опорной поверхности расточки подшипника в окружном направлении;

N — число разбивок опорной поверхности расточки подшипника в осевом направлении;

PTM 108.020.101—76 Стр. 5

NHК — количество значений эксцентриситета NHZ — количество значений коэффициента нагрузки £;

UM — кратность мелкого разбиения в окружном направлении;

FI [1 : 4] — углы (в градусах), ограничивающие опорные поверхности нижней и верхней половин вкладыша;

DL — отношение длины подшипника к диаметру расточки;

LV — отношение длины верхнего опорного пояска к длине подшипника (для подшипников без выбранного верха и для цилиндрических подшипников LV=\);*

Т — отношение толщины прокладки к радиальному зазору в разъеме подшипника (удвоенная степень эллиптичности), для цилиндрических подшипников 7=0;

EPS _ приращение безразмерных координат при определении коэффициентов жесткости (рекомендуется брать Ю“⁵—10"⁷);

EPS 1 —точность получения заданного направления нагрузки на подшипник в градусах (рекомендуется задавать 0,05);

ZX — угол между полной и вертикальной нагрузками на подшипники (для вертикальной нагрузки равен нулю);

THETA — начальный угол между вертикалью и линией, соединяющей центр нижней половины вкладыша, О_н и центр цапфы О_ц, равный 90° для цилиндрических и 0° для эллиптических подшипников;

Q3 — булевская переменная, имеющая значение false;

НХ/\ 1:2 X NHX] — массив, в первых NHX элементах которого

задаются значения эксцентриситетов % в порядке их возрастания, а в последующих NHX — соответствующие этим эксцентриситетам углы наклона линии центров (при неизвестных начальных значениях равные нулю). Для эллиптического подшипника без выбранного верха НХ1[\] ^ ^ 7/2, а с выбранным верхом НХ1[ 1]<; ^7— 1 +EPS;

HZI [1 :NHZ] — коэффициенты нагрузки £ в порядке их возрастания, для которых надо получить статические и динамические характеристики.

Стр. 6 РТМ 108.020.101—76

6. ПРИМЕРЫ РАСЧЕТА

6.1. В качестве примера приведен расчет статических характеристик цилиндрического и эллиптического подшипников, параметры которых приведены в табл. 1. Заполнение исходных данных помещено в приложении 2. Результаты расчета безразмерных коэффициентов даны в табл. 2 и 3.

Таблица I Параметры подшипников (из расчета)

Тип подшипника L/D ЦЬ t Y2 Ъ ?4 град. Цилиндрический 0,8 — 0 7,5 172,5 — — Эллиптический 0,8 0,3 1,375 7,5 172,5 187,5 352,5

Таблица 2 Статические характеристики цилиндрического подшипника LjD = 0,8

X 0 С W Ti Яе Яг t 0 1,57 0 2,88 0,50 0,50 0 5,8 0,2 1,18 0,16 3,12 0,55 0,41 0,14 6,5 0,4 0,95 0,41 3,80 0,56 0,31 0,23 3,7 0,6 0,75 0,94 5,15 0,52 0,22 0,29 14,0 0,8 0,54 2,82 8,39 0,43 0,11 0,29 30,6 0,9 0,41 7,26 13,1 0,35 0,06 0,28 64,6

Таблица 3 Статические характеристики эллиптического подшипника L/D=0,8

X 0 г W 4i Яе Яг t 0,538 0,12 0 12,8 0,59 0,30 0,27 61 0,6 0,33 0,46 12,6 0,59 0,27 0,29 63 0,7 0,46 1,27 13,8 0,57 0,23 0,31 63 0,8 0,48 2,77 14,5 0,54 0,18 0,32 70 0,9 0,40 7,38 18,5 0,49 0,17 0,32 56

⁹ tos Od‘OOSOd'/cOOd* Stood‘SSeOd'HQQd Hi ^l°4 «

i ccc SieilST ! |Н:пэ'»'м'

*«• j'le+Ninse'iViNN.-nsx'x'e iiiP 37555 i ими i«d‘i* liili 37515 I 014 ll«II55id

f 151

: П4Г|«4|Г|вя«Н*:*

lo n 1111л в 3115 и Sill » 1131 •«fr'll.

*if iii

f    ill!    ИИ    йшщ    jT»:«

* e iliilit 57551 f я Iiiil ! iH i#Ii i и'» 3511757 ah'r<w inlll 8 U'ih'n'wiwi* iiiilloP

liii

! _Хи1 1551 ^! (Ninzi'ii 1555* 5151157 i itnmu*

* (Хня*«: i I    yis*

)Si'(wnxWKB:«)«eQ9‘tuts 51£§§ 57*5* ! «#*»!*# !    j    i«»i=s:n^m

f |Й«ЧЙЙ^,ХцН'*>И)»*о1‘* *• iRrt'tiliHilHM'ji'iUEyeoJ ! iirt'tx'Ex'ti'e⁰

'ledi'ein'i'Ai'ie'ti'i'itf'ii'ei'ih'tH'jU'tM 1515

!    [9!I»**'cet:H»9 *lli5 * •!** Rlilool ! W3«‘C^N

'EN'AK'A I 'ZHH'lM'CC'il 'm'l'l'f * I * i 1 ' XHH'N'w 139э7н1

MH10Md31)lVdVX ХИЮЗЬЩУЛЗ VA3f»3Vd VWWVdJOdU

эожэ&рншояэ^ l ЭИИЧЖ01ГИви i 'di3 91—10ГОСО'801 Wld

pio*n i

Р0500(1,(4XMXUM+1)*(Я-U+1,8|NS,C0SS) ;

«:*тяи« г

ГОЯ Hlmt    JTIf    1    UM?|L    И 00    1N(K):sM    {

FOR к!* 1    STtp    I    UMTlk    Нз 00    V(H):=0    ;

for к;* 1    step    1    WMT,L    HH 00    xiKi:«o    i

FOR I:ЖI    trip    1    UNT|L    M 00

posoo(ii,ii-n*(NM+M3) + bx,v) ;

REM! *NVIa I |

J0i*IF UM*< THEM * ELSE o i JJD-.rM-l ; н i: r (1+xxxco»* M)-»«**INS 11П t* i FOR Us JO ST|P 1 Wl 333 00 Beojn iv:brim i

FOR к:«1 ST(P 1 UNTIL NN DO XS|K3:sO FOR k:*1 STEP 1 UNTIL N3 DO «SlK):s0 5 FOR к: «= 1 STEP I UNTIL N 00 XSIKxH-N+K J‘a 1    :

FOR MU:a| STEP I UNTIL £F IsO THEN UM-1 ELSE UM 00

OEcin FOR К: * I STEP 1 UNHL NN DO B(Ki:>0 i на;«и» I nv:*nv+i j cosR:=cossInvj

SlNA:»»lN8tNV) j H5:=-C06RxTS(R+2)t2 i

не:=sinnxt®IR+2)12 ; H7:=(xxxh6-WUxhS)x6 ;