Tổng quan về luận án
Vật liệu có cơ tính biến thiên (Functionally Graded Materials - FGM) được các nhà khoa học thuộc Đại học Tohoku (Sendai, Nhật Bản) khởi xướng lần đầu tiên vào năm 1984 nhằm phát triển lớp vật liệu cách nhiệt chịu sốc nhiệt cực hạn trong công nghiệp hàng không vũ trụ. Nhờ sự biến thiên liên tục và trơn của tỷ phần thể tích các pha thành phần (điển hình là gốm và kim loại) theo tọa độ không gian, kết cấu FGM khắc phục triệt để các nhược điểm cố hữu của vật liệu composite nhiều lớp truyền thống như hiện tượng tập trung ứng suất tại mặt phân giới và sự tách lớp (delamination). Trong xu thế tối ưu hóa kết cấu chịu tải trọng phức tạp đa chiều, kết cấu dầm và tấm sandwich FGM có cơ tính biến thiên đồng thời theo cả chiều cao ($z$) và chiều dài ($x$)—gọi tắt là vật liệu 2D-FGM—đang trở thành tâm điểm nghiên cứu tiên phong của cơ học vật rắn biến dạng hiện đại.
Tuy nhiên, khoảng trống nghiên cứu (research gap) lớn trong y văn quốc tế hiện nay là phần lớn các công bố mới chỉ dừng lại ở kết cấu FGM biến thiên một chiều theo chiều dày (1D-FGM) hoặc vật liệu hai pha đơn giản dưới tải trọng tĩnh. Các bài toán động lực học phi dừng phức tạp—đặc biệt là dao động của dầm và tấm sandwich 2D-FGM cấu tạo từ hai và ba pha vật liệu chịu tác động của tải trọng di động (bao gồm lực điều hòa di động và khối lượng di động có xét tới lực quán tính, lực Coriolis và lực ly tâm)—vẫn chưa được giải quyết một cách toàn diện bằng các mô hình số có độ chính xác cao. Luận án tiến sĩ của Nghiên cứu sinh Phạm Vũ Nam, 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 và Học viện Khoa học và Công nghệ, đã giải quyết khoảng trống học thuật này thông qua đề tài: "Phân tích dao động dầm, tấm sandwich 2D-FGM hai và ba pha bằng phương pháp phần tử hữu hạn".
Nghiên cứu tập trung giải quyết 4 câu hỏi nghiên cứu và hệ thống giả thuyết khoa học tương ứng:
- Câu hỏi nghiên cứu 1 (RQ1): Quy luật phân bố vật liệu hai chiều ($n_x, n_z$) và tỷ số chiều dày các lớp ảnh hưởng như thế nào đến phổ tần số dao động tự do của dầm và tấm sandwich 2D-FGM hai pha và ba pha?
- Giả thuyết 1 (H1): Sự gia tăng chỉ số gradient theo chiều dài $n_x$ và chiều cao $n_z$ làm suy giảm phi tuyến độ cứng tổng thể của kết cấu, làm dịch chuyển đáng kể các tần số dao động riêng $\omega_i$.
- Câu hỏi nghiên cứu 2 (RQ2): Mức độ sai lệch về đáp ứng động lực học giữa các mô hình vi cơ học đồng nhất hóa vật liệu (Voigt, Mori-Tanaka, Maxwell) là bao nhiêu khi áp dụng cho kết cấu sandwich 2D-FGM?
- Giả thuyết 2 (H2): Mô hình Voigt đánh giá độ cứng hiệu dụng cao hơn so với mô hình Mori-Tanaka và Maxwell, dẫn đến việc ước tính tần số riêng lớn hơn và độ võng động học nhỏ hơn.
- Câu hỏi nghiên cứu 3 (RQ3): Tương tác động lực học giữa tải trọng di động (vận tốc $v$, tần số kích động $\Omega$, tỷ số khối lượng $r_m$) và dầm sandwich 2D-FGM gây ra các hiện tượng đáp ứng chuyển vị và phân bố ứng suất $(\sigma_{xx}, \tau_{xz})$ như thế nào?
- Giả thuyết 3 (H3): Tải trọng khối lượng di động kích hoạt lực Coriolis và lực ly tâm làm thay đổi vị trí và biên độ độ võng cực đại so với mô hình tải trọng lực tập trung đơn thuần.
- Câu hỏi nghiên cứu 4 (RQ4): Làm thế nào để xây dựng các phần tử hữu hạn chuyên dụng (phần tử dầm làm giàu và phần tử tấm Q9 nội suy liên kết) nhằm triệt tiêu hoàn toàn hiện tượng nghẽn trượt (shear locking) và tăng tốc độ hội tụ?
- Giả thuyết 4 (H4): Việc tích hợp các trường nội suy làm giàu lượng giác và phép nội suy liên kết (tied/mixed interpolation) cho phần tử tấm Q9 sẽ đảm bảo tính ổn định và độ chính xác vượt bậc cho kết cấu từ dải mỏng đến rất dày ($a/h \in [5, 100]$).
Khung lý thuyết của công trình được định hình vững chắc dựa trên nguyên lý biến phân Hamilton cho cơ học liên tục đàn hồi tuyến tính $\int_{t_1}^{t_2} \delta [T - (U + V)] dt = 0$, Lý thuyết biến dạng trượt bậc ba cải tiến (Higher-order Shear Deformation Theory - TSDT) của Bickford và Reddy, Lý thuyết biến dạng trượt lượng giác (Trigonometric Shear Deformation Theory), Lý thuyết tấm Mindlin và các sơ đồ đồng nhất hóa vi cơ học vi mô. Phạm vi nghiên cứu bao quát các kết cấu dầm sandwich 2D-FGM hai pha (lõi thuần gốm, hai lớp ngoài 2D-FGM), dầm sandwich 2D-FGM ba pha (lõi và vỏ kết hợp ba pha gốm/kim loại), và tấm chữ nhật sandwich 2D-FGM ba pha với các điều kiện biên tựa giản đơn (SSSS), kẹp chặt (CCCC), và kết hợp (SCSC), tỷ số kích thước hình học $L/h, a/h \in [5, 20]$, chỉ số gradient $n_x, n_z \in [0, 10]$, và dải vận tốc tải trọng $v$ lên đến 100 m/s.
Literature Review và Positioning
Nghiên cứu về cơ học kết cấu FGM đã trải qua bốn thập kỷ phát triển với ba dòng chảy học thuật (research streams) chủ đạo:
Dòng nghiên cứu thứ nhất tập trung vào dao động dầm FGM 1D biến đổi theo chiều cao. Các công trình kinh điển của Aydogdu & Taskin (2007) sử dụng nghiệm giải tích Navier cho dầm P-FGM và E-FGM; Thai & Vo (2012) phát triển lý thuyết biến dạng trượt bậc cao loại bỏ nhu cầu sử dụng hệ số điều chỉnh trượt; Sayyad & Ghugal (2015) tổng hợp các biến thể lý thuyết dầm sandwich FGM. Hướng tiếp cận này sau đó được mở rộng sang dầm FGM biến thiên dọc trục (axial gradient) bởi Alshorbagy et al. (2011), Shahba et al. (2011), và Huang & Li (2010). Điểm giới hạn chung của dòng nghiên cứu này là chỉ khảo sát gradient 1D đơn hướng, không phản ánh được sự biến thiên trường cơ tính phức tạp trong các kết cấu thực tế chịu tải trọng cục bộ hoặc nhiệt trường phân bố hai chiều.
Dòng nghiên cứu thứ hai hướng tới vật liệu 2D-FGM và kết cấu tấm sandwich FGM. Wang et al. (2016) và Karamanlı (2017) khảo sát dao động dầm 2D-FGM biến thiên theo hàm số mũ; Lieu et al. (2018) và Li et al. (2019) sử dụng phương pháp đẳng hình học (IGA) phân tích uốn và dao động của tấm 2D-FGM. Đối với tấm sandwich FGM, các nghiên cứu của Zenkour (2005), Neves et al. (2012) với lý thuyết quasi-3D, và Thai et al. (2014) đã đặt nền móng giải tích quan trọng. Tuy nhiên, các lời giải giải tích truyền thống chỉ giới hạn ở biên tựa giản đơn và không thể giải quyết được hệ phương trình vi phân chuyển động có hệ số phụ thuộc biến không gian $x$ phức tạp của kết cấu sandwich 2D-FGM đa pha.
Dòng nghiên cứu thứ ba tập trung vào ứng xử động lực học dưới tải trọng di động. Các công bố đột phá của Şimşek & Kocatürk (2009), Şimşek (2010), Khalili et al. (2010), và Gan et al. (2015) đã chứng minh rằng vận tốc di động, hiệu ứng quán tính và tần số kích động làm thay đổi căn bản cơ chế đáp ứng động lực học của kết cấu composite. Dù vậy, bài toán tải trọng di động trên kết cấu sandwich 2D-FGM hai pha và ba pha hầu như chưa được khai phá.
Trong y văn tồn tại hai cuộc tranh luận lý thuyết chưa ngã ngũ:
- Tranh luận về mô hình đồng nhất hóa: Mô hình Voigt cổ điển dựa trên quy tắc hỗn hợp tuyến tính đơn giản nhưng bỏ qua tương tác vi cấu trúc hạt, trong khi mô hình Mori-Tanaka (1973) và Maxwell (1873) phản ánh chính xác nồng độ ứng suất cục bộ của pha cốt trong nền. Đa số nghiên cứu trước đây áp dụng Voigt do tính thuận tiện toán học mà thiếu sự đánh giá sai số định lượng trên kết cấu 2D gradient.
- Tranh luận về bậc của trường chuyển vị: Lý thuyết cổ điển Euler-Bernoulli và Mindlin (FSDT) dẫn đến hiện tượng nghẽn trượt nghiêm trọng khi tỷ số $L/h$ nhỏ hoặc cần hệ số hiệu chỉnh trượt khó xác định cho vật liệu dị hướng không đồng nhất.
So sánh trực tiếp với hai nghiên cứu quốc tế tiêu biểu:
- Nghiên cứu của Songsuwan et al. (2018) phân tích dầm sandwich 1D-FGM chịu lực điều hòa di động bằng phương pháp Ritz nhưng bị bó hẹp trong phạm vi biến thiên một chiều ($n_x = 0$) và chỉ khảo sát mô hình tải trọng lực thuần túy.
- Nghiên cứu của Karamanlı (2017) sử dụng lý thuyết biến dạng trượt bậc ba phân tích uốn tĩnh dầm sandwich 2D-FGM hai pha bằng đa thức xấp xỉ nhưng không giải quyết bài toán động lực học phi dừng và dao động tự do dưới tải trọng di động.
Luận án này định vị tại điểm giao thoa tiên phong: tích hợp đồng thời cơ tính biến thiên hai chiều ($x, z$), cấu trúc sandwich đa pha (2 pha và 3 pha), kết hợp tải trọng quán tính di động phức tạp thông qua việc xây dựng hệ phần tử hữu hạn chuyên dụng làm giàu hàm nội suy và thuật toán tích phân thời gian Newmark ẩn.
Đóng góp lý thuyết và khung phân tích
Đóng góp cho lý thuyết
Luận án mở rộng và làm sâu sắc thêm ba hệ thống lý thuyết nền tảng trong cơ học kết cấu:
-
Mở rộng lý thuyết biến dạng trượt bậc ba Bickford-Reddy: Ứng dụng thành công trường chuyển vị phân tách uốn - trượt cho kết cấu sandwich 2D-FGM hai pha có mô-men tĩnh và độ cứng biến thiên hai chiều $A_{ij}(x), B_{ij}(x), D_{ij}(x)$. Trường chuyển vị cải tiến triệt tiêu hoàn toàn ứng suất trượt trên các bề mặt tự do trên và dưới dầm ($\tau_{xz}(x, \pm h/2) = 0$), loại bỏ triệt để sự phụ thuộc vào hệ số hiệu chỉnh trượt nhân tạo.
-
Phát triển lý thuyết biến dạng trượt lượng giác cho dầm 2D-FGM ba pha: Thiết lập hệ phương trình vi phân chuyển động phi tuyến động lực học có xét đến ảnh hưởng đồng thời của lực quán tính dịch chuyển, lực Coriolis và lực ly tâm sinh ra bởi khối lượng di động $m$:
$$\mathbf{m}m = m \mathbf{N}^T \mathbf{N}, \quad \mathbf{c}m = 2 m v \mathbf{N}^T \mathbf{N}{,x}, \quad \mathbf{k}m = m v^2 \mathbf{N}^T \mathbf{N}{,xx} + m \ddot{w}(x_m, t) \mathbf{N}^T \mathbf{N}{,x}$$
-
Chuyển dịch mô hình vi cơ học đồng nhất hóa đa pha: Định lượng hóa phạm vi ứng dụng và sai số biên của mô hình quy tắc hỗn hợp Voigt so với mô hình trường tự nhất quán Mori-Tanaka và mô hình Maxwell trong việc xác định các mô-đun đàn hồi tương đương $E_f(x, z)$ và $G_f(x, z)$ của composite 2D-FGM hai pha và ba pha.
Khung phân tích độc đáo
Khung phân tích được xây dựng trên sự tích hợp chặt chẽ của ba trụ cột lý thuyết: Cơ học môi trường liên tục biến phân, Lý thuyết vi cơ học đồng nhất hóa, và Phương pháp phần tử hữu hạn nâng cao.
Mô hình phân bố tỷ phần thể tích lũy thừa 2D của dầm sandwich hai pha được mô tả chính xác qua biểu thức:
$$V_m(x, z) = \begin{cases} \left(\frac{z - z_1}{z_0 - z_1}\right)^{n_z} \left(1 - \frac{x}{2L}\right)^{n_x}, & z \in [z_0, z_1] \ 0, & z \in [z_1, z_2] \ \left(\frac{z - z_2}{z_3 - z_2}\right)^{n_z} \left(1 - \frac{x}{2L}\right)^{n_x}, & z \in [z_2, z_3] \end{cases}$$
với $V_c(x, z) = 1 - V_m(x, z)$.
Mô hình Mori-Tanaka được thiết lập dựa trên mô-đun khối $K_f$ và mô-đun trượt $G_f$ hiệu dụng:
$$E_f(x, z) = \frac{9 K_f G_f}{3 K_f + G_f}, \quad \nu_f(x, z) = \frac{3 K_f - 2 G_f}{6 K_f + 2 G_f}$$
trong đó:
$$\frac{K_f - K_m}{K_c - K_m} = \frac{V_c}{1 + V_m \frac{K_c - K_m}{K_m + \frac{4}{3}G_m}}, \quad \frac{G_f - G_m}{G_c - G_m} = \frac{V_c}{1 + V_m \frac{G_c - G_m}{G_m + G_m \frac{9K_m + 8G_m}{6K_m + 12G_m}}}$$
Điều kiện biên xác định (boundary conditions): Áp dụng cho các cấu hình liên kết thực tế gồm tựa đơn (Simple support - S), ngàm cứng (Clamped - C), và tự do (Free - F), đảm bảo tính đại diện cho kết cấu công trình nhịp lớn và vỏ phương tiện bay.
Phương pháp nghiên cứu tiên tiến
Thiết kế nghiên cứu
Nghiên cứu tuân thủ chặt chẽ lập trường nhận thức luận thực chứng (Positivism) và phương pháp luận diễn dịch toán - số (deductive computational mechanics). Thiết kế nghiên cứu kết hợp mô hình hóa giải tích biến phân với mô phỏng số phần tử hữu hạn nhiều mức (multi-level discrete modeling):
- Cấp độ dầm 2D-FGM 2 pha: Sử dụng phần tử dầm 2 nút với 10 bậc tự do (DOFs):
$$\mathbf{d}e = {u_1, u_2, w{b1}, w_{b1,x}, w_{b2}, w_{b2,x}, w_{s1}, w_{s1,x}, w_{s2}, w_{s2,x}}^T$$
kết hợp hàm nội suy tuyến tính cho chuyển vị dọc trục $u$ và đa thức Hermite bậc ba cho chuyển vị uốn $w_b$ và trượt $w_s$.
- Cấp độ dầm 2D-FGM 3 pha: Ứng dụng trường nội suy làm giàu (enriched trigonometric shape functions) kết hợp hàm dạng Hermite cải tiến, loại bỏ nhu cầu tăng mật độ lưới phần tử mà vẫn đảm bảo độ chính xác vượt trội.
- Cấp độ tấm 2D-FGM 3 pha: Phát triển phần tử tứ giác 9 nút Q9 (Quadrilateral 9-node) với 45 bậc tự do tổng thể, kết hợp kỹ thuật nội suy liên kết (tied/mixed interpolation) cho các trường biến dạng trượt ngang $\gamma_{xz}, \gamma_{yz}$, triệt tiêu hoàn toàn hiện tượng khóa trượt (shear locking) khi tỷ số độ mảnh $a/h$ tiến tới vô cùng.
Quy trình nghiên cứu rigorous
Quy trình tính toán số được xây dựng theo một chu trình nghiêm ngặt:
- Thiết lập toán học: Áp dụng nguyên lý Hamilton thu được hệ phương trình vi phân chuyển động liên tục dạng đạo hàm riêng cấp 4 có hệ số biến thiên theo $x$.
- Rời rạc hóa không gian: Sử dụng phép cầu phương Gauss (Gauss Quadrature) bậc cao để tính toán tích phân các ma trận độ cứng $\mathbf{k}e$, ma trận khối lượng $\mathbf{m}e$, và véc-tơ lực nút $\mathbf{f}e$ do tính không đồng nhất của vật liệu:
$$\mathbf{k}{aa} = \int_0^l \mathbf{N}{,x}^T A{11}(x) \mathbf{N}{,x} dx, \quad \mathbf{k}{bb} = \int_0^l \mathbf{H}{,xx}^T A{22}(x) \mathbf{H}_{,xx} dx$$
- Mô hình hóa cản: Tích hợp ma trận cản Rayleigh tổng thể $\mathbf{C} = \alpha \mathbf{M} + \beta \mathbf{K}$, với các hệ số cản $\alpha, \beta$ được xác định chính xác từ tỷ số cản thực nghiệm $\xi_i = \xi_j = 2.5%$ tương ứng với hai tần số dao động riêng đầu tiên $\omega_1, \omega_2$.
- Tích phân thời gian trực tiếp: Triển khai thuật toán Newmark gia tốc trung bình ẩn (Average Acceleration Method) với tính ổn định vô điều kiện (unconditionally stable):
$$\mathbf{K}{\text{eff}} \mathbf{D}{n+1} = \mathbf{F}{\text{eff}}^{n+1}$$
trong đó ma trận độ cứng hiệu dụng và véc-tơ tải trọng hiệu dụng tại thời điểm $(n+1)\Delta t$ có dạng:
$$\mathbf{K}{\text{eff}} = \frac{4}{\Delta t^2} \mathbf{M} + \frac{2}{\Delta t} \mathbf{C} + \mathbf{K}$$
$$\mathbf{F}{\text{eff}}^{n+1} = \mathbf{F}{n+1} + \mathbf{M}\left(\frac{4}{\Delta t^2}\mathbf{D}_n + \frac{4}{\Delta t}\mathbf{\dot{D}}_n + \mathbf{\ddot{D}}_n\right) + \mathbf{C}\left(\frac{2}{\Delta t}\mathbf{D}_n + \mathbf{\dot{D}}_n\right)$$
Data và phân tích
Toàn bộ hệ thống mã nguồn mô phỏng được lập trình độc lập bằng công cụ tính toán ma trận và đại số biểu thức Symbolic trên nền tảng MATLAB và Maple.
Kiểm tra độ hội tụ và độ ổn định lưới:
- Kết quả khảo sát sự hội tụ trên dầm sandwich 2D-FGM ba pha ($L/h = 20, r_m = 0.5, v = 50\text{ m/s}$) chứng minh rằng phần tử dầm làm giàu đạt trạng thái hội tụ tuyệt đối ở tham số tần số cơ bản $\mu_1 = 2.8452$ và hệ số động lực học $D_d = 1.2418$ chỉ với $NE = 20$ phần tử, vượt trội hơn hẳn so với phần tử Lagrange bậc một truyền thống đòi hỏi $NE > 60$.
- Đối với phần tử tấm Q9, sai số tương đối so với nghiệm giải tích 3D Elasticity của dầm/tấm 1D-FGM hai pha chỉ dao động trong khoảng $0.08% - 0.35%$.
Phát hiện đột phá và implications
Những phát hiện then chốt
-
Quy luật chi phối của gradient hai chiều ($n_x, n_z$): Khi đồng thời tăng các chỉ số gradient từ $n_x = n_z = 0.2$ lên $n_x = n_z = 5$, tham số tần số dao động cơ bản $\mu_1$ của dầm sandwich 2D-FGM hai pha (cấu hình 2-1-2, $L/h = 10$) giảm mạnh từ $2.7841$ xuống $1.9125$ (mức suy giảm $31.3%$). Nguyên nhân vật lý cốt lõi là sự sụt giảm tỷ phần thể tích của pha gốm ($E_c = 380\text{ GPa}, \rho_c = 3800\text{ kg/m}^3$) và sự gia tăng chiếm lĩnh của pha kim loại ($E_m = 70\text{ GPa}, \rho_m = 2707\text{ kg/m}^3$), làm mềm hóa kết cấu một cách rõ rệt.
-
Sự sai lệch mang tính hệ thống giữa các mô hình vi cơ học: Kết quả so sánh định lượng khẳng định mô hình Voigt luôn ước tính tần số dao động riêng cao hơn mô hình Mori-Tanaka từ $3.2%$ đến $8.7%$, và cao hơn mô hình Maxwell khoảng $4.1% - 9.3%$. Tương ứng, hệ số động lực học $D_d = w_{\text{dynamic}}^{\max} / w_{\text{static}}^{\max}$ tính theo mô hình Mori-Tanaka luôn lớn hơn so với mô hình Voigt từ $5.4% - 11.2%$ tại cùng dải vận tốc $v = 50 - 80\text{ m/s}$. Điều này chỉ ra rằng việc sử dụng mô hình Voigt truyền thống có thể dẫn đến các thiết kế thiếu an toàn trong thực tế công trình.
| Cấu hình dầm / Tham số |
Mô hình Voigt ($\mu_1$) |
Mô hình Mori-Tanaka / Maxwell ($\mu_1$) |
Mức sai lệch tương đối (%) |
Hệ số động lực học $D_d$ (Voigt) |
Hệ số động lực học $D_d$ (MT/Maxwell) |
| Dầm 2D-FGM 2P ($L/h = 5, n_x=n_z=0.5$) |
2.5412 |
2.4418 |
+4.07% |
1.1824 |
1.2510 |
| Dầm 2D-FGM 2P ($L/h = 20, n_x=n_z=1.0$) |
2.1145 |
2.0120 |
+5.09% |
1.3412 |
1.4580 |
| Dầm 2D-FGM 3P ($L/h = 10, n_x=n_z=0.5$) |
2.6840 |
2.5210 |
+6.46% |
1.2150 |
1.3320 |
| Tấm 2D-FGM 3P SSSS ($a/h = 10, n_x=n_z=1$) |
1.8420 |
1.7150 |
+7.40% |
-- |
-- |
| Tấm 2D-FGM 3P CCCC ($a/h = 10, n_x=n_z=1$) |
3.4210 |
3.1980 |
+6.97% |
-- |
-- |
-
Hiện tượng dịch chuyển đỉnh đáp ứng và biến dạng ứng suất phi đối xứng: Dưới tác dụng của tải trọng di động, sự biến thiên tính chất theo phương dọc trục ($n_x > 0$) làm phá vỡ tính đối xứng hình học của đường cong độ võng thời gian $w(L/2, t)$. Vị trí độ võng cực đại không còn xuất hiện khi tải trọng đi qua chính giữa nhịp ($t = \Delta T / 2$) mà bị trễ về phía nửa sau dầm tùy thuộc vào vận tốc $v$. Đồng thời, phân bố ứng suất pháp $\sigma_{xx}$ và ứng suất trượt $\tau_{xz}$ theo chiều cao tiết diện biểu hiện tính chất phi tuyến bậc cao phức tạp, với vị trí trục trung hòa động liên tục dịch chuyển theo tọa độ không gian $x$.
-
Hiệu ứng tương tác khối lượng di động (Moving Mass Effect): Khác biệt hoàn toàn với mô hình lực di động thuần túy (Moving Force), sự xuất hiện của tỷ số khối lượng $r_m = m / (\rho A L) > 0$ sinh ra lực quán tính phụ và lực Coriolis làm giảm tần số dao động tức thời của hệ khi khối lượng tiến vào dầm. Khi vận tốc vượt qua ngưỡng tới hạn $v_{\text{critical}} \approx 65\text{ m/s}$ (với $L/h = 20, r_m = 0.5$), biên độ dao động dư sau khi tải trọng rời khỏi dầm tăng vọt, gây ra ứng suất trượt $\tau_{xz}$ cục bộ rất lớn tại gối tựa.
-
Hiệu năng chống nghẽn trượt của phần tử tấm Q9: Phần tử Q9 với kỹ thuật nội suy liên kết duy trì độ chính xác tuyệt đối đối với cả tấm rất dày ($a/h = 5$) và tấm rất mỏng ($a/h = 100$), loại bỏ hoàn toàn hiện tượng khóa cứng ma trận độ cứng thường gặp trong các phần tử tấm Mindlin đẳng tham số bậc thấp tiêu chuẩn.
Implications đa chiều
- Về mặt lý thuyết: Cung cấp hệ thống bằng chứng toán học và số học xác nhận vai trò sống còn của việc mô hình hóa gradient hai chiều và lực quán tính Coriolis trong tính toán động lực học kết cấu composite thế hệ mới.
- Về phương pháp luận: Chuẩn hóa quy trình phần tử hữu hạn làm giàu và phần tử tấm nội suy liên kết, tạo tiền đề mở rộng cho các bài toán phân tích phi tuyến hình học, kết cấu vỏ thoải và vật liệu xốp có lỗ rỗng (porous FGM).
- Về ứng dụng kỹ thuật và thực tiễn: Cung cấp cơ sở khoa học để tối ưu hóa thiết kế phân bố tỷ phần thể tích vật liệu trong các dầm cầu đường sắt tốc độ cao, dầm dẫn hướng tàu đệm từ, tấm sàn chịu tải trọng xe cơ giới hạng nặng, và vỏ bảo vệ động cơ phản lực hàng không, giúp giảm trọng lượng bản thân từ $15% - 25%$ trong khi vẫn kiểm soát hoàn hảo biên độ rung động động học.
Limitations và Future Research
Mặc dù đạt được nhiều kết quả đột phá, công trình vẫn tồn tại một số giới hạn học thuật cần được làm rõ:
- Giới hạn miền tuyến tính: Nghiên cứu giả định ứng xử vật liệu đàn hồi tuyến tính và chuyển vị biến dạng nhỏ, chưa xét đến hiệu ứng phi tuyến hình học Von Kármán khi dầm/tấm chịu tải trọng động lực học cực lớn.
- Điều kiện môi trường nhiệt độ: Chưa tích hợp trường nhiệt độ không dừng và gradient nhiệt độ phức tạp gây biến đổi cơ lý tính phụ thuộc nhiệt độ ($E(T), \rho(T)$).
- Hình thái khuyết tật vật liệu: Mô hình giả định vật liệu liên tục lý tưởng, chưa tính đến ảnh hưởng của lỗ rỗng tế vi ngẫu nhiên (micro-porosities) sinh ra trong quá trình chế tạo thiêu kết gốm - kim loại.
Định hướng nghiên cứu tiếp theo:
- Mở rộng thuật toán PTHH sang phân tích phi tuyến hình học và động lực học phi tuyến của dầm, vỏ sandwich 2D-FGM dưới tác động của tải trọng va đập động năng cao.
- Thiết lập mô hình bài toán nhiệt - cơ - điện từ ghép nối (Thermo-Electro-Mechanical coupled analysis) cho dầm/tấm sandwich 2D-FGM áp điện thông minh.
- Ứng dụng các thuật toán tối ưu hóa đa mục tiêu (như giải thuật di truyền GA, tối ưu hóa bầy đàn PSO) để tìm kiếm quy luật phân bố gradient tối ưu nhằm triệt tiêu dao động cộng hưởng.
- Thực hiện các thử nghiệm thực nghiệm trên các mẫu in 3D FGM đa vật liệu để hiệu chuẩn chính xác hơn các hệ số cản và mô hình đồng nhất hóa.
Tác động và ảnh hưởng
- Tác động học thuật: Luận án thiết lập một hệ thống dữ liệu chuẩn (benchmark numerical data) phong phú và có độ tin cậy cao cho bài toán dao động dầm và tấm sandwich 2D-FGM. Các ấn phẩm khoa học xuất bản từ luận án trên các tạp chí quốc tế chuyên ngành uy tín thuộc danh mục ISI/Scopus (như Composite Structures, Thin-Walled Structures, Journal of Sandwich Structures & Materials) đóng góp trực tiếp vào chỉ số trích dẫn và nâng cao vị thế khoa học của cơ học tính toán Việt Nam.
- Tác động công nghiệp: Cung cấp công cụ phần mềm tính toán chuyên dụng cho các kỹ sư R&D trong ngành công nghiệp quốc phòng, hàng không vũ trụ và giao thông vận tải, hỗ trợ trực tiếp việc thiết kế các kết cấu nhẹ có khả năng chống rung động, chống mỏi và cách nhiệt vượt trội.
- Lợi ích kinh tế - xã hội: Việc tối ưu hóa vật liệu 2D-FGM giúp giảm thiểu việc sử dụng các kim loại quý hiếm, giảm tải trọng bản thân của kết cấu công trình, kéo dài chu kỳ bảo dưỡng và nâng cao mức độ an toàn vận hành của hệ thống giao thông tốc độ cao.
Đối tượng hưởng lợi
- Nghiên cứu sinh và Giảng viên đại học: Tiếp cận được khung lý thuyết biến dạng trượt bậc cao hoàn chỉnh, các công thức giải tích biến phân tường minh và phương pháp lập trình phần tử hữu hạn nâng cao trên MATLAB/Maple để phục vụ công tác giảng dạy và nghiên cứu chuyên sâu.
- Kỹ sư tính toán kết cấu (CAE Engineers): Sở hữu công cụ số và các biểu đồ khảo sát tham số chi tiết để lựa chọn tỷ số hình học, cấu hình lớp sandwich và gradient vật liệu phù hợp cho các cấu kiện chịu tải trọng di động.
- Các nhà hoạch định tiêu chuẩn thiết kế: Có thêm cơ sở dữ liệu thực nghiệm số vững chắc để xây dựng các quy chuẩn kỹ thuật mới cho việc ứng dụng vật liệu composite biến tính chức năng trong xây dựng công trình giao thông và phương tiện vận tải hiện đại.
Câu hỏi chuyên sâu
-
Đóng góp lý thuyết độc đáo nhất của luận án là gì?
- Đóng góp độc đáo nhất là việc mở rộng thành công lý thuyết biến dạng trượt bậc ba Bickford-Reddy và lý thuyết lượng giác cho kết cấu sandwich 2D-FGM hai và ba pha, đồng thời định lượng hóa sự phân kỳ của các mô hình đồng nhất hóa vi cơ học (Voigt, Mori-Tanaka, Maxwell) trong điều kiện gradient hai chiều $(x, z)$.
-
Cải tiến phương pháp luận có gì vượt trội so với các nghiên cứu trước đây?
- So với phương pháp Ritz của Songsuwan et al. (2018) hay giải tích Navier cổ điển vốn chỉ áp dụng cho dầm 1D-FGM biên khớp, luận án đã phát triển thành công phần tử dầm làm giàu hàm lượng giác 10 bậc tự do và phần tử tấm Q9 nội suy liên kết. Kỹ thuật này xử lý hoàn hảo các hệ số biến thiên không gian phức tạp và loại bỏ triệt để hiện tượng nghẽn trượt trên mọi tỷ số độ mảnh $a/h \in [5, 100]$.
-
Phát hiện nào mang tính bất ngờ nhất từ kết quả số?
- Phát hiện bất ngờ nhất là việc mô hình Voigt đánh giá quá mức độ cứng của dầm sandwich 2D-FGM, dẫn đến việc dự báo tần số riêng cao hơn tới $8.7%$ và độ võng động học nhỏ hơn tới $11.2%$ so với mô hình Mori-Tanaka/Maxwell. Điều này chỉ ra rằng các tiêu chuẩn thiết kế dựa trên quy tắc hỗn hợp Voigt có thể dẫn đến nguy cơ mất an toàn nghiêm trọng dưới tải trọng động.
-
Luận án có cung cấp đầy đủ giao thức để tái lập kết quả (Replication Protocol) không?
- Hoàn toàn đầy đủ. Toàn bộ công thức ma trận phần tử $(\mathbf{k}_e, \mathbf{m}_e, \mathbf{f}_e)$, thuật toán tích phân thời gian Newmark ẩn, các bảng thông số vật liệu chi tiết (gốm Alumina, Zirconia; kim loại Nhôm, Thép, Ti-6Al-4V) và mã nguồn tính toán trên MATLAB/Maple đều được trình bày minh bạch trong văn bản và hệ thống phụ lục.
-
Chương trình nghiên cứu 10 năm tiếp theo được định hình như thế nào?
- Lộ trình 10 năm tập trung vào: (1) Mở rộng bài toán sang miền phi tuyến lớn và va đập tốc độ cao; (2) Tích hợp trường nhiệt - điện - từ thông minh; (3) Tối ưu hóa cấu trúc gradient bằng trí tuệ nhân tạo (AI-driven inverse design); và (4) Chế tạo mẫu thực nghiệm bằng công nghệ in 3D kim loại - gốm FGM để thương mại hóa ứng dụng.
Kết luận
Công trình nghiên cứu của Nghiên cứu sinh Phạm Vũ Nam đã hoàn thành xuất sắc các mục tiêu khoa học đề ra với 5 đóng góp cốt lõi:
- Thiết lập thành công mô hình toán học giải tích chặt chẽ dựa trên lý thuyết biến dạng trượt bậc ba cải tiến và lý thuyết lượng giác cho dầm và tấm sandwich 2D-FGM hai và ba pha.
- Phát triển hệ phần tử hữu hạn dầm hai nút làm giàu 10 bậc tự do và phần tử tấm tứ giác Q9 nội suy liên kết, khắc phục triệt để hiện tượng nghẽn trượt và tối ưu hóa tốc độ hội tụ.
- Làm sáng tỏ cơ chế tương tác động lực học phức tạp giữa tải trọng di động (lực điều hòa và khối lượng di động có xét lực Coriolis, ly tâm) với kết cấu sandwich 2D gradient.
- Đánh giá định lượng toàn diện sai số giữa các mô hình đồng nhất hóa vật liệu Voigt, Mori-Tanaka và Maxwell, đưa ra khuyến nghị kỹ thuật then chốt cho tính toán an toàn công trình.
- Cung cấp bộ dữ liệu số chuẩn xác thực cao, mở ra ba hướng nghiên cứu mới về cơ học kết cấu composite tiên tiến, đóng góp quan trọng vào sự phát triển của cơ học tính toán trong nước và quốc tế.
