Optimal Design of Plane Frames under Uncertainty: A Comparison of Two Approaches

Using the first collapse theorem, the necessary and sufficient survival conditions of an elasto-plastic structure consist of the yield condition and the equilibrium condition. In practical applications several random model parameters have to be taken into account. This leads to a stochastic optimization problem which cannot be solved using the traditional methods. Instead of that, appropriate (deterministic) substitute problems must be formulated. Here, the design of plane frames is considered, where the applied load is supposed to be stochastic. In the first approach the recourse problem will be formulated in the standard form of stochastic linear programming (SLP). In order to apply efficient numerical solution procedures (LP-solvers), approximate recourse problems based on discretization (RPD) and the expected value problem (EVP) are introduced. In the second approach – based on the yield condition – a quadratic cost function will be introduced. After the formulation of the stochastic optimization problem, the expected cost based optimization problem (ECBOP) and the minimum expected cost problem (MECP) are formulated as representatives of appropriate substitute problems. Subsequently, comparative numerical results using these methods are presented. (© 2009 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Plastic Structural Analysis Under Stochastic Uncertainty

Problems from limit load or shakedown analysis are based on the convex, linear or linearized yield/strength condition and the linear equilibrium equation for the generic stress vector. Having to take into account, in practice, stochastic variations of the model parameters (e.g., yield stresses, plastic capacities) and external loadings, the basic stochastic plastic analysis problem must be replaced by an appropriate deterministic substitute problem. Instead of calculating approximatively the probability of failure based on a certain choice of failure modes, here, a direct approach is presented based on the costs for missing carrying capacity and the failure costs (e.g., costs for damage, repair, compensation for weakness within the structure, etc.). Based on the basic mechanical survival conditions, the failure costs may be represented by the minimum value of a convex and often linear program. Several mathematical properties of this program are shown. Minimizing then the total expected costs subject to the remaining (simple) deterministic constraints, a stochastic optimization problem is obtained which may be represented by a “Stochastic Convex Program (SCP) with recourse”. Working with linearized yield/strength conditions, a “Stochastic Linear Program (SLP) with complete fixed recourse” is obtained. In case of a discretely distributed probability distribution or after the discretization of a more general probability distribution of the random structural parameters and loadings as well as certain random cost factors one has a linear program (LP) with a so-called “dual decomposition data” structure. For stochastic programs of this type many theoretical results and efficient numerical solution procedures (LP-solver) are available. The mathematical properties of theses substitute problems are considered. Furthermore approximate analytical formulas for the limit load factor are given.     

Topology optimization of plane frames in case of random external loadings

Using the first collapse–theorem, the necessary and sufficient survival conditions of an elasto–plastic structure consist of the yield condition and the equilibrium condition. In practical applications several random model parameters have to be taken into account. This leads to a stochastic optimization problem which cannot be solved using the traditional methods. Instead of that, appropriate (deterministic) substitute problems must be formulated. Here, the topology optimization of frames is considered, where the external load is supposed to be stochastic. The recourse problem will be formulated in general and in the standard form of stochastic linear programming (SLP). After the formulation of the stochastic optimization problem, the recourse problem with discretization and the expected value problem are introduced as representatives of substitute problems. Subsequently, numerical results using these methods are presented. (© 2008 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Optimal design and sensitivity of large scale plane frames for the case of uncertainty

This paper deals with the optimal design and the robustness of large scale plane frames in dependence of their height. Using the first collapse theorem, the necessary and sufficient survival conditions of an elasto-plastic structure consist of the yield condition and the equilibrium condition. The basis for our consideration is provided by a plane n-storey frame which will be increased successively, and which is affected by applied random forces and moments. Taking into account these random applied loads, we get a stochastic structural optimization problem which cannot be solved using the traditional methods. Instead of that, an appropriate (deterministic) substitute problem is formulated. First, the recourse problem will be formulated in general and in the standard form of stochastic linear programming (SLP), and after the formulation of the stochastic optimization problem, the Recourse Problem based on Discretization (RPD) is introduced as a representative of substitute problems. The resulting (large scale) linear program (LP) can be solved efficiently by means of usual LP-solvers. (© 2010 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Production planning with discounting and stochastic demands

A general continuous review production planning problem with stochastic demand is considered. Conditions under which the stochastic problem may be correctly solved using an equivalent deterministic problem are developed. This deterministic problem is known to have the same solution as the stochastic problem. Moreover, conditions are established under which the deterministic equivalent problem differs from a commonly used deterministic approximation to the problem only in the interest rate used in discounting. Thus, solving the stochastic problem is no more difficult than solving a commonly used approximation of the problem.     

A Deterministic Approach To Stochastic Optimal Control With Application To Anticipative Control

Using the decomposition of solution of SDE, we consider the stochastic optimal control problem with anticipative controls as a family of deterministic control problems parametrized by the paths of the driving Wiener process and of a newly introduced Lagrange multiplier stochastic process (nonanticipativity equality constraint). It is shown that the value function of these problems is the unique global solution of a robust equation (random partial differential equation) associated to a linear backward Hamilton-Jacobi-Bellman stochastic partial differential equation (HJB SPDE). This appears as limiting SPDE for a sequence of random HJB PDE's when linear interpolation approximation of the Wiener process is used. Our approach extends the Wong-Zakai type results [20] from SDE to the stochastic dynamic programming equation by showing how this arises as average of the limit of a sequence of deterministic dynamic programming equations. The stochastic characteristics method of Kunita [13] is used to represent the value function. By choosing the Lagrange multiplier equal to its nonanticipative constraint value the usual stochastic (nonanticipative) optimal control and optimal cost are recovered. This suggests a method for solving the anticipative control problems by almost sure deterministic optimal control. We obtain a PDE for the “cost of perfect information” the difference between the cost function of the nonanticipative control problem and the cost of the anticipative problem which satisfies a nonlinear backward HJB SPDE. Poisson bracket conditions are found ensuring this has a global solution. The cost of perfect information is shown to be zero when a Lagrangian submanifold is invariant for the stochastic characteristics. The LQG problem and a nonlinear anticipative control problem are considered as examples in this framework     

Equivalence of laws and null controllability for SPDEs driven by a fractional Brownian motion

We obtain necessary and sufficient conditions for equivalence of law for linear stochastic evolution equations driven by a general Gaussian noise by identifying the suitable space of controls for the corresponding deterministic control problem. This result is applied to semilinear (reaction-diffusion) equations driven by a fractional Brownian motion. We establish the equivalence of continuous dependence of laws of solutions to semilinear equations on the initial datum in the topology of pointwise convergence of measures and null controllability for the corresponding deterministic control problem.     

Stochastic control problems with delay

We consider optimal control problems for systems described by stochastic differential equations with delay. We state conditions for certain classes of such systems under which the stochastic control problems become finite-dimensional. These conditions are illustrated with three applications. First, we solve some linear quadratic problems with delay. Then we find the optimal consumption rate in a financial market with delay. Finally, we solve explicitly a deterministic fluid problem with delay which arises from admission control in ATM communication networks.     

A Chance Constrained Approach to Fractional Programming with Random Numerator

This paper presents a chance constrained programming approach to the problem of maximizing the ratio of two linear functions of decision variables which are subject to linear inequality constraints. The coefficient parameters of the numerator of the objective function are assumed to be random variables with a known multivariate normal probability distribution. A deterministic equivalent of the stochastic linear fractional programming formulation has been obtained and a subsidiary convex program is given to solve the deterministic problem.     

Solution of a class of stochastic linear-convex control problems using deterministic equivalents

In this paper, we identify a new class of stochastic linearconvex optimal control problems, whose solution can be obtained by solving appropriate equivalent deterministic optimal control problems. The term linear-convex is meant to imply that the dynamics is linear and the cost function is convex in the state variables, linear in the control variables, and separable. Moreover, some of the coefficients in the dynamics are allowed to be random and the expectations of the control variables are allowed to be constrained. For any stochastic linear-convex problem, the equivalent deterministic problem is obtained. Furthermore, it is shown that the optimal feedback policy of the stochastic problem is affine in its current state, where the affine transformation depends explicitly on the optimal solution of the equivalent deterministic problem in a simple way. The result is illustrated by its application to a simple stochastic inventory control problem.This research was supported in part by NSERC Grant A4617, by SSHRC Grant 410-83-0888, and by an INRIA Post-Doctoral Fellowship.     

Copositive optimization - Recent developments and applications

Due to its versatility, copositive optimization receives increasing interest in the Operational Research community, and is a rapidly expanding and fertile field of research. It is a special case of conic optimization, which consists of minimizing a linear function over a cone subject to linear constraints. The diversity of copositive formulations in different domains of optimization is impressive, since problem classes both in the continuous and discrete world, as well as both deterministic and stochastic models are covered. Copositivity appears in local and global optimality conditions for quadratic optimization, but can also yield tighter bounds for NP-hard combinatorial optimization problems. Here some of the recent success stories are told, along with principles, algorithms and applications.     

Investigation of a problem of an elastic half space subjected to stochastic temperature distribution at the boundary

The present work investigates the responses of stochastic type temperature distribution applied at the boundary of an elastic medium in the context of thermoelasticity without energy dissipation. We consider an one dimensional problem of half space and assume that the bounding surface of the half space is traction free and is subjected to two types of time dependent temperature distributions which are of stochastic types. In order to compare the results predicted by stochastic temperature distributions with the results of deterministic type temperature distribution, the stochastic type temperature distributions applied at the boundary are taken in such a way that they reduce to the cases of deterministic types as special cases. Integral transform technique along with stochastic calculus is used to solve the problem. The approximated solutions for physical fields like, stress, temperature, displacement etc. are derived for very small values of time where stochastic type boundary conditions are taken to be of white noise type. The problem is further illustrated with graphical representation of numerical solutions of the problem for a particular case. A detailed comparison of the results of stochastic temperature, displacement and stress distributions inside the half space with the corresponding results of deterministic distributions is presented and special features of the effects of stochastic type boundary conditions are highlighted.     

Optimal PID controller of robots under stochastic uncertainty

In practice often it is not possible to specify exact model parameters. Hence, precomputed controller based on some parameter estimates can produce bad results. In this presentation the aim is to combine classical PID control theory and stochastic optimisation methods in order to obtain robust optimal feedback control. The method works with cost functions being minimized and takes into account stochastic parameter varations. After Taylor expansion to calculate expected cost functions and a few transformations an approximate deterministic substitute PID control problem follows. Here, retaining only linear terms, approximation of expectations and variances of the expected cost functions can be calculated explicitly. By means of splines, numerical approximations of the objective function and the differential equations are obtained then. Using stochastic optimization methods, random parameter variations are incorporated into the optimal control process. Hence, robust optimal feedback controls are obtained. (© 2008 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Stochastic MOLP with Incomplete Information: An Interactive Approach with Recourse

Numerous multiobjective linear programming (MOLP) methods have been proposed in the last two decades, but almost all for contexts where the parameters of problems are deterministic. However, in many real situations, parameters of a stochastic nature arise. In this paper, we suppose that the decision-maker is confronted with a situation of partial uncertainty where he possesses incomplete information about the stochastic parameters of the problem, this information allowing him to specify only the limits of variation of these parameters and eventually their central values. For such situations, we propose a multiobjective stochastic linear programming methodology; it implies the transformation of the stochastic objective functions and constraints in order to obtain an equivalent deterministic MOLP problem and the solving of this last problem by an interactive approach derived from the STEM method. Our methodology is illustrated by a didactical example.     

A New Decomposition Technique in Solving Multistage Stochastic Linear Programs by Infeasible Interior Point Methods

Multistage stochastic linear programming (MSLP) is a powerful tool for making decisions under uncertainty. A deterministic equivalent problem of MSLP is a large-scale linear program with nonanticipativity constraints. Recently developed infeasible interior point methods are used to solve the resulting linear program. Technical problems arising from this approach include rank reduction and computation of search directions. The sparsity of the nonanticipativity constraints and the special structure of the problem are exploited by the interior point method. Preliminary numerical results are reported. The study shows that, by combining the infeasible interior point methods and specific decomposition techniques, it is possible to greatly improve the computability of multistage stochastic linear programs.     

Quadratic Control for Stochastic Systems Defined by Evolution Operators and Square Integrable Martingales

The infinite dimensional version of the linear quadratic cost control problem is studied by Curtain and Pritchard [2], Gibson [5] by using Riccati integral equations, instead of differential equations. In the present paper the corresponding stochastic case over a finite horizon is considered. The stochastic perturbations are given by Hilbert valued square integrable martingales and it is shown that the deterministic optimal feedback control is also optimal in the stochastic case. Sufficient conditions are given for the convergence of approximate solutions of optimal control problems.     

Planning logistics operations in the oil industry

In this paper we apply stochastic programming modelling and solution techniques to planning problems for a consortium of oil companies. A multiperiod supply, transformation and distribution scheduling problem—the Depot and Refinery Optimization Problem (DROP)—is formulated for strategic or tactical level planning of the consortium's activities. This deterministic model is used as a basis for implementing a stochastic programming formulation with uncertainty in the product demands and spot supply costs (DROPS), whose solution process utilizes the deterministic equivalent linear programming problem. We employ our STOCHGEN general purpose stochastic problem generator to ‘recreate’ the decision (scenario) tree for the unfolding future as this deterministic equivalent. To project random demands for oil products at different spatial locations into the future and to generate random fluctuations in their future prices/costs a stochastic input data simulator is developed and calibrated to historical industry data. The models are written in the modelling language XPRESS-MP and solved by the XPRESS suite of linear programming solvers. From the viewpoint of implementation of large-scale stochastic programming models this study involves decisions in both space and time and careful revision of the original deterministic formulation. The first part of the paper treats the specification, generation and solution of the deterministic DROP model. The stochastic version of the model (DROPS) and its implementation are studied in detail in the second part and a number of related research questions and implications discussed.     

Pricing problems with a continuum of customers as stochastic Stackelberg games

The pricing problem where a company sells a certain kind of product to a continuum of customers is considered. It is formulated as a stochastic Stackelberg game with nonnested information structure. The inducible region concept, recently developed for deterministic Stackelberg games, is extended to treat the stochastic pricing problem. Necessary and sufficient conditions for a pricing scheme to be optimal are derived, and the pricing problem is solved by first delineating its inducible region, and then solving a constrained optimal control problem.The research work reported here as supported in part by the National Science Foundation under Grant ECS-81-05984, Grant ECS-82-10673, and by the Air Force Office of Scientific Research under AFOSR Grant 80-0098.     

Optimal PID control in case of random initial conditions

The aim of this presentation is to construct a robust optimal PID feedback controller, taking into account stochastic uncertainties in the initial conditions. Usually, a precomputed feedback control is based on exactly known model parameters. However, in practice, often exact information about model parameters and initial values is not given. Hence, having an inital point, which differs from the nominal values, a standard precomputed controller may produce bad results. Supposing now that the probability distribution of the random parameter variations is known, in the following stochastic optimisation methods will be applied in order to obtain robust optimal feedback controls. Taking into account stochastic parameter variations at the initial point, the method works with expected total costs arising from the primary control expenses and the tracking error. Furthermore, the free regulator parameters are selected then such that the expected total costs are minimized. After Taylor expansion to calculate expected cost functions and a few transformations an approximate deterministic substitute control problem follows. Here, retaining only linear terms, approximation of expectations and variances of the expected cost functions can be calculated explicitly. By means of splines, numerical approximations of the objective function and the differential equations are obtained then. Using stochastic optimization methods, random parameter variations are incorporated into the optimal control process. Hence, robust optimal feedback controls are obtained. (© 2009 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

A class of discrete time generalized Riccati equations

In this paper, a class of discrete-time backward non-linear equations defined on some ordered Hilbert spaces of symmetric matrices is considered. The problem of the existence of some global solutions is investigated. The class of considered discrete-time non-linear equations contains, as special cases, a great number of difference Riccati equations both from the deterministic and the stochastic framework. The results proved in the paper provide the sets of necessary and sufficient conditions that guarantee the existence of some special solutions of the considered equations as: the maximal solution, the stabilizing solution and the minimal positive semi-definite solution. These conditions are expressed in terms of the feasibility of some suitable systems of linear matrix inequalities (LMI). One shows that in the case of the equations with periodic coefficients to verify the conditions that guarantee the existence of the maximal or the stabilizing solution, we have to check the solvability of some systems of LMI with a finite number of inequations. The proofs are based on some suitable properties of discrete-time linear equations defined by the positive operators on some ordered Hilbert spaces chosen adequately. The results derived in this paper provide useful conditions that guarantee the existence of the maximal solution or the stabilizing solution for different classes of difference matrix Riccati equations involved in many problems of robust control both in the deterministic and the stochastic framework. The proofs are deterministic and are accessible to the readers less familiarized with the stochastic reasonings.