Nghiên cứu các thông số công nghệ để nâng cao chất lượng và độ chính xác gia công của sản phẩm khi mài thép làm khuôn skd61

ĐẠI HỌC THÁI NGUYÊN TRƯỜNG ĐẠI HỌC KỸ THUẬT CÔNG NGHIỆP ****** THUYẾT MINH LUẬN VĂN THẠC SĨ CÔNG NGHỆ CHẾ TẠO MÁY ĐỀ TÀI: NGHIÊN CỨU CÁC THÔNG SỐ CÔNG NGHỆ ĐỂ NÂNG CAO CHẤT LƯỢNG VÀ ĐỘ CHÍNH XÁC GIA CƠNG CỦA SẢN PHẨM KHI MÀI THÉP LÀM KHN SKD61 Học Viên: Hồng Văn Quyết Lớp: CHK10 CTM Chuyên ngành: Công nghệ Chế tạo máy Người HD Khoa học: TS Nguyễn Văn Hùng NGƯỜI HƯỚNG DẪN KHOA HỌC HỌC VIÊN TS Nguyễn Văn Hùng KHOA ĐÀO TẠO SAU ĐẠI HỌC Hoàng Văn Quyết BGH TRƯỜNG ĐẠI HỌC KỸ THUẬT CƠNG NGHIỆP THÁI NGUN - 2010 Số hóa Trung tâm Học liệu – Đại học Thái Nguyên http://www.lrc-tnu.edu.vn LỜI CAM ĐOAN Việc nghiên cứu lĩnh vực công nghệ thay hàng nhập ngoại góp phần làm tăng lực cạnh tranh nhà sản xuất nước Dựa yêu cầu thực tiễn lĩnh vực khuôn mẫu, luận văn chọn hướng nghiên cứu chế tạo loại khuôn mẫu phải nhập ngoại nhằm tăng hiệu kinh tế - kỹ thuật sản phẩm Mặc dù, có cố gắng nhiều việc nghiên cứu tính tốn tiến hành thực nghiệm luận văn cịn nhiều thiếu sót cịn nhiều điểm cần đề xuất trao đổi, thảo luận thêm Tác giả mong trân trọng đóng góp, phê bình thầy giáo đồng nghiệp luận văn Em xin trân trọng cảm ơn Khoa Đào tạo Sau đại học, Ban Giám hiệu Trường Đại học Kỹ thuật Công nghiệp tạo điều kiện thuận lợi cho em suốt q trình họa tập hồn thành luận văn Xin trân trọng cảm ơn thầy giáo Hội đồng bảo vệ Đề cương luận văn góp ý, chỉnh sửa phê duyệt đề cương để luận văn em hoàn thành với nội dung tốt Đặc biệt, em xin trân trọng cảm ơn thầy TS Nguyễn Văn Hùng tận tình hướng dẫn giúp em hoàn thành luận văn Xin chân thành cảm ơn đồng nghiệp, cộng tác viên giúp đỡ đề xuất giải pháp tốt trình viết luận văn thí nghiệm Xin chân thành cảm ơn gia đình bạn bè động viên, cổ vũ tinh thần vật chất cho thân suốt trình học tập làm luận văn Em xin cam đoan số liệu kết nêu luận văn trung thực chưa công bố công trình khác Trừ phần tham khảo ghi rõ luận văn Học viên Hoàng Văn Quyết Số hóa Trung tâm Học liệu – Đại học Thái Nguyên http://www.lrc-tnu.edu.vn MỤC LỤC Phần mở đầu Tính cấp thiết đề tài Đối tượng, mục đích, nội dung phương pháp nghiên cứu Trang 1 2.1 2.2 Đối tượng nghiên cứu Mục đích nghiên cứu 2 2.3 2.4 Nội dung nghiên cứu Phương pháp nghiên cứu 2 Ý nghĩa đề tài 3.1 Ý nghĩa khoa học 3.2 Ý nghĩa thực tiến Chương 1: Tổng quan mài phẳng mài thép làm khuôn 1.1 Đặc điểm trình mài 4 1.2 1.2.1 1.2.2 1.3 1.3.1 Tổng quan mài thép làm khn Vị trí ngun công mài gia công kim loại Thực trạng ngành công nghiệp sản xuất khuôn mẫu Cơ sở trình cắt mài phẳng Quá trình tạo phoi mài phẳng 5 7 1.3.2 Động hình học trình mài phẳng 1.3.3 1.3.4 1.4 1.4.1 1.4.2 1.4.3 1.5 1.5.1 1.5.2 1.5.3 Tốc độ bóc vật liệu Chiều dày phoi tương đương Động lực học mài Lực cắt mài Công suất cắt mài Rung động mài Mòn tuổi bền đá mài Các dạng mòn đá mài Tuổi bền đá mài phương pháp xác định tuổi bền đá mài Sửa đá mài 10 11 11 11 13 13 15 15 17 19 1.6 1.7 1.7.1 1.7.2 Nhiệt độ mài Chất lượng bề mặt mài Sự hình thành nhám bề mặt mài Sự hình thành sóng bề mặt mài 19 21 21 22 Số hóa Trung tâm Học liệu – Đại học Thái Nguyên http://www.lrc-tnu.edu.vn 1.7.3 1.8 Sự thay đổi cấu trúc lớp bề mặt, hình thành ứng suất dư bề mặt Tính gia cơng vật liệu mài 22 23 1.9 Chế độ cắt ảnh hưởng chế độ cắt đến thông số đặc trưng cho trình cắt mài phẳng 24 1.9.1 Ảnh hưởng chế độ cắt đến lực cắt 24 1.9.2 1.9.3 Ảnh hưởng chế độ cắt đến rung động mài Ảnh hưởng chế độ cắt đến nhiệt cắt 25 25 1.9.4 1.9.5 Ảnh hưởng chế độ cắt đến mòn tuổi bền đá mài Ảnh hưởng chế độ cắt đến độ nhám bề mặt mài 26 26 1.10 Kết luận chương 1.10.1 Tóm tắt cơng nghệ mài, mài vật liệu làm khuôn mẫu 1.10.2 Định hướng nghiên cứu đề tài 27 27 28 1.10.3 Nội dung nghiên cứu đề tài 1.10.4 Định hướng nghiên cứu 28 28 29 Chương 2: Tối ưu hóa chế độ cắt mài vật liệu thép làm khuôn SKD61 2.1 Mô hình trình mài phẳng 2.2 Những định hướng nghiên cứu tối ưu hóa chế độ cắt mài thép làm khn SKD61 2.3 2.4 Mơ hình hóa tốn học trình mài phẳng Kết luận chương Chương 3: Xây dựng sở nghiên cứu thực nghiệm 3.1 Mục đích nghiên cứu thực nghiệm 3.1.1 Chọn thơng số vào 3.1.2 Chọn tiêu đánh giá 3.1.3 Nhiễu mài phẳng 3.2 Xây dựng quy hoạch thực nghiệm 3.2.1 Các nguyên tắc quy hoạch thực nghiệm 3.2.2 Chọn loại kế hoạch thực nghiệm mơ hình hồi quy thực nghiệm 3.2.3 Xác định miền qui hoạch 3.3 Hệ thống thiết bị thí nghiệm 3.3.1 Máy thí nghiệm 3.3.2 Vật liệu thí nghiệm Số hóa Trung tâm Học liệu – Đại học Thái Nguyên 29 30 31 34 35 35 35 37 38 38 38 40 51 52 52 53 http://www.lrc-tnu.edu.vn 3.3.3 3.3.4 Đá mài thí nghiệm Chế độ sửa đá 53 54 3.3.5 3.3.6 Dung dịch trơn nguội Hệ thống đo lường 54 54 3.4 Kết luận chương 56 Chương 4: Nghiên cứu thực nghiệm thảo luận kết 57 4.1 Các bước nghiên cứu tối hóa q trình mài phẳng thép làm khuôn SKD61 57 4.1.1 Các hàm mục tiêu mài tinh thép làm khuôn SKD61 58 4.1.2 Chọn dạng hàm hồi quy 58 4.2 61 4.2.1 Tối ưu hóa trình mài phẳng mài tinh thép làm khn SKD61 Xây dựng thuật toán quy hoạch thực nghiệm 4.2.2 4.2.3 4.2.4 Xác định hàm hồi quy Xác định hàm mục tiêu tối ưu Miền tối ưu hoá 61 64 66 4.3 Nghiên cứu hình thái, cấu trúc tính chất lý lớp bề mặt mài tinh thép làm khn SKD61 Nghiên cứu hình thái bề mặt gia cơng 71 73 75 4.4.2 4.4.3 4.5 4.5.1 4.5.2 Nghiên cứu cấu trúc tính chất lý lớp kim loại bề mặt Thảo luận kết nghiên cứu mài tinh thép làm khuôn SKD61 Ảnh hưởng chế độ cắt đến độ nhám bề mặt độ xác gia cơng Hình thái bề mặt gia cơng Cấu trúc lớp kim loại bề mặt gia công Triển khai ứng dụng kết vào sản phẩm thực tế Sản phẩm ứng dụng Kết ứng dụng 4.6 Kết luận chương 79 80 82 84 4.3.1 4.3.2 4.4 4.4.1 KẾT LUẬN CHUNG TÀI LIỆU THAM KHẢO PHỤ LỤC Số hóa Trung tâm Học liệu – Đại học Thái Nguyên 61 71 75 76 76 77 77 78 http://www.lrc-tnu.edu.vn DANH MỤC CÁC KÍ HIỆU CHÍNH Ký hiệu Ý nghĩa Đơn vị nđ Tốc độ quay đá mài Vb Vận tốc bàn máy (lượng chạy dao dọc) m/ph Sn Lượng chạy dao ngang mm/ht az Chiều sâu cắt hạt mài mm bz Chiều rộng phoi cắt mm t Chiều sâu cắt mài mm b Chiều rộng mài mm B Chiều rộng đá mài mm Dđ Đường kính đá mài mm đ Tốc độ đá mài m/s ct Tốc độ chi tiết gia công m/ph Lc Chiều dài cung tiếp xúc tĩnh mm De Đường kính tương đương mm Chiều dày phoi khơng biến dạng lớn mm htđ Chiều dày phoi tương đương mm Qw Tốc độ bóc vật liệu Q’ w Tốc độ bóc vật liệu đơn vị bề rộng mài Piz Thành phần lực cắt theo phương tiếp tuyến tác dụng lên hạt mài N Piy Thành phần lực cắt theo phương pháp tuyến tác dụng lên hạt mài N Pz Thành phần lực cắt tiếp tuyến N Py Thành phần lực cắt pháp tuyến N hmax Ra, Rz, Rt Thông số đánh giá độ nhám bề mặt gia công Kp = Py/Pz Hệ số lực cắt N Công suất mài u Năng lượng riêng mài Kc Hệ số khả cắt đá mài G Hệ số mài T Tuổi bền đá mài Tm Nhiệt độ mài Ssđ Lượng chạy dao ngang sửa đá Số hóa Trung tâm Học liệu – Đại học Thái Nguyên Vòng/ph mm3/s mm3/s.m m W J/mm3 mm3/s.N phút C m/ph http://www.lrc-tnu.edu.vn tsđ Chiều sâu cắt sửa đá xi Giá trị mã hố thơng số vào Xi g mm Khoảng thay đổi thông số vào Gia tốc Số hóa Trung tâm Học liệu – Đại học Thái Nguyên m/s2 http://www.lrc-tnu.edu.vn DANH MỤC CÁC BẢNG BIỂU TT Bảng 1.1 Nội dung Trang Hệ số truyền nhiệt vật liệu phụ thuộc vào hàm lượng 20 hợp kim 3.1 Ma trận quy hoạch trung tâm hợp thành đối xứng với thông số ảnh hưởng 52 3.2 Ký hiệu tương đương thép SKD61 nước 53 3.3 Thành phần hoá học mẫu thí nghiệm thép SKD61 53 4.1 Giá trị mã hóa giá trị thực mức thông số 60 4.2 Ma trận thí nghiệm mài tinh thép SKD61 60 4.3 Kết thí nghiệm mài tinh thép SKD61 đá 61 Cn60.CV1.G.V1 250x25x75 đ = 30m/s 4.4 Ứng dụng thông số chế độ cắt tối ưu vào mài tinh nửa khuôn dập trục khuỷu xe máy Honda Cơng ty Diezel Sơng Cơng (DISOCO) Số hóa Trung tâm Học liệu – Đại học Thái Nguyên 78 http://www.lrc-tnu.edu.vn DANH MỤC CÁC HÌNH VẼ TT Hình số Nội dung 1.1 Quá trình tách phoi hạt mài 1.2 Sơ đồ cắt hạt mài 1.3 Vùng tiếp xúc đá – chi tiết mài 10 1.4 Lực cắt tác dụng lên hạt mài 11 1.5 Quan hệ Py – Pz 13 1.6 Rung động gây sóng bề mặt gia cơng 14 1.7 Các dạng mòn đá mài 16 1.8 Q trình mịn đá mài 16 1.9 Cấu trúc lớp bề mặt mài 23 10 1.10 Profin bề mặt mài lý tưởng 26 11 2.1 Mơ hình trình mài phẳng 29 12 3.1 Sơ đồ nghiên cứu thực nghiệm 35 13 3.2 Sơ đồ khối xử lý kết đo 41 14 3.3 Sơ đồ khối xác định hàm hồi quy biến 43 15 3.4 Sơ đồ khối xác định hàm hồi quy nhiều biến 48 16 3.5 Sơ đồ khối thuật toán Gradient 51 17 3.6 Máy đo độ nhám bề mặt SJ – 201 54 18 3.7 Máy đo tọa độ CMM 544 – Mitutoyo Trang 55 R 19 3.8 Đầu đo tiếp xúc TP20 – Hãng RENISHAW 56 20 4.1 Sơ đồ cấu trúc chương trình phần mềm hồi quy tối ưu hóa 57 21 4.2a Quan hệ nhám bề mặt với Sn, Vb (t tối ưu) 67 22 4.2b Quan hệ nhám bề mặt với t, Vb (Sn tối ưu) 67 23 4.2c Quan hệ nhám bề mặt với t, Sn (Vb tối ưu) 68 24 4.3a Quan hệ sai số gia công với Sn, Vb (t tối ưu) 69 25 4.3b Quan hệ sai số gia công với t, Vb (Sn tối ưu) 69 26 4.3c Quan hệ sai số gia công với t, Sn (Vb tối ưu) 70 27 4.4 Ảnh SEM bề mặt mài thép SKD61 đá mài Hải 71 Dương Cn60.CV1.G.V1 250x25x75 với đ = 30m/s; t = 0,005 mm; Sn = mm/ht; b = m/ph; dung dịch trơn Số hóa Trung tâm Học liệu – Đại học Thái Nguyên http://www.lrc-tnu.edu.vn nguội nhũ tương hãng TOTAN pha với nước đạt nồng độ 10% 28 4.5 Ảnh SEM bề mặt mài thép SKD61 đá mài Hải 72 Dương Cn60.CV1.G.V1 250x25x75 với đ = 30m/s; t = 0,01 mm ; Sn = mm/ht; b = 11 m/ph; dung dịch trơn nguội nhũ tương hãng TOTAN pha với nước đạt nồng độ 10% 29 4.6 Ảnh SEM bề mặt mài thép SKD61 đá mài Hải 72 Dương Cn60.CV1.G.V1 250x25x75 với đ = 30m/s; t = 0,01 mm ; Sn = mm/ht; b = 11 m/ph; dung dịch trơn nguội nhũ tương hãng TOTAN pha với nước đạt nồng độ 10% 30 4.7 Cấu trúc bề mặt mài thép SKD61 đá mài Hải 73 Dương Cn60.CV1.G.V1 250x25x75 với đ = 30m/s; t = 0,005 mm ; Sn = mm/ht; b = m/ph; dung dịch trơn nguội nhũ tương hãng TOTAN pha với nước đạt nồng độ 10% 31 4.8 Cấu trúc bề mặt mài thép SKD61 đá mài Hải 74 Dương Cn60.CV1.G.V1 250x25x75 với đ = 30m/s; t = 0,01 mm ; Sn = mm/ht; b = 11 m/ph; dung dịch trơn nguội nhũ tương hãng TOTAN pha với nước đạt nồng độ 10% 32 4.9 Cấu trúc bề mặt mài thép SKD61 đá mài Hải 74 Dương Cn60.CV1.G.V1 250x25x75 với đ = 30m/s; t = 0,01 mm ; Sn = mm/ht; b = 11 m/ph; dung dịch trơn nguội nhũ tương hãng TOTAN pha với nước đạt nồng độ 10% 33 4.10 Nửa khuôn dập trục khuỷu xe máy Honda 77 34 4.11 Nửa khn dập trục khuỷu xe máy Honda 78 Số hóa Trung tâm Học liệu – Đại học Thái Nguyên http://www.lrc-tnu.edu.vn -91- % xp0=[0.0075 8]; LB=[0.0045 4.6 4.4]; UB=[0.011 9.4 11.6]; % defunc=inline('exp(kp1+kp2*log(xp(1))+log(xp(2))*kp3+log(xp(3))*kp4+kp5*l og(xp(1))*log(xp(2))+kp6*log(xp(1))*log(xp(3))+kp7*log(xp(2))*log(xp(3))+ kp8*log(xp(1))^2+kp9*log(xp(2))^2+kpA*log(xp(3))^2)','xp','kp1','kp2','kp 3','kp4','kp5','kp6','kp7','kp8','kp9','kpA'); options=optimset('LargeScale','off'); [deoptim,deval,deexit,deoutput,delamda,degrad]=fmincon(defunc,xp0,[],[],[ ],[],LB,UB,[],options,bde(1),bde(2),bde(3),bde(4),bde(5),bde(6),bde(7),bd e(8),bde(9),bde(10)); %================== fprintf(fp3,'\n\nSAI SO GIA CONG DELTA:'); fprintf(fp3,'\n '); fprintf(fp3,'\n\n SAU KHI MU HOA- Sai so gia cong co quan he theo dang:'); fprintf(fp3,'\nDe= '); in(fp3,br); fprintf(fp3,'\n '); save tinhcx.mat Ttinh Stinh Vtinh bde deoptim; if rexit > 0.0 fprintf(fp3,'\n Ham muc tieu hoi tu'); fprintf(fp3,'\n Sai so gia cong Delta dat tai:'); fprintf(fp3,'\nt=%12.4f\nS=%12.4f\nV=%12.3f\n',deoptim); fprintf(fp3,'\n Delta min=%f\n\n',deval); fprintf(fp3,'\nSo buoc lap :%d',deoutput.iterations); fprintf(fp3,'\nThuat toan :%s',deoutput.algorithm); fprintf(fp3,'\nSo lan tinh toan gia tri ham muc tieu:%d',deoutput.funcCount); fprintf(fp3,'\nFinal step size:%d',deoutput.stepsize); fprintf(fp3,'\nFisrt order optimality:%d',deoutput.firstorderopt); fprintf(fp3,'\n Gia tri GRADIENT cua ham muc tieu :%f',degrad); else fprintf(fp3,'\nHam muc tieu khong hoi tu'); end Số hóa Trung tâm Học liệu – Đại học Thái Nguyên http://www.lrc-tnu.edu.vn -92- fprintf(fp3,'\n*********************************************************' ); %%%%%%%%%%%%%% stat1=fclose(fp2); stat2=fclose(fp3); msgbox('Chuong trinh ket thuc tot dep!','MatLab thong bao'); clc; disp('ket qua hoi qui ghi vao file HQB2TINH.OUT'); disp('Ket qua toi uu ghi vao file OPGCTINH.OUT'); %************************************** function [x]= vec2matr(m,n,A) %generate matrix(m,n) from vector A syntax [x]=vec2ma(m,n,A) %If size of matrix X is bigger than size of Vector A then the remain will be filled with x=zeros(m,n); z=m*n; l=length(A); if z >l for t=l+1:z A(t)=0; end end % for j=1:m for k=1:n i=(j-1)*n+k; x(j,k)=A(i); end end %%%%%%%%%%%%%%%%%%%%% function in(fp,bp) fprintf(fp,' e^(%f',bp(1)); if bp(2)>0 fprintf(fp,'+%f*lnt',bp(2)); else fprintf(fp,'%f*lnt',bp(2)); end if bp(3)>0 fprintf(fp,'+%flnS',bp(3)); Số hóa Trung tâm Học liệu – Đại học Thái Nguyên http://www.lrc-tnu.edu.vn -93- else fprintf(fp,'%flnS',bp(3)); end if bp(4)>0 fprintf(fp,'+%flnV',bp(4)); else fprintf(fp,'%flnV',bp(4)); end if bp(5)>0 fprintf(fp,'+%flnt*lnS',bp(5)); else fprintf(fp,'%flnt*lnS',bp(5)); end if bp(6)>0 fprintf(fp,'+%flnt*lnV',bp(6)); else fprintf(fp,'%flnt*lnV',bp(6)); end if bp(7)>0 fprintf(fp,'+%flnS*lnV',bp(7)); else fprintf(fp,'%flnS*lnV',bp(7)); end if bp(8)>0 fprintf(fp,'+%flnt^2',bp(8)); else fprintf(fp,'%flnt^2',bp(8)); end if bp(9)>0 fprintf(fp,'+%flnS^2',bp(9)); else fprintf(fp,'%flnS^2',bp(9)); end if bp(10)>0 fprintf(fp,'+%flnV^2)',bp(10)); else fprintf(fp,'%flnV^2)',bp(10)); end Số hóa Trung tâm Học liệu – Đại học Thái Nguyên http://www.lrc-tnu.edu.vn -94- Số liệu thí nghiệm So lan thi nghiem mai tinh thep SKD61 15 So nhan to So thi nghiem o tam t(mm) S(mm/ht) V(m/ph) Ra(µm) Rz(µm) Delta(µm) 0.005 0.23 1.44 15 0.010 5 0.30 1.78 8.4 0.005 0.19 1.16 17.6 0.010 0.31 1.82 8.4 0.005 11 0.22 1.59 7.4 0.010 11 0.56 3.55 10.5 0.005 11 0.31 2.04 9.6 0.010 11 0.39 1.92 5.2 0.011 0.46 3.18 7.6 0.0045 0.18 1.24 9.5 0.0075 9.4 0.37 2.25 4.9 0.0075 4.6 0.34 1.85 7.2 0.0075 11.6 0.45 2.74 9.8 0.0075 4.4 0.20 1.23 11.5 0.0075 0.26 1.44 10.2 Thi nghiem lap Ra0 Rz0 De0 0.25 1.52 10.5 0.26 1.44 10.3 0.23 1.42 10.2 Số hóa Trung tâm Học liệu – Đại học Thái Nguyên http://www.lrc-tnu.edu.vn -95- Hàm đồ thị nhám bề mặt function dothinham %Ve thi cua %Cac gia tri can de ve thi duoc load tu file %co dinh mot bien de ve thi sau ve contour clc load tinhra [x,y]=meshgrid([min(Ttinh):0.002:max(Ttinh)],[min(Stinh):0.2:max(Stinh)]) ; z=exp(br(1)+br(2)*log(x)+br(3)*log(y)+br(4)*log(roptim(3))+br(5)*log(x).* log(y)+br(6)*log(roptim(3))*log(x)+br(7)*log(roptim(3))*log(y)+br(8)*log( x).^2+br(9)*log(y).^2+br(10)*log(roptim(3))^2); h1=figure('Position' , [120 90 400 460]); subplot(2,1,1); surfc(x,y,z) xlabel('Chieu sau cat t(mm)') ylabel('Luong chay dao S (mm/ht)') zlabel('Nham be mat Ra (\mu m)') set(gca,'GridlineStyle','-') grid on title('Quan he cua Nham be mat Ra voi t,S ({\itV toi uu})','FontSize',14) subplot(2,1,2) [C,h]=contour(x,y,z,5); Clabel(C,h); xlabel('Chieu sau cat t(mm)') ylabel('Luong chay dao S (mm/ht)') set(gca,'GridlineStyle','-') grid on print -djpeg -r200 ra_t_s % h2=figure('Position' , [120 90 400 460]); subplot(2,1,1); [x,y]=meshgrid([min(Stinh):0.2:max(Stinh)],[min(Vtinh):1:max(Vtinh)]); z=exp(br(1)+br(2)*log(roptim(1))+br(3)*log(x)+br(4)*log(y)+br(5)*log(ropt im(1))*log(x)+br(6)*log(roptim(1))*log(y)+br(7)*log(x).*log(y)+br(8)*log( roptim(1))^2+br(9)*log(x).^2+br(10)*log(y).^2); surfc(x,y,z) xlabel('Luong chay dao S(mm/ht)') Số hóa Trung tâm Học liệu – Đại học Thái Nguyên http://www.lrc-tnu.edu.vn -96- ylabel('Van toc ban V (m/ph)') zlabel('Nham be mat Ra (\mu m)') set(gca,'GridlineStyle','-') grid on title('Quan he cua Nhan be mat Ra voi S,V ({\itt toi uu})','FontSize',14) subplot(2,1,2) [C,h]=contour(x,y,z,5); Clabel(C,h); xlabel('Luong chay dao S (mm/ht)') ylabel('Van toc ban V (m/ph)') set(gca,'GridlineStyle','-') grid on print -djpeg ra_s_v %%%%%5 h3=figure('Position' , [120 90 400 460]); subplot(2,1,1); [x,y]=meshgrid([min(Ttinh):0.002:max(Ttinh)],[min(Vtinh):1:max(Vtinh)]); z=exp(br(1)+br(2)*log(x)+br(3)*log(roptim(2))+br(4)*log(y)+br(5)*log(ropt im(2))*log(x)+br(6)*log(x).*log(y)+br(7)*log(roptim(2))*log(y)+br(8)*log( x).^2+br(9)*log(roptim(2))^2+br(10)*log(y).^2); surfc(x,y,z) xlabel('Chieu sau cat t(mm)') ylabel('Van toc ban V (m/ph)') zlabel('Nham be mat Ra (\mu m)') set(gca,'GridlineStyle','-') grid on title('Quan he cua Nham be mat Ra voi t,V ({\itS toi uu})','FontSize',14) subplot(2,1,2) [C,h]=contour(x,y,z,5); Clabel(C,h); xlabel('Chieu sau cat t (mm)') ylabel('Van toc ban V (m/ph)') set(gca,'GridlineStyle','-') grid on print -djpeg ra_t_v Số hóa Trung tâm Học liệu – Đại học Thái Nguyên http://www.lrc-tnu.edu.vn -97- Hàm đồ thị sai số gia công function dtsaisogiacong %Ve thi cua chinh xac Delta %Cac gia tri can de ve thi duoc load tu file %co dinh mot bien de ve thi sau ve contour clc load tinhcx [x,y]=meshgrid([min(Ttinh):0.002:max(Ttinh)],[min(Stinh):0.2:max(Stinh)]) ; z=exp(bde(1)+bde(2)*log(x)+bde(3)*log(y)+bde(4)*log(deoptim(3))+bde(5)*lo g(x).*log(y)+bde(6)*log(deoptim(3))*log(x)+bde(7)*log(deoptim(3))*log(y)+ bde(8)*log(x).^2+bde(9)*log(y).^2+bde(10)*log(deoptim(3))^2); h1=figure('Position' , [120 90 400 460]); subplot(2,1,1); surfc(x,y,z) xlabel('Chieu sau cat t (mm)') ylabel('Luong chay dao S (mm/ht)') zlabel('Sai so Delta (\mu m)') set(gca,'GridlineStyle','-') grid on title('Quan he cua Sai so gia cong voi t,S ({\itV toi uu})','FontSize',14) subplot(2,1,2) [C,h]=contour(x,y,z,5); Clabel(C,h); xlabel('Chieu sau cat t (mm)') ylabel('Luong chay dao S (mm/ht)') set(gca,'GridlineStyle','-') grid on print -djpeg -r200 saiso_t_s % h2=figure('Position' , [120 90 400 460]); subplot(2,1,1); [x,y]=meshgrid([min(Stinh):0.2:max(Stinh)],[min(Vtinh):1:max(Vtinh)]); z=exp(bde(1)+bde(2)*log(deoptim(1))+bde(3)*log(x)+bde(4)*log(y)+bde(5)*lo g(deoptim(1))*log(x)+bde(6)*log(deoptim(1))*log(y)+bde(7)*log(x).*log(y)+ bde(8)*log(deoptim(1))^2+bde(9)*log(x).^2+bde(10)*log(y).^2); surfc(x,y,z) Số hóa Trung tâm Học liệu – Đại học Thái Nguyên http://www.lrc-tnu.edu.vn -98- xlabel('Luong chay dao S (mm/ht)') ylabel('Van toc ban V (m/ph )') zlabel('Sai so Delta (\mu m)') set(gca,'GridlineStyle','-') grid on title('Quan he cua Sai so gia cong voi S,V ({\itt toi uu})','FontSize',14) subplot(2,1,2) [C,h]=contour(x,y,z,5); Clabel(C,h); xlabel('Luong chay dao S (mm/ht)') ylabel('Van toc ban V (m/ph )') set(gca,'GridlineStyle','-') grid on print -djpeg saiso_s_v %%%%%5 h3=figure('Position' , [120 90 400 460]); subplot(2,1,1); [x,y]=meshgrid([min(Ttinh):0.002:max(Ttinh)],[min(Vtinh):1:max(Vtinh)]); z=exp(bde(1)+bde(2)*log(x)+bde(3)*log(deoptim(2))+bde(4)*log(y)+bde(5)*lo g(deoptim(2))*log(x)+bde(6)*log(x).*log(y)+bde(7)*log(deoptim(2))*log(y)+ bde(8)*log(x).^2+bde(9)*log(deoptim(2))^2+bde(10)*log(y).^2); surfc(x,y,z) xlabel('Chieu sau cat t (mm)') ylabel('Van toc ban V (m/ph )') zlabel('Sai so Delta (\mu m)') set(gca,'GridlineStyle','-') grid on title('Quan he cua Sai so gia cong voi t,V ({\itS toi uu})','FontSize',14) subplot(2,1,2) [C,h]=contour(x,y,z,5); Clabel(C,h); xlabel('Chieu sau cat t (mm)') ylabel('Van toc ban V (m/ph )') set(gca,'GridlineStyle','-') grid on print -djpeg saiso_t_v Số hóa Trung tâm Học liệu – Đại học Thái Nguyên http://www.lrc-tnu.edu.vn -99- Kết quy hoạch thực nghiệm 5.1 Hàm hồi quy bậc hai Mai SKD61 voi che gia cong tinh Thoi gian chay chuong trinh:20-Jun-2010 *************************************************************** Ket qua hoi qui voi tin cay 95 phan tram la: Phuong trinh Hoi qui tuyen tinh cap II co dang: y=b0+b1*x1+b2*x2+b3*x3+b12*x1*x2+b13*x1*x3+b23*x2*x3+b11*x1^2+b22*x2^2+b3 3*x3^2 CAC HE SO HOI QUI CUA NHAM BE MAT Ra: b0=4.791821 b1=0.045477 b2=-7.001463 b3=2.141319 b12=-0.510209 b13=0.472273 b23=0.231714 b11=-0.078076 b22=1.051174 b33=0.068081 He so xac dinh R^2= 0.859529 Gia tri ham so thong ke de kiem dinh gia thiet F0=3.399384 va gia tri p-value=0.095605 Kiem dinh ket qua CAC SO THONG KE Phan bo Student voi muc y nghia 0.05 bac tu t_anpha=2.919986 Phan bo F voi muc y nghia 0.01 bac tu va F_anpha=99.415852 Phuong sai du:0.021562 Phuong sai tai sinh (binh phuong)=0.003920 t(1)=296.403886 > t_anpha=2.919986, He so beta(1) co nghia t(2)=13.902202 > t_anpha=2.919986, He so beta(2) co nghia t(3)=836.460008 > t_anpha=2.919986, He so beta(3) co nghia t(4)=271.845192 > t_anpha=2.919986, He so beta(4) co nghia t(5)=301.236364 > t_anpha=2.919986, He so beta(5) co nghia t(6)=296.301305 > t_anpha=2.919986, He so beta(6) co nghia Số hóa Trung tâm Học liệu – Đại học Thái Nguyên http://www.lrc-tnu.edu.vn -100- t(7)=56.815487 > t_anpha=2.919986, He so beta(7) co nghia t(8)=118.842995 > t_anpha=2.919986, He so beta(8) co nghia t(9)=250.363513 > t_anpha=2.919986, He so beta(9) co nghia t(10)=18.627338 > t_anpha=2.919986, He so beta(10) co nghia -KIEM TRA y^: -F^=S^2du/S^2ts =3.399384 t_anpha=2.919986, He so beta(1) co nghia t(2)=3491.356216 > t_anpha=2.919986, He so beta(2) co nghia t(3)=824.975281 > t_anpha=2.919986, He so beta(3) co nghia Số hóa Trung tâm Học liệu – Đại học Thái Nguyên http://www.lrc-tnu.edu.vn -101- t(4)=882.648325 > t_anpha=2.919986, He so beta(4) co nghia t(5)=3586.137076 > t_anpha=2.919986, He so beta(5) co nghia t(6)=2708.152222 > t_anpha=2.919986, He so beta(6) co nghia t(7)=714.802343 > t_anpha=2.919986, He so beta(7) co nghia t(8)=1633.872180 > t_anpha=2.919986, He so beta(8) co nghia t(9)=2004.937716 > t_anpha=2.919986, He so beta(9) co nghia t(10)=1244.742311 > t_anpha=2.919986, He so beta(10) co nghia -KIEM TRA y^: -F^=S^2du/S^2ts =2.781914