Worst case analysis of mechanical structures by interval methods
|
|
- Oswald Phelps
- 5 years ago
- Views:
Transcription
1 Worst case analysis of mechanical structures by interval methods Arnold Neumaier (University of Vienna, Austria) (joint work with Andrzej Pownuk)
2 Safety Safety studies in structural engineering are supposed to guard against failure in all reasonable situations encountered during the lifetime of a structure. The uncertainty present in assessing what is reasonable has been under discussion since the dawn of history. How can we know what will happen to us when the LORD alone decides? (Proverbs 20:24) Experience has shown that the LORD decides not arbitrarily, but in accordance with the laws of Nature,
3 Thus one way of knowing what will happen to us is to compute worst case bounds on critical response variables, given time-proven finite element models together with worst case bounds on the uncertainties of the input variables. This leads to finite element calculations involving interval parameters.
4 Finite element structural analysis The finite element analysis of mechanical structures amounts in many cases to the solution of a large and sparse linear system with a symmetric, positive definite coefficient matrix. Uncertainties in the material parameters or the execution of a given design result in linear systems with uncertain coefficients. However, as Cramer s rule shows, the uncertainties enter nonlinearly into the equations.
5 FEM equations become nonconvex when data are uncertain Several percent errors in elasticiy modules, and much larger errors in forces are not uncommon Current safety regulation laws require worst case analysis
6 A nonconvex toy problem 1.5 quartic image of nested circles quartic image of nested squares
7 The graphs depict the image of a family of concentric circles and squares (indicating increasing amounts of uncertainty) under the harmless looking function y = F (x) defined by only two nonlinear operations (squaring), u = x 2 x 2 1, z 1 = cx 1 su, z 2 = sx 1 + cu, v = z 2 z 2 1, y 1 = cz 1 sv, y 2 = sz 1 + cv, for c = 0.4, s = 1 c 2. Typical finite element computations are far more complex.
8 Large uncertainties imply large nonlinearities and require new methods. Monte Carlo techniques underestimate worst case Monotonicity-based methods underestimate worst case Local optimization methods underestimate worst case interval analysis, global optimization
9 The goal verification problem A goal verification problem is the quest for verifying that for all parameter combinations in some feasible region, a family of constraints are satisfied, or exhibiting a feasible parameter combination for which some constraint is violated.
10 In many applications involving the design of a structure, a factory, or a machine, it is important that certain goals are met under a variety of conditions that cannot be known in advance. Violating the goals may mean physical danger (risk of failure, often life-threatening) or financial danger (loss of money).
11 Generally, these conditions determining the goal can be specified in terms of a vector x of parameters whose components are unknown. If x is known to lie in a box x = [x, x], and all choices of x C are meaningful scenarios, a deterministic worst case analysis is appropriate.
12 Given the conditions x, the goals are assumed to be expressible in the form of vector constraints F (x) int F for all x x, (1) where F is a vector-valued function, int F is an open box of acceptable values of F, and int denotes the interior. We call (1) the safety constraint(s) since, in the majority of applications, their satisfaction implies that it is safe (for both the designing and the using party) to build and use the structure, factory, or machine, while violation of (1) implies potential danger.
13 The goal verification problem is to show that F (x) int F for all x x, (2) or find a counterexample. The fact that there are infinitely many constraints in (2) makes the problem hard and nonstandard. Past practice is to check (2) only for a number of randomly or systematically generated sample cases. This makes it quite possible that the worst case is overlooked. For safety critical applications, a complete search seems imperative, though it is usually regarded as infeasible to do.
14 Condition (2) is equivalent to checking that the range of F over the box x is contained in the interior of F. In principle, this can be checked by a computation of the range. Using interval analysis, one can often get fairly cheaply enclosures for the range. However, the wrapping effect produces often overly pessimistic enclosures. Moreover, if the computation of F (x) involves the solution of linear systems (as in finite element applications), the enclosure algorithms may even fail completely due to overestimation in intermediate results.
15 Linear systems with uncertain coefficients We consider the uncertain linear system B(x)u(x) = b(x), where the coefficient matrix B(x) and/or the right hand side b(x) depend on a parameter vector x x. For simplicity, we assume a single safety constraint, expressed in terms of the displacement vector u(x), F (x, u(x)) < 0.
16 Centered form approach To solve B(x)u(x) = b(x), we choose a center x 0 and write x = x 0 + s, s s = x x 0. For an arbitrary preconditioning matrix J, we compute enclosures JB(x 0 + Ds) B 0 + D B l s l for all s s, Jb(x 0 + Ds) b 0 + D b l s l for all s s, {u 0 B 0 u 0 b 0 for some B 0 B 0 } u 0.
17 If an interval enclosure ( B 0 + D B l s l ) 1 S exists and then X := S[b 1 B 1 u 0,..., b n B n u 0 ], u(x) u 0 + X(x x 0 ) for all x x. This can be used to check the safety constraint F (x, u(x)) < 0 by another centered form. This works well if uncertainties are only in the right hand side (Muhanna & McMullen), but not for uncertainties in the coefficient matrix.
18 A counterexample ą ć B(x) = x 1 + x 2 x 1 x 2 [0.5, 1.5] A, x x 2 x 1 x 2 x 1 + x 2 [0.5, 1.5] A 1 With x 0 = mid x =, we have s [ 0.5, 0.5] [ 0.5, 0.5] A, J = B(x 0 ) 1 = I, B A, B 1 A, B B 0 + B 1 s 1 + B 2 s 2 [0.5, 1.5] [ 0.5, 0.5] [ 0.5, 0.5] [0.5, 1.5] A contains the singular A although det B(x) = x 1 x 2 > 0 for all x x. 1 A,
19 In higher dimensions, the same problem tends to appear already for much smaller uncertainties. Finite element applications therefore call for a modified approach, which exploits the special form of the finite element equations.
20 Representing FEM matrices In many finite element problems, the only uncertainty in the coefficient matrix is in the element stiffness coefficients x l. (Additional uncertainty in the forces = right hand sides is allowed, too.) For example, in a truss structure, x l = E l a l /L l > 0 where l is the element index, E l the Young modulus describing material properties of the lth bar, a l its cross section area and L l its length.
21 In general, the coefficient matrix depends both on the element stiffness coefficients and on lengths and angles, but if the geometry is assumed fixed then the dependence takes a simple form: m B(x) = x l A T l A l (3) l=1 with extremely sparse matrices A l with few rows. An important special case is where each A l has a single row only. This is the case for truss structures, but not for beams and more complex elements.
22 In the case of truss structures, we may rewrite (3) as B(x) = A T D(x)A, where A = A 1., D(x) = Diag(x 1,..., x m ). A m A is a sparse rectangular matrix, and D(x) diagonal with positive diagonal entries.
23 Unfortunately, even this special form of B(x) is not enough to make the traditional methods work. The example ( x1 + x 2 x 2 x 2 x 1 + x 2 ) u = ( ) 60, x 61 ( ) [1 δ, 1 + δ], [5 δ, 5 + δ] has a coefficient matrix of the form KA T D(x)A with diagonal D(x) with positive diagonal entries, although not constructed from finite elements. We display the boundaries of the solution set for δ = 0.05, 0.1, 0.15, 0.2, 0.25
24
25 The monotonicity method and the vertex sensitivity method fail to find the maximum of u 1 since the maximizer is not a corner. A local optimization may get stuck in the nonglobal minimum. when trying to minimize u 1. Similar problems can be seen in (higher-dimensional) finite element problems.
26 However, coefficient matrices of the form A T D(x)A allow one to iteratively improve the centered form, leading to excellent interval enclosures. We consider more generally uncertain linear systems of the form (K + BDA)u = a + F b, with uncertainties in D and b only.
27 Main Theorem Let D 0 R n n be such that C := (K + BD 0 A) 1 exists, and put d = (D 0 D)v, where v = Au. If there are vectors w 0, w > 0 and w such that w w D 0 D ACB w, w D 0 D ACa + ACF b, then d d := [ αw, αw], α = max i w i w i, (4) and the solution u of (K + BDA)u = a + F b is related to v and d by the equations u = Ca + CF b + CBd, (5) v = ACa + ACF b + ACBd. (6)
28 Proof of (5) (6): Since d = (D 0 D)v, and (K + BDA)u = a + F b, equation u = Ca + CF b + CBd follows from CBd = CB(D 0 D)Au = C(K + BD 0 A)u C(K + BDA)u = u C(a + F b). Multiplication with A gives v = Au = ACa + ACF b + ACBd.
29 Proof of (4): Put β = max i d i /w i, so that d βw, with equality in some component i. w The definition of α = max i i w i implies w αw, hence d = (D 0 D)(ACa + ACF b + ACBd) D 0 D ACa + ACF b + D 0 D ACB βw w + β(w w ) αw + β(w w ). Thus βw i = d i αw i + β(w i w i ), hence βw i αw i. Since w > 0, we conclude that β α, and the desired conclusion d d := [ αw, αw] follows.
30 We now assume D D, b b as interval bounds for the data uncertainties. If we take D 0 = mid D and w as the all-one vector then w := w D 0 D ACB w, w := D 0 D ACa + ACF b, are small and satisfy the required conditions if w > 0 (which is a very moderate condition). By the theorem, d d := [ αw, αw], a narrow box. Using the formulas of the theorem one gets narrow enclosures u for u, v for v and a generally improved enclosure for d. The enclosures can be further improved by iterating this and intersecting with the previously computed enclosures.
31 To get realistic bounds on quantities z = Z(u) (such as safety constraints) dependent linearly or nonlinearly on the solution u of the uncertain linear system, one should intersect the simple enclosure z = Z(u) with the enclosure obtained from the centered form z = Z(CF mid b) + (SCF )(b mid b) + (SCB)d, where S = Z[mid u, u] is a slope matrix for Z.
32 We implemented a Matlab interface to the ANSYS finite element modeling system. The set of ANSYS commands which describe a particular truss structure can be created by using a standard ANSYS GUI. Then, using the algorithm presented here it is possible to calculate the interval solution. Finally it is possible to plot the interval solution in the ANSYS GUI.
33 The following calculations were done in Matlab, using the Intlab package of Rump for the interval calculations. The computer used was an AMD Athlon MP with a 1680 MHz CPU and 3GB memory, running under LINUX.
34
35 10-bay 10-floor truss structure used in work by Muhanna & Mullen
36 Test case: 10-bay 10-floor truss 230 equations (FEM equations for displacements of 121 nodes) 420 two-sided bound constraints (stiffness uncertainties, one per element) total of 650 variables (230 linear and 420 nonlinear ones) stiffness uncertainty α = 5% α defines search region Sizes increase proportional to the number of bays and to the number of floors.
37 The monotonicity method is too slow (exponential amount of work). The vertex sensitivity, although faster and sharper than Monte Carlo methods, underestimates the worst case. The monotonicity assumption is no longer valid at the specified range of uncertainty (already for the 5 5 grid). Traditional interval methods already fails for tiny search regions (uncertainty of α = 0.1%).
38 The relative uncertainty in the stiffness (diagonal entries of D) is varied to be able to assess the degradation of the bounds as the uncertainty increases. The following table lists the widths of the enclosures of the displacements in horizontal (u x ) and vertical (u y ) direction of the upper right corner, and compares it with the widths u x and u y of the enclosure obtained by the element-by-element method of Muhanna & Mullen.
39 uncertainty 1% 2% 3% 4% 5% wid u x wid u x wid u y wid u y With increasing uncertainty, our bounds are increasingly better than those from the element-by-element method.
40 The times for the computation of the enclosures of all components of the solution were times the time needed for the computation of the solution with the midpoint stiffness coefficients. This compares very favorably with Monte Carlo simulations, the monotonicity method and the vertex sensitivity method, although all these only give inner enclosures of the solution set. For smaller examples where it was feasible to run the vertex sensitivity method, the comparison of the resulting inner enclosure with our outer enclosure shows that for the range of uncertainties considered here, our bounds are optimal within a few percent of the width.
41 To see how work and results scale with dimension and uncertainty, we consider a n-bay n-floor truss with variable n, with the same other specifications as in the previous example. truss size n n, n = number of rows of A number of columns of A For n = 60, the memory capacity (about 1GB of free memory) of Matlab, with which we did our tests, was insufficient to store the dense matrix ACA T (with entries) needed in the iteration for our bounds.
42 n = 10, uncertainty 1% 5% 10% 15% 20% 25% u x ux CPU time time ratio n = 20, uncertainty 1% 5% 10% 15% 20% 25% u x ux CPU time time ratio n = 30, uncertainty 1% 5% 10% 15% 20% 25% u x ux CPU time time ratio n = 40, uncertainty 1% 5% 10% 15% 20% 25% u x ux CPU time 2m:55 3m:04 3m:09 3m:09 3m:09 3m:09 time ratio n = 50, uncertainty 1% 5% 10% 15% 20% 25% u x ux CPU time 9m:18 8m:56 9m:13 9m:09 9m:10 9m:10 time ratio
43 In the table, lower and upper bounds are outward rounded. CPU times are given in seconds or minutes : seconds. The time ratio (significant only in the leading digits, due to random variations under repetition) indicates the number of trial point evaluations (using the direct sparse linear solver of Matlab) that can be made in the same time. As can be seen, it is much smaller than the dimensions of the problem, leading to a time advantage over all current methods for approximate worst case analysis.
44 For uncertainties up to about 15%, the enclosures are of high quality. On the other hand, one can see that the enclosures get wider as either the problem size or the uncertainty get larger. In extremal cases (25% uncertainty for n 30), not even the sign of the displacement is guaranteed, probably an artifact of the method caused by overestimation of the worst case.
45 Conclusion The new method gives fast, reliable and accurate worst case bounds for uncertain linear systems. The method scales without difficulties to high-dimensional problems, hence appear to be suitable for many previously untractable instances of worst case analysis in structural engineering. The method works not only for truss structures but applies apply to all kinds of finite element equations. Applications to other mechanical structures are in preparation.
46 Further details can be found in the manuscript A. Neumaier and A. Pownuk, Linear systems with large uncertainties, with applications to truss structures, 2004 A preprint is available from neum/papers.html#linunc
Interval Finite Element Methods for Uncertainty Treatment in Structural Engineering Mechanics Rafi L. Muhanna Georgia Institute of Technology USA
Interval Finite Element Methods for Uncertainty Treatment in Structural Engineering Mechanics Rafi L. Muhanna Georgia Institute of Technology USA Second Scandinavian Workshop on INTERVAL METHODS AND THEIR
More informationUncertainty modeling for robust verifiable design. Arnold Neumaier University of Vienna Vienna, Austria
Uncertainty modeling for robust verifiable design Arnold Neumaier University of Vienna Vienna, Austria Safety Safety studies in structural engineering are supposed to guard against failure in all reasonable
More informationComparison study of the computational methods for eigenvalues IFE analysis
Applied and Computational Mechanics 2 (2008) 157 166 Comparison study of the computational methods for eigenvalues IFE analysis M. Vaško a,,m.sága a,m.handrik a a Department of Applied Mechanics, Faculty
More informationUncertainty is unavoidable in engineering system. Probabilistic approach is the traditional approach
Interval Finite Elements as a Basis for Generalized Models of Uncertainty in Engineering Analysis Rafi L. Muhanna Georgia Institute of Technology USA Collaborative Research Center 528 German National Science
More informationarxiv: v1 [math.na] 28 Jun 2013
A New Operator and Method for Solving Interval Linear Equations Milan Hladík July 1, 2013 arxiv:13066739v1 [mathna] 28 Jun 2013 Abstract We deal with interval linear systems of equations We present a new
More informationSharp Interval Estimates for Finite Element Solutions with Fuzzy Material Properties
Sharp Interval Estimates for Finite Element Solutions with Fuzzy Material Properties R. L. Mullen, F.ASCE Department of Civil Engineering, Case Western Reserve University, Cleveland, Ohio 44106-7201 rlm@po.cwru.edu
More informationGeometric Misfitting in Structures An Interval-Based Approach
Geometric Misfitting in Structures An Interval-Based Approach M. V. Rama Rao Vasavi College of Engineering, Hyderabad - 500 031 INDIA Rafi Muhanna School of Civil and Environmental Engineering Georgia
More informationHowever, reliability analysis is not limited to calculation of the probability of failure.
Probabilistic Analysis probabilistic analysis methods, including the first and second-order reliability methods, Monte Carlo simulation, Importance sampling, Latin Hypercube sampling, and stochastic expansions
More informationFinite Element Structural Analysis using Imprecise Probabilities Based on P-Box Representation
Finite Element Structural Analysis using Imprecise Probabilities Based on P-Box Representation Hao Zhang 1, R. L. Mullen 2, and R. L. Muhanna 3 1 University of Sydney, Sydney 2006, Australia, haozhang@usyd.edu.au
More informationDesign of 2D Elastic Structures with the Interval Parameters
Design of 2D Elastic Structures with the Interval Parameters ANDRZEJ POWNUK The University of Texas at El Paso Department of Mathematical Sciences El Paso, Texas USA andrzej@pownuk.com NAVEEN KUMAR GOUD
More informationFast Algorithms for SDPs derived from the Kalman-Yakubovich-Popov Lemma
Fast Algorithms for SDPs derived from the Kalman-Yakubovich-Popov Lemma Venkataramanan (Ragu) Balakrishnan School of ECE, Purdue University 8 September 2003 European Union RTN Summer School on Multi-Agent
More informationReduction of Random Variables in Structural Reliability Analysis
Reduction of Random Variables in Structural Reliability Analysis S. Adhikari and R. S. Langley Department of Engineering University of Cambridge Trumpington Street Cambridge CB2 1PZ (U.K.) February 21,
More informationStochastic optimization - how to improve computational efficiency?
Stochastic optimization - how to improve computational efficiency? Christian Bucher Center of Mechanics and Structural Dynamics Vienna University of Technology & DYNARDO GmbH, Vienna Presentation at Czech
More informationConstraint propagation for univariate quadratic constraints
Constraint propagation for univariate quadratic constraints Arnold Neumaier Institut für Mathematik, Universität Wien Strudlhofgasse 4, A-1090 Wien, Austria email: Arnold.Neumaier@univie.ac.at WWW: http://www.mat.univie.ac.at/
More informationLinear Matrix Inequalities in Control
Linear Matrix Inequalities in Control Delft Center for Systems and Control (DCSC) Delft University of Technology The Netherlands Department of Electrical Engineering Eindhoven University of Technology
More informationChapter 9: Systems of Equations and Inequalities
Chapter 9: Systems of Equations and Inequalities 9. Systems of Equations Solve the system of equations below. By this we mean, find pair(s) of numbers (x, y) (if possible) that satisfy both equations.
More informationCOMPUTING UNCERTAINTIES IN A DELTA WING COMPOSITE STRUCTURE USING INTERVAL-BASED METHODS
COMPUTING UNCERTAINTIES IN A DELTA WING COMPOSITE STRUCTURE USING INTERVAL-BASED METHODS Franck Delcroix 1, Chris Boyer 2, Vincent Braibant 2 1 Laboratoire de Modélisation et Mécanique des Structures,
More informationMaximum-weighted matching strategies and the application to symmetric indefinite systems
Maximum-weighted matching strategies and the application to symmetric indefinite systems by Stefan Röllin, and Olaf Schenk 2 Technical Report CS-24-7 Department of Computer Science, University of Basel
More informationA Probability-Based Approach to Solving Minesweeper. Step 1: Aryeh Drager MATH 60 Prof. Peter Doyle Dartmouth College 5/31/2012
Aryeh Drager MATH 60 Prof. Peter Doyle Dartmouth College 5/31/2012 A Probability-Based Approach to Solving Minesweeper For this project I have trained MATLAB to play Minesweeper. This training required
More informationComparison between Interval and Fuzzy Load Flow Methods Considering Uncertainty
Comparison between Interval and Fuzzy Load Flow Methods Considering Uncertainty T.Srinivasarao, 2 P.Mallikarajunarao Department of Electrical Engineering, College of Engineering (A), Andhra University,
More informationContent. Department of Mathematics University of Oslo
Chapter: 1 MEK4560 The Finite Element Method in Solid Mechanics II (January 25, 2008) (E-post:torgeiru@math.uio.no) Page 1 of 14 Content 1 Introduction to MEK4560 3 1.1 Minimum Potential energy..............................
More informationDesigning Information Devices and Systems I Spring 2018 Lecture Notes Note Introduction to Linear Algebra the EECS Way
EECS 16A Designing Information Devices and Systems I Spring 018 Lecture Notes Note 1 1.1 Introduction to Linear Algebra the EECS Way In this note, we will teach the basics of linear algebra and relate
More informationKey words. numerical linear algebra, global optimization, validated computing, interval analysis
ON PRECONDITIONERS AND SPLITTING IN THE INTERVAL GAUSS SEIDEL METHOD R. BAKER KEARFOTT AND SIRIPORN HONGTHONG Abstract. Finding bounding sets to solutions to systems of algebraic equations with uncertainties
More informationA Robustness Optimization of SRAM Dynamic Stability by Sensitivity-based Reachability Analysis
ASP-DAC 2014 A Robustness Optimization of SRAM Dynamic Stability by Sensitivity-based Reachability Analysis Yang Song, Sai Manoj P. D. and Hao Yu School of Electrical and Electronic Engineering, Nanyang
More informationIn Chapters 3 and 4 we introduced linear programming
SUPPLEMENT The Simplex Method CD3 In Chapters 3 and 4 we introduced linear programming and showed how models with two variables can be solved graphically. We relied on computer programs (WINQSB, Excel,
More informationIntroduction - Motivation. Many phenomena (physical, chemical, biological, etc.) are model by differential equations. f f(x + h) f(x) (x) = lim
Introduction - Motivation Many phenomena (physical, chemical, biological, etc.) are model by differential equations. Recall the definition of the derivative of f(x) f f(x + h) f(x) (x) = lim. h 0 h Its
More informationDesigning Information Devices and Systems I Fall 2018 Lecture Notes Note Introduction to Linear Algebra the EECS Way
EECS 16A Designing Information Devices and Systems I Fall 018 Lecture Notes Note 1 1.1 Introduction to Linear Algebra the EECS Way In this note, we will teach the basics of linear algebra and relate it
More informationMath 354 Transition graphs and subshifts November 26, 2014
Math 54 Transition graphs and subshifts November 6, 04. Transition graphs Let I be a closed interval in the real line. Suppose F : I I is function. Let I 0, I,..., I N be N closed subintervals in I with
More informationx x2 2 + x3 3 x4 3. Use the divided-difference method to find a polynomial of least degree that fits the values shown: (b)
Numerical Methods - PROBLEMS. The Taylor series, about the origin, for log( + x) is x x2 2 + x3 3 x4 4 + Find an upper bound on the magnitude of the truncation error on the interval x.5 when log( + x)
More informationIntroduction to PDEs: Assignment 9: Finite elements in 1D
Institute of Scientific Computing Technical University Braunschweig Dr. Elmar Zander Dr. Noemi Friedman Winter Semester 213/14 Assignment 9 Due date: 6.2.214 Introduction to PDEs: Assignment 9: Finite
More informationA Hybrid Method for the Wave Equation. beilina
A Hybrid Method for the Wave Equation http://www.math.unibas.ch/ beilina 1 The mathematical model The model problem is the wave equation 2 u t 2 = (a 2 u) + f, x Ω R 3, t > 0, (1) u(x, 0) = 0, x Ω, (2)
More information6-1 Study Guide and Intervention Multivariable Linear Systems and Row Operations
6-1 Study Guide and Intervention Multivariable Linear Systems and Row Operations Gaussian Elimination You can solve a system of linear equations using matrices. Solving a system by transforming it into
More informationOptimization of nonlinear trusses using a displacement-based approach
Struct Multidisc Optim 23, 214 221 Springer-Verlag 2002 Digital Object Identifier (DOI) 10.1007/s00158-002-0179-1 Optimization of nonlinear trusses using a displacement-based approach S. Missoum, Z. Gürdal
More information(x 1 +x 2 )(x 1 x 2 )+(x 2 +x 3 )(x 2 x 3 )+(x 3 +x 1 )(x 3 x 1 ).
CMPSCI611: Verifying Polynomial Identities Lecture 13 Here is a problem that has a polynomial-time randomized solution, but so far no poly-time deterministic solution. Let F be any field and let Q(x 1,...,
More informationSpecial cases of linear programming
Special cases of linear programming Infeasible solution Multiple solution (infinitely many solution) Unbounded solution Degenerated solution Notes on the Simplex tableau 1. The intersection of any basic
More informationMath Numerical Analysis Mid-Term Test Solutions
Math 400 - Numerical Analysis Mid-Term Test Solutions. Short Answers (a) A sufficient and necessary condition for the bisection method to find a root of f(x) on the interval [a,b] is f(a)f(b) < 0 or f(a)
More information1. Introduction. We consider the general global optimization problem
CONSTRUCTION OF VALIDATED UNIQUENESS REGIONS FOR NONLINEAR PROGRAMS IN WHICH CONVEX SUBSPACES HAVE BEEN IDENTIFIED R. BAKER KEARFOTT Abstract. In deterministic global optimization algorithms for constrained
More information7.1 Solving Systems of Equations
Date: Precalculus Notes: Unit 7 Systems of Equations and Matrices 7.1 Solving Systems of Equations Syllabus Objectives: 8.1 The student will solve a given system of equations or system of inequalities.
More informationParallelism in Structured Newton Computations
Parallelism in Structured Newton Computations Thomas F Coleman and Wei u Department of Combinatorics and Optimization University of Waterloo Waterloo, Ontario, Canada N2L 3G1 E-mail: tfcoleman@uwaterlooca
More information6.854J / J Advanced Algorithms Fall 2008
MIT OpenCourseWare http://ocw.mit.edu 6.85J / 8.5J Advanced Algorithms Fall 008 For information about citing these materials or our Terms of Use, visit: http://ocw.mit.edu/terms. 8.5/6.85 Advanced Algorithms
More informationLinear programming I João Carlos Lourenço
Decision Support Models Linear programming I João Carlos Lourenço joao.lourenco@ist.utl.pt Academic year 2012/2013 Readings: Hillier, F.S., Lieberman, G.J., 2010. Introduction to Operations Research, 9th
More informationEfficient Method of Solution of Large Scale Engineering Probelms with Interval Parameters
Efficient Method of Solution of Large Scale Engineering Probelms with Interval Parameters ANDRZEJ POWNUK Department of Civil Engineering, Silesian University of Technology,Gliwice, Poland, pownuk@zeus.polsl.gliwice.pl,
More informationA direct solver for elliptic PDEs in three dimensions based on hierarchical merging of Poincaré-Steklov operators
(1) A direct solver for elliptic PDEs in three dimensions based on hierarchical merging of Poincaré-Steklov operators S. Hao 1, P.G. Martinsson 2 Abstract: A numerical method for variable coefficient elliptic
More informationAMS 529: Finite Element Methods: Fundamentals, Applications, and New Trends
AMS 529: Finite Element Methods: Fundamentals, Applications, and New Trends Lecture 3: Finite Elements in 2-D Xiangmin Jiao SUNY Stony Brook Xiangmin Jiao Finite Element Methods 1 / 18 Outline 1 Boundary
More informationMultiple integrals: Sufficient conditions for a local minimum, Jacobi and Weierstrass-type conditions
Multiple integrals: Sufficient conditions for a local minimum, Jacobi and Weierstrass-type conditions March 6, 2013 Contents 1 Wea second variation 2 1.1 Formulas for variation........................
More information5.7 Cramer's Rule 1. Using Determinants to Solve Systems Assumes the system of two equations in two unknowns
5.7 Cramer's Rule 1. Using Determinants to Solve Systems Assumes the system of two equations in two unknowns (1) possesses the solution and provided that.. The numerators and denominators are recognized
More informationMatrix Assembly in FEA
Matrix Assembly in FEA 1 In Chapter 2, we spoke about how the global matrix equations are assembled in the finite element method. We now want to revisit that discussion and add some details. For example,
More informationCan ANSYS Mechanical Handle My Required Modeling Precision?
Can ANSYS Mechanical Handle My Required Modeling Precision? 2/8/2017 www.telescope-optics.net Recently, PADT received the following inquiry from a scientist interested in using ANSYS for telescope application
More informationthrough any three given points if and only if these points are not collinear.
Discover Parabola Time required 45 minutes Teaching Goals: 1. Students verify that a unique parabola with the equation y = ax + bx+ c, a 0, exists through any three given points if and only if these points
More informationIntroduction, basic but important concepts
Introduction, basic but important concepts Felix Kubler 1 1 DBF, University of Zurich and Swiss Finance Institute October 7, 2017 Felix Kubler Comp.Econ. Gerzensee, Ch1 October 7, 2017 1 / 31 Economics
More informationInstitute for Statics und Dynamics of Structures Fuzzy probabilistic safety assessment
Institute for Statics und Dynamics of Structures Fuzzy probabilistic safety assessment Bernd Möller Fuzzy probabilistic safety assessment randomness fuzzy randomness fuzzyness failure Fuzzy- Probabilistik
More informationPreconditioned Parallel Block Jacobi SVD Algorithm
Parallel Numerics 5, 15-24 M. Vajteršic, R. Trobec, P. Zinterhof, A. Uhl (Eds.) Chapter 2: Matrix Algebra ISBN 961-633-67-8 Preconditioned Parallel Block Jacobi SVD Algorithm Gabriel Okša 1, Marián Vajteršic
More informationInferring from data. Theory of estimators
Inferring from data Theory of estimators 1 Estimators Estimator is any function of the data e(x) used to provide an estimate ( a measurement ) of an unknown parameter. Because estimators are functions
More informationON STOCHASTIC STRUCTURAL TOPOLOGY OPTIMIZATION
ON STOCHASTIC STRUCTURAL TOPOLOGY OPTIMIZATION A. EVGRAFOV and M. PATRIKSSON Department of Mathematics, Chalmers University of Technology, SE-4 96 Göteborg, Sweden J. PETERSSON Department of Mechanical
More informationLinear Algebra. Matrices Operations. Consider, for example, a system of equations such as x + 2y z + 4w = 0, 3x 4y + 2z 6w = 0, x 3y 2z + w = 0.
Matrices Operations Linear Algebra Consider, for example, a system of equations such as x + 2y z + 4w = 0, 3x 4y + 2z 6w = 0, x 3y 2z + w = 0 The rectangular array 1 2 1 4 3 4 2 6 1 3 2 1 in which the
More informationPoisson Solvers. William McLean. April 21, Return to Math3301/Math5315 Common Material.
Poisson Solvers William McLean April 21, 2004 Return to Math3301/Math5315 Common Material 1 Introduction Many problems in applied mathematics lead to a partial differential equation of the form a 2 u +
More informationSolving Linear Systems of Equations
Solving Linear Systems of Equations Gerald Recktenwald Portland State University Mechanical Engineering Department gerry@me.pdx.edu These slides are a supplement to the book Numerical Methods with Matlab:
More informationMechanical Models with Interval Parameters
In K. Gürlebeck L. Hempel C. Könke (Eds.) IKM2003: Digital Proceedings of 16th International Conference on the Applications of Computer Science and Mathematics in Architecture and Civil Engineering, ISSN
More informationGeneralized interval arithmetic on compact matrix Lie groups
myjournal manuscript No. (will be inserted by the editor) Generalized interval arithmetic on compact matrix Lie groups Hermann Schichl, Mihály Csaba Markót, Arnold Neumaier Faculty of Mathematics, University
More informationMTH501- Linear Algebra MCQS MIDTERM EXAMINATION ~ LIBRIANSMINE ~
MTH501- Linear Algebra MCQS MIDTERM EXAMINATION ~ LIBRIANSMINE ~ Question No: 1 (Marks: 1) If for a linear transformation the equation T(x) =0 has only the trivial solution then T is One-to-one Onto Question
More informationTopic 17. Analysis of Algorithms
Topic 17 Analysis of Algorithms Analysis of Algorithms- Review Efficiency of an algorithm can be measured in terms of : Time complexity: a measure of the amount of time required to execute an algorithm
More informationHarbor Creek School District
Unit 1 Days 1-9 Evaluate one-sided two-sided limits, given the graph of a function. Limits, Evaluate limits using tables calculators. Continuity Evaluate limits using direct substitution. Differentiability
More informationLINEAR PROGRAMMING: A GEOMETRIC APPROACH. Copyright Cengage Learning. All rights reserved.
3 LINEAR PROGRAMMING: A GEOMETRIC APPROACH Copyright Cengage Learning. All rights reserved. 3.4 Sensitivity Analysis Copyright Cengage Learning. All rights reserved. Sensitivity Analysis In this section,
More informationFast Verified Solutions of Sparse Linear Systems with H-matrices
Fast Verified Solutions of Sparse Linear Systems with H-matrices A. Minamihata Graduate School of Fundamental Science and Engineering, Waseda University, Tokyo, Japan aminamihata@moegi.waseda.jp K. Sekine
More informationOptimization Problems with Probabilistic Constraints
Optimization Problems with Probabilistic Constraints R. Henrion Weierstrass Institute Berlin 10 th International Conference on Stochastic Programming University of Arizona, Tucson Recommended Reading A.
More informationMULTIPLE CHOICE QUESTIONS DECISION SCIENCE
MULTIPLE CHOICE QUESTIONS DECISION SCIENCE 1. Decision Science approach is a. Multi-disciplinary b. Scientific c. Intuitive 2. For analyzing a problem, decision-makers should study a. Its qualitative aspects
More informationSession-Based Queueing Systems
Session-Based Queueing Systems Modelling, Simulation, and Approximation Jeroen Horters Supervisor VU: Sandjai Bhulai Executive Summary Companies often offer services that require multiple steps on the
More informationElementary and Intermediate Algebra
Elementary and Intermediate Algebra 978-1-63545-106-1 To learn more about all our offerings Visit Knewton.com Source Author(s) (Text or Video) Title(s) Link (where applicable) OpenStax Lynn Marecek, Santa
More informationChapter 5 HIGH ACCURACY CUBIC SPLINE APPROXIMATION FOR TWO DIMENSIONAL QUASI-LINEAR ELLIPTIC BOUNDARY VALUE PROBLEMS
Chapter 5 HIGH ACCURACY CUBIC SPLINE APPROXIMATION FOR TWO DIMENSIONAL QUASI-LINEAR ELLIPTIC BOUNDARY VALUE PROBLEMS 5.1 Introduction When a physical system depends on more than one variable a general
More informationNotes for CS542G (Iterative Solvers for Linear Systems)
Notes for CS542G (Iterative Solvers for Linear Systems) Robert Bridson November 20, 2007 1 The Basics We re now looking at efficient ways to solve the linear system of equations Ax = b where in this course,
More information4 Finite Element Method for Trusses
4 Finite Element Method for Trusses To solve the system of linear equations that arises in IPM, it is necessary to assemble the geometric matrix B a. For the sake of simplicity, the applied force vector
More informationPartial Differential Equations (PDEs)
C H A P T E R Partial Differential Equations (PDEs) 5 A PDE is an equation that contains one or more partial derivatives of an unknown function that depends on at least two variables. Usually one of these
More informationA Parametric Simplex Algorithm for Linear Vector Optimization Problems
A Parametric Simplex Algorithm for Linear Vector Optimization Problems Birgit Rudloff Firdevs Ulus Robert Vanderbei July 9, 2015 Abstract In this paper, a parametric simplex algorithm for solving linear
More informationMIT Algebraic techniques and semidefinite optimization February 14, Lecture 3
MI 6.97 Algebraic techniques and semidefinite optimization February 4, 6 Lecture 3 Lecturer: Pablo A. Parrilo Scribe: Pablo A. Parrilo In this lecture, we will discuss one of the most important applications
More informationChapter 3 Variational Formulation & the Galerkin Method
Institute of Structural Engineering Page 1 Chapter 3 Variational Formulation & the Galerkin Method Institute of Structural Engineering Page 2 Today s Lecture Contents: Introduction Differential formulation
More informationIterative Methods for Linear Systems
Iterative Methods for Linear Systems 1. Introduction: Direct solvers versus iterative solvers In many applications we have to solve a linear system Ax = b with A R n n and b R n given. If n is large the
More informationOverview of course. Introduction to Optimization, DIKU Monday 12 November David Pisinger
Introduction to Optimization, DIKU 007-08 Monday November David Pisinger Lecture What is OR, linear models, standard form, slack form, simplex repetition, graphical interpretation, extreme points, basic
More informationCS 6820 Fall 2014 Lectures, October 3-20, 2014
Analysis of Algorithms Linear Programming Notes CS 6820 Fall 2014 Lectures, October 3-20, 2014 1 Linear programming The linear programming (LP) problem is the following optimization problem. We are given
More informationSensitivity and Reliability Analysis of Nonlinear Frame Structures
Sensitivity and Reliability Analysis of Nonlinear Frame Structures Michael H. Scott Associate Professor School of Civil and Construction Engineering Applied Mathematics and Computation Seminar April 8,
More informationSOLVING POWER AND OPTIMAL POWER FLOW PROBLEMS IN THE PRESENCE OF UNCERTAINTY BY AFFINE ARITHMETIC
SOLVING POWER AND OPTIMAL POWER FLOW PROBLEMS IN THE PRESENCE OF UNCERTAINTY BY AFFINE ARITHMETIC Alfredo Vaccaro RTSI 2015 - September 16-18, 2015, Torino, Italy RESEARCH MOTIVATIONS Power Flow (PF) and
More informationRobotics I. February 6, 2014
Robotics I February 6, 214 Exercise 1 A pan-tilt 1 camera sensor, such as the commercial webcams in Fig. 1, is mounted on the fixed base of a robot manipulator and is used for pointing at a (point-wise)
More informationSimple Examples on Rectangular Domains
84 Chapter 5 Simple Examples on Rectangular Domains In this chapter we consider simple elliptic boundary value problems in rectangular domains in R 2 or R 3 ; our prototype example is the Poisson equation
More informationDistributed Optimization. Song Chong EE, KAIST
Distributed Optimization Song Chong EE, KAIST songchong@kaist.edu Dynamic Programming for Path Planning A path-planning problem consists of a weighted directed graph with a set of n nodes N, directed links
More informationVerification of Nonlinear Hybrid Systems with Ariadne
Verification of Nonlinear Hybrid Systems with Ariadne Luca Geretti and Tiziano Villa June 2, 2016 June 2, 2016 Verona, Italy 1 / 1 Outline June 2, 2016 Verona, Italy 2 / 1 Outline June 2, 2016 Verona,
More informationIterative Methods for Ax=b
1 FUNDAMENTALS 1 Iterative Methods for Ax=b 1 Fundamentals consider the solution of the set of simultaneous equations Ax = b where A is a square matrix, n n and b is a right hand vector. We write the iterative
More informationMatrices. A matrix is a method of writing a set of numbers using rows and columns. Cells in a matrix can be referenced in the form.
Matrices A matrix is a method of writing a set of numbers using rows and columns. 1 2 3 4 3 2 1 5 7 2 5 4 2 0 5 10 12 8 4 9 25 30 1 1 Reading Information from a Matrix Cells in a matrix can be referenced
More information15-780: LinearProgramming
15-780: LinearProgramming J. Zico Kolter February 1-3, 2016 1 Outline Introduction Some linear algebra review Linear programming Simplex algorithm Duality and dual simplex 2 Outline Introduction Some linear
More informationCS6999 Probabilistic Methods in Integer Programming Randomized Rounding Andrew D. Smith April 2003
CS6999 Probabilistic Methods in Integer Programming Randomized Rounding April 2003 Overview 2 Background Randomized Rounding Handling Feasibility Derandomization Advanced Techniques Integer Programming
More informationA TOUR OF LINEAR ALGEBRA FOR JDEP 384H
A TOUR OF LINEAR ALGEBRA FOR JDEP 384H Contents Solving Systems 1 Matrix Arithmetic 3 The Basic Rules of Matrix Arithmetic 4 Norms and Dot Products 5 Norms 5 Dot Products 6 Linear Programming 7 Eigenvectors
More informationEconomics 203: Intermediate Microeconomics. Calculus Review. A function f, is a rule assigning a value y for each value x.
Economics 203: Intermediate Microeconomics Calculus Review Functions, Graphs and Coordinates Econ 203 Calculus Review p. 1 Functions: A function f, is a rule assigning a value y for each value x. The following
More informationS. Freitag, B. T. Cao, J. Ninić & G. Meschke
SFB 837 Interaction Modeling Mechanized Tunneling S Freitag, B T Cao, J Ninić & G Meschke Institute for Structural Mechanics Ruhr University Bochum 1 Content 1 Motivation 2 Process-oriented FE model for
More informationLinear Algebra Homework and Study Guide
Linear Algebra Homework and Study Guide Phil R. Smith, Ph.D. February 28, 20 Homework Problem Sets Organized by Learning Outcomes Test I: Systems of Linear Equations; Matrices Lesson. Give examples of
More informationStress analysis of a stepped bar
Stress analysis of a stepped bar Problem Find the stresses induced in the axially loaded stepped bar shown in Figure. The bar has cross-sectional areas of A ) and A ) over the lengths l ) and l ), respectively.
More informationROBUST DESIGN OF A MECHANICAL LATCH
LA-UR-15-21296 Page 1 of 20 ROBUST DESIGN OF A MECHANICAL LATCH François Hemez Technical Staff Member, XTD-Division, Los Alamos National Laboratory Adjunct Professor, University of California San Diego
More informationResearch Collection. Basics of structural reliability and links with structural design codes FBH Herbsttagung November 22nd, 2013.
Research Collection Presentation Basics of structural reliability and links with structural design codes FBH Herbsttagung November 22nd, 2013 Author(s): Sudret, Bruno Publication Date: 2013 Permanent Link:
More informationValidated Solutions of Initial Value Problems for Parametric ODEs
Validated Solutions of Initial Value Problems for Parametric ODEs Youdong Lin and Mark A. Stadtherr Department of Chemical and Biomolecular Engineering University of Notre Dame, Notre Dame, IN 46556, USA
More informationSOLUTIONS FOR ADMISSIONS TEST IN MATHEMATICS, COMPUTER SCIENCE AND JOINT SCHOOLS WEDNESDAY 31 OCTOBER 2018
SOLUTIONS FOR ADMISSIONS TEST IN MATHEMATICS, COMPUTER SCIENCE AND JOINT SCHOOLS WEDNESDAY OCTOBER 8 Mark Scheme: Each part of Question is worth marks which are awarded solely for the correct answer. Each
More informationDeterministic Global Optimization Algorithm and Nonlinear Dynamics
Deterministic Global Optimization Algorithm and Nonlinear Dynamics C. S. Adjiman and I. Papamichail Centre for Process Systems Engineering Department of Chemical Engineering and Chemical Technology Imperial
More informationWeek 8. 1 LP is easy: the Ellipsoid Method
Week 8 1 LP is easy: the Ellipsoid Method In 1979 Khachyan proved that LP is solvable in polynomial time by a method of shrinking ellipsoids. The running time is polynomial in the number of variables n,
More informationRESPONSE SURFACE METHODS FOR STOCHASTIC STRUCTURAL OPTIMIZATION
Meccanica dei Materiali e delle Strutture Vol. VI (2016), no.1, pp. 99-106 ISSN: 2035-679X Dipartimento di Ingegneria Civile, Ambientale, Aerospaziale, Dei Materiali DICAM RESPONSE SURFACE METHODS FOR
More information