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_rigid.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
