Eroding

A projectile strikes a plate at a critical angle.

*CONTACT_ERODING_SURFACE_TO_SURFACE
Projectile Penetrates Plate
LS-DYNA Manual Section: *CONTACT_ERODING_SURFACE_TO_SURFACE

Additional Sections:
*INITIAL_VELOCITY_GENERATION

Example: Projectile Penetrates Plate

Filename: contact.projectile.k

Description:
A projectile strikes a plate at a critical angle.

Model:
The hemispherical projectile has a length of 7.67 cm and a diameter of
0.767 cm. The plate measures 23.01 cm * 23 cm * 0.64 cm. The projectile and
the plate are elastic perfectly plastic with failure strain. The initial
velocity of the projectile is 0.129 cm/ microsecond at an angle of 75 degrees.
The calculation terminates at 110.0 usec.

Input:
The initial velocity (magnitude and direction) of the projectile is set using
*INITIAL_VELOCITY_GENERATION. Eroding contact between the projectile surface
and plate surface is defined so that the contact erodes as the element erodes
(*CONTACT_ERODING_SURFACE_TO_SURFACE). This allows the contact to work
correctly as layers of the parts erode during penetration.

Results:
The projectile fractures into a tip and trailing portion. The trailing portion
punches a hole through the plate while the tip deflects off the plate.

*CONTACT_ERODING_SURFACE_TO_SURFACE
*CONTROL_CONTACT
*CONTROL_ENERGY
*CONTROL_OUTPUT
*CONTROL_TERMINATION
*DATABASE_BINARY_D3PLOT
*DATABASE_BINARY_D3THDT
*DATABASE_EXTENT_BINARY
*DATABASE_GLSTAT
*DATABASE_MATSUM
*DATABASE_RCFORC
*DATABASE_SLEOUT
*ELEMENT_SOLID
*END
*INITIAL_VELOCITY_GENERATION
*KEYWORD
*MAT_PLASTIC_KINEMATIC
*NODE
*PART
*SECTION_SOLID
*SET_PART
*TITLE

  *KEYWORD
  *TITLE
  Projectile Penetrating Plate
  $
  $  LSTC Example
  $
  $  Last Modified: September 8, 1997
  $
  $   Units: gram, cm, microsec, 1e+07 N, Mbar, 1e+07 N-cm
  $
  $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
  $
  $$$$  Control Ouput
  $
  $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
  $
  $...>....1....>....2....>....3....>....4....>....5....>....6....>....7....>....8
  $
  *CONTROL_TERMINATION
  $   endtim    endcyc     dtmin    endneg    endmas
   1.100E+02
  $
  *CONTROL_CONTACT
  $   slsfac    rwpnal    islchk    shlthk    penopt    thkchg     orien
         1.0
  $   usrstr    usrfac     nsbcs    interm    xpenen
  
  $
  *CONTROL_ENERGY
  $     hgen      rwen    slnten     rylen
           2         2
  $
  *CONTROL_OUTPUT
  $    npopt    neecho    nrefup    iaccop     opifs    ipnint    ikedit
           1         3
  $
  $
  *DATABASE_BINARY_D3PLOT
  $       dt      lcdt
   10.000000
  $
  *DATABASE_EXTENT_BINARY
  $    neiph     neips    maxint    strflg    sigflg    epsflg    rltflg    engflg
  
  $          ieverp = 1 put each plot state in separate d3plot files
  $   cmpflg    ieverp    beamip
                     1
  $
  *DATABASE_BINARY_D3THDT
  $       dt      lcdt
      999999
  $
  *DATABASE_GLSTAT
  $       dt
        0.10
  $
  *DATABASE_MATSUM
  $       dt
        0.10
  $
  *DATABASE_RCFORC
  $       dt
        0.10
  $
  *DATABASE_SLEOUT
  $       dt
        0.10
  $
  $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
  $
  $$$$  Define Contacts
  $
  $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
  $
  $...>....1....>....2....>....3....>....4....>....5....>....6....>....7....>....8
  $
  *CONTACT_ERODING_SURFACE_TO_SURFACE
  $     ssid      msid     sstyp     mstyp    sboxid    mboxid       spr       mpr
           1         2         3         3
  $
  $       fs        fd        dc        vc       vdc    penchk        bt        dt
  
  $
  $      sfs       sfm       sst       mst      sfst      sfmt       fsf       vsf
  
  $
  $     isym    erosop      iadj
           1         1
  $
  $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
  $
  $$$$  Initial Conditions
  $
  $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
  $
  $$$$  Assign an initial velocity to the projectile (part 1) angled down
  $$$$  towards the plate.
  $
  $...>....1....>....2....>....3....>....4....>....5....>....6....>....7....>....8
  $
  *INITIAL_VELOCITY_GENERATION
  $      sid      styp     omega        vx        vy        vz
           1         1           1.246E-01 0.000E+00-3.339E-02
  $
  $       xc        yc        zc        nx        ny        nz     phase
   
  $
  *SET_PART
           1
           1
  $
  $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
  $
  $$$$  Define Parts and Materials
  $
  $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
  $
  $...>....1....>....2....>....3....>....4....>....5....>....6....>....7....>....8
  $
  *PART
  $      pid       sid       mid     eosid      hgid    adpopt
  Projectile            
           1         1         1
  Plate
           2         1         2
  $
  $
  $$$$$$  Materials
  $
  $$$$  failure strain for erosion of the projectile and plate elements are
  $$$$  set as:  fs = 0.8 
  $
  *MAT_PLASTIC_KINEMATIC
  $      mid        ro         e        pr      sigy      etan      beta
           1 1.862E+01 1.170E+00      0.22 1.790E-02                 1.0
  $      src       srp        fs
                             0.8
  $
  $
  *MAT_PLASTIC_KINEMATIC
  $      mid        ro         e        pr      sigy      etan      beta
           2 7.896E+00 2.100E+00     0.284 1.000E-02                 1.0
  $      src       srp        fs
                             0.8
  $
  $
  $$$$$$  Sections
  $
  *SECTION_SOLID
  $      sid    elform
           1         0
  $
  $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
  $
  $$$$  Define Nodes and Elements
  $
  $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
  $
  *END