Di kelas FEM, biasanya diterima begitu saja bahwa matriks kekakuan adalah pasti positif, tapi saya tidak bisa mengerti mengapa. Adakah yang bisa memberikan penjelasan?
Sebagai contoh, kita bisa mempertimbangkan masalah Poisson: yang matriks kekakuannya adalah: yang simetris dan pasti positif. Simetri adalah sifat yang jelas, tetapi kepastian positif tidak begitu jelas bagi saya.
finite-element
matrix
stiffness
pengguna123
sumber
sumber
Jawaban:
Properti berikut dari properti persamaan diferensial parsial (bentuk lemah); ini adalah salah satu kelebihan dari metode elemen hingga dibandingkan dengan, misalnya, metode beda hingga.
Untuk melihat itu, pertama-tama ingat bahwa metode elemen hingga dimulai dari bentuk lemah dari persamaan Poisson (Saya mengasumsikan kondisi batas Dirichlet di sini): Cari sedemikian sehingga Properti penting di sini adalah (Ini mengikuti ketidaksetaraan Poincaré.)a ( u , v ) : = ∫ Ω ∇ u ⋅ ∇ vu∈H10(Ω)
Sekarang pendekatan elemen hingga klasik adalah untuk mengganti ruang dimensi tak terbatas dengan subruang dimensi-terbatas dan temukan sedemikian sehingga Properti penting di sini adalah bahwa Anda menggunakan dan subruang yang sama ( diskritisasi yang sesuai ); itu artinya Anda masih memilikiH10( Ω ) Vh⊂ H10( Ω ) kamuh∈ Vh
Sekarang untuk langkah terakhir: Untuk mengubah bentuk variasi menjadi sistem persamaan linear, Anda memilih basis dari , tulis dan masukkan , ke dalam . Matriks kekakuan kemudian memiliki entri (yang bertepatan dengan apa yang Anda tulis).{φ1,…,φN} Vh uh=∑Ni=1uiφi vh=φj 1≤j≤N (2) K Kij=a(φi,φj)
Sekarang ambil vektor sembarang dan set . Kemudian kita miliki dengan dan bilinearity dari (yaitu, Anda dapat memindahkan skalar dan jumlah ke kedua argumen) Karena arbitrer, ini menyiratkan bahwa adalah positif pasti.v⃗ =(v1,…,vN)T∈RN vh:=∑Ni=1viφi∈Vh (3) a → v K
TL; DR: Matriks kekakuan adalah pasti positif karena berasal dari diskritisasi yang sesuai dari persamaan diferensial parsial elliptic (self-adjoint) .
sumber
Jika kekakuan elemen tidak positif, maka sistem tidak stabil. Jadi modelnya kemungkinan besar tidak benar. Lihatlah persamaan osilator harmonik yang paling dasar
Solusinya tidak stabil jika negatif (lihat akar dari persamaan karakteristik). Itu artinya solusinya akan meledak. Kekakuan harus menjadi kekuatan pemulih. Setidaknya untuk pegas fisik. Matriks kekakuan meluas ke sejumlah besar elemen (matriks kekakuan global). Itu semuanya. Tapi itu ide dasar yang sama. Dasar FEM adalah dalam metode matriks kekakuan untuk analisis struktural di mana setiap elemen memiliki kekakuan yang terkait dengannya.k
sumber