Defo2Rigid

Two spheres are connected to wires to form two pendulums. One sphere is in a horizontal position with gravitational acceleration, base acceleration and is given an initial velocity in the vertical direction. The other sphere is in the vertical direction. The spheres are treated as rigid bodies while no contact or deformation occurs (i.e., when the horizontal pendulum swings down towards the vertical pendulum). The spheres are switched to deformable through a restart file so that they become flexible during contact.

*DEFORMABLE_TO_RIGID
Interaction of Pendulums
LS-DYNA Manual Section: *DEFORMABLE_TO_RIGID

Additional Sections:
*BOUNDARY_SPC_NODE
*LOAD_BODY_Y *RIGID_DEFORMABLE_R2D

Example: Interaction of Pendulums

Filenames: deformable_to_ri gid.pendulum.k
deformable_to_rigid.pendulum.res

Execution lines:
ls940 i= deformable_to_rigid.pendulum.k
ls940 i= deformable_to_rigid.pendulum.res r=d3dump01

Description:
Two spheres are connected to wires to form two pendulums. One sphere is in a
horizontal position with gravitational acceleration, base acceleration and is
given an initial velocity in the vertical direction. The other sphere is in
the vertical direction. The spheres are treated as rigid bodies while no
contact or deformation occurs (i.e., when the horizontal pendulum swings down
towards the vertical pendulum). The spheres are switched to deformable through
a restart file so that they become flexible during contact.

Model:
Both spheres are modeled using shell elements. The pendulum wires are modeled
using elastic beams. Automatic single surface contact is used during the impact
phase.

Reference:
Reid, J.D.

deformable_to_rigid.pendulum.k:$ - uses *DEFORMABLE_TO_RIGID option to decrease execution time before impact
deformable_to_rigid.pendulum.k:*BOUNDARY_SPC_NODE
deformable_to_rigid.pendulum.k:*CONTACT_AUTOMATIC_SINGLE_SURFACE
deformable_to_rigid.pendulum.k:*CONTROL_CONTACT
deformable_to_rigid.pendulum.k:*CONTROL_ENERGY
deformable_to_rigid.pendulum.k:*CONTROL_OUTPUT
deformable_to_rigid.pendulum.k:*CONTROL_SHELL
deformable_to_rigid.pendulum.k:*CONTROL_TERMINATION
deformable_to_rigid.pendulum.k:*DATABASE_BINARY_D3PLOT
deformable_to_rigid.pendulum.k:*DATABASE_BINARY_D3THDT
deformable_to_rigid.pendulum.k:*DATABASE_EXTENT_BINARY
deformable_to_rigid.pendulum.k:*DATABASE_GLSTAT
deformable_to_rigid.pendulum.k:*DATABASE_HISTORY_NODE
deformable_to_rigid.pendulum.k:*DATABASE_MATSUM
deformable_to_rigid.pendulum.k:*DATABASE_NODOUT
deformable_to_rigid.pendulum.k:*DATABASE_RBDOUT
deformable_to_rigid.pendulum.k:*DEFINE_BOX
deformable_to_rigid.pendulum.k:*DEFINE_CURVE
deformable_to_rigid.pendulum.k:*DEFORMABLE_TO_RIGID
deformable_to_rigid.pendulum.k:*ELEMENT_BEAM
deformable_to_rigid.pendulum.k:*ELEMENT_SHELL
deformable_to_rigid.pendulum.k:*END
deformable_to_rigid.pendulum.k:*INITIAL_VELOCITY
deformable_to_rigid.pendulum.k:*KEYWORD
deformable_to_rigid.pendulum.k:*LOAD_BODY_Y
deformable_to_rigid.pendulum.k:*MAT_ELASTIC
deformable_to_rigid.pendulum.k:*MAT_PLASTIC_KINEMATIC
deformable_to_rigid.pendulum.k:*NODE
deformable_to_rigid.pendulum.k:*PART
deformable_to_rigid.pendulum.k:*SECTION_BEAM
deformable_to_rigid.pendulum.k:*SECTION_SHELL
deformable_to_rigid.pendulum.k:*TITLE
deformable_to_rigid.pendulum.res:*CONTROL_TERMINATION
deformable_to_rigid.pendulum.res:*DATABASE_BINARY_D3PLOT
deformable_to_rigid.pendulum.res:*END
deformable_to_rigid.pendulum.res:*KEYWORD
deformable_to_rigid.pendulum.res:*RIGID_DEFORMABLE_R2D
deformable_to_rigid.pendulum.res:*TITLE

  *KEYWORD
  *TITLE
  Pendulum with 2 spheres colliding
  $
  $  LSTC Example
  $
  $  - uses *DEFORMABLE_TO_RIGID option to decrease execution time before impact
  $
  $  - one sphere is given an initial velocity (gravity alone just takes
  $      too long for the pendulum to swing)
  $
  $
  $  Last Modified: September 16, 1997
  $
  $  Units: mm, kg, ms, kN, GPa, kN-mm
  $
  $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
  $
  $$$$  Control Ouput
  $
  $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
  $
  $...>....1....>....2....>....3....>....4....>....5....>....6....>....7....>....8
  $
  *CONTROL_TERMINATION 
  $   endtim    endcyc     dtmin    endeng    endmas
        11.0         0       0.0       0.0       0.0
  $
  *CONTROL_CONTACT
  $   slsfac    rwpnal    islchk    shlchk    penopt    thkchg     orien
                                         2   
  $   usrstr    usrfrc     nsbcs    interm     xpene
  
  $
  *CONTROL_ENERGY
  $     hgen      rwen    slnten     rylen
           2         2
  $
  *CONTROL_OUTPUT
  $    npopt    neecho    nrefup    iaccop     opifs    ipnint    ikedit
           1         3
  $
  *CONTROL_SHELL
  $   wrpang    itrist     irnxx    istupd    theory       bwc     miter
                                         1         2
  $
  $
  *DATABASE_BINARY_D3PLOT
  $       dt      lcdt
        1.00
  $
  *DATABASE_EXTENT_BINARY
  $    neiph     neips    maxint    strflg    sigflg    epsflg    rltflg    engflg
  
  $   cmpflg    ieverp    beamip
                     1
  $
  *DATABASE_BINARY_D3THDT
  $       dt      lcdt
      999999
  $
  $
  *DATABASE_GLSTAT
  $       dt
        0.10
  $
  *DATABASE_MATSUM
  $       dt
        0.10
  $      
  *DATABASE_NODOUT
  $       dt
        0.10
  $
  *DATABASE_HISTORY_NODE
  $      id1       id2       id3
         350       374       678       713
  $
  *DATABASE_RBDOUT
  $       dt
        0.10
  $
  $
  $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
  $
  $$$$  Define Contacts - Sliding Interfaces
  $
  $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
  $
  $...>....1....>....2....>....3....>....4....>....5....>....6....>....7....>....8
  $
  *CONTACT_AUTOMATIC_SINGLE_SURFACE
  $     ssid      msid     sstyp     mstyp
           0
  $  Equating ssid to zero means that all segments are included
  $
  $       fs        fd
        0.08      0.08
  
  $
  $
  $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
  $
  $$$$  Gravity
  $
  $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
  $
  *LOAD_BODY_Y
  $     lcid        sf    lciddr        xc        yc        zc
           1   0.00981
  $
  *DEFINE_CURVE
  $     lcid      sidr      scla      sclo      offa      offo
           1
  $           abscissa            ordinate
  $
                  0.00               1.000
              10000.00               1.000
  $
  $
  $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
  $
  $$$$  Boundary and Initial Conditions
  $
  $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
  $
  $...>....1....>....2....>....3....>....4....>....5....>....6....>....7....>....8
  $
  $$$$  Constrain translation of end points of beams
  $ 
  *BOUNDARY_SPC_NODE
  $      nid       cid      dofx      dofy      dofz     dofrx     dofry     dofrz
       45004         0         1         1         1         0         0         0
       45005         0         1         1         1         0         0         0
       45010         0         1         1         1         0         0         0
       45011         0         1         1         1         0         0         0
  $
  $
  $$$$  The nodes within box 5 are given an initial velocity.
  $
  *INITIAL_VELOCITY
  $     nsid    nsidex     boxid
                               5
  $
  $       vx        vy        vz        wx        wy        wz
         0.0     -12.0       0.0
  $
  *DEFINE_BOX
  $    boxid       xmm       xmx       ymn       ymx       zmn       zmx
           5    -120.0     -80.0      80.0     120.0     -30.0      30.0
  $
  $
  $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
  $
  $$$$  Define Parts and Materials
  $
  $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
  $ 
  $...>....1....>....2....>....3....>....4....>....5....>....6....>....7....>....8
  $
  $$$$$   SPHERES
  $
  *PART
  $      pid       sid       mid     eosid      hgid    adpopt
  sphere1
           1         1         1
  sphere2
           2         2         1
  $
  $ 
  $$$$ Materials
  $
  $  Aluminum
  $
  *MAT_PLASTIC_KINEMATIC
  $      mid        ro         e        pr      sigy      etan      beta
           1   2.70e-6      68.9     0.330     0.286   0.00689
  $      src       srp        fs
  
  $ 
  $$$$ Sections
  $
  $
  *SECTION_SHELL
  $ 
  $      sid    elform      shrf       nip     propt   qr/irid     icomp
           1         2
  $       t1        t2        t3        t4      nloc
         1.0       1.0       1.0       1.0
  $
  *SECTION_SHELL
  $ 
  $      sid    elform      shrf       nip     propt   qr/irid     icomp
           2         2
  $       t1        t2        t3        t4      nloc
         1.0       1.0       1.0       1.0
  $
  $
  $$$$$   PENDULUM WIRES - ELASTIC BEAMS
  $
  *PART
  $      pid       sid       mid     eosid      hgid    adpopt
  Pendulum Wires - Elastic Beams
          45        45        45
  $
  $
  *MAT_ELASTIC
  $      mid        ro         e        pr        da        db         k
          45   7.86e-6     210.0      0.30
  $
  $
  *SECTION_BEAM
  $      sid    elform      shrf   qr/irid       cst
          45         3   1.00000       1.0
  $ res:   a       iss       itt       irr        sa
        10.0
  $
  $
  $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
  $
  $$$$  Deformable Switching
  $
  $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
  $
  $...>....1....>....2....>....3....>....4....>....5....>....6....>....7....>....8
  $
  *DEFORMABLE_TO_RIGID
  $      pid       mrb
           1
  $
  *DEFORMABLE_TO_RIGID
  $      pid       mrb
           2
  $
  $
  $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
  $
  $$$$  Define Nodes
  $
  $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
  $
  $...>....1....>....2....>....3....>....4....>....5....>....6....>....7....>....8
  $
  $ 
  *END