---

Loch_Iso.txt

Ich schreibe im Moment meine Bachelorarbeit zum Thema FEM-Modellierung von Sheet Moulding Compound. Dabei geht es vorrangig um die Darstellung der elastischen Materialeigenschaften anhand einer UMAT.

Nachdem meine anfängliche sehr lange UMAT mir stets Konvergenz Probleme bereitet hat habe ich sie durch eine sehr einfache UMAT ersetzt und das Modell auch stark vereinfacht. Totzdem erscheint weiterhin folgendes:

AVERAGE FORCE                      1.298E-03  TIME AVG. FORCE      1.298E-03
LARGEST RESIDUAL FORCE            -0.600      AT NODE        297  DOF  2
  INSTANCE: LOCHPLATTE2-1                                                                 
LARGEST INCREMENT OF DISP.        -5.692E-05  AT NODE        569  DOF  2
  INSTANCE: LOCHPLATTE2-1                                                                 
LARGEST CORRECTION TO DISP.      -5.692E-05  AT NODE        569  DOF  2
  INSTANCE: LOCHPLATTE2-1                                                                 
          FORCE    EQUILIBRIUM NOT ACHIEVED WITHIN TOLERANCE.
.
.
.
***NOTE: THE SOLUTION APPEARS TO BE DIVERGING. CONVERGENCE IS JUDGED UNLIKELY.
***ERROR: TOO MANY ATTEMPTS MADE FOR THIS INCREMENT

Das Modell mit normaler Definition von E und nu (ohne UMAT) funktioniert einwandfrei.
Unten befindet sich die einfache UMAT und im Anhang die .inp datei.
Nach zwei Wochen Fehlersuche hoffe ich nun dass mir hier jemand behilflich sein kann.
Könnte es am Element Typ liegen (C3D8)? oder an den eingestellten Optionen des STEP?

      SUBROUTINE UMAT(STRESS,STATEV,DDSDDE,SSE,SPD,SCD,
    1    RPL,DDSDDT,DRPLDE,DRPLDT,
    2    STRAN,DSTRANm,TIME,DTIME,TEMP,DTEMP,PREDEF,DPRED,CMNAME,
    3    NDI,NSHR,NTENS,NSTATV,PROPS,NPROPS,COORDS,DROT,PNEWDT,
    4    CELENT,DFGRD0,DFGRD1,NOEL,NPT,LAYER,KSPT,KSTEP,KINC)
C
      INCLUDE 'ABA_PARAM.INC'
C
      CHARACTER*80 CMNAME
      DIMENSION STRESS(NTENS),STATEV(NSTATV),
    1    DDSDDE(NTENS,NTENS),DDSDDT(NTENS),DRPLDE(NTENS),
    2    STRAN(NTENS),DSTRAN(NTENS),TIME(2),PREDEF(1),DPRED(1),
    3    PROPS(NPROPS),COORDS(3),DROT(3,3),DFGRD0(3,3),DFGRD1(3,3)
C
C
      PARAMETER (ZERO = 0.D0,ONE = 1.D0,TWO = 2.D0,THREE = 3.D0)
C
C
C  ELASTIC PROPERTIES
      EMOD=PROPS(1)
      ENU=PROPS(2)
      EBULK3=EMOD/(ONE-TWO*ENU)
      EG2=EMOD/(ONE+ENU)
      EG=EG2/TWO
      EG3=THREE*EG
      ELAM=(EBULK3-EG2)/THREE
C
C ELASTIC STIFFNESS
C
      DO K1=1, NDI
        DO K2=1, NDI
          DDSDDE(K2,K1)=ELAM
        END DO
        DDSDDE(K1,K1)=EG2+ELAM
      END DO
      DO K1=NDI+1, NTENS
        DDSDDE(K1,K1)=EG
      END DO
C
C CALCULATE STRESS
C
      DO K1=1, NTENS
        DO K2=1, NTENS
          STRESS(K2)=STRESS(K2)+DDSDDE(K2,K1)*DSTRAN(K1)
        END DO
      END DO
C
      RETURN
      END

Eine Antwort auf diesen Beitrag verfassen (mit Zitat/Zitat des Beitrags) IP

---

Das Problem ist gelöst. Ich wurde darauf hingewiesen dass im code "DSTRANm" statt "DSTRAN" steht.

So schnell kann einem ein Tippfehler 2 Wochen seiner Zeit rauben. Unglaublich.

Eine Antwort auf diesen Beitrag verfassen (mit Zitat/Zitat des Beitrags) IP