Tổng quan luận án

Đề tài "Sóng mặt và sóng trong các cấu trúc mỏng" do nghiên cứu sinh Nguyễn Thị Khánh Linh thực hiện tại Viện Cơ học thuộc Viện Hàn lâm Khoa học và Công nghệ Việt Nam (2013), dưới sự hướng dẫn khoa học của PGS. TS Phạm Chí Vĩnh, thuộc chuyên ngành Cơ học vật rắn (mã số 62 44 21 01).

Sóng mặt Rayleigh truyền trong môi trường đàn hồi là cơ sở lý thuyết cho nhiều ứng dụng trong địa chấn học, âm học, địa vật lý, công nghệ truyền thông, thiết kế các bộ lọc sóng âm bề mặt (SAW) trong thiết bị di động, và đặc biệt là kỹ thuật đánh giá không phá hủy (NDT) các đặc trưng cơ học của vật liệu và kết cấu. Trong kỹ thuật hiện đại, các cấu trúc lớp mỏng phủ trên bán không gian và vật liệu composite phân lớp tuần hoàn (như composite cốt sợi dùng trong vỏ tàu thủy, máy bay, ô-tô) được ứng dụng rộng rãi. Để xác định các tính chất cơ học hoặc ứng suất dư của cấu trúc từ dữ liệu thực nghiệm đo vận tốc sóng (giải bài toán ngược), việc thiết lập các phương trình tán sắc và công thức giải tích xác định vận tốc sóng dạng tường minh là yêu cầu then chốt.

Khoảng trống nghiên cứu được luận án chỉ ra bao gồm:

  1. Đối với bán không gian đẳng hướng không nén được chịu tác dụng của trọng trường, các nghiên cứu trước đây dừng lại ở việc thiết lập phương trình tán sắc mà chưa đưa ra công thức giải tích chính xác và công thức xấp xỉ có độ chính xác cao cho vận tốc sóng.
  2. Đối với bản mỏng đàn hồi trực hướng bán vô hạn, phương trình tán sắc sóng Rayleigh chính trước đây bị chia thành nhiều biểu thức cồng kềnh khác nhau tùy theo nhóm vật liệu; trong khi phương trình tán sắc sóng Rayleigh không chính mới chỉ dừng lại ở dạng ẩn.
  3. Đối với bán không gian đàn hồi nằm dưới lớp nước, bài toán truyền sóng trước đây chỉ được giải quyết dưới giả thiết lớp nước nông và trọng trường nhỏ; bài toán tổng quát với độ dày và trọng trường tùy ý chưa được hoàn thiện.
  4. Đối với bán không gian phủ lớp mỏng, các phương trình tán sắc xấp xỉ hiện có hầu hết dừng lại ở bậc một hoặc bậc hai, và chủ yếu giới hạn trong môi trường đàn hồi đẳng hướng.
  5. Đối với môi trường phân lớp tuần hoàn, các khai triển tiệm cận phương trình tán sắc sóng SH mới chỉ tính đến bậc thấp $\Omega_1, \Omega_3$; chưa có công thức tổng quát bậc cao cho sóng SH và sóng Lamb khi có ứng suất trước.

Luận án xác định 3 mục tiêu nghiên cứu cụ thể:

  • Mục tiêu 1: Áp dụng các công cụ giải tích mới để phát triển kết quả của các bài toán sóng mặt Rayleigh:
    • Bài toán 1: Tìm công thức chính xác và các công thức xấp xỉ độ chính xác cao cho vận tốc sóng Rayleigh trong môi trường đàn hồi đẳng hướng không nén được chịu tác dụng của trọng trường; khảo sát điều kiện tồn tại và tính duy nhất của sóng.
    • Bài toán 2: Thiết lập phương trình tán sắc dạng tường minh thống nhất cho sóng chính Rayleigh trong bản mỏng trực hướng bán vô hạn; xây dựng công thức vận tốc chính xác và xấp xỉ; tìm phương trình tán sắc dạng tường minh cho sóng không chính Rayleigh bằng phương pháp tích phân đầu.
    • Bài toán 3: Nghiên cứu bài toán tổng quát sóng Rayleigh truyền trong bán không gian đàn hồi không nén được nằm dưới lớp nước không nén được với độ dày và trọng trường tùy ý; thiết lập điều kiện tồn tại, phương trình tán sắc và các công thức vận tốc xấp xỉ.
  • Mục tiêu 2: Xây dựng phương trình tán sắc xấp xỉ bậc ba của sóng Rayleigh truyền trong bán không gian phủ lớp mỏng cho hai trường hợp: (i) Bán không gian và lớp đều là đàn hồi trực hướng (nén được hoặc không nén được); (ii) Bán không gian và lớp đều có biến dạng trước.
  • Mục tiêu 3: Mở rộng nghiên cứu của Norris và Santosa (1992), thiết lập công thức truy hồi xác định các hệ số tiệm cận bậc cao $\Omega_{2m+1}$ ($m \ge 2$) cho sóng SH trong môi trường phân lớp tuần hoàn đẳng hướng nén được và sóng Lamb trong môi trường phân lớp tuần hoàn không nén được có ứng suất trước.

Đối tượng nghiên cứu là sóng mặt Rayleigh, sóng SH và sóng Lamb truyền trong các bán không gian đàn hồi, bán không gian phủ lớp mỏng và môi trường phân lớp tuần hoàn. Phạm vi nghiên cứu giới hạn ở việc mô hình hóa giải tích cơ học môi trường liên tục biến dạng phẳng/ứng suất phẳng, dẫn xuất các phương trình tán sắc chính xác, phương trình tán sắc xấp xỉ và công thức vận tốc sóng.

Tổng quan tài liệu và vị trí của luận án

Sự phát triển lý thuyết sóng mặt Rayleigh khởi đầu từ công trình của Rayleigh (1885). Trong hơn hai thập kỷ trước thời điểm thực hiện luận án, lĩnh vực này phát triển với sự ra đời của các công cụ giải tích mới:

Hướng nghiên cứu Các công trình tiêu biểu Kết quả và hạn chế đã có
Phương pháp giải tích sóng mặt Taziev (1989), Ting (2004), Mozhaev (1995), Destrade (2001, 2003) Đề xuất phương pháp véctơ phân cực và phương pháp tích phân đầu để vượt qua hạn chế của phương trình đặc trưng trong môi trường dị hướng.
Công thức chính xác vận tốc sóng Rahman & Barber (1995), Nkemzi (1997), Malischewsky (2000, 2004), Phạm Chí Vĩnh & Ogden (2004), Ogden & Phạm Chí Vĩnh (2004), Phạm Chí Vĩnh (2009, 2010), Phạm Chí Vĩnh & Phạm Thị Hà Giang (2010) Rahman-Barber dùng phương trình bậc ba nhưng công thức tách thành hai nhánh theo dấu biệt thức. Malischewsky dùng Cardan nhưng chưa chứng minh. Phạm Chí Vĩnh và cộng sự chứng minh chặt chẽ và mở rộng cho vật liệu trực hướng, biến dạng trước.
Công thức xấp xỉ vận tốc sóng Bergmann (1948), Achenbach (1973), Brekhovskikh (1960), Briggs (1992), Nesvijski (2000), Li (2002), Phạm Chí Vĩnh & Malischewsky (2004, 2007, 2008) Các công thức ban đầu sai số còn lớn. Phạm Chí Vĩnh & Malischewsky áp dụng phương pháp bình phương tối thiểu để tối ưu hóa đa thức bậc hai cho hàm bậc ba.
Ảnh hưởng của trọng trường và chất lỏng Bromwich (1898), Biot (1940, 1965), Kuipers (1958), Phạm Chí Vĩnh & Seriani (2009), Phạm Chí Vĩnh (2009) Bromwich, Biot, Kuipers đưa ra phương trình tán sắc cho bán không gian không nén được chịu trọng trường nhưng dừng ở nghiệm số. Bromwich giải bài toán lớp nước dưới giả thiết lớp mỏng và trọng trường nhỏ.
Bản mỏng đàn hồi trực hướng Ohyoshi (2000), Cerv (2008), Cerv và cs (2010), Chadwick (1976) Ohyoshi chỉ tìm được phương trình dạng ẩn. Cerv phân tách thành các biểu thức riêng cho vật liệu dị hướng yếu và dị hướng mạnh lập phương; sóng không chính Rayleigh chỉ có phương trình ẩn.
Bán không gian phủ lớp mỏng Achenbach & Keshava (1967), Steigmann (2007), Trần Thanh Tuấn (2007) Xấp xỉ bậc bốn của Achenbach chứa hằng số chưa xác định. Steigmann dừng ở bậc hai cho lớp đẳng hướng ngang trên bán không gian đẳng hướng. Trần Thanh Tuấn dừng ở bậc một.
Cấu trúc mỏng tuần hoàn Norris & Santosa (1992) Thiết lập khai triển tiệm cận cho sóng SH trong môi trường tuần hoàn gồm $N$ lớp đẳng hướng nén được nhưng chỉ tính $\Omega_1, \Omega_3$.

Vị trí của luận án: Luận án kế thừa các phương pháp giải tích hiện đại (lý thuyết phương trình bậc ba, bình phương tối thiểu, tích phân đầu, điều kiện biên hiệu dụng, khai triển tiệm cận), giải quyết trọn vẹn các bài toán còn dang dở, cung cấp các biểu thức giải tích tường minh, hợp nhất và nâng cao bậc chính xác của các mô hình truyền sóng trong cấu trúc lớp mỏng.

Cơ sở lý thuyết và phương pháp nghiên cứu

Luận án dựa trên lý thuyết đàn hồi động lực học tuyến tính, lý thuyết biến dạng hữu hạn (cho môi trường có ứng suất trước) và cơ học môi trường liên tục dưới trạng thái ứng suất phẳng hoặc biến dạng phẳng.

Các khái niệm và mô hình lý thuyết

  • Mô hình đàn hồi chịu trọng trường: Sử dụng hệ phương trình chuyển động của Biot với các ứng suất giả định $s'{11} = s{11} + \rho g u_3$, $s'{33} = s{33} + \rho g u_3$, $s'{13} = s{13}$ và điều kiện không nén được $u_{1,1} + u_{3,3} = 0$.
  • Mô hình bản mỏng trực hướng: Trạng thái ứng suất phẳng ($\sigma_{31} = \sigma_{32} = \sigma_{33} = 0$). Quan hệ ứng suất - biến dạng biểu diễn qua các hệ số độ cứng $B_{ij}$ ($B_{11}, B_{22}, B_{12}, B_{66}$) hoặc ma trận độ cứng quay $Q_{ij}(\theta)$.
  • Điều kiện biên hiệu dụng bậc ba: Biểu diễn tương tác cơ học của lớp phủ mỏng lên bán không gian bằng cách khai triển chuỗi Taylor trường ứng suất và chuyển dịch theo độ dày không thứ nguyên $\varepsilon = kh \ll 1$ đến bậc ba.

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

  1. Phương pháp lý thuyết phương trình bậc ba: Khảo sát đại số phương trình tán sắc sau khi hữu tỷ hóa về phương trình bậc ba dạng $z^3 + 3qz - 2R = 0$, sử dụng công thức nghiệm Cardan để tìm công thức vận tốc chính xác, đồng thời khảo sát dấu của biệt thức $D = q^3 + R^2$ để chứng minh sự tồn tại và duy nhất của nghiệm sóng Rayleigh.
  2. Phương pháp bình phương tối thiểu: Thay thế hàm số bậc ba $x^3$ hoặc $y^3$ bằng đa thức bậc hai xấp xỉ tốt nhất trên các khoảng đóng xác định theo chuẩn không gian $L^2$ hoặc $C[a, b]$:
    • Trên khoảng $[x_0, 1]$: $x^3 \approx 2.73787x^2 - 2.47605x + 0.73818$.
    • Trên khoảng $[0, y_0]$ ($y_0 \approx 0.29560$): $y^3 \approx 0.44315y^2 - 0.04618y + 0.00078$.
    • Trên khoảng $[0, 1]$: $x^3 \approx 1.84375x^2 - 0.89063x + 0.04688$ (hoặc $1.875x^2 - 0.9375x + 0.03125$). Kỹ thuật này chuyển phương trình tán sắc bậc ba thành phương trình bậc hai, giải ra nghiệm vận tốc xấp xỉ tường minh.
  3. Phương pháp tích phân đầu (First Integral Method): Áp dụng cho bài toán sóng không chính Rayleigh trong bản mỏng trực hướng. Thiết lập hệ phương trình vi phân cấp hai dạng ma trận đối với biên độ ứng suất $\boldsymbol{\Sigma}$, đưa vào các tích phân vô hướng dẫn đến hệ ma trận phản đối xứng $D, E, F$, từ đó rút ra phương trình tán sắc tường minh mà không cần tìm nghiệm của phương trình đặc trưng bậc bốn.
  4. Phương pháp khai triển tiệm cận: Biểu diễn vận tốc và tần số dưới dạng chuỗi lũy thừa theo tham số nhỏ $\varepsilon = kh$ (với $c^2 = \Omega_1 + (kh)^2 \Omega_3 + \dots$), sử dụng tính chất hàm chẵn để triệt tiêu các bậc chẵn ($\Omega_{2m} = 0$), thiết lập hệ thức truy hồi xác định các hệ số bậc lẻ $\Omega_{2n+1}$.

Nguồn dữ liệu kiểm chứng

Luận án sử dụng thông số vật liệu thực tế để tính toán số trị và so sánh: Kevlar-Epoxy (Aramid), Wolfram, Fibredux-914-C, $\alpha$-Fe, Silicon, Crystalline Gold, hợp kim SE84LV, Boron-Epoxy. Kết quả giải tích được đối chiếu với kết quả số từ phương pháp phần tử hữu hạn (FEA) của Cerv và cộng sự (2010).

Nội dung chính theo từng chương

Chương 1: Tổng quan

Chương 1 hệ thống hóa lịch sử nghiên cứu sóng mặt Rayleigh và sóng trong cấu trúc mỏng trong hơn một thế kỷ. Tác giả phân tích chi tiết:

  • Sự tiến triển của các phương pháp thiết lập phương trình tán sắc dạng hiện (phương pháp véctơ phân cực của Taziev, Ting, Destrade; phương pháp tích phân đầu của Mozhaev, Destrade).
  • Các giai đoạn phát triển của công thức vận tốc sóng Rayleigh từ công thức Cardan của Rahman-Barber (1995), công thức hàm biến phức của Nkemzi (1997), nghiệm giải tích của Malischewsky (2000, 2004) đến các công thức chính xác và xấp xỉ của Phạm Chí Vĩnh, Ogden, Malischewsky (2004-2010).
  • Tổng quan về đánh giá ứng suất trước bằng sóng âm, phương trình sóng trong bán không gian phủ lớp mỏng và cấu trúc tuần hoàn.
  • Tổng kết tình hình nghiên cứu trong nước và xác định 3 mục tiêu cụ thể của luận án.

Chương 2: Sóng mặt Rayleigh

Chương này tập trung giải quyết hai bài toán lớn:

1. Sóng Rayleigh trong bán không gian đàn hồi đẳng hướng, không nén được, chịu tác dụng của trọng trường

  • Mô hình xét bán không gian $x_3 \ge 0$. Phương trình tán sắc dẫn xuất có dạng: $$\sqrt{(2 - x)^2 - 4\sqrt{1 - x}} - \delta x = 0$$ trong đó $x = c^2/c_2^2$ ($c_2 = \sqrt{\mu/\rho}$), tham số trọng trường $\delta = g/(k c_2^2) \ge 0$.
  • Điều kiện tồn tại và duy nhất: Sóng Rayleigh tồn tại khi và chỉ khi $0 \le \delta < 1$. Khi $\delta \ge 1$, phương trình không có nghiệm thực trong khoảng $(0, 1)$. Khi $0 \le \delta < 1$, nghiệm $x_r(\delta)$ là duy nhất và là hàm đơn điệu tăng thực sự theo $\delta$, biến thiên từ $x_0 \approx 0.91262$ (khi $\delta = 0$) đến $1$ (khi $\delta \to 1$).
  • Công thức chính xác: Luận án thiết lập 2 công thức chính xác qua lý thuyết phương trình bậc ba:
    • Công thức 1: $$x_r = \frac{4 + \delta}{3} + \sqrt[3]{R + \sqrt{D}} + \sqrt[3]{R - \sqrt{D}}$$ với $q = (\delta^2 + 8\delta - 8)/9$, $R = -(\delta^3 + 12\delta^2 + 12\delta + 136)/27$, $D = q^3 + R^2 > 0, \forall \delta \in [0, 1)$.
    • Công thức 2 (qua biến đổi $y = \sqrt{1-x}$): $$x_r = 1 - \left[ -\frac{1}{3} + \sqrt[3]{\frac{26 - 9\delta}{27} + \sqrt{\frac{(\delta + 11)(\delta + 8)^2}{27}}} + \sqrt[3]{\frac{26 - 9\delta}{27} - \sqrt{\frac{(\delta + 11)(\delta + 8)^2}{27}}} \right]^2$$
  • Công thức xấp xỉ: Bằng phương pháp bình phương tối thiểu, luận án dẫn ra 2 công thức xấp xỉ bậc hai:
    • $x_1(\delta) = \frac{B - \sqrt{B^2 - 4AC}}{2A}$ với $A = 2.73787 - 2(4+\delta)$, $B = 2.47605 - 20 - 8\delta - \delta^2$, $C = 16 - 0.73818$.
    • $x_2(\delta) = 1 - y_r^2$ với $y_r = \frac{-(2.04618 + \delta) + \sqrt{(2.04618 + \delta)^2 - 4(1.44315)(\delta - 0.99922)}}{2 \times 1.44315}$.

2. Sóng Rayleigh trong bản mỏng đàn hồi trực hướng, bán vô hạn

  • Sóng chính Rayleigh: Truyền theo trục chính $x_1$ của bản mỏng $x_2 \ge 0$. Phương trình tán sắc thống nhất cho mọi vật liệu trực hướng được thiết lập: $$\sqrt{B_{66} - \rho c^2} [B_{12}^2 - B_{22}(B_{11} - \rho c^2)] + \rho c^2 \sqrt{B_{22} B_{66} (B_{11} - \rho c^2)(B_{66} - \rho c^2)} = 0$$ Biểu thức này đơn giản hơn các phương trình tách rời của Cerv (2008). Vận tốc không thứ nguyên $x = \rho c^2 / B_{66}$ được biểu diễn chính xác và xấp xỉ qua 3 tham số không thứ nguyên $b_1 = E_2/E_1$, $b_2 = 1 - \nu_{12}^2 E_2/E_1$, $b_3 = 2(1 - \nu_{12} E_2/E_1) - E_2/G_{12}$.
  • Sóng không chính Rayleigh: Trục truyền sóng hợp với trục vật liệu chính góc $\theta$. Bằng phương pháp tích phân đầu, luận án thiết lập phương trình tán sắc tường minh bậc bốn theo $X = \rho c^2$: $$\det \begin{bmatrix} \alpha_{11} & \beta_{11} & \gamma_{11} \ \alpha_{12} & \beta_{12} & \gamma_{12} \ \alpha_{22} & \beta_{22} & \gamma_{22} \end{bmatrix} = 0$$ Chứng minh chặt chẽ về mặt toán học tính đối xứng của đường cong vận tốc: $c_R(\theta) = c_R(\pi - \theta)$, và khi $E_1 = E_2$ thì $c_R(\theta) = c_R(\pi/2 - \theta)$ (có trục đối xứng tại $\theta = \pi/4$ và $\theta = 3\pi/4$).

So sánh số trị vận tốc sóng $c_R$ (m/s) giải từ phương trình tán sắc giải tích với kết quả phần tử hữu hạn (FEA) của Cerv và cộng sự (2010):

Vật liệu Góc $\theta$ Kết quả FEA (Cerv et al., 2010) Kết quả giải tích (Luận án)
SE84LV $52^\circ$ 1976.0 1975.2
Fibredux $45^\circ$ 2337.0 2338.2
Boron-Epoxy $49^\circ$ 2320.2 2333.3

Chương 3: Sóng Rayleigh trong bán không gian đàn hồi nằm dưới lớp nước

Chương 3 mở rộng bài toán Bromwich (1898) sang trường hợp tổng quát:

  • Khảo sát bán không gian đàn hồi đẳng hướng không nén được nằm dưới lớp nước lý tưởng không nén được có độ dày $h$ tùy ý, chịu gia tốc trọng trường $g$ tùy ý.
  • Thiết lập phương trình tán sắc chính xác dạng giải tích liên hệ giữa số sóng $k$, vận tốc sóng $c$, độ dày không thứ nguyên $\varepsilon = kh$ và tham số trọng trường $\delta = g/(kc_2^2)$.
  • Xác định miền tồn tại của sóng Rayleigh trong bài toán ghép chất lỏng - chất rắn.
  • Dẫn xuất các phương trình tán sắc xấp xỉ và công thức vận tốc xấp xỉ trong 3 trường hợp:
    1. Khi cả hai đại lượng $\delta$ và $\varepsilon$ đều nhỏ.
    2. Khi chỉ có $\varepsilon$ nhỏ còn $\delta$ tùy ý.
    3. Các công thức xấp xỉ toàn cục trên toàn miền biến thiên của tham số.

Chương 4: Sóng Rayleigh trong bán không gian đàn hồi phủ lớp mỏng

Chương này nghiên cứu phương trình tán sắc xấp xỉ bậc ba bằng phương pháp điều kiện biên hiệu dụng:

  • Bán không gian đàn hồi trực hướng nén được phủ lớp mỏng trực hướng nén được: Thiết lập điều kiện biên hiệu dụng bậc ba tại mặt phân chia $x_3 = 0$. Rút ra phương trình tán sắc xấp xỉ bậc ba theo $\varepsilon = kh$ và suy biến về trường hợp đẳng hướng cũng như công thức vận tốc xấp xỉ bậc hai.
  • Bán không gian đàn hồi trực hướng không nén được phủ lớp mỏng: 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 ba tương ứng.
  • Bán không gian đàn hồi có ứng suất trước phủ lớp mỏng có ứng suất trước: Dẫn xuất phương trình tán sắc xấp xỉ bậc ba cho môi trường biến dạng trước tùy ý; khảo sát các trường hợp riêng như môi trường không có ứng suất trước hoặc môi trường có biến dạng trước đẳng hướng trong mặt phẳng.

Chương 5: Sóng trong cấu trúc mỏng tuần hoàn

Chương 5 mở rộng lý thuyết của Norris và Santosa (1992):

  • Sóng SH trong môi trường vô hạn phân lớp tuần hoàn, các lớp đều mỏng: Mỗi chu kỳ gồm $N$ lớp vật liệu đàn hồi đẳng hướng nén được. Sử dụng khai triển tiệm cận chuỗi lũy thừa theo tham số nhỏ $\varepsilon = kh \ll 1$: $$c^2 = \Omega_1 + (kh)^2 \Omega_3 + (kh)^4 \Omega_5 + \dots$$ Luận án thiết lập công thức tường minh tính $\Omega_1, \Omega_3$ và xây dựng hệ thức truy hồi tổng quát để tính toán các hệ số bậc cao $\Omega_{2n+1}$ $\forall n \ge 2$.
  • Sóng Lamb trong môi trường vô hạn phân lớp tuần hoàn không nén được có biến dạng trước: Thiết lập phương trình tán sắc dạng tiệm cận, công thức tính $\Omega_k$ và công thức truy hồi xác định $\Omega_{2n+1}$ cho sóng Lamb hai thành phần.

Kết quả và những đóng góp mới

Các đóng góp mới của luận án bao gồm:

  1. Tìm ra 2 công thức chính xác và 2 công thức xấp xỉ bậc hai có độ chính xác cao cho vận tốc sóng Rayleigh trong bán không gian đàn hồi đẳng hướng không nén được chịu tác dụng của trọng trường. Chứng minh điều kiện cần và đủ để sóng Rayleigh tồn tại và duy nhất là $0 \le \delta < 1$.
  2. Thiết lập phương trình tán sắc dạng tường minh thống nhất cho sóng chính Rayleigh trong bản mỏng trực hướng bán vô hạn áp dụng cho mọi loại vật liệu trực hướng; dẫn xuất 1 công thức vận tốc chính xác và 2 công thức xấp xỉ.
  3. Rút ra phương trình tán sắc dạng tường minh cho sóng không chính Rayleigh trong bản mỏng trực hướng bằng phương pháp tích phân đầu; chứng minh toán học tính chất đối xứng của đường cong vận tốc sóng theo góc chỉ phương $\theta$.
  4. Giải quyết trọn vẹn bài toán sóng Rayleigh trong bán không gian đàn hồi không nén được nằm dưới lớp nước với độ dày lớp và trọng trường tùy ý; đưa ra điều kiện tồn tại sóng và các công thức vận tốc xấp xỉ.
  5. Xây dựng điều kiện biên hiệu dụng bậc ba và thiết lập các phương trình tán sắc xấp xỉ bậc ba cho sóng Rayleigh trong bán không gian trực hướng (nén được, không nén được) và bán không gian có ứng suất trước phủ lớp mỏng.
  6. Mở rộng công trình của Norris và Santosa (1992), thiết lập các công thức truy hồi xác định các hệ số tiệm cận bậc cao $\Omega_{2n+1}$ ($n \ge 2$) cho sóng SH và sóng Lamb trong cấu trúc phân lớp tuần hoàn có ứng suất trước.

Hạn chế và hướng nghiên cứu tiếp

  • Hạn chế: Các bài toán chủ yếu khảo sát trên mô hình vật liệu đàn hồi tuyến tính lý tưởng, chưa xét đến tính tiêu tán năng lượng (nhớt đàn hồi), hiệu ứng nhiệt đàn hồi hoặc tương tác điện - từ đàn hồi (piezoelectric/piezomagnetic). Đối với cấu trúc phủ lớp mỏng, các nghiệm tiệm cận chỉ áp dụng trong miền bước sóng dài so với độ dày lớp ($kh \ll 1$).
  • Hướng nghiên cứu tiếp:
    • Mở rộng phương pháp tích phân đầu và lý thuyết phương trình bậc ba cho môi trường có tính dị hướng cao hơn (monoclinic, triclinic) chịu tác dụng đồng thời của trường nhiệt độ và điện từ trường.
    • Nghiên cứu bài toán ngược xác định đồng thời độ dày và các hằng số đàn hồi của lớp phủ từ dữ liệu đo phổ tán sắc thực nghiệm.
    • Phát triển mô hình sóng trong cấu trúc phân lớp tuần hoàn có xét đến liên kết không hoàn hảo (imperfect interface) giữa các lớp.

Giá trị tham khảo

  • Đối với nghiên cứu sinh và giảng viên ngành Cơ học vật rắn, Toán ứng dụng, Địa vật lý: Luận án cung cấp hệ thống phương pháp giải tích (phương pháp tích phân đầu, phương pháp phương trình bậc ba, xấp xỉ bình phương tối thiểu, điều kiện biên hiệu dụng). Chương 2 và Chương 4 là tài liệu tham khảo cho việc giải các bài toán sóng mặt trong môi trường dị hướng và cấu trúc đa lớp.
  • Đối với kỹ sư kiểm tra không phá hủy (NDT) và khoa học vật liệu: Các phương trình tán sắc tường minh và công thức vận tốc xấp xỉ cung cấp công cụ tính toán để trích xuất các thông số cơ học, kiểm tra độ bám dính của màng mỏng hoặc đo ứng suất dư trong vật liệu composite cốt sợi mà không cần giải các hệ phương trình vi phân phức tạp.

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

1. Điều kiện tồn tại và tính duy nhất của sóng Rayleigh trong bán không gian đẳng hướng không nén được chịu trọng trường là gì?

Sóng Rayleigh tồn tại và duy nhất khi và chỉ khi tham số trọng trường $\delta = g/(kc_2^2)$ thỏa mãn $0 \le \delta < 1$. Khi $\delta \ge 1$, phương trình tán sắc không có nghiệm thực trong khoảng $(0, 1)$. Khi $0 \le \delta < 1$, vận tốc không thứ nguyên bình phương $x_r(\delta) = c^2/c_2^2$ là hàm đơn điệu tăng thực sự từ $x_0 \approx 0.91262$ (khi $\delta = 0$) đến $1$ (khi $\delta \to 1$).

2. Phương pháp tích phân đầu giải quyết bài toán sóng không chính Rayleigh như thế nào?

Trong bản mỏng trực hướng khi góc truyền sóng $\theta \ne 0$, phương trình vi phân liên hệ chuyển dịch và ứng suất dẫn đến phương trình đặc trưng bậc bốn đầy đủ, khiến phương pháp thông thường không rút ra được nghiệm tường minh. Phương pháp tích phân đầu chuyển bài toán sang hệ vi phân cấp hai đối với biên độ ứng suất $\boldsymbol{\Sigma}$, sử dụng tích vô hướng và tính chất phản đối xứng của các ma trận $D, E, F$ để lập hệ phương trình tuyến tính thuần nhất ba phương trình ba ẩn số. Cho định thức cấp ba này bằng không sẽ thu được phương trình tán sắc dạng tường minh bậc bốn đối với $\rho c^2$ mà không phụ thuộc vào nghiệm của phương trình đặc trưng.

3. Tính chất đối xứng của vận tốc sóng theo góc định hướng sợi $\theta$ được chứng minh ra sao?

Phương trình tán sắc của sóng không chính Rayleigh thỏa mãn đẳng thức hàm $F(X, \theta) = F(X, \pi - \theta)$, suy ra đường cong vận tốc $c_R(\theta)$ luôn nhận đường $\theta = \pi/2$ làm trục đối xứng trong khoảng $[0, \pi]$. Trong trường hợp đặc biệt khi $B_{11} = B_{22}$ (tức $E_1 = E_2$), hàm thỏa mãn thêm $F(X, \theta) = F(X, \pi/2 - \theta)$, tạo thêm hai trục đối xứng tại $\theta = \pi/4$ và $\theta = 3\pi/4$.

4. Bậc xấp xỉ của phương trình tán sắc cho bán không gian phủ lớp mỏng trong luận án đạt đến bậc mấy?

Luận án sử dụng phương pháp điều kiện biên hiệu dụng để thiết lập phương trình tán sắc xấp xỉ đến bậc ba theo độ dày không thứ nguyên $\varepsilon = kh$, nâng cao bậc chính xác so với các xấp xỉ bậc một của Trần Thanh Tuấn (2007) và xấp xỉ bậc hai của Steigmann (2007).

5. Kết quả mở rộng nghiên cứu của Norris và Santosa (1992) trong Chương 5 là gì?

Norris và Santosa (1992) nghiên cứu sóng SH trong môi trường phân lớp tuần hoàn $N$ lớp nhưng chỉ xác định hai hệ số đầu $\Omega_1, \Omega_3$ trong chuỗi tiệm cận $c^2 = \Omega_1 + (kh)^2 \Omega_3 + \dots$. Luận án đã thiết lập công thức truy hồi tổng quát cho phép tính toán tất cả các hệ số bậc cao $\Omega_{2n+1}$ ($\forall n \ge 2$) cho sóng SH, đồng thời mở rộng thành công công thức truy hồi này cho bài toán sóng Lamb trong môi trường phân lớp tuần hoàn không nén được có ứng suất trước.

Kết luận

Luận án "Sóng mặt và sóng trong các cấu trúc mỏng" của tác giả Nguyễn Thị Khánh Linh là công trình nghiên cứu giải tích trong lĩnh vực Cơ học vật rắn. Luận án đã giải quyết trọn vẹn các bài toán về phương trình tán sắc và công thức vận tốc sóng mặt Rayleigh trong bán không gian chịu tác dụng của trọng trường, bản mỏng trực hướng và bán không gian nằm dưới lớp chất lỏng. Đồng thời, công trình đã nâng bậc xấp xỉ điều kiện biên hiệu dụng lên bậc ba cho bán không gian phủ lớp mỏng và thiết lập công thức truy hồi bậc cao cho sóng trong cấu trúc tuần hoàn. Các kết quả giải tích và bảng số liệu trong luận án cung cấp cơ sở lý thuyết phục vụ tính toán âm học, địa chấn học và kỹ thuật đánh giá không phá hủy vật liệu.