本發(fā)明專利技術(shù)提出了一種基于無(wú)網(wǎng)格模型計(jì)算周期結(jié)構(gòu)板聲學(xué)散射系數(shù)的方法，首先推導(dǎo)適合于周期結(jié)構(gòu)板節(jié)點(diǎn)聲壓計(jì)算的系統(tǒng)方程；另一方面，將周期結(jié)構(gòu)板及相同尺寸的參考板用節(jié)點(diǎn)進(jìn)行建模，然后利用移動(dòng)最小二乘法構(gòu)建形函數(shù)，結(jié)合兩者獲得節(jié)點(diǎn)處聲壓差后，即可求得接收點(diǎn)聲壓，進(jìn)而求得方向散射系數(shù)和平均散射系數(shù)。本發(fā)明專利技術(shù)將無(wú)網(wǎng)格法引入到周期結(jié)構(gòu)的聲學(xué)散射系數(shù)數(shù)值計(jì)算之中，避免了傳統(tǒng)數(shù)值方法中因?yàn)榫W(wǎng)格的存在而導(dǎo)致的一系列問(wèn)題。本發(fā)明專利技術(shù)有良好的自適應(yīng)性，在需要提高計(jì)算頻率上限時(shí)，可局部地增加節(jié)點(diǎn)密度，而不需對(duì)模型進(jìn)行重新劃分。通過(guò)與測(cè)量實(shí)驗(yàn)對(duì)比，驗(yàn)證了本發(fā)明專利技術(shù)具有較高的精度，因而在周期結(jié)構(gòu)散射系數(shù)數(shù)值計(jì)算中具有廣闊應(yīng)用前景。
【技術(shù)實(shí)現(xiàn)步驟摘要】
本專利技術(shù)涉及聲散射
，具體為，應(yīng)用于計(jì)算任意形狀周期結(jié)構(gòu)板的聲學(xué)散射系數(shù)。
技術(shù)介紹
周期結(jié)構(gòu)是一種特殊的界面聲散射結(jié)構(gòu)，它是指材料的幾何形狀在宏觀空間上按照周期規(guī)則排列的結(jié)構(gòu)。周期結(jié)構(gòu)在形式上早已存在，且因其良好的聲擴(kuò)散性能而應(yīng)用廣泛，多用于音樂(lè)廳、錄音室、劇院等對(duì)聲場(chǎng)擴(kuò)散要求較高的場(chǎng)所。周期結(jié)構(gòu)應(yīng)用的關(guān)鍵是其 散射系數(shù)的獲取。聲學(xué)散射系數(shù)的獲得通常有兩種方式，一種是測(cè)量實(shí)驗(yàn)方法，一種是數(shù)值計(jì)算方法。對(duì)于測(cè)量實(shí)驗(yàn)，直到2004年才建立國(guó)際標(biāo)準(zhǔn)，即混響室轉(zhuǎn)臺(tái)法，該方法還有諸多需要完善的方面，而且面對(duì)種類繁多的周期結(jié)構(gòu)，實(shí)驗(yàn)測(cè)量需要消耗大量的人力物力。另夕卜，對(duì)于處于設(shè)計(jì)階段的周期結(jié)構(gòu)是無(wú)法進(jìn)行實(shí)驗(yàn)測(cè)量的，因此研究散射系數(shù)的數(shù)值計(jì)算方法是非常必要的。但是由于理論上的限制，鏡面反射的聲能很難通過(guò)計(jì)算得到，因此數(shù)值計(jì)算的發(fā)展一直較為緩慢。2000年E. Mommertz提出計(jì)算鏡面反射聲能時(shí)可以用一個(gè)與周期結(jié)構(gòu)尺寸相同的純平的結(jié)構(gòu)代替周期結(jié)構(gòu)，由此發(fā)展了散射系數(shù)的數(shù)值算法。此后，對(duì)散射系數(shù)的數(shù)值計(jì)算獲得了快速發(fā)展以及越來(lái)越多的關(guān)注。目前關(guān)于散射系數(shù)的數(shù)值計(jì)算主要是采用邊界元法等基于網(wǎng)格的方法，由于網(wǎng)格的存在，這些方法存在一些固有缺陷首先是前處理困難，計(jì)算之前要生成網(wǎng)格模型，使得數(shù)據(jù)準(zhǔn)備的工作量大，尤其對(duì)于比較復(fù)雜的結(jié)構(gòu)模型，容易出現(xiàn)畸變網(wǎng)格；第二，這類方法普遍采用低階多項(xiàng)式作為聲壓函數(shù)插值函數(shù)，不可能對(duì)較高頻率聲波傳播問(wèn)題給出很好的近似；第三，計(jì)算結(jié)果不是光滑連續(xù)的，需要進(jìn)行光順化后處理。
技術(shù)實(shí)現(xiàn)思路
要解決的技術(shù)問(wèn)題為解決現(xiàn)有技術(shù)存在的問(wèn)題，本專利技術(shù)提出了一種。無(wú)網(wǎng)格法是近年來(lái)在力學(xué)等領(lǐng)域發(fā)展起來(lái)的一種新型數(shù)值計(jì)算方法，它采用一組相互獨(dú)立的節(jié)點(diǎn)來(lái)離散求解區(qū)域，直接借助于離散點(diǎn)來(lái)構(gòu)造形函數(shù)，從而可以徹底或部分地消除網(wǎng)格。相對(duì)于有限元法、邊界元法，無(wú)網(wǎng)格法在以下方面具有優(yōu)勢(shì)一是前處理簡(jiǎn)便，因?yàn)閮H需利用節(jié)點(diǎn)描述模型，無(wú)需考慮節(jié)點(diǎn)之間的拓?fù)潢P(guān)系；二是采用緊支函數(shù)的無(wú)網(wǎng)格法可得到帶狀稀疏矩陣，適用于求解大型科學(xué)工程問(wèn)題；三是自適應(yīng)性好，可在高誤差區(qū)域靈活增加節(jié)點(diǎn)數(shù)目或提高插值函數(shù)的階次；四是可以提供連續(xù)性好、形式靈活的形函數(shù)，從而使計(jì)算結(jié)果光滑連續(xù)，無(wú)需光順化后處理。本方法以無(wú)網(wǎng)格法為基礎(chǔ)，結(jié)合邊界元計(jì)算散射系數(shù)的思想，首先推導(dǎo)適合于周期結(jié)構(gòu)板節(jié)點(diǎn)聲壓計(jì)算的系統(tǒng)方程；另一方面，將周期結(jié)構(gòu)板及相同尺寸的參考板用節(jié)點(diǎn)進(jìn)行建模，然后利用移動(dòng)最小二乘法構(gòu)建形函數(shù)，結(jié)合兩者獲得節(jié)點(diǎn)處聲壓差后，即可求得接收點(diǎn)聲壓，進(jìn)而求得方向散射系數(shù)和平均散射系數(shù)。技術(shù)方案本專利技術(shù)的技術(shù)方案為所述一種，其特征在于包括以下步驟步驟I :將待計(jì)算的周期結(jié)構(gòu)板劃分為包含η個(gè)均勻分布節(jié)點(diǎn)的周期結(jié)構(gòu)板節(jié)點(diǎn)模型；建立一個(gè)純平的參考板，參考板與周期結(jié)構(gòu)板的平面投影形狀相同、厚度相同，將所述參考板也劃分為包含η個(gè)均勻分布節(jié)點(diǎn)的參考板節(jié)點(diǎn)模型；步驟2 :在以周期結(jié)構(gòu)板節(jié)點(diǎn)模型幾何中心為球心的一個(gè)半球面上均勻設(shè)置各不少于100個(gè)聲源點(diǎn)及接收點(diǎn)，所述半球面變徑不小于周期結(jié)構(gòu)板幾何中心到周期結(jié)構(gòu)板邊緣最大距離的兩倍；也在以參考板節(jié)點(diǎn)模型幾何中心為球心的一個(gè)半球面上設(shè)置相同位置和數(shù)量的聲源點(diǎn)及接收點(diǎn)；步驟3 :利用移動(dòng)最小二乘法，建立周期結(jié)構(gòu)板節(jié)點(diǎn)模型和參考板節(jié)點(diǎn)模型的形函數(shù)；步驟4:分別計(jì)算周期結(jié)構(gòu)板節(jié)點(diǎn)模型和參考板節(jié)點(diǎn)模型對(duì)應(yīng)的系統(tǒng)方程，得到周期結(jié)構(gòu)板節(jié)點(diǎn)模型及參考板節(jié)點(diǎn)模型上節(jié)點(diǎn)處的聲壓差；所述系統(tǒng)方程通過(guò)離散Helmholtz微分方程得到，系統(tǒng)方程的形式為權(quán)利要求1.一種，其特征在于包括以下步驟 步驟I :將待計(jì)算的周期結(jié)構(gòu)板劃分為包含η個(gè)均勻分布節(jié)點(diǎn)的周期結(jié)構(gòu)板節(jié)點(diǎn)模型；建立一個(gè)純平的參考板，參考板與周期結(jié)構(gòu)板的平面投影形狀相同、厚度相同，將所述參考板也劃分為包含η個(gè)均勻分布節(jié)點(diǎn)的參考板節(jié)點(diǎn)模型； 步驟2 :在以周期結(jié)構(gòu)板節(jié)點(diǎn)模型幾何中心為球心的一個(gè)半球面上均勻設(shè)置各不少于100個(gè)聲源點(diǎn)及接收點(diǎn)，所述半球面變徑不小于周期結(jié)構(gòu)板幾何中心到周期結(jié)構(gòu)板邊緣最大距離的兩倍；也在以參考板節(jié)點(diǎn)模型幾何中心為球心的一個(gè)半球面上設(shè)置相同位置和數(shù)量的聲源點(diǎn)及接收點(diǎn)； 步驟3 :利用移動(dòng)最小二乘法，建立周期結(jié)構(gòu)板節(jié)點(diǎn)模型和參考板節(jié)點(diǎn)模型的形函數(shù)；步驟4 :分別計(jì)算周期結(jié)構(gòu)板節(jié)點(diǎn)模型和參考板節(jié)點(diǎn)模型對(duì)應(yīng)的系統(tǒng)方程，得到周期結(jié)構(gòu)板節(jié)點(diǎn)模型及參考板節(jié)點(diǎn)模型上節(jié)點(diǎn)處的聲壓差； 所述系統(tǒng)方程通過(guò)離散Helmholtz微分方程得到，系統(tǒng)方程的形式為(C+D) · pd = F 式中C、D為nXn階的系數(shù)矩陣，Pd為周期結(jié)構(gòu)板節(jié)點(diǎn)模型或參考板節(jié)點(diǎn)模型上所有節(jié)點(diǎn)的聲壓差、F為ηΧ I階向量的載荷矩陣；其中全文摘要本專利技術(shù)提出了一種，首先推導(dǎo)適合于周期結(jié)構(gòu)板節(jié)點(diǎn)聲壓計(jì)算的系統(tǒng)方程；另一方面，將周期結(jié)構(gòu)板及相同尺寸的參考板用節(jié)點(diǎn)進(jìn)行建模，然后利用移動(dòng)最小二乘法構(gòu)建形函數(shù)，結(jié)合兩者獲得節(jié)點(diǎn)處聲壓差后，即可求得接收點(diǎn)聲壓，進(jìn)而求得方向散射系數(shù)和平均散射系數(shù)。本專利技術(shù)將無(wú)網(wǎng)格法引入到周期結(jié)構(gòu)的聲學(xué)散射系數(shù)數(shù)值計(jì)算之中，避免了傳統(tǒng)數(shù)值方法中因?yàn)榫W(wǎng)格的存在而導(dǎo)致的一系列問(wèn)題。本專利技術(shù)有良好的自適應(yīng)性，在需要提高計(jì)算頻率上限時(shí)，可局部地增加節(jié)點(diǎn)密度，而不需對(duì)模型進(jìn)行重新劃分。通過(guò)與測(cè)量實(shí)驗(yàn)對(duì)比，驗(yàn)證了本專利技術(shù)具有較高的精度，因而在周期結(jié)構(gòu)散射系數(shù)數(shù)值計(jì)算中具有廣闊應(yīng)用前景。文檔編號(hào)G06F19/00GK102938017SQ201210388230公開日2013年2月20日 申請(qǐng)日期2012年10月12日 優(yōu)先權(quán)日2012年10月12日專利技術(shù)者曾向陽(yáng), 王海濤 申請(qǐng)人:西北工業(yè)大學(xué)本文檔來(lái)自技高網(wǎng)...
【技術(shù)保護(hù)點(diǎn)】
一種基于無(wú)網(wǎng)格模型計(jì)算周期結(jié)構(gòu)板聲學(xué)散射系數(shù)的方法，其特征在于：包括以下步驟：步驟1：將待計(jì)算的周期結(jié)構(gòu)板劃分為包含n個(gè)均勻分布節(jié)點(diǎn)的周期結(jié)構(gòu)板節(jié)點(diǎn)模型；建立一個(gè)純平的參考板，參考板與周期結(jié)構(gòu)板的平面投影形狀相同、厚度相同，將所述參考板也劃分為包含n個(gè)均勻分布節(jié)點(diǎn)的參考板節(jié)點(diǎn)模型；步驟2：在以周期結(jié)構(gòu)板節(jié)點(diǎn)模型幾何中心為球心的一個(gè)半球面上均勻設(shè)置各不少于100個(gè)聲源點(diǎn)及接收點(diǎn)，所述半球面變徑不小于周期結(jié)構(gòu)板幾何中心到周期結(jié)構(gòu)板邊緣最大距離的兩倍；也在以參考板節(jié)點(diǎn)模型幾何中心為球心的一個(gè)半球面上設(shè)置相同位置和數(shù)量的聲源點(diǎn)及接收點(diǎn)；步驟3：利用移動(dòng)最小二乘法，建立周期結(jié)構(gòu)板節(jié)點(diǎn)模型和參考板節(jié)點(diǎn)模型的形函數(shù)；步驟4：分別計(jì)算周期結(jié)構(gòu)板節(jié)點(diǎn)模型和參考板節(jié)點(diǎn)模型對(duì)應(yīng)的系統(tǒng)方程，得到周期結(jié)構(gòu)板節(jié)點(diǎn)模型及參考板節(jié)點(diǎn)模型上節(jié)點(diǎn)處的聲壓差；所述系統(tǒng)方程通過(guò)離散Helmholtz微分方程得到，系統(tǒng)方程的形式為：(C+D)·pd＝F式中C、D為n×n階的系數(shù)矩陣，pd為周期結(jié)構(gòu)板節(jié)點(diǎn)模型或參考板節(jié)點(diǎn)模型上所有節(jié)點(diǎn)的聲壓差、F為n×1階向量的載荷矩陣；其中c(Pi)=1P∈Ω,P∉Γα/4πP∈Γ0P∉Ω,α為節(jié)點(diǎn)Pi處表面的立體角；Dij=∫ΓNQj(∂G(Pi,Q)∂nQ+β∂G(Pi,Q)∂nPi∂nQ)dΓ,Fi=jρωqω(r0)[G(Pi,r0)+jk∂G(Pi,r0)∂nPi],為格林函數(shù)，rPQ表示P，Q之間的距離，N為周期結(jié)構(gòu)板節(jié)點(diǎn)模型或參考板節(jié)點(diǎn)模型的形函數(shù)，ρ為空氣密度，ω為圓頻率，k為波數(shù)；步驟5：利用下式和步驟4得到的周期結(jié)構(gòu)板節(jié)點(diǎn)模型及參考板節(jié)點(diǎn)模型上節(jié)點(diǎn)處的聲壓差分別計(jì)算周期結(jié)構(gòu)板節(jié)點(diǎn)模型和參考板節(jié)點(diǎn)模型對(duì)應(yīng)的接收點(diǎn)的聲壓：pω(R)＝？B·pd其中B=∫ΓNQ∂G(R,Q)∂nQdΓ;步驟6：根據(jù)下式計(jì)算每個(gè)聲源點(diǎn)的方向散射系數(shù)：式中：θ和分別代表聲源點(diǎn)相對(duì)于周期結(jié)構(gòu)板節(jié)點(diǎn)模型幾何中心的俯仰角與方位角；n為接收點(diǎn)的數(shù)量；θ“和分別代表第i個(gè)接收點(diǎn)的俯仰角與方位角；p1是周期結(jié)構(gòu)板節(jié)點(diǎn)模型所對(duì)應(yīng)的接收點(diǎn)的聲壓；p0是參考板節(jié)點(diǎn)模型所對(duì)應(yīng)的接收點(diǎn)的聲壓；*代表復(fù)共軛；步驟7：得到所有聲源點(diǎn)的方向散射系數(shù)之后，按照下式計(jì)算平均散射系數(shù)：FDA00002247675200011.jpg,FDA00002247675200021.jpg,FDA00002247675200023.jpg,FDA00002247675200024.jpg,FDA00002247675200025.jpg,FDA00002247675200026.jpg...
【技術(shù)特征摘要】
【專利技術(shù)屬性】
技術(shù)研發(fā)人員：曾向陽(yáng)，王海濤，
申請(qǐng)(專利權(quán))人：西北工業(yè)大學(xué)，
類型：發(fā)明
國(guó)別省市：
還沒(méi)有人留言評(píng)論。發(fā)表了對(duì)其他瀏覽者有用的留言會(huì)獲得科技券。