Sources Esope | Cast3M

vonmis0
C VONMIS0   SOURCE    KICH      18/10/04    21:15:42     9946           
C  VONMIS0
      FUNCTION VONMIS0(S,NSTRS,MFR,IFOUR,EP1,ALFAH)
 
      IMPLICIT INTEGER(I-N)
      IMPLICIT REAL*8(A-H,O-Z)
      DIMENSION S(*)
 
C---------ELEMENTS MASSIFS MFR=1 /  31-----------------------------------
 
      IF(MFR.EQ.1.or.mfr.eq.31) THEN
 
        SOM=0.D0
        DO I=4,NSTRS
          SOM=SOM+S(I)*S(I)
        ENDDO
 
        SEQ=S(1)*S(1)+S(2)*S(2)+S(3)*S(3)-S(1)*S(2)-S(1)*S(3)-S(2)*S(3)
     &  +3.D0*SOM
C PPw   VONMIS0=SQRT(SEQ)
 
C---------ELEMENTS COQUES AVEC OU SANS CT (MFR=3 ET MFR=9)--------------
 
      ELSE IF(MFR.EQ.3.OR.MFR.EQ.9) THEN
 
        AUX1=1.D0/(EP1*EP1)
        AUX2=ALFAH*36.D0*AUX1*AUX1
 
        IF(IFOUR.GE.1) THEN
          SEQ=(S(1)*S(1)+S(2)*S(2)+3.D0*S(3)*S(3)-S(1)*S(2))*AUX1
     &    +(S(4)*S(4)+S(5)*S(5)+3.D0*S(6)*S(6)-S(4)*S(5))*AUX2
C PPw     VONMIS0=SQRT(SEQ)
        ELSE
C PPw     VONMIS0=SQRT((S(1)*S(1)+S(2)*S(2)-S(1)*S(2))*AUX1
C PPw&    +(S(3)*S(3)+S(4)*S(4)-S(3)*S(4))*AUX2)
          SEQ=(S(1)*S(1)+S(2)*S(2)-S(1)*S(2))*AUX1
     &    +(S(3)*S(3)+S(4)*S(4)-S(3)*S(4))*AUX2
        ENDIF
 
C---------ELEMENTS COQUES EPAISSES MFR=5--------------------------------
 
      ELSE IF(MFR.EQ.5) THEN
 
        SOM=0.D0
        DO I=3,5
          SOM=SOM+S(I)*S(I)
        ENDDO
C PPw   VONMIS0=SQRT(S(1)*S(1)+S(2)*S(2)-S(1)*S(2)+3*SOM)
        SEQ=S(1)*S(1)+S(2)*S(2)-S(1)*S(2)+3*SOM
 
C---------ELEMENTS POUTRES MFR=7----------------------------------------
 
      ELSE IF(MFR.EQ.7) THEN
 
C PPw   VONMIS0=SQRT(S(1)*S(1)+S(4)*S(4)+S(5)*S(5)+S(6)*S(6))
        SEQ=S(1)*S(1)+S(4)*S(4)+S(5)*S(5)+S(6)*S(6)
 
      ENDIF
C+PPw
      VONMIS0=SQRT(ABS(SEQ))
C+PPw
      END