在工程中的数学问题

PDF
在工程中的数学问题/2020/文章

研究论文|开放存取

体积 2020 |文章的ID 9752040 | 10 网页 | https://doi.org/10.1155/2020/9752040

建立改进的动态模型体外超声碎石合作医疗机器人

学术编辑:扬Kostavelis
收到 2019年8月8日
修订 2020年1月30日
接受 2020年2月24日
发表 2020年4月10日

摘要

结石是具有高发生率的常见的疾病之一。有效的治疗方法是体外超声碎石。目前,它是低有关碎石机的智能化,自动化水平,并已逐渐不能满足治疗需要。体外超声碎石合作医疗机器人可以有效地解决这些问题,并且它是动态模型的精确建模的研究同样重要。根据以前的研究和实验的基础上,提出了一个修正理论,以提高在协同工作的机器人模型误差的动态模型的准确性。该研究首次建立动态模型和协作机器人的实体模型,然后从实体模型减去动态模型的值,以获得修改后的方程。最后,动态模型的精度是通过修改方程改善。实验结果表明,该动力学模型修正理论是有效和转矩方程的修正后可以改善动态建模的准确性。实验结果表明,该改进的动态模型理论是有效和转矩方程的修正后可以改善动态模型的建模精度。修改后的方程已在第五次多项式的最佳校​​正效果,并且可以被用于体外超声碎石合作医疗机器人控制。

1.介绍

近年来,医疗机器人协作研究一直方兴未艾,并在医疗领域也逐渐取得了良好的应用。张某等人。总结医疗机器人产业的发展特点,分析了63主流医疗机器人企业的特点,全面梳理医疗机器人发展的总体格局[1]。这提供了什么样的机械臂在本文中使用的参考。因此,在该纸中使用的机器人是一个协作机器人。

Mihelj等人详细描述协同机器人的特征及保护,工作条件,工作环境的角度来解释它们与传统工业机器人之间的差异。他们与人[解释的协作机器人工作的特点结合起来2]。Ren等人研究了协同机器人的教与学方法[3.]。这提供了在本文中合作的机器人数据采集的想法。法西奥等。研究协作和传统机器人组件的对比模型并比较手工装配和非合作自动装配[之间的差异4,为协同机器人换探头的研究提供了参考。Zanchettin等人在collaborative robot special report中指出,机器人的灵活性和有限性阻碍了中小型企业对机器人的使用。协作机器人可以成功地解决这类问题[5],这间接支持医疗协作机器人的研究。加图索等人研究了安全运行的合作机器人的阻抗控制,并提出了新的合作机器人阻抗控制技术。控制器通过实现能源和电力限制和验证使用KUKA LWR 4+手臂[阻抗安全人机合作6]。控制算法提供了碎片机器人的障碍物避免和碰撞算法的想法,并且还提供用于高精度动态建模的研究方向。

Xiao G等人研究了一种新的多机器人协同路径规划算法,该算法采用了集中式和分布式相结合的多机器人系统架构。采用融合免疫共进化算法和人工势场法解决了全局路径规划和局部路径规划问题,提高了机器人的自适应水平[7]。这提供了用于体外超声碎石机器人的实验路径规划的方法的参考。毛泽东等人的研究,检测平台上的协作机器人的安全控制。通过软件设计,在工作区的人机合作,实现避​​免事故的发生,具有较高的安全性[8为碎片机器人的安全设计提供了参考。Yan等研究了基于动量观测的全机械臂单点接触信息的实时估计以及人机合作中的优化算法,提出了一种基于机器人控制和运动状态的信息。基于动量观测和优化算法,整个机械手单点接触信息的实时估计方法,无需借助传感器实现机械臂的外力感知[9]。这为体外超声碎石机器人的碰撞算法和路径规划提供了新的思路。Chen研究了双高斯过程的协同机器人自适应策略,使机器人能够适应目标对象的姿态变化,获得平滑的关节运动[10],它提供了碎片机器人的人机协作实验的参考方法。萧研究了轻质模块化协同机器人的碰撞检测算法,并提出了碰撞检测算法无需增加额外的传感器和改变控制系统的复杂性[11]。这提供了用于动态碰撞检测体外超声碎石机器人的一个新的理论参考。何军研究了可变负载双臂机器人的阻抗自适应控制系统,并提出了双回路阻抗可变刚度跟踪控制策略。使用合力分解成内力和外力的方法来控制所述内回路阻抗和外回路阻抗,分别。它解决了稳定性控制问题两臂合作掌握和对象对环境的相互作用力[12],它提供了用于动态建模和体外超声碎石机器人的控制的参考方法。

对于医学协同机器人的运动学和动力学,常用的运动学描述方法有D-H参数法、改进的D-H参数法、指数积方程法(POE) [13- - - - - -15]。对于动力学建模的常用方法包括拉格朗日方法[16,17],牛顿 - 欧拉方法[18- - - - - -20],凯恩方法[21- - - - - -24],旋到偶数方法[25,26]和虚功原理方法[27,28]。拉格朗日方法具有大量当机器人具有高自由度的计算,但避免了线性加速度和角加速度的计算。牛顿 - 欧拉方法太复杂的自由multidegree,但计算的综合,所以惯性的影响往往简化或忽略。凯恩方法避免了推导和计算效率高。它也以并联机器人的动态建模和高精度的控制使用的29- - - - - -31],这是便于计算机控制。它认为,广义的动态之和广义惯性力为零[27,32,33]。Zheng等人建模分析,可靠性和非线性动力学响应于多链路关节和摩擦方面[34- - - - - -37为本文的建模提供了参考。

对于本文的内容,实体模型比传统的动态模型更准确,但不能直接应用于控制系统。传统的动态模型精度较低,但可以直接用于控制系统。利用实体模型和传统动态模型的实验数据,建立了改进的动态模型。通过修正方程,提高了动态模型的精度。利用数学建模和实验方法对修正方程的待定系数矩阵进行求解。证明了体外超声碎石机器人的补体校正方程是五阶线性相关的。

2.合作机器人的组成和参数

医疗机器人协作系统主要由三个部分组成:治疗床,控制系统,和协同臂系统。该图1显示了床和双合作机械臂系统的设计。本文的动态建模研究是围绕在图中所示的6自由度机器人的协作系统中进行2。数据2(一个)- - - - - -2 (c)分别为合作机器人的前视图、侧视图和后视图。

在治疗过程中,医生确定该宝石的位置,并发送一个位置信号给控制系统。相配合的机器人然后接收位置信号,并计划用于机器人致动器的路径积极到达砾石位置并开始碎石治疗。它的动态响应要求较高的速度和准确性。本研究使用由实验室开发,如图协作医疗机器人系统2。它有六个单自由度关节。1号模型用于自底向上和一二接头。3 - 4接头采用2号模型。3号模型用于五六关节。合作机器人和关节的参数如表所示1。本文的实验和实体造型是基于这些参数。


指定 单元 1号 第2号 第3号

额定转矩 纳米 480 180 100
重复精度 ≤0.01 ≤0.01 ≤0.01
峰值扭矩 纳米 1200 700 340
重量 公斤 12.8 7.6 5.6
旋转角度 ≤360 ≤360 ≤360
减速比 1 160 160 120
最大角速度 34 74 74
额定电压 直流 48 48 48
有效载荷 13 kg
总重量 55 kg
工作半径 1400 mm
重复定位精度 ±0.02 mm
自由度 6
TCP最大速度 1.6 m/s
环境温度 −10°C∼50°C
电气接口 USB;RS232;CAN总线
电源 230v交流50 - 60hz
能量消耗 Typical operating conditions are approximately 650 W

3.建立运动学模型

根据图2,机器人的原理图如图所示3.。在图中,Z1∼Z6是第z-每个关节的笛卡尔坐标系的坐标轴。由于末端执行器的碎石机理是刚性的,动力学研究中忽略了用刚体代替该部分。

据的大小关系和坐标系统的位置,协作机器人的d-H参数列于表2。进一步得到了运动学变换矩阵。


链接 (毫米)

1 0 0 420 0
2 Pi / 2相 0 0 Pi / 2相
3. 254 0 Pi / 2相
4 Pi / 2相 0 469 0
5 Pi / 2相 0 111 0
6 -Pi / 2 0 251 0

然后,运动模型方程式如下所示:

4.建立动态模型

凯恩法是一种利用达朗贝尔原理和虚位移原理建立动态方程的方法。建立机器人机构的动力学模型是一种通用的方法。其基本思想是将广义坐标替换为系统的自变量。这意味着广义动力学和广义惯性力之和为零。凯恩方程法在动态建模中的突出优点是只需要计算向量点积和叉乘运算,避免了导数运算。因此,凯恩方法在动力学建模过程中的计算量远远小于其他方法,大大提高了计算效率,便于控制。

根据凯恩方法的计算过程中,操作的初始条件如下。

根据固体模型的参数,每个关节的旋转矩阵如下所示:

重心的位置矢量 构件描述通过杆坐标系 ,其位置矢量如下:

如示于方程每个关节的坐标系中的惯性张量矩阵(4A)- (4F)。单位是 :

根据运动学方程和结构尺寸,广义位移推断如下:

凯恩方法需要推导从接头1的广义速度和加速度值,以关节6.对于旋转接头,它使用以下的公式来计算。广义速度在等式如图所示(6)。质心的广义速度在方程式所示(7)和(8)。一般角速度如式(9)。获得广义角速度,如等式(10)。广义角加速度如式(11)。广义加速度如式(12)。质心的广义加速度如式(13):

由工况确定的初始条件如式(14)。它假定 -绝对参考帧的轴是垂直于地面和向上指向:

凯恩方法是基于所述凯恩动态方程动态方法。主动态方程示于方程(15)和(16):

没有额外的外力负载在端部作为一个简化的条件被应用。代入已知参数代入方程式(15)和(16),最终的转矩可以表达如方程来获得(17A)- (17F)后迭代计算:

在方程(17A)- (17F

整理最终扭矩表达式后,得到动力学模型的总表达式如下:

5.动态模型分析

对体外超声碎石医疗合作机器人进行路径规划,将结果代入方程(17A)- (17F)。取代的主要参数示于表3.。末端执行器的运动位移曲线如图所示4


指定 单元 关节1 关节2 关节3 联合4 联合5 联合6

广义位移 x 毫米 0 0 0 0 0 0
y 毫米 0 0 -254 0 0 0
z 毫米 420 0 0 469 111 251
质量中心 x 毫米 0 127 0 0 0 0
y 毫米 0 0 −58.46 0 58.30 0
z 毫米 -100 147 0 -85 0 0
接头质量 公斤 12.8 12.8 7.6 7.6 5.6 5.6
惯性张量 方程(4A) 方程(4B) 方程(4C) 方程(4D) 方程(4E) 方程(4F)
初始条件 , , , ,

的条件被输入到动力学模型,并且该任务计划下的动态模型的响应特性示于图5

然后,基于图中示出的动态模型参数和所述协作机器人2,在ADAMS实体模型。它的参数是完全一样的协作机器人。代在图中所示的输入条件4到实体模型,在实体模型关节的动态响应曲线示于图6。运动1〜6的运动在图6表示接头1至接头6的响应曲线。

6.改进的动态模型分析与建模

基于动态建模和实体建模,在这部分分析模型误差。首先,存在于所有的动态模型的错误,这是实际的扭矩变化,并在数学模型中的扭矩变化之间的差异。也就是说,误差等于实际值减去数学模型值。误差被定义为 ,如式(19), 是动态模型的实际矩函数吗 为数学模型中弯矩的变化量。

当动态模型的精度可以提高时,考虑修正函数 存在,所以方程(20) 成立。该 定义为动力学模型提高精度后的误差值,如式(20):

校正函数的表达 由方程(19)和(20)在等式所示(21), 为模型误差变化量:

大家都知道,动态模型的真正价值往往很难找到和过于复杂的模型公式,使模型失去了实际应用的可能性。当模型构建以获得误差比一般的数学模型下,我们认为,这是更接近真实价值,而不是为建模和计算方法,更好的模拟结果仍然可以得到的实际值。在校正功能的研究 ,我们使用实际的值 而不是真正的价值 建模和计算。这是因为我们更关心的是模型本身,而不是模型直接到达真值的精度的提高。

校正功能 与时间参数有关吗t。根据泰勒定理, 是一个n时间的八阶函数t。校正函数被定义,如等式(22), 未定系数矩阵, 是时间变化量, 为机器人当前关节自由度变化量, 是函数的顺序,并 是待确定的功能的顺序:

由此可以看出,方程(22)是动力学模型校正功能的隐函数。为了获得该功能的一个明确的溶液,需要的数学模型和实际模型来定义和解决相应的参数。

为了找到实体模型和各个关节的动力学模型之间的差异,我们减去从实体模型的响应值的动态模型的响应值。将得到的模型差异显示在图7,各关节部分误差的差值如表所示4


数量 关节1 关节2 关节3 联合4 联合5 联合6

1 2386.993 6350.648 −6500.23 -285.029 243.4248 30.51686
2 963326.4 1605860 -1.2E+ 7 -4.9E+ 7 −6322412 3456341
3. 1490851 1763153 -1.1E+ 7 -4.7E+ 7 -5900607 3269695
4 1604890 521430.1 −6484790 -4.5E+ 7 -5056274 3164968
5 1553503 -1710918 684622.4 -3.9E+ 7 −2993026 2884595
6 385449.2 -8416316 14053210 −9698177 2759736 1795244
7 -919109 -1.1E+ 7 18790044 6699360 3397692 1746705
8 -1696889 -1.1E+ 7 21712862 17946281 2836632 1935051
9 −1469893 −1E+ 7 21058280 24486418 2035476 2076084
10 −143465 -8416636 15059756 29383435 1141529 2232086

根据表的数据4,发现拟合结果好时j等于5。然后,未确定的系数矩阵 下式为:23)。基质0-5从底部到顶部与电源相关的:

根据公式(22)和(23),在每个关节的动态模型的校正函数在方程中示出(24)。它是在式所示(25),其校正函数的总表达式与接头1至关节6。此外,该改进的动态模型方程是式(26),从而完成对动态模型精度的改进:

根据表中的数据4和图56,与原始模型方程相比这样的点的校正量的改善百分比示于表5和关节模型精度的平均增加百分比示于表6。校正模型精度的综合整治率为34.10%。


数量 关节1 关节2 关节3 联合4 联合5 联合6

点1 (%) 29.84 27.11 28.82 30.33 28.32 30.24
2 (%) 33.62 33.82 30.83 34.19 34.31 34.82
点3(%) 37.89 36.81 36.82 37.18 34.87 35.83
点4(%) 35.61 37.84 37.91 35.42 36.77 36.91
点5(%) 36.31 35.42 38.33 33.03 37.72 37.33
点6(%) 33.23 30.03 35.44 32.17 36.33 34.01
7点(%) 37.22 34.37 34.08 34.17 35.13 33.23
8点(%) 34.32 32.33 33.13 34.33 32.15 33.68
点9(%) 30.82 33.05 34.73 34.43 32.78 32.14


数量 关节1 关节2 关节3 联合4 联合5 联合6

百分比(%) 34.32 33.42 34.45 33.92 34.26 34.24

7.结论

本文提出了一种体外超声碎石合作医疗机器人,它可以有效地提高动态模型的任务执行的准确性改进的动态模型。根据凯恩方法,协作机器人的动态模型被建立,并且在实体模型是基于协作机器人建造。通过比较动力学模型和实体模型之间的差异,改进的动态模型被建立用于协同的机器人。装修后,校正功能被发现是最好的,当订单j待定系数的是等于五。与原来的动态模型相比,精度由34.1%,这提供了使用所述体外超声碎石合作医疗机器人动态模型的理论模型提高。

数据可用性

支持本研究结果的数据可根据要求从通信作者处获得。

利益冲突

作者宣称,有兴趣就本文发表任何冲突。

致谢

这项研究得到了中国的甘肃省自然科学基金(18JR3RA142)和教育研究基金(GJ2019C-57)的支持。

参考

  1. 张,马,周,周,“医疗机器人技术发展综述”,发表于机电一体化和自动化2018年IEEE国际会议论文集(ICMA),第443-448,中国长春,2018月。查看在:出版商网站|谷歌学术
  2. M. Mihelj, T. Bajd, A. Ude等人,“机器人技术”,in协作机器人,施普林格,柏林,德国,2019。查看在:谷歌学术
  3. T.仁,Y栋,D吴和陈K.“的基础上力的互动协作机器人的直接教学行为设计”智能和机器人系统杂志,第96卷,不。1、2019年第83-93页。查看在:出版商网站|谷歌学术
  4. M.法乔,M. Bottin,和G.罗萨蒂,“协作和传统机器人组件:一个对比模型,”国际先进制造技术杂志,第102卷,第2期。5-8,第1355-1372页,2019年。查看在:出版商网站|谷歌学术
  5. A. M. Zanchettin,E.克罗夫特,H丁,和M.李,“在工作场所协作机器人[从特邀编辑],”IEEE机器人与自动化杂志卷。25,没有。2,第16-17,2018。查看在:出版商网站|谷歌学术
  6. G.莱奥拉,C·A·卡德纳斯,T. S. Tadele,T.德弗里斯和S. Stramigioli,“协作机器人安全和能源意识到阻抗控制器的发展,”IEEE机器人和自动化通讯卷。3,没有。2,第1237至1244年,2018。查看在:出版商网站|谷歌学术
  7. 肖国斌,严新辉,“新型多机器人协同路径规划算法”,计算机科学卷。40,没有。4,第217-220,2013。查看在:谷歌学术
  8. M.久镛,“关于检测平台合作机器人的安全控制,”电气工程卷。19,没有。4,第51-53,2018。查看在:谷歌学术
  9. H.燕,王十,刘Z.,和T.李“的基础上,气势观察和优化算法在人类与机器人合作的整个机械臂的单个联系人信息实时估计,”机器人卷。40,没有。4,第393-400 ,, 2018。查看在:谷歌学术
  10. Y.陈,郭J.和J.刘,“基于双高斯过程协同机器人的自适应策略”计算机集成制造系统,第23卷,第3期。9、1898 - 1898年,2017年。查看在:谷歌学术
  11. 肖,田,洪,“光与模组化合作机器人之碰撞侦测演算法”,天津大学学报卷。50,没有。11,第1140至1147年,2017年。查看在:谷歌学术
  12. J.他,阻抗控制器的自适应刚度双臂机器人协作处理变载荷安徽:中国科学,淮南的技术,中国,2016年大学。
  13. A.乔拉,C. Frazelle和I.沃克,“为多分连续操纵恒定曲率正向运动学的比较,”在在机器人计算IEEE国际会议论文集,拉古纳山,CA,USA,2018年1月。查看在:出版商网站|谷歌学术
  14. S. Pedrammehr,B. Danaei,H.阿卜迪,M. T. Masouleh和S. Nahavandi,“Hexarot的动态分析:轴对称平行机械手,”Robotica,第36卷,不。2018年第225-240页。查看在:出版商网站|谷歌学术
  15. A. Schweikard和F恩斯特,“机器人运动学,”在医疗机器人,施普林格国际出版,柏林,德国,2015年。查看在:出版商网站|谷歌学术
  16. B.张C.邵,李Y.,H坦和D.江“的基础上ADAMS仿人机器人步行系统的动态模拟分析,”上海交通大学学报(自然科学版),第24卷,不。1、2019年第58-63页。查看在:出版商网站|谷歌学术
  17. R. Boughdiri,H.纳赛尔,H. Bezine,N. K. M'Sirdi,A. Naamane和A M. Alimi,“在考虑耦合动力学自由运动的多指机器手的拉格朗日动力学建模,”先进材料研究,第588-589卷,第1659-1663页,2017年。查看在:出版商网站|谷歌学术
  18. 梁兵,李天成,陈志成,王玉英,廖玉英,“基于深度学习与物理仿真之机器人手臂动力学控制”,国立中央大学机械工程研究所硕士论文2018年第37届中国控制会议论文集(CCC), 2921-2925页,武汉,2018年7月。查看在:出版商网站|谷歌学术
  19. 杨,姜,何,李志,徐斌,“具有有限时间收敛之机器人控制系统之自适应参数估计与控制设计”,国立中央大学机械工程研究所硕士论文。IEEE交易在工业电子卷。65,没有。10,第8112-8123,2018。查看在:出版商网站|谷歌学术
  20. J.李,C. K.柳,F. C.公园和S. S.斯里尼瓦沙,在“A线性时间变积分器多体系统”,国际研讨会的机器人技术的基础算法程序,梅里达,墨西哥,2018月。查看在:谷歌学术
  21. 杨,“水下自重构机器人的动态模型与CPG网络生成”,先进的机器人,第30卷,不。14、2016年第925-937页。查看在:出版商网站|谷歌学术
  22. 王晓明,许振邦,王兵,“6-UHP并联平台在基底激励下的动态建模与仿真”,机器人卷。38,没有。6,第687-695,2016。查看在:谷歌学术
  23. X.李,H.孙,L.辽,和J.宋,“建立用于7自由度可重构模块化机器人的改进的凯恩动态模型,”在工程中的数学问题, 2017年卷,文章ID 4028602, 13页,2017年。查看在:出版商网站|谷歌学术
  24. k·杨,X.-Y。基于凯恩动态模型和中心模式发生器的水下蛇形机器人游动仿真平台上海交通大学学报(自然科学版)卷。19,没有。3,第294-301,2014。查看在:出版商网站|谷歌学术
  25. X.菊梁,Z.苏雷,G.浩,T.黄,和B.林,“研究为双机器人镜铣削协同加工方法中,”国际期刊先进制造技术,第105卷,不。2018年第4071-4084页。查看在:出版商网站|谷歌学术
  26. 张f,曲j,刘h,傅y,“冗余双臂机器人的位姿/力对称协调方法”,装配自动化卷。38,没有。5,第678-688,2018。查看在:出版商网站|谷歌学术
  27. S. R. Ahmadzadeh,“模态的方法来超冗余机械臂动力学和模糊控制器的系统的设计,建模” IEEE,皮斯卡塔韦,新泽西州,美国,2015年,M.S.论文。查看在:谷歌学术
  28. G.汉,F.解,和X. J.柳,“一个高速并行机器人的电力消耗的评价,”机械工程前沿卷。13,没有。2,2018。查看在:出版商网站|谷歌学术
  29. Y.韵Y.李,“设计和一种新的六自由度的冗余的分析致动具有用于高精度定位符合铰链平行机器人”非线性动力学卷。61,没有。4,第829-845,2010。查看在:出版商网站|谷歌学术
  30. “一类多自由度主动振动控制机器人的一般动力学与控制模型”,机构与机器原理卷。46,没有。10,第1549至1574年,2011。查看在:出版商网站|谷歌学术
  31. Y.云和Y. M.李,“一种基于凯恩的动态主动控制并联机构的一般模式,”在在IEEE亚太区会议上电路与系统的论文集(APCCAS),页1830-1833,中国澳门,2008年11月。查看在:出版商网站|谷歌学术
  32. 沈,吴,陈,“多自由度连续机器人之高效逆动力学建模”,中国机械工程卷。27,没有。1,第20-24页,2016。查看在:谷歌学术
  33. Y.许,Y H.张,W. H.刘,X. G.段,和L.高,“动力学建模流畅灵活的机器人关节的控制,”在在机电工程与自动化2015年IEEE国际会议论文集(ICMA),第631-636,北京,中国,2015年8月。查看在:出版商网站|谷歌学术
  34. E.郑,朱河,S朱,和X路“上的灵活的多连杆机构,包括关节有间隙和润滑的超精密冲床动力学研究”非线性动力学,第83卷,第3期。1-2页137-159,2016。查看在:出版商网站|谷歌学术
  35. Z.王,问天,H.胡和P.弗洛雷斯,“非线性动力学和不确定关节间隙柔性多体系统的混沌控制,”非线性动力学卷。86,没有。3,第1571至1597年,2016。查看在:出版商网站|谷歌学术
  36. 张智雄,“具有多重不确定性之平面并联机械臂之误差建模与运动可靠性分析”,国立中央大学机械工程研究所硕士论文。机构与机器原理,第124卷,第55-72页,2018年。查看在:出版商网站|谷歌学术
  37. 十,王,刘G.,S.马,和R.堂“研究用干回转间隙铰平面多体系统动态响应:数值模拟和实验方法”声音与振动杂志的卷。438,第116-138,2019。查看在:出版商网站|谷歌学术

版权所有©2020小李等人。这是下发布的开放式访问文章知识共享署名许可,其允许在任何介质无限制地使用,分发和再现时,所提供的原始工作正确的引用。

0 查看 | 0 下载 | 0 引文
PDF 下载引用 引文
下载其他格式更多的
订单打印副本订单

我们致力于快速,安全地与COVID-19尽可能共享成果。任何作者提交COVID-19纸应该通知我们的help@hindawi.com以确保他们的研究是快速跟踪和尽快预印本服务器上公布。我们将针对与COVID-19接受的文章中提供的出版费用减免无限。