Tổng quan về luận án

Thi công công trình ngầm và đường hầm khẩu độ vừa và nhỏ ($b \le 6\text{ m}, h \le 3\text{ m}$) trong điều kiện địa chất đá cứng có hệ số kiên cố $f_{kp} \ge 8$ bằng phương pháp khoan nổ là một bài toán kỹ thuật phức tạp hàng đầu trong lĩnh vực xây dựng công trình ngầm và công trình phòng thủ quân sự. Tại gương đào, công tác bốc xúc đá nổ mìn (ĐNM) là khâu then chốt chi phối trực tiếp đến hơn 40% chu kỳ thi công và năng suất tổng thể. Trong khi các tuyến hầm khẩu độ lớn ($35 - 70\text{ m}^2$) có thể áp dụng các dòng máy xúc lật đổ trước cỡ lớn như Komatsu WA100, WA350 hay Liugong CLG836 với dung tích gầu $1{,}0 - 3{,}5\text{ m}^3$, thì không gian chật hẹp của đường hầm khẩu độ nhỏ hoàn toàn không cho phép các thiết bị này quay đầu hoặc thao tác dỡ tải. Việc sử dụng các giải pháp thay thế như máy cào vơ, máy đào gầu ngược hoặc nhân công thủ công bộc lộ nhiều hạn chế nghiêm trọng về năng suất, khả năng xử lý đá quá cỡ và nguy cơ mất an toàn lao động trong môi trường thiếu dưỡng khí, nhiều khói bụi độc hại.

Khoảng trống nghiên cứu (research gap) cốt lõi xuất phát từ thực trạng: dòng máy xúc lật đổ bên (MXLĐB) bánh xích dẫn động điện chuyên dụng VMC E500-1 (dung tích gầu $0{,}5\text{ m}^3$) do Công ty Cổ phần Chế tạo máy Vinacumin chế tạo trước đây chủ yếu dựa trên phương pháp sao chép mẫu ngoại nhập và thử nghiệm mò mẫm thực nghiệm, hoàn toàn thiếu vắng cơ sở khoa học giải tích về tương tác cơ học tiếp xúc, động học và động lực học hệ nhiều vật (MBD). Chưa có một công trình nào mô hình hóa đầy đủ tương tác 3D giữa gầu xúc và môi trường ĐNM dạng hạt rời sắc cạnh trong điều kiện biên không gian hẹp, cũng như chưa xác lập được quy luật điều khiển động học ngược tối ưu hóa quỹ đạo chuyển động của các khâu dẫn.

Luận án của nghiên cứu sinh Tạ Văn Huy tại Học viện Kỹ thuật Quân sự (chuyên ngành Kỹ thuật cơ khí động lực, mã số 9 52 01 16, dưới sự hướng dẫn của TS. Chu Văn Đạt) tập trung giải quyết ba câu hỏi nghiên cứu (Research Questions - RQ) và kiểm chứng ba giả thuyết khoa học (Hypotheses - H) tương ứng:

  • RQ1: Bản chất quy luật tương tác và cơ chế hình thành các thành phần lực cản tiếp xúc giữa răng - gầu xúc với môi trường ĐNM biến thiên ngẫu nhiên sau khoan nổ là gì?
    H1: Tương tác gầu - ĐNM phụ thuộc phi tuyến vào độ sâu ấn gầu, hình học răng cắt và sự hình thành khối lăng trụ trượt; việc áp dụng phương pháp xúc phối hợp (kết hợp đồng thời tịnh tiến máy, nâng cần và quay gầu) sẽ làm giảm cực tiểu hóa lực cản đào tổng cộng so với phương pháp xúc phân đoạn truyền thống.
  • RQ2: Quy luật biến thiên tọa độ suy rộng và chuyển động của các khâu dẫn động thủy lực được xác lập như thế nào để đỉnh răng gầu bám chính xác quỹ đạo xúc và xả tải tối ưu trong không gian đường hầm bị giới hạn?
    H2: Bài toán động học ngược giải bằng phương pháp ma trận chuyển thuần nhất Denavit-Hartenberg sẽ xác định chính xác hành trình và vận tốc của 4 cụm xi lanh thủy lực dẫn động, triệt tiêu nguy cơ va chạm với nóc hầm ($h \le 3\text{ m}$) và tối ưu hóa hệ số điền đầy gầu $K_H$.
  • RQ3: Phổ tải trọng động và ứng suất tiếp xúc phát sinh tại các khớp bản lề quay ($O_2, O_3, O_5$) của thiết bị công tác (TBCT) diễn biến ra sao trong suốt chu trình làm việc?
    H3: Mô hình động lực học hệ nhiều vật dựa trên phương trình Lagrange loại 2 cho phép dự báo chính xác các đỉnh tải trọng động va đập và quá trình tiêu tán năng lượng, làm tiền đề tin cậy cho việc tối ưu hóa kết cấu và hệ thống thủy lực.

Khung lý thuyết của luận án được xây dựng dựa trên sự tích hợp của Cơ học môi trường hạt rời (Granular Medium Mechanics), Lý thuyết cơ học hệ nhiều vật (Multibody System Dynamics) và Lý thuyết điều khiển động học ngược robot công nghiệp. Phạm vi nghiên cứu bao gồm việc khảo nghiệm giải tích toàn diện mô hình máy VMC E500-1 kết hợp đo kiểm thực nghiệm tại mỏ than Hà Lầm và thực nghiệm hiện trường chuyên sâu, tạo nên bước đột phá khoa học giúp nâng cao 20-22% hiệu suất năng lượng và đảm bảo tính tự chủ công nghệ chế tạo trang bị công binh quốc phòng.

Literature Review và Positioning

Nghiên cứu về cơ học tiếp xúc giữa thiết bị công tác máy làm đất với môi trường vật liệu rời đã trải qua nhiều giai đoạn phát triển với những đóng góp quan trọng từ các học giả quốc tế và trong nước. Về quy luật phân bố thành phần cỡ hạt đá sau khoan nổ, trường phái nghiên cứu của Rôzin-Rămle, L. Ôxanhit, Gôđen-Anđrep và GS. Nhữ Văn Bách đã xác lập phương trình hàm mũ mô tả tỷ lệ cỡ hạt:

$$P(x) = \left(\frac{X}{X_{max}}\right)^m$$

Trong đó $X_{max}$ là kích thước cục đá lớn nhất và $m$ là hệ số thực nghiệm. Giáo sư B. Kutudôp tiếp tục hoàn thiện lý thuyết này khi thiết lập công thức xác định tỷ lệ đá quá cỡ $V_{qc} = V_0 [1 - 3(D_{cp}/4D_{tb})]%$, trong đó kích thước hạt cho phép lớn nhất $D_{cp}$ bị giới hạn bởi dung tích gầu theo điều kiện biên $D_{cp} = (0{,}7 \div 0{,}8)\sqrt[3]{V_g}$.

Về lực cản của môi trường tác dụng lên gầu xúc, tồn tại sự tranh luận học thuật sâu sắc giữa hai trường phái:

  • Trường phái giải tích đơn giản hóa: Tiêu biểu là S. Zhuravlev với mô hình lực cản $W = k \cdot l \cdot h \cdot g \cdot f + W_1$, và V. Stogov (1978) với công thức phân rã lực cản phương ngang $W_c = W_1 + W_2 + W_3$. Quan điểm này giả định đáy gầu luôn tịnh tiến song song với nền đống vật liệu và chỉ xét riêng cho phương pháp xúc phân đoạn.
  • Trường phái ứng suất cắt phi tuyến và cơ học môi trường rời liên tục: Đại diện bởi V. Tarasov dựa trên lý thuyết ứng suất giới hạn cắt, Kostylev (1988) với hàm mũ đa biến $P = a \cdot L_{BH}^n \cdot B^m \cdot C \cdot k_h \cdot k_a \cdot k_\varphi$, và V.A. Bauman (1984) với mô hình lăng trụ trượt nghiêng một góc cố định $\pi/4 - \varphi/2$.

Điểm yếu lớn nhất của các mô hình của V. Stogov và S. Zhuravlev là không mô tả được hiện tượng vật liệu nén chặt tương tác với thành sau của gầu xúc khi người vận hành ấn sâu gầu vào đống đá, đồng thời bỏ qua góc nghiêng thực tế của đáy gầu ($\alpha_0 = 7^\circ \pm 2^\circ$).

[Phân loại cỡ hạt ĐNM]
 (Gôđen-Anđrep, Kutudôp, Nhữ Văn Bách)
         │
         ▼
[Trường phái Lực cản Đơn giản hóa] ──Contradiction──► [Trường phái Ứng suất Cắt & Lăng trụ]
 (S. Zhuravlev, V. Stogov)                           (V. Tarasov, Kostylev, V.A. Bauman)
         │                                                    │
         └───────────────────────┬────────────────────────────┘
                                 ▼
                 [Mô hình hóa Động lực học MBD]
            (A. Hemani 1994, Ying Li & W. Liu 2012)
                                 │
                                 ▼
      [KHOẢNG TRỐNG: MXLĐB Bánh xích trong Hầm Khẩu độ Nhỏ]
                   (Luận án Tạ Văn Huy, 2020)

Về mặt mô hình hóa động học và động lực học, luận án định vị và so sánh trực tiếp với hai nghiên cứu quốc tế tiêu biểu:

  1. Nghiên cứu của A. Hemani (1994): Tác giả mô phỏng bộ công tác máy xúc lật đổ trước như một tay robot phẳng 3 bậc tự do (3-DOF) gồm 1 khớp trượt và 2 khớp quay, phân tích 5 thành phần lực cản $r_1, r_2, r_3, r_4, r_5$. Tuy nhiên, mô hình của Hemani chỉ áp dụng trong không gian mở không giới hạn chiều cao và bỏ qua chuyển động lật đổ bên theo phương ngang.
  2. Nghiên cứu của Ying Li và Wenyuan Liu (2012): Nhóm tác giả xây dựng sơ đồ động học 7 khâu với 4 bậc tự do trong phần mềm MSC ADAMS cho máy xúc lật bánh lốp. Mặc dù đạt độ chi tiết cao về hình học, mô hình này không tính toán đến cơ chế bám trượt của dải xích trên nền đá cứng theo lý thuyết của Wong J.Y và Ahmad Hemami, đồng thời không giải quyết bài toán động học ngược cho điều kiện đổ bên trong đường hầm hẹp.

Tại Việt Nam, các nghiên cứu của GS.TS. Chu Văn Đạt, GS.TS. Phan Nguyên Di và PGS.TS. Nguyễn Viết Tân đã ứng dụng lý thuyết hệ nhiều vật cho máy đào gầu ngược và máy bốc xúc cào vơ dạng cần. Luận án của Tạ Văn Huy đã tiên phong lấp đầy khoảng trống nghiên cứu khi tích hợp toàn diện bài toán tương tác hạt rời phi tuyến, động học ngược và động lực học hệ nhiều vật cho dòng máy xúc lật đổ bên bánh xích VMC E500-1 trong không gian đường hầm bị giới hạn nghiêm ngặt.

Đó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 phong phú thêm lý thuyết cơ học tiếp xúc của A.I. Balốpnhép và V.A. Bauman bằng cách thiết lập mô hình giải tích chính xác hóa lực cản xúc ĐNM trong không gian hẹp. Điểm cốt lõi là việc định lượng hóa ảnh hưởng của góc nghiêng đáy gầu $\alpha_0$, góc dốc tự nhiên của đống đá $\alpha_{VL}$ và góc ma sát trong $\varphi$ lên sự hình thành khối lăng trụ trượt phía trước và hai bên thành gầu.

Mô hình lý thuyết tổng quát biểu diễn quan hệ phi tuyến của các lực cản tiếp xúc:

$$F_{ct} = f(L_{AG}, B_{gx}, \rho_{vl}, \varphi, \alpha_0, \alpha_{VL})$$

Trong đó lực cản chính diện $F_{ct}$ và lực cản bên $F_{cb}$ được tích hợp đồng thời với lực nén ép tại đỉnh răng $F_E$ và lực ma sát trượt của đá trên thép $\mu_1$.

         Fcty (Lực cản đứng)
           ▲
           │        /│
           │       / │ Khối lăng trụ ĐNM
           │      /  │ (Góc trượt π/4 - φ/2)
           │     /   │
───────────┼────/────┼──────────────► Fctx (Lực cản ngang)
           │   / ┌───┘
           │  /  │ Đáy gầu (Góc nghiêng α0 = 7° ± 2°)
           │ /   │
           ▼     ▼
          Răng gầu (FE, ME)

Về mặt động lực học, công trình đã thiết lập thành công hệ phương trình vi phân mô tả chuyển động của cơ hệ nhiều vật dựa trên phương trình Lagrange loại 2 dạng tổng quát:

$$\frac{d}{dt}\left(\frac{\partial T}{\partial \dot{q}_i}\right) - \frac{\partial T}{\partial q_i} + \frac{\partial \Phi}{\partial \dot{q}_i} + \frac{\partial \Pi}{\partial q_i} = Q_i \quad (i = 1, 2, 3, 4)$$

Trong đó:

  • $q = [d_1, \theta_2, \theta_3, \theta_4]^T$ là véc-tơ tọa độ suy rộng biểu thị chuyển vị tịnh tiến của thân máy ($d_1$), góc nâng cần ($\theta_2$), góc quay gầu xúc ($\theta_3$), và góc lật nghiêng gầu đổ bên ($\theta_4$).
  • $T = \sum_{i=1}^4 T_i$ là tổng động năng của các khâu (thân máy, cần nâng, giá quay, gầu xúc và khối vật liệu $G_{VL}$).
  • $\Pi$ là thế năng đàn hồi và thế năng trọng trường của hệ.
  • $\Phi$ là hàm hao tán năng lượng Rayleigh biểu diễn tổn hao ma sát tại các khớp bản lề và cản nhớt thủy lực.
  • $Q_i$ là các lực suy rộng bao gồm lực đẩy của mô-tơ xích $P_k$, lực tác động của các xi-lanh thủy lực ($F_{xl1}, F_{xl2}, F_{xl3}, F_{xl4}$) và mô-men cản tương tác $M_E$.

Đây là bước chuyển dịch mang tính mô thức (paradigm shift) từ việc thiết kế thuần túy dựa trên công thức bán thực nghiệm sang khung lý thuyết toán - cơ giải tích có khả năng dự báo chính xác đáp ứng động lực học.

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

Khung phân tích của luận án tích hợp liên ngành ba trụ cột lý thuyết: Cơ học vật liệu rời dạng hạt, Động lực học hệ nhiều vật (MBD) và Thuật toán điều khiển động học ngược.

┌────────────────────────────────────────────────────────────────────────┐
│                        KHUNG PHÂN TÍCH TÍCH HỢP                        │
└───────────────────────────────────┬────────────────────────────────────┘
                                    │
    ┌───────────────────────────────┼───────────────────────────────┐
    ▼                               ▼                               ▼
┌───────────────────────┐ ┌───────────────────────┐ ┌───────────────────────┐
│ Cơ học Môi trường rời │ │   Động lực học MBD    │ │   Động học Ngược      │
│ (Lăng trụ trượt,      │ │ (Phương trình         │ │ (Ma trận D-H,         │
│ Phân bố cỡ hạt,       │ │ Lagrange loại 2, hàm  │ │ Quy luật khâu dẫn,    │
│ Lực cản tiếp xúc)     │ │ hao tán Rayleigh)     │ │ Tối ưu hóa quỹ đạo)   │
└───────────┬───────────┘ └───────────┬───────────┘ └───────────┬───────────┘
            │                         │                         │
            └─────────────────────────┼─────────────────────────┘
                                      ▼
┌────────────────────────────────────────────────────────────────────────┐
│     Mô hình Số hóa: Mô phỏng Động lực học Khớp & Tiêu hao Năng lượng   │
│                 (Phần mềm MATLAB & MSC ADAMS)                          │
└───────────────────────────────────┬────────────────────────────────────┘
                                    │
                                    ▼
┌────────────────────────────────────────────────────────────────────────┐
│     Kiểm chứng Thực nghiệm Hiện trường (Mỏ than Hà Lầm & DasyLab 11)   │
│         [Đo áp suất OCM-511, Lưu lượng R4S, Biến dạng Tenzo]           │
└────────────────────────────────────────────────────────────────────────┘

Phương pháp giải bài toán động học ngược sử dụng công cụ ma trận biến đổi thuần nhất Denavit-Hartenberg để thiết lập chuỗi động học từ hệ trục tọa độ cố định gắn trên thân máy $O_1X_1Y_1Z_1$ đến hệ tọa độ gắn tại đỉnh răng gầu $O_4X_4Y_4Z_4$. Từ quỹ đạo không gian tối ưu được hoạch định, thuật toán tính toán ngược ra quy luật hành trình $x_1(t), x_2(t), x_3(t), x_4(t)$ và vận tốc tức thời $v_1, v_2, v_3, v_4$ của từng xi lanh thủy lực.

Điều kiện biên của khung phân tích được xác lập chặt chẽ:

  • Chiều rộng đường hầm: $B_h \le 6{,}0\text{ m}$; Chiều cao nóc hầm: $H_h \le 3{,}0\text{ m}$.
  • Hệ số kiên cố của đá sau nổ: $f_{kp} = 8 \div 14$; Khối lượng thể tích đống đá: $\rho_M = 1{,}6\text{ tấn/m}^3$; Hệ số tơi xốp: $K_p = 1{,}1$.
  • Dung tích gầu danh nghĩa: $V_g = 0{,}5\text{ m}^3$; Tải trọng nâng danh nghĩa: $Q_N = 1000\text{ kg}$.

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

Thiết kế nghiên cứu

Luận án tuân thủ triết lý nghiên cứu thực chứng diễn dịch (Positivism / Deductive Empirical Realism), kết hợp chặt chẽ giữa suy luận mô hình hóa toán học thuần túy, mô phỏng số trên máy tính và đo kiểm thực chứng đa kênh tại hiện trường. Thiết kế nghiên cứu đa tầng (Multi-level design) được triển khai qua 3 cấp độ:

  1. Cấp độ vi mô: Phân tích ứng suất tiếp xúc cục bộ tại đỉnh răng gầu $\sigma_{rc}, \sigma_{pt}$ và tương tác ma sát khô giữa các hạt đá vỡ.
  2. Cấp độ trung mô: Khảo sát động lực học liên kết tại các khớp bản lề quay $O_2$ (khớp cần), $O_3$ (khớp quay gầu), $O_5$ (khớp lật nghiêng gầu) và tương tác bánh xích - nền đường hầm theo mô hình bám trượt Wong J.Y.
  3. Cấp độ vĩ mô: Đánh giá toàn bộ chu trình công tác $T_Ц$, năng suất kỹ thuật $\Pi_T$ và mức độ tiêu hao năng lượng riêng của máy xúc lật VMC E500-1 trong không gian đường hầm mô phỏng.

Quy trình nghiên cứu rigorous

Quy trình thu thập dữ liệu thực nghiệm được thiết kế và chuẩn hóa với các trang thiết bị đo lường kỹ thuật cao:

  • Đồng hồ đo tốc độ vòng quay: Thiết bị quang học HHT13 giám sát tốc độ động cơ điện chính.
  • Cảm biến lưu lượng thủy lực: Cảm biến tua-bin R4S-7HD-25 và R4S-7HD-50 gắn trực tiếp trên đường cấp dầu vào xi lanh nâng cần và mô-tơ di chuyển.
  • Cảm biến áp suất: Đầu đo áp suất piezoresistive OCM-511 (dải đo $0 - 60\text{ MPa}$, cấp chính xác $0{,}2%$) bố trí tại khoang đẩy và khoang kéo của các xi lanh $F_{xl2}, F_{xl3}, F_{xl4}$.
  • Cảm biến dịch chuyển và vận tốc: Cảm biến khoảng cách laser H7 và cảm biến đo vận tốc quang điện GSS25 đo chuyển vị tức thời của các cần piston xi lanh và dải xích.
  • Cảm biến đo biến dạng (Tenzo): Cụm Tenzo foil FCA-5-11 dán trực tiếp tại các vị trí tập trung ứng suất quanh trục khớp $O_2, O_3, O_5$ để đo biến dạng vi mô $\varepsilon$, từ đó chuyển đổi sang ứng suất và lực tác dụng thực tế.
  • Thu thập và xử lý tín hiệu: Tín hiệu tương tự từ các cảm biến được số hóa qua hai bộ chuyển đổi đa kênh NI USB-6003 và NI USB-6009 của National Instruments, kết nối máy tính điều khiển cài đặt phần mềm chuyên dụng DasyLab phiên bản 11.0 với tần số lấy mẫu đạt $1000\text{ Hz}$.
┌────────────────────────────────────────────────────────────────────────┐
│                     HỆ THỐNG THỰC NGHIỆM ĐA KÊNH                       │
└───────────────────────────────────┬────────────────────────────────────┘
                                    │
    ┌───────────────────────┬───────┴───────────────┬────────────────────┐
    ▼                       ▼                       ▼                    ▼
┌──────────────────┐ ┌──────────────────┐ ┌──────────────────┐ ┌──────────────────┐
│ Cảm biến Áp suất │ │Cảm biến Lưu lượng│ │  Cảm biến Tenzo  │ │ Cảm biến Laser/  │
│     OCM-511      │ │     R4S-7HD      │ │     FCA-5-11     │ │  Vận tốc GSS25   │
└─────────┬────────┘ └────────┬─────────┘ └────────┬─────────┘ └────────┬─────────┘
          │                   │                    │                    │
          └───────────────────┼────────────────────┴────────────────────┘
                              ▼
        ┌──────────────────────────────────────────────┐
        │  Bộ Thu thập & Xử lý Tín hiệu Đa kênh        │
        │           NI USB-6003 / NI USB-6009          │
        └─────────────────────┬────────────────────────┘
                              ▼
        ┌──────────────────────────────────────────────┐
        │     Phần mềm Phân tích DasyLab 11.0          │
        │     & Xử lý Dữ liệu Tính toán MATLAB         │
        └──────────────────────────────────────────────┘

Phương pháp tam giác đạc (triangulation) được áp dụng nghiêm ngặt giữa kết quả giải tích lý thuyết, mô phỏng số và dữ liệu thực nghiệm. Độ tin cậy của phép đo được kiểm chuẩn qua 15 chu kỳ xúc thử nghiệm lặp lại, hệ số tin cậy đạt mức $\alpha > 0{,}95$.

Data và phân tích

Hệ phương trình vi phân chuyển động phi tuyến bậc hai được lập trình giải bằng thuật toán Runge-Kutta bậc 4-5 thích ứng (ode45) trong môi trường MATLAB. Đồng thời, mô hình động lực học không gian 3D hoàn chỉnh của máy VMC E500-1 được thiết lập trong phần mềm MSC ADAMS để kiểm chứng chéo (robustness check).

Dữ liệu thực nghiệm thu thập từ mỏ than Hà Lầm được lọc nhiễu số học qua bộ lọc thông thấp Butterworth bậc 4. Kết quả so sánh giữa đường cong tính toán lý thuyết và đường cong thực nghiệm cho thấy sự tương thích xuất sắc: đồ thị lực tại khớp nâng cần $O_2$ và khớp quay gầu $O_3$ có sai số tương đối cực đại không vượt quá $12{,}8%$, hệ số tương quan $R^2 > 0{,}91$, khẳng định độ tin cậy tuyệt đối của mô hình toán học đã xây dựng.

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

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

  1. Phát hiện 1 (Tối ưu hóa quỹ đạo xúc phối hợp): Luận án chứng minh rằng phương pháp xúc phối hợp (đồng thời tịnh tiến thân máy $d_1$, nâng cần $\theta_2$ và quay gầu $\theta_3$) giúp giảm lực cản đào cực đại tác dụng lên gầu từ $28{,}5\text{ kN}$ xuống còn $19{,}2\text{ kN}$ (giảm $32{,}6%$) so với phương pháp xúc phân đoạn đẩy thẳng, đồng thời nâng hệ số điền đầy gầu $K_H$ từ $0{,}78$ lên $0{,}95$.
  2. Phát hiện 2 (Quy luật phân bố tải trọng động tại các khớp): Khảo sát giải tích và thực nghiệm chỉ ra rằng lực tác dụng lên khớp nâng cần $O_2$ đạt giá trị lớn nhất ($F_{O2_max} \approx 65{,}4\text{ kN}$) ở giai đoạn cuối của quá trình nạy tách tảng đá lớn ra khỏi đống đá nổ mìn, chứ không phải ở giai đoạn nâng tải lên cao. Lực tác dụng lên khớp lật gầu $O_5$ trong quá trình xả tải nghiêng một bên đạt đỉnh $F_{O5_max} \approx 22{,}1\text{ kN}$ trong khoảng thời gian $t = 1{,}2 - 1{,}8\text{ s}$.
  3. Phát hiện 3 (Cơ chế trượt xích và giới hạn bám): Khi gầu xúc bị cản trở bởi các cục đá quá cỡ ($D > D_{cp} = 0{,}55\text{ m}$), lực cản ngang $F_{cn}$ tăng đột biến vượt quá lực bám lớn nhất $P_{1max}$ xác định theo điều kiện bám dải xích ($P_k \le G_m \cdot \varphi_x$), dẫn đến hiện tượng trượt quay dải xích $j \ge 15%$, gây tiêu hao năng lượng vô ích và mài mòn kết cấu.
  4. Phát hiện 4 (Tiết kiệm năng lượng chu trình): Bằng việc phối hợp tối ưu hóa vận tốc của các khâu dẫn động thủy lực theo nghiệm của bài toán động học ngược, thời gian của một chu trình công tác $T_Ц$ giảm từ $42\text{ s}$ xuống còn $33\text{ s}$ (rút ngắn $21{,}4%$), làm tăng năng suất kỹ thuật $\Pi_T$ của máy VMC E500-1 lên mức $28 - 35\text{ m}^3/\text{h}$.
Thông số Khảo sát Phương pháp Xúc Phân đoạn Phương pháp Xúc Phối hợp (Luận án) Mức độ Cải thiện
Lực cản đào cực đại ($W_{c_max}$) $28{,}5\text{ kN}$ $19{,}2\text{ kN}$ Giảm $32{,}6%$
Hệ số điền đầy gầu ($K_H$) $0{,}78$ $0{,}95$ Tăng $21{,}8%$
Thời gian chu trình ($T_Ц$) $42\text{ s}$ $33\text{ s}$ Rút ngắn $21{,}4%$
Năng suất kỹ thuật ($\Pi_T$) $22\text{ m}^3/\text{h}$ $32\text{ m}^3/\text{h}$ Tăng $45{,}5%$
Lực động đỉnh tại khớp $O_2$ $82{,}1\text{ kN}$ $65{,}4\text{ kN}$ Giảm $20{,}3%$

Implications đa chiều

  • Về mặt học thuật và lý thuyết: Luận án đã đặt nền móng lý thuyết vững chắc cho chuyên ngành Cơ khí động lực công trình ngầm tại Việt Nam, đóng góp một khung phương pháp luận hoàn chỉnh từ mô hình hạt rời tiếp xúc đến động lực học hệ nhiều vật cho các máy thi công ngầm chuyên dụng.
  • Về mặt phương pháp luận: Quy trình tích hợp đo kiểm thực nghiệm thời gian thực (DasyLab - Tenzo - Cảm biến áp suất/lưu lượng) kết hợp mô phỏng số đồng thời (MATLAB - ADAMS) có thể chuyển giao trực tiếp để nghiên cứu các dòng máy xây dựng khác như máy đào hầm TBM, máy cào vơ, máy khoan hầm thủy lực.
  • Về mặt thực tiễn kỹ thuật: Các đường đặc tính tải trọng động và quy luật chuyển động khâu dẫn là tài liệu kỹ thuật chuẩn để Công ty Cổ phần Chế tạo máy Vinacumin cải tiến kết cấu thép của cần nâng, tối ưu hóa đường kính xi lanh thủy lực và thiết lập thông số cho cụm van phân phối tỷ lệ trên máy VMC E500-1.
  • Về mặt quốc phòng và chính sách: Cung cấp luận cứ khoa học để Bộ Quốc phòng và Binh chủng Công binh phê duyệt phương án cơ giới hóa đồng bộ công tác đào hầm quân sự bí mật khẩu độ nhỏ, giảm bớt quân số thi công thủ công, nâng cao tốc độ đào hầm chiến lược lên gấp 2,5 - 3 lần.

Limitations và Future Research

Mặc dù đạt được những kết quả xuất sắc, luận án thẳng thắn thừa nhận các giới hạn học thuật cần tiếp tục nghiên cứu:

  • Môi trường ĐNM trong mô hình lý thuyết được giả định là môi trường hạt rời đồng nhất với ma sát khô thuần túy, chưa xét đến ảnh hưởng của độ ẩm bão hòa, nước ngầm chảy tràn hoặc bùn nhão tại gương hầm làm thay đổi hệ số liên kết $C$ và góc ma sát trong $\varphi$.
  • Mô hình động lực học hệ nhiều vật xem các khâu (thân máy, cần, giá quay, gầu) là các vật rắn tuyệt đối (rigid bodies), chưa tích hợp tính biến dạng đàn hồi của kết cấu kim loại (Flexible Multibody Dynamics) khi chịu va đập xung kích mạnh từ đá tảng.
  • Khảo sát thực nghiệm mới chỉ thực hiện trên đối tượng duy nhất là máy xúc lật đổ bên VMC E500-1 với một loại gầu xúc có răng tiêu chuẩn dung tích $0{,}5\text{ m}^3$.

Chương trình nghiên cứu mở rộng trong tương lai bao gồm:

  1. Phát triển mô hình mô phỏng liên kết giữa Phương pháp Phần tử rời (Discrete Element Method - DEM) và Động lực học hệ nhiều vật (MBD) để mô phỏng chân thực từng hạt đá vỡ va chạm với gầu xúc.
  2. Tích hợp hệ thống điều khiển thông minh tự động hóa hoàn toàn (Autonomous Loading) dựa trên mạng nơ-ron sâu và thuật toán điều khiển mờ (Fuzzy Logic) để máy tự nhận diện hình dạng đống đá và tự điều chỉnh quỹ đạo xúc tối ưu theo thời gian thực.
  3. Nghiên cứu độ bền mỏi và tích lũy hư hỏng của kết cấu thép khung cần máy VMC E500-1 dưới tác động của phổ tải trọng động ngẫu nhiên thu thập từ thực nghiệm.

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

Công trình nghiên cứu mang lại tác động sâu rộng trên nhiều phương diện:

  • Tác động học thuật: Dự kiến tạo ra nguồn trích dẫn lớn trong các tạp chí chuyên ngành Cơ khí cơ giới hóa xây dựng, Kỹ thuật khai thác mỏ và Công trình quân sự ngầm; đóng vai trò là tài liệu đào tạo sau đại học tại Học viện Kỹ thuật Quân sự, Trường Đại học Mỏ - Địa chất và Đại học Bách khoa Hà Nội.
  • Chuyển đổi công nghiệp khai khoáng: Giúp Tập đoàn Công nghiệp Than - Khoáng sản Việt Nam (Vinacumin) hoàn thiện quy trình công nghệ chế tạo nội địa hóa máy xúc lật đổ bên, giảm giá thành đầu tư thiết bị từ $40 - 50%$ so với việc nhập khẩu máy chuyên dụng từ châu Âu hoặc Nhật Bản.
  • An ninh - Quốc phòng: Ứng dụng trực tiếp vào nhiệm vụ xây dựng các công trình công sự phòng thủ chiến lược của các Quân khu, hầm sở chỉ huy, hầm cất chứa khí tài quân sự có độ ngụy trang cao và cửa hầm nhỏ hẹp.
  • Lợi ích xã hội: Thay thế hoàn toàn phương pháp xúc đá thủ công nặng nhọc và độc hại, giảm thiểu triệt để nguy cơ tai nạn lao động sập đổ đá cục bộ tại gương đào hầm, không gây ô nhiễm khí thải nhờ sử dụng động cơ điện.

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

  • Nghiên cứu sinh và Học viên cao học: Tiếp cận phương pháp luận mẫu mực về mô hình hóa cơ học hệ nhiều vật và kỹ thuật thực nghiệm đo động lực học đa kênh thời gian thực.
  • Giảng viên và Các nhà khoa học cơ khí: Kế thừa hệ thống phương trình vi phân Lagrange loại 2 và thuật toán động học ngược để phát triển cho các thiết bị công trình ngầm phức tạp hơn.
  • Kỹ sư R&D tại các nhà máy chế tạo máy (Vinacumin, Tổng cục Công nghiệp Quốc phòng): Sử dụng trực tiếp bộ thông số hình học, lực tác dụng tại các khớp và quy luật điều khiển xi lanh để cải tiến thiết kế cơ khí và hệ thống thủy lực.
  • Các đơn vị thi công công trình ngầm (Binh chủng Công binh, Tổng công ty Sông Đà, Lũng Lô): Ứng dụng quy trình vận hành xúc phối hợp để nâng cao năng suất thi công, giảm tiêu hao năng lượng và kéo dài tuổi thọ thiết bị.

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 thiết lập mô hình giải tích tích hợp xác định lực cản xúc ĐNM trong không gian hẹp, mở rộng trực tiếp lý thuyết cân bằng giới hạn lăng trụ trượt của V.A. Bauman (1984) và lý thuyết tương tác môi trường rời của A.I. Balốpnhép. Luận án đã bổ sung thành công các thông số hình học biên của gầu ($\alpha_0 = 7^\circ \pm 2^\circ$), góc dốc tự nhiên $\alpha_{VL}$ và hàm phân bố cỡ hạt phi tuyến vào công thức cân bằng lực, khắc phục triệt để hạn chế của mô hình V. Stogov vốn chỉ áp dụng cho phương pháp xúc phân đoạn đơn giản.

2. Điểm mới về phương pháp luận nghiên cứu so với các công trình quốc tế là gì?

So với nghiên cứu của A. Hemani (1994) và Ying Li & Wenyuan Liu (2012), luận án đã:

  • Thiết lập chuỗi động học không gian đầy đủ của máy xúc lật đổ bên bánh xích 4 bậc tự do có tính đến điều kiện biên giới hạn không gian đường hầm ($b \le 6\text{ m}, h \le 3\text{ m}$).
  • Giải quyết trọn vẹn bài toán động học ngược để tìm ra quy luật biến thiên tức thời của 4 cụm xi lanh thủy lực đáp ứng quỹ đạo đỉnh răng gầu tối ưu.
  • Đồng bộ hóa phương pháp mô phỏng số (MATLAB/ADAMS) với hệ thống đo thực nghiệm đa kênh thời gian thực (DasyLab 11, NI-DAQ, Tenzo FCA-5-11, cảm biến áp suất OCM-511) với sai số kiểm chứng dưới $12{,}8%$.
Quỹ đạo Đỉnh Răng Gầu Hợp Lý (Xúc phối hợp)
  │
  ├─► Giải Bài toán Động học Ngược (Ma trận D-H)
  │     │
  │     ├─► Hành trình xi lanh nâng cần x2(t)
  │     ├─► Hành trình xi lanh quay gầu x3(t)
  │     ├─► Hành trình xi lanh lật gầu x4(t)
  │     └─► Chuyển vị tịnh tiến thân máy d1(t)
  │
  └─► Đầu vào cho Mô hình Động lực học MBD (Lagrange loại 2)
        │
        └─► Xuất Tải trọng Động tại Khớp (O2, O3, O5) ──► Kiểm chứng với Đo Tenzo Thực nghiệm

3. Phát hiện bất ngờ nhất trong quá trình khảo sát động lực học là gì?

Phát hiện bất ngờ nhất là giá trị cực đại của lực tác dụng lên khớp bản lề nâng cần $O_2$ ($F_{O2_max} \approx 65{,}4\text{ kN}$) và áp suất dầu trong xi lanh nâng cần $F_{xl2}$ không xuất hiện ở vị trí nâng gầu lên độ cao xả tải lớn nhất mà xuất hiện tại thời điểm kết thúc hành trình ấn gầu và bắt đầu chuyển động quay nạy gầu để tách khối đá nổ mìn. Điều này chứng minh rằng lực cản quán tính và mô-men liên kết hạt rời tại chân đống đá đóng vai trò áp đảo so với trọng trường tĩnh của khối vật liệu.

4. Giao thức lặp lại (Replication Protocol) của nghiên cứu được mô tả ra sao?

Nghiên cứu cung cấp giao thức thực nghiệm hoàn chỉnh bao gồm:

  • Sơ đồ đấu nối chi tiết các cảm biến đo áp suất OCM-511 trên mạch thủy lực nâng và quay gầu.
  • Bản vẽ định vị và quy trình dán cảm biến biến dạng Tenzo FCA-5-11 trên các tai khớp $O_2, O_3, O_5$.
  • Cấu hình khối chức năng trên phần mềm DasyLab 11.0 kết nối qua cổng USB với card NI 6003/NI 6009.
  • Bộ thông số đầu vào chi tiết của máy VMC E500-1 (khối lượng, tọa độ trọng tâm các khâu, đặc tính kéo bánh xích) được công bố minh bạch trong các phụ lục của luận án, cho phép tái lập hoàn toàn mô phỏng trên MATLAB và ADAMS.

5. Lộ trình phát triển 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 giai đoạn chiến lược:

  • Giai đoạn 1 (2-3 năm): Hoàn thiện mô hình động lực học vật rắn biến dạng (Flexible Multibody System) và tối ưu hóa hệ thống thủy lực tỷ lệ điện - thủy lực (Electro-hydraulic Proportional System) cho máy VMC E500-1.
  • Giai đoạn 2 (3-5 năm): Phát triển hệ thống điều khiển bán tự động, tích hợp cảm biến LiDAR quét 3D biên dạng đống đá tại gương hầm để tự động lập quỹ đạo xúc.
  • Giai đoạn 3 (5-10 năm): Chế tạo thế hệ máy xúc lật đổ bên tự hành thông minh hoàn toàn (AI-driven Autonomous Mining Loader) không người lái, ứng dụng pin nhiên liệu hoặc năng lượng sạch trong công trình ngầm sâu.

Kết luận

  1. Luận án đã giải quyết trọn vẹn một bài toán kỹ thuật cấp thiết có ý nghĩa khoa học và thực tiễn cao: tối ưu hóa quá trình bốc xúc vật liệu của máy xúc lật đổ bên bánh xích VMC E500-1 trong thi công đường hầm khẩu độ vừa và nhỏ bằng phương pháp khoan nổ.
  2. Thiết lập thành công mô hình giải tích tương tác gầu - ĐNM có tính đến các yếu tố hình học phi tuyến, làm sáng tỏ bản chất cơ học tiếp xúc và xác định chính xác các thành phần lực cản đào trong điều kiện biên đường hầm hẹp.
  3. Xác định quỹ đạo chuyển động hợp lý của đỉnh răng gầu theo phương pháp xúc phối hợp và giải quyết triệt để bài toán động học ngược, tìm ra quy luật điều khiển tối ưu cho các khâu dẫn động thủy lực, giúp giảm $32{,}6%$ lực cản đào cực đại và rút ngắn $21{,}4%$ thời gian chu trình làm việc.
  4. Xây dựng hoàn chỉnh hệ phương trình vi phân động lực học hệ nhiều vật mô tả chuyển động không gian của cơ hệ, xác định chính xác phổ tải trọng động và quy luật biến thiên lực tại các khớp quay bản lề chịu lực chính ($O_2, O_3, O_5$).
  5. Thực hiện quy trình nghiên cứu thực nghiệm mẫu mực, tin cậy cao trên thiết bị đo hiện đại (DasyLab 11, NI DAQ, cảm biến Tenzo, áp suất, lưu lượng), kiểm chứng sự tương thích xuất sắc giữa mô hình lý thuyết và thực tế với sai số dưới $12{,}8%$.
  6. Công trình mở ra hướng ứng dụng đột phá trong việc nâng cấp, hoàn thiện kết cấu dòng máy VMC E500-1 của Tập đoàn Vinacumin, đồng thời cung cấp luận cứ khoa học vững chắc phục vụ chiến lược hiện đại hóa trang bị công binh trong xây dựng công trình ngầm quốc phòng tại Việt Nam.