静电除尘器极板振动建模新进展：混合有限元与条带法的实测验证

基于波兰Bielsko-Biala工业大学Nowak团队对SIGMA收尘极板的实验对比研究

关键词
electrostatic precipitator, collecting electrode, hybrid finite element method, strip method, vibration analysis, 静电除尘器, 振打清灰

电除尘器（ESP）在燃煤电厂、水泥、冶金及垃圾焚烧等行业的超低排放改造中已经成为刚性“标配”。在长期运行中，影响静电除尘效率的关键因素并不只在电场区的荷电与迁移过程，机械系统的振打清灰效果同样至关重要。工业现场早有共识：如果收尘极板上的积灰不能被有效清除，再高端的高压电源或流场优化也难以挽回排放恶化的问题。因此，围绕“极板振动”这一传统环节进行精细化建模与实验验证，正逐渐成为静电除尘技术升级的一个隐性风口。

波兰Bielsko-Biala工业大学的Andrzej P. Nowak基于电除尘器收尘极板，完成了两种数值模型的实验验证工作，聚焦于典型的SIGMA波纹极板结构。文章标题为《Experimental Validation of Numerical Models of Collecting Electrodes》，通过将混合有限元方法（HFEM）与条带法（SPL）计算结果与实测振动加速度对比，为工程界提供了一套可推广的极板振动数值建模与验证框架，对提升ESP振打系统设计质量和运行诊断能力具有重要参考价值。

从工程机理看，电除尘器的振打系统通过锤击或电磁锤对刷击梁（brushing bar）施加轴向冲击，冲击力沿振打装置与收尘极板系统传播，在板面形成以切向与法向加速度为特征的复杂振动响应。工业界通常采用“法向最大加速度≥100 g”作为合格与否的经验判据，这是在早期测量手段有限背景下形成的规则[7]。随着三轴加速度传感器、数据记录系统和信号处理技术的发展，简单依赖单一峰值指标已经无法满足精细化判断和设计优化的需求，更系统的振动响应建模与实验验证变得尤为关键。

在建模方法上，Nowak团队提出了两种思路。第一种是混合有限元方法（Hybrid Finite Element Method, HFEM）。其核心思想是：将极板视为由多个刚体单元和弹性连接构成的柔性多体系统，其中弹性变形能采用传统有限元壳单元求解，而质量与惯性特性则用刚体有限元方法（Rigid Finite Element Method）集中表示[11]。在HFEM中，极板被划分为多个矩形壳单元，每个单元在局部坐标系下具有三平动、三转动共六个自由度，弹性应变能可整理成标准二次型，得到24×24的单元刚度矩阵。随后，通过刚体有限元对质量进行集中处理，使得系统总质量矩阵呈对角形式，极大简化了动力学方程的数值积分开销。这种“刚柔耦合”的混合有限元建模思路，尤其适合同时考虑极板、刷击梁与上部吊挂梁三部分联动的整体振动特性[6]。不过，由于对板面进行了相对细致的离散，HFEM模型自由度数量较大，计算时间在大尺度工业系统中仍然不容忽视。

第二种方法是条带法（Strip Method, SPL），可以看作一种半解析的有限条带方法[2]。该方法关注到实际收尘极板在长度方向（通常是垂直方向）几何形状相对恒定，而截面沿长度的变化有限，因此可将极板离散为若干条“纵向条带”，在条带长度方向采用B3样条函数近似挠度分布[14]。具体而言，对每一条带，法向挠度被表达为长度方向函数与宽度方向函数的乘积形式，长度方向函数由B3样条展开，时间相关系数则作为广义坐标。由于B3样条具备C2连续性、插值性好、能灵活施加不同边界条件，条带法不仅大幅降低了系统的自由度，而且刚度与质量矩阵天然呈稀疏结构，更适合进行重复算例和参数扫描。条带法尤其契合SIGMA型极板由多个折线长条组合而成的几何特征，能够在保证主导振型精度的前提下显著缩短计算时间[2]。

为了验证这两种收尘极板振动模型的可靠性，Nowak依托某静电除尘器制造企业搭建了全尺寸试验平台。试验对象为九块SIGMA型收尘极板，悬挂于公共吊梁底部，下端通过刷击梁和砧座联接。试验采用TEAC LX110 16通道记录仪和5只三轴ICP加速度传感器（PCB 356A02），在24 kHz采样率下记录单次锤击激励下的三向加速度响应。由于振打由人工锤击触发，数据采集无法完全自动同步，采集到的多次响应需要通过后处理进行对齐与时间同步，才能开展统计分析。每种传感器布置工况下重复试验10–25次，分别在竖向不同高度（I–V层）以及板面不同横向位置的控制点上获取足够多的样本，为后续的数值模型验证提供基础数据集。

在响应评价指标上，本文没有停留在传统的峰值比较，而引入了更具统计意义的runRMS（运行均方根）、RMS分布、FACk因子和命中率qε等量化指标。其中，runRMS是沿时间滑动窗口计算的瞬时均方根曲线，可用于刻画振动能量随时间的演变并提取最大runRMS值作为对比基准。在数字仿真方面，两种数值模型均采用实测冲击力F(t)作为输入边界，进而求解各控制点的总加速度时程，随后通过同样的runRMS与RMS处理流程与实验数据对齐比对。

FACk因子是一种用于模型评价的统计指标，最初应用于环境微尺度风场模型评估[3][10]。其思想是统计在所有测点中，有多少比例的模拟RMS结果落在实测RMS的某一倍数范围内（例如k=2，即0.5–2倍区间）。FACk越高，说明模型对整体空间分布的预测越可靠。qε命中率则从相对误差角度评价模型精度，即统计有多少测点的模拟RMS相对误差低于给定阈值ε（例如40%）。FACk偏重“合格率”，而qε更敏感于模型精细度，两者综合使用，可以更客观地给出数值模型的工程适用性判断。

在统一设置k=2、ε=0.4、全部35个测点参与统计的前提下，HFEM和SPL两种静电除尘极板模型在FACk和qε指标上呈现出有意思的差异。对于法向加速度，SPL条带模型的FACk约有95%的测点落入0.5–2倍区间，总加速度的FACk则约为90%；HFEM相应指标分别低约10%和5%。在对模型区分度更高的qε指标上，对法向加速度而言，SPL模型约有74%的测点在±40%误差范围内，总加速度约为84%；而HFEM仅有约42%和47%的测点满足该精度要求。这表明，在收尘极板这种长板型构件上，条带法凭借更紧凑的广义坐标表达和针对性几何抽象，在保证计算效率的同时，其对关键加速度响应的预测精度显著优于HFEM。

值得强调的是，作者在文中也明确指出，两种方法的精度均已达到可用于工程预评估的水准，其计算结果与商业有限元软件（如ABAQUS、ANSYS）的对比显示误差处于同一量级，而计算成本却明显更低[1][6]。这对于需要在电厂、冶金或垃圾焚烧项目中快速比选不同极板型式、振打参数或结构改造方案的工程师而言极具现实意义：

第一，这类“轻量级”数值模型可以在设计阶段帮助判断某一振打方案是否能在目标控制截面上实现“全覆盖”≥100 g法向加速度或更合理的三向加速度分布，减少对现场反复试打与调试的依赖。

第二，通过与实测加速度数据对标，当运行过程中监测到某些区域响应明显衰减或偏离模型预期时，可以及时怀疑挂耳松动、刷击梁变形、极板变形积灰异常等机械故障，从而将极板振动分析引入静电除尘器的状态诊断与预测性维护体系。

第三，在超低排放改造和高比集电面积（SCA）工况下，收尘极板高度与列数不断增加，振打系统能量分布的不均匀性问题会更加突出。SPL条带法这种兼具效率与精度的振动建模方法，非常适合作为电除尘器结构优化中“多工况批量计算”的内核，比如优化极板刚度分布、吊挂形式或振打锤质量与落距等关键参数[2][8]。

综上，Nowak团队的这项工作并非仅仅是对两种数值方法的学术比较，而是通过严谨的实验验证，建立起收尘极板实际振动行为与数值模型之间的可信链接。对于正在推进ESP–FF混合除尘、错流电除尘器、新型高频振打装置或智能振打控制的设备制造商与设计院来说，这一技术路线为“从经验振打走向可计算振打”提供了可落地的方法论。随着工业现场在线测振手段的进一步普及，把HFEM和SPL这种基于物理机理的模型与数据驱动的算法结合，有望形成面向下一代静电除尘器的数字孪生基础，为工业烟气治理提供更高效、更可控的振打清灰解决方案。

参考文献
[1] Adamiec-Wójcik I., Nowak A., Wojciech S. Comparison of methods for vibration analysis of electrostatic precipitators. Acta Mechanica Sinica. 2011;27(1):72-79.
[2] Adamiec-Wójcik I., Nowak A., Wojciech S. Dynamic analysis of electrostatic precipitators using finite strip method. In: DSTA 2009 Conference Proceedings. Vol. 2. Łódź, Poland; 2009. p. 889-896.
[3] Britter R., Schatzmann M., editors. Background and justification document to support the model evaluation guidance and protocol. COST Action 732. Brussels: COST Office; 2007.
[4] Long Z., Yao Q., Song Q., Li S. A second-order accurate finite volume method for the computation of electrical conditions inside a wire-plate electrostatic precipitator on unstructured meshes. Journal of Electrostatics. 2009;67(4):597-604.
[5] Neimarlija N., Demirdžić I., Muzaferija S. Finite volume method for calculation of electrostatic fields in electrostatic precipitators. Journal of Electrostatics. 2009;67(1):37-47.
[6] Nowak A.P., Adamiec-Wójcik I. Vibration analysis of collecting electrodes of precipitators by means of the hybrid finite element method. In: Multibody Dynamics, ECCOMAS Thematic Conference. Warsaw, Poland; 2009. p. 214-217.
[7] Nowak A. Measuring evaluation of vibration quality of the collecting electrodes section in ESP (in Polish). Measurement Automation and Monitoring. Warsaw; 2011. (in print).
[8] Nowak A.P., Wojciech S. Optimisation and experimental verification of a dust-removal beater for the electrodes of electrostatic precipitators. Computers & Structures. 2004;82(22):1785-1792.
[9] Talaie M.R. Mathematical modelling of wire-duct single-stage electrostatic precipitators. Journal of Hazardous Materials. 2005;124(1–3):44-52.
[10] VDI 3783 Part 9. Environmental meteorology – Prognostic microscale wind field models – Evaluation for flow around buildings and obstacles. Düsseldorf: VDI; 2005.
[11] Wittbrodt E., Adamiec-Wójcik I., Wojciech S. Dynamics of Flexible Multibody Systems: Rigid Finite Element Method. Berlin–Heidelberg–New York: Springer; 2006.
[12] Yang X.F., Kang Y.M., Zhong K. Effects of geometric parameters and electric indexes on the performance of laboratory-scale electrostatic precipitators. Journal of Hazardous Materials. 2009;169(1–3):941-947.
[13] Zhao X., Luo S. A discussion on the ESP-FF hybrid precipitator. In: Electrostatic Precipitation, 11th International Conference on Electrostatic Precipitation (ICESP XI). Hangzhou, China: Zhejiang University Press–Springer; 2008. p. 472-474.
[14] Zienkiewicz O.C., Taylor R.L. The Finite Element Method. Vol. 2: Solid Mechanics. Oxford: Butterworth-Heinemann; 2003.

获取更多静电除尘相关专业论文，请访问 https://isesp.org/conference-papers/