Search results

The purpose of this paper is to develop an effective numerical approach to assess the nonlinear dynamic responses of a near‐bed submarine pipeline.

A coupled numerical approach is proposed in this paper to assess the nonlinear dynamic responses of this pipeline. The boundary‐element method is first used to get the nonlinear dynamic fluid loading induced by the asymmetric flow. The meshless technique is used to discretize the structure of the pipeline. A numerical example is first presented to verify the effectivity of the present method. Then, the coupled technique is used to simulate the nonlinear dynamic fluid‐structure interaction problem of a near‐bed pipeline. A Newton‐Raphson iteration procedure is used herein to solve the nonlinear system of equations, and the Newmark method is adopted for the time integration.

The presence of seabed results in a large negative lift on a pipeline in a horizontal current. Studies reveal that there exists a critical current velocity, above which the pipeline will become instable, and the critical velocity is significantly affected by the initial gap from the pipeline to the seabed.

The near‐bed submarine pipeline is a widely used structure in marine engineering. This paper originally develops a numerical approach to model this special fluid‐structure interaction problem. It has demonstrated by the examples that the present approach is very effective and has good potential in the practical applications.

This article aims to develop an accurate and efficient meshfree Galerkin method based on the strain smoothing technique for linear elastic continuous and fracture problems.

This paper proposed a generalized linear smoothed meshfree method (LSMM), in which the compatible strain is reconstructed by the linear smoothed strains. Based on the idea of the weighted residual method and employing three linearly independent weight functions, the linear smoothed strains can be created easily in a smoothing domain. Using various types of basic functions, LSMM can solve the linear elastic continuous and fracture problems in a unified way.

On the one hand, the LSMM inherits the properties of high efficiency and stability from the stabilized conforming nodal integration (SCNI). On the other hand, the LSMM is more accurate than the SCNI, because it can produce continuous strains instead of the piece-wise strains obtained by SCNI. Those excellent performances ensure that the LSMM has the capability to precisely track the crack propagation problems. Several numerical examples are investigated to verify the accurate, convergence rate and robustness of the present LSMM.

This study provides an accurate and efficient meshfree method for simulating crack growth.

The purpose of the article is to extend the node-based smoothed point interpolation method (NS-PIM) for soil consolidation analysis based on the Biot’s theory.

The shape functions for displacements and pore pressures are constructed using the PIM separately, leading to the Kronecker delta property and easy implementation of essential boundary conditions. Then, a benchmark problem of 2D consolidation under ramp load is solved to investigate the validity of this application. Meanwhile, convergence features of different solutions are studied. Furthermore, the incompressible and impermeable condition under instant load is investigated. The results calculated by the NS-PIM solution with different orders of shape functions are compared. Finally a 2D consolidation problem in construction period is solved. An error estimation method is applied to check the mesh quality.

The results of the NS-PIM solution show good agreement with those certified results. Useful convergence features are found when comparing the results of the NS-PIM and the FEM solutions. A simple method is introduced to estimate the errors of the model with rough grids. The convergence features and error estimation method can be applied to check the mesh quality and get accurate results. More stable results can be obtained using the NS-PIM solution with lower order of pore pressure shape functions under the incompressible and impermeable condition.

It cannot be denied that the calculation of NS-PIM solution takes more time than that of the FEM solution, and more work needs to be carried out to optimize the NS-PIM solution. Also, in further study, the feasibility of more complicated and practical engineering problems can still be probed in the NS-PIM solution.

This paper introduced a method for the consolidation analysis on the situation of construction loads (“ramp load”) using the NS-PIM which is quite indispensable in many foundation problems. Also, more stable results can be obtained using the NS-PIM solution with lower order of pore pressure shape functions than that with same order of shape functions.

This study first focuses on the situation of construction loads (“ramp load”) in the NS-PIM consolidation analysis which is quite indispensable in many foundation problems. An error estimation method is introduced to evaluate the mesh quality and get accurate values based on the convergence features of the FEM and NS-PIM solutions. Then, the incompressible and impermeable condition under instant load is investigated, and the analysis show that the NS-PIM with lower order of pore pressure shape functions can get stable results in such conditions.

This study aims to propose a new numerical method for solving non-linear partial differential equations on irregular domains.

The main aim of the current paper is to propose a local meshless collocation method to solve the two-dimensional Klein-Kramers equation with a fractional derivative in the Riemann-Liouville sense, in the time term. This equation describes the sub-diffusion in the presence of an external force field in phase space.

First, the authors use two finite difference schemes to discrete temporal variables and then the radial basis function-differential quadrature method has been used to estimate the spatial direction. To discrete the time-variable, the authors use two different strategies with convergence orders O(τ1+γ) and O(τ2−γ) for 0 < γ < 1. Finally, some numerical examples have been presented to show the high accuracy and acceptable results of the proposed technique.

The proposed numerical technique is flexible for different computational domains.

This paper aims to first apply more advanced anisotropic yield criterions as Yld91 and Yld2004 to spherical indentation simulations, and investigate plastic anisotropy identified from indentation simulations following different yield criterions (Hill48, Yld91, Yld2004) to discover laws. It also aims to compare the difference in plastic anisotropy identified from indentation on three yield criterions and evaluate the applicability of plastic anisotropy.

This paper uses indentation simulations on different yield criterions to identify plastic anisotropy. First, the trust-region techniques based on the nonlinear least-squares method are used to determine anisotropy coefficients of Yld91 and Yld2004. Then, Yld91 and Yld2004 are implemented into ABAQUS software using user-defined material (UMAT) subroutines with the proposed universal structure. Finally, through considering comprehensively the key factors, the locations of the optimal data acquisition points in indentation simulations on different yield criterions are determined. And, the identified stress–strain curves are compared with experimental data.

This paper discovers that indentation on Yld2004 is able to fully identify difference in equivalent plastic strain between 0° and 90° directions when indentation depth h_t is relatively smaller. And, this research demonstrates conclusively that plastic anisotropy identified from indentation on Yld2004 and Yld91 is more applicable at larger strains than that on Hill48, and that on Yld2004 is more applicable than that on Yld91, overall. In addition, the method on the determination of the locations of the optimal data acquisition points is demonstrated to be also valid for anisotropic material.

This paper first investigates plastic anisotropic properties and laws identified from indentation simulations following more advanced anisotropic yield criterions and provides reference for later research.

To increase the use of the meshless method, a hybrid stress method is introduced into the meshless method.

The method is based on the radial point interpolation method (RPIM). According to the Hellinger Reissner principle, stress functions are introduced into the solution procedure. Finite elements are used as background cells for integration. All cells are divided into two types – the H cells, which are around the traction-free circular boundary, and the G cells. For the H cells, stress functions in polar coordinates are created. For the G cells, 12-parameter stress functions in Cartesian coordinates are used. Stress functions are based on equilibrium equations and stress compatible equation.

Numerical results show that this method is reliable.

Hybrid stress methods have been applied to finite element methods, but the finite element methods have not been applied into meshless methods.

The purpose of this paper is to study thermal (natural) convection in nine different containers involving the same area (area= 1 sq. unit) and identical heat input at the bottom wall (isothermal/sinusoidal heating). Containers are categorized into three classes based on geometric configurations [Class 1 (square, tilted square and parallelogram), Class 2 (trapezoidal type 1, trapezoidal type 2 and triangle) and Class 3 (convex, concave and triangle with curved hypotenuse)].

The governing equations are solved by using the Galerkin finite element method for various processing fluids (Pr = 0.025 and 155) and Rayleigh numbers (10³ ≤ Ra ≤ 10⁵) involving nine different containers. Finite element-based heat flow visualization via heatlines has been adopted to study heat distribution at various sections. Average Nusselt number at the bottom wall ( Nub¯) and spatially average temperature (θ^) have also been calculated based on finite element basis functions.

Based on enhanced heating criteria (higher Nub¯ and higher θ^), the containers are preferred as follows, Class 1: square and parallelogram, Class 2: trapezoidal type 1 and trapezoidal type 2 and Class 3: convex (higher θ^) and concave (higher Nub¯).

The comparison of heat flow distributions and isotherms in nine containers gives a clear perspective for choosing appropriate containers at various process parameters (Pr and Ra). The results for current work may be useful to obtain enhancement of the thermal processing rate in various process industries.

Heatlines provide a complete understanding of heat flow path and heat distribution within nine containers. Various cold zones and thermal mixing zones have been highlighted and these zones are found to be altered with various shapes of containers. The importance of containers with curved walls for enhanced thermal processing rate is clearly established.

The purpose of this paper is to extend the natural neighbour radial point interpolation method (NNRPIM) to the dynamic analysis (free vibrations and forced vibrations) of two‐dimensional, three‐dimensional and bending plate problems.

The NNRPIM shape‐function construction is briefly presented, as are the dynamic equations and the mode superposition method is used in the forced vibration analysis. Several benchmark examples of two‐dimensional and plate bending problems are solved and compared with the three‐dimensional NNRPIM formulation. The obtained results are compared with the available exact solutions and the finite element method (FEM) solutions.

The developed NNRPIM approach is a good alternative to the FEM for the solution of dynamic problems, once the obtained results with the EFGM shows a high similarity with the obtained FEM results and for the majority of the studied examples the NNRPIM results are more close to the exact solution results.

Comparing the FEM and the NNRPIM, the computational cost of the NNRPIM is higher.

The paper demonstrates extension of the NNRPIM to the dynamic analysis of two‐dimensional, three‐dimensional and bending plate problems. The elimination of the shear‐locking phenomenon in the NNRPIM plate bending formulation. The various solved examples prove a high convergence rate and accuracy of the NNRPIM.

This meshless collocation method is applicable not only to the Helmholtz equation with Dirichlet boundary condition but also mixed boundary conditions. It can calculate not only the high wavenumber problems, but also the variable wave number problems.

In this paper, the authors developed a meshless collocation method by using barycentric Lagrange interpolation basis function based on the Chebyshev nodes to deduce the scheme for solving the three-dimensional Helmholtz equation. First, the spatial variables and their partial derivatives are treated by interpolation basis functions, and the collocation method is established for solving second order differential equations. Then the differential matrix is employed to simplify the differential equations which is on a given test node. Finally, numerical experiments show the accuracy and effectiveness of the proposed method.

The numerical experiments show the advantages of the present method, such as less number of collocation nodes needed, shorter calculation time, higher precision, smaller error and higher efficiency. What is more, the numerical solutions agree well with the exact solutions.

Compared with finite element method, finite difference method and other traditional numerical methods based on grid solution, meshless method can reduce or eliminate the dependence on grid and make the numerical implementation more flexible.

The Helmholtz equation has a wide application background in many fields, such as physics, mechanics, engineering and so on.

This meshless method is first time applied for solving the 3D Helmholtz equation. What is more the present work not only gives the relationship of interpolation nodes but also the test nodes.

The study aims to present a new meshless method based on the Taylor series expansion. The compact supported radial basis functions (CSRBFs) are very attractive, can be considered as a numerical tool for the engineering problems and used to obtain the trial solution and its derivatives without differentiating the basis functions for a meshless method. A meshless based on the CSRBF and Taylor series method has been developed for the solutions of engineering problems.

This paper is devoted to present a truly meshless method which is called a radial basis Taylor series method (RBTSM) based on the CSRBFs and Taylor series expansion (TSE). The basis function and its derivatives are obtained without differentiating CSRBFs.

The RBTSM does not involve differentiation of the approximated function. This property allows us to use a wide range of CSRBF and weight functions including the constant one. By using a different number of terms in the TSE, the global convergence properties of the RBTSM can be improved. The global convergence properties are satisfied by the RBTSM. The computed results based on the RBTSM shows excellent agreement with results given in the open literature. The RBTSM can provide satisfactory results even with the problem domains which have curved boundaries and irregularly distributed nodes.

The CSRBFs have been widely used for the construction of the basic function in the meshless methods. However, the derivative of the basis function is obtained with the differentiation of the CSRBF. In the RBTSM, the derivatives of the basis function are obtained by using the TSE without differentiating the CSRBF.