Classification of explicit three-stage symplectic difference schemes for the numerical solution of natural Hamiltonian systems: A comparative study of the accuracy of high-order schemes on molecular dynamics problems

The natural Hamiltonian systems (systems with separable Hamiltonians) are considered. The variety of explicit three-stage symplectic schemes is described. A classification of the third-order accurate schemes is given. All fourth-order schemes are found (there are seven of them). It is proved that there are no fifth-order schemes. The schemes with improved properties, such as invertibility and optimality with respect to the phase error, are listed. Numerical results that demonstrate the properties of these schemes are presented, and their comparative analysis with respect to the accuracy–efficiency criterion is given. The disbalance of total energy is used as the accuracy criterion.     

Translation Ovoids and Spreads of the Generalized Hexagon H(q)

Ovoids of the finite classical generalized hexagon H(q) that are translation with respect to a point are classified. By duality, translation spreads with respect to a line are classified when the characteristic is three. When the characteristic is not equal to three, it is shown that there are no ovoids that are translation with respect to a flag.     

Infectious disease models with time-varying parameters and general nonlinear incidence rate

Infectious disease models with time-varying parameters and general nonlinear incidence rates are analyzed. The functional form of the nonlinear incidence rate is assumed to change in time, due to, for example, environmental factors or a change in population behavior. More specifically, a new SIR model with time-varying parameters and switched nonlinear incidence rate is studied. The stability of the disease-free equilibrium is investigated, as well as disease persistence in the endemic case. A switched epidemic model with generalized compartments and time-varying parameters is also proposed and analyzed. Pulse vaccination and pulse treatment are applied to the new SIR model with seasonality and switched incidence rate. A control strategy with vaccine failure is applied to the switched epidemic model with generalized compartments. The control strategies are analyzed to determine their success in eradicating the disease. Some examples are given, with simulations, to illustrate the threshold conditions found.     

Two-storage inventory model with lot-size dependent fuzzy lead-time under possibility constraints via genetic algorithm

Multi-item inventory models with stock dependent demand and two storage facilities are developed in a fuzzy environment where processing time of each unit is fuzzy and the processing time of a lot is correlated with its size. These are order-quantity reorder-point models with back-ordering if required. Here possibility and crisp constraints on investment and capacity of the small storehouse respectively are considered. The models are formulated as fuzzy chance constrained programming problem and is solved via generalized reduced gradient (GRG) technique when crisp equivalent of the constraints are available. A genetic algorithm (GA) is developed based on fuzzy simulation and entropy where region of search space gradually decreases to a small neighborhood of the optima and it is used to solve the models whenever the equivalent crisp form of the constraint is not available. The models are illustrated with some numerical examples and some sensitivity analyses have been done. For some particular cases results observed via GRG and GA are compared.     

Wiener processes with random effects for degradation data

This article studies the maximum likelihood inference on a class of Wiener processes with random effects for degradation data. Degradation data are special case of functional data with monotone trend. The setting for degradation data is one on which n independent subjects, each with a Wiener process with random drift and diffusion parameters, are observed at possible different times. Unit-to-unit variability is incorporated into the model by these random effects. EM algorithm is used to obtain the maximum likelihood estimators of the unknown parameters. Asymptotic properties such as consistency and convergence rate are established. Bootstrap method is used for assessing the uncertainties of the estimators. Simulations are used to validate the method. The model is fitted to bridge beam data and corresponding goodness-of-fit tests are carried out. Failure time distributions in terms of degradation level passages are calculated and illustrated.     

单调空间与度量化定理

本文回答了关于MCM空间遗传性的一个问题,讨论了k-MCM空间是k半层空间的条件,得到了一些用g函数刻划的度量化定理.主要结论有:MCM空间是关于Fσ子空间遗传的;在正规空间类中,q空间(ωN空间,k-MCM空间)是关于开Fσ子空间遗传的;如果X是具有Gδ对角线的正则次中紧 k-MCM空间,则X是k半层空间;X是可度量化空间的充要条件是存在X上的g函数满足对X中任意不相交的闭集F与紧集C,都有某个n∈ω,使得(∪x∈F g(n,x))∩(∪y∈C g(n,y))=(?).     

Stable generalization of convex functions

A kind of generalized convex functions is said to be stable with respect to some property (P) if this property is maintaincd during an arbitrary function from this class is disturbed by a linear functional with sufficiently small norm. Unfortunately. known generallzed convexities iike quasicunvexity, explicit quasiconvexity. and pseudoconvexity are not stable with respect to such optimization properties which are expected to be true by these generalizations, even if the domain ol the functions is compact. Therefore, we introduce the notion of s-quasiconvex functions. These functions are quasiconvex, explicitly quasicon vex. and pseudoconvex if they are continuously differentiable. Especially, the s-quasiconvexity is stable with respect to the following important properties: (Pl) all lower level sets are convex, (P2) each local minimum is a global minimum. and (P3) each stationary point is a global minimizer. In this paper, different aspects. of s–quasiconvexity and its stability are investigated.     

Tracking nonperiodic trajectories with the overtaking criterion

A tracking problem on an infinite time interval is studied, where the plant is linear with quadratic cost, and the tracked trajectory is not necessarily period. Optimal solutions with respect to the overtaking criterion are studied. Existence and uniqueness of such optimal solutions are proved and they are shown to be given by a linear feedback law which is the same as in the periodic case. A close relation between the solutions of tracking problems where the tracked trajectories are different only for very large times is established.     

Two storage inventory model in a mixed environment

Multi-item inventory models with two storage facility and bulk release pattern are developed with linearly time dependent demand in a finite time horizon under crisp, stochastic and fuzzy-stochastic environments. Here different inventory parameters—holding costs, ordering costs, purchase costs, etc.—are assumed as probabilistic or fuzzy in nature. In particular cases stochastic and crisp models are derived. Models are formulated as profit maximization principle and three different approaches are proposed for solution. In the first approach, fuzzy extension principle is used to find membership function of the objective function and then it’s Graded Mean Integration Value (GMIV) for different optimistic levels are taken as equivalent stochastic objectives. Then the stochastic model is transformed to a constraint multi-objective programming problem using Stochastic Non-linear Programming (SNLP) technique. The multi-objective problems are transferred to single objective problems using Interactive Fuzzy Satisfising (IFS) technique. Finally, a Region Reducing Genetic Algorithm (RRGA) based on entropy has been developed and implemented to solve the single objective problems. In the second approach, the above GMIV (which is stochastic in nature) is optimized with some degree of probability and using SNLP technique model is transferred to an equivalent single objective crisp problem and solved using RRGA. In the third approach, objective function is optimized with some degree of possibility/necessity and following this approach model is transformed to an equivalent constrained stochastic programming problem. Then it is transformed to an equivalent single objective crisp problem using SNLP technique and solved via RRGA. The models are illustrated with some numerical examples and some sensitivity analyses have been presented.     

Approximations of optimal control problems for semilinear elliptic equations with discontinuous coefficients and solutions and with control in matching boundary conditions

Mathematical formulations of nonlinear optimal control problems for semilinear elliptic equations with discontinuous coefficients and solutions and with control in matching conditions are examined. Finite difference approximations of optimization problems are constructed, and the approximation error is estimated with respect to the state and the cost functional. Weak convergence of the approximations with respect to the control is proved. The approximations are regularized using Tikhonov regularization.     

The simulation of offshore turbulent dispersion using seeded eddies

A way of representing turbulence in a two-dimensional situation is introduced appropriate to depth-independent offshore fluid mechanics. The turbulence is simulated by a collection of eddies, each of which has an analytically simple form but whose size, strength and position is governed by stochastically assigned variables. The problem addressed here is how contaminant is dispersed in such an eddy field. A number of experiments are performed whereby the eddies are seeded with marked particles that move with the fluid. The variance of these particles is monitored as time varies, and the results are compared with an assumed power law distribution. Although not a perfect fit, the results are in general accord with a power law with index between 1.5 and 2.5, which is in agreement with the observed power law of 2.34 due to Okubo, and a marked improvement on random walk models which give a variance directly proportional to time. Some further applications of this technique are discussed, namely the simulation of turbulent boundary layers and the simulation of the cascade of energy up turbulent length scales.     

Splitting and conjoining objects in monotopological categories

Power-sets are defined for any concrete category (over Set) with finite concrete products, and their structure described for monotopological categories. These sets are used to define the notions of splitting object and of conjoining object. Characterizations of the existence of these objects in monotopological categories are given. It is proved that no proper monotopological category can be concretely cartesian closed. Most well-known monotopological categories with splitting objects are topological or are c-categories, but it is shown that there are many proper monotopological categories which are not c-categories, and yet have splitting objects, and may even be cartesian closed. One of the characterizations of the existence of splitting objects is used to prove that a monotopological category with splitting objects is cartesian closed iff the largest initial completion in which it is epireflective is cartesian closed iff its MacNeille completion is cartesian closed.     

A Dynamic Model for a Hybrid Articulated Robot

The dynamic modeling of hybrid systems, consisting of flexible and rigid parts results in large partial differential equation systems (PDE). With the assumption of small deflections and the Ritz expansion the PDE can be approximated by an ordinary differential equation system (ODE) but the number of degrees of freedom is generally high. In this paper a hybrid articulated robot with 2 flexible links and 6 joints is under consideration. The joints are equipped with Harmonic Drive gears with high gear ratio but relative low stiffness. Therefore additionally degrees of freedom are introduced for the elastic deflection of the gears. The links are modeled with flexibility in two bending directions and in torsional sense. To be able to achieve structured equations the projection equation in subsystem representation is used. The projection equation is based on the momentum and the angular momentum equations of each single finite or infinitesimal body which are projected into the space of minimal coordinates and subsequently are summed up. Groups of bodies are collected to the so called subsystems with separated describing velocities. These subsystems are linked together with the kinematical chain. Because the robot is tree structured it is possible to obtain an explicit expression for the second derivatives of the minimal coordinates with a recursive scheme (O(n) efficiency). The robot is controlled with a feed forward controller and a linear PD joint controller. Simulation results and measured data are presented and compared. (© 2008 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

An optimal asymptotic-numerical method for convection dominated systems having exponential boundary layers

The convection dominated diffusion problems are studied. Higher order accurate numerical methods are presented for problems in one and two dimensions. The underlying technique utilizes a superposition of given problem into two independent problems. The first one is the reduced problem that refers to the outer or smooth solution. Stretching transformation is used to obtain the second problem for inner layer solution. The method considered for outer or degenerate problems are based on higher order Runge–Kutta methods and upwind finite differences. However, inner problem is solved analytically or asymptotically. The schemes presented are proved to be consistent and stable. Possible extensions to delay differential equations and to nonlinear problems are outlined. Numerical results for several test examples are illustrated and a comparative analysis is presented. It is observed that the method presented is highly accurate and easy to implement. Moreover, the numerical results obtained are not only comparable with the exact solution but also in agreement with the theoretical estimates.     

Quadratic forms over formally real fields with eight square classes

Complete classification of formally real fields with 8 square classes with respect to the behaviour of quadratic forms is given. Two fields F and K are equivalent with respect to quadratic forms if the quadratic form schemes of the two fields are isomorphic or in other words, if the Witt rings W(F) and W(K) are isomorphic. It is shown here that for formally real fields with 8 square classes there are exactly seven possible quadratic form schemes and for each of the seven schemes a formally real field with 8 square classes and the given scheme is constructed.     

Micromechanical analysis of the stress transfer in single-fiber composite: The influence of the uniform and graded interphase with finite-thickness

A theoretical model is developed to analyze the stress transfer between fiber and matrix through the interphase with finite thickness. The Young's modulus of interphase is assumed to be homogeneous uniform or power-graded along radial direction while other material parameters are constants. The bonds between fiber and interphase as well as between interphase and matrix are perfect. The geometrical equations are strictly satisfied except that the radial displacement gradient with respect to the axial direction is neglected, as its magnitude is much smaller than that of the axial displacement gradient with respect to the radial direction. The equilibrium equations along radial direction are strictly satisfied, while the equilibrium equations along axial direction are satisfied in the integral forms. In addition, both the interfacial displacement and stress continuity conditions as well as stress boundary conditions are enforced exactly. Two coupled 2nd-order ordinary differential equations can be obtained in terms of average axial stresses in fiber and matrix. Finite element analysis (FEA) with refined mesh for single-fiber composite containing uniform interphase with finite thickness is developed to validate the present model. Series of parameter studies are performed to investigate the influence of interphase properties and thickness as well as the fiber volume content and model length on the stress distribution in composites.     

Thermodynamics of multi-temperature fluids with applications to turbulence modelling

The thermodynamics of a fluid with several distinct temperatures is considered. The balance laws are studied and a set of linear constitutive equations is derived; equations of motion are obtained and discussed. The application of the developed theory to turbulence modelling is studied and results are shown to be consistent with the one-equation models of turbulence. The thermodynamics of turbulence is also briefly discussed.     

Nonlinear dynamics of high-dimensional models of in-plane and out-of-plane vibration in an axially moving viscoelastic beam

Nonlinear dynamics of high-dimensional models of an axially moving viscoelastic beam with in-plane and out-of-plane vibration with combined parametric and forcing excitations are investigated by the incremental harmonic balance (IHB) method in this paper. Governing equations of transverse in-plane and out-of-plane and longitudinal vibration are obtained basing on the Hamilton's principle. The Galerkin method is used to separate time variable and spatial variable to obtain a set of multi-order differential equations. The IHB method with the fast Fourier transform (FFT) is used to solve periodic response of high-dimensional models of the beam for which convergent mode is reached. Stability of the steady-state periodic solutions is analyzed using the multivariable Floquet theory. Particular attention is paid to in-plane and out-of-plane vibration on convergent mode of the beam with combined parametric and forcing excitations. Multiple solutions are observed, and jump phenomena between in-plane and out-of-plane vibration with different transverse cross sections are discovered.     

基于改进区间数密度集结算子指标群赋权方法

针对指标数据已知,而权重数据未知的群组赋权问题,给出了一种基于改进的区间数密度集结算子来进行指标群组赋权的决策方法。首先给出了区间数和区间数密度集结算子(IDM)的定义及性质,改进了以前区间数聚类的方法,应用直接法对一维区间数据组进行聚类,并定义了模糊统计量,以确定最为合理的一种聚类方式。然后基于改进的区间数密度集结算子这种数学模型,来解决指标值数据已知,而权重未知的群组赋权问题。最后举例说明该方法的可行性和实用性。     

Control Strategies for a Hybrid Articulated Robot

Due to higher requirements in productivity and cost efficiency of production lines, robots and other manipulators have to move faster. One possibility to fulfill the mentioned goals is to build lightweight constructions having elastic deformations in joints and links. The elastic components tend to vibrations and static deflections. Methods that compensate or minimize these drawbacks are the focus of this paper. An articulated robot with 6 joints and flexibility in joints and links is under consideration. The joints are actuated by DC motors combined with Harmonic Drive gears which offer high gear ratios but undergo elastic deformations. The links are flexible in two bending directions and in torsional sense. To achieve ordinary differential equations, a Ritz approach together with the projection equation is used. The obtained model is used for feedforward and feedback control design. Based on reference trajectories and on a rigid body model, estimations for the elastic deflections are calculated. These deflections are used to alter the reference trajectory in order to minimize the error of the tool center point. For basic active damping, non-local curvature feedback is used. Together with PD joint control and the feedforward control, satisfying results are obtained. Additionally, a sliding control approach is presented. The stiffness of the tool center point is enhanced with the drawback of less active damping. Simulation results and measured data are presented and compared. (© 2010 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)