Luận án Tiến sĩ Nghiên cứu ảnh hưởng của các thông số công nghệ và nhiệt độ đến biến dạng tạo hình khi dập vuốt chi tiết dạng cốc từ vật liệu SPCC

BỘ GIÁO DỤC VÀ ĐÀO TẠO TRƯỜNG ĐẠI HỌC SƯ PHẠM KỸ THUẬT HƯNG YÊN LUYỆN THẾ THẠNH NGHIÊN CỨU ẢNH HƯỞNG CỦA CÁC THÔNG SỐ CÔNG NGHỆ VÀ NHIỆT ĐỘ ĐẾN BIẾN DẠNG TẠO HÌNH KHI DẬP VUỐT CHI TIẾT DẠNG CỐC TỪ VẬT LIỆU SPCC LUẬN ÁN TIẾN SĨ KỸ THUẬT CƠ KHÍ Hưng Yên - 2022 BỘ GIÁO DỤC VÀ ĐÀO TẠO TRƯỜNG ĐẠI HỌC SƯ PHẠM KỸ THUẬT HƯNG YÊN LUYỆN THẾ THẠNH NGHIÊN CỨU ẢNH HƯỞNG CỦA CÁC THÔNG SỐ CÔNG NGHỆ VÀ NHIỆT ĐỘ ĐẾN BIẾN DẠNG TẠO HÌNH KHI DẬP VUỐT CHI TIẾT DẠNG CỐC TỪ VẬT LIỆU SPCC Chuyên ngành: Kỹ thuật khí Mã số: 9520103 LUẬN ÁN TIẾN SĨ KỸ THUẬT CƠ KHÍ NGƯỜI HƯỚNG DẪN KHOA HỌC: GS.TS Nguyễn Đức Toàn GS.TSKH Bành Tiến Long Hưng Yên - 2022 LỜI CÁM ƠN Tôi xin gửi lời cám ơn chân thành sâu sắc tới GS TSKH NGND Bành Tiến Long GS.TS Nguyễn Đức Toàn, người Thầy tận tình hướng dẫn, động viên, giúp đỡ tơi vượt qua khó khăn để hồn thành luận án Tôi xin chân thành cám ơn Công ty TNHH Cơ khí Ơ tơ Đức Hịa tạo điều kiện thuận lợi giúp đỡ tơi hồn thành thực nghiệm luận án Tôi xin gửi lời cám ơn tới Ban lãnh đạo Trường, Khoa Cơ khí, Bộ mơn Tự động hóa thiết kế cơng nghệ khí, Trường Đại học Sư phạm Kỹ thuật Hưng Yên tạo điều kiện giúp đỡ tơi hồn thành luận án Cuối xin gửi lời cám ơn chân thành đến tất đồng nghiệp, bạn bè, gia đình người thân ln bên động viên khích lệ mong muốn tơi hồn thành luận án Nghiên cứu sinh (đã ký) Luyện Thế Thạnh i LỜI CAM ĐOAN Tôi xin cam đoan nội dung luận án cơng trình nghiên cứu riêng tơi Tất số liệu kết nghiên cứu luận án trung thực chưa công bố cơng trình nghiên cứu khác TM TẬP THỂ HƯỚNG DẪN KHOA HỌC (đã ký) NGHIÊN CỨU SINH (đã ký) GS.TS Nguyễn Đức Toàn ii Luyện Thế Thạnh MỤC LỤC MỤC LỤC iii DANH MỤC CÁC KÝ HIỆU VÀ CHỮ VIẾT TẮT vi DANH MỤC CÁC BẢNG ix DANH MỤC CÁC HÌNH ẢNH, ĐỒ THỊ xii MỞ ĐẦU 1 Tính cấp thiết đề tài nghiên cứu Mục đích, đối tượng phạm vi nghiên cứu 2.1 Mục đích nghiên cứu 2.2 Đối tượng phạm vi nghiên cứu Phương pháp nghiên cứu Ý nghĩa khoa học ý nghĩa thực tiễn đề tài Những đóng góp đề tài Cấu trúc nội dung luận án .4 CHƯƠNG TỔNG QUAN VỀ CÔNG NGHỆ DẬP VUỐT CÓ GIA NHIỆT 1.1 Cơng nghệ dập sản xuất khí 1.2 Công nghệ dập vuốt yếu tố ảnh hưởng đến chất lượng sản phẩm dập vuốt……………………………………………………………………………… ……6 1.2.1 Khái quát công nghệ dập vuốt………………………………………………6 1.2.2 Xác định thông số trình dập vuốt…………… ….………………10 1.2.3 Các yếu tố ảnh hưởng đến chất lượng sản phẩm dập vuốt 17 1.3 Các phương pháp gia nhiệt dập vuốt 20 1.3.1 Các mơ hình gia nhiệt phơi dập vuốt……………………………………20 1.3.2 Mơ hình gia nhiệt khn dập vuốt………………………………………24 1.4 Tổng quan tình hình nghiên cứu ngồi nước gia cơng tạo hình vật liệu có gia nhiệt 25 1.4.1 Tình hình nghiên cứu nước………………………………………… 25 1.4.2 Tình hình nghiên cứu nước 26 KẾT LUẬN CHƯƠNG 34 CHƯƠNG CƠ SỞ LÝ THUYẾT VỀ PHÁ HỦY VẬT LIỆU VÀ XÂY DỰNG ĐƯỜNG CONG GIỚI HẠN TẠO HÌNH CỦA VẬT LIỆU SPCC TẠI CÁC NHIỆT ĐỘ KHÁC NHAU 35 2.1 Mơ hình thuộc tính phá huỷ vật liệu 35 iii 2.1.1 Phá huỷ dẻo vật liệu………………………………………………………… 35 2.1.2 Mơ hình phá hủy vật liệu…………………………………………………… 37 2.1.3 Mơ hình thuộc tính vật liệu 44 2.1.4 Xác định tính vật liệu SPCC sử dụng nghiên cứu 45 2.2 Xây dựng FLC vật liệu SPCC nhiệt độ phòng nhiệt độ gia nhiệt độ khác .49 2.2.1 Xây dựng FLC vật liệu SPCC nhiệt độ phòng 50 2.2.2 Xây dựng FLC vật liệu SPCC nhiệt độ gia nhiệt độ khác nhau………………………………………………………………………… …….55 KẾT LUẬN CHƯƠNG 63 CHƯƠNG NGHIÊN CỨU MÔ PHỎNG VÀ THỰC NGHIỆM ẢNH HƯỞNG CỦA MỘT SỐ THƠNG SỐ ĐẾN CHIỀU CAO TẠO HÌNH VÀ CHIỀU DÀY PHÂN BỐ CỦA CHI TIẾT DẠNG CỐC TRỤ KHI DẬP VUỐT 64 3.1 Sơ đồ hệ thống thực nghiệm .64 3.2 Thiết lập thực nghiệm 66 3.2.1 Vật liệu thực nghiệm………………………………………………………….66 3.2.2 Phôi dụng cụ phục vụ thực nghiệm…………………………………66 3.2.3 Máy thực nghiệm…………………………………………………………… 67 3.2.4 Bộ khuôn thực nghiệm……………………………………………………… 68 3.2.5 Thiết bị gia nhiệt thu thập liệu nhiệt độ………………………….69 3.2.6 Thiết bị đo .71 3.3 Mô số gia công dập vuốt 73 3.3.1 Mơ hình phần tử hữu hạn (FEM)…………………………………………… 73 3.3.2 Thiết lập thông số mô thực nghiệm……………………………74 3.4 Nghiên cứu trình gia nhiệt dập vuốt thực nghiệm .76 3.4.1 Sơ đồ thực nghiệm……………………………………………………………76 3.4.2 Xây dựng mơ hình tốn học thể mối quan hệ thời gian gia nhiệt nhiệt độ phôi dập vuốt 77 3.5 Kiểm chứng độ xác FLC vật liệu SPCC nhiệt độ khác thông qua mô thực nghiệm 79 3.5.1 Kiểm chứng FLC nhiệt độ phòng 79 3.5.2 Kiểm chứng FLC nhiệt độ khác nhau…………………………………….81 3.6 Nghiên cứu ảnh hưởng lực chặn phôi, mức độ dập vuốt, bán kính cong chày nhiệt độ đến chiều cao tạo hình chiều dày phân bố chi tiết dạng cốc trụ 82 iv 3.6.1 Nghiên cứu ảnh hưởng lực chặn phôi (FBH) đến chiều cao tạo hình chi tiết dạng cốc trụ 82 3.6.2 Nghiên cứu ảnh hưởng mức độ dập vuốt (Mt) đến chiều cao tạo hình chi tiết dạng cốc trụ .85 3.6.3 Nghiên cứu ảnh hưởng bán kính cong chày Rp đến chiều cao tạo hình chi tiết dạng cốc trụ 87 3.6.4 Nghiên cứu ảnh hưởng nhiệt độ phôi T (0C) đến chiều cao tạo hình chi tiết dạng cốc trụ .89 3.7 Nghiên cứu ảnh hưởng nhiệt độ đến chiều dày phân bố chi tiết dạng cốc trụ dập vuốt 91 3.7.1 Nghiên cứu ảnh hưởng nhiệt độ đến chiều dày phân bố chi tiết dạng cốc trụ nhiệt độ phòng 93 3.7.2 Nghiên cứu ảnh hưởng nhiệt độ đến chiều dày phân bố chi tiết dạng cốc trụ nhiệt độ khác 95 3.7.3 So sánh phân bố chiều dày chi tiết dạng cốc trụ thông qua mô nhiệt độ khác 97 KẾT LUẬN CHƯƠNG 99 CHƯƠNG XÂY DỰNG MƠ HÌNH TỐN THỂ HIỆN ẢNH HƯỞNG CỦA MỘT SỐ THÔNG SỐ ĐẾN CHIỀU CAO TẠO HÌNH VÀ CHIỀU DÀY PHÂN BỐ CHI TIẾT DẠNG CỐC TRỤ KHI DẬP VUỐT 100 4.1 Thiết kế thực nghiệm .100 4.2 Điều kiện thực nghiệm 104 4.3 Nghiên cứu mối quan hệ thông số đầu vào thông số đầu 105 4.3.1 Phân tích mối quan hệ thơng số (FBH, Rp, Mt) đến chiều cao tạo hình (HR) chi tiết dạng cốc trụ dập vuốt 105 4.3.2 Phân tích mối quan hệ thông số (FBH, Rp, T) đến chiều cao tạo hình (HR2) chi tiết dạng cốc trụ dập vuốt 110 4.3.3 Phân tích mối quan hệ thông số (FBH, Rp, T) đến chiều dày phân bố (tP) chi tiết dạng cốc trụ dập vuốt .116 KẾT LUẬN CHƯƠNG 121 KẾT LUẬN VÀ HƯỚNG NGHIÊN CỨU TIẾP THEO 122 TÀI LIỆU THAM KHẢO 123 DANH MỤC CÁC CƠNG TRÌNH ĐÃ CƠNG BỐ CỦA LUẬN ÁN 129 PHỤ LỤC v DANH MỤC CÁC KÝ HIỆU VÀ CHỮ VIẾT TẮT d1 Diễn giải Đơn vị Forming Limit Curve – Đường cong giới hạn tạo hình Forming Limit Diagram – Sơ đồ đường cong giới hạn Fracture Forming Limit Diagram – Sơ đồ giới hạn tạo hình phá hủy Finite Element Method – Mô phương pháp phần tử hữu hạn Modified Maximum Force Criterion – Tiêu chuẩn trở lực tối đa sửa đổi Steel – Plate – Cold rolled – Carbon – Thép cán nguội Mơ hình Johnson – Cook Mơ hình Bao-Wierzbicki Plane Strain – Biến dạng phẳng Uniaxial Tension – Biến dạng kéo đơn trục Equibiaxial Tension – Biến dạng kéo đồng thời theo hai phương Stretch Forming – Biến dạng tạo hình kéo Pure Shear – Biến dạng cắt túy Uniaxial Compression – Biến dạng nén đơn trục Bán kính lượn cối mm Bán kính cong chày mm Khe hở chày cối mm Đường kính phơi ban đầu mm Bán kính phơi mm Bán kính vành chi tiết sau dập vuốt mm Bán kính chày dập vuốt mm Đường kính chày dập vuốt mm Đường kính phơi sau lần dập mm d2 Đường kính bán thành phẩm sau lần dập thứ hai mm dn Đường kính bán thành phẩm cho lần dập mm Mt Mức độ dập vuốt HR Chiều cao tạo hình chi tiết bắt đầu xuất rách (Gọi tắt “chiều cao tạo hình”) Chiều cao tạo hình chi tiết dạng cốc trụ thực nghiệm Chiều cao tạo hình chi tiết dạng cốc trụ mơ Chiều cao tạo hình chi tiết dạng cốc trụ (sử dụng FLC xây dựng dựa mơ hình Voce mơ số) Ký hiệu FLC FLD FFLD FEM MMFC SPCC J-C B-W P.S U.T B.T S.F S.P U.C Rd Rp Wc D0 R0 R rp dp HR-TN HR-MP Hv vi mm mm mm mm mm Ft Chiều cao tạo hình chi tiết dạng cốc trụ (sử dụng FLC xây dựng dựa mơ hình Swift mơ số) Chiều cao tạo hình chi tiết dạng cốc trụ (sử dụng FLC xây dựng dựa mơ hình Kim-Tuan mô số) Sai lệch chiều cao tạo hình chi tiết dạng cốc trụ sử dụng mơ hình tương ứng Voce; Swift; Kim-Tuan Chiều cao tạo hình chi tiết dạng cốc trụ sử dụng FLC khác Chiều cao tạo hình chi tiết dạng cốc trụ sử dụng FLC theo phương pháp tỷ lệ Chiều cao tạo hình chi tiết dạng cốc trụ sử dụng FLC theo phương pháp đề xuất Sai lệch chiều cao tạo hình chi tiết dạng cốc trụ sử dụng FLC theo phương pháp tỷ lệ phương pháp đề xuất Sai lệch chiều cao tạo hình chi tiết dạng cốc trụ thay đổi mức độ dập vuốt Sai lệch chiều cao tạo hình chi tiết dạng cốc trụ thay đổi bán kính cong chày Sai lệch chiều cao tạo hình chi tiết dạng cốc trụ thay đổi lực chặn phôi Sai lệch chiều cao tạo hình chi tiết dạng cốc trụ thay đổi nhiệt độ Chiều cao tạo hình chi tiết dạng cốc trụ xác định từ mơ hình nhiệt độ phịng Chiều cao tạo hình chi tiết dạng cốc trụ xác định từ mơ hình nhiệt độ khác Chiều dày phân bố chi tiết dạng cốc trụ Chiều dày phôi Chiều dày lớn phần vành ổ biến dạng Sai lệch phân bố chiều dày chi tiết dạng cốc trụ Chiều dày chi tiết dạng cốc trụ mô Chiều dày chi tiết dạng cốc trụ thực nghiệm Chiều dày nhỏ chi tiết dạng cốc trụ Chiều dày trung bình chi tiết dạng cốc trụ Lực dập vuốt FBH Lực chặn phơi kN µp Hệ số ma sát chày phơi µh Hệ số ma sát chặn phơi µd Hệ số ma sát cối phôi E Mô đun đàn hồi kim loại vii Hs Hkt ∆𝐻𝑣 ; ∆𝐻𝑆 ; ∆𝐻𝑘𝑡 HFLC HSTL HSDX ∆𝐻𝑆𝑇𝐿 ; ∆𝐻𝑆𝐷𝑋 ∆𝐻𝑀𝑡 ∆𝐻𝑅𝑝 ∆𝐻𝐹 ∆𝐻𝑇 HR1-PT HR2-PT tP t tk ∆𝑡𝑝 tMP tTN tmin tTB mm % mm mm mm % % % % % mm mm mm mm mm % mm mm mm mm kN MPa Vp 𝜀1 , 𝜀2  ;  f ;  eq 𝜀̇  𝑇𝑒𝑚𝑝 𝜀̅𝑓 , 𝜀̅𝑓𝑅𝑜𝑜𝑚  , eq 1; ;  * Tốc độ dụng cụ gây biến dạng Biến dạng phụ Biến dạng tương đương Tốc độ biến dạng Biến dạng hướng tiếp tuyến mm/s s-1 Biến dạng phá hủy tương đương tương ứng nhiệt độ cao nhiệt độ phòng Ứng suất tương đương MPa Ứng suất theo phương MPa Tỷ lệ biến dạng Ứng suất tương đương MPa Ứng suất thủy tĩnh 𝜎𝐻 Ứng suất tương đương Von Mises 𝜎𝑒𝑞 F, G, H, L, M Các tham số bất đẳng hướng Hill N Hệ số dị hướng theo phương cán, phương ngang vng r0 ; r90 ; r45 góc với hướng cán theo phương 450 Tỷ lệ ứng suất theo hai phương  Biểu thị hàm cứng hóa H =H ( ) H' 𝑅𝑚  d , d1 , c0 , c1 Biểu thị độ dốc đường cong cứng hóa Giá trị Lankford Chỉ số ứng suất theo ba phương Các hệ số phương trình D1 , D2 , D3 Các hệ số phương trình T Ts1 Ts2 Ts3 Ts4 Ts5 ANOVA S/N MSD yi y0 Su Nhiệt độ Cảm biến nhiệt độ cối dập vuốt Cảm biến nhiệt độ chặn phôi Cảm biến nhiệt độ đế cối dập vuốt Cảm biến nhiệt độ chày dập vuốt Cảm biến nhiệt độ cối Analysis of Variance – Phân tích phương sai Signal – to – Nooise ratio – Tỷ số tín hiệu nhiễu Bình phương độ lệch Các giá trị kiểm tra thực nghiệm r Giá trị tiêu chuẩn giá trị mục tiêu Phương sai dư Sll Phương sai lặp F Hệ số Fisher viii MPa MPa MPa MPa C C C C C C 0,573 0,554 0,563 0,572 0,588 Bảng Hình ảnh chiều dày vị trí chi tiết dạng cốc trụ dập 1500C Nhiệt độ 1500C Các vị trí đo chiều dày Chiều dày (mm) 0,711 0,662 0,613 0,572 Hình ảnh đo 5 0,544 0,562 0,568 0,573 Bảng Hình ảnh chiều dày vị trí chi tiết dạng cốc trụ dập 2500C Nhiệt độ 2500C Các vị trí đo chiều dày Chiều dày (mm) 0,718 0,672 0,614 0,576 Hình ảnh đo 0,542 0,564 0,571 0,576 Bảng Hình ảnh chiều dày vị trí chi tiết dạng cốc trụ dập với thông số tối ưu Các vị trí đo chiều dày Bộ thơng số tối ưu Chiều dày Hình ảnh đo (mm) 0,686 0,642 0,618 0,584 0,571 0,576 0,583 0,585 10 Hình ảnh mơ số ảnh hưởng nhiệt độ đến chiều cao tạo hình chi tiết dạng cốc trụ T11=250C T12=1500C T13=2500C T14=3000C T15=3500C T16=4000C Hình Chiều cao tạo hình mô nhiệt độ khác Kết mô dập vuốt để đo chiều dày chi tiết dạng cốc trụ L1 11 L2 L3 L4 12 L5 L6 L7 13 L8 L9 Hình Hình ảnh mơ để xác định chiều dày chi tiết dạng cốc trụ Bảng Kết đo chiều dày chi tiết dạng cốc trụ thực nghiệm Vị trí đo chiều dày (mm) Taguchi L9 L1 L2 L3 L4 L5 L6 L7 L8 L9 0,692 0,691 0,675 0,688 0,675 0,697 0,695 0,686 0,682 0,657 0,665 0,647 0,631 0,643 0,651 0,659 0,645 0,643 0,610 0,617 0,608 0,605 0,607 0,615 0,612 0,602 0,604 0,584 0,596 0,597 0,594 0,593 0,598 0,581 0,581 0,596 0,562 0,551 0,557 0,558 0,558 0,561 0,551 0,564 0,563 0,570 0,560 0,569 0,561 0,574 0,568 0,570 0,572 0,567 0,570 0,562 0,581 0,582 0,585 0,568 0,573 0,576 0,571 0,574 0,573 0,582 0,584 0,587 0,569 0,574 0,580 0,579 TB 0,602 0,602 0,602 0,600 0,603 0,603 0,602 0,601 0,601 14 Chương chình Matlab xây dựng hàm chiều cao tạo hình với thơng số (FBH, Rp , M t ) clear all;clc; syms a b c d X Y Z A1=10;A2=12.5;A3=15; B1=4;B2=6;B3=8; C1=2.1;C2=2.25;C3=2.4; A=[A1 A1 A1 A2 A2 A2 A3 A3 A3]; B=[B1 B2 B3 B1 B2 B3 B1 B2 B3]; C=[C1 C2 C3 C2 C3 C1 C3 C1 C2]; %% K=(A3-A1)./25; L=(B3-B1)./25; M=(C3-C1)./25; %% x=meshgrid(A1:K:A3); y=meshgrid(B1:L:B3); [X,Y]=meshgrid(x,y); Z=C1; HZ=2464.3.*X.^(-0.4875).*Y.^0.219862.*Z.^(-4.3682); figure(1) surfc(X,Y,HZ), xlabel('FBH (KN)'), ylabel('Rp (mm)'), zlabel('H (mm)'); hold on grid on shading interp color map(gray); %colorbar %% x= meshgrid(A1:K:A3); z=meshgrid(C1:M:C3); [X,Z]=meshgrid(x,z); Y=B3; HY=2464.3.*X.^(-0.4875).*Y.^0.219862.*Z.^(-4.3682); figure(2) surfc(X,Z,HY), xlabel('FBH (KN)'), ylabel('Mt'), zlabel('H (mm)'); hold on grid on shading interp color map(gray); %% z= meshgrid(C1:M:C3); y=meshgrid(B1:L:B3); 15 [Y,Z]=meshgrid(y,z); X=A1; HX=2464.3.*X.^(-0.4875).*Y.^0.219862.*Z.^(-4.3682); figure(3) surfc(Y,Z,HX), xlabel('Rp (mm)'), ylabel('Mt'), zlabel('H (mm)'); hold on grid on shading interp color map(gray); HTN=[41.8 35.7 28.8 28.9 22.2 42.8 19.4 40.1 29.5] X=A;Y=B;Z=C; HNoisuy=2464.3.*X.^(-0.4875).*Y.^0.219862.*Z.^(-4.3682) Giatritrungbinhthinghiem=(HTN(1)+HTN(2)+HTN(3)+HTN(4)+HTN(5)+H TN(6)+HTN(7)+HTN(8)+HTN(9))/9; Giatritrungnoisuy=HNoisuy ynho=HTN-Giatritrungbinhthinghiem; ynhomu=HNoisuy-Giatritrungbinhthinghiem; enho=HTN-HNoisuy tongenhobinh=enho(1)^2+enho(2)^2+enho(3)^2+enho(4)^2+enho(5)^2 +enho(6)^2+enho(7)^2+enho(8)^2+enho(9)^2; TSS=ynho(1)^2+ynho(2)^2+ynho(3)^2+ynho(4)^2+ynho(5)^2+ynho(6)^ 2+ynho(7)^2+ynho(8)^2+ynho(9)^2; ESS=ynhomu(1)^2+ynhomu(2)^2+ynhomu(3)^2+ynhomu(4)^2+ynhomu(5)^ 2+ynhomu(6)^2+ynhomu(7)^2+ynhomu(8)^2+ynhomu(9)^2; SSE=enho(1)^2+enho(2)^2+enho(3)^2+enho(4)^2+enho(5)^2+enho(6)^ 2+enho(7)^2+enho(8)^2+enho(9); n=9;p=4; MSE=SSE/(n-p); %=SSE/(n-p) n=9;p=4 MSM=ESS/(p-1); %=SSM/(p-1) MST=TSS/(n-1); %=SST/(n-1) Phuong sai tai sinh; f1=n-p; f2=p-1 F(f2,f1) %% R binh phuong Rbinhphuong1=1-(SSE/TSS) Rbinhphuonghieuchinh=1-(MSE/MST) %% Fisher %Fisher1=MSM/MST %% phuong sai du phuongsaidu=SSE/(n-p)% phuong sai tuong thich OK %% Phuong sai lap dd=43.26;ee=44.6;ff=45.5; tb=(dd+ee+ff)/3; phuongsailap=((dd-tb)^2+(ee-tb)^2+(ff-tb)^2)/2 Fisher2=phuongsaidu/phuongsailap %% f(alpha,p-1,n-p) 16 Chương chình Matlab xây dựng hàm chiều cao tạo hình với thông số (FBH, Rp, T) clear all;clc; syms a b c d X Y Z A1=10;A2=12.5;A3=15; B1=4;B2=6;B3=8; C1=25;C2=150;C3=250; A=[A1 A1 A1 A2 A2 A2 A3 A3 A3]; B=[B1 B2 B3 B1 B2 B3 B1 B2 B3]; C=[C1 C2 C3 C2 C3 C1 C3 C1 C2]; %% K=(A3-A1)./25; L=(B3-B1)./25; M=(C3-C1)./25; %% x=meshgrid(A1:K:A3); y=meshgrid(B1:L:B3); [X,Y]=meshgrid(x,y); Z=C1; HZ=84.586.*X.^(-0.795).*Y.^0.229.*Z.^(0.194); figure(1) surfc(X,Y,HZ), xlabel('FBH (KN)'), ylabel('Rp (mm)'), zlabel('H (mm)'); hold on grid on shading interp color map(gray); %colorbar %% x= meshgrid(A1:K:A3); z=meshgrid(C1:M:C3); [X,Z]=meshgrid(x,z); Y=B1; HY=84.586.*X.^(-0.795).*Y.^0.229.*Z.^(0.194); figure(2) surfc(X,Z,HY), xlabel('FBH (KN)'), ylabel('T (C)'), zlabel('H (mm)'); hold on grid on shading interp color map(gray); %% z= meshgrid(C1:M:C3); y=meshgrid(B1:L:B3); 17 [Y,Z]=meshgrid(y,z); X=A1; HX=84.586.*X.^(-0.795).*Y.^0.229.*Z.^(0.194); figure(3) surfc(Y,Z,HX), xlabel('Rp (mm)'), ylabel('T (C)'), zlabel('H (mm)'); hold on grid on shading interp color map(gray); HTN=[33.2 56.8 61.8 40.4 49.5 32.6 38.2 27.3 41.9] X=A;Y=B;Z=C; HNoisuy=84.586.*X.^(-0.795).*Y.^0.229.*Z.^(0.194); Giatritrungbinhthinghiem=(HTN(1)+HTN(2)+HTN(3)+HTN(4)+HTN(5)+H TN(6)+HTN(7)+HTN(8)+HTN(9))/9; Giatritrungnoisuy=(HNoisuy(1)+HNoisuy(2)+HNoisuy(3)+HNoisuy(4) +HNoisuy(5)+HNoisuy(6)+HNoisuy(7)+HNoisuy(8)+HNoisuy(9))/9; ynho=HTN-Giatritrungbinhthinghiem; ynhomu=HNoisuy-Giatritrungbinhthinghiem; enho=HTN-HNoisuy tongenhobinh=enho(1)^2+enho(2)^2+enho(3)^2+enho(4)^2+enho(5)^2 +enho(6)^2+enho(7)^2+enho(8)^2+enho(9)^2; TSS=ynho(1)^2+ynho(2)^2+ynho(3)^2+ynho(4)^2+ynho(5)^2+ynho(6)^ 2+ynho(7)^2+ynho(8)^2+ynho(9)^2; ESS=ynhomu(1)^2+ynhomu(2)^2+ynhomu(3)^2+ynhomu(4)^2+ynhomu(5)^ 2+ynhomu(6)^2+ynhomu(7)^2+ynhomu(8)^2+ynhomu(9)^2; SSE=enho(1)^2+enho(2)^2+enho(3)^2+enho(4)^2+enho(5)^2+enho(6)^ 2+enho(7)^2+enho(8)^2+enho(9); n=9;p=4; MSE=SSE/(n-p); %=SSE/(n-p) n=9;p=4 MSM=ESS/(p-1); %=SSM/(p-1) MST=TSS/(n-1); %=SST/(n-1) Phuong sai tai sinh; f1=n-p; f2=p-1 F(f2,f1) %% R binh phuong Rbinhphuong1=1-(SSE/TSS) Rbinhphuonghieuchinh=1-(MSE/MST) %% Fisher %Fisher1=MSM/MST %% phuong sai du phuongsaidu=SSE/(n-p)% phuong sai tuong thich OK %% Phuong sai lap dd=63.1;ee=62.3;ff=64.9; tb=(dd+ee+ff)/3; phuongsailap=((dd-tb)^2+(ee-tb)^2+(ff-tb)^2)/2 Fisher2=phuongsaidu/phuongsailap %% f(alpha,p-1,n-p) 18