Pendulum I
A pendulum with two colliding blocks is calculated. The pendulum is modeled with beam elements of elastic behavior. The blocks are modeled with solid elements of elastic behavior and one of it is given an initial velocity. The *DEFORMABLE_TO_RIGID option is used to decrease execution time before impact.
https://www.dynaexamples.com/intro-by-k.-weimar/misc/pendulum-i
https://www.dynaexamples.com/@@site-logo/LS-DYNA-Examples-Logo480x80.png
Pendulum I
A pendulum with two colliding blocks is calculated. The pendulum is modeled with beam elements of elastic behavior. The blocks are modeled with solid elements of elastic behavior and one of it is given an initial velocity. The *DEFORMABLE_TO_RIGID option is used to decrease execution time before impact.
*BOUNDARY_SPC_NODE *CONTACT_SURFACE_TO_SURFACE *CONTROL_TERMINATION *CONTROL_TIMESTEP *DATABASE_BINARY_D3PLOT *DATABASE_EXTENT_BINARY *DATABASE_GLSTAT *DATABASE_HISTORY_NODE *DATABASE_MATSUM *DATABASE_RCFORC *DEFINE_CURVE *DEFORMABLE_TO_RIGID *ELEMENT_BEAM *ELEMENT_SOLID *END *HOURGLASS *KEYWORD *LOAD_BODY_Y *MAT_ELASTIC *NODE *PART *RIGID_DEFORMABLE_D2R *RIGID_DEFORMABLE_R2D *SECTION_BEAM *SECTION_SOLID *SET_SEGMENT *TITLE
*KEYWORD
*TITLE
Rigid/Deformable Material Switching
*CONTROL_TERMINATION
$ f i f f f
$ endtim endcyc dtmin endneg endmas
0.170
*CONTROL_TIMESTEP
$ DTINIT TSSFAC ISDO TSLIMT DT2MS LCTM ERODE MS1ST
0.9
*DEFORMABLE_TO_RIGID
$ i i
$ pid mrb
1
2
$other possible options:
$*RIGID_DEFORMABLE_R2D
$$ i
$$ pid
$ 1
$ 2
$*RIGID_DEFORMABLE_D2R
$$ i
$$ pid
$ 1
$ 2
*MAT_ELASTIC
$ f f f f f f f
$ mid ro e pr da db k
1 7.85e-9 210000.0 0.30
*SECTION_SOLID
$ i i
$ sid elform
1
*SECTION_BEAM
$ i i f f f
$ sid elform shrf qr/irid cst
2 3 1.0 1.0
$ f f f f f f
$ int: ts1 ts2 tt1 tt2 nsloc ntloc
$ f f f f f
$ res: a iss itt irr sa
10.0
$ f f f f f
$ disc:vol iner cid ca offset
$ f f f f f
$ t1 t2 t3 t4 nloc
*PART
$ c
$ head
Solid - Elastic
$ i i i i i i
$ pid sid mid eosid hgid adpopt
1 1 1 0 0 0
*PART
$ c
$ head
Solid - Elastic
$ i i i i i i
$ pid sid mid eosid hgid adpopt
2 1 1 0 0 0
*HOURGLASS
$ i f
$ hgid ihq
1 6
*PART
$ c
$ head
Beam - Elastic
$ i i i i i i
$ pid sid mid eosid hgid adpopt
3 2 1 1
*PART
$ c
$ head
Beam - Elastic
$ i i i i i i
$ pid sid mid eosid hgid adpopt
4 2 1 1
*DATABASE_BINARY_D3PLOT
$ f i
$ dt/cycl lcdt
0.01
*DATABASE_EXTENT_BINARY
$ i i i i i i i i
$ neiph neips maxint strflg sigflg epsflg rltflg engflg
1
$ i i i
$ cmpflg ieverp beamip
*DATABASE_HISTORY_NODE
$ i i i i i i i i
$ id1 id2 id3 id4 id5 id6 id7 id8
1 11
*DATABASE_GLSTAT
$ f i
$ dt binary
0.001
*DATABASE_MATSUM
$ f i
$ dt binary
0.001
*DATABASE_RCFORC
$ f i
$ dt binary
0.001
*LOAD_BODY_Y
$ i f i f f f
$ lcid df lciddr xc yc zc
1 10000.0
*DEFINE_CURVE
$ i i f f f f
$ lcid sidr scla sclo offa offo
1
$ f f
$ abscissa ordinate
0.00 1.000
1000.00 1.000
*CONTACT_SURFACE_TO_SURFACE
$ i i i i i i i i
$ ssid msid sstyp mstyp sboxid mboxid spr mpr
1 2
$ f f f f f i f f
$ fs fd dc vc vdc penchk bt dt
$ f f f f f f f f
$ sfs sfm sst mst sfst sfmt fsf vsf
*SET_SEGMENT
$ i f f f f
$ sid da1 da2 da3 da4
1
$ i i i i f f f f
$ n1 n2 n3 n4 a1 a2 a3 a4
2 6 8 4
*SET_SEGMENT
$ i f f f f
$ sid da1 da2 da3 da4
2
$ i i i i f f f f
$ n1 n2 n3 n4 a1 a2 a3 a4
11 15 17 13
*BOUNDARY_SPC_NODE
$ NID/NSID CID DOFX DOFY DOFZ DOFRX DOFRY DOFRZ
*NODE
$ i8 f f f i i
$ node x y z tc rc
*ELEMENT_SOLID
$ i i i i i i i i i i
$ eid pid n1 n2 n3 n4 n5 n6 n7 n8
*ELEMENT_BEAM
$ i i i i i
$ eid pid n1 n2 n3
*END