A meshless local Galerkin integral equation method for solving a type of Darboux problems based on the radial basis functions

The main goal of this paper is to solve a class of Darboux problems by converting them into the two-dimensional nonlinear Volterra integral equation of the second kind. The scheme approximates the solution of these integral equations using the discrete Galerkin method together with local radial basis functions, which use a small set of data instead of all points in the solution domain. We also employ the Gauss–Legendre integration rule on the influence domains of shape functions to compute the local integrals appearing in the method. Since the scheme is constructed on a set of scattered points and does not require any background meshes, it is meshless. The error bound and the convergence rate of the presented method are provided. Some illustrative examples are included to show the validity and efficiency of the new technique. Furthermore, the results obtained demonstrate that this method uses much less computer memory than the method established using global radial basis functions. doi:10.1017/S1446181121000377

Exact finite beam element for open thin walled doubly symmetric members under torsional and warping moments

Starting with total potential energy variational principle, the governing equilibrium coupled equations for the torsional-warping static analysis of open thin-walled beams under various torsional and warping moments are derived. The formulation captures shear deformation effects due to warping. The exact closed form solutions for torsional rotation and warping deformation functions are then developed for the coupled system of two equations. The exact solutions are subsequently used to develop a family of shape functions which exactly satisfy the homogeneous form of the governing coupled equations. A super-convergent finite beam element is then formulated based on the exact shape functions. Key features of the beam element developed include its ability to (a) eliminate spatial discretization arising in commonly used finite elements, and (e) eliminate the need for time discretization. The results based on the present finite element solution are found to be in excellent agreement with those based on exact solution and ABAQUS finite beam element solution at a small fraction of the computational and modelling cost involved.

An Alternative Approach to Obtaining Stiffness and Mass Matrices for Three-Dimensional Bernoulli-Euler and Timoshenko Beam-Columns

In the static analysis of beam-column systems using matrix methods, polynomials are using as the shape functions. The transverse deflections along the beam axis, including the axial- flexural effects in the beam-column element, are not adequately described by polynomials. As an alternative method, the element stiffness matrix is modeling using stability parameters. The shape functions which are obtaining using the stability parameters are more compatible with the system’s behavior. A mass matrix used in the dynamic analysis is evaluated using the same shape functions as those used for derivations of the stiffness coefficients and is called a consistent mass matrix. In this study, the stiffness and consistent mass matrices for prismatic three-dimensional Bernoulli-Euler and Timoshenko beam-columns are proposed with consideration for the axial-flexural interactions and shear deformations associated with transverse deflections along the beam axis. The second-order effects, critical buckling loads, and eigenvalues are determined. According to the author’s knowledge, this study is the first report of the derivations of consistent mass matrices of Bernoulli-Euler and Timoshenko beam-columns under the effect of axially compressive or tensile force.

A MESHLESS LOCAL GALERKIN INTEGRAL EQUATION METHOD FOR SOLVING A TYPE OF DARBOUX PROBLEMS BASED ON RADIAL BASIS FUNCTIONS

The main goal of this paper is to solve a class of Darboux problems by converting them into the two-dimensional nonlinear Volterra integral equation of the second kind. The scheme approximates the solution of these integral equations using the discrete Galerkin method together with local radial basis functions, which use a small set of data instead of all points in the solution domain. We also employ the Gauss–Legendre integration rule on the influence domains of shape functions to compute the local integrals appearing in the method. Since the scheme is constructed on a set of scattered points and does not require any background meshes, it is meshless. The error bound and the convergence rate of the presented method are provided. Some illustrative examples are included to show the validity and efficiency of the new technique. Furthermore, the results obtained demonstrate that this method uses much less computer memory than the method established using global radial basis functions.

Wormhole geometries in f(Q) gravity and the energy conditions

Following the recent theory of f(Q) gravity, we continue to investigate the possible existence of wormhole geometries, where Q is the non-metricity scalar. Recently, the non-metricity scalar and the corresponding field equations have been studied for some spherically symmetric configurations in Mustafa (Phys Lett B 821:136612, 2021) and Lin and Zhai (Phys Rev D 103:124001, 2021). One can note that field equations are different in these two studies. Following Lin and Zhai (2021), we systematically study the field equations for wormhole solutions and found the violation of null energy conditions in the throat neighborhood. More specifically, considering specific choices for the f(Q) form and for constant redshift with different shape functions, we present a class of solutions for static and spherically symmetric wormholes. Our survey indicates that wormhole solutions could not exist for specific form function $$f(Q)= Q+ \alpha Q^2$$ f ( Q ) = Q + α Q 2 . To summarize, exact wormhole models can be constructed with violation of the null energy condition throughout the spacetime while being $$\rho \ge 0$$ ρ ≥ 0 and vice versa.

A simplified variant of the time finite element methods based on the shape functions of an axial finite bar

In the field of structural dynamics, the structural responses in the time domain are of major concern. There already exist many methods proposed previously including widely used direct time integration methods such as ones in the β-Newmark family, Houbolt’s method, and Runge-Kutta method. The time finite element methods (TFEM) that followed the well-posed variational statement for structural dynamics are found to bring about a superior accuracy even with large time steps (element sizes), when compared with the results from methods mentioned above. Some high-order time finite elements were derived with the procedure analogous to the conventional finite element methods. In the formulation of these time finite elements, the shape functions are like the ones for a (spatial) 2-order finite beam. In this article, a simplified variant for the TFEM is proposed where the shape functions similar to the ones for a (spatial) axial bar are used. The accuracy in the obtained results of some numerical examples is found to be comparable with the accuracy in the previous results.