热弹性复合材料层合板的响应分析

利用混合变分原理,考虑了热平衡方程与导热方程中变量的对偶关系,通过增加非齐次正则方程的维数,成功地将热弹性体机一热耦合问题的非齐次Hamilton正则方程转化为能独立求解的齐次Hamilton正则方程．所以大大简化了问题的求解。并通过数值实例研究了在温度栽荷和动力栽荷作用下及考虑阻尼的四边简支复合材料层合板的响应问题。     

蜂窝夹层板动力分析的Hamilton体系与辛算法

将蜂窝夹层板动力分析从Lagrange体系改换为Hamilton体系。通过作者早已提出的一条简单而统一的新途径,建立了蜂窝夹层板动力学的相空间非传统Hamilton变分原理,并从该原理推导出相应的Hamilton正则方程、边界条件与初始条件。然后基于这种相空间非传统Hamilton变分原理,提出蜂窝夹层板动力响应分析的辛空间有限元-时间子域法,算例的数值结果表明,这种新方法的计算精度与效率都明显高于国际上常用的Wilson-θ法和Newmark-β法。     

混合小推力航天器轨道保持高性能滑模控制

针对采用太阳帆、太阳电混合小推力推进的航天器，研究了其在日心悬浮轨道的保持控制问题。为解决已有控制方法中未综合考虑内部未建模动态和外部未知扰动的问题，以及进一步提高系统控制性能，设计了一种高性能滑模控制策略。首先，考虑模型不确定性，建立了混合小推力航天器在日心悬浮轨道柱面坐标系的动力学方程；其次，基于改进型条件积分滑模面和径向基（RBF）神经网络设计了控制律，结合自适应方法在线估计不确定参数；接着，将求取的虚拟控制量在推进剂最优条件下转换成实际控制量，即太阳帆姿态角和太阳电推进力；最后，数值仿真验证了上述设计方法提高了系统鲁棒性，减小了轨道位置超调，并且混合推进相比于单一太阳帆推进，在更短收敛时间内控制精度提高了4个数量级，相比于单一太阳电推进，一年可以节省约89.6%的推进剂。     

压电热弹性体开口壳的响应分析

根据广义的Hamilton变分原理推导出了压电热弹性体非齐次的Hamilton正则方程。结合压电热弹性体平衡方程和热平衡方程,成功地导出了压电热弹性体机、电、热耦合问题的齐次状态方程。将非齐次方程转化为齐次方程不仅使问题变得大为简化,同时也减少了数值计算的工作量。数值实例研究了温度载荷和力载荷作用下压电热弹性材料四边简支层合开口壳的响应问题。     

微纳米压痕测试技术:发展与应用

近年来,具有高精度、高通量的微纳米压痕测试技术,已被广泛应用于研究微/纳米尺度下材料力学性能演化规律和变形行为中。然而,在航空航天材料试验测试领域,令研究人员更感兴趣的往往是如何更好地揭示材料工程性能,更好地理解材料在服役环境下变形损伤机制。因此,接近材料真实服役环境（如高/低温、电/磁场）下的微纳米压痕测试系统更具应用潜力。首先对传统的微纳米压痕测试技术进行回顾总结,涉及测试系统的组成、经典分析理论方法及其面临的尺度/尺寸效应。然后,简要描述典型磁电弹性材料在力-电-热-磁多场耦合环境下接触力学行为的解析模型,并着重阐述面向材料实际服役环境下的压痕测试技术的典型应用,包括高/低温纳米压痕测试和电/磁场耦合条件下的纳米压痕测试应用。最后,讨论了目前发展所面临的主要问题和挑战,这对微纳米压痕测试技术的进一步发展和先进应用具有重要意义。     

振动叶栅中尾迹/激波相干的数值研究

研究了尾迹/激波相干对振动叶栅非定常流动的影响.数值方法利用了先进的谐波方法.该方法基于傅立叶变换,把叶栅非定常流变量分解为一时间平均量和若干谐波分量,继而,时域非定常N-S方程被转化为频域定常傅立叶系数方程.定常傅立叶系数方程的求解采用了伪时间推进方法.为加速数值解收敛,采用了本地时间步长和多重网格技术.通过NASA67转子算例结果表明:上游尾迹对振动叶栅非定常流动具有显著影响.当存在尾迹相干时,叶片通道内的激波在轴向方向上表现出较大的漂移运动,叶片所受的非定常力幅值变大,同时失去了时间周期性.      

Structural Nonlinear Flutter Characteristics Analysis for an Actuator-fin System with Dynamic Stiffness

The flutter characteristics of an actuator-fin system are investigated with structural nonlinearity and dynamic stiffness of the electric motor. The component mode substitution method is used to establish the nonlinear governing equations in time domain and frequency domain based on the fundamental dynamic equations of the electric motor and decelerator. The existing describing function method and a proposed iterative method are used to obtain the flutter characteristics containing preload freeplay nonlinea...     

非线性数学规划在飞行力学中的应用

主要介绍非线性数学规划在飞行器轨迹优化中的应用，处理方法是将每一个状态变量用三次多项式表示，状态微分方程就化成一组代数方程，从而使飞机轨迹优化问题化成了离散的数学规划问题，再用广义Ｌａｇｒａｎｇｅ乘子法求解，该方法对一般的轨迹优化问题易于编程实现。     

电动膨胀循环液氧甲烷变推力发动机动态特性研究

电动泵压式液体火箭发动机受到了广泛的关注，然而电池有限的输出功率和过于沉重的质量成为限制电动泵压式发动机发展的重要因素。为此本文提出了一种电机驱动燃料泵和涡轮驱动氧泵的电动膨胀循环液体火箭发动机方案，并着重研究了该型发动机的动态响应特性。首先给出了20 kN级电动膨胀循环发动机的技术指标和部组件参数，进一步基于AMESim平台建立了全系统动力学模型，验证了方案的可行性和部组件动力学模型的准确性，并深入研究了单点工况和调节工况的动态响应特性。结果表明，针对启动过程而言，涡轮泵调整时间较电动泵长，这降低了系统响应速度，但工况越高，系统响应速度越快；高工况启动时，甲烷在冷却通道内的剧烈相变和跨临界状态的不连续物性相互耦合易引发系统振荡；就调节过程而言，推力调节时普遍存在超调或凹坑现象，且系统在两相同工况之间调节时，正调响应速度快于负调，这也导致阶跃幅值相等条件下的系统调整时间随目标工况升高而缩短。     

倾转旋翼机动态倾转过渡过程的操纵策略优化

利用最优控制方法研究倾转旋翼机的最优动态倾转过渡过程,并得到最优操纵策略,使得由时间、姿态角变化以及驾驶员工作负荷等组成的性能指标达到最小。首先,在基本纵向刚体飞行动力学模型的基础上引入混合操纵方程,并使用杆量位移的一阶导数作为控制量,形成适用于计算倾转旋翼机动态倾转过渡过程的飞行动力学模型,从而能在动态倾转操纵策略优化过程中考虑到操纵系统特性对操纵量变化速度的限制,以及避免操纵量在优化过程中出现跳跃不连续。然后,将倾转旋翼机的最优动态倾转过渡过程转化为非线性动态最优控制问题,建立合理的性能指标,并采用直接转换法和序列二次规划算法进行求解。最后,以XV-15倾转旋翼机为样机,分别计算正向和逆向最优动态倾转过渡过程,并与驾驶员飞行仿真数据进行对比。结果表明:飞行状态量的时间历程与文献吻合地较好,且俯仰姿态角和杆量位移变化更加柔和。最优控制方法可以用于研究倾转旋翼机的最优动态倾转过渡过程。     

基于非线性规划的空空导弹水平转弯轨迹优化

基于空空导弹的质点动力学方程建立了轨迹优化模型,并根据非线性规划优化算法的特点,对优化模型中的积分变量进行替换。采用一种改进的多重打靶法将轨迹优化问题转化为约束非线性规划问题,在此基础上采用广义乘子法和变尺度法求解满足各种约束条件的最优转弯轨迹。仿真结果表明,该轨迹优化方法实现了导弹在约束条件下的最短时间转弯,且所得结果符合运动学规律,对于离轴发射以及越肩发射的转弯问题具有一定意义。     

前飞状态旋翼动稳定性及自由度间相互作用研究

为研究直升机在前飞状态下的稳定性,建立了前飞状态下旋翼的结构模型、气动力模型和平衡方程组。对方程组的周期系数处理以及前飞时的相关参数对旋翼动稳定性的影响等进行了研究。在前飞动稳定性分析中引入动力入流自由度。从分析各个自由度间的相互作用入手揭示旋翼动不稳定性的物理本质及影响因素。      

Earth-moon Trajectory Optimization Using Solar Electric Propulsion

The optimization of the Earth-moon trajectory using solar electric propulsion is presented. A feasible method is proposed to optimize the transfer trajectory starting from a low Earth circular orbit (500 km altitude) to a low lunar circular orbit (200 km altitude). Due to the use of low-thrust solar electric propulsion, the entire transfer trajectory consists of hundreds or even thousands of orbital revolutions around the Earth and the moon. The Earth-orbit ascending (from low Earth orbit to high Earth orbit) and lunar descending (from high lunar orbit to low lunar orbit) trajectories in the presence of J2 perturbations and shadowing effect are computed by an analytic orbital averaging technique. A direct/indirect method is used to optimize the control steering for the trans-lunar trajectory segment, a segment from a high Earth orbit to a high lunar orbit, with a fixed thrust-coast-thrust engine sequence. For the trans-lunar trajectory segment, the equations of motion are expressed in the inertial coordinates about the Earth and the moon using a set of nonsingular equinoctial elements inclusive of the gravitational forces of the sun, the Earth, and the moon. By way of the analytic orbital averaging technique and the direct/indirect method, the Earth-moon transfer problem is converted to a parameter optimization problem, and the entire transfer trajectory is formulated and optimized in the form of a single nonlinear optimization problem with a small number of variables and constraints. Finally, an example of an Earth-moon transfer trajectory using solar electric propulsion is demonstrated.     

基于动态吸振器的高超声速复合材料壁板颤振抑制及其优化设计

研究了动态吸振器对高超声速流中复合材料壁板颤振的抑制作用.将壁板和动态吸振器的相互作用力描述为两者相对位移和相对速度的函数,利用von-Karman非线性应变-位移关系和3阶非线性活塞理论,根据哈密顿原理和牛顿第二定律分别建立了壁板和动态吸振器的运动微分方程.采用模态假设法对系统进行离散并进行数值模拟,对动态吸振器的安装位置进行了优化设计.数值计算结果表明:合理选择动态吸振器参数可使壁板临界颤振动压提高51.7%,同时可降低颤振后极限环振动的幅值.      

基于传递函数方法的局部覆盖环状CLD圆柱壳动力学分析

 采用传递函数方法对局部覆盖环状约束层阻尼（CLD）圆柱壳进行了动力学分析。基于Donnell薄壳简化,由Hamilton原理导出了局部覆盖CLD圆柱壳的运动方程和边界条件,通过Laplace变换并引入状态向量后,建立了系统的状态空间方程,利用传递函数方法求解方程得到了系统的固有频率、〖JP2〗损耗因子和频响曲线。算例分析结果验证了该方法的有效性和精确性,且更适于处理黏弹性材料具有频变特性的问题。讨论了CLD覆盖率和黏弹层厚度对系统固有频率和损耗因子的影响,为局部覆盖CLD圆柱壳的优化设计奠定了基础。     

Characteristics of motorized spindle supported by active magnetic bearings

A motorized spindle supported by active magnetic bearings（AMBs） is generally used for ultra-high-speed machining. Iron loss of radial AMB is very great owing to high rotation speed, and it will cause severe thermal deformation. The problem is particularly serious on the occasion of large power application, such as all electric aero-engine. In this study, a prototype motorized spindle supported by five degree-of-freedom AMBs is developed. Homopolar and heteropolar AMBs are independently adopted as radial bearings. The influences of the two types of radial AMBs on the dynamic characteristics of the motorized spindle are comparatively investigated by theoretical analysis, test modal analysis and actual operation of the system. The iron loss of the two types of radial AMBs is analyzed by finite element software and verified through run-down experiments of the system. The results show that the structures of AMB have less influence on the dynamic characteristics of the motorized spindle. However, the homopolar structure can effectively reduce the iron loss of the radial AMB and it is useful for improving the overall performance of the motorized spindle.     

采用微扰法测试碳泡沫的电磁参数

利用网络分析仪连接矩形谐振腔组成材料电磁参数微扰法测试系统。由谐振腔微扰理论出发,得到以托管装载小样品测试时复介电常数和复磁导率的计算公式。对碳泡沫进行电磁参数测试。结果表明,随着碳泡沫电导率的增加,其介电常数缓慢增加,电损耗先增加后减小,而磁损耗不断增加。碳泡沫具有介电常数小、电损耗大、随电导率增加而增加的非本征磁损耗,这种独特的电磁损耗特征使得碳泡沫有用作宽频带高温结构吸波材料的可能性。     

同步发电机动态性能的数字机研究

本文采用P变换阵将同步发电机的时变动态方程,转换成所谓尸派克方程。由此建立了单台发电机的状态方程。在航空无刷交流发电机中,由于调压器的非线性作用及交流励磁机的延滞,使主同步发电机的激励函数U_F（t）可归纳为带滞后（t_E）的指数曲线(等效时间常数T_ff)。 用数字机对上述状态方程按解析法、四阶龙格库塔法、指数矩阵(e^AT)法、网络拓扑法等四种方法,分别打印出八种同步发电机动态曲线。分析了四种算法的计算精度及稳定区域,确定其应用范围,并得出一定的结论。     

直升机桨叶扬起下坠碰撞动响应计算

针对铰接式旋翼桨叶扬起下坠碰撞问题,采用中等变形梁模型,同时引入有限转角模拟桨叶绕铰的刚性运动。运用Hamilton原理建立桨叶扬起下坠碰撞动力学方程,采用隐式Newmark法计算桨叶在有初始挥舞角下的下坠碰撞响应。计算结果与试验数据的对比表明:如果采用定步长的方法计算响应,计算时间较长,计算结果会产生明显的延迟。为了解决上述问题,在碰撞前后设置了一碰撞缓冲区,采用了变步长的方法强制计算点在跳过此缓冲区时必须进入该区。此法在计算桨叶扬起下坠碰撞动响应时,能较好的解决计算延迟问题,计算时间大大减少,计算结果与试验数据吻合较好。     

基于欧拉方程的二维振荡机翼非定常气动设计反命题方法

建立了基于欧拉方程组的二维振荡机翼非定常气动设计反命题方法的数学模型.通过一系列变换,将物理时空中的求解域转换成映象坐标系中的规范区域;并导出了映象坐标下的欧拉方程,结合非定常反命题的边界条件,便可用有限差分方法求解.本方法引入映象坐标系解决反命题边界形状的不确定性,并能利用欧拉方程现有的各种差分格式于反命题求解.