Force on a Taper Plunger Magnet

David Meeker

Introduction

The objective of this note is to compare calculations of Force versus Position for a solenoid plunger with experimental results from Roter's "Electromagnetic Devices." In particular, the actuator pictured in Figure 1 below will be considered.

Figure 1: Roters’ Ch. IX, Fig. 7, "Taper Plunger Magnet and its force-stroke curves."

Model Construction

A model of this geometry is available as rotor1b.fem. The model is pictured below in Figure 2.

The actuator geometry was previously drawn in FEMM in the usual interactive mode. All components of the moving plunger were defined to be members of Group Number 1 so that the plunger can be easily selected and moved. The BH curve for the ½ hard 1020 steel was also obtained from the Roters text.

Comparison to Published Results

The data was stripped from the figure in Roters using WinDig, a free tool for recovering numerical data from graphs.

A combined plot, shown in Figure 2, compares the force for various plunger positions with at 1.5 A coil current as published in Roters to that predicted by the “weighted stress tensor” volume integral in FEMM. There is a reasonable agreement between the experimental and FEA results.

Figure 1: FEMM model of Roters’ “Taper Plunger Magnet”

Figure 1: Comparison of FEMM results to Roters example for 1.5 A current.

Matlab Script Implementation

A Matlab script, roters1b.m, is included in the regular FEMM distribution, for completeness, it is listed below. The script moves the plunger through its stroke, computing force at a number of intermediate locations.

disp('Roters1b: Simulation of a Tapered Plunger Magnet');
disp('David Meeker');
disp('dmeeker@ieee.org');
disp('This geometry comes from Chap. IX, Figure 7 of Herbert Roters');
disp('classic book, Electromagnetic Devices. The program moves');
disp('the plunger of the magnet over a stroke of 1.5in at 1/10in increments');
disp('solving for the field and evaluating the force on the plunger at');
disp('each position. When all positions have been evaluated, the program');
disp('plots a curve of the finite element force predictions.');

openfemm

opendocument('roters1b.fem');
mi_saveas('temp.fem');
n=16;
stroke=1.5;
x=zeros(n,1);
f=zeros(n,1);

for k=1:n
disp(sprintf('iteration %i of %i',k,n));
mi_analyze;
mi_loadsolution;
mo_groupselectblock(1);
x(k)=stroke*(k-1)/(n-1);
f(k)=mo_blockintegral(19)/4.4481;
mi_selectgroup(1);
mi_movetranslate(0,-stroke/(n-1));
mi_clearselected
end

plot(x,f)
xlabel('Displacement, Inches');
ylabel('Force, Lbf');
title('Plunger Force vs. Displacement');

closefemm

References