Multi-point excitation dynamic load identification method, device and equipment and storage medium

Abstract

The invention discloses a multi-point excitation dynamic load identification method, device and equipment and a computer readable storage medium. According to the method, in a multi-degree-of-freedom motion system model, a preset time domain is divided into a plurality of time periods, every m continuous time periods serve as a dynamic load recognition unit, N degrees of freedom exist in the multi-degree-of-freedom motion system model, and in the preset time domain, a dynamic load is applied to each of the m degrees of freedom, the m dynamic load force values are constant in one identification unit, and m is smaller than or equal to N; and aiming at each identification unit, solving the equation set based on the dynamic load force values, and calculating to obtain each dynamic load force value of the identification unit, so that the technical problem that each independent load is difficult to accurately separate and identify when facing a multi-source excitation problem is solved, and the technical requirement of a modern complex engineering structure on multi-load accurate identification is met.

Description

Technical Field

[0001] This application relates to the field of multi-point excitation dynamic load identification, specifically to a method, apparatus, device, and computer-readable storage medium for multi-point excitation dynamic load identification. Background Technology

[0002] With the increasing complexity of modern engineering structures, such as aerospace vehicles, high-speed trains, and large mechanical equipment, these systems often experience dynamic loads from multiple locations simultaneously during operation. Accurate identification of these dynamic loads is crucial for structural health monitoring, fatigue life assessment, vibration control, and optimization of dynamic characteristics.

[0003] In related technologies, traditional dynamic load identification methods are mostly designed for single excitation sources or simplified to equivalent single-point excitation, and load identification is achieved by constructing system transfer functions or state-space equations.

[0004] However, complex structures in actual engineering are often subjected to multiple dynamic loads with different locations and characteristics simultaneously. These loads may have coupling effects, and the characteristics of each load (such as amplitude, frequency, and phase) vary greatly. Existing technologies can only identify single-point excited dynamic loads and cannot effectively handle multi-point excited situations. This makes it difficult to accurately separate and identify individual loads when facing multi-source excitation problems, and fails to meet the technical requirements of accurate identification of multiple loads in modern complex engineering structures. Summary of the Invention

[0005] This application provides a method, apparatus, device, and computer-readable storage medium for identifying multi-point excitation dynamic loads, which can solve the technical problem in the prior art that multiple dynamic loads cannot be identified.

[0006] In a first aspect, embodiments of this application provide a method for identifying multi-point excitation dynamic loads, the method comprising: In the multi-degree-of-freedom motion system model, the preset time domain is divided into several time periods, and each m consecutive time periods is a dynamic load identification unit. In the multi-degree-of-freedom motion system model, there are N degrees of freedom. In the preset time domain, a dynamic load is applied to each of the m degrees of freedom. The force values ​​of the m dynamic loads are constant within an identification unit, and m is less than or equal to N. For each identification unit, the dynamic load force value of each identification unit is calculated by solving the system of equations based on the dynamic load force value.

[0007] In conjunction with the first aspect, in one implementation, the step of calculating the dynamic load force value for each identification unit by solving a system of equations based on the dynamic load force value includes: In a multi-degree-of-freedom motion system model, the displacement and velocity of each degree of freedom corresponding to each time period are collected. In a recognition unit, there are m+1 time points.

[0008] In conjunction with the first aspect, in one embodiment, the equations for solving the dynamic load force value are as follows:

[0009] in, For the first Within the time period, the first i The work done by a dynamic load. For the first The increase in kinetic energy over a given period of time For the first The increase in potential energy over a given time period No. Damping energy dissipation over a time period, of which k The range of values ​​for is [1, m]. i The value range is [1, m].

[0010] In conjunction with the first aspect, in one embodiment, the kinetic energy increment, the potential energy increment, and the damping dissipation energy during the time period include: The increase in kinetic energy within the time period is the difference in system kinetic energy at the beginning and end of the time period; The potential energy increment within the time period is the difference in system potential energy at the beginning and end of the time period. The damping dissipation energy within the time period refers to the energy consumed by damping within a time period.

[0011] In conjunction with the first aspect, in one embodiment, after calculating each dynamic load force value of the identification unit, the process includes: Based on each dynamic load force value of each identification unit, a time-domain identification curve for each dynamic load within a preset time domain is obtained by fitting.

[0012] Secondly, embodiments of this application provide a multi-point excitation dynamic load identification device, the multi-point excitation dynamic load identification device comprising: Force application module: In the multi-degree-of-freedom motion system model, the preset time domain is divided into several time periods, and each m consecutive time period is a dynamic load identification unit. In the multi-degree-of-freedom motion system model, there are N degrees of freedom. Within the preset time domain, a dynamic load is applied to each of the m degrees of freedom. The force values ​​of the m dynamic loads remain constant within an identification unit, and m is less than or equal to N.

[0013] Calculation module: For each identification unit, the system of equations is solved based on the dynamic load force value to calculate the dynamic load force value of each identification unit.

[0014] In conjunction with the second aspect, in one implementation, the computing module includes:

[0015] in, For the first Within the time period, the first i The work done by a dynamic load. For the first The increase in kinetic energy over a given period of time For the first The increase in potential energy over a given time period No. Damping energy dissipation over a time period, of which k The range of values ​​for is [1, m]. i The value range is [1, m].

[0016] In conjunction with the second aspect, in one embodiment, the kinetic energy increment, the potential energy increment, and the damping dissipation energy during the time period include: The increase in kinetic energy within the time period is the difference in system kinetic energy at the beginning and end of the time period; The potential energy increment within the time period is the difference in system potential energy at the beginning and end of the time period. The damping dissipation energy within the time period refers to the energy consumed by damping within a time period.

[0017] Thirdly, embodiments of this application provide a multi-point excitation dynamic load identification device, the multi-point excitation dynamic load identification device including a processor, a memory, and a multi-point excitation dynamic load identification program stored in the memory and executable by the processor, wherein when the multi-point excitation dynamic load identification program is executed by the processor, it implements the steps of the multi-point excitation dynamic load identification method as described in the first aspect.

[0018] Fourthly, embodiments of this application provide a computer-readable storage medium storing a multi-point excitation dynamic load identification program, wherein when the multi-point excitation dynamic load identification program is executed by a processor, it implements the steps of the multi-point excitation dynamic load identification method as described in the first aspect.

[0019] The beneficial effects of the technical solutions provided in this application include: By dividing the preset time domain into several time periods in a multi-degree-of-freedom motion system model, each m consecutive time period is used as a dynamic load identification unit. In the multi-degree-of-freedom motion system model, there are N degrees of freedom. Within the preset time domain, a dynamic load is applied to each of the m degrees of freedom. The force values ​​of the m dynamic loads remain constant within an identification unit, and m is less than or equal to N. For each identification unit, a system of equations is solved based on the dynamic load force values ​​to calculate the dynamic load force value of each identification unit. This solves the technical problem in related technologies where it is difficult to accurately separate and identify each independent load when facing multi-source excitation problems. Attached Figure Description

[0020] Figure 1 This is a flowchart illustrating an embodiment of the multi-point excitation dynamic load identification method of this application; Figure 2 This is a schematic diagram of a three-mass spring-damped system; Figure 3 for The time-domain recognition result curve; Figure 4 for The time-domain recognition result curve; Figure 5 This is a functional module diagram of an embodiment of the multi-point excitation dynamic load identification device of this application; Figure 6 This is a schematic diagram of the hardware structure of the multi-point excitation dynamic load identification device involved in the embodiments of this application. Detailed Implementation

[0021] To enable those skilled in the art to better understand the present application, the technical solutions in the embodiments of the present application will be clearly and completely described below with reference to the accompanying drawings. Obviously, the described embodiments are only some embodiments of the present application, and not all embodiments. Based on the embodiments in the present application, all other embodiments obtained by those of ordinary skill in the art without creative effort are within the scope of protection of the present application.

[0022] To make the objectives, technical solutions, and advantages of this application clearer, the embodiments of this application will be described in further detail below with reference to the accompanying drawings.

[0023] In a first aspect, embodiments of this application provide a method for identifying multi-point excitation dynamic loads.

[0024] In one embodiment, reference is made to Figure 1 , Figure 1 This is a flowchart illustrating an embodiment of the multi-point excitation dynamic load identification method of this application. Figure 1 As shown, the multi-point excitation dynamic load identification method includes: Step S10: In the multi-degree-of-freedom motion system model, the preset time domain is divided into several time periods, and each m consecutive time periods is a dynamic load identification unit. In the multi-degree-of-freedom motion system model, there are N degrees of freedom. In the preset time domain, a dynamic load is applied to each of the m degrees of freedom. The force values ​​of the m dynamic loads are constant within an identification unit, and m is less than or equal to N. In one embodiment, in a multi-degree-of-freedom motion system model, the time period during which the dynamic load is applied is taken as a preset time domain T, and T is divided into a small time periods. ,but In a system with N degrees of freedom, if a dynamic load is applied to each of the m degrees of freedom, and m time intervals are considered as one identification unit (the number of dynamic loads equals the number of time intervals within the identification unit), then there are b identification units. At this point, the dynamic load within a single identification unit is assumed to remain constant. Theoretically speaking, The smaller the value, the more accurate the calculation.

[0025] For example, refer to Figure 2 , Figure 2 This is a schematic diagram of a three-mass spring-damped system. Figure 2 As shown, in a three-mass spring-damped system, there are three degrees of freedom. , and ,exist Apply dynamic load ,exist Apply dynamic load ,Will and The time period of action is taken as the preset time domain T, and so on. Taking a time interval of 0.01s as an example, assuming T = 10s, T is divided into 1000 time intervals. Since there are two dynamic loads in the system, each identification unit should contain two time intervals, resulting in a total of 500 identification units. In each identification unit, the external load is considered a constant value. For example, in the first identification unit... . and Constant and unchanging.

[0026] Furthermore, before calculating the dynamic load force value of each identification unit by solving the system of equations based on the dynamic load force value, the following steps are included: In a multi-degree-of-freedom motion system model, the displacement and velocity of each degree of freedom corresponding to each time period are collected. In one identification unit, there are m+1 time points.

[0027] In one embodiment, in this three-mass spring damping system, the first identification unit has two time periods and three time points. , exist - The displacement over a time period is ,exist - The displacement over a time period is , exist The velocities at time points are as follows: , , , exist - The displacement at time t is ,exist - The displacement over a time period is , exist The velocities at time points are as follows: , , ,in exist The velocities at time points are as follows: , , .

[0028] Step S20: For each identification unit, solve the equation set based on the dynamic load force value to calculate the dynamic load force value of each identification unit.

[0029] Furthermore, the equations for solving the dynamic load force value are as follows:

[0030] in, For the first Within the time period, the first i The work done by a dynamic load. For the first The increase in kinetic energy over a given period of time For the first The increase in potential energy over a given time period No. Damping energy dissipation over a time period, of which k The range of values ​​for is [1, m]. i The value range is [1, m].

[0031] Furthermore, the kinetic energy increment, potential energy increment, and damping dissipation energy during the time period include: The increase in kinetic energy within the time period is the difference in system kinetic energy at the beginning and end of the time period; The potential energy increment within the time period is the difference in system potential energy at the beginning and end of the time period. The damping dissipation energy within the time period refers to the energy consumed by damping within a time period.

[0032] In one embodiment, the velocity and mass of each degree of freedom corresponding to each time period, as well as the connection stiffness between the degrees of freedom, are utilized. , Based on the formulas for kinetic and potential energy, the kinetic and potential energy increments for each time period in each identification unit are obtained; using the damping coefficient... , , Based on the damping dissipation energy calculation formula, the damping dissipation energy of each identification unit in each time period is obtained; Will and Displacement in the first identification unit , Substituting the kinetic energy increment, potential energy increment, and damping dissipation energy obtained above for each time period in each identification unit into the dynamic load force value equation set, we get:

[0033] Solve the first recognition unit .

[0034] in, For the first The work done by the i-th dynamic load within a time period. For the first The increase in kinetic energy over a given period of time For the first The increase in potential energy over a given time period No. The damping dissipation energy over a time period, where k ranges from [1,2] and i ranges from [1,2]; for example... The first time period in the first identification unit The work done The increase in kinetic energy during the first time period. This represents the potential energy increment during the first time interval. Damping energy dissipation during the first time period.

[0035] Following the same solution method, the results for all recognition units are obtained. .

[0036] Further, after calculating each dynamic load force value of the identification unit, the process includes: Based on each dynamic load force value of each identification unit, a time-domain identification curve for each dynamic load within a preset time domain is obtained by fitting.

[0037] like Figure 3 , Figure 4 As shown, Figure 3 for The time-domain recognition result curve; Figure 4 for The temporal domain recognition result curve.

[0038] In this embodiment, the preset time domain is divided into several time periods in the multi-degree-of-freedom motion system model, and each m consecutive time period is used as a dynamic load identification unit. There are N degrees of freedom in the multi-degree-of-freedom motion system model. Within the preset time domain, a dynamic load is applied to each of the m degrees of freedom. The force values ​​of the m dynamic loads remain constant within an identification unit, and m is less than or equal to N. For each identification unit, the equation system is solved based on the dynamic load force value to calculate the dynamic load force value of each identification unit. This solves the technical problem in related technologies where it is difficult to accurately separate and identify each independent load when facing multi-source excitation problems. Secondly, embodiments of this application also provide a multi-point excitation dynamic load identification device.

[0039] In one embodiment, reference is made to Figure 5 , Figure 5 This is a functional module diagram of an embodiment of the multi-point excitation dynamic load identification device of this application. Figure 5 As shown, the multi-point excitation dynamic load identification device includes: Force application module 10: In the multi-degree-of-freedom motion system model, the preset time domain is divided into several time periods, and each m consecutive time period is a dynamic load identification unit. In the multi-degree-of-freedom motion system model, there are N degrees of freedom. In the preset time domain, a dynamic load is applied to each of the m degrees of freedom. The force values ​​of the m dynamic loads are constant within an identification unit, and m is less than or equal to N.

[0040] Calculation module 20: For each identification unit, solves the equation system based on the dynamic load force value to calculate the dynamic load force value of each identification unit.

[0041] Furthermore, in one embodiment, the multi-point excitation dynamic load identification device further includes a data acquisition module, used for: In a multi-degree-of-freedom motion system model, the displacement and velocity of each degree of freedom corresponding to each time period are collected. In a recognition unit, there are m+1 time points.

[0042] Furthermore, in one embodiment, the computing module 20 is used for:

[0043] in, For the first Within the time period, the first i The work done by a dynamic load. For the first The increase in kinetic energy over a given period of time For the first The increase in potential energy over a given time period No. Damping energy dissipation over a time period, of which k The range of values ​​for is [1, m]. i The value range is [1, m].

[0044] Furthermore, in one embodiment, the kinetic energy increment, the potential energy increment, and the damping dissipation energy during the time period include: The increase in kinetic energy within the time period is the difference in system kinetic energy at the beginning and end of the time period; The potential energy increment within the time period is the difference in system potential energy at the beginning and end of the time period. The damping dissipation energy within the time period refers to the energy consumed by damping within a time period.

[0045] Furthermore, in one embodiment, the multi-point excitation dynamic load identification device further includes a fitting module, used for: Based on each dynamic load force value of each identification unit, a time-domain identification curve for each dynamic load within a preset time domain is obtained by fitting.

[0046] The functions of each module in the multi-point excitation dynamic load identification device correspond to the steps in the multi-point excitation dynamic load identification method embodiment, and their functions and implementation processes will not be described in detail here.

[0047] Thirdly, embodiments of this application provide a multi-point excitation dynamic load identification device, which can be a personal computer (PC), laptop computer, server, or other device with data processing capabilities.

[0048] Reference Figure 6 , Figure 6 This is a schematic diagram of the hardware structure of the multi-point excitation dynamic load identification device involved in the embodiments of this application. In the embodiments of this application, the multi-point excitation dynamic load identification device may include a processor, a memory, a communication interface, and a communication bus.

[0049] The communication bus can be of any type and is used to interconnect the processor, memory, and communication interface.

[0050] The communication interface includes input / output (I / O) interfaces, physical interfaces, and logical interfaces used for interconnecting internal components of the multi-point excitation dynamic load identification device, as well as interfaces used for interconnecting the multi-point excitation dynamic load identification device with other devices (such as other computing devices or user equipment). Physical interfaces can be Ethernet interfaces, fiber optic interfaces, ATM interfaces, etc.; user equipment can be displays, keyboards, etc.

[0051] Memory can be various types of storage media, such as random access memory (RAM), read-only memory (ROM), non-volatile RAM (NVRAM), flash memory, optical storage, hard disk, programmable ROM (PROM), erasable PROM (EPROM), electrically erasable PROM (EEPROM), etc.

[0052] The processor can be a general-purpose processor, which can call the multi-point excitation dynamic load identification program stored in the memory and execute the multi-point excitation dynamic load identification method provided in the embodiments of this application. For example, the general-purpose processor can be a central processing unit (CPU). The method executed when the multi-point excitation dynamic load identification program is called can be referred to in the various embodiments of the multi-point excitation dynamic load identification method of this application, and will not be repeated here.

[0053] Those skilled in the art will understand that Figure 6 The hardware structure shown does not constitute a limitation of this application and may include more or fewer components than shown, or combine certain components, or have different component arrangements.

[0054] Fourthly, embodiments of this application also provide a computer-readable storage medium.

[0055] The present application stores a multi-point excitation dynamic load identification program on a computer-readable storage medium, wherein when the multi-point excitation dynamic load identification program is executed by a processor, it implements the steps of the multi-point excitation dynamic load identification method as described above.

[0056] The method implemented when the multi-point excitation dynamic load identification program is executed can be referred to in various embodiments of the multi-point excitation dynamic load identification method of this application, and will not be repeated here.

[0057] It should be noted that the sequence numbers of the embodiments in this application are for descriptive purposes only and do not represent the superiority or inferiority of the embodiments.

[0058] The terms "comprising" and "having," and any variations thereof, in the specification, claims, and accompanying drawings of this application are intended to cover non-exclusive inclusion. For example, a process, method, system, product, or apparatus that includes a series of steps or units is not limited to the listed steps or units, but may optionally include steps or units not listed, or may optionally include other steps or units inherent to such process, method, product, or apparatus. The terms "first," "second," and "third," etc., are used to distinguish different objects, etc., and do not indicate a sequence, nor do they limit "first," "second," and "third" to different types.

[0059] In the description of the embodiments of this application, terms such as "exemplary," "for example," or "for instance" are used to indicate examples, illustrations, or explanations. Any embodiment or design described as "exemplary," "for example," or "for instance" in the embodiments of this application should not be construed as being more preferred or advantageous than other embodiments or designs. Specifically, the use of terms such as "exemplary," "for example," or "for instance" is intended to present the relevant concepts in a concrete manner.

[0060] In the description of the embodiments of this application, unless otherwise stated, " / " means "or". For example, A / B can mean A or B. The "and / or" in the text is merely a description of the relationship between related objects, indicating that there can be three relationships. For example, A and / or B can mean: A exists alone, A and B exist simultaneously, and B exists alone. In addition, in the description of the embodiments of this application, "multiple" means two or more.

[0061] In some processes described in the embodiments of this application, multiple operations or steps are included in a specific order. However, it should be understood that these operations or steps may not be executed in the order they appear in the embodiments of this application, or they may be executed in parallel. The sequence number of the operation is only used to distinguish different operations, and the sequence number itself does not represent any execution order. In addition, these processes may include more or fewer operations, and these operations or steps may be executed sequentially or in parallel, and these operations or steps may be combined.

[0062] Through the above description of the embodiments, those skilled in the art can clearly understand that the methods of the above embodiments can be implemented by means of software plus necessary general-purpose hardware platforms. Of course, they can also be implemented by hardware, but in many cases the former is a better implementation method. Based on this understanding, the technical solution of this application, in essence, or the part that contributes to the prior art, can be embodied in the form of a software product. This computer software product is stored in a storage medium (such as ROM / RAM, magnetic disk, optical disk) as described above, and includes several instructions to cause a terminal device to execute the methods described in the various embodiments of this application.

[0063] The above are merely preferred embodiments of this application and do not limit the patent scope of this application. Any equivalent structural or procedural transformations made using the content of this application's specification and drawings, or direct or indirect applications in other related technical fields, are similarly included within the patent protection scope of this application.

Claims

1. A method for identifying multi-point excitation dynamic loads, characterized in that, The multi-point excitation dynamic load identification method includes: In the multi-degree-of-freedom motion system model, the preset time domain is divided into several time periods, and each m consecutive time periods is a dynamic load identification unit. In the multi-degree-of-freedom motion system model, there are N degrees of freedom. In the preset time domain, a dynamic load is applied to each of the m degrees of freedom. The force values ​​of the m dynamic loads are constant within an identification unit, and m is less than or equal to N. For each identification unit, the dynamic load force value of each identification unit is calculated by solving the system of equations based on the dynamic load force value.

2. The multi-point excitation dynamic load identification method as described in claim 1, characterized in that, Before calculating the dynamic load force value of each identification unit by solving the system of equations based on the dynamic load force value, the following steps are included: In a multi-degree-of-freedom motion system model, the displacement and velocity of each degree of freedom corresponding to each time period are collected. In a recognition unit, there are m+1 time points.

3. The multi-point excitation dynamic load identification method as described in claim 2, characterized in that, The equations for solving the dynamic load force value are as follows: in, For the first Within the time period, the first i The work done by a dynamic load. For the first The increase in kinetic energy over a given period of time For the first The increase in potential energy over a given time period No. Damping energy dissipation over a time period, of which k The range of values ​​for is [1, m]. i The value range is [1, m].

4. The multi-point excitation dynamic load identification method as described in claim 3, characterized in that, The kinetic energy increment, potential energy increment, and damping dissipation energy during the time period include: The increase in kinetic energy within the time period is the difference in system kinetic energy at the beginning and end of the time period; The potential energy increment within the time period is the difference in system potential energy at the beginning and end of the time period. The damping dissipation energy within the time period refers to the energy consumed by damping within a time period.

5. The multi-point excitation dynamic load identification method as described in claim 3, characterized in that, After calculating the dynamic load force value of each identification unit, the following is included: Based on each dynamic load force value of each identification unit, a time-domain identification curve for each dynamic load within a preset time domain is obtained by fitting.

6. A multi-point excitation dynamic load identification device, characterized in that, The multi-point excitation dynamic load identification device includes: Force application module: In the multi-degree-of-freedom motion system model, the preset time domain is divided into several time periods, and each m consecutive time period is a dynamic load identification unit. In the multi-degree-of-freedom motion system model, there are N degrees of freedom. Within the preset time domain, a dynamic load is applied to each of the m degrees of freedom. The force values ​​of the m dynamic loads remain constant within an identification unit, and m is less than or equal to N. Calculation module: For each identification unit, the system of equations is solved based on the dynamic load force value to calculate the dynamic load force value of each identification unit.

7. The multi-point excitation dynamic load identification device as described in claim 6, characterized in that, The computing module includes: in, For the first Within the time period, the first i The work done by a dynamic load. For the first The increase in kinetic energy over a given period of time For the first The increase in potential energy over a given time period No. Damping energy dissipation over a time period, of which k The range of values ​​for is [1, m]. i The value range is [1, m].

8. The multi-point excitation dynamic load identification device as described in claim 7, characterized in that, The kinetic energy increment, potential energy increment, and damping dissipation energy during the time period include: The increase in kinetic energy within the time period is the difference in system kinetic energy at the beginning and end of the time period; The potential energy increment within the time period is the difference in system potential energy at the beginning and end of the time period. The damping dissipation energy within the time period refers to the energy consumed by damping within a time period.

9. A multi-point excitation dynamic load identification device, characterized in that, The multi-point excitation dynamic load identification device includes a processor, a memory, and a multi-point excitation dynamic load identification program stored in the memory and executable by the processor, wherein when the multi-point excitation dynamic load identification program is executed by the processor, it implements the steps of the multi-point excitation dynamic load identification method as described in any one of claims 1 to 5.

10. A computer-readable storage medium, characterized in that, The computer-readable storage medium stores a multi-point excitation dynamic load identification program, wherein when the multi-point excitation dynamic load identification program is executed by a processor, it implements the steps of the multi-point excitation dynamic load identification method as described in any one of claims 1 to 5.

Images