多层圆柱体锚固结构中超声导波频散特征的研究

超声导波无损检测技术因其高效和快捷的优点成为检测锚杆锚固质量的有效方法。但锚固锚杆结构中超声导波的多模态性、频散性与能量泄露导致完整的锚杆底端反射信号的获得具有不确定性。应用弹性动力学理论并采用全局矩阵法建立超声导波在多层圆柱体锚固结构中的频散方程的通用表达式,然后通过非线性外推法和二分法两步算法求解频散方程的精确解,解决了频散曲线的分类和相交等难题,获得了具有自主知识产权的求解多层圆柱体锚固结构频散曲线的程序。利用该程序计算了不同锚固锚杆结构的频散曲线,并与商用软件Disperse计算的结果吻合较好。同时用本程序计算了实验室有限锚固结构中的频散曲线,验证了低频导波在有限与无限结构中的巨大差异,而这一结构特征的频散特性Disperse并未给出相关算例。     

磁场对电磁弹性板中Lamb波频散特性的影响

研究了静态磁场下电磁弹性结构中Lamb波的传播行为.在确定静态磁场下电磁弹性板中耦合初始广义应力(弹性应力、电位移和磁感应强度的基础上,推导了含初始广义应力时板中Lamb波传播的运动方程,并由此获得了对称模态和反对称模态时的频散方程.以由BaTiO3-CoFe2O4材料构成的电磁弹性板模型作为数值算例,绘制了Lamb波传播的对称模态和反对称模态频散曲线.计算结果表明磁场对Lamb波的频散特性有一定的影响.     

考虑粘弹性材料频变特性的复合结构频率响应分析

论文采用基于Layerwise离散层理论的有限元板单元建模粘弹性阻尼复合结构,发展了一种考虑粘弹性材料的频变特性的复合结构频率响应计算方法.该方法首先计算结构的质量矩阵和各个频率点下的结构复刚度矩阵,然后求解运动方程的若干阶复特征对并由此计算各阶模态的传递函数,最后将各模态传函线性叠加得到近似的总传函.为了提高计算效率,论文采用了一种高效的数值方法,即只计算若干频率点下的特征向量与特征值,并计算这些点处特征向量关于频率的高阶导数,通过泰勒系数展开逼近和Rayleigh商式,可求得附近若干频率点处的特征向量和特征值,从而避免了在各个频率点下求解大自由度结构特征方程的问题,可以极大地提高计算效率.对一端简支的三层约束阻尼梁算例进行了分析,并与文献中的结果作对比,结果验证了方法的有效性和计算效率.     

波导本征问题分析的比例边界有限元方法

引入了一种求解波导本征值问题的高效而精确算法-比例边界有限元方法SBFEM (Scaled Boundary Finite Element Method).该方法的一个特点是只需在边界上进行离散,问题降低一维,使计算工作量大大减少;另一特点是所建立的控制方程为二阶常微分方程,可以解析地求解,使计算精度得到了保证.论文利用变分原理并通过比例边界坐标变换,推导了TE波和TM波波导的比例边界有限元频域方程以及波导动剐度方程,同时给出了波导动刚度矩阵的连分式解形式,通过引入辅助变量进一步得出波导特征值方程并求出波导本征值.以矩形、L形波导和叶型加载矩形波导的本征问题分析为例,通过与解析解及其他数值方法比较,结果表明,此方法具有精度高、计算工作量小的优点,而且随着连分式阶数增加收敛速度快.进一步分析了一类角切四脊正方形波导的传输特性.     

基于谱元法考虑层间接触状态的Rayleigh波频散特征正演计算

为了提升Rayleigh波应用于成层土基勘探的精细化水平,本文基于谱元法原理,通过Goodman模型和Matsui层间滑移系数建立了可考虑层间不同接触状态的Rayleigh波理论频散方程。针对典型土基成层结构,运用谱元法和快速矢量传递解析法对比计算了层间完全连续状态下土基的Rayleigh波多阶模态频散曲线,结果显示谱元法计算结果与解析法相应结果之间的平均相对误差在0.3%以下,具有较高的计算精度。在此基础上,通过改变层间接触状态和敏感性分析,揭示了层间接触状态对Rayleigh波基阶频散特征的影响。最后,结合速度-应力有限差分数值计算,验证了谱元法计算层间不同接触状态下R波频散特征的可靠性。     

A p-TYPE SUPERCONVERGENT RECOVERY METHOD FOR FE ELASTIC STABILITY ANALYSIS OF EULER BEAMS1)

本文将有限元p型超收敛算法应用于欧拉梁弹性稳定分析。该法基于有限元解答中失稳载荷和失稳模态结点位移的超收敛特性,建立了单元上失稳模态近似满足的线性常微分方程边值问题,在每个单元上,对该边值问题采用一个高次元进行求解,获得失稳模态的超收敛解,再将失稳模态的超收敛解代入瑞利商的解析表达式,最终获得失稳载荷的超收敛解。该法思路简明,通过少量计算即可显著提高失稳载荷和失稳模态的精度与收敛阶。数值算例表明,该法高效、可靠,值得进一步研究和推广到各类杆系结构。     

欧拉梁弹性稳定分析的有限元p型超收敛算法

本文将有限元p型超收敛算法应用于欧拉梁弹性稳定分析。该法基于有限元解答中失稳载荷和失稳模态结点位移的超收敛特性,建立了单元上失稳模态近似满足的线性常微分方程边值问题,在每个单元上,对该边值问题采用一个高次元进行求解,获得失稳模态的超收敛解,再将失稳模态的超收敛解代入瑞利商的解析表达式,最终获得失稳载荷的超收敛解。该法思路简明,通过少量计算即可显著提高失稳载荷和失稳模态的精度与收敛阶。数值算例表明,该法高效、可靠,值得进一步研究和推广到各类杆系结构。     

端点位移激励下斜拉索非线性振动计算方法研究

考虑拉索不同阶模态大幅振动之间的耦合效应,根据拉索的振动理论,详细地推导了单根拉索在端点位移激励下发生大幅振动时的非线性振动方程。根据某实际斜拉桥拉索参数,讨论了不同垂跨比对拉索振动特性的影响。使用四阶Runge-Kutta法求解拉索的非线性振动方程,通过对比有限元模型的非线性动力时程积分数值计算结果,验证了理论模型的可靠性与适用性。     

端点位移激励下斜拉索非线性振动计算方法研究Investigation on computing method of nonlinear vibration of the stayed cable with displacement excitation at the end point

考虑拉索不同阶模态大幅振动之间的耦合效应,根据拉索的振动理论,详细地推导了单根拉索在端点位移激励下发生大幅振动时的非线性振动方程。根据某实际斜拉桥拉索参数,讨论了不同垂跨比对拉索振动特性的影响。使用四阶Runge-Kutta法求解拉索的非线性振动方程,通过对比有限元模型的非线性动力时程积分数值计算结果,验证了理论模型的可靠性与适用性。     

考虑频率依赖性的涂层复合结构固有特性求解

粘弹性阻尼材料的力学特性参数会随着频率的变化而改变,即具有频率依赖性,因而传统的动力学建模及分析方法不能满足实际涂层结构优化设计的需要。在简要介绍粘弹性阻尼材料频率依赖性的基础上,本文提出用特征向量增值法来求解涂层复合结构的固有特性,并详细推导了特征向量增值法的求解原理。由此,提出了特征向量增值法的计算流程,包括计算无阻尼系统的固有特性,用Fox and Kapoor或者Nelson方法计算复特征向量增量;用Rayleigh熵法求解复特征值。最后,以涂敷粘弹性阻尼材料的钛基薄板为例,求解了该复合结构的固有特性,并与经典的模态应变能法进行了比较,证明了所提方法的正确性。     

A wideband fast multipole boundary element method for half-space/plane-symmetric acoustic wave problems

This paper presents a novel wideband fast multipole boundary element approach to 3D half-space/planesymmetric acoustic wave problems.The half-space fundamental solution is employed in the boundary integral equations so that the tree structure required in the fast multipole algorithm is constructed for the boundary elements in the real domain only.Moreover,a set of symmetric relations between the multipole expansion coefficients of the real and image domains are derived,and the half-space fundamental solution is modified for the purpose of applying such relations to avoid calculating,translating and saving the multipole/local expansion coefficients of the image domain.The wideband adaptive multilevel fast multipole algorithm associated with the iterative solver GMRES is employed so that the present method is accurate and efficient for both lowand high-frequency acoustic wave problems.As for exterior acoustic problems,the Burton-Miller method is adopted to tackle the fictitious eigenfrequency problem involved in the conventional boundary integral equation method.Details on the implementation of the present method are described,and numerical examples are given to demonstrate its accuracy and efficiency.     

伸缩虚拟边界元法解二维Helmholtz外问题

以位势理论为基础，提出了求解Helmholtz外问题的伸缩虚拟边界元法．给出了该方法在全波数域内获得唯一解的严格数学证明，其核心是通过伸缩虚拟边界使对偶内问题的特征频率(本征值)避开与波数重合，从而保证了解的唯一性，同以往前人提出的几种解法途径相比，该法简单得多；通过诸多边界曲线形状和不同边界量的声辐射算例，从计算精度、稳定性以及克服解的非唯一性等方面，对该方法进行了检验．计算结果表明：对远场或近场辐射声压，该方法都具有非常高的效率和精度．     

A modified Galerkin/finite element method for the numerical solution of the Serre‐Green‐Naghdi system

A new modified Galerkin/finite element method is proposed for the numerical solution of the fully nonlinear shallow water wave equations. The new numerical method allows the use of low‐order Lagrange finite element spaces, despite the fact that the system contains third order spatial partial derivatives for the depth averaged velocity of the fluid. After studying the efficacy and the conservation properties of the new numerical method, we proceed with the validation of the new numerical model and boundary conditions by comparing the numerical solutions with laboratory experiments and with available theoretical asymptotic results. Copyright © 2016 John Wiley & Sons, Ltd.     

Neural networks for BEM analysis of steady viscous flows

This paper presents a new neural network‐boundary integral approach for analysis of steady viscous fluid flows. Indirect radial basis function networks (IRBFNs) which perform better than element‐based methods for function interpolation, are introduced into the BEM scheme to represent the variations of velocity and traction along the boundary from the nodal values. In order to assess the effect of IRBFNs, the other features used in the present work remain the same as those used in the standard BEM. For example, Picard‐type scheme is utilized in the iterative procedure to deal with the non‐linear convective terms while the calculation of volume integrals and velocity gradients are based on the linear finite element‐based method. The proposed IRBFN‐BEM is verified on the driven cavity viscous flow problem and can achieve a moderate Reynolds number of 1400 using a relatively coarse uniform mesh. The results obtained such as the velocity profiles along the horizontal and vertical centrelines as well as the properties of the primary vortex are in very good agreement with the benchmark solution. Furthermore, the secondary vortices are also captured by the present method. Thus, it appears that an ability to represent the boundary solution accurately can significantly improve the overall solution accuracy of the BEM. Copyright © 2003 John Wiley & Sons, Ltd.     

Coupled frequencies of a rectangular hydroelastic system with two fluids

This paper deals with the study of behaviour of an idealized 2D hydroelastic system involving two inviscid liquids with an elastic rectangular container. The main objective is to investigate the influence of the physical parameters on eigenfrequencies and eigenmodes of the system. The study extends the previous results obtained for hydroelastic systems with one fluid. The governing equations describing the behaviour of the system are analyzed by using the concept of normal modes and their solutions presented in the form of infinite series. The expansion coefficients for the velocity potentials are calculated by employing a new inner product that allows the orthogonalisation of the normal modes. An eigenfrequency equation is then derived from the existence condition of a nontrivial solution. The numerical calculations are performed by varying only some relevant parameters.     

Solution of the Navier–Stokes equations in velocity–vorticity form using a Eulerian–Lagrangian boundary element method

This paper describes the Eulerian–Lagrangian boundary element model for the solution of incompressible viscous flow problems using velocity–vorticity variables. A Eulerian–Lagrangian boundary element method (ELBEM) is proposed by the combination of the Eulerian–Lagrangian method and the boundary element method (BEM). ELBEM overcomes the limitation of the traditional BEM, which is incapable of dealing with the arbitrary velocity field in advection‐dominated flow problems. The present ELBEM model involves the solution of the vorticity transport equation for vorticity whose solenoidal vorticity components are obtained iteratively by solving velocity Poisson equations involving the velocity and vorticity components. The velocity Poisson equations are solved using a boundary integral scheme and the vorticity transport equation is solved using the ELBEM. Here the results of two‐dimensional Navier–Stokes problems with low–medium Reynolds numbers in a typical cavity flow are presented and compared with a series solution and other numerical models. The ELBEM model has been found to be feasible and satisfactory. Copyright © 2000 John Wiley & Sons, Ltd.     

边界元特征值分析及Burton-Miller法探究BEM eigenvalue analysis and the investigation of the Burton-Miller method

运用围道积分方法将边界元非线性特征值问题转化为规模很小的广义特征值问题,从而构造出一种边界元特征值分析方法。数值算例验证了该方法的求解精度。针对外声场问题,通过对常规、法向导数和Burton‐Miller边界积分方程的虚假特征频率的计算和比较,揭示了Burton‐Miller法规避虚假特征频率的本质,并对其中的叠加常数的最优取值给出了一种新的解释。     

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

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

ADAPTIVE LINEARIZATION AND GRID ITERATIONS WITH THE TRI-TREE MULTIGRID REFINEMENT–RECOARSEMENT ALGORITHM FOR THE NAVIER–STOKES EQUATIONS

The tri-tree algorithm for refinements and recoarsements of finite element grids is explored. The refinement–recoarsement algorithm not only provides an accurate solution in certain parts of the grid but also has a major influence on the finite element equation system itself. The refinements of the grid lead to a more symmetric and linear equation matrix. The recoarsements will ensure that the grid is not finer than is necessary for preventing divergence in an iterative solution procedure. The refinement–recoarsement algorithm is a dynamic procedure and the grid is adapted to the instant solution. In the tri-tree multigrid algorithm the solution from a coarser grid is scaled relatively to the increase in velocity boundary condition for the finer grid. In order to have a good start vector for the solution of the finer grid, the global Reynolds number or velocity boundary condition should not be subject to large changes. For each grid and velocity solution the element Reynolds number is computed and used as the grid adaption indicator during the refinement–recoarsement procedure. The iterative tri-tree multigrid method includes iterations with respect to the grid. At each Reynolds number the same boundary condition s are applied and the grid is adapted to the solution iteratively until the number of unknowns and elements in the grid becomes constant. In the present paper the following properties of the tri-tree algorithm are explored: the influence of the increase in boundary velocities and the size of the grid adaption indicator on the amount of work for solving the equations, the number of linear iterations and the solution error estimate between grid levels. The present work indicates that in addition to the linear and non-linear iterations, attention should also be given to grid adaption iterations. © 1997 by John Wiley & Sons, Ltd.