Symplectic solution system for reissner plate bending

Based on the Hellinger-Reissner variatonal principle for Reissner plate bendingand introducing dual variables, Hamiltonian dual equations for Reissner plate bending werepresented. Therefore Hamiltonian solution system can also be applied to Reissner platebending problem, and the transformation from Euclidian space to symplectic space and fromLagrangian system to Hamiltonian system was realized. So in the symplectic space whichconsists of the original variables and their dual variables, the problem can be solved viaeffective mathematical physics methods such as the method of separation of variables andeigenfunction-vector expansion. All the eigensolutions and Jordan canonical formeigensolutions for zero eigenvalue of the Hamiltonian operator matrix are solved in detail, and their physical meanings are showed clearly. The adjoint symplectic orthonormal relation of the eigenfunction vectors for zero eigenvalue are formed. It is showed that the alleigensolutions for zero eigenvalue are basic solutions of the Saint-Venant problem and theyform a perfect symplectic subspace for zero eigenvalue. And the eigensolutions for nonzeroeigenvalue are covered by the Saint-Venant theorem. The symplectic solution method is notthe same as the classical semi-inverse method and breaks through the limit of the traditional semi-inverse solution. The symplectic solution method will have vast application.     

含裂纹纳米板振动问题的哈密顿体系方法

基于裂纹处范德华力效应,采用非局部弹性理论构造纳米板模型,并通过导入哈密顿体系建立含裂纹纳米板振动问题的对偶正则控制方程组。在全状态向量表示的哈密顿体系下,将含裂纹纳米板的固有频率和振型问题归结为广义辛本征值和本征解问题。利用哈密顿体系具有的辛共轭正交关系,得到问题解的级数解析表达式。结合边界条件,得到固有频率与辛本征值的代数方程关系式,进而直接给出固有频率的表达式。数值结果表明,非局部尺寸参数和裂纹长度对纳米板振动的各阶固有频率有直接的影响。对比表明,辛方法是准确且可靠的,可为工程应用提供依据。     

陀螺系统与反对称矩阵辛本征解的计算

用状态向量法,引出陀螺线性系统的广义本征问题,证明了本征向量之间的加权共轭辛正交关系,以及用本征向量对任意状态向量的展开定理。运用反对称矩阵胞块组成的LDL~T分解,将本征方程导向辛本征问题的标准型。这套方法适用于陀螺系统K阵不正定的情形。对于辛本征问题用SH变换将矩阵化为半边三对角线胞块阵或三对角线胞块阵,然后再求解其全部本征解。为陀螺系统的模态分析打下了基础。     

Paradox solution on elastic wedge dissimilar materials

According to the Hellinger-Reissner variational principle and introducing proper transformation of variables , the problem on elastic wedge dissimilar materials can be led to Hamiltonian system, so the solution of the problem can be got by employing the separation of variables method and symplectic eigenfunction expansion under symplectic space, which consists of original variables and their dual variables . The eigenvalue - 1 is a special one of all symplectic eigenvalue for Hamiltonian system in polar coordinate . In general, the eigenvalue - 1 is a single eigenvalue, and the classical solution of an elastic wedge dissimilar materials subjected to a unit concentrated couple at the vertex is got directly by solving the eigenfunction vector for eigenvalue - 1. But the eigenvalue - 1 becomes a double eigenvalue when the vertex angles and modulus of the materials satisfy certain definite relationships and the classical solution for the stress distribution becomes infinite at this moment, that is, the para     

空腔内粘性流问题与哈密顿体系方法

以双板驱动空腔粘性流问题为研究对象,根据其特点建立了哈密顿体系下的对偶正则方程,将问题归结为辛体系下的本征值问题.利用辛本征解空间的完备性、正交性和展开理论,形成一套封闭的求解问题方法.算例的数值结果揭示了一些空腔流动的特点.同时这种方法也为研究其他问题提供了一条思路.     

大型辛矩阵本征问题的逆迭代法

基于共轭辛子空间迭代法,求解了大型辛矩阵的主要本征解。随着迭代的进行,可以无限地逼近其精确解。     

大型陀螺特征值问题的广义Arnoldi减缩算法

基于Ａｒｎｏｌｄｉ法，建立陀螺特征值问题的广义Ａｒｎｏｌｄｉ格式，并利用系统矩阵的反对称特性，得到极其简洁的甚至比对称矩阵Ｌａｎｃｚｏｓ法更为简单的递推格式，可称为陀螺Ａｒｎｏｌｄｉ减缩算法。     

Fracture analysis of mode III crack problems for the piezoelectric bimorph

In this paper, a symplectic method based on the Hamiltonian system is proposed to analyze the interfacial fracture in the piezoelectric bimorph under anti-plane deformation. A set of Hamiltonian governing equations is derived from the Hamiltonian function by introducing dual variables of generalized displacements and stresses which can be expanded in series in terms of the symplectic eigensolutions. With the aid of the adjoint symplectic orthogonality, coefficients of the series are determined by the boundary conditions along the crack faces and along the external geometry. The stress\electric displacement intensity factors and energy release rates (G) directly relate to the first few terms of the nonzero eigenvalue solutions. The two ideal crack boundary conditions, namely the electrically impermeable and permeable crack assumptions, are considered. Numerical examples including the complex mixed boundary conditions are considered to show fracture behaviors of the interface crack and discuss the influencing factors.     

圆柱型正交各向异性弹性楔的佯谬解

圆柱型正交各向异性弹性楔体顶端受有集中力偶的经典解,当顶角满足一定关系时,其应力成为无穷大,这是个佯谬．该文在哈密顿体系下将该问题进行重新求解,即利用极坐标各向异性弹性力学哈密顿体系．在原变量和其对偶变量组成的辛几何空间求解特殊本征值的约当型本征解,从而直接给出该佯谬问题的解析解．结果再次表明经典力学中的弹性楔佯谬解对应的是哈密顿体系下辛几何的约当型解．     

Eigenvalue and stability analysis for transverse vibrations of axially moving strings based on Hamiltonian dynamics

The Hamiltonian dynamics is adopted to solve the eigenvalue problem for transverse vibrations of axially moving strings. With the explicit Hamiltonian function the canonical equation of the free vibration is derived. Non-singular modal functions are obtained through a linear, symplectic eigenvalue analysis, and the symplectic-type orthogonality conditions of modes are derived. Stability of the transverse motion is examined by means of analyzing the eigenvalues and their bifurcation, especially for strings transporting with the critical speed. It is pointed out that the motion of the string does not possess divergence instability at the critical speed due to the weak interaction between eigenvalue pairs. The expansion theorem is applied with the non-singular modal functions to solve the displacement response to free and forced vibrations. It is demonstrated that the modal functions can be used as the base functions for solving linear and nonlinear vibration problems. The project supported by the National Natural Science Foundation of China (10472021, 10421002 and 10032030), the NSFC-RFBR Collaboration Project (1031120166/10411120494) and the Scientific Research Foundation for the Returned Overseas Chinese Scholars, State Education Ministry The English text was polished by Keren Wang     

代数黎卡提方程的求解与辛子空间迭代法

综合共轭辛子空间迭代法以及２＾Ｎ消元迭代算法，利用接触变换下黎卡提方程解的变换公式，给出了代数黎卡提方程的有效解法。即使有│μ│＝１的本征根，解法依然有效。数例表明了方法的有效性。     

On the computation of the main eigen-pairs of the continuous-time linear quadratic control problem

The degeneration of the eigenvalue equation of the discrete-time linear quadraticcontrol problem to the continuous-time one when△t→0?is given first.When thecontinuous-time n-dimensional eigenvalue equation,which has all the eigenvalues located inthe left half plane,has been reduced from the original2n-dimensional one,the present paperproposes that several of the eigenvalues nearest to the imaginary axis be obtained by thematrix transformation A_e=e~A.All the eigenvalues of A_e are in the unit circle,with theeigenvectors unchanged and the original eigenvaiues can be obtained by a logarithmoperation.And several of the eigenvalues of A_e nearest to the unit circle can be calculated bythe dual subspace iteration method.     

哈密顿体系与弹性楔体问题

将哈密体系引入到级坐标下的弹性力学楔体问题,利用该体系辛空间的性质,将问题化为本征值和本征向量求解上,得到了完备的解空间,从而改变了弹性力学传统的拉格朗日体系以应力函数为特征的半逆法的讨论去解决该类问题的思路,给出了一条求解该类问题的直接法。     

粘弹性厚壁筒问题的辛本征解方法

将哈密顿体系引进到粘弹性力学厚壁筒问题中,在辛体系下重新描述了基本问题,即建立了正则方程组。借助于积分变换,得到了拉伸、扭转和弯曲等问题的解以及有边界局部效应的解。将原问题归结为辛几何空间中的零本征值本征解和非零本征值本征解问题,从而建立了一种有效的分析问题方法和数值方法。为解决同类问题提供了一条可行的路径。     

弹性圆柱壳在轴向冲击载荷和温度耦合作用下的屈曲

通过引入哈密顿体系,将临界载荷和临界温度及它们所对应的屈曲模态归结为辛体系下的广义本征值和本征解问题。根据辛本征解的正交性和完备性,给出了全部的且独立存在的屈曲模态。数值结果表明,在轴向冲击载荷和温度耦合作用下,弹性圆柱壳的屈曲呈现出复杂的模式,温度直接影响冲击临界载荷的大小。随着温度的增加,冲击临界荷载降低,最后,文中给出各种条件下的屈曲模态。     

Analytical and numerical methods of symplectic system for Stokes flow in two-dimensional rectangular domain

In this paper, a new analytical method of symplectic system, Hamiltonian system, is introduced for solving the problem of the Stokes flow in a two-dimensional rectangular domain. In the system, the fundamental problem is reduced to an eigenvalue and eigensolution problem. The solution and boundary conditions can be expanded by eigensolutions using adjoint relationships of the symplectic ortho-normalization between the eigensolutions. A closed method of the symplectic eigensolution is presented based on completeness of the symplectic eigensolution space. The results show that fundamental flows can be described by zero eigenvalue eigensolutions, and local effects by nonzero eigenvalue eigensolutions. Numerical examples give various flows in a rectangular domain and show effectiveness of the method for solving a variety of problems. Meanwhile, the method can be used in solving other problems.     

Analytical and numerical methods of symplectic system for Stokes flow in two-dimensional rectangular domain

In this paper,a new analytical method of symplectic system.Hamiltonian system,is introduced for solving the problem of the Stokes flow in a two-dimensional rectangular domain.In the system,the fundamental problem is reduced to all eigenvalue and eigensolution problem.The solution and boundary conditions call be expanded by eigensolutions using ad.ioint relationships of the symplectic ortho-normalization between the eigensolutions.A closed method of the symplectic eigensolution is presented based on completeness of the symplectic eigensolution space.The results show that fundamental flows can be described by zero eigenvalue eigensolutions,and local effects by nonzero eigenvalue eigensolutions.Numerical examples give various flows in a rectangular domain and show effectivenees of the method for solving a variety of problems.Meanwhile.the method can be used in solving other problems.     

陀螺效应对转子横向振动的影响分析

举例说明了在动力转子系统中陀螺效应对实际模型的影响。着重分析了转子陀螺效应对进动角速度、振型以及临界角速度的影响。并应用状态空间法求解陀螺系统的本征值问题。数值结果表明，在一些工程问题中，陀螺力对于转子系统振动特性的影响是不能忽略的。     

基于辛本征值方法的受约束介电弹性体中褶皱不稳定性的分析

辛弹性力学已广泛应用于弹性学中各种边值问题的精确解、计算表面波模式以及预测多层超弹性薄膜中的表面褶皱。本文展示了辛分析框架还可应用于受约束介电弹性体中的表面褶皱。机械和电位移向量是两个基本变量来描述介电弹性体中机械变形与电场紧密耦合。褶皱的临界电压可以通过引入基本变量的对偶变量来从辛本征值问题中解决。本文采用扩展的W-W（Wittrick-Williams）算法和精确的积分方法,准确而高效地解决制定的辛本征值问题的本征值。通过将褶皱电压和波数与有无表面能的褶皱基准结果进行比较,验证了辛分析的有效性。辛分析框架简洁且适用于其他不稳定问题,如分层电介质弹性体、磁弹性不稳定性以及层压复合结构的微观和宏观不稳定性。     

ARNOLDI REDUCTION ALGORITHM FOR LARGE SCALE GYROSCOPIC EIGENVALUE PROBLEM

Based on Arnoldi's method, a version of generalized Arnoldi algorithm has been devel-oped for the reduction of gyroscopic eigenvalue problems. By utilizing the skew symmetry of systemmatrix, a very simple recurrence scheme, named gyroscopic Arnoldi reduction algorithm has been ob-tained, which is even simpler than the Lanczos algorithm for symmetric eigenvalue problems. Thecomplex number computation is completely avoided. A restart technique is used to enable the reductionalgorithm to have iterative characteristics. It has been found that the restart technique is not only ef-fective for the convergence of multiple eigenvalues but it also furnishes the reduction algorithm with atechnique to check and compute missed eigenvalues. By combining it with the restart technique, the al-gorithm is made practical for large-scale gyroscopic eigenvalue problems. Numerical examples are giv-en to demonstrate the effectiveness of the method proposed.