(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 LUAN VAN CHAT LUONG download : add luanvanchat@agmail.com 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 LUAN VAN CHAT LUONG download : add luanvanchat@agmail.com 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 tơ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 LUAN VAN CHAT LUONG download : add luanvanchat@agmail.com 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ấ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 Luyện Thế Thạnh ii LUAN VAN CHAT LUONG download : add luanvanchat@agmail.com 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ố 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 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 LUAN VAN CHAT LUONG download : add luanvanchat@agmail.com 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 LUAN VAN CHAT LUONG download : add luanvanchat@agmail.com 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 LUAN VAN CHAT LUONG download : add luanvanchat@agmail.com 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 mm mm mm mm vi LUAN VAN CHAT LUONG download : add luanvanchat@agmail.com 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 LUAN VAN CHAT LUONG download : add luanvanchat@agmail.com 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 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 s-1 MPa MPa MPa MPa C C C C C C viii LUAN VAN CHAT LUONG download : add luanvanchat@agmail.com 0,573 0,554 0,563 0,572 0,588 LUAN VAN CHAT LUONG download : add luanvanchat@agmail.com 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 LUAN VAN CHAT LUONG download : add luanvanchat@agmail.com 0,544 0,562 0,568 0,573 LUAN VAN CHAT LUONG download : add luanvanchat@agmail.com 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 LUAN VAN CHAT LUONG download : add luanvanchat@agmail.com 0,542 0,564 0,571 0,576 LUAN VAN CHAT LUONG download : add luanvanchat@agmail.com 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 LUAN VAN CHAT LUONG download : add luanvanchat@agmail.com 0,571 0,576 0,583 0,585 10 LUAN VAN CHAT LUONG download : add luanvanchat@agmail.com 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 LUAN VAN CHAT LUONG download : add luanvanchat@agmail.com L2 L3 L4 12 LUAN VAN CHAT LUONG download : add luanvanchat@agmail.com L5 L6 L7 13 LUAN VAN CHAT LUONG download : add luanvanchat@agmail.com 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 LUAN VAN CHAT LUONG download : add luanvanchat@agmail.com 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 LUAN VAN CHAT LUONG download : add luanvanchat@agmail.com [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 LUAN VAN CHAT LUONG download : add luanvanchat@agmail.com 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 LUAN VAN CHAT LUONG download : add luanvanchat@agmail.com [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 LUAN VAN CHAT LUONG download : add luanvanchat@agmail.com ... .90 Hình 34 Ảnh hưởng nhiệt độ đến chi? ??u cao tạo hình chi tiết dạng cốc trụ dập vuốt vật liệu SPCC 90 Hình 35 Chi? ??u cao tạo hình chi tiết dạng cốc trụ mô dập vuốt nhiệt độ từ 25 0C đến. .. 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... 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ừ