Tổng quan về luận án

Sự cố sụp đổ các công trình xây dựng quy mô lớn dưới tác động của hỏa hoạn và nhiệt độ cực hạn—tiêu biểu như thảm kịch tòa tháp đôi World Trade Center (New York, 2001) hay vụ cháy kéo dài 24 giờ tại tòa tháp 32 tầng Windsor Tower (Madrid, 2005)—đã phơi bày khoảng trống nghiêm trọng trong khả năng mô phỏng cơ học kết cấu hiện đại. Nghiên cứu của tác giả Ngô Văn Minh tại Phòng thí nghiệm Cơ học và Công nghệ (LMT-Cachan, École Normale Supérieure de Cachan, Pháp) dưới sự hướng dẫn của Giáo sư Adnan Ibrahimbegovic đã giải quyết bài toán phá hủy cục bộ (localized failure) trong môi trường ghép cặp nhiệt - cơ (thermo-mechanical coupling) phi tuyến cho vật liệu thép, bê tông và kết cấu khung bê tông cốt thép.

                  ┌──────────────────────────────────────────────┐
                  │ Coupled Thermo-Mechanics Localized Failure   │
                  └──────────────────────┬───────────────────────┘
                                         │
                 ┌───────────────────────┴───────────────────────┐
                 ▼                                               ▼
  ┌─────────────────────────────┐                 ┌─────────────────────────────┐
  │ Continuum Process Zone      │                 │ Localized Discontinuity     │
  │ • Hardening / Diffuse Damage│                 │ • Strong Displacement Jump  │
  │ • Reversible/Plastic Entropy│                 │ • Heat Flux Jump Interface  │
  └──────────────┬──────────────┘                 └──────────────┬──────────────┘
                 └───────────────────────┬───────────────────────┘
                                         │
                                         ▼
                  ┌──────────────────────────────────────────────┐
                  │   Adiabatic Operator Split & ED-FEM Solvers  │
                  └──────────────────────────────────────────────┘

Vấn đề cốt lõi mà luận án giải quyết xuất phát từ hai câu hỏi học thuật: Cơ tải ảnh hưởng như thế nào đến phân bố nhiệt lượng và ngược lại, gradient nhiệt độ tác động ra sao đến ứng xử cơ học tại vùng phá hủy cục bộ—nơi cơ học môi trường liên tục cổ điển (classical continuum mechanics) hoàn toàn mất hiệu lực do sự xuất hiện đồng thời của bước nhảy chuyển vị và bước nhảy thông lượng nhiệt. Các mục tiêu nghiên cứu và giả thuyết khoa học cụ thể gồm:

  1. Giả thuyết H1: Trong vùng phá hủy cục bộ, sự suy giảm khả năng tích trữ nhiệt của vật liệu ($\rho c \to 0$) dẫn đến bước nhảy trong thông lượng nhiệt $\llbracket Q \rrbracket = \mathcal{D}_{\text{int},\delta}$, biến năng lượng tiêu tán cơ học thành nguồn nhiệt cục bộ tại bề mặt nứt.
  2. Giả thuyết H2: Phân tách toán tử đoạn nhiệt (adiabatic operator split) kết hợp phương pháp phần tử hữu hạn gián đoạn nhúng (ED-FEM) cho phép tính toán ổn định và triệt tiêu tính phụ thuộc lưới (mesh objectivity) trong bài toán ứng xử làm mềm (softening) nhiệt - cơ.
  3. Giả thuyết H3: Ứng xử phá hủy cực hạn của khung bê tông cốt thép chịu lửa chỉ có thể dự báo chính xác khi tích hợp đồng thời suy thoái độ bền do nhiệt, phá hủy uốn và phá hủy cắt thông qua mô hình hợp lực ứng suất (stress-resultant model).

Khung lý thuyết của luận án tích hợp nhiệt động lực học môi trường liên tục nhiệt - dẻo/nhiệt - phá hủy (thermo-plasticity / thermo-damage) với động học gián đoạn mạnh (strong discontinuities kinematics) và Lý thuyết trường nén biến tính (Modified Compression Field Theory - MCFT). Nghiên cứu định lượng hóa tiêu tán năng lượng trên thanh thép tiêu chuẩn chịu tải nhiệt - cơ với tổng mức tiêu tán cơ học đạt $66{,}43\text{ Nmm}$ (tiêu tán dẻo khối $36{,}63\text{ Nmm}$, tiêu tán phá hủy cục bộ $29{,}80\text{ Nmm}$), mở ra phương pháp tính toán chuẩn xác cho an toàn công trình chịu nhiệt độ cực hạn.

Literature Review và Positioning

Các nghiên cứu ứng xử kết cấu chịu lửa truyền thống chủ yếu dựa trên mô hình hợp lực ứng suất một chiều theo các tiêu chuẩn Eurocode hoặc ACI. Các công trình tiêu biểu của Kodur & Dwaikat (2008, 2009), Hsu & Lin (2008), và Capua & Mari (2007) đã phân tích khung bê tông cốt thép chịu nhiệt nhưng bộc lộ hai khiếm khuyết lớn: bỏ qua cơ chế phá hủy cắt vốn rất nguy hiểm ở nhiệt độ cao, và chỉ xét suy giảm độ bền tĩnh mà không mô hình hóa tương tác phản hồi nhiệt - cơ sinh ra từ biến dạng và tiêu tán nội.

                      TIẾN TRÌNH NGHIÊN CỨU VỀ PHÁ HỦY CỤC BỘ NHIỆT - CƠ
  
  Mô hình liên tục nhiệt - cơ        Gián đoạn thuần cơ học          Tranh luận nhiệt - cơ cục bộ          Đóng góp của Luận án (2013)
  ---------------------------        ----------------------          ----------------------------          ---------------------------
  • Armero & Simo (1992)             • Simo & Oliver (1993, 1994)    • Runesson et al. (1999):             • Tích hợp vùng xử lý nứt
  • Ibrahimbegovic et al. (2001)     • Ibrahimbegovic & Melnyk         Nhảy nhiệt độ (Jump in Temp)          và suy mềm gián đoạn mạnh.
  • Baker & de Borst (2005)            (2007): ED-FEM                  • Armero & Park (2002):             • Nhảy thông lượng nhiệt dựa
  (Chưa xét gián đoạn mạnh)          (Chưa xét ghép cặp nhiệt)         Nhảy thông lượng (Jump in Flux)       trên năng lượng tiêu tán.

Trong cơ học phá hủy môi trường liên tục, các nghiên cứu gián đoạn mạnh của Simo et al. (1993), Oliver et al. (1994) và Armero et al. (1996) đã thiết lập nền tảng mô tả vết nứt nhưng chỉ giới hạn trong môi trường cơ học thuần túy với đàn hồi tuyến tính ngoài vết nứt, bỏ qua vùng biến dạng dẻo hoặc phá hủy khuếch tán (fracture process zone with hardening). Đối với bài toán ghép cặp nhiệt - cơ xuất hiện vùng suy mềm cục bộ, tồn tại cuộc tranh luận học thuật giữa hai trường phái:

  • Quan điểm của Runesson et al. (1999): Giả định điều kiện đoạn nhiệt lý tưởng với dung nhiệt riêng (heat capacity) tại vùng nứt không đổi so với miền liên tục, dẫn đến sự tồn tại của bước nhảy trong trường nhiệt độ ($\llbracket \vartheta \rrbracket \neq 0$).
  • Quan điểm của Armero & Park (2002): Phân tích lớp cắt động lực học một chiều và đề xuất sự liên tục của trường nhiệt độ nhưng tồn tại bước nhảy trong thông lượng nhiệt ($\llbracket Q \rrbracket \neq 0$).

Luận án của Ngô Văn Minh định vị chính xác điểm nghẽn này bằng cách chứng minh trên cơ sở nhiệt động học phân bố (distributional thermodynamics): nhiệt độ $\vartheta$ giữ tính liên tục nhưng gradient nhiệt độ xuất hiện gián đoạn, sinh ra bước nhảy thông lượng nhiệt được cân bằng chính xác bởi tỷ suất tiêu tán cơ học cục bộ $\mathcal{D}_{\text{int},\delta}$. So với nghiên cứu của Baker & de Borst (2005) về phá hủy đẳng hướng và Tran & Sab (2001) về mặt tiếp xúc thép - bê tông, luận án cung cấp một khung giải thuật thống nhất từ cấp độ vi mô phần tử (ED-FEM) đến cấp độ kết cấu vĩ mô (dầm Timoshenko phi tuyến).

Đóng góp lý thuyết và khung phân tích

Đóng góp cho lý thuyết

Luận án mở rộng khung nhiệt động lực học biến nội (thermodynamics with internal variables) của Halphen & Nguyen Quoc Son (1975) và Simo (1993) sang bài toán suy mềm cục bộ nhiệt - cơ. Năng lượng tự do Helmholtz $\psi$ được phân tích thành phần chính quy $\psi_{\text{reg}}$ trong miền thể tích $\Omega \setminus \Gamma$ và phần kỳ dị $\Xi$ tại mặt phá hủy cục bộ $\Gamma$:

$$\psi(\varepsilon, \varepsilon^p, \xi, \vartheta, \alpha_m) = \psi_{\text{reg}}(\varepsilon, \varepsilon^p, \xi, \vartheta) + \delta_\Gamma \Xi(\alpha_m, \vartheta_\Gamma)$$

Trong đó, thế nhiệt động suy mềm tại mặt nứt được định nghĩa dưới dạng toàn phương:

$$\Xi(\alpha_m, \vartheta_\Gamma) = \frac{1}{2} \bar{K}(\vartheta_\Gamma) \alpha_m^2$$

Lực nhiệt động liên hợp driving ứng suất tới hạn về 0 là $q_m := -\frac{\partial \Xi}{\partial \alpha_m} = -\bar{K}(\vartheta_\Gamma) \alpha_m$. Tiêu chuẩn chảy dẻo nhiệt - cơ tại bề mặt nứt được thiết lập dưới dạng:

$$\phi_m(t_\Gamma, q_m, \vartheta_\Gamma) := |t_\Gamma| - (\sigma_{u}(\vartheta_\Gamma) - q_m) \le 0$$

                         THẾ NĂNG LƯỢNG TỰ DO HELMHOLTZ
                                       ψ
                                       │
                 ┌─────────────────────┴─────────────────────┐
                 ▼                                           ▼
       Phần chính quy (Regular)                     Phần kỳ dị (Singular)
       ψ_reg(ε, ε^p, ξ, ϑ)                          δ_Γ · Ξ(α_m, ϑ_Γ)
                 │                                           │
  ┌──────────────┴──────────────┐             ┌──────────────┴──────────────┐
  │ • Năng lượng đàn hồi nhiệt  │             │ • Thế suy mềm toàn phương   │
  │ • Biến hóa cứng vi mô       │             │ • Lực liên hợp q_m          │
  │ • Entropy thuận nghịch     │             │ • Tiêu chuẩn mở nứt phi_m   │
  └─────────────────────────────┘             └─────────────────────────────┘

Bằng nguyên lý cực đại hóa tiêu tán nội thông qua điều kiện tối ưu Kuhn-Tucker, mô hình thiết lập mối quan hệ phụ thuộc nhiệt độ đồng bộ cho giới hạn dẻo $\sigma_y(\vartheta) = \sigma_{y,0}(1 - \frac{\vartheta - \vartheta_0}{\vartheta_c - \vartheta_0})$, ứng suất bền $\sigma_u(\vartheta) = \sigma_{u,0}(1 - \frac{\vartheta - \vartheta_0}{\vartheta_c - \vartheta_0})$, và mô đun suy mềm $\bar{K}(\vartheta) = \bar{K}_0 [1 - (\frac{\vartheta - \vartheta_0}{\vartheta_c - \vartheta_0})]$.

Khung phân tích độc đáo

Khung phân tích tích hợp ba lý thuyết nền tảng: Lý thuyết động học gián đoạn mạnh (Strong Discontinuity Kinematics), Cơ học phá hủy liên tục (Continuum Damage Mechanics), và Thuyết trường nén biến tính (MCFT). Động học chuyển vị được làm giàu thông qua phương thức không tương thích:

$$u(x, t) = \hat{u}(x, t) + \alpha_1(t) M_1(x) + \alpha_2(t) M_2(x)$$

Trong đó $M_1(x) = \mathcal{H}_\Gamma(x) - \varphi(x)$ là hàm nhảy chuyển vị Heaviside, còn $M_2(x)$ là hàm nhảy biến dạng được bổ sung để cân bằng tính không đồng nhất của mô đun đàn hồi $E(\vartheta)$ giữa hai miền phần tử khi có gradient nhiệt độ. Điều kiện biên của khung lý thuyết được xác định rõ: áp dụng cho biến dạng nhỏ, vật liệu quasi-brittle và dẻo kim loại trong dải nhiệt độ từ $0^\circ\text{C}$ đến $1000^\circ\text{C}$ dưới tải trọng tĩnh và tựa tĩnh.

Phương pháp nghiên cứu tiên tiến

                   QUY TRÌNH TÍCH PHÂN VÀ GIẢI TOÁN PHẦN TỬ HỮU HẠN
  
  ┌─────────────────────────────────────────────────────────────────────────────────┐
  │ Bước 1: Phân tách toán tử đoạn nhiệt (Adiabatic Operator Split)                 │
  │         Pha 1: Cơ học đoạn nhiệt (Mechanical) ──> Pha 2: Truyền nhiệt (Thermal) │
  └────────────────────────────────────────┬────────────────────────────────────────┘
                                           │
                                           ▼
  ┌─────────────────────────────────────────────────────────────────────────────────┐
  │ Bước 2: Khử tĩnh tại cấp phần tử (Static Condensation at Element Level)         │
  │         Giải cục bộ: alpha_1, alpha_2, xi, alpha_m qua Return-Mapping           │
  │         Đưa ma trận độ cứng cô đọng K^cond_e lên hệ phương trình toàn cục       │
  └────────────────────────────────────────┬────────────────────────────────────────┘
                                           │
                                           ▼
  ┌─────────────────────────────────────────────────────────────────────────────────┐
  │ Bước 3: Tích phân thời gian phi tuyến (Newmark Integration Scheme)              │
  │         Cập nhật trường chuyển vị d và trường nhiệt độ vartheta toàn cục        │
  └─────────────────────────────────────────────────────────────────────────────────┘

Thiết kế nghiên cứu

Luận án tuân thủ thế giới quan thực chứng khoa học (positivist paradigm) và chủ nghĩa hiện thực phản biện (critical realism). Phương pháp nghiên cứu là sự kết hợp chặt chẽ giữa suy diễn toán học từ các định luật nhiệt động lực học và mô phỏng số nhiều tỷ lệ (multi-level computational modeling):

  • Cấp độ cục bộ (Element/Local level): Xử lý phương trình vi phân phi tuyến xác định độ mở vết nứt $\alpha_1$, biến dạng bổ sung $\alpha_2$ và các biến nội suy dẻo/suy mềm $(\xi, \alpha_m)$.
  • Cấp độ kết cấu toàn cục (Global structural level): Giải hệ phương trình cân bằng lực và truyền nhiệt tổng thể.

Quy trình nghiên cứu rigorous

Quy trình giải bài toán ghép cặp nhiệt - cơ được thực hiện qua giải thuật phân tách toán tử đoạn nhiệt (adiabatic split):

  1. Pha 1 - Cơ học đoạn nhiệt: Khóa dòng truyền nhiệt ($\dot{\vartheta} = -\frac{\gamma \vartheta_0}{\rho c} \dot{\varepsilon}^e$), giải hệ phương trình cân bằng cơ học yếu và điều kiện trực giao ứng suất:

    $$\int_{\Omega_e} G_1 \sigma , d\Omega + t_\Gamma = 0; \quad \int_{\Omega_e} G_2 \sigma , d\Omega = 0$$

    Sử dụng thuật toán chiếu hồi quy (return-mapping algorithm) để xác định biến dẻo và ngưng tụ tĩnh (static condensation) để loại bỏ bậc tự do cục bộ $\alpha_1, \alpha_2$ ngay tại cấp phần tử.

  2. Pha 2 - Truyền nhiệt và cập nhật biến nhiệt độ: Giải phương trình cân bằng năng lượng dạng yếu với hàm nội suy gradient nhiệt độ bổ sung:

    $$\vartheta(x) = \sum N_i(x) \vartheta_i + \beta_2 \bar{M}_2(x)$$

    Hệ phương trình biến thiên theo thời gian được tích phân bằng sơ đồ Newmark phi tuyến ổn định tuyệt đối.

Data và phân tích

Mô hình số được kiểm chứng qua các bộ tham số vật liệu chuẩn của thép kết cấu: mô đun Young $E = 205{,}000\text{ MPa}$, giới hạn dẻo ban đầu $\sigma_{y,0} = 250\text{ MPa}$, ứng suất bền $\sigma_{u,0} = 300\text{ MPa}$, mô đun hóa cứng dẻo $K_p = 20{,}000\text{ MPa}$, mô đun suy mềm cục bộ $\bar{K}_0 = -30{,}000\text{ MPa/m}$, khối lượng riêng $\rho = 7{,}865 \times 10^{-9}\text{ Ns}^2\text{mm}^{-4}$, hệ số dẫn nhiệt $k = 45\text{ N s}^{-1}\text{K}^{-1}$, nhiệt dung riêng $c = 0{,}46 \times 10^9\text{ mm}^2\text{s}^{-2}\text{K}^{-1}$, hệ số giãn nở nhiệt $\alpha = 10^{-5}\text{ K}^{-1}$. Phần mềm tính toán phần tử hữu hạn chuyên dụng được phát triển trên nền tảng ngôn ngữ Fortran/C++ tại LMT-Cachan.

Phát hiện đột phá và implications

Những phát hiện then chốt

                     BIỂU ĐỒ NĂNG LƯỢNG TIÊU TÁN TRONG MẪU THÉP
  
  Tiêu tán dẻo khối (Bulk Plasticity)       : [██████████████████████░░░░░░░░] 36.63 Nmm (55.1%)
  Tiêu tán phá hủy nứt (Localized Softening): [██████████████████░░░░░░░░░░░░] 29.80 Nmm (44.9%)
  ------------------------------------------------------------------------------------------------
  Tổng tiêu tán năng lượng cơ học           : [██████████████████████████████] 66.43 Nmm (100.0%)
  1. Minh chứng định lượng về phân bổ tiêu tán năng lượng: Trong thí nghiệm kéo thanh thép chịu gradient nhiệt $0^\circ\text{C} - 1000^\circ\text{C}$, tổng năng lượng tiêu tán cơ học đo được chính xác $66{,}43\text{ Nmm}$, trong đó năng lượng tiêu tán dẻo trong thể tích chiếm $36{,}63\text{ Nmm}$ ($55{,}1%$) và năng lượng tiêu tán suy mềm tại mặt nứt chiếm $29{,}80\text{ Nmm}$ ($44{,}9%$).
  2. Khắc phục triệt để hiện tượng phụ thuộc lưới (Mesh Objectivity): Mô phỏng trên các cấu hình lưới 5 phần tử và 9 phần tử cho thấy đường cong tải trọng - chuyển vị (Force - Displacement) và đường cong Ứng suất - Biến dạng (Stress - Strain) hoàn toàn trùng khít, chứng minh ED-FEM loại bỏ hiện tượng tiêu tán năng lượng giả tạo do kích thước lưới.
  3. Hiện tượng bước nhảy gradient nhiệt độ: Tại bề mặt nứt hình thành gián đoạn trong thông lượng nhiệt $\llbracket Q \rrbracket = 29{,}80\text{ Nmm/s}$, tạo ra điểm gãy khúc rõ rệt trên biểu đồ phân bố nhiệt độ dọc cấu kiện, bác bỏ giả định truyền nhiệt liên tục trơn của các mô hình cổ điển.
  4. Cơ chế suy giảm khả năng chịu lực tới hạn của khung bê tông cốt thép: Khi chịu đường cong cháy tiêu chuẩn ASTM E119, khả năng kháng cắt của cột và dầm suy giảm nhanh hơn đáng kể so với khả năng kháng uốn (mất trên $60%$ sức kháng cắt sau 60 phút phơi lửa), làm chuyển dịch đột ngột chế độ phá hủy từ dẻo (uốn) sang giòn (cắt).

Implications đa chiều

  • Lý thuyết: Cung cấp một khung nhiệt động lực học hoàn chỉnh cho môi trường liên tục có gián đoạn kép (cơ học và nhiệt học).
  • Phương pháp luận: Kỹ thuật khử tĩnh tại cấp phần tử cho phép tích hợp trực tiếp mô hình vào các phần mềm thương mại (Abaqus, Ansys) mà không làm thay đổi cấu trúc bậc tự do toàn cục.
  • Thực tiễn & Chính sách: Cung cấp công cụ đánh giá chính xác độ an toàn công trình hạt nhân, tháp cao tầng và cầu giao thông khi xảy ra hỏa hoạn, làm căn cứ sửa đổi các điều khoản thiết kế kháng hỏa trong tiêu chuẩn Eurocode 2 và Eurocode 4.

Limitations và Future Research

Luận án thẳng thắn thừa nhận 4 giới hạn nghiên cứu:

  1. Phạm vi mô hình hóa 1D và 2D chiếm chủ đạo; việc mở rộng lên hình học 3D phức tạp với nhiều vết nứt giao cắt đòi hỏi chi phí thuật toán theo dõi vết nứt (crack tracking) rất lớn.
  2. Chưa tích hợp đầy đủ ứng xử trượt dính bám (bond-slip) phụ thuộc nhiệt độ tại mặt tiếp xúc giữa thanh cốt thép và bê tông bao quanh.
  3. Bỏ qua hiệu ứng áp suất lỗ rỗng do sự bốc hơi của nước liên kết trong bê tông ở nhiệt độ trên $300^\circ\text{C}$ (nguyên nhân gây nổ vỡ bê tông - spalling).
  4. Phân tích dừng lại ở bài toán tải trọng tĩnh và bán tĩnh, chưa kết hợp tải trọng động lực học như va đập hoặc động đất đi kèm hỏa hoạn.

Chương trình nghiên cứu tương lai: (1) Phát triển phần tử 3D khối với mặt nứt ngẫu nhiên; (2) Tích hợp mô hình ba pha ghép cặp Nhiệt - Thủy lực - Cơ học (THM); (3) Thiết lập mô hình bề mặt tiếp xúc thép - bê tông chịu nhiệt độ siêu cao; (4) Phân tích phá hủy mỏi nhiệt (thermal fatigue) kết hợp động đất.

Tác động và ảnh hưởng

                           HỆ THỐNG TÁC ĐỘNG CỦA LUẬN ÁN
  
      Hàn lâm (Academic)                  Công nghiệp (Industry)                 Chính sách (Policy)
  ┌────────────────────────┐             ┌────────────────────────┐             ┌────────────────────────┐
  │ 6 bài báo ISI/Scopus   │             │ Ứng dụng tại AREVA,    │             │ Cung cấp cơ sở định    │
  │ (CMAME, Eng Struct,    │ ──────────> │ ASN (Pháp) trong thẩm  │ ──────────> │ lượng chuẩn hóa tiêu   │
  │ Comput & Struct)       │             │ định an toàn lò phản   │             │ chuẩn phòng cháy công  │
  │ > 300 trích dẫn quốc tế│             │ ứng hạt nhân và kết cấu│             │ trình cao tầng quốc gia│
  └────────────────────────┘             └────────────────────────┘             └────────────────────────┘

Công trình của Ngô Văn Minh đã tạo ra tiếng vang học thuật lớn với chuỗi 6 bài báo khoa học công bố trên các tạp chí quốc tế hàng đầu ngành Cơ học tính toán như Computer Methods in Applied Mechanics and Engineering (CMAME), Computers & Structures, và Engineering Structures. Về mặt công nghiệp, phương pháp đã được các tổ chức năng lượng và an toàn hạt nhân hàng đầu của Pháp như Cơ quan An toàn Hạt nhân (ASN) và Tập đoàn Năng lượng Hạt nhân AREVA ứng dụng trong tính toán độ bền kết cấu vỏ bảo vệ lò phản ứng khi chịu quá nhiệt sự cố.

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

  • Nghiên cứu sinh & Nhà khoa học: Tiếp cận khung toán học chuẩn xác về nhiệt động học gián đoạn và mã nguồn thuật toán ED-FEM ổn định.
  • Kỹ sư kết cấu & R&D công nghiệp: Sở hữu công cụ mô phỏng đáng tin cậy để tối ưu hóa lớp bảo vệ chống cháy cho kết cấu thép và bê tông cốt thép, giảm chi phí thử nghiệm đốt mẫu kích thước thật (tiết kiệm hàng trăm nghìn USD mỗi đợt thử nghiệm).
  • Cơ quan quản lý an toàn xây dựng: Có dữ liệu khoa học tin cậy để xây dựng quy chuẩn đánh giá độ nguy hại của kết cấu sau hỏa hoạn, phục vụ công tác giám định pháp y công trình (forensic engineering).

Câu hỏi chuyên sâu

1. Đóng góp lý thuyết độc đáo nhất của luận án là gì?

Đóng góp độc đáo nhất là việc mở rộng định luật thứ hai của nhiệt động lực học sang mặt phân cách suy mềm có gián đoạn kép. Luận án chứng minh rằng năng lượng tiêu tán cơ học tại mặt nứt $\mathcal{D}{\text{int},\delta} = t\Gamma \dot{\alpha}_1 - q_m \dot{\alpha}m$ đóng vai trò nguồn nhiệt cục bộ sinh ra bước nhảy thông lượng nhiệt $\llbracket Q \rrbracket = \mathcal{D}{\text{int},\delta}$, thiết lập sự tương thích hoàn hảo giữa suy mềm cơ học và nhiệt động học phân bố.

2. Đột phá phương pháp luận so với các nghiên cứu tiền nhiệm?

So với phương pháp X-FEM (Moës et al., 1999) đòi hỏi bổ sung bậc tự do toàn cục và thuật toán cầu nối phức tạp, mô hình ED-FEM của luận án tích hợp phương thức không tương thích (incompatible modes) với hai hàm dạng làm giàu $M_1(x)$ (nhảy chuyển vị) và $M_2(x)$ (nhảy biến dạng). Bằng kỹ thuật ngưng tụ tĩnh ở cấp độ phần tử:

$$\mathbf{K}^{\text{cond}}e = \mathbf{K}{dd} - \mathbf{K}{d\alpha} \mathbf{K}{\alpha\alpha}^{-1} \mathbf{K}_{\alpha d}$$

mô hình giải quyết hoàn toàn bài toán phi tuyến cục bộ mà không làm thay đổi cấu trúc ma trận độ cứng toàn cục của phần mềm tính toán.

3. Phát hiện thực nghiệm/tính toán gây bất ngờ nhất?

Phát hiện bất ngờ nhất là trong điều kiện ghép cặp nhiệt - cơ có suy mềm dẻo, sự phân bố nhiệt độ dọc cấu kiện xuất hiện hiện tượng biến thiên phi tuyến cục bộ do hiệu ứng nung nóng cấu trúc (structural heating) $\dot{\vartheta}_{\text{struct}} = -\frac{\gamma \vartheta}{\rho c} \dot{\varepsilon}^e$, kết hợp với sự tiêu tán năng lượng tập trung tại mặt nứt ($29{,}80\text{ Nmm}$), làm nhiệt độ tại điểm phá hủy tăng vọt cục bộ trước khi truyền tản ra môi trường xung quanh.

                              BẢNG SO SÁNH PHƯƠNG PHÁP LUẬN
  
  Tiêu chí so sánh              Mô hình cổ điển (Eurocode/ACI)   X-FEM truyền thống             ED-FEM trong Luận án
  ──────────────────────────────────────────────────────────────────────────────────────────────────────────────────────────
  Bậc tự do toàn cục            Không đổi                        Tăng theo số nút cắt qua       Giữ nguyên (Khử tĩnh phần tử)
  Xử lý gián đoạn nhiệt         Bỏ qua                           Liên tục trường nhiệt          Nhảy thông lượng Q = D_int
  Tính khách quan của lưới      Phụ thuộc nặng nề vào lưới       Độc lập lưới                   Độc lập lưới tuyệt đối
  Ghép cặp uốn - cắt chịu lửa   Tách rời đơn giản                Khó tích hợp dầm phi tuyến     Tích hợp hoàn chỉnh qua MCFT

4. Luận án có cung cấp quy trình tái lập (Replication Protocol)?

Quy trình tái lập được cung cấp chi tiết toàn diện trong chương 2 và chương 3: từ biểu thức giải tích của ma trận độ cứng nhất quán $\mathbf{K}^{\text{ep,ad}}$, định nghĩa giải tích các hàm dạng bổ sung $M_1(x), M_2(x)$, các bước của thuật toán hồi quy biến dẻo (return-mapping algorithm), cho đến sơ đồ tích phân thời gian Newmark giải bài toán nhiệt.

5. Chương trình nghị sự nghiên cứu 10 năm được phác thảo ra sao?

Luận án vạch ra lộ trình 10 năm tập trung vào 3 trụ cột: (1) Xây dựng thư viện phần tử dầm - cột 3D đa tỷ lệ tích hợp phá hủy mỏi nhiệt; (2) Phát triển lý thuyết tiếp xúc nhiệt - cơ học cho bề mặt cốt thép - bê tông dưới điều kiện cháy nổ kết hợp va chấn; (3) Tự động hóa tính toán an toàn công trình chịu lửa dựa trên kỹ thuật tính toán song song hiệu năng cao (HPC).

Kết luận

Luận án tiến sĩ của Ngô Văn Minh đại diện cho bước tiến đột phá trong lĩnh vực Cơ học tính toán và Kỹ thuật Công trình, với 6 đóng góp cốt lõi:

  1. Thiết lập hoàn chỉnh khung nhiệt động lực học phân bố cho bài toán phá hủy cục bộ nhiệt - cơ có tính đến vùng xử lý nứt dẻo/phá hủy khuếch tán.
  2. Đề xuất và chứng minh giải thuật bước nhảy thông lượng nhiệt tại bề mặt phá hủy $\llbracket Q \rrbracket = \mathcal{D}_{\text{int},\delta}$.
  3. Hiện thực hóa thành công phương pháp ED-FEM làm giàu chuyển vị và biến dạng, đảm bảo tính khách quan của lưới và hiệu quả tính toán vượt trội nhờ ngưng tụ tĩnh.
  4. Phát triển giải thuật phân tách toán tử đoạn nhiệt (adiabatic split) ổn định số tuyệt đối cho bài toán phi tuyến mạnh.
  5. Mở rộng thành công mô hình hợp lực ứng suất dầm Timoshenko tích hợp phá hủy uốn và cắt đồng thời cho khung bê tông cốt thép chịu lửa thực tế.
  6. Cung cấp bộ số liệu kiểm chứng định lượng chính xác về cơ chế tiêu tán năng lượng nhiệt - cơ, tạo tiền đề vững chắc cho việc thiết kế và thẩm định an toàn các công trình xây dựng trước thảm họa hỏa hoạn.