LOAD_FORCE_INTERACT

Boundary conditions
Attention: This command is in the beta stage and the format may change over time.
*LOAD_FORCE_INTERACT
"Optional title"
coid
entype${}_1$, enid${}_1$, entype${}_2$, enid${}_2$, fid, dtype, dynamic
Parameter definition
VariableDescription
coid Command ID (optional)
entype${}_1$ Entity type of first body
options: P, PS
enid${}_1$ Entity identification number of first body
entype${}_2$ Entity type of second body
options: P, PS
enid${}_2$ Entity identification number of second body
fid Force FUNCTION
dtype Definiton of distance between bodies
options:
0 $\rightarrow$ force based on distance between neareast elements
1 $\rightarrow$ force based on distance between centers of gravity
dynamic Flag to activate dynamic search for nearest elements
options:
0 $\rightarrow$ use nereast elements from time zero
1 $\rightarrow$ use current geometry
Description

This command is used to apply forces beteen two bodies. The force per unit body volume is defined by a FUNCTION with ID=fid. The total force is obtained by integrating the force intensity over the volume of the first body.

The relative location of the bodies defines the direction of the force. The magnitude of the force is allowed to depend on their relative distance (intrinsic variable dist). There are two different relative distance options, dtype=0 and dtype=1.

If dtype=0 then we always find the nearest point on the second body for every point (element center) in the first body. If dtype=1 we do no work with local distances, but with the relative location of the bodies center of gravity. Note that if dtype=1 and the interacting bodies are slender (e.g. cables) then the local centers of gravity of the cross sections are used.

Example
Two magnets

A model with two magnets separated by a board.


*PARAMETER
R = 0.01, "magnet radius"
h = 0.01, "magnet thickness"
t = 0.01, "table thickness"
C = 1000.0, "magnet constant"
tend = 0.11, "end time"
*TIME
[%tend]
#
# --- MESH ---
#
*COMPONENT_CYLINDER
"magnet 1"
1, 1, 1, 4
0, 0, [-%h-%t], 0, 0, [-%t], [%R]
*COMPONENT_CYLINDER
"magnet 2"
2, 2, 1, 4
0, 0, 0, 0, 0, [%h], [%R]
*COMPONENT_BOX
"table"
3, 3, 10, 10, 1
-0.1, -0.1, [-%t], 0.1, 0.1, 0.0
*CHANGE_P-ORDER
ALL, 0, 2
*SMOOTH_MESH
ALL, 0, 45.0 ,
#
# --- MATERIAL ---
#
*MAT_RIGID
"iron"
1, 7800.0
*MAT_RIGID
"wood"
2, 1000.0
*PART
"magnet 1"
1, 1
"magnet 2"
2, 1
"table"
3, 2
#
# --- BC ---
#
*BC_MOTION
"magnet 2"
2
P, 2
V, X, 1
*FUNCTION
1
1.0
*BC_MOTION
"table"
3
P, 3, XYZ, XYZ
#
# --- MAGNETIC FORCES ---
#
*LOAD_FORCE_INTERACT
"magnets"
1
P, 1, P, 2, 1000, 0, 0
*FUNCTION
"magnetic force"
1000
-%C/dist^2
#
# --- GRAVITY ---
#
*LOAD_GRAVITY
Z, 981
*FUNCTION
981
9.81
#
# --- CONTACT ---
#
*CONTACT
"general"
1
ALL, 0, ALL, 0, 0.1, 1.0e11
*END