System identification with reaction-mass devices

Reaction-mass devices are to be used in a dual role: first for the purpose of system identification, and second for the purpose of vibration suppression. Unlike ground-based shakers, reaction-mass devices have a position constraint which must be maintained. A two-part system identification procedure is described which maintains this constraint. In the first part, a new algorithm is used to obtain the impulse response from input-output data. The second part uses the eigensystem realization algorithm to determine linear matrix models of the structure and of the reaction-mass device. This procedure is demonstrated on a small experimental structure. The usefulness of the model for the design of a vibration suppression controller is also demonstrated.     

粘性流体中弹性板振动的有限元耦合问题

流体－结构耦合作用问题是工程中比较常见的问题，具有重要意义，由于流体计算的复杂性，迄今为止，大部分的流体－结构耦合分析都是建立在对流体充分简化的基础上，尤其是将流体视为无粘、无旋的理想流体。该文在近年来前人工作的基础上，发展了一种流体－弹性结构耦合计算模式。流体视为不可压、有粘的介质，流场没有自由表面。该文采用ＳＵ／ＰＧ方法形成流体的有限元方程，采用ＡＬＥ格式处理流体和结构之间的移动界面。采用预估     

Time marching for simulation of fluid–structure interaction problems

Numerical simulation of industrial multi-physics problems is still a challenge. It generally requires large computational resources. It may involve complex code coupling techniques. It also relies on appropriate numerical methods making data transfer possible, quick and accurate. In the framework of partitioned procedures, multi-physics computations require the right choice of code coupling schemes, because several physical mechanisms are involved. Numerical simulation of fluid–structure interactions is one of these issues. It is investigated in this paper. First the computational process involving a code coupling procedure is presented. Then, applications and test cases involving fluid structure interactions are investigated using several examples. A partitioned procedure involves several operators ensuring code coupling. A special attention must be paid to energy conservation at the fluid–structure interface, especially when it is moving and when strong non-linear behaviour occurs in both fluid and structure systems. In the present work, several fluid–structure code-coupling schemes are compared and discussed in terms of stability and energy conservation properties. The criteria are based on the evaluation of the energy that is numerically created at the fluid–structure interface. This is achieved by considering the staggering process due to the time lag between the fluid and structure solvers. Comparisons are made, and finally the article gives recommendations for creating a tool devoted to coupled simulations of fluid structure interactions.     

Numerical form-finding of tensegrity structures

A novel and versatile numerical form-finding procedure that requires only a minimal knowledge of the structure is presented. The procedure only needs the type of each member, i.e. either compression or tension, and the connectivity of the nodes to be known. Both equilibrium geometry and force densities are iteratively calculated. A condition of a maximal rank of the force density matrix and minimal member length, were included in the form-finding procedure to guide the search of a state of self-stress with minimal elastic potential energy. It is indeed able to calculate novel configurations, with no assumptions on cable lengths or cable-to-strut ratios. Moreover, the proposed approach compares favourably with all the leading techniques in the field. This is clearly exemplified through a series of examples.     

非线性方程组的解法:局部弧长法

描述了一个新的非线性方程组的求解方法——局部弧长法．该方法是在弧长法的基础上发展起来的适合于材料非线性有限元分析的数值解法．其约束方程充分利用了结构中破坏区域内的非线性变形信息，有效地解决了材料非线性分析中的稳定性与收敛性问题．数值计算表明，该方法不仅适合于求解结构的极限承载能力，也适合于求解结构达到极限承载参力以后的荷载－变形的全过程     

Investigation of shock/elastic obstacles interactions by means of a coupling technique

Arrays of obstacles have proved to be an efficient way of attenuating shock waves generated by large scale explosions. The present study intends to take into account fluid–structure interactions that may occur when elastic obstacles are used. A tractable coupling tool based on a partitioned procedure is exposed, validated on the supersonic flutter of a panel and applied to a configuration composed of square section cylinders. Several numerical difficulties related to staggering are emphasized and workarounds discussed. A methodical procedure involving one and two-way coupled simulations highlights the influence of material properties as well as the acceleration of waves in the fluid when an initial motion is prescribed to the obstacles. Finally, it is shown that uncoupled simulations may be relevant to investigate shock mitigation in some given cases.     

On Force Stabilization of a Dissipative System

The class of ADGFE systems that can be stabilized by gyroscopic and positional nonpotential forces is defined more exactly. A procedure is proposed for selection of these forces, explicitly stabilizing the motion of a dissipative mechanical system in even- and odd-dimensional cases under some additional constraints on the mutual structure of the matrices of dissipative and potential coefficients     

流固偶合运动的边界元法

本文给出了流固偶合运动(包括物体散射辐射及偶合运动)的边界元法理论和应用.对于散射问题,求出了物体引起的散射势及入射波作用于物体的载荷.对于辐射问题,求出了辐射势及物体在流体中运动的附加质量和附加阻尼.偶合问题包括求其中包含的散射势和辐射势以及作用于物体之上的散射力、物体的附加质量、附加阻尼、物体在入射波作用下的运动.在偶合运动问题中,本文采取了边界积分方程与物体在流体中的运动方程联立求解的方法,并将其运用到边界元法的数值过程中.所编制的程序有较高的精度.最后给出了数值计算结果与理论解的比较.     

Combined interface boundary condition method for coupled thermal simulations

A new procedure for modeling the conjugate heat‐transfer process between fluid and structure subdomains is presented. The procedure relies on higher‐order combined interface boundary conditions (CIBC) for improved accuracy and stability. Traditionally, continuity of temperature and heat flux along interfaces is satisfied through algebraic jump conditions in a staggered fashion. More specifically, Dirichlet temperature conditions are usually imposed on the fluid side and Neumann heat‐flux conditions are imposed on the solid side for the stability of conventional sequential staggered procedure. In this type of treatment, the interface introduces additional stability constraints to the coupled thermal simulations. By utilizing the CIBC technique on the Dirichlet boundary conditions, a staggered procedure can be constructed with the same order of accuracy and stability as those of standalone computations. Using the Godunov–Ryabenkii normal‐mode analysis, a range of values of the coupling parameter is found that yields a stable and accurate interface discretization. The effectiveness of the method is investigated by presenting and discussing performance evaluation data using a 1D finite‐difference formulation for each subdomain. Copyright © 2007 John Wiley & Sons, Ltd.     

On the contact of a spherical membrane enclosing a fluid with rigid parallel planes

The mechanical response of an inflated spherical membrane-fluid structure in contact with rigid parallel planes is studied. The membrane is assumed to be a two-dimensional non-linear elastic and isotropic structure, while no assumption is imposed on the fluid. A numerical procedure is employed to compute the equilibrium configurations of the membrane-fluid structure. This study provides information regarding the contact force, stress distribution and pressure in the membrane and in the enclosed fluid, respectively. It was observed that a transition between unwrinkled to partially wrinkled configurations of the membrane occurs subjected to the loading conditions. Further investigation of the wrinkled configurations is presented.     

基于应变脉冲响应协方差的损伤识别方法研究

基于振动参数的结构损伤识别,是近年来土木工程的热点研究课题,振动参数包括频率、振型、频响函数、模态应变能、应变响应和加速度响应等,当结构损伤时,损伤位置附近将产生应力重分布,从而引起应变的变化,因此对比损伤前后的应变或者应变响应参数,可以用来识别结构损伤.提出了一种应变脉冲响应协方差参数,它是应变脉冲响应在时间区间上的能量积分;推导并证明了该参数是结构模态参数(频率,位移模态,应变模态,阻尼等)的函数,可用来表征结构状态.相比于传统的模态参数识别方法,可以保留更高阶的模态参数,而且避免了模态识别可能引起的误差;基于简支钢梁的多种损伤工况,研究和展示了该参数的特性,通过数值模拟发现,该参数能简单直观地判定损伤发生和识别损伤位置,无需建立结构分析模型,只需比较结构损伤前后的应变脉冲响应协方差参数即可;该参数简便易算,具有较好的抗噪性能,对结构损伤敏感,而且对结构刚度减少呈现一致变化特性,所以适合实际工程结构的健康监测和损伤识别.     

静气动弹性计算方法研究

对基于结构网格的Euler方程及N-S方程求解器和基于非结构网格的Euler方程求解器,采用结构模态分析方法和柔度矩阵方法,对无人机大展弦比机翼在Ma=0.6,α=2?, 飞行高度20km的巡航状态下的静气动弹性特性进行了数值模拟. 验证了两种求解器对静气动弹性模拟的准确性. 同时,对模态分析方法和柔度矩阵方法进行了对比研究,发现柔度矩阵方法更适用于静气动弹性数值模拟. 另外,对应用物面法向偏转方法替代网格变形技术模拟静气动弹性进行了研究,计算表明物面法向偏转方法可以大大提高静气动弹性计算效率和克服机翼结构变形过大时动网格技术无法处理的不足.      

求解非比例阻尼体系复模态的实模态摄动法

根据工程结构的实际情况,建立了非比例阻尼结构体系复模态特性的近似求解方 法------实模态摄动法. 这一方法以复Ritz向量展开原理为基础, 把非比例阻尼结构体系复模态特性的分析过程分解为两个基本步骤,首先以结构体系的实模 态向量构建复Ritz向量的求解子空间,然后通过非线性复代数方程组的求解代替扩阶后的复 特征值方程的求解,从而简化了计算过程.通过两个算例表明: 这一方法不仅计算简便,而且具有较高的计算精度和执行效率, 对于复杂的非比例阻尼系统是很适用的,具有一定的工程应用价值.     

结构静态拓扑重分析的摄动-Pade逼近法

提出一种拓扑修改的静态重分析的方法．该方法将所有新增加的自由度分步、逐次得增加到原结构上．在每一子步中利用过渡矩阵得到摄动基，并对摄动基进行正交化处理来进行Pade逼近，得到的近似解作为下一子步的原始解．循环结束时得到修改结构的近似解．     

Study of the anisotropy of sliding stops of mobile robots

The interaction of mini-robot contact elements with the surface is studied. The values of basic parameters necessary for successful operation of such elements are estimated. A procedure for calculating the value of anisotropic friction typical to adhesive material with haired structure is proposed. The results of this paper can be used to design robots intended for motion in bounded spaces, in particular, in narrow tubes.     

Limit Analysis of Structures with Stochastic Strengths by a Static Approach

A procedure for the evaluation of the conditional probability of collapse of elastic–plastic structures with stochastic strengths is presented. The procedure represents an extension to the stochastic framework of the static approach, coupled with the method of redundant unknowns, well established in the classical deterministic limit analysis. In this paper, a kinematic approach for probabilistic limit analysis, provided in the literature, is also recalled in order to show the duality with the proposed static approach.A beam of elastic perfectly plastic material with correlated stochastic strengths is studied and the influence of the degree of correlation between cross-sections is evidenced. When a low correlation is encountered an accurate discretisation is necessary. Furthermore, investigation of the influence of different collapse events on the conditional probability of collapse is conducted for a frame structure.     

两相介质饱和土三维非轴对称稳态响应分析

用积分变换方法分析饱和土三维非轴对称波动方程，用刚度矩阵方法分析层状饱和土三维非轴对称稳态响应。以数值算例对比分析了单相土介质与两相饱和土介质三维非轴对称稳态动力响应。分析中未引入势函数及传统的Helmholtz分解，推导过程及解答均非常简捷，为边界元方法在三维非轴对称饱和土与基础结构动力相互作用研究中的应用奠定基础。     

基于改进粒子群的薄壁变截面刚架临界载荷优化算法

针对大型变截面薄壁结构的稳定问题,以一类任意约束对称结构受非对称载荷的单跨刚架为研究对象,结构拆分为相关铁木辛柯（Timoshenko）梁,结合差分原理和最优化方法,以每段刚架的每个离散点挠度、临界载荷、轴力、剪力和梁端弯矩为设计变量,建立求解满足边界条件的非线性差分方程模型,提出基于优胜劣汰粒子更新的粒子群（IPSO）临界载荷优化算法。运用JAVA编程语言编制对应优化程序,分析典型算例并核实ABAQUS仿真结果。研究表明,本文提出的优化算法获得了有效的变形位型和高精度的临界载荷计算,能更好地描述刚架受力下位型和载荷的力学关系,进一步为工程设计与分析提供支持。