CN109626044B - 一种单张纸印刷机递纸吸嘴运动轨迹优化设计方法 - Google Patents

一种单张纸印刷机递纸吸嘴运动轨迹优化设计方法 Download PDF

Info

Publication number
CN109626044B
CN109626044B CN201811364889.6A CN201811364889A CN109626044B CN 109626044 B CN109626044 B CN 109626044B CN 201811364889 A CN201811364889 A CN 201811364889A CN 109626044 B CN109626044 B CN 109626044B
Authority
CN
China
Prior art keywords
section
motion
track
suction nozzle
length
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
CN201811364889.6A
Other languages
English (en)
Other versions
CN109626044A (zh
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.)
Fuzhou University
Original Assignee
Fuzhou University
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 Fuzhou University filed Critical Fuzhou University
Priority to CN201811364889.6A priority Critical patent/CN109626044B/zh
Publication of CN109626044A publication Critical patent/CN109626044A/zh
Application granted granted Critical
Publication of CN109626044B publication Critical patent/CN109626044B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • BPERFORMING OPERATIONS; TRANSPORTING
    • B65CONVEYING; PACKING; STORING; HANDLING THIN OR FILAMENTARY MATERIAL
    • B65HHANDLING THIN OR FILAMENTARY MATERIAL, e.g. SHEETS, WEBS, CABLES
    • B65H5/00Feeding articles separated from piles; Feeding articles to machines
    • B65H5/22Feeding articles separated from piles; Feeding articles to machines by air-blast or suction device
    • B65H5/222Feeding articles separated from piles; Feeding articles to machines by air-blast or suction device by suction devices

Landscapes

  • Engineering & Computer Science (AREA)
  • Mechanical Engineering (AREA)

Abstract

本发明涉及一种单张纸印刷机递纸吸嘴运动轨迹优化设计方法,其特征在于轨迹曲线的分段之间采用超椭圆曲线过渡,引入轨迹分段的水平比例系数和垂直比例系数,根据轨迹区域、周期,确定轨迹方程。根据总周期和轨迹分段的长度,确定分段的时间,建立分段运动行程时间多项式。通过轨迹方程和运动行程时间多项式,计算轨迹的位移、速度、加速度。将水平比例系数、垂直比例系数和超椭圆曲线阶次作为设计变量,轨迹整段最大合加速度作为优化设计目标进行设计。该轨迹设计方法设计参数少,轨迹过渡段曲线光滑,吸嘴运动平稳、冲击振动小。

Description

一种单张纸印刷机递纸吸嘴运动轨迹优化设计方法
技术领域
本发明涉及印刷机结构设计技术领域,具体涉及一种单张纸印刷机递纸吸嘴运动轨迹优化设计方法。
背景技术
传统递纸吸嘴机构吸嘴运动轨迹分为四段的设计,虽然能够实现递纸功能,但在水平位移转折点处产生较大的速度波动,在转角过渡段采用圆弧曲线过渡,容易产生较大的加速度,对于机构容易产生一定的振动,降低了机器的使用寿命,增加印刷成本。
发明内容
有鉴于此,本发明的目的在于提供一种单张纸印刷机递纸吸嘴运动轨迹优化设计方法,减小了传统设计带来的机构冲击,降低了机构的振动,得到了形状过渡光滑,平稳高效的运动轨迹。
为实现上述目的,本发明采用如下技术方案:
一种递纸机构吸嘴运动轨迹优化设计方法,包括如下步骤:
步骤S1:根据工作区域长度、周期,确定曲线过度段的轨迹方程;
步骤S2:对过渡段进行曲线积分计算过渡段弧长,求得匀速直线段的长度,并建立运动行程时间多项式;
步骤S3:根据得到的轨迹方程和运动行程时间多项式,计算得到运动的轨迹形状和位移、速度、加速度随时间变化的运动规律图;
步骤S4:将垂直比例系数,水平比例系数和超椭圆曲线阶次作为设计变量,整段最大合加速度作为设计目标进行优化设计,并根据运动的轨迹形状和位移、速度、加速度随时间变化的运动规律图,得到最优运动轨迹和位移、速度、加速度随时间变化的运动规律图。
进一步的,所述步骤S1具体为:
步骤S11:将运动轨迹分为八段;
步骤S12:根据递纸吸嘴工艺要求,确定轨迹总行程水平距离L1、回程上升段垂直距离L2、取纸上升段垂直距离L3、回程下降段后半段水平距离L4、轨迹周期T;
步骤S13:引入垂直比例系数λ1,水平比例系数λ2,超椭圆阶次m,确定曲线过度段的轨迹方程。
进一步的,所述运动分为八段包括:取纸上升段ABC,AB为前半段,BC为后半段;
递纸段CD;回程上升段DEF,前半段DE,后半段EF;回程段FG;回程下降段GHA,前半段GH,后半段HA;其中递纸段和回程段采用直线段。取纸上升段、回程上升段、回程下降段采用超椭圆曲线过渡。
进一步的,所述步骤S2具体为:
步骤S21:对过渡段进行曲线积分计算过渡段弧长,并求得匀速直线段的长度;
步骤S22:根据总周期和每段的长度,确定每段的时间,考虑每段首末位置位移、速度、加速度、跃度的约束,建立行程时间多项式.
进一步的,所述步骤S3具体为:
步骤S31:将周期分为N个微段,确定每段微分时长,根据运动行程时间多项式,对应时长将总行程相应分为N微段,确定每段微分长度;
步骤S32:通过轨迹方程x,y的关系,计算局部坐标下各微段的x,y方向位移、速度、加速度。并通过坐标平移变换;
步骤S33:将局部坐标转化为整体坐标,得到运动轨迹形状,和位移、速度、加速度随时间变化的运动规律图。
进一步的,所述步骤S4具体为:
步骤S41:将垂直比例系数λ1,水平比例系数λ2和超椭圆曲线阶次m作为设计变量,整段最大合加速度作为设计目标;
步骤S42:构建m-λ12-amax图;
步骤S43:根据m-λ12-amax图,计算最小加速度对应的长度比例系数λ1和λ2,超椭圆阶次m,得到最优的运动轨迹形状,和位移、速度、加速度随时间变化的运动规律图。
本发明与现有技术相比具有以下有益效果:
1、本发明通过四个区域长度和一个周期五个参数,对递纸吸嘴运动轨迹进行描述,解决了传统设计参数过多,轨迹描述不佳,设计繁琐的难题;
2、本发明对于递纸吸嘴机构吸嘴运动轨迹进行重新分段,将轨迹几何形状尺寸和运动规律结合进行同步设计,降低了传统设计中二者不协调,导致的运动学性能不佳;
3、本发明在六个拐角处采用超椭圆曲线进行过渡处理,并对超椭圆曲线参数进行优化,降低了传统设计中圆弧过渡产生的过高加速度,有效的降低了机构的振动,增加机器的使用寿命,减小印刷成本。
附图说明
图1是本发明实施例中运动轨迹设计尺寸及过渡段参数示意图;
图2是本发明实施例中微分处理示意图;
图3是本发明实施例中m-λ12-amax图;
图4是本发明实施例中最优行程时间图;
图5是本发明实施例中最优运动轨迹图;
图6是本发明实施例中最优轨迹对应递纸吸嘴运动规律图。
具体实施方式
下面结合附图及实施例对本发明做进一步说明。
请参照图1,本发明提供一种单张纸印刷机递纸吸嘴运动轨迹优化设计方法,按照如下步骤实现:
步骤S1:根据区域长度、周期。将工作区域形状抽象,将运动轨迹分为八段,用超椭圆曲线过渡,引入垂直比例系数λ1,水平比例系数λ2,超椭圆阶次m,确定曲线过度段的轨迹方程。
本实施例,取L1=70mm,L2=11mm,L3=4mm,L4=7mm,T=3s;
超椭圆过渡段阶次m的取值依次为2、4、6;
垂直比例系数取值区间λ1的取值区间为[0.2,0.7];
水平比例系数取值区间λ2的取值区间为[0.01,0.06]。
根据初始设计空间区域大小,使用垂直比例系数λ1,水平比例系数λ2,来确定过渡段的两半轴的大小a、b,利用超椭圆曲线过渡段和直线段将运动轨迹分为八段,对过渡段进行曲线积分计算过渡段弧长,并求得匀速直线段的长度大小,为多项式的建立做准备;
如图1,AB段为第1段,BC段为第2段,CD段为第3段,DE段为第4段,EF段为第5段,FG为第6段,GH为第7段,HA为第8段。除第3、6段为匀速直线段,其余各段为一定运动规律的超椭圆弧段。第1、2段具有相同的横半轴长度,第4、5、6具有相同的横半轴长度,第4、5具有相同纵半轴长度。第8段的横半轴长度为L4的大小。
过渡段参数a、b、m计算方法:
第1段:
Figure GDA0001980608180000051
b1=L3×λ1;m1=m;
第2段:a2=a1;b2=L3×(1-λ1);m2=m;
第4段:
Figure GDA0001980608180000052
m4=m;
第5段:a5=a4
Figure GDA0001980608180000053
m5=m;
第7段:a7=a4;b7=(L2+L3)×(1-λ1);m7=m;
第8段:a8=L4;b8=(L2+L3)×λ1;m8=m;
将各段参数带入
Figure GDA0001980608180000054
即确定曲线过度段的轨迹方程。
步骤S2:对过渡段进行曲线积分计算过渡段弧长,并求得匀速直线段的长度,根据总周期和每段的长度,确定每段的时间,考虑每段首末位置位移、速度、加速度、跃度的约束,建立行程时间多项式;
1)过渡段曲线弧长计算:
超椭圆方程:
Figure GDA0001980608180000061
其中,a为横半轴长度,b为纵半轴长度,m为超椭圆阶次
以局部坐标系第一象限为例:
Figure GDA0001980608180000066
(若为二三四象限曲线,则取相应的正切值)
Figure GDA0001980608180000062
对x(θ)、y(θ)求一阶导数,该段弧长为:
Figure GDA0001980608180000063
2)匀速直线段长度计算:
l3=L1-L4-a4
l6=L1-2×a4
3)计算各段时间:
将各段长度作为权重,根据总周期T,计算出各段时长,以及直线段的速度大小;
Figure GDA0001980608180000064
Figure GDA0001980608180000065
Ti-第i段时长、T-轨迹总周期、li-各段轨迹长度、L-轨迹总长度。
第3、6段的速度大小:
Figure GDA0001980608180000071
Figure GDA0001980608180000072
其中,δ1,δ2为直线段速度系数,可根据实际要求调整往返速度大小,本例中取值为1;
4)运动行程时间多项式:
采用七次多项式的方法,考虑每段运动的首末路程,速度,加速度,跃度的约束,建立每段的运动规律s(t);
s(t)=a0+a1×t+a2×t2+a3×t3+a4×t4+a5×t5+a6×t6+a7×t7
边界条件:
第1段:
s(0)=0;s'(0)=0;s”(0)=0;s”'(0)=0;
s(T1)=l1;s'(T1)=0;s”(T1)=0;s”'(T1)=0;
第2段:
s(0)=0;s'(0)=0;s”(0)=0;s”'(0)=0;
s(T2)=l2;s'(T2)=v3;s”(T2)=0;s”'(T2)=0;
第4段:
s(0)=0;s'(0)=v3;s”(0)=0;s”'(0)=0;
s(T4)=l4;s'(T4)=0;s”(T4)=0;s”'(T4)=0;
第5段:
s(0)=0;s'(0)=0;s”(0)=0;s”'(0)=0;
s(T5)=l5;s'(T5)=v6;s”(T5)=0;s”'(T5)=0;
第7段:
s(0)=0;s'(0)=v6;s”(0)=0;s”'(0)=0;
s(T7)=l7;s'(T7)=0;s”(T7)=0;s”'(T7)=0;
第8段:
s(0)=0;s'(0)=0;s”(0)=0;s”'(0)=0;
s(T8)=l8;s'(T8)=0;s”(T8)=0;s”'(T8)=0;
求得过渡段的行程时间多项式s1(t)、s2(t)、s4(t)、s5(t)、s7(t)、s8(t)。直线段采用匀速直线运动的方式,由其直线长度除以该段时间计算速度,并得到其行程时间多项式s3(t)、s6(t);
步骤S3:将整个周期分为N个微段,确定每段微分时长,根据运动规律多项式,对应时长将总行程相应分为N微段,确定每段微分长度。通过轨迹方程x,y的关系。计算局部坐标下各微段的x,y方向位移、速度、加速度。并通过坐标平移变换,将局部坐标转化为整体坐标,得到运动的轨迹形状,和位移、速度、加速度随时间变化的运动规律图;
如图2,以局部坐标系第一象限为例:
微段时长:
Figure GDA0001980608180000081
微分弧长为:u(i+1)=s(ti+1)-s(ti);
轨迹方程:
Figure GDA0001980608180000082
求导:
Figure GDA0001980608180000083
y方向增量:
Figure GDA0001980608180000084
y方向位移:yi+1=yi+Δyi+1
x方向位移:
Figure GDA0001980608180000085
x方向增量:Δxi+1=xi+1-xi
当|x'(yi)|>|y'(xi)|时
x方向增量:
Figure GDA0001980608180000091
x方向位移:xi+1=xi+Δxi+1
y方向位移:
Figure GDA0001980608180000092
y方向增量:Δyi+1=yi+1-yi
采用以下方式计算各节点x、y方向速度。
Figure GDA0001980608180000093
采用以下方式计算各节点x、y方向加速度。
Figure GDA0001980608180000094
其它过渡段通过循环迭代同理可得局部坐标系下的轨迹形状尺寸。坐标平移变换可得到整体坐标系下x,y的坐标,即轨迹形状,如图5。将位移、速度、加速度与时间对应可得到运动轨迹的随时间变化的运动规律图。
步骤S4:将垂直比例系数λ1,水平比例系数λ2和超椭圆曲线阶次m作为设计变量,整段最大合加速度作为设计目标进行优化设计,在m-λ12-amax图中找到最小合加速度对应的垂直比例系数λ1、水平比例系数λ2和超椭圆曲线阶次m值的大小,得到最优运动轨迹及运动规律图。
1)寻找最优参数:
在每段建立加速度数组,将每次每段循环所得的加速度存入数组,得到m、λ1、λ2、amax对应的数据列表。从数据列表中导出数据,如图3建立λ12-m-amax图,找到最小加速度对应的长度比例系数λ1和λ2,超椭圆阶次m;
如图3a),当m=2时,运动轨迹合加速度显示有最小区域。缩小搜索区域,取m=2,将λ1、λ2作为变量,计算各次参数对应合加速度大小,如图3b),即当m=2,λ1=0.6,λ2=0.02时有最小合加速度,9.05m/s2,如图6d)。即最优设计结果。与其他比例系数下的合加速度进行对比,为最优设计结果。
2)得到设计结果:
输出m=2,λ1=0.6,λ2=0.02对应下的吸嘴运动行程时间多项式,如图4;最优运动轨迹,图5;最优轨迹对应递纸吸嘴位移-速度-加速度曲线,如图6。为机构的反向设计提供重要参考。
以上所述仅为本发明的较佳实施例,凡依本发明申请专利范围所做的均等变化与修饰,皆应属本发明的涵盖范围。

Claims (6)

1.一种递纸吸嘴机构吸嘴运动轨迹优化设计方法,其特征在于,包括如下步骤:
步骤S1:根据工作区域长度、周期,确定曲线过渡段的轨迹方程;
步骤S2:对过渡段进行曲线积分计算过渡段弧长,求得匀速直线段的长度,并建立运动行程时间多项式;
步骤S3:根据得到的轨迹方程和运动行程时间多项式,计算得到运动的轨迹形状和位移、速度、加速度随时间变化的运动规律图;
步骤S4:将垂直比例系数,水平比例系数和超椭圆曲线阶次作为设计变量,整段最大合加速度作为设计目标进行优化设计,并根据运动的轨迹形状和位移、速度、加速度随时间变化的运动规律图,得到最优运动轨迹和位移、速度、加速度随时间变化的运动规律图。
2.根据权利要求1所述的一种递纸吸嘴机构吸嘴运动轨迹优化设计方法,其特征在于:所述步骤S1具体为:
步骤S11:将运动轨迹分为八段;
步骤S12:根据递纸吸嘴工艺要求,确定轨迹总行程水平距离L1、回程上升段垂直距离L2、取纸上升段垂直距离L3、回程下降段后半段水平距离L4、轨迹周期T;
步骤S13:引入垂直比例系数
Figure 761683DEST_PATH_IMAGE002
,水平比例系数
Figure 898266DEST_PATH_IMAGE004
,超椭圆阶次m,确定曲线过渡段的轨迹方程。
3.根据权利要求2所述的一种递纸吸嘴机构吸嘴运动轨迹优化设计方法,其特征在于:所述运动轨迹分为八段包括:取纸上升段ABC, AB为前半段,BC为后半段;
递纸段CD;回程上升段DEF,前半段DE,后半段EF;回程段FG;回程下降段GHA,前半段GH,后半段HA;其中递纸段和回程段采用直线段;
取纸上升段、回程上升段、回程下降段采用超椭圆曲线过渡。
4.根据权利要求1所述的一种递纸吸嘴机构吸嘴运动轨迹优化设计方法,其特征在于:所述步骤S2具体为:
步骤S21:对过渡段进行曲线积分计算过渡段弧长,并求得匀速直线段的长度;
步骤S22:根据总周期和每段的长度,确定每段的时间,考虑每段首末位置位移、速度、加速度、跃度的约束,建立运动行程时间多项式。
5.根据权利要求1所述一种递纸吸嘴机构吸嘴运动轨迹优化设计方法,其特征在于:所述步骤S3具体为:
步骤S31:将周期分为N个微段,确定每段微分时长,根据运动行程时间多项式,对应时长将总行程相应分为N微段,确定每段微分长度;
步骤S32:通过轨迹方程x,y的关系,计算局部坐标下各微段的x,y方向位移、速度、加速度;并通过坐标平移变换;
步骤S33:将局部坐标转化为整体坐标,得到运动轨迹形状,和位移、速度、加速度随时间变化的运动规律图。
6.根据权利要求2所述的一种递纸吸嘴机构吸嘴运动轨迹优化设计方法,其特征在于:所述步骤S4具体为:
步骤S41:将垂直比例系数
Figure DEST_PATH_IMAGE007
,水平比例系数
Figure 934410DEST_PATH_IMAGE008
和超椭圆曲线阶次m作为设计变量,整段最大合加速度作为设计目标;
步骤S42:构建m-
Figure 626423DEST_PATH_IMAGE007
-
Figure 301118DEST_PATH_IMAGE008
-amax图;
步骤S43:根据m-
Figure 824503DEST_PATH_IMAGE007
-
Figure 192030DEST_PATH_IMAGE008
-amax图,计算最小加速度对应的长度比例系数
Figure 869874DEST_PATH_IMAGE007
Figure 82681DEST_PATH_IMAGE008
,超椭圆阶次m,得到最优的运动轨迹形状,和位移、速度、加速度随时间变化的运动规律图。
CN201811364889.6A 2018-11-16 2018-11-16 一种单张纸印刷机递纸吸嘴运动轨迹优化设计方法 Active CN109626044B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201811364889.6A CN109626044B (zh) 2018-11-16 2018-11-16 一种单张纸印刷机递纸吸嘴运动轨迹优化设计方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201811364889.6A CN109626044B (zh) 2018-11-16 2018-11-16 一种单张纸印刷机递纸吸嘴运动轨迹优化设计方法

Publications (2)

Publication Number Publication Date
CN109626044A CN109626044A (zh) 2019-04-16
CN109626044B true CN109626044B (zh) 2020-04-10

Family

ID=66068169

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201811364889.6A Active CN109626044B (zh) 2018-11-16 2018-11-16 一种单张纸印刷机递纸吸嘴运动轨迹优化设计方法

Country Status (1)

Country Link
CN (1) CN109626044B (zh)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114476864B (zh) * 2022-01-13 2024-05-10 北京精密机电控制设备研究所 一种机电伺服柔索驱动装置的控制方法

Family Cites Families (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
DE2861397D1 (en) * 1978-06-26 1982-01-28 Heidelberger Druckmasch Ag Camless drive for an oscillating front gripper
CN107522005B (zh) * 2017-09-22 2024-02-06 天津长荣科技集团股份有限公司 一种下吸式递纸装置
CN207209484U (zh) * 2017-09-27 2018-04-10 河北鑫宏源印刷包装有限责任公司 覆膜机的送纸机构的送纸头

Also Published As

Publication number Publication date
CN109626044A (zh) 2019-04-16

Similar Documents

Publication Publication Date Title
CN107390691B (zh) 一种agv路径跟踪方法
CN110865610B (zh) 一种基于机床振动抑制的刀具轨迹插值和速度规划方法
CN106393106B (zh) 参数自适应密化的机器人nurbs曲线运动插补方法
CN108890644B (zh) 一种多轴同步轨迹规划方法、***及计算机可读存储介质
CN105527927A (zh) 一种基于机床各旋转轴角加速度分治优化的五轴加工刀轴矢量插值方法
CN109626044B (zh) 一种单张纸印刷机递纸吸嘴运动轨迹优化设计方法
CN113741426B (zh) 一种基于局部点云曲线拟合的机器人加工路径规划方法
CN111897290A (zh) 一种轴加加速度光滑的拐角过渡平滑方法
CN106647623A (zh) 一种几何精度及衔接速度最优化的五轴联动平滑插补方法
CN113635301B (zh) 一种六轴机械臂运动速度控制改进方法
CN109657282A (zh) 一种基于拉格朗日动力学的h型运动平台建模方法
CN107436592A (zh) 基于b样条的精密打磨轨迹生成方法
CN110209123B (zh) 一种机床驱动约束下的五轴数控加工刀轴矢量局部优化方法
CN113341719A (zh) 一种融合teb算法和模糊控制的agv轨迹跟踪控制方法
CN111158318A (zh) 一种非对称性四次曲线柔性加减速规划方法
CN111310106A (zh) 一种基于逐步逼近原始轮廓的b样条拟合优化方法
CN106951632A (zh) 一种铲齿凸轮的结构优化设计方法
CN111551178A (zh) 一种基于最短路径的分段轨迹时间规划方法
CN113359721B (zh) 一种结合运动控制的改进a*的agv路径规划方法
CN110068832A (zh) 一种激光导航agv的高精度定位方法
CN110018635A (zh) 一种并联机器人拾放操作轨迹规划方法及***
CN116880209A (zh) 一种机器人时间-能量最优平滑轨迹规划方法
Bai et al. Quadruped free gait generation for straight-line and circular trajectories
CN115659549A (zh) 一种内曲线液压马达凸轮环曲线修形方法
WO2022153843A1 (ja) 移動装置、および移動装置制御方法

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