一种考虑流固耦合的核电厂外层安全壳外挂水箱结构的抗震设计方法技术

本发明专利技术涉及一种考虑流固耦合的核电厂外层安全壳外挂水箱结构的抗震设计方法，该方法根据结构和水的不同材料特性，采用不同的有限元网格描述方式，厂房墙体和水箱结构采用拉格朗日网格方式模拟，水体部分采用欧拉网格方式模拟，采用流固耦合分析方法和显式动力学分析方法，从不同材料网格特性、相互接触关系、流体状态方程与结构动力学方程等方面出发，真实模拟地震作用下外挂水箱中水与结构的相互作用。本发明专利技术可解决形状复杂、液面净空较小的水箱结构的水和结构在地震作用下相互耦合作用的分析问题。

【技术实现步骤摘要】
【专利摘要】本专利技术涉及，该方法根据结构和水的不同材料特性，采用不同的有限元网格描述方式，厂房墙体和水箱结构采用拉格朗日网格方式模拟，水体部分采用欧拉网格方式模拟，采用流固耦合分析方法和显式动力学分析方法，从不同材料网格特性、相互接触关系、流体状态方程与结构动力学方程等方面出发，真实模拟地震作用下外挂水箱中水与结构的相互作用。本专利技术可解决形状复杂、液面净空较小的水箱结构的水和结构在地震作用下相互耦合作用的分析问题。【专利说明】
本专利技术涉及一种核电厂外层安全壳外挂水箱结构的抗震设计方法，尤其涉及形状不规则、液面净空较小、地震作用下水与结构的相互作用较为明显的水箱结构。
技术介绍
核电厂对地震作用下的安全性有很高要求，尤其是核电厂中的水箱结构，在地震作用下，水箱中的水的晃动效应会对结构产生附加的内力响应，准确分析该内力响应，对核电厂水箱的设计具有重要的指导作用。目前国际上已有的核电厂水箱结构设计中，水在地震作用下的晃动效应，通常通过以下两种方法计算:(1)将水箱内的水转化为集中质量固结于水箱墙体上，直接进行整体抗震分析，而忽略了水箱水的晃动参与质量、以及水-结构流固耦合作用的影响，计算结果偏于保守。(2)水对结构的作用分为晃动水和冲动水两部分，根据经验公式计算得到冲动水的参与质量，模拟质量点参与结构整体抗震分析。而局部的晃动力作用，通过经验公式计算液面晃动高度的形式估算。此方法无法考虑晃动水对于水箱顶板的冲击作用，因此不适用于水箱形状复杂、水面净空较小的情况。
技术实现思路
本专利技术的目的是克服目前国际上已有方法比较简化、计算结果偏于保守、无法分析形状复杂和水面净空较小的水箱结构等问题，提供一种能够充分考虑水和结构的相互作用，并且可应用于形状复杂、水面净空较小的核电厂外层安全壳外挂水箱结构的抗震设计方法。本专利技术的技术方案如下:一种考虑流固稱合的核电厂外层安全壳外挂水箱结构的抗震设计方法，包括如下步骤:(I)根据核电厂外层安全壳及防护厂房的设计图纸，建立与实际结构一致的三维有限元分析模型，其中，厂房墙体和水箱结构采用拉格朗日网格方式模拟，水箱中的水采用欧拉网格方式模拟；(2)分别确定结构和水体的有限单元类型、材料特性及本构模型；(3)确定流固耦合分析方法及水体与结构的接触关系；(4)确定地震作用荷载及边界条件；(5)采用显示动力学分析方法进行求解，得到在地震作用下的结构响应；(6)将得到的结构在地震作用下的内力响应用于结构配筋设计。进一步，如上所述的考虑流固耦合的核电厂外层安全壳外挂水箱结构的抗震设计方法，步骤(I)中采用ABAQUS有限元分析软件建立三维有限元分析模型。进一步，如上所述的考虑流固耦合的核电厂外层安全壳外挂水箱结构的抗震设计方法，步骤(2)中，结构的有限单元类型为三节点壳单元和四节点壳单元，水体的有限单元类型为缩减积分欧拉单元；结构的混凝土材料采用线弹性材料模型，水体的材料特性为不可压缩粘性流体。更进一步，如上所述的考虑流固耦合的核电厂外层安全壳外挂水箱结构的抗震设计方法，步骤(2)中，水的材料本构模型采用状态方程的形式描述，状态方程根据单位质量的内能增量等于应力总功率与热功率的和的方式建立。进一步，如上所述的考虑流固耦合的核电厂外层安全壳外挂水箱结构的抗震设计方法，步骤(3)中，采用相容拉格朗日-欧拉法模拟流固耦合作用关系，水体与结构的接触关系通过罚函数法实现，根据流体与结构的接触力为罚函数力的反作用力，在计算接触力的过程中通过不断迭代，平衡流体与结构的相互之间的相对位移，求得每一时刻的流体与结构的接触力。进一步，如上所述的考虑流固耦合的核电厂外层安全壳外挂水箱结构的抗震设计方法，步骤(4)中，所述的边界条件为模型底部施加X，Y，Z三个方向加速度时程，分别代表地震作用在三个方向的分量，峰值加速度为0.3g。进一步，如上所述的考虑流固耦合的核电厂外层安全壳外挂水箱结构的抗震设计方法，步骤(5)中，采用显式动力学求解器ABAQUS/Explicit对核电厂外层安全壳及防护厂房结构的地震响应进行计算分析，获得结构和流体在地震加速度荷载作用下的结构位移时程、结构单元应力、结构单元内力、水的晃动效果时程、晃动高度、结构与流体的相互接触力。进一步，如上所述的考虑流固耦合的核电厂外层安全壳外挂水箱结构的抗震设计方法，步骤(6)中，结合Matlab程序化设计，首先，对结构单元的各个不同内力的时程提取各自的最大值；然后，进行单元内力最不利组合；最后，得到组合内力结果后，参考能源行业标准NB/T20012-2010《压水堆核电厂核安全有关的混凝土结构设计要求》对结构进行配筋计算。与现有的核电厂水箱结构抗震分析方法相比，本专利技术具有如下特点:1、采用国际上得到广泛认可的大型通用有限元分析软件ABAQUS，模拟了整个核电厂外层安全壳及防护厂房结构的三维模型，结构模型具有真实性、准确性、和可靠性。2、对于水箱中水的模拟，采用了更适用于描述变形较大、流动性较强的流体物质的欧拉网格描述方式，能更准确有效的描述在地震作用下的水箱中水的运动状态。3、本专利技术采用流固耦合分析方法和显式动力学分析方法，从不同材料网格特性、相互接触关系、流体状态方程与结构动力学方程等方面出发，真实模拟地震作用下外挂水箱中水与结构的相互作用。4、通过分析求解得到结构内力响应结果后，提取结构各个单元整个地震作用时程中的最大内力响应值，用于结构配筋设计，拓展了时程分析在结构设计中的应用。结合Matlab程序化设计，对结构内力输出结果进行处理、配筋计算等。5、本专利技术方法可以对外层安全壳外挂水箱相关部位结构设计进行优化。通过外挂水箱流固稱合分析找到结构薄弱环节，提出如改变结构断面等相应的优化措施，进一步提高结构的安全性和可靠性。6、本专利技术方法可解决形状复杂、液面净空较小的水箱结构的水和结构在地震作用下相互耦合作用的分析问题。【专利附图】【附图说明】图1为外层安全壳外挂水箱平面不意图；图2为外层安全壳外挂水箱剖面不意图；图3为本专利技术的方法原理示意图；图4为本专利技术中流体与结构接触关系的定义方式示意图。图1、图2中，10.外层安全壳11.水箱A12.水箱B13.水箱C图3中，1.水箱结构(拉格朗日网格)2.欧拉网格3.欧拉网格部分位置填入水材料4.结构与水之间设置接触关系5.结构底部输入地震加速度时程【具体实施方式】下面结合附图和实施例对本专利技术进行详细的描述。( I)建立三维有限元分析模型外挂环形水箱属于核电厂外层安全壳及防护厂房，位于外层安全壳筒体10上部，如图1和图2所示，水箱整体为圆环形，由隔板分三段，如图中的水箱A-11，、水箱B-12、水箱C-13，每段水箱中有部分位置水深较深，水深为11.68m，水底标高为42.3m，其他位置水深为3.38m，水底标高50.6m ;水面标高为53.98m，水面净空1.62m，水量共计约为3000吨。在本专利技术方法分析中，采用国际通用有限元分析软件ABAQUS建立三维有限元模型。根据结构和水的不同材料特性，采用不同的有限元网格描述方式。厂房墙体和水箱结构采用拉格朗日网格方式模拟(如图3`中标号1)，拉格朗日单元网格随材料的变形而变化，材料边界清晰，是描述固体结构的最常用方法。水箱中的水采用欧拉网格方式模拟本文档来自技高网...

【技术保护点】
一种考虑流固耦合的核电厂外层安全壳外挂水箱结构的抗震设计方法，包括如下步骤：（1）根据核电厂外层安全壳及防护厂房的设计图纸，建立与实际结构一致的三维有限元分析模型，其中，厂房墙体和水箱结构采用拉格朗日网格方式模拟，水箱中的水采用欧拉网格方式模拟；（2）分别确定结构和水体的有限单元类型、材料特性及本构模型；（3）确定流固耦合分析方法及水体与结构的接触关系；（4）确定地震作用荷载及边界条件；（5）采用显示动力学分析方法进行求解，得到在地震作用下的结构响应；（6）将得到的结构在地震作用下的内力响应用于结构配筋设计。

【技术特征摘要】

【专利技术属性】
技术研发人员：蔡利建，姚迪，王黎丽，李玉民，张超琦，
申请(专利权)人：中国核电工程有限公司，
类型：发明
国别省市：北京;11