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