直升机的振动一直是非常突出的问题。旋翼作为直升机的主要升力面和操纵面，是直升机主要振源之一。直升机在前飞时，旋翼的转动使桨叶经历周期性的气流环境变化，弹性桨叶在周期气动激振力作用下产生动载荷；各片桨叶动载荷合成到桨毂上形成桨毂六力素。桨叶振动载荷可以看作是由基频为1，2，3 Ω…各阶谐波分量叠加组成，对于N片桨叶的旋翼，桨毂相当于“滤波器”，只允许频率为kN Ω及(kN±1)Ω的谐波载荷通过传递到机身结构上。因此，通过降低桨叶的振动水平能够降低桨叶根部的振动载荷，从而使得直升机机身结构的激励大幅降低。

为了降低直升机振动，直升机大部分采用传统的减振吸振装置来抑制旋翼/机体振动。现如今，通过被动的吸振器和优化的动力学设计，直升机的振动水平已经从20世纪50年代的0.3～0.5g降低到0.1g左右。虽然直升机振动水平得到一定的改善，但相对于固定翼飞行器，上述振动水平还是较大，与直升机民用市场的预期振动水平和NASA提出的0.02g的目标相差较大。因此，需要一个更高效的方法来控制振动。

近二十年来，直升机振动主动控制技术日益受到重视。直升机的主动控制是通过控制系统输入能量从而来降低振动载荷。主要控制方法[1]有高阶谐波控制(Higher-order harmonic control，HHC)、独立桨叶控制(Individual blade control，IBC)、主动后缘襟翼控制和结构响应主动控制(Active control of structure response，ACSR)等。

(1)高阶谐波控制(HHC)

高阶谐波控制采用闭环控制系统，通过监测直升机飞行时的振动反馈信号，以一定控制规律传递给与自动倾斜器不旋转环相连接的作动器产生附加的变距角，从而可以抵消振动响应中的部分或全部谐波分量。

(2)独立桨叶控制(IBC)

独立桨叶控制[2-3]属于高阶谐波控制的一个分支，其通过独立安装在桨叶拉杆上的作动器直接对每片桨叶单独进行高阶谐波控制(图1)。与传统的HHC技术相比，其更具灵活性并可以补偿桨叶间的差异。

(3)IBC高阶谐波输入

从田文镜对黄河的治理可以看出，田文镜能够在中年之后得以雍正帝的重用，主要在于他的大义为民。他的多数主张符合了雍正帝的用人标准：厉行新政，尽职尽责，刚正不阿，清廉节俭。他在河南的政绩确实十分显著，使得官场风气为之一新，雍正便赞称他为“巡抚中第一人”，如此之高的评价可见雍正帝对他的宠爱。

IBC的变距输入信号表达式

如果要追问是什么事物最深刻改变了现代学生的生存方式，又是什么事物集结了大学生的时代意识，那么，这个事物一定非网络技术莫属。网络让学生们找到了一条通向满足各种需求的阳关大道，极大迎合了年轻人寻找心灵家园的迫切感。手机APP等新兴网络技术的出现和应用也给大学生的人格培养造成了种种挑战。较多见的人格障碍有因刻意追求完美而带有不完善感的强迫型人格障碍、以多疑敏感为主要表现的偏执型人格障碍、以情绪低落、多愁善感、缺乏自信、对环境困难估计过高为主要表现的抑郁型人格障碍等。

ϑi,IBC=ϑncos(nψi-φn)

(1)

式中：ϑi,IBC为IBC的变距量，ψi为i片桨叶的方位角；则控制信号有3个参变量：幅值ϑn、谐波阶次n、相位φn；通过控制这3个变量选取合适的谐波分量加载，使得振动控制达到最优。

为考察过境免签政策实施的动态效果，本文引入多期DID模型，将政策虚拟变量与政策实施后的年份的虚拟变量构造交乘项，从而判断政策实施的阶段性效果。修正后的计量模型如下：

  

图1 独立桨叶控制原理图Fig.1 Schematic of IBC

1 旋翼载荷模型

1.1 旋翼气动力模型

定常来流下，翼型的迎角不随时间改变，升力系数随迎角是线性变化的。旋翼在周期变化的气动环境中，迎角随时间不断改变，翼型处于非定常来流状态下，所以气动力响应也是非定常的。

本文采用Leishman-Beddoes[4-6]非定常/动态失速模型，并引入动态入流模型，计算旋翼的气动响应。该气动模型对气动响应模拟主要分为以下3个部分：

(1)附着流状态下，气动力主要由环量响应和非环量响应构成；

(2)通过比较滞后法向力与临界法向力，判定翼型剖面气流是否发生后缘分离，计算分离点并重新校准气动力；

(3)当后缘分离下的滞后法向力大于临界法向力时，深度失速发生，前缘涡产生并严重影响气动力值。

通过离散时间，选定时间步长，由Duhamel积分推导气动力计算公式，通过编程进行迭代计算。

通过计算桨叶旋转状态下各阶固有模态生成MNF文件导入ADAMS，从而计算弹性桨叶的动力学方程。如图3所示，计算前进比μ=0.14时桨毂垂向载荷。

  

图2 L-B气动模型算例验证Fig.2 Validation of L-B aerodynamic model

⋮

根据加藤回想，他必须研究世代传承的祖业即兵学，但“世上喜好西洋兵学的人增多了，并且比起兵学我更喜欢研究哲学、伦理学、法学等学科了……认为对社会也有些益处，故遂改志，决定要从事自己喜欢的研究”[注]加藤弘之：《加藤弘之自叙传》，收入《传记丛书》88，东京：大空社，1991年，第26—27页。此书复刻加藤弘之先生八十岁祝贺会编：《加藤弘之自叙传：附-金婚式记事该略·追远碑建设始末》，东京：编者，1915年。 。1864(元治元)年，29岁时由幕府拔擢为直属家臣，就任开成所教授职并一职。

 

(2)

式中：为桨叶某处的展向坐标；ψ为对应的相位角；λ0为时均入流；λ1c为一阶横向入流分量；λ1s为一阶纵向入流分量。

式中：n表示变距谐波的阶次；An为相应谐波的幅值；φn为相应谐波的相位。

1.2 旋翼动力学模型

集成约束方程，建立系统的动力学方程。

动力学计算时采用拉格朗日方程

 

(3)

本文采用动力学软件ADAMS建立旋翼系统的动力学方程。ADAMS中刚体广义坐标用直角坐标和方位角表示q=[x,y,z,ψ,θ,φ]T。

 

(4)

式中：P为广义动量；H为外加力的坐标转换矩阵。

表1是模型旋翼的相关参数，将其导入ADAMS，建立桨叶与桨毂之间的铰接约束。将桨叶分成若干等分段，每段桨叶所受外加激励即气动载荷通过编译L-B模型的气动子程序定义加载。

建议：男性肝脏的承受能力是每天40克酒精，女性减半。一般40克酒精相当于含酒精6度的啤酒1000毫升，含酒精12度的红酒500毫升，含酒精度50度的白酒100毫升。

 

表1 模型旋翼参数

 

Tab.1 Parameters of modeling rotor

  

参数数值旋翼转速/(r·m-1)387桨叶片数3旋翼半径/m5.25翼型NACA0012弦长/m0.35负扭转/(°)5.65桨叶质量/kg39.64挥舞铰偏置量/m0.11变距铰偏置量/m0.25

采用NACA0012翼型进行算例验证，迎角变化为α=2.1+8.2sinωt。如图2所示，L-B气动模型计算值与试验值基本吻合。

  

图3 μ=0.14桨毂载荷值Fig.3 Hub load value at μ=0.14

2 独立桨叶控制对桨毂载荷的影响

单片桨叶控制[8-12]在输入不同阶次、幅值、相位的高阶谐波成分通过控制桨叶的桨距角来达到减振目的。桨距角可表示为

微课视频的内容要强调理论知识的重难点，对实践部分要强调过程、注重细节；在平台建设方面构建学习平台并上传学习资源；在知识储备方面，熟悉教学内容，预测可能出现的问题；在策略制定方面，教师制定教学策略，对教师的知识储备能力要求较高。

θ=θ0+θccos(ωt)+θssin(ωt)+

A2cos(2ωt+φ2)+

A3cos(3ωt+φ3)+

由于直升机流场的复杂性，本文采用准确度较高、计算较精确的Pitt-Peters[7]动力入流模型。其将诱导速度表示为由桨叶展向和方位角的一阶傅里叶级数数学表达式

苗期田间铺设银灰色膜避蚜，每667平方米铺设银灰膜5千克或将银灰膜剪成10～15厘米宽的膜条，间距15厘米，纵横拉成网眼状。

Ancos(nωt+φn)

(5)

旋翼总入流速度由悬停爬升速度和诱导速度叠加组成。

如图4所示，在算例迎角变化基础上对二维翼型加载高阶谐波分量后，翼型的迎角会发生变化，翼型的动态气动力会随之改变，在特定谐波下，可以使得翼型某一时间范围内不进入动态失速区域。

  

图4 谐波分量对气动力的影响Fig.4 Effects on the aerodynamic load by the higher harmonic pitch

为方便分析，定义桨毂载荷相对比值作为减振衡量指标。载荷比表示：某一方向上未加谐波前载荷最大值与最小值之差与加载谐波后载荷最大值和最小值之差的比值。载荷比越小，表明振动幅度越低，减振效果越好。

2.1 谐波相位对振动载荷的影响

而对于不同旋翼模型在不同的飞行状态下，虽然其有一定的差异，但应当能够找出相应的高阶谐波阶次、幅值、相位组合加载到独立桨叶变距上，使得机体的振动水平得到降低。

由图5可知，加载二阶谐波时，桨毂纵向、横向载荷比大于1，即二阶谐波对桨毂力只有激振作用，两方向力变化趋势相同；在谐波相位180°附近，桨毂载荷比小于1，对桨毂垂向力具有减振作用，当谐波相位处于180°时，减振效果最优。

  

图5 谐波相位对振动载荷比的影响Fig.5 Effects on the vibration of harmonic phase to the pitch

加载三阶谐波时，谐波相位对3个方向力载荷比影响趋势相同，对桨毂纵向力与垂向力影响较大而对横向力影响十分微小。随着谐波相位增大，振动载荷比先减小后增大，垂向力由减振逐渐变为激振；当谐波相位处于120°时，减振效果最优。

对上述计算结果分析，3片桨叶旋翼系统，其桨毂垂向振动载荷主要来源于各片桨叶气动力3n阶分量(n=1,2,3,…)。要通过高阶谐波变距减小各方向振动，就需要谐波相位与振动载荷相位匹配，否则不仅不能达到减振目的，还会变成激振。在不同状态，不同阶次谐波达到减振效果所需要的相位不同。

此诗与乾隆年间的政治斗争有关，所以对该诗主旨的准确把握，必须建立在对相关史料进行分析的基础上。《清史稿·管世铭列传》载：

2.2 谐波幅值对振动载荷的影响

在前飞速度μ=0.14时，由相位影响分析，二阶谐波和三阶谐波分别选择对应的最优减振相位值(谐波相位值分别为φ2=180°，φ3=90°)。

分别选取二阶谐波和三阶谐波独立加载，在固定相位不变的情况下，改变谐波幅值，绘制出载荷比变化曲线(图6)。幅值取值从0.1°开始，每间隔0.1°增加幅值大小，计算桨毂周期内的变化量并算出载荷比，将离散数据拟合成曲线。

由图6可知，加载二阶谐波，二阶谐波幅值较小时，对力载荷有一定的减振效果，垂向力在幅值为0.7时，减振最明显，随着谐波幅值增大，桨毂力载荷增大，会形成激振的效果。

  

图6 三阶谐波幅值对振动载荷比的影响Fig.6 Effects on the vibration of harmonic amplitude to the pitch

加载三阶谐波时，能够有效降低桨毂振动载荷，垂向和纵向载荷的减振效果最为明显，侧向力在幅值为0.2左右时，减振效果最好，之后，减振逐渐转变成激振；垂向力在幅值0.8时，减振达到最优。

总体来看，谐波幅值较大时，一方面变距拉杆激励器难以实现，另一方面高阶谐波变成了振动激励，并不能达到降低振动的目的；当谐波幅值太小时，高阶谐波作用弱化，减振效果也随之减弱。

由于低压变频器具有调速范围广、调速范围宽、运行效率高、操作方便、易于与其他设备的接口相连等一系列的优点在各行各业中得以广泛的应用，具有非常稳定的性能，使用上优于直流调速系统。另外，变频器体积较小、而且低噪音，所以在三相异步电动机上应用具有成本低、易维护的特点，不仅有效的降低了成本，同时也使工艺得以简化。对于变频器的应用，通常可以有效的提高单位生产率、提高设备自动化水平、产品质量及有效的使生活环境得以改善等。同时由于低压变频器的使用，可以有效的节约能源，并使企业的生产成本得以有效的降低。

在中医看来，花椒可以温中止痛、祛湿散寒。用花椒水泡脚和用当归、红花泡脚有异曲同工之效，都能活血通络，使整个机体血脉畅通，浑身暖融融的。此外，花椒还是一种天然的消毒剂，用花椒水泡脚还能帮助治脚气。

3 结 论

对于本文所选用的3片铰接式旋翼模型，在一定前飞状态下进行独立桨叶变距仿真，可以得出如下结论：

(1)二阶、三阶谐波都可以大幅降低桨毂垂向振动，最大降幅能达到80%，同时对桨盘拉力影响很小。

式中：Kh为磁滞损耗系数；对一般硅钢片α取1.6～2.3；f为电机磁场的交变频率；Bm为磁密幅值。建立Ansoft模型，磁滞损耗由Ansoft计算[3]，电机为分数槽时曲线如图4。

(2)不同阶次高阶谐波实现减振所需的谐波幅值不同，谐波幅值过小，减振效果弱，谐波幅值过大，减振转换成激振。

(3)不同状态下独立桨叶控制对不同阶次谐波的相位需求不同，高阶谐波相位对桨毂载荷影响较大，最大振幅与最小振幅对应相位相差180°。

在前进比μ=0.14时，分别选取二阶谐波和三阶谐波独立加载到桨距上，在固定幅值不变的情况下(谐波幅值分别为A2=1，A3=0.2)，改变谐波相位，绘制出载荷比变化曲线，如图5所示。

参考文献：

[1] 张永孝, 李强. 直升机结构响应主动控制技术工程化应用研究[J]. 飞行力学, 2015, 33(4): 371-375.

ZHANG Yongxiao, LI Qiang. Research of Engineering application for active control of helicopter structure responses[J].Flight Dynamics, 2015, 33(4): 371-375.

Design of Efficient Steering Controller for Unmanned Surface Vehicle

[2] DAN P, LIU L, CHANDRASEKAR J, et al. Higher-harmonic-control algorithm for helicopter vibration reduction revisited[J]. Journal of Guidance Control & Dynamics, 2012, 28(5):918-930.

[3] 徐海, 王华明, 杨仁国. 独立桨叶高阶谐波变距对旋翼垂向载荷的影响分析[J]. 南京航空航天大学学报, 2016, 48(2):200-204.

XU Hai,WANG Huaming,YANG Renguo. Analysis on effect of individual high harmonic blade pitch on vertical hub load[J].Journal of Nanjing University of Aeronautics & Astronautics, 2016, 48(2):200-204.

[4] LEISHMAN J G, BEDDOES T S.A semi-empirical model for dynamic stall[J]. Journal of the American Helicopter Society, 1989, 34(3):3-17.

[5] 宋辰瑶, 徐国华. 基于非定常气动模型的翼型动态响应计算[J]. 实验流体力学, 2005, 19(2):28-34.

SONG Chenyao, XU Guohua. Computations of airfoil dynamic responses based on unsteady aerodynamic model[J].Journal of Experiments in Fluid Mechanics, 2005, 19(2):28-34.

[6] LARSEN J W, NIELSEN S R K, KRENK S. Dynamic stall model for wind turbine airfoils[J]. Journal of Fluids & Structures, 2007, 23(7):959-982.

[7] GAONKAR G,PETERS D.Review of dynamic inflow modeling for rotorcraft flight dynamics [J]. Vertica, 1988, 12(3):213-242.

[8] JACKLIN S, LEYLAND J, BLAAS A, et al. Full-scale wind tunnel investigation of a helicopter individual blade control system[C]∥50th Annual Forum of the AHS. USA:AIAA,1993.

[9] ARNOLD U T P, FüRST D. Closed loop IBC results from CH-53G flight tests[J]. Aerospace Science & Technology, 2005, 9(5):421-435.

[10] YEO H, JAIN R, JAYARAMAN B. Investigation of rotor vibratory loads of a UH-60A individual blade control system[J]. Journal of the American Helicopter Society, 2011, 56(4):1-18.

2.死亡时间。据王清兰分析，体大膘肥的母猪易造成临产时胎儿死亡，死亡占出生仔猪总数的9.2%，占哺乳仔猪死亡总数的69.2%。又据河南正阳猪场的分析，3日龄以内死亡的仔猪占死亡总数的26.63%，4～7日龄占29.27%，8～15日龄占20.21%，16～20日龄占9.92%，21～25日龄占7.27%，26～35日龄占2.17%，36～45日龄占1.15%，46～60日龄占3.02%。加拿大对6 890头仔猪的分析，从出生到20日龄，仔猪死亡率为25.6%，其中分娩时死亡占死亡总数的15.3%，7日龄以内死亡占43.7%。可见，仔猪日龄越小，死亡率越高。

[11] JACKLIN S A, HABER A, SIMONE G D, et al. Full-scale wind tunnel test of an individual blade control system for a UH-60 helicopter [C]//AHS 58th Annual Forum. Montreal, Canada: American Helicopter Society Inernational, 2002.

[12] RAVICHANDRAN K, CHOPRA I, WAKE B E, et al. Trailing-edge flaps for rotor performance enhancement and vibration reduction[J]. Journal of the American Helicopter Society, 2013, 58(2):1-13.