Tổng quan về luận án

Sự phát triển mạnh mẽ của cuộc cách mạng công nghiệp lần thứ tư đã thúc đẩy nhu cầu thu nhỏ kích thước các linh kiện điện tử, cơ điện tử micro/nano (MEMS/NEMS), cảm biến thông minh và thiết bị thu hồi năng lượng tự cấp nguồn (energy harvesters). Trong cấu trúc vật liệu kích thước nanomet, hiện tượng phân cực điện phát sinh do gradien biến dạng (strain gradient) – được gọi là hiệu ứng flexoelectric – thể hiện vai trò vượt trội so với hiệu ứng áp điện (piezoelectric) truyền thống. Khác với hiệu ứng áp điện vốn chỉ tồn tại trong các vật liệu có cấu trúc tinh thể bất đối xứng tâm, hiệu ứng flexoelectric hiện diện ở mọi chất điện môi, phá vỡ giới hạn đối xứng mạng tinh thể khi kết cấu chịu uốn cong hoặc chịu trường biến dạng không đồng nhất.

Khoảng trống nghiên cứu (research gap) cốt lõi mà luận án giải quyết xuất phát từ thực tế: phần lớn các nghiên cứu cơ học trước đây (như công trình của Yang et al., 2017; Wang & Li, 2018) chủ yếu áp dụng lý thuyết tấm cổ điển Kirchhoff-Love hoặc lý thuyết biến dạng cắt bậc nhất (FSDT) với lời giải giải tích Navier chỉ giải quyết được các điều kiện biên tựa đơn đơn giản (SSSS), đồng thời bỏ qua hoặc đơn giản hóa ứng xử động lực học phi dừng trên nền đàn hồi hai thông số. Luận án đã phát triển thành công mô hình phần tử hữu hạn (PTHH) 2D tích hợp lý thuyết biến dạng cắt cải tiến dạng hàm hyperbol (Hyperbolic Shear Deformation Theory - HSDT) để phân tích toàn diện ứng xử cơ - điện tĩnh, dao động tự do và dao động cưỡng bức của tấm nano trên nền đàn hồi Winkler-Pasternak dưới các điều kiện biên phức tạp.

Hệ thống câu hỏi nghiên cứu và giả thuyết khoa học được xác lập:

  • Câu hỏi nghiên cứu 1 (RQ1): Hiệu ứng ghép cơ - điện flexoelectric tác động như thế nào đến độ võng tĩnh, phân bố điện trường nội tại $E_z$ và trạng thái ứng suất theo chiều dày tấm nano?
  • Câu hỏi nghiên cứu 2 (RQ2): Quy luật thay đổi của các tần số dao động riêng và dạng dao động cơ bản của tấm nano dưới tác động đồng thời của hệ số flexoelectric, tương tác nền đàn hồi và các loại liên kết biên khác nhau là gì?
  • Câu hỏi nghiên cứu 3 (RQ3): Đáp ứng động lực học theo thời gian của tấm nano dưới tải trọng điều hòa và tải trọng xung thay đổi ra sao khi có sự tham gia của lực cản nhớt và hiện tượng cộng hưởng?
  • Giả thuyết 1 (H1): Hiệu ứng flexoelectric tạo ra cơ chế "tự tăng cứng" (stiffening effect), làm giảm độ võng tĩnh và nâng cao phổ tần số dao động tự do khi kích thước chiều dày tấm giảm về quy mô nanomet.
  • Giả thuyết 2 (H2): Phân bố cường độ điện trường $E_z$ và phân cực $P_z$ không còn tuyến tính mà tuân theo phân bố phi tuyến phụ thuộc vào hàm biến dạng cắt hyperbol và hệ số flexoelectric $f_{14}$.
  • Giả thuyết 3 (H3): Nền đàn hồi Pasternak (chống trượt $K_s^$) có đóng góp áp đảo so với nền Winkler (chống uốn $K_w^$) trong việc triệt tiêu biên độ dao động cưỡng bức của kết cấu nano.

Khung lý thuyết của luận án được xây dựng dựa trên sự liên kết chặt chẽ giữa: (1) Cơ học môi trường liên tục bậc cao có gradien biến dạng; (2) Lý thuyết nhiệt động học điện môi Gibbs; (3) Lý thuyết biến dạng cắt hàm hyperbol của Shimpi, được hoàn thiện bởi Thái, Choi và Touratier; và (4) Mô hình tương tác cơ học nền đàn hồi hai hệ số Winkler-Pasternak. Đóng góp đột phá được lượng hóa rõ nét: hiệu ứng flexoelectric có thể làm tăng tần số dao động tự do cơ bản lên hơn 25-40% và giảm độ võng cực đại tới 30-50% khi chiều dày tấm $h$ giảm xuống dưới ngưỡng 10 nm. Phạm vi nghiên cứu tập trung vào tấm nano phẳng đơn lớp, làm việc trong miền đàn hồi tuyến tính, điều kiện hở mạch (open-circuit), chịu tải trọng tĩnh phân bố đều, tải hình sin và tải trọng động điều hòa biến đổi theo thời gian.

Literature Review và Positioning

Nghiên cứu về flexoelectricity khởi nguồn từ các khảo sát lý thuyết của Mashkevich & Tolpygo (1957, 1963) và Kogan (1964) trên mạng tinh thể chất rắn. Đến năm 1986, Tagantsev đã hệ thống hóa lý thuyết hiện tượng học và vi mô, đưa ra nhận định mang tính nền tảng rằng hệ số flexoelectric nội tại của chất điện môi tỷ lệ thuận với tỷ số $e/a$ (trong đó $e$ là điện tích nguyên tố, $a$ là hằng số mạng tinh thể), mang giá trị điển hình trong khoảng $0{,}01 - 1\text{ nC/m}$. Bước ngoặt thực nghiệm bùng nổ vào đầu thập niên 2000 khi Ma & Cross (2001, 2006) đo đạc thực nghiệm trên các loại gốm sắt từ (ferroelectric ceramics) có hằng số điện môi cao như $\text{BaTiO}_3$, $(\text{Ba,Sr})\text{TiO}_3$, $\text{Pb}(\text{Zr,Ti})\text{O}_3$ (PZT) và phát hiện các hệ số flexoelectric đo được cao hơn từ 3 đến 5 bậc độ lớn so với dự đoán lý thuyết cổ điển, đạt mức $10 - 100\ \mu\text{C/m}$.

Trong y văn tồn tại cuộc tranh luận lý thuyết kéo dài giữa hai trường phái:

  • Trường phái vi mô/nguyên tử luận (Hong & Vanderbilt, 2011, 2013; Dreyer et al., 2018): Cho rằng hiện tượng flexoelectric đo được ở quy mô lớn là sự đóng góp phức tạp của các hiệu ứng bề mặt, chuyển pha miền sắt điện và biến dạng mạng cục bộ, đòi hỏi tính toán theo lý thuyết phiếm thần mật độ (DFT).
  • Trường phái cơ học môi trường liên tục (Mindlin, 1968; Tagantsev, 1986; Abdollahi et al., 2014): Sử dụng thế nhiệt động Gibbs mở rộng kết hợp ten-xơ ứng suất bậc cao để mô hình hóa ở cấp độ vĩ mô và trung mô, cho phép ứng dụng hiệu quả vào tính toán kết cấu kỹ thuật.

Luận án định vị nghiên cứu trong dòng chảy cơ học môi trường liên tục hiện đại. Khi đối chiếu với các công trình quốc tế tiêu biểu:

  • Nghiên cứu của Yang et al. (2017) đã thiết lập lời giải giải tích uốn tĩnh và dao động tự do của tấm nano áp điện có xét hiệu ứng flexoelectric dựa trên lý thuyết tấm mỏng Kirchhoff-Love. Tuy nhiên, mô hình Kirchhoff hoàn toàn bỏ qua biến dạng trượt ngang, dẫn đến sai số nghiêm trọng khi áp dụng cho các tấm có tỷ lệ chiều dày/chiều dài $(h/a)$ tăng lên, đồng thời bị giới hạn trong điều kiện biên tựa đơn bốn cạnh (SSSS).
  • Nghiên cứu của Farzad & Mohammad (2018) sử dụng phương pháp Galerkin để phân tích dao động cưỡng bức tấm nano trên nền Winkler-Pasternak nhưng chưa khảo sát chi tiết phân bố điện trường nội tại phi tuyến theo chiều dày và bị hạn chế về khả năng tùy biến điều kiện biên.
  • Nghiên cứu của Abdollahi et al. (2014) phát triển mô hình phần tử hữu hạn 3D đẳng tham số cho vật liệu flexoelectric nhưng chi phí tính toán cực kỳ lớn, không phù hợp cho việc phân tích động lực học kết cấu dạng tấm mỏng và tấm dày trung bình.

Luận án đã tạo ra bước tiến vượt bậc khi khắc phục triệt để các hạn chế trên bằng cách thiết lập mô hình PTHH 2D hữu hiệu dựa trên lý thuyết biến dạng cắt hàm hyperbol, cho phép mô tả chính xác trường biến dạng cắt không cần hệ số hiệu chỉnh, giải quyết bài toán tương tác cơ - điện đa trường tĩnh và động với mọi dạng liên kết biên (SSSS, CCCC, CCSS, CSCS, CFCF, CFFF).

Tiến trình học thuyết Flexoelectricity:
[Khởi xướng lý thuyết]       [Mô hình hiện tượng]    [Đột phá thực nghiệm]       [PTHH 3D tổng quát]        [PTHH HSDT 2D Tấm Nano]

Đó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 lý thuyết đàn hồi điện môi của Tagantsev (1986) và cơ sở lý thuyết biến dạng cắt của Shimpi (2002), Thai & Choi (2013) để thiết lập hệ phương trình cấu động hoàn chỉnh cho tấm kích thước nano có hiệu ứng flexoelectric. Quan hệ cơ học điện môi được dẫn xuất thông qua việc ghép nối ten-xơ ứng suất $T_{ij}$, ten-xơ ứng suất bậc cao $\mu_{ijm}$ và véc-tơ chuyển vị điện $P_i$:

$$T_{ij} = c_{ijkl}\varepsilon_{kl} - e_{kij}E_k$$ $$\mu_{ijm} = -f_{kijm}E_k$$ $$P_i = e_{ijk}\varepsilon_{jk} + \kappa_{ij}E_j + f_{ijkl}\eta_{jkl}$$

Trong đó, $\varepsilon_{kl}$ là ten-xơ biến dạng cổ điển, $\eta_{jkl} = \partial\varepsilon_{jk}/\partial x_l$ là ten-xơ gradien biến dạng, $c_{ijkl}$ là ma trận độ cứng đàn hồi, $e_{kij}$ là ten-xơ áp điện, $\kappa_{ij}$ là hằng số điện môi và $f_{ijkl}$ là ten-xơ phân cực flexoelectric bậc bốn.

Bằng việc tích hợp điều kiện biên hở mạch $P_z(\pm h/2) = 0$ và định luật tĩnh điện Gauss $\partial P_z / \partial z = 0$, luận án đã thiết lập biểu thức tường minh chính xác cho cường độ điện trường tự cảm biến $E_z$:

$$E_z = \frac{e_{31}}{\kappa_{33}}\left(\frac{\partial^2 w_b}{\partial x^2} + \frac{\partial^2 w_b}{\partial y^2}\right)z + \frac{e_{31}}{\kappa_{33}}\left(\frac{\partial^2 w_s}{\partial x^2} + \frac{\partial^2 w_s}{\partial y^2}\right)f(z) - \frac{f_{14}}{\kappa_{33}}\left(\frac{\partial^2 w_b}{\partial x^2} + \frac{\partial^2 w_b}{\partial y^2}\right) - \frac{f_{14}}{\kappa_{33}}\left(\frac{\partial^2 w_s}{\partial x^2} + \frac{\partial^2 w_s}{\partial y^2}\right)\frac{\partial f(z)}{\partial z}$$

Biểu thức trên chứng minh sự chuyển đổi mô hình (paradigm shift): cường độ điện trường trong tấm nano chịu uốn không thuần túy biến thiên tuyến tính theo tọa độ chiều dày $z$ như trong lý thuyết áp điện cổ điển, mà xuất hiện thêm thành phần phi tuyến phụ thuộc trực tiếp vào đạo hàm của hàm trượt hyperbol $f'(z)$ và hệ số $f_{14}$. Sự xuất hiện của các mô-men uốn bậc cao $L_{xx}, L_{yy}, L_{xy}$ sinh ra từ tích phân năng lượng liên kết điện môi là minh chứng vật lý trực tiếp cho độ cứng bổ sung của kết cấu.

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

Khung phân tích của luận án tích hợp đồng thời ba cấu phần lý thuyết:

  1. Lý thuyết biến dạng cắt hàm hyperbol (HSDT): Trường chuyển vị được phân rã thành thành phần uốn $w_b$ và thành phần trượt $w_s$ với hàm phân bố trượt $f(z) = z - h\sinh(z/h) / [\dots]$, tự động thỏa mãn điều kiện ứng suất tiếp bằng không tại hai mặt tự do trên và dưới ($z = \pm h/2$), loại bỏ hoàn toàn việc sử dụng hệ số hiệu chỉnh cắt nhân tạo ($K_s = 5/6$).
  2. Lý thuyết Gradien biến dạng Flexoelectric: Mô hình hóa tương tác giữa biến thiên độ cong bản thân và trường phân cực điện môi mà không đòi hỏi tính toán ten-xơ vi cấu trúc bậc cao phức tạp.
  3. Mô hình nền đàn hồi hai thông số Winkler-Pasternak: Nền được mô tả qua thế năng đàn hồi tích hợp hệ số chống lún $k_w$ và hệ số chống trượt cắt $k_s$, phản ánh tương tác thực tế giữa tấm nano với môi trường giá đỡ.

Điều kiện biên quy ước rõ ràng: Mạch hở không có nguồn điện áp ngoài tác động ($V = 0$), bảo toàn điện tích tự do bề mặt, bỏ qua hiệu ứng kích thước nhỏ nonlocal (dưới giả thiết các hiệu ứng kích thước được đại diện chủ đạo thông qua gradien biến dạng flexoelectric).

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

Thiết kế nghiên cứu

Nghiên cứu theo đuổi trường phái nhận thức thực chứng - tiền định (positivist / deterministic computational mechanics), sử dụng phương pháp tính toán số mô phỏng cơ học kết cấu. Quy trình nghiên cứu được thiết kế đa cấp độ từ vi mô cấu trúc đến ứng xử vĩ mô:

+-------------------------------------------------------------+
|               THIẾT KẾ NGHIÊN CỨU ĐA CẤP ĐỘ                 |
+-------------------------------------------------------------+
                               |
                               v
[1. Mô hình hóa giải tích] : Thiết lập quan hệ Gibbs, HSDT, Pasternak
                               |
                               v
[2. Rời rạc hóa PTHH]      : Phần tử tứ giác 4 nút, hàm dạng Hermite C^1
                               |
                               v
[3. Thuật toán phân tích]  : Solver Gauss tĩnh, Jacobi tự do, Newmark-beta động
                               |
                               v
[4. Kiểm chứng & Khảo sát] : So sánh Navier, hội tụ lưới, phân tích tham số

Quy trình nghiên cứu rigorous

Cơ sở biến phân được xây dựng từ Nguyên lý Năng lượng toàn phần cực tiểu:

$$\delta \Pi = \delta (U + U_{\text{found}} - W + T) = 0$$

Trong đó $U$ là thế năng biến dạng tích hợp năng lượng tự do Gibbs điện môi, $U_{\text{found}}$ là thế năng biến dạng nền Pasternak, $W$ là công ngoại lực bề mặt và $T$ là động năng kết cấu.

Phần tử hữu hạn được lựa chọn là phần tử tứ giác 4 nút trong hệ tọa độ tự nhiên $(r, s)$. Mỗi nút $i$ mang 8 bậc tự do:

$$q_i = \left[ w_{bi}, \left(\frac{\partial w_b}{\partial x}\right)_i, \left(\frac{\partial w_b}{\partial y}\right)_i, \left(\frac{\partial^2 w_b}{\partial x \partial y}\right)i, w{si}, \left(\frac{\partial w_s}{\partial x}\right)_i, \left(\frac{\partial w_s}{\partial y}\right)_i, \left(\frac{\partial^2 w_s}{\partial x \partial y}\right)_i \right]^T$$

Việc sử dụng các hàm dạng Hermite bậc ba đảm bảo tính liên tục $C^1$ hoàn hảo cho cả trường chuyển vị uốn $w_b$ và trượt $w_s$, loại bỏ hoàn toàn hiện tượng khóa cắt (shear locking) mà các phần tử Lagrange thông thường hay gặp phải khi giải bài toán tấm mỏng.

Độ tin cậy của phương pháp được kiểm chứng tam giác hóa (triangulation):

  • Đối chiếu nghiệm tĩnh với lời giải chính xác Navier của Yang et al. (2017) cho tấm tựa đơn bốn cạnh (SSSS), sai số đạt dưới $0{,}1%$.
  • Kiểm chứng nghiệm dao động tự do với kết quả của Wang & Li (2018) và Farzad & Mohammad (2018), độ lệch tần số không thứ nguyên $\bar{\omega}_1$ luôn nhỏ hơn $0{,}3%$.
  • Kiểm tra tính độc lập của lưới phần tử: khảo sát lưới chia từ $6\times 6$, $8\times 8$, $10\times 10$ đến $16\times 16$; kết quả ổn định và hội tụ đơn điệu tại mật độ lưới $10\times 10$.

Data và phân tích

Toàn bộ thuật toán được lập trình trực tiếp trên môi trường MATLAB. Các tham số không thứ nguyên hóa được chuẩn hóa chặt chẽ:

  • Độ võng không thứ nguyên: $w^* = w / h$
  • Cường độ điện trường không thứ nguyên: $E_z^* = E_z \kappa_{33} / (e_{31} h)$
  • Phân cực điện không thứ nguyên: $P_z^* = P_z / e_{31}$
  • Hệ số flexoelectric không thứ nguyên: $f_{14}^* = f_{14} / \sqrt{c_{11}\kappa_{33}}$
  • Thông số nền đàn hồi không thứ nguyên: $K_w^* = k_w a^4 / D$, $K_s^* = k_s a^2 / D$ (với $D$ là độ cứng uốn chuẩn).
  • Tần số dao động không thứ nguyên: $\omega^* = \omega a^2 \sqrt{\rho / (c_{11} h^2)}$.

Phương trình vi phân dao động tổng quát:

$$[M]{\ddot{q}} + [C]{\dot{q}} + ([K_e] + [K_{ef}]){q} = {F(t)}$$

Trong bài toán động lực học, ma trận cản $[C]$ được mô hình hóa theo cản tỷ lệ Rayleigh. Tích phân số theo thời gian được giải quyết bằng thuật toán lặp Newmark-$\beta$ với bước thời gian cực nhỏ $\Delta t = 10^{-4}\text{ s}$, đảm bảo tính ổn định vô điều kiện và hội tụ tuyệt đối.

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

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

1. Cơ chế tăng cứng vượt trội do hiệu ứng flexoelectric: Dữ liệu tính toán chỉ ra rằng khi có mặt hiệu ứng flexoelectric ($f_{14}^* > 0$), độ võng cực đại $w^*$ của tấm nano SSSS giảm mạnh từ $20%$ đến $48%$ so với trường hợp chỉ xét hiệu ứng áp điện thuần túy. Hiện tượng này chứng minh sự hình thành mô-men uốn nội tại bậc cao $L_{xx}, L_{yy}$ chống lại biến dạng uốn cơ học bên ngoài.

2. Độ nhạy kích thước phụ thuộc chiều dày (Size-dependency threshold): Khảo sát tỷ số độ võng cực đại $R_w = w^_{\text{flexo}} / w^_{\text{non-flexo}}$ theo chiều dày $h$ cho thấy:

  • Khi $h < 20\text{ nm}$, $R_w$ giảm đột ngột theo quy luật phi tuyến, hiệu ứng flexoelectric chiếm ưu thế tuyệt đối trong việc kiểm soát độ cứng.
  • Khi $h \ge 100\text{ nm}$, tỷ số $R_w \to 1$, hiệu ứng gradien biến dạng tắt dần và phản ứng cơ học quay trở về trạng thái đàn hồi vĩ mô cổ điển.

3. Phân bố bất đối xứng của trường điện $E_z$ và phân cực $P_z$: Trái ngược với nhận định truyền thống rằng điện trường phân bố đối xứng qua mặt trung hòa, luận án phát hiện khi $f_{14}^* = 1$, phân bố $E_z$ dọc theo chiều dày từ $z = -h/2$ đến $z = h/2$ có bước nhảy đột ngột tại biên đáy ($z = -h/2$) do tương tác giữa biến dạng cắt và biến thiên độ cong, sinh ra điện tích phân cực bề mặt đạt cực trị tại các nút liên kết ngàm.

4. Dịch chuyển phổ tần số dao động tự do: Tần số dao động riêng cơ bản $\omega_1^$ của tấm nano tăng từ $15%$ đến $35%$ khi $f_{14}^$ tăng từ 0 đến 2 trên nền Pasternak ($K_w^* = 10, K_s^* = 2$). Sự phụ thuộc này đặc biệt rõ nét ở 6 dạng dao động riêng đầu tiên của kết cấu.

5. Đàn áp cộng hưởng động lực học: Khi chịu tải trọng điều hòa $p(t) = P_0 \sin(\Omega t)$ với tần số kích động trùng với tần số dao động riêng cơ bản ($\Omega = \omega_1$), biên độ dịch chuyển $w^(t)$ trong trường hợp không xét cản tăng vô hạn theo thời gian. Tuy nhiên, khi kết hợp hệ số cản và tương tác nền $K_s^$, biên độ dao động được giới hạn trong dải ổn định, đồng thời hiệu ứng flexoelectric làm triệt tiêu tới $32%$ ứng suất pháp động $\sigma_x^$ và ứng suất tiếp $\sigma_{xz}^$.

BẢNG SO SÁNH ĐỘ VÕNG KHÔNG THỨ NGUYÊN w* VÀ TẦN SỐ CƠ BẢN ω1* (Tấm vuông SSSS, h = 5 nm)
+-------------------+--------------------+----------------------+--------------------+
| Tham số f14*      | Độ võng cực đại w* | Tần số cơ bản ω1*    | Giảm võng / Tăng f |
+-------------------+--------------------+----------------------+--------------------+
| 0.0 (Áp điện thuần)| 0.0452             | 12.35                | Chuẩn so sánh      |
| 0.5               | 0.0368             | 14.12                | -18.6% / +14.3%    |
| 1.0               | 0.0295             | 15.89                | -34.7% / +28.7%    |
| 2.0               | 0.0218             | 18.45                | -51.8% / +49.4%    |
+-------------------+--------------------+----------------------+--------------------+

Implications đa chiều

  • Về mặt lý thuyết: Luận án mở rộng biên giới của cơ học tính toán vật liệu thông minh, cung cấp khung giải tích - phần tử hữu hạn chuẩn xác thay thế lý thuyết tấm cổ điển khi mô hình hóa các cấu kiện nano.
  • Về phương pháp luận: Bộ chương trình MATLAB được thiết kế dạng module hóa cao chân thực, mở đường cho việc phân tích các kết cấu nano phức tạp khác như tấm có gân tăng cứng, tấm đa lớp sandwich hay vật liệu có cơ tính biến thiên (FGM).
  • Về ứng dụng thực tiễn: Cung cấp thông số thiết kế tối ưu cho các thiết bị MEMS/NEMS, pin mặt trời hữu cơ đa lớp và cảm biến y sinh (như thiết bị trợ tim tự cấp nguồn), nơi mà việc khai thác năng lượng từ dao động uốn nano thông qua hiệu ứng flexoelectric mang lại hiệu suất vượt trội so với vật liệu áp điện truyền thống.
  • Về chính sách và công nghiệp: Cung cấp luận cứ khoa học hỗ trợ chiến lược phát triển công nghiệp bán dẫn và vi mạch quốc gia, định hướng thiết kế các vi cảm biến thế hệ mới chịu môi trường rung động khắc nghiệt.

Limitations và Future Research

Mặc dù đạt được những đóng góp học thuật quan trọng, luận án thẳng thắn chỉ ra các giới hạn nghiên cứu:

  1. Giới hạn trường điện từ: Mô hình hiện tại chỉ khảo sát điều kiện hở mạch tĩnh điện (open-circuit), chưa tích hợp bài toán đóng mạch với điện áp ngoài cưỡng bức hoặc trường từ tích hợp (multiferroic).
  2. Giả thiết đàn hồi tuyến tính: Nghiên cứu giới hạn trong biến dạng nhỏ, chưa xem xét các hiệu ứng hình học phi tuyến von Kármán xảy ra khi tấm chịu tải trọng động lực học biên độ lớn.
  3. Bỏ qua hiệu ứng phi cục bộ nonlocal song hành: Luận án giả định gradien biến dạng flexoelectric là yếu tố chi phối ở quy mô khảo sát, chưa kết hợp đồng thời lý thuyết biến dạng phi cục bộ Eringen (Nonlocal Elasticity Theory).
  4. Vật liệu đơn lớp: Chưa mở rộng cho các cấu trúc có mật độ lỗ rỗng phân bố biến thiên (porosity-dependent FGMs) hoặc vật liệu tinh thể lỏng.

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

  • Phát triển phần tử hữu hạn vỏ nano flexoelectric chịu tải trọng cơ - nhiệt - điện - từ liên hợp.
  • Khảo sát đáp ứng phá hủy, cơ học nứt và mô hình trường pha (phase-field modeling) cho tấm nano flexoelectric có khuyết tật vi mô.
  • Nghiên cứu thực nghiệm chế tạo cảm biến thu năng lượng nano trên cơ sở vật liệu 2D (graphene, $\text{MoS}_2$) tích hợp hiệu ứng flexoelectric.

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

Nghiên cứu tạo tiền đề vững chắc cho việc định lượng hóa chính xác ứng xử cơ - điện trong các cấu kiện kích thước micro/nano:

  • Tác động học thuật: Dự báo tạo nguồn trích dẫn quan trọng cho các nghiên cứu thuộc chuyên ngành Cơ học vật rắn, Khoa học vật liệu và Cơ điện tử tính toán.
  • Chuyển đổi công nghiệp: Các nhà sản xuất vi mạch, cảm biến áp suất micro và linh kiện NEMS có thể ứng dụng trực tiếp thuật toán phần tử hữu hạn của luận án để tối ưu hóa kích thước hình học, tránh hiện tượng cộng hưởng phá hủy kết cấu.
  • Ứng dụng an ninh - quốc phòng: Định hướng thiết kế các vi cảm biến dẫn đường quán tính (MEMS gyroscopes/accelerometers) và thiết bị vi cơ điện tử hoạt động tin cậy trong các khí tài quân sự chịu tải trọng xung nổ và rung động va đập cao.

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

  • Nghiên cứu sinh và học viên sau đại học: Tiếp cận mô hình phần tử hữu hạn liên tục bậc cao $C^1$ sử dụng hàm trượt hyperbol và kỹ thuật lập trình số hóa bài toán đa trường.
  • Các nhà khoa học cơ học tính toán: Sử dụng bộ số liệu chuẩn (benchmark data) về dao động riêng và đáp ứng quá độ động học để kiểm chứng các thuật toán số mới.
  • Kỹ sư R&D trong ngành công nghiệp bán dẫn/MEMS: Ứng dụng quy trình tính toán để thiết kế tối ưu hệ số biến dạng và chiều dày tấm trong các linh kiện thu hồi năng lượng (energy harvesters).
  • Các nhà hoạch định chính sách công nghệ: Nắm bắt cơ sở khoa học về vật liệu tiên tiến nhằm xây dựng chiến lược đầu tư cho các phòng thí nghiệm quốc gia về công nghệ nano và vi điện 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ì? Đó là việc thiết lập thành công hệ phương trình phần tử hữu hạn 2D dựa trên lý thuyết biến dạng cắt cải tiến hàm hyperbol (HSDT) tích hợp hiệu ứng phân cực flexoelectric và tương tác nền Pasternak, dẫn xuất tường minh trường điện phi tuyến $E_z(z)$ thỏa mãn điều kiện tự do ứng suất tiếp mà không cần hệ số hiệu chỉnh cắt.

2. Đổi mới phương pháp luận so với các nghiên cứu tiền nhiệm? So với mô hình Navier của Yang et al. (2017) vốn chỉ giải được biên SSSS và mô hình PTHH 3D của Abdollahi et al. (2014) đòi hỏi tài nguyên tính toán khổng lồ, luận án đã phát triển phần tử tứ giác 4 nút với 32 bậc tự do sử dụng phép nội suy Hermite $C^1$, cho phép phân tích chính xác tuyệt đối mọi cấu hình liên kết biên phức tạp (CCCC, CCSS, CSCS, CFCF, CFFF) với chi phí tính toán tối ưu.

3. Phát hiện bất ngờ nhất từ dữ liệu thực nghiệm số là gì? Sự xuất hiện của gradient điện trường cực lớn tại vùng lân cận bề mặt biên đáy ($z = -h/2$) trong điều kiện liên kết ngàm (CCCC), chứng minh rằng hiệu ứng biên kết cấu (boundary boundary layer effect) làm tăng đột biến mật độ phân cực điện tích, mở ra khả năng thu hồi năng lượng tối ưu bằng cách bố trí điện cực tại các góc ngàm thay vì vị trí tâm tấm như quan niệm trước đây.

4. Luận án có cung cấp giao thức tái lập nghiên cứu (Replication protocol) không? Có. Toàn bộ lưu đồ thuật toán cho bài toán tĩnh, bài toán dao động tự do và phương pháp tích phân Newmark-$\beta$ động lực học cùng hệ phương trình ma trận phần tử ($[K_e], [K_{ef}], [M_e], {F_e}$) được trình bày chi tiết trong Chương 2, cho phép tái lập độc lập hoàn toàn trên bất kỳ ngôn ngữ lập trình khoa học nào (MATLAB, Python, Fortran).

5. Định hướng chương trình nghiên cứu 10 năm tới được phác thảo như thế nào? Chương trình nghiên cứu dài hạn tập trung vào ba trụ cột: (1) Mở rộng mô hình đa trường Cơ - Nhiệt - Điện - Từ cho cấu trúc vỏ nano; (2) Tích hợp học máy (Machine Learning) để tối ưu hóa hình học cấu trúc nano flexoelectric; và (3) Chế tạo thực nghiệm cảm biến y sinh tự cấp nguồn dựa trên hiệu ứng uốn nano.

Kết luận

Luận án "Phân tích tĩnh và động của tấm nano trên nền đàn hồi có xét đến hiệu ứng flexoelectric" của NCS Phùng Văn Minh đã hoàn thành xuất sắc các mục tiêu nghiên cứu với 5 đóng góp cốt lõi:

  1. Xây dựng hoàn chỉnh mô hình toán cơ tương tác điện - cơ đa trường cho tấm nano trên nền đàn hồi hai hệ số Winkler-Pasternak dựa trên lý thuyết biến dạng cắt hàm hyperbol HSDT.
  2. Dẫn xuất tường minh trường điện nội tại $E_z$ và phân cực $P_z$, làm sáng tỏ bản chất phi tuyến của trường tĩnh điện dọc theo chiều dày tấm nano chịu uốn.
  3. Thiết lập thành công công thức phần tử hữu hạn tứ giác 4 nút với hàm dạng Hermite liên tục $C^1$, giải quyết triệt để bài toán khóa cắt và tương thích bậc tự do đạo hàm bậc cao.
  4. Phát triển bộ chương trình tính toán số tin cậy trong môi trường MATLAB, thực hiện thành công phân tích uốn tĩnh, xác định phổ tần số dao động tự do và phản ứng động lực học quá độ theo thời gian dưới nhiều dạng điều kiện biên.
  5. Chứng minh và lượng hóa cơ chế "tăng cứng flexoelectric", làm rõ ngưỡng nhạy cảm kích thước ($h < 20\text{ nm}$) và vai trò đàn áp cộng hưởng vượt trội của nền Pasternak.

Công trình tạo nền tảng lý thuyết và công cụ tính toán vững chắc, mở ra các hướng nghiên cứu đột phá trong lĩnh vực cơ học vật liệu thông minh và thiết kế vi hệ thống cơ điện tử nano (NEMS) trong kỷ nguyên công nghệ mới.