Contact I — Welcome to LS-DYNA Examples

Home » Knowledge Base » Simple Examples » Examples Manual » Contact » Contact I

A shell element drops and rebounds on an elastic plate.

Introduction

*CONTACT
Shell Rebounds from Plate Using Five Contact Types
LS-DYNA Manual Section: *CONTACT

Additional Sections:
*INITIAL_VELOCITY

Example: Shell Rebounds from Plate Using Five Contact Types

Filename: contact.plates.k

Description:
A shell element drops and rebounds on an elastic plate.

Model:
The plate measures 40*40*1 mm3 and contains 16 shell elements. The dropped
shell element has a side length of 10 mm, a thickness of 2 mm and drop height
of 10 mm. All shell elements are elastic with Belytschko-Tsay formulation.
The dropped shell element has an initial velocity of 100,000 mm/second
vertically towards the plate. The calculations terminate at 0.0002 seconds.

Input:
All four nodes of the dropped shell element have an initial velocity specified
by *INITIAL_VELOCITY. Contact types 3, 5 and 10 use the dropped shell element
as slave side and the four shell elements in the center of the plate as
master side. The example file has type 3 contact activated, while the other
contact types are commented out. To change contact types, simply comment out
type 3 and un-comment the desired contact.

Type 3 contact is a two way surface to surface algorithm. The segments on the
slave side are checked for penetration of the master segment then the opposite
search takes place.

Type 4 is a single surface algorithm. The nodes of all segments are checked
for penetration of all segments.

Type 5 is a node to surface one way algorithm. The program checks that no slave
node penetrates any master segment.

Type 10 converts surface to surface definition into a node to surface
definition.
Type 13 is a more robust version of the single surface algorithm.

Reference:
Schweizerhof, K. and Weimar, K.

Keywords

$$*CONTACT_AUTOMATIC_SINGLE_SURFACE
$$*CONTACT_NODES_TO_SURFACE
$$*CONTACT_ONE_WAY_SURFACE_TO_SURFACE
$$*CONTACT_SINGLE_SURFACE
$$*SET_NODE_LIST
$$*SET_SEGMENT
*CONTACT_SURFACE_TO_SURFACE
*CONTROL_ENERGY
*CONTROL_HOURGLASS
*CONTROL_TERMINATION
*CONTROL_TIMESTEP
*DATABASE_BINARY_D3PLOT
*DATABASE_BINARY_D3THDT
*DATABASE_EXTENT_BINARY
*DATABASE_GLSTAT
*DATABASE_HISTORY_NODE
*DATABASE_MATSUM
*DATABASE_NCFORC
*DATABASE_NODOUT
*DATABASE_RCFORC
*ELEMENT_SHELL
*END
*INITIAL_VELOCITY
*KEYWORD
*MAT_ELASTIC
*NODE
*PART
*SECTION_SHELL
*SET_NODE_LIST
*SET_SEGMENT
*TITLE

Reduced Input

  *KEYWORD
  *TITLE
  Sliding Interface Types 3,4,5,10,13 
  $
  $  LSTC Example
  $
  $  Last Modified: September 5, 1997
  $
  $  Five different contacts are defined for the same problem.  The only one
  $  active is type 3, surface to surface.  The other four are commented out.
  $  To switch contact types, comment out the active one and remove the comments
  $  from the desired one.  
  $
  $  Units: mm, s
  $
  $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
  $
  $$$$  Control Ouput
  $
  $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
  $
  *CONTROL_TERMINATION
  $   endtim    endcyc     dtmin    endneg    endmas
     .200E-3
  $
  *CONTROL_ENERGY
  $     hgen      rwen    slnten     rylen 
           2         1         2
  $
  *CONTROL_HOURGLASS
  $      ihq        qh
           4
  $
  *CONTROL_TIMESTEP
  $   dtinit      scft      isdo    tslimt      dtms      lctm     erode     ms1st
                  0.50         0
  $
  *DATABASE_BINARY_D3PLOT
  $       dt      lcdt
   0.0100E-3
  $
  *DATABASE_BINARY_D3THDT
  $       dt      lcdt
   2.0000E-3
  $
  *DATABASE_EXTENT_BINARY
  $    neiph     neips    maxint    strflg    sigflg    epsflg    rltflg    engflg
                                         1
  $   cmpflg    ieverp    beamip
  
  $
  *DATABASE_GLSTAT
  $       dt
    0.01e-04
  $
  *DATABASE_NCFORC
  $       dt
    0.01e-04
  $
  *DATABASE_NODOUT
  $       dt
    0.01e-04
  $
  *DATABASE_HISTORY_NODE
  $      id1       id2       id3       id4       id5       id6       id7       id8
          12        13       101
  $
  *DATABASE_MATSUM
  $       dt
    0.10e-05
  $
  *DATABASE_RCFORC
  $       dt
    0.01e-04
  $
  $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
  $
  $$$$  Define Contacts
  $
  $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
  $
  $...>....1....>....2....>....3....>....4....>....5....>....6....>....7....>....8
  $
  $$$$$$$$$$  Type 3, surface to surface
  $
  *CONTACT_SURFACE_TO_SURFACE
  $     ssid      msid     sstyp     mstyp    sboxid    mboxid       spr       mpr
           1         2                                                 1         1
  $       fs        fd        dc        vc       vdc    penchk        bt        dt
  
  $      sfs       sfm       sst       mst      sfst      sfmt       fsf       vsf
  
  $
  *SET_SEGMENT
  $      sid
           1
  $       n1        n2        n3        n4
         101       103       104       102
  $
  *SET_SEGMENT
  $      sid
           2
  $       n1        n2        n3        n4
           7         8        13        12
           8         9        14        13
          12        13        18        17
          13        14        19        18
  $
  $...>....1....>....2....>....3....>....4....>....5....>....6....>....7....>....8
  $
  $$$$$$$$$$  Type 4, single surface
  $                     to make active, remove the $$ from the lines below
  $
  $$*CONTACT_SINGLE_SURFACE
  $     ssid      msid     sstyp     mstyp    sboxid    mboxid       spr       mpr
  $$         1         0                                                 1         1
  $       fs        fd        dc        vc       vdc    penchk        bt        dt
  $$
  $      sfs       sfm       sst       mst      sfst      sfmt       fsf       vsf
  $$
  $
  $$*SET_SEGMENT
  $      sid
  $$         1
  $       n1        n2        n3        n4
  $$       101       103       104       102
  $$         7         8        13        12
  $$         8         9        14        13
  $$        12        13        18        17
  $$        13        14        19        18
  $
  $...>....1....>....2....>....3....>....4....>....5....>....6....>....7....>....8
  $
  $$$$$$$$$$  Type 5, node to surface
  $                     to make active, remove the $$ from the lines below
  $
  $$*CONTACT_NODES_TO_SURFACE
  $     ssid      msid     sstyp     mstyp    sboxid    mboxid       spr       mpr
  $$         1         2         4                                       1         1
  $       fs        fd        dc        vc       vdc    penchk        bt        dt
  $$
  $      sfs       sfm       sst       mst      sfst      sfmt       fsf       vsf
  $$
  $
  $$*SET_NODE_LIST
  $      sid
  $$         1
  $     nid1      nid2      nid3      nid4      nid5      nid6      nid7      nid8
  $$       101       103       104       102
  $
  $$*SET_SEGMENT
  $      sid
  $$         2
  $       n1        n2        n3        n4
  $$         7         8        13        12
  $$         8         9        14        13
  $$        12        13        18        17
  $$        13        14        19        18
  $
  $...>....1....>....2....>....3....>....4....>....5....>....6....>....7....>....8
  $
  $$$$$$$$$$  Type 10, surface to surface
  $                     to make active, remove the $$ from the lines below
  $
  $$*CONTACT_ONE_WAY_SURFACE_TO_SURFACE
  $     ssid      msid     sstyp     mstyp    sboxid    mboxid       spr       mpr
  $$         1         2                                                 1         1
  $       fs        fd        dc        vc       vdc    penchk        bt        dt
  $$
  $      sfs       sfm       sst       mst      sfst      sfmt       fsf       vsf
  $$
  $
  $$*SET_SEGMENT
  $      sid
  $$         1
  $       n1        n2        n3        n4
  $$       101       103       104       102
  $
  $$*SET_SEGMENT
  $      sid       da1       da2       da3       da4
  $$         2
  $       n1        n2        n3        n4
  $$         7         8        13        12
  $$         8         9        14        13
  $$        12        13        18        17
  $$        13        14        19        18
  $
  $...>....1....>....2....>....3....>....4....>....5....>....6....>....7....>....8
  $
  $$$$$$$$$$  Type 13, automatic single surface
  $                     to make active, remove the $$ from the lines below
  $
  $$*CONTACT_AUTOMATIC_SINGLE_SURFACE
  $     ssid      msid     sstyp     mstyp    sboxid    mboxid       spr       mpr
  $$         0         0                                                 1         1
  $       fs        fd        dc        vc       vdc    penchk        bt        dt
  $$
  $      sfs       sfm       sst       mst      sfst      sfmt       fsf       vsf
  $$
  $
  $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
  $
  $$$$  Initial Conditions
  $
  $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
  $
  $$$$  Nodes of the Impactor Material are given an initial velocity.
  $
  *INITIAL_VELOCITY
  $     nsid    nsidex     boxid
           1
  $       vx        vy        vz       
         0.0       0.0 -100000.0
  $
  *SET_NODE_LIST
  $      sid
           1
  $     nid1      nid2      nid3      nid4      nid5      nid6      nid7      nid8
         101       102       103       104
  $
  $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
  $
  $$$$  Define Parts and Materials
  $
  $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
  $
  $...>....1....>....2....>....3....>....4....>....5....>....6....>....7....>....8
  $
  *PART
  $      pid       sid       mid     eosid      hgid    adpopt
  Impacted Material 
           1         1         1
  Impactor Material 
           2         2         1
  $
  $
  *MAT_ELASTIC
  $      mid        ro         e        pr        da        db         k
           1   1.00e-8  100000.0     0.300
  $
  $
  *SECTION_SHELL
  $      sid    elform      shrf       nip     propt   qr/irid     icomp
           1             0.83333       2.0       3.0
  $       t1        t2        t3        t4      nloc
         1.0       1.0       1.0       1.0
  $
  *SECTION_SHELL
  $      sid    elform      shrf       nip     propt   qr/irid     icomp
           2             0.83333       2.0       3.0
  $       t1        t2        t3        t4      nloc
         2.0       2.0       2.0       2.0
  $
  $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
  $
  $$$$  Define Nodes and Elements
  $
  $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
  $
  $$$ Outer edge nodes of the Impacted Material are fixed in translation (tc = 7)
  $
  *END

Figures

image01.gif

image02.gif

Animated Result

Result

Download