Tổng quan về luận án

Nghiên cứu ứng xử cơ học của các kết cấu vỏ composite nhiều lớp (laminated composite shells) và composite có cơ tính biến thiên (Functionally Graded Materials - FGM) đóng vai trò then chốt trong kỹ thuật hàng không vũ trụ, công trình biển, tàu ngầm và công nghiệp quốc phòng. Vỏ là kết cấu chịu lực tối ưu nhờ sự kết hợp giữa hiệu ứng màng (membrane) và hiệu ứng uốn (bending). Tuy nhiên, việc mô hình hóa số chính xác vỏ mỏng và vỏ có độ dày trung bình chịu tải trọng tĩnh, dao động tự do, phi tuyến hình học và tải trọng sóng xung kích nổ ngầm dưới nước (underwater shock wave) đặt ra nhiều thách thức tính toán nghiêm ngặt, đặc biệt là hiện tượng "khóa cắt" (shear locking) và "khóa màng" (membrane locking) trong phương pháp phần tử hữu hạn cổ điển (FEM).

Khoảng trống nghiên cứu (Research Gap) xuất phát từ hạn chế của các phương pháp hiện hữu: Các phương pháp giải tích (Reddy, 2004; Đào Huy Bích et al., 2013) chỉ giải quyết được các hình học đối xứng và điều kiện biên lý tưởng; phần tử tam giác 3 nút tiêu chuẩn trong FEM truyền thống bị hiện tượng quá cứng (overly stiff) và khóa cắt nghiêm trọng khi tỷ lệ chiều dài trên chiều dày $a/h > 100$; trong khi các kỹ thuật nội suy hỗn hợp thành phần ten-xơ (MITC3 của Lee & Bathe, 2004) hoặc phương pháp làm trơn trên miền (CS-FEM của Liu et al., 2007) chưa giải quyết triệt để sự cân bằng giữa độ mềm lưới tam giác và triệt tiêu hoàn toàn khóa cắt cho vỏ composite có độ cong phức tạp.

Các câu hỏi nghiên cứu và giả thuyết khoa học của luận án bao gồm:

  1. RQ1: Làm thế nào để xây dựng một phần tử vỏ tam giác 3 nút phẳng mới vừa triệt tiêu hoàn toàn hiện tượng khóa cắt vừa khắc phục độ cứng giả tạo của phần tử tam giác chuẩn? Giả thuyết H1: Việc kết hợp kỹ thuật làm trơn biến dạng trên cạnh (Edge-based Smoothed FEM - ES-FEM) với kỹ thuật nội suy hỗn hợp các thành phần ten-xơ (MITC3) sẽ tạo ra phần tử ES-MITC3 có ma trận độ cứng tối ưu hóa năng lượng biến dạng, đảm bảo độ hội tụ đơn điệu và chính xác cao.
  2. RQ2: Ứng xử phi tuyến hình học chuyển vị lớn và các điểm mất ổn định tới hạn (snap-through, snap-back) của vỏ FGM và vỏ composite lớp được mô phỏng chính xác ra sao? Giả thuyết H2: Tích hợp thuật toán chiều dài cung (Arc-length) với phương pháp lặp Newton-Raphson trên nền tảng phần tử ES-MITC3 sẽ vượt qua các điểm kỳ dị ma trận tiếp tuyến để theo dõi chính xác toàn bộ đường cong cân bằng lực - chuyển vị.
  3. RQ3: Đáp ứng động lực học phi tuyến của vỏ composite 2 độ cong dưới tác dụng của sóng xung kích nổ trong môi trường nước biến thiên như thế nào theo thời gian và tham số kết cấu? Giả thuyết H3: Sự kết hợp giữa thuật toán tích phân thời gian ẩn Newmark, phương pháp lặp Newton-Raphson và mô hình tải trọng sóng xung kích Cole cho phép dự báo ứng xử quá độ chính xác dưới tác động của khối lượng nổ TNT ($M = 10\text{ kg}$) tại khoảng cách nổ ($S = 3\text{ m}$).

Nghiên cứu ứng dụng khung lý thuyết biến dạng cắt bậc nhất (First-order Shear Deformation Theory - FSDT) của Reissner-Mindlin, lý thuyết biến dạng lớn phi tuyến hình học von Kármán, nguyên lý công khả dĩ (Lagrange principle) và cơ học môi trường liên tục. Phạm vi khảo sát tập trung vào vỏ cầu ($R_x = R_y = R$), vỏ trụ ($R_x = R, R_y = \infty$), vỏ hyperbol/yên ngựa ($R_x = -R_y$), với các góc xếp lớp $[0^\circ/90^\circ/0^\circ]$, $[0^\circ/90^\circ/90^\circ/0^\circ]$, $[45^\circ/-45^\circ/0^\circ]$, cùng dải chiều dày $t = 6.35\text{ mm}$ đến $t = 25.4\text{ mm}$ và hệ số mũ thể tích FGM $n \in [0, 10]$.

Literature Review và Positioning

Các nghiên cứu giải tích và số về vỏ composite đã trải qua quá trình phát triển sâu rộng với ba nhánh học thuật chính:

Nhánh thứ nhất tập trung vào lý thuyết giải tích và bán giải tích cho tấm và vỏ composite. Reddy (2004) phát triển lý thuyết biến dạng cắt bậc cao (HSDT) loại bỏ hệ số hiệu chỉnh cắt $\kappa = 5/6$. Tornabene và Viola (2009, 2013) phát triển phương pháp cầu phương vi phân tổng quát (GDQ) giải bài toán tĩnh và dao động vỏ composite 2 độ cong. Tại Việt Nam, Đào Huy Bích và cộng sự (2013, 2015), Nguyễn Đình Đức (2014, 2016) công bố các công trình giải tích xuất sắc về vỏ FGM và vỏ có gân gia cường trên nền đàn hồi Pasternak bằng phương pháp Bubnov-Galerkin và Runge-Kutta. Dù vậy, phương pháp giải tích tồn tại rào cản lớn khi xử lý hình học tùy ý và điều kiện biên phức hợp trong kỹ thuật thực tế.

Nhánh thứ hai nghiên cứu kỹ thuật xử lý hiện tượng khóa cắt trong FEM. Khi tỷ lệ kích thước cạnh trên chiều dày $a/h > 100$, phần tử Mindlin 3 nút và 4 nút cổ điển xuất hiện hiện tượng khóa cắt nghiêm trọng (shear locking). Zienkiewicz et al. (1971) đề xuất tích phân giảm (reduced integration) nhưng làm xuất hiện các mode dao động năng lượng zero (spurious zero-energy modes). Bletzinger et al. (2000) đề xuất phương pháp rời rạc lệch trượt (Discrete Shear Gap - DSG3). Đặc biệt, Bathe và Dvorkin (1985), Lee và Bathe (2004), Jeon et al. (2015) đã thiết lập phương pháp nội suy hỗn hợp các thành phần ten-xơ (MITC3, MITC4, MITC3+) dựa trên việc nội suy độc lập trường biến dạng trượt ngang tại các điểm lấy mẫu hình học trên cạnh phần tử.

Nhánh thứ ba là sự ra đời của phương pháp phần tử hữu hạn trơn (Smoothed Finite Element Method - SFEM) do G. R. Liu và cộng sự khởi xướng (2007, 2009). Bằng cách chuyển đổi tích phân đạo hàm hàm dạng trên diện tích phần tử sang tích phân đường bao trên các miền làm trơn, SFEM mang lại ma trận độ cứng "mềm" hơn, cải thiện đáng kể độ hội tụ nghiệm chuyển vị và trường ứng suất. Các biến thể bao gồm: làm trơn trên miền tế bào (Cell-based CS-FEM), làm trơn trên nút (Node-based NS-FEM) và làm trơn trên cạnh (Edge-based ES-FEM). Tại Việt Nam, Nguyễn Xuân Hùng et al. (2010), Phùng Văn Phúc et al. (2013) đã kết hợp ES-FEM với kỹ thuật DSG3 (tạo ra ES-DSG3) cho tấm phẳng.

Luận án này định vị chính xác tại điểm giao thoa chưa từng được khai phá: Tích hợp kỹ thuật làm trơn trên cạnh (ES-FEM) với phần tử tam giác vỏ phẳng nội suy ten-xơ hỗn hợp (MITC3) để hình thành phần tử vỏ phẳng ES-MITC3. Khi so sánh với công trình quốc tế của Lee & Bathe (2004) trên phần tử MITC3 thuần túy và Cui et al. (2011) trên phần tử CS-DSG3, phương pháp ES-MITC3 của luận án chứng minh tính vượt trội về tốc độ hội tụ năng lượng biến dạng, không phụ thuộc vào hiện tượng méo lưới (mesh distortion), đồng thời triệt tiêu hoàn toàn khóa cắt và khóa màng trên các bài toán vỏ thoải đa độ cong.

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

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

Luận án tạo ra bước tiến quan trọng trong cơ học tính toán kết cấu vỏ:

  1. Mở rộng lý thuyết phần tử hữu hạn trơn sang vỏ composite đa hướng và FGM: Luận án chứng minh rằng việc áp dụng phép toán làm trơn gradient biến dạng màng, uốn và cắt trên các miền trơn gắn với cạnh phần tử ($\Omega_k$) làm giảm độ cứng giả tạo của mô hình biến phân chuẩn, đưa nghiệm chuyển vị số tiệm cận sát hơn với cận trên của nghiệm giải tích chính xác.
  2. Xây dựng mô hình tương tác phi tuyến hình học theo trường biến dạng von Kármán: Mô hình hóa tường minh sự ghép nối giữa ma trận độ cứng màng $\mathbf{A}$, ma trận ghép nối uốn - nén $\mathbf{B}$ và ma trận độ cứng uốn $\mathbf{D}$ trong quan hệ nội lực: $$\begin{Bmatrix} \mathbf{N} \ \mathbf{M} \end{Bmatrix} = \begin{bmatrix} \mathbf{A} & \mathbf{B} \ \mathbf{B} & \mathbf{D} \end{bmatrix} \begin{Bmatrix} \hat{\boldsymbol{\varepsilon}}_m \ \hat{\boldsymbol{\kappa}} \end{Bmatrix}$$
  3. Chuyển dịch hệ hình (Paradigm shift) trong giải thuật mất ổn định phi tuyến: Phối hợp toán tử gradient làm trơn với phương pháp giải phi tuyến bán kính cung lặp (Arc-length method), luận án giải quyết trọn vẹn hiện tượng nhảy bước mất ổn định qua điểm giới hạn (snap-through) và thụt lùi (snap-back) mà các phương pháp kiểm soát tải trọng (load-control) truyền thống không thể hội tụ.

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

Khung phân tích được xây dựng dựa trên sự hợp nhất chặt chẽ của 4 nền tảng lý thuyết:

  • Lý thuyết biến dạng cắt bậc nhất Reissner-Mindlin (FSDT): Trường chuyển vị 6 bậc tự do tại mỗi nút $\mathbf{u} = [u_0, v_0, w_0, \beta_x, \beta_y, \beta_z]^T$, cho phép mô hình hóa chính xác hiệu ứng trượt ngang với hệ số hiệu chỉnh cắt $\kappa = 5/6$.
  • Kỹ thuật nội suy ten-xơ hỗn hợp MITC3: Trường biến dạng cắt tự nhiên $\gamma_{rz}, \gamma_{sz}$ được tính toán bằng phép chiếu nội suy từ các giá trị trượt đo dọc theo 3 cạnh của tam giác, đảm bảo biến dạng cắt không bị suy biến khi chiều dày $h \to 0$.
  • Toán tử làm trơn biến dạng trên cạnh (Edge-based Smoothing Operation): Miền tính toán $\Omega$ được phân chia thành các miền trơn $\Omega_k$ liên kết với từng cạnh $k$. Trường biến dạng màng trơn $\tilde{\boldsymbol{\varepsilon}}m^k$, biến dạng uốn trơn $\tilde{\boldsymbol{\kappa}}^k$ và biến dạng cắt trơn $\tilde{\boldsymbol{\gamma}}^k$ được xác định qua tích phân đường bao: $$\tilde{\boldsymbol{\varepsilon}}^k = \frac{1}{A_k} \int{\Omega_k} \boldsymbol{\varepsilon}(\mathbf{x}) d\Omega = \frac{1}{A_k} \oint_{\Gamma_k} \mathbf{n}_k(\mathbf{x}) \mathbf{u}(\mathbf{x}) d\Gamma$$
  • Lý thuyết tải trọng sóng xung kích ngầm nổ (Underwater Explosion - UNDEX): Tích hợp công thức thực nghiệm của Cole và áp lực truyền sóng xung kích phụ thuộc thời gian $P(t) = P_m e^{-t/\theta}$, phản ánh chính xác tương tác va đập chất lỏng - kết cấu vỏ.

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

Thiết kế nghiên cứu

Nghiên cứu được thiết kế theo trường phái thực chứng tính toán (Computational Positivism) và cơ học môi trường liên tục suy biến (Degenerated Continuum Mechanics). Mô hình sử dụng cách tiếp cận số đa tầng:

[Miền vỏ composite 3D] 
[Rời rạc hóa lưới tam giác phẳng 3 nút (MITC3)]
[Tạo lập hệ thống miền trơn gắn với cạnh phần tử (ES-FEM)]
[Thiết lập ma trận độ cứng trơn K̃, ma trận khối lượng M̃]
[Phân tích Tuyến tính/   [Phân tích Phi tuyến     [Động lực học Phi tuyến
 Dao động riêng]          Tĩnh (Arc-length + NR)]   (Newmark + NR + Sóng nổ)]

Quy trình nghiên cứu rigorous

  1. Rời rạc hóa và chuyển đổi hệ tọa độ: Miền kết cấu vỏ 3 chiều được rời rạc hóa thành lưới các phần tử tam giác phẳng 3 nút. Thiết lập hệ tọa độ cục bộ $\hat{O}\hat{x}\hat{y}\hat{z}$ cho từng phần tử và hệ tọa độ ảo cho từng miền trơn cạnh, liên kết với hệ tọa độ tổng thể $Oxyz$ qua ma trận chuyển đổi hình học Jacobi và ma trận cosin chỉ phương $\boldsymbol{\Lambda}_{0j}^i$.
  2. Tính toán ma trận độ cứng trơn: Ma trận biến dạng trơn $\tilde{\mathbf{B}}^k$ được thiết lập trực tiếp từ tọa độ các nút biên của miền trơn mà không cần đạo hàm hàm dạng cục bộ. Ma trận độ cứng của miền trơn được tính bằng tích phân trực tiếp không cần điểm Gauss bậc cao: $$\tilde{\mathbf{K}}k = \int{\Omega_k} \tilde{\mathbf{B}}_k^T \mathbf{D} \tilde{\mathbf{B}}_k d\Omega = A_k \tilde{\mathbf{B}}_k^T \mathbf{D} \tilde{\mathbf{B}}_k$$
  3. Thuật toán phi tuyến tĩnh và động:
    • Phi tuyến tĩnh: Kết hợp thuật toán lặp Newton-Raphson với kỹ thuật Arc-length kiểm soát độ dài cung $s = \sqrt{\Delta \mathbf{d}^T \Delta \mathbf{d} + \Delta \lambda^2 \psi^2 \mathbf{F}^T \mathbf{F}}$.
    • Động lực học phi tuyến: Sử dụng phương pháp tích phân từng bước Newmark ($\gamma = 0.5, \beta = 0.25$ - sơ đồ gia tốc trung bình không tiêu tán nhân tạo) kết hợp hiệu chỉnh lặp Newton-Raphson tại mỗi bước thời gian $\Delta t = 10^{-5}\text{ s}$ đến $10^{-4}\text{ s}$.

Data và phân tích

Toàn bộ thuật toán được tác giả tự lập trình hoàn toàn trong môi trường ma trận tính toán khoa học MATLAB. Các kiểm chứng độ tin cậy (Verification & Validation) được thực hiện đa chiều:

  • Dữ liệu phân tích tĩnh: Khảo sát chuyển vị không thứ nguyên $\bar{w} = 1000 w(a/2, b/2, 0) h^3 E_2 / (P a^4)$ và ứng suất $\bar{\sigma}_i = \sigma_i h^2 / (P a^2)$ trên vỏ 2 độ cong chịu tải phân bố đều, phân bố dạng sin và tải tập trung.
  • Dữ liệu dao động riêng: Tần số dao động không thứ nguyên $\bar{\omega} = \omega (a^2/h)\sqrt{\rho/E_2}$ cho 6 mode dao động đầu tiên trên các lưới $10 \times 10$, $20 \times 20$, $30 \times 30$.
  • Dữ liệu phi tuyến nổ ngầm: Tải trọng sóng nổ TNT $M = 10\text{ kg}$ đặt tại khoảng cách $S = 3\text{ m}$ trong môi trường nước có khối lượng riêng $\rho_{water} = 1000\text{ kg/m}^3$, vận tốc âm thanh $c_0 = 1460\text{ m/s}$, áp suất đỉnh sóng $P_m = 53.1 \times 10^6 (M^{1/3}/S)^{1.13}\text{ Pa}$, hằng số thời gian suy giảm $\theta = 0.084 M^{1/3} (M^{1/3}/S)^{-0.23}\text{ ms}$.

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

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

  1. Khắc phục triệt để hiện tượng khóa cắt với tốc độ hội tụ vượt bậc: Kết quả tính toán chuyển vị tâm vỏ cầu mỏng ($a/h = 1000$) bằng phần tử ES-MITC3 cho thấy sai số nhỏ hơn $0.12%$ so với nghiệm giải tích Navier, trong khi phần tử tam giác tiêu chuẩn FEM sai số lên đến $68.4%$ do nghẽn cắt. So với phần tử MITC3 chuẩn của Lee & Bathe (2004), ES-MITC3 đạt độ chính xác tương đương trên lưới thưa hơn $40%$ số bậc tự do.
  2. Khả năng nắm bắt hoàn hảo các trạng thái giới hạn phi tuyến phức tạp: Trên mô hình vỏ trụ FGM ($Al/ZrO_2$) chịu tải trọng tập trung tại tâm với hệ số mũ $n = 2$, chiều dày thay đổi từ $t = 6.35\text{ mm}$, $t = 12.7\text{ mm}$ đến $t = 25.4\text{ mm}$, phương pháp Arc-length ES-MITC3 đã mô tả liên tục và mượt mà hiện tượng nhảy lật dạng snap-through và snap-back mà các công bố trước đây sử dụng Newton-Raphson tiêu chuẩn bị phân kỳ tại điểm tới hạn.
  3. Quy luật đáp ứng động lực học phi tuyến dưới sóng nổ ngầm: Khi tăng khoảng cách nổ từ $S = 2\text{ m}$ lên $S = 5\text{ m}$ với lượng nổ $M = 10\text{ kg}$ TNT, biên độ chuyển vị đỉnh của tâm vỏ composite hai độ cong $[0^\circ/90^\circ/0^\circ]$ giảm phi tuyến đến $76.8%$. Đặc biệt, kết cấu vỏ hai độ cong có độ cứng động lực học cao hơn vỏ trụ và vỏ phẳng tương đương lần lượt là $34.2%$ và $58.6%$.
  4. Ảnh hưởng quyết định của cấu trúc lớp và hệ số biến thiên vật liệu: Vỏ FGM có cơ tính chuyển tiếp theo hàm mũ $V_c(z) = (z/h + 1/2)^n$ chuyển dịch từ hoàn toàn gốm ($n=0$) sang hoàn toàn kim loại ($n \to \infty$) làm giảm tần số dao động tự do bậc 1 tới $42.5%$ do sự suy giảm mô đun đàn hồi tương đương. Đối với vỏ composite lớp, góc hướng sợi $[45^\circ/-45^\circ/0^\circ]$ thể hiện khả năng kháng xoắn và phân tán áp lực sóng nổ ngầm vượt trội so với sơ đồ $[0^\circ/90^\circ/0^\circ]$.

Implications đa chiều

  • Về mặt lý thuyết: Chứng minh tính đúng đắn của việc áp dụng toán tử làm trơn trường biến dạng trên biên phần tử tam giác phẳng để mô phỏng hình học vỏ cong 3D, mở rộng giới hạn áp dụng của lý thuyết vỏ Reissner-Mindlin trong miền phi tuyến hình học lớn.
  • Về mặt phương pháp tính: Cung cấp giải thuật ES-MITC3 có tính tổng quát cao, lập trình thuận tiện, không cần khử bậc tự do xoay quanh trục pháp tuyến (drilling degree of freedom) phức tạp, có thể tích hợp trực tiếp vào các phần mềm thương mại CAE.
  • Về mặt thực tiễn kỹ thuật & Quốc phòng: Đưa ra cơ sở dữ liệu số chuẩn xác phục vụ thiết kế vỏ tàu ngầm, thân tàu chiến composite, kết cấu bảo vệ bồn chứa dầu khí ngoài khơi chống lại các mối đe dọa từ vũ khí nổ ngầm và va đập thủy động lực.
  • Về mặt chính sách & Quy chuẩn: Cung cấp công cụ tính toán mô phỏng thay thế các thử nghiệm nổ ngầm thực địa tốn kém và nguy hiểm, phục vụ việc xây dựng tiêu chuẩn thiết kế kết cấu vỏ tàu quân sự thế hệ mới.

Limitations và Future Research

Nghiên cứu có các giới hạn khoa học cần được ghi nhận minh bạch:

  1. Luận án xây dựng trên nền tảng lý thuyết biến dạng cắt bậc nhất (FSDT), trong đó các ứng suất trượt ngang qua chiều dày được giả thiết phân bố không đổi và hiệu chỉnh qua hệ số $\kappa = 5/6$. Với các vỏ composite rất dày ($a/h < 5$), mô hình chưa mô tả được đầy đủ hiệu ứng biến dạng phi tuyến qua chiều dày như lý thuyết HSDT hoặc lý thuyết từng lớp (Layerwise).
  2. Tương tác sóng nổ ngầm được mô hình hóa dưới dạng trường áp lực quy ước tác dụng lên bề mặt kết cấu theo công thức Cole, chưa xét đến tương tác hai chiều đầy đủ giữa chất lỏng và kết cấu (Acoustic Fluid-Structure Interaction - FSI) có kể đến hiện tượng xâm thực (cavitation).
  3. Vật liệu composite được giả thiết làm việc trong miền đàn hồi tuyến tính phi tuyến hình học, chưa xét đến tổn hại tích lũy, phá hủy tách lớp (delamination) hoặc chảy dẻo kim loại trong vật liệu FGM.

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

  • Hướng 1: Mở rộng phần tử ES-MITC3 kết hợp lý thuyết biến dạng cắt bậc cao (HSDT) hoặc lý thuyết Zigzag cải tiến để phân tích chính xác ứng suất cắt cục bộ giữa các lớp vỏ composite dày.
  • Hướng 2: Tích hợp mô hình phá hủy lũy tiến (Progressive Damage Modeling - Hashin/Puck failure criteria) và cơ học nứt vỡ (XFEM/Phase-field) trên nền tảng ES-FEM để dự báo sự phá hủy kết cấu dưới tải nổ.
  • Hướng 3: Phát triển bài toán ghép FSI hai chiều đồng thời giữa ES-FEM và phương pháp hạt mịn trơn (Smoothed Particle Hydrodynamics - SPH) để mô phỏng trọn vẹn hiện tượng nổ ngầm kèm tạo bọt khí xâm thực va đập.
  • Hướng 4: Tối ưu hóa đa mục tiêu góc xếp lớp sợi composite và hàm phân bố vật liệu FGM nhằm cực tiểu hóa khối lượng và cực đại hóa khả năng kháng chấn sóng xung kích.

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

  • Tác động học thuật (Academic Impact): Tạo ra tiền đề học thuật vững chắc trong cộng đồng cơ học tính toán kết cấu tại Việt Nam và khu vực. Phương pháp ES-MITC3 mở ra nhánh nghiên cứu mới về phần tử trơn cho kết cấu thành mỏng, với tiềm năng trích dẫn cao trong các tạp chí quốc tế hàng đầu ISI/Scopus (như Composite Structures, Thin-Walled Structures, Computers & Structures).
  • Chuyển đổi công nghệ công nghiệp (Industry Transformation): Ứng dụng trực tiếp trong ngành đóng tàu cao tốc, hàng không vũ trụ và chế tạo bồn bể áp lực composite. Giúp rút ngắn chu kỳ thử nghiệm ảo (virtual testing) từ vài tháng xuống vài giờ tính toán trên máy tính cá nhân.
  • Ý nghĩa quốc phòng - an ninh: Nâng cao năng lực tự chủ tính toán thiết kế vỏ khí tài quân sự, thân tàu ngầm mini, ống phóng ngư lôi và công trình phòng thủ ven biển chịu tải trọng nổ xung kích trong môi trường nước.

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

  • Nghiên cứu sinh & Giảng viên ngành Cơ kỹ thuật / Kỹ thuật Xây dựng / Hàng không: Tiếp cận mã nguồn thuật toán, phương pháp luận phát triển phần tử hữu hạn trơn tiên tiến và kỹ thuật xử lý phi tuyến hình học phức tạp.
  • Kỹ sư R&D tại các Viện nghiên cứu và Doanh nghiệp đóng tàu: Sử dụng quy trình và dữ liệu tham số để thiết kế tối ưu kết cấu vỏ composite nhẹ, độ bền cao, chống chịu va đập nổ ngầm.
  • Các cơ quan quản lý và hoạch định tiêu chuẩn quân sự: Khai thác các phát hiện về đáp ứng động lực học nổ ngầm để xây dựng quy chuẩn đánh giá sức bền sống còn của trang thiết bị ngầm.

Câu hỏi chuyên sâu

1. Đóng góp lý thuyết độc đáo nhất của luận án là gì và đã mở rộng lý thuyết nào?

Đóng góp lý thuyết độc đáo nhất là việc thiết lập thành công mô hình phần tử vỏ phẳng tam giác ES-MITC3 bằng cách tích hợp toán tử làm trơn biến dạng trên cạnh (Edge-based SFEM) vào kỹ thuật nội suy ten-xơ hỗn hợp (MITC3). Luận án đã mở rộng lý thuyết tấm/vỏ Reissner-Mindlin (FSDT) và lý thuyết phần tử trơn của G. R. Liu sang bài toán vỏ thoải đa độ cong composite và FGM, loại bỏ hoàn toàn hiện tượng khóa cắt mà không cần bổ sung các hàm nội suy bậc cao tốn kém chi phí tính toán.

2. Sự đổi mới về phương pháp luận so với các nghiên cứu quốc tế tương đương?

So với phần tử MITC3 cổ điển của Lee & Bathe (2004) vốn có xu hướng cứng hơn thực tế trên lưới thưa, và so với phần tử trơn CS-DSG3 của Cui et al. (2011) vốn nhạy cảm với hình dạng miền phân chia tế bào, ES-MITC3 tạo ra ma trận độ cứng có năng lượng biến dạng tối ưu tiệm cận nghiệm chính xác liên tục. Đồng thời, việc kết hợp đồng thời ba thuật toán: ES-MITC3 + Arc-length phi tuyến + Newmark ẩn tạo nên một khung phương pháp luận hoàn chỉnh và mạnh mẽ cho cả bài toán tĩnh, mất ổn định giới hạn và động lực học nổ ngầm.

3. Phát hiện bất ngờ và có ý nghĩa cơ học sâu sắc nhất là gì?

Phát hiện nổi bật là sự biến đổi độ cong Gauss ($K = 1/(R_x R_y)$) tạo ra sự thay đổi căn bản trong cơ chế chịu tải của vỏ dưới áp lực sóng nổ ngầm: Vỏ yên ngựa ($K < 0$, hyperbolic paraboloid) phân tán ứng suất dạng sóng biến dạng uốn phân tán nhanh hơn vỏ trụ ($K = 0$), nhưng lại chịu biên độ dao động cục bộ lớn hơn tại các góc tự do. Bên cạnh đó, việc điều chỉnh góc đặt cốt sợi $[45^\circ/-45^\circ/0^\circ]$ giúp giảm ứng suất Von-Mises đỉnh tới $28.4%$ so với cấu trúc trực giao $[0^\circ/90^\circ/0^\circ]$ dưới cùng một tải trọng nổ ngầm.

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

Hoàn toàn có. Toàn bộ công thức giải tích ma trận chuyển đổi tọa độ $\boldsymbol{\Lambda}$, các ma trận chuyển đổi biến dạng $\boldsymbol{\Lambda}{m1}^k, \boldsymbol{\Lambda}{b1}^k, \boldsymbol{\Lambda}_{s1}^k$, công thức tính ma trận độ cứng trơn $\tilde{\mathbf{K}}_k$, sơ đồ khối thuật toán lặp Arc-length Newton-Raphson và sơ đồ tích phân Newmark động lực học đều được chi tiết hóa từng bước trong luận án, cho phép tái lập độc lập 100% trên MATLAB, Fortran hoặc C++.

5. Định hướng chương trình nghiên cứu 10 năm tiếp theo được vạch ra như thế nào?

Lộ trình 10 năm tập trung vào: (1) Nâng cấp ES-MITC3 lên phần tử 3D suy biến phi tuyến vật liệu dẻo - từ biến; (2) Xây dựng mô-đun mô phỏng phá hủy tách lớp đa tỷ lệ (multiscale delamination) cho composite cấu trúc phức tạp; (3) Hoàn thiện phần mềm CAE chuyên dụng tính toán an toàn nổ ngầm cho các kết cấu công trình biển và phương tiện tác chiến ngầm tự hành (UUV).

Kết luận

  1. Phát triển thành công phần tử vỏ tam giác phẳng mới ES-MITC3, kết hợp hoàn hảo kỹ thuật làm trơn trên cạnh và phương pháp nội suy ten-xơ hỗn hợp, triệt tiêu hoàn toàn hiện tượng khóa cắt trên mọi tỷ lệ chiều dày ($a/h$ từ $10$ đến $1000$).
  2. Thiết lập khung phân tích phi tuyến hình học toàn diện bằng cách tích hợp ES-MITC3 với phương pháp Arc-length và Newton-Raphson, mô phỏng xuất sắc các điểm uốn mất ổn định snap-through và snap-back trên vỏ trụ, vỏ cầu và vỏ FGM.
  3. Lần đầu tiên giải quyết bài toán động lực học phi tuyến của vỏ composite chịu sóng xung kích nổ ngầm bằng phương pháp phần tử hữu hạn trơn, làm sáng tỏ định lượng ảnh hưởng của khối lượng nổ TNT ($M$), khoảng cách nổ ($S$), chiều dày và góc xếp lớp sợi.
  4. Xây dựng bộ chương trình tính toán số độc lập hoàn chỉnh trên MATLAB, đạt độ chính xác cao và độ hội tụ vượt trội khi đối sánh với nghiệm giải tích Navier và phần mềm thương mại.
  5. Mở ra 3 hướng nghiên cứu học thuật mới: Cơ học tính toán vỏ trơn bậc cao (HSDT-SFEM), tương tác chất lỏng - kết cấu nổ ngầm (FSI-SFEM) và tối ưu hóa kháng chấn cho kết cấu composite trong môi trường khắc nghiệt.
  6. Đóng góp giá trị thực tiễn quan trọng cho nền khoa học kỹ thuật quân sự và công nghiệp đóng tàu hiện đại, khẳng định tính tiên phong và chiều sâu học thuật của một công trình nghiên cứu tiến sĩ xuất sắc.