Konstruksi dasar elemen hingga hingga / sesuai untuk mesh segitiga atau tetrahedral

9

Dalam makalah Hierarchical Conforming Finite Element Methods untuk Persamaan Biharmonic , P. Oswald mengklaim elemen tipe Clough-Tocher memiliki kontinuitas sambil menjadi polinomial kubik pada setiap segitiga. Dia tidak memberikan satu set fungsi dasar eksplisit hanya derajat standar kebebasan pada titik-titik quadrature.C1

Demikian pula, dalam buku Teori Matematika Metode Elemen Hingga Bab 3, penulis memberi kita konstruksi elemen hingga Hermit kubik, tetapi mereka tidak menyebutkan kontinuitas elemen Hermit kubik.

Namun, dalam makalah Kompleks diferensial dan stabilitas numerik , Doulgas Arnold mengusulkan bahwa untuk ruang diskrit / sesuai, kita harus menggunakan elemen hingga Hermite quintic (atau lebih tepatnya Argyris), yang sangat rumit untuk diekspresikan secara eksplisit.H 2C1H2

Jadi inilah pertanyaanku:

(1) Apakah ada kertas yang muncul dengan formula eksplisit untuk elemen hingga / sesuai pada mesh segitiga atau tetrahedral?H 2C1H2

(2) Apakah sebaiknya kubik menjadi tingkat minimum persyaratan polinomial untuk kontinuitas ?C1

Shuhao Cao
sumber

Jawaban:

5

Elemen hermit kubik memiliki turunan normal kontinu tetapi tidak kontinuitas penuh . Secara khusus, turunan normal mungkin tidak cocok pada batas dua elemen, jauh dari simpul. Jika Anda ingin kontinuitas penuh Anda harus menggunakan elemen Argyris atau Hsieh-Clough-Tucker atau sesuatu. Saya merekomendasikan diskusi di bab 6 buku elemen hingga Ciarlet.C 1C1C1

Tingkat polinomial yang diperlukan untuk kontinuitas akan tergantung pada dimensi spasial Anda, tetapi dalam 2D ​​atau 3D saya tidak berpikir Anda bisa lolos dengan polinomial kurang dari kubik. Anda dapat mempertimbangkan beberapa jenis metode yang tidak sesuai yang memungkinkan ruang elemen hingga yang lebih sederhana.C1

Andrew T. Barker
sumber
Err, jika suatu fungsi kontinu melintasi antarmuka antara dua sel, dan jika fungsi pada setiap sel dalam sebagaimana jika itu adalah polinomial, lalu bagaimana turunan tangensial dapat terputus pada antarmuka sel? Atau maksud Anda bahwa turunan tangensial dapat diskontinyu pada simpul, yaitu titik akhir dari setiap antarmuka ? C
Wolfgang Bangerth
Anda benar sekali, saya mengedit jawabannya.
Andrew T. Barker
3

Saya merujuk Anda ke buku Splines on Triangulations . Saya tidak dapat menemukan salinan saya saat ini untuk memberikan jawaban yang lebih baik, tetapi saya ingat diskusi / teorema pada urutan polinomial yang diperlukan untuk ruang . Jika saya ingat dengan benar, Lai membuktikan bahwa dalam kondisi tertentu tidak apa-apa, tetapi selalu cukup. p = 3 p = 5C1p=3p=5

Sayangnya, saya juga ingat bahwa Lai kemudian tidak menunjukkan bagaimana membangun ruang , hanya membuktikan mereka ada diberikan triangulasi dan ruang spline. Setelah ia memiliki bukti ini, ia menyelesaikan aplikasinya dengan persamaan kendala linear tambahan untuk menegakkan kondisi .C 1C1C1

Nathan
sumber
selamat datang di scicomp Mr. Collier :)
Aron Ahmadia