Bush rigidity-based load secondary decomposition general modeling method, device and equipment

Abstract

The invention relates to a load secondary decomposition general modeling method, device and equipment based on bushing rigidity. Comprising the steps of obtaining a to-be-simulated geometric model and performing geometric processing to obtain a processed model; based on the whole vehicle coordinate system, the sizes of the rubber bushing in three directions and the clearance amount between the rubber bushing and the target structural member are obtained, and the load critical value in each direction is determined by combining the nonlinear stiffness curve. And determining a first load applied to the hard point through the first RBE3 unit and a second load applied to the contact position through the second RBE3 unit according to the total load in the three directions and the critical value. And dividing grids for the processed model, endowing material attributes, constructing a new finite element model in combination with boundary constraints and loads, and outputting a simulation result. Therefore, the problems that a current structural part contains a rubber bushing, but the nonlinear rigidity of the rubber bushing is not considered, and the stress position of the structural part is not correctly determined are solved, and the simulation precision and accuracy of the structural part are greatly improved.

Description

TECHNICAL FIELD

[0001] The present application relates to the technical field of vehicles, in particular to a load secondary decomposition general modeling method based on bushing stiffness, device and equipment. BACKGROUND

[0002] In the field of strength calculation of chassis suspension bracket type structural members, the accuracy of load decomposition and finite element modeling directly affects the reliability of structural design.

[0003] In related technologies, the load decomposition of chassis suspension bracket type structural members is often based on the coordinates of the structural hard points. In the calculation of strength, the conventional finite element modeling directly applies loads to the hard point positions of the structure, and connects all the contact surfaces of the structural members through RBE2 / RBE3 units, including the contact surfaces that may be in contact with the structural members due to the stress deformation of the rubber bushing. The load is directly applied to the hard point positions of the structure.

[0004] However, in this modeling method, all the contact surfaces are under stress, and the stiffness of the rubber is assumed to be linear for modeling, without reflecting the physical properties of the hyperelasticity of the rubber bushing. Its nonlinear stiffness is not reflected in the finite element model, which does not conform to the actual situation. SUMMARY

[0005] The present application provides a load secondary decomposition general modeling method based on bushing stiffness to solve the problem that the nonlinear stiffness of the rubber bushing is not considered and the stress position of the structural member is not correctly determined in the current structural member containing the rubber bushing, greatly improving the simulation accuracy and accuracy of the structural member.

[0006] To achieve the above-mentioned purpose, the first aspect of the present application proposes a load secondary decomposition general modeling method based on bushing stiffness, comprising the following steps: Obtain a simulation geometry model, and perform geometric processing on the simulation geometry model to obtain a processed geometry model; Based on the vehicle coordinate system, obtain the size of the first to third directions of the rubber bushing in the processed geometry model and the gap amount of the first to third directions between the rubber bushing and the target structural member in the processed geometry model, and based on the size of the first to third directions and the gap amount of the first to third directions, determine the load critical value of the first to third directions according to the nonlinear stiffness curve of the rubber bushing; Based on the total load to be calculated in the first to third directions, determine the first load in the first to third directions applied to the hard point position through the first RBE3 unit and the second load in the first to third directions applied to the contact position of the rubber bushing and the target structural member through the second RBE3 unit according to the load critical value of the first to third directions; The processed geometric model is meshed and material properties are assigned, and based on preset boundary constraint conditions, a new finite element model is obtained by processing the meshed geometric model according to the first load in the first to third directions and the second load in the first to third directions, and a simulation result is output based on the new finite element model.

[0007] According to one embodiment of the application, the first to third direction load critical values are determined according to the nonlinear stiffness curve of the rubber bushing based on the first to third direction sizes and the first to third direction gap amounts, comprising: According to the size of each direction and the gap amount of the corresponding direction, the actual displacement parameters of the rubber bushing in the first to third directions are obtained; The displacement parameter set of the rubber bushing in the first to third directions is determined according to the nonlinear stiffness curve of the rubber bushing; The first to third direction load critical values are determined according to the actual displacement parameters in the first to third directions and the displacement parameter set in the first to third directions.

[0008] According to one embodiment of the application, the first to third direction load critical values are determined according to the actual displacement parameters in the first to third directions and the displacement parameter set in the first to third directions, comprising: It is judged whether the actual displacement parameter in the first direction is in the displacement parameter set in the first direction; If the actual displacement parameter in the first direction is in the displacement parameter set in the first direction, the first direction load critical value is obtained based on the displacement parameter set in the first direction, otherwise, the first direction load critical value is obtained by interpolation according to the displacement parameter set in the first direction based on the actual displacement parameter in the first direction; It is judged whether the actual displacement parameter in the second direction is in the displacement parameter set in the second direction; If the actual displacement parameter in the second direction is in the displacement parameter set in the second direction, the second direction load critical value is obtained based on the displacement parameter set in the second direction, otherwise, the second direction load critical value is obtained by interpolation according to the displacement parameter set in the second direction based on the actual displacement parameter in the second direction; It is judged whether the actual displacement parameter in the third direction is in the displacement parameter set in the third direction; If the actual displacement parameter in the third direction is in the displacement parameter set in the third direction, the third direction load critical value is obtained based on the displacement parameter set in the third direction, otherwise, the third direction load critical value is obtained by interpolation according to the displacement parameter set in the third direction based on the actual displacement parameter in the third direction.

[0009] According to one embodiment of the present application, the determining the first load in the first to third directions applied to the hard point position by the first RBE3 unit and the second load in the first to third directions applied to the contact position of the rubber bushing and the target structural member by the second RBE3 unit based on the total load in the first to third directions to be calculated comprises: determining whether the total load in the first direction is greater than the load critical value in the first direction; if the total load in the first direction is greater than the load critical value in the first direction, calculating a first difference between the total load in the first direction and the load critical value in the first direction, taking the load critical value in the first direction as the first load in the first direction applied to the hard point position by the first RBE3 unit, and taking the first difference as the second load in the first direction applied to the contact position by the second RBE3 unit; otherwise, taking the total load in the first direction as the first load in the first direction applied to the hard point position by the first RBE3 unit.

[0010] According to one embodiment of the present application, the determining the first load in the first to third directions applied to the hard point position by the first RBE3 unit and the second load in the first to third directions applied to the contact position of the rubber bushing and the target structural member by the second RBE3 unit based on the total load in the first to third directions to be calculated comprises: determining whether the total load in the second direction is greater than the load critical value in the second direction; if the total load in the second direction is greater than the load critical value in the second direction, calculating a second difference between the total load in the second direction and the load critical value in the second direction, taking the load critical value in the second direction as the first load in the second direction applied to the hard point position by the first RBE3 unit, and taking the second difference as the second load in the second direction applied to the contact position by the second RBE3 unit; otherwise, taking the total load in the second direction as the first load in the second direction applied to the hard point position by the first RBE3 unit.

[0011] According to one embodiment of the present application, the determining the first load in the first to third directions applied to the hard point position by the first RBE3 unit and the second load in the first to third directions applied to the contact position of the rubber bushing and the target structural member by the second RBE3 unit based on the total load in the first to third directions to be calculated comprises: determining whether the total load in the third direction is greater than the load critical value in the third direction; If the total load in the third direction is greater than the load threshold in the third direction, a third difference between the total load in the third direction and the load threshold in the third direction is calculated, the load threshold in the third direction is taken as the first load in the third direction applied to the hard point position by the first RBE3 unit, and the third difference is taken as the second load in the third direction applied to the contact position by the second RBE3 unit.

[0012] Otherwise, the total load in the third direction is taken as the first load in the third direction applied to the hard point position by the first RBE3 unit.

[0013] According to the universal modeling method for load secondary decomposition based on bush stiffness provided in the embodiments of the present application, the size and gap amount of the rubber bush are obtained by processing the to-be-simulated geometric model and combining the vehicle coordinate system, the load threshold is determined according to the nonlinear stiffness curve of the rubber bush, and then the loads at the hard point and the contact position are decomposed from the total load. After grid division, material assignment and constraint setting, the finite element model is constructed and the simulation result is output. In this way, the problem that the nonlinear stiffness of the rubber bush is not considered and the force position of the structural part is not correctly determined in the current structural part containing the rubber bush is solved, and the simulation accuracy and accuracy of the structural part are greatly improved.

[0014] To achieve the above object, the second aspect of the embodiments of the present application provides a universal modeling device for load secondary decomposition based on bush stiffness, comprising: A first obtaining module is configured to obtain a to-be-simulated geometric model and perform geometric processing on the to-be-simulated geometric model to obtain a processed geometric model. A second obtaining module is configured to obtain, based on a vehicle coordinate system, the size of a rubber bush in a first direction to a third direction and the gap amount between the rubber bush and a target structural part in the first direction to the third direction in the processed geometric model, and determine the load threshold in the first direction to the third direction according to the nonlinear stiffness curve of the rubber bush based on the size in the first direction to the third direction and the gap amount in the first direction to the third direction. A determining module is configured to determine, based on the total load in the first direction to the third direction to be calculated, the first load in the first direction to the third direction applied to the hard point position by the first RBE3 unit and the second load in the first direction to the third direction applied to the contact position between the rubber bush and the target structural part by the second RBE3 unit according to the load threshold in the first direction to the third direction. The decomposition module divides the processed geometric model into grids and assigns material properties, and based on preset boundary constraint conditions, processes the geometric model after grid division according to the first load in the first to third directions and the second load in the first to third directions to obtain a new finite element model, and outputs a simulation result based on the new finite element model. According to an embodiment of the present application, According to an embodiment of the present application, the second acquisition module is specifically configured to: According to the size of each direction and the gap amount of the corresponding direction, the actual displacement parameters of the rubber bushing in the first to third directions are obtained; According to the nonlinear stiffness curve of the rubber bushing, the displacement parameter set of the rubber bushing in the first to third directions is determined; According to the actual displacement parameters in the first to third directions and the displacement parameter set in the first to third directions, the load critical values in the first to third directions are determined.

[0015] According to an embodiment of the present application, the second acquisition module is further configured to: determine whether the actual displacement parameter in the first direction is in the displacement parameter set in the first direction; If the actual displacement parameter in the first direction is in the displacement parameter set in the first direction, the load critical value in the first direction is obtained based on the displacement parameter set in the first direction, otherwise, the load critical value in the first direction is obtained by interpolation according to the displacement parameter set in the first direction based on the actual displacement parameter in the first direction; determine whether the actual displacement parameter in the second direction is in the displacement parameter set in the second direction; If the actual displacement parameter in the second direction is in the displacement parameter set in the second direction, the load critical value in the second direction is obtained based on the displacement parameter set in the second direction, otherwise, the load critical value in the second direction is obtained by interpolation according to the displacement parameter set in the second direction based on the actual displacement parameter in the second direction; determine whether the actual displacement parameter in the third direction is in the displacement parameter set in the third direction; If the actual displacement parameter in the third direction is in the displacement parameter set in the third direction, the load critical value in the third direction is obtained based on the displacement parameter set in the third direction, otherwise, the load critical value in the third direction is obtained by interpolation according to the displacement parameter set in the third direction based on the actual displacement parameter in the third direction.

[0016] According to an embodiment of the present application, the determination module is specifically configured to: determine whether the total load in the first direction is greater than the load critical value in the first direction; if the total load in the first direction is greater than the load threshold in the first direction, calculating a first difference between the total load in the first direction and the load threshold in the first direction, and taking the load threshold in the first direction as a first load in the first direction applied to the hard point position by the first RBE3 unit and the first difference as a second load in the first direction applied to the contact position by the second RBE3 unit; otherwise, taking the total load in the first direction as the first load in the first direction applied to the hard point position by the first RBE3 unit.

[0017] According to an embodiment of the present application, the determining module is further configured to: determine whether the total load in the second direction is greater than the load threshold in the second direction; if the total load in the second direction is greater than the load threshold in the second direction, calculating a second difference between the total load in the second direction and the load threshold in the second direction, and taking the load threshold in the second direction as a first load in the second direction applied to the hard point position by the first RBE3 unit and the second difference as a second load in the second direction applied to the contact position by the second RBE3 unit; otherwise, taking the total load in the second direction as the first load in the second direction applied to the hard point position by the first RBE3 unit.

[0018] According to an embodiment of the present application, the determining module is further configured to: determine whether the total load in the third direction is greater than the load threshold in the third direction; if the total load in the third direction is greater than the load threshold in the third direction, calculating a third difference between the total load in the third direction and the load threshold in the third direction, and taking the load threshold in the third direction as a first load in the third direction applied to the hard point position by the first RBE3 unit and the third difference as a second load in the third direction applied to the contact position by the second RBE3 unit.

[0019] otherwise, taking the total load in the third direction as the first load in the third direction applied to the hard point position by the first RBE3 unit.

[0020] According to the embodiment of the present application, the general modeling device for load secondary decomposition based on bushing stiffness is used to process the geometric model to be simulated, obtain the size and gap of the rubber bushing in combination with the whole vehicle coordinate system, determine the load critical value according to the nonlinear stiffness curve, and then decompose the load of the hard point and the contact position according to the total load. After the grid division, material assignment and constraint setting, the finite element model is constructed and the simulation result is output. Thus, the problem that the rubber bushing is contained in the current structure but the nonlinear stiffness of the rubber bushing is not considered and the force position of the structure is not correctly determined is solved, and the simulation precision and accuracy of the structure are greatly improved.

[0021] To achieve the above object, the third aspect of the present application provides an electronic device, comprising: a memory, a processor and a computer program stored in the memory and executable on the processor, wherein the processor executes the program to implement the general modeling method for load secondary decomposition based on bushing stiffness as described in the above embodiments.

[0022] To achieve the above object, the fourth aspect of the present application provides a computer readable storage medium having a computer program stored thereon, wherein the program is executed by a processor to implement the general modeling method for load secondary decomposition based on bushing stiffness as described in the above embodiments.

[0023] To achieve the above object, the fifth aspect of the present application provides a computer program product comprising a computer program, wherein the computer program is executed by a processor to implement the general modeling method for load secondary decomposition based on bushing stiffness as described in the above embodiments.

[0024] The additional aspects and advantages of the present application will be partially given in the following description, partially become obvious from the following description, or be known by the practice of the present application. BRIEF DESCRIPTION OF DRAWINGS

[0025] The above and / or additional aspects and advantages of the present application will become apparent and more readily appreciated from the following description of the embodiments, taken in conjunction with the accompanying drawings, in which: Figure 1 A flow chart of a general modeling method for load secondary decomposition based on bushing stiffness according to an embodiment of the present application; Figure 2 A geometric model structure schematic diagram according to an embodiment of the present application; Figure 3 A rubber bushing stiffness parameter and curve schematic diagram according to an embodiment of the present application; Figure 4 A finite element calculation model structure schematic diagram according to an embodiment of the present application; Figure 5A flow chart of a general modeling method based on bushing stiffness for secondary load decomposition is provided according to an embodiment of the present application. Figure 6 A block schematic diagram of a general modeling device based on bushing stiffness for secondary load decomposition is provided according to an embodiment of the present application. Figure 7 A structural schematic diagram of an electronic device is provided according to an embodiment of the present application. DETAILED DESCRIPTION

[0026] Embodiments of the present application are described in detail below with reference to the accompanying drawings, in which the same or similar components have the same or similar designations and functions throughout, and examples of the embodiments described below are illustrative and are intended to explain the present application, and are not to be understood as limiting the present application.

[0027] A general modeling method based on bushing stiffness for secondary load decomposition is described below with reference to the accompanying drawings according to an embodiment of the present application. In view of the problems that the current structural member contains a rubber bushing but does not consider the nonlinear stiffness of the rubber bushing and does not correctly determine the force position of the structural member in the background art, the present application proposes a general modeling method based on bushing stiffness for secondary load decomposition. The method obtains the size and gap amount of the rubber bushing by processing the simulation geometry model and combining the vehicle coordinate system, determines the load critical value according to the nonlinear stiffness curve, and then decomposes the load at the hard point and the contact position according to the total load. After meshing, material assignment and constraint setting, a finite element model is constructed and simulation results are output. The method solves the problems that the current structural member contains a rubber bushing but does not consider the nonlinear stiffness of the rubber bushing and does not correctly determine the force position of the structural member, and greatly improves the simulation accuracy and accuracy of the structural member.

[0028] First, the general modeling method based on bushing stiffness for secondary load decomposition according to an embodiment of the present application is described with reference to the accompanying drawings.

[0029] Figure 1 A flow chart of a general modeling method based on bushing stiffness for secondary load decomposition according to an embodiment of the present application.

[0030] As shown in Figure 1 the general modeling method based on bushing stiffness for secondary load decomposition includes the following steps: In step S101, a simulation geometry model is obtained, and the simulation geometry model is geometrically processed to obtain a processed geometry model.

[0031] The simulation geometry model refers to an original three-dimensional digitized model used for this simulation analysis, which is the object carrier of the simulation. The geometric processing refers to a series of preprocessing of the simulation geometry model. The processed geometry model refers to an optimized model obtained after geometric processing.

[0032] Specifically, the geometric model to be simulated is usually obtained by computer-aided design software or using a three-dimensional scanning device, and there are many ways to obtain it, which will not be described in detail here. The purpose of geometric processing is to optimize the model to make it more suitable for subsequent meshing, load application and other simulation operations, mainly including: deleting redundant features, repairing geometric defects, geometric simplification and geometric segmentation extraction and other operation modes. Here will not be described in detail.

[0033] In step S102, based on the whole vehicle coordinate system, the sizes of the first to third directions of the rubber bushing in the processed geometric model and the gap amounts of the first to third directions between the rubber bushing and the target structural member in the processed geometric model are obtained, and based on the sizes of the first to third directions and the gap amounts of the first to third directions, the load critical values of the first to third directions are determined according to the nonlinear stiffness curve of the rubber bushing.

[0034] Optionally, in some embodiments, determining the load critical values of the first to third directions according to the nonlinear stiffness curve of the rubber bushing based on the sizes of the first to third directions and the gap amounts of the first to third directions comprises: obtaining actual displacement parameters of the rubber bushing in the first to third directions according to the size of each direction and the gap amount of the corresponding direction; determining a set of displacement parameters of the rubber bushing in the first to third directions according to the nonlinear stiffness curve of the rubber bushing; and determining the load critical values of the first to third directions according to the actual displacement parameters of the first to third directions and the set of displacement parameters of the first to third directions.

[0035] Further, in some embodiments, determining the load critical values of the first to third directions according to the actual displacement parameters of the first to third directions and the set of displacement parameters of the first to third directions comprises: determining whether the actual displacement parameter of the first direction is in the set of displacement parameters of the first direction; if the actual displacement parameter of the first direction is in the set of displacement parameters of the first direction, obtaining the load critical value of the first direction based on the set of displacement parameters of the first direction, otherwise, obtaining the load critical value of the first direction by interpolation according to the set of displacement parameters of the first direction based on the actual displacement parameter of the first direction; determining whether the actual displacement parameter of the second direction is in the set of displacement parameters of the second direction; if the actual displacement parameter of the second direction is in the set of displacement parameters of the second direction, obtaining the load critical value of the second direction based on the set of displacement parameters of the second direction, otherwise, obtaining the load critical value of the second direction by interpolation according to the set of displacement parameters of the second direction based on the actual displacement parameter of the second direction; determining whether the actual displacement parameter of the third direction is in the set of displacement parameters of the third direction; if the actual displacement parameter of the third direction is in the set of displacement parameters of the third direction, obtaining the load critical value of the third direction based on the set of displacement parameters of the third direction, otherwise, obtaining the load critical value of the third direction by interpolation according to the set of displacement parameters of the third direction based on the actual displacement parameter of the third direction.

[0036] Wherein, the whole vehicle coordinate system is a unified three-dimensional coordinate system used in the automobile industry to describe the position and movement of various components of a vehicle; the rubber bushing is a commonly used elastic connecting piece, which is widely used in mechanical structures such as automobiles and plays a role of buffering, vibration isolation, and connection; the first to third directions generally correspond to the X, Y, and Z coordinate axis directions of the whole vehicle coordinate system; the size refers to the geometric measurement of the rubber bushing in the first to third directions; the gap amount is the size of the gap between the rubber bushing and the target structural part in the first to third directions; the nonlinear stiffness curve is a curve of the stiffness of the rubber bushing changing with deformation; the load critical value is the key load value that the rubber bushing can withstand in the first to third directions corresponding to a specific displacement; the actual displacement parameter is the actual displacement size that the rubber bushing may generate under the condition of being loaded, which is calculated according to the size of the rubber bushing in each direction and the gap amount with the target structural part; the displacement parameter set is a set of a series of different displacement values and related parameters such as stiffness or load corresponding to the nonlinear stiffness curve of the rubber bushing; interpolation is a mathematical method for estimating the value of an unknown point between known discrete data points.

[0037] Specifically, as shown in Figure 2 , Figure 2 is a geometric model structure diagram provided according to an embodiment of the present application, wherein, Figure 2 A in the above formula represents the gap amount between the rubber bushing and the structure, Figure 2 B in the above formula represents the size of the rubber bushing itself, Figure 2 E in the above formula represents the contact surface after the rubber bushing is projected in section, Figure 2 F in the above formula represents the initial contact surface, the size of the blue rubber bushing in the X, Y, and Z directions under the whole vehicle coordinate system, and the gap amount with the gray structural part. Since the rubber bushing has very small compressibility, it is specified that the maximum compression amount can only be compressed to 1 / 3 of its own size, i.e., it does not deform after being compressed to 1 / 3 of its own size. According to the size of the rubber bushing in section, the contact area of the rubber bushing with the structural part is determined, and the critical value of load secondary decomposition is determined by combining the nonlinear stiffness of the rubber bushing, the gap amount between the rubber bushing and the structural part, and the maximum compression amount of the rubber bushing to redistribute the X, Y, and Z direction loads. First, it is determined whether the rubber bushing will contact the structural part under the condition of being loaded in the +X, -X, +Y, -Y, +Z, and -Z directions. If it contacts, load secondary decomposition is performed. If it does not contact, there is no need to adjust.

[0038] In the case that the rubber bushing is in contact with the structure under load, then (A+B / 3) is used as the displacement parameter of the bushing. If the rubber bushing displacement parameter does not contain (A+B / 3) point, then the load corresponding to the (A+B / 3) point is obtained by linear interpolation of the upper and lower two points closest to the (A+B / 3) point, so as to ensure that the corresponding load C is as much as possible on the nonlinear curve of the bushing, and then C is the critical value of the load quadratic decomposition, as shown in Figure 3 Figure 3 The rubber bushing stiffness parameter and the curve diagram provided according to an embodiment of the application, wherein the text data on one side presents the parameter values related to the stiffness of the rubber bushing in different directions, which are the basic data support for determining the nonlinear stiffness curve of the rubber bushing and calculating the load critical value, and the other side is the stiffness curve, wherein the horizontal coordinate represents the displacement of the rubber bushing, and the vertical coordinate represents the force borne by the rubber bushing. Different color curves respectively represent the force-displacement relationship of the rubber bushing in different directions. As can be seen from the curve shape, the stiffness of the rubber bushing has nonlinear characteristics, that is, the force and displacement are not in a simple linear proportional relationship, and the change rate of the force will change with the displacement. The bearing capacity change rule of the rubber bushing in different displacement stages is complex and nonlinear.

[0039] Therefore, based on the whole vehicle coordinate system, in combination with the size of the rubber bushing, the gap amount with the target structure and the nonlinear stiffness curve, the load critical value in each direction is determined through the relationship between the actual displacement parameter and the displacement parameter set (including judging whether it is in the set or not, and interpolating if not), which fully considers the nonlinear stiffness characteristics of the rubber bushing and breaks through the limitation of the linear assumption in the conventional modeling. By accurately determining the load critical value and reasonably decomposing the load, the actual stress state of the structure under the limiting action of the rubber bushing can be truly reflected, the unreasonable problem that all contact surfaces are stressed in the traditional model is solved, and the accuracy and reliability of the finite element simulation are significantly improved, thereby providing a more actual working condition basis for the strength analysis of the structure.

[0040] In step S103, based on the total load in the first to third directions to be calculated, the first load in the first to third directions applied to the hard point position through the first RBE3 unit and the second load in the first to third directions applied to the contact position of the rubber bushing and the target structure through the second RBE3 unit are determined according to the load critical values in the first to third directions.

[0041] ​Optionally, in some embodiments, based on the total load of the first to third directions to be calculated, the first to third directions' first loads applied to the hard point position by the first RBE3 unit and the first to third directions' second loads applied to the contact position of the rubber bushing and the target structural member by the second RBE3 unit are determined according to the load critical values of the first to third directions, comprising: judging whether the total load of the first direction is greater than the load critical value of the first direction; if the total load of the first direction is greater than the load critical value of the first direction, calculating a first difference value between the total load of the first direction and the load critical value of the first direction, and taking the load critical value of the first direction as the first load of the first direction applied to the hard point position by the first RBE3 unit, and taking the first difference value as the second load of the first direction applied to the contact position by the second RBE3 unit; otherwise, taking the total load of the first direction as the first load of the first direction applied to the hard point position by the first RBE3 unit.

[0042] Optionally, in some embodiments, based on the total load of the first to third directions to be calculated, the first to third directions' first loads applied to the hard point position by the first RBE3 unit and the first to third directions' second loads applied to the contact position of the rubber bushing and the target structural member by the second RBE3 unit are determined according to the load critical values of the first to third directions, comprising: judging whether the total load of the second direction is greater than the load critical value of the second direction; if the total load of the second direction is greater than the load critical value of the second direction, calculating a second difference value between the total load of the second direction and the load critical value of the second direction, and taking the load critical value of the second direction as the first load of the second direction applied to the hard point position by the first RBE3 unit, and taking the second difference value as the second load of the second direction applied to the contact position by the second RBE3 unit; otherwise, taking the total load of the second direction as the first load of the second direction applied to the hard point position by the first RBE3 unit.

[0043] Optionally, in some embodiments, based on the total load in the first to third directions to be calculated, the first load in the first to third directions applied to the hard point position by the first RBE3 unit and the second load in the first to third directions applied to the contact position of the rubber bushing and the target structure by the second RBE3 unit according to the load threshold value in the first to third directions are determined, including: determining whether the total load in the third direction is greater than the load threshold value in the third direction; if the total load in the third direction is greater than the load threshold value in the third direction, calculating a third difference value of the total load in the third direction and the load threshold value in the third direction, and taking the load threshold value in the third direction as the first load in the third direction applied to the hard point position by the first RBE3 unit, and taking the third difference value as the second load in the third direction applied to the contact position by the second RBE3 unit. Otherwise, the total load in the third direction is taken as the first load in the third direction applied to the hard point position by the first RBE3 unit.

[0044] Wherein, the total load in the first to third directions to be calculated refers to the total force borne by the structure in the X, Y, Z directions of the vehicle coordinate system; the load threshold value in the first to third directions is the maximum load threshold value that the rubber bushing can bear in the X, Y, Z directions; the first RBE3 unit is a rigid element used to simulate force or motion transmission in finite element analysis; the hard point position is a key connection point or force reference point on the structure; the first load is the load applied to the hard point position by the first RBE3 unit, representing the load part that the rubber bushing itself can bear; the second RBE3 unit is another rigid connection unit used to transmit the load exceeding the bearing capacity of the rubber bushing to the contact position, simulating the state that the contact structure participates in force after the rubber bushing reaches the limit; the contact position is the part where the rubber bushing and the target structure may contact, which is the load application position of the second RBE3 unit and bears the load part exceeding the threshold value; the second load is the load applied to the contact position by the second RBE3 unit, which exists only when the total load exceeds the load threshold value, is equal to the difference between the total load and the threshold value, and represents the additional load borne by the contact structure; the first / second / third difference value refers to the difference between the total load and the load threshold value in the corresponding direction, which is used to quantify the load part exceeding the bearing capacity of the rubber bushing and is the source of the value of the second load.

[0045] The load acting on the structure in the vehicle coordinate system is shown in Table 1: Table 1

[0046] Wherein, Fx, Fy, Fz respectively represent the forces (unit: Newton, N) borne by the structural member in the X-axis, Y-axis, and Z-axis directions of the vehicle coordinate system; F1, F2, and F3 are respectively the values of the loads borne by the structural member in the X-axis, Y-axis, and Z-axis directions, for quantifying the magnitude of the forces borne by the structural member in different directions, and providing basic data for performance analysis of the structural member such as strength and stiffness.

[0047] Specifically, as shown in Figure 4 , Figure 4 FIG. 1 is a structural schematic diagram of a finite element calculation model provided according to an embodiment of the present application. For the load F1 in the X direction of the vehicle coordinate system, if F1 is greater than the critical value C, a load with a size of (F1-C) is applied to the position where the compressed rubber bushing contacts the structural member, and a load with a size of C is directly applied to the hard point position and is applied to the position where the uncompressed rubber contacts the structural member through the RBE3 unit; if the load F1 in the X direction of the vehicle coordinate system is less than the critical value C, F1 is directly applied to the hard point position and is applied to the position where the uncompressed rubber bushing contacts the structural member through the RBE3 unit. The same applies to the loads F2 and F3 in the Y and Z directions of the vehicle coordinate system. Details are not described herein.

[0048] Thus, for the X, Y, and Z directions, it is first determined whether the total load exceeds the load critical value in the direction. If it exceeds, the critical value part is taken as the first load applied to the hard point position through the first RBE3 unit, borne by the rubber bushing, and the difference between the total load and the critical value is taken as the second load applied to the contact position through the second RBE3 unit, borne by the contact structure; if it does not exceed, the total load is all taken as the first load applied to the hard point position through the first RBE3 unit. This distribution mode can truly reflect the force state of the rubber bushing under different loads, i.e., only the rubber bushing bears the force when it does not reach the limit, and the contact structure participates in the force bearing when it exceeds the limit.

[0049] In step S104, the processed geometric model is meshed and assigned with material properties, and based on the preset boundary constraint condition, the geometric model after meshing is processed to obtain a new finite element model according to the first load in the first to third directions and the second load in the first to third directions, and the simulation result is output based on the new finite element model.

[0050] The dividing grid refers to a key step in finite element analysis, that is, dividing a continuous geometric model into a large number of small, interconnected units to form a discretized grid structure; the material attribute refers to the physical and mechanical characteristic parameters of the material constituting the model; the preset boundary constraint condition can be a boundary constraint condition preset by the user, a boundary constraint condition obtained through a finite number of experiments, or a boundary constraint condition obtained through a finite number of computer simulations, which is not limited here; the first load in the first to third directions refers to the load applied at the hard point position through the first RBE3 unit, representing the load part that the rubber bushing itself can bear; the second load in the first to third directions refers to the load applied at the contact position of the rubber bushing and the target structure through the second RBE3 unit, which exists only when the total load exceeds the load threshold, representing the load part that exceeds the bearing capacity of the rubber bushing and needs to be borne by the contact structure; the geometric model after grid division is a discretized model after grid division processing, composed of a large number of units, and is the direct object of subsequent load application, constraint and finite element calculation; the new finite element model integrates the grid, material attribute, boundary constraint condition and first and second loads to form a complete finite element analysis model; the simulation result is the calculation result output after running the new finite element model. Specifically, the geometric model is divided into a grid, material attributes are assigned, and new boundary conditions and loads are applied to the finite element model. The finite element analysis software is submitted for calculation, and the simulation result is output. Specifically, the model after geometric processing is divided into a grid, that is, the continuous geometric model is divided into a large number of small units to form a discretized grid structure. The first and second loads in the first to third directions are applied to the grid model according to the corresponding directions and positions. The new finite element model is run, and the stress distribution, deformation, displacement trend, etc. of the structure are calculated by the solver to obtain the simulation results. These results can be used to evaluate the strength performance of the structure, verify the rationality of the design, and provide data support for optimization and improvement.

[0051] In summary, the specific process of the general modeling method based on the load secondary decomposition of the bushing stiffness is as shown in Figure 5 Figure 5 The general modeling method based on the load secondary decomposition of the bushing stiffness according to one embodiment of the present application is provided, and the steps are as follows: S501, start finite element simulation; S502, pre-process the geometric model for simulation, including repairing geometric defects, deleting redundant features, simplifying complex structures, etc. geometric processing; S503, determine the size of the bushing, the gap between the structure and the stiffness; ​S504, determine the critical value of the load secondary decomposition according to the size, gap amount and stiffness of the rubber bushing, and the action position of different parts after the load decomposition; S505, judge whether the original total load is greater than the critical value determined in step S504, if greater than the critical value, execute step S506, otherwise, execute step S508; S506, re-decompose the load; S507, divide the geometric model into grids and assign material properties; S508, apply the preset boundary constraint condition and the decomposed or undecomposed load to the model with divided grids and assigned material properties, so that the model simulates the actual stress and constraint state; S509, submit the finite element software calculation; S510, output the simulation result.

[0052] Therefore, the whole process can effectively simulate the stress and response of the structure containing the rubber bushing by accurate model processing, load analysis and decomposition, combined with the finite element method, to provide a basis for structure performance evaluation and optimization.

[0053] According to the load secondary decomposition general modeling method based on bushing stiffness provided by the embodiment of the application, the size and gap amount of the rubber bushing are obtained by processing the to-be-simulated geometric model and combining the vehicle coordinate system, the load critical value is determined according to the nonlinear stiffness curve of the rubber bushing, and then the load of the hard point and the contact position is decomposed from the total load. After grid division, material assignment and constraint setting, the finite element model is constructed and the simulation result is output. Therefore, the problem that the rubber bushing is contained in the current structure but the nonlinear stiffness of the rubber bushing is not considered and the stress position of the structure is not correctly determined is solved, and the simulation precision and accuracy of the structure are greatly improved.

[0054] Secondly, the load secondary decomposition general modeling device based on bushing stiffness according to the embodiment of the application is described with reference to the accompanying drawings.

[0055] Figure 6 is a block schematic diagram of the load secondary decomposition general modeling device based on bushing stiffness according to an embodiment of the application.

[0056] As shown in Figure 6 , the load secondary decomposition general modeling device based on bushing stiffness 10 comprises a first acquisition module 100, a second acquisition module 200, a determination module 300 and a decomposition module 400.

[0057] The first acquisition module 100 is configured to acquire a to-be-simulated geometric model, and perform geometric processing on the to-be-simulated geometric model to obtain a processed geometric model. The second acquisition module 200 acquires, based on the vehicle coordinate system, the sizes of the rubber bushing in the first to third directions in the processed geometric model and the gap amounts in the first to third directions between the rubber bushing and the target structural member in the processed geometric model, and determines the load critical values in the first to third directions according to the nonlinear stiffness curve of the rubber bushing based on the sizes in the first to third directions and the gap amounts in the first to third directions. The determination module 300 determines, based on the total load in the first to third directions to be calculated, the first load in the first to third directions applied to the hard point position through the first RBE3 unit and the second load in the first to third directions applied to the contact position of the rubber bushing and the target structural member through the second RBE3 unit according to the load critical values in the first to third directions. The decomposition module 400 divides the processed geometric model into grids and assigns material properties, and processes the geometric model after grid division to obtain a new finite element model based on the preset boundary constraint condition according to the first load in the first to third directions and the second load in the first to third directions, and outputs a simulation result based on the new finite element model.

[0058] According to an embodiment of the present application, the second acquisition module 200 is specifically configured to: obtain actual displacement parameters of the rubber bushing in the first to third directions according to the size in each direction and the gap amount in the corresponding direction; determine the displacement parameter set of the rubber bushing in the first to third directions according to the nonlinear stiffness curve of the rubber bushing; and determine the load critical values in the first to third directions according to the actual displacement parameters in the first to third directions and the displacement parameter set in the first to third directions.

[0059] According to an embodiment of the present application, the second acquisition module 200 is further configured to: determine whether the actual displacement parameter in the first direction is in the displacement parameter set in the first direction; if the actual displacement parameter in the first direction is in the displacement parameter set in the first direction, obtain the load critical value in the first direction based on the displacement parameter set in the first direction, otherwise, obtain the load critical value in the first direction by interpolation according to the displacement parameter set in the first direction based on the actual displacement parameter in the first direction; determine whether the actual displacement parameter in the second direction is in the displacement parameter set in the second direction; if the actual displacement parameter in the second direction is in the displacement parameter set in the second direction, obtain the load critical value in the second direction based on the displacement parameter set in the second direction, otherwise, obtain the load critical value in the second direction by interpolation according to the displacement parameter set in the second direction based on the actual displacement parameter in the second direction; determine whether the actual displacement parameter in the third direction is in the displacement parameter set in the third direction; if the actual displacement parameter in the third direction is in the displacement parameter set in the third direction, obtain the load critical value in the third direction based on the displacement parameter set in the third direction, otherwise, obtain the load critical value in the third direction by interpolation according to the displacement parameter set in the third direction based on the actual displacement parameter in the third direction.

[0060] According to one embodiment of the present application, the determining module 300 is specifically configured to: determine whether the total load in the first direction is greater than the load threshold in the first direction; if the total load in the first direction is greater than the load threshold in the first direction, calculate a first difference between the total load in the first direction and the load threshold in the first direction, and take the load threshold in the first direction as the first load in the first direction applied to the hard point position by the first RBE3 unit and take the first difference as the second load in the first direction applied to the contact position by the second RBE3 unit; otherwise, take the total load in the first direction as the first load in the first direction applied to the hard point position by the first RBE3 unit.

[0061] According to one embodiment of the present application, the determining module 300 is further configured to: determine whether the total load in the second direction is greater than the load threshold in the second direction; if the total load in the second direction is greater than the load threshold in the second direction, calculate a second difference between the total load in the second direction and the load threshold in the second direction, and take the load threshold in the second direction as the first load in the second direction applied to the hard point position by the first RBE3 unit and take the second difference as the second load in the second direction applied to the contact position by the second RBE3 unit; otherwise, take the total load in the second direction as the first load in the second direction applied to the hard point position by the first RBE3 unit.

[0062] According to one embodiment of the present application, the determining module 300 is further configured to: determine whether the total load in the third direction is greater than the load threshold in the third direction; if the total load in the third direction is greater than the load threshold in the third direction, calculate a third difference between the total load in the third direction and the load threshold in the third direction, and take the load threshold in the third direction as the first load in the third direction applied to the hard point position by the first RBE3 unit and take the third difference as the second load in the third direction applied to the contact position by the second RBE3 unit; otherwise, take the total load in the third direction as the first load in the third direction applied to the hard point position by the first RBE3 unit.

[0063] It should be noted that the foregoing explanation of the embodiment of the general modeling method for load secondary decomposition based on bushing stiffness is also applicable to the embodiment of the general modeling device for load secondary decomposition based on bushing stiffness, which will not be repeated here.

[0064] According to the load secondary decomposition general modeling device based on the bushing stiffness provided by the embodiment of the application, the size and gap amount of the rubber bushing are obtained by processing a to-be-simulated geometric model in combination with a whole vehicle coordinate system, the load critical value is determined according to the nonlinear stiffness curve, and then the load of the hard point and the contact position is decomposed according to the total load, and the finite element model is constructed and the simulation result is output after grid division, material assignment and constraint setting. Therefore, the problem that the rubber bushing is contained in the current structure but the nonlinear stiffness of the rubber bushing is not considered and the force position of the structure is not correctly determined is solved, and the simulation precision and accuracy of the structure are greatly improved.

[0065] Figure 7 The structure schematic diagram of the electronic equipment provided by the embodiment of the application is provided. The electronic equipment can include: The memory 701, the processor 702 and the computer program stored in the memory 701 and executable on the processor 702.

[0066] The processor 702 implements the load secondary decomposition general modeling method based on the bushing stiffness provided in the above embodiment when executing the program.

[0067] Further, the electronic equipment further includes: The communication interface 703 is used for communication between the memory 701 and the processor 702.

[0068] The memory 701 is used for storing the computer program executable on the processor 702.

[0069] The memory 701 can include a high-speed RAM (Random Access Memory, random access memory) memory, and can also include a non-volatile memory, for example, at least one disk memory.

[0070] If the memory 701, the processor 702 and the communication interface 703 are independently implemented, the communication interface 703, the memory 701 and the processor 702 can be connected to each other through a bus and complete the communication between each other. The bus can be an ISA (Industry Standard Architecture, industry standard architecture) bus, a PCI (Peripheral Component Interconnect, peripheral component interconnect) bus or an EISA (Extended Industry Standard Architecture, extended industry standard architecture) bus, etc. The bus can be divided into an address bus, a data bus, a control bus, etc. For the convenience of representation, Figure 7 Only one thick line is used in the figure, but it does not mean that there is only one bus or one type of bus.

[0071] Optionally, if the memory 701, the processor 702 and the communication interface 703 are integrated on a chip, the memory 701, the processor 702 and the communication interface 703 can complete the communication among each other through an internal interface.

[0072] The processor 702 can be a CPU (Central Processing Unit), or an ASIC (Application Specific Integrated Circuit), or one or more integrated circuits configured to implement one or more embodiments of the application.

[0073] The embodiment of the application further provides a computer readable storage medium, which has a computer program stored thereon, and the computer program is executed by a processor to implement the load secondary decomposition general modeling method based on bushing stiffness as above.

[0074] The embodiment of the application further provides a computer program product, which comprises a computer program, and the computer program is executed by a processor to implement the load secondary decomposition general modeling method based on bushing stiffness as above.

[0075] In addition, the terms "first", "second", "third", etc. are used herein only to describe different instances, and cannot be construed to indicate or imply relative importance or imply the number of the indicated technical features. Therefore, the features defined with "first", "second" can explicitly or implicitly include at least one of the features. In the description of the present application, the meaning of "a plurality of" is at least two, for example, two, three, etc., unless otherwise specifically limited.

[0076] In the description of the present application, the description of the terms "one embodiment", "some embodiments", "example", "specific example", or "some examples" means that the specific features, structures, materials or characteristics described in connection with the embodiment or example are included in at least one embodiment or example of the present application. In the present application, the illustrative description of the above terms is not necessarily for the same embodiment or example. Moreover, the specific features, structures, materials or characteristics described can be combined in any appropriate manner in any one or more embodiments or examples. In addition, the skilled in the art can combine and combine the different embodiments or examples described in the present application and the features of the different embodiments or examples without contradiction.

[0077] Although the embodiments of the present application have been shown and described above, it should be understood that the above embodiments are exemplary and cannot be construed as limiting the present application, and those skilled in the art can make changes, modifications, replacements and variations to the above embodiments within the scope of the present application.

Claims

1. A bushing stiffness-based load secondary path generalized modeling method, characterized by, The method comprises the following steps: acquiring a to-be-simulated geometric model, and performing geometric processing on the to-be-simulated geometric model to obtain a processed geometric model; based on a whole vehicle coordinate system, acquiring sizes of a rubber bushing in first to third directions in the processed geometric model and gap amounts in the first to third directions between the rubber bushing and a target structural member in the processed geometric model, and based on the sizes in the first to third directions and the gap amounts in the first to third directions, determining load critical values in the first to third directions according to a nonlinear stiffness curve of the rubber bushing; based on total loads in the first to third directions to be calculated, determining first loads in the first to third directions applied to a hard point position through a first RBE3 unit and second loads in the first to third directions applied to a contact position of the rubber bushing and the target structural member through a second RBE3 unit according to the load critical values in the first to third directions; performing meshing on the processed geometric model and assigning material properties, and based on preset boundary constraint conditions, processing the geometric model after meshing to obtain a new finite element model according to the first loads in the first to third directions and the second loads in the first to third directions, and outputting a simulation result based on the new finite element model.

2. The method of claim 1, wherein, The method comprises the following steps: obtaining actual displacement parameters of the rubber bushing in the first to third directions according to the size of each direction and the gap amount in the corresponding direction; determining a set of displacement parameters of the rubber bushing in the first to third directions according to the nonlinear stiffness curve of the rubber bushing; determining the load critical values in the first to third directions according to the actual displacement parameters in the first to third directions and the set of displacement parameters in the first to third directions.

3. The method of claim 2, wherein, The method comprises the following steps: determining whether the actual displacement parameter in the first direction is in the set of displacement parameters in the first direction; if the actual displacement parameter in the first direction is in the set of displacement parameters in the first direction, obtaining the load critical value in the first direction based on the set of displacement parameters in the first direction, otherwise, obtaining the load critical value in the first direction by interpolation according to the set of displacement parameters in the first direction based on the actual displacement parameter in the first direction; determining whether the actual displacement parameter in the second direction is in the set of displacement parameters in the second direction; if the actual displacement parameter in the second direction is in the set of displacement parameters in the second direction, obtaining the load critical value in the second direction based on the set of displacement parameters in the second direction, otherwise, obtaining the load critical value in the second direction by interpolation according to the set of displacement parameters in the second direction based on the actual displacement parameter in the second direction; determining whether the actual displacement parameter in the third direction is in the set of displacement parameters in the third direction; If the actual displacement parameter of the third direction is in the displacement parameter set of the third direction, a load critical value of the third direction is obtained based on the displacement parameter set of the third direction, otherwise, the load critical value of the third direction is obtained by interpolation based on the actual displacement parameter of the third direction and the displacement parameter set of the third direction.

4. The method of claim 1, wherein, The first to third loads applied to the hard point position by the first RBE3 unit and the first to third loads applied to the contact position of the rubber bushing and the target structure by the second RBE3 unit are determined based on the total load of the first to third directions to be calculated and the load critical values of the first to third directions, comprising: determining whether the total load of the first direction is greater than the load critical value of the first direction; If the total load of the first direction is greater than the load critical value of the first direction, a first difference between the total load of the first direction and the load critical value of the first direction is calculated, the load critical value of the first direction is taken as the first load of the first direction applied to the hard point position by the first RBE3 unit, and the first difference is taken as the second load of the first direction applied to the contact position by the second RBE3 unit; Otherwise, the total load of the first direction is taken as the first load of the first direction applied to the hard point position by the first RBE3 unit.

5. The method of claim 1, wherein, The first to third loads applied to the hard point position by the first RBE3 unit and the first to third loads applied to the contact position of the rubber bushing and the target structure by the second RBE3 unit are determined based on the total load of the first to third directions to be calculated and the load critical values of the first to third directions, comprising: determining whether the total load of the second direction is greater than the load critical value of the second direction; If the total load of the second direction is greater than the load critical value of the second direction, a second difference between the total load of the second direction and the load critical value of the second direction is calculated, the load critical value of the second direction is taken as the first load of the second direction applied to the hard point position by the first RBE3 unit, and the second difference is taken as the second load of the second direction applied to the contact position by the second RBE3 unit; Otherwise, the total load of the second direction is taken as the first load of the second direction applied to the hard point position by the first RBE3 unit.

6. The method of claim 1, wherein, The first to third loads applied to the hard point position by the first RBE3 unit and the first to third loads applied to the contact position of the rubber bushing and the target structure by the second RBE3 unit are determined based on the total load of the first to third directions to be calculated and the load critical values of the first to third directions, comprising: determining whether the total load of the third direction is greater than the load critical value of the third direction; If the total load in the third direction is greater than the load threshold in the third direction, a third difference between the total load in the third direction and the load threshold in the third direction is calculated, the load threshold in the third direction is taken as a first load in the third direction applied to the hard point position by the first RBE3 unit, and the third difference is taken as a second load in the third direction applied to the contact position by the second RBE3 unit; Otherwise, the total load in the third direction is taken as the first load in the third direction applied to the hard point position by the first RBE3 unit.

7. A bushing stiffness-based load secondary path generalized modeling apparatus, characterized by, The method comprises: The first obtaining module is configured to obtain a to-be-simulated geometric model, and perform geometric processing on the to-be-simulated geometric model to obtain a processed geometric model; The second obtaining module is configured to obtain, based on a whole-vehicle coordinate system, sizes in first to third directions of a rubber bushing in the processed geometric model and first to third direction gap amounts between the rubber bushing and target structural members in the processed geometric model, and determine, based on the sizes in the first to third directions and the first to third direction gap amounts, load thresholds in the first to third directions according to a nonlinear stiffness curve of the rubber bushing; The determining module is configured to determine, based on total loads in the first to third directions to be calculated, first loads in the first to third directions applied to hard point positions by first RBE3 units and second loads in the first to third directions applied to contact positions of the rubber bushing and the target structural members by second RBE3 units according to the load thresholds in the first to third directions. The decomposing module is configured to divide the processed geometric model into grids, assign material properties, and process the geometric model after grid division based on preset boundary constraint conditions according to the first to third loads and the second to third loads to obtain a new finite element model, and output a simulation result based on the new finite element model.

8. An electronic device, comprising: The method comprises: A memory, a processor, and a computer program stored on the memory and executable on the processor, wherein the processor executes the program to implement the general modeling method based on secondary decomposition of bushing stiffness loads according to any one of claims 1-6.

9. A computer readable storage medium having stored thereon a computer program, characterized in that, The program is executed by the processor to implement the general modeling method based on secondary decomposition of bushing stiffness loads according to any one of claims 1-6.

10. A computer program product comprising a computer program, characterized in that, The computer program is executed by the processor to implement the general modeling method based on secondary decomposition of bushing stiffness loads according to any one of claims 1-6.

Images