CN113581201B - Potential field model-based collision avoidance control method and system for unmanned vehicle - Google Patents
Potential field model-based collision avoidance control method and system for unmanned vehicle Download PDFInfo
- Publication number
- CN113581201B CN113581201B CN202110829535.XA CN202110829535A CN113581201B CN 113581201 B CN113581201 B CN 113581201B CN 202110829535 A CN202110829535 A CN 202110829535A CN 113581201 B CN113581201 B CN 113581201B
- Authority
- CN
- China
- Prior art keywords
- vehicle
- control
- model
- collision avoidance
- potential field
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 38
- 238000006073 displacement reaction Methods 0.000 claims abstract description 13
- 230000008569 process Effects 0.000 claims abstract description 5
- 239000013598 vector Substances 0.000 claims description 23
- 238000005381 potential energy Methods 0.000 claims description 18
- 238000004422 calculation algorithm Methods 0.000 claims description 16
- 239000011159 matrix material Substances 0.000 claims description 10
- 238000005457 optimization Methods 0.000 claims description 9
- 230000001133 acceleration Effects 0.000 claims description 7
- 238000012937 correction Methods 0.000 claims description 4
- 206010039203 Road traffic accident Diseases 0.000 claims description 3
- 238000009795 derivation Methods 0.000 claims description 3
- 230000000694 effects Effects 0.000 claims description 3
- 229910052739 hydrogen Inorganic materials 0.000 claims 3
- 239000001257 hydrogen Substances 0.000 claims 3
- 150000002431 hydrogen Chemical class 0.000 claims 2
- UFHFLCQGNIYNRP-UHFFFAOYSA-N Hydrogen Chemical compound [H][H] UFHFLCQGNIYNRP-UHFFFAOYSA-N 0.000 claims 1
- 238000010586 diagram Methods 0.000 description 6
- 230000008859 change Effects 0.000 description 4
- 238000004364 calculation method Methods 0.000 description 3
- 230000003993 interaction Effects 0.000 description 3
- 238000005265 energy consumption Methods 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 230000004044 response Effects 0.000 description 2
- 230000004075 alteration Effects 0.000 description 1
- 239000010426 asphalt Substances 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 239000003795 chemical substances by application Substances 0.000 description 1
- 150000001875 compounds Chemical class 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 238000009499 grossing Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000005096 rolling process Methods 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 230000006641 stabilisation Effects 0.000 description 1
- 238000011105 stabilization Methods 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
- 239000000126 substance Substances 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
Images
Classifications
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B60—VEHICLES IN GENERAL
- B60W—CONJOINT CONTROL OF VEHICLE SUB-UNITS OF DIFFERENT TYPE OR DIFFERENT FUNCTION; CONTROL SYSTEMS SPECIALLY ADAPTED FOR HYBRID VEHICLES; ROAD VEHICLE DRIVE CONTROL SYSTEMS FOR PURPOSES NOT RELATED TO THE CONTROL OF A PARTICULAR SUB-UNIT
- B60W60/00—Drive control systems specially adapted for autonomous road vehicles
- B60W60/001—Planning or execution of driving tasks
- B60W60/0011—Planning or execution of driving tasks involving control alternatives for a single driving scenario, e.g. planning several paths to avoid obstacles
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B60—VEHICLES IN GENERAL
- B60W—CONJOINT CONTROL OF VEHICLE SUB-UNITS OF DIFFERENT TYPE OR DIFFERENT FUNCTION; CONTROL SYSTEMS SPECIALLY ADAPTED FOR HYBRID VEHICLES; ROAD VEHICLE DRIVE CONTROL SYSTEMS FOR PURPOSES NOT RELATED TO THE CONTROL OF A PARTICULAR SUB-UNIT
- B60W30/00—Purposes of road vehicle drive control systems not related to the control of a particular sub-unit, e.g. of systems using conjoint control of vehicle sub-units
- B60W30/08—Active safety systems predicting or avoiding probable or impending collision or attempting to minimise its consequences
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B60—VEHICLES IN GENERAL
- B60W—CONJOINT CONTROL OF VEHICLE SUB-UNITS OF DIFFERENT TYPE OR DIFFERENT FUNCTION; CONTROL SYSTEMS SPECIALLY ADAPTED FOR HYBRID VEHICLES; ROAD VEHICLE DRIVE CONTROL SYSTEMS FOR PURPOSES NOT RELATED TO THE CONTROL OF A PARTICULAR SUB-UNIT
- B60W50/00—Details of control systems for road vehicle drive control not related to the control of a particular sub-unit, e.g. process diagnostic or vehicle driver interfaces
- B60W2050/0001—Details of the control system
- B60W2050/0019—Control system elements or transfer functions
- B60W2050/0028—Mathematical models, e.g. for simulation
-
- Y—GENERAL 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
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
- Y02T10/00—Road transport of goods or passengers
- Y02T10/10—Internal combustion engine [ICE] based vehicles
- Y02T10/40—Engine management systems
Landscapes
- Engineering & Computer Science (AREA)
- Automation & Control Theory (AREA)
- Transportation (AREA)
- Mechanical Engineering (AREA)
- Human Computer Interaction (AREA)
- Control Of Driving Devices And Active Controlling Of Vehicle (AREA)
- Steering Control In Accordance With Driving Conditions (AREA)
- Control Of Position, Course, Altitude, Or Attitude Of Moving Bodies (AREA)
Abstract
The invention relates to the field of automobile collision avoidance control, in particular to a potential field model-based unmanned automobile collision avoidance control method and system; the system comprises a potential field model building module, an MPC collision avoidance path planning module, an MPC control module and a vehicle model module; the potential field model is used for calculating collision avoidance path discrete points, the MPC collision avoidance path planning module is used for optimizing the collision avoidance path discrete points, the MPC control module is used for optimizing and solving the front wheel turning angle of the automobile according to the lateral displacement, the yaw angle and the actual state quantity of the automobile and inputting the front wheel turning angle into the vehicle model module, and the vehicle model module is used for outputting the actual state quantity of the vehicle. The invention adds a path planning on the basis of the traditional tracking control, can optimize the potential field minimum discrete point, plan the path, and then track the path through the tracking control, thereby realizing the obstacle avoidance control under the dynamic traffic environment and improving the smoothness and the same-row efficiency of the path in the obstacle avoidance process.
Description
Technical Field
The invention relates to the field of automobile collision avoidance control, in particular to a potential field model-based unmanned automobile collision avoidance control method and system.
Background
With the continuous increase of automobile holding capacity, the problems of urban congestion, energy consumption, traffic accidents and the like are increasingly prominent, and the automatic driving of automobiles becomes the development trend of the automobile industry. The collision avoidance planning is one of the intelligent automobile active collision avoidance core technologies, so that the burden of a driver is reduced, the commuting efficiency can be effectively improved, and the energy consumption is reduced.
The students at home and abroad have a lot of research results in the aspects of path planning and control methods of autonomous vehicles, and common trajectory planning algorithms comprise a random search method based on trajectory tracking, a trajectory planning method based on a specific function, a model prediction method based on an optimized trajectory, an artificial potential field method and the like. The above-described method is mainly applied to motion planning of autonomous cars and modeling of driver behavior in specific traffic scenarios, such as car following models. The main problems of the existing model are that consideration on risk factors is lacked, complex road conditions and planning and modeling on vehicle-road interaction are not carried out, and the existing evasion model is difficult to adapt to interaction and dynamic change of traffic environment and vehicle state.
Disclosure of Invention
The invention provides a potential field model-based unmanned automobile collision avoidance control method and system, aiming at the problems that the existing obstacle avoidance method does not consider the influence of complex and variable environmental factors and the self state of a vehicle on obstacle avoidance of an automatic driving vehicle, and is difficult to adapt to the interaction and dynamic change of a traffic environment and the state of the vehicle. The method comprises the steps of firstly obtaining obstacle avoidance path discrete points through calculation by using a dynamic potential field model, inputting the obstacle avoidance path discrete points into an MPC controller for optimization, then providing a tire rigidity prediction method based on obstacle avoidance path information, realizing prediction and linearization of tire force in the prediction field by using predicted tire state rigidity, and finally realizing vehicle obstacle avoidance control through the designed MPC obstacle avoidance controller.
In a first aspect of the present invention, the present invention provides a potential field model-based collision avoidance control method for an unmanned vehicle, the method comprising:
establishing a potential field model according to an object vector in a road and the road condition, and calculating collision avoidance path discrete points according to the potential field model;
optimizing the discrete points of the avoidance path, and calculating to obtain the minimum point of a risk field formed by combining the potential energy of each vehicle in the road condition from the discrete points of the avoidance path;
constructing a point quality model for obstacle avoidance path planning based on vehicle dynamics constraint according to the motion vector and the position vector of the vehicle;
constructing a first cost objective function by using the point quality model and taking a system stable point as a control target;
calculating the first generation price target function through an optimization algorithm to obtain an optimal control sequence of the system in any time domain corresponding to each moment;
adding the minimum point of the risk field into a point quality model corresponding to a first cost objective function, and calculating the lateral displacement of an avoidance path and the optimal control variable corresponding to a yaw angle under the target of minimizing the sum of the output state quantity deviation, the control variable and the potential energy of the minimum point of the risk field;
designing a linearized tire model, establishing a monorail vehicle kinematic model, and calculating a state quantity in a prediction time sequence and a control quantity of prediction output from the monorail vehicle kinematic model through a current state quantity of a system and a control increment in a control time domain;
taking the two norms of the deviation between the predicted output state quantity and the actual state quantity as performance indexes of speed control, taking the two norms of control increment as control input smooth indexes, and constructing a second cost objective function;
calculating the second cost objective function by adopting an Active-Set algorithm to obtain the optimal open-loop control increment corresponding to the front wheel corner and the yaw moment of the vehicle;
and inputting the optimal control variable and the first element of the optimal open-loop control increment serving as actual control variables and control increments into a vehicle model to complete path tracking control of the vehicle.
In a second aspect of the present invention, the present invention further provides a collision avoidance control system for an unmanned vehicle based on a potential field model, the system comprising:
the potential field model building module is used for building a potential field model according to object vectors in a road and road conditions and calculating discrete points of a collision avoidance path according to the potential field model;
the MPC collision avoidance path planning module is used for optimizing the discrete points of the collision avoidance path and calculating the minimum point of the risk field formed by combining the potential energy of each vehicle in the road condition from the discrete points of the collision avoidance path; constructing a point quality model for obstacle avoidance path planning based on vehicle dynamics constraint according to the motion vector and the position vector of the vehicle; constructing a first cost objective function by using the point quality model and taking a system stable point as a control target; calculating the first generation price target function through an optimization algorithm to obtain an optimal control sequence of the system in any time domain corresponding to each moment; adding the minimum point of the risk field into a point quality model corresponding to the first cost objective function, and calculating the optimal control variable corresponding to the lateral displacement and the yaw angle of the avoidance path under the target of minimizing the sum of the deviation of the output state quantity, the control variable and the potential energy of the minimum point of the risk field;
the MPC control module is used for designing a linearized tire model, establishing a single-track vehicle kinematic model, and calculating a state quantity in a prediction time sequence and a control quantity of prediction output from the single-track vehicle kinematic model through a current state quantity of a system and a control increment in a control time domain; taking the two norms of the deviation between the predicted output state quantity and the actual state quantity as performance indexes of speed control, taking the two norms of control increment as control input smooth indexes, and constructing a second cost objective function; calculating the second cost objective function by adopting an Active-Set algorithm to obtain the optimal open-loop control increment corresponding to the front wheel corner and the yaw moment of the vehicle;
and the vehicle model module is used for inputting the optimal control variable and the first element of the optimal open-loop control increment as actual control variables and control increments to a vehicle model to complete path tracking control of the vehicle.
The invention has the beneficial effects that:
according to the method, the collision avoidance path points under the dynamic working condition can be predicted by establishing the comprehensive potential field model considering the road environment and the vehicle characteristics, so that the vehicle can autonomously avoid the collision under the dynamic working condition, the collision avoidance problem under the conditions of considering the road environment and the characteristics of the autonomous vehicle can be effectively solved, the collision avoidance control effect under the dynamic working condition can be obviously improved through the MPC controller and the corresponding control method, and the obstacle avoidance efficiency and stability can be improved.
Drawings
FIG. 1 is a flow chart of a collision avoidance control method for an unmanned vehicle according to an embodiment of the present invention;
FIG. 2 is a dynamic potential field map and a static potential field map in an embodiment of the present invention;
FIG. 3 is a schematic view of a kinematic model of a vehicle in an embodiment of the invention;
FIG. 4 is a schematic structural diagram of a collision avoidance control system of an unmanned vehicle according to an embodiment of the present invention;
fig. 5 is a schematic diagram of an obstacle avoidance planning path under a dynamic working condition in the embodiment of the present invention;
FIG. 6 is a diagram illustrating a response of an obstacle avoidance control yaw angle under a dynamic condition in an embodiment of the present invention;
fig. 7 is a response diagram of the obstacle avoidance control sideslip angle under the dynamic condition in the embodiment of the present invention.
Detailed Description
The technical solutions in the embodiments of the present invention will be clearly and completely described below with reference to the drawings in the embodiments of the present invention, and it is obvious that the described embodiments are only a part of the embodiments of the present invention, and not all of the embodiments. All other embodiments, which can be derived by a person skilled in the art from the embodiments given herein without making any creative effort, shall fall within the protection scope of the present invention.
Fig. 1 is a flowchart of a collision avoidance control method for an unmanned vehicle in an embodiment of the present invention, and as shown in fig. 1, the control method includes:
101. establishing a potential field model according to an object vector in a road and the road condition, and calculating collision avoidance path discrete points according to the potential field model;
in the embodiment of the present invention, a potential field model is first established according to the vectors of objects in the road and the condition of the road itself.
In some embodiments, the potential field model may be represented as follows:
wherein, U K_q A potential field model representing an object q; a. The q Representing uncertainty constants, for correction of the model, gamma 1 For calibration constants greater than 1, corrected for different vehicle performance, gamma 2 A correction factor greater than 0, determined from historical traffic accident data for the road segment, for calibrating the effect of speed on driving risk; m is a group of q Representing the equivalent mass, W, of an object q in the road q A road condition factor representing the location of object q; r is q =(x-x q ,y-y q ) Representing an object qFrom the surrounding positions (x, y) to the position (x) of the object q q ,y q ) A distance vector of (d); v. of q Indicating the speed, theta, of a moving object q q Indicating the direction of the object q speed and r q The clockwise direction of the formed included angle is positive; t is c Represents the road curvature factor, ± represents the vehicle turning direction, + represents the vehicle turning direction as left, -represents the vehicle turning direction as right.
In further embodiments, the potential field model may also be represented as follows:
wherein, U K_q A potential field model representing an object q; a. The q Representing an uncertainty constant for correcting the model; m is a group of q Representing the equivalent mass, W, of an object q in the road q A road condition factor representing the location of object q; r is q =(x-x q ,y-y q ) Representing the positions (x, y) around the object q to the position (x) where the object q is located q ,y q ) The distance vector of (a); k is a radical of formula 1 And k 2 Is the gain factor of the potential energy field; v. of q Representing the speed, theta, of a moving object q q Indicating the direction of the object q speed and r q The clockwise direction of the formed included angle is positive; t is c Represents the road curvature factor, ± represents the vehicle turning direction, + represents the vehicle turning direction as left, -represents the vehicle turning direction as right.
For the above model, the road condition factor W q Can be expressed as:
where κ is the visibility coefficient. μ is the adhesion coefficient. I is the road gradient coefficient. In addition, eta 1 And η 2 And correcting the coefficient according to accident data of the road section. Kappa * As a standard value of the visibility coefficient, μ * Is a standard value of road adhesion coefficient. Two of themThis parameter is determined according to the local road conditions and is typically 1.
Equivalent mass M q Can be expressed as:
M q =[16.435·ln(Δv q )-35.472]m q G q H q
wherein m is q Is the mass of object q, G q Is the type of object H q Is the structural shape of the object, Δ v q Is the standard deviation of the vehicle speed from the current road average speed.
Fig. 2 is a dynamic potential field diagram in the embodiment of the present invention, as shown in fig. 2, the center of the potential field is the center of mass of the moving object, and the closer to the front of the vehicle, the more the kinetic energy potential energy is, the more the risk characteristics of actual vehicle driving are met. Through the potential field model, discrete points of the collision avoidance path can be determined, wherein the discrete points are expressed as min (U) K_q ) (ii) a The collision avoidance path discrete points are a series of object discrete points, and the collision avoidance path can be obtained after the discrete points are fitted.
102. Optimizing the discrete points of the avoidance path, and calculating the minimum point of a risk field formed by combining the potential energy of each vehicle in the road condition from the discrete points of the avoidance path;
because the path formed by the discrete points of the avoidance path calculated in the step 101 is not smooth enough and is not suitable for vehicle control, the method needs to optimize the path and obtains the minimum point of the risk field formed by combining the potential energy of each vehicle in the road condition by calculation from the discrete points of the avoidance path;
calculating to obtain a minimum point of a risk field formed by combining potential energy of each vehicle in the road condition, wherein the expression is as follows:
wherein E is the driving risk field U of the position of the vehicle q obs,q I.e. by the smallest n vehicle potential energy fields U K_q And (4) forming.
103. Constructing a point quality model for obstacle avoidance path planning based on vehicle dynamics constraint according to the motion vector and the position vector of the vehicle;
in consideration of vehicle dynamics constraints, a point mass model is employed, whose expression is as follows:
in the formula, xi = [ v ] y ,v x ,φ,X,Y]The total number of the 5 discrete state quantities is 5, and the state quantities are the vehicle speed, the heading angle and the abscissa and the ordinate of the vehicle body in the geodetic coordinate system in turn in the y direction and the x direction of the vehicle. ξ (t) represents the state quantity corresponding to the time t, i.e. in the control input sequence U t A state quantity track acting under the system; u (t) represents a control input amount at time t, and here is a control amount corresponding to a state quantity trajectory composed of a lateral displacement and a yaw angle of the avoidance path. Mu is road adhesion coefficient, g is gravitational acceleration. Superscript denotes derivation.
104. Constructing a first cost objective function by using the point quality model and taking a system stable point as a control target;
in the embodiment of the present invention, the system stabilization point is f (0, 0) =0, and with this point as the system control target, for any time domain N, the optimization target is the first generation price target function J N Expressed as:
wherein, U t =[u t ,...,u t+N-1 ]Is a control input sequence in the time domain N, and u (t) represents the control input quantity at the time t; ξ (t) represents the state quantity corresponding to the time t, i.e. in the control input sequence U t A state quantity track acting under the system; l (-) represents the traceability of the trace of the state quantity of the expected output, and P (-) represents the interruption constraint.
105. Calculating the first generation price target function through an optimization algorithm to obtain an optimal control sequence of the system in any time domain corresponding to each moment;
by optimizing and solving the first generation price target function, the optimal control sequence U at the t moment of the system can be obtained at the t moment of the system t,t =[u t,t ,...u t+N-1,t ]At this time, the first element of the optimal control sequence can be used as the actual input of the controlled object, i.e., u (ξ (t)) = u t,t 。
The optimization algorithm may be any one of the existing algorithms, for example, a convex approximation-based obstacle avoidance principle and an unmanned vehicle path planning model prediction algorithm, which are published in an automated report in 1 month 2020 from korea and the like, may be adopted. 106. Adding the minimum point of the risk field into a point quality model corresponding to a first cost objective function, and calculating the lateral displacement of an avoidance path and the optimal control variable corresponding to a yaw angle under the target of minimizing the sum of the output state quantity deviation, the control variable and the potential energy of the minimum point of the risk field;
taking the point quality model established in the step 104 as a prediction model of an obstacle avoidance path planning module, and adding the potential field minimum point into the point quality model to plan a collision avoidance path, wherein the expression is as follows:
s.tU min ≤U t ≤U max
wherein, U t Representing the corresponding optimal control sequence of the system at the time t; n is a radical of p Is a prediction time domain; e is a driving safety field formed by n obstacles, lambda represents a track point output by the system, and lambda ref Representing the track point of the reference output, wherein (t + i | t) is the prediction of the current time t to the time after i steps; q and R are respectively a weight matrix for adjusting the tracking capability of the reference track point and an adjusting control sequence U i A weight matrix of (a); first itemIndicating the predicted time domain N p Internal systemDeviation of the output from a reference output; second itemThe requirement of the system on the magnitude of the control variable is indicated, the third item E indicates the requirement on the magnitude of the field intensity value of the driving safety field, namely, in the optimization process, the driving risk is ensured to be minimum, and the deviation with the reference track is considered; u shape min And U max Respectively representing the minimum value of the control sequence and the maximum value set of the control sequence in the system.
107. Designing a linearized tire model, establishing a monorail vehicle kinematic model, and calculating a state quantity in a prediction time sequence and a control quantity of prediction output from the monorail vehicle kinematic model through a current state quantity of a system and a control increment in a control time domain;
designing a linear tire model based on the vehicle kinematics model, represented as follows:
1) The state of stiffness of the tire is defined,
defining the tyre state rigidity C as the lateral force F under each slip angle alpha y The ratio of the slip angle α is expressed as follows:
wherein the slip angles alpha of the front and rear tires f And alpha r Are respectively defined as follows:
wherein alpha is f ,α r Is the side slip angle, v, of the front and rear tires x Is the longitudinal speed, v, of the vehicle y Is the lateral speed of the vehicle, omega is the yaw rate of the vehicle, l f ,l r Respectively, the wheelbase and the wheelbase, and delta is the steering angle.
The slip rate s of the tire on the ground can be calculated by the following formula:
wherein, ω is T Is the angular velocity of the wheel, v is the wheel center velocity, and r is the wheel rolling radius.
2) And (3) carrying out linear equation design on the linear tire model, substituting the tire lateral force and the tire slip angle in the nonlinear tire model into the defined tire state rigidity C to obtain the tire rigidity of each tire, wherein based on the obtained tire rigidity, the lateral force of the front tire and the rear tire can be represented in a linear mode as follows:
F x,i =C xi α xi
F y,i =C i α i
wherein, subscripts i = f, r respectively refer to front and rear tires, and xy respectively refer to x-direction and y-direction; i.e. F x,i Represents the lateral force in the x direction received by the front and rear tires i; f y,j Indicating the lateral force in the y-direction experienced by the front and rear tires i.
After the nonlinear tire model and the linear tire model are established, a single-track vehicle kinematic model is established, and the expression of the single-track vehicle kinematic model is as follows:
where m is the vehicle mass, v x ,v y ω is the longitudinal speed, lateral speed and yaw rate of the vehicle, a x And a y Acceleration of the mass center in x-direction and y-direction in the coordinate game of the vehicle, I z For the moment of inertia of the vehicle about the z-axis, F xf ,F xr ,∑M z Longitudinal, transverse and yaw forces respectively, F yf And F yr Lateral forces respectively applied to the front and rear wheels, /) f ,l r Respectively the fore-aft wheelbase. X and Y are positions in a geodetic coordinate system.Is the yaw angle, δ is the steering angle; superscript denotes derivation.
Substituting the linearized representation of the lateral forces of the front and rear tires into the above-described monorail vehicle kinematics model, a predictive model for the MPC controller can be obtained as:
obtaining a vehicle linear kinetic equation by linearizing the predictive model about the operating point of the vehicle:
wherein A (t) and B (t) are Jacobian matrixes of a state equation expression f relative to a state quantity xi and a controlled quantity u respectively, and are selectedSolving the deflection angle of the front wheel for output quantity according to the slip angle and the transverse displacement; then a discrete linear time-varying equation of state can be obtained as:
wherein the content of the first and second substances,
A k,t =I+TA(t),B k,t =TB(t),i is the identity matrix and T is the sampling time. For the Jacobian matrix solution there are:
suppose thatConverting the control quantity u (t) of the discrete linear time-varying state equation into a control increment delta u (t), and performing corresponding transformation on the system state equation to obtain a new spatial state expression as follows:
wherein each matrix is defined as follows:
Δu(k|t)=u(k|t)-u(k-1|t)
wherein n is a state quantity dimension, and m is a control quantity dimension. If the state quantity and the control increment at the time t are acquired from the system, the output quantity of the system at the time t +1 can be predicted by using the formula, and the output quantity of the system at the time k can be acquired through iterative calculation.
Calculating the prediction output, and taking the prediction time domain as N according to the model prediction control theory p Control time domain as N c The state quantity and prediction output in the prediction time domain can be obtained as follows:
and (3) representing the output of the system at the future time in a matrix mode:
in the formula (I), the compound is shown in the specification,
therefore, the state quantity and the output quantity in the prediction time domain can be calculated through the current state quantity of the system and the control increment in the control time domain.
108. Taking the two norms of the deviation between the predicted output state quantity and the actual state quantity as performance indexes of speed control, taking the two norms of control increment as control input smooth indexes, and constructing a second cost objective function;
in order to construct the second generation price target function, firstly some constraints need to be designed:
1) In order to make the tracking process more stable, the invention sets the vehicle limit front wheel deflection angle and the front wheel deflection angle increment as follows:
-25°≤δ≤25°
-0.5°≤Δδ≤0.5°
2) According to the previous researches, on a good dry asphalt pavement, the centroid deviation angle limit of the vehicle for stable running can reach +/-12 degrees, on an ice and snow road, the limit value can reach +/-2 degrees, so the centroid deviation angle constraint condition is set as follows:
beta is more than or equal to 12 degrees and less than or equal to 12 degrees (good road surface)
Beta is more than or equal to-2 degrees and less than or equal to 2 degrees (ice and snow road surface)
3) The longitudinal acceleration and the lateral acceleration are limited by the ground adhesion coefficient, and the following relationship exists:
if the longitudinal speed is not changed, further obtaining:
|a y |≤μg
4) If the road attachment condition is good, the constraint is loose, but the stability of the vehicle body is affected by the overlarge lateral acceleration, the optimal solution is not generated in the solving process due to the undersize, and therefore the constraint is set as a soft constraint condition, namely, each control period controller dynamically adjusts the constraint to the solving condition, and the constraint is as follows:
wherein a is y.min And a y.max Respectively, the minimum and maximum of the lateral acceleration, and epsilon is the relaxation factor.
After the constraint condition is set, the embodiment of the present invention continues to use the two norms of the deviation between the expected trajectory and the actual trajectory as the performance index of the speed control, and use the two norms of the control input increment as the control input smoothing index, and obtain the following second generation price target function:
wherein J represents a second cost objective function; n is a radical of p Is a prediction time domain, N c Is the control time domain; λ (k + i) represents the trace point output by the system at time k + i, λ ref (k) Representing the trace point of the reference output of the system at the moment k; q and R are weighting matrices; Δ u (k + i-1) represents the optimal open-loop control increment corresponding to the front wheel angle and yaw moment of the vehicle at time k + i-1.
Thus, the tracking control problem can be described as:
-25°≤δ≤25°
-0.5°≤Δδ≤0.5°
beta is more than or equal to 12 degrees and less than or equal to 12 degrees (good road surface)
Beta is more than or equal to-2 degrees and less than or equal to 2 degrees (ice and snow road surface)
|a y |≤μg
109. Calculating the second cost objective function by adopting an Active-Set algorithm to obtain the optimal open-loop control increment corresponding to the front wheel corner and the yaw moment of the vehicle;
in the embodiment of the invention, the optimal open-loop control increment corresponding to the front wheel corner and the yaw moment of the vehicle can be obtained by solving the tracking control problem model corresponding to the second cost objective function;
the expression of the optimal open-loop control sequence delta u is as follows:
110. and inputting the optimal control variable and the first element of the optimal open-loop control increment serving as actual control variables and actual control increments into a vehicle model to complete path tracking control of the vehicle.
In the embodiment of the invention, the first element in the control system, namely the first element comprising the optimal control variable and the first element of the optimal open-loop control increment are used as the actual control output increment and input into a Carsim vehicle model to realize the automobile path tracking control. Fig. 4 is a schematic structural diagram of a collision avoidance control system of an unmanned vehicle in an embodiment of the present invention, and as shown in fig. 4, the system structure includes: the potential field model building module is used for building a potential field model according to object vectors in a road and road conditions and calculating discrete points of a collision avoidance path according to the potential field model;
the MPC collision avoidance path planning module is used for optimizing the collision avoidance path discrete points and calculating the minimum point of the risk field formed by combining the potential energy of each vehicle in the road condition from the collision avoidance path discrete points; constructing a point quality model for obstacle avoidance path planning based on vehicle dynamics constraint according to the motion vector and the position vector of the vehicle; constructing a first cost objective function by using the point quality model and taking a system stable point as a control target; calculating the first generation price target function through an optimization algorithm to obtain an optimal control sequence of the system in any time domain corresponding to each moment; adding the minimum point of the risk field into a point quality model corresponding to the first cost objective function, and calculating the optimal control variable corresponding to the lateral displacement and the yaw angle of the avoidance path under the target of minimizing the sum of the deviation of the output state quantity, the control variable and the potential energy of the minimum point of the risk field;
the MPC control module is used for designing a nonlinear tire model and a linearized tire model, establishing a single-track vehicle kinematic model, and calculating state quantity in a prediction time sequence and control quantity of prediction output from the single-track vehicle kinematic model through the current state quantity of a system and control increment in a control time domain; taking the two norms of the deviation between the predicted output state quantity and the actual state quantity as performance indexes of speed control, taking the two norms of control increment as control input smooth indexes, and constructing a second cost objective function; calculating the cost target by adopting an Active-Set algorithm to obtain the optimal open-loop control increment corresponding to the front wheel corner and the yaw moment of the vehicle;
and the vehicle model module is used for inputting the optimal control variable and the first element of the optimal open-loop control increment as actual control variables and control increments to a vehicle model to complete path tracking control of the vehicle. In an embodiment of the invention, the potential field model construction module is arranged to apply the desired lateral displacement Y ref (t) and expectationInputting the data into an MPC collision avoidance path planning module; the MPC collision avoidance path planning module expects the planned result to be the transverse displacement Y ref_local (t) and the desired longitudinal displacement X ref_local (t) input to the MPC control module, which outputs the steering angle δ of the front tire to the vehicle model module f At this time, the vehicle model module may feed back the state quantity ξ in real time to the MPC control module.
In an embodiment, this embodiment specifically describes the method of the present invention with a certain vehicle model and Simulink of Carsim automobile simulation software as a platform, and the main parameters are shown in table 1:
TABLE 1 potential field model parameters and controller parameters
Based on the potential field model-based unmanned automobile collision avoidance control method and system provided by the invention, simulation result graphs shown in fig. 5-7 can be obtained by calculating the potential field model parameters and the controller parameters, fig. 5 shows an obstacle avoidance planning path graph of the invention, fig. 6 shows a change graph of a tracked vehicle yaw angle along with tracking time, and fig. 7 shows a change graph of a tracked vehicle sideslip angle along with tracking time.
In the description of the present invention, it is to be understood that the terms "coaxial", "bottom", "one end", "top", "middle", "other end", "upper", "one side", "top", "inner", "outer", "front", "center", "both ends", and the like are used in the orientations and positional relationships indicated in the drawings, which are for convenience of description and simplicity of description, and do not indicate or imply that the device or element referred to must have a particular orientation, be constructed and operated in a particular orientation, and therefore, are not to be construed as limiting the present invention.
In the present invention, unless otherwise explicitly specified or limited, the terms "mounted," "disposed," "connected," "fixed," "rotated," and the like are to be construed broadly, e.g., as being fixedly connected, detachably connected, or integrated; can be mechanically or electrically connected; the terms may be directly connected or indirectly connected through an intermediate agent, and may be used for communicating the inside of two elements or interacting relation of two elements, unless otherwise specifically defined, and the specific meaning of the terms in the present invention can be understood by those skilled in the art according to specific situations.
Although embodiments of the present invention have been shown and described, it will be appreciated by those skilled in the art that changes, modifications, substitutions and alterations can be made in these embodiments without departing from the principles and spirit of the invention, the scope of which is defined in the appended claims and their equivalents.
Claims (10)
1. A collision avoidance control method of an unmanned automobile based on a potential field model is characterized by comprising the following steps:
establishing a potential field model according to an object vector in a road and the road condition, and calculating collision avoidance path discrete points according to the potential field model;
optimizing the collision avoidance path discrete points, and calculating to obtain a minimum point of a risk field formed by combining the potential energy of each vehicle in the road condition from the collision avoidance path discrete points;
constructing a point quality model for obstacle avoidance path planning based on vehicle dynamics constraint according to the motion vector and the position vector of the vehicle;
constructing a first cost objective function by using the point quality model and taking a system stable point as a control target;
calculating the first price target function through an obstacle avoidance principle based on convex approximation and an unmanned vehicle path planning model prediction algorithm to obtain an optimal control sequence of the system in any time domain corresponding to each moment;
adding the minimum point of the risk field into a point quality model corresponding to a first cost objective function, and calculating the lateral displacement of an avoidance path and the optimal control variable corresponding to a yaw angle under the target of minimizing the sum of the output state quantity deviation, the control variable and the potential energy of the minimum point of the risk field;
designing a linearized tire model, establishing a monorail vehicle kinematic model, and calculating a state quantity in a prediction time sequence and a control quantity of prediction output from the monorail vehicle kinematic model through a current state quantity of a system and a control increment in a control time domain;
taking the two norms of the deviation between the predicted output state quantity and the actual state quantity as performance indexes of speed control, taking the two norms of control increment as control input smooth indexes, and constructing a second cost objective function;
calculating the second cost objective function by adopting an Active-Set algorithm to obtain the optimal open-loop control increment corresponding to the front wheel corner and the yaw moment of the vehicle;
and inputting the optimal control variable and the first element of the optimal open-loop control increment serving as actual control variables and control increments into a vehicle model to complete path tracking control of the vehicle.
2. The potential field model-based collision avoidance control method for the unmanned automobile, according to claim 1, wherein the potential field model comprises:
wherein, U K_q Representing a potential field model; a. The q Representing uncertainty constants, for correction of the model, gamma 1 For calibration constants greater than 1, corrected for different vehicle performance, gamma 2 A correction factor greater than 0, determined from historical traffic accident data for the road segment, for calibrating the effect of speed on driving risk; m q Representing the equivalent mass, W, of an object q in the road q A road condition factor representing the location of object q; r is a radical of hydrogen q =(x-x q ,y-y q ) Indicates the positions (x, y) around the object q to the position (x) of the object q q ,y q ) The distance vector of (a); v. of q Representing the speed, theta, of a moving object q q Representing the direction of the speed of the object q and r q The clockwise direction of the formed included angle is positive; t is q Representing a road curvature factor; + -represents the vehicle turning direction, + represents the vehicle turning direction as left, -represents the vehicle turning direction as right.
3. The potential field model-based collision avoidance control method for the unmanned vehicle according to claim 1, wherein the potential field model further comprises:
wherein, U K_q Representing a potential field model; a. The q Representing an uncertainty constant for correcting the model; m is a group of q Representing the equivalent mass, W, of an object q in the road q A road condition factor representing the location of object q; r is q =(x-x q ,y-y q ) Representing the positions (x, y) around the object q to the position (x) where the object q is located q ,y q ) A distance vector of (d); k is a radical of formula 1 And k 2 A gain factor expressed as a potential energy field; v. of q Representing the speed, theta, of a moving object q q Representing the direction of the speed of the object q and r q The clockwise direction of the formed included angle is positive; t is q Indicates the road curvature factor, + indicates the vehicle cornering direction is left, -indicates the vehicle cornering direction is right.
4. The potential field model-based collision avoidance control method for the unmanned vehicle according to claim 1, wherein the first cost objective function is expressed as:
wherein, J N Representing a first cost objective function; u shape t To indicate any timeControl input sequence, U, within field N t =[u(t),...,u(t+N-1)]U (t) represents a control input amount at time t; ξ (t) represents the state quantity corresponding to time t, i.e. in the control input sequence U t A state quantity track acting under the system; l (-) represents the traceability of the trace of the state quantity of the expected output, and P (-) represents the interruption constraint.
5. The unmanned automobile collision avoidance control method based on the potential field model as claimed in claim 1, wherein the calculating of the optimal control variables corresponding to the lateral displacement and the yaw angle of the avoidance path, that is, the calculating of the optimized path, comprises:
s.t U min ≤U t ≤U max
wherein, U t Representing the optimal control sequence corresponding to the system at the time t; n is a radical of hydrogen p Is a prediction time domain; e is a driving safety field formed by n obstacles, lambda represents a track point output by the system, and lambda is a safety field ref Representing the track point of the reference output, wherein (t + i | t) is the prediction of the current time t to the time after i steps; u shape i The control variable corresponding to the i moment in the prediction time domain; q and R are respectively a weight matrix for adjusting the tracking capability of the reference track point and an adjusting control sequence U i A weight matrix of (a); first itemIndicating the predicted time domain N p Deviation of an internal system output from a reference output; second itemThe requirement of the system on the magnitude of the control variable is shown, the third item | E | represents the requirement on the magnitude of the field intensity value of the driving safety field, namely, in the optimization process, the driving risk is ensured to be minimum, and the deviation from the reference track is considered; u shape min And U max Respectively represent systemAnd the control sequence minimum value and the control sequence maximum value set in the system.
6. The potential field model-based unmanned automobile collision avoidance control method according to claim 5, wherein the minimum point of the risk field is a driving safety field formed by n obstacles, and is represented asWherein E is the risk field U of the location of the vehicle q obs,q I.e. by the smallest n vehicle potential fields U K_q And (4) forming.
7. The potential field model-based unmanned automobile collision avoidance control method according to claim 1, wherein the single-track vehicle kinematic model comprises tire state stiffness defined by a ratio of a lateral force of a tire to a slip angle corresponding to the lateral force; designing a linear tire model based on the tire state stiffness; obtaining a single-track vehicle kinematic model from the non-linear tire model and the linearized tire model, represented as:
where m is the vehicle mass, v x ,v y ω is the longitudinal speed, lateral speed and yaw rate of the vehicle, respectively; a is x And a y The mass center acceleration in the x direction and the y direction in the vehicle coordinate system is obtained; I.C. A z For the moment of inertia of the vehicle about the z-axis, F xf ,F xr ,∑M z The longitudinal force, the transverse force and the yawing moment borne by the vehicle are respectively; f yf And F yr The lateral forces applied to the front wheel and the rear wheel are respectively; l. the f ,l r Respectively the front and rear wheelbases; x and Y are the positions of the vehicles in the geodetic coordinate system;is a yaw angle; delta is a steering angle; superscript denotes derivation.
8. The potential field model-based collision avoidance control method for the unmanned vehicle according to claim 1, wherein the second cost objective function is expressed as:
wherein J represents a second cost objective function; n is a radical of hydrogen p Is a prediction time domain, N c Is the control time domain; λ (k + i) represents the trace point output by the system at time k + i, λ ref (k) Representing the trace point of the reference output of the system at the moment k; q 1 And R 1 Respectively adjusting a weighting matrix of the reference track point tracking capability and a weighting matrix of the control increment; Δ u (k + i-1) represents the optimum open-loop control increment corresponding to the front wheel angle and yaw moment of the vehicle at time k + i-1.
9. The potential field model-based unmanned automobile collision avoidance control method according to claim 8, wherein the optimal open-loop control increment corresponding to the front wheel turning angle and the yaw moment of the vehicle is represented as:
where Δ u represents an optimum open-loop control increment corresponding to a front wheel angle and a yaw moment of the vehicle.
10. An unmanned vehicle collision avoidance control system based on a potential field model, the system comprising:
the potential field model building module is used for building a potential field model according to object vectors in a road and road conditions, and calculating discrete points of a collision avoidance path according to the potential field model;
the MPC collision avoidance path planning module is used for optimizing the collision avoidance path discrete points and calculating the minimum point of a risk field formed by combining the potential energy of each vehicle in the road condition from the collision avoidance path discrete points; constructing a point quality model for obstacle avoidance path planning based on vehicle dynamics constraint according to the motion vector and the position vector of the vehicle; constructing a first cost objective function by using the point quality model and taking a system stable point as a control target; calculating the first price target function through an obstacle avoidance principle based on convex approximation and an unmanned vehicle path planning model prediction algorithm to obtain an optimal control sequence of the system in any time domain corresponding to each moment; adding the minimum point of the risk field into a point quality model corresponding to the first cost objective function, and calculating the optimal control variable corresponding to the lateral displacement and the yaw angle of the avoidance path under the target of minimizing the sum of the deviation of the output state quantity, the control variable and the potential energy of the minimum point of the risk field;
the MPC control module is used for designing a linearized tire model, establishing a single-rail vehicle kinematic model, and calculating a state quantity in a prediction time sequence and a control quantity of prediction output from the single-rail vehicle kinematic model through a current state quantity of a system and a control increment in a control time domain; taking the two norms of the deviation between the predicted output state quantity and the actual state quantity as performance indexes of speed control, taking the two norms of control increment as control input smooth indexes, and constructing a second cost objective function; calculating the second cost target function by adopting an Active-Set algorithm to obtain the optimal open-loop control increment corresponding to the front wheel corner and the yaw moment of the vehicle;
and the vehicle model module is used for inputting the optimal control variable and the first element of the optimal open-loop control increment as actual control variables and control increments to a vehicle model to complete path tracking control of the vehicle.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110829535.XA CN113581201B (en) | 2021-07-22 | 2021-07-22 | Potential field model-based collision avoidance control method and system for unmanned vehicle |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110829535.XA CN113581201B (en) | 2021-07-22 | 2021-07-22 | Potential field model-based collision avoidance control method and system for unmanned vehicle |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113581201A CN113581201A (en) | 2021-11-02 |
CN113581201B true CN113581201B (en) | 2022-11-04 |
Family
ID=78248744
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110829535.XA Active CN113581201B (en) | 2021-07-22 | 2021-07-22 | Potential field model-based collision avoidance control method and system for unmanned vehicle |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113581201B (en) |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114312848B (en) * | 2022-01-17 | 2024-04-12 | 吉林大学 | Intelligent driving automobile track planning and tracking control method based on double-layer MPC |
JP7400912B1 (en) | 2022-09-26 | 2023-12-19 | いすゞ自動車株式会社 | automatic driving device |
JP7416301B1 (en) | 2023-03-10 | 2024-01-17 | いすゞ自動車株式会社 | Dynamic route planning device and dynamic route planning method |
CN116494993A (en) * | 2023-04-04 | 2023-07-28 | 同济大学 | Intelligent driving real-time track planning method considering high-precision vehicle dynamic characteristics |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107885932A (en) * | 2017-11-07 | 2018-04-06 | 长春工业大学 | It is a kind of to consider man-machine harmonious automobile emergency collision avoidance layer-stepping control method |
Family Cites Families (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US9199668B2 (en) * | 2013-10-28 | 2015-12-01 | GM Global Technology Operations LLC | Path planning for evasive steering maneuver employing a virtual potential field technique |
US10725470B2 (en) * | 2017-06-13 | 2020-07-28 | GM Global Technology Operations LLC | Autonomous vehicle driving systems and methods for critical conditions |
CN107561942B (en) * | 2017-09-12 | 2021-02-09 | 重庆邮电大学 | Intelligent vehicle trajectory tracking model prediction control method based on model compensation |
CN110217227B (en) * | 2019-06-25 | 2020-09-18 | 长春工业大学 | Steering and braking combined collision avoidance control method suitable for ice and snow road working conditions |
CN111873992A (en) * | 2020-08-11 | 2020-11-03 | 北京理工大学重庆创新中心 | Artificial potential field method for automatic driving vehicle decision layer path planning |
-
2021
- 2021-07-22 CN CN202110829535.XA patent/CN113581201B/en active Active
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107885932A (en) * | 2017-11-07 | 2018-04-06 | 长春工业大学 | It is a kind of to consider man-machine harmonious automobile emergency collision avoidance layer-stepping control method |
Also Published As
Publication number | Publication date |
---|---|
CN113581201A (en) | 2021-11-02 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN113581201B (en) | Potential field model-based collision avoidance control method and system for unmanned vehicle | |
CN111845774B (en) | Automatic driving automobile dynamic trajectory planning and tracking method based on transverse and longitudinal coordination | |
CN108227491B (en) | Intelligent vehicle track tracking control method based on sliding mode neural network | |
Chen et al. | Path tracking and handling stability control strategy with collision avoidance for the autonomous vehicle under extreme conditions | |
CN113320542B (en) | Tracking control method for automatic driving vehicle | |
Cai et al. | Implementation and development of a trajectory tracking control system for intelligent vehicle | |
CN107561942A (en) | Intelligent vehicle track following model predictive control method based on model compensation | |
CN111007722B (en) | Transverse robust fault-tolerant control system and method for four-wheel steering automatic driving automobile | |
Xu et al. | Model predictive control for lane keeping system in autonomous vehicle | |
Attia et al. | Coupled longitudinal and lateral control strategy improving lateral stability for autonomous vehicle | |
CN114379583B (en) | Automatic driving vehicle track tracking system and method based on neural network dynamics model | |
EL HAJJAMI et al. | Neural network based sliding mode lateral control for autonomous vehicle | |
CN111959500A (en) | Automobile path tracking performance improving method based on tire force distribution | |
CN114572251A (en) | High-speed automatic driving automobile track tracking method based on predictive control | |
CN112578672A (en) | Unmanned vehicle trajectory control system based on chassis nonlinearity and trajectory control method thereof | |
CN114637292A (en) | Vehicle track tracking robust control method and system considering obstacle avoidance | |
CN112606843A (en) | Intelligent vehicle path tracking control method based on Lyapunov-MPC technology | |
CN116048081A (en) | Automatic driving vehicle decision and regulation method considering safety boundary constraint | |
Ko et al. | Integrated path planning and tracking control of autonomous vehicle for collision avoidance based on model predictive control and potential field | |
Peng et al. | Lane-change model and tracking control for autonomous vehicles on curved highway sections in rainy weather | |
Li et al. | Path planning and obstacle avoidance control for autonomous multi-axis distributed vehicle based on dynamic constraints | |
CN112596516B (en) | Dubins curve-based multi-vehicle formation switching method | |
Chen et al. | An overtaking obstacle algorithm for autonomous driving based on dynamic trajectory planning | |
Lee et al. | Robust control based on DOBC for improving lateral string stability of truck platooning | |
CN114044003A (en) | Tracking control method for front-rear double-shaft steering vehicle |
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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
TR01 | Transfer of patent right | ||
TR01 | Transfer of patent right |
Effective date of registration: 20240316 Address after: 401320 2-1, Unit 1, Building 20, District D, the Pearl River River City, Huaxi Industrial Park, Banan District, Chongqing Patentee after: Chongqing Zhiche Qinxing Technology Co.,Ltd. Country or region after: China Address before: 400065 Chongwen Road, Nanshan Street, Nanan District, Chongqing Patentee before: CHONGQING University OF POSTS AND TELECOMMUNICATIONS Country or region before: China |