B-spline solution of non-linear singular boundary value problems arising in physiology

We use B-spline functions to develop a numerical method for computing approximations to the solution of non-linear singular boundary value problems associated with physiology science. The original differential equation is modified at singular point then the boundary value problem is treated by using B-spline approximation. The numerical method is tested for its efficiency by considering three model problems from physiology.     

Extended multiscale FEM for design of beams and frames with complex topology

A numerical method for design of beams and frames with complex topology is proposed. The method is based on extended multi-scale finite element method where beam finite elements are used on coarse scale and continuum elements on fine scale. A procedure for calculation of multi-scale base functions, up-scaling and downscaling techniques is proposed by using a modified version of window method that is used in computational homogenization. Coarse scale finite element is embedded into a frame of a material that is representing surrounding structure in a sense of mechanical properties. Results show that this method can capture displacements, shear deformations and local stress-strain gradients with significantly reduced computational time and memory comparing to full scale continuum model. Moreover, this method includes a special hybrid finite elements for precise modelling of structural joints. Hence, the proposed method has a potential application in large scale 2D and 3D structural analysis of non-standard beams and frames where spatial interaction between structural elements is important.     

The Cubic B-Spline Method for a Class of Caputo-Fabrizio Fractional Differential Equations北大核心CSCD

基于分数阶微积分基本定理和三次B样条理论,构造了求解线性Caputo-Fabrizio型分数阶微分方程数值解的三次B样条方法,利用分数阶微积分基本定理将初值问题转化为关于解函数的表达式,再使用三次B样条函数逼近表达式中积分项的被积函数,进而计算了一类Caputo-Fabrizio型分数阶微分方程的数值解.给出了所构造的三次B样条方法的误差估计、收敛性和稳定性的理论证明.数值实验表明,该文数值方法在求解一类Caputo-Fabrizio型分数阶微分方程数值解时具有一定的可行性和有效性,且计算精度和计算效率优于现有的两种数值方法.     

Three-dimensional complex variable element-free Galerkin method

The complex variable element-free Galerkin (CVEFG) method is an efficient meshless Galerkin method that uses the complex variable moving least squares (CVMLS) approximation to form shape functions. In the past, applications of the CVMLS approximation and the CVEFG method are confined to 2D problems. This paper is devoted to 3D problems. Computational formulas and theoretical analysis of the CVMLS approximation on 3D domains are developed. The approximation of a 3D function is formed with 2D basis functions. Compared with the moving least squares approximation, the CVMLS approximation involves fewer coefficients and thus consumes less computing times. Formulations and error analysis of the CVEFG method to 3D elliptic problems and 3D wave equations are provided. Numerical examples are given to verify the convergence and accuracy of the method. Numerical results reveal that the CVEFG method has better accuracy and higher computational efficiency than other methods such as the element-free Galerkin method.     

Free vibration and buckling analysis of plates using B-spline wavelet on the interval Mindlin element

A finite element method (FEM) of B-spline wavelet on the interval (BSWI) is used in this paper to solve the free vibration and buckling problems of plates based on Reissner–Mindlin theory. By aid of the high accuracy of B-spline functions approximation for structural analysis, the proposed method could obtain a fast convergence and a satisfying numerical accuracy with fewer degrees of freedoms (DOF). The numerical examples demonstrate that the present BSWI method achieves the high accuracy compared to the exact solution and others existing approaches in the literatures. The BSWI finite element has potential to be used as a numerical method in analysis and design.     

Second-order two-scale method for bending behavior analysis of composite plate with 3-D periodic configuration and its approximation

This paper considers the bending behaviors of composite plate with 3-D periodic configuration.A second-order two-scale(SOTS)computational method is designed by means of construction way.First,by 3-D elastic composite plate model,the cell functions which are defined on the reference cell are constructed.Then the effective homogenization parameters of composites are calculated,and the homogenized plate problem on original domain is defined.Based on the Reissner-Mindlin deformation pattern,the homogenization solution is obtained.And then the SOTS’s approximate solution is obtained by the cell functions and the homogenization solution.Second,the approximation of the SOTS’s solution in energy norm is analyzed and the residual of SOTS’s solution for 3-D original in the pointwise sense is investigated.Finally,the procedure of SOTS’s method is given.A set of numerical results are demonstrated for predicting the effective parameters and the displacement and strains of composite plate.It shows that SOTS’s method can capture the 3-D local behaviors caused by3-D micro-structures well.     

A two-level nesting smoothed meshfree method for structural dynamic analysis

This paper presents a novel two-level nesting smoothed meshfree method (NSMM), which significantly improves the computational efficiency of the meshfree Galerkin methods without losing their accuracy, thus facilitates the employment of meshfree methods in applications where background integration cells would be prohibitively expensive. In the NSMM, the system stiffness matrix is calculated using the general smoothing strain technique over the two-level nesting smoothing sub-domains where fewer integration points are used and the costly derivative computation of meshfree shape functions is avoided. The accuracy, efficiency and stability of the present method are assessed by virtue of several numerical examples for problems involving free and forced vibration analysis of the linear elastic continua and dynamic crack response of elastic solid. The results reveal that the NSMM stands out and achieves better performance compared to other existing approaches in the literature.     

High order approximations using spline-based differential quadrature method: Implementation to the multi-dimensional PDEs

A new differential quadrature method based on cubic B-spline is developed for the numerical solution of differential equations. In order to develop the new approach, the B-spline basis functions are used on the grid and midpoints of a uniform partition. Some error bounds are obtained by help of cubic spline collocation, which show that the method in its classic form is second order convergent. In order to derive higher accuracy, high order perturbations of the problem are generated and applied to construct the numerical algorithm. A new fourth order method is developed for the numerical solution of systems of second order ordinary differential equations. By solving some test problems, the performance of the proposed methods is examined. Also the implementation of the method for multi-dimensional time dependent partial differential equations is presented. The stability of the proposed methods is examined via matrix analysis. To demonstrate the applicability of the algorithms, we solve the 2D and 3D coupled Burgers’ equations and 2D sine-Gordon equation as test problems. Also the coefficient matrix of the methods for multi-dimensional problems is described to analyze the stability.     

The construction of plane elastomechanics and Mindlin plate elements of B-spline wavelet on the interval

Based on two-dimensional tensor product B-spline wavelet on the interval (BSWI), a class of C0 type plate elements is constructed to solve plane elastomechanics and moderately thick plate problems. Instead of traditional polynomial interpolation, the scaling functions of two-dimensional tensor product BSWI are employed to form the shape functions and construct BSWI elements. Unlike the process of direct wavelets adding in the previous work, the elemental displacement field represented by the coefficients of wavelets expansions is transformed into edges and internal modes via the constructed transformation matrix in this paper. The method combines the versatility of the conventional finite element method (FEM) with the accuracy of B-spline functions approximation and various basis functions for structural analysis. Some numerical examples are studied to demonstrate the proposed method and the numerical results presented are in good agreement with the closed-form or traditional FEM solutions.     

Numerical solutions of the modified Burgers’ equation by Petrov-Galerkin method

The modified Burgers’ equation (MBE) is solved numerically by the Petrov-Galerkin method using a linear hat function as the trial function and a cubic B-spline function as the test function. Product approximation has been used in this method. A linear stability analysis of the scheme shows it to be unconditionally stable. The accuracy of the presented method is demonstrated by two test problems. The numerical results are found in good agreement with the exact solutions.     

Trigonometric Tension B-Spline Method for the Solution of Problems in Calculus of Variations

In this paper, the tension B-spline collocation method is used for finding the solution of boundary value problems which arise from the problems of calculus of variations. The problems are reduced to an explicit system of algebraic equations by this approximation. We apply some numerical examples to illustrate the accuracy and implementation of the method.     

A Petrov-Galerkin method for equal width equation

The equal width equation is solved numerically by Petrov-Galerkin method using linear hat function and quadratic B-spline function as trial and test functions respectively. Product approximation has been used in this method. A linear stability analysis of the scheme is shown to be conditionally stable. Test problems including the single soliton and the interaction of solitons are used to validate the suggested method, which is found to be accurate and efficient. The Maxwellian initial condition pulse and the development of an undular bore are also studied.     

Numerical Solution of Singularly Perturbed Boundary Value Problems Based on Optimal Control Strategy

In this paper, we develop a numerical technique for singularly perturbed boundary value problems using B-spline functions and least square method. The approximate solution derived in this article is convergent to the exact solution and can be applied both to linear and nonlinear models. The numerical examples and computational results illustrate and guarantee a higher accuracy for this technique.     

Elastodynamic Contact Problems for Interface Cracks under Harmonic Loading

2-D fracture dynamics' problems for elastic bimaterials with cracks located at the bonding interface under the oblique time harmonic wave are considered in the study. The system of boundary integral equations for displacements and tractions is derived from Somigliana identity taking the contact interaction of the opposite crack faces into account. For the numerical solution the collocation method with piecewise constant approximation is used. The numerical results are obtained for various values of the angle of the wave incidence and the wave frequency taking the friction effects into account. (© 2011 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

A Galerkin boundary element method based on interpolatory Hermite trigonometric wavelets

A Galerkin boundary element method based on interpolatory Hermite trigonometric wavelets is presented for solving 2-D potential problems defined inside or outside of a circular boundary in this paper. In this approach, an equivalent variational form of the corresponding boundary integral equation for the potential problem is used; the trigonometric wavelets are employed as trial and test functions of the variational formulation. The analytical formulae of the matrix entries indicate that most of the matrix entries are naturally zero without any truncation technique and the system matrix is a block diagonal matrix. Each block consists of four circular submatrices. Hence the memory spaces and computational complexity of the system matrix are linear scale. This approach could be easily coupled into domain decomposition method based on variational formulation. Finally, the error estimates of the approximation solutions are given and some test examples are presented.     

Meshfree Approximation for Stochastic Optimal Control Problems

In this work,we study the gradient projection method for solving a class of stochastic control problems by using a mesh free approximation ap-proach to implement spatial dimension approximation.Our main contribu-tion is to extend the existing gradient projection method to moderate high-dimensional space.The moving least square method and the general radial basis function interpolation method are introduced as showcase methods to demonstrate our computational framework,and rigorous numerical analysis is provided to prove the convergence of our meshfree approximation approach.We also present several numerical experiments to validate the theoretical re-sults of our approach and demonstrate the performance meshfree approxima-tion in solving stochastic optimal control problems.     

Meshless method with ridge basis functions

Meshless collocation methods for the numerical solutions of PDEs are increasingly adopted due to their advantages including efficiency and flexibility, and radial basis functions are popularly employed to represent the solutions of PDEs. Motivated by the advantages of ridge basis function representation of a given function, such as the connection to neural network, fast convergence as the number of terms is increased, better approximation effects and various applications in engineering problems, a meshless method is developed based on the collocation method and ridge basis function interpolation. This method is a truly meshless technique without mesh discretization: it neither needs the computation of integrals, nor requires a partition of the region and its boundary. Moreover, the method is applied to elliptic equations to examine its appropriateness, numerical results are compared to that obtained from other (meshless) methods, and influence factors of accuracy for numerical solutions are analyzed.     

B-spline curve fitting with invasive weed optimization

B-spline curves and surfaces are generally used in computer aided design (CAD), data visualization, virtual reality, surface modeling and many other fields. Especially, data fitting with B-splines is a challenging problem in reverse engineering. In addition to this, B-splines are the most preferred approximating curve because they are very flexible and have powerful mathematical properties and, can represent a large variety of shapes efficiently [1]. The selection of the knots in B-spline approximation has an important and considerable effect on the behavior of the final approximation. Recently, in literature, there has been a considerable attention paid to employing algorithms inspired by natural processes or events to solve optimization problems such as genetic algorithms, simulated annealing, ant colony optimization and particle swarm optimization. Invasive weed optimization (IWO) is a novel optimization method inspired from ecological events and is a phenomenon used in agriculture. In this paper, optimal knots are selected for B-spline curve fitting through invasive weed optimization method. Test functions which are selected from the literature are used to measure performance. Results are compared with other approaches used in B-spline curve fitting such as Lasso, particle swarm optimization, the improved clustering algorithm, genetic algorithms and artificial immune system. The experimental results illustrate that results from IWO are generally better than results from other methods.     

A new numerical approach for the advective-diffusive transport equation

A new numerical solution procedure is presented for the one-dimensional, transient advective-diffusive transport equation. The new method applies Herrera's algebraic theory of numerical methods to the spatial derivatives to produce a semi-discrete approximation. The semi-discrete system is then solved by standard time marching algorithms. The algebraic theory, which involves careful choice of test functions in a weak form statement of the problem, leads to a numerical approximation that inherently accommodates different degrees of advection domination. Algorithms are presented that provide either nodal values of the unknown function or nodal values of both the function and its spatial derivative. Numerical solution of several test problems demonstrates the behavior of the method.