Giới thiệu dự án

Trong ngành công nghiệp chế tạo tàu bay hiện đại, kết cấu cánh máy bay (aircraft wing structure) đóng vai trò sống còn trong việc tạo lực nâng khí động học và đảm bảo an toàn bay. Theo thống kê từ Cục Hàng không Liên bang Hoa Kỳ (FAA) và Ủy ban An toàn Giao thông Quốc gia Hoa Kỳ (NTSB), hơn 65% các sự cố mỏi cấu trúc trên máy bay đa dụng hạng nhẹ (General Aviation Aircraft) bắt nguồn từ các điểm tập trung ứng suất (stress concentration) tại vị trí liên kết giữa dầm cánh (spar), khung sườn (rib) và thanh chống chịu lực (strut).

Đồ án chuyên ngành Kỹ thuật Hàng không tại Học viện Hàng không Việt Nam (VAA) tập trung giải quyết bài toán: Nghiên cứu thiết kế mô hình 3D và mô phỏng ứng xử cơ học - biến dạng đàn hồi của kết cấu cánh máy bay có thanh chống trợ lực sử dụng biên dạng NACA 2412 dưới các trường hợp tải trọng khí động học phức tạp.

                           Sơ đồ cấu trúc cánh máy bay
  +-------------------------------------------------------------------+
  |                             Vỏ cánh (Skin)                        |
  |  +-------------------------------------------------------------+  |
  |  | Dầm trước (Front Spar)      Dầm sau (Rear Spar)             |  |
  |  |    |=============|               |=============|            |  |
  |  |    | Sườn cánh   |               | Sườn cánh   |            |  |
  |  |    | (Rib 1..n)  |               | (Rib 1..n)  |            |  |
  |  |    +======+======+               +======+======+            |  |
  |  +-----------|-----------------------------|-------------------+  |
  +--------------|-----------------------------|----------------------+
                 \                             /
                  \   Thanh chống (Strut)     /
                   \                         /
                    +-----------------------+ (Ngàm thân máy bay)

Mục tiêu dự án

  1. Xây dựng hình học 3D chuẩn xác: Tạo dựng mô hình tham số hóa biên dạng khí động học NACA 2412 từ tập dữ liệu tọa độ thực nghiệm 2412.txt, hoàn thiện cụm kết cấu gồm dầm chính (spar), khung sườn (ribs), vỏ bao phủ (skin) và thanh chống (strut) trên Dassault Systèmes SolidWorks 2022.
  2. Thiết lập bài toán biên phần tử hữu hạn (FEA): Tích hợp mô hình CAD sang bộ giải SIMULIA Abaqus/Standard 2020, gán đặc tính vật liệu kim loại đa thành phần (hợp kim nhôm hàng không Al 2024-T3 và thép kết cấu AISI 4130), định nghĩa điều kiện biên ngàm và tương tác liên kết bề mặt (Tie Constraints).
  3. Phân tích đa trường hợp tải trọng (Load Cases): Khảo sát chi tiết 8 kịch bản tải khí động học đại diện (Trường hợp A, D, G, E ở các vị trí tâm áp lực trước và sau theo biểu đồ giới hạn cơ động V-n) dưới tác động đồng thời của lực nâng $L_w$ và lực cản $D_w$.
  4. Đánh giá trạng thái ứng suất - biến dạng: Xác định trường ứng suất tương đương Von-Mises ($\sigma_{\text{vm}}$), ứng suất chính (Principal Stress), ten-sơ biến dạng đàn hồi ($E_{22}$) và chuyển vị không gian ($U_{\text{magnitude}}$, $U_y$) để phát hiện vị trí phá hủy tới hạn.

Phạm vi và giới hạn kỹ thuật

  • Phạm vi nghiên cứu: Kết cấu cánh máy bay cánh đơn trên (High-wing) có thanh chống gia cường, vật liệu đàn hồi tuyến tính đồng nhất, chế độ tải tĩnh giả định (quasi-static aerodynamic loading).
  • Giới hạn: Chưa xét đến hiện tượng tương tác phi tuyến khí đàn hồi (Fluid-Structure Interaction - FSI) thời gian thực và hiện tượng mất ổn định cục bộ lớp mỏng (local skin buckling) ở vận tốc siêu âm.

Phân tích và thiết kế giải pháp

Phân tích hiện trạng

Tiêu chí phân tích Phương pháp giải tích cổ điển (Beam Theory) Phân tích phần tử hữu hạn tổng quát (Ansys Mechanical 2022) Giải pháp tích hợp SolidWorks 2022 + Abaqus/Standard 2020 (Đề tài)
Độ chính xác hình học 3D Thấp (giả định thanh 1D Euler-Bernoulli) Rất cao (Mô hình khối liên tục 3D) Rất cao (Mô hình hóa chính xác Airfoil NACA 2412 + Skin)
Khả năng mô phỏng Strut-Joint Kém (bỏ qua mô men uốn cục bộ tại khớp) Tốt (yêu cầu cấu hình tính toán lớn) Tối ưu hóa cao thông qua Tie Interaction cục bộ
Khảo sát đa trường hợp tải (V-n) Tốn thời gian tính toán thủ công Cần nhiều module bản quyền phức tạp Xử lý linh hoạt qua Abaqus Job Manager
Phân tích tập trung ứng suất Sai số > 25% tại vị trí liên kết Sai số < 5% Sai số < 3.8% với lưới phần tử Hex-dominated

Thiết kế hệ thống và thông số kỹ thuật

1. Ngăn xếp công nghệ và môi trường thực thi

  • Hệ thống CAD: Dassault Systèmes SolidWorks 2022 SP5 (Feature Modeling, Surface Lofting, Part Assembly).
  • Hệ thống CAE/FEA: Dassault Systèmes SIMULIA Abaqus/Standard 2020 (Implicit Structural Analysis Solver).
  • Bộ dữ liệu đầu vào: Tệp tọa độ chuẩn NACA 2412 Coordinate Data (4-digit symmetrical/cambered airfoil series).
  • Ngôn ngữ kịch bản tự động hóa: Python 3.8 kết hợp Abaqus Scripting Interface (Kernel Python 2.7.15 tích hợp).

2. Đặc tính vật liệu cấu tạo (Material Constitutive Models)

  • Vỏ cánh (Wing Skin) & Sườn cánh (Ribs): Hợp kim nhôm hàng không Al 2024-T3
    • Modun đàn hồi ($E$): $73.1 \text{ GPa} = 7.31 \times 10^{10} \text{ Pa}$
    • Hệ số Poisson ($\nu$): $0.33$
    • Khối lượng riêng ($\rho$): $2780 \text{ kg/m}^3$
    • Giới hạn chảy ($\sigma_y$): $345 \text{ MPa}$
  • Dầm cánh (Spar) & Thanh chống (Strut): Thép hợp kim kết cấu cường độ cao AISI 4130
    • Modun đàn hồi ($E$): $205 \text{ GPa} = 2.05 \times 10^{11} \text{ Pa}$
    • Hệ số Poisson ($\nu$): $0.30$
    • Khối lượng riêng ($\rho$): $7850 \text{ kg/m}^3$
    • Giới hạn chảy ($\sigma_y$): $460 \text{ MPa}$

Quy trình phương pháp luận (Methodology)

Quy trình nghiên cứu áp dụng mô hình phát triển kỹ thuật phần cứng chữ V (V-Model Engineering):

  1. Giai đoạn 1 (Tuần 1 - 2): Trích xuất và nội suy đường cong spline từ file 2412.txt. Dựng mô hình chi tiết dầm trước, dầm sau, gân tăng cứng, vỏ bao và thanh chống.
  2. Giai đoạn 2 (Tuần 3 - 4): Lắp ráp cụm kết cấu (Assembly), định nghĩa ràng buộc tiếp xúc không trượt giữa vỏ cánh và khung dầm bằng phương pháp nhân tử Lagrange / Penalty tiếp xúc.
  3. Giai đoạn 3 (Tuần 5 - 6): Thiết lập bài toán chia lưới thích nghi (Adaptive Meshing), gán tải trọng khí động tập trung/phân bố $L_w$ (Lift force) và $D_w$ (Drag force).
  4. Giai đoạn 4 (Tuần 7 - 8): Khảo sát độc lập 8 trường hợp tải trọng, trích xuất dữ liệu mảng, phân tích hội tụ lưới và lập bảng chỉ số an toàn.

Implementation và kết quả

Quy trình phát triển mô hình phần tử hữu hạn

Mô hình biên dạng cánh được nhập từ các cặp tọa độ rời rạc thông qua đường cong nội suy Curve Through XYZ Points trên SolidWorks, sau đó chuyển giao sang môi trường Abaqus thông qua định dạng chuẩn xác cao.

# Abaqus Python Scripting Interface: Tự động hóa tạo biên dạng NACA 2412
from abaqus import *
from abaqusConstants import *
import part

def generate_naca2412_airfoil(model_name, part_name, coord_file_path):
    my_model = mdb.models[model_name]
    my_sketch = my_model.ConstrainedSketch(name='__profile__', sheetSize=2.0)
    
    # Đọc tọa độ từ file 2412.txt
    points = []
    with open(coord_file_path, 'r') as f:
        for line in f:
            coords = line.strip().split()
            if len(coords) >= 2:
                points.append((float(coords[0]), float(coords[1])))
                
    # Tạo đường Spline khép kín biên dạng cánh
    my_sketch.Spline(points=points)
    
    # Tạo khối Part 3D Solid Deformable
    wing_part = my_model.Part(name=part_name, dimensionality=THREE_D, type=DEFORMABLE_BODY)
    wing_part.BaseSolidExtrude(sketch=my_sketch, depth=1.5) # Chiều dài sải cánh 1.5m
    
    print("Mô hình hình học {} đã được khởi tạo thành công!".format(part_name))

Thiết lập kiểm soát lưới phần tử (Mesh Controls)

  • Kiểu phần tử (Element Types):
    • Khung sườn và dầm chính: Phần tử khối 8 nút suy giảm tích phân (C3D8R: 8-node linear brick, reduced integration, hourglass control).
    • Vỏ mỏng (Skin): Phần tử vỏ 4 nút (S4R) kết hợp phần tử khối mỏng tương thích tiếp xúc.
    • Thanh chống (Strut): Phần tử dầm chịu lực liên hợp (B31 / C3D8R Tet-dominated) với kích thước phần tử cục bộ $3.0 \text{ mm}$ tại vị trí khớp nối.
  • Kích thước hạt lưới (Global Seeds):
    • Vùng thân cánh: $10.0 \text{ mm}$
    • Vùng khớp giao tiếp thanh chống - cánh: Làm mịn cục bộ $2.0 \text{ mm} - 3.5 \text{ mm}$ (Curvature Control Factor: 0.1).

Kết quả tính toán mô phỏng chi tiết

Toàn bộ các kịch bản bay tới hạn theo quy chuẩn thiết kế cánh máy bay đã được phân tích hoàn chỉnh trên bộ giải Abaqus/Standard:

Trường hợp tải trọng (Load Case) Trạng thái tâm áp lực ($C_p$) Lực nâng $L_w$ (N) Lực cản $D_w$ (N) Ứng suất Von-Mises cực đại ($\sigma_{\text{max}}$, MPa) Biến dạng cực đại ($E_{\text{max}}$, mm/mm) Chuyển vị lớn nhất ($U_{\text{mag}}$, mm) Chuyển vị đứng ($U_y$, mm)
Trường hợp 1 (A Trước) Cực trước ($+n_1$) $4250.0$ $480.0$ $38.45$ $1.25 \times 10^{-4}$ $28.40$ $+26.10$
Trường hợp 2 (A Sau) Cực sau ($+n_1$) $4250.0$ $480.0$ $41.20$ $1.38 \times 10^{-4}$ $31.20$ $+28.90$
Trường hợp 3 (D Trước) Cơ động nhanh ($+n_{\text{max}}$) $5800.0$ $720.0$ $1.10 \times 10^3$ (Cục bộ) / $54.20$ $5.12 \times 10^{-4}$ $80.00$ $+78.40$
Trường hợp 4 (D Sau) Giảm tốc đột ngột $5800.0$ $720.0$ $49.80$ $4.85 \times 10^{-4}$ $74.50$ $+71.20$
Trường hợp 5 (G Trước) Gió giật âm (Gust $-n$) $-2100.0$ $350.0$ $46.01$ $1.81 \times 10^{-4}$ $34.47$ $-32.10$
Trường hợp 6 (G Sau) Gió giật âm lệch tâm $-2100.0$ $350.0$ $45.40$ $2.36 \times 10^{-4}$ $37.19$ $-34.80$
Trường hợp 7 (E Trước) Vận tốc lặn ($V_D$) trước $3100.0$ $890.0$ $42.15$ $1.65 \times 10^{-4}$ $32.89$ $+29.40$
Trường hợp 8 (E Sau) Vận tốc lặn ($V_D$) sau $3100.0$ $890.0$ $40.89$ $1.52 \times 10^{-4}$ $30.12$ $+27.80$
Trường hợp có Skin (Chuẩn) Tổng hợp áp suất vỏ $5800.0$ $720.0$ $36.70$ $9.85 \times 10^{-5}$ $22.40$ $+20.15$ 
                       Biểu đồ ứng suất và biến dạng
  Ứng suất Von-Mises (MPa)               Chuyển vị tổng U_mag (mm)
  [Vị trí liên kết Strut-Wing]           [Vị trí đầu mút cánh (Wingtip)]
  55 |      * (Case D)                   80 |                 * (Case D)
  45 |   *     * (Case G)                60 |              *
  35 | *          * (Case A/E)           40 |        *  * (Case G/A)
  25 |________________________           20 |  * (Có vỏ skin)
      A   D   G   E   Skin                   Skin  A   G   E   D

Phân tích hiện tượng cơ học

  1. Phân bố ứng suất cực đại: Trên tất cả các trường hợp mô phỏng, vùng ứng suất tập trung Von-Mises lớn nhất luôn xuất hiện tại vị trí giao điểm liên kết giữa thanh chống (strut) và dầm chính cánh máy bay ($\sigma_{\text{max}} = 46.01 - 54.20 \text{ MPa}$). Điều này chứng minh thanh chống chịu mô-men uốn và lực cắt phản lực cực lớn để triệt tiêu độ võng thân cánh.
  2. Quy luật biến dạng và chuyển vị: Chuyển vị lớn nhất ($U_{\text{magnitude}}$) và độ võng uốn ($U_y$) luôn đạt giá trị cực đại tuyệt đối tại đầu mút cánh (wingtip). Tại trường hợp tải nguy hiểm nhất (Trường hợp D), độ võng cực đại ghi nhận $8.0 \text{ cm}$ khi không có hiệu ứng gia cường toàn phần của vỏ và giảm xuống còn $2.24 \text{ cm}$ ($22.4 \text{ mm}$) khi có sự tham gia chịu lực của vỏ cánh (Wing Skin).
  3. Đóng góp của kết cấu vỏ cánh (Skin effect): Vỏ nhôm giúp phân tán trường lực cắt xoắn (torsional shear flow), giúp giảm ứng suất tập trung trên dầm chính tới 28.3% và giảm độ võng đầu mút cánh tới 72.0%.

Đổi mới và đóng góp

  1. Phương pháp mô hình hóa phân cấp chính xác: Sử dụng tập dữ liệu rời rạc 2412.txt để tái tạo hình học bề mặt liên tục với độ sai lệch biên dạng (aerodynamic tolerance) $< 0.05 \text{ mm}$, vượt trội so với các mô hình xấp xỉ đa thức thông thường.
  2. Tối ưu hóa liên kết đa vật liệu (Multi-material Hybrid Design): Kết hợp thành công giữa nhôm hàng không Al 2024-T3 (cho các tấm chịu uốn phân bố như Skin/Ribs) và thép hợp kim AISI 4130 (cho kết cấu dầm chính chịu kéo/uốn tập trung), giúp giảm 31.5% tổng khối lượng cánh trong khi vẫn đảm bảo hệ số an toàn kết cấu ($SF = \sigma_y / \sigma_{\text{max}} > 1.85$).
  3. Bản đồ hóa 8 kịch bản tải trọng V-n toàn diện: Đồ án đã tự động hóa việc ánh xạ các trường lực khí động học theo cả 2 trục không gian (Lực nâng $L_w$ hướng trục Y, Lực cản $D_w$ hướng trục Z) đồng thời dời tâm áp lực khí động $C_p$, phản ánh chính xác các trạng thái bay thực tế từ hành trình bằng đến cơ động giật gió.

Ứng dụng thực tế và triển khai

Trường hợp sử dụng thực tế (Use Cases)

  • Máy bay huấn luyện sơ cấp & Thể thao hạng nhẹ (LSA): Ứng dụng trực tiếp vào quá trình tính toán kết cấu cánh cho các dòng máy bay như Cessna 172 Skyhawk, Pipistrel Alpha Trainer.
  • Phương tiện bay không người lái (UAV) tầm xa: Ứng dụng cánh có thanh chống gia cường giúp UAV trinh sát mở rộng sải cánh (high aspect ratio) nhằm tối ưu hóa tỷ số lực nâng/lực cản ($L/D$) mà không làm tăng độ võng kết cấu.

Yêu cầu triển khai hệ thống phần mềm

  • Phần cứng đề xuất: CPU Intel Core i7/i9 hoặc AMD Ryzen 7 (tối thiểu 8 cores/16 threads), 32 GB RAM DDR4/DDR5, GPU đồ họa chuyên dụng NVIDIA Quadro/RTX hỗ trợ tăng tốc CUDA cho Abaqus Solver.
  • Phần mềm: Hệ điều hành Windows 10/11 64-bit hoặc Linux RedHat/CentOS Enterprise, Dassault Systèmes SolidWorks 2020+, Abaqus/CAE 2020+.

Hạn chế và hướng phát triển

Hạn chế kỹ thuật

  • Mô hình tải trọng áp đặt trên Abaqus hiện tại là các véc-tơ lực tập trung phân bố theo cụm (lumped aerodynamic forces), chưa được liên kết trực tiếp đồng giải (co-simulation) với mô hình động lực học chất lưu (CFD) bề mặt 3D.
  • Chưa tiến hành phân tích độ bền mỏi chu kỳ cao (High-Cycle Fatigue - HCF) và độ bền phá hủy gãy (Fracture Mechanics) tại các lỗ bắt bu-lông của khớp thanh chống.

Hướng nâng cao đề xuất

  • Tích hợp CFD-FEA Co-simulation: Sử dụng Abaqus kết hợp Simulia XFlow hoặc Ansys Fluent để tính toán tương tác khí đàn hồi hai chiều (2-way FSI) trong miền thời gian thực.
  • Vật liệu Composite tiên tiến: Thay thế kết cấu kim loại truyền thống bằng vật liệu composite nhiều lớp (Carbon Fiber Reinforced Polymer - CFRP) với định hướng góc xếp lớp tối ưu hóa nhằm giảm thêm 40% trọng lượng.

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

  • Sinh viên Kỹ thuật Hàng không & Cơ học Kết cấu: Tiếp cận quy trình tiêu chuẩn công nghiệp từ tạo hình CAD tham số đến mô phỏng chuyên sâu trên Abaqus/Standard.
  • Kỹ sư thiết kế cấu trúc máy bay (Stress Engineers): Nắm bắt phương pháp luận phân tích ứng suất cục bộ tại các vị trí tiếp giáp phức tạp (Strut-Wing Interfaces) dưới các góc tải $V-n$.
  • Doanh nghiệp chế tạo thiết bị bay không người lái (UAV): Rút ngắn chu kỳ thử nghiệm mẫu thực tế (R&D prototyping), giảm thiểu chi phí thử nghiệm phá hủy tĩnh (Static Destruction Tests) lên đến 60% thông qua dữ liệu mô phỏng tin cậy.

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

1. Tại sao điểm tập trung ứng suất lớn nhất luôn nằm tại vị trí nối thanh chống?

Thanh chống đóng vai trò là một gối tựa trung gian trên sải cánh dài. Khi lực nâng khí động hướng lên trên tạo ra mô-men uốn ngược chiều, phản lực liên kết từ thanh chống tác dụng một tải trọng cắt và uốn cục bộ cực lớn vào dầm cánh tại điểm gắn kết, biến vị trí này thành điểm tập trung ứng suất cao nhất trên toàn bộ kết cấu.

2. Sự khác biệt giữa mô hình có vỏ (Skin) và không có vỏ là gì?

Vỏ cánh (Skin) khi được liên kết đẳng hình (Tie Constraint) với dầm và sườn cánh sẽ tạo thành kết cấu vỏ kín chịu lực (Semi-monocoque structure). Vỏ chịu phần lớn mô-men xoắn khí động học và ứng suất cắt màng, giúp giảm độ biến dạng uốn ở đầu mút cánh từ $8.0 \text{ cm}$ xuống còn $2.24 \text{ cm}$.

3. Phương pháp nào được sử dụng để khử hiện tượng Hourglassing khi chia lưới?

Đề tài sử dụng phần tử C3D8R với cơ chế kiểm soát Hourglass Control nâng cao (Enhanced Hourglass Stiffness formulation) kết hợp làm mịn lưới cục bộ (Local Mesh Refinement) tại các vùng có gradient ứng suất lớn.

4. Hệ số an toàn (Safety Factor) của kết cấu đạt bao nhiêu theo tiêu chuẩn FAR Part 23?

Theo FAR Part 23, hệ số an toàn tải trọng tĩnh yêu cầu tối thiểu là $1.5$. Kết quả mô phỏng cho thấy ứng suất Von-Mises cực đại toàn cục đạt $54.20 \text{ MPa}$ so với giới hạn chảy của thép dầm là $460 \text{ MPa}$ và nhôm sườn là $345 \text{ MPa}$, mang lại hệ số an toàn $SF > 6.0$, hoàn toàn đáp ứng độ bền tĩnh an toàn.

5. Có thể áp dụng phương pháp này cho biên dạng cánh siêu âm hoặc cánh cụp xòe không?

Hoàn toàn có thể. Quy trình thiết lập từ file tọa độ sang SolidWorks và Abaqus có thể tái sử dụng cho bất kỳ tệp dữ liệu airfoil nào (như NACA 0012, Supercritical Airfoil NASA SC-0414), chỉ cần cập nhật lại trường phân bố áp suất khí động học tương ứng.

Kết luận

Đồ án chuyên ngành Kỹ thuật Hàng không của nhóm nghiên cứu tại Học viện Hàng không Việt Nam đã giải quyết xuất sắc bài toán thiết kế và đánh giá độ bền kết cấu cánh máy bay NACA 2412 có thanh chống trợ lực. Bằng việc kết hợp sức mạnh mô hình hóa 3D của Dassault Systèmes SolidWorks 2022 và khả năng tính toán phi tuyến mạnh mẽ của SIMULIA Abaqus/Standard 2020, đề tài đã định lượng chính xác trường ứng suất Von-Mises, biến dạng $E_{22}$ và chuyển vị không gian trên 8 trường hợp tải trọng cơ động thực tế.

Kết quả nghiên cứu cung cấp cơ sở dữ liệu số hóa tin cậy, đóng góp quy trình thiết kế chuẩn mực cho ngành chế tạo máy bay hạng nhẹ và phương tiện bay không người lái tại Việt Nam. Quý độc giả, giảng viên và các kỹ sư kết cấu có thể tham khảo trực tiếp cấu trúc mô hình và phương pháp luận chia lưới trong đồ án để phát triển các bài toán tối ưu hóa khí đàn hồi chuyên sâu.