Товары в корзине: 0 шт Оформить заказ
Стр. 1 

22 страницы

396.00 ₽

Купить ГОСТ 25645.204-83 — бумажный документ с голограммой и синими печатями. подробнее

Распространяем нормативную документацию с 1999 года. Пробиваем чеки, платим налоги, принимаем к оплате все законные формы платежей без дополнительных процентов. Наши клиенты защищены Законом. ООО "ЦНТИ Нормоконтроль"

Наши цены ниже, чем в других местах, потому что мы работаем напрямую с поставщиками документов.

Способы доставки

  • Срочная курьерская доставка (1-3 дня)
  • Курьерская доставка (7 дней)
  • Самовывоз из московского офиса
  • Почта РФ

Устанавливает требования к заданию объекта и алгоритм вычисления функций, характеризующих экранированность точек внутри объекта-фантома с окружающей его защитой.

 Скачать PDF

Оглавление

1 Общие положения

2 Задание объекта

3 Алгоритм расчета функции экранированности

Приложение Способ реализации алгоритма (программа OPTIC)

Стр. 1
стр. 1
Стр. 2
стр. 2
Стр. 3
стр. 3
Стр. 4
стр. 4
Стр. 5
стр. 5
Стр. 6
стр. 6
Стр. 7
стр. 7
Стр. 8
стр. 8
Стр. 9
стр. 9
Стр. 10
стр. 10
Стр. 11
стр. 11
Стр. 12
стр. 12
Стр. 13
стр. 13
Стр. 14
стр. 14
Стр. 15
стр. 15
Стр. 16
стр. 16
Стр. 17
стр. 17
Стр. 18
стр. 18
Стр. 19
стр. 19
Стр. 20
стр. 20
Стр. 21
стр. 21
Стр. 22
стр. 22

ГОСУДАРСТВЕННЫЙ СТАНДАРТ СОЮЗА ССР

БЕЗОПАСНОСТЬ РАДИАЦИОННАЯ ЭКИПАЖА КОСМИЧЕСКОГО АППАРАТА В КОСМИЧЕСКОМ ПОЛЕТЕ

МЕТОДИКА РАСЧЕТА ЭКРАНИРОВАННОСТИ ТОЧЕК ВНУТРИ ФАНТОМА

ГОСТ 25645.204-83

Цена 10 коп.


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

ГОСУДАРСТВЕННЫЙ КОМИТЕТ СССР ПО СТАНДАРТАМ Москва

ИСПОЛНИТЕЛИ

П. А. Барсов; А .И. Григорьев, д-р мед. наук; Е. Е. Ковалев, д-р техн. наук; Л. М. Коварскйй, канд. техн. наук; Е. И. Кудряшов, канд. техн. наук; Е. Н. Лесновский, канд. техн. наук; В. А. Панин; Н. М. Пинчук; И. Я. Ремизов, канд. техн. наук; В. А. Сакович, канд. техн. наук; В. М. Сахаров, канд. техн. наук; В. Б. Хвостов, канд. физ.-мат. наук

УТВЕРЖДЕН И ВВЕДЕН В ДЕЙСТВИЕ Постановлением Государственного комитета СССР по стандартам от 20 декабря 1983 г, № 6360
ГОСТ 25645.204-83 Стр. 9

606 CONTINUE

IF(I0ENTI511.511,512 51? READ 5fXC .Y%Z0fUO,V9.wn

CALL GEOMIN(XOfYO,ZO,UO,vn,WO,Ll) PPINT 70,X0?YP#70,UOfVPtWO

70    P0RM&T ( ?Xf3(t:ln*3) ,5Х*3(С10«3) ) PRINT 71,L1

71    FORMAT(9Х»14)

DO 81 IZ-1 ,1.1

P.PINf 72s IP(IZ ) rSPf 17+1 )?KP{ IZ )

72    FORMAT«5X,I4,5X,E10.3,5X, 14)

81 CONTINUE

GO TO 50 ?

511 READ lb ISTOR 11 FOR MAT(16)

13 CONTINUE

READ l?fX0.Y0.7O 1? FORMAT ( 3E10 «л )

c ----------------------------------

DO 306 IN=1,190 DO PT11N i =л «

DOPT1C TN)-0.

00 PT21 IN?-0«

DOPT3CIN)-0.

D0PT41IN)=0.

DO P T 5(IN 1=0.

306 continue

DO 400 IS=1#I5TDR .

CALI POMEGA(UO.VO.WO 1 CALL GEOMIN(XO*YQ ? 70sUO tVO ?40 * L 1) 01=0*

DL 1=0•

DL2=0*

DO 308 KL=2?L1 NL=KL“1

IF<KP(NU *EQ*CM GO TQ 30?

IF(KP(NL)•GT*1) GO TO 351 OL 1=DL1 +SP < KL l-SPINU GO TO 308

351 DL?=DL?MSP(KL)~5P(NL>r*RO(KP(N! )) 308 CONTINUE 307 CONTINUE DL-OLI+0L 2

IN0=1NDfX(OL)

INDWNDEXfDLl)

IND2-INDEXCDL21

DOPTC INfH = DOPTUNDl*I •

DOpTK I ND1 ? = D0PTH If 01? + 1*

DOPT?U4D?)=DOPT-4 IND?) +1.

400 CONTINUE

PrTnT 220

PRINT 430 »ISTOR »JO s- YC a 70 ^

430 FORMATl/IOXt МИСЯ0 KCToPWf? =• » I6.5X #*КОГГ*ОИНАТЫ T Л'*ХМ» * ?Xf 4 XO= S

Стр. 10 ГОСТ 25645.204-83

DIDO

0101

sue/

01 '3 Ci£4-

0105

*Е1С.З,Ч '5 -5 ,2Х.в YQ-r* ,СЮ.З.ЧГН) ,%2 ХЛ2^* .EI0.3*U *г/* PRINT 42#

РР I N7 4 5*

450 FORMATS'4,Л Ч 3    »9К,а    Р    *    ,    7Х.«    И    f    tlQXt    3    PI    '-Т.,5    V

М,1ПХ.в *?    »9 л * )

PRINT 2'2

222 FORMAT! 12% ? * ' /гм**?) 9 * 6Х. •    •    9<5Х*    Ч    С‘«*    *?/Г    р    f    7у    *    •    «    ?    СХг    5    \*    *■'

*2/Г>Ч8ХЛ    45ХЛН 4**?/Г)Ч

слое

01П7

ОШВ 010 > sup

О i 11 01X2 V i 3 0И4 0115 011 6

0117

0118

0119

0120 0121 012?

0123

0124

0125

0126 П127 0128

0129

0130

0131

0132

0133

0134

0135

0136

0137

0138

PRINT 420

ПР 500 IL=I * 50

0ПРТ(1L1=DOPT<U 1/ISTOR

DOPTHm =ГЧЛТ1 { ILI/ISTCR

ПОРТ? UL ШиГ 5Т2 J !L1/Tf " О?

IFnL.Lr.il 5 'Л ТП 501 miL.LT.lM OP ТП 502 iff iL.ir.3?) е;п тп 50з IL 1 = < IL“37 ) 40+100 И ?=Ш*10

DOPT3I I L ) =ООРТ ( IL ) *0* 1 ПОРТ4( Ш=О0РТ1( IL 1*0.1 ППРТ5 ( lD = Dl)PT2(IL ИЛ. 1 GO ТП 504 501 IL1=1L-1 IL 2= 1st 1 +1

ОПРТЩ ID =ОПРТ( IL)

D0PT41 ID = ООРТ I ( IL 1 ООРТ 5(IL)=DPPT2<П 1 GO ТП 504 50? IL!=(IL-m*2+lO IL2=IL1+2

DOPT3( IU=OOPT(!L 1*0.5 ППРТ4( IL)=00PT! ( ID*0.D0PT5I IL ШПГФТ ' ( IL)*0.5 GO TO 504

503    IL l = CIL-1645+20 IL 2= TL1 + 5

DOPT31 IL)=DOPT(IL)*0•2 DQPT4IIL 1 = D0PT1 ( IU*0.2 DOPT 5(IL)=DOPT2(IL1*0.2

504    CONTINUE

PRINT 4tO.IL * IL 1 .И_?.ПОРТ( IU , DOPT 3 ( ID t DOPTl (IL) * DOPT4IL b ♦ DOPT2T m*00PT5( It)

9

0139

0140

0141

014*?

0143

0001

D0O2 0OSJ 0004 С 005 С 00 6

410 FORMAT! ?X,J3*5X,IM 9-9,14,2X.2!2X»E10.3)*3X,2<2X*E10.*Ь3X.?I?X *E10.3))

500 CONTINUE PRINT 220 220 FORMATC/2X, 110ПН*)/)

GO TO 13 END

SUBROUTINE GEOMINIXO. YOtZ0#Ua.4VO.WOftU

COMMON/AG/UP »VV»WRfAi50.10)*R£H 30)

COMMON/AG1/N.NZON.ITJ 50 ),NC450bMI (30 >. IPZ П8«6) rIDttG.6) «KPN< 30) COMMON / A G2/ fOI( 100) . (P(100 IfSP11081»*Pt 100 J SP( 1)=Э.О IP(11=0

COOT

11*1

0008

Ш1

S5H-i5i?^0;Xc''r3,U0»v5fW!>*Ln

IF f L 1-2 ? *3 j 4*>«* ч

00 £.©

44

L=U“1

ООН

00 46 J=2,L

0012

Al^SPC J1

8013

K1=0

0014

JV=J+l

0015

DO 47 I-JVfL1

0016

IFCAULE.SP(D! GO TO 47

О0| 7

A1=SP(11

001 в

K1 = I

оою

11 = 1

0020

47

CONTINUE

0021

IFfК1)43,46,<8

0022

48

SI=$Pf J)

0023

12=IP(J )

Э-^24

= u>

С025

f P i111-12

ОС ?6

SPC J)=A1

0027

SPiI1)=SI

0028

46

CONTINUE

0029

45

CONTINUE ПП 19 J=2 fLI

0030

003]

1=1 PiJ)

0032

DO 19 J1=1,N7PN

DD33

j^^KPNC J1)

0034

J2=MI(Jl)

0035

DP J3=lfJ?

0036

J4« IPZf J1fJ3I

0937

JFf IDI Я /J3>*7PH Jb))2Q,l9,19

0038

CONTINUE

0039

KPf J-l)=J5

0040

IDI ( П=-1«ЮН I )

0041

GO TO 1

0042

19

CONTINUE

0043

KP(J-15-0

DD44

inim=-i*iDim

0045

1

CONTINUE

0046

18

CGNTINUE

0047

KP(L1)=0

004$

23

CONTINUE

0049

RETURN

0050

END

0001

FUNCTION RANNQCNMB!

С

*** РЧБОР СПУЧДЙНОгО чиспл НД 1

зш

OAT* 1Х/1/

CALL RANDU?IX,IV,YH)

0004

IX= IY

0005

PANNP=YFL

0006

RETURN

0007

END

0001

SU3RnU4NE Rr-:3U(IX,IY,YFU

0002

IY=I v^6ff'.7 39

0003

IFIIY>5,5.6

QO04

5

IY=IY^214748364T+1

0005

6

YFL=IY

0006

YFL=YFL*с 46 56613f~Q

0007

RETURN

0008

END

0001

SUBROUTINE ROMEGAfU9,vr,W0)

0002

DATA NMR/1/

оооч

CALL CSri(rFbSFI)

0004

CTET=1 .-2.*RANNP(NMR)

0005

STET-SQRT< l.-CTET**2)

0006

uo=stft*cfj

0007

VO* ST FT^ SFI

0008

wo=ctft

0009

RETURN

0010

ENO

0001

SUBROUT INE CROSS UJt * VR*WR ,L 1 , SP* IP* I

0002

DIMENSION $ P(l^ ), IP(IDO)

0003

IFlWR)l,??l

0004

2

IF(VF*NE*C *0) ГП 1 J ^

0005

RETURN

0006

1

0=Vft*VK-UP*WR

0007

IF(0) «,9,4

0008

я

S=-VR/W»

00Э9

IF1S.LE*L.C) GO Tn 5

0010

I 1~L 1 H

ООН

SPUD =5

0012

IPCtn=Il

0013

5

RETURN

0014

3

S=-UR/f 2'>VR)

0015

IFCS.LE.O.O) GO TO 6

0016

Ll=Lin

0017

SPC11HS

0018

IP«U) = I1

0019

6

RET URN

002С

4

D-SQRT{D)

0021

5*(-VR-D>/WR

0022

IF C SeLEe^.O) GO TO 7

0023

L1U1 + 1

0024

SP(LI)*S

0025

Ip U I) -11

0026

7

SU-VRO)/WR

0027

IF(S•LE«0*Q) GO TO 8

0028

LI=L1*1

0029

$D(L1) *S

0030

iPfLl)=il

0031

Я

fETURN

0032

FND

0001

SUB ROUT I N*" STS{ XO* YO* 70 *iIC,VOfWO*Ln

00С2

Cn,«MnN/AG/l»R»VV*WR*Al D*10) tROf 30 5

000 3

COMMON/ДО1'NfNZQN*ГП50)*NCF(50)•MIC

0004

COMMON/Г G?/IDH109)*IP( 130 ) , SP ( 100 ) *

0005

ПП 17 11=1,4

0006

j = it cm

0007

GO TO (l*?*U4,5,6,7*Qf9tl0*ll*l??15

0008

1

CALL TIP1(11г XD ?UO* IDD)

0009

GO ТП 15


KPN(30)


ото

?

САМ Т IP ! { u fvn ty t in0)

ООП

0П TO I*

эп?

3

САП Г Г°Н I b ? t* *1 nni

Об 1 з

GO Tn 16

0014

4

CALL TI02( I b 7r fY^fV2 ,W3, ПО)

2*45

GP T0 16

001ft

6

САП HP2( I 1 ,7^*V0,umw"f X OD)

ООН

GO T ° 16

ОС 18

6

CALL I 1 » Y ,* ilh »VMnn|

0010

GP Tn 16

0020

7

ГАЦ HP3f I1,X' , Y' » 7ftfMntV'»,WO»IOr))

??? I

GO 15

0022

3

CALI TI»3< H fY'’f¥Of 70iVPftJ0fvr,.f pf))

0023

gp m is

9

r ALL ^IP3( П f ZOf Yfir Y0»wr fir »vo» iooi

0025

GP TP 15

0026

10

CALL T 1 P 4 ( 11 , Y ■* t 7 9 f Wл , I no )

РС?7

GP Tp 16

0028

1 1

CALL ИРМ 4 ,Yr- f 7%<jr,w'P,IGD)

0029

с*з?

1?

Gp ТП 16

CALL TTP4C I 1 ,XC fYOtUOfV^tTOO)

0031

Г,п T П 16

0032

13

CAl L т 1P6( П , vp # Y% 7 b 1*0. VO f WO t 1001

0033

GP to 16

ост

14

ГАН T IP6( I I. Xf*, YO» 70 »UP, VP, WO f TDD 1

0035

15

CALL rpnSMURtVVrWP,Ll, SP.IPf II)

0^36

mm 1jno

ОбЗТ

17

CONTINUE

0038

PPTUPN

0039

END

0Э31

г

FUNCTION INDEX f GL)

*** ВЫЧИСЛЕНА НОМЕРА Y4ACtkA ДП6 ЗАДАННОГО РАЗЬИЕНКЛ ТОЛЩИН ***

0002

J F 1 DL - 1 0 . ) 40 7,40 2,40 3

0003

402

IN0*OL

Э?*94

1N0=IND+1

0С05

GO TП 4H

С00 6

4Г 3

IF<0L-2 4> 404,404,405

0007

404

XIND*f0t-10.)/2•

Л8

IND = *INDH 1

0009

GO T П 410

0010

405

IF(OL-130.I 406,4(6,407

ООП

406

XIМП = (PL- 20 -)/5.

0012

INH=X1N0+16

0013

GO TO 410

0214

!F<0L-29Г И 408.408,409

0015

408

XINO*(OL-lOf.1/10.

0016

lNP*X?N0+3?

0017

GO TC 410

0018

409

IN0=51

0019

410

CONTINUE

OQ20

INDEX*IND

0021

ftFTllRN

00 22

CNO

от 1

0Г,пз

Щ04

проб 0006 0007 ОЭС 8 000*9

0001

отг2

0001

о*чк
2    -*5

0006

000    I ОЭС 8

0009

0010 ООН

0001 0002

0003

0004

0005

0006

0007

0008

0009

0010 ООН 0012

0013

0014

0015

0016

0001

ООО?

0003

0004

0005

0006

0007

Q008

1)Г"9

0010

ООН

0012


VJ*nn«JTI-*|P TIPI ( II, Х,и,ПН>

COMMON/AG/UP fVV , WR,MS tlO) f^n(V } UR=-A( Т|,ПИ I inr-i

IF!UR.ОТ.r,01 ID0= I VV=U/?

WR-0.0 RETURN END


SUBROUTINE TIP?( Il#7fYfVtW»|09| С'*ЧМПМ/АО/иа ,VV,WR ,AC50,10>,ROC 30)

S=A(T1.3)-A(!1,1)

UP = SMM!!,?)“7}-Sl*fAUl?l)-r)

IDP^- 1

I F ( UR. GT #no 0 ) IPD=1

vv=si/?*v-$/?*w

WR-0.n

RETURN

END

SUBROUTINE НРЗС 1 1 , vvVf 7tu,V, WtlOH) COMMON/AG/UR,VV,WP,At 50,101,RO(30) S-MIl*41*ftni,5)-Afn#M*AUif4)

$ =— c

51- M 11 f5)-A(11»3)

S?=A(Ii,6)-A(11,4)

UR= ( S1*X+S )**?- S?*^*UA(Iirl)-YHr'f ( \ ( H UR = - up

V/ = X*'J*Sl**?+S]*S*U + S?*S?*f V*( A< !1»1)-Y)+;

VV=^VV

IOOs-1

IF(IJR.G^.O.O) ТП9=+1 WR=S1*SI*U*U-S2*S2M V*V+W*W)

WR=-WR

RETURN

ENO

SUBROUT I NF Trp4mfY.T,V,WFI0n)

COMMON/AG/UR, VVrWR, A (s;, 1«;J fpp( 3A)

S*Af I 1 , n*AU1 ,4)

S1*MI1,3MAC 11,2)

52- AU1 »?)*M!1,4)

UR = S* S+$l*$l~S2*    S‘Mn»4)?Y*?*SJTA(H

HAM ) , 4)**?*Y»Y*AC II, ?1 **?*Z47 IDD=-1

IRUR.07.'',0> 100=1 VV=-T*A( I 1,4) *V- $' T A ( I! ,?

WR = M T 1 r 4'** ;» v*\M A (

RETURN

t=ND


, M-7 '+4 м * С A < I 1 , ? ) - / i )


,?>*/

'▼Mil,?'4*"



Л001

по?

оооз

0054

ОООЪ

0106


C3D7

оэоо

DJG9


ООП

0012


г П ' П.х

7U° * Vt WP f AI 5

' Д< < 1 ,б )

у Ъ )

5°- М г* ,41*М П|?! ш*м I и и** с

j^s»* s-:aac n »t )***s*x- ?*ac и , -u*si

?^*S'

mo---!

H4<P* ,T.C*€) ino=i

VV=-A{T1 - 1>*<*S*U-A(П t 3)*Sl*Sl*V-A(11f 5)*r?*S?*W* I S^S^iJ^X + Sl^Sl^V^Y + SZ^S^^W^Z VfR«S*$*U*ll + *1 *S I *V*V+!::?*c?*W*W

PC7URN


I/ »

СГ^ОЬ/ * *


? s ! > f V Г)


О)


< П *3 )**?*S1*SI+A(i1l f 5 )**?* ** 2*RZ-A< ! I * с‘^“ и . л    S1*Y--?*a    (I: г 5) $ "*S?*Z*


ООО?    SUB^in^r TIP6U 1 tXfYtZ*UtVtWf IODJ

0002    COMWW/AO/MPfVVtWR ? A (5^ ,101 »PO( 301

0003    UR=AU 1 V I > +MI b? )*X* A( II ,5HYtMI! f4>*Z*A (!Ь« > *X* У*Д (I Ь91 *¥♦

IZiAIII,)0)*X*/ AIM,5)*X*X+A(I!*6)*Y*Y+A(11♦7)*?*l

ООО7    IDD*-I

0005    \Р(И».9Т*0) I DO* i

0006    VV"4,5MA( I1.2>*U<-A4It3)*V + A4it4»*WUAni»5)'*U*AfMII »6)&V*Y+

IA ( 11 ,?)*W*Z+6. 5*1 A I II ,8)*<V*4H)*Y)*A! 11,9>MW*Y*V*Z>*A(IЬ10>*

?< ^x + u*ni

0007    WP -A ( П » 5 )*U*U*A(Г ?    ,7)WA(!1    ,8)    <'j*    V: A (11* ©

1^(11 ,10)*M*W

Огг:э    ю j°n

ОС *1 ^    с ч и

0331    SUB40UTI ME

C    *** PPW^fcJHJ i ПУЧАйги< СО!>1 (Cl И SINP! C SI ***

00v 2    DATA NMB/1/

CPC 9    FI*RANNf4NH3)    fc*?»5

0034    C~COS(FT i

0D04    S-S IN* F Г5

0006    RETURN

0D07    E^O


4. Пример расчета

Объектом расчета является цилиндрический фантом, задаваемый в соответствии с ГОСТ, экранированный с боков цилиндрическим слоем алюминия высотой (Я) 127 см с внутренним диаметром (D) MX) см и толщиной стенки (S) 10 см

Точка с0 расположена на оси Z на высоте (h) 54 см Геометрия объекта приведена на черт 2

Пакет входных данных для задачи и результаты расчетов по программе OPTIC представлен ниже Время счета данного примера —2 мин на ЭВМ ЕС 1040



Геометрия объект*


1, 2, .... 7—номера, присвоенные поверхностям (V==7); I, II, V-* номера, присвоенные зонам (NZON~5); вещество фантома — в зонах I н II; вакуум — в зонах III и IV; алюминий в зоне V/(NMAT=*s3)

Черт. 2

ГбСТ 25645.204—83 Стр. M



Стр. 18 ГОСТ 25M5.Z04—ЙЛ


1

!

_

*

»

гг

*

}

к

сс

5

j

•"^ гЧ

<1

i

1

5

СО о

1

1

О

I

»

•х

ЬАч ш

!

U

CJ Г*

1

!

W

«-4-

Г

j

Г~

• • # • •

1

1

1

5

QL

сосоо

I

1

!

4

Г*

1

1

2:

<1

*

1

а.

•“* чН (\i r\i fTi

I

1

1

»

1

м

1

i

fe

1

1

У

I

Tfr

1

с

оссов

к

I

к

1

*

!

Й-

1

*

1

j

1

чГ

1

*

1

rt

<1

f

IT

{

!

**f

1

1

LT,

ООСОС

Jk

1

*

1

■н

1

l\

if

1

tr

1

X

«

1

■%

1

«3

1

•X

1

•£.

*

1

f

*r?

.*

I

€5 О""' гчН

г*

«Г

1

***

'*

1

к

{

*

1

Сь

*

tr

?

к

1

►м

н

<1

i

_1

#

{

X

к

!

c

к

1

<

*

\

or

к

\

гЛ

чН put< НИИ

У

*

\

*

1

1 ! 1

*—«

ii

t

u

к

1

W-

■н

1

X

к

1

■л

I

*:

к

1

S

к

CVf^fMfVj }

к

1

н

яг

C COC {

tt

а

1

rvj

о

*

1

X

к

1

Mill

С

к

U UJlbU )

x

к

\

т

м

<5

^CCw 1

Q

Я

1

X

ч

i

f—

<н>

!

О

К

• • • • i

Lb

«

1

1

к

CC C C I

5

Г’

1

1

1

1

I

г

I

C

к

1

с

-ft

1

u

К

1

г

и

1

w

к

1

X

1

й

I

^

я

1

•к

t

•Й

1

соооо

>

Г*'

I

*

1

к

<i

1

я

г

ec c r i

1

• • • • i

г

-?

^ C^C l

j

■*

1

*

1

“ОСССС

л

1

к

1

♦-♦ечкмгч \

S

|

fc ООО t

1

К

i

I

о.

(V

ШШШШ I

1

►ч J-

О Ci* чС sDf^

с

COCO 1

1

КЛИтяБ 1

1

* * • • 1

1

C COO J

1

1

i

1

1

го

it ^

CM(\|Cu J

{

1

C CC 1

ft

1

f~i

uiUjai i

1

1

«а

f

\

см

О-Н^СЛ'^С 1

1

• ••*••* 1

1

COCfCCO 1 I 1

1

t

1

1

1

гН

U

i

I

1

*

С

I

1

1

1

1

1

1

I

|

»м«

чг

h-

P^f^rr-fMNfstfV 1

1

!

»~»

«М«им>И^4 j

1

i

1

1

1

1

1

ич

г

t

и-

1

*

1

)

***!№***$£*****#*******$*** $~*******ф**эМ *******************************^&***************:к***%*****?*«*£


УДК 62—784.7:539.16:629.78:006.354

ГОСУДАРСТВЕННЫЙ

Группа Ф40

ССР

СТАНДАРТ

СОЮЗА

Безопасность радиационная экипажа космического аппарата в космическом полете
МЕТОДИКА РАСЧЕТА ЭКРАНИРОВАННОСТИ ТОЧЕК ВНУТРИ ФАНТОМА

Soacecrew radiation safety during spaceflight. Computation methods of points shielding inside fantom


гост

25645.204-83


ОКП 696800


Постановлением Государственного комитета СССР по стандартам от 20 декабря 1983 г. № 6360 срок введения установлен

с 01.01.85

Настоящий стандарт устанавливает требования к заданию объекта и алгоритм вычисления функций, характеризующих экранированность точек внутри объекта-фантома с окружающей его защитой.

Под защитой в стандарте понимают конструкцию космического аппарата (КА), его оборудование и специальное снаряжение, защищающее (экранирующее) космонавта от ионизирующего излучения.

Стандарт предназначен для подготовки исходных данных, необходимых при расчетах на предприятиях и организациях, занимающихся научно-исследовательскими и опытно-конструкторскими работами, связанными с обеспечением радиационной безопасности экипажа космического аппарата в космическом полете.

1. ОБЩИЕ ПОЛОЖЕНИЯ

1.1.    Экранированность точки г0) расположенной внутри фантома, характеризует функция экра

нированное™ w (|, го) такая, что w (g, Го) d£ представляет вероятность для лучей, изотропно испущенных из точки го, встретить на своем пути суммарное количество вещества фантома и защиты \ в интервале от g до g +    выраженное в массовых единицах длины.

£=£i+i2,

где — количество вещества фантома;

|2— количество вещества защиты.

1.2.    Под массовой единицей длины в веществе понимают произведение линейной единицы длины на плотность вещества.

—►

1.3.    Самоэкранированность точки г0, расположенной внутри фантома, характеризует функция еамоэкранированности w\ (£1, г0), тождественно равная w (g, г0) $8=<ь

1.4. Экранированность защитой точки го, расположенной внутри фантома, характеризует функция экранированности защитой w2 (g, /о), тождественно равная w (£, г01==о .

2. ЗАДАНИЕ ОБЪЕКТА

2.1. Объект, в виде выпуклого тела, задают совокупностью зон с постоянными физическими свойствами вещества в пределах зоны. Каждой зоне присваивают номер К— 1, 2, ..., Ктах> где Ктах — максимальное количество зон, необходимое для задания объекта.

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

Издание официально*    Перепечатка    воспрещена

© Издательство стандартов, 1984


ЧИСЛО ИСТОРИЙ * 10000    КООРОПИАТЫ 'точки    Х0= 0.0 <СМ), YD= 0.0 ССМ), Z0= 0.210?    0?(СМ)


N


I

I

I

I

I


I

I

*


I

I

I

I

I

I

I

I

i

I

I

i

i

i

I

i

i

i

i

i

\

i

i

i

i

i

*

i

i

i

i

i

i

i

i

i

i

i

i

I

i

i

i

i

i

i

i

i

i

i

i

!

I

I

V


i

I

i

I

I

l


i

<

i


•я

•*


ft


(Vi

#

II

U


Cl


U

>4

N


XT

О


a



o.


ГЧ


*

*

~>T

О'ч


rvCv C\l <\ <\j <\ f\i <V <\ СЧ (\ rv *— >~* *-, f\

еалосссссооеос^гс" c<-*

YTi i i i i i i i i i i i i I I i i

LLU*U-L-UuLL LULUUjIUU. ILUilULL

вс. г* С о С С О О С С С СИЛ IT tf\njO O' СО ОС ^ S ^ N ^ m о * К О 4- <Ni С Г-14» О

С С'СССС twr


Cc<P0CCCC0CC0w60u0DCc.CCCCt сое


< с с С с «л


• * •

«С Cl


LtCctU-tCrCCl t'f. ct ос. сс ^coCcr ©а


Cf\ ГЧ:ГУ^Г\(\,ЛчГЧ|ГЧГч|Р\.Гч. ГЧ/ОО*-»*-*

лг^О^СОСООСООССс СОСО I Т ! I I I I I I I I I I I I II

oCq/JOCCOCCCCCOOCO^CM^ f-nS-f ^rtl'CfMfsar^oet^C* очлгсгс    ^

i-4^fsjfsi<sjprif^<Mf4tOC-4-f4jS Счсес С СС С С * СС COv

ccVooococcССССОССССС сссссссс. tc^x


с с «ее* cot ос СЧ Ос С с ^С" -С


CCOCCCCCwCCCС СПС I I I II I I М II II I м

LLUjUJILUJUJU U-UjUJU UJUjU-U LLi fCCCN l^^Cir^t^fSurv OGMCCtC r* PKlCCCrVCfMU'CCTvOr-b fr ^

Co<r OOCOCC О О h-LOUMTfT'^s04j-rvfv<t HRHh Cccc.e £ - ’ c    wCf    u    c

ОС 1    «Л.    Ct    CLCCt    I    О    ■    Xv    .    *^'4    ^    ^


CCCC СС^-н->^-.М»тол occooooccocoxcc .    I    I    II I I I I II

iLlU'lbU U- la LL LI Ц|ИДЫ1 HULL IT *.1^ «tf'^vChfVf-HO'fvCCC C*4 «"нООК^ГСГч 'tf^cnsOCOff

-r с c CCCC    «tf'Of^ccc.    Cc

.•••••••••••••••••••«••••••••a

CC^CCCC 050C Cf-orcc Cc CtCCO- CwO c c


c, С СССч. ,CC' C C OCC c c ••«•••••••••••••

»    /    C    С    ОСГ    C    c. " c


t

ft

a

к

ft

к

A

*

-A

ft

ft

V

*-

h

*

A

\

ft

#

Ц

•ar

I-

ft

ft-

*

•Ip

*4

•я

n

*

1*

'ft

ft

ft

ft

Л

4

a-

*K

4

я

•ft

•ж

ft

к

■*


ГС


r-ч r-4 •— «—t f4! (\l (\|

C CoerC COIcccco I I I I II I II I I M I

ULLLia; LLtLlI u.UUUllUUU.

eoir or^ oc^ca-Hf^i сеема

h* ^sT'COCCM< —чггсог^ОС-Гч с CCCCOC ОС CCVtnxCTCCs0ei4^^^'J‘f-<—c c


CO< C - cc Ct k-c 04, СС С С СС Сс CCC СС ОС


(М*ин*-»*Ч»"СС<иОИ#ч(ЧЛ

ccc cccvc* eoccc I II I I I    I I I I

шша lilU,UJLU|jjlLU UJUjUJILI

o<ro^t(M<rccctra octree siC'Ncrrc 'tet't с »4' CCCCCCCir^H^^fTrtMM^inMeCrir. 1C

0CJ30 30C0CCC0CCC OCCCCOCWOCCOC4 cc


^ C * Г ОС    ч*с CC"

< Cr.C.CiC-C    wc.t c


^    ,    C~OOrcC    C^CC'XC

‘С    СО    ЮОООСЮССОО


;-Г<ГкГ vor^ cc0sОГО0CC.LTOieiClAOinOlTiC LT'OLTr 1ЛОСС t ^^^rH^fMrVfTirr»^ st ITilO vOC-N r^OCttC O' Or-nrvjfT


/?Cf сое cot ooececc

Y\ (vrv Гч Г'-.СЧ f4i<Vj(Ni


I I I I I I I I I I I I , I I I I I I I I I I i I I I I

c^oiroiniiroinoor^'?:^    CrjCOi.    CC


S I I I I I I I I I I I I I I I I I I I I

A^r^sf ir^<£^c3<w>occairicir©u-\c- ___ ^    „.     Л    ^    +tr.    _

j -|Д || » . r>-<    r~4 r~j f\ P4J f\J f\J (\J <\    f\,CV(V


4f 1Л CO


t^^iVfVfsjrgrvfMf^ rMisifS^frr^ тглгг^г0


•Dh ССесмЛ fT> IT sOh-ODO'O fT fT |П гг ^ -d* 'rf- s?    ^ ^ ^    1Л


»

*

•4


A

-*

*

•ft-

*

a

ft

#

#

*

*

•K

*

ft

■Jr

•ft

•ft

*

•ft

«

*

ft

ft

ft

«

*

#

•ft

#


2.2. Каждая зона объекта должна быть задана вектором поверхностей /*■, вектором неопределен-

ности а к! (г к ), индексом, характеризующим принадлежность вещества к фантому или защите, и плотностью вещества в зоне рк.

Вид поверхности

Уравнение поверхности

Тип

поверхности

Максимальное число коэффициентов

Вводимые коэффи циенты

Плоскость, перпендикулярная оси:

X

х-с

1

1

С

У

У=С

2

1

с

Z

z=c

3

1

с

Плоскость, параллельная оой:

X

r-Ki y,-v\ Z—Zx~ z3—zx

4

4

Ух» ^1» ^«1 ^2

У

х—хг X3—Xx Z—Zx Z3—Zx

5

4

Xj, Zlf Х3, z3

Z

X—Xx Xa—Xx

6

4

Хх, ^1» ^*1 Ki

У-УГ уз уi

Конус, параллельный оси: X

V(y-yi)2+(Z-z^-R,__

X—Xx

7

6

^1. *i. Дх,

■^1* ^2» ^2

-

У

V (X-Xx)2+(z-Zx)2- R1 _

8

6

•^1» RXt Ylt *lf Уа

R'i Ri

Z

у (Х-Х1)2+(У-К1)2- Rr _

9

6

Хх, Ух, /?х»

Zx, /?3, z3

Z—Zx

Rl

Za— Zx

Цилиндр, параллельный оси:

X

(K-Kx)2 (Z-Zx)> .

a 2 + 62 ~1

10

4

У1( я» Zx, ^

У

(X-Xi)2. (Z-Zx)2 .

a2 + 62

11

4

Xi » , Zx, b

Z

(X-Xx)\ (У-Ух)2 ,

a2 + _1

12

4

Xlt a, Ylt b

Эллипсоид

(X-Xx)2 (У-Ух)2. (Z-Zx)2 . a% 4- + C21

13

6

Xx t » Уг j b» »

Поверхность второго порядка об-

Й1-]-СЗ'аХ-|-А[зУ-4-Й4 Z-|-Д5Х 3-|-(2ьУ 2-}-

14

10

«1, fla i • • ■ i

щего вида

H-a7Z2+a8xy+^9yZ+a10XZ =0

2.2.1.    Поверхности задают в виде уравнений 1 и 2-го порядков в декартовой системе координат

r={X, Y, Z} в общем и (г')=0 или каноническом виде в соответствии с таблицей. Каждой поверхности присваивают номер t= 1, 2,/тах, где /тах — максимальное количество поверхностей, необходимое для задания объекта.

2.2.2.    Совокупность номеров поверхностей, ограничивающих К-ю зону {1}к> из множества но-

—►

меров поверхностей {i} (i—1, 2, /гаах) образует вектор поверхностей

2.2.3.    Каждая поверхность ut (г') = 0 разделяет два объема: внутренний — щ(г ')<0 и внешний — щ (г')>0. Принадлежность точки г* к внутреннему или внешнему объему характеризу-

ГОСТ 25645.204—83 Стр. 3

ют признаком, именуемым индексом неопределенности б/ (г*), значение которого определяется выражением

(1)

2.2.4.    Все точки зоны должны иметь одинаковые индексы неопределенности относительно поверхностей, ограничивающих ее.

2.2.5.    Совокупность индексов неопределенности произвольной точки г* для вектора /* образует

—> —>    —*    *•

вектор неопределенности ак (г*). Вектор неопределенности для точек /С-й зоны записывают как а я(г)/

3. АЛГОРИТМ РАСЧЕТА ФУНКЦИИ ЭКРАНИРОВАННОСТИ

3.1. Функцию экранированное™ w (£, г0) вычисляют в виде функции wn~> (г ) кусочнопостоянной на отрезке (|; , |/+i )

wW (

где 1= 1, Lmax — ном:ер отрезка;

о»^)

], если ?|<6(г0> о, если g(rT, 2)<E/;g(r^    ;

I (r0, Й) — количество вещества, встреченного на пути луча из точки г0 в направлении й.

3.2. Для определения функции w (£, г0) необходимо задать расчетную сетку {£,} в диапазоне 0<l<lmaxt причем ширину интервала Д| следует выбирать исходя из требований к погрешности функционала, вычисляемого с использованием w (£, г0).

3.3. Для вычисления величины g (го, й) необходимо определить расстояние, пройденное лучом в зонах объекта, что требует выполнения ряда операций, изложенных в пп. 3.3.1—3.3.7.

—> —^    —►    —►

3.3.1. Вычисляют расстояния 5 (г0, й, i) от точки г0 до пересечения луча в направлении Й со всеми поверхностями, решив для этого относительно S совместно систему уравнений, описывающих

поверхность и прямую в направлении Й, проходящую через точку г0

(3)

U0=_o

1 Р=7>+SQ ,

где 0<S<oo — расстояние от точки г0 по лучу Й до пересечения с i-й поверхностью.

Система уравнений (3) для каждой поверхности может иметь одно, два или ни одного решения, что соответственно означает однократное, двукратное или отсутствие пересечения i-й поверхности лучом.

Полученным решениям присваивают номер п (п= 1, ..., N, где N — максимальное количество

пересечений лучом поверхностей объекта).

*    —►    —►

3.3.2.    Располагают полученный массив значений Sn0, й) (я = 1, ..., N) в порядке возрастания,

—> —>

формируя при этом последовательность соответствующих номеров поверхностей i п0, Й).

—>    —у

3.3.3.    Вычисляют длины отрезков t п (го, й) между последовательными пересечениями

^п(го» 2) = Sn(rо, Й) Sn-\(rо, 2),    (4)

^ __ —>

положив S0 (ль П) — 0 (пересечение лучом точки г0).

3.3.4.    Вычисляют 6, (гп ) в произвольной точке гп каждого из отрезков tn0, Й) (п =

■=1, 2, ..., N) относительно всех поверхностей ut (г')=0 (i= 1, 2, 1тах), используя соотношение (1) и рекуррентные соотношения:

2 Зак. 1378

Стр. 4 ГОСТ 25645.204-83

(.V)

8i(rn+l) =

i 1)-")Апах> ^    1,

где in — номер пересекаемой лучом поверхности.

3.3.5.    Из полученных индексов неопределенности для точки гп отбирают относящиеся к /(-й зоне и формируют совокупность векторов неопределенности а*(г„) (/С=1, .... Ктйх).

3.3.6.    Определяют последовательность номеров зон Кп (r0, £2) (л=1, ..., N0), в которых располо

жены отрезки луча tn0, £2), путем отыскания одинаковых (равных) векторов среди а* (г^) и ак(гп) (К= 1..... /Стах)-    Отсутствие таких векторов для некоторой точки rNa (N0*cN) свидетельст

вует о ее расположении вне объекта и процесс идентификации отрезков t п (r0, Q) для ti>N0 прекращают.

3.3.7.    Вычисляют количество вещества на пути луча | (r0t £2) по формуле

«Л ®= 2 /л(^,2)-ркЙ, 2).    (б)

1 4 7

3.4.    Функции самоэкранированности w\ (gi, г0) и экранированности защитой w2 (g2, г0) долж-

ны быть рассчитаны аналогично w (g, г0), причем для вычисления igi (r0, £2) и g2 (r0, £2) осуществляют раздельное суммирование расстояний, пройденных лучом в зонах фантома и защиты, умноженных на плотность вещества в соответствующих зонах.

3.5.    Возможный способ реализации алгоритма приведен в рекомендуемом приложении.

ПРИЛОЖЕНИЕ

Рекомендуемое

СПОСОБ РЕАЛИЗАЦИИ АЛГОРИТМА (ПРОГРАММА OPTIC)

1. Описание программы OPTIC

1JA Программа OPTIC предназначена для^расчета функций экранированности w (g, r0)t самоэкранированности

оч (lu Го) и экранированности защитой w2 (ga, г0) точек г0 в объектах сложной геометрической конфигурации с непостоянными физическими свойствами по объему. Вычисление функции экранированности, определяемой выражением (2) настоящего стандарта, осуществляется методом Монте-Карло. Программа написана на языке Фортран и ориентирована на ЭВМ типа ЕС или БЭСМ-6. Время счета одного варианта зависит от сложности объекта и требуемой точности вычисляемого функционала. Блок-схема вычисления функционалов и хю2 представлена на черт. 1 (в левом углу блоков указаны номера, соответствующие пунктам описания программы).

Передача информации между подпрограммами и связь их с управляющей программой осуществляется в виде описания COMMON-областей и путем задания формальных параметров.

1.2. Описание COMM ON - областей

1.2.1.    COMMON /AG/ UR, VV, WR, A (50, 10), RO (30), где UR, VV, WR —* рабочие ячейки;

A (50, 10)—массив коэффициентов, описывающих поверхности (задается в соответствии с таблицей настоящего стандарта;

RO (30) — плотность вещества в зоне

1.2.2.    COMMON AG1/N, NZON, IT (50), NCF (50), MI (30), IPZ (30,6), ID (30,6), KPN (30), где N<50 — количество поверхностей, применяемое для задания объекта;

NZON<30 — количество зон, применяемое для задания объекта (включая пустоты);

IT (1)<14, 1=1, . . . , N — тип поверхности;

NCF (I) <10, 1= 1,    .    .    ,    N    —максимальное количество коэффициентов, необходимое для задания поверхности I;

М1(К)<6, К= 1, . . , NZON — количество поверхностей, ограничивающих зону К;

IPZ (К, J)<50— порядковый номер поверхности, ограничивающей К-ю зону (K=l,..., NZON;    MI(K));

KPN (К) — индекс материала в зоне К (предполагается, что индекс KPN=li, имеет вещество фантома).

1.2.3.    COMMON /AG2/ IDI (100), IP (100), SP (100), КР (100), где IDI (100) — массив рабочих ячеек;

Г6СТ 25645.204—S3 Стр. 5

IP (ЮО) — массив порядковых номеров поверхностей, пересекаемых прямой в направлении Q, в порядке очередности;

—►    —^

SP (100) —массив расстояний от точки г0, расположенной в объекте в направлении Й, до пересечения с поверхностями в порядке возрастания (SP (1) = 0)i;    '

Блок-схема вычисления функций экранированности wt самоэкранированности и экранированности защитой w2

Черт 1

КР (100) — массив индексов материалов, пересекаемых лучом, в порядке очередности (КР(1) — индекс материала в зоне, содержащей точку г0);

КР = 0— признак выхода из объекта.

I31 Входная информация о геометрии объекта

Входная информация о геометрии объекта считывается с перфокарт и содержится COMMON-областях, описанных в пп. 1.2.1 и 1.2.2.

В данной версии программы предусмотрено использование не более 50 поверхностей 1 и 2-го порядка (задаваемых в соответствии с таблицей настоящего стандарта) для описания геометрии объекта. Максимальное количество зон не превышает 30, причем, каждая зона должна быть ограничена не более, чем шестью поверхностями. Все расстояния задаются в сантиметрах, плотность вещества в зоне — в граммах на кубический сантиметр. При необходимости расширить число зон и поверхностей для описания объекта необходимо изменить соответствующие размерности в COMMON-областях

1.4.    ISTOR — число истории, необходимое для расчетафункций экранированности (рекомендуемое значение ISTOR >10000).

1.5.    г0={Х0, Y0, Z0} — декартовы координаты точки г0.

1.6.    Подпрограмма ROMEGA (U0, V0,    W0)    —    подпрограмма    для    розыгрыша    случайного    направления    вектора

^ ^

Q, имеющего изотропное распределение; ПО, V0, W0 — направляющие косинусы вектора Q в декартовой системе координат. Подпрограмма использует датчик случайных чисел, равномерно распределенных на участке (0,1).

1.7.    Подпрограмма GEOMIN (Х0, Y0, ZO, U0,    V0,    W0,    L1) — основной модуль программы, предназначенный для

, —*• вычисления расстояний от внутренней точки объекта г0= {Х0,    Y0,    Z0}    в    направлении    Q    =    |U0,    V0,    W0)

до пересечения с поверхностями, описывающими объект, а также идентификация материалов, пересекаемых при этом лучом.

Стр. 6 ГОСТ 25645.204-83

Выходная информация содержится в COMMON-области, описанной в п. 1.2 3, и параметре L1.

L1 — максимальное количество пересечений (плюс 1) луча с поверхностями до выхода из объекта (КР (Ll)=0).

1.8. Вычисление толщин вещества фантома ^ и защиты £2 осуществляется раздельным суммированием расстояний,

вещества в соответствующих зонах

цая сетка

разбиения для

{£/} :

A£/= 1 -

0<£< 10

1=1...

.,10

де,-2.

10ч<£<20

/=11,.

..,15

Ag/=5.

20<g<100

/=16,.

..,31

о

7

10 <£<290

/=32,.

..,50

Все случаи, когда £>290, фиксируются в накопителе / = 51. При попадании £ (г0, Q) в соответствующий интервал А£/ 1.10. Конечные функционалы получаются делением величин рину интервала Д£/.

На печать выдаются распределения    и

ности:


в накопитель информации добавляется 1. р(П


на число историй (ISTOR) и , 50


соответствующую ши-а также соответствующие величины вероят-


пройденных лучом в фантоме и защите в направлении Q, умноженным на плотность 19. Анализируется попадание величин |ь £3 и £=£[ + £3 в заданные интервал В программе используется еле;

Д8/;    Р^==и|[^-Д£/ и Р(^=и^г* ■ Д£;.

1.11. Пакет программ содержит все необходимые для проведения расчетов вспомогательные подпрограммы, включая датчик случайных чисел для ЕС ЭВМ (подпрограмма RANDU). Для проведения расчетов на ЭВМ БЭСМ-6 рекомендуется использовать генератор случайных чисел RNDN (библиотечная программа мониторной системы «Дубна»). В ©том случае необходимо заменить функцию RANNO на следующую:

FUNCTION RANNO (NMB)

RANNO=RNDM (—1)

RETURN

END

ГОСТ 25645,204—83 Стр. 7

2. Инструкция по вводу исходных данных

№ и/к

Считываемый символ

Формат

Назначение символа

1

/

N, NZON, NMAT

313

N — число поверхностей;

NZON — число зон;

NMAT — число разных материалов

2—1

2—2

2-3

(IT (I), 1 = 1, N)

2413

1Т *—тип г-й поверхности (в соответствии с таблицей настоящего стандарта) в порядке нумерации поверхностей. При N>24 заносить данные на карты 2—2 и 2—3

3—1

3-2

3-3

(NCF (I), 1*1, N)

2413

NCF,- — число вводимых коэффициентов (в соответствии с таблицей настоящего стандарта) в порядке нумерации поверхностей. При N>24 заносить данные на карты 3—2

и 3—3'

4—1

4—2

(MI(K), K=l, NZON)

2413

Число поверхностей, ограничивающих К-ю зону в порядке принятой нумерации зон. При NZON >24 заносить данные на карту 4—2

5-1

5—NZON

{(IPZ, (K, J), J—1,6), K=l, NZON)

613

Векторы поверхности — номера поверхностей (в принятой нумерации), ограничивающих К-ю зону Требуется ввести п/к с 5—] до 5—NZON

—^ ^

6—1

((ID (K, J), J= 1,6), K=l, NZON)

613

Векторы неопределенности. а^( гк) — индексы неопределенности внутренней точки зоны К относительно ограничивающих ее поверхностей.

Последовательность номеров поверхностей при описании

зоны задается вектором .

Требуется ввести п/к с 6—1 до 6^NZON

6— NZON

7- 1

7—N

((A (1, J), J=l, NCF (I)), 1=1, N)

6E10.0

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

Требуется ввести п/к с 7—1 до 7—N.

Для поверхности общего вида, содержащей более 6 коэффициентов, информация вводится на 2-х п/к, т. е. 7—i—1 и 7—i—2.

8—1

8—2

(KPN (K), K=l, NZON)

2413

Номер материала, расположенного в /С-й зоне, в соответствии с принятой нумерацией.

Значение KPN=1 принято для вещества фантома При NZON>24 заносить данные на п/к 8—2

9-1

9-5

(RO (I), 1=1, NMAT)

6E10.0

Плотность вещества в соответствии с принятой нумерацией (см п/к 8).

При числе различных веществ, большем 6, данные заносить на п/к 9—1, 9—2, . . .

10

ISTOR

16

Число историй, необходимое для оценки интегралов ш, wu w2 (рекомендуется > 10000)

11-1

X0, YO, ZO

3E10.0

—►

Координаты точки г0, в которой производится расчет функций Wt и WТребуется столько п/к, сколько вариантов расчета для

разных точек г0 в данном объекте