基于惯性释放的细节有限元模型分析

本文简要地介绍了惯性释放的概念与原理,通过两个算例说明在静力学求解问题中将惯性释放应用于细节模型中的具体方法,将算例的计算结果与试验结果对比说明基于惯性释放的细节有限元模型分析方法的正确性,并且阐述了该方法的后续发展方向。     

高速发动机油膜惯性对活塞裙润滑影响

利用一种新的迭代算法 ,分析了高速发动机油膜惯性对活塞裙润滑特性的影响 .该迭代方法是利用有限元法和差分法交替求解Navier Stocks方程和雷诺方程 ,据此导出了适用于高速发动机活塞裙的润滑计算的混和润滑雷诺模型 .新的模型借助惯性系数 ,引入了油膜惯性项 ;同时给出求解含有油膜惯性项的迭代步骤和有限元表达式。计算结果表明 :随着惯性系数和活塞裙的长径比的不同 ,油膜惯性会对油膜的摩擦力、压力和承载力产生不同的影响 ,这种影响对承载力尤为明显 .该润滑模型也可用于中、低速发动机的活塞裙润滑计算以及不计入惯性项 (惯性系数置为零 )的某些润滑问题求解 .     

基于Broyden法的旋翼多体系统气动弹性分析

建立了旋翼多体系统气动弹性模型并给出了一种适合于该模型响应计算的数值计算方法。采用柔性多体系统动力学方法建立旋翼气动弹性模型,利用驱动约束显著简化约束方程形式,集成大变形桨叶模型,准确考虑变形的非线性,适合于对采用柔性结构的先进旋翼进行气动弹性分析。基于Broyden法改进隐式积分法积分一步中非线性方程的求解,避免求取切线矩阵和矩阵求逆运算,保持隐式积分法具有较好稳定性的同时提高计算效率,解决了旋翼多体系统气动弹性力学方程隐式表达且具有较强非线性和较高刚性比造成的响应计算困难。通过模型旋翼桨叶响应计算验证了结构模型与气动弹性响应求解方法。采用建立的气动弹性模型计算悬停和前飞状态旋翼气动弹性稳定性,与试验结果对比验证了模型的正确性。研究了不同的稳定性计算方法、桨叶结构模型和入流模型等对悬停和前飞稳定性计算的影响,结果表明本文所采用的结构、气动模型及气动弹性稳定性计算方法提高了气动弹性稳定性分析精度。     

状态空间气动弹性模型及其在颤振计算中的应用

本文评述了已有的状态空间气动弹性模型,并对它们提出了一些改进意见。本文还提出了一种新的状态空间气动弹性模型。用两种机翼为例进行了颤振计算。算例表明:这种新的状态空间气动弹性模型是一种精度较高而阶数较低的模型。本文还提出了用状态空间气动弹性模型进行颤振计算时自动识别模态及自动确定颤振点的方法。     

飞行器气动弹性风洞试验技术综述

研究气动弹性问题的主要手段有数值计算、风洞试验和飞行试验三种,气动弹性风洞试验具有可靠性高(相比数值计算)和代价低(相比飞行试验)等优势,已成为航空航天飞行器气动弹性性能评估和校核的重要手段。以亚声速和跨声速气动弹性风洞试验技术为主,分别从静气动弹性、颤振和阵风试验三个角度,阐述了国内外在模型设计、数据采集处理、模型支撑、阵风发生装置、阵风载荷减缓等方面开展的主要研究工作,总结了气动弹性风洞试验在飞行器研制中的重要意义,并对我国未来气动弹性试验能力的发展提出几点建议。     

挠性捷联惯性系统惯性器件两种模型的研究

本文提出了挠性捷联惯性系统惯性器件的两种常用的稳态模型，并指出它们之间的换算关系。这两种模型在标定，对准，导航和参数补偿过程都有一定的区别，文章根据实际经验指出了模型应用的范围。     

应用惯性释放方法的静气动弹性发散分析

翼面发散作为静气动弹性的重点分析对象,传统的分析方法通常将翼面根部固支约束,不考虑刚体模态的影响。根据咨询通告AC 25.629-1B 的要求,刚体模态或短周期模态也可能会造成发散。因此,在民用飞机静气动弹性分析及适航符合性验证中,要计及刚体模态对发散的影响。通过对惯性释放方法的研究,在发散分析中计及刚体模态,给出考虑刚体模态的模态法发散分析方法;以某型民用飞机为例,应用本文方法分别计算空机状态、不同燃油和商载的静气动弹性发散结果,并和不考虑刚体模态的模态法、考虑/不考虑刚体模态的颤振法的发散计算结果进行对比分析。结果表明：考虑刚体模态后的发散速度有所提高,临界发散模态不变,发散速度随燃油和商载的增加而降低。     

飞行员跳伞模拟训练系统人伞惯性参数研究

针对飞行员跳伞模拟训练系统研发需要,对人伞惯性参数进行了研究.首先,建立了人伞系统体固定坐标系,并在此基础上给出了伞衣、伞绳的惯性参数计算方法;其次,基于Hanavan人体模型,给出了人体模型各环节的重量分配、质心计算、惯性参数的计算方法;再次,基于转动惯量平行轴定理,给出了人伞系统体固定坐标系下惯性参数的计算方法;最...     

基于非定常气动力辨识技术的气动弹性数值模拟

选择离散型输入输出差分模型,运用最小二乘方法进行非定常气动力建模,并将辨识得到的降阶模型用于气动弹性的数值模拟。1个马赫数下的颤振临界点的计算仅需调用一次非定常流场求解器。计算精度保持与非定常欧拉方程计算方法相当的同时计算效率提高了1~2个量级。计算了跨声速具有S型颤振边界的气动弹性标准算例－Isogaiwing和三维气动弹性标模算例AGARD445.6,辨识模型计算边界与非定常Euler方程计算结果吻合。证明非定常气动力辨识技术可以提供高效的高精度的气动弹性分析。     

车载惯性定位标定的理论推导及误差分析

为了解决车载惯性定位系统误差常用的起点纬度标定计算法在安装偏差角90°附近时误差较大已不适用的问题,以及充实车载惯性定位系统标定计算的理论,首先进行了车载惯性定位系统误差标定计算的推导,证明了误差的可标定性。借鉴等角航线反解理论,得到车载惯性定位系统标定理论依据。参考中分纬度理论,推导出使用平均纬度进行简化计算的航程及航向角误差公式,将其引入惯性定位系统标定计算,提出了一种使用平均纬度代替起点纬度进行标定计算的方法,并对两种标定计算方法的误差进行了公式推导及仿真分析。仿真结果表明,使用平均纬度进行标定计算的误差明显小于使用起点纬度,在安装偏差角较大时,使用平均纬度进行标定计算可使标定计算误差降低90%。证明了提出的标定计算方法在南北纬55°、行程不超过90km时可用,高纬度、更长行程下可参照文中公式进行误差估计。     

一种基于光固化快速成型的飞机静弹性风洞试验模型设计方法

高速风洞静弹性模型设计和制造是静弹性风洞试验的一个关键。为解决模型设计周期长、制造费用高等问题,提出了一种基于立体光固化快速成型面向高速风洞大展弦比机翼静弹性模型研制方法。基于机翼刚度分布相似参数,采用机翼钢梁骨架和树脂蒙皮组合结构,通过优化结构尺寸完成静弹性模型结构设计;使用机械加工和快速成型技术完成模型制造,并通过地面刚度试验对加工模型进行了刚度分布验证。风洞试验结果表明：基于立体光固化成型技术设计和制造的静弹性风洞试验模型工程实用、可行,与传统静弹性模型研制过程相比,具有研制周期短、成本低而且不存在因填充物带来附加刚度的显著优势。     

基于非线性试验气动力的飞机静气动弹性响应分析

分别使用线性和非线性静气动弹性分析方法,对某飞机纵向静气动弹性响应特性随飞行动压、攻角、平尾偏度和纵向过载变化的趋势进行了分析,并将分析结果和飞行试验进行了比较。线性方法的气动力计算使用亚音速偶极子格网法。非线性方法由风洞试验提供刚体气动力,并使用线性气动力影响系数矩阵对其进行弹性化处理。通过对两种方法的计算结果进行对比分析可以看出:①使用刚体试验气动力的非线性分析方法能够获得和飞行试验比较一致的结果;②在线性方法所提供的结果中,翼面的弯曲变形及剪力、弯矩和扭矩等部分参数能为初步设计提供大致参考,但气动力系数的弹性增量、翼面的扭转变形、翼面的剪力、弯矩和扭矩的弹性部分可能不能为型号设计提供正确的指导;③不论刚体气动力为线性还是非线性,气动力系数随动压变化的趋势均呈线性。     

大展弦比机翼的静气动弹性分析

大型飞机普遍采用的大展弦比机翼,其气动弹性问题显得尤为突出。通过求解基于三维结构网格的Euler/边界层耦合方程组,采用结构模态分析法,对HIRENASD机翼的静气动弹性进行了数值模拟。计算得到了机翼静弹变形前和变形后的压力分布以及不同攻角下翼尖变形值。计算结果和实验值对比分析,验证了Euler/边界层耦合方程组求解器对静气动弹性数值模拟的准确性,表明所发展的数值模拟方法可以用于大展弦比机翼静气动特性的分析研究。     

复材机翼气弹特性工程化设计分析技术

本文在静气弹的基本原理分析基础上,发展了机翼气弹特性工程化设计方法：基于参数化的机翼气动弹性模型,研究了复合材料机翼整体蒙皮重要设计参数对舵效、静变形和静强度的影响,采用变参分析方法进行了灵敏度分析,获取影响舵效的关键参数变化趋势,用以指导舵效气动弹性剪裁优化设计。对优化设计结果进行了工程化的处理并进行校核,获得了工程可用的蒙皮铺层设计方案,在改善舵效使其满足飞机性能要求的同时,还获得了满意的减重105 kg（总质量的14%）效益。     

高精度的机翼静气动弹性数值模拟方法研究

气动结构耦合是研究静气动弹性问题的关键所在。采用三维Euler方程为控制方程,以有限元方法计算结构弹性变形,用细平板样条法(TPS)实现气动结构数据间的转换,进行机翼静气动弹性数值的模拟。算例结果表明,本方法是一种具有很高计算精度和工程实用性的气动弹性模拟方法。     

活塞理论气动力静压修正方法及其在曲壁板颤振分析中的应用

提出了一种采用计算流体力学(CFD)计算的压力分布对活塞理论气动力进行静压修正的方法,将该方法应用到曲壁板的静气动弹性变形及颤振稳定性分析中,并与采用曲率修正活塞理论气动力的计算结果进行了对比.分析结果表明,采用本文提出的活塞理论气动力静压修正方法进行曲壁板的气动弹性分析,在圆柱曲壁板曲率较小的情况下,与采用曲率修正活塞理论气动力方法得到的静气动弹性变形、稳定性边界差别不大;而在曲率较大时,采用本文方法计算得到的曲壁板静气动弹性变形,其曲壁板靠近前缘部分被压的更低,而曲壁板的颤振稳定性边界更小,且这种差别随着圆柱曲壁板曲率的增加而不断增大.该方法突破了曲率修正活塞理论的小曲率限制,扩大了活塞理论气动力在曲壁板颤振分析中的适用范围.     

充气机翼的褶皱与失效行为研究

在经典工程梁理论的基础上,结合张力薄膜的应力状态分析,提出充气机翼褶皱失稳的判据。计入表面薄膜褶皱引起的刚度退化效应,将机翼等效处理为一个变截面刚度的梁,建立了充气悬臂机翼的等效梁模型,并采用微分求积法进行充气机翼弯曲变形分析。计算结果与充气机翼的静力弯曲试验结果相吻合,验证了充气机翼弯曲变形分析方法的有效性。应用片条理论引入气动力模型,并与所建立的等效梁模型相耦合,建立充气机翼的静气动弹性耦合模型,并用迭代算法进行求解。考虑起皱和失稳两种判据,并计算获取了试验机翼的起皱动压和皱褶失稳动压形式,计算结果与风洞试验结果一致。根据所建立的充气机翼静气动弹性分析方法,可以预测充气机翼表面褶皱区的扩展和弯曲变形,进而绘制充气机翼的静气弹许用包线,为充气机翼的设计提供必要的安全边界参考依据。     

跨声速风扇叶片的静态气动弹性问题

使用时域的流固耦合数值计算方法,研究了跨声速风扇叶片在气动力和离心力共同作用下的静态气动弹性问题,分析了叶片在不同工况下的变形规律及叶片变形对整体气动性能的影响.NASA rotor 67的静态气动弹性计算说明气动力对叶片最大变形的贡献达13.07%, 而且叶片变形明显地改变了通道激波的位置和强度.宽弦空心跨声速风扇叶片的静态气动弹性计算说明叶片变形对总体气动效率的影响为0.15%~ 0.5%,其中气动力对变形贡献在叶片尖部的前缘可达41%,考虑气动力引起的变形使得该风扇的流量增大,气动特性线整体向右偏移.计算结果说明:气动力的非线性对跨声速风扇叶片静态变形问题有显著的影响,工程实践中从设计叶型到制造叶型的反扭过程应该采用流固耦合方法以得到更准确的叶型.      

NNW-FSI软件静气动弹性耦合加速策略设计与实现

在国家数值风洞（NNW）工程的资助下,依托NNW-FSI流固耦合模拟软件平台,从气动载荷作用下飞行器结构静变形大小与收敛过程无关的物理机制出发,基于变形增量叠加的方式,设计和实现了一种静气动弹性耦合加速策略,通过松弛因子对耦合迭代的收敛过程进行调整。结合超大展弦比无人机和CHN-T1模型两种不同外形,开展了不同松弛因子下的静气动弹性耦合数值模拟,对耦合加速策略的参数影响和加速效果进行了测试和评估。从计算误差控制角度对松弛因子加速耦合迭代收敛的作用机制进行了理论分析,弄清3种类型静气动弹性耦合模拟过程中松弛因子发挥的作用,并给出了松弛因子选取范围的建议。静气动弹性耦合模拟和理论分析结果表明,针对不同类型的静气动弹性耦合问题,选取合适的松弛因子,能够达到抑制振荡并加速收敛的效果。     

基于CFD和CSM耦合的通用静气弹分析方法

提出了一种适用于有限元精细化建模的流固耦合插值点选择方法,通过RBF(径向基函数)方法实现流固耦合面的数据交换,实现了基于CFD/CSM(computational fluid dynamics/computational structural mechanics)耦合的通用非线性静气弹分析方法。以HIRENASD(high Reynolds number aero-structural dynamics)风洞试验模型为验证对象,数值结果很好地与风洞试验结构变形、气动压力分布吻合,验证了所发展非线性CFD/CSM耦合静气弹求解器的精度。详细研究了HIRENASD模型在大迎角(AOA)流动下的静气动弹性特性,以及该模型弹性变形对机翼气动特性影响规律。研究表明:HIRENASD弹性模型变形后其升力小于刚性模型;在小迎角范围内刚性、弹性模型升力差随迎角增大呈线性增长;当迎角大于4°后,升力差先减小后基本保持不变,呈非线性关系。