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.
*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
*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
