Solving the homogeneous isotropic linear elastodynamics equations using potentials and finite elements. The case of the rigid boundary condition

In this article, elastic wave propagation in a homogeneous isotropic elastic medium with a rigid boundary is considered. A method based on the decoupling of pressure and shear waves via the use of scalar potentials is proposed. This method is adapted to a finite element discretization, which is discussed. A stable, energy preserving numerical scheme is presented, as well as 2D numerical results.     

Coupled FEM-BEM for elastoplastic contact problems using Lagrange multipliers

The coupling of the elastoplastic finite element and elastic boundary element methods for two-dimensional frictionless contact stress analysis is presented. Interface traction matching (boundary element approach), which involves the force terms in the finite element analysis being transformed to tractions, is chosen for the coupling method. The analysis at the contact region is performed by the finite element method, and the Lagrange multiplier approach is used to apply the contact constraints. Since the analyses of elastoplastic problems are non-linear and involve iterative solution, the reduced size of the final system of equations introduced by combining the two methods is very advantageous, especially for contact problems where the nature of the problem also involves an iterative scheme.     

A modified scaled boundary finite element method for dynamic response of a discontinuous layered half-space

In this paper, a modified scaled boundary finite element method is proposed to deal with the dynamic analysis of a discontinuous layered half-space. In order to describe the geometry of discontinuous layered half-space exactly, splicing lines, rather than a point, are chosen as the scaling center. Based on the modified scaled boundary transformation of the geometry, the Galerkin's weighted residual technique is applied to obtain the corresponding scaled boundary finite element equations in displacement. Then a modified version of dimensionless frequency is defined, and the governing first-order partial differential equations in dynamic stiffness with respect to the excitation frequency are obtained. The global stiffness is obtained by adding the dynamic stiffness of the interior domain calculated by a standard finite element method, and the dynamic stiffness of far field is calculated by the proposed method. The comparison of two existing solutions for a horizontal layered half-space confirms the accuracy and efficiency of the proposed approach. Finally, the dynamic response of a discontinuous layered half-space due to vertical uniform strip loadings is investigated.     

Three-dimensional simulation of flat rolling through a combined finite element–boundary element approach

This paper presents an innovative approach for analysing three-dimensional flat rolling. The proposed approach is based on a solution resulting from the combination of the finite element method with the boundary element method. The finite element method is used to perform the rigid–plastic numerical modelling of the workpiece allowing the estimation of the roll separating force, rolling torque and contact pressure along the surface of the rolls. The boundary element method is applied for computing the elastic deformation of the rolls. The combination of the two numerical methods is made using the finite element solution of the contact pressure along the surface of the rolls to define the boundary conditions to be applied on the elastic analysis of the rolls. The validity of the proposed approach is discussed by comparing the theoretical predictions with experimental data found in the literature.     

Extension of the scaled boundary finite element method to geometrical and physical nonlinearity

The scaled boundary finite element method (SBFEM) has been used in many fields of engineering to solve the governing equations in bounded and unbounded 2D as well as 3D domains. In solid mechanics, the semi-analytical solution strategy of the SBFE formulation (numerical in circumferential direction, analytical in radial direction) is based on the assumption of linear elastic material behavior and only small geometrical changes. However, a large group of materials (e.g. rubber) shows geometrical and physical nonlinearity at mechanical loading. In this contribution, the extension of the SBFEM to geometrical and physical nonlinearity is examined. A plane finite element is developed which uses the concept of shape functions constructed by the SBFEM in the framework of a nonlinear finite element analysis. (© 2014 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Analysis of free-edge effects by boundary finite element method

The scaled boundary finite element method is a novel semi-analytical analysis technique that combines the advantages of the finite element method and the boundary element method. Only a part of the boundary of the considered domain has to be discretized but nevertheless the method is solely finite element based. The governing equations are solved in the so-called scaling direction analytically, whereas a finite element approximation of the solution is performed in the circumferential directions, which form the boundary of the considered domain. Thus, the numerical effort can be reduced considerably when handling stress concentration problems such as e.g. the free-edge effect in laminated plates. In order to analyze the free-edge effect in a semi-infinite half plane, some kinematic coupling equations have to be introduced, that not only couple the degrees of freedom on the boundary, but also within the non-discretized domain. The implementation of kinematic coupling equations within the method is presented. Finally, the efficiency of the new approach is shown in some benchmark examples. (© 2006 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

A numerical analysis of the Cahn–Hilliard equation with non-permeable walls

In this article we consider the numerical analysis of the Cahn–Hilliard equation in a bounded domain with non-permeable walls, endowed with dynamic-type boundary conditions. The dynamic-type boundary conditions that we consider here have been recently proposed in Ruiz Goldstein et al. (Phys D 240(8):754–766, 2011) in order to describe the interactions of a binary material with the wall. The equation is semi-discretized using a finite element method for the space variables and error estimates between the exact and the approximate solution are obtained. We also prove the stability of a fully discrete scheme based on the backward Euler scheme for the time discretization. Numerical simulations sustaining the theoretical results are presented.     

Analysis of notches and cracks in circular Kirchhoff plates using the scaled boundary finite element method

A new formulation of the scaled boundary finite element method (SBFEM) is presented for the analysis of circular plates in the framework of Kirchhoff's plate theory. Essential for the SBFEM is, that a domain is described by the mapping of its boundary with respect to a scaling centre. The governing partial differential equations are transformed into scaled boundary coordinates and are reduced to a set of ordinary differential equations, which can be solved in a closed-form analytical manner. If the scaling centre is selected at the root of an existent crack or notch, the SBFEM enables the effective and precise calculation of singularity orders of cracked and notched structures. The element stiffness matrices for bounded and unbounded media are derived. Numerical examples show the performance and efficiency of the method, applied to plate bending problems. (© 2012 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Second Order Unconditionally Stable and Convergent Linearized Scheme for a Fluid-Fluid Interaction Model

In this paper, a fully discrete finite element scheme with second-order temporal accuracy is proposed for a fluid-fluid interaction model, which consists of two Navier-Stokes equations coupled by a linear interface condition. The proposed fully discrete scheme is a combination of a mixed finite element approximation for spatial discretization, the second-order backward differentiation formula for temporal discretization, the second-order Gear's extrapolation approach for the interface terms and extrapolated treatments in linearization for the nonlinear terms. Moreover, the unconditional stability is established by rigorous analysis and error estimate for the fully discrete scheme is also derived. Finally, some numerical experiments are carried out to verify the theoretical results and illustrate the accuracy and efficiency of the proposed scheme.     

A partitioned finite element scheme based on Gauge-Uzawa method for time-dependent MHD equations

In this paper, we mainly introduce a partitioned scheme based on Gauge-Uzawa finite element method for the 2D time-dependent incompressible magnetohydrodynamics (MHD) equations. It is a fully decoupled projection method which combines the Gauge and Uzawa methods within a variational formulation. Firstly, the temporal discretization is based on backward Euler technique for the linear term and semi-implicit scheme for the nonlinear term. Secondly, the spatial approximation of fluid velocity, hydrodynamic pressure, and magnetic field apply the mixed element method. Finally, the validity, reliability, and accuracy of the proposed algorithms are supported by numerical experiments.     

Elastic wave radiation and scattering in heterogeneous soil using the scaled boundary finite element method

A time-domain approach for the simulation of elastic waves in heterogeneous soil domains is presented. It is based on modelling both near and far field by the scaled boundary finite element method (SBFEM). The SBFEM facilitates the use of a structured mesh in the near field region without the need to circumvent hanging nodes. The quadtree mesh is obtained automatically from image data. Radiation damping in the far field is modelled accurately by means of a displacement unit-impulse-response-based formulation. An example analysis of wave radiation by an alluvial basin illustrates the potential of the proposed methodology. (© 2017 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

应用边界元法模拟纤维增强复合材料平面弹性问题

将含有随机分布多种夹杂相复合材料的二维弹性力学问题归结为复连通区域的边界积分方程,进而转化成矩阵方程进行求解和分析．根据同类夹杂相外在边界上的面力与位移之间关系矩阵完全相同的特点,使得最后的矩阵方程阶数得到大规模减少,这正是此处提出改进的边界元方法的主要思路．数值算例表明,对于此类问题,与常规的边界元分域解法相比更加有效．以该方法为基础,可以详细给出纤维增强复合材料二维条件下的宏观等效力学性质．     

The comparison between the DRBEM and DQM solution of nonlinear reaction–diffusion equation

In this study, both the dual reciprocity boundary element method and the differential quadrature method are used to discretize spatially, initial and boundary value problems defined by single and system of nonlinear reaction–diffusion equations. The aim is to compare boundary only and a domain discretization method in terms of accuracy of solutions and computational cost. As the time integration scheme, the finite element method is used achieving solution in terms of time block with considerably large time steps. The comparison between the dual reciprocity boundary element method and the differential quadrature method solutions are made on some test problems. The results show that both methods achieve almost the same accuracy when they are combined with finite element method time discretization. However, as a method providing very good accuracy with considerably small number of grid points differential quadrature method is preferrable.     

Semi-analytical hybrid approach for the simulation of layered waveguide with a partially debonded piezoelectric structure

This paper presents an efficient semi-analytical hybrid approach for simulating the dynamic interaction of perfectly bonded or damaged piezoelectric structures with a layered elastic waveguide. In the proposed approach, the frequency domain spectral element method is utilized for the discretization of the finite-sized surface mounted piezoelectric structure, and the semi-analytical boundary integral method is employed for the evaluation of wave phenomena in the host laminate structure. While the spectral element method allows cost-effective simulation of dynamics of a complex-shaped transducer (e.g. curvilinear or with wrapped electrodes), the analytically-based technique reliably describes wave excitation and propagation in multi-layered structures. The coupling of these methods is achieved through the rigorous fulfillment of the boundary conditions at the area of waveguide-transducer contact. Three various combinations of approximation polynomials and surface-load interpolation functions are applied in order to obtain the solution in a frequency domain. The time-domain solution is evaluated employing the inverse Laplace transform. Convergence of the method is confirmed for different bonding conditions. The paper demonstrates the efficiency of the proposed method for the multi-parameter analysis of the dependence of the resonance characteristics on the debonding parameters and contact conditions. The approach can be used for such a crucial task as diagnosing failures of piezoelectric devices incorporated into a structural health monitoring system based on guided waves.     

ERROR ESTIMATES FOR TWO-SCALE COMPOSITE FINITE ELEMENT APPROXIMATIONS OF NONLINEAR PARABOLIC EQUATIONS

We study spatially semidiscrete and fully discrete two-scale composite finite element method for approximations of the nonlinear parabolic equations with homogeneous Dirich-let boundary conditions in a convex polygonal domain in the plane.This new class of finite elements,which is called composite finite elements,was first introduced by Hackbusch and Sauter[Numer.Math.,75(1997),pp.447-472]for the approximation of partial differential equations on domains with complicated geometry.The aim of this paper is to introduce an efficient numerical method which gives a lower dimensional approach for solving par-tial differential equations by domain discretization method.The composite finite element method introduces two-scale grid for discretization of the domain,the coarse-scale and the fine-scale grid with the degrees of freedom lies on the coarse-scale grid only.While the fine-scale grid is used to resolve the Dirichlet boundary condition,the dimension of the finite element space depends only on the coarse-scale grid.As a consequence,the resulting linear system will have a fewer number of unknowns.A continuous,piecewise linear composite finite element space is employed for the space discretization whereas the time discretization is based on both the backward Euler and the Crank-Nicolson methods.We have derived the error estimates in the L∞(L2)-norm for both semidiscrete and fully discrete schemes.Moreover,numerical simulations show that the proposed method is an efficient method to provide a good approximate solution.     

A High-Order Finite Difference Scheme for 3D Unsteady Convection Diffusion Reaction Equations北大核心CSCD

针对三维非稳态对流扩散反应方程,构造了一种高精度紧致有限差分格式,对空间的离散采用四阶紧致差分方法,对时间的离散采用Taylor级数展开和余项修正技术,所提格式在时间上的精度为二阶、在空间上的精度为四阶。利用Fourier稳定性分析法证明了该格式是无条件稳定的。最后给出数值算例验证了理论结果。     

Dynamical analysis of Kirchhoff plates by an explicit time integration Morley element method

An explicit time integration finite element method is proposed to investigate dynamical analysis of Kirchhoff plates, where the Morley element is used for spatial discretization and the second-order central scheme for time discretization. Certain error estimates in the energy norm are achieved. A number of numerical results are included to show computational performance of the method.     

Implicit finite element schemes for stationary compressible particle-laden gas flows

The derivation of macroscopic models for particle-laden gas flows is reviewed. Semi-implicit and Newton-like finite element methods are developed for the stationary two-fluid model governing compressible particle-laden gas flows. The Galerkin discretization of the inviscid fluxes is potentially oscillatory and unstable. To suppress numerical oscillations, the spatial discretization is performed by a high-resolution finite element scheme based on algebraic flux correction. A multidimensional limiter of TVD type is employed. An important goal is the efficient computation of stationary solutions in a wide range of Mach numbers. This is a challenging task due to oscillatory correction factors associated with TVD-type flux limiters and the additional strong nonlinearity caused by interfacial coupling terms. A semi-implicit scheme is derived by a time-lagged linearization of the nonlinear residual, and a Newton-like method is obtained in the limit of infinite CFL numbers. The original Jacobian is replaced by a low-order approximation. Special emphasis is laid on the numerical treatment of weakly imposed boundary conditions. It is shown that the proposed approach offers unconditional stability and faster convergence rates for increasing CFL numbers. The strongly coupled solver is compared to operator splitting techniques, which are shown to be less robust.     

CONVERGENCE OF THE CRANK-NICOLSON/NEWTON SCHEME FOR NONLINEAR PARABOLIC PROBLEM

In this paper, the Crank-Nicolson/Newton scheme for solving numerically secondorder nonlinear parabolic problem is proposed. The standard Galerkin finite element method based on P2 conforming elements is used to the spatial discretization of the problem and the Crank-Nicolson/Newton scheme is applied to the time discretization of the resulted finite element equations. Moreover, assuming the appropriate regularity of the exact solution and the finite element solution, we obtain optimal error estimates of the fully discrete CrankNicolson/Newton scheme of nonlinear parabolic problem. Finally, numerical experiments are presented to show the efficient performance of the proposed scheme.     

区间运算和静力区间有限元

用均值和离差两参数表征区间变量的不确定性,根据区间运算规则,论证了区间变量的运算特性.将区间分析和有限元方法相结合,提出了非概率不确定结构的一种区间有限元分析方法.将区间有限元静力控制方程中n自由度不确定位移场特征参数的求解归结为求解一2n阶线性方程组.实例分析表明文中方法是有效和可行的.