(12)发明专利申请
(10)申请公布号 CN 112668072 A(43)申请公布日 2021.04.16
(21)申请号 202011279107.6(22)申请日 2020.11.16
(71)申请人 西南交通大学
地址 610000 四川省成都市二环路北一段(72)发明人 李鹏 代成栋 张德春 周跃民
郭严 杨翊仁 刘欢 杨钰莹 李伟才 (74)专利代理机构 成都慕川专利代理事务所
(普通合伙) 51278
代理人 谢芳(51)Int.Cl.
G06F 30/13(2020.01)G06F 30/23(2020.01)G06F 30/28(2020.01)G06F 119/14(2020.01)
权利要求书2页 说明书9页 附图4页
(54)发明名称
用于核反应堆燃料组件抗震分析的建模方法
(57)摘要
本发明具体涉及一种用于核反应堆燃料组件抗震分析的建模方法，其特征在于包括如下步骤：步骤1.进行运动特征进行分类；步骤2.将多种结构类型根据分类进行组合；步骤3.对步骤2中得到的组合进行分类设计；步骤4.总结出组件建模的特点，梳理出其中两类重复的两类基本子结构，发展相应的两类组合单元类型；步骤5.将步骤4中涉及到的两类组合单元类型分为用于模拟可以发生横向振动组件各根组件模型和用于模拟计算组件外侧的围板的有限元模型；步骤6.将步骤5中的单元模型进行组合方式设计；步骤7.选择合单元模型子结构将其进行拼接形成整体；步骤8.将步骤7中的子结构拼接耦合成一个整体，完成对核反应堆燃料组件抗震分析建模。
CN 112668072 ACN 112668072 A
权 利 要 求 书
1/2页
1.用于核反应堆燃料组件抗震分析的建模方法，其特征在于包括如下步骤：步骤1.对组件整体依据结构的局部和各部件的运动特征进行分类；步骤2.将步骤1中得到的多种结构类型根据分类进行组合；步骤3.对步骤2中得到的组合进行分类设计；步骤4.总结出组件建模的特点，梳理出其中两类重复的两类基本子结构，分别针对这两类结构发展相应的两类组合单元类型；
步骤5.将步骤4中涉及到的两类组合单元类型分为用于模拟可以发生横向振动组件各根组件模型的第一类组合单元有限元模型和用于模拟就算组件外侧的围板的第二类组合单元有限元模型；
步骤6.将步骤5中第一类组合单元模型和第二类组合的单元模型进行拼接；
并步骤7.按照组件个数任意选择第一类组合单元模型及第二类组合单元模型子结构，
将其进行拼接形成整体的全结构；
步骤8.采用共节点力传递的方式将步骤7中的各个子结构之间耦合成一个整体，完成对核反应堆燃料组件抗震分析的建模。
2.根据权利要求1所述用于核反应堆燃料组件抗震分析的建模方法，其特征在于：所述子结构类型包括固定不可发生横向振动的围板结构和可发生横向振动单根组件结构。
3.根据权利要求1所述用于核反应堆燃料组件抗震分析的建模方法，其特征在于：所述第一类复合型单元模型包括：梁单元模型、非线性弹性单元模型、线性弹性单元模型和，所述梁单元模型和非线性弹性单元模型为线性弹性单元和间隙单元组合而成。
4.根据权利要求1所述用于核反应堆燃料组件抗震分析的建模方法，其特征在于：所述第二类组合单元模型包括：非线性弹性单元模型，其由线性弹性单元和间隙单元组合而成；附加流体质量单元模型。
5.根据权利要求1所述用于核反应堆燃料组件抗震分析的建模方法，其特征在于：所述子结构类型包括固定不可发生横向振动的围板结构和可发生横向振动的单根组件结构。
6.根据权利要求3所述用于核反应堆燃料组件抗震分析的建模方法，其特征在于：所述梁单元模型、弹簧单元模型、非线性弹簧单元模型和附加流体质量单元的模型基础参数可以根据需求进行修改。
7.根据权利要求3所述用于核反应堆燃料组件抗震分析的建模方法，其特征在于：所述梁单元模型包含二维欧拉伯努利有限元梁模型、考虑剪切效应的有限元梁模型和轴向变形的有限元梁模型。
8.根据权利要求3所述用于核反应堆燃料组件抗震分析的建模方法，其特征在于：所述弹簧单元为线性弹簧约束的线性弹簧单元模型，其特征在于所述线性弹簧单元用于约束相邻两个节点的位移。
9.根据权利要求3所述用于核反应堆燃料组件抗震分析的建模方法，其特征在于：所述非线性弹簧单元由线性弹簧单元模型及间隙单元模型依靠端部节点的位移相耦合组合而成，所述非线性弹簧单元模型同时具备线性弹簧单元和间隙单元的作用，所述附加流体质量单元模型包括流体附加质量影响系数。
10.根据权利要求1所述用于核反应堆燃料组件抗震分析的建模方法，其特征在于：所述步骤8中的详细步骤为：采用自子结构建模会建立多个单根组件的振动方程，振动方程如
2
CN 112668072 A
权 利 要 求 书
2/2页
下：
……
将上述方程将依靠右端力进行耦合，依据力耦合及预估校正格式构建耦合迭代求解格式。
3
CN 112668072 A
说 明 书
用于核反应堆燃料组件抗震分析的建模方法
1/9页
技术领域
[0001]本发明属于新型能源电力生产领域，涉及核反应燃料组件抗震设计分析领域，具体涉及一种用于核反应堆燃料组件抗震分析的新型建模方法。
背景技术
[0002]核能作为清洁、高效和安全的绿色能源，具有储量丰富、能量密集度高、低碳无污染的优点，受到各国的极大重视。实现核电技术的完全自主研发，必定需要有完备的基础理论研究予以支撑。针对核电结构中的基础理论问题尤其是基础的力学问题进行研究，并以
对新一代反应堆的自主设计和发展至关重要。此储备核电设计技术，
[0003]燃料组件是反应堆内的重要部件。压水堆燃料组件的设计准则要求：在极限的事故工况下，如地震及冷却剂丧失工况，燃料组件应保持结构完整，确保堆芯可冷却；燃料组件的横向变形应不影响控制棒组件的插入和反应堆的安全停堆。因此燃料组件在正常运行和事故工况下的结构动力学特性将影响燃料组件的结构完整性，并直接影响反应堆运行的安全性和可靠性。反应堆燃料组件事故工况下的安全性已被各国列为重点安全审查项目，各核能设计机构都需要具备有相应自主开发的计算软件。针对燃料组件事故工况安全评估主要包括计算分析及实验验证两个方面，两者相辅相成互为支撑。横向计算分析在整个设计阶段起着重要作用，可以为燃料组件模型的设计、修改提供所需要的计算数据支持。现阶段针对燃料组件的计算分析主要借助于专用软件及商业软件进行，主要的计算模型可以分为：单个燃料组件横向详细模型、横向简化模型、横向碰撞模型及一排多个燃料组件的抗震计算分析模型，以下简称组件抗震模型，如图1所示。抗震模型的计算分析效率则强烈依赖于多个组件的建模方式。
[0004]目前组件抗震模型主要是将多个组件依次排列成一排而获得。各个组件之间依靠上下管座、围板边界条件及非线性碰撞力相互耦合。采用有限元方法对全部 15个组件进行“全结构”式的离散，如图2所示，通过施加边界条件而获得系统的运动微分方程组。考虑地震等极限工况下的等效外部激励，采用比例阻尼或者模态阻尼的施加系统阻尼，通过直接积分求解该微分方程组获得系统的响应。由于系统中含有碰撞等强非线性环节，需要采用较小的计算时间步进行内外交替迭代求解才能保证结果的收敛性，这显然要以牺牲计算效率为代价。
[0005]目前组件抗震模型主要采用繁琐的“全结构”式建模方法，该建模方式存在以下明显的缺陷和不足：[0006]①“全结构”的建模方式未考虑燃料组件结构的自相似特征，建模过程复杂，建模成本高；[0007]②“全结构”的建模方式会导致高维强非线性的求解效率和计算精度低的问题；[0008]③“全结构”的建模方式需要一次性建立全部结构的模型，不适合复杂工况的计算，复杂工况下的计算效率普遍偏低。
4
CN 112668072 A
说 明 书
2/9页
发明内容
[0009]针对上述现有技术中忽略组件结构存在的“自相似”特征而普遍采用“全结构”式的建模方式而导致建模过程繁复、建模效率低、适用性差、计算效率和计算精度差的问题，本发明充分利用了组件结构的局部特征，发展“子结构”式的建模策略，提出了核反应堆燃料组件抗震分析的建模新方法。
[0010]本发明提出的用于核反应堆燃料组件抗震分析的建模方法，其特征在于包括如下步骤：
[0011]步骤1.对组件整体依据结构的局部和各部件的运动特征进行分类；[0012]步骤2.将步骤1中得到的多种结构类型根据分类进行组合；[0013]步骤3.对步骤2中得到的组合进行分类设计；[0014]步骤4.总结出组件建模的特点，梳理出其中两类重复的两类基本子结构，分别针对这两类结构发展相应的两类组合单元类型；
[0015]步骤5.将步骤4中涉及到的两类组合单元类型分为用于模拟可以发生横向振动组件各根组件模型的第一类组合单元有限元模型和用于模拟计算组件外侧的围板的第二类组合单元有限元模型；
[0016]步骤6.将步骤5中第一类组合单元模型和第二类组合的单元模型进行组合方式设计；
[0017]步骤7.按照组件个数任意选择第一类组合单元模型及第二类组合单元模型子结构，并将其进行拼接形成整体的全结构；
[0018]步骤8.采用共节点力传递的方式将步骤7中的各个子结构拼接耦合成一个整体，完成对核反应堆燃料组件抗震分析的建模。[0019]所述第一类复合型单元模型包括：梁单元模型、非线性弹性单元模型、线性弹性单元模型和，所述梁单元模型和非线性弹性单元模型为线性弹性单元和间隙单元组合而成。[0020]所述第二类组合单元模型包括：非线性弹性单元模型，其由线性弹性单元和间隙单元组合而成；附加流体质量单元模型。
[0021]所述子结构类型包括固定不可发生横向振动的围板结构和可发生横向振动的单根组件结构。
[0022]所述梁单元模型、弹簧单元模型、非线性弹簧单元模型和附加流体质量单元的模型基础参数可以根据需求进行修改。
[0023]所述梁单元模型包含二维欧拉伯努利有限元梁模型、考虑剪切效应的有限元梁模型和轴向变形的有限元梁模型。
[0024]所述弹簧单元为线性弹簧约束的线性弹簧单元模型，其特征在于所述线性弹簧单元用于约束相邻两个节点的位移。
[0025]所述非线性弹簧单元由线性弹簧单元模型及间隙单元模型依靠端部节点的位移相耦合组合而成，所述非线性弹簧单元模型同时具备线性弹簧单元和间隙单元的作用。[0026]所述附加流体质量单元模型包括流体附加质量影响系数。[0027]所述系数为外部输入系数，程序内部不计算，只是读取后用于生成流体单元，该系数由其他软件计算而得到。
[0028]所述步骤8中的详细步骤为：采用子结构建模获得多个单根组件的振动方程，如
5
CN 112668072 A
说 明 书
3/9页
下：
[0029][0030][0031][0032]
……
将上述方程将依靠右端的非线性力进行耦合，依据力耦合及预估校正格式构建耦
合迭代求解格式。
[0034]本发明的有益效果：[0035]1.本发明由于采用“子结构”式的建模方式，会大大降低计算量，提高计算效率。[0036]2.本发明建模方式考虑结构局燃料组件的自相似特征，建模过程简单，建模成本低。
[0037]3.本发明不会导致高维强非线性的求解效率问题，提升了求解效率，增加了计算精度。
[0038]4.本发明不需要一次性建立全部结构的模型，适合复杂工况的计算，复杂工况下的计算效率高。
附图说明
[0039]图1为压水堆及堆内燃料组件示意图；[0040]图2为“全结构”建模方式建立的燃料组件横向排列模型示意图；[0041]图3为两类基本子结构单元‑‑‑组合类示意图；[0042]图4为有限元派生类示意图；
[0043]图5为各个子结构之间的力传递示意图；[0044]图6为各个子结构之间的拼接方式示意图；
[0045]图7为子结构振动之间的与预估校正迭代格式的流程图。
具体实施方式：[0046]实施例1：[0047]用于核反应堆燃料组件抗震分析的建模方法，其特征在于包括如下步骤：[0048]步骤1.对组件整体依据结构的局部和各部件的运动特征进行分类；[0049]步骤2.将步骤1中得到的多种结构类型根据分类进行组合；[0050]步骤3.对步骤2中得到的组合进行分类设计；[0051]步骤4.总结出组件建模的特点，梳理出其中两类重复的两类基本子结构，分别针对这两类结构发展相应的两类组合单元类型；
[0052]步骤5.将步骤4中涉及到的两类组合单元类型分为用于模拟可以发生横向振动组件各根组件模型的第一类组合单元有限元模型和用于模拟就算组件外侧的围板的第二类组合单元有限元模型；
[0053]步骤6.将步骤5中第一类组合单元模型和第二类组合的单元模型进行组合设计；[0054]步骤7.按照组件个数任意选择第一类组合单元模型及第二类组合单元模型子结构，并将其进行拼接形成整体的全结构；
6
[0033]
CN 112668072 A[0055]
说 明 书
4/9页
步骤8.采用共节点力传递的方式将步骤7中的各个子结构之间耦合成一个整体，
完成对核反应堆燃料组件抗震分析的建模。[0056]所述第一类复合型单元模型包括：梁单元模型、非线性弹性单元模型、线性弹性单元模型和，所述梁单元模型和非线性弹性单元模型为线性弹性单元和间隙单元组合而成。[0057]所述第二类组合单元模型包括：非线性弹性单元模型，其由线性弹性单元和间隙单元组合而成；附加流体质量单元模型。
[0058]所述子结构类型包括固定不可发生横向振动的围板结构和可发生横向振动的单根组件结构。
[0059]所述梁单元模型、弹簧单元模型、非线性弹簧单元模型和附加流体质量单元的模型基础参数可以根据需求进行修改。
[0060]所述梁单元模型包含二维欧拉伯努利有限元梁模型、考虑剪切效应的有限元梁模型和轴向变形的有限元梁模型。
[0061]所述弹簧单元为线性弹簧约束的线性弹簧单元模型，其特征在于所述线性弹簧单元用于约束相邻两个节点的位移。
[0062]所述非线性弹簧单元由线性弹簧单元模型及间隙单元模型依靠端部节点的位移相耦合组合而成，所述非线性弹簧单元模型同时具备线性弹簧单元和间隙单元的作用。[0063]所述附加流体质量单元模型包括流体附加质量影响系数。[0064]所述系数为外部输入系数，程序内部不计算，只是读取后用于生成流体单元，该系数由其他软件计算而得到。
[0065]所述步骤8中的详细步骤为：采用子结构建模获得多个单根组件的振动方程，如下：
[0066][0067][0068][0069]
……
将上述方程将依靠右端的非线性力进行耦合，依据力耦合及预估校正格式构建耦
合迭代求解格式。[0071]实施例2：
[0072]下面以具有n个组件的系统的建模为例说明计算量的变化，其中每个组件还有m个自由度。
[0073]当采用全结构建模方式，建立的方程为
[0070][0074][0075][0076][0077][0078][0079][0080]
矩阵M,C,K的维数均为(m×n)×(m×n)；而求解方程的计算量为
T0＝(m×n)×(m×n)当采用子结构建模时，需建立n个子结构方程，
……
7
CN 112668072 A[0081]
说 明 书
5/9页
这时，矩阵Ms,Cs,Ks的维数均为(m×m),而求解的整个方程的计算量为
[0083]T＝(m×m)
[0084]两种计算量比较[0085]T/T＝1/n0
[0086]由此可知，采用现有的“子结构”建模方式求解计算量仅为“全结构”建模计算的1/n。当建立包含15根组件的组件抗震工程计算模型时，“子结构”式的计算量仅为原来的1/15。采用相同的硬件计算可知，“子结构”式计算时间也为原来的1/15。因此，采用“子结构”式的建模方式，会大大降低计算量，提高计算效率。[0087]实施例3：[0088]用于核反应堆燃料组件抗震分析的建模方法，其特征在于包括如下步骤：[0089]步骤1.对组件整体依据结构的局部和各部件的运动特征进行分类；[0090]步骤2.将步骤1中得到的多种结构类型根据分类进行组合；[0091]步骤3.对步骤2中得到的组合进行分类设计；[0092]步骤4.总结出组件建模的特点，梳理出其中两类重复的两类基本子结构，分别针对这两类结构发展相应的两类组合单元类型；
[0093]步骤5.将步骤4中涉及到的两类组合单元类型分为用于模拟可以发生横向振动组件各根组件模型的第一类组合单元有限元模型和用于模拟就算组件外侧的围板的第二类组合单元有限元模型；
[0094]步骤6.将步骤5中第一类组合单元模型和第二类组合的单元模型进行拼接；
[0095]步骤7.按照组件个数任意选择第一类组合单元模型及第二类组合单元模型子结构，并将其进行拼接形成整体的全结构；
[0096]步骤8.采用共节点力传递的方式将步骤7中的各个子结构之间耦合成一个整体，完成对核反应堆燃料组件抗震分析的建模。[0097]所述第一类复合型单元模型包括：梁单元模型、非线性弹性单元模型、线性弹性单元模型和，所述梁单元模型和非线性弹性单元模型为线性弹性单元和间隙单元组合而成。[0098]所述第二类组合单元模型包括：非线性弹性单元模型，其由线性弹性单元和间隙单元组合而成；附加流体质量单元模型。
[0099]所述子结构类型包括固定不可发生横向振动的围板结构和可发生横向振动的单根组件结构。
[0100]所述梁单元模型、弹簧单元模型、非线性弹簧单元模型和附加流体质量单元的模型基础参数可以根据需求进行修改。
[0101]所述梁单元模型包含二维欧拉伯努利有限元梁模型、考虑剪切效应的有限元梁模型和轴向变形的有限元梁模型。
[0102]所述弹簧单元为线性弹簧约束的线性弹簧单元模型，其特征在于所述线性弹簧单元用于约束相邻两个节点的位移。
[0103]所述附加流体质量单元模型包括流体附加质量影响系数。[0104]所述系数为外部输入系数，程序内部不计算，只是读取后用于生成流体单元，该系数由其他软件计算而得到。
[0082]
8
CN 112668072 A[0105]
说 明 书
6/9页
所述步骤8中的详细步骤为：采用子结构建模获得多个单根组件的振动方程，如
下：
[0106][0107][0108][0109]
……
依据力耦合及预估校正格式构建耦将上述方程将依靠右端的非线性力进行耦合，
合迭代求解格式。[0111]实施例4：[0112]用于核反应堆燃料组件抗震分析的建模方法，其特征在于包括如下步骤：[0113]步骤1.对组件整体依据结构的局部和各部件的运动特征进行分类；[0114]步骤2.将步骤1中得到的多种结构类型根据分类进行组合；[0115]步骤3.对步骤2中得到的组合进行分类设计；[0116]步骤4.总结出组件建模的特点，梳理出其中两类重复的两类基本子结构，分别针对这两类结构发展相应的两类组合单元类型；
[0117]步骤5.将步骤4中涉及到的两类组合单元类型分为用于模拟可以发生横向振动组件各根组件模型的第一类组合单元有限元模型和用于模拟就算组件外侧的围板的第二类组合单元有限元模型；
[0118]步骤6.将步骤5中第一类组合单元模型和第二类组合的单元模型进行拼接；
[0119]步骤7.按照组件个数任意选择第一类组合单元模型及第二类组合单元模型子结构，并将其进行拼接形成整体的全结构；
[0120]步骤8.采用共节点力传递的方式将步骤7中的各个子结构之间耦合成一个整体，完成对核反应堆燃料组件抗震分析的建模。[0121]所述第一类复合型单元模型包括：梁单元模型、非线性弹性单元模型、线性弹性单元模型和，所述梁单元模型和非线性弹性单元模型为线性弹性单元和间隙单元组合而成。[0122]所述第二类组合单元模型包括：非线性弹性单元模型，其由线性弹性单元和间隙单元组合而成；附加流体质量单元模型。
[0123]所述子结构类型包括固定不可发生横向振动的围板结构和可发生横向振动的单根组件结构。
[0124]所述梁单元模型、弹簧单元模型、非线性弹簧单元模型和附加流体质量单元的模型基础参数可以根据需求进行修改。
[0125]所述梁单元模型包含二维欧拉伯努利有限元梁模型、考虑剪切效应的有限元梁模型和轴向变形的有限元梁模型。
[0126]所述弹簧单元为线性弹簧约束的线性弹簧单元模型，其特征在于所述线性弹簧单元用于约束相邻两个节点的位移。
[0127]所述非线性弹簧单元由线性弹簧单元模型及间隙单元模型依靠端部节点的位移相耦合组合而成，所述非线性弹簧单元模型同时具备线性弹簧单元和间隙单元的作用。[0128]所述附加流体质量单元模型包括流体附加质量影响系数。[0129]所述步骤8中的详细步骤为：采用子结构建模获得多个单根组件的振动方程，如
[0110]
9
CN 112668072 A
说 明 书
7/9页
下：
[0130][0131][0132][0133]
……
将上述方程将依靠右端的非线性力进行耦合，依据力耦合及预估校正格式构建耦
合迭代求解格式。[0135]实施例5：[0136]用于核反应堆燃料组件抗震分析的建模方法，其特征在于包括如下步骤：[0137]步骤1.对组件整体依据结构的局部和各部件的运动特征进行分类；[0138]步骤2.将步骤1中得到的多种结构类型根据分类进行组合；[0139]步骤3.对步骤2中得到的组合进行分类设计；[0140]步骤4.总结出组件建模的特点，梳理出其中两类重复的两类基本子结构，分别针对这两类结构发展相应的两类组合单元类型；
[0141]步骤5.将步骤4中涉及到的两类组合单元类型分为用于模拟可以发生横向振动组件各根组件模型的第一类组合单元有限元模型和用于模拟就算组件外侧的围板的第二类组合单元有限元模型；
[0142]步骤6.将步骤5中第一类组合单元模型和第二类组合的单元模型进行拼接；
[0143]步骤7.按照组件个数任意选择第一类组合单元模型及第二类组合单元模型子结构，并将其进行拼接形成整体的全结构；
[0144]步骤8.采用共节点力传递的方式将步骤7中的各个子结构之间耦合成一个整体，完成对核反应堆燃料组件抗震分析的建模。[0145]所述第一类复合型单元模型包括：梁单元模型、非线性弹性单元模型、线性弹性单元模型和，所述梁单元模型和非线性弹性单元模型为线性弹性单元和间隙单元组合而成。[0146]所述第二类组合单元模型包括：非线性弹性单元模型，其由线性弹性单元和间隙单元组合而成；附加流体质量单元模型。
[0147]所述子结构类型包括固定不可发生横向振动的围板结构和可发生横向振动的单根组件结构。
[0148]所述梁单元模型、弹簧单元模型、非线性弹簧单元模型和附加流体质量单元的模型基础参数可以根据需求进行修改。
[0149]所述梁单元模型包含二维欧拉伯努利有限元梁模型、考虑剪切效应的有限元梁模型和轴向变形的有限元梁模型。
[0150]所述弹簧单元为线性弹簧约束的线性弹簧单元模型，其特征在于所述线性弹簧单元用于约束相邻两个节点的位移。
[0151]所述非线性弹簧单元由线性弹簧单元模型及间隙单元模型依靠端部节点的位移相耦合组合而成，所述非线性弹簧单元模型同时具备线性弹簧单元和间隙单元的作用。[0152]所述附加流体质量单元模型包括流体附加质量影响系数。[0153]所述系数为外部输入系数，程序内部不计算，只是读取后用于生成流体单元，该系
[0134]
10
CN 112668072 A
说 明 书
8/9页
数由其他软件计算而得到。[0154]实施例6：[0155]用于核反应堆燃料组件抗震分析的建模方法，其特征在于包括如下步骤：[0156]步骤1.对组件整体依据结构的局部和各部件的运动特征进行分类；[0157]步骤2.将步骤1中得到的多种结构类型根据分类进行组合；[0158]步骤3.对步骤2中得到的组合进行分类设计；[0159]步骤4.总结出组件建模的特点，梳理出其中两类重复的两类基本子结构，分别针对这两类结构发展相应的两类组合单元类型；
[0160]步骤5.将步骤4中涉及到的两类组合单元类型分为用于模拟可以发生横向振动组件各根组件模型的第一类组合单元有限元模型和用于模拟就算组件外侧的围板的第二类组合单元有限元模型；
[0161]步骤6.将步骤5中第一类组合单元模型和第二类组合的单元模型进行拼接；
[0162]步骤7.按照组件个数任意选择第一类组合单元模型及第二类组合单元模型子结构，并将其进行拼接形成整体的全结构；
[0163]步骤8.采用共节点力传递的方式将步骤7中的各个子结构之间耦合成一个整体，完成对核反应堆燃料组件抗震分析的建模。[0164]所述第一类复合型单元模型包括：梁单元模型、非线性弹性单元模型、线性弹性单元模型和，所述梁单元模型和非线性弹性单元模型为线性弹性单元和间隙单元组合而成。[0165]所述第二类组合单元模型包括：非线性弹性单元模型，其由线性弹性单元和间隙单元组合而成；附加流体质量单元模型。
[0166]所述子结构类型包括固定不可发生横向振动的围板结构和可发生横向振动的单根组件结构。
[0167]所述梁单元模型、弹簧单元模型、非线性弹簧单元模型和附加流体质量单元的模型基础参数可以根据需求进行修改。
[0168]所述梁单元模型包含二维欧拉伯努利有限元梁模型、考虑剪切效应的有限元梁模型和轴向变形的有限元梁模型。
[0169]所述弹簧单元为线性弹簧约束的线性弹簧单元模型，其特征在于所述线性弹簧单元用于约束相邻两个节点的位移。
[0170]所述非线性弹簧单元由线性弹簧单元模型及间隙单元模型依靠端部节点的位移相耦合组合而成，所述非线性弹簧单元模型同时具备线性弹簧单元和间隙单元的作用。[0171]所述附加流体质量单元模型包括流体附加质量影响系数。[0172]所述系数为外部输入系数，程序内部不计算，只是读取后用于生成流体单元，该系数由其他软件计算而得到。
[0173]所述步骤8中的详细步骤为：采用子结构建模获得多个单根组件的振动方程，如下：
[0174][0175][0176][0177]
11
……
CN 112668072 A[0178]
说 明 书
9/9页
将上述方程将依靠右端的非线性力进行耦合，依据力耦合及预估校正格式构建耦
合迭代求解格式。
12
CN 112668072 A
说 明 书 附 图
1/4页
图1
图2
13
CN 112668072 A
说 明 书 附 图
2/4页
图3
图4
14
CN 112668072 A
说 明 书 附 图
3/4页
图5
图6
15
CN 112668072 A
说 明 书 附 图
4/4页
图7
16
因篇幅问题不能全部显示，请点此查看更多更全内容