Ansys hyperelastic + viscoelastic material trouble
I'm having trouble when using the combined hyperelastic material (neo-hookean) with viscoelastic (one term prony) in Ansys. The problem is a cone being pulled by uniform pressure at the top and fixed at the
bottom. Plane183 (8 node plane element with uniform reduced integration) is used with axisymmetric option. This element empolyed mixed u-p formation. The hyperelastic material is assumed to be incompressible (u=0.5) or nearly incompressible (u=0.4995). (actually if use u=0.4995, Ansys will not run at all, the element will not form at the first substep)
The model runs fine if only the neo-hookean hyperelastic
material (u=0.5) is used but will crash at about 30% of completion when viscoelastic material is combined. By looking at the hydrostatic pressure distribution in post-processing, there is obvious checkerboard pattern along the loading edge, so it seems that volumetric locking causes the analysis to stop.
My question is why this happens when viscoelastic material is combined but not happens when only hyperelastic material is used.
My guess is our version of Ansys is too limited (university intermediate version)?
So if anyone in this forum has higher version Ansys in his/her hand, can he/she kindly run this analysis for me and see if it's due to the version problem?
The attached input file ?scripts_e02r1pull? combines hyperelastic and viscoelastic material. The run will take only a couple of minutes. Thanks a lot.
/com >> define element type and material
/com >> define keypoint
/com >> generate quater circle
CIRCLE,P51X, , , ,90, ,
/com >> generate straight lines
LSTR, 9, 1
LSTR, 1, 10
LSTR, 8, 4
LSTR, 6, 7
/com >> generate small circle
LARC, 4, 6, 5
/com >> divide lines
LDIV, 1, 0.70
/com >> copy contact zone lines
LGEN,2,7, , ,0,-0.04, , ,0
/com >> shorten copied lines
LDIV, 8, 0.4539919309218844E-01
LSTR, 11, 14
LSTR, 10, 13
/com >> create area
LSBL, 3, 9
/com >> manual meshing the contact zone
/com >> transverse lines
LESIZE,8, , ,5, , , , ,1
LESIZE,7, , ,80, , , , ,1
/com >> mesh the area
/com >> apply surface load
D,all,UY,0, ! APPLY DISPL
/COM, CONTACT PAIR CREATION - START
/COM, CONTACT PAIR PROPERTIES - END
/com >>analysis parameters
/com >>message display control
/com >>nonlinear diagnose
/com >>run the job
/com 2nd step
/com time, 2000
I would be happy to run your simulation is you provide an input file for ABAQUS :wink:
Perhaps someone else can run it through ANSYS...