Loss mass effect calculation method for nuclear island plant anti-seismic analysis and related components

Abstract

The invention relates to the technical field of earthquake resistance, and discloses a nuclear island plant earthquake resistance analysis loss mass effect calculation method and related components, and the method comprises the steps: calculating static effects Rx, Ry and Rz of a nuclear island plant structure under the condition that unit acceleration excitation is applied in x, y and z directions; performing modal analysis on the nuclear island plant structure to obtain all vibration modes with natural vibration frequencies smaller than a preset rigid reaction frequency, and extracting natural vibration circular frequencies omega i corresponding to the vibration modes, coefficients gamma ix, gamma iy and gamma iz in vibration mode participation coefficients, static reaction # imgabs0 # corresponding to vibration mode vectors {phi i} and other parameters; and respectively calculating corresponding lost mass effects when unit acceleration excitation is applied in the x, y and z directions according to the parameters. According to the method, on the basis of the principle of structural dynamics, the corresponding'lost mass' effect is calculated when unit acceleration excitation is applied in the x, y and z directions, and the method has the advantage of being accurate in calculation result.

Description

Technical Field

[0001] The embodiments of the present invention relate to the field of earthquake resistance technology, and in particular to a method for calculating lost mass effect of earthquake resistance analysis of a nuclear island plant and related components. Background Art

[0002] Earthquake action is one of the control conditions for the structural design of nuclear island buildings. The earthquake-resistant analysis of nuclear island buildings usually uses the vibration mode decomposition response spectrum method to calculate its earthquake effect (such as internal force, deformation, etc.); the modal superposition time history analysis method is used to calculate the earthquake floor response spectrum to provide input data for the earthquake-resistant analysis of equipment, pipelines, etc. The basis of these two analysis methods is the results of structural modal analysis.

[0003] Nuclear island buildings usually use reinforced concrete structures composed of thick walls and thick plates. Compared with the softer civil building structures, their natural vibration frequency is greater than the rigid response frequency f ZPA (usually 33Hz) is an important vibration mode, and there are many local vibrations. In actual calculations, it is impossible to cover all vibration modes of the nuclear island structure, which inevitably ignores the response corresponding to the high-order vibration modes. Even if the seismic analysis method of controlling the number of vibration modes of civil building structures is borrowed (that is, the cumulative vibration mode participation mass in the three orthogonal translation directions is required to reach more than 90% respectively), the local seismic response of the nuclear island plant structure will be underestimated, because the unconsidered high-order vibration modes may be the vibration of a local area. Therefore, the seismic analysis of the nuclear island structure must take into account the influence of "lost mass".

[0004] The calculation theory of the "lost mass" effect aims to obtain the difference terms of the equivalent seismic force and apply them to the corresponding nodes of the structural Ansys calculation model to ensure the accuracy of the seismic calculation of the nuclear island structure.

[0005] The existing technology uses the equivalent nodal force method specified in the design specification to calculate the "lost mass" effect under unit acceleration excitation in the x, y, and z directions. To calculate the equivalent nodal force, it is necessary to extract the moment of inertia Jx of each node in the Ansys calculation model. j , Jy j , Jz j , but it is difficult to obtain the exact value at present, and it is usually assumed to be 0, resulting in inaccurate calculation results of the equivalent nodal force. Summary of the invention

[0006] The embodiment of the present invention provides a method for calculating the lost mass effect of seismic analysis of a nuclear island plant and related components, aiming to solve the problem that the existing technology uses an equivalent nodal force method to calculate the "lost mass" effect with uncertain assumptions, resulting in an inaccurate calculation structure.

[0007] In a first aspect, an embodiment of the present invention provides a method for calculating the missing mass effect of a nuclear island plant in seismic analysis, which is characterized by comprising:

[0008] Calculate the static effect R of the nuclear island building structure under unit acceleration excitation in the x, y, and z directions x , R y and R z ;

[0009] Perform modal analysis on the nuclear island plant structure to obtain all vibration modes whose natural frequencies are less than the preset rigid response frequency, and extract the natural circular frequency ω corresponding to each vibration mode. i , the coefficient γ in the mode participation coefficient ix , γ iy , γ iz and vibration mode vector {φ i The corresponding static response

[0010] According to the static effect R x , R y and R z , and the natural circular frequency ω corresponding to each vibration mode i , the coefficient γ in the mode participation coefficient ix , γ iy , γ iz and vibration mode vector {φ i The corresponding static response The corresponding missing mass effect when unit acceleration excitation is applied in the x, y, and z directions is calculated respectively.

[0011] In a second aspect, an embodiment of the present invention provides a device for calculating lost mass effect of seismic analysis of nuclear island plant, characterized by comprising:

[0012] The static effect calculation unit is used to calculate the static effect R of the nuclear island plant structure under unit acceleration excitation in the x, y, and z directions. x , R y and R z ;

[0013] The parameter acquisition unit is used to perform modal analysis on the nuclear island plant structure to obtain all vibration modes whose natural vibration frequencies are less than the preset rigid response frequency, and to extract the natural vibration circular frequency ω corresponding to each vibration mode. i , the coefficient γ in the mode participation coefficient ix , γ iy , γ iz and vibration mode vector {φ i The corresponding static response

[0014] The lost mass effect calculation unit is used to calculate the static effect R x , R y and R z, and the natural circular frequency ω corresponding to each vibration mode i , the coefficient γ in the mode participation coefficient ix , γ iy , γ iz and vibration mode vector {φ i The corresponding static response The corresponding missing mass effect when unit acceleration excitation is applied in the x, y, and z directions is calculated respectively.

[0015] In a third aspect, an embodiment of the present invention provides a computer device, which includes a memory, a processor, and a computer program stored in the memory and executable on the processor, wherein when the processor executes the computer program, the method for calculating the lost mass effect of seismic analysis of a nuclear island building as described in the first aspect above is implemented.

[0016] In a fourth aspect, an embodiment of the present invention provides a computer-readable storage medium, wherein the computer-readable storage medium stores a computer program, and when the computer program is executed by a processor, the processor executes the method for calculating the lost mass effect of seismic analysis of a nuclear island building as described in the first aspect above.

[0017] The beneficial effects of the embodiments of the present invention are:

[0018] The present invention is based on the principle of structural dynamics and realizes the calculation of the "lost mass" effect corresponding to the application of unit acceleration excitation in the x, y and z directions, and has the advantage of accurate calculation results.

[0019] The present invention calculates the "lost mass" effect when unit acceleration excitation is applied in the x, y and z directions respectively, and on this basis, can conveniently calculate the total "lost mass" effect based on the vibration mode decomposition response spectrum method and the modal superposition time history analysis method commonly used in the seismic analysis of the nuclear island plant structure under any excitation size in each direction. BRIEF DESCRIPTION OF THE DRAWINGS

[0020] In order to more clearly illustrate the technical solutions of the embodiments of the present invention, the accompanying drawings required for use in the description of the embodiments will be briefly introduced below. Obviously, the accompanying drawings described below are some embodiments of the present invention. For ordinary technicians in this field, other accompanying drawings can be obtained based on these accompanying drawings without paying any creative work.

[0021] Figure 1 A schematic diagram of deformation of a multi-degree-of-freedom structural system under earthquake action provided by an embodiment of the present invention;

[0022] Figure 2 A schematic flow chart of a method for calculating the lost mass effect of a nuclear island plant building seismic analysis provided by an embodiment of the present invention;

[0023] Figure 3 Another schematic flow chart of a method for calculating the lost mass effect of a nuclear island plant building seismic analysis provided by an embodiment of the present invention;

[0024] Figure 4 A schematic block diagram of a device for calculating lost mass effect of seismic analysis of a nuclear island plant provided by an embodiment of the present invention;

[0025] Figure 5 A schematic block diagram of a computer device provided in an embodiment of the present invention. DETAILED DESCRIPTION

[0026] The following will be combined with the drawings in the embodiments of the present invention to clearly and completely describe the technical solutions in the embodiments of the present invention. Obviously, the described embodiments are part of the embodiments of the present invention, not all of the embodiments. Based on the embodiments of the present invention, all other embodiments obtained by ordinary technicians in this field without creative work are within the scope of protection of the present invention.

[0027] It should be understood that when used in this specification and the appended claims, the terms "include" and "comprises" indicate the presence of described features, integers, steps, operations, elements and / or components, but do not exclude the presence or addition of one or more other features, integers, steps, operations, elements, components and / or combinations thereof.

[0028] It should also be understood that the terms used in this specification of the present invention are only for the purpose of describing specific embodiments and are not intended to limit the present invention. As used in the specification of the present invention and the appended claims, unless the context clearly indicates otherwise, the singular forms "a", "an" and "the" are intended to include plural forms.

[0029] It should be further understood that the term "and / or" used in the present description and the appended claims refers to any and all possible combinations of one or more of the associated listed items, and includes these combinations.

[0030] As used in this specification and the appended claims, the term "if" may be interpreted as "when" or "upon" or "in response to determining" or "in response to detecting," depending on the context. Similarly, the phrases "if it is determined" or "if [described condition or event] is detected" may be interpreted as meaning "upon determination" or "in response to determining" or "upon detection of [described condition or event]" or "in response to detecting [described condition or event]," depending on the context.

[0031] To facilitate understanding of the invention of this application, the calculation principle of the existing "missing mass" effect is first introduced as follows:

[0032] Under the action of ground motion caused by earthquake, assuming that the ground motion is the same at all support points of the structure, the degree of freedom is the deformation of the structural system (n is the number of nodes in the structural system, each node has horizontal x, y and vertical z translational freedom and rotational freedom around the x, y, z axis) as shown in Figure 1 As shown, its differential equation of motion is as follows:

[0033]

[0034] In formula (1), [M], [C] and [K] are the mass matrix, damping matrix and stiffness matrix of the structural system of N×N order, respectively; the mass matrix [M] is a diagonal matrix, which is expressed in block matrix form as follows:

[0035]

[0036] In formula (1), and {u(t)} are the N×1-order relative acceleration (including translational and rotational acceleration), relative velocity (including translational and rotational velocity) and relative displacement (including translational and rotational displacement) reaction vectors of the structural system, which are functions of time t; corresponding to the mass matrix [M], the first 1~n, n+1~2n, 2n+1~3n elements correspond to the translational relative acceleration, translational relative velocity and translational relative displacement reactions of the structural system in the x, y, and z directions, respectively; the 3n+1~4n, 4n+1~5n, and 5n+1~6n elements correspond to the relative angular acceleration, relative angular velocity and relative angular displacement reactions of the structural system around the x, y, and z axes, respectively.

[0037] In formula (1), {R} is the N×1-order earthquake excitation indicator vector of the structural system. Corresponding to the mass matrix, it can be expressed as {R}={{D x} T ,{D y} T ,{D z} T , {0} T ,{0} T ,{0} T}, where {D x} T ,{D y} T ,{D z} T are all N-order column vectors. For the assumption that the ground motion is the same at all support points of the structural system, the elements are the same, that is, D jx =D x j∈(1,n), D jy =D y j∈(1,n), Djz =D z j∈(1,n); when there is earthquake excitation in only one direction (assuming x direction), D x =1, the rest D y =D z =0,{0} T is an n-order column vector, that is, there is no seismic excitation in the direction of rotation around the x, y, and z axes

[0038] In formula (1), is the time history of earthquake ground motion acceleration.

[0039] According to the theory of linear algebra, an N-dimensional vector can always be expressed as a linear combination of N independent vectors. The seismic relative displacement response vector {u(t)} of the structural system can be expressed as a vibration mode vector {φ i}, (i=1,2,……,N) linear combination:

[0040]

[0041] In formula (2), the vibration mode vector {φ i} is the displacement offset vector of the undamped structural system corresponding to a certain natural frequency. It is the ratio of displacement, not the actual displacement. By setting a suitable ratio, the following formula M can be obtained: i * ={φ i} T [M]{φ i}=1(i=j);{φ i The element φ in niUx ,φ niUy ,φ niUz are the translational displacement coefficients of the nth node in the x, y, and z directions under the ith vibration mode; {φ i} satisfies the following formula (3):

[0042] [K]{φ i}=ω i 2 [M]{φ i}(3);

[0043] In formula (3), ω i is the natural circular frequency of the i-th vibration mode of the structural system; q i (t)(i=1,2,……,N) is the vibration mode canonical coordinate, which is a function of time t.

[0044] Substituting formula (2) into formula (1), we can obtain:

[0045]

[0046] Multiply both sides of equation (4) by {φi} T Then divide by {φ i} T [M]{φ i}, and use the orthogonality of the vibration mode with respect to the mass matrix [M] and the stiffness matrix [K], and assume that the vibration mode is also orthogonal with respect to the damping matrix [C], that is:

[0047]

[0048] And order

[0049]

[0050] We can get:

[0051]

[0052] It can be seen that the original N-dimensional simultaneous differential motion equations reflecting the motion characteristics of the structural system with N degrees of freedom are decomposed into N independent canonical coordinates q i The differential equation of motion of the single-degree-of-freedom system; the natural circular frequency ω of each single-degree-of-freedom system i and damping ratio ζ i is the circular frequency and damping ratio corresponding to each vibration mode of the multi-degree-of-freedom system, γ i is the participation coefficient of the i-th mode of the multi-freedom system.

[0053] By using Duhamel integral, the solution of equation (6) is:

[0054]

[0055] In formula (7), Δ i (t) is the self-oscillating circular frequency ω i , the damping ratio is ζ i A single degree of freedom system under earthquake excitation The seismic displacement response under the condition of ; therefore, equation (2) can be converted into:

[0056]

[0057] The seismic effect of the structural system at each moment can be determined by the structural static analysis [K]{u(t)}, which is the [K]{u i (t)}=[K]γ i {φ i}Δ i (t), the algebraic sum of i∈(1,N).

[0058] The vector {R} can also be expressed as N independent vectors {φ i}, that is:

[0059]

[0060] Multiply both sides of equation (9) by {φ i} T [M] then divided by {φ i} T [M]{φ i}, and using the orthogonality of vibration modes, we get:

[0061]

[0062] Visible expansion coefficient Γ i and mode participation coefficient γ i The same, that is It shows that the sum of the products of each order mode vector and the corresponding mode participation coefficient is a constant vector, namely the earthquake excitation indicator vector {R}; therefore, equation (4) can be converted into:

[0063]

[0064] From equation (11) and the orthogonality of the above vibration modes with respect to the mass matrix [M], the damping matrix [C] and the stiffness matrix [K], it can be seen that the force vector Only the corresponding i-th vibration mode response is generated and the i-th vibration mode response is completely determined by the force vector When only the first m (m<<N) order vibration modes are calculated, the difference term of the equivalent seismic force shown on the right side of equation (11) is

[0065]

[0066] Formula (12) is the theoretical basis for the calculation method of the "lost mass" effect in the seismic analysis of the existing nuclear island plant. At the same time, it is required that the natural frequencies of the uncalculated modes, i.e., the m+1 to N-th order modes, are all greater than the rigid response frequency f ZPA (usually 33Hz), the responses corresponding to these vibration modes only contain the rigid response components of the structure.

[0067] Based on this, the existing technology can calculate the difference term of the equivalent seismic force based on formula (12) and apply it to the corresponding nodes of the structure:

[0068] When there is only earthquake excitation in the x direction When (Dx=1, Dy=Dz=0), the equivalent nodal force of the unit acceleration seismic vector error shown in equation (12) is:

[0069]

[0070] The corresponding "missing mass" effect is the effect under the above equivalent nodal force multiplied by the corresponding seismic excitation

[0071] When there is only seismic excitation in the y direction, (Dy=1, Dx=Dz=0), the equivalent nodal force of the unit acceleration seismic vector error shown in equation (12) is:

[0072]

[0073] The corresponding "missing mass" effect is the effect under the above equivalent nodal force multiplied by the corresponding seismic excitation

[0074] When there is only earthquake excitation in the z direction When (Dz=1, Dx=Dy=0), the equivalent nodal force of the unit acceleration seismic vector error shown in equation (12) is:

[0075]

[0076] The corresponding "missing mass" effect is the effect under the above equivalent nodal force multiplied by the corresponding seismic excitation

[0077] To calculate the equivalent nodal forces, it is necessary to extract the structural mass matrix [M] formed by the calculation software to obtain the value mx of any element in the mass matrix [M] j 、my j 、mz j , Jx j , Jy j , Jz j , to calculate the forces and moments applied to the nodes. However, the structural mass matrix [M] formed by the calculation software is related to the calculation assumptions of the software, and it is often difficult to obtain the data directly and accurately. When the software calculates, the degrees of freedom of all nodes of the structure are usually constrained, and unit acceleration is applied in the three translational x, y, and z directions respectively, and the constraint reaction force of each constraint node in the corresponding direction is extracted as the mass size in that direction. That is, when a unit acceleration is applied in the x direction, the constraint reaction force of each constraint node in the x direction is extracted, and the mass size mx of each node in the x direction can be obtained. j , similarly we can find my j 、mz j . However, for the moment of inertia Jx of each node j , Jy j , Jz j Difficult to obtain, usually assumed to be 0.

[0078] The above is the calculation of the "lost mass" effect under unit acceleration excitation in the x, y, and z directions using the equivalent nodal force method specified in the existing design specifications. To calculate the equivalent nodal force, it is necessary to extract the moment of inertia Jx of each node in the Ansys calculation model. j , Jy j , Jzj , but it is difficult to obtain the exact value at present, and it is usually assumed to be 0, resulting in inaccurate calculation results of the equivalent nodal force.

[0079] Based on this, this application develops a calculation method for the "lost mass" effect of nuclear island plant seismic analysis based on the principle of structural dynamics to solve the above-mentioned node moment of inertia Jx j , Jy j , Jz j It is difficult to obtain and is usually assumed to be 0, which leads to inaccurate calculation results.

[0080] The principle of calculating the "lost mass" effect of nuclear island buildings in the present invention is briefly described as follows:

[0081] when When the nuclear island building structure is subjected to unit acceleration excitation static load, its response relative to time t is a constant, that is, {u(t)}={u}; that is, the differential equation of motion (1) can be "degenerated" to

[0082]

[0083] In this way, the response of the nuclear island building structure can be regarded as and The sum of these two parts of excitation; when only the first m (m<<N) order vibration modes are calculated, the structure is The following reaction is the "lost mass" effect, namely:

[0084]

[0085] The nuclear island plant structure is The response to the stimulus can be seen as The sum of the reactions under i The reaction under} is r i mode , obviously The reaction coefficient is With r i mode The product of The response to the stimulus is

[0086] According to the linear superposition principle, the "missing mass" effect of the nuclear island structure can be obtained by subtracting its effect under -[M]{R} from its effect under The effect of the acceleration is obtained; therefore, the response under the unit acceleration excitation in the x direction can be obtained. The response under the unit acceleration excitation in the x direction is recorded as R x The response under unit acceleration excitation in the y direction is recorded as R yand the response under unit acceleration excitation in the z direction only is denoted as R z .

[0087] When there is only unit acceleration earthquake excitation in the x direction (Dx=1, Dy=Dz=0), the corresponding "missing mass" effect is

[0088] When there is only unit acceleration earthquake excitation in the y direction (Dz = 1, Dx = Dy = Rx = Ry = Rz = 0), the corresponding "missing mass" effect is

[0089] When there is only unit acceleration earthquake excitation in the z direction (Dz = 1, Dx = Dy = Rx = Ry = Rz = 0), the corresponding "missing mass" effect is

[0090] Based on the brief description of the principle of calculating the "missing mass" effect of the nuclear island plant in the present invention, the process of calculating the "missing mass" effect of the nuclear island plant using the Ansys platform is as follows:

[0091] See also Figure 2 , Figure 2 A schematic flow chart of a method for calculating the lost mass effect of a nuclear island plant building seismic analysis provided by an embodiment of the present invention;

[0092] like Figure 2 As shown, the method includes steps S201 to S103.

[0093] S201. Calculate the static effect R of the nuclear island building structure under unit acceleration excitation in the x, y and z directions respectively. x , R y and R z .

[0094] S202, perform modal analysis on the nuclear island plant structure to obtain all vibration modes whose natural vibration frequencies are less than the preset rigid response frequency, and extract the natural vibration circular frequency ω corresponding to each vibration mode i , the coefficient γ in the mode participation coefficient ix , γ iy , γ iz and vibration mode vector {φ i The corresponding static response

[0095] S203, according to the static effect R x , R y and R z , and the natural circular frequency ω corresponding to each vibration mode i , the coefficient γ in the mode participation coefficient ix , γ iy , γ izand vibration mode vector {φ i The corresponding static response The corresponding missing mass effect when unit acceleration excitation is applied in the x, y, and z directions is calculated respectively.

[0096] In this embodiment, the finite element software Ansys can be used to perform the static effect R of step S101. x , R y and R z , and the extraction of each parameter in step S202.

[0097] In this embodiment, based on the process of steps S201-S203, it is not necessary to obtain the rotational inertia Jx of each node in the calculation model of the software Ansys. j , Jy j , Jz j , there is no need to assume the moment of inertia to be 0. This embodiment is based on the principle of structural dynamics and realizes the calculation of the corresponding missing mass effect when a unit acceleration excitation is applied in the x, y, and z directions, and has the advantage of accurate calculation results.

[0098] In one embodiment, step S203 includes:

[0099] When there is a unit acceleration earthquake excitation in the x direction, the corresponding missing mass effect is calculated as follows:

[0100]

[0101] When there is a unit acceleration earthquake excitation in the y direction, the corresponding missing mass effect is calculated according to the following formula:

[0102]

[0103] When there is a unit acceleration earthquake excitation in the z direction, the corresponding missing mass effect is calculated as follows:

[0104]

[0105] In this embodiment, the parameter data obtained in steps S201-S202 are substituted into the above corresponding formulas for calculation, and the corresponding lost mass effect when unit acceleration excitation is applied in the x, y, and z directions can be obtained.

[0106] In one embodiment, see Figure 3 The method also includes steps S301 to S302.

[0107] S301, based on the calculation of the missing mass effect corresponding to the unit acceleration excitation applied in the x, y, and z directions, applying the seismic excitation in the x, y, and z directions to obtain the total missing mass effect based on the modal superposition time history analysis method;

[0108] S302: Calculate the total missing mass effect based on the mode decomposition response spectrum method based on the absolute maximum value of the time history of the seismic excitation in the x, y, and z directions.

[0109] In this embodiment, steps S201-S203 are used to calculate the "lost mass" effect when unit acceleration excitation is applied in the x, y, and z directions respectively. Based on this, the total "lost mass" effect based on the modal superposition time-history analysis method and the vibration mode decomposition response spectrum method commonly used in the seismic analysis of the nuclear island plant structure under any excitation size in each direction can be conveniently calculated.

[0110] Specifically, in step S301, the total lost mass effect R based on the modal superposition time history analysis method is calculated according to the following formula: missingmass :

[0111]

[0112] in, represents the earthquake excitation in the x direction, represents the seismic excitation in the y direction, represents the earthquake excitation in the z direction, R x-missingmass represents the missing mass effect corresponding to the unit acceleration earthquake excitation in the x direction, R y-missingmass represents the missing mass effect corresponding to the unit acceleration earthquake excitation in the y direction, R z-missingmass It represents the missing mass effect corresponding to the unit acceleration seismic excitation in the z direction.

[0113] Specifically, in step S302, the total lost mass effect R based on the vibration mode decomposition response spectrum method is calculated according to the following formula: missingmass :

[0114] R missingmass =R x-missingmass ×ZPA x +R y-missingmass ×ZPA y +R z-missingmass ×ZPA z ;

[0115] in, ZPAx represents the absolute maximum value of the time history of the earthquake excitation in the x direction; ZPAy represents the absolute maximum value of the time history of the earthquake excitation in the y direction; ZPAz represents the absolute maximum value of the time history of the earthquake excitation in the z direction.

[0116] The embodiment of the present invention also provides a device for calculating the missing mass effect of a nuclear island plant in seismic analysis, which is used to execute any embodiment of the method for calculating the missing mass effect of a nuclear island plant in seismic analysis. Figure 4 , Figure 4 It is a schematic block diagram of a device for calculating lost mass effect of seismic analysis of a nuclear island plant provided by an embodiment of the present invention.

[0117] like Figure 4 As shown, the missing mass effect calculation device 400 for seismic analysis of nuclear island buildings includes: a static effect calculation unit 401, a parameter acquisition unit 402 and a missing mass effect calculation unit 403.

[0118] The static effect calculation unit 401 is used to calculate the static effect R of the nuclear island building structure under unit acceleration excitation in the x, y and z directions respectively. x , R y and R z ;

[0119] The parameter acquisition unit 402 is used to perform modal analysis on the nuclear island plant structure to obtain all vibration modes whose natural vibration frequencies are less than the preset rigid response frequency, and extract the natural vibration circular frequency ω corresponding to each vibration mode. i , the coefficient γ in the mode participation coefficient ix , γ iy , γ iz The static response corresponding to the vibration mode vector;

[0120] The lost mass effect calculation unit 403 is used to calculate the static effect R x , R y and R z , and the natural circular frequency ω corresponding to each vibration mode i , the coefficient γ in the mode participation coefficient ix , γ iy , γ iz The static response corresponding to the vibration mode vector is calculated by respectively calculating the missing mass effect when unit acceleration excitation is applied in the x, y, and z directions.

[0121] The device can use the finite element software Ansys to calculate the static effects Rx, Ry and Rz in step S101 and extract the parameters in step S202.

[0122] In this embodiment, based on the process of steps S201-S203, it is not necessary to obtain the rotational inertia Jx of each node in the calculation model of the software Ansys. j , Jy j , Jz j, there is no need to assume the moment of inertia to be 0. This embodiment is based on the principle of structural dynamics and realizes the calculation of the corresponding missing mass effect when a unit acceleration excitation is applied in the x, y, and z directions, and has the advantage of accurate calculation results.

[0123] Those skilled in the art can clearly understand that, for the convenience and brevity of description, the specific working processes of the above-described devices and units can refer to the corresponding processes in the aforementioned method embodiments and will not be repeated here.

[0124] The above-mentioned calculation device for lost mass effect of nuclear island building seismic analysis can be realized in the form of a computer program. The computer program can be used in Figure 5 Runs on the computer device shown.

[0125] See also Figure 5 , Figure 5 5 is a schematic block diagram of a computer device provided by an embodiment of the present invention. The computer device 500 is a server, which can be an independent server or a server cluster composed of multiple servers.

[0126] See also Figure 5 The computer device 500 includes a processor 502 , a memory and a network interface 505 connected via a system bus 501 , wherein the memory may include a non-volatile storage medium 503 and an internal memory 504 .

[0127] The non-volatile storage medium 503 can store an operating system 5031 and a computer program 5032. When the computer program 5032 is executed, the processor 502 can execute a method for calculating the missing mass effect of nuclear island plant seismic analysis.

[0128] The processor 502 is used to provide computing and control capabilities to support the operation of the entire computer device 500 .

[0129] The internal memory 504 provides an environment for the operation of the computer program 5032 in the non-volatile storage medium 503. When the computer program 5032 is executed by the processor 502, the processor 502 can execute the lost mass effect calculation method for the seismic analysis of the nuclear island plant.

[0130] The network interface 505 is used for network communication, such as providing data information transmission, etc. Those skilled in the art will understand that Figure 5 The structure shown in the figure is only a block diagram of a part of the structure related to the solution of the present invention, and does not constitute a limitation on the computer device 500 to which the solution of the present invention is applied. The specific computer device 500 may include more or fewer components than those shown in the figure, or combine certain components, or have a different arrangement of components.

[0131] Those skilled in the art will understand that Figure 5 The embodiments of the computer device shown in the figure do not constitute a limitation on the specific composition of the computer device. In other embodiments, the computer device may include more or fewer components than shown in the figure, or combine certain components, or arrange the components differently. For example, in some embodiments, the computer device may only include a memory and a processor. In such embodiments, the structure and function of the memory and the processor are the same as those of the embodiment of the present invention. Figure 5 The embodiments shown are consistent and will not be described again here.

[0132] It should be understood that in the embodiment of the present invention, the processor 502 may be a central processing unit (CPU), and the processor 502 may also be other general-purpose processors, digital signal processors (DSP), application-specific integrated circuits (ASIC), field-programmable gate arrays (FPGA) or other programmable logic devices, discrete gate or transistor logic devices, discrete hardware components, etc. Among them, the general-purpose processor may be a microprocessor or the processor may also be any conventional processor, etc.

[0133] It is understood by those skilled in the art that all or part of the processes in the method for implementing the above embodiment can be completed by instructing the relevant hardware through a computer program. The computer program can be stored in a storage medium, which is a computer-readable storage medium. The computer program is executed by at least one processor in the computer system to implement the process steps of the embodiment of the above method.

[0134] In another embodiment of the present invention, a computer-readable storage medium is provided. The computer-readable storage medium may be a non-volatile computer-readable storage medium. The computer-readable storage medium stores a computer program, which, when executed by a processor, causes the processor to execute any embodiment of the method for calculating the lost mass effect of the seismic analysis of a nuclear island plant.

[0135] The storage medium is a physical, non-transient storage medium, for example, it can be a USB flash drive, a mobile hard disk, a read-only memory (ROM), a magnetic disk or an optical disk, etc., which can store program codes.

[0136] Those of ordinary skill in the art will appreciate that the units and algorithm steps of each example described in conjunction with the embodiments disclosed herein can be implemented in electronic hardware, computer software, or a combination of the two. In order to clearly illustrate the interchangeability of hardware and software, the composition and steps of each example have been generally described in terms of function in the above description. Whether these functions are performed in hardware or software depends on the specific application and design constraints of the technical solution. Professional and technical personnel can use different methods to implement the described functions for each specific application, but such implementation should not be considered to be beyond the scope of the present invention.

[0137] In the several embodiments provided by the present invention, it should be understood that the disclosed devices and methods can be implemented in other ways. For example, the device embodiments described above are only schematic. For example, the division of each unit is only a logical function division, and there may be other division methods in actual implementation. For example, multiple units or components can be combined or integrated into another system, or some features can be ignored or not executed.

[0138] The steps in the method of the embodiment of the present invention can be adjusted in order, combined and deleted according to actual needs. The units in the device of the embodiment of the present invention can be combined, divided and deleted according to actual needs. In addition, the functional units in the various embodiments of the present invention can be integrated into one processing unit, or each unit can exist physically separately, or two or more units can be integrated into one unit.

[0139] If the integrated unit is implemented in the form of a software functional unit and sold or used as an independent product, it can be stored in a storage medium. Based on this understanding, the technical solution of the present invention is essentially or the part that contributes to the prior art, or all or part of the technical solution can be embodied in the form of a software product, and the computer software product is stored in a storage medium, including a number of instructions for an air conditioner to perform all or part of the steps of the method described in each embodiment of the present invention.

[0140] In the above embodiments, the description of each embodiment has its own emphasis. For parts that are not described in detail in a certain embodiment, reference can be made to the relevant descriptions of other embodiments.

[0141] Obviously, those skilled in the art can make various changes and modifications to the present invention without departing from the spirit and scope of the present invention. Thus, if these modifications and variations of the present invention fall within the scope of the claims of the present invention and their equivalents, the present invention is also intended to include these modifications and variations.

[0142] The above is only a specific embodiment of the present invention, but the protection scope of the present invention is not limited thereto. Any technician familiar with the technical field can easily think of various equivalent modifications or replacements within the technical scope disclosed by the present invention, and these modifications or replacements should be included in the protection scope of the present invention. Therefore, the protection scope of the present invention shall be based on the protection scope of the claims.

Claims

1. A method for calculating the lost mass effect of seismic analysis of nuclear island buildings, characterized in that: include: Calculate the static effect R of the nuclear island building structure under unit acceleration excitation in the x, y, and z directions x , R y and R z ; Perform modal analysis on the nuclear island plant structure to obtain all vibration modes whose natural frequencies are less than the preset rigid response frequency, and extract the natural circular frequency ω corresponding to each vibration mode. i , the coefficient γ in the mode participation coefficient ix , γ iy , γ iz and vibration mode vector {φ i The corresponding static reaction r i mode ; According to the static effect R x , R y and R z , and the natural circular frequency ω corresponding to each vibration mode i , the coefficient γ in the mode participation coefficient ix , γ iy , γ iz and vibration mode vector {φ i The corresponding static reaction r i mode , calculate the corresponding lost mass effect when unit acceleration excitation is applied in the x, y, and z directions respectively.

2. The method for calculating the lost mass effect of seismic analysis of nuclear island buildings according to claim 1 is characterized in that: The step of calculating the corresponding lost mass effect when a unit acceleration excitation is applied in the x, y, and z directions comprises: When there is a unit acceleration earthquake excitation in the x direction, the corresponding missing mass effect is calculated as follows:

3. The method for calculating the lost mass effect of seismic analysis of nuclear island buildings according to claim 1 is characterized in that: The calculation of the corresponding lost mass effect when a unit acceleration excitation is applied in the x, y, and z directions includes: When there is a unit acceleration earthquake excitation in the y direction, the corresponding missing mass effect is calculated according to the following formula:

4. The method for calculating the lost mass effect of seismic analysis of nuclear island buildings according to claim 1 is characterized in that: The calculation of the corresponding lost mass effect when a unit acceleration excitation is applied in the x, y, and z directions includes: When there is a unit acceleration earthquake excitation in the z direction, the corresponding missing mass effect is calculated as follows:

5. The method for calculating the lost mass effect of seismic analysis of nuclear island buildings according to claim 1, characterized in that: Also includes: On the basis of calculating the missing mass effect corresponding to the unit acceleration excitation applied in the x, y and z directions, the seismic excitation in the x, y and z directions is applied to obtain the total missing mass effect based on the modal superposition time history analysis method; Based on the absolute maximum values ​​of the seismic excitations in the x, y and z directions during the time history, the total missing mass effect based on the mode decomposition response spectrum method is calculated.

6. The method for calculating the lost mass effect of seismic analysis of nuclear island buildings according to claim 5 is characterized in that: The method of applying seismic excitations in the x, y, and z directions based on the calculation of the corresponding missing mass effects when unit acceleration excitations are applied in the x, y, and z directions to obtain the total missing mass effect based on the modal superposition time history analysis method includes: The total missing mass effect R based on the modal superposition time history analysis method is calculated as follows: missingmass : in, represents the earthquake excitation in the x direction, represents the seismic excitation in the y direction, represents the earthquake excitation in the z direction, R x-missingmass represents the missing mass effect corresponding to the unit acceleration earthquake excitation in the x direction, R y-missingmass represents the missing mass effect corresponding to the unit acceleration earthquake excitation in the y direction, R z-missingmass It represents the missing mass effect corresponding to the unit acceleration seismic excitation in the z direction.

7. The method for calculating the lost mass effect of seismic analysis of nuclear island buildings according to claim 5 is characterized in that: The absolute maximum value in the time history process of the seismic excitation in the x, y, and z directions is used to calculate the total lost mass effect based on the vibration mode decomposition response spectrum method, including: The total missing mass effect R based on the modal decomposition response spectrum method is calculated as follows missingmass : R missingmass =R x-missingmass ×ZPA x +R y-missingmass ×ZPA y +R z-missingmass ×ZPA z ; in, ZPAx represents the absolute maximum value of the time history of the earthquake excitation in the x direction; ZPAy represents the absolute maximum value of the time history of the earthquake excitation in the y direction; ZPAz represents the absolute maximum value of the time history of the earthquake excitation in the z direction; R x-missingmass represents the missing mass effect corresponding to the unit acceleration earthquake excitation in the x direction, R y-missingmass represents the missing mass effect corresponding to the unit acceleration earthquake excitation in the y direction, R z-missingmass It represents the missing mass effect corresponding to the unit acceleration seismic excitation in the z direction.

8. A device for calculating the lost mass effect of seismic analysis of nuclear island buildings, characterized in that: include: The static effect calculation unit is used to calculate the static effect R of the nuclear island plant structure under unit acceleration excitation in the x, y, and z directions. x , R y and R z ; The parameter acquisition unit is used to perform modal analysis on the nuclear island plant structure to obtain all vibration modes whose natural vibration frequencies are less than the preset rigid response frequency, and to extract the natural vibration circular frequency ω corresponding to each vibration mode. i , the coefficient γ in the mode participation coefficient ix , γ iy , γ iz and vibration mode vector {φ i The corresponding static response The lost mass effect calculation unit is used to calculate the static effect R x , R y and R z , and the natural circular frequency ω corresponding to each vibration mode i , the coefficient γ in the mode participation coefficient ix , γ iy , γ iz and vibration mode vector {φ i The corresponding static response The corresponding missing mass effect when unit acceleration excitation is applied in the x, y, and z directions is calculated respectively.

9. A computer device comprising a memory, a processor, and a computer program stored in the memory and executable on the processor, characterized in that: When the processor executes the computer program, the method for calculating the lost mass effect of seismic analysis of a nuclear island plant according to any one of claims 1 to 7 is implemented.

10. A computer-readable storage medium, characterized in that: The storage medium stores a computer program, which, when executed by a processor, can implement the method for calculating lost mass effect of seismic analysis of a nuclear island plant according to any one of claims 1 to 7.

Images