基于ROM技术的阵风响应分析方法

阵风响应分析是大型飞机设计过程中必不可少的环节. 现有的阵风响应分析主要采用基于线化升力面理论的气动力模型,不能考虑到各种非线性效应,不适合于跨音速气动弹性的分析. 基于CFD技术,采用系统辨识方法,在状态空间内建立了降阶的非定常气动力模型(reduced order model, ROM). 耦合结构运动方程、非定常气动力模型(结构运动)、外激阵风的气动力模型,建立了基于CFD技术的阵风响应分析模型.算例研究了某一典型机翼在方波激励下的阵风响应问题,对比了各阶模态位移的响应以及翼根弯矩的响应. 基于ROM技术的计算结果与CFD/CSD直接耦合仿真结果吻合,证明了该方法的正确性和精度.      

CFD结合降阶模型预测阵风响应

传统的阵风响应主要在频域内进行分析,气动载荷基于线性方法计算,不能考虑黏性和跨声速流动影响. 飞机设计需考虑不同频率和不同形状阵风的响应,基于CFD的阵风响应预测由于计算工况太多,工作量巨大. 本文发展了一种CFD结合非定常气动力ARMA(autoregressive-moving-averagemodel)降阶模型的阵风响应分析方法,CFD只要针对给定频率和形状的一种阵风响应进行计算,对获得的气动力时间历程运用线性最小二乘法参数辨识ARMA降阶模型的系数,则对任意频率和形状的阵风,代入降阶模型即可确定该阵风的响应,大大提高了计算效率. 为效验发展的方法,先计算NACA0012翼型在低马赫数0.11的阵风响应,通过对比CFD、ARMA降阶模型及早期发展的不可压阵风响应预测方法的结果,验证了方法的有效性. 再对比CFD、ARMA在跨声速马赫数0.8的阵风响应预测结果,证实所发展的方法对跨声速阵风响应预测亦是有效的.      

基于传递函数法的大展弦比机翼阵风响应分析

将传递函数法应用于大展弦比机翼的阵风响应分析。首先,基于二元机翼的运动方程和准定常片条理论建立机翼的阵风响应微分方程,对其进行Laplace变换,并转换为状态空间方程形式。然后,运用传递函数方法,获得机翼响应在频域的解析解,通过Laplace数值逆变换求得机翼在时域内的响应。通过与已有文献结果对比,验证了本文方法的正确性。最后,采用该方法求解了"1-cos"型阵风和连续大气湍流作用下的机翼响应,并对结果进行了分析讨论。     

基于非定常N-S方程的三维机翼颤振特性研究

本研究采用N-S方程为控制方程,求解三维机翼的非定常气动力,并与颤振方程耦合迭代计算,求解机翼广义坐标响应的时间历程,根据广义坐标的时间历程的衰减、等幅和发散振荡等情况确定机翼颤振临界条件.以某导弹翼面作为算例,计算结果与实验结果和理论分析相吻合.     

ANALYSIS OF THE HIGH-ASPECT-RATIO WING GUST RESPONSE BY TRANSFER FUNCTION METHOD1)

将传递函数法应用于大展弦比机翼的阵风响应分析。首先,基于二元机翼的运动方程和准定常片条理论建立机翼的阵风响应微分方程,对其进行Laplace变换,并转换为状态空间方程形式。然后,运用传递函数方法,获得机翼响应在频域的解析解,通过Laplace数值逆变换求得机翼在时域内的响应。通过与已有文献结果对比,验证了本文方法的正确性。最后,采用该方法求解了“1-cos”型阵风和连续大气湍流作用下的机翼响应,并对结果进行了分析讨论。     

考虑粘弹性材料频变特性的复合结构频率响应分析

论文采用基于Layerwise离散层理论的有限元板单元建模粘弹性阻尼复合结构,发展了一种考虑粘弹性材料的频变特性的复合结构频率响应计算方法.该方法首先计算结构的质量矩阵和各个频率点下的结构复刚度矩阵,然后求解运动方程的若干阶复特征对并由此计算各阶模态的传递函数,最后将各模态传函线性叠加得到近似的总传函.为了提高计算效率,论文采用了一种高效的数值方法,即只计算若干频率点下的特征向量与特征值,并计算这些点处特征向量关于频率的高阶导数,通过泰勒系数展开逼近和Rayleigh商式,可求得附近若干频率点处的特征向量和特征值,从而避免了在各个频率点下求解大自由度结构特征方程的问题,可以极大地提高计算效率.对一端简支的三层约束阻尼梁算例进行了分析,并与文献中的结果作对比,结果验证了方法的有效性和计算效率.     

弹性连接旋转柔性梁动力学分析

采用Chebyshev 谱方法对考虑根部连接弹性的平面内旋转柔性梁动力学特性进行研究. 基于Gauss-Lobatto 节点与Chebyshev 多项式方法对柔性梁变形场进行离散,通过投影矩阵法施加固定及弹性连接边界条件. 利用Chebyshev 谱方法获得了系统固有频率和模态振型数值解,通过与有限元方法及加权残余法的比较,验证了方法的有效性. 分析了弹性连接刚度、角速度比率、系统径长比及梁的长细比等参数对系统固有频率及模态振型的影响. 研究发现:由于系统弯曲模态、拉伸模态的频率随各参数的变化规律不一致,将出现频率转向与振型转换现象;随着弹性连接刚度、角速度比率及系统径长比的增大,低阶弯曲模态频率增大并超过高阶拉伸模态频率,随着梁的长细比的增大,低阶拉伸模态频率增大并超过高阶弯曲模态频率.      

某飞机机翼颤振模型模态测试及分析

以某飞机机翼颤振缩比模型为研究对象,通过纯模态测试得到试验件前九阶模态参数。建立机翼梁架有限元模型并进行计算,得到各阶模态振型及频率。将各阶模态的计算结果与试验结果进行对比,频率的计算结果与纯模态试验值的误差均在5%以内,主要模态的振型计算结果与试验结果也基本一致,验证了有限元计算模型的准确性。在此基础上可以进行下一步的颤振分析,且试验结果可以作为机翼颤振模型风洞试验的参考依据。     

基于动网格技术的柔性后缘自适应机翼气动特性分析

研究了带柔性后缘的可变弯度自适应机翼在自适应变弯度过程中的气动特性.自适应变弯度过程中的气动力计算采用了基于弹簧理论的非结构动网格技术,求解NS方程时采用有限体积的二阶迎风格式离散,时间推进为隐格式双时间推进方法.通过计算柔性后缘机翼的升力特性、阻力特性及升阻比特性,并与带刚性后缘机翼的气动特性进行比较,发现柔性后缘机翼在后缘偏转时,其最大升阻比对应的迎角随着偏转角增大而降低.在中等迎角及接近失速迎角情况下,柔性后缘机翼升力系数明显优于刚性后缘机翼,并且其升力线变化较为平缓,有效迎角范围更大.     

黏弹性地基上功能梯度材料板的振动分析

为研究黏弹性地基上功能梯度材料板的自由和强迫振动特性,基于Reddy高阶剪切变形理论以及由Shen导得的广义Karman型方程,用双重Fourier级数法推导了三参数黏弹性地基上四边简支功能梯度材料板自由振动和动力响应的解析解,计算了各模态自振频率和半波冲击载荷作用下的动力响应,讨论了材料组分指数、黏弹性地基参数、边厚比等因素对自由振动和动力响应的影响.结果表明,黏弹性地基的剪切和压缩刚度显著提升了功能梯度材料板的振动频率,减小了动力响应;另外,地基的黏性对振动频率和动力响应也有一定的影响.     

The 7th International Conference on Fluid Mechanics May 24-27,2015,Qingdao,China

正http://www.icfm7.org First Announcement and Call for PapersThe objective of International Conference on Fluid Mechanics(ICFM)is to provide a forum for researchers to exchange new ideas and recent advances in the fields of theoretical,experimental,computational Fluid Mechanics as well as interdisciplinary subjects.It was successfully convened by the Chinese Society of Theoretical and Applied Mechanics(CSTAM)in Beijing(1987,     

INFORMATION FOR CONTRIBUTORS

Contributions： The Journal, Acta Mechanica Solida Sinica, is pleased to receive papers from engineers and scientists working in various aspects of solid mechanics. All contributions are subject to critical review prior to acceptance and publication.     

Preface

This special issue of PARTICUOLOGY is devoted to the first UK-China Particle Technology Forum taking place in Leeds, UK, on 1-3 April 2007. The forum was initiated by a number of UK and Chinese leading academics and organised by the University of Leeds in collaboration with Chinese Society of Particuology, Particle Technology Subject Group （PTSG） of the Institution of Chemical Engineers （IChemE）, Particle Characterisation Interest Group （PCIG） of the Royal Society of Chemistry （RSC） and International Fine Particle Research Institute （IFPRI）. The forum was supported financially by the Engineering and Physics Sciences Research Council （EPSRC） of United Kingdom,     

基于鲁棒滤波的捷联导引头视线角速度估计方法

针对捷联导引头无法直接获取视线角速度等信息的问题,研究了鲁棒滤波在大气层外飞行器捷联导引头视线角速度估计中的应用。为了建立非线性滤波估计模型,考虑目标视线角速度的慢变特性,采用一阶马尔科夫模型建立了状态方程;推导了视线角速度的解耦模型,并建立了量测方程;考虑到实际应用中存在系统噪声统计特性失准的问题,基于Huber-Based鲁棒滤波方法,设计了视线角速度滤波器,并完成了基于Huber-Based滤波方法和扩展卡尔曼滤波方法的数学仿真。仿真结果表明Huber-Based滤波方法的视线角、视线角速度及视线角加速度估计精度分别达到0.1140'、0.1423'/s、0.0203'/s2,而扩展卡尔曼滤波方法的视线角、视线角速度及视线角加速度估计精度仅分别为0.6577'、0.6415'/s、0.0979'/s~2。仿真结果证明了该方法可以有效地估计出相对视线角速度等信息,并且在非高斯噪声的条件下,依然可获得较高的估计精度,具有一定的鲁棒性。     

Guide for authors

正Each of the sections below provides essential information for authors.We recommend that you take the time to read them before submitting a contribution to Acta Mechanica Sinica.We hope our guide to authors may help you navigate to the appropriate section.How to prepare a submission This document provides an outline of the editorial process involved in publishing a scientific paper in Acta Mechanica     

Use specification of multiscale materials for life spanned over macro-, micro-, nano-, and pico-scale

Multiscale material intends to enhance the strength and life of mechanical systems by matching the transmitted spatiotemporal energy distribution to the constituents at the different scale, say—macro, micro, nano, and pico,—, depending on the needs. Lower scale entities are, particularly, critical to small size systems. Large structures are less sensitive to microscopic effects. Scale shifting laws will be developed for relating test data from nano-, micro-, and macro-specimens. The benefit of reinforcement at the lower scale constituents needs to be justified at the macroscopic scale. Filling the void and space in regions of high energy density is considered.Material inhomogeneity interacts with specimen size. Their combined effect is non-equilibrium. Energy exchange between the environment and specimen becomes increasingly more significant as the specimen size is reduced. Perturbation of the operational conditions can further aggravate the situation. Scale transitional functions and/or f_j/j+1 are introduced to quantify these characteristics. They are represented, respectively, by , and (f_mi/ma,f_na/mi,f_pi/na). The abbreviations pi, na, mi, and ma refer to pico, nano, micro and macro.Local damage is assumed to initiate at a small scale, grows to a larger scale, and terminate at an even larger scale. The mechanism of energy absorption and dissipation will be introduced to develop a consistent book keeping system. Compaction of mass density for constituents of size 10⁻¹², 10⁻⁹, 10⁻⁶, 10⁻³ m, will be considered. Energy dissipation at all scales must be accounted for. Dissipations at the smaller scale must not only be included but they must abide by the same physical and mathematical interpretation, in order to avoid inconsistencies when making connections with those at the larger scale where dissipations are eminent.Three fundamental Problems I, II, and III are stated. They correspond to the commonly used service conditions. Reference is made to a Representative Tip (RT), the location where energy absorption and dissipation takes place. The RT can be a crack tip or a particle. At the larger size scales, RT can refer to a region. Scale shifting of results from the very small to the very large is needed to identify the benefit of using multiscale materials.