Luận án TS Trần Thị Thanh Vân: Vỏ composite tròn xoay chịu áp lực trong, công nghệ quấn

Tổng quan về luận án

Nghiên cứu về kết cấu vỏ mỏng tròn xoay chế tạo từ vật liệu composite cốt sợi liên tục chịu áp lực trong đóng vai trò cốt lõi trong kỹ thuật hàng không vũ trụ, công nghiệp quốc phòng và hệ thống lưu trữ năng lượng hiện đại. Vỏ chịu áp lực truyền thống chế tạo từ hợp kim kim loại độ bền cao (như thép Cr-Ni-Mo, thép không gỉ austenit hoặc hợp kim titan Ti-6Al-4V) bộc lộ những giới hạn nghiêm trọng về tỷ trọng lớn, nguy cơ tập trung ứng suất nhiệt tại mối hàn và hiện tượng suy giảm độ dai va đập tại vùng chuyển tiếp liên kết. Trong bối cảnh đó, vật liệu composite cốt sợi thủy tinh, sợi aramid hoặc sợi cacbon kết hợp nền polymer nhiệt rắn mang lại bước nhảy vọt nhờ độ bền riêng ($R_m/\rho$) và mô đun đàn hồi riêng ($E/\rho$) vượt trội. Minh chứng lịch sử từ việc ứng dụng vỏ động cơ tên lửa nhiên liệu rắn Polaris A2 (1960) hay động cơ tên lửa "Titan-III" đã chứng minh khả năng cắt giảm khối lượng kết cấu từ 40% đến 80% so với kết cấu kim loại, nâng cao đáng kể hệ số hoàn thiện kết cấu $W = p_b V / m$ (với $p_b$ là áp suất nổ phá hủy, $V$ là dung tích bên trong và $m$ là khối lượng tổng thể).

Khoảng trống học thuật (Research Gap) then chốt mà luận án của NCS. Trần Thị Thanh Vân (2021) tập trung giải quyết xuất phát từ tính tương tác phức tạp của mối quan hệ ba chiều "Vật liệu – Kết cấu – Công nghệ". Trong công nghệ chế tạo bình áp lực composite, sơ đồ quấn phẳng (Polar Winding) – một trường hợp riêng của quấn phi trắc địa (Non-geodesic Winding) – có ưu thế vượt trội về năng suất nhờ đầu quấn chuyển động quay tròn liên tục quanh phôi mà không cần cơ cấu đảo chiều đổi hướng phức tạp như quấn xoắn (Helical Winding). Tuy nhiên, các công trình kinh điển của Hartung (1963) và Vydrin cùng cộng sự (1978) đã chỉ ra rằng đường cong kinh tuyến của đáy vỏ trụ kín quấn phẳng luôn xuất hiện điểm uốn hình học, khiến quỹ đạo đáy không thể hội tụ về lỗ cực nếu không có giải pháp hiệu chỉnh. Đồng thời, nghiên cứu của Wang và cộng sự (2011) khi xấp xỉ đáy bằng hai cung tròn đã làm xuất hiện bước nhảy đột ngột của hệ số trượt sợi, gây mất tính liên tục cơ học. Cho đến trước công trình này, cơ sở giải tích xác định miền giới hạn hình học tường minh ($\bar{e}, \alpha$) nhằm triệt tiêu hoàn toàn nguy cơ trượt sợi ($\lambda \le [\lambda]$) cho bình trụ kín quấn phẳng vẫn chưa được thiết lập hoàn chỉnh.

Để giải quyết triệt để khoảng trống trên, luận án thiết lập 4 câu hỏi nghiên cứu (Research Questions) và 4 giả thuyết tương ứng:

• RQ1: Làm thế nào để thiết lập mô hình toán tổng quát xác định biên dạng đáy cân bằng của vỏ trụ composite quấn phẳng nhằm khắc phục hiện tượng điểm uốn hình học?
  Giả thuyết H1: Việc tích hợp điều kiện nắp đóng kín ($q = p r_p / 2$) cho phần đường cong cơ sở và điều kiện lỗ cực hở ($q = 0$) tại vùng lân cận cực ($r_f \ge r_i$) sẽ tạo ra đường cong kinh tuyến đáy liên tục, hội tụ trơn tru về cực.
• RQ2: Mối quan hệ hàm số giữa các tham số hình học ban đầu (độ lệch tâm không thứ nguyên $\bar{e}$, góc quấn ban đầu $\alpha$) và hệ số trượt sợi $\lambda$ dọc theo trục tọa độ $\bar{z}$ diễn ra theo quy luật nào?
  Giả thuyết H2: Phân bố hệ số trượt cực đại $\lambda_{max}$ phụ thuộc phi tuyến vào cặp tham số $(\bar{e}, \alpha)$, cho phép thiết lập bản đồ đường đồng mức xác định miền tham số an toàn không trượt.
• RQ3: Giới hạn hình học tỷ lệ chiều dài trên đường kính ($L/D < 2$) trong quấn phẳng có bản chất cơ học giải tích như thế nào?
  Giả thuyết H3: Giới hạn $L/D$ là hệ quả trực tiếp của điều kiện cân bằng lực ma sát tiếp tuyến - pháp tuyến và sức căng sợi dọc trục trên phần thân trụ.
• RQ4: Mức độ tương thích giữa mô hình giải tích dựa trên lý thuyết lưới (Netting Theory) và kết quả thực nghiệm chế tạo vỏ bình 1 lít chịu áp lực 150 at là bao nhiêu?
  Giả thuyết H4: Khi độ bền kéo dọc sợi vượt trội so với nền polyme, sai số giữa áp suất nổ thực nghiệm và áp suất phá hủy lý thuyết sẽ nằm trong giới hạn kiểm soát dưới 10%.

Khung lý thuyết của luận án tích hợp Lý thuyết lưới mở rộng cho quỹ đạo phi trắc địa kết hợp Hình học vi phân mặt cong tròn xoay và Cơ học tiếp xúc ma sát Coulomb. Đóng góp đột phá của công trình được định lượng thông qua việc xác định chính xác miền thông số công nghệ khả thi với hệ số trượt cho phép $[\lambda] \in [0,1; 0,4]$, chế tạo thành công mẫu bình composite sợi thủy tinh E/epoxy biến tính dung tích $V = 1\text{ lít}$, chiều dày tính toán tối ưu, đạt áp suất phá hủy nổ $p_b \ge 150\text{ at}$ ($15\text{ MPa}$), mở ra lộ trình tự chủ công nghệ quấn bình áp lực cao trong công nghiệp dân sự và quân sự.

Literature Review và Positioning

Lịch sử nghiên cứu thiết kế và chế tạo vỏ composite tròn xoay chịu áp lực trong trên thế giới phát triển qua ba dòng học thuật (Research Streams) chính:

[Dòng 1: Lý thuyết lưới & Quấn trắc địa]
Read (1963), Vasiliev (2009)
(Định luật Clairaut: r.sin β = const; Triệt tiêu lực tiếp tuyến: λ = 0)
                      │
                      ▼
[Dòng 2: Quấn phi trắc địa & Cơ học trượt sợi]
Hartung (1963), Vydrin (1978), Bunakov (1982), Wang (2011)
(Cân bằng vi phân sợi trên mặt cong; Điều kiện ổn định ma sát: λ ≤ [λ])
                      │
                      ▼
[Dòng 3: Cơ học môi trường liên tục dị hướng]
De-Jong (1983), Fukunaga & Uemura (1983), Zu (2010, 2014)
(Tiêu chuẩn phá hủy Tsai-Hill/Tsai-Wu; Cực tiểu hóa năng lượng biến dạng U)

Dòng thứ nhất bắt đầu từ công trình nền tảng của Read (1963) về lý thuyết lưới, đặt giả thiết toàn bộ tải trọng áp lực trong do hệ cốt sợi đơn hướng đảm nhận, còn nhựa nền chỉ đóng vai trò truyền ứng suất cắt. Quỹ đạo sợi được rải theo đường trắc địa tuân thủ định luật Clairaut: $$r \sin \beta = \text{const}$$ ở đó $r$ là bán kính hướng tâm và $\beta$ là góc quấn hợp bởi tiếp tuyến sợi và tiếp tuyến kinh tuyến. Vasiliev (2009) đã hệ thống hóa các phương trình vi phân biên dạng đáy cho bình có lỗ cực hở và lỗ cực đóng kín, đồng thời chứng minh rằng quấn trắc địa triệt tiêu hoàn toàn lực cắt ngang ($f_t = 0, \lambda = 0$), mang lại độ ổn định hình học tuyệt đối khi quấn nhưng lại hạn chế nghiêm ngặt không gian hình học đáy và gây phình dày cục bộ tại vùng cực.

Dòng thứ hai tập trung vào bài toán quấn phi trắc địa nhằm tạo ra các hình học đáy linh hoạt. Khi sợi bị làm lệch khỏi đường trắc địa, xuất hiện lực tiếp tuyến $f_t$ có xu hướng làm trượt sợi khỏi bề mặt khuôn quấn. Để sợi không bị trượt, tỷ số giữa lực tiếp tuyến và lực pháp tuyến $f_n$ phải thỏa mãn điều kiện ổn định ma sát: $$\lambda = \left|\frac{f_t}{f_n}\right| \le [\lambda]$$ với $[\lambda]$ là hệ số trượt cho phép (tương đương hệ số ma sát tĩnh $\mu$ giữa sợi tẩm nhựa và bề mặt khuôn). Hartung (1963) áp dụng lý thuyết cân bằng màng để thiết kế đáy quấn phẳng nhưng gặp phải hiện tượng uốn cong đường kinh tuyến trước khi tới trục quay mà không đưa ra giải pháp hiệu chỉnh. Bunakov và cộng sự (1982) đề xuất thuật toán hiệu chỉnh hai vùng nhưng vướng phải bài toán lặp hội tụ phức tạp do tham số bán kính cực $r_p$ không được cố định trước. Wang và cộng sự (2011) cố gắng khắc phục bằng mô hình hình học ghép hai cung tròn, song lại tạo ra điểm kỳ dị gián đoạn về hệ số trượt tại vị trí chuyển tiếp cung tròn, gây suy giảm nghiêm trọng độ bền chịu tải cục bộ.

Dòng thứ ba sử dụng Lý thuyết cơ học vật liệu composite dị hướng nhiều lớp kết hợp các tiêu chuẩn phá hủy như Tsai-Hill, Tsai-Wu hoặc cực tiểu hóa năng lượng biến dạng (De-Jong, 1983; Zu và cộng sự, 2010, 2014). Mặc dù mô hình dị hướng mô tả chính xác trạng thái ứng suất phẳng đa trục của vật liệu, nó đòi hỏi phải xác định trước 5 hằng số đàn hồi độc lập ($E_1, E_2, \nu_{12}, G_{12}, \dots$) và giải hệ phương trình vi phân phi tuyến phức tạp bậc cao.

Tranh luận học thuật cốt lõi tồn tại giữa quan điểm của Lý thuyết lưới (coi mô đun ngang $E_2 = 0$) và Lý thuyết môi trường liên tục dị hướng ($E_2 \neq 0$). Phân tích dữ liệu thực nghiệm cho thấy với các hệ vật liệu composite độ bền cao có mức độ dị hướng vượt trội như sợi thủy tinh/epoxy ($E_1 = 60\text{ GPa}, E_2 = 13\text{ GPa}$), sợi carbon/epoxy ($E_1 = 140\text{ GPa}, E_2 = 11\text{ GPa}$) hoặc sợi aramid/epoxy ($E_1 = 95\text{ GPa}, E_2 = 5,1\text{ GPa}$), tham số vật liệu $k = \frac{E_2 (1 - \nu_{21})}{E_1 (1 - \nu_{12})}$ rất nhỏ ($k \approx 0,07 - 0,16$). Khi $k \to 0$, mô hình cơ học dị hướng tiệm cận hoàn toàn về mô hình lý thuyết lưới. Do đó, việc sử dụng lý thuyết lưới mở rộng để tối ưu hóa kết cấu vỏ quấn phẳng không chỉ đảm bảo độ chính xác kỹ thuật mà còn đơn giản hóa cấu trúc giải tích của mô hình thiết kế.

Luận án định vị chính xác điểm giao thoa học thuật: Kế thừa tính tường minh của lý thuyết lưới, phát triển mô hình giải tích quấn phẳng giải quyết triệt để bài toán điểm uốn hình học bằng kỹ thuật hiệu chỉnh bán kính trơn $r_f = 1,5 r_p$, và lần đầu tiên lượng hóa thành công không gian tham số $(\bar{e}, \alpha, \bar{L}_{max})$ dưới ràng buộc khắt khe của hệ số trượt $\lambda \le [\lambda]$.

Đó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 Lý thuyết lưới cổ điển của Read (1963) sang bài toán quấn phẳng phi trắc địa có xét đến ma sát biên giới hạt. Trên một phân tố diện tích vỏ tròn xoay chịu áp suất trong $p$, phương trình cân bằng màng liên kết nội lực kinh tuyến $N_1$ và nội lực vĩ tuyến $N_2$ với các bán kính cong chính $R_1$ và $R_2$: $$\frac{N_1}{R_1} + \frac{N_2}{R_2} = p$$ Theo mô hình mạng lưới sợi chịu kéo thuần túy với lực kéo dọc trục $F$, các thành phần nội lực được phân rã theo góc quấn $\beta$: $$N_1 = F \cos^2 \beta, \quad N_2 = F \sin^2 \beta$$ Kết hợp với phương trình cân bằng vi phân sợi trên bề mặt có độ cong pháp tuyến $k_n$ và độ cong trắc địa $k_g$, hệ số trượt sợi không thứ nguyên $\lambda$ được thiết lập dưới dạng giải tích: $$\lambda = \frac{k_g}{k_n \sin \beta + \cos \beta \frac{d\phi}{ds_m}}$$ ở đó $s_m$ là chiều dài cung kinh tuyến và $\phi$ là góc nghiêng tiếp tuyến kinh tuyến.

                    ┌──────────────────────────────────────────────┐
                    │      Phương trình cân bằng màng vỏ mỏng      │
                    │         N₁/R₁ + N₂/R₂ = p ; N₂/N₁ = tan²β    │
                    └──────────────────────┬───────────────────────┘
                                           │
                                           ▼
                    ┌──────────────────────────────────────────────┐
                    │    Hình học vi phân & Cơ học trượt sợi       │
                    │   λ = |ft / fn| = |kg / (kn sinβ + ...)|     │
                    └──────────────────────┬───────────────────────┘
                                           │
                                           ▼
                    ┌──────────────────────────────────────────────┐
                    │      Mô hình quấn phẳng tham số hóa          │
                    │    r̄ = r/R, z̄ = z/R, ē = e/R, α = góc ban đầu │
                    └──────────────────────┬───────────────────────┘
                                           │
                                           ▼
                    ┌──────────────────────────────────────────────┐
                    │   Hiệu chỉnh điểm uốn biên dạng đáy          │
                    │   Vùng 1 (r̄ ≥ r̄f): q = p·rp/2 (Đáy đóng)     │
                    │   Vùng 2 (r̄ < r̄f): q = 0       (Đáy hở)      │
                    └──────────────────────────────────────────────┘

Mô hình lý thuyết đưa ra bước chuyển đổi hệ hình (Paradigm Shift) khi chuyển từ việc giải phương trình vi phân quỹ đạo sợi đơn lẻ sang khảo sát trường tương quan đồng mức giữa độ lệch tâm không thứ nguyên $\bar{e} = e/R$ và góc quấn ban đầu $\alpha$. Luận án chứng minh rằng biên dạng đáy không bị suy biến khi áp dụng giải pháp phân vùng tải trọng: vùng cơ sở ($\bar{r} \ge \bar{r}_f$) ứng với nắp đóng kín chịu lực dọc trục phân bố $q = p r_p / 2$, và vùng chuyển tiếp ($\bar{r} < \bar{r}_f$) ứng với điều kiện mép tự do $q = 0$ tại bán kính hiệu chỉnh $\bar{r}_f \ge 1,5 \bar{r}_p$.

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

Khung phân tích của luận án được xây dựng từ sự tích hợp của ba trụ cột lý thuyết:

1. Lý thuyết hình học vi phân mặt cong tròn xoay: Mô tả mối quan hệ giữa bán kính cực $\bar{r}_p$, bán kính trụ $R$, và các độ cong kinh tuyến $R_1$, vĩ tuyến $R_2$.
2. Lý thuyết lưới sợi đơn hướng: Thiết lập quan hệ phân bố chiều dày lớp composite trên phần trụ $h_c$, tại xích đạo $h_{eq}$, và phân bố liên tục trên nắp đáy $h(\bar{z})$: $$h(\bar{z}) = \frac{h_{eq}}{\bar{r} \cos \beta \sqrt{1 + (r')^2}}$$
3. Cơ học ma sát tiếp xúc phi trắc địa Coulomb: Xác lập ranh giới ổn định động học của sợi gia cường trên bề mặt khuôn tẩm nhựa ướt.

Điều kiện biên (Boundary Conditions) được xác định tường minh:

• Tại xích đạo đáy bình ($\bar{z} = 0, \bar{r} = 1$): Góc quấn đạt cực tiểu $\beta = \beta_{eq}$, chiều dày đạt giá trị danh nghĩa $h_{eq}$.
• Tại điểm uốn hình học $\bar{r}_i$: Đạo hàm bậc hai $r''(\bar{z}) = 0$, chuyển đổi trạng thái tiếp xúc sang vùng hiệu chỉnh.
• Tại lỗ cực ($\bar{r} = \bar{r}_p$): Quỹ đạo sợi tiếp tuyến với đường kính bích kim loại, lực kéo dọc trục sợi triệt tiêu tải trọng áp lực tác dụng lên nắp đóng kín.

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ủ lập trường nhận thức luận Thực chứng kết hợp Hiện thực phê phán (Positivism & Critical Realism), sử dụng thiết kế đa mức (Multi-level Design) kết hợp chặt chẽ giữa mô hình hóa toán học giải tích, mô phỏng số tích phân Runge-Kutta và thực nghiệm kiểm chứng phá hủy áp lực thủy tĩnh.

[MỨC 1: Vi mô / Trung mô] ──► Xác định cơ tính đơn lớp cốt sợi E-Glass/Epoxy
                               (Thử kéo mẫu vòng Split-Disk: σ_C = 1501 MPa, E_C = 47 GPa)
                                      │
                                      ▼
[MỨC 2: Vĩ mô / Động học] ──► Tích phân số hệ phương trình vi phân biên dạng & trượt sợi
                               (Giải thuật Runge-Kutta trên Matlab/Mathematica: ē, α, λ_max)
                                      │
                                      ▼
[MỨC 3: Cấu trúc hệ thống] ──► Chế tạo máy quấn phẳng & Thử nghiệm nổ bình 1L
                               (Hydrostatic Burst Test: Đo áp suất nổ p_b, kiểm chứng sai số)

Quy mô mẫu thử nghiệm bao gồm:

• Mẫu xác định cơ tính vật liệu: Các mẫu vòng composite quấn đơn hướng chế tạo trên khuôn quấn chuyên dụng, thử kéo theo tiêu chuẩn ASTM D2290 (Nolan split-disk method) với $N = 10$ mẫu để đảm bảo độ tin cậy thống kê.
• Mẫu cấu trúc bình hoàn chỉnh: Bình composite hình trụ có hai đáy kín, dung tích $V = 1\text{ lít}$, đường kính $D = 80\text{ mm}$, chiều dày tính toán $h_c \approx 1,2\text{ mm}$, được quấn phẳng ướt và đóng rắn nhiệt.

Quy trình nghiên cứu rigorous

Quy trình nghiên cứu thực nghiệm và thu thập dữ liệu được kiểm soát nghiêm ngặt qua các công đoạn:

1. Phối chế hệ vật liệu nền: Nhựa epoxy biến tính được kiểm soát độ nhớt động học trong dải tối ưu cho quá trình quấn ướt, duy trì độ ổn định độ nhớt lên đến $48\text{ giờ}$ ở dải nhiệt độ $25 - 30^\circ\text{C}$, đảm bảo khả năng thấm ướt kiệt sợi thủy tinh E-glass roving trước khi dải lên khuôn.
2. Hệ thống tạo sức căng và chia sợi: Sợi từ các bôbin được dẫn qua hệ thống bạc dẫn hướng, lược chia sợi và khối tạo sức căng sơ bộ chủ động, đi qua lô tẩm nhựa nóng nhằm triệt tiêu bọt khí và rỗ xốp.
3. Chế tạo thiết bị quấn phẳng chuyên dụng: Thiết kế và chế tạo máy quấn phẳng vận hành theo nguyên lý đầu quấn quay quanh trục phôi với cơ cấu dẫn động điều khiển số, đảm bảo góc quấn $\alpha$ và độ lệch tâm $e$ chính xác theo tọa độ cài đặt.
4. Quy trình thử nghiệm nổ thủy tĩnh (Hydrostatic Burst Testing): Mẫu bình sau đóng rắn được nạp môi chất lỏng (nước), gia áp cưỡng bức bằng bơm cao áp thủy lực với tốc độ tăng áp ổn định $0,5\text{ MPa/s}$, trang bị cảm biến áp suất kỹ thuật số độ chính xác cao và hệ thống ghi dữ liệu thời gian thực cho đến thời điểm phá hủy hoàn toàn.

Tam giác đạc phương pháp luận (Triangulation) được thiết lập vững chắc giữa: (1) Nghiệm giải tích vi phân màng $\leftrightarrow$ (2) Mô phỏng tích phân số $\leftrightarrow$ (3) Dữ liệu biến dạng và áp suất phá hủy thực nghiệm.

Data và phân tích

Đặc trưng cơ lý tính của vật liệu composite đơn hướng sợi thủy tinh E/epoxy biến tính xác định từ thực nghiệm thử kéo mẫu vòng đạt các giá trị cụ thể:

• Độ bền kéo dọc trục sợi: $[\sigma_c] = 1501\text{ MPa}$
• Mô đun đàn hồi dọc trục: $E_C = 47\text{ GPa}$
• Độ dãn dài tương đối khi đứt: $\epsilon_C = 2,8%$
• Tỷ trọng vật liệu: $\rho = 2,053\text{ g/cm}^3$

Hệ thống thuật toán số tích phân đường cong kinh tuyến và trường hệ số trượt được lập trình trên phần mềm tính toán kỹ thuật chuyên dụng. Khảo sát tham số được thực hiện toàn diện với $\bar{e}$ biến thiên từ $0,0$ đến $0,5$ và góc $\alpha$ từ $\pi/20$ đến $\pi/4$. Kết quả kiểm tra độ bền vững (Robustness Checks) khẳng định nghiệm số hội tụ ổn định, không phát sinh dao động số tại lân cận điểm uốn sau khi áp dụng thuật toán hiệu chỉnh bán kính $r_f$.

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

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

Nghiên cứu mang lại 4 phát hiện khoa học mang tính đột phá:

                            ┌────────────────────────────────────────────────────────┐
                            │              4 PHÁT HIỆN ĐỘT PHÁ CHÍNH                 │
                            └───────────────────────────┬────────────────────────────┘
                                                        │
         ┌───────────────────────┬──────────────────────┴────────────────┬────────────────────────┐
         ▼                       ▼                                       ▼                        ▼
┌──────────────────┐   ┌──────────────────┐                    ┌──────────────────┐     ┌──────────────────┐
│  Phát hiện 1     │   │   Phát hiện 2    │                    │   Phát hiện 3    │     │   Phát hiện 4    │
│Xác lập bán kính  │   │Bản đồ đường đồng │                    │Giải mã giới hạn  │     │Thực nghiệm nổ    │
│hiệu chỉnh tối ưu │   │mức (ē, α) - λ_max│                    │hình học L/D < 2  │     │bình 1L đạt 150 at│
│r_f = 1,5 r_p     │   │loại bỏ trượt sợi │                    │dưới góc độ cơ học│     │sai số lý thuyết  │
│khử uốn biên dạng │   │với λ ≤ [λ]       │                    │cân bằng ma sát   │     │dưới 6%           │
└──────────────────┘   └──────────────────┘                    └──────────────────┘     └──────────────────┘

1. Khắc phục hoàn toàn hiện tượng uốn biên dạng đáy: Luận án chỉ ra rằng điểm uốn hình học $r_i$ là thuộc tính tự nhiên của phương trình vi phân nắp đóng kín. Bằng cách thiết lập điểm bắt đầu đường cong hiệu chỉnh tại $r_f = 1,5 r_p$, đường cong kinh tuyến đáy chuyển đổi trơn tru sang mô hình mép tự do ($q = 0$), loại bỏ hoàn toàn hiện tượng phân kỳ toán học mà Hartung (1963) chưa xử lý được.
2. Xác lập bản đồ đường đồng mức tham số $(\bar{e}, \alpha) - \lambda_{max}$: Đã phát hiện quy luật phân bố phi tuyến của hệ số trượt cực đại $\lambda_{max}$. Khi độ lệch tâm $\bar{e}$ tăng từ $0,0$ lên $0,5$, phạm vi góc quấn ban đầu $\alpha$ cho phép để sợi không trượt ($\lambda \le [\lambda]$) bị thu hẹp đáng kể. Cụ thể, với hệ số ma sát thực tế $[\lambda] = 0,2$, việc lựa chọn cặp tham số $(\bar{e} = 0,2; \alpha = \pi/20)$ đảm bảo phân bố hệ số trượt $\lambda(\bar{z})$ luôn nằm sâu trong miền an toàn ($\lambda \le 0,15$).
3. Giải mã bản chất cơ học của giới hạn $L/D < 2$: Nghiên cứu đã chứng minh bằng giải tích rằng quy tắc kinh nghiệm $L/D < 2$ trong quấn phẳng (từng được Eckold, 1994 và Peters, 2011 nêu ra) thực chất xuất phát từ giới hạn tích lũy góc xoắn sợi trên phần thân trụ. Khi chiều dài phần trụ không thứ nguyên $\bar{L} = L/R$ vượt quá ngưỡng giới hạn $\bar{L}_{max}(\bar{e}, \alpha)$, lực kéo dọc sợi tạo ra thành phần lực tiếp tuyến vượt quá giới hạn ma sát nghỉ $[\lambda]$, dẫn đến hiện tượng trượt sợi hàng loạt trên thân trụ.
4. Kiểm chứng thực nghiệm áp suất nổ chính xác: Mẫu bình 1 lít chế tạo bằng công nghệ quấn phẳng ướt đạt áp suất phá hủy thực tế $p_b = 150\text{ at}$ ($15\text{ MPa}$). Phân tích dạng phá hủy cho thấy toàn bộ kết cấu bị phá hủy đồng thời do đứt sợi dọc trục tại phần thân trụ và nắp đáy, không xuất hiện hiện tượng tách lớp sớm hay trượt bóc sợi tại vùng cực, chứng minh tính tối ưu hóa ứng suất đồng đều của biên dạng thiết kế với sai số so với mô hình lý thuyết dưới $6%$.

Implications đa chiều

• Về mặt học thuật: Công trình khẳng định tính đúng đắn và hiệu quả của Lý thuyết lưới mở rộng khi ứng dụng cho các hệ composite có độ bất đẳng hướng cơ học cao ($E_1/E_2 > 4,5$), cung cấp công cụ toán học tường minh thay thế cho các thuật toán giải lặp phức tạp.
• Về mặt phương pháp luận: Thiết lập phương pháp chuẩn hóa không thứ nguyên ($\bar{r}, \bar{z}, \bar{e}, \bar{L}$) kết hợp biểu đồ đồng mức, cho phép các kỹ sư thiết kế nhanh chóng tra cứu thông số kết cấu - công nghệ tối ưu mà không cần chạy lại toàn bộ quy trình mô phỏng phần tử hữu hạn đa tỷ lệ.
• Về mặt ứng dụng thực tiễn: Công nghệ quấn phẳng cho phép phát triển các dòng máy quấn kết cấu đơn giản, vận hành ở vận tốc góc cao và liên tục, rút ngắn chu kỳ gia công bình áp lực dung tích nhỏ ($< 5\text{ lít}$) xuống $30 - 40%$ so với máy quấn xoắn 4 trục hoặc 5 trục CNC đắt tiền.
• Về mặt chính sách và công nghiệp quốc phòng: Tạo nền tảng làm chủ độc lập công nghệ chế tạo vỏ động cơ tên lửa cỡ nhỏ, bình khí tài dưỡng khí cứu hộ, bình chứa khí nén cao áp cho phương tiện không người lái (UAV) và bình chứa nhiên liệu khí hóa lỏng CNG/LPG, giảm thiểu phụ thuộc vào trang thiết bị nhập khẩu.

Limitations và Future Research

Mặc dù đạt được những kết quả đột phá, nghiên cứu vẫn tồn tại một số giới hạn học thuật cần được thừa nhận khách quan:

1. Giả thiết màng mỏng bỏ qua ứng suất uốn cục bộ: Mô hình lý thuyết lưới thuần túy giả định trạng thái ứng suất màng phân bố đều qua chiều dày, chưa tính đến ảnh hưởng của mômen uốn cục bộ và ứng suất tiếp cắt lớp tại vùng chuyển tiếp hình học giữa thân trụ và nắp đáy ($r = R$).
2. Mô hình ma sát đẳng hướng đơn giản hóa: Giả thiết hệ số trượt cho phép $[\lambda]$ là hằng số ma sát tĩnh đẳng hướng chưa phản ánh đầy đủ bản chất phi tuyến phụ thuộc vào nhiệt độ đóng rắn, độ nhớt biến thiên của nhựa nền epoxy và vận tốc rải sợi động trong quá trình quấn.
3. Phạm vi thực nghiệm giới hạn ở quy mô phòng thí nghiệm: Kiểm chứng thực nghiệm mới dừng lại ở mẫu bình dung tích $1\text{ lít}$, chịu tải trọng áp lực tĩnh đơn kỳ ở nhiệt độ phòng, chưa khảo sát độ bền mỏi chu kỳ cao (áp suất nén - xả tuần hoàn) và tác động môi trường nhiệt độ khắc nghiệt.

Chương trình nghiên cứu tương lai (Future Research Agenda) cần tập trung vào các hướng:

• Mở rộng mô hình sang bài toán vỏ dày nhiều lớp có xét đến hiệu ứng cắt lớp ngang và suy giảm độ cứng suy rộng (Degradation continuum damage mechanics).
• Phát triển hệ thống điều khiển tự động thích nghi sức căng sợi theo thời gian thực (Closed-loop real-time tension control) nhằm tối ưu hóa độ nén ép lớp và mật độ sợi thể tích $V_f > 65%$.
• Nghiên cứu ứng dụng sợi cacbon cường độ siêu cao (T700/T800) kết hợp vỏ làm kín lót polyme nhiệt dẻo HDPE chế tạo bình áp lực Type IV siêu nhẹ cho hệ thống lưu trữ nhiên liệu Hydro áp suất cực cao ($350 - 700\text{ bar}$).

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

Nghiên cứu tạo tiền đề vững chắc cho việc thiết lập tiêu chuẩn thiết kế và chế tạo bình chịu áp composite tại Việt Nam:

┌───────────────────────────────────────────────────────────────────────────┐
│                    MA TRẬN TÁC ĐỘNG VÀ ẢNH HƯỞNG                          │
├──────────────────────┬────────────────────────────────────────────────────┤
│ Lĩnh vực             │ Kết quả lượng hóa & Tác động thực tiễn             │
├──────────────────────┼────────────────────────────────────────────────────┤
│ Học thuật & Nghiên cứu│ Đóng góp hệ phương trình giải tích biên dạng đáy;  │
│                      │ Cung cấp phương pháp luận chuẩn hóa không thứ      │
│                      │ nguyên cho các nghiên cứu sinh cơ khí động lực.    │
├──────────────────────┼────────────────────────────────────────────────────┤
│ Công nghiệp Quốc phòng│ Chế tạo vỏ động cơ tên lửa nhiên liệu rắn;         │
│                      │ Cắt giảm 40-80% khối lượng so với vỏ kim loại,     │
│                      │ tăng cự ly bay và tải trọng chiến thuật.           │
├──────────────────────┼────────────────────────────────────────────────────┤
│ Công nghiệp Dân sinh │ Bình khí nén cứu hỏa, khí tài dưỡng khí thợ lặn,   │
│                      │ bình composite chứa LPG/CNG nhẹ hơn 75% so với kim │
│                      │ loại, loại bỏ nguy cơ nổ vỡ mảnh nguy hiểm.        │
├──────────────────────┼────────────────────────────────────────────────────┤
│ Tự chủ Công nghệ     │ Tự chủ chế tạo máy quấn phẳng tốc độ cao giá rẻ,   │
│                      │ tiết kiệm 40% chi phí đầu tư thiết bị CNC ngoại.   │
└──────────────────────┴────────────────────────────────────────────────────┘

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

• Nghiên cứu sinh và Giảng viên ngành Kỹ thuật Cơ khí - Động lực: Tiếp cận hệ thống cơ sở lý thuyết lưới phi trắc địa và mô hình toán học giải tích mẫu mực để phát triển các cấu trúc composite phức tạp.
• Kỹ sư R&D tại các Viện nghiên cứu và Nhà máy Quốc phòng: Ứng dụng trực tiếp quy trình tính toán kích thước hình học, góc quấn ban đầu $\alpha$, độ lệch tâm $e$, và phân bố chiều dày $h(z)$ trong thiết kế vỏ chịu lực khí tài quân sự.
• Doanh nghiệp chế tạo thiết bị áp lực và vật liệu mới: Sở hữu giải pháp công nghệ quấn phẳng chi phí thấp, năng suất cao để sản xuất các dòng bình tích áp thủy lực, bình lọc nước composite công nghiệp và bình chứa khí nén chất lượng cao.
• Cơ quan quản lý tiêu chuẩn kỹ thuật: Có cơ sở khoa học để xây dựng bộ tiêu chuẩn đánh giá an toàn, kiểm tra không phá hủy (NDT) và thử nghiệm thủy tĩnh cho các sản phẩm bình áp lực composite nội địa.

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 độc đáo nhất là việc mở rộng thành công Lý thuyết lưới (Netting Theory) cổ điển sang sơ đồ quấn phẳng phi trắc địa với mô hình hiệu chỉnh biên dạng đáy hai vùng trơn ($r_f = 1,5 r_p$). Công trình đã khắc phục triệt để nghịch lý điểm uốn hình học của Hartung (1963) và xử lý dứt điểm điểm kỳ dị bước nhảy hệ số trượt trong mô hình hai cung tròn của Wang và cộng sự (2011), tạo ra hệ phương trình giải tích thuần nhất xác định biên dạng đáy cân bằng hoàn chỉnh cho bình composite kín.

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

So sánh với nghiên cứu của Bunakov và cộng sự (1982) phải sử dụng thuật toán giải lặp tham số bán kính cực $r_p$ phức tạp và Wang và cộng sự (2011) dựa trên hình học ghép cung tròn gượng ép, luận án đề xuất phương pháp chuẩn hóa tham số không thứ nguyên toàn phần ($\bar{r}, \bar{z}, \bar{e}, \bar{L}$). Từ đó, luận án thiết lập biểu đồ đường đồng mức xác định miền giới hạn hình học tường minh $(\bar{e}, \alpha) \to \lambda_{max}$, cho phép tối ưu hóa kết cấu trực tiếp mà không cần giải lặp số đa biến.

3. Phát hiện thực nghiệm bất ngờ nhất được hỗ trợ bởi dữ liệu là gì?

Phát hiện bất ngờ nhất là tính chất suy giảm đột ngột của miền tham số góc quấn ban đầu khả thi $\alpha$ khi độ lệch tâm không thứ nguyên $\bar{e}$ vượt quá giá trị $0,3$. Dữ liệu số chỉ ra rằng tại $\bar{e} = 0,5$, giá trị $\lambda_{max}$ tăng vọt lên trên $0,4$ ngay cả ở góc quấn nhỏ $\alpha = \pi/20$, giải thích chính xác nguyên nhân gây trượt sợi hàng loạt trong các thử nghiệm công nghệ trước đây khi cố gắng tăng độ lệch tâm để rút ngắn chiều dài nắp đáy.

4. Quy trình lặp lại thực nghiệm (Replication Protocol) có được cung cấp đầy đủ không?

Quy trình được cung cấp đầy đủ với độ chi tiết cao:

• Công thức đơn liệu nhựa nền epoxy biến tính ổn định độ nhớt 48 giờ ở $25 - 30^\circ\text{C}$.
• Quy trình quấn mẫu vòng và đồ gá thử kéo split-disk theo chuẩn xác định độ bền kéo $[\sigma_c] = 1501\text{ MPa}$ và $E_C = 47\text{ GPa}$.
• Sơ đồ động học máy quấn phẳng và lưu đồ thuật toán tính toán 2 bước xác định thông số công nghệ cho bình dung tích $1\text{ lít}$, áp suất nổ $150\text{ at}$.

5. Chương trình nghị sự 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 3 trọng tâm:

1. Nâng cấp mô hình giải tích sang vật liệu composite nhiệt dẻo (Thermoplastic composites) kết hợp công nghệ quấn dán nhiệt tại chỗ (In-situ consolidation).
2. Phát triển bình composite Type IV và Type V không lót kim loại phục vụ mạng lưới lưu trữ và phân phối Hydro lỏng/khí nén $700\text{ bar}$.
3. Ứng dụng trí tuệ nhân tạo (AI) và mạng nơ-ron vật lý (PINNs) để dự đoán quá trình tích tụ hư hỏng mỏi đa trục và dự báo tuổi thọ còn lại của bình chịu áp trong điều kiện vận hành thực tế.

Kết luận

Luận án tiến sĩ kỹ thuật của NCS. Trần Thị Thanh Vân đã hoàn thành xuất sắc các mục tiêu nghiên cứu đề ra với 5 đóng góp học thuật và thực tiễn cốt lõi:

1. Thiết lập thành công hệ phương trình toán học giải tích xác định biên dạng đáy cân bằng của vỏ composite hình trụ kín chế tạo theo sơ đồ quấn phẳng, giải quyết dứt điểm hiện tượng điểm uốn hình học bằng giải pháp hiệu chỉnh bán kính trơn $r_f = 1,5 r_p$.
2. Xây dựng mô hình quan hệ giải tích giữa các tham số hình học ban đầu ($\bar{e}, \alpha$) với trường phân bố hệ số trượt sợi $\lambda(\bar{z})$, lượng hóa chính xác biểu đồ đường đồng mức giới hạn không trượt sợi ứng với các cấp ma sát $[\lambda] \in [0,1; 0,4]$.
3. Luận giải khoa học bản chất cơ học của giới hạn tỷ lệ chiều dài trên đường kính $L/D < 2$ trong công nghệ quấn phẳng dựa trên điều kiện cân bằng ma sát tiếp xúc.
4. Chế tạo thành công thiết bị quấn phẳng chuyên dụng quy mô phòng thí nghiệm và làm chủ quy trình công nghệ quấn ướt sử dụng hệ vật liệu sợi thủy tinh E/epoxy biến tính có độ bền kéo đạt $1501\text{ MPa}$.
5. Kiểm chứng thực nghiệm thành công trên mẫu bình composite 1 lít, đạt áp suất nổ phá hủy $150\text{ at}$ với độ sai lệch so với mô hình lý thuyết màng mỏng dưới $6%$.

Công trình mở ra 3 hướng nghiên cứu mới về cơ học kết cấu composite: tối ưu hóa đa mục tiêu vỏ chịu áp phi trắc địa, phát triển bình áp lực cao Type IV/V cho năng lượng Hydro, và tự động hóa điều khiển thích nghi quá trình rải sợi composite trong nền công nghiệp sản xuất thông minh.