300 likes | 482 Views
János Lógó Department of Structural Mechanics Budapest University of Technology and Economics Hungary. Optimal topologies in case of probabilistic loading. Introduction , Motivation Mathematical background Assumptions , Mechanical models Parametric Study Conclusions.
E N D
János Lógó Department of Structural Mechanics Budapest University of Technology and Economics Hungary Optimal topologies in case of probabilisticloading IFIP/IIASA/GAMM, 10-12 Dec, 2007, Laxenburg, Austria
Introduction, Motivation • Mathematical background • Assumptions, Mechanical models • Parametric Study • Conclusions IFIP/IIASA/GAMM, 10-12 Dec, 2007, Laxenburg, Austria
Introduction, Motivation K. Marti “Stochastic Optimization Methods”, Springer-Verlag, Berlin-Heidelberg, 2005. K. Marti, “Reliability Analysis of Technical Systems/Structures by means of polyhedral Approximation of the Safe/Unsafe Domain”, GAMM-Mitteilungen, 30, 2, 211-254, 2007. G. Kharmada, N. Olhoff, A. Mohamed, M. Lemaire “Reliability-based Topology Optimization”, Structural and Multidisplinary Optimization, 26, 295-307, 2004. A. Prékopa ”Stochastic Programming”, Akadémia Kiadó and Kluwer, Budapest, Dordrecht, 1995. J. Logo „New Type of Optimality Criteria Method in Case of Probabilistic Loading Conditions”, Mechanics Based Design of Structures and Machines, 35(2), 147-162, 2007. IFIP/IIASA/GAMM, 10-12 Dec, 2007, Laxenburg, Austria
Mathematical background Joint normal distribution Prekopa (1995) -Kataoka (1963): where IFIP/IIASA/GAMM, 10-12 Dec, 2007, Laxenburg, Austria
Mechanical models, Assumptions (3.a) subject to (3.b-d) Stochastically linearized form: IFIP/IIASA/GAMM, 10-12 Dec, 2007, Laxenburg, Austria
Probabilistic compliance constraint Prekopa model: IFIP/IIASA/GAMM, 10-12 Dec, 2007, Laxenburg, Austria
Mechanical model (5.a) subject to (5.b-d) IFIP/IIASA/GAMM, 10-12 Dec, 2007, Laxenburg, Austria
Minimum weight design withstochastically calculated compliance (6.a) subject to (6.b-d) IFIP/IIASA/GAMM, 10-12 Dec, 2007, Laxenburg, Austria
Iterative formulation Determination of the active and passive sets if if if IFIP/IIASA/GAMM, 10-12 Dec, 2007, Laxenburg, Austria
Calculation of the Lagrange-multipliern IFIP/IIASA/GAMM, 10-12 Dec, 2007, Laxenburg, Austria
42 f1=50 f2=50 30 30 30 30 Example 1. 20160 FEs, Poisson’s ratio is 0. The compliance limit is C=410000.q=0.9 IFIP/IIASA/GAMM, 10-12 Dec, 2007, Laxenburg, Austria
Stochastic optimal topology with covariances: 0.1, 0.1, 0, 0 and expected probability value 0.9 IFIP/IIASA/GAMM, 10-12 Dec, 2007, Laxenburg, Austria
Stochastic optimal topology with covariances: 0.1, 0., 0, 0 and expected probability value 0.9 IFIP/IIASA/GAMM, 10-12 Dec, 2007, Laxenburg, Austria
Stochastic optimal topology with covariances: 0, 0.1, 0, 0 and expected probability value 0.9 IFIP/IIASA/GAMM, 10-12 Dec, 2007, Laxenburg, Austria
Stochastic optimal topology with covariances: 0.5, 0.5, 0.0, 0.0 and expected probability value 0.9 IFIP/IIASA/GAMM, 10-12 Dec, 2007, Laxenburg, Austria
Stochastic optimal topology with covariances: 5, 5, 0, 0 and expected probability value 0.9 IFIP/IIASA/GAMM, 10-12 Dec, 2007, Laxenburg, Austria
Probabilistic topologies with variable expected probability with fix mean and covariance values: 0.4, 0.4, 0.01, 0.01 expected probability value: 0.60 IFIP/IIASA/GAMM, 10-12 Dec, 2007, Laxenburg, Austria
Probabilistic topologies with variable expected probability with fix mean and covariance values: 0.4, 0.4, 0.01, 0.01 expected probability value:0.65 IFIP/IIASA/GAMM, 10-12 Dec, 2007, Laxenburg, Austria
Probabilistic topologies with variable expected probability with fix mean and covariance values: 0.4, 0.4, 0.01, 0.01 expected probability value:0.70 IFIP/IIASA/GAMM, 10-12 Dec, 2007, Laxenburg, Austria
Probabilistic topologies with variable expected probability with fix mean and covariance values: 0.4, 0.4, 0.01, 0.01 expected probability value:0.75 IFIP/IIASA/GAMM, 10-12 Dec, 2007, Laxenburg, Austria
Probabilistic topologies with variable expected probability with fix mean and covariance values: 0.4, 0.4, 0.01, 0.01 expected probability value: 0.80 IFIP/IIASA/GAMM, 10-12 Dec, 2007, Laxenburg, Austria
Probabilistic topologies with variable expected probability with fix mean and covariance values: 0.4, 0.4, 0.01, 0.01 expected probability value:0.85 IFIP/IIASA/GAMM, 10-12 Dec, 2007, Laxenburg, Austria
Probabilistic topologies with variable expected probability with fix mean and covariance values: 0.4, 0.4, 0.01, 0.01 expected probability value:0.90 IFIP/IIASA/GAMM, 10-12 Dec, 2007, Laxenburg, Austria
Probabilistic topologies with variable expected probability with fix mean and covariance values: 0.4, 0.4, 0.01, 0.01 expected probability value:0.95 IFIP/IIASA/GAMM, 10-12 Dec, 2007, Laxenburg, Austria
Minimum volumes in function of the expected probability value IFIP/IIASA/GAMM, 10-12 Dec, 2007, Laxenburg, Austria
40 f1=50 f2=50 40 40 Example 2. Cantilever with two forces 24200 FEs, Poisson’s ratio is 0. The compliance limit is C=320000q=0.9 The covariances: , , , IFIP/IIASA/GAMM, 10-12 Dec, 2007, Laxenburg, Austria
Probabilistic topologies with variable expected probability expected probability value:0.75 IFIP/IIASA/GAMM, 10-12 Dec, 2007, Laxenburg, Austria
Probabilistic topologies with variable expected probability expected probability value:0.95 IFIP/IIASA/GAMM, 10-12 Dec, 2007, Laxenburg, Austria
Minimum volumes in function of the expected probability value IFIP/IIASA/GAMM, 10-12 Dec, 2007, Laxenburg, Austria
Conclusions • The probabilistically constrained topology optimization problem was solved • The introduced algorithm provides an iterative tool which allows to use thousands of design variables • The algorithm is rather stable and provides the convergence to reach the optimum. • Needs rather simple computer programming • The covariance values have significant effect for the optimal topology IFIP/IIASA/GAMM, 10-12 Dec, 2007, Laxenburg, Austria