Linear Dipole Model

The linear dipole model is mechanically the same as the Linear Model, comprised of both linear and dashpot components. In addition to unbonded mechanical interactions, forces and moments based on discrete magnetic dipoles may be included to explicitly model permanent and induced magnetization.

Introduction

The linear dipole model provides mechanical linear and dashpot components that act in parallel with one another. The linear component provides linear elastic (no-tension) frictional behavior, while the dashpot component provides viscous behavior. Both components act over a vanishingly small area, and thus transmit only a force. The mechanical model is thoroughly elucidated in the Linear Model document, and the user is referred there for further details of the mechanical characteristics of the model.

In order to model magnetic interactions, each piece is abstracted as a magnetic dipole based on [Thomaszewski2008]. Forces and moments due to the dipole interactions are computed and also applied to the pieces. The moments applied to each piece due to the dipole interactions are, in general, not the same. One must also use the proximity concept (see the contact cmat proximity command) so that contacts may be deemed active despite the lack of physical overlap. Inheritance of the dipole_mo surface property is also required along with orientation tracking (see the model orientation-tracking command).

Behavior Summary

The linear dipole model provides the behavior of an infinitesimal interface that does not resist relative rotation. The mechanical contact force is resolved into linear and dashpot components. The linear component provides linear elastic (no-tension), frictional behavior. The dashpot component provides viscous behavior. Review the Linear Model document for further details of the mechanical response. Additional forces and moments due to interacting magnetic dipoles are also added. The total contact force is Fc=Fl+Fd+Fdi, where Fl is the linear force, Fd is the dashpot force, and Fdi is the magnetic dipole force. The moments applied to each piece from the magnetic dipoles are generally different from each other and different from zero (i.e., M1M20, where M1 is the moment applied to piece 1 and M2 is the moment applied to piece 2). These forces and moments are based on magnetostatic computations assuming dipoles as outlined below.

Activity-Deletion Criteria

A contact with the linear dipole model is active if the surface gap is less than or equal to the dipole distance dd (0 by default). The force-displacement law is skipped for inactive contacts. Mechanical forces may not arise even if dipole forces/moments exist.

Force-Displacement Law

In addition to the force displacements computed in the linear model, this model computes forces and moments due to magnetic dipoles interacting at a distance. These forces are added to the linear and dashpot forces.

Given the magnetic field B, the magnetizing field H, the current density J, and the magnetization M, Maxwell’s equations can be written as

(1)B=0
(2)×H=J
(3)H=1μ0BM

where μ0 is the permeability of free space and is the gradient.

Suppose there is a magnetic moment m positioned at the center of a rigid piece given a specified current density. Assuming this can be represented as a vanishingly small magnetic dipole, the force exerted on the rigid body in the presence of a magnetic field is given by

(4)Fdi=(mB)

while the moment is

(5)Mdi=m×B

The magnetic field as a function of distance from a vanishingly small magnetic dipole is

(6)B(r)=μ04π[3n(nm)m(|rr0|)3]

where n=(rr0)/|rr0| and r0 is the center of the rigid piece.

As noted in [Thomaszewski2008], for an arbitrarily shaped object with uniform magnetization, one can decompose it into N magnetic dipoles so that the total magnetic field is

(7)B(r)=μ04πNi=1[3ni(nimi)mi(|rri|)3]

For an infinite number of dipoles, the solution converges to the total magnetic field if their spatial distribution discretizes the piece shape. As the shapes of the pieces/bodies are not incorporated explicitly in the model, one can discretize a body into multiple clump pebbles to simulate non-spherical magnet geometries using this contact model.

The forces/moments represent the interactions of two magnetic dipoles with magnetic moments m1 and m2 at positions r1 and r2; these are the piece positions. Given the definition of the magnetic field due to a dipole, the resulting force is

(8)Fdi=3μ04π|r1r2|4[n12(m1m2)+m1(m2n12)+m2(m1n12)5n12(m1n12)(m2n12)]

where n12=(r1r2)/|r1r2|. This is the force applied to piece 2. The force applied to piece 1 is equal and opposite, replacing n12 with n21. As the distance between pieces shrinks to a vanishingly small value the force becomes infinite. The dcap is used to place a cap on the minimum distance used for force and moment computations.

The moments are

(9)M2=μ04π|r1r2|3[3(m2×n12)(m1n12)m2×m1]

and

(10)M1=μ04π|r1r2|3[3(m1×n21)(m2n21)m1×m2]

These moments are not equal and opposite. Also notice that the moments decay less rapidly than the forces, meaning that moments dominate at larger distances.

The above forces and moments presume the magnetic moments are in the local coordinate system. They should be specified as a piece surface property and inherited. Since these are vectors, though, they must be rotated to account for the rigid body rotation of the pieces. Thus, at the point in the simulation when they are specified as piece surface properties, it is presumed no rotations have occurred and orientation tracking is active. Subsequently they are rotated as the pieces spin.

Energy Partitions

The energies computed are the same as for the linear energies.

Properties

The properties defined by the linear dipole contact model are listed in the table below as a concise reference. See the Contact Properties section for a description of the information in the table columns. The mapping from the surface inheritable properties to the contact model properties is also discussed below.

Table 1: Linear Dipole Model Properties
Keyword Symbol Description Type Range Default Modifiable Inheritable
linear Model name
Linear Group:
kn kn Normal stiffness [force/length] FLT [0.0,+) 0.0 YES YES
ks ks Shear stiffness [force/length] FLT [0.0,+) 0.0 YES YES
fric μ Friction coefficient [-] FLT [0.0,+) 0.0 YES YES
rgap gr Reference gap [length] FLT R 0.0 YES NO
lin_mode Ml Normal-force update mode [-] INT {0,1} 0 YES NO
    {0: update is absolute1: update is incremental          
emod E Effective modulus [force/area] FLT [0.0,+) 0.0 NO N/A
kratio κ Normal-to-shear stiffness ratio [-] FLT [0.0,+) 0.0 NO N/A
    κknks          
lin_slip s Slip state [-] BOOL {false,true} false NO N/A
    {true: slippingfalse: not slipping          
lin_force Fl Linear force (contact plane coord. system) VEC R3 0 YES NO
    (Fln,Flss,Flst)(2D model: Flss0)          
user_area A Constant area [length*length] FLT (0.0,+) 0.0 YES NO
Dashpot Group:
dp_nratio βn Normal critical damping ratio [-] FLT [0.0,1.0] 0.0 YES NO
dp_sratio βs Shear critical damping ratio [-] FLT [0.0,1.0] 0.0 YES NO
dp_mode Md Dashpot mode [-] INT {0,1,2,3} 0 YES NO
    {0: full normal & full shear1: no-tension normal & full shear2: full normal & slip-cut shear3: no-tension normal & slip-cut shear          
dp_force Fd Dashpot force (contact plane coord. system) VEC R3 0 NO NO
    (Fdn,Fdss,Fdst)(2D model: Fdss0)          
Dipole Group:
dipole_cap dcap Capping distance for force/moment [length] FLT (0.0,+) 0.0 YES NO
dipole_d dd Activity distance [length] FLT (0.0,+) 0.0 YES NO
dipole_f Fdi Dipole force (contact plane coord. system) VEC R3 0 YES NO
dipole_m1 M1 Moment applied to piece 1 (contact plane coord. system) VEC R3 0 YES NO
dipole_m2 M2 Moment applied to piece 2 (contact plane coord. system) VEC R3 0 YES NO
dipole_mo1 m1 Dipole moment of piece 1 VEC R3 0 YES YES
dipole_mo2 m2 Dipole moment of piece 2 VEC R3 0 YES YES
By convention, κ=0.0 if either the normal or the shear stiffness is 0.

Note

Out-of-balance forces acting on bodies are accumulated during force-displacement calculations. Modifying the forces stored in the contact models will not alter them instantaneously. Therefore, any change to Fl may only be effective during the next force-displacement calculation.

Surface Property Inheritance

The linear stiffnesses, kn and ks, and the friction coefficient, μ, may be inherited from the contacting pieces, along with the original dipole moments m1 and m2. Remember that for a property to be inherited, the inheritance flag for this property must be set to true (default), and both contacting pieces must hold a property with the exact same name.

The linear stiffnesses are inherited, assuming that both pieces’ stiffnesses act in series (Figure 5):

(11)1kn=1k(1)n+1k(2)n1ks=1k(1)s+1k(2)s

where (1) and (2) denote the properties of piece 1 and 2, respectively.

The friction coefficient is inherited using the minimum of the values set for the pieces:

(12)μ=min(μ(1),μ(2)).

The dipole moments are taken from each piece and stored.

Methods

Table 2: Linear Dipole Model Methods
Method Arguments Symbol Type Range Default Description
area Set user_area to the area
deformability Set deformability
  emod E FLT [0.0,+) N/A Effective modulus
  kratio κ FLT [0.0,+) N/A Normal-to-shear stiffness ratio
By convention, setting κ=0.0 sets ks=0.0 but does not alter kn.

See here for a description of the methods.

Callback Events

Table 3: Linear Dipole Model Callback Events
Event Array Slot Value Type Range Description
contact_activated Contact has become active
  1 C_PNT N/A Contact pointer
slip_change Slip state has changed
  1 C_PNT N/A Contact pointer
  2 INT {0,1} Slip change mode
        {0: slip has initiated1: slip has ended

Usage and Verification Examples

The verification “Simulating Square Permanent Magnets” compares the analytic response of the attractive forces/moments of a square magnet with increasing levels of discretization of a cube with dipoles.

Model Summary

An alphabetical list of the linear dipole model methods is given here. An alphabetical list of the linear dipole model properties is given here.

References

[Thomaszewski2008](1, 2) Thomaszewski, B., A. Gunmann, S. Pabst and W. Straßer. “Magnets in Motion,” in ACM Transactions on Graphics, 27, pp. 1-9, 2008.