CN112149308A - Method for quickly and automatically adjusting excitation force - Google Patents
Method for quickly and automatically adjusting excitation force Download PDFInfo
- Publication number
- CN112149308A CN112149308A CN202011059206.3A CN202011059206A CN112149308A CN 112149308 A CN112149308 A CN 112149308A CN 202011059206 A CN202011059206 A CN 202011059206A CN 112149308 A CN112149308 A CN 112149308A
- Authority
- CN
- China
- Prior art keywords
- excitation force
- exciting
- force
- excitation
- force vector
- 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/20—Design optimisation, verification or simulation
-
- 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
- G06F30/25—Design optimisation, verification or simulation using particle-based methods
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/06—Multi-objective optimisation, e.g. Pareto optimisation using simulated annealing [SA], ant colony algorithms or genetic algorithms [GA]
-
- 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
- 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)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
Abstract
The invention discloses a method for quickly and automatically adjusting excitation force, which comprises the following steps: firstly, according to the characteristics of a tested structure or by using an impact force hammering method, determining the number of exciting forces and the position of an exciting point by enabling a vibration exciter to be positioned at a point where each interested order mode has enough displacement; step two, carrying out multipoint sine frequency sweep or multipoint random excitation test on the tested structure to obtain a frequency response function matrix and modal frequency of the structure; thirdly, a purity indication function is used as a fitness function in the algorithm, and a group of optimal excitation force amplitudes can be obtained after the particle swarm optimization iteration; and step four, sending the optimal excitation force amplitude to the vibration exciter to enable the vibration exciter to quickly and automatically adjust the excitation force. The invention provides a new method for quickly and automatically adjusting exciting force in a modal test, which uses an optimization algorithm to realize the automatic adjustment of exciting force, the force adjustment process is simple, and the test time is greatly shortened.
Description
Technical Field
The application relates to the technical field of automatic adaptive exciting force in pure modal tests, in particular to a method for quickly and automatically adapting exciting force.
Background
The automatic adjustment of the exciting force in the pure modal test is always a difficult problem, and the traditional Dark method and the like cannot solve the problem of the automatic adjustment of the exciting force in the pure modal test of a large complex structure because the optimal distribution of the exciting force cannot be searched globally.
The method for adjusting the vibration force developed based on the idea that the modal force of the adaptive mode is not zero and the modal force of the non-adaptive mode is zero is to obtain the approximate distribution of the modal shape by a frequency domain modal parameter identification method.
Disclosure of Invention
In order to overcome the problems in the prior art, the invention provides a novel method for quickly and automatically adjusting the exciting force in a modal test, the exciting force is automatically adjusted by using an optimization algorithm, the force adjusting process is simple, and the test time is greatly shortened.
The invention relates to a method for quickly and automatically adjusting excitation force, which comprises the following steps:
firstly, according to the characteristics of a tested structure or by using an impact force hammering method, determining the number of exciting forces and the position of an exciting point by enabling a vibration exciter to be positioned at a point where each interested order mode has enough displacement;
step two, carrying out multipoint sine frequency sweep or multipoint random excitation test on the tested structure to obtain a frequency response function matrix and modal frequency of the structure;
thirdly, a purity indication function is used as a fitness function in the algorithm, and a group of optimal excitation force amplitudes can be obtained after the particle swarm optimization iteration;
and step four, sending the optimal excitation force amplitude to the vibration exciter to enable the vibration exciter to quickly and automatically adjust the excitation force.
Wherein the differential equation of vibration of the structure under test is expressed as follows
Wherein M, C and K are a mass matrix, a damping matrix and a rigidity matrix of the tested structure respectively; { x (t) } and { f (t) } are the displacement vector of the structure under test and the external force vector acting on the structure, respectively.
Further, the second step specifically includes the following steps:
step 2.1, the system is excited by simple harmonic as
{f(t)}=Fejωt
In the formula: f is an excitation force amplitude array, and the steady-state response of the system can be expressed as
{x(t)}=X(ω)ejωt
In the formula: x (omega) is a steady-state response amplitude array;
step 2.2, substituting the steady state response into the vibration differential equation to obtain a system response:
X(ω)=(K-ω2M+jωC)-1·F=H(ω)·F
wherein the frequency response function is H (omega) ═ K-omega2M+jωC)-1。
Further, the purity indication function in the third step is as follows:
as a preferable selection, the particle swarm algorithm in step three has the characteristics of intuitive background, simplicity and easy implementation, and wide adaptability to different types of functions, and is often used for optimizing complex nonlinear functions, combinatorial optimization and the like.
The particle swarm optimization is used for the automatic excitation force optimization in the patent and specifically comprises the following steps: setting the number of exciting forces in a test as d, and setting the number of exciting force vectors of each group in an optimization algorithm as m; each group of exciting force vector is F ═ F1,F2,…,Fm) (ii) a Single excitation force vector is Fi=(fi,1,fi,2,…,fi,d) (ii) a Velocity vector corresponding to excitation force vector is Vi=(vi,1,vi,2,…,vi,d) (ii) a Each excitation force vector uniquely corresponds to a pure mode indication function value P ═ P1,P2,…,Pm) Wherein the optimal excitation force vector is Pi=(pi,1,pi,2,…,pi,d) (ii) a The global optimal excitation force vector is Pg=(pg,1,pg,2,…,pg,d) (ii) a Each generation of excitation force vector update and corresponding velocity update is as follows:
vi,j(t+1)=vi,j(t)+c1r1[Pi,j-fi,j(t)]+c2r2[Pg,j-fi,j(t)]
fi,j(t+1)=fi,j(t)+vi,j(t+1),j=1,…,d
fi,j(t) is the exciting force vector at the moment t; v. ofi,j(t) is the speed corresponding to the exciting force vector at the moment t; f. ofi,j(t +1) is an exciting force vector at the moment (t + 1); v. ofi,j(t +1) is the speed corresponding to the excitation force vector at the moment (t + 1);
c1and c2A learning factor that is a non-negative constant; r is1And r2Are mutually independent pseudo-random numbers; d is the number of exciting forces in the modal test; pi,jThe optimal excitation force vector at the moment t; pg,jAnd (4) setting the optimal excitation force vector for all the excitation force vectors until the moment t. And obtaining the optimal excitation force vector through the iteration of an automatic optimization algorithm.
Compared with the prior art, the method for quickly and automatically adjusting the excitation force has the following beneficial effects: the force adjusting process in the pure mode test is simpler, and the exciting force amplitude which enables the pure mode indicating function to be very high can be obtained in a short time.
Drawings
In order to more clearly illustrate the technical solution of the present invention, the drawings needed to be used in the present invention will be briefly described below, and it is apparent that the drawings in the following description are only some embodiments of the present invention, and it is obvious for those skilled in the art that other drawings can be obtained based on these drawings without inventive labor.
FIG. 1 is a schematic diagram of simulation structure meshing.
Detailed Description
In order to make the objects, technical solutions and advantages of the present invention more apparent, embodiments of the present invention will be described in detail with reference to the accompanying drawings.
The method is explained in detail by the simulation of the double U-shaped beam and the attached drawings as follows:
when the number of the exciting forces of the tested structure and the positions of the exciting points are determined, a group of exciting force amplitudes must exist in the exciting force amplitude combination in the range of the vibration exciter, so that the purity indicating function reaches the maximum value. The method provided by the invention is to use an optimization algorithm to perform optimization in the amplitude and phase combination of the exciting force so as to realize automatic adjustment of the exciting force.
The invention relates to a method for quickly and automatically adjusting excitation force, which comprises the following steps:
firstly, according to the characteristics of a tested structure or by using an impact force hammering method, determining the number of exciting forces and the position of an exciting point by enabling a vibration exciter to be positioned at a point where each interested order mode has enough displacement;
step two, carrying out multipoint sine frequency sweep or multipoint random excitation test on the tested structure to obtain a frequency response function matrix and modal frequency of the structure;
thirdly, a purity indication function is used as a fitness function in the algorithm, and a group of optimal excitation force amplitudes can be obtained after the particle swarm optimization iteration;
and step four, sending the optimal excitation force amplitude to the vibration exciter to enable the vibration exciter to quickly and automatically adjust the excitation force.
Wherein the differential equation of vibration of the structure under test is expressed as follows
Wherein M, C and K are a mass matrix, a damping matrix and a rigidity matrix of the tested structure respectively; { x (t) } and { f (t) } are the displacement vector of the structure under test and the external force vector acting on the structure, respectively.
Further, the second step specifically includes the following steps:
step 2.1, the system is excited by simple harmonic as
{f(t)}=Fejωt
In the formula: f is an excitation force amplitude array, and the steady-state response of the system can be expressed as
{x(t)}=X(ω)ejωt
In the formula: x (omega) is a steady-state response amplitude array;
step 2.2, substituting the steady state response into the vibration differential equation to obtain a system response:
X(ω)=(K-ω2M+jωC)-1·F=H(ω)·F
wherein the frequency response function is H (omega) ═ K-omega2M+jωC)-1。
Further, the purity indication function in the third step is as follows:
as a preferable selection, the particle swarm algorithm in step three has the characteristics of intuitive background, simplicity and easy implementation, and wide adaptability to different types of functions, and is often used for optimizing complex nonlinear functions, combinatorial optimization and the like.
The particle swarm optimization is used for the automatic excitation force optimization in the patent and specifically comprises the following steps: setting the number of exciting forces in a test as d, and setting the number of exciting force vectors of each group in an optimization algorithm as m; each group of exciting force vector is F ═ F1,F2,…,Fm) (ii) a Single excitation force vector is Fi=(fi,1,fi,2,…,fi,d) (ii) a Velocity vector corresponding to excitation force vector is Vi=(vi,1,vi,2,…,vi,d) (ii) a The excitation force vector only corresponds to a pure mode indication function value P ═1,P2,…,Pm) Wherein the optimal excitation force vector is Pi=(pi,1,pi,2,…,pi,d) (ii) a Global maximumThe vector of the excellent exciting force is Pg=(pg,1,pg,2,…,pg,d) (ii) a Each generation of excitation force vector update and corresponding velocity update is as follows:
vi,j(t+1)=vi,j(t)+c1r1[Pi,j-fi,j(t)]+c2r2[Pg,j-fi,j(t)]
fi,j(t+1)=fi,j(t)+vi,j(t+1),j=1,…,d
fi,j(t) is the exciting force vector at the moment t; v. ofi,j(t) is the speed corresponding to the exciting force vector at the moment t; f. ofi,j(t +1) is an exciting force vector at the moment (t + 1); v. ofi,j(t +1) is the speed corresponding to the excitation force vector at the moment (t + 1);
c1and c2A learning factor that is a non-negative constant; r is1And r2Are mutually independent pseudo-random numbers; d is the number of exciting forces in the modal test; pi,jThe optimal excitation force vector at the moment t; pg,jAnd (4) setting the optimal excitation force vector for all the excitation force vectors until the moment t. And obtaining the optimal excitation force vector through the iteration of an automatic optimization algorithm.
Example 1
As shown in fig. 1, the measured structure provided in this example is a double U-shaped beam, which is made of steel, and has a density of 7.8, an elastic modulus of 210Gpa, and a poisson's ratio of 0.31. Its total length is 2500, height is 100, width is 100, and thickness is 10. And (3) deriving an overall mass matrix and an overall stiffness matrix of the model by using workbench, and introducing matrix data of the model and a modal damping matrix with a modal damping ratio of 0.05 into an algorithm for calculation. Since the first six-order mode of the simulation structure is a rigid body mode, we take the first five-order bending mode to verify the method.
Since we only need to obtain the ratio of the excitation forces of the excitation points, one of the excitation force amplitudes is set to 1 or-1. The results obtained with the method of this patent are as follows:
order of mode | 1 | 2 | 3 | 4 | 5 |
Mode frequency Hz | 7.54 | 7.68 | 10.87 | 22.22 | 33.86 |
Indicating function T | 0.94 | 0.99 | 0.99 | 0.99 | 0.99 |
Excitation force 1 | 1 | 1 | -1 | 1 | -1 |
|
1.02 | -0.99 | -0.99 | -0.94 | 0.98 |
Excitation point 1 | 7 | 7 | 4 | 15 | 7 |
|
19 | 19 | 40 | 21 | 19 |
It can be seen that the method of the present patent can be used to quickly and automatically adjust the exciting force. And through inspection, the normalized excitation array is basically consistent with the theoretical array.
The above description is only exemplary of the present application and should not be taken as limiting the present application, as any modification, equivalent replacement, or improvement made within the spirit and principle of the present application should be included in the protection scope of the present application.
Claims (6)
1. A method for fast and automatic tuning of excitation forces, characterized in that the method comprises the steps of:
firstly, according to the characteristics of a tested structure or by using an impact force hammering method, determining the number of exciting forces and the position of an exciting point by enabling a vibration exciter to be positioned at a point where each interested order mode has enough displacement;
step two, carrying out multipoint sine frequency sweep or multipoint random excitation test on the tested structure to obtain a frequency response function matrix and modal frequency of the structure;
thirdly, a purity indication function is used as a fitness function in the algorithm, and a group of optimal excitation force amplitudes can be obtained after the particle swarm optimization iteration;
and step four, sending the optimal excitation force amplitude to the vibration exciter to enable the vibration exciter to quickly and automatically adjust the excitation force.
2. The method of claim 1, wherein the differential equation of the vibration of the structure under test is expressed as follows
Wherein M, C and K are a mass matrix, a damping matrix and a rigidity matrix of the tested structure respectively; { x (t) } and { f (t) } are the displacement vector of the structure under test and the external force vector acting on the structure, respectively.
3. The method for rapidly and automatically adjusting excitation force according to claim 2, wherein the second step specifically comprises the following steps:
step 2.1, the system is excited by simple harmonic as
{f(t)}=Fejωt
In the formula: f is an excitation force amplitude array, and the steady-state response of the system can be expressed as
{x(t)}=X(ω)ejωt
In the formula: x (omega) is a steady-state response amplitude array;
step 2.2, substituting the steady state response into the vibration differential equation to obtain a system response:
X(ω)=(K-ω2M+jωC)-1·F=H(ω)·F
wherein the frequency response function is H (omega) ═ K-omega2M+jωC)-1。
5. the method according to claim 4, wherein the particle swarm algorithm in the third step is to generate a group of particles at random, that is, an excitation force vector and a velocity vector corresponding to the particle, and after a system frequency response function and a purity indication function are programmed into the particle swarm algorithm, each particle uniquely corresponds to a system response value and a purity indication function value, and an optimal particle is determined by the purity indication function value; the next generation of particles can be generated after the position and speed of the particles are updated in the algorithm.
6. The method for rapidly and automatically adjusting excitation force according to claim 5, wherein the excitation force vector updating in the combination of the particle swarm optimization and the excitation force automatic optimization is specifically as follows:
vi,j(t+1)=vi,j(t)+c1r1[Pi,j-fi,j(t)]+c2r2[Pg,j-fi,j(t)]
fi,j(t+1)=fi,j(t)+vi,j(t+1),j=1,…,d
fi,j(t) is the exciting force vector at the moment t; v. ofi,j(t) is the speed corresponding to the exciting force vector at the moment t; f. ofi,j(t +1) is an exciting force vector at the moment (t + 1); v. ofi,j(t +1) is the speed corresponding to the excitation force vector at the moment (t + 1);
c1and c2A learning factor that is a non-negative constant; r is1And r2Are mutually independent pseudo-random numbers; d is the number of exciting forces in the modal test; pi,jThe optimal excitation force vector at the moment t; pg,jAnd (4) setting the optimal excitation force vector for all the excitation force vectors until the moment t. And obtaining the optimal excitation force vector through the iteration of an automatic optimization algorithm.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011059206.3A CN112149308B (en) | 2020-09-30 | 2020-09-30 | Method for quickly and automatically adjusting excitation force |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011059206.3A CN112149308B (en) | 2020-09-30 | 2020-09-30 | Method for quickly and automatically adjusting excitation force |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112149308A true CN112149308A (en) | 2020-12-29 |
CN112149308B CN112149308B (en) | 2021-10-01 |
Family
ID=73894443
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202011059206.3A Active CN112149308B (en) | 2020-09-30 | 2020-09-30 | Method for quickly and automatically adjusting excitation force |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112149308B (en) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113032910A (en) * | 2021-04-13 | 2021-06-25 | 哈尔滨工业大学 | Layout optimization method for actuator of flexible spacecraft |
CN118468924A (en) * | 2024-06-27 | 2024-08-09 | 南京航空航天大学 | Mixed intelligent force adjusting method in structural pure mode test |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20050155431A1 (en) * | 2004-01-19 | 2005-07-21 | Mayumi Fukuyama | Vibration test system and method for structures |
CN101561342A (en) * | 2009-05-15 | 2009-10-21 | 北京工业大学 | System and method for measuring time-sharing quick steady-state sine sweep excitation frequency response function |
CN101598632A (en) * | 2009-07-07 | 2009-12-09 | 华中科技大学 | The measuring method of strong point axial rigidity and damping in the screw rod drive feeding system |
CN106570264A (en) * | 2016-10-28 | 2017-04-19 | 南京航空航天大学 | Method for quickly adjusting excitation force in pure modal test |
-
2020
- 2020-09-30 CN CN202011059206.3A patent/CN112149308B/en active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20050155431A1 (en) * | 2004-01-19 | 2005-07-21 | Mayumi Fukuyama | Vibration test system and method for structures |
CN101561342A (en) * | 2009-05-15 | 2009-10-21 | 北京工业大学 | System and method for measuring time-sharing quick steady-state sine sweep excitation frequency response function |
CN101598632A (en) * | 2009-07-07 | 2009-12-09 | 华中科技大学 | The measuring method of strong point axial rigidity and damping in the screw rod drive feeding system |
CN106570264A (en) * | 2016-10-28 | 2017-04-19 | 南京航空航天大学 | Method for quickly adjusting excitation force in pure modal test |
Non-Patent Citations (2)
Title |
---|
孟霞等: "复合材料板架结构中纯模态试验技术的应用", 《噪声与振动控制》 * |
郭家骅: "多点相位协调及纯模态试验研究", 《中国优秀硕士学位论文全文数据库 工程科技Ⅱ辑》 * |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113032910A (en) * | 2021-04-13 | 2021-06-25 | 哈尔滨工业大学 | Layout optimization method for actuator of flexible spacecraft |
CN113032910B (en) * | 2021-04-13 | 2022-11-11 | 哈尔滨工业大学 | Layout optimization method for flexible spacecraft actuator |
CN118468924A (en) * | 2024-06-27 | 2024-08-09 | 南京航空航天大学 | Mixed intelligent force adjusting method in structural pure mode test |
Also Published As
Publication number | Publication date |
---|---|
CN112149308B (en) | 2021-10-01 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN112149308B (en) | Method for quickly and automatically adjusting excitation force | |
Yam et al. | Theoretical and experimental study of modal strain analysis | |
Yamaguchi et al. | Chaotic vibrations of a cylindrical shell-panel with an in-plane elastic-support at boundary | |
CN110188512B (en) | Method for modeling single-duplex clamp pipe system | |
Bertolini et al. | Seismic noise filters, vertical resonance frequency reduction with geometric anti-springs: a feasibility study | |
CN112182697B (en) | High-precision dynamic measurement method for system tension of damping suspender | |
Mohd Zin et al. | Modal and FRF based updating methods for the investigation of the dynamic behaviour of a plate | |
CN110096779B (en) | Servo mechanism dynamic characteristic analysis method | |
CN108560619B (en) | The analysis method and system of dynamic test for pile foundation bearing capacity | |
CN102680214B (en) | Vibration isolation performance test device and method for corrugated pipe | |
CN102818686B (en) | Grid-control TWT metal grid mesh Modal Experimental Method | |
CN106570264A (en) | Method for quickly adjusting excitation force in pure modal test | |
CN116187053A (en) | Method and device for measuring uncertainty of arm support stress based on Monte Carlo method | |
Sujatha | Basics of experimental modal analysis | |
Froli et al. | Experimental static and dynamic tests on a large-scale free-form Voronoi grid shell mock-up in comparison with finite-element method results | |
Rezaiee-Pajand et al. | Vibration analysis of plane frames by customized stiffness and diagonal mass matrices | |
CN112231916A (en) | Modal quality measuring method | |
Zhang | Effect of damping coefficients on structural damage identification | |
Chen et al. | Experiments on fluidelastic vibrations of tube arrays | |
Miao et al. | A load identification application technology based on regularization method and finite element modified model | |
CN111783316B (en) | Modified Biot composite material frequency-dependent damping behavior prediction method | |
Miao et al. | Research Article A Load Identification Application Technology Based on Regularization Method and Finite Element Modified Model | |
CN108304653A (en) | Crystal cultivates process environments vibration effect and quantifies rapid evaluation technology | |
Hong et al. | A Novel Structural Modification Method for Vibration Reduction: Stiffness Sensitivity Analysis with Principal Strain Application | |
Yin et al. | Optimal Selection of Exciting and Measuring Points for Modal Test |
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 |