CN102279101A - 六维力高频疲劳试验机及使用方法 - Google Patents

六维力高频疲劳试验机及使用方法 Download PDF

Info

Publication number
CN102279101A
CN102279101A CN2011101955232A CN201110195523A CN102279101A CN 102279101 A CN102279101 A CN 102279101A CN 2011101955232 A CN2011101955232 A CN 2011101955232A CN 201110195523 A CN201110195523 A CN 201110195523A CN 102279101 A CN102279101 A CN 102279101A
Authority
CN
China
Prior art keywords
side chain
control
theta
power
centerdot
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.)
Granted
Application number
CN2011101955232A
Other languages
English (en)
Other versions
CN102279101B (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.)
Beihang University
Original Assignee
Beihang 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 Beihang University filed Critical Beihang University
Priority to CN201110195523.2A priority Critical patent/CN102279101B/zh
Publication of CN102279101A publication Critical patent/CN102279101A/zh
Application granted granted Critical
Publication of CN102279101B publication Critical patent/CN102279101B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)

Abstract

本发明一种六维力高频疲劳试验机,它包括6-UPS并联机构和控制系统两个部分,所述6-UPS并联机构由加载平台、拉压力传感器、筒式直线电机、虎克铰链、静平台及球铰链组成;该控制系统由系统硬件、系统控制方法及软件两部分组成。该六维力高频疲劳试验机以动载荷机械零件的疲劳试样为测试对象的六维力高频疲劳试验机及使用方法,可用于模拟动载荷机械零件的实际受力情况,测定其在静载荷、动载荷及交变疲劳试验中的变形情况,以解决现有疲劳试验机不能施加高频多维力的不足,以接近于实际工作状态下的受力状况,准确的测定试件的疲劳承载能力。

Description

六维力高频疲劳试验机及使用方法
(一)技术领域
本发明提供一种六维力高频疲劳试验机及使用方法,它特指一类具有三维移动和三维转动的六自由度疲劳加载并联机构,属于试验机技术领域。
(二)背景技术
在实际的生产生活中,很多装备及设备都受到疲劳破坏的影响,以飞机为例。飞机是一种大型、复杂的先进装备。飞机机体结构特别是机翼等关键部件的寿命决定了飞机的总寿命,因此,飞机结构的使用寿命是决定飞机使用寿命的基础,飞机结构使用寿命的评定对飞机使用寿命评定起着决定性作用。现代飞机的结构设计要满足多种要求,而其中飞机结构强度要求和疲劳寿命是飞机结构设计的关键参数。据不完全统计,自2001年以来,世界上因机械故障造成机毁人亡的恶性事故超过15起,人员死亡超过1600人。
从结构强度和寿命设计技术的角度看,更新传统试验设备、探索新试验方法的研究亟需加强。在进行动载荷机械零件设计时,一般都是根据试件的标准力学性能测试结果、在理论计算的基础上而展开的。迄今为止,标准力学性能测试的基本载荷形式是单纯拉伸、压缩、弯曲以及扭转等。而试样的实际受力情况远非如此简单,往往是多种载荷的组合,也就是说,试样工作环境与设计依据之间存在较大差别。其结果将直接影响设备的使用寿命及安全系数,这也是目前制约结构合理化设计的一个方面。
目前,虽然可以通过实体三维造型,采用有限元仿真分析软件进行试样的受力状况进行分析,由于仿真软件自身的局限性(各种仿真软件均有其对某领域的专长)以及约束条件建立困难等原因,一般只用来对所设计试样进行校核或为改进提供参考,不可能完全替代实验。同时,由于各种动载荷机械零件需要进行大量严格的实物测试试验,因此急需一种能够模拟实际受力情况的测试装置,建立尽量接近实际工作条件并测评试样的工作状态,为动载荷机械零件的结构优化、安全可靠性设计提供依据。
(三)发明内容
本发明提出一种以动载荷机械零件的疲劳试样为测试对象的六维力高频疲劳试验机及使用方法,可用于模拟动载荷机械零件的实际受力情况,测定其在静载荷、动载荷及交变疲劳试验中的变形情况,以解决现有疲劳试验机不能施加高频多维力的不足,以接近于实际工作状态下的受力状况,准确的测定试件的疲劳承载能力。
本发明一种六维力高频疲劳试验机,它包括6-UPS并联机构和控制系统两个部分(如图1),其特征在于:
1.6-UPS并联机构
该六维力高频疲劳试验机使用6-UPS并联机构,其是6-SPS并联机构的变形。6-SPS并联机构由6个运动链连接动平台和静平台组成,并以移动副作为输入,通常称为Stewart平台。六维力高频疲劳试验机使用6-SPS的变形机构6-UPS机构,其与6-SPS机构相同,由动平台,静平台及加载支链三部分组成,不同是6-SPS并联机构的加载支链两端由球铰链分别与动平台和静平台连接。而6-UPS并联机构的加载支链两端分别由球铰链和虎克铰链与动平台和静平台连接。6-SPS并联机构的每一根杆要有六自由度,因为与动平台相连的关节轴承提供了三个转动自由度,移动副提供了一个移动自由度,故静平台一侧的关节轴承只需提供两个转动自由度即可。因此静平台采用了两个转动自由度的虎克铰链。这样即形成了6-UPS并联机构。
该6-UPS并联机构由加载平台1、拉压力传感器2、筒式直线电机3、虎克铰链4、静平台5及球铰链6组成,如图2;
其中,所述加载平台1,即6-UPS并联机构的动平台,其由六个球铰链6支撑。高刚度的加载平台在安装夹具后可以对动载荷机械零件件进行加载。
由球铰链6、拉压力传感器2及筒式直线电机3和虎克铰链4组成加载支链,通过筒式直线电机3的伸缩运动来完成加载支链的长度变化,通过筒式直线电机3与球铰链6和虎克铰链4共同作用下改变加载平台1的位姿和加载力的大小。所述筒式直线电机3作为试验机的驱动部件,也是移动副。筒式直线电机3具有高频,高加速度,低惯量的优点,是疲劳加载试验优秀性能的前提。
所述静平台5作为所述6-UPS并联机构的基座,通过虎克铰链4与加载支链连接,具有高强度,高刚度的性能。
所述静平台5和加载平台1以6个加载支链相联,每个加载支链两端分别是球铰链6和虎克铰链4,中间是移动副(筒式直线电机3)。筒式直线电机3在驱动器作用下作相对移动,改变加载支链的长度,使加载平台1的位置和姿态发生变化。
该六维力高频疲劳试验机使用高性能筒式直线电机3作为移动副,使6-UPS并联机构实现位置变化及加载运动。通过六根加载支链的运动组合使得加载平台(动平台)产生六自由度的加载运动,可以对被测试样施加六维载荷(三维力和三维力矩)。六根加载支链的运动由计算机控制,可以根据所需要载荷谱产生给定的运动组合及其变化规律。
动载荷机械零件的疲劳试样通过夹具接受载荷,可以实现三维力和三维力矩的高频力学性能测试。动平台的移动由六个加载支链带动实现。加载支链部件在自身伺服主机的驱动下,实现精确的轴向伸缩。六根加载支链的伸缩组合,实现对上夹具的各种运动驱动,进而对动载荷机械零件的疲劳试样施加各种单个或组合载荷。
该六维力高频疲劳试验机的工作原理是:每根加载支链通过筒式直线电机3产生伸缩运动,六根加载支链的运动组合使得加载平台1产生六自由度的加载运动,可以对被测试件施加六维载荷(三维力和三维力矩)。六根加载支链的运动由计算机控制,可以根据所需的载荷谱产生给定的运动组合及其变化规律。
2.控制系统
该控制系统由系统硬件、系统控制方法及软件两部分组成;
(1).系统硬件
该系统硬件包括工控机(PC)、运动控制卡、筒式直线电机3及拉压力传感器2,如图3。工控机和运动控制卡组成六维力高频疲劳试验机控制部分,筒式直线电机3和拉压力传感器2安装在6-UPS并联机构上,如图2。
该工控机采用全钢机箱、无源底板,工业电源及全长主板(2.0G Pentium IVCPU,2G内存)。
该运动控制卡采用6轴以上可编程多轴运动卡。其单轴最快的采样周期为60μs(微秒)。该运动控制卡提供运动控制、离散控制、内务处理、同主机的交互等数控的基本功能。该运动控制卡可以同步控制多个驱动轴,实现复杂的多轴协调运动。
该筒式直线电机3采用大推力高加速度筒式直线电机。该筒式直线电机3安装在6-UPS并联机构的加载支链上,一端连接拉压力传感器2,另外一端连接虎克铰链4,为六维力高频疲劳试验机提供动力。该筒式直线电机3是三相无轴承水冷同步伺服直线电机,其最大推力不小于850N(牛),持续推力不小于550N(牛),最大功率不小于1400W(瓦),电机空载最大加速度不小于350m/s2(米/秒2)。
该拉压力传感器2采用柱式高精度拉压力传感器。该拉压力传感器2一端与球铰链6连接,另外一端与筒式直线电机3连接。该拉压力传感器2的弹性体为柱式结构,用于拉伸力和压缩力测量。全密封结构,其输出灵敏度为1.5mV/V(毫伏/伏),直线度为±0.05%F.S(满量程),重复性为±0.05%F.S(满量程),最大量程为±5000N(牛)。
控制系统的核心为运动控制卡,工业控制计算机负责信息流和数据流的管理,以及从传感器、编码器读取力和位置数据,并经过计算后发送控制指令。驱动器负责直线电机的功率驱动,实现位置、速度和力的控制。
(2).系统控制方法及软件
系统控制方法及软件由系统控制、反馈系统控制及系统软件三部分组成。
6-UPS并联机构的控制包括运动学和动力学控制。运动学控制指的是位置控制,动力学控制是指力控制或力和位置的混合控制。对于运动学方法,我们需要知道机构的运动学正解和运动学逆解。当控制机构运动时,首先根据末端执行器的位置通过运动学逆解计算出各关节的位置,然后控制关节到达指定位置。动力学控制是指综合考虑机构的运动学和动力学模型生成控制指令。由于通过简单的运动学逆解运算即可生成运动指令,运动学控制简单、可靠易实现,因此应用广泛。并联机床的运动呈现高度的非线性和强耦合,使得并联机床的运动尤其在高速运动时会对运动控制精度产生较大影响,因此需要建立动力学控制系统来对并联机床进行控制。
a.系统控制
该六维力高频疲劳试验机的控制使用力伺服控制。该力伺服控制通过模糊PID控制器来实现。目前机器人的力伺服控制主要有以下几类:阻抗控制,力/位混合控制,自适应控制和智能控制。模糊控制是智能控制中的一种,模糊控制的优点在于对非线性系统有良好的控制作用。图4即为模糊逻辑控制器流程图,其是模糊控制的核心,它是根据各种边界条件和一些重要的专家规则确定的。
在工程实际中,应用最为广泛的调节器控制规律为比例、积分、微分控制,简称PID控制。它以其结构简单、稳定性好、工作可靠、调整方便而成为工业控制的主要技术之一。当被控对象的结构和参数不能完全掌握,或得不到精确的数学模型时,控制理论的其它技术难以采用时,系统控制器的结构和参数必须依靠经验和现场调试来确定,这时应用PID控制技术最为方便。而PID控制的不足是当系统存在非线性,噪声及强耦合时无法解决稳定性及准确性之间的矛盾,并且PID控制较难取得令人满意的动态响应。
模糊PID控制器结合传统PID控制和模糊控制的优点,稳定性强,在线控制方便且能适应模型参数的变化,因此该六维力高频疲劳试验机的控制——力伺服控制采用模糊PID控制器来实现。模糊PID控制器对各个加载支链的驱动力进行实时跟踪,从而实现对6-UPS并联机构的力伺服控制。
具体地说,首先基于6-UPS并联机构各加载支链的动力学方程建立其控制模型,然后采用模糊PID控制器对加载支链驱动力进行跟踪,从而保证6-UPS并联机构在整个运动过程中具有稳定而准确的驱动力输出。
要进行力伺服控制,所以先介绍一下6-UPS并联机构的受力分析。在静力分析时,可以应用螺旋理论建立起6-UPS并联机构的静力平衡方程,并转化为用影响系数表示的矩阵方程。对于该试验机,机构的位形是一定的,六根杆长l1、l2、l3、l4、l5、l6是确定的,它的动平台经由6个加载支链与静平台相连,每杆两端分别为球铰链和虎克铰链,此6-UPS并联机构为稳定结构,结构示意图如图5。在理想状况下,六根加载支链上沿着加载支链的方向产生6个力,在动平台上合成为一个六维力矢。
当以力旋量表示6个力时,考虑动平台的平衡,6个杆的力螺旋之和应与平台的6维力相平衡,因此可列出螺旋方程。
f1$1+f2$2+…+fi$i+…=F+∈M    (1)
其中,i=1、2、3、4、5或6;
式中,fi为第i个杆受到的轴力;$i为第i杆轴线对固定坐标系的单位线矢,$i=Si+∈SOi,Si·Si=1,Si·SOi=0;F及M分别为平台上作用力的主矢和对坐标原点主矩。上述螺旋方程可以写为矩阵形式的平衡方程:
F = [ G f F ] f
其中F={Fx,Fy,Fz,Mx,My,Mz}T;Fx,Fy,Fz,Mx,My,Mz分别为沿着轴x,y,z三个方向的力和扭矩;f={f1,f2…f6}T
Figure BDA0000075468050000062
为静力影响系数矩阵。
[ G f F ] = S 1 S 2 . . . S 6 S 01 S 02 . . . S 06 - - - ( 2 )
若动平台6个球铰链分别记为b1,b2…b6,它们对O-XYZ坐标系的空间位置以矢量b1,b2…b6表示,静平台的6个虎克铰链以B1,B2……B6表示,而对O-XYZ坐标系的位置,以矢量B1,B2,…,B6表示,则:
S i = b i - B i | b i - B i | ; S Oi = B i × S i = B i × b i | b i - B i | 其中,i=1、2、3、4、5或6;
则影响系数矩阵为
[ G f F ] = b 1 - B 1 | b 1 - B 1 | b 2 - B 2 | b 2 - B 2 | . . . b 6 - B 6 | b 6 - B 6 | B 1 × b 1 | b 1 - B 1 | B 2 × b 2 | b 2 - B 2 | . . . B 6 × b 6 | b 6 - B 6 | - - - ( 3 )
根据6-UPS并联机构影响系数矩阵(力雅克比矩阵)可知,力雅克比矩阵是满秩矩阵,所以当动平台上产生合力及合力矩已知,在各加载支链的位置、速度、加速度即已经确定并且唯一。
由牛顿-欧拉法可得加载支链的动力学方程为
τ ri = m ru l · · ri + 3 4 m ru l ri α · i 2 + m ru g sin α i + F · l ri - - - ( 4 )
其中,i=1、2、3、4、5或6;
式中,τri为加载支链要求驱动力;mrui为加载支链前段(与动平台连接)质量;
Figure BDA0000075468050000068
为加载支链上端构件沿自身矢量方向上的移动加速度;lri为第i个加载支链的长度;αi为第i个加载支链与动平台平面的夹角,由几何关系可知,
Figure BDA0000075468050000069
g为单位重力矢;Fi为第i个球铰链所受外载荷力;lri为第i个加载支链矢量。根据拉格朗日法可将上式写成如下形式:
τ ri = H ri ( q i ) q · · i + C ri ( q i , q · i ) q · i + G ri ( q i ) - - - ( 5 )
其中,i=1、2、3、4、5或6;
式中,Hri(qi)为第i个加载支链的惯性矩阵;
Figure BDA0000075468050000072
为第i个加载支链阻尼项系数矩阵;Gri(qi)为第i个加载支链重力和动平台外载荷力相关的矩阵;
Figure BDA0000075468050000073
为第i个的末端加速度向量;
Figure BDA0000075468050000074
为速度向量;qi为机构的末端位姿,如图6为6-UPS并联机构加载支链受力图。
以上建立了6-UPS并联机构各个加载支链的控制模型,但该6-UPS并联机构在高频加载过程中,其控制模型参数随时间变化,因此需采用模糊PID控制器对6-UPS并联机构进行控制,可以解决疲劳加载试验中试验载荷的幅值漂移问题,从而降低了反馈系统的误差,同时也提高了试验精度。正如前文所述,该六维力高频疲劳试验机将使用模糊控制和传统PID控制相结合,设计一种模糊PID控制器,以误差e和误差变化
Figure BDA0000075468050000075
作为输入,可以满足不同时刻对e和
Figure BDA0000075468050000076
对控制参数自整定的要求。根据模糊控制理论对控制参数实行在线修改,使得被控对象品质指标始终保持在最佳范围之内。
设模糊PID控制器的输入输出集合为
R = { ( e , e · , u ) | ( e , e · , u ) ∈ { NB , NM , NS , ZO , PS , PM , PB } } - - - ( 6 )
式中,NB表示数值为负大,NM为负中,NS为负小,ZO为零,PS为正小,PM为正中,PB为正大,u(u∈R)为模糊PID控制器输出。设i∈I为误差e的个数,将式中各数值分别量化,则对于e∈[ei,ei+1],有
A i ( e ) = e i + 1 - e e i + 1 - e i A i + 1 ( e ) = e - e i e i + 1 - e i A k ( e ) = 0 ( k ≠ ( i , i + 1 ) ∈ I ) - - - ( 7 )
设j∈J为误差变化
Figure BDA0000075468050000082
的个数,则对于
Figure BDA0000075468050000083
B j ( e · ) = e · j + 1 - e · e · j + 1 - e · j B j + 1 ( e · ) = e · - e · j e · j + 1 - e · j B t ( e · ) = 0 ( t ≠ ( j , j + 1 ) ∈ J ) - - - ( 8 )
模糊推理函数为
f ij = A i ( e ) B j ( e · ) (i∈I,j∈J)(9)
式中Ai(e)和(i∈I,j∈J)分别表示e和的相应输入量模糊集合。
采用重心法对输出量模糊集进行去模糊化,则可得模糊推理机的实际输出为:
u = Σ i , j f ij u ij Σ i , j f ij = Σ k = ( i , i + 1 ) t = ( j , j + 1 ) ( A k ( e ) B t ( e · ) ) u kt Σ k = ( i , i + 1 ) t = ( j , j + 1 ) ( A k ( e ) B t ( e · ) ) - - - ( 10 )
式中,ukt∈{-3,-2,-1,0,1,2,3}为推理机模糊输出值的数值化表示,u∈R为模糊PID控制器输出,其他参数定义见公式(9)。
定义模糊PID控制器参数kP、kI和kD的调整算式如下:
k P = k P 0 + f { e , e · } k P = k P 0 + Δ k P k I = k I 0 + f { e , e · } k I = k I 0 + Δ k I k D = k D 0 + f { e , e · } k D = k D 0 + Δ k D - - - ( 11 )
其中,kP为模糊PID控制器的输出参数中的比例系数;kI为模糊PID控制器的输出参数中的积分系数;kD为模糊PID控制器的输出参数中的微分系数;
比例系数kP增加,则系统的响应速度快,但超调量增加;积分系数kI增加,则系统的稳态误差消除快,但过渡过程超调大;微分系数kD增加,则系统动态特性好,但抗干扰能力变差;
kP0、kI0和kD0为相应参数的初始值;ΔkP、ΔkI和ΔkD为在线控制过程中,由模糊PID控制器输出的参数调整量。在机构运动过程中,在线检测加载支链实际的驱动力和指令值之间的差值e和
Figure BDA0000075468050000091
由上式即可完成模糊PID控制器的参数调整。整个控制的流程图如图7。
b.反馈系统控制
该六维力高频疲劳试验机的反馈系统也应用力伺服控制。力伺服控制是指将压力传感器安装在6-UPS并联机构加载支链上,通过六个加载支链上的压力传感器得到力反馈信号,通过对六个加载支链的反馈信号的处理,进而对6-UPS并联机构的动平台进行力伺服控制,如图8。力伺服控制的优点在于单个力输出不存在相互干扰的问题,通过6-UPS并联机构力的反解很容易得到加载平台合力的结果,再通过模糊PID控制器就可以对当前信号进行对比、分析、控制。并且应用力伺服控制时,各个加载支链的互换性较好。
c.系统软件
该六维力高频疲劳试验机软件部分采用分层结构,按照功能分成4个部分,分别完成人机界面层、数据/参数管理层、通讯层及控制层(如图9)。
人机界面层主要完成运动模式选择、运动参数输入、运动启动/停止控制等操作任务,同时需要与数据库连接,将系统控制设置参数以及运动过程中的状态参数显示,以便对运动过程进行人工干预,因此它的关联部分为数据库。其具体任务描述为图10。
数据/参数管理层实际上就是一个数据库,需要管理的数据主要包括:各类系统控制参数、默认参数、运动状态等等(如图11)。
考虑到6-UPS并联结构控制的特殊性,拟将各轴运动规划也放在该部分内,便于计算形成各种参数的管理。
通讯层相对简单,主要任务就是根据具体的硬件进行时序协调管理。例如通讯接口、网口等。
在本软件系统中,控制层即对运动控制卡(PMAC)及筒式直线电机3驱动器的编程及调试。为了使软件系统达到良好的控制效果,要保证运动的稳定性、准确性及快速性,控制层的主要工作包括筒式直线电机3驱动器与运动控制卡连接,系统反馈控制与运动控制卡连接以及6-UPS并联机构本身位置与运动控制卡的连接,如图12。
为了增强系统软件的实时性和系统软件适应高频加载试验的需要,利用多线程技术对软件进行编写。因为系统软件运行时循环模块不停地采集数据、显示数据、存储数据、调整、判断保护,无限循环下去直到程序退出为止,因此软件会耗费大量的运算时间,同时由于控制系统资源被大量占用,采集到的数据的同步性并不好,例如应力到峰值时应变却不是峰值,这就造成实验数据的精度不够高。利用多线程技术能很好地解决这些问题,使程序运行得更快并且更稳定。
利用多线程来优化系统软件,就是把程序的循环模块单独放在一个线程中来执行,把它与主线程隔离开,让主线程来完成循环模块以外的所有其他功能。循环模块所在的线程是一个高速运行的线程,依靠循环完成不停采集数据和控制等功能。而窗体的主线程除了负责循环模块以外的其他功能外,还负责响应用户输入,例如鼠标的点击、键盘的输入等。这样可以很好的提高系统软件的实时性。
本发明一种六维力高频疲劳试验机的使用方法,其步骤如下:
①分析试样在实际工作中的受力状况和建立试样的时变载荷谱;
分析确定试样在工作中的受力状况,包括集中载荷作用点及方向,分布载荷的分布规律,试样的约束情况,以及上述两类载荷状况随时间的变化规律,根据载荷的空间、时间分布规律,建立试样的载荷谱:F=F(x,y,z,t),即载荷F与载荷作用点x,y,z和时间t之间的函数关系;
②求出试样应力应变分布及其时间函数;
依据载荷谱采用有限元方法计算出结构件的应力、应变在试样中的分布,求得这种分布状况随时间的变化规律:σ=σ(x,y,z,t),ε=ε(x,y,z,t),即应力σ和应变ε的空间位置x,y,z分布与时间t之间的函数关系。
③采用有限元仿真技术确定这种应力状态下的加载点位移及其时间函数;
根据上一步求出的应力应变分布及其时间函数,在有限元仿真环境下,采用单点集中力给试样加载(可以是多维力载荷),产生与工作状态尽量相似的应力应变分布及其随时间的变化,确定加载点位置、载荷大小及方向、产生的位移大小及方向,以及位移的时间函数。
④按照加载位移的时间函数进行多维力加载,根据加载位移的时间函数和模糊PID控制,确定加载平台运动要求,由工控机根据加载平台的运动求解6-UPS并联机构反解,得出每根并联加载支链的运动形式,筒式直线电机3系统产生加载运动。
为求运动反解,在加载平台建立坐标系(O′-X′Y′Z′),在静平台建立坐标系(O-XYZ),每一根杆用一个矢量Li表示(i=1,2,...6.).,如图5。
图中,H是加载平台坐标原点在静平台坐标系中的位置矢量,bi是加载平台铰链点位置矢量,Bi是静平台铰链点位置矢。由图示的矢量关系可得:
Bi=H+biRT-Li    (12)
式中,R是加载平台座标系相对于静平台坐标系的姿态矩阵:
R = c θ 2 c θ 3 s θ 1 s θ 2 c θ 3 - c θ 1 s θ 3 c θ 1 s θ 2 c θ 3 + s θ 1 s θ 3 c θ 2 s θ 3 s θ 1 s θ 2 s θ 3 + c θ 1 c θ 3 c θ 1 s θ 2 s θ 3 - s θ 1 c θ 3 - s θ 2 sθ 1 c θ 2 c θ 1 c θ 2 - - - ( 13 )
矩阵中sin和cos简记为s和c,θ1、θ2、θ3分别表示加载平台绕x、y、z轴的转角。
根据加载平台加载时所需的位置H和姿态θ1、θ2、θ3,由式12和13,即可算出静平台铰链点Bi位置,实现对6-UPS并联机构的控制。
控制系统根据所算出的一系列杆长值,控制各杆产生相应位移,使加载平台实现加载运动,然后通过拉压力传感器2的反馈,应用力伺服控制,对加载运动所产生的加载力和力矩与试样在工作中的受力状况进行对比,经过误差计算最后确定加载支链杆长的变化值,使得加载力的大小和方向与试样在工作中的受力状况相同。
⑤分析试样的实际应力分布,作为结构设计的参考;
进行多维力加载的同时,采用粘贴应变片等措施在试样关键位置监测应力应变,采样系统采集应力应变值,与载荷谱的理论值对比;分析试样的实际应力分布,作为试样结构设计的参考。
对于小型动载荷机械零件,一个并联试验机构即可实现对结构件的多维力加载(如图13)。
本发明一种六维力高频疲劳试验机及使用方法,其优点和积极效果是:其采用并联试验机构对被测试样施加三维力和三维力矩,在尽可能符合实际受理环境下对试样施加多维变频力或者疲劳测试,得到更加准确的应力应变分布数据、更加准确的找到危险位置和更加真实的破坏形式,从而更真实的模拟试样在工作状态下的受力状况,更准确的测定试样的承载能力。为各种飞机等精密武器结构的轻量化、优化、提高可靠性,具有重要的意义。
(四)附图说明
图1为本发明六维力高频疲劳试验机的结构框图
图2为本发明6-UPS六维力高频疲劳试验机构指示图
图3为本发明控制系统组成示意图
图4为本发明应用模糊逻辑控制器流程图
图5为本发明6-UPS结构示意图
图6为本发明加载支链受力图
图7为本发明控制系统流程图
图8为本发明系统反馈的力伺服控制流程图
图9为本发明软件分层结构
图10为本发明界面层任务
图11为本发明管理层任务
图12为本发明控制层任务
图13为本发明对动载荷机械部件进行疲劳加载示意图
图中的符号说明如下:
图5中bi为加载平台铰链点;Bi为静平台铰链点;xyz为静坐标系;x′y′z′为动坐标系;
图6中Ai和Bi分别为加载平台和静平台铰链点矢量;τr为加载支链要求驱动力;mrui为加载支链前段(与加载平台连接)质量;α为加载支链与加载平台平面的夹角,由几何关系可知,
Figure BDA0000075468050000121
g为单位重力矢;Fi为球铰链所受外载荷力;
图7中F为加载平台要求驱动力;
Figure BDA0000075468050000122
为力雅克比矩阵;fi为各加载支链要求分力;kP、kI和kD为模糊PID控制器的输出参数;kPO、kIO和kDO为相应参数的初始值;ΔkP、ΔkI和ΔkD为在线控制过程中,由模糊运算器输出的参数调整量;
Hri(qi)为第i个加载支链的惯性矩阵;
为第i个加载支链阻尼项系数矩阵;
Gr1(qi)为第i个加载支链重力和加载平台外载荷力相关的矩阵;
为第i个的末端加速度向量;
Figure BDA0000075468050000133
为速度向量;qi为机构的末端位姿;
F*为加载平台最终输出力;
图8中F为加载平台要求驱动力;Fe为反解后单加载支链驱动力;
Fg为电机单杆实际输出力;F*为加载平台最终输出力。
(五)具体实施方式
本发明提出一种以动载荷机械零件的疲劳试样为测试对象的六维力高频疲劳试验机及使用方法,可用于模拟动载荷机械零件的实际受力情况,测定其在静载荷、动载荷及交变疲劳试验中的变形情况,以解决现有疲劳试验机不能施加高频多维力的不足,以接近于实际工作状态下的受力状况,准确的测定试件的疲劳承载能力。
本发明一种六维力高频疲劳试验机,它包括6-UPS并联机构和控制系统两个部分(如图1),其特征在于:
1.6-UPS并联机构
该六维力高频疲劳试验机使用6-UPS并联机构,其是6-SPS并联机构的变形。6-SPS并联机构由6个运动链连接动平台和静平台组成,并以移动副作为输入,通常称为Stewart平台。六维力高频疲劳试验机使用6-SPS的变形机构6-UPS机构,其与6-SPS机构相同,由动平台,静平台及加载支链三部分组成,不同是6-SPS并联机构的加载支链两端由球铰链分别与动平台和静平台连接。而6-UPS并联机构的加载支链两端分别由球铰链和虎克铰链与动平台和静平台连接。6-SPS并联机构的每一根杆要有六自由度,因为与动平台相连的关节轴承提供了三个转动自由度,移动副提供了一个移动自由度,故静平台一侧的关节轴承只需提供两个转动自由度即可。因此静平台采用了两个转动自由度的虎克铰链。这样即形成了6-UPS并联机构。
该6-UPS并联机构由加载平台1、拉压力传感器2、筒式直线电机3、虎克铰链4、静平台5及球铰链6组成,如图2;
其中,所述加载平台1,即6-UPS并联机构的动平台,其由六个球铰链6支撑。高刚度的加载平台在安装夹具后可以对动载荷机械零件件进行加载。
由球铰链6、拉压力传感器2及筒式直线电机3和虎克铰链4组成加载支链,通过筒式直线电机3的伸缩运动来完成加载支链的长度变化,通过筒式直线电机3与球铰链6和虎克铰链4共同作用下改变加载平台1的位姿和加载力的大小。所述筒式直线电机3作为试验机的驱动部件,也是移动副。筒式直线电机3具有高频,高加速度,低惯量的优点,是疲劳加载试验优秀性能的前提。
所述静平台5作为所述6-UPS并联机构的基座,通过虎克铰链4与加载支链连接,具有高强度,高刚度的性能。
所述静平台5和加载平台1以6个加载支链相联,每个加载支链两端分别是球铰链6和虎克铰链4,中间是移动副(筒式直线电机3)。筒式直线电机3在驱动器作用下作相对移动,改变加载支链的长度,使加载平台1的位置和姿态发生变化。
该六维力高频疲劳试验机使用高性能筒式直线电机3作为移动副,使6-UPS并联机构实现位置变化及加载运动。通过六根加载支链的运动组合使得加载平台(动平台)产生六自由度的加载运动,可以对被测试样施加六维载荷(三维力和三维力矩)。六根加载支链的运动由计算机控制,可以根据所需要载荷谱产生给定的运动组合及其变化规律。
动载荷机械零件的疲劳试样通过夹具接受载荷,可以实现三维力和三维力矩的高频力学性能测试。动平台的移动由六个加载支链带动实现。加载支链部件在自身伺服主机的驱动下,实现精确的轴向伸缩。六根加载支链的伸缩组合,实现对上夹具的各种运动驱动,进而对动载荷机械零件的疲劳试样施加各种单个或组合载荷。
该六维力高频疲劳试验机的工作原理是:每根加载支链通过筒式直线电机3产生伸缩运动,六根加载支链的运动组合使得加载平台1产生六自由度的加载运动,可以对被测试件施加六维载荷(三维力和三维力矩)。六根加载支链的运动由计算机控制,可以根据所需的载荷谱产生给定的运动组合及其变化规律。
2.控制系统
该控制系统由系统硬件、系统控制方法及软件两部分组成;
(1).系统硬件
该系统硬件包括工控机(PC)、运动控制卡、筒式直线电机3、拉压力传感器2构成,如图3。工控机和运动控制卡组成六维力高频疲劳试验机控制部分,筒式直线电机3和拉压力传感器2安装在6-UPS并联机构上,如图2。
该工控机为全钢机箱(19″,4U),全长主板、无源底板,工业电源(2.0GPentium IV CPU,2G内存)。工控机作为控制卡的载体,实现与控制卡的通讯,完成系统参数管理、用户界面的操作、非实时性任务的实现(例如运动规划和运算等)等功能。
该运动控制卡采用6轴以上可编程多轴运动卡。该运动控制卡可以同步控制8个驱动轴,实现复杂的多轴协调运动。它能够对存储在它内部的程序进行单独的运算,执行运动程序、PLC程序、进行伺服环更新,及以串口、总线两种方式与主计算机进行通讯,而且它还可以自动对任务进行优先等级判别,从而进行实时的多任务处理,这使得它在处理时间和任务切换这两方面大大减轻主机和编程器的负担,提高了整个控制系统的运行速度和控制精度。其单轴最快的采样周期为60μs(微秒)。
该筒式直线电机3采用大推力高加速度筒式直线电机。该筒式直线电机3安装在6-UPS并联机构的加载支链上,一端连接拉压力传感器2,另外一端连接虎克铰链,为六维力高频疲劳试验机提供动力。该筒式直线电机3是三相无轴承水冷同步伺服直线电机。该筒式直线电机3是三相无轴承水冷同步伺服直线电机,其最大推力不小于850N(牛),持续推力不小于550N(牛),最大功率不小于1400W(瓦),电机空载最大加速度不小于350m/s2(米/秒2)。
该拉压力传感器2采用柱式高精度拉压力传感器。该拉压力传感器2一端与球铰链6连接,另外一端与筒式直线电机3连接。该拉压力传感器2的弹性体为柱式结构,用于拉伸力和压缩力测量。全密封结构,其输出灵敏度为1.5mV/V(毫伏/伏),直线度为±0.05%F.S(满量程),重复性为±0.05%F.S(满量程),最大量程为±5000N(牛)。
控制系统的核心为运动控制卡,工业控制计算机负责信息流和数据流的管理,以及从传感器、编码器读取力和位置数据,并经过计算后发送控制指令。驱动器负责直线电机的功率驱动,实现位置、速度和力的控制。
(2).系统控制方法及软件
系统控制方法及软件由系统控制、反馈系统控制及系统软件三部分组成。
6-UPS并联机构的控制包括运动学和动力学控制。运动学控制指的是位置控制,动力学控制是指力控制或力和位置的混合控制。对于运动学方法,我们需要知道机构的运动学正解和运动学逆解。当控制机构运动时,首先根据末端执行器的位置通过运动学逆解计算出各关节的位置,然后控制关节到达指定位置。动力学控制是指综合考虑机构的运动学和动力学模型生成控制指令。由于通过简单的运动学逆解运算即可生成运动指令,运动学控制简单、可靠易实现,因此应用广泛。并联机床的运动呈现高度的非线性和强耦合,使得并联机床的运动尤其在高速运动时会对运动控制精度产生较大影响,因此需要建立动力学控制系统来对并联机床进行控制。
a.系统控制
该六维力高频疲劳试验机的控制使用力伺服控制。该力伺服控制通过模糊PID控制器来实现。目前机器人的力伺服控制主要有以下几类:阻抗控制,力/位混合控制,自适应控制和智能控制。模糊控制是智能控制中的一种,模糊控制的优点在于对非线性系统有良好的控制作用。图4即为模糊逻辑控制器流程图,其是模糊控制的核心,它是根据各种边界条件和一些重要的专家规则确定的。
在工程实际中,应用最为广泛的调节器控制规律为比例、积分、微分控制,简称PID控制。它以其结构简单、稳定性好、工作可靠、调整方便而成为工业控制的主要技术之一。当被控对象的结构和参数不能完全掌握,或得不到精确的数学模型时,控制理论的其它技术难以采用时,系统控制器的结构和参数必须依靠经验和现场调试来确定,这时应用PID控制技术最为方便。而PID控制的不足是当系统存在非线性,噪声及强耦合时无法解决稳定性及准确性之间的矛盾,并且PID控制较难取得令人满意的动态响应。
模糊PID控制器结合传统PID控制和模糊控制的优点,稳定性强,在线控制方便且能适应模型参数的变化,因此该六维力高频疲劳试验机的控制——力伺服控制采用模糊PID控制器来实现。模糊PID控制器对各个加载支链的驱动力进行实时跟踪,从而实现对6-UPS并联机构的力伺服控制。
具体地说,首先基于6-UPS并联机构各加载支链的动力学方程建立其控制模型,然后采用模糊PID控制器对加载支链驱动力进行跟踪,从而保证6-UPS并联机构在整个运动过程中具有稳定而准确的驱动力输出。
要进行力伺服控制,所以先介绍一下6-UPS并联机构的受力分析。在静力分析时,可以应用螺旋理论建立起6-UPS并联机构的静力平衡方程,并转化为用影响系数表示的矩阵方程。对于该试验机,机构的位形是一定的,六根杆长l1、l2、l3、l4、l5、l6是确定的,它的动平台经由6个加载支链与静平台相连,每杆两端分别为球铰链和虎克铰链,此6-UPS并联机构为稳定结构,结构示意图如图5。在理想状况下,六根加载支链上沿着加载支链的方向产生6个力,在动平台上合成为一个六维力矢。
当以力旋量表示6个力时,考虑动平台的平衡,6个杆的力螺旋之和应与平台的6维力相平衡,因此可列出螺旋方程。
f1$1+f2$2+…+fi$i+…=F+∈M  (i=1...6)(1)
式中,fi为第i个杆受到的轴力;$i为第i杆轴线对固定坐标系的单位线矢,$i=Si+∈S0i,Si·Si=1,Si·S0i=0;F及M分别为平台上作用力的主矢和对坐标原点主矩。上述螺旋方程可以写为矩阵形式的平衡方程:
F = [ G f F ] f
其中F={Fx,Fy,Fz,Mx,My,Mz}T;Fx,Fy,Fz,Mx,My,Mz分别为沿着轴x,y,z三个方向的力和扭矩;f={f1,f2…f6}T
Figure BDA0000075468050000172
为静力影响系数矩阵。
[ G f F ] = S 1 S 2 . . . S 6 S 01 S 02 . . . S 06 - - - ( 2 )
若动平台6个球铰链分别记为b1,b2…b6,它们对O-XYZ坐标系的空间位置以矢量b1,b2…b6表示,静平台的6个虎克铰链以B1,B2……B6表示,而对O-XYZ坐标系的位置,以矢量B1,B2,…,B6表示,则:
S i = b i - B i | b i - B i | ; S 0 i = B i × S i = B i × b i | b i - B i | (i=1...6)
则影响系数矩阵为
[ G f F ] = b 1 - B 1 | b 1 - B 1 | b 2 - B 2 | b 2 - B 2 | . . . b 6 - B 6 | b 6 - B 6 | B 1 × b 1 | b 1 - B 1 | B 2 × b 2 | b 2 - B 2 | . . . B 6 × b 6 | b 6 - B 6 | - - - ( 3 )
根据6-UPS并联机构影响系数矩阵(力雅克比矩阵)可知,力雅克比矩阵是满秩矩阵,所以当动平台上产生合力及合力矩已知,在各加载支链的位置、速度、加速度即已经确定并且唯一。
由牛顿-欧拉法可得加载支链的动力学方程为
τ ri = m ru l · · ri + 3 4 m ru l ri α · i 2 + m ru g sin α i + F · l ri (i=1...6)(4)
式中,τri为加载支链要求驱动力;mrui为加载支链前段(与动平台连接)质里;
Figure BDA0000075468050000182
为加载支链上端构件沿自身矢量方向上的移动加速度;lri为第i个加载支链的长度;αi为第i个加载支链与动平台平面的夹角,由几何关系可知,g为单位重力矢;Fi为第i个球铰链所受外载荷力;lr1为第i个加载支链矢量。根据拉格朗日法可将上式写成如下形式:
τ ri = H ri ( q i ) q · · i + C ri ( q i , q · i ) q · i + G ri ( q i ) (i=1...6)(5)
式中,Hr1(qi)为第i个加载支链的惯性矩阵;
Figure BDA0000075468050000185
为第i个驱动支链阻尼项系数矩阵;Gri(qi)为第i个加载支链重力和动平台外载荷力相关的矩阵;
Figure BDA0000075468050000186
为第i个的末端加速度向量;
Figure BDA0000075468050000187
为速度向量;qi为机构的末端位姿,如图6为6-UPS并联机构加载支链受力图。
以上建立了6-UPS并联机构各个加载支链的控制模型,但该6-UPS并联机构在高频加载过程中,其控制模型参数随时间变化,因此需采用模糊PID控制器对6-UPS并联机构进行控制,可以解决疲劳加载试验中试验载荷的幅值漂移问题,从而降低了反馈系统的误差,同时也提高了试验精度。正如前文所述,该六维力高频疲劳试验机将使用模糊控制和传统PID控制相结合,设计一种模糊PID控制器,以误差e和误差变化
Figure BDA0000075468050000188
作为输入,可以满足不同时刻对e和
Figure BDA0000075468050000189
对控制参数自整定的要求。利用模糊控制对控制参数实行在线修改,使得被控对象品质指标始终保持在最佳范围之内。
设模糊PID控制器的输入输出集合为
R = { ( e , e · , u ) | ( e , e · , u ) ∈ { NB , NM , NS , ZO , PS , PM , PB } } - - - ( 6 )
式中,NB表示数值为负大,NM为负中,NS为负小,ZO为零,PS为正小,PM为正中,PB为正大,u(u∈R)为模糊PID控制器输出。设i∈I为误差e的个数,将式中各数值分别量化,则对于e∈[ei,ei+1],有
A i ( e ) = e i + 1 - e e i + 1 - e i A i + 1 ( e ) = e - e i e i + 1 - e i A k ( e ) = 0 ( k ≠ ( i , i + 1 ) ∈ I ) - - - ( 7 )
设j∈J为误差变化的个数,则对于
Figure BDA0000075468050000194
B j ( e · ) = e · j + 1 - e · e · j + 1 - e · j B j + 1 ( e · ) = e · - e · j e · j + 1 - e · j B t ( e · ) = 0 ( t ≠ ( j , j + 1 ) ∈ J ) - - - ( 8 )
模糊推理函数为
f ij = A i ( e ) B j ( e · ) (i∈I,j ∈J)(9)
式中Ai(e)和(i∈I,j∈J)分别表示e和
Figure BDA0000075468050000198
的相应输入量模糊集合。
采用重心法对输出量模糊集进行去模糊化,则可得模糊推理机的实际输出为:
u = Σ i , j f ij u ij Σ i , j f ij = Σ k = ( i , i + 1 ) t = ( j , j + 1 ) ( A k ( e ) B t ( e · ) ) u kt Σ k = ( i , i + 1 ) t = ( j , j + 1 ) ( A k ( e ) B t ( e · ) ) - - - ( 10 )
式中,ukt∈{-3,-2,-1,0,1,2,3}为推理机模糊输出值的数值化表示,u∈R为模糊PID控制器输出。
众所周知,模糊PID控制器的参数变化影响被控系统的过渡特性,如:比例系数kP增加,则系统的响应速度快,但超调量增加;积分系数kI增加,则系统的稳态误差消除快,但过渡过程超调大;微分系数kD增加,则系统动态特性好,但抗干扰能力变差。
定义模糊PID控制器参数kP、kI和kD的调整算式如下:
k P = k P 0 + f { e , e · } k P = k P 0 + Δ k P k I = k I 0 + f { e , e · } k I = k I 0 + Δ k I k D = k D 0 + f { e , e · } k D = k D 0 + Δ k D - - - ( 11 )
其中,kP、kI和kD为模糊PID控制器的输出参数;kP0、kI0和kD0为相应参数的初始值;ΔkP、ΔkI和ΔkD为在线控制过程中,由模糊PID控制器输出的参数调整量。在机构运动过程中,在线检测加载支链实际的驱动力和指令值之间的差值e和
Figure BDA0000075468050000202
由上式即可完成模糊PID控制器的参数调整。整个控制的流程图如图7。
b.反馈系统控制
该六维力高频疲劳试验机的反馈系统也应用力伺服控制。力伺服控制是指将压力传感器安装在6-UPS并联机构加载支链上,通过六个加载支链上的压力传感器得到力反馈信号,通过对六个加载支链的反馈信号的处理,进而对6-UPS并联机构的动平台进行力伺服控制,如图8。力伺服控制的优点在于单个力输出不存在相互干扰的问题,通过6-UPS并联机构力的反解很容易得到加载平台合力的结果,再通过模糊PID控制器就可以对当前信号进行对比、分析、控制。并且应用力伺服控制时,各个加载支链的互换性较好。
c.系统软件
该六维力高频疲劳试验机软件部分采用分层结构,按照功能分成4个部分,分别完成人机界面层、数据/参数管理层、通讯层及控制层(如图9)。
人机界面层主要完成运动模式选择、运动参数输入、运动启动/停止控制等操作任务,同时需要与数据库连接,将系统控制设置参数以及运动过程中的状态参数显示,以便对运动过程进行人工干预,因此它的关联部分为数据库。其具体任务描述为图10。
数据/参数管理层实际上就是一个数据库,需要管理的数据主要包括:各类系统控制参数、默认参数、运动状态等等(如图11)。
考虑到6-UPS并联结构控制的特殊性,拟将各轴运动规划也放在该部分内,便于计算形成各种参数的管理。
通讯层相对简单,主要任务就是根据具体的硬件进行时序协调管理。例如通讯接口、网口等。
在本软件系统中,控制层即对运动控制卡(PMAC)及筒式直线电机3驱动器的编程及调试。为了使软件系统达到良好的控制效果,要保证运动的稳定性、准确性及快速性,控制层的主要工作包括筒式直线电机3驱动器与运动控制卡连接,系统反馈控制与运动控制卡连接以及6-UPS并联机构本身位置与运动控制卡的连接,如图12。
为了增强系统软件的实时性和系统软件适应高频加载试验的需要,利用多线程技术对软件进行编写。因为系统软件运行时循环模块不停地采集数据、显示数据、存储数据、调整、判断保护,无限循环下去直到程序退出为止,因此软件会耗费大量的运算时间,同时由于控制系统资源被大量占用,采集到的数据的同步性并不好,例如应力到峰值时应变却不是峰值,这就造成实验数据的精度不够高。利用多线程技术能很好地解决这些问题,使程序运行得更快并且更稳定。
利用多线程来优化系统软件,就是把程序的循环模块单独放在一个线程中来执行,把它与主线程隔离开,让主线程来完成循环模块以外的所有其他功能。循环模块所在的线程是一个高速运行的线程,依靠循环完成不停采集数据和控制等功能。而窗体的主线程除了负责循环模块以外的其他功能外,还负责响应用户输入,例如鼠标的点击、键盘的输入等。这样可以很好的提高系统软件的实时性。
本发明一种六维力高频疲劳试验机的使用方法,其步骤如下:
①分析试样在实际工作中的受力状况和建立试样的时变载荷谱;
分析确定试样在工作中的受力状况,包括集中载荷作用点及方向,分布载荷的分布规律,试样的约束情况,以及上述两类载荷状况随时间的变化规律,根据载荷的空间、时间分布规律,建立试样的载荷谱:F=F(x,y,z,t),即载荷F与载荷作用点x,y,z和时间t之间的函数关系;
②求出试样应力应变分布及其时间函数;
依据载荷谱采用有限元方法计算出结构件的应力、应变在试样中的分布,求得这种分布状况随时间的变化规律:σ=σ(x,y,z,t),ε=ε(x,y,z,t),即应力σ和应变ε的空间位置x,y,z分布与时间t之间的函数关系。
③采用有限元仿真技术确定这种应力状态下的加载点位移及其时间函数;
根据上一步求出的应力应变分布及其时间函数,在有限元仿真环境下,采用单点集中力给试样加载(可以是多维力载荷),产生与工作状态尽量相似的应力应变分布及其随时间的变化,确定加载点位置、载荷大小及方向、产生的位移大小及方向,以及位移的时间函数。
④按照加载位移的时间函数进行多维力加载,根据加载位移的时间函数和模糊PID控制,确定加载平台运动要求,由工控机根据加载平台的运动求解6-UPS并联机构反解,得出每根并联加载支链的运动形式,筒式直线电机3系统产生加载运动。
为求运动反解,在加载平台建立坐标系(O′-X′Y′Z′),在静平台建立坐标系(O-XYZ),每一根杆用一个矢量Li表示(i=1,2,...6.).,如图5。
图中,H是加载平台坐标原点在静平台坐标系中的位置矢量,bi是加载平台铰链点位置矢量,Bi是静平台铰链点位置矢。由图示的矢量关系可得:
Bi=H+biRT-Li    (12)
式中,R是加载平台座标系相对于静平台坐标系的姿态矩阵:
R = c θ 2 c θ 3 s θ 1 s θ 2 c θ 3 - c θ 1 s θ 3 c θ 1 s θ 2 c θ 3 + s θ 1 s θ 3 c θ 2 s θ 3 s θ 1 s θ 2 s θ 3 + c θ 1 c θ 3 c θ 1 s θ 2 s θ 3 - s θ 1 c θ 3 - s θ 2 sθ 1 c θ 2 c θ 1 c θ 2 - - - ( 13 )
矩阵中sin和cos简记为s和c,θ1、θ2、θ3分别表示加载平台绕x、y、z轴的转角。
根据加载平台加载时所需的位置H和姿态θ1、θ2、θ3,由式12和13,即可算出静平台铰链点Bi位置,实现对6-UPS并联机构的控制。
控制系统根据所算出的一系列杆长值,控制各杆产生相应位移,使加载平台实现加载运动,然后通过拉压力传感器2的反馈,应用力伺服控制,对加载运动所产生的加载力和力矩与试样在工作中的受力状况进行对比,经过误差计算最后确定加载支链杆长的变化值,使得加载力的大小和方向与试样在工作中的受力状况相同。
⑤分析试样的实际应力分布,作为结构设计的参考;
进行多维力加载的同时,采用粘贴应变片等措施在试样关键位置监测应力应变,采样系统采集应力应变值,与载荷谱的理论值对比;分析试样的实际应力分布,作为试样结构设计的参考。
对于小型动载荷机械零件,一个并联试验机构即可实现对结构件的多维力加载(如图13)。

Claims (9)

1.一种六维力高频疲劳试验机,它包括6-UPS并联机构和控制系统两个部分,其特征在于:
6-UPS并联机构:
该6-UPS并联机构由加载平台、拉压力传感器、筒式直线电机、虎克铰链、静平台及球铰链组成;
所述加载平台,即6-UPS并联机构的动平台,由六个球铰链支撑;加载平台在安装夹具后对动载荷机械零件件进行加载;
所述球铰链、拉压力传感器、筒式直线电机及虎克铰链组成加载支链;通过筒式直线电机的伸缩运动来完成加载支链的长度变化,通过筒式直线电机与球铰链和虎克铰链共同作用下改变加载平台的位姿和加载力的大小;所述筒式直线电机作为该六维力高频疲劳试验机的驱动部件,也是移动副;
所述静平台作为所述6-UPS并联机构的基座,通过虎克铰链与加载支链;
所述静平台和加载平台以6个加载支链相联,每个加载支链两端分别是球铰链和虎克铰链,中间是该筒式直线电机;该筒式直线电机在驱动器作用下作相对移动,改变加载支链的长度,使加载平台的位置和姿态发生变化;
所述六维力高频疲劳试验机使用高性能筒式直线电机作为移动副,使6-UPS并联机构实现位置变化及加载运动;通过6个加载支链的运动组合使得加载平台产生六自由度的加载运动,对被测试样施加六维载荷;6个加载支链的运动由计算机控制,根据所需要载荷谱产生给定的运动组合及其变化规律;
控制系统:
该控制系统由系统硬件、系统控制方法及软件两部分组成;
(1).系统硬件
该系统硬件包括工控机、运动控制卡、筒式直线电机及拉压力传感器;该工控机和运动控制卡组成该六维力高频疲劳试验机控制部分,筒式直线电机和拉压力传感器安装在6-UPS并联机构上;
控制系统的核心为运动控制卡,工业控制计算机负责信息流和数据流的管理,以及从传感器、编码器读取力和位置数据,并经过计算后发送控制指令;驱动器负责直线电机的功率驱动,实现位置、速度和力的控制;
(2).系统控制方法及软件
系统控制方法及软件由系统控制、反馈系统控制及系统软件三部分组成;
a.系统控制
该六维力高频疲劳试验机的控制使用力伺服控制;该力伺服控制通过模糊PID控制器来实现;模糊PID控制器对各个加载支链的驱动力进行实时跟踪,从而实现对6-UPS并联机构的力伺服控制;
具体地说,首先基于6-UPS并联机构各加载支链的动力学方程建立其控制模型,然后采用模糊PID控制器对加载支链驱动力进行跟踪,从而保证6-UPS并联机构在整个运动过程中具有稳定而准确的驱动力输出;
要进行力伺服控制,先要进行6-UPS并联机构的受力分析;在静力分析时,应用螺旋理论建立起6-UPS并联机构的静力平衡方程,并转化为用影响系数表示的矩阵方程;对于该试验机,机构的位形是一定的,六根杆长l1、l2、l3、l4、l5、l6是确定的,它的动平台经由6个加载支链与静平台相连,每杆两端分别为球铰链和虎克铰链,此6-UPS并联机构为稳定结构,6个加载支链上沿着加载支链的方向产生6个力,在动平台上合成为一个六维力矢;
当以力旋量表示6个力时,考虑动平台的平衡,6个杆的力螺旋之和应与平台的六维力相平衡,因此列出螺旋方程;
f1$1+f2$2+…+fi$i+…=F+∈M    (1)
其中,i=1、2、3、4、5或6;
式中,fi为第i个杆受到的轴力;$i为第i杆轴线对固定坐标系的单位线矢,$i=Si+∈S0i,Si·Si=1,Si·S0i=0;F及M分别为平台上作用力的主矢和对坐标原点主矩;上述螺旋方程写为矩阵形式的平衡方程:
F = [ G f F ] f
其中F={Fx,Fy,Fz,Mx,My,Mz}T;Fx,Fy,Fz,Mx,My,Mz分别为沿着轴x,y,z三个方向的力和扭矩;f={f1,f2…f6}T
Figure FDA0000075468040000022
为静力影响系数矩阵;
[ G f F ] = S 1 S 2 . . . S 6 S 01 S 02 . . . S 06 - - - ( 2 )
若动平台6个球铰链分别记为b1,b2…b6,它们对O-XYZ坐标系的空间位置以矢量b1,b2…b6表示,静平台的6个虎克铰链以B1,B2……B6表示,而对O-XYZ坐标系的位置,以矢量B1,B2,…,B6表示,则:
S i = b i - B i | b i - B i | ; S 0 i = B i × S i = B i × b i | b i - B i | 其中,i=1、2、3、4、5或6;
则影响系数矩阵为
[ G f F ] = b 1 - B 1 | b 1 - B 1 | b 2 - B 2 | b 2 - B 2 | . . . b 6 - B 6 | b 6 - B 6 | B 1 × b 1 | b 1 - B 1 | B 2 × b 2 | b 2 - B 2 | . . . B 6 × b 6 | b 6 - B 6 | - - - ( 3 )
根据6-UPS并联机构影响系数矩阵可知,力雅克比矩阵是满秩矩阵,所以当动平台上产生合力及合力矩已知,在各加载支链的位置、速度及加速度即已经确定并且唯一;
由牛顿-欧拉法得到加载支链的动力学方程为
τ ri = m ru l · · ri + 3 4 m ru l ri α · i 2 + m ru g sin α i + F · l ri - - - ( 4 )
其中,i=1、2、3、4、5或6;
式中,τri为加载支链要求驱动力;mrui为加载支链前段质量;
Figure FDA0000075468040000035
为加载支链上端构件沿自身矢量方向上的移动加速度;lri为第i个加载支链的长度;αi为第i个加载支链与动平台平面的夹角,由几何关系可知,
Figure FDA0000075468040000036
g为单位重力矢;Fi为第i个球铰链所受外载荷力;lri为第i个加载支链矢量;根据拉格朗日法将上式写成如下形式:
τ ri = H ri ( q i ) q · · i + C ri ( q i , q · i ) q · i + G ri ( q i ) - - - ( 5 )
其中,i=1、2、3、4、5或6;
式中,Hr1(qi)为第i个加载支链的惯性矩阵;
Figure FDA0000075468040000038
为第i个加载支链阻尼项系数矩阵;Gri(qi)为第i个加载支链重力和动平台外载荷力相关的矩阵;
Figure FDA0000075468040000039
为第i个的末端加速度向量;为速度向量;qi为机构的末端位姿;
以上建立了6-UPS并联机构各个加载支链的控制模型,但该6-UPS并联机构在高频加载过程中,其控制模型参数随时间变化,因此需采用模糊PID控制器对6-UPS并联机构进行控制;以误差e和误差变化
Figure FDA0000075468040000041
作为输入,满足不同时刻对e和
Figure FDA0000075468040000042
对控制参数自整定的要求;根据模糊控制理论对控制参数实行在线修改,使得被控对象品质指标始终保持在最佳范围之内;
设模糊PID控制器的输入输出集合为
R = { ( e , e · , u ) | ( e , e · , u ) ∈ { NB , NM , NS , ZO , PS , PM , PB } } - - - ( 6 )
式中,NB表示数值为负大,NM为负中,NS为负小,ZO为零,PS为正小,PM为正中,PB为正大,u(u∈R)为模糊PID控制器输出;设i∈I为误差e的个数,将式中各数值分别量化,则对于e∈[ei,ei+1],有
A i ( e ) = e i + 1 - e e i + 1 - e i A i + 1 ( e ) = e - e i e i + 1 - e i A k ( e ) = 0 ( k ≠ ( i , i + 1 ) ∈ I ) - - - ( 7 )
设j∈J为误差变化
Figure FDA0000075468040000045
的个数,则对于
Figure FDA0000075468040000046
B j ( e · ) = e · j + 1 - e · e · j + 1 - e · j B j + 1 ( e · ) = e · - e · j e · j + 1 - e · j B t ( e · ) = 0 ( t ≠ ( j , j + 1 ) ∈ J ) - - - ( 8 )
模糊推理函数为
f ij = A i ( e ) B j ( e · ) (i∈I,j∈J)(9)
式中Ai(e)和(i∈I,j∈J)分别表示e和
Figure FDA00000754680400000410
的相应输入量模糊集合;
采用重心法对输出量模糊集进行去模糊化,则得到模糊推理机的实际输出为:
u = Σ i , j f ij u ij Σ i , j f ij = Σ k = ( i , i + 1 ) t = ( j , j + 1 ) ( A k ( e ) B t ( e · ) ) u kt Σ k = ( i , i + 1 ) t = ( j , j + 1 ) ( A k ( e ) B t ( e · ) ) - - - ( 10 )
式中,ukt∈{-3,-2,-1,0,1,2,3}为推理机模糊输出值的数值化表示,u∈R为模糊PID控制器输出,其他参数定义见公式(9);
定义模糊PID控制器参数kP、kI和kD的调整算式如下:
k P = k P 0 + f { e , e · } k P = k P 0 + Δ k P k I = k I 0 + f { e , e · } k I = k I 0 + Δ k I k D = k D 0 + f { e , e · } k D = k D 0 + Δ k D - - - ( 11 )
其中,kP为模糊PID控制器的输出参数中的比例系数;kI为模糊PID控制器的输出参数中的积分系数;kD为模糊PID控制器的输出参数中的微分系数;
比例系数kP增加,则系统的响应速度快,但超调量增加;积分系数kI增加,则系统的稳态误差消除快,但过渡过程超调大;微分系数kD增加,则系统动态特性好,但抗干扰能力变差;
kP0、kI0和kD0为相应参数的初始值;ΔkP、ΔkI和ΔkD为在线控制过程中,由模糊PID控制器输出的参数调整量;在机构运动过程中,在线检测加载支链实际的驱动力和指令值之间的差值e和
Figure FDA0000075468040000053
由上式完成模糊PID控制器的参数调整;整个控制的流程图如图7;
b.反馈系统控制
该六维力高频疲劳试验机的反馈系统也应用力伺服控制;该力伺服控制将压力传感器安装在6-UPS并联机构加载支链上,通过6个加载支链上的压力传感器得到力反馈信号,通过对6个加载的反馈信号的处理,进而对6-UPS并联机构的动平台进行力伺服控制,该力伺服控制的优点在于单个力输出不存在相互干扰的问题,通过6-UPS并联机构力的反解很容易得到加载平台合力的结果,再通过模糊PID控制器就对当前信号进行对比、分析、控制;并且应用力伺服控制时,各个加载支链的互换性较好;
c.系统软件
该六维力高频疲劳试验机软件部分采用分层结构,按照功能分成4个部分,分别完成人机界面层、数据/参数管理层、通讯层及控制层;
人机界面层主要完成运动模式选择、运动参数输入、运动启动/停止控制的操作任务,同时需要与数据库连接,将系统控制设置参数以及运动过程中的状态参数显示,以便对运动过程进行人工干预,因此它的关联部分为数据库;
数据/参数管理层实际上就是一个数据库,需要管理的数据主要包括:各类系统控制参数、默认参数及运动状态;
通讯层主要任务就是根据具体的硬件进行时序协调管理;
在该软件系统中,控制层即对运动控制卡及筒式直线电机驱动器的编程及调试;为了使软件系统达到良好的控制效果,要保证运动的稳定性、准确性及快速性,控制层的主要工作包括筒式直线电机驱动器与运动控制卡连接,系统反馈控制与运动控制卡连接以及6-UPS并联机构本身位置与运动控制卡的连接;
为了增强系统软件的实时性和系统软件适应高频加载试验的需要,利用多线程技术对软件进行编写;利用多线程来优化系统软件,就是把程序的循环模块单独放在一个线程中来执行,把它与主线程隔离开,让主线程来完成循环模块以外的所有其他功能;循环模块所在的线程是一个高速运行的线程,依靠循环完成不停采集数据和控制等功能;而窗体的主线程除了负责循环模块以外的其他功能外,还负责响应用户输入;这样可以很好的提高系统软件的实时性;
2.根据权利要求1所述的六维力高频疲劳试验机,其特征在于:该工控机采用全钢机箱、全长主板、无源底板及工业电源;
3.根据权利要求1所述的六维力高频疲劳试验机,其特征在于:该运动控制卡采用6轴以上可编程多轴运动卡,其单轴最快的采样周期为60μs;
4.根据权利要求1所述的六维力高频疲劳试验机,其特征在于:该筒式直线电机采用筒式直线电机;该筒式直线电机安装在6-UPS并联机构的加载支链上,一端连接拉压力传感器,另外一端连接虎克铰链,为该六维力高频疲劳试验机提供动力;
5.根据权利要求1所述的六维力高频疲劳试验机,其特征在于:该筒式直线电机是三相无轴承水冷同步伺服直线电机,其最大推力不小于850N,持续推力不小于550N,最大功率不小于1400W,电机空载最大加速度不小于350m/s2
6.根据权利要求1所述的六维力高频疲劳试验机,其特征在于:该拉压力传感器采用柱式高精度拉压力传感器;该拉压力传感器一端与球铰链连接,另外一端与筒式直线电机连接;
7.根据权利要求1所述的六维力高频疲劳试验机,其特征在于:该拉压力传感器的弹性体为柱式结构,用于拉伸力和压缩力测量;
8.根据权利要求1所述的六维力高频疲劳试验机,其特征在于:该拉压力传感器采用全密封结构,其要求输出灵敏度为1.5mV/V,直线度为±0.05%F.S,重复性为±0.05%F.S,最大量程为±5000N;
9.一种六维力高频疲劳试验机的使用方法,其特征在于,包括下述步骤:
步骤一:分析试样在实际工作中的受力状况和建立试样的时变载荷谱;
分析确定试样在工作中的受力状况,包括集中载荷作用点及方向,分布载荷的分布规律,试样的约束情况,以及上述两类载荷状况随时间的变化规律,根据载荷的空间、时间分布规律,建立试样的载荷谱:F=F(x,y,z,t),即载荷F与载荷作用点x,y,z和时间t之间的函数关系;
步骤二:求出试样应力应变分布及其时间函数;
依据载荷谱采用有限元方法计算出结构件的应力、应变在试样中的分布,求得这种分布状况随时间的变化规律:σ=σ(x,y,z,t),ε=ε(x,y,z,t),即应力σ和应变ε的空间位置x,y,z分布与时间t之间的函数关系;
步骤三:采用有限元仿真技术确定这种应力状态下的加载点位移及其时间函数;
根据步骤二求出的应力应变分布及其时间函数,在有限元仿真环境下,采用单点集中力给试样加载,产生与工作状态尽量相似的应力应变分布及其随时间的变化,确定加载点位置、载荷大小及方向、产生的位移大小及方向,以及位移的时间函数;
步骤四:按照加载位移的时间函数进行多维力加载,根据加载位移的时间函数和模糊PID控制,确定加载平台运动要求,由工控机根据加载平台的运动求解6-UPS并联机构反解,得出每根并联加载支链的运动形式,筒式直线电机系统产生加载运动;
为求运动反解,在加载平台建立坐标系(O′-X′Y′Z′),在静平台建立坐标系(O-XYZ),每一根杆用一个矢量Li表示,其中,i=1、2、3、4、5或6;
H是加载平台坐标原点在静平台坐标系中的位置矢量,bi是加载平台铰链点位置矢量,Bi是静平台铰链点位置矢;由矢量关系得到:
Bi=H+biRT-Li    (12)
式中,R是加载平台座标系相对于静平台坐标系的姿态矩阵:
R = c θ 2 c θ 3 s θ 1 s θ 2 c θ 3 - c θ 1 s θ 3 c θ 1 s θ 2 c θ 3 + s θ 1 s θ 3 c θ 2 s θ 3 s θ 1 s θ 2 s θ 3 + c θ 1 c θ 3 c θ 1 s θ 2 s θ 3 - s θ 1 c θ 3 - s θ 2 sθ 1 c θ 2 c θ 1 c θ 2 - - - ( 13 )
矩阵中sin和cos简记为s和c,θ1、θ2、θ3分别表示加载平台绕x、y、z轴的转角;
根据加载平台加载时所需的位置H和姿态θ1、θ2、θ3,由式12和13,算出静平台铰链点Bi位置,实现对6-UPS并联机构的控制;
控制系统根据所算出的一系列杆长值,控制各杆产生相应位移,使加载平台实现加载运动,然后通过拉压力传感器的反馈,应用力伺服控制,对加载运动所产生的加载力和力矩与试样在工作中的受力状况进行对比,经过误差计算最后确定加载支链杆长的变化值,使得加载力的大小和方向与试样在工作中的受力状况相同;
步骤五:分析试样的实际应力分布,作为结构设计的参考;
进行多维力加载的同时,采用粘贴应变片等措施在试样关键位置监测应力应变,采样系统采集应力应变值,与载荷谱的理论值对比;分析试样的实际应力分布,作为试样结构设计的参考。
CN201110195523.2A 2011-07-13 2011-07-13 六维力高频疲劳试验机的使用方法 Expired - Fee Related CN102279101B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201110195523.2A CN102279101B (zh) 2011-07-13 2011-07-13 六维力高频疲劳试验机的使用方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201110195523.2A CN102279101B (zh) 2011-07-13 2011-07-13 六维力高频疲劳试验机的使用方法

Publications (2)

Publication Number Publication Date
CN102279101A true CN102279101A (zh) 2011-12-14
CN102279101B CN102279101B (zh) 2014-10-08

Family

ID=45104672

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201110195523.2A Expired - Fee Related CN102279101B (zh) 2011-07-13 2011-07-13 六维力高频疲劳试验机的使用方法

Country Status (1)

Country Link
CN (1) CN102279101B (zh)

Cited By (18)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103020351A (zh) * 2012-12-10 2013-04-03 中国飞机强度研究所 一种飞机位姿的三维实时显示方法
CN103020350A (zh) * 2012-12-10 2013-04-03 中国飞机强度研究所 一种紧固件载荷变形曲线数值仿真方法
CN103543029A (zh) * 2013-11-05 2014-01-29 中国航空工业集团公司西安飞机设计研究所 一种用于飞机全机主操纵系统疲劳试验的监控分析方法
CN104006959A (zh) * 2014-06-18 2014-08-27 苏州东菱科技有限公司 一种复合疲劳试验的测量方法
CN104564913A (zh) * 2014-12-10 2015-04-29 黑龙江真美广播通讯器材有限公司 并联六自由度转台液压控制系统
CN104748959A (zh) * 2015-04-09 2015-07-01 上海理工大学 三轴多维度加载力学性能测试试验台
CN105159137A (zh) * 2015-07-27 2015-12-16 江苏大学 一种混联式汽车电泳涂装输送机构动力学建模方法
CN106066903A (zh) * 2016-01-26 2016-11-02 西北工业大学 一种面向机器人末端执行器的结构强度校核方法
CN106570199A (zh) * 2015-10-10 2017-04-19 北汽福田汽车股份有限公司 用于滑动轴承的轴承座或轴颈的结构强度计算方法和装置
CN109129158A (zh) * 2018-10-30 2019-01-04 吉林大学 基于新型并联工具系统的精密铣磨成形机床及其控制方法
CN109799150A (zh) * 2019-02-19 2019-05-24 北京建筑大学 铁路i型弹条高频疲劳测试方法
CN111215648A (zh) * 2020-01-17 2020-06-02 清华大学 电主轴可靠性快速实验加载方法及加载系统
CN111409103A (zh) * 2020-03-27 2020-07-14 伯朗特机器人股份有限公司 工业机器人位姿特性中互换性的测量方法
CN111999046A (zh) * 2020-07-09 2020-11-27 东风汽车集团有限公司 一种疲劳耐久试验台架及试验台
CN112613159A (zh) * 2020-11-27 2021-04-06 北华航天工业学院 一种Stewart平台6-UHU结构运动学求解及误差标定方法
CN114608817A (zh) * 2022-03-01 2022-06-10 中航航空模拟系统有限公司 六自由度运动平台电动缸故障诊断方法
CN115256353A (zh) * 2022-09-05 2022-11-01 中国科学院光电技术研究所 一种外满秩、冗余多平台复合型串并联闭环式运动机构
CN118482855A (zh) * 2024-05-20 2024-08-13 湖北省计量测试技术研究院 一种机器人六维力传感器在线校准装置及主动驱动校准方法

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN2075574U (zh) * 1990-10-27 1991-04-24 东北重型机械学院秦皇岛分校 一种并联多环结构机器人
CN101451933A (zh) * 2008-12-18 2009-06-10 北京林业大学 六维加载材料力学试验机
US20090321631A1 (en) * 2008-06-25 2009-12-31 Axcelis Technologies, Inc. Low-inertia multi-axis multi-directional mechanically scanned ion implantation system
CN102063122A (zh) * 2010-11-10 2011-05-18 哈尔滨工业大学 一种空间六自由度运动台模态控制方法
CN102095341A (zh) * 2009-12-15 2011-06-15 上海宇航系统工程研究所 一种具有测量基准和固定坐标系的六自由度并联机构

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN2075574U (zh) * 1990-10-27 1991-04-24 东北重型机械学院秦皇岛分校 一种并联多环结构机器人
US20090321631A1 (en) * 2008-06-25 2009-12-31 Axcelis Technologies, Inc. Low-inertia multi-axis multi-directional mechanically scanned ion implantation system
CN101451933A (zh) * 2008-12-18 2009-06-10 北京林业大学 六维加载材料力学试验机
CN102095341A (zh) * 2009-12-15 2011-06-15 上海宇航系统工程研究所 一种具有测量基准和固定坐标系的六自由度并联机构
CN102063122A (zh) * 2010-11-10 2011-05-18 哈尔滨工业大学 一种空间六自由度运动台模态控制方法

Cited By (27)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103020350A (zh) * 2012-12-10 2013-04-03 中国飞机强度研究所 一种紧固件载荷变形曲线数值仿真方法
CN103020351A (zh) * 2012-12-10 2013-04-03 中国飞机强度研究所 一种飞机位姿的三维实时显示方法
CN103020351B (zh) * 2012-12-10 2016-04-20 中国飞机强度研究所 一种飞机位姿的三维实时显示方法
CN103543029A (zh) * 2013-11-05 2014-01-29 中国航空工业集团公司西安飞机设计研究所 一种用于飞机全机主操纵系统疲劳试验的监控分析方法
CN103543029B (zh) * 2013-11-05 2016-06-01 中国航空工业集团公司西安飞机设计研究所 一种用于飞机全机主操纵系统疲劳试验的监控分析方法
CN104006959A (zh) * 2014-06-18 2014-08-27 苏州东菱科技有限公司 一种复合疲劳试验的测量方法
CN104564913A (zh) * 2014-12-10 2015-04-29 黑龙江真美广播通讯器材有限公司 并联六自由度转台液压控制系统
CN104748959B (zh) * 2015-04-09 2017-03-22 上海理工大学 三轴多维度加载力学性能测试试验台
CN104748959A (zh) * 2015-04-09 2015-07-01 上海理工大学 三轴多维度加载力学性能测试试验台
CN105159137B (zh) * 2015-07-27 2018-11-06 江苏大学 一种混联式汽车电泳涂装输送机构动力学建模方法
CN105159137A (zh) * 2015-07-27 2015-12-16 江苏大学 一种混联式汽车电泳涂装输送机构动力学建模方法
CN106570199A (zh) * 2015-10-10 2017-04-19 北汽福田汽车股份有限公司 用于滑动轴承的轴承座或轴颈的结构强度计算方法和装置
CN106570199B (zh) * 2015-10-10 2020-03-24 北京宝沃汽车有限公司 用于滑动轴承的轴承座或轴颈的结构强度计算方法和装置
CN106066903B (zh) * 2016-01-26 2019-03-29 西北工业大学 一种面向机器人末端执行器的结构强度校核方法
CN106066903A (zh) * 2016-01-26 2016-11-02 西北工业大学 一种面向机器人末端执行器的结构强度校核方法
CN109129158A (zh) * 2018-10-30 2019-01-04 吉林大学 基于新型并联工具系统的精密铣磨成形机床及其控制方法
CN109799150B (zh) * 2019-02-19 2021-04-20 北京建筑大学 铁路i型弹条高频疲劳测试方法
CN109799150A (zh) * 2019-02-19 2019-05-24 北京建筑大学 铁路i型弹条高频疲劳测试方法
CN111215648A (zh) * 2020-01-17 2020-06-02 清华大学 电主轴可靠性快速实验加载方法及加载系统
CN111409103A (zh) * 2020-03-27 2020-07-14 伯朗特机器人股份有限公司 工业机器人位姿特性中互换性的测量方法
CN111999046A (zh) * 2020-07-09 2020-11-27 东风汽车集团有限公司 一种疲劳耐久试验台架及试验台
CN111999046B (zh) * 2020-07-09 2021-05-11 东风汽车集团有限公司 一种疲劳耐久试验台架及试验台
CN112613159A (zh) * 2020-11-27 2021-04-06 北华航天工业学院 一种Stewart平台6-UHU结构运动学求解及误差标定方法
CN112613159B (zh) * 2020-11-27 2023-10-03 北华航天工业学院 一种Stewart平台6-UHU结构运动学求解及误差标定方法
CN114608817A (zh) * 2022-03-01 2022-06-10 中航航空模拟系统有限公司 六自由度运动平台电动缸故障诊断方法
CN115256353A (zh) * 2022-09-05 2022-11-01 中国科学院光电技术研究所 一种外满秩、冗余多平台复合型串并联闭环式运动机构
CN118482855A (zh) * 2024-05-20 2024-08-13 湖北省计量测试技术研究院 一种机器人六维力传感器在线校准装置及主动驱动校准方法

Also Published As

Publication number Publication date
CN102279101B (zh) 2014-10-08

Similar Documents

Publication Publication Date Title
CN102279101B (zh) 六维力高频疲劳试验机的使用方法
CN103558079B (zh) 基于并联机构驱动力闭环的多自由度加载方法
Dongsu et al. Adaptive sliding control of six-DOF flight simulator motion platform
Duchaine et al. Computationally efficient predictive robot control
Lu et al. Real-time simulation system for UAV based on Matlab/Simulink
She et al. On the impact force of human-robot interaction: Joint compliance vs. link compliance
EP2105263A2 (en) Real time control method and device for robots in virtual projection
Liégeois Performance and computer-aided design
Chen et al. Design and dynamics modeling of a novel 2R1T 3-DOF parallel motion simulator
Huerzeler et al. Configurable real-time simulation suite for coaxial rotor uavs
Antonya et al. Co-simulation environment for the analysis of the driving simulator’s actuation
CN109241553B (zh) 一种机械臂关节实时动力学建模方法
Tang et al. Development of a laboratory HILs testbed system for small UAV helicopters
Post Robust state estimation for the control of flexible robotic manipulators
Moormann et al. The Modelica flight dynamics library
Aponso et al. Identification of higher order helicopter dynamics using linear modeling methods
Fan et al. Dynamic modelling and model predictive control of flexible-link manipulators
Aghili A mechatronic testbed for revolute-joint prototypes of a manipulator
Chen et al. Flight dynamics modelling and experimental validation for unmanned aerial vehicles
Nye et al. Design and modeling of a computer-controlled planar manipulator
Xu et al. A Novel Cable-Driven Joint Module for Space Manipulators: Design, Modeling, and Characterization
Pascault Finite-time tracking of trajectories in task space for stewart platform
Guida et al. Effect of Strain Wave Gear Teeth Wear on the Behaviour of an Electro-Mechanical Flight Control Actuator for eVTOL Aircrafts
Da Silva et al. Joint dynamics modeling and parameter identification for space robot applications
Wang et al. Friction Identification of Heavy-Duty Manipulator CMOR Using Markov Chain Monte Carlo Method

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20141008

Termination date: 20160713

CF01 Termination of patent right due to non-payment of annual fee