极坐标哈密顿体系约当型与弹性楔的佯谬解

讨论了极坐标弹性平面哈密顿体系的当型,并通过约当型的求解,直接给出了相关弹性楔体佯谬问题的解,从理论上阐明了经典弹性力学中某些佯谬问题的出现是由于其对应的是哈密顿体系中特殊的约当型解,同时也很自然地为该类问题提供了一个通用,有效的求解方法。     

平面粘性流体扰动与哈密顿体系

通过变分原理,将哈密顿体系的理论引入到平面粘性流体扰动的问题中,导出一套哈密顿算子矩阵的本征函数向量展开求解问题的方法。基于直接法求解流体力学基本方程,导出流场一般特征关系,通过本征值的求解及本征向量的叠加,得到波扰动解,继可分析流场端部效应。从而在该领域用在哈密顿体系下辛几何空间中研究问题的方法代替了传统在拉格朗日体系欧几里德空间分析问题的方法。为流体力学的研究提供一条新途径。     

非线性轨迹优化问题的保辛自适应求解方法

非线性轨迹优化问题一般是一个非线性最优控制问题。将非线性系统的最优控制问题导入到哈密顿体系的辛几何空间当中,基于对偶变量变分原理提出了求解非线性最优控制问题的一种保辛自适应方法。以时间区段两端协态作为独立变量,在时间区段内采用拉格朗日插值近似状态和协态变量,并利用对偶变量变分原理将非线性最优控制问题转化为非线性方程组的求解,保持了哈密顿系统的辛几何结构。并进一步,提出了基于多层次迭代的自适应算法,提高了非线性最优控制问题的求解效率。数值实验验证了该算法在求解非线性轨迹优化问题中的有效性。     

受r^n分布载荷的楔：佯谬的解决

对表面受与ｒｎ（ｎ≥０）成正比的分布载荷的楔，当楔顶角２α与ｎ之间满足一定关系时，经典解为无穷大，这是一个佯谬．本文采用复变函数方法，研究了这个佯谬的所有情形，并发现存在二次佯谬，即对某些特定的（ｎ，α），佯谬解仍为无穷大，对此本文也予以解决     

受一般载荷的楔:佯谬的解决

对于一个楔,当其表面受有与r~n(n≥0)成正比的外载荷时,按照经典弹性力学的解,对于顶角2α为π或2π的楔,其应力为无限大。这个佯谬,对于n=0的情形,已由Dempscy和T.C.T.Ting解决。本文研究n>0的一般情形。     

固支矩形Mindlin板弯曲问题的辛求解体系

本文利用二类变量广义变分原理推出了Mindlin板弯曲问题的Hamilton体系,利用辛几何方法对全状态向量进行分离变量,得到相应的横向本征问题,在求出其本征值后,按本征函数展开法导出了原问题的辛本征通解。给出了一个承受集中载荷的四边固支矩形薄板的算例,按本文求解体系得到的解与经典解吻合较好。本文直接从Mindlin板弯曲问题出发,在其Hamilton体系内使用辛几何方法给出了的一套新的求解体系,突破了传统解法的局限性,具有一般性及较高的理论推广价值。     

顶端受集中力偶或表面受均匀载荷的圆柱型正交各向异性楔:佯谬的解决

对楔顶角为２α的圆柱型正交各向异性楔，其顶端受集中力偶的经典解当α满足ｔｇλα＝λα（λ＝（ａ６６＋２ａ１２＋２ａ１１）／ａ１１，ａ６６、ａ１２、ａ１１为弹性常数）时，应力成为无穷大；其表面受均匀载荷的经典解当α满足ｓｉｎλα＝０（对称变形）或ｔｇλα＝λα（反对称变形）时，也有同样问题。这二个佯谬均由文中利用叠加齐次解的方法解决     

预端受集中力偶或表面受均匀载荷的圆柱型正交各向异性楔：佯谬…

对楔顶角为２α的圆柱型正交各向异性楔，其顶端受集中力偶的经典解当α满足ｔｇλα＝λα（λ＝√（ａ６６＋２ａ１２＋２ａ１１）／ａ１１，ａ６６、ａ１２、ａ１１为弹性常数）时，应力成为无穷大；其表面受均匀载荷的经典解当ａ满足ｓｉｎλα＝０（对称变形）或ｔｇλα＝λα（反对称变形）时，也有同样问题。这二个佯谬均由文中利用叠加齐次解的方法解决。     

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

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

四边任意支承条件下弹性矩形薄板弯曲问题的解析解

利用辛几何法推导出了四边为任意支承条件下矩形薄板弯曲的解析解。在分析过程中首先把矩形薄板弯曲问题表示成Hamilton正则方程，然后利用辛几何方法对全状态相变量进行分离变量，求出其本征值后，再按本征函数展开的方法求出四边为任意支承条件下矩形薄板弯曲的解析解。由于在求解过程中并不需要人为的事先选取挠度函数，而是从弹性矩形薄板弯曲的基本方程出发，直接利用数学的方法求出问题的解析解，使得这类问题的求解更加理论化和合理化。文中的最后还给出了计算实例来验证本文方法的正确性。     

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     

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.     

SYMPLECTIC DUALITY SYSTEM ON PLANE MAGNETOELECTROELASTIC SOLIDS

By means of the generalized variable principle of magnetoelectroelastic solids, the plane magnetoelectroelastic solids problem was derived to Hamiltonian system. In symplectic geometry space, which consists of original variables, displacements, electric potential and magnetic potential, and their duality variables, lengthways stress, electric displacement and magnetic induction, the effective methods of separation of variables and symplectic eigenfunction expansion were applied to solve the problem. Then all the eigen-solutions and the eigen-solutions in Jordan form on eigenvalue zero can be given, and their specific physical significations were shown clearly. At last, the special solutions were presented with uniform loader, constant electric displacement and constant magnetic induction on two sides of the rectangle domain.     

HAMILTONIAN SYSTEM AND THE SAINT VENANT PROBLEM IN ELASTICITY

ＨＡＭＩＬＴＯＮＩＡＮＳＹＳＴＥＭＡＮＤＴＨＥＳＡＩＮＴＶＥＮＡＮＴＰＲＯＢＬＥＭＩＮＥＬＡＳＴＩＣＩＴＹＺｈｏｎｇＷａｎｘｉｅ（钟万勰）；ＸｕＸｉｎｓｈｅｎｇ（徐新生）；ＺｈａｎｇＨｏｎｇｗｕ（张洪武）（ＲｅｃｅｉｖｅｄＪｕｎｅ５，１９９５）Ａｂｓ...     

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

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

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

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

STRESS SINGULARITY ANALYSIS AT CRACK TIP ON BI-MATERIAL INTERFACES BASED ON HAMILTONIAN PRINCIPLE

In this paper, the stress singularity analysis at the crack tip on elastic bi-material inter-faces is considered. The governing equations of plane elasticity in sectorial domain are derived to be inHamiltonian form via variable substitution and variational principle. The methods of separation ofvariables and conjugate symplectic eigen-function expansion are developed to solve the equations insectorial domain. The general formulae for the solution of stress singularities at the crack tip on bi-ma-terial interfaces are put forward, and a new solution technique for fracture problems is presented.     

Hamiltonian principle based stress singularity analysis near crack corners of multi-material junctions

This paper presents a new method for the stress singularity analysis near the crack corners of a multi-material junctions. The stress singularities near the crack corners of multi-dissimilar isotropic elastic material junctions are studied analytically in terms of the methods developed in Hamiltonian system. The governing equations of plane elasticity in a sectorial domain are derived in Hamiltonian form via variable substitution and variational principle respectively. Both of the methods of global state variable separation and symplectic eigenfunction expansion are used to find the analytical solution of the problem. The relationships among the state vectors in different material spaces are obtained by means of coordinate transformation and consistent conditions between the two adjacent domains. The expression of the original problem is thus changed into a new form where the solutions of symplectic generalized eigenvalues and eigenvectors are needed. The closed form of expressions is established for the stress singularity analysis near the corner with arbitrary vertex angles. Numerical results are presented with several chosen angles and multi-material constants. To show the potential of the new method proposed, a semi-analytical finite element is furthermore developed for the numerical analysis of crack problems.     

Hamilton体系下介电弹性体圆形薄膜的动力学建模与辛求解

采用辛算法研究了Hamilton体系下介电弹性体圆形薄膜的动力学响应。首先，将该问题引入Hamilton对偶变量体系，借助Legendre变换，给出系统的广义动量和Hamilton函数，通过对Hamilton函数作用量的变分，得到Hamilton体系下的正则方程。其次，对于得到的正则方程给出了辛Runge-Kutta的计算格式。最后，采用二级四阶辛Runge-Kutta算法对动力学系统进行了数值求解，和四级四阶经典Runge-Kutta算法进行对比，结果表明，二级四阶辛Runge-Kutta算法具有保能量以及长时间数值稳定的优势，同时说明四级四阶经典Runge-Kutta算法对于步长依赖的局限性。