The mode III stress/electric intensity factors and singularities analysis for edge-cracked circular piezoelectric shafts

This paper introduces a new analytical method to determine the stress and electric intensity factors for edge-cracked circular piezoelectric shafts using Hamiltonian formalism. The singularities near the crack tip are represented in terms of exponential series that can show the boundary layer effects effectively. A symplectic system is established directly by introducing dual vectors in terms of the symplectic eigenfunctions. The coefficients of the series are determined from the lateral boundary conditions along the crack faces and the outer boundary conditions along the exterior geometric domain. The intensity factors are determined by the first two coefficients of non-zero eigenvalue solutions. Numerical examples for various boundary conditions are given. The influencing parameters on the intensity factors are investigated.     

短波近似的保辛算法

WKBJ短波近似是最常用的有效求解方法之一。保守体系的微分方程可用Hamilton体系的方法描述,其特点是保辛。保辛给出保守体系结构最重要的特性。但WKBJ短波近似却未曾考虑保辛的问题。WKBJ近似可用自变量坐标变换,然后再给出其保辛摄动。数值例题展示了本文变换保辛算法的有效性。     

Closed form stress distribution in 2D elasticity for all boundary conditions

This paper applies a Hamiltonian method to study analytically the stress dis- tributions of orthotropic two-dimensional elasticity in(x,z)plane for arbitrary boundary conditions without beam assumptions.It is a method of separable variables for partial differential equations using displacements and their conjugate stresses as unknowns.Since coordinates(x,z)can not be easily separated,an alternative symplectic expansion is used. Similar to the Hamiltonian formulation in classical dynamics,we treat the x coordinate as time variable so that z becomes the only independent coordinate in the Hamiltonian ma- trix differential operator.The exponential of the Hamiltonian matrix is symplectic.There are homogenous solutions with constants to be determined by the boundary conditions and particular integrals satisfying the loading conditions.The homogenous solutions consist of the eigen-solutions of the derogatory zero eigenvalues(zero eigen-solutions) and that of the well-behaved nonzero eigenvalues(nonzero eigen-solutions).The Jordan chains at zero eigenvalues give the classical Saint-Venant solutions associated with aver- aged global behaviors such as rigid-body translation,rigid-body rotation or bending.On the other hand,the nonzero eigen-solutions describe the exponentially decaying localized solutions usually ignored by Saint-Venant's principle.Completed numerical examples are newly given to compare with established results.     

哈密顿系统正则变换在时变最优控制中的应用

利用哈密顿系统正则变换和生成函数理论求解线性时变最优控制问题,构造了新的最优控制律形式并提出了控制增益计算的保结构算法. 利用生成函数求解最优控制导出的哈密顿系统两端边值问题,并构造线性时变系统的最优控制律,由第2类生成函数所构造的最优控制律避免了末端时刻出现无穷大反馈增益. 控制系统设计中需求解生成函数满足的时变矩阵微分方程组. 根据生成函数与哈密顿系统状态转移矩阵之间的关系,从正则变换的辛矩阵描述出发,导出了求解这组微分方程组的保结构递推算法.为了保持递推计算中的辛矩阵结构,哈密顿系统状态转移矩阵的计算中利用了Magnus级数.      

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

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

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.     

多层层合板圣维南问题的解析解

将哈密尔顿体系理论引入到多层层合板问题之中，建立了一套求解该问题的横向哈密尔顿算子矩阵的本征函数向量展开解法，并成功地求解出圣维南问题的解析解．进一步显示了弹性力学新求解体系的有效性及其应用潜力     

多自由度参激系统稳定性的数值分析

提出多自由度周期参激系统稳定性的数值直接法。通过将扰动方程表示成状态方程形式,再根据Flo-quet理论将扰动解表示成指数特征分量与周期分量之积,并将其周期分量与系统周期系数展成Fourier级数,导出一系列代数方程,建立矩阵特征值问题,从而由数值求解特征值可直接确定参激系统的稳定性。该方法可用于一般周期参激阻尼系统,特征值矩阵不含逆子阵。应用于斜拉索在支座周期运动激励下的参激振动不稳定性分析,数值结果表明该方法的有效性。     

线性哈密顿系统的本征摄动法

利用哈密顿算子辛自共轭的特点讨论了保守哈密顿体系的摄动问题,给出了哈密顿矩阵的本征值与本征向量的二阶摄动分析方法。即当系统在哈密顿框架下进行较小修改时,不重复求解大型哈密顿矩阵的本征问题,只需在原系统的模态参数基础上进行模态分析即可,这种矩阵摄动法给出了修改后矩阵的二阶本征值和本征向量,为一般线性保守体系的本征摄动求解提出了一个新方法。     

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

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

Determination of stress intensity factors by the finite element discretized symplectic method

When rewriting the governing equations in Hamiltonian form, analytical solutions in the form of symplectic series can be obtained by the method of separation of variable satisfying the crack face conditions. In theory, there exists sufficient number of coefficients of the symplectic series to satisfy any outer boundary conditions. In practice, the matrix relating the coefficients to the outer boundary conditions is ill-conditioned unless the boundary is very simple, e.g., circular. In this paper, a new two-level finite element method using the symplectic series as global functions while using the conventional finite element shape functions as local functions is developed. With the available classical finite elements and symplectic series, the main unknowns are no longer the nodal displacements but are the coefficients of the symplectic series. Since the first few coefficients are the stress intensity factors, post-processing is not required. A number of numerical examples as well as convergence studies are given.     

Numerical solutions of linear quadratic control for time-varying systems via symplectic conservative perturbation

Optimal control system of state space is a conservative system, whose approximate method should be symplectic conservation. Based on the precise integration method, an algorithm of symplectic conservative perturbation is presented. It gives a uniform way to solve the linear quadratic control (LQ control) problems for linear time-varying systems accurately and efficiently, whose key points are solutions of differential Riccati equation (DRE) with variable coefficients and the state feedback equation. The method is symplectic conservative and has a good numerical stability and high precision. Numerical examples demonstrate the effectiveness of the proposed method.     

Numerical solutions of linear quadratic control for time-varying systems via symplectic par conservative perturbation

Optimal control system of state space is a conservative system, whose approximate method should be symplectic conservation. Based on the precise integration method, an algorithm of symplectic conservative perturbation is presented.It gives a uniform way to solve the linear quadratic control (LQ control) problems for linear time-varying systems accurately and efficiently, whose key points are solutions of differential Riccati equation (DRE) with variable coefficients and the state feedback equation.The method is symplectic conservative and has a good numerical stability and high precision. Numerical examples demonstrate the effectiveness of the proposed method.     

An analytical symplectic approach to the vibration analysis of orthotropic graphene sheets

A nonlocal continuum orthotropic plate model is proposed to study the vibration behavior of single-layer graphene sheets (SLGSs) using an analytical symplectic approach.A Hamiltonian system is established by introduc-ing a total unknown vector consisting of the displacement amplitude,rotation angle,shear force,and bending moment. The high-order governing differential equation of the vibra-tion of SLGSs is transformed into a set of ordinary differential equations in symplectic space.Exact solutions for free vibra-tion are obtianed by the method of separation of variables without any trial shape functions and can be expanded in series of symplectic eigenfunctions. Analytical frequency equations are derived for all six possible boundary con-ditions. Vibration modes are expressed in terms of the symplectic eigenfunctions.In the numerical examples,com-parison is presented to verify the accuracy of the proposed method. Comprehensive numerical examples for graphene sheets with Levy-type boundary conditions are given.A para-metric study of the natural frequency is also included.     

Mixed-mode thermal stress intensity factors from the finite element discretized symplectic method

A finite element discretized symplectic method is introduced to find the thermal stress intensity factors (TSIFs) under steady-state thermal loading by symplectic expansion. The cracked body is modeled by the conventional finite elements and divided into two regions: near and far fields. In the near field, Hamiltonian systems are established for the heat conduction and thermoelasticity problems respectively. Closed form temperature and displacement functions are expressed by symplectic eigen-solutions in polar coordinates. Combined with the analytic symplectic series and the classical finite elements for arbitrary boundary conditions, the main unknowns are no longer the nodal temperature and displacements but are the coefficients of the symplectic series after matrix transformation. The TSIFs, temperatures, displacements and stresses at the singular region are obtained simultaneously without any post-processing. A number of numerical examples as well as convergence studies are given and are found to be in good agreement with the existing solutions.     

HAMILTONIAN SYSTEM AND THE SAINT VENANT PROBLEM IN ELASTICITY

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

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.     

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

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

薄板理论的正交关系及其变分原理

利用平面弹性与板弯曲的相似性理论，将弹性力学新正交关系中构造对偶向量的思路推广到 各向同性薄板弹性弯曲问题，由混合变量求解法直接得到对偶微分方程并推导了对应的变分 原理. 所导出的对偶微分矩阵具有主对角子矩阵为零矩阵的特点. 发现了两个独立的、对称 的正交关系，利用薄板弹性弯曲理论的积分形式证明了这种正交关系的成立. 在恰当选择对 偶向量后，弹性力学的新正交关系可以推广到各向同性薄板弹性弯曲理论.     

随机和区间非齐次线性哈密顿系统的比较研究及其应用

随着计算机技术的飞速发展,更高效、更稳定和长时间模拟能力更强的数值算法需求迫切.哈密顿系统辛算法与传统算法相比在稳定性和长期模拟方面具有显著优越性.但动力系统中不可避免地存在大量不同程度的不确定性,动力学分析中需要考虑这些不确定性的影响以确保合理有效性. 然而,目前考虑参数不确定性的哈密顿系统响应分析的研究基础还比较薄弱. 为此,本文考虑随机和区间参数不确定性,对两种不确定性非齐次线性哈密顿系统分析计算结果进行了比较研究,从而突破了传统哈密顿系统的局限性, 并应用于结构动力响应评估中. 首先,针对确定性非齐次线性哈密顿系统, 提出了考虑确定性扰动的参数摄动法;在此基础上, 分别提出了随机、区间非齐次线性哈密顿系统的参数摄动法,得到了它们响应界限的数学表达; 随后,用数学理论推导得到了区间响应范围包含随机响应范围的相容性结论; 最后,两个数值算例在较小时间步长下验证了所提方法在结构动力响应中的可行性和有效性,体现了随机、区间哈密顿系统响应结果之间的包络关系,并在较大时间步长下与传统方法相比较凸显了哈密顿系统辛算法的数值计算优势、与蒙特卡洛模拟方法相比较验证了所提方法的精度.