Phương pháp phần tử hữu hạn so với Phương pháp phần tử hữu hạn mở rộng (FEM so với XFEM)


9

Sự khác biệt chính giữa FEM và XFEM là gì? Khi nào chúng ta (không) nên sử dụng XFEM intead của FEM và ngược lại? Nói cách khác, khi tôi gặp một vấn đề mới, làm thế nào tôi có thể biết sử dụng cái nào trong số chúng?


1
Hầu hết các lần tôi gặp XFEM, đó là để đối phó với sự không liên tục liên quan đến lan truyền vết nứt và gãy trong cơ học rắn. Tôi chưa thực sự thấy nó được sử dụng bên ngoài ứng dụng này.
Paul

Trên thực tế, vẫn còn nhiều lĩnh vực khác cũng sử dụng XFEM để giải quyết. Đó là lý do tại sao tôi cần biết cách nhận ra phương pháp này bất cứ khi nào tôi bắt đầu giải quyết vấn đề.
Anh-Thị DINH

Câu trả lời:


11

Phương pháp phần tử hữu hạn (FEM) là phương thức cha mẹ đã truyền cảm hứng cho nhiều, nhiều phương thức và phương thức khác thực sự là FEM nhưng giả vờ là không.

Trong phương pháp phần tử hữu hạn, "các hàm hình dạng" được sử dụng để cung cấp một không gian gần đúng để giải pháp có thể được biểu diễn bằng một vectơ. Trong FEM cổ điển, các hàm hình dạng này là đa thức.

Trong Phương pháp phần tử hữu hạn mở rộng (XFEM), các hàm "làm giàu" bổ sung được sử dụng để tính gần đúng giải pháp bên cạnh các hàm hình dạng đa thức. Các hàm làm giàu này được chọn để có các thuộc tính mà giải pháp được biết là tuân theo.

Các hàm làm giàu XFEM rõ ràng nhất là các hàm năng lượng được giới thiệu tại các vết nứt góc nhọn để biểu thị các điểm kỳ dị trong gradient giải pháp (nghĩa là điểm kỳ dị trong ứng suất cho các vấn đề cơ học vững chắc). XFEM có thể được sử dụng cho các chức năng làm giàu khác và các miền giải pháp khác (đáng chú ý là truyền nhiệt), nhưng tên đồng nghĩa với phân tích gãy xương.

Sự khác biệt giữa các phương pháp khác nhau - đây có phải là XFEM hay không?, V.v. - là khó khăn và tinh tế và không quan trọng.

Đối với việc sử dụng, XFEM thấy rất ít sử dụng thực tế. Có một số ứng dụng trong mã phần tử hữu hạn thực sự, đáng chú ý nhất là Abaqus, nhưng chúng chưa được chấp nhận rộng rãi.

Đối với hầu hết các vấn đề thực tế, FEM cổ điển sẽ được sử dụng. Đối với hầu hết các vấn đề phân tích gãy xương, FEM cổ điển vẫn có thể được sử dụng với sàng lọc lưới phù hợp và / hoặc sàng lọc p trong khu vực của đầu vết nứt. Các mô hình gãy xương khác, ít nghiêm ngặt hơn cũng có thể được sử dụng.


Không muốn bỏ qua câu trả lời (xuất sắc) này, các hàm số ít đại diện cho các thành phần giải pháp ở các góc reentrant thực tế thuộc loại trong đó là khoảng cách đến góc, cho chuyển vị (giải pháp) và cho ứng suất (đạo hàm của nó). r 0 < α < 1 - 1 < α < 0rαr0<α<1-1<α<0
Wolfgang Bangerth

@WolfgangBangerth, Cảm ơn! Tôi đã chỉnh sửa câu trả lời của mình để nói "các hàm năng lượng", đó là những gì tôi muốn đặt ở vị trí đầu tiên, mặc dù nó vẫn không chính xác. Tôi gần như đặt sqrt (r) (cho một vết nứt kín) để vẽ một bức tranh rõ ràng hơn, nhưng tôi không chắc liệu nó có gây mất tập trung không. Tôi biết có rất nhiều chi tiết, tôi biết, về việc triển khai XFEM nghiêm túc (một số tôi đã nghiên cứu và những người khác thì không).
Mike

@Mike: một câu hỏi ít liên quan hơn là sự khác biệt giữa FEM bong bóng P1 và XFEM là gì? Bạn có thể cho tôi xem?
Anh-Thị DINH

@PoBo, Có chút tương đồng. Cả hai phương pháp đều liên quan đến việc thêm các hàm hình dạng mà không thay đổi lưới và cả hai đều dựa trên cùng một phép toán cơ bản chung cho toàn bộ họ FEM, nhưng đó là nơi tương tự kết thúc.
Mike

Nếu bạn không hiểu rõ về phiên bản p hoặc phương pháp tiếp cận chức năng hình dạng bong bóng P1, bạn có thể thử một câu hỏi cấp cao khác hoặc chọn một trong những cuốn sách về nó (nói chung của Szabo và Babuska khá nghiêm ngặt, nhưng ít hơn nhiều so với những phiên bản khác bao gồm phiên bản p.)
Mike

7

Cả hai của Mike câu trả lời và Jed của một mô tả tốt sự phân đôi XFEM / FEM và chính xác chỉ ra rằng khu vực quan trọng nhất của ứng dụng là 3D Gãy Mechanics, nơi bạn có một vết nứt, tức là trọng lượng rẽ gián đoạn trên một bề mặt bên trong tên miền của bạn.

Các vết nứt khó mô hình hóa trong FEM cổ điển vì hai lý do:

  1. Lưới phải đồng nhất trên vết nứt: chính xác hơn là vết nứt phải nằm ở ranh giới của một tên miền phụ của FE. Các vết nứt không thể nằm bên trong (vượt qua) một yếu tố hữu hạn.

  2. Trường ứng suất đơn ở đầu vết nứt đòi hỏi các yếu tố đặc biệt và / hoặc kỹ thuật chia lưới (yếu tố điểm quý, lưới tập trung) để được mô hình hóa với độ chính xác tốt.

Từ quan điểm kỹ thuật trong cơ học gãy xương, bạn có hai loại vấn đề chính:

  1. Tính toán cường độ ứng suất ,

  2. phân tích lan truyền vết nứt, ví dụ trong phân tích dung sai hoặc mỏi.

Đối với các loại đầu tiên của vấn đề FEM cổ điển là nhiều hơn cần thiết và là những công cụ kỹ thuật tiêu chuẩn. (Điều này là do, may mắn thay, có các phương pháp năng lượng để đánh giá các yếu tố cường độ ứng suất không nhạy cảm với các lỗi số gần đầu vết nứt.)

Phân tích lan truyền crack là một câu chuyện hoàn toàn khác: trong hầu hết các trường hợp bạn không biết trước đường dẫn crack, và do đó việc làm lại thường xuyên là cần thiết. Lời hứa chính của XFEM là cho phép lan truyền vết nứt bên trong lưới FEM cố định , vết nứt cắt theo cách của anh ta không chỉ ở ranh giới giữa các tên miền phụ, mà bên trong chính FE.

XFEM là một kỹ thuật tương đối mới, vẫn còn lâu mới là một công cụ kỹ thuật tiêu chuẩn. Câu trả lời của tôi cho câu hỏi OP, ít nhất là trong phân tích cơ học và kỹ thuật rắn, là XFEM có một lĩnh vực ứng dụng rất hẹp và chuyên biệt trong phân tích lan truyền vết nứt và thiệt hại, đối với hình học 3D phức tạp, khi đường nứt có thể được ước tính là tiên nghiệm .

Tuy nhiên, để tôi nhấn mạnh rằng cơ học gãy xương là một lĩnh vực rất quan trọng trong kỹ thuật: ví dụ, máy bay ngày nay cũng an toàn vì có thể dự đoán số lượng thiệt hại và lan truyền vết nứt giữa các khoảng bảo trì. XFEM, hoặc các kỹ thuật mới tương tự, sẽ trở thành công cụ quan trọng trong tương lai gần.


tầm quan trọng của XFEM trong cơ học gãy xương được thể hiện bởi tất cả các bạn nhưng vẫn còn các lĩnh vực khác cần sử dụng XFEM thay vì FEM cổ điển? Ví dụ, trong sự tăng trưởng của màng sinh học, giao diện của màng sinh học trong chất nền thay đổi theo thời gian. Ranh giới có thể thay đổi (di chuyển ranh giới). Nếu chúng ta sử dụng FEM cổ điển, lưới phải được tạo ra ở mỗi bước thời gian, điều đó có đúng không? Điều đó thực sự không tốt, đặc biệt là trong trường hợp 3D. Hoặc nếu chúng ta xem xét 2 pha chất lỏng với độ dốc nồng độ khác nhau, có vẻ như cũng cần sử dụng XFEM?
Anh-Thị DINH

Có rất nhiều vấn đề trong đó bạn có bề mặt tự do hoặc ranh giới di chuyển, rất khó với các phương pháp Lagrangian thuần túy (do thường xuyên làm lại). XFEM nói thêm về mô hình hóa sự không liên tục trong các miền. Tôi biết các quy trình khớp nối sử dụng sự gián đoạn để thể hiện ranh giới di chuyển ... nhưng tôi không phải là chuyên gia trong các lĩnh vực này.
Stefano M

một câu hỏi ít liên quan hơn là sự khác biệt giữa FEM bong bóng P1 và XFEM là gì? Bạn có thể cho tôi xem?
Anh-Thị DINH

Tôi sẽ đề nghị mở một câu hỏi mới. Tóm lại P1-bubble / P1 là một yếu tố hữu hạn cụ thể (đối với giải pháp của phương trình Stokes) trong khi XFEM là một khái niệm tổng quát hơn, liên quan đến việc sử dụng các hàm làm giàu để mô hình hóa sự không liên tục, khai thác cách tiếp cận Phân vùng Thống nhất.
Stefano M

6

FEM là tập con của XFEM. XFEM là một phương pháp để làm phong phú các không gian phần tử hữu hạn để xử lý các vấn đề với sự không liên tục (chẳng hạn như gãy xương). Với FEM cổ điển, để đạt được độ chính xác tương tự thường đòi hỏi phải chia lưới phù hợp và tinh chỉnh thích ứng phức tạp, trong khi XFEM thực hiện với một lưới duy nhất, di chuyển độ phức tạp hình học đó vào các yếu tố (XFEM rất phức tạp để thực hiện, đặc biệt là trong 3D). Trong khi đó, XFEM dẫn đến các ma trận cực kỳ điều hòa đòi hỏi người giải quyết trực tiếp hoặc phương pháp đa trị rất chuyên biệt (ví dụ: Gerstenberger và Tuminaro (2012) ).


Liệu nỗ lực trong việc chuyển sự phức tạp từ lưới sang các chức năng hình dạng thực sự phải trả giá cuối cùng? Cả hai dường như phức tạp theo cùng một cách.
shuhalo

Như thường xảy ra trong khoa học tính toán, nó phụ thuộc vào người bạn hỏi và vấn đề bạn giải quyết. Nhiều học viên XFEM đã chơi khăm bằng cách sử dụng một phương trình bậc hai thô thay vì một phần thích nghi với sự không liên tục của phần tử.
Jed Brown
Khi sử dụng trang web của chúng tôi, bạn xác nhận rằng bạn đã đọc và hiểu Chính sách cookieChính sách bảo mật của chúng tôi.
Licensed under cc by-sa 3.0 with attribution required.