Скачиваний:
8
Добавлен:
15.03.2019
Размер:
296.5 Кб
Скачать

 

я

ц

А

Ы

Ч

ц «Э

я

»

Ч

 

№4

:

. . - 15-1-3___________________

.

:

___________________ Ч

. .

Ц :

.

: -

; -

.

я

я

:

 

 

 

I 0=δI 0 –

 

 

 

 

 

 

 

 

 

,

/ 2; δ –

,

(T0)^4

 

 

 

5.67 10−8 ,

 

 

0.95; –

-

,

,

/( 2 · 4 ); T0 –

, .

 

,

 

 

 

 

 

 

M=−LN( I0/S0 ) ,

M –

 

; S0 –

,

1366 / 2 .

A

 

,

z

 

 

 

 

A= A NA Δz

 

 

 

 

; A –

, .

 

, 2 ; NA –

 

,

-3; Δz –

 

 

Э A=2 π

RA^2

r A –

, .

= M + A

I=S0 exp(− )

/

T=

.

ΔT=T−T0

,

Python,

 

,

 

1.

#!/*sr/bin/env python

# coding: utf8

#from __future__ import unicode_literals *

from numpy import *

from math import *

from scipy import *

from matplotlib.pyplot import *

ra = 1E-7

T0 = 300.0

Na = 10**8

sigma = 5.67E-8

S0 = 1366.0

b = 0.95

z = 2000

Namin = 1E+8

Namax = 1E+9

dNa = 1E+8

m = int((Namax-Namin)/dNa+1)

Na = zeros(m)

ta = zeros((len(Na)), dtype=float)

t = zeros((len(Na)), dtype=float)

I = zeros((len(Na)), dtype=float)

T = zeros((len(Na)), dtype=float)

dT =zeros((len(Na)), dtype=float)

for j in range(m):

Na[j] = Namax - j*dNa

I0 = b*sigma*(T0**4)

tm = -log(I0 / S0)

qa = 2*pi*(ra**2)

ta[j] = qa*Na[j]*z

t[j] = ta[j]+tm

I[j] = S0*exp(-t[j])

T[j] = (I[j]/(b*sigma))**(1.0/4.0)

dT[j] = T[j]-T0

print I0, tm, t, Na

semilogx(Na, t, color="blue", lw=1)

xlabel("Na,m", fontsize=12)

ylabel("t", fontsize=12)

grid(True)

show()

,

:

>>>

436.3065 1.14129706211 [ 1.26696077 1.2543944 1.24182803 1.22926166 1.21669529 1.20412892

1.19156254 1.17899617 1.1664298 1.15386343] [ 1.00000000e+09 9.00000000e+08 8.00000000e+08 7.00000000e+08

6.00000000e+08 5.00000000e+08 4.00000000e+08 3.00000000e+08

2.00000000e+08 1.00000000e+08] 6.28318530718e-14 [ 0.12566371 0.11309734 0.10053096 0.08796459 0.07539822 0.06283185

0.05026548 0.03769911 0.02513274 0.01256637] [ 290.72172789 291.63649329 292.55413704 293.47466817 294.39809579

295.32442901 296.25367696 297.18584882 298.12095379 299.0590011 ] [-9.27827211 -8.36350671 -7.44586296 -6.52533183 -5.60190421 -4.67557099

-3.74632304 -2.81415118 -1.87904621 -0.9409989 ]

>>>

:

,

:

 

 

>>>

436.3065 1.14129706211 [ 1.26696077 1.2543944 1.24182803 1.22926166 1.21669529 1.20412892

1.19156254 1.17899617 1.1664298 1.15386343] [ 1.00000000e+09 9.00000000e+08 8.00000000e+08 7.00000000e+08

6.00000000e+08 5.00000000e+08 4.00000000e+08 3.00000000e+08

2.00000000e+08 1.00000000e+08] 6.28318530718e-14 [ 0.12566371 0.11309734 0.10053096 0.08796459 0.07539822 0.06283185

0.05026548 0.03769911 0.02513274 0.01256637] [ 290.72172789 291.63649329 292.55413704 293.47466817 294.39809579

295.32442901 296.25367696 297.18584882 298.12095379 299.0590011 ] [-9.27827211 - 8.36350671 -7.44586296 -6.52533183 -5.60190421 -4.67557099

-3.74632304 -2.81415118 -1.87904621 -0.9409989 ]

>>>

:

,

.

:

,

№4

.

;

,

;

, z,

;

,

.

: -

; -

.

,

,

.

Соседние файлы в папке Attachments_ivan.skachkoff@yandex.ru_2017-12-24_20-40-50