A nonlinear truck model and its treatment as a multibody system

A planar vertical truck model with nonlinear suspension and its multibody system formulation are presented. The equations of motion of the model form a system of differential-algebraic equations (DAEs). All equations are given explicitly, including a complete set of parameter values, consistent initial values, and a sample road excitation. Thus the truck model allows various investigations of the specific DAE effects and represents a test problem for algorithms in control theory, mechanics of multibody systems, and numerical analysis. Several numerical tests show the properties of the model.     

A compact formulation of an elastoplastic analysis problem

Two important problems in the area of engineering plasticity are limit load analysis and elastoplastic analysis. It is well known that these two problems can be formulated as linear and quadratic programming problems, respectively (Refs. 1–2). In applications, the number of variables in each of these mathematical programming problems tends to be large. Consequently, it is important to have efficient numerical methods for their solution. The purpose of this paper is to present a method which allows the quadratic programming formulation of the elastoplastic analysis to be reformulated as an equivalent quadratic programming problem which has significantly fewer variables than the original formulation. Indeed, in Section 4, we will present details of an example for which the original quadratic programming formulation required 297 variables and for which the equivalent formulation presented here required only two variables. The method is based on a characterization of the entire family of optimal solutions for a linear programming problem.This research was supported by the Natural Science and Engineering Council of Canada under Grant No. A8189 and by a Leave Fellowship from the Social Sciences and Humanities Research Council of Canada. The author takes pleasure in acknowledging many stimulating discussions with Professor D. E. Grierson.     

Smooth boundary integration method in a linear two-dimensional problem of incompressible fluid and solid

This paper presents a new application of a theoretical and computational method of smooth boundary integration which belongs to the methods of boundary integral equations. Smooth integration is not a method of approximation. In its final analytical form, a smooth-kernel integral equation is computerized easily and accurately.

Smooth integration is associated with a “pressure-vorticity” formulation which covers linear problems in elasticity and fluid mechanics. The solution presented herein is essentially the same as that reported in an earlier paper for regular elasticity. The constraint of incompressibility does not cause difficulties in the pressure-vorticity formulation.

The linear fluid mechanics problem formulated and solved in this paper covers Stokes' problem of a slow viscous flow, and has a wider interpretation. The translational inertia forces are incorporated in the linear problem, as in Euler's dynamic theory of inviscid flow. The centrifugal inertia forces are left for the non-linear problem. The linear problem is perceived as a step in solution of the non-linear problems.     

An Outer Approximations Approach to Reliability-Based Optimal Design of Structures

We present a new formulation of the problem of minimizing the initial cost of a structure subject to a minimum reliability requirement, expressed in terms of the so-called design points of the first-order reliability theory, i.e., points on limit-state surfaces that are nearest to the origin in a transformed standard normal space, as well as other deterministic constraints. Our formulation makes it possible to use outer approximations algorithms for the solution of such optimal design problems, eliminating some of the major objections associated with treating them as bilevel optimization problems. A numerical example is presented that illustrates the reliability and efficiency of the algorithm.     

Plastic flow in a conical channel: Qualitative features of the solutions under different yield conditions

The problem of the convergence of the solutions of problems of plasticity theory, with a yield condition which depends on the hydrostatic stress, to solutions based on classical plasticity theory with von Mises or Tresea conditions is considered, with a particular choice of the parameters of the material model. For the case of axisymmetric flow of material in a channel with converging and diverging walls, solutions according to two plasticity theories with a yield condition which depends on the hydrostatic stress are compared with the classical solution. It is shown that only the solution using Spencer's model shows all the main features of the classical solution. As the internal criterion of the choice of the preferred plasticity theory when examining a special class of problems, it is suggested that the criterion of the convergence of the solutions to the solutions of classical plasticity theory should be used.     

Two formulations of elastoplastic problems and the theoretical determination of the location of neck formation in samples under tension

Two formulations of elastoplastic problems in the mechanics of deformable solids with finite displacements and deformations are investigated. The first of these is formulated starting from the classical geometrically non-linear equations of the theory of elasticity and plasticity, in which the components of the Cauchy–Green strain tensor, associated with the components of the conditional stress tensor by physically non-linear relations according to flow theory in the simplest version of their representation, are taken as a measure of the deformations. The second formulation is based on the introduction of the true tensile and shear strains which, according to Novoshilov, are associated with the components of the true stresses by physical relations of the above-mentioned form. It is shown that, in the second version of the formulation of the problem, the use of the corresponding equations, complied taking account of the elastoplastic properties of the material with correct modelling of the ends of cylindrical samples and the method of loading (stretching) them, enables the location of the formation of a neck to be determined theoretically and enables the initial stage of its formation to be described without making any assumptions regarding the existence of initial irregularities in the geometry of the samples.     

A new numerical approach for the advective-diffusive transport equation

A new numerical solution procedure is presented for the one-dimensional, transient advective-diffusive transport equation. The new method applies Herrera's algebraic theory of numerical methods to the spatial derivatives to produce a semi-discrete approximation. The semi-discrete system is then solved by standard time marching algorithms. The algebraic theory, which involves careful choice of test functions in a weak form statement of the problem, leads to a numerical approximation that inherently accommodates different degrees of advection domination. Algorithms are presented that provide either nodal values of the unknown function or nodal values of both the function and its spatial derivative. Numerical solution of several test problems demonstrates the behavior of the method.     

Economical and ecological utility oriented analysis of the process of anaerobic digestion of waste waters

A problem which has been constantly emphasized is the creation of criteria adequate to characterize the complexity of ecological analysis. The objective of the present paper is to demonstrate the capabilities of multiattribute utility theory in difficult-to-formalize problems. The multiattribute utility and the proposed algorithms provide a logically and operationally tested method which includes value in complex ecological problems. The results obtained and the constructed utility functions should be accepted as an iterative stage in real investigations, rather than as complete research that offer a final decision. The value estimations of the decision maker are the basis for interest in a given ecological problem. But they are often not explicitly or consistently addressed in the real investigations. The proposed methods account for otherwise uninterpretable information. The constructed value function can be used for automatic computer control and monitoring of anaerobic waste water digestion, which could reveal a new potential from the practical point of view.     

A New Formulation for Spanning Trees

Spanning trees are fundamental structures in graph theory. Furthermore, computing them is a central part in many relevant algorithms, used in either practical or theoretical applications. The classical Minimum Spanning Tree problem is solvable in polynomial time but almost all of its variants are NP-Hard. In this paper, a novel polynomial size mixed integer linear programming formulation is introduced for spanning trees. This formulation is based on a new characterization we propose for acyclic graphs. Preliminary computational results show that this formulation is capable of solving small instances of the diameter constrained minimum spanning tree problem. It should be possible to strengthen the formulation to tackle larger instances of that problem. Additionally, our spanning tree formulation may prove to be a more effective model for some related applications.     

Fractional multi-commodity flow problem: Duality and optimality conditions

This paper deals with multi-commodity flow problem with fractional objective function. The optimality conditions and the duality concepts of this problem are given. For this aim, the fractional linear programming formulation of this problem is considered and the weak duality, the strong direct duality and the weak complementary slackness theorems are proved applying the traditional duality theory of linear programming problems which is different from same results in Chadha and Chadha (2007) [1]. In addition, a strong (strict) complementary slackness theorem is derived which is firstly presented based on the best of our knowledge. These theorems are transformed in order to find the new reduced costs for fractional multi-commodity flow problem. These parameters can be used to construct some algorithms for considered multi-commodity flow problem in a direct manner. Throughout the paper, the boundedness of the primal feasible set is reduced to a weaker assumption about solvability of primal problem which is another contribution of this paper. Finally, a real world application of the fractional multi-commodity flow problem is presented.     

Approximations of thermoelastic and viscoelastic control systems

This paper deals with the development and analysis of well-posed models and computational algorithms for control of a class of partial differential equations that descrive the motions of thermo-viscoelastic structures. We first present an abstract “state space” framework and general well-posedness result that can be applied to a large class of thermo-elastic and thermo-viscoelastic models. This state space framework is used in the development of a computational scheme to be used in the solution of an LQR control problem. A detailed convergence proof is provided for the viscoelastic model, and several numerical results are presented to illustrate the theory and to analyze problems for which the theory is incomplete.     

Variational formulation of a modified couple stress theory and its application to a simple shear problem

A variational formulation is provided for the modified couple stress theory of Yang et al. by using the principle of minimum total potential energy. This leads to the simultaneous determination of the equilibrium equations and the boundary conditions, thereby complementing the original work of Yang et al. where the boundary conditions were not derived. Also, the displacement form of the modified couple stress theory, which is desired for solving many problems, is obtained to supplement the existing stress-based formulation. All equations/expressions are presented in tensorial forms that are coordinate-invariant. As a direct application of the newly obtained displacement form of the theory, a simple shear problem is analytically solved. The solution contains a material length scale parameter and can capture the boundary layer effect, which differs from that based on classical elasticity. The numerical results reveal that the length scale parameter (related to material microstructures) can have a significant effect on material responses.      

Variational formulation of a modified couple stress theory and its application to a simple shear problem

A variational formulation is provided for the modified couple stress theory of Yang et al. by using the principle of minimum total potential energy. This leads to the simultaneous determination of the equilibrium equations and the boundary conditions, thereby complementing the original work of Yang et al. where the boundary conditions were not derived. Also, the displacement form of the modified couple stress theory, which is desired for solving many problems, is obtained to supplement the existing stress-based formulation. All equations/expressions are presented in tensorial forms that are coordinate-invariant. As a direct application of the newly obtained displacement form of the theory, a simple shear problem is analytically solved. The solution contains a material length scale parameter and can capture the boundary layer effect, which differs from that based on classical elasticity. The numerical results reveal that the length scale parameter (related to material microstructures) can have a significant effect on material responses.     

Berth allocation at indented berths for mega-containerships

This paper addresses the berth allocation problem at a multi-user container terminal with indented berths for fast handling of mega-containerships. In a previous research conducted by the authors, the berth allocation problem at a conventional form of the multi-user terminal was formulated as a nonlinear mathematical programming, where more than one ship are allowed to be moored at a specific berth if the berth and ship lengths restriction is satisfied. In this paper, we first construct a new integer linear programming formulation for easier calculation and then the formulation is extended to model the berth allocation problem at a terminal with indented berths, where both mega-containerships and feeder ships are to be served for higher berth productivity. The berth allocation problem at the indented berths is solved by genetic algorithms. A wide variety of numerical experiments were conducted and interesting findings were explored.     

Strain gradient plasticity solution for an internally pressurized thick-walled cylinder of an elastic linear-hardening material

An elastic-plastic solution is presented for an internally pressurized thick-walled plane strain cylinder of an elastic linear-hardening plastic material. The solution is derived in a closed form using a strain gradient plasticity theory. The inner radius of the cylinder enters the solution not only in non-dimensional forms but also with its own dimensional identity, which differs from that in classical plasticity based solutions and makes it possible to capture the size effect at the micron scale. The classical plasticity solution of the same problem is recovered as a special case of the current solution. To further illustrate the newly derived solution, formulas and numerical results for the plastic limit pressure are provided. These results reveal that the load-carrying capacity of the cylinder increases with decreasing inner radius at the micron scale. It is also seen that the macroscopic behavior of the pressurized cylinder can be well described by using classical plasticity based solutions.     

A new adaptive approach of the Metropolis-Hastings algorithm applied to structural damage identification using time domain data

In the present work, the formulation and solution of the inverse problem of structural damage identification is presented based on the Bayesian inference, a powerful approach that has been widely used for the formulation of inverse problems in a statistical framework. The structural damage is continuously described by a cohesion field, which is spatially discretized by the finite element method, and the solution of the inverse problem of damage identification, from the Bayesian point of view, is the posterior probability densities of the nodal cohesion parameters. In this approach, prior information about the parameters of interest and the quantification of the uncertainties related to the magnitudes measured can be used to estimate the sought parameters. Markov Chain Monte Carlo (MCMC) method, implemented via the Metropolis-Hastings (MH) algorithm, is commonly used to sample such densities. However, the conventional MH algorithm may present some difficulties, for instance, in high dimensional problems or when the parameters of interest are highly correlated or the posterior probability density is very peaked. In order to overcome these difficulties, a new adaptive MH algorithm (P-AMH) is proposed in the present work. Numerical results related to an inverse problem of damage identification in a simply supported Euler-Bernoulli beam are presented. Synthetic experimental time domain data, obtained with different damage scenarios, and noise levels, were addressed with the aim at assessing the proposed damage identification approach. An adaptive MH algorithm (H-AMH) and the conventional MH algorithm, already consolidated in the literature, were also considered for comparison purposes. The numerical results show that both adaptive algorithms outperformed the conventional MH. Besides, the P-AMH provided Markov chains with faster convergence and better mixing than the ones provided by the H-AMH.     

Strain gradient plasticity solution for an internally pressurized thick-walled cylinder of an elastic linear-hardening material

An elastic-plastic solution is presented for an internally pressurized thick-walled plane strain cylinder of an elastic linear-hardening plastic material. The solution is derived in a closed form using a strain gradient plasticity theory. The inner radius of the cylinder enters the solution not only in non-dimensional forms but also with its own dimensional identity, which differs from that in classical plasticity based solutions and makes it possible to capture the size effect at the micron scale. The classical plasticity solution of the same problem is recovered as a special case of the current solution. To further illustrate the newly derived solution, formulas and numerical results for the plastic limit pressure are provided. These results reveal that the load-carrying capacity of the cylinder increases with decreasing inner radius at the micron scale. It is also seen that the macroscopic behavior of the pressurized cylinder can be well described by using classical plasticity based solutions.     

Lower Bounds for Fixed Spectrum Frequency Assignment

Determining lower bounds for the sum of weighted constraint violations in fixed spectrum frequency assignment problems is important in order to evaluate the performance of heuristic algorithms. It is well known that, when adopting a binary constraints model, clique and near-clique subproblems have a dominant role in the theory of lower bounds for minimum span problems. In this paper we highlight their importance for fixed spectrum problems. We present a method based on the linear relaxation of an integer programming formulation of the problem, reinforced with constraints derived from clique-like subproblems. The results obtained are encouraging both in terms of quality and in terms of computation time.     

High resolution multi-moment finite volume method for supersonic combustion on unstructured grids

In this study, we present a novel numerical model for simulating detonation waves on unstructured grids. In contrast to the conventional finite volume method (FVM), two types of moment comprising the volume-integrated average (VIA) and the point value (PV) at the cell vertex are treated as the evolution variables for the reacting Euler equations. The VIA is computed based on a finite volume formulation of the flux form where the conventional Riemann problem is solved by the HLLC Riemann solver. The PV is updated in a point-wise manner by using the differential formulation where the Roe solver is used to compute the differential Riemann problems. In order to increase the accuracy around discontinuities, numerical oscillations and dissipations are reduced using the boundary variation diminishing algorithm. Convergence tests demonstrated that the proposed model could achieve third-order accuracy with unstructured grids for reacting Euler equations. The high resolution property of the proposed method was verified based on simulations of several detonation wave propagation problems in two and three dimensions. In particular, the current model could resolve the cellular structures with fewer degrees of freedom for the unstable oblique detonation wave problem. These fine structures may be smoothed out by the conventional FVM due to the excessive amount of numerical dissipation errors. Importantly, a simulation of stiff detonation waves showed that the proposed method could capture the correct position of the reaction front whereas the conventional FVMs produced spurious phenomena. Thus, the proposed model can obtain highly accurate solutions for detonation problems on unstructured grids, which is highly advantageous for real applications involving complex geometrical configurations.     

The matrix sign function and computations in systems

The matrix sign function has several interesting properties which form the basis of new solution algorithms for problems which occur frequently in systems and control theory applications. Presented in this paper are new algorithms, based on the matrix sign function, for the solution of algebraic matrix Riccati equations, Lyapunov equations, coupled Riccati equations, spectral factorization, matrix square roots, pole assignment, and the algebraic eigenvalue-eigenvector problem. Examples of the application of each algorithm are also presented.