基于极值选取疲劳应力筛选算法研究.docVIP

  • 6
  • 0
  • 约6.55千字
  • 约 12页
  • 2018-08-30 发布于福建
  • 举报

基于极值选取疲劳应力筛选算法研究.doc

基于极值选取疲劳应力筛选算法研究

基于极值选取疲劳应力筛选算法研究   摘 要:核电站关键设备在线疲劳监测系统需要将重要的中间输出结果,即应力数据回写到实时数据库PI系统中。对未经筛选的应力写入PI数据库非常耗时,并且很多应力对重现应力趋势并没有显著贡献。因此,对应力数据进行筛选尤为必要。提出了一种基于极值选取的应力筛选算法,在保证筛选前后应力趋势基本一致的前提下,极大缩短了应力回写时间,从而有效控制了应力计算步骤后处理的总体耗时。   关键词:疲劳应力;极值选取;筛选算法   中图分类号:TP312 文献标识码:A 文章编号:1672-7800(2016)003-0045-03   作者简介:甘国华(1982-),男,湖北当阳人,中核武汉核电运行技术股份有限公司工程师,研究方向为数据库管理、核电厂工程应用软件开发;陈银强(1981-),男,甘肃庄浪人,中核武汉核电运行技术股份有限公司高级工程师,研究方向为核电厂关键设备疲劳监测系统开发、核电厂疲劳分析评价;张娜妮(1988-),女,湖北荆门人,中核武汉核电运行技术股份有限公司助理工程师,研究方向为核电厂工程应用软件开发及测试。   0 引言   目前,国内许多电厂都在建立厂级监控信息系统(SIS),而实时数据库是电厂SIS体系结构的基础。国内电厂较常用的是美国OSI软件公司开发的PI(Plant Information)实时数据库系统[1]。PI实时数据库采用了“旋转门”压缩专利技术,使进入到PI数据库的数据经过高效的压缩算法处理,极大地节省了硬盘空间[2,3]。除此之外,PI数据库还提供了功能强大的二次开发接口,如PI-API、PI-SDK,使得电厂各专业领域的应用系统能够访问PI数据库实现特定的业务需求。   疲劳是造成核电站关键设备失效的主要原因之一,是长期承受扰动应力并在足够多的循环扰动作用下造成的[4]。国外许多核电站都通过疲劳监测系统(FMS)来监测关键设备的疲劳损伤状态,以应对核电站热疲劳相关事故[5]。田湾核电站一回路设备在线疲劳监测系统(JNPC_FMS)是集热工数据抽取、应力计算、疲劳寿命评价于一体的在线FMS软件,实时数据库采用PI,关系数据库采用Oracle,后台疲劳寿命计算程序采用Windows服务方式开发并连续运行, 集成至SIS平台中的前台B/S页面选用VS 2010开发。   为了间接验证疲劳损伤系数的合理性,开发团队需要将每月的应力计算结果保存至PI数据库中,以供客户借助SIS平台查看各监测点的月度应力趋势。经过简单的性能测试、分析,发现面对每月30万量级的应力数据(每条应力记录又包含6个应力分量),即便采用PI-API调用方式加上应力分量异步并行写入方式,耗时长达5分钟(CPU八核且单核频率2.4GHz),这已经接近月度应力计算耗时,意味着占用大量的硬件资源;另一方面,很多应力记录中存在大量并非特征值的非极值点,这类应力对重现应力趋势并没有显著贡献,且很多应力会被PI数据库的旋转门算法在归档过程中给过滤掉。因此,对应力数据进行筛选尤为必要。   本文先提出了应力筛选算法的基本原则、初步算法描述,接着用假设的测试应力数据验证该算法的正确性,然后用真实的应力数据验证正确性时发现了该算法的不足,并对该算法进行改进和再次验证,在保证筛选前后应力趋势基本一致的前提下,极大缩短了应力回写时间,从而有效控制了应力计算步骤后处理的总体耗时。   1 应力筛选算法初步设计及验证   1.1 问题提出及分析   由于非法的热工运行数据在执行应力计算之前,就通过热工数据筛选算法(包括为空、数值转换异常、超出热工参数值上下边界范围、超出温度变化速率等)被筛选掉了,因此,可以初步假定应力计算结果不会出现异常应力数据。   应力筛选处理的基本原则:一个月的应力计算结果,针对每个时刻的应力分量,除了首、尾应力分量保留外,其它时刻的应力分量只保留处于极值(峰值或谷值)的数据点,非极值点一律滤掉。   基于上述原则,通过理论分析可以得出:在机组满功率运行时,由于各工艺参数测量值几乎稳定不变,则得到的月度应力趋势应该是近乎于6条直线(每条直线代表一个应力分量的趋势),这种场景只需保留月初首个时刻的应力、月末最后一个时刻的应力;在机组降功率停堆(或升功率启堆)时,除了满功率运行的那一时段只需保留首尾时刻的两个应力外,存在应力波动的其它时刻除了保留较大的极值应力外,还将保留应力阶跃过程中出现的小幅波动的极值应力;在出现异常工况时,应力显著波动时的极值应力也应能够保留下来。   1.2 算法初步设计   为方便研究,需将应力分量描述成含有两个元素的结构体,形如struct SC {double T,double V}。其中,T表示监测时刻,V表示T时刻对应的应力分量(单位Pa)。在

文档评论(0)

1亿VIP精品文档

相关文档