Купить ГОСТ 25645.204-83 — бумажный документ с голограммой и синими печатями. подробнее
Распространяем нормативную документацию с 1999 года. Пробиваем чеки, платим налоги, принимаем к оплате все законные формы платежей без дополнительных процентов. Наши клиенты защищены Законом. ООО "ЦНТИ Нормоконтроль"
Наши цены ниже, чем в других местах, потому что мы работаем напрямую с поставщиками документов.
Устанавливает требования к заданию объекта и алгоритм вычисления функций, характеризующих экранированность точек внутри объекта-фантома с окружающей его защитой.
1 Общие положения
2 Задание объекта
3 Алгоритм расчета функции экранированности
Приложение Способ реализации алгоритма (программа OPTIC)
Дата введения | 01.01.1985 |
---|---|
Добавлен в базу | 01.09.2013 |
Актуализация | 01.01.2021 |
20.12.1983 | Утвержден | Госстандарт СССР | 6360 |
---|---|---|---|
Издан | Издательство стандартов | 1984 г. |
Чтобы бесплатно скачать этот документ в формате PDF, поддержите наш сайт и нажмите кнопку:
ГОСУДАРСТВЕННЫЙ СТАНДАРТ СОЮЗА ССР
БЕЗОПАСНОСТЬ РАДИАЦИОННАЯ ЭКИПАЖА КОСМИЧЕСКОГО АППАРАТА В КОСМИЧЕСКОМ ПОЛЕТЕ
Цена 10 коп.
Издание официальное
ГОСУДАРСТВЕННЫЙ КОМИТЕТ СССР ПО СТАНДАРТАМ Москва
П. А. Барсов; А .И. Григорьев, д-р мед. наук; Е. Е. Ковалев, д-р техн. наук; Л. М. Коварскйй, канд. техн. наук; Е. И. Кудряшов, канд. техн. наук; Е. Н. Лесновский, канд. техн. наук; В. А. Панин; Н. М. Пинчук; И. Я. Ремизов, канд. техн. наук; В. А. Сакович, канд. техн. наук; В. М. Сахаров, канд. техн. наук; В. Б. Хвостов, канд. физ.-мат. наук
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.c 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• |
0СЛ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 |
0Г,пз
Щ04
проб 0006 0007 ОЭС 8 000*9
0001
0001
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
у Ъ )
5°- М г* ,41*М П|?! ш*м I и и** с
j^s»* s-:aac n »t )***s*x- ?*ac и , -u*si
?^*S'
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>*
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—ЙЛ
|
***!№***$£*****#*******$*** $~*******ф**эМ *******************************^&***************:к***%*****?*«*£ |
УДК 62—784.7:539.16:629.78:006.354
Группа Ф40
Soacecrew radiation safety during spaceflight. Computation methods of points shielding inside fantom
ОКП 696800
Постановлением Государственного комитета СССР по стандартам от 20 декабря 1983 г. № 6360 срок введения установлен
с 01.01.85
Настоящий стандарт устанавливает требования к заданию объекта и алгоритм вычисления функций, характеризующих экранированность точек внутри объекта-фантома с окружающей его защитой.
Под защитой в стандарте понимают конструкцию космического аппарата (КА), его оборудование и специальное снаряжение, защищающее (экранирующее) космонавта от ионизирующего излучения.
Стандарт предназначен для подготовки исходных данных, необходимых при расчетах на предприятиях и организациях, занимающихся научно-исследовательскими и опытно-конструкторскими работами, связанными с обеспечением радиационной безопасности экипажа космического аппарата в космическом полете.
1.1. Экранированность точки г0) расположенной внутри фантома, характеризует функция экра
нированное™ w (|, го) такая, что w (g, Го) d£ представляет вероятность для лучей, изотропно испущенных из точки го, встретить на своем пути суммарное количество вещества фантома и защиты \ в интервале от g до g + выраженное в массовых единицах длины.
где — количество вещества фантома;
|2— количество вещества защиты.
1.2. Под массовой единицей длины в веществе понимают произведение линейной единицы длины на плотность вещества.
—►
1.3. Самоэкранированность точки г0, расположенной внутри фантома, характеризует функция еамоэкранированности w\ (£1, г0), тождественно равная w (g, г0) $8=<ь
1.4. Экранированность защитой точки го, расположенной внутри фантома, характеризует функция экранированности защитой w2 (g, /о), тождественно равная w (£, г0)е1==о .
2.1. Объект, в виде выпуклого тела, задают совокупностью зон с постоянными физическими свойствами вещества в пределах зоны. Каждой зоне присваивают номер К— 1, 2, ..., Ктах> где Ктах — максимальное количество зон, необходимое для задания объекта.
Примечание. Если исходный объект представляет собой вогнутое тело, то его следует дополнить пустыми зонами.
Издание официально* Перепечатка воспрещена
★
© Издательство стандартов, 1984
ГОСТ 25645.204-83 Стр. 19
ЧИСЛО ИСТОРИЙ * 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- + C2 “1 |
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
ют признаком, именуемым индексом неопределенности б/ (г*), значение которого определяется выражением
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. Располагают полученный массив значений Sn (г0, й) (я = 1, ..., N) в порядке возрастания,
—> —>
формируя при этом последовательность соответствующих номеров поверхностей i п (г0, Й).
—> —у
3.3.3. Вычисляют длины отрезков t п (го, й) между последовательными пересечениями
^п(го» 2) = Sn(rо, Й) Sn-\(rо, 2), (4)
^ __ —>
положив S0 (ль П) — 0 (пересечение лучом точки г0).
3.3.4. Вычисляют 6, (гп ) в произвольной точке гп каждого из отрезков tn (г0, Й) (п =
■=1, 2, ..., N) относительно всех поверхностей ut (г')=0 (i= 1, 2, 1тах), используя соотношение (1) и рекуррентные соотношения:
2 Зак. 1378
(.V)
8i(rn+l) =
i 1)-")Апах> ^ 1,
где in — номер пересекаемой лучом поверхности.
3.3.5. Из полученных индексов неопределенности для точки гп отбирают относящиеся к /(-й зоне и формируют совокупность векторов неопределенности а*(г„) (/С=1, .... Ктйх).
3.3.6. Определяют последовательность номеров зон Кп (r0, £2) (л=1, ..., N0), в которых располо
жены отрезки луча tn (г0, £2), путем отыскания одинаковых (равных) векторов среди а* (г^) и ак(гп) (К= 1..... /Стах)- Отсутствие таких векторов для некоторой точки rNa (N0*cN) свидетельст
вует о ее расположении вне объекта и процесс идентификации отрезков t п (r0, Q) для ti>N0 прекращают.
3.3.7. Вычисляют количество вещества на пути луча | (r0t £2) по формуле
3.4. Функции самоэкранированности w\ (gi, г0) и экранированности защитой w2 (g2, г0) долж-
ны быть рассчитаны аналогично w (g, г0), причем для вычисления igi (r0, £2) и g2 (r0, £2) осуществляют раздельное суммирование расстояний, пройденных лучом в зонах фантома и защиты, умноженных на плотность вещества в соответствующих зонах.
3.5. Возможный способ реализации алгоритма приведен в рекомендуемом приложении.
ПРИЛОЖЕНИЕ
Рекомендуемое
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) — массив рабочих ячеек;
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)
до пересечения с поверхностями, описывающими объект, а также идентификация материалов, пересекаемых при этом лучом.
Выходная информация содержится в 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
2. Инструкция по вводу исходных данных | ||||||||||||||||||||||||||||||||||||||||||||||||
|