x
Our website uses cookies. By using the website you agree ot its use. More information can be found in our privacy policy.

Reduced Input

$-----------------------------------------------------------------------------
$
$ Example provided by Pierre (LSTC)
$
$ E-Mail: info@dynamore.de
$ Web: http://www.dynamore.de
$
$ Copyright, 2015 DYNAmore GmbH
$ Copying for non-commercial usage allowed if
$ copy bears this notice completely.
$
$X------------------------------------------------------------------------------
$X
$X 1. Run file as is.
$X    Requires LS-DYNA Beta version 119000 (or higher) with double precision 
$X
$X------------------------------------------------------------------------------
$# UNITS: (kg/m/s) 
$X------------------------------------------------------------------------------
$X
*KEYWORD
*TITLE
EM: Fall of a sphere on a 10 cells battery
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
$                                                                              $
$                            ELECTROMAGNETISM                                  $
$                                                                              $
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
*EM_CONTROL
$--------1---------2---------3---------4---------5---------6---------7---------8
$    emsol     numls  emdtinit   emdtmax   emtinit    emtend  ncyclFem  ncyclBem
         3              &em_dt                                    5000      5000
*EM_MAT
$--------1---------2---------3---------4---------5---------6---------7---------8
$   em_mid     mtype     sigma     eosId                              randletype
        11         2      6.e7                                                 1
        12         1                                                           2
        13         1                                                           3
        14         1                                                           4
        15         2      3.e7                                                 5
         1         2      1.e7
         2         2      1.e7
*SET_PART
1
1,2,3,4,5,6,7,8
9,10
*EM_RANDLES_LAYERED
$--------1---------2---------3---------4---------5---------6---------7---------8
$ randleId randlType partSetId   rdlArea
         1         1         1         2
$        Q        cQ   SOCinit    SOCtoU
      150.  2.777e-2      100.      -444
$    r0cha     r0dis    r10cha    r10dis    c10cha    c10dis
      0.02      0.02     0.008     0.008      110.      110.
$     temp fromTherm r0ToTherm      dUdT
       25.         0         1         0
$  useSocS   tauSocS  lcidSocS
         0        0.         0
*EM_ISOPOTENTIAL
$--------1---------2---------3---------4---------5---------6---------7---------8
$    isoId   setType     setId  randType
         1         2         1         1
*EM_ISOPOTENTIAL
$    isoId   setType     setId  randType
         2         2         2         1
*EM_ISOPOTENTIAL
$    isoId   setType     setId  randType
         3         2         3         1
*EM_ISOPOTENTIAL
$    isoId   setType     setId  randType
         4         2         4         1
*EM_ISOPOTENTIAL
$    isoId   setType     setId  randType
         5         2         5         1
*EM_ISOPOTENTIAL
$    isoId   setType     setId  randType
         6         2         6         1
*EM_ISOPOTENTIAL
$    isoId   setType     setId  randType
         7         2         7         1
*EM_ISOPOTENTIAL
$    isoId   setType     setId  randType
         8         2         8         1
*EM_ISOPOTENTIAL
$    isoId   setType     setId  randType
         9         2         9         1
*EM_ISOPOTENTIAL
$    isoId   setType     setId  randType
        10         2        10         1
*EM_ISOPOTENTIAL
$    isoId   setType     setId  randType
        11         2        11         5
*EM_ISOPOTENTIAL
$    isoId   setType     setId  randType
        12         2        12         5
*EM_ISOPOTENTIAL
$    isoId   setType     setId  randType
        13         2        13         5
*EM_ISOPOTENTIAL
$    isoId   setType     setId  randType
        14         2        14         5
*EM_ISOPOTENTIAL
$    isoId   setType     setId  randType
        15         2        15         5
*EM_ISOPOTENTIAL
$    isoId   setType     setId  randType
        16         2        16         5
*EM_ISOPOTENTIAL
$    isoId   setType     setId  randType
        17         2        17         5
*EM_ISOPOTENTIAL
$    isoId   setType     setId  randType
        18         2        18         5
*EM_ISOPOTENTIAL
$    isoId   setType     setId  randType
        19         2        19         5
*EM_ISOPOTENTIAL
$    isoId   setType     setId  randType
        20         2        20         5
*EM_ISOPOTENTIAL
$    isoId   setType     setId  randType
        21         2        21
*EM_ISOPOTENTIAL
$    isoId   setType     setId  randType
        22         2        22
*EM_ISOPOTENTIAL_CONNECT
$--------1---------2---------3---------4---------5---------6---------7---------8
$   connid  connType isoPotId1 isoPotId2     R,V,I      lcid
         1         3        22                  0.
*EM_ISOPOTENTIAL_CONNECT
$   connid  connType isoPotId1 isoPotId2     R,V,I      lcid
         2         4        21        22                 555
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
$                                                                              $
$                             SHORT CIRCUIT                                    $
$                                                                              $
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
*EM_RANDLE_SHORT
$--------1---------2---------3---------4---------5---------6---------7---------8
$ areaType   functId
         2       501
*DEFINE_FUNCTION_TABULATED
$#     fid     definition
       502     (thick,res) pair data
$# title
$--------1---------2---------3---------4---------5---------6---------7---------8
resistanceVsThickSep
                  0.               5.e-5
               1.e-3               5.e-5
               3.e-3               5.e-5
                  2.               5.e-5
                1.e2               5.e-5
*DEFINE_FUNCTION
501
float resistance_short_randle(float time,
                             float x_sep, float y_sep, float z_sep,
                             float x_sem, float y_sem, float z_sem,
                             float x_ccp, float y_ccp, float z_ccp,
                             float x_ccm, float y_ccm, float z_ccm)
{ 
  float distCC;
  float distSEP;
  distCC=sqrt(
     pow((x_ccp-x_ccm),2.)+pow((y_ccp-y_ccm),2.)+pow((z_ccp-z_ccm),2.));
  distSEP=sqrt(
     pow((x_sep-x_sem),2.)+pow((y_sep-y_sem),2.)+pow((z_sep-z_sem),2.));
  if (distCC < 0.000117) {
     return resistanceVsThickSep(distSEP) ;
  } else {
    return -1. ;
  }
}
*EM_OUTPUT
$--------1---------2---------3---------4---------5---------6---------7---------8
$     matS      matF      solS      solF      mesh    memory    timing    d3plot
         2         2         2         2                                       0
$      mf2       gmv                                            randle
                   1                                                 0
*END