1.阀门简介
某核二级电动闸阀（型号：NB200Z960Y-100R）安全级别为二级，抗震类型为1A类，其主要设计参数见表1。
表1 阀门固有频率
模态阶数
频率/Hz
模态振型说明
1
43.66
水平X方向（轴向）的振动
2
63.77
水平Z方向（侧向）的振动
3
262.21
垂直Y方向的振动
阀体、阀盖、闸板材料：SA351MCF8M，200℃下基本许用应力值为132 MPa；
支架材料：CF8制成，200℃下基本许用应力值为133 MPa。
阀体与管道通过焊接连接，阀体通过螺栓M33与阀盖连接，阀盖与支架通过螺栓M30连接，电动装置的重量约为180kg, 其重心坐标（-120,50,1021）。
核二级电动闸阀三维模型
2.振动模态计算
目的：确定第一阶固有频率是大于还是小于33Hz，从而确定是利用等效静力方法还是动力分析方法来进行后继计算。
对阀门进出口法兰的端面进行固定约束，进行模态计算，提取三阶固有频率数，其计
算结果见表1。
阀门前三阶模态振型
根据该阀门在中国核动力研究院进行的试验报告，其测得的前三阶固有频率分别为41.45Hz，56.38Hz，第三阶＞100Hz,计算基频与试验结果基本一致。由于该阀门最低固有频率都超过33Hz，可视为刚性结构，因此，可用等效静力分析方法进行抗震分析。
3.抗震分析模型
地震载荷输入
通常当设备固有频率大于33HZ时，可采用等效静力分析方法进行计算，即SSE地震载荷取三方向的加速度值做为地震输入。这里根据《研制任务书》规定,两个水平方向作用5g地震加速度值，垂直方向作用3.5g加速度值。
载荷组合与应力限值
抗震分析一般只对第二类和第四类工况进行，第二类工况考虑OBE地震，第四类工况考虑SSE地震。这里为了计算方便和结果保守起见，在进行建模计算和限值评定时，使用最严重的载荷组合（即事故工况D级载荷）和最严格的应力限值(即A级的应力限值)，见表2。因此，如果这种情况下的应力限值得到满足，则其它工况的计算结果也必小于相应的应力限值，即设计和A、B、C 级工况均可满足。
本次核阀抗震分析充分考虑设计压力、设计温度、设备自重、螺栓预紧力等使用载荷和地震载荷。阀门的设计压力为6.4Mpa，为计算保守，事故工况下作用压力取值为1.5倍设计压力，阀门的设计温度主要用来确定此温度下材料的许用应力，而不考虑由于温度作用产生的热应力（为二次应力，这里相对由压力引起的应力来说，数值很小，基本可忽略不计），螺栓预紧力值则根据满足紧密要求的设计计算获得，设备自重考虑为在设备施加1g的重力加速度。
表2 载荷组合及所应用的应力限值
工况
载荷组合
应力限值
设计和正常工况
（A级使用限制）
设计压力
设计温度
设备自重
其它运行载荷
σm≤1.0S
σm（或σL）+σb≤1.5S
事故工况
（D级使用限制）
事故工况下压力
事故工况下温度
设备自重
SSE载荷
其它运行载荷
σm≤2.0S
σm（或σL）+σb≤2.4S
注：表中的应力限值取自ASME规范第III 篇表NC-3521-1，表中标识符意义如下：
σm:总体薄膜应力强度 σb:一次弯曲应力强度 σL:局部薄膜应力 S:最大许用应力值
抗震分析有限元模型
在建模时,在不影响分析结果的前提下对阀门结构部件做一定的简化。同时考虑到电动装置的整体刚性很强, 它对阀门的主要作用是施加惯性载荷, 因此在计算中将电动装置简化为在其质心处的一个集中质量点, 通过多点约束的形式连接到支架的上端面上。
利用强大的前处理软件HyperMesh建立三维有限元模型，模型采用SOLID95实体单元和TARGE170、CONTA174接触单元及PRETS179预紧单元、MASS21质量单元, 一共划分129627个单元，373086个结点。阀门整体模型各部件之间的联接主要通过网格节点协调连续、MPC约束、建立接触对的方法进行处理。
考虑最恶劣的工况承受最不利的载荷组合，计算事故工况（1.5×设计内压+螺栓预紧力+设备自重+SSE载荷），在阀门进出口法兰的端面进行固定约束，作用均布压力为6.4*1.5=9.6MPa，螺栓预紧力值根据设计计算书获得，SSE地震载荷按照水平方向5g、垂直方向3.5g加速度设置，自重影响只在垂直方向设置为1g加速度。
阀门计算有限元模型
4.承压边界计算结果及评定
阀门承压边界的分析对象包括阀体、阀盖，利用ASME规范进行应力强度评定。而对于支架，由于采用不锈钢塑性材料，为了保守起见，采用第三强度理论进行评定。
阀体、阀盖应力评定
图4为阀门整机在内压、自重、螺栓预紧力和SSE地震载荷组合作用下的应力强度计算结果云图。选取的阀体、阀盖薄弱部位的应力评定截面（穿过壁厚的应力评定线）如图5所示，利用ANSYS后处理器的Path功能计算评定截面的薄膜应力σm和薄膜加弯曲应力σm +σb，最大薄膜应力发生在B-B截面，值为73.64MPa，最大薄膜加弯曲应力也产生在B-B截面，值为113.85MPa，本分析中用应力限值保守地取A级的限值，即σm≤1.0S,σm +σb≤1.5S，S为阀体、阀盖的许用应力，值为132MPa,薄弱部位的应力评定结果见表3。
图4 阀门总体应力强度云图 图5 应力评定截面示意图
表3 阀体、阀盖各薄弱部位的应力评定结果
应力评定线
分析内容
A-A
B-B
C-C
D-D
薄膜应力σm
43.10
73.64
38.85
57.07
薄膜应力限值1.0S
132
132
132
132
膜加弯曲应力σm +σb
97.20(I)
113.85(I)
78.80(I)
97.53(O)
膜加弯应力限值1.5S
198
198
198
198
结论
合格
合格
合格
合格
注：①应力单位均为MPa，I是指壁厚方向的内表面，O是指壁厚方向的外表面。
支架强度校核
在内压、自重、螺栓预紧力和SSE地震载荷组合作用下，支架的最大应力强度值为112.7 MPa，支架材料为SA-351 CF8M，根据ASME规范第二卷Part D Table 1A查得，其基本许用应力值为138MPa，满足强度要求。
5.螺栓与法兰计算及评定
中法兰联接螺栓
在内压、自重、螺栓预紧力和地震加速度联合作用下,法兰螺栓应力评定见表3，中法兰应力评定见表4。
中法兰联接螺栓的拉应力与剪应力
法兰纵向、径向、切向应力
6.结论
本抗震分析采用最恶劣的工况、最不利的载荷组合，计算阀门设备的各重要部件的应力，然后利用最严格的应力限值（评定准则）来进行评定。核二级电动闸阀经过抗震分析表明:
1）阀门自振频率高于33Hz，可近似认为是刚性结构，具有良好的抗震性能。
2）承压边界阀体、阀盖在各种规定载荷综合作用下产生的薄膜应力(σm)、薄膜加弯曲应力(σm+σb)都在适当的许用值范围内，支架的等效应力值小于许可应力，不会产生屈服破坏，满足强度要求。
3）中法兰及其联接螺栓的各类应力均在要求的限值内，能保持结构边界的连接。
因此，在内压、螺栓预紧力和地震载荷联合作用下，该阀门能保持结构的完整性。