发明名称 一种航空遥感惯性稳定平台摩擦参数辨识方法
摘要 一种航空遥感惯性稳定平台摩擦参数辨识方法。在分析系统正弦响应曲线基础上,建立惯性稳定平台的LuGre摩擦模型,通过对惯性稳定平台三个框架系统分别设计参数辨识方案,利用两步辨识及动态参数优化的摩擦参数辨识方法得到LuGre摩擦模型参数。本发明能够较准确的得到LuGre摩擦模型参数,实现摩擦力矩补偿,提高惯性稳定平台系统的控制精度和稳定精度,可广泛适用于各类航空遥感惯性稳定平台系统。
申请公布号 CN103344243B 申请公布日期 2015.12.09
申请号 CN201310275237.6 申请日期 2013.07.02
申请人 北京航空航天大学 发明人 周向阳;宫国浩;刘炜;房建成;赵强;李永;岳海潇;俞瑞霞
分类号 G01C21/18(2006.01)I;G05B13/04(2006.01)I 主分类号 G01C21/18(2006.01)I
代理机构 北京科迪生专利代理有限责任公司 11251 代理人 成金玉;贾玉忠
主权项 一种航空遥感惯性稳定平台摩擦参数辨识方法,其特征在于:首先进行摩擦模型静态参数辨识,构造航空遥感惯性稳定平台三框架伺服系统PI闭环,使三框架伺服电机以不同转速正反转相同角度,低速运行,获得Stribeck负斜率段曲线,测得对应转速下的电流环电流,然后拟合Stribeck曲线,进行静态参数辨识,同时,对稳定平台三框架伺服系统输入适当电流,获得稳态角位移;其次,进行动态参数辨识,根据动力学模型计算动态参数初始值,利用动态参数初始值以稳态条件下实测摩擦力矩随速率变化曲线和摩擦模型仿真曲线的相似程度为目标函数进行动态参数优化,具体包括以下步骤:(1)采用K<sub>p</sub>+K<sub>i</sub>/S控制器对航空遥感惯性稳定平台三框架伺服电机系统构造PI控制速率、电流双闭环系统,K<sub>p</sub>+K<sub>i</sub>/S为拉普拉斯表示法下的控制器结构,K<sub>p</sub>为比例环节放大倍数,K<sub>i</sub>/S表示积分环节,K<sub>i</sub>为积分环节放大倍数;(2)在静基座下,三个框架伺服电机在负载及速率闭环情况下,分别给定不同的横滚框电机转速W<sub>L_R</sub>、俯仰框电机转速W<sub>L_P</sub>、方位框电机转速W<sub>L_A</sub>,在±0.5~±6°/s的转速下运行,使三个框架伺服电机正向或反向恒速转动同样角度,运行于Stribeck负斜率段曲线,测量不同转速恒速运行下对应的横滚框电流环电流I<sub>R</sub>、俯仰框电流环电流I<sub>P</sub>、方位框电流环电流I<sub>A</sub>,其中,Stribeck曲线在这里是框架转速与对应摩擦力矩的曲线;(3)采用非线性最小二乘方法,结合步骤(2)测得的速率和电流值,对三个框架的稳态LuGre摩擦模型<img file="FDA0000772228910000011.GIF" wi="830" he="180" />进行拟合,获得三个框架的稳态LuGre摩擦模型的静态参数T<sub>s</sub>,T<sub>c</sub>,ω<sub>s</sub>和σ<sub>2</sub>;T<sub>fs</sub>表示静基座稳态时的摩擦力矩,T<sub>s</sub>,T<sub>c</sub>,ω<sub>s</sub>和σ<sub>2</sub>为摩擦模型的四个静态参数,T<sub>s</sub>表示最大静摩擦力矩,T<sub>c</sub>表示库仑摩擦力矩,ω<sub>s</sub>表示临界Stribeck角速度,σ<sub>2</sub>表示粘性摩擦系数,ω<sub>l</sub>表示负载端角速度;(4)根据步骤(3)拟合得到的静态参数,俯仰轴、横滚轴电机输入电流使其伺服电机输出力矩介于T<sub>s</sub>+T<sub>i</sub>和T<sub>c</sub>+T<sub>i</sub>之间,其中T<sub>i</sub>为干扰力矩,测得横滚框、俯仰框静摩擦区域的对应稳态角位移θ<sub>S_R</sub>和θ<sub>S_P</sub>,方位轴电机输入电流使伺服电机输出力矩介于T<sub>s</sub>和T<sub>c</sub>之间,测得对应稳态角位移θ<sub>s_A</sub>;(5)动态参数辨识,首先根据公式σ<sub>0</sub>≈T<sub>c</sub>sgn(ω<sub>l</sub>)/θ<sub>s</sub>,ω<sub>l</sub>为负载端角速度,由步骤(3)、(4)所得静态参数,求出运动停止时三框架模型的动态参数初值σ<sub>0</sub>,然后根据<img file="FDA0000772228910000021.GIF" wi="424" he="87" />求得动态参数初值σ<sub>1</sub>,其中J为等效到负载端的转动惯量;(6)动态参数优化,以稳态条件下实测摩擦力矩随速率变化曲线和摩擦模型仿真曲线的相似程度目标函数<img file="FDA0000772228910000022.GIF" wi="460" he="116" />为动态参数优化的依据,利用Matlab优化工具箱,得到动态参数σ<sub>0</sub>与σ<sub>1</sub>优化值,式中T<sub>f</sub>(i)和T<sub>f</sub>(ω<sub>l</sub>(i)分别为稳态条件下实测摩擦力矩和稳态条件下仿真摩擦力矩。
地址 100191 北京市海淀区学院路37号