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 Iñaki (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 MPP Dev 118745  (or higher) with double precision 
$X
$X------------------------------------------------------------------------------
$# UNITS: (S.I) 
$X------------------------------------------------------------------------------
$X
*KEYWORD
*TITLE
EM Resistive spot welding problem (RSW)
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
$                                                                              $
$                             PARAMETERS                                       $
$                                                                              $
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
*PARAMETER
R    T_end       0.2
R  dt_plot     0.002
$
$--- EM
$
R    em_dt     0.001
Iem_bemmtx        10
Iem_femmtx        10
R  em_cond      1.e6
R  em_curr      845. 
$--- STRUC
$
R struc_dt     0.001
R  rho_el1  8.89e+03
R  rho_el2  8.89e+03
R  rho_pl1    8.02e3
R  rho_pl2    7.87e3
R  struc_E   6.89e10
R struc_nu      0.33
$--- THERMAL
$
R therm_dt     0.001
R therm_hc       480
R therm_tc        52 
R therm_ti       25.
*PARAMETER_EXPRESSION
$#    prmr                                                            expression
rtr       293.15
rt0       273.15
rti       293.15
rdt0      0.001
rtdt0     0.01
*INCLUDE
mesh.k
*INCLUDE
struc.k
*include
thermal.k 
*DATABASE_BINARY_D3PLOT
&dt_plot
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
$                                                                              $
$                             EM CONTROL                                       $
$                                                                              $
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
*EM_CONTROL
$    emsol     numls  emdtinit   emdtmax   emtinit    emtend  ncyclFem  ncyclBem
         3                                                  &em_femmtx&em_bemmtx 
*EM_CONTROL_CONTACT
         1         1
*EM_CONTROL_TIMESTEP
$   tstype    dtcons      lcid    factor
         1    &em_dt
*EM_OUTPUT
$     matS      matF      solS      solF      mesh    memory    timing d3plotAsc
         2         2         2         2         0         0         0         0
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
$                                                                              $
$                             EM CONTACT                                       $
$                                                                              $
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
*EM_CONTACT
$   CONTID    COTYPE     PSIDM     PSIDS
         1         1     10001     10002       0.3       0.3                5e-4
*EM_CONTACT
$   CONTID    COTYPE     PSIDM     PSIDS
         2         1     10002     10003       0.3       0.3                5e-4
*EM_CONTACT
$   CONTID    COTYPE     PSIDM     PSIDS
         3         1     10003     10004       0.3       0.3                5e-4
*EM_CONTACT_RESISTANCE
$#    crid    contid     ctype    circid   jhrtype       
         1         1         1         0         1
$#    lcid    
       900
*EM_CONTACT_RESISTANCE
$#    crid    contid     ctype    circid   jhrtype       
         2         2         1         0         1
$#    lcid    
       900
*EM_CONTACT_RESISTANCE
$#    crid    contid     ctype    circid   jhrtype       
         3         3         1         0         1
$#    lcid    
       900
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
$                                                                              $
$                             EM ISOPOTENTIAL                                  $
$                                                                              $
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
*EM_ISOPOTENTIAL
$#   isoid   settype     setid     
         1         2         1
*EM_ISOPOTENTIAL
$#   isoid   settype     setid     
         2         2         2
*EM_ISOPOTENTIAL_CONNECT
$#   conid   contype    isoid1    isoid2       val      lcid    
         1         4         1         2       0.0         1
*EM_ISOPOTENTIAL_CONNECT
$#   conid   contype    isoid1    isoid2       val      lcid    
         2         3         2         0       0.0         0
*DEFINE_CURVE_TITLE
Load-Time-Function for Current 
$#    lcid      sidr       sfa       sfo      offa      offo    dattyp     lcint
         1         0       1.0  &em_curr       0.0       0.0         0         0
$#                a1                  o1  
                 0.0                 1.0
            900000.0                 1.0
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
$                                                                              $
$                             EM MAT and EOS                                   $
$                                                                              $
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
*EM_MAT_001
$   em_mid     mtype     sigma       eos              
         4         2                   1
*EM_MAT_001
$   em_mid     mtype     sigma       eos              
         7         2                   1    
$
*EM_MAT_001
$   em_mid     mtype     sigma       eos              
         5         2                   1
*EM_MAT_001
$   em_mid     mtype     sigma       eos              
         6         2                   1
*EM_EOS_TABULATED1
$#   EOSID      LCID
         1         2
*DEFINE_CURVE_TITLE
Electrical Conductivity function of temperature
$#    LCID      SIDR       SFA       SFO      OFFA      OFFO    DATTYP
         2         0       1.0  &em_cond       0.0       0.0         0
0.0,1.
25,1.
100,0.1
200,0.05
2000,0.05  
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
$                                                                              $
$                             DEFINE_FUNCTIONS                                 $
$                                                                              $
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
*DEFINE_FUNCTION
$#     fid                                                               heading
       901                                                                      
$#                                                                      function
float r901(float time,float area_glo,float area_loc,float dist_ct,
float curr_loc,float curr_glo,float JHrate_loc,
float ctpress_mst,float ctpress_slv,float temp_mst,float temp_slv,
float cond_mst,float cond_slv,float mass_mst,float mass_slv,
float vmstress_mst,float vmstress_slv)
{
float f,t;
f=1./(ctpress_mst+ctpress_slv);  
return f;
}
*DEFINE_FUNCTION
$#     fid                                                               heading
       900                                                                      
$#                                                                      function
float r900(float time,float area_glo,float area_loc,float dist_ct,
float curr_loc,float curr_glo,float JHrate_loc,
float ctpress_mst,float ctpress_slv,float temp_mst,float temp_slv,
float cond_mst,float cond_slv,float mass_mst,float mass_slv,
float vmstress_mst,float vmstress_slv)
{
float f,t;
f=0.000001;  
return f;
}
*DEFINE_FUNCTION
$#     fid                                                               heading
       902                                                                      
$#                                                                      function
float r902(float time,float area_glo,float area_loc,float dist_ct,
float curr_loc,float curr_glo,float JHrate_loc,
float ctpress_mst,float ctpress_slv,float temp_mst,float temp_slv,
float cond_mst,float cond_slv,float mass_mst,float mass_slv,
float vmstress_mst,float vmstress_slv)
{
float f,t;
if (ctpress_slv < 0.0) {
ctpress_slv = 0.0;
}
t=(temp_mst+temp_slv)/2;
f=1.0/area_glo*0.56*
(((ctpress_slv)-(-5))/(20-(-5)))**-0.3333333333*
((t-1000+(293.15-t)*2**(-1/-0.3))
/(293.15-1000))**-0.3;
if (f <= 0.0) {
f = 1.0E-12;
}
return f;
}
*DEFINE_FUNCTION
$#     fid                                                               heading
       903                                                                      
$#                                                                      function
float r903(float time,float area_glo,float area_loc,float dist_ct,
float curr_loc,float curr_glo,float JHrate_loc,
float ctpress_mst,float ctpress_slv,float temp_mst,float temp_slv,
float cond_mst,float cond_slv,float mass_mst,float mass_slv,
float vmstress_mst,float vmstress_slv)
{
float f,t;
if (ctpress_slv < 0.0) {
if (ctpress_mst < 0.0) {
f = 1.0e-12;
return f;
}
}
f=3*vmstress_slv/(ctpress_mst/2+ctpress_slv/2)*
((1/cond_mst)+(1/cond_slv))*0.5;
if (f <= 0.0) {
f = 1.0E-12;
}
return f;
}
*END