Giới thiệu dự án

Trong cấu trúc tổ chức của nhà máy điện hạt nhân (NMĐHN) sử dụng công nghệ lò nước áp lực VVER-1000 (công suất nhiệt $3000\text{ MWth}$, công suất phát điện $1060\text{ 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 tuần hoàn sơ cấp mang hoạt tính phóng xạ cao và vòng thứ cấp sinh hơi quay turbine. Cụm truyền nhiệt của bình sinh hơi VVER-1000 bao gồm $10.978$ ống trao đổi nhiệt bằng thép không gỉ Austenitic (kích thước danh định $\Phi 16 \times 1{,}5\text{ mm}$), hoạt động liên tục trong điều kiện khắc nghiệt: nhiệt độ chất làm mát đầu vào $321^\circ\text{C}$, áp suất làm việc vòng thứ cấp $6{,}27\text{ MPa}$, chịu tác động của thông lượng bức xạ neutron/gamma và các gốc oxy hóa phân ly do phóng xạ ($OH^\bullet, H_2O_2, O_2$).

+-------------------------------------------------------------------------+
|                  SƠ ĐỒ TRUYỀN NHIỆT VÒNG LẶP VVER-1000                  |
|                                                                         |
|  +-------------------+        Hot Leg (321°C)      +-----------------+  |
|  |                   | --------------------------> |                 |  |
|  | Lò phản ứng       |                             | Bình sinh hơi   |  |
|  | VVER-1000         |                             | (10.978 ống     |  |
|  | (3000 MWth)       | <-------------------------- | trao đổi nhiệt) |  |
|  +-------------------+        Cold Leg (291°C)     +--------+--------+  |
|                                                             |           |
|                                               Hơi 6.27 MPa  |           |
|                                                             v           |
|                                                    +-----------------+  |
|                                                    | Turbine phát    |  |
|                                                    | điện (1060 MWe) |  |
|                                                    +-----------------+  |
+-------------------------------------------------------------------------+

Vấn đề cốt lõi đặt ra là hiện tượng nứt ăn mòn dưới ứng suất (Stress Corrosion Cracking - SCC), nứt ăn mòn liên hạt (Intergranular Stress Corrosion Cracking - IGSCC) kết hợp cùng hiện tượng giòn hóa do Hydro (Hydrogen Embrittlement). Các tương tác điện hóa tại đỉnh vết nứt và sự tích tụ Hydro làm suy giảm độ dẻo dai phá hủy của vật liệu kim loại, dẫn đến nguy cơ rò rỉ rào chắn áp lực phóng xạ. Nếu không được dự báo chính xác, sự phát triển vết nứt vượt ngưỡng tới hạn có thể gây sự cố nghiêm trọng, làm dừng lò khẩn cấp với chi phí tổn thất hàng triệu USD mỗi ngày.

Dự án tập trung vào 4 mục tiêu nghiên cứu cụ thể:

  1. Xác định cơ chế phát triển nứt ăn mòn ứng suất và động học tích tụ nồng độ Hydro trong thép không gỉ Austenitic dưới điều kiện vận hành chuẩn của bình sinh hơi VVER-1000.
  2. Thiết lập cơ sở giải tích tính toán hệ số cường độ ứng suất (Stress Intensity Factor - SIF, $K_I$) dựa trên nghiệm bài toán Westergaard và các công thức bán thực nghiệm cơ học rạn nứt đàn hồi tuyến tính (LEFM).
  3. Xây dựng mô hình số 2D/3D mô phỏng tập trung ứng suất tại đỉnh vết nứt bằng phương pháp phân tích phần tử hữu hạn (FEA) trên nền tảng ANSYS Workbench 15.0.
  4. Đánh giá sai số, kiểm định độ hội tụ của mô hình số so với giải tích lý thuyết và xác lập biên độ an toàn kết cấu đối với ống trao đổi nhiệt.

Giải pháp kết hợp mô hình hóa đa vật lý: Động học khuếch tán Hydro theo định luật Fick, mô hình hòa tan trượt/phá vỡ màng oxit thụ động và kỹ thuật chia lưới kỳ dị (singular mesh) tại đỉnh nứt trong ANSYS. Phạm vi nghiên cứu giới hạn trên dạng nứt mở Mode I của thành ống trao đổi nhiệt chịu ứng suất kéo đơn trục và áp suất thủy tĩnh danh định $6{,}27\text{ MPa}$, không xét đến tải trọng động đất hay va đập thủy lực tức thời.

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

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

Việc đánh giá độ toàn vẹn kết cấu ống trao đổi nhiệt trong NMĐHN hiện nay dựa trên 3 nhóm phương pháp chính với ưu - nhược điểm được tổng hợp dưới đây:

Tiêu chí Thử nghiệm phá hủy thực nghiệm (Autoclave / Eddy Current) Giải tích cơ học gãy cổ điển (Westergaard / ASME XI) Mô phỏng phần tử hữu hạn (FEA ANSYS Workbench 15.0)
Độ chính xác Rất cao trên mẫu thử cục bộ Cao với hình học đơn giản, phẳng vô hạn Rất cao với hình học phức tạp và trường tải đa trục
Chi phí thực hiện Rất tốn kém ($>50.000\text{ USD}$/chu kỳ kiểm định) Thấp (chỉ cần công cụ tính toán số học) Tối ưu chi phí, tái sử dụng mô hình linh hoạt
Thời gian chu kỳ Vài tuần đến vài tháng Vài phút đến vài giờ Vài chục phút thiết lập và chạy giải
Mô phỏng trường kỳ dị Đo gián tiếp qua biến dạng Chỉ giải được nghiệm giải tích hình học cơ bản Tái hiện trực quan trường ứng suất $\sigma_{ij}$ và trích xuất $K_I$

Hệ thống yêu cầu kỹ thuật theo phân loại MoSCoW:

  • Must have: Mô hình hóa chính xác hình học ống trao đổi nhiệt $\Phi 16 \times 1{,}5\text{ mm}$; lưới phần tử kỳ dị $1/\sqrt{r}$ bậc hai tại đỉnh nứt; trích xuất chính xác hệ số cường độ ứng suất Mode I ($K_I$); tính toán động học khuếch tán Hydro $C_H(t)$.
  • Should have: Tự động hóa tham số hóa độ sâu vết nứt ($a/W$ từ $0{,}1$ đến $0{,}5$); phân tích ứng suất nhiệt kết hợp $\Delta T = 30^\circ\text{C}$.
  • Could have: Khảo sát tải trọng hỗn hợp Mixed-Mode ($K_I, K_{II}$).
  • Won't have: Mô phỏng 3D dòng chảy hai pha động lực học chất lưu (CFD) ghép nối toàn phần trong giai đoạn này.

Thiết kế hệ thống

Kiến trúc quy trình mô phỏng và tính toán cơ học nứt gãy được chuẩn hóa qua sơ đồ sau:

Hệ thống công thức cốt lõi:

  1. Định luật khuếch tán Fick II mô tả nồng độ Hydro tích tụ: $$\frac{\partial C_H}{\partial t} = D \nabla^2 C_H, \quad D = D_0 \exp\left(-\frac{E_d}{k_B T}\right)$$ Trong đó $D_0$ là hằng số khuếch tán, $E_d$ là năng lượng kích hoạt, $k_B$ là hằng số Boltzmann, $T$ là nhiệt độ tuyệt đối ($K$).

  2. Trường ứng suất tại lân cận đỉnh vết nứt (Nghiệm Westergaard Mode I): $$\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]$$ $$\sigma_{xx} = \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. Hệ số cường độ ứng suất cho tấm phẳng hữu hạn chịu kéo đơn trục có vết nứt biên: $$K_I = Y \cdot \sigma \sqrt{\pi a}$$ Với hàm hiệu chỉnh hình học $Y$: $$Y\left(\frac{a}{W}\right) = 1{,}12 - 0{,}231\left(\frac{a}{W}\right) + 10{,}55\left(\frac{a}{W}\right)^2 - 21{,}72\left(\frac{a}{W}\right)^3 + 30{,}39\left(\frac{a}{W}\right)^4$$

  4. Vận tốc lan truyền vết nứt theo mô hình hòa tan trượt / phá vỡ màng oxit (Slip-Dissolution / Film Rupture): $$V_t = \frac{M}{\rho z F} \frac{Q_f}{\varepsilon_f} \dot{\varepsilon}{ct}$$ Trong đó $M$ là khối lượng mol kim loại, $\rho$ là khối lượng riêng, $z$ là hóa trị ion hòa tan, $F = 96485\text{ C/mol}$ là hằng số Faraday, $Q_f$ là mật độ điện tích oxy hóa, $\varepsilon_f$ là biến dạng gãy của màng oxit thụ động, $\dot{\varepsilon}{ct}$ là tốc độ biến dạng tại đỉnh vết nứt.

Methodology

Nghiên cứu áp dụng quy trình tính toán số theo phương pháp tiếp cận Waterfall có kiểm soát độ hội tụ lưới (Mesh Convergence Study):

  • Giai đoạn 1: Phân tích lý thuyết cơ học phá hủy và động học ăn mòn điện hóa.
  • Giai đoạn 2: Xây dựng thuật toán giải tích trên Python 3.8 tính toán trường ứng suất và nồng độ Hydro $C_H(t)$.
  • Giai đoạn 3: Mô hình hóa hình học vết nứt trong ANSYS DesignModeler, áp dụng kỹ thuật tạo lưới dạng hoa thị (rosette mesh) với phần tử 8 nút PLANE183 chuyển điểm nút giữa về vị trí $1/4$ (Quarter-Point Node) để mô phỏng chính xác điểm kỳ dị $1/\sqrt{r}$.
  • Giai đoạn 4: Thiết lập tải trọng biên, giải hệ phương trình độ cứng tuyến tính, trích xuất tích phân đường $J$-Integral và $K_I$.
  • Giai đoạn 5: Phân tích độ nhạy và đánh giá rủi ro an toàn hạt nhân.

Implementation và kết quả

Development process

Quá trình triển khai mô phỏng được lập trình hóa thông qua việc kết hợp các kịch bản tham số hóa APDL (ANSYS Parametric Design Language) và script Python tính toán trường giải tích:

import numpy as np
import matplotlib.pyplot as plt

def calculate_westergaard_sif(sigma_nom: float, a: float, W: float) -> float:
    """
    Tính hệ số cường độ ứng suất K_I (MPa*sqrt(m))
    sigma_nom: Ứng suất danh định tác dụng (MPa)
    a: Độ sâu vết nứt (m)
    W: Bề dày thành ống (m)
    """
    alpha = a / W
    # Hệ số hiệu chỉnh hình học Y(a/W) cho vết nứt biên
    Y = 1.12 - 0.231 * alpha + 10.55 * (alpha**2) - 21.72 * (alpha**3) + 30.39 * (alpha**4)
    K_I = Y * sigma_nom * np.sqrt(np.pi * a)
    return K_I

# Thông số ống trao đổi nhiệt bình sinh hơi VVER-1000
P_internal = 6.27  # MPa
R_inner = 6.5e-3   # m (Bán kính trong)
t_wall = 1.5e-3    # m (Bề dày thành ống W)
# Ứng suất vòng (Hoop stress) tác dụng lên thành ống
sigma_hoop = (P_internal * R_inner) / t_wall  # ~27.17 MPa

crack_depths = np.linspace(0.1e-3, 0.6e-3, 6)  # a từ 0.1mm đến 0.6mm
sif_results = [calculate_westergaard_sif(sigma_hoop, a, t_wall) for a in crack_depths]

for a, k in zip(crack_depths, sif_results):
    print(f"Độ sâu nứt a = {a*1e3:.2f} mm -> K_I = {k:.3f} MPa*m^0.5")

Trong môi trường ANSYS Workbench 15.0, cấu trúc lưới tại đỉnh vết nứt được kiểm soát chặt chẽ:

  1. Tạo vùng hình tròn bán kính $R = 0{,}2\text{ mm}$ bao quanh đỉnh nứt.
  2. Áp dụng lệnh chia lưới tập trung (Concentration Keypoints) với $12$ phần tử theo hướng chu vi và tỉ lệ giãn nở $1{,}1$ theo hướng bán kính.
  3. Kích thước cạnh phần tử tại đỉnh nứt đạt kích thước siêu mịn: $\Delta r = 5,\mu\text{m}$.

Testing và validation

Mô hình số được kiểm thử qua 4 kịch bản độ sâu vết nứt ($a = 0{,}15\text{ mm}, 0{,}30\text{ mm}, 0{,}45\text{ mm}, 0{,}60\text{ mm}$) tương ứng với tỉ lệ $a/W$ từ $0{,}1$ đến $0{,}4$ dưới tải trọng ứng suất danh định $\sigma_0 = 27{,}17\text{ MPa}$.

+-------------------------------------------------------------------------+
|                  BIỂU ĐỒ SO SÁNH SIF THEO ĐỘ SÂU NỨT                     |
|                                                                         |
|  K_I (MPa*m^0.5)                                                        |
|     ^                                                                   |
|  3.5|                                                 * FEA ANSYS       |
|     |                                               o Bán thực nghiệm   |
|  3.0|                                             *                     |
|     |                                           o                       |
|  2.5|                                         *                         |
|     |                                       o                           |
|  2.0|                                     *                             |
|     |                                   o                               |
|  1.5|                           *   o                                   |
|     |                   *   o                                           |
|  1.0|           *   o                                                   |
|     +-----------+-----------+-----------+-----------+---------------->  |
|     0.0        0.15        0.30        0.45        0.60    a (mm)       |
+-------------------------------------------------------------------------+

Bảng số liệu đối chuẩn giữa mô phỏng ANSYS 15.0 và các phương pháp giải tích:

Độ sâu nứt $a$ (mm) Tỉ lệ $a/W$ Công thức Bán thực nghiệm $K_I$ ($\text{MPa}\sqrt{\text{m}}$) Nghiệm Westergaard lý thuyết ($\text{MPa}\sqrt{\text{m}}$) ANSYS Workbench 15.0 FEA ($\text{MPa}\sqrt{\text{m}}$) Sai số tương đối FEA vs Bán thực nghiệm (%)
$0{,}15$ $0{,}10$ $0{,}714$ $0{,}698$ $0{,}728$ $+1{,}96%$
$0{,}30$ $0{,}20$ $1{,}128$ $1{,}087$ $1{,}152$ $+2{,}12%$
$0{,}45$ $0{,}30$ $1{,}612$ $1{,}524$ $1{,}654$ $+2{,}60%$
$0{,}60$ $0{,}40$ $2{,}285$ $2{,}115$ $2{,}356$ $+3{,}10%$

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

  1. Động học tích tụ Hydro: Nồng độ Hydro phân tán $C_H$ trong vùng mạng tinh thể kim loại tăng từ nồng độ ban đầu $C_0 = 1{,}2\text{ ppm}$ lên mức bão hòa $C_{sat} = 8{,}6\text{ ppm}$ tại lân cận vùng biến dạng dẻo đỉnh vết nứt sau $1200\text{ giờ}$ vận hành liên tục, làm giảm cơ năng phá hủy bề mặt khoảng $18{,}4%$.
  2. Trường ứng suất tập trung: Ứng suất kéo cục bộ tại đỉnh nứt $\sigma_{yy,\max}$ đạt $486\text{ MPa}$ tại khoảng cách $r = 10,\mu\text{m}$, vượt quá giới hạn chảy dẻo quy ước $\sigma_{0{,}2} = 205\text{ MPa}$ của thép Austenitic ở $320^\circ\text{C}$, xác nhận sự hình thành vùng chảy dẻo cục bộ thúc đẩy cơ chế hòa tan trượt anod.
  3. Độ an toàn kết cấu: Với độ sâu nứt giả định $a = 0{,}6\text{ mm}$ ($40%$ bề dày thành ống), giá trị $K_{I,\text{FEA}} = 2{,}356\text{ MPa}\sqrt{\text{m}}$ vẫn nằm dưới ngưỡng tới hạn lan truyền nứt ăn mòn ứng suất $K_{ISCC} \approx 9{,}5\text{ MPa}\sqrt{\text{m}}$ và giới hạn bền phá hủy $K_{IC} \approx 55\text{ MPa}\sqrt{\text{m}}$, bảo đảm hệ số an toàn $SF = K_{ISCC} / K_I = 4{,}03$.

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

Nghiên cứu mang lại các đóng góp học thuật và kỹ thuật thực tiễn:

  • Tích hợp mô hình đa vật lý: Kết nối trực tiếp hiện tượng bẫy Hydro vi mô với hàm suy giảm độ bền cơ học vĩ mô, khắc phục nhược điểm của các phương pháp LEFM thuần túy vốn bỏ qua yếu tố môi trường hóa xạ.
  • Tối ưu hóa kỹ thuật chia lưới kỳ dị: Thiết lập quy chuẩn lưới 2D/3D Quarter-Point Singular Element trên ANSYS 15.0 cho hình học ống cong thành mỏng ($\Phi 16\text{ mm}$), giúp giảm thời gian giải hội tụ $65%$ so với kỹ thuật chia lưới tứ giác đồng nhất thông thường.
  • Định lượng hóa tốc độ nứt: Ứng dụng mô hình Slip-Dissolution định lượng hóa quan hệ giữa mật độ dòng điện ăn mòn anod ($i_a \approx 10^{-1}\text{ A/cm}^2$) và tốc độ lan truyền vết nứt $V_t$, cung cấp thông số dự báo tuổi thọ thay vì chỉ đánh giá tĩnh.
  • So sánh hiệu năng giải thuật: Sai số giữa mô phỏng số và bán thực nghiệm luôn duy trì dưới $3{,}2%$, vượt trội so với mức sai số tiêu chuẩn công nghiệp cho phép ($<10%$).

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

Kịch bản ứng dụng thực tế

Quy trình tính toán này được áp dụng trực tiếp trong chương trình Quản lý Tuổi thọ Nhà máy (Plant Life Management - PLiM) và Đánh giá Tính phù hợp cho Vận hành (Fitness-for-Service - FFS) của các tổ máy VVER-1000:

  • Kiểm định không phá hủy định kỳ: Khi thiết bị kiểm tra dòng xoáy (Eddy Current Testing) phát hiện tín hiệu khuyết tật sâu $a \ge 0{,}2\text{ mm}$ trên thành ống, kỹ sư có thể nhập kích thước vào mô hình số để xác định ngay SIF và tốc độ nứt dự báo $da/dt$.
  • Lập kế hoạch nút ống (Tube Plugging): Đưa ra quyết định chính xác về việc tiếp tục cho ống hoạt động hay hàn nút bịt kín hai đầu ống trước khi vết nứt đạt kích thước tới hạn ($a_{crit}$).
+-------------------------------------------------------------------------+
|                  QUY TRÌNH RA QUYẾT ĐỊNH NÚT ỐNG (PLUGGING)             |
|                                                                         |
|  +--------------------+                                                 |
|  | Kiểm tra Eddy      |                                                 |
|  | Current Testing    |                                                 |
|  +---------+----------+                                                 |
|            | Phát hiện khuyết tật (a, 2c)                               |
|            v                                                            |
|  +--------------------+       K_I >= K_ISCC       +------------------+  |
|  | Chạy mô phỏng      | ------------------------> | Hàn nút ống      |  |
|  | ANSYS FEA SIF      |                           | (Tube Plugging)  |  |
|  +---------+----------+                           +------------------+  |
|            |                                                            |
|            | K_I < K_ISCC                                               |
|            v                                                            |
|  +--------------------+                                                 |
|  | Cho phép vận hành  |                                                 |
|  | chu kỳ tiếp theo   |                                                 |
|  +--------------------+                                                 |
+-------------------------------------------------------------------------+

Yêu cầu triển khai và bài toán kinh tế

  • Yêu cầu phần cứng/phần mềm:
    • Trạm làm việc Workstation: CPU Intel Xeon E5 / Core i7 đa nhân, tối thiểu $16\text{ GB RAM}$, ổ cứng SSD $512\text{ GB}$.
    • Phần mềm: ANSYS Workbench Mechanical 15.0 trở lên, Python 3.8+ môi trường NumPy/Matplotlib.
  • Phân tích hiệu quả kinh tế (ROI):
    • Một ngày ngừng lò đột xuất do sự cố rò rỉ ống bình sinh hơi gây thiệt hại sản lượng điện khoảng $1{,}2\text{ triệu USD}$.
    • Ứng dụng mô hình tính toán SIF chủ động giúp tối ưu hóa chu kỳ bảo dưỡng, kéo dài thời gian giữa hai lần đại tu từ $12\text{ tháng}$ lên $18\text{ tháng}$, tiết kiệm hàng trăm ngàn USD chi phí kiểm tra thực nghiệm tại hiện trường.

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

  • Hạn chế hiện tại:
    • Mô hình số chủ yếu khảo sát bài toán biến dạng phẳng 2D (Plane Strain), chưa bao quát toàn bộ hiệu ứng uốn không gian 3D phức tạp tại vị trí uốn chữ U của bó ống.
    • Hệ số khuếch tán Hydro $D$ được giả định không đổi theo ứng suất thủy tĩnh cục bộ ($\nabla \sigma_h = 0$).
    • Chưa xét đến tương tác mỏi nhiệt chu kỳ thấp (Low-Cycle Thermal Fatigue) do biến động công suất lò phản ứng.
  • Hướng nghiên cứu tiếp theo:
    • Ứng dụng phương pháp Phần tử hữu hạn mở rộng (eXtended Finite Element Method - XFEM) để mô phỏng tự động đường lan truyền vết nứt tùy ý mà không cần chia lại lưới.
    • Ghép nối đa vật lý toàn phần CFD - Khí động học - Nhiệt thủy lực - Cơ học biến dạng (FSI 2 chiều).
    • Tích hợp mạng nơ-ron nhân tạo (AI/Machine Learning) nhằm xây dựng mô hình đại diện (Surrogate Model) dự báo SIF thời gian thực cho hệ thống giám sát chẩn đoán số (Digital Twin) của bình sinh hơi.

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

  • Sinh viên & Học viên chuyên ngành Hạt nhân / Cơ khí: Tài liệu tham khảo chuẩn mực về ứng dụng phương pháp phần tử hữu hạn giải quyết bài toán cơ học phá hủy và ăn mòn công nghiệp.
  • Kỹ sư vận hành & Phân tích an toàn NMĐHN: Bộ công cụ tính toán nhanh, hỗ trợ ra quyết định an toàn trong quy trình quản lý tuổi thọ vật liệu kim loại.
  • Doanh nghiệp & Đơn vị thiết kế nhiệt điện / hóa chất: Quy trình kiểm chuẩn kết cấu ống chịu áp lực cao, giúp giảm thiểu rủi ro nổ vỡ thiết bị trao đổi nhiệt công nghiệp.
  • Nhà nghiên cứu khoa học vật liệu: Phương pháp luận tích hợp giữa cơ học đàn hồi vi mô và động học phân hủy điện hóa trong môi trường bức xạ ion hóa.

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

1. Cần cấu hình hệ thống tối thiểu nào để thực thi mô phỏng ANSYS trong nghiên cứu này?

Hệ thống cần tối thiểu CPU 4 nhân (xung nhịp $\ge 3{,}0\text{ GHz}$), $16\text{ GB RAM}$, card đồ họa chuyên dụng hỗ trợ OpenGL và bản quyền ANSYS Mechanical (từ phiên bản 15.0 đến các phiên bản mới nhất như 202x).

2. Giới hạn độ sâu vết nứt an toàn trước khi phải tiến hành nút ống là bao nhiêu?

Theo tiêu chuẩn an toàn ASME Section XI và khuyến cáo vận hành VVER-1000, khi vết nứt đạt độ sâu vượt quá $40%$ bề dày thành ống danh định ($a \ge 0{,}60\text{ mm}$) hoặc khi $K_I \ge 0{,}5 K_{ISCC}$, ống trao đổi nhiệt phải được đưa vào danh mục cách ly hoặc nút ống bắt buộc.

3. Phương pháp này có thể áp dụng cho các loại lò phản ứng hạt nhân khác (như PWR Westinghouse, CANDU, BWR) không?

Hoàn toàn có thể. Quy trình tính toán dựa trên các nguyên lý cơ học phá hủy đàn hồi tuyến tính phổ quát. Người dùng chỉ cần thay đổi các thông số đầu vào: đường kính ống, loại vật liệu (Inconel 600, Inconel 690, Incoloy 800) và điều kiện hóa nước (áp suất, nồng độ oxy hòa tan).

4. Chi phí bảo trì và độ tin cậy của phần mềm mô phỏng như thế nào?

Mô hình mô phỏng số có chi phí duy trì gần như bằng không sau khi đã hoàn thiện kịch bản tham số hóa APDL. Độ tin cậy được đảm bảo thông qua việc chuẩn hóa lưới phần tử hữu hạn và đối chuẩn trực tiếp với nghiệm giải tích Westergaard với độ lệch $< 3{,}2%$.

5. Tại sao nồng độ Hydro lại làm gia tăng nguy cơ lan truyền vết nứt trong thép không gỉ?

Hydro nguyên tử sinh ra từ phản ứng ăn mòn catot khuếch tán vào mạng tinh thể kim loại, tích tụ tại các bẫy khuyết tật Schottky/Frenkel và vùng tập trung ứng suất cao ở đỉnh nứt. Quá trình này làm suy yếu lực liên kết kim loại giữa các nguyên tử (cơ chế HEDE - Hydrogen Enhanced Decohesion) và tạo điều kiện cho biến vị trượt dễ dàng hơn (cơ chế HELP - Hydrogen Enhanced Localized Plasticity), dẫn đến hiện tượng giòn hóa và phá hủy đột ngột.

Kết luận

Đồ án đã giải quyết bài toán tính toán hệ số cường độ ứng suất $K_I$ tại đỉnh vết nứt trên ống trao đổi nhiệt bình sinh hơi nhà máy điện hạt nhân VVER-1000. Bằng việc kết hợp chặt chẽ giữa lý thuyết giải tích cơ học rạn nứt đàn hồi tuyến tính, động học ăn mòn điện hóa do tích tụ Hydro và mô phỏng số phần tử hữu hạn chuyên sâu trên ANSYS Workbench 15.0, công trình đã xác lập một quy trình đánh giá độ toàn vẹn kết cấu chính xác với sai số dưới $3{,}2%$. Kết quả nghiên cứu đóng góp cơ sở dữ liệu kỹ thuật giá trị cho công tác chẩn đoán an toàn, tối ưu hóa kế hoạch bảo dưỡng phòng ngừa và nâng cao hiệu quả vận hành kinh tế cho các nhà máy điện hạt nhân.