CN113305809B - Force sense simulation control method of fully-constrained space rope driving parallel mechanism - Google Patents

Force sense simulation control method of fully-constrained space rope driving parallel mechanism Download PDF

Info

Publication number
CN113305809B
CN113305809B CN202110586409.6A CN202110586409A CN113305809B CN 113305809 B CN113305809 B CN 113305809B CN 202110586409 A CN202110586409 A CN 202110586409A CN 113305809 B CN113305809 B CN 113305809B
Authority
CN
China
Prior art keywords
rope
force
vector
platform
guide pulley
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN202110586409.6A
Other languages
Chinese (zh)
Other versions
CN113305809A (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.)
Harbin Institute of Technology Shenzhen
Original Assignee
Harbin Institute of Technology Shenzhen
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 Harbin Institute of Technology Shenzhen filed Critical Harbin Institute of Technology Shenzhen
Priority to CN202110586409.6A priority Critical patent/CN113305809B/en
Publication of CN113305809A publication Critical patent/CN113305809A/en
Application granted granted Critical
Publication of CN113305809B publication Critical patent/CN113305809B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

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
    • BPERFORMING OPERATIONS; TRANSPORTING
    • B25HAND TOOLS; PORTABLE POWER-DRIVEN TOOLS; MANIPULATORS
    • B25JMANIPULATORS; CHAMBERS PROVIDED WITH MANIPULATION DEVICES
    • B25J9/00Programme-controlled manipulators
    • B25J9/003Programme-controlled manipulators having parallel kinematics
    • B25J9/0078Programme-controlled manipulators having parallel kinematics actuated by cables
    • 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/1605Simulation of manipulator lay-out, design, modelling of manipulator

Landscapes

  • Engineering & Computer Science (AREA)
  • Robotics (AREA)
  • Mechanical Engineering (AREA)
  • Automation & Control Theory (AREA)
  • Physics & Mathematics (AREA)
  • Mathematical Physics (AREA)
  • Manipulator (AREA)

Abstract

The invention discloses a force sense simulation control method of a fully constrained space rope driven parallel mechanism, which comprises the following steps: step one, determining a motion calculation method; step two, solving a Jacobian matrix; step three, determining a tension distribution method; step four, establishing a dynamic model of the mechanism when the elasticity of the rope is considered; step five, building a closed-loop motion control frame; step six, determining a force sense simulation control scheme; and step seven, completing the construction of the overall control framework. The control frame of the invention is composed of three parts, which can ensure the good position tracking of the mobile platform, ensure the small and adjustable tension on the rope, and increase the man-machine interaction function of the rope parallel mechanism. In the specific control scheme, the dynamic performance of the mechanism is improved by introducing the internal force coordination ring, and the interaction with people is realized by introducing the force sense simulation ring. The force sensing simulation ring can be used or not according to the use requirement, and the rest two parts still have a good control effect.

Description

一种全约束空间绳索驱动并联机构的力感模拟控制方法A force-sensing simulation control method for a fully constrained space rope-driven parallel mechanism

技术领域technical field

本发明属于控制技术领域,涉及一种全约束空间绳索驱动并联机构的力感模拟控制方法。The invention belongs to the technical field of control, and relates to a force-sensing simulation control method of a fully constrained space rope-driven parallel mechanism.

背景技术Background technique

绳索驱动并联机构具有机构质量小、运动空间大、运动惯量小、刚度可调节、负载能力强等优势,逐渐受到众多学者的关注,而我国在绳索并联机构方面的研究起步较晚,目前可应用于众多领域,例如:射电望远镜、风洞试验、康复机器人、工业相机以及搬运装配等。由于绳索单向受力的特点,机构要实现n自由度的灵活控制,绳索数目m应满足m≥n+1,即实现六个自由度的控制至少需要7根绳索,为了保证机构的对称性,大多采用8根绳索实现全约束空间的控制。在控制绳索驱动并联机构时涉及机构的动力学特性,同时机构为冗余约束,调节绳索内力也至关重要。其次,使机构能够与人进行交互也是当今时代发展的趋势,这对绳索驱动并联机构的控制框架提出了更高的要求,研究绳索驱动并联机构的控制框架十分重要。The rope-driven parallel mechanism has the advantages of small mechanism mass, large motion space, small motion inertia, adjustable stiffness, and strong load capacity. It is used in many fields, such as radio telescopes, wind tunnel tests, rehabilitation robots, industrial cameras, and handling and assembly. Due to the unidirectional force of the rope, the mechanism needs to realize flexible control of n degrees of freedom, and the number m of ropes should satisfy m ≥ n+1, that is, at least 7 ropes are required to realize the control of six degrees of freedom. In order to ensure the symmetry of the mechanism , most of which use 8 ropes to realize the control of the fully constrained space. The dynamic characteristics of the mechanism are involved in the control of the rope-driven parallel mechanism. At the same time, the mechanism is a redundant constraint, and it is also crucial to adjust the internal force of the rope. Secondly, it is also the development trend of the current era to enable the mechanism to interact with people, which puts forward higher requirements for the control framework of the rope-driven parallel mechanism. It is very important to study the control framework of the rope-driven parallel mechanism.

发明内容SUMMARY OF THE INVENTION

为了解决上述问题,本发明提供了一种全约束空间绳索驱动并联机构的力感模拟控制方法。本发明结构简单,运动范围大,绳索具有本质安全的特点,机构一方面能够稳定运行,实现对于给定轨迹的运动;另一方面能够与人进行交互,可以对外界施加的力做出相应的反应,实现了对机构移动平台位置和绳索张力两个被控量的良好控制,同时能够用来模拟人与虚拟物体交互的力感。In order to solve the above problems, the present invention provides a force-sensing simulation control method of a fully constrained space rope-driven parallel mechanism. The invention has a simple structure, a large range of motion, and the rope has the characteristics of intrinsic safety. On the one hand, the mechanism can run stably and realize the movement of a given trajectory; on the other hand, it can interact with people and make corresponding responses to external forces It can realize the good control of the two controlled quantities of the position of the mechanism moving platform and the tension of the rope, and can be used to simulate the force feeling of the interaction between people and virtual objects.

本发明的目的是通过以下技术方案实现的:The purpose of this invention is to realize through the following technical solutions:

一种全约束空间绳索驱动并联机构的力感模拟控制方法,包括如下步骤:A force-sensing simulation control method for a fully constrained space rope-driven parallel mechanism, comprising the following steps:

步骤一、确定运动解算方法:Step 1. Determine the motion solution method:

所述运动解算方法一方面考虑导向滑轮的影响,另一方面求解所述卷扬结构时电机转角和绳长变化量的关系,因此将绳索分成两部分进行分析:工作空间内部分和驱动系统内部分,对于工作空间内部分,求解摆动角βi、包角αi、出绳点坐标BSi,最后得到工作空间部分绳长lW,i;对于驱动系统内部分,通过两个时刻绳长状态得到恒等关系求解电机转角和绳长变化量关系,其中:The motion calculation method considers the influence of the guide pulley on the one hand, and solves the relationship between the motor rotation angle and the change of the rope length when solving the hoisting structure, so the rope is divided into two parts for analysis: the part in the working space and the driving system In the inner part, for the inner part of the working space, solve the swing angle β i , the wrapping angle α i , the coordinates of the rope exit point B S i , and finally get the rope length l W,i of the working space part; for the inner part of the driving system, through two moments The rope length state is obtained by obtaining the identity relationship to solve the relationship between the motor rotation angle and the rope length variation, where:

(1)摆动角βi表示为:(1) The swing angle β i is expressed as:

Figure BDA0003087631480000021
Figure BDA0003087631480000021

式中,Bai,x为动铰点向量在全局坐标系x轴的分量;Bai,y为动铰点向量在全局坐标系y轴的分量;Bbi,x为导向滑轮摆动轴与导向滑轮边缘切点向量在全局坐标系x轴的分量;Bbi,y为导向滑轮摆动轴与导向滑轮边缘切点向量在全局坐标系y轴的分量;In the formula, B a i,x is the component of the moving hinge point vector on the x-axis of the global coordinate system; B a i,y is the component of the moving hinge point vector on the y-axis of the global coordinate system; B b i,x is the swing of the guide pulley The component of the tangent point vector between the axis and the edge of the guide pulley on the x-axis of the global coordinate system; B b i,y is the component of the tangent point vector between the swing axis of the guide pulley and the edge of the guide pulley on the y-axis of the global coordinate system;

(2)包角αi表示为:(2) The wrap angle α i is expressed as:

αi=π-(αi,1i,2);α i =π-(α i,1i,2 );

Figure BDA0003087631480000031
Figure BDA0003087631480000031

Figure BDA0003087631480000032
Figure BDA0003087631480000032

式中,Bmi为导向滑轮中心在全局坐标系中的向量;Bmi,z为导向滑轮中心向量在全局坐标系z轴的分量;Bai为动铰点在全局坐标系中的向量;Bai,z为动铰点向量在全局坐标系z轴的分量;ρ为导向滑轮半径。In the formula, B m i is the vector of the center of the guide pulley in the global coordinate system; B m i,z is the component of the center vector of the guide pulley in the z-axis of the global coordinate system; B a i is the dynamic hinge point in the global coordinate system. vector; B a i,z is the component of the dynamic hinge point vector on the z-axis of the global coordinate system; ρ is the radius of the guide pulley.

(3)出绳点坐标BSi表示为:(3) The coordinate B S i of the rope exit point is expressed as:

Figure BDA0003087631480000033
Figure BDA0003087631480000033

式中,Bmi,x为导向滑轮中心向量在全局坐标系x轴的分量;Bmi,y为导向滑轮中心向量在全局坐标系y轴的分量;In the formula, B m i,x is the component of the center vector of the guide pulley on the x-axis of the global coordinate system; B m i,y is the component of the center vector of the guide pulley on the y-axis of the global coordinate system;

(4)工作空间部分绳长lW,i表示为:(4) The rope length l W,i of the working space part is expressed as:

lW,i=qP+lAS=ρ·αi+||Bai-BSi||2l W,i =q P +l AS =ρ·α i +|| B a i - B S i || 2 ;

式中,qP为导向滑轮上包裹的绳长;lAS为导向滑轮出绳点到动铰点间的绳长;In the formula, q P is the length of the rope wrapped on the guide pulley; l AS is the length of the rope between the rope exit point of the guide pulley and the moving hinge point;

(5)电机转角和绳长变化量关系表示为:(5) The relationship between motor rotation angle and rope length variation is expressed as:

Figure BDA0003087631480000034
Figure BDA0003087631480000034

式中,ΔlW为绳索有效变化量;h为绳索缠绕螺距;

Figure BDA0003087631480000036
为电机转角;vPG为电机和卷扬之间的传动比;r为卷扬半径;In the formula, Δl W is the effective variation of the rope; h is the rope winding pitch;
Figure BDA0003087631480000036
is the motor rotation angle; v PG is the transmission ratio between the motor and the winch; r is the winch radius;

步骤二、求解雅可比矩阵:Step 2. Solve the Jacobian matrix:

Figure BDA0003087631480000035
Figure BDA0003087631480000035

式中,J为机构的雅可比矩阵;BRP为全局坐标系和体坐标系之间位姿变化的旋转矩阵;ln为单位绳长向量;Pai为动铰点在体坐标系中的向量,BaiBRP PaiIn the formula, J is the Jacobian matrix of the mechanism; B R P is the rotation matrix of the pose change between the global coordinate system and the body coordinate system; ln is the unit rope length vector; P a i is the dynamic hinge point in the body coordinate system. A vector in , B a i = B R P P a i ;

步骤三、确定张力分配方法:Step 3. Determine the tension distribution method:

采用Closed-form理论对绳索进行张力分配:Tension distribution of the rope using the Closed-form theory:

fis=-J+Tw+fm-J+TJTfmf is =-J +T w+f m -J +T J T f m ;

式中,fis为绳索张力;J+T为雅可比矩阵伪逆;w为外力;fm为张力范围最大最小值的平均值;In the formula, f is the rope tension; J + T is the pseudo-inverse of the Jacobian matrix; w is the external force; f m is the average value of the maximum and minimum values of the tension range;

步骤四、建立考虑绳索弹性时机构的动力学模型:Step 4. Establish the dynamic model of the mechanism when considering the elasticity of the rope:

假设绳索与移动平台、绳索与导向滑轮、绳索与卷扬之间的连接均为理想状态,基于以上假设,动力学建模包括移动平台、卷扬和导向滑轮以及绳索,其中:Assuming that the connections between the rope and the mobile platform, the rope and the guide pulley, and the rope and the hoist are all ideal, based on the above assumptions, the dynamic modeling includes the mobile platform, the hoist and the guide pulley, and the rope, where:

(1)采用牛顿欧拉法建模移动平台:(1) Using Newton's Euler method to model the mobile platform:

Figure BDA0003087631480000041
Figure BDA0003087631480000041

式中,L为各绳的单位方向向量所构成的矩阵;PA为动平台铰点坐标矩阵;E为3×3单位矩阵;I为动平台相对于定平台的惯量矩阵;Ω为反对称矩阵;m为移动平台质量;g为重力加速度向量;f2i为平台侧绳索张力;t为动平台相对于定平台的移动位移;ω和

Figure BDA0003087631480000043
为动平台相对于定平台的转动角速度和角加速度;In the formula, L is the matrix formed by the unit direction vector of each rope; P A is the coordinate matrix of the hinge point of the moving platform; E is the 3×3 unit matrix; I is the inertia matrix of the moving platform relative to the fixed platform; Ω is the antisymmetric matrix; m is the mass of the moving platform; g is the gravitational acceleration vector; f 2i is the rope tension on the platform side; t is the displacement of the moving platform relative to the fixed platform; ω and
Figure BDA0003087631480000043
is the rotational angular velocity and angular acceleration of the moving platform relative to the fixed platform;

(2)采用力矩平衡法建模卷扬和导向滑轮:(2) Using the moment balance method to model the hoist and guide pulley:

Figure BDA0003087631480000042
Figure BDA0003087631480000042

式中,Jr为卷扬和滑轮模型的等效转动惯量;

Figure BDA0003087631480000044
为第i个电机的转角矢量;Cr为卷扬和滑轮模型的等效阻尼系数矩阵;λs为卷扬和滑轮模型的等效半径;f1i为第i根绳索驱动侧的绳索拉力;Tact,i为驱动器的转矩;where J r is the equivalent moment of inertia of the hoist and pulley models;
Figure BDA0003087631480000044
is the rotation angle vector of the ith motor; C r is the equivalent damping coefficient matrix of the hoist and pulley models; λ s is the equivalent radius of the hoist and pulley models; f 1i is the rope tension on the driving side of the ith rope; T act,i is the torque of the driver;

(3)将绳索建模为弹簧阻尼模型;(3) Model the rope as a spring damping model;

Figure BDA0003087631480000051
Figure BDA0003087631480000051

式中,ci为绳索的刚度,di为绳索的阻尼;Δli为绳索形变量;

Figure BDA0003087631480000053
为绳索形变速度;In the formula, ci is the stiffness of the rope, d i is the damping of the rope; Δl i is the deformation of the rope;
Figure BDA0003087631480000053
is the rope deformation speed;

(4)动力学数学模型为:(4) The dynamic mathematical model is:

Figure BDA0003087631480000052
Figure BDA0003087631480000052

式中,F为作用在移动平台上的外力;In the formula, F is the external force acting on the mobile platform;

步骤五、搭建闭环运动控制框架:Step 5. Build a closed-loop motion control framework:

所述闭环运动控制框架包括基于关节反馈的运动控制框架和基于位姿反馈的运动控制框架,其中:The closed-loop motion control framework includes a joint feedback-based motion control framework and a pose feedback-based motion control framework, wherein:

所述基于关节反馈的运动控制框架和基于位姿反馈的运动控制框架均包括期望位姿、移动平台位姿反馈、卷扬转角反馈、内力协调、电机驱动系统、带有弹性绳索的动力学模型;The joint feedback-based motion control framework and the pose feedback-based motion control framework both include desired pose, mobile platform pose feedback, hoisting angle feedback, internal force coordination, motor drive system, and a dynamic model with elastic ropes ;

所述基于关节反馈的运动控制框架中,由期望位姿通过运动学逆解解算得到期望的各绳绳长变化量,将结果转换为期望电机转角,再与反馈的实际转角比较得到关节偏差,输入到PID控制器中实现位置控制;In the motion control framework based on joint feedback, the expected change in rope length of each rope is obtained through the inverse kinematics solution from the desired pose, the result is converted into the expected motor rotation angle, and then the joint deviation is obtained by comparing with the actual rotation angle of the feedback. , input into the PID controller to realize position control;

所述基于位姿反馈的运动控制框架中,由期望位姿和模型测量位姿比较得到位姿误差,根据绳索长度误差和位姿误差之间相差一个雅可比矩阵可以得到绳索长度误差,实现由任务空间误差到关节空间误差的转换,进而求得电机的转角误差并输入到PD控制器中实现位置控制;In the motion control framework based on pose feedback, the pose error is obtained by comparing the desired pose and the model measured pose, and the rope length error can be obtained according to the difference between the rope length error and the pose error by a Jacobian matrix. The conversion of the task space error to the joint space error, and then the angle error of the motor is obtained and input into the PD controller to realize the position control;

所述内力协调环中,根据张力分配方法求解期望张力,与模型测量得到的绳索张力构成力内环反馈。In the inner force coordination loop, the expected tension is solved according to the tension distribution method, and the force inner loop feedback is formed with the rope tension measured by the model.

步骤六、确定力感模拟控制方案:Step 6: Determine the force-sensing simulation control scheme:

采用力感模拟环将整个机构等效为质量-弹簧-阻尼系统,从而得到力和运动间的传递函数,使外力转换为平台的位置修正量,通过设置不同的参数组合改变机构的物理特性,即机构的移动平台在外力F作用下产生不同的运动状态输出,其中:The whole mechanism is equivalent to a mass-spring-damping system by using a force-sensing simulation ring, so as to obtain the transfer function between force and motion, so that the external force can be converted into the position correction value of the platform, and the physical characteristics of the mechanism can be changed by setting different parameter combinations. That is, the mobile platform of the mechanism produces different motion state outputs under the action of the external force F, among which:

力和运动间的传递函数表示为:The transfer function between force and motion is expressed as:

Figure BDA0003087631480000061
Figure BDA0003087631480000061

式中,Md为等效质量;Dd为等效阻尼;Kd为等效刚度;F为作用在移动平台上的外力;E为移动平台的位移修正量;where M d is the equivalent mass; D d is the equivalent damping; K d is the equivalent stiffness; F is the external force acting on the mobile platform; E is the displacement correction of the mobile platform;

步骤七、完成整体控制框架搭建:Step 7. Complete the construction of the overall control framework:

所述整体控制框架至少包括两个部分:移动平台位置闭环(两种框架)、内力协调环;如果并联机构需要实现力感模拟,则整体控制框架还包括力感模拟环,执行步骤八;The overall control framework includes at least two parts: a closed-loop position of the mobile platform (two kinds of frameworks) and an internal force coordination loop; if the parallel mechanism needs to realize force-sensing simulation, the overall control framework also includes a force-sensing simulation loop, and step 8 is performed;

步骤八、对不同的力感进行模拟,验证人机交互的可行性。Step 8: Simulate different force senses to verify the feasibility of human-computer interaction.

相比于现有技术,本发明具有如下优点:Compared with the prior art, the present invention has the following advantages:

本发明的控制框架采用三个部分组成,能够保证移动平台的位置跟踪良好,也能够保证绳上张力对抗小且可调,同时增加了绳索并联机构的人机交互功能。具体控制方案中通过引入内力协调环提高机构的动态性能,通过引入力感模拟环实现了与人的交互。其中,力感模拟环可以根据使用需要选择是否采用,其余两部分仍然具备较好的控制效果。The control frame of the present invention is composed of three parts, which can ensure the position tracking of the mobile platform is good, and can also ensure that the tension on the rope is small and adjustable, and at the same time, the human-computer interaction function of the rope parallel mechanism is increased. In the specific control scheme, the dynamic performance of the mechanism is improved by introducing an internal force coordination loop, and the interaction with humans is realized by introducing a force-sensing simulation loop. Among them, the force-sensing simulation ring can be selected according to the needs of use, and the other two parts still have a good control effect.

附图说明Description of drawings

图1为整体机构简图;Figure 1 is a schematic diagram of the overall mechanism;

图2为图1的局部放大图;Fig. 2 is a partial enlarged view of Fig. 1;

图3为工作空间内绳索封闭向量示意图;Figure 3 is a schematic diagram of a rope closure vector in the workspace;

图4为绳索传动示意图;Figure 4 is a schematic diagram of a rope drive;

图5为卷扬机构安装示意图;Figure 5 is a schematic diagram of the installation of the hoisting mechanism;

图6为动力学模型框图;Figure 6 is a block diagram of a dynamic model;

图7为基于位姿反馈的整体控制框架;Figure 7 shows the overall control framework based on pose feedback;

图8为基于关节反馈的整体控制框架。Figure 8 shows the overall control framework based on joint feedback.

具体实施方式Detailed ways

下面结合附图对本发明的技术方案作进一步的说明,但并不局限于此,凡是对本发明技术方案进行修改或者等同替换,而不脱离本发明技术方案的精神和范围,均应涵盖在本发明的保护范围中。The technical solutions of the present invention will be further described below in conjunction with the accompanying drawings, but are not limited thereto. Any modification or equivalent replacement of the technical solutions of the present invention without departing from the spirit and scope of the technical solutions of the present invention shall be included in the present invention. within the scope of protection.

本发明为八根绳索驱动的并联机构,并联机构通过电机连接卷扬,卷扬上缠绕绳索,绳索通过导向滑轮导向最终连接在末端执行机构即移动平台上,通过控制电机的转角完成末端执行机构的运动控制,所述并联机构的移动平台上安装有手柄,方便人与机构进行交互。图1和图2为整体机构简图,人们可以通过手握手柄移动平台与机构进行交互。图4为一根绳索完整的传动简图,图5所示为卷扬机构的安装示意图。具体步骤如下:The invention is a parallel mechanism driven by eight ropes. The parallel mechanism is connected to the hoist through the motor, the rope is wound on the hoist, and the rope is guided by the guide pulley and finally connected to the end effector, that is, the mobile platform, and the end effector is completed by controlling the rotation angle of the motor. A handle is installed on the mobile platform of the parallel mechanism, which is convenient for people to interact with the mechanism. Figures 1 and 2 are schematic diagrams of the overall mechanism, and people can interact with the mechanism by holding the handle to move the platform. Figure 4 is a schematic diagram of a complete transmission of a rope, and Figure 5 is a schematic diagram of the installation of the hoisting mechanism. Specific steps are as follows:

步骤1:确定运动解算方法。Step 1: Determine the motion solution method.

并联机构的导向滑轮会对绳长产生影响,同时并联机构的传动结构会造成卷扬松开的绳长和末端释放的绳长不等,因此本发明一方面将考虑导向滑轮的影响,另一方面求解该结构时电机转角和绳长变化量的关系。The guide pulley of the parallel mechanism will affect the length of the rope, and the transmission structure of the parallel mechanism will cause the length of the rope loosened by the hoist and the length of the rope released at the end to be different. Therefore, the present invention will consider the influence of the guide pulley on the one hand, and the other The relationship between the motor rotation angle and the rope length variation when the structure is solved.

将图4所示的绳索传动示意图分成两部分进行分析:工作空间内部分(卷扬出绳点C-导向滑轮摆动轴与导向滑轮边缘的切点B)和驱动系统内部分(导向滑轮摆动轴与导向滑轮边缘的切点B-动铰点A)。The schematic diagram of the rope transmission shown in Figure 4 is divided into two parts for analysis: the inner part of the working space (the hoisting rope point C - the tangent point B between the swing shaft of the guide pulley and the edge of the guide pulley) and the inner part of the drive system (the swing shaft of the guide pulley) Tangent point B to the edge of the guide pulley - moving hinge point A).

一、工作空间内部分:如图3所示,已知工作空间的滑轮半径为ρ,考虑滑轮对绳长的影响,将工作空间内绳索长度分解为两部分来计算:第一部分为包裹在滑轮上的长度,即绳索与滑轮两个切点之间的弧长,第二部分为岀绳点与动平台连接点之间的长度。1. The inner part of the working space: As shown in Figure 3, the radius of the pulley in the known working space is ρ. Considering the influence of the pulley on the rope length, the length of the rope in the working space is decomposed into two parts to calculate: the first part is wrapped in the pulley. The length above is the arc length between the two tangent points of the rope and the pulley, and the second part is the length between the rope point and the connection point of the moving platform.

(1)求解摆动角βi(1) Solve the swing angle β i :

Figure BDA0003087631480000081
Figure BDA0003087631480000081

式中,Bai,x为动铰点向量在全局坐标系x轴的分量;Bai,y为动铰点向量在全局坐标系y轴的分量;Bbi,x为导向滑轮摆动轴与导向滑轮边缘切点向量在全局坐标系x轴的分量;Bbi,y为导向滑轮摆动轴与导向滑轮边缘切点向量在全局坐标系y轴的分量。In the formula, B a i,x is the component of the moving hinge point vector on the x-axis of the global coordinate system; B a i,y is the component of the moving hinge point vector on the y-axis of the global coordinate system; B b i,x is the swing of the guide pulley The component of the tangent point vector between the axis and the edge of the guide pulley on the x-axis of the global coordinate system; B b i,y is the component of the tangent point vector between the swing axis of the guide pulley and the edge of the guide pulley on the y-axis of the global coordinate system.

(2)求解包角αi(2) Solve the wrap angle α i :

Figure BDA0003087631480000091
Figure BDA0003087631480000091

式中,

Figure BDA0003087631480000098
为导向滑轮关于摆动角的旋转矩阵;ex为H坐标系中沿x轴正向的单位向量,ex=[1 0 0]T。In the formula,
Figure BDA0003087631480000098
is the rotation matrix of the guide pulley about the swing angle; e x is the unit vector along the positive x-axis in the H coordinate system, e x =[1 0 0] T .

则包角αi可表示为αi=π-(αi,1i,2)。Then the wrap angle α i can be expressed as α i =π-(α i,1i,2 ).

其中,

Figure BDA0003087631480000092
in,
Figure BDA0003087631480000092

(3)求解出绳点坐标BSi(3) Solve the rope point coordinates B S i :

Figure BDA0003087631480000093
Figure BDA0003087631480000093

(4)求解工作空间部分绳长lW,i(4) Solve the part of the rope length l W,i in the workspace:

将两部分长度相加得,lW,i=qP+lAS=ρ·αi+||Bai-BSi||2Adding the lengths of the two parts, l W,i =q P +l AS =ρ·α i +|| B a i - B S i || 2 .

式中,qP为导向滑轮上包裹的绳长;lAS为导向滑轮出绳点到动铰点间的绳长。In the formula, q P is the length of the rope wrapped on the guide pulley; l AS is the length of the rope between the rope exit point of the guide pulley and the moving hinge point.

二、驱动系统内部分:如图4所示,s为驱动系统内滑轮与工作空间内导向滑轮间的直线绳长,可由卷扬转角

Figure BDA0003087631480000097
和绳索缠绕螺距h求得:2. Internal part of the drive system: As shown in Figure 4, s is the length of the straight rope between the pulley in the drive system and the guide pulley in the working space, which can be adjusted by the hoisting angle
Figure BDA0003087631480000097
And the rope winding pitch h is obtained:

Figure BDA0003087631480000094
Figure BDA0003087631480000094

卷扬松开的绳长为:The length of the rope loosened by the winch is:

Figure BDA0003087631480000095
Figure BDA0003087631480000095

该传动形式时松开后驱动系统内部长度

Figure BDA0003087631480000096
加绳索有效变化量ΔlW等于松开前驱动系统内部长度加即将松开的长度,如下式:In this transmission form, the internal length of the rear drive system is released
Figure BDA0003087631480000096
The effective change Δl W of the rope is equal to the internal length of the drive system before loosening plus the length to be loosened, as follows:

Figure BDA0003087631480000101
Figure BDA0003087631480000101

式中,cA和cB见图4所示(cA为卷扬出绳点到绳索与滑轮γ左侧切点间的绳长,cB为包裹在滑轮γ上的绳索长度);lW0,i为工作空间部分绳索原长。In the formula, c A and c B are shown in Figure 4 (c A is the rope length from the point of hoisting the rope to the tangent point between the rope and the left side of the pulley γ, c B is the length of the rope wrapped on the pulley γ); l W0,i is the original length of the rope in the workspace.

整理得:Arranged:

Figure BDA0003087631480000102
Figure BDA0003087631480000102

步骤2:求解雅可比矩阵。Step 2: Solve the Jacobian matrix.

雅可比矩阵可以用来描述机构关节变量速度和移动平台速度之间的关系,同时也描述了关节力和末端移动平台受力的关系,因此求解雅可比矩阵为绳索张力分配和关节变量到末端变量转换奠定基础。The Jacobian matrix can be used to describe the relationship between the speed of the joint variable of the mechanism and the speed of the moving platform, as well as the relationship between the joint force and the force of the end moving platform, so solve the Jacobian matrix to assign the rope tension and the joint variable to the end variable Conversion lays the foundation.

对式

Figure BDA0003087631480000103
两边同时求导并化简得:
Figure BDA0003087631480000105
pair
Figure BDA0003087631480000103
Differentiate both sides and simplify to get:
Figure BDA0003087631480000105

忽略影响较小的项,对上式化简得:Ignoring less influential terms, the above formula can be simplified to:

Figure BDA0003087631480000106
Figure BDA0003087631480000106

步骤3:确定张力分配方法。Step 3: Determine the tension distribution method.

所述并联机构上的绳索具有弹性和单向受力的特点,在移动平台运动时应保证绳索始终处于张紧状态,对于全约束机构来说应对绳索进行张力分配来保证绳上张力在期望范围内。The rope on the parallel mechanism has the characteristics of elasticity and one-way force. When the mobile platform moves, the rope should always be in a tensioned state. For a full restraint mechanism, the rope should be tensioned to ensure that the tension on the rope is within the desired range. Inside.

采用Closed-form理论对绳索进行张力分配:Tension distribution of the rope using the Closed-form theory:

f=-J+Tw+fm-J+TJTfmf=-J + Tw + fm -J + TJTfm .

式中,f为绳索张力;J+T为雅可比矩阵伪逆;w为外力;fm为张力范围最大最小值的平均值。In the formula, f is the rope tension; J + T is the pseudo-inverse of the Jacobian matrix; w is the external force; f m is the average value of the maximum and minimum values of the tension range.

步骤4:建立考虑绳索弹性时机构的动力学模型。Step 4: Build a dynamic model of the mechanism considering the elasticity of the rope.

本发明所研究的绳索驱动并联机构不需要承受较大的负载,移动平台的质量较轻。假设绳索与移动平台、绳索与导向滑轮、绳索与卷扬之间的连接均为理想状态,即机构中不产生摩擦力。基于以上假设,动力学建模包括移动平台、卷扬和导向滑轮以及绳索。The rope-driven parallel mechanism studied in the present invention does not need to bear a large load, and the mass of the mobile platform is relatively light. It is assumed that the connections between the rope and the moving platform, the rope and the guide pulley, and the rope and the hoist are all ideal, that is, there is no friction in the mechanism. Based on the above assumptions, the dynamic modeling includes the moving platform, hoisting and guiding pulleys, and ropes.

移动平台动力学:基于牛顿定律建立末端平台平动过程的动力学模型,基于欧拉方程建立移动平台转动过程的动力学模型。Dynamics of mobile platform: establish the dynamic model of the translation process of the end platform based on Newton's law, and establish the dynamic model of the rotation process of the mobile platform based on Euler's equation.

Figure BDA0003087631480000111
Figure BDA0003087631480000111

Figure BDA0003087631480000112
Figure BDA0003087631480000112

整理得:Arranged:

Figure BDA0003087631480000113
Figure BDA0003087631480000113

卷扬和导向滑轮动力学:假设Jr为卷扬和滑轮模型的等效转动惯量;

Figure BDA0003087631480000114
第i个电机的转角矢量,规定使绳索伸长的方向为正方向;Cr为卷扬和滑轮模型的等效阻尼系数矩阵;λs为卷扬和滑轮模型的等效半径;f1i为第i根绳索驱动侧的绳索拉力;Tact,i为驱动器的转矩。则由力矩平衡得:Hoist and guide sheave dynamics: Assume J r to be the equivalent moment of inertia of the hoist and sheave models;
Figure BDA0003087631480000114
The rotation angle vector of the i-th motor specifies that the direction of rope elongation is the positive direction; C r is the equivalent damping coefficient matrix of the hoist and pulley models; λ s is the equivalent radius of the hoist and pulley models; f 1i is the The rope tension on the drive side of the i-th rope; T act,i is the torque of the drive. Then from the moment balance:

Figure BDA0003087631480000115
Figure BDA0003087631480000115

绳索动力学:将绳索等效为一个弹簧阻尼模型,设定绳索的刚度为ci和绳索的阻尼为di,则有:Rope dynamics: The rope is equivalent to a spring damping model, and the stiffness of the rope is set as ci and the damping of the rope as d i , there are:

Figure BDA0003087631480000121
Figure BDA0003087631480000121

完整动力学:当机构与用户进行交互时,假设用户施加的作用力为F,此时整体的动力学数学模型为,Complete dynamics: When the mechanism interacts with the user, assuming that the force exerted by the user is F, the overall dynamic mathematical model at this time is,

Figure BDA0003087631480000122
Figure BDA0003087631480000122

搭建的动力学仿真框架如图6所示。The dynamic simulation framework built is shown in Figure 6.

步骤5:搭建闭环运动控制框架。Step 5: Build a closed-loop motion control framework.

本发明提供两种控制框架:一种基于位姿反馈的运动控制框架,一种是基于关节反馈的运动控制框架。所述两种框架包括期望位姿、移动平台位姿反馈、卷扬转角反馈、内力协调、电机驱动系统、带有弹性绳索的动力学模型,其中:The present invention provides two control frameworks: a motion control framework based on position and attitude feedback, and a motion control framework based on joint feedback. The two frameworks include desired pose, mobile platform pose feedback, hoisting angle feedback, internal force coordination, motor drive system, and dynamic models with elastic ropes, wherein:

关节空间控制方法由期望位姿通过运动学逆解解算得到期望的各绳绳长变化量,将结果转换为期望电机转角,再与反馈的实际转角比较得到关节偏差,输入到PID控制器中实现位置控制;The joint space control method obtains the expected change of each rope length from the expected pose and kinematic inverse solution, converts the result into the expected motor rotation angle, and then compares it with the feedback actual rotation angle to obtain the joint deviation, which is input into the PID controller. Realize position control;

位姿闭环为任务空间控制方法,由期望位姿和模型测量位姿比较得到位姿误差,根据绳索长度误差和位姿误差之间相差一个雅可比矩阵可以得到绳索长度误差,实现由任务空间误差到关节空间误差的转换,进而求得电机的转角误差并输入到PD控制器中实现位置控制;The pose closed loop is a task space control method. The pose error is obtained by comparing the desired pose and the model measured pose. According to the difference between the rope length error and the pose error by a Jacobian matrix, the rope length error can be obtained. The conversion to the joint space error, and then the angle error of the motor is obtained and input into the PD controller to realize the position control;

控制框架中还包括绳索张力控制方法,即内力协调环,根据张力分配方法求解期望张力,与模型测量得到的绳索张力构成力内环反馈。The control framework also includes a rope tension control method, that is, an internal force coordination loop, which solves the expected tension according to the tension distribution method, and forms a force inner loop feedback with the rope tension measured by the model.

本发明中,基于关节反馈的运动控制框架如图7所示,期望位姿通过运动学逆解解算得到期望的各绳绳长变化量Δl,根据步骤1得到的结果将其转换为期望电机转角,再与反馈的实际转角比较得到关节偏差Δθ,输入到PID控制器中,张力控制采用步骤3所述的张力分配方法构成力内环反馈。In the present invention, the motion control framework based on joint feedback is shown in Figure 7. The desired pose is obtained through the inverse kinematics solution to obtain the desired change Δl of each rope length, and is converted into the desired motor according to the result obtained in step 1. The rotation angle is compared with the actual rotation angle of the feedback to obtain the joint deviation Δθ, which is input into the PID controller, and the tension control adopts the tension distribution method described in step 3 to form a force inner loop feedback.

本发明中,基于位姿反馈的运动控制框架如图8所示,由期望位姿和测量位姿比较可以得到位姿偏差Δq,根据绳索长度误差和位姿误差之间相差一个雅可比矩阵可以得到绳索长度误差Δl=J·Δq,完成了由任务空间误差到关节空间误差的转换,进而可以求得电机的转角误差并输入到PD控制器中实现位置控制,张力控制采用步骤3所述的张力分配方法构成力内环反馈。In the present invention, the motion control framework based on pose feedback is shown in Figure 8. The pose deviation Δq can be obtained by comparing the expected pose and the measured pose. According to the difference between the rope length error and the pose error, a Jacobian matrix can be The rope length error Δl=J·Δq is obtained, and the conversion from the task space error to the joint space error is completed, and then the angle error of the motor can be obtained and input into the PD controller to realize the position control, and the tension control adopts the step 3. The tension distribution method constitutes an inner loop feedback of force.

本发明中,上述两种控制均通过电机驱动系统输出驱动力矩控制机构运动。In the present invention, the above two kinds of control both control the movement of the mechanism by outputting the driving torque of the motor driving system.

步骤6:确定力感模拟控制方案。Step 6: Determine the force-sensing simulation control scheme.

在实现运动控制基础上,本发明还提供了一种用于人机交互的控制方案,采用力感模拟环将整个机构等效为质量-弹簧-阻尼系统,通过设置不同的参数组合来改变机构的物理特性,即机构的移动平台可以在外力F作用下产生不同的运动状态输出。力和运动间的传递函数为:On the basis of realizing motion control, the present invention also provides a control scheme for human-computer interaction. The whole mechanism is equivalent to a mass-spring-damping system by using a force-sensing simulation ring, and the mechanism is changed by setting different parameter combinations. The physical characteristics of the mechanism, that is, the mobile platform of the mechanism can produce different motion state outputs under the action of the external force F. The transfer function between force and motion is:

Figure BDA0003087631480000131
Figure BDA0003087631480000131

能够使外力转换为平台的位置修正量。The external force can be converted into the position correction amount of the platform.

本发明中,采用的是绳索张力传感器,当移动平台受到外力作用后,此时必然会导致绳上张力发生变化,通过雅可比矩阵和受力分析将测量的绳张力转换为外力输入到力感模拟环的控制器中。In the present invention, a rope tension sensor is used. When the mobile platform is subjected to external force, the tension on the rope will inevitably change. Through Jacobian matrix and force analysis, the measured rope tension is converted into external force and input to the force sense. in the controller of the analog loop.

除绳索张力之外的力的合力wis可以表示为:The resultant force w is of the forces other than the rope tension can be expressed as:

wis=-JTfisw is =-J T f is .

式中,fis为绳索张力。where f is the rope tension.

经过受力分析抵消掉惯性力

Figure BDA0003087631480000141
以及重力winit,外力wadm表示为:After force analysis, the inertial force is canceled
Figure BDA0003087631480000141
and the gravity w init , the external force w adm is expressed as:

Figure BDA0003087631480000142
Figure BDA0003087631480000142

步骤7:完成整体控制框架搭建。Step 7: Complete the construction of the overall control framework.

整体控制框架包括移动平台位置闭环(两种框架)、内力协调环、力感模拟环三个部分,如图7和图8所示。如果机构不需要实现力感模拟,只需要去掉控制框架中的力感模拟环,剩余部分仍然可以实现对机构的良好控制。The overall control framework includes three parts of the mobile platform position closed loop (two kinds of frameworks), the internal force coordination loop, and the force sense simulation loop, as shown in Figure 7 and Figure 8. If the mechanism does not need to achieve force-sensing simulation, it is only necessary to remove the force-sensing simulation ring in the control frame, and the remaining parts can still achieve good control of the mechanism.

步骤8:对不同的力感进行模拟,验证人机交互的可行性。Step 8: Simulate different force senses to verify the feasibility of human-computer interaction.

通过仿真模拟几种参数设置情形,可以类比为人与现实中的某种物体进行交互时的力感受,例如:可以模拟推墙壁、推箱子、拍皮球等情形的力感。具体模拟情形如下:By simulating several parameter setting situations, it can be compared to the force feeling when a person interacts with a certain object in reality, for example, it can simulate the force feeling of pushing a wall, pushing a box, shooting a ball, etc. The specific simulation situation is as follows:

当刚度、阻尼和惯量三者均存在时,移动平台可以在一个瞬时推力作用下呈现衰减震荡最后稳定在平衡位置;当阻尼为0时,移动平台可以呈现弹性的效果,人们可以感受到像触摸弹簧时的力感受,此时平台可以用来模拟皮球;当刚度和惯量均为0时,移动平台可以呈现对于外力的跟随运动,当在某一时刻撤去外力,移动平台可以停止在力消失的位置,人们可以感受到像将一个质量很轻的物体从一个位置推到另一个位置时的力感受,此时平台可以用来模拟箱子;当刚度为0时,移动平台在一个瞬时推力的作用下会逐渐走远直到停止,人们可以感受到像推远水面上小船时的力感受,此时平台可以用来模拟水面上的小船。当在某个位置处开始设置无限大的刚度和阻尼,此时在人们持续增加外力的作用下,平台也不会被移动,人们可以感受到像碰撞到墙面时的力感受,此时平台可以用来模拟墙壁。When the stiffness, damping and inertia all exist, the mobile platform can present a damped oscillation under the action of an instantaneous thrust and finally stabilize at the equilibrium position; when the damping is 0, the mobile platform can exhibit an elastic effect, and people can feel like a touch The force feeling of the spring, the platform can be used to simulate a ball; when the stiffness and inertia are both 0, the mobile platform can follow the external force, and when the external force is removed at a certain time, the mobile platform can stop at the point where the force disappears. position, people can feel the force feeling like pushing a very light object from one position to another, at this time the platform can be used to simulate the box; when the stiffness is 0, the moving platform acts on an instantaneous thrust The bottom will gradually go further until it stops, and people can feel the force feeling like pushing a boat on the water away. At this time, the platform can be used to simulate the boat on the water. When the infinite stiffness and damping are set at a certain position, the platform will not be moved under the action of people continuously increasing the external force, and people can feel the force when it hits the wall. At this time, the platform will not be moved. Can be used to simulate walls.

Claims (6)

1. A force sense simulation control method of a rope-driven parallel mechanism in a fully-constrained space is characterized by comprising the following steps of:
step one, determining a motion calculation method:
the rope was divided into two parts for analysis: the working space inner part and the driving system inner part solve the swing angle beta for the working space inner partiAngle of wrap alphaiCoordinates of rope outlet pointBSiFinally, the length l of the rope of the working space part is obtainedW,i(ii) a For the part in the driving system, obtaining an identity relation through the rope length states at two moments to solve the relation between the motor rotation angle and the rope length variation;
step two, solving a Jacobian matrix:
the jacobian matrix is represented as:
J=[ln T (BRP Pai×ln)T];
wherein J is a Jacobian matrix of the mechanism;BRPa rotation matrix of pose change between a global coordinate system and a body coordinate system; lnIs a unit rope length vector;Paiis the vector of the movable hinge point in the body coordinate system,BaiBRP Pai
step three, determining a tension distribution method:
the ropes were tension distributed using Closed-form theory:
fis=-J+Tw+fm-J+TJTfm
in the formula (f)isIs the rope tension; j. the design is a square+TIs the pseudo-inverse of the Jacobian matrix; w is an external force; f. ofmThe average value of the maximum value and the minimum value of the tension range;
step four, establishing a dynamic model of the mechanism in consideration of the elasticity of the rope:
assuming that the connections between the ropes and the moving platform, between the ropes and the guide pulleys, and between the ropes and the winches are all in an ideal state, the mathematical model of the dynamics is expressed as:
Figure FDA0003526729300000021
wherein L is a matrix formed by unit direction vectors of each rope;Pa is a movable platform hinge point coordinate matrix; c. CiIs the stiffness of the rope; diDamping for the rope; Δ liIs the deformation amount of the rope;
Figure FDA0003526729300000023
the rope deformation speed; e is a 3 × 3 identity matrix; i is an inertia matrix of the movable platform relative to the fixed platform; omega is an antisymmetric matrix; m is the mobile platform mass; g is a gravity acceleration vector; t is the moving displacement of the movable platform relative to the fixed platform; ω and
Figure FDA0003526729300000022
the rotating angular velocity and the angular acceleration of the movable platform relative to the fixed platform are obtained; f is an external force acting on the mobile platform;
step five, building a closed-loop motion control framework:
the closed-loop motion control framework comprises a joint feedback-based motion control framework and a pose feedback-based motion control framework, wherein: the motion control frame based on joint feedback and the motion control frame based on pose feedback both comprise an expected pose, moving platform pose feedback, winch corner feedback, internal force coordination, a motor driving system and a dynamic model with an elastic rope;
step six, determining a force sense simulation control scheme:
the whole mechanism is equivalent to a mass-spring-damping system by adopting a force sense simulation ring, so that a transfer function between force and motion is obtained, external force is converted into position correction of a platform, and the physical characteristics of the mechanism are changed by setting different parameter combinations, namely the moving platform of the mechanism is subjected to external force FextDifferent motion state outputs are generated under the action;
and seventhly, completing the construction of an integral control framework:
the overall control frame comprises at least two parts: a mobile platform position closed loop and an internal force coordination loop.
2. The method for simulating and controlling the force sense of a rope-driven parallel mechanism in a fully constrained space as claimed in claim 1, wherein in the first step, the swing angle β isiAngle of wrap alphaiCoordinates of rope outlet pointBSiLength of rope in working spaceW,iAnd the relation between the motor rotation angle and the rope length variation is expressed as follows:
(1) angle of oscillation betai
Figure FDA0003526729300000031
In the formula,Bai,xthe component of the dynamic hinge point vector on the x axis of the global coordinate system is shown;Bai,ythe component of the dynamic hinge point vector on the y axis of the global coordinate system is shown;Bbi,xthe component of the vector of the tangent point of the swing axis of the guide pulley and the edge of the guide pulley on the x axis of the global coordinate system is shown;Bbi,ythe component of the vector of the oscillating axis of the guide pulley and the tangent point of the edge of the guide pulley on the y axis of the global coordinate system;
(2) wrap angle alphai
αi=π-(αi,1i,2);
Figure FDA0003526729300000032
Figure FDA0003526729300000033
In the formula,Bmithe vector of the center of the guide pulley in the global coordinate system is taken as the vector;Bmi,zis the component of the guide pulley center vector on the z axis of the global coordinate system;Baithe vector of the movable hinge point in the global coordinate system is obtained;Bai,zfor the vector of the pivot point in the z-axis of the global coordinate systemA component; rho is the radius of the guide pulley;
(3) coordinates of rope outlet pointBSi
Figure FDA0003526729300000034
In the formula,Bmi,xis the component of the guide pulley center vector on the x axis of the global coordinate system;Bmi,ythe component of the guide pulley center vector on the y axis of the global coordinate system is shown;
(4) rope length l of working space partW,i
lW,i=qP+lAS=ρ·αi+||Bai-BSi||2
In the formula, qPThe length of the rope wrapped on the guide pulley; lASThe length of the rope between the rope outlet point of the guide pulley and the movable hinge point is determined;
(5) the relationship between the motor rotation angle and the rope length variation is as follows:
Figure FDA0003526729300000041
in the formula,. DELTA.lWIs the effective variation of the rope; h is the winding pitch of the rope;
Figure FDA0003526729300000042
is the motor corner; v. ofPGThe transmission ratio between the motor and the winch is set; r is the hoisting radius.
3. The force sense simulation control method of the fully constrained space rope driven parallel mechanism according to claim 1, wherein in the fourth step, the dynamic modeling comprises a moving platform, a hoisting and guiding pulley and a rope, wherein:
(1) adopting a Newton Euler method to model a mobile platform:
Figure FDA0003526729300000043
in the formula (f)2iIs the platform side rope tension;
(2) a moment balance method is adopted to model a hoisting and guide pulley:
Figure FDA0003526729300000044
in the formula, JrEquivalent moment of inertia of the winch model and the pulley model;
Figure FDA0003526729300000045
the rotation angle vector of the ith motor; crAn equivalent damping coefficient matrix of a winch model and a pulley model; lambda [ alpha ]sThe equivalent radius of the winch model and the pulley model; f. of1iThe rope tension of the ith rope driving side; t isact,iIs the torque of the drive;
(3) modeling the rope as a spring damping model;
Figure FDA0003526729300000051
4. the force sense analog control method of the fully-constrained space rope-driven parallel mechanism according to claim 1, characterized in that in the fifth step, in a motion control frame based on joint feedback, the expected rope length variation of each rope is obtained by inverse kinematics solution from the expected pose, the result is converted into an expected motor rotation angle, and then the expected motor rotation angle is compared with the fed-back actual rotation angle to obtain joint deviation, and the joint deviation is input into a PID controller to realize position control; in a motion control frame based on pose feedback, comparing an expected pose with a model measurement pose to obtain a pose error, obtaining a rope length error according to a Jacobian matrix which is different from the rope length error and the pose error, realizing the conversion from a task space error to a joint space error, further solving a corner error of a motor and inputting the corner error into a PD controller to realize position control; in the internal force coordination loop, expected tension is solved according to a tension distribution method, and the expected tension and rope tension obtained by model measurement form force internal loop feedback.
5. The method for simulating and controlling the force sense of a rope-driven parallel mechanism in the fully constrained space according to claim 1, wherein in the sixth step, the transfer function between the force and the motion is expressed as:
Figure FDA0003526729300000052
in the formula, MdIs equivalent mass; ddEquivalent damping is adopted; kdIs the equivalent stiffness; and E is the displacement correction quantity of the movable platform.
6. The force sense simulation control method of the fully-constrained space rope-driven parallel mechanism according to claim 1, wherein in the seventh step, if the parallel mechanism needs to realize force sense simulation, the overall control frame further comprises a force sense simulation ring, and the eighth step is executed;
and step eight, simulating different force senses and verifying the feasibility of human-computer interaction.
CN202110586409.6A 2021-05-27 2021-05-27 Force sense simulation control method of fully-constrained space rope driving parallel mechanism Active CN113305809B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110586409.6A CN113305809B (en) 2021-05-27 2021-05-27 Force sense simulation control method of fully-constrained space rope driving parallel mechanism

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110586409.6A CN113305809B (en) 2021-05-27 2021-05-27 Force sense simulation control method of fully-constrained space rope driving parallel mechanism

Publications (2)

Publication Number Publication Date
CN113305809A CN113305809A (en) 2021-08-27
CN113305809B true CN113305809B (en) 2022-04-05

Family

ID=77375679

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110586409.6A Active CN113305809B (en) 2021-05-27 2021-05-27 Force sense simulation control method of fully-constrained space rope driving parallel mechanism

Country Status (1)

Country Link
CN (1) CN113305809B (en)

Families Citing this family (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114347005B (en) * 2022-03-18 2022-07-15 中国科学技术大学 A continuous reconfiguration planning method for a rope-pulled parallel robot capable of avoiding obstacles
CN114367964B (en) * 2022-03-22 2022-07-15 中国科学技术大学 A method for reconfiguration planning of rope-pulled parallel robots
CN114643584B (en) * 2022-05-17 2022-09-30 中国科学技术大学 A fast terminal sliding mode synchronization control method for a rope-pulled parallel robot
CN114706312B (en) * 2022-06-06 2022-09-30 中国科学技术大学 A high-precision adaptive cooperative control method for a rope-pulled parallel robot
CN114933031B (en) * 2022-06-17 2025-04-01 北京航空航天大学 A multi-degree-of-freedom platform
CN115096496B (en) * 2022-06-23 2023-03-17 大连理工大学 A method for decoupling measurement of spatial six-dimensional force using cable drive mechanism
CN115096594B (en) * 2022-06-23 2023-09-19 大连理工大学 A device that uses a cable drive mechanism to load and measure six-dimensional forces in space
CN115862407A (en) * 2022-12-05 2023-03-28 中仿智能科技(上海)股份有限公司 Rope traction simulator
CN119347787B (en) * 2024-12-23 2025-03-25 启元实验室 Control method, system and storage medium for adaptive anti-relaxation of rope-driven joints

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107182271B (en) * 2011-12-16 2014-09-03 总装备部工程设计研究总院 Rope drive system pulling force and displacement self-adaptation control method in parallel
CN105974797A (en) * 2016-07-01 2016-09-28 厦门大学 Rope-driven parallel robot motion control method in consideration of elasticity effect and compensation
CN107009348A (en) * 2017-04-18 2017-08-04 中国科学技术大学 A kind of multi-configuration rope driving parallel robot and its spatial pose method for solving
CN107414811A (en) * 2017-06-30 2017-12-01 深圳市波心幻海科技有限公司 Two phase flow parallel connection impedance driver and robot and control method
CN108748091A (en) * 2018-05-21 2018-11-06 哈尔滨工业大学 Rope driving gravity compensation system in parallel towards big reachable tree and rope driving unit via Self-reconfiguration Method
CN109724482A (en) * 2019-01-07 2019-05-07 哈尔滨工业大学 A simulating equipment for recoverable rocket landing conditions based on a rope-driven parallel robot
CN110450135A (en) * 2019-07-30 2019-11-15 华中科技大学鄂州工业技术研究院 A kind of dynamic suspension type weight support system
CN111400917A (en) * 2020-03-18 2020-07-10 中国科学技术大学 A kinematic optimization solution method for a rope-pulled parallel robot with variable structure

Family Cites Families (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102088933B (en) * 2008-05-20 2015-03-18 伯克利仿生技术公司 Device and method for decreasing energy consumption of a person by use of a lower extremity exoskeleton
CN104647367B (en) * 2014-12-29 2016-05-25 合肥工业大学 The robot palletizer in parallel of the compound driving of a kind of rope bar
US10639785B2 (en) * 2015-05-18 2020-05-05 The Regents Of The University Of California Variable force generators for arm supporting exoskeletons

Patent Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107182271B (en) * 2011-12-16 2014-09-03 总装备部工程设计研究总院 Rope drive system pulling force and displacement self-adaptation control method in parallel
CN105974797A (en) * 2016-07-01 2016-09-28 厦门大学 Rope-driven parallel robot motion control method in consideration of elasticity effect and compensation
CN107009348A (en) * 2017-04-18 2017-08-04 中国科学技术大学 A kind of multi-configuration rope driving parallel robot and its spatial pose method for solving
CN107414811A (en) * 2017-06-30 2017-12-01 深圳市波心幻海科技有限公司 Two phase flow parallel connection impedance driver and robot and control method
CN108748091A (en) * 2018-05-21 2018-11-06 哈尔滨工业大学 Rope driving gravity compensation system in parallel towards big reachable tree and rope driving unit via Self-reconfiguration Method
CN109724482A (en) * 2019-01-07 2019-05-07 哈尔滨工业大学 A simulating equipment for recoverable rocket landing conditions based on a rope-driven parallel robot
CN110450135A (en) * 2019-07-30 2019-11-15 华中科技大学鄂州工业技术研究院 A kind of dynamic suspension type weight support system
CN111400917A (en) * 2020-03-18 2020-07-10 中国科学技术大学 A kinematic optimization solution method for a rope-pulled parallel robot with variable structure

Non-Patent Citations (6)

* Cited by examiner, † Cited by third party
Title
A Portable Intuitive Haptic Device on a Desk for User-Friendly Teleoperation of a Cable-Driven Parallel Robot;Park, JH;《APPLIED SCIENCES-BASEL》;20210523;全文 *
Kinematics Modeling of Rigid and Flexible Parallel Joint Mechanisms Driven by Rope;YiTao Pan;《2018 IEEE International Conference on Information and Automation (ICIA)》;20190826;全文 *
Modal space decoupled controller for hydraulically driven six degree of freedom parallel robot;Peter O. Ogbobe;《2010 2nd International Conference on Mechanical and Electronics Engineering》;20100826;全文 *
基于Lagrange方程的3RPS/UPS冗余驱动并联机构柔顺控制;崔学良;《计算机集成制造系统》;20161031;全文 *
绳系并联支撑机构的绳迟滞效应及影响实验研究;冀洋锋;《工程力学》;20191130;全文 *
考虑被动关节阻尼的液压驱动Stewart平台模态空间控制策略;田体先;《华南理工大学学报(自然科学版)》;20150630;全文 *

Also Published As

Publication number Publication date
CN113305809A (en) 2021-08-27

Similar Documents

Publication Publication Date Title
CN113305809B (en) Force sense simulation control method of fully-constrained space rope driving parallel mechanism
CN104477420B (en) Nine rope driven machine people devices of a kind of zero-gravity simulation and low gravity
Yeo et al. Design and analysis of cable-driven manipulators with variable stiffness
Yoshikawa Multifingered robot hands: Control for grasping and manipulation
Kurazume et al. A new index of serial-link manipulator performance combining dynamic manipulability and manipulating force ellipsoids
Picard et al. Pulleys and force sensors influence on payload estimation of cable-driven parallel robots
Yoo et al. Position-tracking control of dual-rope winch robot with rope slip compensation
Picard et al. Control strategies for a cable-driven parallel robot with varying payload information
Seo et al. Vibration reduction of flexible rope-driven mobile robot for safe facade operation
Picard et al. Control solution for a cable-driven parallel robot with highly variable payload
Filipovic et al. The significance of adopted Lagrange’s principle of virtual work used for modeling aerial robots
Phuoc Tho et al. Evaluating cable tension distributions of CDPR for virtual reality motion simulator
Lee et al. Façade operation robot on convex surface using an embedded rope ascender
You et al. Dynamic control of a 3-DOF cable-driven robot based on backstepping technique
Xu et al. A novel cable-driven joint module for space manipulators: Design, modeling, and characterization
CN118579173A (en) Bionic biped robot balance control method and humanoid robot system
Picard Modeling and robust control of cable-driven parallel robots for industrial applications
Collard et al. Dynamics modelling of large suspended parallel cable-driven robots
Torres-Mendez et al. Analysis of a high stiffness warehousing cable-based robot
CN116968057A (en) Parallel rope drives balanced control system
Zhu et al. Optimization design and dynamic stability analysis of 3-DOF cable-driven parallel robot with an elastic telescopic rod
CN112276945B (en) External active gravity compensation system of robot and simulation verification method
JPH0938877A (en) Direct teaching control device
Rashvand et al. Kinematic and dynamic analysis of CPR-A: An aerial cable suspended parallel robot
Zhao et al. Dynamic analysis and trajectory solution of multi-robot coordinated towing system

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