CN104392067A - Method for calculating turbulence buffeting responses of heat exchanger - Google Patents

Abstract

The invention discloses a method for calculating turbulence buffeting responses of a heat exchanger. The method comprises steps as follows: a flow field of a secondary side of the heat exchanger is calculated by thermotechnical hydraulic system analysis software to obtain parameters of flow fields, distributed along a heat exchange pipe, of flow on the secondary side; on the basis of the parameters of the flow fields and the structure parameters of the heat exchange pipe, whether the parameters of the flow fields belong to a first preset range and a second preset range is judged; if the parameters of the flow fields belong to the first preset range and the second preset range, a reference equivalent power spectral density function of exciting force is calculated by means of formula (1)-(6); a root-mean-square value of the turbulence buffeting vibration responses is calculated by means of the random vibration theory based on a density function, flow-induced vibration evaluation is performed by means of the root-mean-square value, and the technical effects of more reasonable input and more accurate analysis result during calculation analysis of turbulence buffeting responses of the heat exchanger are achieved.

Description

A kind of computing method of heat interchanger turbulent buffeting response

Technical field

The present invention relates to reactor structure mechanics field, particularly relate to the computing method of a kind of heat interchanger turbulent buffeting response.

Background technology

A large amount of heat-exchanger rigs is had in nuclear energy, thermoelectricity, petrochemical complex equipment.The Flow vibration problem of these heat-exchanger rigs in operational process annoyings engineering circles always.Up to the present, academia compares the Flow vibration mechanism of approval is " vortex shedding ", " turbulent buffeting ", " fluid-elastic instability ", " acoustic resonance " etc.Wherein, the effect due to turbulent buffeting is long-term and cannot avoids, and needs to calculate turbulent buffeting response, for wearing and tearing and analysis of fatigue, to predict the security of heat interchanger long-time running in the design phase of the heat interchanger such as nuclear boiler.

The turbulent buffeting phenomenon of heat-transfer pipe tube bank refers to: produce turbulent flow when heat interchanger shell fluid flows through tube bank, and in turbulent flow, the pressure of pulsatile change and velocity field constantly supply pipe energy, and pipe absorbs energy and also produces vibration.The frequency range of turbulent fluctuation is wider and have very strong randomness, and turbulent buffeting RESPONSE CALCULATION needs to use random seismic response analysis method.

When the Flow vibration of heat interchanger is analyzed, the tube bundle vibration amplitude caused due to turbulent flow is little, usually turbulent buffeting is considered as the solid coupled problem of weak current, namely first determines the power spectrum density (PSD) that fluid encourages, then calculate the PSD response of tube bank.Here the fluid force that fluid matasomatism is structurally irrelevant with structure motion is called exciting force.The key of turbulent buffeting analysis is the PSD determining exciting force.If the PSD of exciting force can be obtained, just obtain the excitation input parameter that tube bundle vibration is analyzed, then in conjunction with certain assumed condition, just can according to random vibration theory, calculate bundle vibration response, thus the wearing and tearing can evaluated further because turbulent buffeting causes and fatigue damage.

In single-phase flow, the experimental study of turbulent exciting force is carried out comparatively abundant.By a large amount of experiments, the exciting force in different fluid media (medium)s and tube bank arrangement situation can be obtained.These experimental datas are compared by desirable dimensionless method, just very naturally can obtain the envelope spectrum of exciting force.Dimensionless method acknowledged at present reduces exciting force with dynamic head, reduces frequency, namely with the ratio of pitch flow velocity and caliber

$\begin{array}{c}{\stackrel{\‾}{\mathrm{\Φ}}}_E^0(f/f_0)=\frac{f_0}{{\left(p_{0}D\right)}^2}{\mathrm{\Φ}}_E^0\left(f\right)\\ p_0=\frac{1}{2}\mathrm{\ρ}{U_p}^2\\ f_0=U_p/D\end{array}---\left(7\right)$

In formula, for the equivalence of exciting force on unit pipe range is with reference to PSD, and for dimensionless number.The conversion of the experimental result of exciting force under a large amount of various flows field condition can be the dimensionless PSD that can conveniently compare by through type (7), thus reasonably can determine the enveloping curve of exciting force PSD.

For two-phase flow, similar research work does not have single-phase flow smooth like that.Void fraction lower than 10% time, random force and the single-phase flow of two-phase flow are similar, can use the envelope spectrum of single-phase flow.When void fraction is higher, the impact of two-phase mixture starts to play a leading role.Because two-phase flow has the different of a lot of essence compared with single-phase flow, the method for normalizing of single-phase flow exciting force and PSD envelope spectrum are not also suitable for two-phase flow, carry out Calculation of Vibration Response do not guard yet with the PSD envelope spectrum of single-phase flow exciting force to the tube bank under two-phase flow effect.

In in the past 30 years, the experimental result in a large number about two-phase flow exciting force is published at crowd, and particularly after the nineties in last century, this aspect research has a large amount of growths.Many researchers also propose method for normalizing for some experimental data.Although these researchs obtain much useful conclusion, but still leave over many problems and do not solved, the such as impact on exciting force such as viscosity, surface tension, density ratio and flow pattern.Although the problems referred to above do not solve, in order to the devisers solving heat interchanger are badly in need of two-phase flow exciting force as the requirement of engineering calculating input, some scholars have studied the method for normalizing of the two-phase random excitation power of Tube Sheet of Heat Exchanger.So far, comparatively comprehensively consider various experimental data obtain two-phase flow random excitation power PSD envelope spectrum for de Langre & Villard, they are on the basis of numerous scholar's experimental data, define the dimensionless formula of two-phase flow exciting force, and give the envelope spectrum of two-phase flow exciting force PSD on this basis.The contracting that they provide than parameter is:

f ₀＝V _p/D _w

(8)

p ₀＝ρ _lgD _w

In formula, V _pfor the pitch flow velocity of two-phase mixture, ρ _lfor the density of liquid phase, g is acceleration of gravity, length dimension D _wbe defined as:

$D_w=0.1D/\sqrt{1-{\mathrm{\α}}_H}---\left(9\right)$

D is pipe diameter, α _hfor homogeneous flow void fraction.

De Langre and Villard considers the contracting of a lot of pressure than parameter, finds that the contracting only based on gravity is the most effective than parameter.Author points out to contract and represents the impact of dynamic pressure than the gravity in parameter, and this dynamic pressure result in the generation of drift velocity between gas-liquid two-phase.Length dimension is the simplification form of distortion of the cavity characteristic length that Taylor etc. introduces.Utilize these two to contract than parameter and the definition with reference to equivalent PSD, the two-phase flow experimental data of multidigit scholar is shown in Fig. 2 by de Langre etc. simultaneously.Corresponding envelope spectrum is:

$\begin{array}{c}{\stackrel{\‾}{\mathrm{\Φ}}}_E^0\left(f_R\right)=10{f_R}^{-0.5},0.001\≤f_R\≤0.06\\ {\stackrel{\‾}{\mathrm{\Φ}}}_E^0\left(f_R\right)=2\×{10}^{-3}{f_R}^{-3.5},0.06\≤f_R\≤1\end{array}---\left(10\right)$

The fundamental purpose of de Langre and Villard is to solve requirement of engineering, but owing to have ignored some influence factors, it is too conservative that they provide exciting force PSD envelope spectrum, limits its application in engineering.Be compelled against one's will, in current engineering, the analysis of steam generator Flow vibration still mainly uses the PSD envelope spectrum based on single-phase flow.

In sum, present inventor, in the process realizing invention technical scheme in the embodiment of the present application, finds that above-mentioned technology at least exists following technical matters:

In the computing method of existing heat interchanger turbulent buffeting response, due to the exciting force method for normalizing based on single-phase flow PSD envelope spectrum and be not suitable for two-phase flow; And only PSD envelope spectrum based on two-phase flow is owing to have ignored some influence factors, provide exciting force PSD envelope spectrum too conservative.So it is unreasonable too high with conservative property to there is analysis input in existing heat interchanger Flow vibration analytical approach, causes the inaccurate technical matters of analysis result.

Summary of the invention

The invention provides the computing method of a kind of heat interchanger turbulent buffeting response, solve existing heat interchanger Flow vibration analytical approach existence analysis input unreasonable too high with conservative property, cause the inaccurate technical matters of analysis result, obtain the PSD envelope spectrum of new two-phase flow exciting force, the two-phase flow that new envelope spectrum is more suitable for tube bank than the envelope spectrum based on single-phase flow buffets calculating, reduce again the too high conservative property of existing two-phase flow envelope spectrum to a certain extent simultaneously, achieve heat interchanger turbulent buffeting RESPONSE CALCULATION analysis input more reasonable, analysis result is technique effect more accurately.

For solving the problems of the technologies described above, the embodiment of the present application provides the computing method of a kind of heat interchanger turbulent buffeting response, and described method comprises:

Utilize the flow field of thermal-hydraulic systematic analysis software heat exchanger secondary side to calculate, obtain the flow field parameter that described secondary side fluid distributes along heat-transfer pipe;

Based on the structural parameters of described flow field parameter and described heat-transfer pipe, judge whether described flow field parameter belongs to the first preset range and the second preset range;

If described flow field parameter belongs to described first preset range and described second preset range, then formula (1)-formula (6) is utilized to calculate the reference equivalent power spectral density function of exciting force

Based on described utilize random vibration theory, calculate the root-mean-square value of the vibratory response of turbulent buffeting, utilize described root-mean-square value to carry out Flow vibration evaluation;

Wherein, the envelope spectrum of dimensionless with reference to equivalent power spectral density of the two-phase flow exciting force in tube bank is acted on namely dimensionless with reference to the upper limit of equivalent power spectral density is:

When 10 ^-3≤ f _rduring <0.04, ${\left[{\stackrel{\&OverBar;}{\mathrm{\&Phi;}}}_E^0\right]}_U=4{f_R}^{-0.7}---\left(1\right)$

As 0.04≤f _rwhen≤1, ${\left[{\stackrel{\&OverBar;}{\mathrm{\&Phi;}}}_E^0\right]}_U=5\&times;{10}^{-4}{f_R}^{-3.5}---\left(2\right)$

Wherein:

${\left[{\stackrel{\&OverBar;}{\mathrm{\&Phi;}}}_E^0\right]}_U=\frac{V_i}{{\left({\mathrm{\&rho;}}_{l}g{D}_{w}D\right)}^2{D}_w}{\mathrm{\&Phi;}}_E^0\left(f\right)---\left(3\right)$

$f_R=\frac{{\mathrm{fD}}_w}{V_i}---\left(4\right)$

$D_w=\frac{0.1D}{\sqrt{1-{\mathrm{\&alpha;}}_H}}---\left(5\right)$

Interface flow velocity V _ibe defined as:

$V_i=0.73(J_g+J_l)+\sqrt{g{D}_e{\mathrm{\&rho;}}_l-{\mathrm{\&rho;}}_g/{\mathrm{\&rho;}}_l}---\left(6\right)$

ρ in formula _lfor the density of liquid phase, g is acceleration of gravity, and D is pipe diameter, effective diameter of pipe D _ebe defined as D _e=2 (P-D), P are the centre distance (pitch) of tube bank adjacent pipe, and f is the natural frequency of pipe, α _hfor homogeneous flow void fraction, J _gand J _lbe respectively the apparent velocity in Guan Zhen, for the equivalent reference power spectral density of exciting force on unit pipe range, for the dimensionless number of 0, D _wfor length dimension, f _rfor reduced frequency, i.e. a kind of dimensionless number of f.

Further, described first preset range is specially: homogeneous flow pitch flow velocity is between 0.2 ~ 14m/s, and described second preset range is specially: homogeneous flow void fraction is between 10 ~ 95%.

Further, described flow field parameter specifically comprises: void fraction, density and flow velocity.

Further, described structural parameters are specially: each first order mode of described heat-transfer pipe, frequency.

Further, described heat transfer tube configuration parameter can be carried out calculating and be obtained by the computing method provided in GB151 and TEMA Specification, or is calculated by commercial finite element software.

The one or more technical schemes provided in the embodiment of the present application, at least have following technique effect or advantage:

First the flow field of thermal-hydraulic systematic analysis software heat exchanger secondary side is utilized to calculate owing to have employed, obtain the flow field parameter that described secondary side fluid distributes along heat-transfer pipe, then theoretical analysis or finite element method is utilized, obtain the structural parameters of described heat-transfer pipe, then based on described flow field parameter, judge whether described flow field parameter belongs in the first preset range and described second preset range, if then described flow field parameter belongs to described first preset range and described structural parameters belong to the second preset range, formula (1)-formula (6) is then utilized to calculate the equivalent PSD of reference of exciting force, obtain last based on described utilize random vibration theory, calculate the root-mean-square value of the vibratory response of turbulent buffeting, described root-mean-square value is utilized to carry out the technological means of Flow vibration assay, namely the PSD envelope spectrum of new two-phase flow exciting force is obtained, the two-phase flow that new envelope spectrum is more suitable for tube bank than the envelope spectrum based on single-phase flow buffets calculating, reduce again the too high conservative property of existing two-phase flow envelope spectrum to a certain extent simultaneously, so, efficiently solve existing heat interchanger Flow vibration analytical approach existence analysis input unreasonable too high with conservative property, cause the inaccurate technical matters of analysis result, and then it is more reasonable to achieve heat interchanger turbulent buffeting RESPONSE CALCULATION analysis input, analysis result is technique effect more accurately.

Accompanying drawing explanation

Fig. 1 is the process flow diagram of Flow vibration analytical approach in the embodiment of the present application one;

Fig. 2 is the envelope spectrum of the two-phase flow exciting force PSD that de Langre and Villard provides;

Fig. 3 is that the two-phase flow exciting force dimensionless that proposes of analytical approach in the embodiment of the present application one is with reference to equivalent PSD envelope spectrum;

Fig. 4 is the illustraton of model of the steam generator heat-transfer pipe in the embodiment of the present application one;

Fig. 5 is the secondary side fluid transverse velocity acting on heat-transfer pipe bending section in the embodiment of the present application one;

Fig. 6 is the distribution of vibration displacement Root mean square response on heat-transfer pipe that in the embodiment of the present application one, three kinds of envelope spectrums calculate.

Embodiment

The invention provides the computing method of a kind of heat interchanger turbulent buffeting response, solve existing heat interchanger Flow vibration analytical approach existence analysis input unreasonable too high with conservative property, cause the inaccurate technical matters of analysis result, obtain the PSD envelope spectrum of new two-phase flow exciting force, the two-phase flow that new envelope spectrum is more suitable for tube bank than the envelope spectrum based on single-phase flow buffets calculating, reduce again the too high conservative property of existing two-phase flow envelope spectrum to a certain extent simultaneously, achieve heat interchanger turbulent buffeting RESPONSE CALCULATION analysis input more reasonable, analysis result is technique effect more accurately.

Technical scheme during the application implements is for solving the problems of the technologies described above.General thought is as follows:

Have employed and first utilize the flow field of thermal-hydraulic systematic analysis software heat exchanger secondary side to calculate, obtain the flow field parameter that described secondary side fluid distributes along heat-transfer pipe, then theoretical analysis or finite element method is utilized, obtain the structural parameters of described heat-transfer pipe, then based on described flow field parameter, judge whether described flow field parameter belongs in the first preset range and described second preset range, if then described flow field parameter belongs to described first preset range and described structural parameters belong to the second preset range, formula (1)-formula (6) is then utilized to calculate the equivalent PSD of reference of exciting force, obtain last based on described utilize random vibration theory, calculate the root-mean-square value of the vibratory response of turbulent buffeting, described root-mean-square value is utilized to carry out the technological means of Flow vibration assay, namely the PSD envelope spectrum of new two-phase flow exciting force is obtained, the two-phase flow that new envelope spectrum is more suitable for tube bank than the envelope spectrum based on single-phase flow buffets calculating, reduce again the too high conservative property of existing two-phase flow envelope spectrum to a certain extent simultaneously, so, efficiently solve existing heat interchanger Flow vibration analytical approach existence analysis input unreasonable too high with conservative property, cause the inaccurate technical matters of analysis result, and then it is more reasonable to achieve heat interchanger turbulent buffeting RESPONSE CALCULATION analysis input, analysis result is technique effect more accurately.

In order to better understand technique scheme, below in conjunction with Figure of description and concrete embodiment, technique scheme is described in detail.

Embodiment one:

In embodiment one, provide the computing method of a kind of heat interchanger turbulent buffeting response, please refer to Fig. 1-Fig. 6, described method comprises:

S10, utilizes the flow field of thermal-hydraulic systematic analysis software heat exchanger secondary side to calculate, obtains the flow field parameter that described secondary side fluid distributes along heat-transfer pipe;

After step slo, the method for the embodiment of the present application just enters step S20, that is: based on the structural parameters of described flow field parameter and described heat-transfer pipe, judge whether described flow field parameter belongs to the first preset range and the second preset range;

After step S20, the method of the embodiment of the present application just enters step S30, that is: if described flow field parameter belongs to described first preset range and described second preset range, then formula (1)-formula (6) is utilized to calculate the reference equivalent power spectral density function of exciting force

After step S30, the method for the embodiment of the present application just enters step S40, that is: based on described utilize random vibration theory, calculate the root-mean-square value of the vibratory response of turbulent buffeting, utilize described root-mean-square value to carry out Flow vibration evaluation;

Wherein, the envelope spectrum of dimensionless with reference to equivalent power spectral density of the two-phase flow exciting force in tube bank is acted on namely dimensionless with reference to the upper limit of equivalent power spectral density is:

When 10 ^-3≤ f _rduring <0.04, ${\left[{\stackrel{\&OverBar;}{\mathrm{\&Phi;}}}_E^0\right]}_U=4{f_R}^{-0.7}---\left(1\right)$

As 0.04≤f _rwhen≤1, ${\left[{\stackrel{\&OverBar;}{\mathrm{\&Phi;}}}_E^0\right]}_U=5\&times;{10}^{-4}{f_R}^{-3.5}---\left(2\right)$

Wherein:

${\left[{\stackrel{\&OverBar;}{\mathrm{\&Phi;}}}_E^0\right]}_U=\frac{V_i}{{\left({\mathrm{\&rho;}}_{l}g{D}_{w}D\right)}^2{D}_w}{\mathrm{\&Phi;}}_E^0\left(f\right)---\left(3\right)$

$f_R=\frac{{\mathrm{fD}}_w}{V_i}---\left(4\right)$

$D_w=\frac{0.1D}{\sqrt{1-{\mathrm{\&alpha;}}_H}}---\left(5\right)$

Interface flow velocity V _ibe defined as:

$V_i=0.73(J_g+J_l)+\sqrt{g{D}_e{\mathrm{\&rho;}}_l-{\mathrm{\&rho;}}_g/{\mathrm{\&rho;}}_l}---\left(6\right)$

ρ in formula _lfor the density of liquid phase, g is acceleration of gravity, and D is pipe diameter, effective diameter of pipe D _ebe defined as D _e=2 (P-D), P are the centre distance (pitch) of tube bank adjacent pipe, and f is the natural frequency of pipe, α _hfor homogeneous flow void fraction, J _gand J _lbe respectively the apparent velocity in Guan Zhen, for the equivalent reference power spectral density of exciting force on unit pipe range, for dimensionless number, D _wfor length dimension, f _rfor reduced frequency, i.e. a kind of dimensionless number of f.

Wherein, in the embodiment of the present application, described first preset range is specially: homogeneous flow pitch flow velocity is between 0.2 ~ 14m/s, and described second preset range is specially: homogeneous flow void fraction is between 10 ~ 95%.

Wherein, in the embodiment of the present application, described flow field parameter specifically comprises: void fraction, density and flow velocity.

Wherein, in the embodiment of the present application, described structural parameters are specially: each first order mode of described heat-transfer pipe, frequency.

Wherein, in the embodiment of the present application, described heat transfer tube configuration parameter can be carried out calculating and be obtained by the computing method provided in GB151 and TEMA Specification, or is calculated by commercial finite element software.

Wherein, in actual applications, commercial finite element software comprises: ANSYS, ABAQUS etc.

Wherein, in the embodiment of the present application, the dimensionless of two-phase flow exciting force is shown in formula (1) and formula (2) with reference to the concrete form of equivalent PSD envelope spectrum, and wherein the flex point of envelope spectrum corresponds to f _r=0.04; The coefficient of low-frequency range PSD envelope spectrum is 4, and index is-0.7; The coefficient of high band PSD envelope spectrum is 5 × 10 ^-4, index is-3.5.

Illustrate below and the Flow vibration analytical approach in the embodiment of the present application is specifically introduced, the PSD response of the U-shaped heat-transfer pipe of steam generator is researched and analysed.

The physical dimension of table 1 heat-transfer pipe

Fig. 4 is the illustraton of model of steam generator heat-transfer pipe, wherein, being node and block plan in the left-half of the horizontal direction of Fig. 4, is joint constraint condition at the right half part of the horizontal direction of Fig. 4, and the physical dimension that heat-transfer pipe model is corresponding and material parameter are as shown in table 1-table 2.Then utilize the two-phase flow flow field of thermal-hydraulic systematic analysis software to heat-transfer pipe present position to calculate, the basic parameter of Flow Field Calculation is shown in Table 3.Calculate the secondary side fluid transverse velocity of heat-transfer pipe bending section, as shown in Figure 5.By judging, the data of the steam generator heat-transfer pipe of this example meet the scope of application of envelope spectrum in the embodiment of the present application: homogeneous flow void fraction 10 ~ 95%; Homogeneous flow pitch flow velocity (i.e. V _p) 0.2 ~ 14m/s.Then utilize formula 1-formula 6 to calculate the equivalent PSD of reference of exciting force, obtain recycling random vibration theory, calculates the root-mean-square value of the vibratory response of turbulent buffeting.The distribution of the vibration displacement Root mean square response that three kinds of envelope spectrums calculate on heat-transfer pipe as shown in Figure 6, as can be seen from Figure 6 between the dynamic respond that the analytical approach that the dynamic respond that the Flow vibration analytical approach in application the embodiment of the present application obtains provides between dynamic respond and application de Langre and Villard based on single-phase flow envelope spectrum obtains, reduce the conservative property that existing envelope spectrum is too high, dynamic respond can be further used for wearing and tearing and the analysis of fatigue of heat-transfer pipe.

Technical scheme in above-mentioned the embodiment of the present application, at least has following technique effect or advantage:

First the flow field of thermal-hydraulic systematic analysis software heat exchanger secondary side is utilized to calculate owing to have employed, obtain the flow field parameter that described secondary side fluid distributes along heat-transfer pipe, then theoretical analysis or finite element method is utilized, obtain the structural parameters of described heat-transfer pipe, then based on described flow field parameter, judge whether described flow field parameter belongs in the first preset range and described second preset range, if then described flow field parameter belongs to described first preset range and described structural parameters belong to the second preset range, formula (1)-formula (6) is then utilized to calculate the equivalent PSD of reference of exciting force, obtain last based on described utilize random vibration theory, calculate the root-mean-square value of the vibratory response of turbulent buffeting, described root-mean-square value is utilized to carry out the technological means of Flow vibration assay, namely the PSD envelope spectrum of new two-phase flow exciting force is obtained, the two-phase flow that new envelope spectrum is more suitable for tube bank than the envelope spectrum based on single-phase flow buffets calculating, reduce again the too high conservative property of existing two-phase flow envelope spectrum to a certain extent simultaneously, so, efficiently solve existing heat interchanger Flow vibration analytical approach existence analysis input unreasonable too high with conservative property, cause the inaccurate technical matters of analysis result, and then it is more reasonable to achieve heat interchanger turbulent buffeting RESPONSE CALCULATION analysis input, analysis result is technique effect more accurately.

Although describe the preferred embodiments of the present invention, those skilled in the art once obtain the basic creative concept of cicada, then can make other change and amendment to these embodiments.So claims are intended to be interpreted as comprising preferred embodiment and falling into all changes and the amendment of the scope of the invention.

Obviously, those skilled in the art can carry out various change and modification to the present invention and not depart from the spirit and scope of the present invention.Like this, if these amendments of the present invention and modification belong within the scope of the claims in the present invention and equivalent technologies thereof, then the present invention is also intended to comprise these change and modification.

Claims (5)

1. computing method for heat interchanger turbulent buffeting response, it is characterized in that, described method comprises:

Utilize the flow field of thermal-hydraulic systematic analysis software heat exchanger secondary side to calculate, obtain the flow field parameter that described secondary side fluid distributes along heat-transfer pipe;

Based on the structural parameters of described flow field parameter and described heat-transfer pipe, judge whether described flow field parameter belongs to the first preset range and the second preset range;

If described flow field parameter belongs to described first preset range and described second preset range, then formula (1)-formula (6) is utilized to calculate the reference equivalent power spectral density function of exciting force

Based on described utilize random vibration theory, calculate the root-mean-square value of the vibratory response of turbulent buffeting, utilize described root-mean-square value to carry out Flow vibration evaluation;

Wherein, the envelope spectrum of dimensionless with reference to equivalent power spectral density of the two-phase flow exciting force in tube bank is acted on namely dimensionless with reference to the upper limit of equivalent power spectral density is:

When 10 ^-3≤ f _rduring <0.04, ${\left[{\stackrel{\&OverBar;}{\mathrm{\&Phi;}}}_E^0\right]}_U=4{f_R}^{-0.7}---\left(1\right)$

As 0.04≤f _rwhen≤1, ${\left[{\stackrel{\&OverBar;}{\mathrm{\&Phi;}}}_E^0\right]}_U=5\&times;{10}^{-4}{f_R}^{-3.5}---\left(2\right)$

Wherein:

${\left[{\stackrel{\&OverBar;}{\mathrm{\&Phi;}}}_E^0\right]}_U=\frac{V_i}{{\left({\mathrm{\&rho;}}_{l}g{D}_{w}D\right)}^2{D}_w}{\mathrm{\&Phi;}}_E^0\left(f\right)---\left(3\right)$

$f_R=\frac{f{D}_w}{V_i}---\left(4\right)$

$D_w=\frac{0.1D}{\sqrt{1-{\mathrm{\&alpha;}}_H}}---\left(5\right)$

Interface flow velocity V _ibe defined as:

$V_i=0.73(J_g+J_l)+\sqrt{g{D}_e({\mathrm{\&rho;}}_l-{\mathrm{\&rho;}}_g)/{\mathrm{\&rho;}}_l}---\left(6\right)$

ρ in formula _lfor the density of liquid phase, g is acceleration of gravity, and D is pipe diameter, effective diameter of pipe D _ebe defined as D _e=2 (P-D), P are the centre distance (pitch) of tube bank adjacent pipe, and f is the natural frequency of pipe, α _hfor homogeneous flow void fraction, J _gand J _lbe respectively the apparent velocity in Guan Zhen, for the equivalent reference power spectral density of exciting force on unit pipe range, for dimensionless number, D _wfor length dimension, f _rfor reduced frequency, i.e. a kind of dimensionless number of f.

2. method according to claim 1, is characterized in that, described first preset range is specially: homogeneous flow pitch flow velocity is between 0.2 ~ 14m/s, and described second preset range is specially: homogeneous flow void fraction is between 10 ~ 95%.

3. method according to claim 1, is characterized in that, described flow field parameter specifically comprises: void fraction, density and flow velocity.

4. method according to claim 1, is characterized in that, described structural parameters are specially: each first order mode of described heat-transfer pipe, frequency.

5. method according to claim 1, is characterized in that, described heat transfer tube configuration parameter can be carried out calculating and be obtained by the computing method provided in GB151 and TEMA Specification, or is calculated by commercial finite element software.