*HEADING
Verification of shell and C3D8R element kinematics under large
number of rigid body rotations (50) using SECOND ORDER ACCURACY=YES
on *SECTION CONTROLS -- Rotating Propeller.
**
*NODE
**
**
*ELEMENT, TYPE=C3D8R, ELSET=solid1
*ELEMENT, TYPE=S4R, ELSET=shell1
**
** hub
**
*NSET, NSET=HUB
**
** hub
**
*ELSET, ELSET=HUB
**
** solid
**
*NSET, NSET=SOLID
**
** solid
**
*ELSET, ELSET=SOLID
**
** fem
**
*ELSET, ELSET=FEM
**
** shlfan
**
*NSET, NSET=SHLFAN
**
** shlfan
**
*ELSET, ELSET=SHLFAN
**
** edges
**
*NSET, NSET=EDGES
**
** solfan
**
*NSET, NSET=SOLFAN
**
** solfan
**
*ELSET, ELSET=SOLFAN
**
*MATERIAL,NAME=MAT1
*ELASTIC
3.0e7,0.3
*DENSITY
0.73e-3,
*SECTION CONTROLS, NAME=CONT, SECOND ORDER ACCURACY=YES
*SOLID SECTION,ELSET=SOLID1,MATERIAL=MAT1, CONTROLS=CONT
*SHELL SECTION,ELSET=SHELL1,MATERIAL=MAT1, CONTROLS=CONT
0.5
**
*NODE
*NSET,NSET=NREF
*NSET,NSET=NOUT
** nodes on the outer edges of the blades at its axis.
*ELSET,ELSET=EOUT
** elements on the outer edges of the blades at its axis.
*NSET,NSET=NROOT
*RIGID BODY,ELSET=HUB,REFNODE=90001,TIE NSET=NROOT
*ELSET,ELSET=EREF
**
*RESTART, WRITE, NUM=1
*** 3000 rpm.
*INITIAL CONDITIONS, TYPE=VEL
90001,6,314.16
*INITIAL CONDITIONS, TYPE=ROTATING VEL
SHLFAN,314.16,0.,0.,0.
0.,0.,0.,0.,0.,5.
SOLFAN,314.16,0.,0.,0.
0.,0.,0.,0.,0.,5.
*STEP
*DYNAMIC,EXPLICIT
,1.0
*BOUNDARY,TYPE=VELOCITY
90001,1,5,0.0
*******
*OUTPUT, HISTORY, TIME INTERVAL=0.001
*ENERGY OUTPUT, VAR=ALL, ELSET=SHLFAN
*ENERGY OUTPUT, VAR=ALL, ELSET=SOLFAN
*NODE OUTPUT,NSET=NREF
U,RF,UR3,VR3
*ELEMENT OUTPUT,ELSET=EREF
S,MISES
*NODE OUTPUT,NSET=NOUT
U,UR3,V,VR3, A, AR3
*ELEMENT OUTPUT,ELSET=EOUT
S,MISES
*******
*OUTPUT, FIELD, NUMBER INTERVAL=50, VAR=PRESELECT
*******
**
** Output requests for qa testing
**
*NSET,NSET=QA_TEST
NREF,
*FILE OUTPUT, NUM=9, TIME MARKS=YES
** frequency 9 used to avoid output time points
** matching full/half revolutions.
*NODE FILE,NSET=NREF
U,RF
*ENERGY FILE
*OUTPUT,FIELD,NUMBER INTERVAL=9, TIMEMARKS=YES
*NODE OUTPUT, NSET=QA_TEST
U,RF
*OUTPUT,HISTORY,VAR=PRESELECT,TIME INTERVAL=0.25
*******
*END STEP