Laminate sheet closure using magnets — Welcome to LS-DYNA Examples

Home » Knowledge Base » EM / Electromagnetics » Tutorial Eddy Currents » Laminate sheet closure using magnets

In this simulation, several magnets are embedded in a soft laminate structure. The objective is to study the latching of the two structure together, based on the positioning of the different magnets. It uses the FEM-BEM monolithic magnetostatics solver (used for zero-conductivity ferromagnetic materials) and a AMS preconditionner.

Description

In this simulation, several magnets are embedded in a soft laminate structure. The objective is to study the latching of the two structure together, based on the positioning of the different magnets. It uses the FEM-BEM monolithic magnetostatics solver (used for zero-conductivity ferromagnetic materials) and a AMS preconditionner.

Animated Result

EM surface magnetic forces :

https://youtu.be/U7pyZlqxqy8

Keywords

*CONTROL_TERMINATION
*CONTROL_TIMESTEP
*DATABASE_BINARY_D3PLOT
*DEFINE_CURVE_TITLE
*EM_CONTROL
*EM_CONTROL_COUPLING
*EM_CONTROL_TIMESTEP
*EM_MAT_001
*EM_MAT_002
*EM_OUTPUT
*EM_PERMANENT_MAGNET
*EM_SOLVER_BEM
*EM_SOLVER_BEMMAT
*EM_SOLVER_FEM
*EM_SOLVER_FEMBEM_MONOLITHIC
*END
*KEYWORD
*MAT_RIGID
*PARAMETER
*PART
*SECTION_SOLID
*TITLE

Reduced Input

*KEYWORD
$ ****************************************************************************
$ electromagnetism
$ ****************************************************************************
*EM_CONTROL
$--------1---------2---------3---------4---------5---------6---------7---------8
$    emsol     numls        dt
         1                                                  &em_femmtx&em_bemmtx
*EM_CONTROL_TIMESTEP
1,&emdt
*EM_SOLVER_FEM
$--------1---------2---------3---------4---------5---------6---------7---------8
$   reltol   maxiter     stype    precon   uselast  ncyclfem
     1.e-6     10000         1         1         1
*EM_SOLVER_BEM
$--------1---------2---------3---------4---------5---------6---------7---------8
$   reltol   maxiter solvetype    precon   uselast  ncyclbem
     1.e-6     10000         2         2         1
$*EM_SOLVER_FEMBEM
$$--------1---------2---------3---------4---------5---------6---------7---------8
$$   reltol   maxiter    forcon
$     1.e-6        50                             3
*EM_SOLVER_BEMMAT
$--------1---------2---------3---------4---------5---------6---------7---------8
$    matId                                                                relTol
         1                                                                1.e-12
*EM_SOLVER_BEMMAT
$--------1---------2---------3---------4---------5---------6---------7---------8
$    matId                                                                relTol
         2                                                                1.e-12
*EM_SOLVER_BEMMAT
$--------1---------2---------3---------4---------5---------6---------7---------8
$    matId                                                                relTol
         3                                                                1.e-12
*EM_MAT_001
$--------1---------2---------3---------4---------5---------6---------7---------8
$      mid     mtype     sigma     eosid
         1         1
*EM_MAT_001
$--------1---------2---------3---------4---------5---------6---------7---------8
$      mid     mtype     sigma     eosid
         2         1
*EM_MAT_002
$ Magnetic material 1
$--------1---------2---------3---------4---------5---------6---------7---------8
$      mid     mtype     sigma     eosid     murel     eosmu
         3         4    &sigma         0    &murel
*EM_MAT_002
$ Magnetic material 1
$--------1---------2---------3---------4---------5---------6---------7---------8
$      mid     mtype     sigma     eosid     murel     eosmu
         4         4    &sigma         0    &murel
*EM_PERMANENT_MAGNET
$--------1---------2---------3---------4---------5---------6---------7---------8
$       ID   Part ID         -     North     South         M
         1         1                   2         3      9.e5
*EM_PERMANENT_MAGNET
$--------1---------2---------3---------4---------5---------6---------7---------8
$       ID   Part ID         -     North     South         M
         2         4                   5         6      9.e5
*EM_PERMANENT_MAGNET
$--------1---------2---------3---------4---------5---------6---------7---------8
$       ID   Part ID         -     North     South         M
         3         7                   8         9      9.e5
*EM_PERMANENT_MAGNET
$--------1---------2---------3---------4---------5---------6---------7---------8
$       ID   Part ID         -     North     South         M
         4        10                  11        12      9.e5
*EM_PERMANENT_MAGNET
$--------1---------2---------3---------4---------5---------6---------7---------8
$       ID   Part ID         -     North     South         M
         5        13                  14        15      9.e5
*EM_PERMANENT_MAGNET
$--------1---------2---------3---------4---------5---------6---------7---------8
$       ID   Part ID         -     North     South         M
         6        16                  17        18      9.e5
*EM_PERMANENT_MAGNET
$--------1---------2---------3---------4---------5---------6---------7---------8
$       ID   Part ID         -     North     South         M
         7        19                  20        21      9.e5
*EM_PERMANENT_MAGNET
$--------1---------2---------3---------4---------5---------6---------7---------8
$       ID   Part ID         -     North     South         M
         8        22                  23        24      9.e5
*EM_PERMANENT_MAGNET
$--------1---------2---------3---------4---------5---------6---------7---------8
$       ID   Part ID         -     North     South         M
         9        25                  26        27      9.e5
*EM_PERMANENT_MAGNET
$--------1---------2---------3---------4---------5---------6---------7---------8
$       ID   Part ID         -     North     South         M
        10        28                  29        30      9.e5
*EM_OUTPUT
$--------1---------2---------3---------4---------5---------6---------7---------8
$     matS      matF      solS      solF      mesh    memory    timing d3plotAsc
         3         2         3         2         0         0         0         0
$      mf2       gmv d3plotFor  timeHist
         0         0         0
*EM_SOLVER_FEMBEM_MONOLITHIC
$--------1---------2---------3---------4---------5---------6---------7---------8
$ COUPLING KSP_TYPE    KSP_ATOL KSP_RTOL KSP_MAXIT GMRES_RST               DEBUG
         0              1.e-10     1.e-5     10000                   1
$--------1---------2---------3---------4---------5---------6---------7---------8
$    FORCE NEWT_STOL NEWT_ATOL NEWT_RTOL NEWT_MAXI
               1.e-4     1.e-4     1.e-4
$--------1---------2---------3---------4---------5---------6---------7---------8
$    LS_ON   LS_FTOL   LS_GTOL   LS_RTOL  LS_SAMPL LS_NUM_SP LS_MAXFUN
         0     1.e-3       0.9     1.e-4                            30
$--------1---------2---------3---------4---------5---------6---------7---------8
$USE_ENERG  STOP_ERR WIDTH_ERR    XTRAPF
        
*EM_CONTROL_COUPLING
,2
*END

Figures

 

bookclose.png

Download