高级检索

留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

木材曲线锯送料平台动力学解析及补偿控制策略研究

孟兆新 曹甲甲 朱莉 马婧尧 石晋菘

孟兆新, 曹甲甲, 朱莉, 马婧尧, 石晋菘. 木材曲线锯送料平台动力学解析及补偿控制策略研究[J]. 北京林业大学学报, 2020, 42(2): 159-166. doi: 10.12171/j.1000-1522.20190234
引用本文: 孟兆新, 曹甲甲, 朱莉, 马婧尧, 石晋菘. 木材曲线锯送料平台动力学解析及补偿控制策略研究[J]. 北京林业大学学报, 2020, 42(2): 159-166. doi: 10.12171/j.1000-1522.20190234
Meng Zhaoxin, Cao Jiajia, Zhu Li, Ma Jingyao, Shi Jinsong. Kinetics analysis and strategy of compensation control study for feeding platform of curve saw for wood[J]. Journal of Beijing Forestry University, 2020, 42(2): 159-166. doi: 10.12171/j.1000-1522.20190234
Citation: Meng Zhaoxin, Cao Jiajia, Zhu Li, Ma Jingyao, Shi Jinsong. Kinetics analysis and strategy of compensation control study for feeding platform of curve saw for wood[J]. Journal of Beijing Forestry University, 2020, 42(2): 159-166. doi: 10.12171/j.1000-1522.20190234

木材曲线锯送料平台动力学解析及补偿控制策略研究

doi: 10.12171/j.1000-1522.20190234
基金项目: 中央高校基本科研业务费专项(2572018CP08)
详细信息
    作者简介:

    孟兆新,博士,教授。主要研究方向:机械工程。Email:349325368@qq.com  地址:150040 黑龙江省哈尔滨市香坊区和兴路26号东北林业大学机电工程学院

  • 中图分类号: S777

Kinetics analysis and strategy of compensation control study for feeding platform of curve saw for wood

  • 摘要: 目的送料平台是模仿木材锯切时人工送料研发的新型机构,平台末端姿态的准确性是保障锯切质量的关键。但是平台驱动支链多,误差从硬件方面进行控制很难实现,为提高板材曲线锯切时的送料精确度,本文提出一种补偿控制策略。方法运用拉格朗日法建立平台的动力学传递函数,通过动力学解析,分析平台关键部件的运动特征。通过数学解算、模拟仿真等方法,研究末端姿态的误差及其特性,从而建立相应的误差模型;以传统PID控制为基础,融合RBF神经网络进行参数优化,设计适用于送料平台的单神经元PID控制器,用以补偿各支链的驱动位移误差,实现实时补偿控制;运用Matlab和Adams联合仿真的方法对补偿控制策略进行验证和分析,并将成功的算法移植到试验台控制器中进行板材锯切试验。结果经文中设计的单神经元PI经文中设计的单神经元PID算法补偿后,送料平台末端轨迹曲线在XY方向偏移误差从补偿前的3 mm降至小于等于1.5 mm,倾斜角误差从补偿前的3.5°降至1.5°,平台末端轨迹曲线与指令曲线相比,在大部分曲线段完全重合。结论文中所设计的单神经元PID补偿控制策略可以有效提高木工锯送料平台末端姿态的准确性,经过补偿控制的送料平台可以实现板材的精密曲线锯切任务。
  • 图  1  一种送料平台结构示意图

    D、E、F、M、N为丝杠传动组;LD、LE、LF为连杆;Oj为等效锯切点;l1l2为丝杠长度;l5为连杆LE和夹具铰接点与丝杠E的距离;a为丝杠间距;b为相邻铰接点间距;hY向丝杠M与光杆间距;θ'为夹具结构角;θ31为丝杠D和连杆LD夹角;θ32为丝杠E和连杆LE夹角;θ33为丝杠F和连杆LF夹角。D, E, F, M and N are screw drive modules,LD, LE and LF are connecting-rods,Oj is equivalent sawing point, l1 and l2 are length of screw,l5 is distance from joint that clamp and rod LE to screw E,a is distance between adjacent screw,b is the distance between adjacent joint,h is distance between screw M and polish rod of direction Yθ' is structural angle of clamp,θ31 is included angle between screw D and rod LDθ32 is included angle between screw E and rod LEθ33 is included angle between screw F and rod LF.

    Figure  1.  Diagrammatic sketch of feeding platform

    图  2  支链D结构解析图

    OjXjYj为等效锯切点坐标系;lD4为螺母D初始位置到等效锯切点距离;θ为夹具转角;βde夹角;d为连杆LD和夹具的铰接点与等效锯切点间距离;e为夹具底边中心点与等效锯切点间距离;yD为螺母D初始位置到丝杠D右侧的距离;LD为连杆。OjXjYj is coordinate system of equivalent sawing point,lD4 is distance from nut D initial position to equivalent sawing point,θ is rotation angle of clamp,β is include angle of d and ed is distance from joint of rod LD and clamp to equivalent sawing point,e is distance from center point of clamp’s bottom line to equivalent sawing point,yD is distance from nut D initial position to right of screw D,LD is connection-rod.

    Figure  2.  Analysis of branched chain D

    图  3  LD、LE、LF连杆伸缩−时间曲线图

    Figure  3.  Graphs between time and connecting-rod LD, LE, LF

    图  4  角度调节机构各连杆形变对末端姿态的影响曲线

    Figure  4.  The curve effected to end posture from deformation of each connect-rod of angle adjusting mechanism

    图  5  丝杠传动误差对连杆转角的影响曲线

    Figure  5.  Curves effected to angle of joint rotation from transfer motion error of screw

    图  6  支链D单神经元PID控制器方案

    yi. 输入信号Input signal;yd. 输出信号Output signal;kp. 比例系数Scale factor;ki. 积分系数Integral factor;kd. 微分系数Differential factor;e. 误差Error;u. 补偿后输入信号Input signal after compensation;f. 干扰信号Disturb signal

    Figure  6.  Control strategy of SN-PID of branch D

    图  7  Matlab与Adams联合仿真Simulink模型

    error_x、error_y、error_1、error_2和error_3为各螺母位移误差;xylDlElF为各螺母输入位移;output_x、output_y、output_1、output_2和output_3为各螺母输出位移;output_j为夹具输出转角。error_x, error_y, error_1, error_2 and error_3 are errors of displacement of each nut. x, y, lD, lE and lF are input displacement of each nut. output_x, output_y, output_1, output_2 and output_3 are output displacement of each nut. output_j is rotation angle of clamp.

    Figure  7.  Simulink model for simulation through Matlab and Adams

    图  8  各轴螺母补偿后位移曲线

    Figure  8.  Curves of each screw nut’s displacement after compensation

    图  9  夹具上一点运动轨迹

    Figure  9.  The trajectory of one point on the clamp

    图  10  送料平台试验样机及加工轨迹曲线

    Figure  10.  Model machine of feeding platform and the curve of trajectory

    表  1  夹具转角$\theta $l5关系

    Table  1.   Relationship about rotation angle θ of clamp and l5

    θ/rad0.5060.4540.4010.3670.3490.3310.3140.297
    ${l_5}$/mm 147 138 134 123 123 117 110 104
    θ/rad 0.262 0.244 0.227 0.209 0.175 0.140 0.122
    ${l_5}$/mm 99 95 85 76 65 48 30
    下载: 导出CSV
  • [1] 朱城伟, 朱大昌, 陈健伟, 等. 3-RPS并联机构SMC控制系统仿真[J]. 机械设计与制造, 2018(9):119−122. doi:  10.3969/j.issn.1001-3997.2018.09.033

    Zhu C W, Zhu D C, Chen J W, et al. SMC control system simulation of 3-RPS parallel mechanism[J]. Machinery Design & Manufacture, 2018(9): 119−122. doi:  10.3969/j.issn.1001-3997.2018.09.033
    [2] 柴保明, 王远东. 3TPT并联机构的误差补偿方法[J]. 机床与液压, 2013, 41(17):85−90. doi:  10.3969/j.issn.1001-3881.2013.17.022

    Chai B M, Wang Y D. An error compensating method for 3TPT parallel mechanism[J]. Machine Tool & Hydraulics, 2013, 41(17): 85−90. doi:  10.3969/j.issn.1001-3881.2013.17.022
    [3] Ni Y B, Wu N, Zhong X Y, et al. Dimensional synthesis of a 3-DOF parallel manipulator with full circle rotation[J]. Chinese Journal of Mechanical Engineering, 2015, 28(4): 830−840. doi:  10.3901/CJME.2015.0122.051
    [4] 秦海宁, 姚单, 杜渝. 3-PRR并联机器人空间位置分析及误差研究[J]. 激光杂志, 2017, 38(12):92−95.

    Qing H N, Yao D, Du Y. The analysis and error research of 3-PRR parallel robot space position[J]. Laser Journal, 2017, 38(12): 92−95.
    [5] Vinoth V, Singh Y, Santhakumar M. Indirect disturbance compensation control of a planar parallel (2-PRP and 1-PPR) robotic manipulator[J]. Robotics and Computer-Integrated Manufacturing, 2014, 30(5): 556−564. doi:  10.1016/j.rcim.2014.03.010
    [6] Castillo-Castaneda E, Takeda Y. Improving path accuracy of a crank-type 6-dof parallel mechanism by stiction compensation[J]. Mechanism and Machine Theory, 2008, 43(1): 104−114. doi:  10.1016/j.mechmachtheory.2006.12.002
    [7] 夏筱筠, 林浒. 基于双模糊控制算法的数控机床过象限误差补偿方法研究[J]. 小型微型计算机系统, 2018, 39(5):241−245.

    Xia X J, Lin H. Study on the compensation for the over quadrant error of CNC based on the dual fuzzy control algorithm[J]. Journal of Chinese Computer Systems, 2018, 39(5): 241−245.
    [8] 彭志文, 高宏力, 梁超, 等. RBF神经网络补偿的并联机器人控制研究[J]. 机械设计与制造, 2018(3):252−259. doi:  10.3969/j.issn.1001-3997.2018.03.074

    Peng Z W, Gao H L, Liang C, et al. Research on parallel robot control based on RBF neural network compensation[J]. Machinery Design & Manufacture, 2018(3): 252−259. doi:  10.3969/j.issn.1001-3997.2018.03.074
    [9] Behrad D, Sasan T, Brian S, et al. A novel adaptive neural network compensator as applied to position control of a pneumatic system[J]. Intelligent Control & Automation, 2011, 2(4): 388−395.
    [10] Niu K, Wang J, Ting K L, et al. Output error analysis of planar 2-dof five-bar mechanism[J]. IOP Conference Series: Materials Science and Engineering, 2018, 324(1): 1−5.
    [11] Si H, Wang L. Error compensation in the five-axis flank milling of thin-walled workpieces[J]. Proceedings of the Institution of Mechanical Engineers, 2019, 233(4): 1224−1234. doi:  10.1177/0954405418780163
    [12] Schneider H, Drust M, Puzik A, et al. Compensation of errors in robot machining with a parallel 3D-piezo compensation mechanism[J]. Procedia Cirp, 2013, 7(12): 305−310.
    [13] 丁建文, 郝宁仲, 徐兆军, 等. 数控木工曲线带锯机的研制[J]. 木材加工机械, 2005(3):16−19. doi:  10.3969/j.issn.1001-036X.2005.03.004

    Ding J W, Hao N Z, Xu Z J, et al. Development of the numerical control woodworking curve band saw[J]. Wood Processing Machinery, 2005(3): 16−19. doi:  10.3969/j.issn.1001-036X.2005.03.004
    [14] Niu K J, Wang J, Kwun-Lon T, et al. Output error analysis of planar 2-DOF five-bar mechanism[J]. Applied Mechanics and Materials, 2013, 324(1): 1−5.
  • [1] 王雷, 徐家琛, 甄雅星, 张恒.  基于热重分析的呼和浩特主要园林树种热解动力学研究 . 北京林业大学学报, 2020, 42(2): 87-95. doi: 10.12171/j.1000-1522.20190280
    [2] 蒋利斌, 张会兰, 杨坪坪, 刘文剑.  降雨和坡面流共同作用下的坡面薄层流水动力学特性 . 北京林业大学学报, 2017, 39(8): 77-86. doi: 10.13332/j.1000-1522.20170126
    [3] 孟兆新, 肖定福, 尹航, 李尚, 张聪.  并联式曲线送料平台机构误差补偿方法的研究 . 北京林业大学学报, 2016, 38(9): 95-101. doi: 10.13332/j.1000-1522.20160026
    [4] 张晓滕, 薛磊, 张宇, 储德淼, 母军.  复合NP阻燃剂处理杨木的热解特性与动力学分析 . 北京林业大学学报, 2016, 38(1): 112-117. doi: 10.13332/j.1000--1522.20150173
    [5] 韩婧, 吴益, 赵琳, 艾云苾, 刘燕, 高荣孚.  光周期对促成栽培芍药生长开花和叶绿素荧光动力学影响 . 北京林业大学学报, 2015, 37(9): 62-69. doi: 10.13332/j.1000-1522.20140418
    [6] 黄河山, 李婷, 刘晋浩.  6-HUS并联运动平台动力学建模与仿真分析 . 北京林业大学学报, 2015, 37(4): 143-150. doi: DOI:10.13332/j.1000-1522.20140349
    [7] 张宇, 母军, 李思锦, 赵阳.  阻燃剂硼酸--硼砂对杨木定向刨花板热解特性的影响 . 北京林业大学学报, 2015, 37(1): 127-138. doi: 10.13332/j.cnki.jbfu.2015.01.003
    [8] 钱晶晶, 黎昕, 杨月月, 林健, 迟德富.  匍枝筋骨草悬浮培养中β-蜕皮甾酮的积累与调控 . 北京林业大学学报, 2015, 37(9): 91-100. doi: 10.13332/j.1000-1522.20150159
    [9] 郭庆启, 张娜, 孙瑶, 王振宇.  超声波辅助提取落叶松树皮多酚的动力学和热力学研究 . 北京林业大学学报, 2014, 36(3): 115-120. doi: 10.13332/j.cnki.jbfu.2014.03.018
    [10] 何正斌, 郭月红, 伊松林, 张璧光
    木材超声波-真空协同干燥的动力学研究 . 北京林业大学学报, 2012, 34(2): 133-136.
    [11] 彭金根, 刘燕, 郭翎, 高荣孚.  越冬返青期间小叶黄杨叶色与叶绿素荧光动力学变化 . 北京林业大学学报, 2011, 33(6): 27-34.
    [12] 苏雅静, 孙爱东, 孙红男, 高雪娟.  竹笋壳黄酮提取物对酪氨酸酶活性影响的动力学研究 . 北京林业大学学报, 2010, 32(4): 276-279.
    [13] 傅深渊, 程书娜, 赵广杰, .  不同催化剂条件下毛竹材液化反应动力学研究 . 北京林业大学学报, 2008, 30(6): 119-123.
    [14] 杨新兵, 余新晓, 鲁绍伟, 李永慈.  基于土壤水动力学的树木根系吸水模型构建(Ⅰ)——单株 . 北京林业大学学报, 2008, 30(supp.2): 197-200.
    [15] 杨新兵, 余新晓, 鲁绍伟, 李永慈, .  基于土壤水动力学的树木根系吸水模型构建(Ⅱ)——林分 . 北京林业大学学报, 2008, 30(supp.2): 201-205.
    [16] 王乐乐, 聂立水, 戴伟.  北京山地不同针叶林型下土壤脲酶动力学和热力学特征研究 . 北京林业大学学报, 2008, 30(2): 89-94.
    [17] 周永学, 毛俊娟, 魏潇潇, 胡胜华, 刘杏娥, 高黎, 张洪江, 吴彩燕, 袁怀文, 颜绍馗, 殷亚方, 杨平, 白岗栓, 黄荣凤, 郑小贤, 王费新, 胡万良, 何亚平, 张莉俊, 李瑞, 王芳, 张璧光, 邓小文, 秦爱光, 戴思兰, 王胜华, 王兆印, 樊军锋, 费世民, 崔赛华, 李猛, 王正, 汪思龙, 谭学仁, 王小青, 孙向阳, 乔建平, 张克斌, 罗晓芳, NagaoHirofumi, 赵天忠, 常旭, 杜社妮, 张岩, 刘燕, 王晓欢, KatoHideo, 李昀, 高荣孚, 张双保, 陈放, 龚月桦, 范冰, 江玉林, , 王海燕, 张占雄, 韩士杰, 李华, 徐嘉, 张旭, 刘云芳, 江泽慧, 孔祥文, 刘秀英, 侯喜录, 常亮, 陈宗伟, 陈秀明, 杨培华, , 李媛良, 丁磊, 李晓峰, 郭树花, IdoHirofumi, 任海青, 薛岩, , 高建社, 张代贵, 陈学平, 张桂兰, 李考学, 徐庆祥, 蒋俊明, 费本华, , 金鑫, , 刘永红, 李雪峰, 续九如, 王晓东, 涂代伦, 丁国权, , 张红丽, , .  非线性植被-侵蚀动力学模型初探 . 北京林业大学学报, 2007, 29(6): 123-128.
    [18] 贺窑青, 熊丹, 陈佳, 范丙友, 欧阳杰, 孙青, 胡海英, 周章义, 李艳华, 刘丽, 曲红, 胡晓丹, 段旭良, 李在留, 雷庆哲, 王莉, 刘美芹, 乔海莉, 郝晨, 李莉, 孙月琴, 金莹, 隋金玲, 张玲, 王丰俊, 姚娜, 石娟, 程堂仁, 路端正, 骆有庆, 陈晓阳, 尹伟伦, 孙爱东, 王建中, 张艳霞, 康向阳, 赵亚美, 续九如, 郭锐, 陈发菊, 李凤兰, 张香, 张志毅, 阎伟, 沈昕, 陆海, 尹伟伦, 田呈明, 冯菁, 骆有庆, 张德权, 武彦文, 冯秀兰, 周燕, 郑彩霞, 李云, 卢存福, 胡晓丹, 蒋湘宁, 高述民, 史玲玲, 梁宏伟, 王百田, 吴晓成, 骆有庆, 赵蕾, 武海卫, 郑永唐, 王晓东, 郝俊, 胡德夫, 姜金仲, 安新民, 孙爱东, 李忠秋, 阎晓磊, 王华芳, 骆有庆, 马钦彦, 沈繁宜, 梁华军, 王冬梅, 严晓素, 崔彬彬
    , 张志翔, 李凯, 刘玉军, 冯晓峰, 骈瑞琪, 赵兵, 谢磊, 高荣孚, 王瑛, 王玉兵, 王建中, 吴坚, 温秀凤3, 王晓楠, 王华芳, 邹坤, 郭晓萍, 尹伟伦, 于京民2, 冯仲科, 呼晓姝, 陶凤杰, 杨伟光, 刘艳, 王民中, 林善枝, 丁霞, 孙建华, 李镇宇, 张庆, 陈卫平, 刘玉军, 张兴杰, 王玉春, 沈应柏, 李凤兰, 马建海, 汪植, 赵新丽, 付瑞海, 蒋平.  肉苁蓉细胞悬浮培养动力学研究 . 北京林业大学学报, 2007, 29(5): 132-136.
    [19] 贾彩凤, 何恒斌, 李长洪, 张春雨, 于占源, 李秀芬, 何利娟, 马履一, 钟健, 张力平, 杨永福, 李笑吟, 史军辉, 赵广亮, 林峰, 王献溥, 姜春宁, 于格, 郭小平, 王勇, 王华, 赖巧玲, 吕兆林, 张亚利, 许月卿, 曹金珍, 杨明嘉, 曾德慧, 习宝田, 黄忠良, 尚晓倩, 郭惠红, 鲁春霞, 朱教君, 王骏, 李鸿琦, 贾桂霞, 朱清科, 毕华兴, 王希群, 郝玉光, 杨培岭, 赵博光, 郑彩霞, D.PascalKamdem, 孙长霞, 赵秀海, 邵晓梅, 胥辉, 于顺利, 王继兴, 李悦, 甘敬, 朱金兆, 陈宏伟, 丁琼, 郑景明, 欧阳学军, 周金池, 张志2, 崔小鹏, 丁琼, 李黎, 费孛, 包仁艳, 王秀珍, 刘燕, 任树梅, 谢高地, 姜凤岐, 尚宇, 张榕, 王庆礼, 杨为民, 范志平, 贾桂霞, 沈应柏, 刘足根, 刘鑫, 刘艳, 贾昆锋, 何晓青, 张中南, 沈应柏, , 张池, , 蔡宝军, 毛志宏, 唐小明, 鹿振友, 李林, 李凤兰, 纳磊, 陈伏生, 周金池, 张方秋, 周小勇, 赵琼, 马玲, , 申世杰, .  植物单宁改性树脂吸附机理的研究 . 北京林业大学学报, 2006, 28(4): 6-11.
    [20] 李世东, 时尽书, 孙晓梅, 徐双民, 潘存德, 张宇清, 金小娟, 窦军霞, 张冰玉, 谭伟, 翟明普, 王玉杰, 颜容, 吕建雄, 胡晓丽, 范丙友, 李绍才, 肖生春, 南海龙, 朱教君, 陈文汇, 李发东, 杨振德, 高峻, 肖洪浪, 师瑞峰, 王云琦, 胡诗宇, 张守攻, 朱清科, 刘红霞, 三乃, 张一平, 李建章, 康宏樟, 田小青, 孙海龙, 苏晓华, 谢益民, 周春江, 韩海荣, 冯仲科, 孟平, 刘俊昌, 宋献方, 骆秀琴, 齐实, 赵博光, 李义良, 李智辉, 杨志荣, 周文瑞, 蒋佳荔, 王笑山, 蔡怀, 陆海, 吴斌, 赵双菊, 张雁, 齐实, 马钦彦, 姜伟, 刘昌明, 岳良松, 朱金兆, 姚山, 蒲俊文, 于静洁, 伊力塔, 宋清海, 葛颂, 张岩, 齐力旺, 张劲松, 张永安, 蒋湘宁, 赵有科, 何磊, 张德荣, 吴庆利, 康峰峰, 刘元, 杨聪, 褚建民, 石丽萍, 吕守芳, 马超德, 崔保山, 曲良建, 刘鑫宇, 朱林峰, 王建华, 王玉珠, 刘相超, 田颖川, 胡堃, 唐常源.  小钻杨和刺槐根系养分吸收的动力学研究 . 北京林业大学学报, 2006, 28(2): 29-33.
  • 加载中
图(10) / 表 (1)
计量
  • 文章访问数:  1138
  • HTML全文浏览量:  657
  • PDF下载量:  28
  • 被引次数: 0
出版历程
  • 收稿日期:  2019-05-28
  • 修回日期:  2019-06-24
  • 网络出版日期:  2019-10-19
  • 刊出日期:  2020-03-03

木材曲线锯送料平台动力学解析及补偿控制策略研究

doi: 10.12171/j.1000-1522.20190234
    基金项目:  中央高校基本科研业务费专项(2572018CP08)
    作者简介:

    孟兆新,博士,教授。主要研究方向:机械工程。Email:349325368@qq.com  地址:150040 黑龙江省哈尔滨市香坊区和兴路26号东北林业大学机电工程学院

  • 中图分类号: S777

摘要: 目的送料平台是模仿木材锯切时人工送料研发的新型机构,平台末端姿态的准确性是保障锯切质量的关键。但是平台驱动支链多,误差从硬件方面进行控制很难实现,为提高板材曲线锯切时的送料精确度,本文提出一种补偿控制策略。方法运用拉格朗日法建立平台的动力学传递函数,通过动力学解析,分析平台关键部件的运动特征。通过数学解算、模拟仿真等方法,研究末端姿态的误差及其特性,从而建立相应的误差模型;以传统PID控制为基础,融合RBF神经网络进行参数优化,设计适用于送料平台的单神经元PID控制器,用以补偿各支链的驱动位移误差,实现实时补偿控制;运用Matlab和Adams联合仿真的方法对补偿控制策略进行验证和分析,并将成功的算法移植到试验台控制器中进行板材锯切试验。结果经文中设计的单神经元PI经文中设计的单神经元PID算法补偿后,送料平台末端轨迹曲线在XY方向偏移误差从补偿前的3 mm降至小于等于1.5 mm,倾斜角误差从补偿前的3.5°降至1.5°,平台末端轨迹曲线与指令曲线相比,在大部分曲线段完全重合。结论文中所设计的单神经元PID补偿控制策略可以有效提高木工锯送料平台末端姿态的准确性,经过补偿控制的送料平台可以实现板材的精密曲线锯切任务。

English Abstract

孟兆新, 曹甲甲, 朱莉, 马婧尧, 石晋菘. 木材曲线锯送料平台动力学解析及补偿控制策略研究[J]. 北京林业大学学报, 2020, 42(2): 159-166. doi: 10.12171/j.1000-1522.20190234
引用本文: 孟兆新, 曹甲甲, 朱莉, 马婧尧, 石晋菘. 木材曲线锯送料平台动力学解析及补偿控制策略研究[J]. 北京林业大学学报, 2020, 42(2): 159-166. doi: 10.12171/j.1000-1522.20190234
Meng Zhaoxin, Cao Jiajia, Zhu Li, Ma Jingyao, Shi Jinsong. Kinetics analysis and strategy of compensation control study for feeding platform of curve saw for wood[J]. Journal of Beijing Forestry University, 2020, 42(2): 159-166. doi: 10.12171/j.1000-1522.20190234
Citation: Meng Zhaoxin, Cao Jiajia, Zhu Li, Ma Jingyao, Shi Jinsong. Kinetics analysis and strategy of compensation control study for feeding platform of curve saw for wood[J]. Journal of Beijing Forestry University, 2020, 42(2): 159-166. doi: 10.12171/j.1000-1522.20190234
  • 对机构进行软件层面的补偿控制是提高机械加工准确性快速且有效的方法。国内外学者在误差补偿控制的研究中已取得了一定的成果,如朱城伟等[1]设计出SMC-PID控制器,有效降低机构误差;柴保明等[2]对3TPT并联机构建立误差模型并对模型研究提出一种误差补偿方法;Ni等[3]利用一阶摄动建立了误差矩阵对三自由度并联机械手实现机构补偿控制;秦海宁等[4]对3-PRR并联机构的误差进行了定性分析,建立了误差测量实验平台,进行了误差消除试验,提出补偿控制方法;Vinoth等[5]针对自己设计的三自由度平面并联机构建立了基于EKF的间接自适应补偿控制;Castillo-Castaneda等[6]对一种6-DOF并联机构设计了基于PMAC运动控制器的轨迹误差补偿算法;夏筱筠等[7]运用双模糊控制算法对数控机床过象限误差进行补偿;彭志文等[8]提出了基于计算力矩控制的RBF神经网络在线补偿控制策略;Behrad等[9]研究出一种新型自适应在线神经网络补偿器应用于气动龙门机器人的位置控制中;Niu等[10]研究了平面二自由度五杆运动副背隙误差,并建立了误差模型和控制策略;Si等[11]提出了一种迭代补偿策略,减小五轴侧铣薄壁件时刀具或工件形变引起的误差;Schneider等[12]提出了一种三维压电补偿控制单元以及匹配的控制算法用以补偿机器人加工系统的加工精度。

    为实现木工板材曲线锯切,模仿人工送料研发机械装置代替人工,本文提出一种新型五自由度送料平台。相较于现有木工曲线锯切[13],该平台避免了锯转动、工件插补运动时锯条承受转矩过大的问题。送料机构由三自由度并联机构和十字滑台串联而成,一方面,虽然许多专家对相关类型三自由度并联机构的控制提出了补偿策略,但已有研究成果均具有局限性,并不适用于本文提出的机构;另一方面,送料平台的误差模型较难建立,且不具备可移植性。为满足木材板料的精确加工要求,实现木工带锯机的精确送料,需要分析送料平台误差来源,建立补偿控制策略。平台姿态的准确性直接决定了木工板料锯切的质量,因此分析机构的误差来源,消除传动误差的影响,保证瞬时姿态的准确性,建立补偿控制策略是必要的。

    • 建立机构的数学模型,求解传递函数是动力学解析的关键。为简化建模对机构作如下假设:其一,各个杆件的质量均匀分布,约束为理想约束;其二,忽略杆件弹性变形和对机构产生负面影响的摩擦力;其三,杆件的质量集中于杆件质心。

    • 图1所示,送料平台为五自由度混联机构,其角度调节机构底座下方为十字丝杠传动模组X-YY轴一边为丝杠传动模组,另一边为起辅助支撑作用的光杆;角度调节机构底座上方平行固定3组角度调节丝杠传动模组D、E、F,各自螺母与连杆铰接,三连杆又分别与夹具铰接。为求解各电机补偿控制输入量,建立其逆解数学模型。对机构进行分析,以支链D(丝杠传动模组D所在运动传动链,简称支链D,下同)为例,分析几何模型如图2所示。

      图  1  一种送料平台结构示意图

      Figure 1.  Diagrammatic sketch of feeding platform

      图  2  支链D结构解析图

      Figure 2.  Analysis of branched chain D

      假设锯切轨迹曲线为$y = f(x)$,则某一时刻等效锯切点坐标为$({x_0},f({x_0}))$,各轴螺母的相对位移量lxlylDlElF可运用多次正余弦定理求得,求解如式(1),式(1)中与机构动态相关的变量θβedL5求解如式(2)。

      $$\left\{ \begin{array}{l} {l_x} = {x_0} + {l_5} \\ {l_y} = {y_0} + {l_{\rm{E}}} - \sqrt {l_{\rm{E}}^2 - l_5^2} \\ {l_{\rm{D}}} = {l_{{\rm{D}}4}} - d \cos \left( {\beta + \theta } \right) - \sqrt {{l_{\rm{D}}}^2 - {{\left[ {a - d \sin \left( {\beta + \theta } \right)} \right]}^2}} \\ {l_{\rm{E}}} = {l_{{\rm{E}}4}} - e \cos \theta - \sqrt {{l_{\rm{E}}}^2 - {{\left( {e \sin \theta } \right)}^2}} \\ {l_{\rm{F}}} = {l_{{\rm{F}}4}} - d \cos \left( {\beta - \theta } \right) - \sqrt {{l_{\rm{F}}}^2 - {{\left[ {a - d \sin \left( {\beta - \theta } \right)} \right]}^2}} \\ \end{array} \right.$$ (1)

      式中:lD4lE4lF4分别为螺母D(丝杠传动模组D的螺母,简称螺母D,下同)、螺母E和螺母F初始位置到等效锯切点间距(单位:mm);lDlElF分别为连杆LD、LE和LF杆长(单位:mm);l5为连杆LE和夹具铰接点与丝杠E的距离(单位:mm);a为并联机构相邻丝杠间距离(单位:mm);d为连杆LD和夹具的铰接点与等效锯切点间距离;e为夹具底边中心点与等效锯切点间距离(单位:mm);b为相邻铰接点间距(单位:mm);θ为夹具转角(单位:rad);βde夹角(单位:rad)。

      $$\left\{ \begin{array}{l} \theta =\arctan (f'({x_0})),\;{l_5} = {f_1}(\theta ) \\ e = {f_2}(\theta ) = {f_1}(\theta ) \csc (\theta ) + b \sin (\theta ') \\ d = {f_3}(\theta ) = {(b \cos \theta ')^2} + {(e + b \sin \theta ')^2} \\ \beta = {f_4}(\theta ) = \arccos [(b \cos \theta )/d] \\ \end{array} \right.$$ (2)

      式中:$\theta '$为夹具结构角(单位:rad);f1θ)为拟合曲线,拟合数据见表1b为相邻铰接点间距(单位:mm)。

      表 1  夹具转角$\theta $l5关系

      Table 1.  Relationship about rotation angle θ of clamp and l5

      θ/rad0.5060.4540.4010.3670.3490.3310.3140.297
      ${l_5}$/mm 147 138 134 123 123 117 110 104
      θ/rad 0.262 0.244 0.227 0.209 0.175 0.140 0.122
      ${l_5}$/mm 99 95 85 76 65 48 30
    • 送料平台各杆件只有绕Z轴的旋转,没有平移,为完整有势系统,其拉格朗日方程可简化为式(3)。定义机构中5个螺母的位移为广义坐标,即${q_1} = {l_x}, $${q_2} = {l_y},{q_3} = {l_{\rm{D}}},{q_4} = {l_{\rm{E}}},{q_5} = {l_{\rm{F}}}$

      $$\frac{{\rm{d}}}{{{\rm{d}}t}}\left( {\frac{{\partial T}}{{\partial {{\dot q}_\alpha }}}} \right) - \frac{{\partial T}}{{\partial {q_\alpha }}} = {{{Q}}_\alpha }\;\;(\alpha = 1,2, 3,4,5)$$ (3)

      式中:T为系统动能(单位:J);${ Q_\alpha }$为输入广义驱动力矩矢量和(单位:N·mm);${q_\alpha }$为广义坐标。

      求解各杆件转动惯量如式(4)所示:

      $${J_{\rm{s}}} = \frac{{{m_{\rm{s}}}{{\phi}^2}}}{8},{J_{\rm{l}}} = \frac{{{m_{\rm{l}}}{l^2}}}{3},{J_{\rm{j}}} = \frac{{m({h^2} + {w^2})}}{{12}} + ml_{\rm{d}}^2 $$ (4)

      式中:Js为丝杠转动惯量(单位:kg·m2);Jl为连杆转动惯量(单位:kg·m2);Jj为夹具转动惯量(单位:kg·m2);m为夹具质量,ms为丝杠质量,ml为连杆质量(单位:kg);$\phi $为直径(单位:mm);l为杆长(单位:mm);hw分别为夹具长、宽(单位:mm);${l_{\rm{d}}}$为夹具转动中心到夹具型心距离(单位:mm)。

      支链D、E、F中连杆与丝杠的夹角分别为θ31θ32θ33,其值如式(5):

      $$\left\{ \begin{array}{l} {\theta _{31}}=\arccos [({l_{\rm{D}}} \cos \gamma - {q_3})/{l_{\rm{D}}}] \\ {\theta _{32}}=\arccos [({l_{\rm{E}}} - {q_4})/{l_{\rm{E}}}] \\ {\theta _{33}}=\arccos [({l_{\rm{F}}} \cos \gamma - {q_5})/{l_{\rm{F}}}] \\ \end{array} \right.$$ (5)

      式中:γ为初始状态时支链D中连杆与丝杠夹角(单位:rad)。

      设夹具质地均匀,几何相对质心坐标为$({l_5},{l_{{\rm{D}}4}} - $${l_{\rm{E}}} \cos {\theta _{32}})$,转角$\theta = {\rm{arctan[}}{l_5}{\rm{/(}}{l_{{\rm{D}}4}}{\rm{ - }} $${l_{\rm{E}}} \cos {\theta _{32}}{\rm{)]}}$,所以机构的动能为$T = {T_{\rm{s}}} + {T_{\rm{m}}} + {T_{\rm{l}}} + {T_{\rm{j}}}$,其中${T_{\rm{s}}}$为所有丝杠总动能(单位:J);${T_{\rm{m}}}$为所有螺母总动能(单位:J),${T_{\rm{l}}}$为3个连杆总动能(单位:J),${T_{\rm{j}}}$为夹具和被加工件动能(单位:J)。其值如式(6)所示:

      $$\left\{ \begin{array}{l} {T_{\rm{s}}} = \dfrac{1}{2}\displaystyle \mathop \sum \limits_{\alpha = 1}^5 {J_{{\rm{s}}\alpha }}{\left( {\dfrac{{{q_\alpha }}}{{{P_\alpha }}}} \right)^2},{T_{\rm{m}}} = \dfrac{1}{2}\mathop \sum \limits_{\alpha = 1}^5 {m_\alpha }{\left( {{{\dot q}_\alpha }} \right)^2}\\ {T_{\rm{l}}} = \dfrac{1}{2}\displaystyle \mathop \sum \limits_{\alpha = 3}^5 {J_{{\rm{l}}\alpha }}{\left( {\dfrac{{{{\dot q}_\alpha }}}{{{k_\alpha }}}} \right)^2},{T_{\rm{j}}} = \dfrac{1}{2}{J_{\rm{j}}}{\left( {\frac{{k_4^{1.5} {L_5} {{\dot q}_4}}}{{L_5^2 + k_4^{1.5}}}} \right)^2} \end{array} \right.$$ (6)

      式中:Pα为丝杠的导程(单位:mm);qα为螺母相对位移量(单位:mm);连杆动能系数${k_1} = $$ \sqrt {{l_{\rm{D}}}^2 - {{\left( {{l_{\rm{D}}}\; \cos \gamma - {q_3}} \right)}^2}} $${k_2} = \sqrt {{l_{\rm{E}}}^2 - {{\left( {{l_{\rm{E}}} - {q_4}} \right)}^2}} $${k_3} = $$\sqrt {{l_{\rm{F}}}^2 - {{\left( {{l_{\rm{F}}}\; \cos \gamma - {q_5}} \right)}^2}} $${k_4} = {l_{{\rm{D}}4}} - {l_{\rm{E}}} + {q_4}$

      对机构总动能求导得:

      $$\frac{{\partial T}}{{\partial {q_\alpha }}} = \mathop \sum \limits_{\alpha = 1}^5 {J_{{\rm{s}}\alpha }} \cdot \frac{{{q_\alpha }}}{{{P_\alpha }}}$$ (7)
      $$\frac{{\partial T}}{{\partial {{\dot q}_\alpha }}} = \mathop \sum \limits_{\alpha = 1}^5 {m_\alpha }{\dot q_\alpha } + \sum\limits_{\alpha =3}^5 {{J_{{\rm{l}}\alpha }}\frac{{{{\dot q}_\alpha }}}{{{k_\alpha }}}} + {J_{\rm{j}}}\frac{{k_4^{1.5} {l_5}}}{{l_5^2 + k_4^{1.5}}}{\dot q_4}$$ (8)

      将式(7)、(8)带入式(3)得送料平台转角与各丝杠进给量之间的关系式:

      $$\mathop \sum \limits_{\alpha = 1}^5 {m_\alpha }{\ddot q_\alpha } + \mathop \sum \limits_{\alpha = 1}^3 {J_{{\rm{l}}\alpha }}\frac{{{{\ddot q}_{\alpha + 2}}}}{{{k_\alpha }}} + {J_{\rm{j}}}\frac{{k_4^{1.5} {l_5}}}{{l_5^2 + k_4^{1.5}}}{\ddot q_4} - \mathop \sum \limits_{\alpha = 1}^5 {J_{{\rm{s}}\alpha }} \cdot \frac{{{q_\alpha }}}{{{P_\alpha }}} = {{{Q}}_\alpha }$$ (9)

      于是机构转角与广义坐标的关系式可以表示为式(10):

      $$\theta {\rm{ = arctan}}\frac{{{J_{\rm j}} {{\ddot q}_4} k_4^{1.5}}}{{{\rm{2}} \! \left[ {{{{Q}}_\alpha } \!-\! \left( {\displaystyle\sum \limits_{\alpha = 1}^5 {m_\alpha }{{\ddot q}_\alpha }\! +\! \displaystyle\sum \limits_{\alpha = 1}^3 {J_{{\rm{l}}\alpha }}\dfrac{{{{\ddot q}_{\alpha + 2}}}}{{{k_\alpha }}} \!-\! \displaystyle\sum \limits_{\alpha = 1}^5 {J_{{\rm{s}}\alpha }} \dfrac{{{q_\alpha }}}{{{P_\alpha }}}} \right)} \right]}}$$ (10)
    • 送料平台位移瞬态的准确性,取决于各传动件瞬态位置的准确性,运用ADAMS解算送料平台模型,测量连杆长度随时间变化的关系,并结合虚功原理求解连杆误差模型;测量螺母的移动规律,分析丝杠的传动误差。

    • 图3所示,LD、LE、LF支链连杆(依次为Rod-1、Rod-2、Rod-3)在运动起始阶段杆长变化频繁,出现反复拉压的现象,但经过0.15 s后归于稳定。杆件伸缩量均不超过2.5 mm,属于微小变形。假设伸缩量在[− 2 mm,2 mm]之间变化,测量形变对送料平台末端姿态的变化影响如图4所示,其中Rod-1、Rod-3的误差仅在局部时间段内使平台末端轨迹偏离理想曲线,而Rod-2误差的影响覆盖了末端轨迹全程。但不难发现,杆的伸缩量与末端轨迹误差呈线性关系。

      图  3  LD、LE、LF连杆伸缩−时间曲线图

      Figure 3.  Graphs between time and connecting-rod LD, LE, LF

      图  4  角度调节机构各连杆形变对末端姿态的影响曲线

      Figure 4.  The curve effected to end posture from deformation of each connect-rod of angle adjusting mechanism

      为了消除连杆形变对末端姿态的影响,依据虚功原理对连杆与夹具组成系统进行解析。已知3个连杆与夹具共18个坐标分量,15个约束,定义广义坐标分别为${\theta _{31}}$${\theta _{32}}$${\theta _{33}}$,列受力平衡方程式求解得式(11)。

      $$\begin{split} {{{F}}_{\rm{D}}} = - \frac{{{J_{{\rm{L}}1}}{\rm{ - }}{J_{\rm{j}}}}}{{{L_{\rm{D}}}\cos {{\theta _{31}}} }}、{{{F}}_{\rm{E}}} = - \frac{{{J_{{\rm{L}}2}} + {J_{\rm{j}}}}}{{{L_{\rm{E}}}\cos {{\theta _{32}}} }}、{{{F}}_{\rm{F}}} = - \frac{{{J_{{\rm{L}}3}} - 2{J_{\rm{j}}}}}{{{L_{\rm{F}}}\cos {{\theta _{33}}} }}\\ \end{split}$$ (11)

      式中:${F_{\rm{D}}}$${F_{\rm{E}}}$${F_{\rm{F}}}$分别为连杆LD、LE、LF所受外力(单位:N);${J_{{\rm{L}}1}}$${J_{{\rm{L}}2}}$${J_{{\rm{L}}3}}$分别为连杆LD、LE、LF的转动惯量(单位:kg·m2);${J_{\rm{j}}}$为夹具和板料转动惯量(单位:kg·m2)。

      连杆变形引起的机构末端姿态误差等于连杆伸长量在丝杠方向的投影长度,求其误差得式(12),可看出连杆引起的误差与时间无关。

      $${e_{{\rm{LD}}}} = - \frac{{{J_{{\rm{L}}1}}{\rm{ - }}{J_{\rm{j}}}}}{{{l_{\rm{DS}}}}}、{e_{{\rm{LE}}}} = - \frac{{{J_{{\rm{L}}2}} + {J_{\rm{j}}}}}{{{l_{\rm{ES}}}}}、{e_{{\rm{LF}}}} = - \frac{{{J_{{\rm{L}}3}} - 2{J_{\rm{j}}}}}{{{l_{\rm{FS}}}}}$$ (12)

      式中:${e_{{\rm{LD}}}}$${e_{{\rm{LE}}}}$${e_{{\rm{LF}}}}$分别为连杆LD、LE、LF形变引起的误差(单位:mm);lDSlESlFS分别为连杆LD、LE和LF的形变量(单位:mm)。

    • 假设螺母位移误差绝对值小于0.5 mm,分析该误差对连杆转角影响如图5(其中支链D、E和F的丝杠简称为丝杠D、丝杠E和丝杠F),图5a5b5c依次为丝杠D、丝杠E、丝杠F的螺母位移误差对连杆转角的影响。很明显较小的位移便能产生较大的转角误差,甚至影响连杆转向,而连杆转角直接影响了末端姿态。依据文献[14]对关节转角误差的分析,考虑到螺旋传动受背隙、摩擦等影响,很难建立一个准确的数学模型,因此为每组丝杠安装尺光栅传感器,用以反馈螺母位置信息,用ls表示,伺服电机自带编码盘作为其位置检测装置,角度用γ表示。于是,丝杠螺母的误差可以表示为${e_{{\rm{s}}\alpha }} = {l_{\rm{s}}} - {\gamma _\alpha } {P_\alpha }/2\pi , $$\left( {\alpha = 1,2, 3,4,5} \right)$

      图  5  丝杠传动误差对连杆转角的影响曲线

      Figure 5.  Curves effected to angle of joint rotation from transfer motion error of screw

    • 设计单神经元PID算法对各支链进行控制层面的误差补偿。为验证该方法的可行性,运用MATLAB和ADAMS仿真进行分析和检验。

    • 应用PID算法对送料平台的各误差进行补偿,是一种有效的方法。但在变输入参数的模型中误差较大,于是引入神经网络算法对PID各系数进行实时调节,以提高PID控制精度。算法结构如图6

      图  6  支链D单神经元PID控制器方案

      Figure 6.  Control strategy of SN-PID of branch D

      依据电机的驱动形式,选用增量式PID控制器,设计单神经元PID,补偿后输入量u可表示为式(13)。

      $$ \begin{aligned} {{u}}(k) =& {k_{\rm{p}}}(e(k) - e(k - 1)) + {k_{\rm{i}}}e(k) + {k_{\rm{d}}}(e(k) - \\ &2e(k - 1) + e(k - 2)) \end{aligned}$$ (13)

      式中:${k_{\rm{p}}}$${k_{\rm{i}}}$${k_{\rm{d}}}$分别为比例、积分和微分系数;$e(k)$为误差,$e(k)={y_{\rm{i}}} - {y_{\rm{d}}}$,其中yiyd分别为控制系统的输入和输出量。

      RBF神经网络结构采用三层前向网络。第一层为丝杠位移误差矩阵,即 ${{X}}{\rm{ = [}}{{{e}}_{\alpha} }{{\rm{]}}^T},(\alpha = 1,2, 3,4, 5)$;第二层为隐含层,径向基函数为${{H}} = {[{h_{\alpha}}]^T}$,其中${h_{\alpha}}$为高斯基函数,${h_{\alpha}} = g(\left\| {{{X}} - {{{C}}_{\alpha}}} \right\|/{b_{\alpha}}^2)$${{{C}}_{\alpha}} = [{c_{{\alpha}1}},{c_{{\alpha}2}}, \cdots, {c_{{\alpha}m}}]$表示网络节点α的高斯基函数中心点,${{b}} = [{b_1},{b_2}, \cdots, $${b_m}]$表示某个节点的高斯基函数宽度;第三层为输出层,${{{K}}_\alpha } = [{k_{{\rm{p}}\alpha }},{k_{{\rm{i}}\alpha }},{k_{{\rm{d}}\alpha }}]$,由于PID的3个系数具有不确定性,设不确定项 ${f_\alpha }(x) = {{{W}}_\alpha }^T{{{K}}_\alpha }(x) + {\varepsilon _\alpha }$,其中W为神经网络权值,${{{W}}_\alpha } = {[{W_{\alpha 1}},{W_{\alpha 2}}, \cdots,{W_{\alpha m}}]^T}$,其中α = 1,2,3,4,5,$\varepsilon $为神经网络逼近误差。所以神经网络的输出为 ${\bar f_\alpha }(x) = $${{{W}}_\alpha }^T{h_\alpha }(x)$,其中${\bar f_\alpha }(x)$为不确定项${f_\alpha }(x)$的神经网络系统代替值。

      采用RBF对PID的参数进行训练时,其误差指标按式(14)进行对比。根据梯度下降法,神经网络权值按式(15)、(16)调节。

      $${{E}}(t) = \frac{1}{2}{[{y_{\rm{i}}}(t) - {y_{\rm{d}}}(t)]^2}$$ (14)

      式中:yi为输入信号;t为时间;yd为输出信号;E为误差。

      $$\Delta {{{w}}_i}(t) = - \eta \frac{{\partial {{E}}}}{{\partial {w_i}}} = \eta ({y_{\rm{i}}}(t) - {y_{\rm{d}}}(t)){h_{\alpha}}\;\;(i=1,2,3,4,5)$$ (15)

      式中:wi为神经网络;权值$\eta \in (0,1)$为学习速率;hα为径向基函数。

      $${{{w}}_{\alpha}}(t) = {{{w}}_{\alpha}}(t - 1) + \Delta {{{w}}_{\alpha}}(t) + \lambda ({{{w}}_{\alpha}}(t - 1) - {{{w}}_{\alpha}}(t - 2))$$ (16)

      式中:$\lambda \in (0,1)$为动量因子。

    • 5个电机的驱动函数由Matlab实时计算,经补偿算法在线修正后输出。在Adams中设置输入变量为丝杠位移量,具体数值由补偿算法实时求解。设置输出变量,即螺母实际位移和夹具转角,输出变量为反馈量,与机构反解位移作比对,差值由神经网络PID控制调节整定并存储在error中。在Simulink中建立控制模型如图7所示,PID的各系数由RBF神经网络在线调节,实现送料平台的精确运动,xylDlElF为运动模型的初始输入量,由逆解计算得出。

      图  7  Matlab与Adams联合仿真Simulink模型

      Figure 7.  Simulink model for simulation through Matlab and Adams

    • 假设板材加工轨迹可表示为曲线${{y}} = 15\cos (0.1x),$为了验证补偿控制的可行性,运用Matlab和Adams联合仿真的方法求解机构各输入、输出曲线。补偿后各螺母位移−时间关系曲线如图8所示,可以看出各位移都在预设行程内。补偿前后,板材加工轨迹曲线如图9所示。从图9中可以看出:经补偿控制后加工轨迹曲线与指令曲线在XY方向最大偏移量不超过1.5 mm。计算其转角误差不大于1.5°,而补偿控制前该偏移量最大为3 mm,角度误差达4.8°。在部分区段,对加工轨迹曲线补偿后,可与指令曲线重合。其中,锯切轨迹在曲线顶峰和低谷区间的误差较大。这是由于在这一曲线段及其类似部分,为使锯与板料达到相切的目的,需要机构夹具的转角达到90°,即夹具需与X轴所在直线垂直。然而机构自身的结构决定该转角无法达到90°,于是锯切该范围曲线时,只能做到近似相切,补偿后,误差虽然较大,但已能满足木工加工精度要求。

      图  8  各轴螺母补偿后位移曲线

      Figure 8.  Curves of each screw nut’s displacement after compensation

      图  9  夹具上一点运动轨迹

      Figure 9.  The trajectory of one point on the clamp

    • 选用实心细木工板(大芯板)作为试验对象,板材规格为500 mm × 300 mm × 12 mm,采用带锯对其进行锯切试验,带锯机选用MJ344,锯条规格为3 035 mm × 27 mm × 0.9 mm × R4/3。图10a为送料平台样机,其中a = 240 mm,b = 190 mm,LD4 = 240 mm,LE4 = 335 mm,LF4 = 240 mm,c = 300 mm。各丝杠螺距皆为5 mm,L1 = 300 mm,L2 = 550 mm。电机转速ɷ = 30 r/min。送料平台初始位置如图10a所示。测量并绘制锯切轨迹曲线图10b,与指令函数对比可知其最大误差不超过1.5 mm,与仿真结果一致。

      图  10  送料平台试验样机及加工轨迹曲线

      Figure 10.  Model machine of feeding platform and the curve of trajectory

    • 为提高送料平台末端姿态精度,本文按照拉格朗日法与虚功原理建立机构传递函数,通过对机构进行动力学仿真,分析机构误差来源以及各误差对平台末端姿态准确性的影响,建立相应的误差模型。设计了一种单神经元PID控制算法对平台各驱动位移进行在线补偿。通过仿真分析的方法验证了补偿策略的可行性,并用该方法对试验台进行补偿控制试验,结果表明:文中设计的单神经元PID补偿控制可有效提高送料平台末端姿态的准确性,使其轨迹在XY方向的偏移误差从补偿前的3 mm,降至小于等于1.5 mm,倾斜角误差从补偿前的3.5°降低至1.5°,该误差在木材加工领域属于允许范围。说明文中设计的单神经元PID补偿控制可以提高送料平台送料精准度,为送料平台在木材加工机床的使用方面提供了一种实用的补偿控制方法。

参考文献 (14)

目录

    /

    返回文章
    返回