Giáo án Phương pháp phần tử hữu hạn

ử dụng lưới tương tự, trong khi Courant chia những miền liên tục thành những miền có hình tam giác cho cách giải thứ hai của phương trình vi phân từng phần elliptic, xuất hiện từ các bài toán về xoắn của phần tử thanh hình trụ. Sự đóng góp của Courant là phát triển, thu hút một số người nhanh chóng đưa ra kết quả cho PPVPTP elliptic được phát triển bởi Rayleigh, Ritz, và Galerkin. Sự phát triển chính thức của PPPTHH được bắt đầu vào nửa sau những năm 1950 trong việc phân tích kết cấu khung máy bay và công trình xây dựng, và đã thu được nhiều kết quả ở Berkeley (xem Early Finite Element Research at Berkeley) trong những năm 1960 trong ngành xây dựng. Phương pháp này được cung cấp nền tảng toán học chặt chẽ vào năm 1973 với việc xuất bản cuốn Strang và tổng kết Phương pháp phần tử hữu hạn 2 trong An Analysis of The Finite element Method và kể từ đó PPPTHH được tổng quát hóa thành một ngành của toán ứng dụng, một mô hình số học cho các hệ thống tự nhiên, được ứng dụng rộng rãi trong kĩ thuật, ví dụ như điện từ học và động lực học chất lỏng. Sự phát triển của PPPTHH trong cơ học kết cấu đặt cơ sở cho nguyên lý năng lượng, ví dụ như: nguyên lý công khả dĩ, PPPTHH cung cấp một cơ sở tổng quát mang tính trực quan theo quy luật tự nhiên, đó là một yêu cầu lớn đối với những kỹ sư kết cấu. Bài toán minh họa Chúng ta sẽ minh họa việc sử dụng PPPTHH từ hai ví dụ mà phương pháp chung có thể là ngoại suy. Chúng ta xem như người đọc đã quen thuộc với tính toán và đại số tuyến tính. Chúng ta sẽ sử dụng bài toán một chiều, tại đây, hàm f được xác định bởi u và u một hàm ẩn của x, u’’ là đạo hàm cấp 2 của u theo x Ví dụ cho bài toán hai chiều là bài toán Dirichlet Ở đây, miền Ω là một miền đơn liên mở trong mặt phẳng (x,y), có biên ∂Ω rất “đẹp” (ví dụ: một đa tạp trơn hoặc một đa giác), và là đạo hàm riêng cấp hai theo biến x và y. Ở ví dụ P1, có thể giải trực tiếp bằng cách lấy nguyên hàm. Tuy nhiên, phương pháp này chỉ thực hiện được trong không gian một chiều và không thể giải được trong trường hợp không gian có hơn hai chiều hoặc trong bài toán u + u’’ = f. Chính vì lí do này mà chúng ta sẽ phát triển phát triển PPPTHH cho trường hợp P1 và phác họa tổng quát của PPPTHH cho trường hợp P2. Lời giải sẽ bao gồm hai bước, nó phản ánh hai bước chủ yếu phải thực hiện để giải một bài toán biên bằng PPPTHH. Ở bước đầu tiên, chúng ta sẽ biểu diễn lại bài toán biên trong dạng gần đúng của nó hoặc dạng biến phân. Rất it hoặc không có máy tính được dùng để thực hiện bước này, việc này được làm bằng tay ở trên giấy. Bước thứ hai là rời rạc hóa, dạng gần đúng được rời rạc trong một không gian hữu hạn chiều. Sau bước thứ hai này, chúng ta sẽ có biểu thức cụ thể cho toàn bộ bài toán nhưng lời giải của bài toán trong không gian hữu hạn chiều tuyến tính chỉ là lời giải gần đúng của bài toán biên. Bài toán trong không gian hữu hạn chiều này sau đó được giải bằng máy tính. So sánh PPPTHH với phương pháp sai phân hữu hạn (PPSPHH) PPSPHH là một phương pháp khác để giải phương trình vi phân từng phần. Sự khác nhau giữa PPPTHH và PPSPHH là: • PPSPHH xấp xỉ bài toán phương trình vi phân; còn PPPTHH thì xấp xỉ lời giải của bài toán này • Điểm đặc trưng nhất của PPPTHH là nó có khả năng áp dụng cho những bài toán hình học và những bài toán biên phức tạp với mối quan hệ rời rạc. Trong khi đó PPSPHH về căn bản chỉ áp dụng được trong dạng hình chữ nhật với mối quan hệ đơn giản, việc vận dụng kiến thức hình học trong PPPTHH là đơn giản về lý thuyết. • Điểm đặc trưng của phương pháp sai phân hữu hạn là có thể dễ dàng thực hiện được. • Trong một vài trường hợp, PPSPHH có thể xem như là một tập con của PPPTHH xấp xỉ. Việc lựa chọn hàm cơ sở là hàm không đổi từng phần hoặc là hàm delta Dirac. Trong cả hai phương pháp xấp xỉ, việc xấp xỉ được tiến hành trên toàn miền, nhưng miền đó không cần liên tục. Như một sự lựa chọn, nó có thể xác định một hàm trên một miền rời rạc, với kết quả là toán tử vi phân liên tục không sinh ra chiều dài hơn, tuy nhiên việc xấp xỉ này không phải là PPPTHH. • Có những lập luận để lưu ý đến cơ sở toán học của việc xấp xỉ phần tử hữu hạn trở lên đúng đắn hơn, ví dụ, bởi vì trong PPSPHH đặc điểm của việc xấp xỉ những điểm lưới còn hạn chế. Phương pháp phần tử hữu hạn 3 • Kết quả của việc xấp xỉ bằng PPPTHH thường chính xác hơn PPSPHH, nhưng điều này còn phụ thuộc vào nhiều vấn đề khác và một số trường hợp đã cho kết quả trái ngược. Nói chung, PPPTHH là một phương pháp thích hợp để phân tích các bài toán về kết cấu (giải các bài toán về biến dạng và ứng suất của vật thể dạng khối hoặc động lực học kết cấu), trong khi đó phương pháp tính trong động lực học chất lỏng có khuynh hướng sử dụng PPSPHH hoặc những phương pháp khác (như phương pháp khối lượng hữu hạn).Những bài toán của động lực học chất lỏng thường yêu cầu phải rời rạc hóa bài toán thành một số lượng lớn những “ô vuông” hoặc những điểm lưới (hàng triệu hoặc hơn), vì vậy mà nó đòi hỏi cách giải phải đơn giản hơn để xấp xỉ các “ô vuông”. Điều này đặc biệt đúng cho các bài toán về dòng chảy ngoài, giống như dòng không khí bao quanh xe hơi hoặc máy bay, hoặc việc mô phỏng thời tiết ở một vùng rộng lớn. Có rất nhiều bộ phần mềm về phương pháp phần tử hữu hạn, một số miễn phí và một số được bán. Tài liệu tham khảo • Wikipedia tiếng Anh [1] • Chu Quốc Thắng, Phương pháp phần tử hữu hạn, NXB khoa học và kỹ thuật, 1997. Các liên kết ngoài Các phần mềm mã nguồn mở cho Phần tử hữu hạn bao gồm CAST3M, Z88, SLFFEA, YADE, FEniCS, deal.II, getFEM, libMesh, freeFEM, Elmer and Code-Aster. • CAST3M [2] • Diffpack [3] • Z88 [4] • SLFFEA [5] • YADE [6] • FEniCS [7] • deal.II [8] • getFEM [9] • libMesh [10] • freeFEM [11] • Code-Aster [12] • Impact [13] • IMTEK Mathematica Supplement (IMS) [14] • Calculix [15] • Elmer [16] • OOFEM [17]—a free, object oriented, general purpose FEM code • IFER [18]—Internet Finite Element Resources - an annotated list of FEA links and programs • Workshop "The Finite Element Method in Biomedical Engineering, Biomechanics and Related Fields" [19] Các phần mềm thương mại cho Phương pháp phần tử hữu hạn bao gồm ABAQUS, ANSYS, LS-DYNA, Nastran, Marc, and COMSOL Multiphysics, SAP2000, MIDAS, STAAP PRO, ETABS • ABAQUS, Inc. [20] • Analysis for Windows [21] • ANSYS, Inc. [22] • COMSOL, Inc. [23] - Provide COMSOL Multiphysics (formerly FEMLAB) • COSMOS/M and COSMOSWorks [24] - Integrated Design Validation in SolidWorks • CUBIT [25] - Hex/Tet meshing tool (Sandia National Laboratories) • ESI Pam-Crash [26] - Advanced crash simultion software Phương pháp phần tử hữu hạn 4 • FEMPRO [27] - Provide Algor • GLView Inova [28] • Infolytica Corp. [29] • Livermore Software Technology Corporation, LSTC [30] - Provide LS-DYNA and LS-PrePost [31] • MEGA [32] - Electromagnetics • Moldflow [33] - Provide Moldflow Plastics Insight [34] flow analysis software • MSC Software [35] - Provide Nastran and Marc • NAFEMS [36] - The International Association for the Engineering Analysis Community • NEiNastran [37] • NISA [38] Cranes Software International Ltd. • Plaxis [39] - Finite element code for soil and rock analysis • GEO5 [40] - Modelling of various geotechnical problems by the finite element method • Pro/MECHANICA [41] • Quickfield [42] - Free Demo (limited to 255 mesh nodes) • SOFiSTiK [43] - Structural FEA software based on AutoCAD • Strand7 [44] (aka Straus7 in Europe outside the UK) • StressCheck [45] - P-version FEA software • UGS NX Nastran [46] • VisualFEA [47] • Visual-Mesh [26] - Advanced meshing tool for crash, thermal and magnetic applications • SAMCEF [48] Chú thích [1] http:/ / en. wikipedia. org/ wiki/ Finite_element_method [2] http:/ / www-cast3m. cea. fr/ [3] http:/ / www. diffpack. com/ [4] http:/ / z88. org/ [5] http:/ / slffea. sourceforge. net/ [6] http:/ / yade. berlios. de/ [7] http:/ / www. fenics. org [8] http:/ / www. dealii. org [9] http:/ / www-gmm. insa-toulouse. fr/ getfem/ [10] http:/ / libmesh. sourceforge. net [11] http:/ / www. freefem. org [12] http:/ / www. code-aster. org [13] http:/ / impact. sourceforge. net/ [14] http:/ / www. imtek. uni-freiburg. de/ simulation/ mathematica/ IMSweb/ [15] http:/ / www. calculix. de [16] http:/ / www. csc. fi/ elmer [17] http:/ / www. oofem. org [18] http:/ / homepage. usask. ca/ ~ijm451/ finite/ fe_resources/ fe_resources. html [19] http:/ / www. uni-ulm. de/ uni/ intgruppen/ fem/ [20] http:/ / www. abaqus. com [21] http:/ / www. cuylaerts. net [22] http:/ / www. ansys. com [23] http:/ / www. comsol. com [24] http:/ / www. cosmosm. com [25] http:/ / cubit. sandia. gov [26] http:/ / www. esi-group. com/ [27] http:/ / www. algor. com [28] http:/ / www. fe-design. de/ en/ glview/ glview. html [29] http:/ / www. infolytica. com [30] http:/ / www. lstc. com Phương pháp phần tử hữu hạn 5 [31] http:/ / www. lstc. com/ lspp [32] http:/ / www. bath. ac. uk/ elec-eng/ aerc/ [33] http:/ / www. moldflow. com [34] http:/ / www. moldflow. com/ stp/ english/ products/ mpi. htm [35] http:/ / www. mscsoftware. com/ products/ products. cfm?Q=132& Z=149 [36] http:/ / www. nafems. org [37] http:/ / www. nenastran. com [38] http:/ / www. nisasoftware. com [39] http:/ / www. plaxis. nl [40] http:/ / www. finesoftware. eu [41] http:/ / www. ptc. com/ appserver/ mkt/ products/ home. jsp?k=1183 [42] http:/ / www. quickfield. com [43] http:/ / www. sofistik. com. [44] http:/ / www. strand7. com [45] http:/ / www. esrd. com [46] http:/ / www. ugs. com/ products/ nx/ simulation/ advanced/ nastran/ [47] http:/ / www. visualfea. com [48] http:/ / www. samcef. com Nguồn và người đóng góp vào bài 6 Nguồn và người đóng góp vào bài Phương pháp phần tử hữu hạn  Nguồn:  Người đóng góp: Micavn, Nguyenthephuc, Nguyễn Thanh Quang, Nguyễn Văn Minh, Pathétique, VnPenguin, 10 sửa đổi vô danh Giấy phép Creative Commons Attribution-Share Alike 3.0 //creativecommons.org/licenses/by-sa/3.0/