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 PDFInfo
- 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
Links
Classifications
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B25—HAND TOOLS; PORTABLE POWER-DRIVEN TOOLS; MANIPULATORS
- B25J—MANIPULATORS; CHAMBERS PROVIDED WITH MANIPULATION DEVICES
- B25J9/00—Programme-controlled manipulators
- B25J9/16—Programme controls
- B25J9/1602—Programme controls characterised by the control system, structure, architecture
- B25J9/1607—Calculation of inertia, jacobian matrixes and inverses
-
- G—PHYSICS
- G05—CONTROLLING; REGULATING
- G05B—CONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
- G05B2219/00—Program-control systems
- G05B2219/30—Nc systems
- G05B2219/39—Robotics, robotics to robotics hand
- G05B2219/39073—Solve 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
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:
θ23=θ2+θ3
nx=c θ1cθ23cθ4cθ5+sθ1sθ4cθ5+cθ1sθ23sθ5
ny=s θ1sθ23sθ5+sθ1cθ23cθ4cθ5-cθ1sθ4cθ5
nz=c θ23Sθ5-sθ23cθ4cθ5
ox=s θ1cθ4-cθ1cθ23sθ4
oy=-c θ1cθ4-sθ1cθ23sθ4
oz=s θ23sθ4
ax=c θ1cθ23cθ4sθ5-cθ1sθ23cθ5+sθ1sθ4Sθ5
ay=s θ1cθ23cθ4sθ5-sθ1sθ23cθ5-cθ1sθ4sθ5
az=-s θ23cθ4sθ5-cθ23cθ5
px=-d4cθ1sθ23+l2cθ1cθ2+d3sθ1-d2sθ1+d6ax+l6ox
py=-d4sθ1sθ23+l2sθ1cθ2+d2cθ1-d3cθ1+d6ay+l6oy
pz=-d4cθ23-l2sθ2+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 θ23=θ2+θ3It 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 θ1cθ23ax+sθ1cθ23ay-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 θ1cθ23cθ4+sθ1sθ4)ax+(sθ1cθ23cθ4-cθ1sθ4)ay-sθ23cθ4az
M2=-c θ1sθ23ax-sθ1sθ23ay-cθ23az
For the 6th joint angles, it can obtain:
Wherein:
N=c θ1cθ23cθ4cθ5+sθ1sθ4cθ5+cθ1sθ23sθ5
O=s θ1cθ4-cθ1cθ23sθ4
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:
-oxsθ1+oycθ1=c θ4sθ6±sθ4cθ6
θ4±θ6(the o of=A tan 2xsθ1-oycθ1,-oxcθ1cθ23-oysθ1cθ23+ozsθ23)
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:
θ23=θ2+θ3
nx=c θ1cθ23cθ4cθ5+sθ1sθ4cθ5+cθ1Sθ23Sθ5
ny=s θ1sθ23Sθ5+sθ1cθ23cθ4cθ5-cθ1sθ4cθ5
nz=c θ23Sθ5-Sθ23cθ4cθ5
ox=S θ1cθ4-cθ1cθ23sθ4
oy=-c θ1cθ4-sθ1cθ23sθ4
oz=S θ23Sθ4
ax=c θ1cθ23cθ4sθ5-cθ1sθ23cθ5+sθ1sθ4Sθ5
ay=s θ1cθ23cθ4sθ5-sθ1sθ23cθ5-cθ1sθ4sθ5
az=-s θ23cθ4Sθ5-cθ23cθ5
px=-d4cθ1sθ23+l2cθ1cθ2+d3sθ1-d2sθ1+d6ax+l6ox
py=-d4Sθ1sθ23+l2sθ1cθ2+d2cθ1-d3cθ1+d6ay+l6oy
pz=-d4cθ23-l2sθ2+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:
sθ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:
pxcθ1+pysθ1=-d4sθ23+l2cθ2
pz=-d4cθ23-l2sθ2
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:
cθ1cθ23px+sθ1cθ23py-sθ23pz-l2cθ3=0
-cθ1sθ23px-sθ1sθ23py-cθ23pz+l2sθ3=d4
To sum up, θ2Solution can be write as:
So, by θ23=θ2+θ3It 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:
cθ1cθ23ax+sθ1cθ23ay-sθ23az=c θ4sθ5
-sθ1ax+cθ1ay=-s θ4sθ5
:θ4(the s θ of=A tan 21ax-cθ1ay, c θ1cθ23ax+sθ1cθ23ay-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θ1cθ23cθ4+sθ1sθ4)ax+(sθ1cθ23cθ4-cθ1sθ4)ay-sθ23cθ4az=s θ5
-cθ1sθ23ax-sθ1Sθ23ay-cθ23az=c θ5
:
θ5(± the M of=A tan 21, M2)
Wherein:
M1=(c θ1cθ23cθ4+sθ1sθ4)ax+(sθ1cθ23cθ4-cθ1sθ4)ay-sθ23cθ4az
M2=-c θ1sθ23ax-sθ1sθ23ay-cθ23az
Now, θ1、θ2、θ3、θ4、θ5It is known.In formula before substitution, it can obtain:
N=c θ1cθ23cθ4cθ5+sθ1sθ4cθ5+cθ1sθ23sθ5
O=s θ1cθ4-cθ1cθ23sθ4
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 2xsθ1-oycθ1,-oxcθ1cθ23-oysθ1cθ23+ozsθ23)
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&theta;</mi>
<mn>6</mn>
</msub>
<mo>-</mo>
<msub>
<mi>n</mi>
<mi>x</mi>
</msub>
<msub>
<mi>s&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&theta;</mi>
<mn>6</mn>
</msub>
<mo>+</mo>
<msub>
<mi>o</mi>
<mi>x</mi>
</msub>
<msub>
<mi>s&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&theta;</mi>
<mn>6</mn>
</msub>
<mo>-</mo>
<msub>
<mi>n</mi>
<mi>y</mi>
</msub>
<msub>
<mi>s&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&theta;</mi>
<mn>6</mn>
</msub>
<mo>+</mo>
<msub>
<mi>o</mi>
<mi>y</mi>
</msub>
<msub>
<mi>s&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&theta;</mi>
<mn>6</mn>
</msub>
<mo>-</mo>
<msub>
<mi>n</mi>
<mi>z</mi>
</msub>
<msub>
<mi>s&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&theta;</mi>
<mn>6</mn>
</msub>
<mo>+</mo>
<msub>
<mi>o</mi>
<mi>z</mi>
</msub>
<msub>
<mi>s&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:
θ23=θ2+θ3
nx=c θ1cθ23cθ4cθ5+sθ1sθ4cθ5+cθ1sθ23sθ5
ny=s θ1sθ23sθ5+sθ1cθ23cθ4cθ5-cθ1sθ4cθ5
nz=c θ23sθ5-sθ23cθ4cθ5
ox=s θ1cθ4-cθ1cθ23sθ4
oy=-c θ1cθ4-sθ1cθ23sθ4
oz=s θ23sθ4
ax=c θ1cθ23cθ4sθ5-cθ1sθ23cθ5+sθ1sθ4sθ5
ay=s θ1cθ23cθ4sθ5-sθ1sθ23cθ5-cθ1sθ4sθ5
az=-s θ23cθ4sθ5-cθ23cθ5
px=-d4cθ1sθ23+l2cθ1cθ2+d3sθ1-d2sθ1+d6ax+l6ox
py=-d4sθ1sθ23+l2sθ1cθ2+d2cθ1-d3cθ1+d6ay+l6oy
pz=-d4cθ23-l2sθ2+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>&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>&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>&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&theta;</mi>
<mn>1</mn>
</msub>
</mrow>
</mtd>
<mtd>
<mrow>
<msub>
<mi>s&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&theta;</mi>
<mn>1</mn>
</msub>
</mrow>
</mtd>
<mtd>
<mrow>
<msub>
<mi>c&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>&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>&rho;</mi>
</mfrac>
<mo>,</mo>
<mo>&PlusMinus;</mo>
<mfrac>
<msqrt>
<mrow>
<msup>
<mi>&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>&rho;</mi>
</mfrac>
<mo>)</mo>
</mrow>
</mrow>
In formula:
<mrow>
<mi>&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>&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>&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&theta;</mi>
<mn>1</mn>
</msub>
<msub>
<mi>c&theta;</mi>
<mn>23</mn>
</msub>
</mrow>
</mtd>
<mtd>
<mrow>
<msub>
<mi>s&theta;</mi>
<mn>1</mn>
</msub>
<msub>
<mi>c&theta;</mi>
<mn>23</mn>
</msub>
</mrow>
</mtd>
<mtd>
<mrow>
<mo>-</mo>
<msub>
<mi>s&theta;</mi>
<mn>23</mn>
</msub>
</mrow>
</mtd>
<mtd>
<mrow>
<mo>-</mo>
<msub>
<mi>l</mi>
<mn>2</mn>
</msub>
<msub>
<mi>c&theta;</mi>
<mn>3</mn>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>-</mo>
<msub>
<mi>c&theta;</mi>
<mn>1</mn>
</msub>
<msub>
<mi>s&theta;</mi>
<mn>23</mn>
</msub>
</mrow>
</mtd>
<mtd>
<mrow>
<mo>-</mo>
<msub>
<mi>s&theta;</mi>
<mn>1</mn>
</msub>
<msub>
<mi>s&theta;</mi>
<mn>23</mn>
</msub>
</mrow>
</mtd>
<mtd>
<mrow>
<mo>-</mo>
<msub>
<mi>c&theta;</mi>
<mn>23</mn>
</msub>
</mrow>
</mtd>
<mtd>
<mrow>
<msub>
<mi>l</mi>
<mn>2</mn>
</msub>
<msub>
<mi>s&theta;</mi>
<mn>3</mn>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>-</mo>
<msub>
<mi>s&theta;</mi>
<mn>1</mn>
</msub>
</mrow>
</mtd>
<mtd>
<mrow>
<msub>
<mi>c&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>&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&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&theta;</mi>
<mn>1</mn>
</msub>
<mo>+</mo>
<msub>
<mi>p</mi>
<mi>y</mi>
</msub>
<msub>
<mi>s&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&theta;</mi>
<mn>3</mn>
</msub>
</mrow>
<msqrt>
<mrow>
<msup>
<mrow>
<mo>(</mo>
<mi>p</mi>
<mi>x</mi>
<mi>c</mi>
<mi>&theta;</mi>
<mn>1</mn>
<mo>+</mo>
<mi>p</mi>
<mi>y</mi>
<mi>s</mi>
<mi>&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&theta;</mi>
<mn>1</mn>
</msub>
<mo>+</mo>
<msub>
<mi>p</mi>
<mi>y</mi>
</msub>
<msub>
<mi>s&theta;</mi>
<mn>1</mn>
</msub>
<mo>)</mo>
</mrow>
</mrow>
So, by θ23=θ2+θ3It can obtain:
<mrow>
<msub>
<mi>&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&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&theta;</mi>
<mn>1</mn>
</msub>
<mo>+</mo>
<msub>
<mi>p</mi>
<mi>y</mi>
</msub>
<msub>
<mi>s&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&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&theta;</mi>
<mn>1</mn>
</msub>
<mo>+</mo>
<msub>
<mi>p</mi>
<mi>y</mi>
</msub>
<msub>
<mi>s&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&theta;</mi>
<mn>1</mn>
</msub>
<mo>+</mo>
<msub>
<mi>p</mi>
<mi>y</mi>
</msub>
<msub>
<mi>s&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>&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 θ1cθ23ax+sθ1cθ23ay-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&theta;</mi>
<mn>1</mn>
</msub>
<msub>
<mi>c&theta;</mi>
<mn>23</mn>
</msub>
<msub>
<mi>c&theta;</mi>
<mn>4</mn>
</msub>
<mo>+</mo>
<msub>
<mi>s&theta;</mi>
<mn>1</mn>
</msub>
<msub>
<mi>s&theta;</mi>
<mn>4</mn>
</msub>
</mrow>
</mtd>
<mtd>
<mrow>
<msub>
<mi>s&theta;</mi>
<mn>1</mn>
</msub>
<msub>
<mi>c&theta;</mi>
<mn>23</mn>
</msub>
<msub>
<mi>c&theta;</mi>
<mn>4</mn>
</msub>
<mo>-</mo>
<msub>
<mi>c&theta;</mi>
<mn>1</mn>
</msub>
<msub>
<mi>s&theta;</mi>
<mn>4</mn>
</msub>
</mrow>
</mtd>
<mtd>
<mrow>
<mo>-</mo>
<msub>
<mi>s&theta;</mi>
<mn>23</mn>
</msub>
<msub>
<mi>c&theta;</mi>
<mn>4</mn>
</msub>
</mrow>
</mtd>
<mtd>
<mrow>
<mo>-</mo>
<msub>
<mi>d</mi>
<mn>4</mn>
</msub>
<msub>
<mi>c&theta;</mi>
<mn>1</mn>
</msub>
<msub>
<mi>s&theta;</mi>
<mn>23</mn>
</msub>
<mo>+</mo>
<msub>
<mi>l</mi>
<mn>2</mn>
</msub>
<msub>
<mi>c&theta;</mi>
<mn>1</mn>
</msub>
<msub>
<mi>c&theta;</mi>
<mn>2</mn>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>-</mo>
<msub>
<mi>c&theta;</mi>
<mn>1</mn>
</msub>
<msub>
<mi>c&theta;</mi>
<mn>23</mn>
</msub>
<msub>
<mi>s&theta;</mi>
<mn>4</mn>
</msub>
<mo>+</mo>
<msub>
<mi>s&theta;</mi>
<mn>1</mn>
</msub>
<msub>
<mi>c&theta;</mi>
<mn>4</mn>
</msub>
</mrow>
</mtd>
<mtd>
<mrow>
<mo>-</mo>
<msub>
<mi>s&theta;</mi>
<mn>1</mn>
</msub>
<msub>
<mi>c&theta;</mi>
<mn>23</mn>
</msub>
<msub>
<mi>s&theta;</mi>
<mn>4</mn>
</msub>
<mo>-</mo>
<msub>
<mi>c&theta;</mi>
<mn>1</mn>
</msub>
<msub>
<mi>c&theta;</mi>
<mn>4</mn>
</msub>
</mrow>
</mtd>
<mtd>
<mrow>
<msub>
<mi>s&theta;</mi>
<mn>23</mn>
</msub>
<msub>
<mi>s&theta;</mi>
<mn>4</mn>
</msub>
</mrow>
</mtd>
<mtd>
<mrow>
<mo>-</mo>
<msub>
<mi>d</mi>
<mn>4</mn>
</msub>
<msub>
<mi>s&theta;</mi>
<mn>1</mn>
</msub>
<msub>
<mi>s&theta;</mi>
<mn>23</mn>
</msub>
<mo>+</mo>
<msub>
<mi>l</mi>
<mn>2</mn>
</msub>
<msub>
<mi>s&theta;</mi>
<mn>1</mn>
</msub>
<msub>
<mi>c&theta;</mi>
<mn>2</mn>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>-</mo>
<msub>
<mi>c&theta;</mi>
<mn>1</mn>
</msub>
<msub>
<mi>s&theta;</mi>
<mn>23</mn>
</msub>
</mrow>
</mtd>
<mtd>
<mrow>
<mo>-</mo>
<msub>
<mi>s&theta;</mi>
<mn>1</mn>
</msub>
<msub>
<mi>s&theta;</mi>
<mn>23</mn>
</msub>
</mrow>
</mtd>
<mtd>
<mrow>
<mo>-</mo>
<msub>
<mi>c&theta;</mi>
<mn>23</mn>
</msub>
</mrow>
</mtd>
<mtd>
<mrow>
<mo>-</mo>
<msub>
<mi>d</mi>
<mn>4</mn>
</msub>
<msub>
<mi>c&theta;</mi>
<mn>23</mn>
</msub>
<mo>-</mo>
<msub>
<mi>l</mi>
<mn>2</mn>
</msub>
<msub>
<mi>s&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 θ1cθ23cθ4+sθ1sθ4)ax+(sθ1cθ23cθ4-cθ1sθ4)ay-sθ23cθ4az
M2=-c θ1sθ23ax-sθ1sθ23ay-cθ23az
For the 6th joint angles, it can obtain:
<mrow>
<msub>
<mi>&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 θ1cθ23cθ4cθ5+sθ1sθ4cθ5+cθ1sθ23sθ5
O=s θ1cθ4-cθ1cθ23sθ4
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&theta;</mi>
<mn>1</mn>
</msub>
<msub>
<mi>c&theta;</mi>
<mn>23</mn>
</msub>
</mrow>
</mtd>
<mtd>
<mrow>
<msub>
<mi>s&theta;</mi>
<mn>1</mn>
</msub>
<msub>
<mi>c&theta;</mi>
<mn>23</mn>
</msub>
</mrow>
</mtd>
<mtd>
<mrow>
<mo>-</mo>
<msub>
<mi>s&theta;</mi>
<mn>23</mn>
</msub>
</mrow>
</mtd>
<mtd>
<mrow>
<mo>-</mo>
<msub>
<mi>l</mi>
<mn>2</mn>
</msub>
<msub>
<mi>c&theta;</mi>
<mn>3</mn>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>-</mo>
<msub>
<mi>c&theta;</mi>
<mn>1</mn>
</msub>
<msub>
<mi>s&theta;</mi>
<mn>23</mn>
</msub>
</mrow>
</mtd>
<mtd>
<mrow>
<mo>-</mo>
<msub>
<mi>s&theta;</mi>
<mn>1</mn>
</msub>
<msub>
<mi>s&theta;</mi>
<mn>23</mn>
</msub>
</mrow>
</mtd>
<mtd>
<mrow>
<mo>-</mo>
<msub>
<mi>c&theta;</mi>
<mn>23</mn>
</msub>
</mrow>
</mtd>
<mtd>
<mrow>
<msub>
<mi>l</mi>
<mn>2</mn>
</msub>
<msub>
<mi>s&theta;</mi>
<mn>3</mn>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>-</mo>
<msub>
<mi>s&theta;</mi>
<mn>1</mn>
</msub>
</mrow>
</mtd>
<mtd>
<mrow>
<msub>
<mi>c&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:
-oxsθ1+oycθ1=c θ4sθ6±sθ4cθ6
θ4±θ6=Atan2 (oxsθ1-oycθ1,-oxcθ1cθ23-oysθ1cθ23+ozsθ23)
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.
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)
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)
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 |
-
2017
- 2017-09-28 CN CN201710902038.1A patent/CN107791248B/en active Active
Patent Citations (7)
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)
Title |
---|
于凌涛等: "一类不满足Pieper准则的机器人逆运动学解析解获取方法", 《机器人 ROBOT》 * |
刘松国等: "6R机器人实时逆运动学算法研究", 《控制理论与应用》 * |
唐奥林: "面向主从式微创外科手术机器人的遥操作运动控制策略研究", 《中国博士学位论文全文数据库 信息科技辑》 * |
Cited By (8)
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 |