Method for constructing virtual vibration system and related equipment

Abstract

The invention provides a method for constructing a virtual vibration system and related equipment. The method comprises the following steps: constructing a first model; the first model is constructed by using a finite element and statistical energy analysis hybrid method; the first model comprises at least one measuring point; based on the obtained response of each measuring point, correcting at least two internal loss factors in the first model to obtain a second model; based on the second model, determining the relationship between the coil force and the vibration control spectrum in each frequency band range so as to determine a target coil force; and based on the target coil force and the second model, determining a third model to construct a virtual vibration system. According to the embodiment of the invention, the internal loss factor is corrected based on the measured data, a mathematical method without cross terms in Taylor expansion is ingeniously applied to the correction of the internal loss factor, so that an accurate model is established, and the excitation condition can be reversely deduced according to the function relationship between the fitting force and the response to construct the virtual vibration system.

Description

Technical Field

[0001] This application relates to the technical field of constructing virtual vibration systems, and particularly to a method for constructing a virtual vibration system and related devices. Background Art

[0002] The hybrid method of finite element and statistical energy analysis (Finite Element-Statistical Energy Analysis, FE-SEA) can make up for the deficiencies of the finite element analysis method (Finite Element, FEA) and the statistical energy analysis method (Statistical Energy Analysis, SEA). However, the parameters of the statistical energy subsystem (SEA) in the hybrid method of finite element and statistical energy analysis (FE-SEA) (such as coupling loss factor, self-loss factor, etc.) are not easy to determine and still have a great impact on the prediction results. It is also difficult to accurately apply vibration excitation at the control points for complex hybrid models of finite element and statistical energy analysis.

[0003] For the vibration system of a shaker-fixture-product composed of a shaker, a fixture, and a product, an FE-SEA model is established. A part of the shaker, fixture, and product is modeled by finite element and divided into FE subsystems. The FE model needs to be calibrated through modal tests. Another part of the product is modeled by the statistical energy method and divided into SEA subsystems. The SEA subsystems need to determine the coupling loss factor and self-loss factor to finally determine the vibration response.

[0004] To determine the coupling loss factor and self-loss factor of the SEA subsystem, in addition to the theoretical calculation method, there is also an experimental method. The theoretical calculation method requires the product object to be regular, and the calculation results for irregular objects are also inaccurate. The experimental methods include the power input method, the energy ratio method, etc. Among them, the power input method is a commonly used experimental method that can obtain both the internal loss factor and the coupling loss factor at the same time. Generally, the experimental method requires the product divided into SEA subsystems to be suspended separately by rubber ropes, sensors to be arranged on the product, and the product to be excited by an exciter equipped with an impedance head, and the input parameters of the relevant exciter and the output parameters of the sensors to be measured, and then the coupling loss factor and self-loss factor of the product are calculated with certain theoretical support. However, the internal loss factor and the coupling loss factor cannot be directly obtained by using the vibration test data of the shaker-fixture-product vibration system. Summary of the Invention

[0005] In view of this, the purpose of this application is to propose a method for constructing a virtual vibration system and related devices.

[0006] Based on the above purpose, this application provides a method for constructing a virtual vibration system, including:

[0007] Build a first model; the first model is built using a hybrid method of finite element and statistical energy analysis; the first model includes at least one measurement point;

[0008] Based on the responses of each of the measurement points obtained, perform correction processing on at least two internal loss factors in the first model to obtain a second model;

[0009] Based on the second model, determine the relationship between the coil force magnitude and the vibration control spectrum within each frequency band range to determine the target coil force;

[0010] Based on the target coil force and the second model, determine a third model to construct a virtual vibration system.

[0011] In a possible implementation manner, the performing correction processing on at least two internal loss factors in the first model based on the responses of each of the measurement points obtained to obtain a second model includes:

[0012] Based on the responses of each of the measurement points obtained and the internal loss factors, build a first correction model;

[0013] Perform Taylor expansion on the first correction model to obtain a second correction model;

[0014] Based on the second correction model, perform correction processing on the internal loss factors to obtain the second model.

[0015] In a possible implementation manner, the first model includes at least one subsystem; each of the internal loss factors corresponds one-to-one to each of the subsystems;

[0016] The performing correction processing on the internal loss factors based on the second correction model to obtain the second model includes:

[0017] Determine the center frequency and the first parameter corresponding to the center frequency;

[0018] Change the internal loss factor corresponding to one of the subsystems, keep the internal loss factors corresponding to other subsystems as the first parameter, record the responses corresponding to the measurement points, and fit the first relationship between the responses corresponding to the measurement points and the changed internal loss factor;

[0019] Change the internal loss factor corresponding to another different subsystem, keep the internal loss factors corresponding to other subsystems as the first parameter, record the responses corresponding to the measurement points, and fit the second relationship between the responses corresponding to the measurement points and the internal loss factor corresponding to the other different subsystem;

[0020] Based on the first relationship and the second relationship, perform correction processing on the internal loss factors to obtain the second model.

[0021] In a possible implementation, the method of correcting the internal loss factor based on the first relationship and the second relationship to obtain the second model includes:

[0022] Derive the first relationship to obtain a first result;

[0023] Derive the second relationship to obtain a second result;

[0024] Determine a correction value of the internal loss factor based on the first result, the second result, and the second correction model;

[0025] Determine the second model based on the correction value.

[0026] In a possible implementation, the method of determining the relationship between the coil force magnitude and the vibration control spectrum in each frequency band range based on the second model to determine the target coil force includes:

[0027] Based on the second model, gradually change the magnitude of the input coil force and record the vibration control spectrum corresponding to each coil force;

[0028] Perform fitting processing on the data points formed by the coil force and the vibration control spectrum to obtain a fitting result;

[0029] Calculate the target coil force based on the obtained target vibration control spectrum and the fitting result.

[0030] In a possible implementation, the method of calculating the target coil force based on the obtained target vibration control spectrum and the fitting result includes:

[0031] Calculate the value of the sum of squares of errors and the square of the correlation coefficient between the target vibration control spectrum and the fitting result to determine the target coil force.

[0032] Based on the same inventive concept, an embodiment of the present application further provides a device for constructing a virtual vibration system, including:

[0033] A construction module configured to construct a first model; the first model is constructed by using a hybrid method of finite element and statistical energy analysis; at least one measurement point is included in the first model;

[0034] A correction module configured to correct at least two internal loss factors in the first model based on the response obtained for each measurement point to obtain a second model;

[0035] A determination module configured to determine the relationship between the coil force magnitude and the vibration control spectrum in each frequency band range based on the second model to determine the target coil force;

[0036] A system module, configured to determine a third model based on the target coil force and the second model to construct a virtual vibration system.

[0037] Based on the same inventive concept, an embodiment of the present application further provides an electronic device, including a memory, a processor, and a computer program stored on the memory and executable on the processor. When the processor executes the program, the method for constructing a virtual vibration system as described in any one of the above is implemented.

[0038] Based on the same inventive concept, an embodiment of the present application further provides a non-transitory computer-readable storage medium. The non-transitory computer-readable storage medium stores computer instructions, and the computer instructions are used to cause the computer to execute the method for constructing a virtual vibration system as described in any one of the above.

[0039] Based on the same inventive concept, an embodiment of the present application further provides a computer program product, which includes computer program instructions, and the computer instructions are used to cause the computer program product to execute the method for constructing a virtual vibration system as described in any one of the above.

[0040] As can be seen from the above, the method and related devices for constructing a virtual vibration system provided by the present application construct a first model; the first model is constructed by using a hybrid method of finite element and statistical energy analysis; at least one measurement point is included in the first model; based on the response of each of the obtained measurement points, at least two internal loss factors in the first model are corrected to obtain a second model; based on the second model, the relationship between the coil force magnitude and the vibration control spectrum in each frequency band range is determined to determine the target coil force; based on the target coil force and the second model, a third model is determined to construct a virtual vibration system. The embodiment of the present application makes up for the deficiencies of a pure finite element model or a pure statistical energy model by establishing an FE-SEA vibration system composed of a vibration table - fixture - product. After establishing the FE-SEA model of the vibration table - fixture - product, the internal loss factor is corrected based on the measured data, so as to establish an accurate model. Specifically, the internal loss factor is determined by combining the actual response, and the consistency between the simulation result and the physical result is good, and the power spectral density error is very small. In order to facilitate the comparison between the simulation result and the physical test result, it is necessary to make the control spectrum of the control point after applying excitation on the FE-SEA model consistent with the control of the physical test. That is, if the physical test uses multi-point average control, then the average value of multiple control points on the model after applying excitation on the model is consistent with the control spectrum of the physical test, which is convenient for predicting the response under any vibration level. Description of the Drawings

[0041] To more clearly illustrate the technical solutions in the present application or related technologies, the following will briefly introduce the accompanying drawings required for use in the embodiments or related technology descriptions. Obviously, the accompanying drawings in the following descriptions are only embodiments of the present application. For those of ordinary skill in the art, without creative efforts, other accompanying drawings can also be obtained based on these drawings.

[0042] Figure 1 Schematic flowchart of the method for constructing a virtual vibration system according to an embodiment of the present application;

[0043] Figure 2 Schematic structural diagram of the FE-SEA model according to an embodiment of the present application;

[0044] Figure 3 Schematic structural diagram of the device for constructing a virtual vibration system according to an embodiment of the present application;

[0045] Figure 4 Schematic structural diagram of an electronic device according to an embodiment of the present application. Detailed implementation manners

[0046] To make the objectives, technical solutions, and advantages of the present application clearer and more understandable, the following further elaborates on the present application in detail with reference to specific embodiments and the accompanying drawings.

[0047] It should be noted that unless otherwise defined, the technical terms or scientific terms used in the embodiments of the present application should have the ordinary meaning understood by those of ordinary skill in the art to which the present application belongs. The "first", "second", and similar terms used in the embodiments of the present application do not indicate any order, quantity, or importance, but are only used to distinguish different components. The terms such as "including" or "comprising" mean that the elements or objects appearing before this word cover the elements or objects listed after this word and their equivalents, without excluding other elements or objects. The terms such as "connected" or "coupled" are not limited to physical or mechanical connections, but may include electrical connections, whether direct or indirect. The terms such as "upper", "lower", "left", "right", etc. are only used to represent relative positional relationships, and when the absolute position of the object being described changes, the relative positional relationship may also change accordingly.

[0048] It can be understood that before using the technical solutions of the various embodiments of the present disclosure, the types, usage scopes, usage scenarios, etc. of the personal information involved will be informed to the user in an appropriate manner and the user's authorization will be obtained.

[0049] For example, when a user's active request is received, a prompt message is sent to the user to clearly prompt the user that the operation requested to be executed will require obtaining and using the user's personal information. Thus, the user can autonomously choose whether to provide personal information to software or hardware such as an electronic device, an application program, a server, or a storage medium that executes the operation of the technical solution of the present disclosure according to the prompt message.

[0050] As an optional but non-limiting implementation manner, the manner of sending a prompt message to the user in response to receiving the user's active request may be, for example, in the form of a pop-up window, and the prompt message may be presented in text in the pop-up window. In addition, the pop-up window may also carry selection controls for the user to select "agree" or "disagree" to provide personal information to the electronic device.

[0051] It can be understood that the above process of notifying and obtaining user authorization is only illustrative and does not limit the implementation manner of the present disclosure. Other manners that comply with relevant laws and regulations can also be applied to the implementation manner of the present disclosure.

[0052] As described in the background art section, the hybrid method of finite element and statistical energy analysis (Finite Element-Statistical Energy Analysis, FE-SEA) can make up for the deficiencies of the finite element analysis method (Finite Element, FEA) and the statistical energy analysis method (Statistical Energy Analysis, SEA). However, the parameters (coupling loss factor, self-loss factor, etc.) of the statistical energy subsystem (SEA) in the hybrid method of finite element and statistical energy analysis (FE-SEA) are not easy to determine and still have a great impact on the prediction result. It is also difficult to accurately apply a vibration excitation at the control point for a complex hybrid model of finite element and statistical energy analysis.

[0053] An FE-SEA model is established for the vibration table-fixture-product vibration system composed of a vibration table, a fixture, and a product. A part of the vibration table, the fixture, and the product is modeled by finite element and divided into FE subsystems. The FE model needs to be calibrated through a modal test. Another part of the product is modeled by the statistical energy method and divided into SEA subsystems. The SEA subsystems need to determine the coupling loss factor and the self-loss factor to finally determine the vibration response.

[0054] In addition to theoretical calculation methods, experimental methods can also be used to determine the coupling loss factor and self-loss factor of the SEA subsystem. Theoretical calculation methods require regular product objects, and the calculation results for irregular objects are inaccurate. Experimental methods include the power input method, the energy ratio method, etc. Among them, the power input method is a commonly used experimental method that can obtain both the internal loss factor and the coupling loss factor simultaneously. Generally, for products divided into SEA subsystems, they are individually suspended with rubber ropes, sensors are arranged on the products, excited by an exciter equipped with an impedance head, and the input parameters of the relevant exciter and the output parameters of the sensors are measured. Based on certain theoretical support, the coupling loss factor and self-loss factor of the products are calculated. However, the internal loss factor and the coupling loss factor cannot be directly obtained using the vibration table-fixture-product vibration system combined with existing vibration test data.

[0055] Considering the above, the embodiments of this application propose a method for constructing a virtual vibration system by constructing a first model. The first model is constructed using the finite element and statistical energy analysis hybrid method. The first model includes at least one measurement point. Based on the response obtained for each measurement point, correction processing is performed on at least two internal loss factors in the first model to obtain a second model. Based on the second model, the relationship between the coil force magnitude and the vibration control spectrum within each frequency band range is determined to determine the target coil force. Based on the target coil force and the second model, a third model is determined to construct a virtual vibration system. The embodiments of this application make up for the deficiencies of a pure finite element model or a pure statistical energy model by establishing an FE-SEA vibration system composed of a vibration table-fixture-product. After establishing the FE-SEA model of the vibration table-fixture-product, the internal loss factor is corrected based on the measured data to establish an accurate model. Specifically, the internal loss factor is determined by combining the actual response, and the consistency between the simulation results and the physical results is good, and the power spectral density error is very small. To facilitate the comparison between the simulation results and the physical test results, the control spectrum of the control points after applying excitation on the FE-SEA model needs to be consistent with that of the physical test. That is, if the physical test uses multi-point average control, then the average value of multiple control points on the model after applying excitation on the model is consistent with the physical test control spectrum, which is convenient for predicting the response at any vibration level.

[0056] Hereinafter, the technical solutions of the embodiments of this application will be described in detail through specific embodiments.

[0057] Refer to Figure 1 , the method for constructing a virtual vibration system according to the embodiments of this application includes the following steps:

[0058] Step S101, construct a first model. The first model is constructed using the finite element and statistical energy analysis hybrid method. The first model includes at least one measurement point.

[0059] Step S102: Based on the responses of each of the measurement points obtained, perform correction processing on at least two internal loss factors in the first model to obtain a second model.

[0060] Step S103: Based on the second model, determine the relationship between the coil force magnitude and the vibration control spectrum within each frequency band range to determine the target coil force.

[0061] Step S104: Based on the target coil force and the second model, determine a third model to construct a virtual vibration system.

[0062] Regarding step S101, first, it is necessary to construct the first model using the hybrid method of finite element and statistical energy analysis.

[0063] Reference Figure 2 , which is a schematic diagram of the FE-SEA model structure of the embodiment of the present application.

[0064] The structure consists of three parts: a vibration table, a fixture, and a product. The product is a cylindrical shell-like structure with three layers of circular flanges distributed in the upper, middle, and lower parts. The vibration table, fixture, and product flanges are modeled using finite elements, that is, the FE subsystem, also called the deterministic subsystem. The cylindrical shell surface of the product is modeled using statistical energy and is divided into upper and lower two SEA subsystems by the middle flange, also called the stochastic subsystems.

[0065] Only by using finite element modeling can the dynamic responses of the corresponding measurement points be obtained. The statistical energy method can only obtain the responses of the entire SEA subsystem and cannot obtain the response of a specific point. In order to compare with the results of physical tests, finite element modeling is used for the upper, middle, and lower flanges of the product (sensors are arranged on the flanges in physical tests).

[0066] After constructing the basic virtual vibration system, for step S102, perform correction processing on the internal loss factors.

[0067] In some embodiments, the performing correction processing on at least two internal loss factors in the first model based on the responses of each of the measurement points obtained to obtain a second model includes: constructing a first correction model based on the responses of each of the measurement points obtained and the internal loss factors; performing Taylor expansion on the first correction model to obtain a second correction model; and performing correction processing on the internal loss factors based on the second correction model to obtain the second model.

[0068] In some embodiments, the first model includes at least one subsystem; each of the internal loss factors corresponds to each of the subsystems one by one; the process of correcting the internal loss factors based on the second correction model to obtain the second model includes: determining a center frequency and a first parameter corresponding to the center frequency; changing the internal loss factor corresponding to one of the subsystems, keeping the internal loss factors corresponding to the other subsystems as the first parameter, recording the response corresponding to the measurement point, and fitting a first relationship between the response corresponding to the measurement point and the changed internal loss factor; changing the internal loss factor corresponding to another different subsystem, keeping the internal loss factors corresponding to the other subsystems as the first parameter, recording the response corresponding to the measurement point, and fitting a second relationship between the response corresponding to the measurement point and the internal loss factor corresponding to the other different subsystem; correcting the internal loss factors based on the first relationship and the second relationship to obtain the second model.

[0069] In some embodiments, the process of correcting the internal loss factors based on the first relationship and the second relationship to obtain the second model includes: taking the derivative of the first relationship to obtain a first result; taking the derivative of the second relationship to obtain a second result; determining a correction value of the internal loss factor based on the first result, the second result, and the second correction model; determining the second model based on the correction value.

[0070] In some embodiments, based on the second model, determining the relationship between the coil force magnitude and the vibration control spectrum within each frequency band range to determine the target coil force includes: based on the second model, gradually changing the input coil force magnitude and recording the vibration control spectrum corresponding to each coil force; performing a fitting process on the data points formed by the coil force and the vibration control spectrum to obtain a fitting result; calculating the target coil force based on the obtained target vibration control spectrum and the fitting result.

[0071] Specifically, apply the Taylor expansion mathematical method to the correction of the internal loss factor. For two internal loss factors η1 and η2, it is speculated through theoretical analysis that there is no cross term after Taylor expansion. Fix one parameter, for example, take η1 = η 1,k , select a series of η2, and cleverly fit the relationship between the response and η2 according to the measured data. Taking the derivative of it can obtain That is A series of η1 and η2 can be obtained according to this idea. In addition to the parameter of the internal loss factor, there are also parameters such as the coupling loss factor and the modal density. The determination of these parameters can follow the existing methods. Details are as follows:

[0072] First, a model for correcting the internal loss factor based on measured data was proposed. According to the number of statistical energy analysis subsystems divided, a mathematical model (the first correction model) between the measured point response and the internal loss factor was established (where f represents the center frequency, k represents the measured point number, η i represents the internal loss factor of subsystem i, n represents the number of subsystems, Φ(η1, η2, …, η i , …, η n ) represents a function of η1, η2, …, η i , …, η n ), and then according to the power flow balance equation of the statistical energy analysis subsystem:

[0073]

[0074] where, represents a function of η i , and L represents a constant term.

[0075] Then

[0076]

[0077] Performing a Taylor expansion of Φ(η1, η2, …, η i , …, η n ) at (η 1,k , η 2,k , …, η i,k , …, η n,k ) gives (the second correction model):

[0078]

[0079] Based on the above formula, it can be seen that to obtain the internal loss factor correction model, it is necessary to obtain

[0080] Therefore, in the subsequent calculation process, only η i is changed, and the remaining η j (i ≠ j) is kept unchanged (the theoretical reference value given by commercial software can be taken). After correcting η i multiple times, the corresponding can be obtained, and then the least squares fitting is used to obtain It is considered that Through the above method,

[0081] According to the deviation range of ±3dB of the measured data of the vibration test, combined with the internal loss factor correction model Φ(η1, η2, …, η i , …, ηn ) η i (i = 1, 2, …, n) ranges, and the system model is accurately established.

[0082] If there are multiple measurement points (k = 1, 2, …, m), the internal loss factor η of the jth subsystem j can take the average value of multiple measurement points:

[0083]

[0084] In a feasible embodiment, the established hybrid model includes two statistical energy analysis subsystems. The internal loss factor of the statistical energy analysis subsystem affects the final calculation result and needs to be corrected. The relationship between the responses at the positions of each measurement point and the loss factor can be expressed as:

[0085]

[0086] where k represents the measurement point number, f represents the center frequency, and η1 and η2 respectively represent the internal loss factors of the two statistical energy analysis subsystems.

[0087] When the center frequency f is determined, the above formula can be expressed as:

[0088]

[0089] Perform Taylor expansion on the above formula at (η 1,k , η 2,k ), and it can be expressed as:

[0090]

[0091] According to the internal loss factor and the hybrid finite element - statistical energy analysis theory, it can be known that there is no f term (uv ≠ 0) in psd (η1, η2), that is Therefore, psd f (η1, η2) can be obtained through the following steps:

[0092] Determine the center frequency f and the corresponding (η 1,k , η 2,k ) at this center frequency;

[0093] Ensure η1 = η 1,k , select a series of η2, calculate the corresponding relationship between the response and η2 and fit psd f (η 1,k , η2), and take the derivative of it to obtain That is

[0094] Ensure η2 = η 2,k , select a series of η1, calculate the corresponding relationship between the response and η1 and fit psd f (η1, η 2,k ), take the derivative of it to get That is

[0095] Calculate psd f (η 1,k , η 2,k ), obtain psd according to the above Taylor expansion f (η1, η2);

[0096] Change the center frequency f, repeat the above steps, and obtain psd(f, η1, η2).

[0097] According to the above steps, use the vertical table target control spectrum 1 vibration test results to obtain the internal loss factor and apply it to other test conditions.

[0098] Taking the 2000 Hz at the 6th measurement point as an example, after establishing the model in VAone, VAone will give a theoretically calculated internal loss factor. Take it as (η 1,k , η 2,k ), where η 1,k = 0.0004509, η 2,k = 0.002215. The values of η1 and η2 are 0.001, 0.005, 0.008, and 0.01. Calculate according to the above steps to get:

[0099]

[0100] Taking the derivative to get:

[0101]

[0102] The psd can be calculated through the hybrid model f (η 1,k , η 2,k ) = 0.002319, according to the above Taylor expansion, we can get:

[0103] psd f (η1, η2)

[0104] = 0.002319 - 0.4822(η1 - 0.0004509)

[0105] - 0.4655(η2 - 0.002215) = 0.0036 - 0.4822η1 - 0.4655η2

[0106] According to the test results and considering the 3dB error, we can get psdf (η1, η2) ∈ [1.3636e - 5, 5.4285e - 5]. Since η1, η2 > 0, the values of η1 and η2 can be determined from the above range. Take η1 = 0.00341 and η2 = 0.004143.

[0107] By performing the above - mentioned analysis on the remaining center frequencies, the corrected results of the internal loss factor (as shown in Table 11) can be obtained and used in the calculations for other test conditions to verify the accuracy of the hybrid model and the corrected results.

[0108] Table 1 Internal loss factor

[0109]

[0110]

[0111] Furthermore, for step S103, determine the target coil force.

[0112] In some embodiments, calculating the target coil force based on the obtained target vibration control spectrum and the fitting result includes: calculating the value of the sum of squares of the errors and the square of the correlation coefficient between the target vibration control spectrum and the fitting result, and determining the target coil force.

[0113] In this embodiment, during the vibration test, closed - loop control is performed by monitoring the dynamic response of the tabletop through control points to make the tabletop reach the preset vibration conditions. While in the virtual test system, the error of the given vibration conditions is directly used for response prediction, and different vibration conditions need to be indirectly simulated by changing the magnitude of the coil force. Therefore, when using the hybrid FE - SEA model for response prediction, the functional relationship between the magnitude of the coil force and the vibration control spectrum within each frequency band needs to be clearly given:

[0114]

[0115] Then, judge the accuracy of the fitting result according to the sum of squares of the errors of the fitting result and the square of the correlation coefficient between the measured data and the inferred data. When the sum of squares of the errors is close to 0 and the correlation coefficient is closer to 1, it indicates a better fitting effect. According to the fitting result and combined with the actual input control spectrum, the magnitude of the input force for each frequency band can be finally obtained.

[0116] In a feasible embodiment, the table shows the vibration test conditions. When conducting tests on the vibration table, the calculation conditions need to be determined for response prediction.

[0117] For the convenience of comparison, the points in the finite element - statistical energy model corresponding to the control points in the actual test are used as the control points in the model. After obtaining the transfer functions of these points in the finite element - statistical energy model, the transfer functions of these points with respect to the entire control system are further obtained. According to the random vibration control algorithm, the control spectra of these points are obtained, and this control spectrum is used as the input of the calculation conditions.

[0118] Table 2 Vibration test conditions

[0119] Frequency Magnitude 20 to 80 Hz +3 dB / Oct 80 to 350 Hz 0.04 g2 / Hz 350 to 2000 Hz -3 dB / Oct

[0120] When performing the hybrid finite element - statistical energy analysis calculation, the 1 / 3 octave is used for analysis, and the calculation conditions also need to be converted into the 1 / 3Oct form. The conversion can be carried out according to the actual calculation process of η1 and η2 in the above formula:

[0121]

[0122] Δf can be calculated according to the frequency range of 1 / 3Oct within the range of 200 - 2000Hz in Table 3.

[0123] Table 3 Center frequencies and frequency ranges of 1 / 3 octave

[0124] Center frequency Frequency range Center frequency Frequency range 200 180～224 800 710～900 250 224～280 1000 900～1120 310 280～355 1250 1120～1400 400 355～560 1600 1400～1800 500 450～560 2000 1800～2240 630 560～710

[0125] During the calculation, a force is applied to the coil, and whether the calculation conditions are met is controlled through five control points on the fixture. Therefore, it is necessary to clarify the relationship between the force applied to the coil and the responses of the measurement points on the fixture. By changing the force magnitudes in each frequency band, the responses of the five measurement points are calculated through the hybrid model, and the least - squares fitting is performed to obtain the relationship function between the force and the response:

[0126]

[0127] Taking the vibration table test conditions as an example, by changing the force magnitudes on the coil, 10N, 100N, 1000N, 10000N, and 100000N are respectively selected as the output conditions, and the hybrid model is used for calculation to obtain the responses of the five control points. The cftool toolbox of MATLAB is used for fitting to obtain the fitting results. Table 4 shows the sum of squared errors of the fitting results and the square value of the correlation coefficient between the measured data and the inferred data. Among them, the sum of squared errors is close to 0, and the correlation coefficient is 1, indicating that the fitting effect is good. According to the fitting results and the 1 / 3Oct - form PSD spectrum of the actual input conditions, the input force magnitudes in each frequency band can be obtained, as shown in

[0128] the table.

[0129] Table 4 Sum of squared errors of the fitting results and square value of the correlation coefficient

[0130]

[0131]

[0132] Table 5 Spectrum of coil force under each input condition

[0133]

[0134] In the embodiment of the present application, a finite element statistical energy model of a vibration system composed of a vibration table, a fixture, and a product is established, accurate parameters are determined, and accurate excitation is applied to achieve a control effect consistent with the physical test. The simulation results of the response points are consistent with the physical test results. The internal loss factor can be determined by combining the actual response, and the consistency between the simulation results and the physical results is good, and the power spectral density is within ±3 dB. The excitation application part in the model is the same as the actual application part of the vibration table, and the simulation model is correspondingly applied to the moving coil of the finite element part. The achieved control effect is consistent with the actual situation. If the actual test is multi-point average value control, the functional relationship between the force (excitation) applied to the moving coil at different frequency bands and the average value of the responses of each control point (response) is fitted, and based on this functional relationship, the average value of each control point (response) applied to the FE-SEA model is adjusted by adjusting the moving coil force (excitation) to be consistent with the physical control spectrum. Similarly, the average value, maximum value, and minimum value control methods can all be realized. The internal loss factor and coupling loss factor can be directly obtained by using the previous vibration test data in combination with software, without the need to prepare impedance heads, exciters, etc. like the power input method and energy ratio method, and use rubber ropes to suspend the SEA subsystem, arrange sensors, etc. to complete another test, and obtain the corresponding data and process it under certain theoretical guidance to obtain the results, but directly use it. After obtaining a model with good consistency with the physical test under a control spectrum, the target control spectrum can be obtained by adjusting the excitation, and the prediction results under different vibration levels can be obtained. It is convenient to perform response prediction under any vibration level.

[0135] In addition, after establishing the FE-SEA model of the vibration table-fixture-product, the internal loss factor is corrected based on the measured data, and the Taylor expansion without cross terms is cleverly used, and the mathematical method is applied to the correction of the internal loss factor to establish an accurate model. After the Taylor expansion, the derivative of the Taylor expansion is deduced backwards, and finally the internal loss factor is determined. The present application can obtain parameters according to the actual vibration test data without the need to design another test separately to obtain parameters. Since the parameters are directly obtained from the physical test results, the consistency with the physical test results is better.

[0136] In order to facilitate the comparison between the simulation results and the physical test results, the control spectrum after applying excitation to the model needs to be consistent with that of the physical test. That is, if the multi-point average value control is adopted in the physical test, the average value of multiple control points on the model after applying excitation to the model is consistent with the control spectrum of the physical test. By changing the magnitude of the force on the coil, a series of forces such as 10N, 100N, 1000N, 10000N, 100000N, etc. are selected as the excitation conditions respectively, and the responses of a series of control points are obtained respectively. By fitting the functional relationship between the force and the response, for any desired control effect, the excitation conditions can be inversely deduced according to this fitted functional relationship. This application combines the least squares method with the physical process of the shaking table to obtain a vibration control effect consistent with that of the physical test, and corresponding creative work has been done.

[0137] After applying the least squares method to the fitting between the force of the moving coil of the shaking table and the response of the control point, a very good control effect consistent with the response of the control point in the physical test can be obtained, thus facilitating the comparison between the simulation results and the physical test results at the response point. On this basis, the excitation conditions can be adjusted according to this fitting relationship to obtain any target control spectrum, so as to predict the vibration response under any target control spectrum.

[0138] Finally, after constructing the FE-SEA model of the vibration system, the measured data is corrected for the internal loss factor, the vibration control excitation conditions are accurately determined, and the methods for determining other parameters commonly used (such as the coupling loss factor and modal density of the SEA subsystem, and verifying the parameters of the FE subsystem through modal tests, etc.) are used to construct the entire virtual vibration system for vibration prediction.

[0139] As can be seen from the above embodiments, the method for constructing a virtual vibration system according to the embodiments of the present application includes constructing a first model; the first model is constructed by using a hybrid method of finite element and statistical energy analysis; at least one measurement point is included in the first model; based on the response of each of the measurement points obtained, at least two internal loss factors in the first model are corrected to obtain a second model; based on the second model, the relationship between the coil force magnitude and the vibration control spectrum in each frequency band range is determined to determine the target coil force; based on the target coil force and the second model, a third model is determined to construct a virtual vibration system. The embodiments of the present application make up for the deficiencies of a pure finite element model or a pure statistical energy model by establishing an FE-SEA vibration system composed of a vibration table-fixture-product. After establishing the FE-SEA model of the vibration table-fixture-product, the internal loss factor is corrected based on the measured data, so as to establish an accurate model. Specifically, the internal loss factor is determined in combination with the actual response, and the consistency between the simulation result and the physical result is good, and the power spectral density error is very small. In order to facilitate the comparison between the simulation result and the physical test result, it is necessary to make the control spectrum of the control point after applying excitation on the FE-SEA model consistent with the control of the physical test. That is, if the physical test uses multi-point average control, the average value of multiple control points on the model after applying excitation on the model is consistent with the control spectrum of the physical test, which is convenient for predicting the response under any vibration level.

[0140] It should be noted that the method of the embodiments of the present application can be executed by a single device, such as a computer or a server. The method of this embodiment can also be applied to a distributed scenario and completed by multiple devices cooperating with each other. In this case of a distributed scenario, one of the multiple devices can only execute one or more steps of the method of the embodiments of the present application, and these multiple devices will interact with each other to complete the described method.

[0141] It should be noted that some embodiments of the present application have been described above. Other embodiments are within the scope of the appended claims. In some cases, the actions or steps recited in the claims can be executed in a different order from that in the above embodiments and still achieve the desired result. Additionally, the processes depicted in the drawings do not necessarily require the specific order or continuous order shown to achieve the desired result. In certain embodiments, multi-tasking and parallel processing are also possible or may be advantageous.

[0142] Based on the same inventive concept, corresponding to the method of any of the above embodiments, the present application also provides an apparatus for constructing a virtual vibration system.

[0143] Referring to Figure 3 , the apparatus for constructing a virtual vibration system includes:

[0144] The building module 31 is configured to build a first model; the first model is built by using a hybrid method of finite element and statistical energy analysis; at least one measuring point is included in the first model;

[0145] The correction module 32 is configured to perform correction processing on at least two internal loss factors in the first model based on the responses of each of the measuring points obtained, so as to obtain a second model;

[0146] The determination module 33 is configured to determine the relationship between the coil force magnitude and the vibration control spectrum within each frequency band range based on the second model, so as to determine the target coil force;

[0147] The system module 34 is configured to determine a third model based on the target coil force and the second model, so as to build a virtual vibration system.

[0148] For the convenience of description, when describing the above device, various modules are described separately according to their functions. Of course, when implementing the present application, the functions of each module can be implemented in one or more software and / or hardware.

[0149] The device in the above embodiment is used to implement the corresponding method for building a virtual vibration system in any of the foregoing embodiments, and has the beneficial effects of the corresponding method embodiments, which will not be elaborated here.

[0150] Based on the same inventive concept, corresponding to the method in any of the above embodiments, the present application further provides an electronic device, including a memory, a processor, and a computer program stored on the memory and executable on the processor. When the processor executes the program, it implements the method for building a virtual vibration system described in any of the above embodiments.

[0151] Figure 4 FIG. shows a more specific schematic diagram of the hardware structure of the electronic device provided in this embodiment. The device may include: a processor 1010, a memory 1020, an input / output interface 1030, a communication interface 1040, and a bus 1050. Among them, the processor 1010, the memory 1020, the input / output interface 1030, and the communication interface 1040 are communicatively connected to each other inside the device through the bus 1050.

[0152] The processor 1010 may be implemented in a general-purpose CPU (Central Processing Unit), a microprocessor, an application-specific integrated circuit (ASIC), or one or more integrated circuits, etc., and is used to execute relevant programs to implement the technical solutions provided in the embodiments of this specification.

[0153] The memory 1020 can be implemented in the form of ROM (Read Only Memory), RAM (Random Access Memory), static storage devices, dynamic storage devices, etc. The memory 1020 can store the operating system and other application programs. When implementing the technical solutions provided in the embodiments of this specification through software or firmware, the relevant program codes are stored in the memory 1020 and called and executed by the processor 1010.

[0154] The input / output interface 1030 is used to connect to the input / output module to achieve information input and output. The input / output module can be configured as a component in the device (not shown in the figure) or externally connected to the device to provide corresponding functions. Among them, the input devices can include keyboards, mice, touchscreens, microphones, various sensors, etc., and the output devices can include displays, speakers, vibrators, indicator lights, etc.

[0155] The communication interface 1040 is used to connect to the communication module (not shown in the figure) to achieve communication interaction between this device and other devices. Among them, the communication module can achieve communication through wired means (such as USB, network cable, etc.) or through wireless means (such as mobile network, WIFI, Bluetooth, etc.).

[0156] The bus 1050 includes a path for transmitting information between various components of the device (such as the processor 1010, the memory 1020, the input / output interface 1030, and the communication interface 1040).

[0157] It should be noted that although the above device only shows the processor 1010, the memory 1020, the input / output interface 1030, the communication interface 1040, and the bus 1050, in the specific implementation process, the device may also include other components necessary for normal operation. In addition, those skilled in the art can understand that the above device may also only include the components necessary for implementing the solutions of the embodiments of this specification, and does not necessarily include all the components shown in the figure.

[0158] The electronic device in the above embodiments is used to implement the corresponding method for constructing a virtual vibration system in any of the foregoing embodiments, and has the beneficial effects of the corresponding method embodiments, which will not be elaborated here.

[0159] Based on the same inventive concept, corresponding to the method in any of the above embodiments, the present application also provides a non-transitory computer-readable storage medium storing computer instructions for causing the computer to execute the method for constructing a virtual vibration system as described in any of the foregoing embodiments.

[0160] The computer-readable medium of this embodiment includes permanent and non-permanent, removable and non-removable media, and information storage can be implemented by any method or technology. The information can be computer-readable instructions, data structures, program modules, or other data. Examples of computer storage media include, but are not limited to, phase change memory (PRAM), static random access memory (SRAM), dynamic random access memory (DRAM), other types of random access memory (RAM), read-only memory (ROM), electrically erasable programmable read-only memory (EEPROM), flash memory, or other memory technologies, compact disc read-only memory (CD-ROM), digital versatile disc (DVD), or other optical storage, magnetic cassette tapes, magnetic tape magnetic disk storage, or other magnetic storage devices, or any other non-transmission medium that can be used to store information accessible by a computing device.

[0161] The computer instructions stored in the storage medium of the above embodiment are used to cause the computer to execute the method of constructing a virtual vibration system as described in any one of the above embodiments, and have the beneficial effects of the corresponding method embodiments, which will not be elaborated here.

[0162] Based on the same inventive concept, corresponding to the method of constructing a virtual vibration system described in any of the above embodiments, the present disclosure also provides a computer program product, which includes computer program instructions. In some embodiments, the computer program instructions can be executed by one or more processors of the computer to cause the computer and / or the processor to execute the method of constructing a virtual vibration system. Corresponding to the execution subjects corresponding to the steps in each embodiment of the method of constructing a virtual vibration system, the processor executing the corresponding steps can belong to the corresponding execution subject.

[0163] The computer program product of the above embodiment is used to cause the computer and / or the processor to execute the method of constructing a virtual vibration system as described in any one of the above embodiments, and has the beneficial effects of the corresponding method embodiments, which will not be elaborated here.

[0164] Those of ordinary skill in the art should understand that the discussion of any of the above embodiments is only exemplary and is not intended to imply that the scope of the present application (including the claims) is limited to these examples; under the concept of the present application, the technical features in the above embodiments or different embodiments can also be combined, the steps can be implemented in any order, and there are many other variations in different aspects of the embodiments of the present application as described above, which are not provided in detail for the sake of brevity.

[0165] In addition, for simplicity of explanation and discussion, and so as not to make the embodiments of the present application difficult to understand, well-known power / ground connections to integrated circuit (IC) chips and other components may or may not be shown in the provided drawings. Further, the devices may be shown in block diagram form in order to avoid making the embodiments of the present application difficult to understand, and this also takes into account the fact that details of the implementation of such block diagram devices are highly dependent on the platform on which the embodiments of the present application are to be implemented (i.e., these details should be fully within the understanding of those skilled in the art). In cases where specific details (such as circuits) are set forth to describe exemplary embodiments of the present application, it will be apparent to those skilled in the art that the embodiments of the present application may be practiced without these specific details or with variations of these specific details. Accordingly, these descriptions are to be regarded as illustrative rather than restrictive.

[0166] Although the present application has been described in connection with specific embodiments thereof, many alternatives, modifications, and variations of these embodiments will be apparent to those of ordinary skill in the art in light of the foregoing description. For example, other memory architectures (such as dynamic RAM (DRAM)) may be used with the embodiments discussed.

[0167] Embodiments of the present application are intended to cover all such alternatives, modifications, and variations that fall within the broad scope of the appended claims. Therefore, any omissions, modifications, equivalent substitutions, improvements, etc., made within the spirit and principle of the embodiments of the present application shall be included within the protection scope of the present application.

Claims

1. A method for constructing a virtual vibration system, characterized in that Including: Constructing a first model; The first model is constructed by using a hybrid method of finite element and statistical energy analysis; The first model includes at least one measurement point; Based on the responses of each of the measurement points obtained, correcting at least two internal loss factors in the first model to obtain a second model; Based on the second model, determining the relationship between the coil force magnitude and the vibration control spectrum within each frequency band range to determine the target coil force; Based on the target coil force and the second model, determining a third model to construct a virtual vibration system.

2. The method according to claim 1, wherein The step of, based on the responses of each of the measurement points obtained, correcting at least two internal loss factors in the first model to obtain a second model, includes: Based on the responses of each of the measurement points obtained and the internal loss factors, constructing a first correction model; Performing a Taylor expansion on the first correction model to obtain a second correction model; Based on the second correction model, correcting the internal loss factors to obtain the second model.

3. The method according to claim 2, wherein The first model includes at least one subsystem; each of the internal loss factors corresponds to each of the subsystems one by one; The step of, based on the second correction model, correcting the internal loss factors to obtain the second model, includes: Determining a center frequency and a first parameter corresponding to the center frequency; Changing the internal loss factor corresponding to one of the subsystems, keeping the internal loss factors corresponding to the other subsystems as the first parameter, recording the response corresponding to the measurement point, and fitting a first relationship between the response corresponding to the measurement point and the changed internal loss factor; Changing the internal loss factor corresponding to another different subsystem, keeping the internal loss factors corresponding to the other subsystems as the first parameter, recording the response corresponding to the measurement point, and fitting a second relationship between the response corresponding to the measurement point and the internal loss factor corresponding to the other different subsystem; Based on the first relationship and the second relationship, correcting the internal loss factors to obtain the second model.

4. The method according to claim 3, characterized in that, The step of, based on the first relationship and the second relationship, correcting the internal loss factors to obtain the second model, includes: Taking the derivative of the first relationship to obtain a first result; Taking the derivative of the second relationship to obtain a second result; Based on the first result, the second result, and the second correction model, determining a correction value of the internal loss factor; Based on the correction value, determining the second model.

5. The method according to claim 1, characterized in that The step of, based on the second model, determining the relationship between the coil force magnitude and the vibration control spectrum within each frequency band range to determine the target coil force, includes: Based on the second model, gradually changing the magnitude of the input coil force and recording the vibration control spectrum corresponding to each coil force; Performing a fitting process on the data points formed by the coil force and the vibration control spectrum to obtain a fitting result; Based on the obtained target vibration control spectrum and the fitting result, calculating to obtain the target coil force.

6. The method according to claim 5, characterized in that The step of, based on the obtained target vibration control spectrum and the fitting result, calculating to obtain the target coil force, includes: Calculate the value of the square of the correlation coefficient of the sum of squared errors between the target vibration control spectrum and the fitting result, and determine the target coil force.

7. An apparatus for constructing a virtual vibration system, characterized in that, It includes: A construction module configured to construct a first model; The first model is constructed by using the hybrid method of finite element and statistical energy analysis; The first model includes at least one measurement point; A correction module configured to perform correction processing on at least two internal loss factors in the first model based on the responses of each of the measurement points obtained, to obtain a second model; A determination module configured to determine the relationship between the coil force magnitude and the vibration control spectrum within each frequency band range based on the second model, so as to determine the target coil force; A system module configured to determine a third model based on the target coil force and the second model, so as to construct a virtual vibration system.

8. An electronic device, comprising a memory, a processor, and a computer program stored on the memory and executable on the processor, characterized in that, When the processor executes the program, it implements the method according to any one of claims 1 to 6.

9. A non-transitory computer-readable storage medium storing computer instructions, characterized in that, The computer instructions are used to cause the computer to execute the method according to any one of claims 1 to 6.

10. A computer program product, including computer program instructions, when the computer program instructions run on a computer, cause the computer to execute the method according to any one of claims 1 to 6.

Images