CN113919104A - Method and system for acquiring nonlinear dynamic response of rotary drum - Google Patents
Method and system for acquiring nonlinear dynamic response of rotary drum Download PDFInfo
- Publication number
- CN113919104A CN113919104A CN202111513866.9A CN202111513866A CN113919104A CN 113919104 A CN113919104 A CN 113919104A CN 202111513866 A CN202111513866 A CN 202111513866A CN 113919104 A CN113919104 A CN 113919104A
- Authority
- CN
- China
- Prior art keywords
- rotary drum
- equation set
- nonlinear
- cylindrical shell
- representing
- 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
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/17—Mechanical parametric or variational design
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/11—Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
- G06F17/12—Simultaneous equations, e.g. systems of linear equations
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/11—Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
- G06F17/13—Differential equations
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Data Mining & Analysis (AREA)
- General Engineering & Computer Science (AREA)
- Geometry (AREA)
- Operations Research (AREA)
- Software Systems (AREA)
- Databases & Information Systems (AREA)
- Algebra (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Feedback Control In General (AREA)
Abstract
The invention provides a method and a system for acquiring nonlinear dynamic response of a rotary drum, wherein the method comprises the following steps: preprocessing a target rotary drum barrel to obtain parameter information of the target rotary drum barrel; respectively acquiring strain energy, kinetic energy and external force acting corresponding to the target rotary drum based on the parameter information; on the basis of strain energy, kinetic energy and external force acting, a nonlinear forced vibration partial differential control equation set corresponding to the target rotary drum is established; carrying out weighted integral processing on the control equation set, and converting the control equation set into an ordinary differential equation set with multiple coupling degrees of freedom; and carrying out numerical solution on the ordinary differential equation set to obtain a nonlinear amplitude-frequency curve and a force amplitude response curve of the target rotary drum. The invention can effectively solve the problem of nonlinear vibration response when the rotary drum of the engine is excited by multiple harmonics, and has the characteristics of high precision, strong universality and the like.
Description
Technical Field
The invention relates to the technical field of mechanical dynamics, in particular to a nonlinear dynamic response obtaining method and system of a rotary drum.
Background
At present, for an engine with a complex structure, a rotary drum in a rotor system is required to have enough structural strength and rigidity to bear axial thrust, internal hydraulic pressure and strong vibration; and it is required to be as light as possible in weight to improve the thrust-weight ratio of the engine, so that the thin-walled cylindrical shell is widely used in various engines with its excellent structure. However, the rotating drum is affected by manufacturing and installation errors or thermal deformation in the operation process, so that unbalance exists inevitably, the rotating drum is easily excited by multiple harmonics of complex surrounding environment, beat vibration, resonance or more complex phenomena occur, the dynamic behavior of the rotating drum is more complex, if the dynamic response of the engine cannot be accurately obtained, the normal use of a rotor system can be greatly damaged, the whole engine can be stopped in serious cases, and huge economic loss and severe social influence are caused.
Therefore, the research on the nonlinear dynamic response of the rotary thin-wall cylindrical shell structure plays an important role in ensuring the structural vibration characteristic and stable operation after the rotor system is subjected to complex multi-harmonic excitation.
At present, the dynamic research scheme of the rotary thin-wall cylindrical shell under multi-harmonic excitation is less, and most of the rotary thin-wall cylindrical shell adopts a finite element method and has no universality; in addition, when the structural model is too complex or nonlinear vibration characteristics occur, the finite element method cannot accurately acquire amplitude deformation, and the limitation is also large. In addition, some traditional analytic or numerical methods simplify the research process of dynamic response, only consider the linear vibration generated by the system, and are incapable of accurately predicting the system prediction error of large-amplitude vibration, especially the complex phenomena such as nonlinear beat vibration and resonance. For example, such errors are unacceptable for precision structures such as aircraft engine rotor systems.
Therefore, the research and establishment of the dynamic response analysis method of the rotary thin-wall cylindrical shell, which has universality and high precision and can be used in a complex multi-harmonic excitation environment, have very important significance.
Disclosure of Invention
In view of the above problems, an object of the present invention is to provide a method and a system for obtaining nonlinear dynamic response of a rotating drum, so as to solve the problems of low precision, no universality and the like of the existing method for studying vibration characteristics of a rotating thin-walled cylindrical shell excited by multiple harmonics.
The invention provides a nonlinear dynamic response obtaining method of a rotary drum, which comprises the following steps: preprocessing a target rotary drum barrel to obtain parameter information of the target rotary drum barrel; respectively acquiring strain energy, kinetic energy and external force acting corresponding to the target rotary drum based on the parameter information; on the basis of strain energy, kinetic energy and external force acting, a nonlinear forced vibration partial differential control equation set corresponding to the target rotary drum is established; carrying out weighted integral processing on the control equation set, and converting the control equation set into an ordinary differential equation set with multiple coupling degrees of freedom; and carrying out numerical solution on the ordinary differential equation set to obtain a nonlinear amplitude-frequency curve and a force amplitude response curve of the target rotary drum.
According to another aspect of the present invention, there is provided a non-linear dynamic response acquiring system of a rotary drum, including: the parameter information acquisition unit is used for preprocessing the target rotary drum barrel and acquiring the parameter information of the target rotary drum barrel; an energy and work obtaining unit for respectively obtaining strain energy, kinetic energy and external force acting corresponding to the target rotary drum based on the parameter information; the control equation set acquisition unit is used for doing work based on strain energy, kinetic energy and external force and establishing a nonlinear forced vibration partial differential control equation set corresponding to the target rotary drum; the system comprises an ordinary differential equation set conversion unit, a control equation set calculation unit and a control equation set calculation unit, wherein the ordinary differential equation set conversion unit is used for carrying out weighted integral processing on the control equation set and converting the control equation set into an ordinary differential equation set with multiple coupling degrees of freedom; and the response curve acquisition unit is used for carrying out numerical solution on the ordinary differential equation set to acquire a nonlinear amplitude-frequency curve and a force amplitude response curve of the target rotary drum.
By utilizing the nonlinear dynamic response obtaining method and the system of the rotary drum, the parameter information of the target rotary drum is obtained by preprocessing the target rotary drum, then the strain energy, the kinetic energy and the external force which correspond to the target rotary drum are respectively obtained to do work, a nonlinear forced vibration partial differential control equation set corresponding to the target rotary drum is established, further the weighted integral processing is carried out on the control equation set, the control equation set is converted into a coupled multi-degree-of-freedom ordinary differential equation set, the ordinary differential equation set is numerically solved, a nonlinear amplitude-frequency curve and a force amplitude response curve of the target rotary drum are obtained, the method and the system can be used for solving the complicated nonlinear dynamic response of the rotary drum, the calculation is convenient, the readability is strong, and the accuracy is high.
To the accomplishment of the foregoing and related ends, one or more aspects of the invention comprise the features hereinafter fully described. The following description and the annexed drawings set forth in detail certain illustrative aspects of the invention. These aspects are indicative, however, of but a few of the various ways in which the principles of the invention may be employed. Further, the present invention is intended to include all such aspects and their equivalents.
Drawings
Other objects and results of the present invention will become more apparent and more readily appreciated as the same becomes better understood by reference to the following description taken in conjunction with the accompanying drawings. In the drawings:
FIG. 1 is a flow chart of a method of obtaining a nonlinear dynamic response of a rotating drum according to an embodiment of the present invention;
FIG. 2 is a block diagram of an equivalent rotating thin-walled cylindrical shell model according to an embodiment of the invention;
fig. 3 is a detailed flowchart of a nonlinear dynamic response acquisition method of a rotary drum according to an embodiment of the present invention;
fig. 4 shows a block schematic diagram of a non-linear dynamic response acquisition system for a rotating drum according to an embodiment of the present invention.
The same reference numbers in all figures indicate similar or corresponding features or functions.
Detailed Description
In the following description, for purposes of explanation, numerous specific details are set forth in order to provide a thorough understanding of one or more embodiments. It may be evident, however, that such embodiment(s) may be practiced without these specific details. In other instances, well-known structures and devices are shown in block diagram form in order to facilitate describing one or more embodiments.
In the description of the present invention, it is to be understood that the terms "central," "longitudinal," "lateral," "length," "width," "thickness," "upper," "lower," "front," "rear," "left," "right," "vertical," "horizontal," "top," "bottom," "inner," "outer," "clockwise," "counterclockwise," "axial," "radial," "circumferential," and the like are used in the orientations and positional relationships indicated in the drawings for convenience in describing the invention and to simplify the description, and are not intended to indicate or imply that the referenced devices or elements must have a particular orientation, be constructed and operated in a particular orientation, and are therefore not to be considered limiting of the invention.
To describe in detail the present invention's method and system for obtaining nonlinear dynamic response of a rotating drum, a detailed description of embodiments of the present invention will be given below with reference to the accompanying drawings.
Fig. 1 shows a schematic flow of a non-linear dynamic response acquisition method of a rotating drum according to an embodiment of the present invention.
As shown in fig. 1, the method for acquiring nonlinear dynamic response of a rotating drum according to an embodiment of the present invention mainly includes the following steps:
s110: and preprocessing the target rotary drum barrel to obtain the parameter information of the target rotary drum barrel.
In this step, the target rotating drum can be equated to a rotating thin-walled cylindrical shell model; then, acquiring parameter information of the rotary thin-wall cylindrical shell model as parameter information of a target rotary drum barrel; wherein the parameter information comprises a geometric parameter and a physical parameter.
Specifically, FIG. 2 illustrates a schematic structure of an equivalent rotating thin-walled cylindrical shell model according to an embodiment of the present invention, as shown in FIG. 2, for rotating a target drum, e.g., a aeronautical instrumentThe drum-type rotor structure of the air engine is equivalent to a rotary thin-wall cylindrical shell model, and then the radius of the rotary thin-wall cylindrical shell model is measuredRIs long and longLAnd thicknesshModulus of elasticityEPoisson ratiovAnd mass densityρThe parameter information may be used as parameter information of the target rotary drum.
S120: and respectively acquiring strain energy, kinetic energy and external force acting corresponding to the target rotary drum based on the parameter information.
Wherein the process of acquiring strain energy corresponding to the target rotating drum based on the parameter information may further include the steps of:
s121: establishing a relation model between nonlinear strain and displacement of the rotary thin-wall cylindrical shell model based on the parameter information;
specifically, the expression of the relational model between the nonlinear strain and the displacement is:
wherein the content of the first and second substances,、andrepresenting the strain of any point on the rotating thin-wall cylindrical shell model;、andrepresenting the mid-plane strain of the rotating thin-wall cylindrical shell model;、andrespectively representing the mid-plane curvature and twist in the thin-walled cylindrical shell model of revolution.
Wherein, the expression of the middle surface strain is as follows:
the expressions for mid-plane curvature and twist are:
wherein the content of the first and second substances,u,v,wrespectively a middle surface edge of a rotary thin-wall cylindrical shell modelx、、zDisplacement in the axial direction.
S122: and acquiring a stress and strain relation model under the thermal effect based on the relation model.
Specifically, the expression of the stress-strain relationship model is:
wherein the content of the first and second substances,andrepresenting the coefficient of thermal expansion in the corresponding direction;、andrepresents the strain of any point on the rotating thin-wall cylindrical shell model, deltaTExpressed as a change in temperature relative to the initial state;C ij the component of the elastic stiffness matrix is represented by the following specific expression:
wherein the content of the first and second substances,Eit means the modulus of elasticity of the polymer,vrepresenting the poisson ratio.
S123: acquiring strain energy of the rotary thin-wall cylindrical shell model based on a stress and strain relation model; wherein the strain energy comprises elastic strain energy and additional strain energy.
Specifically, the expression formula of the elastic strain energy is as follows:
wherein the content of the first and second substances,srepresenting the infinitesimal area of the rotating thin-walled cylindrical shell model,hrepresenting the thickness of the rotating thin-walled cylindrical shell model,、、the relationship between the stress and the strain is shown,、andrepresenting the strain of any point on the rotating thin-wall cylindrical shell model;
the expression formula of the additional strain energy is:
wherein the content of the first and second substances,srepresenting the infinitesimal area of the rotating thin-walled cylindrical shell model,,which represents the angular velocity of rotation of the rotating body,Rrepresenting the radius of the rotating thin-walled cylindrical shell model,ρthe mass density is expressed in terms of,u,v,wrespectively the middle surface of the rotary thin-wall cylindrical shell model is arranged alongx、、zDisplacement in the axial direction.
Further, the expression of the kinetic energy corresponding to the target rotary drum is:
wherein the content of the first and second substances,Lrepresenting the length of the rotating thin-walled cylindrical shell model,,ρthe mass density is expressed in terms of,which represents the angular velocity of rotation of the rotating body,u,v,wrespectively the middle surface of the rotary thin-wall cylindrical shell model is arranged alongx、、zDisplacement in the axial direction;
the expression of the external force acting corresponding to the target rotary drum is as follows:
wherein the content of the first and second substances,representing the radially distributed force acting on a unit area, the expression:
wherein the content of the first and second substances,fin terms of radial force amplitude (x 0, ) Is the location of the point of applied radial force;representing a dirac function;ω f is the excitation frequency;Jrepresents the total number of harmonics;andrespectively representing different phases not equal to the fundamental harmonicThe magnitude coefficient of the bit.
S130: and establishing a nonlinear forced vibration partial differential control equation set corresponding to the target rotary drum based on the strain energy, the kinetic energy and the external force acting.
Based on the strain energy, the kinetic energy and the external force to do work, the process of establishing a nonlinear forced vibration partial differential control equation set corresponding to the target rotary drum comprises the following steps:
constructing a variation expression based on a Hamilton variation principle, the strain energy, the kinetic energy and the external force acting:
wherein the content of the first and second substances,tthe time is represented by the time of day,the kinetic energy is represented by a representation of the kinetic energy,which is indicative of the elastic strain energy of the material,which represents the additional strain energy, is,the external force is shown to do work,representing a dirac function;
using Hamilton principle variation, namely:then, based on the variational expressions, the information aboutu,v,wVariation of direction to obtain the damping coefficient containing structurec d The system of nonlinear forced vibration partial differential control equations:
wherein the content of the first and second substances,Rrepresenting the radius of the rotating thin-walled cylindrical shell model,hrepresenting the thickness of the rotating thin-walled cylindrical shell model,which represents the angular velocity of rotation of the rotating body,,,the moment of inertia is represented as a function of,which is indicative of an external force,x、andzare all indicative of the coordinate axes,u,v,wrespectively the middle surface of the rotary thin-wall cylindrical shell model is arranged alongx、、zDisplacement in the axial direction.
S140: and performing weighted integral processing on the control equation set, and converting the control equation set into an ordinary differential equation set with multiple coupling degrees of freedom.
This step may include: establishing a displacement function meeting boundary conditions; and then, converting the control equation set into a multi-degree-of-freedom coupled ordinary differential equation set by utilizing a multi-mode Galerkin method and a displacement function.
Specifically, by establishing that the boundary condition(s) is satisfied,,,) The displacement function of (2) simplifies the nonlinear system into finite dimensions, and the expression is as follows:
wherein the content of the first and second substances,u m,n (t),v m,n (t)and w m,n (t) Representing a displacement amplitude component;M 1,N 1representing a modal truncation coefficient;min the order of the axial half-wave number,nis the circumferential wave number; subscriptcIndicating driving mode, subscriptsRepresenting the accompanying modality. In nonlinear vibration, these two modes interact through one-to-one internal resonance.
Then, a multi-mode Galerkin method is utilized to convert the control equation set into an ordinary differential equation set coupled with multiple degrees of freedom, namely, the original equation set is weighted by proper functions in sequence and is integrated on the middle surface of the rotary thin-wall cylindrical shell model, wherein the weighting functions are as follows:
then, the system of ordinary differential equations is expressed in matrix form as:
wherein the content of the first and second substances,m, C and K are denoted as stiffness matrix, damping matrix and mass matrix, respectively; g represents a gyro matrix generated by rotation;represent other non-linear portions; f represents the vector of the external force.
S150: and carrying out numerical solution on the ordinary differential equation set to obtain a nonlinear amplitude-frequency curve and a force amplitude response curve of the target rotary drum.
Wherein, because the system is non-autonomous bivalent and can not be directly solved, virtual variables can be introducedAs the speed, a first-order non-autonomous system is obtained as follows:
it can be seen that the above system has been reduced to first order, but is still not autonomous, and therefore can adopt the Hopf paradigm:
without loss of generality, the first harmonic excitation and the second harmonic excitation are chosen here for elucidation, namely:let (1) , ) = (1,0), = cos(ωt)、 (t) = sin(ωt) (ii) a Then, let cos (2) in the code of the MATLAB program definition modelωt)= -Andthus, it is possible to obtain:
neglecting the above system of ordinary differential equationsuAndvthe inertia term and the rotation term of the direction can be obtained only in relation to the radial directionwThe ordinary differential equation of (2); then, on the basis of MATLAB software, a numerical algorithm is written by adopting a quasi-arc length continuation method to solve the nonlinear coupling ordinary differential equation set, the result is compared with the experimental result, and if the verification is successful, an amplitude-frequency response curve and a force-amplitude response curve can be obtained. If the verification fails, the displacement function is reconstructed, and the coupled multi-degree-of-freedom ordinary differential equation set is solved again. The above operation is then carried out.
As a specific example, fig. 3 shows a detailed flow of a nonlinear dynamic response acquisition method of a rotary drum according to an embodiment of the present invention.
As shown in fig. 3, before processing a target rotary drum or a rotary drum, the target rotary drum or the rotary drum is equivalent to a rotary thin-wall cylindrical shell, and then strain energy, kinetic energy and external force acting under coriolis force, centrifugal force and thermal effect are obtained based on multi-harmonic excitation; then constructing a nonlinear forced vibration partial differential equation set on the basis of the Hamilton principle on the basis of acting of strain energy, kinetic energy and external force; further constructing an ordinary differential equation set coupling multiple degrees of freedom based on a multi-mode Galerkin method; further based on a quasi-arc length continuation method, obtaining a dynamic characteristic calculation result of the rotary thin-wall cylindrical shell, verifying the dynamic characteristic calculation result, and outputting a nonlinear amplitude-frequency and force-amplitude response curve corresponding to the rotary drum if the verification is passed; otherwise, solving the ordinary differential equation set of the coupling multiple degrees of freedom again.
The nonlinear dynamic response obtaining method of the rotary drum barrel can be used for solving the complex nonlinear dynamic response of the rotary drum barrel, such as the complex problems of nonlinear beat vibration, nonlinear multiple internal resonance and the like which are often found by engineering, and has the following advantages:
1. by adopting the modified Donnell nonlinear shell theory, the defect that the traditional Donnell nonlinear shell theory is inaccurate in wavelet number can be overcome;
2. the nonlinear dynamic response analysis of different rotary drum barrels under multi-harmonic excitation can be met only by controlling the geometric physical parameters of the rotary cylindrical shell, external excitation and other parameters, and modeling and any modification on a calculation program are not needed;
3. and the multi-mode Galerkin method and the quasi-arc length continuation method are adopted to carry out frequency solution on the nonlinear forced vibration of the rotating thin-wall cylindrical shell, so that the calculation is convenient and the readability is strong.
Corresponding to the nonlinear dynamic response acquisition method of the rotary drum, the invention also provides a nonlinear dynamic response acquisition system of the rotary drum.
In particular, fig. 4 shows a schematic logic of a non-linear dynamic response acquisition system for a rotating drum according to an embodiment of the present invention.
As shown in fig. 4, the non-linear dynamic response acquisition system 100 for a rotating drum of the present invention comprises:
a parameter information obtaining unit 101, configured to perform preprocessing on a target rotary drum to obtain parameter information of the target rotary drum;
an energy and work obtaining unit 102 configured to obtain strain energy, kinetic energy, and external force work corresponding to the target rotary drum, respectively, based on the parameter information;
a control equation set obtaining unit 103, configured to apply work based on strain energy, kinetic energy, and external force, and establish a nonlinear forced vibration partial differential control equation set corresponding to the target rotary drum;
an ordinary differential equation set conversion unit 104, configured to perform weighted integration processing on the control equation set, and convert the control equation set into an ordinary differential equation set with multiple degrees of freedom in coupling;
and the response curve acquiring unit 105 is configured to perform numerical solution on the system of ordinary differential equations to acquire a nonlinear amplitude-frequency curve and a force amplitude response curve of the target rotary drum.
By utilizing the nonlinear dynamic response obtaining method and the nonlinear dynamic response obtaining system of the rotary drum, after the target rotary drum is equivalent to a rotary thin-wall cylindrical shell structure, the nonlinear control equation set of the structure under the multi-harmonic excitation is obtained by utilizing the improved Donnell nonlinear thin-shell theory and the Hamilton variation principle, then further converting the nonlinear control equation set into a coupled multi-degree-of-freedom ordinary differential equation set by a multi-mode Galerkin method, and based on MATLAB software, the ordinary differential equation set is numerically solved by adopting a quasi-arc length continuation method to obtain nonlinear dynamic response of the rotary thin-wall cylindrical shell under the multi-harmonic excitation, the problem of nonlinear vibration response when the rotary drum of the aero-engine is subjected to multi-harmonic excitation can be effectively solved, the method has the advantages of being high in precision, strong in universality and the like, and has great significance in the aspects of safety protection of the aero-engine, structural resonance avoidance and the like.
The non-linear dynamic response acquisition method and system of a rotating drum according to the present invention are described above by way of example with reference to the accompanying drawings. However, it should be understood by those skilled in the art that various modifications can be made to the non-linear dynamic response obtaining method and system for a rotating drum proposed by the present invention without departing from the scope of the present invention. Therefore, the scope of the present invention should be determined by the contents of the appended claims.
Claims (10)
1. A method of obtaining a nonlinear dynamic response of a rotating drum, comprising:
preprocessing a target rotary drum barrel to obtain parameter information of the target rotary drum barrel;
respectively acquiring strain energy, kinetic energy and external force acting corresponding to the target rotary drum based on the parameter information;
based on the strain energy, the kinetic energy and the external force acting, establishing a nonlinear forced vibration partial differential control equation set corresponding to the target rotary drum;
carrying out weighted integral processing on the control equation set, and converting the control equation set into an ordinary differential equation set with multiple coupling degrees of freedom;
and carrying out numerical solution on the ordinary differential equation set to obtain a nonlinear amplitude-frequency curve and a force amplitude response curve of the target rotary drum.
2. The method for acquiring nonlinear dynamic response of a rotating drum according to claim 1, wherein the preprocessing of the target rotating drum and the acquiring of the parameter information of the target rotating drum comprise:
the target rotary drum barrel is equivalent to a rotary thin-wall cylindrical shell model;
acquiring parameter information of the rotating thin-wall cylindrical shell model as parameter information of the target rotating drum; wherein the parameter information comprises a geometric parameter and a physical parameter.
3. The method of acquiring a nonlinear dynamic response of a rotary drum according to claim 2, wherein the process of acquiring the strain energy corresponding to the target rotary drum based on the parameter information includes:
establishing a relation model between nonlinear strain and displacement of the rotary thin-wall cylindrical shell model based on the parameter information;
obtaining a stress and strain relation model under the thermal effect based on the relation model;
acquiring strain energy of the rotary thin-wall cylindrical shell model based on the stress and strain relation model; wherein the strain energy comprises elastic strain energy and additional strain energy.
4. The method of obtaining a nonlinear dynamic response of a rotating drum of claim 3, wherein the expression of the model of the relationship between nonlinear strain and displacement is:
wherein the content of the first and second substances,、andrepresenting the strain of any point on the rotating thin-wall cylindrical shell model;、andrepresenting the mid-plane strain of the rotating thin-walled cylindrical shell model;、andrespectively representing the surface curvature and torsion in the rotating thin-wall cylindrical shell model;
the expression of the mid-plane strain is as follows:
the expressions of the mid-plane curvature and twist are:
5. The method of claim 4, wherein the stress versus strain model is expressed as:
wherein the content of the first and second substances,andrepresenting the coefficient of thermal expansion in the corresponding direction;、andrepresents the strain of any point on the rotating thin-wall cylindrical shell model, deltaTExpressed as a change in temperature relative to the initial state;C ij the component of the elastic stiffness matrix is represented by the following specific expression:
wherein the content of the first and second substances,Eit means the modulus of elasticity of the polymer,vrepresenting the poisson ratio.
6. The method of obtaining nonlinear dynamic response of a rotating drum of claim 5, wherein the elastic strain energy is expressed by the formula:
wherein the content of the first and second substances,srepresenting the infinitesimal area of the rotating thin-walled cylindrical shell model,hrepresenting the thickness of the rotating thin-walled cylindrical shell model,、、the relationship between the stress and the strain is shown,、andrepresenting the strain of any point on the rotating thin-wall cylindrical shell model;
the expression formula of the additional strain energy is:
wherein the content of the first and second substances,srepresenting the infinitesimal area of the rotating thin-walled cylindrical shell model,,which represents the angular velocity of rotation of the rotating body,Rrepresenting the radius of the rotating thin-walled cylindrical shell model,ρthe mass density is expressed in terms of,u,v,wrespectively the middle surface of the rotary thin-wall cylindrical shell model is arranged alongx、、zDisplacement in the axial direction.
7. The method of obtaining a nonlinear dynamic response of a rotating drum of claim 2, wherein the expression of the kinetic energy corresponding to the target rotating drum is:
wherein the content of the first and second substances,Lrepresenting the length of the rotating thin-walled cylindrical shell model,,ρthe mass density is expressed in terms of,which represents the angular velocity of rotation of the rotating body,u,v,wrespectively the middle surface of the rotary thin-wall cylindrical shell model is arranged alongx、、zDisplacement in the axial direction;
the expression of the external force acting corresponding to the target rotary drum is as follows:
wherein the content of the first and second substances,representing the radially distributed force acting on a unit area, the expression:
wherein the content of the first and second substances,fin terms of radial force amplitude (x 0,) Is the location of the point of applied radial force;representing a dirac function;ω f is the excitation frequency;Jrepresents the total number of harmonics;andrespectively representing amplitude coefficients not equal to different phases of the fundamental harmonic.
8. The method of obtaining a nonlinear dynamic response of a rotary drum of claim 3, wherein the process of establishing a system of nonlinear forced vibration partial differential control equations corresponding to the target rotary drum based on the strain energy, the kinetic energy, and the external force to do work comprises:
constructing a variation expression based on a Hamilton variation principle, the strain energy, the kinetic energy and the external force acting:
wherein the content of the first and second substances,tthe time is represented by the time of day,the kinetic energy is represented by a representation of the kinetic energy,which is indicative of the elastic strain energy of the material,which represents the additional strain energy, is,the external force is shown to do work,representing a dirac function;
based on the variational expression, obtaining a damping coefficient containing a structurec d The system of nonlinear forced vibration partial differential control equations:
wherein the content of the first and second substances,Rrepresenting the radius of the rotating thin-walled cylindrical shell model,hrepresenting the thickness of the rotating thin-walled cylindrical shell model,which represents the angular velocity of rotation of the rotating body,,,the moment of inertia is represented as a function of,which is indicative of an external force,x、andzare all indicative of the coordinate axes,u,v,wrespectively the middle surface of the rotary thin-wall cylindrical shell model is arranged alongx、、zDisplacement in the axial direction.
9. The method of obtaining nonlinear dynamic response of a rotating drum according to claim 4, wherein the process of performing weighted integration on the control equation set to convert the control equation set into an ordinary differential equation set coupling multiple degrees of freedom comprises:
establishing a displacement function meeting boundary conditions;
converting the control equation set into a multi-degree-of-freedom coupled ordinary differential equation set by utilizing a multi-mode Galerkin method and the displacement function; wherein the system of control equations is weighted in turn with a suitable weighting function, the weighting function being:
the matrix form of the system of ordinary differential equations is expressed as:
wherein the content of the first and second substances,m, C and K are denoted as stiffness matrix, damping matrix and mass matrix, respectively; g denotes a gyro matrix generated by the rotation,representing other non-linear portions, F representing a vector of external forces;
performing numerical solution on the ordinary differential equation set to obtain a nonlinear amplitude-frequency curve and a force amplitude response curve of the target rotary drum, wherein the numerical solution comprises the following steps:
introducing virtual variablesAs the speed, a first order non-autonomy is obtainedThe system is as follows:
using the Hopf paradigm:
wherein the content of the first and second substances,let (1) , ) = (1,0), = cos(ωt)、 (t) = sin(ωt) (ii) a Then, let cos (2) in the code of the MATLAB program definition modelωt)= -Andtherefore, the deformed ordinary differential equation system can be obtained as follows:
and acquiring a nonlinear amplitude-frequency curve and a force amplitude response curve of the target rotary drum based on the deformed ordinary differential equation set.
10. A non-linear dynamic response acquisition system for a rotating drum, comprising:
the parameter information acquisition unit is used for preprocessing a target rotary drum barrel and acquiring the parameter information of the target rotary drum barrel;
an energy and work obtaining unit, configured to obtain strain energy, kinetic energy, and external force work corresponding to the target rotary drum, respectively, based on the parameter information;
a control equation set obtaining unit, configured to apply work based on the strain energy, the kinetic energy, and the external force, and establish a nonlinear forced vibration partial differential control equation set corresponding to the target rotary drum;
the system comprises an ordinary differential equation set conversion unit, a control equation set calculation unit and a control equation set calculation unit, wherein the ordinary differential equation set conversion unit is used for carrying out weighted integral processing on the control equation set and converting the control equation set into an ordinary differential equation set with multiple coupling degrees of freedom;
and the response curve acquisition unit is used for carrying out numerical solution on the ordinary differential equation set to acquire a nonlinear amplitude-frequency curve and a force amplitude response curve 2 of the target rotary drum.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111513866.9A CN113919104B (en) | 2021-12-13 | 2021-12-13 | Method and system for acquiring nonlinear dynamic response of rotary drum |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111513866.9A CN113919104B (en) | 2021-12-13 | 2021-12-13 | Method and system for acquiring nonlinear dynamic response of rotary drum |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113919104A true CN113919104A (en) | 2022-01-11 |
CN113919104B CN113919104B (en) | 2022-04-26 |
Family
ID=79248634
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202111513866.9A Active CN113919104B (en) | 2021-12-13 | 2021-12-13 | Method and system for acquiring nonlinear dynamic response of rotary drum |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113919104B (en) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115826407A (en) * | 2022-11-29 | 2023-03-21 | 中国航发沈阳发动机研究所 | Control method for reducing rotation inertia excitation of drum shaft |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
GB2069727A (en) * | 1980-02-15 | 1981-08-26 | Kloeckner Humboldt Deutz Ag | Controlling the difference in speed between two rotating components |
US20030019285A1 (en) * | 1998-06-29 | 2003-01-30 | Juranitch James C. | Isolation and loading arrangement for differential gear system under test |
CN109800512A (en) * | 2019-01-23 | 2019-05-24 | 东北大学 | Rotating cylindrical shell-variable cross-section disk-pretwist blade system dynamic modeling method |
CN110442971A (en) * | 2019-08-06 | 2019-11-12 | 东北大学 | A kind of rotating cylindrical shell kinetic characteristics Uncertainty Analysis Method |
CN110457823A (en) * | 2019-08-13 | 2019-11-15 | 大连民族大学 | The MLP method of super-elasticity cylindrical thin shell strong nonlinear vibration |
-
2021
- 2021-12-13 CN CN202111513866.9A patent/CN113919104B/en active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
GB2069727A (en) * | 1980-02-15 | 1981-08-26 | Kloeckner Humboldt Deutz Ag | Controlling the difference in speed between two rotating components |
US20030019285A1 (en) * | 1998-06-29 | 2003-01-30 | Juranitch James C. | Isolation and loading arrangement for differential gear system under test |
CN109800512A (en) * | 2019-01-23 | 2019-05-24 | 东北大学 | Rotating cylindrical shell-variable cross-section disk-pretwist blade system dynamic modeling method |
CN110442971A (en) * | 2019-08-06 | 2019-11-12 | 东北大学 | A kind of rotating cylindrical shell kinetic characteristics Uncertainty Analysis Method |
CN110457823A (en) * | 2019-08-13 | 2019-11-15 | 大连民族大学 | The MLP method of super-elasticity cylindrical thin shell strong nonlinear vibration |
Non-Patent Citations (1)
Title |
---|
陈科文: ""转动圆柱壳非线性动力学建模与数值仿真"", 《中国优秀硕士学位论文全文数据库(电子期刊)基础科学》 * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115826407A (en) * | 2022-11-29 | 2023-03-21 | 中国航发沈阳发动机研究所 | Control method for reducing rotation inertia excitation of drum shaft |
CN115826407B (en) * | 2022-11-29 | 2024-04-09 | 中国航发沈阳发动机研究所 | Control method for reducing drum shaft rotation inertia excitation |
Also Published As
Publication number | Publication date |
---|---|
CN113919104B (en) | 2022-04-26 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Bazilevs et al. | Isogeometric fatigue damage prediction in large-scale composite structures driven by dynamic sensor data | |
Gu et al. | Free vibration of rotating cantilever pre-twisted panel with initial exponential function type geometric imperfection | |
Hashemi et al. | Free vibration analysis of rotating thick plates | |
CN109800512B (en) | Dynamic modeling method for rotating cylindrical shell-variable cross-section disc-pre-twisted blade system | |
Kerboua et al. | Numerical model to analyze the aerodynamic behavior of a combined conical–cylindrical shell | |
Hussain et al. | Vibration characteristics of fluid-filled functionally graded cylindrical material with ring supports | |
CN113919104B (en) | Method and system for acquiring nonlinear dynamic response of rotary drum | |
CN110020460A (en) | It is bolted flanged (FLGD) cylindrical shell structure frequency response function Uncertainty Analysis Method | |
Rawat et al. | Free vibration analysis of thin circular cylindrical shell with closure using finite element method | |
Patil et al. | A scalable time-parallel solution of periodic rotor dynamics in X3D | |
Chai et al. | Frequency response of spinning cylindrical shells with discontinuous boundary conditions: a semi-analytical method | |
CN107291988B (en) | Method for acquiring equivalent excitation force of momentum wheel mounting interface | |
Mavriplis et al. | Recent advances in high-fidelity multidisciplinary adjoint-based optimization with the NSU3D flow solver framework | |
Campagnari et al. | Estimation of axial load in tie-rods using experimental and operational modal analysis | |
Sabri et al. | Finite element method applied to supersonic flutter of circular cylindrical shells | |
Zeng et al. | Dynamic modeling and coupling characteristics of rotating inclined beams with twisted-shape sections | |
Liu et al. | Travelling-modes of a spinning cylindrical shell with elastic boundaries | |
Swami et al. | Vibration analysis of a tire under static loading using flexible ring-based model | |
Pešek et al. | Dry friction interblade damping by 3d fem modelling of bladed disk: Hpc calculations compared with experiment | |
Jiammeepreecha et al. | Nonlinear axisymmetric free vibration analysis of liquid-filled spherical shell with volume constraint | |
Dhar et al. | Transient stress analysis and fatigue life estimation of turbine blades | |
JP4513776B2 (en) | Earthquake response analysis method | |
Seraj | Free vibration and dynamic instability analyses of doubly-tapered rotating laminated composite beams | |
Thresher et al. | Structural dynamic behavior of wind turbines | |
Liu et al. | An efficient dynamic modeling technique for a central tie rod rotor |
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 |