Giới thiệu dự án
Trong các nhà máy điện hạt nhân (NMĐHN) thế hệ III sử dụng công nghệ lò phản ứng nước áp lực VVER-1000 (công suất nhiệt 3000 MWth, công suất điện 1060 MWe), bình sinh hơi (Steam Generator - SG) đóng vai trò là ranh giới an toàn then chốt ngăn cách giữa vòng tuần hoàn sơ cấp mang phóng xạ và vòng thứ cấp sinh hơi quay tua bin. Bình sinh hơi VVER-1000 nằm ngang chứa 10.978 ống trao đổi nhiệt bằng thép không gỉ austenitic với đường kính 16 × 1,5 mm. Các ống này phải vận hành liên tục dưới điều kiện cực kỳ khắc nghiệt: áp suất vòng sơ cấp 15,7 MPa, áp suất hơi vòng thứ cấp 6,27 MPa, nhiệt độ dòng chất làm mát đầu vào lên tới 321 °C (594 K), cùng với bức xạ neutron và môi trường hóa nước sinh ra các gốc oxy hóa nguy hiểm ($H_2O_2$, $OH^-$, $O_2$, $e_{aq}^-$).
Tác nhân ăn mòn (SCC/IGSCC) Tích tụ Hydro & Ứng suất
Khởi sinh và phát triển vết nứt tại đỉnh
Tính toán SIF (K_I) qua ANSYS 15.0 & Giải tích
Vấn đề thực tiễn và bài toán kỹ thuật
Hiện tượng ăn mòn dưới tác dụng của ứng suất và môi trường (Stress Corrosion Cracking - SCC), kết hợp với sự giòn hóa do tích tụ hydro (Hydrogen Embrittlement), là nguyên nhân hàng đầu dẫn đến sự khởi sinh và lan truyền vết nứt tế vi tại bề mặt và liên kết hạt (IGSCC) trên thành ống trao đổi nhiệt. Một sự cố nứt vỡ ống trao đổi nhiệt không chỉ gây dừng lò khẩn cấp (thiệt hại ước tính từ 800.000 đến 1.200.000 USD/ngày vận hành) mà còn làm gia tăng nguy cơ rò rỉ phóng xạ từ vòng sơ cấp sang vòng thứ cấp. Việc xác định chính xác trường ứng suất và Hệ số cường độ ứng suất (Stress Intensity Factor - SIF, ký hiệu $K_I$) tại đỉnh vết nứt là yêu cầu bắt buộc để dự đoán tuổi thọ còn lại và thiết lập chu kỳ bảo dưỡng định kỳ.
Mục tiêu của dự án
- Phân tích cơ chế ăn mòn hóa lý: Thiết lập mô hình động học quá trình phân hủy phóng xạ của nước, cơ chế hòa tan trượt (Slip Dissolution Model) và quy luật khuếch tán, tích tụ hydro trong thép austenitic theo định luật Fick.
- Xây dựng mô hình giải tích vết nứt: Thiết lập phương trình trường ứng suất và hệ số cường độ ứng suất $K_I$ dạng giải tích theo bài toán Westergaard và các công thức bán thực nghiệm cho tấm phẳng hữu hạn chứa vết nứt biên.
- Mô phỏng phần tử hữu hạn (FEA): Triển khai mô hình 2D/3D chi tiết hình học vết nứt trên phần mềm ANSYS Workbench phiên bản 15.0 để trích xuất giá trị $K_I$ tại đỉnh vết nứt với kỹ thuật chia lưới kỳ dị (Singular crack tip meshing).
- Đối chuẩn và đánh giá sai số: Định lượng mối quan hệ tương quan giữa nồng độ tích tụ hydro, chiều sâu vết nứt ($a$) và giá trị $K_I$, đồng thời so sánh độ chính xác giữa phương pháp số (FEA) và phương pháp lý thuyết/bán thực nghiệm.
Phân tích và thiết kế giải pháp
Phân tích hiện trạng
| Phương pháp đánh giá |
Ưu điểm |
Nhược điểm |
Sai số ước tính |
Chi phí triển khai |
| Thực nghiệm kéo nứt & Đo dòng Anode |
Phản ánh chính xác điều kiện vi cấu trúc thực tế |
Chi phí cực cao, khó tái lập môi trường bức xạ/nhiệt áp cao |
Cơ sở chuẩn ($\pm 1,5%$) |
Rất cao (>50.000 USD/mẫu) |
| Giải tích Westergaard cổ điển |
Tính toán tức thời, công thức đóng đơn giản |
Chỉ đúng cho tấm phẳng vô hạn, bỏ qua hiệu ứng biên và phi tuyến |
$8% - 15%$ |
Rất thấp (thực hiện thủ công) |
| Công thức bán thực nghiệm Newman-Raju |
Có xét đến hệ số hiệu chỉnh hình học $\beta(a/w)$ |
Bị giới hạn bởi các dải tỷ lệ hình học chuẩn |
$4% - 7%$ |
Thấp (Lập trình tính toán) |
| Mô phỏng FEA (ANSYS Workbench 15.0) |
Mô phỏng chính xác hình học phức tạp, phân bố ứng suất đa trục |
Đòi hỏi kỹ năng chia lưới đỉnh vết nứt và tài nguyên tính toán |
$1,5% - 3,5%$ |
Trung bình (Bản quyền phần mềm) |
Yêu cầu hệ thống và ràng buộc kỹ thuật (MoSCoW)
- Must-have: Mô hình hóa chính xác trạng thái biến dạng phẳng (Plane Strain) tại đỉnh vết nứt trên ống thép austenitic; tính toán được phân bố $K_I$ theo chiều dài nứt từ $a/w = 0,1$ đến $a/w = 0,6$.
- Should-have: Tích hợp hàm suy giảm đặc tính cơ học do nồng độ Hydro $C_H(t)$ theo thời gian vận hành lò phản ứng.
- Could-have: Tự động hóa quá trình chia lưới vùng lân cận vết nứt bằng APDL (ANSYS Parametric Design Language).
- Won't-have: Mô phỏng tương tác đa trường đầy đủ 3D Fluid-Structure Interaction (FSI) tức thời trong khóa luận này.
Thiết kế hệ thống tính toán và cơ sở phương pháp luận
Ngăn xếp công nghệ và công cụ kỹ thuật
- Nền tảng mô phỏng chính: ANSYS Workbench 15.0 (Module Mechanical / Fracture Mechanics).
- Bộ giải số: Preconditioned Conjugate Gradient (PCG) Solver với phần tử 8-nút bậc hai (PLANE183).
- Thuật toán trích xuất vết nứt: Tích phân J-Integral và Phương pháp ngoại suy độ mở chuyển vị đỉnh vết nứt (Displacement Extrapolation Method).
- Công cụ tính toán bổ trợ: MATLAB R2016a / Python 3.8 (NumPy, SciPy) xử lý tích phân phương trình khuếch tán Hydro và định luật Faraday.
Cơ sở toán học và cơ học rạn nứt
-
Trường ứng suất vùng đỉnh vết nứt theo Westergaard (Dạng mở Mode I):
$$\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]$$
$$\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]$$
$$\tau_{xy} = \frac{K_I}{\sqrt{2\pi r}} \cos\left(\frac{\theta}{2}\right)\sin\left(\frac{\theta}{2}\right)\cos\left(\frac{3\theta}{2}\right)$$
-
Hệ số cường độ ứng suất cho tấm phẳng hữu hạn có vết nứt biên:
$$K_I = Y \cdot \sigma \sqrt{\pi a}$$
Trong đó hàm hình học $Y = f(a/w)$ được xác định theo chuỗi đa thức bậc 4:
$$Y = 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$$
-
Mô hình hòa tan trượt (Slip Dissolution Model theo Định luật Faraday):
Tốc độ lan truyền trung bình của vết nứt $V_t$:
$$V_t = \frac{M}{\rho \cdot z \cdot F} \frac{Q_f}{\epsilon_f} \dot{\epsilon}{ct}$$
Trong đó $M$ là khối lượng mol nguyên tử kim loại, $\rho$ là khối lượng riêng, $F = 96485 \text{ C/mol}$ là hằng số Faraday, $z$ là số electron hóa trị giải phóng, $Q_f$ là mật độ điện tích oxy hóa, và $\dot{\epsilon}{ct}$ là tốc độ biến dạng tại đỉnh vết nứt.
Implementation và kết quả
Quy trình tính toán và mô hình hóa chi tiết
Quy trình mô phỏng trên ANSYS Workbench 15.0 được chuẩn hóa qua 5 giai đoạn:
- Khởi tạo thông số vật liệu: Thép không gỉ austenitic với mô đun đàn hồi $E = 195 \times 10^3 \text{ MPa}$, hệ số Poisson $\nu = 0,3$, giới hạn chảy $\sigma_{ys} = 240 \text{ MPa}$, mật độ $\rho = 7900 \text{ kg/m}^3$.
- Xây dựng hình học vết nứt: Khởi tạo mô hình thành ống trao đổi nhiệt bề dày $w = 1,5 \text{ mm}$, chiều dài mẫu chuẩn $2H = 20 \text{ mm}$, vết nứt biên dạng nứt đơn Mode I với chiều sâu $a$ thay đổi từ $0,15 \text{ mm}$ đến $0,90 \text{ mm}$ ($a/w = 0,1 \div 0,6$).
- Kỹ thuật chia lưới kỳ dị (Quarter-Point Singular Elements): Tại vùng đỉnh vết nứt (Crack Tip), các phần tử 2D PLANE183 được chuyển đổi thành phần tử kỳ dị bằng cách dời các nút giữa cạnh về vị trí $1/4$ khoảng cách tính từ đỉnh vết nứt. Kỹ thuật này giúp tái lập chính xác đặc tính kỳ dị bậc $1/\sqrt{r}$ của trường ứng suất cơ học gãy đàn hồi tuyến tính (LEFM).
! APDL Script Snippet for Crack Tip Meshing in ANSYS Mechanical APDL
/PREP7
ET,1,PLANE183,1 ! Element type 8-node plane strain
MP,EX,1,195E3 ! Young's Modulus: 195 GPa
MP,PRXY,1,0.3 ! Poisson's Ratio: 0.3
K,1,0,0,0 ! Crack Tip Node
K,2,0.15,0,0 ! Crack Mouth Node
L,1,2
KSCON,1,0.005,1,12 ! Concentrate elements at Crack tip (Quarter-point node shift)
AMESH,ALL
/SOLU
ANTYPE,0 ! Static structural analysis
SFL,ALL,PRES,-9.43 ! Apply net differential pressure (MPa)
SOLVE
- Khuếch tán và tích tụ Hydro: Tính toán theo phương trình vi phân Fick 2 với hệ số khuếch tán $D = D_0 \exp\left(-\frac{E_d}{k_B T}\right)$ với $D_0 = 1,2 \times 10^{-7} \text{ m}^2/\text{s}$, $E_d = 0,54 \text{ eV}$. Nồng độ tích tụ hydro $C_H(t)$ gia tăng theo chu kỳ làm việc, làm giảm năng lượng bề mặt liên kết $\gamma_s$ và tăng độ nhạy cảm ứng suất.
# Script tính toán nồng độ Hydro tích tụ và suy giảm giới hạn bền
import numpy as np
def calculate_hydrogen_accumulation(time_hours, D_coeff, C_0):
"""
Tính nồng độ Hydro tích tụ tại đỉnh vết nứt theo thời gian vận hành
"""
t_seconds = time_hours * 3600.0
# Nghiệm hàm sai số bù (erfc) cho bài toán bán vô hạn
x = 1.5e-3 # Độ dày thành ống (m)
c_profile = C_0 * (1.0 - (x / (2.0 * np.sqrt(D_coeff * t_seconds + 1e-12))))
return np.maximum(c_profile, C_0)
time_steps = np.array([1000, 5000, 10000, 20000, 40000]) # Giờ vận hành
D_austenitic = 1.2e-7 * np.exp(-0.54 * 1.602e-19 / (1.38e-23 * 594.0)) # m^2/s tại 321°C
print("D_eff at 321C:", D_austenitic)
Kết quả thử nghiệm và kiểm chứng mô hình
Bảng so sánh giá trị Hệ số cường độ ứng suất $K_I$ ($\text{MPa}\cdot\sqrt{\text{m}}$) thu được giữa Phương pháp Giải tích Westergaard, Bán thực nghiệm và Mô phỏng ANSYS FEA 15.0 dưới ứng suất danh nghĩa $\sigma = 120 \text{ MPa}$:
| Chiều sâu nứt $a$ (mm) |
Tỷ lệ $a/w$ |
$K_I$ Westergaard ($\text{MPa}\sqrt{\text{m}}$) |
$K_I$ Bán thực nghiệm ($\text{MPa}\sqrt{\text{m}}$) |
$K_I$ ANSYS FEA 15.0 ($\text{MPa}\sqrt{\text{m}}$) |
Sai số FEA vs Bán thực nghiệm (%) |
| 0,15 |
0,10 |
2,605 |
2,918 |
2,845 |
$-2,50%$ |
| 0,30 |
0,20 |
3,684 |
4,237 |
4,120 |
$-2,76%$ |
| 0,45 |
0,30 |
4,512 |
5,460 |
5,310 |
$-2,74%$ |
| 0,60 |
0,40 |
5,210 |
6,825 |
6,630 |
$-2,85%$ |
| 0,75 |
0,50 |
5,825 |
8,563 |
8,290 |
$-3,18%$ |
| 0,90 |
0,60 |
6,381 |
10,912 |
10,480 |
$-3,95%$ |
K_I (MPa√m)
0.10 0.20 0.30 0.40 0.50 0.60 a/w
Đánh giá định lượng kết quả đạt được
- Độ chính xác mô hình số: Kết quả trích xuất $K_I$ từ ANSYS 15.0 cho độ lệch chỉ từ $2,50%$ đến $3,95%$ so với công thức bán thực nghiệm Newman-Raju, chứng minh tính tin cậy tuyệt đối của phương pháp phần tử hữu hạn có hiệu chỉnh phần tử kỳ dị.
- Tác động của tỷ lệ $a/w$: Khi vết nứt phát triển từ $a/w = 0,1$ lên $a/w = 0,6$, hệ số $K_I$ tăng phi tuyến tới $268%$ (từ $2,845 \text{ MPa}\sqrt{\text{m}}$ lên $10,480 \text{ MPa}\sqrt{\text{m}}$).
- Mối đe dọa từ tích tụ Hydro: Dữ liệu tính toán cho thấy khi nồng độ hydro nguyên tử đạt ngưỡng bão hòa trong mạng tinh thể FCC ($> 5 \text{ ppm}$ sau 30.000 giờ làm việc), giá trị ngưỡng chịu nứt ăn mòn $K_{ISCC}$ của thép không gỉ bị suy giảm xấp xỉ $35% - 42%$, thúc đẩy vết nứt chuyển nhanh từ giai đoạn 1 sang giai đoạn 2 (ổn định) với tốc độ tăng trưởng vết nứt $V_t \approx 10^{-9} \div 10^{-8} \text{ m/s}$.
Đổi mới và đóng góp
- Chuẩn hóa quy trình liên kết Động học ăn mòn - Cơ học rạn nứt: Tích hợp thành công mô hình giải phóng điện tích anode $Q_f$ từ mô hình phân rã màng thụ động (Film Rupture) với trường tập trung ứng suất cơ học gãy đàn hồi tuyến tính (LEFM).
- Kỹ thuật chia lưới tối ưu hóa tài nguyên: Thiết lập cấu hình lưới đa vùng (Multi-zone Mesh) kết hợp nút dịch chuyển $1/4$ tại đỉnh vết nứt trên ANSYS Workbench 15.0, giúp giảm $60%$ số lượng phần tử cần thiết mà vẫn đảm bảo độ hội tụ và độ chính xác của tích phân J-integral trong phạm vi sai số $< 3%$.
- Bộ dữ liệu tham chiếu chuyên biệt cho VVER-1000: Cung cấp đường cong tương quan chuẩn giữa chiều sâu vết nứt ($a/w$), nồng độ Hydro tích tụ và hệ số $K_I$ trực tiếp cho ống trao đổi nhiệt vật liệu thép không gỉ austenitic trong môi trường lò áp lực VVER-1000, thay thế cho các giả định gần đúng từ các dòng lò PWR phương Tây (Westinghouse/Framatome).
Ứng dụng thực tế và triển khai
Bối cảnh ứng dụng thực tế
- Quản lý tuổi thọ nhà máy (Plant Life Management - PLiM): Cung cấp tiêu chuẩn tính toán độ an toàn kết cấu cho phép cơ quan kiểm định hạt nhân đưa ra quyết định nút ống (tube plugging) hoặc thay thế từng cụm ống trao đổi nhiệt trong các kỳ đại tu (Refueling Outage).
- Thiết lập chu kỳ kiểm tra không phá hủy (NDT/ISI): Dựa trên tốc độ lan truyền vết nứt tính từ giá trị $K_I$, nhóm vận hành có thể xác định chính xác tần suất kiểm tra dòng điện xoáy (Eddy Current Testing - ECT) định kỳ cho 10.978 ống trong bình sinh hơi.
[Kiểm tra ECT định kỳ]
[So sánh K_I với ngưỡng K_ISCC]
K_I < 0.6 * K_ISCC K_I >= 0.6 * K_ISCC
[Cho phép vận hành tiếp tục] [Nút ống bảo vệ (Tube Plugging)]
Phân tích hiệu quả kinh tế và độ khả thi
- Ngăn chặn dừng lò đột xuất: Một sự cố rò rỉ ống bình sinh hơi buộc phải dừng lò mất tối thiểu 3-5 ngày để xác định vị trí và nút ống, tiêu tốn khoảng $3.000.000 - 5.000.000 \text{ USD}$. Ứng dụng mô hình tính toán SIF giúp chuyển đổi từ bảo trì bị động sang bảo trì dự đoán (Predictive Maintenance), tiết kiệm hàng triệu USD chi phí vận hành hàng năm.
- Tiêu chuẩn triển khai: Phương pháp luận hoàn toàn tương thích với các bộ tiêu chuẩn kiểm tra an toàn hạt nhân quốc tế như ASME Boiler and Pressure Vessel Code Section XI (Rules for Inservice Inspection of Nuclear Power Plant Components) và tiêu chuẩn PNAE G-7 của Nga.
Hạn chế và hướng phát triển
Hạn chế kỹ thuật hiện tại
- Giới hạn đàn hồi tuyến tính (LEFM): Nghiên cứu chủ yếu áp dụng cơ học gãy đàn hồi tuyến tính, chưa mở rộng sang cơ học gãy đàn dẻo (Elastic-Plastic Fracture Mechanics - EPFM / Đường cong kháng nứt J-R curve) khi vật liệu chịu quá tải dẻo cục bộ lớn tại đỉnh vết nứt.
- Mô hình 2D phẳng: Mô hình mô phỏng xem xét trạng thái biến dạng phẳng 2D, chưa xét tới hình học vết nứt bán elip 3D không đối xứng theo chu vi ống trao đổi nhiệt.
- Bỏ qua hiệu ứng phân tách do bức xạ (RIS): Chưa ghép nối trực tiếp trường bức xạ neutron neutron flux $\phi$ với sự biến đổi thành phần hợp kim (giảm Crom cục bộ tại biên hạt do bức xạ).
Định hướng nghiên cứu mở rộng
- Phát triển mô hình vết nứt 3D bề mặt (Surface Semi-Elliptical Crack): Ứng dụng công nghệ XFEM (Extended Finite Element Method) trên ANSYS để mô phỏng sự lan truyền vết nứt 3D động dưới tác động của tải rung động thủy nhiệt (Flow-Induced Vibration - FIV).
- Mô hình đa trường vật lý tương tác (Coupled Multi-Physics): Tích hợp liên kết hai chiều giữa module CFD (mô phỏng dòng chảy hai pha sôi ở vòng thứ cấp) và module Cơ học cấu trúc (Transient Structural Mechanics) để khảo sát ảnh hưởng của dao động nhiệt (Thermal fatigue).
Đối tượng hưởng lợi
- Kỹ sư cơ học hạt nhân & Nhà nghiên cứu vật liệu: Tiếp cận phương pháp luận chuẩn xác kết hợp giữa mô phỏng số FEA bằng ANSYS 15.0 và lý thuyết gãy cơ học, làm nền tảng phát triển các mô hình đánh giá an toàn kết cấu nâng cao.
- Nhà máy điện hạt nhân và Cơ quan quản lý an toàn: Có cơ sở khoa học để đánh giá tuổi thọ bình sinh hơi VVER-1000, tối ưu hóa kế hoạch bảo trì ngăn chặn sự cố nứt vỡ ống trao đổi nhiệt, nâng cao hệ số khả dụng (Capacity Factor) của nhà máy.
- Sinh viên và học viên cao học ngành Kỹ thuật Hạt nhân/Cơ khí: Nguồn tài liệu học thuật hoàn chỉnh về cơ chế ăn mòn SCC, tương tác Hydro trong kim loại và ứng dụng công cụ mô phỏng phần tử hữu hạn trong công nghiệp năng lượng cao.
Câu hỏi thường gặp
1. Yêu cầu cấu hình hệ thống tối thiểu để chạy mô phỏng ANSYS Workbench 15.0 cho bài toán này là gì?
Mô hình 2D biến dạng phẳng tập trung phần tử kỳ dị yêu cầu cấu hình tối thiểu: CPU Intel Core i5 (4 nhân), RAM 8 GB, card đồ họa hỗ trợ OpenGL, dung lượng ổ cứng khả dụng 20 GB. Đối với mô hình mở rộng 3D dạng vết nứt bán elip (Semi-elliptical crack), cấu hình khuyến nghị là CPU Intel Xeon hoặc Core i7/i9 (8-16 nhân), RAM tối thiểu 32 GB và card đồ họa chuyên dụng NVIDIA Quadro.
2. Giới hạn độ sâu vết nứt an toàn trên ống trao đổi nhiệt bình sinh hơi VVER-1000 là bao nhiêu?
Theo tiêu chuẩn an toàn kỹ thuật của lò VVER-1000, giới hạn chiều sâu vết nứt tối đa cho phép trước khi bắt buộc phải nút ống (tube plugging) là $40%$ bề dày danh định của thành ống (tương đương $a = 0,6 \text{ mm}$ trên ống $16 \times 1,5 \text{ mm}$). Khi $a/w \ge 0,4$, hệ số $K_I$ bắt đầu tiến sát ngưỡng tới hạn gây giòn mỏi $K_{ISCC}$ dưới điều kiện tích tụ Hydro.
3. Tại sao trong thép không gỉ austenitic, hàm lượng Crom $\ge 11%$ lại có khả năng chống ăn mòn và hiện tượng IGSCC xảy ra do đâu?
Hàm lượng Crom từ $11%$ trở lên (thường là $17 - 18%$ trong thép austenitic) tạo ra một lớp màng oxit thụ động $(Cr_2O_3)$ siêu mỏng (từ $10$ đến $50 \text{ \AA}$) ngăn cách kim loại với môi trường oxy hóa. Tuy nhiên, khi kim loại bị nung nhiệt trong khoảng $450 - 850 \text{ °C}$ hoặc do ứng suất nhiệt kéo dài, các nguyên tử Carbon liên kết với Crom tạo kết tủa cacbit crom ($Cr_{23}C_6$) dọc theo ranh giới hạt, tạo ra "vùng nghèo Crom" ($< 12% \text{ Cr}$). Vùng này mất tính thụ động và trở thành anode điện hóa bị hòa tan nhanh chóng, gây ra hiện tượng ăn mòn nứt liên kết hạt (IGSCC).
4. Phương pháp trích xuất hệ số cường độ ứng suất $K_I$ trong ANSYS Workbench hoạt động như thế nào?
ANSYS Workbench sử dụng hai phương pháp chính: (1) Tích phân J-Integral: Tính toán tốc độ giải phóng năng lượng biến dạng trên các đường bao khép kín quanh đỉnh vết nứt, sau đó chuyển đổi sang SIF qua công thức $K_I = \sqrt{J \cdot E / (1 - \nu^2)}$ đối với trạng thái biến dạng phẳng; (2) Ngoại suy chuyển vị (Displacement Extrapolation): Đo độ mở chuyển vị tương đối $u_y$ của các nút tại bề mặt vết nứt gần đỉnh và sử dụng phương trình giải tích đàn hồi để ngoại suy về vị trí $r \to 0$.
5. Hydro làm thế nào thâm nhập vào cấu trúc thép và gây suy giảm SIF?
Trong môi trường nước nhiệt độ cao bị phân hủy phóng xạ, các phản ứng catot sinh ra hydro nguyên tử ($H^+$ nhận electron trở thành $H$). Do kích thước nguyên tử cực nhỏ, hydro khuếch tán vào mạng tinh thể thép austenitic theo định luật Fick và bị bẫy tại các vị trí khuyết tật mạng (lỗ trống Schottky/Frenkel, biến vị, biên hạt). Sự tích tụ hydro làm suy giảm năng lượng liên kết nguyên tử kim loại (cơ chế HEDE - Hydrogen-Enhanced Decohesion) và tăng độ linh động của biến vị (cơ chế HELP - Hydrogen-Enhanced Localized Plasticity), làm giảm giá trị năng lượng phá hủy bề mặt và hạ thấp ngưỡng ứng suất tới hạn $K_{ISCC}$.
Kết luận
Đồ án đã giải quyết bài toán cơ học rạn nứt kết cấu trong bình sinh hơi nhà máy điện hạt nhân VVER-1000:
- Làm sáng tỏ cơ chế nứt phá hủy: Xác lập mối liên hệ giữa các quá trình hóa lý vùng đỉnh vết nứt (ăn mòn SCC, phân hủy màng thụ động, khuếch tán Hydro) với sự suy giảm độ bền cơ học của thép austenitic.
- Làm chủ phương pháp số hiện đại: Ứng dụng thành công phần mềm ANSYS Workbench 15.0 với kỹ thuật phần tử kỳ dị bậc hai để tính toán chính xác Hệ số cường độ ứng suất $K_I$, đạt độ tương thích cao với sai số chỉ dưới $4%$ so với công thức bán thực nghiệm.
- Giá trị thực tiễn vượt trội: Cung cấp cơ sở định lượng quan trọng cho công tác kiểm tra an toàn định kỳ, thiết lập tiêu chuẩn kỹ thuật loại bỏ khuyết tật và tối ưu hóa chi phí vận hành an toàn cho các nhà máy điện hạt nhân.
