The equations of motion of a system under the action of the impulsive constraints

In order to solve the problem of motion for the system with n degrees of freedom under the action of p impulsive constraints, we must solve the simultaneous equations consisting of n+p equations. In this paper, it has been shown that the undetermined multipliers in the equations of impact can be cancelled for the cases of both the generalized coordinates and the quasi-coordinates. Thus there are only n-p equations of impact. Combining these equations with p impulsive constraint equations, we have simultaneous equations consisting of n equations. Therefore, only n equations are necessary to solve the problem of impact for the system subjected to impulsive constraints. The method proposed in this paper is simpler than ordinary methods.     

Multigrid solution of compressible turbulent flow on unstructured meshes using a two-equation model

The steady state solution of the system of equations consisting of the full Navier-Stokes equations and two turbulence equations has been obtained using a multigrid strategy on unstructured meshes. The flow equations and turbulence equations are solved in a loosely coupled manner. The flow equations are advanced in time using a multistage Runge-Kutta time-stepping scheme with a stability-bound local time step, while the turbulence equations are advanced in a point-implicit scheme with a time step which guarantees stability and positivity. Low-Reynolds-number modifications to the original two-equation model are incorporated in a manner which results in well-behaved equations for arbitrarily small wall distances. A variety of aerodynamic flows are solved, initializing all quantities with uniform freestream values. Rapid and uniform convergence rates for the flow and turbulence equations are observed.     

考虑质量耦合效应的流体饱和孔隙介质波动方程

利用混合物理论和连续介质力学的基本原理，推导了考虑质量耦合效应的流体饱和弹性孔隙介质的波动方程，并与经典的Biot波动方程进行了对比．结果表明：该文得到的方程包含了Biot波动方程的所有要素，且形式与后者基本相同．比较而言，该文推导过程具有更明确的物理意义，概念也更完整．     

多体系统指标2运动方程HHT方法违约校正1）

采用Cartesian绝对坐标建模方法,完整约束多体系统运动方程是指标3的微分--代数方程（differentialalgebraic equations,DAEs）,数值求解指标3的DAEs属于高指标问题,通过对位置约束方程求导,可使运动方程的指标降为2.位置约束方程求导得到的是速度约束方程.直接求解指标3的运动方程,速度约束方程得不到满足,而且高指标DAEs的数值求解存在一些问题.论文首先采用HHT（Hilber--Hughes--Taylor）直接积分方法求解降指标得到的指标2运动方程,此时速度约束方程参与离散计算,从机器精度上讲速度约束自然得到满足,而位置约束方程没有参与计算,存在“违约”.针对违约问题,采用基于Moore--Penrose广义逆理论的违约校正方法,消除位置约束方程的违约.指标2运动方程HHT方法违约校正,将HHT方法和违约校正方法很好地结合,在数值求解指标2运动方程的过程中,位置约束方程和速度约束方程都不存在违约问题,而且新方法没有引入新的未知数向量,离散得到的非线性方程组的方程数量与原指标2运动方程的方程数量相同,求解规模没有扩大.新方法的实用和有效性通过算例的数值实验得到验证,数值实验也说明新方法保持了HHT方法本身具有的数值阻尼可以控制和二阶精度的特性.最后从非线性方程组的求解规模和计算速度上与其他方法进行了比较分析,说明新方法的优势所在.     

高雷诺数流动的基本控制方程体系和扩散跑物化Navier-Stokes方程组的意义和用途

在计算机发达的时代, 高雷诺($Re$)数绕流计算中有无必要使用简化NS方程组, 本文讨论这个问题. 主要内容如下: (1)高$Re$数绕流包含3种基本流动: 所有方向对流占优流动、所有方向对流扩散竞争流动和部分方向对流占优部分方向对流扩散竞争流动(简称干扰剪切流动), 3个基本流动的特征彼此不同且在流场中所占领域大小彼此相差悬殊, NS方程区域很小,它们的最简单控制方程组Euler、Navier-Stokes (NS)和扩散抛物化(DP) NS方程组的数学性质彼此不同, 因此利用Euler-DPNS-NS方程组体系分析计算高$Re$数绕流流动就是一个合乎逻辑的选择, 该法与利用单一NS方程组的常用方法可以彼此检验和补充. (2)流体之间以及流体与外界的动量、能量和质量交换, 流态从层流到湍流的演化主要发生在干扰剪切流动中, 干扰剪切流及其最简单控制方程------DPNS方程组具有基础意义; DPNS方程组笔者在1967年已提出. (3)诸简化NS方程组: DPNS、抛物化(P)NS、薄层(TL)NS、黏性层(VL)NS方程组的发展、相互关系, 它们的历史贡献和今后的用途; 它们的数学性质均为扩散抛物型, 但它们包含的黏性项彼此有所不同; 从流体力学角度来看, 它们中只有DPNS方程组能够准确描述干扰剪切流动. 提出把诸简化NS方程组统一为DPNS方程组的建议. (4)干扰剪切流------DPNS方程组与无干扰剪切流------边界层方程组之间的关系以及进一步研究干扰剪切流的意义.      

Unsymmetrical nonlinear bending problem of circular thin plate with variable thickness

Introduction Manystructuralelements(pole,plate,shell)withunevenandvariablethicknessarewidely usedinallkindsofengineeringfields.Engineerscansavematerialswhentheydesignbecause theseelementshavebetteroptimizedshapeofstructuralfeature,butthisaddsdifficultytotheanalysisoftheirmechanicalcapability.Manypreviouspapers[1-4]havesolvedtheproblemof symmetricalaxis,butnobodyhassolvedtheunsymmetricalnonlineardeformationproblemof circularthinplatewithvariablethicknessandunsymmetricalaxisuptonow,afewworkonly …     

多体动力学超定运动方程广义-α求解新算法

完整约束多体系统第一类Lagrange方程建模得到的运动方程是指标-3形式的微分-代数方程(differental-algebraic equations,DAEs).如果同时考虑速度约束,将得到超定运动方程,该方程是指标-2的超定微分-代数方程(over-determined differential-algebraic equations,ODAEs).基于结构动力学中常用的广义-α方法,将其拓展,求解包含速度约束的超定运动方程,相对于其他求解指标-2 ODAEs的算法,新的算法没有增加离散得到的非线性方程组方程的数目.通过数值实验验证算法,并说明其求解ODAEs不存在精度降阶的现象,仍然具有二阶精度,同时算法的数值耗散也是可以控制的.最后新方法与其他求解多体系统ODAEs形式运动方程算法的CPU时间进行了比较分析.     

多体系统动力学方程违约修正的数值计算方法

多体系统动力学方程为微分代数方程,一般将其转化成常微分方程组进行数值计算,在数值积分的过程中约束方程的违约会逐渐增大.本文对具有完整、定常约束的多体系统,在修改的带乘子Lagrange正则形式的方程的基础上,根据Baumgarte提出的违约修正的方法,给出了一种多体系统微分代数方程违约修正法和系统的动力学方程的矩阵表达式.通过对曲柄-滑块机构的数值仿真,计算结果表明本文给出的方法在计算精度和计算效率上好于Baumgarte提出的两种违约修正的方法.     

A new multigrid algorithm for non-linear equations in conjunction with time-marching procedures

Full approximate storage (FAS) multigrid algorithm is the most commonly used multigrid algorithm for non-linear equations. The algorithm initially developed for steady-state equations was later extended to obtain steady-state solutions employing unsteady equations. In extending the FAS algorithm for the steady-state non-linear equations to unsteady non-linear equations, the FAS algorithm does not to take into account that the governing equations are typically linearized in time before they are solved. Thus, there is a scope to develop a new multigrid algorithm to apply the multigrid technique to the equations linearized in time. In the present work, such an algorithm is developed exploring this possibility and is implemented for two-dimensional incompressible and compressible flows coupled with explicit time marching procedures. The results of the new algorithm compare favourably with those of the FAS multigrid method and single grid.     

Another form of equations of motion for constrained multibody systems

This paper presents a new and simplified set of explicit equations of motion for constrained mechanical systems. The equations are applicable with both holonomic and nonholonomic systems and the constraints may, or may not, be ideal. It is shown that this set of equations is equivalent to governing equations developed earlier by others. The connection of these equations with Kane's equations is discussed. It is shown that the developed equations are directly applicable with controlled systems where the controlling forces and moments may be subject to constraints. Finally, a procedure is presented for determining which control force systems are equivalent. Examples are presented to demonstrate the advantages, features, and range of application of the equations.     

Representation of solutions to equations of hyperbolic type

The general solutions to hyperbolic equations of fourth and sixth order are obtained using Vekua’s method for the representation of the general solutions to elliptic equations of order 2n with the aid of n analytic functions. It is assumed that the right-hand sides of the hyperbolic equations can be expanded in time series of sines. The systems of equations of various approximations for a prismatic thin body in terms of moments with respect to the system of Legendre polynomials can be reduced to these equations and to some hyperbolic-type equations of higher order.     

Motion of strings, embedding problem and soliton equations

The motion of a flexible string of constant length in E ³ in interaction, corresponding to a variety of physical situations, is considered. It is pointed out that such a system could be studied in terms of the embedding problem in differential geometry, either as a moving helical space curve in E ³ or by the embedding equations of two dimensional surfaces in E ³. The resulting integrability equations are identifiable with standard soliton equations such as the non-linear Schrödinger, modified K-dV, sine-Gordon, Lund-Regge equations, etc. On appropriate reductions the embedding equations in conjunction with suitable local space-time and/or gauge symmetries reproduce the AKNS-type eigenvalue equations and Riccati equations associated with soliton equations. The group theoretical properties follow naturally from these studies. Thus the above procedure gives a simple geometric interpretation to a large class of the soliton possessing nonlinear evolution equations and at the same time solves the underlying string equations.     

Burnett simulations of gas flow in microchannels

The Burnett equations with slip boundary conditions are used to model the gas flow in microchannels in transition flow regime. As the Navier-Stokes equations are not appropriate to model the gas flow in this regime, the higher-order Burnett equations are adopted in the present study. In earlier studies, convergent solutions of the Burnett equations of microPoiseuille flow could only be obtained when Knudsen number is less than 0.2. By using a relaxation method on the boundary values, convergent solutions of the Burnett equations can be obtained even when Knudsen number reaches 0.4. The solutions of Burnett equations agree very well with experimental data and direct simulation Monte Carlo (DSMC) results. The pressure distributions and velocity profiles are then discussed in detail.     

Dynamic simulation of the hydrodynamic interaction among immersed particles in Stokes flow

A numerical method for the dynamic simulation of the hydrodynamic interaction among particles in Stokes flow is developed. The method couples the quasi-static Stokes equations for the fluid with the equilibrium equations for the particles. The boundary element method is used to represent the velocity at a general field point in terms of surface velocities and stresses. However, neither the stresses nor the velocities are assumed to be known on the surface of the particles. Kinematic equations relating the linear and angular velocities at the centroids of the particles to the surface velocities are combined with the discretized boundary element equations and the equilibrium equations to generate a system of linear equations. The associated coefficient matrix is correspondent to the grand resistance matrix which relates the velocities of the particles to a given geometry.     

带约束非线性多体系统动力学方程数值分析方法

Lagrange方法是建立带约束多体系统动力学方程的普遍方法之一 ,其方程的形式为微分 代数方程组 ,数值计算与数值分析是研究多体系统动力学特性的重要方法。本文利用缩并法给出了带约束多体系统动力学方程的隐式数值计算方法和Lyapunov指数的计算方法。将数值仿真、Lya punov指数计算和Poincare映射有机结合 ,分析非线性多体系统动力学行为。通过一个算例 ,说明该方法的有效性     

Solutions of matrix equations in problems of mechanics and control

The paper considers algorithms for solving linear matrix equations related to problems of mechanics and control, namely, the Lyapunov and Sylvester matrix equations and Riccati-type nonlinear matrix equations. These algorithms are capable of solving both linear equations and linear matrix inequalities. Algorithms based on the Bass relations are used to solve Riccati-type nonlinear matrix equations in so-called special cases where some eigenvalues of the matrix pencil are on a unit circle. These algorithms are compared with those of other authors by way of examples. It is shown that the algorithms can be implemented in symbolic computing routings, which allows solving these equations with high accuracy     

多体系统Lagrange方程数值算法的研究进展

Lagrange方法是建立多体系统动力学方程的普遍方法之一, 其方程的形式为常微分方程组或微分-代数方程组,数值计算与数 值分析是研究多体系统动力学特性的重要方法。本文简要介绍了多 体系统动 力学方程的第一、二类Lagrange方程和修正的Lagrange方 程的基本形式及这些方程的正则形式,着重介绍了正则方程在数值 计算中的特点,就多体系统Lagrange方程的隐式算法、辛算法和多 体系统动力学特性的数值分析方法(包括数值仿真、 Poincarè映射 和Lyapunov指数的计算方法)的研究现状进行了综述。     

双面约束多点摩擦多体系统的建模和数值方法

提出了一种建立具有固定双面约束多点摩擦的多体系统动力学方程的方法. 用笛卡尔坐标阵描述系统的位形,根据局部方法的递推关系建立系统的约束方程,应用第一类Lagrange方程建立该系统的动力学方程,使得具有摩擦的约束面的法向力与Lagrange乘子一一对应,便于摩擦力的分析与计算,并用矩阵形式给出了摩擦力的广义力的一般表达式. 应用增广法将微分-代数方程组转化为常微分方程组,并用分块矩阵的形式给出,以便于方程的编程与计算.给出了一种改进的试算法,可提高计算效率. 最后给出了一个算例,应用试算法和RK法对算例进行了数值仿真.      

Bulk equations and Knudsen layers for the regularized 13 moment equations

The order of magnitude method offers an alternative to the Chapman-Enskog and Grad methods to derive macroscopic transport equations for rarefied gas flows. This method yields the regularized 13 moment equations (R13) and a generalization of Grad’s 13 moment equations for non-Maxwellian molecules. Both sets of equations are presented and discussed. Solutions of these systems of equations are considered for steady state Couette flow. The order of magnitude method is used to further reduce the generalized Grad equations to the non-linear bulk equations, which are of second order in the Knudsen number. Knudsen layers result from the linearized R13 equations, which are of the third order. Superpositions of bulk solutions and Knudsen layers show good agreement with DSMC calculations for Knudsen numbers up to 0.5.      

基于群桩振动方程组的群桩基础阻抗函数求解及比较

本文首先在位移场叠加的基础上,对经典的单桩振动方程进行了修正得到群桩振动方程组;然后根据群桩振动方程组的近似求解,导出了群桩阻抗函数求解的常用方法一群桩阻抗方程方法。本文中给出的群桩阻抗方程方法在具体的表达形式上与常见的阻抗方程计算方法略有差别;而后根据振动方程组系数矩阵在复空间为正规阵的性质,提出了该四阶耦合微分方程组解耦求解的详细步骤,在直接求解微分方程组的基础上给出了群桩基础阻抗计算的一种新方法。与阻抗方程方法相比较,新方法未采用柔度因子来考虑桩与桩之间的相互影响,而是通过耦合的运动方程组对桩与桩之间的相互作用作了新的表述,并且近似地考虑了屏蔽效应。本文最后给出了上述两种计算方法的数值算例。