On the high-frequency response of unsteady lift and circulation: A dynamical systems perspective

In this paper, we consider the unsteady aerodynamics of a two-dimensional airfoil as a dynamical system whose input is the angle of attack (or airfoil motion) and output is the lift force. Based on this view, we discuss the evolution of lift and circulation from a purely dynamical perspective through step response, frequency response, transfer function, etc. In particular, we point to the relation between the high-frequency gain of the transfer function and the physics of the development of lift and circulation. Based on this view, we show that the circulatory lift dynamics is different from the circulation dynamics. That is, we show that the circulatory lift is not lift due to circulation. In fact, we show that the circulatory–non-circulatory classification is arbitrary. By comparing the steady and unsteady thin airfoil theory, we show that the circulatory lift possesses some acceleration (added-mass) effects. Finally, we perform simulations of Navier–Stokes equations to show that a non-circulatory maneuver in the absence of a free stream induces viscous circulation over the airfoil.     

运动副摩擦参数的识别方法研究

采用粘性 库仑摩擦模型建立了单自由度回转机械系统运动微分方程,分别从时域和频域角度提出了两种识别运动副摩擦参数的方法。仿真结果表明,两种方法在较高信噪比下均能有效地辨识出库仑摩擦和粘性摩擦参数,其中频域法比时域法具有更高的参数估计精度。对电机驱动的单自由度转动机械系统实验装置的运动副建立了相应的摩擦模型,实验结果表明:在不同试验条件下两种方法计算得到的库仑和粘性摩擦参数具有较好重复性,从而验证了本文参数辨识方法的正确性。     

变截面梁动力响应的频域传递矩阵计算方法

为提高变截面梁振动分析的计算效率,提出了基于频域传递矩阵法的动力计算算法．首先,选择线速度、角速度、弯矩和剪力作为求解变量,通过Laplace变换将变截面梁的动力响应时域偏微分方程转换为频域常微分方程;然后,通过求解频域方程并结合协调和边界条件建立变截面梁的频域传递矩阵;通过构造傅里叶级数展开形式的时域响应函数,对变截面梁传递矩阵方法求解的频响函数进行Laplace逆变换,建立了变截面梁的固有特性计算和时域瞬态响应计算方法,最后,借助数值仿真软件,开发了变截面梁动力响应分析的计算程序．完成对算例的仿真计算和分析,并与有限元计算结果进行对比,数值结果验证了该方法的正确性和有效性．     

On the analyses of vibration and stability for orthotropic circular cylindrical shells conveying fluid

Based on flügge's thin orthotropic shell equations and a linear potential flow theory, buckling instability and flutter-type instability of orthotropic circular cylindrical shells conveying fluid are studied. By means of the expansion in beam-mode functions and Galerkin's method they can be reduced to solving a generalized complex eigenvalue problem. In calculating the generalized forces defined as Eq. (23), a direct numerical integration technique, which has proven to be more efficient than others, is used. For clamped-clamped orthotropic shells several numerical examples show that the effects of axial elastic modulus on the dynamic behaviors are quite different from those of circumferential elastic modulus.     

Time-domain Impedance Boundary Conditions for Computational Acoustics and Aeroacoustics

This paper reviews the short history, motivation, numerical and theoretical issues, and development of methods for treating a boundary as a reflective/absorptive surface for the time-domain computation of waves in general and acoustic waves in particular. It begins with the extension and implementation of the frequency-domain impedance to a time-domain impedance-equivalent boundary condition (TDIBC), and illustrates how the theoretical, numerical, and implementation issues are addressed and resolved for acoustic/aeroacoustic applications. Comments are also made on the extendibility and applicability of the concept and TDIBC to other fields and types of problems.     

Determination of relaxation and retardation spectrum by inverse functional filtering

The article is devoted for the determination of the relaxation and retardation spectrum (RRS) from monotonic time- and frequency-domain material functions by the inverse functional filters executing discrete convolution algorithms for geometrically spaced data. It is shown that the problem of RRS determination from a wide variety of material functions leads to the three inverse filtering tasks on a logarithmic time or frequency scale with the three specific frequency responses concerning: (i) the time-domain functions, (ii) the real parts and (iii) the imaginary parts of the frequency-domain functions, and three algorithms (having the versions with even and odd number of coefficients) are to be applied to: (i) time-domain compliance and modulus functions, (ii) their derivatives, and (iii) frequency-domain functions. It is demonstrated that ill-posedness of an inverse filter manifests as large sampling-rate-dependent noise amplification coefficients. A novel regularization strategy allowing to ensure the desired noise immunity is proposed based on choosing sampling rate for geometrically spaced data. The performance of the algorithms is investigated. Optimal sampling rates are disclosed for specific material functions. The frequency range of 2–3 decades is established to be optimal for the recovery of a single RRS point estimate ensuring maximum accuracy with reasonable noise immunity. Practical algorithms are proposed for recovering RRS from the real and imaginary parts of frequency-domain functions. Some known non-parametric methods are compared with the suggested functional filters.     

平稳高斯激励下线性结构随机振动分析的辅助简谐激励广义法

相比于时域法,频域法是更为高效、易行的随机振动分析方法,但对于平稳激励下的随机振动分析,现有频域方法常需振型截断或功率谱矩阵分解,将会影响计算精度和效率.为此,本文在频域法的框架下,针对平稳高斯激励下线性结构的随机振动分析提出了一种精确且高效的辅助简谐激励广义法.首先,引入广义脉冲响应函数和广义频响函数的概念,推导了与...     

随机游走对罗经方位对准的精度影响分析

分析了随机游走对罗经方位对准的影响,通过罗经方位对准的频域误差传递函数,求出时域的单位脉冲响应函数,进而推导出该系统在白噪声驱动下的时域表达式,最后求出该表达式的统计特性,其标准差即为随机游走对罗经方位对准的误差公式。通过分析误差公式可知,罗经方位对准的精度与随机游走系数成正比,还与设置的半个二阶阻尼振荡周期的平方根成反比。当陀螺随机游走一定时,单纯的通过延长对准时间并不能提高对准精度,还需要延长二阶阻尼振荡周期,数值仿真验证了所提观点的正确性。     

基于高次广义柔度灵敏度的结构损伤识别

柔度矩阵可以由结构的低价模态近似计算获得,因此被广泛用于结构的模型修正和损伤识别中.由普通柔度派生而来的广义柔度,可以由低价模态数据更加精确的获得,且随着广义柔度次数的增高其精度越高,往往只需要第一或二阶模态数据即可获得很准确的高次广义柔度.因此,广义柔度灵敏度方法自提出以来受到广泛关注.论文详细研究了基于高次广义柔度灵敏度的损伤识别计算方法,研究中发现,利用广义柔度灵敏度进行损伤识别计算时,并非越高次的广义柔度其识别结果越准确,随着广义柔度次数的增加,损伤识别结果精度呈现出先提高但随后显著降低的趋势.究其原因在于,虽然随着广义柔度次数的增加,广义柔度本身的精度更高,但与之相应的灵敏度方程组系数矩阵的条件数却也显著增大了,即方程组的病态性反而更加严重了,这导致了基于高次广义柔度计算所得的损伤参数的精度反而不如低次广义柔度的情况.因此,论文的研究表明,工程中利用广义柔度进行模型修正或损伤识别时,一般采用一次广义柔度或二次广义柔度即可,且计算中为了克服方程组的病态性和数据噪声的不利影响,论文提出了一种反馈奇异值截断法,能够明显提高计算精度,获得较准确的识别结果.     

PATH INTEGRAL SOLUTION OF NONLINEAR DYNAMIC BEHAVIOR OF STRUCTURE UNDER WIND EXCITATION

IntroductionThe dynamic behavior of the nonlinear structure under wind excitation has beenobserved very complicated.Taking guyed masts as an example,only a few collapsingaccidents occurred under extreme atmospheric conditions[1],many took place under mild…     

STUDY OF PARAMETRIC VIBRATION OF GYROSCOPIC CONTINUA BASED ON NONLINEAR NORMAL MODES AND A NUMERICAL ITERATIVE APPROACH

以脉动流输流管为例,利用非线性模态技术和一种数值迭代法研究陀螺连续体的非线性参数振动响应问题. 通过谐波平衡法将系统非线性非自治控制方程转化为拟自治方程,并在状态空间上利用不变流形法构造系统的非线性模态. 以对应自治系统的解为初值,采用一种数值迭代法来求解拟自治控制方程的模态系数,结果证明了该迭代法的快速收敛性. 在频域分析中得到了幅频响应和相空间上的不变流形,而在时域复模态分析中则发现了参激陀螺系统的正交相位差和行波振动现象.     

基于非线性模态和一种数值迭代法的陀螺连续体参数振动研究

以脉动流输流管为例,利用非线性模态技术和一种数值迭代法研究陀螺连续体的非线性参数振动响应问题. 通过谐波平衡法将系统非线性非自治控制方程转化为拟自治方程,并在状态空间上利用不变流形法构造系统的非线性模态. 以对应自治系统的解为初值,采用一种数值迭代法来求解拟自治控制方程的模态系数,结果证明了该迭代法的快速收敛性. 在频域分析中得到了幅频响应和相空间上的不变流形,而在时域复模态分析中则发现了参激陀螺系统的正交相位差和行波振动现象.     

基于有限条带的厚/薄板矩形通用单元

基于两广义位移梁理论，利用解析试函数来建立厚／薄梁单元的横向位移、转角位移、曲率、剪应变等位移模式，构造出厚／薄梁通用单元．应用有限条带，将厚／薄梁单元的位移模式应用于厚／薄板矩形弯曲单元，直接构造出单元的横向位移、转角、曲率、剪应变，导出了单元的刚度矩阵和结点等效力，编制了计算程序，进行了数值计算和比较，结果表明，所研究的单元不出现剪切闭锁且精度较好．     

A direct finite element method study of generalized thermoelastic problems

By using the principle of virtual work, the finite element equations corresponding to the generalized thermoelasticity with two relaxation times (i.e., the G–L theory) are derived. In order to improve the accuracy of integral transform method, especially in two/three-dimension, the equations are solved directly in time-domain. As numerical examples, the developed method is used to investigate the generalized thermoelastic behavior of a slim strip and a semi-infinite plate subject to thermal shock loading. The results demonstrate that the developed method can faithfully predict the deformation of the structure and most importantly the delicate second sound effect of heat conduction in both one and two-dimensional solids whilst it is usually difficult to model, especially in two-dimensional, by using the transform method. This method can also be used to study generalized piezothermoelastic and magnetothermoelastic problems.     

成层饱和介质平面波斜入射问题的一维化时域方法

地震波斜入射下自由场的输入是大型结构抗震分析中亟待解决的问题之一,尤其是成层饱和多孔介质自由场问题,由于问题的复杂性,目前研究甚少. 本文基于Biot提出的饱和多孔介质动力方程,建立了一种新的求解平面波斜入射下基岩上覆饱和多孔介质成层场地自由场分析的一维化时域计算方法. 该方法首先根据Snell定律将饱和多孔介质二维空间问题转化为一维时域问题,通过对深度方向的有限元离散,得到饱和多孔介质波动问题的一维化有限元方程,然后采用单相弹性介质精确人工边界条件模拟基岩半空间的波动辐射和输入特征,通过考虑基岩与饱和多孔介质间透水或不透水边界条件以及不同饱和多孔介质交界面边界条件,形成基岩上覆成层饱和介质系统的整体有限元方程,最后采用中心差分法与Newmark平均加速度近似格式相结合的方法对时间进行离散,得到节点的动力时程的显式表达. 典型场地的地震反应分析表明,本文方法的计算结果与传递矩阵法结合傅里叶变换的计算结果完全吻合,证明了其有效性.