I am trying to simulate a rotating disk consisting of C3D20(R) elements. I have tried both with and without R, the results were not affected. The disk is rotated around the x-axis which is outside of the solid elements.
In a nutshell, here are the control cards of the input deck:
*TIE, NAME = TIE_CYCLIC, CYCLIC SYMMETRY S_CYC_SLAVE,S_CYC_MASTER *CYCLIC SYMMETRY MODEL , N = 50, TIE = TIE_CYCLIC 0.,0.,0.,1.,0.,0. *TRANSFORM,TYPE=C,NSET=N_SUPPRESS 0.,0.,0.,1.,0.,0. *BOUNDARY N_SUPPRESS,2,3 *STEP, NLGEOM = NO *STATIC *DLOAD Disk,CENTRIF,1000000.,0.,0.,0.,1.,0.,0.
You can download my complete model on filetransfer.io
I have cyclic symmetry conditions on the cut faces of the disk so it radially carries loads itself. Axial and circumferential rigid body modes are prevented by prescribing Dirichlet boundary conditions to the respective DOF of the nodal set N_SUPPRESS. I have tried varying N_SUPPRESS, e.g. a single node or multiple adjacent nodes, but the result didn’t change.
The result is unphysical because the smallest displacement is not on the lower end and also the displacements vary at the top of the holes.
I have falsified this agains commercial FE solvers. However, I am not able to fix this input deck. Does anyone know how to fix this?
I am using calculix 2.19 and have tried both spooles and pardiso solver on multiple and on a single core, which didn’t seem to affect anything (except runtime).
Thank you very much!