Nghiên cứu mô phỏng và tính toán trường nhiệt độ khi hàn tấm mỏng

BỘ GIÁO DỤC VÀ ĐÀO TẠO TRƯỜNG ĐẠI HỌC BÁCH KHOA HÀ NỘI - NGÔ THỊ THẢO NGHIÊN CỨU, MÔ PHỎNG VÀ TÍNH TOÁN TRƯỜNG NHIỆT ĐỘ KHI HÀN TẤM MỎNG 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 NGUYỄN TIẾN DƯƠNG Hà Nội - Năm 2011 LỜI CAM ĐOAN Tôi xin cam đoan số liệu kết nêu Luận văn trung thực chưa công bố công trình khác Trừ phần tham khảo nêu rõ Luận văn Tác giả Ngô Thị Thảo LỜI CẢM ƠN Tác giả xin chân thành cảm ơn TS Nguyễn Tiến Dương, người hướng dẫn giúp đỡ tận tình từ định hướng đề tài, tổ chức thực nghiệm đến trình viết hoàn chỉnh Luận văn Tác giả bày tỏ lòng biết ơn thầy cô môn Hàn Công nghệ kim loại - viện Cơ khí trường Đại hoc Bách Khoa Hà Nội Xin cảm ơn Ban lãnh đạo Viện đào tạo Sau đại học, Viện Cơ khí trường Đại học Bách Khoa Hà Nội tạo điều kiện thuận lợi để hoàn thành Luận văn Tác giả chân thành cảm ơn môn Công nghệ hàn Khoa Cơ khí trường ĐHSPKT Hưng Yên tạo điều kiện giúp đỡ tác giả thực đề tài Do lực thân nhiều hạn chế nên Luận văn không tránh khỏi sai sót, tác giả mong nhận đóng góp ý kiến Thầy, Cô giáo, nhà khoa học bạn đồng nghiệp Tác giả Ngô Thị Thảo DANH MỤC BẢNG Bảng 5.1 -Bảng kết thực nghiệm Bảng 5.2 –So sánh kết DANH MỤC HÌNH VẼ Hình 2.1 Hàn hồ quang tay(SMAW) Hình 2.2-Hàn TIG Hình 2.3-Hàn MAG Hình 2.4-Hàn lớp thuốc Hình 2.5-Hàn điện xỉ Hình 2.6 - Nguồn nhiệt điểm tức thời vật thể bán vô hạn Hình 2.7- Nguồn nhiệt đường tức thời vật thể bán vô hạn Hình 2.8- Nguồn nhiệt phẳng tức thời vật thể bán vô hạn Hình 3.1- Nguồn nhiệt đường mỏng Hình 3.2 Nguồn nhiệt đường di chuyển mỏng Hình 3.3- Đồ thị hàm Bessel K (u ) K1 (u ) Hình 3.4- Đồ thị lời giải mỏng Rosenthal Hình 5.1- Chuẩn bị phôi Hình 5.2- Đo chiều rộng mối hàn Hình 5.3- Đo chiều rộng đường bao Hình 5.4- Đo chiều rộng vùng ảnh hưởng nhiệt Hình 5.5- Đo nhiệt độ hàn PHẦN MỞ ĐẦU * Lý lựa chọn đề tài, lĩnh vực nghiên cứu - Hầu hết trình hàn liên quan đến nhiệt biến thiên nhiệt độ vật hàn Việc nghiên cứu trường nhiệt độ hàn đưa nghiệm khả thi cho toán hàn khác vấn đề quan tâm nhiều năm qua - Trường nhiệt độ hàn có ảnh hưởng định đến đặc tính liên kết hàn như: khả xuất ứng suất dư biến dạng dư hàn, hình dạng kích thước mối hàn, cấu trúc kim loại mối hàn vùng ảnh hưởng nhiệt… Những yếu tố ảnh hưởng đến chất lượng kết cấu hàn khả vận hành chúng điều kiện thực tế - Việc nghiên cứu vấn đề tiêu biểu công nghệ hàn ứng suất biến dạng, mô hình hóa cấu trúc tế vi kim loại mối hàn… đòi hỏi phải có thông tin chi tiết trình biến thiên nhiệt độ phần tử kêt cấu hàn Các nguồn nhiệt đóng vai trò chủ đạo phần lớn trình hàn, chi tiết nung nguồn nhiệt cục dòng nhiệt có xu hướng chuyển động từ nơi có nhiệt độ cao tới nơi có nhiệt độ thấp chi tiết Sự thay đổi nhiệt độ vị trí vật thể phụ thuộc vào tốc độ nguội vị trí *Mục đích nghiên cứu - Mô tính toán hình dạng kích thước vũng hàn, trường nhiệt độ hàn hồ quang liên kết giáp mối - Tính toán thể tích vũng hàn, tiết diện ngang đường hàn, chiều rộng vùng ảnh hưởng nhiệt - Tính toán thời gian nguội tốc độ nguội khoảng nhiệt độ từ 8000C đến 5000C * Lịch sử nghiên cứu Việc nghiên cứu vấn đề liên quan đến hàn ứng suất dư, biến dạng dư cấu trúc tế vi yêu cầu thong tin chi tiết lịch sử nhiệt phần tử kết cấu hàn Việc nghiên cứu thực nhiều năm qua với mục đích mô tính toán hình dạng kích thước vũng hàn, tiết diện ngang đường hàn, chiều rộng vùng ảnh hưởng nhiệt tính toán thời gian nguội, tốc độ nguội hàn * Tóm tắt luận văn đóng góp Với yêu cầu đòi hỏi tạo sản phẩm hàn có chất lượng đảm bảo, điều quan trọng cần nghiên cứu yếu tố ảnh hưởng đến chất lượng sản phẩm Trường nhiệt độ hàn định đến khả hình thành ứng suất dư, biến dạng dư, kích thước, hình dạng mối hàn Từ lựa chọn phương pháp hàn, chế độ hàn hợp lý cho vật liệu hàn cụ nhằm đáp ứng yêu cầu chất lượng kinh tế Ngoài việc nghiên cứu vấn đề tiêu biểu công nghệ hàn ứng suất biến dạng, mô hình hóa cấu trúc tế vi kim loại mối hàn… đòi hỏi phải có thông tin chi tiết trình biến thiên nhiệt độ phần tử kêt cấu hàn Các nguồn nhiệt đóng vai trò chủ đạo phần lớn trình hàn.Sự thay đổi nhiệt độ vị trí vật thể phụ thuộc vào tốc độ nguội vị trí Từ vấn đề đó, nội dung đề tài nghiên cứu Nội dung gồm chương Chương 1: Trình bày tổng quan phương trình truyền nhiệt Tác giả đưa khái niệm truyền nhiệt, công thức tính trường nhiệt cho trường hợp điển hình: vật thể khối, vật thể bán vô hạn, phẳng… Chương 2: Khái quát trình hàn nguồn nhiệt hàn Ở chương tác giả giới thiệu số trình hàn nóng chảy dùng phổ biến đưa số nguồn nhiệt hàn thường gặp Chương 3: Trường nhiệt hàn mỏng Tác giả giới thiệu mô hình tính toán, phân bố nhiệt công thức tính toán trạng thái giả ổn định Chương 4: Tính toán Tác giả lựa chọn trình hàn với chế độ hàn xác định cho vật liệu hàn cụ thể Sử dụng phần mềm Excel để tính toán Chương 5: Thực nghiệm Tác giả tiến hành thí nghiệm mẫu thật với thông số phần tính toán lý thuyết, sau lấy kết so sánh với lý thuyết Chương 6: Mô trường nhiệt độ Tác giả sử dụng phần mềm Ansys để mô trường nhiệt độ hàn mỏng Phần kết luận tác giả trình bày tóm tắt toàn nội dung luận văn kiến nghị Phần phụ lục tác giả trình bày chương trình tính toán Excel chương trình mô trường nhiệt độ Ansy * Phương pháp nghiên cứu - Thu thập phương pháp tính toán mối hàn sách khoa học; - Tiến hành nghiên cứu thực nghiệm để đo kích thước vũng hàn, thể tích vũng hàn, chiều rộng vùng ảnh hưởng nhiệt, thời gian nguội khoảng nhiệt độ từ 8000C đến 5000C - Xây dựng bảng biểu, đồ thị Dùng phần mềm mô trường nhiệt độ hàn - Xây dựng công thức tính kích thước, thể tích vũng hàn, chiều rộng vùng ảnh hưởng nhiệt tốc độ nguội * Ý nghĩa khoa học đề tài Việc nghiên cứu, tính toán trường nhiệt độ giúp ta tính toán kiểm soát nhiệt độ trình hàn, từ có biện pháp cải thiện, nâng suất chất lượng hàn * Hướng giải quyết, nơi thực đề tài nghiên cứu luận án - Hướng giải quyết: mặt lý luận dựa theo lý truyền nhiệt hàn ứng dụng - Về mặt thí nghiệm : tiến hành thí nghiệm hàn hồ quang xưởng thực tập trường ĐHSPKT Hưng Yên + Vật liệu thép cacbon thấp: CT38 + Phương pháp hàn: GMAW + Vị trí hàn : 1G + Thí nghiệm hàn đắp hàn giáp mối tấm, tiến hành đo vùng ảnh hưởng nhiệt dựa thay đổi màu sắc vùng mẫu, đo kích thước lượt hàn, cắt mẫu dùng hoá chất tẩm thực để đo vùng nóng chảy kim loại lớp mối hàn, đo thời gian nguội từ 8000C đến 5000C Chương I: TỔNG QUAN VỀ CÁC PHƯƠNG TRÌNH TRUYỀN NHIỆT 1.1 Khái niệm truyền nhiệt Hầu hết phương pháp hàn dùng nguồn nhiệt tập trung vào vị trí mối hàn đến nhiệt độ biến dạng dẻo Tuỳ thuộc vào tích chất vật liệu hàn kết cấu mà truyền nhiệt nhanh hay chậm khác vùng xung quanh Chế độ nung nóng làm nguội không nguyên nhân làm xuất khuyết tật hàn (nứt nóng, nứt nguội, mối hàn không ngấu, làm liên kết cong vênh, tạo ứng suất dư) làm giảm khả làm việc kết cấu Trạng thái nhiệt kim loại hàn, xỉ hàn phần tử khác tương tác trình hình thành mối hàn, trình hoá lý, trình luyện kim, làm thay đổi tổ chức liên kết hàn Đặc biệt hàn nóng chảy trình luyện kim xảy ngắn ảnh hưởng tới tổ chức tính chất liên kết hàn Những giả thuyết đặc trưng cho ứng suất biến dạng phát sinh kết cấu hàn chủ yếu phụ thuộc vào chu trình nung nóng làm nguội, phụ thuộc lan truyền nhiệt liên kết, tổ chức kim loại mối hàn Kim loại gần mối hàn, vùng ảnh hưởng nhiệt xác định đặc điểm phân bố lan truyền nhiệt, nguội tốc độ lan truyền nhiệt định vùng ảnh hưởng nhiệt Nếu không tính toán đến trạng thái nhiệt kim loại không giải thích sâu sắc đầy đủ tượng xuất hàn Muốn nghiên cứu trình hàn phải nghiên cứu đầy đủ trình xảy ra, phải đặt giả thuyết gần định luật nung nóng vật thể, kết cấu lan truyền nhiệt kết cấu Muốn nghiên cứu cách tỉ mỉ phải hiểu đại lượng vật lý vật liệu kết cấu • Số lượng nhiệt nguồn sinh : Q (cal) • Nhiệt dung khối lượng riêng : C (cal/g 0C) • Nhiệt dung thể tích : C γ ( cal /cm3.0C) Với chất rắn nhiệt dung thay đổi phụ thuộc vào nhiệt độ, áp suất, muốn nung nóng , vật thể lên 10C khoảng nhiệt độ hỗn hợp tiêu hao 1lượng nhiệt khác Vật thể nhân khối lượng nhiệt Q để nung nóng khối lượng nhiệt G từ T10C Æ T20C xác định giá trị trung bình nhiệt dung khoảng nhiệt độ cho trước C= Q (cal / g.0 C ) G (T2 − T1 ) - Nhiệt dung nhiệt độ T : CT = dQ (cal / g.0 C ) GdT - Số lượng nhiệt để nung nóng 1g vật chất từ 00C đên T0C cho trước S (cal/g) S =CT Æ C = S T - Sự phân bố nhiệt bề mặt hay thể tích kết cấu, đặc trưng môi trường nhiệt độ ¾ Trường nhiệt độ : tập hợp nhiệt độ phân bố theo thời gian tất điểm vật thể hay bề mặt mà có nhiệt độ Trường nhiệt độ đặc trưng đường đẳng nhiệt mặt đẳng nhiệt - Đường đẳng nhiệt: đường nhiệt độ điểm - Mặt đẳng nhiệt : mặt hình học có nhiệt độ điểm bề mặt Giao mặt đẳng nhiệt mặt thuộc vật thể ta có đường đẳng nhiệt - Biểu diễn trường nhiệt độ hệ phương trình hệ tọa độ không gian T = T(x,y,z,t) = f (x,y,z,t) với hệ tọa độ trụ T =T (r,ϕ,z) - Sự thay đổi nhiệt độ hệ phương trình hai đường đẳng nhiệt T1 T2 tính có LIM ( - T1 − T1 Khi hai đường đẳng nhiệt gần ∆s → O ta ∆S ∂T T1 − T1 ) ∆s → O = : Gradien nhiệt độ theo phương Sz : ∆S ∂S Gradien nhiệt độ theo phương tiếp tuyến với đường đẳng nhiệt =0, theo phương pháp tuyến có giá trị cực đại - Gradien nhiệt độ điểm cho trước có véctơ phương biến đổi nhiệt độ lớn trùng với pháp tuyến mặt đẳng nhiệt - Khi T0 tăng gradien > ngược lại Nguồn nhiệt tác dụng tâm t2 < t1→ gradien < - Nếu vật hàn có phân bố nhiệt độ không phần nhiệt độ cao truyền sang phần nhiệt độ thấp đến lúc điểm xét, nhiệt độ truyền qua có cân nhiệt - Trong vật rắn nhiệt độ truyền từ nơi cao đến nơi thấp gọi truyền nhiệt nội vật thể - Nhiệt truyền từ nội kết cấu hàn vật bên môi trường lỏng, khí gọi nhiệt truyền đối lưu xạ • Trong trình truyền nhiệt xảy ba điều kiện biên : Quá trình đoạn nhiệt: nhiệt độ bề mặt vật thể phụ thuộc vào tọa độ bề mặt thời gian, nhiệt độ phân bố bề mặt mà không phân bố bên vật thể, tức nhiệt độ bề mặt không thay đổi suốt thời gian t TF = T( F ) ( x, y , z , t ) + Dòng nhiệt phân bố qua bề mặt vật thể thuộc tọa độ bề mặt thời gian q F = q( F ) ( x, y , z , t ) KẾT LUẬN VÀ KIẾN NGHỊ Trong trình nghiên cứu thực đề tài tác giả tính toán thể tích vũng hàn, tiết diện ngang vũng hàn, tốc độ nguội thời gian nguội từ 8000C đến 5000C Sử dụng thiết bị sẵn có sở tiến hành hàn mối hàn mỏng, thu kết so sánh với kết tính toán lý thuyêt Áp dụng phần mềm Ansys để mô trường nhiệt độ hàn Nếu có điều kiện tác giả nghiên cứu mô trường nhiệt độ khối 3D số phần mềm khác 76 TÀI KIỆU THAM KHẢO [1]- Nguyễn Thế Ninh, Phân tích truyền nhiệt hàn ứng dụng, Nhà xuất Bách Khoa Hà Nội, 2008 [2]- Ngô Lê Thông, Công nghệ hàn điện nóng chảy (1&2), Nhà xuất Khoa học kỹ thuật Hà Nội, 2004 [3]- Đinh Bá Trụ, Hướng dẫn sử dụng phần mềm Ansys (1&2), Học viện kỹ thuật quân sự, 2003 [4]- Nghiêm Hùng, Vật liệu học sở, Nhà xuất Khoa học kỹ thuật, 2007 [5]- AWS D1.1/D1.1:2010 [6]- N.T Nguyen, Thermal analysis of welds, WIT Pres, Southhamton, 2004 77 PHỤ LỤC Phụ lục 1: Chương trình tính toán kích thước theo lý thuyết Excel 78 79 80 Phụ lục 2: Chương trình mô trường nhiệt độ hàn mỏng phần mếm Ansys !================ /Filname,giapmoi !Ten file /Units,SI !He don vi chuan *AFUN,DEG !Don vi goc la "do" /view,,-1,2,2 !Diem nhin 3D /Title,Luan van tot nghiep_Ngo Thi Thao_CNH2009 !===============Dinh nghia cac thong so: dai=250e-3 ![m] ! -Cac thong so qua trinh han: U=20 ![V]-Dien ap han I=130 ![A]-Dong dien han Vh=3.5e-3 nuy=0.8 ![m/s]-Toc han !He so nhiet hieu dung q=U*I*nuy/Vh ![W/m]-Nang luong duong cua Nguon nhiet han Lt=8e-3 inc=Lt/Vh conv1=1000 conv=20 Tstop=500 ![m]-Vet nguon nhiet ![s]-Buoc thoi gian ![W/(m2.K)]-truyen nhiet vao ga ![W/(m2.K)]-doi luu khong ![oC]-nhiet Stop tinh toan ! -Cac doan chia luoi: d1=1 d2=2 d3=10 d4=25 !==============Pre-processor: /PREP7 !====Xay dung cac diem keypoint: 81 K,1 K,2,-120e-3 K,3,-120e-3,3e-3 K,4,,3e-3 K,5,120e-3,3e-3 K,6,120e-3, K,7,-4e-3,3e-3 K,8,4e-3,3e-3 K,9,,5e-3 ! -Ve cac Line, Arc: L,1,2 L,2,3 L,3,7 L,7,4 L,4,1 L,4,8 L,8,5 L,5,6 L,6,1 Spline,7,9,8 ! -Ve cac AREA: A,1,2,3,7,4 !A1 A,1,4,8,5,6 !A2 AL,4,6,11,10 !A3-moi han ! -Dinh doan chia luoi cho cac AREA: Lesize,4,,,d2 Lesize,6,,,d2 Lesize,10,,,d2 Lesize,11,,,d2 82 Lesize,5,,,d1 Lesize,2,,,d1 Lesize,8,,,d1 Lesize,9,,,d3,0.5 !chia khong deu Lesize,7,,,d3,2 !chia khong deu Lesize,3,,,d3,0.5 !chia khong deu Lesize,1,,,d3,2 !chia khong deu ! -Keo dai cac Areas VOLUME: Vext,all,,,,,dai !Keo dai Area Volume ! -Dinh doan chia luoi cho cac canh cua Volume: Lesize,18,,,d3 Lesize,19,,,d3 Lesize,27,,,d3 Lesize,28,,,d3 Lesize,20,,,d4 Lesize,17,,,d4 Lesize,21,,,d4 Lesize,31,,,d4 Lesize,26,,,d4 !============Cac thuoc tinh phi tuyen cua Steel: !===1-Elastic Modul -[Pa]- Modun dan hoi MPTEMP,,,,,,, !xoa bang nhiet MPTEMP,1,20,200,400,600,800,1000 !bang nhiet moi MPTEMP,7,1100,1200,1300,1500 MPDATA,EX,1,1,210000E+6,200000E+6,175000E+6,135000E+6,78000E+6,1500 0E+6 MPDATA,EX,1,7,7000E+6,3000E+6,1000E+6,1000E+6 ! -!===2-Poisson - He so poat xong 83 MPTEMP,,,,,,, !xoa bang nhiet MPTEMP,1,20,1500 !bang nhiet moi MPDATA,PRXY,1,1,0.33,0.33 !bang gia tri moi ! -!===3-DEnsity -[kg/m3]- Ty = f(temp) MPTEMP,,,,,,, !xoa bang nhiet MPTEMP,1,20,1400 !bang nhiet moi MPDATA,DENS,1,1,7800,7300 !bang gia tri moi ! -!===4-Specific Heat -[J/kg.K]- Nhiet dung rieng MPTEMP,,,,,,, !xoa bang nhiet MPTEMP,1,20,100,300,500,800,900 !bang nhiet moi MPTEMP,7,1000,1500 MPDATA,C,1,1,400,500,600,650,850,650 !bang gia tri moi MPDATA,C,1,7,660,680 ! -!===5-Thermal conductivity -[W/m.K]- He so dan nhiet MPTEMP,,,,,,, !xoa bang nhiet MPTEMP,1,20,100,200,300,400,500 !bang nhiet moi MPTEMP,7,600,700,800,1500 MPDATA,KXX,1,1,46,46,45,43,40,37 !bang gia tri moi MPDATA,KXX,1,7,33,28,24,32 ! -!===6-Thermal Expansion -[%]- He so dan no nhiet MPTEMP,,,,,,, MPTEMP,1,20,1300 !xoa bang nhiet !bang nhiet moi MPDATA,ALPX,1,1,0,20e-6 !bang gia tri moi ! -!============Chon element va chia luoi cho cac VOLUME: 84 ET,1,Solid98 !Loai phan tu Solid 98 VMESH,ALL !Chia luoi cho tat ca cac Volume FINISH !Ket thuc Pre-processor !=====================Main Processor: /SOLU !Solution Antype,Trans !Transient solution Rescontrol,,last,last !Cho phep chay noi tiep - Restart solution NROPT,full !full Newton-Raphson method ESTIF,1e-10 !Gan gia tri cho Matrix se Kill OutRES,nsol,last Solcontrol,on !Dieu khien viec ghi du lieu tinh toan !Dieu khien tu dong viec noi suy Pstres,on !Ung suat truoc > Pre-Stress ! -Dieu kien ban dau: Toffset,273 !Lay 273 K lam C TREF,20 !Ban dau = 20 C (TUNIF=20 C) !==============Dieu kien bien ga kep han: Asel,s,area,,6 !Chon mat so 12 de dat dieu kien bien Asel,a,area,,13 !Chon them mat so 15 Nsla,s,1 !Chon cac nut thuoc mat vua chon o tren D,all,Ux,0 !Dieu kien bien Lien ket theo phuong X D,all,Uy,0 D,all,Uz,0 Allsel !Chon tat ca mo hinh !===========Face film - Truyen nhiet vao Do ga: Asel,s,area,,6 Asel,a,area,,13 Nsla,s,1 SF,all,CONV,conv1,20 ![W/(m2.K)],20[oC] !===========Face film - Doi luu Khong khi: 85 Asel,s,area,,7 Asel,a,area,,12 Asel,a,area,,14 Asel,a,area,,5 Asel,a,area,,4 Asel,a,area,,10 Asel,a,area,,1 Asel,a,area,,2 Nsla,s,1 SF,all,CONV,conv,20 ![W/(m2.K)],20[oC] Allsel !===============Tieu diet het cac phan tu thuoc moi han: Vsel,s,volu,,3 !Chon moi han NSLV,s,0 !Chon nut nam o phia cua khoi EslN,s,0,all !Chon Element gan voi cac nut vua chon EKILL,all !Kill tat ca cac Element vua chon !############################################### !=============Tien hanh han:================= !############################################### !===============Xac dinh vung tai sinh elements: tim=inc !End time of the first Load step z2=Vh*tim z1=z2-Lt !Xac dinh diem dat tai !Chieu dai dat tai !===============Tai sinh cac phan tu han: Vsel,s,volu,,3 !Chon khoi so - moi han NSLV,s,0 !Chon nut nam o phia cua khoi Nsel,r,loc,z,z1,z2 !Chon lai cac nut se duoc sinh Han EslN,s,0,all EALIVE,all !Chon Element gan voi cac nut vua chon !Tai sinh cac Element vua chon 86 !===============Dat tai nguon nhiet vao cac nut tai sinh: Vsel,s,volu,,3 !Chon khoi so - moi han NSLV,s,1 !Chon nut thuoc khoi (o + tren bien cua khoi) Nsel,r,loc,z,z1,z2 !Chon lai cac nut se duoc sinh Han *GET,sonut,NODE,,COUNT !Xac dinh so nut vua chon F,all,HEAT,q*Lt/sonut !Dat tai nguon nhiet han CHO TUNG NUT !===============Giai: Time,tim !Time of transient solution AutoTS,on !Tu dong tim buoc noi suy Allsel !Chon tat ca Model Solve !Giai he PT ma tran !===============Xoa tai nguon nhiet de chuyen sang vi tri moi: Vsel,s,volu,,3 !Chon khoi so - moi han NSLV,s,1 !Chon nut thuoc khoi (o + tren bien cua khoi) Nsel,r,loc,z,z1,z2 !Chon nut se duoc sinh Han Fdele,all,HEAT !Xoa tai nguon nhiet han ! -Luu ket qua: Save !Save to restart solution FINISH !Ket thuc /Solution !==============================Hien thi ket qua: /POST1 !Xuat ket qua /CVal,,600,727,850,911,1200,1300,1500,20000 !Thang nhiet