Radial Point Interpolation Method Research Articles

Nonlinear analysis of three-dimensional hyperelastic problems using radial point interpolation method

Hyperelastic materials are primarily common in real life as well as in industry applications, and studying this kind of material is still an active research area. Naturally, the characteristic of hyperelastic material will be expressed when it undergoes large deformation, so the geometrical nonlinear effect should be considered. To analyze the behavior of hyperelastic material, the Neo-Hookean model is imposed in this study because of its simplicity. The model shows the nonlinear behavior when the deformation becomes large due to the nonlinear displacement-strain relation. This constitutive relation also gives a good correlation with experimental data. This study performs a meshless method, namely the radial point interpolation method (RPIM), to analyze the nonlinear behavior of Neo-Hookean hyperelastic material under a finite deformation state in three-dimensional space. The standard Newton-Raphson technique is applied to obtain the nonlinear solutions. Unlike mesh-based approaches, the meshless method shows its advantages in large deformation problems due to its mesh independence. In this paper, the numerical results of three-dimensional problems that undergo large deformation will be calculated and validated with solutions derived from previous studies. By investigating the obtained results, the superior ability of RPIM in hyperelastic problems can be proved.

Meshfree method for bending and free vibration analysis of laminated plates using the Reissner’s mixed variational theorem

Free vibration and bending behavior of laminated plates are investigated using a combination of Reissner’s mixed variational theorem and the meshless radial point interpolation method. By introducing Reissner’s mixed variational theorem, a refined transverse shear stress field is derived based on three-dimensional (3D) elasticity equations. The efficiency and accuracy of this model are demonstrated through numerical examples. It is found that accurate prediction of shear stresses significantly improves the prediction of natural frequency for composite laminated plates. The difference between shear stresses derived from 3D equilibrium equations and that from constitutive equations is significant.

Extending the Natural Neighbour Radial Point Interpolation Meshless Method to the Multiscale Analysis of Sandwich Beams with Polyurethane Foam Core

This work investigates the mechanical behaviour of sandwich beams with cellular cores using a multiscale approach combined with a meshless method, the Natural Neighbour Radial Point Interpolation Method (NNRPIM). The analysis is divided into two steps, aiming to analyse the efficiency of NNRPIM formulation when combined with homogenisation techniques for a multiscale computational framework of large-scale sandwich beam problems. In the first step, the cellular core material undergoes a controlled modification process in which circular holes are introduced into bulk polyurethane foam (PUF) to create materials with varying volume fractions. Subsequently, a homogenisation technique is combined with NNRPIM to determine the homogenised mechanical properties of these PUF materials with different porosities. In this step, NNRPIM solutions are compared with high-order FEM simulations. While the results demonstrate that RPIM can approximate high-order FEM solutions, it is observed that the computational cost increases significantly when aiming for comparable smoothness in the approximations. The second step applies the homogenised mechanical properties obtained in the first step to analyse large-scale sandwich beam problems with both homogeneous and functionally graded cores. The results reveal the capability of NNRPIM to closely replicate the solutions obtained from FEM analyses. Furthermore, an analysis of stress distributions along the beam thickness highlights a tendency for some NNRPIM formulations to yield slightly lower stress values near the domain boundaries. However, convergence towards agreement among different formulations is observed with mesh refinement. The findings of this study show that NNRPIM can be used as an alternative numerical method to FEM for analysing sandwich structures.

A Meshless Radial Point Interpolation Method for Solving Fractional Navier–Stokes Equations

This paper aims to develop a meshless radial point interpolation (RPI) method for obtaining the numerical solution of fractional Navier–Stokes equations. The proposed RPI method discretizes differential equations into highly nonlinear algebraic equations, which are subsequently solved using a fixed-point method. Furthermore, a comprehensive analysis regarding the effects of spatial and temporal discretization, polynomial order, and fractional order is conducted. These factors’ impacts on the accuracy and efficiency of the solutions are discussed in detail. It can be shown that the meshless RPI method works quite well for solving some benchmark problems accurately.

Enhanced meshfree method with nodal integration for analysis of functionally graded material sandwich curved shells

This paper presents a nodal integration technique, the sub-domain stabilized conforming integration (SSCI), for the meshfree radial point interpolation method (RPIM) applied to the static and modal analysis of functionally graded material (FGM) sandwich curved shells. FGM sandwich shells with different kinds of core and face sheets are considered in this work while the interested curved shell is formulated by the first-order shear deformation theory. The numerical integration technique to compute the stiffness and mass matrices in the equilibrium equation is the SSCI, which is a stabilized nodal integration with strain smoothing to preserve the accuracy and stability of the numerical results. The RPIM shape functions are utilized in this study for interpolating both the field variables and the geometry of the curved shell due to their ability to satisfy the Kronecker delta property, a rare advantage among meshfree methods. The static and modal analysis of different geometry curved shells with various sandwich FGMs are conducted. Through several numerical examples, the accuracy and efficiency of the SSCI technique in the meshfree RPIM are demonstrated and discussed.

Deformation monitoring for fixed-wing UAS through Inverse Mesh-free Method

This study proposes the Inverse Local Radial Point Interpolation Method (iLRPIM) for reconstructing the deformation of irregular plate structures. The proposed methodology eliminates the drawbacks of the existing inverse finite-element method (iFEM), where the Jacobian matrix results in reduced computational accuracy for irregular plate-shape sensing. To this end, this paper released from the element and used the Mesh-Free method to describe the calculated strains. Subsequently, the system equation was formulated by minimizing the error between the measured and calculated strains. The deformation parameters were then computed by solving this system equation. Additionally, compared to iFEM, iLRPIM allows for easier adjustment of the interpolation function, making it more suitable for real-world engineering applications. The effectiveness of iLRPIM has been validated through simulations and practical experiments using four typical models. In simulations of the UAS wing box equivalent model, peak errors were below 0.8 mm for a deformation range of 25 mm. In practical scenarios, peak errors were below 3.1 mm for a deformation range of 85 mm. These results underscore the system’s capability to accurately sense the shapes of irregular plate structures using iLRPIM.

An in-depth analysis of the Radial Point Interpolation Method’s parameters in relation to the accuracy of fracture problem

The fracture problem has been modeled using Radial Point Interpolation Method (RPIM) to examine the effects of support domain size and shape parameter value on the calculation of Stress Intensity Factor (SIF). The study focuses on an edge crack problem subjected to normal load and shear loads, determining that the support domain size of ‘ α ≥ 3.5 and shape parameter value between ‘ n = 1.9 -2.1’ yield results of accuracy with error % ≤ 3 % . Additionally, research extended by varying the problem domain size and crack location to verify consistency. The findings indicate that selecting the optimal shape parameter value will ensure better accuracy over the enrichment techniques with lesser computational cost.

Multiscale Analysis of Sandwich Beams with Polyurethane Foam Core: A Comparative Study of Finite Element Methods and Radial Point Interpolation Method.

This study presents a comprehensive multiscale analysis of sandwich beams with a polyurethane foam (PUF) core, delivering a numerical comparison between finite element methods (FEMs) and a meshless method: the radial point interpolation method (RPIM). This work aims to combine RPIM with homogenisation techniques for multiscale analysis, being divided in two phases. In the first phase, bulk PUF material was modified by incorporating circular holes to create PUFs with varying volume fractions. Then, using a homogenisation technique coupled with FEM and four versions of RPIM, the homogenised mechanical properties of distinct PUF with different volume fractions were determined. It was observed that RPIM formulations, with higher-order integration schemes, are capable of approximating the solution and field smoothness of high-order FEM formulations. However, seeking a comparable field smoothness represents prohibitive computational costs for RPIM formulations. In a second phase, the obtained homogenised mechanical properties were applied to large-scale sandwich beam problems with homogeneous and approximately functionally graded cores, showing RPIM's capability to closely approximate FEM results. The analysis of stress distributions along the thickness of the beam highlighted RPIM's tendency to yield lower stress values near domain edges, albeit with convergence towards agreement among different formulations. It was found that RPIM formulations with lower nodal connectivity are very efficient, balancing computational cost and accuracy. Overall, this study shows RPIM's viability as an alternative to FEM for addressing practical elasticity applications.

Buckling analysis of functionally graded sandwich thin plates using a meshfree Hermite Radial Point Interpolation Method

Buckling analysis of functionally graded sandwich thin plates using a meshfree Hermite Radial Point Interpolation Method

Three-dimensional elastoplastic post-buckling analysis of functionally graded plates using a novel meshfree Tchebychev-radial point interpolation approach

This paper presents a three-dimensional (3D) analysis of the post-buckling behavior of functionally graded (FG) plates for the first time using an elastoplasticity-based meshfree formulation. A novel 3D Tchebychev-radial point interpolation approach, which combines radial basis functions with Tchebychev polynomials, is developed to solve the nonlinear post-buckling problem. The incremental plastic deformation is modeled by the Prandtl–Reuss flow rule along the isotropic hardening von Mises criterion. The governing equations are derived using the principle of virtual work by considering the 3D full Green–Lagrange strain components. The Newton–Raphson technique along with the arc-length method is used to determine post-buckling equilibrium paths of FG plates. The effective elastoplastic properties of the functionally graded material are assessed by exploiting the homogenization method, named Tamura–Tomota–Ozawa (TTO) model. It has been demonstrated that the accuracy of the solution is improved by incorporating Tchebychev polynomials into the radial point interpolation method (RPIM) shape function, and the stability and robustness of the results are independent of variations in the shape parameter. Furthermore, it is confirmed that TRPIM, with a slightly higher CPU time compared to RPIM, exhibits a higher convergence rate. The excellent agreement of the results with those existing in the literature shows that the proposed meshfree method can be used as a robust and accurate numerical tool to predict the elastoplastic post-buckling behavior of FG plates. Further numerical assessments indicate that post-buckling paths are significantly affected by factors such as geometric parameters, material gradient, loading ratio, and boundary conditions (BCs).

A robust radial point interpolation method empowered with neural network solvers (RPIM-NNS) for nonlinear solid mechanics

In this work, we proposed a robust radial point interpolation method empowered with neural network solvers (RPIM-NNS) for solving highly nonlinear solid mechanics problems. It is enabled by neural network solvers via minimizing an energy-based functional loss. The RPIM-NNS has the following key ingredients: (1) It uses radial basis functions (RBFs) for displacement interpolation at arbitrary points in the problem domain, permitting irregular node distributions. (2) Nodes are placed also beyond the domain boundary, allowing the convenient implementation of boundary conditions of both Dirichlet and Neumann types. (3) It uses strain energy in an integral form as a part of the loss function, ensuring solution stability. (4) A well-developed gradient descendant algorithm in machine learning is employed to find the optimal solution, enabling robustness and ease in handling material and geometrical nonlinearities. (5) The proposed RPIM-NNS is compatible with parallel computing schemes. The performance of this method is tested using nonlinear problems including Cook's membrane and 3D twisting rubber problems, demonstrating its remarkable stability and robustness. This work, which seamlessly integrates the neural network solvers with mechanics governing equations and computational mechanics techniques, offers an excellent alternative for nonlinear solid mechanics problems. MATLAB codes are made available at https://github.com/JinshuaiBai/RPIM_NNS for free downloading.

Predicting the mixed-mode fracture propagation in single-leg bending using a radial point interpolation meshless method

Adhesive bonding is a widely used joining technique with increasing applications in the automotive and aircraft industries. However, the numerical simulation of adhesively bonded joints is challenging due to the intricate mechanical behavior of the bonding layer, especially when structures are subjected to mixed-mode loading. This work proposes a meshless fracture propagation algorithm to simulate mixed-mode fracture propagation in single-leg bending (SLB) adhesive joints. The proposed algorithm combines the radial point interpolation method (RPIM) with a mixed-mode strain-based criterion to predict crack growth. The RPIM permits a flexible discretization of the fracture region and eases the geometric remeshing of the integration cells to account for the crack tip propagation. Additionally, the RPIM provides accurate and smooth strain/stress fields, allowing implementation of the mixed-mode fracture criterion. The numerical model was validated against experimental data for SLB adhesive joints with a brittle adhesive. The results showed that the proposed algorithm can accurately predict the resistance curves and critical energy release rates of the adhesive joints.

Fractional Boundary Element Solution for Nonlinear Nonlocal Thermoelastic Problems of Anisotropic Fibrous Polymer Nanomaterials

This paper provides a new fractional boundary element method (BEM) solution for nonlinear nonlocal thermoelastic problems with anisotropic fibrous polymer nanoparticles. This comprehensive BEM solution comprises two solutions: the anisotropic fibrous polymer nanoparticles problem solution and the nonlinear nonlocal thermoelasticity problem. The nonlinear nonlocal thermoelasticity problem solution separates the displacement field into complimentary and specific components. The overall displacement is obtained using the boundary element methodology, which solves a Navier-type problem, and the specific displacement is derived using the local radial point interpolation method (LRPIM). The new modified shift-splitting (NMSS) technique, which minimizes memory and processing time requirements, was utilized to solve BEM-created linear systems. The performance of NMSS was evaluated. The numerical results show how fractional and graded parameters influence the thermal stresses of nonlinear nonlocal thermoelastic issues involving anisotropic fibrous polymer nanoparticles. The numerical findings further reveal that the BEM results correlate very well with the finite element method (FEM) and analytical results, demonstrating the validity and correctness of the proposed methodology.

A meshless computational framework for a modified dynamic system of vehicle coupled with plate structure

In previous simulations of train–bridge interaction systems (TBIS), the supporting system for the train are commonly treated as beam structures, leading to less accurate results, particularly for small-span cases. To address this limitation, a modified vertical TBIS is proposed. In the presented TBIS, the supporting system is considered as a Reissner–Mindlin plate, and the displacement field is described by first-order shear deformation theory (FSDT). To establish the model, radial point interpolation method (RPIM), a meshless method, is employed. Finally, a coupled dynamic equation is established to calculate various responses of the system. Several numerical examples are presented to illustrate the disparities between the system based on plate model and traditional beam model. The results indicate that the beam model yields higher estimates of the mid-span vertical displacement of the bridge, while the peak of the mid-span vertical acceleration is smaller compared to the plate model; additionally, it is observed that the carbody is primarily influenced by rail irregularities. Consequently, the proposed plate model offers distinct advantages over the beam model in providing comprehensive structural response information, thereby offering novel insights into bridge design and analysis. Additionally, this marks the inaugural application of the meshless method in the field of TBIS, which further extends the application scope of meshless methods.

Radial basis point interpolation for strain field calculation in digital image correlation

In order to extract smooth and accurate strain fields from the noisy displacement fields obtained by digital image correlation (DIC), a point interpolation meshless (PIM) method with a radial basis function (RBF) is introduced for full-field strain calculation, which overcomes the problems of slow calculation speed and unstable matrix inverse calculation of the element-free Galerkin method (EFG). The radial basis point interpolation method (RPIM) with three different radial basis functions and the moving least squares (MLS) and pointwise least squares (PLS) methods are compared by analyzing and validating the strain fields with high-strain gradients in simulation experiments. The results indicate that the RPIM is nearly 80% more computationally efficient than the MLS method when a larger support domain is used, and the efficiency of the RPIM is nearly 26% higher than that of the MLS method when a smaller support domain is used; the strain calculation accuracy is slightly lower than that of the MLS method by 0.3–0.5%, but the stability of the calculation is significantly improved. In contrast with the PLS method, which is easily affected by the noise and the size of the strain calculation window, the RPIM is insensitive to the displacement noise and the size of the support domain and can obtain a similar calculation accuracy. The RPIM with multiquadric (MQ) radial basis functions performs well in balancing the computational accuracy and efficiency and is insensitive to shape parameters. The application cases show that the method can effectively compute the strain field at the crack tip, validating its applicability to the study of the plastic region at the crack tip. In conclusion, the proposed RPIM-based method provides an accurate, practical, and robust approach for full-field strain measurements.

Dynamic Analysis of an FGM Beam Undergoing Large Overall Motion in Temperature Field

The dynamic analysis of a functionally gradient material (FGM) beam undergoing large overall motion based on meshless methods is studied in the variable temperature field environment. Considering the high-order terms of nonlinear coupling deformation which were often ignored in previous studies, the high-order rigid–flexible coupled (HOC) dynamic model of a rotating FGM beam with a lumped mass in a temperature field is established by employing Lagrange’s equations of the second kind. Compared with the first-order approximate coupled (FOAC) dynamic model, the HOC dynamic model is more accurate and has a wider scope of application. The validity of the point interpolation (PIM) method and radial point interpolation method (RPIM) in this paper is verified by comparison with simulation results of the assumed mode method (AMM) and the finite element method (FEM). On this basis, the factors affecting the dynamic characteristics of the FGM beam such as the form of temperature field, the functionally gradient index, and the location of added lumped mass are discussed. In this paper, the meshless method with higher computational efficiency and the HOC dynamic model with more accurate results are adopted. This study will provide guidance and reference for dealing with problems in engineering, machinery, aerospace, and other fields in complex environment.

Nonlinear dynamic analysis of aircraft CFRP sandwich wings under explosive blast loading: Introducing SVM-DNN algorithm to predict dynamical information

The aircraft wings are a crucial component, particularly during mechanical shock excitation. Thus, for the first time, an inventive adaptable model to simulate the nonlinear dynamics of the wings under explosive blast loading at the moment of the collision is described in this article. It is crucial to increase the structure's stability because of explosive blast loading. The mentioned structure is made of two composite face sheets and a carbon fiber-reinforced polymer (CFRP). Functionally graded (FG) graphene origami (GOri)-enabled auxetic metallic metamaterial (GOEAM) is introduced to improve the stability of this kind of applicable structure, especially in air maneuvers. The analysis is conducted utilizing Reddy's third-order shear deformation of Reddy theory (TSDRT), incorporating Hamilton's principle and Von Kármán nonlinear geometric assumptions. The governing and boundary equations are accurately computed at the domain and boundary edges of the doubly curved structure. The coupled mesh-free radial point interpolation method (MRPIM) and Newton-Raphson method are applied for numerically simulating the current doubly curved structure by employing Newmark's temporal integration method with a constant-average acceleration approach. The present findings are cross-checked against the machine learning method and open-source results from the literature. Based on physical data, machine learning integrates supervised machine learning to predict the nonlinear vibration of the system. In this case, nonlinear transient bending properties are found using hybrid machine learning methods and solutions. The findings demonstrate the significance of GOEAM for nonlinear vibrations of the composite system. The current study's conclusions could be used as a benchmark and useful recommendations for future structural design techniques with enhanced features.

Dynamic analysis of sandwich conical shell with magneto/electro rheological core using a truly meshfree radial point interpolation method

Dynamic analysis of sandwich conical shell with magneto/electro rheological core using a truly meshfree radial point interpolation method

Transient analyses of underwater acoustic propagation with the modified meshfree radial point interpolation method and newmark time integration techniques

Underwater acoustic propagation plays a very important role in underwater communication and ocean engineering field. The radial point interpolation method (RPIM) is an effective numerical approach to solve underwater acoustic propagations though considerable numerical errors can always be caused due to the discontinuously differentiable nodal interpolation function in one integration cell. The origin of this numerical error primarily comes from the mismatch between the background integration cell and the nodal interpolation function support. To address this issue, a modified meshfree RPIM is developed in this work for transient analyses of underwater acoustic propagations. In the proposed modified RPIM, the construction of the interpolation function supports are dependent on the integration cells. In consequence, in one background cell all the quadrature points share one mutual support domain and the continuously differentiable interpolation functions can be obtained. Therefore, the numerical integration errors can be evidently decreased and the related numerical dispersion effects in wave analysis can be largely suppressed. Owing to the adequately small spatial dispersion errors from the proposed modified RPIM, it is observed that the modified RPIM with the Newmark time domain discretization scheme can monotonically improve the solution precision by directly reducing the temporal discretization interval. This property makes the modified RPIM clearly outperform the original RPIM and is especially suitable for handling very complex underwater acoustic propagations.

A novel train–bridge interaction computational framework based on a meshless box girder model

In traditional train–bridge coupled system (TBCS), simply supported box girder bridges are often modeled using Euler beam elements, neglecting their spatial structure. This simplification may yield inaccurate results, impacting the running safety analysis of high-speed railway. To address this issue, a novel train–bridge interaction computational framework based on first-order shear deformation theory (FSDT) and radial point interpolation method (RPIM) for the box girder bridge model is proposed. In this model, the displacement fields of top, bottom, and web plates are represented using FSDT and numerically discretized by RPIM. The traditional TBCS is upgraded by replacing the Euler beam model with the novel model. This is the first time that the framework has been applied to TBCS field. Several numerical examples are presented to highlight the accuracy of this novel model, and illustrate the differences, and advantages of it over the traditional model. The results indicate that the proposed model closely matches the accuracy of the solid element (C3D20) model; it can provide a comprehensive bridge response compared to traditional TBCS, and the latter may underestimate the dynamic response of the structure. The proposed model holds significant potential for the simulation of box structures and widespread application in the field of TBCS.