Computationally efficient modelling of pattern dependencies in the micro-embossing of thermoplastic polymers

Similar documents
Modeling of MEMS Fabrication Processes

Outline: droplet-dispensed NIL simulation

Towards nanoimprint lithography-aware layout design checking. 25 February 2010 Hayden Taylor and Duane Boning Massachusetts Institute of Technology

process dependencies in nanoimprint

A continuum theory of amorphous solids undergoing large deformations, with application to polymeric glasses

Tool- and pattern-dependent spatial variations in silicon deep reactive ion etching

UNLOADING OF AN ELASTIC-PLASTIC LOADED SPHERICAL CONTACT

STANDARD SAMPLE. Reduced section " Diameter. Diameter. 2" Gauge length. Radius

STUDIES ON NANO-INDENTATION OF POLYMERIC THIN FILMS USING FINITE ELEMENT METHODS

Effect of embedment depth and stress anisotropy on expansion and contraction of cylindrical cavities

6.37 Determine the modulus of resilience for each of the following alloys:

Supplementary Information: Nanoscale heterogeneity promotes energy dissipation in bone

MECHANICS OF MATERIALS

Prediction of the bilinear stress-strain curve of engineering material by nanoindentation test

Mechanical Properties of Polymers. Scope. MSE 383, Unit 3-1. Joshua U. Otaigbe Iowa State University Materials Science & Engineering Dept.

Final Project: Indentation Simulation Mohak Patel ENGN-2340 Fall 13

Lecture 7, Foams, 3.054

A multiscale framework for lubrication analysis of bearings with textured surface

A CONTACT-MECHANICS BASED MODEL FOR DISHING AND EROSION IN

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

NANOINDENTATION STUDIES OF PAPER

Elements of Polymer Structure and Viscoelasticity. David M. Parks Mechanics and Materials II February 18, 2004

Agricultural Science 1B Principles & Processes in Agriculture. Mike Wheatland

Analysis of contact deformation between a coated flat plate and a sphere and its practical application

Back Calculation of Rock Mass Modulus using Finite Element Code (COMSOL)

Constitutive Model for High Density Polyethylene to Capture Strain Reversal

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

Supplementary Figures

Flexural properties of polymers

On Springback Prediction In Stamping Of AHSS BIW Components Utilizing Advanced Material Models

Modelling the behaviour of plastics for design under impact

Example-3. Title. Description. Cylindrical Hole in an Infinite Mohr-Coulomb Medium

Mechanics of Earthquakes and Faulting

Contact Modeling of Rough Surfaces. Robert L. Jackson Mechanical Engineering Department Auburn University

Ratcheting deformation in thin film structures

Supporting Information. Interfacial Shear Strength of Multilayer Graphene Oxide Films

A FINITE ELEMENT STUDY OF ELASTIC-PLASTIC HEMISPHERICAL CONTACT BEHAVIOR AGAINST A RIGID FLAT UNDER VARYING MODULUS OF ELASTICITY AND SPHERE RADIUS

Measurement and modelling of contact stiffness

Massachusetts Institute of Technology Department of Aeronautics and Astronautics Cambridge, MA Problem Set 9

Stress-Strain Behavior

Application of nanoindentation technique to extract properties of thin films through experimental and numerical analysis

C.J. Bennett, W. Sun Department of Mechanical, Materials and Manufacturing Engineering, University of Nottingham, Nottingham NG7 2RD, UK

MICROMECHANICAL ANALYSIS OF FRP COMPOSITES SUBJECTED TO LONGITUDINAL LOADING

CHAPTER 17 MECHANICAL PROPERTIES OF THIN POLYMER FILMS WITHIN CONTACTS

Modelling the nonlinear shear stress-strain response of glass fibrereinforced composites. Part II: Model development and finite element simulations

Supporting Information

UNIVERSITY OF SASKATCHEWAN ME MECHANICS OF MATERIALS I FINAL EXAM DECEMBER 13, 2008 Professor A. Dolovich

Optimization of blank dimensions to reduce springback in the flexforming process

FINITE ELEMENT ANALYSIS OF SLIDING CONTACT BETWEEN A CIRCULAR ASPERITY AND AN ELASTIC URFACE IN PLANE STRAIN CONDITION

Theory at a Glance (for IES, GATE, PSU)

FCP Short Course. Ductile and Brittle Fracture. Stephen D. Downing. Mechanical Science and Engineering

Nanoindentation of Fibrous Composite Microstructures: Experimentation and Finite Element Investigation. Mark Hardiman

Determining the Elastic Modulus and Hardness of an Ultrathin Film on a Substrate Using Nanoindentation

6.4 A cylindrical specimen of a titanium alloy having an elastic modulus of 107 GPa ( psi) and

ME 2570 MECHANICS OF MATERIALS

Prediction of geometric dimensions for cold forgings using the finite element method

EFFECT OF PILE-UP ON THE MECHANICAL CHARACTERISTICS OF STEEL WITH DIFFERENT STRAIN HISTORY BY DEPTH SENSING INDENTATION

Nonlinear Finite Element Modeling of Nano- Indentation Group Members: Shuaifang Zhang, Kangning Su. ME 563: Nonlinear Finite Element Analysis.

Effect of Strain Hardening on Unloading of a Deformable Sphere Loaded against a Rigid Flat A Finite Element Study

NORMAL STRESS. The simplest form of stress is normal stress/direct stress, which is the stress perpendicular to the surface on which it acts.

Exercise: concepts from chapter 8

(MPa) compute (a) The traction vector acting on an internal material plane with normal n ( e1 e

EMA 3702 Mechanics & Materials Science (Mechanics of Materials) Chapter 2 Stress & Strain - Axial Loading

The plastic behaviour of silicon subjected to micro-indentation

Presentation at IEEE, SC2 Mtg., Shanghai 16th October Marek Tengler

EE C245 ME C218 Introduction to MEMS Design Fall 2007

MECE 3321 MECHANICS OF SOLIDS CHAPTER 3

Mechanical Properties of Materials

Finite Element Method in Geotechnical Engineering

Mechanics of Irregular Honeycomb Structures

Prediction of Elastic Constants on 3D Four-directional Braided

Design of a fastener based on negative Poisson's ratio foam adapted from

Outline. 4 Mechanical Sensors Introduction General Mechanical properties Piezoresistivity Piezoresistive Sensors Capacitive sensors Applications

Characterisation of the mechanical behaviour of a polyurethane elastomer based on indentation and tensile creep experiments

Calculation of Thick-walled Shell Taking into Account Nonlinearity and Temperature Inhomogeneity of the Concrete at Subzero Temperatures

MASSACHUSETTS INSTITUTE OF TECHNOLOGY DEPARTMENT OF MATERIALS SCIENCE AND ENGINEERING CAMBRIDGE, MASSACHUSETTS 02139

Multiscale modeling of failure in ABS materials

Module-4. Mechanical Properties of Metals

International Communications in Heat and Mass Transfer 36 (2009) Contents lists available at ScienceDirect

Keysight Technologies Instrumented Indentation Testing with the Keysight Nano Indenter G200. Application Note

20. Rheology & Linear Elasticity

Ratcheting and Rolling Contact Fatigue Crack Initiation Life of Rails under Service Loading. Wenyi YAN Monash University, Australia

Mechanics of Solids. Mechanics Of Solids. Suraj kr. Ray Department of Civil Engineering

SEMM Mechanics PhD Preliminary Exam Spring Consider a two-dimensional rigid motion, whose displacement field is given by

Application of Discrete Element Method to Study Mechanical Behaviors of Ceramic Breeder Pebble Beds. Zhiyong An, Alice Ying, and Mohamed Abdou UCLA

TABLE OF CONTENTS CHAPTER TITLE PAGE TITLE PAGE DECLARATION DEDIDATION ACKNOWLEDGEMENTS ABSTRACT ABSTRAK

A Finite Element Study of Elastic-Plastic Hemispherical Contact Behavior against a Rigid Flat under Varying Modulus of Elasticity and Sphere Radius

Determining thermal noise limiting properties of thin films

Chapter Two: Mechanical Properties of materials

ME 207 Material Science I

3.032 Problem Set 4 Fall 2007 Due: Start of Lecture,

DERIVATIVE OF STRESS STRAIN, DEVIATORIC STRESS AND UNDRAINED COHESION MODELS BASED ON SOIL MODULUS OF COHESIVE SOILS

Ultrasonic Anisotropic Conductive Films (ACFs) Bonding of Flexible Substrates on Organic Rigid Boards at Room Temperature

Mechanical properties of polymers: an overview. Suryasarathi Bose Dept. of Materials Engineering, IISc, Bangalore

Frontiers of Fracture Mechanics. Adhesion and Interfacial Fracture Contact Damage

PURE BENDING. If a simply supported beam carries two point loads of 10 kn as shown in the following figure, pure bending occurs at segment BC.

Fig. 1. Circular fiber and interphase between the fiber and the matrix.

Introduction to Engineering Materials ENGR2000. Dr. Coates

L Young modulus = gradient L/A B1

HELICAL BUCKLING OF DRILL-STRINGS

Transcription:

Computationally efficient modelling of pattern dependencies in the micro-embossing of thermoplastic polymers Hayden Taylor and Duane Boning Microsystems Technology Laboratories Massachusetts Institute of Technology Ciprian Iliescu and Bangtao Chen Institute of Bioengineering and Nanotechnology, Singapore 24 September 2007

Hot micro- and nano-embossing stamp Glass-transition temperature temperature polymer load time t load t hold To choose an optimal process, we need to assign values to Heat Time Our load and temperature are constrained by Equipment Stamp and substrate properties 2

PMMA in compression N.M. Ames, Ph.D. thesis, MIT, 2007 3

PMMA in compression, 140 C using model of N.M. Ames, Ph.D. thesis, MIT, 2007 4

PMMA in compression (T g = 105 C) using model of N.M. Ames, Ph.D. thesis, MIT, 2007 5

Starting point: linear elastic material model E(T) Embossing done at high temperature, with low elastic modulus Deformation frozen in place by cooling before unloading Wish to compute deformation of a layer when embossed with an arbitrarily patterned stamp Take discretized representations of stamp and substrate 6

Response of material to unit pressure at one location General load response: 2 1 ν w( x, y) = πe p( ξ, η) 2 ( x ξ ) + ( y η) 2 dξdη w load radius, r Point load response wr = constant Response to unit pressure in a single element of the mesh: 1 ν = π E 2 F i, j 2 2 1 2 2 1 + 1, [ f ( x, y ) f ( x, y ) f ( x, y ) f ( x y )] ( ) ( ) ( ) 2 2 2 2 x, y = y ln x + x + y + x ln y + x y f + 1 F i,j defined here x 1,y 1 x 2,y 2 Unit pressure here 7

1-D verification of approach for PMMA at 130 C Iteratively find distribution of pressure consistent with stamp remaining rigid while polymer deforms Fit elastic modulus that is consistent with observed deformations Extracted Young s modulus ~ 5 MPa at 130 C 8

2-D linear elastic model succeeds with PMMA at 125 C Si stamp 1 2 cavity Simulation protrusion 1 mm 15 μm Topography (micron) 3 5 4 6 7 8 1 2 3 4 5 6 7 8 0 Lateral position (mm) Lateral position (mm) Thick, linear-elastic material model Experimental data 9

Linear-elastic model succeeds at 125 C, p ave = 0.5 MPa stamp penetration w polymer p 10

Linear-elastic model succeeds at 125 C, p ave = 1 MPa Features filled, 1MPa 11

Linear elastic model succeeds below yielding at other temperatures 12

Extracted PMMA Young s moduli from 110 to 140 C 13

Material flows under an average pressure of 8 MPa at 110 C stamp polymer 14

Yielding at 110 C stamp penetration polymer w Simple estimates of strain rate: penetration w t hold 2 10-3 to 10-1 during loading 10-4 to 10-3 during hold Local contact pressure at feature corners > 8 MPa 15

Yielding at 110 C stamp penetration polymer w Simple estimates of strain rate: penetration w t hold 2 10-3 to 10-1 during loading 10-4 to 10-3 during hold Local contact pressure at feature corners > 8 MPa N.M. Ames, Ph.D. Thesis, MIT, 2007 16

Modelling combined elastic/plastic behavior Compressive stress Yield stress 0.4 Compressive strain Plastic flow Deborah number De = t material /t load, hold De << 1 De ~ 1 De >> 1 Consider plastic deformation instantaneous Consider flow to be measurable but not to modify the pressure distribution substantially during hold 17

Modelling combined elastic/plastic behavior Elastic: E(T) Plastic flow De << 1 De ~ 1 De >> 1 Plastic flow w + e yield ( ) ( ) ( ) ( ) ( ) x,y = p x,y f x,y + p x,y p A Bt f ( x,y) hold p Existing linear-elastic component Tuned to represent cases from capillary filling to non-slip Poiseuille flow f e Material compressed f p Volume conserved radius radius 18

Summary and future directions The merits of a linear-elastic embossing polymer model have been probed This simulation approach completes an 800x800-element simulation in: ~ 45 s (without filling) ~ 4 min (with some filling) Our computational approach can be extended to capture yielding and plastic flow Is a single pressure distribution solution sufficient to model viscoelasto-plastic behaviour? Abstract further: mesh elements containing many features 19

Acknowledgements Nici Ames, Matthew Dirckx, David Hardt, Yee Cheong Lam and Lallit Anand The Singapore-MIT Alliance 20