A Lagrange–Eulerian formulation of an axially moving beam based on the absolute nodal coordinate formulation

In the scope of this paper, a finite-element formulation for an axially moving beam is presented. The beam element is based on the absolute nodal coordinate formulation, where position and slope vectors are used as degrees of freedom instead of rotational parameters. The equations of motion for an axially moving beam are derived from generalized Lagrange equations in a Lagrange–Eulerian sense. This procedure yields equations which can be implemented as a straightforward augmentation to the standard equations of motion for a Bernoulli–Euler beam. Moreover, a contact model for frictional contact between an axially moving strip and rotating rolls is presented. To show the efficiency of the method, simulations of a belt drive are presented.     

A comparison of finite elements for nonlinear beams: the absolute nodal coordinate and geometrically exact formulations

Two of the most popular finite element formulations for solving nonlinear beams are the absolute nodal coordinate and the geometrically exact approaches. Both can be applied to problems with very large deformations and strains, but they differ substantially at the continuous and the discrete levels. In addition, implementation and run-time computational costs also vary significantly. In the current work, we summarize the main features of the two formulations, highlighting their differences and similarities, and perform numerical benchmarks to assess their accuracy and robustness. The article concludes with recommendations for the choice of one formulation over the other.     

Clamped end conditions and cross section deformation in the finite element absolute nodal coordinate formulation

In the finite element absolute nodal coordinate formulation (ANCF), the elimination of the relative translations and rotations at a point does not necessarily define a fully clamped joint, particularly in the case of fully parameterized ANCF finite elements that allow for the deformation of the cross section. In this investigation, the formulations and results of two different sets of clamped end conditions that define two different joints are compared. The first joint, called the partially clamped joint, eliminates only the translations and rotations at a point on the cross section. The second joint, called the fully clamped joint, eliminates all the translation, rotation and deformation degrees of freedom at a point on the cross section. The kinematic equations that define the partially and fully clamped joints are developed, and the dynamic equations used in the comparative numerical study presented in this paper are shown. As discussed in this investigation, the fully clamped joint does not allow for the deformation of the cross section at the joint node since the gradient vectors remain orthogonal unit vectors. The partially clamped joint, on the other hand, allows for the deformation of the cross section. Nanson’s formula is used as a measure of the deformation of the cross section in the case of the partially clamped joint. A very flexible pendulum that has a rigid body attached to its free end is used to compare the results of the partially and fully clamped joints. The numerical results obtained using this very flexible pendulum example show that, while the type of joint (partially or fully clamped) does not significantly affect the gross reference motion of the system, there are fundamental differences between the two joints since the partially clamped joint allows for the cross section deformations at the joint node.     

Simulation of a viscoelastic flexible multibody system using absolute nodal coordinate and fractional derivative methods

This paper employs a new finite element formulation for dynamics analysis of a viscoelastic flexible multibody system. The viscoelastic constitutive equation used to describe the behavior of the system is a three-parameter fractional derivative model. Based on continuum mechanics, the three-parameter fractional derivative model is modified and the proposed new fractional derivative model can reduce to the widely used elastic constitutive model, which meets the continuum mechanics law strictly for pure elastic materials. The system equations of motion are derived based on the absolute nodal coordinate formulation (ANCF) and the principle of virtual work, which can relax the small deformation assumption in the traditional finite element implementation. In order to implement the viscoelastic model into the absolute nodal coordinate, the Grünwald definition of the fractional derivative is employed. Based on a comparison of the HHT-I3 method and the Newmark method, the HHT-I3 method is used to solve the equations of motion. Another particularity of the proposed method based on the ANCF method lies in the storage of displacement history only during the integration process, reducing the numerical computation considerably. Numerical examples are presented in order to analyze the effects of the truncation number of the Grünwald series (fading memory phenomena) and the value of several fractional model parameters and solution convergence aspects. An erratum to this article can be found at     

Three-dimensional beam propagation method based on the variable transformed Galerkin''''s method

A novel three-dimensional beam propagation method (BPM) based on the variable transformed Galerkin's method is introduced for simulating optical field propagation in three-dimensional dielectric structures. The infinite Cartesian x-y plane is mapped into a unit square by a tangent-type function transformation. Consequently, the infinite region problem is converted into the finite region problem. Thus, the boundary truncation is eliminated and the calculation accuracy is promoted. The three-dimensional BPM basic equation is reduced to a set of first-order ordinary differential equations through sinusoidal basis function, which fits arbitrary cladding optical waveguide, then direct solution of the resulting equations by means of the Runge-Kutta method. In addition, the calculation is efficient due to the small matrix derived from the present technique. Both z-invariant and z-variant examples are considered to test both the accuracy and utility of this approach.     

An analysis of the time integration algorithms for the finite element solutions of incompressible Navier–Stokes equations based on a stabilised formulation

This work is concerned with the analysis of time integration procedures for the stabilised finite element formulation of unsteady incompressible fluid flows governed by the Navier–Stokes equations. The stabilisation technique is combined with several different implicit time integration procedures including both finite difference and finite element schemes. Particular attention is given to the generalised-α method and the linear discontinuous in time finite element scheme. The time integration schemes are first applied to two model problems, represented by a first order differential equation in time and the one dimensional advection–diffusion equation, and subjected to a detailed mathematical analysis based on the Fourier series expansion. In order to establish the accuracy and efficiency of the time integration schemes for the Navier–Stokes equations, a detailed computational study is performed of two standard numerical examples: unsteady flow around a cylinder and flow across a backward facing step. It is concluded that the semi-discrete generalised-α method provides a viable alternative to the more sophisticated and expensive space–time methods for simulations of unsteady flows of incompressible fluids governed by the Navier–Stokes equations.     

Recursive formulations for multibody systems with frictional joints based on the interaction between bodies

In practice, the clearances of joints in a great number of mechanical systems are well under control. In these cases, some of the existing methods become unpractical because of the little differences in the order of magnitude between relative movements and computational errors. Assuming that the effects of impacts are negligible, we proved that both locations and forces of contacts in joints can be fully determined by parts of joint reaction forces. Based on this fact, a method particularly suited for multibody systems possessing frictional joints with tiny clearances is presented. In order to improve the efficiency of computation, recursive formulations are proposed based on the interactions between bodies. The proposed recursive formulations can improve the computation of joint reaction forces. With the methodology presented in this paper, not only the motion of bodies in a multibody system but also the details about the contacts in joints, such as forces of contacts and locations of contact points, can be obtained. Even with the assumption of impact free, the instants of possible impacts can be detected without relying upon any ambiguous parameters, as indicated by numerical examples in this paper.     

Extracting exact answers from large-scale corpus based on hybrid strategy

This paper provides a novel and efficient method for extracting exact textual answers from the returned documents that are retrieved by traditional IR system in large-scale collection of texts. The main intended contribution of this paper is to propose System Similarity Model （SSM）, which can be considered as an extension of vector space model （VSM） to rank passages, It presents a method of formalized answer extraction based on pattern learning and applies binary logistic regression model （LRM）, which seldom be used in IE to extract special information from candidate data sets. The parameters estimated for the data gathers with serious problem of data sparse, therefore we take stratified sampling method, and improve traditional logistic regression model parameters estimated methods. The series of experimental results show that the overall performance of our system is good and our approach is effective. Our system, lnsun05QAl, which participated in QA track of TREC 2005 obtained excellent results.     

A Prolate-Element Method for Nonlinear PDEs on the Sphere

A p-type spectral-element method using prolate spheroidal wave functions (PSWFs) as basis functions, termed as the prolate-element method, is developed for solving partial differential equations (PDEs) on the sphere. The gridding on the sphere is based on a projection of the prolate-Gauss-Lobatto points by using the cube-sphere transform, which is free of singularity and leads to quasi-uniform grids. Various numerical results demonstrate that the proposed prolate-element method enjoys some remarkable advantages over the polynomial-based element method: (i) it can significantly relax the time step size constraint of an explicit time-marching scheme, and (ii) it can increase the accuracy and enhance the resolution.     

A C° continuous linear beam/bilinear plate flexure element

Design of a simple C°-continuous beam/plate flexure element based on a shear-deformable theory is attempted. Emphasis is placed on development of a low order linear/bilinear element. However, past experience shows that such elements become very “stiff” especially when thickness is reduced (conforming to Kirchhoff mode). This phenomenon is called “locking.” Attempts have been made in the last several years by a few investigators to overcome this problem.The locking problem is resolved here through a new approach. The total strain energy is split into bending and shear energies and an antiparameter in the shear energy term is introduced to avoid locking. By numerical experimentation on beam and plate problems it is shown that the present approach gives good results in the thin limit. It is also shown that the additional/spurious zero energy modes do not arise here because reduced/selective integration is avoided.     

The χ-schemes: a new consistent high-resolution formulation based on the normalized variable methodology

This paper deals with the formulation and testing of a new class of consistent high-resolution schemes, denoted as the χ-schemes. These schemes, combine consistency, accuracy and boundedness across systems of equations and are suitable for use in the simulation of multi-phase and multi-component flows. The consistency feature refers to the capability of these schemes to implicitly satisfy the additional algebraic constraint representing a global conservation relation governing certain sets of equations (e.g., species mass fraction, volume fraction, etc.). Four χ-schemes are implemented within an unstructured grid finite-volume framework, tested by solving four multi-component pure-advection test problems, and shown to be consistent.     

Biologically-inspired robot spatial cognition based on rat neurophysiological studies

This paper presents a robot architecture with spatial cognition and navigation capabilities that captures some properties of the rat brain structures involved in learning and memory. This architecture relies on the integration of kinesthetic and visual information derived from artificial landmarks, as well as on Hebbian learning, to build a holistic topological-metric spatial representation during exploration, and employs reinforcement learning by means of an Actor-Critic architecture to enable learning and unlearning of goal locations. From a robotics perspective, this work can be placed in the gap between mapping and map exploitation currently existent in the SLAM literature. The exploitation of the cognitive map allows the robot to recognize places already visited and to find a target from any given departure location, thus enabling goal-directed navigation. From a biological perspective, this study aims at initiating a contribution to experimental neuroscience by providing the system as a tool to test with robots hypotheses concerned with the underlying mechanisms of rats’ spatial cognition. Results from different experiments with a mobile AIBO robot inspired on classical spatial tasks with rats are described, and a comparative analysis is provided in reference to the reversal task devised by O’Keefe in 1983.     

AILU preconditioning for the finite element formulation of the incompressible Navier–Stokes equations

In this paper, the effect of a variable reordering method on the performance of “adapted incomplete LU (AILU)” preconditioners applied to the P2P1 mixed finite element discretization of the three-dimensional unsteady incompressible Navier–Stokes equations has been studied through numerical experiments, where eigenvalue distribution and convergence histories are examined. It has been revealed that the performance of an AILU preconditioner is improved by adopting a variable reordering method which minimizes the bandwidth of a globally assembled saddle-point type matrix. Furthermore, variants of the existing AILU(1) preconditioner have been suggested and tested for some three-dimensional flow problems. It is observed that the AILU(2) outperforms the existing AILU(1) with a little extra computing time and memory.     

Geometric nonlinear effects on the planar dynamics of a pivoted flexible beam encountering a point-surface impact

Flexible-body modeling with geometric nonlinearities remains a hot topic of research by applications in multibody system dynamics undergoing large overall motions. However, the geometric nonlinear effects on the impact dynamics of flexible multibody systems have attracted significantly less attention. In this paper, a point-surface impact problem between a rigid ball and a pivoted flexible beam is investigated. The Hertzian contact law is used to describe the impact process, and the dynamic equations are formulated in the floating frame of reference using the assumed mode method. The two important geometric nonlinear effects of the flexible beam are taken into account, i.e., the longitudinal foreshortening effect due to the transverse deformation, and the stress stiffness effect due to the axial force. The simulation results show that good consistency can be obtained with the nonlinear finite element program ABAQUS/Explicit if proper geometric nonlinearities are included in the floating frame formulation. Specifically, only the foreshortening effect should be considered in a pure transverse impact for efficiency, while the stress stiffness effect should be further considered in an oblique case with much more computational effort. It also implies that the geometric nonlinear effects should be considered properly in the impact dynamic analysis of more general flexible multibody systems.     

A mixed formulation for the direct approximation of the control of minimal $$L^2$$-norm for linear type wave equations

This paper deals with the numerical computation of null controls for the wave equation with a potential. The goal is to compute approximations of controls that drive the solution from a prescribed initial state to zero at a large enough controllability time. In [Cîndea, Fernández-Cara & Münch, Numerical controllability of the wave equation through primal methods and Carleman estimates, 2013], a so called primal method is described leading to a strongly convergent approximation of boundary controls: the controls minimize quadratic weighted functionals involving both the control and the state and are obtained by solving the corresponding optimality condition. In this work, we adapt the method to approximate the control of minimal square-integrable norm. The optimality conditions of the problem are reformulated as a mixed formulation involving both the state and its adjoint. We prove the well-posedeness of the mixed formulation (in particular the inf-sup condition) then discuss several numerical experiments. The approach covers both the boundary and the inner controllability. For simplicity, we present the approach in the one dimensional case.     

Unsupervised classification for polarimetric SAR images based on the improved CFSFDP algorithm

In polarimetric synthetic aperture radar (PolSAR) image processing, the number of classes is an important factor for PolSAR image classification. Therefore, how to accurately estimate the number of PolSAR image classes is an important issue. In this article, we propose a novel unsupervised classification method which can accurately estimate the number of classes for PolSAR images. First, the PolSAR image is initialized into many small clusters by using the complementary information from Yamaguchi decomposition and distribution characteristics of data. Second, the improved clustering by fast search and find of density peaks, named as improved CFSFDP algorithm, is introduced to select the appropriate category number. Finally, to improve the representation of each category, the PolSAR data set is classified by an iterative fine-tuning process based on a complex K-Wishart function. The performance of the proposed classification approach is presented and analysed on three real data sets. The experimental results show that the proposed classification method can accurately estimate the category number and enhance the classification accuracy in comparison with other traditional methods. It is also shown that the data distribution characteristic has the additional information beyond the target scattering decomposition, and this information is important for the initialization.     

A new model for the integrated vehicle-crew-rostering problem and a computational study on rosters

Operational planning within public transit companies has been extensively tackled but still remains a challenging area for operations research models and techniques. This phase of the planning process comprises vehicle-scheduling, crew-scheduling and rostering problems. In this paper, a new integer mathematical formulation to describe the integrated vehicle-crew-rostering problem is presented. The method proposed to obtain feasible solutions for this binary non-linear multi-objective optimization problem is a sequential algorithm considered within a preemptive goal programming framework that gives a higher priority to the integrated vehicle-crew-scheduling goal and a lower priority to the driver rostering goals. A heuristic approach is developed where the decision maker can choose from different vehicle-crew schedules and rosters, while respecting as much as possible management’s interests and drivers’ preferences. An application to real data of a Portuguese bus company shows the influence of vehicle-crew-scheduling optimization on rostering solutions.     

OpenMP based parallel normalized direct methods for sparse finite element linear systems

A new parallel normalized exact inverse algorithm is presented for solving sparse symmetric finite element linear systems on symmetric multiprocessor systems (SMP), based upon an antidiagonal motion approach (“wave”-like pattern) for overcoming the data dependencies. The proposed algorithm was implemented using OpenMP directives. Numerical results, such as speedups and efficiency, are presented illustrating the efficient performance on a symmetric multiprocessor computer system, where the proposed algorithmic solution method achieves good speedups.

The space of EDF deadlines: the exact region and a convex approximation

It is well known that the performance of computer controlled systems is heavily affected by delays and jitter occurring in the control loops, which are mainly caused by the interference introduced by other concurrent activities. A common approach adopted to reduce delay and jitter in periodic task systems is to decrease relative deadlines as much as possible, but without jeopardizing the schedulability of the task set. In this paper, we formally characterize the region of admissible deadlines so that the system designer can appropriately select the desired values to maximize a given performance index defined over the task set. Finally we also provide a sufficient region of feasible deadlines which is proved to be convex.