CN109635500A - Aviation pipeline three-dimensional flow consolidates coupling parameter resonance response characteristic prediction method and device - Google Patents
Aviation pipeline three-dimensional flow consolidates coupling parameter resonance response characteristic prediction method and device Download PDFInfo
- Publication number
- CN109635500A CN109635500A CN201910001794.6A CN201910001794A CN109635500A CN 109635500 A CN109635500 A CN 109635500A CN 201910001794 A CN201910001794 A CN 201910001794A CN 109635500 A CN109635500 A CN 109635500A
- Authority
- CN
- China
- Prior art keywords
- aviation
- conveying pipe
- response characteristic
- equation
- fluid conveying
- 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.)
- Granted
Links
- 230000004044 response Effects 0.000 title claims abstract description 93
- 230000008878 coupling Effects 0.000 title claims abstract description 43
- 238000010168 coupling process Methods 0.000 title claims abstract description 43
- 238000005859 coupling reaction Methods 0.000 title claims abstract description 43
- 238000000034 method Methods 0.000 title claims abstract description 39
- 239000012530 fluid Substances 0.000 claims abstract description 136
- 230000033001 locomotion Effects 0.000 claims abstract description 41
- 239000011159 matrix material Substances 0.000 claims description 11
- 230000010349 pulsation Effects 0.000 claims description 8
- 238000006243 chemical reaction Methods 0.000 claims description 3
- 238000010586 diagram Methods 0.000 description 8
- 238000006073 displacement reaction Methods 0.000 description 6
- 230000005284 excitation Effects 0.000 description 6
- 238000013016 damping Methods 0.000 description 5
- 230000009466 transformation Effects 0.000 description 4
- 238000005452 bending Methods 0.000 description 3
- 230000008901 benefit Effects 0.000 description 3
- 239000000463 material Substances 0.000 description 3
- 229910000831 Steel Inorganic materials 0.000 description 2
- 230000005540 biological transmission Effects 0.000 description 2
- 230000003068 static effect Effects 0.000 description 2
- 239000010959 steel Substances 0.000 description 2
- 238000006467 substitution reaction Methods 0.000 description 2
- 238000004026 adhesive bonding Methods 0.000 description 1
- 210000001367 artery Anatomy 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 230000006378 damage Effects 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 238000009792 diffusion process Methods 0.000 description 1
- 238000005315 distribution function Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 239000000446 fuel Substances 0.000 description 1
- 230000005484 gravity Effects 0.000 description 1
- 239000004615 ingredient Substances 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 238000005086 pumping Methods 0.000 description 1
- 238000013139 quantization Methods 0.000 description 1
- 238000010791 quenching Methods 0.000 description 1
- 230000003014 reinforcing effect Effects 0.000 description 1
- 238000000926 separation method Methods 0.000 description 1
- 210000003462 vein Anatomy 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Measuring Volume Flow (AREA)
- Pipeline Systems (AREA)
Abstract
The present invention relates to the dynamics Design technical field of aviation fluid conveying pipe, propose that a kind of aviation pipeline three-dimensional flow consolidates coupling parameter resonance response characteristic prediction method and device.The aviation pipeline three-dimensional flow consolidates coupling parameter resonance response characteristic prediction method and consolidates Coupling Dynamic Model according to the three-dimensional flow that relevant parameter establishes aviation fluid conveying pipe；Consolidate the Coupled Dynamics equation of motion according to the three-dimensional flow that kinetic model establishes aviation fluid conveying pipe；Finite dimensioned equation is obtained according to dynamical motion equation；Parameter resonance response characteristic of the aviation fluid conveying pipe within the scope of default ripple frequency is obtained according to finite dimensioned equation.Parameter resonance response characteristic has been obtained, interval of resonance can have effectively been avoided, prevent parameter resonance, has improved the service life and safety of aviation fluid conveying pipe.
Description
Technical field
The present invention relates to the dynamics Design technical field of aviation fluid conveying pipe more particularly to a kind of aviation pipeline three-dimensional flows
Gu coupling parameter resonance response characteristic prediction method and device.
Background technique
Hydraulic air, fuel passage system have the important function of medium transmission and energy transmission.Its load environment is compared
It is more special for traditional industry fluid conveying pipe.Make in complicated internal pressure pulsation and the comprehensive of external structural vibrations excitation
With.The variation of tube fluid motion state can interact with structural pipe wall, can generate fluid structurecoupling phenomenon, induce pipeline configuration
It vibrates.When especially motivating trandfer fluid by pumping source, the flow velocity of pipeline is generally not ideal stability, and aviation pipeline is pumped
Alternation suction oil movement, tube fluid is the flow velocity of pulses in most cases.When flow velocity, ripple frequency and pulsation amplitude
When meeting some requirements, parameter resonance will occur for fluid conveying pipe.Be embodied in, the negative stiffness item of system be time-varying and
Its frequency is at twice of system frequency or covibration that the sum of certain two rank intrinsic frequency place occurs.
This parameter resonance can cause tube noise diffusion, and pipeline life reduces and the destruction of pipe fitting, by serious shadow
Ring the safety for arriving aviation Pipe Conveying Fluid.
Therefore, it is necessary to propose a kind of new aviation pipeline three-dimensional flow consolidate coupling parameter resonance response characteristic prediction method and
Device.
Above- mentioned information disclosed in the background technology part are only used for reinforcing the understanding to background of the invention, therefore it can
To include the information not constituted to the prior art known to persons of ordinary skill in the art.
Summary of the invention
It is an object of the invention to overcome the unpredictable parameter resonance response characteristic of above-mentioned prior art aviation fluid conveying pipe
And lead to the lower deficiency of the safety of aviation Pipe Conveying Fluid, it is total to provide a kind of parameter that can predict aviation fluid conveying pipe
It is pre- that aviation pipeline three-dimensional flow of the vibration response characteristic to improve aviation fluid conveying pipe safety consolidates coupling parameter resonance response characteristic
Survey method and apparatus.
Additional aspect and advantage of the invention will be set forth in part in the description, and partly will be from description
It is apparent from, or can practice through the invention and acquistion.
According to an aspect of the present invention, a kind of aviation pipeline three-dimensional flow consolidates coupling parameter resonance response Predicting Performance Characteristics side
Method, comprising:
Consolidate Coupling Dynamic Model according to the three-dimensional flow that relevant parameter establishes aviation fluid conveying pipe；
Consolidate the Coupled Dynamics equation of motion according to the three-dimensional flow that the kinetic model establishes the aviation fluid conveying pipe；
Finite dimensioned equation is obtained according to the dynamical motion equation；
Parameter of the aviation fluid conveying pipe within the scope of default ripple frequency is obtained according to the finite dimensioned equation
Resonance response characteristic.
In a kind of exemplary embodiment of the disclosure, the relevant parameter includes:
Geometric parameter, internal diameter, outer diameter and total length including the fluid conveying pipe；
Physical parameter, elasticity modulus and density including the fluid conveying pipe；
Duty parameter, including the mean flow rate of fluid, pulsation amplitude and ripple frequency in the fluid conveying pipe.
In a kind of exemplary embodiment of the disclosure, finite dimensioned side is obtained according to the dynamical motion equation
Journey, comprising:
The finite dimensioned equation is obtained using the golden method of gal the Liao Dynasty according to the dynamical motion equation.
In a kind of exemplary embodiment of the disclosure, the aviation pipe conveying fluid is obtained according to the finite dimensioned equation
The parameter resonance response characteristic in road, comprising:
The parameter resonance response characteristic of fluid conveying pipe is calculated using Runge Kutta method according to the finite dimensioned equation.
In a kind of exemplary embodiment of the disclosure, the parameter resonance response characteristic includes:
Aviation fluid conveying pipe fluid velocity inside three-dimensional flow of aviation fluid conveying pipe when being lower than critical flow velocity consolidates coupling parameter
The three-dimensional fluid structurecoupling of aviation fluid conveying pipe when resonance response characteristic and aviation fluid conveying pipe fluid velocity inside are more than critical flow velocity
Parameter resonance response characteristic.
In a kind of exemplary embodiment of the disclosure, the parameter resonance response characteristic is within the scope of default ripple frequency
The vibration amplitude of aviation fluid conveying pipe corresponding to different ripple frequencies.
In a kind of exemplary embodiment of the disclosure, the aviation pipeline three-dimensional flow consolidates coupling parameter resonance response characteristic
Prediction technique further include:
Eigenvalues analysis is carried out to the dynamical matrix equation and obtains intrinsic frequency.
In a kind of exemplary embodiment of the disclosure, the intrinsic frequency includes:
The aviation fluid conveying pipe fluid velocity inside intrinsic frequency of aviation fluid conveying pipe and aviation when being lower than critical flow velocity is defeated
The intrinsic frequency of aviation fluid conveying pipe when flow tube road fluid velocity inside is more than critical flow velocity.
In a kind of exemplary embodiment of the disclosure, the parameter resonance response characteristic is detected according to the intrinsic frequency
Accuracy.
According to one aspect of the disclosure, a kind of aviation pipeline three-dimensional flow is provided and consolidates the pre- of coupling parameter resonance response characteristic
Survey device, comprising:
Model building module, the three-dimensional flow for establishing aviation fluid conveying pipe consolidate Coupling Dynamic Model；
Establishing equation module consolidates Coupled Dynamics matrix for establishing aviation pipeline three-dimensional flow according to the kinetic model
Equation；
Equation conversion module, for obtaining finite dimensioned equation according to dynamical matrix equation；
Computing module, the parameter resonance for obtaining the aviation fluid conveying pipe according to the finite dimensioned equation respond
Characteristic.
As shown from the above technical solution, the present invention has at least one of following advantages and good effect:
Aviation pipeline three-dimensional flow of the present invention consolidates coupling parameter resonance response characteristic prediction method, is established and is navigated according to relevant parameter
The three-dimensional flow of empty fluid conveying pipe consolidates Coupling Dynamic Model；Three-dimensional flow, which is established, according to kinetic model consolidates Coupled Dynamics movement side
Journey, and then finite dimensioned equation is obtained, aviation fluid conveying pipe is then obtained in default pulsation frequency by finite dimensioned equation
Parameter resonance response characteristic within the scope of rate；Parameter resonance response characteristic has been obtained, interval of resonance can be effectively avoided, has prevented from sending out
Raw parameter resonance, improves the service life and safety of aviation fluid conveying pipe.
Detailed description of the invention
Its example embodiment is described in detail by referring to accompanying drawing, above and other feature of the invention and advantage will become
It is more obvious.
Fig. 1 is the flow chart that aviation pipeline three-dimensional flow of the present invention consolidates coupling parameter resonance response characteristic prediction method；
Fig. 2 is the schematic diagram of relevant parameter；
Fig. 3 is the solution procedure schematic diagram of intrinsic frequency；
Fig. 4 is the classification chart of parameter resonance response characteristic；
Fig. 5 is that the three-dimensional flow of aviation fluid conveying pipe consolidates Coupling Dynamic Model schematic diagram；
Fig. 6 is the schematic diagram of intrinsic frequency before and after aviation fluid conveying pipe plane buckling；
Fig. 7 is the schematic diagram of intrinsic frequency before and after aviation fluid conveying pipe out-of-plane buckling；
Fig. 8 is amplitude-frequency response before buckling；
Fig. 9 to Figure 19 is time and space response cloud atlas；
Figure 20 to Figure 28 is that ripple frequency is ωp=15 response diagram；
Figure 29 to Figure 36 is that ripple frequency is ωp=45 response diagram；
Figure 37 to Figure 41 is that ripple frequency is ωp=45 motion profile figure；
Figure 42 to Figure 46 is that ripple frequency is ωp=400 motion profile figure；
Figure 47 is amplitude-frequency response after buckling；
Figure 48 to Figure 56 is that ripple frequency is ωp=14 response diagram；
Figure 57 to Figure 65 is that ripple frequency is ωp=395 response.
Specific embodiment
Example embodiment is described more fully with reference to the drawings.However, example embodiment can be with a variety of shapes
Formula is implemented, and is not understood as limited to embodiment set forth herein；On the contrary, thesing embodiments are provided so that the present invention will
Fully and completely, and by the design of example embodiment comprehensively it is communicated to those skilled in the art.Identical attached drawing in figure
Label indicates same or similar structure, thus the detailed description that will omit them.
Present invention firstly provides a kind of aviation pipeline three-dimensional flows to consolidate coupling parameter resonance response characteristic prediction method.Reference
Shown in Fig. 1, which, which consolidates coupling parameter resonance response characteristic prediction method, be may comprise steps of:
Step S110 consolidates Coupling Dynamic Model according to the three-dimensional flow that relevant parameter establishes aviation fluid conveying pipe；
Step S120 consolidates Coupled Dynamics fortune according to the three-dimensional flow that the kinetic model establishes the aviation fluid conveying pipe
Dynamic equation；
Step S130 obtains finite dimensioned equation according to the dynamical motion equation；
Step S140 obtains the aviation fluid conveying pipe in default ripple frequency range according to the finite dimensioned equation
Interior parameter resonance response characteristic.
Above-mentioned steps are described in detail below.
In step s 110, Coupling Dynamic Model is consolidated according to the three-dimensional flow that relevant parameter establishes aviation fluid conveying pipe.
In the step s 120, the three-dimensional fluid structurecoupling power of the aviation fluid conveying pipe is established according to the kinetic model
Learn the equation of motion.
Referring to shown in Fig. 2, relevant parameter may include that geometric parameter, physical parameter, duty parameter geometric parameter can wrap
Include the internal diameter, outer diameter and total length of fluid conveying pipe；Physical parameter may include the elasticity modulus and density of the fluid conveying pipe；Work
Condition parameter may include the mean flow rate of fluid in the fluid conveying pipe, pulsation amplitude and ripple frequency.
Duct length is L, linear mass mp, Young's modulus E, cross sectional moment of inertia I, pipeline section product is Ap,
Fluid units linear mass is mf, fluid velocity U.The equilibrium state of pipeline is along x-axis, it is assumed that and pipeline is Euler's beam model,
Two, y-o-z plane orthogonal displacements are respectively W (x, t) and V (x, t).
In this example embodiment, it then follows it is assumed hereinafter that: pipe diameter relative duct size length is very small, meets euler beam
Model application conditions, pipeline is shear-deformable and rotary inertia is ignored；Fluid is incompressible without gluing；Gravity influence is neglected
Slightly；Pipeline external axial force and Fluid pressure are not considered.
Pipeline considers the influence of viscoelastic damping, and the viscoelastic damping coefficient of pipeline is E*.Axial force in pipeline is
Wherein ε is the axial strain of pipeline,It is axial strain to the first derivative of time.Consider the shadow of geometrical non-linearity
It rings, the axial strain of pipeline can be expressed as
It is by the expression formula that (2) formula substitution (1) formula obtains axial force
(3) are brought into the non-linear three-dimensional fortune for removing to obtain fluid conveying pipe in the equation of motion of classical fluid conveying pipe
Dynamic equation:
Fluid flow rate is fluctuating velocity in system
U=U0(1+μcos(Ωpt)) (6)
Referring to shown in Fig. 3, for convenience, following dimensionless group is introduced:
ξ representation space coordinate, w and v indicate nondimensionalization position, and τ indicates the nondimensionalization time, and u is nondimensionalization flow velocity,
β is mass ratio, and ω indicates nondimensionalization frequency, α nondimensionalization damped coefficient, a kind of quantization of duct length and system when S,
Above-mentioned dimensionless group is definite value in one timing of pipeline.
Dimensionless group is updated to (16), the nondimensionalization equation of motion can be obtained in (17) formula:
Both ends simple boundary condition is
η may include w and v in formula, be displacement parameter.
In step s 130, finite dimensioned equation is obtained according to the dynamical motion equation.
In step S140, the aviation fluid conveying pipe is obtained in default ripple frequency according to the finite dimensioned equation
Parameter resonance response characteristic in range.
Referring to shown in Fig. 3, aviation pipeline three-dimensional flow of the present invention consolidates coupling parameter resonance response characteristic prediction method can be with
Aviation fluid conveying pipe fluid velocity inside is obtained lower than critical flow including carrying out Eigenvalues analysis to the dynamical matrix equation
Aviation fluid conveying pipe when the intrinsic frequency of aviation fluid conveying pipe and aviation fluid conveying pipe fluid velocity inside are more than critical flow velocity when fast
Intrinsic frequency.
Referring to shown in Fig. 4, the aviation fluid conveying pipe is obtained in default ripple frequency according to the finite dimensioned equation
Parameter resonance response characteristic in range may include obtaining aviation fluid conveying pipe fluid velocity inside boat when lower than critical flow velocity
The three-dimensional flow of empty fluid conveying pipe consolidates coupling parameter resonance response characteristic and aviation fluid conveying pipe fluid velocity inside is more than critical flow velocity
When aviation fluid conveying pipe three-dimensional flow consolidate coupling parameter resonance response characteristic.
Finite dimensioned equation, equation (16), (17) are obtained using the golden method of gal the Liao Dynasty according to the dynamical matrix equation
The three-dimensional motion model of Pipe Conveying Fluid is described, Garlerkin method is used to the discretization model, two, pipeline lateral positions
Shifting can set solution as
Wherein wr(τ) and vr(τ) is the generalized coordinates of discrete system, and N is that number, φ is truncated in Garlerkinr(ξ) is both ends
The mode function of freely-supported beam model bending vibration
(19), (20) and (21) are substituted into (16), (17) formula, second nonlinear is obtained using Garlerkin method and simplifies
Model
Wherein
It is generalized coordinates to the second dervative of nondimensional time coordinate,For broad sense seat
Mark the first derivative to nondimensional time coordinate, wj(τ)、vj(τ) is generalized coordinates.
In this example embodiment, Garlerkin truncation number takes 6, and generalized coordinates can be set to following column vector
Form
W=[w1 w2 w3 w4 w5 w6]T (28)
V=[v1 v2 v3 v4 v5 v6]T (29)
Equation (22), (23) can be written as following matrix equation
Equation obtains after arranging
Introduce state variable u1, u2, u3, u4,
It will be in (34) substitution (32), (33)
Quadravalence Long Gekuta Iteration:
H is iteration step length, T0It is iteration starting point, therefore has
t1The first derivative of the state variable at moment
t2The state variable at moment
t2The first derivative of the state variable at moment
t3The state variable at moment
t3The first derivative of the state variable at moment
t4The state variable at moment
State variable after an iteration
The interative computation for carrying out 1,000,000 times obtains stablizing u2、u4, and then obtain the vibration of stable aviation fluid conveying pipe
Amplitude, i.e. u2、u4The value of vector sum.
Below by way of a specific embodiment, the present invention will be described:
In this example embodiment, iteration step length=0.0005 is taken, iteration total degree is 1000000 times.
During numerical value calculates, steel pipe is chosen as object, pipeline Young's modulus is E=210GPa, channel density ρp=
7850kg/m3, outer diameter tube Do=12mm, internal diameter Di=10mm, duct length L=1m, fluid density ρf=
870kg/m3。
To solve (fluid flow rate be less than critical flow velocity) fluid conveying pipe intrinsic frequency before buckling, ignore geometrical non-linearity item and
Viscoelastic damping item, equation (16) and (17) can simplify for
Referring to shown in Fig. 3, aviation pipeline three-dimensional flow of the present invention consolidates coupling parameter resonance response characteristic prediction method can be with
Aviation fluid conveying pipe fluid velocity inside is obtained lower than critical flow including carrying out Eigenvalues analysis to the dynamical matrix equation
Aviation fluid conveying pipe when the intrinsic frequency of aviation fluid conveying pipe and aviation fluid conveying pipe fluid velocity inside are more than critical flow velocity when fast
Intrinsic frequency.Equation (44) and (45) are handled according to Galerkin method, then carries out Eigenvalues analysis and acquires system
Here iteration step length h=0.0005 is taken, iteration total degree is 1000000 times.
During numerical value calculates, steel pipe is chosen as object, pipeline Young's modulus is E=210GPa, channel density ρp=
7850kg/m3, outer diameter tube Do=12mm, internal diameter Di=10mm, duct length L=1m, fluid density ρP=
870kg/m3。
To solve fluid conveying pipe intrinsic frequency before buckling, ignore geometrical non-linearity item and viscoelastic damping item, by equation
(44) it is handled with (45) according to Galerkin method, then carries out the intrinsic frequency that Eigenvalues analysis acquires system.
When fluid flow rate is more than critical flow velocity (first natural frequency is equal to zero), bending deformation can occur for fluid conveying pipe, by
In the presence of geometrical non-linearity, system does not occur unstability and breaks ring, but balances in a curved configuration states, herein only
Single order buckling is discussed, the displacement of buckling can be set to
The displacement amplitude in the direction wherein y, z is a1、a2, Displacements Distribution function is first-order modal function
The buckling configuration of (fluid flow rate is less than critical flow velocity) fluid conveying pipe after solving buckling, omit equation (16) and
(17) the time correlation item in obtains following equation:
(46) and (47) are updated in (48) and (49) and are acquired
In order to solve the intrinsic frequency of fluid conveying pipe after buckling, the displacement of fluid conveying pipe assume that comprising static part and
Dynamic part
Wherein w0(ξ)、v0(ξ) indicates static part,WithIndicate dynamic part.
Formula (51) is updated to the linear oscillator equation of fluid conveying pipe after going to obtain buckling in equation (16) and (17)
Garlerkin method discretization equation (52) and (53) then carry out Eigenvalues analysis and obtain pipe conveying fluid after buckling
The intrinsic frequency in road.
Shown in referring to figure 6 and figure 7, after bending deformation occurs for pipeline, there are two kinds of feelings outside in face and face for the vibration of pipeline
Shape, when velocity of flow of zero dimension is zero, first three rank dimensionless intrinsic frequency is ωn1=9.87, ωn2=39.48 and ωn3=88.83.
It is preceding 6 rank dimensionless intrinsic frequency under u=2 and u=3.5 situation that table 1, which gives velocity of flow of zero dimension,.
Table 1: fluid velocity is the dimensionless intrinsic frequency of u=2 and u=3.5
Fluid velocity is fluctuation velocity
U=u0(1+μcos(ωpτ)) (54)
Wherein u0It is mean flow rate, μ is pulsation amplitude, ωpFor ripple frequency.Runge Kutta method solves fluid conveying pipe
Response, primary condition areV (1)=[0, -0.0001,0,0,0,0] ', damping system
Number is α1=0.0001, dimensionless fluid velocity is u0=2, pulsation amplitude is μ=0.2.
Referring to shown in Fig. 8, ripple frequency is 13 < ωp<17(r1-1), 43 < ωp<47 (r1-2), 71 < ωp<79(r2-2),
121<ωp<128(r2-3), 166 < ωp<182 (r3-3), 237 < ωp<248(r3-4), 298 < ωp<327(r4-4), 393 < ωp<
408(r4-5), 469 < ωp<512(r5-5), 591 < ωp<607 (r5-6), 681 < ωp<740(r6-6), fluid conveying pipe can join
Number resonance, as ripple frequency increases continuously, subharmonic resonance and combination resonance are alternately present.
Referring to shown in Fig. 9 to Figure 19, r1-1, r2-2, r3-3, r4-4, r5-5, r6-6, r1-2, r2-3, r3-4, r4-5, r5-6Subharmonic
Resonance and combination resonance are alternately present, it can be seen that subharmonic resonance is the vibration of single mode, and combination resonance is continuous
The superposition of two rank mode.
The selecting frequency point ω on amplitude-frequency responsep=15, the response of the total strong point in recording distance end 3/10.
It is ω referring to ripple frequencyp=15 response: Figure 20 and Figure 21 is steady-state response, Fu that Figure 22 and Figure 23 are in response to
In leaf transformation, Figure 24 and Figure 25 are phasor and Poincare section, and Figure 26 is motion profile, and Figure 27 and Figure 28 are two transverse directions
Snapshot plotting.
Referring to shown in Figure 20 to Figure 28, ripple frequency ωpWhen=15, excitation has played the first step mode period vibration of system
It is dynamic, and motion profile is in one plane, the tilt angle of plane is related with primary condition.
The transverse movement snapshot plotting of the above subharmonic resonance show fluid conveying pipe axially each point response exist difference,
That is the stationary wave vibration that is not standard of the vibration of pipeline at this time, vibration here further includes a small amount of traveling wave ingredient.
Referring to shown in Figure 29 to Figure 41, the selecting frequency point ω on amplitude-frequency responsep=45 and ωp=400, record away from
Response from the total strong point in end 3/10.
Ripple frequency is ωp=45 response: Figure 29 and Figure 30 is steady-state response, in Fu that Figure 31 and Figure 32 are in response to
Leaf transformation, Figure 33 and Figure 34 are phasor and Poincare section, and Figure 35 and Figure 36 are two transverse direction snapshot plottings.
Ripple frequency is ωp=45 motion profile figure: in Tu37Zhong, ξ=1/10, Tu38Zhong, ξ=3/10, Fig. 9, ξ=
5/10, Tu40Zhong, ξ=7/10, Tu41Zhong, ξ=9/10.
Ripple frequency is ωpWhen=45, excitation has played the first rank and the second-order modal vibration of system, and ripple frequency is about etc.
In the sum of response frequency ωp≈ω1+ω2≈ 8.795+36.44, response frequency are approximately equal to the first rank and second-order intrinsic frequency
ωn1=7.55 and ωn)=37.51.Two transverse direction movement snapshot plotting display pipes movements are in the first first order mode and second-order
It is mutually converted between the vibration shape.
Motion profile is spatial form, and the motion profile at each position is all to have 5 cusps, this is that have with response frequency
It closes, the ratio between response frequency is ω1/ω)=0.241, it is approximately equal to 1/4, the number of cusp is equal to the sum of molecule and denominator of ratio.
Referring to shown in Figure 42 to Figure 46, ripple frequency ωp=400 motion profile figure: Tu42Zhong, ξ=1/10, figure
In 43, ξ=3/10, Tu44Zhong, ξ=5/10, Tu45Zhong, ξ=7/10, Tu46Zhong, ξ=9/10.
Ripple frequency is ωpWhen=400, excitation has played the fourth order and the 5th rank modal vibration of system, and ripple frequency is about
Equal to the sum of response frequency ωp≈ω4+ω5≈ 158.3+242.5, response frequency are approximately equal to fourth order and the 5th intrinsic frequency of rank
Rate ωn4=156.06 and ωn5=244.91.
Motion profile is spatial form, and the motion profile at each position is all to have 5 cusps, this is that have with response frequency
It closes, the ratio between response frequency is ω4/ω5=0.653, it is approximately equal to 2/3, the number of cusp is equal to the sum of molecule and denominator of ratio.
The parameter resonance response characteristic after fluid conveying pipe buckling is studied below, and dimensionless fluid velocity is u0=3.5, arteries and veins
Dynamic amplitude is μ=0.2.
Referring to shown in Figure 47, fluid flow rate u0=3.5 be more than critical flow velocity m, and theoretical single order buckling amplitude isThe buckling amplitude of numerical solution is 0.004, and the two is very close.
Here it can be observed that after buckling fluid conveying pipe respond " self-locking " phenomenon, in specific frequency separation, pipe conveying fluid
Subharmonic resonance or combination resonance phenomenon can occur for road.There is similar " quench in the parameter resonance frequency zones more than critical flow velocity
Go out " phenomenon, i.e., there is the jump parameter resonance responded response on the boundary of resonance and off-resonance and is dropped under buckling amplitude suddenly.
Ripple frequency is ωp=14 response: the Fourier that Figure 48 and Figure 49 steady-state response, Figure 50 and Figure 51 are in response to
Transformation, Figure 52 and Figure 53 are phasor and Poincare section, and Figure 54 is motion profile, and Figure 55 and Figure 56 are two transverse direction snapshots
Figure.
Referring to shown in Figure 48 to Figure 56, ripple frequency ωpWhen=14, excitation has played the first step mode period vibration of system
It is dynamic, and motion profile is in one plane, the tilt angle of plane of movement is related with primary condition.
Ripple frequency is ωp=395 response: Figure 57 and Figure 58 is steady-state response, in Fu that Figure 59 and Figure 60 are in response to
Leaf transformation, Figure 61 and Figure 62 are phasor and Poincare section, and Figure 63 is motion profile, and Figure 64 and Figure 65 are that two transverse directions are fast
According to figure.
Referring to shown in Figure 57 to Figure 65, ripple frequency ωpWhen=395, excitation has played fourth order and the 5th rank of system
Modal vibration, ripple frequency are approximately equal to the sum of response frequency ωp≈ω4+ω5≈ 158.3+236.2, response frequency are approximately equal to
Fourth order and the 5th rank natural frequency ωn4=153 and ωn5=242.Two transverse direction movement snapshot plotting display pipes movements
It is mutually converted between the fourth order vibration shape and the 5th first order mode.Motion profile is steric configuration, the motion profile at each position
All be to have 5 cusps, this be it is related with response frequency, the ratio between response frequency be ω4/ω5=0.67, it is approximately equal to 2/3, cusp
Number is equal to the sum of molecule and denominator of ratio.
Further invention also provides a kind of aviation pipeline three-dimensional flow and consolidate coupling parameter resonance response Predicting Performance Characteristics device, this
Device may include a lower module
Model building module, the three-dimensional flow for establishing aviation fluid conveying pipe consolidate Coupling Dynamic Model；
Establishing equation module consolidates Coupled Dynamics matrix for establishing aviation pipeline three-dimensional flow according to the kinetic model
Equation；
Equation conversion module, for obtaining finite dimensioned equation according to dynamical matrix equation；
Computing module, the parameter resonance for obtaining the aviation fluid conveying pipe according to the finite dimensioned equation respond
Characteristic
Above-mentioned described feature, structure or characteristic can be incorporated in one or more embodiment party in any suitable manner
In formula, if possible, it is characterized in discussed in each embodiment interchangeable.In the above description, it provides many specific thin
Section fully understands embodiments of the present invention to provide.It will be appreciated, however, by one skilled in the art that this can be practiced
The technical solution of invention, or can be using other methods, component, material without one or more in the specific detail
Material etc..In other cases, known features, material or operation are not shown in detail or describe to avoid each side of the invention is obscured
Face.
The term of " about " " about " is used to be generally represented within the 20% of a given value or range in this specification, preferably
It is within 10%, and is more preferably within 5%.Given quantity is quantity about herein, implies that the feelings in not certain illustrated
Under condition, " about " " about " " substantially " meaning of " general " can be still implied.
Although the term of relativity, such as "upper" "lower" is used to describe a component of icon for another in this specification
The relativeness of one component, but these terms are in this manual merely for convenient, for example, with reference to the accompanying drawings described in show
The direction of example.It is appreciated that, if making it turn upside down the device overturning of icon, the component described in "upper" will
As the component in "lower".Term of other relativities, such as "high" " low " "top" "bottom" " preceding " " rear " " left side " " right side " etc. are also made
With similar meaning.When certain structure is at other structures "upper", it is possible to refer to that certain structural integrity is formed in other structures, or
Refer to that certain structure is " direct " to be arranged in other structures, or refers to that certain structure is arranged in other structures by the way that another structure is " indirect ".
In this specification, term "one", " one ", "the", " described " and "at least one" indicating there are one or
Multiple element/component parts/etc.；Term "comprising", " comprising " and " having " are to indicate the open meaning being included
And refer to the element in addition to listing/component part/also may be present other than waiting other element/component part/etc.；Term " the
One ", " second " and " third " etc. only use as label, are not the quantity limitations to its object.
It should be appreciated that the present invention is not limited in its application to the detailed construction and arrangement of the component of this specification proposition
Mode.The present invention can have other embodiments, and can realize and execute in many ways.Aforesaid deformation form and
Modification is fallen within the scope of the present invention.It should be appreciated that this disclosure and the present invention of restriction extend in text
And/or it is mentioned in attached drawing or all alternative combinations of two or more apparent independent features.It is all these different
Combination constitutes multiple alternative aspects of the invention.Embodiment described in this specification illustrates to become known for realizing the present invention
Best mode, and will enable those skilled in the art using the present invention.
Claims (10)
1. a kind of aviation pipeline three-dimensional flow consolidates coupling parameter resonance response characteristic prediction method characterized by comprising
Consolidate Coupling Dynamic Model according to the three-dimensional flow that relevant parameter establishes aviation fluid conveying pipe；
Consolidate the Coupled Dynamics equation of motion according to the three-dimensional flow that the kinetic model establishes the aviation fluid conveying pipe；
Finite dimensioned equation is obtained according to the dynamical motion equation；
Parameter resonance of the aviation fluid conveying pipe within the scope of default ripple frequency is obtained according to the finite dimensioned equation
Response characteristic.
2. aviation pipeline three-dimensional flow according to claim 1 consolidates coupling parameter resonance response characteristic prediction method, feature
It is, the relevant parameter includes:
Geometric parameter, internal diameter, outer diameter and total length including the fluid conveying pipe；
Physical parameter, elasticity modulus and density including the fluid conveying pipe；
Duty parameter, including the mean flow rate of fluid, pulsation amplitude and ripple frequency in the fluid conveying pipe.
3. aviation pipeline three-dimensional flow according to claim 1 consolidates coupling parameter resonance response characteristic prediction method, feature
It is, finite dimensioned equation is obtained according to the dynamical motion equation, comprising:
The finite dimensioned equation is obtained using the golden method of gal the Liao Dynasty according to the dynamical motion equation.
4. aviation pipeline three-dimensional flow according to claim 1 consolidates coupling parameter resonance response characteristic prediction method, feature
It is, the parameter resonance response characteristic of the aviation fluid conveying pipe is obtained according to the finite dimensioned equation, comprising:
The parameter resonance response characteristic of fluid conveying pipe is calculated using Runge Kutta method according to the finite dimensioned equation.
5. aviation pipeline three-dimensional flow according to claim 4 consolidates coupling parameter resonance response characteristic prediction method, feature
It is, the parameter resonance response characteristic includes:
Aviation fluid conveying pipe fluid velocity inside three-dimensional flow of aviation fluid conveying pipe when being lower than critical flow velocity consolidates coupling parameter resonance
The three-dimensional flow of aviation fluid conveying pipe consolidates coupling parameter when response characteristic and aviation fluid conveying pipe fluid velocity inside are more than critical flow velocity
Resonance response characteristic.
6. aviation pipeline three-dimensional flow according to claim 4 consolidates coupling parameter resonance response characteristic prediction method, feature
It is, the parameter resonance response characteristic is aviation fluid conveying pipe corresponding to different ripple frequencies within the scope of default ripple frequency
Vibration amplitude.
7. aviation pipeline three-dimensional flow according to claim 1 consolidates coupling parameter resonance response characteristic prediction method, feature
It is, the aviation pipeline three-dimensional flow consolidates coupling parameter resonance response characteristic prediction method further include:
Eigenvalues analysis is carried out to the dynamical motion equation and obtains intrinsic frequency.
8. aviation pipeline three-dimensional flow according to claim 7 consolidates coupling parameter resonance response characteristic prediction method, feature
It is, the intrinsic frequency includes:
The intrinsic frequency and aviation pipe conveying fluid of aviation fluid conveying pipe fluid velocity inside aviation fluid conveying pipe when being lower than critical flow velocity
The intrinsic frequency of aviation fluid conveying pipe when road fluid velocity inside is more than critical flow velocity.
9. aviation pipeline three-dimensional flow according to claim 8 consolidates coupling parameter resonance response characteristic prediction method, feature
It is, the accuracy of the parameter resonance response characteristic is detected according to the intrinsic frequency.
10. the prediction meanss that a kind of aviation pipeline three-dimensional flow consolidates coupling parameter resonance response characteristic characterized by comprising
Model building module, the three-dimensional flow for establishing aviation fluid conveying pipe consolidate Coupling Dynamic Model；
Establishing equation module consolidates Coupled Dynamics matrix side for establishing aviation pipeline three-dimensional flow according to the kinetic model
Journey；
Equation conversion module, for obtaining finite dimensioned equation according to dynamical matrix equation；
Computing module, the parameter resonance for obtaining the aviation fluid conveying pipe according to the finite dimensioned equation respond special
Property.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910001794.6A CN109635500B (en) | 2019-01-02 | 2019-01-02 | Aviation pipeline three-dimensional fluid-solid coupling parameter resonance response characteristic prediction method and device |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910001794.6A CN109635500B (en) | 2019-01-02 | 2019-01-02 | Aviation pipeline three-dimensional fluid-solid coupling parameter resonance response characteristic prediction method and device |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN109635500A true CN109635500A (en) | 2019-04-16 |
| CN109635500B CN109635500B (en) | 2024-02-02 |
Family
ID=66056669
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910001794.6A Active CN109635500B (en) | 2019-01-02 | 2019-01-02 | Aviation pipeline three-dimensional fluid-solid coupling parameter resonance response characteristic prediction method and device |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109635500B (en) |
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111967204A (en) * | 2020-08-12 | 2020-11-20 | 西北工业大学 | Anti-resonance reliability and sensitivity analysis method for flow transmission pipeline |
| CN112100892A (en) * | 2020-09-18 | 2020-12-18 | 哈尔滨工业大学（威海） | Prediction method for natural frequency of in-band flow flexible pipeline under different boundary conditions |
| CN114459714A (en) * | 2022-02-17 | 2022-05-10 | 中国核电工程有限公司 | Method capable of reducing pipeline flow-induced vibration, pipeline and verification device |
| CN114818292A (en) * | 2022-04-15 | 2022-07-29 | 中国核动力研究设计院 | Energy band structure analysis method for phononic crystal current-carrying pipeline |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20100139417A1 (en) * | 2008-12-02 | 2010-06-10 | Krohne Messtechnik Gmbh | Method for operating a resonance-measuring system and a resonance-measuring system |
| CN105205258A (en) * | 2015-09-22 | 2015-12-30 | 中国核动力研究设计院 | Method for analyzing vortex shedding induced vibration of heat exchanger heat transfer pipe |
| CN107063295A (en) * | 2016-10-26 | 2017-08-18 | 中国矿业大学(北京) | A kind of method for analyzing stability of resonant mode gyro |
| CN107688703A (en) * | 2017-08-25 | 2018-02-13 | 西北工业大学 | The length ratio design method of the soft or hard pipeline that connects |
| CN108763628A (en) * | 2018-04-13 | 2018-11-06 | 西北工业大学 | The design method and device of multispan functionally gradient fluid conveying pipe |
| CN109063305A (en) * | 2018-07-25 | 2018-12-21 | 中南大学 | Defeated stream straight pipeline Vibration Absorption Designing method under random vibration environment |
-
2019
- 2019-01-02 CN CN201910001794.6A patent/CN109635500B/en active Active
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20100139417A1 (en) * | 2008-12-02 | 2010-06-10 | Krohne Messtechnik Gmbh | Method for operating a resonance-measuring system and a resonance-measuring system |
| CN105205258A (en) * | 2015-09-22 | 2015-12-30 | 中国核动力研究设计院 | Method for analyzing vortex shedding induced vibration of heat exchanger heat transfer pipe |
| CN107063295A (en) * | 2016-10-26 | 2017-08-18 | 中国矿业大学(北京) | A kind of method for analyzing stability of resonant mode gyro |
| CN107688703A (en) * | 2017-08-25 | 2018-02-13 | 西北工业大学 | The length ratio design method of the soft or hard pipeline that connects |
| CN108763628A (en) * | 2018-04-13 | 2018-11-06 | 西北工业大学 | The design method and device of multispan functionally gradient fluid conveying pipe |
| CN109063305A (en) * | 2018-07-25 | 2018-12-21 | 中南大学 | Defeated stream straight pipeline Vibration Absorption Designing method under random vibration environment |
Non-Patent Citations (3)
| Title |
|---|
| L. WANG 等: "Three-dimensional dynamics of supported pipes conveying fluid", 《ACTA MECHANICA SINICA》, vol. 33, no. 6, 19 October 2017 (2017-10-19), pages 1065 - 1074, XP036376835, DOI: 10.1007/s10409-017-0718-z * |
| YUN-DONG LI 等: "Nonlinear vibration of slightly curved pipe with conveying pulsating fluid", 《NONLINEAR DYNAMICS》, vol. 88, 23 February 2017 (2017-02-23), pages 2513 - 2529, XP036247444, DOI: 10.1007/s11071-017-3393-5 * |
| 王乙坤 等: "分布式运动约束下悬臂输液管的参数共振研究", 《HTTP://KNS.CNKI.NET/KCMS/DETAIL/11.2062.O3.20181128.1639.004.HTML》, 28 November 2018 (2018-11-28), pages 1 - 12 * |
Cited By (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111967204A (en) * | 2020-08-12 | 2020-11-20 | 西北工业大学 | Anti-resonance reliability and sensitivity analysis method for flow transmission pipeline |
| CN112100892A (en) * | 2020-09-18 | 2020-12-18 | 哈尔滨工业大学（威海） | Prediction method for natural frequency of in-band flow flexible pipeline under different boundary conditions |
| CN112100892B (en) * | 2020-09-18 | 2023-08-18 | 哈尔滨工业大学（威海） | Method for predicting natural frequency of in-band flow flexible pipeline under different boundary conditions |
| CN114459714A (en) * | 2022-02-17 | 2022-05-10 | 中国核电工程有限公司 | Method capable of reducing pipeline flow-induced vibration, pipeline and verification device |
| CN114818292A (en) * | 2022-04-15 | 2022-07-29 | 中国核动力研究设计院 | Energy band structure analysis method for phononic crystal current-carrying pipeline |
| CN114818292B (en) * | 2022-04-15 | 2023-08-25 | 中国核动力研究设计院 | Phonon crystal current carrying pipeline energy band structure analysis method |
Also Published As
| Publication number | Publication date |
|---|---|
| CN109635500B (en) | 2024-02-02 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN109635500A (en) | Aviation pipeline three-dimensional flow consolidates coupling parameter resonance response characteristic prediction method and device | |
| Kolahchi et al. | Visco-nonlocal-refined Zigzag theories for dynamic buckling of laminated nanoplates using differential cubature-Bolotin methods | |
| Wang et al. | Free vibration analysis of metal foam core sandwich beams on elastic foundation using Chebyshev collocation method | |
| Aydogdu | Vibration analysis of cross-ply laminated beams with general boundary conditions by Ritz method | |
| Ahmadian et al. | Modelling and updating of large surface-to-surface joints in the AWE-MACE structure | |
| Li et al. | Aeroelastic flutter analysis for 2D Kirchhoff and Mindlin panels with different boundary conditions in supersonic airflow | |
| Asgari et al. | Aeroelastic characteristics of magneto-rheological fluid sandwich beams in supersonic airflow | |
| Arani et al. | Nonlinear vibration of embedded smart composite microtube conveying fluid based on modified couple stress theory | |
| Almitani et al. | Influence of the perforation configuration on dynamic behaviors of multilayered beam structure | |
| Abbasnejad et al. | Stability analysis of a piezoelectrically actuated micro-pipe conveying fluid | |
| Zou et al. | Fluid-induced vibration of composite natural gas pipelines | |
| Xia et al. | Application of particle damping for vibration attenuation in brake drum | |
| Kuo et al. | An efficient model for the dynamic behaviour of a single pile in viscoelastic soil | |
| Dai et al. | Dynamics and stability of magnetically actuated pipes conveying fluid | |
| CN107169220A (en) | A kind of passive vibration control device design and optimization method for pipe conveying fluid vibration control | |
| Esfahanian et al. | Fluid-Structure Interaction in microchannel using Lattice Boltzmann method and size-dependent beam element | |
| Li et al. | Dynamic response of pipe conveying fluid with lateral moving supports | |
| Sánchez-Merino et al. | Simplified longitudinal seismic response of tunnels linings subjected to surface waves | |
| Qu et al. | Stress response of the hydraulic composite pipe subjected to random vibration | |
| Ma et al. | A planar shock isolation system with high-static-low-dynamic-stiffness characteristic based on cables | |
| Wang et al. | Vibration of slender structures subjected to axial flow or axially towed in quiescent fluid | |
| Chakka et al. | Modeling and reduction of shocks on electronic components within a projectile | |
| Fu et al. | Boundary element based discontinuous deformation analysis | |
| Bakhshandeh et al. | Benchmark solution for free vibration analysis of transversely isotropic thick rectangular plates | |
| CN102353516A (en) | Method for determining time-domain pulsating drag force of deepwater riser |
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 | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |