CN112836269A - Method for splicing nuclear reactor fuel assembly anti-seismic analysis modeling substructure types - Google Patents

Method for splicing nuclear reactor fuel assembly anti-seismic analysis modeling substructure types Download PDF

Info

Publication number
CN112836269A
CN112836269A CN202011279106.1A CN202011279106A CN112836269A CN 112836269 A CN112836269 A CN 112836269A CN 202011279106 A CN202011279106 A CN 202011279106A CN 112836269 A CN112836269 A CN 112836269A
Authority
CN
China
Prior art keywords
model
types
unit
fuel assembly
modeling
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Pending
Application number
CN202011279106.1A
Other languages
Chinese (zh)
Inventor
李鹏
代成栋
张德春
周跃民
郭严
杨翊仁
刘欢
杨钰莹
李伟才
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Southwest Jiaotong University
Original Assignee
Southwest Jiaotong University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Southwest Jiaotong University filed Critical Southwest Jiaotong University
Priority to CN202011279106.1A priority Critical patent/CN112836269A/en
Publication of CN112836269A publication Critical patent/CN112836269A/en
Pending legal-status Critical Current

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/10Geometric CAD
    • G06F30/13Architectural design, e.g. computer-aided architectural design [CAAD] related to design of buildings, bridges, landscapes, production plants or roads
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • G06F30/23Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • G06F30/28Design optimisation, verification or simulation using fluid dynamics, e.g. using Navier-Stokes equations or computational fluid dynamics [CFD]
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/14Force analysis or force optimisation, e.g. static or dynamic forces

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Geometry (AREA)
  • General Physics & Mathematics (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • General Engineering & Computer Science (AREA)
  • Mathematical Analysis (AREA)
  • Pure & Applied Mathematics (AREA)
  • Mathematical Optimization (AREA)
  • Fluid Mechanics (AREA)
  • Computing Systems (AREA)
  • Mathematical Physics (AREA)
  • Algebra (AREA)
  • Architecture (AREA)
  • Civil Engineering (AREA)
  • Structural Engineering (AREA)
  • Computational Mathematics (AREA)
  • Monitoring And Testing Of Nuclear Reactors (AREA)

Abstract

The invention particularly relates to a method for splicing nuclear reactor fuel assembly anti-seismic analysis modeling substructure types, which is characterized by comprising the following steps of: step 1, summarizing the characteristics of component modeling, combing two repeated types of basic substructures, and developing two corresponding types of combined unit types aiming at the two types of structures respectively; step 2, dividing the types of the two types of combined units into a model for simulating each component of the component capable of generating transverse vibration and a coaming for simulating and calculating the outer side of the component; step 3, designing the two types of combined unit models in the step 2 in a combined mode; step 4, randomly selecting a combined unit model type substructure according to the number of the components, and splicing the combined unit model type substructure to form an integral whole structure; and 5, splicing and coupling the sub-structures in the step 4 into a whole in a mode of common-node force transmission to complete modeling of the nuclear reactor fuel assembly anti-seismic analysis.

Description

Method for splicing nuclear reactor fuel assembly anti-seismic analysis modeling substructure types
Technical Field
The invention belongs to the field of novel energy power production, relates to the field of nuclear reaction fuel assembly anti-seismic design analysis, and particularly relates to a method for splicing nuclear reactor fuel assembly anti-seismic analysis modeling substructure types.
Background
The nuclear energy is used as clean, efficient and safe green energy, has the advantages of abundant reserves, high energy intensity, low carbon and no pollution, and is greatly valued by various countries. The complete independent research and development of the nuclear power technology are realized, and complete basic theory research is required to support the technology. The method aims at the basic theoretical problem in a nuclear power structure, particularly the basic mechanical problem, researches are carried out, and a nuclear power design technology is reserved according to the basic theoretical problem, so that the method is important for the independent design and development of a new-generation reactor.
Fuel assemblies are important components within the reactor. Design criteria for pressurized water reactor fuel assemblies require: under extreme accident conditions, such as earthquakes and coolant loss conditions, the fuel assemblies should maintain structural integrity to ensure that the reactor core can be cooled; lateral deformation of the fuel assembly should not affect the insertion of the control rod assembly and safe shutdown of the reactor. The structural dynamics of the fuel assembly under normal operation and accident conditions will therefore affect the structural integrity of the fuel assembly and directly affect the safety and reliability of the reactor operation. The safety of the reactor fuel assembly under the accident condition is already listed as a key safety examination project by various countries, and each nuclear energy design mechanism needs to be provided with corresponding independently developed computing software. The safety assessment method for the accident condition of the fuel assembly mainly comprises two aspects of calculation analysis and experimental verification, wherein the two aspects supplement each other to support each other. The transverse calculation analysis plays an important role in the whole design stage, and can provide required calculation data support for the design and modification of the fuel assembly model. At present, the computational analysis of the fuel assembly is mainly performed by means of special software and commercial software, and the main computational models can be divided into: the method comprises the following steps of transversely detailed models, transversely simplified models, transversely collision models and seismic calculation analysis models of a row of multiple fuel assemblies, which are hereinafter referred to as assembly seismic models, and is shown in FIG. 1. The efficiency of the computational analysis of the seismic model is strongly dependent on the way the components are modeled.
At present, the component anti-seismic model is mainly obtained by arranging a plurality of components in a row in sequence. The components are mutually coupled by the boundary conditions of the upper and lower pipe seats and the enclosing plate and the nonlinear impact force. The finite element method is adopted to carry out 'full structure' type dispersion on all 15 components, and as shown in figure 2, a motion differential equation system of the system is obtained by applying boundary conditions. And (3) considering equivalent external excitation under extreme working conditions such as earthquake and the like, applying system damping by adopting proportional damping or modal damping, and solving the differential equation set through direct integration to obtain the response of the system. Because the system contains strong nonlinear links such as collision, the convergence of the result can be ensured only by adopting a smaller calculation time step to carry out inner and outer alternative iterative solution, which obviously takes the cost of the calculation efficiency as a price.
At present, the anti-seismic component model mainly adopts a fussy 'full structure' modeling method, and the modeling method has the following obvious defects and shortcomings:
firstly, the self-similarity characteristic of the fuel assembly structure is not considered in the modeling mode of the 'full structure', the modeling process is complex, and the modeling cost is high;
secondly, the modeling mode of the 'full structure' can cause the problems of low solving efficiency and low calculation precision of high-dimensional strong nonlinearity;
and thirdly, the modeling mode of the 'full structure' needs to establish models of all structures at one time, is not suitable for the calculation of complex working conditions, and has low calculation efficiency under the complex working conditions.
Disclosure of Invention
Aiming at the problems of complex modeling process, low modeling efficiency, poor applicability, and poor calculation efficiency and calculation accuracy caused by neglecting self-similarity characteristics of the component structure and generally adopting a full-structure modeling mode in the prior art, the invention fully utilizes local characteristics of the component structure, develops a substructure-type modeling strategy, and provides a new modeling method for nuclear reactor fuel component seismic analysis.
The invention provides a method for splicing nuclear reactor fuel assembly anti-seismic analysis modeling substructure types, which is characterized by comprising the following steps of:
step 1, summarizing the characteristics of component modeling, combing two repeated types of basic substructures, and developing two corresponding types of combined unit types aiming at the two types of structures respectively;
step 2, dividing the two types of combination unit types involved in the step 1 into a first type of combination unit finite element model for simulating each component model of the component capable of generating transverse vibration and a second type of combination unit finite element model for simulating and calculating a coaming outside the component;
step 3, designing a combination mode of the first-type combination unit model and the second-type combination unit model in the step 2;
step 4, randomly selecting substructures of the first type of combination unit model and the second type of combination unit model according to the number of the components, and splicing the substructures to form an integral whole structure;
and 5, splicing and coupling the sub-structures in the step 4 into a whole in a mode of common-node force transmission to complete modeling of the nuclear reactor fuel assembly anti-seismic analysis.
The first type composite unit model comprises: the beam unit model and the nonlinear elastic unit model are formed by combining a linear elastic unit and a gap unit.
The second type of combined unit model comprises: a nonlinear elastic cell model composed of a linear elastic cell and a gap cell; a fluid mass cell model is appended.
The substructure types include a structure for fixing enclosing plates which cannot generate transverse vibration and a structure for generating single assembly which can generate transverse vibration.
The model base parameters of the beam element model, the spring element model, the nonlinear spring element model and the additional fluid mass element can be modified as required.
The beam element model comprises a two-dimensional Euler Bernoulli finite element beam model, a finite element beam model considering a shearing effect and an axially deformed finite element beam model.
The spring unit is a linear spring unit model constrained by a linear spring, and is characterized in that the linear spring unit is used for constraining the displacement of two adjacent nodes.
The nonlinear spring unit is formed by coupling and combining a linear spring unit model and a gap unit model by means of displacement of end node, and the nonlinear spring unit model has the functions of the linear spring unit and the gap unit at the same time.
The additional fluid mass cell model includes a fluid additional mass influence coefficient.
The coefficients are externally input coefficients, are not calculated inside the program, but are read and used for generating the fluid unit, and are calculated by other software.
The detailed steps in the step 5 are as follows: obtaining vibration equations of a plurality of single assemblies by adopting substructure modeling, wherein the vibration equations comprise the following steps:
Figure RE-GDA0002959275620000031
Figure RE-GDA0002959275620000032
……
Figure RE-GDA0002959275620000033
the equations are coupled by means of nonlinear force at the right end, and a coupling iterative solution format is constructed according to the force coupling and the estimation correction format.
The invention has the beneficial effects that:
1. the invention adopts a sub-structure type modeling mode, so that the calculation amount can be greatly reduced, and the calculation efficiency can be improved.
2. The modeling mode of the invention considers the self-similarity characteristic of the fuel assembly of the structural office, the modeling process is simple, and the modeling cost is low.
3. The invention can not cause the problem of solving efficiency of high-dimensional strong nonlinearity, improves the solving efficiency and increases the calculation precision.
4. The method does not need to establish models of all structures at one time, is suitable for calculation under complex working conditions, and has high calculation efficiency under the complex working conditions.
Drawings
FIG. 1 is a schematic diagram of a pressurized water reactor and in-reactor fuel assemblies;
FIG. 2 is a schematic diagram of a fuel assembly transverse arrangement model established by a full structure modeling manner;
FIG. 3 is a schematic diagram of two types of basic sub-structural units- -combination type;
FIG. 4 is a schematic diagram of finite element derivation classes;
FIG. 5 is a schematic diagram of the force transfer between the various substructures;
FIG. 6 is a schematic diagram of the manner of stitching between the various substructures;
FIG. 7 is a flow chart of an iterative format of AND estimation correction between substructure vibrations.
The specific implementation mode is as follows:
example 1:
a method for splicing nuclear reactor fuel assembly seismic analysis modeling substructure types, comprising the steps of:
step 1, summarizing the characteristics of component modeling, combing two repeated types of basic substructures, and developing two corresponding types of combined unit types aiming at the two types of structures respectively;
step 2, dividing the two types of combination unit types involved in the step 1 into a first type of combination unit finite element model for simulating each component model of the component capable of generating transverse vibration and a second type of combination unit finite element model for simulating and calculating a coaming outside the component;
step 3, designing a combination mode of the first-type combination unit model and the second-type combination unit model in the step 2;
step 4, randomly selecting substructures of the first type of combination unit model and the second type of combination unit model according to the number of the components, and splicing the substructures to form an integral whole structure;
and 5, splicing and coupling the sub-structures in the step 4 into a whole in a mode of common-node force transmission to complete modeling of the nuclear reactor fuel assembly anti-seismic analysis.
The first type composite unit model comprises: a two-dimensional Euler Bernoulli beam finite element model, a finite element beam model considering a shearing effect, and a finite element beam model deforming in the axial direction; a nonlinear elastic cell model composed of a linear elastic cell and a gap cell; linear elastic cell model, additional fluid mass cell model.
The second type of combined unit model comprises a nonlinear elastic unit model, a second type of combined unit model and a third type of combined unit model, wherein the nonlinear elastic unit model is formed by combining linear elastic units and gap units; a fluid mass cell model is appended.
The substructure types include a structure for fixing enclosing plates which cannot generate transverse vibration and a structure for generating single assembly which can generate transverse vibration.
The model base parameters of the beam element model, the spring element model, the nonlinear spring element model and the additional fluid mass element can be modified as required.
The beam element includes a two-dimensional Euler Bernoulli finite element beam model, a finite element beam model considering shear effect, and an axially deformed finite element beam model.
The spring unit is a linear spring unit model constrained by a linear spring, and is characterized in that the linear spring unit is used for constraining the displacement of two adjacent nodes.
The nonlinear spring unit is formed by coupling and combining a linear spring unit model and a gap unit model by means of displacement of end node; the nonlinear spring element model may function as both a linear spring element and a gap element.
The additional fluid mass cell model includes a fluid additional mass influence coefficient.
The coefficients are externally input coefficients, are not calculated inside the program, but are read and used for generating the fluid unit, and are calculated by other software.
The detailed steps in the step 5 are as follows: obtaining vibration equations of a plurality of single assemblies by adopting substructure modeling, wherein the vibration equations comprise the following steps:
Figure RE-GDA0002959275620000051
Figure RE-GDA0002959275620000052
……
Figure RE-GDA0002959275620000053
the equations are coupled by means of nonlinear force at the right end, and a coupling iterative solution format is constructed according to the force coupling and the estimation correction format.
Example 2:
the variation of the calculated quantities is illustrated below by way of example of a modeling of a system with n components, each component having m degrees of freedom.
When a full-structure modeling mode is adopted, the established equation is
Figure RE-GDA0002959275620000061
The dimensions of the matrices M, C, K are all (M × n) × (M × n); and the amount of calculation to solve the equation is
T0=(m×n)×(m×n)
When adopting the substructure modeling, n substructure equations need to be established,
Figure RE-GDA0002959275620000062
Figure RE-GDA0002959275620000063
……
Figure RE-GDA0002959275620000064
at this time, the matrix Ms,Cs,KsAre all (m × m), and the whole equation to be solved is calculated as
T=(m×m)
Comparison of two calculated quantities
T/T0=1/n
Therefore, the solving calculation amount is only 1/n of the modeling calculation of the 'full structure' by adopting the existing 'substructure' modeling mode. When a component earthquake-resistant engineering calculation model containing 15 components is built, the calculation amount of the 'substructure' formula is only 1/15 of the original formula. The same hardware calculation shows that the calculation time of the "substructure" formula is also 1/15. Therefore, the adoption of the sub-structure modeling mode can greatly reduce the calculation amount and improve the calculation efficiency.
Example 3:
a method for splicing nuclear reactor fuel assembly seismic analysis modeling substructure types, comprising the steps of:
step 1, summarizing the characteristics of component modeling, combing two repeated types of basic substructures, and developing two corresponding types of combined unit types aiming at the two types of structures respectively;
step 2, dividing the two types of combination unit types involved in the step 1 into a first type of combination unit finite element model for simulating each component model of the component capable of generating transverse vibration and a second type of combination unit finite element model for simulating and calculating a coaming outside the component;
step 3, designing a combination mode of the first-type combination unit model and the second-type combination unit model in the step 2;
step 4, randomly selecting substructures of the first type of combination unit model and the second type of combination unit model according to the number of the components, and splicing the substructures to form an integral whole structure;
and 5, splicing and coupling the sub-structures in the step 4 into a whole in a mode of common-node force transmission to complete modeling of the nuclear reactor fuel assembly anti-seismic analysis.
The first type composite unit model comprises: a two-dimensional Euler Bernoulli beam finite element model, a finite element beam model considering a shearing effect, and a finite element beam model deforming in the axial direction; a nonlinear elastic cell model composed of a linear elastic cell and a gap cell; linear elastic cell model, additional fluid mass cell model.
The second type of combined unit model comprises: a nonlinear elastic cell model composed of a linear elastic cell and a gap cell; a fluid mass cell model is appended.
The model base parameters of the beam element model, the spring element model, the nonlinear spring element model and the additional fluid mass element can be modified as required.
The beam element includes a two-dimensional Euler Bernoulli finite element beam model, a finite element beam model considering shear effect, and an axially deformed finite element beam model.
The additional fluid mass cell model includes a fluid additional mass influence coefficient.
The coefficients are externally input coefficients, are not calculated inside the program, but are read and used for generating the fluid unit, and are calculated by other software.
The detailed steps in the step 5 are as follows: obtaining vibration equations of a plurality of single assemblies by adopting substructure modeling, wherein the vibration equations comprise the following steps:
Figure RE-GDA0002959275620000071
Figure RE-GDA0002959275620000072
……
Figure RE-GDA0002959275620000073
the equations are coupled by means of nonlinear force at the right end, and a coupling iterative solution format is constructed according to the force coupling and the estimation correction format.
Example 4:
a method for splicing nuclear reactor fuel assembly seismic analysis modeling substructure types, comprising the steps of:
step 1, summarizing the characteristics of component modeling, combing two repeated types of basic substructures, and developing two corresponding types of combined unit types aiming at the two types of structures respectively;
step 2, dividing the two types of combination unit types involved in the step 1 into a first type of combination unit finite element model for simulating each component model of the component capable of generating transverse vibration and a second type of combination unit finite element model for simulating and calculating a coaming outside the component;
step 3, designing a combination mode of the first-type combination unit model and the second-type combination unit model in the step 2;
step 4, randomly selecting substructures of the first type of combination unit model and the second type of combination unit model according to the number of the components, and splicing the substructures to form an integral whole structure;
and 5, splicing and coupling the sub-structures in the step 4 into a whole in a mode of common-node force transmission to complete modeling of the nuclear reactor fuel assembly anti-seismic analysis.
The first type composite unit model comprises: a two-dimensional Euler Bernoulli beam finite element model, a finite element beam model considering a shearing effect, and a finite element beam model deforming in the axial direction; a nonlinear elastic cell model composed of a linear elastic cell and a gap cell; linear elastic cell model, additional fluid mass cell model.
The second type of combined unit model comprises: a nonlinear elastic cell model composed of a linear elastic cell and a gap cell; a fluid mass cell model is appended.
The substructure types include a structure for fixing enclosing plates which cannot generate transverse vibration and a structure for generating single assembly which can generate transverse vibration.
The model base parameters of the beam element model, the spring element model, the nonlinear spring element model and the additional fluid mass element can be modified as required.
The beam element includes a two-dimensional Euler Bernoulli finite element beam model, a finite element beam model considering shear effect, and an axially deformed finite element beam model.
The spring unit is a linear spring unit model constrained by a linear spring, and is characterized in that the linear spring unit is used for constraining the displacement of two adjacent nodes.
The nonlinear spring unit is formed by coupling and combining a linear spring unit model and a gap unit model by means of displacement of end node; the nonlinear spring element model may function as both a linear spring element and a gap element.
The additional fluid mass cell model includes a fluid additional mass influence coefficient.
The detailed steps in the step 5 are as follows: obtaining vibration equations of a plurality of single assemblies by adopting substructure modeling, wherein the vibration equations comprise the following steps:
Figure RE-GDA0002959275620000091
Figure RE-GDA0002959275620000092
……
Figure RE-GDA0002959275620000093
the equations are coupled by means of nonlinear force at the right end, and a coupling iterative solution format is constructed according to the force coupling and the estimation correction format.
Example 5:
step 1, summarizing the characteristics of component modeling, combing two repeated types of basic substructures, and developing two corresponding types of combined unit types aiming at the two types of structures respectively;
step 2, dividing the two types of combination unit types involved in the step 1 into a first type of combination unit finite element model for simulating each component model of the component capable of generating transverse vibration and a second type of combination unit finite element model for simulating and calculating a coaming outside the component;
step 3, designing a combination mode of the first-type combination unit model and the second-type combination unit model in the step 2;
step 4, randomly selecting substructures of the first type of combination unit model and the second type of combination unit model according to the number of the components, and splicing the substructures to form an integral whole structure;
and 5, splicing and coupling the sub-structures in the step 4 into a whole in a mode of common-node force transmission to complete modeling of the nuclear reactor fuel assembly anti-seismic analysis.
The first type composite unit model comprises: a two-dimensional Euler Bernoulli beam finite element model, a finite element beam model considering a shearing effect, and a finite element beam model deforming in the axial direction; a nonlinear elastic cell model composed of a linear elastic cell and a gap cell; linear elastic cell model, additional fluid mass cell model.
The second type of combined unit model comprises: a nonlinear elastic cell model composed of a linear elastic cell and a gap cell; a fluid mass cell model is appended.
The substructure types include a structure for fixing enclosing plates which cannot generate transverse vibration and a structure for generating single assembly which can generate transverse vibration.
The model base parameters of the beam element model, the spring element model, the nonlinear spring element model and the additional fluid mass element can be modified as required.
The beam element includes a two-dimensional Euler Bernoulli finite element beam model, a finite element beam model considering shear effect, and an axially deformed finite element beam model.
The spring unit is a linear spring unit model constrained by a linear spring, and is characterized in that the linear spring unit is used for constraining the displacement of two adjacent nodes.
The nonlinear spring unit is formed by coupling and combining a linear spring unit model and a gap unit model by means of displacement of end node; the nonlinear spring element model may function as both a linear spring element and a gap element.
The additional fluid mass cell model includes a fluid additional mass influence coefficient.
The coefficients are externally input coefficients, are not calculated inside the program, but are read and used for generating the fluid unit, and are calculated by other software.
Example 6:
a method for splicing nuclear reactor fuel assembly seismic analysis modeling substructure types, comprising the steps of:
step 1, summarizing the characteristics of component modeling, combing two repeated types of basic substructures, and developing two corresponding types of combined unit types aiming at the two types of structures respectively;
step 2, dividing the two types of combination unit types involved in the step 1 into a first type of combination unit finite element model for simulating each component model of the component capable of generating transverse vibration and a second type of combination unit finite element model for simulating and calculating a coaming outside the component;
step 3, designing a combination mode of the first-type combination unit model and the second-type combination unit model in the step 2;
step 4, randomly selecting substructures of the first type of combination unit model and the second type of combination unit model according to the number of the components, and splicing the substructures to form an integral whole structure;
and 5, splicing and coupling the sub-structures in the step 4 into a whole in a mode of common-node force transmission to complete modeling of the nuclear reactor fuel assembly anti-seismic analysis.
The first type composite unit model comprises: a two-dimensional Euler Bernoulli beam finite element model, a finite element beam model considering a shearing effect, and a finite element beam model deforming in the axial direction; a nonlinear elastic cell model composed of a linear elastic cell and a gap cell; linear elastic cell model, additional fluid mass cell model.
The second type of combined unit model comprises a nonlinear elastic unit model, a second type of combined unit model and a third type of combined unit model, wherein the nonlinear elastic unit model is formed by combining linear elastic units and gap units; a fluid mass cell model is appended.
The substructure types include a structure for fixing enclosing plates which cannot generate transverse vibration and a structure for generating single assembly which can generate transverse vibration.
The model base parameters of the beam element model, the spring element model, the nonlinear spring element model and the additional fluid mass element can be modified as required.
The beam element includes a two-dimensional Euler Bernoulli finite element beam model, a finite element beam model considering shear effect, and an axially deformed finite element beam model.
The spring unit is a linear spring unit model constrained by a linear spring, and is characterized in that the linear spring unit is used for constraining the displacement of two adjacent nodes.
The nonlinear spring unit is formed by coupling and combining a linear spring unit model and a gap unit model by means of displacement of end node; the nonlinear spring element model may function as both a linear spring element and a gap element.
The additional fluid mass cell model includes a fluid additional mass influence coefficient.
The coefficients are externally input coefficients, are not calculated inside the program, but are read and used for generating the fluid unit, and are calculated by other software.
The detailed steps in the step 5 are as follows: obtaining vibration equations of a plurality of single assemblies by adopting substructure modeling, wherein the vibration equations comprise the following steps:
Figure RE-GDA0002959275620000111
Figure RE-GDA0002959275620000112
……
Figure RE-GDA0002959275620000113
the equations are coupled by means of nonlinear force at the right end, and a coupling iterative solution format is constructed according to the force coupling and the estimation correction format.

Claims (10)

1. A method for splicing nuclear reactor fuel assembly seismic analysis modeling substructure types, comprising the steps of:
step 1, summarizing the characteristics of component modeling, combing two repeated types of basic substructures, and developing two corresponding types of combined unit types aiming at the two types of structures respectively;
step 2, dividing the two types of combination unit types involved in the step 1 into a first type of combination unit finite element model for simulating each component model of the component capable of generating transverse vibration and a second type of combination unit finite element model for simulating and calculating a coaming outside the component;
step 3, designing a combination mode of the first-type combination unit model and the second-type combination unit model in the step 2;
step 4, randomly selecting substructures of the first type of combination unit model and the second type of combination unit model according to the number of the components, and splicing the substructures to form an integral whole structure;
and 5, splicing and coupling the sub-structures in the step 4 into a whole in a mode of common-node force transmission to complete modeling of the nuclear reactor fuel assembly anti-seismic analysis.
2. The method for splicing nuclear reactor fuel assembly seismic analysis modeling substructure types according to claim 1, characterized in that: the first type composite unit model comprises: the beam unit model and the nonlinear elastic unit model are formed by combining a linear elastic unit and a gap unit.
3. The method for splicing nuclear reactor fuel assembly seismic analysis modeling substructure types according to claim 1, characterized in that: the second type of combined unit model comprises: a nonlinear elastic cell model composed of a linear elastic cell and a gap cell; a fluid mass cell model is appended.
4. The method for splicing nuclear reactor fuel assembly seismic analysis modeling substructure types according to claim 1, characterized in that: the substructure types include a structure for fixing enclosing plates which cannot generate transverse vibration and a structure for generating single assembly which can generate transverse vibration.
5. The method for splicing nuclear reactor fuel assembly seismic analysis modeling substructure types according to claim 3, characterized in that: the model base parameters of the beam element model, the spring element model, the nonlinear spring element model and the additional fluid mass element can be modified as required.
6. The method for splicing nuclear reactor fuel assembly seismic analysis modeling substructure types according to claim 3, characterized in that: the beam element model comprises a two-dimensional Euler Bernoulli finite element beam model, a finite element beam model considering a shearing effect and an axially deformed finite element beam model.
7. The method for splicing nuclear reactor fuel assembly seismic analysis modeling substructure types according to claim 3, characterized in that: the spring unit is a linear spring unit model constrained by a linear spring, and is characterized in that the linear spring unit is used for constraining the displacement of two adjacent nodes.
8. The method for splicing nuclear reactor fuel assembly seismic analysis modeling substructure types according to claim 3, characterized in that: the nonlinear spring unit is formed by coupling and combining a linear spring unit model and a gap unit model by means of displacement of end node, and the nonlinear spring unit model has the functions of the linear spring unit and the gap unit at the same time.
9. The method for splicing nuclear reactor fuel assembly seismic analysis modeling substructure types according to claim 3, characterized in that: the additional fluid mass cell model includes a fluid additional mass influence coefficient.
10. The method for splicing nuclear reactor fuel assembly seismic analysis modeling substructure types according to claim 1, characterized in that: the detailed steps in the step 5 are as follows: obtaining vibration equations of a plurality of single assemblies by adopting substructure modeling, wherein the vibration equations comprise the following steps:
Figure RE-FDA0002998196330000021
the equations are coupled by means of nonlinear force at the right end, and a coupling iterative solution format is constructed according to the force coupling and the estimation correction format.
CN202011279106.1A 2020-11-16 2020-11-16 Method for splicing nuclear reactor fuel assembly anti-seismic analysis modeling substructure types Pending CN112836269A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202011279106.1A CN112836269A (en) 2020-11-16 2020-11-16 Method for splicing nuclear reactor fuel assembly anti-seismic analysis modeling substructure types

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202011279106.1A CN112836269A (en) 2020-11-16 2020-11-16 Method for splicing nuclear reactor fuel assembly anti-seismic analysis modeling substructure types

Publications (1)

Publication Number Publication Date
CN112836269A true CN112836269A (en) 2021-05-25

Family

ID=75923113

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202011279106.1A Pending CN112836269A (en) 2020-11-16 2020-11-16 Method for splicing nuclear reactor fuel assembly anti-seismic analysis modeling substructure types

Country Status (1)

Country Link
CN (1) CN112836269A (en)

Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS6239789A (en) * 1985-08-14 1987-02-20 ウエスチングハウス・エレクトリック・コ−ポレ−ション Neutron reflector
CN101905340A (en) * 2010-07-23 2010-12-08 西安交通大学 Method for rapidly judging stability of high-speed milling
CN104636556A (en) * 2015-02-09 2015-05-20 武汉理工大学 Vibration response calculating method of limited baseboard structure connected at any angle
CN107044895A (en) * 2017-02-22 2017-08-15 中国水利水电科学研究院 Force test method is collided when spentnuclear fuel screen work is with component seismic test
CN107195341A (en) * 2017-05-17 2017-09-22 华北电力大学 The on-line monitoring device that nuclear power plant containment shell fragment influences on fuel assembly pressure drop
CN108899098A (en) * 2018-07-20 2018-11-27 中广核研究院有限公司 A kind of in-pile component structure suitable for presurized water reactor
CN109378096A (en) * 2018-11-12 2019-02-22 中国原子能科学研究院 Reactor core refueling method for deep pool type low-temperature heat supply reactor
CN110598324A (en) * 2019-09-12 2019-12-20 西安交通大学 Nuclear reactor dispersion plate type fuel element core fluid-solid coupling calculation method
CN110853784A (en) * 2019-11-19 2020-02-28 中国核动力研究设计院 Pressurized water reactor nuclear power plant stabiliser bearing structure

Patent Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS6239789A (en) * 1985-08-14 1987-02-20 ウエスチングハウス・エレクトリック・コ−ポレ−ション Neutron reflector
CN101905340A (en) * 2010-07-23 2010-12-08 西安交通大学 Method for rapidly judging stability of high-speed milling
CN104636556A (en) * 2015-02-09 2015-05-20 武汉理工大学 Vibration response calculating method of limited baseboard structure connected at any angle
CN107044895A (en) * 2017-02-22 2017-08-15 中国水利水电科学研究院 Force test method is collided when spentnuclear fuel screen work is with component seismic test
CN107195341A (en) * 2017-05-17 2017-09-22 华北电力大学 The on-line monitoring device that nuclear power plant containment shell fragment influences on fuel assembly pressure drop
CN108899098A (en) * 2018-07-20 2018-11-27 中广核研究院有限公司 A kind of in-pile component structure suitable for presurized water reactor
CN109378096A (en) * 2018-11-12 2019-02-22 中国原子能科学研究院 Reactor core refueling method for deep pool type low-temperature heat supply reactor
CN110598324A (en) * 2019-09-12 2019-12-20 西安交通大学 Nuclear reactor dispersion plate type fuel element core fluid-solid coupling calculation method
CN110853784A (en) * 2019-11-19 2020-02-28 中国核动力研究设计院 Pressurized water reactor nuclear power plant stabiliser bearing structure

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
齐欢欢等: "基于ANSYS的燃料组件事故动力分析程序", 《核动力工程》 *

Similar Documents

Publication Publication Date Title
CN111291494B (en) Multi-scale multi-physical field coupling simulation method for TRISO fuel particles of nuclear reactor
Lin et al. A high efficient assembly technique for large PEMFC stacks: Part I. Theory
CN105653869B (en) A kind of supercritical water reactor reactor core Three dimensional transient method for analyzing performance
CN110619145A (en) Self-adaptive modeling method for flexible support gear transmission device
CN112668072A (en) Modeling method for nuclear reactor fuel assembly seismic analysis
Subhan et al. Development of nuclear fuel assembly 3D reduced model and modal analysis
Leturcq et al. A new reduced order model to represent the creep induced fuel assembly bow in PWR cores
CN112836269A (en) Method for splicing nuclear reactor fuel assembly anti-seismic analysis modeling substructure types
CN112699433A (en) Method for classifying nuclear reactor fuel assembly seismic analysis mold building structure types
Lee et al. Design improvement of an OPT-H type nuclear fuel rod support grid by using an axiomatic design and an optimization
CN108038262A (en) A kind of floor response spectrum simplified calculation method of consideration SSSI effects
Aldebie et al. Neutronic and thermal analysis of heat pipe cooled graphite moderated micro reactor
Kępisty et al. SFR mechanical scenarios and neutron transport transients with CAST3M code
Sezen et al. Mechanistic and probabilistic seismic assessment of structures and components in nuclear power plants
Stimpson et al. Demonstration of Coupled Tiamat Single Assembly Calculations
Stimpson et al. Demonstration of Coupled Fuel Performance Calculations in VERA on Watts Bar Unit 1, Cycle 1
Hussain et al. CNPP fuel rod vibration analysis using finite element method
Choia et al. Modeling and Simulation Needs and Capabilities for Artificial Intelligence Based Plant Reload Optimization Platform
Anjum et al. Multiphysics analysis of mechanical responses in micro-reactors under varied operating conditions
Roberts et al. Accelerated Thermal Battery Design Through Digital Engineering Tools.
XiCheng et al. A Finite Element Model For NHR200-? 9x9 Fuel Assembly In Dynamic Analysis
CN115186603A (en) Method for calculating interaction between fuel rod and clamping mechanism
CN118312762A (en) Polycondensation method and system based on data driving
Brandhorst Jr New 5 Kilowatt Free‐Piston Stirling Space Converter Developments
Chaube et al. A Preliminary Fluid-Structure Coupling of NekRS and MOOSE via Cardinal

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
WD01 Invention patent application deemed withdrawn after publication

Application publication date: 20210525

WD01 Invention patent application deemed withdrawn after publication