Contact Overview
Example showing various contact cases in LS-DYNA
https://www.dynaexamples.com/show-cases/contact-overview
https://www.dynaexamples.com/@@site-logo/LS-DYNA-Examples-Logo480x80.png
Contact Overview
Example showing various contact cases in LS-DYNA
Watch the animation on YouTube
*CONTACT_AUTOMATIC_GENERAL_ID
*CONTACT_AUTOMATIC_SINGLE_SURFACE_ID
*CONTACT_AUTOMATIC_SINGLE_SURFACE_MORTAR_ID
*SECTION_BEAM_TITLE
*SECTION_SHELL_TITLE
*SECTION_SOLID_TITLE
$-----------------------------------------------------------------------------
$
$ Example provided by Dr. Nils Karajan (DYNAmore Corporation)
$
$ Date: Mon Jun 17 10:42:53 CET 2019
$
$ E-Mail: info@dynamore.de
$ Web: http://www.dynamore.de
$
$ Copyright, 2019 DYNAmore Corp.
$ Copying for non-commercial usage allowed if
$ copy bears this notice completely.
$
$-------------------------------------------------------------------------------
$
$ Run main file as is. It was tested with LS-DYNA R9.2
$
$-------------------------------------------------------------------------------
$# Units: kg / mm / ms / kN / GPa / kN-mm
$-------------------------------------------------------------------------------
$
$ Part Numbering Scheme is Defined as:
$
$ | ET | T | CE | N |
$
$ Where ET, T, CE, N are Digits Defined by:
$
$ Element Type 'ET': Beam 10
$ Shell 20
$ Solid 30
$
$ Top/Bottom Flag 'T': Top 1
$ Bottom 0
$
$ Contact Entity 'CE': Beam Node 10
$ Beam 11
$ Shell Corner 20
$ Shell Edge 21
$ Shell Segment 22
$ Solid Corner 30
$ Solid Edge 31
$ Solid Segment 32
$
$ Number 'N': Three-Digit Number to Define Contacting Pairs
$
$ e.g.: Shell Edge to Solid Segment
$
$ Shell Part Number: 20121001
$ Solid Part Number: 30032001
$
*KEYWORD
$
*INCLUDE
$mesh.k
mesh_largerSolids.k
$
*TITLE
$CONTACT_AUTOMATIC_SINGLE_SURFACE
CONTACT_AUTOMATIC_GENERAL
$CONTACT_AUTOMATIC_SINGLE_SURFACE_MORTAR
$
$--- CONTROL
$
*CONTROL_TERMINATION
$# endtim endcyc dtmin endeng endmas nosol
30.0
*CONTROL_TIMESTEP
$# dtinit tssfac isdo tslimt dt2ms lctm erode ms1st
0.9
$# dt2msf dt2mslc imscl unused unused rmscl
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
$
$ CONTROL cards - recommendations for explicit crash applications
$
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
$
*CONTROL_ACCURACY
$# osu inn pidosu iacc
4
*CONTROL_SOLUTION
$# soln nlq isnan lcint lcacc ncdcf
1 1001
*CONTROL_ENERGY
$# hgen rwen slnten rylen
2 2 2 2
*CONTROL_HOURGLASS
$# ihq qh
8 0.1
*CONTROL_BULK_VISCOSITY
$# q1 q2 type btype
1.5 0.06 -1
*CONTROL_CONTACT
$# slsfac rwpnal islchk shlthk penopt thkchg orien enmass
0.1 1.0 2 1
$# usrstr usrfrc nsbcs interm xpene ssthk ecdt tiedprj
$# sfric dfric edc vfc th th_sf pen_sf
$# ignore frceng skiprwg outseg spotstp spotdel spothin
2 1 2 1
$# isym nserod rwgaps rwgdth rwksf icov swradf ithoff
1 1.0
$# shledg pstiff ithcnt tdcnof ftall unused shltrw
*CONTROL_RIGID
$# lmf jntf orthmd partm sparse metalf plotel rbsms
*CONTROL_SHELL
$# wrpang esort irnxx istupd theory bwc miter proj
30.0 1 1 1
$# rotascl intgrd lamsht cstyp6 thshel
$# psstupd sidt4tu cntco itsflg irquad w-mode stretch icrq
$# nfail1 nfail4 psnfail keepcs delfr drcpsid drcprm intperr
1 1.0
*CONTROL_SOLID
$# esort fmatrix niptets swlocl psfail t10jtol icohed tet13k
1 1
$# pm1 pm2 pm3 pm4 pm5 pm6 pm7 pm8 pm9 pm10
*CONTROL_OUTPUT
$# npopt neecho nrefup iaccop opifs ipnint ikedit iflush
1 3 1
$# iprtf ierode tet10s8 msgmax ipcurv gmdt ip1dblt eocs
1 10000
$# tolev newleg frfreq minfo solsig msgflg cdetol
*CONTROL_MPP_IO_NODUMP
$
$
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
$
$ CONTACT cards
$
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
$
$*CONTACT_AUTOMATIC_SINGLE_SURFACE_ID
*CONTACT_AUTOMATIC_GENERAL_ID
$*CONTACT_AUTOMATIC_GENERAL_INTERIOR_ID
$*CONTACT_AUTOMATIC_GENERAL_EDGEONLY_ID
$*CONTACT_AUTOMATIC_SINGLE_SURFACE_MORTAR_ID
$# cid title
1
$# ONLY _MPP
$# ignore bckt lcbckt ns2trk inititr parmax cparm8
$ 2
$# ssid msid sstyp mstyp sboxid mboxid spr mpr
42 2
$# fs fd dc vc vdc penchk bt dt
0.2 20.0
$# sfs sfm sst mst sfst sfmt fsf vsf
$# soft sofscl lcidab maxpar sbopt depth bsort frcfrq
1 2
$# penmax thkopt shlthk snlog isym i2d3d sldthk sldstf
$# igap ignore dprfac dtstif unused unused flangl cid_rcf
2
$# q2tri dtpchk sfnbr fnlscl dnlscl tcso tiedid shledg
$# sharec cparm8 ipback srnde fricsf icor ftorq region
$# pstiff ignroff fstol 2dbinr
1
*SET_PART_LIST_GENERATE
$# sid da1 da2 da3 da4 solver
42
$# b1beg b1end b2beg b2end b3beg b3end b4beg b4end
10000000 30199999
$
$
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
$
$ BOUNDARY / INITIAL CONDITIONS
$
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
$
*BOUNDARY_SPC_SET_ID
$# id heading
1 Fix lower elements
$# nsid cid dofx dofy dofz dofrx dofry dofrz
100200300 1 1 1 1 1 1
*SET_NODE_ADD
$# sid da1 da2 da3 da4 solver
100200300 MECH
$# nsid1 nsid2 nsid3 nsid4 nsid5 nsid6 nsid7 nsid8
100200301 100200302
*SET_NODE_LIST_GENERATE
$# sid da1 da2 da3 da4 solver
100200301 MECH
$# b1beg b1end b2beg b2end b3beg b3end b4beg b4end
1000000 1009999 2000000 2009999 3000000 3009999
*SET_NODE_GENERAL
$# sid da1 da2 da3 da4 solver
100200302 MECH
$# option e1 e2 e3 e4 e5 e6 e7
PART 10110000 10111000 20120000
PART 20120000 20121000 20122000
PART 30130000 30131000 30132000
$
$---
$
*INITIAL_VELOCITY_GENERATION
$#nsid/pid styp omega vx vy vz ivatn icid
101201301 1 -0.2
$# xc yc zc nx ny nz phase irigid
*SET_PART_LIST_GENERATE
$# sid da1 da2 da3 da4 solver
101201301 MECH
$# b1beg b1end b2beg b2end b3beg b3end b4beg b4end
10100000 10199999 20100000 20199999 30100000 30199999
$
$
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
$
$ DATABASE
$
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
$
*DATABASE_GLSTAT
$# dt binary lcur ioopt
0.25 1
*DATABASE_MATSUM
$# dt binary lcur ioopt
0.25 1
*DATABASE_RCFORC
$# dt binary lcur ioopt
0.25 1
*DATABASE_BINARY_D3PLOT
$# dt lcdt beam npltc psetid
0.25
$# ioopt
$
$
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
$
$ SECTIONS
$
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
$
*SECTION_BEAM_TITLE
Beams
$# secid elform shrf qr/irid cst scoor nsm
1 1 0.833 1
$# ts1 ts2 tt1 tt2 nsloc ntloc
0.5 0.5
*SECTION_BEAM_TITLE
Spotweld Beams
$# secid elform shrf qr/irid cst scoor nsm
11 9 1.0 2 1
$# ts1 ts2 tt1 tt2 print
0.5 0.5
*SECTION_SHELL_TITLE
Shells
$# secid elform shrf nip propt qr/irid icomp setyp
2 16 0.833 5
$# t1 t2 t3 t4 nloc marea idof edgset
0.5
*SECTION_SOLID_TITLE
Solids
$# secid elform aet
3 -1
$
$
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
$
$ MATERIALS
$
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
$
*MAT_PLASTIC_KINEMATIC_TITLE
Mild Steel
$# mid ro e pr sigy etan beta
1 7.830E-6 210.0 0.3 0.3 50.0
$# src srp fs vp
*MAT_SPOTWELD_TITLE
Mild Steel
$# mid ro e pr sigy eh dt tfail
11 7.830e-6 210.0 0.3 0.3 50.0
$# efail nrr nrs nrt mrr mss mtt nf
*MAT_NULL
$# mid ro pc mu terod cerod ym pr
111 7.830e-6 210.0 0.3
*END
All Contact Scenarios
Individual Contact Scenario
