State vector solutions for nonaxisymmetric problem of multilayered half space piezoelectric medium

A state space formulation is established for the nonaxisymmetric space problem of transversely isotropic piezoelectric media in a system of cylindrical coordinate by introducing the state vector. Using the Hankel transform and the Fourier series, the state vector equations are transformed into ordinary differential equations. By the use of the matrix methods, the analytical solutions of a single piezoelectric layer are presented in the form of the product of initial state variables and transfer matrix. The applications of state vector solutions are discussed. An analytical solution for a semiinfinite piezoelectric medium subjected to the vertical point forceP _z, horizontal point forceP _x along x-direction and point electric charge Q at the origin of the surface is presented. According to the continuity conditions at the interfaces, the general solution formulation forN-layered transversely isotropic piezoelectric media is given. Project supported by the National Natural Science Foundation of China (Grant No. 59648001).     

Model of motion of a controlled body in the base attached trihedron

The position of a controlled body in space is described in a moving Cartesian system of coordinates, which we call the base attached trihedron. A model of motion of the body in these axes is derived in the form of a system of differential equations allowing for the Earth's rotation and the variable curvature of the Earth's surface.Translated from Vychislitel'naya i Prikladnaya Matematika, No. 61, pp. 112–118, 1987.     

Systematic dynamic modelling of mechanical systems containing kinematic loops

In this tutorial paper a systematic procedure is presented to obtain the dynamic models of mechanical systems containing kinematic loops, with a main emphasis on efficiency and with particular regard to robotic systems. The procedure retains a minimal set of generalized coordinates for the corresponding open loop structure, obtained by removing some additional constraints closing loops in the original structure, while adopting an efficient Newton-Euler formulation of the equations of motion. Two methods for including the loop closure equations are then discussed: the Lagrange multipliers method and the method based on an explicit solution of the constraint equations. In the first case the dynamic model is given in the form of an implicit Differential Algebraic Equations (DAE) system, while in the second case an Ordinary Differential Equations (ODE) system could be obtained.     

Expansion of solution of an inverted pendulum system with time delay

In this paper, the solution expansion for an inverted pendulum system with time delay is studied. The linearized model of this nonlinear system near its equilibrium is derived on the assumption that a unique equilibrium exists in it. Then the asymptotic expressions of its eigenvalues and the eigenvalues’ corresponding eigenvectors are obtained. Moreover, although the set of these eigenvectors does not form a Schauder basis for the state space, the solution of this model still can be expressed by these eigenvectors in the form of infinite series under certain conditions. Finally, a simulation is provided to support these results.     

From Generalized Clifford Algebras to nambu’s formulation of dynamics

We consider a commutative part of the Generalized Clifford Algebras, denominated asalgebra of multicomplex numbers. By using the multicomplex algebra as underlying algebraic structure we construct oscillator model for the Nambu’s formulation of dynamics. We propose a new dynamicals principle which gives rise to two kinds of Hamilton-Nambu equations inD≥2-dimensional phase space. The first one is formulated with (D−1)-evolution parameter and a single Hamiltonian. The Haniltonian of the oscillator model in such dynamics is given byD-degree homogeneous form. In the second formulation, vice versa, the evolution of the system along a single evolution parameter is generated by (D−1) Hamiltonian.     

A two-grid algorithm based on Newton iteration for the stream function form of the Navier-Stokes equations

In this paper, we propose a two-grid algorithm for solving the stream function formulation of the stationary Navier-Stokes equations. The algorithm is constructed by reducing the original system to one small, nonlinear system on the coarse mesh space and two similar linear systems (with same stiffness matrix but different right-hand side) on the fine mesh space. The convergence analysis and error estimation of the algorithm are given for the case of conforming elements. Furthermore, the algorithm produces a numerical solution with the optimal asymptotic H ²-error. Finally, we give a numerical illustration to demonstrate the effectiveness of the two-grid algorithm for solving the Navier-Stokes equations.     

A clustering approach to the planar hub location problem

The problem tackled in this paper is as follows: consider a set ofn interacting points in a two-dimensional space. The levels of interactions between the observations are given exogenously. It is required to cluster then observations intop groups, so that the sum of squared deviations from the cluster means is as small as possible. Further, assume that the cluster means are adjusted to reflect the interaction between the entities. (It is this latter consideration which makes the problem interesting.) A useful property of the problem is that the use of a squared distance term yields a linear system of equations for the coordinates of the cluster centroids. These equations are derived and solved repeatedly for a given set of cluster allocations. A sequential reallocation of the observations between the clusters is then performed. One possible application of this problem is to the planar hub location problem, where the interacting observations are a system of cities and the interaction effects represent the levels of flow or movement between the entities. The planar hub location problem has been limited so far to problems with fewer than 100 nodes. The use of the squared distance formulation, and a powerful supercomputer (Cray Y-MP) has enabled quick solution of large systems with 250 points and four groups. The paper includes both small illustrative examples and computational results using systems with up to 500 observations and 9 clusters.     

Numerical integration of differential-algebraic equations with mixed holonomic and control constraints

The present work aims at the incorporation of control (or servo) constraints into finite–dimensional mechanical systems subject to holonomic constraints. In particular, we focus on underactuated systems, defined as systems in which the number of degrees of freedom exceeds the number of inputs. The corresponding equations of motion can be written in the form of differential–algebraic equations (DAEs) with a mixed set of holonomic and control constraints. Apart from closed–loop multibody systems, the present formulation accommodates the so–called rotationless formulation of multibody dynamics. To this end, we apply a specific projection method to the DAEs in terms of redundant coordinates. A similar projection approach has been previously developed in the framework of generalized coordinates by Blajer & Kołodziejczyk [1]. (© 2008 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

On boundary-volume element discretization of inhomogenous elastodynamic problems

This paper presents an integral equation formulation and its discretization scheme for the elastodynamic problem in which the material properties are prescribed as arbitrary, continuous and differentiable functions of the spatial coordinates. The formulation is made by using the Green's function for the corresponding problem in homogenous elasticity. From a weighted residual statement of the problem, the governing differential equation is transformed into a set of the integral equations in the inner domain as well as on the boundary. These integral equations are discretized by introducing a finite number of the boundary-volume-time elements, and the solution for the system of linear equations thus obtained is discussed.     

Inverse dynamics of underactuated multibody systems

The present work deals with controlled mechanical systems subject to holonomic constraints. In particular, we focus on underactuated systems, defined as systems in which the number of degrees of freedom exceeds the number of inputs. The governing equations of motion can be written in the form of differential-algebraic equations (DAEs) with a mixed set of holonomic and control constraints. The rotationless formulation of multibody dynamics will be considered [1]. To this end, we apply a specific projection method to the DAEs in terms of redundant coordinates. A similar projection approach has been previously developed in the framework of generalized coordinates by Blajer & Kołodziejczyk [2]. (© 2009 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Finding strong defining hyperplanes of PPS using multiplier form

The production possibility set (PPS) is defined as the set of all inputs and outputs of a system in which inputs can produce outputs. In this paper, we deal with the problem of finding the strong defining hyperplanes of the PPS. These hyperplanes are equations that form efficient surfaces. It is well known that the optimal solutions of the envelopment formulation for extreme efficient units are often highly degenerate and, therefore, may have alternate optima for the multiplier form. Every optimal solution of the multiplier form yields a hyperplane which is supporting at the PPS. We will show that the hyperplane which corresponds to an extreme optimal solution of the multiplier form (in evaluating an efficient DMU), and whose components corresponding to inputs and outputs are non zero is a strong defining hyperplane of the PPS. This will be discussed in details in this paper. These hyperplanes are useful in sensitivity and stability analysis, the status of returns to scale of a DMU, incorporating performance into the efficient frontier analysis, and so on. Using numerical examples, we will demonstrate how to use the results.     

Controllability and stabilization of the programmed motions of a transport robot

A non-linear mathematical model for the motion of a transport robot (TR) with a caterpillar chassis and with drives based on DC motors, which is a non-holonomic electromechanical system, is considered. Non-linear canonical transformations of the coordinates of the state and control space are constructed, which reduce the initial equations of motion of the TR to a simpler canonical form, which is convenient for analysing and synthesizing control systems for the TR. The conditions for the TR to be controllable as a controlled object are found. Algorithms are given for constructing programmed motions (PMs) of the TR. Stabilizing control laws are synthesized under which the PMs of the TR are asymptotically stable and transients of a specified nature are ensured.     

A kinetic formulation for a model coupling free surface and pressurised flows in closed pipes

The aim of this paper is to present a kinetic formulation of a model for the coupling of transient free surface and pressurised flows. Firstly, we revisit the system of Saint-Venant equations for free surface flow: we state some properties of Saint-Venant equations, we propose a kinetic formulation and we verify that this kinetic formulation leads to a Gibbs equilibrium that minimises (in some general case) an energy and preserves the still water steady state. Secondly, we propose a model for pressurised flows in a Saint-Venant-like conservative formulation. We then propose a kinetic formulation and we verify that this kinetic formulation leads to a Gibbs equilibrium that minimises in any case an energy and preserves the still water steady state. Finally, we propose a dual model that couples these two types of flow.     

On the construction and investigation of solutions of boundary value problems of thermoviscoelasticity

The uncoupled mixed boundary value problem of thermoviscoelasticity is considered in a quasistatic formulation. The temperature distribution is assumed nonstationary and inhomogeneous. The influence of the temperature on the viscoelastic properties of the material is taken into account by the introduction of a reduced time. The equations of state of the material are written in differential form as a system of kinetic equations in some tensor-type strain parameters. The system mentioned is equivalent to a Volterra integral equation with kernel in the form of a sum of exponents. The differential approach used is apparently more convenient for numerical realization /1/ (especially in nonuniform problems) and results in a substantially different mathematical formulation as compared with that based on the integral form of writing the equations of state investigated in /2,3/. Precisely for going over to the boundary value problem are the kinetic differential equations converted into an operator differential equation in Hubert space. The existence, uniqueness, and stability of the solution of the problem formulated are established, and conditions for the convergence of the Galerkin approximations and the stability of the difference approximations in time are formulated.     

极坐标系下非分裂PML及时域有限元实现

在弹性波传播的数值模拟中,吸收边界被广泛应用于截取有限空间进行无限空间问题的分析.完全匹配层（perfect matched layer, PML）吸收边界较其他吸收边界条件具有更优越的吸收性能,已被成功应用于直角坐标系下的弹性波方程正演模拟.考虑极坐标系下二阶弹性波动方程,通过采用辅助函数的方法,提出了一种非分裂格式的完全匹配层吸收边界条件.并且基于Galerkin近似技术,给出了非对称以及轴对称条件下的时域有限元计算格式.通过数值算例分析了该极坐标系下分裂格式的完全匹配层吸收边界的有效性.     

A Primal-Dual Formulation for Certifiable Computations in Schubert Calculus

Formulating a Schubert problem as solutions to a system of equations in either Plücker space or local coordinates of a Schubert cell typically involves more equations than variables. We present a novel primal-dual formulation of any Schubert problem on a Grassmannian or flag manifold as a system of bilinear equations with the same number of equations as variables. This formulation enables numerical computations in the Schubert calculus to be certified using algorithms based on Smale’s \(\alpha \)-theory.     

An Optimal Control Problem for a Singular System of Solid–Liquid Phase Transition

In this article, we deal with a control problem for a singular system regarding a phase-field model which describes a solid–liquid transition by the Ginzburg–Landau theory. The purpose is to control the system by the means of the heat supply r able to guide it into a certain state with a solid (or liquid) part in a prescribed subset Ω₀ of the space domain Ω, and maintain it in this state during a period of time. The transition is described by a nonlinear differential system of two equations for the phase field and temperature. The control problem is set for some expressions of the cost functional which might reveal cases of physical interest. An approximating control problem is introduced and the existence of at least an optimal pair is proved. The first-order optimality conditions for the approximating problem are determined and a convergence result is given.