Nghiên cứu ảnh hưởng của mn và si trong thành phần thuốc hàn thiêu kết hệ

BỘ GIÁO DỤC VÀ ĐÀO TẠO TRƢỜNG ĐẠI HỌC BÁCH KHOA HÀ NỘI - NGUYỄN CÔNG VƢỢNG NGHIÊN CỨU ẢNH HƢỞNG CỦA Mn VÀ Si TRONG THÀNH PHẦN THUỐC HÀN THIÊU KẾT HỆ BAZƠ TƢƠNG ĐƢƠNG LOẠI F7A(P)4 THEO AWS A5.17-80 ĐẾN CƠ TÍNH MỐI HÀN LUẬN VĂN THẠC SĨ KHOA HỌC CHUYÊN NGÀNH: CÔNG NGHỆ HÀN Hà Nội – Năm 2013 BỘ GIÁO DỤC VÀ ĐÀO TẠO TRƢỜNG ĐẠI HỌC BÁCH KHOA HÀ NỘI - NGUYỄN CÔNG VƢỢNG NGHIÊN CỨU ẢNH HƢỞNG CỦA Mn VÀ Si TRONG THÀNH PHẦN THUỐC HÀN THIÊU KẾT HỆ BAZƠ TƢƠNG ĐƢƠNG LOẠI F7A(P)4 THEO AWS A5.17-80 ĐẾN CƠ TÍNH MỐI HÀN Chuyên ngành : CÔNG NGHỆ HÀN LUẬN VĂN THẠC SĨ KHOA HỌC CHUYÊN NGÀNH: CÔNG NGHỆ HÀN NGƯỜI HƯỚNG DẪN KHOA HỌC : TS: Vũ Huy Lân Hà Nội – Năm 2013 LỜI CAM ĐOAN Tôi Nguyễn Công Vượng, học viên lớp Cao học Công nghệ hàn – Khoá 2011B, Trường Đại học Bách Khoa Hà Nội, thực đề tài: “Nghiên cứu ảnh hƣởng Mn Si thành phần thuốc hàn thiêu kết hệ bazơ tƣơng đƣơng loại F7A(P)4 theo AWS A5.17-80 đến tính mối hàn ” Tác giả tham gia Đề tài Khoa học Công nghệ Trọng điểm cấp Nhà nước mã số KC.02.04/11-15 TS Vũ Huy Lân làm Chủ nhiệm Đề tài Tác giả xin cam đoan rằng: Ngoại trừ số liệu, bảng biểu, đồ thị,… trích dẫn từ tài liệu tham khảo số liệu, nội dung lại công bố Luận văn tác giả nhóm tác giả tham gia Đề tài đưa Nếu sai xin hoàn toàn chịu trách nhiệm Hà Nội, ngày 28 tháng năm 2013 Học viên Nguyễn Công Vƣợng DANH MỤC CÁC KÝ HIỆU VÀ CHỮ VIẾT TẮT AC Alternating Curent AWS American Welding Society ASTM American Society for Testing of Materials ASME American Society of Machine Engineers DC Direct Current DCEN, DC - Direct Current Electrode Negative DCEP, DC + Direct Current Electrode Positive DT Destructive Testing EN European Standards HAZ Heat affected zone IIW International Institute Welding ISO International Standard Organization JIS Japanese Industrial Standards KLCB Kim loại NDT Nondestructive Testing PWHT Post Weld Heat Treatment SAW Submerged Arc Welding SMAW Shielded Metal Arc Welding TCVN Tiêu chuẩn Việt Nam VAHN Vùng ảnh hưởng nhiệt WPS Welding Procedure Specification WPQR Welding Procedure Qualiffication Recode ak Độ dai va đập σch Giới hạn chảy σk Độ bền kéo δ Độ dãn dài tương đối Ψmh Hệ số hình dạng mối hàn Ψn Hệ số ngấu mối hàn DANH MỤC CÁC BẢNG Bảng 1.1 Các nguồn nguyên liệu nước dự kiến sử dụng [1] 23 Bảng 2.1 Phân loại kí hiệu thuốc hàn theo IIW – 545 –78 31 Bảng 2.2 yêu cầu tính kim loại mối hàn theo AWS A5.17 – 80 34 Bảng 2.3 Công thử độ dai va đập kim loại mối hàn 34 Bảng 2.4 Thành phần hoá học số thông số chủ yếu dây hàn 34 tự động lớp thuốc theo AWS A5.17 – 80 34 Bảng 2.5.Các loại xỉ hàn theo thành phần hoá học chất theo IIW 46 Bảng 2.6: Thành phần sơ mẻ liệu thuốc hàn (mẻ liệu khô) 53 Bảng 2.7: Quy đổi thành phần mẻ liệu có tính đến nước thủy tinh 53 Bảng 3.1: Giá trị khoảng biến thiên yếu tố đầu vào 71 Bảng 3.2: Kế hoạch thực nghiệm bậc hai trực giao (với biến mã hóa) 72 Bảng 3.3: Thành phần hoá học dây hàn SM400B, (%): 79 Bảng 3.4: Thành phần hoá học dây hàn EM12, (%): 79 Bảng 3.5: Các tiêu tính dây hàn: 80 Bảng 3.6: Kích thước mẫu thử độ bền kéo 89 Bảng 4.1: Giá trị khoảng biến thiên yếu tố 92 Bảng 4.2: Kế hoạch thực nghiệm thí nghiệm trực giao bậc 92 Bảng 4.3: Các thông số chế độ hàn 93 Bảng 4.4: Kết thí nghiệm 95 Bảng 5.1 Hàm lượng Fe-Mn, Fe-Si đưa vào mẻ liệu thuốc hàn F7A4-1 dự kiến thành phần hóa học tính toán kim loại mối hàn 103 Bảng5.2 Kết thử nghiệm thuốc hàn F7A4-1 thành phần hóa học 103 kim loại mối hàn 103 Bảng 5.3 Hàm lượng Fe-Mn, Fe-Si đưa vào mẻ liệu thuốc hàn F7A4-2 dự kiến thành phần hóa học tính toán kim loại mối hàn 104 Bảng 5.4 Kết thử nghiệm thuốc hàn F7A4-2 thành phần hóa học 104 kim loại mối hàn (Mẫu thử theo ANSI/AWS B4.0) 104 DANH MỤC CÁC HÌNH VẼ VÀ ĐỒ THỊ Hình 1.1 Sơ đồ nguyên lý trình hàn tự động lớp thuốc 14 Hình 1.2 Một số thiết bị hàn hồ quang tự động lớp thuốc 14 Hình 1.3 Một số hình ảnh hàn hồ quang tự động lớp thuốc 15 Hình 1.4 Sơ đồ nguyên lý hàn hồ quang lớp thuốc điện cực băng 15 Hình 1.5 Hàn hồ quang lớp thuốc điện cực băng 16 Hình 2.1 Sơ đồ kí hiệu thuốc hàn – dây hàn AWS A5.17 - 80 33 Hình 2.2 Hàm lượng oxi kim loại đắp phụ thuộc vào hệ số hoạt tính 39 hoá học thuốc-xỉ hàn hàn chế độ hàn trung bình 39 thép bon thép hợp kim thấp 39 Hình 2.3 Độ sệt xỉ nhiệt độ 1400°C, poa 41 Hình 2.4 Sự thay đổi độ nhớt xỉ phụ thuộc vào nhiệt độ 43 1-xỉ dài; 2- xỉ ngắn 43 Hình 2.5 Giản đồ trạng thái nhiệt độ hệ xỉ hàn, °C 49 Hình 2.6 Độ sệt xỉ nhiệt độ 1600°C, Pa.s 50 Hình 2.7 Sức căng bề mặt hệ xỉ hàn 1600°C , J/m2 50 Hình 2.8 Sơ đồ nghiên cứu tổng quát thuốc hàn 54 Hình 2.9 Sơ đồ nghiên cứu đề tài 55 Hình 3.1 Ảnh hưởng hàm lượng ôxit mangan thuốc hàn 59 đến hàm lượng gia tăng S kim loại mối hàn 59 Hình 3.2 Hàm lượng S kim loại mối hàn phụ thuộc 59 vào tính bazơ thuốc hàn 59 Hình 3.3 Hàm lượng P kim loại mối hàn phụ thuộc vào 60 tính bazơ thuốc hàn; 1– (P)th = 0,03%; 2– (P)th = 0,05 60 Hình 3.4: Máy hàn tự động Dosun DZ1000 80 Hình 3.5: Kích thước mẫu thử phân tích thành phần hóa học kim loại mối hàn 82 Hình 3.6: Phôi chuẩn bị hàn mẫu thử phân tích thành phần hóa học kim loại 83 mối hàn 83 Hìn 3.7: Máy phân tích quang phổ 84 Hình 3.8: Các kích thước mầu hàn để chuẩn bị mẫu thử tính 86 Hình 3.9: Vị trí lấy mẫu thử độ bền kéo độ dai va đập 87 Hình 3.10 : Mẫu thử kéo kim loại mối hàn: 88 Hình 3.11: Cắt mẫu thử độ bền kéo độ dãn dài 88 Hình 3.12: Ví dụ hình dạng kích thước mẫu thử độ bền kéo độ dãn dài 88 Hình3.13: Máy kéo nén H011- Matest- Italia 89 Hình3.14: Hình dạng kích thước mẫu chuẩn đầy đủ thử dai va đập 90 theo tiêu chuẩn AWS B4.0 90 Hình 3.15: Mẫu chuẩn đầy đủ thử dai va đập theo tiêu chuẩn AWS (có dung sai) 90 Hình 3.16: Máy thử va đập vị trí búa đập 91 Hình 4.1: Các phôi mẫu sau hàn xong tiến hành cắt để lấy mẫu thử 94 thành phần hóa học kim loại mối hàn 94 Hình 4.2 : Mẫu hàn cắt để chuẩn bị phân tích thành phần hóa học 94 Hình 5.1 Sự phụ thuộc hàm lượng Mn kim loại mối hàn vào 98 hàm lượng Fe-Mn mẻ liệu thuốc hàn 98 Hình 5.2 Sự phụ thuộc hàm lượng Mn kim loại mối hàn vào hàm 98 lượng Fe-Si mẻ liệu thuốc hàn 99 Hình 5.3 Sự phụ thuộc hàm lượng Mn kim loại mối hàn vào 99 hàm lượng Fe-Mn Fe-Si mẻ liệu thuốc hàn 99 Hình 5.4 Sự phụ thuộc hàm lượng Si,(%) kim loại mối hàn vào 99 hàm lượng Fe-Si, (%) mẻ liệu thuốc hàn 99 Hình 5.5 Sự phụ thuộc hàm lượng Si kim loại mối hàn vào hàm lượng Fe-Mn mẻ liệu thuốc hàn 100 Hình 5.6 Sự phụ thuộc hàm lượng Si kim loại mối hàn vào hàm lượng Fe-Mn Fe-Si mẻ liệu thuốc hàn 100 Hình 5.7 Gá lắp phôi để hàn mẫu thử tính kim loại mối hàn 106 Hình 5.8 Hàn mẫu thử tính kim loại mối hàn 106 Hình 5.9 Quá trình hàn mẫu thử tính kim loại mối hàn 107 Hình 5.10 Mẫu sau thử kéokim loại mối hàn 107 Hình 5.11 Các mẫu sau thử độ dai va đậpkim loại mối hàn 108 MỞ ĐẦU Trong nghiệp công nghiệp hoá đất nước ngành khoa học công nghệ phát triển mạnh, ngành công nghệ hàn không ngừng phát triển bề rộng có công nghệ hàn đa dạng giải pháp nâng cao suất chất lượng hàn Do ngành hàn ngày ứng dụng rộng rãi, số công nghệ hàn tự động lớp thuốc ứng dụng phổ biến giới Ở Việt Nam công nghệ ngày ứng dụng rộng rãi, đặc biệt lĩnh vực đóng tàu chế tạo kết cấu thép, Tuy nhiên, thuốc hàn chủ yếu phải nhập từ nước ngoài, để khắc phục tình trạng số sở sản xuất vật liệu hàn lớn nước Công ty sản xuất que hàn Việt Đức, công ty sản xuất vật liệu hàn Nam Triệu, nghiên cứu chế tạo thử Nhưng việc ứng dụng hạn chế, chưa nghiên cứu đầy đủ tiêu cần thiết để nâng cao chất lượng thuốc hàn Trong số tiêu quan trọng này, hàm lượng Mn Si kim loại mối hàn giải pháp để hợp kim mối Nội dung cụ thể luận văn tốt nghiệp: “ Nghiên cứu ảnh hƣởng Mn Si thành phần thuốc hàn thiêu kết hệ bazơ tƣơng đƣơng loại F7A(P)4 theo AWS A5.17-80 đến tính mối hàn ” Có vấn đề cần giải sau: - Giới thiệu tổng quan việc nghiên cứu giới Việt Nam thuốc hàn thiêu kết Xác định hướng nghiên cứu đề tài - Các vấn đề cần giải quyết: + Phân tích, lựa chọn nguyên tố hợp kim phổ biến từ fero hợp kim dùng thành phần thuốc hàn chế tạo thuốc hàn thiêu kết tương đương với loại F7A(P)4 theo AWS A5.17-80 + Nghiên cứu ảnh hưởng fero hợp kim (Fe-Mn, Fe-Si) thành phần thuốc hàn đến thành phần hóa học tính kim loại mối hàn Nhận ét Các kết thử tính cho thấy hàm lượng Mn Si kim loại mối hàn đáp ứng tốt yêu cầu đề Các kết làm sở cho việc thiết kế mẻ liệu thuốc hàn thiêu kết mác F7A4 kết hợp với dây hàn EM12K 5.5 Giới thiệu thuốc hàn để sản uất thử nghiệm Với hàm lượng ferô xác định kết hợp với dây hàn EL8 tiến hành hàn thử nghiệm phân tích thành phần hóa học, thử tính kim loại mối hàn đáp ứng tốt tiêu chí yêu cầu Do vậy, sử dụng kết nghiên cứu trên, đưa hàm lượng ferô Fe-Mn Fe-Si vào mẻ liệu thuốc hàn F7A4 để sản xuất thử nghiệm sau: Bảng5.5 Hàm lượng Fe-Mn, Fe-Si đưa vào mẻ liệu thuốc hàn để sản xuất thử nghiệm thuốc hàn F7A2 Hàm lượng ferô hợp Hàm lượng nguyên tố kim loại mối hàn, % kim, % Fe-Mn Fe-Si C Mn Si S P 7,5 2,0 0,0874 1,0179 0,4886 0,0073 0,0239 Thành phần hóa học tính kim loại mối hàn hàn thuốc hàn F7A2 thử nghiệm đạt yêu cầu 110 KẾT LUẬN VÀ KIẾN NGHỊ Kết luận: Đề 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 toá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 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ế Trong phạm vi tác giả đề cập đến nội dung sau:  Tổng quan công nghệ hàn tự động lớp thuốc  Ảnh hưởng chế độ hàn đến thành phần hóa học ( Mn, Si) đến tính mối hàn  Nghiên cứu dịch chuyển Mn Si từ ferô Fe-Mn Fe-Si mẻ liệu thuốc hàn vào mối hàn  Xác định hàm lượng Fe-Mn Fe-Si cho mẻ liệu thuốc hàn đảm bảo thành phần hóa học tính kim loại mối hàn theo yêu cầu Kết chủ yếu đề tài gồm nội dung sau đây:  Đề 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 111 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 toá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 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 Kiến nghị:  Các kết cấu hàn phương pháp hàn tự động lớp thuốc nên hàn vị trí hàn  Để nhận mối hàn mong muốn thay đổi chế độ hàn cường độ dòng điện hàn, hiệu điện hàn, vận tốc hàn ta cần nghiên cứu thêm ảnh hưởng tầm với điện cực, kích thước điện cực góc nghiêng điện cực, cần nghiên cứu nhiều chế độ hàn khác có ảnh hưởng đến mối hàn SAW, từ có đánh giá toàn diện điều chỉnh trình hàn SAW hiệu 112 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] Bùi Minh Trí (1980), Các phương pháp tính toán tối ưu hóa T1, T2 [4] Quách Tuần Ngọc (1996), Giáo trình pascal, ĐHBK Hà Nội [5] - TCVN 3223 : 2000, Que hàn thép cacbon thép hợp kim thấp [6] - TCVN 3223 : 2000, Que hàn thép cacbon thép hợp kim thấp [7] - TCVN 3223 : 2000, Que hàn thép cacbon thép hợp kim thấp [8] - Петров Г Л Сварочные материалы Машиностроение, Ленинград 1972 [9] - Пoтaпова Н Н Сварочные материалы для дуговой сварки Машиностроение, Москва 1989 [10] - Cпецэлектрод Covered electrodes for manual arc welding, depositing and cutting (Catalô que hàn hãng Spetselektrode, 2000) [11] - AWS Welding Handbook, 9th Edition, 2001 [12]- AWS D1.1/D1.1M (2006), Structural Welding Code - Steel [13] - Esab Welding Consumables, Esab Welding Co., LTD., 2006 [14] - Huyndai Welding Consumables, Huyndai Welding Co., LTD., Korea 2004 [15] - Kobelco Welding Handbook , Kobe Steel LTD., Japan 2004 [16] - Lincoln Welding Handbook , Lincoln Welding Co., LTD., USA 2005 113 PHỤ LỤC Phụ lục Danh sách nguyên vật liệu để sản xuất thuốc hàn Việt Nam Phụ lục Các tiêu k thuật nguyên vật liệu sản xuất thuốc hàn STT Tên nguyên liệu Kí hiệu Thành phần hoá học Độ hạt Dây hàn H-08A Fe - Mn C :  0,1 Mn : 0,35 - 0,6 Si :  0,03 P :  0,03 S :  0,03 Mn : 78 - 85 C :1 Si :  1,5 P :  0,2 S :  0,03 SiO2 : 63 - 67 Al2O3 : 14 - 19 Fenspat Fe2O3 : 0,1 - 0,5 K2O : - 12 P :  0,03 S :  0,03 SiO2 : 50 - 56 Al2O3 : 28 - 35 Cao lanh Fe2O3 : max P :  0,03 S :  0,03 CaCO3 :  97 114 Đá vôi P :  0,03 S :  0,03 CaCO3 : 41 - 48 Đôlômit MgO3 : 43 - 76 SiO2 : 55 - 56 MgO : 28 - 31 Al2O3 : max Bột tan Fe2O3 : max P :  0,06 S :  0,03 Rutil 85 (Rutil nhân tạo) Rutil 10 TiO2 :  85 FeO :  3,5 P :  0,03 S :  0,03 TiO2 :  92 FeO :  1,8 P :  0,03 S :  0,03 TiO2 Ilmenit 11 : 50 - 52 FeO : 38 - 40 P :  0,03 S :  0,03 SiO2 : 45 - 55 AL2O3 : 29 - 36 Mica Fe2O3 : max 115 P :  0,03 S :  0,06 12 Ôxít Titan (Ti2O) TiO2 13 Si C Fe - Si 14 Mn : 98  99 : 42 - 45 :1 :  1,5 P :  0,2 S :  0,03 Tỷ trọng : 150g/1 Nước thủy tinh 116 Phụ lục Phần mềm 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 begin for j:=1 to n 117 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 begin d2[i,j]:=0; for j:=1 to m d2[i,k]:=d2[i,k]+a1[i,j]*b2[j,k]; end; 118 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]; b[i,l]:=b[m,l]; b[m,l]:=c; end; if m=i then begin t:=false; 119 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; 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 120 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; 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); 121 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 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 122 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; 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; 123 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