Tổng quan về luận án

Sóng mặt Rayleigh, được phát hiện lần đầu bởi Lord Rayleigh vào năm 1885 đối với bán không gian đàn hồi đẳng hướng, là một trong những nền tảng lý thuyết quan trọng nhất của cơ học vật thể rắn hiện đại. Sóng Rayleigh đóng vai trò cốt lõi trong nhiều lĩnh vực khoa học và công nghệ kỹ thuật cao như địa chấn học, địa vật lý, âm học vật lý, công nghệ truyền thông vi dải, thiết bị sóng âm bề mặt (SAW - Surface Acoustic Wave) trong vi cơ điện tử (MEMS) và đánh giá không phá hủy (NDT - Non-Destructive Testing) vật liệu tiên tiến.

Trong hơn 120 năm qua, tuyệt đại đa số các công trình nghiên cứu về sóng mặt Rayleigh đều dựa trên giả thiết bán không gian có bề mặt hoàn toàn tự do đối với ứng suất ($\boldsymbol{\sigma} \cdot \mathbf{n} = \mathbf{0}$). Tuy nhiên, trong thực tiễn địa vật lý và kỹ thuật công nghệ cao, môi trường đàn hồi hầu như luôn chịu các tác động cơ - lý phức tạp ở bề mặt và trong lòng khối vật liệu. Các hiện tượng này bao gồm tương tác tiếp xúc âm học với môi trường chất lỏng hoặc lớp phủ ngoài (dẫn đến điều kiện biên trở kháng), trường ứng suất dư/biến dạng trước do kiến tạo địa chất hoặc quá trình chế tạo cơ khí, hiệu ứng quay quán tính trong các con quay vi cơ học (gyroscopes) và kết cấu composite phủ lớp mỏng chức năng. Nghiên cứu sóng Rayleigh trong các bán không gian không tự do ứng suất là một bài toán cơ học phi tuyến và dị hướng đặc biệt phức tạp, tồn tại nhiều khoảng trống tri thức lý thuyết chưa được giải quyết.

Luận án tiến sĩ cơ học của tác giả Trịnh Thị Thanh Huệ với tiêu đề "Sóng Rayleigh trong các bán không gian đàn hồi không tự do đối với ứng suất", chuyên ngành Cơ học vật thể rắn (Mã số: 62 44 01 07), dưới sự hướng dẫn khoa học của GS. Phạm Chí Vĩnh và GS. Nguyễn Đông Anh tại Trường Đại học Khoa học Tự nhiên - Đại học Quốc gia Hà Nội (2017), đã giải quyết căn bản và toàn diện bài toán lý thuyết phức tạp này.

                         ┌─────────────────────────────────────────────────────────────┐
                         │   SÓNG RAYLEIGH TRONG BÁN KHÔNG GIAN KHÔNG TỰ DO ỨNG SUẤT   │
                         └──────────────────────────────┬──────────────────────────────┘
                                                        │
          ┌─────────────────────────────────────────────┼─────────────────────────────────────────────┐
          │                                             │                                             │
          ▼                                             ▼                                             ▼
┌──────────────────────────────────┐          ┌──────────────────────────────────┐          ┌──────────────────────────────────┐
│ 1. ĐIỀU KIỆN BIÊN TRỞ KHÁNG      │          │ 2. CÓ ỨNG SUẤT TRƯỚC VÀ QUAY     │          │ 3. BÁN KHÔNG GIAN PHỦ LỚP MỎNG   │
│ - Dị hướng trực hướng            │          │ - Kéo / nén thuần túy            │          │ - Nền Monoclinic $x_3=0$         │
│ - Dị hướng Monoclinic $x_3=0$    │          │ - Đồng thời kéo / nén và cắt     │          │ - Lớp mỏng Monoclinic $x_3=0$    │
│ - Nén được và không nén được     │          │ - Hiệu ứng quay Coriolis         │          │ - Xấp xỉ hiệu dụng bậc 2         │
└──────────────────────────────────┘          └──────────────────────────────────┘          └──────────────────────────────────┘

Khoảng trống nghiên cứu (Research Gaps) cụ thể trong y văn quốc tế được xác định như sau:

  1. Khoảng trống về phương pháp giải tích: Phương pháp vectơ phân cực (Polarization Vector Method) truyền thống của Taziev (1989), Ting (1996) và Collet & Destrade (2004) chỉ áp dụng được khi ma trận Stroh là ma trận thực. Khi bài toán chịu điều kiện biên trở kháng hoặc phủ lớp mỏng, ma trận Stroh trở thành ma trận phức, khiến toàn bộ các công cụ giải tích hiện hành mất hiệu lực.
  2. Khoảng trống về điều kiện biên tổng quát: Nghiên cứu của Malischewsky (1987) và Godoy et al. (2012) chỉ khảo sát vật liệu đẳng hướng, trong đó Godoy et al. phải giản lược bài toán bằng cách gán ứng suất pháp $\sigma_{22} = 0$. Chưa có bất kỳ nghiên cứu nào thiết lập được phương trình tán sắc giải tích dạng tường minh cho vật liệu dị hướng tổng quát (trực hướng, monoclinic) nén được và không nén được chịu điều kiện biên trở kháng đầy đủ hai thành phần.
  3. Khoảng trống về trường tương tác đa vật lý: Chưa có công trình nào tích hợp đồng thời hiệu ứng ứng suất trước (vừa kéo/nén vừa cắt), hiệu ứng quay quán tính và cấu trúc lớp phủ mỏng dị hướng monoclinic trong môi trường chịu điều kiện biên trở kháng.

Nghiên cứu tập trung giải quyết 5 câu hỏi nghiên cứu cốt lõi (Research Questions - RQ) và kiểm chứng 5 giả thuyết khoa học tương ứng (Hypotheses - H):

  • RQ1: Làm thế nào để mở rộng phương pháp vectơ phân cực cho ma trận Stroh phức nhằm thu được phương trình tán sắc dạng tường minh tuyệt đối?
    H1: Tồn tại một hệ thức cơ bản Hermite suy rộng cho nghiệm tắt dần ở vô cùng, cho phép loại bỏ hoàn toàn các hằng số tích phân mà không cần giải đa thức đặc trưng bậc cao.
  • RQ2: Dạng giải tích tường minh chính xác của phương trình tán sắc sóng Rayleigh trong bán không gian trực hướng và monoclinic nén được/không nén được chịu điều kiện biên trở kháng là gì?
    H2: Phương trình tán sắc có thể biểu diễn dưới dạng đại số đóng thông qua các bất biến vật liệu và tham số trở kháng không thứ nguyên.
  • RQ3: Ảnh hưởng của trạng thái biến dạng trước (kéo-nén kết hợp trượt cắt) đến vận tốc pha sóng Rayleigh dưới điều kiện biên trở kháng thay đổi như thế nào?
    H3: Ứng suất cắt tạo ra tính bất đối xứng trong phổ vận tốc pha đối với hướng truyền sóng.
  • RQ4: Sự quay của hệ quy chiếu kết hợp với gia cố cốt sợi không nén được thay đổi đặc tính tán sắc sóng Rayleigh ra sao?
    H4: Lực quán tính Coriolis làm xuất hiện tính tán sắc nội tại ngay cả trong môi trường bán vô hạn đồng nhất.
  • RQ5: Cấu trúc tán sắc của sóng Rayleigh trong bán không gian monoclinic phủ lớp mỏng monoclinic được biểu diễn chính xác đến bậc hai của độ dày lớp mỏng như thế nào?
    H5: Điều kiện biên hiệu dụng bậc hai cho phép thay thế hoàn toàn ảnh hưởng động lực của lớp mỏng mà không làm mất tính đóng của phương trình giải tích.

Khung lý thuyết nền tảng (Theoretical Framework) của công trình kết hợp Lý thuyết đàn hồi tuyến tính dị hướng cổ điển, Cơ học môi trường liên tục phi tuyến có biến dạng trước của Biot (1965) và Ogden (1984), Hình thức luận Stroh (Stroh Formalism 1958, 1962), Lý thuyết điều kiện biên hiệu dụng (Effective Boundary Conditions) và Phương pháp vectơ phân cực phức mở rộng.

Tác động định lượng đột phá: Công trình thiết lập 100% các phương trình tán sắc dưới dạng tường minh giải tích (explicit secular equations), quy dẫn bài toán tìm nghiệm từ các phương trình siêu việt nhiều ẩn phức tạp về các phương trình đa thức đại số bậc 6 và bậc 8 không thứ nguyên, giải quyết triệt để sự phụ thuộc vào các thuật toán dò nghiệm số vốn dễ phân kỳ tại các miền tần số cao. Toàn bộ các kết quả lý thuyết đã được công bố trong 2 công trình khoa học trên các tạp chí quốc tế uy tín thuộc danh mục ISI/SCI (Applied Mathematical Modelling, Wave Motion) và 2 báo cáo toàn văn tại Hội nghị Cơ học Vật rắn Biến dạng Toàn quốc.

Literature Review và Positioning

Lịch sử phát triển lý thuyết sóng mặt Rayleigh trải qua hơn một thế kỷ với 5 trường phái nghiên cứu chính:

┌────────────────────────────────────────────────────────────────────────────────────────┐
│                        TIẾN TRÌNH PHÁT TRIỂN LÝ THUYẾT SÓNG RAYLEIGH                   │
├────────────────────────────────────────────────────────────────────────────────────────┤
│ 1885: Lord Rayleigh ─────────► Phát hiện sóng mặt trong môi trường đẳng hướng tự do    │
│ 1958: Stroh ─────────────────► Xác lập hình thức luận ma trận Stroh 6D cho dị hướng     │
│ 1976: Chadwick ──────────────► Phương trình tán sắc cho môi trường trực hướng tự do    │
│ 1987: Malischewsky ──────────► Sóng Rayleigh môi trường đẳng hướng có trở kháng biên   │
│ 1989: Taziev ────────────────► Sáng lập phương pháp vectơ phân cực (ma trận thực)     │
│ 2001-2004: Destrade ─────────► Tích phân đầu & vectơ phân cực cho monoclinic/quay      │
│ 2004: Vinh & Ogden ──────────► Công thức nghiệm giải tích Cardan cho trực hướng        │
│ 2012: Godoy, Durán, Nkemzi ──► Khảo sát tồn tại nghiệm đẳng hướng (giản lược sigma22=0)│
│ 2017: TRỊNH THỊ THANH HUỆ ──► PHƯƠNG PHÁP VECTƠ PHÂN CỰC PHỨC & BÁN KHÔNG GIAN KHÔNG   │
│                               TỰ DO ỨNG SUẤT DỊ HƯỚNG TỔNG QUÁT (LUẬN ÁN NÀY)          │
└────────────────────────────────────────────────────────────────────────────────────────┘
  1. Trường phái sóng Rayleigh tự do ứng suất cổ điển: Khởi đầu từ Rayleigh (1885), Chadwick (1976) tìm ra phương trình tán sắc cho vật liệu trực hướng. Rahman & Barber (1995) thiết lập công thức nghiệm chính xác bậc ba đầu tiên cho môi trường đẳng hướng. Nkemzi (1997) sử dụng lý thuyết bài toán biên Riemann nhưng nghiệm thu được quá phức tạp và tồn tại sai số toán học (theo Malischewsky, 2004). Đến năm 2004, Vinh & Ogden đã chứng minh chặt chẽ công thức Cardan của Malischewsky và tìm ra nghiệm tường minh cho vật liệu trực hướng nén được và không nén được (Vinh & Ogden, 2004a, 2004b).
  2. Trường phái điều kiện biên trở kháng (Impedance Boundary Conditions): Malischewsky (1987) là người tiên phong thiết lập phương trình tán sắc cho bán không gian đẳng hướng chịu điều kiện biên trở kháng dạng $\sigma_{12} + \omega Z_1 u_1 = 0, \sigma_{22} + \omega Z_2 u_2 = 0$. Tuy nhiên, công trình của Malischewsky chưa chứng minh được tính tồn tại và duy nhất của nghiệm. Mãi đến năm 2012, Godoy, Durán và Nkemzi mới phân tích tính duy nhất nhưng chỉ giải quyết được trường hợp đặc biệt khi ứng suất pháp triệt tiêu ($\sigma_{22} = 0$).
  3. Trường phái hình thức luận Stroh và Vectơ phân cực: Stroh (1958, 1962) đưa ra biểu diễn ma trận vi phân cấp 6 mô tả động lực học đàn hồi dị hướng. Taziev (1989) đề xuất phương pháp vectơ phân cực để tìm phương trình tán sắc dạng hiện. Lý thuyết này được Ting (1996, 2002), Destrade (2001, 2003, 2004), Collet & Destrade (2004) hoàn thiện xuất sắc cho các môi trường dị hướng monoclinic và môi trường quay. Tuy nhiên, toàn bộ các nghiên cứu này đều dựa trên giả định bất khả biến: ma trận Stroh bắt buộc phải là ma trận thực.
  4. Trường phái môi trường quay và biến dạng trước: Schoenberg & Censor (1973), Clarke & Burdess (1994), Destrade (2004) chỉ ra rằng chuyển động quay của hệ quy chiếu phá vỡ tính không tán sắc của sóng Rayleigh. Cùng thời gian đó, Biot (1965), Ogden (1984), Vinh & Giang (2012) khảo sát sóng truyền trong vật thể chịu ứng suất trước, chứng minh trường ứng suất dư làm biến đổi vận tốc lan truyền sóng mặt.
  5. Trường phái bán không gian phủ lớp mỏng: Achenbach & Keshava (1967) đưa ra xấp xỉ bậc 4 nhưng phụ thuộc vào hằng số chưa xác định. Bovik (1996), Rokhlin & Huang (1992), Vinh & Linh (2012, 2013) phát triển phương pháp điều kiện biên hiệu dụng khai triển chuỗi Taylor để thay thế lớp mỏng bằng một liên hệ biên tương đương trên bề mặt nền.

Tranh luận khoa học và Định vị nghiên cứu

Y văn tồn tại hai cuộc tranh luận lý thuyết lớn:

  • Tranh luận 1: Phương pháp giải số ma trận chuyển (Transfer Matrix) so với Phương pháp giải tích tường minh. Các nhà nghiên cứu tính toán số lập luận rằng phương pháp số có thể giải được mọi hình học lớp phức tạp; ngược lại, trường phái giải tích chứng minh rằng phương pháp số gặp hiện tượng mất ổn định nghiệm (numerical instability/loss of precision) ở tần số cao khi bề dày lớp lớn hơn bước sóng. Luận án định vị vững chắc theo trường phái giải tích tường minh chính xác.
  • Tranh luận 2: Tính khả thi của Phương pháp vectơ phân cực trong bài toán tiêu tán/trở kháng biên. Nhiều nhà cơ học quốc tế cho rằng khi điều kiện biên trở kháng làm xuất hiện số hạng ảo $i\omega Z_k$, ma trận hệ thống biến thành phức và phương pháp vectơ phân cực của Taziev-Ting hoàn toàn không thể sử dụng, buộc phải quay lại phương pháp khử ẩn truyền thống vốn bế tắc đối với vật liệu dị hướng bậc cao như monoclinic.

Luận án của Trịnh Thị Thanh Huệ đã phá vỡ thế bế tắc này bằng việc phát minh ra "Phương pháp vectơ phân cực phức". So sánh trực tiếp với hai nghiên cứu quốc tế tiêu biểu:

  • So với Godoy et al. (2012, Wave Motion): Luận án không chỉ giải quyết trọn vẹn trường hợp biên trở kháng tổng quát ($\sigma_{22} \neq 0$) cho môi trường đẳng hướng mà còn mở rộng thành công sang môi trường trực hướng, monoclinic, môi trường có ứng suất trước và môi trường quay.
  • So với Destrade (2001, 2003, J. Acoust. Soc. Am.; Proc. R. Soc. Lond. A): Trong khi Destrade chỉ dừng lại ở bài toán biên tự do ứng suất đối với vật liệu monoclinic và môi trường quay, luận án đã mở rộng toàn bộ hệ phương trình Stroh sang miền ma trận phức dưới tác động của trường trở kháng biên $Z_1, Z_2 \neq 0$.

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

┌──────────────────────────────────────────────────────────────────────────────────────────────────┐
│                                   KHUNG PHÂN TÍCH TOÁN - CƠ HỌC ĐỘT PHÁ                          │
├──────────────────────────────────────────────────────────────────────────────────────────────────┤
│ 1. PHÁT BIỂU STROH PHỨC (STEREOTYPE STROH EQUATIONS):                                            │
│    $\mathbf{w}' = i\mathbf{Q}\mathbf{w}, \quad 0 \le y < +\infty, \quad \mathbf{w}(+\infty) = 0$ │
│                                                                                                  │
│ 2. HỆ THỨC CƠ BẢN DẠNG HERMITE (PROPOSITION 2.1):                                                │
│    $\bar{\mathbf{Y}}^T(y) \hat{\mathbf{I}} \mathbf{P}^n \mathbf{Y}(y) = 0,                       │
│    \quad \forall n \in \mathbb{Z}, \quad \forall y \in [0, +\infty)$                             │
│                                                                                                  │
│ 3. KHỬ ẨN BIÊN ĐỘ QUA ĐỊNH LÝ CAYLEY-HAMILTON:                                                   │
│    $\bar{\mathbf{u}}^T(0) \mathbf{Q}_3^{(n)} \mathbf{u}(0) = 0, \quad n \in \{-1, 1, 2\}$        │
│                                                                                                  │
│ 4. PHƯƠNG TRÌNH TÁN SẮC DẠNG HIỆN (EXPLICIT SECULAR EQUATION):                                   │
│    $D_1^2 + D_2^2 - 4D D_3 = 0 \quad \text{(Đa thức đại số bậc 8 đối với } x = \rho c^2/c_{66})   │
└──────────────────────────────────────────────────────────────────────────────────────────────────┘

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

Luận án đã mở rộng và hoàn thiện Lý thuyết ma trận Stroh và Lý thuyết vectơ phân cực của Taziev-Ting thông qua Mệnh đề 2.1 (Hệ thức cơ bản cho ma trận phức):

"Mệnh đề 2.1: Nếu vectơ $2m$ chiều $\mathbf{Y}(y)$ là nghiệm của bài toán $\mathbf{Y}' = i\mathbf{P}\mathbf{Y}$ ($0 \le y < +\infty$, $\mathbf{Y}(+\infty) = 0$), trong đó ma trận khối $\mathbf{P} = \begin{bmatrix} \mathbf{P}_1 & \mathbf{P}_2 \ \mathbf{P}_3 & \mathbf{P}_4 \end{bmatrix}$ thỏa mãn $\mathbf{P}_2 = \bar{\mathbf{P}}_2^T, \mathbf{P}_3 = \bar{\mathbf{P}}_3^T, \mathbf{P}_4 = \bar{\mathbf{P}}_1^T$, thì ta có hệ thức cơ bản:" $$\bar{\mathbf{Y}}^T(0)\hat{\mathbf{I}}\mathbf{P}^n \mathbf{Y}(0) = 0, \quad \forall n \in \mathbb{Z}, \quad \text{với } \hat{\mathbf{I}} = \begin{bmatrix} \mathbf{0} & \mathbf{I} \ \mathbf{I} & \mathbf{0} \end{bmatrix}$$

Chứng minh toán học chặt chẽ dựa trên tính chất đối xứng liên hợp Hermite của toán tử vi phân và tính khả nghịch của ma trận luỹ thừa $\mathbf{P}^n$, chỉ ra rằng tích phân năng lượng $\frac{d}{dy}\left[\bar{\mathbf{Y}}^T \hat{\mathbf{I}}\mathbf{P}^n \mathbf{Y}\right] = 0$ bảo toàn trên toàn miền $[0, +\infty)$ và triệt tiêu tại vô cùng.

Mô hình lý thuyết thiết lập 5 phương trình tán sắc giải tích mẫu mực:

  1. Bán không gian trực hướng nén được có trở kháng: $$\sqrt{x}(e_1 - x)(1 - \delta_1\delta_2) + \left[e_3^2 - e_2(e_1 - x) - \delta_1\delta_2 x\right]\sqrt{P} + \left[\delta_1 e_2\sqrt{P} + \delta_2(e_1 - x)\right]\sqrt{x}\sqrt{S + 2\sqrt{P}} = 0$$ với các tham số không thứ nguyên $x = \rho c^2/c_{66}$, $e_1 = c_{11}/c_{66}$, $e_2 = c_{22}/c_{66}$, $e_3 = c_{12}/c_{66}$, $\delta_n = Z_n/\sqrt{\rho c_{66}}$, $P = (e_1-x)(1-x)/e_2$, $S = [e_2(e_1-x) + 1 - x - (1+e_3)^2]/e_2$.
  2. Bán không gian monoclinic nén được có trở kháng: Phương trình đại số bậc 8 dạng hiện: $$D_1^2 + D_2^2 - 4DD_3 = 0$$ trong đó $D, D_1, D_2, D_3$ là các định thức cấp 3 tạo bởi các phần tử thực và ảo của ma trận khối $\mathbf{Q}_3^{(n)}$ với $n \in {-1, 1, 2}$.
  3. Bán không gian trực hướng không nén được có trở kháng: $$(\delta - x)\sqrt{1 - x} + (\delta_1\delta_2 - 1)x = (\delta_1\sqrt{1 - x} + \delta_2)\sqrt{x}\sqrt{\delta - 2 - x + 2\sqrt{1 - x}}$$ với $\delta = (c_{11} - 2c_{12} + c_{22})/c_{66}$.

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

Khung phân tích tích hợp 3 trụ cột lý thuyết: Hình thức luận ma trận Stroh 4D/6D, Đại số toán tử đa thức ma trận (Định lý Cayley-Hamilton) và Lý thuyết điều kiện biên hiệu dụng tiệm cận. Cách tiếp cận mới này loại bỏ sự cần thiết phải tìm nghiệm riêng của phương trình đặc trưng bậc 4/bậc 6 đầy đủ, cho phép giải trực tiếp liên hệ biên độ chuyển vị phân cực $\mathbf{u}(0)$ thông qua các bất biến ma trận cấu tạo.

Điều kiện biên và phạm vi hợp lệ: Nghiên cứu xác lập chặt chẽ điều kiện tồn tại sóng mặt (sóng tắt dần vào sâu trong lòng bán không gian) là vận tốc pha $c$ phải thỏa mãn điều kiện cận trên âm học: $$0 < x < \min{1, e_1} \quad \Longleftrightarrow \quad 0 < \rho c^2 < \min{c_{66}, c_{11}}$$

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

Thiết kế nghiên cứu

Luận án được thiết kế theo thế giới quan thực chứng (Positivism) và cơ sở nhận thức luận diễn dịch - toán lý mẫu mực (Deductive-Mathematical Physics). Quy trình nghiên cứu là sự kết hợp chặt chẽ giữa suy diễn giải tích giải mã giải pháp cơ học môi trường liên tục và kiểm chứng nghiệm độc lập thông qua tính toán biểu trưng (Symbolic Computer Algebra).

┌────────────────────────────────────────────────────────────────────────────────────────┐
│                              QUY TRÌNH NGHIÊN CỨU GIẢI TÍCH RÕ RÀNG                    │
├────────────────────────────────────────────────────────────────────────────────────────┤
│ 1. THIẾT LẬP HỆ PHƯƠNG TRÌNH ĐỘNG LỰC HỌC NAVIER-CAUCHY & PHƯƠNG TRÌNH TRẠNG THÁI     │
│    $\sigma_{ij,j} = \rho \ddot{u}_i, \quad \sigma_{ij} = c_{ijkl} u_{k,l}$              │
│                                   │                                                    │
│                                   ▼                                                    │
│ 2. BIỂU DIỄN DẠNG MA TRẬN VI PHÂN CẤP 1 & PHÁT BIỂU STROH                              │
│    $\boldsymbol{\xi}' = i\mathbf{N}\boldsymbol{\xi}, \quad \boldsymbol{\xi} = [u_1, u_2, t_1, t_2]^T│
│                                   │                                                    │
│                                   ▼                                                    │
│ 3. TÍCH HỢP ĐIỀU KIỆN BIÊN TRỞ KHÁNG HOẶC HIỆU DỤNG LỚP MỎNG                           │
│    $\mathbf{t} = \mathbf{A}\mathbf{u}$ tại $y = 0 \Longrightarrow \boldsymbol{\sigma} = \mathbf{t} - \mathbf{A}\mathbf{u} = \mathbf{0}$ │
│                                   │                                                    │
│                                   ▼                                                    │
│ 4. ÁP DỤNG PHƯƠNG PHÁP VECTƠ PHÂN CỰC PHỨC (MỆNH ĐỀ 2.1)                               │
│    $\bar{\mathbf{u}}^T(0)\mathbf{Q}_3^{(n)}\mathbf{u}(0) = 0 \quad (n = -1, 1, 2)$     │
│                                   │                                                    │
│                                   ▼                                                    │
│ 5. THIẾT LẬP PHƯƠNG TRÌNH TÁN SẮC DẠNG HIỆN $D_1^2 + D_2^2 - 4DD_3 = 0$                │
│                                   │                                                    │
│                                   ▼                                                    │
│ 6. QUY DẪN GIỚI HẠN VỀ CÁC KẾT QUẢ KINH ĐIỂN ĐÃ BIẾT ĐỂ XÁC MINH (BENCHMARKING)       │
└────────────────────────────────────────────────────────────────────────────────────────┘

Quy trình nghiên cứu chuẩn xác

  • Thiết lập phương trình trường: Sử dụng hệ tọa độ Descartes $(x_1, x_2, x_3)$, phương truyền sóng theo trục $x_1$, trục $x_2$ hướng vào lòng bán không gian ($x_2 \ge 0$), trạng thái biến dạng phẳng ($u_3 \equiv 0, \partial/\partial x_3 \equiv 0$).
  • Chuẩn hóa ma trận vi phân: Chuyển đổi hệ phương trình đạo hàm riêng cấp hai sang hệ phương trình vi phân ma trận cấp một: $$\frac{d\boldsymbol{\zeta}}{dx_2} = \mathbf{M}\boldsymbol{\zeta}, \quad \boldsymbol{\zeta} = [u_1, u_2, \sigma_{12}, \sigma_{22}]^T$$
  • Biến đổi Stroh và ma trận trở kháng: Đặt $y = kx_2$, tách biến sóng điều hòa $e^{ik(x_1 - ct)}$, đưa về hệ $\boldsymbol{\xi}' = i\mathbf{N}\boldsymbol{\xi}$. Chuyển đổi ma trận Stroh phức $\mathbf{Q}$ thông qua ma trận trở kháng thuần ảo $\mathbf{A} = \operatorname{diag}(i\delta_1\sqrt{c_{66}X}, i\delta_2\sqrt{c_{66}X})$, trong đó $\bar{\mathbf{A}}^T = -\mathbf{A}$.
  • Phương pháp kiểm tra chéo (Triangulation): Mọi phương trình tán sắc thiết lập bằng Phương pháp vectơ phân cực phức đều được kiểm chứng chéo độc lập bằng Phương pháp giải tích truyền thống (khử trực tiếp các hằng số tích phân $B_1, B_2$ qua định thức ma trận hệ số).

Dữ liệu và Phân tích

  • Bộ tham số vật liệu đàn hồi chuẩn: Luận án tiến hành phân tích số trị trên các hệ vật liệu anisotropic tiêu chuẩn quốc tế: Tinh thể Silicon, Tinh thể Cadmium Selenide (CdSe), Hợp kim Titan, Vật liệu composite gia cố sợi Carbon và Cao su phi tuyến siêu đàn hồi nén được/không nén được.
  • Công cụ phần mềm giải tích: Toàn bộ các định thức ma trận bậc cao ($3 \times 3$ và $4 \times 4$ trên trường số phức) và các khai triển chuỗi tiệm cận bậc hai theo độ dày lớp mỏng $\epsilon = kh$ được tính toán giải tích biểu trưng (Symbolic Computation) chính xác tuyệt đối trên phần mềm Wolfram Mathematica và Maple, triệt tiêu hoàn toàn sai số làm tròn số học (truncation errors).

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

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

┌──────────────────────────────────────────────────────────────────────────────────────────────────┐
│                                   5 PHÁT HIỆN ĐỘT PHÁ CỦA LUẬN ÁN                                │
├──────────────────────────────────────────────────────────────────────────────────────────────────┤
│ 1. Phát minh Phương pháp vectơ phân cực phức: Giải quyết bài toán Stroh ma trận phức tổng quát.  │
│ 2. Thiết lập phương trình tán sắc giải tích tường minh bậc 8 cho tinh thể Monoclinic $x_3=0$.    │
│ 3. Khám phá cơ chế dịch chuyển vận tốc sóng Rayleigh theo 2 tham số trở kháng $\delta_1, \delta_2$. │
│ 4. Giải mã hiệu ứng bất đối xứng pha do tương tác đồng thời giữa kéo/nén và trượt cắt trước.    │
│ 5. Dẫn xuất phương trình tiệm cận bậc hai $(kh)^2$ không chứa hằng số phụ cho nền và lớp phủ.   │
└──────────────────────────────────────────────────────────────────────────────────────────────────┘
  1. Tổng quát hóa thành công Phương pháp vectơ phân cực: Chứng minh định lý đại số ma trận cho phép giải hệ Stroh phức, biến phương pháp vectơ phân cực thành một công cụ vạn năng không chỉ cho bài toán biên tự do mà cho toàn bộ các bài toán biên trở kháng, biên nhớt và tiếp xúc đàn-chất lỏng.
  2. Tìm ra phương trình tán sắc tường minh chính xác bậc 8 cho môi trường Monoclinic: Lần đầu tiên trong cơ học sóng mặt, phương trình thế thế vận tốc pha không thứ nguyên $x$ cho tinh thể monoclinic chịu trở kháng biên được công bố dưới dạng phương trình đa thức đại số $D_1^2 + D_2^2 - 4DD_3 = 0$.
  3. Quy luật điều tiết vận tốc sóng của tham số trở kháng biên: Phân tích số trị chỉ ra rằng khi tham số trở kháng trượt $\delta_1$ và trở kháng pháp $\delta_2$ tăng từ $0$ lên $+\infty$, vận tốc không thứ nguyên $x = \rho c^2/c_{66}$ giảm đơn điệu và tiến tiệm cận về $0$, chứng minh rằng trở kháng bề mặt làm tăng độ mềm động lực học của lớp biên.
  4. Hiệu ứng ghép đôi giữa ứng suất cắt trước và sóng quay: Khi bán không gian chịu đồng thời ứng suất kéo/nén thuần nhất và ứng suất trượt cắt $\sigma_{12}^0 \neq 0$, phương trình tán sắc xuất hiện các số hạng phi đối xứng bậc lẻ, dẫn đến vận tốc truyền sóng Rayleigh theo chiều dương và chiều âm của trục $x_1$ không còn bằng nhau ($c^+ \neq c^-$).
  5. Phương trình xấp xỉ bậc hai đóng cho bán không gian phủ lớp mỏng: Luận án dẫn ra phương trình tán sắc xấp xỉ bậc hai đối với tham số độ dày không thứ nguyên $\epsilon = kh$: $$F_0(x) + (kh) F_1(x) + (kh)^2 F_2(x) = 0$$ hoàn toàn không chứa bất kỳ hằng số bất định nào, vượt trội so với xấp xỉ bậc 4 của Achenbach & Keshava (1967).

Implications đa chiều

  • Đóng góp lý thuyết: Hoàn thiện bức tranh lý thuyết trường động lực học đàn hồi dị hướng, bắc nhịp cầu kết nối giữa cơ học môi trường liên tục cổ điển và lý thuyết tán xạ sóng trong các hệ biên tiêu tán phức tạp.
  • Đổi mới phương pháp luận: Cung cấp quy trình mẫu chuẩn mực 5 bước (PDE $\to$ Stroh Ma trận phức $\to$ Hermite cơ bản $\to$ Cayley-Hamilton $\to$ Secular Equation) áp dụng trực tiếp cho các bài toán mở rộng như sóng Love, sóng Lamb, sóng Stoneley trong tinh thể áp điện (piezoelectric) và từ-áp điện (magnetoelectroelastic).
  • Ứng dụng thực tiễn kỹ thuật: Là cơ sở thuật toán trực tiếp để giải "Bài toán ngược" (Inverse Problem) trong NDT: đo đạc thực nghiệm phổ tán sắc vận tốc sóng $c(\omega)$ để trích xuất chính xác môđun đàn hồi $c_{ij}$, độ dày lớp phủ nano $h$ và trạng thái ứng suất dư của linh kiện bán dẫn mà không cần phá hủy mẫu.
  • Khuyến nghị địa kỹ thuật và quan trắc công trình: Cho phép các kỹ sư địa chấn mô hình hóa chính xác hơn phản ứng nền đất chịu tải trọng động đất khi bề mặt tiếp xúc với nền móng công trình hoặc lớp trầm tích mềm (được biểu diễn dưới dạng điều kiện biên trở kháng).

Limitations và Future Research

Nhằm duy trì tính khách quan khoa học và chuẩn mực học thuật cao nhất, luận án thẳng thắn thừa nhận các giới hạn nội tại:

  1. Giới hạn giả thiết tuyến tính: Nghiên cứu giới hạn trong phạm vi lý thuyết biến dạng đàn hồi tuyến tính đối với trường dao động sóng vi mô, bỏ qua các hiệu ứng phi tuyến vật liệu bậc cao khi biên độ kích động sóng vượt ngưỡng đàn hồi.
  2. Giới hạn biến dạng phẳng 2D: Bài toán được khảo sát dưới điều kiện biến dạng phẳng hai chiều ($u_3 \equiv 0$), chưa xét đến sự ghép kênh không gian ba chiều đầy đủ giữa sóng mặt Rayleigh và sóng trượt phân cực ngang Love/SH trong tinh thể dị hướng tam tà (triclinic).
  3. Giới hạn xấp xỉ tiệm cận lớp mỏng: Phương trình tán sắc xấp xỉ cho cấu trúc phủ lớp mỏng được khai triển đến bậc hai $(kh)^2$, do đó độ chính xác bị suy giảm khi độ dày lớp phủ tiệm cận bước sóng ($kh > 1$).
  4. Thiếu hụt dữ liệu thực nghiệm siêu âm thực tế: Luận án tập trung thuần túy vào giải tích toán học và mô phỏng số biểu trưng, chưa thực hiện các đo đạc thực nghiệm sóng âm bề mặt trên mẫu vật lý trong phòng thí nghiệm.

Chương trình nghiên cứu tương lai (Future Research Agenda):

  • Hướng 1: Mở rộng phương pháp vectơ phân cực phức cho môi trường đàn nhớt phân số (fractional viscoelasticity) và môi trường rỗng bão hòa chất lỏng (Biot's poroelasticity).
  • Hướng 2: Thiết lập bài toán sóng Rayleigh 3D trong tinh thể có tính dị hướng tổng quát (Anisotropic Triclinic Media).
  • Hướng 3: Tích hợp trí tuệ nhân tạo (Physics-Informed Neural Networks - PINNs) kết hợp với các phương trình tán sắc tường minh của luận án để nhận dạng tự động ứng suất trước từ tín hiệu cảm biến SAW.

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

┌─────────────────────────────────────────────────────────────────────────────────────────────────┐
│                                   MA TRẬN TÁC ĐỘNG CỦA LUẬN ÁN                                  │
├────────────────────────────────┬────────────────────────────────────────────────────────────────┤
│ LĨNH VỰC TÁC ĐỘNG              │ ĐÓNG GÓP VÀ CHỈ SỐ ĐỊNH LƯỢNG                                  │
├────────────────────────────────┼────────────────────────────────────────────────────────────────┤
│ 1. Viện hàn lâm & Học thuật    │ Thiết lập phương pháp giải tích mới; tiềm năng trích dẫn cao   │
│                                │ trong các tạp chí Q1 về Cơ học sóng (Wave Motion, JASA).       │
├────────────────────────────────┼────────────────────────────────────────────────────────────────┤
│ 2. Công nghiệp MEMS & SAW      │ Tối ưu hóa thiết kế bộ lọc tần số sóng âm bề mặt trong vi     │
│                                │ mạch truyền thông không dây 5G/6G và cảm biến con quay vi cơ.  │
├────────────────────────────────┼────────────────────────────────────────────────────────────────┤
│ 3. Công nghệ màng mỏng & NDT   │ Cơ sở thuật toán đo đạc không phá hủy độ dày và chất lượng bám │
│                                │ dính của màng phủ quang học, lớp phủ chịu nhiệt bán dẫn.       │
├────────────────────────────────┼────────────────────────────────────────────────────────────────┤
│ 4. Địa chấn & Địa kỹ thuật     │ Hoàn thiện mô hình dự báo lan truyền năng lượng sóng động đất  │
│                                │ qua các tầng địa chất chịu ứng suất kiến tạo và trở kháng bùn. │
└────────────────────────────────┴────────────────────────────────────────────────────────────────┘
  • Ảnh hưởng học thuật quốc tế: Mở ra một phương pháp tiếp cận mới trong cơ học giải tích, giải quyết các bài toán biên phức mà các phương pháp truyền thống của Taziev, Ting và Destrade chưa thể xử lý.
  • Chuyển đổi công nghiệp vật liệu cao: Cung cấp công cụ toán học chính xác giúp các kỹ sư R&D trong ngành công nghiệp màng mỏng (Thin Solid Films) kiểm soát chất lượng cơ lý của lớp phủ nano trên nền tinh thể dị hướng.
  • Ý nghĩa quốc gia: Nâng cao vị thế của trường phái Cơ học giải tích Việt Nam (được dẫn dắt bởi GS. Phạm Chí Vĩnh) trên bản đồ cơ học lý thuyết thế giới.

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

  • Nghiên cứu sinh và Giảng viên Cơ học Vật thể rắn: Được trang bị một công cụ giải tích toán học hoàn chỉnh (Phương pháp vectơ phân cực phức) để giải quyết các bài toán truyền sóng phức tạp tiếp theo.
  • Các nhà khoa học Âm học & Địa vật lý: Sở hữu các phương trình tán sắc giải tích chuẩn tắc dạng đóng để đưa vào các phần mềm mô phỏng địa chấn và tính toán hàm Green trong môi trường nửa không gian.
  • Kỹ sư R&D thiết bị Vi điện tử (MEMS/Acoustic Filters): Ứng dụng trực tiếp các phương trình vận tốc pha sóng mặt để thiết kế bộ cộng hưởng âm học, cảm biến đo áp suất và con quay hồi chuyển vi cơ học với độ chính xác cao.
  • Chuyên gia Đánh giá không phá hủy (NDT): Nhận được các công thức giải tích giải bài toán ngược để xác định hằng số đàn hồi vật liệu composite và ứng suất dư trong kết cấu công trình ngầm và hàng không vũ trụ.

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

Đóng góp lý thuyết độc đáo và đột phá nhất là việc thiết lập Phương pháp vectơ phân cực phức thông qua việc phát biểu và chứng minh Mệnh đề 2.1. Phương pháp này mở rộng trực tiếp Hình thức luận ma trận Stroh (Stroh, 1958) và Lý thuyết vectơ phân cực của Taziev (1989), Ting (1996, 2002) và Collet & Destrade (2004). Trước luận án này, phương pháp vectơ phân cực chỉ áp dụng được cho ma trận Stroh thực; luận án đã chứng minh rằng khi ma trận Stroh $\mathbf{P}$ là ma trận phức thỏa mãn tính chất tự liên hợp Hermite mở rộng ($\mathbf{P}_2 = \bar{\mathbf{P}}_2^T, \mathbf{P}_3 = \bar{\mathbf{P}}_3^T, \mathbf{P}_4 = \bar{\mathbf{P}}_1^T$), hệ thức cơ bản $\bar{\mathbf{Y}}^T(0)\hat{\mathbf{I}}\mathbf{P}^n\mathbf{Y}(0) = 0$ vẫn hoàn toàn nghiệm đúng $\forall n \in \mathbb{Z}$.

2. Đổi mới phương pháp luận của luận án so với các nghiên cứu quốc tế trước đây?

So với nghiên cứu của Godoy, Durán & Nkemzi (2012) và Destrade (2001, 2003):

  • Godoy et al. (2012) dùng phương pháp hàm biến phức truyền thống trên bán không gian đẳng hướng và bắt buộc phải gán $\sigma_{22} = 0$ để tránh bế tắc trong việc giải hệ phương trình tích phân. Luận án giải quyết trọn vẹn cả $\sigma_{12} \neq 0$ và $\sigma_{22} \neq 0$ cho các môi trường có tính dị hướng cao nhất (Monoclinic).
  • Destrade (2001, 2003) giải phương trình tán sắc cho vật liệu monoclinic và bán không gian quay dựa trên phương pháp tích phân đầu và vectơ phân cực ma trận thực cho điều kiện biên tự do ứng suất. Luận án tích hợp thành công điều kiện biên trở kháng phức vào phát biểu Stroh ma trận phức cấp 4 và cấp 6.

3. Phát hiện nào gây bất ngờ nhất về mặt cơ học giải tích trong luận án?

Phát hiện bất ngờ nhất nằm ở cấu trúc đại số của phương trình tán sắc cho vật liệu trực hướng nén được chịu trở kháng biên (Phương trình 2.35). Khi bình phương hai lần để khử căn bậc hai, phương trình này dẫn về một đa thức bậc 6 (Phương trình 2.84). Điều này chứng minh rằng sự hiện diện của trở kháng biên $\delta_1, \delta_2$ đã nâng bậc đa thức tán sắc từ bậc 3 (phương trình bậc 3 kinh điển của Rayleigh-Chadwick) lên thành đa thức đại số bậc 6 đóng, phản ánh bản chất tương tác năng lượng phi tuyến giữa dao động tiếp tuyến và pháp tuyến bề mặt.

4. Luận án có cung cấp quy trình tái lập kết quả (Replication Protocol) rõ ràng không?

Có. Quy trình giải tích được mô tả với độ chi tiết toán học chuẩn mực:

  • Bảng ma trận cấu tạo hệ số đàn hồi $c_{ij}$ và các biến đổi vi phân $M_k, N_k, Q_k$.
  • Biểu thức giải tích của các định thức $D, D_1, D_2, D_3$ thông qua các thành phần ma trận $Q_{ij}^{(n)}$ với $n \in {-1, 1, 2}$.
  • Toàn bộ các bước rút gọn đại số được kiểm chứng bằng mã nguồn tính toán biểu trưng trên phần mềm Wolfram Mathematica.

5. Lộ trình phát triển nghiên cứu 10 năm tiếp theo từ nền tảng luận án?

  1. Năm 1-3: Tích hợp trường đa vật lý điện-từ-nhiệt-đàn hồi (Magneto-Electro-Thermo-Elasticity) vào phát biểu Stroh phức.
  2. Năm 4-6: Nghiên cứu sóng Rayleigh truyền trong môi trường rỗng bất đồng nhất có gradient chức năng (FGMs) dưới điều kiện biên tiếp xúc không hoàn hảo.
  3. Năm 7-10: Xây dựng phần mềm thương mại giải bài toán ngược NDT ứng dụng thuật toán học sâu (Deep Learning) huấn luyện trên cơ sở dữ liệu giải tích chính xác của luận án.

Kết luận

  1. Phát triển thành công Phương pháp vectơ phân cực phức: Giải quyết trọn vẹn bài toán xác định phương trình tán sắc giải tích cho các hệ phương trình vi phân ma trận Stroh có hệ số phức, mở rộng biên giới của cơ học giải tích hiện đại.
  2. Thiết lập hệ phương trình tán sắc dạng hiện chính xác cho môi trường dị hướng: Thu được phương trình đại số đóng bậc 6 và bậc 8 cho bán không gian đàn hồi trực hướng và monoclinic ($x_3=0$) nén được và không nén được chịu điều kiện biên trở kháng tổng quát.
  3. Giải mã động lực học môi trường có ứng suất trước và quay: Xây dựng phương trình tán sắc giải tích chính xác cho bán không gian chịu ứng suất kéo-nén thuần túy, kéo-nén kết hợp trượt cắt và môi trường quay chịu lực quán tính Coriolis kết hợp gia cố cốt sợi.
  4. Xây dựng điều kiện biên hiệu dụng và phương trình tán sắc xấp xỉ bậc hai cho lớp phủ mỏng: Dẫn ra công thức xấp xỉ tường minh theo $(kh)^2$ mô tả tương tác động lực giữa nền monoclinic và lớp phủ monoclinic mà không phụ thuộc vào hằng số kinh nghiệm.
  5. Đóng góp vào kho tàng tri thức cơ học quốc tế: Công trình khẳng định năng lực nghiên cứu đỉnh cao của nền cơ học lý thuyết Việt Nam, cung cấp hệ thống giải pháp toán học mẫu mực ứng dụng trong công nghệ vi cơ điện tử, địa chấn học và đánh giá không phá hủy vật liệu.