A symmetric BEM–FEM method for an axisymmetric eddy current problem

This paper deals with the numerical solution of the time-harmonic eddy current model in an axisymmetric unbounded domain. To this end, a new symmetric BEM–FEM formulation is derived and also analyzed. Moreover, error estimates for the corresponding discretization are proven.     

高分辨率空间相机的工程分析

以大型结构设计软件SolidEdge建立了高分辨率空间相机的精确的几何模型, 在此基础上用大型有限元分析软件MSC.Patran建立了相机的有限元模型, 并利用MSC/Nastran对该相机整机进行了工程分析. 其计算结果与试验结果相差10%以内.     

深基坑工程渗流场特性分析

深基坑工程的施工特点决定了其渗流场水力条件的复杂性。大量的工程实践显示，渗流问题是许多基坑工程事故的主要或重要原因。本文采用有限单元法分析了基坑渗流场分布特性。比较了不同水力条件下渗流作用对基坑土体渗透稳定性的影响，探讨了工程中可能出现的不利因素，特别是防渗体破坏情况对基坑安全造成的危害性。本文的工作也表明，利用数值方法分析基坑渗流场分布特性，对于深基坑工程的安全性评价和工程施工都具有重要的指导意义。     

Blended kernel approximation in the ℋ︁‐matrix techniques

Several types of ??‐matrices were shown to provide a data‐sparse approximation of non‐local (integral) operators in FEM and BEM applications. The general construction is applied to the operators with asymptotically smooth kernel function provided that the Galerkin ansatz space has a hierarchical structure. The new class of ??‐matrices is based on the so‐called blended FE and polynomial approximations of the kernel function and leads to matrix blocks with a tensor‐product of block‐Toeplitz (block‐circulant) and rank‐k matrices. This requires the translation (rotation) invariance of the kernel combined with the corresponding tensor‐product grids. The approach allows for the fast evaluation of volume/boundary integral operators with possibly non‐smooth kernels defined on canonical domains/manifolds in the FEM/BEM applications. (Here and in the following, we call domains canonical if they are obtained by translation or rotation of one of their parts, e.g. parallelepiped, cylinder, sphere, etc.) In particular, we provide the error and complexity analysis for blended expansions to the Helmholtz kernel. Copyright © 2002 John Wiley & Sons, Ltd.     

粘性流动的局部隐式有限元解法

粘性流动的局部隐式有限元解法依凤鸣，吴遇春（哈尔滨工业大学动力系哈尔滨１５Ｏ００１）关键词粘性流动，局部隐式有限元法，Ｎ．Ｓ．方程１引言有限元法具有对复杂几何边界适应性强，对流场内变量梯度大的区域计算精度高以及易于编制大型通用程序等特点，已经成了叶轮...     

Effect of a characteristic length on crack spacing in a reinforced concrete bar under tension

The paper presents quasi-static FE-simulations of the crack formation in a reinforced concrete bar without stirrups subject to tension. The material was modeled with a continuum smeared crack model using an elasto-plastic constitutive law. A linear Rankine criterion with isotropic softening and associated flow rule was adopted in a tensile regime. To ensure the mesh-independency, the softening parameter was enhanced by a characteristic length of micro-structure by means of a non-local theory. Attention was laid to the effect of a different characteristic length of micro-structure and initial bond-slip stiffness on the spacing of localized zones.     

Lyapunov spectrum determination from the FEM simulation of a chaotic advecting flow

The problem of the determination of the Lyapunov spectrum in chaotic advection using approximated velocity fields resulting from a standard FEM method is investigated. A fourth order Runge–Kutta scheme for trajectory integration is combined with a third order Jacobian matrix method with QR ‐factorization. After checking the algorithm on the standard Lorenz and coupled quartic oscillator systems, the method is applied to a model 3‐D steady flow for which an analytical expression is known. Both linear and quadratic approximated velocity fields succeed in predicting the Lyapunov exponents as well as describing the chaotic or regular regions inside the flow with satisfactory accuracy. A more realistic flow is then studied in order to delineate the possible limitations of the approach. Copyright © 2005 John Wiley & Sons, Ltd.     

Time-harmonic acoustic propagation in the presence of a shear flow

This work deals with the numerical simulation, by means of a finite element method, of the time-harmonic propagation of acoustic waves in a moving fluid, using the Galbrun equation instead of the classical linearized Euler equations. This work extends a previous study in the case of a uniform flow to the case of a shear flow. The additional difficulty comes from the interaction between the propagation of acoustic waves and the convection of vortices by the fluid. We have developed a numerical method based on the regularization of the equation which takes these two phenomena into account. Since it leads to a partially full matrix, we use an iterative algorithm to solve the linear system.     

芯片封装焊球连接疲劳寿命预测分析--能量法和有效应变法之比较

研究了目前常用的两种芯片封装焊球连接的寿命预测方法一能量法和有效应变法。对焊球连接的材料本构关系进行了对比分析，并利用FORTRAN语言编制了相应的材料本构模型子程序，将其与有限元仿真工具ANSYS耦合，实现了将焊球连接的材料本构模型用户子程序导入到ANSYS的材料库。在此基础上模拟了三维芯片叠层球栅阵列尺寸封装的焊球结构，在热循环条件下(-40～ 125℃的工作状态，并分别利用能量法和有效应变法对焊球连接的寿命预测进行了比较分析，最后对两种方法作出了评价。     

Study on Combined Method Based on 3-D ESPI

The finite element method (FEM),whether the calculation is accurate or not,depends closely on object boundary condition.If the three dimensional displacement of the object obtained in experiment is regarded as its boundary condition,a new method combining the results of experiment and calculation,called combined method (CM),is formed.The combined method possess advantages of experiment and calculation.It can correct calculation and improve the accuracy of FEM.Accordingly it has more practicability.In this paper,the three dimensional displacement fields of a typical beam loaded at three points are tested by using 3-D electric speckle pattern interferometry (ESPI).Using the experimental results as boundary condition the whole three-dimensional displacement fields can be calculated by FEM.The beam′s three-dimensional displacement fields obtained by FEM agree very well with those obtained by experiment.This proves that the combined method is effective and practicable.