Self-Piercing Riveting
Joining aluminum sheets with steel rivets
https://www.dynaexamples.com/introduction/process_simulation/self-piercing-riveting
https://www.dynaexamples.com/@@site-logo/LS-DYNA-Examples-Logo480x80.png
Self-Piercing Riveting
Joining aluminum sheets with steel rivets
*BOUNDARY_PRESCRIBED_MOTION_RIGID_ID *CONTACT_2D_AUTOMATIC_SINGLE_SURFACE_ID *CONTROL_ADAPSTEP *CONTROL_ADAPTIVE *CONTROL_TERMINATION *CONTROL_TIMESTEP *DATABASE_BINARY_D3PLOT *DATABASE_BNDOUT *DATABASE_GLSTAT *DATABASE_RCFORC *DATABASE_SLEOUT *DATABASE_SPCFORC *DEFINE_CURVE_TITLE *ELEMENT_SHELL *END *KEYWORD *MAT_RIGID_TITLE *MAT_PIECEWISE_LINEAR_PLASTICITY_TITLE *NODE *PARAMETER *PART *PART_ADAPTIVE_FAILURE *SECTION_SHELL *SET_PART_LIST_TITLE *TITLE
$-----------------------------------------------------------------------------
$
$ Example provided by Dr. Daniel Hilding (DYNAmore Nordic AB) and Dr. Nils Karajan (DYNAmore)
$
$ E-Mail: info@dynamore.de
$ Web: http://www.dynamore.de
$
$ Copyright, 2015 DYNAmore Nordic AB
$ Copying for non-commercial usage allowed if
$ copy bears this notice completely.
$
$X------------------------------------------------------------------------------
$X
$X 1. Run file as is.
$X Requires LS-DYNA SMP R7.1.2 (or higher) with single precision
$X
$X------------------------------------------------------------------------------
$# UNITS: kg / mm / ms / kN / GPa / kNmm
$X------------------------------------------------------------------------------
$X
*KEYWORD
*TITLE
$# title
$#
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
$ $
$ CONTROL CARDS (GENERAL) $
$ $
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
*PARAMETER
R dtcurv 0.001
R dtplot 0.02
*CONTROL_TERMINATION
$# endtim endcyc dtmin endeng endmas
10.0
*CONTROL_TIMESTEP
$# dtinit tssfac isdo tslimt dt2ms lctm erode ms1st
0.7 -1.0E-5
$# dt2msf dt2mslc imscl
$#
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
$ $
$ CONTROL CARDS (ADAPTIVITY) $
$ $
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
*CONTROL_ADAPSTEP
$# factin dfactr
1.0 0.010
*CONTROL_ADAPTIVE
$# adpfreq adptol adpopt maxlvl tbirth tdeath lcadp ioflag
0.30 0.1 8 3
$# adpsize adpass ireflg adpene adpth memory orient maxel
$# ladpn90 ladpgh ncfred ladpcl adpctl cbirth cdeath lclvl
1
$# cnla
*PART_ADAPTIVE_FAILURE
$# pid t
4 0.10
$#
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
$ $
$ DATABASE CARDS $
$ $
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
*DATABASE_BNDOUT
$# dt binary lcur ioopt
&dtcurv
*DATABASE_RCFORC
$# dt binary lcur ioopt
&dtcurv
*DATABASE_SPCFORC
$# dt binary lcur ioopt
&dtcurv
*DATABASE_SLEOUT
$# dt binary lcur ioopt
&dtcurv
*DATABASE_GLSTAT
$# dt binary lcur ioopt
&dtcurv
*DATABASE_FORMAT
$# iform ibinary
1
*DATABASE_BINARY_D3PLOT
$# dt lcdt beam npltc psetid
&dtplot
$# ioopt
$#
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
$ $
$ CONTACT CARDS $
$ $
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
*CONTACT_2D_AUTOMATIC_SINGLE_SURFACE_ID
$# cid title
1All to all
$# sids sidm sfact freq fs fd dc
1 0.200
$# tbirth tdeath sos som nds ndm cof init
*SET_PART_LIST_TITLE
Contacting parts
$# sid da1 da2 da3 da4 solver
1 MECH
$# pid1 pid2 pid3 pid4 pid5 pid6 pid7 pid8
1 2 3 4 5 6
$#
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
$ $
$ DEFORMABLE PARTS $
$ $
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
*PART
$# title
Steel Rivet
$# pid secid mid eosid hgid grav adpopt tmid
1 1 3 2
*PART
$# title
Top sheet (aluminum)
$# pid secid mid eosid hgid grav adpopt tmid
4 1 4 2
*PART
$# title
Lower sheet (aluminum)
$# pid secid mid eosid hgid grav adpopt tmid
6 1 4 2
*SECTION_SHELL
$# secid elform shrf nip propt qr/irid icomp setyp
1 15 4 1
$# t1 t2 t3 t4 nloc marea idof edgset
1.0
*MAT_PIECEWISE_LINEAR_PLASTICITY_TITLE
Steel
$# mid ro e pr sigy etan fail tdel
3 7.800E-6 210.00 0.300
$# c p lcss lcsr vp
1
$# eps1 eps2 eps3 eps4 eps5 eps6 eps7 eps8
$# es1 es2 es3 es4 es5 es6 es7 es8
*DEFINE_CURVE_TITLE
Steel
$# lcid sidr sfa sfo offa offo dattyp
1
$# a1 o1
0.000 1.500
0.067 1.600
0.100 1.620
*MAT_PIECEWISE_LINEAR_PLASTICITY_TITLE
Aluminum
$# mid ro e pr sigy etan fail tdel
4 2.800E-6 70.000 0.300
$# c p lcss lcsr vp
2
$# eps1 eps2 eps3 eps4 eps5 eps6 eps7 eps8
$# es1 es2 es3 es4 es5 es6 es7 es8
*DEFINE_CURVE_TITLE
Aluminum
$# lcid sidr sfa sfo offa offo dattyp
2
$# a1 o1
0.000 0.073
0.040 0.130
0.100 0.170
0.200 0.200
$#
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
$ $
$ RIGID PARTS $
$ $
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
*PART
$# title
Plunger
$# pid secid mid eosid hgid grav adpopt tmid
2 2 1
*PART
$# title
Blank holder
$# pid secid mid eosid hgid grav adpopt tmid
5 2 1
*PART
$# title
Base tool
$# pid secid mid eosid hgid grav adpopt tmid
3 2 2
*SECTION_SHELL
$# secid elform shrf nip propt qr/irid icomp setyp
2 15 1 1
$# t1 t2 t3 t4 nloc marea idof edgset
1.0
*MAT_RIGID_TITLE
Sliding steel part
$# mid ro e pr n couple m alias
1 7.800E-6 210.00 0.300
$# cmo con1 con2
1.0 6 7
$# lco or a1 a2 a3 v1 v2 v3
*MAT_RIGID_TITLE
Fixed rigid steel part
$# mid ro e pr n couple m alias
2 7.800E-6 210.00 0.300
$# cmo con1 con2
1.0 7 7
$# lco or a1 a2 a3 v1 v2 v3
$#
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
$ $
$ BC CARDS $
$ $
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
*BOUNDARY_PRESCRIBED_MOTION_RIGID_ID
$# id heading
1Plunger
$# pid dof vad lcid sf vid death birth
2 2 2 3 -8.500
$# id heading
2Blank holder
$# pid dof vad lcid sf vid death birth
5 2 2 4 -0.500
*DEFINE_CURVE_TITLE
Tool motion plunger
$# lcid sidr sfa sfo offa offo dattyp
3
$# a1 o1
0.000 0.000
10.000 1.000
100.000 1.000
*DEFINE_CURVE_TITLE
Tool motion holder
$# lcid sidr sfa sfo offa offo dattyp
4
$# a1 o1
0.000 0.000
1.000 1.000
100.000 1.000
*END