EXPERIMENT 4: AN ELECTRICAL-THERMAL ACTUATOR

Similar documents
Module 2: Thermal Stresses in a 1D Beam Fixed at Both Ends

Module 10: Free Vibration of an Undampened 1D Cantilever Beam

Transient Thermal Analysis of a Fin

Due Monday, September 14 th, 12:00 midnight

MATERIAL MECHANICS, SE2126 COMPUTER LAB 2 PLASTICITY

Workshop 8. Lateral Buckling

Leaf Spring (Material, Contact, geometric nonlinearity)

MATERIAL MECHANICS, SE2126 COMPUTER LAB 3 VISCOELASTICITY. k a. N t

Linear Static Analysis of a Cantilever Beam (SI Units)

Creating Axisymmetric Models in FEMAP

Benchmarking COMSOL Multiphysics 3.4

Composite FEM Lab-work

Linear Static Analysis of a Cantilever Beam (CBAR Problem)

Thermal Stress Analysis of a Bi- Metallic Plate

Cable Tension APPENDIX D. Objectives: Demonstrate the use of elastic-plastic material properties. Create an enforced displacement on the model.

Coupled Field Analysis using the ANSYS/Multiphysics Commercial FEA Code

Lab 1: Numerical Solution of Laplace s Equation

CHAPTER 4 DESIGN AND ANALYSIS OF CANTILEVER BEAM ELECTROSTATIC ACTUATORS

Module I: Two-dimensional linear elasticity. application notes and tutorial. Problems

Finite Element Analysis of Piezoelectric Cantilever

DESIGN AND SIMULATION OF ACCELEROMETER SPRINGS

MATERIAL MECHANICS, SE2126 COMPUTER LAB 4 MICRO MECHANICS. E E v E E E E E v E E + + = m f f. f f

Response Spectrum Analysis Shock and Seismic. FEMAP & NX Nastran

Non-linear and time-dependent material models in Mentat & MARC. Tutorial with Background and Exercises

Plane and axisymmetric models in Mentat & MARC. Tutorial with some Background

Geometric Nonlinear Analysis of a Cantilever Beam

Chapter 5. Vibration Analysis. Workbench - Mechanical Introduction ANSYS, Inc. Proprietary 2009 ANSYS, Inc. All rights reserved.

D && 9.0 DYNAMIC ANALYSIS

Settlement and Bearing Capacity of a Strip Footing. Nonlinear Analyses

Project. First Saved Monday, June 27, 2011 Last Saved Wednesday, June 29, 2011 Product Version 13.0 Release

Coupling Physics. Tomasz Stelmach Senior Application Engineer

Depletion Analysis of an Oilfield

Institute of Structural Engineering Page 1. Method of Finite Elements I. Chapter 2. The Direct Stiffness Method. Method of Finite Elements I

Elastic Properties of Solids Exercises I II for one weight Exercises III and IV give the second weight

WSEAS TRANSACTIONS on APPLIED and THEORETICAL MECHANICS. John K. Sakellaris

Multiphysics Modeling

ENSC387: Introduction to Electromechanical Sensors and Actuators LAB 3: USING STRAIN GAUGES TO FIND POISSON S RATIO AND YOUNG S MODULUS

Axisymmetric Modeling. This tutorial gives an overview of axisymmetric modeling. Learn how to:

Project Engineer: Wesley Kinkler Project Number: 4.14 Submission Date: 11/15/2003. TAMUK Truss Company Trusses Made Simple

Analysis of a Casted Control Surface using Bi-Linear Kinematic Hardening

2D Liquefaction Analysis for Bridge Abutment

EE C247B / ME C218 INTRODUCTION TO MEMS DESIGN SPRING 2014 C. Nguyen PROBLEM SET #4

Procedure for Performing Stress Analysis by Means of Finite Element Method (FEM)

Example 37 - Analytical Beam

Magneto-Mechanical Modeling and Simulation of MEMS Sensors Based on Electroactive Polymers

Force vs time. IMPULSE AND MOMENTUM Pre Lab Exercise: Turn in with your lab report

Quintic beam closed form matrices (revised 2/21, 2/23/12) General elastic beam with an elastic foundation

Activity P27: Speed of Sound in Air (Sound Sensor)

CHAPTER 5 FIXED GUIDED BEAM ANALYSIS

ELEC 1908 The Electric Potential (V) March 28, 2013

ANSYS problem. Situation:

Mechanical Properties Rev 5.0

Modeling a Composite Slot Cross-Section for Torsional Analysis

Obtaining the Center of Gravity, Volume, Weight and More

Effects of Chrome Pattern Characteristics on Image Placement due to the Thermomechanical Distortion of Optical Reticles During Exposure

Material Property Definition

Tutorial Number 18: Heat transfer analysis of a teapot

ANSYS Coupled-Field Analysis Guide

Due Monday, October 19 th, 12:00 midnight. Problem 1 2D heat conduction in a rectangular domain (hand calculation) . A constant heat source is given:

CHEMDRAW ULTRA ITEC107 - Introduction to Computing for Pharmacy. ITEC107 - Introduction to Computing for Pharmacy 1

Software Verification

Contents. Feature Articles. On the Web. Resources. A Publication for ANSYS Users

Software BioScout-Calibrator June 2013

Analysis of Planar Truss

An example of panel solution in the elastic-plastic regime

ME 475 Modal Analysis of a Tapered Beam

Application Note. Capacitive Pressure sensor design

BUILDING BASICS WITH HYPERCHEM LITE

ENGN 2340 Final Project Report. Optimization of Mechanical Isotropy of Soft Network Material

APPENDIX 1 MATLAB AND ANSYS PROGRAMS

Finite Element Modelling with Plastic Hinges

13 Dewatered Construction of a Braced Excavation

Thermo-Mechanical Analysis of a Multi-Layer MEMS Membrane

MET 487 Instrumentation and Automatic Controls. Lecture 13 Sensors

Stress analysis of a stepped bar

St art. rp m. Km /h 1: : : : : : : : : : : : :5 2.5.

Design and Simulation of Micro-cantilever

Steps in the Finite Element Method. Chung Hua University Department of Mechanical Engineering Dr. Ching I Chen

Two problems to be solved. Example Use of SITATION. Here is the main menu. First step. Now. To load the data.

Institute of Structural Engineering Page 1. Method of Finite Elements I. Chapter 2. The Direct Stiffness Method. Method of Finite Elements I

Modelling and numerical simulation of the wrinkling evolution for thermo-mechanical loading cases

UNLOADING OF AN ELASTIC-PLASTIC LOADED SPHERICAL CONTACT

General elastic beam with an elastic foundation

Lab 1: Dynamic Simulation Using Simulink and Matlab

NMR Predictor. Introduction

Advanced Simulation of Sealings CADFEM GmbH Rainer Rauch

1106. Numerical investigation of dynamical properties of vibroactive pad during hot imprint process

Design and Analysis of Various Microcantilever Shapes for MEMS Based Sensing

Hot Spot / Point Density Analysis: Kernel Smoothing

Femtet What s New. Murata Software Co., Ltd. All Rights Reserved, Copyright c Murata Manufacturing Co., Ltd.

TRANSVERSE STRESSES IN SHEAR LAG OF BOX-GIRDER BRIDGES. Wang Yuan

Chapter 10 AC Analysis Using Phasors

Finite Element Analysis Prof. Dr. B. N. Rao Department of Civil Engineering Indian Institute of Technology, Madras. Module - 01 Lecture - 11

ISSN: [HARIHARAN * et al., 7(3): March, 2018] Impact Factor: 5.164

SECOORA Data Portal Exercises

DESIGN AND SIMULATION OF UNDER WATER ACOUSTIC MEMS SENSOR

OPTIMIZATION OF THE COMPOSITE MATERIALS OF TANKS USING FINITE ELEMENT METHOD AND STRAIN GAUGES

SIMULATION AND OPTIMIZATION OF MEMS PIEZOELECTRIC ENERGY HARVESTER WITH A NON-TRADITIONAL GEOMETRY

Load Cell Design Using COMSOL Multiphysics

Example Resistive Heating

Transcription:

EXPERIMENT 4: AN ELECTRICAL-THERMAL ACTUATOR 1. OBJECTIVE: 1.1 To analyze an electrical-thermal actuator used in a micro-electromechanical system (MEMS). 2. INTRODUCTION 2.1 Introduction to Thermal Actuator This experiment demonstrates how to analyze an electrical-thermal actuator used in a micro-electromechanical system (MEMS). The thermal actuator is fabricated from polysilicon and is shown below. Figure 1: The thermal actuator used in MEMS The thermal actuator works on the basis of a differential thermal expansion between the thin arm and blade. The required analysis is a coupled-field multiphysics analysis that accounts for the interaction (coupling) between thermal, electric, and structural fields. A potential difference applied across the electrical connection pads induces a current to flow through the arm and blade. The current flow and the resistivity of the polysilicon produce Joule heating (I 2 R) in the arm blade. The Joule heating causes the arm and the blade to heat up. Temperatures in the range of 700-1300 o K are generated. These temperatures produce thermal strain and thermally induced deflections. The current in the thin arm is greater than the current in the blade. Therefore, the thin arm heats up more than the blade, which causes the actuator to bend towards the blade. The maximum deformation occurs at the actuator tip. The amount of tip deflection (or force applied if the tip is restrained) is a direct function of the applied potential difference. 1

Therefore, the amount of tip deflection (or applied force) can be accurately calibrated as a function of applied voltage. These thermal actuators are used to move micro devices, such as ratchets and gear trains. Arrays of thermal actuators can be connected together at their blade tips to multiply the effective force. The main objective of the analysis is to compute the blade tip deflection for an applied potential difference across the electrical connection pads. Additional objectives are to: Obtain temperature, voltage, and displacement plots Animate displacement results Determine total current and heat flow. To define material properties for this analysis, you must convert the given units for Young's modulus, resistivity, and thermal conductivity to µmksv units. The units have been converted to µmksv for you, and are shown in the following table. Material Properties for Polysilicon (µmksv units) Young's modulus Poisson's ratio 0.22 Resistivity Coefficient of thermal expansion Thermal conductivity 2.2 Import Geometry Step 1: Import IGES file. 169e3 MPa 2.3e-11 ohm-µm 2.9e-6/ o K 150e6 pw/µm o K You will begin by importing the model of the actuator. It is provided for you in the form of an IGES file. 1. Utility Menu> File> Import> IGES 2. No defeaturing 3. [OK] 4. File name: actuator.iges \Program Files\Ansys Inc\V81\ANSYS\data\models\actuator.iges 5. [OK] 2

Figure 2: Imported.iges file 2.3 Define Materials Step 2: Define element type. You will now specify the element type as coupled-field element SOLID98, and the default degrees of freedom [KEYOPT(1) = 0]: UX, UY, UZ, TEMP, VOLT, MAG. This problem makes use of all but the MAG degree of freedom. 1. Main Menu> Preprocessor> Element Type> Add/Edit/Delete 2. [Add...] 3. Coupled Field (left column) 4. Scalar Tet 98 (right column) 5. [OK] 6. [Close] 2.4 Define Material Properties Step 3: Define material properties. You will now enter the material property values for polysilicon. These values are for Young's modulus, Poisson's ratio, thermal expansion coefficient, thermal conductivity, and resistivity. The units for these values are µmksv. 1. Main Menu> Preprocessor> Material Props> Material Models 2. (double-click) Structural, then Linear, then Elastic, then Isotropic 3. EX = 169e3 4. PRXY = 0.22 5. [OK] 6. (double-click) Thermal Expansion Coef, then Isotropic 7. ALPX = 2.9e-6 8. [OK] 9. (double-click) Thermal, then Conductivity, then Isotropic 10. KXX = 150e6 11. [OK] 12. (double-click) Electromagnetics, then Resistivity, then Constant 13. RSVX = 2.3e-11 14. [OK] 15. Material> Exit 3

2.5 Generate Mesh Step 4: Mesh the model. You will now obtain a course volume mesh using tetrahedral shapes. Setting the SmartSizing control to 10 (maximum) achieves the course mesh. 1. Main Menu> Preprocessor> Meshing> MeshTool 2. (check) SmartSize 3. (slide) Course = 10 4. [Mesh] 5. [Pick All] 6. [Close] MeshTool 2.6 Apply Loads Step 5: Plot areas. Figure 3: Meshed actuator You will now plot individual areas in preparation for when you apply the boundary conditions to the electrical connection pads, which you will do in the next two steps. You will first set ANSYS to display each of the areas in a distinguishing color and number. 1. Utility Menu> PlotCtrls> Numbering 2. (check) Area numbers to On 3. [OK] 4. Utility Menu> Plot> Areas Figure 4: Area numbering 4

Step 6: Apply boundary conditions to electrical connection pad 1. You will now apply the voltage, temperature, and displacement boundary conditions to electrical connection pad 1. You will first select this pad area so you only have to perform the picking function once there, even though you are applying three different types of boundary conditions. 1. Utility Menu> Select> Entities 2. (first drop down) Areas 3. (second drop down) By Num/Pick 4. [OK] 5. Pick electrical connection pad 1 (the upper pad). Ensure that you have picked the correct area by holding the mouse button down and dragging the mouse until ONLY the pad area highlights, then release the button. Figure 5: Selecting the electrical connection pad 6. [OK] Now apply the voltage boundary condition to pad 1. 7. Main Menu> Preprocessor> Loads> Define Loads> Apply> Electric> Boundary> Voltage> On Areas 8. [Pick All] By choosing Pick All, you pick only the area representing pad 1 because that is the only entity you currently have selected. 9. Load VOLT Value = 5 10. [OK] Notice the voltage boundary condition symbols added to pad 1. Now apply the temperature boundary condition to pad 1. Figure 6: Voltage boundary condition is applied to pad 1 11. Main Menu> Preprocessor> Loads> Define Loads> Apply> Thermal> Temperature> On Areas 12. [Pick All] 13. Load TEMP value = 30 14. [OK] Notice the temperature boundary condition symbols added to pad 1. Now apply the displacement boundary conditions to pad 1. 5

Figure 7: Temperature boundary condition applied to pad 1 15. Main Menu> Preprocessor> Loads> Define Loads> Apply> Structural> Displacement> On Areas 16. [Pick All] 17. DOF to be constrained = UX 18. Displacement value = 0 19. [Apply] 20. [Pick All] 21. DOF to be constrained = UY 22. [Apply] 23. [Pick All] 24. DOF to be constrained = UZ 25. [OK] Notice that ANSYS cumulatively added the displacement boundary condition symbols to the pad after you applied them (i.e., X constraint, then Y constraint, then Z constraint). Figure 8: The addition of displacement boundary condition symbols to the pad 1 Step 7: Apply boundary conditions to electrical connection pad 2. You will now apply the voltage, temperature, and displacement boundary conditions to electrical connection pad 2. The procedure you follow is identical to the one you just performed to add boundary conditions to pad 1. 1. Utility Menu> Select> Entities 2. (first drop down) Areas 3. (second drop down) By Num/Pick 4. [OK] 5. Pick electrical connection pad 2 (the lower pad). Ensure that you have picked the correct area by holding the mouse button down and dragging the mouse until ONLY the pad area highlights, then release the button. 6

Figure 9: Picking pad 2 for the electrical connection 6. [OK] Now apply the voltage boundary condition to pad 2. 7. Main Menu> Preprocessor> Loads> Define Loads> Apply> Electric> Boundary> Voltage> On Areas 8. [Pick All] By choosing Pick All, you pick only the area representing pad 2 because that is the only entity you currently have selected. 9. Load VOLT Value = 0 10. [OK] Notice the voltage boundary condition symbols added to pad 2. Now apply the temperature boundary condition to pad 2. Figure 10: Temperature boundary condition is applied to pad 2 11. Main Menu> Preprocessor> Loads> Define Loads> Apply> Thermal> Temperature> On Areas 12. [Pick All] 13. Load TEMP Value = 30 14. [OK] Notice the temperature boundary condition symbols added to pad 2. Now apply the displacement boundary conditions to pad 2. Figure 11: Displacement boundary condition is applied to pad 2 15. Main Menu> Preprocessor> Loads> Define Loads> Apply> Structural> Displacement> On Areas 16. [Pick All] 17. DOF to be constrained = UX 18. Displacement value = 0 19. [Apply] 20. [Pick All] 21. DOF to be constrained = UY 22. [Apply] 7

23. [Pick All] 24. DOF to be constrained = UZ 25. [OK] Notice that ANSYS cumulatively added the displacement boundary condition symbols to the pad after you applied them (i.e., X constraint, then Y constraint, then Z constraint). Before solving the problem, you must select the entire finite element model. Figure 12: The addition of displacement boundary condition symbols to pad 2 26. Utility Menu> Select> Everything 2.7 Obtain Solution Step 8: Solve. You will now initiate the ANSYS solution. 1. Main Menu> Solution> Solve> Current LS 2. Review the information in the status window, then choose: File> Close (Windows platforms), 3. [OK] 4. [Yes] to warning message. 5. [Close] the information window when the solution is done. 2.8 Review results Step 9: Plot temperature results. You will now plot the temperature results. This is one of the objectives of this analysis. 1. Main Menu> General Postproc> Read Results> Last Set 2. Main Menu> General Postproc> PlotResults> Contour Plot> Nodal Solu 3. (left column) DOF solution 4. (right column) Temperature TEMP 5. [OK] Discussion a. Interpret your results by referring to the legend beneath your plot. b. Explain on the colours of the electrical connection pad. c. Which part of the actuator is at higher temperature? 8

Step 10: Plot voltage results. You will now plot the voltage results. This is one of the objectives of this analysis. 1. Main Menu> General Postproc> Contour Plot> Plot Results> Nodal Solu 2. (left column) DOF solution 3. (right column) Elec poten VOLT 4. [OK] Discussion a. Interpret your results by referring to the legend beneath your plot. b. Why the electrical connection pads are distinctly two different colours? c. Explain the colour changing from the pads end to the blade tip end. Step 11: Plot displacement results and animate. You will now plot and animate the displacement results. These are two of the objectives of this analysis. 1. Main Menu> General Postproc> PlotResults> Contour Plot> Nodal Solu 2. (left column) DOF solution 3. (right column) Translation UY 4. [OK]. Answer the discussion question before you continue to next section. 5. Utility Menu> PlotCtrls> Animate> Deformed Results 6. (left column) DOF solution 7. (right column) Translation UY 8. [OK] 9. [Close] after making desired choices in the Animation Controller. Discussion a. Interpret your result by referring to the legend beneath your plot. Step 12: List total heat flow and current. You will now obtain a listing of results including the total heat flow and current. These are two of the objectives for this analysis. 1. Utility Menu> Plot> Areas 2. Utility Menu> Select> Entities 3. (first drop down) Areas 4. (second drop down) By Num/Pick 5. [OK] 6. Pick the electrical connection pad on the thin side of the actuator (lower pad as shown below). 9

7. [OK] 8. Utility Menu> Select> Entities 9. (first drop down) Nodes 10. (second drop down) Attached to 11. Areas, all 12. [OK] 13. Utility Menu> List> Results> Reaction Solution 14. 1 st 10 items 15. [OK] 16. When you are done viewing the listing, choose: File> Close (Windows platforms), 17. Utility Menu> Select> Everything Discussion a. List all your total heat flow and current. Step 13: Exit the ANSYS program. 1. Toolbar: Quit 2. Quit - No Save! 3. [OK] 10

4.6 Discussions Step 9: Plot temperature results. a. Interpret your results by referring to the legend beneath your plot. b. Explain on the colours of the electrical connection pad. c. Which part of the actuator is at higher temperature? 11

Step 10: Plot voltage results. a. Interpret your results by referring to the legend beneath your plot. b. Why the electrical connection pads are distinctly two different colours? c. Explain the colour changing from the pads end to the blade tip end. 12

Step 11: Plot displacement results and animate. a. Interpret your result by referring to the legend beneath your plot. Step 12: List total heat flow and current. a. List all your total heat flow and current. 13

In your own words, describe an electrical-thermal actuator. 5.0 CONCLUSION - In your own words, conclude the Experiment 4: An electrical-thermal actuator 14