数控机床是装备制造业的工作母机，其可靠性技术目前已成为制约行业发展的关键共性技术[1]。现阶段检测数控机床可靠性的主要方法是“经时试验”，即使用该机床长时间切削金属坯料，使机床主轴承受切削力，以此发现机床长期运行时所发生的故障。但是，该方法耗费大量的原料和刀具，显著增加了测试成本。所以，需要研制一种可向机床主轴施加多维载荷的加载装置，以改进传统测试方法。

并联机构由于其结构紧凑、承载能力强、积累误差小和动态性能好等特点在工业领域被广泛使用[2]。其动平台具备空间多自由度，故可以将其应用于多维力加载领域。关于并联加载装置的研究，最早可追溯到1962年Gough 和Whitehall发明的6自由度并联机器人[3],将其应用于对轮胎的加载中，测试轮胎在不同载荷下的性能。徐彬[4]采用Stewart平台结构的形式，对大型旋挖钻机动臂加载性能参数进行了分析。刘少欣等[5]对6-UPS并联加载装置的动力学模型进行了简化。王博和黄其涛[6]提出了一种正交并联机构的多轴加载试验系统，用以模拟对接机构组合体实际应用中受到的复杂载荷，是并联机构一种新的应用形式。Stokes等[7]基于6自由度并联机构设计一种针对脊柱的加载装置，在末端安装多维力传感器，对动物脊柱进行了刚度测量。Nierenberger和Flohic等[8-9]将Stewart平台用作多维力材料加载试验机，并对混凝土试验件进行拉力-剪力-剪切力二维加载试验。Guo等[10-12]同样基于Stewart平台，设计多维力加载试验机，对试验机的输出载荷进行检测。但是，上述研究采用的Stewart平台的驱动电机位于运动支链中，随连杆一起运动，使机构运动部件的质量增大，对机构的动态特性产生不利影响。相比之下，6-PSS并联机构由于其驱动电机固定，可有效减小运动部件的质量，机构的动态特性好，便于在多维加载过程中实时跟随目标的空间运动。王洪瑞等[13]建立了6-PSS并联机器人的动力学模型；孙小勇等[14]对6-PSS并联机器人的参数进行了优化设计；Wang等[15]基于Hexaglide型并联机构设计多维力加载装置，并对其刚度分布进行分析，但是并未对加载装置进行多维力加载试验。

应当把这种做好事的动力和传统儒家文化中提倡的在自律基础上的弘仁、达人联系起来。做什么事情通达无愧于心的境界就好了。做好事是非功利的，目的在于心安。不如此就难以解释南通当代有许多人自觉长期照顾那些与自己毫无血缘关系的人员的奉献之举，如沈亚秋长期照顾三对与自己并无任何血缘关系的老人。

本文提出一种基于6-PUS并联机构的多维力加载装置，可对目标提供多维载荷。首先根据向量叠加原理和牛顿迭代法分别介绍了6-PUS并联机构的逆、正解运算，并由旋量理论建立了机构的运动学和静力学模型；然后根据虚功原理建立了机构的动力学模型并进行仿真验证；最后对多维力加载装置进行多维力加载试验。

1 6-PUS并联机构的运动学建模

1.1 运动学逆解

如图1所示，6-PUS并联机构由动平台、静平台及6根完全相同的支链组成。其中，动平台复合球铰中心Bi(i=1,2,…,6)的分布半径为rp，6个球铰链分为3组，沿圆周均布，每组之间夹角为120°；6个直线驱动及静平台复合虎克铰Ci(i=1,2,…,6)同样分为3组在静平台上沿圆周均布，每组驱动的安装平面倾角为β。将滑块置于零位，每组两驱动的中点的分布半径为 rb，各支链中的连杆长度为L。

选任一支链,使用向量叠加描述该运动链得

(1)

式中：为各运动支链中移动副的驱动量；Li为定长连杆的长度；si为静坐标系中移动副的方向向量；ki为静坐标系中沿连杆方向的单位向量；p和R分别为动坐标系O′{x′,y′,z′}相

图1 多维力加载装置结构及坐标系定义示意图 Fig.1 Schematic of structure and coordinate system definition of tested multi-dimensional loading device

对于静坐标系O{x,y,z}的原点位置和坐标系旋转矩阵；在动坐标系中定义

令根据机构构型和运动的连续性求解式(1)得

(2)

1.2 运动学正解

与运动学逆解相反，多维力加载装置的运动学正解是已知6个运动支链的电机驱动量，求解机构动平台位姿的过程，即已知驱动量di(i=1,2,…,6)，求解动平台中心点O′的位置向量p=[px py pz]T和姿态向量Θ=[θx θy θz]T，其中向量p和Θ的3个分量分别是动平台沿x、y、z轴的位移和绕其轴线的转角。

将动平台的位姿表示为Q=[pT ΘT]T,可以根据姿态向量求解出旋转矩阵R=f(Θ)，由运动学反解可知

(3)

i6

(4)

则Q=[p A]T，A=[a0 a1 a2 a3]T。

将R代入式(3)，再与式(4)相加可得7个方程，最后利用牛顿迭代法即可求解得出位置向量Q的7个分量未知数。求解过程如图2所示，首先确定动平台初始位姿Q0，然后用牛顿迭代法求出迭代位姿Qn(n为迭代次数，n≥0),当由末端迭代位姿Qn反解计算出的杆长与已知杆长的差值向量Δf=[f1 f2 … f7]的模满足f(Qn)<ξ=10-10 m时，即可认为此时的迭代位姿Qn可以作为当前动平台的实际位姿向量Q，迭代结束。

图2 机构运动学正解流程图 Fig.2 Flowchart of forward kinematics of mechanism

2 运动学与静力学分析

2.1 速度及静力雅可比矩阵分析

任取一条运动支链作为研究对象，如图3所示。将动平台中心点的速度旋量用Plücker 坐标表示为T=[vT ωT]T，其中v和ω分别表示动平台质心点相对静坐标系的平移速度和旋转角速度。

由旋量法可得

一身空姐职业装的程凤萍今天接受人民交通网专访的时候，给记者的印象是“干练、热情、泼辣、大方”，难怪南航湖南分公司的同事们都亲切地称她“凤姐”。

式中：旋转矩阵R用四元数表示，相比欧拉角，可有效避免奇异问题。通常，可将四元数定义为A=a0+a1i+a2j+a3k,(ai∈R,i=0,1,2,3)，且

二是建立岗位业务提升机制。整合可利用的培训资源，结合业务工作实际加强党员业务知识的学习培训，开阔党员同志的视野，充实、调整和完善知识结构，提升业务岗位技能和工作能力。通过建立和落实党员业务服务调研制度，进一步提升岗位服务的实效性。重视培育良好的工作态度、工作作风，打造一支特别能奉献、特别能战斗的高素质党员队伍，使党组织的服务能力长效可行。

(5)

将驱动量整理为矩阵形式为

将式(5)的两端同时与单位速度旋量取互易积，且 i∘ ij=0(j=2,3,…,6)，整理为矩阵形式可得

(6)

式中：动平台质心速度滑块移动速速度雅可比矩阵(i=1,2,…,6)。

由速度和静力学对偶关系可知：

缓斜井断面型式采用马蹄型断面，断面净尺寸5.0m×5.958m（宽×高，半圆城门洞型）。缓坡段设置错车道，断面净尺寸7.1m×5.958m。纵断面组成为200m（12.9%）陡坡段+20m（3%）缓坡段，综合纵坡12%，井底布置30m长纵坡为零的洞段，缓斜井长640m。洞身内布置有行车道、2条直径1.0m的风管、各种管路和供电线路等设施。电缆及管子布置在洞壁上，电缆采用挂钩固定在井壁上，位于管子上方，管子布置在缓斜井底板的管座上。缓斜井横断面如图3。

(7)

式中：F=[fT mT]T为施加在动平台质心的外力载荷f为外力，m为力矩；τ=[τ1 τ2 … τ6]T为移动副滑块驱动力向量；Jf为力雅可比矩阵。

图3 支链坐标系及速度旋量分布示意图 Fig.3 Schematic of branched-chain coordinate system and velocity screw distribution

2.2 速度分析

选取任一定长连杆为研究对象，以Ci为原点,建立如图3所示的支链坐标系Oi{xi,yi,zi},其初始状态x、y和z轴分别与静坐标系的3个坐标轴平行，故连杆方向向量可使用ZYZ欧拉角表示，即支链坐标系Oi先绕静坐标系z轴旋转角φi，再绕新坐标系的轴旋转角度θi，此时iki=zi。

对式(1)的右侧求导可知：

(8)

式中：

将vbi通过坐标变换转换至支链坐标系得

(9)

式中：为支链坐标系相对静坐标系的旋转矩阵；iJbi=[iJbixT iJbiyT iJbizT]T。

[4] 徐彬.基于并联机构的多维力加载试验系统研究[D].哈尔滨： 哈尔滨工业大学,2013:10-40.

(10)

对式(10)的两端同时点乘iki，可得

式中：和分别为转动副旋转轴方向向量和旋转速率，且ki6=ki，si和分别为移动副的方向向量和移动速率。

(11)

式中：

对式(10)两端同时叉乘iki，由iωi·iki=0可得连杆角速度在支链坐标系中的表示为

在曲线面板的左上角，我们可以看到一个靶心形状的小图标。单击图标即可以选择目标调整工具，使用该工具在画面任意位置上下拖动，就等于直接在曲线上上下调整该区域所对应的影调范围。这个方法相比起直接编辑曲线，又更简单一些。

木材热处理改性后，其粗糙度的变化规律一直得到各界广泛的关注。在前人的大量文献中就可以看出，木材热处理后其粗糙度的变化规律若得到更细致更明确的总结，对板材胶合、涂饰工艺有着一定的参考价值。本文将在前人研究的基础上，重新对4种木材进行热处理实验，并将实验过程进行记录，对木材试样热处理前后粗糙度进行测量及统计，并结合扫描电镜图，以便得出木材试样的热处理粗糙度变化规律。

式中：iΩki=iki。

由滑块的质心向量rci=ai+disi可得其质心速度为

1.1.2 排除标准 ①Ⅰ型糖尿病患者；②合并低血压、酮症酸中毒等较为严重的糖尿病并发症患者；③合并重症感染、外伤、心肺肝肾功能损伤的患者；④临床研究依从性较差的患者；⑤合并精神疾病、意识障碍、心理障碍、听力障碍或由于其他病理、客观因素影响而无法进行正常沟通的患者[2]。

(12)

同理，由连杆的质心向量rli=ai+disi+Liki/2可得连杆的质心速度，将其转换至支链坐标系：

(13)

使用Plucker坐标表示滑块的质心速度和连杆质心速

(14)

(15)

式中：

Jli2=

首先，选拔优秀合同管理人才。可依照合同管理人员应具有的素质条件，选择本企业优秀人才担任合同管理人员，也可以通过公开考评和竞争招聘方式引进、选拔人员。在三、四级要推行专兼职合同管理员制度，在使用过程中坚持优胜劣汰的原则，把优秀人才放在合同管理岗位上。

2.3 加速度分析

对式(10)求导且两端同时点乘iki可得滑块的加速度在支链坐标系中为

(16)

对式(10)求导且两端同时叉乘iki可得连杆角加速度在支链坐标系中为

(17)

对式(12)求导得滑块质心加速度向量为

一方面，利用多种形式宣传造势企业管理制度。移动网络的普遍，给中小企业领导者开辟了新的管理方式。可以利用移动网络形式开展管理制度学习的比赛，或收集对企业管理制度的意见和对企业组织发展的想法建议等，以求企业员工理解管理制度，明确制度的意义，认可制度，使得企业管理工作正常进行，管理制度能够得以贯彻落实。

(18)

对式(13)求导得连杆质心加速度向量，并将其转换至支链坐标系得

(19)

3 动力学分析

3.1 动力学方程

由虚功原理可知，多维力加载装置各部件所受广义主动力及广义惯性力合力所做的功为零。

静坐标系下动平台中心点广义主动力和广义惯性力的合力(包含力fp及力矩np)为

式中：fe和ne分别为动平台中心点外力和外力矩；mp为动平台质量；AIp为动平台在静坐标系下的转动惯量。

同理，静坐标系下各滑块质心点的广义主动力和广义惯性力合力(包含力fci及力矩nci)为

2) 6-PUS并联机构具备对目标进行多维力加载的能力，加载试验结果表明，该机构可对目标施加六维三角波和正弦载荷信号，最大跟随误差均小于10%，证明其可以对目标进行匀速和变加速多维力加载，为数控机床加载试验提供新的加载手段和理论基础。

式中：mli为连杆质量；iIli为支链坐标系下连杆的转动惯量。

由虚功原理得

(20)

当各构件运动时间相同时，δq=Jpδxp, δxci=Jciδxp, δixli=iJliδxp , 将此式(3)代入式(20)得

(21)

3.2 动力学仿真验证

首先，建立Adams理论模型，并结合模型的质量参数在MATLAB中推导多维力加载装置的动力学方程，建立动力学方程；其次，同时设定Adams和MATLAB中的机构模型的动平台质心运动轨迹如式(22)所示：

xp=[40cos(πt),40sin(πt),5t,20cos(πt),

1403 自制腹膜反折悬吊装置在腹侧肾肿瘤后腹腔镜下肾部分切除术中的应用 王 磊，储传敏，干思舜，叶剑青，曲发军，杨 炜，田毅君，潘秀武，杨启维，崔心刚

20sin(πt),30sin(πt)]T

(22)

则其质心速度和加速度向量也可求出；再次，同时设定Adams和MATLAB中的机构模型的动平台质心点所受外力为沿z轴负方向1 000 N；最后，在Adams中测量6个滑块的驱动力，与MATLAB计算得出的驱动力进行比较，并计算二者之间的误差，如图4所示。由图4可知，Adams机构理论模型与MATLAB动力学模型误差维持在2 N以内，小于驱动力的0.5%，支链2、3和5在0.04 s时误差达到4 N，但仍小于驱动力的1%，验证了本文运动学及动力学分析的正确性。

图4 Adams与MATLAB模型所得各支链驱动力的误差 Fig.4 Deviation of branched-chain actuation forces derived from Adams and MATLAB model

4 多维力加载试验

进行多维力加载试验时，将多维力加载装置放置在三轴联动龙门式机床试验平台上，动平台中心与机床主轴连接，主轴处于静止状态，改变动平台的位姿，其与主轴的相对运动会对应产生作用力，实现对主轴的多维力加载。具体加载过程如下：第一,设定期望动平台输出的多维力值；第二,由光栅尺反馈各支链的驱动量，正解得机构动平台位姿；第三,根据动平台位姿和动力学方程求解得出施加期望多维力所需要的各支链驱动力，并由此推导得出各传感器期望力；第四,使用模糊PID控制算法[16]对各支链进行力闭环控制，使动平台输出期望多维力。

多维力加载试验分为2部分：六维三角波载荷和六维正弦载荷，分别代表对机床主轴匀速加载和变加速度加载2种形式。试验采用“双十指标”评估加载性能，即实际加载力与参考加载力误差小于10%，相位误差小于10°。通过比较各方向实际载荷和参考载荷间的跟随误差，分析和评估多维力加载装置的实时性和准确性。

水库枢纽由大坝、溢洪道、输水洞、泄洪洞和水电站等五部分组成。大坝为浆砌石重力墙式堆石坝，最大坝高69.5 m，坝顶高程1 031.5 m，坝顶长206 m，宽4 m；是岚漪河上的一座以防洪为主，兼顾灌溉、防洪等任务的中型水库。

1) 六维三角波载荷加载

同时在主轴上施加幅值为200 N的三维三角波力和幅值为10 N·m三维三角波力矩，得到理论力、理论力矩及加载误差如图5所示。

观察图5可知，多维力加载装置对主轴施加的六维三角波载荷中，沿x、y和z轴的加载力误差最大值分别为4.680、2.636和5.264 N，基本处于参考加载力最大幅值的5%以内。另外，沿三轴加载力矩误差最大值分别为0.652、0.515和0.321 N·m，小于参考加载力矩最大幅值的10%，满足加载需要。

2) 六维正弦载荷加载

随着时代的发展，信息技术的进步，也为了满足现代化教学的需求，一些多媒体工具被运用到高中数学教学中来弥补传统教学的不足，替代一些落后的教学方式供教师所用.它与传统教学方式不同，能够丰富教师的教学内容，使学生接触到更多教材之外的知识点，同时也有效地达到了吸引学生的作用，使得学生更加积极主动地参与到学习活动中.但同时，多媒体教学也并不一定就能够发挥积极正向的作用，最重要的是看我们教师如何对其进行运用.

同时在主轴上施加幅值为200 N的三维正弦力和幅值为10 N·m的三维正弦力矩，得到理论力、理论力矩及加载误差如图6所示。

由图6可知，多维力加载装置对主轴施加的六维正弦载荷中，沿x、y和z轴加载力误差最大值分别为18.291、22.417和19.201 N，基本维持在参考载荷的10%以内。另外，沿三轴的加载力矩误差最大值分别为1.583、2.062和1.607 N·m，在参考加载力矩的10%以内，满足加载需要。

图5 六维三角波载荷加载误差 Fig.5 Loading errors along six axes with triangular wave load

图6 六维正弦载荷加载误差 Fig.6 Loading errors along six axes with sine load

5 结 论

本文提出一种基于6-PUS并联机构的多维力加载装置，可对机床主轴施加多维载荷。建立了加载装置的运动学、静力学以及动力学模型，并进行了动力学仿真验证和多维力加载试验，得出以下结论：

1) 采用虚功原理可对6-PUS并联机构进行准确的动力学建模。用MATLAB建立动力学模型，并在Adams中建立结构相同的虚拟样机，通过给定相同的动平台运动轨迹和外载荷，对比Adams和MATLAB仿真所得的驱动力，误差小于1%，验证了动力学建模的正确性。

式中：mci为滑块质量。同理，支链坐标系下各连杆质心点的广义主动力和广义惯性力的合力(包含力fli及力矩nli)表示为

根据吉林省6所典型高职院校旅游管理专业核心课程设置情况以及对省内其他高职院校负责人进行面访调研，吉林省高职院校总体上兼顾了“高等技术”与“应用型”这两个对人才需求的“关键词”，但是也存在一些诸如核心课程设置宽泛、地域性不突出、核心课程授课教师自身实践技能匮乏等问题。

参考文献 (References)

[1] 杨兆军,陈传海,陈菲,等.数控机床可靠性技术的研究进展[J].机械工程学报,2013,49(20):130-139.

YANF Z J,CHEN C H,CHEN F,et al.NC machine tool reliability technology database applications[J].Journal of Mechanical Engineering,2013,49(20):130-139(in Chinese).

[2] 郭江真,王丹,樊锐,等.3PRS/UPS冗余驱动并联机器人刚度特性分布[J].北京航空航天大学学报,2014,40(4):500-506.

GUO J Z,WANG D,FAN R,et al.Stiffness characteristics distribution of 3PRS/UPS parallel manipulator with actuation redundancy[J].Journal of Beijing University of Aeronautics and Astronautics,2014,40(4):500-506(in Chinese).

[3] GOUGH V E,WHITEHALL S G.Universal tyre test machine[C]∥Proceedings of the FISITA Ninth International Technical Congress.London:Institution of Mechanical Engineers,1962:117-137.

对式(1)等号左侧求导且将其转化至支链坐标系中得

XU B.Research on multi-DOF loading system based on parallel mechanism[D].Harbin：Harbin Institute of Technology,2013:10-40(in Chinese).

[5] 刘少欣,王丹,陈五一.6-UPS并联多维力加载装置的动力学建模及简化[J].机床与液压,2014,42(19):5-8.

LIU S X,WANG D,CHEN W Y.Dynamics modeling and simplification of a 6-UPS parallel multi-dimensional loading device[J].Machine Tool and Hydraulics,2014,42(19):5-8(in Chinese).

[6] 王博,黄其涛.正交并联多自由度加载试验系统研究[J].机电工程,2016,33(11):1334-1338.

WANG B,HUANG Q T.Loading test system of orthogonal parallel multi degree of freedom[J].Mechanical and Electrical Engineering,2016,33(11):1334-1338(in Chinese).

[7] STOKES I A,GARDMEMORSE M,CHURCHILL D,et al.Measurement of a spinal motion segment stiffness matrix[J].Journal of Biomechanics,2002,35(4):517-521.

[8] NIERENBERGER M,PONCELET M,PATTOFATTO S,et al.Multiaxial testing of materials using a stewart platform:Case study of the Nooru-Mohamed test[J].Experimental Techniques,2012,38(2):74-83.

[9] FLOHIC J L,PARPOIL V,BOUISSOU S,et al.A 3D displacement control by digital image correlation for the multiaxial testing of materials with a Stewart platform[J].Experimental Mechanics,2014,54(5):817-828.

[10] GUO J Z,WANG D,FAN R,et al.Development of a material testing machine with multi-dimensional loading capability[J].Journal of Advanced Mechanical Design Systems & Manufacturing,2016,10(2):JAMDSM0017.

[11] GUO J Z,FAN R,WANG D,et al.Development of multidimensional loading material testing machine[J].Materials Science Forum,2013,770:164-169.

[12] 樊锐,李茜,王丹.6PUS并联机构的运动学整机标定[J].北京航空航天大学学报,2016,42(5):871-877.

FAN R,LI X,WANG D.Integral kinematics calibration of 6PUS parallel mechanism[J].Journal of Beijing University of Aeronautics and Astronautics,2016,42(5):871-877(in Chinese).

[13] 王洪瑞,马丽,肖金壮.一种6-PSS并联机器人动力学模型的建立[J].中国高新技术企业,2010(10):19-21.

WANG H R,MA L,XIAO J Z.A 6-PSS parallel robot dynamic model[J].China High-Tech Enterprises,2010(10):19-21(in Chinese).

[14] 孙小勇,郑彬,鲍捷,等.高速6-PSS并联机器人参数优化设计[J].农业机械学报,2015,46(5):372-378.

SUN X Y,ZHENG B,BAO J,et al.Parameter optimization design of high speed 6-PSS parallel robot[J].Transactions of the Chinese Society of Agricultural Machinery,2015,46(5):372-378(in Chinese).

[15] WANG D,FAN R,CHEN W.Stiffness analysis of a Hexaglide parallel loading mechanism[J].Mechanism & Machine Theory,2013,70(6):454-473.

[16] 樊锐,刘力军,王丹,等.6-PUS并联加载机构模糊PID力控制系统设计[J].航空制造技术,2016,513(18):34-39.

FAN R,LIU L J,WANG D,et al.Design of fuzzy PID force control system based on 6-PUS parallel machine[J].Aeronautical Manufacturing Technology,2016,513(18):34-39(in Chinese).