Bird
This deck shows calculation of the impact of an idealized bird against a turbine blade. The bird is modeled using SPH method, the turbine is discretized using Lagrangian shell elements. Contact is specified with *CONTACT_AUTOMATIC_NODES_TO_SURFACE keyword.
https://www.dynaexamples.com/sph/bird
https://www.dynaexamples.com/@@site-logo/LS-DYNA-Examples-Logo480x80.png
Bird
This deck shows calculation of the impact of an idealized bird against a turbine blade. The bird is modeled using SPH method, the turbine is discretized using Lagrangian shell elements. Contact is specified with *CONTACT_AUTOMATIC_NODES_TO_SURFACE keyword.
*BOUNDARY_PRESCRIBED_MOTION_RIGID *CONSTRAINED_EXTRA_NODES_SET *CONTACT_AUTOMATIC_NODES_TO_SURFACE_TITLE *CONTROL_CONTACT *CONTROL_ENERGY *CONTROL_HOURGLASS *CONTROL_SHELL *CONTROL_SPH *CONTROL_TERMINATION *CONTROL_TIMESTEP *DATABASE_BINARY_D3PLOT *DATABASE_EXTENT_BINARY *DATABASE_GLSTAT *DATABASE_MATSUM *DATABASE_RBDOUT *DEFINE_BOX *DEFINE_CURVE *ELEMENT_SHELL *ELEMENT_SPH *END *EOS_GRUNEISEN *INITIAL_VELOCITY_NODE *KEYWORD *MAT_NULL *MAT_PIECEWISE_LINEAR_PLASTICITY *MAT_RIGID *NODE *PART *PART_INERTIA *SECTION_SHELL *SECTION_SPH *SET_NODE_LIST *SET_NODE_LIST_GENERATE *SET_PART_LIST *TITLE
*KEYWORD
$* LS-DYNA Keyword file: LSTC LS-PRE/POST Ver 1.0(Beta)- 31OCT2002(1553)
$* Total Number of Nodes = 5185,Largest node ID = 1004160
$* Total Number of Shells = 960,Largest Shell ID = 960
$* Total Number of Parts = 3,Largest Part ID = 101
$* Total Number of Mats = 3 Largest Mat ID = 3
*TITLE
SPH / TP 3 Bird Impact
*CONTROL_CONTACT
$* slsfac rwpnal islchk shlthk penopt thkchg orien enmass
0.0 0.0 0 2 0 1
$* usrstr usrfrc nsbcs interm xpene ssthk ecdt tiedprj
0 0 0 0 0.0
$* sfric dfric edc vfc th th_sf pen_sf
0.0 0.0 0.0 0.0 0.0 0.0 0.0
*CONTROL_ENERGY
$* hgen rwen slnten rylen
2 2 2 1
*CONTROL_HOURGLASS
$* ihq qh
4 0.400000
*CONTROL_SHELL
$* wrpang esort irnxx istupd theory bwc miter proj
20.00000 2 -2 1 11 2 1
*CONTROL_SPH
$
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ *CONTROL_SPH $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$
$* ncbs boxid dt dim memory form start maxv
1 101 0.0 0 0 0 0.0 0.0
*CONTROL_TERMINATION
$* endtim endcyc dtmin endeng endmas
0.001000 0 0.0 0.0 0.0
*CONTROL_TIMESTEP
$* dtinit tssfac isdo tslimt dt2ms lctm erode ms1st
0.0 0.600000 0 0.0 0.0
*DATABASE_GLSTAT
9.9998e-06
*DATABASE_MATSUM
9.9998e-06
*DATABASE_RBDOUT
9.9998e-06
*DATABASE_BINARY_D3PLOT
1.0000e-05
*DATABASE_EXTENT_BINARY
0 0 5
0 0 0 0 0 0 0
*CONTACT_AUTOMATIC_NODES_TO_SURFACE_TITLE
$* cidtitle
1
$* ssid msid sstyp mstyp sboxid mboxid spr mpr
101 1 4 2
$* fs fd dc vc vdc penchk bt dt
0.0 0.0 0.0 0.0 0.0 0 0.0 0.0
$* sfs sfm sst mst sfst sfmt fsf vsf
0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0
$* soft sofscl lcidab maxpar edge depth bsort frcfrq
1 0.100000 0 0.0 0.0 0 3
*SET_NODE_LIST_GENERATE
101 0.0 0.0 0.0 0.0
1000001 1004160
*DEFINE_BOX
101 0.560000 0.750000 -0.220000 -0.068000 -0.049000 0.250000
*DEFINE_CURVE
$
*ELEMENT_SHELL
$
*ELEMENT_SPH
$
*NODE
$
*INITIAL_VELOCITY_NODE
$
*CONSTRAINED_EXTRA_NODES_SET
2 1
*SET_NODE_LIST
1 0.0 0.0 0.0 0.0
1 2 3 4 5 6 7 8
9 10 11 12 13 274 278 281
284 287 290 293 296 299 302 305
308
*BOUNDARY_PRESCRIBED_MOTION_RIGID
$
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$
$$$$ LOADING
$
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$
$* pid dof vad lcid sf vid death birth
2 7 0 1 1.000000 0 0.0 0.0
*PART
$* title
$* pid secid mid eosid hgid grav adpopt tmid
1 1 1
*SECTION_SHELL
$* secid elform shrf nip propt qr/irid icomp setyp
1 11 0.0 5 0 0.0
$* t1 t2 t3 t4 nloc marea
0.003000 0.003000 0.003000 0.003000 0 0.0
*MAT_PIECEWISE_LINEAR_PLASTICITY
$* mid ro e pr sigy etan fail tdel
1 4000.0002.0000e+11 0.300000 1.000000 0.0 0.0 0.0
$* c p lcss lcsr vp
0.0 0.0 2 0 0.0
$* eps1 eps2 eps3 eps4 eps5 eps6 eps7 eps8
0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0
$* es1 es2 es3 es4 es5 es6 es7 es8
0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0
*PART_INERTIA
$* title
$* pid secid mid eosid hgid grav adpopt tmid
2 2 2
$* xc yc zc tm ircs nodeid
0.0 0.0 0.125947 0.130000
$* ixx ixy ixz iyy iyz izz
0.0013833.3715e-041.1187e-079.3520e-04-1.491e-076.9611e-04
$* vtx vty vtz vrx vry vrz
0.0 0.0 0.0 0.0 0.0 -546.6000
*SECTION_SHELL
$* secid elform shrf nip propt qr/irid icomp setyp
2 11 0.0 5 0 0.0
$* t1 t2 t3 t4 nloc marea
0.003000 0.003000 0.003000 0.003000 0 0.0
*MAT_RIGID
$* mid ro e pr n couple malias
2 4000.0002.0000e+11 0.300000 0.0 0.0 0.0
$* cmo con1 con2
1.000000 7.000000 4.000000
$lco or a1 a2 a3 v1 v2 v3
0.0 0.0 0.0 0.0 0.0 0.0
*PART
$
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ SPH PART $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$
$* title
$* pid secid mid eosid hgid grav adpopt tmid
101 101 3 1
*SECTION_SPH
$
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ *SECTION_SPH $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$
$* secid cslh hmin hmax sphini death start
101 1.200000 0.200000 2.000000 0.0 0.0 0.0
*MAT_NULL
$* mid ro pc mu terod cerod ym pr
3 1000.000-1.000e+06 0.001000 0.0 0.0 0.0 0.0
*EOS_GRUNEISEN
$* eosid c s1 s2 s3 gamao a e0
1 1647.000 2.480000 0.0 0.0 0.0 0.0 0.0
$* v0
0.0
*SET_PART_LIST
1 0.0 0.0 0.0 0.0
1
*END