电动力绳系周期解求解及稳定性控制

电动力绳系姿态摆动动力学方程具有非线性特点,且姿态摆动状态直接影响着空间任务能否正常进行。首先,基于哑铃模型及偶极子地磁场模型假设,推导了在倾斜圆轨道上运行的电动力绳系姿态动力学方程。其次,给出了利用摄动法和不动点迭代法相互配合的数值方法求解该非线性方程周期解的过程及解的稳定性判断准则,并针对不稳定周期解提出了相应的稳定性控制策略。最后,通过数值仿真分析,得到了姿态摆动运动的周期解图像,验证了周期解求解算法的有效性,并对不稳定周期解实现了稳定性控制。     

发生风雨激振拉索表面水膜形态数值分析和实测初探

结合Lemaitre斜拉索表面水膜风致形变解析算法,对数学模型中拉索运动方程进行修正,使其能计入拉索风致运动对表面附着水膜形态的影响,通过差分方法对方程进行数值求解,并进行参数分析:研究拉索表面水膜在重力、表面风压、表面摩擦力、表面张力作用下的形态变化,得到拉索发生风雨激振时表面水线的运动规律及形态变化.采用超声波探测仪对拉索表面水膜进行动态测量,实时获得拉索风雨激振形态下表面水膜动态响应数据,通过与水膜运动数值算法结果进行对比一定程度上探讨了圆截面斜拉索风雨振发生机理.     

一种惯导测试转台摩擦力模型及其参数辨识方法

针对转台运动系统中存在的摩擦力环节及其特性提出了一种静态摩擦力模型及其参数辨识方法,该模型克服了库伦摩擦模型的不连续性,并利用该模型对摩擦力在转台伺服系统中的影响,进行了仿真分析及试验对比。结果表明,该摩擦力模型可充分描述摩擦力对转台伺服系统运动性能的影响,为进一步提高系统的运动性能提供了一个研究基础。     

含间隙球轴承-转子系统的主共振研究

 为了深刻揭示球轴承-转子系统的共振机理,考虑了重力、装配等因素引起的非对称轴承间隙以及滚珠与滚道之间的Hertz接触,建立了球轴承-转子系统的运动方程,首次应用平均法求得了系统在主共振情况下的解析解,并讨论了系统参数对主共振的影响。根据奇异性理论方法计算得到了系统参数平面上的转迁集,研究了不同参数域内共振解的分岔模式,最后通过数值仿真验证了理论分析的正确性,从而为球轴承-转子系统的振动控制提供了理论依据。     

三维弹塑性接触问题的边界元法

将边界元法发展用于解决三维弹塑性接触问题,推导了无摩擦的三维弹塑性接触问题的积分方程,讨论了其数值解技术。数值结果与有限无法的数值结果对比表明,所发展的方法是成功的。     

在运动突风作用下亚声速机翼非定常气动力数值计算

本文应用运动涡有限基本解方法求得亚声速机翼在运动突风作用下非定常气动力的数值解。对于不同平面形状的机翼进行了计算。计算结果表明,该方法解类似的问题是有效的。     

六阶精度的群速度直接控制紧致格式及其应用

差分数值解产生非物理振荡的直接原因是在于非均一的波群的群速度,因此本文利用直接群速度控制的方法重新构造具有六阶精度的紧致型差分格式,以达到改善激波数值解的目的。对所构造的格式的精度及数值解随控制参数变化的行为进行了分析。最后文中给出了一些典型算例,证明了文中所构造的格式具有方法简便,物理含义清楚,精度高,网格基架点小,和捕捉激波能力较强的优点。      

径向简谐激励下航空发动机机匣的振动分析

根据航空发动机机匣薄壁结构的几何特征,将机匣简化为一端简支一端固支的圆柱壳结构,建立了系统的动力学方程,并采用多尺度法对动力学方程进行了求解,进而研究了径向简谐激励下机匣动力学方程的零解及非零解的稳定性,并讨论了系统单模态运动发展到复合模态运动以及能量在模态之间的转移过程,最后采用数值方法研究了径向简谐激励下机匣的振动特性.结果表明:随着径向简谐激励的增加,系统表现为周期运动、倍周期运动及混沌运动交替出现.      

多层手征介质圆波导的传输特性

阐述多层手征介质圆波导传输特性的递推特征函数解 ,导出了 N层手征介质圆波导传播常数的特征方程 ,提供了若干例子的数值结果     

轴向流中圆柱体的非线性动力学分析

通过引入流体摩擦力的二次阻尼项，推导出了轴向流中两端简支圆柱体的新的运动微分方程。对添加二次阻尼项的运动微分方程进行了无量纲化和Ritz-Galerkin离散。通过数值模拟，对两端铰支圆柱体的稳定性进行了系统的分析，并得出颤振发生的临界流速条件。     

适用于混合网格的约束最小二乘重构方法

为强化边界条件对计算流体力学(CFD)计算中流场重构过程的影响,基于三维任意多面体混合网格,对一种约束最小二乘法进行了研究。通过对不同边界类型的透明化处理,统一了约束方程组的形式,从而简化了该方法的应用。为了充分吸收消去法和加权法各自求解约束方程组的优点,提出了混合采用两种方法的思路。针对层流平板研究了加权法中加权系数如何取值问题,数值实验表明加权系数取到5基本已经足够大。最后通过亚声速层流平板和跨声速湍流ONERA-M6算例对比了混合法、加权法与原始最小二乘法,计算结果表明:约束最小二乘法相比原始最小二乘法的流场重构更加准确,特别是对边界质量较差网格的计算优势更为突出;混合消去加权法相比单独采用加权法的计算结果也有所改善。     

Robust wavelet neuro control for linear brushless motors

Design, simulation and experimental implementation of a wavelet basis function network learning controller for linear brushless dc motors (LBDCM) are considered. Stability robustness with position tracking is the primary concern. The proposed controller deals mainly with external disturbances, e.g. nonlinear friction force and payload variation in motion control of linear motors. It consists of two parts, one is a state feedback component, and the other one is a learning feedback component. The state feedback controller is designed on the basis of a simple linear model, and the learning feedback component is a wavelet neural controller. The attenuation effect of wavelet neural networks on friction force is first verified by the numerical method. The learning effect of wavelet neural networks on friction force is also shown in the numerical results. Then, a wavelet neural network is applied on a real LBDCM to on-line suppress the friction force, which may be variable due to the different lubrication. The effectiveness of the proposed control schemes is demonstrated by simulated and experimental results.     

二自由度扭转干摩擦系统的频域响应计算方法

针对一类作扭转振动的干摩擦系统,基于库仑阻尼模型和谐波平衡思想给出了其频域响应的求解方法,考虑了干摩擦界面上可能出现的周期内全滑移运动及阻滞-滑移运动.首先,建立了扭转干摩擦系统的二自由度集中质量模型,并进行了无量纲化.其次,建立了以谐波平衡思想处理该动力学模型的方法,并给出了系统的接触面作周期内全滑移运动及阻滞-滑移运动的判定条件.最后,对于这两种接触运动状态,分别给出了其摩擦力的处理方法和谐波展开形式,从而建立了扭转干摩擦系统的频域求解方法.以此方法计算了系统在一组不同的正压力下的响应,其结果与数值积分法所得的结果取得了相当的一致.此半解析方法较数值积分方法更能揭示扭转干摩擦系统的工作机理,并可以大幅提高计算时间,其推导和求解过程可用于类似的非线性系统.      

Evaluation of thermal effects on temperature-sensitive operating force of flow servo valve for fuel metering unit

The temperature-induced variation in operating force of flow control valves may result in performance degradation or even jam faults of fuel metering unit (FMU), which significantly affects the safety of aircrafts. In this work, an analytical modeling approach of temperature-sensitive operating-force of servo valve is proposed to investigate the temperature characteristics in varying temperature conditions. Considering the temperature effects, a new extended model of flow force is built and an analytical model of valve friction is also derived theoretically based on the dynamic clearance induced by thermal effects. The extremum condition of friction is obtained to analyze the characteristic-temperature points where jam faults occur easily. The numerical results show that flow force increases firstly and then decreases as temperature increases under a constant valve opening. The maximum friction of flow servo valve can be uniquely determined when the structural parameters and ambient temperature are given. The worst situation just happens at the characteristic-temperature points, which are linearly related to the axial temperature gradients of valve spool. Such evaluations may give an explanation for the temperature-induced jam faults of vulnerable valves and provide a reference for designers to determine a suitable working-temperature range of valves in practice.     

端柱面组合气膜密封系统稳态特性数值模拟与分析

针对端柱面组合气膜密封系统建立稳态特性的分析模型,给出基于有限元数值方法求解气膜压力分布的数值模拟方法和流程,实现端面、柱面区域的耦合求解,进而给出了端柱面组合气膜密封系统稳态特性数值模拟计算的方法.通过与Fluent软件仿真结果的比较验证数值模拟方法的正确性.重点研究膜厚与旋转轴转速对端柱面组合气膜密封系统稳态特性(端面气膜反力、柱面气膜反力、泄漏量、摩擦转矩)的影响,计算分析表明:随着膜厚的增加,泄漏量增加、摩擦转矩减小;随着端面膜厚的增加,端面气膜反力缓慢减小、柱面气膜反力出现极值;随着柱面膜厚的增加,端面气膜反力和柱面气膜反力都在减小;柱面膜厚对稳态特性的影响大于端面膜厚.端柱面结构中密封气膜压力的耦合对密封特性影响显著,稳态特性数值模拟方法的实现为动态特性分析提供了基础.      

基于库伦摩擦效应的箔片轴承极限承载力理论与试验

基于有限单元法建立了考虑库伦摩擦的波箔型径向气体箔片轴承的箔片结构模型,采用有限差分法和有限单元法耦合求解Reynolds方程和气膜厚度方程,通过求解轴颈达到极限偏心率时的轴承极限承载力,研究了箔片结构库伦摩擦效应对轴承极限承载力的影响规律,并搭建了轴承极限承载力测试试验台,利用温度法测量了两个具有不同轴承壳内表面粗糙度的波箔型径向气体箔片轴承的轴承极限承载力.通过对比分析仿真结果与试验结果表明:轴承壳圆柱孔内表面粗糙度为0.4μm的轴承在10000r/min和20000r/min下,轴承极限承载力分别为15.5N和42.3N;而表面粗糙度为1.6μm的轴承极限承载力为10.9N和29.6N,这是由于波纹箔片和轴承壳体之间的库伦摩擦力增大了波纹箔片的刚度,因此增大箔片结构摩擦因数使得轴承极限承载力降低,并且仿真结果变化趋势和试验结果变化趋势吻合.      

有摩擦弹性接触问题的增量数值解法

 指出应用全量变分法于有摩擦弹性接触问题是不能普遍适用的，提出了这一类问题增量变分法的数值解法和适用于增量法的库伦摩擦力的正则条件。对可转变的边界条件限于较简单的有限个离散点的情况，可采用状态试凑法来确定可变约束条件的状态。数值计算结果表明了本文提出的增量方法的正确性。     

基于接触刚度分布特性的叶片振动响应分析

根据接触刚度在锯齿形叶冠结合面上非均匀分布的特性,提出一种基于定义和有限元计算相结合的接触刚度计算方法。在此基础上,将微-宏滑动摩擦模型作为叶冠结合面处的摩擦模型,推导微滑动和宏滑动状态下摩擦力表达式。利用谐波平衡法将非线性摩擦力转化为等效刚度和等效阻尼进行振动特性分析。针对叶冠结合面相对位移幅值动态变化的特点,提出了一种迭代求解的振动响应分析方法。与文献提供的带摩擦阻尼结构悬臂梁振动响应实验数据相比,共振峰附近的振幅误差为1.99 mm,相对误差为3.9%。振幅的最大误差为5.53 mm,出现在远离共振峰的位置,验证了响应分析方法的可行性。将该方法应用于带冠叶片上,结果表明当激振力频率为812.3 Hz时,振幅为0.56 mm。     

具有非线性摩擦阻尼随机失谐的叶盘系统响应特性

针对目前含非线性摩擦阻尼失谐叶盘系统振动局部化机理研究存在的诸多困难,尤其是其非线性力学模型的建立和求解等难题.建立了具有滞迟干摩擦阻尼随机失谐的叶盘系统力学模型,并利用增量谐波平衡法分析了叶盘系统干摩擦力失谐度、耦合强度、黏性阻尼比和摩擦力强度等系统参数对叶盘系统受迫响应的影响规律.研究表明:含非线性摩擦阻尼的谐调叶...     

箔片结构库伦摩擦效应对径向箔片轴承特性的影响

基于有限元法建立波箔型气体径向轴承箔片结构的库伦摩擦模型,通过改变平箔片与波纹箔片之间以及波纹箔片与轴承座之间的摩擦因数,对比分析了在各种载荷分布条件下波纹箔片库伦摩擦模型与文献中线性弹簧模型的刚度特性,研究了库伦摩擦效应对波纹箔片刚度特性的影响规律.在此基础上,运用有限单元法和有限差分法求解雷诺方程和气膜厚度方程,研究了在两个工作转速下气体波箔片轴承中截面处最小气膜厚度随轴承承载力的变化规律以及承载力随偏心率的变化规律.通过数值仿真对该模型、文献中线性弹簧模型和刚性表面气体轴承进行对比分析,并把气膜厚度分布与文献结果进行了对比.结果表明:箔片的库伦摩擦力在一定程度上增大了波纹箔片的刚度,并且随着摩擦因数的增大其刚度以及两端固定的波纹箔片个数也增加,使得箔片轴承表面变“刚”,因此轴承静特性更趋于刚性表面轴承,此外当轴承承载力一定时,箔片摩擦因数越大轴承的最小气膜厚度越小.