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

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

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

传统的阵风响应主要在频域内进行分析,气动载荷基于线性方法计算,不能考虑黏性和跨声速流动影响.飞机设计需考虑不同频率和不同形状阵风的响应,基于CFD的阵风响应预测由于计算工况太多,工作量巨大.本文发展了一种CFD结合非定常气动力ARMA(autoregressive-moving-average model)降阶模型的阵风响应分析方法,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级数法推导了三参数黏弹性地基上四边简支功能梯度材料板自由振动和动力响应的解析解,计算了各模态自振频率和半波冲击载荷作用下的动力响应,讨论了材料组分指数、黏弹性地基参数、边厚比等因素对自由振动和动力响应的影响.结果表明,黏弹性地基的剪切和压缩刚度显著提升了功能梯度材料板的振动频率,减小了动力响应;另外,地基的黏性对振动频率和动力响应也有一定的影响.     

一种基于直接计算高阶奇异积分的断裂力学双边界积分方程分析法

相对于有限元法,边界单元法在求解断裂问题上有着独特的优势,现有的边界单元法中主要有子区域法和双边界积分方程法.采用一种改进的双边界积分方程法求解二维、三维断裂问题的应力强度因子,对非裂纹边界采用传统的位移边界积分方程,只需对裂纹面中的一面采用面力边界积分方程,并以裂纹间断位移为未知量直接用于计算应力强度因子.采用一种高阶奇异积分的直接法计算面力边界积分方程中的超强奇异积分;对于裂纹尖端单元,提供了三种不同形式的间断位移插值函数,采用两点公式计算应力强度因子.给出了多个具体的算例,与现存的精确解或参考解对比,可得到高精度的计算结果.      

汶川地震诱发文家沟巨型滑坡 碎屑流基本特征及成因机制初步分析


2008年5月12日,汶川M80地震在四川省绵竹市清平乡文家沟内诱发一巨型滑坡。通过现场调查得知,滑坡前后缘高差455m,厚度20～30m,滑面为基岩层面,初始方量2750×107m3。滑体在运动中转化为碎屑流。滑坡-碎屑流总的水平运动距离为4022m,垂直运动距离为1443m,遗留的堆积物体积达5×107m3。滑坡距映秀—北川断裂仅36km,位于其下盘,地震烈度达XI度。滑坡导致文家沟中48人遇害,并形成一条完整的地震次生地质灾害链。初步分析表明滑坡启动速度快,滑坡向碎屑流转化过程明显、地点明确。碎屑流运动过程复杂,伴有强烈的“气垫效应”和“前缘气浪冲击效应”。作者认为,文家沟滑坡的高启动速度是长持时强烈地震动作用的结果,与山体的猛烈碰撞是导致滑体解体并转化为碎屑流的原因。     

Preface: Orbits around asteroid

One of the core issues in modern celestial mechanics is the orbital dynamics in the near-regime gravitational field of as- teroids, which provides deep insights into the mathematical nature of a class of nonlinear systems, and plays as a critical basis for in situ explorations of different science goals. Lots of efforts have been made to reveal the characteristics of orbital motion in the vicinity of asteroids, and to improve the skills of asteroid research in methodology.     

Effet de l'anisotropie élastique cristalline sur la distribution des facteurs de Schmid à la surface des polycristauxInfluence of crystalline elasticity anisotropy on Schmid factor distribution at the free surface of polycrystals

Experimental studies of the plasticity mechanisms of polycrystals are usually based on the Schmid factor distribution supposing crystalline elasticity isotropy. A numerical evaluation of the effect of crystalline elasticity anisotropy on the apparent Schmid factor distribution at the free surface of polycrystals is presented. Cubic elasticity is considered. Order II stresses (averaged on all grains with the same crystallographic orientation) as well as variations between averages computed on grains with the same crystallographic orientation but with different neighbour grains are computed. The Finite Element Method is used. Commonly studied metals presenting an increasing anisotropy degree are considered (aluminium, nickel, austenite, copper). Concerning order II stresses in strongly anisotropic metals, the apparent Schmid factor distribution is drifted towards small Schmid factor values (the maximum Schmid factor is equal to 0.43 instead of 0.5) and the slip activation order between characteristic orientations of the crystallographic standard triangle is modified. The computed square deviations of the stresses averaged on grains with the same crystallographic orientation but with different neighbour grains are a bit higher than the second order ones (inter-orientation scatter). Our numerical evaluations agree quantitatively with several observations and measures of the literature concerning stress and strain distribution in copper and austenite polycrystals submitted to low amplitude loadings. Hopefully, the given apparent Schmid factor distributions could help to better understand the observations of the plasticity mechanisms taking place at the free surface of polycrystals. To cite this article: M. Sauzay, C. R. Mecanique 334 (2006).     

Nonlinear behavior of matrix-inclusion composites under high confining pressure: application to concrete and mortar

This paper is devoted to a micromechanics-based simulation of the response of concrete to hydrostatic and oedometric compressions. Concrete is described as a composite made up of a cement matrix in which rigid inclusions are embedded. The focus is put on the role of the interface between matrix and inclusion which represent the interfacial transition zone (ITZ). A plastic behavior is considered for both the matrix and the interfaces. The effective response of the composite is derived from the modified secant method adapted to the situation of imperfect interfaces. To cite this article: T.H. Le et al., C. R. Mecanique 336 (2008).     

Dynamics and Breakup of Viscoelastic Liquids (A Review)

The phenomena of hydrodynamic breakup of liquid jets, drops, films, bridges, and filaments are reviewed for liquids with viscoelastic properties. The reasons for breakup are capillary instabilities, collisions with rigid obstacles, and other forms of dynamic action. The relationship between the properties of the liquids and the features of the breakup process is discussed.     

Flow stability analysis and excitation using pulsating jets

Classical flow stability applied to transition from laminar to turbulent flow may also describe the behavior of vorticity fluctuations created by a pulsating jet placed along a solid boundary. A numerical laminar flow experiment involving a pulsating jet placed along the surface of a duct with flow separation downstream, resulted in eliminating most part of the separated flow region. Applying the same approach to a turbulent flow, it was possible to develop a turbulent stability flow formulation and apply successfully turbulent pulsating jet flow separation control. To cite this article: D. Skamnakis, K. Papailiou, C. R. Mecanique 333 (2005).