Eigenfrequency analysis of cable structures with inclined cables

The approximate eigenfrequencies for the in-plane vibrations of a cable struc- ture consisting of inclined cables,together with point masses at various points were com- puted.It was discovered that the classical transfer matrix method was inadequate for this task,and hence the larger exterior matrices were used to determine the eigenfrequency equation.Then predictions of the dynamics of the general cable structure based on the asymptotic estimates of the exterior matrices were made.     

Hyperbolicity of Velocity-Stress Equations for Waves in Anisotropic Elastic Solids

This paper reports mathematical properties of the three-dimensional, first-order, velocity-stress equations for propagating waves in anisotropic, linear elastic solids. The velocity-stress equations are useful for numerical solution. The original equations include the equation of motion and the elasticity relation differentiated by time. The result is a set of nine, first-order partial differential equations (PDEs) of which the velocity and stress components are the unknowns. Cast into a vector-matrix form, the equations can be characterized by three Jacobian matrices. Hyperbolicity of the equations is formally proved by analyzing (i) the spectrum of a linear combination of the three Jacobian matrices, and (ii) the eigenvector matrix for diagonalizing the linearly combined Jacobian matrices. In the three-dimensional space, linearly combined Jacobian matrices are shown to be connected to the classic Christoffel matrix, leading to a simpler derivation for the eigenvalues and eigenvectors. The results in the present paper provide critical information for applying modern numerical methods, originally developed for solving conservation laws, to elastodynamics.     

Bifurcations and chaos of an inclined cable

The nonlinear behavior of an inclined cable subjected to a harmonic excitation is investigated in this paper. The Galerkin’s method is applied to the partial differential governing equations to obtain a two-degree-of-freedom nonlinear system subjected to harmonic excitation. The nonlinear systems in the presence of both external and 1:1 internal resonances are transformed to the averaged equations by using the method of averaging. The averaged equations are numerically examined to obtain the steady-state responses and chaotic solutions. Five cascades of period-doubling bifurcations leading to chaotic solutions, 3-periodic solutions leading to chaotic solution, boundary crisis phenomena, as well as the Shilnikov mechanism for chaos, are observed. In order to study the global dynamics of an inclined cable, after determining the averaged equations of motion in a suitable form, a new global perturbation technique developed by Kova?i? and Wiggins is used. This technique provides analytical results for the critical parameter values at which the dynamical system, through the Shilnikov type homoclinic orbits, possesses a Smale horseshoe type of chaos.     

层状弹性半空间非轴对称动力问题的奇异解

在柱坐标系下，利用关于方位角的Ｆｏｕｒｉｅｒ变换及关于径向的Ｈａｎｋｅｌ变换，将弹性力学基本方程组转化为非齐次的一阶常微分方程组的标准形式．采用求解微分方程组的矩阵法，建立了介质层的传递矩阵．由层间完全接触条件，导出了在任意埋藏源作用下层状弹性半空间频域奇异解，时域奇异解可通过关于频率的Ｆｏｕｒｉｅｒ积分得到．该方法可应用到固体、流体层的情况     

中心刚体-楔形梁-质点刚柔耦合系统动力学分析

研究了中心刚体-楔形梁-质点系统的固有特性和动力学响应.楔形梁为Euler-Bernoulli梁,高度和宽度都沿着梁的长度方向线性变化.利用广义Hamilton原理和一阶近似耦合模型得到了含有楔形梁完全耦合且时变的微分/代数控制方程.考虑了离心刚化效应,利用有限元得到了系统完全耦合的有限维方程.忽略轴向与横向位移的相互作用,得到了系统的一致质量、阻尼和刚度矩阵.最后对楔形梁和等截面梁在有无端部质点的四种结构进行仿真,结果表明存在显著差异,重点比较了同等条件下楔形梁与等截面梁的差异指数,说明均匀梁和楔形梁的截面细微的差别能够导致系统频率和动力学响应的明显差别.指出实际系统中使用楔形梁模型能够得到更为精确的仿真结果.     

三维声学特征频率灵敏度分析的边界元法

声系统特征频率的灵敏度分析为其优化设计提供了基础,具有重要意义。边界元法在声学问题的求解中具有独特优势,但因其系统方程系数矩阵的频率相关性导致的非线性特征值问题给声学特征频率的灵敏度分析带来了很大困难。为此,本文首先对非线性特征值问题进行了线性化处理,利用围道积分投影方法将非线性特征方程转换为小规模广义特征方程,然后对其关于设计变量直接求导,并引入左特征向量和转换矩阵构造了一种适用于内外声场的三维声学单/重特征频率灵敏度分析的边界元法。数值算例验证了该方法的适用性,以及对单/重特征频率灵敏度的计算精度。     

Nonlinear dynamic response and active vibration control of the viscoelastic cable with small sag

IntroductionCablesareveryefficientstructuralmembersandhencehavebeenwidelyusedinmanylong_spanstructures,includingcable_supportbridges,guyedtowersandcable_supportroofs.Sincecablesarelight,veryflexibleandlightlydamped ,structuresutilizingcables,i.e .,cable_structuresystems,usuallyhavevariousdynamicproblems.Theirmodelsarethereforeverimportantinpredictingandcontrollingtheirresponses.Inthelastdecade,thenonlineardynamicvibrationandstabilitybehaviorofcablesandcable_structureshavedrawntheattentionofman…     

非半单分叉的Normal Form计算

在文［１］基础上，提出一种仅知道派生线性系统零实部特征值时求解非线性系统非半单分叉ＮｏｒｍａｌＦｏｒｍ的方法．通过适当的分类，将要求解的线性代数方程组分为若干相互独立的方程组．将所求系数向量按字典序列排列后，各独立方程组的系数矩阵是上三角矩阵．在非共振情形，各系数向量可按反字典序列递推求出．在共振情形，根据文中的二个定理，巧妙地由一简单的常数矩阵的最大秩子矩阵，定位其系数矩阵的满秩子矩阵，解决了这类方程组的降维简化．通过消元法，把简化后的方程化成类似于半单分叉ＮｏｒｍａｌＦｏｒｍ求解过程中方程的形式，其解法也类似．该方法非常易于在计算机代数软件平台上程序化．     

Forces associated with non-linear non-holonomic constraint equations

A concise method has been formulated for identifying a set of forces needed to constrain the behavior of a mechanical system, modeled as a set of particles and rigid bodies, when it is subject to motion constraints described by non-holonomic equations that are inherently non-linear in velocity. An expression in vector form is obtained for each force; a direction is determined, together with the point of application. This result is a consequence of expressing constraint equations in terms of dot products of vectors rather than in the usual way, which is entirely in terms of scalars and matrices. The constraint forces in vector form are used together with two new analytical approaches for deriving equations governing motion of a system subject to such constraints. If constraint forces are of interest they can be brought into evidence in explicit dynamical equations by employing the well-known non-holonomic partial velocities associated with Kane's method; if they are not of interest, equations can be formed instead with the aid of vectors introduced here as non-holonomic partial accelerations. When the analyst requires only the latter, smaller set of equations, they can be formed directly; it is not necessary to expend the labor first to form the former, larger set and subsequently perform matrix multiplications.     

考虑弹簧阻尼作动器解析雅可比矩阵的多刚体动力学分析

弹簧-阻尼-作动器(spring-damper-actuator,SDA)是多体系统中常见的力元,在工程领域中有着广泛的应用.采用绝对坐标方法建立的多体系统动力学控制方程通常是复杂的非线性微分-代数方程组.为了保证数值解的精度和稳定性,通常需要采用隐式算法求解动力学方程,而雅可比矩阵的计算在隐式数值求解过程中至关重要.对于含有SDA的多体系统,SDA造成的附加雅可比矩阵是与广义坐标和广义速度相关的高度非线性函数.目前的很多研究工作专注于广义力向量的计算,然而对附加雅克比矩阵的计算则少有关注.针对含SDA的多刚体系统进行动力学分析,首先基于Newmark算法研究其在动力学方程求解中的雅可比矩阵的构成形式;然后推导SDA的广义力向量对应的附加雅可比矩阵,其中包括广义力向量对广义坐标和对广义速度的偏导数矩阵.最后通过两个数值算例研究附加雅可比矩阵对动力学分析收敛性的影响;数值分析表明:当SDA的刚度、阻尼和作动力数值较大时,SDA导致的附加雅可比矩阵对数值解的收敛性有重要影响;当考虑SDA对应的附加雅可比矩阵时,动力学分析可以以较少的迭代步实现收敛,从而减少分析时间.     

Dynamic instability of inclined cables under combined wind flow and support motion

In this paper an inclined nearly taut stay, belonging to a cable-stayed bridge, is considered. It is subject to a prescribed motion at one end, caused by traveling vehicles, and embedded in a wind flow blowing simultaneously with rain. The cable is modeled as a non-planar, nonlinear, one-dimensional continuum, possessing torsional and flexural stiffness. The lower end of the cable is assumed to undergo a vertical sinusoidal motion of given amplitude and frequency. The wind flow is assumed uniform in space and constant in time, acting on the cable along which flows a rain rivulet. The imposed motion is responsible for both external and parametric excitations, while the wind flow produces aeroelastic instability. The relevant equations of motion are discretized via the Galerkin method, by taking one in-plane and one out-of-plane symmetric modes as trial functions. The two resulting second-order, non-homogeneous, time-periodic, ordinary differential equations are coupled and contain quadratic and cubic nonlinearities, both in the displacements and velocities. They are tackled by the Multiple Scale perturbation method, which leads to first-order amplitude-phase modulation equations, governing the slow dynamics of the cable. The wind speed, the amplitude of the support motion and the internal and external frequency detunings are set as control parameters. Numerical path-following techniques provide bifurcation diagrams as functions of the control parameters, able to highlight the interactions between in-plane and out-of-plane motions, as well as the effects of the simultaneous presence of the three sources of excitation.     

基于Hamilton体系的弹性力学问题的比例边界有限元方法

比例边界有限元方法是求解偏微分方程的一种半解析半数值解法。对于弹性力学问题,可采用基于力学相似性、基于比例坐标相似变换的加权余量法和虚功原理得到以位移为未知量的系统控制方程,属于Lagrange体系。但在求解时,又引入了表面力为未知量,控制方程属于Hamilton体系。因而,本文提出在比例边界有限元离散方法的基础上,利...     

病态代数方程的精细积分解法

基于精细积分思想,提出了一种有效的病态代数方程组求解方法。类似于稳态热传导方程可视为瞬态热传导方程的极限形式,将具有正定对称实系数矩阵的病态代数方程组归结为一个常微分方程组初值问题的极限形式,并在此基础上建立了病态代数方程组的精细积分解法。该方法不仅精度高,而且能以指数速度收敛,具有较高的效率。本文还讨论了病态代数方程...     

?????????????????????????

 H.B.Jayaraman在20世纪80年代推导的悬链线索元有限元 法计算精度高，特别适用于精度要求比较高的大型索结构. 但是， 当索原长Lu的取值与悬索两节点之间的直线长度相近时，迭 代不易收敛，甚至发散. 提出了当该迭代不收敛时，应采用的迭 代策略. 计算结果表明，该方法准确，计算精度高，可供悬索结 构设计、施工时参考.     

奇异Hamilton系统矩阵的精细积分法

常规位移有限元的结构振动方程是n个二阶常微分方程组.采用一般交分原理推导,将结构振动问题引入Hamiltoil体系,将得到2n个一阶常微分方程组.精细积分法宜于处理一阶方程,应用于线性定常结构动力问题求解,可以得到在数值上逼近精确解的结果.对于非齐次动力方程,当结构具有刚体位移时,系统矩阵将出现奇异.本文借鉴全元选大元高斯-约当法求解线性方程组的经验,提出全元选大元法求奇异矩阵零本征解的方法,该方法可以简便快速地寻求奇异矩阵零本征值对应的子空间.利用Hamiltoil体系已有研究成果及Hamilton系统的共轭辛正交归一关系,迅速将零本征值对应的子空间分离出来,通过投影排除奇异部分,然后用精细积分法求得问题的解.数值算例表明,该方法对Hamilton系统奇异问题,处理方便,计算量小,易于实现,同时保持了精细算法的优点.     

Response of Uncertain Nonlinear Vibration Systems with 2D Matrix Functions

This paper extends the response of uncertain nonlinear vibration systems to vector-valued and matrix-valued functions. Random variables and system derivatives are conveniently arranged into 2D matrices. The method is based on a second order expansion of the governing equations and matrix calculus, Kronecker algebra are used in the mathematical development. The results derived are easily amenable to computational procedures.     

Nonlinear Parameter Identification of a Mechanical Interface Based on Primary Wave Scattering

We study stress-wave propagation in an impulsively forced split Hopkinson bar system incorporating a threaded interface. We first consider only primary transmission and reflection and reduce the problem to a first-order, strongly nonlinear ordinary differential equation governing the displacement across the interface, called the primary-pulse model. The interface is modeled as an adjusted-Iwan element, which is characterized by matching experimental and numerical eigenfrequencies as well as primary pulse amplitudes. We find that the adjusted-Iwan element parameters are dependent on preload and impact velocity (input force). A high-order finite element model paired with the identified adjusted-Iwan element is used to simulate multiple transmissions and reflections across the interface. We find that the finite element simulation reproduces the experimental results in both the wavelet and Fourier domains, validating the identification method. Our findings demonstrate that the primary-pulse model can be used for experimental parameter identification of nonlinear interfaces in waveguides.     

多自由度参激系统稳定性的数值分析

提出多自由度周期参激系统稳定性的数值直接法。通过将扰动方程表示成状态方程形式,再根据Flo-quet理论将扰动解表示成指数特征分量与周期分量之积,并将其周期分量与系统周期系数展成Fourier级数,导出一系列代数方程,建立矩阵特征值问题,从而由数值求解特征值可直接确定参激系统的稳定性。该方法可用于一般周期参激阻尼系统,特征值矩阵不含逆子阵。应用于斜拉索在支座周期运动激励下的参激振动不稳定性分析,数值结果表明该方法的有效性。     

Dynamics Investigation of Three Coupled Rods with a Horizontal Barrier

This report is a part of the larger project of non-linear dynamics investigation of three coupled physical pendulums with damping and with arbitrary situated barriers, and externally driven. The set of differential equations and the set of algebraic inequalities (representing a barrier) governing the motion of three coupled rods are presented in the non-dimensional form. The system of governing equations is integrated between two successive impacts, and the discontinuity points are detected (by halving time step until a required precision is obtained). In each impact time, the state of the system is transformed using the extended restitution coefficient rule. The theory of Aizerman and Gantmakher is used to calculate the fundamental solution matrices in the analyzed system exhibiting discontinuities. The fundamental matrices are used during calculation of Lyapunov exponents, during stability analysis of periodic solutions (Floquet multipliers) and in shooting method applied to detect and trace periodic orbits. Some examples for three coupled identical rods with horizontal barrier are reported.     

A composite multigrid method for calculating unsteady incompressible flows in geometrically complex domains

A time-accurate, finite volume method for solving the three-dimensional, incompressible Navier-Stokes equations on a composite grid with arbitrary subgrid overlapping is presented. The governing equations are written in a non-orthogonal curvilinear co-ordinate system and are discretized on a non-staggered grid. A semi-implicit, fractional step method with approximate factorization is employed for time advancement. Multigrid combined with intergrid iteration is used to solve the pressure Poisson equation. Inter-grid communication is facilitated by an iterative boundary velocity scheme which ensures that the governing equations are well-posed on each subdomain. Mass conservation on each subdomain is preserved by using a mass imbalance correction scheme which is secondorder-accurate. Three test cases are used to demonstrate the method's consistency, accuracy and efficiency.