CN108446528A - Front suspension optimum design method, device and computer readable storage medium - Google Patents
Front suspension optimum design method, device and computer readable storage medium Download PDFInfo
- Publication number
- CN108446528A CN108446528A CN201810560437.9A CN201810560437A CN108446528A CN 108446528 A CN108446528 A CN 108446528A CN 201810560437 A CN201810560437 A CN 201810560437A CN 108446528 A CN108446528 A CN 108446528A
- Authority
- CN
- China
- Prior art keywords
- front suspension
- coordinates
- optimum design
- optimization
- surface model
- 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.)
- Pending
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/15—Vehicle, aircraft or watercraft design
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Geometry (AREA)
- General Physics & Mathematics (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Computer Hardware Design (AREA)
- Automation & Control Theory (AREA)
- Aviation & Aerospace Engineering (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
The invention discloses a kind of front suspension optimum design methods, include the following steps：Obtain the hard spot coordinate of Simulation Calculation；Sensitivity analysis is carried out to the coordinate parameters of the hard spot coordinate, and obtains coordinates of targets parameter；The second-order response surface model of the Simulation Calculation is built according to the coordinates of targets parameter；According to the second-order response surface model foundation multi-goal optimizing function；The Optimum Design Results of optimization aim are obtained by the multi-goal optimizing function.Invention additionally discloses a kind of front suspension optimization design device and computer readable storage mediums.The present invention is designed a model by second-order response surface model instead of simulation optimization, solves the technical problem that Optimal design and calculation amount is big, and the time-consuming period is grown.
Description
Technical field
The present invention relates to automobile technical field more particularly to front suspension optimum design method, device and computer-readable deposit
Storage media.
Background technology
With the development of society, automobile is more and more universal, therefore the design technology of automobile is also increasingly perfect therewith.Its
Middle change suspension frame structure hard spot reduces variable quantity of wheel alignment parameter during wheel hop, is effectively improved as one kind
The method of suspension property is studied extensively.
The prior art generally by choosing and improving the position coordinates of hard spot, carries out l-G simulation test, and then judge current shape
Whether the parameter value under condition meets the requirements, to realize the optimization design to suspension.Due to above-mentioned position coordinates variation range very
Greatly, cause this design process operand very big and the time-consuming period is long, while effect that can not be to positional parameter and significance level
Ad hoc analysis is carried out, causes the specific aim of optimum results poor.
The above is only used to facilitate the understanding of the technical scheme, and is not represented and is recognized that the above is existing skill
Art.
Invention content
The main purpose of the present invention is to provide a kind of front suspension optimum design method, device and computer-readable storage mediums
Matter, it is intended to realize the calculation amount for reducing front suspension optimization design and carry out the mesh of specific aim analysis to the significance level of positional parameter
's.
To achieve the above object, the present invention provides a kind of front suspension optimum design method, the front suspension optimization design side
Method includes the following steps：
Obtain the hard spot coordinate of Simulation Calculation；
Sensitivity analysis is carried out to the coordinate parameters of the hard spot coordinate, and obtains coordinates of targets parameter；
The second-order response surface model of the Simulation Calculation is built according to the coordinates of targets parameter；
According to the second-order response surface model foundation multi-goal optimizing function；
The Optimum Design Results of optimization aim are obtained by the multi-goal optimizing function.
Preferably, the coordinate parameters to the hard spot coordinate carry out sensitivity analysis, and obtain coordinates of targets parameter
The step of include：
Determine the variation range of the coordinate parameters of the hard spot coordinate；
According to the variation range, the coordinate parameters are emulated by dynamics analysis software；
The coordinate parameters for meeting sensitivity condition are chosen according to simulation result, as coordinates of targets parameter.
Preferably, described the step of meeting the coordinate parameters of sensitivity condition according to simulation result selection, includes：
The influence value to the Simulation Calculation of the coordinate parameters is determined according to the simulation result；
Judge whether the coordinate parameters meet the sensitivity condition according to the influence value, and chooses and meet the spirit
The above-mentioned coordinate parameters of sensitivity condition are as coordinates of targets parameter.
Preferably, the second-order response surface model that the Simulation Calculation is built according to the coordinates of targets parameter
Step includes：
Central Composite is carried out to the coordinates of targets parameter to test to obtain test result；
The second-order response surface model of the optimization aim is determined according to the test result.
Preferably, the front suspension optimum design method further includes：
Obtain the regression coefficient of the second-order response surface model；
The multiple correlation coefficient for calculating the regression coefficient, to determine the second-order response face mould according to the multiple correlation coefficient
The accuracy of type；
When the accuracy is unsatisfactory for precise requirements, the second-order response surface model is corrected.
Preferably, described to include according to the step of second-order response surface model foundation multi-goal optimizing function：
The correspondence weighted factor of the optimization aim is calculated by the second-order response surface model；
It is weighted combination according to the weighted factor and the optimization aim, to establish multi-goal optimizing function.
Preferably, the step of Optimum Design Results that optimization aim is obtained by the multi-goal optimizing function wrap
It includes：
Obtain the variation range of the coordinates of targets parameter；
The coordinates of targets parameter in the variation range is optimized by multi-objective optimization algorithm；
The coordinates of targets parameter after the multi-objective optimization algorithm is optimized inputs the multi-goal optimizing function, with
Obtain the Optimum Design Results of optimization aim.
In addition, to achieve the above object, the present invention also provides a kind of front suspension optimization design devices, which is characterized in that institute
Stating front suspension optimization design device includes：It memory, processor and is stored on the memory and can be on the processor
The Optimized Program of operation, the Optimized Program realize front suspension optimization as described above when being executed by the processor
The step of design method.
In addition, to achieve the above object, the present invention also provides a kind of computer readable storage mediums, which is characterized in that institute
It states and is stored with Optimized Program on computer readable storage medium, realized such as when the Optimized Program is executed by processor
Above the step of front suspension optimum design method.
A kind of front suspension optimum design method, device and the computer readable storage medium that the embodiment of the present invention proposes, institute
Front suspension optimum design method is stated, the hard spot coordinate of Simulation Calculation is first obtained, then the coordinate of the hard spot coordinate is joined
Number carries out sensitivity analysis, and obtains coordinates of targets parameter so that can build the emulation according to the coordinates of targets parameter
The second-order response surface model of computation model, and then according to the second-order response surface model foundation multi-goal optimizing function, finally lead to
Cross the Optimum Design Results that the multi-goal optimizing function obtains optimization aim.Since the present invention passes through second-order response surface model generation
It is optimized for actual emulation computation model, and optimized variable chooses the highest optimized variable of susceptibility, so that contracting
Subtract the calculation amount in process of optimization, and makes optimum results more targeted.
Description of the drawings
Fig. 1 is the terminal structure schematic diagram for the hardware running environment that the embodiment of the present invention is related to；
Fig. 2 is the flow diagram of front suspension optimum design method first embodiment of the present invention；
Fig. 3 is the flow diagram of front suspension optimum design method second embodiment of the present invention；
Fig. 4 is the flow diagram of front suspension optimum design method 3rd embodiment of the present invention；
Fig. 5 is the flow diagram of front suspension optimum design method fourth embodiment of the present invention；
Fig. 6 is the flow diagram of the 5th embodiment of front suspension optimum design method of the present invention；
Fig. 7 is the flow diagram of front suspension optimum design method sixth embodiment of the present invention；
Fig. 8 is the flow diagram of the 7th embodiment of front suspension optimum design method of the present invention.
The embodiments will be further described with reference to the accompanying drawings for the realization, the function and the advantages of the object of the present invention.
Specific implementation mode
It should be appreciated that the specific embodiments described herein are merely illustrative of the present invention, it is not intended to limit the present invention.
The primary solutions of the embodiment of the present invention are：
Obtain the hard spot coordinate of Simulation Calculation；
Sensitivity analysis is carried out to the coordinate parameters of the hard spot coordinate, and obtains coordinates of targets parameter；
The second-order response surface model of the Simulation Calculation is built according to the coordinates of targets parameter；
According to the second-order response surface model foundation multi-goal optimizing function；
The Optimum Design Results of optimization aim are obtained by the multi-goal optimizing function.
Due in the prior art, generally by choosing and improving the position coordinates of hard spot, carrying out l-G simulation test, and then sentence
Whether the parameter value under disconnected the present situation meets the requirements, to realize the optimization design to suspension.Due to the change of above-mentioned position coordinates
It is very big to change range, causes this design process operand very big and the time-consuming period is long, at the same can not to the effect of positional parameter and
Significance level carries out ad hoc analysis, causes the specific aim of optimum results poor.
A kind of front suspension optimum design method, device and the computer readable storage medium that the embodiment of the present invention proposes, institute
Front suspension optimum design method is stated, the hard spot coordinate of Simulation Calculation is first obtained, then the coordinate of the hard spot coordinate is joined
Number carries out sensitivity analysis, and obtains coordinates of targets parameter so that can build the emulation according to the coordinates of targets parameter
The second-order response surface model of computation model, and then according to the second-order response surface model foundation multi-goal optimizing function, finally lead to
Cross the Optimum Design Results that the multi-goal optimizing function obtains optimization aim.Since the present invention passes through second-order response surface model generation
It is optimized for actual emulation computation model, and optimized variable chooses the highest optimized variable of susceptibility, so that contracting
Subtract the calculation amount in process of optimization, and makes optimum results more targeted.
As shown in Figure 1, the terminal structure schematic diagram for the hardware running environment that Fig. 1, which is the embodiment of the present invention, to be related to.
Terminal of the embodiment of the present invention can be PC, can also be that the terminals such as pocket computer, tablet computer or server are set
It is standby.
As shown in Figure 1, the terminal may include：Processor 1001, such as CPU, network interface 1004, user interface
1003, memory 1005, communication bus 1002.Wherein, communication bus 1002 is for realizing the connection communication between these components.
User interface 1003 may include display screen (Display), input unit such as keyboard (Keyboard), mouse etc., can be selected
Family interface 1003 can also include standard wireline interface and wireless interface.Network interface 1004 may include optionally standard
Wireline interface, wireless interface (such as WI-FI interfaces).Memory 1005 can be high-speed RAM memory, can also be stable deposit
Reservoir (non-volatile memory), such as magnetic disk storage.Memory 1005 optionally can also be independently of aforementioned place
Manage the storage device of device 1001.
It will be understood by those skilled in the art that the restriction of the not structure paired terminal of terminal structure shown in Fig. 1, can wrap
It includes than illustrating more or fewer components, either combines certain components or different components arrangement.
As shown in Figure 1, as may include that operating system, network are logical in a kind of memory 1005 of computer storage media
Believe module, Subscriber Interface Module SIM and Optimized Program.
In terminal shown in Fig. 1, network interface 1004 is mainly used for connecting background server, is carried out with background server
Data communicate；User interface 1003 is mainly used for connecting client (user terminal), with client into row data communication；And processor
1001 can be used for calling the Optimized Program stored in memory 1005, and execute following operation：
Obtain the hard spot coordinate of Simulation Calculation；
Sensitivity analysis is carried out to the coordinate parameters of the hard spot coordinate, and obtains coordinates of targets parameter；
The second-order response surface model of the Simulation Calculation is built according to the coordinates of targets parameter；
According to the second-order response surface model foundation multi-goal optimizing function；
The Optimum Design Results of optimization aim are obtained by the multi-goal optimizing function.
Further, processor 1001 can call the Optimized Program stored in memory 1005, also execute following
Operation：
Determine the variation range of the coordinate parameters of the hard spot coordinate；
According to the variation range, the coordinate parameters are emulated by dynamics analysis software；
The coordinate parameters for meeting sensitivity condition are chosen according to simulation result, as coordinates of targets parameter.
Further, processor 1001 can call the Optimized Program stored in memory 1005, also execute following
Operation：
The influence value to the Simulation Calculation of the coordinate parameters is determined according to the simulation result；
Judge whether the coordinate parameters meet the sensitivity condition according to the influence value, and chooses and meet the spirit
The above-mentioned coordinate parameters of sensitivity condition are as coordinates of targets parameter.
Further, processor 1001 can call the Optimized Program stored in memory 1005, also execute following
Operation：
Central Composite is carried out to the coordinates of targets parameter to test to obtain test result；
The second-order response surface model of the optimization aim is determined according to the test result.
Further, processor 1001 can call the Optimized Program stored in memory 1005, also execute following
Operation：
Obtain the regression coefficient of the second-order response surface model；
The multiple correlation coefficient for calculating the regression coefficient, to determine the second-order response face mould according to the multiple correlation coefficient
The accuracy of type；
When the accuracy is unsatisfactory for precise requirements, the second-order response surface model is corrected.
Further, processor 1001 can call the Optimized Program stored in memory 1005, also execute following
Operation：
The correspondence weighted factor of the optimization aim is calculated by the second-order response surface model；
It is weighted combination according to the weighted factor and the optimization aim, to establish multi-goal optimizing function.
Further, processor 1001 can call the Optimized Program stored in memory 1005, also execute following
Operation：
Obtain the variation range of the coordinates of targets parameter；
The coordinates of targets parameter in the variation range is optimized by multi-objective optimization algorithm；
The coordinates of targets parameter after the multi-objective optimization algorithm is optimized inputs the multi-goal optimizing function, with
Obtain the Optimum Design Results of optimization aim.
With reference to Fig. 2, front suspension optimum design method first embodiment of the present invention, the front suspension optimum design method packet
It includes：
Step S10 obtains the hard spot coordinate of Simulation Calculation；
In the present embodiment, the Simulation Calculation of front suspension is first established, and then is obtained according to the Simulation Calculation
Front suspension hard spot coordinate.
Specifically, for example, when the front suspension to certain vehicle optimizes practical, for more body power may be used
It learns software ADAMS/Car and establishes dynamics simulation model.In above-mentioned Simulation Calculation, except elastic element and rubber member
Outside part, all parts can be considered as rigid body；All connections between parts can be reduced to rigid hinge.
Further, simplified Suspension Model may include：Lower control arm, knuckle assembly, damper and the horizontal drawing of steering
Bar.Other mass property parameters, the mechanics parameters of spring, bushing and damper are measured by testing.
Front suspension hard spot parameter is more, but suspension main pin axis and suspension hanging point and the space geometry and fortune of guiding mechanism
Correlation in close relations is moved, is to determine the core element of wheel alignment parameter.By adjusting under main pin axis transversal swinging arm and steering
The length of hinge joint and track rod is saved to adjust wheel alignment parameter, therefore some insensitive hard spots ginsengs can be saved
Number, so choose in track rod under tie point, damper tie point before installation point, lower swing arm, lower swing arm rear tie points, under
Hard spot of the 5 crucial hard spots such as outer tie point of swing arm as computation model, and obtain the initial coordinate (x, y, z) of each hard spot.
In addition, when establishing Simulation Calculation, first front suspension simplify and it is assumed that its content may include：It is outstanding
Rod piece in frame does not deform with the bounce of wheel, it is assumed that rod piece is rigid body；The deformation for not considering wheel, wheel letter
Turn to rigid body；Internal clearance between parts is disregarded, and connection type is all reduced to hinge connection；All parts are all recognized
To be rigid body, each kinematic pair is rigid connection, ignores the frictional force between each kinematic pair；Vehicle body is motionless with respect to ground, only
Study the characteristic of suspension.
Step S20 carries out sensitivity analysis to the coordinate parameters of the hard spot coordinate, and obtains coordinates of targets parameter；
In the present embodiment, sensitivity refers to that design parameter changes the qualitative assessment influenced on design performance index, sensitive
Degree analysis be research and analyse system (or model) state or output variation systematic parameter or ambient conditions are changed it is quick
The method of sense degree.Sensitivity analysis is frequently utilized that in optimal method to study initial data inaccuracy or change
When optimal solution stability.It can determine which parameter has large effect to system or model by sensitivity analysis.This reality
Input variable of the example using the coordinate parameters of hard spot coordinate as sensitivity analysis is applied, then by sensitivity analysis, chooses and meets
The input variable of sensitivity condition is as coordinates of targets parameter.
Specifically, the experimental design of use carries out sensitivity analysis, and experimental design discussion is multiple variables while occurring
When variation, influence of each design variable for model machine performance.The present embodiment can carry out sensitivity analysis according to following steps：
Test objective is specified, the purpose of experimental design is found out in numerous input variables is affected to model machine performance
Parameter, and determine object function as evaluation goal；
Testing program is formulated, limits a range to each design factor, test factor can only become within this range
Change, all combined situations are then arranged in experimental design matrix；
It is tested, is emulated in corresponding dynamics analysis software, system can be directed to each variable-value automatically
Situation is all once emulated, and recording simulation results.
Test result is analyzed by simulation software, and impact factor contribution degree is obtained according to test result, and then selects tribute
Degree of offering meets susceptibility condition, as coordinates of targets parameter, wherein it refers to its influence value (impact factor to meet susceptibility condition
Contribution degree) it is more than default influence value (can be arranged by User Defined).
For example, choosing to tie point P in track rod1, installation point P under damper2, tie point P before lower swing arm3, the bottom
Arm rear tie points P4, the outer tie point P of lower swing arm5Coordinate parameters Deng 5 hard spots as the input variable sensitively analyzed, i.e., with
P1、P2、P3、P4And P5Change of the coordinate parameters on three directions of x, y, z of space coordinates be turned to the defeated of sensitivity analysis
Enter variable.
Further, the sensitivity analysis carried out in ADAMS/Insight, by P1、P2、P3、P4And P5Deng 5 hard spots
Totally 15 coordinate parameters obtain highest 6 coordinate parameters of sensitivity, i.e., installation point y-coordinate x under damper as input quantity1、
Tie point z coordinate x in drag link2, tie point y-coordinate x before lower swing arm3, tie point z coordinate x before lower swing arm4, connect outside lower swing arm
Point y-coordinate x5, the outer tie point z coordinate x of lower swing arm6, using these coordinates as coordinates of targets parameter.
Step S30 builds the second-order response surface model of the Simulation Calculation according to the coordinates of targets parameter；
In the present embodiment, center experimental design is first carried out according to coordinates of targets parameter, is then designed according to Central Composite
As a result second-order response surface model is established.
Specifically, it establishes response surface approximate model to need to consider a certain number of experimental design data, using experimental design
When need to make the entire design space of experimental design sample instinct reflection.Central Composite experimental design is used to establish second-order response surface model,
It is capable of providing the ability of abundant information and fitting second order relationship about experiment variable and test error.
For example, Central Composite experiment can be carried out according to Central Composite experimental design scheme table as follows：
Central Composite experimental design scheme table：
Further, second-order response surface model is the multinomial with clear expression-form by approximation one, right
Response target is influenced to carry out modeling analysis by multiple variables, according to the second-order response for the generation that Central Composite test result obtains
Surface model is as follows：
Wherein, xi、xjFor independent variable, β0、βi、βii、βijFor regression coefficient, the number N=(n+2) (n+1) of regression coefficient/
2；ε is approximate error, in the case where meeting requirement of engineering precision it is believed that ε=0, n are cumulative peak value.
Step S40, according to the second-order response surface model foundation multi-goal optimizing function；
Specifically, toe-in angle sets to coordinate front-wheel camber, to mitigate tire wear；Kingpin inclination and
Castor is then to generate aligning torque and reduce steering force when turning to.According to the above analysis, before the suspension
Wheel positional parameter is divided into two groups.Wherein, toe-in of front wheel angle and front-wheel camber are one group, and kingpin inclination and castor are
One group.Every group of parameter is set as an object function f (x) by the way that set of weights is legal, form is as follows：
f1(x)=w1y1(x)+w2y2(x)
f2(x)=w3y3(x)+w4y4(x)
Wherein, wiFor each subhead scalar functions yi(x) weighted factor, value are decided by the order of magnitude of the objectives and important
Degree.Influence of the point counting target to vehicle control stability and tire wear for convenience is established unified using direct weighting method
Object function, weighted factor wiChoosing method is as follows：
If known certain subhead scalar functions yi(x) mobility scale is：
αi≤yi(x)≤βi, i=1,2,3,4
Then claim：
For the tolerance of the target, the weighted factor that then can use the target is：
This follow the example of is based on requiring each subhead scalar functions in unified target function to tend to reach system on the order of magnitude
One balance.When the numerical value change of a certain object function is wider, the tolerance of target is bigger, and weighted factor is with regard to smaller；And it counts
When value variation is narrower, for the tolerance of target with regard to smaller, weighted factor is bigger, to reach the work for balancing each object function magnitude
With.
It is analyzed according to l-G simulation test, determines the value range of each target variable, the appearance of each target is calculated according still further to formula
Limit and weighted factor can obtain the weighted factor table of the weighted factor of each positional parameter of wheel each target variable as follows:
The weighted factor table of each target variable:
Can obtain multi-goal optimizing function according to the weighted factor of each subhead scalar functions determined above is：
Step S50 obtains the Optimum Design Results of optimization aim by the multi-goal optimizing function.
In the present embodiment, it after being optimized to coordinates of targets parameter according to the multi-goal optimizing function, obtains and achieves
Coordinate parameters after optimization.
In the present embodiment, simplified mathematical model is screened into line sensitivity and is analyzed, and then utilizes Central Composite experimental design,
Second-order response surface model is established, to substitute computation model, finally constructs multiple-objection optimization letter with the method for weighted array
Number, which simplify calculate solution procedure, hence it is evident that shortens and calculates the time, improves optimization design efficiency.
Further, with reference to Fig. 3, front suspension optimum design method second embodiment of the present invention is implemented based on above-mentioned first
Example, the step S20 include：
Step S21 determines the variation range of the coordinate parameters of the hard spot coordinate；
Step S22 emulates the coordinate parameters by dynamics analysis software according to the variation range；
In the present embodiment, each coordinate parameters is given to limit a range, coordinate parameters can only become within this range
Change, all combined situations are then arranged in experimental design matrix.By the experimental design matrix in corresponding dynamic analysis
It is emulated in software, system all can be emulated once for each variable-value situation automatically, and the result of emulation all can
It is recorded.
Step S23 chooses the coordinate parameters for meeting sensitivity condition according to simulation result, joins as coordinates of targets
Number.
In the present embodiment, simulation software voluntarily analyzes test result, and chooses the seat that influence value is more than default influence value
Parameter is marked as coordinates of targets parameter.
It should be noted that the simulation calculation software can be ADAMS/Insight.
In the present embodiment, by sensitivity analysis selection target coordinate parameters, which simplify calculate solution procedure, contracting
The short calculating time.
Further, with reference to Fig. 4, front suspension optimum design method 3rd embodiment of the present invention, based on above-mentioned first to the
Two embodiments, the step S23 include：
Step S231 determines the influence to the Simulation Calculation of the coordinate parameters according to the simulation result
Value；
Step S232 judges whether the coordinate parameters meet the sensitivity condition according to the influence value, and chooses
Meet the above-mentioned coordinate parameters of the sensitivity condition as coordinates of targets parameter.
In the present embodiment, simulation calculation software (can be ADAMS/Insight) is determined according to the simulation experiment result and is sat
Influence value of the parameter to simulation result is marked, the coordinate parameters that influence value is more than default influence value is then selected to join as coordinates of targets
Number.
It should be noted that the default influence value is obtained by user interface, user can be according to input equipments such as keyboards
Influence value is preset in input.
In the present embodiment, according to the influence value of simulation result coordinates computed parameter, influence value is then selected to be more than default shadow
The coordinate parameters of value are rung as coordinates of targets parameter, determine the selection criteria of coordinates of targets parameter in this way.
Further, with reference to Fig. 5, front suspension optimum design method fourth embodiment of the present invention, based on above-mentioned first to the
Three embodiments, the step S30 include：
Step S31 carries out Central Composite to the coordinates of targets parameter and tests to obtain test result；
Step S32 determines the second-order response surface model of the optimization aim according to the test result.
In the present embodiment, composite testing in first being carried out according to the coordinates of targets parameter, is then tried according to Central Composite
The test result tested obtains the information of experiment variable and test error, and obtains the ability of experiment variable fitting second order relationship.
Further, it is closed by software (can be MATLAB), the information of experiment variable and test error and fitting second order
The ability of system determines second-order response surface model.
In the present embodiment, second-order response surface model is established using Central Composite experimental design, to substitute simulation calculation
Model, which further simplifies multifactor, multivariable Simulation Calculations, improve optimization design efficiency.
Further, with reference to Fig. 6, the 5th embodiment of front suspension optimum design method of the present invention, based on above-mentioned first to the
Four embodiments, the front suspension optimum design method further include：
Step S60 obtains the regression coefficient of the second-order response surface model；
Step S70 calculates the multiple correlation coefficient of the regression coefficient, to determine the second order according to the multiple correlation coefficient
The accuracy of response surface model；
Step S80 corrects the second-order response surface model when the accuracy is unsatisfactory for precise requirements.
In the present embodiment, the multiple correlation coefficient of regression coefficient can be calculated according to following formula：
Wherein, SSE is the quadratic sum of response and response estimation value difference；SSY is square of response and the equal value difference of response
With；N is assessment number of test points；yiFor simulation data value；For Response Face Function value；R2It is one between [0,1] to change
Value, value indicate that second-order response surface model is more accurate closer to 1.
Further, (can be R when the accuracy of second-order response surface model is unsatisfactory for precise requirements2Less than 0.9
When, judgement is unsatisfactory for precise requirements), correct the second-order response surface model.
It should be noted that the precise requirements can be set according to different design requirements, the present embodiment is accurate
Degree requires to be not construed as limiting.
In the present embodiment, the multiple correlation coefficient for calculating regression coefficient, to judge the accuracy of second-order response surface model,
When being unsatisfactory for precise requirements, the second-order response surface model is corrected, has ensured the accuracy of optimum results in this way.
Further, with reference to Fig. 7, front suspension optimum design method sixth embodiment of the present invention, based on above-mentioned first to the
Five embodiments, the step S40 include：
Step S41 calculates the correspondence weighted factor of the optimization aim by the second-order response surface model；
Specifically, if known certain subhead scalar functions yi(x) mobility scale is：
αi≤yi(x)≤βi, i=1,2,3,4
Then claim：
For the tolerance of the target, the weighted factor that then can use the target is：
This follow the example of is based on requiring each subhead scalar functions in unified target function to tend to reach system on the order of magnitude
One balance.When the numerical value change of a certain object function is wider, the tolerance of target is bigger, and weighted factor is with regard to smaller；And it counts
When value variation is narrower, for the tolerance of target with regard to smaller, weighted factor is bigger, to reach the work for balancing each object function magnitude
With.
Step S42 is weighted combination according to the weighted factor and the optimization aim, to establish multiple-objection optimization letter
Number.
Specifically, toe-in angle sets to coordinate front-wheel camber, to mitigate tire wear；Kingpin inclination and
Castor is then to generate aligning torque and reduce steering force when turning to.According to the above analysis, before the suspension
Wheel positional parameter is divided into two groups.Wherein, toe-in of front wheel angle and front-wheel camber are one group, and kingpin inclination and castor are
One group.Every group of parameter is set as an object function by the way that set of weights is legal, form is as follows：
f1(x)=w1y1(x)+w2y2(x)
f2(x)=w3y3(x)+w4y4(x)
Wherein, wiFor each subhead scalar functions yi(x) weighted factor, value are decided by the order of magnitude of the objectives and important
Degree.
Can obtain multi-goal optimizing function according to the weighted factor of each subhead scalar functions determined above is：
In the present embodiment, the correspondence weighted factor that the second-order response surface model calculates the optimization aim is first passed through,
Then combination is weighted according to the weighted factor and the optimization aim, to establish multi-goal optimizing function, will determined in this way
Position Parametric optimization problem converts for the optimization problem of two targets such as knuckle positioning and stub positioning.
Further, with reference to Fig. 8, the 7th embodiment of front suspension optimum design method of the present invention, based on above-mentioned first to the
Six embodiments, the step S50 include：
Step S51 obtains the variation range of the coordinates of targets parameter；
In the present embodiment, multi-goal optimizing function is：
Obtain above formula xiVariation range
Step S52 optimizes the coordinates of targets parameter in the variation range by multi-objective optimization algorithm；
Step S53, it is excellent that the coordinates of targets parameter after the multi-objective optimization algorithm is optimized inputs the multiple target
Change function, to obtain the Optimum Design Results of optimization aim.
In the present embodiment, willX in rangeiValue is used using NCGA algorithms (multi-objective optimization algorithm)
After vicinity Crossover Strategy optimizes, according to multi-goal optimizing function calculation optimization result.
In the present embodiment, the variation range for first obtaining coordinates of targets parameter, by multi-objective optimization algorithm to changing model
Coordinates of targets parameter in enclosing optimizes, then the coordinates of targets parameter input multiple target after multi-objective optimization algorithm is optimized is excellent
Change function, to obtain the Optimum Design Results of optimization aim, in this way ensure elitism strategy, density value valuation strategy, quickly it is non-bad
While ordering strategy and raising calculating speed, the exploration performance of algorithm is significantly improved, optimum results are more reasonable.
In addition, the embodiment of the present invention also proposes a kind of front suspension optimization design device, the front suspension optimization design device
Including：Memory, processor and it is stored in the Optimized Program that can be run on the memory and on the processor, institute
It states and realizes the front suspension optimum design method described in as above each embodiment when Optimized Program is executed by the processor
Step.
In addition, the embodiment of the present invention also proposes a kind of storage medium, it is stored with Optimized Program on the storage medium,
The step of the front suspension optimum design method described in as above each embodiment is realized when the Optimized Program is executed by processor
Suddenly.
It should be noted that herein, the terms "include", "comprise" or its any other variant are intended to non-row
His property includes, so that process, method, article or system including a series of elements include not only those elements, and
And further include other elements that are not explicitly listed, or further include for this process, method, article or system institute it is intrinsic
Element.In the absence of more restrictions, the element limited by sentence "including a ...", it is not excluded that including this
There is also other identical elements in the process of element, method, article or system.
The embodiments of the present invention are for illustration only, can not represent the quality of embodiment.
Through the above description of the embodiments, those skilled in the art can be understood that above-described embodiment side
Method can add the mode of required general hardware platform to realize by software, naturally it is also possible to by hardware, but in many cases
The former is more preferably embodiment.Based on this understanding, technical scheme of the present invention substantially in other words does the prior art
Going out the part of contribution can be expressed in the form of software products, which is stored in one as described above
In storage medium (such as ROM/RAM, magnetic disc, CD), including some instructions use so that a station terminal equipment (can be mobile phone,
Computer, server or network equipment etc.) execute method described in each embodiment of the present invention.
It these are only the preferred embodiment of the present invention, be not intended to limit the scope of the invention, it is every to utilize this hair
Equivalent structure or equivalent flow shift made by bright specification and accompanying drawing content is applied directly or indirectly in other relevant skills
Art field, is included within the scope of the present invention.
Claims (9)
1. a kind of front suspension optimum design method, which is characterized in that the front suspension optimum design method includes the following steps：
Obtain the hard spot coordinate of Simulation Calculation；
Sensitivity analysis is carried out to the coordinate parameters of the hard spot coordinate, and obtains coordinates of targets parameter；
The second-order response surface model of the Simulation Calculation is built according to the coordinates of targets parameter；
According to the second-order response surface model foundation multi-goal optimizing function；
The Optimum Design Results of optimization aim are obtained by the multi-goal optimizing function.
2. front suspension optimum design method as described in claim 1, which is characterized in that the coordinate to the hard spot coordinate
Parameter carries out sensitivity analysis, and the step of obtaining coordinates of targets parameter includes：
Determine the variation range of the coordinate parameters of the hard spot coordinate；
According to the variation range, the coordinate parameters are emulated by dynamics analysis software；
The coordinate parameters for meeting sensitivity condition are chosen according to simulation result, as coordinates of targets parameter.
3. front suspension optimum design method as described in claim 1, which is characterized in that described chosen according to simulation result meets
The step of coordinate parameters of sensitivity condition includes：
The influence value to the Simulation Calculation of the coordinate parameters is determined according to the simulation result；
Judge whether the coordinate parameters meet the sensitivity condition according to the influence value, and chooses and meet the sensitivity
The above-mentioned coordinate parameters of condition are as coordinates of targets parameter.
4. front suspension optimum design method as described in claim 1, which is characterized in that described according to the coordinates of targets parameter
The step of second-order response surface model for building the Simulation Calculation includes：
Central Composite is carried out to the coordinates of targets parameter to test to obtain test result；
The second-order response surface model of the optimization aim is determined according to the test result.
5. front suspension optimum design method as described in claim 1, which is characterized in that the front suspension optimum design method is also
Including：
Obtain the regression coefficient of the second-order response surface model；
The multiple correlation coefficient for calculating the regression coefficient, to determine the second-order response surface model according to the multiple correlation coefficient
Accuracy；
When the accuracy is unsatisfactory for precise requirements, the second-order response surface model is corrected.
6. front suspension optimum design method as described in claim 1, which is characterized in that described according to the second-order response face mould
Type establishes the step of multi-goal optimizing function and includes：
The correspondence weighted factor of the optimization aim is calculated by the second-order response surface model；
It is weighted combination according to the weighted factor and the optimization aim, to establish multi-goal optimizing function.
7. front suspension optimum design method as described in claim 1, which is characterized in that described to pass through the multiple-objection optimization letter
Number obtain optimization aims Optimum Design Results the step of include：
Obtain the variation range of the coordinates of targets parameter；
The coordinates of targets parameter in the variation range is optimized by multi-objective optimization algorithm；
The coordinates of targets parameter after the multi-objective optimization algorithm is optimized inputs the multi-goal optimizing function, to obtain
The Optimum Design Results of optimization aim.
8. a kind of front suspension optimization design device, which is characterized in that the front suspension optimization design device includes：Memory, place
It manages device and is stored in the Optimized Program that can be run on the memory and on the processor, the Optimized Program
The step of front suspension optimum design method as described in any one of claim 1 to 7 is realized when being executed by the processor.
9. a kind of computer readable storage medium, which is characterized in that be stored with optimization on the computer readable storage medium and set
It has the records of distance by the log sequence, realizes that the front suspension as described in any one of claim 1 to 7 is excellent when the Optimized Program is executed by processor
The step of changing design method.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810560437.9A CN108446528A (en) | 2018-06-01 | 2018-06-01 | Front suspension optimum design method, device and computer readable storage medium |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810560437.9A CN108446528A (en) | 2018-06-01 | 2018-06-01 | Front suspension optimum design method, device and computer readable storage medium |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN108446528A true CN108446528A (en) | 2018-08-24 |
Family
ID=63206543
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201810560437.9A Pending CN108446528A (en) | 2018-06-01 | 2018-06-01 | Front suspension optimum design method, device and computer readable storage medium |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN108446528A (en) |
Cited By (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109808467A (en) * | 2019-03-04 | 2019-05-28 | 南京航空航天大学 | A kind of cactus biomimetic features anticollision car door and its optimum design method |
| CN111241628A (en) * | 2020-01-05 | 2020-06-05 | 开沃新能源汽车集团有限公司 | Method for optimizing suspension positioning parameters |
| CN111639388A (en) * | 2020-04-28 | 2020-09-08 | 东风汽车集团有限公司 | Method and system for simulating parameters of elastic element of automobile |
| CN112339884A (en) * | 2020-09-25 | 2021-02-09 | 广州汽车集团股份有限公司 | Stabilizer bar arrangement position determining method and device and readable storage medium |
| CN117725765A (en) * | 2024-02-07 | 2024-03-19 | 中汽研汽车检验中心（天津）有限公司 | Vehicle suspension multi-objective optimization method, device and medium based on response analysis |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US8504331B2 (en) * | 2009-04-03 | 2013-08-06 | Airbus Operations (S.A.S.) | Aerodynamic model identification process for aircraft simulation process |
| CN107590312A (en) * | 2017-08-11 | 2018-01-16 | 中国第汽车股份有限公司 | A kind of check method of the front suspension divertical motion based on skeleton pattern |
| CN107609303A (en) * | 2017-09-28 | 2018-01-19 | 北京新能源汽车股份有限公司 | Suspension optimization method and system for vehicle |
| CN107784133A (en) * | 2016-08-25 | 2018-03-09 | 上汽通用汽车有限公司 | The collocation method of McPherson front suspension hard spot |
| CN107871036A (en) * | 2017-09-28 | 2018-04-03 | 北京新能源汽车股份有限公司 | Sensitivity analysis method and system for suspension hard point |
-
2018
- 2018-06-01 CN CN201810560437.9A patent/CN108446528A/en active Pending
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US8504331B2 (en) * | 2009-04-03 | 2013-08-06 | Airbus Operations (S.A.S.) | Aerodynamic model identification process for aircraft simulation process |
| CN107784133A (en) * | 2016-08-25 | 2018-03-09 | 上汽通用汽车有限公司 | The collocation method of McPherson front suspension hard spot |
| CN107590312A (en) * | 2017-08-11 | 2018-01-16 | 中国第汽车股份有限公司 | A kind of check method of the front suspension divertical motion based on skeleton pattern |
| CN107609303A (en) * | 2017-09-28 | 2018-01-19 | 北京新能源汽车股份有限公司 | Suspension optimization method and system for vehicle |
| CN107871036A (en) * | 2017-09-28 | 2018-04-03 | 北京新能源汽车股份有限公司 | Sensitivity analysis method and system for suspension hard point |
Non-Patent Citations (1)
| Title |
|---|
| 钮彦磊 等: "ADAMS/Car与Insight在双横臂悬架优化设计中的应用", 《制造业自动化》 * |
Cited By (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109808467A (en) * | 2019-03-04 | 2019-05-28 | 南京航空航天大学 | A kind of cactus biomimetic features anticollision car door and its optimum design method |
| CN111241628A (en) * | 2020-01-05 | 2020-06-05 | 开沃新能源汽车集团有限公司 | Method for optimizing suspension positioning parameters |
| CN111639388A (en) * | 2020-04-28 | 2020-09-08 | 东风汽车集团有限公司 | Method and system for simulating parameters of elastic element of automobile |
| CN112339884A (en) * | 2020-09-25 | 2021-02-09 | 广州汽车集团股份有限公司 | Stabilizer bar arrangement position determining method and device and readable storage medium |
| CN112339884B (en) * | 2020-09-25 | 2021-10-29 | 广州汽车集团股份有限公司 | Stabilizer bar arrangement position determining method and device and readable storage medium |
| CN117725765A (en) * | 2024-02-07 | 2024-03-19 | 中汽研汽车检验中心（天津）有限公司 | Vehicle suspension multi-objective optimization method, device and medium based on response analysis |
| CN117725765B (en) * | 2024-02-07 | 2024-05-31 | 中汽研汽车检验中心（天津）有限公司 | Vehicle suspension multi-objective optimization method, device and medium based on response analysis |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN108446528A (en) | Front suspension optimum design method, device and computer readable storage medium | |
| Gao et al. | Robust sampled-data $ H_ {\infty} $ control for vehicle active suspension systems | |
| US6741957B1 (en) | Analytical tire model for vehicle durability and ride comfort analysis | |
| CN106909743B (en) | Macpherson suspension hard point coordinate optimization method based on multi-target particle swarm optimization | |
| CN108959686A (en) | A kind of correction method for finite element model based on sensitivity analysis | |
| CN116090111B (en) | Automobile leaf spring fatigue life prediction method based on deep learning model | |
| CN111931285A (en) | Vehicle lightweight forward design method based on multidisciplinary multi-objective optimization | |
| CN105846925A (en) | Three-dimensional MIMO OAT channel modeling method and apparatus | |
| CN115563714A (en) | Method for verifying vibration characteristic parameters of key positions of complex equipment | |
| CN116861553A (en) | Suspension hard point optimization method and device, electronic equipment and storage medium | |
| CN114065387A (en) | Frame optimization method and device, computer equipment and storage medium | |
| CN111639388B (en) | Method and system for simulating parameters of elastic element of automobile | |
| CN118036183A (en) | Vehicle suspension parameter identification method, device and equipment based on genetic algorithm | |
| CN113095513A (en) | Double-layer fair federal learning method, device and storage medium | |
| CN115935761B (en) | Reliability simulation method and device for equipment, computer equipment and storage medium | |
| CN111371604A (en) | Equipment system network reconstruction method, system, machine readable medium and equipment | |
| CN114993604B (en) | Wind tunnel balance static calibration and measurement method based on deep learning | |
| CN109614658A (en) | A kind of the multidisciplinary optimization platform and optimization method of vehicle air suspension system | |
| CN113627074A (en) | Ground wave propagation delay prediction method based on transfer learning | |
| Li et al. | An Improved Genetic Algorithm to Optimize Spatial Locations for Double‐Wishbone Type Suspension System with Time Delay | |
| CN117725668A (en) | Data-driven vehicle suspension K & C system and construction method | |
| Paluskar et al. | Taguchi method (DOE) based performance optimization of a three link rigid axle passenger car suspension using MBD simulations | |
| CN114861335B (en) | Calibration method of automobile dynamics calculation model and related equipment | |
| CN114970228B (en) | Lightweight design analysis method and related device for cross arm of extra-heavy chassis suspension | |
| Zhu et al. | [Retracted] Research on Finite Element Structure of Vehicle Suspension Control Arm Based on Neural Network Sensing Control |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| RJ01 | Rejection of invention patent application after publication |
Application publication date: 20180824 |
|
| RJ01 | Rejection of invention patent application after publication |