Optimization analysis of modified axial force element of sting balance
-
摘要: 针对天平轴向力受法向力干扰较大的问题,提出了一种竖直梁轴向力元件的改进结构。通过理论分析和有限元仿真,获得了轴向力/法向力作用下的轴向力元件变形情况;对结构进行受力分析,得到法向力作用下的轴向力测量梁的变形特点及变形原因。提出了一种改进的轴向力元件测量梁结构,测量梁与主梁的连接夹角为不等于90°的优化夹角,该夹角通过仿真结果中的法向力对轴向力的干扰输出与轴向力有效输出之比和该夹角的线性函数关系获得。对改进的天平轴向力元件的有限元仿真结果表明:与传统轴向力元件相比,改进结构的轴向力有效输出仅减小2.77%,但法向力对轴向力的干扰输出减小了99.32%。改进的轴向力元件具有良好的抗法向力干扰的效果,适用于具有大升阻比特性的飞行器风洞模型气动力试验要求。Abstract: Aiming at the problem that the measuring of axial force is greatly disturbed by the normal force, a modified axial force element of sting balance is proposed. By means of theoretical analysis and finite element simulation, the deformation of axial force element under axial force and normal force is obtained. Meanwhile, mechanical analysis has been carried out to get the characteristics and reasons of the deformation of the measuring beam under normal force. The deflection angle between measuring beams and main beams of the modified structure is specific, and it is not 90°, which is obtained by linear fitting between the ratio of the interference of normal force on axial force to the output of axial force and the deflection angle. Simulation result of the modified structure shows that the axial force output is only reduced by 2.77%, moreover, the interference of normal force is reduced by 99.32%, compared to the conventional axial force element. The modified axial force element can obviously reduce normal force interference on normal force, so it can be used for wind tunnel aerodynamic test with large lift-to-drag ratio requirement.
-
Keywords:
- sting balance /
- axial force element /
- normal force interference /
- measuring beam /
- linear fitting
-
0 引言
FL-12风洞作为主力风洞承担着大量重要试验任务,其中包括进气道、吹气流动控制、引射短舱和TPS等多项试验,这些试验任务都需要供气控制系统输送稳定、可控的气流到试验模型,并要达到优于±3g/s的绝对控制精度,同时要求流量调节范围大、稳定时间短,因此实现流量宽范围快速精确控制是供气试验最为基础和关键的环节[1-3]。
目前模拟调节阀的高精度阀门定位器死区0.1%,综合误差0.5%[4],达不到试验要求的流量控制精度,而且在实际工况中其工作流量特性随着试验模型的变化而变化,同时存在死区、调节滞后严重、调节范围不宽、密封件易磨损等问题,因此无法满足当前供气试验要求。
与模拟调节阀相比,数字阀具有可调范围宽、控制精度高、响应速度快、结构简单、鲁棒性强、可靠性高等特点[5-7],目前数字阀主要分为PWM(脉冲宽度调制)数字阀和PCM(脉冲编码调制)数字阀[8-9]。PWM数字阀也称高速开关式数字阀,通过调节PWM信号的占空比来控制阀的平均流量,但是数字阀的流量和得失电时间存在一定的矛盾。PCM数字阀是将开关阀按照一定规律的排列组合(二进制、斐波那契数列等)实现对流量的控制,但是离散化的形式会导致流量输出的不连续性,同时小口径喷嘴的加工难度较高,进退阶时波动较大,频繁切换阀会造成冲击和噪声,缩短低位开关阀的寿命。
为解决PWM数字阀流量和得失电时间的矛盾和PCM数字阀的流量不连续性和低位开关阀的寿命问题,本文设计了PCM+PWM数字阀。该数字阀通过PCM数字阀提供基准流量,PWM数字阀在基准流量的基础上再进行调节,其中PCM数字阀采用多路文丘里喷嘴与开关阀按照二进制规律排列组合而成,PWM数字阀由单路小口径文丘里喷嘴与高频开关阀组成。
1 数字阀的设计指标和设计方案
1.1 设计指标
依据试验任务需求,FL-12风洞供气控制系统中数字阀设计的主要技术指标:(1) 阀前压力范围:6.0~8.0MPa;(2) 阀后压力范围:0.1~6.0MPa;(3) 流量控制范围:0.1~6.0kg/s;(4) 绝对控制精度:≤±3g/s;(5) 稳定控制时间:≤30s。
1.2 设计方案
本文中设计的数字阀由13位半PCM数字阀+ PWM数字阀,共计16路喷嘴和16路开关阀组成。图 1为数字阀支路示意图,气流依次通过喷嘴、高压金属软管、电磁阀后流出,系统先通过开关不同位的电磁阀来控制每路管道的通断以达到组合出不同流量的目的,当流量输出值与给定值的差值进入预先设定的误差带内,则停止动作PCM数字阀,由PWM数字阀通过脉宽调制控制。在实际流量控制中,PCM数字阀起到量程调节作用,PWM数字阀起到精细调节作用。
2 数字阀设计的关键技术
2.1 喷嘴的选择与设计
为使喷嘴流量特性不受后端试验模型状态变化及背压影响,采用了临界流喷嘴。临界流喷嘴主要包括临界流标准喷嘴和临界流文丘里喷嘴,其中临界流标准喷嘴结构如图 2所示。如果保持喷嘴入口压力p1和温度T1不变,使其出口压力p2逐渐减小,则通过喷嘴的气体流量qm将逐渐增加。当出口压力下降至pc时,通过喷嘴的流量将达到最大值,气体流速为当地声速,再继续降低出口压力,通过喷嘴的流量将不再增加,流速也保持声速不变。将喷嘴出口的流速达到声速的压力pc称为临界压力,pc/p1称为临界压力比,通过喷嘴的流量(qm)max称为临界流量[10]。
对于空气(γ=1.4),可以通过公式1计算得到临界压力比为0.528。由于数字阀前压力一般稳定在8MPa,因此要保证“恒流”特性不变,出口压力p2必须低于临界压力4.224MPa,显然无法满足最大出口压力为6MPa的设计要求,而且压力损耗太大。
(1) 为使出口压力得到恢复,采用喉道后面带扩压管的临界流文丘里喷嘴,其出口压力高低与其扩散段尺寸有关[11]。根据JJG620-2008《临界流文丘里喷嘴》,可以计算出最大背压比p2/p1与扩散段面积比A2/A*的对应参数表(见表 1),其中p1为入口压力,p2为出口压力,A*为喉道面积,A2为扩散段出口面积。
表 1 临界流文丘里喷嘴的最大允许背压比表(γ=1.4)Table 1 Maximum allowable back pressure ratio of the critical flow Venturi-type nozzleA2/A* p2/p1 A2/A* p2/p1 1.1 0.6863 1.9 0.8495 1.2 0.7376 2.0 0.8554 1.3 0.7701 2.5 0.8743 1.4 0.7930 3.0 0.8842 1.5 0.8101 3.5 0.8900 1.6 0.8233 4.0 0.8937 1.7 0.8339 5.0 0.8981 1.8 0.8424 10.0 0.9038 依据表 1,选择喷嘴扩散段面积比为4.0,在理想条件下,当出口压力低于0.8937倍入口压力,即可达到临界状态,喷嘴喉道流速为声速,流量由入口压力、入口温度与喉道面积唯一确定。依据圆环形喉道临界流文丘里喷嘴设计准则[12],喷嘴的结构设计如图 3所示。
2.2 喉道面积的分配
PCM数字阀的喷嘴喉道面积按照二进制规律排列,通过这些不同的喉道可以组合不同喉道面积,组合喉道面积范围在所有喉道面积之和与最小喉道面积之间,最小分辨率是最小喉道面积,因此设计足够小的喉道面积就可以解决定位精度的问题,设计足够多的位数就可以解决流量调节范围的问题。
PCM数字阀的设计工况是入口压力为8.0MPa,入口温度为20℃,流量为6.0kg/s,在最大允许背压比范围内,根据公式(2) 可得PCM数字阀需要的喉道总面积为3.176cm2。
(2) 式中:p1为喷嘴入口压力(Pa);T1为喷嘴入口温度(K);A*为喉道面积(m2)。
经计算,将喉道总面积按照13位半二进制分配,各喷嘴喉道面积分配值及工作流量分配值如表 2所示。
表 2 PCM数字阀各喷嘴喉道面积及对应流量计算表Table 2 Nozzle throat area and calculated flow value of the PCM digital valveBit位 临界流量/(kg·s-1) 喷嘴喉道面积/cm2 喉道通径/mm 1 0.00059332 0.00031416 0.2000 2 0.0011866 0.00062832 0.2828 3 0.0023733 0.0012566 0.4000 4 0.0047466 0.0025133 0.5657 5 0.0094932 0.0050265 0.8000 6 0.018986 0.010053 1.1314 7 0.037973 0.020106 1.6000 8 0.075945 0.040212 2.2627 9 0.15189 0.080425 3.2000 10 0.30378 0.16085 4.5255 11 0.60756 0.3217 6.4000 12 1.2151 0.6434 9.0510 13 2.4303 1.2868 12.8000 13半 1.2151 0.6434 9.0510 合计 6.075 3.2170 实际设计PCM数字阀总喉道面积为3.217cm2,在不考虑喷嘴加工误差的情况下,其喉道面积分辨率(定位精度)为1/10239,在设计工况下对应的最小流量为0.6g/s。
PWM喷嘴喉道面积的选取需依据试验流量精度要求(±3g/s),并充分考虑裕量,预设误差带为±6g/s,可在表 2中选取喉道口径0.8mm的喷嘴作为PWM数字阀的喷嘴。
2.3 电磁阀参数计算与型号选择
流量系数CV是电磁阀主要的计算参数,表征为阀门流通能力大小,根据设计工况,并按照公式(3) 进行计算,结果如表 3所示。
(3) 表 3 电磁阀流量系数计算与型号选择Table 3 Magnetic valve flow coefficient calculation and model selectionBit位 计算流量系数/(m3·h-1) 电磁阀流量系数/(m3·h-1) 电磁阀阀座口径/mm 电磁阀型号 PWM 0.02892 0.05 8.0 2/943 1 0.00195 0.13 2.0 A5232 2 0.00390 0.13 2.0 A5232 3 0.00781 0.13 2.0 A5232 4 0.01563 0.13 2.0 A5232 5 0.03125 2.0 8.0 B4623 6 0.0625 2.0 8.0 B4623 7 0.12501 2.0 8.0 B4623 8 0.25003 2.0 8.0 B4623 9 0.50005 2.0 8.0 B4623 10 1.00010 2.0 8.0 B4623 11 2.00021 8.0 25.0 A4625 12 4.00042 8.0 25.0 A4625 13 7.99895 20.0 32.0 A4827 13半 4.00042 8.0 25.0 A4625 式中:Q为标准状态下最大流量,Nm3/h;G为空气比重,为1.1kg/Nm3; T为气体温度, ℃; p1为入口压力, kg/cm2; p2为出口压力, kg/cm2; Δp=p1-p2, kg/cm2。
PCM数字阀的开关阀采用德国GSR二位二通系列的电磁阀, PWM数字阀的高频开关阀采用德国GSR 2/943,该阀动作频率可达1kHz,具有脉宽调制功能,线性流量控制。为使各喷嘴入口压力、温度等条件一致,入口压力稳定,出口压力满足喷嘴喉道气流为声速流要求,电磁阀的流量系数尽量选大,以减小工作压力降,同时需要一定容积的阀体/管道作为稳定段,增加喷嘴入口压力的稳定性[13]。
由于电磁阀不同口径的开关时间存在一定差异,其口径不大于25mm时阀门开关时间约为0.1s,而口径为32mm时阀门开关时间长达0.44s。为了减少各阀门在时间上的差异,增加阀门切换过程中的稳定性,电磁阀的口径都选用25mm以内的。口径32mm的电磁阀则由2个口径25mm的电磁阀组合替代,因此13位半的PCM数字阀具有15个支路,再加上PWM数字阀,共计16个支路。
3 流量控制策略
3.1 供气控制系统
供气控制系统的高压气源来自12~22MPa的高压储气罐,经减压阀降低至稳定的8MPa,再通过数字阀输出试验需要的流量值。在数字阀后端有高精度的文丘里流量计实时监测和反馈输出的流量值,从而组成闭环控制系统,如图 4所示。
3.2 流量调节运算方法
供气控制系统采用了微分前置非线性PID流量调节运算方法,其控制量通过数字阀完成,实现主管路流量的控制。该运算方法首先根据主管路流量给定与流量反馈完成PID运算,然后通过PID运算结果和数字阀气流入口总压总温,计算数字阀阀位偏差,将此偏差作为数字阀阀位补偿值重新调整阀位,从而调整主管路流量输出值,当流量输出值与给定值的差值进入预先设定的误差带内(±6g/s),则停止动作PCM数字阀,由PWM数字阀通过脉宽调制控制。该运算的微分前置环节可对供气管道及管道设备(数字阀后的供气管道、空气换热器、空气过滤器等)容积造成的流量滞后进行补偿,缩短调节过渡过程时间,增加系统稳定性;比例和积分环节的运算律按照大误差小比例、小积分增益,以及小误差大比例、大积分增益的控制策略构建。
4 调试与应用
4.1 调试
FL-12风洞供气控制系统的数字阀安装示意图如图 5所示。
在调试中,需要在不同流量、不同试验模型状态的情况下,多次调整PID运算的比例增益、积分增益和调节微分增益等参数,使流量控制达到最佳状态。表 4给出了流量闭环控制的测试结果,当给定流量1.2、2.4、3.55和4.67kg/s时,系统流量控制的最大绝对误差分别为1.9、2.2、2.5和2.6g/s,最大相对误差分别为0.15%、0.09%、0.07%、和0.06%。调试数据表明,该套数字阀的流量控制精度能够满足设计指标±3g/s,同时能够实现宽范围的流量调节,且效果平稳,重复性好,可靠性高。
表 4 流量闭环控制的测试结果(单位:kg/s)Table 4 Test results of the flow closed loop control给定流量 1.2 2.4 3.55 4.67 均值流量 1.2012 2.3984 3.5518 4.6688 最大流量 1.2019 2.4022 3.5525 4.6720 最小流量 1.1996 2.3979 3.5480 4.6674 极差 0.0023 0.0043 0.0045 0.0046 最大绝对误差 0.0019 0.0022 0.0025 0.0026 最大相对误差 0.15% 0.09% 0.07% 0.05% 4.2 应用
数字阀在供气控制系统中调试完成后,开展了数十项供气试验。其中以某型飞机进气道试验为例,本次试验的供气流量为1.0和3.2kg/s,流量控制调节过程曲线和稳态过程曲线分别如图 6和7所示。
当供气流量为1.0kg/s时,流量稳定时间为15s,最大绝对误差为1.8g/s,最大相对误差为0.15%;当供气流量为3.2kg/s时,流量稳定时间为16s,最大绝对误差为2.7g/s,最大相对误差为0.08%。
5 结论
大量的供气试验结果表明,采用13位半PCM数字阀和PWM数字阀组合的数字阀在FL-12风洞供气控制系统中的应用是成功的,全面达到了设计指标。其控制精度高、响应时间短、调节范围宽和可靠稳定等特点,为风洞供气试验提供了强有力的支撑。
-
表 1 六分量载荷作用下轴向力元件测量梁应变
Table 1 Strain of axial force measuring element under six loads
应变计编号 A=200 N N=2000 N C=200 N M=50 N·m Nb=20 N·m L=20 N·m 应变输出/(10−6) 1 –218.10 –714.20 –11.80 –32.70 –205.90 835.00 2 234.80 –342.70 4.20 –41.40 70.20 130.80 3 235.20 –342.70 2.40 39.90 69.00 129.50 4 –217.10 –712.20 13.40 32.90 –210.30 829.50 5 –216.80 –711.30 12.00 –32.80 204.80 –830.60 6 234.70 –342.00 –3.90 –41.10 –70.10 –130.80 7 234.60 –344.00 –2.60 40.20 –68.00 –132.00 8 –217.40 –713.20 –13.60 32.80 210.70 –830.50 表 2 六分量载荷作用下轴向力元件输出
Table 2 Bridge outputs of axial force measuring element under six loads
桥路 A=200 N N=2000 N C=200 N M=50 N·m Nb=20 N·m L=20 N·m 电桥输出/(10−3) ΔUA1/U0 0.453 0.371 0.003 –0.001 0.278 –0.702 ΔUA2/U0 0.452 0.369 –0.002 0 –0.277 0.699 ΔUA 0.904 0.740 0.001 –0.001 0.001 –0.003 表 3 测量梁与支撑梁的结构尺寸与网格尺寸
Table 3 Structure and simulation grid size of measuring sheets and supporting sheets
测量梁 支撑梁 网格参数/mm 数量 y向高度/mm 宽度/mm 厚度/mm 数量 y向高度/mm 宽度/mm 厚度/mm 应变计仿真区域 测量梁和支撑梁 主梁 2 10 5 2 12 10 5 1 0.25 0.50 2.00 表 4 不同夹角的轴向力元件在6个分量载荷作用下的输出
Table 4 Axial force output with different deflection angle under six loads
夹角/(° ) 轴向力输出ΔU/(10−3) ΔUN/ΔUA A=200 N N=2000 N C=200 N M=50 N·m Nb=20 N·m L=20 N·m 85.0 0.855 −0.134 0 −0.001 0 0 −0.157 87.5 0.897 0.293 0 −0.003 0.001 −0.004 0.327 90.0 0.904 0.740 0 −0.001 0.001 −0.003 0.818 92.5 0.880 1.159 0.006 −0.003 −0.001 0.001 1.317 95.0 0.824 1.500 0 −0.001 0 0.001 1.820 表 5 优化夹角的轴向力元件在轴向力/法向力作用下的输出
Table 5 Axial force output with optimal deflection angle under axial force or normal force
夹角/(°) 轴向力输出ΔU/(10−3) ΔUN/ΔUA A=200 N N=2000 N 85.7 0.866 −0.017 −0.019 85.8 0.879 −0.005 −0.005 -
[1] 贺徳馨. 风洞天平[M]. 北京: 国防工业出版社, 2001年. [2] NOURI N M,MOSTAFAPOUR K,KAMRAN M,et al. Design methodology of a six-component balance for measuring forces and moments in water tunnel tests[J]. Measurement,2014,58:544-555. doi: 10.1016/j.measurement.2014.09.011
[3] RINGEL M,LEVIN D. Improved drag element for wind tunnel sting balances[J]. Journal of Aircraft,1985,22(11):979-982. doi: 10.2514/3.45234
[4] ZHANG S X,LI X,MA H Q,et al. Mechanical analysis of normal force interference on axial force measurement for internal sting balance[J]. Aerospace Science and Technology,2016,58:351-357. doi: 10.1016/j.ast.2016.08.028
[5] 陆文祥,彭超. 高精度大升阻比民机测力天平研制[J]. 流体力学实验与测量,1999,13(3):75-79,96. LU W X,PENG C. Development of the force measuring balances for civil planes with high accuracy and lift drag ratio[J]. Experiments and Measurements in Fluid Mechanics,1999,13(3):75-79,96.
[6] 史玉杰,黄勇,田正波. 小展弦比飞翼布局高速标模测力天平研制[J]. 实验流体力学,2015,29(5):50-54. DOI: 10.11729/syltlx20150015 SHI Y J,HUANG Y,TIAN Z B. Strain gauge balance development for force test on small aspect ratio flying wing high speed standard model[J]. Journal of Experiments in Fluid Mechanics,2015,29(5):50-54. doi: 10.11729/syltlx20150015
[7] 史玉杰,陈竹,田正波. 横Π型梁在风洞应变天平阻力结构上的应用[J]. 实验流体力学,2012,26(4):83-86. DOI: 10.3969/j.issn.1672-9897.2012.04.017 SHI Y J,CHEN Z,TIAN Z B. The application of thwart Π beam to axial force structure of wind tunnel strain gauge balance[J]. Journal of Experiments in Fluid Mechanics,2012,26(4):83-86. doi: 10.3969/j.issn.1672-9897.2012.04.017
[8] 熊琳,宫建,王金印,等. 小直径杆式应变天平轴向力元件设计问题的探讨[J]. 实验流体力学,2013,27(5):75-78. DOI: 10.3969/j.issn.1672-9897.2013.05.014 XIONG L,GONG J,WANG J Y,et al. Discussion about axial force element design of bending-beam strain-gauge balance with small diameter[J]. Journal of Experiments in Fluid Mechanics,2013,27(5):75-78. doi: 10.3969/j.issn.1672-9897.2013.05.014
[9] 闫万方,蒋坤,张江. 基于减小轴向力测量干扰的高精度测力天平研制[J]. 实验流体力学,2018,32(6):61-67. DOI: 10.11729/syltlx20180082 YAN W F,JIANG K,ZHANG J. Development of a six-component wind tunnel balance with lower interference on axial force measurement[J]. Journal of Experiments in Fluid Mechanics,2018,32(6):61-67. doi: 10.11729/syltlx20180082
[10] 顾岩,虞伟建. 有限元分析法在风洞天平中的应用[J]. 流体力学实验与测量,1999,13(4):82-86. GU Y,YU W J. The application of finite element analysis in wind-tunnel balance[J]. Experiments and Measurements in Fluid Mechanics,1999,13(4):82-86.
[11] 姚裕,张召明. 整体式盒式应变天平有限元设计[J]. 南京航空航天大学学报,2010,42(1):58-61. DOI: 10.3969/j.issn.1005-2615.2010.01.011 YAO Y,ZHANG Z M. Finite element design on integrated box balance[J]. Journal of Nanjing University of Aeronautics & Astronautics,2010,42(1):58-61. doi: 10.3969/j.issn.1005-2615.2010.01.011
[12] PARKER P, DELOACH R. Structural optimization of a force balance using a computational experiment design[R]. AIAA-2002-0540, 2002. doi: 10.2514/6.2002-540
[13] RHEW R. Strain-gage balance axial section design optimiza-tion using design of experiments[R]. AIAA 2005-7600, 2005. doi: 10.2514/6.2005-7600