首页 / 专利库 / 地热能 / 地热循环 / 闭环系统 / 一种基于地面频响试验的气动伺服弹性稳定性预测方法

一种基于地面频响试验的气动伺服弹性稳定性预测方法

阅读:186发布:2020-05-14

专利汇可以提供一种基于地面频响试验的气动伺服弹性稳定性预测方法专利检索,专利查询,专利分析的服务。并且本 发明 提出了一种基于地面频响试验的 气动 伺服弹性 稳定性 预测方法。该方法以真实 飞行器 为试验对象,通过地面频响试验测得飞行器结构‑控制环节相关传递函数,通过气动 力 理论计算得到气动力影响系数矩阵,结合两者结果预测气动伺服弹性稳定性。该方法将气动伺服弹性系统中结构‑控制环节和非定常气动力环节分开考虑,避免了模拟加载的诸多 缺陷 ;在控制回路开环状态下进行地面频响试验,避免了控制回路闭环失稳的 风 险,试验安全性高;通过地面频响试验和非定常气动力理论计算相结合的方式进行气动伺服弹性稳定性预测,预测 精度 高,具有一定的工程应用价值。,下面是一种基于地面频响试验的气动伺服弹性稳定性预测方法专利的具体信息内容。

1.一种适用于真实飞行器的基于地面频响试验的气动伺服弹性稳定性预测方法,其特征在于包括:
A)将飞行器的气动伺服弹性系统分为结构自由度和控制自由度,结构自由度由各气动段下洗控制点的转zθ和纵向速度zv组成的广义坐标z=[zθ zv]T描述,控制自由度由偏角信号δ描述,
B)把飞行器的机体划分为M个气动段,其中舵面被视为刚体并作为单独的气动段,C)对于N阶振动模态的飞行器,用以下参数表征所述气动伺服弹性系统:
舵面耦合惯性质量Mδ,其维数为N×1,
舵面气动影响系数Aδ,其维数为1×1,
全机气动力影响系数矩阵Az,其维数为(M+1)×2(M+1),
舵面惯性力fMδ到广义坐标z的传递函数P11,其维数为2(M+1)×N,
由舵偏产生的非定常气动力fAδ到广义坐标z的传递函数P12,其维数为2(M+1)×1,由结构振动产生的非定常气动力fAz到广义坐标z的传递函数P13,其维数为2(M+1)×(M+
1),
舵面惯性力fMδ到惯性测量元件输入角速度 的传递函数P21,其维数为1×N,由舵偏产生的非定常气动力fAδ到惯性测量元件输入角速度 的传递函数P22,其维数为
1×1,
由结构振动产生的非定常气动力fAz到惯性测量元件输入角速度 的传递函数P23,其维数为1×(M+1),
舵机传递函数KS,其维数为1×1,
惯性测量元件传递函数KI,其维数为1×1,
飞行控制系统传递函数KC,其维数为1×1,
D)进行地面频响试验,包括在气动伺服弹性系统开环情况下执行以下操作:
D1)以输入舵机的舵偏指令u作为输入信号,采集舵机输出的舵偏角信号δ,测得舵机传递函数KS;采集飞行控制系统解算指令uc,测得舵偏指令到解算指令的传递函数Tcs=-s2KCKIP21MδKS,其维数为1×1;采集各气动段下洗控制点的转角zθ和纵向速度zv,测得舵偏指令到广义坐标z的传递函数Tzs=-s2P11MδKS,其维数为2(M+1)×1,其中s表示拉氏变量,D2)用激振器在舵面气动段压心处进行激励,采集飞行控制系统解算指令uc,测得由舵偏产生的非定常气动力到解算指令的传递函数Tcδ=KCKIP22,其维数为1×1;采集各气动段下洗控制点的转角zθ和纵向速度zv,测得由舵偏产生的非定常气动力到广义坐标z的传递函数P12,其维数为2(M+1)×1,
D3)用激振器在每个气动段压心处进行激励,采集飞行控制系统解算指令uc,测得由结构振动产生的非定常气动力到解算指令的传递函数Tcz=KCKIP23,其维数为1×(M+1);采集各气动段下洗控制点的转角zθ和纵向速度zv,测得由结构振动产生的非定常气动力到广义坐标z的传递函数P13,其维数为2(M+1)×(M+1),
E)进行基于非定常气动力理论的算法处理,其中所述非定常气动力理论是从偶极子格网法和气动导数法中选出的一种,
定义全机气动力影响系数矩阵Az为:
其中,S表示各气动段气动力计算参考面积组成的对角矩阵,其维数为(M+1)×(M+1),D表示相关工程非定常气动力理论求得的气动力影响系数矩阵,其维数为(M+1)×(M+1);w表示各气动段下洗控制点位移,其维数为(M+1)×1,ρ表示大气密度,V表示飞行速度,按照下式确定舵面气动力影响系数Aδ:
其中,Sδ表示舵面气动力计算参考面积,Dδ表示相关工程非定常气动力理论求得的舵面气动力影响系数,δ(t)表示舵偏角信号且δ(t)=δ,
F)结合地面频响试验测得的传递函数和气动力影响系数矩阵Az和Aδ,令拉氏变量s=iω,得到气动伺服弹性系统的开环频响函数,
G(iω)=Tcs+TczAz(I-P13Az)-1Tcs+TcδAδKS+TczAz(I-P13Az)-1P12AδKS (3)其中I表示单位矩阵,
根据奈奎斯特稳定性判据,通过开环频响函数的幅相曲线,判断闭环系统的稳定性,包括:
对于气动伺服弹性系统,飞行器的弹性系统是稳定的,控制回路开环是稳定的,当奈奎斯特曲线包围临界点(-1,0)时,则判定气动伺服弹性系统是稳定的,否则判定气动伺服弹性系统是不稳定的,
判断不同速度下气动伺服弹性系统开环频响函数G(iω)的奈奎斯特曲线是否包围临界点(-1,0),从而得到飞行器气动伺服弹性临界稳定速度,通过气动伺服弹性系统开环频响函数G(iω)的伯德图,得到飞行器在不同速度下的幅值裕度和相位裕度。

说明书全文

一种基于地面频响试验的气动伺服弹性稳定性预测方法

技术领域

[0001] 本发明涉及一种基于地面频响试验对真实飞行器进行气动伺服弹性稳定性预测的方法。

背景技术

[0002] 气动伺服弹性问题是飞行器弹性结构、非定常气动以及飞行控制系统相互耦合的动气动弹性问题。飞行器气动伺服弹性不稳定会造成结构疲劳损伤、降低操纵系统性能、甚至导致严重的结构破坏,因此在飞行器设计中必须重视气动伺服弹性稳定性。
[0003] 目前研究气动伺服弹性问题的主要途径有三种:一是数值计算,这种方法分析对象是飞行器的气动弹性数学模型,引入了较多理论假设,不能完全反映真实飞行器对象的动态特性;二是洞试验,该方法采用真实飞行器的缩比模型,且受到风洞性能条件的诸多限制,与真实飞行器的真实情况还是存在一定差别;三是飞行试验,通过进行真实飞行器飞行包线内的飞行试验,验证其气动伺服弹性稳定性,结果真实可靠,但试验成本高、风险大,因此需要先进行飞行器气动弹性相关的数值计算、地面试验和风洞试验。

发明内容

[0004] 近年来提出了一种气动伺服弹性地面模拟试验方法,通过飞行器非定常气动力的模拟加载,从而以真实飞行器为试验对象完成气动伺服弹性地面模拟试验。该方法与传统方法相比具有低成本、低风险的优势,可以成为传统方法的有力补充。该方法的技术关键在于非定常气动力的模拟加载,相关研究机构采用了诸多方法进行非定常气动力的计算及模拟加载,其结果具有一定的工程应用价值,但其受到加载系统的局限,非定常气动力的模拟加载都难以准确实现,也影响了气动伺服弹性稳定性预测的精度
[0005] 针对现有技术的上述问题,本发明人提出了一种基于地面频响试验的气动伺服弹性稳定性预测方法,其将非定常气动力环节与结构-控制环节分开考虑,结合真实飞行器的地面频响试验和非定常气动力理论计算进行气动伺服弹性稳定性预测。
[0006] 基于根据本发明的一种基于地面频响试验的气动伺服弹性稳定性预测方法,可以将气动伺服弹性系统结构-控制环节和非定常气动力环节分开考虑,以真实飞行器为试验对象,开展开环系统地面频响试验,结合试验测得的相关结构-控制环节传递函数和非定常气动力理论计算得到的气动力影响系数矩阵,进行气动伺服弹性稳定性预测。
[0007] 本发明提供了一种适用于真实飞行器的基于地面频响试验的气动伺服弹性稳定性预测方法,其特征在于包括:
[0008] A)将飞行器的气动伺服弹性系统分为结构自由度和控制自由度,结构自由度由各气动段下洗控制点的转zθ和纵向速度zv组成的广义坐标z=[zθ zv]T描述,控制自由度由偏角信号δ描述,
[0009] B)把飞行器的机体划分为M个气动段,其中舵面被视为刚体并作为单独的气动段,[0010] C)对于N阶振动模态的飞行器,用以下参数表征所述气动伺服弹性系统:
[0011] 舵面耦合惯性质量Mδ,其维数为N×1,
[0012] 舵面气动力影响系数Aδ,其维数为1×1,
[0013] 全机气动力影响系数矩阵Az,其维数为(M+1)×2(M+1),
[0014] 舵面惯性力fMδ到广义坐标z的传递函数P11,其维数为2(M+1)×N,[0015] 由舵偏产生的非定常气动力fAδ到广义坐标z的传递函数P12,其维数为2(M+1)×1,[0016] 由结构振动产生的非定常气动力fAz到广义坐标z的传递函数P13,其维数为2(M+1)×(M+1),
[0017] 舵面惯性力fMδ到惯性测量元件输入角速度 的传递函数P21,其维数为1×N,[0018] 由舵偏产生的非定常气动力fAδ到惯性测量元件输入角速度 的传递函数P22,其维数为1×1,
[0019] 由结构振动产生的非定常气动力fAz到惯性测量元件输入角速度 的传递函数P23,其维数为1×(M+1),
[0020] 舵机传递函数KS,其维数为1×1,
[0021] 惯性测量元件传递函数KI,其维数为1×1,
[0022] 飞行控制系统传递函数KC,其维数为1×1,
[0023] D)进行地面频响试验,包括在气动伺服弹性系统开环情况下执行以下操作:
[0024] D1)以输入舵机的舵偏指令u作为输入信号,采集舵机输出的舵偏角信号δ,测得舵机传递函数KS;采集飞行控制系统解算指令uc,测得舵偏指令到解算指令的传递函数Tcs=-s2KCKIP21MδKS,其维数为1×1;采集各气动段下洗控制点的转角zθ和纵向速度zv,测得舵偏指令到广义坐标z的传递函数Tzs=-s2P11MδKS,其维数为2(M+1)×1,其中s表示拉氏变量,[0025] D2)用激振器在舵面压心处进行激励,采集飞行控制系统解算指令uc,测得由舵偏产生的非定常气动力到解算指令的传递函数Tcδ=KCKIP22,其维数为1×1;采集各气动段下洗控制点的转角zθ和纵向速度zv,测得由舵偏产生的非定常气动力到广义坐标z的传递函数P12,其维数为2(M+1)×1,
[0026] D3)用激振器在每个气动段压心处进行激励,采集飞行控制系统解算指令uc,测得由结构振动产生的非定常气动力到解算指令的传递函数Tcz=KCKIP23,其维数为1×(M+1);采集各气动段下洗控制点的转角zθ和纵向速度zv,测得由结构振动产生的非定常气动力到广义坐标z的传递函数P13,其维数为2(M+1)×(M+1)。
[0027] E)进行基于非定常气动力理论的算法处理,其中所述非定常气动力理论是从偶极子格网法和气动导数法中选出的一种,
[0028] 定义全机气动力影响系数矩阵Az为:
[0029]
[0030] 其中,S表示各气动段气动力计算参考面积组成的对角矩阵,其维数为(M+1)×(M+1),D表示相关工程非定常气动力理论求得的气动力影响系数矩阵,其维数为(M+1)×(M+
1);w表示各气动段下洗控制点位移,其维数为(M+1)×1,ρ表示大气密度,V表示飞行速度,[0031] 按照下式确定舵面气动力影响系数Aδ:
[0032]
[0033] 其中,Sδ表示舵面气动力计算参考面积,Dδ表示相关工程非定常气动力理论求得的舵面气动力影响系数,δ(t)表示舵偏角信号且δ(t)=δ,
[0034] F)结合地面频响试验测得的传递函数和气动力影响系数矩阵Az和Aδ,令拉氏变量s=iω,得到气动伺服弹性系统的开环频响函数,
[0035] G(iω)=Tcs+TczAz(I-P13Az)-1Tcs+TcδAδKS+TczAz(I-P13Az)-1P12AδKS   (3)[0036] 其中I表示单位矩阵,
[0037] 根据奈奎斯特稳定性判据,通过开环频响函数的幅相曲线,判断闭环系统的稳定性,包括:
[0038] 对于气动伺服弹性系统,飞行器的弹性系统是稳定的,控制回路开环是稳定的,当奈奎斯特曲线包围临界点(-1,0)时,则判定气动伺服弹性系统是稳定的,否则判定气动伺服弹性系统是不稳定的,
[0039] 判断不同速度下气动伺服弹性系统开环频响函数G(iω)的奈奎斯特曲线是否包围临界点(-1,0),从而得到飞行器气动伺服弹性临界稳定速度,
[0040] 通过气动伺服弹性系统开环频响函数G(iω)的伯德图,得到飞行器在不同速度下的幅值裕度和相位裕度。附图说明
[0041] 图1为根据本发明的一个实施例的弹箭类飞行器示意图;
[0042] 图2为根据本发明的一个实施例的气动伺服弹性系统的框图
[0043] 图3为根据本发明的一个实施例的气动伺服弹性稳定性预测方法流程图
[0044] 图4A-图4C为运用根据本发明的实施例进行气动伺服弹性稳定性预测结果与标称系统的对比,其中,图4A为幅频特性曲线;图4B为相频特性曲线;图4C为乃奎斯特曲线。
[0045] 附图标记:
[0046] 1—弹体 2—舵面 3—舵机 4—激振器
[0047] 5—飞行控制系统 6—惯性测量元件 7—加速传感器
[0048] 具体实施方法
[0049] 图3是根据本发明的一个实施例的一种基于地面频响试验的气动伺服弹性稳定性预测方法的流程图,其以图1所示的弹箭类飞行器为应用对象,包括以下步骤:
[0050] (1)地面频响试验
[0051] 将弹体1沿轴向均匀划分为M个气动段,舵面2被视为刚体并作为单独气动段,舵面2下洗控制点和压心被等效投影到弹体1上。
[0052] 在各气动段下洗控制点两侧近距离布置两个加速度传感器7。通过一次积分,根据几何关系,测得下洗控制点的纵向速度zv;两次积分,测得下洗控制点的转角zθ,记为广义坐标z=[zθzv]T。
[0053] 图2是该弹箭类飞行器的气动伺服弹性系统框图,在气动伺服弹性系统开环情况下,从201处向舵机3输入舵偏指令u,采集202处舵机3的输出舵偏角信号δ,测得舵机3的传递函数KS;采集206处飞行控制系统5的解算指令uc,测得舵偏指令到解算指令的传递函数Tcs=-s2KCKIP21MδKS;采集205处各气动段下洗控制点的转角zθ和纵向速度zv,测得舵偏指令到下洗控制点广义坐标的传递函数Tzs=-s2P11MδKS。
[0054] 在203表示的弹体1上舵面2压心投影点处用激振器进行激励,采集206处飞行控制系统5的解算指令uc,测得由舵偏产生的非定常气动力到解算指令的传递函数Tcδ=KCKIP22;采集205处各气动段下洗控制点的转角zθ和纵向速度zv,测得由舵偏产生的非定常气动力到广义坐标z的传递函数P12。
[0055] 在204表示的每个气动段压心处用激振器分别进行激励,采集206处飞行控制系统解算指令uc,测得由结构振动产生的非定常气动力到解算指令的传递函数Tcz=KCKIP23;采集205处各气动段下洗控制点的转角zθ和纵向速度zv,测得由结构振动产生的非定常气动力到广义坐标z的传递函数P13。
[0056] (2)基于非定常气动力理论的算法处理
[0057] 该实施例弹箭类飞行器采用气动导数法进行基于非定常气动力理论的算法处理,[0058] 全机气动力影响系数矩阵Az由以下式子确定:
[0059]
[0060] 其中,S表示各气动段气动力计算参考面积组成的对角矩阵,CLα表示各气动段气动导数组成的对角矩阵;
[0061] 舵面气动力影响系数Aδ由以下式子确定:
[0062]
[0063] 其中,Sδ表示舵面气动力计算参考面积, 表示舵面气动导数。
[0064] (3)气动伺服弹性稳定性预测
[0065] 结合地面频响试验传递函数和非定常气动力理论计算得到的气动力影响系数矩阵,令拉氏变量s=iω,得到气动伺服弹性系统开环频响函数,
[0066] G(iω)=Tcs+TczAz(I-P13Az)-1Tcs+TcδAδKS+TczAz(I-P13Az)-1P12AδKS   (3)[0067] 根据奈奎斯特稳定性判据,通过开环频响函数的幅相曲线,判断闭环系统的稳定性。对于气动伺服弹性系统,飞行器的弹性系统是稳定的,控制回路开环是稳定的,因此只要奈奎斯特曲线包围临界点(-1,0),则判定气动伺服弹性系统是稳定的,否则判定气动伺服弹性系统是不稳定的。
[0068] 通过判断不同速度下气动伺服弹性系统开环频响函数G(iω)的奈奎斯特曲线是否包围临界点(-1,0),得到飞行器气动伺服弹性临界稳定速度。
[0069] 通过气动伺服弹性系统开环频响函数G(iω)的伯德图,得到飞行器在不同速度下的幅值裕度和相位裕度。
[0070] 对于图1所示弹箭类飞行器,在速度V=450m/s时,应用根据本发明的上述实施例的方法,得到了如图4A所示该实施例幅频特性曲线、如图4B所示相频特性曲线、以及图4C所示乃奎斯特曲线,可以看出根据本发明的该方法得到的气动伺服弹性系统开环频响函数与标称系统计算结果吻合,具有高的精度。
[0071] 本发明的优点包括:
[0072] 将气动伺服弹性系统中结构-控制环节和非定常气动力环节分开考虑,采用非定常气动力理论方法进行非定常气动力计算,避免了模拟加载的诸多缺陷;在控制回路开环状态下进行地面频响试验,避免了控制回路闭环失稳的风险,试验安全性高;通过地面频响试验和非定常气动力理论计算相结合的方式进行气动伺服弹性稳定性预测,预测精度高,具有一定的工程应用价值。
高效检索全球专利

专利汇是专利免费检索,专利查询,专利分析-国家发明专利查询检索分析平台,是提供专利分析,专利查询,专利检索等数据服务功能的知识产权数据服务商。

我们的产品包含105个国家的1.26亿组数据,免费查、免费专利分析。

申请试用

分析报告

专利汇分析报告产品可以对行业情报数据进行梳理分析,涉及维度包括行业专利基本状况分析、地域分析、技术分析、发明人分析、申请人分析、专利权人分析、失效分析、核心专利分析、法律分析、研发重点分析、企业专利处境分析、技术处境分析、专利寿命分析、企业定位分析、引证分析等超过60个分析角度,系统通过AI智能系统对图表进行解读,只需1分钟,一键生成行业专利分析报告。

申请试用

QQ群二维码
意见反馈