三角形单元协调与非协调位移的能量正交关系

基于组合稳定化变分原理,周天孝提出的组合杂交法是绝对收敛和稳定的,它给出了一种系统化的增强应力/应变方法,并建立了一簇低阶仿射等价的n-cube(n=2,3)单元。本文论证了单元上应力插值为线性,位移插值为协调线性部分和非协调二次部分之和的三角形组合杂交单元其协调部分与非协调部分的能量正交关系,进而得到此三角形单元刚度矩阵等同于协调的三角形线性元刚度矩阵,即非协调部分无应变增强特性。     

基于混合网格的低声爆反设计方法研究

低声爆设计方法已成为新一代军民用超声速飞机研制过程中必须解决的关键难题之一。针对传统SGD低声爆外形反设计方法无法对声爆近场非线性效应进行描述和分析的缺点,提出了利用CFD方法求解得到的声爆近场压力分布代替F函数进行低声爆反设计的方法。声爆近场预测采用点-点对接的结构/非结构混合网格,充分利用非结构网格对复杂外形适应性强和结构化网格计算效率高的优点。结果分析表明,基于改进后的低声爆反设计方法得到的方案在声爆超压以及感觉噪声级等方面都比基于原始SGD方法得到的方案有较大改善。     

不同润湿性纳米通道内库埃特流动的模拟

利用非平衡分子动力学模拟方法, 模拟了两无限大平行平板组成的纳米通道内的库埃特流动, 并给出了壁面润湿性和速度对流场密度、速度分布及壁面滑移的影响规律.数值模拟中, 统计系综采用微正则系综, 势能函数选用LJ/126模型, 壁面设为刚性原子壁面, 温度校正使用速度定标法, 牛顿运动方程的求解则采用文莱特算法.结果表明, 纳米通道内流体密度呈对称的衰减振荡分布, 且随壁面润湿性的降低, 振荡幅度减小, 振荡周期保持不变;滑移量随壁面润湿性的提高而降低, 甚至在亲水壁面时出现负滑移现象;随壁面速度的增加滑移速度逐渐增大, 且在流体呈现非线性流动阶段其增幅显著加大.另外, 还发现当壁面设置为超疏水性时, 壁面滑移呈现出随润湿性降低而减小的反常现象, 并基于杨氏方程对其进行了解释.     

基本变量区域重要性测度及其稀疏网格解

为了提高现有基本变量对样本均值贡献的区域重要性测度指标的稳定性和收敛性, 提出了一个新的衡量基本变量内部各个区域对输出均值影响的重要性测度指标.并将其进一步扩展提出了一个衡量基本变量内部各个区域对输出总方差分解式中一阶方差影响的区域重要性测度指标.分析了所提指标的性质, 并探讨了它们与现有基本变量对样本均值贡献区域重要性测度指标和对样本方差贡献的区域重要性测度 指标之间的关系. 另外, 针对所提指标的特点, 还建立了其求解高效的稀疏网格积分法.算例结果表明, 所提新的基本变量对输出均值贡献的区域重要性测度指标不仅继承了现有指标的优点, 而且比现有指标具有更高的收敛性和稳定性.所提基本变量对一阶方差贡献的区域重要性指标能够在基本变量对样本方差贡献区域重要性测度的基础上, 进一步提供基本变量内部各个区域对总方差的一阶分量的影响信息.而所建稀疏网格积分法可以在保证计算精度的同时大幅度提高基本变量区域重要性分析的效率.     

基于损伤力学方法的铆接结构疲劳寿命预估

将损伤力学-有限元方法应用于飞机结构疲劳寿命预估中,对机翼蒙皮压窝形式的铆接结构进行了寿命计算和分析,得到了结构的疲劳寿命和疲劳损伤演化过程,并对铆接试验件进行了疲劳试验,将损伤力学方法的预测结果与试验结果进行了比较。分析结果反映出结构疲劳损伤累积初期较为缓慢、后期速率急剧上升的过程,符合损伤演化规律。在180MPa载荷水平下,计算寿命与试验平均寿命的相对误差为11.3%,满足工程精度要求。从而验证了损伤力学方法应用于实际工程结构寿命预测的可行性和可靠性。     

DH36钢拉伸塑性流动特性及本构关系

对DH36钢在温度从293~800 K、应变率为0.001和0.1 s-1的拉伸塑性流动特性进行实验研究,通过端口形貌图对变形前后的试样进行了微观分析,结果表明：(1)在实验温度范围内,0.001和0.1 s-1的应变率下,第三型应变时效现象出现,随应变率的增加,时效发生的温度区域移向更高温度;(2)第三型应变时效的发生与合金原子在晶界和晶粒中大量的第二相析出强化有关联;(3)建立包含第三型应变时效现象的统一本构模型,通过比较该模型能够较好的预测DH36的塑性拉伸流动应力。     

THE EVOLUTION CHARACTERIZATION OF LIBRATIN POINT ORBITS AND APPLICATION RESEARCH

三体问题中, 轨道的受力和运动规律非常复杂. 对于特定的任务, 如何选择轨道的初始解是一大难题.针对平面三体问题, 利 用近拱点庞加莱映射, 对平动点顺行轨道和逆行轨道的长期和短期演化性质进行分析.根据轨道的初始状态将其分为逃逸轨道和捕获轨道.对于逃逸轨道, 给出了同宿轨道和异宿轨道的设计方法, 并利用两级微分修正法消除了拼接点处的位置不连续问题.对于捕获轨道, 得到了几类典型的周期和准周期轨道.对逆行轨道的演化性质进行分析时发现, 逆行轨道通常为准周期轨道, 比顺行轨道更加稳定.利用近拱点庞加莱映射可以快速确定不同类型轨道对应的初始状态, 为特定任务需求下的轨道设计提供了一种快速而有效的选择方案.     

机翼双梁模型的动力学修正及应用

由于全尺寸机翼有限元模型太过复杂,不能用于大型客机翼吊发动机振动载荷的传递特性研究,本文提出了一个全新的机翼双梁动力学模型概念。基于机翼双梁模型的主要模态对其结构刚度不确定参数的灵敏度分析,选取翼梁及翼肋的水平、垂直、扭转刚度参数作为有效设计变量对双梁模型进行了动力学修正,使其与全尺寸机翼有限元参考模型前6阶模态频率的最大误差小于11.11%,大大提升了机翼双梁模型的动力符合性;集成建立了大型客机“吊架-机翼-机身”全机动力有限元模型,并计算了巡航状态下发动机振动载荷经机翼传递到机身各框段的动载荷分布特性。分析结果表明：机翼双梁模型较好地满足了发动机振动载荷传递及振动传递路径分析研究需求,为民机舱内声学预计及发动机隔振安装设计提供了基础数据。     

气泡弹性对同心筒水下发射影响研究

潜射导弹自同心简装置发射过程中在筒口处会形成燃气泡,由于气泡在水下的弹性效应,水气耦合作用导致气泡周期性地膨胀压缩,会对导弹弹道参数、受载特性产生一定影响.以CFD为技术手段,建立了导弹水下发射动态仿真模型,通过数值模拟获得了发射过程中多相流场、导弹弹道参数及受力的时间历程曲线,分析了气体弹性效应对参数变化的影响;并针对发射深度及其发射艇速对气泡弹性效应的影响进行了仿真分析,给出了其影响规律.仿真方法和结果对工程研究有一定参考价值.     

两种湍流模型时域颤振计算方法研究

采用时域计算分析方法进行了机翼跨音速颤振特性研究。在结构运动网格的基础上,采用格点格式有限体积方法进行空间离散和双时间全隐式方法进行时间推进求解雷诺平均N-S方程。针对流动粘性分别应用了SST湍流模型和SSG雷诺应力模型,通过对跨音速标模算例AGARD445.6机翼的计算结果与实验值的对比分析,其中应用SST湍流模型得到的颤振速度与实验值最为接近,特别是在跨音速段平均相对误差在3%以内;并且计算结果整体上反映了跨音速颤振"凹坑"物理特性,验证了方法的有效性。