Tổng quan về luận án

Thi công các công trình ngầm và đường hầm quân sự (ĐHQS) khẩu độ vừa và nhỏ ($S_{tđ} < 16\text{ m}^2$, bề rộng $B < 4\text{ m}$) trong điều kiện địa hình chia cắt, địa chất phức tạp đặt ra yêu cầu cấp thiết về cơ giới hóa đồng bộ công tác khoan nổ mìn. Luận án tiến sĩ kỹ thuật chuyên ngành Kỹ thuật Cơ khí Động lực (Mã số: 9 52 01 16) của nghiên cứu sinh Bùi Văn Trầm, dưới sự hướng dẫn khoa học của TS. Nguyễn Lâm Khánh và NGƯT. Chu Văn Đạt tại Học viện Kỹ thuật Quân sự (2019), mang tiêu đề "Nghiên cứu động lực học thiết bị khoan xoay đập lắp trên máy đào phục vụ thi công hầm khẩu độ vừa và nhỏ", đã tiên phong giải quyết bài toán cốt lõi: tối ưu hóa không gian thao tác và kiểm soát rung động động lực học của hệ giá đỡ máy khoan nội địa MK-Z49 lắp trên máy cơ sở Yanmar.

Trong thực tiễn thi công công trình ngầm, việc nhập khẩu các tổ hợp khoan hầm tự hành cỡ lớn từ các nước công nghiệp phát triển gặp rào cản nghiêm trọng về kích thước cồng kềnh, chi phí đầu tư đắt đỏ và sự phụ thuộc hoàn toàn vào chuỗi cung ứng phụ tùng thay thế. Ngược lại, các thiết bị giá khoan chế tạo trong nước trước đây chủ yếu dựa trên phương pháp sao chép cơ học thuần túy hoặc kinh nghiệm gia công rời rạc, dẫn đến hai khoảng trống nghiên cứu (research gaps) mang tính then chốt:

  1. Thiếu mô hình động học không gian tổng quát dẫn đến hệ giá đỡ bị "điểm mù", không thể bao phủ toàn bộ các vị trí lỗ khoan theo hộ chiếu gương hầm kích thước $2,8\text{ m} \times 1,8\text{ m}$ (đặc biệt là các lỗ viền vách có góc nghiêng $\alpha_h = 5^\circ - 10^\circ$ và lỗ đột khẩu $\alpha_h = 12^\circ - 15^\circ$).
  2. Thiếu mô hình động lực học hệ nhiều vật liên kết mạch vòng có xét đến tính nén đàn hồi phi tuyến của dầu thủy lực và tương tác đá - mũi khoan, khiến giá dẫn hướng bị rung lắc dữ dội, làm sai lệch đường tâm lỗ khoan vượt ngưỡng tiêu chuẩn kỹ thuật Việt Nam $\Delta x \le 1,5 L_k$ (cm) và gây hiện tượng bó kẹt cần khoan nghiêm trọng.

Để giải quyết triệt để vấn đề này, luận án xác định hai câu hỏi nghiên cứu và giả thuyết khoa học tường minh:

  • Câu hỏi nghiên cứu 1 (RQ1): Cần điều chỉnh cấu trúc động học không gian của hệ giá đỡ như thế nào để khâu thao tác cuối đạt góc xoay linh hoạt, phủ kín $100%$ diện tích gương hầm $2,8\text{ m} \times 1,8\text{ m}$ mà không gây va chạm hình học với vách hầm?
  • Giả thuyết 1 (H1): Việc bổ sung khớp quay $q_4$ tại chụp đầu cần với dải biến thiên góc quay $q_4 \in [90^\circ, 270^\circ]$ sẽ mở rộng không gian thao tác trong hệ tọa độ Descartes $OXYZ$, cho phép điều chỉnh độc lập hướng véc-tơ trục khoan mà không phụ thuộc vào tư thế thân máy cơ sở.
  • Câu hỏi nghiên cứu 2 (RQ2): Mức độ rung lắc của điểm tỳ giá dẫn hướng chịu tác động ra sao dưới các chế độ tải trọng va đập tần số cao của búa thủy lực kết hợp với áp suất và dung tích dầu trong các khoang xi lanh?
  • Giả thuyết 2 (H2): Mô hình hóa xi lanh thủy lực dẫn động dưới dạng phần tử đàn nhớt phi tuyến có độ cứng tương đương $c_{d(tđ)} = c_{d(A)} + c_{d(B)}$ tích hợp trong hệ phương trình vi phân - đại số (DAE) Lagrange dạng nhân tử sẽ phản ánh chính xác đáp ứng dịch chuyển rung chấn của điểm $S$ trên giá dẫn hướng, cho phép xác định kích thước đường kính xi lanh tối ưu nhằm triệt tiêu rung chấn.

Khung lý thuyết của luận án tích hợp ma trận truyền vị trí Denavit - Hartenberg (D-H), cơ học hệ nhiều vật mạch vòng kín (closed-loop multibody dynamics) và lý thuyết tương tác cơ học đất đá Hsieh - Reissner. Đóng góp đột phá được lượng hóa thông qua việc tái cấu trúc thành công máy khoan MK-Z49 (trọng lượng $3160\text{ kg}$, búa YYTZ28 công suất $22\text{ kW}$, tần số đập $60\text{ Hz}$, năng lượng va đập $65\text{ J}$), giảm biên độ rung lắc điểm dẫn hướng về dưới $0,5\text{ mm}$, khống chế sai lệch lỗ khoan $\Delta x \le 1,5\text{ cm}$ trên toàn bộ chiều sâu khoan $L_k = 1,4\text{ m}$, nâng cao năng suất đào hầm lên $25 - 30%$ và đảm bảo an toàn tuyệt đối cho kết cấu chống đỡ.

Literature Review và Positioning

Tổng quan y văn quốc tế và trong nước về lĩnh vực thiết bị thi công ngầm ghi nhận ba dòng nghiên cứu chính với những luận điểm và sự đối lập sâu sắc:

                      ┌──────────────────────────────────────────────────────────┐
                      │             DÒNG CHẢY NGHIÊN CỨU Y VĂN                   │
                      └──────────────────────────┬───────────────────────────────┘
                                                 │
         ┌───────────────────────────────────────┼───────────────────────────────────────┐
         ▼                                       ▼                                       ▼
┌────────────────────────────────┐ ┌───────────────────────────┐ ┌────────────────────────────────┐
│   Tương tác Mũi khoan - Đá     │ │  Động lực học Xi lanh &   │ │   Cơ học Hệ nhiều vật &        │
│                                │ │     Rung động Thủy lực    │ │     Không gian Thao tác        │
├────────────────────────────────┤ ├───────────────────────────┤ ├────────────────────────────────┤
│ • Hustrulid & Fairhurst (1971) │ │ • P. Korchagin (2010):    │ │ • Denavit & Hartenberg (1955)  │
│   Sóng ứng suất 1D búa đập     │ │   Lò xo tuyến tính 10-DOF │ │ • Nguyen Van Khang (2007-2017) │
│ • Lundberg (1973): Tải xoay rời│ │ • Tran Xuan Tuy (2004):   │ │   DAE sang ODE mạch vòng       │
│ • Chiang & Elias (2000):       │ │   Hệ thủy lực 1 bậc tự do │ │ • Chu Anh My (2013-2018):      │
│   Mô men cắt thuần túy         │ │ • Luận án:                │ │   Robot lai nối tiếp - song    │
│ • Luận án: Tích hợp hàm hóa    │ │   Độ cứng phi tuyến 3D    │ │   song 6 bậc tự do             │
│   tải trọng thực nghiệm 60 Hz  │ │   $c_{d(tđ)} = f(x, p, V)$│ │ • Luận án: Mở rộng khớp $q_4$  │
└────────────────────────────────┘ └───────────────────────────┘ └────────────────────────────────┘
  1. Dòng nghiên cứu tương tác mũi khoan và đất đá: Hustrulid và Fairhurst (1971), Lundberg (1973) cùng các nghiên cứu mở rộng của Chiang và Elias (2000) đã phân tích sâu cơ chế truyền sóng ứng suất trong choòng khoan đập đỉnh. Tuy nhiên, các tác giả chủ yếu khảo sát tải trọng va đập thuần túy dọc trục hoặc tải trọng cắt xoay tĩnh riêng biệt. Tranh luận học thuật lớn tồn tại giữa quan điểm của Whittaker (1992) – cho rằng lực cắt tế vi do ma sát đóng vai trò quyết định độ mài mòn mũi khoan – và quan điểm của Lundberg (1973) – khẳng định ứng suất nén xung kích từ sóng phản xạ mới là nhân tố chính gây phá hủy thể tích đá. Luận án đã định vị khoảng trống: chưa có nghiên cứu nào kết hợp đồng thời mô men cản cắt $M_C(t)$ và lực đập xung kích tần số cao $F_G(t)$ ở tần số $60\text{ Hz}$ tác động trực tiếp lên hệ giá đỡ không gian trong điều kiện địa chất đa dạng ($f_{kp} = 3 - 7$).
  2. Dòng nghiên cứu động lực học máy công trình và xi lanh thủy lực: Nghiên cứu kinh điển của P. Korchagin (2010) khảo sát mô hình máy đào 10 bậc tự do (10-DOF) nhưng đã đơn giản hóa xi lanh thủy lực thành các phần tử đàn nhớt tuyến tính cố định; trong khi đó, công trình của Trần Xuân Tùy (2004) chỉ xem xét tính nén của dầu trong chuyển động tịnh tiến một chiều tuyến tính. Điểm nghẽn là việc quy dẫn tuyến tính cố định bỏ qua sự biến thiên thể tích khoang dầu $V_A(t), V_B(t)$ theo hành trình công tác của piston, dẫn đến sai số động lực học đáng kể khi cơ hệ chịu tải chấn động mạnh.
  3. Dòng nghiên cứu cơ cấu mạch vòng và robot công nghiệp: Các công trình của Nguyễn Văn Khang (2007, 2011, 2017) và Chu Anh Mỳ (2013, 2018) đã phát triển mạnh phương pháp giải hệ nhiều vật mạch vòng kín và robot song song bằng phương trình vi phân - đại số (DAE). Tuy nhiên, các ứng dụng chủ yếu dừng lại ở tay máy công nghiệp hoặc thiết bị hàn phẳng.

So sánh với hai nghiên cứu quốc tế điển hình:

  • Nghiên cứu của P. Korchagin (2010) về động lực học máy công trình: Chỉ áp dụng mô hình phẳng, bỏ qua tương tác không gian ba chiều và không gian hẹp của công trình ngầm. Luận án của Bùi Văn Trầm vượt trội nhờ thiết lập mô hình không gian 6 khâu hoàn chỉnh với 5 xi lanh thủy lực điều khiển phi tuyến.
  • Nghiên cứu của Chiang & Elias (2000) về cơ chế phá hủy đất đá: Chỉ khảo sát thiết bị khoan cố định trên bệ thử, thiếu tính tương tác ngược giữa phản lực đất đá và độ biến dạng đàn hồi của hệ thống đỡ di động.

Luận án đã lấp đầy các khoảng trống trên bằng cách thiết lập một khung phân tích thống nhất từ động học D-H, động lực học DAE hệ nhiều vật đến tương tác đàn nhớt đá - mũi khoan Hsieh - Reissner.

Đó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 động lực học hệ nhiều vật (Multibody System Dynamics) của Schiehlen (1997) và Shabana (2013), đồng thời kế thừa trường phái cơ học giải tích Lagrange tại Việt Nam (Nguyễn Văn Khang, Đỗ Sanh), mang lại những bước tiến lý thuyết cụ thể:

  • Phát triển mô hình động học không gian 6 bậc tự do: Thiết lập hệ thức ma trận biến đổi thuần nhất Denavit - Hartenberg (D-H) với việc bổ sung biến khớp quay chủ động $q_4$ tại khâu chụp đầu cần:

    $$\mathbf{T}6^0 = \prod{i=1}^6 \mathbf{H}_{i}^{i-1}(q_i)$$

    Mô hình này giải quyết triệt để bài toán động học ngược (Inverse Kinematics), xác định chuẩn xác bộ thông số cấu hình khớp $\mathbf{q} = [q_1, q_2, q_3, q_4, q_5, q_6]^T$ tương ứng với mọi tọa độ điểm khoan và véc-tơ chỉ phương mũi khoan trên gương hầm $2,8\text{ m} \times 1,8\text{ m}$.

  • Thiết lập hệ phương trình vi phân - đại số (DAE) chỉ số 3: Mô tả chính xác động lực học hệ giá đỡ dạng mạch vòng kín chịu liên kết hình học phi tuyến bằng phương trình Lagrange dạng nhân tử:

    $$\mathbf{M}(\mathbf{q})\ddot{\mathbf{q}} + \boldsymbol{\Phi}_{\mathbf{q}}^T \boldsymbol{\lambda} = \mathbf{Q}(\mathbf{q}, \dot{\mathbf{q}}, t)$$

    $$\boldsymbol{\Phi}(\mathbf{q}, t) = \mathbf{0}$$

    Trong đó, $\mathbf{M}(\mathbf{q})$ là ma trận khối lượng suy rộng kích thước $6 \times 6$, $\boldsymbol{\Phi}_{\mathbf{q}}$ là ma trận Jacobi của các phương trình liên kết, $\boldsymbol{\lambda}$ là véc-tơ nhân tử Lagrange biểu diễn nội lực liên kết tại các khớp trượt và khớp quay, $\mathbf{Q}$ là véc-tơ lực suy rộng bao gồm ngoại lực khoan và phản lực thủy lực.

+---------------------------------------------------------------------------------------+
|                 KHUNG ĐỘNG LỰC HỌC HỆ NHIỀU VẬT KHÔNG GIAN CỦA MK-Z49                 |
+---------------------------------------------------------------------------------------+
|                                                                                       |
|   [ Ngoại lực đập - xoay ]       [ Đàn nhớt chất lỏng ]        [ Phản lực đá Hsieh ]  |
|      F_G(t) & M_G(t) 60Hz       c_d(tđ) = K'(S_A^2/V_A + ...)    c_đ & k_đ (C30-C70)  |
|             │                               │                           │             |
|             └───────────────────────┬───────┴───────────────────────────┘             |
|                                     ▼                                                 |
|               ┌───────────────────────────────────────────────┐                       |
|               │  HỆ PHƯƠNG TRÌNH VI PHÂN - ĐẠI SỐ (DAE CHỈ SỐ 3) │                  |
|               │                                               │                       |
|               │   M(q) * q'' + Φ_q^T * λ = Q(q, q', t)        │                       |
|               │   Φ(q, t) = 0                                 │                       |
|               └──────────────────────┬────────────────────────┘                       |
|                                      │                                                |
|                                      ▼ (Thuật toán biến đổi Khang - Sanh)             |
|               ┌───────────────────────────────────────────────┐                       |
|               │       HỆ PHƯƠNG TRÌNH ODE CỰC TIỂU TRÊN C30   │                       |
|               │   Mô phỏng Runge-Kutta / MAPLE - MATLAB       │                       |
|               └──────────────────────┬────────────────────────┘                       |
|                                      │                                                |
|                                      ▼                                                |
|               [ Phân tích sai lệch quỹ đạo rung chấn điểm S ]                         |
|               [ Khống chế dung sai: Δx <= 1.5 cm theo TCVN  ]                         |
|                                                                                       |
+---------------------------------------------------------------------------------------+

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

Sự tích hợp liên ngành giữa ba nền tảng lý thuyết:

  1. Lý thuyết thủy khí động lực học: Mô hình hóa tính nén đàn nhớt của dầu thủy lực công tác thông qua mô đun đàn hồi thể tích $K'$ và dung tích thay đổi theo thời gian:

    $$c_{d(tđ)} = c_{d(A)} + c_{d(B)} = K' \left( \frac{S_A^2}{V_{A0} + S_A x} + \frac{S_B^2}{V_{B0} - S_B x} \right)$$

    Công thức này phản ánh chính xác trạng thái độ cứng tương đương đạt giá trị cực tiểu tại vị trí cân bằng giữa hành trình và tăng phi tuyến về hai phía đầu mút xi lanh.

  2. Lý thuyết tương tác cơ học môi trường tiếp xúc Hsieh - Reissner: Tính toán độ cứng tương đương $c_đ$ và độ cản nhớt $k_đ$ của vách gương hầm tỳ vào mũi định vị giá khoan:

    $$c_đ = G \cdot r \cdot \frac{f_1}{f_1^2 + f_2^2}, \quad k_đ = \frac{r^2 \sqrt{G \rho}}{a_0} \cdot \frac{-f_2}{f_1^2 + f_2^2}$$

    Với $G$ là mô đun đàn hồi trượt, $\rho$ là khối lượng riêng của đá, $r$ là bán kính tiếp xúc thực tế và $f_1, f_2$ là các hàm dịch chuyển Reissner không thứ nguyên.

  3. Kỹ thuật biến đổi DAE sang ODE: Ứng dụng phương pháp tọa độ suy rộng cực tiểu nhằm loại bỏ các nhân tử Lagrange $\boldsymbol{\lambda}$, biến đổi hệ phương trình DAE chỉ số cao phức tạp thành hệ phương trình vi phân thường (ODE) cấp 2 hoàn toàn khả tích, đảm bảo tính ổn định số và triệt tiêu sai số trôi dạt liên kết (constraint drift).

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 học thực chứng (positivism) với thiết kế nghiên cứu kết hợp chặt chẽ giữa mô hình hóa giải tích toán học, mô phỏng số trên máy tính và kiểm chứng thực nghiệm hiện trường đa cấp độ (multi-level validation).

                      ┌────────────────────────────────────────────────────────┐
                      │              QUY TRÌNH THỰC NGHIỆM ĐA CẤP ĐỘ           │
                      └───────────────────────────┬────────────────────────────┘
                                                  │
         ┌────────────────────────────────────────┴────────────────────────────┐
         ▼                                                                     ▼
┌───────────────────────────────────────────────┐     ┌────────────────────────────────────────────────┐
│   CẤP ĐỘ 1: PHÒNG THÍ NGHIỆM LAS-XD72         │     │    CẤP ĐỘ 2: HIỆN TRƯỜNG NHÀ MÁY Z49           │
├───────────────────────────────────────────────┤     ├────────────────────────────────────────────────┤
│ • Mẫu bê tông chuẩn C30, C40, C50, C60, C70   │     │ • Máy khoan nguyên mẫu MK-Z49 trên Yanmar      │
│ • Kiểm chuẩn cường độ nén và độ kiên cố f_kp  │     │ • Đo tải áp suất p_1(t), p_2(t) tại xi lanh    │
│ • Cảm biến áp lực strain-gauge Kyowa          │     │ • Cảm biến gia tốc rung chấn 3D bệ dẫn hướng   │
│ • Đo biến dạng tấm thép bệ gá động cơ         │     │ • Thử nghiệm khoan gương nổ mìn thực địa       │
└───────────────────────┬───────────────────────┘     └────────────────────────┬───────────────────────┘
                        │                                                      │
                        └──────────────────────────┬───────────────────────────┘
                                                   ▼
                      ┌────────────────────────────────────────────────────────┐
                      │    ĐỐI SOÁT DỮ LIỆU SỐ (MAPLE / MATLAB Runge-Kutta)    │
                      │    Độ lệch pha cực đại < 8.5%, R^2 = 0.934 - 0.958     │
                      └────────────────────────────────────────────────────────┘
  • Mức độ 1 (Cơ cấu chấp hành): Khảo sát độc lập từng xi lanh dẫn động thủy lực nâng hạ cần ($x_1$) và nâng hạ giá trượt ($x_2$).
  • Mức độ 2 (Hệ thống không gian): Khảo sát liên kết toàn bộ 6 khâu của hệ giá đỡ máy khoan MK-Z49 trong không gian ba chiều.
  • Mức độ 3 (Tương tác hiện trường): Khảo sát động lực học đồng thời giữa máy cơ sở Yanmar ($2600\text{ kg}$), hệ giá đỡ ($560\text{ kg}$), cụm đầu khoan YYTZ28 ($30\text{ kg}$) và khối đá hầm có độ kiên cố Protodyakonov $f_{kp} = 3 - 10$.

Quy trình nghiên cứu rigorous

Quy trình thu thập và xử lý dữ liệu được thiết kế nghiêm ngặt:

  • Chiến lược lấy mẫu và đối tượng thí nghiệm: Sử dụng các khối mẫu bê tông chuẩn hóa có cấp độ bền C30, C40, C50, C60, C70 tương ứng với hệ số độ kiên cố đất đá $f_{kp} = 3, 4, 5, 6, 7$, mô phỏng chính xác các loại đá vôi mềm, sa thạch, cát kết và granit thường gặp trong các công trình công binh của Lữ đoàn 131 Hải quân và Quân khu 3.
  • Thiết bị đo lường và giao thức thu thập:
    • Cảm biến biến dạng điện trở (strain-gauge) dán trực tiếp lên tấm thép gá cụm động cơ nhằm đo tải trọng nén và mô men uốn $M_x$. Hệ số đầu đo $K$, cầu đo điện áp nuôi $V_0$ truyền dữ liệu qua bộ biến đổi tín hiệu tốc độ cao.
    • Cảm biến áp suất thủy lực lắp tại khoang nâng hạ cần và khoang giá trượt, đo dải áp suất công tác $0 - 25\text{ MPa}$.
    • Cảm biến đo chuyển vị và gia tốc rung chấn đa trục gắn tại điểm $S$ trên đầu giá dẫn hướng khoan.
  • Xử lý và làm trơn số liệu: Chuỗi dữ liệu thời gian thực được lọc nhiễu dải thông (band-pass filter), sau đó khai triển chuỗi Fourier và hồi quy phi tuyến theo phương pháp bình phương tối thiểu (least-squares method), xác định chính xác các hàm lực $F_G(t)$ và mô men $M_G(t)$ có tần số cơ bản $60\text{ Hz}$ trùng khớp với tần suất va đập của búa khoan thủy lực YYTZ28.

Data và phân tích

Toàn bộ thuật toán mô phỏng được lập trình trên phần mềm chuyên dụng MAPLE và MATLAB, sử dụng thuật toán tích phân số Runge-Kutta bậc 4-5 thích nghi (ode45/ode15s) chuyên trị hệ phương trình vi phân cứng.

Bảng tổng hợp tham số cơ học và thông số động lực học đầu vào:

Tham số đại lượng Ký hiệu Giá trị đơn vị Ý nghĩa cơ học
Trọng lượng máy cơ sở Yanmar $M_0$ $2600\text{ kg}$ Khối lượng nền cơ sở cố định
Áp suất thủy lực cực đại $p_{max}$ $17,5\text{ MPa}$ Giới hạn van an toàn hệ thống
Tần số đập của búa YYTZ28 $f_d$ $60\text{ Hz}$ Tần suất xung kích phá hủy đá
Năng lượng va đập búa $E_d$ $65\text{ J}$ Năng lượng một nhát đập piston
Tốc độ quay cụm đầu búa $n_k$ $230\text{ vg/ph}$ Tốc độ cắt xoay nghiền phoi
Độ cứng đá mẫu C30 ($f_{kp}=3$) $c_đ$ $3,89 \cdot 10^7\text{ N/m}$ Khả năng đàn hồi chống lún đá
Hệ số cản nhớt đá C30 $k_đ$ $1,12 \cdot 10^5\text{ N.s/m}$ Hệ số suy giảm năng lượng sóng
Độ cứng đá mẫu C70 ($f_{kp}=7$) $c_đ$ $8,50 \cdot 10^7\text{ N/m}$ Độ cứng mẫu đá kiên cố cao
Hệ số cản nhớt đá C70 $k_đ$ $2,45 \cdot 10^5\text{ N.s/m}$ Độ tiêu tán dao động đá cứng

Độ tin cậy của mô hình số được khẳng định qua hệ số xác định $R^2 \ge 0,934$ trên tất cả các kênh đo áp suất và chuyển vị, với khoảng tin cậy $95%$ ($p < 0,001$).

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

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

  1. Khám phá và triệt tiêu "vùng mù động học": Phân tích không gian thao tác nguyên bản của máy MK-Z49 chỉ ra rằng cơ cấu cũ chỉ có 5 bậc tự do, khiến đầu khoan không thể tiếp cận các lỗ viền vách đáy và vòm hầm $2,8\text{ m} \times 1,8\text{ m}$ dưới góc nghiêng tiêu chuẩn $5^\circ - 10^\circ$. Sau khi bổ sung khớp quay $q_4$ cho phép quay quanh trục cần từ $90^\circ$ đến $270^\circ$, không gian thao tác đã bao phủ $100%$ diện tích gương hầm, giải quyết hoàn toàn bài toán công nghệ thi công lỗ khoan viền và đột khẩu.
  2. Quy luật phi tuyến về đáp ứng rung chấn của điểm tỳ dẫn hướng $S$: Quỹ đạo mặt phẳng pha $YOZ$ của điểm $S$ cho thấy rung động dao động có tính chu kỳ biên độ nhỏ nhưng phụ thuộc phi tuyến vào độ kiên cố của đất đá. Khi độ cứng đất đá tăng từ mẫu C30 ($f_{kp} = 3$) lên C70 ($f_{kp} = 7$), biên độ dịch chuyển chuyển vị cưỡng bức của điểm $S$ giảm từ $0,68\text{ mm}$ xuống $0,24\text{ mm}$ do năng lượng phản xạ từ vách đá cứng hạn chế biến dạng trượt của mũi tỳ.
  3. Mối quan hệ định lượng giữa đường kính xi lanh và độ ổn định: Mô phỏng số trên dải đường kính xi lanh từ $60\text{ mm}$ đến $100\text{ mm}$ chỉ ra rằng đường kính xi lanh nâng hạ cần ($d_{XL1}$) và xi lanh nâng hạ giá trượt ($d_{XL2}$) có ảnh hưởng chi phối đến độ cứng thủy lực tương đương $c_{d(tđ)}$. Khi tăng đường kính piston xi lanh thêm $20%$, độ cứng vững toàn hệ thống tăng $44%$, làm giảm biên độ rung lắc điểm $S$ tới $36,2%$.
  4. Kiểm chứng độ chính xác viền lỗ khoan: Dữ liệu thực nghiệm khẳng định sai lệch viền lỗ khoan cực đại đo được trên mẫu C30 là $\Delta x_{max} = 1,28\text{ cm}$ và trên mẫu C70 là $\Delta x_{max} = 0,86\text{ cm}$ (đo tại chiều sâu khoan $L_k = 1,2\text{ m}$), hoàn toàn thỏa mãn tiêu chuẩn kỹ thuật Việt Nam $\Delta x \le 1,5 L_k = 1,5 \times 1,2 = 1,8\text{ cm}$.
  Sai lệch lỗ khoan Δx (cm)
    2.0 ┬──────────────────────────────────────────────────────────────────
        │                                        [Ngưỡng TCVN: Δx = 1.8 cm]
    1.5 ┼ - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
        │      ● C30 (Δx = 1.28 cm)
    1.0 ┼                             ● C50 (Δx = 1.02 cm)
        │                                                    ● C70 (Δx = 0.86 cm)
    0.5 ┼
        │
    0.0 ┴──────────┬──────────────────────────┬──────────────────────────┬─
                  C30                        C50                        C70
                                     Độ cứng mẫu bê tông / Đá
  1. Sự tương đồng giữa lý thuyết và thực nghiệm hiện trường: Sai số tương đối giữa áp suất tính toán lý thuyết và áp suất đo thực nghiệm tại xi lanh nâng hạ cần ($x_1$) không vượt quá $8,4%$, và tại xi lanh giá trượt ($x_2$) không vượt quá $7,6%$. Độ lệch pha dao động thời gian cực đại $< 0,005\text{ s}$, chứng minh tính đúng đắn tuyệt đối của mô hình toán DAE tích hợp độ cứng thủy lực phi tuyến.

Implications đa chiều

  • Tiến bộ học thuật: Cung cấp phương pháp luận tổng quát để giải các bài toán động lực học không gian hệ nhiều vật mạch vòng kín có chứa các khâu đàn nhớt phi tuyến chịu kích động ngẫu nhiên tần số cao ($60\text{ Hz}$), mở rộng ứng dụng cho các dòng máy xây dựng và tay máy thủy lực phức tạp.
  • Đổi mới phương pháp luận: Xây dựng thành công thuật toán kết nối hai chiều giữa bài toán động học ngược (D-H) và bài toán động lực học thuận (DAE), cho phép tính toán tức thời ứng suất và phản lực khớp tại từng vị trí lỗ khoan bất kỳ trên gương hầm.
  • Ứng dụng công nghiệp quốc phòng và dân dụng: Cung cấp tài liệu thiết kế và hồ sơ công nghệ hoàn chỉnh cho Nhà máy Z49 - Binh chủng Công binh để nâng cấp dây chuyền sản xuất máy khoan MK-Z49 loạt số 0, thay thế hoàn toàn máy khoan nhập khẩu trong thi công hầm thủy điện nhỏ và hầm phòng thủ chiến lược.
  • Chính sách và tiêu chuẩn hóa: Đề xuất bổ sung tiêu chuẩn nghiệm thu độ ổn định động lực học và rung chấn cho thiết bị khoan ngầm chế tạo trong nước vào hệ thống tiêu chuẩn kỹ thuật của Bộ Quốc phòng và Bộ Giao thông Vận tải.

Limitations và Future Research

Mặc dù đạt được những kết quả đột phá, công trình vẫn tồn tại một số giới hạn nghiên cứu mang tính biên định:

  1. Giả thiết về độ cứng khâu: Luận án coi các khâu cơ học và khớp bản lề là cứng tuyệt đối (rigid bodies), rung động toàn hệ quy tụ về tính nén đàn nhớt của cột dầu thủy lực. Trong thực tế làm việc lâu dài, hiện tượng biến dạng đàn hồi của kết cấu cần thép hộp và độ rơ lỏng mài mòn tại các chốt khớp bản lề sau hàng nghìn chu kỳ va đập có thể gây ra các thành phần dao động ký sinh bậc cao.
  2. Hiện tượng nhiệt - thủy lực: Bỏ qua sự suy giảm độ nhớt và mô đun đàn hồi thể tích $K'$ của dầu thủy lực khi nhiệt độ dầu tăng cao ($> 65^\circ\text{C}$) trong ca làm việc liên tục kéo dài.
  3. Môi trường địa chất ngẫu nhiên: Mô hình tương tác đá dựa trên các mẫu bê tông đồng nhất (C30 - C70), chưa bao quát hết các đới đứt gãy kiến tạo ngẫu nhiên, hang karst hoặc đất đá ngậm nước áp lực cao.

Chương trình nghiên cứu tương lai (5 - 10 năm):

  • Phát triển mô hình động lực học hệ nhiều vật đàn hồi (Flexible Multibody System Dynamics) kết hợp phương pháp phần tử hữu hạn (FEM) để khảo sát hiện tượng mỏi kết cấu hệ giá đỡ dưới tải trọng va đập $60\text{ Hz}$.
  • Tích hợp hệ thống điều khiển thích nghi thông minh (Adaptive Electro-hydraulic Control) sử dụng mạng nơ-ron để tự động điều chỉnh áp suất và lưu lượng cấp cho xi lanh tùy theo độ cứng thực tế của từng lớp đá khi mũi khoan tiến sâu.
  • Nghiên cứu cơ chế tự động hóa định vị đường tâm lỗ khoan bằng cảm biến laser và thị giác máy tính (Computer Vision) trong điều kiện bụi hầm nồng độ cao.

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

  • Tác động học thuật: Dự báo công trình sẽ thu hút từ 50 - 80 trích dẫn học thuật chuyên ngành trong các tạp chí cơ học, máy xây dựng và kỹ thuật công trình ngầm; trở thành tài liệu tham khảo mẫu mực tại Học viện Kỹ thuật Quân sự, Trường Đại học Bách khoa Hà Nội và Trường Đại học Giao thông Vận tải.
  • Chuyển dịch công nghiệp cơ khí chế tạo: Giúp các nhà máy công binh nội địa hóa $100%$ quy trình thiết kế và chế tạo hệ giá đỡ máy khoan tự hành, giảm chi phí sản xuất thiết bị xuống chỉ bằng $35 - 40%$ so với việc nhập khẩu máy khoan đồng bộ từ nước ngoài (tiết kiệm hàng chục tỷ đồng cho mỗi dự án hầm).
  • Ý nghĩa quốc phòng - an ninh: Nâng cao vượt bậc tốc độ thi công công trình chiến đấu, công trình hầm ngầm cất giấu vũ khí trang bị kỹ thuật tại các địa bàn chiến lược hiểm trở, bảo đảm tính bí mật, cơ động và sẵn sàng chiến đấu cao của Quân đội Nhân dân Việt Nam.

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

  • Nghiên cứu sinh và học viên sau đại học: Tiếp cận phương pháp luận chuẩn xác về giải toán DAE hệ nhiều vật mạch vòng và kỹ thuật đo lường động lực học thủy lực thực nghiệm.
  • Kỹ sư thiết kế R&D tại các nhà máy cơ khí (Z49, Z751, Viện Cơ ký chính xác): Ứng dụng trực tiếp bộ thông số kết cấu tối ưu về đường kính xi lanh, hành trình khớp và cấu hình giá trượt để chế tạo các dòng máy khoan hầm thế hệ mới.
  • Chỉ huy các đơn vị công binh công trình và nhà thầu xây dựng ngầm: Sở hữu cẩm nang công nghệ vận hành chuẩn xác, điều chỉnh áp suất khoan tối ưu theo từng cấp đất đá, kéo dài tuổi thọ choòng khoan và bảo vệ an toàn công trình.
  • Cơ quan quản lý kỹ thuật quân sự: Có cơ sở khoa học tin cậy để thẩm định, phê duyệt và tiêu chuẩn hóa các dự án cải tiến vũ khí trang bị kỹ thuật công binh.

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 tích hợp thành công mô hình độ cứng thủy lực phi tuyến phụ thuộc hành trình $c_{d(tđ)} = K'(S_A^2/V_A + S_B^2/V_B)$ vào hệ phương trình vi phân - đại số (DAE) Lagrange dạng nhân tử cho cơ hệ không gian 6 bậc tự do. Công trình này mở rộng trực tiếp lý thuyết động lực học hệ nhiều vật của Schiehlen và Shabana sang lĩnh vực máy thi công ngầm chịu tải chấn động va đập chu kỳ cao ($60\text{ Hz}$).

2. Đổi mới phương pháp luận so với hai nghiên cứu quốc tế tiêu biểu?

So với P. Korchagin (2010) (chỉ nghiên cứu mô hình phẳng với xi lanh tuyến tính hóa) và Chiang & Elias (2000) (chỉ đo lực cắt phá hủy đá tĩnh), luận án đã:

  1. Xây dựng mô hình động học không gian 3D đầy đủ với 6 biến khớp D-H, bổ sung khớp quay $q_4 \in [90^\circ, 270^\circ]$.
  2. Thiết lập hàm hóa ngoại lực ngẫu nhiên thực tế qua chuỗi Fourier thực nghiệm kết hợp lý thuyết tương tác Hsieh - Reissner.

3. Phát hiện bất ngờ nhất có số liệu thực nghiệm chứng minh là gì?

Hiện tượng biên độ rung lắc của điểm tỳ giá dẫn hướng $S$ giảm mạnh khi độ cứng của đối tượng khoan tăng lên: biên độ dao động tại mẫu C30 ($f_{kp} = 3$) đạt $0,68\text{ mm}$, trong khi tại mẫu đá cứng C70 ($f_{kp} = 7$) chỉ còn $0,24\text{ mm}$. Kết quả này chứng minh rằng phản lực đàn hồi tiếp xúc $c_đ$ từ đá cứng đóng vai trò như một gối tựa thủy tĩnh hỗ trợ dập tắt dao động của giá dẫn hướng, trái ngược với giả định trực giác ban đầu rằng đá cứng sẽ làm tăng rung lắc bệ máy.

4. Giao thức tái lập nghiên cứu (Replication Protocol) có được cung cấp đầy đủ?

Luận án cung cấp tường minh $100%$ các thông số hình học ma trận D-H, bảng khối lượng tâm quán tính 6 khâu, thông số kích thước 5 xi lanh thủy lực, các câu lệnh thiết lập DAE trong phần mềm MAPLE/MATLAB và sơ đồ bố trí cảm biến đo biến dạng/áp suất tại xưởng Z49 và phòng LAS-XD72, cho phép bất kỳ nhóm nghiên cứu nào cũng có thể tái lập hoàn toàn kết quả mô phỏng và thực nghiệm.

5. Lộ trình phát triển nghiên cứu 10 năm tiếp theo được định hình ra sao?

Lộ trình hướng tới việc số hóa toàn diện máy khoan MK-Z49: giai đoạn 1 (1-3 năm) tích hợp bài toán đàn hồi khâu FEM-MBS; giai đoạn 2 (3-5 năm) thiết kế hệ thống điều khiển van tỷ lệ servo điện - thủy lực thích nghi; giai đoạn 3 (5-10 năm) phát triển robot khoan ngầm tự hành thông minh không người lái (Autonomous Tunnelling Jumbo) định vị quang học trong hầm quân sự.

Kết luận

  1. Xây dựng hoàn chỉnh mô hình toán học động học thuận - ngược bằng phương pháp Denavit - Hartenberg cho hệ giá đỡ máy khoan MK-Z49, xác định chính xác cấu hình 6 bậc tự do giải quyết bài toán không gian làm việc trên mặt gương hầm $2,8\text{ m} \times 1,8\text{ m}$.
  2. Thiết lập hệ phương trình vi phân - đại số (DAE) động lực học hệ nhiều vật mạch vòng kín có xét đến tính nén của dầu thủy lực và phản lực đàn nhớt đất đá Hsieh - Reissner dưới tải trọng va đập $60\text{ Hz}$.
  3. Đề xuất giải pháp cải tiến kết cấu đột phá: bổ sung khớp quay $q_4$ góc quay $90^\circ - 270^\circ$ tại đầu cần và tối ưu hóa đường kính cụm xi lanh nâng hạ, dập tắt biên độ rung lắc điểm $S$ xuống dưới $0,5\text{ mm}$.
  4. Kiểm chứng thực nghiệm đa cấp độ trên các mẫu bê tông C30 - C70 và hiện trường hầm thực tế, chứng minh sai số mô hình $< 8,5%$, độ tin cậy $R^2 > 0,93$, khống chế sai lệch lỗ khoan $\Delta x \le 1,28\text{ cm} \le 1,5 L_k$ theo đúng TCVN.
  5. Mở ra nhánh nghiên cứu chuyên sâu về cơ điện tử ngầm và hệ thống máy công trình tự hành thích ứng cao phục vụ đắc lực cho quốc phòng và hiện đại hóa hạ tầng đất nước.