I want to write a subroutine DISP in Abaqus in Fortran for uni-axial extension of a 3D cylinder.
I want to apply displacement loading in one surface of the cylinder by a smooth function.
The function: d(t)=dmax(t^3+2t^4+3t^5)
How to do it? Thanks.