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.     

辛体系下平面热黏弹性圣维南问题的解析解

借助积分变换,将辛体系引入平面热黏弹性问题,建立了基本问题的对偶方程,并将全 部圣维南问题归结为满足共轭辛正交关系的零本征值本征解问题. 同时,利用变量代换和本 征解展开等技术给出了一套求解边界条件问题的具体方法. 算例讨论了几种典型边界条件问 题,描述了热黏弹性材料的蠕变和松弛特征,体现了这种辛方法的有效性.     

Stokes流问题中的辛本征解方法

通过引入哈密顿体系,将二维Stokes流问题归结为哈密顿体系下的本 征值和本征解问题. 利用辛本征解空间的完备性,建立一套封闭的求解问题方法. 研究结果 表明零本征值本征解描述了基本的流动,而非零本征值本征解则显示着端部效应影响特点. 数值算例给出了辛本征值和本征解的一些规律和具体例子. 这些数值例子说明了端部非规则 流动的衰减规律. 为研究其它问题提供了一条路径.     

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.     

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

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

EIGENVALUE PROBLEM OF A LARGE SCALE INDEFINITE GYROSCOPIC DYNAMIC SYSTEM

Gyroscopic dynamic system can be introduced to Hamiltonian system.Based on an adjoint symplectic subspace iteration method of Hamiltonian gyroscopic system, an adjoint symplectic subspace iteration method of indefinite Hamiltonian function gy- roscopic system was proposed to solve the eigenvalue problem of indefinite Hamiltonian function gyroscopic system.The character that the eigenvalues of Hamiltonian gyroscopic system are only pure imaginary or zero was used.The eigenvalues that Hamiltonian function is negative can be separated so that the eigenvalue problem of positive definite Hamiltonian function system was presented,and an adjoint symplectic subspace iteration method of positive definite Hamiltonian function system was used to solve the separated eigenvalue problem.Therefore,the eigenvalue problem of indefinite Hamiltonian function gyroscopic system was solved,and two numerical examples were given to demonstrate that the eigensolutions converge exactly.     

双层石墨烯系统稳态受迫振动问题中的辛方法

针对悬臂石墨烯系统提出一种求解其稳态受迫振动问题的辛解析方法。基于Eringen非局部理论,将石墨烯层板受迫振动问题导入哈密顿体系。采用边界条件分解技术,将问题化为三种边界条件的子问题。通过辛解析方法,得到由辛本征值和辛本征解表示的双层石墨烯系统受迫振动问题的解析解表达式。数值结果表明,辛本征解级数具有很好的收敛性和精度,并与文献结果吻合;在一定的外载激励下可发生同向振动模式和反向振动模式;在一定的参数下,得到一些新的现象和结论。     

基于辛体系的正交各向异性地基梁解析解

基于二维弹性理论, 利用Hellinger-Reissner变分原理, 通过引入对偶变量, 推导 了双参数地基上正交各向异性梁平面应力问题的辛对偶方程组; 采用分离变量法和本征展 开方法, 将原问题归结为求解零本征值本征解和非零本征值本征解, 得到了适用于任意横纵 比的梁的解析解. 由于在求解过程中不需要事先人为地选取试函数, 而是从梁的基本方程出 发, 直接利用数学方法求出问题的解, 使得问题的求解更加合理化. 其中, 地基对梁的力学 行为的影响看作是侧边边界条件, 类似于外载, 可通过零本征解的线性展开来评价, 非零本 征值本征解对应圣维南原理覆盖的部分. 还利用哈密顿变分原理, 给出了两端固支梁的 一种新的改进边界条件. 编程计算了细梁和深梁等算例, 研究了地基上梁的变形沿着厚度方 向的变化特性, 验证了辛方法的有效性.     

辛方法和弹性圆柱壳在内外压和轴向冲击下的动态屈曲

针对有内压或外压的弹性圆柱壳在轴向冲击载荷耦合作用下的动态屈曲问题,构造哈密顿体系,在辛空间中将临界载荷和动态屈曲模态归结为辛本征值和本征解问题,从而形成一种辛方法。该方法直接得到非轴对称的屈曲模态。数值结果给出了圆柱壳问题的临界载荷和屈曲模态以及一些规律。     

矩形梁受分布荷载的辛解答

在弹性力学平面直角坐标辛体系中,采用分离变量法,放弃齐次边界条件,得到了矩形梁侧边受幂函数形式分布荷载问题的辛解答,给出了这类问题在辛体系中的一般解法,分别对矩形梁受法向和切向分布荷载的问题进行了求解,显示了此方法的有效性.辛解法采用对偶的二类变量进行求解,可同时给出位移和应力;由于辛解法能较好地处理各种边界条件,因此不仅能求解静定问题,也能直接求解静不定问题.     

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     

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

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

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

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

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

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

四边固支矩形薄板自由振动的哈密顿解析解

在哈密顿体系中利用辛几何方法求解了四边固支矩形薄板自由振动问题的解析解。首先,从基本方程出发,将问题表示成Hamilton正则方程,然后利用辛几何方法导出本征值问题,从而得到本征函数解,使之满足边界条件;再由方程组有非零解的条件,最终推导出四边固支矩形薄板的自振频率方程,得到频率的解析解。计算了不同长宽比情况下四边固支矩形薄板的频率,结果与已有文献完全一致。该解法有望推广至更多尚未得到解析解的矩形板的振动问题。     

电磁波导的半解析辛分析

根据电磁波导的Hamilton体系，辛几何可用于任意各向异性材料，而且便于处理不同区段的界面条件，横向的电场和磁场构成了对偶向量．基于Hamilton变分原理用半解析法进行横向离散应当保持体系的辛结构．离散后可以运用应用力学的有效算法，求解其辛本征值问题．每段波导可以引入两端Riccati矩阵，用精细积分法求解其方程组．     

Symplectic eigenfunction expansion theorem for elasticity of rectangular planes with two simply-supported opposite sides

The eigenvalue problem of the Hamiltonian operator associated with plane elasticity problems is investigated.The eigenfunctions of the operator are directly solved with mixed boundary conditions for the displacement and stress in a rectangular region.The completeness of the eigenfunctions is then proved,providing the feasibility of using separation of variables to solve the problems.A general solution is obtained with the symplectic eigenfunction expansion theorem.     

New analytic solutions to 2D transient heat conduction problems with/without heat sources in the symplectic space

The two-dimensional (2D) transient heat conduction problems with/without heat sources in a rectangular domain under different combinations of temperature and heat flux boundary conditions are studied by a novel symplectic superposition method (SSM). The solution process is within the Hamiltonian system framework such that the mathematical procedures in the symplectic space can be implemented, which provides an exceptional direct rigorous derivation without any assumptions or predetermination of the solution forms compared with the conventional inverse/semi-inverse methods. The distinctive advantage of the SSM offers an access to new analytic heat conduction solutions. The results obtained by the SSM agree well with those obtained from the finite element method (FEM), which confirms the accuracy of the SSM.     

STRESS SINGULARITY ANALYSIS OF MULTI-MATERIAL WEDGES UNDER ANTIPLANE DEFORMATION

With the help of the coordinate transformation technique, the symplectic dual solving system is developed for multi-material wedges under antiplane deformation. A virtue of present method is that the compatibility conditions at interfaces of a multi-material wedge are expressed directly by the dual variables, therefore the governing equation of eigenvalue can be derived easily even with the increase of the material number. Then, stress singularity on multi-material wedges under antiplane deformation is investigated, and some solutions can be presented to show the validity of the method. Simultaneously, an interesting phenomenon is found and proved strictly that one of the singularities of a special five-material wedge is independent of the crack direction.     

平面粘弹性问题的辛求解方法

将弹性力学辛对偶求解方法与Laplace变换相结合,提出了一个求解粘弹性平面问题的新方法。首先利用Laplace变换,将粘弹性平面问题转化为一个准弹性问题,在辛弹性力学的框架下,利用分离变量和辛本征展开法对其进行求解,然后由逆变换得到原问题的解。为证明方法的有效性,求解分析了矩形域平面粘弹性圣维南问题,得到了令人满意的结果。