Nonlinear dynamic analysis of a complex dual rotor-bearing system based on a novel model reduction method

In this paper, a dynamic model of a complex dual rotor-bearing system of an aero-engine is established based on the finite element method with three types of beam elements (rigid disc, cylindrical beam element and conical beam element), as well as taking into account the nonlinearities of all of the supporting rolling element bearings. To rapidly and accurately analyze dynamic behaviors of the complex dual rotor-bearing system, a two-level model order reduction (MOR) method is proposed by combining component mode synthesis (CMS) method and proper orthogonal decomposition (POD) technique. The first-level reduced-order model (ROM) of the dual rotors is obtained by CMS method with a high precision for the original system. Then, the POD method is applied to second-level model order reduction to further decrease the degrees of freedom (DOFs) of first-level ROM. Second-level ROM with mode expansion and direct second-level ROM are obtained, and the nonlinear displacement responses of the two ROMs are compared with the first-level ROM. The numerical results demonstrate that the proposed method has a higher computational efficiency and accuracy in terms of mode expansion than the direct model reduction by using POD method. In addition, the nonlinear vibration responses of the dual rotor-bearing system are studied by this second-level ROM in the case of different clearances of the inter-shaft bearing. The results indicate that the dynamic characteristics of the dual rotor-bearing system are very complicated for a large clearance.     

超支化烯丙基聚硅氮烷-巯基化合物体系紫外光固化研究

采用等温差示光量热扫描技术(DPC)研究了超支化烯丙基聚硅氮烷-巯基化合物体系的紫外光固化.对比了超支化聚合物-巯基化合物体系、二官能烯丙基硅氮烷-巯基化合物体系和超支化自聚体系的紫外光固化特性,结果表明超支化-巯基化合物体系可在引发剂浓度低(0.1wt%)、辐照强度低(5mW/cm2)的情况下迅速发生光固化反应;与超支化自聚反应相比固化速率高、双键转化率高;而与低官能-巯基化合物体系相比,由于超支化分子结构的独特性,导致固化速率快,双键转化率偏低.对引发剂浓度、辐照强度和反应温度对固化行为的影响表明,在引发剂浓度不超过0.5wt%和辐照强度不大于30.50mW/cm2时,反应速率分别与引发剂浓度和辐照强度的1/2次方成线性关系.运用带扩散因子的自催化固化动力学模型研究了不同温度下的固化行为,计算出特定条件下的光固化动力学参数,表观总反应级数约为8.76,表观活化能为13.97kJ/mol.     

Dynamic symmetry breaking and structure-preserving analysis on the longitudinal wave in an elastic rod with a variable cross-section

The longitudinal wave propagating in an elastic rod with a variable cross-section owns wide engineering background, in which the longitudinal wave dissipation determines some important performances of the slender structure. To reproduce the longitudinal wave dissipation effects on an elastic rod with a variable cross-section, a structure-preserving approach is developed based on the dynamic symmetry breaking theory. For the dynamic model controlling the longitudinal wave propagating in the elast...     

Accurate and straightforward symplectic approach for fracture analysis of fractional viscoelastic media

An accurate and straightforward symplectic method is presented for the fracture analysis of fractional two-dimensional(2D) viscoelastic media. The fractional Kelvin-Zener constitutive model is used to describe the time-dependent behavior of viscoelastic materials. Within the framework of symplectic elasticity, the governing equations in the Hamiltonian form for the frequency domain(s-domain) can be directly and rigorously calculated. In the s-domain, the analytical solutions of the displacement ...     

Dynamic analysis of geared transmission system for wind turbines with mixed aleatory and epistemic uncertainties

This paper deals with the co-existence of mixed aleatory and epistemic uncertainties in a wind turbine geared system for more reliable and robust vibration analyses. To this end, the regression-based polynomial chaos expansion(PCE) is used to track aleatory uncertainties, and the polynomial surrogate approach(PSA) is developed to treat the epistemic uncertainties. This non-intrusive dual-layer framework shares the same collocation pool, which is extracted from the Legendre series. Moreover, the ...     

Research on data assimilation strategy of turbulent separated flow over airfoil

In order to increase the accuracy of turbulence field reconstruction, this paper combines experimental observation and numerical simulation to develop and establish a data assimilation framework, and apply it to the study of S809 low-speed and high-angle airfoil flow. The method is based on the ensemble transform Kalman filter(ETKF) algorithm, which improves the disturbance strategy of the ensemble members and enhances the richness of the initial members by screening high flow field sensitivity ...     

基于双向流固耦合的超空泡射弹入水研究

超空泡射弹通过超空泡减阻技术在水下高速长距离航行, 是对抗水下近距离威胁的有效手段. 为了扩大防御范围、增加杀伤力, 超空泡射弹具有很高的发射速度. 高速超空泡射弹在入水时中受到极大的冲击载荷, 发生显著的结构变形, 结构变形与流场之间存在相互影响和作用, 常规的基于刚体假设的仿真研究方法不再适用. 为了研究高速超空泡射弹入水过程中的结构变形及其对流体动力特性的影响, 通过耦合流体力学求解器和结构动力学求解器, 建立了射弹高速入水双向流固耦合仿真模型, 并通过与文献中的试验结果进行对比验证了该模型空泡形态计算方法和耦合方法的准确性. 使用双向流固耦合的方法对高速射弹在不同初始攻角入水过程中的超空泡流动特性及结构变形特性进行了数值模拟研究, 通过对比流固耦合模型与刚体模型的计算结果, 得到了超空泡射弹的结构弯曲变形对流体动力载荷的影响. 研究结果表明: 高速射弹入水过程中流固耦合效应对超空泡流型及流体动力载荷的计算结果有显著影响; 本文所研究的射弹在考虑流固耦合效应, 带攻角垂直入水两倍弹长的范围内, 超空泡射弹的流体动力载荷与弯曲变形之间形成正反馈; 高速超空泡射弹在入水过程中受到的流体动力载荷及弹体应力应变随入水初始攻角的增加显著增大, 研究对象在初速1400m/s的条件下入水时, 当初始攻角不超过2°时不存在结构安全性问题.      

基于准点缺陷理论探索非晶合金蠕变机制

作为典型多体相互作用非平衡体系,如何明晰非晶合金多场耦合激励下变形机制,建立非晶合金变形行为、流动特性与微观结构特征本征关联始终是非晶合金力学性能的重要研究内容.本文以具有显著β弛豫行为的La_56.16Ce_14.04Ni_19.8Al₁₀非晶合金作为研究载体,通过开展宽温度应力窗口蠕变实验,着重考查了蠕变柔量、准稳态蠕变速率、特征弛豫时间、蠕变应力指数及蠕变激活能演变规律,系统研究了非晶合金蠕变行为与蠕变机制.基于准点缺陷理论分析了非晶合金蠕变行为由弹性向黏弹性及黏塑性逐步转变的过程,从微观结构演化角度构建了非晶合金蠕变行为完整物理图像.研究结果表明,非晶合金高温蠕变行为是一典型热力耦合、非线性过程,其潜在蠕变机制受控于温度、应力与加载时间.应力较低时,非晶合金蠕变机制对应于热激活单粒子流动.应力较高时,蠕变机制则对应于应力诱导局部剪切变形增强与温度诱导原子扩散等复杂耦合过程.非晶合金蠕变变形过程所涉及弹塑性转变源于非晶合金准点缺陷激活、微剪切畴经热力耦合激励形核长大、扩展与不可逆融合.     

k-ω SST湍流模式三维激波分离流动修正

“激波?边界层分离”是航空气动领域的典型湍流非平衡流动问题, 准确模拟激波分离对于跨声速飞行器气动性能评估和优化设计具有重要意义. 然而传统涡黏性湍流模式中涡黏性系数的定义方式并不适用于非平衡流动, k-ω SST湍流模式为此引入的Bradshaw假设在应用于三维强逆压梯度和较大分离流动时反而限制了雷诺应力的生成, 导致包括k-ω SST在内的常用涡黏性湍流模式均无法对此类流动进行准确模拟. 同时, 现有的非线性雷诺应力本构关系也并不能有效提高模拟精度. 为此, 针对k-ω SST模式分别提出了基于Bradshaw假设和基于长度尺度的两种激波分离流动修正方法. 前者通过提高Bradshaw常数的方式放宽了对雷诺应力生成的限制, 后者则从湍流长度尺度概念出发, 利用混合长度理论、湍动能生成/耗散之比和一种新定义的长度尺度之比构造了ω方程耗散项修正函数, 提高了模式在三维激波分离流动中的建模长度尺度. 两种方法对ONERA M6机翼跨声速大攻角流动均能得到较雷诺应力模式更好的模拟结果. 进一步的雷诺应力分析表明, 三维激波分离流动中“主雷诺应力分量”的概念不再成立, 各雷诺应力分量大小接近. 网格收敛性分析、对其他攻角状态的验证以及湍流平板边界层壁面律验证进一步确认了所提出的两种修正方法的合理性、有效性和通用性.      

Walsh函数有限体积法的多重网格特征研究

Walsh函数有限体积法(FVM-WBF)是一种能够在网格内部捕捉间断的新型数值方法. 持续增加Walsh基函数数目能够稳步提高FVM-WBF方法的求解分辨率, 但计算量暴发式增长和收敛速度下降的问题也会同步出现. 针对Walsh基函数数目增加而引起的计算效率问题, 本文分析了Walsh基函数及其系数所能影响的网格单元局部均值区域尺度, 发现其中隐含类似多重网格的尺度特征, 据此提出一种结合多重网格策略的FVM-WBF方法. 在定常流场计算中根据各级Walsh基函数影响尺度的不同, 对每级Walsh基函数设置满足其稳定性约束的时间步长, 在时间推进求解的过程中快速消除不同波长的数值误差, 实现多重网格的加速收敛效果. 选取NACA0012翼型和二维圆柱的定常无黏绕流问题作为算例, 对引入多重网格策略的FVM-WBF方法和不考虑多重网格策略的FVM-WBF方法进行对比测试. 数值结果证实: 新发展的FVM-WBF方法具备多重网格的关键特征, 在不增加任何特殊处理和计算量的情况下, 只需通过时间步长的调整, 就能够达到多重网格的加速效果, 显著提升计算效率.