Tổng quan về luận án

Nghiên cứu về sự lan truyền sóng và ứng xử cơ học tại các bề mặt phân chia có hình học phức tạp giữ vai trò then chốt trong cơ học vật rắn biến dạng, địa vật lý thăm dò và kỹ thuật vật liệu tiên tiến. Luận án tiến sĩ cơ học với đề tài "Thuần nhất hóa biên phân chia độ nhám cao của lý thuyết đàn hồi xốp và đàn hồi micropolar" do NCS. Nguyễn Thị Kiều thực hiện (chuyên ngành Cơ học vật rắn, mã số 9440109.02, dưới sự hướng dẫn khoa học của GS. TSKH. Phạm Chí Vĩnh và TS. Đỗ Xuân Tùng tại Trường Đại học Khoa học Tự nhiên – Đại học Quốc gia Hà Nội, 2021) là một công trình tiên phong giải quyết trọn vẹn bài toán cơ học tiếp xúc vi mô - vĩ mô tại biên có độ gồ ghề lớn.

Trong thực tế kỹ thuật, biên phân chia độ nhám cao (highly rough interface) – nơi biên độ nhám $A$ lớn hơn rất nhiều so với chu kỳ tuần hoàn $\epsilon$ ($A \gg \epsilon$, hay tỉ số $A/\epsilon \gg 1$) – xuất hiện phổ biến trong cấu trúc địa tầng nứt nẻ, liên kết vật liệu composite phân lớp, xương xốp sinh học và các siêu bề mặt nhân tạo. Phương pháp số truyền thống (FEM, BEM) khi xử lý biên nhám cao thường gặp khủng hoảng về số lượng phần tử lưới vi mô gần biên, dẫn đến nghiệm số mất ổn định và chi phí tính toán bùng nổ. Phương pháp thuần nhất hóa tiệm cận (asymptotic homogenization) khắc phục triệt để nhược điểm này bằng cách thay thế lớp không gian chứa biên nhám $-A < x_3 < 0$ bằng một lớp vật liệu tương đương có biên phẳng hoàn toàn, được chi phối bởi hệ phương trình vi phân đạo hàm riêng hiệu dụng (governing homogenized equations).

   Miền thực tế (Độ nhám cao A >> ε)              Mô hình thuần nhất hóa tương đương
   
      Ω(+) (Bán không gian trên)                      Ω(+) (Bán không gian trên)
 ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~          -------------------------------------- x3 = 0
 /\    /\    /\    /\    /\   Biên nhám               Lớp vật liệu thuần nhất hóa hiệu dụng
/  \  /  \  /  \  /  \  /  \  L: x3=h(x1/ε)           - Hệ số tường minh <·>(x3)
\  /  \  /  \  /  \  /  \  /  (-A < x3 < 0)           - Thay thế chính xác biên nhám L
 \/    \/    \/    \/    \/                      -------------------------------------- x3 = -A
 ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
      Ω(-) (Bán không gian dưới)                      Ω(-) (Bán không gian dưới)

Khoảng trống nghiên cứu then chốt (Research Gap): Trước công trình này, các nghiên cứu thuần nhất hóa biên nhám cao chủ yếu dừng lại ở các bài toán vô hướng đơn phương trình (toán tử Laplace, Poisson với điều kiện Dirichlet, Robin, Signorini như Amirat et al., 2004; De Maio et al., 2005). Đối với các hệ đàn hồi nhiều trường chuyển vị, Nevard & Keller (1997) và Gilbert & Ou (2003) chỉ thu được hệ phương trình thuần nhất hóa dưới dạng ẩn (implicit form), đòi hỏi phải giải số đồng thời 27 phương trình vi phân đạo hàm riêng trên ô cơ sở để xác định hệ số. Mặc dù Vinh & Tùng (2010–2015) đã mở đường cho phương trình dạng hiện (explicit form) trong môi trường đàn hồi cổ điển, đàn-nhiệt và đàn-điện, nhưng đối với các môi trường dị hướng phức tạp như môi trường đàn hồi xốp (poroelasticity) và môi trường đàn hồi vi cực (micropolar elasticity), bài toán dạng hiện vẫn hoàn toàn bỏ ngỏ do sự xuất hiện của các số hạng chuyển vị trực tiếp trong điều kiện liên tục ứng suất vi mô.

Luận án tập trung giải quyết 3 câu hỏi nghiên cứu (Research Questions - RQ) và kiểm chứng 3 giả thuyết khoa học (Hypotheses - H):

  • RQ1: Làm thế nào để thiết lập hệ phương trình đạo hàm riêng thuần nhất hóa dạng giải tích hiển (explicit closed-form PDEs) cho lớp biên nhám cao trong môi trường đàn hồi xốp 2D (mô hình Auriault) và 3D (mô hình Biot)?
  • RQ2: Cấu trúc ma trận nào cho phép tương thích hóa điều kiện liên tục biên chứa cả đạo hàm bậc nhất lẫn bản thân vector chuyển dịch trong lý thuyết đàn hồi Micropolar của Eringen, nhằm triệt tiêu sự cồng kềnh của toán tử vi mô?
  • RQ3: Quy luật phản xạ, khúc xạ của sóng ngang SH và sóng chuyển dịch dọc qua lớp biên nhám tương đương diễn ra như thế nào, và liệu cấu trúc nhám tuần hoàn có tạo ra các dải cấm sóng (band gaps) chặn hoàn toàn năng lượng sóng hay không?
  • H1: Tồn tại một biểu diễn ma trận mở rộng bậc 4 (chứa các khối $A_{14}, A_{34}$) đồng dạng cấu trúc với điều kiện nhảy liên tục tiếp xúc trên đường cong $L$, cho phép đưa nghiệm ô cơ sở $N^1, N^{11}, N^{13}$ về các tích phân trung bình giải tích $\langle \cdot \rangle(x_3)$.
  • H2: Hệ số phản xạ ($R$) và khúc xạ ($T$) của sóng đàn hồi qua biên nhám cao có thể tính toán chính xác bằng công thức dạng đóng (closed-form solutions) thông qua ma trận chuyển (transfer matrix) của lớp vật liệu không thuần nhất.
  • H3: Cấu trúc biên nhám hình lược (comb-like), hình răng cưa (sawtooth) và hình sin (sinusoidal) đóng vai trò như các bộ lọc siêu vật liệu cơ học (mechanical metamaterial filters), tạo ra các dải cấm sóng phụ thuộc trực tiếp vào tỉ số độ cứng $\mu_0/\mu_1$, góc tới $\theta$ và tham số hình học $A/a$.

Nghiên cứu áp dụng khung lý thuyết liên ngành tích hợp: Lý thuyết thuần nhất hóa hai thang đo của Auriault (1980, 1985), Lý thuyết cơ học môi trường xốp Biot (1956, 1962), Lý thuyết môi trường vi cực Cosserat-Eringen (1966) kết hợp với kỹ thuật ma trận vi mô - vĩ mô của Vinh & Tùng (2012, 2015). Phạm vi nghiên cứu bao quát từ biến dạng phẳng suy rộng hai chiều đến không gian ba chiều, khảo sát định lượng toàn diện trên các cấu trúc nhám tuần hoàn với các góc tới từ $0^\circ$ đến $90^\circ$, tần số không thứ nguyên $\omega H/(2\pi b_1)$ biến thiên rộng, cung cấp một bước đột phá có tính định lượng cao cho cơ học tính toán hiện đại.

Literature Review và Positioning

Lịch sử phát triển của cơ học tiếp xúc và tán xạ sóng tại biên gồ ghề được định hình qua hai nhánh phương pháp luận chính: phương pháp nhiễu loạn cho độ nhám thấp và phương pháp thuần nhất hóa tiệm cận cho độ nhám cao.

1. Nhánh nghiên cứu biên phân chia độ nhám thấp ($A \ll \epsilon$)

Khi biên độ gồ ghề nhỏ hơn rất nhiều so với chu kỳ và bước sóng, phương pháp nhiễu loạn giải tích (perturbation method) là công cụ kinh điển được áp dụng rộng rãi. Hawwa & Asfar (1996) đã nghiên cứu quá trình truyền sóng SH trong lớp có hai mặt hình sin yếu bằng cách khai triển chuỗi lũy thừa theo tham số bé đặc trưng cho độ nhấp nhô biên, đưa bài toán về chuỗi các phương trình xấp xỉ liên tiếp. Singh & Tomar (2008) tiếp tục mở rộng phương pháp nhiễu loạn bậc nhất để tính hệ số phản xạ, khúc xạ của sóng quasi-P qua mặt phân chia nhám giữa hai bán không gian đàn hồi có ứng suất trước, chứng minh sự phụ thuộc của trường tán xạ vào góc tới và gia số tính chất cơ học. Tuy nhiên, các kỹ thuật nhiễu loạn hoàn toàn suy biến và mất tính hội tụ khi tỉ số độ nhám $A/\epsilon \ge 1$, khiến chúng không thể áp dụng cho các bề mặt nhám sâu.

2. Nhánh nghiên cứu biên phân chia độ nhám cao ($A \gg \epsilon$)

Để xử lý các biên dao động mạnh, các hướng tiếp cận ban đầu dựa trên mô hình hóa hình học và phương pháp số:

  • Zaki & Neureuther (1971) giải số hệ phương trình tích phân để khảo sát sự tán xạ sóng điện từ trên bề mặt hình sin độ nhám cao, nhưng gặp trở ngại nghiêm trọng về thời gian tính toán.
  • Kohn & Vogelius (1984) đặt nền móng giải tích khi sử dụng lý thuyết đàn hồi tuyến tính 3D và khai triển tiệm cận để nghiên cứu sự uốn của bản mỏng có độ dày biến thiên nhanh (như bản có gân gia cường dày đặc), thu được phương trình vi phân bậc bốn cho chuyển dịch mặt trung hòa.
  • Talbot et al. (1990) lần đầu tiên đề xuất ý tưởng thay thế miền chứa biên phân chia nhám cao tuần hoàn tùy ý giữa hai môi trường điện từ bằng một lớp chuyển tiếp dị hướng có tính chất biến thiên liên tục theo chiều dày, giúp giảm thiểu đáng kể hệ số phản xạ sóng.
  • Trong giải tích thuần túy, Amirat et al. (2004, 2007) và De Maio et al. (2005) đã thiết lập các đánh giá sai số nghiêm ngặt giữa nghiệm xấp xỉ thuần nhất hóa và nghiệm chính xác cho toán tử Laplace và Poisson với biên răng lược mang điều kiện biên Robin hoặc Signorini (Kazimerchuk & Melnyk, 2009).
Tiến trình phát triển của lý thuyết thuần nhất hóa biên nhám cao:

[Trước 2010] Khủng hoảng dạng ẩn (Implicit Barrier)
  - Nevard & Keller (1997): Đàn hồi 3D dị hướng general -> Hệ ẩn 27 phương trình vi phân
  - Gilbert & Ou (2003): Đàn hồi xốp 3D (Auriault) -> Hệ phương trình dạng ẩn, không khả thi giải tích
        │
        ▼
[2010 - 2015] Đột phá dạng hiện cấp 1 (Vinh & Tùng)
  - Biểu diễn ma trận phương trình cơ bản + điều kiện liên tục
  - Thu được phương trình dạng hiện (Explicit) cho Đàn hồi cổ điển, Đàn - Nhiệt, Đàn - Điện
  - Hạn chế: Chưa xử lý được hệ có chuyển vị trong điều kiện liên tục (Bị cồng kềnh hoặc bế tắc)
        │
        ▼
[Luận án - 2021] Đột phá dạng hiện cấp 2 (Nguyễn Thị Kiều)
  - Phát hiện "Dạng tương đồng ma trận" (Isomorphic Matrix Formulation)
  - Giải quyết trọn vẹn: Đàn hồi xốp 2D (Auriault), Đàn hồi xốp 3D (Biot), Đàn hồi Micropolar (Eringen)
  - Thiết lập công thức dạng đóng cho hệ số R, T và phát hiện dải cấm sóng (Band Gaps)

3. Vị thế học thuật và đối sánh quốc tế của luận án

Khoảng cách lớn nhất của giai đoạn trước năm 2010 là "rào cản dạng ẩn" (implicit barrier). Khi Nevard & Keller (1997) thuần nhất hóa biên nhám cao 3D giữa hai vật thể đàn hồi dị hướng, các hệ số hiệu dụng chỉ được định nghĩa một cách trừu tượng thông qua hệ 27 phương trình đạo hàm riêng trên ô nhân tuần hoàn. Tương tự, Gilbert & Ou (2003) nghiên cứu sự truyền sóng âm qua biên nhám cao giữa hai vật liệu đàn hồi xốp theo mô hình Auriault nhưng cũng dừng lại ở hệ phương trình thuần nhất hóa dạng ẩn không thể ứng dụng trực tiếp.

Luận án của NCS. Nguyễn Thị Kiều đã định vị một bước nhảy vọt về phương pháp luận khi kế thừa kỹ thuật ma trận của Vinh & Tùng (2010, 2012) nhưng vượt qua được giới hạn cố hữu của phương pháp này. So sánh trực tiếp với công trình của Gilbert & Ou (2003) và nghiên cứu của Kumar & Sharma (2008) về tán xạ sóng vi cực, luận án lần đầu tiên cung cấp hệ phương trình thuần nhất hóa dạng hiện tường minh 100% cho cả hai môi trường phức tạp bậc nhất: Đàn hồi xốp (Poroelasticity) và Đàn hồi vi cực (Micropolar elasticity). Các hệ số thuần nhất hóa thu được là hàm giải tích trực tiếp của các hằng số đàn hồi vật liệu ($\mathbf{C}$, Darcy permeability $\mathbf{K}$, Biot coupling $\alpha$, Micropolar parameters $\mu, \kappa, \gamma$) và hàm hình học biên $\langle \cdot \rangle(x_3)$, mở ra công cụ phân tích giải tích hoàn toàn thay thế việc giải số 27 phương trình vi mô phức tạp.

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

Đóng góp cho lý thuyết (Theoretical Advancements)

Luận án mở rộng và làm phong phú hệ thống lý thuyết cơ học môi trường liên tục thông qua ba đóng góp cốt lõi:

  1. Mở rộng lý thuyết đàn hồi xốp vi mô Auriault và liên tục Biot: Cầu nối giải tích giữa cấu trúc vi mô chất lỏng nhớt Newton không nén được trong mạng lỗ rỗng với đáp ứng vĩ mô của sóng âm được thiết lập dưới dạng giải tích đóng.
  2. Nâng tầm lý thuyết đàn hồi Micropolar của Eringen: Luận án đã giải quyết bài toán biên nhám vi cực – nơi trường biến dạng được mô tả qua 6 bậc tự do (3 chuyển vị tịnh tiến $u_i$ và 3 chuyển vị vi quay $\phi_i$ của phần tử vật chất) – vốn chứa ten-sơ vi ứng suất bất đối xứng và ten-sơ ngẫu lực ứng suất (couple stress tensor).
  3. Phát hiện quy luật tương thích cấu trúc ma trận (Matrix Isomorphism): Đây là đóng góp đột phá về mặt toán - cơ. Tác giả phát hiện rằng nguyên nhân gây ra sự bùng nổ cồng kềnh trong các nghiên cứu trước đây là do sự bất đối xứng giữa phương trình cân bằng vi mô và điều kiện nhảy biên. Bằng cách tái cấu trúc phương trình trạng thái ma trận chứa thêm các toán tử chuyển vị trực tiếp, phương trình thuần nhất hóa thu được đạt độ tinh gọn tối ưu, chỉ chứa duy nhất các toán tử lấy trung bình một chiều $\langle \cdot \rangle(x_3)$.
                      MÔ HÌNH TOÁN - CƠ ĐỒNG BỘ CỦA LUẬN ÁN
                      
  [Lý thuyết Đàn hồi xốp 2D/3D]               [Lý thuyết Đàn hồi Micropolar]
   (Auriault, 1980; Biot, 1956)                     (Eringen, 1966)
   - Chuyển vị khung rắn u                          - Chuyển vị tịnh tiến u
   - Dịch chuyển tương đối lỏng w                   - Vector vi quay phi
   - Áp suất lỗ rỗng p                              - Ứng suất & ngẫu lực bất đối xứng
                │                                           │
                └─────────────────────┬─────────────────────┘
                                      │
                                      ▼
             [TÁI CẤU TRÚC PHƯƠNG TRÌNH TRẠNG THÁI DẠNG MA TRẬN]
                 (A_11 v_{,1} + A_13 v_{,3} + A_14 v)_{,1} + ...
                  Khớp hoàn hảo với điều kiện liên tục [Σ_i n_i] = 0
                                      │
                                      ▼
             [KHAI TRIỂN TIỆM CẬN ĐA THANG ĐO VI MÔ - VĨ MÔ]
                 U = V + ε(N^1 V + N^11 V_{,1} + N^13 V_{,3}) + ...
                                      │
                                      ▼
             [PHƯƠNG TRÌNH THUẦN NHẤT HÓA DẠNG HIỆN TƯỜNG MINH]
                 Σ̂_1,1 + Σ̂_3,3 + B̂ V_{,1} + D̂ V_{,3} + Ê V + F̂ = 0

Khung phân tích độc đáo (Novel Analytical Framework)

Hệ phương trình cơ bản dạng ma trận được luận án thiết lập dưới dạng tổng quát: $$(A_{11}\mathbf{v}{,1} + A{13}\mathbf{v}{,3} + A{14}\mathbf{v}){,1} + (A{31}\mathbf{v}{,1} + A{33}\mathbf{v}{,3} + A{34}\mathbf{v}){,3} + B\mathbf{v}{,1} + D\mathbf{v}_{,3} + E\mathbf{v} + F = 0$$

Trong đó vector ẩn trạng thái liên hợp được chuẩn hóa thành $\mathbf{v} = [u_1, u_2, u_3, p]^T$ (đối với đàn hồi xốp) hoặc chứa bổ sung thành phần vi quay $\mathbf{v} = [u_1, u_2, \phi_3]^T$ (đối với đàn hồi micropolar biến dạng phẳng).

Toán tử lấy trung bình giải tích dọc theo đường $x_3 = \text{const}$ (với $-A < x_3 < 0$) qua ô chu kỳ $[0, 1]$ được định nghĩa chính xác: $$\langle \Phi \rangle(x_3) = \int_0^1 \Phi(y, x_3),dy = (y_2 - y_1)\Phi^{(1)} + (1 - y_2 + y_1)\Phi^{(2)}$$ với $y_1, y_2$ ($0 \le y_1 < y_2 \le 1$) là hai nghiệm của phương trình biên hình học $h(y) = x_3$.

Nhờ khung phân tích này, nghiệm của ba bài toán ô cơ sở địa phương (local cell problems) đối với các ma trận hiệu chỉnh bậc một $N^1, N^{11}, N^{13}$ được xác định bằng các biểu thức giải tích đóng:

  • $N^1_{,y} = A_{11}^{-1} \langle A_{11}^{-1} \rangle^{-1} \langle A_{11}^{-1} A_{14} \rangle - A_{11}^{-1} A_{14}$
  • $I + N^{11}{,y} = A{11}^{-1} \langle A_{11}^{-1} \rangle^{-1}$
  • $N^{13}{,y} = A{11}^{-1} \langle A_{11}^{-1} \rangle^{-1} \langle A_{11}^{-1} A_{13} \rangle - A_{11}^{-1} A_{13}$

Kết quả dẫn đến hệ phương trình chuyển động thuần nhất hóa tương đương của lớp vật liệu trung gian: $$\hat{\Sigma}{1,1} + \hat{\Sigma}{3,3} + \langle B A_{11}^{-1} \rangle \langle A_{11}^{-1} \rangle^{-1} \mathbf{V}{,1} + \Big[ \langle D \rangle + \langle B A{11}^{-1} \rangle \langle A_{11}^{-1} \rangle^{-1} \langle A_{11}^{-1} A_{13} \rangle - \langle B A_{11}^{-1} A_{13} \rangle \Big] \mathbf{V}_{,3} + \dots = 0$$ Điều này chứng minh tính chuẩn xác tuyệt đối về mặt toán học và bảo toàn định luật biến suy năng lượng cơ học.

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

Thiết kế nghiên cứu (Research Design)

Luận án tuân thủ thế giới quan thực chứng suy diễn toán học (Positivist / Mathematical Deductive Framework), kết hợp giải tích tiệm cận vi mô - vĩ mô (two-scale asymptotic expansion) với phương pháp tính số trên nền tảng đại số ma trận Stroh. Thiết kế đa quy mô được chia thành:

  • Quy mô vi mô (Micro-scale $y = x_1/\epsilon$): Nơi biên độ nhám $h(y)$ biến thiên tuần hoàn nhanh trong ô cơ sở $Y = [0, 1]$.
  • Quy mô vĩ mô (Macro-scale $x_1, x_3$): Nơi trường sóng trung bình truyền qua môi trường tương đương dưới dạng các mode sóng điều hòa theo thời gian $\exp(i\omega t)$.

Quy trình nghiên cứu chuẩn mực (Rigorous Research Protocol)

  1. Thiết lập phương trình trạng thái ma trận: Chuyển đổi toàn bộ hệ phương trình vi phân cân bằng động lực học, định luật thấm Darcy suy rộng và hệ phương trình cấu tạo vật liệu dị hướng sang dạng ma trận khối cấp $4 \times 4$ hoặc $6 \times 6$.
  2. Đồng nhất hóa tiệm cận các bậc $\epsilon$: Khai triển trường chuyển vị $U(x_1, y, x_3) = \mathbf{V} + \epsilon(N^1\mathbf{V} + N^{11}\mathbf{V}{,1} + N^{13}\mathbf{V}{,3}) + \epsilon^2(\dots) + O(\epsilon^3)$. Thay thế vào phương trình cơ bản và điều kiện biên tiếp xúc dạng nhảy $[v]_L = 0, [\Sigma_1 n_1 + \Sigma_3 n_3]_L = 0$, sau đó đồng nhất hóa triệt tiêu các hệ số của $\epsilon^{-1}, \epsilon^0, \epsilon^1$.
  3. Phân tích giải tích ô cơ sở: Giải chính xác hệ phương trình vi phân thường đối với biến vi mô $y$ trên khoảng $[0, 1]$, bảo đảm điều kiện tuần hoàn $N(0) = N(1)$.
  4. Kỹ thuật chia lớp liên tục (Stratified Layer Discretization): Đối với lớp vật liệu không thuần nhất $-A < x_3 < 0$, chia thành $N$ lớp con thuần nhất mỏng ($N$ từ 50 đến 200 lớp) để xây dựng ma trận chuyển truyền sóng (Propagator Transfer Matrix), bảo đảm tính ổn định nghiệm số tuyệt đối so với các phương pháp phần tử hữu hạn thông thường.
       QUY TRÌNH TÍNH TOÁN SỐ BẰNG MA TRẬN CHUYỂN TẦNG (TRANSFER MATRIX)
       
  [Lớp biên phẳng x3 = 0]    v(+) = [u1, u2, u3, p]^T,  Σ3(+)
  ==============================================================
    Lớp con 1 (Thuần nhất)   Matran T1 = exp(A(x3_1) * Δz)
  --------------------------------------------------------------
    Lớp con 2 (Thuần nhất)   Matran T2 = exp(A(x3_2) * Δz)
  --------------------------------------------------------------
    ...                      ...
  --------------------------------------------------------------
    Lớp con N (Thuần nhất)   Matran TN = exp(A(x3_N) * Δz)
  ==============================================================
  [Lớp biên phẳng x3 = -A]   v(-) = [u1, u2, u3, p]^T,  Σ3(-)
  
  ==> Ma trận chuyển toàn phần: M = TN * TN-1 * ... * T2 * T1
  ==> Thiết lập hệ phương trình đại số tuyến tính tìm trực tiếp R và T!

Dữ liệu và phân tích số (Data & Computation)

Toàn bộ thuật toán được lập trình giải tích và mô phỏng số bằng phần mềm khoa học chuyên dụng (MATLAB/Mathematica với sự hỗ trợ lập trình từ TS. Trần Thanh Tuấn). Độ chính xác của nghiệm số được kiểm chứng chéo (robustness check) qua việc so sánh các trường hợp giới hạn: khi độ cứng hai môi trường đồng nhất ($\mu^{(1)} = \mu^{(2)}$), hệ số phản xạ triệt tiêu $|R| \to 0$ và hệ số truyền qua đạt cực đại $|T| \to 1$, khớp chính xác 100% với lý thuyết truyền sóng cổ điển.

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

Những phát hiện then chốt (Key Breakthrough Findings)

1. Thiết lập thành công phương trình dạng hiện cho môi trường xốp Auriault (2D) và Biot (3D)

Luận án đã thu được trọn vẹn biểu thức giải tích của các toán tử ma trận thuần nhất hóa. Đối với vật liệu đàn hồi xốp trực hướng, các phương trình dạng thành phần được rút gọn hoàn toàn về các tích phân hình học của hàm $h(y)$, xóa bỏ hoàn toàn nhu cầu giải số 27 phương trình vi mô như các công bố quốc tế trước đó (Nevard & Keller, 1997; Gilbert & Ou, 2003).

2. Giải mã hiện tượng Dải cấm sóng (Band Gaps) do biên nhám cao tuần hoàn

Khi khảo sát sự phản xạ và khúc xạ của sóng ngang SH và sóng chuyển dịch dọc, luận án phát hiện: Biên phân chia độ nhám cao tuần hoàn hoạt động như một cấu trúc tinh thể âm học (phononic crystal), có khả năng ngăn chặn hoàn toàn sự truyền qua của sóng trong các dải tần số xác định ($|T| = 0, |R| = 1$).

Đồ thị biểu diễn cơ chế dải cấm sóng (Band Gap) theo tần số:
Hệ số truyền qua |T|
 1.0 ┤──────\                  /──────\                  /──────
     │       \                /        \                /       
 0.5 ┤        \              /          \              /        
     │         \            /            \            /         
 0.0 ┼──────────\__________/──────────────\__________/──────────
     0         ω_L1      ω_R1            ω_L2      ω_R2         Tần số ωH/(2π b1)
                 [ DẢI CẤM SÓNG 1 ]        [ DẢI CẤM SÓNG 2 ]
                 (|T| = 0, |R| = 1)        (|T| = 0, |R| = 1)
  • Ảnh hưởng của tỉ số mô đun trượt $\mu_0/\mu_1$: Khi độ tương phản cơ học giữa hai bán không gian tăng từ $\mu_0 = 2\mu_1$ lên $\mu_0 = 20\mu_1$, độ rộng dải cấm sóng mở rộng đáng kể; các dải cấm xuất hiện dày đặc hơn ở miền tần số thấp.
  • Ảnh hưởng của hình học biên nhám:
    • Biên dạng hình lược (comb-like) tạo ra các dải cấm sâu và sắc nét nhất.
    • Biên dạng răng cưa (sawtooth) và hình sin (sinusoidal) làm dịch chuyển biên dải cấm về vùng tần số cao hơn và làm mịn các đỉnh cộng hưởng tán xạ.
  • Ảnh hưởng của góc tới $\theta$: Khi góc tới $\theta$ tiến gần đến góc tới hạn $\theta_c$, xuất hiện hiện tượng phản xạ toàn phần và biến đổi pha phức tạp giữa sóng phản xạ và sóng khúc xạ.

3. Bảng tổng hợp các đặc tính sóng và dải cấm sóng qua các dạng biên nhám

Dưới đây là bảng tổng hợp định lượng các thông số cơ học và đặc tính tán xạ sóng qua ba cấu trúc biên nhám tiêu biểu được khảo sát trong luận án:

Thông số & Đặc tính tán xạ Biên dạng Hình lược (Comb-like) Biên dạng Răng cưa (Sawtooth) Biên dạng Hình sin (Sinusoidal)
Hàm mô tả biên $h(y)$ Phân đoạn hằng số ($y_1, y_2 = \text{const}$) Bậc nhất từng đoạn ($h(y) = ky + b$) Hàm lượng giác $h(y) = -\frac{A}{2}(1 + \cos 2\pi y)$
Độ phức tạp hệ số $\langle \cdot \rangle(x_3)$ Hằng số độc lập với $x_3$ (rất gọn) Hàm đại số tuyến tính theo $x_3$ Hàm phi tuyến chứa $\arccos(x_3)$
Độ rộng dải cấm sóng ($\Delta \omega$) Rộng nhất (Khả năng chặn sóng tối đa) Trung bình (Bị dịch về tần số cao) Hẹp và phụ thuộc mạnh vào biên độ $A$
**Hệ số phản xạ cực đại $ R _{\max}$** Đạt xấp xỉ $1.0$ trên toàn dải cấm
Độ nhạy với góc tới $\theta$ Duy trì dải cấm ổn định ($0^\circ \le \theta \le 60^\circ$) Dải cấm biến dạng nhanh khi $\theta > 45^\circ$ Biến thiên liên tục theo hàm trơn

Implications đa chiều

                       CÁC ỨNG DỤNG ĐA NGÀNH CỦA LUẬN ÁN
                       
 [Địa chấn & Địa vật lý]      [Kỹ thuật Siêu âm NDT]       [Siêu vật liệu Âm học]
  - Giải đoán sóng địa chấn    - Kiểm tra khuyết tật sâu    - Thiết kế tấm chắn sóng
  - Thăm dò bồn chứa dầu khí   - Đánh giá liên kết hàn nhám - Bộ lọc tần số cơ học
  - Mô hình hóa đứt gãy xốp    - Đo đặc trưng lớp composite - Triệt tiêu rung chấn
  • Ý nghĩa học thuật: Cung cấp chuẩn mực toán học mới để giải quyết các bài toán biên dao động nhanh trong cơ học lượng tử, âm học và điện từ học cấu trúc nano.
  • Thăm dò địa vật lý & Dầu khí: Cho phép giải đoán chính xác dữ liệu địa chấn tại các tầng chứa dầu khí có cấu trúc đá trầm tích xốp nứt nẻ bão hòa chất lỏng, nơi biên tiếp xúc địa tầng luôn có độ gồ ghề phức tạp.
  • Kiểm tra không phá hủy (NDT): Cung cấp thuật toán giải ngược (inverse problem) giúp xác định chiều sâu và hình dạng nhám của vết nứt ngầm thông qua phổ tần số phản xạ của sóng siêu âm.
  • Kỹ thuật vật liệu & Chống chấn: Đặt cơ sở lý thuyết cho việc thiết kế các rào cản địa chấn nhân tạo (metamaterial seismic barriers) bằng cách tạo rãnh nhám tuần hoàn dưới lòng đất để ngăn chặn sóng động đất phá hủy công trình.

Limitations và Future Research

Mặc dù đạt được những kết quả đột phá, luận án duy trì góc nhìn khoa học khách quan khi chỉ rõ 4 giới hạn nghiên cứu (Boundary Conditions):

  1. Giả thiết tính tuần hoàn lý tưởng: Biên nhám được giả thiết tuần hoàn hoàn hảo với chu kỳ $\epsilon$. Trong thực tế địa chất, bề mặt nhám thường mang tính ngẫu nhiên (random roughness) hoặc cấu trúc tự đồng dạng phân dạng (fractal surfaces).
  2. Miền biến dạng đàn hồi tuyến tính: Các phương trình cấu tạo giả định biến dạng nhỏ và chuyển động dao động điều hòa theo thời gian $\exp(i\omega t)$; chưa xét đến hiệu ứng phi tuyến vật liệu hoặc phi tuyến hình học khi biên độ sóng quá lớn.
  3. Mô hình chất lỏng lỗ rỗng đơn pha: Nghiên cứu áp dụng cho chất lỏng Newton không nén được bão hòa hoàn toàn; chưa tính đến trường hợp bão hòa từng phần (hỗn hợp khí - dầu - nước đa pha).
  4. Giới hạn tiếp xúc hoàn hảo: Giả thiết hai miền dính kết chặt tại biên vi mô ($[u]_L = 0$); chưa xét đến sự trượt ma sát (slip boundary) hoặc bong tách tiếp xúc vi mô (debonding).

Chương trình nghiên cứu 10 năm tiếp theo (Future Research Agenda):

  • Mở rộng 1: Phát triển lý thuyết thuần nhất hóa ngẫu nhiên (stochastic homogenization) cho biên nhám phân dạng (fractal interfaces).
  • Mở rộng 2: Khảo sát môi trường xốp phi tuyến chịu tải trọng va đập động lực học cường độ cao (blast loading).
  • Mở rộng 3: Tích hợp dòng chảy đa pha phi-Newton trong môi trường xốp đàn hồi nhiệt (thermo-poroelasticity).
  • Mở rộng 4: Thực nghiệm kiểm chứng trên hệ thống đo siêu âm laser đa kênh tại các mẫu vật liệu composite in 3D bề mặt nhám vi mô.

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

  • Tác động học thuật (Academic Impact): Công trình là bước tiến vượt bậc của trường phái Cơ học tính toán Việt Nam (nhóm nghiên cứu GS. TSKH. Phạm Chí Vĩnh), được kỳ vọng tạo ra chỉ số trích dẫn cao trong các tạp chí cơ học quốc tế hàng đầu ISI/Scopus (như Journal of Mechanics and Physics of Solids, Wave Motion, International Journal of Solids and Structures).
  • Chuyển đổi công nghiệp & Công nghệ vật liệu: Mở ra phương pháp chế tạo các lớp phủ triệt tiêu sóng âm (anechoic coatings) cho tàu ngầm và siêu bề mặt tiêu âm trong kỹ thuật hàng không vũ trụ.
  • Giá trị kinh tế - xã hội: Cải thiện độ chính xác trong định vị túi dầu/khí ngầm tại thềm lục địa Việt Nam, giảm thiểu rủi ro khoan thăm dò tốn kém hàng triệu USD; đồng thời đóng góp giải pháp bảo vệ công trình hạt nhân và đập thủy điện trước hiểm họa động đất.

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

                              BẢN ĐỒ ĐỐI TƯỢNG HƯỞNG LỢI
                              
  [Nghiên cứu sinh & Nhà toán - cơ] ──► Tiếp nhận khung giải tích ma trận thuần nhất hóa
  [Chuyên gia Địa vật lý & NDT]      ──► Sử dụng công thức đóng R, T để giải bài toán ngược
  [Kỹ sư R&D Siêu vật liệu]          ──► Ứng dụng dải cấm sóng (Band Gaps) để chặn dao động
  [Cơ quan Quản lý An toàn Động đất] ──► Căn cứ khoa học thiết kế hào chắn sóng địa chấn
  • Nghiên cứu sinh & Nhà khoa học chuyên sâu: Sở hữu một công cụ toán học mẫu mực để giải mã các phương trình vi phân nhiều trường trên miền biên phức tạp.
  • Chuyên gia Địa chấn & Dầu khí: Tiếp cận thuật toán tính toán nhanh hệ số sóng mà không cần xây dựng lưới vi mô nặng nề.
  • Kỹ sư R&D Vật liệu & Thiết kế âm học: Nắm bắt quy luật hình học biên để chủ động điều khiển dải cấm sóng theo yêu cầu kỹ thuật.

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ì và đã mở rộng lý thuyết nào?

Trả lời: Đóng góp độc đáo nhất là việc thiết lập hệ phương trình thuần nhất hóa dạng hiện (tường minh 100%) cho môi trường đàn hồi xốp (theo cả mô hình vi mô Auriault và mô hình liên tục Biot) và đàn hồi vi cực Micropolar. Công trình mở rộng trực tiếp lý thuyết thuần nhất hóa tiệm cận của Bensoussan et al. (1978), phá vỡ bế tắc "dạng ẩn" tồn tại suốt nhiều thập kỷ của Nevard & Keller (1997) và Gilbert & Ou (2003), đồng thời mở rộng khung phương pháp ma trận của Vinh & Tùng (2012, 2015) sang các hệ môi trường liên tục phi cổ điển.

2. Đột phá phương pháp luận của luận án khi so sánh với các công trình quốc tế?

Trả lời: So với cách tiếp cận dạng thành phần truyền thống (dẫn đến 27 phương trình vi mô ẩn như Nevard & Keller), luận án phát hiện ra "Dạng tương đồng cấu trúc ma trận" (Matrix Isomorphism). Bằng cách bổ sung các khối ma trận liên hợp chuyển vị ($A_{14}, A_{34}$) vào phương trình cơ bản, điều kiện liên tục vi mô trên biên nhám cao $L$ khớp hoàn toàn với cấu trúc đạo hàm không gian, biến các bài toán ô cơ sở trên chu kỳ $[0, 1]$ thành các tích phân giải tích một chiều đơn giản $\langle \cdot \rangle(x_3)$, triệt tiêu toàn bộ sự cồng kềnh đại số.

3. Phát hiện bất ngờ nhất từ dữ liệu số là gì?

Trả lời: Phát hiện bất ngờ nhất là sự hình thành dải cấm sóng tuyệt đối (Band Gaps) của sóng đàn hồi chỉ nhờ vào sự nhám hóa tuần hoàn của biên tiếp xúc, ngay cả khi cả hai bán không gian đều là các vật liệu đàn hồi thông thường không có cấu trúc tuần hoàn khối. Đặc biệt, với biên nhám hình lược, khi tỉ số độ cứng $\mu_0/\mu_1 = 20$, dải cấm sóng mở rộng đột biến và triệt tiêu hoàn toàn năng lượng sóng khúc xạ ($|T| = 0, |R| = 1$) trên nhiều khoảng tần số rộng.

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

Trả lời: Có. Luận án trình bày chi tiết từng bước: từ thiết lập hệ ma trận khối $A_{mn}, B, D, E, F$, giải thuật ba bài toán nhân $N^1, N^{11}, N^{13}$, định nghĩa toán tử tích phân trung bình hình học $\langle \cdot \rangle(x_3)$, đến thuật toán chia $N$ lớp con thuần nhất mỏng để xây dựng ma trận truyền sóng (Transfer Matrix) tính toán hệ số $R$ và $T$. Toàn bộ hệ phương trình đã được tường minh hóa thành phần trong các phụ lục và chương chuyên sâu.

5. Khung chương trình nghiên cứu 10 năm tới được phác thảo như thế nào?

Trả lời: Định hướng mở rộng bao gồm: (1) Chuyển đổi từ biên tuần hoàn sang bề mặt nhám ngẫu nhiên/phân dạng; (2) Tích hợp hiệu ứng tiếp xúc ma sát/trượt phi tuyến vi mô; (3) Khảo sát trường sóng bão hòa đa pha trong công nghệ thu hồi dầu tăng cường (EOR); (4) Thiết kế và thực nghiệm siêu vật liệu chặn sóng địa chấn trên mô hình bàn rung địa chấn thực tế.

Kết luận

Luận án tiến sĩ của NCS. Nguyễn Thị Kiều là một công trình khoa học xuất sắc, mẫu mực về sự kết hợp giữa tư duy toán học giải tích chặt chẽ và ứng dụng cơ học kỹ thuật hiện đại, đóng góp 5 giá trị nền tảng:

  1. Thiết lập thành công hệ phương trình thuần nhất hóa dạng hiện dạng ma trận và dạng thành phần cho biên phân chia độ nhám cao trong môi trường đàn hồi xốp 2D (mô hình Auriault) và 3D (mô hình Biot).
  2. Xây dựng hệ phương trình thuần nhất hóa dạng hiện cho lý thuyết đàn hồi Micropolar vi cực, mở rộng năng lực mô tả biến dạng vi xoay của vật liệu hạt và cấu trúc sinh học.
  3. Phát hiện quy luật tương đồng cấu trúc ma trận giữa phương trình vi phân và điều kiện tiếp xúc, giải quyết triệt để rào cản cồng kềnh đại số tồn tại nhiều thập kỷ.
  4. Tìm ra công thức dạng đóng tính toán chính xác hệ số phản xạ ($R$) và khúc xạ ($T$) của sóng ngang SH và sóng chuyển dịch dọc qua các cấu trúc biên nhám hình lược, răng cưa và hình sin.
  5. Chứng minh và định lượng hóa hiện tượng dải cấm sóng (Band Gaps), mở ra hướng ứng dụng mới trong chế tạo siêu vật liệu cơ học và rào chắn triệt tiêu sóng địa chấn.
  6. Mở ra 3 nhánh nghiên cứu liên ngành đột phá: Âm học siêu bề mặt nhám vi cấu trúc, thăm dò địa chấn tầng chứa nứt nẻ nâng cao và công nghệ kiểm tra không phá hủy siêu âm độ phân giải cao.

Công trình khẳng định bước tiến vượt bậc của Cơ học vật rắn Việt Nam, hội nhập và giải quyết trực tiếp những bài toán cốt lõi của cơ học môi trường liên tục thế giới.