Bushing dynamic stiffness test method, device, test equipment, system and medium

Abstract

The invention discloses a bushing dynamic stiffness test method, device, test equipment, system and medium. The method comprises the following steps: exciting two subsystems, and determining an initial transfer function matrix; reducing the initial transfer function matrix, and determining an original point transfer function matrix corresponding to the two subsystems and an original point transfer function matrix corresponding to the coupling system; and determining the dynamic stiffness of the bushing based on the original point transfer function matrixes corresponding to the two subsystems and the original point transfer function matrix corresponding to the coupling system. According to the method, the dynamic stiffness of the bushing is tested through the coupling system connected with the bushing, so that the problem that the dynamic stiffness of the bushing is inaccurately tested in a high-frequency range due to resonance and other problems of existing test equipment can be effectively avoided, and the relatively accurate and stable dynamic stiffness test can be performed on the bushing in a relatively large frequency range.

Description

Technical Field

[0001] The present invention relates to the field of bushing testing, and in particular to a method, device, test equipment, system and medium for testing the dynamic stiffness of a bushing. Background Art

[0002] In the prior art, the dynamic stiffness of a bushing is mainly tested by relying on an MTS (Mechanical Test Systems) test system. Due to problems such as resonance of the test equipment, the frequency range for the MTS to test the dynamic stiffness of the bushing is small, and it is impossible to accurately test the dynamic stiffness of the bushing within a large frequency range (especially within a high-frequency range). Therefore, the technical problem to be solved currently is how to accurately test the dynamic stiffness of the bushing within a large frequency range. Summary of the Invention

[0003] Embodiments of the present invention provide a method, device, test equipment, system and medium for testing the dynamic stiffness of a bushing to solve the technical problem of how to accurately test the dynamic stiffness of the bushing within a large frequency range.

[0004] A method for testing the dynamic stiffness of a bushing is applicable to a coupled system, where the coupled system includes two subsystems, and the two subsystems are respectively connected to both ends of the bushing. The method for testing the dynamic stiffness of the bushing includes: Exciting the two subsystems to determine an initial transfer function matrix; the initial transfer function matrix is a matrix formed by K*K initial transfer function values determined based on K excitation data and K response data corresponding to each excitation data. Reducing the initial transfer function matrix to determine the origin transfer function matrix corresponding to the two subsystems and the origin transfer function matrix corresponding to the coupled system; the origin transfer function matrix corresponding to each subsystem is a matrix formed by L*L origin transfer function values determined based on L excitation data corresponding to the subsystem and L response data corresponding to each excitation data in the same subsystem; the origin transfer function matrix corresponding to the coupled system is a matrix formed by L*L origin transfer function values determined based on L excitation data corresponding to any one of the subsystems and L response data corresponding to each excitation data in the other subsystem, where 2 < L < K. Determining the dynamic stiffness of the bushing based on the origin transfer function matrix corresponding to the two subsystems and the origin transfer function matrix corresponding to the coupled system.

[0005] The above-mentioned bushing dynamic stiffness testing method determines the initial transfer function matrix by exciting two subsystems of the coupled system, reduces the initial transfer function matrix generated by the excitation, obtains the reduced origin transfer function matrices corresponding to the two subsystems and the reduced origin transfer function matrix corresponding to the coupled system, and determines the dynamic stiffness of the bushing connecting the coupled system according to the reduced origin transfer function matrices corresponding to the two subsystems and the reduced origin transfer function matrix corresponding to the coupled system. This method determines the dynamic stiffness of the bushing according to the reduced origin transfer function matrix, can effectively reduce the amount of data processing, and can relatively quickly realize the dynamic stiffness testing of the bushing connecting the coupled system. This method tests the dynamic stiffness of the bushing through the coupled system connected to the bushing, can effectively avoid the problem that the existing test equipment causes inaccurate dynamic stiffness testing of the bushing in the high-frequency range due to resonance and other problems, and can perform relatively accurate and stable dynamic stiffness testing of the bushing in a large frequency range.

[0006] Preferably, each of the subsystems includes N indicating points, and the indicating points are excitation points and response points; The exciting the two subsystems to determine the initial transfer function matrix includes: exciting the two subsystems to determine K excitation data and K response data corresponding to each excitation data. The K excitation data are data formed by exciting 2*N excitation points along M excitation directions respectively, and the K response data are data formed by the 2*N response points responding in the M response directions for each excitation. K = 2*N*M, N≥2, M≥1; based on an excitation data and a response data corresponding to the excitation data, determining an initial transfer function value, and based on K*K initial transfer function values, determining the initial transfer function matrix. In this embodiment, the initial transfer function matrix is determined to determine the dynamic stiffness of the bushing connected to the two subsystems according to the initial transfer function matrix, so as to realize the dynamic stiffness testing of the bushing in a large frequency range.

[0007] Preferably, the reducing the initial transfer function matrix to determine the origin transfer function matrices corresponding to the two subsystems and the origin transfer function matrix corresponding to the coupled system includes: reducing the initial transfer function matrix to determine the target transfer function matrices corresponding to the two subsystems and the target transfer function matrix corresponding to the coupled system; the target transfer function matrix corresponding to each subsystem is a matrix formed by R*R initial transfer function values determined based on R excitation data corresponding to the subsystem and R response data corresponding to each excitation data in the same subsystem; the origin transfer function matrix corresponding to the coupled system is a matrix formed by R*R initial transfer function values determined based on R excitation data corresponding to any one of the subsystems and R response data corresponding to each excitation data in the other subsystem, L < R = Reduce the target transfer function matrices corresponding to the two subsystems and the target transfer function matrix corresponding to the coupled system based on the excitation point coordinates corresponding to the excitation data and the response point coordinates corresponding to the response data, and respectively determine the origin transfer function matrices corresponding to the two subsystems and the origin transfer function matrix corresponding to the coupled system.

[0008] In this embodiment, obtain the origin transfer function matrices corresponding to the two subsystems and the origin transfer function matrix corresponding to the coupled system, so as to accurately and stably determine the dynamic stiffness of the bushing according to the origin transfer function matrix subsequently.

[0009] Preferably, the reducing the target transfer function matrices corresponding to the two subsystems and the target transfer function matrix corresponding to the coupled system based on the excitation point coordinates corresponding to the excitation data and the response point coordinates corresponding to the response data, and respectively determining the origin transfer function matrices corresponding to the two subsystems and the origin transfer function matrix corresponding to the coupled system includes: Determine the excitation reduction matrix L* corresponding to each subsystem based on the excitation point coordinates corresponding to the excitation data of each subsystem and the origin coordinates corresponding to each subsystem; Based on the response point coordinates corresponding to the response data of each subsystem and the origin coordinates corresponding to each subsystem, determine the response reduction matrix *L corresponding to each subsystem; Use the excitation reduction matrices corresponding to the two subsystems and the response reduction matrices corresponding to the two subsystems to reduce the target transfer function matrix of the target system, and determine the origin transfer function matrix corresponding to the target system; the target system is the two subsystems or the coupled system.

[0010] In this embodiment, reduce the excitation points and response points in the target system to the origin corresponding to the target system, and reduce the target initial transfer function matrix to the origin transfer function matrix corresponding to the origin, so as to make it feasible to determine the dynamic stiffness of the bushing according to the origin transfer function matrix.

[0011] Preferably, the origin transfer function matrix corresponding to the target system is the product of the inverse matrix of the response reduction matrix corresponding to the target system, the target transfer function matrix, and the inverse matrix of the excitation reduction matrix. In this embodiment, obtain the origin transfer function matrices corresponding to the first subsystem, the second subsystem, and the coupled system, which are used to accurately and stably determine the dynamic stiffness of the bushing within a large frequency range (especially within the high-frequency range).

[0012] Preferably, the dynamic stiffness of the bushing is the quotient of the origin transfer function matrix corresponding to the coupling system and the target determinant; the target determinant is the difference between the first matrix product and the second matrix product, the first matrix product is the product of the origin transfer function matrices corresponding to two subsystems, and the second matrix product is the product of the origin transfer function matrix corresponding to the coupling system and the origin transfer function matrix corresponding to the coupling system. In this embodiment, according to the relationship between the deformation of the bushing relative to each of the two subsystems and the relationship between the equivalent excitation of each subsystem and the deformation of the bushing, the origin transfer function matrix is processed to accurately determine the dynamic stiffness of the bushing and realize the test of the dynamic stiffness of the bushing.

[0013] A bushing dynamic stiffness testing device includes: An initial transfer function matrix determination module, configured to excite two subsystems to determine an initial transfer function matrix; the initial transfer function matrix is a matrix formed by K*K initial transfer function values determined based on K excitation data and K response data corresponding to each of the excitation data. An origin transfer function matrix determination module, configured to reduce the initial transfer function matrix to determine the origin transfer function matrices corresponding to two subsystems and the origin transfer function matrix corresponding to the coupling system; each origin transfer function matrix corresponding to a subsystem is a matrix formed by L*L origin transfer function values determined based on L excitation data corresponding to the subsystem and L response data corresponding to each of the excitation data in the same subsystem; the origin transfer function matrix corresponding to the coupling system is a matrix formed by L*L origin transfer function values determined based on L excitation data corresponding to any one of the subsystems and L response data corresponding to each of the excitation data in the other subsystem, where L < K. A dynamic stiffness determination module, based on the origin transfer function matrices corresponding to two subsystems and the origin transfer function matrix corresponding to the coupling system, determines the dynamic stiffness of the bushing.

[0014] A testing device includes a memory, a processor, and a computer program stored in the memory and executable on the processor. When the processor executes the computer program, the above-mentioned bushing dynamic stiffness testing method is implemented.

[0015] A bushing dynamic stiffness testing system includes the above-mentioned testing device and a coupling system connected to the testing device. The coupling system includes two subsystems, and the two subsystems are respectively connected to both ends of the bushing.

[0016] A computer-readable storage medium stores a computer program, and when the computer program is executed by a processor, the above-mentioned bushing dynamic stiffness testing method is implemented. Description of the Drawings

[0017] To more clearly illustrate the technical solutions of the embodiments of the present invention, the following briefly introduces the drawings required for the description of the embodiments of the present invention. Obviously, the drawings in the following description are only some embodiments of the present invention. For those of ordinary skill in the art, without creative efforts, other drawings can be obtained based on these drawings.

[0018] Figure 1 is a flowchart of a bushing dynamic stiffness testing method in an embodiment of the present invention; Figure 2 is a schematic diagram of a bushing dynamic stiffness testing device in an embodiment of the present invention; Figure 3 is a schematic diagram of the excitation point and response point of a coupling system connected to a bushing in an embodiment of the present invention; Figure 4 is an image of the initial transfer function values of a response point in the X-axis, Y-axis, and Z-axis directions within a partial frequency range when an excitation is applied to an excitation point in the Z-axis direction in an embodiment of the present invention; Figure 5 is an image of the test result of the bushing dynamic stiffness within the frequency range of [0, 3000 Hz] in an embodiment of the present invention; In the figure, A, the first subsystem; B, the second subsystem; C, the bushing. Detailed implementation manners

[0019] The following will clearly and completely describe the technical solutions in the embodiments of the present invention with reference to the drawings in the embodiments of the present invention. Obviously, the described embodiments are some, but not all, of the embodiments of the present invention. All other embodiments obtained by those of ordinary skill in the art based on the embodiments of the present invention without creative efforts fall within the protection scope of the present invention.

[0020] The bushing dynamic stiffness testing method provided by the embodiments of the present invention is used to accurately test the dynamic stiffness of the bushing within a relatively large frequency range.

[0021] In one embodiment, the bushing dynamic stiffness test method is applicable to a coupled system. The coupled system includes two subsystems, and the two subsystems are respectively connected to both ends of the bushing for dynamically testing the bushing through the two subsystems. In this example, the two subsystems are respectively defined as the first subsystem and the second subsystem. The first subsystem and the second subsystem are the two subsystems of the coupled system, and the first subsystem and the second subsystem are respectively connected to the bushing. Each of the subsystems includes N indicating points, and the two subsystems are provided with 2*N indicating points. These 2*N indicating points can be either excitation points or response points. In this example, when any indicating point is excited, this indicating point is an excitation point, and for this excitation, all indicating points can be their corresponding response points.

[0022] In one embodiment, as Figure 1 shown, a bushing dynamic stiffness test method is provided. Taking the application of this method in a test device as an example, the method includes the following steps: S101: Excite the two subsystems to determine the initial transfer function matrix; the initial transfer function matrix is a matrix formed by K*K initial transfer function values determined based on K excitation data and K response data corresponding to each excitation data; S102: Reduce the initial transfer function matrix to determine the origin transfer function matrix corresponding to the two subsystems and the origin transfer function matrix corresponding to the coupled system; the origin transfer function matrix corresponding to each subsystem is a matrix formed by L*L origin transfer function values determined based on L excitation data corresponding to the subsystem and L response data corresponding to each excitation data in the same subsystem; the origin transfer function matrix corresponding to the coupled system is a matrix formed by L*L origin transfer function values determined based on L excitation data corresponding to any one subsystem and L response data corresponding to each excitation data in the other subsystem, where 2 < L < K; S103: Determine the dynamic stiffness of the bushing based on the origin transfer function matrix corresponding to the two subsystems and the origin transfer function matrix corresponding to the coupled system.

[0023] Among them, the initial transfer function matrix is a matrix formed by K*K initial transfer function values determined based on K excitation data and K response data corresponding to each excitation data, where K > 2. The excitation data refers to the data used to excite the excitation points of the subsystem. The response data refers to the data detected at the response points when the excitation points are excited. The excitation data in this example is the excitation force applied to the excitation points by a force hammer, and the response data is the acceleration detected at the response points. The initial transfer function value refers to the transfer function value between the excitation point and the response point when the excitation point on the subsystem is excited and a response is received at the response point.

[0024] As an example, in step S101, during the dynamic stiffness test of the bushing, the two subsystems can be excited separately. When the excitation points of each subsystem are excited according to the excitation data, the response data are obtained at the response points of each subsystem in sequence. According to the excitation data of each excitation point and the response data of each response point, the initial transfer function values between each excitation point and each response point are determined, and the initial transfer function values are arranged in matrix form to obtain the initial transfer function matrix. In this example, at each frequency in a preset large frequency range, a force hammer is used to excite one or more excitation directions of each excitation point in the two subsystems according to the excitation data, and the response data corresponding to each excitation data are obtained at each response point in the two subsystems. According to the response data and the excitation data, K*K initial transfer function values are determined, and the K*K initial transfer function values are arranged with K excitation data as columns and the K response data corresponding to each excitation data as rows to obtain the K*K initial transfer function matrix corresponding to each frequency. At this time, in the initial transfer function matrix, the initial transfer function value in the i-th row and the j-th column represents the transfer function value determined by the response data of the j-th response point and the excitation data of the i-th excitation point when the i-th excitation point is excited by the excitation data corresponding thereto. It can be understood that it can also be arranged with K excitation data as rows and the K response data corresponding to each excitation data as columns, and its setting method can be determined independently according to the actual situation. At this time, in the initial transfer function matrix, the initial transfer function value in the i-th row and the j-th column represents the transfer function value determined by the response data of the i-th response point and the excitation data of the j-th excitation point when the j-th excitation point is excited by the excitation data corresponding thereto.

[0025] Among them, the origin transfer function matrix is the matrix after reducing the initial transfer function matrix. The origin transfer function value refers to the transfer function value between the equivalent excitation and the equivalent response. The equivalent excitation refers to the excitation when the excitation data acts equivalently on the origin, and the equivalent response refers to the response when the response data acts equivalently on the origin. The origin refers to the equivalent point after reducing the excitation points and response points in the two subsystems. In this example, the origin is the center point of the bushing connected to the two subsystems.

[0026] As an example, in step S102, the steps for the test device to determine the origin transfer function matrix corresponding to the two subsystems and the origin transfer function matrix corresponding to the coupled system include: (1) The steps for determining the origin transfer function matrix corresponding to any one of the two subsystems include: First, the test device reduces the initial transfer function matrix to determine the transfer function matrix corresponding to the subsystem. Among them, the transfer function matrix corresponding to the subsystem refers to the matrix formed by the transfer function values between the excitation point and the response point when the excitation point and the response point are in the same subsystem.

[0027] Secondly, the test equipment equivalentizes the excitation points and response points in the subsystem to the corresponding origin of the subsystem, reduces the excitation data of multiple excitation points to L excitation data at the origin, and reduces the response data of multiple response points to L response data at the origin. The operation relationship among multiple excitation points, multiple response points, and the transfer function matrix corresponding to the subsystem is represented by using the L excitation data at the origin, the L response data at the origin, and the transfer function matrix corresponding to the subsystem. In this example, the origin of the subsystem is the center point of the bushing connected to the subsystem.

[0028] Finally, process the above operation relationship to determine the processing result corresponding to the ratio of the L excitation data at the origin and the L response data at the origin, and this processing result is the origin transfer function matrix corresponding to the subsystem.

[0029] (2) Determining the origin transfer function matrix corresponding to the coupled system includes the following steps: First, the test equipment reduces the initial transfer function matrix to determine the transfer function matrix corresponding to the coupled system. Among them, the transfer function matrix corresponding to the coupled system refers to the matrix formed by the transfer function values between the excitation point and the response point when the excitation point and the response point are respectively in two subsystems.

[0030] Secondly, the test equipment equivalentizes the excitation points in one subsystem and the response points in another subsystem to the origin corresponding to the coupled system, reduces the excitation data of multiple excitation points to L excitation data at the origin, and reduces the response data of multiple response points to L response data at the origin. The operation relationship among multiple excitation points, multiple response points, and the transfer function matrix corresponding to the coupled system is represented by using the L excitation data at the origin, the L response data at the origin, and the transfer function matrix corresponding to the subsystem; in this example, the origin of the coupled system is the center point of the bushing connected to the coupled system.

[0031] Finally, process the above operation relationship to determine the processing result corresponding to the ratio of the L excitation data at the origin and the L response data at the origin, and this processing result is the origin transfer function matrix corresponding to the subsystem.

[0032] In this example, by reducing the initial transfer function matrices corresponding to the coupling system and the two subsystems in the coupling system, the K*K initial transfer function values are reduced to the L*L origin transfer function matrices corresponding to the two subsystems and the L*L origin transfer function matrix corresponding to the coupling system. The reduced L*L origin transfer function matrix can be directly used to determine the dynamic stiffness of the bushings connected to the two subsystems. The testing device reduces the initial transfer function matrix corresponding to each frequency in the preset large frequency range through the above method, and determines the origin transfer function matrices corresponding to the two subsystems and the origin transfer function matrix corresponding to the coupling system corresponding to each frequency, so as to determine the dynamic stiffness of the bushings corresponding to each frequency in the preset large frequency range according to the reduced origin transfer function matrix.

[0033] As an example, in step S103, for each frequency in the preset frequency range, the testing device can use the dynamic stiffness algorithm to process the several inputs of the L*L origin transfer function matrices corresponding to the two subsystems and the L*L origin transfer function matrix corresponding to the coupling system, and determine the dynamic stiffness of the bushings connected to the two subsystems at each frequency.

[0034] In this example, the preset frequency range can be [0, 3000Hz]. Figure 5 Show the test result image of the bushing dynamic stiffness in the frequency range of [0, 3000Hz]. Figure 5 The first curve image in shows the response data of a response direction of the origin in the coupling system at different frequencies at different frequencies. Figure 5 The second curve image in is the dynamic stiffness of the bushing at different frequencies. From Figure 5 It can be seen that this dynamic stiffness testing method can perform dynamic stiffness testing on bushings in a large frequency range, especially in the frequency range greater than 1000Hz. Compared with the prior art that can only perform dynamic stiffness testing on bushings in the frequency range below 1000Hz, this method can perform dynamic stiffness testing on bushings in a large frequency range, especially in the high-frequency range above 1000Hz, and has a wide frequency testing range.

[0035] In this embodiment, by exciting two subsystems of the coupling system, an initial transfer function matrix is determined, and the initial transfer function matrix generated by the excitation is reduced to obtain a reduced origin transfer function matrix corresponding to the two subsystems and a reduced origin transfer function matrix corresponding to the coupling system. Then, based on the reduced origin transfer function matrices corresponding to the two subsystems and the reduced origin transfer function matrix corresponding to the coupling system, the dynamic stiffness of the bushing connecting the coupling system is determined. This method determines the dynamic stiffness of the bushing based on the reduced origin transfer function matrix, which can effectively reduce the amount of data processing and relatively quickly realize the dynamic stiffness test of the bushing connecting the coupling system. By using the coupling system connected to the bushing to perform the dynamic stiffness test on the bushing, this method can effectively avoid the problem that the existing test equipment causes inaccurate dynamic stiffness test of the bushing in the high-frequency range due to resonance and other problems, and can perform relatively accurate and stable dynamic stiffness test on the bushing in a large frequency range.

[0036] In one embodiment, each subsystem includes N indicating points, and the indicating points are excitation points and response points. Herein, the indicating point refers to a point in the coupling system used for the dynamic stiffness test of the bushing connecting the coupling system. As an example, each subsystem in the coupling system includes N indicating points, and the same indicating point can be used as the excitation point of the coupling system or as the response point of the coupling system. For example, as Figure 3 shown, in the first subsystem A, the indicating points numbered 1 to 4 can be used as excitation points or as response points. In the second subsystem B, the indicating points numbered 5 to 8 can be used as excitation points or as response points.

[0037] In one embodiment, step S101, that is, exciting the two subsystems to determine the initial transfer function matrix, includes: S201: Excite the two subsystems to determine K excitation data and K response data corresponding to each excitation data. The K excitation data are data formed by exciting 2*N excitation points along M excitation directions respectively, and the K response data are data formed by the 2*N response points responding in the M response directions for each excitation. K = 2*N*M, N≥2, M≥1; S202: Based on each excitation data and a response data corresponding to the excitation data, determine an initial transfer function value, and based on the K*K initial transfer function values, determine the initial transfer function matrix.

[0038] As an example, in step S201, the test device excites N excitation points of each subsystem in M directions according to the method of exciting one excitation point in one excitation direction each time, determines K excitation data for excitation and K response data corresponding to each excitation data, where K = 2 * N * M, N ≥ 2, and M ≥ 1. In this example, each subsystem includes N indication points, which are both excitation points and response points. Each excitation point and each response point correspond to M excitation directions. The coupled system includes two subsystems, and the two subsystems include 2 * N excitation points and 2 * N response points. The test device excites one excitation direction of one excitation point each time according to the excitation data, and obtains the response data of the 2 * N response points in the two subsystems in M response directions for each excitation data, that is, obtains K response data corresponding to each excitation data, where K = 2 * N * M. The test device excites the 2 * N excitation points in M directions respectively according to the above method, determines the response data corresponding to each excitation data in the M directions of the 2 * N response points, and determines K excitation data and K * K response data corresponding to the K excitation data.

[0039] For example Figure 3 For the coupled system shown, the first subsystem A includes 4 indication points numbered from 1 to 4, and the second subsystem B includes 4 indication points numbered from 5 to 8. For the first subsystem A and the second subsystem B, N = 4. If each excitation point is excited in the X-axis direction, Y-axis direction, and Z-axis direction respectively, and the response data of each response point in the X-axis direction, Y-axis direction, and Z-axis direction is obtained, at this time, M = 3. In the coupled system corresponding to the first subsystem A and the second subsystem B, the excitation data is K = 2 * N * M = 24, and the response data corresponding to each excitation data is K = 2 * N * M = 24. The test device excites the 8 excitation points in the first subsystem A and the second subsystem B in the X-axis direction, Y-axis direction, and Z-axis direction respectively according to the K excitation data, and determines the response data of the 8 response points in the first subsystem A and the second subsystem B in the X-axis direction, Y-axis direction, and Z-axis direction relative to the response data corresponding to each excitation data, obtaining 24 * 24 response data. As shown in Table 1 below, the columns in Table 1 represent the response points and the corresponding response directions of the first subsystem A and the second subsystem B. For example, 1X in the column means the response point is numbered 1 and the response direction is the X-axis direction. The rows in Table 1 represent the excitation points and the excitation directions, and 1X in the row means the excitation point is numbered 1 and the excitation direction is the X-axis direction.

[0040] Table 1 As an example, in step S202, after the test device determines each excitation data and a response data corresponding to the excitation data, it can process the two input parameters of the excitation data and the corresponding response data based on a preset transfer function formula to determine the initial transfer function value between the two. For example, the transfer function formula can be set as a formula for calculating the ratio of the response data to its corresponding excitation data. Then, by processing the ratio of the K response data to each excitation data, K initial transfer function values corresponding to each excitation data are determined. The test device processes the K excitation data and the K*K response data corresponding to the K excitation data in the above manner to determine K*K initial transfer function values, arranges the initial transfer function values in K rows and K columns, and determines a K*K initial transfer function matrix.

[0041] As Figure 3 shown, M = 3, N = 4, K = 2*N*M = 24. For example, if the excitation point is No. 1 and the excitation direction is the X-axis direction, the response points can be No. 1 to No. 8. The response direction of each response point is the X-axis direction, the Y-axis direction, and the Z-axis direction. For one excitation point, each response point corresponds to a response data in the X-axis direction, the Y-axis direction, and the Z-axis direction respectively. Therefore, for the excitation data corresponding to the excitation point with the label 1 and the excitation direction being the X-axis direction, there are 24 response data. The ratio of these 24 response data to their corresponding excitation data can be processed to obtain 24 initial transfer function values corresponding to the 24 response data when the excitation point is No. 1 and the excitation direction is the X-axis direction (1X). , , , , , , . The test device determines 24*24 initial transfer function values corresponding to the excitation points and response points in the X-axis direction, the Y-axis direction, and the Z-axis direction when the labels are respectively from No. 1 to No. 8 and the response points are respectively from No. 1 to No. 8 according to the above method, and obtains the initial transfer function values between 24 excitation data and 24 response data in Table 1. , , . Then, as shown in Table 1, determine a 24*24 initial transfer function matrix corresponding to No. 1 to No. 8 in the X-axis, Y-axis, and Z-axis directions. , where the row represents the excitation point and the excitation direction of the excitation point, the column represents the response point and the response direction of the response point, and the initial transfer function matrix is: = where, Denotes the initial transfer function value determined by the excitation data with the excitation point labeled 1 and the excitation direction along the X-axis, and the response data with the response point labeled 1 and the response direction along the X-axis. Denotes the initial transfer function value determined by the excitation data with the excitation point labeled 1 and the excitation direction along the X-axis, and the response data with the response point labeled 8 and the response direction along the Z-axis. Denotes the initial transfer function value determined by the excitation data with the excitation point labeled 8 and the excitation direction along the Z-axis, and the response data with the response point labeled 1 and the response direction along the X-axis. Denotes the initial transfer function value determined by the excitation data with the excitation point labeled 8 and the excitation direction along the Z-axis, and the response data with the response point labeled 8 and the response direction along the Z-axis. As Figure 4 shown, it shows the initial transfer function values corresponding to each frequency in a certain frequency range when the excitation point is labeled 1 and the excitation direction is along the Z-axis, and the response point is 2 and the response directions are along the X-axis direction , the initial transfer function values in the Y-axis direction and the initial transfer function values in the Z-axis direction .

[0042] In this example, the test device can, at each frequency within a preset frequency range, excite N excitation points of each subsystem in two subsystems along M excitation directions respectively, determine the K*K initial transfer function matrix corresponding to each frequency, and determine the dynamic stiffness of the bushing connected to the two subsystems based on the initial transfer function matrix, so as to realize the dynamic stiffness test of the bushing within a relatively large frequency range.

[0043] In one embodiment, step S102, that is, reducing the initial transfer function matrix to determine the origin transfer function matrix corresponding to the two subsystems and the origin transfer function matrix corresponding to the coupled system, includes: S301: Reduce the initial transfer function matrix to determine the target transfer function matrix corresponding to the two subsystems and the target transfer function matrix corresponding to the coupled system; the target transfer function matrix corresponding to each subsystem is a matrix formed by R*R initial transfer function values determined based on R excitation data corresponding to the subsystem and R response data corresponding to each excitation data in the same subsystem; the origin transfer function matrix corresponding to the coupled system is a matrix formed by R*R initial transfer function values determined based on R excitation data corresponding to any one subsystem and R response data corresponding to each excitation data in the other subsystem, L<R= ; S302: Based on the excitation point coordinates corresponding to the excitation data and the response point coordinates corresponding to the response data, reduce the target transfer function matrices corresponding to the two subsystems and the target transfer function matrix corresponding to the coupled system, and respectively determine the origin transfer function matrices corresponding to the two subsystems and the origin transfer function matrix corresponding to the coupled system.

[0044] As an example, in step S301, after the test device obtains the initial transfer function matrix, it can reduce the initial transfer function matrix according to whether the excitation point and the response point are in the same system in the coupled system, and obtain the target transfer function matrices corresponding to the two subsystems and the target transfer function matrix corresponding to the coupled system.

[0045] In this example, the target transfer function matrix of the first subsystem is a matrix corresponding to the first transfer function values determined by the excitation data of N excitation points in M excitation directions in the first subsystem and the response data of N response points in M response directions. The first transfer function value refers to the initial transfer function value between the response data of N response points in M response directions in the first subsystem when the N excitation points in the first subsystem are excited by the excitation data corresponding to M excitation directions. It can be understood that there are R excitation data for N excitation points in M excitation directions in the first subsystem, and there are R response data for N response points in M response directions in the first subsystem, where R = N * M = , therefore, the target transfer function matrix of the first subsystem is a matrix formed by R * R initial transfer function values determined based on the R excitation data corresponding to the first subsystem and each excitation data in the R response data corresponding to the first subsystem.

[0046] As an example, as Figure 3 shown, the bushing C is respectively connected to the first subsystem A and the second subsystem B. The 12 excitation data of the excitation points 1 to 4 (N = 4) in the first subsystem A in the X-axis direction, Y-axis direction, and Z-axis direction (M = 3) and the 12 response data of the response points 1 to 4 in the first subsystem A in the X-axis direction, Y-axis direction, and Z-axis direction determine 12 * 12 first transfer function values, R = 12, and the first transfer function values are sorted in 12 rows and 12 columns to obtain the 12 * 12 target transfer function matrix of the first subsystem A , the target transfer function matrix of the first subsystem A is: = where, represents the first transfer function value determined by the excitation data with the excitation point numbered 1 and the excitation direction in the X-axis direction in the first subsystem A and the response data with the response point numbered 1 and the response direction in the X-axis direction. It represents the first transfer function value determined by the excitation data with the excitation point labeled 1 and the excitation direction along the X-axis in the first subsystem A, and the response data with the response point labeled 4 and the response direction along the Z-axis. It represents the first transfer function value determined by the excitation data with the excitation point labeled 4 and the excitation direction along the Z-axis in the first subsystem A, and the response data with the response point labeled 1 and the response direction along the X-axis. It represents the first transfer function value determined by the excitation data with the excitation point labeled 4 and the excitation direction along the Z-axis in the first subsystem A, and the response data with the response point labeled 4 and the response direction along the Z-axis.

[0047] The target transfer function matrix of the second subsystem is a matrix corresponding to the second transfer function values determined by the excitation data of N excitation points in M excitation directions and the response data of N response points in M response directions in the second subsystem. The second transfer function value refers to the initial transfer function value between the response data of N response points in M response directions in the second subsystem when the N excitation points in the second subsystem are excited by the excitation data corresponding to M excitation directions. It can be understood that there are R excitation data for N excitation points in M excitation directions in the second subsystem, and there are R response data for N response points in M response directions in the second subsystem, where R = N * M = Therefore, the target transfer function matrix of the second subsystem is a matrix formed by R * R initial transfer function values based on the corresponding R excitation data of the second subsystem and the corresponding R response data for each excitation data in the second subsystem.

[0048] As an example, as Figure 3 shown, 12 excitation data of excitation points 5 to 8 (N = 4) in the X-axis direction, Y-axis direction, and Z-axis direction (M = 3) in the second subsystem B and 12 response data of response points 5 to 8 in the X-axis direction, Y-axis direction, and Z-axis direction in the second subsystem B determine 12 * 12 second transfer function values. R = 12, and the second transfer function values are sorted in 12 rows and 12 columns to obtain the 12 * 12 target transfer function matrix of the second subsystem B The target transfer function matrix of the second subsystem B is: = Among them, It represents the second transfer function value determined by the excitation data with the excitation point labeled 5 and the excitation direction along the X-axis in the second subsystem B, and the response data with the response point labeled 5 and the response direction along the Z-axis. It represents the second transfer function value determined by the excitation data with the excitation point labeled 5 and the excitation direction along the X-axis in the second subsystem B, and the response data with the response point labeled 8 and the response direction along the Z-axis. It represents the second transfer function value determined by the excitation data with the excitation point labeled 8 and the excitation direction along the Z-axis in the second subsystem B, and the response data with the response point labeled 5 and the response direction along the X-axis. It represents the second transfer function value determined by the excitation data with the excitation point labeled 8 and the excitation direction along the Z-axis in the second subsystem B, and the response data with the response point labeled 8 and the response direction along the Z-axis.

[0049] The target transfer function matrix of the coupled system is the matrix corresponding to the first transfer function values determined by the excitation data of N excitation points in M excitation directions in the first subsystem and the response data of N response points in M response directions in the second subsystem. The third transfer function value refers to the initial transfer function value between the response data of N response points in M response directions in the second subsystem and the excitation data corresponding to the excitation of N excitation points in M excitation directions in the first subsystem when the N excitation points in the first subsystem are excited by the excitation data corresponding to M excitation directions. Alternatively, the third transfer function value refers to the initial transfer function value between the response data of N response points in M response directions in the first subsystem and the excitation data corresponding to the excitation of N excitation points in M excitation directions in the second subsystem when the N excitation points in the second subsystem are excited by the excitation data corresponding to M excitation directions. It can be understood that there are R excitation data for N excitation points in M excitation directions in the first subsystem, R response data for N response points in M response directions in the first subsystem, R excitation data for N excitation points in M excitation directions in the second subsystem, and R response data for N response points in M response directions in the second subsystem, where R = N * M = , therefore, the target transfer function matrix of the coupled system is a matrix formed by R * R initial transfer function values based on the R excitation data corresponding to the first subsystem and the R response data corresponding to each excitation data in the second subsystem, or the target transfer function matrix of the coupled system is a matrix formed by R * R initial transfer function values based on the R excitation data corresponding to the second subsystem and the R response data corresponding to each excitation data in the first subsystem.

[0050] As an example, as Figure 3 shown, the 12 excitation data of excitation points 1 to 4 (N = 4) in the X-axis direction, Y-axis direction, and Z-axis direction (M = 3) in the first subsystem A and the 12 response data of response points 5 to 8 in the X-axis direction, Y-axis direction, and Z-axis direction in the second subsystem B determine 12 * 12 third transfer function values. R = 12, and the third transfer function values are sorted in 12 rows and 12 columns to obtain the 12 * 12 target transfer function matrix of the coupled system , in this case, the target transfer function matrix of the coupling system is: = , where represents the third transfer function value determined by the excitation data with the excitation point numbered 1 and the excitation direction along the X-axis in the first subsystem A, and the response data with the response point numbered 5 and the response direction along the X-axis in the second subsystem B. represents the third transfer function value determined by the excitation data with the excitation point numbered 1 and the excitation direction along the X-axis in the first subsystem A, and the response data with the response point numbered 8 and the response direction along the Z-axis in the second subsystem B. represents the third transfer function value determined by the excitation data with the excitation point numbered 4 and the excitation direction along the Z-axis in the first subsystem A, and the response data with the response point numbered 5 and the response direction along the X-axis in the second subsystem B. represents the third transfer function value determined by the excitation data with the excitation point numbered 4 and the excitation direction along the Z-axis in the first subsystem A, and the response data with the response point numbered 8 and the response direction along the Z-axis in the first subsystem A.

[0051] As another example, as Figure 3 shown, 12 excitation data of the excitation points 5 to 8 (N = 4) in the second subsystem B along the X-axis, Y-axis, and Z-axis directions (M = 3) and 12 response data of the response points 1 to 4 in the first subsystem A along the X-axis, Y-axis, and Z-axis directions determine 12 * 12 third transfer function values. R = 12, and the third transfer function values are sorted in 12 rows and 12 columns to obtain the target transfer function matrix of the 12 * 12 coupling system , in this case, the target transfer function matrix of the coupling system is: = , where represents the third transfer function value determined by the excitation data with the excitation point numbered 5 and the excitation direction along the X-axis in the second subsystem B, and the response data with the response point numbered 1 and the response direction along the X-axis in the first subsystem A. represents the third transfer function value determined by the excitation data with the excitation point numbered 5 and the excitation direction along the X-axis in the second subsystem B, and the response data with the response point numbered 4 and the response direction along the Z-axis in the first subsystem A. represents the third transfer function value determined by the excitation data with the excitation point numbered 8 and the excitation direction along the Z-axis in the second subsystem B, and the response data with the response point numbered 1 and the response direction along the X-axis in the first subsystem A. It represents the third transfer function value determined by the excitation data with the excitation point numbered 8 and the excitation direction along the Z-axis in the second subsystem B, and the response data with the response point numbered 4 and the response direction along the Z-axis in the first subsystem A.

[0052] As an example, in step S302, the test device determines the excitation reduction matrix corresponding to each subsystem according to the excitation point coordinates corresponding to the excitation points in each subsystem; determines the response reduction matrix corresponding to each subsystem according to the response point coordinates corresponding to the response points in each subsystem; uses the excitation reduction matrix corresponding to each subsystem and the response reduction matrix corresponding to each subsystem to reduce the R*R target transfer function matrix corresponding to each subsystem, and obtains the L*L origin transfer function matrix corresponding to each subsystem in the two subsystems. The R*R target transfer function matrix corresponding to the coupled system is reduced by using the excitation reduction matrix corresponding to one subsystem and the response reduction matrix corresponding to the other subsystem to obtain the L*L origin transfer function matrix corresponding to the coupled system.

[0053] As Figure 3 shown, the test device determines the excitation reduction matrix corresponding to the first subsystem A according to the excitation point coordinates corresponding to the excitation points in the first subsystem A; determines the excitation reduction matrix corresponding to the second subsystem B according to the excitation point coordinates corresponding to the excitation points in the second subsystem B; determines the response reduction matrix corresponding to the first subsystem A according to the response point coordinates corresponding to the response points in the first subsystem A; determines the response reduction matrix corresponding to the second subsystem B according to the response point coordinates corresponding to the response points in the second subsystem B.

[0054] The test device uses the excitation reduction matrix corresponding to the first subsystem A and the response reduction matrix corresponding to the first subsystem A to reduce the target transfer function matrix corresponding to the first subsystem A to obtain the origin transfer function matrix corresponding to the first subsystem A.

[0055] The test device uses the excitation reduction matrix corresponding to the second subsystem B and the response reduction matrix corresponding to the second subsystem B to reduce the target transfer function matrix corresponding to the second subsystem B to obtain the origin transfer function matrix corresponding to the second subsystem B.

[0056] The test device uses the excitation reduction matrix corresponding to the first subsystem A and the excitation reduction matrix corresponding to the second subsystem B to reduce the target transfer function matrix corresponding to the coupled system to obtain the origin transfer function matrix Alternatively, the test device uses the excitation reduction matrix corresponding to the second subsystem B and the response reduction matrix corresponding to the first subsystem A to reduce the target transfer function matrix corresponding to the coupled system to obtain the origin transfer function matrix corresponding to the coupled system . Among them, = , , , and are all 6*6 matrices.

[0057] In this embodiment, the target transfer function matrix corresponding to each subsystem is used to reflect the transfer status of each subsystem, and the target transfer function matrix corresponding to the coupled system is used to reflect the transfer status between the two subsystems. To obtain the dynamic stiffness of the bushing, it is necessary to reduce the target transfer function matrices corresponding to the two subsystems and the target transfer function matrix corresponding to the coupled system to the origin between the two subsystems and the bushing, and obtain the origin transfer function matrices corresponding to the two subsystems and the origin transfer function matrix corresponding to the coupled system, so as to accurately and stably determine the dynamic stiffness of the bushing connected to the two subsystems based on the origin transfer function matrices corresponding to the two subsystems and the origin transfer function matrix corresponding to the coupled system.

[0058] In one embodiment, step S302, that is, based on the excitation point coordinates corresponding to the excitation data and the response point coordinates corresponding to the response data, reducing the target transfer function matrices corresponding to the two subsystems and the target transfer function matrix corresponding to the coupled system, and respectively determining the origin transfer function matrices corresponding to the two subsystems and the origin transfer function matrix corresponding to the coupled system, includes: S401: Based on the excitation point coordinates corresponding to the excitation data of each subsystem and the origin coordinates corresponding to each subsystem, determine the excitation reduction matrix of L* corresponding to each subsystem; S402: Based on the response point coordinates corresponding to the response data of each subsystem and the origin coordinates corresponding to each subsystem, determine *L response reduction matrix corresponding to each subsystem; S403: Use the excitation reduction matrices corresponding to the two subsystems and the response reduction matrices corresponding to the two subsystems to reduce the target transfer function matrix of the target system, and determine the origin transfer function matrix of the target system. The target system is a subsystem or a coupled system.

[0059] As an example, in step S401, the subsystem is the first subsystem or the second subsystem. Determining the excitation reduction matrix of L* corresponding to each subsystem includes the following steps: When the subsystem is the first subsystem, the origin coordinates corresponding to the first subsystem are (0, 0, 0). The test device determines the excitation reduction matrix corresponding to the first subsystem based on the excitation point coordinates corresponding to the excitation data of the first subsystem and the origin coordinates of the first subsystem. , = where, = where i is the number of excitation points, represents the excitation data of the i-th excitation point in the first subsystem, represents the equivalent excitation obtained by equivalenting the excitation data of the i-th excitation point in the first subsystem to the origin corresponding to the first subsystem, is the X-axis coordinate of the i-th excitation point in the first subsystem, is the Y-axis coordinate of the i-th excitation point in the first subsystem, is the Z-axis coordinate of the i-th excitation point in the first subsystem. is the reduction matrix for equivalenting the excitation data of the i-th excitation point in the first subsystem in M directions to the equivalent excitation corresponding to the origin, where, is a matrix with L rows and K columns. As Figure 3 shown, in the first subsystem A, N = 4, M = 3, i = 1, 2, 3, and 4, L = 6, and K = 24. The excitation reduction matrix corresponding to the first subsystem is = where, is a 6-row and 12-column excitation reduction matrix.

[0060] When the subsystem is the second subsystem, the origin coordinates corresponding to the second subsystem are (0, 0, 0). The test device determines the excitation reduction matrix corresponding to the second subsystem based on the excitation point coordinates corresponding to the excitation data of the second subsystem and the origin coordinates of the second subsystem. is: = where, = where i is the number of excitation points, represents the excitation data of the i-th excitation point in the second subsystem, represents the equivalent excitation obtained by equivalenting the excitation data of the i-th excitation point in the second subsystem to the origin corresponding to the second subsystem, is the X-axis coordinate of the i-th excitation point in the second subsystem, is the Y-axis coordinate of the i-th excitation point in the second subsystem, is the Z-axis coordinate of the i-th excitation point in the second subsystem. The reduction matrix for equivalent excitation corresponding to the origin of the excitation data of the $i$-th excitation point in the second subsystem in $M$ directions, where is an excitation reduction matrix with $L$ rows and $K$ columns. As Figure 3 shown, in the second subsystem $B$, $N = 4$, $M = 3$, $i = 5, 6, 7$ and $i = 8$, $L = 6$ and $K = 24$, the excitation reduction matrix corresponding to the second subsystem is = , where is an excitation reduction matrix with 6 rows and 12 columns.

[0061] As an example, in step S402, the subsystem is the first subsystem or the second subsystem, and the steps for determining the *$L$ response reduction matrix corresponding to each subsystem are as follows: When the subsystem is the first subsystem, the origin coordinates corresponding to the first subsystem are $(0, 0, 0)$. The test equipment determines the response reduction matrix corresponding to the first subsystem based on the response point coordinates corresponding to the response data of the first subsystem and the origin coordinates corresponding to the first subsystem as: = , = , where $j$ is the number of response points, represents the response data of the $j$-th response point of the first subsystem, represents the conversion of the response data of the origin corresponding to the first subsystem to the response data of the $j$-th response point in the first subsystem, is the $X$-axis coordinate of the $j$-th response point in the first subsystem, is the $Y$-axis coordinate of the $j$-th response point in the first subsystem, is the $Z$-axis coordinate of the $j$-th response point in the first subsystem. The reduction matrix for equivalent response corresponding to the origin of the response data of the $j$-th response point in the first subsystem in $M$ directions, where is rows and $L$ columns of the response reduction matrix. As Figure 3 shown, in the first subsystem $A$, $N = 4$, $M = 3$, $i = 1, 2, 3$ and $i = 4$, $L = 6$ and $K = 24$, the response reduction matrix corresponding to the first subsystem is = , where is a response reduction matrix with 12 rows and 6 columns.

[0062] When the subsystem is the second subsystem, the origin coordinates corresponding to the second subsystem are (0, 0, 0). The test device determines the response reduction matrix corresponding to the second subsystem based on the response point coordinates corresponding to the response data of the second subsystem and the origin coordinates corresponding to the second subsystem. It is: = , = , where j is the number of response points. represents the response data of the j-th response point in the second subsystem. represents converting the response data of the origin corresponding to the second subsystem to the response data of the j-th response point in the second subsystem. represents the equivalent response of the origin corresponding to the second subsystem. is the X-axis coordinate of the j-th response point in the second subsystem. is the Y-axis coordinate of the j-th response point in the second subsystem. is the Z-axis coordinate of the j-th response point in the second subsystem. is the reduction matrix for equivalent the response data of the j-th response point in the second subsystem in M directions to the equivalent response corresponding to the origin. Among them, is a response reduction matrix with L columns and K rows. As Figure 3 shown, in the second subsystem B, N = 4, M = 3, i = 5, 6, 7, and 8, L = 6, and K = 24. The response reduction matrix corresponding to the second subsystem is = , where is a 12-row and 6-column response reduction matrix.

[0063] As an example, in step S403, the test system reduces the target transfer function matrix corresponding to the subsystem based on the excitation reduction matrix and response reduction matrix corresponding to the same subsystem to determine the origin transfer function matrix corresponding to the subsystem; and reduces the target transfer function matrix corresponding to the coupled system based on the excitation reduction matrix corresponding to any one subsystem and the response reduction matrix corresponding to another subsystem to determine the origin transfer function matrix corresponding to the coupled system. For example, the test device uses the excitation reduction matrix corresponding to the first subsystem and the response reduction matrix corresponding to the first subsystem to reduce the target transfer function matrix corresponding to the first subsystem to determine the origin transfer function matrix corresponding to the first subsystem. The test device uses the excitation reduction matrix corresponding to the second subsystem and the response reduction matrix corresponding to the second subsystem to reduce the target transfer function matrix corresponding to the second subsystem to determine the origin transfer function matrix corresponding to the second subsystem. The test device uses the excitation reduction matrix of the first subsystem and the response reduction matrix of the second subsystem to reduce the target transfer function matrix of the coupled system to determine the origin transfer function matrix corresponding to the coupled system; or uses the excitation reduction matrix of the second subsystem and the response reduction matrix of the first subsystem to reduce the target transfer function matrix of the coupled system to determine the origin transfer function matrix corresponding to the coupled system.

[0064] In this embodiment, according to the excitation reduction matrix corresponding to the target system and the response reduction matrix corresponding to the target system, the target transfer function matrix of the target system is reduced to determine the origin transfer function matrix corresponding to the target system, the excitation points and response points in the target system are reduced to the origin corresponding to the target system, and the target initial transfer function matrix is reduced to the origin transfer function matrix corresponding to the origin, so as to make it feasible to determine the dynamic stiffness of the bushing according to the origin transfer function matrix.

[0065] In one embodiment, the origin transfer function matrix corresponding to the target system is the product of the inverse matrix of the response reduction matrix corresponding to the target system, the target transfer function matrix, and the inverse matrix of the excitation reduction matrix.

[0066] In this example, when the target system is the first subsystem, the origin transfer function matrix corresponding to the first subsystem is the product of the inverse matrix of the response reduction matrix corresponding to the first subsystem, the target transfer function matrix corresponding to the first subsystem, and the inverse matrix of the excitation reduction matrix corresponding to the first subsystem.

[0067] Among them, the origin transfer function matrix corresponding to the first subsystem refers to the matrix of L*L origin transfer function values between the equivalent excitations in L directions and the equivalent responses in L directions of the origin corresponding to the first subsystem after reducing the excitation points and response points of the first subsystem to the origin corresponding to the first subsystem. The equivalent excitation in the first subsystem refers to the excitation data of the excitation points of the first subsystem being equivalent to the excitation at the origin. The equivalent response in the first subsystem refers to the response data of the response points of the first subsystem being equivalent to the response at the origin. The origin transfer function value corresponding to the first subsystem refers to the transfer function value between the equivalent excitation and the equivalent response of the origin corresponding to the first subsystem.

[0068] The origin transfer function matrix corresponding to the first subsystem is: ( ), where is the excitation reduction matrix corresponding to the first subsystem, is the response reduction matrix corresponding to the first subsystem, represents the equivalent excitation of the origin corresponding to the first subsystem, represents the equivalent response of the origin corresponding to the first subsystem, () represents the generalized inverse matrix.

[0069] Understandably, as Figure 3 shown, the reduction relationship between the equivalent excitation of the origin corresponding to the first subsystem A and the excitation data of the 4 excitation points in the first subsystem A is: = , then = . Among them, = , respectively represent the equivalent excitations of the origin corresponding to the first subsystem A in 6 excitation directions of the X-axis direction, Y-axis direction, Z-axis direction, rotation direction along the X-axis, rotation direction along the Y-axis, and rotation direction along the Z-axis, that is, L = 6, is the excitation data of the excitation point numbered 1 in the first subsystem A, = , where is the excitation data of the excitation point numbered 1 in the X-axis direction in the first subsystem A, is the excitation data of the excitation point numbered 1 in the Y-axis direction in the first subsystem A, is the excitation data of the excitation point numbered 1 in the Z-axis direction in the first subsystem A. Similarly, = , = , = 。

[0070] The equivalent response of the origin corresponding to the first subsystem A and the response data of 4 response points in the first subsystem A The reduction relationship between them is: = That is = ( ) . Among them, = respectively represent the equivalent responses of the origin of the first subsystem A in 6 response directions: the X-axis direction, the Y-axis direction, the Z-axis direction, the rotation direction along the X-axis, the rotation direction along the Y-axis, and the rotation direction along the Z-axis. That is, L = 6, is the response data of the response point numbered 1 in the first subsystem A, = . Among them, is the response data of the response point numbered 1 in the first subsystem A in the X-axis direction, is the response data of the response point numbered 1 in the first subsystem A in the Y-axis direction, is the response data of the response point numbered 1 in the first subsystem A in the Z-axis direction. Similarly, = , = , = .

[0071] According to the equivalent response and equivalent excitation of the origin corresponding to the first subsystem A, determine the origin transfer function matrix corresponding to the first subsystem A is: = That is . Therefore, = ( ). From the target transfer function matrix of the first subsystem = it can be seen that = ( ).

[0072] In this example, when the target system is the second subsystem, the origin transfer function matrix corresponding to the second subsystem is the product of the inverse matrix of the response reduction matrix corresponding to the second subsystem, the target transfer function matrix corresponding to the second subsystem, and the inverse matrix of the excitation reduction matrix corresponding to the second subsystem.

[0073] Among them, the origin transfer function matrix corresponding to the second subsystem refers to the matrix of L*L origin transfer function values between the equivalent excitations in L directions and the equivalent responses in L directions of the origin corresponding to the second subsystem after reducing the excitation points and response points of the second subsystem to the origin corresponding to the second subsystem. The equivalent excitation in the second subsystem refers to the excitation data of the excitation point of the second subsystem being equivalent to the excitation at the origin. The equivalent response in the second subsystem refers to the response data of the response point of the second subsystem being equivalent to the response at the origin. The origin transfer function value corresponding to the second subsystem refers to the transfer function value between the equivalent excitation and the equivalent response of the origin corresponding to the second subsystem.

[0074] As an example, the origin transfer function matrix corresponding to the second subsystem is: ( ) ; where, is the excitation reduction matrix corresponding to the second subsystem, is the response reduction matrix corresponding to the second subsystem. represents the equivalent excitation of the origin corresponding to the second subsystem, represents the equivalent response of the origin corresponding to the second subsystem.

[0075] Understandably, as Figure 3 shown, the relationship between the equivalent excitation of the origin corresponding to the second subsystem B and the excitation data of the 4 excitation points in the second subsystem B is: = , that is = ( ) . Among them, = , respectively represent the equivalent excitations of the origin corresponding to the second subsystem B in 6 excitation directions of the X-axis direction, Y-axis direction, Z-axis direction, rotation direction along the X-axis, rotation direction along the Y-axis, and rotation direction along the Z-axis, that is, L = 6, is the excitation data of the excitation point numbered 5 in the second subsystem B, = , where, is the excitation data of the excitation point numbered 5 in the second subsystem B in the X-axis direction, is the excitation data of the excitation point numbered 5 in the second subsystem B in the Y-axis direction, is the excitation data of the excitation point numbered 5 in the second subsystem B in the Z-axis direction. Similarly, = , = , = 。

[0076] The equivalent response of the origin corresponding to the second subsystem B and the response data of 4 response points in the second subsystem B The reduction relationship between them is: = , that is = ( ) . Among them, = , respectively represent the equivalent responses of the origin corresponding to the second subsystem B in 6 response directions of the X-axis direction, Y-axis direction, Z-axis direction, rotation direction along the X-axis, rotation direction along the Y-axis, and rotation direction along the Z-axis, that is, L = 6, is the response data of the response point numbered 5 in the second subsystem B, = , among which, is the response data of the response point numbered 5 in the second subsystem B in the X-axis direction, is the response data of the response point numbered 5 in the second subsystem B in the Y-axis direction, is the response data of the response point numbered 5 in the second subsystem B in the Z-axis direction. Similarly, , = , = .

[0077] According to the equivalent response and equivalent excitation of the origin corresponding to the second subsystem B, determine the origin transfer function matrix corresponding to the second subsystem as: = , that is , therefore, = ( ), from the target transfer function matrix of the second subsystem = it can be seen that ( ).

[0078] In this example, when the target system is a coupled system, the origin transfer function matrix corresponding to the coupled system is the product of the inverse matrix of the response reduction matrix corresponding to the coupled system, the target transfer function matrix corresponding to the coupled system, and the inverse matrix of the excitation reduction matrix corresponding to the coupled system.

[0079] Among them, the origin transfer function matrix corresponding to the coupling system refers to the matrix of L*L origin transfer function values between the equivalent excitations of the origin of the first subsystem in L directions and the equivalent responses of the origin of the second subsystem in L directions after reducing the excitation point of the first subsystem to the origin corresponding to the first subsystem and reducing the response point of the second subsystem to the origin corresponding to the second subsystem, or the matrix of L*L origin transfer function values between the equivalent excitations of the origin of the first subsystem in L directions and the equivalent responses of the origin of the second subsystem in L directions.

[0080] As an example, the origin transfer function matrix corresponding to the coupling system is: ( ); Or, the origin transfer function matrix corresponding to the coupling system is: ( ), and = .

[0081] Understandably, according to the equivalent excitation of the origin corresponding to the first subsystem A and the equivalent response of the origin corresponding to the second subsystem B, the origin transfer function matrix corresponding to the coupling system is: = ∙ That is, ( ) = ∙ That is, = ( ) From the target transfer function matrix of the coupling system = it can be seen that ( ). Or, according to the equivalent excitation of the origin corresponding to the second subsystem B and the equivalent response of the origin corresponding to the first subsystem A, the origin transfer function matrix corresponding to the coupling system is: = ∙ That is, ( ) = ∙ That is, = ( ) , from the target transfer function matrix of the coupling system = it can be seen that .

[0082] In this embodiment, for the dynamic stiffness test of the bushing, the bushing is respectively connected to the first subsystem and the second subsystem. After testing multiple excitation points and multiple response points in the first subsystem and the second subsystem, it is necessary to equivalent multiple excitation points and multiple response points to the origin corresponding to the first subsystem and the origin corresponding to the second subsystem (i.e., the center point of the bushing connected to the first subsystem and the second subsystem). When obtaining the origin corresponding to the first subsystem as both the excitation point and the response point, the origin transfer function matrix corresponding to the first subsystem formed by the origin transfer function values in the L excitation directions and response directions, the origin transfer function matrix corresponding to the second subsystem formed by the origin transfer function values in the L excitation directions and response directions when the origin corresponding to the second subsystem is used as both the excitation point and the response point, and the origin transfer function matrix corresponding to the coupling system when the origin corresponding to the first subsystem is used as the excitation point and the origin corresponding to the second subsystem is used as the response point, or the origin transfer function matrix corresponding to the coupling system when the origin corresponding to the second subsystem is used as the excitation point and the origin corresponding to the first subsystem is used as the response point, are used to accurately and stably determine the dynamic stiffness of the bushing within a relatively large frequency range (especially within the high-frequency range) according to the origin transfer function matrix corresponding to the first subsystem, the origin transfer function matrix corresponding to the second subsystem, and the origin transfer function matrix corresponding to the coupling system.

[0083] In one embodiment, the excitation data is the excitation force and the response data is the acceleration.

[0084] As an example, as Figure 3 shown, the excitation data of the 4 excitation points in the first subsystem A in the X-axis direction, Y-axis direction, and Z-axis direction are all excitation forces, and the equivalent excitation of the origin corresponding to the first subsystem A includes the excitation forces of the origin corresponding to the first subsystem A in 6 excitation directions: the X-axis direction, Y-axis direction, Z-axis direction, rotation direction along the X-axis, rotation direction along the Y-axis, and rotation direction along the Z-axis. The excitation data of the 4 excitation points in the second subsystem B in the X-axis direction, Y-axis direction, and Z-axis direction are all excitation forces, and the equivalent excitation of the origin corresponding to the second subsystem B includes the excitation forces of the origin corresponding to the second subsystem B in 6 excitation directions: the X-axis direction, Y-axis direction, Z-axis direction, rotation direction along the X-axis, rotation direction along the Y-axis, and rotation direction along the Z-axis.

[0085] The response data of the 4 response points in the first subsystem A in the X-axis direction, Y-axis direction, and Z-axis direction are all accelerations, and the equivalent response of the origin corresponding to the first subsystem A Including the accelerations in six response directions of the origin corresponding to the first subsystem A in the X-axis direction, Y-axis direction, Z-axis direction, rotation direction along the X-axis, rotation direction along the Y-axis, and rotation direction along the Z-axis. The response data of the four response points in the second subsystem B in the X-axis direction, Y-axis direction, and Z-axis direction are all accelerations, and the equivalent response of the origin corresponding to the second subsystem B Including the accelerations in six response directions of the origin corresponding to the second subsystem B in the X-axis direction, Y-axis direction, Z-axis direction, rotation direction along the X-axis, rotation direction along the Y-axis, and rotation direction along the Z-axis.

[0086] In one embodiment, the dynamic stiffness of the bushing is the quotient of the transfer function matrix of the origin corresponding to the coupled system and the target determinant; The target determinant is the difference between the first matrix product and the second matrix product. The first matrix product is the product of the transfer function matrices of the origins corresponding to the two subsystems, and the second matrix product is the product of the transfer function matrix of the origin corresponding to the coupled system and the transfer function matrix of the origin corresponding to the coupled system.

[0087] As an example, the dynamic stiffness of the bushing . Among them, the target determinant is a 2×2 determinant: . In this example, through the coupled system formed by the two subsystems and the bushing connected to the coupled system, the dynamic stiffness test of the bushing is carried out. There is no strict requirement for the test frequency, and the dynamic stiffness test of the bushing can be carried out within a relatively large frequency range. Taking Figure 3 the first subsystem A, the second subsystem B in and the bushing C respectively connected to the first subsystem A and the second subsystem B as an example, the equivalent excitation of the origin corresponding to the first subsystem A is , The relationship with the deformation of the bushing relative to the first subsystem A and the second subsystem B is: = ( ), The relationship with the deformation of the bushing relative to the first subsystem A and the second subsystem B is: = ( ), where is the dynamic stiffness of the first subsystem A, is the dynamic stiffness of the second subsystem B, is the deformation of the bushing relative to the first subsystem A, is the deformation of the bushing relative to the second subsystem B. The above relational expressions can also be written as the following relational expression (1): In one example, there is the following relationship between the equivalent excitations of the first subsystem A and the second subsystem B and the deformation of the bushing: That is, there is the following relational expression (2): The elements containing the dynamic stiffness in the above relational expression (1) are corresponding to the ( ) at the same position in the relational expression (2), and we get = , that is, the dynamic stiffness of the bushing . Among them, is the target determinant, is the first matrix product, is the second matrix product, is the origin transfer function matrix corresponding to the first subsystem, is the origin transfer function matrix corresponding to the second subsystem, is the origin transfer function matrix corresponding to the coupled system.

[0088] In another example, there is the following relationship between the equivalent excitations of the first subsystem and the second subsystem and the deformation of the bushing: That is, there is the following relational expression (3): The elements containing the dynamic stiffness in the above relational expression (1) are corresponding to the ( ) at the same position in the relational expression (3), and we get = , that is, the dynamic stiffness of the bushing . Among them, is the target determinant, is the first matrix product, is the second matrix product, is the origin transfer function matrix corresponding to the first subsystem, is the origin transfer function matrix corresponding to the second subsystem, is the origin transfer function matrix corresponding to the coupled system.

[0089] In this embodiment, according to the relationship between the deformation of the bushing relative to each subsystem in the two subsystems, and the relationship between the equivalent excitation of each subsystem in the two subsystems and the deformation of the bushing, the origin transfer function matrices corresponding to the two subsystems and the origin transfer function matrix corresponding to the coupled system are processed to accurately determine the dynamic stiffness of the bushing connected to the two subsystems, realizing the test of the dynamic stiffness of the bushing. The above dynamic stiffness test method is realized through the coupled system formed by the two subsystems and the bushing connected to the coupled system. This method has no strict limit on the frequency and can test the dynamic stiffness of the bushing within a large frequency range, with strong applicability.

[0090] It should be understood that the magnitudes of the sequence numbers of the steps in the above embodiments do not mean the order of execution. The execution order of each process should be determined according to its function and internal logic, and should not constitute any limitation to the implementation process of the embodiments of the present invention.

[0091] In one embodiment, a bushing dynamic stiffness test device is provided, which corresponds one-to-one with the bushing dynamic stiffness test method in the above embodiment. As Figure 2 shown, the bushing dynamic stiffness test device includes an initial transfer function matrix determination module 21, an origin transfer function matrix determination module 22, and a dynamic stiffness determination module 23. The detailed description of each functional module is as follows: The initial transfer function matrix determination module 21 is used to excite the two subsystems to determine the initial transfer function matrix; the initial transfer function matrix is a matrix formed by K*K initial transfer function values determined based on K excitation data and K response data corresponding to each excitation data; The origin transfer function matrix determination module 22 is used to reduce the initial transfer function matrix to determine the origin transfer function matrices corresponding to the two subsystems and the origin transfer function matrix corresponding to the coupled system; the origin transfer function matrix corresponding to each subsystem is a matrix formed by L*L origin transfer function values determined based on L excitation data corresponding to the subsystem and L response data corresponding to each excitation data in the same subsystem; the origin transfer function matrix corresponding to the coupled system is a matrix formed by L*L origin transfer function values determined based on L excitation data corresponding to any one of the subsystems and L response data corresponding to each excitation data in the other subsystem, where L < K; The dynamic stiffness determination module 23 determines the dynamic stiffness of the bushing based on the origin transfer function matrices corresponding to the two subsystems and the origin transfer function matrix corresponding to the coupled system.

[0092] In one embodiment, the initial transfer function matrix determination module 21 includes: A data acquisition sub-module, which is used to stimulate two subsystems, determine K excitation data and K response data corresponding to each excitation data. The K excitation data are data formed by stimulating 2*N excitation points along M excitation directions respectively, and the K response data are data formed by the 2*N response points responding in the M response directions for each stimulation. K = 2*N*M, N≥2, M≥1; An initial transfer function matrix determination sub-module, which determines an initial transfer function value based on each excitation data and a response data corresponding to the excitation data, and determines an initial transfer function matrix based on K*K initial transfer function values.

[0093] In one embodiment, the origin transfer function matrix determination module 22 includes: A target transfer function matrix determination sub-module, which is used to reduce the initial transfer function matrix to determine the target transfer function matrix corresponding to the two subsystems and the target transfer function matrix corresponding to the coupled system; the target transfer function matrix corresponding to each subsystem is a matrix formed by R*R initial transfer function values determined based on R excitation data corresponding to the subsystem and R response data corresponding to each excitation data in the same subsystem; the origin transfer function matrix corresponding to the coupled system is a matrix formed by R*R initial transfer function values determined based on R excitation data corresponding to any one subsystem and R response data corresponding to each excitation data in the other subsystem, L<R= ; An origin transfer function matrix determination sub-module, which reduces the target transfer function matrix corresponding to the two subsystems and the target transfer function matrix corresponding to the coupled system based on the excitation point coordinates corresponding to the excitation data and the response point coordinates corresponding to the response data, and respectively determines the origin transfer function matrix corresponding to the two subsystems and the origin transfer function matrix corresponding to the coupled system.

[0094] In one embodiment, the origin transfer function matrix determination sub-module includes: An excitation reduction matrix determination unit, which determines the excitation reduction matrix of L* corresponding to each subsystem based on the excitation point coordinates corresponding to the excitation data of each subsystem and the origin coordinates corresponding to each subsystem; A response reduction matrix determination unit, which determines the response reduction matrix of *L corresponding to each subsystem based on the response point coordinates corresponding to the response data of each subsystem and the origin coordinates corresponding to each subsystem; An origin transfer function matrix determination unit, which is used to reduce the target transfer function matrix of the target system by using the excitation reduction matrices corresponding to the two subsystems and the response reduction matrices corresponding to the two subsystems, and determine the origin transfer function matrix corresponding to the target system. The target system is a subsystem or a coupled system.

[0095] For the specific limitations of the bushing dynamic stiffness testing device, reference can be made to the limitations of the bushing dynamic stiffness testing method in the above text, which will not be elaborated here. Each module in the above bushing dynamic stiffness testing device can be implemented in whole or in part by software, hardware, and their combination. Each of the above modules can be embedded in the processor of the testing device in hardware form or be independent of it, or can be stored in the memory of the testing device in software form, so that the processor can call and execute the operations corresponding to each of the above modules.

[0096] In one embodiment, a testing device is provided, including a memory, a processor, and a computer program stored on the memory and executable on the processor. When the processor executes the computer program, the bushing dynamic stiffness testing method in the above embodiment is implemented. For example Figure 1 as shown in S101 - S103. To avoid repetition, it will not be elaborated here. Alternatively, when the processor executes the computer program, the functions of each module / unit in this embodiment of the bushing dynamic stiffness testing device are implemented. For example Figure 2 the functions of the initial transfer function matrix determination module 21, the origin transfer function matrix determination module 22, and the dynamic stiffness determination module 23 as shown. To avoid repetition, it will not be elaborated here.

[0097] In one embodiment, a bushing dynamic stiffness testing system is provided, including the above testing device and a coupling system connected to the testing device. The coupling system includes two subsystems, and the two subsystems are respectively connected to both ends of the bushing. The testing device is respectively connected to the two subsystems, and is used to excite the two subsystems in a relatively large frequency range, determine the initial transfer function matrix, reduce the initial transfer function matrix, determine the origin transfer function matrices corresponding to the two subsystems and the origin transfer function matrix corresponding to the coupling system, and determine the dynamic stiffness of the bushing based on the origin transfer function matrices corresponding to the two subsystems and the origin transfer function matrix corresponding to the coupling system, so as to realize the dynamic stiffness testing of the bushing in a relatively large frequency range.

[0098] In one embodiment, a computer-readable storage medium is provided. A computer program is stored on the computer-readable storage medium. When the computer program is executed by the processor, the bushing dynamic stiffness testing method in the above embodiment is implemented. For example Figure 1 as shown in S101 - S103. To avoid repetition, it will not be elaborated here. Alternatively, when the computer program is executed by the processor, the functions of each module / unit in this embodiment of the bushing dynamic stiffness testing device are implemented. For example Figure 2 the functions of the initial transfer function matrix determination module 21, the origin transfer function matrix determination module 22, and the dynamic stiffness determination module 23 as shown. To avoid repetition, it will not be elaborated here.

[0099] Those skilled in the art can clearly understand that, for the convenience and brevity of description, only the above-mentioned division of each functional unit and module is used as an example. In actual applications, the above-mentioned functions can be allocated to different functional units and modules according to needs, that is, the internal structure of the device is divided into different functional units or modules to complete all or part of the functions described above.

[0100] The above-described embodiments are only used to illustrate the technical solutions of the present invention, rather than to limit it; although the present invention has been described in detail with reference to the foregoing embodiments, those of ordinary skill in the art should understand that they can still modify the technical solutions recorded in the foregoing embodiments, or perform equivalent replacements on some of the technical features; and these modifications or replacements do not cause the essence of the corresponding technical solutions to deviate from the spirit and scope of the technical solutions of the embodiments of the present invention, and should all be included in the protection scope of the present invention.

Claims

1. A method for testing the dynamic stiffness of a bushing, characterized in that, Applicable to a coupling system, the coupling system includes two subsystems, and the two subsystems are respectively connected to both ends of a bushing. The dynamic stiffness test method of the bushing includes: Exciting the two subsystems to determine an initial transfer function matrix; the initial transfer function matrix is a matrix formed by K*K initial transfer function values determined based on K excitation data and K response data corresponding to each excitation data. Reducing the initial transfer function matrix to determine the origin transfer function matrix corresponding to the two subsystems and the origin transfer function matrix corresponding to the coupling system; the origin transfer function matrix corresponding to each subsystem is a matrix formed by L*L origin transfer function values determined based on L excitation data corresponding to the subsystem and L response data corresponding to each excitation data in the same subsystem; the origin transfer function matrix corresponding to the coupling system is a matrix formed by L*L origin transfer function values determined based on L excitation data corresponding to any one of the subsystems and L response data corresponding to each excitation data in the other subsystem, where 2 < L < K. Determining the dynamic stiffness of the bushing based on the origin transfer function matrix corresponding to the two subsystems and the origin transfer function matrix corresponding to the coupling system.

2. The bushing dynamic stiffness test method according to claim 1, wherein Each subsystem includes N indicating points, and the indicating points are excitation points and response points. The exciting the two subsystems to determine the initial transfer function matrix includes: Exciting the two subsystems to determine K excitation data and K response data corresponding to each excitation data. The K excitation data are data formed by exciting 2*N excitation points along M excitation directions respectively, and the K response data are data formed by the 2*N response points responding in M response directions for each excitation, where K = 2*N*M, N ≥ 2, and M ≥ 1. Determining an initial transfer function value based on an excitation data and a response data corresponding to the excitation data, and determining the initial transfer function matrix based on K*K of the initial transfer function values.

3. The bushing dynamic stiffness testing method according to claim 1, characterized in that The reducing the initial transfer function matrix to determine the origin transfer function matrix corresponding to the two subsystems and the origin transfer function matrix corresponding to the coupling system includes: Reduce the initial transfer function matrix to determine the target transfer function matrix corresponding to the two subsystems and the target transfer function matrix corresponding to the coupling system; the target transfer function matrix corresponding to each of the subsystems is a matrix formed by R×R initial transfer function values determined based on the R excitation data corresponding to the subsystem and the R response data corresponding to each of the excitation data in the same subsystem; the origin transfer function matrix corresponding to the coupling system is a matrix formed by R×R initial transfer function values determined based on the R excitation data corresponding to any one of the subsystems and the R response data corresponding to each of the excitation data in the other subsystem, where L < R = ; Reducing the target transfer function matrix corresponding to the two subsystems and the target transfer function matrix corresponding to the coupling system based on the excitation point coordinates corresponding to the excitation data and the response point coordinates corresponding to the response data, and respectively determining the origin transfer function matrix corresponding to the two subsystems and the origin transfer function matrix corresponding to the coupling system.

4. The bushing dynamic stiffness testing method according to claim 3, characterized in that, The reducing the target transfer function matrix corresponding to the two subsystems and the target transfer function matrix corresponding to the coupling system based on the excitation point coordinates corresponding to the excitation data and the response point coordinates corresponding to the response data, and respectively determining the origin transfer function matrix corresponding to the two subsystems and the origin transfer function matrix corresponding to the coupling system includes: Determine the excitation reduction matrix corresponding to each subsystem based on the excitation point coordinates corresponding to the excitation data of each subsystem and the origin coordinates corresponding to each subsystem ; Based on the response point coordinates corresponding to the response data of each subsystem and the origin coordinates corresponding to each subsystem, determine the response reduction matrix of *L for each subsystem; Using the excitation reduction matrix corresponding to the two subsystems and the response reduction matrix corresponding to the two subsystems to reduce the target transfer function matrix of the target system to determine the origin transfer function matrix corresponding to the target system. The target system is two subsystems or the coupled system.

5. The bushing dynamic stiffness testing method according to claim 4, characterized in that, The origin transfer function matrix corresponding to the target system is the product of the inverse matrix of the response reduction matrix corresponding to the target system, the target transfer function matrix, and the inverse matrix of the excitation reduction matrix.

6. The bushing dynamic stiffness testing method according to claim 1, characterized in that The dynamic stiffness of the bushing is the quotient of the origin transfer function matrix corresponding to the coupled system and the target determinant; The target determinant is the difference between a first matrix product and a second matrix product. The first matrix product is the product of the origin transfer function matrices corresponding to the two subsystems, and the second matrix product is the product of the origin transfer function matrix corresponding to the coupled system and the origin transfer function matrix corresponding to the coupled system.

7. A bushing dynamic stiffness test device, characterized in that, Including: An initial transfer function matrix determination module, configured to excite two subsystems to determine an initial transfer function matrix; The initial transfer function matrix is a matrix formed by K*K initial transfer function values determined based on K excitation data and K response data corresponding to each of the excitation data. An origin transfer function matrix determination module, configured to reduce the initial transfer function matrix to determine the origin transfer function matrix corresponding to the two subsystems and the origin transfer function matrix corresponding to the coupled system; the origin transfer function matrix corresponding to each subsystem is a matrix formed by L*L origin transfer function values determined based on L excitation data corresponding to the subsystem and L response data corresponding to each of the excitation data in the same subsystem; the origin transfer function matrix corresponding to the coupled system is a matrix formed by L*L origin transfer function values determined based on L excitation data corresponding to any one of the subsystems and L response data corresponding to each of the excitation data in the other subsystem, where L < K; A dynamic stiffness determination module, based on the origin transfer function matrices corresponding to the two subsystems and the origin transfer function matrix corresponding to the coupled system, determines the dynamic stiffness of the bushing.

8. A testing device, comprising a memory, a processor, and a computer program stored in the memory and executable on the processor, characterized in that, When the processor executes the computer program, it implements the bushing dynamic stiffness test method according to any one of claims 1 to 6.

9. A bushing dynamic stiffness test system, characterized in that Including the test device according to claim 8 and a coupled system connected to the test device. The coupled system includes two subsystems, and the two subsystems are respectively connected to both ends of the bushing.

10. A computer-readable storage medium storing a computer program, characterized in that, When the computer program is executed by a processor, it implements the bushing dynamic stiffness test method according to any one of claims 1 to 6.

Images