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.     

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.     

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.     

Diffusional growth of cloud particles: existence and uniqueness of solutions

Diffusional growth of cloud particles is commonly described by a coupled system of parabolic equations and ordinary differential equations. The Dirichlet boundary condition for the parabolic equation is obtained from the solution of the ordinary differential equations, but this solution itself depends on the solution of the parabolic equations. We first present the governing equations describing diffusional growth of cloud particles. In a second step, we consider a simplified model problem, motivated by the diffusional growth equations. The main difference between the simplified model problem and the diffusional growth equations consists in neglecting the dependence of the domain for the parabolic equations on the solution. For the model problem, we show unique solvability using a fixed point method. Finally, we discuss application of the main result for the model problem to the diffusional growth equations and illustrate these equations with the help of a numerical solution.     

Equations of cylindrical bending of orthotropic plates with arbitrary conditions on their front surfaces

Based on approximations of solutions of elasticity theory equations by Legendre polynomial segments, differential equations for bending of orthotropic plates are constructed. In contrast to equations constructed with the use of kinematic and force hypotheses, the order of these differential equations is independent of the type of conditions on front surfaces. The matrices of the constructed equations depend on the type of boundary conditions. An analytical solution is given for the system of equations in the case with normal and shear stresses being specified on the upper and lower front surfaces.     

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

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

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.     

A tensor method for the derivation of the equations of rigid body dynamics

A tensor method for the derivation of the equations of rigid body dynamics,based onthe concepts of continuum mechanics,is presented.The formula of time derivative of theinertia tensor with zero corotational rate is used to prove the equivalences of five methods,namely,Lagrange’s equations,Nielsen’s equations,Gibbs-Appell’s equations,Kane’sequations and the generalized momentum type of Kane’s equations.Some differentialidentities on angular velocity and angular acceleration are given.     

Oscillatory properties of the solutions of nonlinear delay hyperbolic differential equations of neutral type

Ｔｈｅｎｅｗｌｙ—ｄｅｖｅｌｏｐｅｄｔｈｅｏｒｙｏｆｎｏｎｌｉｎｅａｒｄｅｌａｙｐａｒｔｉａｌｆｕｎｃｔｉｏｎａｌｄｉｆｆｅｒｅｎｔｉａｌｅｑｕａｔｉｏｎｓａｒｅａｐｐｌｉｅｄｉｎｍａｎｙｆｉｅｌｄｓ,ｓｕｃｈａｓｉｎｅｎｇｉｎｅｅｒｉｎｇ ,ｂｉｏｌｏｇｙ,ｍｅｄｉｃｉｎｅ,ｐｈｙｓｉｃｓａｎｄｃｈｅｍｉｓｔｒｙ .Ａｓｅｒｉｅｓｏｆｓｕｆｆｉｃｉｅｎｔｃｏｎｄｉｔｉｏｎｓ,ｎｅｃｅｓｓａｒｙａｎｄｓｕｆｆｉｃｉｅｎｔｃｏｎｄｉｔｉｏｎｓｆｏｒｏｓｃｉｌｌａｔｉｏｎｓｏｆｔｈｅｅｑｕａｔｉｏｎｓｗｅｒ…     

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     

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 …     

Modified equations of finite-size layered plates made of orthotropic material. Comparison of the results of numerical calculations with analytical solutions

This paper describes the modified bending equations of layered orthotropic plates in the first approximation. The approximation of the solution of the equation of the three-dimensional theory of elasticity by the Legendre polynomial segments is used to obtain differential equations of the elastic layer. For the approximation of equilibrium equations and boundary conditions of three-dimensional theory of elasticity, several approximations of each desired function (stresses and displacements) are used. The stresses at the internal points of the plate are determined from the defining equations for the orthotropic material, averaged with respect to the plate thickness. The construction of the bending equations of layered plates for each layer is carried out with the help of the elastic layer equations and the conjugation conditions on the boundaries between layers, which are conditions for the continuity of normal stresses and displacements. The numerical solution of the problem of bending of the rectangular layered plate obtained with the help of modified equations is compared with an analytical solution. It is determined that the maximum error in determining the stresses does not exceed 3 %.     

Nonlinear vibration of corrugated shallow shells under uniform load

Based on the large deflection dynamic equations of axisymmetric shallow shells of revolution, the nonlinear forced vibration of a corrugated shallow shell under uniform load is investigated. The nonlinear partial differential equations of shallow shell are reduced to the nonlinear integral-differential equations by the method of Green's function. To solve the integral-differential equations, expansion method is used to obtain Green's function. Then the integral-differential equations are reduced to the form with degenerate core by expanding Green's function as series of characteristic function. Therefore, the integral-differential equations become nonlinear ordinary differential equations with regard to time. The amplitude-frequency response under harmonic force is obtained by considering single mode vibration. As a numerical example, forced vibration phenomena of shallow spherical shells with sinusoidal corrugation are studied. The obtained solutions are available for reference to design of corrugated shells     

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.     

Invariants of a Remarkable Family of Nonlinear Equations

In classical literature, invariants of families of differentialequations were considered for linear equations only, e.g. the renownedLaplace invariants for linear hyperbolic partial differential equationsand invariants of linear ordinary differential equations with variablecoefficients. The restriction to linear equations was essential inpioneering works of Cockle, Laguerre, Halphen, andForsyth for tackling the problem of invariants of differentialequations. Lie regretted that these authors did not use advantagesprovided by his theory of infinite continuous groups, but he himself didnot undertake further developments in this direction.Recently, the present author considered the possibility hinted byLie's remark and introduced the infinitesimal technique in thetheory of invariants of families of differential equations thatwas lacking in old methods. In consequence, a simple unifiedapproach was developed for calculation of invariants of algebraicand differential equations independent on the assumption oflinearity of equations. It was employed recently for calculationof Laplace type invariants for parabolic equations. Here, themethod is applied to calculation of invariants for the family ofnonlinear equations appearing in the problem on linearization ofnonlinear ordinary differential equations.     

The Effective Elastic Constants of Porous Polycrystals of Trigonal Symmetry

A theory behind the effective properties of piezoelectric crystals is considered. It is based on the stochastic electroelastic equations of a microinhomogeneous medium. The initial equations are reduced to integral equations, which, after conditional averaging, are reduced to a system of linear algebraic equations in component-average parameters. The effective electroelastic properties of porous polycrystals with trigonal symmetry are evaluated     

半圆环截面细环壳的屈曲

本文由Sanders非线性平衡方程和Koiter小应变协调方程推导出细环壳的非线性微分方程和稳定方程。用伽辽金法求解了静水压或边界载荷作用下的半园环截面细环壳的稳定方程。对于不同的边界条件及一系列几何参数,计算得到了临界载荷及屈曲模态。     

分数因子与分数阶完整力学系统的运动方程和循环积分

引入分数因子和分数增量,给出了分数阶微积分的定义和性质;基于分数阶导数的定义,证明了含有分数因子的等时变分与分数阶算子的交换关系;提出了分数阶完整保守和非保守系统的Hamilton原理;建立了分数阶完整保守系统和非保守系统的运动微分方程;得到了分数阶完整保守系统的循环积分;并利用分数阶循环积分导出分数阶罗兹方程．最后给出了两个例子．研究表明利用分数因子给出的分数阶微分方程是一个含有分数因子的通常的微分方程,那么分数阶系统运动微分方程的求解都可以采用通常微分方程的求解方法．     

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

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

层合圆柱三维温度场分析的半解析-精细积分法

针对圆柱体的三维温度场分析,提出了一种高效的半解析-精细积分法。将温度场展开为环向坐标的Fourier级数,并对径向坐标进行差分离散,从而把三维热传导方程简化为一系列二阶常微分方程;将这些二阶常微分方程转化为哈密顿体系下的一阶状态方程,并利用两点边值问题的精细积分法求解。由于该方法仅对径向坐标进行差分离散,故相对于传统的数值方法离散规模大幅度减少,不仅提高了计算效率、降低了存贮量,而且缓解了代数方程的病态问题。此外,针对Fourier半解析解,根据热平衡原理推导出了两种材料衔接面的半解析差分方程,从而为求解复合材料层合柱问题打下了基础。算例结果表明,即使对于细长比高达400的圆柱杆件,此方法仍然可以给出精度较高的解答。