为了提高飞机起降安全性,高升力装置在现代大型飞机设计中得以不断发展和广泛应用。在型号研制的初期阶段,一般通过借助计算流体力学(CFD)[1-8]或低速风洞测压试验数据[9-10],并借助有限元求解或多体动力学仿真等方法,确定襟翼结构、运动机构及其操纵系统的设计载荷。襟翼构型、缝道参数、扰流板打开对襟翼气动流场影响非常复杂,目前的CFD及风洞试验技术还难以精确捕捉和模拟流场环境及分布载荷情况[11]。适航标准要求对类似气动操纵面的特殊装置,其受载情况必须由试验数据确定[12]。
当前,国内外飞机制造商一般通过在襟翼翼面进行打孔测压等方式,进行襟翼气动压力分布试飞验证[13]。襟翼打孔测压的试飞改装,需新设计和生产特殊的试飞襟翼试验件,满足测压孔布置、测压管路敷设及结构强度相关要求。特殊构型的打孔测压襟翼试验件具备研制成本高、周期长和机上测压管路改装复杂等特点。结构强度和改装空间等诸多限制,致使测压孔有限,单块翼面一般包含2~3个测压剖面和40~50个测压点。基于有限的测压数据和仿真结果,通过经验插值和面积积分方式,确定的翼面总体压力分布及操纵系统关键构件操纵载荷仍存在较大的不确定性。
光纤应变传感测量,具有灵敏度高、响应速度快、机上改装简单、分布式、重量轻、抗电磁干扰强等优点,可以实现折射率、应变、曲率和温度等参数的探测[14]。可以较好适应襟翼运动机构复杂的振动、温度及电磁工作环境[15],开展飞行过程中的襟翼运动机构内力测定,进而完成襟翼翼面操纵系统驱动/约束载荷和翼面总载荷试飞测量与验证。
1 后缘襟翼运动机构某型号后缘襟翼分为内襟翼及外襟翼,内、外襟翼分别通过2套运动机构连接到机翼盒段、机身框等结构。后缘襟翼运动机构采用的机构布置如图 1所示,1号运动机构布置于内襟翼内侧,位于机翼根部整流罩内部。2~4号运动机构悬挂于机翼后梁及辅助梁之上,从翼根至翼尖分别布置2~4号运动机构。整流罩随动机构位于2~4号整流罩内部。襟翼驱动支座安装在机翼盒段辅助梁及后梁上,通过襟翼操纵系统驱动和控制襟翼的角度及运动。
![]() |
图 1 某后缘襟翼机构布置 Fig. 1 Arrangement of a certain tailing-edge flap mechanism |
在襟翼卡位位置上和运动过程中,后缘襟翼翼面总载荷P,与滑轨约束力F和驱动连杆约束力R构成一组平衡力系。4套运动机构的传感器布置相似,以3#运动机构为例进行传感器布置及组网说明,如图 2所示,在驱动连杆两个截面上布置8个应变传感器和1个温度传感器,襟翼滑轨布置2个应变传感器和1个温度传感器。所有传感器串联组网,形成一个测试通道,确定一端为激光输入端即主测试端口,另一端为激光输出端即备用测试端口。
![]() |
图 2 后缘襟翼机构传力路径 Fig. 2 Load path of trailing-edge flap mechanism |
通过测量试飞过程中驱动连杆的载荷应变,并扣除温度应变影响,确定襟翼摇臂作用于驱动连杆的载荷R。通过测量试飞过程中滑轨应变,并扣除温度应变影响,确定襟翼滑轮架作用于滑轨的载荷F。通过分别对内襟翼和外襟翼两套驱动机构,襟翼摇臂作用于驱动连杆的载荷R和襟翼滑轮架作用于滑轨的载荷F求和,即可获得内襟翼或外襟翼的总载荷P。襟翼摇臂作用于驱动连杆的载荷R,相对与操纵系统驱动点求矩,即为操纵系统所需克服的驱动/约束力矩T。
2.2 光纤传感测量系统为了确保测量系统及其机上改装的高可靠性,采用预封装式布拉格光栅传感器(Fiber Brag Grating, FBG)[16]。其基本原理是将光纤特定位置纤芯,制成折射率周期分布的光栅区。将一宽带光源发出的光注入光纤后,特定波长(布拉格反射光)的光波在这个区域内将被反射,如图 3所示。
![]() |
图 3 光纤布拉格光栅应变传感原理 Fig. 3 Theory of grating strain sensor for fiber brag |
反射光中心波长信号跟光栅周期和纤芯的有效折射率有关,即
$ \lambda_{\mathrm{B}}=2 n_{\mathrm{eff}} \varLambda $ | (1) |
式中:neff为纤芯的有效折射率;Λ为布拉格光栅栅距;λB为反射光中心波长。
当光纤轴向受拉力时, 将使光纤产生拉伸应变, 其布拉格光栅栅距Λ亦相应改变, 从而导致反射光波长变化[17]。其波长改变量ΔλB与应变ε的关系为
$ \Delta \lambda_{\mathrm{B}} / \lambda_{\mathrm{B}}=\left(1-P_{\mathrm{e}}\right) \varepsilon $ | (2) |
式中:ε为光栅应变;Pe为有效光弹系数,
$ P_{{\rm e}}=n_{\mathrm{eff}}^{2}\left[P_{12}-k\left(P_{11}+P_{12}\right)\right] / 2 $ | (3) |
式中:k为泊松比;P11和P12为光弹系数。因此应变与波长解调精度及光纤布拉格光栅栅距的关系为
$ \varepsilon=\left(1-P_{\mathrm{e}}\right) \cdot \Delta \lambda_{\mathrm{B}} /\left(2 n_{\mathrm{eff}} \cdot \varLambda\right) $ | (4) |
当布拉格光栅栅距Λ越小和激光波长解调仪越灵敏时,应变精度越高。某型号光纤布拉格光栅传感器的应变测量精度为1 με,温度测量精度为0.1 ℃。
根据测量点空间位置,将多个具有不同栅距的布拉格光栅按照不同的间隔,定制在同一根光纤上。光栅栅距需根据测量应变值范围定制,确保其反射光谱波长范围互不重叠,以实现对待测结构的准分布式测量[18]。光纤传感测量系统如图 4所示。
![]() |
图 4 光纤传感测量系统 Fig. 4 Optic fiber sensor measure system |
4个运动机构分别采用一根测试光纤,降低测试光纤机上敷设难度和提高测试系统可靠性。测试总线为一根8芯集成光缆,其中4根纤芯分别与4套运动机构测试光纤的激光输入端及反射光对接,另外一根纤芯分别与4套运动机构测试光纤的输出端对接,以用于输出光数据采集,作为数据备份。
2.3 测试系统校准通过有限元仿真虚拟标定,将构件应变响应大、应力梯度小和无运动碰触位置,确定为应变测量点[19]。然后进行机构构件的温度-应变地面和载荷-应变标定试验,进而建立应变-载荷相关方程[20],具体的地面校准试验流程如图 5所示。
![]() |
图 5 地面校准试验流程 Fig. 5 Flow chart of ground calibration test |
将完成应变传感器和温度传感器改装的试验件,放入地面高低温环境试验箱,开展温度-应变地面标定试验。温度变化值Δt与[εj]之间可以表示为多元线性回归方程形式:
$ \boldsymbol{a}_{j} \Delta t=\boldsymbol{\varepsilon}_{j} $ | (5) |
式中:j=1, 2, …, n为应变传感器个数。通过计算,求解温度-应变的相关系数aj。
通过地面载荷应变标定试验,建立载荷应变相关方程。驱动连杆输出端承受任意方向力R,可以表示成图 2中参考坐标系下的fix、fiy和fiz(可忽略相对小量),fix、fiy和fiz分别驱动连杆输出端承受力R的x、y和z向分量,与应变εij之间可以表示为多元线性回归方程形式:
$ \left[\begin{array}{cccc} \varepsilon_{11} & \varepsilon_{12} & \cdots & \varepsilon_{1 n} \\ \varepsilon_{21} & \varepsilon_{22} & \cdots & \varepsilon_{2 n} \\ \vdots & \vdots & & \vdots \\ \varepsilon_{m 1} & \varepsilon_{m 2} & \cdots & \varepsilon_{m n} \end{array}\right]\left[\begin{array}{ccc} \beta_{11} & \beta_{12} & \beta_{13} \\ \beta_{21} & \beta_{22} & \beta_{23} \\ \vdots & \vdots & \vdots \\ \beta_{n 1} & \beta_{n 2} & \beta_{n 3} \end{array}\right]=\\ \left[\begin{array}{ccc} f_{1 x} & f_{1 y} & f_{1 z} \\ f_{2 x} & f_{2 y} & f_{2 z} \\ \vdots & \vdots & \vdots \\ f_{m x} & f_{m y} & f_{m z} \end{array}\right] $ | (6) |
式中:m为地面标定试验工况数;n为应变传感器个数。通过矩阵运算,求解相关系数β。
对于不好拆卸,进行单独的载荷应变地面物理标定试验的襟翼滑轨,通过有限元方法,在滑轨位置xi施加垂向力fiy和侧向力fiz(可忽略相对小量),与应变εij之间可以表示为多元线性回归方程形式:
$ \left[\begin{array}{cccc} \varepsilon_{11} & \varepsilon_{12} & \cdots & \varepsilon_{1 n} \\ \varepsilon_{21} & \varepsilon_{22} & \cdots & \varepsilon_{2 n} \\ \vdots & \vdots & & \vdots \\ \varepsilon_{m 1} & \varepsilon_{m 2} & \cdots & \varepsilon_{m n} \end{array}\right]\left[\begin{array}{cc} \beta_{11} & \beta_{12} \\ \beta_{21} & \beta_{22} \\ \vdots & \vdots \\ \beta_{n 1} & \beta_{n 2} \end{array}\right]=\\ \left[\begin{array}{cc} x_{1} f_{1 x} & x_{1} f_{1 z} \\ x_{2} f_{2 x} & x_{2} f_{2 z} \\ \vdots & \vdots \\ x_{m} f_{m x} & x_{m} f_{m z} \end{array}\right] $ | (7) |
应变响应传感器的优选,遵循应变响应大且线性度好的原则,剔除遵循应变响应小且线性度差的应变通道,以提高载荷识别精度[21]。某型号襟翼运动机构驱动连杆,应变-载荷相关方程的载荷识别误差如表 1所示,Fx为驱动连杆输出端承受任意方向R洽驱动连杆绕驱动端旋转的切向载荷;Fy为R绕驱动端旋转径向载荷。载荷识别误差最大不超过0.77%,小于2%的一般载荷试飞校准误差要求。
驱动连杆名称 | 校验载荷 | 识别载荷 | 识别误差% | |||
Fx/N | Fy/N | Fx/N | Fy/N | |||
1# | 5 000 | 3 000 | 5 023 | 2 941 | -0.18 | |
8 000 | 2 000 | 8 007 | 1 960 | -0.03 | ||
2# | 5 000 | 3 000 | 5 048 | 2 982 | 0.55 | |
8 000 | 2 000 | 8 036 | 1 960 | 0.31 | ||
3# | 5 000 | 4 000 | 4 964 | 3 966 | -0.77 | |
8 000 | 2 000 | 7 984 | 1 960 | -0.30 | ||
4# | 5 000 | 5 000 | 5 099 | 4 976 | 0.76 | |
8 000 | 3 000 | 8 050 | 2 996 | 0.53 |
在完成基于光纤传感的试飞测试系统改装和地面标定试验后,在特定飞行状态下(高度15 000 ft (1 ft=0.304 8 m)、飞行速度160 ktIAS~230 ktIAS(1 kt=1.852 km/h)、减速板打开),开展了襟翼操纵载荷试飞测试验证。在试飞过程中,以起飞地面的应变和温度为零值,测量和记录了应变ε′ij和温度变化Δt(应变和温度作为零值)。采用式(8)消除温度应变影响,获得载荷作用下的应变值。
$ \varepsilon_{i j}=\varepsilon_{i j}^{\prime}-\Delta t a_{j} $ | (8) |
然后,将驱动连杆对应的应变数据εij代入式(6),获得驱动连杆输出端承受任意方向力R。将襟翼滑轨对应的应变数据εij和通过翼面偏度及运动几何参数确定的xi,代入式(7),获得襟翼滑轨承受滑轮架作用的任意方向力F。电阻应变与光纤光栅传感对3#襟翼驱动连杆驱动载荷P的识别结果对比如图 6所示。
![]() |
图 6 电阻应变与光纤传感载荷识别结果对比 Fig. 6 Comparison of identifing loads between resistance strain and optical fiber sensor |
通过对比可以看出光纤光栅传感对3#襟翼驱动连杆驱动载荷峰值识别结果与电阻应变载荷峰值识别结果误差不大于5%。载荷识别结果的时间历程曲线形态及跟随性基本一致,表明光纤光栅传感的载荷识别结果有效性。
减速板角度对襟翼驱动连杆的驱动载荷具有直接的影响,打开相同的减速板角度,10°襟翼3#驱动连杆驱动载荷变化及其载荷峰值,大于19°襟翼3#驱动连杆驱动载荷变化及其载荷峰值。驱动连杆驱动载荷R相对于驱动点取矩即为襟翼翼面操纵系统驱动/约束载荷试飞,基于光纤应变的襟翼驱动连杆驱动载荷识别结果如表 2所示。
襟翼角度/(°) 数据来源 |
减速板角度/(°) | 内襟翼驱动载荷矩/(N·m) | 外襟翼驱动载荷矩/(N·m) | |||
1# | 2# | 3# | 4# | |||
10实测值 | 0 | 550 | 1 600 | 2 000 | 1 100 | |
10 | 540 | 1 650 | 3 000 | 1 500 | ||
(-2%) | (3%) | (50%) | (36%) | |||
25 | 530 | 1 800 | 3 800 | 1 800 | ||
(-4%) | (13%) | (90%) | (64%) | |||
10理论值 | 25 | 600 | 2 050 | 4 200 | 2 050 | |
20实测值 | 0 | 850 | 1 500 | 2 200 | 1 400 | |
10 | 845 | 1 600 | 2 900 | 1 600 | ||
(-1%) | (7%) | (32%) | (14%) | |||
25 | 880 | 1 800 | 3 100 | 1 800 | ||
(4%) | (20%) | (41%) | (29%) | |||
20理论值 | 25 | 900 | 1 800 | 3 400 | 1 850 | |
注:()内的值为相对误差。 |
通过对外襟翼2套驱动机构的驱动连杆的载荷R和襟翼滑轮架作用于滑轨的载荷F求和,获得外襟翼的总载荷P如表 3所示。
数据来源 | 减速板角度/(°) | X坐标/mm | Y坐标/mm | 总载荷/N |
理论值 | 0 | 24 248 | 9 465 | 49 000 |
实测值 | 0 | 24 150 | 9 331 | 46 000 |
理论值 | 25 | 24 173 | 9 357 | 53 000 |
实测值 | 25 | 24 081 | 9 208 | 51 000 |
总体上来看,试飞实测结果验证了理论计算结果的合理保守性和可靠性。试飞实测结果表明,当襟翼10°时,加速板打开情况下的操纵载荷影响(最大增幅90%)大于襟翼20°的影响(最大增幅32%)。位于外襟翼上表面的减速板打开,不仅影响外襟翼3#和4#作动器操纵载荷,同时也影响内襟翼靠外的2#作动器操纵载荷(增幅约20%)。因此,相较于减速板关闭时的机翼构型,减速板打开情况下的襟翼作动器操纵载荷具有显著的增加量,必须谨慎分析相关影响。
4 结论通过对比,可以看出光纤光栅传感载荷识别结果与电阻应变载荷识别结果误差较小,满足工程精度需求。两者之间载荷识别结果的时间历程曲线基本一致,表明了光纤光栅传感载荷识别技术的有效性。光纤光栅传感技术,为某型号襟翼操纵载荷飞行试验提供了一种相较于打孔测压及电阻应变襟翼载荷试飞,更高精度和更低成本的技术途径和方法。其独特的灵敏度高、响应速度快、机上改装简单、分布式、重量轻、抗电磁干扰强等优点,对民用飞机后缘襟翼运动机构复杂的振动、温度及电磁工作环境具有很好的适应性。可以预见,在其他复杂受载环境的民用飞机运动机构载荷试飞测定和结构健康检查等领域,光纤光栅应变传感技术的应用前景广阔,值得深入探讨和分析研究。
[1] | LU W S, TIAN Y, LIU P Q. Aerodynamic optimization and mechanism design of fiexible variable camber trailing-edge flap[J]. Chinese Journal of Aeronautics, 2017, 30(3): 998-1003. |
Click to display the text | |
[2] |
王运涛, 李松, 孟德虹, 等. 不同襟翼偏角梯形翼构型气动特性数值模拟[J]. 航空学报, 2015, 36(6): 1823-1829. WANG Y T, LI S, MENG D H, et al. Numerical simulation of aerodynamic characteristics of trapezoidal wing configuration at different flap angles[J]. Acta Aeronautica et Astronautica Sinica, 2015, 36(6): 1823-1829. (in Chinese) |
Cited By in Cnki | Click to display the text | |
[3] |
米百刚, 詹浩, 朱军. 二维干净翼型、增升构型地面效应的数值模拟研究[J]. 应用力学学报, 2013, 30(6): 822-827. MI B G, ZHAN H, ZHU J. Numerical simulation on aerodynamic characteristic of clean and multi-element airfoils in ground effect[J]. Chinese Journal of Applied Mechanics, 2013, 30(6): 822-827. (in Chinese) |
Cited By in Cnki (1) | Click to display the text | |
[4] |
秦绪国, 刘沛清, 屈秋林, 等. 多段翼型地面效应数值模拟与分析[J]. 航空动力学报, 2011, 26(4): 890-896. QIN X G, LIU P Q, QU Q L, et al. Numerical simulation on stall of wing in ground effect[J]. Journal of Aerospace Power, 2011, 26(4): 890-896. (in Chinese) |
Cited By in Cnki (5) | Click to display the text | |
[5] |
秦绪国, 刘沛清, 屈秋林, 等. 三维多段机翼地面效应数值模拟[J]. 航空学报, 2011, 32(2): 257-264. QIN X G, LIU P Q, QU Q L, et al. Numercial simulation on 3Dmulti-element wings in ground effect[J]. Acta Aeronautica et Astronautica Sinica, 2011, 32(2): 257-264. (in Chinese) |
Cited By in Cnki | Click to display the text | |
[6] |
李孝伟, 乔志德. 多段翼型大迎角下主翼、襟翼上的分离流及缝道流动[J]. 航空学报, 1999, 20(1): 55-57. LI X W, QIAO Z D. Seperated flows over main element and flap of multielement airfoil at high angles of attack[J]. Acta Aeronautica et Astronautica Sinica, 1999, 20(1): 55-57. (in Chinese) |
Cited By in Cnki (20) | Click to display the text | |
[7] |
朱自强, 陈迎春, 吴宗成, 等. 高升力系统外形的数值模拟计算[J]. 航空学报, 2005, 26(6): 257-262. ZHU Z Q, CHEN Y C, WU Z C, et al. Numerical simulation of high lift system configuration[J]. Acta Aeronautica et Astronautica Sinica, 2005, 26(6): 257-262. (in Chinese) |
Cited By in Cnki (71) | Click to display the text | |
[8] |
李孝伟, 乔志德. 带襟、副翼三维机翼粘性绕流计算[J]. 航空学报, 2001, 22(3): 55-57. LI X W, QIAO Z D. Navier-stokes computations for the wings with flaps or aileron[J]. Acta Aeronautica et Astronautica Sinica, 2001, 22(3): 55-57. (in Chinese) |
Cited By in Cnki | Click to display the text | |
[9] |
孙智伟, 白俊强, 高正红, 等. 现代超临界翼型设计及其风洞试验[J]. 航空学报, 2015, 36(3): 804-818. SUN Z W, QIAO J Q, GAO Z H, et al. Design and wind tunnel test investigation of the modern supercritical airfoil[J]. Acta Aeronautica et Astronautica Sinica, 2015, 36(3): 804-818. (in Chinese) |
Cited By in Cnki (8) | Click to display the text | |
[10] |
肖涛, 代钦. Gurney襟翼对机翼地面效应气动特性和流动结构的影响实验研究[J]. 空气动力学学报, 2013, 31(5): 572-578. XIAO T, DAI Q. The experimental study of aerodynamics and flow structures of a wing with Gurney flap in ground effect[J]. Acta Aerodynamica Sinica, 2013, 31(5): 572-578. (in Chinese) |
Cited By in Cnki (6) | Click to display the text | |
[11] |
熊磊, 刘洋, 毛俊. 大型运输类飞机后缘襟翼气动载荷特性分析[J]. 空气动力学学报, 2017, 35(3): 400-403. XIONG L, LIU Y, MAO J. Analysis of flight load characteristics of transport aircraft flap[J]. Acta Aerodynamica Sinica, 2017, 35(3): 400-403. (in Chinese) |
Cited By in Cnki | Click to display the text | |
[12] |
中国民用航空局.中国民用航空规章第25部: 运输类飞机适航标准: CCAR25-R4[S].北京: 中国民用航空局, 2011: 43. Civil Aviation Administration of China China. Civil Aviation Regulations Part 25 Airworthiness Standard: Transport Category Airplanes: CCAR25-R4[S]. Beijing: Civil Aviation Authority of China, 2011: 43(in Chinese). |
[13] |
王兆东. 民机飞行载荷的试飞验证[J]. 民用飞机设计与研究, 2011(2): 16-20. WANG Z D. The test flight demonstration of civil airplane flight load[J]. Civil Aircraft Design and Research, 2011(2): 16-20. (in Chinese) |
Cited By in Cnki (7) | Click to display the text | |
[14] |
刘冲冲, 邹翔, 周正仙. 高精度动态光纤应变传感研究[J]. 光子学报, 2016, 45(11): 1106007-1. LIU C C, ZOU X, ZHOU Z X. Study of high-precision dynamic optical fiber strain sensing[J]. Acta Photonica Sinica, 2016, 45(11): 1106007-1. (in Chinese) |
Cited By in Cnki (1) | Click to display the text | |
[15] |
卢卫涛, 张拴民, 陈安锋, 等. 高灵敏度光纤应变传感器[J]. 应用光学, 2017, 38(5): 848-851. LU W T, ZHANG S M, CHEN A F, et al. High sensitive fiber strain sensor[J]. Journal of Applied Optic, 2017, 38(5): 848-851. (in Chinese) |
Cited By in Cnki | Click to display the text | |
[16] |
徐先东.光纤光栅温敏与温度补偿式封装技术的研究[D].武汉: 武汉理工大学, 2003: 12-23. XU X D. Study of fiber bragg grating packaging technology on temperature sensitivity and temperature compensation[D]. Wuhan: Wuhan University of Technology, 2003: 12-23(in Chinese). |
Cited By in Cnki | Click to display the text | |
[17] |
张春光, 黄民双, 陶宝祺. 光纤应变传感技术研究[J]. 连云港化工高专学报, 1998, 11(4): 7-10. ZHANG C G, HUANG M S, TAO B Q. Study of optical fiber strain technology[J]. Journal of Lianyungang College of Chemical Technology, 1998, 11(4): 7-10. (in Chinese) |
Cited By in Cnki (2) | Click to display the text | |
[18] |
陈池, 周宜红.大量程分布式光纤传感技术研究及工程应用[D].武汉: 武汉大学, 2013: 5-10. CHEN C, ZHOU Y H. The sensing technology and application of broad dynamic range distributed optical fiber sensor[D]. Wuhan: Wuhan University, 2013: 5-10(in Chinese). |
Cited By in Cnki | Click to display the text | |
[19] |
杨全伟, 舒成辉, 赵华. 飞机机翼载荷校准试验的有限元仿真方法[J]. 系统仿真学报, 2008, 20(23): 6569-6571. YANG Q W, SHU C H, ZHAO H. Technology of FEM simulation in aircraft wing loads calibration test[J]. Journal of System Simulation, 2008, 20(23): 6569-6571. (in Chinese) |
Cited By in Cnki (8) | Click to display the text | |
[20] |
范华飞, 李俊, 李志蕊. 短连杆飞行载荷实测技术[J]. 航空工程进展, 2015, 6(4): 447-451. FAN H F, LI J, LI Z R. Flight loads measurement technique for short rods[J]. Advances in Aeronautical Science and Eingineering, 2015, 6(4): 447-451. (in Chinese) |
Cited By in Cnki (2) | Click to display the text | |
[21] |
金秀芬, 余建虎, 于秀伟. 民用飞机载荷校准试验关键技术浅析[J]. 航空工程进展, 2017, 8(1): 105-108. JIN X F, YU J F, YU X W. Initial analysis on the key technologies of civil aircraft load calibration test[J]. Advances in Aeronautical Science and Engineering, 2017, 8(1): 105-108. (in Chinese) |
Cited By in Cnki | Click to display the text |