ABAQUS error while solving constant triaxiality RVE analysis

classic Classic list List threaded Threaded
1 message Options
Reply | Threaded
Open this post in threaded view
|

ABAQUS error while solving constant triaxiality RVE analysis

Akash Bhuwal
Dear All,
I am trying to model RVE with constant triaxiality with/without void using
ABAQUS. I am following the paper
Article Representative volume element calculations under constant st...
<https://www.researchgate.net/publication/272743177_Representative_volume_element_calculations_under_constant_stress_triaxiality_Lode_parameter_and_shear_ratio>
The model is for simplified 2D uniaxial tension and I am using MPC
subroutine and .INP file attached.
I am getting an error saying "ELEMENT PART-1-1.1 IS DISTORTING
EXCESSIVELY." I am not sure how to proceed. Please help.
Regards,
Akash

  ----------

      SUBROUTINE MPC(UE,A,JDOF,MDOF,N,JTYPE,X,U,UINIT,MAXDOF,LMPC,

     * KSTEP,KINC,TIME,NT,NF,TEMP,FIELD,LTRAN,TRAN)

C

      INCLUDE 'ABA_PARAM.INC'

C

      DIMENSION UE(MDOF),A(MDOF,MDOF,N),JDOF(MDOF,N),X(6,N),

     * U(MAXDOF,N),UINIT(MAXDOF,N),TIME(2),TEMP(NT,N),

     * FIELD(NF,NT,N),LTRAN(N),TRAN(3,3,N)

C

      IF (JTYPE .EQ. 1) THEN

C  

        rho=1

C

        u1m=U(1,2)

        u2m=U(2,2)

        u2n2=U(2,3)

C

        UE(1)=(u1m+2*rho*(1+u2m)/(1+u1m)*(u2n2-u2m))

       

        A(1,1,1)=1.0

        A(1,1,2)=-1+2*rho*(1+u2m)*(u2n2-u2m)/(1+u1m)**2

        A(1,2,2)=-2*rho*(u2n2-2*u2m-1)/(1+u1m)

        A(1,1,3)=-2*rho*(1+u2m)/(u1m)

C

        JDOF(1,1) = 1

        JDOF(1,2) = 1

        JDOF(2,2) = 2

        JDOF(1,3) = 2

      END IF

C

      RETURN

      END

[Non-text portions of this message have been removed]