在哈密顿体系下分析非线性动力学问题

首先将ｎ维未知向量ｑ的二阶非线性力系统Ｍｑ＋Ｇｑ＋Ｋｑ＝Ｆ（ｑ,ｑ,ｔ）转化为与其等价的２ｎ维未知向量ｖ的一阶微分方程ｖ＝Ｈｖ＋ｆ（ｖ,ｔ）,其中非线性部分ｊｉ（ｖ,ｔ）＝０（ｉ＝１,．．．ｎ）,ｆｉ（ｖ,ｔ）＝Ｆｉ－ｎ（ｑ,ｑ,ｔ）（ｉ＝ｎ＋１,．．．２ｎ）;然后给出一种求解ｖ的逐步积分公式,从而将精细积分法进一步推广应用到非线性动力学问题。算例表明本方法的计算量较小且结果合理可靠。     

多自由度非线性动力学方程的能量校准算法

针对非线性动力学方程,通过Taylor展开和Duhamel积分,得到一个具有待定参数的逐步积分求解公式;通过数学变换,将原动力学方程转换为一个能确定待定参数的能量校准方程;最后将该参数回代入逐步积分公式,得到数值解.数值算例的结果说明了该方法的有效性,可以消除算法阻尼和抑制数值解发散,同时,在大步长的条件下也得到了非常准确而且稳定的结果,可以对系统长期性态进行仿真.     

弹性混合板的振动和稳定问题

通过将非线性应变位移关系引入Hellinger-Reissner变分原理,推导出了各向异性板基于弹性理论下的振动和屈曲控制方程.用精细积分法研究了四边简支混合矩形板.精细积分法与传统的有限差分法相比,可以给出计算机精度允许的很精确的数值结果.所以,给出的混合板的振动和稳定的结果可以看作是近似解析的.这些结果可以作为衡量各种板理论准确性的一个标准.而且,非对称层合板中出现的各种耦合影响,比如弯扭耦合,拉弯耦合等在同一组控制方程里被考虑了.     

大转动梁的几何非线性分析讨论

本文借助Ｌａｇｒａｎｇｅ（Ｔ．Ｌ．）法、修正的Ｌａｇｒａｎｇｅ（Ｕ．Ｌ．）法及带有动坐标的迭代法求解梁的几何非线性问题，说明了各自的特点，澄清了若干基本概念，指出动坐标方法实质上就是Ｕ．Ｌ．法，它适合于分析具有大转动梁的问题，并可方便地推广到大转动的板壳问题，同时指出对于几何非线性问题，可以不必区分Ｃａｕｃｈｙ应力和Ｋｉｒｃｈｈｏｆｆ应力。     

反应、扩散化学动力学方程的有效解法

对扩散、化学反应或瞬态温度场问题，给出了具有4阶精度、自起步的隐式时间积分算法．算例显示，其精度和稳定性都好于四阶Runge-Kutta法，并且保留了原系数矩阵的稀疏存储方式和稀疏矩阵的运算规则，使紧缩存储技术和减少计算时间有效的结合．以旋转填充床内的竞争串联反应为算例，表明该算法是有效的．     

THE CURVED BEAM ELEMENT AND ITS CONVERGENCE RATE

The quasi-conforming element of the curved beam and shallow curved beam is given inthis paper.Numerical examples illustrate that the quasi-conforming elements of the curvedbeam and shallow curved beam Which is used to aooroximate the curved beam have betteraccuracy than the straight beam element.The curved beam element constructed bydisplacement method can not satisly rigid body motion condition and the verv fine gridshave to be used in order to satisfy rigid body motion condition approximately.In this paper it is proved that the straight beam element and the quasi-conformingelement of the curved beam and shallow curved beam.when element size is reducedinfinitely,have convergence rate with the same order O(l~2)and when regulur elements areused,I is the element length.     

三维接触问题的拟二维序列解法

在每个接触点引入一系列测动坐标系，使三维接触问题转化为一系列可直接用线性互补的Ｌｅｍｋｅ算法求解的拟二维接触问题序列，逼近原三维接触问题的解，给出了加速收敛的算法。算例表明提出的方法是行之有效的。     

一种黏弹塑性统一本构模型

经过对大量有关统一本构模型的文献资料分析，指出了现有统一本构模型存在的问题，并通过对材料实验数据的分析，指出了黏弹性现象在实验中的表现，并据此将黏弹性引入到弹性黏塑性统一本构模型之中，建立了黏弹塑性统一本构模型，通过模型的数值模拟证明：模型计算结果无论在变形趋势上，还是在数值精度上都与实验数据符合得很好，克服了此前统一本构模型存在的问题。黏弹塑性统一本构模型的产生将统一本构模型的产生将统一本构理论的内涵扩大到黏弹性范围，进而构造了一个黏弹塑性理论的新框架。     

EFFECTIVE SOLUTION METHOD OF CHEMICAL REACTION KINETICS WITH DIFFUSE

The time integration method with four-order accuracy, self-starting and implicit for the diffuse chemical reaction kinetics equation or the transient instantaneous temperature filed equation was presented. The examples show that both accuracy and stability are better than Runge-Kutta method with four-order. The coefficients of the equation are stored with sparse matrix pattern, so an algorithm is presented which combines a compact storage scheme with reduced computation cost. The computation of the competitive and consecutive reaction in the rotating packed bed, taken as examples, shows that the method is effective.     

NOVEL REGULARIZED BOUNDARY INTEGRAL EQUATIONS FOR POTENTIAL PLANE PROBLEMS

The universal practices have been centralizing on the research of regulariza-tion to the direct boundary integal equations (DBIEs). The character is elimination of singularities by using the simple solutions. However, up to now the research of regular ization to the first kind integral equations for plane potential problems has never been found in previous literatures. The presentation is mainly devoted to the research on the regularization of the singular boundary integral equations with indirect unknowns. A novel view and idea is presented herein, in which the regularized boundary integral equations with indirect unknowns without including the Cauchy principal value (CPV) and Hadamard-finite-part (HFP) integrals are established for the plane potential problems. With some numerical results, it is shown that the better accuracy and higher efficiency, especially on the boundary, can be achieved by the present system.