CN105021351A - 一种基于光纤光栅振动传感系统的直升机桨叶质量不平衡检测方法 - Google Patents
一种基于光纤光栅振动传感系统的直升机桨叶质量不平衡检测方法 Download PDFInfo
- Publication number
- CN105021351A CN105021351A CN201510502911.9A CN201510502911A CN105021351A CN 105021351 A CN105021351 A CN 105021351A CN 201510502911 A CN201510502911 A CN 201510502911A CN 105021351 A CN105021351 A CN 105021351A
- Authority
- CN
- China
- Prior art keywords
- blade
- fiber
- frequency
- vibration
- 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.)
- Granted
Links
Landscapes
- Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
Abstract
本发明提供了一种基于光纤光栅振动传感系统的直升机桨叶质量不平衡检测方法,将光纤光栅传感器布置于直升机的桨叶表面(1)和桨叶底座表面(2)上,通过光纤滑环和光缆连接到后部的解调仪和上位机,在桨叶转动时感受桨叶表面的应变和桨叶底座的振动(3)并采集数据(4),后部的数据处理系统通过时域分析(5)和频域分析(6),提取数据特征,与正常桨叶的数据进行对比,从而判断直升机桨叶的质量是否不平衡以及不平衡程度(7)。本发明克服了人工纸筒式标杆法费时费力、精度低、具有一定危险性和频闪仪靶标法依靠人眼不客观、定标时间长的缺点,能对直升机桨叶的质量不平衡进行客观、精确的检测和定量分析。
Description
技术领域
本发明属于结构健康监测领域,特别涉及一种基于光纤光栅振动传感系统的直升机桨叶质量不平衡检测方法。
背景技术
直升机旋翼作为向直升机飞行提供升力和操纵力的核心部件,直接关系到直升机的安全及其他各项重要性能的优劣。直升机旋翼高速旋转产生的椎体是直升机旋翼系统技术状态的综合反映,包括桨叶质量变化、桨叶安装和固定情况等,均反映到旋翼椎体的变化。其中桨叶的质量不平衡会引起飞机机身的振动,影响飞行安全,严重时有可能会引起桨叶的折断,造成机毁人亡的后果,所以为了飞行安全,很有必要对这一常见的问题进行检测。
测量桨叶质量不平衡的传统方法有两种:人工纸筒式标杆法和频闪仪靶标法。人工纸筒式标杆法是在各根桨叶的叶尖上涂不同的颜色,当旋翼旋转到规定转速时,将吊挂白色纸筒的标杆逐渐靠近桨叶叶尖,通过人工观察纸筒上不同颜色的痕迹来粗略判断桨叶质量的平衡状况。人工纸筒式标杆法较为简单,但是有几点弊端:1、人工靠杆需要空、地勤人员的协同,且靠杆有一定危险性;2、重复性差,无法记录数据;3、比较依赖熟练技术工的经验,精度低。频闪仪靶标法则是利用人眼视觉暂留现象发展出的一种测量方法,在各桨叶上粘贴特定的反光靶标,在直升机地面开车后,用频闪仪固定照射桨叶某一位置,反光靶标随之发光,由于人眼的迟滞性,会记录下靶标位置,若桨叶质量不平衡,转动轨迹变化,则靶标位置也会变化。此方法相对于人工纸筒式靶标法有所进步,但仍存在缺点:1、无法记录数据,依赖人眼判断,不够客观也不能进行定量分析;2、首次锁定靶标时间较长,有可能对直升机机体结构造成危害。
光纤光栅传感器作为一种新兴的传感器,具有小巧轻便、抗电磁干扰、耐腐蚀、灵敏度高、成本低及易于实现复用等优点,在传感领域发挥着重要作用,并在航空航天、桥梁、矿山、地质勘探、铁路检测等领域中得到广泛应用。目前将光纤光栅传感器应用在直升机桨叶质量不平衡检测上的实例和专利暂时还没有,但是光纤光栅传感器所拥有的这些优点,令其具有良好的应用前景,而此专利的发明内容,也是富有实践意义的。
发明内容
本发明的目的在于,克服已有的技术局限,利用光纤布拉格光栅组建振动传感系统,提供了一种基于光纤光栅振动传感系统的直升机桨叶质量不平衡检测方法,该方法能够有效地检测到直升机桨叶的质量不平衡以及质量不平衡的程度。
本发明采用的技术方案为:一种基于光纤光栅振动传感系统的直升机桨叶质量不平衡检测方法,将光纤光栅传感器粘贴在直升机桨叶的表面和桨叶底座上,直升机桨叶表面的首个光纤光栅传感器连接至光纤滑环的动端,桨叶底座的首个光纤光栅传感器连接至光纤滑环的静端,桨叶底座的最后一个光纤光栅传感器连接至光纤光栅振动解调系统;在桨叶转动时感受桨叶表面的应变和桨叶底座的振动并采集数据,处理数据得到振动的时域参数:振幅、均值、方差、均方根值、峰值、峰值因子、峭度系数;然后通过快速傅里叶变换(FFT)计算频域参数:基频、重心频率、均方频率、均方根频率、频率方差;将时域、频域参数与正常桨叶的振动数据对比,判断桨叶质量是否平衡以及不平衡程度。
进一步的,光纤光栅传感器的粘贴方式采取栅区两侧用3M胶带粘贴,栅区用α-氰基丙烯酸乙酯胶粘贴,各个光纤光栅传感器熔接为一串,直升机桨叶表面的首个光纤光栅传感器连接至光纤滑环的动端,桨叶底座的首个光纤光栅传感器连接至光纤滑环的静端,桨叶底座的最后一个光纤光栅传感器连接至光纤光栅振动解调系统。
进一步的,采集到的原始数据为光纤光栅的中心波长,先将其转换为应变值,得到振动幅值的时域数据,计算均值、方差、均方根值、峰值、峰值因子、峭度系数时域参数,然后通过快速傅里叶变换(FFT)转换成频域数据,计算重心频率、均方频率、均方根频率、频率方差频域参数,通过将以上数据与质量平衡的桨叶振动数据相对比来判断所测桨叶是否质量不平衡以及质量不平衡的程度。
其中,由振动幅值计算得到均值、方差、均方根值、峰值、峰值因子、峭度系数时域参数的方法为时域分析方法。由时域信号进行快速傅里叶变换并计算得到基频、重心频率、均方频率、均方根频率、频率方差频域参数的方法为频域分析方法。
本发明与现有技术相比的优点在于:可以准确测量直升机桨叶旋转时的数据,并可保存下来量化分析;多个特征参数对比得出结论,对桨叶质量不平衡的判断成功率高;数据客观,检测过程重复性好,且操作简单、安全;目前检测直升机桨叶质量不平衡的人工纸筒式标杆法和频闪仪靶标法都依赖于人工的技术,判断不够客观且无法记录数据,而本发明利用光纤光栅振动传感系统进行测量和记录数据,能够避免人工观测的主观性,而且更加准确、可靠,此为创新性之一;本发明在布置传感器时,根据不同部位的不同数据特征在桨叶表面和桨叶底座表面布置了传感器,通过光纤滑环连为一串,在全面测量数据的同时又解决了桨叶转动带来的连接线扰动问题,此为创新性之二;本发明在数据处理时,从时域参数和频域参数两方面进行计算,通过多个参数的对比来判断桨叶质量的不平衡,成功率高,而且能判断出桨叶质量不平衡的程度,此为创新性之三。
附图说明
图1为本专利的方法示意图;
图2为光纤光栅振动传感系统的结构连接示意图;
图3为光纤光栅传感器在直升机桨叶表面粘贴的示意图。
具体实施方式
下面结合附图对本发明的具体实施方式进行描述,以便更好地理解本发明。需要特别提醒注意的是,在以下的描述中,当采用已知功能和设计的详细描述也许会淡化本发明的主要内容时,这些描述在这里将被忽略。
图1为本专利的整体思路与检测流程图,1为在直升机桨叶表面上粘贴光纤光栅传感器,2为在桨叶底座表面上粘贴光纤光栅传感器,3是直升机开大车状态下开机检测,4是选定合适的采样频率来记录数据,5是对采集到的原始数据进行时域分析,包括求取峰值、峰值因子、峭度系数、均方根值等参数,6是对采集到的数据进行频域分析,包括快速傅里叶变换(FFT)、求取基频、重心频率、均方频率等参数,7是将求取的时域、频域参数与正常桨叶运转的相应参数进行对比,判断所测桨叶是否存在质量不平衡问题以及质量不平衡的严重程度。
图2为光纤光栅振动传感系统的结构图,1为粘贴在直升机桨叶表面的光纤光栅传感器,这些传感器通过熔接接为一串,通过光纤滑环的转动端2连接到光纤滑环3,粘贴在桨叶底座表面上的光纤光栅传感器5,也是通过熔接接为一串,一端连接到光纤滑环的固定端4上,另一端连接到光纤光栅振动解调仪6上,6连接到PC 7上进行数据存储与后期处理。这套系统具有体积小巧,灵敏度高,抗电磁干扰等优点,尤其适合在直升机开大车时进行检测,而且通过创造性地引入光纤滑环,使得在直升机开大车状态下,桨叶的转动不会造成线缆扰动和弯曲,从而降低了发生危险的概率并减小了信号的衰减。
图3为光纤光栅传感器在直升机桨叶表面上的布置粘贴图,在粘贴光纤光栅传感器时,采取的方式为光纤光栅的栅区两侧用3M胶带粘贴,栅区用α-氰基丙烯酸乙酯粘贴,在粘贴时预加一部分应力在光栅两端,这样可以保证光纤光栅传感器能紧贴材料表面,从而充分传递材料本身的应变;光纤光栅传感器布置在桨叶的叶尖处,一个横向布置,一个纵向布置,这样可以感受不同方向的应变。
光纤光栅传感器是用光栅栅距的变化来反映外界物理量变化的,当直升机桨叶转动时,粘贴在桨叶表面和桨叶底座表面的光纤光栅传感器感受到的振动,实质上是动态应变,光纤光栅传感器测振动的原理如下:
光纤布拉格光栅(FBG)的中心波长为:
λB=2neffΛ (1)
式中的neff即代表纤芯的有效折射率,Λ则代表光纤光栅的周期(栅距)。当一束宽带光入射进布拉格光栅时,符合上述光栅谐振条件的窄带频谱将被反射回来,当外界的物理量如压力、温度、声场作用在光纤光栅上时,有效折射率neff和光栅周期Λ都会受到影响而发生改变,引起后向反射光中心波长的漂移。于是,可以得到:
式中:△neff为光纤纤芯有效折射率变化,△Λ为光栅栅距变化。
其中光纤产生应变时的折射率变化为:
式中:ε是轴向应变,μ是泊松比。对于典型的石英光纤:neff=1.46,P11=0.12,P12=0.27,μ=0.16,则P=0.22。若设则式(2)可以写为:
式(4)就是FBG应变测量的一般公式,可以看出当FBG的材料一旦确定,其对应的传感特性基本上是与材料系数相关的常数,这也说明FBG具有很好的线性应变特性。
光纤光栅振动传感系统测得的原始数据为各光纤光栅传感器的中心波长,利用(4)式将其转化为应变,按照时间点依次作图,即可得到应变的时域图,同时得到应变的幅值x(t)。时域参数的计算公式如下:
均值一个周期内函数绝对值的平均值,
其中,T为振动信号的样本长度,x(t)为振动的时域信号;
方差σx 2:各个数据分别与其平均数之差的平方的和的平均数,
其中,N为离散化信号的样本点数,x(ti)为i点处的样本值;
峰值XP:信号波形中的最大幅值,
XP=max[|x(ti)|] (7)
峰值因子C:信号峰值与其有效值的比值,用来描述信号波形尖峰度的一个指标,
其中,Xrms为信号的有效值;
峭度系数K:用来表示样本的密度函数p(x)图形顶峰的凸平度,
其中,x~N(μ,σ2),N为样本长度,μ为样本均值,σ为样本标准差;
频域参数的计算公式如下:
重心频率FC:信号的频率谱中能量最集中的频率,
其中,fi为i时刻的功率谱所对应的频率值;pi为i时刻的功率谱幅值,N为时间长度;
频率方差VF:各个频率分别与中心频率作差的平方的和的平均数,
其中,fc为中心频率;
均方频率MSF:单位功率内的频率平方值,
其中:fi为i时刻的功率谱所对应的频率值;pi为i时刻的功率谱幅值,N为时间长度;
均方根频率RMSF:单位功率内的频率值,
尽管上面对本发明说明性的具体实施方式进行了描述,以便于本技术领的技术人员理解本发明,但应该清楚,本发明不限于具体实施方式的范围,对本技术领域的普通技术人员来讲,只要各种变化在所附的权利要求限定和确定的本发明的精神和范围内,这些变化是显而易见的,一切利用本发明构思的发明创造均在保护之列。
Claims (3)
1.一种基于光纤光栅振动传感系统的直升机桨叶质量不平衡检测方法,其特征在于:将光纤光栅传感器粘贴在直升机桨叶的表面和桨叶底座上,直升机桨叶表面的首个光纤光栅传感器连接至光纤滑环的动端,桨叶底座的首个光纤光栅传感器连接至光纤滑环的静端,桨叶底座的最后一个光纤光栅传感器连接至光纤光栅振动解调系统;在桨叶转动时感受桨叶表面的应变和桨叶底座的振动并采集数据,处理数据得到振动的时域参数:振幅、均值、方差、均方根值、峰值、峰值因子、峭度系数;然后通过快速傅里叶变换(FFT)计算频域参数:基频、重心频率、均方频率、均方根频率、频率方差;将时域、频域参数与正常桨叶的振动数据对比,判断桨叶质量是否平衡以及不平衡程度。
2.根据权利要求1所述的一种基于光纤光栅振动传感系统的直升机桨叶质量不平衡检测方法,其特征在于:光纤光栅传感器的粘贴方式采取栅区两侧用3M胶带粘贴,栅区用α-氰基丙烯酸乙酯胶粘贴,各个光纤光栅传感器熔接为一串,直升机桨叶表面的首个光纤光栅传感器连接至光纤滑环的动端,桨叶底座的首个光纤光栅传感器连接至光纤滑环的静端,桨叶底座的最后一个光纤光栅传感器连接至光纤光栅振动解调系统。
3.根据权利要求1所述的一种基于光纤光栅振动传感系统的直升机桨叶质量不平衡检测方法,其特征在于:采集到的原始数据为光纤光栅的中心波长,先将其转换为应变值,得到振动幅值的时域数据,计算均值、方差、均方根值、峰值、峰值因子、峭度系数时域参数,然后通过快速傅里叶变换(FFT)转换成频域数据,计算重心频率、均方频率、均方根频率、频率方差频域参数,通过将以上数据与质量平衡的桨叶振动数据相对比来判断所测桨叶是否质量不平衡以及质量不平衡的程度。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510502911.9A CN105021351B (zh) | 2015-08-16 | 2015-08-16 | 一种基于光纤光栅振动传感系统的直升机桨叶质量不平衡检测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510502911.9A CN105021351B (zh) | 2015-08-16 | 2015-08-16 | 一种基于光纤光栅振动传感系统的直升机桨叶质量不平衡检测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105021351A true CN105021351A (zh) | 2015-11-04 |
CN105021351B CN105021351B (zh) | 2017-09-26 |
Family
ID=54411488
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510502911.9A Active CN105021351B (zh) | 2015-08-16 | 2015-08-16 | 一种基于光纤光栅振动传感系统的直升机桨叶质量不平衡检测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105021351B (zh) |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107314840A (zh) * | 2017-07-26 | 2017-11-03 | 北京航空航天大学 | 一种基于光纤光栅实测应变的机翼剪力实时监测方法 |
CN107782438A (zh) * | 2017-11-28 | 2018-03-09 | 广西师范大学 | 基于脉冲编码和边沿滤波法的光纤光栅振动测量系统 |
CN109376449A (zh) * | 2018-11-09 | 2019-02-22 | 中国直升机设计研究所 | 直升机机体振动水平评价方法及装置 |
CN112683189A (zh) * | 2020-12-10 | 2021-04-20 | 武汉理工大学 | 基于fbg的船用螺旋桨叶片动应变测量装置 |
CN113915077A (zh) * | 2021-11-26 | 2022-01-11 | 龙源(北京)风电工程技术有限公司 | 光纤光栅式风电机组变桨法兰、叶片智能诊断方法及系统 |
CN116448319A (zh) * | 2023-06-16 | 2023-07-18 | 天津赛象科技股份有限公司 | 一种轮胎动平衡检测方法、装置、介质、设备及终端 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
DE19754321A1 (de) * | 1997-12-08 | 1999-06-10 | Schenck Rotec Gmbh | Verfahren zur Unwuchtbestimmung und Unwuchtmeßeinrichtung |
CN101285452A (zh) * | 2008-04-29 | 2008-10-15 | 唐德尧 | 一种主动控制风力发电机桨叶动平衡的方法 |
CN101706356A (zh) * | 2009-12-04 | 2010-05-12 | 北京英华达软件工程有限公司 | 一种自动显示汽轮发电机组振动信号特征的装置及方法 |
KR20110116459A (ko) * | 2010-04-19 | 2011-10-26 | 한국항공우주산업 주식회사 | 밸런스 웨이트를 갖는 항공기 로터용 회전 시험체 모듈 |
CN104075795A (zh) * | 2014-06-18 | 2014-10-01 | 北京金风科创风电设备有限公司 | 用于监测风力发电机组的叶轮振动状态的方法和系统 |
CN104215395A (zh) * | 2014-09-09 | 2014-12-17 | 中国石油大学(北京) | 一种转子不平衡故障的检测方法及装置 |
-
2015
- 2015-08-16 CN CN201510502911.9A patent/CN105021351B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
DE19754321A1 (de) * | 1997-12-08 | 1999-06-10 | Schenck Rotec Gmbh | Verfahren zur Unwuchtbestimmung und Unwuchtmeßeinrichtung |
CN101285452A (zh) * | 2008-04-29 | 2008-10-15 | 唐德尧 | 一种主动控制风力发电机桨叶动平衡的方法 |
CN101706356A (zh) * | 2009-12-04 | 2010-05-12 | 北京英华达软件工程有限公司 | 一种自动显示汽轮发电机组振动信号特征的装置及方法 |
KR20110116459A (ko) * | 2010-04-19 | 2011-10-26 | 한국항공우주산업 주식회사 | 밸런스 웨이트를 갖는 항공기 로터용 회전 시험체 모듈 |
CN104075795A (zh) * | 2014-06-18 | 2014-10-01 | 北京金风科创风电设备有限公司 | 用于监测风力发电机组的叶轮振动状态的方法和系统 |
CN104215395A (zh) * | 2014-09-09 | 2014-12-17 | 中国石油大学(北京) | 一种转子不平衡故障的检测方法及装置 |
Non-Patent Citations (1)
Title |
---|
张保钦等: "风机叶片故障预测的振动方法研究", 《电子测量与仪器学报》 * |
Cited By (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107314840A (zh) * | 2017-07-26 | 2017-11-03 | 北京航空航天大学 | 一种基于光纤光栅实测应变的机翼剪力实时监测方法 |
CN107314840B (zh) * | 2017-07-26 | 2019-06-25 | 北京航空航天大学 | 一种基于光纤光栅实测应变的机翼剪力实时监测方法 |
CN107782438A (zh) * | 2017-11-28 | 2018-03-09 | 广西师范大学 | 基于脉冲编码和边沿滤波法的光纤光栅振动测量系统 |
CN109376449A (zh) * | 2018-11-09 | 2019-02-22 | 中国直升机设计研究所 | 直升机机体振动水平评价方法及装置 |
CN109376449B (zh) * | 2018-11-09 | 2022-07-01 | 中国直升机设计研究所 | 直升机机体振动水平评价方法及装置 |
CN112683189A (zh) * | 2020-12-10 | 2021-04-20 | 武汉理工大学 | 基于fbg的船用螺旋桨叶片动应变测量装置 |
CN112683189B (zh) * | 2020-12-10 | 2022-07-08 | 武汉理工大学 | 基于fbg的船用螺旋桨叶片动应变测量装置 |
CN113915077A (zh) * | 2021-11-26 | 2022-01-11 | 龙源(北京)风电工程技术有限公司 | 光纤光栅式风电机组变桨法兰、叶片智能诊断方法及系统 |
CN116448319A (zh) * | 2023-06-16 | 2023-07-18 | 天津赛象科技股份有限公司 | 一种轮胎动平衡检测方法、装置、介质、设备及终端 |
CN116448319B (zh) * | 2023-06-16 | 2023-09-15 | 天津赛象科技股份有限公司 | 一种轮胎动平衡检测方法、装置、介质、设备及终端 |
Also Published As
Publication number | Publication date |
---|---|
CN105021351B (zh) | 2017-09-26 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105021351A (zh) | 一种基于光纤光栅振动传感系统的直升机桨叶质量不平衡检测方法 | |
CN107271090A (zh) | 一种基于光纤光栅的飞机机翼弯矩实时监测方法 | |
CN105115438B (zh) | 一种光纤传感系统温度补偿方法 | |
Frieden et al. | High-speed internal strain measurements in composite structures under dynamic load using embedded FBG sensors | |
CN101339093B (zh) | 光纤陀螺用光纤环质量的测量方法及其装置 | |
CN104007175B (zh) | 一种悬臂柔性梁多裂缝损伤识别装置及方法 | |
Dewynter-Marty et al. | Embedded fiber Bragg grating sensors for industrial composite cure monitoring | |
CN105784237A (zh) | 一种微推力测试系统及方法 | |
Jin et al. | A fibre-optic grating sensor for the study of flow-induced vibrations | |
CN103575664B (zh) | 基于集成式复合传感器的结构多尺度健康监测装置与方法 | |
CN107727483A (zh) | 一种基于光纤光栅用于地基原位测试的贯入剪切装置及方法 | |
CN106053010A (zh) | 多分量光纤天平及其测量方法 | |
CN108180976A (zh) | 基于光纤光栅的振动状况和摩擦系数同时测量的方法 | |
CN110108340A (zh) | 一种汽车动态称重装置 | |
CN206311247U (zh) | 一种基于分布式光纤感测的力和位移测量的传感器装置 | |
CN106124123B (zh) | 一种扭矩传感器定度实验装置 | |
CN206019923U (zh) | 一种扭矩传感器定度实验装置 | |
CN207798532U (zh) | 一种基于光纤光栅用于地基原位测试的贯入剪切装置 | |
CN108152130A (zh) | 三边法风电叶片静力测试方法 | |
CN109708586A (zh) | 一种光纤布拉格光栅应变传感器的封装方法 | |
CN107044828A (zh) | 一种光纤光栅结构变形检测装置、网络与方法 | |
CN104121983B (zh) | 基于双折射光纤环镜的振动动态监测系统及方法 | |
CN106441653B (zh) | 一种地质力学模型试验隧道明洞模型受力应变监测方法 | |
CN110424362A (zh) | 一种光纤式温度自补偿静力触探传感器 | |
CN108180977A (zh) | 基于光纤光栅的振动状况和摩擦系数同时测量的装置 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
TR01 | Transfer of patent right | ||
TR01 | Transfer of patent right |
Effective date of registration: 20210419 Address after: 261057 room 17, building 23, Wuzhou garden, 1166 Weifang Economic Development Zone, Shandong Province Patentee after: Shandong Shuangshi Security Information Technology Industry Research Institute Co., Ltd Address before: 100191 Haidian District, Xueyuan Road, No. 37, Patentee before: BEIHANG University |