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.     

基于哈密顿体系的平面无限解析元

在有限元法中,无限域的问题不便于处理求解。但无限域往往可以由规则的无限外域再加上有限的局部域组成。将无限域问题中的有限局部域用有限元法处理,在规则的无限外域中建立极坐标系,将规则无限域问题导向哈密顿体系,利用本征向量展开的方法,推导出一种新的半解析无限解析元,其刚度阵是精确的。该单元可用常规方法作为一个超级有限单元与有限的局部域连接。数值计算结果表明,该单元具有精度高,应用方便,数据处理非常简单的特点。对无限域问题的数值求解有重要意义。该方法可推广到三维无限域问题中。     

标准有限元方法模拟大气压流注等离子体中数值振荡与伪扩散问题研究与改进

基于商用软件的等离子体数值模拟方案大多是基于标准Galerkin格式的有限元方法,在计算大气压流注放电时,存在数值振荡和伪扩散的问题,计算结果可靠性不高.针对针-板体流注放电结构,构建了修正的COMSOL Multiphysics等离子体模型.修正模型与COMSOL Multiphysics内置等离子体模块的区别主要在于增加了基于弱形式的人工稳定项.利用修正模型复现了经典的雪崩-流注转捩算例,并与经典纳秒秒冲过电压针板放电实验结果进行了对比验证.结果表明,引入修正后的商用软件模型能够克服数值振荡和伪扩散缺陷,确保等离子体数值计算的准确性和可复现性.      

Superconvergent estimates of conforming finite element method for nonlinear time‐dependent Joule heating equations

In this article, we study the superconvergence analysis of conforming bilinear finite element method (FEM) for nonlinear Joule heating equations. Based on the rigorous estimates together with high accuracy analysis of this element, mean value technique and interpolation postprocessing approach, the superclose and superconvergent estimates about the related variables in H¹‐norm are derived for semidiscrete and a linearized backward Euler fully discrete schemes, which extends the results of optimal estimates obtained for conforming FEMs in the previous literature. At last, a numerical experiment is performed to verify the theoretical analysis.     

The coupling of bem and fem for the two-dimensional viscous flow problem

In this paper we propose a numerical scheme for treating the problem of sJow viscous flow past an obstacle in the plane. This scheme is a combination of boundary element and finite element methods. By introducing an auxiliary boundary curve, we divide the region under consideration into two subregions, an inner and an outer region. In the inner region, we employ a finite element method (FEM) for solving a system of simplified field equations with proper natural boundary conditions. In the outer region, the solution is expressed in the form of a simple-layer potential with density function satisfying a system of modified integral equations of the first kind. The latter are solved by a boundary element method (BEM). Both solutions are matched on the common auxiliary boundary curve. Error estimates in suitable function spaces are derived in terms of the mesh widths as well as the small parameters, the Reynolds numbers     

Laser forming of a bowl shaped surface with a stationary laser beam

Despite a lot of research done in the field of laser forming, generation of a symmetric bowl shaped surface by this process is still a challenge mainly because only a portion of the sheet is momentarily deformed in this process, unlike conventional sheet metal forming like deep drawing where the entire blank undergoes forming simultaneously reducing asymmetry to a minimum. The motion of laser beam also makes the process asymmetric. To counter these limitations this work proposes a new approach for laser forming of a bowl shaped surface by irradiating the centre of a flat circular blank with a stationary laser beam. With high power lasers, power density sufficient for laser forming, can be availed at reasonably large spot sizes. This advantage is exploited in this technique. Effects of duration of laser irradiation and beam spot diameter on the amount of bending and asymmetry in the formed surface were investigated. Laser power was kept constant while varying irradiation time. While varying laser spot diameter laser power was chosen so as to keep the surface temperature nearly constant at just below melting. Experimental conditions promoted almost uniform heating through sheet thickness. The amount of bending increased with irradiation time and spot diameter. It was interesting to observe that blanks bent towards the laser beam for smaller laser beam diameters and the reverse happened for larger spot diameters (~10 times of the sheet thickness). Effect of spot diameter variation has been explained with the help of coupled thermal-structural finite element simulations.     

Acoustic properties of micro-perforated panel absorber having arbitrary cross-sectional perforations

Micro-perforated panel absorber is used in many noise control applications as a next-generation absorbing material. Perforation shapes of micro-perforated panel studied are usually circular in the past. However, in practice, the perforations are often non-circular or irregular shape due to manufacturing techniques. Sound absorption coefficient and absorption bandwidth of the micro-perforated panel absorber may be further improved, when the perforations in shape are changed. In view of the existing exact solutions of sound propagation in tubes, the simple formulas of specific acoustic impedances of the tubes for triangle and square cross-sectional perforations are derived. Mass reactance end correction of the micro-perforated panel is obtained based on the sound radiation of a shaped piston. The specific acoustic impedance ratio of the micro-perforated panel absorber is calculated and analyzed, which can predict its sound absorption bandwidth. Finally, for closed perforations, the influences of the perforations in shape (including triangle, circle, square and irregular circle) on sound absorption of the MPP absorber are discussed in collaboration with FE simulations.     

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

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

固体材料脆性断裂的多尺度时空模型(英文)

通过实例分析了动态断裂过程所固有的典型动态作用,基于断裂孵化时间理论的时空结构方法提出了固体材料断裂的统一解释.此外,还提出了计及尺度效应的广义断裂模型,该模型可以用于预测准脆性非均质材料的多尺度断裂,并证明此模型基于实验室尺度的实验数据可以预测更高尺度(真实尺度)的宏观断裂.     

A refined numerical model for determining inflation-burst behavior of composite membrane structures

The airship structures made of multi-layer composite fabrics or membranes can offer the platform for earth observations, wireless communications and space research due to light weight and good mechanical performance. The structural safety and serviceability strongly depend on material properties and working conditions. Available studies are limited within service stress limits or are lack of suitable biaxial tensile constitutive models for understanding structural behavior. This paper thus focuses on a refined numerical model for determining inflation-burst behavior of composite airship structures considering new biaxial constitutive equations, novel failure criteria and manufacture factors.The differences between ideal and real forms of airship structures, e.g. volume difference, demonstrate the necessity for incorporating cutting-pattern effects in the initial numerical model. For structural analysis, stress distributions on real structural forms are different from those on ideal forms because of welding parts that can enhance local stiffness. The ultimate pressures are 56.7 kPa and 59.5 kPa for ideal and real structural forms. Structural breaking initiated at the maximum diameter of ideal structural forms propagates fast while welding parts can prevent breaking propagation for real structural forms. Therefore, the refined numerical model can reveal basic structural behavior and safety performance of airship structures in the inflation-burst processes.