本發(fā)明屬于風(fēng)工程技術(shù)領(lǐng)域，具體涉及一種風(fēng)驅(qū)雨量的目標(biāo)軌跡追蹤計(jì)算方法。

背景技術(shù)：

風(fēng)驅(qū)雨是風(fēng)工程研究的一個(gè)熱點(diǎn)問題，風(fēng)驅(qū)雨的研究在建筑、交通等領(lǐng)域均具有重要意義。風(fēng)驅(qū)雨會(huì)影響建筑墻面的保溫性能，導(dǎo)致墻面發(fā)生霉變，甚至影響墻面使用壽命；風(fēng)驅(qū)雨還可能影響高速列車、飛機(jī)等交通工具的運(yùn)行安全。為評(píng)估風(fēng)驅(qū)雨的影響，就需要確定風(fēng)驅(qū)雨量的大小。

目前，采用數(shù)值模擬方法定量計(jì)算風(fēng)驅(qū)雨量大小是國內(nèi)外的研究熱點(diǎn)。風(fēng)驅(qū)雨數(shù)值模擬方法通常基于歐拉-拉格朗日和歐拉-歐拉兩類模型。由于在實(shí)際降雨條件下雨滴所占體積分?jǐn)?shù)遠(yuǎn)小于10％，不滿足通常認(rèn)為的采用歐拉-歐拉模型時(shí)各相體積分?jǐn)?shù)應(yīng)大于10％的要求，因此基于歐拉-拉格朗日模型的計(jì)算方法應(yīng)用更為廣泛，且通常認(rèn)為該模型能夠更準(zhǔn)確地模擬風(fēng)驅(qū)雨現(xiàn)象中雨滴的運(yùn)動(dòng)過程。在基于歐拉-拉格朗日模型的風(fēng)驅(qū)雨計(jì)算方法中，首先需要獲得建筑等研究對(duì)象的繞流流場(chǎng)。基于該流場(chǎng)，計(jì)算一定量離散雨滴的軌跡，其中一些軌跡應(yīng)終止于待研究的物體表面。最后根據(jù)終止于物面的雨滴軌跡計(jì)算該面上的風(fēng)驅(qū)雨量。表征風(fēng)驅(qū)雨量大小的參數(shù)主要有特定收集率和收集率。其中，特定收集率表示某一直徑雨滴在物面的風(fēng)驅(qū)雨強(qiáng)度與在未受擾動(dòng)水平面降雨強(qiáng)度之比；收集率表示將特定收集率按不同直徑雨滴體積分?jǐn)?shù)加權(quán)平均后的結(jié)果。

當(dāng)前對(duì)建筑表面風(fēng)驅(qū)雨量的研究主要考慮建筑迎風(fēng)面的收集率分布，僅有少量研究對(duì)建筑側(cè)面收集率分布進(jìn)行了分析，對(duì)建筑背風(fēng)面收集率分布則尚未見到研究報(bào)道。分析建筑繞流流場(chǎng)的速度分布特點(diǎn)可知，撞擊到建筑側(cè)面和背風(fēng)面的雨滴數(shù)量所占比例小于撞擊到建筑迎風(fēng)面的雨滴。由于當(dāng)前基于歐拉-拉格朗日模型的計(jì)算方法需要待計(jì)算表面具有一定數(shù)量的雨滴以保證收集率計(jì)算精度，因此若要準(zhǔn)確計(jì)算建筑側(cè)面和背風(fēng)面的收集率分布，需要計(jì)算的雨滴數(shù)量比僅計(jì)算建筑迎風(fēng)面時(shí)所需雨滴數(shù)量更大，導(dǎo)致計(jì)算量巨大。

另外，在建筑迎風(fēng)面的收集率分布在邊界和角點(diǎn)處具有較大梯度，為準(zhǔn)確計(jì)算邊界、角點(diǎn)等區(qū)域的收集率，當(dāng)前方法需要在這些區(qū)域計(jì)算更多雨滴軌跡，以減小根據(jù)雨滴終點(diǎn)所劃分區(qū)域的面積，提高收集率計(jì)算精度。然而由于建筑繞流流場(chǎng)復(fù)雜，當(dāng)前的雨滴軌跡計(jì)算方法無法直接確定合適的雨滴初始條件來保證雨滴終點(diǎn)位于邊界和角點(diǎn)區(qū)域，只能通過在可能的范圍內(nèi)計(jì)算大量雨滴軌跡，以提高終點(diǎn)落入關(guān)心區(qū)域的雨滴軌跡數(shù)量。

由以上分析可見，若采用當(dāng)前基于歐拉-拉格朗日模型的計(jì)算方法研究建筑側(cè)面和背風(fēng)面以及邊界和角點(diǎn)區(qū)域的收集率分布，則需要計(jì)算巨量的雨滴軌跡，且難以保證收集率計(jì)算的準(zhǔn)確性。

技術(shù)實(shí)現(xiàn)要素：

本發(fā)明提出一種風(fēng)驅(qū)雨量的目標(biāo)軌跡追蹤計(jì)算方法，通過建立新的雨滴軌跡計(jì)算方法，實(shí)現(xiàn)對(duì)建筑側(cè)面、背風(fēng)面及邊界和角點(diǎn)區(qū)域收集率的準(zhǔn)確計(jì)算，同時(shí)所需的雨滴軌跡計(jì)算量相比現(xiàn)有方法顯著下降。

對(duì)于建筑側(cè)面、背風(fēng)面、邊界、角點(diǎn)等區(qū)域而言，采用現(xiàn)有方法計(jì)算雨滴軌跡難以保證有足夠數(shù)量的雨滴落在所關(guān)心區(qū)域內(nèi)。本發(fā)明提供的目標(biāo)軌跡追蹤計(jì)算方法則預(yù)先給定雨滴終點(diǎn)，并通過建立描述雨滴起點(diǎn)和終點(diǎn)位置關(guān)系的模型來計(jì)算恰好終止于給定終點(diǎn)的雨滴軌跡。采用本發(fā)明的方法，只需確定合理的雨滴終點(diǎn)位置，即能得到準(zhǔn)確的收集率分布。

雨滴在建筑繞流流場(chǎng)中的運(yùn)動(dòng)受自身重力以及阻力的作用。其中，阻力與流場(chǎng)風(fēng)速、雨滴本身運(yùn)動(dòng)速度、雨滴直徑等因素有關(guān)。由于流場(chǎng)風(fēng)速在空間的分布是連續(xù)的，雨滴運(yùn)動(dòng)速度隨時(shí)間也是連續(xù)變化，因此根據(jù)雨滴運(yùn)動(dòng)方程，相鄰時(shí)刻雨滴的坐標(biāo)滿足連續(xù)函數(shù)關(guān)系；又由于雨滴軌跡是對(duì)雨滴運(yùn)動(dòng)方程進(jìn)行時(shí)間積分得到，因此雨滴終點(diǎn)坐標(biāo)與起點(diǎn)坐標(biāo)也應(yīng)滿足一個(gè)連續(xù)的函數(shù)關(guān)系。然而，由于建筑繞流流場(chǎng)十分復(fù)雜，雨滴在運(yùn)動(dòng)過程中的受力隨位置變化而不斷變化，因此這些函數(shù)關(guān)系難以直接根據(jù)雨滴運(yùn)動(dòng)方程解析得到。由于雨滴運(yùn)動(dòng)過程的復(fù)雜性，針對(duì)上述問題，可采用代理模型方法建立描述雨滴起點(diǎn)坐標(biāo)和終點(diǎn)坐標(biāo)的函數(shù)關(guān)系。代理模型方法并不考慮雨滴的實(shí)際運(yùn)動(dòng)過程，而是根據(jù)一定數(shù)量樣本軌跡的起點(diǎn)和終點(diǎn)坐標(biāo)的對(duì)應(yīng)關(guān)系，建立函數(shù)表達(dá)式。在此基礎(chǔ)上，根據(jù)預(yù)先給定的雨滴終點(diǎn)坐標(biāo)即可給出對(duì)應(yīng)的起點(diǎn)坐標(biāo)，之后再根據(jù)雨滴運(yùn)動(dòng)方程計(jì)算得到雨滴軌跡，若軌跡終點(diǎn)與給定終點(diǎn)的誤差滿足精度要求，則認(rèn)為所得軌跡是有效的，否則對(duì)函數(shù)表達(dá)式進(jìn)行修正，直到所有雨滴軌跡均滿足精度要求為止。

在應(yīng)用代理模型方法時(shí)，需要一定數(shù)量的雨滴樣本軌跡。研究發(fā)現(xiàn)，當(dāng)給定終點(diǎn)在雨滴樣本軌跡的終點(diǎn)覆蓋范圍之內(nèi)時(shí)，采用所建立的函數(shù)關(guān)系預(yù)測(cè)得到的起點(diǎn)坐標(biāo)更為準(zhǔn)確。另外，在計(jì)算邊界、角點(diǎn)區(qū)域的收集率時(shí)，需要準(zhǔn)確計(jì)算終點(diǎn)位于邊界和角點(diǎn)上的雨滴樣本軌跡。因此，在給定雨滴樣本軌跡時(shí)，應(yīng)使得雨滴樣本軌跡的終點(diǎn)的覆蓋范圍超過待計(jì)算區(qū)域的范圍，從而使得邊界、角點(diǎn)上的點(diǎn)位于雨滴樣本軌跡的終點(diǎn)覆蓋范圍內(nèi)，以提高雨滴樣本軌跡計(jì)算的精度。為獲得超過待計(jì)算區(qū)域范圍的軌跡終點(diǎn)，需要將待計(jì)算區(qū)域進(jìn)行延展，建立一個(gè)空間中的虛平面，當(dāng)雨滴運(yùn)動(dòng)至該虛平面時(shí)，記錄當(dāng)前雨滴終點(diǎn)坐標(biāo)作為該雨滴軌跡的終點(diǎn)坐標(biāo)。將終止于虛平面和待計(jì)算區(qū)域的雨滴軌跡共同作為雨滴樣本軌跡，采用代理模型方法建立函數(shù)關(guān)系式，則對(duì)待計(jì)算區(qū)域內(nèi)部、邊界和角點(diǎn)處的雨滴軌跡均可進(jìn)行準(zhǔn)確計(jì)算，從而獲得準(zhǔn)確的收集率分布。

對(duì)于建筑側(cè)面和背風(fēng)面而言，只需存在可終止于待計(jì)算區(qū)域的雨滴軌跡，即可應(yīng)用目標(biāo)軌跡追蹤計(jì)算方法獲得該待計(jì)算區(qū)域內(nèi)其它給定終點(diǎn)的雨滴軌跡，進(jìn)而得到待計(jì)算區(qū)域內(nèi)的收集率分布。

本發(fā)明提供的風(fēng)驅(qū)雨量的目標(biāo)軌跡追蹤計(jì)算方法，分以下步驟：

步驟1：給定風(fēng)場(chǎng)來流條件，求解雷諾平均n-s方程獲得建筑繞流流場(chǎng)；

步驟2：根據(jù)待計(jì)算區(qū)域建立虛平面，并計(jì)算樣本軌跡在待計(jì)算區(qū)域或虛平面上的樣本軌跡終點(diǎn)；

步驟3：給定待計(jì)算區(qū)域上的雨滴軌跡終點(diǎn)坐標(biāo)；

步驟4：利用代理模型方法，根據(jù)樣本軌跡起點(diǎn)和樣本軌跡終點(diǎn)建立關(guān)于雨滴起點(diǎn)和終點(diǎn)坐標(biāo)關(guān)系的函數(shù)關(guān)系式；

步驟5：對(duì)步驟3給定的雨滴軌跡終點(diǎn)坐標(biāo)，利用步驟4所建立的函數(shù)關(guān)系式計(jì)算雨滴軌跡起點(diǎn)坐標(biāo)；

步驟6：根據(jù)步驟5計(jì)算得到的雨滴軌跡起點(diǎn)坐標(biāo)求解雨滴運(yùn)動(dòng)方程，將求解得到的預(yù)測(cè)軌跡和預(yù)測(cè)軌跡終點(diǎn)坐標(biāo)與之前給定的雨滴軌跡終點(diǎn)坐標(biāo)進(jìn)行對(duì)比，若誤差滿足精度要求，則記錄該預(yù)測(cè)軌跡為有效軌跡，否則將該預(yù)測(cè)軌跡加入樣本軌跡，返回步驟4。

步驟7：重復(fù)步驟4～6，得到所有直徑雨滴在給定雨滴軌跡終點(diǎn)的有效軌跡，然后根據(jù)有效軌跡計(jì)算該待計(jì)算區(qū)域的收集率分布。

本發(fā)明的優(yōu)點(diǎn)在于：

1.在現(xiàn)有技術(shù)中的基于歐拉-拉格朗日模型的雨滴軌跡計(jì)算方法中，雨滴軌跡是根據(jù)給定的初始條件積分雨滴運(yùn)動(dòng)方程得到的，在計(jì)算結(jié)束前無法確定雨滴的終點(diǎn)位置，這就導(dǎo)致有大量雨滴軌跡最終未終止于待計(jì)算的表面，而這些軌跡對(duì)收集率計(jì)算而言是無效的。而在本發(fā)明提供的目標(biāo)軌跡追蹤計(jì)算方法中，則是先給定雨滴軌跡終點(diǎn)，再根據(jù)該終點(diǎn)確定相應(yīng)的初始條件。這樣就能保證所計(jì)算的雨滴軌跡均可用于所關(guān)心區(qū)域的收集率計(jì)算，而不需額外計(jì)算大量無用的雨滴軌跡，因而可有效降低雨滴軌跡計(jì)算量。

2.本發(fā)明的目標(biāo)軌跡追蹤計(jì)算方法預(yù)先給定了待計(jì)算區(qū)域的終點(diǎn)位置，因此可根據(jù)流動(dòng)條件和物面幾何形狀特征等在不同區(qū)域給定不同密度的軌跡終點(diǎn)分布，即對(duì)待計(jì)算區(qū)域進(jìn)行加密，從而更準(zhǔn)確地獲得該區(qū)域的收集率分布。

3.本發(fā)明的目標(biāo)軌跡追蹤計(jì)算方法能夠快速獲得建筑側(cè)面、背風(fēng)面等區(qū)域的收集率分布，并能通過區(qū)域加密，準(zhǔn)確計(jì)算邊界、角點(diǎn)等梯度較大區(qū)域的收集率分布。

附圖說明

圖1為本發(fā)明的目標(biāo)軌跡追蹤技術(shù)方法中計(jì)算給定終點(diǎn)軌跡示意圖；

圖2為三維方形建筑計(jì)算區(qū)域示意圖；

圖3為三維方形建筑虛平面示意圖；

圖4為現(xiàn)有參考文獻(xiàn)給出的三維方形建筑迎風(fēng)面上半?yún)^(qū)收集率等值線圖；

圖5為采用本發(fā)明的目標(biāo)軌跡追蹤計(jì)算方法得到的三維方形建筑迎風(fēng)面上半?yún)^(qū)收集率等值線圖。

圖中：

1.樣本軌跡釋放平面；2.待計(jì)算區(qū)域；3.樣本軌跡起點(diǎn)；4.樣本軌跡；5.樣本軌跡終點(diǎn)；6.雨滴軌跡終點(diǎn)；7.雨滴軌跡起點(diǎn)；8.預(yù)測(cè)軌跡；9.預(yù)測(cè)軌跡終點(diǎn)；10.計(jì)算區(qū)域入口；11.方形建筑迎風(fēng)面；12.虛平面；13.虛平面上的樣本軌跡終點(diǎn)；14.待計(jì)算平面上的目標(biāo)軌跡終點(diǎn)。

具體實(shí)施方式

下面結(jié)合附圖和實(shí)施例對(duì)本發(fā)明進(jìn)行詳細(xì)說明。

本發(fā)明提供一種風(fēng)驅(qū)雨量的目標(biāo)軌跡追蹤計(jì)算方法，在待計(jì)算區(qū)域預(yù)先給定雨滴軌跡終點(diǎn)，通過建立描述雨滴起點(diǎn)和雨滴終點(diǎn)坐標(biāo)關(guān)系的函數(shù)表達(dá)式，獲得給定雨滴軌跡終點(diǎn)對(duì)應(yīng)的雨滴軌跡起點(diǎn)，再根據(jù)雨滴軌跡起點(diǎn)求解雨滴運(yùn)動(dòng)方程得到預(yù)測(cè)軌跡和預(yù)測(cè)軌跡終點(diǎn)，最后根據(jù)預(yù)測(cè)軌跡計(jì)算相應(yīng)的待計(jì)算區(qū)域的收集率。在建立函數(shù)關(guān)系式時(shí)，采用了代理模型方法，所需雨滴樣本軌跡數(shù)量遠(yuǎn)小于現(xiàn)有技術(shù)中的基于歐拉-拉格朗日模型方法所需計(jì)算的雨滴數(shù)量，因此本發(fā)明的目標(biāo)軌跡追蹤計(jì)算方法相比現(xiàn)有方法可有效降低計(jì)算量。同時(shí)，本發(fā)明的目標(biāo)軌跡追蹤計(jì)算方法能準(zhǔn)確獲得建筑側(cè)面、背風(fēng)面、邊界、角點(diǎn)等區(qū)域的收集率分布，而現(xiàn)有方法難以給出準(zhǔn)確結(jié)果。

本發(fā)明提供的風(fēng)驅(qū)雨量的目標(biāo)軌跡追蹤計(jì)算方法，結(jié)合圖1和圖3，具體包括以下步驟：

步驟1：建立建筑繞流計(jì)算網(wǎng)格，給定風(fēng)場(chǎng)來流條件(所述的來流條件是指來流速度或速度型)，求解雷諾平均n-s方程，獲得建筑繞流流場(chǎng)。

步驟2：在建筑上選取待計(jì)算區(qū)域2，根據(jù)待計(jì)算區(qū)域2建立虛平面12。

如圖3所示，當(dāng)待計(jì)算區(qū)域2為單一平面時(shí)，直接將該待計(jì)算區(qū)域2進(jìn)行平面延展，得到與待計(jì)算區(qū)域2共面且包含待計(jì)算區(qū)域2的無限大空間平面，作為虛平面12；當(dāng)待計(jì)算區(qū)域2由多個(gè)平面組成時(shí)，分別延展各個(gè)平面，獲得不同平面對(duì)應(yīng)的虛平面12；當(dāng)待計(jì)算區(qū)域2為曲面時(shí)，需根據(jù)待計(jì)算區(qū)域2上給定的雨滴軌跡終點(diǎn)6的分布，分別對(duì)每個(gè)給定的雨滴軌跡終點(diǎn)6處建立曲面的切平面，作為該給定的雨滴軌跡終點(diǎn)6對(duì)應(yīng)的虛平面12。

根據(jù)步驟1中計(jì)算得到的建筑繞流流場(chǎng)，確定樣本軌跡釋放平面1的高度，在該高度處流場(chǎng)風(fēng)速與來流相同。根據(jù)待計(jì)算區(qū)域2的范圍，確定樣本軌跡釋放平面1的范圍，應(yīng)使得從該樣本軌跡釋放平面1釋放的雨滴，沿樣本軌跡在待計(jì)算區(qū)域2對(duì)應(yīng)虛平面12上的樣本軌跡終點(diǎn)所覆蓋范圍包含待計(jì)算區(qū)域2。

在樣本軌跡釋放平面1上均勻分布樣本軌跡起點(diǎn)3。確定待計(jì)算的樣本雨滴直徑，在釋放平面1內(nèi)給定各直徑樣本雨滴初始位置(樣本軌跡起點(diǎn)3坐標(biāo))和初始速度。其中，初始速度的水平分量與當(dāng)?shù)仫L(fēng)場(chǎng)速度的水平分量相等，豎直分量設(shè)為該直徑樣本雨滴的下落終速。針對(duì)每一個(gè)樣本軌跡起點(diǎn)3，求解雨滴運(yùn)動(dòng)方程，方程形式如下：

其中，ρp、ρa(bǔ)分別為水和空氣的密度，x、y、z為雨滴坐標(biāo)，u、v、w分別為x、y、z方向的風(fēng)速，up、vp、wp分別為雨滴在x、y、z方向的速度，ψ、γ分別為雨滴與風(fēng)場(chǎng)相對(duì)速度與x、y、z方向的夾角，μ為空氣粘性系數(shù)，re為雨滴相對(duì)于風(fēng)場(chǎng)速度的雷諾數(shù)，g為重力加速度，cd為雨滴阻力系數(shù)，d為雨滴直徑，t為時(shí)間。

根據(jù)式(1)計(jì)算樣本軌跡4，獲得在待計(jì)算區(qū)域2或虛平面12上的樣本軌跡終點(diǎn)5。

步驟3：在待計(jì)算區(qū)域2內(nèi)給定雨滴軌跡終點(diǎn)6的坐標(biāo)，為保證收集率計(jì)算精度，應(yīng)使得待計(jì)算區(qū)域2的邊界和角點(diǎn)區(qū)域終點(diǎn)分布較密。

步驟4：選擇代理模型方法，將步驟2中的樣本軌跡起點(diǎn)3和樣本軌跡終點(diǎn)5的坐標(biāo)代入代理模型，計(jì)算代理模型中的系數(shù)，建立描述雨滴起點(diǎn)和終點(diǎn)的坐標(biāo)關(guān)系的函數(shù)關(guān)系式。在本步驟中，應(yīng)針對(duì)不同直徑樣本雨滴分別建立函數(shù)關(guān)系式。

步驟5：將步驟3中給定的雨滴軌跡終點(diǎn)6的坐標(biāo)帶入步驟4中所建立的函數(shù)關(guān)系式，計(jì)算得到相應(yīng)直徑雨滴的雨滴軌跡起點(diǎn)7的坐標(biāo)。

步驟6：將步驟5計(jì)算得到的雨滴軌跡起點(diǎn)7的坐標(biāo)作為雨滴初始位置，將相同直徑樣本雨滴的初始速度作為初始速度，求解雨滴運(yùn)動(dòng)方程得到雨滴運(yùn)動(dòng)的預(yù)測(cè)軌跡8和預(yù)測(cè)軌跡終點(diǎn)9，并計(jì)算該預(yù)測(cè)軌跡終點(diǎn)9與給定的雨滴軌跡終點(diǎn)6的坐標(biāo)的誤差，若滿足精度要求則記錄該預(yù)測(cè)軌跡8為有效軌跡。否則，將該預(yù)測(cè)軌跡8作為樣本軌跡，返回步驟4。

步驟7：重復(fù)步驟4～6，得到所有直徑雨滴在給定雨滴軌跡終點(diǎn)的有效軌跡，然后根據(jù)有效軌跡計(jì)算該待計(jì)算區(qū)域的收集率分布。

所述的代理模型方法中的代理模型可根據(jù)實(shí)際計(jì)算問題選擇多項(xiàng)式響應(yīng)面模型、徑向基函數(shù)模型、kriging模型等，本發(fā)明實(shí)施例以徑向基函數(shù)模型為例進(jìn)行計(jì)算。

實(shí)施例：如圖2所示的三維方形建筑迎風(fēng)面收集率計(jì)算。建筑長(zhǎng)、寬、高均為10m。計(jì)算域入口邊界10給定來流速度型，表達(dá)式如下：

其中，u代表當(dāng)?shù)厮斤L(fēng)速；h代表當(dāng)?shù)馗叨?；u0代表參考速度，取10m/s；h0代表參考高度，取10m；α為指數(shù)，取0.15。

采用徑向基函數(shù)模型建立雨滴起點(diǎn)坐標(biāo)和終點(diǎn)坐標(biāo)的函數(shù)關(guān)系式。徑向基函數(shù)模型的一般形式為：

其中，yk為k點(diǎn)的輸出值；wi為系數(shù)；n為樣本點(diǎn)數(shù)；rik為i點(diǎn)與k點(diǎn)的歐式距離；為基函數(shù)。在本實(shí)施例中，基函數(shù)選擇multiquadric函數(shù)，其形式為：

其中，c為常數(shù)。本研究中取c為樣本點(diǎn)距離函數(shù)值的標(biāo)準(zhǔn)差。

應(yīng)用本發(fā)明提供的目標(biāo)軌跡追蹤計(jì)算方法和參考文獻(xiàn)1(參考文獻(xiàn)1：b.blocken,j.carmeliet.onthevalidityofthecosineprojectioninwind-drivenraincalculationsonbuildings[j].buildingandenvironment,2006,41:1182-1189.)提供的方法分別計(jì)算所述三維方形建筑迎風(fēng)面收集率，結(jié)果如圖5和圖4所示，圖中給出了迎風(fēng)面11收集率等值線圖，降雨強(qiáng)度為1mm/h?？梢钥吹剑繕?biāo)軌跡追蹤計(jì)算方法準(zhǔn)確得到了平面上的收集率分布，尤其是在平面的邊界和角點(diǎn)區(qū)域，目標(biāo)軌跡追蹤計(jì)算方法與參考文獻(xiàn)結(jié)果十分接近；但目標(biāo)軌跡追蹤計(jì)算方法所需計(jì)算量不到參考文獻(xiàn)中方法的10％?？梢?，本發(fā)明能夠顯著降低計(jì)算量，同時(shí)還可獲得準(zhǔn)確的收集率分布結(jié)果。

本發(fā)明的上述實(shí)施方案只是舉例說明，并不是僅有的。本發(fā)明中收集率計(jì)算平面的虛平面以及代理模型形式根據(jù)實(shí)際研究對(duì)象和條件確定，并適用于任何外形和來流及降雨條件下的風(fēng)驅(qū)雨量計(jì)算，所有在本發(fā)明范圍內(nèi)或等同于本發(fā)明的范圍內(nèi)的改變均被本發(fā)明包含。