CN102663199A - 一种基于响应曲面分析的声屏障优化设计方法 - Google Patents

一种基于响应曲面分析的声屏障优化设计方法 Download PDF

Info

Publication number
CN102663199A
CN102663199A CN201210116552XA CN201210116552A CN102663199A CN 102663199 A CN102663199 A CN 102663199A CN 201210116552X A CN201210116552X A CN 201210116552XA CN 201210116552 A CN201210116552 A CN 201210116552A CN 102663199 A CN102663199 A CN 102663199A
Authority
CN
China
Prior art keywords
sound barrier
sound
design
barrier
analysis
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN201210116552XA
Other languages
English (en)
Other versions
CN102663199B (zh
Inventor
张�林
殷承启
许雪记
仝凯
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
China Design Group Co Ltd
Original Assignee
Jiangsu Provincial Communication Planning and Design Institute Co Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Jiangsu Provincial Communication Planning and Design Institute Co Ltd filed Critical Jiangsu Provincial Communication Planning and Design Institute Co Ltd
Priority to CN201210116552.XA priority Critical patent/CN102663199B/zh
Publication of CN102663199A publication Critical patent/CN102663199A/zh
Application granted granted Critical
Publication of CN102663199B publication Critical patent/CN102663199B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

本发明涉及一种基于响应曲面分析的声屏障优化设计方法，其特征在于，包括如下步骤：（1）计算无限长声屏障绕射声衰减量，得到有限长声屏障理论插入损失；（2）选取影响声屏障对各受声点插入损失的三个重要变量；（3）计算不同方案下声屏障对各受声点的插入损失和声屏障造价；（4）得到二次多元回归模型；（5）计算方差分析结果和显著性检验结果；（6）对各设计参数进行显著性影响分析；（7）对声屏障设计参数进行优化，从优化区域内分析得到最佳设计参数。本发明公开的基于响应曲面分析的声屏障优化设计方法实现对多个敏感目标进行降噪优化研究，有效控制声屏障的降噪效果和工程造价。

Description

一种基于响应曲面分析的声屏障优化设计方法
技术领域
本发明涉及一种声屏障优化设计的方法，更具体的说涉及利用Minitab软件进行声屏障优化设计的方法。
背景技术
近年来，随着我国环城高架、高速公路的快速发展，以及公众对声环境质量要求的日益提高，交通噪声的污染与控制受到越来越多的关注。声屏障作为最实用、有效的交通噪声防治措施之一，在交通环保领域得到广泛的应用。在同一段声屏障保护范围内，由于各处房屋高度、路肩高差、至路肩距离等众多因素不尽相同，声屏障会对各保护目标产生不同的降噪效果；而对于同一处受声点，声屏障设计高度、长度和起始位置也会对插入损失及工程造价带来显著差异。研究表明，部分声屏障设计长度、高度不足，或者安置位置不合理，保护区域内声环境质量不能完全达标。
国内对声屏障工程的优化设计研究已取得若干积极成果。结合工程实践，利用Cadna/A、SYSNOISE等声学分析软件模拟声波辐射、散射和传播过程，可以实现对声屏障降噪效果和等声级线分布预测。然而这些研究大多以声屏障设置路段内特定位置的声学响应为单一优化指标，对保护范围内多处敏感目标的降噪优化研究，以及声屏障工程造价的控制却少有关注，因而难以实现声屏障工程的最佳综合环境效益。
发明内容
为解决现有技术的不足，本发明的目的在于提供一种对多个敏感目标进行降噪优化研究，有效控制声屏障的降噪效果和工程造价的基于响应曲面分析的声屏障优化设计方法。
为达到上述目的，本发明是通过以下的技术方案来实现的：
一种基于响应曲面分析的声屏障优化设计方法，包括如下步骤：
（1）根据惠更斯-菲涅尔原理计算无限长声屏障绕射声衰减量，并进行物理量修正，得到有限长声屏障理论插入损失；
（2）基于理论插入损失计算方法，选取影响声屏障对各受声点插入损失的三个重要变量：高度、长度、起点桩号；
（3）根据Box-Behnken的中心组合设计原理与方法，利用Minitab软件自动生成声屏障设计方案，计算不同方案下声屏障对各受声点的插入损失和声屏障造价；
（4）利用Minitab设计软件对步骤（3）中计算的插入损失和声屏障造价数据整理分析拟合，得到二次多元回归模型：                                                
Figure 487905DEST_PATH_IMAGE001
，式中，Yi为因变量响应值，β0为常数项，βi为线性系数，βii为交互项系数，βij为二次项系数，Xi为自变量；
（5）对各回归方程系数及方程模型进行F检验和P检验，得出方差分析结果和显著性检验结果；
（6）根据模型生成的响应曲面及响应等高线图，对各设计参数进行显著性影响分析；
（7）根据步骤（5）中得到的结果和步骤（6）的曲面分析结果，利用Minitab软件中值线功能对声屏障设计参数进行优化，并使用Minitab优化器绘制优化图，从优化区域内分析得到最佳设计参数。
对本发明进一步限定的技术方案为：步骤（1）中所述的修正的物理量包括遮蔽角百分率和地面吸声衰减。
进一步的：步骤（3）中所述的计算声屏障造价的公式为：Y=（0.09+0.066X1）X2，其中，X1、X2分别为声屏障设计高度和长度。
本发明的有益效果是：本发明采用的响应曲面分析法是一种综合实验设计与数学建模的新型优化方法，在多因素优化分析研究中表现出显著的优越性。本发明通过在给定区域上构建因素与响应值间的明确函数关系，得到满足多个响应值的各因素的最佳组合，实现对多个敏感目标进行降噪优化研究，有效控制声屏障的降噪效果和工程造价。
附图说明
图1为本发明的基于响应曲面分析的声屏障优化设计方法的流程图；
图2为实施例1中声屏障高度、长度对A点插入损失的影响的响应曲面图；
图3为实施例1中声屏障高度、长度对A点插入损失的影响的等高线图；
图4为实施例1中声屏障高度、起点桩号对B点插入损失的影响的响应曲面图；
图5为实施例1中声屏障高度、起点桩号对B点插入损失的影响的等高线图；
图6为实施例1中声屏障长度、起点桩号对C点插入损失的影响的响应曲面图；
图7为实施例1中声屏障长度、起点桩号对C点插入损失的影响的等高线图；
图8为实施例1中A、B、C三点插入损失及声屏障造价等值线图；
图9为实施例1中A、B、C三点的重叠等值线和响应优化器求解图。
具体实施方式
以下结合附图和具体实施例对本发明作具体的介绍。
图1为本发明方法的流程图。
本实施例中，选取一段道路，在道路营运期间，部分房屋会受到一定程度的交通噪声影响，选取A（K3+880）、B(K3+795)和C(K3+760)处首排房屋为代表受声点，进行数据分析。
为避免构建数学模型过于繁冗难以分析，在建模过程中做出以下合理假设：
（1）交通噪声源为无限长线声源，取500Hz为声源等效频率；
（2） 声屏障型式为直立型，不考虑屏体吸声作用、透射声、反射声对声屏障插入损失的影响；
（3） 受声点至道路隔离栅之间为平坦疏松地面。
根据假设条件，一种基于响应曲面分析的声屏障优化设计方法，包括如下步骤：
（1）根据惠更斯-菲涅尔原理计算无限长声屏障绕射声衰减量，并进行物理量修正，修正的物理量包括遮蔽角百分率和地面吸声衰减，得到有限长声屏障理论插入损失。
（2）基于理论插入损失计算方法，选取影响声屏障对各受声点插入损失的三个重要变量：声屏障设计高度、长度、起点桩号，分别以X1 、X2 、X3代表，每个自变量的低、中、高设计水平分别以-1、0、+1编码。
响应曲面设计因素与设计水平的关系如表1所示：
表1：
Figure 883114DEST_PATH_IMAGE002
（3）根据Box-Behnken的中心组合设计原理与方法，利用Minitab软件自动生成声屏障设计方案，计算不同方案下声屏障对各受声点的插入损失和声屏障造价。计算声屏障造价的公式为：Y=（0.09+0.066X1）X2，其中，X1、X2分别为声屏障设计高度和长度。
响应曲面设计方案与结果如表2所示：
表2
Figure 572853DEST_PATH_IMAGE003
（4）利用Minitab设计软件对步骤（3）中计算的插入损失和声屏障造价数据整理分析拟合，得到二次多元回归模型：：
Figure 457632DEST_PATH_IMAGE004
，式中，Yi为因变量响应值，β0为常数项，βi为线性系数，βii为交互项系数，βij为二次项系数，Xi为自变量。
以声屏障设计高度、长度和起点桩号为自变量，A、B、C三处受声点插入损失为响应值的二次多元回归方程分别为：
Figure 36250DEST_PATH_IMAGE005
（5）对各回归方程系数及方程模型进行F检验和P检验，得出方差分析结果和显著性检验结果。
方差分析结构和显著性检验结果如表3所示：
表3：
Figure 235150DEST_PATH_IMAGE006
对表3的数据进行软件分析，结果显示：对各受声点声屏障插入损失建立的二次回归模型相关系数R2≥0.9736，校正相关系数R2adj≥0.9261，P值≤0.0009；拟合不足部分P值≥0.1061。这说明各回归模型拟合度较高，可以较好地反映声屏障设计参数对各受声点插入损失的影响规律。
（6）根据模型生成的响应曲面及响应等高线图，对各设计参数进行显著性影响分析，曲面斜率反应各因素的影响显著性水平，等高线弯曲程度反应因素间交互作用的强弱。
其中，图2为实施例1中声屏障高度、长度对A点插入损失的影响的响应曲面图。图3为实施例1中声屏障高度、长度对A点插入损失的影响的等高线图。根据图2和图3可以得出，A点插入损失受声屏障设计高度、长度和起点桩号等因素影响均较为显著。这是由于A点距路肩很近，声屏障高度的变化对绕射声程差影响较大；同时A点靠近声屏障设置终点，声屏障长度和起始桩号位置是决定遮蔽角百分率的主要因素。
图4为实施例1中声屏障高度、起点桩号对B点插入损失的影响的响应曲面图。图5为实施例1中声屏障高度、起点桩号对B点插入损失的影响的等高线图。根据图4和图5可以得出，B点位于声屏障保护范围的中部位置，遮蔽角百分率受声屏障设置路段范围的制约程度较小，故声屏障起始桩号位置对B点插入损失影响不显著。
图6为实施例1中声屏障长度、起点桩号对C点插入损失的影响的响应曲面图；图7为实施例1中声屏障长度、起点桩号对C点插入损失的影响的等高线图。根据图6和图7可以得出，C点距路肩较远，交通线声源有效长度较A点、B点增加，延长声屏障对增大遮蔽角百分率意义不大，故C点插入损失受声屏障设计长度影响不显著。
以上分析结果与方差分析中一次项系数P值检验结果是一致的。
（7）根据步骤（5）中得到的结果和步骤（6）的曲面分析结果，利用Minitab软件中值线功能对声屏障设计参数进行优化，并使用Minitab优化器绘制优化图，从优化区域内分析得到最佳设计参数。
声屏障工程设计实例中，其高度、长度和设置路段范围不仅取决于声敏感点的目标降噪值，同时受到工程造价和公路沿线设备设施的制约。为实现设计声屏障实用性和经济性的统一，利用Minitab软件中重叠等值线功能对声屏障设计参数进行优化。声屏障技术经济指标设定为：各处受声点夜间噪声达标；工程造价控制在理想的范围内。
使用Minitab优化器绘制优化图，如图8和图9所示，从优化区域内分析得到最佳设计参数：声屏障设计高度3.7m，设计长度214m，起点、终点桩号分别为K3+728和K3+942。在此设计条件下，预测声屏障工程对A、B、C处居民房屋的插入损失分别为：11.1dB、9.1dB和4.4dB，可以满足各受声点声环境质量标准；声屏障工程造价约为71.7万元。Minitab优化器绘制的重叠等值线和响应优化器求解图如图9所示。
上述实施例不以任何形式限制本发明，凡采用等同替换或等效变换的方式所获得的技术方案，均落在本发明的保护范围内。

Claims (3)

1.一种基于响应曲面分析的声屏障优化设计方法，其特征在于，包括如下步骤：
（1）根据惠更斯-菲涅尔原理计算无限长声屏障绕射声衰减量，并进行物理量修正，得到有限长声屏障理论插入损失；
（2）基于理论插入损失计算方法，选取影响声屏障对各受声点插入损失的三个重要变量：高度、长度、起点桩号；
（3）根据Box-Behnken的中心组合设计原理与方法，利用Minitab软件自动生成声屏障设计方案，计算不同方案下声屏障对各受声点的插入损失和声屏障造价；
（4）利用Minitab设计软件对步骤（3）中计算的插入损失和声屏障造价数据整理分析拟合，得到二次多元回归模型：                                                
Figure 201210116552X100001DEST_PATH_IMAGE001
，式中，Yi为因变量响应值，β0为常数项，βi为线性系数，βii为交互项系数，βij为二次项系数，Xi为自变量；
（5）对各回归方程系数及方程模型进行F检验和P检验，得出方差分析结果和显著性检验结果；
（6）根据模型生成的响应曲面及响应等高线图，对各设计参数进行显著性影响分析；
（7）根据步骤（5）中得到的结果和步骤（6）的曲面分析结果，利用Minitab软件中值线功能对声屏障设计参数进行优化，并使用Minitab优化器绘制优化图，从优化区域内分析得到最佳设计参数。
2.根据权利要求1所述的一种基于响应曲面分析的声屏障优化设计方法，其特征在于，步骤（1）中所述的修正的物理量包括遮蔽角百分率和地面吸声衰减。
3.根据权利要求1所述的一种基于响应曲面分析的声屏障优化设计方法，其特征在于，步骤（3）中所述的计算声屏障造价的公式为： Y=（0.09+0.066X1）X2，其中，X1、X2分别为声屏障设计高度和长度。
CN201210116552.XA 2012-04-20 2012-04-20 一种基于响应曲面分析的声屏障优化设计方法 Expired - Fee Related CN102663199B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201210116552.XA CN102663199B (zh) 2012-04-20 2012-04-20 一种基于响应曲面分析的声屏障优化设计方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201210116552.XA CN102663199B (zh) 2012-04-20 2012-04-20 一种基于响应曲面分析的声屏障优化设计方法

Publications (2)

Publication Number Publication Date
CN102663199A true CN102663199A (zh) 2012-09-12
CN102663199B CN102663199B (zh) 2015-05-13

Family

ID=46772690

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201210116552.XA Expired - Fee Related CN102663199B (zh) 2012-04-20 2012-04-20 一种基于响应曲面分析的声屏障优化设计方法

Country Status (1)

Country Link
CN (1) CN102663199B (zh)

Cited By (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103335709A (zh) * 2013-06-26 2013-10-02 中山大学 基于指向性声源的声屏障过渡区域绕射衰减室外测量方法
CN103970954A (zh) * 2014-05-13 2014-08-06 温州大学 一种基于鲁棒设计的高密度集成电路封装的优化方法
CN109271753A (zh) * 2018-11-29 2019-01-25 东北大学 基于响应面分析的沉钒废水氨氮回收的工艺参数优化方法
CN110357676A (zh) * 2019-05-30 2019-10-22 中山大学 基于响应曲面法确定鸟粪石回收氮磷工艺参数的方法
CN111612235A (zh) * 2020-05-13 2020-09-01 吉林省农业科学院 一种养殖污水处理磁分离工艺参数组合优化方法
CN111814358A (zh) * 2020-08-06 2020-10-23 中国电子科技集团公司第四十六研究所 一种微波复合介质基板综合性能的多因子共优化设计方法
CN112726858A (zh) * 2020-12-26 2021-04-30 广东电网有限责任公司电力科学研究院 一种基于声屏障的噪声治理优化方法
CN113652983A (zh) * 2021-08-20 2021-11-16 浙江绍兴零贝环保科技有限公司 一种基于机器学习的声屏障设计方法
CN118395843A (zh) * 2024-04-15 2024-07-26 中科海智(青岛)轨道交通研究院有限公司 一种道路声屏障声学设计模拟计算方法及系统
CN118395843B (zh) * 2024-04-15 2024-09-27 中科海智(青岛)轨道交通研究院有限公司 一种道路声屏障声学设计模拟计算方法及系统

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101515302A (zh) * 2008-05-16 2009-08-26 上海船舶运输科学研究所 道路声屏障声学设计方法
CN101582093A (zh) * 2008-05-16 2009-11-18 上海船舶运输科学研究所 一种操作简单、功能全面的声屏障结构设计方法

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101515302A (zh) * 2008-05-16 2009-08-26 上海船舶运输科学研究所 道路声屏障声学设计方法
CN101582093A (zh) * 2008-05-16 2009-11-18 上海船舶运输科学研究所 一种操作简单、功能全面的声屏障结构设计方法

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
卢向明: "道路声屏障声学特性与声学设计研究", 《浙江大学硕士学位论文》 *
沙学锋等: "基于SYSNOISE的环境", 《长春理工大学学报》 *
赵健植等: "基于响应曲面法的除雾器", 《中国电机工程学报》 *
魏显威等: "高速公路声屏障设计中几个问题", 《公路交通科技》 *

Cited By (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103335709A (zh) * 2013-06-26 2013-10-02 中山大学 基于指向性声源的声屏障过渡区域绕射衰减室外测量方法
CN103335709B (zh) * 2013-06-26 2014-12-17 中山大学 基于指向性声源的声屏障过渡区域绕射衰减室外测量方法
CN103970954A (zh) * 2014-05-13 2014-08-06 温州大学 一种基于鲁棒设计的高密度集成电路封装的优化方法
CN103970954B (zh) * 2014-05-13 2017-01-25 温州大学 一种基于鲁棒设计的高密度集成电路封装的优化方法
CN109271753A (zh) * 2018-11-29 2019-01-25 东北大学 基于响应面分析的沉钒废水氨氮回收的工艺参数优化方法
CN110357676A (zh) * 2019-05-30 2019-10-22 中山大学 基于响应曲面法确定鸟粪石回收氮磷工艺参数的方法
CN111612235A (zh) * 2020-05-13 2020-09-01 吉林省农业科学院 一种养殖污水处理磁分离工艺参数组合优化方法
CN111814358A (zh) * 2020-08-06 2020-10-23 中国电子科技集团公司第四十六研究所 一种微波复合介质基板综合性能的多因子共优化设计方法
CN112726858A (zh) * 2020-12-26 2021-04-30 广东电网有限责任公司电力科学研究院 一种基于声屏障的噪声治理优化方法
CN113652983A (zh) * 2021-08-20 2021-11-16 浙江绍兴零贝环保科技有限公司 一种基于机器学习的声屏障设计方法
CN113652983B (zh) * 2021-08-20 2022-11-29 浙江绍兴零贝环保科技有限公司 一种基于机器学习的声屏障设计方法
CN118395843A (zh) * 2024-04-15 2024-07-26 中科海智(青岛)轨道交通研究院有限公司 一种道路声屏障声学设计模拟计算方法及系统
CN118395843B (zh) * 2024-04-15 2024-09-27 中科海智(青岛)轨道交通研究院有限公司 一种道路声屏障声学设计模拟计算方法及系统

Also Published As

Publication number Publication date
CN102663199B (zh) 2015-05-13

Similar Documents

Publication Publication Date Title
CN102663199B (zh) 一种基于响应曲面分析的声屏障优化设计方法
Kragh et al. Development of two Danish building typologies for residential buildings
Gremmelspacher et al. Historical building renovation and PV optimisation towards NetZEB in Sweden
CN103020336B (zh) 一种等效led光源创建方法
CN108805321A (zh) 一种电动汽车充电站规划方法
CN103258235B (zh) 一种基于改进人工鱼群算法的供水管网改扩建优化方法
CN105426619A (zh) 一种确定盾构隧道纵向弯曲刚度的分析方法
CN113312686B (zh) 基于cfd环境数值模拟的立面垂直绿化排布方法及系统
CN104378769A (zh) 基于覆盖预测的td-scdma基站规划点自动选点方法
Baulac et al. Optimization of multiple edge barriers with genetic algorithms coupled with a Nelder–Mead local search
CN106707035B (zh) 无线电环境图场强参数估计算法
CN101873605B (zh) 一种网络规划中自适应传播环境分类方法
Zhang et al. Investigation of typical residential block typologies and their impact on pedestrian-level microclimate in summers in Nanjing, China
Van Renterghem et al. Landscaping for road traffic noise abatement: Model validation
CN105787642A (zh) 一种油气勘探中潜在商业区优选方法
CN101998421A (zh) 一种实现室内综合覆盖的室外天线放置方法
CN103774605B (zh) 一种提高环抱式港池水体交换能力的设计方法
Loikkanen et al. Multicriteria evaluation of sustainable energy solutions for Colosseum
Diz-Mellado et al. Unravelling the impact of courtyard geometry on cooling energy consumption in buildings
Wang et al. Urban network noise control based on road grade optimization considering comprehensive traffic environment benefit
Li et al. Optimizing of architectural geometry and tubular daylight guidance system based on genetic algorithm to enhance daylighting and energy performance in underground office buildings
Wang et al. Low-carbon-oriented commercial district urban form optimization and impact assessment analysis
CN103391557B (zh) 一种基于网格自动取定室内天线需求及布置的方法
CN105208566A (zh) 一种智能自组网无线通讯的物联网终端信息采集器
Yang et al. Morpho-based study on urban street spatial configuration: The case of Nanjing city, China

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant
C56 Change in the name or address of the patentee
CP03 Change of name, title or address

Address after: 210014 Ziyun Road, Qinhuai District, Jiangsu, No. 9, No.

Patentee after: CHINA DESIGN GROUP Co.,Ltd.

Address before: 210005 Zhongshan South Road, Jiangsu, No. 342, No.

Patentee before: Jiangsu Provincial Communications Planning and Design Institute Co.,Ltd.

CP01 Change in the name or title of a patent holder
CP01 Change in the name or title of a patent holder

Address after: 210014 No. 9, Ziyun Avenue, Qinhuai District, Jiangsu, Nanjing

Patentee after: China Design Group Co.,Ltd.

Address before: 210014 No. 9, Ziyun Avenue, Qinhuai District, Jiangsu, Nanjing

Patentee before: CHINA DESIGN GROUP Co.,Ltd.

CF01 Termination of patent right due to non-payment of annual fee
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20150513