基于多项式混沌法的翼型不确定性分析及梯度优化设计

耦合伴随方法和非嵌入式多项式混沌法，发展了高效、可靠的不确定性梯度优化设计方法。利用伴随方程法求解目标函数对不确定性变量的导数，发展了一种梯度增强型多项式混沌法。通过亚声速和跨声速下等多种算例可以证明该方法可以提高不确定性分析的效率和精度。同时，利用基于方差分解的全局敏感性分析方法对不确定性变量的敏感性进行了量化。建立了多项式混沌耦合伴随方程的统计矩梯度求解方法，并结合梯度增强型多项式混沌法搭建不确定性梯度优化设计系统。基于该优化设计系统对二维低亚声速和跨声速翼型开展确定性及不确定性优化设计研究。优化结果显示，相比于确定性优化设计，不确定性优化设计通过合理权衡确定性性能和不确定性性能，可提高抵抗马赫数和迎角不确定性扰动的能力，同时优化性能均值和标准差。其中阻力系数均值最大可降低17%，阻力系数标准差最大可降低80%。而确定性优化设计可能导致性能鲁棒性的降低。     

DE历表的结构、计算与比较

在详细介绍美国JPL(Jet Propulsion Laboratory,喷气推进实验室)发布的DE (Development Ephemerides,行星/月球展开历表)系列历表的基本情况、头文件和数据文件的结构、历表计算原理等内容的基础上,以DE405、DE421、DE431为例,分别计算了1900-2200年间月球、太阳、行星等天体在地心天球参考系中的坐标及速度,并由此得到3个历表间的差异,结果表明DE431与DE421之间的差异,整体而言要小于二者同DE405之间的差异.最后以引潮力的计算为例,表明历袁更新对于引潮力的影响基本可忽略不计,验证了IERS(International Earth Rotation and reference systems Service,国际地球自转与参考系服务)规范中关于历表更新对相关地球科学影响十分有限的结论.     

航空发动机增益调度控制的多项式平方和规划方法

针对现有的线性变参数(linear parameter varying,LPV)控制器设计方法都是关于仿射参数依赖系统而没有专门针对多项式描述的LPV系统这一现状,提出了一种基于多项式平方和(sum of squares,SOS)规划的增益调度控制设计方法,并将其用于转速大范围变化时的航空发动机高压转子转速及压比控制.根据发动机非线性模型获取不同转速下的状态空间模型,并利用多项式拟合的方法建立发动机线性变参数模型.给出能够保证无静差的增益调度控制结构,利用有界实定理和多项式平方和理论推导出能够保证闭环系统鲁棒稳定的SOS约束条件,并形成控制器求解的SOS规划问题,通过求解获得多项式描述的增益调度控制器.分别以LPV模型和发动机非线性模型为对象做阶跃仿真,结果表明:高压转子转速/发动机压比控制系统的调节时间在2s以内,稳态误差不超过0.1%.      

涡轮叶片温度场插值的坐标变换法

研究了利用坐标变换求涡轮叶片温度场的插值方法。选择等叶高面间的叶片体为研究对象,利用20个节点包络此叶片体,在物理坐标与虚拟坐标之间建立坐标变换函数,确立物理坐标节点与虚拟坐标节点的映射关系,通过坐标变换把形状复杂的叶片体变换为规则的正方体,在虚拟坐标系中、规则的正方体上插值求得涡轮叶片体任意点的温度值。通过算例验证了该方法的可行性和正确性。     

机载雷达罩热载荷仿真分析关键技术研究

研究了机载雷达罩热载荷仿真分析的3个关键技术,一是雷达罩及天线的CAD建模技术,并应用Fluent计算罩体外表面温度;二是建立雷达罩的多腔体热辐射模型,以得到雷达罩内表面的温度分布;三是温度插值技术,用以将流场计算得到的温度值插值到对应的热应力分析模型上。     

基于直接配点法的滑翔轨迹快速优化设计

介绍了基于五次Gauss-Lobatto多项式的直接配点法在再入飞行器三维轨迹最优化问题中的应用。首先给出了再入飞行器轨迹优化问题模型,其中运动方程为三自由度模型,性能指标选为到达指定地点飞行时间最短,控制变量则为无量纲升力系数和倾侧角。再入飞行过程中受到加热率、过载和动压约束,终端状态受到目标位置约束。然后,应用直接配点法将最优控制问题离散化为非线性规划问题,将动态优化问题转化为静态参数最优化问题。选取各节点和配点上的状态量和控制量作为优化参数。最后应用SNOPT软件包对参数最优化问题进行求解。仿真结果表明直接配点法对于再入飞行器轨迹初始参数取值不敏感,且求解过程具有一定的实时性。     

稳态温度场作用下涡轮叶片振动特性的研究

根据一维线性分段插值方法对一维温度场进行分段插值拟合,为有限元方法分析温度场作用下叶片振动问题中温度载荷的施加提供了简化可行的方案,利用有限元方法验证该线性拟合方案的可行性,并采用多种方法对计算结果进行验证分析,对一维稳态温度场的分段插值进行了误差分析,得到稳态温度场作用下叶片振动问题的有限元分析计算方法和温度载荷施加方案。     

带旋转修正的弹簧-TFI混合动网格方法

网格存在大变形时,传统的超限插值(TFI)动网格方法易于引起网格正交性问题.依据几何关系并考虑其插值特点,提出了一种带旋转修正的TFI动网格方法.对结构网格块进行分块,在子块角点间建立弹簧元,并基于弹簧类比法计算角点位移;子网格内部结点位移采用修正后的TFI动网格方法进行计算,从而形成具有弹簧-TFI混合特征的动网格新...     

多场耦合计算平台与高超声速热防护结构传热问题研究

从有限元法(FEM)和有限体积法(FVM)的单元特性出发,提出一种具有局部守恒特性的界面载荷插值方法.采用共享内存技术开发适用于通用有限元和计算流体力学(CFD)软件的多场耦合计算平台,并基于分区耦合方式实现流固耦合传热计算.作为验证,分别将ANSYS与Fluent和CFD-FASTRAN软件进行耦合,计算外壁冷却的喷...     

静气弹中非线性气动力建模方法与分析

大展弦比低雷诺数气动布局容易较早出现气流分离,会带来明显的非线性气动力问题。针对此类布局提出了一种建立基于Kriging插值的非线性压力系数分布模型的方法。从Navier-Stokes方程计算的不同状态下飞机的压力系数中提取不同坐标的系数。利用Kriging插值函数建立CFD压力系数对迎角导数的响应面,将插值结果代入偶极子网格法(Double-Lattice Method,DLM)修正其线性方法。利用无限板样条(IPS)方法进行气动结构耦合,实现了有限元结构的非线性气弹响应分析。算例结果验证了方法对于静气弹分析的有效性,同时能准确地反映弹性带来的气动效率的降低和非线性力矩特征。