Luận án TS Nguyễn Văn Mùi: Mô phỏng số, tối ưu hệ thống cấp phôi rung động

Tổng quan về luận án

Trong kỷ nguyên tự động hóa công nghiệp và sản xuất thông minh, các nguyên công lắp ráp, đóng gói và chiết rót đòi hỏi sự đồng bộ hóa cao độ về năng suất và độ chính xác. Đặc biệt, trong ngành sản xuất dược phẩm vô trùng, việc đóng nắp lọ vaccine yêu cầu môi trường vô khuẩn nghiêm ngặt cấp độ A, loại bỏ hoàn toàn sự can thiệp trực tiếp của con người theo tiêu chuẩn GMP-WHO. Thiết bị cấp phôi tự động theo nguyên lý rung động (Vibratory Bowl Feeder - VBF) với máng dẫn xoắn ốc Acsimet đóng vai trò then chốt trong việc phân loại, định hướng và vận chuyển các chi tiết rời rạc có kích thước nhỏ. Tuy nhiên, phương pháp thiết kế truyền thống chủ yếu dựa trên thử nghiệm - sai lỗi (trial-and-error), dẫn đến chi phí chế tạo đơn chiếc tốn kém, thời gian căn chỉnh kéo dài và khó đạt được trạng thái vận hành tối ưu ở chế độ cộng hưởng.

Luận án tiến sĩ kỹ thuật cơ khí với đề tài "Nghiên cứu ứng dụng mô phỏng số để khảo sát động lực học và tối ưu kết cấu điều hướng trong hệ thống cấp phôi tự động theo nguyên lý rung động" do nghiên cứu sinh Nguyễn Văn Mùi thực hiện dưới sự hướng dẫn khoa học của PGS. Lê Giang Nam và TS. Bùi Quí Lực tại Trường Đại học Bách khoa Hà Nội (2019, Mã số ngành: 9520103) đã giải quyết triệt để rào cản kỹ thuật này.

Khoảng trống nghiên cứu (Research Gap) then chốt được xác định: Trong khi các công trình quốc tế thường khảo sát phôi kim loại cứng hoặc bán dẫn, các nghiên cứu chuyên sâu về động lực học tiếp xúc ma sát đàn hồi của phôi phi kim (cụ thể là nắp chai vaccine bằng cao su với hệ số ma sát cao) di chuyển trên máng Inox SUS 304 trong phễu xoắn ốc cộng hưởng còn rất hạn chế. Đồng thời, sự thiếu vắng của một quy trình tích hợp khép kín từ tham số hóa hình học (Parametric CAD), phân tích dao động riêng (FEA Modal Analysis), mô phỏng động lực học đa vật thể (Multibody Dynamics - MBD) đến quy hoạch thực nghiệm tối ưu hóa (Design of Experiments - DOE) là nguyên nhân chính khiến công nghiệp chế tạo máy trong nước phụ thuộc vào nhập khẩu.

Hệ thống câu hỏi nghiên cứu và giả thuyết khoa học:

1. RQ1: Bản chất quy luật động lực học tiếp xúc giữa phôi cao su đàn hồi và máng xoắn Inox SUS 304 thay đổi như thế nào dưới tác động đồng thời của góc nghiêng rãnh xoắn $\theta$, góc rung động $\psi$ và điện áp kích thích $U$?
   • H1: Tồn tại một ngưỡng góc nghiêng rãnh xoắn tới hạn $\theta_{opt}$ mà tại đó vận tốc cấp phôi đạt cực trị trước khi chuyển sang trạng thái trượt lùi hoặc va đập hỗn loạn.
2. RQ2: Làm thế nào để thiết lập mô hình tính toán số có độ tin cậy cao, thay thế các thử nghiệm vật lý tốn kém nhằm xác định tần số dao động riêng $f_r$ trùng với tần số kích từ $f = 100\text{ Hz}$?
   • H2: Quy trình thiết kế tham số hóa kết hợp phương pháp phần tử hữu hạn (FEM) cho phép điều chỉnh chủ động độ cứng của hệ thống 03 cụm lò xo lá bố trí lệch nhau $120^\circ$, đảm bảo sai số tần số dao động riêng dưới $3%$ so với thực tế.
3. RQ3: Cấu hình và thứ tự bố trí các bẫy cơ học điều hướng ảnh hưởng như thế nào đến tỷ lệ cấp phôi đúng hướng và tổng thời gian cấp nắp vaccine?
   • H3: Tối ưu hóa hình học chuỗi 03 bẫy định hướng trên môi trường mô phỏng ảo giúp loại bỏ $100%$ các trạng thái phôi sai hướng (nằm ngửa, nằm nghiêng) và đạt năng suất cấp liệu ổn định tối đa.

Khung lý thuyết nền tảng tích hợp: Động lực học hệ nhiều vật (MBD), lý thuyết ma sát tiếp xúc Coulomb cải biên, động lực học rung phi tuyến, phương pháp phân tích năng lượng Rayleigh-Ritz và phương trình vi phân biến phân Lagrange.

Đóng góp đột phá định lượng: Xây dựng thành công phương trình hồi quy bậc hai mô tả tương tác đa biến giữa điện áp và góc nghiêng máng; xác định chính xác góc nghiêng tối ưu $\theta = 1{,}256^\circ$ tại điện áp $U = 170\text{V}$ cho vận tốc di chuyển cực đại; giảm thiểu hơn $70%$ thời gian thiết kế - thử nghiệm chế tạo thực tế. Phạm vi nghiên cứu bao quát từ kiểm chuẩn tham số vật liệu thực tế trên máy đo ma sát VF, đo gia tốc rung 3 trục bằng thiết bị chuyên dụng FAG Detector III, NACENTECH NI, VIBROPOST 80 ReX, đến kiểm nghiệm mô phỏng động lực học 3D với dải điện áp $100\text{V} - 170\text{V}$ và dải góc nghiêng $1{,}1^\circ - 4{,}5^\circ$.

Literature Review và Positioning

Lịch sử nghiên cứu cơ cấu cấp phôi rung động ghi nhận sự phát triển liên tục từ các mô hình giải tích một bậc tự do đến các mô phỏng động lực học tương tác phi tuyến phức tạp:

graph TD
    A["Mô hình giải tích cơ bản<br>Trần Văn Địch (2000), Boothroyd (1982)"] --> B["Khảo sát thông số động lực học đơn lẻ<br>Lee (2002), Berkowitz & Canny (1996)"]
    B --> C["Mô hình hóa số FEM/DEM & MBD<br>Ziaei-Rad (2013), Mucchi (2011), Chen (2011)"]
    C --> D["Quy trình thiết kế chủ động cộng hưởng<br>& Tối ưu hóa tiếp xúc phi kim<br><b>Nguyễn Văn Mùi (2019)</b>"]

Tổng quan các trường phái nghiên cứu quốc tế và trong nước

1. Trường phái giải tích và cơ học thực nghiệm truyền thống: Trần Văn Địch (2000) và Nguyễn Phương & Nguyễn Thị Phương Giang (2003) đặt nền móng lý thuyết cấp phôi rung tại Việt Nam, thiết lập các điều kiện cân bằng lực quán tính và ma sát $F_{qt} > F_{ms} + F_{ng}$. Tuy nhiên, các phân tích này chủ yếu dừng lại ở bài toán phẳng 2D, coi phôi là chất điểm hoặc vật rắn tuyệt đối, chưa tính đến biến dạng đàn hồi cục bộ và ma sát động thay đổi.
2. Trường phái khảo sát động lực học thực nghiệm và mô phỏng số:
   • Lee (2002) và Han & Lee (2004) khảo sát ảnh hưởng của biên độ rung đến vận tốc di chuyển với tần số $f = 120\text{ Hz}$, hệ số ma sát $\mu = 0{,}5$, góc nghiêng lò xo $\alpha = 18^\circ$, góc nghiêng máng $\theta = 5^\circ$. Các tác giả nhận thấy khi biên độ đạt $0{,}061\text{ mm}$, vận tốc phôi đạt $0{,}0937\text{ m/s}$ (sai lệch $2{,}7%$ so với thực nghiệm). Khi tăng số lượng phôi lên 200, 400, 800 phôi, trạng thái chuyển động hỗn loạn xuất hiện làm suy giảm tính quy luật của hệ số điều khiển.
   • Berkowitz & John Canny (1996) đánh giá xác suất định hướng của 350 phôi kim loại qua 6 trạng thái không gian với góc máng $12{,}8^\circ$ tại biên độ $0{,}6\text{ mm}$ và $1{,}0\text{ mm}$. M. Jiang et al. (2003) và Daniel Hofmann et al. (2013) tiếp tục phân tích tỷ lệ kích thước bẫy $b/w$, chứng minh khi tỷ lệ $b/w > 1{,}1$ thì xác suất phôi qua bẫy đạt $100%$.
   • Paul C. Chao & Chien-Yu Shen (2007) áp dụng phương pháp phân ly Rayleigh-Ritz và phương trình Lagrange mô tả 3 mode dao động đầu tiên của băng tải rung thẳng, xác định hiện tượng cộng hưởng ở tần số $330\text{ Hz}$.
   • Emiliano Mucchi et al. (2011) khảo sát động học đàn hồi của phễu rung với 2 cấu hình bố trí nam châm điện (phương tiếp tuyến và phương thẳng đứng), chỉ ra rằng gia tốc tại vành phễu trên luôn đạt cực đại.
   • H. Ziaei-Rad (2013) ứng dụng phương pháp phần tử rời rạc (DEM) 2D để mô phỏng hạt thép ($\rho = 7800\text{ kg/m}^3, f = 60\text{ Hz}, \mu = 0{,}5, \theta = 5^\circ, \psi = 18^\circ$), cho thấy mức sai lệch giữa mô phỏng và lý thuyết là $10%$.
   • J. Leberle (2011) nghiên cứu băng trượt rung liên tục, ghi nhận sự sai lệch vận tốc lên tới $20%$ khi hệ số ma sát $\mu = 0{,}25$ và góc pha thay đổi từ $0^\circ - 360^\circ$. Đáng chú ý, Leberle nhấn mạnh khi $\mu > 0{,}8$, mô hình giải tích hoàn toàn mất tính tương đồng với thực tế do hiệu ứng kẹt trượt đàn hồi (stick-slip).

Tranh luận khoa học và định vị của luận án

Tồn tại hai quan điểm trái ngược trong y văn: Một trường phái cho rằng việc tăng biên độ rung và góc nghiêng máng luôn tỷ lệ thuận với việc gia tăng năng suất cấp phôi (Maul, 1997; Chen, 2011). Ngược lại, trường phái thứ hai (Hao-Wen Ma & Gang Fang, 2015; Leberle, 2011) chứng minh rằng khi gia tốc pháp tuyến vượt ngưỡng $A_n/g_n > 1{,}0$, phôi sẽ chuyển sang chế độ "bay" (ném tự do), gây va đập hỗn loạn làm mất ổn định dòng phôi và phá vỡ cấu trúc định hướng tại các bẫy cơ học.

Luận án của Nguyễn Văn Mùi định vị tại điểm giao thoa giải quyết tranh luận này: Đối với phôi cao su đàn hồi có hệ số ma sát rất cao tiếp xúc Inox SUS 304, không thể áp dụng các tham số của phôi kim loại. Tác giả tiên phong phát triển phương pháp "Thiết kế chủ động cộng hưởng" (Active Resonance Design), kết hợp tối ưu hóa đa biến thực nghiệm - số, giúp duy trì gia tốc máng dẫn trong vùng cận giới hạn ổn định, tối đa hóa vận tốc mà không làm mất kiểm soát quỹ đạo phôi qua bẫy.

Đó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 cơ cấu cấp phôi rung của Boothroyd và Trần Văn Địch thông qua việc mô hình hóa chính xác tương tác tiếp xúc phi tuyến của vật liệu mềm đàn hồi cao.

                  ┌────────────────────────────────────────────────────────┐
                  │             ĐIỀU KIỆN ĐỘNG LỰC HỌC TỔNG QUÁT           │
                  │             Trượt lên: F_qt > F_ms + m_p * g * sinθ    │
                  └───────────────────────────┬────────────────────────────┘
                                              │
                     ┌────────────────────────┴────────────────────────┐
                     ▼                                                 ▼
     ┌───────────────────────────────┐                 ┌───────────────────────────────┐
     │       CHẾ ĐỘ TRƯỢT (SLIDING)  │                 │        CHẾ ĐỘ NHẢY (HOPPING)  │
     │      N > 0 ; A_n / g_n < 1,0  │                 │        N = 0 ; A_n / g_n ≥ 1,0│
     │  Vận tốc phụ thuộc ma sát μ_s │                 │  Phôi bay rời khỏi mặt máng   │
     └───────────────┬───────────────┘                 └───────────────┬───────────────┘
                     │                                                 │
                     └────────────────────────┬────────────────────────┘
                                              │
                                              ▼
                  ┌────────────────────────────────────────────────────────┐
                  │          VÙNG VẬN HÀNH TỐI ƯU CẬN CỘNG HƯỞNG           │
                  │   f_r = 100 Hz ; θ_opt = 1,256° ; U = 170V ; v_max    │
                  └────────────────────────────────────────────────────────┘

1. Chuẩn hóa hệ phương trình vi phân chuyển động 3D: Lực quán tính kích động được phân rã chính xác theo phương pháp tuyến ($n$) và phương tiếp tuyến ($p$) với máng xoắn Acsimet: $$F_{qt_p} = m_p \cdot a_0 \omega^2 \cos\psi$$ $$F_{qt_n} = m_p \cdot a_0 \omega^2 \sin\psi$$ Điều kiện để phôi trượt tịnh tiến đi lên không gián đoạn: $$m_p \cdot a_0 \omega^2 \cos\psi \ge \mu_s \cdot N + m_p g \sin\theta$$ với phản lực pháp tuyến tức thời: $$N = m_p g \cos\theta - m_p \cdot a_0 \omega^2 \sin\psi$$ Điều kiện xảy ra trạng thái phôi rời khỏi lòng máng (chế độ nhảy đàn hồi): $$N \le 0 \iff a_0 \omega^2 \sin\psi \ge g \cos\theta \iff \frac{A_n}{g_n} \ge 1{,}0$$

2. Mô hình hóa tham số kích từ điện từ kép: Nam châm điện xoay chiều cấp nguồn $50\text{ Hz}$ tạo ra xung lực hút điện từ biến thiên tuần hoàn với tần số gấp đôi: $f = 2 \times 50\text{ Hz} = 100\text{ Hz}$. Cơ hệ được chuẩn hóa thành hệ dao động cưỡng bức phi tuyến có cản, với độ cứng tương đương được kiểm soát bởi 3 cụm nhíp lò xo lá đặt nghiêng góc $\alpha$.

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

Khung phân tích của luận án tích hợp liên hoàn 4 công cụ chuyên sâu:

• CAD Tham số hóa (CATIA V5): Thiết lập bảng liên kết hàm số hình học từ MS Excel đến mô hình 3D, cho phép tự động cập nhật toàn bộ kích thước hình học của phễu, máng xoắn, góc nghiêng và giá đỡ khi thay đổi thông số đầu vào.
• Phân tích phần tử hữu hạn FEM (ANSYS Workbench): Khảo sát dạng dao động (Mode Shapes) và trích xuất chính xác 6 tần số dao động riêng đầu tiên của toàn bộ cơ hệ phễu - nhíp đàn hồi.
• Mô phỏng động lực học đa vật thể MBD (MSC ADAMS): Gán hàm động học tương đương xác định từ thực nghiệm, thiết lập mô hình tiếp xúc va đập phi tuyến Hertz - Coulomb với hệ số ma sát thực đo.
• Quy hoạch thực nghiệm và hồi quy phi tuyến (MATLAB Optimization Toolbox): Xây dựng ma trận thực nghiệm trực giao, xác định phương trình bề mặt đáp ứng và tối ưu hóa hàm mục tiêu năng suất cấp phôi.

Điều kiện biên giới hạn (Boundary Conditions): Góc nghiêng rãnh xoắn $\theta \in [1{,}1^\circ; 4{,}5^\circ]$, điện áp kích từ $U \in [100\text{V}; 170\text{V}]$, tần số kích thích cố định $100\text{ Hz}$, phôi nắp vaccine bằng cao su lưu hóa có kích thước danh nghĩa tiêu chuẩn y tế, vật liệu phễu Inox SUS 304 xử lý bề mặt độ bóng cao.

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ủ chặt chẽ thế giới quan thực chứng (Positivism) kết hợp phương pháp luận kỹ thuật hiện thực phê phán, sử dụng thiết kế đa phương pháp tích hợp (Integrated Multi-method Design) bao gồm: Tính toán giải tích $\rightarrow$ Mô hình hóa tham số $\rightarrow$ Mô phỏng phần tử hữu hạn $\rightarrow$ Thực nghiệm kiểm chuẩn ma sát/dao động $\rightarrow$ Mô phỏng động lực học MBD $\rightarrow$ Tối ưu hóa bề mặt đáp ứng.

flowchart LR
    subgraph G1["Giai đoạn 1: Tham số hóa & FEM"]
        A1["CATIA + Excel<br>Parametric Modeling"] --> A2["ANSYS Workbench<br>Modal Analysis (fr = 100 Hz)"]
    end
    subgraph G2["Giai đoạn 2: Đo lường thực nghiệm"]
        B1["Máy đo VF<br>Hệ số ma sát thực μ"] --> B3["Đo rung 3 trục<br>FAG, NI, Vibropost 80"]
        B2["Mẫu phôi cao su<br>+ Máng SUS 304"] --> B3
    end
    subgraph G3["Giai đoạn 3: MBD & Tối ưu"]
        A2 & B3 --> C1["MSC ADAMS<br>Multibody Simulation"]
        C1 --> C2["MATLAB Response Surface<br>Tối ưu θ=1,256°, U=170V"]
        C2 --> C3["Kiểm nghiệm thực tế<br>Hệ thống nắp vaccine"]
    end

Quy trình nghiên cứu rigorous

1. Giao thức đo hệ số ma sát thực tế: Sử dụng máy đo ma sát chuyên dụng Model VF. Mẫu thử gồm vật liệu cao su của nắp chai vaccine trượt trên đế tấm Inox SUS 304 trong điều kiện độ ẩm và nhiệt độ phòng chuẩn. Dữ liệu lực ma sát tiếp xúc được ghi nhận theo thời gian thực để trích xuất chính xác hệ số ma sát tĩnh $\mu_s$ và ma sát động $\mu_k$.
2. Giao thức đo lường dao động và trích xuất hàm chuyển động:
   • Cảm biến gia tốc áp điện đa trục kết nối với bộ phân tích Detector III (FAG) và phần mềm phân tích Trendline.
   • Hệ thống thu thập dữ liệu chuyên dụng của Viện Ứng dụng Công nghệ (NACENTECH) tích hợp phần mềm LabVIEW/NI.
   • Thiết bị phân tích dao động VIBROPOST 80 kết hợp phần mềm xử lý tín hiệu ReX.
   • Kết quả đo thực nghiệm tại các mức điện áp từ $100\text{V}$ đến $170\text{V}$ cung cấp chính xác biên độ dao động tịnh tiến thẳng đứng $z(t)$ và biên độ dao động góc xoay quanh trục $\varphi(t)$, làm cơ sở xây dựng hàm động học tương đương (Function Builder) nạp vào MSC ADAMS.
3. Phân tích dạng dao động trên ANSYS Workbench:
   • Thiết lập mô hình vật liệu Inox SUS 304: Mô đun đàn hồi $E = 193\text{ GPa}$, hệ số Poisson $\nu = 0{,}31$, khối lượng riêng $\rho = 7930\text{ kg/m}^3$.
   • Gán điều kiện biên ngàm cố định tại mặt dưới của 3 chân đế lò xo lá.
   • Khảo sát các Mode dao động từ Mode 1 đến Mode 6. Kết quả xác định Mode dao động xoắn - tịnh tiến đồng thời có tần số dao động riêng $f_r \approx 100\text{ Hz}$, trùng khớp hoàn toàn với tần số lực kích từ của nam châm điện $50\text{ Hz}$ hai nhịp.

Data và phân tích

Quá trình quy hoạch thực nghiệm được tiến hành nhằm thiết lập mối quan hệ giữa hai biến độc lập: Điện áp cấp cho nam châm điện $X_1 = U\text{ (V)}$ và Góc nghiêng rãnh xoắn $X_2 = \theta\text{ (độ)}$ đối với biến mục tiêu là Vận tốc cấp phôi $Y = v\text{ (mm/s)}$.

Ma trận dữ liệu thực nghiệm và mô phỏng trên MSC ADAMS được xử lý bằng phần mềm MATLAB, áp dụng mô hình hồi quy phi tuyến bậc hai dạng: $$v = \beta_0 + \beta_1 U + \beta_2 \theta + \beta_{12} U \theta + \beta_{11} U^2 + \beta_{22} \theta^2$$

Hệ số xác định tương quan đạt $R^2 > 0{,}96$, kiểm định Fisher ($F$-test) và giá trị $p < 0{,}001$ khẳng định mô hình hồi quy có ý nghĩa thống kê vượt trội, độ tin cậy tuyệt đối trong dải thông số khảo sát.

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

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

1. Xác lập chính xác điểm cực trị hình học và động lực học: Mô phỏng và thực nghiệm chứng minh vận tốc di chuyển của nắp vaccine cao su đạt giá trị cực đại khi góc nghiêng rãnh xoắn $\theta = 1{,}256^\circ$ tại điện áp cực hạn $U = 170\text{V}$. Khi góc $\theta < 1{,}1^\circ$, hành trình nâng phôi bị kéo dài không cần thiết; ngược lại khi $\theta > 2{,}5^\circ$, vận tốc tụt dốc nhanh chóng do thành phần lực cản trọng trường $m_p g \sin\theta$ vượt quá lực ma sát truyền động quán tính.
2. Khám phá hiện tượng mất ổn định định hướng khi vượt áp: Tại mức điện áp $U > 160\text{V} - 170\text{V}$, biên độ dao động thẳng đứng tăng cao khiến tỷ số gia tốc $A_n/g_n \ge 1{,}0$, phôi chuyển từ chế độ trượt êm sang chế độ ném nhảy liên tục. Đối với phôi cao su có độ đàn hồi cao, va đập này tạo ra các dao động xoay ngẫu nhiên, khiến phôi bị văng khỏi máng dẫn hoặc lật sai hướng trước khi đi vào cụm bẫy định hướng.
3. Quy luật phân loại của chuỗi 03 bẫy cơ học điều hướng: Luận án thiết kế và thử nghiệm hệ thống 03 bẫy định hướng chuyên biệt cho nắp chai vaccine:
   • Bẫy 1 (Bẫy loại phôi chồng/nghiêng): Sử dụng thanh gạt giới hạn chiều cao, phôi nằm sai tư thế bị đẩy rơi tự do trở lại lòng đáy phễu.
   • Bẫy 2 (Bẫy định hướng đáy nắp): Lợi dụng sự lệch trọng tâm hình học của nắp chai cao su rỗng, giữ lại phôi có hướng đáy quay xuống dưới và lật/loại phôi quay đáy lên trên.
   • Bẫy 3 (Bẫy thu hẹp rãnh thoát): Chuẩn hóa dòng phôi đơn chiếc đi vào máng dẫn thẳng tiếp nhận. Kết quả mô phỏng sau 30 giây đầu tiên và đánh giá tổng thời gian cấp phôi khẳng định cấu hình bẫy tối ưu cho phép $100%$ phôi đầu ra đạt chuẩn tư thế lắp ráp với thời gian cấp liệu ngắn nhất.
4. Mức độ tương đồng cao giữa mô phỏng MBD và thực nghiệm: Sai số tương đối giữa vận tốc tính toán bằng mô phỏng số MSC ADAMS và kết quả đo thực nghiệm trên mô hình thực tế dao động trong khoảng hẹp từ $3{,}8%$ đến $7{,}2%$. Con số này thấp hơn đáng kể so với mức sai lệch $10% - 20%$ trong các nghiên cứu quốc tế trước đây của Ziaei-Rad (2013) hay Leberle (2011).

Implications đa chiều

• Ý nghĩa lý thuyết: Cung cấp cơ sở khoa học hoàn chỉnh về động lực học tiếp xúc ma sát dị hướng của vật liệu phi kim đàn hồi trong trường kích rung xoắn ốc 3D; bổ sung dữ liệu thực nghiệm quan trọng cho lý thuyết dao động phi tuyến máy cơ khí.
• Đổi mới phương pháp luận: Chuẩn hóa quy trình "Digital Twin" sơ khai cho thiết bị cấp phôi rung: Tham số hóa CATIA $\rightarrow$ Modal ANSYS $\rightarrow$ Dynamics ADAMS $\rightarrow$ MATLAB Optimization. Phương pháp này có khả năng chuyển giao và áp dụng trực tiếp cho nhiều dạng phôi công nghiệp phức tạp khác (vi mạch điện tử, nút bấm, linh kiện nhựa y tế).
• Ứng dụng thực tiễn và chính sách: Giúp các doanh nghiệp chế tạo máy trong nước làm chủ công nghệ thiết kế phễu rung cấp phôi chính xác "đúng ngay từ lần đầu tiên" (Right First Time), cắt giảm triệt để chi phí gia công thử nghiệm lặp lại, thúc đẩy nội địa hóa các dây chuyền đóng gói dược phẩm tự động đạt chuẩn quốc tế.

Limitations và Future Research

Mặc dù đạt được những kết quả đột phá, luận án thẳng thắn chỉ ra các giới hạn nghiên cứu:

1. Giới hạn về mô hình hóa biến dạng đàn hồi: Mô phỏng trong MSC ADAMS chủ yếu coi phôi nắp vaccine là vật thể rắn có gán tiếp xúc đàn hồi cục bộ (Contact Force Function), chưa tích hợp đầy đủ biến dạng phi tuyến phi hình học dạng phần tử hữu hạn phi tuyến (FE-Solid Non-linear) trong thời gian thực do hạn chế về tài nguyên tính toán.
2. Điều kiện môi trường lý tưởng: Các thử nghiệm thực nghiệm được thực hiện trong môi trường phòng thí nghiệm có nhiệt độ và độ ẩm tương đối ổn định. Trong thực tế nhà máy dược phẩm vô trùng, sự thay đổi độ ẩm không khí và hiện tượng tích tụ tĩnh điện trên bề mặt Inox SUS 304 có thể làm biến thiên nhẹ hệ số ma sát $\mu_s$.
3. Phạm vi hình học của đối tượng phôi: Luận án tập trung chuyên sâu vào đối tượng nắp chai vaccine dạng trụ tròn rỗng bậc một đầu kín.

Chương trình nghiên cứu 10 năm tiếp theo (Future Research Agenda)

• Phát triển mô hình mô phỏng đồng thời biến dạng - động lực học (Co-simulation giữa ANSYS Explicit Dynamics và MSC ADAMS) để mô tả chính xác 100% hiện tượng biến dạng đàn nhớt (viscoelasticity) của cao su khi va đập ở tần số cao.
• Nghiên cứu tích hợp hệ thống điều khiển thích nghi thông minh (Adaptive Resonance Control) sử dụng cảm biến phản hồi quang học và biến tần tự động dò quét tần số cộng hưởng theo thời gian thực khi khối lượng phôi trong phễu suy giảm dần.
• Mở rộng ngân hàng dữ liệu tối ưu hóa cho các dạng hình học phôi phi đối xứng, phôi tổ hợp kim loại - nhựa trong ngành công nghiệp bán dẫn và sản xuất thiết bị y tế tinh vi.

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

Tác động học thuật và chỉ số trích dẫn

Công trình đóng góp một tài liệu tham khảo mẫu mực về ứng dụng mô phỏng số kỹ thuật cao trong cơ học máy tại Việt Nam. Phương pháp luận kết hợp ANSYS - ADAMS - MATLAB mở ra hướng nghiên cứu mới cho các nghiên cứu sinh, học viên cao học chuyên ngành Cơ kỹ thuật, Cơ điện tử và Kỹ thuật Chế tạo máy, với tiềm năng trích dẫn cao trong các tạp chí chuyên ngành cơ khí và tự động hóa.

Chuyển đổi công nghiệp và giá trị kinh tế - xã hội

• Ngành công nghiệp Dược phẩm & Y tế: Cung cấp giải pháp công nghệ tự chủ để chế tạo các hệ thống cấp nắp chai tự động cho dây chuyền đóng lọ vaccine vô trùng, giảm giá thành đầu tư thiết bị nhập ngoại từ Đức, Ý, Nhật Bản tới $50% - 60%$.
• Ngành Công nghiệp Phụ trợ & Lắp ráp Điện tử: Ứng dụng quy trình thiết kế chủ động cộng hưởng giúp các nhà sản xuất máy tự động hóa tại Việt Nam rút ngắn chu kỳ nghiên cứu - phát triển sản phẩm (R&D) từ vài tháng xuống còn vài ngày, nâng cao năng lực cạnh tranh quốc gia trong chuỗi cung ứng toàn cầu.

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

mindmap
  root((Hệ sinh thái hưởng lợi))
    Nghiên cứu sinh & Giảng viên
      Cơ sở lý thuyết tiếp xúc ma sát phi tuyến
      Quy trình nghiên cứu kết hợp FEM MBD DOE
      Tài liệu đào tạo chuyên sâu Cơ điện tử
    Kỹ sư R&D Doanh nghiệp
      Quy trình thiết kế tham số hóa CATIA-Excel
      Rút ngắn hơn 70% thời gian căn chỉnh
      Đảm bảo máy chạy đúng dải cộng hưởng 100 Hz
    Nhà máy Dược phẩm
      Làm chủ thiết bị chiết rót vô trùng GMP
      Cắt giảm 50-60% chi phí đầu tư máy ngoại
      Vận hành ổn định năng suất cao
    Cơ quan Quản lý Khoa học
      Cơ sở hoạch định chính sách tự chủ công nghệ
      Tiêu chuẩn hóa thiết bị tự động hóa y tế

• Nghiên cứu sinh và Giới học thuật: Tiếp cận phương pháp luận nghiên cứu thực nghiệm kết hợp mô phỏng số tiên tiến; kế thừa khung phân tích động lực học tiếp xúc phi tuyến để phát triển các đề tài chuyên sâu.
• Kỹ sư R&D và Thiết kế Cơ khí: Nắm vững công cụ tham số hóa và tối ưu hóa kết cấu máng xoắn, làm chủ kỹ thuật điều chỉnh nhíp đàn hồi để hệ thống luôn vận hành tại điểm cộng hưởng cơ học bền vững.
• Doanh nghiệp Sản xuất Dược phẩm: Sở hữu dây chuyền tự động cấp nắp vaccine vận hành tin cậy, đạt tiêu chuẩn phòng sạch vô khuẩn nghiêm ngặt, loại bỏ hoàn toàn nguy cơ lây nhiễm chéo do thao tác thủ công.

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

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

Đóng góp lý thuyết độc đáo nhất là việc mô hình hóa chính xác quy luật động lực học tiếp xúc ma sát dị hướng của phôi đàn hồi phi kim (cao su) trên máng xoắn Inox SUS 304, mở rộng trực tiếp lý thuyết cấp phôi rung động cổ điển của Geoffrey Boothroyd (1982) và Trần Văn Địch (2000). Thay vì coi ma sát là hằng số đơn giản trên mặt phẳng nghiêng, luận án đã tích hợp trường gia tốc dao động xoắn 3D với hệ số ma sát thực nghiệm biến thiên, xác lập điều kiện biên phân tách giữa trạng thái trượt ổn định và trạng thái ném nhảy hỗn loạn.

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

So với nghiên cứu của H. Ziaei-Rad (2013) (chỉ mô phỏng 2D rời rạc bằng DEM) và J. Leberle (2011) (khảo sát băng trượt thẳng với độ sai lệch thực nghiệm lên tới $20%$), luận án đã thiết lập một chuỗi quy trình số hóa khép kín 3D hoàn chỉnh (Seamless Digital Simulation Chain): Tham số hóa hình học từ Excel sang CATIA $\rightarrow$ Phân tích tần số dao động riêng trên ANSYS Workbench $\rightarrow$ Mô phỏng động lực học tiếp xúc đa vật thể trên MSC ADAMS dựa trên hàm kích động thực nghiệm $\rightarrow$ Tối ưu hóa hàm hồi quy phi tuyến đa biến trên MATLAB. Kết quả đưa sai số mô phỏng so với thực tế xuống mức kỷ lục: $3{,}8% - 7{,}2%$.

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

Phát hiện bất ngờ nhất là vận tốc di chuyển cực đại của phôi cao su đạt được tại góc nghiêng rãnh xoắn rất nhỏ $\theta = 1{,}256^\circ$, hoàn toàn trái ngược với khuyến nghị thiết kế truyền thống dành cho chi tiết kim loại thông thường ($\theta = 3^\circ - 5^\circ$). Dữ liệu thực nghiệm và mô phỏng chỉ ra rằng do phôi cao su có hệ số ma sát lớn và tính đàn hồi cao, khi góc nghiêng $\theta > 2{,}5^\circ$, thành phần phản lực pháp tuyến suy giảm kết hợp với lực cản dốc khiến phôi nhanh chóng bị trượt ngược, làm giảm sút năng suất nghiêm trọng.

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

Có. Luận án cung cấp chi tiết toàn bộ quy trình tái lập bao gồm:

• Thông số vật liệu: Inox SUS 304 ($E = 193\text{ GPa}, \rho = 7930\text{ kg/m}^3, \nu = 0{,}31$).
• Kích thước kết cấu phễu rung, góc nghiêng lò xo lá $\alpha$, sơ đồ bố trí 3 cụm nhíp $120^\circ$.
• Giao thức đo kiểm trên máy đo ma sát VF, thiết bị FAG Detector III, NACENTECH NI, VIBROPOST 80.
• Các bước gán liên kết tiếp xúc (Impact Contact Parameters: Stiffness, Force Exponent, Damping, Friction Coefficients) trong môi trường MSC ADAMS.

5. Định hướng nghiên cứu 10 năm tới được vạch ra như thế nào?

Định hướng tập trung vào việc xây dựng hệ thống cấp phôi rung thông minh thế hệ mới (Smart Cyber-Physical Vibratory Feeder): Ứng dụng cảm biến đo gia tốc và thị giác máy tính (Computer Vision) để giám sát lưu lượng phôi; tích hợp mạch điều khiển biến tần tự động thích nghi (Auto-tuning Inverter) giúp duy trì tần số cộng hưởng $f_r$ ổn định tuyệt đối khi khối lượng phôi thay đổi liên tục trong quá trình sản xuất; mở rộng mô phỏng tương tác đa trường vật lý cho các linh kiện siêu nhỏ (Micro/Nano Assembly).

Kết luận

Luận án tiến sĩ của tác giả Nguyễn Văn Mùi là công trình nghiên cứu khoa học công phu, nghiêm túc và mang tính đột phá cao trong lĩnh vực Kỹ thuật Cơ khí và Tự động hóa chế tạo máy. Toàn bộ nội dung công trình được đúc kết qua các thành tựu cốt lõi:

1. Xây dựng thành công phương pháp thiết kế chủ động cộng hưởng: Xác định chính xác kết cấu và thông số đàn hồi của hệ thống lò xo lá để tần số dao động riêng của cơ hệ phễu đạt đúng $f_r = 100\text{ Hz}$, trùng với tần số lực kích từ nam châm điện hai nhịp, triệt tiêu hoàn toàn hiện tượng lệch pha năng lượng.
2. Làm chủ quy trình mô phỏng số tích hợp đa nền tảng: Kết hợp thành công CATIA - ANSYS Workbench - MSC ADAMS - MATLAB, thiết lập mô hình thử nghiệm ảo với độ chính xác cao (sai số kiểm chứng thực nghiệm chỉ từ $3{,}8% - 7{,}2%$).
3. Tìm ra bộ thông số công nghệ tối ưu cho nắp chai vaccine cao su: Xác định góc nghiêng rãnh xoắn tối ưu $\theta = 1{,}256^\circ$ tại điện áp $U = 170\text{V}$, thiết lập phương trình hồi quy toán học mô tả chính xác tương tác vận tốc cấp phôi.
4. Tối ưu hóa kết cấu chuỗi 03 bẫy cơ học điều hướng: Định hướng chính xác $100%$ nắp chai vaccine cao su ở đầu ra, loại bỏ hoàn toàn các tư thế phôi lỗi, đạt năng suất cấp phôi ổn định cho dây chuyền đóng gói y tế vô trùng.
5. Mở ra hướng nghiên cứu mới trong nước: Xóa bỏ sự phụ thuộc vào phương pháp chế tạo thử - sai truyền thống, cung cấp cẩm nang thiết kế chuẩn mực "đúng ngay từ lần đầu tiên", tiết kiệm hàng trăm triệu đồng chi phí R&D cho mỗi chủng loại sản phẩm.
6. Giá trị thực tiễn chuyển giao xuất sắc: Đóng góp trực tiếp vào công cuộc hiện đại hóa ngành cơ khí chế tạo máy và công nghiệp dược phẩm tự động hóa cao, khẳng định năng lực nghiên cứu độc lập, chuyên sâu và hội nhập trình độ khoa học quốc tế của đội ngũ trí thức Việt Nam.