CN109635389B - 一种电动舵机刚度试验数据处理方法 - Google Patents
一种电动舵机刚度试验数据处理方法 Download PDFInfo
- Publication number
- CN109635389B CN109635389B CN201811445168.8A CN201811445168A CN109635389B CN 109635389 B CN109635389 B CN 109635389B CN 201811445168 A CN201811445168 A CN 201811445168A CN 109635389 B CN109635389 B CN 109635389B
- Authority
- CN
- China
- Prior art keywords
- load
- steering engine
- displacement
- electric steering
- data
- 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.)
- Active
Links
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/15—Vehicle, aircraft or watercraft design
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B64—AIRCRAFT; AVIATION; COSMONAUTICS
- B64F—GROUND OR AIRCRAFT-CARRIER-DECK INSTALLATIONS SPECIALLY ADAPTED FOR USE IN CONNECTION WITH AIRCRAFT; DESIGNING, MANUFACTURING, ASSEMBLING, CLEANING, MAINTAINING OR REPAIRING AIRCRAFT, NOT OTHERWISE PROVIDED FOR; HANDLING, TRANSPORTING, TESTING OR INSPECTING AIRCRAFT COMPONENTS, NOT OTHERWISE PROVIDED FOR
- B64F5/00—Designing, manufacturing, assembling, cleaning, maintaining or repairing aircraft, not otherwise provided for; Handling, transporting, testing or inspecting aircraft components, not otherwise provided for
- B64F5/60—Testing or inspecting aircraft components or systems
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Geometry (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Aviation & Aerospace Engineering (AREA)
- Mathematical Optimization (AREA)
- Mathematical Analysis (AREA)
- Computational Mathematics (AREA)
- Pure & Applied Mathematics (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Automation & Control Theory (AREA)
- Manufacturing & Machinery (AREA)
- Transportation (AREA)
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本申请公开了一种电动舵机刚度试验数据处理方法,其包括:获取电动舵机在系统载荷作用下与载荷呈线性关系的载荷‑位移数据;建立包含载荷变量、位移变量与随机误差变量的回归模型,根据载荷‑位移数据对回归模型进行一次回归拟合;选取不超过门限值的随机误差变量,并根据载荷‑位移数据二次求解回归模型获得电动舵机的载荷‑位移曲线。本申请通过采用线性回归分析法及引入误差消除因子,对电动舵机刚度试验数据进行分析处理,从而获得其刚度数据,并可以根据加载、卸载力‑位移曲线推算出舵机间隙。本申请可以预估舵机刚度,提高刚度试验效率及试验结果精度,从而提高了飞机设计仿真精度。
Description
技术领域
本申请属于飞行器试验技术领域,特别涉及一种电动舵机刚度试验数据处理方法。
背景技术
随着无人机技术在近二十年的快速发展,从微型飞行器、旋翼飞行器、常规布局无人机到隐身布局无人机,无人机平台呈现爆发式发展。在其中,对于中小型无人机平台,许多是采用电动舵机对操纵面进行操控的。然而无人机平台使用的电动舵机,其性能特性同大型飞行器通常使用的液压舵机还是存在着一定的区别的。因此,在无人机设计的过程中有必要对其刚度、间隙等性能参数进行精确测定。
发明内容
本申请的目的是提供了一种电动舵机刚度试验数据处理方法,以解决上述任一问题。
本申请的技术方案是:一种电动舵机刚度试验数据处理方法,其包括:获取电动舵机在系统载荷作用下与所述载荷呈线性关系的载荷-位移数据;建立包含载荷变量、位移变量与随机误差变量的回归模型,根据载荷-位移数据对回归模型进行一次回归拟合;选取不超过门限值的所述随机误差变量,并根据所述载荷-位移数据二次求解所述回归模型获得电动舵机的载荷-位移曲线。
本申请中,所述回归模型为:y=α+βx+ε
式中,x为测量点位移,y为系统的载荷,α为截距,β为斜率,即测量点当地的刚度系数,ε为随机误差项。
本申请中,所述门限值为所述随机误差变量的20%。
本申请中,所述载荷-位移数据求解所述回归模型采用最小二乘法求解。
本申请通过采用线性回归分析法及引入误差消除因子,对电动舵机刚度试验数据进行分析处理,从而获得其刚度数据,并可以根据加载、卸载力-位移曲线推算出舵机间隙。本申请的方法简单,数据处理方便,可对飞机电动舵机刚度试验进行数据分析处理,且可以预估舵机刚度,提高刚度试验效率及试验结果精度,从而提高了飞机设计仿真精度。
附图说明
为了更清楚地说明本申请提供的技术方案,下面将对附图作简单地介绍。显而易见地,下面描述的附图仅仅是本申请的一些实施例。
图1为本申请的电动舵机刚度试验数据处理方法流程图。
图2为本申请一实施例的电动舵机力-位移曲线图。
图3为本申请一实施例的电动舵机力-位移曲线拟合图。
具体实施方式
为使本申请实施的目的、技术方案和优点更加清楚,下面将结合本申请实施例中的附图,对本申请实施例中的技术方案进行更加详细的描述。
无人机电动舵机刚度数据是一种加载力与位移值的大量散点数据,为了能够精准获得电动舵的机刚度数据,本申请提出了一种电动舵机刚度试验数据处理方法。本方法可以从大量观测的散点数据中寻找到能反映事物内部的统计规律,并可以按数学模型形式表达出来,其主要是根据最小二乘法原理,通过对微分方程组求偏导数,接触各个常数项,从而最终得到定量公式。从而可以较精确地得到电动舵机在线性工作区间的刚度及评估舵机初始间隙。
如图1,对于飞机结构或部件的刚度特性,通常是通过刚度试验得到的。通过刚度试验测得加载点载荷以及各测量点位移,再通过数据处理,得到试验件的刚度数据。因此,试验数据测量的准确性及数据处理的可靠性对于试验结果的准确性十分重要。本申请主要是采用线性回归对试验数据进行初步分析,引入误差消除因子对误差较大的数据点进行去除,最后使用去除误差点的试验数据再次进行线性回归分析最终得到试验件的刚度数据。
在本申请用于对刚度试验的一组实测数据进行分析处理。对于刚度试验中用力与位移间线性关系表示的单变量关系,可以称为一元回归,其回归方程(回归模型)为:
y=α+βx+ε (1)
式中,x为自变量,表示测量点位移;y为因变量,表示系统的载荷;α为方程的截距;β为方程的斜率,即测量点当地的刚度系数;ε为随机误差项。
对于回归方程中的随机误差项ε,其具有如下性质:
1)随机误差项是随机变量,其期望值为0;
2)对于所有的位移值x,随机误差项的方差σ2为常数;
3)随机误差项之间相互独立,即与一个值相联系的误差对与另一个值相联系的误差没有影响;
4)随机误差项服从正态分布。
在线性回归分析中采用最小二乘法,即通过对微分方程组求偏导数,解出各常数项,从而得到定量公式。
对一组刚度试验测量值(yi,xi)(i=1,2,…,n),满足:
yi=α+βxi+εi (2)
由于在实际的试验条件下很可能出现一些误差较大的测量点,从而影响整体的试验结果。为了解决这一问题,本申请中引入了误差消除因子η对试验数据进行二次处理。对于误差消除因子η应综合考虑试验对象刚度的线性特征、试验数据的样本数量以及试验设备的测量精度等因素确定,本申请中,误差消除因子η不超过随机误差项的20%。即对于每一组测量值(yi,xi)与参数(α,β)的估计量之间都存在一个随机误差项εi。对误差消除因子η定义一个门限值,例如η=0.1,即将随机误差项εi最大的10%测量值剔除,从而保证了试验数据整体的精确度。
根据剔除误差较大数据后的刚度试验测量值(力-位移)进行二次线性回归分析,最终得到试验测量点的力-位移曲线,即该点的刚度值。
最后分别依据加载段、卸载段测量数据得到试验测量点的加载、卸载力-位移曲线。计算加载曲线和卸载曲线在位移轴上截距点的距离即得到电动舵机的间隙。
以某大展弦比缩比验证机的电动舵机刚度试验为例,在试验室条件下进行了方向舵舵机的刚度试验。将舵机固支于试验台上,在舵机通电的情况下,从正、反两个方向上施加不同的静态试验载荷,测量舵机的位移。其中加载过程为从正方向加载、卸载,然后反方向加载、卸载,最后正方向加载、卸载,测量加载和卸载时的位移。刚度试验结果如图1所示的力-位移曲线。
根据图2所示的试验数据,本申请的数据处理方法进行试验试验处理。首先,分别选取加载段和卸载段数据进行一次线性回归拟合;然后,根据试验的实际情况取误差消除因子η=0.125,即分别剔除加载段和卸载段试验数据中误差最大的一个点;之后将剔除误差较大点后的两组试验数据进行二次回归分析;最后得到加载段和卸载段的刚度拟合曲线。两条刚度拟合曲线在位移轴上的截距差即为电动舵机的间隙,拟合结果及间隙评估如图3所示。
本申请通过采用线性回归分析法及引入误差消除因子,对电动舵机刚度试验数据进行分析处理,从而获得其刚度数据,并可以根据加载、卸载力-位移曲线推算出舵机间隙。本申请的方法简单,数据处理方便,可对飞机电动舵机刚度试验进行数据分析处理,且可以预估舵机刚度,提高刚度试验效率及试验结果精度,从而提高了飞机设计仿真精度。
以上所述,仅为本申请的具体实施方式,但本申请的保护范围并不局限于此,任何熟悉本技术领域的技术人员在本申请揭露的技术范围内,可轻易想到的变化或替换,都应涵盖在本申请的保护范围之内。因此,本申请的保护范围应以所述权利要求的保护范围为准。
Claims (1)
1.一种电动舵机刚度试验数据处理方法,其特征在于,包括:
获取电动舵机在系统载荷作用下与所述载荷呈线性关系的载荷-位移数据;
建立包含载荷变量、位移变量与随机误差变量的回归模型,根据载荷-位移数据对回归模型进行一次回归拟合,其中,在线性回归分析中采用最小二乘法,即通过对微分方程组求偏导数,解出各常数项,从而得到定量公式,过程包括:
对一组刚度试验测量值(yi,xi)(i=1,2,…,n),满足回归模型:
yi=α+βxi+εi
式中,xi为测量点位移,yi为系统的载荷,α为截距,β为斜率,即测量点当地的刚度系数,εi为随机误差项;
选取不超过门限值数量的所述随机误差变量的数据,并根据所述载荷-位移数据二次求解所述回归模型获得电动舵机的载荷-位移曲线。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811445168.8A CN109635389B (zh) | 2018-11-29 | 2018-11-29 | 一种电动舵机刚度试验数据处理方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811445168.8A CN109635389B (zh) | 2018-11-29 | 2018-11-29 | 一种电动舵机刚度试验数据处理方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109635389A CN109635389A (zh) | 2019-04-16 |
CN109635389B true CN109635389B (zh) | 2022-12-20 |
Family
ID=66069935
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811445168.8A Active CN109635389B (zh) | 2018-11-29 | 2018-11-29 | 一种电动舵机刚度试验数据处理方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109635389B (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111003204B (zh) * | 2019-12-06 | 2022-12-02 | 江西洪都航空工业集团有限责任公司 | 一种飞机平尾伺服作动机构动刚度测试系统及测试方法 |
CN111003207B (zh) * | 2019-12-06 | 2022-12-02 | 江西洪都航空工业集团有限责任公司 | 一种飞机舵机动刚度试验台稳定裕度测试系统及方法 |
CN112014217A (zh) * | 2020-09-08 | 2020-12-01 | 湖南省力宇燃气动力有限公司 | 弹性体刚度确定方法 |
Citations (18)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2007094567A (ja) * | 2005-09-27 | 2007-04-12 | Fuji Heavy Ind Ltd | 剛性管理装置、及び、剛性管理方法 |
CN101105423A (zh) * | 2007-08-03 | 2008-01-16 | 北京理工大学 | 变刚度舵机仿真加载装置 |
CN102609578A (zh) * | 2012-01-20 | 2012-07-25 | 重庆长安汽车股份有限公司 | 一种汽车钢板弹簧多体模型的构建方法 |
CN202485892U (zh) * | 2011-12-20 | 2012-10-10 | 中国航空工业集团公司沈阳飞机设计研究所 | 一种颤振风洞试验模型舵面操纵刚度模拟装置 |
CN102829747A (zh) * | 2012-08-03 | 2012-12-19 | 北京机电工程研究所 | 一种适用于非线性舵系统的间隙测量方法 |
CN105004620A (zh) * | 2015-07-16 | 2015-10-28 | 浙江工业大学 | 一种高频疲劳试验机的动态载荷误差补偿方法 |
CN105092519A (zh) * | 2015-07-10 | 2015-11-25 | 东北大学 | 基于增量偏最小二乘法的样品成份测定方法 |
CN106096178A (zh) * | 2016-06-24 | 2016-11-09 | 哈尔滨大金工程试验检测有限公司 | 一种桥梁拉索抗弯刚度识别方法 |
CN106407575A (zh) * | 2016-09-23 | 2017-02-15 | 南京航空航天大学 | 一种复合材料柔性零件装配偏差分析方法 |
CN106828972A (zh) * | 2016-12-15 | 2017-06-13 | 中国航空工业集团公司西安飞机设计研究所 | 一种飞机舵面驱动作动器性能测试装置 |
JP2017187404A (ja) * | 2016-04-06 | 2017-10-12 | 住友ゴム工業株式会社 | ゴムの歪み予測方法 |
CN107403021A (zh) * | 2016-05-18 | 2017-11-28 | 北京自动化控制设备研究所 | 一种电动伺服系统刚度建模方法 |
CN107817081A (zh) * | 2017-08-31 | 2018-03-20 | 中国航空工业集团公司沈阳飞机设计研究所 | 一种座舱盖静力试验台刚度分布模拟方法 |
CN107844627A (zh) * | 2017-09-25 | 2018-03-27 | 北京理工大学 | 一种仅输出时变结构模态参数贝叶斯估计方法 |
CN108108522A (zh) * | 2017-11-29 | 2018-06-01 | 中国飞机强度研究所 | 一种加筋壁板极限载荷计算修正方法 |
CN108133075A (zh) * | 2017-11-29 | 2018-06-08 | 中国飞机强度研究所 | 一种航空发动机轴类部件支撑刚度模拟器 |
CN108520089A (zh) * | 2018-02-28 | 2018-09-11 | 哈尔滨工程大学 | 一种柔性连接的波浪滑翔器运动预测方法 |
CN108763656A (zh) * | 2018-05-03 | 2018-11-06 | 东南大学 | 一种基于完全二阶多项式响应面模型的含铰结构区间刚度识别方法 |
Family Cites Families (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
GB2467184A (en) * | 2009-01-27 | 2010-07-28 | Illinois Tool Works | Load testing apparatus |
WO2014158185A1 (en) * | 2013-03-29 | 2014-10-02 | Compagnie Generale Des Etablissements Michelin | Tire uniformity improvement using estimates based on convolution/deconvolution |
EP3044564A1 (en) * | 2013-09-09 | 2016-07-20 | MTS Systems Corporation | Test system having a compliant actuator assembly and iteratively obtained drive |
US9540121B2 (en) * | 2015-02-25 | 2017-01-10 | Cisco Technology, Inc. | Pre-flight self test for unmanned aerial vehicles (UAVs) |
DE102015014770A1 (de) * | 2015-11-13 | 2017-05-18 | GM Global Technology Operations LLC (n. d. Ges. d. Staates Delaware) | Verfahren zur Bestimmung elastischer Eigenschaften eines Kraftfahrzeug-Kabelbaums |
-
2018
- 2018-11-29 CN CN201811445168.8A patent/CN109635389B/zh active Active
Patent Citations (18)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2007094567A (ja) * | 2005-09-27 | 2007-04-12 | Fuji Heavy Ind Ltd | 剛性管理装置、及び、剛性管理方法 |
CN101105423A (zh) * | 2007-08-03 | 2008-01-16 | 北京理工大学 | 变刚度舵机仿真加载装置 |
CN202485892U (zh) * | 2011-12-20 | 2012-10-10 | 中国航空工业集团公司沈阳飞机设计研究所 | 一种颤振风洞试验模型舵面操纵刚度模拟装置 |
CN102609578A (zh) * | 2012-01-20 | 2012-07-25 | 重庆长安汽车股份有限公司 | 一种汽车钢板弹簧多体模型的构建方法 |
CN102829747A (zh) * | 2012-08-03 | 2012-12-19 | 北京机电工程研究所 | 一种适用于非线性舵系统的间隙测量方法 |
CN105092519A (zh) * | 2015-07-10 | 2015-11-25 | 东北大学 | 基于增量偏最小二乘法的样品成份测定方法 |
CN105004620A (zh) * | 2015-07-16 | 2015-10-28 | 浙江工业大学 | 一种高频疲劳试验机的动态载荷误差补偿方法 |
JP2017187404A (ja) * | 2016-04-06 | 2017-10-12 | 住友ゴム工業株式会社 | ゴムの歪み予測方法 |
CN107403021A (zh) * | 2016-05-18 | 2017-11-28 | 北京自动化控制设备研究所 | 一种电动伺服系统刚度建模方法 |
CN106096178A (zh) * | 2016-06-24 | 2016-11-09 | 哈尔滨大金工程试验检测有限公司 | 一种桥梁拉索抗弯刚度识别方法 |
CN106407575A (zh) * | 2016-09-23 | 2017-02-15 | 南京航空航天大学 | 一种复合材料柔性零件装配偏差分析方法 |
CN106828972A (zh) * | 2016-12-15 | 2017-06-13 | 中国航空工业集团公司西安飞机设计研究所 | 一种飞机舵面驱动作动器性能测试装置 |
CN107817081A (zh) * | 2017-08-31 | 2018-03-20 | 中国航空工业集团公司沈阳飞机设计研究所 | 一种座舱盖静力试验台刚度分布模拟方法 |
CN107844627A (zh) * | 2017-09-25 | 2018-03-27 | 北京理工大学 | 一种仅输出时变结构模态参数贝叶斯估计方法 |
CN108108522A (zh) * | 2017-11-29 | 2018-06-01 | 中国飞机强度研究所 | 一种加筋壁板极限载荷计算修正方法 |
CN108133075A (zh) * | 2017-11-29 | 2018-06-08 | 中国飞机强度研究所 | 一种航空发动机轴类部件支撑刚度模拟器 |
CN108520089A (zh) * | 2018-02-28 | 2018-09-11 | 哈尔滨工程大学 | 一种柔性连接的波浪滑翔器运动预测方法 |
CN108763656A (zh) * | 2018-05-03 | 2018-11-06 | 东南大学 | 一种基于完全二阶多项式响应面模型的含铰结构区间刚度识别方法 |
Non-Patent Citations (5)
Title |
---|
Variable stiffness composite material design by using support vector regression assisted efficient global optimization method;Fan Ye et,al;《Structural and Multidisciplinary Optimization》;20171231;第56卷(第1期);第203–219页 * |
一种全地形车橡胶扭力轴套扭转刚度建模及优化方法;姚寿文等;《北京汽车》;20170425(第02期);第8-13+45页 * |
基于遗传算法的电动负载模拟器ITAE控制器设计和仿真;黄勇等;《宇航学报》;20080930(第05期);第61-65页 * |
柔性悬臂梁的振动特性与等效线性化方法的局限性;李斌等;《航空学报》;20130308(第09期);第149-159页 * |
舵机动刚度测试系统设计;吴娟等;《弹箭与制导学报》;20090615(第03期);第42-44页 * |
Also Published As
Publication number | Publication date |
---|---|
CN109635389A (zh) | 2019-04-16 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109635389B (zh) | 一种电动舵机刚度试验数据处理方法 | |
CN107862134B (zh) | 一种考虑自相关测量误差的Wiener过程可靠性分析方法 | |
CN107436983B (zh) | 一种基于多元样本差异的o型橡胶密封圈寿命预测方法 | |
CN107885928B (zh) | 考虑测量误差的步进应力加速性能退化可靠性分析方法 | |
CN111859658A (zh) | 一种产品贮存寿命与可靠性评估方法 | |
CN111079329B (zh) | 一种基于相似结构试验的疲劳寿命评估方法 | |
CN112800616B (zh) | 基于比例加速退化建模的设备剩余寿命自适应预测方法 | |
CN112710690A (zh) | 一种确定加速寿命试验中加速因子的方法及应用 | |
CN110457731B (zh) | 一种应用于直升机动部件振动趋势分析的方法 | |
CN109583039A (zh) | 一种极小子样可靠性增长试验设计方法 | |
CN114018832B (zh) | 一种钢表面涂层防护等级的评估方法 | |
CN114001654B (zh) | 工件端面位姿评价方法 | |
CN109766518B (zh) | 考虑样本个体差异的不确定加速退化建模和分析方法 | |
CN113094640A (zh) | 一种频域下宽带多轴随机振动寿命预测方法 | |
CN110567487A (zh) | 一种激光陀螺可靠性测试系统及方法 | |
CN110736595B (zh) | 轨道交通舱室内电子设备寿命评估方法 | |
RU2419775C1 (ru) | Способ идентификации аэродинамических характеристик летательного аппарата по результатам испытаний | |
CN107732940B (zh) | 一种基于adpss的电力系统稳定器参数优化试验方法 | |
CN113984833B (zh) | 一种环境温度等效及加速试验方法 | |
CN112685912B (zh) | 一种多元广义Wiener过程性能退化可靠性分析方法 | |
Yang et al. | Accelerated life reliability evaluation of grating ruler for CNC machine tools based on competing risk model and incomplete data | |
CN114492074A (zh) | 一种概率损伤容限评估分析方法 | |
CN110920932B (zh) | 一种确定直升机安定面安装结构的载荷谱的方法及装置 | |
CN115186393A (zh) | 一种低样本量的航空发动机结构件c-p-s-n曲线获取方法 | |
CN109388829A (zh) | 一种电子产品寿命测算方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |