1. 引言
核电厂设备运行的安全性极为重要。安全级设备或安全重要非安全级设备都有抗震要求，必须经过抗震鉴定 [1] [2] [3] 。抗震鉴定的目的是证明设备在承受安全停堆地震(SSE)作用力期间和(或)之后能够维持结构完整性和可运行性，并执行其安全功能的能力。设备的抗震鉴定方法通常有四种：试验法、分析法、分析与试验组合法、经验法 [4] [5] [6] [7] 。机械设备大多采用分析法鉴定，电气设备因不确定因素太多一般采用试验法鉴定。对于大型设备，因无法直接安装在试验台上进行试验鉴定，故一般采用分析鉴定。抗震分析主要验证结构完整性，对能动设备难以验证其功能完整性(即可运行性)。
风冷散热器是核电厂应急或备用柴油发电机组的重要辅助设备之一，其主要功能是为柴油发电机组提供循环冷却水，以保证机组的安全运行。本文采用分析与试验组合法对某核电厂风冷散热器进行抗震鉴定，先进行抗震分析验证结构完整性，再对设备的关键部件进行抗震试验，验证其能动部件(风机)的可运行性。通过对计算结果和试验结果的对比分析，验证了计算分析模型和试验样机设计的合理性。鉴定结果表明：风冷散热器的抗震性能满足规范要求。
Figure 1. Structure of the air cooler radiator
图1. 风冷散热器结构
2. 设备描述
某核电厂风冷散热器的设备等级为安全重要非安全级(即F-SC3级)，抗震类别为I类，必须经过抗震鉴定。该风冷散热器主要由3组换热模块及钢结构框架、端部封板、上平台及连接管道等部件组成，其外形尺寸为：长5188 mm × 宽2600 mm × 高2684 mm，总质量约11,750 kg，结构如图1所示。每一组换热模块主要由一台风机、一个高温水散热芯组、一个低温水散热芯组组成。钢结构框架由型钢、钢板焊接而成。换热模块通过连接螺栓固定在框架上部横梁上。框架立柱底部与基础预埋板焊接固定。
3. 抗震分析
3.1. 力学模型
采用Solidworks三维设计软件 [8] 建立风冷散热器的几何模型，将几何模型导入ANSYS Workbench有限元分析软件 [9] ，进行材料设置、网格划分、分析参数设置、载荷施加，最后进行求解及结果后处理。
采用SOLID186及SOLID187实体单元进行划分网格，有限元整体网格模型如图2所示。风机的电机采用质量块模拟。框架的8根立柱底部钢板与基础预埋板焊接固定，设置为固定边界条件，约束节点的6个方向的自由度。
Figure 2. Finite elements model of the air cooler radiator
图2. 风冷散热器有限元模型
3.2. 材料参数
风冷散热器各主要部件的材料及其力学性能参数见表1。
Table 1. Materials property parameters
表1. 材料的性能参数
3.3. 载荷条件
依据RCC-M规范 [10] ，考虑两种载荷工况分别对应A级和C级两种载荷准则：
1) 正常工况(A级准则)，载荷：DW + P；
2) 事故工况(C级准则)，载荷：DW + P + SSE。
其中，DW、P、SSE分别表示重力载荷、设计压力载荷、地震载荷。地震载荷采用风冷散热器安装位置的楼层反应谱的包络谱(见表2，阻尼比取4%)，三个坐标轴方向同时施加。
Table 2. Envelope of floor response spectrum
表2. 楼层反应谱的包络谱
3.4. 应力限值
根据RCC-M规范要求，采用线弹性分析方法，按最大主应力进行应力评定。A级准则和C级准则的应力限值见表3，其中，Pm为一次总体薄膜应力，Pb为一次弯曲应力，S为材料基本许用应力。对于实体单元，通过应力线性化方法获得薄膜应力和弯曲应力。
3.5. 计算结果
3.5.1. 固有频率
采用兰索士法(Lanczos method)求得结构100 Hz以内的所有频率。结构前4阶固有频率见表4，对应的模态振型见图3。
3.5.2. 结构应力
正常工况下的结构应力直接由静力分析计算得到，最大应力值79.554 MPa；地震响应采用反应谱分析计算得到；事故工况下结构应力由反应谱分析结果与静力分析结果叠加得到，最大应力值111.36 MPa。两种工况下结构最大主应力分布见图4，正常工况下的最大应力出现在框架中间一个立柱的顶部位置，事故工况下最大应力出现在框架四角一个立柱的顶部位置。
3.6. 应力评定
正常工况下结构最大主应力为79.554 MPa，小于薄膜应力限值94 MPa (A级准则)；事故工况下结构最大主应力为111.36 MPa，小于薄膜应力限值112 MPa (C级准则)。因此，无需应力线性化即可确定：结构的薄膜应力、薄膜 + 弯曲应力均小于各自的应力限值，强度满足规范要求。
4. 抗震试验
4.1. 试验样机设计
因地震模拟振动台尺寸限制，原风冷散热器无法直接安装在振动台上进行试验。取原设备的三分之一结构进行抗震试验，即取一组换热模块(包括一台风机)及缩短的钢结构框架(长度约为原尺寸的1/3，宽度、高度尺寸保持不变)组成试验样机，并配置了水箱、手动加压泵、阀门、压力表、管道、电源、电缆等辅助设施，满足样机在抗震试验时的运行工作条件。抗震试验布置见图5，图中的MA1、MD1、MS1等分别表示加速度、位移、应变传感器的测点编号。
4.2. 试验内容与要求
1) 抗震试验前的功能试验(即基准性能试验)，已在制造厂完成。现场进行密封性试验，检查确认各密封点没有泄漏。
2) 动态特性探查试验，测定结构的自振频率、振型。用加速度幅值为0.2 g的白噪声随机波在三个方向进行扫描，扫描时间不低于3分钟。
3) 抗震鉴定试验，验证设备的可运行性。按试验要求响应谱构造人造地震波(加速度时间历程)作为试验台的实际输入激励。人造地震波对应的试验反应谱应包络试验要求反应谱，且有10%的裕度。先进行5 次运行基准地震(OBE)鉴定试验，为考核设备在地震过程中执行安全功能的能力，其中二次OBE试验期间进行风机启动运行，其余三次保持运行；最后进行一次SSE 地震鉴定试验，期间也要求风机启动，验证电机和风扇的可运行性。
4) 抗震试验后的功能试验和最终检验，返厂进行并与基准性能试验进行对比。现场采用目测方法检查试验样机的结构完整性。
4.3. 试验结果
经过5次OBE地震和1次SSE地震考核试验，结果表明：风机正常启动并正常运行，结构件未见损伤及变形，螺栓、螺母没有松动、脱落，焊缝也无撕裂。因此，风冷散热器样机满足抗震性能要求。最后返厂进行功能试验，并与基准性能试验数据进行对比，结果表明：风冷散热器性能没有异常变化。
以下简单介绍抗震试验的主要测量结果。
1) 测点的自振频率
试验测得的结构整体振动的自振频率见表5。
Table 5. Natural frequencies by test
表5. 测点的自振频率
2) 测点的应变峰值
5次OBE试验和1次SSE试验得到的应变峰值见表6 (仅列出最大值)。两个测点位置都靠近框架立柱底部。
5. 计算与试验结果对比分析
5.1. 固有频率对比
比较表4和表5，样机的纵向和垂向的固有频率明显降低，横向固有频率变化很小，说明样机长度尺寸变短对结构纵向和垂向的固有特性影响较大，对结构横向固有特性影响很小，这与实际情况相符。比较结构横向固有频率，试验值为18.2 Hz，计算值为19.832 Hz，两者十分吻合。
5.2. 结构应力对比
根据表6测得的应变值，乘以材料弹性模量后可以得到结构的应力值。只考虑事故工况(SSE地震)条件，MS1测点(位于立柱根部)的应力为110.42 MPa，与计算模型的结果111.36 MPa (位于立柱顶部)已很接近。计算模型结果表明：立柱顶部的应力比底部应力大，由此推断，样机立柱顶部的应力应该比计算模型顶部应力还要大些，也就说明试验样机承受的地震载荷比原实际结构更大些。
5.3. 结构模型对比
通过上述分析可知，试验样机(一组换热模块)的固有频率比原结构(三组换热模块)低，结构应力比原结构大，说明试验样机在地震工况下的受载条件比原结构更为苛刻，因此，试验样机的抗震试验鉴定结果可以覆盖原结构的抗震鉴定要求，试验样机的设计是合理的。另一方面，计算模型的横向固有频率与试验样机的实测结果比较接近，验证了计算分析模型的准确性。
6. 结论
采用分析与试验组合法对某核电厂风冷散热器进行抗震鉴定。首先，采用有限元法进行抗震分析，依据RCC-M规范进行应力评定，验证了设备在地震时的结构完整性；其次，对风冷散热器(取一组换热模块作为试验样机)进行抗震试验，验证了设备在地震时的可运行性。最后，通过对试验结果和计算结果进行分析，验证了计算模型和试验样机设计的合理性。抗震鉴定结果表明，该风冷散热器结构设计满足规范要求，在地震载荷作用下能够保持结构完整性并执行其安全功能。本文采用的抗震鉴定方法对核电厂大型设备结构的抗震设计和抗震安全评估有一定的参考价值。