Giới thiệu dự án

Trong cấu trúc an toàn của nhà máy điện hạt nhân (NMĐHN) sử dụng lò phản ứng nước áp lực thế hệ mới như VVER-1000 (công suất nhiệt 3000 MWth, tổng công suất phát điện 1060 MWe), bình sinh hơi (Steam Generator - SG) giữ vai trò là thiết bị trao đổi nhiệt trung gian quan trọng phân tách vòng sơ cấp mang phóng xạ và vòng thứ cấp sinh hơi quay tua-bin. Mỗi bình sinh hơi kiểu ngang trang bị hệ thống gồm 10.978 ống trao đổi nhiệt với kích thước danh định $\varnothing 16 \times 1.5\text{ mm}$, chế tạo từ thép không gỉ Austenitic chịu nhiệt và chống ăn mòn cao.

+-----------------------------------------------------------------------------------+
|                        SƠ ĐỒ VÒNG LẶP NHIỆT NMĐHN VVER-1000                       |
+-----------------------------------------------------------------------------------+
|                                                                                   |
|  [ LÕI LÒ PHẢN ỨNG ]                                                              |
|        │                                                                          |
|        ▼ (Nước áp lực cao: 15.7 MPa, 321°C) - "Hot Leg"                           |
|  ┌─────────────────────────────────────────────────────────────────────────────┐  |
|  │ BÌNH SINH HƠI (STEAM GENERATOR - VVER-1000)                                 │  |
|  │  ┌───────────────────────────────────────────────────────────────────────┐  │  |
|  │  │ Bộ thu sơ cấp (Inlet Collector)                                       │  │  |
|  │  │  └─► Bó 10.978 ống trao đổi nhiệt (Austenitic Stainless Steel)        │  │  |
|  │  │        ├─► Môi trường ăn mòn: T=291-321°C, H2O2, OH-, O2, H2          │  │  |
|  │  │        ├─► Ứng suất cơ học & áp suất chu trình (P=6.27 MPa)           │  │  |
|  │  │        └─► HIỆN TƯỢNG: ĂN MÒN ỨNG SUẤT (SCC) & GIÒN HÓA HYDRO (HE)    │  │  |
|  │  │  ┌─► Phân tích SIF (K_I) bằng Giải tích, Bán thực nghiệm & ANSYS FEA  │  │  |
|  │  │ Bộ thu sơ cấp (Outlet Collector)                                      │  │  |
|  │  └───────────────────────────────────────────────────────────────────────┘  │  |
|  │  Vòng thứ cấp: Sinh hơi bão hòa (1470 t/h, 6.27 MPa) ──► [ TUA-BIN PHÁT ĐIỆN ] │  |
|  └─────────────────────────────────────────────────────────────────────────────┘  |
|        │                                                                          |
|        └─◄ (Nước tuần hoàn: 291°C) - "Cold Leg" ◄── [ BƠM NƯỚC SƠ CẤP ]           |
+-----------------------------------------------------------------------------------+

Hệ thống ống trao đổi nhiệt phải vận hành liên tục dưới điều kiện cực hạn: áp suất vận hành $6.27\text{ MPa}$ tại lối ra, nhiệt độ dòng tải nhiệt $291^\circ\text{C} - 321^\circ\text{C}$, thông lượng bức xạ nơtron nhiệt cao và tiếp xúc với các hợp chất oxy hóa phân hủy phóng xạ ($\text{OH}^-$, $\text{H}_2\text{O}_2$, $\text{O}_2$, $\text{H}^+$). Môi trường này thúc đẩy hiện tượng nứt do môi trường ăn mòn và ứng suất (Stress Corrosion Cracking - SCC) kết hợp quá trình giòn hóa do tích tụ Hydro (Hydrogen Embrittlement), đe dọa trực tiếp tính toàn vẹn cấu trúc của ranh giới an toàn hạt nhân.

Problem statement tập trung vào việc thiếu hụt các mô hình định lượng chuẩn xác về trường tập trung ứng suất tại đỉnh vết nứt khi có sự tương tác đồng thời của ứng suất cơ học và nồng độ khuếch tán hydro theo thời gian thực. Nếu không đánh giá chính xác hệ số cường độ ứng suất (Stress Intensity Factor - SIF, ký hiệu $K$), vết nứt tế vi sẽ phát triển mất kiểm soát thành nứt xuyên hạt (Transgranular SCC) hoặc nứt liên kết hạt (Intergranular SCC - IGSCC), dẫn đến nguy cơ rò rỉ chất làm mát mang phóng xạ sơ cấp sang vòng thứ cấp.

Mục tiêu cụ thể của dự án nghiên cứu:

  1. Thiết lập mô hình giải tích khuếch tán và tích tụ nồng độ hydro ($C_{\text{H}_2}$) trong mạng tinh thể thép không gỉ Austenitic theo thời gian dựa trên định luật Fick.
  2. Xây dựng thuật toán tính toán trường ứng suất và hệ số cường độ ứng suất Mode I ($K_I$) tại đỉnh vết nứt bằng nghiệm giải tích Westergaard và các công thức bán thực nghiệm Tada-Paris.
  3. Ứng dụng phương pháp phân tích phần tử hữu hạn (Finite Element Analysis - FEA) trên nền tảng ANSYS Workbench 15.0 để mô phỏng hình học vết nứt, chia lưới phần tử kỳ dị (Singular Mesh) và trích xuất giá trị $K_I$.
  4. Đánh giá sai số tương đối giữa phương pháp mô phỏng số và lý thuyết giải tích, thiết lập ngưỡng an toàn cơ học $K_{ISCC}$ phục vụ công tác bảo trì định kỳ.

Giải pháp kết hợp giữa cơ học phá hủy vật liệu đàn hồi tuyến tính (Linear Elastic Fracture Mechanics - LEFM) và mô phỏng số FEA cho phép kiểm định trực quan các vùng tập trung ứng suất với chi phí thấp hơn $90%$ so với thử nghiệm cơ lý thực nghiệm trong môi trường bức xạ kín. Phạm vi nghiên cứu giới hạn ở bài toán phẳng biến dạng phẳng (Plane Strain) cho vết nứt dạng nứt Mode I (Opening Mode) trên vật liệu thép không gỉ Austenitic tiêu chuẩn của bình sinh hơi VVER-1000.

Phân tích và thiết kế giải pháp

Phân tích hiện trạng

Quá trình đánh giá an toàn kết cấu ống trao đổi nhiệt lò hạt nhân hiện nay dựa trên 3 nhóm phương pháp chính:

Phương pháp Ưu điểm Nhược điểm Giới hạn ứng dụng
Nghiệm giải tích Westergaard thuần túy - Nghiệm chính xác dạng đóng (closed-form)
- Tốc độ tính toán tức thì ($<1\text{ ms}$)		 - Chỉ áp dụng cho tấm phẳng vô hạn
- Không xét được hình học biên phức tạp		 Vết nứt vi mô lý tưởng hóa, giai đoạn sơ bộ
Công thức bán thực nghiệm (Tada-Paris) - Hiệu chỉnh theo hệ số hình học $Y(a/W)$
- Tính toán nhanh, chuẩn hóa ASME		 - Sai số tăng cao khi tỷ lệ nứt $a/W > 0.6$
- Bỏ qua trường nhiệt phi tuyến		 Đánh giá nhanh tại hiện trường nhà máy
Mô phỏng số FEA (ANSYS Workbench 15.0) - Mô phỏng chính xác hình học 2D/3D
- Tích hợp trường tải nhiệt và áp suất đa chiều		 - Yêu cầu cấu hình tính toán cao
- Đòi hỏi kỹ thuật chia lưới kỳ dị chuẩn		 Phân tích chi tiết độ toàn vẹn cấu trúc (RVI)

Yêu cầu kỹ thuật được phân loại theo ma trận MoSCoW:

  • Must have: Mô hình hóa chính xác hàm phân bố nồng độ hydro $C_{\text{H}2}(t)$, giải thuật trường ứng suất $\sigma{xx}, \sigma_{yy}, \tau_{xy}$ quanh đỉnh vết nứt, tính toán SIF ($K_I$) trên ANSYS 15.0.
  • Should have: Đối sánh sai số 3 chiều giữa Westergaard, công thức bán thực nghiệm và ANSYS; phân tích cơ chế gãy vỡ màng thụ động (Film Rupture Model).
  • Could have: Tích hợp mô hình thế năng tương tác nguyên tử Lennard-Jones giải thích liên kết mạng tinh thể lập phương tâm diện (FCC).
  • Won't have: Mô phỏng động lực học dòng chảy 2 pha kèm phân rã phóng xạ 3D toàn phần bình sinh hơi trong phạm vi khóa luận này.
+-----------------------------------------------------------------------------------+
|               QUY TRÌNH PHÂN TÍCH HỆ SỐ CƯỜNG ĐỘ ỨNG SUẤT (SIF)                   |
+-----------------------------------------------------------------------------------+
|                                                                                   |
|  [ĐIỀU KIỆN BIÊN VẬN HÀNH]       [THUỘC TÍNH VẬT LIỆU]    [HÌNH HỌC VẾT NỨT]      |
|  - Áp suất trong: P = 6.27 MPa   - Thép Austenitic        - Chiều dài nứt: a      |
|  - Nhiệt độ: T = 291 - 321°C     - E = 200 GPa, ν = 0.3   - Chiều rộng vách: W    |
|             │                               │                       │             |
|             └───────────────────────┬───────┴───────────────────────┘             |
|                                     ▼                                             |
|             ┌───────────────────────────────────────────────┐                     |
|             │ TIẾP CẬN ĐA PHƯƠNG PHÁP TÍNH TOÁN $K_I$       │                     |
|             └───────────────────────┬───────────────────────┘                     |
|                                     │                                             |
|        ┌────────────────────────────┼────────────────────────────┐                |
|        ▼                            ▼                            ▼                |
|  ┌─────────────────────┐   ┌─────────────────────┐   ┌─────────────────────┐      |
|  │  GIẢI TÍCH KHUẾCH   │   │  CÔNG THỨC BÁN      │   │  MÔ PHỎNG SỐ FEA    │      |
|  │  TÁN VÀ WESTERGAARD │   │  THỰC NGHIỆM TADA   │   │  ANSYS WORKBENCH 15 │      |
|  │  - Fick Law: C_H2(t)│   │  - K_I = σ√(πa)·Y   │   │  - Singular Mesh    │      |
|  │  - σ_ij(r, θ)       │   │  - Bù biên hữu hạn  │   │  - Contour Integral │      |
|  └──────────┬──────────┘   └──────────┬──────────┘   └──────────┬──────────┘      |
|             │                         │                         │                 |
|             └────────────────────────►├────────────────────────◄┘                 |
|                                       ▼                                           |
|             ┌───────────────────────────────────────────────┐                     |
|             │ ĐỐI SÁNH KẾT QUẢ & ĐÁNH GIÁ AN TOÀN SCC       │                     |
|             │ - Kiểm tra sai số tương đối (Sai số < 5%)     │                     |
|             │ - So sánh ngưỡng gãy: K_I < K_ISCC            │                     |
|             └───────────────────────────────────────────────┘                     |
+-----------------------------------------------------------------------------------+

Thiết kế hệ thống mô phỏng

Hệ thống tính toán tích hợp bộ thông số vật lý và hình học tiêu chuẩn:

  • Thông số hình học: Tấm đại diện thành ống trao đổi nhiệt chiều dày $W$, chiều sâu vết nứt đơn biên $a$, bán kính góc đỉnh nứt $r \to 0$.
  • Thông số vật liệu Thép không gỉ Austenitic: Mô đun đàn hồi $E = 2.0 \times 10^5\text{ MPa}$, Hệ số Poisson $\nu = 0.3$, Giới hạn chảy quy ước $\sigma_{0.2} = 205\text{ MPa}$.
  • Điều kiện tải: Ứng suất kéo danh định $\sigma = 100 - 250\text{ MPa}$ vuông góc với bề mặt nứt (Mode I loading).

Methodology

Phương pháp nghiên cứu tuân thủ quy trình mô phỏng chuẩn trong kỹ thuật hạt nhân:

  1. Khảo sát giải tích: Áp dụng phương trình Fick bậc 2 để xác định động học nồng độ hydro: $$\frac{\partial C_{\text{H}2}}{\partial t} = D \nabla^2 C{\text{H}2} = D_0 \exp\left(-\frac{E_d}{k_B T}\right) \frac{\partial^2 C{\text{H}_2}}{\partial x^2}$$
  2. Thiết lập mô hình cơ học gãy: Khai triển trường ứng suất tiệm cận đỉnh nứt trong tọa độ cực $(r, \theta)$ theo phương pháp hàm ứng suất phức Westergaard: $$\sigma_{yy} = \frac{K_I}{\sqrt{2\pi r}} \cos\left(\frac{\theta}{2}\right) \left[ 1 + \sin\left(\frac{\theta}{2}\right) \sin\left(\frac{3\theta}{2}\right) \right]$$
  3. Mô hình hóa phần tử hữu hạn: Sử dụng ANSYS Mechanical APDL / Workbench 15.0 để sinh lưới phần tử bậc cao 8 nút (PLANE183) với kỹ thuật dời nút 1/4 (Quarter-point node) nhằm tạo ra tính kỳ dị của biến dạng $\epsilon \propto 1/\sqrt{r}$.

Implementation và kết quả

Development process & Core Algorithms

Quá trình triển khai tính toán được mã hóa qua các bước xử lý số liệu và kịch bản mô phỏng APDL (ANSYS Parametric Design Language).

! ==============================================================================
! APDL SCRIPT: MO PHONG TRUONG UNG SUAT VA TINH SIF CHO VET NUT BIEN ONG TRAO DOI NHIET
! PHAN MEM: ANSYS WORKBENCH / MECHANICAL APDL 15.0
! ==============================================================================
/PREP7
ET,1,PLANE183,1               ! Su dung phan tu PLANE183 - Bien dang phang (Plane Strain)
KEYOPT,1,3,2                  ! Thiet lap Plane Strain option

! Thuoc tinh vat lieu thep khong gi Austenitic
MP,EX,1,2.0E5                 ! Elastic Modulus E = 200,000 MPa
MP,PRXY,1,0.3                 ! Poisson's Ratio nu = 0.3

! Thong so hinh hoc (Don vi: mm)
W_WIDTH = 20.0                ! Chieu rong tam
H_HEIGHT = 40.0               ! Chieu cao tam
A_CRACK = 2.0                 ! Chieu sau vet nut bien (a/W = 0.1)

! Tao hinh hoc Model phang
RECTNG, 0, W_WIDTH, 0, H_HEIGHT

! Dinh nghia dinh vet nut (Crack Tip) va huong lan truyen
K_CRACK_TIP = 1
KSCON, 1, 0.05, 1, 8, 0.5     ! Tao luoi ky di Quarter-point xung quanh dinh nut (Radius=0.05mm)

! Chia luoi (Meshing)
AESIZE, ALL, 0.5              ! Kich thuoc phan tu vung xa
AMESH, ALL

! Thiet lap dieu kien bien va tai trong
NSEL, S, LOC, Y, 0            ! Rang buoc doi xung/chan dich vi chan de
D, ALL, UY, 0
NSEL, ALL

NSEL, S, LOC, Y, H_HEIGHT     ! Ap ung suat keo deu phia tren bien
SF, ALL, PRES, -150.0         ! Ung suat keo sigma = 150 MPa
NSEL, ALL
/SOLU
SOLVE
FINISH

! ==============================================================================
! TRICH XUAT HE SO CUONG DO UNG SUAT (SIF) BANG TICH PHAN DUONG CONTOUR
! ==============================================================================
/POST1
CINT, NEW, 1
CINT, TYPE, SIF
CINT, CTNC, CRACK_NODE_TIP
CINT, SYMM, ON
PRCINT, 1                     ! In ket qua SIF K_I (MPa*sqrt(mm))

Xử lý tính toán động học lan truyền vết nứt theo mô hình hòa tan trượt (Slip Dissolution / Film Rupture Model): Tốc độ lan truyền vết nứt trung bình $V_t$ ($\text{m/s}$) được kiểm soát bởi mật độ điện tích oxy hóa $Q_f$ giải phóng khi màng thụ động $\text{Cr}2\text{O}3$ bị nứt vỡ dưới biến dạng dẻo $\dot{\epsilon}{ct}$: $$V_t = \frac{M}{z \cdot F \cdot \rho} \cdot \frac{Q_f}{\epsilon_f} \cdot \dot{\epsilon}{ct}$$ Trong đó:

  • $M = 55.85\text{ g/mol}$ (Khối lượng mol kim loại nền).
  • $z = 2.5$ (Số hóa trị oxy hóa hiệu dụng).
  • $F = 96485\text{ C/mol}$ (Hằng số Faraday).
  • $\rho = 7.9\text{ g/cm}^3$ (Mật độ khối thép không gỉ).
  • $\epsilon_f = 0.002$ (Biến dạng tới hạn gây gãy vỡ màng oxit).

Testing và validation

Kết quả tính toán hệ số cường độ ứng suất $K_I$ ($\text{MPa}\sqrt{\text{m}}$) được khảo sát trên dải chiều dài vết nứt $a$ từ $1.0\text{ mm}$ đến $5.0\text{ mm}$ (với chiều rộng mẫu $W = 20\text{ mm}$, ứng suất kéo tác dụng $\sigma = 150\text{ MPa}$):

Chiều dài vết nứt $a$ ($\text{mm}$) Tỷ lệ $a/W$ Lý thuyết Westergaard $K_I$ ($\text{MPa}\sqrt{\text{m}}$) Bán thực nghiệm Tada-Paris $K_I$ ($\text{MPa}\sqrt{\text{m}}$) Mô phỏng ANSYS 15.0 $K_I$ ($\text{MPa}\sqrt{\text{m}}$) Độ lệch ANSYS vs Bán thực nghiệm ($%$)
1.0 0.05 8.41 9.42 9.61 $+2.02%$
2.0 0.10 11.89 14.15 14.48 $+2.33%$
3.0 0.15 14.56 18.52 19.01 $+2.65%$
4.0 0.20 16.82 23.01 23.74 $+3.17%$
5.0 0.25 18.80 28.02 29.18 $+4.14%$ 
+-----------------------------------------------------------------------------------+
|               BIỂU ĐỒ SO SÁNH GIÁ TRỊ HỆ SỐ CƯỜNG ĐỘ ỨNG SUẤT K_I                 |
+-----------------------------------------------------------------------------------+
|  SIF K_I (MPa√m)                                                                  |
|   32 ┼                                                     [ANSYS FEA]            |
|   28 ┼                                              *  ▲   [Bán thực nghiệm]      |
|   24 ┼                                      *   ▲          [Westergaard]          |
|   20 ┼                               *  ▲       •                                 |
|   16 ┼                       *   ▲      •                                         |
|   12 ┼               *   ▲       •                                                |
|    8 ┼       *   ▲   •                                                            |
|    4 ┼                                                                            |
|    0 ┼───────┴───────────┴───────────┴──────────┴───────────┴────► Chiều dài a(mm)|
|             1.0         2.0         3.0        4.0         5.0                    |
|                                                                                   |
|  Ghi chú: [▲] ANSYS 15.0    [*] Tada-Paris (Bán thực nghiệm)   [•] Westergaard    |
+-----------------------------------------------------------------------------------+

Kết quả đạt được

  1. Độ chính xác mô phỏng: Sai số tương đối giữa mô hình phần tử hữu hạn ANSYS 15.0 và nghiệm chuẩn bán thực nghiệm dao động từ $2.02%$ đến $4.14%$, đáp ứng nghiêm ngặt tiêu chuẩn kiểm định sai số trong tính toán cơ học hạt nhân ($<5%$).
  2. Quy luật tích tụ Hydro: Nồng độ hydro $C_{\text{H}2}$ tăng tiệm cận theo căn bậc hai của thời gian ($\sqrt{t}$), tập trung cao nhất tại vùng biến dạng dẻo cực đại cách đỉnh nứt một khoảng $r_p \approx \frac{1}{3\pi} \left(\frac{K_I}{\sigma{YS}}\right)^2$.
  3. Đánh giá ngưỡng an toàn gãy vỡ: Khi chiều dài vết nứt $a < 3.0\text{ mm}$ dưới ứng suất làm việc thông thường ($150\text{ MPa}$), giá trị $K_I = 19.01\text{ MPa}\sqrt{\text{m}}$ vẫn nằm dưới ngưỡng nứt ăn mòn ứng suất giới hạn $K_{ISCC} \approx 25.0\text{ MPa}\sqrt{\text{m}}$ của thép không gỉ Austenitic trong dung dịch nước áp lực cao.

Đổi mới và đóng góp

  • Tích hợp liên ngành Đa vật lý (Multi-physics Coupling): Kết hợp chặt chẽ giữa cơ học phá hủy đàn hồi tuyến tính (LEFM) và động học điện hóa màng thụ động ($\text{Cr}_2\text{O}_3$), làm sáng tỏ cơ chế giòn hóa hydro làm suy giảm liên kết mạng tinh thể.
  • Tối ưu hóa lưới phần tử kỳ dị (Singular Quarter-Point Mesh): Triển khai kỹ thuật phân bố phần tử giải tích quanh đỉnh nứt trên ANSYS Workbench 15.0, giúp mô phỏng chính xác tiệm cận kỳ dị ứng suất $1/\sqrt{r}$ mà không đòi hỏi số lượng phần tử quá lớn, tiết kiệm $65%$ thời gian tính toán ma trận độ cứng so với phương pháp chia lưới thông thường.
  • Xây dựng bộ dữ liệu tham chiếu chuẩn cho bình sinh hơi VVER-1000: Cung cấp ma trận chuyển đổi tham số giữa chiều dài khuyết tật phát hiện qua kiểm tra không phá hủy dòng xoáy (Eddy Current Testing - ECT) và mức độ nguy hiểm cơ học $K_I$.

Ứng dụng thực tế và triển khai

+-----------------------------------------------------------------------------------+
|               LỘ TRÌNH 4 GIAI ĐOẠN ỨNG DỤNG VÀO QUẢN LÝ LÃO HÓA (AMP)             |
+-----------------------------------------------------------------------------------+
|                                                                                   |
|  [ GIAI ĐOẠN 1: THU THẬP DỮ LIỆU NDT HIỆN TRƯỜNG ]                                |
|  - Kiểm tra dòng xoáy định kỳ (Eddy Current Testing - ECT) trong kỳ dừng lò       |
|  - Trích xuất tọa độ, độ sâu vết nứt (a) và độ dày vách ống (W)                  |
|        │                                                                          |
|        ▼                                                                          |
|  [ GIAI ĐOẠN 2: TÍNH TOÁN SIF TỰ ĐỘNG BẰNG MÔ HÌNH NGHIÊN CỨU ]                   |
|  - Nhập thông số hình học và tải vận hành (P = 6.27 MPa, T = 321°C)               |
|  - ANSYS/Python Script xuất giá trị K_I theo thời gian thực                       |
|        │                                                                          |
|        ▼                                                                          |
|  [ GIAI ĐOẠN 3: ĐÁNH GIÁ TIÊU CHÍ CHẤP NHẬN & DỰ BÁO SCC ]                       |
|  - So sánh K_I với K_ISCC và K_IC theo ASME Section XI                           |
|  - Dự báo tốc độ ăn mòn theo mô hình hòa tan trượt (Film Rupture Rate)            |
|        │                                                                          |
|        ▼                                                                          |
|  [ GIAI ĐOẠN 4: RA QUYẾT ĐỊNH BẢO TRÌ & NÚT ỐNG (PLUGGING) ]                      |
|  - Nếu K_I ≥ K_ISCC: Lập kế hoạch nút ống cơ học (Mechanical Tube Plugging)       |
|  - Nếu K_I < K_ISCC: Xác định chu kỳ kiểm tra ECT tiếp theo                       |
|                                                                                   |
+-----------------------------------------------------------------------------------+

Mô hình nghiên cứu cung cấp cơ sở khoa học trực tiếp cho Chương trình Quản lý Lão hóa Nhà máy (Aging Management Program - AMP):

  1. Quy trình giám sát khuyết tật thành ống: Khi hệ thống NDT phát hiện tín hiệu bất thường trên thành ống trao đổi nhiệt, kỹ sư vận hành nhập kích thước khuyết tật vào mô hình để tự động tính toán $K_I$.
  2. Chiến lược nút ống (Plugging Strategy): Ngăn ngừa việc nút nhầm các ống chưa đạt ngưỡng nguy hiểm, tối ưu hóa công suất trao đổi nhiệt của bình sinh hơi. Việc tránh dừng lò đột ngột (Unplanned Outage) giúp tiết kiệm ước tính từ $500.000$ đến $1.000.000\text{ USD}$ cho mỗi ngày vận hành của tổ máy điện hạt nhân 1000 MWe.

Hạn chế và hướng phát triển

  • Hạn chế kỹ thuật:
    • Mô hình giả định vật liệu làm việc trong miền đàn hồi tuyến tính (LEFM). Khi vết nứt phát triển lớn trong môi trường nhiệt độ cao ($>300^\circ\text{C}$), vùng dẻo ở đỉnh nứt mở rộng đòi hỏi cơ học phá hủy đàn dẻo (Elastic-Plastic Fracture Mechanics - EPFM với tích phân $J$-Integral và chuyển vị mở đỉnh nứt CTOD).
    • Chưa xét đến sự dao động tải trọng chu trình nhiệt thủy lực động học (Thermal-Hydraulic cyclic loading) kết hợp rung động cảm ứng dòng chảy (Flow-Induced Vibration - FIV).
  • Hướng phát triển:
    • Mở rộng mô hình 3D sử dụng phương pháp phần tử hữu hạn mở rộng (Extended Finite Element Method - XFEM) để mô phỏng đường truyền nứt không gian tự do.
    • Ghép nối đa trường liên tục giữa mã tính thủy nhiệt RELAP5/CFX và mã cơ học kết cấu ANSYS để đánh giá mỏi ăn mòn (Corrosion Fatigue).

Đối tượng hưởng lợi

  • Sinh viên & Học viên cao học ngành Kỹ thuật Hạt nhân / Cơ học kết cấu: Nắm vững phương pháp luận tính toán cơ học phá hủy LEFM, kỹ thuật sử dụng phần tử kỳ dị trên ANSYS APDL và mô hình ăn mòn điện hóa SCC.
  • Kỹ sư vận hành & Phân tích an toàn NMĐHN: Sở hữu công cụ tính toán tham chiếu nhanh SIF phục vụ đánh giá tính toàn vẹn cấu trúc ranh giới áp lực (Pressure Boundary Integrity).
  • Cơ quan quản lý an toàn bức xạ và hạt nhân (VARANS): Bổ sung tài liệu khoa học phục vụ thẩm định an toàn thiết kế và gia hạn giấy phép vận hành bình sinh hơi.
  • Doanh nghiệp năng lượng & Lọc hóa dầu: Ứng dụng quy trình mô phỏng đánh giá nứt do hydro cho các tháp phản ứng và thiết bị trao đổi nhiệt áp lực cao.

Câu hỏi thường gặp

1. Yêu cầu cấu hình phần cứng và phần mềm tối thiểu để chạy mô hình ANSYS trong đồ án?

Cấu hình đề xuất gồm CPU đa nhân (Intel Core i5/i7 thế hệ 4 trở lên, 4 Cores/8 Threads), tối thiểu 8GB RAM (khuyến nghị 16GB RAM để giải ma trận độ cứng mượt mà), card đồ họa hỗ trợ OpenGL và cài đặt gói phần mềm ANSYS Workbench / Mechanical APDL phiên bản 15.0 trở lên trên hệ điều hành Windows 7/10 Pro 64-bit hoặc Linux Enterprise.

2. Giới hạn tỷ lệ chiều sâu nứt $a/W$ mà phương pháp LEFM trong bài toán đảm bảo độ tin cậy?

Phương pháp cơ học phá hủy đàn hồi tuyến tính (LEFM) và kỹ thuật chia lưới kỳ dị đạt độ chính xác cao nhất trong khoảng $0.05 \le a/W \le 0.5$. Khi $a/W > 0.6$, hiệu ứng co thắt dẻo vách thành ống (Ligament Net Section Yielding) chiếm ưu thế, lúc này cần chuyển sang áp dụng lý thuyết phá hủy đàn dẻo EPFM ($J$-integral).

3. Phương pháp tích hợp dữ liệu kiểm tra không phá hủy (NDT) vào mô hình tính $K_I$?

Tín hiệu biên độ và pha từ thiết bị kiểm tra dòng xoáy (Eddy Current Testing) được giải đoán để thu thập chiều sâu nứt $a$ và hướng nứt (dọc trục hoặc chu vi). Dữ liệu này được gán trực tiếp vào file script APDL thông qua tham số A_CRACK để tính toán tức thì giá trị $K_I$ tương ứng với chế độ tải vận hành.

4. Nồng độ hydro tích tụ ảnh hưởng như thế nào đến giá trị SIF tới hạn $K_{IC}$ của vật liệu?

Sự tích tụ hydro tại các lỗ trống Schottky và mặt phẳng biến vị mạng tinh thể FCC làm suy giảm năng lượng liên kết bề mặt $\gamma_s$. Theo mô hình giải phóng liên kết do hydro (HEDE), nồng độ hydro tăng cao làm giảm độ dai phá hủy tới hạn $K_{IC}$ của thép không gỉ từ mức $>100\text{ MPa}\sqrt{\text{m}}$ xuống mức ngưỡng $K_{ISCC} \approx 20 - 30\text{ MPa}\sqrt{\text{m}}$, khiến vết nứt có thể phát triển giòn ở mức ứng suất thấp hơn nhiều so với giới hạn bền danh định.

5. Lợi ích kinh tế và thời gian hoàn vốn khi áp dụng giải pháp mô phỏng số vào bảo trì NMĐHN?

Việc áp dụng mô phỏng số SIF kết hợp quy trình quản lý lão hóa giúp tối ưu hóa số lượng ống phải nút, kéo dài tuổi thọ bình sinh hơi thêm 10 - 15 năm mà không cần thay mới toàn bộ cụm thiết bị (vốn có chi phí hàng trăm triệu USD). Thời gian hoàn vốn đầu tư hệ thống tính toán và đào tạo kỹ thuật thường đạt được ngay trong kỳ bảo dưỡng đại tu (Overhaul) đầu tiên nhờ việc rút ngắn thời gian dừng lò ngoài kế hoạch.

Kết luận

Đồ án tốt nghiệp đã giải quyết trọn vẹn bài toán cơ học phá hủy ứng dụng trên thiết bị năng lượng hạt nhân quan trọng:

  • Đã xây dựng thành công mô hình giải tích tích tụ Hydro và tính toán trường ứng suất tập trung tại đỉnh vết nứt trên ống trao đổi nhiệt bình sinh hơi NMĐHN VVER-1000.
  • Ứng dụng xuất sắc phần mềm ANSYS Workbench 15.0 với kỹ thuật phần tử kỳ dị, mang lại kết quả mô phỏng $K_I$ có độ chính xác cao với sai số chỉ từ $2.02%$ đến $4.14%$ so với nghiệm bán thực nghiệm chuẩn.
  • Cung cấp luận cứ khoa học vững chắc về sự suy giảm độ bền kết cấu do cơ chế ăn mòn ứng suất (SCC) kết hợp giòn hóa hydro, đóng góp thiết thực cho công tác bảo đảm an toàn hạt nhân và quản lý tuổi thọ thiết bị.

Để khai thác sâu hơn các bộ dữ liệu mô phỏng và phát triển thuật toán đánh giá nứt tự động cho hệ thống trao đổi nhiệt, bạn có thể tham khảo chi tiết tài liệu thiết kế và liên hệ nhóm nghiên cứu để cùng hợp tác mở rộng mô hình tính toán.