RE: what's wrong with the UMAT subroutine

  1.  Note that ABA_PARAM.INC contains the line ‘IMPLICIT REAL*8 (A-H,O-Z)’.  This means that variables starting with I through N are assumed to be integers.  Thus, your ‘Lambda’ variable would be an integer.
  2.  At least as posted in the e-mail, there seem to be some extraneous characters, for example:

      SUBROUTINE UMAT(STRESS,STATEV,DDSDDE,SSE,SPD,SCD,
    ! ; 1 RPL,DDSDDT,DRPLDE,DRPLDT,

     DDSDDE(3,3) = (1! .-poiss)*Lamda/TWO

Regards,

Dave Lindeman
Staff Scientist
Corporate Research Systems Laboratory
3M Center 235-3G-08
St. Paul, MN 55144
651-733-6383

From: [email protected] [mailto:[email protected]]
Sent: Thursday, December 06, 2018 8:01 PM
To: [email protected]
Subject: [EXTERNAL] [Abaqus] what's wrong with the UMAT subroutine


hi,

    Thanks for your time and patience reading this email.

     Recently I've been trying to define a coordinate-dependent Young's modulus in a 2D plane stress model. I've write a user subroutine as below, but there's a error message ('Problem during compilation') when running the subroutine with ABAQUS/CAE.

    Could anyone tell me what's going wrong with the subroutine?  Thanks again!  Frankey

      SUBROUTINE UMAT(STRESS,STATEV,DDSDDE,SSE,SPD,SCD,
    ! ; 1 RPL,DDSDDT,DRPLDE,DRPLDT,
     2 STRAN,DSTRAN,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)
     4 JSTEP(4)
C     VARIED&! nbsp;YOUNG'S MODULUS WITH COORDINATE
 &! nbsp;    PARAMETER (poiss = .3D0, ONE = 1.D0, TWO= 2.D0, YY=20000.D0)
      X=COORDS(2)
      young = YY*X
      Lamda = young/(ONE - poiss*poiss)
      DO I=1,NTENS
       DO J=1,NTENS
        DDSDDE(I,J) = 0.
       ENDDO
      ENDDO
      DDSDDE(1,1) = Lamda
      DDSDDE(2,2) = Lamda
      DDSDDE(1,2) = Lamda*poiss
      DDSDDE(2,1) = Lamda*poiss
      DDSDDE(3,3) = (1! .-poiss)*Lamda/TWO

      DO K=1,NTENS
       DO L=1,NTENS
          STRESS(K)=STRESS(K)+DDSDDE(K,L)*DSTRAIN(L)
       ENDDO
      ENDDO

      RETURN
      END


________________________________
[email protected]





________________________________



3M Notice: This communication is from an [EXTERNAL] sender.
If this email looks suspicious, do NOT click or open any links or attachments in the email. To report a suspicious email, click on the Report Phishing – PhishMe icon in the Outlook ribbon or forward this email using the report email as spam link in the text below.

Click here<https://spam.mmm.com/pem/pages/digestProcess/digestProcess.jsf?content=aedaaa864ecbae94140305c08d43986d739753923c2cddf483049b16f02976bac6c4ff5fb22097b592aab29dac82584f30f60c161a797f8d8d8a621478a2f25250f75f3a4b570e12125b08ef0dd8e5d1dbc7a045714a146255659761471276b9c37b570f20a98a734732948327931bff93d08ba9b6d561fece99c327acde977017d051a1b719785d> to report this email as spam





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