Learning Forum

[Changping Yi]

Hi,

I downloaded input files for battery modelling from https://lsdyna.ansys.com/battery/#download. Some of them cannot be run due to different problems. 

For the input file of Basic_intshort, it shows *** Error 10246 (KEY+246)
     line contains improperly formatted data
      *** Error reading *EM_MAT card.
          At line# 89 of file
 --------------------------------------------------------------------------
          4         1
 --------------------------------------------------------------------------
 
 *** Error 10450 (KEY+450)
     in keyword command:
         5         2    3.70e7
          At line# 90 of file
For the case of  Basic_socshift,   it shows

 *** Error 10246 (KEY+246)
     line contains improperly formatted data
      *** Error reading *EM_MAT card.
          At line# 67 of file

 --------------------------------------------------------------------------
          4         1
 --------------------------------------------------------------------------
 
 *** Error 10450 (KEY+450)
     in keyword command:
         5         2    3.70e7
          At line# 68 of file
 For the case of Meshless, it shows

 *** Error 10307 (KEY+307)
     CHECKING MATERIAL INPUT Part ID= 1
     PART ID 1 has TMID= 0 which does not exist.
     This is PART 22 in the order of input.
     Coupled structural thermal analysis needs to have a thermal material ID defined.

There are different errors with other cases. I am using R16 solver, I wonder if these input files are old or other reasons. Thank you!

[ErKo]

 

 

 

 

 

 

 

Hello 

Thanks for letting us know (indeed these examples are old)– much appreciated.

 

For the first (basic_intshort) and second (basic_socshift) example you mention (https://lsdyna.ansys.com/wp-content/uploads/2022/12/basic_intshort.zip — https://lsdyna.ansys.com/wp-content/uploads/2022/12/basic_socshift.zip) – same issue is for the meshless example.

the em.k file should be like below (change that in that em file – format it as needed):
For ref., the main issue is the *EM_MAT_001 keyword. We need to have 5 individual *EM_MAT_001 keywords like shown below (something has changed in recent releases so we can not have all 5 properties under one *EM_MAT_001 card) – I have made a note for us internally about this issue – many thanks for bringing it up.

So the em.k should be like the one here: https://lsdyna.ansys.com/wp-content/uploads/2022/12/basic.zip (at least with respect to using multiple em_mat cards as below)

$–EM.k–

*KEYWORD
*INCLUDE
soc.k
$—+—-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                                                        5000      5000
*EM_CONTROL_TIMESTEP
1,&em_dt
*EM_RANDLES_SOLID
$——–1———2———3———4———5———6———7———8
$#   rdlid   rdltype   rdlarea   ccppart   ccnpart   seppart poselpart negelpart
         1         0         2         1         5         3         2         4
$#       q        cq   socinit    soctou      
       20.   0.02777     100.0      -444
$#   r0cha     r0dis    r10cha    r10dis    c10cha    c10dis      
      0.02      0.02   0.00001   0.00001     500.0     500.0
$#    temp   frtherm    r0toth      dudt     tempu     
      25.0         0         1       0.0         0
$# usesocs   tausocs  sicslcid        
         0       0.0         0
*EM_RANDLES_SHORT
2,501
*DEFINE_FUNCTION
501
float resistance_short_randle(float time,
                             float x_sep, float y_sep, float z_sep,
                             float x_sen, float y_sen, float z_sen,
                             float x_ccp, float y_ccp, float z_ccp,
                             float x_ccn, float y_ccn, float z_ccn)
{ 
 if (time > 600.) {
  if (y_sep < 0.12) {
    if (y_sep > 0.08) {
     return 0.002 ;
                      }
    else {
        return -1.;
         }
                    } 
    else {
    return -1. ;
         }
}
}
*DEFINE_CURVE
$uVsSoc
444
0,2.8
10,3.655
100,4.15
200,4.15
*EM_ISOPOTENTIAL
$——–1———2———3———4———5———6———7———8
$#   isoid   settype     setid     
         1         2         1
*EM_ISOPOTENTIAL
$#   isoid   settype     setid     
         2         2         2
*EM_ISOPOTENTIAL
$#   isoid   settype     setid     
         3         2         3
*EM_ISOPOTENTIAL
$#   isoid   settype     setid     
         4         2         4
*EM_ISOPOTENTIAL_CONNECT
$——–1———2———3———4———5———6———7———8
$#   conid   contype    isoid1    isoid2       val      lcid    
         1         2         1         2      0.05        
$*EM_ISOPOTENTIAL_CONNECT
$$#   conid   contype    isoid1    isoid2       val      lcid    
$         2         3         2         0       0.0         
*DEFINE_CURVE
$currentVsTime
555                              
0.,-200
20.,-200
*EM_MAT_001
$——–1———2———3———4———5———6———7———8
$   em_mid     mtype     sigma     eosId          
         1         2    5.88e7
*EM_MAT_001
$——–1———2———3———4———5———6———7———8
$   em_mid     mtype     sigma     eosId          
         2         1
*EM_MAT_001
$——–1———2———3———4———5———6———7———8
$   em_mid     mtype     sigma     eosId          
         3         1
*EM_MAT_001
$——–1———2———3———4———5———6———7———8
$   em_mid     mtype     sigma     eosId          
         4         1
*EM_MAT_001
$——–1———2———3———4———5———6———7———8
$   em_mid     mtype     sigma     eosId          
         5         2     3.7e7
*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
$                                                                     0
*END

 

 

 

 

 

 

 

 

[ErKo]

 

For the third meshless example (https://lsdyna.ansys.com/wp-content/uploads/2022/12/meshless.zip) it should be (em.k file):

$—————————————————————————–
$
$
$ 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                                                           2         2
*EM_CONTROL_CONTACT
1,1,1,1,,,,0.0002
*EM_CONTROL_TIMESTEP
1,&em_dt
*EM_MAT_001
$——–1———2———3———4———5———6———7———8
$   em_mid     mtype     sigma     eosId                              randletype
        11         1      6.e7                                                 1
*EM_MAT_001
$——–1———2———3———4———5———6———7———8
$   em_mid     mtype     sigma     eosId                              randletype
        12         1                                                           2
*EM_MAT_001
$——–1———2———3———4———5———6———7———8
$   em_mid     mtype     sigma     eosId                              randletype
        13         1                                                           3
*EM_MAT_001
$——–1———2———3———4———5———6———7———8
$   em_mid     mtype     sigma     eosId                              randletype
        14         1                                                           4
*EM_MAT_001
$——–1———2———3———4———5———6———7———8
$   em_mid     mtype     sigma     eosId                              randletype
        15         1      3.e7                                                 1
*EM_MAT_001
$——–1———2———3———4———5———6———7———8
$   em_mid     mtype     sigma     eosId                              randletype
         1         2      1.e7
*EM_MAT_001
$——–1———2———3———4———5———6———7———8
$   em_mid     mtype     sigma     eosId                              randletype
         2         2      1.e7
*EM_MAT_001
$——–1———2———3———4———5———6———7———8
$   em_mid     mtype     sigma     eosId                              randletype
         3         4      1.e7
*SET_PART
1
1
*SET_PART
2
2
*SET_PART
3
3
*SET_PART
4
4
*SET_PART
5
5
*SET_PART
6
6
*SET_PART
7
7
*SET_PART
8
8
*SET_PART
9
9
*SET_PART
10
10
*SET_PART
20
11,12,13,14,15,16,17,18
19,20,21,22,23,24,25,26
27,28,29,30
*SET_PART
30
31
*EM_CONTACT
1,1,20,30,,,,0.0002
*EM_RANDLES_MESHLESS
$——–1———2———3———4———5———6———7———8
$ randleId randlType partSetId   rdlArea
         1         0        
$        Q        cQ   SOCinit    SOCtoU
       20.  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_RANDLES_MESHLESS
$——–1———2———3———4———5———6———7———8
$ randleId randlType partSetId   rdlArea
         2         0        
$        Q        cQ   SOCinit    SOCtoU
       20.  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_RANDLES_MESHLESS
$——–1———2———3———4———5———6———7———8
$ randleId randlType partSetId   rdlArea
         3         0  
$        Q        cQ   SOCinit    SOCtoU
       20.  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_RANDLES_MESHLESS
$——–1———2———3———4———5———6———7———8
$ randleId randlType partSetId   rdlArea
         4         0        
$        Q        cQ   SOCinit    SOCtoU
       20.  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_RANDLES_MESHLESS
$——–1———2———3———4———5———6———7———8
$ randleId randlType partSetId   rdlArea
         5         0       
$        Q        cQ   SOCinit    SOCtoU
       20.  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_RANDLES_MESHLESS
$——–1———2———3———4———5———6———7———8
$ randleId randlType partSetId   rdlArea
         6         0       
$        Q        cQ   SOCinit    SOCtoU
       20.  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_RANDLES_MESHLESS
$——–1———2———3———4———5———6———7———8
$ randleId randlType partSetId   rdlArea
         7         0       
$        Q        cQ   SOCinit    SOCtoU
       20.  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_RANDLES_MESHLESS
$——–1———2———3———4———5———6———7———8
$ randleId randlType partSetId   rdlArea
         8         0       
$        Q        cQ   SOCinit    SOCtoU
       20.  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_RANDLES_MESHLESS
$——–1———2———3———4———5———6———7———8
$ randleId randlType partSetId   rdlArea
         9         0        
$        Q        cQ   SOCinit    SOCtoU
       20.  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_RANDLES_MESHLESS
$——–1———2———3———4———5———6———7———8
$ randleId randlType partSetId   rdlArea
        10         0       
$        Q        cQ   SOCinit    SOCtoU
       20.  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         
*EM_ISOPOTENTIAL
$    isoId   setType     setId  randType
         2         2         2        
*EM_ISOPOTENTIAL
$    isoId   setType     setId  randType
         3         2         3         
*EM_ISOPOTENTIAL
$    isoId   setType     setId  randType
         4         2         4         
*EM_ISOPOTENTIAL
$    isoId   setType     setId  randType
         5         2         5         
*EM_ISOPOTENTIAL
$    isoId   setType     setId  randType
         6         2         6         
*EM_ISOPOTENTIAL
$    isoId   setType     setId  randType
         7         2         7         
*EM_ISOPOTENTIAL
$    isoId   setType     setId  randType
         8         2         8         
*EM_ISOPOTENTIAL
$    isoId   setType     setId  randType
         9         2         9         
*EM_ISOPOTENTIAL
$    isoId   setType     setId  randType
        10         2        10         
*EM_ISOPOTENTIAL
$    isoId   setType     setId  randType
        11         2        11         
*EM_ISOPOTENTIAL
$    isoId   setType     setId  randType
        12         2        12         
*EM_ISOPOTENTIAL
$    isoId   setType     setId  randType
        13         2        13         
*EM_ISOPOTENTIAL
$    isoId   setType     setId  randType
        14         2        14         
*EM_ISOPOTENTIAL
$    isoId   setType     setId  randType
        15         2        15         
*EM_ISOPOTENTIAL
$    isoId   setType     setId  randType
        16         2        16         
*EM_ISOPOTENTIAL
$    isoId   setType     setId  randType
        17         2        17         
*EM_ISOPOTENTIAL
$    isoId   setType     setId  randType
        18         2        18         
*EM_ISOPOTENTIAL
$    isoId   setType     setId  randType
        19         2        19         
*EM_ISOPOTENTIAL
$    isoId   setType     setId  randType
        20         2        20        
*EM_ISOPOTENTIAL
$    isoId   setType     setId  randType
        37         2        37
$*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        37                  0.
*EM_ISOPOTENTIAL_CONNECT
$——–1———2———3———4———5———6———7———8
$   connid  connType isoPotId1 isoPotId2     R,V,I      lcid
         2         5         1        11                   1  
*EM_ISOPOTENTIAL_CONNECT
$——–1———2———3———4———5———6———7———8
$   connid  connType isoPotId1 isoPotId2     R,V,I      lcid
         3         5         2        12                   2
*EM_ISOPOTENTIAL_CONNECT
$——–1———2———3———4———5———6———7———8
$   connid  connType isoPotId1 isoPotId2     R,V,I      lcid
         4         5         3        13                   3 
*EM_ISOPOTENTIAL_CONNECT
$——–1———2———3———4———5———6———7———8
$   connid  connType isoPotId1 isoPotId2     R,V,I      lcid
         5         5         4        14                   4 
*EM_ISOPOTENTIAL_CONNECT
$——–1———2———3———4———5———6———7———8
$   connid  connType isoPotId1 isoPotId2     R,V,I      lcid
         6         5         5        15                   5 
*EM_ISOPOTENTIAL_CONNECT
$——–1———2———3———4———5———6———7———8
$   connid  connType isoPotId1 isoPotId2     R,V,I      lcid
         7         5         6        16                   6 
*EM_ISOPOTENTIAL_CONNECT
$——–1———2———3———4———5———6———7———8
$   connid  connType isoPotId1 isoPotId2     R,V,I      lcid
         8         5         7        17                   7 
*EM_ISOPOTENTIAL_CONNECT
$——–1———2———3———4———5———6———7———8
$   connid  connType isoPotId1 isoPotId2     R,V,I      lcid
         9         5         8        18                   8 
*EM_ISOPOTENTIAL_CONNECT
$——–1———2———3———4———5———6———7———8
$   connid  connType isoPotId1 isoPotId2     R,V,I      lcid
        10         5         9        19                   9
*EM_ISOPOTENTIAL_CONNECT
$——–1———2———3———4———5———6———7———8
$   connid  connType isoPotId1 isoPotId2     R,V,I      lcid
        11         5        10        20                  10
$*EM_ISOPOTENTIAL_CONNECT
$$   connid  connType isoPotId1 isoPotId2     R,V,I      lcid
$         2         4        21        22                 
$—+—-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_sen, float y_sen, float z_sen,
                             float x_ccp, float y_ccp, float z_ccp,
                             float x_ccn, float y_ccn, float z_ccn)
{ 
  float distCC;
  float distSEP;
  distCC=sqrt(
     pow((x_ccp-x_ccn),2.)+pow((y_ccp-y_ccn),2.)+pow((z_ccp-z_ccn),2.));
  distSEP=sqrt(
     pow((x_sep-x_sen),2.)+pow((y_sep-y_sen),2.)+pow((z_sep-z_sen),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
         4         4         4         4                                       0
$      mf2       gmv                                            randle
$                   1                                                 0
*END

 

[Changping Yi]

Hello ErKo,

   Thank you so much for the prompt response! basic_socshift and basic_intshort work well now. For Meshless, the error message(error 10307) shows a thermal material ID is missing for many parts. 

*** Error 10307 (KEY+307)
     CHECKING MATERIAL INPUT Part ID= 33
     PART ID 33 has TMID= 0 which does not exist.
     This is PART 24 in the order of input.
     Coupled structural thermal analysis needs to have a thermal material ID defined.

[ErKo]

 

 

 

 

Hi – more than welcome – this works for me (MPP dp R16.1):

yes TMID is not defined – one can use tmid for all parts/composite_shell without one (say tmid=1) – you need also to uncomment the thermal.k include in i.k.

Easier alternative one can have a thermal material card with id TMID=0 and that will be assigned to all parts without any TMID assigned :).
(bit of a trick)

So use in thermal.k (tmid=0 is the last *MAT_THERMAL_ISOTROPIC):

*MAT_THERMAL_ISOTROPIC
$——–1———2———3———4———5———6———7———8
$     tmid       tro     tgrlc    tgmult      tlat      hlat    
         1 7860.0000     0.000     0.000     0.000     0.000
$       hc        tc  
 460.00000 40.000000
*MAT_THERMAL_ISOTROPIC
$——–1———2———3———4———5———6———7———8
$#    tmid       tro     tgrlc    tgmult      tlat      hlat    
         2 7860.0000     0.000     0.000     0.000     0.000
$#      hc        tc  
 460.00000 40.000000
*MAT_THERMAL_ISOTROPIC
$——–1———2———3———4———5———6———7———8
$#    tmid       tro     tgrlc    tgmult      tlat      hlat    
         3 7860.0000     0.000     0.000     0.000     0.000
$#      hc        tc  
 460.00000 40.000000
*MAT_THERMAL_ISOTROPIC
$——–1———2———3———4———5———6———7———8
$#    tmid       tro     tgrlc    tgmult      tlat      hlat    
         0 7860.0000     0.000     0.000     0.000     0.000
$#      hc        tc  
 460.00000 40.000000

All the best

Erko

 

 

 

 

 

[Changping Yi]

 

Hi Erko,

yes, it works if it is downloaded from https://lsdyna.ansys.com/wp-content/uploads/2022/12/meshless.zip. it does not work if  it is downloaded from https://www.dynaexamples.com/em/battery.   After I ran it, I could not find the temperature in LSPP(v9,11, and13) while it is shown in the presentation at 30:20 (

&t=1298s). I wonder if I need to define any term for output.

 

[ErKo]

 

Hi that is great – dynaexamples I can not comment – the one I sent is the official one .

 

for thermal: you need also to uncomment the thermal.k and *include in the i.k.

Then in the thermal.k file add the below instead of what is there now:
so a thermal material card with id TMID=0 and that will be assigned to all parts without any TMID assigned :).
(bit of a trick)

So use in thermal.k (tmid=0 is the last mat with *MAT_THERMAL_ISOTROPIC below):

*KEYWORD
*TITLE
EM: Fall of a sphere on a 10 cells battery
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
$                                                                              $
$                              THERMAL                                         $
$                                                                              $
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
*CONTROL_SOLUTION
$--------1---------2---------3---------4---------5---------6---------7---------8
$     soln
         2
*CONTROL_THERMAL_SOLVER
$--------1---------2---------3---------4---------5---------6---------7---------8
$    atype     ptype    solver     cgtol       gpt    eqheat     fwork       sbc
         1         0        11    1.e-06         8       1.0       1.0       0.0
$   msglvl    maxitr    abstol    reltol     omega                           tsf
         0       500   1.0e-10    1.0e-4       1.0                           1.0
*CONTROL_THERMAL_TIMESTEP
$--------1---------2---------3---------4---------5---------6---------7---------8
$       ts       tip       its      tmin      tmax     dtemp      tscp      lcts
         0        1. &therm_dt &therm_dt &therm_dt        1.       0.5         0
*INITIAL_TEMPERATURE_SET
$--------1---------2---------3---------4---------5---------6---------7---------8
$      nid      temp 
         0       25.
*MAT_THERMAL_ISOTROPIC
$--------1---------2---------3---------4---------5---------6---------7---------8
$     tmid       tro     tgrlc    tgmult      tlat      hlat    
         1 7860.0000     0.000     0.000     0.000     0.000
$       hc        tc  
 460.00000 40.000000
*MAT_THERMAL_ISOTROPIC
$--------1---------2---------3---------4---------5---------6---------7---------8
$#    tmid       tro     tgrlc    tgmult      tlat      hlat    
         2 7860.0000     0.000     0.000     0.000     0.000
$#      hc        tc  
 460.00000 40.000000
*MAT_THERMAL_ISOTROPIC
$--------1---------2---------3---------4---------5---------6---------7---------8
$#    tmid       tro     tgrlc    tgmult      tlat      hlat    
         3 7860.0000     0.000     0.000     0.000     0.000
$#      hc        tc  
 460.00000 40.000000
*MAT_THERMAL_ISOTROPIC
$--------1---------2---------3---------4---------5---------6---------7---------8
$#    tmid       tro     tgrlc    tgmult      tlat      hlat    
         0 7860.0000     0.000     0.000     0.000     0.000
$#      hc        tc  
 460.00000 40.000000
*END


All the best

Erko

 

[Changping Yi]

Hi Erko,

  Thank you for your help! I can get temperature now. However, it does not look like the temperature evolution in the video.  I wonder if there is a short circuit in the model. Battery modelling is new for me. 

[ErKo]

Hello - for sure there is I think some diff. between our model and video :) - me too so not sure why it is.

Erko