Next Article in Journal
A Symmetrical Terahertz Triple-Band Metamaterial Absorber Using a Four-Capacitance Loaded Complementary Circular Split Ring Resonator and an Ultra-Thin ZnSe Substrate
Previous Article in Journal
Symmetric Mass Generation
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Hybrid RBF Collocation Method and Its Application in the Elastostatic Symmetric Problems

1
School of Rail Transportation, Soochow University, Suzhou 215131, China
2
School of Transportation and Civil Engineering, Nantong University, Nantong 226019, China
*
Author to whom correspondence should be addressed.
Symmetry 2022, 14(7), 1476; https://doi.org/10.3390/sym14071476
Submission received: 14 June 2022 / Revised: 12 July 2022 / Accepted: 15 July 2022 / Published: 19 July 2022
(This article belongs to the Section Engineering and Materials)

Abstract

:
In this paper, a new hybrid radial basis function collocation method (HRBF-CM) is proposed to help resolve two-dimensional elastostatic symmetric problems. In the new approach, the hybrid radial basis function (HRBF) combines the infinitely smooth RBF and piecewise smooth RBF, containing two parameters (the shape parameter and the weight parameter). Discretization schemes are presented in detail. We use MATLAB to implement the HRBF-CM and produce numerical results which demonstrate the potential of this method. The new method’s accuracy is higher than that of the traditional methods, especially in the case of a more significant number of nodes. We discuss the new method’s effectiveness compared to the widely used traditional RBF and also investigate the effect of parameters on the method’s performance under the new method.

1. Introduction

The radial basis function was proposed by Hardy in the 1970s [1] and has made outstanding contributions to surface fitting and solving interpolation problems. Several authors discretized the partial differential equations (PDEs) using the collocation method [2,3]. The RBF method has the following three advantages [4]: (i) it can be estimated without using a mesh; (ii) it is highly accurate; and (iii) it has sufficient flexibility in choosing basis functions. Generally, a radial function ϕ ( r j ) is a function of the Euclidean norm r j = x x j 2 , where x n is the center point and x n is a point in the influence domain of x .
Most radial basis functions can be divided into two categories (see some typical RBFs in Table 1). The first category uses infinitely smooth radial basis functions, including multiquadric (MQ), inverse multiquadric (IMQ), and Gaussian (EXP) functions. This type of radial basis function contains a shape parameter c, which needs to be defined by the user and controls the stability and accuracy of the RBF approximation. Therefore, choosing an optimal shape parameter value is imperative, which has also spurred discussion among many scholars. Generally, we prefer a constant optimal shape parameter (CSP) c throughout the computation, and the shape parameter c is the same for each node.
Some famous scholars proposed an experienced formula to select the optimal shape parameters [5,6,7]. Rippa chooses the best shape parameter c based on the minimization of the estimation of the error function and proposes the leave-one-out cross-validation (LOOCV) method [8]. However, it was assumed that CSP is not suitable. In this case, the approximation accuracy of the interpolation method was expected to be reduced, even more so than the approximation method of the general function, and the condition number of the matrix was large, resulting in an ill-conditioned linear system [9,10]. Kansa and Sarra used the variable shape parameter (VSP) form to solve this problem [11,12]. VSP is used to assign different values to the shape parameters which correspond to each center. This shows that each row of the RBF matrix has other shape parameters, which makes the condition number reduction in the RBF matrix more stable and improves the approximation of the interpolation function. Another method was used to determine the interval of the shape parameter, rather than a value, without any minimization and estimation error [13]. The purpose of VSP is to convert several CSPs into an interval range. Some other VSP strategies are summarized in [14,15,16].
It can be seen that choosing optimal shape parameter values in infinitely smooth RBFs is essential and has been an ongoing challenge. In contrast, the second set of piecewise smooth RBFs is not infinitely differentiable. They have no shape parameter and are easier to implement than the former group when used. However, they only lead to algebraic convergence rates and are rarely used alone in applications [17]. Therefore, many researchers prefer using infinitely smooth RBFs in scientific computing and engineering applications [18,19,20,21].
Some scholars bypassed the process of seeking optimal shape parameters and proposed a nontraditional RBF. For instance, Karimi [22] introduces a new Gaussian RBF (GRBF), which is both stable and accurate and is free of instability issues for small values of the shape parameter. Similarly, Zhang [17] proposed a new global RBF based on the coupling of infinitely smooth RBF with the conical spines, known as the coupled radial basis function (CRBF). He used the CRBF method in the Kansa method to produce a relatively good linear system. The condition number is moderate, and the error changes slightly with the change in the shape parameter value. Cao [23] uses this CRBF method to reconstruct the elasticity problem. However, there are still some deficiencies in the calculation accuracy, and the operation time is too long. Recently, Manzoor [24] used this similar idea to propose a hybrid radial basis function (HRBF) method, which combines two different RBFs: an infinite smooth RBF defined by shape parameters and a piecewise smooth RBF that is independent of shape parameters. This method has better accuracy when solving Burgers’ equations. In this paper, we use this method to solve elastostatic symmetric problems.
The predecessors have also made great efforts to solve elastostatic problems [25,26,27,28,29,30,31,32,33,34,35]. This paper mainly studies two well-known examples: cantilever beams and thin plates with circular holes. Zhang [25] applied the globally supported RBFs, i.e., MQ, IMQ, and Wendland’s functions [26]. Tolstykh and Shirobokov [27] applied the local RBF method. Liu et al. proposed a combination of Galerkin weak form and RBFs, namely the radial point interpolation method (RPIM) [28,29,30,31]. Simonenko [32] analyzes the applicability of both the global and the local versions of the method for elastostatic problems. They use multiquadric functions as RBFs and describe how to select an optimal shape parameter value to minimize approximation errors. Based on Simonenko’s article, we use the HRBF instead of MQ-RBF and the global Kansa method to conduct numerical experiments.
This paper is structured as follows. In Section 2, we first introduce the origin of the hybrid radial basis function, list the several types of HRBF, and describe the HRBF-CM’s formulation to find a solution to elastostatic problems. In Section 3, three numerical examples are calculated to verify the validity and accuracy of the HRBF collocation method, namely a patch test, a cantilever beam, and an infinite plate with a central hole. We compare the accuracy of the HRBF-CM with the infinitely smooth RBFs and the piecewise smooth RBFs. The distribution of two kinds of regular and irregular nodes is considered, and the influence of shape and weight parameters on the calculation results is discussed. The last section provides a summary of this method.

2. The HRBF-CM for Elastostatic Symmetric Problems

In this section, we describe the origin of the HRBF, the HRBF used in this article, and the process of solving elastostatic equations using the HRBF-CM in detail.

2.1. The Construction of HRBFs

As presented in [24], by coupling the infinite smoothing RBFs ( ϕ 1 ), such as MQ, IMQ, and GA, with the piecewise smooth RBFs Sn and TPS RBFs ( ϕ 2 ), the HRBF is constructed and can be defined as
R ( r j ) = γ ϕ 1 ( r j ) + ϵ ϕ 2 ( r j ) , r j = x x j 2 ,
where the parameters γ and ϵ are positive real numbers that control the contributions of ϕ 1 and ϕ 2 , respectively. ϕ 1 is an infinite smooth RBF with a user-defined shape parameter c , and ϕ 2 is a piecewise smooth RBF free from shape parameter.
Here, we find that the HRBF has three parameters, which seem to increase the algorithm’s complexity. At this point, we can scale the HRBF by a constant, which will not affect the algorithm. Therefore, the HRBF (1) can be normalized by the constant w = ϵ / γ and has been redefined as
R w ( r j ) = ϕ 1 ( r j ) + w ϕ 2 ( r j ) .
Note that the reconstructed HRBF (2) now contains two parameters (the shape parameter c and the weight parameter w ). These two parameters control the accuracy and stability of the HRBF-CM. The HRBFs to be used in this paper are listed in Table 2. It has been dictated in Mishra’s article [33] that a small value of the weight parameter w usually gives acceptable accuracy with a stable spectrum (when all the eigenvalues of the differentiation matrix lie in the left half-plane). In Manzoor’s article [24], the authors calculated many numerical examples and concluded that in most of the test examples, for the three types of hybrid radial basis functions, MQ + S3, GS + S3, and IMQ + S3, the weight parameter w = 10−3 works best. On the other hand, for the three types of hybrid radial basis functions, MQ + TPS, GS + TPS, and IMQ + TPS, the weight parameter w = 10−9. However, in this work, we use w = [100, 10−1, 10−2, 10−3, …, 10−15] to find the optimal parameter w.

2.2. The Elastostatic Equation

Considering the two-dimensional elastostatic symmetric problem, Ω and Γ represent the problem domain and the boundary, respectively. The equations of elasticity, written in terms of displacements, are defined as
E 1 μ 2 ( 2 u x x 2 + 1 μ 2 2 u x y 2 + 1 + μ 2 2 u y x y ) = f x , E 1 μ 2 ( 2 u y y 2 + 1 μ 2 2 u y x 2 + 1 + μ 2 2 u x x y ) = f y ,       i n         Ω
with the displacement boundary condition
u x = u ¯ x , u y = u ¯ y , o n Γ u
and the stress boundary condition
l σ x x + m τ x y = t ¯ x , l τ y x + m σ y y = t ¯ y     o n       Γ t
where ux and uy represent the displacement components on the boundary; fx and fy are the given body forces; u ¯ and v ¯ are the known functions of nodes on the displacement boundary Γ u ; l and m represent the cosine of the normal direction outside the slope; and σ x x ,   σ y y ,   τ x y represent the stress constraints on the stress boundary Γ t . E and μ stand for the Young’s modulus and the Poisson’s ratio, respectively. The stresses σ x x , σ y y , and τ x y and the displacements ux and uy present the following relationship.
{ σ x x = E 1 μ 2 ( u x x + μ u y y ) , σ y y = E 1 μ 2 ( μ u x x + u y y ) , τ x y = E 2 ( 1 + μ ) ( u x y + u y x ) .

2.3. The HRBF Formulation

Assuming that there are N distinct nodes { x i } i = 1 N , the HRBF-CM is used to solve the elastostatic symmetric problem (3)–(5). In the HRBF-CM, the displacement approximate solutions u x ( x ) and u y ( x ) are separately spanned by a set of translated HRBFs, i.e.,
{ u x ( x ) = j = 1 N a j R w ( r j ) , u y ( x ) = j = 1 N b j R w ( r j ) .
R w ( r j ) with r j = x x j 2 is the HRBF defined in Table 2. The unknown coefficients, aj and bj, are computed by a collocation of Equation (3) at a set of interior nodes and a collocation of the boundary conditions at boundary nodes. Let z be defined as a vector of length 2N containing the unknowns.
{ z j = a j z j + N = b j , j = 1 , , N ,
Then, a vector f is defined,
{ f j = f x ( x j ) f j + N = f y ( x j ) j = 1 , , N .
By substituting Equation (7) into the elastostatic Equation (3) and boundary conditions (4)–(5), we obtain the following linear system:
G z = f .
For i = 1, …, N, the elements of j internal nodes corresponding to the matrix G are
{ G i , j ( i j ) = E 1 μ 2 [ 2 R w ( r i , j ) x 2 + 1 μ 2 2 R w ( r i , j ) y 2 ] , G i , j ( i = j ) = E 1 μ 2 [ 2 ϕ ( r i , j ) x 2 + 1 μ 2 2 ϕ ( r i , j ) y 2 ] .
{ G i , j + N ( i j + N ) = E 2 ( 1 μ ) 2 R w ( r i , j ) x y , G i , j + N ( i = j + N ) = E 2 ( 1 μ ) 2 ϕ ( r i , j ) x y .      
{ G i + N , j ( i + N j ) = E 2 ( 1 μ ) 2 R w ( r i , j ) x y , G i + N , j ( i + N = j ) = E 2 ( 1 μ ) 2 ϕ ( r i , j ) x y .      
{ G i + N , j ( i + N j + N ) = E 1 μ 2 [ 2 R w ( r i , j ) y 2 + 1 μ 2 2 R w ( r i , j ) x 2 ] , G i + N , j ( i + N = j + N ) = E 1 μ 2 [ 2 ϕ ( r i , j ) y 2 + 1 μ 2 2 ϕ ( r i , j ) x 2 ] .
The elements of the boundary nodes of the matrix G depend on the boundary conditions. For instance, in the case of Dirichlet boundary conditions,
{ G i , j = R w ( r i , j ) , f j = u ¯ x ( x j ) ,                   G i + N , j + N = R w ( r i , j ) , f j + N = u ¯ y ( x j ) ,        
where u ¯ x and u ¯ y are the exact solutions of problem (3), and Equation (10) can be written as
( G i , j G i , j + N G i + N , j G i + N , j + N ) ( a j b j ) = ( u ¯ x ( x j ) u ¯ y ( x j ) ) .
Therefore, the coefficients aj and bj of Equation (7) are used to compute the horizontal and vertical displacements ux and uy, respectively. In addition, the stresses are computed from Equation (5) with
u x x = i = 1 N a j R w ( r j ) x ,
u x y = i = 1 N a j R w ( r j ) y ,
u y x = i = 1 N b j R w ( r j ) x ,
u y y = i = 1 N b j R w ( r j ) y .

3. Numerical Examples

In this section, we use three numerical examples to verify the validity and the applicability of the hybrid radial basis function collocation method. The numerical results for these examples are compared with the analytical solutions and the traditional RBF solutions. The effect of different basis functions on computational accuracy is discussed.
To judge the validity and accuracy of the HRBF-CM method, we need a measure which can be measured by the relative error concerning the displacement and the stress:
{ r u = i = 1 N ( u i n u m u i a n a l ) 2 / i = 1 N u i a n a l 2 r σ = i = 1 N ( σ i n u m σ i a n a l ) 2 / i = 1 N σ i a n a l 2
where u i n u m and σ i n u m are the numerical solutions for the displacement and stress of the ith node, respectively, and u i a n a l and σ i a n a l are the corresponding analytical solutions for the ith node.

3.1. Patch Test

The first numerical example is the standard patch test containing 13 nodes, with 5 irregular interior nodes, as shown in Figure 1. In this patch test, the essential boundary is the linear displacement assigned along the boundary. The linear displacement functions are ux = 0.6x and uy = 0.6y. Satisfying the patch test requires the same linear function to give the displacement of any interior node and requires the strains and stresses to be constant in the patch. In this problem, we randomly choose the shape parameter c = 1, the weight parameter referring to Manzoor’s article [24], the variant related to TPS, the weight parameter w = 10−9, the variant associated with S3, and the weight parameter w = 10−3. The HRBF passes the patch test with precision, as well as the traditional RBF. Different radial basis functions have different accuracies. The exact results of the internal node calculations for the MQ-HRBF and the traditional MQ-RBF are shown in Table 3 when the parameters are taken as E = 1 and v = 0.3.

3.2. Cantilever Beam

For the second example, we consider a cantilever beam problem, as presented in Figure 2. A beam of length L and height H subjected to traction at the free end is considered. The beam has a unit thickness; hence, a plane stress problem is also considered here. The closed-form solution is available for the parabolic traction of force P:
{ u x a n a l = P 6 E I ( y H 2 ) [ ( 6 L 3 x ) x + ( 2 + μ ) ( y 2 H y ) ] , u y a n a l = P 6 E I [ 3 μ ( y H 2 ) 2 ( L x ) + ( 4 + 5 μ ) H 2 x 4 + ( 3 L x ) x 2 ] ,
where the moment of inertia I of the beam is given as I = H3/12.
The stresses corresponding to the above displacements are
{ σ x x a n a l = P ( L x ) I ( y H 2 ) , σ y y a n a l = 0 , τ x y a n a l = P 2 I ( y 2 H y ) .
We refer to the study by Simonenko et al. [32] which uses the Dirichlet boundary conditions given by (22) in all the boundaries of the beam. The relevant parameters are L = 12; H = 2; E = 1000; μ = 0.3; p = −10.

3.2.1. Effect of Irregular Node Distribution

To investigate the stability of the HRBF-CM and its comparison with the traditional radial basis function method, both regular and irregular node distributions are used. Four kinds of regular node distributions, i.e., 13 × 3 (13 nodes in the x direction, 3 nodes in the y direction), 25 × 5, 49 × 9, and 73 × 13, are considered. Four irregular distributions are used with the same number of nodes, i.e., 39, 125, 441, and 949. We discuss the effect of the node distributions on the numerical accuracy, and the results of the calculations are presented in Table 4 and Table 5. The two node distributions for the 441 nodes are shown in Figure 3a,b. The deflection distributions for these two node distribution methods are given in Figure 4.
The HRBF-CM can obtain very accurate approximate solutions for the displacements, indicating its high applicability and accuracy in the case of irregular node distribution methods, which perform better than the traditional radial basis function. The accuracy of the IMQ + S3 method reaches up to 10−9 for displacements and 10−16 for stresses. As the number of nodes increases, the results of the HRBF-CM are in better agreement with the analytical solution. Nevertheless, the number of nodes is not as high as possible, as too many nodes can lead to higher computational costs and lower efficiency. As can be seen from Table 4 and Table 5, the deformations of the GA radial basis functions (GA + TPS, GA + S3), especially GA + S3, have better computational accuracy and stability for different node distribution patterns. Table 4 and Table 5 also list the results of sharded smooth TPS and S3. It can be seen that they do not include the choice of parameters, and the calculation time is short, but the accuracy is not high.

3.2.2. Effect of Parameters

When the meshless method analyzes nonlinear elastic problems, different radial basis function shape parameters impact the accuracy of the calculation. In addition to the shape parameter, the HRBF-CM includes a weight parameter w, so this section goes into detail. The influence of the shape parameter c and the weight parameter w on the calculation results is studied. The relationships between r u ,     r σ and the shape parameter c for the regular 73 × 13 (949 in total) nodes and the irregular 949 nodes are plotted in Figure 5. In particular, Figure 5a,b plot the relationship between the shape parameter c (varying from 10−4∼20) and the relative error of displacement r u and Figure 5c,d plot the corresponding relative error of r σ . To explore the variation of other nodes, Figure 6 plots ru and r σ relative to the weight parameter w for the regular 49 × 9 (441 in total) nodes and the irregular 441 nodes, where Figure 6a,b plot the relationship between the weight parameter w ([100, 10−1, 10−2, …, 10−15]) and the relative error of displacement ru, and Figure 6c,d plot the corresponding relative error of r σ .
We can first see from these figures that the HRBF-CM has high accuracy. Figure 5 shows that the shape parameter c has a particular influence on the calculation progress of the conventional and hybrid radial basis functions. However, when the shape parameter takes an enormous value, the relative error of the hybrid radial basis function accuracy is still excellent. The trend shows that as the shape parameter increases, the relative error is smaller and more accurate. The hybrid radial basis function also contains a weighting parameter w. This article examines the effect of the weight parameter on accuracy in detail. We select 16 weighting parameters, and Figure 6 shows that the accuracy of all hybrid radial basis functions becomes higher as the weighting parameters are reduced, showing a convergence trend. The parameters for the TPS hybrid radial basis functions are smoother and not as volatile as those for the S3 hybrid radial basis functions. Table 4 and Table 5 give the optimal parameters for the regular and irregular nodes, and the optimal parameters are chosen by minimizing the relative error in the displacement ru. These numerical results show that the HRBF-CM has better accuracy and better stability than the traditional RBF method when solving cantilever beam problems.

3.2.3. Convergence Analysis and Simulation

In this subsection, we analyze the convergence of the HRBF-CM using the cantilever beam as an example, and also give simulated plots of the computational results. In Table 4 and Table 5, we give the minimum relative error ru and r σ for regular and irregular nodes, respectively, using different calculation methods. We draw the curves of the convergence diagram, as shown in Figure 7. Figure 7a,b show the convergence curves for stresses and displacements with regular node distribution. Figure 7c,d show the same items under irregular nodal distribution. h is the maximum size of the nodal arrangement.
The calculation results in Table 4 and Table 5 show that the HRBF-CM has excellent convergence properties. TPS and S3 converge slowly, and MQ, IMQ, and GA do not show good convergence. The relative errors ru and r σ do not decrease as the number of nodes increases. For the small number of nodes, the accuracy of the traditional infinitely smooth RBF method is similar to the HRBF-CM. However, in the case of a large number of nodes, the HRBF-CM is more advantageous.
We plot the analytical and numerical solutions of the HRBF-CM (MQ + TPS) in the case of 73 × 13 regular nodes, as shown in Figure 8. Figure 8a represents the normal stress σ x x and Figure 8b shows the shear stress τ x y . It can be seen that the numerical solution and the analytical solution are highly consistent, and the calculation accuracy is very high.

3.3. Plate with a Hole

For the third test problem, we consider a center circular hole plate with uniform tensile load σ 0 at infinity in the x direction, as shown in Figure 9a. In actual computations, the geometric parameters and the material parameters are as follows: σ 0 = 1 , Lx = 5, Ly = 5, a = 1, E = 1000, and μ = 0.3. Due to its symmetry, only a quarter of the model is considered in the calculation, as shown in Figure 9b.
We discuss the third numerical problem in Cartesian coordinates. This problem can also be discussed in polar coordinates [34]. Here, we only consider the Cartesian coordinate system. The analytical solution for displacement and stresses is as follows:
{ u x a n a l = σ 0 a 8 G [ x a ( κ + 1 ) + 2 a ( 1 + κ ) x x 2 + y 2 + 2 a x 3 3 x y 2 ( x 2 + y 2 ) 2 ( 1 a 2 x 2 + y 2 ) ] , u y a n a l = σ 0 a 8 G [ y a ( κ 3 ) + 2 a ( 1 κ ) y x 2 + y 2 + 2 a 3 y x 2 y 3 ( x 2 + y 2 ) 2 ( 1 a 2 x 2 + y 2 ) ] ,
{ σ x x a n a l = σ 0 [ 1 3 2 a 2 ( x 2 y 2 ) ( x 2 + y 2 ) 2 a 2 ( x 2 + y 2 ) 3 ( x 4 + y 4 6 x 2 y 2 ) ( 1 3 2 a 2 x 2 + y 2 ) ] , σ y y a n l a = σ 0 [ 1 2 a 2 ( x 2 y 2 ) ( x 2 + y 2 ) 2 + a 2 ( x 2 + y 2 ) 3 ( x 4 + y 4 6 x 2 y 2 ) ( 1 3 2 a 2 x 2 + y 2 ) ] , τ x y a n a l = σ 0 [ a 2 x y ( x 2 + y 2 ) 2 + 4 a 2 x y x 2 y 2 ( x 2 + y 2 ) 3 ( 1 3 2 a 2 x 2 + y 2 ) ] .
where G = E 2 ( 1 + μ ) is the shear modulus. κ = 3 μ 1 + μ . is the Kolosov constant.

3.3.1. Effect of Parameters

We calculate three-node distributions for the plate with a hole, including 289, 626, and 1089. We plot the relationship between the shape parameter c and the relative errors ru and r σ for the 1089 node distribution, respectively, in Figure 10. Figure 11 plots the relationship between the weight parameter w and the relative errors ru and r σ for the same-sized nodes. Due to the difference in the order of magnitude of the errors, the relative error variations for MQ + TPS, MQ + S3, IMQ + TPS, and IMQ + S3 are plotted in Figure 11a,b to provide a more intuitive view. The relative error variations for GA + TPS and GA + S3 are plotted in Figure 11c,d.
As can be seen in Figure 10, the traditional RBF is not as accurate as HRBF when the shape parameter gradually increases. The GA radial basis function is not as accurate, and its variants GA + TPS, as well as GA + S3, are not advantageous either. As seen in Figure 11a,b, the relative error is more affected by the variation of the weighting parameter w, with better results for weighting parameters between [10−5, 10−6, 10−7, 10−8, 10−9]. The GA is less accurate, as seen in Figure 11c,d.

3.3.2. Convergence Analysis and Simulation

We investigate the convergence of the HRBF-CM in solving the plate with a hole by considering three-node distributions, namely 289, 626, and 1089. Figure 12 plots the results of the displacement calculations for nodes 625 and 1089. It can be seen that the settlement accuracy is very high. The calculation results for the other nodes are shown in Table 6. The optimal parameters listed in Table 6 are obtained by minimising the relative error ru of the displacements. The calculation results show that the relative errors ru and r σ decrease as the nodes increase. Figure 13 effectively shows the convergence characteristics, plotting the convergence curves for the nine kinds of radial basis functions. From these results, we find that the deformation HRBFs for all types of conventional radial basis functions are more accurate than the traditional radial basis functions, i.e., MQ + TPS and MQ + S3 are better than the MQ method, IMQ + TPS and IMQ + S3 are better than the MQ method, and GA + TPS and GA + S3 are more accurate than GA. When the number of nodes is 1089, IMQ + S3 has a displacement accuracy of up to 10−5 and a stress accuracy of 10−3. The HRBF-CM also has the highest accuracy for the other node distributions.
Figure 14a,b plot the analytical solution and the HRBF (IMQ + S3) solution for stress σ x x . Figure 14c,d plot the analytical solution and the HRBF (IMQ + S3) solution for the stress σ y y . Figure 14e,f plot the analytical and HRBF (IMQ + S3) solutions for the shear stress τ x y . These figures also show the excellent agreement between the analytical and numerical results.
It should be pointed out that the present study only considers the global collocation method to solve small-scale problems. We can use the local collocation method to solve large-scale problems in the future. The numerical examples are mainly focused on plane elastostatic symmetry problems. However, it can be extended to three-dimensional elastostatic problems to verify the effectiveness of the HRBF-CM.

4. Conclusions

The traditional RBF methods are widely used in many areas of scientific computing, such as interpolating multi-dimensional scattered data and solving PDEs. The HRBF is constructed based on infinite smooth RBFs and piecewise smooth RBFs. It does inherit the advantages of both RBFs, making HRBFs both spectrally convergent and stable. Inspired by Manzoor, [24] used the HRBF to solve the Burgers’ equation; we study his HRBF in detail in this paper, using three elastostatic symmetric problems. The sensitivity analysis of the HBRF concerning the number of nodes N, the weight parameter w, and the shape parameter c is analyzed. The ability of the HRBF-CM to fit the data under different shape parameters and weight parameters is investigated. The numerical results of the HRBF are compared with those of the analytical solutions and other numerical methods, and excellent agreements can be found.

Author Contributions

Conceptualization, C.L. and Y.C.; methodology, C.L., Y.C. and Y.-T.C.; software, Y.-T.C.; validation, C.L., Y.C. and Y.-T.C.; formal analysis, L.-Q.Y.; investigation, Y.-T.C.; resources, C.L. and Y.C.; data curation, L.-Q.Y. and Y.-T.C.; writing—original draft preparation, C.L. and Y.-T.C.; writing—review and editing, L.-Q.Y. and Y.C.; visualization, Y.-T.C.; supervision, C.L., L.-Q.Y. and Y.C.; project administration, C.L.; funding acquisition, C.L. All authors have read and agreed to the published version of the manuscript.

Funding

The work was supported by the National Natural Science Foundation of China (No. 11972240).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data is contained within the present article. The data used in this article are available in each figure or table.

Acknowledgments

The authors express their gratitude to the anonymous reviewers for the constructive suggestions.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Hardy, R.L. Multiquadric equations of topography and other irregular surfaces. J. Geophys. Res. 1971, 76, 1905–1915. [Google Scholar] [CrossRef]
  2. Buhmann, M.D. Radial basis functions. Acta Numer. 2000, 9, 1–38. [Google Scholar] [CrossRef] [Green Version]
  3. Kansa, E.J. Multiquadrics—A scattered data approximation scheme with applications to computational fluiddynamics—II solutions to parabolic, hyperbolic and elliptic partial differential equations. Comput. Math. Appl. 1990, 19, 147–161. [Google Scholar] [CrossRef] [Green Version]
  4. Lee, M.B.; Lee, Y.L.; Sunwoo, Y.L.; Yoon, J. Some issues on interpolation matrices of locally scaled radial basis functions. Appl. Math. Comput. 2011, 217, 5011–5014. [Google Scholar] [CrossRef]
  5. Fasshauer, G.E.; Zhang, J.G. On choosing “optimal” shape parameters for RBF approximation. Numer. Algorithms 2007, 45, 345–368. [Google Scholar] [CrossRef]
  6. Franke, C.; Schaback, R. Solving partial differential equations by collocation using radial basis functions. Appl. Math. Comput. 1998, 93, 73–82. [Google Scholar] [CrossRef]
  7. Hardy, R.L. A contribution of the multiquadric method: Interpolation of potential inside the earth. Comput. Math. Appl. 1992, 24, 81–97. [Google Scholar] [CrossRef] [Green Version]
  8. Rippa, S. An algorithm for selecting a good value for the parameter c in radial basis function interpolation. Adv. Comput. Math. 1999, 11, 193–210. [Google Scholar] [CrossRef]
  9. Diederichs, B.; Iske, A. Improved estimates for condition numbers of radial basis function interpolation matrices. J. Approx. Theory 2019, 238, 38–51. [Google Scholar] [CrossRef]
  10. Majdisova, Z.; Skala, V. Radial basis function approximations: Comparison and applications. Appl. Math. Model. 2017, 51, 728–743. [Google Scholar] [CrossRef] [Green Version]
  11. Kansa, E.J.; Carlson, R. Improved accuracy of multi-quadric interpolation using variable shape parameters. Comput. Math. Appl. 1992, 24, 20–99. [Google Scholar] [CrossRef] [Green Version]
  12. Scotta, A.S.; Derek, S. A random variable shape parameter strategy for radial basis function approximation methods. Eng. Anal. Bound. Elem. 2009, 33, 1239–1245. [Google Scholar]
  13. Biazar, J.; Hosami, M. An interval for the shape parameter in radial basis function approximation. Appl. Math. Comput. 2017, 315, 131–149. [Google Scholar] [CrossRef]
  14. Ahmad, G.; Ehsan, M.; Hamed, R. On the new variable shape parameter strategies for radial basis function. Comput. Appl. Math. 2015, 34, 691–704. [Google Scholar]
  15. Stolbunov, V.; Nair, P.B. Sparse radial basis function approximation with spatially variable shape parameters. Appl. Math. Comput. 2018, 330, 170–184. [Google Scholar] [CrossRef]
  16. Bayona, V.; Moscoso, M.; Kindelan, M. Optimal variable shape parameter for multiquadric based RBF-FD method. J. Comput. Phys. 2012, 231, 2466–2481. [Google Scholar] [CrossRef]
  17. Zhang, Y.M. An accurate and stable RBF method for solving partial differential equations. Appl. Math. Lett. 2019, 97, 93–98. [Google Scholar] [CrossRef]
  18. Nadir, B.; Mohammed, O.; Minh-Tuan, N.; Abderrezak, S. Optimal trajectory generation method to find a smooth robot joint trajectory based on multiquadric radial basis functions. Int. J. Adv. Manuf. Technol. 2022, 120, 297–312. [Google Scholar] [CrossRef]
  19. Feng, L.T.; Katupitiya, J. Radial basis function-based vector field algorithm for wildfire boundary tracking with UAVs. Comput. Appl. Math. 2022, 41, 124. [Google Scholar] [CrossRef]
  20. Swetha, K.; Eldho, T.I.; Singh, L.G.; Kumar, A.V. Groundwater flow simulation in a confined aquifer using local radial point interpolation meshless method (LRPIM). Eng. Anal. Bound. Elem. 2022, 134, 637–649. [Google Scholar] [CrossRef]
  21. Parand, K.; Hemami, M. Numerical study of astrophysics equations by meshless collocation method based on compactly supported radial basis function. Int. J. APPL. Comput. Math. 2017, 3, 1053–1075. [Google Scholar] [CrossRef] [Green Version]
  22. Karimi, N.; Kazem, S.; Ahmadian, D.; Adibi, H.; Ballestra, L.V. On a generalized Gaussian radial basis function: Analysis and applications. Eng. Anal. Bound. Elem. 2020, 112, 46–57. [Google Scholar] [CrossRef]
  23. Chen, Y.T.; Cao, Y. A coupled RBF method for the solution of elastostatic problems. Math. Probl. Eng. 2021, 2121, 6623273. [Google Scholar] [CrossRef]
  24. Hussain, M. Hybrid radial basis function methods of lines for the numerical solution of viscous Burgers’ equation. J. Comput. Appl. Math. 2021, 40, 107. [Google Scholar] [CrossRef]
  25. Zhang, X.; Song, K.Z.; Lu, M.W.; Liu, X. Meshless methods based on collocation with radial basis functions. Comput. Mech. 2000, 26, 333–343. [Google Scholar] [CrossRef]
  26. Wendland, H. Piecewise polynomial positive definite and compactly supported radial basis functions of minimal degree. Adv. Comput. Math. 1995, 4, 389–396. [Google Scholar] [CrossRef]
  27. Tolstykh, A.I.; Shirobokov, D.A. On using radial basis functions in a “finite difference mode” with applications to elasticity problems. Comput. Mech. 2003, 33, 68–79. [Google Scholar] [CrossRef]
  28. Liu, G.R.; Gu, Y.T. A local radial point interpolation method (LRPIM) for free vibration analyses of 2-D solids. J. Sound Vib. 2001, 246, 29–46. [Google Scholar] [CrossRef] [Green Version]
  29. Wang, J.G.; 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]
  30. Wang, J.G.; Liu, G.R. 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]
  31. Cao, Y.; Yao, L.Q.; Yi, S.C. A weighted nodal-radial point interpolation meshless method for 2D solid problems. Eng. Anal. Bound. Elem. 2014, 39, 88–100. [Google Scholar] [CrossRef]
  32. Simonenko, S.; Bayona, V.; Kindelan, M. Optimal shape parameter for the solution of elastostatic problems with the RBF method. J. Eng. Math. 2014, 85, 115–129. [Google Scholar] [CrossRef]
  33. Mishra, P.K.; Nath, S.K.; Kosec, G.; Sen, M.K. An improved radial basis-pseudospectral scheme with hybrid Gaussian-cubic kernels. Eng. Anal. Bound. Elem. 2017, 80, 162–171. [Google Scholar] [CrossRef] [Green Version]
  34. Liu, Z.; Gao, H.F.; Wei, G.F.; Wang, Z.M. The meshfree analysis of elasticity problem utilizing radial basis reproducing kernel particle method. Results Phys. 2020, 17, 103037. [Google Scholar] [CrossRef]
  35. Chen, Z.Y.; Zhao, Q.F.; Lim, C.W. A new static-dynamic equivalence beam bending approach for the stability of a vibrating beam. Mech. Adv. Mater. Struct. 2021, 28, 999–1009. [Google Scholar] [CrossRef]
Figure 1. Example 1: a patch with 13 nodes.
Figure 1. Example 1: a patch with 13 nodes.
Symmetry 14 01476 g001
Figure 2. Example 2: cantilever beam.
Figure 2. Example 2: cantilever beam.
Symmetry 14 01476 g002
Figure 3. Example 2: (a) 441 regular node distribution; (b): 441 irregular node distribution.
Figure 3. Example 2: (a) 441 regular node distribution; (b): 441 irregular node distribution.
Symmetry 14 01476 g003
Figure 4. Example 2: displacement computational results using HRBF (GA + S3)−(a): 441 regular nodes distribution; (b): 441 irregular node distribution.
Figure 4. Example 2: displacement computational results using HRBF (GA + S3)−(a): 441 regular nodes distribution; (b): 441 irregular node distribution.
Symmetry 14 01476 g004
Figure 5. Example 2: relationship between relative errors ru and r σ and shape parameter c—(a,c): 949 regular node distribution; (b,d): 949 irregular node distribution.
Figure 5. Example 2: relationship between relative errors ru and r σ and shape parameter c—(a,c): 949 regular node distribution; (b,d): 949 irregular node distribution.
Symmetry 14 01476 g005
Figure 6. Example 2: relationship between relative errors and weight parameter—(a,c): 441 regular nodes; (b,d): 441 irregular nodes.
Figure 6. Example 2: relationship between relative errors and weight parameter—(a,c): 441 regular nodes; (b,d): 441 irregular nodes.
Symmetry 14 01476 g006
Figure 7. Example 2: convergence study−(a): regular node distribution, convergence curves of ru; (b): regular node distribution, convergence curves of r σ ; (c): irregular node distribution, convergence curves of ru; (d): irregular node distribution, convergence curves of r σ .
Figure 7. Example 2: convergence study−(a): regular node distribution, convergence curves of ru; (b): regular node distribution, convergence curves of r σ ; (c): irregular node distribution, convergence curves of ru; (d): irregular node distribution, convergence curves of r σ .
Symmetry 14 01476 g007
Figure 8. Example 2: stress computational results using HRBF (MQ + TPS, regular 949 nodes). (a) Normal stress σ x x ; (b): shear stress τ x y .
Figure 8. Example 2: stress computational results using HRBF (MQ + TPS, regular 949 nodes). (a) Normal stress σ x x ; (b): shear stress τ x y .
Symmetry 14 01476 g008
Figure 9. Example 3: (a): plate with a hole; (b): a quarter of the model.
Figure 9. Example 3: (a): plate with a hole; (b): a quarter of the model.
Symmetry 14 01476 g009
Figure 10. Example 3: 1089 node distribution−(a): relationship between relative errors ru and shape parameter c; (b): relationship between relative errors r σ and shape parameter c.
Figure 10. Example 3: 1089 node distribution−(a): relationship between relative errors ru and shape parameter c; (b): relationship between relative errors r σ and shape parameter c.
Symmetry 14 01476 g010
Figure 11. Example 3: 1089 node distribution—(a,c): relationship between relative errors ru and shape parameter w; (b,d): relationship between relative errors r σ and shape parameter w.
Figure 11. Example 3: 1089 node distribution—(a,c): relationship between relative errors ru and shape parameter w; (b,d): relationship between relative errors r σ and shape parameter w.
Symmetry 14 01476 g011
Figure 12. Example 3: displacement computational results using the HRBF (IMQ + S3)—(a): 625 nodes; (b): 1089 nodes.
Figure 12. Example 3: displacement computational results using the HRBF (IMQ + S3)—(a): 625 nodes; (b): 1089 nodes.
Symmetry 14 01476 g012
Figure 13. Example 3: convergence study—(a): convergence curves of ru; (b): convergence curves of r σ .
Figure 13. Example 3: convergence study—(a): convergence curves of ru; (b): convergence curves of r σ .
Symmetry 14 01476 g013
Figure 14. Example 3: stress computational results using HRBF (IMQ + S3, 1089 nodes)—(a): analytical solution of normal stress σ x x ; (b): numerical solution of normal stress σ x x ; (c): analytical solution of normal stress σ y y ; (d): numerical solution of normal stress σ y y ; (e): analytical solution of shear stress τ x y ; (f): numerical solution of shear stress τ x y .
Figure 14. Example 3: stress computational results using HRBF (IMQ + S3, 1089 nodes)—(a): analytical solution of normal stress σ x x ; (b): numerical solution of normal stress σ x x ; (c): analytical solution of normal stress σ y y ; (d): numerical solution of normal stress σ y y ; (e): analytical solution of shear stress τ x y ; (f): numerical solution of shear stress τ x y .
Symmetry 14 01476 g014
Table 1. Some typical RBFs.
Table 1. Some typical RBFs.
Infinitely Smooth RBFsPiecewise Smooth RBFs
NameExpressionNameExpression
Multiquadric (MQ) c 2 + r j 2 Splines of degree n (Sn) r j 2 n 1 , n
Inverse multiquadric (IMQ) 1 / c 2 + r j 2 Thin plate spline (TPS) r j 2 n log ( r j )
Gaussian (GA) e c 2 r j 2 Linear r j
Note: r j = x x j 2 , c represents the shape parameter.
Table 2. The hybrid RBFs.
Table 2. The hybrid RBFs.
NameExpression
MQ + TPS r j 2 + c 2 + w r j 4 ln ( r j )
IMQ + TPS 1 / r j 2 + c 2 + w r j 4 ln ( r j )
GA + TPS e c 2 r j 2 + w r j 4 ln ( r j )
MQ + S3 r j 2 + c 2 + w r j 3
IMQ + S3 1 / r j 2 + c 2 + w r j 3
GA + S3 e c 2 r j 2 + w r j 3
Note: r j = x x j 2 , c represents the shape parameter, and w represents the weight parameter.
Table 3. Example 1: computational results of different methods (MQ and its variants).
Table 3. Example 1: computational results of different methods (MQ and its variants).
Internal NodeCoordinatesMethod
ExactMQMQ + TPSMQ + S3
5(1,1)
(ux, uy)(0.600,0.600)(0.593,0.589)(0.593,0.589)(0.593,0.590)
( σ x x , σ y y , τ x y ) (0.857,0.857,0)(0.896,0.890,0.010)(0.896,0.890,0.010)(0.895,0.889,0.010)
10(0.65,1)
(ux, uy)(0.390,0.600)(0.374,0.588)(0.374,0.588)(0.374,0.588)
( σ x x , σ y y , τ x y ) (0.857,0.857,0)(0.893,0.900,0.009)(0.893,0.900,0.009)(0.892,0.899,0.009)
11(0.7,1.5)
(ux, uy)(0.420,0.900)(0.409,0.901)(0.409,0.901)(0.409,0.901)
( σ x x , σ y y , τ x y ) (0.857,0.857,0)(0.906,0.893,−0.006)(0.906,0.893,−0.006)(0.904,0.892,−0.006)
12(1.3,1.2)
(ux, uy)(0.780,0.720)(0.784,0.718)(0.784,0.718)(0.784,0.718)
( σ x x , σ y y , τ x y ) (0.857,0.857,0)(0.892,0.896,0.012)(0.892,0.896,0.012)(0.891,0.895,0.011)
13(1.2,0.6)
(ux, uy)(0.720,0.360)(0.718,0.343)(0.718,0.343)(0.718,0.343)
( σ x x , σ y y , τ x y ) (0.857,0.857,0)(0.909,0.889,−0.003)(0.909,0.889,−0.003)(0.908,0.888,−0.003)
Table 4. Example 2: computational results of different methods for the irregular node distribution.
Table 4. Example 2: computational results of different methods for the irregular node distribution.
Method N = 39N = 125N = 441N = 949
MQcopt201031.6
ru2.5513 × 10−67.1523 × 10−72.5843 × 10−73.7019 × 10−7
r σ 1.4577 × 10−38.6395 × 10−51.7324 × 10−43.3451 × 10−4
CPU (s)0.04450.04510.20920.7934
IMQcopt18931.7
ru2.6001 × 10−69.4806 × 10−79.5554 × 10−72.6342 × 10−6
r σ 1.9224 × 10−36.2694 × 10−46.7735 × 10−42.4746 × 10−3
CPU (s)0.02700.03810.22880.9858
GAcopt14.56.594.5
ru2.7955 × 10−64.1935 × 10−75.3383 × 10−63.8875 × 10−6
r σ 1.2317 × 10−35.5071 × 10−52.3780 × 10−41.9450 × 10−4
CPU (s)0.00710.02690.11280.4085
MQ + TPScopt(19.5,10−15)(14,10−13)(12.5,10−11)(11.5,10−10)
ru3.3220 × 10−62.2796 × 10−72.1210 × 10−81.3474 × 10−8
r σ 2.8319 × 10−37.7843 × 10−59.3717 × 10−68.6273 × 10−6
CPU (s)0.01140.07230.22711.0846
IMQ + TPScopt(19,10−15)(13,10−15)(16,10−13)(11.5,10−12)
ru2.2627 × 10−63.4743 × 10−74.7540 × 10−82.9863 × 10−8
r σ 1.6277 × 10−31.6833 × 10−42.8371 × 10−52.0773 × 10−5
CPU (s)0.01080.03960.27621.3265
GA + TPScopt(19.5,10−15)(9.5,10−14)(10,10−12)(9,10−11)
ru2.7870 × 10−61.7480 × 10−73.5633 × 10−82.7927 × 10−8
r σ 1.5181 × 10−38.5198 × 10−52.1890 × 10−52.1679 × 10−5
CPU (s)0.00930.06860.20680.9901
MQ + S3copt(20,10−12)(13,10−15)(10,10−11)(4.5,10−10)
ru5.1602 × 10−63.3904 × 10−74.5543 × 10−85.0247 × 10−8
r σ 2.8665 × 10−33.0302 × 10−52.1685 × 10−52.5899 × 10−5
CPU (s)0.01980.05850.19820.8878
IMQ + S3copt(20,10−15)(13,10−15)(14,10−13)(15.5,10−13)
ru2.6737 × 10−61.3149 × 10−71.9911 × 10−81.1097 × 10−8
r σ 1.7319 × 10−35.7103 × 10−51.0639 × 10−55.6919 × 10−6
CPU (s)0.01270.06540.24891.1624
GA + S3copt(20,10−12)(11,10−14)(10.5,10−12)(9,10−11)
ru3.0658 × 10−61.1867 × 10−71.3768 × 10−81.3672 × 10−8
r σ 2.9404 × 10−33.4692 × 10−59.9027 × 10−61.0190 × 10−5
CPU (s)0.01180.05730.18180.7992
TPSru2.0007 × 10−47.5755 × 10−51.6952 × 10−58.0060 × 10−6
r σ 1.3620 × 10−12.7900 × 10−26.6000 × 10−32.9000 × 10−3
CPU (s)0.00510.02970.13140.5920
S3ru3.3979 × 10−41.8154 × 10−45.0393 × 10−52.7450 × 10−5
r σ 2.0560 × 10−18.5200 × 10−23.0800 × 10−21.7500 × 10−2
CPU (s)0.02790.02240.09780.4239
Note: The bold numbers represent the best calculated data.
Table 5. Example 2: computational results of different methods for regular node distribution.
Table 5. Example 2: computational results of different methods for regular node distribution.
Method N =39 (3 × 13)N =125 (5 × 25)N =441 (9 × 49)N =949 (13 × 73)
MQcopt18143.51.7
ru8.0866 × 10−62.9725 × 10−64.0491 × 10−75.3096 × 10−7
r σ 8.3123 × 10−26.4158 × 10−41.2788 × 10−43.8119 × 10−4
CPU (s)0.03380.04370.18690.7948
IMQcopt18.516.53.52.1
ru1.2985 × 10−52.1955 × 10−61.4602 × 10−61.1814 × 10−6
r σ 8.0496 × 10−26.9578 × 10−46.7189 × 10−47.5492 × 10−4
CPU (s)0.08520.03900.21820.9871
GAcopt181346
ru7.3142 × 10−62.5473 × 10−63.7281 × 10−65.8424 × 10−6
r σ 6.1412 × 10−23.6249 × 10−42.1470 × 10−42.2683 × 10−4
CPU (s)0.00730.02370.21140.4151
MQ + TPScopt(20,10−15)(20,10−10)(11,10−11)(10.5,10−10)
ru7.1497 × 10−67.4186 × 10−71.4188 × 10−89.3300 × 10−9
r σ 8.5075 × 10−22.9660 × 10−48.1700 × 10−66.4800 × 10−6
CPU (s)0.04670.05720.23211.0553
IMQ + TPScopt(20,10−15)(20,10−14)(11.5,10−13)(10.5,10−12)
ru1.2108 × 10−51.2125 × 10−62.9799 × 10−81.5791 × 10−8
r σ 7.6433 × 10−24.5468 × 10−42.1833 × 10−51.5639 × 10−5
CPU (s)0.04440.06960.28021.3881
GA + TPScopt(20,10−10)(20,10−13)(10,10−12)(9,10−11)
ru3.4546 × 10−63.7997 × 10−72.6746 × 10−82.0292 × 10−8
r σ 5.8012 × 10−21.5217 × 10−41.9862 × 10−51.6821 × 10−5
CPU (s)0.00690.03190.21910.9833
MQ + S3copt(1910−12)(20,10−12)(11,10−11)(9,10−10)
ru5.3438 × 10−67.6381 × 10−73.1988 × 10−82.9318 × 10−8
r σ 7.5212 × 10−22.7886 × 10−41.8097 × 10−53.0168 × 10−5
CPU (s)0.03980.03290.20590.8991
IMQ + S3copt(20,10−15)(20,10−13)(14.5,10−13)(16,10−13)
ru1.3095 × 10−51.1944 × 10−61.2919 × 10−86.9444 × 10−9
r σ 7.7167 × 10−24.3883 × 10−48.8861 × 10−64.4397 × 10−6
CPU (s)0.03090.03720.25851.1514
GA + S3copt(20,10−12)(20,10−13)(11,10−12)(9.5,10−11)
ru4.5154 × 10−63.4575 × 10−71.2244 × 10−89.7053 × 10−9
r σ 5.8101 × 10−21.2654 × 10−48.6607 × 10−68.6815 × 10−6
CPU (s)0.00530.02970.19340.7901
TPSru1.7521 × 10−46.2527 × 10−51.0582 × 10−53.4544 × 10−6
r σ 1.3640 × 10−12.7100 × 10−26.0000 × 10−32.4000 × 10−3
CPU (s)0.00420.02900.13490.5985
S3ru3.3830 × 10−41.7256 × 10−44.1090 × 10−51.6038 × 10−5
r σ 2.0250 × 10−11.3640 × 10−12.7400 × 10−21.3900 × 10−2
CPU (s)0.00340.02710.10250.4247
Note: The bold numbers represent the best calculated data.
Table 6. Example 3: computational results of different methods.
Table 6. Example 3: computational results of different methods.
Method N = 289N = 625N = 1089
MQcopt0.90.80.8
ru4.2960 × 10−41.7965 × 10−47.1667 × 10−5
r σ 1.3975 × 10−25.4578 × 10−32.4684 × 10−3
CPU (s)0.1746 0.35311.1068
IMQcopt1.31.21.2
ru3.1431 × 10−31.0973 × 10−36.0918 × 10−4
r σ 6.8418 × 10−23.0396 × 10−22.0173 × 10−2
CPU (s)0.16580.20701.3669
GAcopt1.822.1
ru5.6372 × 10−22.8088 × 10−21.0097 × 10−2
r σ 9.4539 × 10−16.0751 × 10−14.9868 × 10−1
CPU (s)0.05990.10480.5760
MQ + TPScopt(1.7,10−5)(1.6,10−7)(1.4,10−8)
ru1.8100 × 10−43.1703 × 10−57.6593 × 10−6
r σ 1.3936 × 10−22.2995 × 10−34.9781 × 10−4
CPU (s)0.15350.47091.5308
IMQ + TPScopt(20,10−7)(20,10−7)(20,10−7)
ru6.2954 × 10−42.2173 × 10−41.0269 × 10−4
r σ 3.4004 × 10−21.5125 × 10−28.2553 × 10−3
CPU (s)0.13710.58291.8389
GA + TPScopt(0.01,10−9)(0.01,10−9)(0.01,10−7)
ru5.9965 × 10−42.1045 × 10−49.9975 × 10−5
r σ 3.3397 × 10−21.4878 × 10−28.1967 × 10−3
CPU (s)0.09750.42651.3619
MQ + S3copt(0.9,10−8)(0.9,10−12)(0.8,10−9)
ru4.1938 × 10−49.7755 × 10−56.4772 × 10−5
r σ 1.4432 × 10−24.5325 × 10−32.2693 × 10−3
CPU (s)0.12640.41271.2019
IMQ + S3copt(2.5,10−6)(2.5,10−8)(2.2,10−9)
ru4.6816 × 10−49.0696 × 10−51.8502 × 10−5
r σ 2.6994 × 10−27.4137 × 10−31.8960 × 10−3
CPU (s)0.14950.53711.6242
GA+S3copt(0.4,10−5)(0.8,10−5)(0.8,10−5)
ru1.0062 × 10−33.1934 × 10−41.2175 × 10−4
r σ 4.9156 × 10−21.6263 × 10−29.4336 × 10−3
CPU (s)0.11730.37051.1097
TPSru7.6646 × 10−42.6404 × 10−41.2074 × 10−4
r σ 3.6900 × 10−21.6100 × 10−28.7000 × 10−3
CPU (s)0.12150.27060.7952
S3ru1.2000 × 10−34.9997 × 10−42.5632 × 10−4
r σ 5.5700 × 10−22.8800 × 10−21.7600 × 10−2
CPU (s)0.12090.20320.5977
Note: The bold numbers represent the best calculated data.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Chen, Y.-T.; Li, C.; Yao, L.-Q.; Cao, Y. A Hybrid RBF Collocation Method and Its Application in the Elastostatic Symmetric Problems. Symmetry 2022, 14, 1476. https://doi.org/10.3390/sym14071476

AMA Style

Chen Y-T, Li C, Yao L-Q, Cao Y. A Hybrid RBF Collocation Method and Its Application in the Elastostatic Symmetric Problems. Symmetry. 2022; 14(7):1476. https://doi.org/10.3390/sym14071476

Chicago/Turabian Style

Chen, Ying-Ting, Cheng Li, Lin-Quan Yao, and Yang Cao. 2022. "A Hybrid RBF Collocation Method and Its Application in the Elastostatic Symmetric Problems" Symmetry 14, no. 7: 1476. https://doi.org/10.3390/sym14071476

APA Style

Chen, Y. -T., Li, C., Yao, L. -Q., & Cao, Y. (2022). A Hybrid RBF Collocation Method and Its Application in the Elastostatic Symmetric Problems. Symmetry, 14(7), 1476. https://doi.org/10.3390/sym14071476

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