Đang tải... (xem toàn văn)
‐‐‐‐ 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