x

# Reduced Input

```\$-----------------------------------------------------------------------------
\$
\$ Example provided by Iñaki (LSTC)
\$
\$ E-Mail: info@dynamore.de
\$ Web: http://www.dynamore.de
\$
\$ 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
\$#    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
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
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
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
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```