Study on the Calculation of Magnetic Force Based on the Equivalent Magnetic Charge Method

Similar documents
Studies on Nuclear Fuel Rod Thermal Performance

13.4 Work done by Constant Forces

KINEMATICS OF RIGID BODIES

Contact Analysis on Large Negative Clearance Four-point Contact Ball Bearing

ADVANCEMENT OF THE CLOSELY COUPLED PROBES POTENTIAL DROP TECHNIQUE FOR NDE OF SURFACE CRACKS

CBE 291b - Computation And Optimization For Engineers

State space systems analysis (continued) Stability. A. Definitions A system is said to be Asymptotically Stable (AS) when it satisfies

DESIGN OF FERRITE MAGNET CFSM TO IMPROVE SPEED RESPONSE OF EPS SYSTEM

ISSN. on Ansys. School of. of Shanghai. and. for Science. Abstract This paper. shield bracket. excavation begins, the. the tunneling.

Measuring Electron Work Function in Metal

NUMERICAL INTEGRATION. The inverse process to differentiation in calculus is integration. Mathematically, integration is represented by.

ELE B7 Power Systems Engineering. Power System Components Modeling

LABYRINTH SEALS DYNAMIC COEFFICIENTS AND CRITICAL SPEED RESEARCH BASED ON CHILDS MODEL

#6A&B Magnetic Field Mapping

Higher Checklist (Unit 3) Higher Checklist (Unit 3) Vectors

Reference. Vector Analysis Chapter 2

A REVIEW OF CALCULUS CONCEPTS FOR JDEP 384H. Thomas Shores Department of Mathematics University of Nebraska Spring 2007

Fig. 1. Open-Loop and Closed-Loop Systems with Plant Variations

Thermal Performance of Electrocaloric Refrigeration using Thermal Switches of Fluid Motion and Changing Contact Conductance

Entropy ISSN

ENERGY-BASED METHOD FOR GAS TURBINE ENGINE DISK BURST SPEED CALCULATION

Lesson 8. Thermomechanical Measurements for Energy Systems (MENR) Measurements for Mechanical Systems and Production (MMER)

Freely propagating jet

Physics 201 Lab 3: Measurement of Earth s local gravitational field I Data Acquisition and Preliminary Analysis Dr. Timothy C. Black Summer I, 2018

Estimation of the particle concentration in hydraulic liquid by the in-line automatic particle counter based on the CMOS image sensor

Multi-objective optimization of dielectric layer photonic crystal filter

Goals: Determine how to calculate the area described by a function. Define the definite integral. Explore the relationship between the definite

(See Notes on Spontaneous Emission)

63. Representation of functions as power series Consider a power series. ( 1) n x 2n for all 1 < x < 1

Math 8 Winter 2015 Applications of Integration

The momentum of a body of constant mass m moving with velocity u is, by definition, equal to the product of mass and velocity, that is

INTRODUCTION. The three general approaches to the solution of kinetics problems are:

Driving Cycle Construction of City Road for Hybrid Bus Based on Markov Process Deng Pan1, a, Fengchun Sun1,b*, Hongwen He1, c, Jiankun Peng1, d

Predict Global Earth Temperature using Linier Regression

WELCOME TO THE LECTURE

Numerical Analysis: Trapezoidal and Simpson s Rule

On the Uncertainty of Sensors Based on Magnetic Effects. E. Hristoforou, E. Kayafas, A. Ktena, DM Kepaptsoglou

Pressure Wave Analysis of a Cylindrical Drum

New Expansion and Infinite Series

Chapter 6 Electrostatic Boundary Value Problems. Dr. Talal Skaik

Flutter frequency based on bending - torsion coupling theory

Summary of equations chapters 7. To make current flow you have to push on the charges. For most materials:

First Law of Thermodynamics. Control Mass (Closed System) Conservation of Mass. Conservation of Energy

Simulation of Eclipsing Binary Star Systems. Abstract

Math 1B, lecture 4: Error bounds for numerical methods

Physics 3323, Fall 2016 Problem Set 7 due Oct 14, 2016

Flow in porous media

THE INTERVAL LATTICE BOLTZMANN METHOD FOR TRANSIENT HEAT TRANSFER IN A SILICON THIN FILM

Using air lines as references for VNA phase measurements

Forces from Strings Under Tension A string under tension medites force: the mgnitude of the force from section of string is the tension T nd the direc

Scientific notation is a way of expressing really big numbers or really small numbers.

QUADRATURE is an old-fashioned word that refers to

Calculus of Variations

This final is a three hour open book, open notes exam. Do all four problems.

Multiple Integrals. Review of Single Integrals. Planar Area. Volume of Solid of Revolution

CAPACITORS AND DIELECTRICS

Magnetic forces on a moving charge. EE Lecture 26. Lorentz Force Law and forces on currents. Laws of magnetostatics

Terminal Velocity and Raindrop Growth

CONTRIBUTION TO THE EXTENDED DYNAMIC PLANE SOURCE METHOD

Conservation Law. Chapter Goal. 5.2 Theory

Purpose of the experiment

Chapter 4 Contravariance, Covariance, and Spacetime Diagrams

A formula sheet and table of physical constants is attached to this paper. Linear graph paper is available.

APPROXIMATE INTEGRATION

Simulated Performance of Packed Bed Solar Energy Storage System having Storage Material Elements of Large Size - Part I

Vibrational Relaxation of HF (v=3) + CO

Energy creation in a moving solenoid? Abstract

FEM ANALYSIS OF ROGOWSKI COILS COUPLED WITH BAR CONDUCTORS

P 3 (x) = f(0) + f (0)x + f (0) 2. x 2 + f (0) . In the problem set, you are asked to show, in general, the n th order term is a n = f (n) (0)

Section 4.8. D v(t j 1 ) t. (4.8.1) j=1

Section 4: Integration ECO4112F 2011

Set up Invariable Axiom of Force Equilibrium and Solve Problems about Transformation of Force and Gravitational Mass

Chapters 4 & 5 Integrals & Applications

1 Bending of a beam with a rectangular section

KINEMATICS OF RIGID BODIES

7.1 Integral as Net Change and 7.2 Areas in the Plane Calculus

in a uniform magnetic flux density B = Boa z. (a) Show that the electron moves in a circular path. (b) Find the radius r o

Rolling Contact Bearings (pg 599)

In the diagram below, the rotation continues until N-S alignment, resulting in lock-up that is, if nothing is done to prevent it.

Abstract. Keywords. 1. Introduction. Eric S. Sanches, José A. Santisteban

5.7 Improper Integrals

Review of Calculus, cont d

Calculus - Activity 1 Rate of change of a function at a point.

Unit #9 : Definite Integral Properties; Fundamental Theorem of Calculus


MAT187H1F Lec0101 Burbulla

CHM Physical Chemistry I Chapter 1 - Supplementary Material

x = b a N. (13-1) The set of points used to subdivide the range [a, b] (see Fig. 13.1) is

The heat budget of the atmosphere and the greenhouse effect

Lecture 13 - Linking E, ϕ, and ρ

Chapter 0. What is the Lebesgue integral about?

Physics 202, Lecture 14

Design Data 1M. Highway Live Loads on Concrete Pipe

INVESTIGATION ON THE MODEL OF VORTEX-INDUCED

Properties of Integrals, Indefinite Integrals. Goals: Definition of the Definite Integral Integral Calculations using Antiderivatives

INTRODUCTION TO INTEGRATION

Physics 24 Exam 1 February 18, 2014

Vorticity. curvature: shear: fluid elements moving in a straight line but at different speeds. t 1 t 2. ATM60, Shu-Hua Chen

Applications of Bernoulli s theorem. Lecture - 7

Math& 152 Section Integration by Parts

Transcription:

Avilble online t www.sciencedirect.com Physics Procedi 4 () 9 97 Interntionl Conference on Applied Physics nd Industril Engineering Study on the Clcultion of Mgnetic Force Bsed on the Equivlent Mgnetic Chrge Method Jingng Li,,Qingchng Tn,Yongqi Zhng,Kuo Zhng College of Technology, Ynbin University, Ynji Chin College of Mechnicl Science nd Engineering,, Jilin University, Chngchun Chin Abstrct Mgnetic drivers hve been used widely in the phrmceuticl, chemicl, petroleum, food nd other industries with its perfect seling without contct. Common method of clculting of the mgnetic force re the Mxwell equtions, empiricl formuls, nd he equivlent mgnetic chrge method s well. The Mxwell equtions method is the most complicted nd the empiricl formuls method is the simplest with low ccurcy. The equivlent mgnetic chrge method is simpler thn the Mxwell equtions method nd more ccurte thn the empiricl formuls method. In this pper, the mgnetic force of the mgnetic driver of reciprocte in line is clculted with the equivlent mgnetic chrge method nd ws compred with the experiment. Published by Elsevier by Elsevier B.V. Ltd. Selection Selection nd/or nd/or peer-review peer-review under responsibility under responsibility of ICAPIE of Orgniztion [nme orgnizer] Committee. Open ccess under CC BY-NC-ND license. Keywords: Mgnetic driver; equivlent mgnetic chrge; mgnetic sleeve; fitting curve. Introduction The mgnetic drive is new trnsfer technology of non-contct force or torque using mgnetic force of permnent mgnet mterils or electromgnetic bsed on the theory of modern mgnetic. The device to chieve this technology clled mgnetic driver or mgnetic coupling. The mgnetic drive chieves the trnsmission of force nd torque using the feture of mgnetic field goes through the work spce or the thin distnce sleeve which is different from the mechnicl drive. Such s mgnetic drive pump or gittor, chnges dynmic sel to sttic sel through bolishing trditionl contct sels fter dopts the sel device of mgnetic drive. It elimintes the lekge of shft sel fundmentlly, so the run, drip, lek nd mechnicl sel het, wer of industril pumps, cldron nd other equipment is solved. It cretes fvorble conditions to trnsport of mixing of toxic, flmmble, explosive nd corrosive mterils sfely nd relibly []. The mgnetic driver of reciprocte in line is kind of mchine to trnsfer the motor s rottion to reciprocting motion in line []. It is very complicted to clculte the mgnetic force. Equivlent mgnetic 875-389 Published by Elsevier B.V. Selection nd/or peer-review under responsibility of ICAPIE Orgniztion Committee. Open ccess under CC BY-NC-ND license. doi:.6/j.phpro...9

Jingng Li et l. / Physics Procedi 4 () 9 97 9 chrge method is n effective wy for clcultion of mgnetic force. Tn Qingchng etc. hve done some effective work in this re [3][4]. Using Mtlb softwre, this pper clcultes mgnetic force of the mgnetic driver of reciprocte in line in detil nd derives chrcteristic curve between the xis mgnetic force nd the reltive displcement of inner mgnetic sleeve nd outer mgnetic sleeve. At lst, the result of clcultion is compred with the experiment.. Model of mgnetic force The mgnetic force s clcultion model of mgnetic driver of reciprocte in line is shown in Fig.. Mgnetiztion direction of outer mgnetic sleeve 3 4 R 4 3 y R 3 R y z Mgnetiztion direction of inner mgnetic sleeve R o o x + shown s positive mgnetic chrge - shown s negtive mgnetic chrge Figure. The clcultion model of mgnetic force Becuse of the xil symmetry of inner nd outer mgnetic sleeves, the rdil mgnetic force is zero. Therefore, this pper only studies the xis mgnetic force. In ccordnce with the theory of equivlent mgnetic chrge, the xis mgnetic force between inner nd outer mgnetic sleeves cn be expressed s [] Where A is given by: B B r r A R R4 F AArrdrdrdd R R3 Br is the residul mgnetic flux density of inner mgnetic sleeve, r density of outer mgnetic sleeve, is spce permebility which is defined s A is expressed s A, 3 3 3 A A ( ) A ( ) A is expressed s A r sin r sin r cos r cos. B is the residul mgnetic flux 7 4 N/A, 3. Clcultion of mgnetic force It needs to solve fourfold integrl for clcultion of mgnetic force. With the id of Mtlb softwre, this pper clcultes the mgnetic force. Becuse of there is only up to triple integrl for solving equtions

9 Jingng Li et l. / Physics Procedi 4 () 9 97 in Mtlb, this pper nlyzes eqution () bsed on the principle of differentil clculus, reduces eqution () to triple integrl then clcultes the mgnetic force with the use of Mtlb. This pper prtitions surfces nd (shown s fig.) into n equl prts long the dr direction, so r R ( i )dr, i n. i According to eqution (), when the reltive displcement between inner nd outer mgnetic sleeves (shown s ) is constnt numericl vlue, there is dr dr, so the mgnetic force (shown s F i ) between the i th differentil surfce nd the outer mgnetic sleeve cn be expressed s triple integrl: R4 F Ardr Ardrdd () i i R3 Where r r R ( i )dr in A i. y R d r Figure. Digrm of the differentition R o r i x According to reference [], B B.98 KGs 34 mm, R 5 mm, R 39.4 mm, r r R 45.4 mm, R 8 mm. Mgnetic force F 3 4 i cn be clculted with the commnd of triplequd, Mtlb procedure is s follows: zzi inline(' A ri * d r * r* A',' ',' ',' r'); Fi triplequd( zzi,,*pi,,*pi, R3, R4); Therefore, when the reltive displcement between inner nd outer mgnetic sleeves is constnt numericl vlue, the mgnetic force between inner nd outer mgnetic sleeves cn be expressed s follow: F F F... F F (3) n i i n According to eqution (3), this pper hs clculted mgnetic force of the mgnetic driver for different, shown s Fig.3. It cn be seen in Fig.3 tht the reltionship between xis force F nd, shown s follow: F when. Forwrd long the positive x direction, F is incresing rpidly with the increse of. Axis force F up to the mximum 97 N when. 7 mm. Then F decresed rpidly with the increse of, F when 64. mm. Forwrd long the negtive x direction, the lw of F nd is quite the opposite to the positive x direction.

Jingng Li et l. / Physics Procedi 4 () 9 97 93 Axil mgnetic force/n Reltive displcement between inner nd outer mgnetic sleeves /mm Figure 3. Performnce curve of reltive displcement between inner nd outer mgnetic sleeves nd xil force 4. Mgnetic force s fitting of the mgnetic driver Axis force is more big in the intervl [-5, +5] thn others, shown s the performnce curve in Fig.3. Otherwise, xis force is very smll in the intervl (-, -5] nd [+5, +). Axis force is symmetricl bout the origin of coordintes, so the fitting intervl of xis force cn be selected in (+5]. This pper obtins fitting curve through fourth-order fitting on xis force using Mtlb commnd of polyfit, the fitting eqution is: 4 3 Ff =.8 -.3-3. + 4 + 3.6 (4) The fitting curve of xis force is shown s Fig.4. Axis force/n Theoreticl vlue Clculted vlue Figure 4. Performnce curve of reltive displcement between inner nd outer mgnetic sleeves nd xil force Fitting error nlysis is shown in Tble. It cn be seen tht the xil force fitting errors 46. % in the intervl (, +4]. The reltive displcement between inner nd outer mgnetic sleeves should be controlled in the intervl [-, +] in order to get stble working stte for the mgnetic driver of reciprocte in line. The xil force showing good liner chrcter nd the fitting errors re less thn %, so the fitting result is perfect. TABLE THE FITTING ERRORS OF AXIS FORCE Items The reltive displcement between inner nd outer mgnetic sleeves /mm Theoreticl xis force F /N 66. 6.8 35. 394.6 345.9 58.8 857.4 56.4 The reltive displcement between inner nd outer mgnetic sleeves/mm

94 Jingng Li et l. / Physics Procedi 4 () 9 97 Fitting xis force F /N 3.6 6. 59.8 37. 398.8 33.9 3.8 84. 534.8 f Fitting errors.9%.%.9%.3%.6% 3.%.8% 4.6% 5. Experiments A double-cting hydrulic cylinder is used to provide the lod. The test tble is shown in Fig.5 (), the structure digrm is show in Fig.5 (b). 3 9 4 (b) Structure digrm of the test tble. Double-cting hydrulic cylinder. Restrictive vlve 3. Pull pressure sensor 4. distnce sleeve 5. inner mgnetic sleeve 6. outer mgnetic sleeve 7. Slider-crnk mechnism 8. Pin wheel Cycloid motor 9. Displcement sensor. Bridge box Figure 5. Photo nd structure digrm of test tble 5 6 S N N S () Photo of the test tble N S S N 7 DH5937/8 Dynmic Acquisition System 8 Mitsubishi Frequency Converter Computer 5. Sttic experiment Sttic experiment is pplying lod while the mgnetic driver is rest, then strt mgnetic driver t certin frequency nd mesure the mximum sttic lod when outer mgnetic sleeve prts from the inner. The results re shown in Fig.6, The sttic mgnetic crrying cpcity of the driver (xil force) into liner decrese with the strt frequency increses. 5. Dynmic experiment The crnk rdius is 6 mm, the frequency rnge of Mitsubishi Frequency Converter is limited from to 5 Hz (Becuse of serious crwling t low speed nd lrge vibrtion t high speed, ccording to the working condition of the test tble, so the frequency is limited), the opening of the restrictive vlve is controlled to keep full open or hlf open. Chrcteristic of displcement nd xis force between inner nd outer mgnetic sleeves t different frequency is shown in Fig.7. It cn be seen tht the reltionship between

Jingng Li et l. / Physics Procedi 4 () 9 97 95 Experimentl vlues Fitting line Axis Force/N Strt frequency/hz Figure 6. Sttic lod chrcteristic of the mgnetic driver Axis force/n Reltive displcement between inner nd outer mgnetic sleeves/mm Figure 7. Chrcteristic of displcement nd xis force between inner nd outer mgnetic sleeves t different frequency Axis force/n Restrictive vlve is hlf open Restrictive vlve is full open Reltive displcement between inner nd outer mgnetic sleeves/mm Figure 8. Chrcteristic of displcement nd xis force between inner nd outer mgnetic sleeves t different openning of restrictive vlve reltive displcement nd xis force is independent of the frequency.

96 Jingng Li et l. / Physics Procedi 4 () 9 97 Chrcteristic of displcement nd xis force between inner nd outer mgnetic sleeves t different open of restrictive vlve is shown in Fig.8. It cn be seen tht the reltionship between xis force nd reltive displcement hs nothing to do with the opening of restrictive vlve Axis force/n Experimentl vlues Theoreticl curve Reltive displcement between inner nd outer mgnetic sleeves/mm Figure 9. Comprison of theoreticl nd experimentl Fig.9 shows the comprison of theoreticl nd experimentl. It cn be seen tht the experimentl vlues nd theoreticl curve re consistent with the trend, but there is some discrepncy, minly due to the following points: Friction between the inner mgnetic sleeve nd the distnce sleeve. Although the theoreticl clcultion considers no rdil mgnetic force on the ssumption tht the xes of inner nd outer sleeves hve common xis, but in the mnufcturing, instlltion process hve inevitbly resulted in different xes, or xis devition hs occurred, resulting in some of the rdil mgnetic force which leding to incresed friction. Permnent mgnet demgnetiztion. The residul mgnetic flux density B r used in theoreticl clcultion is the fctory clibrtion vlue. Due to demgnetiztion chrcteristic of permnent mgnet, so the residul mgnetic flux density should be less thn the fctory clibrtion vlues. Eddy current losses. It will produce induced current if the distnce sleeve is cut by the lternting mgnetic field, nd the energy will turn into het. This will cuse eddy current heting to the mgnetic driver nd decrese its efficient, the energy trnsforms into lot of het which cn probbly cuse the distnce sleeve overheting, mgnetic sleeves demgnetiztion or other filure [5][6]. Other resons which cn decrese the efficiency of the mgnetic driver, including temperture, electromgnetic fields in the surrounding, etc. Conclusions In this pper, the method to clculte the xis force of mgnetic driver of reciprocte in line is studied, nd some experiments is done. The conclusions re s follows: Using equivlent mgnetic chrge theory to clculte the mgnetic force of mgnetic driver of reciprocte in line is very effective. In certin rnge of the reltive displcement between inner nd outer sleeves, mgnetic force cn be reduced to simpler clcultion model, shown s eqution (4). The reltionship between the reltive displcement of inner nd outer mgnetic sleeves nd xis force re independent of the frequency nd the opening of restrictive vlve.

Jingng Li et l. / Physics Procedi 4 () 9 97 97 References [] [] [3] [4] [5] [6] Zhng Suolong. STRUCTURE IMPROVEMENT OF MAGNETIC ACTUATION DEVICE FOR KETTLEPETRO- CHEMICAL EQUIPMENT TECHNOLOGY. 999, (5), pp.7-8. ZHAO Guoto, TAN Qingchng, LI Wei. MAGNETIC FORCE TRANSMISSION TECHNOLOGY OF RECIPROCATE IN LINE. CHINESE JOURNAL OF MECHANICAL ENGINEERING. 7, 43(), pp.-4. Tn Qingchng, Meng Huiqin, Xiu shicho. AN INVESTIGATION ON THE DIRECTIONS OF MAGNETIZATION OF PERMANENT-MAGNET BEARINGS. ACTA ARMAMENTARII. 993, 5(), PP.59-63. Tn Qingchng, Xin Di, Li Wei. Study on trnsmitting torque nd synchronism of mgnet couplings. Journl of Mechnicl Engineering Science. 99, 6(), PP.38-384 LI Zhipeng, LIU Xu. Energy Loss Anlysis nd Design Study of Mgnetic Pump. MECHANICAL AND ELECTRICAL EQUIPMENT. (), pp.8-3. FENG Zhongming, CHEN Cundong. Clcultion nlysis of turbulence loss in mgnetic pump & its ppliction. Generl Mchinery. 5(4), pp.8-8387.