CN107860548B - 一种近似在线的颤振边界预测方法 - Google Patents
一种近似在线的颤振边界预测方法 Download PDFInfo
- Publication number
- CN107860548B CN107860548B CN201710816564.6A CN201710816564A CN107860548B CN 107860548 B CN107860548 B CN 107860548B CN 201710816564 A CN201710816564 A CN 201710816564A CN 107860548 B CN107860548 B CN 107860548B
- Authority
- CN
- China
- Prior art keywords
- speed
- frequency
- data
- flutter
- test
- 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.)
- Expired - Fee Related
Links
- 238000000034 method Methods 0.000 title claims abstract description 54
- 238000012360 testing method Methods 0.000 claims abstract description 39
- 230000004044 response Effects 0.000 claims abstract description 18
- 238000012545 processing Methods 0.000 claims abstract description 15
- 238000007664 blowing Methods 0.000 claims abstract description 14
- 238000013213 extrapolation Methods 0.000 claims abstract description 12
- 238000012886 linear function Methods 0.000 claims abstract description 4
- 230000008569 process Effects 0.000 claims description 4
- 238000013016 damping Methods 0.000 description 18
- 230000005284 excitation Effects 0.000 description 7
- 241001123248 Arma Species 0.000 description 5
- 238000004458 analytical method Methods 0.000 description 5
- 238000010586 diagram Methods 0.000 description 4
- 238000001228 spectrum Methods 0.000 description 4
- RYGMFSIKBFXOCR-UHFFFAOYSA-N Copper Chemical compound [Cu] RYGMFSIKBFXOCR-UHFFFAOYSA-N 0.000 description 2
- 238000003491 array Methods 0.000 description 2
- 230000008859 change Effects 0.000 description 2
- 229910052802 copper Inorganic materials 0.000 description 2
- 239000010949 copper Substances 0.000 description 2
- 230000007547 defect Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 230000035945 sensitivity Effects 0.000 description 2
- 241001156002 Anthonomus pomorum Species 0.000 description 1
- 238000012935 Averaging Methods 0.000 description 1
- 230000002238 attenuated effect Effects 0.000 description 1
- 238000005452 bending Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 230000000295 complement effect Effects 0.000 description 1
- 238000012790 confirmation Methods 0.000 description 1
- 230000008878 coupling Effects 0.000 description 1
- 238000010168 coupling process Methods 0.000 description 1
- 238000005859 coupling reaction Methods 0.000 description 1
- 230000003247 decreasing effect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000006073 displacement reaction Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 238000010998 test method Methods 0.000 description 1
- 238000010200 validation analysis Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01M—TESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
- G01M9/00—Aerodynamic testing; Arrangements in or on wind tunnels
- G01M9/02—Wind tunnels
- G01M9/04—Details
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01M—TESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
- G01M7/00—Vibration-testing of structures; Shock-testing of structures
- G01M7/02—Vibration-testing by means of a shake table
- G01M7/025—Measuring arrangements
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Fluid Mechanics (AREA)
- Aerodynamic Tests, Hydrodynamic Tests, Wind Tunnels, And Water Tanks (AREA)
Abstract
本发明公开一种近似在线的颤振边界预测方法,在低速风洞中进行连续风速下的吹风试验,测量出飞机关键部位的振动响应;画出各个模态频率峰值倒数图,下降趋势明显的即为危险模态;继续进行连续风速下的吹风试验,采集数据,只在危险模态的频率附近读取出最大值,记为频率峰值,每个数据点有一一对应的测量风速,该点数据所对应的速度取平均值,即为与峰值匹配的速度;画出峰值倒数与速度的关系图,并用一次函数进行拟合外推得到颤振边界速度;继续采集并处理新的试验数据,不断得到新的颤振边界速度,当所预测的颤振边界速度与实际吹风速度之间的相对误差不足3%时,停止试验。本发明实现在连续风速下进行颤振边界预测,是一种近似在线的预测方法。
Description
技术领域:
本发明涉及一种近似在线的颤振边界预测方法,其属于信息传输与处理技术领域。
背景技术:
颤振是由于气动力、弹性力和惯性力耦合而发生的振幅不衰减的自激振动,它是气动弹性力学中最重要的问题之一。
经过几十年的发展,颤振风行试验技术已经有了很大提高。但依然由三个部分组成:结构激励、响应测量和颤振边界的预测。其中,颤振边界预测是最为关键的一部分。包括阻尼外推法、颤振余度法、包线函数法、ARMA法、NG方法和鲁棒颤振分析法。
阻尼外推法通过观测各阶模态的阻尼随亚临界速度的变化直接外推到阻尼为零的点来求颤振临界速度。对于大型复杂系统,阻尼外推法是最效的,目前在实际的颤振飞行试验中使用也最广泛。但是,阻尼识别精度相比频率识别精度要差很多,为了减小阻尼识别的误差,需要在同一个速度下进行多次测试;同时阻尼是飞行速度的非线性函数,在外推时,需要尽量多的测试速度点才能保证结果的准确性,还需要颤振飞行试验点尽量接近颤振临界点,并选择合适的插值方法才能获得较高的预测精度。
颤振余度法由Zimmerman和Weissenburger提出,其基于经典的弯扭二自由度机翼的运动方程,利用Routh判据确定系统的稳定条件,并通过采样数据获得极点信息,求解描述颤振余度的与动压有关的二次多项式得到颤振临界点。与阻尼外推法一样,颤振余度法也需要进行参数识别,提取系统的阻尼和频率信息。由于颤振余度对阻尼的敏感程度小于对频率的敏感程度,而频率的识别精度相比阻尼要高,可以说,颤振余度法在一定程度上弥补了阻尼外推法的阻尼提取精度要求高、外推方法选择严和突发型颤振预测难等缺陷。但对于多自由度模型,须提前知道哪两阶模态发生耦合,但由于未考虑到其他模态的影响,必然存在一定的误差。
包线函数法由Copper提出,其理论基础是:对于任何稳定的有阻尼系统,脉冲响应是随着时间衰减的。在时域内,响应的衰减所围成的形状称为衰减包线,其大小和形状的不同意味着阻尼的变化,其本质内涵与阻尼外推法类似。包线函数法并不需要识别模态参数,而是通过计算不同速度下包线函数的形状参数来预测系统的稳定性。该方法对测试速度下的每一个监测位置的响应都可计算出一个包线函数,通过平均对提高预测精度有一定帮助。在包线函数法的基础上,Copper等人基于在线差分方程,以随机激励作为输入发展了在线的包线函数法,然而由于没有确定性的输入,该方法需要花费大量的时间收敛到一个稳定的估计,因此对于在线估计,其应用是有限的。
ARMA方法由Matsuzaki和Ando提出,所需的信号由大气紊流对结构产生随机激励而获得。其采用由包含响应项的自回归(AR)模型和代表着白噪声激励的滑动平均(MA)模型“混合”而成的ARMA模型来描述气动系统,通过系统辨识得到自回归项系数后,由Jury判据确定系统的稳定性。该方法对于二阶及二阶以上模态都适用,对于高阶模态系统,具有更高的预测精度。但在测试速度较小的时候,曲线拟合时可能会把缓和型颤振误作突发性颤振。利用ARMA方法进行颤振边界预测,飞行试验时无需额外的激振设备,操作简单,易于实现,但大气紊流激励能量主要集中在低频段,高频模态难以激起,并且响应信号的信噪比一般较低,这是ARMA方法的一大弊端。
NG法由Nissim和Gilyard提出,是一种基于气动弹性系统两个不同亚临界速度下的模态位移响应,通过辨识结构运动方程在亚临界来流动压下运动方程的系数,然后求解不同速度下已被辨识的系统特性获得颤振临界点的方法。NG法只需两个不同测试速度点即可预测出颤振临界特性。当测试动压离颤振临界动压较远时,依然能够准确地预测出颤振临界特性。增加激励力列阵的数目对提高该方法的鲁棒性和预测精度有很大作用,但当出现重频模态或者频率值十分接近时,需提供多个不同的激励力列阵。另外,若刚体的大幅运动掩盖了弹性响应,会使得辨识过程中高频模态被噪声所覆盖。
鲁棒颤振分析法是一种新的颤振边界预测方法,由Lind和Brenner提出,与传统颤振飞行试验方法不同,该方法通过结构奇异值理论将颤振理论模型和飞行试验数据结合起来,以结构奇异值来衡量气动弹性系统的鲁棒稳定性。模型确认是鲁棒颤振分析中非常关键的环节,不确定模型需能够合理反映标称模型与真实系统之间的差异,既不能过雨保守,也不能过于乐观。鲁棒颤振分析方法将理论模型与试验数据结合起来,使得理论与试验优势互补,且抗噪能力更强,可避免由于系统辨识失真造成稳定飞行边界预测的失败。但须注意,该方法具有一定保守性,在模型确认中,不确定性的来源,描述形式的选择和不确定性幅值的大小等都会影响鲁棒颤振分析结果。
发明内容:
目前的颤振边界预测方法都是在固定风速下进行预测的,本发明提供一种近似在线的颤振边界预测方法,是在连续风速下预测颤振边界,是一种近似在线的方法。
本发明所采用的技术方案有:一种近似在线的颤振边界预测方法,其特征在于:包括如下步骤:
(1).首先,在低速风洞中进行连续风速下的吹风试验,测量出飞机关键部位的振动响应;
在其中取1至n个数据点进行处理,利用频谱图识别出关键部位的主要模态,其频率记为f11,f12,f13,每个振动响应数据点有一一对应的测量风速,对该1-n点数据所对应的速度取平均值,记为速度v1,利用下面公式将时域的数据转换为频域:
然后在f11,f12,f13的频率附近读取出频域下数据的最大值,记为频率峰值m11,m12,m13;再取1+L到n+L的试验数据点,用同样的方法得到主要模态频率,记为f21,f22,f23,速度平均值记为v2,频率峰值m21,m22,m23;依次类推对1+2L至n+2L的试验数据进行处理,得到主要模态频率f31,f32,f33,计算出与之对应的平均速度v3,频率峰值m31,m32,m33;
(2).将f11,f21,f31记为第一阶模态频率,f12,f22,f32记为第二阶模态频率,f13,f23,f33记为第三阶模态频率,依次类推;
画出各个模态频率峰值倒数图1/m-v,其中下降趋势明显的即为危险模态;
(3).继续进行连续风速下的吹风试验,采集数据,取N个数据点进行处理,利用下面公式将时域的数据转换为频域:
只在危险模态的频率附近读取出最大值,记为频率峰值,每个数据点有一一对应的测量风速,该N点数据所对应的速度取平均值,即为与峰值匹配的速度;
(4).每隔M个数据,再取N个数据点按照步骤(2)进行处理,依次得到峰值-速度关系,画出峰值倒数与速度的关系图,并用一次函数进行拟合外推得到一个颤振边界速度;
(5).按照步骤(3)和(4)继续采集并处理新的试验数据,不断得到新的颤振边界速度,当所预测的颤振边界速度与实际吹风速度之间的相对误差不足3%时,停止试验。
本发明具有如下有益效果:频率的识别精度相比阻尼要高,本方法避免了利用阻尼进行颤振边界预测时精度要求高,实现难度大的缺陷。更重要的是,本发明可以实现在连续风速下进行颤振边界预测,是一种近似在线的预测方法。
附图说明:
图1为1-2000点的频谱图。
图2为2001-4000点的频谱图。
图3为4001-6000点的频谱图。
图4为6001-8000点的频谱图。
图5为风速与标定值之间成一次函数的关系图。
图6为第一阶模态的频率峰值倒数图。
图7为第二阶模态的频率峰值倒数图。
图8为根据峰值倒数—速度图像拟合外推颤振边界速度。
具体实施方式:
下面结合附图对本发明作进一步的说明。
本发明近似在线的颤振边界预测方法,包括如下步骤:
(1).首先,在低速风洞中进行连续风速下的吹风试验,测量出飞机关键部位的振动响应(如机翼或机身的弯曲、扭转等)。
在其中取1至n个数据点进行处理,利用频谱图识别出关键部位的主要模态,其频率记为f11,f12,f13;每个振动响应数据点有一一对应的测量风速,对该1-n点数据所对应的速度取平均值,记为速度v1,利用下面公式将时域的数据转换为频域:
然后在f11,f12,f13的频率附近读取出频域下数据的最大值,记为频率峰值m11,m12,m13。再取1+L到n+L的试验数据点,用同样的方法得到主要模态频率,记为f21,f22,f23;速度平均值记为v2,频率峰值m21,m22,m23。依次类推对1+2L至n+2L的试验数据进行处理,得到主要模态频率f31,f32,f33;计算出与之对应的平均速度v3,频率峰值m31,m32,m33。
(2).由于吹风试验中速度增长幅度不会很大,不同速度下得到的各阶模态频率相近,f11,f21,f31记为第一阶模态频率,f12,f22,f32记为第二阶模态频率,f13,f23,f33记为第三阶模态频率,依次类推。
另外,每个速度下的识别出的模态数可能是任意正整数,用同样的方法对数据进行处理时也可以进行多次。
因此,根据得到的数据可以整理如下所示的表格:
第一阶模态 | 第二阶模态 | 第三阶模态 | …… | |
v<sub>1</sub> | m<sub>11</sub> | m<sub>12</sub> | m<sub>13</sub> | …… |
v<sub>2</sub> | m<sub>21</sub> | m<sub>22</sub> | m<sub>23</sub> | …… |
v<sub>3</sub> | m<sub>31</sub> | m<sub>32</sub> | m<sub>33</sub> | …… |
…… | …… | …… | …… | …… |
v<sub>k</sub> | m<sub>k1</sub> | m<sub>k2</sub> | m<sub>k3</sub> | …… |
根据上面的表格,可以画出各个模态的频率峰值倒数图1/m-v。其中下降趋势明显的即为危险模态。
(3).继续进行连续风速下的吹风试验,采集数据,取N个数据点进行处理,利用下面公式将时域的数据转换为频域:
只在危险模态的频率附近读取出最大值,记为频率峰值;每个数据点有一一对应的测量风速,该N点数据所对应的速度取平均值,即为与峰值匹配的速度。
(4).每隔M个数据,再取N个数据点按照步骤(2)进行处理,依次得到峰值-速度关系。画出峰值倒数与速度的关系图,并用一次函数进行拟合外推得到一个颤振边界速度。M取的越小,拟合结果越精确。
(5).按照步骤(3)和(4)继续采集并处理新的试验数据,不断得到新的颤振边界速度。随着数据增加,结果将越来越精确。当所预测的颤振边界速度与实际吹风速度之间的相对误差不足3%时,停止试验。
下面通过具体实施例来说明本发明近似在线的颤振边界预测方法。
首先,在低速风洞中进行连续风速下的吹风试验。测量出飞机左翼扭转的振动响应。分别取1至2000、2001至4000、4001至6000和6001至8000个数据点进行处理,画出频谱图,如图1至图4所示:
根据频谱图读取出主要模态,并根据图5所示的风速与标定值之间的关系,找到振动响应数据点相对应的测量风速。
将得到的数据,整理为如下的表格:
根据表格,画出两个模态的频率下的峰值倒数图1/m-v,如图6和图7所示。可以看出:第一阶模态有更明显的下降趋势,为危险模态。
因此,在连续风速下左翼扭转的试验中,取第一阶模态进行近似在线的颤振边界预测。
首先取1-2000点,识别出危险模态下的峰值和及该段响应所对应的速度;每隔5个数据点,再取2000点,即:2006-4006、2011-4011……依次得到危险模态下的峰值和对应的速度。根据得到的峰值和速度来拟合图像,外推出一个颤振边界速度。
随着试验的进行,会得到更多的峰值及速度,便可以时时更新拟合图像和颤振边界速度,用一种近似在线的方法来进行颤振边界预测。
本次试验,在识别了1000组峰值和速度时,得到如图8,即为在连续风速下进行近似在线的颤振边界预测。
以上所述仅是本发明的优选实施方式,应当指出,对于本技术领域的普通技术人员来说,在不脱离本发明原理的前提下还可以作出若干改进,这些改进也应视为本发明的保护范围。
Claims (1)
1.一种近似在线的颤振边界预测方法,其特征在于:包括如下步骤:
(1).首先,在低速风洞中进行连续风速下的吹风试验,测量出飞机关键部位的振动响应;
在其中取1至n个数据点进行处理,利用频谱图识别出关键部位的主要模态,其频率记为f11,f12,f13,每个振动响应数据点有一一对应的测量风速,对该1-n点数据所对应的速度取平均值,记为速度v1,利用下面公式将时域的数据转换为频域:
然后在f11,f12,f13的频率附近读取出频域下数据的最大值,记为频率峰值m11,m12,m13;再取1+L到n+L的试验数据点,用同样的方法得到主要模态频率,记为f21,f22,f23,速度平均值记为v2,频率峰值m21,m22,m23;依次类推对1+2L至n+2L的试验数据进行处理,得到主要模态频率f31,f32,f33,计算出与之对应的平均速度v3,频率峰值m31,m32,m33;
(2).将f11,f21,f31记为第一阶模态频率,f12,f22,f32记为第二阶模态频率,f13,f23,f33记为第三阶模态频率,依次类推;
画出各个模态频率峰值倒数图1/m-v,其中下降趋势明显的即为危险模态;
(3).继续进行连续风速下的吹风试验,采集数据,取N个数据点进行处理,利用下面公式将时域的数据转换为频域:
只在危险模态的频率附近读取出最大值,记为频率峰值,每个数据点有一一对应的测量风速,该N点数据所对应的速度取平均值,即为与峰值匹配的速度;
(4).每隔M个数据,再取N个数据点按照步骤(2)进行处理,依次得到峰值-速度关系,画出峰值倒数与速度的关系图,并用一次函数进行拟合外推得到一个颤振边界速度;
(5).按照步骤(3)和(4)继续采集并处理新的试验数据,不断得到新的颤振边界速度,当所预测的颤振边界速度与实际吹风速度之间的相对误差不足3%时,停止试验。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710816564.6A CN107860548B (zh) | 2017-09-12 | 2017-09-12 | 一种近似在线的颤振边界预测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710816564.6A CN107860548B (zh) | 2017-09-12 | 2017-09-12 | 一种近似在线的颤振边界预测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107860548A CN107860548A (zh) | 2018-03-30 |
CN107860548B true CN107860548B (zh) | 2020-07-07 |
Family
ID=61699360
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710816564.6A Expired - Fee Related CN107860548B (zh) | 2017-09-12 | 2017-09-12 | 一种近似在线的颤振边界预测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107860548B (zh) |
Families Citing this family (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109063290A (zh) * | 2018-07-20 | 2018-12-21 | 中国航空工业集团公司沈阳飞机设计研究所 | 一种基于神经网络技术的颤振预测方法 |
CN109086501A (zh) * | 2018-07-20 | 2018-12-25 | 中国航空工业集团公司沈阳飞机设计研究所 | 一种颤振预测方法 |
CN109886464B (zh) * | 2019-01-20 | 2022-03-18 | 东北电力大学 | 基于优化奇异值分解生成特征集的低信息损失短期风速预测方法 |
CN110657939A (zh) * | 2019-08-30 | 2020-01-07 | 中国空气动力研究与发展中心高速空气动力研究所 | 颤振临界预测方法及装置 |
CN113267304A (zh) * | 2021-04-25 | 2021-08-17 | 上海机电工程研究所 | 导弹伺服振动亚临界试验及稳定边界预示系统及方法 |
CN113408357B (zh) * | 2021-05-21 | 2024-02-13 | 西北工业大学 | 一种颤振紊流响应信号到脉冲响应信号的生成方法 |
CN113945357B (zh) * | 2021-10-28 | 2022-09-06 | 南京航空航天大学 | 一种风力机叶片颤振临界风速预测方法及装置 |
CN114608795B (zh) * | 2022-05-11 | 2022-07-22 | 中国飞机强度研究所 | 一种飞机吹风测试用风洞系统共振边界确定方法 |
Family Cites Families (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103530511A (zh) * | 2013-10-10 | 2014-01-22 | 南京航空航天大学 | 紊流激励条件下的风洞颤振试验颤振边界预测方法 |
CN103617323A (zh) * | 2013-11-27 | 2014-03-05 | 中国飞行试验研究院 | 一种基于飞行试验的颤振数据处理架构 |
CN105509992A (zh) * | 2015-11-27 | 2016-04-20 | 中国航空工业集团公司沈阳飞机设计研究所 | 一种颤振风洞试验外挂物模型系统激励方法 |
US10073002B2 (en) * | 2016-03-03 | 2018-09-11 | United Technologies Corporation | Flutter detection sensor |
-
2017
- 2017-09-12 CN CN201710816564.6A patent/CN107860548B/zh not_active Expired - Fee Related
Also Published As
Publication number | Publication date |
---|---|
CN107860548A (zh) | 2018-03-30 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107860548B (zh) | 一种近似在线的颤振边界预测方法 | |
CN104443427B (zh) | 飞行器颤振预测系统及方法 | |
CN110823537B (zh) | 一种超高层建筑风振控制阻尼器性能的检测评价方法 | |
CN107590317A (zh) | 一种计及模型参数不确定性的发电机动态估计方法 | |
US10634580B2 (en) | Systems and methods for analyzing flutter test data using damped sine curve fitting with the closed form shape fit | |
CN102713263A (zh) | 现场校准风力涡轮机叶片的负载传感器的方法 | |
CN101859146A (zh) | 一种基于预测滤波和经验模态分解的卫星故障预测方法 | |
CN107391818A (zh) | 一种基于状态观测器的振动模态参数识别方法 | |
CN107688554A (zh) | 基于自适应傅里叶分解的频域辨识方法 | |
CN103530511A (zh) | 紊流激励条件下的风洞颤振试验颤振边界预测方法 | |
CN112199793B (zh) | 一种非高斯载荷的损伤等效试验谱优化方法 | |
CN104298870A (zh) | 一种移动荷载下简支梁损伤和移动力同时识别方法 | |
CN113901379B (zh) | 一种边缘端的实时数据动态在线快速处理方法 | |
CN104881585A (zh) | 一种两自由度机翼的颤振边界预测方法 | |
US20190376874A1 (en) | A method of estimating the number of modes for the sparse component analysis based modal identification | |
CN103310060A (zh) | 一种跨音速极限环颤振分析方法 | |
CN105628189A (zh) | 一种直升机机体振动信号的自适应处理方法 | |
CN110263762B (zh) | 一种基于输出的半潜式海洋平台能量传递路径分析方法 | |
CN115688288A (zh) | 飞行器气动参数辨识方法、装置、计算机设备及存储介质 | |
CN109062051B (zh) | 一种提高机器人动力学参数辨识精度的方法 | |
CN103559414A (zh) | 一种基于阿伦尼斯模型的环境因子点估计方法 | |
CN106406290B (zh) | 一种车辆横向动力远程测控系统的故障检测方法 | |
CN104535082A (zh) | 一种基于飞行试验和理论计算判断惯导元件性能的方法 | |
CN105446320B (zh) | 一种基于有限频域的高速列车垂向悬挂系统的故障检测方法 | |
CN115753002A (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 | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20200707 |