考虑液体晃动的反应堆厂房结构确定方法及装置与流程

1.本发明涉及核电厂技术领域，尤其涉及一种考虑液体晃动的反应堆厂房结构确定方法及装置。


背景技术：

2.随着全球能源需求的不断增长，化石燃料等不可再生能源的日益枯竭，核电作为一种经济、高效的清洁能源，在世界范围内得到了广泛应用。然而一旦发生核事故，会对人们的生命财产安全造成毁灭性的灾难。福岛核电站受地震引发的核泄漏事故重新引起了世界范围内对核电站在地震灾害下安全性的重视，影响核电站抗震性能的重要因素之一是考虑液体晃动的反应堆厂房结构的设计。
3.现有技术中采用欧拉网格模拟液体，无法真实模拟反应堆厂房中敞口容器内液体运动，计算结果不够精确，导致反应堆厂房结构设计不够合理。


技术实现要素：

4.本发明实施例提供一种考虑液体晃动的反应堆厂房结构确定方法及装置，解决了耦合欧拉-拉格朗日法计算结果不够精确，导致反应堆厂房结构设计不够合理的问题。
5.为达到上述目的，本发明实施例第一方面提供一种考虑液体晃动的反应堆厂房结构确定方法，包括如下步骤：
6.获取核电厂布置信息和荷载信息；
7.根据核电厂布置信息和荷载信息，建立反应堆厂房结构的几何模型，所述几何模型包括固体部件和液体部件；
8.基于所述几何模型，确定所述固体部件和所述液体部件的接触属性，且将所述液体部件转换为粒子，以得到反应堆厂房仿真模型；
9.根据所述几何模型、所述反应堆厂房仿真模型、核电厂布置信息和荷载信息，确定反应堆厂房的结构。
10.可选的，所述基于所述几何模型，确定所述固体部件和所述液体部件的接触属性，且将所述液体部件转换为粒子，以得到反应堆厂房仿真模型，还包括：
11.对所述几何模型中的所述固体部件和所述液体部件进行网格划分，获得液体部件网格；
12.基于所述液体部件网格，按等参方向和预设粒子数目将所述液体部件转换为粒子，得到反应堆厂房仿真模型。
13.可选的，所述根据所述几何模型、所述反应堆厂房仿真模型、核电厂布置信息和荷载信息，确定反应堆厂房的结构，包括：
14.基于所述几何模型，选取预设位置；
15.根据核电厂布置信息、荷载信息、所述几何模型和所述反应堆厂房仿真模型，获得所述预设位置的地震剪力时程曲线、所述固体部件的位移场、应力场和内力；
16.根据所述预设位置的地震剪力时程曲线、所述固体部件的位移云图、应力云图和内力，确定反应堆厂房结构设计结果。
17.可选的，所述基于所述液体部件网格，按等参方向和预设粒子数目将所述液体部件转换为粒子，得到反应堆厂房仿真模型，包括：
18.基于所述液体部件网格，按等参方向和预设粒子数目将所述液体部件转换为粒子，得到初始模型；
19.基于初始模型，进行粒子数量敏感性分析，获得分析结果；
20.根据所述分析结果，调整所述初始模型，获得中间模型；
21.将所述中间模型作为初始模型，执行利用初始模型，进行粒子数量敏感性分析，获得分析结果以及根据所述分析结果，调整所述初始模型，获得中间模型的步骤，直到所述分析结果符合预设范围；
22.将最后获得的中间模型确定为所述反应堆厂房仿真模型。
23.本发明实施例第二方面提供一种考虑液体晃动的反应堆厂房结构确定装置，包括：
24.获取模块，用于获取核电厂布置信息和荷载信息；
25.第一确定模块，用于根据核电厂布置信息和荷载信息，建立反应堆厂房结构的几何模型，所述几何模型包括固体部件和液体部件；
26.第二确定模块，基于所述几何模型，确定所述固体部件和所述液体部件的接触属性，且将所述液体部件转换为粒子，以得到反应堆厂房仿真模型；
27.第三确定模块，用于根据所述几何模型、所述反应堆厂房仿真模型、核电厂布置信息和荷载信息，确定反应堆厂房的结构。
28.可选的，所述第二确定模块，包括：
29.第一确定子模块，用于对所述几何模型中的所述固体部件和所述液体部件进行网格划分，获得液体部件网格；
30.第二确定子模块，用于基于所述液体部件网格，按等参方向和预设粒子数目将所述液体部件转换为粒子，得到反应堆厂房仿真模型。
31.可选的，所述第三确定模块，包括：
32.基于所述几何模型，选取预设位置；
33.根据核电厂布置信息、荷载信息、所述几何模型和所述反应堆厂房仿真模型，获得所述预设位置的地震剪力时程曲线、所述固体部件的位移场、应力场和内力；
34.根据所述预设位置的地震剪力时程曲线、所述固体部件的位移云图、应力云图和内力，确定反应堆厂房结构设计结果。
35.可选的，所述第二确定子模块，包括：
36.基于所述液体部件网格，按等参方向和预设粒子数目将所述液体部件转换为粒子，得到初始模型；
37.基于初始模型，进行粒子数量敏感性分析，获得分析结果；
38.根据所述分析结果，调整所述初始模型，获得中间模型；
39.将所述中间模型作为初始模型，执行利用初始模型，进行粒子数量敏感性分析，获得分析结果以及根据所述分析结果，调整所述初始模型，获得中间模型的步骤，直到所述分
析结果符合预设范围；
40.将最后获得的中间模型确定为所述反应堆厂房仿真模型。
41.本发明实施例第三方面提供一种电子设备，包括处理器、存储器及存储在所述存储器上并可在所述处理器上运行的计算机程序，所述计算机程序被所述处理器执行时实现如第一方面所述的考虑液体晃动的反应堆厂房结构确定方法中的步骤。
42.本发明实施例第四方面提供一种可读存储介质，所述可读存储介质上存储有程序，所述程序被处理器执行时实现如第一方面所述的考虑液体晃动的反应堆厂房结构确定方法中的步骤。
43.本发明实施例中，将液体部件转换为粒子，得到反应堆厂房仿真模型，可以模拟出敞口水箱内液体晃动会溅溢出水箱外的情形。反应堆厂房中有敞口水箱，对比现有技术中的耦合欧拉-拉格朗日法因为需定义液体的计算域，液体只能在计算域范围内流动，通常耦合欧拉-拉格朗日法定义液体的计算域为水箱的容积范围，即指定了液体只能在水箱内流动，无法模拟出敞口水箱液体晃动时可能会溅溢出水箱外的情形，解决了耦合欧拉-拉格朗日法计算结果不够精确，导致反应堆厂房结构设计不够合理的问题。
附图说明
44.为了更清楚的说明本发明实施例中的技术方案，现对说明书附图作如下说明，显而易见地，下述附图仅是本发明的实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据所列附图获得其他附图。
45.图1是本发明实施例提供的考虑液体晃动的反应堆厂房结构确定方法的流程图之一；
46.图2是本发明实施例提供的考虑液体晃动的反应堆厂房结构确定装置的结构图；
47.图3是本发明实施例提供的电子设备的结构图；
48.图4是本发明实施例提供的考虑液体晃动的反应堆厂房结构确定方法的流程图之二。
具体实施方式
49.下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅是本发明一部分实施例，而不是全部的实施例。在本发明中的实施例的基础上，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
50.参见图1和图4，本发明实施例提供一种考虑液体晃动的反应堆厂房结构确定方法，该方法可应用于电子设备，包括：
51.步骤101，获取核电厂布置信息和荷载信息；
52.步骤102，根据核电厂布置信息和荷载信息，建立反应堆厂房结构的几何模型，所述几何模型包括固体部件和液体部件；
53.步骤103，基于所述几何模型，确定所述固体部件和所述液体部件的接触属性，且将所述液体部件转换为粒子，以得到反应堆厂房仿真模型；
54.步骤104，根据所述几何模型、所述反应堆厂房仿真模型、核电厂布置信息和地震
荷载信息，确定反应堆厂房的结构。
55.在步骤101中，核电厂布置信息包括反应堆厂房的位置信息以及对反应堆厂房的使用要求。荷载，指的是使结构或构件产生内力和变形的外力及其它因素。对于反应堆厂房来讲，荷载包括反应堆厂房自重和地震作用下的荷载。
56.在步骤102中，可以首先根据核电厂布置信息和荷载信息，进行结构模型简化，再根据结构模型简化的结果建立几何模型。结构模型简化包括将反应堆厂房墙板简化为壳，反应堆厂房内部容器中的液体简化为实体。根据简化后的结果建立几何模型。建立几何模型时，应理解，反应堆厂房结构包括固体部件和液体部件，几何模型建模可以首先输入核电厂布置信息，然后根据核电厂布置信息和荷载信息在cad(计算机辅助设计)软件中建立几何模型。建立几何模型过程中，需要对几何模型定义固体部件、液体部件的材料属性和截面属性，具体的，定义固体部件的材料属性为各向同性材料，定义液体部件的材料属性为不可压缩理想液体。定义固体部件的截面属性为壳截面，定义液体部件的截面属性为实体截面。在整体坐标系中，将定义好材料属性和截面属性的固体部件和液体部件组装在一起，得到几何模型。
57.具体的，引入人工压缩性，将液体近似为弱可压缩液体，由以下可压缩液体状态方程，通过粒子的密度和内能计算粒子的压力。
[0058][0059][0060]
式中，ρ0为参考密度，ρ为密度，c0为对应于ρ0的人工声速，g为重力加速度。p为各向同性压力。
[0061]
在步骤103中，基于所述几何模型，确定所述固体部件和所述液体部件的接触属性，且将所述液体部件转换为粒子，以得到反应堆厂房仿真模型。
[0062]
应理解，接触属性包括固体部件、液体部件的接触表面之间的切向作用和法向作用。在本发明实施例中，不考虑摩擦导致的切向作用，只考虑传递正向压力的法向作用。
[0063]
将所述液体部件转换为粒子可以先将液体部件划分为网格单元，对液体部件进行拉格朗日六面体网格划分，基于划分后的拉格朗日六面体网格单元将液体部件转换为粒子。并对固体部件进行拉格朗日四边形网格划分，得到反应堆厂房仿真模型。
[0064]
具体的，对固体部件和液体部件进行拉格朗日网格划分，并在计算初始时刻(t＝0)对液体单元进行粒子离散，在固体与液体的交界面处的粒子与拉格朗日单元在保持接触的同时又防止其相互侵入。算法如下：
[0065]
初始时刻(t＝0)读入所定义的固体与液体的交界面处的拉格朗日网格单元节点的位置坐标和速度然后根据每个网格单元的节点坐标与速度平均值将其等效为一个空间坐标与速度为的粒子，即
[0066][0067]
式中，n为每个拉格朗日网格单元节点总数，β代表坐标方向，β均取区间[1,3]中的正整数，同时遵循张量运算的字母指标法。
[0068]
确定了粒子的空间坐标与速度之后，根据初始时刻该部分液体的材料属性赋予粒子密度、压力、内能等计算参数。耦合界面处生成的粒子按以下方程近似计算：
[0069][0070]
式中，ρ,v(v
α
,v
β
),e,x(x
α
,x
β
),σ
αβ
,t,f
α
,p,δ
αβ
,τ
αβ
,μ,ε
αβ
分别代表密度、速度、内能、空间坐标、液体粒子总应力张量、时间、作用在单位质量液体上的体力、各向同性压力、一个二元函数、黏性剪应力、黏性系数、剪应变率。
[0071]
为耦合界面上粒子的速度；为为耦合界面上粒子的压力。
[0072]
在固液耦合计算过程中，每个时间步开始时刻均更新耦合交界面处的位置与速度。同时为了正确实现耦合界面处力的传递，进行液体与固体之间的接触设置：耦合界面处的拉格朗日网格单元上的单元压力由各单元中心附近的粒子压力平均值求得。
[0073]
引入边界惩罚力以防止液体与固体交界面的穿透，对于流固耦合交界面的液体粒子，其受到的惩罚力f方向指向壁面的法向。
[0074][0075][0076][0077]
式中，n是液体与固体交界面处固体壁面的法线方向；ci为粒子i的声速；ψ为液体粒子距离固体壁面的距离；ξ为位于两固体壁边界粒子之间的液体粒子的位置ξi在两固体
壁边界粒子间连线上的投影；δb为任意两个相邻的边界粒子的间距；u
⊥
为液体粒子速度在固体壁面法线方向上的投影；z为测点液体粒子距离自由液面的深度。
[0078]
在步骤104中，根据所述几何模型、所述反应堆厂房仿真模型、核电厂布置信息和地震荷载信息，确定反应堆厂房的结构。
[0079]
首先根据所述根据几何模型、所述反应堆厂房仿真模型、核电厂布置信息和地震荷载信息计算求得反应堆厂房在所受的所有荷载作用下的位移场，再由反应堆厂房在所受的所有荷载作用下的位移场计算求得应力场，再根据求得的应力场求得内力和配筋面积，最后根据核电厂设计规范，参考内力和配筋面积，进行结构配筋设计，确定反应堆厂房的结构。因由位移场最终确定反应堆厂房的结构的过程为现有结构设计中常用求解手段，故在此不对其过程进行详细阐述。
[0080]
其中，求得反应堆厂房在所受的所有荷载作用下的位移场，需判断求解结果是否收敛，若求解结果收敛，再根据个人工程经验判断收敛的分析结果是否合理，如果不合理，则调整反应堆厂房仿真模型，直到收敛的求解结果正常。
[0081]
其中，求反应堆厂房在所受的所有荷载作用下的位移场时，需要施加荷载作用，施加荷载作用可以首先选择显示动力计算，分别建立显示动力计算的自重载荷步和地震作用载荷步，新建自重载荷步，x向地震作用载荷步或y向地震作用载荷步或z向地震作用载荷步。然后施加边界条件、反应堆厂房自重和地震作用下的荷载。
[0082]
具体的，计算过程中，被转换为粒子的液体部件与被划分为拉格朗日单元的固体部件之间的相互作用满足以下拉格朗日形式的偏微分方程，它是基于质量守恒、动量守恒、能量守恒这三条基本物理守恒定律，同时在运动的拉格朗日框架中引入对总时间的导数而得到的。
[0083][0084]
σ
αβ
＝-pδ
αβ
τ
αβ
[0085][0086][0087]
τ
αβ
＝με
αβ
[0088]
式中，α、β代表坐标方向，α和β均取区间[1,3]中的正整数，同时遵循张量运算的字母指标法；
[0089]
ρ,v(v
α
,v
β
),e,x(x
α
,x
β
),σ
αβ
,t,f
α
,p,δ
αβ
,τ
αβ
,μ,ε
αβ
分别代表密度、速度、内能、空间坐标、液体粒子总应力张量、时间、作用在单位质量液体上的体力、各向同性压力、一个二元函数、黏性剪应力、黏性系数、剪应变率。
[0090]
为方便计算，采用以下粒子形式近似上述基本控制方程：
[0091][0092][0093]
式中，i和j代表粒子编号，∏
ij
是人工黏性；w
ij
是光滑函数；h是光滑长度，初始值取为初始粒子间距的1.2倍；r是与粒子空间坐标和光滑长度有关的一个系数。
[0094]
本发明实施例中，将液体部件转换为粒子，得到反应堆厂房仿真模型，可以模拟出敞口水箱内液体晃动会溅溢出水箱外的情形。反应堆厂房中有敞口水箱，对比现有技术中的耦合欧拉-拉格朗日法因为需定义液体的计算域，液体只能在计算域范围内流动，通常耦合欧拉-拉格朗日法定义液体的计算域为水箱的容积范围，即指定了液体只能在水箱内流动，无法模拟出敞口水箱液体晃动时可能会溅溢出水箱外的情形，解决了耦合欧拉-拉格朗日法计算结果不够精确，导致反应堆厂房结构设计不够合理的问题。
[0095]
可选的，所述基于所述几何模型，确定所述固体部件和所述液体部件的接触属性，且将所述液体部件转换为粒子，以得到反应堆厂房仿真模型，还包括：
[0096]
对所述几何模型中的所述固体部件和所述液体部件进行网格划分，获得液体部件网格；
[0097]
基于所述液体部件网格，按等参方向和预设粒子数目将所述液体部件转换为粒子，得到反应堆厂房仿真模型。
[0098]
具体的，对固体部件进行拉格朗日四边形网格划分，对液体部件进行拉格朗日六面体网格划分。指定在求解初始t＝0时，液体的拉格朗日六面体网格按等参方向n个粒子的原则转换，n可以取区间[2,7]中的正整数，比如取n＝2。
[0099]
对液体部件进行拉格朗日六面体网格划分后再将液体的拉格朗日六面体网格按等参方向n个粒子的原则转换，可以提高计算结果的精确性。
[0100]
可选的，所述根据所述几何模型、所述反应堆厂房仿真模型、核电厂布置信息和荷载信息，确定反应堆厂房的结构，包括：
[0101]
基于所述几何模型，选取预设位置；
[0102]
根据核电厂布置信息、荷载信息、所述几何模型和所述反应堆厂房仿真模型，获得所述预设位置的地震剪力时程曲线、所述固体部件的位移场、应力场和内力；
[0103]
根据所述预设位置的地震剪力时程曲线、所述固体部件的位移云图、应力云图和内力，确定反应堆厂房结构设计结果。
[0104]
选取反应堆厂房固体部件中的筏板底面形心处为预设位置，在预设位置处施加边
界条件和地震作用求解筏板的位移场、应力场，得到单元内力，进行筏板配筋。应理解，筏板为基础工程中的一块混凝土板，板下是地基，板上面有柱，墙等。在预设位置与筏板底板之间建立平动和转动方向的运动耦合，并在预设位置处施加边界条件地震作用下的荷载。具体的，在自重载荷步，预设位置处的边界条件可视实际情况选择简化为固接荷载，或者简化为考虑土结相互作用的土弹簧淤泥；对固体部件中的筏板、安全壳和内部结构施加自重荷载。在x向或y向或z向地震作用载荷步，释放自重载荷步中x向或y向或z向的平动，并在预设位置处输入x向或y向或z向的地震时程曲线。经计算获得所述预设位置的地震剪力时程曲线、所述固体部件的位移场、应力场和内力。
[0105]
应理解，一个面通过运动耦合，耦合到一个控制点即预设位置上，那么这个面跟控制点之间就是一种刚性接触，面和控制点成为一个刚体，即相对位置不变，计算运动也是作为一个整体运动的。一般情况下，面上的应变和应力都是0。只有筏板底板需要施加边界条件和地震作用下的荷载。固体部件中的筏板、安全壳和内部结构都只需施加自重荷载。然后再求解位移场、应力场。因为假定固体部件中的筏板、安全壳和内部结构是一体的，共同作用的，虽然地震作用下的荷载只是加到筏板底板，但地震作用会通过筏板传给安全壳和内部结构的。所以，安全壳、内部结构不需要加边界和地震作用，地震作用会自动传递。
[0106]
求解预设位置处的地震剪力时程曲线，可以确定地震作用对结构产生的最大剪力，然后将这个最大剪力按结构构件刚度进行分配，得到具体某一结构构件的具体剪力值，然后进行抗剪配筋。求解预设位置处的地震剪力时程曲线可以提高计算结果的精确性。此外，设置预设位置，还便于施加边界条件和地震作用，否则需选择整个筏板底部进行加边界条件和地震作用。若选择整个筏板底部施加边界条件和地震作用，计算中的换算比较复杂，容易出错，效率更低。设置预设位置求解，计算效率更高。
[0107]
可选的，所述基于所述液体部件网格，按等参方向和预设粒子数目将所述液体部件转换为粒子，得到反应堆厂房仿真模型，包括：
[0108]
基于所述液体部件网格，按等参方向和预设粒子数目将所述液体部件转换为粒子，得到初始模型；
[0109]
基于初始模型，进行粒子数量敏感性分析，获得分析结果；
[0110]
根据所述分析结果，调整所述初始模型，获得中间模型；
[0111]
将所述中间模型作为初始模型，执行利用初始模型，进行粒子数量敏感性分析，获得分析结果以及根据所述分析结果，调整所述初始模型，获得中间模型的步骤，直到所述分析结果符合预设范围；
[0112]
将最后获得的中间模型确定为所述反应堆厂房仿真模型。
[0113]
具体的，液体的六面体网格等参方向粒子数目n取值2，施加荷载求取固体部件的第一位移场结果，液体的六面体网格等参方向粒子数目n取值3得到初始模型，施加荷载求取固体部件的第二位移场结果，分析第一位移场结果与第二位移场结果的误差，比如第一位移场结果减去第二位移场结果的差除以第二位移场结果，根据误差分析结果，调整初始模型中粒子的数目，获得中间模型，将所述中间模型作为初始模型，执行利用初始模型，进行粒子数量敏感性分析，获得分析结果以及根据所述分析结果，调整所述初始模型，获得中间模型的步骤，直到第(n 1)次与第n次的计算结果误差小于预设值，该预设值根据实际需要设置，本发明实施例中可以取该预设值为5％，可取液体的六面体网格等参方向粒子数目
为n的中间模型确定为反应堆厂房仿真模型。
[0114]
经过如上所述的粒子数量敏感性分析，可以确定更准确选取液体粒子数，以确保计算精度和提高计算效率。
[0115]
参见图2，图2是本发明实施例提供的考虑液体晃动的反应堆厂房结构确定装置的结构图，如图2所示，考虑液体晃动的反应堆厂房结构确定装置200包括：
[0116]
获取模块201，用于获取核电厂布置信息和荷载信息；
[0117]
第一确定模块202，用于根据核电厂布置信息和荷载信息，建立反应堆厂房结构的几何模型，所述几何模型包括固体部件和液体部件；
[0118]
第二确定模块203，基于所述几何模型，确定所述固体部件和所述液体部件的接触属性，且将所述液体部件转换为粒子，以得到反应堆厂房仿真模型；
[0119]
第三确定模块204，用于根据所述几何模型、所述反应堆厂房仿真模型、核电厂布置信息和荷载信息，确定反应堆厂房的结构。
[0120]
可选的，所述第二确定模块202，包括：
[0121]
第一确定子模块，用于对所述几何模型中的所述固体部件和所述液体部件进行网格划分，获得液体部件网格；
[0122]
第二确定子模块，用于基于所述液体部件网格，按等参方向和预设粒子数目将所述液体部件转换为粒子，得到反应堆厂房仿真模型。
[0123]
可选的，所述第三确定模块204，包括：
[0124]
基于所述几何模型，选取预设位置；
[0125]
根据核电厂布置信息、荷载信息、所述几何模型和所述反应堆厂房仿真模型，获得所述预设位置的地震剪力时程曲线、所述固体部件的位移场、应力场和内力；
[0126]
根据所述预设位置的地震剪力时程曲线、所述固体部件的位移云图、应力云图和内力，确定反应堆厂房结构设计结果。
[0127]
可选的，所述第二确定子模块，包括：
[0128]
基于所述液体部件网格，按等参方向和预设粒子数目将所述液体部件转换为粒子，得到初始模型；
[0129]
基于初始模型，进行粒子数量敏感性分析，获得分析结果；
[0130]
根据所述分析结果，调整所述初始模型，获得中间模型；
[0131]
将所述中间模型作为初始模型，执行利用初始模型，进行粒子数量敏感性分析，获得分析结果以及根据所述分析结果，调整所述初始模型，获得中间模型的步骤，直到所述分析结果符合预设范围；
[0132]
将最后获得的中间模型确定为所述反应堆厂房仿真模型。
[0133]
考虑液体晃动的反应堆厂房结构确定装置200能够实现图1的方法实施例中电子设备实现的各个过程，为避免重复，这里不再赘述。
[0134]
参见图3，图3是本发明实施例提供的另一种电子设备的结构示意图。如图3所示，电子设备300包括：处理器301、存储器302及存储在所述存储器302上并可在所述处理器上运行的计算机程序，电子设备300中的各个组件通过总线系统303耦合在一起。可理解，总线系统303用于实现这些组件之间的连接通信。
[0135]
其中，处理器301，用于获取核电厂布置信息和荷载信息；
[0136]
根据核电厂布置信息和荷载信息，建立反应堆厂房结构的几何模型，所述几何模型包括固体部件和液体部件；
[0137]
基于所述几何模型，确定所述固体部件和所述液体部件的接触属性，且将所述液体部件转换为粒子，以得到反应堆厂房仿真模型；
[0138]
根据所述几何模型、所述反应堆厂房仿真模型、核电厂布置信息和荷载信息，确定反应堆厂房的结构。
[0139]
进一步的，处理器301，还用于对所述几何模型中的所述固体部件和所述液体部件进行网格划分，获得液体部件网格；
[0140]
基于所述液体部件网格，按等参方向和预设粒子数目将所述液体部件转换为粒子，得到反应堆厂房仿真模型。
[0141]
进一步的，处理器301，还用于基于所述几何模型，选取预设位置；
[0142]
根据核电厂布置信息、荷载信息、所述几何模型和所述反应堆厂房仿真模型，获得所述预设位置的地震剪力时程曲线、所述固体部件的位移场、应力场和内力；
[0143]
根据所述预设位置的地震剪力时程曲线、所述固体部件的位移云图、应力云图和内力，确定反应堆厂房结构设计结果。
[0144]
进一步的，处理器301，还用于基于所述液体部件网格，按等参方向和预设粒子数目将所述液体部件转换为粒子，得到初始模型；
[0145]
基于初始模型，进行粒子数量敏感性分析，获得分析结果；
[0146]
根据所述分析结果，调整所述初始模型，获得中间模型；
[0147]
将所述中间模型作为初始模型，执行利用初始模型，进行粒子数量敏感性分析，获得分析结果以及根据所述分析结果，调整所述初始模型，获得中间模型的步骤，直到所述分析结果符合预设范围；
[0148]
将所述中间模型确定为所述反应堆厂房仿真模型。
[0149]
电子设备300能够实现前述实施例中电子设备实现的各个过程，为避免重复，这里不再赘述。
[0150]
本发明实施例还提供一种计算机可读存储介质，计算机可读存储介质上存储有计算机程序，该计算机程序被处理器执行时实现上述考虑液体晃动的反应堆厂房结构确定方法实施例的各个过程，且能达到相同的技术效果，为避免重复，这里不再赘述。其中，所述的计算机可读存储介质，如只读存储器(read-only memory，简称rom)、随机存取存储器(random access memory，简称ram)、磁碟或者光盘等。
[0151]
需要说明的是，在本文中，术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含，从而使得包括一系列要素的过程、方法、物品或者装置不仅包括那些要素，而且还包括没有明确列出的其他要素，或者是还包括为这种过程、方法、物品或者装置所固有的要素。在没有更多限制的情况下，由语句“包括一个
……”
限定的要素，并不排除在包括该要素的过程、方法、物品或者装置中还存在另外的相同要素。
[0152]
通过以上的实施方式的描述，本领域的技术人员可以清楚地了解到上述实施例方法可借助软件加必需的通用硬件平台的方式来实现，当然也可以通过硬件，但很多情况下前者是更佳的实施方式。基于这样的理解，本发明的技术方案本质上或者说对现有技术做出贡献的部分可以以软件产品的形式体现出来，该计算机软件产品存储在一个存储介质
(如rom/ram、磁碟、光盘)中，包括若干指令用以使得一台终端(可以是手机，计算机，服务器，或者网络设备等)执行本发明各个实施例所述的方法。
[0153]
以上所述，仅为本发明的具体实施方式，但本发明的保护范围并不局限于此，任何熟悉本技术领域的技术人员在本发明揭露的技术范围内，可轻易想到变化或替换，都应涵盖在本发明的保护范围之内。因此，本发明的保护范围应以权利要求的保护范围为准。