Pulsed Magnet Crimping

Similar documents
Introduction to Laplace Transform Techniques in Circuit Analysis

FUNDAMENTALS OF POWER SYSTEMS

ELECTROMAGNETIC WAVES AND PHOTONS

PHYS 110B - HW #2 Spring 2004, Solutions by David Pace Any referenced equations are from Griffiths Problem statements are paraphrased

Tuning of High-Power Antenna Resonances by Appropriately Reactive Sources

s much time does it take for the dog to run a distance of 10.0m

Given the following circuit with unknown initial capacitor voltage v(0): X(s) Immediately, we know that the transfer function H(s) is

PHYS 110B - HW #6 Spring 2004, Solutions by David Pace Any referenced equations are from Griffiths Problem statements are paraphrased

Lecture 23 Date:

SIMON FRASER UNIVERSITY School of Engineering Science ENSC 320 Electric Circuits II. Solutions to Assignment 3 February 2005.

Lecture 15 - Current. A Puzzle... Advanced Section: Image Charge for Spheres. Image Charge for a Grounded Spherical Shell

Solving Differential Equations by the Laplace Transform and by Numerical Methods

Final Comprehensive Exam Physical Mechanics Friday December 15, Total 100 Points Time to complete the test: 120 minutes

Chapter 4. The Laplace Transform Method

Bogoliubov Transformation in Classical Mechanics

ECE Linear Circuit Analysis II

Laplace Transformation

Math 273 Solutions to Review Problems for Exam 1

two equations that govern the motion of the fluid through some medium, like a pipe. These two equations are the

EELE 3332 Electromagnetic II Chapter 10

Department of Mechanical Engineering Massachusetts Institute of Technology Modeling, Dynamics and Control III Spring 2002

Chapter 13. Root Locus Introduction

Chapter K - Problems

Bernoulli s equation may be developed as a special form of the momentum or energy equation.

into a discrete time function. Recall that the table of Laplace/z-transforms is constructed by (i) selecting to get

Social Studies 201 Notes for March 18, 2005

Magnetism Lecture Series. Electromagnetic Induction

1. /25 2. /30 3. /25 4. /20 Total /100

Phy 213: General Physics III 6/14/2007 Chapter 30 Worksheet 1

Lecture 6: Resonance II. Announcements

Figure 1 Siemens PSSE Web Site

EE C128 / ME C134 Problem Set 1 Solution (Fall 2010) Wenjie Chen and Jansen Sheng, UC Berkeley

Lecture 17: Analytic Functions and Integrals (See Chapter 14 in Boas)

ECE382/ME482 Spring 2004 Homework 4 Solution November 14,

1. Basic introduction to electromagnetic field. wave properties and particulate properties.

NCAAPMT Calculus Challenge Challenge #3 Due: October 26, 2011

Online supplementary information

III.9. THE HYSTERESIS CYCLE OF FERROELECTRIC SUBSTANCES

v 2,p = v 3,p. The total energy at P is then mv 2 p = 6.68mv 2 p 4.49Gm2 d. (3) P 2 O 3 r o Gm = v2 p d2 P 3

ECE 325 Electric Energy System Components 6- Three-Phase Induction Motors. Instructor: Kai Sun Fall 2015

Chapter 2 Sampling and Quantization. In order to investigate sampling and quantization, the difference between analog

μ + = σ = D 4 σ = D 3 σ = σ = All units in parts (a) and (b) are in V. (1) x chart: Center = μ = 0.75 UCL =

R L R L L sl C L 1 sc

Question 1 Equivalent Circuits

ME 375 EXAM #1 Tuesday February 21, 2006

7.2 INVERSE TRANSFORMS AND TRANSFORMS OF DERIVATIVES 281

Notes on Phase Space Fall 2007, Physics 233B, Hitoshi Murayama

V = 4 3 πr3. d dt V = d ( 4 dv dt. = 4 3 π d dt r3 dv π 3r2 dv. dt = 4πr 2 dr

Social Studies 201 Notes for November 14, 2003

Digital Control System

Constant Force: Projectile Motion

A novel protocol for linearization of the Poisson-Boltzmann equation

MICRO-HYDRO INSTALLATION SIZING CALCULATIONS Jacques Chaurette eng. January 17, 2008

Discover the answer to this question in this chapter.

R. W. Erickson. Department of Electrical, Computer, and Energy Engineering University of Colorado, Boulder

Physics 41 Homework Set 3 Chapter 17 Serway 7 th Edition

Physics 2. Angular Momentum. Prepared by Vince Zaccone For Campus Learning Assistance Services at UCSB

BASIC INDUCTION MOTOR CONCEPTS

AP Physics Charge Wrap up

ON THE APPROXIMATION ERROR IN HIGH DIMENSIONAL MODEL REPRESENTATION. Xiaoqun Wang

Fluid-structure coupling analysis and simulation of viscosity effect. on Coriolis mass flowmeter

AMS 212B Perturbation Methods Lecture 20 Part 1 Copyright by Hongyun Wang, UCSC. is the kinematic viscosity and ˆp = p ρ 0

FRTN10 Exercise 3. Specifications and Disturbance Models

Lecture 21. The Lovasz splitting-off lemma Topics in Combinatorial Optimization April 29th, 2004

Physics 218: Exam 1. Class of 2:20pm. February 14th, You have the full class period to complete the exam.

1.1. Curves Curves

2015 PhysicsBowl Solutions Ans Ans Ans Ans Ans B 2. C METHOD #1: METHOD #2: 3. A 4.

Moment of Inertia of an Equilateral Triangle with Pivot at one Vertex

4.6 Principal trajectories in terms of amplitude and phase function

Physics 741 Graduate Quantum Mechanics 1 Solutions to Final Exam, Fall 2014

Chapter 19. Capacitors, Resistors and Batteries

Equivalent POG block schemes

online learning Unit Workbook 4 RLC Transients

arxiv: v1 [math.mg] 25 Aug 2011

The continuous time random walk (CTRW) was introduced by Montroll and Weiss 1.

EECS2200 Electric Circuits. RLC Circuit Natural and Step Responses

A Constraint Propagation Algorithm for Determining the Stability Margin. The paper addresses the stability margin assessment for linear systems

Lecture 10 Filtering: Applied Concepts

time? How will changes in vertical drop of the course affect race time? How will changes in the distance between turns affect race time?

Linear Motion, Speed & Velocity

11.2 Stability. A gain element is an active device. One potential problem with every active circuit is its stability

Suggested Answers To Exercises. estimates variability in a sampling distribution of random means. About 68% of means fall

Problem Set 8 Solutions

The statistical properties of the primordial fluctuations

Several schematic symbols for a capacitor are shown below. The symbol resembles the two conducting surfaces separated with a dielectric.

ME 375 FINAL EXAM Wednesday, May 6, 2009

Unit I Review Worksheet Key

Preemptive scheduling on a small number of hierarchical machines

Modeling in the Frequency Domain

March 18, 2014 Academic Year 2013/14

Poynting Vector and Energy Flow W14D1

a = f s,max /m = s g. 4. We first analyze the forces on the pig of mass m. The incline angle is.

Spring 2014 EE 445S Real-Time Digital Signal Processing Laboratory. Homework #0 Solutions on Review of Signals and Systems Material

Fermi Distribution Function. n(e) T = 0 T > 0 E F

3pt3pt 3pt3pt0pt 1.5pt3pt3pt Honors Physics Impulse-Momentum Theorem. Name: Answer Key Mr. Leonard

MATEMATIK Datum: Tid: eftermiddag. A.Heintz Telefonvakt: Anders Martinsson Tel.:

New bounds for Morse clusters

At the end of this lesson, the students should be able to understand:

Section Induction motor drives

Convective Heat Transfer

Transcription:

Puled Magnet Crimping Fred Niell 4/5/00 1 Magnetic Crimping Magnetoforming i a metal fabrication technique that ha been in ue for everal decade. A large capacitor bank i ued to tore energy that i ued to bend, weld, or otherwie work a piece of metal. In the following example, we will look at a magnetoforming apparatu. The apparatu ue a capacitor and a park gap to tranfer a large amount of energy into a work coil. A thin tube of Aluminum i placed in the work coil. The work coil exert ome force on the tube, and the tube i crimped radially around a plug placed in it open end (ee diagram). Below i an electrical chematic of the apparatu, and a decription of what we are going to do to analyze thi machine. Firt, we will decribe the action of the crimper, Second, we will derive an expreion for the current a a function of time in the coil, the reonant frequency, the peak current, and the kin depth of the Aluminum tube with the emitted EM radiation, Lat, our tak remain to derive an expreion for the preure a a funtion of time, the peak preure, and a imple expreion for the peak preure baed on the phyical parameter of the apparatu. 1. Crimping effect If we conider the Al tube a a um of infiniteimal look tacked up inide a larger um of current loop (the inductor), we can ue Faraday law of induction and Lenz law to explain the crimping effect. If we create a large magnetic induction B with the olenoid over a relatively hort time period with a large I, we create a large B. Faraday how that E = B (1) So an E field i created whoe direction i uch that an emf i induced in the tube. Auming the tube i roughly the ame ize a the olenoid (all of the field goe through the tube), the induced emf in the wall of the tube will have the folowing form: 1

E = M I where the I i the changing current in the inductor, and M the mutual inductance of the ytem. Thi emf (E) give rie to the induced current. Note, however, that from E, effectively Lenz Law, the induced current flow in the oppoite direction of the original current. Thi induced current lead to a new B field, oppoing the original B field. Thee two oppoing induction exert a repulive force on each other. The force exerted by the olenoid on the tube will be much greater than the force exerted by the tube on the olenoid. Thu, the tube feel a force radially inward, giving rie to the crimping effect.. Current, Peak current, Frequency, Skin Depth, etc. The differential equation for thi circuit ha the following form: V o R b i 1 1 C t 0 () i 1+ dt L i R i = 0 (3) Where V o i the initial charging voltage of the capacitor, i 1 i the current in loop 1, and i i the current through loop. We have ued Kirchoff law and Ohm law implicitly. At the peak value, the current through R b i only 0.5mA, and the current through R i imilarly derived to be everal order of magnitude higher, o we will neglect the effect of R b in the differential equation to maintain anity. Now, we have With the caviat V o R 1 C t 0 idt L i = 0 (4) L i + R i + i C = 0 (5) The general olution i found to be i(t) = e Rt V o ω 0 L inh R t 1 4L RC (6) But, we have 4L > 1, o our general olution reduce to R C i(t) = e Rt V o ω o L in ω ot (7) Note that thi ha the form we would expect- a ine function with an exponential decay. Alo notice that ω o i the reonant frequency of the ytem. The reonant frequency ha the following form:

ω o = 1 LC R 4L (8) With the familiar relation f = ωo. Plugging in the number given, we find the frequency π of ocillation to be ω o =181811.8 rad =8.94 khz. For the i max we mut take the derivative of i(t) and et it to zero- giving u the time when the current reache a maximum. Thu, we have i(t) = 0 = R e Rt V 0 ω 0 L in ω 0t + e Rt V 0 L co ω 0t (9) Thi give u a time of 8.1 10 6 econd, which correpond to a maximum current of 1.6 10 4 Amp, or 16kA. From the damping factor of a tranmitted EM wave in a diperive medium, i.e. e βz with z the direction of propagation, we find that the ditance the wave travel before falling to 1 e of it original value i given by z = β 1, which i generally referred to a the kin depth, δ. δ = c πµωo σ (10) Above, c i the peed of light, µ i permiability of the material (µµ 0, which in our cae i imply the permiability of free pace, and σ i the conductivity of the medium, or 1 ρ. Plugging in the given number, we get δ =4.97x10 4 meter, or 0.5 mm. 3. Preure a a function of time The magnetic preure can be derived in a number of way. I will outline two, a the anwer wa upriing to me (independent of the phyical dimenion). Firt we will look at virtual work method. I will work in MKS unit: B U magnetic dτ = B (volume) = B πr l (11) µ 0 µ 0 µ 0 with l, and r being the phyical dimention of the urface we want to calculate the energy inide of. Let u call the ˆr direction the radial direction. We know that force i defined a F U. Therefore, we can define a magnetic force in the following manner: F magnetic U magnetic = U magnetic ˆr = B πrl (1) r µ 0 Notice that the force i in the ˆr direction. Thi mean that the force i directed radially inward. So, the force will be puhing the tube in an inward, cruhing manner. We can define a preure 3

which i an intereting reult. f magnetic F magnetic area = F magnetic πrl Not believing it, I tarted from the tre tenor T ij. = B µ 0 (13) T ij ɛ(e i E j 1 E ij ) + 1 µ 0 (B i B j 1 B ij ) (14) For our ytem i = j and we are only intereted in the B field, o T ij reduce to: T i=j = B µ 0 (15) Force i the time derivitive of momentum of the field, and i given a p i=j Applying the divergence theorem = j v x j T i=j d 3 x (16) Thi give u p i=j = x j T i=j ˆnda (17) F magnetic = T i=j ˆnda = B µ 0 (area) (18) Which i the exact ame reult we got earlier. The magnetic preure i independent of the actual ize of the urface. The only place the phyical dimention of the urface the preure i exerted on i preent in the B field. Getting back to our problem, the magnetic field from an ideal olenoid i given a B olenoid = µ 0 ni(t)ẑ (19) where n i the number of turn per unit length, n = N. Here, the field i only dependent l on the number of turn and the total length of the coil. Now we have the preure a a function of time: f magnetic = µ 0 n I (t max ) = µ 0n e Rt Vo L ωol (in ω ot) (0) And the peak preure i calculated uing the peak current a calculated above: f peak = 40 MPa = 5805 pi. Thi i plenty of preure. It i certainly enough to cruh Coke can, if not more. Quarter placed in thi olenoid would alo probably hrink. 4

In order to get a good etimate for the actual preure, we can make a few approximation. We know that the in term of the preure will reach maximum when the t = t max. So, we let that equal 1. Alo, we know that the exponential term will be very cloe to 1, ince the t max i o mall. We are left with: f magnetic µ 0n Vo (1) ωol Now, ince we are till making approximation we can ay that ω 0 ( LC) 1. Plugging thi in we have: f magnetic = µ 0n V o C = µ ( ) 0C N ( ) l V C V0 o = () πl R for a olenoid. Notice that the radiu of the coil R i in the denominator. Any error in that meaurment could yield very different preure. Perhap the better method of predicting the preure i with the inductance meaured directly, perhap the next-tolat relation. Alo, thi i independent of the number of turn in the coil. To maximize preure, one would have to minimize the radiu firt, then the length. 5