Search results

This paper aims to review the existing bond-based peridynamic (PD) and state-based PD heat conduction models, and further propose a refined bond-based PD thermal conduction model by using the PD differential operator.

The general refined bond-based PD is established by replacing the local spatial derivatives in the classical heat conduction equations with their corresponding nonlocal integral forms obtained by the PD differential operator. This modeling approach is representative of the state-based PD models, whereas the resulting governing equations appear as the bond-based PD models.

The refined model can be reduced to the existing bond-based PD heat conduction models by specifying particular influence functions. Also, the refined model does not require any calibration procedure unlike the bond-based PD. A systematic explicit dynamic solver is introduced to validate 1 D, 2 D and 3 D heat conduction in domains with and without a crack subjected to a combination of Dirichlet, Neumann and convection boundary conditions. All of the PD predictions are in excellent agreement with the classical solutions and demonstrate the nonlocal feature and advantage of PD in dealing with heat conduction in discontinuous domains.

The existing PD heat conduction models are reviewed. A refined bond-based PD thermal conduction model by using PD differential operator is proposed and 3 D thermal conduction in intact or cracked structures is simulated.

The purpose of this paper is to apply the Peridynamic differential operator (PDDO) to incompressible inviscid fluid flow with moving boundaries. Based on the potential flow theory, a Lagrangian formulation is used to cope with non-linear free-surface waves of sloshing water in 2D and 3D rectangular and square tanks.

In fact, PDDO recasts the local differentiation operator through a nonlocal integration scheme. This makes the method capable of determining the derivatives of a field variable, more precisely than direct differentiation, when jump discontinuities or gradient singularities come into the picture. The issue of gradient singularity can be found in tanks containing vertical/horizontal baffles.

The application of PDDO helps to obtain the velocity field with a high accuracy at each time step that leads to a suitable geometry updating for the procedure. Domain/boundary nodes are updated by using a second-order finite difference time algorithm. The method is applied to the solution of different examples including tanks with baffles. The accuracy of the method is scrutinized by comparing the numerical results with analytical, numerical and experimental results available in the literature.

Based on the investigations, PDDO can be considered a reliable and suitable approach to cope with sloshing problems in tanks. The paper paves the way to apply the method for a wider range of problems such as compressible fluid flow.

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.

The paper aims to use a switching technique which allows to couple a nonlocal bond-based Peridynamic approach to the Meshless Local Exponential Basis Functions (MLEBF) method, based on classical continuum mechanics, to solve planar problems.

The coupling has been achieved in a completely meshless scheme. The domain is divided in three zones: one in which only Peridynamics is applied, one in which only the meshless method is applied and a transition zone where a gradual transition between the two approaches takes place.

The new coupling technique generates overall grids that are not affected by ghost forces. Moreover, the use of the meshless approach can be limited to a narrow boundary region of the domain, and in this way, it can be used to remove the “surface effect” from the Peridynamic solution applied to all internal points.

The current study paves the road for future studies on dynamic and static crack propagation problems.

The present work involves analytical and experimental investigation of sloshing in a two-dimensional rectangular tank including the effect of porous baffles to control and/or reduce the wave motion in the sloshing tank. The purpose of this study is to assess the analytical solutions of the drag coefficient effect on porous baffles performance to track free surface motion variation in the sloshing tank by comparison with experimental shake table tests under a range of sway excitation.

The linear second-order ordinary differential equations for liquid sloshing in the rectangular tank were solved using Newmark’s beta method and obtained the analytical solutions for liquid sloshing with dual vertical porous baffles of full submergence depths in a sway-oscillated rectangular tank following the methodology similar to Warnitchai and Pinkaew (1998) and Tait (2008).

The porous baffles significantly reduce wave elevation in the varying filled levels of the tank compared to the baffle-free tank under the range of excitation frequencies. It is observed that the Reynolds number-dependent drag coefficient for porous baffles in the tank can significantly reduce the sloshing elevations and is found to be effective to achieve higher damping compared to the porosity-dependent drag coefficient for porous baffles in the sloshing tank. The analytical model’s response to free surface elevation variations in the sloshing tank was compared with the experiment’s test results. The analytical results matched with shake table test results with a quantitative difference near the first resonant frequency.

The scope of the study is limited to porous baffles performance under range sway motion and three different filling levels in the tank. The porous baffle performance includes Reynolds number dependent drag coefficient to explore the damping effect in the sloshing tank.

The porous baffles with low-level porosities in the sloshing tank have many engineering applications where the first resonant mode of sloshing in the tank is more important. The porous baffle drag coefficient is an important parameter to study the baffle’s damping effect in sloshing tanks. Hence, obtained analytical solution for liquid sloshing in the rectangular tank with Reynolds number as well as porosity-dependent drag coefficient (model 1) and porosity-dependent drag coefficient porous baffles (model 2) performance is discussed. The model’s test results were validated using a series of shake table sloshing experiments for three fill levels in the tank with sway motion at various excitation frequencies covering the first four sloshing resonant modes.

The purpose of this paper is to extend the meshless local exponential basis functions (MLEBF) method to the case of nonlinear and linear, variable coefficient partial differential equations (PDEs).

The original version of MLEBF method is limited to linear, constant coefficient PDEs. The reason is that exponential bases which satisfy the homogeneous operator can only be determined for this class of problems. To extend this method to the general case of linear PDEs, the variable coefficients along with all involved derivatives are first expanded. This expanded form is evaluated at the center of each cloud, and is assumed to be constant over the entire cloud. The solution procedure is followed as in the former version. Nonlinear problems are first converted to a succession of linear, variable coefficient PDEs using the Newton-Kantorovich scheme and are subsequently solved using the aforementioned approach until convergence is achieved.

The results obtained show good performance of the method as solution to a wide range of problems. The results are compared with the well-known methods in the literature such as the finite element method, high-order finite difference method or variants of the boundary element method.

The MLEBF method is a simple yet effective tool for analyzing various kinds of problems. It is easy to implement with high parallelization potential. The proposed method addresses the biggest limitation of the method, and extends it to linear, variable coefficient PDEs as well as nonlinear ones.

The purpose of this paper is to describe a novel implementation of a multispatial method, multitime-scheme subdomain differential algebraic equation (DAE) framework allowing a mix of different space discretization methods and different time schemes by a robust generalized single step single solve (GS4) family of linear multistep (LMS) algorithms on a single body analysis for the first-order nonlinear transient systems.

This proposed method allows the coupling of different numerical methods, such as the finite element method and particle methods, and different implicit and/or explicit algorithms in each subdomain into a single analysis with the GS4 framework. The DAE, which constrains both space and time in multi-subdomain analysis, combined with the GS4 framework ensures the second-order time accuracy in all primary variables and Lagrange multiplier. With the appropriate GS4 parameters, the algorithmic temperature rate variable shift can be matched for all time steps using the DAE. The proposed method is used to solve various combinations of spatial methods and time schemes between subdomains in a single analysis of nonlinear first-order system problems.

The proposed method is capable of coupling different spatial methods for multiple subdomains and different implicit/explicit time integration schemes in the GS4 framework while sustaining second-order time accuracy.

Traditional approaches do not permit such robust and flexible coupling features. The proposed framework encompasses most of the LMS methods that are second-order time accurate and unconditionally stable.

The purpose of this paper is to study large deviations for the solution processes of a stochastic equation incorporated with the effects of nonlocal condition.

A weak convergence approach is adopted to establish the Laplace principle, which is same as the large deviation principle in a Polish space. The sufficient condition for any family of solutions to satisfy the Laplace principle formulated by Budhiraja and Dupuis is used in this work.

Freidlin–Wentzell type large deviation principle holds good for the solution processes of the stochastic functional integral equation with nonlocal condition.

The asymptotic exponential decay rate of the solution processes of the considered equation towards its deterministic counterpart can be estimated using the established 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].

The purpose of this paper is to propose an improved proportional topology optimization (IPTO) algorithm for tackling the stress-constrained minimum volume optimization problem, which can meet the requirements that are to get rid of the problems of numerical derivation and sensitivity calculation involved in the process of obtaining sensitivity information and overcome the drawbacks of the original proportional topology optimization (PTO) algorithm.

The IPTO algorithm is designed by using the new target material volume update scheme and the new density variable update scheme and by introducing the improved density filter (considering the weighting function based on the Gaussian distribution) and Heaviside-type projection operator on the basis of the PTO algorithm. The effectiveness of the IPTO algorithm is demonstrated by solving the stress-constrained minimum volume optimization problems for two numerical examples and being compared with the PTO algorithm.

The results of this paper show that the uses of the proposed strategies contribute to improving the optimized results and the performance (such as the ability to obtain accurate solutions, robustness and convergence speed) of the IPTO algorithm. Compared with the PTO algorithm, the IPTO algorithm has the advantages of fast convergence speed, enhancing the ability to obtain accurate solutions and improving the optimized results.

This paper achieved the author’s intended purpose and provided a new idea for solving the stress-constrained optimization problem under the premise of avoiding obtaining sensitivity information.