CN114970237B - Method for improving fatigue durability simulation precision of stabilizer bar system - Google Patents
Method for improving fatigue durability simulation precision of stabilizer bar system Download PDFInfo
- Publication number
- CN114970237B CN114970237B CN202210424742.1A CN202210424742A CN114970237B CN 114970237 B CN114970237 B CN 114970237B CN 202210424742 A CN202210424742 A CN 202210424742A CN 114970237 B CN114970237 B CN 114970237B
- Authority
- CN
- China
- Prior art keywords
- stabilizer bar
- rubber
- simulation
- bar system
- cover plate
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 239000003381 stabilizer Substances 0.000 title claims abstract description 78
- 238000004088 simulation Methods 0.000 title claims abstract description 48
- 238000000034 method Methods 0.000 title claims abstract description 43
- 238000006073 displacement reaction Methods 0.000 claims abstract description 28
- 239000000463 material Substances 0.000 claims abstract description 23
- 238000012360 testing method Methods 0.000 claims abstract description 22
- 238000004364 calculation method Methods 0.000 claims abstract description 21
- 238000003825 pressing Methods 0.000 claims abstract description 8
- 230000000452 restraining effect Effects 0.000 claims abstract description 3
- 238000005457 optimization Methods 0.000 claims description 16
- 238000013519 translation Methods 0.000 claims description 9
- 230000035515 penetration Effects 0.000 claims description 6
- 238000006243 chemical reaction Methods 0.000 claims description 3
- 238000011156 evaluation Methods 0.000 claims description 3
- 238000012545 processing Methods 0.000 claims description 3
- 230000005540 biological transmission Effects 0.000 abstract description 4
- 238000011161 development Methods 0.000 description 2
- 229910000831 Steel Inorganic materials 0.000 description 1
- 230000002159 abnormal effect Effects 0.000 description 1
- 239000006096 absorbing agent Substances 0.000 description 1
- 239000000853 adhesive Substances 0.000 description 1
- 230000001070 adhesive effect Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 239000011248 coating agent Substances 0.000 description 1
- 238000000576 coating method Methods 0.000 description 1
- 238000005094 computer simulation Methods 0.000 description 1
- 230000005577 local transmission Effects 0.000 description 1
- 238000010422 painting Methods 0.000 description 1
- 230000008707 rearrangement Effects 0.000 description 1
- 238000005096 rolling process Methods 0.000 description 1
- 230000035939 shock Effects 0.000 description 1
- 239000010959 steel Substances 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 239000002699 waste material Substances 0.000 description 1
Abstract
The invention relates to a method for improving fatigue durability simulation precision of a stabilizer bar system, which comprises finite element modeling; calibrating parameters of the rubber material; applying a bolt pre-tightening load, pressing and connecting the cover plate to a local structure of the auxiliary frame, and simultaneously pressing and assembling the cover plate and the rubber bushing together; performing simulation calculation of external load fatigue conditions, restraining local structures of auxiliary frames, adding relative displacement at two ends of a stabilizer bar, and solving to obtain two stress field results; and (5) leading in two stress fields to finish the calculation of the constant-amplitude fatigue endurance damage value between the two stress fields. According to the invention, the contact constraint among the rubber bushing, the stabilizer bar and the cover plate is established, so that the load transmission is ensured to be accurate; setting a simulation assembly process through bolt pre-tightening and contact constraint initial interference; the real rubber bushing model is used for replacing a virtual rubber bushing, so that the problem that the rigidity test data of the bushing does not meet the finite element simulation requirement is solved; through the measures, three difficult problems in the simulation process of the stabilizer bar system are solved, and the simulation precision is improved.
Description
Technical Field
The invention belongs to the technical field of stabilizer bar systems, and particularly relates to a method for improving fatigue durability simulation precision of a stabilizer bar system.
Background
The transverse stabilizer bar system is used for preventing the vehicle body from excessively transversely rolling when turning so as to ensure the stability and smoothness of the vehicle. In the stabilizer bar system, the both ends of stabilizer bar are connected on control arm or bumper shock absorber support through the connecting rod, and the straight section in centre has controlling two points to pass through bushing connection on the sub vehicle frame, and the bush is built-in rubber spare, and the coating adhesive is and vulcanizes into an organic whole between stabilizer bar and the stabilizer bar bush, avoids producing wearing and tearing and abnormal sound in the use. When the automobile turns, the side tilting of the automobile body is caused by the action of centrifugal force, so that the bent inner wheel is suspended and stretched, the bent outer wheel is suspended and compressed, the stabilizer bar body is twisted, and great torsion is generated. A stabilizer bar of an automobile in the whole life cycle can be twisted for 30 ten thousand times approximately, so that the stabilizer bar is designed to meet the fatigue life of 30 ten thousand times twisting, and meanwhile, the stabilizer bar cannot exceed a target value too much, otherwise, the stabilizer bar is over-designed, and material waste and cost increase are caused.
Fatigue durability of stabilizer bar systems is an important performance indicator for vehicle development. There are two methods for predicting the fatigue life of the stabilizer bar at present, one is through bench test and the other is through finite element simulation calculation. The former has high precision, but has longer cost and period; the latter is low in cost and short in period, but the simulation accuracy is greatly different due to the difference of the simulation methods. The finite element modeling of the existing stabilizer bar system has three difficulties: firstly, bushing assembly simulation needs to consider interference assembly stress of the bushing, and the assembly process has a large influence on the strength of the stabilizer bar cover plate; secondly, the simulation of the position connection relation of the bushing needs to consider the contact of the bushing, the stabilizer bar and the cover plate, so that the accuracy of local load transmission is ensured; thirdly, the rigidity data of the bushing does not meet the finite element simulation requirement, and the nonlinear interval of the rigidity test data of the general bushing is shorter, so that the problem of poor convergence exists.
At present, when a stabilizer bar system is simulated, the positions of rubber bushings are connected by virtual bushing units, one ends of the bushings are connected to a stabilizer bar through rb3 units, the other ends of the bushings are connected to a stabilizer bar cover plate through rb3 units, the bushings are endowed with various rigidity to simulate real bushings, when loads are transmitted through rb3, the load transmission of each connecting point follows a certain weight coefficient, the weight coefficient is difficult to determine, and therefore the local transmission of the loads is inaccurate. The method cannot evaluate the connecting position of the stabilizer bar, and in the past experience, the failure position is often concentrated at the connecting position of the bushing, so that the simulation precision is low. Meanwhile, the interference assembly process of the bushing cannot be considered, and the assembly stress cannot be substituted into calculation, so that the simulation precision is further reduced. Based on the above, there is a strong need to develop a method capable of effectively improving the fatigue endurance simulation precision of the stabilizer bar system.
Disclosure of Invention
The invention aims to provide a method for improving fatigue endurance simulation precision of a stabilizer bar system, which is based on real system modeling and considers the influence of an assembly process at the same time so as to solve the problem of lower finite element fatigue simulation precision of the existing stabilizer bar system. The method can improve the simulation precision and save the development cost.
The invention aims at realizing the following technical scheme:
a method for improving fatigue durability simulation precision of a stabilizer bar system comprises the following steps:
A. finite element modeling
A1, sequentially drawing grids on parts of a stabilizer bar system, wherein the parts of the stabilizer bar system consist of stabilizer bars, a plurality of rubber bushings, a plurality of stabilizer bar cover plates and a plurality of bolts;
a2, establishing material properties;
A3, endowing the unit with attributes, wherein all structures use entity units;
A4, assembling a model;
A5, establishing constraint and load working conditions; B. the rubber material constitutive parameters are calibrated by calculating initial constitutive parameters of rubber through rubber hardness, and then the actual constitutive parameters of the rubber are calibrated through parameter optimization means;
C. Assembly process simulation
Applying a bolt pre-tightening load, pressing and connecting the cover plate to a local structure of the auxiliary frame, and simultaneously pressing and assembling the cover plate and the rubber bushing together;
D. Durable condition calculation
C, based on the step, performing simulation calculation of external load fatigue working conditions, restraining local structures of auxiliary frames, adding certain relative displacement at two ends of a stabilizer bar, and solving to obtain two stress field results;
E. Durable damage value calculation
And D, importing the results of the two stress fields calculated in the step D, and finishing calculation of the constant-amplitude fatigue endurance damage value between the two stress fields.
Further, in the step A1, the number of the rubber bushings and the stabilizer bar cover plates is 2, and the number of bolts is 4.
Further, in the step A1, the rubber bushing is divided into hexahedral grids, the aspect ratio is ensured to be more than or equal to 5, and other part grids are divided into second-order tetrahedral grids.
Further, step A2 specifically includes: the rubber body structure material adopts a two-term Mooney-Rivlin model, namely a formula (1), initial values of rubber constitutive parameters C 10 and C 01 are calculated through the input of rubber hardness of a material constitutive parameter empirical formula (2) -a formula (4), and the density of the material is endowed;
W＝C10(I1-3)+C01(I2-3) (1)
lgE0＝0.0184Hr-0.4575 (2)
G＝E0/3＝2(C10+C01) (3)
C01/C10＝0.05 (4)
wherein C 10,C01 is a material constant; e 0 is the elastic modulus of the rubber, H r is the Rockwell hardness of the rubber, and G is the shear modulus.
Further, step A4 is specifically: connecting bolt holes at two ends of the stabilizer bar by using a rigid unit, wherein a main point is a center point of the bolt holes; the bush and the stabilizer bar are used for tie constraint, and the bush and the cover plate, the bolt and the cover plate and the auxiliary frame are used for local structure contact constraint, wherein the bush and the cover plate structure have initial penetration, and the penetration is interference.
Further, step A5 specifically includes: the auxiliary frame partial structure is restrained, the load working condition comprises two load steps, and the first step is the assembly working condition which comprises bolt pre-tightening and interference fit of a bushing and a cover plate; and the second step is an external load working condition, wherein the external load is the relative displacement quantity of the two ends of the stabilizer bar, the relative displacement quantity is related to the left-right relative runout quantity of the wheel center, and the two ends of the stabilizer bar are converted through the lever ratio.
Further, step B specifically comprises:
b1, ensuring that the displacement output interval of a calculated stiffness curve is consistent with that of a test stiffness curve through result output and data processing, restricting the loading of an outer tube and an inner tube of a bushing under the same simulation working condition as that of a rubber bushing rack stiffness test, and solving a load and displacement curve;
B2, calculating the minimum value of the area S surrounded by the stiffness curve and the test stiffness curve by taking C 10 and C 01 as variables as targets, and calibrating the rubber constitutive parameters;
B3, parameter optimization: firstly preparing a test stiffness curve, then preparing a stiffness calculation finite element model, constructing parameter optimization flow iterative calculation through parameterization C 10 and C 01 and target setting S, executing parameter optimization iterative calculation, and finally determining parameters C 10 and C 01.
Further, the specific steps of calculating the area S surrounded by the stiffness curve and the test stiffness curve are as follows:
S1＝(a1-b1)2+(a2-b2)2+…(an-bn)2 (5)
S2＝(c1-d1)2+(c2-d2)2+…(cn-dn)2 (6)
S3＝(e1-f1)2+(e2-f2)2+…(en-fn)2 (7)
S4＝(g1-h1)2+(g2-h2)2+…(gn-hn)2 (8)
S5＝(i1-j1)2+(i2-j2)2+…(in-jn)2 (9)
S6＝(k1-l1)2+(k2-l2)2+…(kn-ln)2 (11)
S＝S1+S2+S3+S4+S5+S6 (11)
Wherein S 1、S2、S3、S4、S5、S6 is the area of the surrounding city of the rigidity curve and the test rigidity curve calculated in X, Y, Z three translation and three rotation directions respectively, and S is the sum of all areas; a n、cn、en、gn、in、kn is the reaction force of the displacement of the stiffness curve at the nth point calculated in X, Y, Z three translation and three rotation directions respectively; b n、dn、fn、hn、jn、ln are respectively X, Y, Z counter forces of displacement of the stiffness curve at the nth point in three translation and three rotation directions.
Further, in the step D, one stress field result is a relative displacement plus 50mm stress result, and the other stress field result is a relative displacement minus 50mm stress result;
further, in the step E, the cycle number is 30 ten thousand, the evaluation index is 1, the damage value is larger than 1, the durability requirement is not met, and the damage value is smaller than 1.
Compared with the prior art, the invention has the beneficial effects that:
According to the invention, the contact constraint among the rubber bushing, the stabilizer bar and the cover plate is established, so that the load transmission is ensured to be accurate; setting a simulation assembly process through bolt pre-tightening and contact constraint initial interference; the real rubber bushing model is used for replacing a virtual rubber bushing, so that the problem that the rigidity test data of the bushing does not meet the finite element simulation requirement is solved; through the measures, three difficult problems in the simulation process of the stabilizer bar system are solved, and the simulation precision is improved.
Drawings
In order to more clearly illustrate the technical solutions of the embodiments of the present invention, the drawings that are needed in the embodiments will be briefly described below, it being understood that the following drawings only illustrate some embodiments of the present invention and therefore should not be considered as limiting the scope, and other related drawings may be obtained according to these drawings without inventive effort for a person skilled in the art.
FIG. 1 is a flow chart of steps of a method of the present invention for improving fatigue durability simulation accuracy of a stabilizer bar system;
FIG. 2 is a schematic structural view of a stabilizer bar system;
FIG. 3 is a flow chart for parameter optimization.
Detailed Description
The invention is further illustrated by the following examples:
to achieve the object of the invention, a stabilizer bar system for a passenger car is taken as an example, and the process of the invention is described in detail. The invention adopts hypermesh as pretreatment, abaqus as solver, isight as parameter optimization process building platform, femfat software is used for calculating durable damage value, and the flow chart is shown in figure 1.
1. Finite element modeling
As shown in fig. 2, the stabilizer bar system part is composed of four types of stabilizer bars, a plurality of rubber bushings, a plurality of stabilizer bar cover plates and a plurality of bolts. In the embodiment, the number of the rubber bushings and the stabilizer bar cover plates is 2, and the number of the bolts is 4. Finite element model building is divided into 4 steps:
1. Grid painting
The grids are drawn and divided sequentially for the four parts, and as the rubber elastic modulus is lower, larger grids deform in the deformation process, the rubber bushing is drawn and divided into hexahedral grids, and meanwhile, the length-width ratio is ensured to be more than or equal to 5, and the grids of other parts are drawn and divided into second-order tetrahedral grids.
2. Establishing material properties
The rubber body structure material adopts a two-term Mooney-Rivlin model, namely a formula (1), the initial values of the rubber constitutive parameters C 10 and C 01 are calculated by inputting the rubber hardness through an empirical formula (2) -a formula (4), and the material density is given. In this example, the rubber hardness H r =60, and the material density was 1.3e-9ton/mm 3 according to the formula C 10＝0.48,C01 =0.12. Other structures are steel materials, and the elastic modulus E=210000 MPa and the Poisson ratio v=0.3.
Two-term Mooney-Rivlin model
W＝C10(I1-3)+C01(I2-3) (1)
Wherein, C 10,C01 is a material constant;
the model not only uses a small deformation range than classical, but is also applicable to a larger range of deformation. Due to simplicity and practicality. The model is widely used in FEA.
The empirical formula of the constitutive parameters of the material is as follows:
lgE0＝0.0184Hr-0.4575 (2)
G＝E0/3＝2(C10+C01) (3)
C01/C10＝0.05 (4)
Wherein E 0 is the elastic modulus of the rubber, H r is the Rockwell hardness of the rubber, G is the shear modulus, and C 10 and C 01 are the determined material constants.
3. Imparting cell attributes
All structures use physical units.
4. Model assembly
Connecting bolt holes at two ends of the stabilizer bar by using a rigid unit, wherein a main point is a center point of the bolt holes; the bush and the stabilizer bar are used for tie constraint, and the bush and the cover plate, the bolt and the cover plate and the auxiliary frame local structure (used for fixing the stabilizer bar) are used for establishing contact constraint, wherein the bush and the cover plate structure have initial penetration, and the penetration is interference.
5. Establishing constraint and load conditions
The auxiliary frame partial structure is restrained, the load working condition comprises two load steps, and the first step is the assembly working condition which comprises bolt pre-tightening and interference fit of a bushing and a cover plate; and the second step is an external load working condition, wherein the external load is the relative displacement quantity of the two ends of the stabilizer bar, the relative displacement quantity is related to the left-right relative runout quantity of the wheel center, and the two ends of the stabilizer bar are converted through the lever ratio.
2. Calibration of rubber material constitutive parameters
The initial rubber constitutive parameters are calculated through rubber hardness, and have certain differences with the actual rubber constitutive parameters. The actual rubber constitutive parameters need to be calibrated through parameter optimization means.
Calculating the area S enclosed by the stiffness curve and the test stiffness curve by taking C 10 and C 01 as variables, and calibrating the rubber constitutive parameters by taking the minimum value of the formula (5) -the formula (11) as a target. Before the test, the displacement output interval of the calculated stiffness curve and the test stiffness curve is ensured to be consistent through the result output and data processing means, the simulation working condition is the same as that of the rubber bushing rack stiffness test, the bushing outer tube and the inner tube are restrained to be loaded, and the load and displacement curve is obtained. And setting up an optimization flow by isight software, and calling an abaqus solver to solve, wherein the optimization flow is shown in fig. 3. Specifically, a test stiffness curve is prepared, then a stiffness calculation finite element model is prepared, parameter optimization flow iterative computation is built through parameterization C 10 and C 01 and target setting S, parameter optimization iterative computation is executed, and finally parameters C 10 and C 01 are determined.
In the present embodiment of the present invention,
S1＝(a1-b1)2+(a2-b2)2+…(an-bn)2 (5)
S2＝(c1-d1)2+(c2-d2)2+…(cn-dn)2 (6)
S3＝(e1-f1)2+(e2-f2)2+…(en-fn)2 (7)
S4＝(g1-h1)2+(g2-h2)2+…(gn-hn)2 (8)
S5＝(i1-j1)2+(i2-j2)2+…(in-jn)2 (9)
S6＝(k1-l1)2+(k2-l2)2+…(kn-ln)2 (10)
S＝S1+S2+S3+S4+S5+S6 (11)
Wherein S 1、S2、S3、S4、S5、S6 is the area of the surrounding city of the rigidity curve and the test rigidity curve calculated in X, Y, Z three translation and three rotation directions respectively, and S is the sum of all areas; a n、cn、en、gn、in、kn is the reaction force of the displacement of the stiffness curve at the nth point calculated in X, Y, Z three translation and three rotation directions respectively; b n、dn、fn、hn、jn、ln are respectively X, Y, Z counter forces of displacement of the stiffness curve at the nth point in three translation and three rotation directions.
Through the parameter optimization flow, the rubber constitutive parameter is calibrated to be C 10＝0.52,,C01 =0.13.
3. Assembly process simulation
The assembly process mainly comprises the steps of applying bolt pre-tightening load, pressing and connecting the cover plate to a sub-frame local structure, and simultaneously pressing and assembling the cover plate and the rubber bushing together, wherein the initial interference of the bushing and the cover plate is also taken into consideration in the process.
4. Durable condition calculation
Based on the pre-tightening working condition in the previous step, the simulation calculation of the external load fatigue working condition is completed, the local structure of the auxiliary frame is restrained, the relative displacement of plus or minus 50mm is added at the two ends of the stabilizer bar, two stress field results are obtained by using abaqus to solve, one is the stress result of the relative displacement plus 50mm, and the other is the stress result of the relative displacement minus 50 mm.
5. Durable damage value calculation
And step four, the results of the two stress fields calculated in the step four are imported into femfat software or other endurance calculation software to finish the calculation of the constant-amplitude fatigue endurance damage value between the two stress fields, the cycle number is 30 ten thousand times, the evaluation index is 1, the damage value result is larger than 1, the endurance performance requirement is not met, and the damage value result is smaller than 1. In the embodiment, the damage value of the stabilizer bar is 0.96, the failure positions of the stabilizer bar in 34.5 ten thousand times are the same, the simulation precision is over 90 percent, the damage value of the stabilizer bar calculated by the original method is 0.42, and the precision is greatly improved.
Through the five steps, the invention process of improving the fatigue endurance simulation precision of the stabilizer bar is completed, and compared with the traditional simulation method, the simulation precision is greatly improved.
Note that the above is only a preferred embodiment of the present invention and the technical principle applied. It will be understood by those skilled in the art that the present invention is not limited to the particular embodiments described herein, but is capable of various obvious changes, rearrangements and substitutions as will now become apparent to those skilled in the art without departing from the scope of the invention. Therefore, while the invention has been described in connection with the above embodiments, the invention is not limited to the embodiments, but may be embodied in many other equivalent forms without departing from the spirit or scope of the invention, which is set forth in the following claims.
Claims (8)
1. The method for improving the fatigue durability simulation precision of the stabilizer bar system is characterized by comprising the following steps of:
A. finite element modeling
A1, meshing stabilizer bar system parts in sequence, wherein the stabilizer bar system parts consist of stabilizer bars, a plurality of rubber bushings, a plurality of stabilizer bar cover plates and a plurality of bolts;
a2, establishing material properties;
The method comprises the following steps: the rubber body structure material adopts a two-term Mooney-Rivlin model, namely a formula (1), initial values of rubber constitutive parameters C 10 and C 01 are calculated through the input of rubber hardness of a material constitutive parameter empirical formula (2) -a formula (4), and the density of the material is endowed;
W＝C10(I1-3)+C01(I2-3) (1)
lgE0＝0.0184Hr-0.4575 (2)
G＝E0/3＝2(C10+C01) (3)
C01/C10＝0.05 (4)
wherein C 10,C01 is a material constant; e 0 is the elastic modulus of the rubber, H r is the Rockwell hardness of the rubber, and G is the shear modulus;
A3, endowing the unit with attributes, wherein all structures use entity units;
A4, assembling a model;
A5, establishing constraint and load working conditions;
B. The rubber material constitutive parameters are calibrated by calculating initial constitutive parameters of rubber through rubber hardness, and then the actual constitutive parameters of the rubber are calibrated through parameter optimization means;
b1, ensuring that the displacement output interval of a calculated stiffness curve is consistent with that of a test stiffness curve through result output and data processing, restricting the loading of an outer tube and an inner tube of a bushing under the same simulation working condition as that of a rubber bushing rack stiffness test, and solving a load and displacement curve;
B2, calculating the minimum value of the area S surrounded by the stiffness curve and the test stiffness curve by taking C 10 and C 01 as variables as targets, and calibrating the rubber constitutive parameters;
B3, parameter optimization: firstly preparing a test stiffness curve, then preparing a stiffness calculation finite element model, constructing parameter optimization flow iterative computation through parameterization C 10 and C 01 and target setting S, executing parameter optimization iterative computation, and finally determining parameters C 10 and C 01;
C. Assembly process simulation
Applying a bolt pre-tightening load, pressing and connecting the cover plate to a local structure of the auxiliary frame, and simultaneously pressing and assembling the cover plate and the rubber bushing together;
D. Durable condition calculation
C, based on the step, performing simulation calculation of external load fatigue working conditions, restraining local structures of auxiliary frames, adding certain relative displacement at two ends of a stabilizer bar, and solving to obtain two stress field results;
E. Durable damage value calculation
And D, importing the results of the two stress fields calculated in the step D, and finishing calculation of the constant-amplitude fatigue endurance damage value between the two stress fields.
2. The method for improving fatigue durability simulation accuracy of the stabilizer bar system according to claim 1, wherein the method comprises the following steps: and A1, the number of the rubber bushings and the stabilizer bar cover plates is 2, and the number of bolts is 4.
3. The method for improving fatigue durability simulation accuracy of the stabilizer bar system according to claim 1, wherein the method comprises the following steps: and A1, dividing the rubber bushing picture into hexahedral grids, and dividing the rubber bushing picture into second-order tetrahedral grids, wherein the length-width ratio of the rubber bushing picture is more than or equal to 5, and the other part grids are ensured to be divided into second-order tetrahedral grids.
4. The method for improving fatigue durability simulation precision of a stabilizer bar system according to claim 1, wherein step A4 specifically comprises: connecting bolt holes at two ends of the stabilizer bar by using a rigid unit, wherein a main point is a center point of the bolt holes; the bush and the stabilizer bar are used for tie constraint, and the bush and the cover plate, the bolt and the cover plate and the auxiliary frame are used for local structure contact constraint, wherein the bush and the cover plate structure have initial penetration, and the penetration is interference.
5. The method for improving fatigue durability simulation precision of a stabilizer bar system according to claim 1, wherein step A5 specifically comprises: the auxiliary frame partial structure is restrained, the load working condition comprises two load steps, and the first step is the assembly working condition which comprises bolt pre-tightening and interference fit of a bushing and a cover plate; and the second step is an external load working condition, wherein the external load is the relative displacement quantity of the two ends of the stabilizer bar, the relative displacement quantity is related to the left-right relative runout quantity of the wheel center, and the two ends of the stabilizer bar are converted through the lever ratio.
6. The method for improving fatigue durability simulation precision of a stabilizer bar system according to claim 1, wherein the specific steps of calculating the area S enclosed by the stiffness curve and the test stiffness curve are as follows:
S1＝(a1-b1)2+(a2-b2)2+…(an-bn)2 (5)
S2＝(c1-d1)2+(c2-d2)2+…(cn-dn)2 (6)
S3＝(e1-f1)2+(e2-f2)2+…(en-fn)2 (7)
S4＝(g1-h1)2+(g2-h2)2+…(gn-hn)2 (8)
S5＝(i1-j1)2+(i2-j2)2+…(in-jn)2 (9)
S6＝(k1-l1)2+(k2-l2)2+…(kn-ln)2 (10)
S＝S1+S2+S3+S4+S5+S6 (11)
Wherein S 1、S2、S3、S4、S5、S6 is the area of the surrounding city of the rigidity curve and the test rigidity curve calculated in X, Y, Z three translation and three rotation directions respectively, and S is the sum of all areas; a n、cn、en、gn、in、kn is the reaction force of the displacement of the stiffness curve at the nth point calculated in X, Y, Z three translation and three rotation directions respectively; b n、dn、fn、hn、jn、ln are respectively X, Y, Z counter forces of displacement of the stiffness curve at the nth point in three translation and three rotation directions.
7. The method for improving fatigue durability simulation accuracy of the stabilizer bar system according to claim 1, wherein the method comprises the following steps: and D, the two stress field results are a relative displacement plus 50mm stress result and a relative displacement minus 50mm stress result.
8. The method for improving fatigue durability simulation accuracy of the stabilizer bar system according to claim 1, wherein the method comprises the following steps: and E, the cycle times are 30 ten thousand times, the evaluation index is 1, the damage value result is larger than 1, the durability requirement is not met, and the damage value result is smaller than 1.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210424742.1A CN114970237B (en) | 2022-04-21 | Method for improving fatigue durability simulation precision of stabilizer bar system |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210424742.1A CN114970237B (en) | 2022-04-21 | Method for improving fatigue durability simulation precision of stabilizer bar system |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN114970237A CN114970237A (en) | 2022-08-30 |
| CN114970237B true CN114970237B (en) | 2024-07-05 |
Family
ID=
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105699205A (en) * | 2016-03-24 | 2016-06-22 | 安徽机电职业技术学院 | Stabilization rod rigidity testing device |
| CN110309598A (en) * | 2019-07-03 | 2019-10-08 | 中汽研（天津）汽车工程研究院有限公司 | A kind of creation method of Virtual Test Field |
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105699205A (en) * | 2016-03-24 | 2016-06-22 | 安徽机电职业技术学院 | Stabilization rod rigidity testing device |
| CN110309598A (en) * | 2019-07-03 | 2019-10-08 | 中汽研（天津）汽车工程研究院有限公司 | A kind of creation method of Virtual Test Field |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN105095542B (en) | A kind of automotive suspension key structure piece optimization design method | |
| CN111950170B (en) | Method for obtaining high-precision Macpherson front suspension knuckle rack test load | |
| CN116090111B (en) | Automobile leaf spring fatigue life prediction method based on deep learning model | |
| CN113343374B (en) | Automobile plate spring fatigue testing method | |
| CN110414052A (en) | A kind of vehicle structure fatigue life calculation method composed based on digital sample car and virtual road | |
| CN105844039A (en) | Method for evaluating local rigidity contribution degrees of body structure | |
| CN115169167A (en) | Method and system for optimizing and matching motion stroke parameters of automobile plate spring | |
| CN108875188B (en) | Method and device for optimizing a body joint of a motor vehicle | |
| CN114970237B (en) | Method for improving fatigue durability simulation precision of stabilizer bar system | |
| WO2024001902A1 (en) | Method and system for calculating maximum stress distribution of parts of vehicle chassis and storage medium | |
| CN113850024A (en) | Machine learning-based reinforced concrete member anti-collision performance prediction method | |
| CN111859721A (en) | Automobile lower swing arm optimization design method | |
| CN110955933B (en) | Mechanical structure fuzzy fatigue reliability calculation method based on response surface method | |
| CN109033726B (en) | Finite element method-based chassis suspension system strength and rack durability analysis method | |
| CN114970237A (en) | Method for improving fatigue endurance simulation precision of stabilizer bar system | |
| CN112685831B (en) | Optimization method for suspension and auxiliary frame part structures | |
| CN114880904A (en) | Large-scale finite element analysis method for rubber bushing | |
| CN115081116A (en) | Method and system for calculating maximum stress of vehicle chassis part and storage medium | |
| CN112069700A (en) | Method and device for calculating strength of torsion-resistant pull rod of power assembly and computer equipment | |
| CN113190915A (en) | Vehicle strength evaluation method and system | |
| CN116992714A (en) | Method for simulating and analyzing strength and durability of auxiliary frame of passenger car | |
| CN114647897B (en) | Cab structure optimization method based on hierarchical entropy-gray correlation method | |
| CN112651089B (en) | Bending shear metal component hybrid reinforcement restoring force model and parameter determination method | |
| CN116839783B (en) | Method for measuring stress value and deformation of automobile leaf spring based on machine learning | |
| CN117034694A (en) | Vibration response calculation method based on whole vehicle simulation model |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant |