Ask Question

Name:
Title:
Your Question:

Answer Question

Name:
Your Answer:
User Submitted Source Code!


Description:
  mahdimahdimahdi
Language: C/C++
Code:
#include <stdio.h>
int main(int argc, char *argv[])
{
     printf("Hello World!");
     getchar();
     return 0;
}
        SUBROUTINE USDFLD(FIELD,STATEV,PNEWDT,DIRECT,T,CELENT,
    1 TIME,DTIME,CMNAME,ORNAME,NFIELD,NSTATV,NOEL,NPT,LAYER,
    2 KSPT,KSTEP,KINC,NDI,NSHR,COORD,JMAC,JMATYP,MATLAYO,LACCFLA)
C

C
      CHARACTER*80 CMNAME,ORNAME
      CHARACTER*3  FLGRAY(15)
      DIMENSION FIELD(NFIELD),STATEV(NSTATV),DIRECT(3,3),
    1 T(3,3),TIME(2)
      DIMENSION ARRAY(15),JARRAY(15),JMAC(*),JMATYP(*),COORD(*)
C     
C       ***********THIS IS FOR SIMPLE TENSION LOADING****************
C          ***********************************************************
C     ENU IS THE POISSON'S RATIO
      PARAMETER (Y0=25.06D0, Yc=50.18D0, dc=0.508D0)
       PARAMETER (E0=3566D0, ENU0=0.37D0)
       
C
       CALL GETVRM(‘S’, ARRAY, JARRAY, FLGRAY, JRCD,
    1 JMAC, JMATYP, MATLAYO, LACCFLA)
       S11=ARRAY(1)
C
C
       SH=(S11/3.0D0)
       SVM=S11
C       
       IF (SH .GE. 0.0D0) THEN
          Yd0=(-0.00042703D0)*(S11**2.0D0/2.0D0)
          d=STATEV(1)
         IF (d .NE. 0.0D0) THEN
             Yd0=(S11**2.0D0/2.0D0)*((1.507D0-1.028D0*d)/
    1        (3566.0D0*(0.514D0*d**2.0D0-1.507D0*d+0.9948D0)**2.0D0))
              STATEV(3)=Yd0
          ENDIF
C          
       ELSE
          Yd0=(-0.00105812D0)*(S11**2.0D0/6.0D0)
          d=STATEV(1)
          IF (d .NE. 0.0D0) THEN
               Yd0=(S11**2.0D0/6.0D0)*((30.993D0*d**2.0D0+4500.777D0*d
    1         -6657.913D0)/(6358178.0D0
    2         *(0.514d**2.0D0-1.507D0*d+0.9948)**2.0D0))
               STATEV(3)=Yd0
          ENDIF     
       ENDIF
C       
C
       STATEV(3)=Yd0
C     
       IF (Yd0 .GE. Yc) THEN
            d=dc
            STATEV(1)=d
            FIELD(1)=STATEV(1)
            STATEV(2)=FIELD(1)
       ELSEIF (Yd0 .LT. Y0) THEN
            d=0.0D0
            STATEV(1)=d
            FIELD(1)=STATEV(1)
            STATEV(2)=FIELD(1)
       ELSE
           d=dc*((Yd0**2.0D0-Y0**2.0D0)/(Yc**2.0D0-Y0**2.0D0))**0.75D0
           STATEV(1)=d
           FIELD(1)=STATEV(1)
           STATEV(2)=FIELD(1)
       ENDIF
C
      RETURN
      END     
Comments: