Tải bản đầy đủ - 0 (trang)
Tài liệu tham khảo

Tài liệu tham khảo

Tải bản đầy đủ - 0trang

Đồ án tốt nghiệp



KTHN & VLMT – K54



PHỤ LỤC

Phụ lục 1: Bảng dữ liệu đầu vào của chương trình tính tốn hiệu

suất:



STT



Ngày



1

2

3

4

5

6

7

8

9

10

11

12

13

14

15

16

17

18

19

20

21

22

23

24

25

26

27

28

29

30



14/2

18/2

19/2

20/2

21/2

25/2

26/2

27/2

28/2

4/3

5/3

6/3

7/3

11/3

12/3

13/3

14/3

18/3

19/3

20/3

25/3

26/3

27/3

28/3

1/4

2/4

3/4

4/4

8/4

15/4



Cường đợ chùm

tia

(µA)

34

36

34

35

35

35

34

36

20

35

30

34

31

34

24

33

28

35

33

35

36

34

34

36

34

32

34

33

31

35



Đinh Quang Huy



Cường đợ tích

phân

(µA)

11.5

25.5

10.1

24.1

12.5

25.6

19.7

21.8

10.2

19.0

11.5

25.6

21.5

25.1

12.2

25.0

10.0

18.1

11.0

21.7

25.1

16.0

24.2

17.2

24.9

11.9

23.6

11.1

17.5

19.0



Thời gian

tổng hợp

(phút)

37

43

32

33

33

39

35

35

37

33

43

35

31

33

50

35

37

39

44

38

38

36

44

34

32

41

36

40

38

41



Hoạt độ

18F-FDG

(mCi)

467.621

831.081

451.554

831.081

548.513

903.108

781.216

692.567

383.783

770.135

447.108

991.756

830.135

775.675

455.986

880.945

298.081

635.000

311.918

736.891

858.783

410.540

831.081

506.405

692.567

468.729

747.973

445.459

753.513

714.729

Trang 45



Đồ án tốt nghiệp



KTHN & VLMT – K54



Phụ lục 2: Code chương trình tính tốn hiệu suất:

Imports Microsoft.VisualBasic.FileIO

Imports System.IO

Public Class Form1

Private Sub Form1_Load(sender As Object, evt As EventArgs) Handles

MyBase.Load

AddHandler BeamCurrentTbx.TextChanged, AddressOf

BeamCurrentTbx_TextChanged

AddHandler IrradiationTimeTbx.TextChanged, AddressOf

IrradiationTimeTbx_TextChanged

AddHandler IntergratedCurrentTbx.TextChanged, AddressOf

IntergratedCurrentTbx_TextChanged

End Sub

Private Sub BeamCurrentTbx_TextChanged(sender As Object, evt As

EventArgs) 'Handles BeamCurrentTbx.TextChanged

'Dim IntergratedCurrent As Double

'If (BeamCurrentTbx.Text = "") Then

' IntergratedCurrent = 0

'Else

' IntergratedCurrent = CDbl(Beam.CurrentTbx.Text) *

CDbl(IrradiationTimeTbx.Text) / 60.0#

'End If

'IntergratedCurrentTbx.Text = Format(IntergratedCurrent, "###0.0")

Dim IrradiationTime As Double

If BeamCurrentTbx.Text = "" Then

IrradiationTime = 0



Đinh Quang Huy



Trang 46



Đồ án tốt nghiệp



KTHN & VLMT – K54



Else

IrradiationTime = 60.0 * CDbl(IntergratedCurrentTbx.Text) /

CDbl(BeamCurrentTbx.Text)

End If

IrradiationTimeTbx.Text = Format(IrradiationTime, "###0.0")

End Sub

Private Sub IrradiationTimeTbx_TextChanged(sender As Object, evt As

EventArgs) 'Handles IrradiationTimeTbx.TextChanged

'Dim IntergratedCurrent As Double

'If (IrradiationTimeTbx.Text = "") Then

' IntergratedCurrent = 0

'Else

' IntergratedCurrent = CDbl(BeamCurrentTbx.Text) *

CDbl(IrradiationTimeTbx.Text) / 60.0#

'End If

'IntergratedCurrentTbx.Text = Format(IntergratedCurrent, "###0.0")

End Sub

Private Sub IntergratedCurrentTbx_TextChanged(ByVal sender As

System.Object, ByVal e As System.EventArgs) 'Handles

IntergratedCurrentTbx.TextChanged

Dim IrradiationTime As Double

If IntergratedCurrentTbx.Text = "" Then

IrradiationTime = 0

Else

IrradiationTime = 60.0 * CDbl(IntergratedCurrentTbx.Text) /

CDbl(BeamCurrentTbx.Text)

End If

IrradiationTimeTbx.Text = Format(IrradiationTime, "###0.0")

End Sub



Đinh Quang Huy



Trang 47



Đồ án tốt nghiệp



KTHN & VLMT – K54



Private Sub CloseBtn_Click(sender As Object, evt As EventArgs) Handles

CloseBtn.Click

Close()

End Sub

Private Sub CalculateBtn_Click(sender As Object, evt As EventArgs) Handles

CalculateBtn.Click

Const Na = 6.02E+23 ''Avogadro's Number

Const M = 20.0# 'M(H2O)

Const T12 = 110.0# 'm

Const v = 0.95 '%H2O

Const rho = 1.0# 'g/cm3

Dim E0 As Double = BeamEnergyTbx.Text 'TxtBeaEne.Text

Dim Ib As Double = BeamCurrentTbx.Text 'TxtBeaInt.Text

Dim TEOB As Double = IrradiationTimeTbx.Text 'm TxtIrrTime

Dim Tsyn As Double = SynthizedTimeTbx.Text 'TxtSynTime.Text

Dim AFDG As Double = Activity18FTbx.Text 'Txt18FDG

Dim lamda = Math.Log(2) / T12

Dim E() As Double = Nothing

Dim Cs() As Double = Nothing

Dim CsCount As Integer

If Not Read_Cs(E0, E, Cs, CsCount) Then Exit Sub

Dim F(CsCount) As Double 'Hoạt đợ sau quá trình tổng hợp ( khơng tính đến

hiệu suất tổng hợp )

For i = 0 To CsCount - 1

F(i) = Ib / 0.0000000000001601 * rho * v * Na / M * Cs(i) * 1.0E-24 '1barn

= 1.0E-24 cm2

Next

Dim R = 0.0#



Đinh Quang Huy



Trang 48



Đồ án tốt nghiệp



KTHN & VLMT – K54



''Integrated

'x0=Range(E(0))

'xn=Range(E(CsCount-1))=Range (17.3MeV)

For i = CsCount - 1 To 1 Step -1

R = R + (Range(E(i)) - Range(E(i - 1))) * (F(i) + F(i - 1)) / 2 'average value

Next

Dim AEOB = R * (1 - Math.Exp(-lamda * TEOB)) / 37000000000.0 * 1000.0

'mCi '37E9

Dim Af = AEOB * Math.Exp(-lamda * Tsyn)

Dim Eff = AFDG / Af * 100.0#

SynthesisEfficiencyTbx.Text = Format(Eff, "00.0") 'LblEff

CalculatedActivityTbx.Text = Format(Af, "####0.0") 'LblAct

End Sub

Private Shared Function Read_Cs(Eb As Double, ByRef E() As Double, ByRef

Cs() As Double, ByRef CsCount As Integer) As Boolean

Dim capacity = 100

ReDim E(capacity)

ReDim Cs(capacity)

Dim itemCount = 0

Dim E1 As Double

Dim E2 As Double

Dim Cs1 As Double

Dim Cs2 As Double

Dim dum As String = String.Empty

FileOpen(1, "Cross.txt", _

OpenMode.Input, OpenAccess.Read, OpenShare.LockWrite)

Input(1, E1)

Input(1, Cs1)

Input(1, dum)



Đinh Quang Huy



Trang 49



Đồ án tốt nghiệp



KTHN & VLMT – K54



Input(1, dum)

If (Eb < E1) Then

E = Nothing

Cs = Nothing

CsCount = 0

FileClose(1)

MsgBox("Data is out of range")

Return False

End If

Do While Not EOF(1)

Input(1, E2)

Input(1, Cs2)

Input(1, dum)

Input(1, dum)

If (Eb < E2) Then

E(itemCount) = Eb

Cs(itemCount) = LinearInterpolate(Eb, E1, E2, Cs1, Cs2)

itemCount = itemCount + 1

Exit Do

Else

E(itemCount) = E2

Cs(itemCount) = Cs2

itemCount = itemCount + 1

If Eb = E2 Then Exit Do

End If

E1 = E2

Cs1 = Cs2

If itemCount = capacity Then

capacity = capacity + 100



Đinh Quang Huy



Trang 50



Đồ án tốt nghiệp



KTHN & VLMT – K54



ReDim Preserve E(capacity)

ReDim Preserve Cs(capacity)

End If

Loop

FileClose(1)

CsCount = itemCount

ReDim Preserve E(CsCount)

ReDim Preserve Cs(CsCount)

Return True

End Function

'Nội suy tuyến tính

Private Shared Function LinearInterpolate(x As Double, x0 As Double, x1 As

Double, y0 As Double, y1 As Double) As Double

Return (x - x0) * (y1 - y0) / (x1 - x0) + y0

End Function

Private Shared Function Range(ByVal x As Double) As Double 'cm

Const A = 0.0018830001

Const b = 0.012243578

Const c = 0.015855191

Const d = -0.0016302767

Const er = -0.0019644942

'(A+bx+cx^2+d(x^2)Ln(x)+er(e^-x))/10

Return (A + b * x + c * x * x + d * x * x * Math.Log(x) + er * Math.Exp(-x)) /

10.0

End Function

'Private Function Eff_Thick(Ex As Double, dEx As Double) As Double

' Return Range(Ex) - Range(Ex - dEx)

'End Function



Đinh Quang Huy



Trang 51



Đồ án tốt nghiệp



KTHN & VLMT – K54



End Class



Phụ lục 3: Tiết diện phản ứng (theo thực nghiệm của S.Takacs,2003)

Năng

lượng

(MeV)

2.5

2.6

2.7

2.8

2.9

3

3.1

3.2

3.3

3.4

3.5

3.6

3.7

3.8

3.9

4

4.1

4.2

4.3

4.4

4.5

4.6

4.7

4.8

4.9

5

5.1

5.2

5.3

5.4

5.5

5.6

5.7

5.8

5.9

6



Tiết

diện

(barn)

0.0027

0.0523

0.0066

0.0079

0.0059

0.088

0.0464

0.0623

0.0618

0.0405

0.0469

0.043

0.1161

0.1171

0.1446

0.1568

0.2481

0.1832

0.3628

0.2843

0.2202

0.1348

0.1621

0.2616

0.4224

0.5024

0.498

0.465

0.4063

0.3523

0.3441

0.3753

0.2849

0.2358

0.2688

0.4654



Đinh Quang Huy



6.1

6.2

Năng

lượng

(MeV)

6.3

6.4

6.5

6.6

6.7

6.8

6.9

7

7.1

7.2

7.3

7.4

7.5

7.6

7.7

7.8

7.9

8

8.1

8.2

8.3

8.4

8.5

8.6

8.7

8.8

8.9

9

9.1

9.2

9.3

9.4

9.5

9.6



0.3027

0.2276

Tiết

diện

(barn)

0.2183

0.2474

0.2615

0.2657

0.2628

0.2343

0.2216

0.2749

0.3423

0.3404

0.3251

0.3057

0.2496

0.2255

0.2365

0.2435

0.2478

0.2466

0.272

0.2814

0.2507

0.2373

0.242

0.2575

0.2666

0.2481

0.2302

0.2176

0.2042

0.2038

0.2052

0.1899

0.1772

0.1862



9.7

9.8

9.9

10

Năng

lượng

(MeV)

10.1

10.2

10.3

10.4

10.5

10.6

10.7

10.8

10.9

11

11.1

11.2

11.3

11.4

11.5

11.6

11.7

11.8

11.9

12

12.1

12.2

12.3

12.4

12.5

12.6

12.7

12.8

12.9

13

13.1

13.2



0.2003

0.1976

0.1949

0.1921

Tiết

diện

(barn)

0.1894

0.1867

0.184

0.1812

0.1788

0.1766

0.1744

0.1722

0.17

0.1678

0.1656

0.1635

0.1613

0.1591

0.1569

0.1548

0.1526

0.1505

0.1483

0.1462

0.144

0.1419

0.1398

0.1377

0.1356

0.1335

0.1314

0.1293

0.1273

0.1252

0.1232

0.1212

Trang 52



Đồ án tốt nghiệp



13.3

13.4

13.5

13.6

13.7

13.8

13.9

14

14.1

14.2

14.3

14.4

14.5

14.6

14.7

14.8

14.9

15

15.1

15.2

15.3

15.4

15.5

15.6

15.7

15.8

15.9

16

16.1

16.2

16.3

16.4

16.5

16.6

16.7

16.8

16.9

17

17.1

17.2

17.3

17.4



0.1191

0.1171

0.1152

0.1132

0.1112

0.1093

0.1074

0.1054

0.1035

0.1017

0.0998

0.098

0.0961

0.0943

0.0926

0.0908

0.0891

0.0873

0.0856

0.084

0.0823

0.0807

0.0791

0.0775

0.0759

0.0744

0.0729

0.0714

0.0699

0.0685

0.0671

0.0657

0.0643

0.063

0.0617

0.0604

0.0592

0.0579

0.0567

0.0556

0.0544

0.0533



Đinh Quang Huy



KTHN & VLMT – K54



17.5

17.6

17.7

17.8

17.9

18

18.1

18.2

18.3

18.4

18.5

18.6

18.7

18.8

18.9

19

19.1

19.2

19.3

19.4

19.5

19.6

19.7

19.8

19.9

20

20.1

20.2

20.3

20.4

20.5

20.6

20.7

20.8

20.9

21

21.1

21.2

21.3

21.4

21.5

21.6



0.0522

0.0512

0.0501

0.0491

0.0482

0.0472

0.0463

0.0454

0.0445

0.0437

0.0428

0.0421

0.0413

0.0405

0.0398

0.0391

0.0384

0.0378

0.0371

0.0365

0.036

0.0354

0.0349

0.0344

0.0339

0.0335

0.0331

0.0327

0.0323

0.0319

0.0315

0.0312

0.0309

0.0306

0.0303

0.03

0.0297

0.0295

0.0292

0.029

0.0287

0.0285



21.7

21.8

21.9

22

22.1

22.2

22.3

22.4

22.5

22.6

22.7

22.8

22.9

23

23.1

23.2

23.3

23.4

23.5

23.6

23.7

23.8

23.9

24

24.1

24.2

24.3

24.4

24.5

24.6

24.7

24.8

24.9

25

25.1

25.2

25.3

25.4

25.5

25.6

25.7

25.8



0.0283

0.0281

0.0279

0.0277

0.0276

0.0274

0.0272

0.027

0.0269

0.0267

0.0266

0.0264

0.0263

0.0261

0.026

0.0258

0.0257

0.0256

0.0254

0.0253

0.0251

0.025

0.0249

0.0248

0.0246

0.0245

0.0244

0.0243

0.0242

0.024

0.0239

0.0238

0.0237

0.0236

0.0235

0.0234

0.0233

0.0233

0.0232

0.0231

0.023

0.0229



Trang 53



Đồ án tốt nghiệp



KTHN & VLMT – K54



25.9



26



0.0228



Đinh Quang Huy



0.0228



Trang 54



Đồ án tốt nghiệp



26.1

26.2

26.3

26.4

26.5

26.6

26.7

26.8

26.9

27

27.1

27.2

27.3

27.4

27.5

27.6

27.7

27.8

27.9

28

28.1

28.2

28.3

28.4

28.5

28.6

28.7

28.8

28.9

29

29.1

29.2

29.3

29.4

29.5

29.6

29.7

29.8



KTHN & VLMT – K54



0.0227

0.0226

0.0225

0.0225

0.0224

0.0223

0.0222

0.0222

0.0221

0.0221

0.022

0.022

0.0219

0.0219

0.0219

0.0218

0.0218

0.0218

0.0218

0.0218

0.0219

0.0219

0.022

0.022

0.0221

0.0222

0.0223

0.0224

0.0225

0.0226

0.0227

0.0228

0.023

0.0231

0.0233

0.0234

0.0236

0.0237



Đinh Quang Huy



Trang 55



Tài liệu bạn tìm kiếm đã sẵn sàng tải về

Tài liệu tham khảo

Tải bản đầy đủ ngay(0 tr)

×