Finite strip buckling analysis of curved plate assemblies under biaxial loading

A finite strip method is presented for calculating the linear buckling stresses of structural assemblies of long, thin plate components which, in general, are curved and which are rigidly joined together at their longitudinal edges. It is assumed that on buckling under the action of a biaxial direct stress field the perturbation forces and displacements all vary sinusoidally in the longitudinal direction. A stiffness matrix relating the amplitudes of the perturbation forces and displacements is developed for the curved strip on the further assumption of relatively high-order polynomial variations of the displacement components around the plate width. Numerical results are presented of the application of the curved strip in calculating the buckling stresses of plates, cylinders, panels and formed sections.     

Investigation of the Nonlinear Behavior of Beams*

ABSTRACT A system of equations that describes the nonlinear behavior of teams is presented. The differential equations are solved by using the Galerkin method, and a system of nonlinear algebric equations is obtained. A method to deal with discontinuities in structural properties and load distribution along a beam is presented. The derivation includes general expressions for gravity-type loads. Two methods of calculating the resultant moment along the beam are described. One method incorporates differentiation of the expressions for displacements, while the second method is based on integration of loads along the beam. The numerical model is used in order to investigate displacements and moment resultants along a cantilevered beam. Different problems that are associated with nonlinear behavior are presented and discussed. Good agreement is obtained between theoretical and experimental results.     

特征值摄动法在利用动响应结构小损伤检测中的应用

目前在使用遗传算法或神经网络方法进行结构动力学损伤检测,需要基于少量的在线测量损伤结构数据和大量的数值仿真数据来实现,其中通过有限元方法来获得仿真数据的巨大计算量是动力学结构损伤检测方法发展中所面临的一个重要问题。本文在建模方面应用近年来提出的调整单元刚度模拟损伤的先进方法,以保证在损伤前后结构自由度数目不变;在此基础上应用特征值摄动法来减少损伤检测中计算量,并通过对复合材料层合板响应信号的小波分析验证了使用一阶矩阵摄动在有效降低计算量的同时,可以获得对损伤检测而言足够准确的响应信号。     

动态结构系统可靠性的频率灵敏度分析

基于频率可靠性理论提出了动态结构系统的频率可靠性灵敏度的计算方法. 根据动态结构系统的固有频率与激振频率差的绝对值不超过规定值的关系准则,定义了随机结构振动问题的可靠性模式和系统的可靠度,进而提出了频率可靠性灵敏度分析方法,应用随机摄动技术、可靠性理论和灵敏度方法,推导出参数为正态随机变量时的准失效概率灵敏度的表达公式,对动态结构系统共振问题的可靠性灵敏度分析方法进行了探索. 并在此基础上,以一简化的车身-车轮-路面关系的系统模型为例,由可靠性灵敏度矩阵可以看出,各随机参数的变化对动态结构系统的可靠性及失效性的影响程度也不同,其计算结果与通常的定性分析结果完全一致,验证了所提方法的有效性. 该方法为动态结构系统的优化和设计提供了依据.      

An analysis of perturbation based methods for the treatment of parameter uncertainty in numerical groundwater models

Monte Carlo methods are robust approaches for the estimation of prediction uncertainty in groundwater flow and transport modelling under uncertain model parameters. However Monte Carlo procedures estimate the prediction statistics by generating a population of solutions from random realisations of the model parameters which are consistent with the parameter statistics, and as a result are computationally demanding. Taylor series based procedures offer an alternative to Monte Carlo methods for calculating the prediction statistics. Two such approaches, the first-order second moment and McLaughlin and Wood's perturbation method, are based on using a Taylor series to derive approximate expressions for the model predictions first and second statistical moments. In this paper the perturbation method presented by McLaughlin and Wood is rederived using Vetter matrix notation. This is compared with the first-order second moment (FOSM) method and while the steady state expressions for these two approaches are shown to be equivalent, the transient forms are considerably different. A new form of the FOSM is derived, which is simpler and has a lower computational burden. However, the transient McLaughlin and Wood expression is found to have a significantly lower computational overhead than either of the FOSM methods presented.     

结构优化半解析灵敏度分析的改进算法

提出了结构半解析灵敏度分析的改进算法,该算法实现简便,对于设计变量摄动步长具有极佳的数值稳定特性。首先,从总体角度推导静力问题半解析法灵敏度分析新算法,提出了位移与应力灵敏度列式,并给出了算法实施途径;然后,将此思路推广于自振频率、屈曲临界荷载和瞬态响应等多种问题,提出了相应的计算步骤。以梁单元与壳单元等典型结构为例,开展了多个算例测试。测试表明,改进算法计算精度和效率均有提升,特别是设计变量步长有更大的数值稳定区域,为复杂工程结构形状优化的灵敏度分析提供了新途径。     

多自由度系统复模态理论的摄动方法 (三)非对称系统复特征值问题的摄动迭代法

基于实模态理论和相应特征值计算方法,本文在位形空间中用摄动方法求解非对称系统的复特征值和复模态。为提高精度,增加了一个迭代格式。最后给出算例并讨论了计算精度和收敛性。     

多柔体系统响应计算的子循环计算方法研究

过去近30年中,柔性多体系统动力学研究取得了巨大的进展,人们的兴趣集中在柔性多体系统建模、计算及实验研究等3个方面. Belytschko等于1979年提出的子循环算法已经成功地应用于结构动力响应的有限元计算中,然而有关柔性多体动力学的子循环算法研究尚未见报道. 该文提出了一种适合于柔性多体系统响应计算的中心差分类子循环算法,在将非线性微分-代数混合方程组(DAEs)缩并为纯微分方程组(ODE)的基础上,推导出快、慢变分量的同步更新公式和子步更新公式;在变量的数值积分过程中,采用能量平衡计算检查算法的稳定性;算例结果表明该算法可以在保持合适的精度要求下,有效地提高响应的计算效率;对积分步长进行摄动修正可以保持算法的稳定性.      

基于不完备频响函数的结构损伤统计识别研究

基于不完备频响函数数据,结合概率统计方法,提出了一种能同时考虑模型参数不确定性和测试噪声影响的结构损伤统计识别方法。首先,基于频响函数某一行向量在不同频率下的幅值数据,利用矩阵拉直运算建立了关于损伤系数的确定性识别方程。其次,假设模型参数误差和测试噪声为零均值的高斯随机变量,根据摄动理论,推导了损伤后结构刚度参数的前二阶矩。随后,利用结构损伤前后的概率分布得到了结构损伤存在概率。最后,通过一个平面桁架结构模型验证了本文方法的有效性。数值算例的研究结果表明,损伤单元的损伤存在概率远大于非损伤单元;测试噪声对识别结果的影响比模型参数不确定性的影响更为显著。     

一种Birkhoff形式下结构动响应问题的保辛中点格式

结构动响应预测是结构设计的基础,是结构振动控制、载荷识别的前提。本文在辛体系下针对结构动响应问题,提出了一种Birkhoff形式下的保辛中点格式。首先引入状态变量,并基于摄动方法将结构动响应方程转化为线性自治Birkhoff方程的形式,进一步利用中心差分推导出线性自治Birkhoff方程的中点格式,其证明是保辛的。该格式不要求Birkhoff方程系数矩阵非奇异,因此适用于奇数维系统。两个不同数值算例的结果充分验证了本文方法的卓越性,也凸显了相对于传统算法在计算精确度和稳定性方面的明显优势。     

Stability for manifolds of equilibrium states of generalized Hamiltonian systems

For a generalized Hamiltonian system, stability for the manifolds of equilibrium states is presented based on Lyapunov’s stability theories. Equilibrium equations, perturbation equations and first approximate equations of the system are given. A theorem for the stability of manifolds of equilibrium states of general autonomous system is used to the generalized Hamiltonian system, and three propositions on the stability of manifolds of equilibrium states of the system are obtained. Two examples are given to illustrate application of the method and results.     

Sensitivity analysis based on Lanczos algorithm in structural dynamics

IntroductionTheproblemforcomputingeigenpairderivativeshasoccupiedmanyresearchersinthepastseveraldecades.SurveypaperbyAdelmanandHaftka[1]isanexcellentintroductiontothetopicandalsoindicatevariousdisciplinesofapplications.Thereasonwhysomanypeopleareinterestedinthisproblemisthatsensitivitiesplayaveryimportantroleinstructuraldynamics,optimalanalysisandcontrolsystemdesign ,andmayalsobeanexpensiveandtime_consumingtask.Researchonsensitivityofmultipleeigenvalueshasbeenanareaofsignificantinterestinrece…     

振动陀螺仪单向时滞耦合系统的摄动分析

针对振动陀螺仪单向时滞耦合系统，研究了角速度、耦合强度、耦合时滞与系统输出信号之间的关系，着重关注耦合时滞对待检测角速度测量结果的影响．本文采用摄动法对振动陀螺仪单向时滞耦合系统进行理论分析，得到了系统输出信号、系统振动边界和待检测角速度的近似表达式，并将数值和近似计算结果进行了对比．理论分析结果与数值结果比较吻合，验证了运用摄动法分析振动陀螺仪单向时滞耦合系统的有效性，也表明了在运用摄动法分析系统时具有局限性．     

随机杆系结构几何非线性分析的递推求解方法

建立了随机静力作用下考虑几何非线性的随机杆系结构的随机非线性平衡方程. 将和 位移耦合的随机割线弹性模量以及随机响应量表示为非正交多项式展开式，运用传统的摄动方法获 得了关于非正交多项式展式的待定系数的确定性的递推方程. 在求解了待定系数后，利用非 正交多项式展开式和正交多项式展开式的关系矩阵，可以很方便地得到未知响应量的二阶统计矩. 两杆结构和平面桁架拱的算例结果表明，当随机量涨落较大时，递推随机有限元方法比基于 二阶泰勒展开的摄动随机有限元方法更逼近蒙特卡洛模拟结果，显示了该方法对几何非线性 随机问题求解的有效性.     

外压加筋圆柱壳总体屈曲Koiter——边界层奇异摄动法

将Ｋｏｉｔｅｒ理论和奇异摄动理论中的边界层法相结合处理加筋圆柱壳无因次化非线性边界层型Ｋａｒｍａｎ－Ｄｏｎｎｅｌ方程由分支点和边界层导致的双重奇异性，提出外压加筋圆柱壳总体屈曲Ｋｏｉｔｅｒ—边界层奇异摄动法。从摄动意义上分析边界条件，前屈曲非线性和初始几何缺陷对外压加筋圆柱壳屈曲载荷的影响。算例表明，本方法具有良好的计算效率和计算精度，与数值解相比更能揭示内在影响规律。     

结构优化半解析灵敏度及误差修正改进算法

提出结构半解析灵敏度分析及其针对刚体位移的误差修正方法的改进算法, 构建灵敏度分析与误差修正项可分离形式. 该方法实现简便, 数值精度不受摄动步长与单元数目的影响. 首先从总体角度推得静力问题的误差修正半解析灵敏度分析方法, 提出了位移误差修正灵敏度列式, 并给出算法实施途径; 然后将此思路推广于自振频率、屈曲临界载荷问题, 提出了相应的计算步骤. 随后, 给出梁单元与壳单元误差修正项的具体推导方法, 并分别使用两种单元构建有限元模型进行算例测试. 结果表明, 该方法适用于多种分析类型, 数值精度不受单元数目与摄动步长的影响. 由于灵敏度分析与误差修正项可以分开计算, 该方法支持将误差修正项直接叠加于灵敏度求解结果进行误差修正, 使已有灵敏度分析程序得到充分利用. 尤其对于复杂工程结构的优化设计, 特别是形状优化设计以及尺寸、形状混合优化设计, 相比于原误差修正方法, 实现更为简便, 效率有所提升, 能为半解析灵敏度分析方法及其程序实现提供新的思路.      

Iterative harmonic balance for period-one rotating solution of parametric pendulum

In this study, an iterative method based on harmonic balance for the period-one rotation of parametrically excited pendulum is proposed. Based on the definition of the period-one rotating orbit, the exact form of the solution can be obtained using the Fourier series. An iterative harmonic balance process is proposed to estimate the coefficients in the exact solution form. The general formula for each iteration step is presented. The method is evaluated using two criteria, which are the system energy error and the global residual error. The performance of the proposed method is compared with the results from multiscale method and perturbation method. The numerical results obtained with the Dormand?CPrince method (ODE45 in MATLAB?) are used as the baseline of the evaluation.     

Simulation of the aerodynamics of buildings and structures by means of the closed vortex loop method

Problems of the numerical simulation of the air flow past buildings and structures are considered using the closed vortex loop method. A mathematical model, based on the vortex approach, of the time-dependent ideal incompressible fluid flow past a system of bodies is proposed. A numerical scheme for solving the problem and an algorithm for calculating the distributed wind loads over the body surface are outlined. An example of calculating the aerodynamic loads is given for a real building and the results are compared with the known results of testing a model of the building in a wind tunnel. An example of the calculation and analysis of the wind distribution over a system of several buildings is also presented.     

THE STUDY ON A KIND OF CONTROL SYSTEM WITH NONLINEAR PARABOLIC DISTRIBUTED PARAMETERS

ＩｎｔｒｏｄｕｃｔｉｏｎＩｎｔｈｅａｃｔｉｖｉｔｉｅｓｓｕｃｈａｓｔｈｅａｅｒｏｓｐａｃｅｅｎｇｉｎｅｅｒｉｎｇｈｅａｔｐｒｏｔｅｃｔｉｏｎ[1,2 ],ｔｈｅｌａｒｇｅ＿ｓｃａｌｅｃｏｎｃｒｅｔｅｐｏｕｒｉｎｇ ,ｔｈｅｃｈｅｍｉｃａｌｐｒｏｄｕｃｔｉｏｎ ,ｔｈｅｐｏｗｄｅｒｍｅｔａｌｌｕｒｇｙａｎｄｔｈｅｆｉｒｅｓａｆｅｔｙｅｎｇｉｎｅｅｒｉｎｇ[3,4],ｔｈｅｒｅｉｓｏｎｅｋｉｎｄｏｆｎｏｎｌｉｎｅａｒｐａｒａｂｏｌｉｃｄｉｓｔｒｉｂｕｔｅｄｐａｒａｍｅｔｅｒｃｏｎｔｒｏｌｓｙｓｔｅｍｗｉｔｈｍｏｖｉ…     

Spectral Response of a Stochastic Oscillator under Impacts

An approximate analytical procedure is presented to estimate theresponse spectrum of an oscillator with elastic impacts under a Gaussian whitenoise excitation. The proposed approach is based on a perturbation analysis ofthe problem and on the use of the stochastic averaging principle. The basicidea is to replace the initial system by a more regular system obtained byapproximating the nonlinear restoring force by a Chebychev polynomial, and thento construct for this regular system two approximations: one for the flowand one for the stationary distribution of the response amplitude. Ananalytical approximation of the response spectrum can then be derived fromthese results. Predictions from this analytical approximation are compared with corresponding digital simulation estimates and with the ones obtained from theconventional equivalent linearization method.