Ứng dụng phương pháp song song giải bài toán cơ học

ĐẠI HỌC QUỐC GIA TP HCM TRƯỜNG ĐẠI HỌC BÁCH KHOA BÙI HOÀNG GIANG ỨNG DỤNG PHƯƠNG PHÁP SONG SONG GIẢI BÀI TỐN CƠ HỌC Chun ngành :Tốn ứng dụng Mã số: 604636 LUẬN VĂN THẠC SĨ Thành phố HỒ CHÍ MINH - 2011 CƠNG TRÌNH ĐƯỢC HỒN THÀNH TẠI TRƯỜNG ĐẠI HỌC BÁCH KHOA –ĐHQG -HCM Cán hướng dẫn khoa học : TS Đặng Văn Vinh Cán chấm nhận xét :TS Nguyễn Bá Thi Cán chấm nhận xét : PGS.TSKH Bùi Tá Long Luận văn thạc sĩ bảo vệ Trường Đại học Bách Khoa, ĐHQG Tp HCM ngày tháng năm 2011 Thành phần Hội đồng đánh giá luận văn thạc sĩ gồm: (Ghi rõ họ, tên, học hàm, học vị Hội đồng chấm bảo vệ luận văn thạc sĩ) Xác nhận Chủ tịch Hội đồng đánh giá LV Trưởng Khoa quản lý chuyên ngành sau luận văn sửa chữa (nếu có) CHỦ TỊCH HỘI ĐỒNG TRƯỞNG KHOA………… ĐẠI HỌC QUỐC GIA TP.HCM TRƯỜNG ĐẠI HỌC BÁCH KHOA CỘNG HÒA XÃ HỘI CHỦ NGHĨA VIỆT NAM Độc lập - Tự - Hạnh phúc NHIỆM VỤ LUẬN VĂN THẠC SĨ Họ tên học viên: Bùi Hoàng Giang Ngày, tháng, năm sinh: Chuyên ngành: 13 – 05– 1986 Toán ứng dụng MSHV: 02408180 Nơi sinh: Kiên Giang Mã số : 604636 I TÊN ĐỀ TÀI: ỨNG DỤNG PHƯƠNG PHÁP SONG SONG GIẢI BÀI TOÁN CƠ HỌC II NHIỆM VỤ VÀ NỘI DUNG : Khảo sát phương pháp song song giải hệ phương trình tuyến tính áp dụng vào giải tốn học Xây dựng chương trình phần tử hữu hạn tính tốn học kiểm nghiệm kết III NGÀY GIAO NHIỆM VỤ: ngày 14 tháng năm 2011 IV NGÀY HOÀN THÀNH NHIỆM VỤ: ngày 15 tháng năm 2011 V CÁN BỘ HƯỚNG DẪN: TS Đặng Văn Vinh Tp HCM, ngày tháng năm 2011 CÁN BỘ HƯỚNG DẪN CHỦ NHIỆM BỘ MÔN ĐÀO TẠO (Họ tên chữ ký) (Họ tên chữ ký) TS Đặng Văn Vinh PGS.TS Nguyễn Đình Huy TRƯỞNG KHOA….……… (Họ tên chữ ký) Lời cảm ơn Trước hết xin cảm ơn TS Đặng Văn Vinh nhiệt tình hướng dẫn, tạo điều kiện, đóng góp nhiều ý kiến q báu để tơi hồn thành luận văn Tôi chân thành cảm ơn thầy giáo mơn Tốn Ứng Dụng nhiệt tình giảng dạy truyền đạt lại kiến thức kinh nghiệm chun mơn Tơi xin bày tỏ lịng biết ơn tới gia đình, bạn bè động viên khích lệ tơi hồn thành luận văn Thành phố Hồ Chí Minh, Ngày 11 tháng năm 2011 Học viên Bùi Hồng Giang Tóm tắt luận văn Luận văn gồm phần phụ lục Phần I phần lí thuyết tốn học thuật tốn, bao gồm chương: • Chương I trình bày qui tắc phân tích độ nhạy nghiệm toán, ảnh hưởng số điều kiện lên độ xác nghiệm hệ, qui tắc phân tích sai số tính tốn nói chung tính tốn máy tính nói riêng Chương I đề cập đến giải thuật LU để giải trực tiếp hệ phương trình tuyến tính hai cải tiến phương pháp LU partial pivoting & complete pivoting • Chương II trình bày cách biểu diễn ma trận thưa cách lưu trữ ma trận thưa máy tính Trong chương trình bày khái niệm thuật tốn phương pháp chiếu không gian Krylov làm tảng cho phương pháp lặp hiệu lực Conjugate Gradient, GMRES Lí thuyết thuật tốn đề cập • Chương III đề cập đến khái niệm & mơ hình tính tốn song song cluster cách thiết kế chương trình song song dùng MPI OpenMP • Chương IV trình bày lớp thuật tốn preconditioning ma trận ILU ứng dụng preconditioning cho phương pháp Conjugate Gradient GMRES Phần Ứng dụng bao gồm chương: • Chương V trình bày sở lí thuyết học mơi trường liên tục phương trình tính tốn biến dạng vật rắn Khung chương trình (framework) phương pháp phần tử hữu hạn dùng để rời rạc hóa phương trình lượng biến dạng mơ tả nhằm cung cấp sở lí thuyết cho chương trình tính • Chương VI mơ tả chương trình tính tốn kết tính tốn với kết cấu thơng dụng Phụ lục gồm mã chương trình tính phần tử hữu hạn Phụ lục trình bày mã chương trình Petsc solver Summary This thesis contains two parts and two appendices Part I mention the mathematical theory, including four chapters: • Chapter I mainly discuss about sensitivity analysis of the linear system of equations and the effect of round off errors on to the accuracies of the solution Apart from that, theory of computation error on computer is developed and presented Chapter I also mention the LU decomposition method to solve the linear system of equations directly and two variants which is partial pivoting LU and complete pivoting LU • Chapter II present the sparse matrix compression method and storage format in computer memory In this chapter, concept of projection method and Krylov subspace method is mentioned as backbone for renown iterative method "Conjugate Gradient" or "GMRES" The algorithm of these method were also regarded • Chapter III mention the concept of parallel computing model on cluster system and guideline to design program using MPI or OpenMP • Chapter IV present class of preconditioning technique like ILU and application of preconditioner to Conjugate Gradient method and GMRES method The second part is the application part comprises of two chapters • Chapter V is the theory of continuum mechanics and the governing equation of rigid body deformation The framework of Finite Element Method is developed as the basis to discretize the blance energy equation which serve as the mathematical kernel of the computation program • Chapter VI describe the finite element program to compute deformation of rigid body in linear sense and some results for popular cases Appendix I is core part of finite element program code Appendix is code for solver utilising Petsc framwork Lời cam đoan Tơi, Bùi Hồng Giang, xin cam đoan luận văn viết trình bày kinh nghiệm kiến thức cá nhân mà không chép sửa đổi tài liệu khác mà/hoặc khơng có đồng ý tác giả sử dụng tài liệu không rõ nguồn gốc Luận văn viết nhằm hồn tất học trình cao học ngành Toán Ứng dụng Đại học Bách khoa Tp HCM với chủ đích khơng sử dụng đề cập hình thức khác Thành phố Hồ Chí Minh, Ngày 11 tháng năm 2011 Học viên Bùi Hoàng Giang Lời mở đầu Mục tiêu nội dung nghiên cứu Mục tiêu luận văn bao gồm: Tìm hiểu sở lí thuyết tốn học, thuật tốn giải hệ phương trình tuyến tính với ma trận thưa Viết chương trình giải toán biến dạng đàn hồi kết cấu Áp dụng chương trình viết giải số kết cấu thơng dụng có Dùng Ansys để kiểm tra kết tính tốn Hệ phương trình tuyến tính đề cập luận văn có dạng Ax = b (1) Trong đó, A ∈ Mn [R], x ∈ Rn , b ∈ Rn Như biết, phương trình (1) kết phân tích sau nhiều vấn đề lí thuyết thực tiễn sau phép rời rạc hóa lắp ghép ma trận độ cứng phần tử (VD: FEM, BEM, ) Đặc biệt tốn học tính biến dạng đàn hồi kết cấu Việc tìm cách giải hiệu thời gian nhớ điều cần thiết, nhiên khó để tìm phương pháp chung cho tất kiểu tốn thực tế khơng thể tìm phương pháp Sau nhiều nỗ lực nhà toán học, số phương pháp chứng nhận hiệu đề cho kiểu ma trận riêng ma trận thưa, ma trận khối có cấu trúc đặc biệt, Việc tìm phương pháp hiệu cho toán cần khảo sát đặc thù nhiều tính tốn thử nghiệm Trong luận văn này, tác giả khảo sát phương pháp hiệu lực sẵn có dành cho hệ phương trình tuyến tính với ma trận thưa tìm hiểu cách thức giải song song cluster dựa phương pháp Cụ thể hơn, tác giả giải tốn đặc thù tốn tính biến dạng đàn hồi tuyến tính kết cấu sử dụng phương pháp lặp (Conjugate Gradient & GMRES) sử dụng phương pháp trực tiếp để so sánh kết Dựa kết thu được, toán tính lại Ansys để kiểm nghiệm độ xác Ansys chọn phần mềm chứng nhận giới hàn lâm kĩ thuật Tình hình phạm vi nghiên cứu Tính đến thời điểm nay, phương pháp cổ điển nghiên cứu đầy đủ cơng trình tác giả lớn (như Golub & van Loan [6], Ke Chen [4]) Rất nhiều code phát triển để giải hệ phương trình tuyến tính theo phương pháp trực tiếp Lapack, UmfPack có hiệu cao với ma trận có kích cỡ trung bình Các code tinh chỉnh số hệ thống phần cứng cụ thể để đạt hiệu vượt trội Các phương pháp lặp nghiên cứu sâu (chẳng hạn Saad [1]), với phương pháp lặp (CG, BiCG, ) preconditioning ma trận, tạo nhiều công cụ cho kĩ sư nhà khoa học thực tính tốn Hơn nữa, ma trận thưa có cấu trúc đặc biệt (ba đường chéo) có tính chất đặc biệt (đối xứng xác định dương) tồn phương pháp hiệu để giải Mặt khác, có cách chọn phương pháp hiệu thống kê tối ưu hóa ma trận thưa có cấu trúc đặc thù gần suy biến Tuy nhiên chưa có phương pháp thực hiệu cho loại toán Những đóng góp luận văn Luận văn thực nghiên cứu giải thuật lặp có sẵn cách song song hóa thuật tốn cluster để giải tốn tính biến dạng đàn hồi Một thực tế chưa có nhiều nghiên cứu áp dụng tính tốn song song cho tốn dạng Các phần mềm có sẵn Ansys hỗ trợ tính tốn song song cho hệ thống đặc thù mà khơng có hỗ trợ khả chuyển rộng rãi Điểm mạnh luận văn thiết lập gói phần mềm tính toán phần tử hữu hạn cho toán học C++ giải song song toán tương ứng với 310464 phương trình 310464 ẩn Lưu ý với tinh vi phức tạp phương pháp phần tử hữu hạn, khơng có nhiều code mã nguồn mở hỗ trợ Mục đích tác giả viết code nhằm phục vụ nghiên cứu có khả mở rộng cho loại tốn khác Hơn tính khả chuyển dễ hiểu cân nhắc nhằm hướng tới mục tiêu mã nguồn mở tương lai Mục lục Danh sách kí hiệu I Phần lí thuyết Hệ phương trình tuyến tính phương pháp giải trực tiếp 1.1 Phân tích độ nhạy nghiệm 1.2 Phân tích sai số thuật toán 1.2.1 Phép tính dấu chấm động 1.2.2 Sai số làm tròn phép tính 1.3 Phân tích LU 1.3.1 Phân tích LU có pivot 1.3.2 Sai số thuật toán LU Ma trận thưa phương pháp giải lặp 2.1 Biểu diễn đồ thị ma trận thưa 2.2 Hoán vị lại 2.2.1 Các phép lại điển hình 2.3 Lưu trữ ma trận thưa 2.4 Toán tử chiếu 2.4.1 Biểu diễn ma trận toán tử chiếu 2.4.2 Toán tử chiếu vng góc 2.4.3 Tính chất tốn tử chiếu 2.5 Trực giao hóa Householder 2.6 Phương pháp chiếu 2.6.1 Một mơ tả hình học phương pháp chiếu 2.6.2 Phương pháp đường dốc (steepest descent) 2.6.3 Phương pháp số dư nhỏ 2.6.4 Phương pháp đường dốc cho số dư 2.7 Phương pháp không gian Krylov 2.7.1 Một số tính chất không gian Krylov 2.7.2 Xây dựng sở cho không gian Krylov 2.7.3 Phương pháp GMRES 2.7.4 Phương pháp Conjugate Gradient 9 11 11 11 13 15 17 19 20 20 22 22 24 26 26 27 28 30 34 35 35 36 36 37 38 40 42 Phụ lục A 541 FemElementLE.cpp T( ) = N( , ) ; T( ) = N( , ) ; T( ) = N( , ) ; // c o u t