Next Article in Journal
A Macroeconomic SIR Model for COVID-19
Previous Article in Journal
A Remark on the Change of Variable Theorem for the Riemann Integral
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Mesh Free Radial Point Interpolation Based Displacement Recovery Techniques for Elastic Finite Element Analysis

1
Civil Engineering Department, College of Engineering, King Khalid University, Abha 61421, Saudi Arabia
2
State Insurance Corporation, New Delhi 110077, India
3
Mechanical Engineering Department, Royal Commission of Jubail & Yanbu, Yanbu 41912, Saudi Arabia
*
Author to whom correspondence should be addressed.
Mathematics 2021, 9(16), 1900; https://doi.org/10.3390/math9161900
Submission received: 9 July 2021 / Revised: 2 August 2021 / Accepted: 3 August 2021 / Published: 10 August 2021

Abstract

:
The study develops the displacement error recovery method in a mesh free environment for the finite element solution employing the radial point interpolation (RPI) technique. The RPI technique uses the radial basis functions (RBF), along with polynomials basis functions to interpolate the displacement fields in a node patch and recovers the error in displacement field. The global and local errors are quantified in both energy and L2 norms from the post-processed displacement field. The RPI technique considers multi-quadrics/gaussian/thin plate splint RBF in combination with linear basis function for displacement error recovery analysis. The elastic plate examples are analyzed to demonstrate the error convergence and effectivity of the RPI displacement recovery procedures employing mesh free and mesh dependent patches. The performance of a RPI-based error estimators is also compared with the mesh dependent least square based error estimator. The triangular and quadrilateral elements are used for the discretization of plates domains. It is verified that RBF with their shape parameters, choice of elements, and errors norms influence considerably on the RPI-based displacement error recovery of finite element solution. The numerical results show that the mesh free RPI-based displacement recovery technique is more effective and achieve target accuracy in adaptive analysis with the smaller number of elements as compared to mesh dependent RPI and mesh dependent least square. It is also concluded that proposed mesh free recovery technique may prove to be most suitable for error recovery and adaptive analysis of problems dealing with large domain changes and domain discontinuities.

1. Introduction

The finite element method (FEM) has become a widely accepted method for the solution of solid mechanics problems. However, the error in finite element solutions is introduced by the very process of discretization of the problem domain. The causes of errors include error due to approximation made during mathematical modeling of the physical situation, error due to rounding off, and error due to discretization. The discretization error connected with the finite element mesh can be controlled with the choice of the type and size of elements. The discretization errors are introduced when the displacement interpolation polynomial does not accurately represent the behavior of the continuum. The current research developments in the finite element method include the improvement of the accuracy, effectivity, and reliability of FEM results of industrial problems [1]. The recent developments in high-performance finite element methods are reviewed by Cen et al. [2]. In the recovery-based error estimation, the difference between the values of recovered more accurate displacements/stress and displacements/stress obtained by the finite element solution provides a measure of the local error. The point-wise definitions are not only difficult to estimate and become meaningless at singularities. Therefore, alternative scale measures or norms are preferable. A particularly useful one is given by the energy norm, which represents the error in the rate of energy dissipation. The error norms of the displacements, such as the L2 norm, representing straightforward physical meaning, show for many finite elements super-convergence properties [3]. Zienkiewicz and Zhu [4] presented a local projection technique, a super convergent patch recovery technique, to recover the errors in stresses using the least square fit of the local polynomial to the super convergent value of the derivatives over node patch. The post-processing technique to recover the stresses and displacements from the finite element solution is presented by Niu and Shephard [5]. The error post-processing technique for displacements in the L2 norm is due to Li and Wiberg [6]. The recovery by equilibrium and compatibility in patches is proposed by Ubertini [7] to recover the error in the FEM solution for stresses. The traditional recovery techniques use the least square fitting of displacement/stress over the mesh dependent patches. The mesh free methods are recently introduced for recovery of FEM solutions errors through various mesh free support domains based interpolating techniques. The mesh free approaches are developed for situations in which distortion of elements occurs due to large domain changes and domain discontinuities. A coupled element-free Galerkin method (EFG) Method and finite element method (FEM) error recovery procedure based on local maximum entropy shape functions, for material and geometrical nonlinearities problems, is proposed by Ullah et al. [8]. The Chung and Belytschko [9] error estimations are used in the EFG regions and Zienkiewicz and Zhu (ZZ) super convergent patch recovery is used for strains and stresses in the FE region of the problem domain.
Nadal et al. [10] propose the explicit-type recovery error estimator in energy norm for the linear elasticity problem using a smooth solution. The discontinuous Galerkin error estimator-based hp-adaptive finite element analyses were carried out by Bird et al. [11]. The performances of the hp-adaptive scheme are compared with uniform and h-refinement adaptive scheme achieving the exponential rates of convergence. Liu et al. [12] propose a linear conforming point interpolation methodology-based meshfree Galerkin method for two-dimensional solid mechanics problems. Liu and Zhang [13] have proposed a strain-constructed point interpolation method (SC-PIM) for the analysis of static, free, and forced vibration structures problems. A thorough analysis of error estimation considering moving least square interpolation is presented by Mirzaei [14]. The study of factors affecting the accuracy and reliability of smoothed radial point interpolation method (RPIM) is carried out by Hamrani et al. [15]. A hybrid approach based on smoothed radial point interpolation method (CS-RPIM) and FEM is proposed by Zhang et al. [16] for solving fracture problems. Cao et al. [17] proposed a hybrid mesh free Galerkin RPIM and FEM to impose the essential boundary conditions. The mesh free recovery method using support domains, i.e., zone of influence of a node, was recently introduced by Ahmed et al. [18] to recover the FEM solution errors. The comparison of error estimation behavior using meshfree methods, RPI error recovery method, and MLS interpolation, and mesh dependent ZZ recovery method was presented by Ahmed [19] and concluded that the quality of error estimation of meshfree recovery methods was better than the mesh dependent error recovery methods. A two-dimensional finite element model employing mesh free nodes interface model based radial point interpolation method to simulate the interaction between soil and structure is proposed by Gong et al. [20].
From the literature review, it is evident that most of the techniques proposed to recover the discretization error introduced in the finite element solution are mesh-based, i.e., dependent on element mesh connectivity, and the mesh free recovery technique in finite element method is a recent interest. The existing mesh-based techniques may not perform well or are difficult to implement in situations to deal with large domain changes and domain discontinuities. Moreover, more attention has been paid to recover errors in the field variable derivatives. Therefore, there is a need to develop mesh free recovery techniques using the different mesh free procedures with special attention to recover errors in the field variable. The present study aimed to investigate the reliability and effectiveness of RPIM-based displacement error recovery techniques in the mesh free environment. The comparison of mesh free RPIM-based displacement error recovery techniques with mesh dependent RPIM and least square based recovery techniques is also carried out. The study also compares the performance of RPI-based error recovery techniques with Zienkiewicz-Zhu (ZZ) recovery techniques presented in the literature [18]. The benchmark elastic plates examples are analyzed to compare the rate of error convergence, effectivity, and adaptively improved meshes with different displacement recovery procedures. The linear triangular, quadratic triangular, and quadrilateral elements are used for problems domain discretization. The errors of finite element solution are computed directly from the recovered displacement in terms of energy and L2 norm. The influence of radial basis functions in RPI-based recovery techniques is also assessed by considering three radial basis functions, namely multi-quadrics, gaussian (exponential), and thin plate splint for the RPI technique. The effect of radial basis functions with their shape parameters and shape of node zones is also assessed on the quality of the error estimation obtained from the RPI-based error recovery techniques.

2. Elastic Problem Statement

The elastic problems are considered with stress field, “σ” and unknown displacement field “u”, defined over a domain Ω that is bounded by Γ = Γt ∪ Γu. The prescribed tractions ( t ¯ ) and displacements ( u ¯ ) are imposed over boundary Γt and Γu, respectively.
The equilibrium conditions and boundary conditions are satisfied to find (σ,u).
LT σ + f = 0   in Ω,
where f is the force vector, LT is the derivative operator.
σ · n   =   t ¯ on   Γ t ,
u   =   u ¯ on   Γ u ,
Constitutive relation:
σ = D ε, with ε = L u,
where D is the elasticity matrix of linear isotropic material and ε is the strain vector.

3. Least Square Interpolation Technique for Displacement Recovery

The recovery of field variable (displacement) is obtained by least squares fit of the computed nodal field variable using a higher order polynomial over an element neighborhood patch consists of a union of the elements surrounding an element. To perform least square fitting, the following functional is minimized.
πf (a) = 1/2∑i=1np [dih(xi, yi) − di (xi, yi)]2
where ui (xi, yi) = Pi (xi, yi) . a
ui = [ui vi ]T , a = [au av ]T
where ui and vi are the nodal parameters of field variables in x and y direction and a is the vector of unknown parameters au and av.
pi = [1, xi, yi, xi2, xiyi, yi2, . . . . . . .],
where (xi, yi) are the sampling points (np) coordinates.
Minimization condition of πf (a) implies that “a” satisfies the following relation.
i=1np Pi T(xi, yi). Pi (xi, yi). a = ∑i=1np Pi T(xi, yi). di h(xi, yi)
Solving for “a”, the following relation is obtained.
a = A−1 b,
where
A = ∑i=1np Pi T(xi, yi). Pi (xi, yi).a, and b = ∑i=1np Pi T(xi, yi). di h(xi, yi)

4. Radial Point Interpolation Method (RPIM) for Displacement Recovery

The radial point interpolation method (RPIM) considers both the radial basis functions (RBF) and polynomials basis function for interpolation [21]. The approximation uh(x) for the field variables u(x) at a point xq, assuming only the influencing nodes of a point xq have an effect on u(x), is given by,
uh (x, xq) = ∑i=1n Ri ai + ∑j=1m pj bi = RT(x) a + pT(x) b,
where Ri is the RBF computed at xi, the ith node in the support domain of xq, pj is the monomial in the polynomial basis, ai and bj are the coefficient for the radial basis function (Ri) and polynomial basis (pj) respectively, n is the number of nodes in the mesh free nodes zones of xq, m is the number of monomials in the polynomial basis.
a = [a1, a2,…,……..an]
b = [b1, b2,…,……..bn]
R(x) = [R1(x), R2(x),……,…… Rn(x)]T,
p(x) = [p1(x), p2(x),……,…… pn(x)]T,
the text following an equation need not be a new paragraph. Please punctuate equations
The number of the radial basis ‘n’ and the order of polynomial basis ‘m’ is chosen based on the reproduction requirement. Minimum terms of polynomial basis are often adopted for better stability. The coefficients ai, and bj are determined in this way such that Equation (12) passes through ‘n’ data nodes in the mesh free nodes zones and imposing constraints on the polynomial basis functions to guarantee a unique solution [22]. Thus, two intermediate matrices Sa and Sb result as follows,
Rb = [PmT Rq−1 Ri]−1 PmT Rq−1,
Ra = PmTRq−1 Pm Sb,
where Rq is the moment matrix associated to the radial basis function and Pm is the moment matrix related to the polynomial terms. The solution of Equation (12) is only possible if the number of unknown parameters a is smaller than, or at the most equal to, the number of independent equations.
The ‘n’ number of equations can be obtained using Equations (11), (16) and (17) as,
uh = Rq a + Pm b,
The following constraints are imposed to obtain unique solutions of the equation,
pmT a = 0,
Combining Equations (18) and (19) yields.
a = Sa ue,
b = Sb ue,
Substituting the Equations (20) and (21) into Equation (18), the interpolation function ψ(x) can be obtained as.
uh(x, xq) = (RT SapT Sb) ue = ψ(x) ue,

4.1. Multi-Quadrics Radial Function (MQ)

The radial basis function performs an influential role on the behavior of RPI methods. The various types of radial basis functions can be used for RPI technique. The formulation of multi quadrics (MQ) radial basis function with dimensionless shape parameters is as follows.
Ri (xj) = [(xj − xi)2 + (yj − yi)2 + (α0_dc)2]q, α0 ≥ 0
where α0 is a dimensionless coefficient and dc represents characteristic length that relates to the nodal spacing in meshfree nodes zones or nodes patch of the point of interest x, which is the shortest distance between the node i and neighborhood nodes.

4.2. Gaussian (Exponential) Radial Basis Function (Exp)

The formulation of Gaussian (Exponential) radial basis function with dimensionless shape parameters is as follows.
Ri (xj) = exp {−α0 [(xjxi )2 + (yjyi)2/dc2]},

4.3. Thin Plate Splint (TPS) Radial Basis Function (Exp)

The formulation of thin plate splint (TPS) radial basis function is as follows.
Ri (xj) = {√[(xjxi)2 + (yjyi)2]}q,
The α0, q and η are the shape parameters of the RBFs. In utilizing RBFs, the shape parameters need to be determined for good performance. Effects of α0 and q on radial basis functions have been widely investigated [23]. It was found that for 2D square plate problem using MQ radial basis function and errors quantified in energy norm, α0 = 5.0 and q = 1.03 lead to good results while α0 should be less 1 for Gaussian RBF. Hence, these values were used in the study. However, for L2 norm, α0 = 1 and q should be less than 0.5 for best performance with MQ. It was also observed that with increasing node zone size, the value of α0 should also be increased for the best performance of RPI-based recovery analysis. The value of α0 = 0.65 and q = 1.03 is taken in plate problem with opening for MQ RBF. The value of α0 = 0.65 is taken with Gaussian RBF and η = 5 is taken with TPS RBF for error quantified in the energy norm.
The second term in Equation (11) consists of polynomials. The linear polynomial added into the RBF can ensure linear consistency and improve interpolation accuracy [23]. The linear polynomial basis with m as 3 is used in the RPI technique, which is given as
pT(xi) = [1, xi, yi,],
The accuracy of interpolation for the point of interest depends on the nodes in the mesh free nodes zones. The circular form of mesh free nodes zones is constructed using the distance d = ║x − xi║/dm. The (xxi) is the distance from node x to point xi and dm is the size of the influence domain of the point xi, the support size of the ith node, dmi, is computed by dmi = dmax ci, in which dmax is a scaling factor called as dilation parameter, the distance cI is determined by searching for sufficient neighbor nodes distance. The dilation parameter (dmax) is taken as 3.0 in the present study. For the construction of the rectangular form of meshfree nodes zones, the distance along two cartesian directions are rx = ║x − xi║/dmx and rx = ║y − yi║/dmx, where dmx = dmax cxi and dmy = dmax cyi. For uniformly distributed nodes, ci is simply the distance between two neighboring nodes. For nonuniformly distributed nodes, ci can be taken as an average nodal spacing in the support domain of xi.

5. Measurement of Errors

The error in computed state variable or state variable derivative i.e., displacement (u) or stress (σ), (eu* or eσ*) is defined as the difference between the exact (or numerical) values of u (or σ) and respective computed values, uh (or σh).
eu = u* − uh,
eσ = σ* − σh,
The specification of local error in the above manner is generally not convenient and may be occasionally misleading. For this reason, ‘norm’ representing integral scalar quantities are introduced to measure the error. For example, the energy norm is evaluated by computing the strain energy contained in the difference between the discontinuous value σh and the recovered values σ*. The ‘energy norm’ for the error, also called as error estimator can be written as in Equation (29).
║x − xiE = [∫Ω eσ*T D−1 eσ* dΩ]0.5,
where “D” is the elasticity matrix.
A more direct quantification measure is the so-called L2 norm, which can be associated with the error in any quantity. For example, for the displacement u, the L2 norm of the error is obtained as given in Equation (30).
║x − xiL2 = [ ∫Ω eT e* dΩ]0.5,
An estimator is asymptotically exact for a problem if the problem global and local (element) effectivity index (θ), i.e., ratio of estimated error and actual error, converges to one when the mesh size approaches to zero [24].
θ = ║e║/║eex║,
where ║eex║ and ║e║ are the actual error and the estimated error in the energy norm.
The accuracy (η) of a finite element solution may be defined as follows.
η = ║e*║/║σ*║,
║σ*2 = ║σh2 + ║e║E2,
║σh2 = ∑i=1n ║σh2i,
The solution is acceptable if ηηallow where ηallow is the allowable accuracy. If η > ηallow, modification of element size is needed.
A 2-dimensional finite element code is developed incorporating the above described mesh free/mesh dependent error recovery techniques and error estimators guided adaptive techniques. The program is run on intel core i7 with 2.6 GHz processor and 16 GB RAM to obtain the numerical results of the following plate problems analysis. The flow chart to show the adaptive finite element analysis using the error recovery techniques implementation is shown in Figure 1.

6. Elastic Plate Problems

The quality of recovery method is obtained by adaptive finite element analysis of two elastic plate problems, for which analytical solution is available [4], in terms of error convergence rate, effectivity, and adaptive meshes employing above-discussed mesh free and mesh dependent recovery techniques. The problem domains are discretized with triangular (three nodes) and quadrilateral (four nodes) mesh. The RPI recovery-based analysis considers multi-quadrics radial basis function and quadratic polynomial basis function. The errors are quantified in both energy and L2 norms. The mesh free nodes zones for RPI technique to interpolate the displacement/stress is shown in Figure 2. The least square and RPI approach uses the nodal values in the mesh dependent element patches consist of the union of the elements surrounding the element under consideration to interpolate the displacement/stress as shown in Figure 3.

6.1. Square Plate Problem

A 1 × 1 plane square plate problem was analyzed under the action of body forces (bx, by) over the domain using various recovery techniques-based FEM in an adaptive environment. The body forces in the form of polynomials are given in Equations (35) and (36). The analytical solutions for displacements (u, v) are given in Equation (37). The discretized domain of the plate using triangular and quadrilateral elements is depicted in Figure 4. The analysis results for error convergence and effectivity for displacement recovery in different norms are tabled in Table 1, Table 2, Table 3, Table 4, Table 5 and Table 6.
u = 0; v = −x y(1 − x) x (1 − y),
bx = (α + β) . (1 − 2x) . (1 − 2y),
by = −2 β y.(1 − y) − (α + 2 β) 2x .(1 − x),
The constants α and β are given as
α = E ν/[(1 − 2ν) (1 + ν)]; β = E/[2(1 + ν)]
where E, ν are Modulus of elasticity and Poisson’s ratio with a value of 1.0 N/mm2 and 0.3.

6.1.1. Mesh Free and Mesh Dependent Displacement Recovery Techniques

The finite element errors quantified in terms of L2 norms and effectivity of error estimation, using mesh free radial point interpolation (RPI) approach and mesh dependent least square (LS) and RPI approach with triangular (3/6 Nodes) and quadrilateral (4 Nodes) mesh for displacements recovery, are presented in Table 1, Table 2 and Table 3. The RPI approach considers the multi-quadratics radial basis function for interpolation purpose, and structured mesh are employed in the analysis.

6.1.2. RPI Recovery Technique and Radial Basis Function Type

The characteristics of error estimation using RPI-based recovery technique are obtained considering Multi-quadrics (MQ), Gaussian or Exponential (Exp), and Thin plate splint (TSP) radial basis functions. Table 1, Table 2 and Table 3 present Multi-quadrics (MQ) RBF-based RPI recovery analysis results wherein errors are quantified in L2 norms. The finite element errors in energy norms and effectivity of error estimation, using radial point interpolation (RPI) approach considering Multi-quadrics (MQ), Gaussian (Exponential), and Thin plate splint (TSP) RBFs, with triangular (3/6 Nodes) and quadrilateral (4 Nodes) structured mesh are tabulated in Table 4, Table 5 and Table 6.

6.1.3. RPI Recovery Technique and Patch Configuration

In almost all adaptive finite element analyses, unstructured meshes are employed. Therefore, it is of interest to test the robustness of RPI-based displacement recovery technique also in unstructured meshes. The finite element errors quantified in terms of L2 norm and effectivity of error estimation, using mesh free radial point interpolation (RPI) approach and mesh dependent RPI approach with multi-quadratic radial basis function (MQ) and least square (LS) approach considering circular and rectangular mesh free node zones of triangular (3/6 Nodes) and quadrilateral (4 Nodes) unstructured mesh for displacements recovery are presented in Table 7, Table 8 and Table 9. The computation results error estimation in energy norm with circular and rectangular mesh free node zones of quadrilateral (4 Nodes) unstructured mesh is shown in Table 10.

6.2. Square Plate with Opening Problem

A problem involving a square portion from an infinite elastic plate having a circular opening with a radius (a) of 1 unit is also analyzed to demonstrate the effectiveness of RPI-based displacement recovery technique in stress concentration conditions. The one-quarter of the plate domain is modeled because of the symmetry of the plate. The plate domain with an opening is discretized with triangular/quadrilateral elements, as shown in Figure 5. Along the symmetry line, the normal displacement component and shear stress are zero. A unit in-plane traction is applied in the x-direction. The close form solution of stress field for an infinite plate with opening is given in Equations (38)–(40) [4].
σx = σ {1 − [(a2/r2)(1.5 cos2θ + cos4θ)] + (1.5 a4/r4 cos4θ)
σy = σ {0 − [(a2/r2)(0.5 cos2θ − cos4θ)] + (1.5 a4/r4 cos4θ)
σxy = σ {0 − [(a2/r2)(0.5 sin2θ + sin4θ)] + (1.5 a4/r4 sin4θ)
where r2 = y2 + x2 and σ is the uniaxial traction applied at infinity.
The finite element analysis of the problem is carried out incorporating mesh dependent and mesh independent RPI-based displacement error recovery techniques. The error estimation obtained with RPI recovery techniques is compared with least square displacement recovery techniques. The analyses results are also obtained employing Zienkiewicz-Zhu (ZZ) recovery techniques presented in the literature [18] for comparison with RPI-based error recovery techniques. The recovery results are compared in terms of convergence errors and global effectivity. The computational results obtained are tabulated in Table 11, Table 12 and Table 13.
The efficiency and reliability of recovery technique under different mesh free and mesh dependent interpolation schemes is demonstrated through finite element analysis in an adaptive environment. The initial meshes are adaptively refined to bring the solution error below the target error limit. Table 14 shows the global error, number of elements, and degree of freedom DOF in refined meshes at prescribed error limit of 2% using the RPI error recovery employing Multi-quadrics (MQ), Gaussian (GS), and thin plate splint radial basis functions, circular and rectangular support domain shapes. The adaptively refined meshes at 2% target error using mesh free and mesh dependent RPIM, and mesh dependent least square based error estimation employing different radial basis function and node zone shapes are displayed in Figure 6, Figure 7 and Figure 8.

7. Discussion

The reliability and effectiveness of RPIM-based displacement error recovery techniques in a mesh free environment is studied. The comparison of mesh free RPIM-based displacement error recovery techniques with mesh dependent RPIM and least square based recovery techniques is also carried out. The performance of Zienkiewicz-Zhu ZZ recovery techniques [18] is also compared with RPI-based error recovery techniques. The RPI technique uses the nodal values in the local support domain to interpolate the field variable or field variable derivatives. The elastic square plate problems under the action of body forces and unit in-plane traction are analyzed using finite element method coupled with RPI-based error recovery considering mesh free circular/rectangular support domains and mesh dependent node patches. The RPIM-based displacement error recovery techniques employ Multi-Quadrics (MQ) and Gaussian or Exponential (Exp) radial basis functions. The problem domain is discretized using linear/quadratic triangular and quadrilateral shape elements. The error analysis is further utilized for adaptive improvement of domain mesh. The quality of discretization error obtained using recovery techniques is compared in terms of error convergence characteristics, error effectivity, and adaptively refined element meshes with specified target error. The errors are computed in L2 norm as well as in energy norm. The computation results obtained with increasing fineness order, radial basis functions, mesh free node zones shapes, and error norms, for elastic plate problem under the action of body forces, are depicted in Table 1, Table 2, Table 3, Table 4, Table 5, Table 6, Table 7, Table 8, Table 9 and Table 10. The results of error convergence rate, effectivity index, and adaptive refined meshes for 2% target accuracy obtained for an infinite plate with an opening are given in Table 11, Table 12, Table 13 and Table 14 and Figure 6, Figure 7 and Figure 8.
The finite element analysis results of benchmark problems using the mesh free RPI, mesh-based RPI and least square recovery techniques with increasing order of mesh fineness, show that the RPI recovery techniques predict lower order of error and higher convergence rate of the error, than that of the finite element solution and mesh dependent least square based recovery techniques with better effectivity. The mesh free RPI recovery techniques predict lower order of error with higher convergence rate than that of the mesh dependent RPI error recovery techniques. It may be due to the possibility of less accurate recovery for boundary nodes in mesh dependent recovery technique since a smaller number of nodes are available. Such difficulty does not arise with mesh free RPI-based displacement recovery technique. The effect of choice of element for domain discretization on error convergence of RPI recovery techniques and least square recovery techniques is clearly observed in the recovery analysis results. The performance of linear quadrilateral elements discretization used in mesh free and mesh dependent RPI-based recovery techniques and mesh dependent least square recovery techniques is found better than the linear and quadratic triangular element discretization as the higher convergence rate, lesser order of error and effectivity index, in general, near to one. From the error quantified in energy and L2 norms, it is found that it is better to employ the L2 norm for error representation for higher order elements as effectivity converge to nearly one and with a higher convergence rate as compared to lower order element. The convergence rate for linear triangular, quadratic triangular, and linear quadrilateral structured mesh, for finite element solution of plate problem using mesh free RPI (MQ), mesh-based RPI (MQ), and least square recovery techniques, in L2 norm are obtained as (1.93939, 2.01685, 2.06494, 2.15207), (3.04733, 4.16168, 3.86833, 3.721499), and (2.01223, 2.26171, 2.00883, 2.01729), respectively. The corresponding error convergence rate for 3/6/4 node elements in energy norm are found as (0.97875, 1.89163, 1.98282, 1.89225), (1.97848, 3.19100, 2.769632, 2.88449), and (1.00169, 2.01497, 2.05275, 2.02723), respectively. There is a pronounced effect of the kind of radial basis function on the RPI recovery-based error estimation. It is found that the characteristics of error recovery considering the Multi-Quadrics radial basis function are much superior to the Gaussian or Exponential radial basis functions in both mesh free and mesh dependent RPI recovery analysis. The convergence rate with Multi-Quadrics, Exponential, and Thin plate splint radial basis functions in mesh free RPI recovery analysis for 3/6/4 node elements in energy norm are, respectively, obtained as (1.89163, 1.48183, 1.91524), (3.19100, 3.48633, 3.09017), and (2.01497, 1.45472, 1.82359). The convergence rate with Multi-Quadrics and Exponential RBF in mesh dependent RPI recovery analysis for 3/6/4 node elements in energy norm are respectively obtained as (1.98282, 1.24389), (2.769632, 2.39122), and (2.05275, 1.44822). The shape of the support domains, i.e., circular, or rectangular, has no significant influence on the RPI-based error estimation provided that enough nodes are enclosed in the support domain. Similar results are also confirmed from the mesh free and mesh dependent RPI-based displacement error recovery analysis of an infinite plate with an opening as shown in Table 11, Table 12, Table 13 and Table 14, giving rate of convergence and effectivity index. It implies from the discussion that the type of element, the norm used to quantify the errors, and radial basis functions influence considerably on the RPI-based displacement error recovery of finite element solution.
The performance of error estimation using mesh free and mesh dependent RPI displacement error recovery technique and least square recovery technique is also demonstrated in terms of adaptively improved meshes with target error limits. The adaptive finite element results for the plate problem with an opening are presented in Table 14 and Figure 6, Figure 7 and Figure 8 at a prescribed error limit of 2% for different radial basis functions and node zones. The adaptively improved meshes indicate the error distribution pattern in the domain as the meshes become finer in areas of high errors to obtain a uniform accuracy throughout the domain. It is concluded that the RPI-based mesh free error estimation technique can be used to predict the zones of high computational errors. The number of elements required to achieve target accuracy, in general, is smaller in RPI error estimations as compared to mesh dependent RPI and mesh dependent least square error estimation for both triangular and quadrilateral elements discretization. The characteristics of adaptive mesh obtained using the mesh free RPI-based error recovery technique using Multi-Quadrics radial basis function are superior to mesh free RPI-based error recovery technique using Gaussian and thin plate splint radial basis function. It is also observed that the performance of RPI-based displacement recovery technique is superior to ZZ stress recovery technique [18]. The numerical results show that the mesh free RPI-based displacement recovery technique is more effective and achieves target accuracy in adaptive analysis with a smaller number of elements as compared to mesh dependent RPI and mesh dependent least square. It is also concluded that the proposed mesh free recovery technique may prove to be most suitable for error recovery and adaptive analysis of problems dealing with large domain changes and domain discontinuities.

8. Present Study Limitations and Future Research Work

The accuracy and the performance of the RPIM techniques depend on suitable and optimized values of shape parameters of radial basis functions. Therefore, the performance is limited to the shape parameters taken in the study. The present study is limited to the performance of mesh free RPI-based error recovery technique for two-dimensional linear elastic problems. The study could be extended the application of mesh free RPI-based error recovery to incompressible elastic problems. The study could also be conducted by applying mesh free recovery techniques to large deformation problems or plastic problems. The different adaptive procedures could also be tested using RPI-based error estimators.

9. Conclusions

The present study is a contribution towards the development of a reliable and cost-effective displacement recovery-based error estimation technique using RPIM in mesh free and mesh dependent environments. The study considers three radial basis functions, namely Multi-Quadrics, Gaussian, or Exponential and thin plate splint RBF, and circular and rectangular support domain shape for error recovery analysis. Finite element analysis on test examples is carried out using linear and higher order triangular and quadrilateral elements for discretization and employing RPI-based recovery technique. The finite element solution errors in energy/L2 norms are calculated directly from the recovered displacement. The characteristics of RPI-based error recovery with different radial basis function is compared in terms of error convergence properties, effectivity, and adaptively refined meshes. It is found from the results that the L2 norm for error representation is better for higher order elements, while the energy norm is equally good for all types of elements. It is concluded that the multi-quadrics radial basis function found to perform better in RPI-based error recovery of finite element solution and radial basis functions shape parameters can be optimized for better performance. It is verified that radial basis functions with their shape parameters, choice of elements for discretization, and norm used to quantify the errors influence considerably on the RPI-based displacement error recovery of finite element solution. The study concludes that mesh free RPI-based displacement error recovery technique with multi-quadrics radial basis function is more effective and efficient than the mesh dependent least square based displacement recovery technique. The numerical results show that the mesh free RPI-based displacement recovery technique is more effective and achieves target accuracy in adaptive analysis with a smaller number of elements as compared to mesh dependent RPI and mesh dependent least square. The study concludes that mesh free RPI-based displacement recovery technique is more effective and efficient, and the number of elements required to achieve target accuracy is smaller in RPI error estimations as compared to mesh dependent RPI and mesh dependent least square error estimation. It is also concluded that the proposed mesh free recovery technique may prove to be most suitable for error recovery and adaptive analysis of problems dealing with large domain changes and domain discontinuities.

Author Contributions

Conceptualization, M.A. and D.S.; methodology, M.A.; software, D.S.; validation, S.A. and M.A.A.; formal analysis, M.A. and S.A.; investigation, M.A.; resources, S.A.; data curation, M.A.A.; writing—original draft preparation, M.A.; writing—review and editing, S.A. and M.A.A.; visualization, S.A.; supervision, D.S.; project administration, M.A.A.; funding acquisition, M.A. All authors have read and agreed to the published version of the manuscript.

Funding

The authors extend their appreciation to the Deanship of Scientific Research at King Khalid University for funding this work through General Research Project under grant number [R.G.P2/73/41].

Institutional Review Board Statement

Not Applicable.

Informed Consent Statement

Not Applicable.

Data Availability Statement

Not Applicable.

Acknowledgments

The authors acknowledge to the Dean, Faculty of Engineering for his valuable support and help.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

References

  1. Taus, M.; Rodin, G.J.; Hughes, T.J.R.; Scott, M.A. Isogeometric boundary element methods and patch tests for linear elastic problems: Formulation, numerical integration, and applications. Comput. Methods Appl. Mech. Eng. 2019, 357, 112591. [Google Scholar] [CrossRef]
  2. Cen, S.; Wu, C.J.; Li, Z.; Shang, Y.; Li, C. Some advances in high-performance finite element methods. Eng. Comput. 2019, 36, 2811–2834. [Google Scholar] [CrossRef]
  3. Grosse, I.; Katragadda, P.; Benoit, J. An adaptive accuracy-based a posteriori error estimator. Finite Elem. Anal. Des. 1992, 12, 75–90. [Google Scholar] [CrossRef]
  4. Zienkiewicz, O.C.; Zhu, J.Z. The superconvergent patch recovery anda posteriori error estimates. Part 1: The recovery technique. Int. J. Numer. Methods Eng. 1992, 33, 1331–1364. [Google Scholar] [CrossRef]
  5. Niu, Q.; Shephard, M.S. Super-convergent Extraction Techniques for Finite Element Analysis. Int. J. Num. Meth. Eng. 1993, 36, 811–836. [Google Scholar] [CrossRef]
  6. Li, X.; Wiberg, N.-E. A posteriori error estimate by element patch post-processing, adaptive analysis in energy and L2 norms. Comput. Struct. 1994, 53, 907–919. [Google Scholar] [CrossRef]
  7. Ubertini, F. Patch recovery based on complementary energy. Int. J. Numer. Methods Eng. 2004, 59, 1501–1538. [Google Scholar] [CrossRef]
  8. Ullah, Z.; Coombs, W.; Augarde, C. An adaptive finite element/meshless coupled method based on local maximum entropy shape functions for linear and nonlinear problems. Comput. Methods Appl. Mech. Eng. 2013, 267, 111–132. [Google Scholar] [CrossRef] [Green Version]
  9. Chung, H.-J.; Belytschko, T. An error estimate in the EFG method. Comput. Mech. 1998, 21, 91–100. [Google Scholar] [CrossRef]
  10. Yang, C.S.; Lee, F.B.; Kao, S.-P.; Hung, P.-S. Twelve different interpolation methods: A case study of Surfer 8.0. In Proceedings of the XXth ISPRS Congress, Istanbul, Turkey, 12–23 July 2004; pp. 778–785. [Google Scholar]
  11. Chen, X.L.; Liu, G.R.; Lim, S.P. The effects of the enforcement of compatibility in the radial point interpolation method for analyzing mindlin plates. In Advances in Meshfree and X-FEM Methods; World Scientific: Singapore, 2002; pp. 84–89. [Google Scholar]
  12. Liu, G.R.; Zhang, G.Y.; Dai, K.Y.; Wang, Y.Y.; Zhong, Z.H.; Li, G.Y.; Han, X. A linearly conforming point interpolation method (lc-pim) for 2d solid mechanics problems. Int. J. Comput. Methods 2005, 2, 645–665. [Google Scholar] [CrossRef] [Green Version]
  13. Liu, G.R.; Zhang, G.Y. A novel scheme of strain-constructed point interpolation method for static and dynamic mechanics problems. Int. J. Appl. Mech. 2009, 1, 233–258. [Google Scholar] [CrossRef] [Green Version]
  14. Mirzaei, D. Analysis of moving least squares approximation revisited. J. Comput. Appl. Math. 2015, 282, 237–250. [Google Scholar] [CrossRef]
  15. Hamrani, A.; Belaidi, I.; Monteiro, E.; Lorong, P. On the Factors Affecting the Accuracy and Robustness of Smoothed-Radial Point Interpolation Method. Adv. Appl. Math. Mech. 2016, 9, 43–72. [Google Scholar] [CrossRef] [Green Version]
  16. Zhang, G.; Wang, Y.; Jiang, Y.; Jiang, Y.; Zong, Z. A combination of singular cell-based smoothed radial point inter-polation method and FEM in solving fracture problem. Int. J. Comput. Methods 2018, 15, 1850079. [Google Scholar] [CrossRef]
  17. Ahmed, M.; Singh, D.; Desmukh, M.N. Interpolation type stress recovery technique based error estimator for elasticity problems. Mechanika 2018, 24, 672–679. [Google Scholar] [CrossRef] [Green Version]
  18. Ahmed, M.; El Ouni, M.H.; Singh, D.; Kahla, N.B. A parametric study of meshfree interpolation based recovery techniques in finite element elastic analysis. CMES Comput. Model. Engg. Sci. 2019, 121, 687–786. [Google Scholar]
  19. Ahmed, M. A Comparative study of mesh-free radial point interpolation method and moving least squares method-based error estimation in elastic finite element analysis. Arab. J. Sci. Eng. 2020, 45, 3541–3557. [Google Scholar] [CrossRef]
  20. Gong, J.; Zou, D.; Kong, X.; Qu, Y.; Zhou, Y. A Non-Matching Nodes Interface Model with Radial Interpolation Function for Simulating 2D Soil–Structure Interface Behaviors. Int. J. Comput. Methods 2020, 18, 2050023. [Google Scholar] [CrossRef]
  21. Liu, G.R.; Gu, Y.T. A local radial point interpolation method (LR-PIM) for free vibration analyses of 2-D solids. J. Sound Vib. 2001, 246, 29–46. [Google Scholar] [CrossRef] [Green Version]
  22. Wang, J.; Liu, G. On the optimal shape parameters of radial basis functions used for 2-D meshless methods. Comput. Methods Appl. Mech. Eng. 2002, 191, 2611–2630. [Google Scholar] [CrossRef]
  23. Wang, J.; Liu, G.R. A point interpolation meshless method based on radial basis functions. Int. J. Numer. Methods Eng. 2002, 54, 1623–1648. [Google Scholar] [CrossRef]
  24. Zienkiewicz, O.C.; Zhu, J.Z. A simple error estimator and adaptive procedure for practical engineering analysis. Int. J. Numer. Methods Eng. 1987, 24, 337–357. [Google Scholar] [CrossRef]
Figure 1. Least square and RPI recovery techniques coupled adaptive elastic finite element analysis flow chart.
Figure 1. Least square and RPI recovery techniques coupled adaptive elastic finite element analysis flow chart.
Mathematics 09 01900 g001
Figure 2. Mesh free node zones (Circular/Rectangular) for RPI bases recovery technique.
Figure 2. Mesh free node zones (Circular/Rectangular) for RPI bases recovery technique.
Mathematics 09 01900 g002
Figure 3. Element patches for least square (Displacement) and RPI (Displacement) recovery technique (mesh dependent).
Figure 3. Element patches for least square (Displacement) and RPI (Displacement) recovery technique (mesh dependent).
Mathematics 09 01900 g003
Figure 4. Subdivided elements mesh (structured/unstructured) for square plate domain.
Figure 4. Subdivided elements mesh (structured/unstructured) for square plate domain.
Mathematics 09 01900 g004
Figure 5. Triangular/quadrilateral elements mesh for plate with opening domain.
Figure 5. Triangular/quadrilateral elements mesh for plate with opening domain.
Mathematics 09 01900 g005
Figure 6. Adaptively improved mesh in elastic plate with opening using mesh free (MF) RPI having different node zone shape and least square error recovery analysis (initial triangular mesh elements = 548, 2% target error).
Figure 6. Adaptively improved mesh in elastic plate with opening using mesh free (MF) RPI having different node zone shape and least square error recovery analysis (initial triangular mesh elements = 548, 2% target error).
Mathematics 09 01900 g006aMathematics 09 01900 g006b
Figure 7. Adaptively improved mesh in elastic plate with opening using meshfree (MF) RPI, mesh dependent (MD) RPI having different radial basis function (initial triangular mesh elements = 548, 2% target error).
Figure 7. Adaptively improved mesh in elastic plate with opening using meshfree (MF) RPI, mesh dependent (MD) RPI having different radial basis function (initial triangular mesh elements = 548, 2% target error).
Mathematics 09 01900 g007
Figure 8. Adaptively improved mesh in elastic plate with opening using mesh free RPI, mesh dependent RPI and least square error recovery analysis (initial quadrilateral mesh elements = 179, 2% target error).
Figure 8. Adaptively improved mesh in elastic plate with opening using mesh free RPI, mesh dependent RPI and least square error recovery analysis (initial quadrilateral mesh elements = 179, 2% target error).
Mathematics 09 01900 g008
Table 1. Error convergence and global effectivity (θ) obtained for RPI (MQ) and least square based displacement recovery techniques (linear triangular, L2 norm).
Table 1. Error convergence and global effectivity (θ) obtained for RPI (MQ) and least square based displacement recovery techniques (linear triangular, L2 norm).
Mesh Size (1/h)FEM Error
(×10−3)
RPI (Mesh Free)RPI (Mesh Dependent) LS (Mesh Dependent)
Error (×10−3)EffectivityError (×10−3)EffectivityError (×10−3)Effectivity
1/45.3653.0840.630042.4290.762093.7910.70817
1/160.3780.1920.647540.1400.833180.1810.67755
1/320.09510.0470.646150.0340.840880.0430.67375
Rate of Conv.1.939392.01685 2.06494 2.15207
Table 2. Error convergence and global effectivity (θ) obtained for RPI (MQ) and least square based displacement recovery techniques (quadratic triangular, L2 norm).
Table 2. Error convergence and global effectivity (θ) obtained for RPI (MQ) and least square based displacement recovery techniques (quadratic triangular, L2 norm).
Mesh Size (1/h)FEM Error
(×10−3)
RPI (Mesh Free)RPI (Mesh Dependent) LS (Mesh Dependent)
Error (×10−3)EffectivityError (×10−3)EffectivityError (×10−3)Effectivity
1/40.245010.247551.300500.146120.953021.030154.24214
1/120.008420.002241.021950.002021.000230.019302.44485
1/240.001040.000141.010170.000141.004940.001321.58496
Rate of Conv.3.047334.16168 3.86833 3.721499
Table 3. Error convergence and global effectivity (θ) obtained for RPI (MQ) and least square based displacement recovery techniques (quadrilateral, L2 norm).
Table 3. Error convergence and global effectivity (θ) obtained for RPI (MQ) and least square based displacement recovery techniques (quadrilateral, L2 norm).
Mesh Size (1/h)FEM Error
(×10−3)
RPI (Mesh Free)RPI (Mesh Dependent) LS (Mesh Dependent)
Error (×10−3)EffectivityError (×10−3)EffectivityError (×10−3)Effectivity
1/42.0322.2801.643552.8392.024551.8401.83209
1/160.1240.0971.555650.1762.288360.1111.82651
1/320.0310.0281.511510.0442.310240.0281.82654
Rate of Conv.2.012232.26171 2.00883 2.01729
Table 4. Error convergence and global effectivity (θ) obtained for RPI with radial basis function types and least square based displacement recovery techniques (linear triangular, energy norm).
Table 4. Error convergence and global effectivity (θ) obtained for RPI with radial basis function types and least square based displacement recovery techniques (linear triangular, energy norm).
Mesh Size (1/h)FEM Error
(×10−3)
RPI (Mesh Free)RPI (Mesh Dependent) LS (Mesh Dependent)
MQ Exp TSP MQ Exp Error
(×10−3)
θ
Error (×10−3)θError (×10−3)θError (×10−3)θError (×10−3)θError (×10−3)θ
1/493.74917.1140.946958.4760.856830.0790.9227626.2330.893351.1060.909658.4780.9300
1/1624.4411.3070.99317.3850.93891.5630.984961.5861.00388.2981.08604.4620.9896
1/3212.2480.3350.99752.6840.93100.5610.989240.4251.00923.8471.11951.1580.9970
Rate of Conv.0.978751.89163 1.48183 1.91524 1.98282 1.24389 1.89225
Table 5. Error convergence and global effectivity (θ) obtained for RPI with radial basis function types and least square based displacement recovery techniques (quadratic triangular, energy norm).
Table 5. Error convergence and global effectivity (θ) obtained for RPI with radial basis function types and least square based displacement recovery techniques (quadratic triangular, energy norm).
Mesh Size (1/h)FEM Error
(×10−3)
RPI (Mesh Free)RPI (Mesh Dependent) LS (Mesh Dependent)
MQExpTSPMQExp Error
(×10−3)
θ
Error (×10−3)θError (×10−3)θError (×10−3)θError (×10−3)θError (×10−3)θ
1/413.1687.4931.159635.0322.86375.7821.106553.9411.02944.4451.054812.2581.2556
1/121.5150.1841.01350.5601.08530.1991.018360.1741.00930.2571.01620.5291.0405
1/240.3800.0251.00570.06781.02750.0231.002920.0261.00530.06131.01160.0701.0127
Rate of Conv.1.978483.1910 3.48633 3.09017 2.76963 2.39122 2.88449
Table 6. Error convergence and global effectivity (θ) obtained for RPI with radial basis function types and least square based displacement recovery techniques (quadrilateral, energy norm).
Table 6. Error convergence and global effectivity (θ) obtained for RPI with radial basis function types and least square based displacement recovery techniques (quadrilateral, energy norm).
Mesh Size (1/h)FEM Error
(×10−3)
RPI (Mesh Free)RPI (Mesh Dependent) LS (Mesh Dependent)
MQExpTSPMQExpError
(×10−3)
θ
Error (×10−3)θError (×10−3)θError (×10−3)θError (×10−3)θError (×10−3)θ
1/460.26213.2491.054152.2220.817327.5951.012123.0531.026346.7920.856013.6671.0810
1/1615.0150.7970.99677.2721.00631.2720.98051.2941.01306.5930.96690.8131.0051
1/327.5060.2010.99672.5361.01910.6220.96150.3231.01642.3030.97820.2011.0013
Rate of Conv.1.001692.01497 1.45472 1.82359 2.05275 1.44822 2.02723
Table 7. Error convergence and global effectivity (θ) obtained for RPI (MQ) and least square based displacement recovery techniques (linear triangular, L2 norm).
Table 7. Error convergence and global effectivity (θ) obtained for RPI (MQ) and least square based displacement recovery techniques (linear triangular, L2 norm).
Mesh Size FEM
(Exact Error)
(×10−3)
RPI (MQ) RecoveryLS Recovery
Circular Mesh Free PatchRectangular Mesh Free PatchMesh Dependent PatchMesh Dependent Patch
Elem.DOF
Error
(×10−3)
θError
(×10−3)
θError
(×10−3)
θError
(×10−3)
θ
1031361.5340.6300.838650.6350.823290.6030.832640.9150.88585
4705300.2820.1170.934730.1160.910670.1150.971210.1410.93646
188720040.0700.0290.945030.0290.912260.0280.990680.0290.91340
Table 8. Error convergence and global effectivity (θ) obtained for RPI (MQ) and least square based displacement recovery techniques (quadratic triangular, L2 norm).
Table 8. Error convergence and global effectivity (θ) obtained for RPI (MQ) and least square based displacement recovery techniques (quadratic triangular, L2 norm).
Mesh Size FEM
(Exact Error)
(×10−3)
RPI (MQ) RecoveryLS Recovery
Circular Mesh Free PatchRectangular Mesh Free PatchMesh Dependent PatchMesh Dependent Patch
Elem.DOF
Error
(×10−3)
θError
(×10−3)
θError
(×10−3)
θError
(×10−3)
θ
1034760.0560.0390.952360.0340.947960.0320.879510.1212.38461
47019980.2470.00280.962880.00250.966050.00250.901080.00731.66566
188777800.0590.000330.986170.000300.980790.000310.930460.000541.27775
Table 9. Error convergence and global effectivity (θ) obtained for RPI (MQ) and least square based displacement recovery techniques (quadrilateral, L2 norm).
Table 9. Error convergence and global effectivity (θ) obtained for RPI (MQ) and least square based displacement recovery techniques (quadrilateral, L2 norm).
Mesh Size FEM
(Exact Error)
(×10−3)
RPI (MQ) RecoveryLS Recovery
Circular Mesh Free PatchRectangular Mesh Free PatchMesh Dependent PatchMesh Dependent Patch
Elem.DOF
Error
(×10−3)
θError
(×10−3)
θError
(×10−3)
θError
(×10−3)
θ
1864320.3560.3031.435640.2981.434440.2991.482730.2491.40829
59312820.1330.0941.383430.0921.368570.0931.381330.0751.28740
133328060.0540.0401.372380.0401.367110.0411.388470.0341.30001
Table 10. Error convergence and global effectivity (θ) obtained for RPI (MQ) and least square based displacement recovery techniques (quadrilateral energy norm).
Table 10. Error convergence and global effectivity (θ) obtained for RPI (MQ) and least square based displacement recovery techniques (quadrilateral energy norm).
Mesh Size FEM
(Exact Error)
(×10−3)
RPI (MQ) RecoveryLS Recovery
Circular Mesh Free PatchRectangular Mesh Free PatchMesh Dependent PatchMesh Dependent Patch
Elem.DOF
Error
(×10−3)
θError
(×10−3)
θError
(×10−3)
θError
(×10−3)
θ
18643222.7648.2691.049258.3491.042678.5791.088696.3051.06260
593128213.9173.8591.071864.3481.053534.0891.082021.8211.01123
181337907.9942.4101.054542.4521.031802.6351.088251.0111.00801
Table 11. Error convergence and global effectivity (θ) for plate problem with hole using RPI (MQ) with various mesh free node zones and least square based recovery techniques (linear triangular, energy norm).
Table 11. Error convergence and global effectivity (θ) for plate problem with hole using RPI (MQ) with various mesh free node zones and least square based recovery techniques (linear triangular, energy norm).
Mesh Size FEM
(Exact Error)
(×10−3)
RPI (MQ) RecoveryLS Recovery (Mesh Dependent)
Circular ZoneRectangular ZoneMesh Dependent PatchElement PatchNode Patch (ZZ)
[18]
Elem.DOF
Error
(×10−3)
θError
(×10−3)
θError
(×10−3)
θError
(×10−3)
θError
(×10−3)
θ
15519412.9057.6360.785527.7680.784858.2530.7312211.2060.8054113.2311.01807
5486188.5044.4420.918934.4020.940744. 8140.898515.4950.906487.9131.12593
158517005.1712.2510.935672.4500.967812.2440.934472.5720.900424.6721.14923
Table 12. Error convergence and global effectivity (θ) for plate problem with hole using exponential (Exp) and thin plate splint (TPS) radial basis function-based RPI error recovery techniques (linear triangular, energy norm).
Table 12. Error convergence and global effectivity (θ) for plate problem with hole using exponential (Exp) and thin plate splint (TPS) radial basis function-based RPI error recovery techniques (linear triangular, energy norm).
Mesh SizeFEM Error
(×10−3)
RPI (Exp)RPI (TPS)
Mesh FreeMesh DependentMesh FreeMesh Dependent
Elem.DOF
Error (×10−3)θError (×10−3)θError (×10−3)θError (×10−3)θ
15519412.90511.7391.0606110.8451.203908.2560.7911813.6110.92793
5486188.5047.2970.899076.5871.299554.6530.918355.9260.88131
158517005.1718.1081.553894.7830.957502.1570.924394.0950.95173
Table 13. Error convergence and global effectivity (θ) for plate problem with hole using RPI (MQ) and least square based recovery techniques (quadrilateral, energy norm).
Table 13. Error convergence and global effectivity (θ) for plate problem with hole using RPI (MQ) and least square based recovery techniques (quadrilateral, energy norm).
Mesh Size FEM
(Exact Error)
(×10−3)
RPI (MQ) RecoveryLS Recovery (Mesh Dependent)
Circular Mesh free Mesh Dependent PatchElement PatchNode Patch (ZZ)
[18]
Elem.DOF
Error
(×10−3)
θError
(×10−3)
θError
(×10−3)
θError
(×10−3)
θ
1794148.0855.3771.011835.3770.926327.2860.940189.2621.40135
57812544.2912.1960.942012.4200.941672.9490.999332.4040.89920
184238542.8821.6040.945511.6960.896991.7030.9455216460.91160
Table 14. Global errors using various recovery-based finite element analysis and number of elements (N) with DOF in adaptively improved meshes for 2% target error.
Table 14. Global errors using various recovery-based finite element analysis and number of elements (N) with DOF in adaptively improved meshes for 2% target error.
Recovery TypeLinear Triangle (Uniform Mesh Having 548 Elements and 618 DOF)Linear Quadrilateral (Uniform Mesh Having 179 Elements and 414 DOF)
FEM
Error
Proj. ErrorAdaptive Mesh (2%)FEM
Error
Proj. ErrorAdaptive Mesh (2%)
NDOFNDOF
LS (Mesh Dependent)6.846.227818666.506.1414443032
LS (Mesh Dependent-ZZ)6.847.6692210126.509.2518163784
RPI (Mesh Free-MQ-Cir.)6.846.32100410986.506.4214683082
RPI (Mesh Free-MQ-Rec.)6.846.47105011406.505.9614222990
RPI (Mesh Dependent-MQ)6.846.18117512746.506.3914523028
RPI (Mesh Free-Exp.-Cir.)6.846.14417486
RPI (Mesh Dependent-Exp.)6.848.9011551242
RPI (Mesh Free-TSP-Cir.)6.846.319451038
RPI (Mesh Dependent-TSP)6.846.0611001200
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Ahmed, M.; Singh, D.; AlQadhi, S.; Alrefae, M.A. Mesh Free Radial Point Interpolation Based Displacement Recovery Techniques for Elastic Finite Element Analysis. Mathematics 2021, 9, 1900. https://doi.org/10.3390/math9161900

AMA Style

Ahmed M, Singh D, AlQadhi S, Alrefae MA. Mesh Free Radial Point Interpolation Based Displacement Recovery Techniques for Elastic Finite Element Analysis. Mathematics. 2021; 9(16):1900. https://doi.org/10.3390/math9161900

Chicago/Turabian Style

Ahmed, Mohd., Devinder Singh, Saeed AlQadhi, and Majed A. Alrefae. 2021. "Mesh Free Radial Point Interpolation Based Displacement Recovery Techniques for Elastic Finite Element Analysis" Mathematics 9, no. 16: 1900. https://doi.org/10.3390/math9161900

APA Style

Ahmed, M., Singh, D., AlQadhi, S., & Alrefae, M. A. (2021). Mesh Free Radial Point Interpolation Based Displacement Recovery Techniques for Elastic Finite Element Analysis. Mathematics, 9(16), 1900. https://doi.org/10.3390/math9161900

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop