Nghiên cứu hàm lượng hidro trong kim loại mối hàn khi hàn tự động dưới lớp thuốc hàn gốm hệ ba zơ

MẪU BÌA LUẬN VĂN CĨ IN CHỮ NHŨ VÀNG Khổ 210 x 297 mm BỘ GIÁO DỤC VÀ ĐÀO TẠO TRƯỜNG ĐẠI HỌC BÁCH KHOA HÀ NỘI DƯƠNG VĂN HỘI HỌ VÀ TÊN: DƯƠNG VĂN HỘI CHUYÊN NGÀNH: CÔNG NGHỆ HÀN NGHIÊN CỨU HÀM LƯỢNG HIĐRÔ TRONG KIM LOẠI MỐI HÀN KHI HÀN TỰ ĐỘNG DƯỚI LỚP THUỐC HÀN GỐM HỆ BA ZƠ LUẬN VĂN THẠC SĨ KHOA HỌC: CÔNG NGHỆ HÀN KHOÁ: 2009 Hà Nội năm 2011 MẪU TRANG PHỤ BÌA LUẬN VĂN BỘ GIÁO DỤC VÀ ĐÀO TẠO TRƯỜNG ĐẠI HỌC BÁCH KHOA HÀ NỘI Họ tên: DƯƠNG VĂN HỘI NGHIÊN CỨU HÀM LƯỢNG HIĐRÔ TRONG KIM LOẠI MỐI HÀN KHI HÀN TỰ ĐỘNG DƯỚI LỚP THUỐC HÀN GỐM HỆ BA ZƠ Chuyên ngành : Công nghệ hàn LUẬN VĂN THẠC SĨ KHOA HỌC: CÔNG NGHỆ HÀN NGƯỜI HƯỚNG DẪN KHOA HỌC: TS VŨ HUY LÂN Hà Nội – Năm 2011 LUẬN VĂN TỐT NGHIỆP MỤC LỤC Trang Trang phụ bìa Lời cam đoan Danh mục ký hiệu, chữ viết tắt Danh mục bảng Danh mục hình vẽ, đồ thị MỞ ĐẦU Chương TỔNG QUAN NGHIÊN CỨU VỀ HÀM LƯỢNG CỦA HIĐRÔ TRONG KIM LOẠI MỐI HÀN KHI CHẾ TẠO VÀ SỬ DỤNG THUỐC HÀN GỐM BAZƠ 1.1 Tổng quan hàn tự động lớp thuốc 1.2 Ảnh hưởng hiđrô kim loại mối hàn hàn lớp thuốc bảo vệ 1.3 Tổng quan nghiên cứu hàm lượng hiđrô kim loại mối hàn chế tạo sử dụng thuốc hàn gốm bazơ giới 1.4 Việc nghiên cứu hàm lượng hiđrô kim loại mối hàn 10 chế tạo sử dụng thuốc hàn gốm bazơ Việt Nam Chương NGHIÊN CỨU ẢNH HƯỞNG CỦA HIĐRÔ TRONG KIM 11 LOẠI MỐI HÀN 2.1 Nguồn gốc hiđrô 11 2.1.1 Nguồn hiđrô từ môi trường xung quanh 11 2.1.2 Nguồn hiđrô từ vật hàn vật liệu hàn 11 2.1.3 Nguồn hiđrô tồn thuốc hàn 11 0B 1B 2B 2.2 Cơ chế ảnh hưởng 12 Người hướng dẫn: TS Vũ Huy Lân Học viên: Dương Văn Hội I LUẬN VĂN TỐT NGHIỆP Nghiên cứu giải pháp chủ yếu giảm hàm lượng hiđrô 2.3 19 kim loại mối hàn chế tạo thuốc hàn 2.3.1 Chuẩn bị vật liệu hàn 19 2.3.2 Sử dụng CaF2 thành phần thuốc hàn 19 2.3.3 Chế độ sấy sơ thiêu kết chế tạo thuốc hàn 20 2.3.4 Chế độ sấy trước sử dụng thuốc hàn 21 3B 4B 5B 6B Chương GIỚI THIỆU VỀ THÀNH PHẦN THUỐC HÀN GỐM VÀ ẢNH 22 HƯỞNG CỦA HUỲNH THẠCH ĐẾN HÀM LƯỢNG HIĐRÔ TRONG KIM LOẠI MỐI HÀN Giới thiệu thành phần mẻ liệu thuốc hàn gốm 3.1 22 3.1.1 Nhóm tạo khí 22 3.1.2 Nhóm tạo xỉ 22 3.1.3 Nhóm chất khử, hợp kim hóa biến tính kim loại mối hàn 22 3.1.4 Nhóm ổn định hồ quang 23 3.1.5 Nhóm chất tạo hình 23 3.1.6 Nhóm chất dính kết 23 Ảnh hưởng huỳnh thạch đến hàm lượng hiđrô kim 3.2 24 loại mối hàn với tạo xỉ 3.2.1 Các yếu tố ảnh hưởng đến hàm lượng hiđrô kim loại 24 mối hàn 3.2.2 Ảnh hưởng huỳnh thạch 24 Chương TÍNH TỐN SƠ BỘ THÀNH PHẦN MẺ LIỆU THUỐC 25 HÀN GỐM 4.1 Giới thiệu sơ lược xỉ hàn lựa chọn tạo xỉ hàn 25 4.1.1 Khái niệm xỉ hàn 25 4.1.2 Đặc tính xỉ hàn 25 4.1.2.1 Tính axit( bazơ) xỉ hàn Người hướng dẫn: TS Vũ Huy Lân 25 Học viên: Dương Văn Hội II LUẬN VĂN TỐT NGHIỆP 4.1.2.2 Hoạt tính hố học xỉ hàn 26 4.1.2.3 Độ nhớt 27 4.1.3 Một số đặc tính khác xỉ hàn 29 4.1.3.1 Sức căng bề mặt phân pha kim loại - xỉ hàn 29 4.1.3.2 Độ thẩm thấu khí xỉ hàn 29 4.1.4 Phân loại xỉ hàn 29 7B 4.1.4.1 Phân theo tính axit hay bazơ 29 4.1.4.2 Phân loại theo độ hoạt tính 30 4.1.4.3 Phân loại theo thay đổi độ nhớt 30 4.2 Cơ sở lựa chọn tạo xỉ cho mẻ liệu thuốc hàn gốm 30 4.2.1 Yêu cầu tạo xỉ 30 4.2.2 Một số tạo xỉ thông dụng 30 4.2.3 Chọn tạo xỉ cho thuốc hàn gốm 31 T T T 28T 8B 9B Lựa chọn nguyên liệu cho mẻ liệu thuốc hàn 32 4.3.1 Sơ lược nguồn nguyên liệu khoáng chất nước 32 4.3 10B Lựa chọn sơ thành phần mẻ liệu thuốc hàn gốm 33 Tính tốn thành phần mẻ liệu thuốc hàn 34 4.3.2 4.4 4.4.1 Giới thiệu thành phần hoá học tính thép làm dây hàn H- 34 08A 4.4.2 Tính tốn lượng ơxy sinh 34 4.4.3 Tính tốn lượng chất khử cần dùng 36 4.4.4 Tính tốn hàm lượng chất hợp kim hố 36 4.4.5 Tính tốn chất xỉ hàn 38 12B 13B 14B T 28T Người hướng dẫn: TS Vũ Huy Lân Học viên: Dương Văn Hội III LUẬN VĂN TỐT NGHIỆP 4.4.5.1 Hàm lượng TiO2 38 4.4.5.2 Từ Đá vôi (98% CaCO3, 1% SiO2, 1% Al2O3) 38 4.4.5.3 Từ quặng Bôxit (52% Al2O3, 5% SiO2) 39 4.4.5.4 Từ cao lanh (52% SiO2, 35% Al2O3, 1,25% Fe2O3) 39 4.4.5.5 Từ Fero – Si (75% Si ) 39 4.4.5.6 Từ Fero – Mn (80% Mn) 40 4.4.5.7 Từ trường thạch (Fenspat) (51,1%SiO2, 31%Al2O3, 40 4,3%K2O+Na2O) 4.4.5.8 Từ Huỳnh thạch (96% CaF , 2% SiO , 2%CaCO ) 40 4.4.5.9 Từ nước thủy tinh (25% SiO2, 11% K2O) 41 R Chương R R R R R QUY TRÌNH CƠNG NGHỆ CHẾ TẠO THUỐC HÀN GỐM 43 5.1 Sơ đồ ngun lý quy trình cơng nghệ chế tạo thuốc hàn gốm 43 5.2 Các bước chế tạo thuốc hàn gốm hệ bazơ 44 5.2.1 Tuyển chọn nguyên vật liệu 44 5.2.2 Nghiền, sàng, chuẩn bị nguyên liệu 44 5.2.3 Xử lý nguyên liệu 45 5.2.4 Phối liệu thành phần mẻ liệu thuốc hàn 47 5.2.5 Trộn nguyên liệu 47 5.2.5.1 Trộn khô 47 5.2.5.2 Trộn ướt 47 5.2.6 Tạo hạt 48 5.2.7 Sấy thuốc hàn (sấy sơ bộ, thiêu kết) 48 5.2.8 Sàng tuyển hạt, kiểm tra kích thích thước hạt, độ ẩm 49 5.2.9 Bao gói, kiểm tra thành phẩm nhập kho 49 Chương XÂY DỰNG KẾ HOẠCH THỰC NGHIỆM, TIẾN HÀNH 50 THÍ NGHIỆM VÀ XỬ LÝ SỐ LIỆU 6.1 Lựa chọn mô hình lập trình Người hướng dẫn: TS Vũ Huy Lân 50 Học viên: Dương Văn Hội IV LUẬN VĂN TỐT NGHIỆP 6.1.1 Lựa chọn mơ hình biến đầu vào, hàm mục tiêu 50 6.1.2 Xây dựng kế hoạch thực nghiệm phương pháp giải 53 6.1.2.1 Loại kế hoạch thực nghiệm 53 6.1.2.2 Xây dựng kế hoạch thực nghiệm cho trường hợp 54 nghiên cứu 6.1.3 Lập trình phần mềm xác định hệ số phương trình hồi quy 55 6.1.3.1 Phương pháp xác định hệ số mơ hình 55 6.1.3.2 Phần mềm lập trình xác định hệ số mơ hình 59 6.2 Xây dựng kế hoạch thực nghiệm tiến hành thí nghiệm 59 6.2.1 Mơ tả phương pháp thí nghiệm 59 6.2.2 Tiến hành thí nghiệm kết 64 15B 16B 6.3 Xử lý kết thực nghiệm 65 6.3.1 Kết xử lý số liệu 65 6.3.2 Biểu diễn đường đặc tính kết luận khoa học 65 6.3.3 Xác định giá trị biến đầu vào đảm bảo hàm lượng hiđrô yêu 67 17B 18B 19B cầu KẾT LUẬN VÀ KIẾN NGHỊ 69 TÀI LIỆU THAM KHAO 71 PHỤ LỤC 72 B B B Người hướng dẫn: TS Vũ Huy Lân Học viên: Dương Văn Hội V LUẬN VĂN TỐT NGHIỆP DANH MỤC CÁC KÍ HIỆU, CÁC CỤM TỪ VIẾT TẮT SỬ DỤNG TRONG LUẬN VĂN Ký hiệu Đơn vị Ý nghĩa H2 - Hiđrô PH Mpa Áp suất hyđrơ phân tử pha khí KH2 , H - Hằng số cân hỗn hợp hyđrô phân tử nguyên tử m [g] Khối lượng Tncxi o C Nhiệt độ nóng chảy xỉ Tnc o C Nhiệt độ nóng chảy kim loại R kimloai P P - H P P P Hoạt tính hóa học Ha - Hoạt tính hóa học của ơxit axit Hb - Hoạt tính hố học ơxit bazơ η kl - Độ nhớt kim loại nóng chảy η xi - Độ nhớt xỉ R - Hệ số tương quan mơ hình [R] - Hệ số tương quan cho phép R R Người hướng dẫn: TS Vũ Huy Lân Học viên: Dương Văn Hội VI LUẬN VĂN TỐT NGHIỆP DANH MỤC CÁC BẢNG BIỂU TT Nội dung Trang Bảng:2 Hàm lượng giới hạn hiđrô kim loại đắp kim 25 loại mối hàn Bảng 4.1 thành phần hóa học chủ yếu xỉ hàn sử dụng thuốc hàn 31 ESAB Bảng 4.2 Bảng thành phần hóa học vật liệu dùng để sản xuất thuốc 32 hàn Bảng 4.3 Thành phần hoá học dây hàn H – 08A 34 Bảng 4.4 Cơ tính dây hàn H – 08A 34 Bảng 4.5 Thành phần sơ mẻ liệu thuốc hàn gốm 37 Bảng 4.6 Quy đổi thành phần mẻ liệu kể phần khô nước thuỷ tinh 38 Bảng 4.7 Hàm lượng (%) chất mẻ liệu thuốc hàn 41 Bảng 5.1 Danh mục quặng, sa khoáng, fero hợp kim lựa chọn 44 Bảng 5.2 Tỷ lệ chất mẻ liệu qui đổi 45 Bảng 5.3 Tỷ lệ chất mẻ liệu khô 46 Bảng 6.1 Kế hoạch thực nghiệm bậc hai trực giao 53 Bảng 6.2 Giá trị khoảng biến thiên yếu tố 54 Bảng 6.3 Kế hoạch thực nghiệm kết thí nghiệm 54 Bảng 6.4 Bảng giá trị tổng quát 55 Bảng 6.5 Kính thước mẫu kiểm tra đệm 59 Bảng 6.6 Kế hoạch thực nghiệm kết thí nghiệm 64 Bảng 6.7 Kết tối ưu biến đầu vào hàm lượng hiđrô 68 Người hướng dẫn: TS Vũ Huy Lân Học viên: Dương Văn Hội VII LUẬN VĂN TỐT NGHIỆP DANH MỤC HÌNH VẼ, ĐỒ THỊ ,SƠ ĐỒ Nội dung TT Trang Hình 1.1 Hàn lớp thuốc bảo vệ Hình 1.2 Hàn dầm chữ I phương pháp hàn tự động lớp thuốc bảo vệ Hình 1.3 Sự phân bố hiđrơ thép Hình 1.4 Các trường hợp nứt kim loại mối hàn Hình 2.1 Mức độ phân ly hiđrô phân tử phụ thuộc vào nhiệt độ 14 Hình 2.2 Độ hồ tan hyđrô sắt phụ thuộc vào nhiệt độ, áp 15 suất pha khí ≤ 1kg / cm Hình 2.3 Ảnh hưởng nhiệt độ áp suất hyđrơ pha khí 16 đến độ hồ tan sắt lỏng Hình 2.4 Hàm lượng hiđrô kim loại mối hàn phụ thuộc vào hàm 18 lượng thuốc hàn Hình 4.1 Sự phụ thuộc độ nhớt xỉ hàn vào nhiệt độ 28 Hình 5.1 Sơ đồ nguyên lý quy trình cơng nghệ chế tạo thuốc hàn gốm 43 Hình 6.1 Sơ đồ tổng quát 52 Hình 6.2 Sơ đồ nghiên cứu đề tài 52 Hình 6.3 Sơ đồ thuật tốn tính hệ số phương trình hồi quy 58 Hình 6.4 Hình dạng kích thước mẫu hàn thí nghiệm 59 Hình 6.5 Cân mẫu thí nghiệm trước hàn 60 Hình 6.6 Thuốc hàn mẫu thí nghiệm 60 Hình 6.7 Tiến hành hàn mẫu thí nghiệm 61 Hình 6.8 Mẫu thí nghiệm sau hàn 61 Người hướng dẫn: TS Vũ Huy Lân Học viên: Dương Văn Hội VIII LUẬN VĂN TỐT NGHIỆP Hình 6.14 Hàm lượng hiđrô mối hàn phụ thuộc vào nhiệt độ sấy thiêu kết Hình 6.15 Hàm lượng hiđrơ mối hàn phụ thuộc vào thời gian sấy thiêu kết Nhận xét: Trên sở quan hệ đường đặc tính cho phép rút số kết luận quan đây: Người hướng dẫn: TS Vũ Huy Lân Học viên: Dương Văn Hội 66 LUẬN VĂN TỐT NGHIỆP - Từ phương trình hồi quy cho phép xây dựng đường đặc tính phản ánh phụ thuộc hàm lượng hiđrô mối hàn vào hàm lượng huỳnh thạch (CaF ) R R mẻ liệu thuốc hàn, nhiệt độ sấy thiêu kết thời gian thiêu kết - Đặc tính đường biểu diễn phản ánh rõ ảnh hưởng yếu tố đến hàm lượng hiđrô kim loại mối hàn vùng ảnh hưởng nhiệt, cho thấy tác dụng yếu tố - Mức độ ảnh hưởng yếu tố có khác nhau, hàm lượng huỳnh thạch (CaF ) R R mẻ liệu thuốc hàn nhiệt độ sấy thiêu kết ảnh hưởng mạnh đến hàm lượng hiđrơ mối hàn Cịn mức độ ảnh hưởng thời gian thiêu kết yếu yếu tố - Khi yếu tố đạt đến giá trị tới hạn, mức độ ảnh hưởng thay đổi Các kết phù hợp với lý thuyết 6.3.3 Xác định giá trị biến đầu vào đảm bảo hàm lượng hiđrô yêu cầu Nguyên tắc xác định giá trị yếu tố hợp lý: Xác định giá trị yếu tố hợp lý (tối ưu) từ mơ hình đảm bảo hàm lượng hiđrơ theo u cầu – tức đảm bảo tiêu chí kinh tế - kỹ thuật Do vậy, việc giải toán phải đáp ứng điều kiện ràng buộc Bài toán quy hoạch đa mục tiêu trường hợp tổng quát với n biến, m buộc, p mục tiêu là: Max(min) j k ( x1 , x , , x n ), k = 1,2, , p  g ( x1 , x , , x n ) ≤ bi  x j ≥ i = 1,2,…,m; j = 1,2,…,n Ta ký hiệu miền buộc miền D Trong toán tối ưu mục tiêu, việc giải tốn tìm giá trị biến đầu vào thỏa mãn giá trị hàm mục tiêu cho trước Người hướng dẫn: TS Vũ Huy Lân Học viên: Dương Văn Hội 67 LUẬN VĂN TỐT NGHIỆP 1) Các ràng buộc biến đầu vào: Là giới hạn biến đầu vào: ≤ X ≤ 14 ; R R 500 ≤ X ≤ 700 ; R R 40 ≤ X ≤ 80 R R 2) Xác định hàm mục tiêu cần tối ưu: Trong trường hợp hàm mục tiêu hàm lượng hiđrô mối hàn, nhiên phụ thuộc vào tiêu tính, mà chủ yếu tiêu độ dai va đập độ giãn dài tương đối kim loại mối hàn mà lựa chọn hàm lượng hiđrô thấp đến mức độ yêu cầu (rất thấp, thấp, trung bình,…) theo giới hạn Trong trường hợp nguyên tắc tối ưu vừa đảm bảo tính kỹ thuật (hàm lượng hiđrơ thấp ≤ cm3/100g) tính kinh tế : P P Y ≤ cm3/100g ; P P Và nhiệt độ thiêu kết thời gian sấy thiêu kết tối thiểu Sau chạy chương trình tối ưu xác định kết sau: Bảng 6.7 Kết tối ưu biến đầu vào hàm lượng hiđrô Thời gian sấy thiêu kết, (phút) X3 60 Hàm lượng huỳnh Nhiệt độ sấy, (°C) X2 thạch, (%) X R R R 10 700 Hàm lượng hiđrô mối hàn, cm3/100g 2,93 P P Với giá trị yếu tố đầu vào tìm tiến hành thí nghiệm kiểm tra, kết kiểm tra cho thấy giá trị tính tốn hàm lượng hiđrơ mối hàn so với thực nghiệm đủ xác, sai lệch 5% Người hướng dẫn: TS Vũ Huy Lân Học viên: Dương Văn Hội 68 LUẬN VĂN TỐT NGHIỆP KẾT LUẬN VÀ KIẾN NGHỊ Đề tài nghiên cứu ảnh hưởng yếu tố gồm có thành phần thuốc hàn chế độ sấy thiêu kết đến hàm lượng hiđrô mối hàn vùng ảnh hưởng nhiệt chế tạo thuốc hàn thiêu kết hệ bazơ để hàn tự động lớp thuốc kết cấu thép cacbon thấp Trên sở kết tính tốn thí nghiệm tiên nghiệm chọn cấu tử chủ yếu thành phần mẻ liệu thuốc hàn Sau giữ cố định cấu tử chuyển sang nghiên cứu thay đổi huỳnh thạch thông số chế độ thiêu kết ảnh hưởng đến hàm lượng hiđrơ mối hàn Trong q trình nghiên cứu ứng dụng phương pháp quy hoạch thực nghiệm để xây dựng mơ hình, xác lập mối quan hệ yếu tố đầu vào hàm lượng hiđrơ mối hàn Từ phương trình hồi quy xây dựng được, cho phép xác định giá trị tối ưu biến đầu vào đảm bảo hàm lượng hiđrô theo yêu cầu có xét đến tiêu kinh tế Kết chủ yếu đề tài gồm nội dung sau đây: - Đề tài nêu rõ nguồn gốc hiđrô mối hàn dạng tồn nước chứa vật liệu hàn điều kiện nhiệt độ sấy, nhiệt độ thiêu kết cần thiết khử nước (độ ẩm) trong loại thuốc hàn khác Trên sở lựa chọn khoảng nhiệt độ đầu vào cần thiết công đoạn sấy thiêu kết chế tạo thuốc hàn hệ bazơ - Đề tài tính toán lựa chọn tạo xỉ hợp lý cho thuốc hàn tương đương với loại OK Flux10.80 theo tiêu chuẩn hãng ESAB 1904-2004 (Thụy Điển) với hệ số bazơ B = 1,1 - Đề tài ứng dụng quy hoạch thực nghiệm, xác lập quan hệ hàm lượng hiđrô mối hàn phụ thuộc vào hàm lượng huỳnh thạch mẻ liệu thuốc hàn chế độ sấy thiêu kết Kết ý nghĩa khoa học cao tính mới, mà có giá trị thực tiễn việc sản xuất thuốc hàn vật liệu nước - Các biểu đồ xây dựng từ phương trình tốn học, phản ánh rõ nét ảnh hưởng huỳnh thạch thành phần thuốc hàn chế độ sấy thiêu kết đến hàm Người hướng dẫn: TS Vũ Huy Lân Học viên: Dương Văn Hội 69 LUẬN VĂN TỐT NGHIỆP lượng hiđrô, đồng thời giúp cho chuyên gia kỹ thuật ứng dụng thuận tiện vào việc tính toán điều chỉnh thành phần huỳnh thạch chế độ thiêu kết yêu cầu hàm lượng hiđrô thay đổi Kết nghiên cứu đề tài áp dụng với hệ xỉ bazơ B = 1,1, có điều kiện nên nghiên cứu mở rộng cho số loại thuốc hàn với hệ số bazơ phổ biến, để giúp cho sở sản xuất vật liệu hàn có đầy đủ tài liệu tham khảo Tuy nhiên, để kết đề tài sớm ứng dụng thực tiễn, cần nghiên cứu kết hợp với nội dung khác tiêu công nghệ hàn tính kỹ thuật thuốc hàn để chế tạo thuốc hàn đảm bảo tiêu kinh tế - kỹ thuật đề Người hướng dẫn: TS Vũ Huy Lân Học viên: Dương Văn Hội 70 LUẬN VĂN TỐT NGHIỆP TÀI LIỆU THAM KHẢO [1] - TS Vũ Huy Lân, TS Bùi Văn Hạnh, Giáo trình Vật liệu hàn, NXB Bách Khoa Hà Nội, 2010 [2] - TS Ngơ Lê Thơng, Cơng nghệ hàn điện nóng chảy (Tập 1&2), NXB Khoa học Kỹ thuật, Hà Nội 2004 [3] - TCVN 3223 : 2000, Que hàn thép cacbon thép hợp kim thấp [4] - Петров Г Л Сварочные материалы Машиностроение, Ленинград 1972 [5] - Пoтaпова Н Н Сварочные материалы для дуговой сварки Машиностроение, Москва 1989 [6] - Cпецэлектрод Covered electrodes for manual arc welding, depositing and cutting (Catalô que hàn hãng Spetselektrode, 2000) [7] - AWS Welding Handbook, 9th Edition, 2001 [8] - Kobelco Welding Handbook , Kobe Steel LTD., Japan 2004 [9] - Huyndai Welding Consumables, Huyndai Welding Co., LTD., Korea 2004 Người hướng dẫn: TS Vũ Huy Lân Học viên: Dương Văn Hội 71 LUẬN VĂN TỐT NGHIỆP PHỤ LỤC Phần mềm xác địch hệ số hệ phương trình hồi quy PROGRAM TTQHTN; uses crt; type matran = array[1 25,1 25] of real; var mt:matran; a:matran; b:matran; X: matran; Y: matran; cv1:matran; d:matran; d1:matran; d2:matran; yi:matran; ytb:matran; c,max:real; t:boolean; ck:char; m,n,i,j,k,q,l,a1,b1,k1: integer; f:file of matran; procedure matranX_Y;{ lap ma tran Y va X_ma tran thong tin } begin writeln(' Ma tran Y '); writeln; for i:=1 to m begin Y[i,1]:=mt[i,n]; write(Y[i,1]:8:2);writeln; end; writeln; Writeln(' Ma tran thong tin X'); writeln; for i:=1 to m d[i,1]:=1; for i:=1 to m for j:=1 to n d[i,j+1]:=mt[i,j]; for i:=1 to m Người hướng dẫn: TS Vũ Huy Lân Học viên: Dương Văn Hội 72 LUẬN VĂN TỐT NGHIỆP begin for j:=1 to n write(d[i,j]:8:2); writeln; end; end; (*******************************************************) Procedure matran_chuyenvi;{thanh lap ma tran chuyen vi X' } begin for i:=1 to n for j:=1 to m cv1[i,j]:=d[j,i]; writeln(' Ma tran chuyen vi '); writeln; for i:=1 to n begin for j:=1 to m write(cv1[i,j]:8:2); writeln; end; end; (*******************************************************) procedure tinh_tich_hai_ma_tran(a1,b1,b2:matran; var d1,d2:matran); { nhan ma tran A=X'.X va B=X'.Y } begin for i:=1 to n for k:=1 to begin d1[i,j]:=0; for j:=1 to m d1[i,k]:=d1[i,k]+a1[i,j]*b1[j,k]; end; writeln; Writeln(' Ma tran B '); writeln; for i:=1 to n begin for k:=1 to write(d1[i,k]:8:4); writeln; end; writeln; for i:=1 to n for k:=1 to n Người hướng dẫn: TS Vũ Huy Lân Học viên: Dương Văn Hội 73 LUẬN VĂN TỐT NGHIỆP begin d2[i,j]:=0; for j:=1 to m d2[i,k]:=d2[i,k]+a1[i,j]*b2[j,k]; end; writeln; writeln(' Ma tran A '); writeln; for i:=1 to n begin for k:=1 to n write(d2[i,k]:10:4); writeln; end; writeln; end; (******************************************************) procedure Gauss(a,b:matran;var x:matran);{ giai pt chinh tac =pp Gauss } begin t:=true; i:=1; while t begin if a[i,i]=0 then begin max:=0; m:=i ; for k:=i+1 to n if max< abs(a[k,i]) then begin m:=k; max:=abs(a[k,i]); end; if mi then begin for j:=i to n for l:=1 to begin c:=a[i,j]; a[i,j]:=a[m,j]; a[m,j]:=c; end; c:=b[i,l]; Người hướng dẫn: TS Vũ Huy Lân Học viên: Dương Văn Hội 74 LUẬN VĂN TỐT NGHIỆP b[i,l]:=b[m,l]; b[m,l]:=c; end; if m=i then begin t:=false; writeln(' Ma tran A suy bien '); exit; end; end; if a[i,i]0 then begin c:=1/a[i,i]; for j:=i+1 to n for l:=1 to a[i,j]:=a[i,j]*c; b[i,l]:=b[i,l]*c; for k:=i+1 to n begin for j:=i+1 to n for l:=1 to a[k,j]:=a[k,j]-a[i,j]*a[k,i]; b[k,l]:=b[k,l]-b[i,l]*a[k,i]; end; writeln; i:=i+1; if i=n+1 then t:=false; end; if i=n+1 then begin for i:=n-1 downto for j:=i+1 to n for l:= to b[i,l]:=b[i,l]-a[i,j]*b[j,l]; writeln(' Ma tran nghiem X '); writeln; Người hướng dẫn: TS Vũ Huy Lân Học viên: Dương Văn Hội 75 LUẬN VĂN TỐT NGHIỆP for i:=1 to n for l:=1 to begin x[i,l]:=b[i,l]; writeln(x[i,l]:8:4);writeln; end; end; end; end; (**********************************************************) procedure tinh_tich_ma_tran(ai,bi:matran; var yi:matran); {tinh gia tri cua ham so theo cong thuc thuc nghiem } begin for i:=1 to m for k:=1 to n begin yi[i,j]:=0; for j:=1 to n yi[i,k]:=yi[i,k]+ai[i,j]*bi[j,k]; end; writeln; Writeln(' Gia tri cua ham so theo cong thuc thuc nghiem '); writeln; for i:=1 to m begin for k:=1 to writeln(yi[i,k]:8:4); writeln end; writeln; end; procedure he_so_tuong_quan_R; var R1,s1,s2,s3:real; begin s1:=0; for i:=1 to m for j:=1 to s1:=s1+Y[i,j]; ytb[i,j]:=s1/m; writeln( 'Ytb=',ytb[i,j]:8:4); writeln; s2:=0; s3:=0; Người hướng dẫn: TS Vũ Huy Lân Học viên: Dương Văn Hội 76 LUẬN VĂN TỐT NGHIỆP for i:=1 to m for j:=1 to begin s2:= s2+ sqr(yi[i,j]- ytb[i,j]); s3:= s3+ sqr(Y[i,j]-ytb[i,j]); R1:=sqrt(s2/s3); end; Writeln('R=',R1:8:4); writeln; if (R1>=0) and(1>=R1) then Writeln('==> Thoa man ') else Writeln ( '==> Khong thoa man '); end; Begin {Chuong trinh chinh} clrscr; assign(f,'matran.dat'); rewrite(f); Write('Cho so hang cua bang so lieu:');readln(m); t:=true; while t begin Write('Cho bac k cua da thuc:'); readln(k1); if k1=1 then begin n:=4 ; write(' So cot cua bang so lieu:',n); t:=false; end; if k1=2 then begin n:=10; write('So cot cua bang so lieu:',n); t:=false; end; if k1=3 then begin n:=20; write('So cot cua bang so lieu:',n); t:=false; end; if (k11) and(k12) and(k13) then Người hướng dẫn: TS Vũ Huy Lân Học viên: Dương Văn Hội 77 LUẬN VĂN TỐT NGHIỆP begin t:=false; Write('So cot cua bang so lieu :');readln(n); end; end; writeln; for i:=1 to m for j:=1 to n begin write('a[',i,',',j,']='); readln(mt[i,j]); end; repeat begin writeln('Ban co muon sua bang so lieu vua nhap ko?'); writeln; Writeln('Co thi nhan y or Y, Ko thi nhan phim bat ki'); readln(ck); if ((ck='y')or(ck='Y')) then begin write(' Can sua hang:');readln(a1); write(' Can sua cot:' );readln(b1); write('a[',a1,',',b1,']='); readln(mt[a1,b1]); end; end; until (ck'y')and(ck'Y'); write(f,mt); close(f); assign(f,'matran.dat'); reset (f); read(f,mt); writeln('*********** Bang so lieu ************'); writeln; for i:=1 to m begin for j:=1 to n write(mt[i,j]:8:2); writeln; end; writeln; matranX_Y;readln; Người hướng dẫn: TS Vũ Huy Lân Học viên: Dương Văn Hội 78 LUẬN VĂN TỐT NGHIỆP writeln; matran_chuyenvi; readln; writeln; tinh_tich_hai_ma_tran(cv1,Y,d,d1,d2);readln; writeln; Gauss(d2,d1,x);readln; writeln; tinh_tich_ma_tran(d,x,yi);readln; writeln; he_so_tuong_quan_R; close(f); readln; end Phần mềm tối ưu hoá program toiuu; var i,j,k:integer; Y,x1,x2,x3,x1kq,x2kq,x3kq:real; Ykq:real; begin Y:=6; for i:=1 to 100 for j:=1 to 100 for k:=1 to 100 begin x1:=6+0.1*i; x2:=500+0.1*j; x3:=40+0.1*k; if (x1