本发明提出一种线源等效的不规则形状面声源几何发散计算方法。方法中，把面声源对接收点的声音辐射能量等价于面声源的各个轮廓线对该接收点的声音辐射能量之和。首先根据声源总声功率相等的原则，将面源总声功率分配到面源的边界轮廓线上；然后针对任一线声源对接收点进行声能计算，整个线声源在接收点的声功率为单位长度源强下接收点的声功率在线源长度上的积分；最后将各轮廓线对接收点的声能进行叠加，得到整个面源对接收点的声音辐射。本发明应用于不规则形状面声源声音几何发散的计算，解决了复杂边界面声源几何发散计算中源强处理复杂、计算精度过低和适用性差的问题。
1.一种线源等效的不规则面声源几何发散计算方法，其特征在于，所述方法包含以下步骤：(a)获取单位面声源强度、面声源面积、面声源轮廓线周长及面声源总功率参数，基于声源总功率相等原则，将面声源总功率分配到面声源的边界轮廓线上；
(b)根据面声源的几何形状，对面声源的边界轮廓线进行离散化，得到若干直线段轮廓线，并将每条直线段轮廓线均视为有限长线声源；
(c)获取接收点和声源之间的环境参数，所述环境参数包括：单位面声源强度、单位等效线声源强度、面声源的面积、面声源轮廓线的周长、线声源的两个端点、接收点到线声源两端点的夹角、接收点到线声源的垂直距离、离散后线声源的个数和线声源的长度；
(d)基于获取的环境参数和声功率积分方法，计算每一条线声源在接收点处的声音辐射能量；
(e)通过每一条线声源在接收点处的声音辐射能量叠加计算得到接收点处的面声源对接收点的声音辐射能量；
所述步骤(a)中，声源总功率相等的原则，能表示为面声源总功率等于各轮廓线声功率总和，具体为： 式中，StrengthA表示单位面声源强度，StrengthL表示单位等效线声源强度，AA表示面声源的面积，CL表示面声源轮廓线的周长；
所述步骤(a)中，将面声源总功率分配到面声源的边界轮廓线上，单位等效线声源强度具体为： 式中，StrengthA表示单位面声源强度，StrengthL表示单位等效线声源强度，AA表示面声源的面积，CL表示面声源轮廓线的周长。
2.根据权利要求1所述的一种线源等效的不规则面声源几何发散计算方法，其特征在于，面声源对接收点的声音辐射能量等价于面声源各个轮廓线对该接收点的声音辐射能量之和。
3.根据权利要求1所述的一种线源等效的不规则面声源几何发散计算方法，其特征在于，所述步骤(b)中，面声源的边界轮廓线，包括直线段线和弧形线，具体离散化过程为：对于直线段线，直接形成对应的直线段线声源；对于弧形线，形成能够用于计算的多条有限长直线段线声源。
4.根据权利要求1所述的一种线源等效的不规则面声源几何发散计算方法，其特征在于，所述步骤(d)中，计算每一条线声源对接收点的声音辐射能量值具体为：LL＝10lg(IL)，式中，LL为一条线声源对接收点的声音辐射能量值，IL为此条线声源对接收点的声功率贡献量。
5.根据权利要求4所述的一种线源等效的不规则面声源几何发散计算方法，其特征在于，所述此条线声源对接收点的声功率贡献量，采用声功率积分的方法进行计算，具体为：式中，IL为此条线声源对接收点的声功率贡献量，x1、x2分别
表示线声源的两个端点，Ix表示任意单位线源对接收点的声功率贡献，StrengthL表示单位等效线声源强度，θ为接收点到线源两端点的夹角，d为接收点到线声源的垂直距离。
6.根据权利要求5所述的一种线源等效的不规则面声源几何发散计算方法，其特征在于，所述任意单位线源对接收点的声功率贡献，具体为： 式中，Ix表示任意单位线源对接收点的噪声声功率贡献，StrengthL表示单位等效线声源强度，x表示线声源的长度，d为接收点到线声源的垂直距离。
7.根据权利要求1所述的一种线源等效的不规则面声源几何发散计算方法，其特征在于，所述步骤(e)中，面声源对接收点的声音辐射能量，通过各轮廓线对接收点的声能叠加得到，具体为： 式中，LA表示面声源对接收点的声音辐射能量，LLi为第i个线声源对接收点的声能贡献值，StrengthL表示单位等效线声源强度，n表示离散后线声源的个数，θi为接收点到第i个线源两端点的夹角，di为接收点到第i个线声源的垂直距离。
一种线源等效的不规则面声源几何发散计算方法
技术领域
[0001] 本发明涉及面声源声音排放和平面几何技术领域，更具体地，涉及一种线源等效的不规则面声源几何发散计算方法。
背景技术
[0002] 近年来，随着人们环境意识的增加，城市声环境质量成为人们越发关注的问题。声源的振动形成声波并向外传播，声源的简化、特性和声场在空间的分布随着声源类型的不同而有很大的差异，因此，不同类型声源声音在几何空间中的发散规律是学者们在研究和工程上都必须要面对的问题。实际生活中声源的形状复杂多样，很难通过严格的控制按照数学方法对复杂边界声源的声场进行精确的计算，因此，在实际工程应用中，通常根据声源的形状、接收点和声源的相对位置以及声音本身的特性对声源进行简化。
[0003] 对于面声源的声衰减计算，在HJ2.4-2009《环境影响评价技术导则声环境》中，对尺寸为a×b长方形(b>a)的面声源声衰减计算做出了规定，但导则的面声源形状局限于矩形，适用性较差。其他面声源声几何发散方法包括了重心替代法、网格化点源替代法和面源单元积分法等。在重心替代法中，面声源用位于面源几何重心的点声源替代，这种方法虽然极大简化了计算过程，但计算结果过于粗糙，不能适用于接收点距离面声源较近的情况。对于网格化点源替代法，则是把面声源网格化，用一系列位于网格中心的点声源阵列来代替面声源，这种方法适用于声强分布平均的复杂形状的面声源声发散计算，但当面声源各部分声辐射功率差异较大时计算出现较大模拟误差，且该方法的声源处理和声发散计算复杂。而面源单元积分法则通过精确积分法、数值积分法或有限分割法等竖直方法进行求解，这类方法以任意点为研究对象，通过单元积分的形式得到面源的辐射情况，方法需要大量的计算，而且在遇到不规则形状的面声源时，需要结合普适性公式等近似方法来克服精确解无法求得的情况。
发明内容
[0004] 本发明目的在于克服现有技术的缺点与不足，提供一种计算精度高、处理方法简单的一种线源等效的不规则面声源几何发散计算方法，用于不规则几何边界的面声源对接收点的声能贡献计算。
[0005] 为解决上述技术问题，本发明采用的技术方案是：
[0006] 一种线源等效的不规则面声源几何发散计算方法，包含以下步骤：
[0007] (a)获取单位面声源强度、面声源面积、面声源轮廓线周长及面声源总功率等参数，基于声源总功率相等原则，将面声源总功率分配到面声源的边界轮廓线上；
[0008] (b)根据面声源的几何形状，对面声源的边界轮廓线进行离散化，得到若干直线段轮廓线，并将每条直线段轮廓线均视为有限长线声源；
[0009] (c)获取接收点和声源之间的环境参数，所述环境参数包括：单位面声源强度、单位等效线声源强度、面声源的面积、面声源轮廓线的周长、线声源的两个端点、接收点到线声源两端点的夹角、接收点到线声源的垂直距离、离散后线声源的个数和线声源的长度；
[0010] (d)基于获取的环境参数和声功率积分方法，计算每一条线声源在接收点处的声音辐射能量(声能贡献)；
[0011] (e)通过每一条线声源在接收点处的声音辐射能量(声能贡献)叠加计算得到接收点处的面声源对接收点的声音辐射能量(声能贡献)。
[0012] 优选地，所述面声源对接收点的声音辐射能量等价于面声源各个轮廓线对该接收点的声音辐射能量之和。
[0013] 优选地，所述声源总功率相等的原则，可表示为面源总声功率等于各轮廓线声功率总和，具体为： 式中，StrengthA表示单位面声源强度，StrengthL表示单位等效线声源强度，AA表示面声源的面积，CL表示面声源轮廓线的周长。
[0014] 优选地，所述步骤(a)中，将面声源总功率分配到面声源的边界轮廓线上，单位等效线声源强度具体为： 式中，StrengthA表示单位面声源强度，StrengthL表示单位等效线声源强度，AA表示面声源的面积，CL表示面声源轮廓线的周长。
[0015] 优选地，所述步骤(b)中，面声源的边界轮廓线，包括直线段线和弧形线，具体离散化过程为：对于直线段线，直接形成对应的直线段线声源；对于弧形线，形成能够用于计算的多条有限长直线段线声源。
[0016] 优选地，所述步骤(d)中，计算每一条线声源对接收点的声音辐射能量值具体为：LL＝10lg(IL)，式中，LL为一条线声源对接收点的声音辐射能量值，IL为此条线声源对接收点的声功率贡献量。
[0017] 优选地，所述此条线声源对接收点的声功率贡献量，采用声功率积分的方法进行计算，具体为： 式中，IL为此条线声源对接收点的声功率贡献量，x1、x2分别表示线声源的两个端点，Ix表示任意单位线源对接收点的声功率贡献，StrengthL表示单位等效线声源强度，θ为接收点到线源两端点的夹角，d为接收点到线声源的垂直距离。
[0018] 优选地，所述任意单位线源对接收点的声功率贡献，具体为：式中，Ix表示任意单位线源对接收点的噪声声功率贡献，StrengthL表示单位等效线声源强度，x表示线声源的长度，d为接收点到线声源的垂直距离。
[0019] 优选地，所述步骤(e)中，面声源对接收点的声音辐射能量，通过各轮廓线对接收点的声能叠加得到，具体为： 式中，LA表示面声源对接收点的声音辐射能量，LLi为第i个线声源对接收点的声能贡献值，StrengthL表示单位等效线声源强度，n表示离散后线声源的个数，θi为接收点到第i个线源两端点的夹角，di为接收点到第i个线声源的垂直距离。
[0020] 与现有技术相比，本发明具有以下优点：
[0021] 本方法从考虑面声源几何形状和面积出发，通过对面声源声能辐射的等效替代，把面声源声能计算转化为易操作且有成熟计算方法的有限长线声源计算。本方法能够应用于不规则形状面声源声音几何发散计算，有效的解决了工程应用中不规则面声源声音几何发散计算精度低、处理复杂和适用性差等问题，其具体表现为：
[0022] (1)考虑面声源几何形状的计算方法更加精确合理。相对于重心替代法，本发明充分考虑了面声源的几何形状，而不是简单的用面声源质心作点声源替代，能够精确的计算复杂边界对面声源声能发散的影响，更贴合实际情况，计算精度更高；
[0023] (2)本发明在计算面声源几何发散过程中，只需要对面声源的几何边界进行离散化，离散的单元即直线段边界单元，而不需要对面声源进行网格化分割，对源强的处理简单，能大大提高面声源声发散计算的效率。
附图说明
[0024] 图1是本发明流程示意图。
[0025] 图2是本发明中面源离散化示意图。
[0026] 图3是本发明中等效线声源和接收点示意图。
具体实施方式
[0027] 下面结合附图和具体实施方式对本发明作进一步的详细说明，应当理解，此处所描述的具体实施例仅仅用于解释本发明，并不用于限定本发明。
[0028] 如图1所示，本发明提供一种线源等效的不规则面声源几何发散计算方法。通过以下步骤实现：
[0029] S100.获取单位面声源强度、面声源面积、面声源轮廓线周长及面声源总功率参数，基于声源总功率相等原则，将面声源总功率分配到面声源的边界轮廓线上；
[0030] S110.根据面声源的几何形状，对面声源的边界轮廓线进行离散化，得到若干直线段轮廓线，并将每条直线段轮廓线均视为有限长线声源；
[0031] 面声源的边界轮廓，可以包括直线段和弧形，具体离散化过程为：对于直线段边界，直接形成对应的直线段线声源；对于弧形边界，形成可供现实计算的多条有限长直线段线声源来代替；
[0032] S120.获取接收点和声源之间的环境参数，所述环境参数包括：单位面声源强度、单位等效线声源强度、面声源的面积、面声源轮廓线的周长、线声源的两个端点、接收点到线源两端点的夹角、接收点到线声源的垂直距离、离散后线声源的个数和线声源的长度；
[0033] S130.根据线声源的长度、接收点到线源两端点的夹角等参数，基于声功率积分方法，计算每一条线声源在接收点处的声音辐射能量(声能贡献)；
[0034] 根据线声源对接收点的声功率贡献计算任意线声源对接收点的声音辐射能量(声能贡献)。
[0035] 根据线声源的两个端点区间、贡献任意单位线源对接收点的声功率、单位等效线声源强度、接收点到线源两端点的夹角和接收点到线声源的垂直距离等参数，应用声功率积分的方法，计算单个线声源对接收点的声功率贡献量。
[0036] 根据线声源的长度、单位等效线声源强度和接收点到线声源的垂直距离等参数，计算任意单位线源对接收点的声功率贡献。
[0037] S140.通过每一条线声源在接收点处的声音辐射能量叠加计算得到接收点处面声源对接收点的声音辐射能量；
[0038] 根据线声源对接收点的声音辐射能量、单位等效线声源强度、离散后线声源的个数、接收点到线源两端点的夹角和接收点到线声源的垂直距离等参数，计算面声源对接收点的声音辐射能量。
[0039] 实施例1
[0040] 本发明适用于不规则面声源声音几何发散的计算。面声源对接收点的声音辐射能量近似的看作面声源各个轮廓线对该接收点的声音辐射能量之和。
[0041] 首先，根据单位面声源强度、面声源面积、面声源轮廓线周长及面声源总功率参数，基于声源总功率相等原则，将面声源总功率分配到面声源的边界轮廓线上。
[0042] 面声源总功率等于各轮廓线声功率总和，具体如下式所示：
[0043]
[0044] 式中，StrengthA表示单位面声源强度，StrengthL表示单位等效线声源强度，AA表示面声源的面积，CL表示面声源轮廓线的周长。
[0045] 将面源总声功率分配到面源的边界轮廓线上，单位等效线声源强度StrengthL可按下式进行计算：
[0046]
[0047] 式中，StrengthA表示单位面声源强度，AA表示面声源的面积，CL表示面声源轮廓线的周长。
[0048] 然后，如附图2所示，根据面声源的几何形状，对面声源的边界轮廓线进行离散化，得到若干直线段轮廓线，并将所有直线段轮廓线均视为有限长线声源。面声源的边界轮廓，可以包括直线段和弧形，具体离散化过程为：对于直线段边界，直接形成对应的直线段线声源；对于弧形边界，形成可供现实计算的多条有限长直线段线声源来代替。
[0049] 然后，获取接收点和声源之间的环境参数，所述环境参数包括：单位面声源强度StrengthA、单位等效线声源强度StrengthL、面声源的面积AA、面声源轮廓线的周长CL、线声源的两个端点x1和x2、接收点到线源两端点的夹角θ、接收点到线声源的垂直距离d、离散后线声源的个数n和线声源的长度x。
[0050] 然后，如附图3所示，根据线声源的长度、接收点到线源两端点的夹角等参数，基于声功率积分方法，计算每一条线声源在接收点处的声能贡献。
[0051] 根据线声源对接收点的声功率贡献量，任意线声源对接收点的声能贡献值LL可按下式进行计算：
[0052] LL＝10lg(IL)
[0053] 式中，IL为此条线声源对接收点的声功率贡献量。
[0054] 根据线声源的两个端点区间、任意单位线源对接收点的声功率贡献、单位等效线声源强度、接收点到线源两端点的夹角和接收点到线声源的垂直距离等参数，单个线声源对接收点的声功率贡献量IL可按下式进行计算：
[0055]
[0056] 式中，x1、x2分别表示线声源的两个端点，Ix表示任意单位线源对接收点的声功率贡献，StrengthL表示单位等效线声源强度，θ为接收点到线源两端点的夹角，d为接收点到线声源的垂直距离。
[0057] 根据线声源的长度、单位等效线声源强度和接收点到线声源的垂直距离等参数，任意单位线源对接收点的声功率贡献Ix可按下式进行计算：
[0058]
[0059] 式中，StrengthL表示单位等效线声源强度，x表示线声源的长度，d为接收点到线声源的垂直距离。
[0060] 最后，通过每一条线声源在接收点处的声能贡献叠加计算得到接收点处面声源对接收点的声能贡献。
[0061] 根据线声源对接收点的声能贡献值、单位等效线声源强度、离散后线声源的个数、接收点到线源两端点的夹角和接收点到线声源的垂直距离等参数，可按下式对面声源对接收点的声能贡献LA进行计算：
[0062]
[0063] 式中，LLi为第i个线声源对接收点的声能贡献值，StrengthL表示单位等效线声源强度，n表示离散后线声源的个数，θi为接收点到第i个线源两端点的夹角，di为接收点到第i个线声源的垂直距离。
[0064] 以上所述的本发明的实施方式，并不构成对本发明保护范围的限定。任何在本发明的精神原则之内所作出的修改、等同替换和改进等，均应包含在本发明的权利要求保护范围之内。
