Tổng quan về luận án

Trong kỹ thuật cơ khí động lực và ma sát học (Tribology), ổ đầu to thanh truyền trong động cơ đốt trong là chi tiết bản lề chịu tải trọng biến thiên chu kỳ phức tạp (720° góc quay trục khuỷu đối với động cơ bốn kỳ), vận tốc góc thay đổi liên tục cùng sự tương tác đa trường vật lý. Độ tin cậy và tuổi thọ của cụm ổ trượt này quyết định trực tiếp đến hiệu suất cơ học, mức tiêu hao nhiên liệu và độ bền nhiệt động của toàn bộ hệ thống động cơ. Khi vận hành dưới áp lực nổ cực đại và tốc độ cao, màng dầu bôi trơn giữa trục khuỷu và bạc lót bị trượt cắt mãnh liệt, sinh nhiệt nội tại làm suy giảm độ nhớt động lực học, dẫn tới nguy cơ sụp đổ màng dầu thủy động (hydrodynamic lubrication film collapse), gây ma sát nửa ướt, mài mòn ôvan và cào xước kim loại nghiêm trọng.

Tuy nhiên, phần lớn các công trình nghiên cứu cơ học tiếp xúc trong nước trước đây chủ yếu dừng lại ở bài toán thủy động đàn hồi đẳng nhiệt (Isothermal EHD) hoặc bôi trơn ổ chịu tải tĩnh, hoàn toàn bỏ qua tương tác nhiệt - thủy động - đàn hồi phi tuyến. Nhằm giải quyết khoảng trống học thuật này, luận án tiến sĩ mang tên "Nghiên cứu đặc tính bôi trơn nhiệt thủy động của ổ có dạng đầu to thanh truyền" do NCS. Phạm Trung Thiên thực hiện tại Bộ môn Máy và Ma sát học - Viện Cơ khí, Trường Đại học Bách khoa Hà Nội, dưới sự hướng dẫn của TS. Trần Thị Thanh Hải và TS. Phạm Minh Hải, đã tạo lập bước đột phá khi tích hợp mô phỏng số nhiệt thủy động (Thermal Hydrodynamic - THD) với thực nghiệm quang đàn hồi trên băng thử chuyên dụng B2016-BKA-20.

Cụ thể, luận án giải quyết các câu hỏi nghiên cứu cốt lõi:

  1. RQ1: Phân bố trường áp suất thủy động $p(\theta, z)$ và trường nhiệt độ màng dầu $T(\theta, y, z)$ biến thiên thế nào theo chu kỳ làm việc 720° của động cơ đốt trong dưới tác động đồng thời của tải trọng chu kỳ và tốc độ quay?
  2. RQ2: Mức độ ảnh hưởng của đặc tính lưu chất bôi trơn (độ nhớt danh định, chỉ số độ nhớt, hệ số dẫn nhiệt) đến sự gia tăng nhiệt độ cực đại $\Delta T_{max}$ và khả năng mang tải của màng dầu là bao nhiêu?
  3. RQ3: Độ chính xác và tương thích giữa mô hình số giải thuật Galerkin phần tử hữu hạn 3 chiều với dữ liệu đo đạc thực nghiệm trên thanh truyền quang đàn hồi đạt ngưỡng tin cậy nào?

Giả thuyết khoa học tương ứng bao gồm:

  • H1: Nhiệt độ màng dầu đạt giá trị cực đại đồng pha với thời điểm góc chất tải nổ cực đại (xung quanh góc quay $0^\circ / 360^\circ$ của trục khuỷu) và tạo ra gradient nhiệt độ lớn nhất tại tiết diện giữa ổ theo phương chiều dài ($z=0$).
  • H2: Sự gia tăng vận tốc góc từ 100 vg/ph đến 200 vg/ph làm tăng tổn thất ma sát nhớt, thúc đẩy tăng chênh lệch nhiệt độ $\Delta T$ màng dầu phi tuyến tính theo hàm số mũ của vận tốc trượt.

Khung lý thuyết của công trình đặt nền tảng trên hệ phương trình Reynolds mở rộng hai chiều kết hợp phương trình vi phân năng lượng nhiệt ba chiều không thứ nguyên và quy luật biến thiên độ nhớt theo nhiệt độ dạng hàm mũ của Vogel/Reynolds. Đối tượng nghiên cứu thực nghiệm là cụm ổ đầu to thanh truyền có đường kính trong $d = 97\text{ mm}$, đường kính ngoài $D = 118\text{ mm}$, chiều dài $L = 20\text{ mm}$, khe hở bán kính $C = 0.038\text{ mm}$, chế tạo bằng vật liệu quang đàn hồi epoxy PLM4/PLMH4, được khảo sát trên dải vận tốc $100\text{ vg/ph}$, $150\text{ vg/ph}$, $200\text{ vg/ph}$ (tương ứng $1.67\text{ vg/s}$ đến $3.34\text{ vg/s}$) cùng hai loại dầu chuẩn công nghiệp Besil F100 và Atox 320 qua 3000 chu kỳ thử nghiệm liên tục.

Literature Review và Positioning

Nghiên cứu về bôi trơn ổ trục động lực học đã trải qua hơn năm thập kỷ phát triển với hai trường phái chính: mô hình hóa lý thuyết số trị và thử nghiệm thực nghiệm chuyên sâu. Dòng nghiên cứu mô phỏng tải trọng tĩnh khởi xướng từ Reddy et al. (1969), Oh & Huebner (1973) khi lần đầu ứng dụng phần tử hữu hạn (FEM) để giải phương trình Reynolds cho ổ đàn hồi với giả thiết đẳng nhớt, sử dụng thuật toán Newton-Raphson nghịch đảo ma trận Jacobin nhưng bị hạn chế nghiêm trọng do tốc độ hội tụ chậm khi tỷ số lệch tâm $\varepsilon$ lớn. Tiếp đó, Fantino & Frêne (1979) phân tích ảnh hưởng của độ nhớt dầu bôi trơn đến biến dạng ổ đỡ thanh truyền bằng sai phân hữu hạn kết hợp lặp Gauss-Seidel.

Dòng nghiên cứu mô phỏng tải trọng động chứng kiến sự cạnh tranh học thuật mạnh mẽ giữa các nhóm tác giả quốc tế. Trường phái thứ nhất tập trung vào hiệu ứng đàn hồi và quán tính cơ học: Booker & Shu (1984), Goenka & Oh (1986), Aitken & McCallion (1991), Guines et al. (1995), Bonneau et al. (1995) đã tích hợp phần tử 8 nút và 20 nút 3D kết hợp thuật toán Murty để mô hình hóa vùng gián đoạn màng dầu (cavitation). Ngược lại, trường phái thứ hai nhấn mạnh hiệu ứng nhiệt thủy động (THD/TEHD) và biến dạng biên dạng: Ozasa et al. (1997), Okamoto et al. (1999), Piffeteau & Souchet (2000), Chamani et al. (2012) và Norbert Lorenz (2017) đã chứng minh rằng việc bỏ qua sinh nhiệt ma sát sẽ dẫn tới sai số dự báo áp suất cực đại lên tới $30-45%$. Đặc biệt, Bonneau & Hajjam (2001) đã ứng dụng mô hình bảo toàn dòng JFO (Jakobson-Floberg-Olsson) để phân tích ổ thanh truyền động cơ đua F1 chịu điều kiện cực hạn.

Tại Việt Nam, các công trình của Phan Thạch Hổ (2001, 2002) trên động cơ NT855, Trần Thị Thanh Hải (1999, 2014, 2018), Phạm Văn Hùng (2005) và Nguyễn Đình Tân (2018) đã giải quyết bài toán thủy động học ổ đỡ, ảnh hưởng của lực siết bu-lông và khảo sát trường áp suất thủy động. Tuy nhiên, một khoảng trống nghiên cứu cốt tử vẫn tồn tại: chưa có bất kỳ công bố nội địa nào xây dựng hoàn chỉnh thuật toán số 3D dự báo đồng thời trường nhiệt màng dầu kết hợp thiết kế hệ thống cảm biến vi cơ điện tử đo nhiệt - áp - lực thời gian thực trên thanh truyền mô hình.

                      TIẾN TRÌNH NGHIÊN CỨU BÔI TRƠN Ổ ĐẦU TO THANH TRUYỀN
  
  [1969 - 1983: Tải trọng tĩnh & Đẳng nhiệt]
  Reddy (1969), Oh & Huebner (1973), Fantino & Frêne (1979), Pierre-Eugene (1983)
  [1984 - 2005: Tải trọng động & Thủy động đàn hồi (EHD)]
  Booker & Shu (1984), Goenka (1986), Aitken (1991), Bonneau & Hajjam (2001)
  [2005 - Nay: Nhiệt thủy động đàn hồi (TEHD) & Đo đạc đa kênh]
  Fatu & Michaud (Megapascal Bench, 2005), Chamani (2012), Norbert Lorenz (2017)
  [VỊ TRÍ LUẬN ÁN: NCS. PHẠM TRUNG THIÊN (HUST - 2018/2021)]
  Tích hợp giải thuật FEM-Galerkin THD 3D + Băng thử B2016-BKA-20 cảm biến đo vô tuyến

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

  1. Băng thử Megapascal của Michaud & Fatu (2005) tại Đại học Poitiers (Pháp): Khảo sát thanh truyền thật động cơ F1 ở tốc độ lên tới $20.000\text{ vg/ph}$, tải nén $90\text{ kN}$ và tải kéo $60\text{ kN}$. Dù tái hiện chân thực điều kiện khắc nghiệt, phương pháp này đòi hỏi chi phí thiết bị hàng triệu USD và không thể trực quan hóa biến dạng quang đàn hồi bên trong thân bạc lót.
  2. Mô hình quang đàn hồi của Pierre-Eugene et al. (1983) và Optasanu (2000) tại LMS Poitiers: Sử dụng nhựa epoxy PLM4 chịu tải $60-300\text{ N}$ ở tốc độ $50-200\text{ vg/ph}$. Luận án của NCS. Phạm Trung Thiên kế thừa ưu điểm trực quan quang học nhưng vượt trội hơn nhờ tích hợp module đo nhiệt độ đa kênh NI9213, bộ truyền tín hiệu vô tuyến RF cho cảm biến áp suất màng siêu nhỏ XCQ-062 và chương trình mô phỏng số 3D hoàn chỉnh viết trên nền tảng Fortran95.

Đó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 khung lý thuyết bôi trơn thủy động cổ điển của Osborne Reynolds sang hệ điều kiện bôi trơn nhiệt thủy động (THD) phi tuyến chịu tải trọng chu kỳ động học. Mô hình lý thuyết của tác giả thiết lập sự tương tác hai chiều (two-way coupling) giữa trường thủy động và trường nhiệt độ: áp suất $p$ sinh ra trường vận tốc cắt chất lưu $(u, v, w)$, công ma sát nhớt biến thành nhiệt lượng làm tăng nhiệt độ $T$, nhiệt độ gia tăng làm suy giảm độ nhớt $\mu(T)$ theo quan hệ phi tuyến:

$$\mu(T) = \mu_0 e^{-\beta(T - T_0)}$$

Sự sụt giảm độ nhớt quay ngược lại làm thay đổi căn bản phân bố áp suất màng dầu và chiều dày màng cực tiểu $h_{min}$.

Khung mệnh đề lý thuyết được chứng minh bao gồm:

  • Proposition 1: Gradient nhiệt độ theo phương hướng kính ($\partial T / \partial y$) đóng vai trò quyết định đến ứng suất cắt ma sát nhớt tại bề mặt tiếp xúc kim loại - dầu, trong khi gradient nhiệt theo phương dọc trục ($\partial T / \partial z$) chi phối lượng nhiệt truyền qua đối lưu biên.
  • Proposition 2: Sự chuyển dịch pha nhiệt động học: đỉnh nhiệt độ màng dầu cực đại không xảy ra ở điểm có khe hở màng dầu nhỏ nhất $h_{min}$ mà định vị tại vùng hội tụ có tích số áp suất và gradient vận tốc cắt lớn nhất ($\tau = \mu \frac{\partial u}{\partial y} = \max$), tương ứng góc quay $0^\circ / 360^\circ$ tại đỉnh chu trình nổ.

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

Khung phân tích của luận án tích hợp đồng thời ba cấu trúc toán học - vật lý:

  1. Phương trình vi phân Reynolds mở rộng cho màng dầu ổ đỡ:

$$\frac{\partial}{\partial x}\left(h^3 \frac{\partial p}{\partial x}\right) + \frac{\partial}{\partial z}\left(h^3 \frac{\partial p}{\partial z}\right) = 6\mu (U_1 - U_2)\frac{\partial h}{\partial x} + 12\mu \frac{\partial h}{\partial t}$$

  1. Phương trình năng lượng nhiệt 3 chiều quy chuẩn sang hệ tọa độ không thứ nguyên $(\bar{x}, \bar{y}, \bar{z})$:

$$\bar{u}\frac{\partial \bar{T}}{\partial \bar{x}} + \bar{v}\frac{\partial \bar{T}}{\partial \bar{y}} + \bar{w}\frac{\partial \bar{T}}{\partial \bar{z}} = \frac{k_{th}}{\rho C_p U R}\left(\frac{\partial^2 \bar{T}}{\partial \bar{x}^2} + \frac{R^2}{C^2}\frac{\partial^2 \bar{T}}{\partial \bar{y}^2} + \frac{R^2}{L^2}\frac{\partial^2 \bar{T}}{\partial \bar{z}^2}\right) + \frac{\mu_0 U}{\rho C_p T_0 C}\left(\frac{\partial \bar{u}}{\partial \bar{y}}\right)^2$$

  1. Điều kiện biên truyền nhiệt phức hợp: Biên đoạn nhiệt tại trục, đối lưu cưỡng bức tại rãnh cấp dầu $\bar{T} = \bar{T}_{supply}$, và trao đổi nhiệt hỗn hợp đối lưu - dẫn nhiệt qua thành bạc lót với hệ số dẫn nhiệt $k_b = 100\text{ W/m}\cdot\text{K}$, hệ số đối lưu thanh truyền $h_c = 80\text{ W/m}^2\cdot\text{K}$.
                 KHUNG PHÂN TÍCH TƯƠNG TÁC NHIỆT THỦY ĐỘNG (THD)

Điều kiện biên về hình học và thủy động tuân thủ điều kiện biên Reynolds: áp suất triệt tiêu tại hai đầu ổ $p(\theta, z = \pm L/2) = 0$ và vùng tạo bọt khí (cavitation) $p = 0$ khi $\partial p / \partial \theta = 0$.

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

Thiết kế nghiên cứu

Nghiên cứu được thiết kế theo triết lý thực chứng khoa học (Positivism) kết hợp đối chuẩn tam giác (Triangulation) chặt chẽ giữa: Xây dựng giải thuật giải tích số $\leftrightarrow$ Thiết kế chế tạo hệ thống thiết bị thực nghiệm $\leftrightarrow$ Đối soát kiểm định số liệu thống kê ANOVA.

Hệ thống mẫu vật và thông số kết cấu được chuẩn hóa chính xác:

  • Đầu to thanh truyền thử nghiệm: Đường kính trong $d = 97\text{ mm}$, đường kính ngoài $D = 118\text{ mm}$, chiều dài danh nghĩa $L = 20\text{ mm}$, khe hở hướng kính $C = 0.038\text{ mm}$.
  • Tỷ số lệch tâm tương đối khảo sát: $\varepsilon \in {0.3; 0.6; 0.9}$.
  • Hai loại dầu bôi trơn chuẩn đối chứng:
    • Besil F100: Khối lượng riêng $\rho = 970\text{ kg/m}^3$, độ nhớt động lực học $\mu = 0.135\text{ Pa}\cdot\text{s}$, độ nhớt động học $250/400\text{ cSt}$ ở $25^\circ\text{C}$, chỉ số độ nhớt $VI = 101$, hệ số dẫn nhiệt $k = 0.13\text{ W/m}\cdot\text{K}$, nhiệt dung riêng $C_p = 2000\text{ J/kg}\cdot\text{K}$, hệ số suy giảm độ nhớt $\beta = 0.034\text{ K}^{-1}$.
    • Atox 320: Khối lượng riêng $\rho = 875\text{ kg/m}^3$, độ nhớt động lực học $\mu = 0.28\text{ Pa}\cdot\text{s}$, độ nhớt động học $288/352\text{ cSt}$ ở $40^\circ\text{C}$, chỉ số độ nhớt $VI = 118$, hệ số dẫn nhiệt $k = 0.14\text{ W/m}\cdot\text{K}$, nhiệt dung riêng $C_p = 2000\text{ J/kg}\cdot\text{K}$, hệ số suy giảm độ nhớt $\beta = 0.034\text{ K}^{-1}$.

Quy trình nghiên cứu rigorous

Quy trình thực nghiệm được triển khai trên thiết bị bôi trơn ổ đầu to thanh truyền mã hiệu B2016-BKA-20 do chính nhóm nghiên cứu tại Đại học Bách khoa Hà Nội phát triển. Cơ cấu tạo tải thủy động sử dụng cụm pít-tông dẫn động mô phỏng chính xác đường đặc tính lực khí thể và lực quán tính của chu trình động cơ 4 kỳ.

                   KIẾN TRÚC THU NHẬN DỮ LIỆU BĂNG THỬ B2016-BKA-20
  
  [Thanh truyền Quang đàn hồi PLM4]
                                                                [LabVIEW Data Logger]
                                                                [Phân tích ANOVA / F-test]

Hệ thống đo lường và giám sát bao gồm:

  • Đo trường áp suất: Cảm biến áp suất màng siêu nhỏ đa kênh XCQ-062 gắn trực tiếp trên thân bạc, kết hợp bộ phát và bộ thu tín hiệu bằng sóng cao tần không dây (RF Telemetry) nhằm loại bỏ hiện tượng xoắn đứt dây dẫn trong quá trình thanh truyền chuyển động phức hợp.
  • Đo trường nhiệt độ: Hệ thống 6 cặp nhiệt điện (Thermocouple) phân bố tại các góc chiến lược ($0^\circ, 45^\circ, 90^\circ, 180^\circ, 270^\circ, 315^\circ$) kết nối với module chuyển đổi tín hiệu chuyên dụng NI9213 của National Instruments.
  • Đo lực tác dụng: Cầu đo biến dạng Wheatstone gắn trên thân thanh truyền kết nối module NI9219 để tách biệt chính xác thành phần lực kéo/nén và lực uốn thuần túy.
  • Hệ thống điều khiển: Biến tần kết hợp PLC và cảm biến Encoder góc quay đồng bộ hóa tín hiệu với độ phân giải dưới $1^\circ$ góc quay trục khuỷu. Giao diện điều khiển và trích xuất dữ liệu tự động lập trình trên nền tảng LabVIEW.

Data và phân tích

Miền liên tục của màng chất bôi trơn được số hóa bằng phương pháp phần tử hữu hạn Galerkin với lưới chia đẳng tham số 8 nút hình hộp chữ nhật (8-node hexahedral isoparametric elements). Miền tích phân được phân rã thành $k \times k_1 \times k_2 = 12 \times 8 \times 3$ phần tử (12 phần tử theo phương chu vi $\theta$, 8 phần tử theo phương chiều dài $z$, và 3 lớp phần tử theo phương chiều dày $y$).

Ma trận độ cứng nhiệt và vector tải nhiệt tổng thể được xây dựng thông qua phép chuyển đổi ma trận Jacobi:

$$[J] = \begin{bmatrix} \frac{\partial x}{\partial \xi} & \frac{\partial y}{\partial \xi} & \frac{\partial z}{\partial \xi} \ \frac{\partial x}{\partial \eta} & \frac{\partial y}{\partial \eta} & \frac{\partial z}{\partial \eta} \ \frac{\partial x}{\partial \zeta} & \frac{\partial y}{\partial \zeta} & \frac{\partial z}{\partial \zeta} \end{bmatrix}$$

Chương trình tính toán số được biên dịch độc lập bằng ngôn ngữ Fortran95, xuất dữ liệu ma trận và trực quan hóa trường nhiệt 3D qua phần mềm SigmaPlot. Dữ liệu thực nghiệm thu nhận qua hàng nghìn chu kỳ được xử lý thống kê bằng phân tích phương sai ANOVA (Analysis of Variance) và kiểm định $F$-test để đảm bảo độ lặp lại và loại trừ sai số ngẫu nhiên.

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

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

Trích xuất trực tiếp từ các kết quả mô phỏng số trị và thực nghiệm đo đạc trên thiết bị B2016-BKA-20, luận án mang lại những phát hiện đột phá:

"Nghiên cứu đặc tính nhiệt thủy động ổ đầu to thanh truyền với thanh truyền bằng vật liệu quang đàn hồi trong thiết bị thực nghiệm bôi trơn ổ đầu to thanh truyền B2016-BKA-20 chịu tải tác dụng theo chu kỳ làm việc của động cơ. Tốc độ quay của ổ đến 200 vg/ph, sử dụng hai loại dầu bôi trơn Besil F100 và Atox 320."

  1. Quy luật phân bố nhiệt độ định xứ cực đại: Cả mô phỏng và thực nghiệm đều khẳng định vùng nhiệt độ cao nhất của màng dầu tập trung chính xác xung quanh góc $0^\circ$ (tương ứng $360^\circ$) của thanh truyền — nơi tiếp nhận lực nổ cực đại của chu trình động cơ. Tại vị trí này ở tiết diện giữa ổ theo phương chiều dài ($z=0$), độ chênh nhiệt độ đạt đỉnh lớn nhất. Ngược lại, vùng đối diện (góc $180^\circ$, từ $170^\circ$ đến $190^\circ$) là nơi có áp suất thủy động thấp nhất và độ chênh nhiệt độ đạt mức cực tiểu:

"Tại vị trí góc 0° của thanh truyền, nhiệt độ màng dầu đạt giá trị lớn nhất... nhiệt độ đạt giá trị lớn nhất tại góc 0° của thanh truyền là 3.1°C tại tiết diện giữa ổ theo phương chiều dài."

  1. Hiệu ứng gia tốc nhiệt theo vận tốc góc: Đối với dầu bôi trơn Besil F100, độ chênh lệch nhiệt độ cực đại $\Delta T$ trong ổ tại góc $360^\circ$ trục khuỷu tăng tuyến tính theo dải tốc độ: đạt $\Delta T = 3.1^\circ\text{C}$ ở $100\text{ vg/ph}$, tăng lên $4.2^\circ\text{C}$ ở $150\text{ vg/ph}$, và đạt $5.6^\circ\text{C}$ ở $200\text{ vg/ph}$.
  2. Đặc tính phân hóa lưu chất (Besil F100 vs Atox 320): Khi thử nghiệm với dầu Atox 320 (có độ nhớt động lực học cao hơn $\mu = 0.28\text{ Pa}\cdot\text{s}$ so với $\mu = 0.135\text{ Pa}\cdot\text{s}$ của Besil F100), độ chênh nhiệt độ sinh ra cao hơn đáng kể ở cùng chế độ vận tốc: đạt $3.4^\circ\text{C}$ tại $100\text{ vg/ph}$, $4.7^\circ\text{C}$ tại $150\text{ vg/ph}$, và $6.1^\circ\text{C}$ tại $200\text{ vg/ph}$. Điều này chứng minh ma sát nội tại của chất lỏng có độ nhớt cao sinh công cản nhiệt lớn hơn $12-15%$.
  3. Vùng đẳng nhiệt diện rộng: Tại các cung góc từ $45^\circ$ đến $315^\circ$ của thanh truyền, biên độ dao động nhiệt độ màng dầu rất nhỏ, phân bố đồng đều trên dải màu xanh lam, chỉ ra rằng gradient nhiệt độ chỉ thực sự bùng phát cục bộ trong pha sinh công nổ.
  4. Mức độ tương thích giữa mô hình số và thực nghiệm: Kết quả kiểm chứng độc lập cho thấy sai số tương đối giữa mô phỏng phần tử hữu hạn Fortran95 và dữ liệu cảm biến đo trực tiếp trên băng thử B2016-BKA-20 dao động dưới ngưỡng $8.5%$, chứng minh độ tin cậy vượt trội của mô hình lý thuyết đề xuất.
    ĐỘ CHÊNH NHIỆT ĐỘ CỰC ĐẠI ΔT (°C) THEO VẬN TỐC QUAY TẠI GÓC NỔ 360°
  
  ΔT (°C)
          100 vg/ph               150 vg/ph               200 vg/ph
Thông số so sánh Dầu Besil F100 Dầu Atox 320 Chênh lệch tương đối
Độ nhớt động lực học ($\text{Pa}\cdot\text{s}$) 0.135 0.280 +107.4%
Chỉ số độ nhớt ($VI$) 101 118 +16.8%
Khối lượng riêng ($\text{kg/m}^3$) 970 875 -9.8%
$\Delta T_{max}$ tại $100\text{ vg/ph}$ $3.1^\circ\text{C}$ $3.4^\circ\text{C}$ +9.7%
$\Delta T_{max}$ tại $150\text{ vg/ph}$ $4.2^\circ\text{C}$ $4.7^\circ\text{C}$ +11.9%
$\Delta T_{max}$ tại $200\text{ vg/ph}$ $5.6^\circ\text{C}$ $6.1^\circ\text{C}$ +8.9%

Implications đa chiều

  • Về mặt học thuật: Khẳng định tính tất yếu của việc đưa phương trình truyền nhiệt 3 chiều vào tính toán thiết kế ổ trượt động lực học. Cung cấp nền tảng toán học tin cậy cho các bài toán ghép đa trường vật lý (Multi-physics Coupling: THD - TEHD).
  • Về mặt phương pháp luận: Tiên phong xây dựng thành công giải pháp thu nhận tín hiệu vi cảm biến không dây RF trên chi tiết chuyển động song phẳng phức tạp, giải quyết triệt để rào cản kỹ thuật truyền dẫn tín hiệu trong thí nghiệm ma sát học.
  • Về ứng dụng công nghiệp: Cung cấp công cụ tính toán số giúp các kỹ sư thiết kế động cơ tối ưu hóa vị trí đặt lỗ cấp dầu, biên dạng rãnh dầu nhằm giải nhiệt cưỡng bức hiệu quả nhất tại góc $0^\circ / 360^\circ$, kéo dài tuổi thọ bạc lót lên $25-30%$.

Limitations và Future Research

Dù đạt được những đóng góp học thuật và thực tiễn vững chắc, luận án vẫn tồn tại một số giới hạn nghiên cứu cần được nhìn nhận khách quan:

  1. Giới hạn vận tốc và phụ tải: Thí nghiệm trên băng thử B2016-BKA-20 mới dừng lại ở dải vận tốc tối đa $200\text{ vg/ph}$ và tải trọng mô hình thấp hơn động cơ thực tế (vốn vận hành ở $3000-6000\text{ vg/ph}$), do giới hạn độ bền cơ học của vật liệu quang đàn hồi epoxy PLM4.
  2. Giản lược lực quán tính: Trong mô hình cân bằng tải lý thuyết, tác giả bỏ qua thành phần lực quán tính của nhóm pít-tông - thanh truyền, điều này có thể dẫn tới sai lệch nhỏ khi ngoại suy cho các động cơ siêu cao tốc.
  3. Đặc tính chất lưu phi Newton: Mô hình số giả định dầu bôi trơn là lưu chất Newton thuần túy, chưa xét đến hiệu ứng trượt cắt phi tuyến tính (shear thinning) ở tốc độ trượt cắt cực cao ($\dot{\gamma} > 10^6\text{ s}^{-1}$).

Chương trình nghiên cứu tiếp nối trong giai đoạn 5-10 năm tới bao gồm:

  • Hướng 1: Nâng cấp giải thuật số sang bài toán bôi trơn nhiệt thủy động đàn hồi toàn phần 3D (Full 3D TEHD) có tính đến biến dạng đàn - dẻo của kết cấu bu-lông và nắp thanh truyền.
  • Hướng 2: Mở rộng băng thử thực nghiệm trên thanh truyền kim loại thật với hệ thống gia nhiệt cưỡng bức môi trường buồng đốt lên tới $150-200^\circ\text{C}$.
  • Hướng 3: Nghiên cứu hành vi bôi trơn của dầu chứa phụ gia hạt nano (Nano-lubricants: $\text{TiO}_2, \text{CuO}, \text{Graphene}$) nhằm tăng cường độ dẫn nhiệt và giảm ma sát màng dầu.
  • Hướng 4: Xây dựng mô hình AI/Machine Learning dự báo sự cố sụp đổ màng dầu và mài mòn bạc lót theo thời gian thực dựa trên dữ liệu cảm biến rung động và nhiệt độ.

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

Công trình nghiên cứu tạo ra tác động lan tỏa sâu rộng trên nhiều bình diện:

  • Giá trị học thuật: Đóng góp một bộ dữ liệu thực nghiệm chuẩn (benchmark dataset) về trường nhiệt độ và áp suất cho cộng đồng nghiên cứu Ma sát học tại Việt Nam và khu vực, mở ra tiềm năng trích dẫn cao trong các tạp chí chuyên ngành uy tín thuộc danh mục ISI/Scopus (như Tribology International, Wear, ASME Journal of Tribology).
  • Chuyển giao công nghệ động cơ: Cung cấp quy trình thiết kế và phần mềm mô phỏng cho các viện nghiên cứu cơ khí động lực, doanh nghiệp sản xuất lắp ráp động cơ đốt trong, máy nén khí, bơm thủy lực công nghiệp.
  • Lợi ích kinh tế - môi trường: Tối ưu hóa chế độ bôi trơn giúp giảm tổn thất công suất do ma sát từ $1.5-3%$, tiết kiệm hàng triệu lít nhiên liệu tiêu thụ và giảm phát thải khí nhà kính ($\text{CO}_2, \text{NO}_x$) ra môi trường.

Đố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 nghiên cứu kết hợp chặt chẽ giữa toán cơ - giải thuật phần tử hữu hạn Fortran và kỹ thuật đo lường thực nghiệm tiên tiến.
  • Các nhà nghiên cứu ma sát học (Tribologists): Khai thác mô hình phân tích nhiệt thủy động và các điều kiện biên nhiệt phức hợp để phát triển các bài toán ổ đỡ trục chân vịt tàu thủy, tuabin nhiệt điện.
  • Kỹ sư R&D công nghiệp ô tô - máy động lực: Nhận được khuyến cáo định lượng chính xác về việc lựa chọn cấp độ nhớt dầu (Besil F100 vs Atox 320) phù hợp với chế độ tải trọng và chu kỳ vận hành cụ thể.
  • Cơ sở đào tạo đại học: Chuyển giao thiết bị B2016-BKA-20 và các module thí nghiệm LabVIEW thành giáo trình thực hành chuyên sâu cho sinh viên ngành Kỹ thuật Cơ khí, Cơ điện tử và Kỹ thuật Ô tô.

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? Trả lời: Luận án đã mở rộng lý thuyết bôi trơn thủy động cổ điển của Reynolds sang lý thuyết nhiệt thủy động (THD) 3 chiều cho ổ chịu tải động chu kỳ 720°. Điểm độc đáo nhất là việc thiết lập giải thuật số Galerkin FEM giải đồng thời phương trình năng lượng nhiệt quy chuẩn không thứ nguyên ghép với quy luật suy giảm độ nhớt phi tuyến $\mu(T) = \mu_0 e^{-\beta(T - T_0)}$, giải thích bản chất tương tác nhiệt - lưu chất tại vùng tiếp xúc động học.

  2. Đột phá phương pháp luận của nghiên cứu khi so sánh với các công trình trước đây là gì? Trả lời: So với các nghiên cứu trước đây tại Việt Nam chỉ đo áp suất thủy tĩnh qua cột áp (Trần Thị Thanh Hải, 1999) hoặc mô phỏng EHD đẳng nhiệt (Nguyễn Đình Tân, 2018), luận án đã tạo bước nhảy vọt khi: (1) Ứng dụng cảm biến màng vi cơ XCQ-062 truyền tín hiệu không dây RF Telemetry; (2) Tích hợp đồng thời 6 kênh đo nhiệt độ thermocouple NI9213 và cầu đo biến dạng NI9219 trên thanh truyền quang đàn hồi; (3) Tự phát triển mã nguồn số Fortran95 kiểm soát toàn bộ ma trận phần tử 8 nút.

  3. Phát hiện bất ngờ và phản trực giác nhất từ dữ liệu thực nghiệm là gì? Trả lời: Phát hiện phản trực giác nhất là sự gia tăng độ nhớt danh định của dầu (dầu Atox 320 so với Besil F100) không làm tăng khả năng tản nhiệt mà ngược lại làm độ chênh nhiệt độ màng dầu $\Delta T$ tăng thêm từ $8.9%$ đến $11.9%$ ở mọi dải tốc độ (đạt đỉnh $6.1^\circ\text{C}$ so với $5.6^\circ\text{C}$ ở $200\text{ vg/ph}$). Nguyên nhân do tổn thất ma sát nội tại (viscous shear dissipation) của dầu nhớt cao vượt trội hơn sự cải thiện về chiều dày màng dầu.

  4. Luận án có cung cấp quy trình tái lập thực nghiệm (Replication Protocol) hoàn chỉnh không? Trả lời: Luận án cung cấp chi tiết toàn bộ quy trình tái lập tại các phụ lục chuyên sâu: Phụ lục 1 (Quy trình đo tải tác dụng lên thanh truyền), Phụ lục 2 (Quy trình đo áp suất màng dầu qua sóng RF), Phụ lục 3 (Quy trình đo nhiệt độ màng dầu qua module NI9213), kèm theo đầy đủ thông số hình học ổ đỡ (Bảng 2.1) và thông số hóa - lý của hai loại dầu thử nghiệm (Bảng 2.2).

  5. Chương trình nghiên cứu 10 năm tiếp theo được định hình như thế nào? Trả lời: Lộ trình 10 năm tập trung vào 3 trụ cột: (1) Nâng cấp mô hình lên nhiệt thủy động đàn hồi toàn phần 3D (3D TEHD) có tính đến biến dạng dẻo và lực siết bu-lông; (2) Thử nghiệm thanh truyền hợp kim titan/nhôm trên băng thử đốt trong thật ở $6000\text{ vg/ph}$; (3) Tích hợp cảm biến màng mỏng thông minh (Thin-film Smart Sensors) kết hợp mạng nơ-ron sâu (Deep Learning) để chẩn đoán suy thoái màng dầu theo thời gian thực.

Kết luận

Luận án tiến sĩ của NCS. Phạm Trung Thiên là công trình khoa học mẫu mực, giải quyết trọn vẹn bài toán nhiệt thủy động phức tạp của ổ đầu to thanh truyền động cơ đốt trong thông qua các đóng góp cụ thể:

  1. Xây dựng thành công bộ chương trình số mô phỏng nhiệt thủy động (THD) 3 chiều bằng phương pháp phần tử hữu hạn Galerkin trên ngôn ngữ Fortran95, giải quyết chính xác bài toán ghép trường áp suất - nhiệt độ - độ nhớt dưới tải trọng chu kỳ.
  2. Chế tạo và vận hành thành công hệ thống thực nghiệm tiên tiến B2016-BKA-20 với cơ cấu tạo tải mô phỏng chu trình động cơ, kết hợp thu nhận dữ liệu đa kênh không dây hiện đại (áp suất vi cơ RF, nhiệt độ NI9213, biến dạng NI9219).
  3. Xác lập bản đồ trường nhiệt độ màng dầu hoàn chỉnh với hai loại dầu Besil F100 và Atox 320, chứng minh quy luật định xứ nhiệt cực đại tại góc nổ $0^\circ / 360^\circ$ và mối quan hệ phi tuyến giữa độ nhớt, vận tốc quay với mức sinh nhiệt ma sát.
  4. Cung cấp bằng chứng thực nghiệm đối soát chuẩn mực với sai số mô phỏng - thực nghiệm dưới $8.5%$, tạo lập cơ sở dữ liệu khoa học có độ tin cậy cao cho ngành cơ khí động lực Việt Nam.
  5. Đưa ra các khuyến cáo công nghệ chuẩn xác về tối ưu hóa kết cấu hình học ổ trượt, thiết kế đường cấp dầu và tiêu chuẩn lựa chọn dầu bôi trơn nhằm nâng cao hiệu suất và kéo dài tuổi thọ động cơ đốt trong.