‐‐‐‐ This following program solves Reynolds equation for the bearings ‐‐‐‐‐‐‐‐‐‐‐ by the numerical method to determine its quantities ‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐ Chương trình sau giải phương trình Reynold viết cho ổ đỡ ‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐‐ bằng phương pháp số để xác định các đại lượng làm việc của ổ ‐‐
Simulation Program.cpp // Họ và tên: Đặng Xuân Hải // MSSV: 20131217 // Mã lớp: 98135 // Học phần: Công nghệ bôi trơn ‐ ME4972 // Ngày thi: 28/12/2017 #include #include #include //‐‐‐‐ This following program solves Reynold's equation for the bearings ‐‐ //‐‐‐‐‐‐‐‐‐ by the numerical method to determine its quantities ‐‐‐‐‐‐‐‐‐‐‐ //‐‐‐‐‐ Chương trình sau giải phương trình Reynold viết cho ổ đỡ ‐‐‐‐‐‐‐‐‐‐ //‐‐‐‐‐‐‐‐‐ bằng phương pháp số để xác định các đại lượng làm việc của ổ ‐‐ main(void) { int i,j,x,m,n,dem=1; float c,L,exp,R; //input float a1,a2,a3,dh,theta,l,k,H; float A1,B,C,D,s,a; //các biến tính ma trận áp suất char tieptuc; //lap lai chuong trinh float coxi, muy, N1,N2,N3,N4,w=0; float W; //các biến tính khả năng tải float A[300][301]; //ma tran he so float P[50][50]; //ma tran ap suat float pg[4]; //ap suat tai cac diem Gauss float h[50][50], h1[50][50], h2[50][50], h1g[4], h2g[4]; //h matrix float hs;//h min float I1=0,I2=0; //Numerical Intergration float Ca; //frictional moment float omega, Muy; int speed; //Các thông số làm việc của ổ và dầu bôi trơn const float pi=3.141593; //Main Prog printf("‐‐Chuong trinh tinh ma tran ap suat cua o do va kha nang tai‐‐"); while(1) { printf("\n\nLan chay thu %d!",dem); dem++; printf("\nHay nhap cac thong so dau vao cua o do!"); //Input printf("\nNhap chieu dai o L(mm)= "); scanf("%f",&L); printf("Nhap ban kinh o R(mm)= "); scanf("%f", &R); printf("Nhap khe ho ban kinh C(mm)= "); scanf("%f",&c); printf("Nhap do lech tam tuong doi epsilon= "); scanf("%f",&exp); printf("Nhap so khoang chia chieu dai n= "); scanf("%d",&n); printf("Nhap so khoang chia chu vi m= "); scanf("%d", &m); c=c/1000; R=R/1000; L=L/1000; //đổi đơn vị sang mét Page 1 Simulation Program.cpp //Chapter I: Building of Presure Field //Set unchange variable k=2*pi/(m+1); l=1./(n+1); a1=1/(k*k); //Tính các hệ số của hệ for(j=1;j