Filter techniques in shape optimization with CAD-free parametrization
|
|
- Jade Williams
- 6 years ago
- Views:
Transcription
1 6 th World Congresses of Structural and Multidisciplinary Optimization Rio de Janeiro, 30 May - 03 June 2005, Brazil Filter techniques in shape optimization with CAD-free parametrization F. Daoud 1, M. Firl 2, K.-U. Bletzinger 3 daoud@bv.tum.de 1, firl@bv.tum.de 2, kub@bv.tum.de 3, Lehrstuhl für Statik, TU München, Arcisstr. 21, D München, Germany 1,2,3 1. Abstract In this contribution an innovative method for shape optimization with parametrization directly on FE mesh is proposed. The major shortcomings of CAD-free parametrization are discussed and global filter techniques are presented in this context. A global filter for shape optimization with FE-parametrization is developed and applied to numerical examples of basic shapes. The influence of locking phenomena in the context of shape optimization is discussed. Finally a local filter for bead optimization is proposed. 2. Keywords: shape optimization, CAD-free parametrization, FE-parametrization, filter techniques, bead optimization, locking. 3. Introduction Parametrization of the design for shape optimization can be based on CAD models or directly on the FE mesh (CAD-free parametrization). In CAD based design optimization, design variables are chosen among the parameters related to the underlying CAD model. During the optimization process, the update of the finite element mesh is due to changes of the CAD model. In FE based design optimization, nodal coordinates are taken as design variables. This parametrization gives more freedom to the optimization process, since the result of the optimization problem is not restricted to the preselected CAD design space. However, there are many technical facts which in the past made this natural choice of optimization model to appear as prohibitive, as there are: the large number of design variables, the wiggly shape and the inherent problems of FE discretization (e.g. locking) which may lead to mechanically wrong answers. This contribution presents an innovative approach to overcome these shortcomings by a combination of several techniques which results in a powerful and flexible tool for the preliminary design of free form shells. The problem of the large number of variables is overcome by the use of adjoint sensitivity analysis. From the mechanical point of view locking-free elements (shell elements) are indispensable. Not only the correct system response depends on such element formulation, but also the feasible domain of optimization problem is severely affected by locking phenomena, c.f. Camprubi et al. [7], Bischoff et al. [5]. Moreover, it can be shown, that locking phenomena are responsible for shape oscillations on element level, which is the smallest scale that can be represented by a chosen discretization. The main problem of the FE-based design is that wiggly shapes may be obtained as optimal designs. These wiggles are mesh dependent and obviously it is not desirable that the optimum design depends on the coarseness of the mesh used. Therefore, the aim of this research is to provide a tool to control the smoothness of the designs obtained as result of shape optimization and to avoid mesh dependency. Based on the filter proposed by Sigmund in [12] to overcome checkerboard and mesh dependency problems in topology optimization, a filter for shape optimization is proposed. Through this sensitivity filter, not only high frequency waves are smoothed, but also the global wave length of the final design can be controlled by adjusting the filter radius. 4. Filter techniques In general, a shape can be understood as the superposition of several principal shape functions each characterized by a certain wave length. Results of shape optimization in FE based design often of local nature are generated by arbitrary combinations of principal shape functions. As a consequence, successful optimization methods must be able to keep control of the desired mesh independent wave lengths of the resulting shapes by means of filter techniques. Nowadays filter techniques are well established procedures in topology optimization and in many other disciplines, c.f. Bendsoe and Sigmund [3], Bängtsson et al. [2]. These methods are by far more than just pragmatical tools, which do a good job in practice. The research in the last two decades especially in mathematics gives us many technically 1
2 mature methods (for more details see Tikhonov et al. [13]). The above mentioned problem of high-frequency waves in the optimization results arises from the many possible scales of the analyzed design. On the one hand many local minima exist, on the other hand during the optimization process long waves of possible solutions compete with divers short waves or wiggles (of other scales). Exactly at this point the filter techniques are engaged to selectively affect the desired wave length and to filter it out. The aim is to develop a procedure, by which smoothness and curvature of the optimized structures are controlled, involving some user-defined geometrical measurement of the modeled structure. The influence of the high-frequency waves, also called disturbance or noise, can be demonstrated on the basis of the one-dimensional case. A cut-out (Figure 1) shows a long wave as well as other superposed short ones. Figure 1: Influence of high-frequency waves on gradients If the oscillating part of the solution is undesirable, then it is obvious that the wiggles (the superposed short waves) influence the sensitivity severely, e.g. at Point P in Figure 1 we notice even a change of sign. It is obvious that the sensitivity of the undisturbed curve deviates from the value of the disturbed one, which definitely affects the whole optimization procedure, exceedingly in gradient based algorithms. The task now is to get rid of the noisy wiggles systematically. An idea how to control waves in various scales in a noisy system response can be borrowed from signal technologies, where the waves, before transmission, are transformed by means of Fourier transformation or wavelet transformation. It is common technique to approximate a function f by superposing periodical shape functions N i (x) and the corresponding amplitudes or function values f i, respectively f(x) f(x) = i N i (x) f i (1) Beside the shape functions, there are also the dual functions D i (x), which do the opposite job, namely filtering the values f i out. The dual function is defined as follows 1 N i (ξ) D j (ξ) dω = δ ij (2) Ω Ω(x) The integration domain Ω depends on the used shape and dual functions. It usually represents the smallest wave length described in the approximation. If we apply the dual function on the approximation f(x), the evaluations of the original function at the sampling points f i are exactly filtered out. On the other hand, if the dual function is applied on the original (noisy) function f(x), the oscillations, which cannot be represented by the shape functions, will be eliminated, as their wave length is smaller than the radius of Ω. A smoothed function evaluation f m is produced. 1 f(ξ) D(ξ) dω = f m (x) (3) Ω Ω(x) 2
3 The integration, Eq (3), over the domain Ω can be interpreted as a way to determine the weighted mean value at each sampling point, which is usually located in the center of the integration domain. This issue motivates a smoothing method for shape optimization. The choice of a suitable dual function and the discretization of Eq (3), where the sensitivities are filtered instead of the function itself, yields f m x k e = n i=1 D i f x k i n (4) D i n number of nodal points inside the integration domain Ω (including the center e) D i = D(R R ei ) weighting function concerning spacial and topological distance R smoothing radius representing Ω R ei Distance between sampling point e (the center) and node i. k iterator for variable coordinates at the node (e.g. X, Y and Z). i=1 This approach is applied in a similar manner to topology optimization, c.f. Bendsoe and Sigmund [3]. Figure 2: spacial (geometrical) and topological distance. Figure 3: element and node levels arround the center. The radius of Ω controls the smallest desired wave length (size of wiggles) in the optimized structure. An additional difficulty in shape optimization concerns the evaluation of distances on the curved surfaces, as the spacial distance only is not enough to determine the correct weighting factors in the sphere of influence of the filter radius. This fact can be demonstrated by means of an intersection curve on a curved surface (Figure 2). Taking point A as a sampling point (point at which the sensitivity is filtered), we notice that node 1, from the topological point of view, is closer to node A than node 2, and accordingly node 1 should possess a greater weighting factor. However, if only the direct distance is considered, then node 2 will get the bigger weight. One possibility is to measure the distances on the curved surface, which is difficult to perform, particularly for discretized structures. To avoid measuring on the curved surface, we introduce a hierarchical weighting system. The FE nodes around the sampling point are classified in levels. Each level gets a special topological weighting factor, as depicted in Figure 3. The final weighting factor is a combination including both information the topological distance and the spacial one. In this way we keep the algorithmic effort low, by scanning the a priori created node map (in the input phase) within the filter radius. As a result of the filtering process the smallest (shortest) waves in the optimized structures are not controlled by discretization parameters (the size of the finite elements), but rather by user defined measurements like the filter radius. 5. Numerical examples The effects of the presented filter will be demonstrated on several exemplary but meaningful benchmarks of minimizing strain energy (compliance). The examples deal with some basic shapes to show the shortcomings of CAD-free shape optimization without filter techniques. The comparison with the results after applying the proposed filter shows the benefits of the method. The following examples are restricted to the static and linear case. Furthermore a linear elastic material law is applied. 3
4 (a) Example 1 (b) Example 2 (c) Example 3 Figure 4: Start geometry (a) Example 1 (b) Example 2 (c) Example 3 Figure 5: Optimization results without filter (a) Example 1 (b) Example 2 (c) Example 3 Figure 6: Optimization results with filter radius equals half the span length In the first example a Navier-supported circular plate under self-weight is optimized. Since the exact plane (horizontal) position as a start geometry results in a maximum of strain energy, as the whole energy is produced in terms of bending, some disturbance in form of start deflection is necessary. The problems involves 294 design variables which are the Z-coordinates (vertical movement) of all nodes except for the supported ones. As the state of minimal strain energy is reached when the load is carried by membrane forces, we expect a shape which has a tendency towards a minimal surface. Locking-free DSG linear shell elements are used, c.f. Bletzinger et al. [6]. The sensitivity analysis is performed by the adjoint method which reduces the calculation cost severely [10]. In the second example a Navier-supported square plate under point loading in the center is optimized. The same remark about start geometry applies also in this case. For parametrization of the geometry 436 design variables are used. In the third example the start geometry is a rectangular container under internal pressure, where we expect a cylindrical shape as a solution of optimization. In this example 438 design variables are used. The optimization results in Figure 6 show the performance of the proposed method. The wiggles of wave lengths smaller than the filter ra4
5 dius (half the span length) are eliminated, as they cannot be represented by the chosen shape functions. The smoothness of the resulting structures is controlled by means of a user-defined parameter derived from the optimized model (filter radius) and is mesh independent which is the main aim of this approach. 6. Locking in the context of shape optimization Locking phenomena are practically as old as the finite element method itself. Already in the early days of FEM the existance of locking and its influence on the results of structural analysis was noticed and investigated [1, 4]. However, the effect of the various locking phenomena on the optimization process is subject of actual research. feasible domain (Locking) feasible domain (Locking-free) VF Optimum (Locking) Optimum (Locking-free) Figure 7: Influence of locking on the feasible domain. The problems of displacement elements, or any other finite element formulation suffering from locking are presented in the following. In [5] the influence of locking on the formulation of the optimization problem is studied, with the focus on the consequences for the feasible domain. It can be observed, that every function and functional depending on the strain energy and/or the vector of nodal displacements, contain a defect due to locking. As a consequence, any approach based on such an element formulation will include the same phenomena. It can be shown, that locking modifies the functions involved in optimization process, especially stress constraints, which restrict the search space, so that the feasible domain is severely affected [8]. Changes in the feasible domain can cause infeasible optimization results as can be seen exemplarily in Figure 7. The irregular influence of locking on the structural response, particularly on the sensitivities, is responsible for the highfrequency oscillations the wiggles of smallest scale on element level in the results of shape optimization of freeform shells. Due to locking the nodal displacements (as state variables of structural analysis) are always underestimated, however the derived fields as strains and stresses show the typical oscillations on element level (see Figure 8). These oscillations are reproduced in the sensitivity analysis, that provides the gradients for gradient-based optimization algorithms. This inherent waviness enhances the wiggly shapes during optimization. The result of shape optimization of a Naviersupported square plate under uniform surface load is depicted in Figure 9. The result of optimization with displacement elements (Figure 9(a)) manifests the fact, that the oscillating structural response on element level due to locking leads to wiggles in the same scale. The usage of locking-free DSG shell elements however, cancels the wiggles on element level (Figure 9(b)). The wave control on system level The main locking phenomenon investigated in this work is the transverse shear locking. Objective function: minimization of strain energy. Discrete Strain Gap [6, 11]. 5
6 must be performed additionally e.g. by means of the presented filter techniques. (a) Distribution of shear forces for displacement elements. (b) Distribution of shear forces for DSG elements. Figure 8: Influence of locking on shape optimization. (a) Optimization with displacement elements (without filter techniques). (b) Optimization with DSG elements (without filter techniques). Figure 9: Influence of locking on waviness of the solution. 7. Bead optimization The proposed filter technique represents a global filter, which operates structure-wide to homogenize the structural response. It also defines an interface between the designer and the optimized structure, where the scales of waves, which should remain in the optimization process, can be determined independent of mesh coarseness. The basic idea of the bead optimization algorithm is to localize the influence of filter to special regions, the so-called hot spots, where local shape modifications improve the mechanical properties of the structure significantly. These special regions form the bead bodies during optimization. To determine the position and orientation of the hot spots a suitable bead criterion is needed. In [9] Emmrich suggests to place the beads along the trajectories of principle moments. In this work a generic approach basing on an appropriate Sensitivity analysis is applied to achieve the desired improvement by means of local shape alterations. In the case of the minimization of strain energy as objective function the projection of the their sensitivities on the director vector at each FE node is taken as a mechanichal quantity to determine, if 6
7 a FE node belongs to the bead body or not. After sensitivity analysis the gradients are denoised with a global filter and sorted in descending order to choose the candidates for building the bead centers. At each candidate for a bead center the magnitude of the projected sensitivity at the node must exceed a specified limit. By experience the follwoing value s f lim provides a good reference: f s f lim s i s f lim = s f mean + 1 n s f mean = 1 n n i=1 n i=1 ( f s i s f mean ) 2 f s i (5) (a) 1. Iteration. (b) 2. Iteration. (c) After bead fusion. Figure 10: Bead optimization of a cantilever plate (clamped on the right hand side). From the sorted list of candidates the node on top (the highest sensitivity) is marked as the first bead center. The remaining candidates must additionally satisfy the following condition, bevor they can declared as addtional centers: Dist(Cen j, Cand i ) 2αR b (6) with Cenj FE nodes, which are already marked as bead centers Cand i analysed candidate Dist(Cen j, Cand i ) distance between the node Cand i and all known centers Cen j α interaction factor R b influence radius of the bead. The bead optimization of a cantilever plate (Figure 10) shows the functionality of the proposed local filter. The convergence in this case is achieved, when the bead pattern converges, i.e. a stable pattern with no more bead centers is maintained. In the example of the clamped plate convergence is gotton after two iterations. After convergence the structure consists of many separate contiguous beads (Figure 10(b)). The question arises, which of the fugues between the beads should better removed and which should stay. In all presented numerical results the obective is the minimization of strain energy. 7
8 (a) Start Geometry. (b) Stresses in the Start Geometry. (c) Stresses in bead structure. Figure 11: Bead optimization of an uniaxial plate. The answer is attained by an additional sensitivity analysis for the FE nodes lying on the boundary between to beads. The projection of the sensitivity information (of fugue nodes) on the director vector and the comparison with the sensitivities at the involved bead centers provides a mechanical quatity to decide, if the boundary is to be removed or not, as can be seen in Figure 10(c). In the Figures 11 and 12 some numerical examples show the application and effect of the proposed bead optimization algorithm based on filter techniques and model paramertization directly on the FE mesh. (a) Sideview. (b) Topview. (c) Perspective view. 8. References Figure 12: Bead optimization of a container subjected to internal pressure. [1] I. Babuška and M. Suri. On locking and robustness in the finite element method. Technical Report BN-1112, Inst. for Physical Sci and Tech., University of Maryland, College Park Campus, May [2] Erik Bängtsson, Daniel Noreland, and Martin Berggren. Shape optimization of an acoustic horn. Computer methods in applied mechanics and engineering, (192): , [3] M. P. Bendsøe and O. Sigmund. Topology optimization. Theory, methods and applications. Springer Verlag Berlin, Heidelberg, New York,
9 [4] M. Bischoff. Theorie und Numerik einer dreidimensionalen Schalenformulierung. PhD thesis, Institut für Baustatik, Universität Stuttgart, [5] Manfred Bischoff, Natalia Camprubi, and Kai-Uwe Bletzinger. Shape optimization of shells and locking. Computers & Structures, (82): , [6] K.-U. Bletzinger, B. Bischoff, and E. Ramm. A unified approach for shear-locking free triangular and rectangular shell finite elements. Computers & Structures, (75): , [7] N. Camprubi, M. Bischoff, and F. Koschnick. High quality structural response and sensitivity analysis in shape optimal design of shells. In WCSMO 5th World Congress of Structural and Multidisciplinary Optimization, Lidodi Vesolo - Venice, 19 May [8] Natalia Camprubi. Design and Analysis in Shape Optimization of Shells. PhD thesis, Lehrstuhl für Statik, Fakultät für Bauingenieur- und Vermessungswesen, Technische Universität München, [9] Dieter Emmrich. Methodische Sickenentwicklung. PhD thesis, Institut für Fahrzeugtechnik,Technische Hochschule Karlsruhe, [10] S. Kimmich. Strukturoptimierung und Sensibilitätsanalyse mit finiten Elementen. PhD thesis, Institut für Baustatik der Universität Stuttgart, [11] Frank Koschnick. Geometrische Locking-Effekte bei Finiten Elementen und ein allgemeines Konzept zu ihrer Vermeidung. PhD thesis, Lehrstuhl für Statik, Fakultät für Bau- und Vermessungswesen der Technischen Universität München, [12] O. Sigmund and P. Peterson. Numerical instabilities in topology optimization: A survey on procedures dealing with checkerboards, mesh-dependencies and local minima. Structural Optimization, (16):68 75, [13] A. N. Tikhonov, A. V. Stepanov, and A. G. Yagola. Numerical methods for the solution of ill-posed problems. Kluwer Academic Publishers,
Free-Form Shape Optimization using CAD Models
Free-Form Shape Optimization using CAD Models D. Baumgärtner 1, M. Breitenberger 1, K.-U. Bletzinger 1 1 Lehrstuhl für Statik, Technische Universität München (TUM), Arcisstraße 21, D-80333 München 1 Motivation
More informationModeling of Shells with Three-dimensional Finite Elements
Universität Stuttgart Fakultät für Bau und- Umweltingenieurwissenschaften Baustatik und Baudynamik Modeling of Shells with Three-dimensional Finite Elements Manfred Bischoff Institute of Structural Mechanics
More informationGuidelines for proper use of Plate elements
Guidelines for proper use of Plate elements In structural analysis using finite element method, the analysis model is created by dividing the entire structure into finite elements. This procedure is known
More informationLOCAL STRESS ANALYSIS OF STIFFENED SHELLS USING MSC/NASTRAN S SHELL AND BEAM p-elements
LOCAL STRESS ANALYSIS OF STIFFENED SHELLS USING MSC/NASTRAN S SHELL AND BEAM p-elements Sanjay Patel, Claus Hoff, Mark Gwillim The MacNeal-Schwendler Corporation Abstract In large finite element models
More informationA NURBS-BASED APPROACH FOR SHAPE AND TOPOLOGY OPTIMIZATION OF FLOW DOMAINS
6th European Conference on Computational Mechanics (ECCM 6) 7th European Conference on Computational Fluid Dynamics (ECFD 7) 11 15 June 2018, Glasgow, UK A NURBS-BASED APPROACH FOR SHAPE AND TOPOLOGY OPTIMIZATION
More informationFinite Element Method. Chapter 7. Practical considerations in FEM modeling
Finite Element Method Chapter 7 Practical considerations in FEM modeling Finite Element Modeling General Consideration The following are some of the difficult tasks (or decisions) that face the engineer
More informationA nodal based evolutionary structural optimisation algorithm
Computer Aided Optimum Design in Engineering IX 55 A dal based evolutionary structural optimisation algorithm Y.-M. Chen 1, A. J. Keane 2 & C. Hsiao 1 1 ational Space Program Office (SPO), Taiwan 2 Computational
More informationINTERNATIONAL JOURNAL OF CIVIL AND STRUCTURAL ENGINEERING Volume 1, No 2, 2010
Thickness Optimization of Vibrating Shells for Minimum Volume Mallika.A 1,Ramana Rao.N.V 2 1 Associate Professor, Dept. of Civil Engineering, VNR Vignana Jyothi Institute of Engineering and Technology,Hyderabad,Andhra
More informationRegenerating CAD Models with OpenCASCADE and pythonocc from Numerical Models with Application to Shape Optimization
Proceedings of the 7th GACM Colloquium on Computational Mechanics for Young Scientists from Academia and Industry October 11-13, 2017 in Stuttgart, Germany Regenerating CAD Models with OpenCASCADE and
More informationApplication of Finite Volume Method for Structural Analysis
Application of Finite Volume Method for Structural Analysis Saeed-Reza Sabbagh-Yazdi and Milad Bayatlou Associate Professor, Civil Engineering Department of KNToosi University of Technology, PostGraduate
More informationPATCH TEST OF HEXAHEDRAL ELEMENT
Annual Report of ADVENTURE Project ADV-99- (999) PATCH TEST OF HEXAHEDRAL ELEMENT Yoshikazu ISHIHARA * and Hirohisa NOGUCHI * * Mitsubishi Research Institute, Inc. e-mail: y-ishi@mri.co.jp * Department
More informationA Multiple Constraint Approach for Finite Element Analysis of Moment Frames with Radius-cut RBS Connections
A Multiple Constraint Approach for Finite Element Analysis of Moment Frames with Radius-cut RBS Connections Dawit Hailu +, Adil Zekaria ++, Samuel Kinde +++ ABSTRACT After the 1994 Northridge earthquake
More informationFinite Element Analysis Prof. Dr. B. N. Rao Department of Civil Engineering Indian Institute of Technology, Madras. Lecture - 36
Finite Element Analysis Prof. Dr. B. N. Rao Department of Civil Engineering Indian Institute of Technology, Madras Lecture - 36 In last class, we have derived element equations for two d elasticity problems
More informationKEYWORDS Non-parametric optimization, Parametric Optimization, Design of Experiments, Response Surface Modelling, Multidisciplinary Optimization
Session H2.5 OVERVIEW ON OPTIMIZATION METHODS Ioannis Nitsopoulos *, Boris Lauber FE-DESIGN GmbH, Germany KEYWORDS Non-parametric optimization, Parametric Optimization, Design of Experiments, Response
More informationThe Level Set Method applied to Structural Topology Optimization
The Level Set Method applied to Structural Topology Optimization Dr Peter Dunning 22-Jan-2013 Structural Optimization Sizing Optimization Shape Optimization Increasing: No. design variables Opportunity
More informationSizing Optimization for Industrial Applications
11 th World Congress on Structural and Multidisciplinary Optimisation 07 th -12 th, June 2015, Sydney Australia Sizing Optimization for Industrial Applications Miguel A.A.S: Matos 1, Peter M. Clausen 2,
More informationModeling Strategies for Dynamic Finite Element Cask Analyses
Session A Package Analysis: Structural Analysis - Modeling Modeling Strategies for Dynamic Finite Element Cask Analyses Uwe Zencker, Günter Wieser, Linan Qiao, Christian Protz BAM Federal Institute for
More informationCHAPTER 4. Numerical Models. descriptions of the boundary conditions, element types, validation, and the force
CHAPTER 4 Numerical Models This chapter presents the development of numerical models for sandwich beams/plates subjected to four-point bending and the hydromat test system. Detailed descriptions of the
More informationCritical study of design parameterization in topology optimization; The influence of design parameterization on local minima
2 nd International Conference on Engineering Optimization September 6-9, 21, Lisbon, Portugal Critical study of design parameterization in topology optimization; The influence of design parameterization
More informationFirst Order Analysis for Automotive Body Structure Design Using Excel
Special Issue First Order Analysis 1 Research Report First Order Analysis for Automotive Body Structure Design Using Excel Hidekazu Nishigaki CAE numerically estimates the performance of automobiles and
More informationEfficient Robust Shape Optimization for Crashworthiness
10 th World Congress on Structural and Multidisciplinary Optimization May 19-24, 2013, Orlando, Florida, USA Efficient Robust Shape Optimization for Crashworthiness Milan Rayamajhi 1, Stephan Hunkeler
More informationAn explicit feature control approach in structural topology optimization
th World Congress on Structural and Multidisciplinary Optimisation 07 th -2 th, June 205, Sydney Australia An explicit feature control approach in structural topology optimization Weisheng Zhang, Xu Guo
More informationDESIGN AND ANALYSIS OF MEMBRANE STRUCTURES IN FEM-BASED SOFTWARE MASTER THESIS
DESIGN AND ANALYSIS OF MEMBRANE STRUCTURES IN FEM-BASED SOFTWARE MASTER THESIS ARCHINEER INSTITUTES FOR MEMBRANE AND SHELL TECHNOLOGIES, BUILDING AND REAL ESTATE e.v. ANHALT UNIVERSITY OF APPLIED SCIENCES
More informationGlobal and clustered approaches for stress constrained topology optimization and deactivation of design variables
th World Congress on Structural and Multidisciplinary Optimization May 9-24, 23, Orlando, Florida, USA Global and clustered approaches for stress constrained topology optimization and deactivation of design
More informationLevel-set and ALE Based Topology Optimization Using Nonlinear Programming
10 th World Congress on Structural and Multidisciplinary Optimization May 19-24, 2013, Orlando, Florida, USA Level-set and ALE Based Topology Optimization Using Nonlinear Programming Shintaro Yamasaki
More informationElastic Bands: Connecting Path Planning and Control
Elastic Bands: Connecting Path Planning and Control Sean Quinlan and Oussama Khatib Robotics Laboratory Computer Science Department Stanford University Abstract Elastic bands are proposed as the basis
More informationTowards lifetime optimization of hanger connection plates for steel arch bridges
1213 Towards lifetime optimization of hanger connection plates for steel arch bridges M. Baitsch*, D. Hartmann Ruhr-University of Bochum, Germany Abstract Deterioration and damage in hanger connection
More informationSIMULATION OF CARBON-ROVING-STRUCTURES EXTREME LIGHT AND STRONG BY FILAMENT WOUND REINFORCEMENT
SIMULATION OF CARBON-ROVING-STRUCTURES EXTREME LIGHT AND STRONG BY FILAMENT WOUND REINFORCEMENT 1 Dirk Dreißig *, 2 Peter Faßbänder, 1 Ulrich Hindenlang 1 Lasso Ingenieurgesellschaft mbh, Germany, 2 FS
More informationSDC. Engineering Analysis with COSMOSWorks. Paul M. Kurowski Ph.D., P.Eng. SolidWorks 2003 / COSMOSWorks 2003
Engineering Analysis with COSMOSWorks SolidWorks 2003 / COSMOSWorks 2003 Paul M. Kurowski Ph.D., P.Eng. SDC PUBLICATIONS Design Generator, Inc. Schroff Development Corporation www.schroff.com www.schroff-europe.com
More informationStep Change in Design: Exploring Sixty Stent Design Variations Overnight
Step Change in Design: Exploring Sixty Stent Design Variations Overnight Frank Harewood, Ronan Thornton Medtronic Ireland (Galway) Parkmore Business Park West, Ballybrit, Galway, Ireland frank.harewood@medtronic.com
More informationScientific Manual FEM-Design 17.0
Scientific Manual FEM-Design 17. 1.4.6 Calculations considering diaphragms All of the available calculation in FEM-Design can be performed with diaphragms or without diaphragms if the diaphragms were defined
More informationTopology optimization of continuum structures with ɛ-relaxed stress constraints
Topology optimization of continuum structures with ɛ-relaxed stress constraints C.E.M. Guilherme and J.S.O. Fonseca Federal University of Rio Grande do Sul, RS Brazil Abstract Topology optimization of
More informationRevision of the SolidWorks Variable Pressure Simulation Tutorial J.E. Akin, Rice University, Mechanical Engineering. Introduction
Revision of the SolidWorks Variable Pressure Simulation Tutorial J.E. Akin, Rice University, Mechanical Engineering Introduction A SolidWorks simulation tutorial is just intended to illustrate where to
More informationAvailable online at ScienceDirect. Procedia Engineering 155 (2016 )
Available online at www.sciencedirect.com ScienceDirect Procedia Engineering 155 (2016 ) 332 341 International Symposium on "Novel structural skins - Improving sustainability and efficiency through new
More informationINTERNATIONAL JOURNAL OF CIVIL AND STRUCTURAL ENGINEERING Volume 2, No 1, 2011
INTERNATIONAL JOURNAL OF CIVIL AND STRUCTURAL ENGINEERING Volume 2, No 1, 2011 Copyright 2010 All rights reserved Integrated Publishing services Research article ISSN 0976 4399 Topology optimization of
More informationApplication of shape optimization method to artificial leaf design
Design and Nature VI 157 Application of shape optimization method to artificial leaf design M. Shimoda1 & Y. Nakata2 1 Department of Advanced Science and Technology, Toyota Technological Institute, Japan
More informationTopology Optimization of Two Linear Elastic Bodies in Unilateral Contact
2 nd International Conference on Engineering Optimization September 6-9, 2010, Lisbon, Portugal Topology Optimization of Two Linear Elastic Bodies in Unilateral Contact Niclas Strömberg Department of Mechanical
More informationRecent Advances on Higher Order 27-node Hexahedral Element in LS-DYNA
14 th International LS-DYNA Users Conference Session: Simulation Recent Advances on Higher Order 27-node Hexahedral Element in LS-DYNA Hailong Teng Livermore Software Technology Corp. Abstract This paper
More informationCURVILINEAR MESH GENERATION IN 3D
CURVILINEAR MESH GENERATION IN 3D Saikat Dey, Robert M. O'Bara 2 and Mark S. Shephard 2 SFA Inc. / Naval Research Laboratory, Largo, MD., U.S.A., dey@cosmic.nrl.navy.mil 2 Scientific Computation Research
More informationComparative Study of Topological Optimization of Beam and Ring Type Structures under static Loading Condition
Comparative Study of Topological Optimization of Beam and Ring Type Structures under static Loading Condition Vani Taklikar 1, Anadi Misra 2 P.G. Student, Department of Mechanical Engineering, G.B.P.U.A.T,
More informationLarge Rotation Analysis of thin-walled shells
Large Rotation Analysis of thin-walled shells Y. Başar, O. Kintzel Ruhr-University Bochum Bochum, Germany Summary The objective of this contribution is to develop an isoparametric -node shell element by
More informationBehaviour of cold bent glass plates during the shaping process
Behaviour of cold bent glass plates during the shaping process Kyriaki G. DATSIOU *, Mauro OVEREND a * Department of Engineering, University of Cambridge Trumpington Street, Cambridge, CB2 1PZ, UK kd365@cam.ac.uk
More informationOptimization in the Abaqus Environment Using TOSCA
Optimization in the Abaqus Environment Using TOSCA Luca Furbatto, Giovanni Di Lorenzo and Claus B.W. Pedersen Luca Furbatto and Giovanni Di Lorenzo (McLaren Racing), Claus B.W. Pedersen (FE-Design GmbH)
More informationISSN: ISO 9001:2008 Certified International Journal of Engineering and Innovative Technology (IJEIT) Volume 2, Issue 3, September 2012
Mitigation Curves for Determination of Relief Holes to Mitigate Concentration Factor in Thin Plates Loaded Axially for Different Discontinuities Shubhrata Nagpal, S.Sanyal, Nitin Jain Abstract In many
More informationA Numerical Form Finding Method for Minimal Surface of Membrane Structure
10 th World Congress on Structural and Multidisciplinary Optimization May 19-24, 2013, Orlando, Florida, US Numerical Form Finding Method for Minimal Surface of Membrane Structure Koichi Yamane 1, Masatoshi
More informationTOPOLOGY OPTIMIZATION: AN INTRODUCTION. Pierre DUYSINX LTAS Automotive Engineering Academic year
TOPOLOGY OPTIMIZATION: AN INTRODUCTION Pierre DUYSINX LTAS Automotive Engineering Academic year 2017-2018 1 LAY-OUT Introduction Topology problem formulation Problem statement Compliance minimization Homogenization
More informationMesh Processing Pipeline
Mesh Smoothing 1 Mesh Processing Pipeline... Scan Reconstruct Clean Remesh 2 Mesh Quality Visual inspection of sensitive attributes Specular shading Flat Shading Gouraud Shading Phong Shading 3 Mesh Quality
More informationVisualization of errors of finite element solutions P. Beckers, H.G. Zhong, Ph. Andry Aerospace Department, University of Liege, B-4000 Liege, Belgium
Visualization of errors of finite element solutions P. Beckers, H.G. Zhong, Ph. Andry Aerospace Department, University of Liege, B-4000 Liege, Belgium Abstract The aim of this paper is to show how to use
More informationWhat makes Bolt Self-loosening Predictable?
What makes Bolt Self-loosening Predictable? Abstract Dr.-Ing. R. Helfrich, Dr.-Ing. M. Klein (INTES GmbH, Germany) In mechanical engineering, bolts are frequently used as standard fastening elements, which
More informationSTRUCTURAL TOPOLOGY OPTIMIZATION SUBJECTED TO RELAXED STRESS AND DESIGN VARIABLES
STRUCTURAL TOPOLOGY OPTIMIZATION SUBJECTED TO RELAXED STRESS AND DESIGN VARIABLES Hailu Shimels Gebremedhen, Dereje Engida Woldemichael and Fakhruldin M. Hashim Mechanical Engineering Department, Universiti
More informationAcquisition of Line Heating Information for Automatic Plate Forming
Acquisition of Line Heating Information for Automatic Plate Forming Chang Doo Jang 1, Sung Choon Moon 2 and Dae Eun Ko 3 Currently, the promotion of productivity is a significant topic in shipbuilding.
More informationNEW WAVE OF CAD SYSTEMS AND ITS APPLICATION IN DESIGN
Vol 4 No 3 NEW WAVE OF CAD SYSTEMS AND ITS APPLICATION IN DESIGN Ass Lecturer Mahmoud A Hassan Al-Qadisiyah University College of Engineering hasaaneng@yahoocom ABSTRACT This paper provides some lighting
More informationMesh Based Interpolative Coding (MBIC)
Mesh Based Interpolative Coding (MBIC) Eckhart Baum, Joachim Speidel Institut für Nachrichtenübertragung, University of Stuttgart An alternative method to H.6 encoding of moving images at bit rates below
More informationCFD Topology Optimization of Automotive Components
CFD Topology Optimization of Automotive Components Dr.-Ing. Markus Stephan, Dr.-Ing. Dipl.-Phys. Pascal Häußler, Dipl.-Math. Michael Böhm FE-DESIGN GmbH, Karlsruhe, Germany Synopsis Automatic CFD optimization
More informationEnabling Efficient Optimization / Sensitivity and Robustness Analysis for Crashworthiness, NVH, and Multi-disciplinary Concept Assessments
Parametric Modeling of Car Body Structures Enabling Efficient Optimization / Sensitivity and Robustness Analysis for Crashworthiness, NVH, and Multi-disciplinary Concept Assessments White Paper by Dr.
More informationHierarchical Multi-Level-Optimization of crashworthy structures using automatic generated submodels
Hierarchical Multi-Level-Optimization of crashworthy structures using automatic generated submodels Harman Singh 1, Axel Schumacher 1, Carlos J. Falconi D. 2, Alexander F. Walser 2, Sven Trentmann 3, Leyre
More informationInvestigation of the behaviour of single span reinforced concrete historic bridges by using the finite element method
Structural Studies, Repairs and Maintenance of Heritage Architecture XI 279 Investigation of the behaviour of single span reinforced concrete historic bridges by using the finite element method S. B. Yuksel
More informationFluid-structure Interaction by the mixed SPH-FE Method with Application to Aircraft Ditching
Fluid-structure Interaction by the mixed SPH-FE Method with Application to Aircraft Ditching Paul Groenenboom ESI Group Delft, Netherlands Martin Siemann German Aerospace Center (DLR) Stuttgart, Germany
More informationHierarchical topology and shape optimization of crash-loaded profile structures
10 th World Congress on Structural and Multidisciplinary Optimization May 19-24, 2013, Orlando, Florida, USA Hierarchical topology and shape optimization of crash-loaded profile structures Christopher
More informationFinite Element Analysis Using Creo Simulate 4.0
Introduction to Finite Element Analysis Using Creo Simulate 4.0 Randy H. Shih SDC PUBLICATIONS Better Textbooks. Lower Prices. www.sdcpublications.com Powered by TCPDF (www.tcpdf.org) Visit the following
More informationON GRADIENT BASED STRUCTURAL OPTIMIZATION OF A WIND TURBINE BLADE
ON GRADIENT BASED STRUCTURAL OPTIMIZATION OF A WIND TURBINE BLADE E. Lund 1 and J.H. Sjølund 2 1 Department of Materials and Production, Aalborg University Fibigerstræde 16, DK-9220 Aalborg East, Denmark
More informationEngineering Effects of Boundary Conditions (Fixtures and Temperatures) J.E. Akin, Rice University, Mechanical Engineering
Engineering Effects of Boundary Conditions (Fixtures and Temperatures) J.E. Akin, Rice University, Mechanical Engineering Here SolidWorks stress simulation tutorials will be re-visited to show how they
More informationDesign parameterization for topology optimization by intersection of an implicit function
Design parameterization for topology optimization by intersection of an implicit function Peter D. Dunning a, a School of Engineering, University of Aberdeen, Aberdeen AB24 3UE, UK Abstract This paper
More informationCHECKERBOARD PROBLEM IN FINITE ELEMENT BASED TOPOLOGY OPTIMIZATION
CHECKERBOARD PROBLEM IN FINITE ELEMENT BASED TOPOLOGY OPTIMIZATION Avinash Shukla, Anadi Misra, Sunil Kumar Department of Mechanical Engineering, G. B. Pant University of Agriculture and Technology Pantnagar,
More informationAutomatic registration of terrestrial laser scans for geological deformation monitoring
Automatic registration of terrestrial laser scans for geological deformation monitoring Daniel Wujanz 1, Michael Avian 2, Daniel Krueger 1, Frank Neitzel 1 1 Chair of Geodesy and Adjustment Theory, Technische
More informationCombination of Simulation and Optimization in the Design of Electromechanical Systems
Combination of Simulation and Optimization in the Design of Electromechanical Systems Hermann Landes 1, Frank Vogel 2, Wigand Rathmann 2 1 WisSoft, Buckenhof, Germany 2 inutech GmbH, Nuremberg, Germany
More informationarxiv: v1 [math.na] 20 Sep 2016
arxiv:1609.06236v1 [math.na] 20 Sep 2016 A Local Mesh Modification Strategy for Interface Problems with Application to Shape and Topology Optimization P. Gangl 1,2 and U. Langer 3 1 Doctoral Program Comp.
More informationModule 1.6: Distributed Loading of a 2D Cantilever Beam
Module 1.6: Distributed Loading of a 2D Cantilever Beam Table of Contents Page Number Problem Description 2 Theory 2 Geometry 4 Preprocessor 7 Element Type 7 Real Constants and Material Properties 8 Meshing
More informationLoad Balancing for Problems with Good Bisectors, and Applications in Finite Element Simulations
Load Balancing for Problems with Good Bisectors, and Applications in Finite Element Simulations Stefan Bischof, Ralf Ebner, and Thomas Erlebach Institut für Informatik Technische Universität München D-80290
More informationIsogeometric Analysis Application to Car Crash Simulation
Isogeometric Analysis Application to Car Crash Simulation S. Bouabdallah 2, C. Adam 2,3, M. Zarroug 2, H. Maitournam 3 De Vinci Engineering Lab, École Supérieure d Ingénieurs Léonard de Vinci 2 Direction
More informationME 442. Marc/Mentat-2011 Tutorial-1
ME 442 Overview Marc/Mentat-2011 Tutorial-1 The purpose of this tutorial is to introduce the new user to the MSC/MARC/MENTAT finite element program. It should take about one hour to complete. The MARC/MENTAT
More informationModule 1 Lecture Notes 2. Optimization Problem and Model Formulation
Optimization Methods: Introduction and Basic concepts 1 Module 1 Lecture Notes 2 Optimization Problem and Model Formulation Introduction In the previous lecture we studied the evolution of optimization
More informationEXACT BUCKLING SOLUTION OF COMPOSITE WEB/FLANGE ASSEMBLY
EXACT BUCKLING SOLUTION OF COMPOSITE WEB/FLANGE ASSEMBLY J. Sauvé 1*, M. Dubé 1, F. Dervault 2, G. Corriveau 2 1 Ecole de technologie superieure, Montreal, Canada 2 Airframe stress, Advanced Structures,
More informationSmooth finite elements
Smooth finite elements seamless handling of incompressibility, distorted and polygonal meshes; links with equilibrium methods Stéphane Bordas * Nguyen-Xuan Hung ** Nguyen-Dang Hung *** * University of
More informationSHAPE OPTIMIZATION THE EASY WAY: THE METHOD OF TENSILE TRIANGLES
C. Mattheck et al., Int. Journal of Design & Nature. Vol. 2, No. 4 (2007) 301 309 SHAPE OPTIMIZATION THE EASY WAY: THE METHOD OF TENSILE TRIANGLES C. MATTHECK, R. KAPPEL & A. SAUER Forschungszentrum Karlsruhe
More informationTable 1. Results of the questionnaire: Frequency of analysis types (Doubles were possible) [4]
INTERNATIONAL DESIGN CONFERENCE - DESIGN 2004 Dubrovnik, May 18-21, 2004. SUPPORT FOR DESIGNERS USING FEA F. Rieg and F. Koch Keywords: FEA 1. The Objectives In the past, modern computationally based design-
More informationProblem-Adapted Mesh Generation With FEM-Features
INTERNATIONAL DESIGN CONFERENCE - DESIGN 2000 Dubrovnik, May 23-26, 2000. Problem-Adapted Mesh Generation With FEM-Features Dipl.-Ing. Horst Werner, Prof. Dr.-Ing. Christian Weber, cand. ing. Martin Schilke
More informationTopological Machining Fixture Layout Synthesis Using Genetic Algorithms
Topological Machining Fixture Layout Synthesis Using Genetic Algorithms Necmettin Kaya Uludag University, Mechanical Eng. Department, Bursa, Turkey Ferruh Öztürk Uludag University, Mechanical Eng. Department,
More informationCHAPTER 5 FINITE ELEMENT METHOD
CHAPTER 5 FINITE ELEMENT METHOD 5.1 Introduction to Finite Element Method Finite element analysis is a computer based numerical method to deduce engineering structures strength and behaviour. Its use can
More informationSome Aspects for the Simulation of a Non-Linear Problem with Plasticity and Contact
Some Aspects for the Simulation of a Non-Linear Problem with Plasticity and Contact Eduardo Luís Gaertner Marcos Giovani Dropa de Bortoli EMBRACO S.A. Abstract A linear elastic model is often not appropriate
More informationElement Order: Element order refers to the interpolation of an element s nodal results to the interior of the element. This determines how results can
TIPS www.ansys.belcan.com 鲁班人 (http://www.lubanren.com/weblog/) Picking an Element Type For Structural Analysis: by Paul Dufour Picking an element type from the large library of elements in ANSYS can be
More informationFigure 30. Degrees of freedom of flat shell elements
Shell finite elements There are three types of shell finite element; 1) flat elements, 2) elements based on the Sanders-Koiter equations and 3) elements based on reduction of a solid element. Flat elements
More informationINSTRUCTIONAL PLAN L( 3 ) T ( ) P ( ) Instruction Plan Details: DELHI COLLEGE OF TECHNOLOGY & MANAGEMENT(DCTM), PALWAL
DELHI COLLEGE OF TECHNOLOGY & MANAGEMENT(DCTM), PALWAL INSTRUCTIONAL PLAN RECORD NO.: QF/ACD/009 Revision No.: 00 Name of Faculty: Course Title: Theory of elasticity L( 3 ) T ( ) P ( ) Department: Mechanical
More informationThe multi-level hp-method for three-dimensional problems: dynamically changing high-order mesh refinement with arbitrary hanging nodes
The multi-level hp-method for three-dimensional problems: dynamically changing high-order mesh refinement with arbitrary hanging nodes Nils Zander 1, Tino Bog 1, Mohamed Elhaddad 1, Felix Frischmann 1,
More informationStress Analysis of Bolted Joints Part II. Contact and Slip Analysis of a Four Bolt Joint
Modern Mechanical Engineering, 2014, 4, 46-55 Published Online February 2014 (http://www.scirp.org/journal/mme) http://dx.doi.org/10.4236/mme.2014.41006 Stress Analysis of Bolted Joints Part II. Contact
More informationRevised Sheet Metal Simulation, J.E. Akin, Rice University
Revised Sheet Metal Simulation, J.E. Akin, Rice University A SolidWorks simulation tutorial is just intended to illustrate where to find various icons that you would need in a real engineering analysis.
More informationSIMULATION OF METAL FORMING PROCESSES. Konstantin SOLOMONOV a, Victor SVIRIN b
SIMULATION OF METAL FORMING PROCESSES Konstantin SOLOMONOV a, Victor SVIRIN b a Moscow State University of Railway Engineering (Voronezh branch), 75а, Uritskogo street, 394026, Voronezh, Russia, E-mail
More informationSTRUCTURAL & MULTIDISCIPLINARY OPTIMIZATION
STRUCTURAL & MULTIDISCIPLINARY OPTIMIZATION Pierre DUYSINX Patricia TOSSINGS Department of Aerospace and Mechanical Engineering Academic year 2018-2019 1 Course objectives To become familiar with the introduction
More informationModelling Flat Spring Performance Using FEA
Modelling Flat Spring Performance Using FEA Blessing O Fatola, Patrick Keogh and Ben Hicks Department of Mechanical Engineering, University of Corresponding author bf223@bath.ac.uk Abstract. This paper
More informationAlmost Curvature Continuous Fitting of B-Spline Surfaces
Journal for Geometry and Graphics Volume 2 (1998), No. 1, 33 43 Almost Curvature Continuous Fitting of B-Spline Surfaces Márta Szilvási-Nagy Department of Geometry, Mathematical Institute, Technical University
More informationIN-PLANE MATERIAL CONTINUITY FOR THE DISCRETE MATERIAL OPTIMIZATION METHOD
IN-PLANE MATERIAL CONTINUITY FOR THE DISCRETE MATERIAL OPTIMIZATION METHOD René Sørensen1 and Erik Lund2 1,2 Department of Mechanical and Manufacturing Engineering, Aalborg University Fibigerstraede 16,
More informationSheet Metal Forming: Spring-back of hydro mechanical deep drawn parts
4 th European LS-DYNA Users Conference Metal Forming I Sheet Metal Forming: Spring-back of hydro mechanical deep drawn parts Authors: Jens Buchert, University of Applied Sciences, Aalen, Germany David
More informationRobust Design Methodology of Topologically optimized components under the effect of uncertainties
Robust Design Methodology of Topologically optimized components under the effect of uncertainties Joshua Amrith Raj and Arshad Javed Department of Mechanical Engineering, BITS-Pilani Hyderabad Campus,
More informationElastic registration of medical images using finite element meshes
Elastic registration of medical images using finite element meshes Hartwig Grabowski Institute of Real-Time Computer Systems & Robotics, University of Karlsruhe, D-76128 Karlsruhe, Germany. Email: grabow@ira.uka.de
More informationChapter 7 Practical Considerations in Modeling. Chapter 7 Practical Considerations in Modeling
CIVL 7/8117 1/43 Chapter 7 Learning Objectives To present concepts that should be considered when modeling for a situation by the finite element method, such as aspect ratio, symmetry, natural subdivisions,
More informationFB-MULTIPIER vs ADINA VALIDATION MODELING
FB-MULTIPIER vs ADINA VALIDATION MODELING 1. INTRODUCTION 1.1 Purpose of FB-MultiPier Validation testing Performing validation of structural analysis software delineates the capabilities and limitations
More information17th World Conference on Nondestructive Testing, Oct 2008, Shanghai, China
7th World Conference on Nondestructive Testing, 25-28 Oct 2008, Shanghai, China Image Registration Combining Digital Radiography and Computer-Tomography Image Data Frank HEROLD YXLON International X-Ray
More informationModel-based segmentation and recognition from range data
Model-based segmentation and recognition from range data Jan Boehm Institute for Photogrammetry Universität Stuttgart Germany Keywords: range image, segmentation, object recognition, CAD ABSTRACT This
More informationCHAPTER 1. Introduction
ME 475: Computer-Aided Design of Structures 1-1 CHAPTER 1 Introduction 1.1 Analysis versus Design 1.2 Basic Steps in Analysis 1.3 What is the Finite Element Method? 1.4 Geometrical Representation, Discretization
More informationOlivier Brüls. Department of Aerospace and Mechanical Engineering University of Liège
Fully coupled simulation of mechatronic and flexible multibody systems: An extended finite element approach Olivier Brüls Department of Aerospace and Mechanical Engineering University of Liège o.bruls@ulg.ac.be
More information