A combined programming and iteration algorithm for finite element analysis of three-dimensional contact problems

Comparing with two-dimensional contact problems, three-dimensional frictional contact problems are more difficult to deal with, because of the unknown slip direction of the tangential force and enormous computing time. In order to overcome these difficulties, a combined PQP (Parametric Quadratic Programming) and iteration method is derived in this paper. The iteration algorithm, which alleviates the difficulty of unknown slip direction, is used along with the PQP method to cut down computing costs. Numerical example is given to demonstrate the validity of the present algorithm. The project supported by the Machinary and Electronics Ministry of China     

New solution system for circular sector plate bending and its application

Instead of the biharmonic type equation, a set of new governing equations and solving method for circular sector plate bending is presented based on the analogy between plate bending and plane elasticity problems. So the Hamiltonian system can also be applied to plate bending problems by introducing bending moment functions. The new method presents the analytical solutions for the circular sector plate. The results show that the new method is effective. Project supported by National Natural Science Foundation (No. 19732020) and the Doctoral Research Foundation of China.     

电磁波导的半解析辛分析

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

Adaptive selection of parameters for precise computation of matrix exponential based on padé approximation

讨论了基于Pad\'{e}逼近的矩阵指数精细积分方法中加权系数N 和展开项数q的自适应选择问题. 参数(N,q)的选择直接影响到矩阵指数计算的精度和效 率. 采用矩阵函数逼近理论，研究了参数N和q的增加对精度的影响程度，据此，提出了 参数(N,q)优化组合的递推自适应选择方法. 该方法可以根据矩阵本身的性态选择合适的参 数(N,q),而参数选择的计算量与矩阵指数的计算量相比几乎可以忽略，这对于增强矩阵指 数精细积分方法的适应性和提高计算效率是很有益处的. 算例验证了该方法的正确性和有效性.     

摩擦接触弹塑性分析的数学规划法

本文提出摩擦接触弹塑性分析的数学规划法,是弹性接触问题解法的推广,使之能处理比例加载时的塑性效应。利用所导出接触问题总刚度阵的逆,结合参变量线性互补公式,仍可消去由接触单元而引入的惩罚因子。于是,可以获得一个弹塑性问题、弹性接触和弹塑性接触问题的有限元参数二次规划法模型。所选的几个实例说明了方法的实用性。     

多材料交接点裂纹无摩擦奇性应力场

在文献「１」基础上，针对工程中难于求解的多材料交接点裂纹尖应力奇性分析问题，在于哈密顿原理，通过分离变量与共轭辛本征函数地求解，利用材料间的界面连接条件与坐标变换关系，建立了应力奇性与本征函数求解的解析表达式。由于采取裂纹面接触区模型，因而不再发生振荡奇异性。     

关于离散热传导物理模型的探讨

本文从热传导离散物理模型上论证了集中质量热容矩阵模型可以在离散点上满足热量守恒定律，采用它就可以避免许多热传导时间积分中的不合理现象。几个算例表明了该模型具有良好的精度。     

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

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

矩阵黎卡提方程的精细积分法

选择恰当的参数，将２￣Ｎ类算法用于代数与微分黎卡提方程。证明了算得的解是如此精确，几乎是计算机上的精确解。数例验证了该结论。     

Method of separation of variables and Hamiltonian system

In the theory of mechanics and/or mathematical physics problems in a prismatic domain, the method of separation of variables ususally leads to the Sturm–Liouville-type eigenproblems of self-adjoint operators, and then the eigenfunction expansion method can be used in equation solving. However, a number of important application problems cannot lead to self-adjoint operator for the transverse coordinate. From the minimum potential energy variational principle, by selection of the state and its dual variables, the generalized variational principle is deduced. Then, based on the analogy between the theory of structural mechanics and optimal control, the present article leads the problem to the Hamiltonian system. The finite-dimensional theory for the Hamiltonian system is extended to the corresponding theory of the Hamiltonian operator matrix and adjoint symplectic spaces. The adjoint symplectic orthonormality relation is proved for the whole state eigenfunction vectors, and then the expansion of an arbitrary whole state function vector by the eigenfunction vectors is established. Thus the range of classical method of separation of variables is considerably extended. The eigenproblem derived from a plate bending problem in a strip domain is used for illustration. © 1993 John Wiley & Sons, Inc.