Secara matematis, mengapa matriks massa / pemuatan penggumpalan vektor bekerja?

13

Saya tahu bahwa orang sering mengganti matriks massa yang konsisten dengan matriks diagonal yang disatukan. Di masa lalu, saya juga menerapkan kode di mana vektor beban dirakit dalam mode yang disatukan daripada mode yang konsisten dengan FEM. Tapi saya tidak pernah melihat mengapa kita diizinkan melakukan ini sejak awal.

Apa intuisi di balik lumping yang memungkinkan seseorang untuk menerapkannya pada massa dan memuat vektor? Apa pembenaran matematis untuk itu? Dalam situasi apa penggumpalan tidak diizinkan / bukan perkiraan yang baik untuk vektor massa dan beban?

Paul
sumber

Jawaban:

16

Dalam metode elemen hingga, entri matriks dan entri sisi kanan didefinisikan sebagai integral. Kita dapat, secara umum, tidak menghitung ini dengan tepat dan menerapkan quadrature. Tetapi ada banyak rumus quadrature yang bisa dipilih, dan sering kali dipilih dengan cara sedemikian rupa sehingga (i) kesalahan yang diperkenalkan oleh quadrature memiliki urutan yang sama dengan yang disebabkan oleh diskritisasi, atau setidaknya tidak jauh lebih buruk, dan (ii) matriks memiliki sifat-sifat tertentu yang ternyata nyaman.

Mass lumping adalah contoh kerja ini: Jika seseorang memilih rumus quadrature tertentu (yaitu, yang dengan titik quadrature terletak di titik interpolasi dari elemen hingga), maka matriks massa yang dihasilkan terjadi diagonal. Itu cukup nyaman untuk implementasi komputasi, dan alasan mengapa orang menggunakan rumus quadrature ini. Itu juga alasan mengapa itu "bekerja": Pilihan formula quadrature ini masih memiliki urutan yang cukup tinggi.

Wolfgang Bangerth
sumber
Jawaban yang luar biasa, seperti biasa. Saya juga akan sangat tertarik dengan pendapat Anda pada bagian kedua dari pertanyaan, ketika lumping tidak diperbolehkan / perkiraan buruk , jika ada sesuatu yang terlintas dalam pikiran.
Anton Menshov
2
@AntonMenshov: Tampaknya akan sulit (mungkin tidak mungkin?) Untuk mendapatkan perkiraan yang baik melalui penggumpalan untuk elemen tingkat tinggi, karena (misalnya diagonal) penggumpalan dalam kasus itu akan menjadi setara dengan kuadratur orde rendah yang diterapkan pada orde lebih tinggi polinomial.
Paul
@ WolfgangBangerth: Saya pikir saya mengerti sekarang. Jadi, ini seperti menggunakan aturan newton-cotes untuk integrasi alih-alih gaussian quadrature. Karena setiap fungsi interpolasi lagrange memiliki nilai satuan pada satu node tertentu, migrasi titik kuadratur ke node hanya menghasilkan istilah diagonal yang menjadi nol (setidaknya, untuk elemen linier).
Paul
1
@ Paul, saya telah menggunakan lumping massal untuk elemen tingkat tinggi (hingga 14) dan berfungsi dengan baik. Merupakan kebiasaan untuk Metode Elemen Spectral, di mana Anda menggunakan node sebagai titik quadrature. Dalam hal ini, urutan aproksimasi . 2n3
nicoguaro
1
Poin penting adalah bahwa untuk elemen orde tinggi, kita perlu mendefinisikan matriks massa "dikelompokkan" melalui rumus quadrature tertentu. Bentuk asli - dari mana istilah "lumped" berasal - menambahkan entri off-diagonal ke diagonal, tetapi itu hanya berfungsi jika semuanya positif. Jika Anda menerapkan Gaussian quadrature, ini berlaku untuk elemen tingkat terendah, tetapi tidak untuk elemen tingkat tinggi.
Wolfgang Bangerth
5

Matriks diagonal memiliki keuntungan yang jelas dalam mempercepat perhitungan numerik, dan jawaban Wolfgang Bangerth adalah penjelasan yang baik tentang bagaimana menghitung matriks massa diagonal, tetapi itu tidak menjawab pertanyaan OP "mengapa ini bekerja " dalam arti "mengapa itu perkiraan yang baik untuk fisika yang Anda modelkan ".

Secara konseptual, Anda dapat memisahkan respons elemen menjadi tiga bagian: gerak translasi benda tegar, rotasi kaku tentang pusat elemen massa, dan deformasi elemen.

12vTMvv

aa3a5

Oleh karena itu, Anda hanya benar-benar memerlukan pendekatan "baik" ke bagian tubuh yang kaku dari gerakan, yaitu 6 DOF, dan pada kenyataannya pendekatan yang baik untuk hanya KE dari terjemahan tubuh yang kaku , yaitu 3 DOF, akan menyatu ketika ukuran elemen adalah berkurang.

Istilah diagonal dari matriks elemen mengandung lebih dari cukup parameter independen untuk mewakili 3 atau 6 syarat KE dengan akurasi yang cukup. Bahkan untuk elemen tingkat tinggi, Anda dapat menggunakan matriks massa diagonal massa di mana istilah diagonal untuk node sisi tengah adalah nol.

Perhatikan bahwa ini adalah situasi yang sangat berbeda dari energi potensial elemen, di mana kontribusi dari translasi dan rotasi benda tegar adalah nol, dan satu-satunya hal yang penting adalah mewakili energi regangan yang sesuai dengan deformasi elemen . Matriks kekakuan diagonal karena itu tidak akan menjadi ide yang layak!

alephzero
sumber
5

Selain jawaban lain, ada skenario di mana kesalahan dalam matriks massa tidak mempengaruhi hasil yang diinginkan.

K(u) u=f(u)u^K(u) u+C(u) u˙+M u¨=f(u)MCu˙=u¨=0M

MM1

1 Walaupun pemikiran tentang perilaku fisik yang dinamis tentu saja lebih mudah dengan matriks massa "yang benar" - misalnya momentum sudut mungkin secara tidak tepat dilestarikan oleh matriks massa yang disatukan.

Max Langhof
sumber
1
Cu˙Cu˙Mu¨Cu˙C
u¨