1. Thiết lập mô hình toán học
1.1 Rời rạc hóa phương trình Reynolds
Phương trình cốt lõi trong phân tích bôi trơn ổ trượt hai chiều được biểu diễn như sau:
Trong đó:
- \(h = c - e\cosθ - δ(x,z)\) là chiều dày màng dầu
- \(c\): khe hở bán kính của ổ
- \(e\): độ lệch tâm
- \(δ(x,z)): biến dạng đàn hồi của cổ trục
Phương pháp sai phân hữu hạn được áp dụng để rời rạc hóa phương trình điều khiển, sử dụng phần tử đẳng tham số bậc hai với 9 nút nhằm nâng cao độ chính xác.
1.2 Xử lý điều kiện biên
- Điều kiện biên áp suất: tại vị trí rãnh cấp dầu, áp suất bằng áp suất cấp dầu; hai bên của segment chịu tải có áp suất bằng 0 (cửa xả dầu).
- Hiệu chỉnh hiện tượng co rỗng: áp dụng điều kiện biên Reynolds, khi áp suất giảm xuống dưới áp suất môi trường, màng dầu sẽ bị phá vỡ.
- Hiệu chỉnh độ nghiêng của cổ trục: chiều dày màng dầu được bổ sung tham số góc nghiêng:
Trong đó \(α\) là góc nghiêng của cổ trục, \(β\) là góc vị trí của segment chịu tải.
2. Quy trình thực hiện phương pháp sai phân hữu hạn
2.1 Phân chia lưới và khởi tạo
- Lưới theo phương chu vi: sử dụng lưới không đều, vùng có độ lệch tâm cao được làm mịn hơn (ví dụ N = 101 nút).
- Lưới theo phương dọc trục: phân chia đều (ví dụ M = 41 nút), phương chiều dày sử dụng tọa độ logarithm để tăng cường hội tụ.
- Giá trị khởi đoán: trường áp suất được khởi tạo với phân phối đều, chiều dày màng dầu tính theo giả thiết ổ cứng.
2.2 Giải tìm trường áp suất bằng phương pháp lặp
% Xây dựng ma trận hệ số phương trình rời rạc
for i = 2:N-1
for j = 2:M-1
C1 = 0.5*(H(i+1,j)+H(i,j))^3;
C2 = 0.5*(H(i,j+1)+H(i,j))^3;
C3 = ALPHA*gamma*(0.5*(H(i,j+1)+H(i,j))^3);
C4 = ALPHA*gamma*(0.5*(H(i,j-1)+H(i,j))^3);
P(i,j) = (-DX*(H(i+1,j)-H(i,j))/2 + C1*P(i,j+1) + C2*P(i,j-1) + ...
C3*P(i+1,j) + C4*P(i-1,j)) / (C1 + C2 + C3 + C4);
end
end
% Hiệu chỉnh áp suất và kiểm tra hội tụ
P = 0.6*P_cu + 0.4*P_moi;
err = max(abs(P - P_cu));
Phương pháp lặp Newton-Raphson được sử dụng với ngưỡng hội tụ là 10⁻⁶.
2.3 Tính toán chiều dày màng dầu
Có tính đến sự kết hợp giữa nghiêng cổ trục và biến dạng đàn hồi:
h(x,z) = c - e*cosθ - δ_dh(x,z) - δ_nhiet(x,z)
- Biến dạng đàn hồi: được xác định thông qua phương pháp phần tử hữu hạn bằng cách giải ma trận độ cứng của segment.
- Biến dạng nhiệt: tính toán trường nhiệt độ dựa trên phương trình năng lượng, sử dụng mô hình đoạn nhiệt để đơn giản hóa.
3. Phân tích ảnh hưởng của các tham số quan trọng
3.1 Độ lệch tâm và tốc độ quay
- Phân phối áp suất: khi độ lệch tâm tăng, đỉnh áp suất dịch chuyển về phía vùng vào, mối quan hệ giữa áp suất cực đại và độ lệch tâm mang tính phi tuyến (xem hình 1).
- Chiều dày màng dầu: ở tốc độ cao, chiều dày màng dầu tăng lên, nhưng tại một số vùng cục bộ có thể giảm do hiệu ứng ly tâm.
3.2 Góc nghiêng của cổ trục
- Tập trung áp suất: khi góc nghiêng α > 0.01, đỉnh áp suất tăng đáng kể (ví dụ α = 0.02, áp suất tăng 15%).
- Khả năng chịu tải: tối ưu hóa góc nghiêng có thể nâng cao khả năng chịu tải, nhưng cần tránh hiện tượng phá vỡ màng dầu.
3.3 Đặc tính của chất bôi trơn
- Ảnh hưởng của nhiệt độ lên độ nhớt: sử dụng mô hình Walther để mô tả sự thay đổi độ nhớt:
Sự giảm độ nhớt dẫn đến khả năng chịu tải của màng dầu giảm, cần hiệu chỉnh lặp trường nhiệt độ.
4. Kiểm chứng số và tối ưu hóa
4.1 Các trường hợp kiểm chứng
- So sánh với thực nghiệm: với một ổ trục có ε = 0.8, áp suất cực đại tính toán được là 33.07 MPa, sai số so với giá trị trong tài liệu 33.06 MPa nhỏ hơn 0.3%.
- Kiểm chứng với độ nghiêng cổ trục: khi góc nghiêng α = 0.02, phân phối áp suất tính toán khớp với kết quả trong bài báo (sai số 1.2%).
4.2 Tối ưu hóa hội tụ
- Làm mịn lưới: vùng có gradient áp suất lớn (như vùng vào) sử dụng lưới thích ứng (xem hình 2).
- Hệ số thư giãn: trong quá trình lặp áp suất, áp dụng hệ số thư giãn động ω = 0.5 - 0.1·ln(err) để đẩy nhanh hội tụ.
4.3 Ứng dụng trong thực tế
- Chẩn đoán sự cố: phát hiện mài mòn ổ trục thông qua bất thường trong phân phối áp suất (ví dụ đỉnh áp suất lệch trên 5% báo hiệu mài mòn không đều).
- Tối ưu hóa thiết kế: khi tỷ lệ chiều rộng/đường kính B/D > 1, tính đồng đều của phân phối áp suất được cải thiện 20%.
5. Cài đặt với MATLAB (phần cốt lõi)
%% Định nghĩa tham số
L = 0.066; % Chiều rộng ổ (m)
R = 0.030; % Bán kính ổ (m)
Cgap = 0.003; % Khe hở bán kính (m)
n = 3000; % Tốc độ quay (vòng/phút)
V = n*2*pi/60; % Vận tốc dài (m/s)
mu = 0.009; % Độ nhớt động lực (Pa·s)
%% Phân chia lưới
N = 101; % Số nút theo chu vi
M = 41; % Số nút theo dọc trục
DX = 2*pi/(N-1);
DY = L/(M-1);
%% Điều kiện ban đầu
H = zeros(N,M);
P = zeros(N,M);
for i = 1:N
for j = 1:M
H(i,j) = 1 + Cgap/R * cos(2*pi*(i-1)/N) + (j-1)*DY*(tan(0.02*pi/180))/R;
end
end
%% Giải tìm áp suất bằng phương pháp lặp
so_lan_lap_toi_da = 1e3;
nguong_hoi_tu = 1e-6;
for lap = 1:so_lan_lap_toi_da
P_cu = P;
for i = 2:N-1
for j = 2:M-1
% Hệ số phương trình rời rạc
C1 = 0.5*(H(i+1,j)+H(i,j))^3;
C2 = 0.5*(H(i,j+1)+H(i,j))^3;
C3 = mu*0.009*(0.5*(H(i,j+1)+H(i,j))^3);
C4 = mu*0.009*(0.5*(H(i,j-1)+H(i,j))^3);
P(i,j) = (-DX*(H(i+1,j)-H(i,j))/2 + C1*P(i,j+1) + C2*P(i,j-1) + ...
C3*P(i+1,j) + C4*P(i-1,j)) / (C1 + C2 + C3 + C4);
end
end
% Xử lý điều kiện biên
P(:,1) = 0; P(:,M) = 0; P(1,:) = 0; P(N,:) = 0;
% Hiệu chỉnh bằng phương pháp lặp hỗn hợp
P = 0.6*P_cu + 0.4*P;
% Kiểm tra hội tụ
sai_so = max(max(abs(P - P_cu)));
if sai_so < nguong_hoi_tu, break; end
end
%% Tính chiều dày màng dầu
h_dau = H - P/(6*mu*V/R);
6. Trực quan hóa và phân tích kết quả
- Bản đồ phân phối áp suất: hiển thị dạng phân phối áp suất, xác định vùng áp cao và vùng xuất hiện hiện tượng co rỗng.
- Mặt cắt chiều dày màng dầu: phân tích vị trí chiều dày màng dầu nhỏ nhất, đánh giá an toàn bôi trơn.
- Đồ thị độ nhạy của tham số: vẽ đường cong ảnh hưởng của độ lệch tâm, tốc độ quay lên khả năng chịu tải.
Phương pháp này cho phép mô phỏng chính xác đặc tính bôi trơn của ổ trượt trong các điều kiện làm việc phức tạp, cung cấp cơ sở lý thuyết cho việc thiết kế độ tin cậy của máy móc quay.