CN107791248A - Control method based on the six degree of freedom serial manipulator for being unsatisfactory for pipper criterions - Google Patents

Control method based on the six degree of freedom serial manipulator for being unsatisfactory for pipper criterions Download PDF

Info

Publication number
CN107791248A
CN107791248A CN201710902038.1A CN201710902038A CN107791248A CN 107791248 A CN107791248 A CN 107791248A CN 201710902038 A CN201710902038 A CN 201710902038A CN 107791248 A CN107791248 A CN 107791248A
Authority
CN
China
Prior art keywords
msub
mtd
mrow
theta
mtr
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
CN201710902038.1A
Other languages
Chinese (zh)
Other versions
CN107791248B (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 University of Technology ZJUT
Original Assignee
Zhejiang University of Technology ZJUT
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 University of Technology ZJUT filed Critical Zhejiang University of Technology ZJUT
Priority to CN201710902038.1A priority Critical patent/CN107791248B/en
Publication of CN107791248A publication Critical patent/CN107791248A/en
Application granted granted Critical
Publication of CN107791248B publication Critical patent/CN107791248B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • BPERFORMING OPERATIONS; TRANSPORTING
    • B25HAND TOOLS; PORTABLE POWER-DRIVEN TOOLS; MANIPULATORS
    • B25JMANIPULATORS; CHAMBERS PROVIDED WITH MANIPULATION DEVICES
    • B25J9/00Programme-controlled manipulators
    • B25J9/16Programme controls
    • B25J9/1602Programme controls characterised by the control system, structure, architecture
    • B25J9/1607Calculation of inertia, jacobian matrixes and inverses
    • GPHYSICS
    • G05CONTROLLING; REGULATING
    • G05BCONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
    • G05B2219/00Program-control systems
    • G05B2219/30Nc systems
    • G05B2219/39Robotics, robotics to robotics hand
    • G05B2219/39073Solve inverse kinematics by fuzzy algorithm

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Mathematical Physics (AREA)
  • Automation & Control Theory (AREA)
  • Robotics (AREA)
  • Mechanical Engineering (AREA)
  • Pharmaceuticals Containing Other Organic And Inorganic Compounds (AREA)
  • Manipulator (AREA)
  • Numerical Control (AREA)

Abstract

The invention discloses the control method based on the six degree of freedom serial manipulator for being unsatisfactory for pipper criterions, belongs to motion planning and robot control method and technology field.The motion control method of six axle six degree of freedom serial manipulators of prior art, it is only applicable to the six axle six degree of freedom serial manipulators for meeting pipper criterions, not high to the control accuracy that is unsatisfactory for the axle six degree of freedom serial manipulator of pipper criterions six, control process is cumbersome.The present invention proposes a kind of parsing control method of new foundation in identical joint coordinate system, when establishing coordinate system to six axle six degree of freedom serial manipulators, all joint coordinate systems are consistent with robot basis coordinates system coordinate direction, and the six axle six degree of freedom serial manipulators for being unsatisfactory for pipper criterions are handled by the way of " standardization " and " transformation approach ", control process is simply orderly, amount of calculation can effectively be reduced, and the six axle six degree of freedom serial manipulators suitable for being unsatisfactory for pipper criterions, control accuracy are high.

Description

Control method based on the six degree of freedom serial manipulator for being unsatisfactory for pipper criterions
Technical field
The present invention relates to the control method based on the six degree of freedom serial manipulator for being unsatisfactory for pipper criterions, belong to machine People's motion control method technical field.
Background technology
Conventional robot kinematics' method for solving mainly has 3 kinds:DH parametric methods, geometric method and analytic method.DH parametric methods It is a kind of wide variety of method, its foundation to joint coordinate system requires very harsh, it is necessary to performed according to specified principle, And cumbersome matrix operation, be unfavorable for embody algorithm real-time.Geometric method typically solution is more difficult, and versatility is bad, only Suitable for the less mechanical mechanism of the free degree.Conventional analytic method is adapted to robot that is simple in construction, being easy to analysis, but computing The solution of arcsine and anticosine often occurs in journey, more solutions just occur in such Robotic inverse kinematics when calculating.
But traditional parsing control method is only applicable to the six axle six degree of freedom serial machines for meeting pipper criterions People-i.e. its be used for the joints axes in last three joints of control machine people's terminal angle and intersect at a point.If these three are closed The joints axes of section are not if intersecting at a point, then the shifted relative on some position is just certainly existed between them. Also, the position and posture of this biasing can also change and change with the motion state in last three joints.Because the biasing is grown The location status of degree has certain uncertainty, and its position and posture to robot end also has significant impact. So we can not be connected by the parsing control method of routine to this six axle six degree of freedoms for being unsatisfactory for pipper criterions Robot carries out motion control, causes that the robot control process to the type is cumbersome, and computationally intensive, control accuracy is not high.
The content of the invention
The defects of for prior art, it can accurately be controlled it is an object of the invention to provide one kind and be unsatisfactory for pipper standards The control process of six degree of freedom serial manipulator motion then is simple, the control method of the few robot of amount of calculation.
To achieve the above object, the technical scheme is that:
Based on the control method for the six degree of freedom serial manipulator for being unsatisfactory for pipper criterions, comprise the following steps:
Step 1:The six axle six degree of freedom serial manipulators for being unsatisfactory for pipper criterions are changed on an equal basis, make its change Change to and meet on six axle six degree of freedom serial manipulator models of pipper criterions;Step 2:To the six axles six degree of freedom When serial manipulator carries out Forward Kinematics Analysis, the direction vector parameter of biasing is reserved;Step 3:According to trying to achieve before Biasing direction vector parameter, to the six axles six degree of freedom serial manipulator carry out changes in coordinates, transformed to satisfaction Six axle six degree of freedom serial manipulators of pipper criterions;Step 4:Analysis of Inverse Kinematics and meter are carried out to it based on analytic method Calculate;Step 5:The problem of progress during Analysis of Inverse Kinematics is carried out be present to the six axles six degree of freedom serial manipulator Related supplemental and explanation;Step 6:The kinematic parameter solved, motion control card is dealt into, to being unsatisfactory for pipper criterions six Axle six degree of freedom serial manipulator is controlled.The present invention proposes that a kind of new foundation controls in the parsing of identical joint coordinate system Method, control process is simply orderly, can effectively reduce amount of calculation, and suitable for being unsatisfactory for six axles six of pipper criterions certainly It is high by degree serial manipulator, control accuracy.
As optimization technique measure, step 1:The six axles six degree of freedom serial manipulator is changed on an equal basis, makes it It is transformed into the six axle six degree of freedom serial manipulator models for meeting pipper criterions.
Although the type serial manipulator is unsatisfactory for pieper criterions, by after corresponding deformation, can still pass through DH parameter models calculate, and the connect end effector arm of mechanical robot of six axles have been carried out into conversion to a certain extent so that end The extended line of last three cradle head axis is held to intersect at a point.
As optimization technique measure, step 2:Positive kinematics are being carried out to the axle six degree of freedom serial manipulator of the type six When analysis, the direction vector parameter of biasing is reserved.
Robot Forward Kinematics Analysis is exactly known rod member geometric parameter and joint angle vector, seeks end effector of robot Position and posture relative to fixed reference coordinate.
Before robot Forward Kinematics Analysis is carried out, it is desirable to establish change of the coordinate system { i } relative to coordinate system { i-1 } Change.This general conversion is made up of four parameters.To any given robot, this conversion is an only variable Function, the other three parameter are determined by mechanical system.Ginseng is converted between the connecting rod of the six axles six degree of freedom serial manipulator It is as shown in the table for number.
Parameter list is converted between the connecting rod of the axle six degree of freedom serial manipulator of table six
Wherein:
aiFor along the X under the i-th coordinate systemiAxle, the Z under the i-th coordinate systemiDirection of principal axis is moved to Z under i+1 coordinate systemi+1 The distance of axle;
αiFor the X under the i-th coordinate systemiAxle, the Z under the i-th coordinate systemiThe Z that axle is rotated under i+1 coordinate systemi+1Axle Angle;
diFor along the Z under the i-th coordinate systemiAxle, the X under the i-th -1 coordinate systemi-1The X that axle is moved under the i-th coordinate systemiAxle Distance;
θiFor the Z under the i-th coordinate systemiAxle, the X under the i-th -1 coordinate systemi-1The X that axle is rotated under the i-th coordinate systemiAxle Angle.
With homogeneous coordinate transformation matrix come represent the normal solution computing of robot (Represent it must is it from i coordinates to j coordinate systems Between 4 × 4 transfer matrix;In formula:By s θiWhat is represented is sin θi;cθiThat represent is cos θ1.And i=1,2,3...6):
Wherein:
θ2323
nx=c θ12345+sθ145+cθ1235
ny=s θ1235+sθ12345-cθ145
nz=c θ235-sθ2345
ox=s θ14-cθ1234
oy=-c θ14-sθ1234
oz=s θ234
ax=c θ12345-cθ1235+sθ1sθ4Sθ5
ay=s θ12345-sθ1235-cθ145
az=-s θ2345-cθ235
px=-d4123+l212+d31-d21+d6ax+l6ox
py=-d4123+l212+d21-d31+d6ay+l6oy
pz=-d423-l22+d6az+l6oz
It is worth noting that:Here the o out calculatedx、oy、ozWith the o hereinafter quotedx5、oy5、oz5Parameter values are Consistent.By θ6=0 brings into above-mentioned equation, it is seen that:The direction vector at the exactly bigoted place of these three parameters description.
As optimization technique measure, step 3:According to the direction vector parameter for the biasing tried to achieve before, to six axle six Free degree serial manipulator carries out coordinate transform, is transformed to the six axle six degree of freedom serial manipulators for meeting pipper criterions On.
Analysis of Inverse Kinematics is the geometric parameter of known machine people's rod member, and given end effector is sat relative to fixed reference Mark the desired locations and posture (wherein, p of systemx′、py′、pz' represent be end effector coordinate position, nx、ny、nz、ox、 oy、oz、ax、ay、azWhat is represented is the coordinate posture of end effector), can analysis robot make its end effector reach this Individual expected pose.
For Analysis of Inverse Kinematics, the coordinate of ending coordinates system posture is known, that is to say, thatParameters It is known.But due to θ6Uncertainty, can not pass through, it is known that byTo drawTransformation matrix. Herein, " standardization " can be taken such to solve the problems, such as.
Above-mentioned Forward Kinematics Analysis is solved what is comeThe o drawnx5、oy5、oz5It is added in whole Analysis of Inverse Kinematics Come, and given data is pre-processed:
Wherein:px′、py′、pz' it is known parameters, by the p after changex、py、pzBring into matrix.
Why given data is pre-processed, be because the serial manipulator is unsatisfactory for pieper criterions.With biography The multiaxis serial manipulator of system is compared, and the joints axes in three joints in end of the type multiaxis serial manipulator are prolonged Long line fails to intersect at a point.One relative is had more so between them on its position of some on end effector Lateral offset.Also, the position and posture of this biasing can also change and change with the motion state in last three joints.Due to The location status of the offset lengths has certain uncertainty, and its position and posture to robot end also has significantly Influence.Only take to its ending coordinates parameter px′、py′、pz' pretreated method is carried out, it could eliminate by above-mentioned inclined Put and do not know to the influence caused by Solve problems, derived so as to carry out further Analysis of Inverse Kinematics and process.In fact, This processing mode be by it is original, be unsatisfactory for the series connection that pieper criterion serial manipulators are converted into meeting pieper criterions Robot.
As optimization technique measure, step 4:Analysis of Inverse Kinematics and calculating are carried out to it based on analytic method.
By formulaIt can obtain:
It is equal by the element (2,4) on above-mentioned equation both sides for first axis joint angle, obtain:
In formula:
Pay attention to:Sign, θ in above formula be present1There can be multiple solutions.Now, θ1It is known.
It is equal by the element (Isosorbide-5-Nitrae) on above-mentioned equation both sides, (3,4) for the 3rd joint angles, obtain:
Wherein:
It is equal by the element (Isosorbide-5-Nitrae) (2,4) on above-mentioned equation both sides for second joint angles, obtain:
So, by θ2323It can obtain:
According to θ1And θ3Four kinds of solution may combine, and can obtain corresponding four kinds of probable values by above formula, then can obtain θ2 Four kinds may solution.
For the 4th joint angles, by the element (1,3) on above-mentioned equation both sides, (3,3) are equal obtains:
θ4(the s θ of=A tan 21ax-cθ1ay, c θ123ax+sθ123ay-sθ23az)
As s θ5When=0, robot is in Singularity, meets degenerative conditions.Now, joint shaft 4 and 6 overlaps, and can only solve Go out θ6With θ4And with difference.In Singularity, θ can be arbitrarily chosen4Angle, calculating corresponding θ6Value.But this feelings Condition is existed only among the serial manipulator for meeting Pepeer criterions.This problem will be begged in detail in next operation By.
For the 5th axis joint angle, by the element (1,3) on above-mentioned equation both sides, (3,3) are equal obtains:
θ5(± the M of=A tan 21, M2)
Wherein:
M1=(c θ1234+sθ14)ax+(sθ1234-cθ14)ay-sθ234az
M2=-c θ123ax-sθ123ay-cθ23az
For the 6th joint angles, it can obtain:
Wherein:
N=c θ12345+sθ145+cθ1235
O=s θ14-cθ1234
As optimization technique measure, step 5:1) to the analysis of generally robot degenerate problem
As s θ5When=0, robot is in Singularity, meets degenerative conditions.Now, joint shaft 4 and 6 overlaps, and can only solve Go out θ6With θ4And with difference.In Singularity, θ can be arbitrarily chosen4Angle, calculating corresponding θ6Value.
But such case is existed only among the serial manipulator for meeting Pepeer criterions, in the six axles six degree of freedom It is that degenerate problem is not present in serial manipulator.Because under the type serial manipulator, the robot being connected with each other is not Conllinear problem can occur.
When Robotic inverse kinematics are analyzed, first the given data of Analysis of Inverse Kinematics is pre-processed.And And this processing mode be exactly be by it is original, be unsatisfactory for pieper criterion serial manipulators and be converted into meeting pieper criterions Serial manipulator.Therefore, during actually solving, or exist and work as s θ5=0 robot meets the situation of degenerative conditions. Below, this degenerate case is carried out to detailed analysis.
Based on joint shaft 4 for coordinate, by s θ5=0, c θ5=± 1 tries to achieveTransition matrix:
Make (1,2), (3,2) are equal obtains:
-ox1+oy1=c θ46±sθ46
θ4±θ6(the o of=A tan 2x1-oy1,-ox123-oy123+oz23)
By previously mentioned ox5、oy5、oz5Parameter, θ can be obtained6
θ6(the o of=A tan 2x, ox5)
It can to sum up obtain:
Wherein:The selection of the positive negativity of above formula is by c θ5Positive and negative judge.
2) analysis that special circumstances are under degenerate state
Above the Singularity of the six axles six degree of freedom serial manipulator is generally analyzed.It is actual On, when the robot motion is to a certain particular point, unusual situation also occurs.Due to such case only exist and certain The state of a kind of location point.Therefore it is referred to as Singularity in particular cases.
Situation 1:Work as d2=d3, and the axis of joint shaft 1 and 4 (or can not solve θ when coincidence, under this state1 And θ4.Such case there is more than among the robot model after change, in original robot model and exist strange Dystopy shape problem.
Such case can occur in both cases:First, such location point is the transit point of movement locus.Second, such Location point is movement locus starting point.
If such location point is the transit point of movement locus, can solve this problem using " inheritance act ". That is, when there is such case, θ1Value acquiescence inherit θ under neighbouring position dotted state1Value, it is and then right θ4Solved.
If such location point is the starting point of movement locus, can solve this problem using " Reconstruction Method ". That is, when there is such case, θ1Value be defaulted as 0, so as to afterwards to θ4Solved.
Situation 2:Work as d2=d3, and when the axis coincidence of joint shaft 1 and 6, θ can not be solved under this state1And θ6。 Such case is existed only among the robot model after conversion, in original robot model is asked in the absence of Singularity Topic.But in above-mentioned Analysis of Inverse Kinematics, the state of this kind of Singularity still occurs.
In this case, can solve this problem using " sciagraphy ".There is such case that is, working as When, θ can be obtained by the projection that robot is formed in xoy planes1Numerical value.
Compared with prior art, the invention has the advantages that:
The present invention is according to the six axle six degree of freedom serial manipulators control feature for being unsatisfactory for pipper criterions, to six axles six When free degree serial manipulator establishes coordinate system, all joint coordinate systems are sat with robot basis coordinates system (global coordinate system) It is consistent to mark direction, and come the six axle six degree of freedom strings to being unsatisfactory for pipper criterions by the way of " standardization " and " transformation approach " Connection robot is handled;Bivariate tan is used when solving joint angles, so passes through the symbol can of independent variable The quadrant where joint angles is determined, avoids the possibility of Lou solution;Occurs the situation of multigroup solution for some joints, using " just Approximately principle ", a solution nearest from current angular is chosen, using the method rejected step by step, is greatly reduced because certain joint occurs The huge amount of calculation that more solutions are brought.
Solution hardly possible, redundancy for being unsatisfactory for the appearance of the axle six degree of freedom serial manipulator kinematics solution of pipper criterions six The problems such as solution, propose that a kind of new foundation in the Analytic Method method of identical joint coordinate system, overcomes six axle six degree of freedoms Serial manipulator is because mass motion caused by itself being unsatisfactory for pipper criterions (inverse solution) solves failure the shortcomings that.In redundancy During solution selection, the method rejected using nearby principle and step by step greatly reduces amount of calculation, raising arithmetic speed, while also avoid Driving phenomenon caused by choosing solution is improper, and then can accurately control and be unsatisfactory for the axle six degree of freedom cascade machine of pipper criterions six Device people.
The present invention proposes a kind of parsing control method of new foundation in identical joint coordinate system, and control process simply has Sequence, amount of calculation can be effectively reduced, and suitable for being unsatisfactory for six axle six degree of freedom serial manipulators of pipper criterions, control Precision is high.
Brief description of the drawings
Fig. 1 is the mechanical arm, three-D model schematic of such six degree of freedom series connection of the present invention;
Fig. 2 is such six axles six degree of freedom serial manipulator machinery local deformation of the invention and establishment of coordinate system schematic diagram;
Fig. 3 is the robot that the present invention is under degenerate state to such six axles six degree of freedom serial manipulator ordinary circumstance Schematic diagram;
Fig. 4 is in Singularity robot model's schematic diagram 1 time for the present invention to special circumstances;
Fig. 5 is in Singularity robot model's schematic diagram 2 times for the present invention to special circumstances.
Embodiment
In order to make the purpose , technical scheme and advantage of the present invention be clearer, it is right below in conjunction with drawings and Examples The present invention is further elaborated.It should be appreciated that specific embodiment described herein is only to explain the present invention, not For limiting the present invention.
On the contrary, the present invention covers any replacement done in the spirit and scope of the present invention being defined by the claims, repaiied Change, equivalent method and scheme.Further, in order that the public has a better understanding to the present invention, below to the thin of the present invention It is detailed to describe some specific detail sections in section description.Part without these details for a person skilled in the art Description can also understand the present invention completely.
The threedimensional model of six degree of freedom series connection mechanical arm includes pedestal, large arm, forearm and the part of end palm four composition, such as Shown in Fig. 1.Large arm can be around pedestal axial-rotation in model, and forearm can be around forearm axial direction around large arm axial-rotation, end's platform Rotation.
As shown in left in Figure 2.Coordinate system { 0 } is located at pedestal lower surface position, and coordinate system { 1 } is located at pedestal and connected with large arm Place centre position is connect, coordinate system { 2 } is located at the centre position of large arm rear end, and coordinate system { 3 } is located at the centre position of forearm rear end, Coordinate system { 4 } and { 5 } are respectively positioned on the centre position of forearm front end, and coordinate system { 6* } is located at end's platform surface parallel and coordinate It is the focal position on { 4 } z-axis extended line, coordinate system { 6 } is located at the surface location in end's platform.D in Fig. 21For coordinate system { 0 } distance between coordinate system { 1 } along the z-axis direction, d2For coordinate system { 1 } and the distance of coordinate system { 2 } along the z-axis direction, d3To sit The distance that mark system { 2 } moves with coordinate system { 3 } along z3 direction of principal axis, d4Moved for coordinate system { 5 } and coordinate system { 6 } along z5 direction of principal axis Dynamic distance;l1For big arm lengths, l2For forearm lengths, l3For end's platform length.
The present invention relates to the control method based on the six degree of freedom serial manipulator for being unsatisfactory for pipper criterions, described side Method comprises the following steps:
Step 1:The axle six degree of freedom serial manipulator of the type six is changed on an equal basis, it is transformed into satisfaction On six axle six degree of freedom serial manipulator models of pipper criterions.
This six axles six degree of freedom serial manipulator structure still has closing solution.Although the serial manipulator is unsatisfactory for Pieper criterions, but by after corresponding deformation, can still be calculated by DH parameter models.As shown in Fig. 2 by left figure The end effector arm of six axles series connection mechanical robot has carried out conversion to a certain extent, and result of variations is as shown in Fig. 2 right figures.Make The rotary shaft for obtaining last three mechanical arms in end intersects at a point.6* coordinate systems are can be passed through to translate by 6 coordinate systems in right figure Conversion obtains.
Step 2:When Forward Kinematics Analysis is carried out to the axle six degree of freedom serial manipulator of the type six, reserve The direction vector parameter of biasing.
Before robot Forward Kinematics Analysis is carried out, it is desirable to establish change of the coordinate system { i } relative to coordinate system { i-1 } Change.This general conversion is made up of four parameters.To any given robot, this conversion is an only variable Function, the other three parameter are determined by mechanical system.As shown in above-mentioned Fig. 2, the six axles six degree of freedom serial manipulator Connecting rod between transformation parameter it is as shown in the table.
Wherein:
aiFor along the X under the i-th coordinate systemiAxle, the Z under the i-th coordinate systemiDirection of principal axis is moved to Z under i+1 coordinate systemi+1 The distance of axle;
αiFor the X under the i-th coordinate systemiAxle, the Z under the i-th coordinate systemiThe Z that axle is rotated under i+1 coordinate systemi+1Axle Angle;
diFor along the Z under the i-th coordinate systemiAxle, the X under the i-th -1 coordinate systemi-1The X that axle is moved under the i-th coordinate systemiAxle Distance;
θiFor the Z under the i-th coordinate systemiAxle, the X under the i-th -1 coordinate systemi-1The X that axle is rotated under the i-th coordinate systemiAxle Angle.
With homogeneous coordinate transformation matrix come represent the normal solution computing of robot (Represent it must is it from i coordinates to j coordinate systems Between 4 × 4 transfer matrix;In formula:By s θiWhat is represented is sin θi;cθiThat represent is cos θ1.And i=1,2,3...6):
Wherein:
θ2323
nx=c θ12345+sθ145+cθ1235
ny=s θ1235+sθ12345-cθ145
nz=c θ235-Sθ2345
ox=S θ14-cθ1234
oy=-c θ14-sθ1234
oz=S θ234
ax=c θ12345-cθ1235+sθ145
ay=s θ12345-sθ1235-cθ145
az=-s θ2345-cθ235
px=-d4123+l212+d31-d21+d6ax+l6ox
py=-d4123+l212+d21-d31+d6ay+l6oy
pz=-d423-l22+d6az+l6oz+h
It is worth noting that:Here the o out calculatedx、oy、ozWith the o hereinafter quotedx5、oy5、oz5Parameter values are Consistent.By 06=0 brings into above-mentioned equation, it is seen that:The vector side at the exactly bigoted place of these three parameters description To.
Step 3:According to the direction vector parameter for the biasing tried to achieve before, to the axle six degree of freedom serial machine of the type six People carries out changes in coordinates, is transformed to the six axle six degree of freedom serial manipulators for meeting pipper criterions.
Analysis of Inverse Kinematics is the geometric parameter of known machine people's rod member, and given end effector is sat relative to fixed reference Mark the desired locations and posture (wherein, p of systemx’、py’、pz' represent be end effector coordinate position, nx、ny、nz、ox、 oy、oz、ax、ay、azWhat is represented is the coordinate posture of end effector), can analysis robot make its end effector reach this Individual expected pose.
For Analysis of Inverse Kinematics, the coordinate of ending coordinates system posture is known, that is to say, thatParameters It is known.But due to θ6Uncertainty, can not pass through, it is known that byTo drawTransformation matrix. Here it is possible to take " standardization " such to solve the problems, such as.
Above-mentioned Forward Kinematics Analysis is solved what is comeThe o drawnx5、oy5、oz5It is added in whole Analysis of Inverse Kinematics Come, and given data is pre-processed:
Wherein:px′、py′、pz' it is known parameters, by the p after changex、py、pzBring into matrix.
Why given data is pre-processed, be because the serial manipulator is unsatisfactory for pieper criterions.With biography The multiaxis serial manipulator of system is compared, and the joints axes in three joints in end of the type multiaxis serial manipulator are prolonged Long line fails to intersect at a point.One relative is had more so between them on its position of some on end effector Lateral offset.Also, the position and posture of this biasing can also change and change with the motion state in last three joints.Due to The location status of the offset lengths has certain uncertainty, and its position and posture to robot end also has significantly Influence.Only take to its ending coordinates parameter px′、py′、pz' pretreated method is carried out, it could eliminate by above-mentioned inclined Put and do not know to the influence caused by Solve problems, so as to carry out further inverse kinematics derivation and analysis.It is in fact, this Processing mode be by it is original, be unsatisfactory for the serial machine that pieper criterion serial manipulators are converted into meeting pieper criterions People, as shown in Figure 2.The serial manipulator model transformed to after pretreatment is basic herein as shown in the right figure in Fig. 2 The coordinate system of upper foundation is as shown in the right figure in Fig. 2.
Step 4:Analysis of Inverse Kinematics and calculating are carried out to it based on analytic method.
By formulaIt can obtain:
It is equal by the element (2,4) on above-mentioned equation both sides for first joint angles, obtain:
1px-cθ1py=d2-d3
To sum up, θ1Solution can be write as:
Pay attention to:Sign, θ in above formula be present1There can be multiple solutions.Now, θ1It is known.
It is equal by the element (Isosorbide-5-Nitrae) on above-mentioned equation both sides, (3,4) for the 3rd axis joint angle, obtain:
px1+py1=-d423+l22
pz=-d423-l22
To sum up, θ3Solution can be write as:
Wherein:
It is equal by the element (Isosorbide-5-Nitrae) (2,4) on above-mentioned equation both sides for second axis joint angle, obtain:
123px+sθ123py-sθ23pz-l23=0
-cθ123px-sθ123py-cθ23pz+l23=d4
To sum up, θ2Solution can be write as:
So, by θ2323It can obtain:
According to θ1And θ3Four kinds of solution may combine, and can obtain corresponding four kinds of probable values by above formula, then can obtain θ2 Four kinds may solution.
For the 4th joint angles, by the element (1,3) on above-mentioned equation both sides, (3,3) are equal obtains:
123ax+sθ123ay-sθ23az=c θ45
-sθ1ax+cθ1ay=-s θ45
:θ4(the s θ of=A tan 21ax-cθ1ay, c θ123ax+sθ123ay-sθ23az)
As s θ5When=0, robot is in Singularity, meets degenerative conditions.Now, joint shaft 4 and 6 overlaps, and can only solve Go out θ6With θ4And with difference.In Singularity, θ can be arbitrarily chosen4Angle, calculating corresponding θ6Value.But this feelings Condition is existed only among the serial manipulator for meeting Pepeer criterions.This problem will be begged in detail in next operation By.
For the 5th, the 6th joint angles, by the element (1,3) on above-mentioned equation both sides, (3,3) are equal obtains:
(cθ1234+sθ14)ax+(sθ1234-cθ14)ay-sθ234az=s θ5
-cθ123ax-sθ123ay-cθ23az=c θ5
θ5(± the M of=A tan 21, M2)
Wherein:
M1=(c θ1234+sθ14)ax+(sθ1234-cθ14)ay-sθ234az
M2=-c θ123ax-sθ123ay-cθ23az
Now, θ1、θ2、θ3、θ4、θ5It is known.In formula before substitution, it can obtain:
N=c θ12345+sθ145+cθ1235
O=s θ14-cθ1234
And then:
Wherein:
Step 5:There is robot degeneration in above-mentioned Analysis of Inverse Kinematics to the six axles six degree of freedom serial manipulator to ask Topic carries out related supplemental and explanation.
Because the kinematics analysis mode of the serial manipulator to the type is different from conventional Calculation Method, and in method In also related to the conversion between two kinds of different type serial manipulators, conventional robot kinematics' Singularity Analysis The problem of being not enough to tackle in the presence of above-mentioned analysis.Therefore, to using the analysis to degenerate case present in motion process To solve the multiresolution issue occurred in Robotic inverse kinematics analysis.
When robot loses 1 free degree, and therefore do not moved back by Ji Cheng robots during desired state motion Change.Robot can degenerate under two conditions:
(1) joint of robot reaches its physics limit and can not further moved;
(2) if the z-axis in two similar joints becomes conllinear, robot may be changed into degeneration shape in its working space State.
This means now whichever joint motions will all produce same motion, controller, which can not also determine bottom, is Which moving in joint.No matter any situation, the available number of degrees of freedom, of robot is less than 6, thus the equation of robot without Solution.
1. pair generally analysis of robot degenerate problem
As s θ5When=0, robot is in Singularity, meets degenerative conditions.Now, joint shaft 4 and 6 overlaps, and can only solve Go out θ6With θ4And with difference.In Singularity, θ can be arbitrarily chosen4Angle, calculating corresponding θ6Value.
But such case is existed only among the serial manipulator for meeting Pepeer criterions, in the six axles six degree of freedom It is that degenerate problem is not present in serial manipulator, as shown in Figure 2.Because under the type serial manipulator, it is connected with each other Conllinear problem will not occur for robot.
When Robotic inverse kinematics are analyzed, first the given data of Analysis of Inverse Kinematics is pre-processed.And And this processing mode be exactly be by it is original, be unsatisfactory for pieper criterion serial manipulators and be converted into meeting pieper criterions Serial manipulator.Therefore, during actually solving, or exist and work as s θ5=0 robot meets the situation of degenerative conditions. Below, this degenerate case is carried out to detailed analysis.
Based on axle 4 for coordinate, by s θ5=0, c θ5=± 1 tries to achieveTransition matrix:
Make (1,2), (3,2) are equal obtains:
θ4±θ6(the o of=A tan 2x1-oy1,-ox123-oy123+oz23)
By previously mentioned ox5、oy5、oz5Parameter, θ can be obtained6
θ6(the o of=A tan 2x, ox5)
It can to sum up obtain:
Wherein:The selection of the positive negativity of above formula is by c θ5Positive and negative judge.
2. the analysis that pair special circumstances are under degenerate state
Above the Singularity of the axle six degree of freedom serial manipulator of the type six is generally analyzed.It is real On border, when the robot motion is to a certain particular point, unusual situation also occurs.Due to such case only exist with The state of certain a kind of location point.Therefore it is referred to as Singularity in particular cases.
Situation 1:As shown in figure 4, work as d2=d3, and the axis of joint shaft 1 and 4 is (or when coincidence, under this state θ can not be solved1And θ4.Such case there is more than among the robot model after change, in original robot model Singularity be present.
Such case can occur in both cases:First, such location point is the transit point of movement locus;Second, such Location point is movement locus starting point.
If such location point is the transit point of movement locus, can solve this problem using " inheritance act ". That is, when there is such case, θ1Value acquiescence inherit θ under neighbouring position dotted state1Value, it is and then right θ4Solved.
If such location point is the starting point of movement locus, can solve this problem using " Reconstruction Method ". That is, when there is such case, θ1Value be defaulted as 0, so as to afterwards to θ4Solved.
Situation 2:As shown in figure 5, work as d2=d3, and when the axis coincidence of joint shaft 1 and 6, can not be solved under this state Go out θ1And θ6.Such case is existed only among the robot model after conversion, is to be not present very in original robot model Dystopy shape problem.But in above-mentioned Analysis of Inverse Kinematics, the state of this kind of Singularity still occurs.
In this case, can solve this problem using " sciagraphy ".There is such case that is, working as When, θ can be obtained by the projection that robot is formed in xoy planes1Numerical value.
It is worth it is to be noted that:It is this to ask method to draw two different results.

Claims (6)

1. the control method based on the six degree of freedom serial manipulator for being unsatisfactory for pipper criterions, it is characterised in that:Including following Step:
Step 1:The six axle six degree of freedom serial manipulators for being unsatisfactory for pipper criterions are changed on an equal basis, transform to it Meet on six axle six degree of freedom serial manipulator models of pipper criterions;Step 2:Connected to the six axles six degree of freedom When robot carries out Forward Kinematics Analysis, the direction vector parameter of biasing is reserved;Step 3:It is inclined according to what is tried to achieve before The direction vector parameter put, changes in coordinates is carried out to the six axles six degree of freedom serial manipulator, is transformed to satisfaction Six axle six degree of freedom serial manipulators of pipper criterions;Step 4:Analysis of Inverse Kinematics and meter are carried out to it based on analytic method Calculate;Step 5:The problem of progress during Analysis of Inverse Kinematics is carried out be present to the six axles six degree of freedom serial manipulator Related supplemental and explanation;Step 6:The kinematic parameter solved, motion control card is dealt into, to being unsatisfactory for pipper criterions six Axle six degree of freedom serial manipulator is controlled.
2. the control method as claimed in claim 1 based on the six degree of freedom serial manipulator for being unsatisfactory for pipper criterions, its It is characterised by:
Step 1:The six axles six degree of freedom serial manipulator is changed on an equal basis, it is transformed into and meets pipper criterions Six axle six degree of freedom serial manipulator models;
Although the type serial manipulator is unsatisfactory for pieper criterions, by after corresponding deformation, can still be joined by DH Exponential model calculates, and the connect end effector arm of mechanical robot of six axles has been carried out into conversion to a certain extent so that end is most The extended line of three cradle head axis intersects at a point afterwards.
3. the control method as claimed in claim 2 based on the six degree of freedom serial manipulator for being unsatisfactory for pipper criterions, its It is characterised by:Step 2:Robot Forward Kinematics Analysis is exactly known rod member geometric parameter and joint angle vector, asks robot last Hold position and posture of the actuator relative to fixed reference coordinate;
Before robot Forward Kinematics Analysis is carried out, conversion of the coordinate system { i } relative to coordinate system { i-1 } is established;It is general this Individual conversion is made up of four parameters;To any given robot, this conversion is the function of only one variable, in addition Three parameters are determined by mechanical system;Transformation parameter such as following table institute between the connecting rod of the six axles six degree of freedom serial manipulator Show;
Wherein:
aiFor along the X under the i-th coordinate systemiAxle, the Z under the i-th coordinate systemiDirection of principal axis is moved to Z under i+1 coordinate systemi+1Axle Distance;
αiFor the X under the i-th coordinate systemiAxle, the Z under the i-th coordinate systemiThe Z that axle is rotated under i+1 coordinate systemi+1The angle of axle Degree;
diFor along the Z under the i-th coordinate systemiAxle, the X under the i-th -1 coordinate systemi-1The X that axle is moved under the i-th coordinate systemiAxle away from From;
θiFor the Z under the i-th coordinate systemiAxle, the X under the i-th -1 coordinate systemi-1The X that axle is rotated under the i-th coordinate systemiThe angle of axle Degree;
The normal solution computing of robot is represented with homogeneous coordinate transformation matrix:(Represent it must is from i coordinates to j coordinate systems 4 × 4 transfer matrix;In formula:By s θiWhat is represented is sin θi;cθiThat represent is cos θ1.And i=1,2,3...6)
<mrow> <mmultiscripts> <mi>T</mi> <mn>6</mn> <mn>0</mn> </mmultiscripts> <mo>=</mo> <mmultiscripts> <mi>T</mi> <mn>3</mn> <mn>0</mn> </mmultiscripts> <mmultiscripts> <mi>T</mi> <mn>5</mn> <mn>3</mn> </mmultiscripts> <mmultiscripts> <mi>T</mi> <mrow> <mn>6</mn> <mo>*</mo> </mrow> <mn>5</mn> </mmultiscripts> <mmultiscripts> <mi>T</mi> <mn>6</mn> <mrow> <mn>6</mn> <mo>*</mo> </mrow> </mmultiscripts> <mo>=</mo> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <mrow> <msub> <mi>o</mi> <mi>x</mi> </msub> <msub> <mi>c&amp;theta;</mi> <mn>6</mn> </msub> <mo>-</mo> <msub> <mi>n</mi> <mi>x</mi> </msub> <msub> <mi>s&amp;theta;</mi> <mn>6</mn> </msub> </mrow> </mtd> <mtd> <mrow> <mo>-</mo> <mrow> <mo>(</mo> <mrow> <msub> <mi>n</mi> <mi>x</mi> </msub> <msub> <mi>c&amp;theta;</mi> <mn>6</mn> </msub> <mo>+</mo> <msub> <mi>o</mi> <mi>x</mi> </msub> <msub> <mi>s&amp;theta;</mi> <mn>6</mn> </msub> </mrow> <mo>)</mo> </mrow> </mrow> </mtd> <mtd> <msub> <mi>a</mi> <mi>x</mi> </msub> </mtd> <mtd> <msub> <mi>p</mi> <mi>x</mi> </msub> </mtd> </mtr> <mtr> <mtd> <mrow> <msub> <mi>o</mi> <mi>y</mi> </msub> <msub> <mi>c&amp;theta;</mi> <mn>6</mn> </msub> <mo>-</mo> <msub> <mi>n</mi> <mi>y</mi> </msub> <msub> <mi>s&amp;theta;</mi> <mn>6</mn> </msub> </mrow> </mtd> <mtd> <mrow> <mo>-</mo> <mrow> <mo>(</mo> <mrow> <msub> <mi>n</mi> <mi>y</mi> </msub> <msub> <mi>c&amp;theta;</mi> <mn>6</mn> </msub> <mo>+</mo> <msub> <mi>o</mi> <mi>y</mi> </msub> <msub> <mi>s&amp;theta;</mi> <mn>6</mn> </msub> </mrow> <mo>)</mo> </mrow> </mrow> </mtd> <mtd> <msub> <mi>a</mi> <mi>y</mi> </msub> </mtd> <mtd> <msub> <mi>p</mi> <mi>y</mi> </msub> </mtd> </mtr> <mtr> <mtd> <mrow> <msub> <mi>o</mi> <mi>z</mi> </msub> <msub> <mi>c&amp;theta;</mi> <mn>6</mn> </msub> <mo>-</mo> <msub> <mi>n</mi> <mi>z</mi> </msub> <msub> <mi>s&amp;theta;</mi> <mn>6</mn> </msub> </mrow> </mtd> <mtd> <mrow> <mo>-</mo> <mrow> <mo>(</mo> <mrow> <msub> <mi>n</mi> <mi>z</mi> </msub> <msub> <mi>c&amp;theta;</mi> <mn>6</mn> </msub> <mo>+</mo> <msub> <mi>o</mi> <mi>z</mi> </msub> <msub> <mi>s&amp;theta;</mi> <mn>6</mn> </msub> </mrow> <mo>)</mo> </mrow> </mrow> </mtd> <mtd> <msub> <mi>a</mi> <mi>z</mi> </msub> </mtd> <mtd> <msub> <mi>p</mi> <mi>z</mi> </msub> </mtd> </mtr> <mtr> <mtd> <mn>0</mn> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mn>1</mn> </mtd> </mtr> </mtable> </mfenced> </mrow>
Wherein:
θ2323
nx=c θ12345+sθ145+cθ1235
ny=s θ1235+sθ12345-cθ145
nz=c θ235-sθ2345
ox=s θ14-cθ1234
oy=-c θ14-sθ1234
oz=s θ234
ax=c θ12345-cθ1235+sθ145
ay=s θ12345-sθ1235-cθ145
az=-s θ2345-cθ235
px=-d4123+l212+d31-d21+d6ax+l6ox
py=-d4123+l212+d21-d31+d6ay+l6oy
pz=-d423-l22+d6az+l6oz
It is worth noting that:Here the o out calculatedx、oy、ozWith the o hereinafter quotedx5、oy5、oz5Parameter values are consistent 's;By θ6=0 brings into above-mentioned equation, it is seen that:The direction vector at the exactly bigoted place of these three parameters description;Its In, px’、py’、pz' represent be end effector coordinate position, nx、ny、nz、ox、oy、oz、ax、ay、azWhat is represented is end The coordinate posture of actuator.
4. the control method as claimed in claim 3 based on the six degree of freedom serial manipulator for being unsatisfactory for pipper criterions, its It is characterised by:
Step 3:According to the direction vector parameter for the biasing tried to achieve before, the six axles six degree of freedom serial manipulator is carried out Coordinate transform, transformed to and met on six axle six degree of freedom serial manipulators of pipper criterions;
Analysis of Inverse Kinematics is the geometric parameter of known machine people's rod member, gives end effector relative to fixed reference frame Desired locations and posture (wherein, px’、py’、pz' represent be end effector coordinate position, nx、ny、nz、ox、oy、oz、 ax、ay、azWhat is represented is the coordinate posture of end effector), can analysis robot make its end effector reach this expection Pose;
For Analysis of Inverse Kinematics, the coordinate of ending coordinates system posture is known, that is to say, thatKnown to parameters; But due to θ6Uncertainty, can not pass through, it is known that byTo drawTransformation matrix;Herein, " standardization " can be taken such to solve the problems, such as;
Above-mentioned Forward Kinematics Analysis is solved what is comeThe o drawnx5、oy5、oz5Whole Analysis of Inverse Kinematics is added to, And given data is pre-processed:
<mfenced open = "{" close = ""> <mtable> <mtr> <mtd> <mrow> <msub> <mi>p</mi> <mi>x</mi> </msub> <mo>=</mo> <msup> <msub> <mi>p</mi> <mi>x</mi> </msub> <mo>&amp;prime;</mo> </msup> <mo>-</mo> <msub> <mi>d</mi> <mn>6</mn> </msub> <msub> <mi>a</mi> <mi>x</mi> </msub> <mo>-</mo> <msub> <mi>l</mi> <mn>6</mn> </msub> <msub> <mi>o</mi> <mrow> <mi>x</mi> <mn>5</mn> </mrow> </msub> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msub> <mi>p</mi> <mi>y</mi> </msub> <mo>=</mo> <msup> <msub> <mi>p</mi> <mi>y</mi> </msub> <mo>&amp;prime;</mo> </msup> <mo>-</mo> <msub> <mi>d</mi> <mn>6</mn> </msub> <msub> <mi>a</mi> <mi>y</mi> </msub> <mo>-</mo> <msub> <mi>l</mi> <mn>6</mn> </msub> <msub> <mi>o</mi> <mrow> <mi>y</mi> <mn>5</mn> </mrow> </msub> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msub> <mi>p</mi> <mi>z</mi> </msub> <mo>=</mo> <msup> <msub> <mi>p</mi> <mi>z</mi> </msub> <mo>&amp;prime;</mo> </msup> <mo>-</mo> <msub> <mi>d</mi> <mn>6</mn> </msub> <msub> <mi>a</mi> <mi>z</mi> </msub> <mo>-</mo> <msub> <mi>l</mi> <mn>6</mn> </msub> <msub> <mi>o</mi> <mrow> <mi>z</mi> <mn>5</mn> </mrow> </msub> </mrow> </mtd> </mtr> </mtable> </mfenced>
Wherein:px′、py′、pz' it is known parameters, by the p after changex、py、pzBring into matrix;
Why given data is pre-processed, be because the serial manipulator is unsatisfactory for pieper criterions;With it is traditional Multiaxis serial manipulator is compared, the extended line of the joints axes in three joints in end of the type multiaxis serial manipulator Fail to intersect at a point;The relative transverse direction of one is had more so between them on its position of some on end effector Biasing;Also, the position and posture of this biasing can also change and change with the motion state in last three joints;Because this is inclined Putting the location status of length has certain uncertainty, and its position and posture to robot end also has significant shadow Ring;Only take to its ending coordinates parameter px′、py′、pz' pretreated method is carried out, it could eliminate by above-mentioned biasing not It is determined that to the influence caused by Solve problems, derived so as to carry out further Analysis of Inverse Kinematics and process;It is in fact, this Processing mode be by it is original, be unsatisfactory for the serial machine that pieper criterion serial manipulators are converted into meeting pieper criterions People.
5. the control method as claimed in claim 4 based on the six degree of freedom serial manipulator for being unsatisfactory for pipper criterions, its It is characterised by:Step 4:Analysis of Inverse Kinematics and calculating are carried out to it based on analytic method;
By formulaIt can obtain:
<mrow> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <mrow> <msub> <mi>c&amp;theta;</mi> <mn>1</mn> </msub> </mrow> </mtd> <mtd> <mrow> <msub> <mi>s&amp;theta;</mi> <mn>1</mn> </msub> </mrow> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mn>0</mn> </mtd> </mtr> <mtr> <mtd> <mrow> <mo>-</mo> <msub> <mi>s&amp;theta;</mi> <mn>1</mn> </msub> </mrow> </mtd> <mtd> <mrow> <msub> <mi>c&amp;theta;</mi> <mn>1</mn> </msub> </mrow> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mn>0</mn> </mtd> </mtr> <mtr> <mtd> <mn>0</mn> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mn>1</mn> </mtd> <mtd> <mn>0</mn> </mtd> </mtr> <mtr> <mtd> <mn>0</mn> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mn>1</mn> </mtd> </mtr> </mtable> </mfenced> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <msub> <mi>n</mi> <mi>x</mi> </msub> </mtd> <mtd> <msub> <mi>o</mi> <mi>x</mi> </msub> </mtd> <mtd> <msub> <mi>a</mi> <mi>x</mi> </msub> </mtd> <mtd> <msub> <mi>p</mi> <mi>x</mi> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>n</mi> <mi>y</mi> </msub> </mtd> <mtd> <msub> <mi>o</mi> <mi>y</mi> </msub> </mtd> <mtd> <msub> <mi>a</mi> <mi>y</mi> </msub> </mtd> <mtd> <msub> <mi>p</mi> <mi>y</mi> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>n</mi> <mi>z</mi> </msub> </mtd> <mtd> <msub> <mi>o</mi> <mi>z</mi> </msub> </mtd> <mtd> <msub> <mi>a</mi> <mi>z</mi> </msub> </mtd> <mtd> <msub> <mi>p</mi> <mi>z</mi> </msub> </mtd> </mtr> <mtr> <mtd> <mn>0</mn> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mn>1</mn> </mtd> </mtr> </mtable> </mfenced> <mo>=</mo> <mmultiscripts> <mi>T</mi> <mn>6</mn> <mn>1</mn> </mmultiscripts> </mrow>
It is equal by the element (2,4) on above-mentioned equation both sides for first axis joint angle, obtain:
<mrow> <msub> <mi>&amp;theta;</mi> <mn>1</mn> </msub> <mo>=</mo> <mi>A</mi> <mi> </mi> <mi>t</mi> <mi>a</mi> <mi>n</mi> <mn>2</mn> <mrow> <mo>(</mo> <msub> <mi>p</mi> <mi>y</mi> </msub> <mo>,</mo> <msub> <mi>p</mi> <mi>x</mi> </msub> <mo>)</mo> </mrow> <mo>+</mo> <mi>A</mi> <mi> </mi> <mi>tan</mi> <mn>2</mn> <mrow> <mo>(</mo> <mfrac> <mrow> <msub> <mi>d</mi> <mn>2</mn> </msub> <mo>-</mo> <msub> <mi>d</mi> <mn>3</mn> </msub> </mrow> <mi>&amp;rho;</mi> </mfrac> <mo>,</mo> <mo>&amp;PlusMinus;</mo> <mfrac> <msqrt> <mrow> <msup> <mi>&amp;rho;</mi> <mn>2</mn> </msup> <mo>-</mo> <msup> <mrow> <mo>(</mo> <msub> <mi>d</mi> <mn>2</mn> </msub> <mo>-</mo> <msub> <mi>d</mi> <mn>3</mn> </msub> <mo>)</mo> </mrow> <mn>2</mn> </msup> </mrow> </msqrt> <mi>&amp;rho;</mi> </mfrac> <mo>)</mo> </mrow> </mrow>
In formula:
<mrow> <mi>&amp;rho;</mi> <mo>=</mo> <msqrt> <mrow> <msup> <msub> <mi>p</mi> <mi>x</mi> </msub> <mn>2</mn> </msup> <mo>+</mo> <msup> <msub> <mi>p</mi> <mi>y</mi> </msub> <mn>2</mn> </msup> </mrow> </msqrt> </mrow>
Pay attention to:Sign, θ in above formula be present1There can be multiple solutions;Now, θ1It is known;
For the 3rd joint angles, by the element (Isosorbide-5-Nitrae) on above-mentioned equation both sides, (3,4), equal (wherein element (Isosorbide-5-Nitrae) refers to The row parameter of equation both sides the 1st row of matrix the 4th), obtain:
<mrow> <msub> <mi>&amp;theta;</mi> <mn>3</mn> </msub> <mo>=</mo> <mi>A</mi> <mi> </mi> <mi>tan</mi> <mn>2</mn> <mrow> <mo>(</mo> <mi>k</mi> <mo>,</mo> <mo>&amp;PlusMinus;</mo> <msqrt> <mrow> <mn>1</mn> <mo>-</mo> <msup> <mi>k</mi> <mn>2</mn> </msup> </mrow> </msqrt> <mo>)</mo> </mrow> </mrow>
Wherein:
<mrow> <mi>k</mi> <mo>=</mo> <mfrac> <mrow> <msup> <msub> <mi>d</mi> <mn>4</mn> </msub> <mn>2</mn> </msup> <mo>+</mo> <msup> <msub> <mi>l</mi> <mn>2</mn> </msub> <mn>2</mn> </msup> <mo>-</mo> <msup> <msub> <mi>p</mi> <mi>x</mi> </msub> <mn>2</mn> </msup> <mo>-</mo> <msup> <msub> <mi>p</mi> <mi>y</mi> </msub> <mn>2</mn> </msup> <mo>-</mo> <msup> <msub> <mi>p</mi> <mi>z</mi> </msub> <mn>2</mn> </msup> </mrow> <mrow> <mn>2</mn> <msub> <mi>l</mi> <mn>2</mn> </msub> <msub> <mi>d</mi> <mn>4</mn> </msub> </mrow> </mfrac> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <mrow> <msub> <mi>c&amp;theta;</mi> <mn>1</mn> </msub> <msub> <mi>c&amp;theta;</mi> <mn>23</mn> </msub> </mrow> </mtd> <mtd> <mrow> <msub> <mi>s&amp;theta;</mi> <mn>1</mn> </msub> <msub> <mi>c&amp;theta;</mi> <mn>23</mn> </msub> </mrow> </mtd> <mtd> <mrow> <mo>-</mo> <msub> <mi>s&amp;theta;</mi> <mn>23</mn> </msub> </mrow> </mtd> <mtd> <mrow> <mo>-</mo> <msub> <mi>l</mi> <mn>2</mn> </msub> <msub> <mi>c&amp;theta;</mi> <mn>3</mn> </msub> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mo>-</mo> <msub> <mi>c&amp;theta;</mi> <mn>1</mn> </msub> <msub> <mi>s&amp;theta;</mi> <mn>23</mn> </msub> </mrow> </mtd> <mtd> <mrow> <mo>-</mo> <msub> <mi>s&amp;theta;</mi> <mn>1</mn> </msub> <msub> <mi>s&amp;theta;</mi> <mn>23</mn> </msub> </mrow> </mtd> <mtd> <mrow> <mo>-</mo> <msub> <mi>c&amp;theta;</mi> <mn>23</mn> </msub> </mrow> </mtd> <mtd> <mrow> <msub> <mi>l</mi> <mn>2</mn> </msub> <msub> <mi>s&amp;theta;</mi> <mn>3</mn> </msub> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mo>-</mo> <msub> <mi>s&amp;theta;</mi> <mn>1</mn> </msub> </mrow> </mtd> <mtd> <mrow> <msub> <mi>c&amp;theta;</mi> <mn>1</mn> </msub> </mrow> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mn>0</mn> </mtd> </mtr> <mtr> <mtd> <mn>0</mn> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mn>1</mn> </mtd> </mtr> </mtable> </mfenced> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <msub> <mi>n</mi> <mi>x</mi> </msub> </mtd> <mtd> <msub> <mi>o</mi> <mi>x</mi> </msub> </mtd> <mtd> <msub> <mi>a</mi> <mi>x</mi> </msub> </mtd> <mtd> <msub> <mi>p</mi> <mi>x</mi> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>n</mi> <mi>y</mi> </msub> </mtd> <mtd> <msub> <mi>o</mi> <mi>y</mi> </msub> </mtd> <mtd> <msub> <mi>a</mi> <mi>y</mi> </msub> </mtd> <mtd> <msub> <mi>p</mi> <mi>y</mi> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>n</mi> <mi>z</mi> </msub> </mtd> <mtd> <msub> <mi>o</mi> <mi>z</mi> </msub> </mtd> <mtd> <msub> <mi>a</mi> <mi>z</mi> </msub> </mtd> <mtd> <msub> <mi>p</mi> <mi>z</mi> </msub> </mtd> </mtr> <mtr> <mtd> <mn>0</mn> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mn>1</mn> </mtd> </mtr> </mtable> </mfenced> <mo>=</mo> <mmultiscripts> <mi>T</mi> <mn>6</mn> <mn>3</mn> </mmultiscripts> </mrow>
It is equal by the element (Isosorbide-5-Nitrae) (2,4) on above-mentioned equation both sides for second joint angles, obtain:
<mrow> <msub> <mi>&amp;theta;</mi> <mn>23</mn> </msub> <mo>=</mo> <mi>A</mi> <mi> </mi> <mi>t</mi> <mi>a</mi> <mi>n</mi> <mn>2</mn> <mrow> <mo>(</mo> <mfrac> <mrow> <msub> <mi>l</mi> <mn>2</mn> </msub> <msub> <mi>s&amp;theta;</mi> <mn>3</mn> </msub> <mo>-</mo> <msub> <mi>d</mi> <mn>4</mn> </msub> </mrow> <msqrt> <mrow> <msup> <mrow> <mo>(</mo> <msub> <mi>p</mi> <mi>x</mi> </msub> <msub> <mi>c&amp;theta;</mi> <mn>1</mn> </msub> <mo>+</mo> <msub> <mi>p</mi> <mi>y</mi> </msub> <msub> <mi>s&amp;theta;</mi> <mn>1</mn> </msub> <mo>)</mo> </mrow> <mn>2</mn> </msup> <mo>+</mo> <msup> <msub> <mi>p</mi> <mi>z</mi> </msub> <mn>2</mn> </msup> </mrow> </msqrt> </mfrac> <mo>,</mo> <mfrac> <mrow> <msub> <mi>l</mi> <mn>2</mn> </msub> <msub> <mi>c&amp;theta;</mi> <mn>3</mn> </msub> </mrow> <msqrt> <mrow> <msup> <mrow> <mo>(</mo> <mi>p</mi> <mi>x</mi> <mi>c</mi> <mi>&amp;theta;</mi> <mn>1</mn> <mo>+</mo> <mi>p</mi> <mi>y</mi> <mi>s</mi> <mi>&amp;theta;</mi> <mn>1</mn> <mo>)</mo> </mrow> <mn>2</mn> </msup> <mo>+</mo> <msup> <mi>pz</mi> <mn>2</mn> </msup> </mrow> </msqrt> </mfrac> <mo>)</mo> </mrow> <mo>-</mo> <mi>A</mi> <mi> </mi> <mi>tan</mi> <mn>2</mn> <mrow> <mo>(</mo> <msub> <mi>p</mi> <mi>z</mi> </msub> <mo>,</mo> <msub> <mi>p</mi> <mi>x</mi> </msub> <msub> <mi>c&amp;theta;</mi> <mn>1</mn> </msub> <mo>+</mo> <msub> <mi>p</mi> <mi>y</mi> </msub> <msub> <mi>s&amp;theta;</mi> <mn>1</mn> </msub> <mo>)</mo> </mrow> </mrow>
So, by θ2323It can obtain:
<mrow> <msub> <mi>&amp;theta;</mi> <mn>2</mn> </msub> <mo>=</mo> <mi>A</mi> <mi> </mi> <mi>tan</mi> <mn>2</mn> <mrow> <mo>(</mo> <mrow> <mfrac> <mrow> <msub> <mi>l</mi> <mn>2</mn> </msub> <msub> <mi>s&amp;theta;</mi> <mn>3</mn> </msub> <mo>-</mo> <msub> <mi>d</mi> <mn>4</mn> </msub> </mrow> <msqrt> <mrow> <msup> <mrow> <mo>(</mo> <mrow> <msub> <mi>p</mi> <mi>x</mi> </msub> <msub> <mi>c&amp;theta;</mi> <mn>1</mn> </msub> <mo>+</mo> <msub> <mi>p</mi> <mi>y</mi> </msub> <msub> <mi>s&amp;theta;</mi> <mn>1</mn> </msub> </mrow> <mo>)</mo> </mrow> <mn>2</mn> </msup> <mo>+</mo> <msup> <msub> <mi>p</mi> <mi>z</mi> </msub> <mn>2</mn> </msup> </mrow> </msqrt> </mfrac> <mo>,</mo> <mfrac> <mrow> <msub> <mi>l</mi> <mn>2</mn> </msub> <msub> <mi>c&amp;theta;</mi> <mn>3</mn> </msub> </mrow> <msqrt> <mrow> <msup> <mrow> <mo>(</mo> <mrow> <msub> <mi>p</mi> <mi>x</mi> </msub> <msub> <mi>c&amp;theta;</mi> <mn>1</mn> </msub> <mo>+</mo> <msub> <mi>p</mi> <mi>y</mi> </msub> <msub> <mi>s&amp;theta;</mi> <mn>1</mn> </msub> </mrow> <mo>)</mo> </mrow> <mn>2</mn> </msup> <mo>+</mo> <msup> <msub> <mi>p</mi> <mi>z</mi> </msub> <mn>2</mn> </msup> </mrow> </msqrt> </mfrac> </mrow> <mo>)</mo> </mrow> <mo>-</mo> <mi>A</mi> <mi> </mi> <mi>tan</mi> <mn>2</mn> <mrow> <mo>(</mo> <mrow> <msub> <mi>p</mi> <mi>z</mi> </msub> <mo>,</mo> <msub> <mi>p</mi> <mi>x</mi> </msub> <msub> <mi>c&amp;theta;</mi> <mn>1</mn> </msub> <mo>+</mo> <msub> <mi>p</mi> <mi>y</mi> </msub> <msub> <mi>s&amp;theta;</mi> <mn>1</mn> </msub> </mrow> <mo>)</mo> </mrow> <mo>-</mo> <mi>A</mi> <mi> </mi> <mi>tan</mi> <mn>2</mn> <mrow> <mo>(</mo> <mrow> <mi>k</mi> <mo>,</mo> <mo>&amp;PlusMinus;</mo> <msqrt> <mrow> <mn>1</mn> <mo>-</mo> <msup> <mi>k</mi> <mn>2</mn> </msup> </mrow> </msqrt> </mrow> <mo>)</mo> </mrow> </mrow>
According to θ1And θ3Four kinds of solution may combine, and can obtain corresponding four kinds of probable values by above formula, then can obtain θ2Four Kind may solve;
For the 4th joint angles, by the element (1,3) on above-mentioned equation both sides, (3,3) are equal obtains:
θ4=Atan2 (s θ1ax-cθ1ay, c θ123ax+sθ123ay-sθ23az)
As s θ5When=0, robot is in Singularity, meets degenerative conditions;Now, joint shaft 4 and 6 overlaps, and can only solve θ6 With θ4And with difference;In Singularity, θ can be arbitrarily chosen4Angle, calculating corresponding θ6Value;But such case is only It is present among the serial manipulator for meeting Pepeer criterions;This problem will be discussed in detail in next operation:
<mrow> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <mrow> <msub> <mi>c&amp;theta;</mi> <mn>1</mn> </msub> <msub> <mi>c&amp;theta;</mi> <mn>23</mn> </msub> <msub> <mi>c&amp;theta;</mi> <mn>4</mn> </msub> <mo>+</mo> <msub> <mi>s&amp;theta;</mi> <mn>1</mn> </msub> <msub> <mi>s&amp;theta;</mi> <mn>4</mn> </msub> </mrow> </mtd> <mtd> <mrow> <msub> <mi>s&amp;theta;</mi> <mn>1</mn> </msub> <msub> <mi>c&amp;theta;</mi> <mn>23</mn> </msub> <msub> <mi>c&amp;theta;</mi> <mn>4</mn> </msub> <mo>-</mo> <msub> <mi>c&amp;theta;</mi> <mn>1</mn> </msub> <msub> <mi>s&amp;theta;</mi> <mn>4</mn> </msub> </mrow> </mtd> <mtd> <mrow> <mo>-</mo> <msub> <mi>s&amp;theta;</mi> <mn>23</mn> </msub> <msub> <mi>c&amp;theta;</mi> <mn>4</mn> </msub> </mrow> </mtd> <mtd> <mrow> <mo>-</mo> <msub> <mi>d</mi> <mn>4</mn> </msub> <msub> <mi>c&amp;theta;</mi> <mn>1</mn> </msub> <msub> <mi>s&amp;theta;</mi> <mn>23</mn> </msub> <mo>+</mo> <msub> <mi>l</mi> <mn>2</mn> </msub> <msub> <mi>c&amp;theta;</mi> <mn>1</mn> </msub> <msub> <mi>c&amp;theta;</mi> <mn>2</mn> </msub> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mo>-</mo> <msub> <mi>c&amp;theta;</mi> <mn>1</mn> </msub> <msub> <mi>c&amp;theta;</mi> <mn>23</mn> </msub> <msub> <mi>s&amp;theta;</mi> <mn>4</mn> </msub> <mo>+</mo> <msub> <mi>s&amp;theta;</mi> <mn>1</mn> </msub> <msub> <mi>c&amp;theta;</mi> <mn>4</mn> </msub> </mrow> </mtd> <mtd> <mrow> <mo>-</mo> <msub> <mi>s&amp;theta;</mi> <mn>1</mn> </msub> <msub> <mi>c&amp;theta;</mi> <mn>23</mn> </msub> <msub> <mi>s&amp;theta;</mi> <mn>4</mn> </msub> <mo>-</mo> <msub> <mi>c&amp;theta;</mi> <mn>1</mn> </msub> <msub> <mi>c&amp;theta;</mi> <mn>4</mn> </msub> </mrow> </mtd> <mtd> <mrow> <msub> <mi>s&amp;theta;</mi> <mn>23</mn> </msub> <msub> <mi>s&amp;theta;</mi> <mn>4</mn> </msub> </mrow> </mtd> <mtd> <mrow> <mo>-</mo> <msub> <mi>d</mi> <mn>4</mn> </msub> <msub> <mi>s&amp;theta;</mi> <mn>1</mn> </msub> <msub> <mi>s&amp;theta;</mi> <mn>23</mn> </msub> <mo>+</mo> <msub> <mi>l</mi> <mn>2</mn> </msub> <msub> <mi>s&amp;theta;</mi> <mn>1</mn> </msub> <msub> <mi>c&amp;theta;</mi> <mn>2</mn> </msub> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mo>-</mo> <msub> <mi>c&amp;theta;</mi> <mn>1</mn> </msub> <msub> <mi>s&amp;theta;</mi> <mn>23</mn> </msub> </mrow> </mtd> <mtd> <mrow> <mo>-</mo> <msub> <mi>s&amp;theta;</mi> <mn>1</mn> </msub> <msub> <mi>s&amp;theta;</mi> <mn>23</mn> </msub> </mrow> </mtd> <mtd> <mrow> <mo>-</mo> <msub> <mi>c&amp;theta;</mi> <mn>23</mn> </msub> </mrow> </mtd> <mtd> <mrow> <mo>-</mo> <msub> <mi>d</mi> <mn>4</mn> </msub> <msub> <mi>c&amp;theta;</mi> <mn>23</mn> </msub> <mo>-</mo> <msub> <mi>l</mi> <mn>2</mn> </msub> <msub> <mi>s&amp;theta;</mi> <mn>2</mn> </msub> </mrow> </mtd> </mtr> <mtr> <mtd> <mn>0</mn> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mn>1</mn> </mtd> </mtr> </mtable> </mfenced> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <msub> <mi>n</mi> <mi>x</mi> </msub> </mtd> <mtd> <msub> <mi>o</mi> <mi>x</mi> </msub> </mtd> <mtd> <msub> <mi>a</mi> <mi>x</mi> </msub> </mtd> <mtd> <msub> <mi>p</mi> <mi>x</mi> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>n</mi> <mi>y</mi> </msub> </mtd> <mtd> <msub> <mi>o</mi> <mi>y</mi> </msub> </mtd> <mtd> <msub> <mi>a</mi> <mi>y</mi> </msub> </mtd> <mtd> <msub> <mi>p</mi> <mi>y</mi> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>n</mi> <mi>z</mi> </msub> </mtd> <mtd> <msub> <mi>o</mi> <mi>z</mi> </msub> </mtd> <mtd> <msub> <mi>a</mi> <mi>z</mi> </msub> </mtd> <mtd> <msub> <mi>p</mi> <mi>z</mi> </msub> </mtd> </mtr> <mtr> <mtd> <mn>0</mn> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mn>1</mn> </mtd> </mtr> </mtable> </mfenced> <mo>=</mo> <mmultiscripts> <mi>T</mi> <mn>5</mn> <mn>4</mn> </mmultiscripts> </mrow>
For the 5th axis joint angle, by the element (1,3) at above-mentioned equation both ends, (3,3) are equal obtains:
θ5=Atan2 (± M1, M2)
Wherein:
M1=(c θ1234+sθ14)ax+(sθ1234-cθ14)ay-sθ234az
M2=-c θ123ax-sθ123ay-cθ23az
For the 6th joint angles, it can obtain:
<mrow> <msub> <mi>&amp;theta;</mi> <mn>6</mn> </msub> <mo>=</mo> <mi>A</mi> <mi> </mi> <mi>tan</mi> <mn>2</mn> <mrow> <mo>(</mo> <mfrac> <mrow> <mo>-</mo> <msub> <mi>o</mi> <mi>x</mi> </msub> </mrow> <msqrt> <mrow> <msup> <mi>n</mi> <mn>2</mn> </msup> <mo>+</mo> <msup> <mi>o</mi> <mn>2</mn> </msup> </mrow> </msqrt> </mfrac> <mo>,</mo> <mfrac> <msub> <mi>n</mi> <mi>x</mi> </msub> <msqrt> <mrow> <msup> <mi>n</mi> <mn>2</mn> </msup> <mo>+</mo> <msup> <mi>o</mi> <mn>2</mn> </msup> </mrow> </msqrt> </mfrac> <mo>)</mo> </mrow> <mo>-</mo> <mi>A</mi> <mi> </mi> <mi>t</mi> <mi>a</mi> <mi>n</mi> <mn>2</mn> <mrow> <mo>(</mo> <mfrac> <mi>n</mi> <msqrt> <mrow> <msup> <mi>n</mi> <mn>2</mn> </msup> <mo>+</mo> <msup> <mi>o</mi> <mn>2</mn> </msup> </mrow> </msqrt> </mfrac> <mo>,</mo> <mfrac> <mi>o</mi> <msqrt> <mrow> <msup> <mi>n</mi> <mn>2</mn> </msup> <mo>+</mo> <msup> <mi>o</mi> <mn>2</mn> </msup> </mrow> </msqrt> </mfrac> <mo>)</mo> </mrow> </mrow>
Wherein:
N=c θ12345+sθ145+cθ1235
O=s θ14-cθ1234
6. the control method as claimed in claim 5 based on the six degree of freedom serial manipulator for being unsatisfactory for pipper criterions, its It is characterised by:
Step 5:1) to the analysis of generally robot degenerate problem
As s θ5When=0, robot is in Singularity, meets degenerative conditions;Now, joint shaft 4 and 6 overlaps, and can only solve θ6 With θ4And with difference;In Singularity, θ can be arbitrarily chosen4Angle, calculating corresponding θ6Value;
But such case is existed only among the serial manipulator for meeting Pepeer criterions, connected in the six axles six degree of freedom It is that degenerate problem is not present in robot;Because under the type serial manipulator, the robot being connected with each other is to send out Raw conllinear problem;
When Robotic inverse kinematics are analyzed, first the given data of Analysis of Inverse Kinematics is pre-processed;And this Kind of processing mode be exactly be by it is original, be unsatisfactory for the series connection that pieper criterion serial manipulators are converted into meeting pieper criterions Robot;Therefore, during actually solving, or exist and work as s θ5=0 robot meets the situation of degenerative conditions;Below, This degenerate case is carried out to detailed analysis;
<mrow> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <mrow> <msub> <mi>c&amp;theta;</mi> <mn>1</mn> </msub> <msub> <mi>c&amp;theta;</mi> <mn>23</mn> </msub> </mrow> </mtd> <mtd> <mrow> <msub> <mi>s&amp;theta;</mi> <mn>1</mn> </msub> <msub> <mi>c&amp;theta;</mi> <mn>23</mn> </msub> </mrow> </mtd> <mtd> <mrow> <mo>-</mo> <msub> <mi>s&amp;theta;</mi> <mn>23</mn> </msub> </mrow> </mtd> <mtd> <mrow> <mo>-</mo> <msub> <mi>l</mi> <mn>2</mn> </msub> <msub> <mi>c&amp;theta;</mi> <mn>3</mn> </msub> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mo>-</mo> <msub> <mi>c&amp;theta;</mi> <mn>1</mn> </msub> <msub> <mi>s&amp;theta;</mi> <mn>23</mn> </msub> </mrow> </mtd> <mtd> <mrow> <mo>-</mo> <msub> <mi>s&amp;theta;</mi> <mn>1</mn> </msub> <msub> <mi>s&amp;theta;</mi> <mn>23</mn> </msub> </mrow> </mtd> <mtd> <mrow> <mo>-</mo> <msub> <mi>c&amp;theta;</mi> <mn>23</mn> </msub> </mrow> </mtd> <mtd> <mrow> <msub> <mi>l</mi> <mn>2</mn> </msub> <msub> <mi>s&amp;theta;</mi> <mn>3</mn> </msub> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mo>-</mo> <msub> <mi>s&amp;theta;</mi> <mn>1</mn> </msub> </mrow> </mtd> <mtd> <mrow> <msub> <mi>c&amp;theta;</mi> <mn>1</mn> </msub> </mrow> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mn>0</mn> </mtd> </mtr> <mtr> <mtd> <mn>0</mn> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mn>1</mn> </mtd> </mtr> </mtable> </mfenced> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <msub> <mi>n</mi> <mi>x</mi> </msub> </mtd> <mtd> <msub> <mi>o</mi> <mi>x</mi> </msub> </mtd> <mtd> <msub> <mi>a</mi> <mi>x</mi> </msub> </mtd> <mtd> <msub> <mi>p</mi> <mi>x</mi> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>n</mi> <mi>y</mi> </msub> </mtd> <mtd> <msub> <mi>o</mi> <mi>y</mi> </msub> </mtd> <mtd> <msub> <mi>a</mi> <mi>y</mi> </msub> </mtd> <mtd> <msub> <mi>p</mi> <mi>y</mi> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>n</mi> <mi>z</mi> </msub> </mtd> <mtd> <msub> <mi>o</mi> <mi>z</mi> </msub> </mtd> <mtd> <msub> <mi>a</mi> <mi>z</mi> </msub> </mtd> <mtd> <msub> <mi>p</mi> <mi>z</mi> </msub> </mtd> </mtr> <mtr> <mtd> <mn>0</mn> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mn>1</mn> </mtd> </mtr> </mtable> </mfenced> <mo>=</mo> <mmultiscripts> <mi>T</mi> <mn>6</mn> <mn>3</mn> </mmultiscripts> </mrow>
Based on joint shaft 4 for coordinate, by s θs=O, c θ5=± 1 tries to achieveTransition matrix:
Make (1,2), (3,2) are equal obtains:
-ox1+oy1=c θ46±sθ46
θ4±θ6=Atan2 (ox1-oy1,-ox123-oy1cθ23+oz23)
By previously mentioned ox5、oy5、oz5Parameter, θ can be obtained6
θ6=Atan2 (ox, ox5)
It can to sum up obtain:
Wherein:The selection of the positive negativity of above formula is by c θ5Positive and negative judge;
2) analysis that special circumstances are under degenerate state
Above the Singularity of the six axles six degree of freedom serial manipulator is generally analyzed;In fact, When the robot motion is to a certain particular point, unusual situation also occurs;Due to such case only exist with it is a certain The state of class location point;Therefore it is referred to as Singularity in particular cases;
Situation 1:Work as d2=d3, and the axis of joint shaft 1 and 4 (or can not solve θ when coincidence, under this state1And θ4; Such case there is more than among the robot model after change, in original robot model and Singularity be present Problem;
Such case can occur in both cases:First, such location point is the transit point of movement locus;Second, such position Point is movement locus starting point;
If such location point is the transit point of movement locus, can solve this problem using " inheritance act ";Also It is to say, when there is such case, θ1Value acquiescence inherit θ under neighbouring position dotted state1Value, and then to θ4Enter Row solves;
If such location point is the starting point of movement locus, can solve this problem using " Reconstruction Method ";Also It is to say, when there is such case, θ1Value be defaulted as 0, so as to afterwards to θ4Solved;
Situation 2:Work as d2=d3, and when the axis coincidence of joint shaft 1 and 6, θ can not be solved under this state1And θ6;This feelings Condition is existed only among the robot model after conversion, is that Singularity problem is not present in original robot model; But in above-mentioned Analysis of Inverse Kinematics, the state of this kind of Singularity still occurs;
In this case, can solve this problem using " sciagraphy ";That is, when there is such case Wait, θ can be obtained by the projection that robot is formed in xoy planes1Numerical value.
CN201710902038.1A 2017-09-28 2017-09-28 Control method of six-degree-of-freedom series robot based on criterion of not meeting Pieper Active CN107791248B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710902038.1A CN107791248B (en) 2017-09-28 2017-09-28 Control method of six-degree-of-freedom series robot based on criterion of not meeting Pieper

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710902038.1A CN107791248B (en) 2017-09-28 2017-09-28 Control method of six-degree-of-freedom series robot based on criterion of not meeting Pieper

Publications (2)

Publication Number Publication Date
CN107791248A true CN107791248A (en) 2018-03-13
CN107791248B CN107791248B (en) 2021-04-30

Family

ID=61532943

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710902038.1A Active CN107791248B (en) 2017-09-28 2017-09-28 Control method of six-degree-of-freedom series robot based on criterion of not meeting Pieper

Country Status (1)

Country Link
CN (1) CN107791248B (en)

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108372499A (en) * 2018-05-15 2018-08-07 哈尔滨理工大学 Spatial cell robot via Self-reconfiguration cradle head module
CN108381540A (en) * 2018-05-15 2018-08-10 哈尔滨理工大学 A kind of space self-reorganization robot elbow turns cell
CN108942943A (en) * 2018-08-16 2018-12-07 居鹤华 Multi-axis robot positive kinematics calculation method based on axis invariant
CN112157653A (en) * 2020-09-11 2021-01-01 北京如影智能科技有限公司 Shielding detection method and device
CN112776341A (en) * 2020-12-15 2021-05-11 同济大学 Shaft action optimization method of six-shaft robot in three-dimensional printing equipment
CN113580135A (en) * 2021-08-09 2021-11-02 华中科技大学 Real-time inverse solution method for seven-axis robot with offset

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102243620A (en) * 2011-06-02 2011-11-16 安凯 Rapid solving method for inverse kinematics problem of six-joint mechanical arm
CN102637158A (en) * 2012-04-28 2012-08-15 谷菲 Inverse kinematics solution method for six-degree-of-freedom serial robot
CN103942427A (en) * 2014-04-11 2014-07-23 哈尔滨工程大学 Quick and simple method for solving inverse kinematics of six-degree-of-freedom mechanical arm
US20150283704A1 (en) * 2014-04-07 2015-10-08 Canon Kabushiki Kaisha Information processing apparatus and information processing method
CN105573143A (en) * 2015-11-30 2016-05-11 珞石(北京)科技有限公司 Inverse kinematics solving method for 6-DOF (degree of freedom) industrial robot
CN105975795A (en) * 2016-05-23 2016-09-28 湖北工业大学 High-precision multi-joint serial connection mechanical arm anti-kinematics solution
CN106363607A (en) * 2016-10-31 2017-02-01 北京控制工程研究所 Space manipulator system for capturing with strong bearing ability and large motion range

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102243620A (en) * 2011-06-02 2011-11-16 安凯 Rapid solving method for inverse kinematics problem of six-joint mechanical arm
CN102637158A (en) * 2012-04-28 2012-08-15 谷菲 Inverse kinematics solution method for six-degree-of-freedom serial robot
US20150283704A1 (en) * 2014-04-07 2015-10-08 Canon Kabushiki Kaisha Information processing apparatus and information processing method
CN103942427A (en) * 2014-04-11 2014-07-23 哈尔滨工程大学 Quick and simple method for solving inverse kinematics of six-degree-of-freedom mechanical arm
CN105573143A (en) * 2015-11-30 2016-05-11 珞石(北京)科技有限公司 Inverse kinematics solving method for 6-DOF (degree of freedom) industrial robot
CN105975795A (en) * 2016-05-23 2016-09-28 湖北工业大学 High-precision multi-joint serial connection mechanical arm anti-kinematics solution
CN106363607A (en) * 2016-10-31 2017-02-01 北京控制工程研究所 Space manipulator system for capturing with strong bearing ability and large motion range

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
于凌涛等: "一类不满足Pieper准则的机器人逆运动学解析解获取方法", 《机器人 ROBOT》 *
刘松国等: "6R机器人实时逆运动学算法研究", 《控制理论与应用》 *
唐奥林: "面向主从式微创外科手术机器人的遥操作运动控制策略研究", 《中国博士学位论文全文数据库 信息科技辑》 *

Cited By (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108372499A (en) * 2018-05-15 2018-08-07 哈尔滨理工大学 Spatial cell robot via Self-reconfiguration cradle head module
CN108381540A (en) * 2018-05-15 2018-08-10 哈尔滨理工大学 A kind of space self-reorganization robot elbow turns cell
CN108372499B (en) * 2018-05-15 2024-02-27 哈尔滨理工大学 Space cell robot self-reconstruction rotary joint module
CN108942943A (en) * 2018-08-16 2018-12-07 居鹤华 Multi-axis robot positive kinematics calculation method based on axis invariant
CN112157653A (en) * 2020-09-11 2021-01-01 北京如影智能科技有限公司 Shielding detection method and device
CN112157653B (en) * 2020-09-11 2022-02-01 北京如影智能科技有限公司 Shielding detection method and device
CN112776341A (en) * 2020-12-15 2021-05-11 同济大学 Shaft action optimization method of six-shaft robot in three-dimensional printing equipment
CN113580135A (en) * 2021-08-09 2021-11-02 华中科技大学 Real-time inverse solution method for seven-axis robot with offset

Also Published As

Publication number Publication date
CN107791248B (en) 2021-04-30

Similar Documents

Publication Publication Date Title
CN107791248A (en) Control method based on the six degree of freedom serial manipulator for being unsatisfactory for pipper criterions
CN105014677B (en) Vision Mechanical arm control method based on Camshift visual tracking and D-H modeling algorithm
CN109859275A (en) A kind of monocular vision hand and eye calibrating method of the rehabilitation mechanical arm based on S-R-S structure
CN107363813A (en) A kind of desktop industrial robot teaching system and method based on wearable device
CN107589934A (en) A kind of acquiring method of articulated manipulator inverse kinematics parsing solution
CN108241339A (en) The movement solution of apery mechanical arm and configuration control method
CN103085069B (en) Novel robot kinematics modeling method
CN106671079A (en) Motion control method for welding robot in coordination with positioner
CN105643619B (en) A kind of industrial robot instrument posture control method of use framework description
CN107756400A (en) A kind of 6R Robotic inverse kinematics geometry solving methods based on spinor theory
CN107584474A (en) A kind of spherical joint Dual-Arm Robot Coordinate Motion method based on geometric projection
CN106003034A (en) Master-slave robot control system and control method
CN108406769A (en) The joint rotation angle of serial manipulator determines method and device
CN107169196A (en) Dynamic modeling method of the robot for space from end effector to pedestal
CN101244561A (en) Kinematic inverse method for installation process of steam generator examining and repairing mechanical arm
CN109834706A (en) The method and device of kinematicsingularities is avoided in robot motion planning
CN111791234A (en) Anti-collision control algorithm for working positions of multiple robots in narrow space
CN111590567A (en) Space manipulator teleoperation planning method based on Omega handle
CN111496783A (en) Inverse kinematics solving method for 6R industrial robot
Bratchikov et al. Development of digital twin for robotic arm
CN116330267A (en) Control method based on industrial robot wrist singular point calculation
CN111515954A (en) Method for generating high-quality motion path of mechanical arm
Zhang et al. Kuka youBot arm shortest path planning based on geodesics
CN109866224A (en) A kind of robot Jacobian matrix calculation method, device and storage medium
Xuan et al. Reverse-driving trajectory planning and simulation of joint robot

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