Topology Optimization Using the SIMP Method
|
|
- Dustin Briggs
- 6 years ago
- Views:
Transcription
1 Fabian Wein
2 Introduction Introduction About this document This is a fragment of a talk given interally Intended for engineers and mathematicians SIMP basics Detailed introduction (based on linear elasticity) Optimization vs. simulation Do it yourself with Sigmund s 99-line code Piezoelectric Loupspeaker Extended SIMP model Results
3 Optimization Introduction General optimization Objective function (scalar) Design variable Structural optimization Topology design of truss structures Shape optimization Parametrization Level set metheod (eventually topology gradient) Topology optimization Homogenization SIMP
4 Linear elasticity Classical SIMP Background Hooke s law [σ] = [c 0 ][S] (in Voigt notation: σ = [c 0 ]Bu) with [σ], σ : Cauchy stress tensor [c 0 ] : tensor of elastic modului [S],S : linear strain tensor u : displacement x z y B = 0 y 0 z 0 x 0 0 z y x 0 T : differential operator
5 Strong formulation Classical SIMP Background PDE Find u : Ω R 3 fulfilling B T [c 0 ]Bu = f in Ω with the boundary conditions u = 0 n T [σ] = 0 on Γ s on ΩΓ s
6 Classical SIMP Discrete FEM formulation Background Solve Global System Ku = f with Assembly K = n e e=1 K e ; K e = [k pq ]; k pq = (B) T [c 0 ]B dω Ω e
7 Classical SIMP Proportional stiffness model Incredients Parametrization by design variable Model structure by local stiffness (full and void). Define local stiffness (finite) element wise: ρ = (ρ 1 ρ ne ) T Continuous interpolation with ρ min ρ e 1. Introduce pseudo density ρ [ c e ](ρ) = ρ e [c 0 ]; Ke (ρ) = ρ e K e ; K(ρ)u(ρ) = f
8 Classical SIMP Minimal mean compliance Incredients Different interpretations Maximize stiffness Minimize mean compliance Minimize stored mechanical energy Minimize compliance min ρ J(u(ρ)) = min ρ ft u(ρ) = min ρ u(ρ)t K(ρ)u(ρ)
9 Find derivative Classical SIMP Incredients General optimization procedure Evaluate objective function Find descent direction δ (e.g. gradient) Find step length along δ (line search) Techniques to find descent direction Use gradient free methods Use finite differences Analytical first derivative Analytical second derivative
10 Sensitvity analysis Classical SIMP Incredients Sensitivity analysis provides analytical derivatives Abbreviate ( ) ρ e by ( ) Derive mean compliance f T u J = f T u + f T u = f T u Find J by deriving state condition Ku = f Solve for every u Ku = K u
11 Adjoint method Classical SIMP Incredients The adjoint method is based on the fixed vector λ J = f T u + λ T ( Ku f) J = f T u + λ T ( K u + Ku ) = (f T λ T K)u + λ T K u Solve: Kλ = f = J u J = u T K u The compliance problem is self-adjoint The general adjoint problem can be efficiently solved by (incomplete) LU decomposition
12 Naive approach Classical SIMP Application of the SIMP method Minimize compliance: straight forward, initial design 0.5 min ρ ft u s.th.: Ku = f ρ e [ρ min : 1] note: K e = ρ e K e, K e = K e,
13 Naive approach Classical SIMP Application of the SIMP method Minimize compliance: straight forward, initial design 0.5 min ρ ft u s.th.: Ku = f ρ e [ρ min : 1] note: K e = ρ e K e, K e = K e, The optimal topology is the trivial solution full material
14 Add constraint Classical SIMP Application of the SIMP method Minimize compliance: volume constraint 50% min ρ ft u s.th.: Ω ρ 1 2 V 0
15 Add constraint Classical SIMP Application of the SIMP method Minimize compliance: volume constraint 50% min ρ ft u s.th.: Ω ρ 1 2 V 0 Grey material has no physical interpretation
16 Third try Classical SIMP Application of the SIMP method Minimize compliance: penalize ρ by ρ p with p = 3 min ρ ft u note: K e = ρ 3 e K e, K e = 3ρ 2 e K e,
17 Third try Classical SIMP Application of the SIMP method Minimize compliance: penalize ρ by ρ p with p = 3 min ρ ft u note: K e = ρ 3 e K e, K e = 3ρ 2 e K e, We have a desired 0-1 pattern but checkerboard structure
18 Forth try Classical SIMP Application of the SIMP method Minimize compliance: use averaged gradients min ρ ft u note: K e = ρ i H i i ρ e 3ρe 2 K e i H i with H i = r min dist(e,i)
19 Forth try Classical SIMP Application of the SIMP method Minimize compliance: use averaged gradients min ρ ft u note: K e = ρ i H i i ρ e 3ρe 2 K e i H i with H i = r min dist(e,i) No checkerboards and no mesh dependency
20 Classical SIMP Optimizers Comparison of different optimizers SCPIP (MMA implementation by Ch. Zillober) Optimality Condition (heuristic for SIMP) IPOPT (general second order optimizer)
21 Performance Classical SIMP Optimizers
22 Optimality Condition Classical SIMP Optimizers Optimality Condition: fix-point type update scheme With max{(1 ζ )ρ ek,ρ min } if ρ ek Be η k max{(1 η)ρ ek,ρ min } ρ ek+1 = min{(1 + ζ )ρ ek,1} if min{(1 + ζ )ρ ek,1} ρ ek Be η k ρ ek Be η k else B ek = Λ 1 K e Λ to be found by bisection Step width ζ e.g. 0.2 Damping η e.g. 0.5
23 Extensions to SIMP Multiple loads Complex load vs. multiple load cases For multiple loadcases several problems are averaged Figure: Two loads applied simultaniously (left) and seperatly (right) The left case is instable if the loads are not applied simultaniously
24 Extensions to SIMP Problem specific optimization Multiple loads Now only the left load is applied to the optimized structures Figure: The scaling of the displacement is the same
25 Extensions to SIMP Optimization for arbitrary nodes Synthesis of compliant mechanisms - aka no title Generalizing the compliance to J = l T u with l = ( ) T.
26 Extensions to SIMP Optimization for arbitrary nodes Synthesis of compliant mechanisms - aka no title Generalizing the compliance to J = l T u with l = ( ) T. For this case one has to apply springs to the load and output nodes
27 Harmonic optimization Extensions to SIMP Harmonic optimization Two common approaches Optimize for eigenvalues Perform SIMP with forced vibrations Harmonic excitation Excite with a single frequency Gain steady-state solution in one step Complex numbers Complex FEM system ( K + jω C ω 2 M)u = f S(ω)u = f ST = S
28 Extensions to SIMP Harmonic optimization Harmonic objective functions: J(u(ρ)) R Compliance J = u T f J = R(sign(J)u T S u) J = (u T f) 2 J = 2(u T f)u T S u J = u T R f I u T I f R J = 2R(λ T S u) Sλ = 2 f j J = u T ū J = 2R(λ T S u) Sλ = ū Optimize for output J = u T Lū J = 2R(λ T S u) Sλ = L T ū Optimize for velocity Optimize for coupled quantities
29 Extensions to SIMP Harmonic tranfer functions Harmonic optimization Classical SIMP converges faster than mass to zero ρe 3 if ρ > 0.1 ρ e µ Pedersen (ρ e ) = ρ e µ RAMP (ρ e ) = if ρ q(1 ρ e ) 100 1e+000 8e-001 SIMP idendity RAMP Contribution 6e-001 4e-001 2e-001 0e Design variable
30 Extensions to SIMP (Global) dynamic compliance Harmonic optimization We optimize for u T ū and u T Lū with L selecting f This illustrates general optimization problems One has to know what one wants One might not want what one gets
31 Piezoelectric Model Piezoelectricity Model We couple linear elasticity with electrostatic Material law [σ] = [c E 0 ][S] [e 0] T E, D = [e 0 ][S] + [ε S 0 ]E. E : electric field intensity in V/m D : electric displacement field C/m 2 [σ] : Cauchy stress tensor [S] : linear strain tensor [c E 0 ] : tensor of elastic moduli [c m ] : tensor of elastic moduli [ε S 0 ] : tensor of dielectric permittivities [e 0 ] : tensor of piezoelectric moduli
32 Strong formulation Piezoelectricity Model Loaded Electrode Grounded Electrode Piezoelectric Material Plate Strong formulation Support Find u p : Ω p R 3, u m : Ω m R 3, φ : Ω p R fulfilling ( ) B T [c E 0 ]Bu p + [e 0 ] T Bφ = 0 in Ω p, B T ( [e 0 ]Bu p [ε S 0 ] Bφ ) = 0 in Ω p, B T [c m ]Bu m = 0 in Ω m
33 Strong formulation Piezoelectricity Model Boundary conditions u m = 0 on Γ s, n T p [σ p ] = 0 on Ω p \ Γ g, n T m[σ m ] = 0 on Ω m \ (Γ g Γ s ), n T p [σ p ] = n T m[σ m ] on Γ g, n p = n m on Γ g, u p = u m on Γ g, φ = 0 on Γ g, φ = φ l on Γ l, n T p D = 0 on Ω p \ (Γ l Γ g )
34 FEM formulation Piezoelectricity Model Bilinearforms K e uu = [k uu pq ]; K e uφ = [k uφ pq ]; K e φφ = [k φφ pq ]; kpq uu = k uφ pq = k φφ pq = Ω e ( B u p ) T [c E 0 ]B u q dω, B T p [e 0 ] B q dω, Ω e BT p [ε S 0 ] B q dω. Ω e Global system K um u m K um u p 0 u m 0 K T u m u p K up u p K up φ u p = q u 0 K T u p φ K φφ φ q φ
35 SIMP optimization of piezoelectric devices Piezoelectric SIMP model Model Design variables pseudo density ρ pseudo polarization ρ p Extended material tensors [ c E e ] = µ c (ρ e )[c E 0 ], [ẽ e ] = µ e (ρ e )[e 0 ], [ ε S e ] = µ ε (ρ e )µ p (ρ p e )[ε S 0 ]. ρ e [ρ min ;1] ρ p e [ 1;1] Adjoint PDE for inhomogeneous forward problem No electric excitation but load from objective Gradient reduces with first order elements to J = w T e K e uuu e + w T e K e uφ φ e
36 SIMP optimization of piezoelectric devices Transduction basics Transduction Reciprocal theorem in elasticity t T a u b dγ = t T b u adγ Γ ta Γ tb... by knowing the body response for one load case, we can calculate the displacement at any point of the body caused by another load case.
37 SIMP optimization of piezoelectric devices Transduction basics Transduction Reciprocal theorem in elasticity t T a u b dγ = t T b u adγ Γ ta Γ tb... by knowing the body response for one load case, we can calculate the displacement at any point of the body caused by another load case. Extension to piezoelectricity promisses (Kögl and Silva, 2005):... the conversion of electrical into elastic energy and vice versa.
38 SIMP optimization of piezoelectric devices Transduction in piezoelectricity Transduction Loadcase a: f a 0 Q a = 0 Loadcase b: f b = 0 Q b 0 Loadcase c: f c 0 grounded electrodes Extension to piezoelectricity L ab = φ T b KT uφ u a + φ T a K φφ φ b L ab = û T a K û b. Difference to Kögl and Silva (2005) We have fixed supporting mechanical plate Kögl and Silva have loadcase c and volume constraint Original objective J(ρ) = w lnl ab (1 w)lnl cc, 0 w 1,
39 SIMP optimization of piezoelectric devices Results of mean transduction Transduction The result is the trivial result vanishing material We excite with fixed charge and fixed nodal force Loadcase a (Q a = 0): minimal stiffening (u) and maximal bending (φ) Loadcase b (f b = 0): surface charge density and E tend to infinity Displacement (m) 5.0e e e e e e e e e e-006 mechanical excitation electric excitation Area void/piezo (%) Voltage (V) 1.0e e e e e e-001 mechanical excitation electric excitation Area void/piezo (%) Figure: Parameter study with varying area covered by piezoelectric material
40 SIMP optimization of piezoelectric devices Results Optimize for maximum displacement of ground plate Thickness of layers Piezoelectric layer: 50 µm Supporting layer: 10, 20, 50, 100, 200, 500 µm No volume constraint and no filtering required!
41 SIMP optimization of piezoelectric devices Locking effects Results Very thin elements Shear locking occures when using linear finite elements Straight forward extension of incompatible modes by SIMP Solve local system of higher order ( K e uu K e αu K e )( uα ue K e αα α e ) = ( ) fe 0 Displacement 3.0e e e e e e e+000 linear, with locking linear, locking free quadratic Discretization of edge Volume fraction linear, with locking linear, locking free quadratic Discretization of edge
42 The End End Last comments The optimal solution lays inside the PDE (plus adjoint RHS) Optimization helps to understand systems better Optimization is the next step after simulation Thanks for your time!
C6.1. Topology Optimization of a Piezoelectric Energy Harvester. 1 Introduction. 2 The Composite Model
C6.1 Topology Optimization of a Piezoelectric Energy Harvester Wein, Fabian; Weller, Elias; Albach, Thorsten; Sutor, Alexander; Lerch, Reinhard Department for Sensor Technology, University Erlangen-Nuremberg
More informationTopology Optimization Mathematics for Design
Topology Optimization Mathematics for Design Homogenization Design Method (HMD) Why topology? Change in shape & size may not lead our design criterion for reduction of structural weight. 1 Structural Design
More informationImproved Two-Phase Projection Topology Optimization
10 th World Congress on Structural and Multidisciplinary Optimization May 19-24, 2013, Orlando, Florida, USA Improved Two-Phase Projection Topology Optimization Josephine V. Carstensen and James K. Guest
More informationDesign of in-plane piezoelectric sensors for static response by simultaneously optimizing the host structure and the electrode profile
10 th World Congress on Structural and Multidisciplinary Optimization May 19-24, 2013, Orlando, Florida, USA Design of in-plane piezoelectric sensors for static response by simultaneously optimizing the
More informationFinite Element Method in Geotechnical Engineering
Finite Element Method in Geotechnical Engineering Short Course on + Dynamics Boulder, Colorado January 5-8, 2004 Stein Sture Professor of Civil Engineering University of Colorado at Boulder Contents Steps
More informationContinuum Mechanics and the Finite Element Method
Continuum Mechanics and the Finite Element Method 1 Assignment 2 Due on March 2 nd @ midnight 2 Suppose you want to simulate this The familiar mass-spring system l 0 l y i X y i x Spring length before/after
More informationDesign of piezoelectric microtransducers based on the topology optimization method
Microsyst Technol (2016) 22:1733 1740 DOI 10.1007/s00542-016-2828-0 TECHNICAL PAPER Design of piezoelectric microtransducers based on the topology optimization method David Ruiz 1 Alberto Donoso 1 José
More informationTwo-material Topology Optimization for the Design of Passive Thermal Control Structures
ICAST2014 # 017 Two-material Topology Optimization for the Design of Passive Thermal Control Structures Pierre F. Thurier 1, George A. Lesieutre 2, Mary I. Frecker 3, James H. Adair 4 1 M.S. student, Aerospace
More informationSIMULATION AND OPTIMIZATION OF MEMS PIEZOELECTRIC ENERGY HARVESTER WITH A NON-TRADITIONAL GEOMETRY
SIMULATION AND OPTIMIZATION OF MEMS PIEZOELECTRIC ENERGY HARVESTER WITH A NON-TRADITIONAL GEOMETRY S. Sunithamani 1, P. Lakshmi 1, E. Eba Flora 1 1 Department of EEE, College of Engineering, Anna University,
More information3. Numerical integration
3. Numerical integration... 3. One-dimensional quadratures... 3. Two- and three-dimensional quadratures... 3.3 Exact Integrals for Straight Sided Triangles... 5 3.4 Reduced and Selected Integration...
More informationTheoretical Manual Theoretical background to the Strand7 finite element analysis system
Theoretical Manual Theoretical background to the Strand7 finite element analysis system Edition 1 January 2005 Strand7 Release 2.3 2004-2005 Strand7 Pty Limited All rights reserved Contents Preface Chapter
More informationMODELING AND FEM ANALYSIS OF DYNAMIC PROPERTIES OF THERMALLY OPTIMAL COMPOSITE MATERIALS
11th World Congress on Computational Mechanics (WCCM XI) 5th European Conference on Computational Mechanics (ECCM V) 6th European Conference on Computational Fluid Dynamics (ECFD VI) E. Oñate, J. Oliver
More informationPIEZOELECTRIC TECHNOLOGY PRIMER
PIEZOELECTRIC TECHNOLOGY PRIMER James R. Phillips Sr. Member of Technical Staff CTS Wireless Components 4800 Alameda Blvd. N.E. Albuquerque, New Mexico 87113 Piezoelectricity The piezoelectric effect is
More informationCSMA Towards optimal design of multiscale nonlinear structures and reduced-order modeling approaches. 1 Introduction.
CSMA 2017 13ème Colloque National en Calcul des Structures 15-19 Mai 2017, Presqu île de Giens (Var) Towards optimal design of multiscale nonlinear structures and reduced-order modeling approaches Liang
More informationFinite Element Analysis of Dynamic Properties of Thermally Optimal Two-phase Composite Structure
Vibrations in Physical Systems Vol.26 (2014) Finite Element Analysis of Dynamic Properties of Thermally Optimal Two-phase Composite Structure Abstract Maria NIENARTOWICZ Institute of Applied Mechanics,
More informationA new approach for stress-based topology optimization: Internal stress penalization
10 th World Congress on Structural and Multidisciplinary Optimization May 19-24, 2013, Orlando, Florida, USA A new approach for stress-based topology optimization: Internal stress penalization Alexander
More informationHeaviside projection based aggregation in stress constrained topology optimization
Heaviside projection based aggregation in stress constrained topology optimization Cunfu Wang Xiaoping Qian Department of Mechanical Engineering, University of Wisconsin-Madison, 1513 University Avenue,
More informationBasic Principles of Weak Galerkin Finite Element Methods for PDEs
Basic Principles of Weak Galerkin Finite Element Methods for PDEs Junping Wang Computational Mathematics Division of Mathematical Sciences National Science Foundation Arlington, VA 22230 Polytopal Element
More informationTopology optimization design of functionally graded bimorph-type piezoelectric actuators
IOP PUBLISHING Smart Mater. Struct. 16 (2007) 2605 2620 Topology optimization design of functionally graded bimorph-type piezoelectric actuators SMART MATERIALS AND STRUCTURES doi:10.1088/0964-1726/16/6/065
More informationModal Analysis: What it is and is not Gerrit Visser
Modal Analysis: What it is and is not Gerrit Visser What is a Modal Analysis? What answers do we get out of it? How is it useful? What does it not tell us? In this article, we ll discuss where a modal
More informationTopology Optimization for the Microstructure Design of Plasmonic Composites
10 th World Congress on Structural and Multidisciplinary Optimization May 19-24, 2013, Orlando, Florida, USA Topology Optimization for the Microstructure Design of Plasmonic Composites Masaki Otomori 1,
More informationChapter 2 Finite Element Formulations
Chapter 2 Finite Element Formulations The governing equations for problems solved by the finite element method are typically formulated by partial differential equations in their original form. These are
More informationJEPPIAAR ENGINEERING COLLEGE
JEPPIAAR ENGINEERING COLLEGE Jeppiaar Nagar, Rajiv Gandhi Salai 600 119 DEPARTMENT OFMECHANICAL ENGINEERING QUESTION BANK VI SEMESTER ME6603 FINITE ELEMENT ANALYSIS Regulation 013 SUBJECT YEAR /SEM: III
More informationInverse Design (and a lightweight introduction to the Finite Element Method) Stelian Coros
Inverse Design (and a lightweight introduction to the Finite Element Method) Stelian Coros Computational Design Forward design: direct manipulation of design parameters Level of abstraction Exploration
More informationThickness Optimization of a Piezoelectric Converter for Energy Harvesting
Excerpt from the Proceedings of the COMSOL Conference 29 Milan Thickness Optimization of a Piezoelectric Converter for Energy Harvesting M. Guizzetti* 1, V. Ferrari 1, D. Marioli 1 and T. Zawada 2 1 Dept.
More informationOPTIMIZED DESIGN OF SONOTRODES FOR PIEZOELECTRIC TRANSDUCERS USING TOPOLOGY OPTIMIZATION
OPTIMIZED DESIGN OF SONOTRODES FOR PIEZOELECTRIC TRANSDUCERS USING TOPOLOGY OPTIMIZATION César Yukishigue Kiyono, ckiyono@gmail.com Ronny Calixto Carbonari, ronny@usp.br Emílio Carlos Nelli Silva, ecnsilva@usp.br
More informationDesign of Functionally Graded Piezoelectric Actuators Using Topology Optimization
Design of Functionally Graded Piezoelectric Actuators Using Topology Optimization Ronny C. Carbonari a,emílio C. N. Silva a, and Glaucio H. Paulino b a Department of Mechatronics ad Mechanical Systems
More informationEigenvalue Problems in Surface Acoustic Wave Filter Simulations
Eigenvalue Problems in Surface Acoustic Wave Filter Simulations S. Zaglmayr, J. Schöberl, U. Langer FWF Start Project Y-192 3D hp-finite Elements: Fast Solvers and Adaptivity JKU + RICAM, Linz in collaboration
More informationPiezoactuator Design Considering the Optimum Placement of FGM Piezoelectric Material
Piezoactuator Design Considering the Optimum Placement of FGM Piezoelectric Material Ronny C. Carbonari a, Shinji Nishiwaki b, Glaucio H. Paulino c,andemílio C. Nelli Silva a a Department of Mechatronics
More informationTopology optimization of compliant circular-path mechanisms. based on an aggregated linear system and singular value.
INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN ENGINEERING Int. J. Numer. Meth. Engng 2000; 00:1 6 [Version: 2002/09/18 v2.02] Topology optimization of compliant circular-path mechanisms based on an aggregated
More informationPerformance Evaluation of Various Smoothed Finite Element Methods with Tetrahedral Elements in Large Deformation Dynamic Analysis
Performance Evaluation of Various Smoothed Finite Element Methods with Tetrahedral Elements in Large Deformation Dynamic Analysis Ryoya IIDA, Yuki ONISHI, Kenji AMAYA Tokyo Institute of Technology, Japan
More informationME FINITE ELEMENT ANALYSIS FORMULAS
ME 2353 - FINITE ELEMENT ANALYSIS FORMULAS UNIT I FINITE ELEMENT FORMULATION OF BOUNDARY VALUE PROBLEMS 01. Global Equation for Force Vector, {F} = [K] {u} {F} = Global Force Vector [K] = Global Stiffness
More informationMulti-actuated Functionally Graded Piezoelectric Micro-Tools Design Using Topology Optimization
Multi-actuated Functionally Graded Piezoelectric Micro-Tools Design Using Topology Optimization Ronny C. Carbonari a,emílio C. N. Silva a, and Glaucio H. Paulino b a Department of Mechatronics ad Mechanical
More informationSoft Bodies. Good approximation for hard ones. approximation breaks when objects break, or deform. Generalization: soft (deformable) bodies
Soft-Body Physics Soft Bodies Realistic objects are not purely rigid. Good approximation for hard ones. approximation breaks when objects break, or deform. Generalization: soft (deformable) bodies Deformed
More informationContents as of 12/8/2017. Preface. 1. Overview...1
Contents as of 12/8/2017 Preface 1. Overview...1 1.1 Introduction...1 1.2 Finite element data...1 1.3 Matrix notation...3 1.4 Matrix partitions...8 1.5 Special finite element matrix notations...9 1.6 Finite
More informationNUMERICAL EVALUATION OF A TEFLON BASED PIEZOELECTRIC SENSOR EFFECTIVITY FOR THE MONITORING OF EARLY AGE COCRETE STRENGTHING
NUMERICAL EVALUATION OF A TEFLON BASED PIEZOELECTRIC SENSOR EFFECTIVITY FOR THE MONITORING OF EARLY AGE COCRETE STRENGTHING Evangelos V. Liarakos Postdoctoral researcher School of Architecture, Technical
More informationLaminated Composite Plates and Shells
Jianqiao Ye Laminated Composite Plates and Shells 3D Modelling With 62 Figures Springer Table of Contents 1. Introduction to Composite Materials 1 1.1 Introduction 1 1.2 Classification of Composite Materials
More informationCode No: RT41033 R13 Set No. 1 IV B.Tech I Semester Regular Examinations, November - 2016 FINITE ELEMENT METHODS (Common to Mechanical Engineering, Aeronautical Engineering and Automobile Engineering)
More informationVIBRATION CONTROL OF RECTANGULAR CROSS-PLY FRP PLATES USING PZT MATERIALS
Journal of Engineering Science and Technology Vol. 12, No. 12 (217) 3398-3411 School of Engineering, Taylor s University VIBRATION CONTROL OF RECTANGULAR CROSS-PLY FRP PLATES USING PZT MATERIALS DILEEP
More information1 Exercise: Linear, incompressible Stokes flow with FE
Figure 1: Pressure and velocity solution for a sinking, fluid slab impinging on viscosity contrast problem. 1 Exercise: Linear, incompressible Stokes flow with FE Reading Hughes (2000), sec. 4.2-4.4 Dabrowski
More informationCheckerboard instabilities in topological shape optimization algorithms
Checkerboard instabilities in topological shape optimization algorithms Eric Bonnetier, François Jouve Abstract Checkerboards instabilities for an algorithm of topological design are studied on a simple
More informationElimination of void element influence on optimization for nonlinear compliance with a buckling constraint using moving iso-surface threshold method
th World Congress on Structural and Multidisciplinary Optimisation 7 th - th, June 5, Sydney Australia Elimination of void element influence on optimization for nonlinear compliance with a bucling constraint
More informationGeneric Strategies to Implement Material Grading in Finite Element Methods for Isotropic and Anisotropic Materials
University of Nebraska - Lincoln DigitalCommons@University of Nebraska - Lincoln Engineering Mechanics Dissertations & Theses Mechanical & Materials Engineering, Department of Winter 12-9-2011 Generic
More informationTopology Optimization of Compliant Mechanism with Geometrical Advantage
610 Topology Optimization of Compliant Mechanism with Geometrical Advantage Seungjae MIN and Younggi KIM A compliant mechanism is a mechanism that produces its motion by the flexibility of some or all
More informationFINITE ELEMENT APPROXIMATION OF A TIME-DEPENDENT TOPOLOGY OPTIMIZATION PROBLEM
ECCOMAS Congress 2016 VII European Congress on Computational Methods in Applied Sciences and Engineering M. Papadrakakis, V. Papadopoulos, G. Stefanou, V. Plevris (eds.) Crete Island, Greece, 5 10 June
More information16.20 HANDOUT #2 Fall, 2002 Review of General Elasticity
6.20 HANDOUT #2 Fall, 2002 Review of General Elasticity NOTATION REVIEW (e.g., for strain) Engineering Contracted Engineering Tensor Tensor ε x = ε = ε xx = ε ε y = ε 2 = ε yy = ε 22 ε z = ε 3 = ε zz =
More informationA 3 D finite element model for sound transmission through a double plate system with isotropic elastic porous materials
Acoustics and Vibrations Group Université de Sherbrooke, QC CANADA Département génie mécanique Université de Sherbrooke Sherbrooke, QC CANADA Tel.: (819) 821-7157 Fax: (819) 821-7163 A 3 D finite element
More informationUNCONVENTIONAL FINITE ELEMENT MODELS FOR NONLINEAR ANALYSIS OF BEAMS AND PLATES
UNCONVENTIONAL FINITE ELEMENT MODELS FOR NONLINEAR ANALYSIS OF BEAMS AND PLATES A Thesis by WOORAM KIM Submitted to the Office of Graduate Studies of Texas A&M University in partial fulfillment of the
More informationMechanics PhD Preliminary Spring 2017
Mechanics PhD Preliminary Spring 2017 1. (10 points) Consider a body Ω that is assembled by gluing together two separate bodies along a flat interface. The normal vector to the interface is given by n
More information. D CR Nomenclature D 1
. D CR Nomenclature D 1 Appendix D: CR NOMENCLATURE D 2 The notation used by different investigators working in CR formulations has not coalesced, since the topic is in flux. This Appendix identifies the
More informationMeasurement of deformation. Measurement of elastic force. Constitutive law. Finite element method
Deformable Bodies Deformation x p(x) Given a rest shape x and its deformed configuration p(x), how large is the internal restoring force f(p)? To answer this question, we need a way to measure deformation
More informationBack Matter Index The McGraw Hill Companies, 2004
INDEX A Absolute viscosity, 294 Active zone, 468 Adjoint, 452 Admissible functions, 132 Air, 294 ALGOR, 12 Amplitude, 389, 391 Amplitude ratio, 396 ANSYS, 12 Applications fluid mechanics, 293 326. See
More informationContinuum Topology Optimization of Buckling-Sensitive Structures
Continuum Topology Optimization of Buckling-Sensitive Structures Salam Rahmatalla, Ph.D Student Colby C. Swan, Associate Professor Civil & Environmental Engineering Center for Computer-Aided Design The
More informationSparse monolithic compliant mechanisms using continuum structural topology optimization
INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN ENGINEERING Int. J. Numer. Meth. Engng 2005; 62:1579 1605 Published online 3 December 2004 in Wiley InterScience (www.interscience.wiley.com). DOI: 10.1002/nme.1224
More informationNonconservative Loading: Overview
35 Nonconservative Loading: Overview 35 Chapter 35: NONCONSERVATIVE LOADING: OVERVIEW TABLE OF CONTENTS Page 35. Introduction..................... 35 3 35.2 Sources...................... 35 3 35.3 Three
More informationOptimal thickness of a cylindrical shell under dynamical loading
Optimal thickness of a cylindrical shell under dynamical loading Paul Ziemann Institute of Mathematics and Computer Science, E.-M.-A. University Greifswald, Germany e-mail paul.ziemann@uni-greifswald.de
More informationBilinear Quadrilateral (Q4): CQUAD4 in GENESIS
Bilinear Quadrilateral (Q4): CQUAD4 in GENESIS The Q4 element has four nodes and eight nodal dof. The shape can be any quadrilateral; we ll concentrate on a rectangle now. The displacement field in terms
More informationGame Physics. Game and Media Technology Master Program - Utrecht University. Dr. Nicolas Pronost
Game and Media Technology Master Program - Utrecht University Dr. Nicolas Pronost Soft body physics Soft bodies In reality, objects are not purely rigid for some it is a good approximation but if you hit
More informationFinite-element modeling of the transient response of MEMS sensors
Smart Mater. Struct. 6 (1997) 53 61. Printed in the UK PII: S0964-1726(97)79231-0 Finite-element modeling of the transient response of MEMS sensors Young-Hun Lim, Vasundara V Varadan and Vijay K Varadan
More informationStructural Health Monitoring using Shaped Sensors
Structural Health Monitoring using Shaped Sensors Michael I. Friswell and Sondipon Adhikari Swansea University, UK This paper is concerned with distributed sensors to measure the response of beam and plate
More informationChapter 6 2D Elements Plate Elements
Institute of Structural Engineering Page 1 Chapter 6 2D Elements Plate Elements Method of Finite Elements I Institute of Structural Engineering Page 2 Continuum Elements Plane Stress Plane Strain Toda
More informationA HIGHER-ORDER BEAM THEORY FOR COMPOSITE BOX BEAMS
A HIGHER-ORDER BEAM THEORY FOR COMPOSITE BOX BEAMS A. Kroker, W. Becker TU Darmstadt, Department of Mechanical Engineering, Chair of Structural Mechanics Hochschulstr. 1, D-64289 Darmstadt, Germany kroker@mechanik.tu-darmstadt.de,
More informationNONLINEAR CONTINUUM FORMULATIONS CONTENTS
NONLINEAR CONTINUUM FORMULATIONS CONTENTS Introduction to nonlinear continuum mechanics Descriptions of motion Measures of stresses and strains Updated and Total Lagrangian formulations Continuum shell
More informationNon-linear and time-dependent material models in Mentat & MARC. Tutorial with Background and Exercises
Non-linear and time-dependent material models in Mentat & MARC Tutorial with Background and Exercises Eindhoven University of Technology Department of Mechanical Engineering Piet Schreurs July 7, 2009
More informationThe use of second-order information in structural topology optimization. Susana Rojas Labanda, PhD student Mathias Stolpe, Senior researcher
The use of second-order information in structural topology optimization Susana Rojas Labanda, PhD student Mathias Stolpe, Senior researcher What is Topology Optimization? Optimize the design of a structure
More informationCalculating Mechanical Transfer Functions with COMSOL Multiphysics. Yoichi Aso Department of Physics, University of Tokyo
Calculating Mechanical Transfer Functions with COMSOL Multiphysics Yoichi Aso Department of Physics, University of Tokyo Objective Suspension Point You have a pendulum like the one shown on the right.
More information1106. Numerical investigation of dynamical properties of vibroactive pad during hot imprint process
1106. Numerical investigation of dynamical properties of vibroactive pad during hot imprint process B. Narijauskaitė 1, A. Palevičius 2, G. Janušas 3, R. Šakalys 4 International Studies Centre, Kaunas
More informationChapter 5 Structural Elements: The truss & beam elements
Institute of Structural Engineering Page 1 Chapter 5 Structural Elements: The truss & beam elements Institute of Structural Engineering Page 2 Chapter Goals Learn how to formulate the Finite Element Equations
More informationGeneral elastic beam with an elastic foundation
General elastic beam with an elastic foundation Figure 1 shows a beam-column on an elastic foundation. The beam is connected to a continuous series of foundation springs. The other end of the foundation
More informationDevelopment of discontinuous Galerkin method for linear strain gradient elasticity
Development of discontinuous Galerkin method for linear strain gradient elasticity R Bala Chandran Computation for Design and Optimizaton Massachusetts Institute of Technology Cambridge, MA L. Noels* Aerospace
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 informationStructural Dynamics Lecture 4. Outline of Lecture 4. Multi-Degree-of-Freedom Systems. Formulation of Equations of Motions. Undamped Eigenvibrations.
Outline of Multi-Degree-of-Freedom Systems Formulation of Equations of Motions. Newton s 2 nd Law Applied to Free Masses. D Alembert s Principle. Basic Equations of Motion for Forced Vibrations of Linear
More informationMULTIMATERIAL TOPOLOGY OPTIMIZATION BY VOLUME CONSTRAINED ALLEN-CAHN SYSTEM AND REGULARIZED PROJECTED STEEPEST DESCENT METHOD
MULTIMATERIAL TOPOLOGY OPTIMIZATION BY VOLUME CONSTRAINED ALLEN-CAHN SYSTEM AND REGULARIZED PROJECTED STEEPEST DESCENT METHOD R. TAVAKOLI Abstract. A new computational algorithm is introduced in the present
More informationPost Graduate Diploma in Mechanical Engineering Computational mechanics using finite element method
9210-220 Post Graduate Diploma in Mechanical Engineering Computational mechanics using finite element method You should have the following for this examination one answer book scientific calculator No
More informationLECTURE # 0 BASIC NOTATIONS AND CONCEPTS IN THE THEORY OF PARTIAL DIFFERENTIAL EQUATIONS (PDES)
LECTURE # 0 BASIC NOTATIONS AND CONCEPTS IN THE THEORY OF PARTIAL DIFFERENTIAL EQUATIONS (PDES) RAYTCHO LAZAROV 1 Notations and Basic Functional Spaces Scalar function in R d, d 1 will be denoted by u,
More informationTunnel Reinforcement Optimization for Nonlinear Material
November 25-27, 2012, Gold Coast, Australia www.iccm-2012.org Tunnel Reinforcement Optimization for Nonlinear Material T. Nguyen* 1,2, K. Ghabraie 1,2, T. Tran-Cong 1,2 1 Computational Engineering and
More informationComparison of Models for Finite Plasticity
Comparison of Models for Finite Plasticity A numerical study Patrizio Neff and Christian Wieners California Institute of Technology (Universität Darmstadt) Universität Augsburg (Universität Heidelberg)
More informationME 475 Modal Analysis of a Tapered Beam
ME 475 Modal Analysis of a Tapered Beam Objectives: 1. To find the natural frequencies and mode shapes of a tapered beam using FEA.. To compare the FE solution to analytical solutions of the vibratory
More informationNumerical analyses of cement-based piezoelectric smart composites
Numerical analyses of cement-based piezoelectric smart composites *Jan Sladek 1, Pavol Novak 2, Peter L. Bishay 3, and Vladimir Sladek 1 1 Institute of Construction and Architecture, Slovak Academy of
More informationBake, shake or break - and other applications for the FEM. 5: Do real-life experimentation using your FEM code
Bake, shake or break - and other applications for the FEM Programming project in TMA4220 - part 2 by Kjetil André Johannessen TMA4220 - Numerical solution of partial differential equations using the finite
More informationFREE VIBRATION OF A THERMO-PIEZOELECTRIC PLATE
Inter national Journal of Pure and Applied Mathematics Volume 113 No. 11 217, 217 225 ISSN: 1311-88 (printed version); ISSN: 1314-3395 (on-line version) url: http://www.ijpam.eu ijpam.eu FREE VIBRATION
More informationFundamentals of Linear Elasticity
Fundamentals of Linear Elasticity Introductory Course on Multiphysics Modelling TOMASZ G. ZIELIŃSKI bluebox.ippt.pan.pl/ tzielins/ Institute of Fundamental Technological Research of the Polish Academy
More informationAnalytical Mechanics: Elastic Deformation
Analytical Mechanics: Elastic Deformation Shinichi Hirai Dept. Robotics, Ritsumeikan Univ. Shinichi Hirai (Dept. Robotics, Ritsumeikan Univ.) Analytical Mechanics: Elastic Deformation 1 / 60 Agenda Agenda
More informationEshan V. Dave, Secretary of M&FGM2006 (Hawaii) Research Assistant and Ph.D. Candidate. Glaucio H. Paulino, Chairman of M&FGM2006 (Hawaii)
Asphalt Pavement Aging and Temperature Dependent Properties through a Functionally Graded Viscoelastic Model I: Development, Implementation and Verification Eshan V. Dave, Secretary of M&FGM2006 (Hawaii)
More informationCoupled electrostatic-elastic analysis for topology optimization using material interpolation
Institute of Physics Publishing Journal of Physics: Conference Series 34 (2006) 264 270 doi:0.088/742-6596/34//044 International MEMS Conference 2006 Coupled electrostatic-elastic analysis for topology
More informationSensors & Transducers 2016 by IFSA Publishing, S. L.
Sensors & Transducers, Vol. 96, Issue, January 206, pp. 52-56 Sensors & Transducers 206 by IFSA Publishing, S. L. http://www.sensorsportal.com Collapse Mode Characteristics of Parallel Plate Ultrasonic
More informationRobust topology optimization of slender structures with geometric imperfections
10 th World Congress on Structural and Multidisciplinary Optimization May 19-24, 2013, Orlando, Florida, USA Robust topology optimization of slender structures with geometric imperfections Miche Jansen
More informationControl of Earthquake Induced Vibrations in Asymmetric Buildings Using Passive Damping
Control of Earthquake Induced Vibrations in Asymmetric Buildings Using Passive Damping Rakesh K. Goel, California Polytechnic State University, San Luis Obispo Abstract This paper summarizes the results
More informationVIBRATION ENERGY FLOW IN WELDED CONNECTION OF PLATES. 1. Introduction
ARCHIVES OF ACOUSTICS 31, 4 (Supplement), 53 58 (2006) VIBRATION ENERGY FLOW IN WELDED CONNECTION OF PLATES J. CIEŚLIK, W. BOCHNIAK AGH University of Science and Technology Department of Robotics and Mechatronics
More informationTopology optimization of passive constrained layer damping with partial coverage on plate
Shock and Vibration 20 (2013) 199 211 199 DOI 10.3233/SAV-2012-00738 IOS Press Topology optimization of passive constrained layer damping with partial coverage on plate Weiguang Zheng a,b, Yingfeng Lei
More informationShape Optimization for Brake Squeal
10 th World Congress on Structural and Multidisciplinary Optimization May 19-24, 2013, Orlando, Florida, USA Shape Optimization for Brake Squeal Kohei Shintani 1 and Hideyuki Azegami 1 1 Nagoya University,
More informationThe Finite Element Method for the Analysis of Non-Linear and Dynamic Systems: Thermomechanics
The Finite Element Method for the Analysis of Non-Linear and Dynamic Systems: Thermomechanics Prof. Dr. Eleni Chatzi Dr. Giuseppe Abbiati, Dr. Konstantinos Agathos Lecture 13-14 December, 2017 1 / 30 Forewords
More informationLASER GENERATED THERMOELASTIC WAVES IN AN ANISOTROPIC INFINITE PLATE
LASER GENERATED THERMOELASTIC WAVES IN AN ANISOTROPIC INFINITE PLATE H. M. Al-Qahtani and S. K. Datta University of Colorado Boulder CO 839-7 ABSTRACT. An analysis of the propagation of thermoelastic waves
More informationCHAPTER 5 SIMULATION OF A PAYLOAD FAIRING
CHAPTER 5 SIMULATION OF A PAYLOAD FAIRING In the preceding chapters, a model of a PZT actuator exciting a SS cylinder has been presented. The structural model is based on a modal expansion formulation
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 informationReview of Strain Energy Methods and Introduction to Stiffness Matrix Methods of Structural Analysis
uke University epartment of Civil and Environmental Engineering CEE 42L. Matrix Structural Analysis Henri P. Gavin Fall, 22 Review of Strain Energy Methods and Introduction to Stiffness Matrix Methods
More informationLecture 12: Finite Elements
Materials Science & Metallurgy Part III Course M6 Computation of Phase Diagrams H. K. D. H. Bhadeshia Lecture 2: Finite Elements In finite element analysis, functions of continuous quantities such as temperature
More informationMechanics of Structures (CE130N) Lab 6. 2 Principle of stationary potential energy
UNIVERSITY OF CALIFORNIA BERKELEY Department of Civil Engineering Spring 0 Structural Engineering, Mechanics and Materials Professor: S. Govindjee Mechanics of Structures (CE30N) Lab 6 Objective The objective
More informationMAE 323: Chapter 6. Structural Models
Common element types for structural analyis: oplane stress/strain, Axisymmetric obeam, truss,spring oplate/shell elements o3d solid ospecial: Usually used for contact or other constraints What you need
More informationPART I. Basic Concepts
PART I. Basic Concepts. Introduction. Basic Terminology of Structural Vibration.. Common Vibration Sources.. Forms of Vibration.3 Structural Vibration . Basic Terminology of Structural Vibration The term
More information