一种用于噪声图绘制的城市广场声传播的辐射度建模方法技术

一种用于噪声图绘制的城市广场声传播的辐射度建模方法，属于广场型城市空间的噪声图绘制领域。本发明专利技术是为了解决现有声传播模型的建模方法在对建筑群进行建模时由于忽略了建筑群的边界面和地面的不规则性而导致的建模精度低的问题。首先在城市广场中布置点声源和接收点，然后对城市广场进行三维立体建模并对建模后的城市广场的每一个面进行划分，点声源发出声能并分配给每一个面元，获得1阶面元声源，再根据1阶面元声源及形状因子获得k阶面元声源及k阶面元声源的能量，接收点根据接收的k阶面元声源的能量的大小绘制辐射度模型。本发明专利技术适用于广场型城市空间声传播的计算。

【技术实现步骤摘要】
一种用于噪声图绘制的城市广场声传播的辐射度建模方法
本专利技术属于广场型城市空间的噪声图绘制领域。
技术介绍
目前噪声图的绘制是欧洲国家的法规要求，我国也正在积极推进城市噪声图的绘制，但目前噪声图的计算多数基于经验公式，通常的方法是声线法或虚声源法，这两种方法在对建筑进行建模计算时，通常忽略建筑边界面和地面的不规则性，也就是说这两种方法只适用于光滑的表面。但在实际情况下，城市广场、街道或建筑群的边界面和地面总是存在一些不规则性。在利用现有方法对建筑群进行建模时，由于忽略了建筑群的边界面和地面的不规则性，导致现有方法的建模精度很低，采用现有方法建立的声传播模型精度仅为8～10dB，因此研发新的建模方法提高建模精度是非常必要的。
技术实现思路
本专利技术是为了解决现有声传播模型的建模方法在对建筑群进行建模时由于忽略了建筑群的边界面和地面的不规则性而导致的建模精度低的问题。现提供计算城市广场声传播的辐射度模型的建模方法。一种用于噪声图绘制的城市广场声传播的辐射度建模方法，该方法包括下述步骤：步骤一、在城市广场中布置点声源S和接收点R；步骤二、对城市广场进行三维立体建模，并对构成城市广场的每一个面划分面元；步骤三、点声源S发出声能分配给步骤二中的面元，并获得1阶面元声源；步骤四、根据1阶面元声源及形状因子获得k阶面元声源及k阶面元声源的能量；步骤五、接收点R接收k阶面元声源的能量并获得接收点R的能量；步骤六、接收点R根据接收的k阶面元声源的能量的大小绘制辐射度模型。本专利技术所述的一种用于噪声图绘制的城市广场声传播的辐射度建模方法，在对城市广场的各个面进行面元划分的时候，将各个面划分成若干个面元，并且越靠近边缘或者城市广场某一面不规则时，面元划分的越小，然后再进行将声源的声能分配给面元、求形状因子、面元之间的能量交换及能量从面元到接收点的传递，得到的形状因子及能量分配精确到小数三位，从而提高了建模精度，相比现有的声传播模型的建模方法，精度提高了30％以上。本专利技术适用于广场型城市空间声传播的计算。附图说明图1为三面围合的城市广场面元划分示意图；图2为五面围合的城市广场面元划分示意图；图3为具体实施方式三中从点源到达地面上一个面元的能量分布图；图4为具体实施方式四中从发射面元到垂直接收面元的形状因子的确定图；图5为具体实施方式七的实施例绘制出的城市广场的噪声图；其中，图1中A表示A面，B表示B面，G表示地面，H表示高，W表示宽，L表示长；图2中A表示A面，B表示B面，C表示C面，D表示D面，S表示点声源，R表示接收点；图4中“1”表示单位圆的半径是1；图5中“2”表示墙体，“3”表示绿地，“4”表示玻璃。具体实施方式具体实施方式一：参照图1和图2具体说明本实施方式，本实施方式所述的一种用于噪声图绘制的城市广场声传播的辐射度建模方法，该方法包括下述步骤：步骤一、在城市广场中布置点声源S和接收点R；步骤二、对城市广场进行三维立体建模，并对构成城市广场的每一个面划分面元；步骤三、点声源S发出声能分配给步骤二中的面元，并获得1阶面元声源；步骤四、根据1阶面元声源及形状因子获得k阶面元声源及k阶面元声源的能量；步骤五、接收点R接收k阶面元声源的能量并获得接收点R的能量；步骤六、接收点R根据接收的k阶面元声源的能量的大小绘制辐射度模型。具体实施方式二：本实施方式是对具体实施方式一所述的一种用于噪声图绘制的城市广场声传播的辐射度建模方法的进一步说明，本实施方式中，步骤二中所述的对城市广场进行三维立体建模，并对构成城市广场的每一个面划分面元的方法如下：建立空间坐标系，分别为x轴，y轴和z轴，将长为L、宽为W、高为H的城市广场所在的面分别投射到x＝0、z＝0、y＝0的面；沿长度方向划分为l个面元，l＝1...Nx，Nx为正整数，且Nx能被4整除；沿宽度方向划分为m个面元，m＝1...Ny，Ny为正整数；沿高度方向划分为n个面元，n＝1...Nz，Nz为正整数，且越靠近边缘处，划分的面元越小；如果沿某个比如沿宽度方向的尺寸小于50m时，那么该方向的面元尺寸ddm的划分方法为：式中qy(qy>1)是沿宽度方向的两个相邻面元之间的比值，Ny为偶数，km为：如果沿某个比如沿长度方向的尺寸大于50m时，那么该方向的面元尺寸ddl的划分方法为：式中qx(qx>1)是沿长度方向的两个相邻面元的比值，Nx是4的整倍数，kl为：具体实施方式三：参照图3具体说明本实施方式，本实施方式是对具体实施方式一所述的一种用于噪声图绘制的城市广场声传播的辐射度建模方法的进一步说明，本实施方式中，步骤三中所述的点声源S发出声能分配给步骤二中的面元，并获得1阶面元声源的方法如下：设点声源S的坐标为(Sx,Sy,Sz)；面元G1(t)l,m计算如下：其中，t＝Sl,m/c，K是与声源功率有关的常数，αGl,m是面元Gl,m的与入射角无关的吸收系数，吸收系数为1，Sl,m是声源到面元的平均声束程长，近似地取为声源到面元中心的距离；和βl,m是面元Gl,m对于声源的位置决定的角度，分别计算如下：在式(b)中，用于考虑声源位于两个面元Gl,m之间的情况，当dl-ddl/2≤Sx≤dl+ddl/2时，否则，同理，在式(c)中，当dl-ddl/2≤Sy≤dl+ddl/2时，kβ＝-1，否则，kβ＝1；如果声源有指向性，则式(d)应加上表示声源在Gl,m方向的辐射强度；此时，面元被视为声源，即1阶面元声源。声源分配的基本原则是，每个面元的声能的百分比，等于该面元的立体角与点声源总的立体角的比值；将面元视为声源，即获得1阶面元声源。具体实施方式四：参照图4具体说明本实施方式，本实施方式是对具体实施方式一所述的一种用于噪声图绘制的城市广场声传播的辐射度建模方法的进一步说明，本实施方式中，步骤四中所述的根据1阶面元声源及形状因子获得k阶面元声源的方法如下：步骤A、求形状因子：对于正交的两个面元，形状因子采用Nusselt方法计算；Nusselt方法中形状因子等于，将接收面元投射到以辐射面元中心为圆心的单位半球上，将投射部分的面积正交投射到半球的单位圆上，得到的投影面积除以圆的面积；对于平行的两个面元，形状因子采用Cohen和Greenberg的方法计算，形状因子的计算方法为，将接收面元投射到辐射面元中心立方体的上半部分；步骤B、每个1阶面元声源接收点声源S的声能，并将点声源S与1阶面元声源之间的形状因子再分配给其它面元，形成二阶面元声源，二阶面元声源将点声源S与二阶面元声源之间的形状因子再分配给其它面元，形成三阶面元声源，即将k-1阶面元声源的声能分配给其他面元，形成k阶面元声源，其中k＝1...∞。形状因子的获得方法：在图1所示的城市广场中，结合图4说明从辐射面元Al',n'(l'＝1,2,3,...,Nx)到接收面元Gl,m的计算图。考虑了面元Al',n'吸收点声源S发出的声能及空气吸收点声源S发出的声能后，从Al',n'辐射到Gl,m的能量AG(l',n'),(l,m)计算如下：式中，d(l',n'),(l,m)是面元Al',n'到面元Gl,m的平均程长，用两个面元中心之间的距离近似。γ(l',n'),(l,m)、Δγ(l',n'),(l,m)和是决定两个面元相对位置的本文档来自技高网...

【技术保护点】
一种用于噪声图绘制的城市广场声传播的辐射度建模方法，其特征在于，该方法包括下述步骤：步骤一、在城市广场中布置点声源S和接收点R；步骤二、对城市广场进行三维立体建模，并对构成城市广场的每一个面划分面元；步骤三、点声源S发出声能分配给步骤二中的面元，并获得1阶面元声源；步骤四、根据1阶面元声源及形状因子获得k阶面元声源及k阶面元声源的能量；步骤五、接收点R接收k阶面元声源的能量并获得接收点R的能量；步骤六、接收点R根据接收的k阶面元声源的能量的大小绘制辐射度模型。

【技术特征摘要】
1.一种用于噪声图绘制的城市广场声传播的辐射度建模方法，该方法包括下述步骤：步骤一、在城市广场中布置点声源S和接收点R；步骤二、对城市广场进行三维立体建模，并对构成城市广场的每一个面划分面元；步骤三、点声源S发出声能分配给步骤二中的面元，并获得1阶面元声源；步骤四、根据1阶面元声源及形状因子获得k阶面元声源及k阶面元声源的能量；步骤五、接收点R接收k阶面元声源的能量并获得接收点R的能量；步骤六、接收点R根据接收的k阶面元声源的能量的大小绘制辐射度模型；其特征在于，步骤三中所述的点声源S发出声能分配给步骤二中的面元，并获得1阶面元声源的方法如下：设点声源S的坐标为(Sx,Sy,Sz)；面元G1(t)l,m计算如下：其中，t＝Sl,m/c，K是与声源功率有关的常数，αGl,m是面元Gl,m的与入射角无关的吸收系数，吸收系数为1，Sl,m是声源到面元的平均声束程长，近似地取为声源到面元中心的距离；和βl,m是面元Gl,m对于声源的位置决定的角度，分别计算如下：在式(b)中，用于考虑声源位于两个面元Gl,m之间的情况，当dl-ddl/2≤Sx≤dl+ddl/2时，否则，同理，在式(c)中，当dl-ddl/2≤Sy≤dl+ddl/2时，kβ＝-1，否则，kβ＝1；如果声源有指向性，则式(d)应加上表示声源在Gl,m方向的辐射强度；此时，面元被视为声源，即1阶面元声源。2.根据权利要求1所述的一种用于噪声图绘制的城市广场声传播的辐射度建模方法，其特征在于，步骤二中所述的对城市广场进行三维立体建模，并对构成城市广场的每一个面划分面元的方法如下：建立空间坐标系，分别为x轴，y轴和z轴，将长为L、宽为W、高为H的城市广场所在的面分别投射到x＝0、z＝0、y＝0的面；沿长度方向划分为l个面元，l＝1...Nx，Nx为正整数，且Nx能被4整除；沿宽度方向划分为m个面元，m＝1...Ny，Ny为正整数；沿高度方向划分为n个面元，n＝1...Nz，Nz为正整数，且越靠近边缘处，划分的面元越小；如果沿某个比如沿宽度方向的尺寸小于50m时，那么该方向的面元尺寸ddm的划分方法为：式中qy(qy>1)是沿宽度方向的两个相邻面元之间的比值，Ny为偶数，km为：如果沿某个比如沿长度方向的尺寸大于50m时，那么该方向的面元尺寸ddl的划分方法为：式中qx(qx>1)是沿长度方向的两个相邻面元的比值，Nx是4的整倍数，kl为...

【专利技术属性】
技术研发人员：康健，
申请(专利权)人：哈尔滨工业大学，
类型：发明
国别省市：黑龙江;23