CN113239451A - Matlab program-based passenger vehicle drive shaft arrangement checking method - Google Patents

Matlab program-based passenger vehicle drive shaft arrangement checking method Download PDF

Info

Publication number
CN113239451A
CN113239451A CN202110316753.3A CN202110316753A CN113239451A CN 113239451 A CN113239451 A CN 113239451A CN 202110316753 A CN202110316753 A CN 202110316753A CN 113239451 A CN113239451 A CN 113239451A
Authority
CN
China
Prior art keywords
driving shaft
displacement
value range
angle
coordinates
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.)
Granted
Application number
CN202110316753.3A
Other languages
Chinese (zh)
Other versions
CN113239451B (en
Inventor
孟帅
方超
孙波
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Zhejiang Hozon New Energy Automobile Co Ltd
Original Assignee
Zhejiang Hozon New Energy Automobile Co Ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Zhejiang Hozon New Energy Automobile Co Ltd filed Critical Zhejiang Hozon New Energy Automobile Co Ltd
Priority to CN202110316753.3A priority Critical patent/CN113239451B/en
Publication of CN113239451A publication Critical patent/CN113239451A/en
Application granted granted Critical
Publication of CN113239451B publication Critical patent/CN113239451B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/10Geometric CAD
    • G06F30/15Vehicle, aircraft or watercraft design
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/10Geometric CAD
    • G06F30/17Mechanical parametric or variational design
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/20Design reuse, reusability analysis or reusability optimisation
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02TCLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
    • Y02T90/00Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation

Landscapes

  • Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Geometry (AREA)
  • Theoretical Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Computer Hardware Design (AREA)
  • General Engineering & Computer Science (AREA)
  • Evolutionary Computation (AREA)
  • Mathematical Analysis (AREA)
  • Pure & Applied Mathematics (AREA)
  • Mathematical Optimization (AREA)
  • Computational Mathematics (AREA)
  • Automation & Control Theory (AREA)
  • Aviation & Aerospace Engineering (AREA)
  • Automobile Manufacture Line, Endless Track Vehicle, Trailer (AREA)
  • Shafts, Cranks, Connecting Bars, And Related Bearings (AREA)

Abstract

The invention belongs to the technical field of vehicle manufacturing, and particularly relates to a method for checking arrangement of a driving shaft of a passenger vehicle based on a Matlab program. The method comprises the following steps: s1, establishing a mathematical model of the limit values of the inner node displacement and the inner node rotation angle, and obtaining the value range of the length L of the driving shaft; s2, establishing a mathematical model of the limit value of the outer joint rotation angle, and correcting the value range of the driving shaft length L obtained in the step S1 according to the obtained limit value range of the included angle between the outer joint axis and the shaft lever; s3, establishing a mathematical model of the offset condition of the contact point displacement of the inner joint bearing, and correcting the value range of the length L of the drive shaft obtained in the step S2 according to the limit value range of the contact point displacement of the bearing; s4, selecting the optimal solution of the driving shaft length L according to the result of the step S3. The method has the characteristics of better calculation precision and capability of effectively improving the checking efficiency of the arrangement of the driving shafts.

Description

Matlab program-based passenger vehicle drive shaft arrangement checking method
Technical Field
The invention belongs to the technical field of vehicle manufacturing, and particularly relates to a method for checking arrangement of a driving shaft of a passenger vehicle based on a Matlab program.
Background
At present, in the general arrangement process of a power assembly of a passenger vehicle with a front driving shaft, the position of the power assembly and the motion rule of a chassis part need to be considered to meet the motion requirement of the driving shaft. The main problems involved are the confirmation of the pitch and the axial length L of the drive shaft. The feasibility analysis process of the traditional method is as follows:
the method comprises the steps of determining theoretical inner nodes of left and right driving shafts after the position of a power assembly is determined, estimating the length L of the driving shafts by using the theoretical inner nodes, finishing data of an initial driving shaft at a plurality of characteristic positions by the L, measuring values of inner joint displacement, rotation angle and outer joint rotation angle, filling the values into a checking template, calculating the possible ideal length of the driving shafts according to the distribution condition of each curve in the template, correcting the L according to experience to obtain the corrected length L1, repeating the checking process, judging whether the L1 is the optimal solution according to the position of each limit value model of the characteristic curve in the template, continuing correcting to L2 if the L1 is not the optimal solution, and then checking and comparing until the optimal solution of the length is found, wherein the process usually needs to consume great efforts and is low in efficiency.
Therefore, it is necessary to design a method for checking the driving shaft arrangement of a passenger car, which has better calculation accuracy and can effectively improve the efficiency of checking the driving shaft arrangement.
For example, the method for checking the slip deflection angle of the driving shaft and the system for checking the slip deflection angle of the driving shaft described in chinese patent application No. CN202010865995.3 include: acquiring a hard point of the whole vehicle; acquiring the position of the mass center of the power assembly; acquiring adjustment parameter information of a suspension component; and acquiring a driving shaft checking model according to the whole vehicle hard point, the mass center position and the adjustment parameter information, and acquiring a driving shaft enveloping body and a sliding swing angle curve according to the driving shaft checking model. Although the influence parameters such as the displacement of the power assembly, the four-wheel parameter adjustment and the like are established in the drive shaft checking model, the contribution and the influence of each parameter variable on the slip quantity, the swing angle and the envelope of the drive shaft can be effectively analyzed, and meanwhile, the actual whole vehicle layout and parameters can be combined to carry out specific analysis checking, so that the problems of insufficient design or excessive design caused by design checking of all vehicle types only by empirical values are solved, and the obtained drive shaft envelope and the slip swing angle curve are more accurate.
Disclosure of Invention
The invention provides a method for checking arrangement of a driving shaft of a passenger vehicle based on a Matlab program, which has better calculation accuracy and can effectively improve the arrangement and checking efficiency of the driving shaft, and aims to solve the problem of low working efficiency caused by the fact that great energy is consumed to carry out arrangement and checking on the driving shaft and the calculation of the shaft length of the driving shaft in the process of arranging a power assembly of the passenger vehicle with a front driving shaft in the prior art.
In order to achieve the purpose, the invention adopts the following technical scheme:
a passenger vehicle driving shaft arrangement checking method based on a Matlab program comprises the following steps:
s1, establishing a mathematical model of the limit values of the inner node displacement and the inner node rotation angle, and obtaining the value range of the length L of the driving shaft;
s2, establishing a mathematical model of the limit value of the outer joint rotation angle, and correcting the value range of the length L of the driving shaft obtained in the step S1 according to the obtained limit value range of the included angle between the outer joint axis and the shaft rod;
s3, establishing a mathematical model of the offset condition of the contact point displacement of the inner joint bearing, and correcting the value-taking range of the length L of the driving shaft obtained in the step S2 according to the limit value range of the contact point displacement of the bearing;
s4, selecting the optimal solution of the driving shaft length L according to the result of the step S3.
Preferably, step S1 includes the steps of:
establishing a coordinate system, wherein the abscissa is displacement, the ordinate is an inner pitch angle, the displacement is sequentially connected into straight lines according to known points 1, 2, 3, 4 and 5 provided by a supplier to form a limit value mathematical model, and linear equations of the straight lines L1, L2 and L3 are solved through the coordinates of the known points; the straight lines L1, L2 and L3 are obtained according to design experience;
setting a theoretical inner node of a left driving shaft as O1, an actual inner node of the driving shaft as Ot, an outer node of the driving shaft as Mli and any point on the outer node of the driving shaft as Mlo;
setting the axial length of a driving shaft, namely the distance between Mli and Ot as L, making an MliH perpendicular O1Ot extension line passing through Mli and at a point H, setting an angle MliOtH as theta, O1Ot as t, namely an inner node rotation angle position theta, and setting the displacement of an inner node as t;
according to the known coordinates of O1Ot and Y axis in the space coordinate system in parallel and the coordinates of O1(x1, y1, z1), Mli (xi, yi, zi) and Mlo (xo, yo, zo), the coordinates of the H point are (x1, yi, z1) according to the space geometrical relation, and the parametric equation of t and theta is solved:
Figure BDA0002990622570000031
Figure BDA0002990622570000032
let the equation of the straight line L be: theta is equal to kt + b,
Figure BDA0002990622570000033
substituting equations (1) and (2) into the equation of line L yields the following equation:
Figure BDA0002990622570000034
order to
Figure BDA0002990622570000035
Wherein f (L) represents a series of straight line groups that are parallel to the model boundary straight lines; l obtained when f (L) is 0 is a boundary point.
Preferably, step S1 further includes the steps of:
derivation of f (l) yields:
Figure BDA0002990622570000036
when k < 0, F' (L) < 0, F (L) is a decreasing function; if the root of f (L) ═ 0 is Lx3, substituting the coordinates (t, theta) corresponding to all the points on the left side of the straight line L3 into f (L), wherein f (L) < 0, and obtaining the value range of L which is L > Lx 3;
when k is more than 0, because theta is less than 30 degrees, k is more than 0.1 and L is more than 300 in the actual engineering, the method obtains
Figure BDA0002990622570000041
I.e. F' (L) > 0, to obtain F (L) as an increasing function; if the calculated roots of f (L) ═ 0 are Lxb1 and Lxb2 by substituting k and b in the equations of the straight lines L1 and L2, respectively, the coordinates (t, θ) corresponding to all points on the right side of the straight lines L1 and L2 are substituted into f (L), and if f (L) < 0, the value range of L is L < min (Lxb1, lxb 2);
finally, the value range of the obtained L is Lx3 < L < min (Lxb1, lxb 2).
Preferably, step S2 includes the steps of:
the pivot angle of the outer section of the driving shaft is set to be phi, and phi is less than phi0Said psi0To allow for the pendulum angle, the coordinates of Ot from step S1 are (x1, y1-t, z1), and the vector angle ψ is found from the space vectors MliMlo and OtMli:
Figure BDA0002990622570000042
conversion according to formula (3) gives:
Figure BDA0002990622570000043
wherein, in the formula (3), the | OtMli | is L; the value range of L is obtained by using an inequality (4).
Preferably, step S3 includes the steps of:
angular conversion displacement tθTan (θ) × λ/2, λ being a conversion coefficient determined by the inner node structure of the drive shaft;
the geometric relationship of the inner joint structure of the driving shaft is as follows:
when the angular conversion displacement is larger than 0, the bearing contact point displacement tz=t+tθ
When the angular conversion displacement is less than 0, the bearing contact point displacement tz=t-tθ
If the design requirement tzShould be at tmAnd tnIn between, let tz>tmAnd t isz<tn
According to design requirements tzAnd selecting the value range of the L meeting the requirement.
Preferably, the straight lines L1, L2 and L3 in step S1 form a safety boundary, and the point coordinates (t, θ) of the inner node of the driving shaft during the slippage process of each working condition are within the safety boundary.
Preferably, the driving shaft outer pitch pivot angle psi is an included angle between an outer pitch axis and a driving shaft lever axis, and the allowable working pivot angle psi0Is 50 deg.
Preferably, the bearing contact point displacement tzThe displacement t converted from the inner joint displacement t and the inner joint rotation angleθIs calculated to obtain tzThe value range of (A) is from-20 mm to +20 mm.
Compared with the prior art, the invention has the beneficial effects that: the checking method of the invention establishes three mathematical models according to the space motion rule of the driving shaft: the mathematical model of inner node displacement and inner pitch corner, the mathematical model of outer pitch corner and the mathematical model of bearing contact point displacement to according to above mathematical model can be through the coordinate of some points in the known three-dimensional data, the value range of once only calculating the driving shaft length accurately, compare in traditional check-up process, work efficiency improves and is showing, at the in-process of whole car modularized design, has more obvious advantage and meaning.
Drawings
FIG. 1 is a schematic diagram of a mathematical model of the limits of displacement and angle of rotation of an inner node according to the present invention;
FIG. 2 is a schematic diagram of a mathematical model of the limits of the outer joint rotation angle of the present invention;
FIG. 3 is a schematic view of a drive shaft according to the present invention;
FIG. 4 is a schematic main interface diagram of a design program for optimizing the axial length of the drive shaft according to the embodiment 1 of the present invention;
fig. 5 is a schematic main interface diagram of a drive shaft axial length fixed parameter optimization design program in an actual application process in embodiment 1 of the present invention;
FIG. 6 is a graph showing the effect of outputting a characteristic value when the middle value of the range of values of the length of the left driving shaft in FIG. 5 is 450 mm;
FIG. 7 is a graph showing the effect of outputting a characteristic value when the middle value of the range of values of the length of the left driving shaft in FIG. 5 is 452 mm;
fig. 8 is a diagram showing an effect of outputting a characteristic value when the middle value of the range of the length of the right drive shaft in fig. 5 is 361 mm.
Detailed Description
In order to more clearly illustrate the embodiments of the present invention, the following description will explain the embodiments of the present invention with reference to the accompanying drawings. It is obvious that the drawings in the following description are only some examples of the invention, and that for a person skilled in the art, other drawings and embodiments can be derived from them without inventive effort.
Example 1:
a passenger vehicle driving shaft arrangement checking method based on a Matlab program comprises the following steps:
s1, establishing a mathematical model of the limit values of the inner node displacement and the inner node rotation angle, and obtaining the value range of the length L of the driving shaft;
s2, establishing a mathematical model of the limit value of the outer joint rotation angle, and correcting the value range of the length L of the driving shaft obtained in the step S1 according to the obtained limit value range of the included angle between the outer joint axis and the shaft rod;
s3, establishing a mathematical model of the offset condition of the contact point displacement of the inner joint bearing, and correcting the value-taking range of the length L of the driving shaft obtained in the step S2 according to the limit value range of the contact point displacement of the bearing;
s4, selecting the optimal solution of the driving shaft length L according to the result of the step S3.
Further, step S1 includes the following steps:
as shown in fig. 1, a coordinate system is established, the abscissa is displacement, the ordinate is inner pitch angle, the displacement is sequentially connected into straight lines according to known points 1, 2, 3, 4 and 5 provided by a supplier to form a limit value mathematical model, and the linear equations of the straight lines L1, L2 and L3 are obtained through the coordinates of the known points; the straight lines L1, L2 and L3 are obtained according to design experience;
the limit boundary is a set of (t, theta) coordinates that can meet the technical requirements, and is generally provided by a supplier; the straight lines L1, L2 and L3 form a safety boundary, the safety boundary is generally obtained by offsetting a limit boundary by a certain safety amount, and point coordinates (t, theta) of the inner joint of the driving shaft in the sliding process of each working condition are all in the safety boundary.
Further, as shown in fig. 3, a theoretical inner node of the left driving shaft is O1, an actual inner node of the driving shaft is Ot, an outer node of the driving shaft is Mli, and any point on the outer node of the driving shaft (an end point of the outer node in fig. 4) is Mlo;
setting the axial length of a driving shaft, namely the distance between Mli and Ot as L, making an MliH perpendicular O1Ot extension line passing through Mli and at a point H, setting an angle MliOtH as theta, O1Ot as t, namely an inner node rotation angle position theta, and setting the displacement of an inner node as t;
according to the known coordinates of O1Ot and Y axis in the space coordinate system in parallel and the coordinates of O1(x1, y1, z1), Mli (xi, yi, zi) and Mlo (xo, yo, zo), the coordinates of the H point are (x1, yi, z1) according to the space geometrical relation, and the parametric equation of t and theta is solved:
Figure BDA0002990622570000071
Figure BDA0002990622570000072
let the equation of the straight line L be: theta is equal to kt + b,
Figure BDA0002990622570000073
substituting equations (1) and (2) into the equation of line L yields the following equation:
Figure BDA0002990622570000074
order to
Figure BDA0002990622570000075
Wherein f (L) represents a series of straight line groups that are parallel to the model boundary straight lines; l obtained when f (L) is 0 is a boundary point.
Further, step S1 includes the following steps:
derivation of f (l) yields:
Figure BDA0002990622570000081
when k is less than 0 (corresponding to the equation L3), F' (L) < 0, F (L) is a decreasing function; if the root of f (L) ═ 0 is Lx3, substituting the coordinates (t, theta) corresponding to all the points on the left side of the straight line L3 into f (L), wherein f (L) < 0, and obtaining the value range of L which is L > Lx 3;
when k is greater than 0 (corresponding to the linear equations L1 and L2), since theta is less than 30 degrees in engineering practice, k is greater than 0.1, and L is greater than 300, the method is obtained
Figure BDA0002990622570000082
Figure BDA0002990622570000083
I.e. F' (L) > 0, to obtain F (L) as an increasing function; if the roots obtained by substituting k and b in the equations of the straight lines L1 and L2 into f (L) ═ 0 are Lxb1 and Lxb2, respectively, the coordinates (t, θ) corresponding to all points on the right side of the straight lines L1 and L2 are substituted into f (L), and if f (L) < 0, the roots obtained from L are substituted into f (L) < 0, and the roots obtained from L are substituted into LThe value range is L < min (Lxb1, lxb 2);
finally, the value range of the obtained L is Lx3 < L < min (Lxb1, lxb2)
Further, as shown in fig. 2, step S2 includes the following steps:
the pivot angle of the outer section of the driving shaft is set to be phi, and phi is less than phi0Said psi0To allow for the pendulum angle, the coordinates of Ot from step S1 are (x1, y1-t, z1), and the vector angle ψ is found from the space vectors MliMlo and OtMli:
Figure BDA0002990622570000084
conversion according to formula (3) gives:
Figure BDA0002990622570000085
wherein, in the formula (3), the | OtMli | is L; the value range of L is obtained by using an inequality (4).
The driving shaft outer pitch pivot angle psi is a included angle between an outer pitch axis and a driving shaft lever axis, and the outer pitch pivot angle psi is smaller than an allowable working pivot angle psi0Said allowable working pivot angle psi0Is generally taken to be about 50 ° or 50 °. In FIG. 2, the abscissa represents the wheel center oscillation position corresponding to the outer section of the driving shaft, and the ordinate represents the pivot angle psi under the corresponding operating conditioniAll psiiShould be at psi0Below the line.
Further, step S3 includes the following steps:
angular conversion displacement tθTan (θ) × λ/2, λ being a conversion coefficient determined by the inner node structure of the drive shaft;
the geometric relationship of the inner joint structure of the driving shaft is as follows:
when the angular conversion displacement is larger than 0, the bearing contact point displacement tz=t+tθ
When the angular conversion displacement is less than 0, the bearing contact point displacement tz=t-tθ
If the design requirement tzShould be at tmAnd tnIn between, let tz>tmAnd t isz<tn
According to design requirements tzAnd selecting the value range of the L meeting the requirement.
The bearing contact point displacement tzThe displacement t converted from the inner joint displacement t and the inner joint rotation angleθIs calculated to obtain tzThe value range of (A) is-20 mm to +20mm, and the specific requirements are given by suppliers.
And according to the system modeling, analysis and optimization objective function, optimizing the drive shaft design by adopting an MATLAB/GUI function and an optimization function in an optimization tool box. MATLAB is very powerful to matrix data processing ability, and GUI is one kind in MATALB and establishes the graphical user interface, can utilize its own controlling part, realizes the visualization of drive shaft design process through the language c, has good human-computer interaction function.
The invention utilizes the GUI function and the optimization tool box of MATLAB to establish a drive shaft axial length fixed parameter optimization design program, and the main interface is shown in figure 4. The position of the outer node is determined along with wheel jump, coordinates of the outer node are generated to excel by three-dimensional software at one time and then are imported into design software, namely, the value range of the left driving shaft and the right driving shaft can be output, a user can select any solution in the value range, and the software calculates according to the solution and outputs the solution to a model corresponding to the excel to form a judgment image.
And (4) realizing a corresponding solving strategy through software according to the mathematical modeling of the driving shaft length optimization design.
Solving strategies of the mathematical models of the inner node displacement and the inner node rotation angle:
coordinates of Mli (xi, yi, zi) and Mlo (xo, yo, zo) of characteristic positions are generated in a cata suspension motion DMU model, corresponding t and theta matrixes can be obtained by substituting equations (1) and (2) into which coordinate matrixes are substituted, a matrix of a solution of L can be obtained by substituting data in the matrixes into three linear equations f (L) and 0, and a solution of L is obtained as [ Lx3] < L < min ([ Lxb1], [ lxb2 ]).
Solving strategy of the mathematical model of the outer joint rotation angle:
as known from design experience, the allowable conditions of the outer pitch angle are easy to realize, so that the value range of L calculated by solving an equation is very large, the final value is not significant, and a strategy of a trial-and-error method is adopted in calculation: and (3) in the mathematical model of the inner node displacement and the inner node rotation angle, taking an integer in the value range of the solved L to test the value, and taking out the value meeting the inequality (4) to obtain the inequality solution meeting the requirement.
Solving strategy of the mathematical model of the contact point displacement of the inner joint bearing:
in the solving process of the mathematical model of the contact point displacement of the inner bearing, a trial value method is needed to solve, and the trial value strategy is the same as the solving strategy of the mathematical model of the outer joint rotation angle.
The trial value process of the mathematical model of the outer joint corner and the mathematical model of the contact point displacement of the inner joint bearing is easy to realize in computer calculation, so the trial value range can be properly enlarged, and the condition of no solution can be solved for satisfying the safety solution of the mathematical model of the limit value inner joint displacement and the inner joint corner due to the solution of the L in the mathematical model of the inner joint displacement and the inner joint corner, and the safety value can be chosen according to the actual engineering relaxation.
The specific operation steps of the MATLAB software program are as follows:
step 1, filling input conditions required by drive shaft design into an EXCEL template, wherein the EXCEL template comprises the following contents:
a, the theoretical inner node space coordinate of the driving shaft, as shown in Table 1, is defined as left side O1 and right side O2;
b, in the suspension motion process, the rack stroke and the tire jumping stroke of the steering gear under each typical working condition, the corresponding coordinates of the node center point Mli of the outer joint of the left driving shaft and any point Mlo on the axis of the outer joint are shown in tables 2 and 3, and the spatial motion position of the table 2 and the coordinates of the relevant points of the outer joint of the table 3 are in one-to-one correspondence;
c, driving shaft inner pitch slip curve boundary points, as shown in Table 4, wherein safety variables are filled in according to design experience values;
d, driving shaft outer section allowable working swing angle limit value, as shown in table 5;
e, coefficients of the driving shaft angle conversion displacement formula, as shown in Table 6
TABLE 1
Figure BDA0002990622570000111
TABLE 2
Time of exercise Steering (rack stroke) Jumping (wheel jumping stroke)
20
30
……
120
130
……
230
……
TABLE 3
Figure BDA0002990622570000121
TABLE 4
Figure BDA0002990622570000131
TABLE 5
Figure BDA0002990622570000132
TABLE 6
Figure BDA0002990622570000133
Figure BDA0002990622570000141
Wherein, table 1 is determined by the overall arrangement speciality of the whole vehicle, tables 2 and 3 are input conditions of the speciality of the chassis before the arrangement of the engine room, and the information in tables 4, 5 and 6 is determined by the self-attribute of the driving shaft and is generally provided by a supplier.
Step 2, importing the template information into a driving shaft design and analysis system, wherein the importing process is an import file module for running a program, the program automatically reads all information in the template, and the processing and calculation are carried out, and the process is as follows:
s21, reading the A, D, E information in the step 1, and checking the information with a theoretical range, if the information does not meet the requirement, suspending operation and reporting an error;
s22, reading A, B information in step 1: the coordinates of the outer node center and the point on the outer node axis are in one-to-one correspondence with the typical working conditions of suspension motion, the working conditions are numbered in the template, and the program correspondingly stores the serial numbers and the coordinates into fixed arrays in the data processing process. In the process of storage, the program judges invalid data, and the array finally used for calculation eliminates the invalid data. And the space coordinates of the right driving shaft under each working condition are obtained by utilizing the characteristic that the positions of the left outer joint and the right outer joint of the driving shaft are symmetrical.
S23, defining an inner joint rotation angle a, an inner joint displacement t and a driving shaft axial length L, setting a vertical line to a straight line O1Ot through O1, setting a foot to be H, setting an inner node to be Ot in the actual motion process of the driving shaft, and solving a parameter equation of a and t relative to L through a space triangle formed by space points Mli, O1 and H and Mli, Ot and H.
S24, reading the information in the step C to find out the limit range of the inner joint slip curve, wherein the range consists of three straight lines, and the equations of the three straight lines can be expressed by a and t;
and S25, combining the parameter equation in the step S23 and the linear equation set in the step S24 to form three equation sets related to L, finding a set of L values of L in the limit range of the slip curve through monotonicity judgment of a function related to the L equation, and defining the set as S.
S26, adding 20 to the maximum value of L in the set S, subtracting 20 from the minimum value to form a set Sa, obtaining an expression of the outer pitch swing angle relative to L through a linear vector formed by Mli and Mlo and a linear vector formed by Ot and Mli, substituting the set S1 into the expression to obtain a swing angle set N corresponding to the set Sa, comparing the set N with allowable swing angles, and excluding all L corresponding to the swing angles which do not meet the requirements, so that the rest parts in the set Sa are L sets meeting the requirements of the allowable working swing angles of the driving shaft, and the set is defined as S1.
S27, the bearing contact point displacement z is converted to an expression for L, and the relationship between z and L is provided by the drive spindle supplier, and is generally the relationship shown in table 6. Substituting the data in the Sa into the expression to obtain a set Z of the pivot angle corresponding to the set Sa, comparing the set Z with the allowable bearing contact point displacement, and excluding all L corresponding to the displacement which does not meet the requirement, so that the rest part in the set Z is the set L meeting the allowable working pivot angle requirement of the driving shaft, and the set is defined as S2.
And S28, respectively outputting the sets S, S1 and S2, wherein the intersection of the sets is a set of the axial length of the driving shaft, and the set is a feasible scheme in design. Repeating the above steps can also find a feasible solution for the right drive shaft. During the output of S1 and S2, the computational analysis system will make decisions on two sets: if the result of S1 is that the set is empty, it is indicated that the L value does not meet the requirement of the swing angle of the outer joint, the system calculates the corresponding working swing angle output as the reference by the self-output reference value, namely the maximum value and the minimum value of the extension range of the set S (the maximum value and the minimum value in the set are extended by 8 units); if the result of S2 is finally judged to be an empty set, which indicates that the L value of the bearing contact point displacement requirement is not met, the system calculates the corresponding bearing contact point displacement output as the reference by the self-output reference value, namely the maximum value and the minimum value of the extension range (extending the maximum value and the minimum value in the set by 8 units) of the set S.
And S29, respectively storing working condition numbers, inner joint rotation angles, outer joint swing angles and bearing contact point displacement numbers corresponding to Sa, wherein the three numbers are defined as A, T, Zt.
And 3, filling the inner joint rotation angle, the outer joint swing angle and the bearing contact displacement corresponding to the selected driving shaft length scheme into a template, observing the position of the corresponding characteristic curve in the check model, and outputting a plurality of groups of images corresponding to L for comparison to find an optimal solution. The specific process is as follows:
s31, finding a reasonable L value through the set S, S1 and the set S2, inputting the value into a driving shaft design and analysis system, addressing the position of L in A, T, Zt by the system, sequentially outputting the numerical values of L in A, T, Zt corresponding to each working condition to corresponding positions in a corresponding EXCEL template, completing the process by a writing algorithm of MATLAB to EXCEL, and finally automatically displaying a characteristic curve in the EXCEL;
and S32, if the characteristic curve displayed above is not the optimal solution, respectively inputting the L of the multiple feasible schemes into the system, outputting the sheet in the EXCEL after the system processes, and comparing the sheet images corresponding to the L to find the optimal solution, namely the final scheme of the driving shaft length.
Taking the driving shaft parameters of the H7 four-wheel drive vehicle as an example, the specific data are as follows:
inner segment slip limiting model
Figure BDA0002990622570000161
Figure BDA0002990622570000171
(II) outer joint swing angle requirement
Allowable working swing angle 48
(III) Angle conversion Displacement requirement
Figure BDA0002990622570000172
The data obtained by importing the above templates and known condition data into software and solving the data is shown in fig. 5.
According to the result shown in fig. 5, the middle value of the range of the length of the left driving shaft is 450mm, and the characteristic value is output, as shown in fig. 6.
It can be seen from fig. 6 that the bearing contact point displacement is mostly positive, so it is considered to lengthen the drive shaft to 452mm and output a characteristic value, as shown in fig. 7.
Then, according to the result shown in fig. 5, the middle value 361mm of the range of the length of the right driving shaft is selected, and the characteristic value is output, as shown in fig. 8 (when the shaft length is 362mm, the image in the checking model is equivalent to 361mm, which can be judged by the image shift rule).
Finally, according to the result, the shaft length parameters of the four-wheel drive left and right driving shafts are selected to be 452mm and 361mm respectively.
The invention can complete the operation and check of 33 characteristic position working conditions required by the check specification at one time, and input the characteristic data into the check template to form an image. The adopted software program can calculate the value range of the length of the driving shaft, and the sharing of the driving shaft can be realized by taking an intersection in the modularized design, so that the cost of the whole vehicle is reduced; completing the determination time of the axle length parameters of the left front drive shaft and the right front drive shaft of a single vehicle type: the operation time of a designer is only 3 minutes, the operation time of a system background is about 7 minutes when an output file is formed, and compared with the traditional parameter setting process, the design efficiency is greatly improved.
The checking method of the invention establishes three mathematical models according to the space motion rule of the driving shaft: the mathematical model of interior nodal point displacement and interior pitch corner, the mathematical model of outer pitch corner and the mathematical model of bearing contact displacement to can once only calculate the value range of driving shaft length through the coordinate of some points in the known three-dimensional data according to above mathematical model, compare in traditional check process, work efficiency improves and is showing, at the in-process of whole car modularized design, has more obvious advantage and meaning.
The foregoing has outlined rather broadly the preferred embodiments and principles of the present invention and it will be appreciated that those skilled in the art may devise variations of the present invention that are within the spirit and scope of the appended claims.

Claims (8)

1. A passenger vehicle driving shaft arrangement checking method based on a Matlab program is characterized by comprising the following steps:
s1, establishing a mathematical model of the limit values of the inner node displacement and the inner node rotation angle, and obtaining the value range of the length L of the driving shaft;
s2, establishing a mathematical model of the limit value of the outer joint rotation angle, and correcting the value range of the driving shaft length L obtained in the step S1 according to the obtained limit value range of the included angle between the outer joint axis and the shaft lever;
s3, establishing a mathematical model of the offset condition of the contact point displacement of the inner joint bearing, and correcting the value range of the length L of the drive shaft obtained in the step S2 according to the limit value range of the contact point displacement of the bearing;
s4, selecting the optimal solution of the driving shaft length L according to the result of the step S3.
2. The Matlab program based passenger car drive shaft arrangement checking method of claim 1, wherein the step S1 comprises the steps of:
establishing a coordinate system, wherein the abscissa is displacement, the ordinate is an inner pitch angle, the displacement is sequentially connected into straight lines according to known points 1, 2, 3, 4 and 5 provided by a supplier to form a limit value mathematical model, and linear equations of the straight lines L1, L2 and L3 are solved through the coordinates of the known points; the straight lines L1, L2 and L3 are obtained according to design experience;
setting a theoretical inner node of a left driving shaft as O1, an actual inner node of the driving shaft as Ot, an outer node of the driving shaft as Mli and any point on the outer node of the driving shaft as Mlo;
setting the axial length of a driving shaft, namely the distance between Mli and Ot as L, drawing MliH through Mli, and perpendicular to an O1Ot extension line at a point H, setting & lt MliOtH & gt & theta, O1Ot & gt & lt t, namely an inner node rotation angle position theta, and setting the displacement of an inner node as t;
according to the known coordinates of O1Ot and Y axis in the space coordinate system in parallel and the coordinates of O1(x1, y1, z1), Mli (xi, yi, zi) and Mlo (xo, yo, zo), the coordinates of the H point are (x1, yi, z1) according to the space geometrical relation, and the parametric equation of t and theta is solved:
Figure FDA0002990622560000021
Figure FDA0002990622560000022
let the equation of the straight line L be: theta is equal to kt + b,
Figure FDA0002990622560000023
substituting equations (1) and (2) into the equation of line L yields the following equation:
Figure FDA0002990622560000024
order to
Figure FDA0002990622560000025
Wherein f (L) represents a series of straight line groups that are parallel to the model boundary straight lines; l obtained when f (L) is 0 is a boundary point.
3. The Matlab program based passenger car drive axle layout checking method of claim 2, wherein the step S1 further comprises the steps of:
derivation of f (l) yields:
Figure FDA0002990622560000026
when k < 0, F' (L) < 0, F (L) is a decreasing function; if the root of f (L) ═ 0 is Lx3, substituting the coordinates (t, theta) corresponding to all the points on the left side of the straight line L3 into f (L), wherein f (L) < 0, and obtaining the value range of L which is L > Lx 3;
when k is more than 0, because theta is less than 30 degrees, k is more than 0.1 and L is more than 300 in the actual engineering, the method obtains
Figure FDA0002990622560000027
I.e. F' (L) > 0, to obtain F (L) as an increasing function; if the calculated roots of f (L) ═ 0 are Lxb1 and Lxb2 by substituting k and b in the equations of the straight lines L1 and L2, respectively, the coordinates (t, θ) corresponding to all points on the right side of the straight lines L1 and L2 are substituted into f (L), and if f (L) < 0, the value range of L is L < min (Lxb1, lxb 2);
finally, the value range of the obtained L is Lx3 < L < min (Lxb1, lxb 2).
4. The Matlab program based passenger car drive axle layout checking method of claim 3, wherein the step S2 comprises the steps of:
the pivot angle of the outer section of the driving shaft is set to be phi, and phi is less than phi0Said psi0To allow for the pendulum angle, the coordinates of Ot from step S1 are (x1, y1-t, z1), and the vector angle ψ is found from the space vectors MliMlo and OtMli:
Figure FDA0002990622560000031
conversion according to formula (3) gives:
Figure FDA0002990622560000032
wherein, in the formula (3), the | OtMli | is L; the value range of L is obtained by using an inequality (4).
5. The Matlab program based passenger car drive shaft arrangement checking method of claim 4, wherein the step S3 comprises the steps of:
angular conversion displacement tθTan (θ) × λ/2, λ being a conversion factor determined by the inner node structure of the drive shaft;
the geometric relationship of the inner joint structure of the driving shaft is as follows:
when the angular conversion displacement is larger than 0, the bearing contact point displacement tz=t+tθ
When the angular conversion displacement is less than 0, the bearing contact point displacement tz=t-tθ
If the design requirement tzShould be at tmAnd tnIn between, let tz>tmAnd t isz<tn
According to design requirements tzAnd selecting the value range of the L meeting the requirement.
6. The Matlab program based passenger car driving axle layout checking method of claim 2, wherein the straight lines L1, L2 and L3 in step S1 form a safety boundary, and the point coordinates (t, θ) of the driving axle inner node during each operating condition slip are within the safety boundary.
7. The Matlab program based passenger vehicle axle arrangement verification method of claim 4, wherein the axle shaft outer pitch pivot angle ψ is an angle between an outer pitch axis and an axle shaft axis of the axle shaft, and the allowable working pivot angle ψ0Is 50 deg.
8. Matlab program based passenger vehicle drive shaft arrangement verification method according to claim 5, wherein the bearing contact point displacement tzThe displacement t converted from the inner joint displacement t and the inner joint rotation angleθIs calculated to obtain tzThe value range of (A) is from-20 mm to +20 mm.
CN202110316753.3A 2021-03-24 2021-03-24 Matlab program-based passenger car drive shaft arrangement checking method Active CN113239451B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110316753.3A CN113239451B (en) 2021-03-24 2021-03-24 Matlab program-based passenger car drive shaft arrangement checking method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110316753.3A CN113239451B (en) 2021-03-24 2021-03-24 Matlab program-based passenger car drive shaft arrangement checking method

Publications (2)

Publication Number Publication Date
CN113239451A true CN113239451A (en) 2021-08-10
CN113239451B CN113239451B (en) 2023-06-16

Family

ID=77130720

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110316753.3A Active CN113239451B (en) 2021-03-24 2021-03-24 Matlab program-based passenger car drive shaft arrangement checking method

Country Status (1)

Country Link
CN (1) CN113239451B (en)

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
GB2017260A (en) * 1978-03-27 1979-10-03 Celanese Corp Carbon fibre reinforced composite drive shaft
JP2009097896A (en) * 2007-10-15 2009-05-07 Ntn Corp Shaft torque measuring device and measurement method of drive shaft
WO2017079263A1 (en) * 2015-11-02 2017-05-11 Amanda Bent-Bolt Company D/B/A Amanda Manufacturing Floating differential suspension system
CN107590322A (en) * 2017-08-25 2018-01-16 北京新能源汽车股份有限公司 Drive shaft bounce check method
CN210059397U (en) * 2019-06-14 2020-02-14 钟祥市金祥汽车半轴有限公司 Semi-axis correcting unit
CN111209649A (en) * 2019-08-26 2020-05-29 中国第一汽车股份有限公司 Drive shaft jumping checking method
CN112380620A (en) * 2020-11-09 2021-02-19 中国第一汽车股份有限公司 Power assembly positioning method based on CATIA hard spot parameter, computer and storage medium
CN112417609A (en) * 2020-12-15 2021-02-26 中国第一汽车股份有限公司 Steering transmission shaft optimization design method, computer equipment and storage medium

Patent Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
GB2017260A (en) * 1978-03-27 1979-10-03 Celanese Corp Carbon fibre reinforced composite drive shaft
JP2009097896A (en) * 2007-10-15 2009-05-07 Ntn Corp Shaft torque measuring device and measurement method of drive shaft
WO2017079263A1 (en) * 2015-11-02 2017-05-11 Amanda Bent-Bolt Company D/B/A Amanda Manufacturing Floating differential suspension system
CN107590322A (en) * 2017-08-25 2018-01-16 北京新能源汽车股份有限公司 Drive shaft bounce check method
CN210059397U (en) * 2019-06-14 2020-02-14 钟祥市金祥汽车半轴有限公司 Semi-axis correcting unit
CN111209649A (en) * 2019-08-26 2020-05-29 中国第一汽车股份有限公司 Drive shaft jumping checking method
CN112380620A (en) * 2020-11-09 2021-02-19 中国第一汽车股份有限公司 Power assembly positioning method based on CATIA hard spot parameter, computer and storage medium
CN112417609A (en) * 2020-12-15 2021-02-26 中国第一汽车股份有限公司 Steering transmission shaft optimization design method, computer equipment and storage medium

Non-Patent Citations (5)

* Cited by examiner, † Cited by third party
Title
MOHAMMAD CHOWDHURI,R. HOSSAIN: "Design Analysis of an Automotive Composite Drive Shaft", 《INTERNATIONAL JOURNAL OF ENGINEERING AND TECHNOLOGY》 *
张伟,陈忠廷等: "基于CATIA的乘用车前驱动轴布置优化", 《汽车工程师》 *
潘正东: "电动轿车等速万向节驱动轴设计与验证", 《农业装备与车辆工程》 *
董伟佳等: "经济型轿车传动轴布置校核设计方法", 《汽车工程师》 *
齐晓旭,郭健等: "一种驱动轴长度与整车的匹配设计方法", 《汽车技术》 *

Also Published As

Publication number Publication date
CN113239451B (en) 2023-06-16

Similar Documents

Publication Publication Date Title
CN105773620B (en) The trajectory planning control method of industrial robot free curve based on Double quaternions
JP5506617B2 (en) Robot control device
CN113341728B (en) Anti-noise type return-to-zero neural network four-wheel mobile mechanical arm trajectory tracking control method
CN112214916A (en) Method and system for restoring physical test bench simulation process based on virtual model
CN113580148B (en) Parallel robot kinematics calibration method based on equivalent kinematic chain
CN112109084A (en) Terminal position compensation method based on robot joint angle compensation and application thereof
CN112947298A (en) Method, system and terminal for optimally generating curved surface machining track of robot
CN111666643B (en) Method for determining contact performance of complex tooth surface
CN109255141B (en) Optimization method for cross section shape of forward conceptual design of automobile body
CN107609303A (en) The suspension optimization method and system of vehicle
Agrawal et al. Designing variable ackerman steering geometry for formula student race car
CN113239451A (en) Matlab program-based passenger vehicle drive shaft arrangement checking method
CN107832554A (en) A kind of McPherson suspension four-wheel parameters validation method based on CATIA parameterized modules
CN112100747B (en) Vehicle body framework topology optimization method, device, equipment and medium
Raghavan Suspension design for linear toe curves: a case study in mechanism synthesis
CN109614658B (en) Multidisciplinary optimization platform and multidisciplinary optimization method for automobile air suspension system
Zhang et al. Computation and optimization of rack and pinion steering mechanism considering kingpin parameters and tire side slip angle
CN115455599A (en) Method and device for estimating dynamic size of leaf spring
CN113146630B (en) Industrial robot milling error compensation method, system, device and medium
Chen et al. Parameter Optimization of Steering Trapezoid Mechanism Based on Hybrid Genetic Algorithm
Kodati et al. Kinematic analysis of the double wishbone suspension system
Dang et al. Compensation of springback error based on comprehensive displacement method in cold rotary forging for hypoid gear
Blanco et al. Methodology for the design of a double-wishbone suspension system by means of target cascading: problem statement
CN106777605B (en) Suspension side-view geometric motion analysis method and system
He et al. Volumetric error prediction and compensation of NC machine tool based on least square support vector machine

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
CB02 Change of applicant information
CB02 Change of applicant information

Address after: 314500 988 Tong Tong Road, Wu Tong Street, Tongxiang, Jiaxing, Zhejiang

Applicant after: Hezhong New Energy Vehicle Co.,Ltd.

Address before: 314500 988 Tong Tong Road, Wu Tong Street, Tongxiang, Jiaxing, Zhejiang

Applicant before: HOZON NEW ENERGY AUTOMOBILE Co.,Ltd.

Address after: 314500 988 Tong Tong Road, Wu Tong Street, Tongxiang, Jiaxing, Zhejiang

Applicant after: United New Energy Automobile Co.,Ltd.

Address before: 314500 988 Tong Tong Road, Wu Tong Street, Tongxiang, Jiaxing, Zhejiang

Applicant before: Hezhong New Energy Vehicle Co.,Ltd.

GR01 Patent grant
GR01 Patent grant
PE01 Entry into force of the registration of the contract for pledge of patent right

Denomination of invention: A Method for Checking the Layout of Passenger Car Drive Shafts Based on Matlab Program

Granted publication date: 20230616

Pledgee: Tongxiang Science and Technology Entrepreneurship Service Center Co.,Ltd.

Pledgor: United New Energy Automobile Co.,Ltd.

Registration number: Y2024330000709