Realtime 3D stress measurement in curing epoxy packaging

Similar documents
Determination of packaging induced 3D stress utilizing a piezocoefficient mapping device

FISH & CHIPS: Four Electrode Conductivity / Salinity Sensor on a Silicon Multi-sensor chip for Fisheries Research

The Wien Bridge Oscillator Family

The colpitts oscillator family

Diffusion of sulfuric acid in protective organic coatings

The dipole moment of a wall-charged void in a bulk dielectric

Modelling Salmonella transfer during grinding of meat

Determination of the minimum size of a statistical representative volume element from a fibre-reinforced composite based on point pattern statistics

Aalborg Universitet. Lecture 14 - Introduction to experimental work Kramer, Morten Mejlhede. Publication date: 2015

A Note on Powers in Finite Fields

Numerical modeling of the conduction and radiation heating in precision glass moulding

Matrix representation of a Neural Network

Optimizing Reliability using BECAS - an Open-Source Cross Section Analysis Tool

An Equivalent Source Method for Modelling the Lithospheric Magnetic Field Using Satellite and Airborne Magnetic Data

Generating the optimal magnetic field for magnetic refrigeration

Lidar calibration What s the problem?

Sound absorption properties of a perforated plate and membrane ceiling element Nilsson, Anders C.; Rasmussen, Birgit

Reducing Uncertainty of Near-shore wind resource Estimates (RUNE) using wind lidars and mesoscale models

Design of Robust AMB Controllers for Rotors Subjected to Varying and Uncertain Seal Forces

Radioactivity in the Risø District July-December 2013

Radioactivity in the Risø District January-June 2014

On the Einstein-Stern model of rotational heat capacities

Exposure Buildup Factors for Heavy Metal Oxide Glass: A Radiation Shield

Measurements of the angular distribution of diffuse irradiance

Strongly 2-connected orientations of graphs

Space charge fields in DC cables

Validation of Boundary Layer Winds from WRF Mesoscale Forecasts over Denmark

Optimal acid digestion for multi-element analysis of different waste matrices

Accelerating interior point methods with GPUs for smart grid systems

Kommerciel arbejde med WAsP

Feasibility of non-linear simulation for Field II using an angular spectrum approach

A Procedure for Classification of Cup-Anemometers

Reflection and absorption coefficients for use in room acoustic simulations

Dynamic Effects of Diabatization in Distillation Columns

Crack Tip Flipping: A New Phenomenon yet to be Resolved in Ductile Plate Tearing

Steady State Comparisons HAWC2 v12.2 vs HAWCStab2 v2.12

Multiple Scenario Inversion of Reflection Seismic Prestack Data

A practical implementation of microphone free-field comparison calibration according to the standard IEC

Validity of PEC Approximation for On-Body Propagation

Global Wind Atlas validation and uncertainty

A comparison study of the two-bladed partial pitch turbine during normal operation and an extreme gust conditions

Unknown input observer based scheme for detecting faults in a wind turbine converter Odgaard, Peter Fogh; Stoustrup, Jakob

A Trust-region-based Sequential Quadratic Programming Algorithm

A systematic methodology for controller tuning in wastewater treatment plants

Geometry of power flows and convex-relaxed power flows in distribution networks with high penetration of renewables

Alkali Release from Typical Danish Aggregates to Potential ASR Reactive Concrete

On numerical evaluation of two-dimensional phase integrals

The Smallest Transistor-Based Nonautonomous Chaotic Circuit

Application of mesoscale models with wind farm parametrisations in EERA-DTOC

On machine measurements of electrode wear in micro EDM milling

Higher-order spectral modelling of the diffraction force around a vertical circular cylinder

Room Allocation Optimisation at the Technical University of Denmark

Plane Wave Medical Ultrasound Imaging Using Adaptive Beamforming

Comparisons of winds from satellite SAR, WRF and SCADA to characterize coastal gradients and wind farm wake effects at Anholt wind farm

Vindkraftmeteorologi - metoder, målinger og data, kort, WAsP, WaSPEngineering,

Surface composition of ordered alloys: An off-stoichiometric effect

Potential solution for rain erosion of wind turbine blades

Berøringsløs optisk måling af gassammensætning ved UV og IR spektroskopi

Concept Testing of a Simple Floating Offshore Vertical Axis Wind Turbine

Design possibilities for impact noise insulation in lightweight floors A parameter study

Single Shooting and ESDIRK Methods for adjoint-based optimization of an oil reservoir

Simulating wind energy resources with mesoscale models: Intercomparison of stateof-the-art

Document Version Publisher s PDF, also known as Version of Record (includes final page, issue and volume numbers)

Resonant Inerter Based Absorbers for a Selected Global Mode

A beam theory fracture mechanics approach for strength analysis of beams with a hole

Published by: PIONEER RESEARCH & DEVELOPMENT GROUP(

Citation for published version (APA): Burcharth, H. F. (1988). Stress Analysis. Aalborg: Department of Civil Engineering, Aalborg University.

Ion exchange model for phase proton exchange waveguide in LiNbO3.

Reflection of sound from finite-size plane and curved surfaces

Benchmarking of optimization methods for topology optimization problems

Validation and comparison of numerical wind atlas methods: the South African example

Phase Equilibrium in Amino Acid Salt Systems for CO2 Capture

Aalborg Universitet. Land-use modelling by cellular automata Hansen, Henning Sten. Published in: NORDGI : Nordic Geographic Information

Vectorial analysis of dielectric photonic crystal VCSEL

Prediction of airfoil performance at high Reynolds numbers.

Neutronics calculations for the ITER Collective Thomson Scattering Diagnostics

A zero-free interval for chromatic polynomials of graphs with 3-leaf spanning trees

GIS Readiness Survey 2014 Schrøder, Anne Lise; Hvingel, Line Træholt; Hansen, Henning Sten; Skovdal, Jesper Høi

Geometry explains the large difference in the elastic properties of fcc and hcp crystals of hard spheres Sushko, N.; van der Schoot, P.P.A.M.

Multi (building)physics modeling

Proving in the Isabelle Proof Assistant that the Set of Real Numbers is not Countable

Wake decay constant for the infinite wind turbine array Application of asymptotic speed deficit concept to existing engineering wake model

Graphs with given diameter maximizing the spectral radius van Dam, Edwin

Impact of the interfaces for wind and wave modeling - interpretation using COAWST, SAR and point measurements

Battling Bluetongue and Schmallenberg virus Local scale behavior of transmitting vectors

Solar radiation and thermal performance of solar collectors for Denmark

Tectonic Thinking Beim, Anne; Bech-Danielsen, Claus; Bundgaard, Charlotte; Madsen, Ulrik Stylsvig

Thermodynamics of monomer partitioning in polymer latices: effect of molar volume of the monomers

2D PIM Simulation Based on COMSOL

Aalborg Universitet. Published in: Proceedings of Roomvent Publication date: 2007

ELECTROMECHANICAL RELIABILITY TESTING OF THREE-AXIAL SILICON FORCE SENSORS

COMPARISON OF THERMAL HAND MODELS AND MEASUREMENT SYSTEMS

Application of the MCNPX-McStas interface for shielding calculations and guide design at ESS

Aalborg Universitet. Seawall Overtopping Tests sea erosion study : Sohar Cornice Frigaard, Peter. Publication date: 1990

The structure of a nucleolytic ribozyme that employs a catalytic metal ion Liu, Yijin; Wilson, Timothy; Lilley, David

Analysis of Biaxially Stressed Bridge Deck Plates

De Prisco, Renato; Comunian, Michele; Grespan, Franscesco; Piscent, Andrea; Karlsson, Anders

Polydiagnostic study on a surfatron plasma at atmospheric pressure

Transcription:

Downloaded from orbit.dtu.dk on: Dec 17, 2017 Realtime 3D stress measurement in curing epoxy packaging Richter, Jacob; Hyldgård, A.; Birkelund, Karen; Hansen, Ole; Thomsen, Erik Vilain Published in: International Solid-State Sensors, Actuators and Microsystems Conference, 2007. TRANSDUCERS 2007. Link to article, DOI: 10.1109/SENSOR.2007.4300575 Publication date: 2007 Document Version Publisher's PDF, also known as Version of record Link back to DTU Orbit Citation (APA): Richter, J., Hyldgård, A., Birkelund, K., Hansen, O., & Thomsen, E. V. (2007). Realtime 3D stress measurement in curing epoxy packaging. In International Solid-State Sensors, Actuators and Microsystems Conference, 2007. TRANSDUCERS 2007. (pp. 2083-2084). IEEE. DOI: 10.1109/SENSOR.2007.4300575 General rights Copyright and moral rights for the publications made accessible in the public portal are retained by the authors and/or other copyright owners and it is a condition of accessing publications that users recognise and abide by the legal requirements associated with these rights. Users may download and print one copy of any publication from the public portal for the purpose of private study or research. You may not further distribute the material or use it for any profit-making activity or commercial gain You may freely distribute the URL identifying the publication in the public portal If you believe that this document breaches copyright please contact us providing details, and we will remove access to the work immediately and investigate your claim.

REALTIME 3D STRESS MEASUREMENT IN CURING EPOXY PACKAGING J. Richter 1, A. Hyldgård 1, K. Birkelund 1, O. Hansen 1 and E.V. Thomsen 1 1 Department of Micro and Nanotechnology, Technical University of Denmark, Kgs Lyngby DK-2800, DENMARK (Tel : +45-4525-5700; E-mail: jar@mic.dtu.dk) Abstract: This paper presents a novel method to characterize stress in microsystem packaging. A circular p-type piezoresistor is implemented on a (001) silicon chip. We use the circular stress sensor to determine the packaging induced stress in a polystyrene tube filled with epoxy. The epoxy curing process is monitored by stress measurements. From the stress measurements we conclude that the epoxy cures in 8 hours at room temperature. We find the difference in in-plane normal stresses to be σ xx -σ yy =6.7 MPa and (σ xx +σ yy -0.4σ zz )=232 MPa. Keywords: stress sensor, epoxy, packaging. 1. INTRODUCTION Packaging and protection of a MEMS device is the last step before commercialization. It is a difficult task to protect the device and not affect the performance of the device at the same time since the electrical performance can be highly influenced by packaging induced stress. This induced stress can be measured by a dedicated stress sensor. In ref. [1] a CMOS integrated stress sensor for packaged integrated circuit dies is presented, and ref. [2] shows a van der Pauw structure used for stress sensing. In this paper we present a method to determine stress in microsystem packaging. We measure the stress distribution in an epoxy package while the epoxy is curing. The packaging concept is otherwise used to protect chip and electronics in a multisensor system developed for fisheries research [3]. This method makes complex stress analysis possible in microscale packaging. We have designed a (001) p-type circular piezoresistor to determine the stress distribution in the package. The design is inspired by the work of Bartholomeyczik [4] and Steiner [5]. The device is shown in Figure 1. 2. THEORY The chip consists of one circular piezoresistor. A conceptual drawing of the device is shown in Figure 2. We direct the current, I(φ), at an angle φ Figure 1. Presentation of stress sensor. The chip is cut along the [110] direction. In the centre region of the chip a single circular p-type resistor is implemented. A current is forced through the resistor by contacts near the perimeter and 8 contacts placed in a radius of 100 µm with respect to the resistor centre measure 4 potential drops. with respect to the [110] crystal direction by using four outer contacts placed near the resistor perimeter. Thus, in the centre of the circular piezoresistor we spin the current 360 o. The inner contacts are placed in a radius of 100 µm with respect to the centre of the circular piezoresistor. The radius of the resistor is 1700 µm, thus we assume a uniform current distribution in the region of the inner contacts. The potential drop, V i =V ia -V ib, is measured across the inner contacts where i=1,2,3,4. Thus, for each current direction four potential drops are measured. To a first order approximation the relative change, V i /V i,0, is equal to a linear combination of the directional stress components. Defining two fictive voltage drops, V p and V t in Figure 2, as linear combinations of the four 2083

obtain an equation system with three independent equations and we can extract σ xy, (σ xx -σ yy ), and (σ xx +σ yy +a σ zz ), where a=2π 12 /(π 11 +π 12 ). In p-type silicon this ratio is approximately a=2π 12 /(π 11 +π 12 ) -0.4 [6]. 3. EXPERIMENTAL Figure 2. Schematic of the circular piezoresistor. The illustration is not to scale. The resistor radius is 1700 µm and the inner contacts are placed in a radius of 100 µm. The four outer contacts spin the current 360 o in the centre of the resistor. The fictive voltage drops, V p and V t, are linear combinations of all measured potential drops, V i = V ia - V ib where i=1,2,3,4 as defined in Eq. (1). The potential drop V 1 is showed as an example. potential drops we obtain a condensed linear system where the number of current directions determines the number of equations. The fictive voltage drops are defined as. (1) In terms of the piezocoefficients and stress components we find The fabrication of the stress sensor chip is performed by conventional silicon microtechnology fabrication. The doping concentration of the piezoresistor layer is approximately N A =10 18 cm -3. The electrical measurements are performed with four point measurements. This is illustrated in Figure 3. The current is spun in angular steps of 5. Thus we obtain 72 measurements of each of the four voltage drops, V i. The sensors are calibrated in a state-of-the-art four point bending fixture prior to packaging. This calibration step is included in order to determine the piezocoefficients of the sensor. The four point bending bridge applies a uniaxial and uniform stress to the centre region of the chip where the circular piezoresistor is located. With this prior calibration step we are able to extract the piezocoefficients (π 11 +π 12 ) and π 44. The four point bending fixture is illustrated in Figure 4. Next, the chip is placed in a polystyrene tube which is subsequently filled with EpoTek 302-3M two-component epoxy. Figure 5 is a photograph of the packaged chip after curing. I V, (2), (3) where V p =V p -V p,0 is the change in the fictive voltage drop V p and V p,0 is the potential drop at zero stress. The piezocoefficients are represented by π 11, π 12, and π 44 and σ xx, σ xy, and σ yy are the inplane stresses and σ zz is the out-of-plane normal stress as sketched in Figure 1. Using different values of φ in Eqs. (2)~(3) we Figure 3. The p-type piezoresistors (white) are implanted on an n-type (001) Si substrate (dark grey). A thermal oxide (light grey) is used for isolation between substrate and metal (black). The resistance of the piezoresistor is measured by a four point measurement. A constant current is forced through the resistor from the outer contacts while a high impedance voltage measurement is performed on the inner contacts. 2084

Figure 4. Illustration of state-of-the-art four point bending fixture used for calibration of the stress sensor. The chip is placed in between the two sets of blades and an actuator pushes the two parts together. The force on the chip is measured by a highly sensitive force sensor. The fixture applies a uniaxial and uniform stress to the centre region of the chip. With this calibration step we determine the piezocoefficients. Figure 5. Photograph of the packaged device after curing of the epoxy. The circular piezoresistor is placed in the centre region of the chip. A twocomponent epoxy is filled in a polystyrene tube. The packaging concept allows for direct exposure of one end of a chip while protecting the other end (along with electronics, battery, etc.) from the surroundings. We perform the same electrical measurements on the chip while the epoxy is curing as performed in the four point bending fixture. The epoxy is cured at room temperature. Figure 7 shows that the epoxy has cured after approximately 8 hours. From the calculated fictive voltage drops we are able to find the stress distribution in the package. Figure 8 shows a measurement of the fictive voltage drop Vp/Vp,0 both when a unaxial stress is applied to the device and when the device is in the package. Previous stress measurements in epoxy packaging [7] have shown a very low σxy value. Thus we assume this value to be equal to zero. From the slope and offset of the linear fit in Figure 8 we find the stress distribution in the package to be (σxx-σyy)=6.7 MPa, and (σxx+σyy0.4σzz)=232 MPa. The value for (σxx-σyy) can also be found by using the extracted values of Vt/Vp,0. For all equations σzz is coupled to π12. In p-type silicon this piezocoefficient is very small compared to the π44 coefficient [6] resulting in a large uncertainty on the estimated σzz value. This problem can be eliminated by using a n-type piezoresistor, where the magnitudes of the piezocofficients are more comparable. n-type piezoresistors are currently being developed expecting realtime 3D stress determination of the epoxy curing process with a minimum uncertainty. Additionally, in n-type piezoresistors the σzz pre- 4. RESULTS Figure 6 shows the voltage change of the four measured potential drops Vi=Vi-Vi,0 as a function of current direction, φ after the epoxy has cured. This figure also includes measurements during calibration where a uniform and uniaxial stress is applied. A linear combination of these four potential drops give us Vp/Vp,0 and Vt/Vp,0 as defined in Eqs. (2)~(3). During the curing of the epoxy we perform measurements on the stress sensor. Figure 7 shows the change in potential drop V3/V3,0 when the current density vector is directed an angle of 90ż with respect to the [110] direction. The equation for V3/V3,0 for this particular angle is. (4) The change in voltage is caused by a change in stress in the chip. This stress is a result of the contraction of epoxy during the curing process. 2085

0.2 0.1 Calibration: Slope = -82.14 10-3 In package: Slope = 7.11 10-3 V p /V p,0 0.0-0.1 Figure 6. The change in voltage drop, V i, where i=1,2,3,4, as a function of angle, φ, for both a uniaxial stress of σ xx =11 MPa (thin lines) and in package after curing (cross). The measured values of V 1, V 2, V 3, and V 4 are used in Eq. (1) to calculate V p and V t. V 3 /V 3,0 0.08 0.06 0.04 0.02 0.00-0.02 0 5 10 15 20 Time (h) Figure 7. The relative change in voltage drop, V 3 / V 3,0, as a function of time at the current angle φ=90. V 3 /V 3,0 for this angle is described by Eq. (4). factor is a=2π 12 /(π 11 +π 12 ) -2 and thus we are able to extract σ zz -σ yy and equivalently σ zz -σ xx with the use of a n-type piezoresistor. 5. CONCLUSION In this paper we determine the stress distribution during an epoxy packaging curing process, adding -0.2-1.0-0.5 0.0 0.5 1.0 Cos(2φ) Figure 8. The fictive voltage drop V p /V p,0 for both the calibration at a uniaxial stress of 77 MPa (circular symbols) and in the package (triangle). a new dimension to characterization of microsystem packaging. We expect that this method can be used in a wide range of packaging induced stress measurement applications. We can determine curing time with an uncalibrated sensor and using a calibrated sensor we can extract the stress distribution. For the presented packaging concept using an Epotek 302-3M epoxy we find the stress in the package to be (σ xx -σ yy )=6.7 MPa, and (σ xx +σ yy - 0.4σ zz )=232 MPa and a curing time of 8 hours. Using n-type piezoresistors will lower the uncertainty of the extracted stress values and enable an extraction of σ zz -σ xx and σ zz -σ yy. REFERENCES [1] P. Ruther, et al, Proc. IEEE Sensors 2006, pp. 314. [2] A. Mian, J.C. Suhling, R.C. Jaeger, IEEE Sensors Journal, 6-2 (2006), pp. 340-356. [3] A. Hyldgård, O. Hansen and E. V. Thomsen, Proc. IEEE MEMS 2005, pp. 303-306. [4] J. Bartholomeyczik, et al, Proc. IEEE MEMS 2004, pp. 263-266. [5] R. Steiner, et al, Sensors and Actuators A, 66 (1998), pp.167-172. [6] C.S. Smith, Physical Review, 94 (1954), pp. 42-49. [7] J. Richter, et al, Proc. IEEE MEMS 2007, pp. 69-72. 2086