一种基于新型插值单元的稳态传热边界元法

为了提高边界元法在求解稳态热问题时的计算精度,通过使用一种新型单元插值方法(称为扩展单元插值法),实现对稳态传热问题的求解。扩展单元是在传统不连续单元的边界配置虚拟节点,把原非连续单元变成高阶的连续单元,并将其作为新型的插值单元。利用虚拟节点和内部源节点构造出的插值函数,可以精确插值边界上的连续和不连续物理场,插值精度要比原始不连续单元高两阶。另外,边界积分方程只在传统的不连续单元的内部节点处建立,只包含内部源节点的自由度,而虚拟节点的自由度可通过与内部源节点之间的关系消除掉,因此最终系统方程的求解规模不会增加。这种新型的插值单元继承了传统连续和不连续单元的优点,克服了它们的缺点。数值结果表明,此种单元插值方法用于求解稳态传热问题时可获得较高的计算精度和收敛性。     

Finite element and boundary element contact stress analysis with remeshing technique

The fundamental part of the contact stress problem solution using a finite element method is to locate possible contact areas reliably and efficiently. In this research, a remeshing technique is introduced to determine the contact region in a given accuracy. In the proposed iterative method, the meshes near the contact surface are modified so that the edge of the contact region is also an element’s edge. This approach overcomes the problem of surface representation at the transition point from contact to non-contact region. The remeshing technique is efficiently employed to adapt the mesh for more precise representation of the contact region. The method is applied to both finite element and boundary element methods. Overlapping of the meshes in the contact region is prevented by the inclusion of displacement and force constraints using the Lagrange multipliers technique. Since the method modifies the mesh only on the contacting and neighbouring region, the solution to the matrix system is very close to the previous one in each iteration. Both direct and iterative solver performances on BEM and FEM analyses are also investigated for the proposed incremental technique. The biconjugate gradient method and LU with Cholesky decomposition are used for solving the equation systems. Two numerical examples whose analytical solutions exist are used to illustrate the advantages of the proposed method. They show a significant improvement in accuracy compared to the solutions with fixed meshes.     

Integral representation of slowly growing equidistant splines

In this paper we consider polynomial splines S(x) with equidistant nodes which may grow as O (|x|^s). We present an integral representation of such splines with a distribution kernel. This representation is related to the Fourier integral of slowly growing functions. The part of the Fourier exponentials herewith play the so called exponential splines by Schoenberg. The integral representation provides a flexible tool for dealing with the growing equidistant splines. First, it allows us to construct a rich library of splines possessing the property that translations of any such spline form a basis of corresponding spline space. It is shown that any such spline is associated with a dual spline whose translations form a biorthogonal basis. As examples we present solutions of the problems of projection of a growing function onto spline spaces and of spline interpolation of a growing function. We derive formulas for approximate evaluation of splines projecting a function onto the spline space and establish therewith exact estimations of the approximation errors.     

Local mass conservative Hermite interpolation for the solution of flow problems by a multi-domain boundary element approach

This paper presents a local Hermite radial basis function interpolation scheme for the velocity and pressure fields. The interpolation for velocity satisfies the continuity equation (mass conservative interpolation) while the pressure interpolation obeys the pressure equation. Additionally, the Dual Reciprocity Boundary Element method (DRBEM) is applied to obtain an integral representation of the Navier-Stokes equations. Then, the proposed local interpolation is used to obtain the values of the field variables and their partial derivatives at the boundary of the sub-domains. This interpolation allows one to obtain the boundary values needed for the integral formulas for velocity and pressure at some nodes within the sub-domains. In the proposed approach the boundary elements are merely used to parameterize the geometry, but not for the evaluation of the integrals as it is usually done. The presented multi-domain approach is different from the traditional ones in boundary elements because the resulting integral equations are non singular and the boundary data needed for the boundary integrals are approximated using a local interpolation. Some accurate results for simple Stokes problems and for the Navier-Stokes equations at low Reynolds numbers up to Re = 400 were obtained.     

A higher order finite element including transverse normal strain for linear elastic composite plates with general lamination configurations

This paper describes a higher-order global-local theory for thermal/mechanical response of moderately thick laminated composites with general lamination configurations. In-plane displacement fields are constructed by superimposing the third-order local displacement field to the global cubic displacement field. To eliminate layer-dependent variables, interlaminar shear stress compatibility conditions have been employed, so that the number of variables involved in the proposed model is independent of the number of layers of laminates. Imposing shear stress free condition at the top and the bottom surfaces, derivatives of transverse displacement are eliminated from the displacement field, so that C⁰ interpolation functions are only required for the finite element implementation. To assess the proposed model, the quadratic six-node C⁰ triangular element is employed for the interpolation of all the displacement parameters defined at each nodal point on the composite plate. Comparing to various existing laminated plate models, it is found that simple C⁰ finite elements with non-zero normal strain could produce more accurate displacement and stresses for thick multilayer composite plates subjected to thermal and mechanical loads. Finally, it is remarked that the proposed model is quite robust, such that the finite element results are not sensitive to the mesh configuration and can rapidly converge to 3-D elasticity solutions using regular or irregular meshes.     

Spline Smoothing on Surfaces

This article presents a method for estimating functions on topologically and/or geometrically complex surfaces from possibly noisy observations. Our approach is an extension of spline smoothing, using a finite element method. The article has a substantial tutorial component: we start by reviewing smoothness measures for functions defined on surfaces, simplicial surfaces and differentiable structures on such surfaces, subdivison functions, and subdivision surfaces. After describing our method, we show results of an experiment comparing finite element approximations to exact smoothing splines on the sphere, and we give examples suggesting that generalized cross-validation is an effective way of determining the optimal degree of smoothing for function estimation on surfaces.     

Radially symmetric weighted extended b-spline model

In this paper, the weighted extended basis splines approach in the finite element method is applied to the electrostatic, electromagnetic wave and bioheat problems for inhomogeneous boundary conditions and radially symmetric structures. This new method, which does not need mesh generation, overcomes some of the drawbacks of using meshes and piecewise-uniform or linear trial functions. Two-dimensional radially symmetric electrostatic and electromagnetic wave equations are evaluated. We also attempt to propose a three-dimensional radially symmetric unexposed human eye model for simulating changes in corneal temperature using these new finite elements in conjunction with linear, quadratic and cubic b-splines. Our findings indicate that weighted extended basis spline solutions improve the standard finite element method. The simulation results which are verified using the values reported in the literature, point out to better efficiency in terms of the accuracy level.     

The bidimensional interproximation problem and mixed splines

Interproximation methods for surfaces can be used to construct a smooth surface interpolating some data points and passing through specified regions. In this paper we study the use of mixed splines, that is smoothing splines with additional interpolation constraints, to solve the interproximation problem for surfaces in the case of scattered data. The solution is obtained by solving a linear system whose structure can be improved by using “bell-shaped” thin plate splines.     

An efficient algorithm for periodic Hermite spline interpolation with shifted nodes

Generalized Hermite spline interpolation with periodic splines of defect 2 on an equidistant lattice is considered. Then the classic periodic Hermite spline interpolation with shifted interpolation nodes is obtained as a special case.By means of a new generalization of Euler-Frobenius polynomials the symbol of the considered interpolation problem is defined. Using this symbol, a simple representation of the fundamental splines can be given. Furthermore, an efficient algorithm for the computation of the Hermite spline interpolant is obtained, which is mainly based on the fast Fourier transform.     

On the Excitation of Rolling Tires from the Road Surface Texture

Rolling tires are excited from the contact with the rough road surface to vibrations, which cause rolling noise. A two scale approach is suggested, where at the macro–scale the vibration of the rolling tire structure is modeled by quite detailed finite element methods. The road surface is described using measured textures. A fine resolution finite element discretization of the tread rubber is performed in order to resolve the asperity contact. The material properties are described by a non–linear viscoelastic rubber model. The tread patch is enforced to approach the rough surface in a transient dynamics manner. From these investigations an enveloping surface profile is reconstructed to be used for the excitation. (© 2008 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

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.     

Thermo-mechanical modeling of turbine blades taking into account structural defects

In turbine blades of aero-engines typical defects are cracks due to high mechanical and thermal loads. The extended finite element method (XFEM) is used for simulations of fracture mechanics problems with cracks. Discontinuities in the displacement and temperature field are allowed and the crack opening displacement and crack tip stress field are reproduced accurately. Since crack closure and non-physical penetration of the crack surfaces may occur under certain load conditions, it becomes necessary to enforce the non-penetration condition for crack surfaces. This contact formulation is assumed to be frictionless. The node-to-segment approach proposed in [3] is extended to ten-node tetrahedral elements with quadratic shape functions. (© 2017 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

On an implicit particle method for simulation of forming processes

Task is the simulation of forming processes using particle methods. We implemented some mesh-free methods (the element free Galerkin method [1] and others) and the finite element method in one programme system which permits a direct comparison. For the mesh-free methods a moving least squares approximation is applied. The shape functions are not zero or one at the nodes, thus essential boundary conditions cannot be imposed directly [2]. We use a penalty method to enforce essential boundary conditions and contact conditions. The contact algorithm (normal contact of nodes to C¹-continuous surfaces) is checked by means of the element free Galerkin method and the FEM on the basis of numerical examples which deal with forming processes. (© 2005 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

A new finite element of C cubic splines

In this paper, we give a new finite element with dimension 16 of C¹ cubic splines which have interpolation schemes on the Morgan-Scott construction ₀ of a triangle.     

Compactly supported fundamental functions for spline interpolation

Summary In this paper various ways of constructing locally supported fundamental splines leading to highly accurate local interpolation schemes are proposed and analyzed.This work was partially supported by NATO grant     

A double-layer interpolation method for implementation of BEM analysis of problems in potential theory

A double-layer interpolation method (DLIM) is proposed to improve the performance of the boundary element method (BEM). In the DLIM, the nodes of an element are sorted into two groups: (i) nodes inside the element, called source nodes, and (ii) nodes on the vertices and edges of the element, called virtual nodes. With only source nodes, the element becomes a conventional discontinuous element. Taking into account both source and virtual nodes, the element becomes a standard continuous element. The physical variables are interpolated by continuous elements (first-layer interpolation), while the boundary integral equations are collocated at the source nodes only. We further established additional constraint equations between source and virtual nodes using a moving least-squares (MLS) approximation (second-layer interpolation). Using these constraints, a square coefficient matrix of the overall system of linear equations was finally achieved. The DLIM keeps the main advantages of MLS, such as significantly alleviating the meshing task, while providing much better accuracy than the traditional BEM. The method has been used successfully for solving potential problems in two dimensions. Several numerical examples in comparison with other methods have demonstrated the accuracy and efficiency of our method.     

A $\theta$-$L$ Approach for Solving Solid-State Dewetting Problems

We propose a $\theta$-$L$ approach for solving a sharp-interface model about simulating solid-state dewetting of thin films with isotropic/weakly anisotropic surface energies. The sharp-interface model is governed by surface diffusion and contact line migration. For solving the model, traditional numerical methods usually suffer from the severe stability constraint and/or the mesh distribution trouble. In the $\theta$-$L$ approach, we introduce a useful tangential velocity along the evolving interface and utilize a new set of variables (i.e., the tangential angle $\theta$ and the total length $L$ of the interface curve), so that it not only could reduce the stiffness resulted from the surface tension, but also could ensure the mesh equidistribution property during the evolution. Furthermore, it can achieve second-order accuracy when implemented by a semi-implicit linear finite element method. Numerical results are reported to demonstrate that the proposed $\theta$-$L$ approach is efficient and accurate.     

A Two-Dimensional Computational Droplet Contact Model

The interaction between capillary fluid films and micro-structural rough surfaces is one of the main challenges in studying self-cleaning mechanisms. The surface behavior of the deformable fluid film is governed by the Young-Laplace equation, which is highly non-linear. Therefore, a numerical solution is introduced using the finite element method, based on a continuum mechanical formulation. Surface and line contact at the fluid-structure interface are modeled by enforcing a contact constraint, and a contact angle, respectively. The numerical solution is validated against the analytical solution of a test case. (© 2011 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)