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.     

Elasticity solutions for functionally graded rectangular plates with two opposite edges simply supported

England (2006) [13] proposed a novel method to study the bending of isotropic functionally graded plates subject to transverse biharmonic loads. His method is extended here to functionally graded plates with materials characterizing transverse isotropy. Using the complex variable method, the governing equations of three plate displacements appearing in the expansions of displacement field are formulated based on the three-dimensional theory of elasticity for a transverse load satisfying the biharmonic equation. The solution may be expressed in terms of four analytic functions of the complex variable, in which the unknown constants can be determined from the boundary conditions similar to that in the classical plate theory. The elasticity solutions of an FGM rectangular plate with opposite edges simply supported under 12 types of biharmonic polynomial loads are derived as appropriate sums of the general and particular solutions of the governing equations. A comparison of the present results for a uniform load with existing solutions is made and good agreement is observed. The influence of boundary conditions, material inhomogeneity, and thickness to length ratio on the plate deflection and stresses for the load x²yq are studied numerically.     

Simulating thin plate bending problems by a family of two-parameter homogenization functions

This study develops a simple and effective numerical technique, which aims to accurately and quickly address thin plate bending problems. Based on the given boundary conditions, the thin plate homogenization function is constructed and a family of two-parameter homogenization functions are derived. Then, the superposition of homogenization functions method for the thin plate, the clamped plate, and the simply supported plate is obtained, which is meshless without numerical integration and iteration with the merits of easy-to-program and easy-to-implement. Six numerical experiments are employed to verify the effectiveness, accuracy and convergence of the proposed novel strategy. The proposed method is evaluated by the comparisons with the analytical solutions and the referenced solutions. It can be observed that the proposed method is quite accurate for the thin plate, the clamped plate, and the simply supported plate problems.     

On new analytic free vibration solutions of rectangular thin cantilever plates in the symplectic space

In this paper, we obtain accurate analytic free vibration solutions of rectangular thin cantilever plates by using an up-to-date rational superposition method in the symplectic space. To the authors’ knowledge, these solutions were not available in the literature due to the difficulty in handling the complex mathematical model. The Hamiltonian system-based governing equation is first constructed. The eigenvalue problems of two fundamental vibration problems are formed for a cantilever plate. By symplectic expansion, the fundamental solutions are obtained. Superposition of these solutions are equal to that of the cantilever plate, which yields the analytic frequency equation. The mode shapes are then readily obtained. The developed method yields the benchmark analytic solutions with fast convergence and satisfactory accuracy by rigorous derivation, without assuming any trial solutions; thus, it is regarded as rational, and its applicability to more boundary value problems of partial differential equations represented by plates’ vibration, bending and buckling may be expected.     

The resultant approach to computing vector characteristics of multiparameter polynomial matrices

Known types of resultant matrices corresponding to one-parameter matrix polynomials are generalized to the multiparameter case. Based on the resultant approach suggested, methods for solving the following problems for multiparameter polynomial matrices are developed: computing a basis of the matrix range, computing a minimal basis of the right null-space, and constructing the Jordan chains and semilattices of vectors associated with a multiple spectrum point. In solving these problems, the original polynomial matrix is not transformed. Methods for solving other parametric problems of algebra can be developed on the basis of the method for computing a minimal basis of the null-space of a polynomial matrix. Issues concerning the optimality of computing the null-spaces of sparse resultant matrices and numerical precision are not considered. Bibliography: 19 titles. __________ Translated from Zapiski Nauchnykh Seminarov POMI, Vol. 323, 2005, pp. 182–214.     

To solving problems of algebra for two-parameter matrices. VII

The paper discusses the method of rank factorization for solving spectral problems for two-parameter polynomial matrices. New forms of rank factorization, which are computed using unimodular matrices only, are suggested. Applications of these factorizations to solving spectral problems for two-parameter polynomial matrices of both general and special forms are presented. In particular, matrices free of the singular spectrum are considered. Conditions sufficient for a matrix to be free of the singular spectrum and also conditions sufficient for a basis matrix of the null-space to have neither the finite regular nor the finite singular spectrum are provided. Bibliography: 3 titles.     

Application of Daubechies wavelets approximation to plate bending

In this paper a new wavelet-based approach is presented for solving two-dimensional boundary-value mechanical problems on the example of plate bending. The deflection equation of a bending plate is approximated by two-dimensional Daubechies wavelets using a least-squares Galerkin method. Due to the order of the differential equation in mechanics of plate structures is four, a way to perform the calculations of high order connection coefficients (that is, integrals of products of basis functions with their high order derivatives) is suggested. The implementation of two-dimensional Daubechies scaling functions approximation to plate bending is exhibited numerically in some examples. The results show that this method has good precision and reliability. (© 2006 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

To Solving Multiparameter Problems of Algebra. 4. The AB-Algorithm and its Applications

The paper continues the investigation of methods for factorizing q-parameter polynomial matrices and considers their applications to solving multiparameter problems of algebra. An extension of the AB-algorithm, suggested earlier as a method for solving spectral problems for matrix pencils of the form A - λB, to the case of q-parameter (q ≥ 1) polynomial matrices of full rank is proposed. In accordance with the AB-algorithm, a finite sequence of q-parameter polynomial matrices such that every subsequent matrix provides a basis of the null-space of polynomial solutions of its transposed predecessor is constructed. A certain rule for selecting specific basis matrices is described. Applications of the AB-algorithm to computing complete polynomials of a q-parameter polynomial matrix and exhausting them from the regular spectrum of the matrix, to constructing irreducible factorizations of rational matrices satisfying certain assumptions, and to computing “free” bases of the null-spaces of polynomial solutions of an arbitrary q-parameter polynomial matrix are considered. Bibliography: 7 titles. __________ Translated from Zapiski Nauchnykh Seminarov POMI, Vol. 309, 2004, pp. 127–143.     

A new method of fundamental solutions applied to nonhomogeneous elliptic problems

The classical method of fundamental solutions (MFS) has only been used to approximate the solution of homogeneous PDE problems. Coupled with other numerical schemes such as domain integration, dual reciprocity method (with polynomial or radial basis functions interpolation), the MFS can be extended to solve the nonhomogeneous problems. This paper presents an extension of the MFS for the direct approximation of Poisson and nonhomogeneous Helmholtz problems. This can be done by using the fundamental solutions of the associated eigenvalue equations as a basis to approximate the nonhomogeneous term. The particular solution of the PDE can then be evaluated. An advantage of this mesh-free method is that the resolution of both homogeneous and nonhomogeneous equations can be combined in a unified way and it can be used for multiscale problems. Numerical simulations are presented and show the quality of the approximations for several test examples. AMS subject classification 35J25, 65N38, 65R20, 74J20     

Analysis of Mindlin micro plates with a modified couple stress theory and a meshless method

A modified couple stress theory and a meshless method is used to study the bending of simply supported micro isotropic plates according to the first-order shear deformation plate theory, also known as the Mindlin plate theory. The modified couples tress theory involves only one length scale parameter and thus simplifies the theory, since experimentally it is easier to determine the single scale parameter. The equations governing bending of the first-order shear deformation theory are implemented using a meshless method based on collocation with radial basis functions. The numerical method is easy to implement, and it provides accurate results that are in excellent agreement with the analytical solutions.     

Numerical solution to the deflection of thin plates using the two-dimensional Berger equation with a meshless method based on multiple-scale Pascal polynomials

In this study, we employ Pascal polynomial basis in the two-dimensional Berger equation, which is a fourth order partial differential equation with applications to thin elastic plates. The polynomial approximation method based on Pascal polynomial basis can be readily adapted to obtain the numerical solutions of partial differential equations. However, a drawback with the polynomial basis is that the resulting coefficient matrix for the problem considered may be ill-conditioned. Due to this ill-conditioned behavior, we use a multiple-scale Pascal polynomial method for the Berger equation. The ill-conditioned numbers can be mitigated using this approach. Multiple scales are established automatically by selecting the collocation points in the multiple-scale Pascal polynomial method. This method is also a meshless method because there is no requirement to establish complex grids or for numerical integration. We present the solutions of six linear and nonlinear benchmark problems obtained with the proposed method on complexly shaped domains. The results obtained demonstrate the accuracy and effectiveness of the proposed method, as well showing its stability against large noise effects.     

The equilibrium boundary element method in boundary-value problems of elasticity theory

An equilibrium boundary element method is proposed for solving boundary-value problems in the theory of elasticity, thermo-elasticity, the dynamical theory of elasticity, bar torsion calculations, and the bending of a plate. The idea is to use simultaneously the method of constructing bundles of functions which exactly satisfy the equilibrium equations, the boundary variational equations of mechanics, and the methods of discrete finite-element approximation. The variational method of constructing the resolving boundary equations ensures that the linear system is symmetric and easily coupled to the finite-element method. Since volume integrals are eliminated the dimensions of the problem are reduced by one, but, unlike the boundary element method, there is no need to know the fundamental solutions. The solution of some bar torsion and plate bending problems confirms the high numerical efficiency of the method.     

Numerical solution of the Poisson equation over hypercubes using reduced Chebyshev polynomial bases

We describe how to use new reduced size polynomial approximations for the numerical solution of the Poisson equation over hypercubes. Our method is based on a non-standard Galerkin method which allows test functions which do not verify the boundary conditions. Numerical examples are given in dimensions up to 8 on solutions with different smoothness using the same approximation basis for both situations. A special attention is paid on conditioning problems.     

离散算子差分法的单元函数

给出了弹性力学离散算子差分法的离散格式,并给出了该方法的几个板弯曲单元和平面四边形单元,通过对它们的考察,分析了离散算子差分方法中的离散格式对单元位移函数的反映能力。在离散算子差分方法中,无论单元位移函数是否协调,其位移函数均能在离散格式中得到十分好的再现,说明了离散算子差分方法的离散格式是一种性能很优良的离散格式。     

A modified method of approximate particular solutions for solving linear and nonlinear PDEs

The method of approximate particular solutions (MAPS) was first proposed by Chen et al. in Chen, Fan, and Wen, Numer Methods Partial Differential Equations, 28 (2012), 506–522. using multiquadric (MQ) and inverse multiquadric radial basis functions (RBFs). Since then, the closed form particular solutions for many commonly used RBFs and differential operators have been derived. As a result, MAPS was extended to Matérn and Gaussian RBFs. Polyharmonic splines (PS) has rarely been used in MAPS due to its conditional positive definiteness and low accuracy. One advantage of PS is that there is no shape parameter to be taken care of. In this article, MAPS is modified so PS can be used more effectively. In the original MAPS, integrated RBFs, so called particular solutions, are used. An additional integrated polynomial basis is added when PS is used. In the modified MAPS, an additional polynomial basis is directly added to the integrated RBFs without integration. The results from the modified MAPS with PS can be improved by increasing the order of PS to a certain degree or by increasing the number of collocation points. A polynomial of degree 15 or less appeared to be working well in most of our examples. Other RBFs such as MQ can be utilized in the modified MAPS as well. The performance of the proposed method is tested on a number of examples including linear and nonlinear problems in 2D and 3D. We demonstrate that the modified MAPS with PS is, in general, more accurate than other RBFs for solving general elliptic equations.© 2017 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 33: 1839–1858, 2017     

Pseudospectra of Matrix Polynomials that Are Expressed in Alternative Bases

Spectra and pseudospectra of matrix polynomials are of interest in geometric intersection problems, vibration problems, and analysis of dynamical systems. In this note we consider the effect of the choice of polynomial basis on the pseudospectrum and on the conditioning of the spectrum of regular matrix polynomials. In particular, we consider the direct use of the Lagrange basis on distinct interpolation nodes, and give a geometric characterization of “good” nodes. We also give some tools for computation of roots at infinity via a new, natural, reversal. The principal achievement of the paper is to connect pseudospectra to the well-established theory of Lebesgue functions and Lebesgue constants, by separating the influence of the scalar basis from the natural scale of the matrix polynomial, which allows many results from interpolation theory to be applied. This work was partially funded by the Natural Sciences and Engineering Research Council of Canada, and by the MITACS Network of Centres of Excellence.     

The application of trigonal curve to the Mikhailov–Shabat–Sokolov flows

Resorting to the characteristic polynomial of Lax matrix for the Mikhailov–Shabat–Sokolov hierarchy associated with a \({3 \times 3}\) matrix spectral problem, we introduce a trigonal curve, from which we deduce the associated Baker–Akhiezer function, meromorphic functions and Dubrovin-type equations. The straightening out of the Mikhailov–Shabat–Sokolov flows is exactly given through the Abel map. On the basis of these results and the theory of trigonal curve, we obtain the explicit theta function representations of the Baker–Akhiezer function, the meromorphic functions, and in particular, that of solutions for the entire Mikhailov–Shabat–Sokolov hierarchy.     

应用虚拟功的互等定理求解均载悬臂厚矩形板的弯曲

悬臂矩形板的弯曲问题一直是平板经典理论中的著名难题,利用中厚板虚拟功的互等定理,借助付宝连提出的角点静力边界条件,得到了均布载荷作用下悬臂厚矩形板弯曲的封闭解析解,并采用现代数值方法和计算软件对所得解析解进行了数值计算.结果表明功的互等法是求解中厚板弯曲问题的一个简明有效的方法.     

Reissner板中的夹杂和缺陷问题研究

基于保角变换技术和Faber级数展开,研究了含任意形状夹杂或缺陷的无限大Reissner板弯曲问题．将变换域中单位圆内、外解析函数分别展开成Faber级数,并将波动函数展开成第一类和第二类修正的n阶Bessel函数;利用边界位移、剪力和弯矩连续性条件得到问题的高阶线性方程组．以含椭圆形夹杂和缺陷的无限大Reissner板柱面弯曲为例,进一步给出了数值算例和理论分析．结果表明,对于软夹杂,板内力矩随夹杂与板厚尺寸比a/h变化非常敏感;在含硬夹杂条件下,板内力矩随夹杂尺寸变化相对不敏感．     

Global optimization of polynomial-expressed nonlinear optimal control problems with semidefinite programming relaxation

In this paper, we propose a new deterministic global optimization method for solving nonlinear optimal control problems in which the constraint conditions of differential equations and the performance index are expressed as polynomials of the state and control functions. The nonlinear optimal control problem is transformed into a relaxed optimal control problem with linear constraint conditions of differential equations, a linear performance index, and a matrix inequality condition with semidefinite programming relaxation. In the process of introducing the relaxed optimal control problem, we discuss the duality theory of optimal control problems, polynomial expression of the approximated value function, and sum-of-squares representation of a non-negative polynomial. By solving the relaxed optimal control problem, we can obtain the approximated global optimal solutions of the control and state functions based on the degree of relaxation. Finally, the proposed global optimization method is explained, and its efficacy is proved using an example of its application.