Helmholtz方程的有限元方法.pdf

1、Helmholtz方程常出現(xiàn)在聲學(xué)、彈性力學(xué)、電磁學(xué)等領(lǐng)域中，在科學(xué)研究中有重要的實(shí)用價(jià)值。對(duì)于它的研究，大量論文資料給出了數(shù)值求解方法的理論分析，本論文將在部分文獻(xiàn)已推導(dǎo)出的理論基礎(chǔ)上，對(duì)Helmholtz方程的變分形式進(jìn)行分析總結(jié)，并將側(cè)重點(diǎn)放在對(duì)Helmholtz方程的有限元方法數(shù)值求解實(shí)驗(yàn)上，并通過(guò)數(shù)值實(shí)驗(yàn)展示各種方法的數(shù)值解精度。
　　對(duì)于一維問(wèn)題，本文采用了線性插值、Hermite插值和三角插值構(gòu)造單元上的形函數(shù)，并應(yīng)

2、用數(shù)學(xué)軟件MATLAB編制計(jì)算程序進(jìn)行了數(shù)值實(shí)驗(yàn)，針對(duì)波數(shù)k的不同取值，通過(guò)誤差數(shù)據(jù)展示了數(shù)值解的精度。
　　對(duì)于二維Helmholtz方程，文中分別采用了三角形單元剖分和矩形單元剖分。對(duì)于三角形單元，文中構(gòu)造了線性形函數(shù)；對(duì)于矩形單元，文中采用雙線性插值、Hermite插值和三角插值構(gòu)造了單元形函數(shù)。另外，還利用齊次Helmholtz方程的解作為基函數(shù)，導(dǎo)出了相應(yīng)的形函數(shù)。最后編制了程序進(jìn)行數(shù)值實(shí)驗(yàn)，針對(duì)波數(shù)k的不同取值，通過(guò)誤

3、差數(shù)據(jù)展示了各種情況下數(shù)值解的精度。
　　對(duì)于三維Helmholtz方程，本文先對(duì)求解區(qū)域進(jìn)行了三棱柱單元剖分，并構(gòu)造了單元上的線性形函數(shù)，接著又采用長(zhǎng)方體單元對(duì)求解區(qū)域進(jìn)行了剖分，分別采用線性插值、Hermite插值和三角插值構(gòu)造了單元形函數(shù)。并編制程序進(jìn)行了數(shù)值實(shí)驗(yàn)，針對(duì)波數(shù)k的不同取值，通過(guò)誤差數(shù)據(jù)展示了數(shù)值解的精度。
　　同時(shí)，本文繪制了一維問(wèn)題數(shù)值解的誤差曲線和二維問(wèn)題數(shù)值解的誤差曲面，清晰地反映了數(shù)值解與解析解之