The first-kind and the second-kind boundary integral equation systems for solution of frictionless contact problems

An original approach to the numerical solution of displacement boundary integral equation (BIE) and traction hypersingular boundary integral equation (HBIE) by the boundary element method (BEM) for contact problems is given. The main point is to show, how the contact conditions are used to formulate the first-kind and the second-kind BIE systems in the case of frictionless two-body elastic contact. The solution of the first-kind BIE is performed by symmetric Galerkin BEM; the second-kind BIE is solved by an appropriate collocation BEM. The contact problem in itself is solved by the method of subsequent approximations of contact region. Both forms of BIE system are compared in several numerical examples. This comparison is made for different kinds of contact problem. The major emphasis is put on the evaluation of contact pressure. The obtained results are compared with referenced numerical and with the analytical ones.     

Effects of phase lags on wave propagation in an infinite solid due to a continuous line heat source

The present work is concerned with the wave propagation in a homogeneous, isotropic and unbounded solid due to a continuous line heat source under the theory of thermoelasticity with three phase-lags (Roychoudhari in J Therm Stress 30:231?C238, 2007). For the solution of the problem, we employ a potential function approach together with Laplace and Hankel transform method. Analytical expressions for the distributions of different fields like temperature, displacement and stresses inside the medium are derived by inverting Laplace transforms in an approximate manner for small values of time. The problem is illustrated by computing numerical values of the field variables for a particular material. The theoretical as well as numerical results are compared with the corresponding results for other theories of thermoelasticity reported earlier.     

A Thick Plate Problem in the Theory of Generalized Thermoelastic Diffusion

In this work, the problem of a thermoelastic thick plate with a permeating substance in contact with one of the bounding planes is considered in the context of the theory of generalized thermoelastic diffusion with one relaxation time. The bounding surface of the half-space is taken to be traction free and is subjected to a time-dependent thermal shock. The chemical potential is also assumed to be a known function of time on the bounding plane. Laplace transform techniques are used. The solution is obtained in the Laplace transform domain by using a direct approach. The solution of the problem in the physical domain is obtained numerically using a numerical method for the inversion of the Laplace transform based on Fourier expansion techniques. The temperature, displacement, stress, and concentration as well as the chemical potential are obtained. Numerical computations are carried out and represented graphically.     

Generalized thermoelastic diffusion problem for an infinitely long hollow cylinder for short times

In this work, we consider the problem of a thermoelastic infinitely long hollow cylinder in the context of the theory of generalized thermoelastic diffusion with one relaxation time. The outer surface of the cylinder is taken traction free and subjected to a thermal shock, while the inner surface is taken to be in contact with a rigid surface and is thermally insulated. Laplace transform techniques are used. The solution is obtained in the Laplace transform domain by using a direct approach. The solution of the problem in the physical domain is obtained numerically using a numerical method for the inversion of the Laplace transform based on Fourier expansion techniques. The temperature, displacement, stress and concentration as well as the chemical potential are obtained. Numerical computations are carried out and represented graphically.     

Numerical solutions for unsteady flow in unconfined aquifers

Two numerical methods for solving the problem of unsteady flow in unconfined aquifers are studied. They are an explicit finite difference method (FDM), and the finite element method (FEM). The FEM is further subdivided into three schemes: vertical displacement approach, explicit scheme (FEM1), normal velocity approach, explicit scheme (FEM2), and vertical displacement approach, implicit scheme (FEM3). Results from the above methods are compared with experimental results from a sand box model. Various factors affecting the accuracy and numerical stability are investigated. Results indicate that, for a similar degree of accuracy, the FEM requires less computational effort than the explicit FDM. Amongst the three FEM schemes, FEM3 appears to be most attractive as it is the most stable and economical of the three schemes compared.     

齿轮啮合过程中安定状态残余应力的数值方法研究

针对齿轮接触表面层变曲率的特点,引入局部坐标系,构造出局部坐标下残余应力-应变场的分布状态,建立了齿轮连续啮合过程中安定状态残余应力的计算方法.该数值方法将弹塑性问题分解为弹性问题和特征应变决定的残余问题,并采用增量映射方法求解特征应变决定的残余问题,可直接得到齿轮安定状态下的接触残余应力.采用该数值方法分别计算了标准渐开线齿轮和齿顶修缘齿轮啮合过程中安定状态下残余应力的分布状态,并与有关试验结果相比较,验证了该数值方法的有效性.     

Evaluation of the light line displacement location for object shape detection

Abstract

The ability for detecting the location of the displacement of a projected light line with subpixel resolution is examined. This evaluation is carried out for measuring a three-dimensional object shape. The technique for measuring the object shape is based on the light line projection. To determine the height data, displacement of the light line is measured on a reference point and on the object surface. This line displacement is measured by approximating the line intensity to a continuous function by the least-squares method, the Gaussian approximation and the Bezier curves method. The light line position is defined by calculating a maximum from this continuous function. To evaluate these three methods it is necessary to know the best accuracy of the experimental results. To achieve this, the results extracted by these methods are compared with a contact method. This contact method is made with a coordinate-measuring machine. The rms error is calculated using data from these three methods and the data from the contact method to carry out this evaluation. Height data are obtained with great accuracy from the method, which has a better rms value. The time necessary for processing a light line using each of these three methods is also presented. In addition to this information, a correlation coefficient is presented giving the deviation of each function generated by these methods. This line projection technique has a great potential because a very simple experimental set-up is used which is inexpensive. This evaluation is tested with complicated objects and its experimental results are presented.     

An integral method for extremely low frequency magnetic shielding

A novel approach for analyzing conducting shields of extremely low frequency magnetic fields in linear media is presented. It consists of an integral formulation based on the cell method, expressed in terms of network-like loop currents and magnetic vector potential line integrals on the shield surface. This formulation leads to a considerable reduction of field problem variables, thus limiting the amount of allocated memory and speeding-up the numerical procedure compared to other differential and integral techniques. Eddy currents are computed first, then the magnetic vector potential and the magnetic flux density distributions are evaluated by applying the superimposition principle. A detailed comparison between this method and a three-dimensional finite element method code demonstrates the accuracy of the results and the advantages of the method.     

Two-dimensional problem for a half-space with axi-symmetric distribution in the theory of generalized thermoelastic diffusion

In this work, we study a two-dimensional problem of axi-symmetric distribution of temperatures in a half-space with a permeating substance in contact with the bounding plane in the context of the theory of generalized thermoelastic diffusion with one relaxation time. The surface of the half-space is taken as traction free and subjected to axi-symmetric time-dependent thermal shock. The chemical potential is also assumed to be a known function of time on the bounding plane. The Laplace and Hankel transform techniques are used. The analytical solution in the transform domain is obtained by using a direct approach. The inverse of the double transform is obtained by using a numerical method based on Fourier expansion techniques. Numerical results for the temperature, displacement, stress, concentration, and chemical potential are carried out and represented graphically.     

Optimization under adaptive error control for finite element based simulations

Optimization problems constrained by complex dynamics can lead to computationally challenging problems especially when high accuracy and efficiency are required. We present an approach to adaptively control numerical errors in optimization problems approximated using the finite element method. The discrete adjoint equation serves as a key tool to efficiently compute both parameter sensitivities and goal-oriented error estimates at the same discretized levels. By using a recovery method for the error estimators, we avoid expensive higher order adjoint calculations. We nest the adaptivity of the mesh within the optimization algorithm, which is responsible for converging both the state and optimization algorithms and thereby allowing the reuse of state, parameters, and reduced Hessian in subsequent optimization iterations. Our approach is demonstrated on a parameter estimation problem for contamination transport in a contact tank reactor. Significant efficiency and accuracy improvements are realized in comparison to uniform grid refinement strategies and black-box optimization methods. A flexible and maintainable software interface was developed to provide access between the underlying linear algebra of a production simulator and advanced numerical algorithms such as optimization and error estimation.     

考虑异质材料的线接触性能建模与分析

提出了一种考虑材料非均匀特性的线接触弹性场数值建模方法。该方法基于数值化等效夹杂方法（NEIM）,以矩形异质材料作为基本单元,并在此基础上,结合镜像法,可求解得到线接触载荷作用下的材料弹性场。利用对比验证、参数化分析及算例验证的方法,将本文方法计算结果与有限元法（FEM）和传统等效夹杂方法（EIM）计算结果进行对比,并分析研究了异质材料分布参数对材料线接触性能的影响。结果表明:本文方法相较FEM和传统EIM具有优越性,可处理和简化涂层问题及任意分布的异质材料线接触问题。不同位置异质材料将引起基体最大von Mises应力的增大或减小;异质材料的尺寸以及异质材料和基体之间的材料差异将影响基体应力集中程度。      

铆接界面端三维奇异性应力场的研究

针对铆接结构的特点,应用特征函数扩展技术分析柱坐标下接触界面端的应力奇异性问题。建立了柱坐标下圆柱体端面接触边缘附近的三维渐近位移场和应力场渐近表达式,并根据铆钉/被铆接件接触界面端的位移和应力边界条件,建立一个非线性特征方程组。据此方程组可求解界面端邻域的应力奇异性指数、位移和应力角分布函数的数值解。通过与有限元方法计算结果相对比,验证了该方法的有效性。分析了平头、沉头以及半圆头铆钉构成的铆接结构的应力奇异性问题,考察了铆钉材料、几何形式和摩擦系数对接触界面端应力奇异性指数和应力场角分布的影响。     

Free vibration analysis of thin arbitrarily laminated anisotropic plates using boundary-continuous displacement Fourier approach

A boundary continuous displacement based Fourier series solution to the boundary-value problem of free vibration of an arbitrarily laminated thin rectangular plate is presented. This powerful approach is employed to solve a system of three highly coupled partial differential equations arising from the Kirchhoff hypothesis as applied to an anisotropic laminate, with the SS2-type simply supported boundary conditions prescribed at all four edges. The accuracy of the computed eigenvalues (natural frequencies) is ascertained by studying the convergence characteristics of the lowest seven natural frequencies, and also by comparison with the computed degenerate FEM (finite element methods) results. Other important numerical results presented include variation of the response quantities of interest with geometric and material parameters, such as fiber orientation angle and longitudinal-to-transverse modulus ratio.     

Hybrid‐mixed discretization of elasto‐dynamic contact problems using consistent singular mass matrices

An alternative spatial semi‐discretization of dynamic contact based on a modified Hamilton's principle is proposed. The modified Hamilton's principle uses displacement, velocity and momentum as variables, which allows their independent spatial discretization. Along with a local static condensation for velocity and momentum, it leads to an approach with a hybrid‐mixed consistent mass matrix. An attractive feature of such a formulation is the possibility to construct hybrid singular mass matrices with zero components at those nodes where contact is collocated. This improves numerical stability of the semi‐discrete problem: the differential index of the underlying differential‐algebraic system is reduced from 3 to 1, and spurious oscillations in the contact pressure, which are commonly reported for formulations with Lagrange multipliers, are significantly reduced. Results of numerical experiments for truss and Timoshenko beam elements are discussed. In addition, the properties of the novel discretization scheme for an unconstrained dynamic problem are assessed by a dispersion analysis.Copyright © 2013 John Wiley & Sons, Ltd.     

A contact analysis approach based on linear complementarity formulation using smoothed finite element methods

Based on the subdomain parametric variational principle (SPVP), a contact analysis approach is formulated in the incremental form for 2D solid mechanics problems discretized using only triangular elements. The present approach is implemented for the newly developed node-based smoothed finite element method (NS-FEM), the edge-based smoothed finite element method (ES-FEM) as well as standard FEM models. In the approach, the contact interface equations are discretized by contact point-pairs using a modified Coulomb frictional contact model. For strictly imposing the contact constraints, the global discretized system equations are transformed into a standard linear complementarity problem (LCP), which can be readily solved using the Lemke method. This approach can simulate different contact behaviors including bonding/debonding, contacting/departing, and sticking/slipping. An intensive numerical study is conducted to investigate the effects of various parameters and validate the proposed method. The numerical results have demonstrated the validity and efficiency of the present contact analysis approach as well as the good performance of the ES-FEM method, which provides solutions of about 10 times better accuracy and higher convergence rate than the FEM and NS-FEM methods. The results also indicate that the NS-FEM provides upper-bound solutions in energy norm, relative to the fact that FEM provides lower-bound solutions.     

Fracture analysis of cracks in magneto-electro-elastic solids by the MLPG

A meshless method based on the local Petrov–Galerkin approach is proposed for crack analysis in two-dimensional (2-D) and three-dimensional (3-D) axisymmetric magneto-electric-elastic solids with continuously varying material properties. Axial symmetry of geometry and boundary conditions reduces the original 3-D boundary value problem into a 2-D problem in axial cross section. Stationary and transient dynamic problems are considered in this paper. The local weak formulation is employed on circular subdomains where surrounding nodes randomly spread over the analyzed domain. The test functions are taken as unit step functions in derivation of the local integral equations (LIEs). The moving least-squares (MLS) method is adopted for the approximation of the physical quantities in the LIEs. The accuracy of the present method for computing the stress intensity factors (SIF), electrical displacement intensity factors (EDIF) and magnetic induction intensity factors (MIIF) are discussed by comparison with numerical solutions for homogeneous materials.     

A piezoelectric screw dislocation interacting with a nonuniformly coated circular inclusion

This paper investigated the interaction of a piezoelectric screw dislocation with a nonuniformly coated circular inclusion in an unbounded piezoelectric matrix subjected to remote antiplane shear and electric fields. In addition to having a discontinuous displacement and a discontinuous electric potential across the slip plane, the dislocation is subjected to a line force and a line charge at the core. The alternating technique in conjunction with the method of analytical continuation is applied to derive the general solutions in an explicit series form. This approach has a clear advantage in deriving the solution to the heterogeneous problem in terms of the solution for the corresponding homogeneous problem. The presented series solutions have rapid convergence which is guaranteed numerically. The image force acting on the piezoelectric screw dislocation is calculated by using the generalized Peach–Koehler formula. The numerical results show that the varying thickness of the interphase layer will exert a significant influence on the shear stress and electric field within the circular inclusion, and on the direction and magnitude of the image force. This solution can be used as Green’s function for the analysis of the corresponding piezoelectric matrix cracking problem.     

On contact zone models for an electrically impermeable interface crack in a piezoelectric bimaterial

An electrically impermeable interface crack between two semi-infinite piezoelectric planes under remote mechanical tension-shear and electrical loading is studied. Assuming the stresses, strains and displacements are independent on the coordinate x ₂ the expressions for the elastic displacement and potential jumps as well as for the stresses and electrical displacement along the interface via a sectionally holomorphic vector function are found. Introducing an artificial contact zone at the right crack tip and assuming the materials possess the symmetry class 6 mm the problem is reduced for a wide range of bimaterial compounds to a combination of combined Dirichlet–Riemann and Hilbert boundary value problems which are solved analytically. From these solutions clear analytical expressions for characteristic mechanical and electrical parameters are derived. As particular cases of the above mentioned solution the classical (oscillating) and contact zone solutions are obtained. Further, a comparison with an associated solution for an electrically permeable crack has been performed. The fracture mechanical parameters for all models via the remote loads are found analytically and important relationships between these parameters are obtained. Due to these relationships an important algorithm of a numerical method applicable for the investigation of an interface crack in a finite sized piezoelectric bimaterial is suggested.     

Effect of numerical artificial corners appearing when using BEM on contact stresses

It has long been established that when a small-displacement problem involving contact between two or more bodies is solved, conforming algorithm approaches based on BEM are reliable and robust. Nevertheless, the authors have shown in a previous paper that a relative sliding between the bodies in contact, even if small, can originate a fictitious corner in the numerical model, which may alter the stress and the displacement predictions. This problem disappears if a non-conforming algorithm is employed. In this paper a new problem with different features but also presenting artificial alteration of the results is addressed. The compression of a rounded punch on a foundation is analysed for many fillet radii. Conforming and non-conforming algorithms are used. Comparison of results shows that effects of an artificial ‘numerical corner’ may alter even the traction distribution relatively far from it. Consequently, the assumed connection between small displacements and initially conforming discretizations with the use of a conforming approach disappears. An adequate knowledge of the contact problem under consideration is then essential to determine whether the singularity is inherent to the problem or is artificial, a non-conforming algorithm then being recommended in this second case.     

The mixed finite element method applied to two-dimensional elastic contact problems

The application of the mixed finite element method to two-dimensional elastic contact problems is investigated. Since in the mixed method, both displacements and stresses are retained as variables, it is found that all the contact conditions—displacement as well as stress—can be approximated directly. A significant novelty is that some of the displacement variables are treated as natural boundary conditions in the contact region. In cases where the contact region is independent of the applied loading, an iterative procedure is used to establish the sliding and adhering portions of the contact region. In cases where the contact region is a function of the applied loading, for example progressive contact, an incremental formulation is employed to describe the discretized contact stages before invoking the former iterations. Several numerical examples are presented and the results are compared with those from the conventional potential energy or displacement finite element method.