Nghiên cứu ảnh hưởng của các thông số công nghệ và nhiệt độ 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 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 TẬP THỂ HƯỚNG DẪN KHOA HỌC NGHIÊN CỨU SINH (Đã ký) (Đã ký) (Đã ký) GS.TS Nguyễn Đức Toàn GS.TSKH Bành Tiến Long Luyện Thế Thạnh ii 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 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ố q 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 ngồi nước 26 KẾT LUẬN CHƯƠNG 33 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 34 iii 2.1 Mơ hình thuộc tính phá huỷ vật liệu………………………………………………34 2.1.1 Phá huỷ dẻo vật liệu…………………………………………………………… 34 2.1.2 Mơ hình phá hủy vật liệu……………………………………………………… 36 2.1.3 Mơ hình thuộc tính vật liệu 43 2.1.4 Xác định tính vật liệu SPCC sử dụng nghiên cứu 44 2.2 Xây dựng FLC vật liệu SPCC nhiệt độ phòng nhiệt độ khác nhau……48 2.2.1 Xây dựng FLC vật liệu SPCC nhiệt độ phòng 49 2.2.2 Xây dựng FLC vật liệu SPCC nhiệt độ khác nhau………………….54 KẾT LUẬN CHƯƠNG 62 CHƯƠNG NGHIÊN CỨU Ả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 63 3.1 Sơ đồ hệ thống thực nghiệm 63 3.2 Thiết lập thực nghiệm……………………………………………………………… 65 3.2.1 Vật liệu thực nghiệm……………………………………………………………65 3.2.2 Phôi dụng cụ phục vụ thực nghiệm……………………………………65 3.2.3 Máy thực nghiệm……………………………………………………………… 66 3.2.4 Bộ khuôn thực nghiệm……………………………………………………… 67 3.2.5 Thiết bị gia nhiệt thu thập liệu nhiệt……………………………….68 3.2.6 Thiết bị đo .70 3.3 Mô số gia công dập vuốt 72 3.3.1 Mơ hình phần tử hữu hạn (FEM)……………………………………………… 72 3.3.2 Thiết lập thông số mô thực nghiệm………………………………73 3.4 Nghiên cứu trình gia nhiệt dập vuốt thực nghiệm……………… 75 3.4.1 Sơ đồ thực nghiệm………………………………………………………………75 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 76 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 78 3.5.1 Kiểm chứng FLC nhiệt độ phòng 78 3.5.2 Kiểm chứng FLC nhiệt độ khác nhau……………………………………… 80 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ụ .81 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ụ 81 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ụ 84 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ụ .86 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ụ 88 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 90 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 .92 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ụ .94 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 96 KẾT LUẬN CHƯƠNG 98 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 99 4.1 Thiết kế thực nghiệm 99 4.2 Điều kiện thực nghiệm .103 4.3 Nghiên cứu mối quan hệ thông số đầu vào thơng số đầu 104 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 (HR1) chi tiết dạng cốc trụ dập vuốt 104 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 109 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 115 KẾT LUẬN CHƯƠNG 120 KẾT LUẬN VÀ HƯỚNG NGHIÊN CỨU TIẾP THEO 121 TÀI LIỆU THAM KHẢO 122 DANH MỤC CÁC CƠNG TRÌNH ĐÃ CƠNG BỐ CỦA LUẬN ÁN 128 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 (Tỷ lệ đường kính phơi ban đầu đường kính chi tiết sau dập vuốt) 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 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 𝑇𝑒𝑚𝑝 mm/s 𝜀̇ Tốc độ biến dạng Biến dạng hướng tiếp tuyến , 𝜀̅𝑓𝑅𝑜𝑜𝑚 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  𝜀̅𝑓 Tốc độ dụng cụ gây biến dạng Biến dạng phụ Biến dạng tương đương  , eq 1; ;  * Ứ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 vuông 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 s-1 MPa MPa MPa MPa C C C C C C 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 dày (mm) 0,711 0,662 0,613 0,572 Chiều 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 chiều dày (mm) Hình ảnh đo dày 0,718 0,672 0,614 0,576 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 Bộ thơng số tối ưu Chiều chiều dày (mm) 0,686 0,642 0,618 0,584 Hình ảnh đo dày 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 Hình ảnh mơ số ảnh hưởng tốc độ dụng cụ gây biến dạng ( tốc độ ép) đến chiều cao tạo hình chi tiết dạng cốc trụ V= 5mm/s V= 10 mm/s 11 V= 15mm/s V= 20 mm/s Hình Chiều cao tạo hình mơ với tốc độ ép khác Kết mô dập vuốt để đo chiều dày chi tiết dạng cốc trụ L1 L2 12 L3 L4 L5 13 L6 L7 L8 14 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 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; 15 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); [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) 16 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) 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]; %% 17 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); [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); 18 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) 19