弹性力学的混合方程和Hamilton正则方程

本文指出,在弹性力学基本方程中,按变量分类的位移方程,应力方程以外的第三种混合方程,以及按运算子分类的微分方程、变分原理以外的第三种Hamilton方程,它们正好是对应的。本文讨论了静力的和动力的情况以及它们可能的应用。     

On the numerical solution of tracked vehicle dynamic equations

In this investigation, the solution of the nonlinear dynamic equations of the multibody tracked vehicle systems are obtained using different procedures. In the first technique, which is based on the augmented formulation that employes the absolute Cartesian coordinates and Lagrange multipliers, the generalized coordinate partitioning of the constraint Jacobian matrix is used to determine the independent coordinates and the associated independent differential equations. An iterative Newton-Raphson algorithm is used to solve the nonlinear constraint equations for the dependent variables. The numerical problems encountered when one set of independent coordinates is used during the simulation of large scale tracked vehicle systems are demonstrated and their relationship to the track dynamics is discussed. The second approach employed in this investigation is the velocity transformation technique. One of the versions of this technique is discussed in this paper and the numerical problems that arise from the use of inconsistent system of kinematic equations are reported. In the velocity transformation technique, the tracked vehicle system is assumed to consist of two kinematically decoupled subsystems; the first subsystem consists of the chassis, the rollers, the sprocket and the idler, while the second subsystem consists of the track which is represented as a closed kinematic chain that consists of rigid links connected by revolute joints. It is demonstrated that the use of one set of recursive equations leads to numerical difficulties because of the change in the track configuration. Singular configurations can be avoided by repeated changes in the recursive equations. The sensitivity of the predictor-corrector multistep numerical integration schemes to the method of formulating the state equations is demonstrated. The numerical results presented in this investigation are obtained using a planner tracked vehicle model that consists of fifty four rigid bodies.     

Plane Problem of Three-Dimensional Stability for a Hinged Plate with Two Symmetric End Cracks

The plane stability problem for a rectangular plate with two symmetric end cracks is solved in three-dimensional formulation. The three-dimensional linearized theory of stability and the finite-difference method are used. The effect of the crack parameter on the critical load is examined__________Translated from Prikladnaya Mekhanika, Vol. 41, No. 4, pp. 47–52, April 2005.     

On the finite displacement analysis of quadrangular plates

In this paper, a numerical method for the linear and geometrically non-linear static analysis of thin plates is presented. The method begins with the elasticity equations pertaining to strain components, stresses, displacement components, strain energy and work due to externally applied loads. The plate geometry is defined by a quadrangular boundary with four straight edges and the natural coordinates in conjunction with the Cartesian coordinates are used to map the geometry. The matrix equation of equilibrium is derived using the work-energy principle with the displacement fields expressed by algebraic polynomials, the coefficients of which are then manipulated to satisfy the kinematic boundary conditions. To validate the results from the present method, square plates having all sides fully fixed and all sides simply supported without in-plane movement are analysed. Comparison is made for the uniformly loaded square plate with the results obtained by Levy who solved the non-linear plate bending problem using the Th.von Ka’rma’n’s equations. Rhombic plates are examined and numerical results corresponding to these cases are presented in this paper. Very good comparison of the results regarding deflection and bending stresses with other sources available in the literature is found.     

充液系统液体-多体耦合动力响应分析

提出了充液系统的液体-多体耦合力学模型，基于ALE有限元法和多体系统动力学理论，发展了液体-多体耦合动力响应分析的一种有效方法. 对于液体子系统，将其运动分解为随同贮箱的大位移运动和相对贮箱的大幅晃动，引入贮箱固连参考系中的任意拉格朗日-欧拉(ALE)运动学描述，建立了贮箱固连非惯性参考系中液体的ALE有限元方程，对液体有限元方程的缩聚大大减少了液体子系统的计算规模. 为了计及液体阻尼的影响，引入了液体修正的Rayleigh阻尼，避免了伪阻尼力的出现. 对于多体子系统，应用多体系统动力学理论建立动力学方程. 在此基础上详细导出了液体-多体耦合动力学方程，并采用预估-多重校正算法(PMA)和时间步长控制算法进行迭代求解，既保证了迭代收敛，又提高了计算效率. 所给算例成功求解了液体运输车辆系统的液体-多体耦合动力响应，深入分析了有关参数对系统动力响应的影响，获得了一些结论.     

Nodeless variable finite element method for heat transfer analysis by means of flux-based formulation and mesh adaptation

Based on flux-based formulation, a nodeless variable element method is developed to analyze two-dimensional steady-state and transient heat transfer problems. The nodeless variable element employs quadratic interpolation functions to provide higher solution accuracy without necessity to actually generate additional nodes. The flux-based formulation is applied to reduce the complexity in deriving the finite element equations as compared to the conventional finite element method. The solution accuracy is further improved by implementing an adaptive meshing technique to generate finite element mesh that can adapt and move along corresponding to the solution behavior. The technique generates small elements in the regions of steep solution gradients to provide accurate solution, and meanwhile it generates larger elements in the other regions where the solution gradients are slight to reduce the computational time and the computer memory. The effectiveness of the combined procedure is demonstrated by heat transfer problems that have exact solutions. These problems are: (a) a steady-state heat conduction analysis in a square plate subjected to a highly localized surface heating, and (b) a transient heat conduction analysis in a long plate subjected to a moving heat source. The English text was polished by Yunming Chen.     

Time-harmonic BEM for 2-D piezoelectricity applied to eigenvalue problems

We will derive the fundamental generalized displacement solution, using the Radon transform, and present the direct formulation of the time-harmonic boundary element method (BEM) for the two-dimensional general piezoelectric solids. The fundamental solution consists of the static singular and the dynamics regular parts; the former, evaluated analytically, is the fundamental solution for the static problem and the latter is given by a line integral along the unit circle. The static BEM is a component of the time-harmonic BEM, which is formulated following the physical interpretation of Somigliana’s identity in terms of the fundamental generalized line force and dislocation solutions obtained through the Stroh–Lekhnitskii (SL) formalism. The time-harmonic BEM is obtained by adding the boundary integrals for the dynamic regular part which, from the original double integral representation over the boundary element and the unit circle, are reduced to simple line integrals along the unit circle.The BEM will be applied to the determination of the eigen frequencies of piezoelectric resonators. The eigenvalue problem deals with full non-symmetric complex-valued matrices whose components depend non-linearly on the frequency. A comparative study will be made of non-linear eigenvalue solvers: QZ algorithm and the implicitly restarted Arnoldi method (IRAM). The FEM results whose accuracy is well established serve as the basis of the comparison. It is found that the IRAM is faster and has more control over the solution procedure than the QZ algorithm. The use of the time-harmonic fundamental solution provides a clean boundary only formulation of the BEM and, when applied to the eigenvalue problems with IRAM, provides eigen frequencies accurate enough to be used for industrial applications. It supersedes the dual reciprocity BEM and challenges to replace the FEM designed for the eigenvalue problems for piezoelectricity.     

新型Poisson括号意义下的无穷维Lie代数

本文首先针对ＫｄＶ方程的Ｈａｍｉｌｔｏｎ形式，建立一种比较容易验证的新型Ｐｏｉｓｓｏｎ括号和无穷维Ｌｉｅ代数．其次，研究ＫｄＶ方程的Ｈａｍｉｌｔｏｎ形式的第一积分与新型Ｐｏｉｓｏｎ括号的关系，得到判定第一积分的充分必要条件．最后，构造ＫｄＶ方程的第一积分．     

DYNAMICS ANALYSIS OF STOCHASTIC SPATIAL FLEXIBLE MULTIBODY SYSTEM 1)

轻质、高精度的柔性多体系统被广泛应用于实际工程领域中.由于实际设计公差、制造误差及环境温度等多种不确定因素的存在,使得柔性多体系统的结构参数(物理参数和几何参数)表现出随机性.具有随机结构参数的动力学模型能够客观地反映出真实系统的动力学行为,且结构参数的不确定性对空间柔性多体系统动力学响应的影响是不容忽视的.针对具有多个随机参数的空间柔性多体系统,提出了一种基于广义alpha算法的非侵入式随机柔性多体系统动力学计算方法.采用绝对节点坐标公式(absolute node coordinate formulation, ANCF)来描述柔性体, 推导建立多体系统动力学模型.利用混沌多项式展开(polynomial chaos expansion, PCE)法构建系统随机动力学方程的代理模型,然后将随机响应面法(stochastic response surface method, SRSM)嵌入广义-alpha方法中,分别采用改进抽样的回归方法(regression method of improved sampling, RMIS)和单项求容积法则(Monte Carlo simulation, MCR)来确定样本点.将数值计算结果与蒙特卡洛模拟(Monte Carlo simulation, MCS)结果进行对比, 验证了所提算法的有效性.在相同的定积分精度的条件下,根据单项求容积法则确定的样本点的计算结果稳定性更强, 且其计算效率更高.     

成型充填过程的ALE有限元模拟

在ALE框架中提出了一个用于成型充填过程有限元数值模拟的模型。应用ALE参考构形及ALE参考粒子速度描写充填过程中的熔体质量运动。摒弃了Hele-Shaw近似假定,因而所提出的模型能用于非薄壁型腔中高分子材料充填过程的数值模拟。应用基于时域分步算法的Taylor-Galerkin方法,对控制成型充填过程的守恒方程建立了弱形式。对移动自由面附近的充填材料区构造了网格生成算法与网格重划分方案。给出了在几种不同形状的典型腔体中充填过程的数值模拟结果,表明了所提出的ALE有限元模型模拟充填过程的有效性。