Main >>
Applications >>
Sample problems >>
Relay dynamics simulation using GNU Octave and QuickField Object Interface

# Relay dynamics simulation using GNU Octave and ActiveField (QuickField Object Model)

The relay consists of the solenoid with moving core, which disconnects the electric switch contacts when the control current is running in the solenoid. The spring keeps the core in the pull out position with switch contacts connected. When the current is turned on, the magnetic field acts on the ferromagnetic core, overcomes the resistance of the spring, and pulls the core inside the solenoid to the pull in position where it is stopped by damper ring, which absorbs the shock. Operating time of this relay and the plunger motion function should be calculated.

To combine the electromagnetic field analysis with the moving core dynamics both GNU Octave and QuickField are used. Interaction between QuickField and GNU Octave is performed using ActiveField application programming interface.

**Problem Type:**

Axisymmetric problem of DC magnetics.

**Geometry:**

Number of turns *N* = 2000;

Current *I* = 0.2 A;

Plunger pull out position *x*_{out} = 10 mm.

Plunger pull in position *x*_{in} = 6 mm.

Plunger weight *m* = 4.5 g;

Spring constant *k* = 4 N/m

Spring free position *x*_{spring.free} = 15 mm.

**Problem:**

The electromagnetic and spring forces act on the plunger. Both forces depends on the plunger position. Calculate plunger motion function.

**Solution:**

The multi-turn winding is replaced with the equivalent total current.

The motion function can be found from second-order differential equation

*m* · *d*^{2}*x*/*d**t*^{2} = *f*(x),

where *m* - is a plunger weight (kg),

*x* - is a plunger position (m)

*f*(x) - is the force acting on the plunger (N).

The force acting on the plunger is a sum of spring force *f*_{spring}(*x*) = *k*·(*x*_{spring.free} - *x*) and electromagnetic force.

The equations are solved in GNU Octave. The dynamic link is used to invoke QuickField and calculate the electromagnetic force at each step.

The calculations are stopped when *x*=*x*_{in} (pull in position, plunger hits damper).

View Octave script in the separate window: relay_dynamics_octave.m, getQfForce_X.m.

**Results:**

The plunger hits initial position after 55 ms.

Video.

Download simulation files.