Search results

1) The OE-MLPG penalty meshfree method is developed to solve cracked structure.(2) Smartening the numerical meshfree method by combining the particle swarm optimization (PSO) optimization algorithms and Voronoi computational geometric algorithm. (3). Selection of base functions, finding optimal penalty factor and distribution of appropriate nodal points to the accuracy of calculation in the meshless local Petrov–Galekrin (MLPG) meshless method.

Using appropriate shape functions and distribution of nodal points in local domains and sub-domains and choosing an approximation or interpolation method has an effective role in the application of meshless methods for the analysis of computational fracture mechanics problems, especially problems with geometric discontinuity and cracks. In this research, computational geometry technique, based on the Voronoi diagram (VD) and Delaunay triangulation and PSO algorithm, are used to distribute nodal points in the sub-domain of analysis (crack line and around it on the crack plane).

By doing this process, the problems caused by too closeness of nodal points in computationally sensitive areas that exist in general methods of nodal point distribution are also solved. Comparing the effect of the number of sentences of basic functions and their order in the definition of shape functions, performing the mono-objective PSO algorithm to find the penalty factor, the coefficient, convergence, arrangement of nodal points during the three stages of VD implementation and the accuracy of the answers found indicates, the efficiency of V-E-MLPG method with Ns = 7 and ß = 0.0037–0.0075 to estimation of 3D-stress intensity factors (3D-SIFs) in computational fracture mechanics.

The present manuscript is a continuation of the studies (Ref. [33]) carried out by the authors, about; feasibility assessment, improvement and solution of challenges, introduction of more capacities and capabilities of the numerical MLPG method have been used. In order to validate the modeling and accuracy of calculations, the results have been compared with the findings of reference article [34] and [35].

This article offered two meshfree algorithms, namely the local radial basis functions-finite difference (LRBF-FD) approximation and local radial basis functions-differential quadrature method (LRBF-DQM) to simulate the multidimensional hyperbolic wave models and work is an extension of Jiwari (2015).

In the evolvement of the first algorithm, the time derivative is discretized by the forward FD scheme and the Crank-Nicolson scheme is used for the rest of the terms. After that, the LRBF-FD approximation is used for spatial discretization and quasi-linearization process for linearization of the problem. Finally, the obtained linear system is solved by the LU decomposition method. In the development of the second algorithm, semi-discretization in space is done via LRBF-DQM and then an explicit RK4 is used for fully discretization in time.

For simulation purposes, some 1D and 2D wave models are pondered to instigate the chastity and competence of the developed algorithms.

The developed algorithms are novel for the multidimensional hyperbolic wave models. Also, the stability analysis of the second algorithm is a new work for these types of model.

In this paper, the smoothed point interpolation method (SPIM) is used to model the slope deformation. However, the computational efficiency of SPIM is not satisfying when modeling the large-scale nonlinear deformation problems of geological bodies.

In this paper, the SPIM is used to model the slope deformation. However, the computational efficiency of SPIM is not satisfying when modeling the large-scale nonlinear deformation problems of geological bodies.

A simple slope model with different mesh sizes is used to verify the performance of the efficient face-based SPIM. The first accelerating strategy greatly enhances the computational efficiency of solving the large-scale slope deformation. The second accelerating strategy effectively improves the convergence of nonlinear behavior that occurred in the slope deformation.

The designed efficient face-based SPIM can enhance the computational efficiency when analyzing large-scale nonlinear slope deformation problems, which can help to predict and prevent potential geological hazards.

Fused Filament Fabrication (FFF) is an extrusion-based manufacturing process using fused thermoplastics. Despite its low cost, the FFF is not extensively used in high-value industrial sectors mainly due to parts' anisotropy (related to the deposition strategy) and residual stresses (caused by successive heating cycles). Thus, this study aims to investigate the process improvement and the optimization of the printed parts.

In this work, a meshless technique – the Radial Point Interpolation Method (RPIM) – is used to numerically simulate the viscoplastic extrusion process – the initial phase of the FFF. Unlike the FEM, in meshless methods, there is no pre-established relationship between the nodes so the nodal mesh will not face mesh distortions and the discretization can easily be modified by adding or removing nodes from the initial nodal mesh. The accuracy of the obtained results highlights the importance of using meshless techniques in this field.

Meshless methods show particular relevance in this topic since the nodes can be distributed to match the layer-by-layer growing condition of the printing process.

Using the flow formulation combined with the heat transfer formulation presented here for the first time within an in-house RPIM code, an algorithm is proposed, implemented and validated for benchmark examples.

This study implements the fourth-order phase field method (PFM) for modeling fracture in brittle materials. The weak form of the fourth-order PFM requires C¹ basis functions for the crack evolution scalar field in a finite element framework. To address this, non-Sibsonian type shape functions that are nonpolynomial types based on distance measures, are used in the context of natural neighbor shape functions. The capability and efficiency of this method are studied for modeling cracks.

The weak form of the fourth-order PFM is derived from two governing equations for finite element modeling. C⁰ non-Sibsonian shape functions are derived using distance measures on a generalized quad element. Then these shape functions are degree elevated with Bernstein-Bezier (BB) patch to get higher-order continuity (C¹) in the shape function. The quad element is divided into several background triangular elements to apply the Gauss-quadrature rule for numerical integration. Both fourth-order and second-order PFMs are implemented in a finite element framework. The efficiency of the interpolation function is studied in terms of convergence and accuracy for capturing crack topology in the fourth-order PFM.

It is observed that fourth-order PFM has higher accuracy and convergence than second-order PFM using non-Sibsonian type interpolants. The former predicts higher failure loads and failure displacements compared to the second-order model due to the addition of higher-order terms in the energy equation. The fracture pattern is realistic when only the tensile part of the strain energy is taken for fracture evolution. The fracture pattern is also observed in the compressive region when both tensile and compressive energy for crack evolution are taken into account, which is unrealistic. Length scale has a certain specific effect on the failure load of the specimen.

Fourth-order PFM is implemented using C¹ non-Sibsonian type of shape functions. The derivation and implementation are carried out for both the second-order and fourth-order PFM. The length scale effect on both models is shown. The better accuracy and convergence rate of the fourth-order PFM over second-order PFM are studied using the current approach. The critical difference between the isotropic phase field and the hybrid phase field approach is also presented to showcase the importance of strain energy decomposition in PFM.

This study aims to simulate the dendritic growth in Stokes flow by iteratively coupling a domain and boundary type meshless method.

A preconditioned phase-field model for dendritic solidification of a pure supercooled melt is solved by the strong-form space-time adaptive approach based on dynamic quadtree domain decomposition. The domain-type space discretisation relies on monomial augmented polyharmonic splines interpolation. The forward Euler scheme is used for time evolution. The boundary-type meshless method solves the Stokes flow around the dendrite based on the collocation of the moving and fixed flow boundaries with the regularised Stokes flow fundamental solution. Both approaches are iteratively coupled at the moving solid–liquid interface. The solution procedure ensures computationally efficient and accurate calculations. The novel approach is numerically implemented for a 2D case.

The solution procedure reflects the advantages of both meshless methods. Domain one is not sensitive to the dendrite orientation and boundary one reduces the dimensionality of the flow field solution. The procedure results agree well with the reference results obtained by the classical numerical methods. Directions for selecting the appropriate free parameters which yield the highest accuracy and computational efficiency are presented.

A combination of boundary- and domain-type meshless methods is used to simulate dendritic solidification with the influence of fluid flow efficiently.

This paper aims to introduce a new numerical approach based on the local weak form and the Petrov–Galerkin idea to numerically simulation of a predator–prey system with two-species, two chemicals and an additional chemotactic influence.

In the first proceeding, the space derivatives are discretized by using the direct meshless local Petrov–Galerkin method. This generates a nonlinear algebraic system of equations. The mentioned system is solved by using the Broyden’s method which this technique is not related to compute the Jacobian matrix.

This current work tries to bring forward a trustworthy and flexible numerical algorithm to simulate the system of predator–prey on the nonrectangular geometries.

The proposed numerical results confirm that the numerical procedure has acceptable results for the system of partial differential equations.

This paper aims to present a methodology for the detection and categorisation of metal powder particles that are partially attached to additively manufactured lattice structures. It proposes a software algorithm to process micro computed tomography (µCT) image data, thereby providing a systematic and formal basis for the design and certification of powder bed fusion lattice structures, as is required for the certification of medical implants.

This paper details the design and development of a software algorithm for the analysis of µCT image data. The algorithm was designed to allow statistical probability of results based on key independent variables. Three data sets with a single unique parameter were input through the algorithm to allow for characterisation and analysis of like data sets.

This paper demonstrates the application of the proposed algorithm with three data sets, presenting a detailed visual rendering derived from the input image data, with the partially attached particles highlighted. Histograms for various geometric attributes are output, and a continuous trend between the three different data sets is highlighted based on the single unique parameter.

This paper presents a novel methodology for non-destructive algorithmic detection and categorisation of partially attached metal powder particles, of which no formal methods exist. This material is available to download as a part of a provided GitHub repository.

This paper aims to investigate the performance of two novel numerical methods, the face-based smoothed finite element method (FS-FEM) and the edge-based smoothed finite element method (ES-FEM), which employ linear tetrahedral elements, for the purpose of strength assessment of a high-speed train hollow axle.

The calculation of stress for the wheelset, comprising an axle and two wheels, is facilitated through the application of the European axle strength design standard. This standard assists in the implementation of loading and boundary conditions and is exemplified by the typical CRH2 high-speed train wheelset. To evaluate the performance of these two methods, a hollow cylinder cantilever beam is first used as a benchmark to compare the present methods with other existing methods. Then, the strength analysis of a real wheelset model with a hollow axle is performed using different numerical methods.

The results of deflection and stress show that FS-FEM and ES-FEM offer higher accuracy and better convergence than FEM using linear tetrahedral elements. ES-FEM exhibits a superior performance to that of FS-FEM using linear tetrahedral elements, showing accuracy and convergence close to FEM using hexahedral elements.

This study channels the novel methods (FS-FEM and ES-FEM) in the static stress analysis of a railway wheelset. Based on the careful testing of FS-FEM and ES-FEM, both methods hold promise as more efficient tools for the strength analysis of complex railway structures.

In this paper, the standard Peridynamic Timoshenko beam model accounting for the shear deformation is chosen to describe the thick beam kinematics. Unfortunately, when applied to very thin beam structures, the standard Peridynamics (PD) encounters the shear locking phenomenon, leading to incorrect solutions.

PD differs from classical continuum mechanics and other nonlocal theories that do not involve spatial derivatives of the displacement field. PD is based on the integral equation instead of differential equations to handle discontinuities and other singularities.

The shear locking can be successfully alleviated using the developed selective integration method. In particular, this technique has been implemented in the standard PD, which allows an accurate result for a wide range of slenderness from very thin to thick (10 < L/t < 10³) structures. It can also accelerate the computational time for particular dynamic problems using fewer neighboring integration particles. Several numerical examples are solved to demonstrate the effectiveness of the proposed method for modeling beam structures.

The paper highlights the severe shear locking phenomenon in the Peridynamic Timoshenko beam available in the literature, especially for very thin structures. A new alternative for the alleviation of shear locking in the Peridynamic Timoshenko beam, using selective integration. Hence the developed Peridynamic Timoshenko beam model is effective for thin and thick structures. A new peridynamic formulation for the low-velocity impact beam models is presented and validated.

1. The paper highlights the severe shear locking phenomenon in the Peridynamic Timoshenko beam proposed in the literature, especially for very thin structures.

2. The developed Peridynamic Timoshenko beam model based on selective integration is effective for thin and thick structures.

3. A new peridynamic formulation for the low-velocity impact beam models is presented and validated.

The paper highlights the severe shear locking phenomenon in the Peridynamic Timoshenko beam proposed in the literature, especially for very thin structures.

The developed Peridynamic Timoshenko beam model based on selective integration is effective for thin and thick structures.

A new peridynamic formulation for the low-velocity impact beam models is presented and validated.