Local and parallel finite element method for solving the biharmonic eigenvalue problem of plate vibration

In this paper, we establish a new local and parallel finite element discrete scheme based on the shifted‐inverse power method for solving the biharmonic eigenvalue problem of plate vibration. We prove the local error estimation of finite element solution for the biharmonic equation/eigenvalue problem and prove the error estimation of approximate solution obtained by the local and parallel scheme. When the diameters of three grids satisfy H⁴ = ?(w²) = ?(h), the approximate solutions obtained by our schemes can achieve the asymptotically optimal accuracy. The numerical experiments show that the computational schemes proposed in this paper are effective to solve the biharmonic eigenvalue problem of plate vibration.     

Approximation of sign-indefinite spectral problems

A variational sign-indefinite eigenvalue problem in an infinite-dimensional Hilbert space is approximated by a problem in a finite-dimensional subspace. We analyze the convergence and accuracy of approximate eigenvalues and eigenelements. The general results are illustrated by a sample scheme of the finite-element method with numerical integration for a one-dimensional sign-indefinite second-order differential eigenvalue problem.     

Fehlerabschätzungen für Eigenwertnäherungen nach der Ersatzkernmethode bei Integralgleichungen

Summary Approximate solutions of the linear integral equation eigenvalue problem can be obtained by replacing the original kernel by an approximate kernel. This procedure results in a linear algebraic eigenvalue problem. In this paper we investigate the order of convergence of this method for simple eigenvalues and corresponding eigenfunctions. Our results confirm some numerically observed superconvergence phenomena.

     

A two-grid discretization scheme for the Steklov eigenvalue problem

In the paper, a two-grid discretization scheme is discussed for the Steklov eigenvalue problem. With the scheme, the solution of the Steklov eigenvalue problem on a fine grid is reduced to the solution of the Steklov eigenvalue problem on a much coarser grid and the solution of a linear algebraic system on the fine grid. Using spectral approximation theory, it is shown theoretically that the two-scale scheme is efficient and the approximate solution obtained by the scheme maintains the asymptotically optimal accuracy. Finally, numerical experiments are carried out to confirm the considered theory.     

Asymptotics of eigenvalues of the nonlinear eigenvalue problem arising from the near mixed-mode crack-tip stress-strain field problems

In the present paper, approximate analytical and numerical solutions to nonlinear eigenvalue problems arising in nonlinear fracture mechanics in studying stress-strain fields near a crack tip under mixed-mode loading are presented. Asymptotic solutions are obtained by the perturbation method (the artificial small parameter method). The artificial small parameter is the difference between the eigenvalue corresponding to the nonlinear eigenvalue problem and the eigenvalue related to the linear “undisturbed” problem. It is shown that the perturbation technique is an effective method of solving nonlinear eigenvalue problems in nonlinear fracture mechanics. A comparison of numerical and asymptotic results for different values of the mixity parameter and hardening exponent shows good agreement. Thus, the perturbation theory technique for studying nonlinear eigenvalue problems is offered and applied to eigenvalue problems arising in fracture mechanics analysis in the case of mixed-mode loading.     

Finite Element Approximation of Eigenvalue Problem for a Coupled Vibration Between Acoustic Field and Plate

1.IntroductionInthispaper,westudyanumericalmethodtocalculateeigen-frequenciesofacoupledvibrationbetweenacousticfieldandplate.Atypicalapplicationofthisresearchistoreduceanoiseinsideacarcausedbyanengineorothersourcesofthesound.OurstudywasmotivatedbytheworkofHagiwaraetal.15].Thebackgroundoftheresearchandsomeapplicationscanbeseenin[5].Werestrictourresearchtotheproblemswhereexactsolutionscanbegiveninaspecialcase.Themainfeatureofourresearchisthemathematicallyrigorousapproachtotheproblem.Weformulat…     

粘性阻尼线性振动系统的复模态特征值问题的一种新的矩阵摄动分析解法

本文提出了对粘性阻尼线性振动系统的复模态二次广义特征值问题进行高效近似求解的一种新的矩阵摄动分析方法,即先将阻尼矩阵分解为比例阻尼部分和非比例阻尼部分之和,并求得系统的比例阻尼实模态特征解;然后以此为初始值,将阻尼矩阵的非比例部分作为对其比例部分的小量修改,利用摄动分析方法简捷地得到系统的复模态特征值问题的近似解.这一新方法适用于振系阻尼分布不十分偏离比例阻尼情况的问题,因此对大阻尼(非过阻尼)振动系统也有效.这是它优于以前提出的基于无阻尼实模态特征解的类似摄动分析方法的重要特点.文中建立了复模态特征值和特征向量的二阶摄动解式,并通过算例证实了其有效性.此外还讨论了利用比例阻尼假定估计阻尼系统固有振动的复特征值的可行性.     

Boundary layer analysis for the eigenvalues of hanging rod

The motion of a naturally straight inextensible flexible elastic hanging rod is formulated and then linearized about the straight solution. To solve this equation by separation of variables, an eigenvalue problem is derived. When the stiffness of the rod is small, the eigenvalue equation is a singular perturbation problem. This paper is devoted to solving this eigenvalue problem by boundary layer analysis when the stiffness is suitably small, especially on the analytic approximate solutions of the first several eigenvalues and eigenfunctions. The first three eigenvalues are also compared with the numerical results computed by a finite difference method. The excellent agreement shows the efficiency of the boundary layer analysis.     

Pseudospectral approximation of eigenvalues of derivative operators with non-local boundary conditions

By taking as a “prototype problem” a one-delay linear autonomous system of delay differential equations we present the problem of computing the characteristic roots of a retarded functional differential equation as an eigenvalue problem for a derivative operator with non-local boundary conditions given by the particular system considered. This theory can be enlarged to more general classes of functional equations such as neutral delay equations, age-structured population models and mixed-type functional differential equations.It is thus relevant to have a numerical technique to approximate the eigenvalues of derivative operators under non-local boundary conditions. In this paper we propose to discretize such operators by pseudospectral techniques and turn the original eigenvalue problem into a matrix eigenvalue problem. This approach is shown to be particularly efficient due to the well-known “spectral accuracy” convergence of pseudospectral methods. Numerical examples are given.     

密封容器组合壳自由振动的精确解

给出了一类密封容器组合壳自由振动问题的精确解,基于Love经典薄壳理论,导出了具有任意经线形状的旋转壳体在轴对称振动时的基本方程,组合壳结构中球壳与柱壳的连接条件是通过连接处的变形连续性和内力平衡关系得出的。问题的数学模型被归结为常微分方程组在球壳和 壳两个区间上的特征值问题。振动模态函数是由Legendre和三角函数构造出来,并且得到了精确的频率方程。所有的计算都是在Maple程序下运行的,无论是精确的符号运算还是具有所需有效数学精度的数值计算,都表明该文所编译的Maple程序是简单而有效的。固有频率的数值结果同文献中有限元法和其它数值方法的结果作了比较。作为一个标准,该文给出的精确解对于检验各种近似方法的精密度是有价值的。     

An efficient spectral‐Galerkin method based on a dimension reduction scheme for eigenvalue problems of Schrödinger equations

In this paper, we propose an efficient spectral‐Galerkin method based on a dimension reduction scheme for eigenvalue problems of Schrödinger equations. Firstly, we carry out a truncation from a three‐dimensional unbounded domain to a bounded spherical domain. By using spherical coordinate transformation and spherical harmonic expansion, we transform the original problem into a series of one‐dimensional eigenvalue problem that can be solved effectively. Secondly, we introduce a weighted Sobolev space to treat the singularity in the effective potential. Using the property of orthogonal polynomials in weighted Sobolev space, the error estimate for the approximate eigenvalues and corresponding eigenfunctions are proved. Error estimates show that our numerical method can achieve spectral accuracy for approximate eigenvalues and eigenfunctions. Finally, we give some numerical examples to demonstrate the efficiency of our algorithms and the correctness of the theoretical results.     

An efficient algorithm for the Schrödinger–Poisson eigenvalue problem

We present a new implementation of the two-grid method for computing extremum eigenpairs of self-adjoint partial differential operators with periodic boundary conditions. A novel two-grid centered difference method is proposed for the numerical solutions of the nonlinear Schrödinger–Poisson (SP) eigenvalue problem.We solve the Poisson equation to obtain the nonlinear potential for the nonlinear Schrödinger eigenvalue problem, and use the block Lanczos method to compute the first k   eigenpairs of the Schrödinger eigenvalue problem until they converge on the coarse grid. Then we perform a few conjugate gradient iterations to solve each symmetric positive definite linear system for the approximate eigenvector on the fine grid. The Rayleigh quotient iteration is exploited to improve the accuracy of the eigenpairs on the fine grid. Our numerical results show how the first few eigenpairs of the Schrödinger eigenvalue problem are affected by the dopant in the Schrödinger–Poisson (SP) system. Moreover, the convergence rate of eigenvalue computations on the fine grid is O(h³)$\mathrm{O}(h^3)$.     

HGAH矩阵的逆特征值问题

In this paper, the inverse eigenvalue problem of Hermitian generalized anti-Hamihonian matrices and relevant optimal approximate problem are considered. The necessary and sufficient conditions of the solvability for inverse eigenvalue problem and an expression of the general solution of the problem are derived. The solution of the relevant optimal approximate problem is given.     

INVERSE EIGENVALUE PROBLEM OF HERMITIAN GENERALIZED ANTI-HAMILTONIAN MATRICES

§1　IntroductionWe considerthe following inverse eigenvalue problem offinding an n-by-n matrix A∈S such thatAxi =λixi,i =1,2 ,...,m,where S is a given set of n-by-n matrices,x1 ,...,xm(m≤n) are given n-vectors andλ1 ,...,λmare given constants.Let X=(x1 ,...,xm) ,Λ=(λ1 ,λ2 ,...,λm) ,then the above inverse eigenvalue problemcan be written as followsProblem Given X∈Cn×m,Λ=(λ1 ,...,λm) ,find A∈S such thatAX =XΛ,where S is a given matrix set.We also discuss the so-called opti…     

解Stokes特征值问题的一种两水平稳定化有限元方法

基于局部Gauss积分,研究了解Stokes特征值问题的一种两水平稳定化有限元方法．该方法涉及在网格步长为H的粗网格上解一个Stokes特征值问题,在网格步长为h=O(H2)的细网格上解一个Stokes问题．这样使其能够仍旧保持最优的逼近精度,求得的解和一般的稳定化有限元解具有相同的收敛阶,即直接在网格步长为h的细网格上解一个Stokes特征值问题．因此,该方法能够节省大量的计算时间．数值试验验证了理论结果．     

Nonlocal boundary-value problems for abstract parabolic equations: well-posedness in Bochner spaces

The well-posedness of the nonlocal boundary-value problem for abstract parabolic differential equations in Bochner spaces is established. The first and second order of accuracy difference schemes for the approximate solutions of this problem are considered. The coercive inequalities for the solutions of these difference schemes are established. In applications, the almost coercive stability and coercive stability estimates for the solutions of difference schemes for the approximate solutions of the nonlocal boundary-value problem for parabolic equation are obtained.     

Backward perturbation analysis of certain characteristic subspaces

Summary This paper gives optimal backward perturbation bounds and the accuracy of approximate solutions for subspaces associated with certain eigenvalue problems such as the eigenvalue problemAx=x, the generalized eigenvalue problem Ax=Bx, and the singular value decomposition of a matrixA. This paper also gives residual bounds for certain eigenvalues, generalized eigenvalues and singular values.This subject was supported by the Swedish Natural Science Research Council and the Institute of Information Processing of the University of Umeå.     

不同理论下圆板特征值之间的解析关系

基于经典板理论(CPT)、一阶剪切变形板理论(FPT)以及Reddy三阶剪切变形板理论(RPT)之间,圆板轴对称特征值问题在数学上的相似性,研究了不同理论之间圆板特征值间的解析关系．将特征值问题的求解转化为代数方程的求解,并导出了不同理论之间圆板特征值的显式精确解析关系．从而,只要已知圆板特征值（临界屈曲载荷和固有频率）的经典结果,便很容易从这些解析关系中得到一阶和三阶理论下圆板特征值的相应结果,这便于工程应用,同时也可检验一阶和三阶理论下板特征值的数值结果的有效性、收敛性以及精确性等问题．     

Adaptive computation of smallest eigenvalues of self‐adjoint elliptic partial differential equations

We consider a new adaptive finite element (AFEM) algorithm for self‐adjoint elliptic PDE eigenvalue problems. In contrast to other approaches we incorporate the inexact solutions of the resulting finite‐dimensional algebraic eigenvalue problems into the adaptation process. In this way we can balance the costs of the adaptive refinement of the mesh with the costs for the iterative eigenvalue method. We present error estimates that incorporate the discretization errors, approximation errors in the eigenvalue solver and roundoff errors, and use these for the adaptation process. We show that it is also possible to restrict to very few iterations of a Krylov subspace solver for the eigenvalue problem on coarse meshes. Several examples are presented to show that this new approach achieves much better complexity than the previous AFEM approaches which assume that the algebraic eigenvalue problem is solved to full accuracy. Copyright © 2010 John Wiley & Sons, Ltd.     

Analytical approximations to nonlinear vibration of an electrostatically actuated microbeam

This paper employs the homotopy analysis method (HAM) to derive analytical approximate solutions for the nonlinear problem with high-order nonlinearity. Such a problem corresponds to the large-amplitude vibration of electrostatically actuated microbeams. The HAM is also optimized to accelerate the convergence of approximate solutions. To verify the accuracy of the present approach, illustrative examples are provided and compared with other analytical and exact solutions.