Télécharger vonmis0.eso

Retour à la liste

Numérotation des lignes :

  1. C VONMIS0 SOURCE KICH 18/10/04 21:15:42 9946
  2. C VONMIS0
  3. FUNCTION VONMIS0(S,NSTRS,MFR,IFOUR,EP1,ALFAH)
  4.  
  5. IMPLICIT INTEGER(I-N)
  6. IMPLICIT REAL*8(A-H,O-Z)
  7. DIMENSION S(*)
  8.  
  9. C---------ELEMENTS MASSIFS MFR=1 / 31-----------------------------------
  10.  
  11. IF(MFR.EQ.1.or.mfr.eq.31) THEN
  12.  
  13. SOM=0.D0
  14. DO I=4,NSTRS
  15. SOM=SOM+S(I)*S(I)
  16. ENDDO
  17.  
  18. 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)
  19. & +3.D0*SOM
  20. C PPw VONMIS0=SQRT(SEQ)
  21.  
  22. C---------ELEMENTS COQUES AVEC OU SANS CT (MFR=3 ET MFR=9)--------------
  23.  
  24. ELSE IF(MFR.EQ.3.OR.MFR.EQ.9) THEN
  25.  
  26. AUX1=1.D0/(EP1*EP1)
  27. AUX2=ALFAH*36.D0*AUX1*AUX1
  28.  
  29. IF(IFOUR.GE.1) THEN
  30. SEQ=(S(1)*S(1)+S(2)*S(2)+3.D0*S(3)*S(3)-S(1)*S(2))*AUX1
  31. & +(S(4)*S(4)+S(5)*S(5)+3.D0*S(6)*S(6)-S(4)*S(5))*AUX2
  32. C PPw VONMIS0=SQRT(SEQ)
  33. ELSE
  34. C PPw VONMIS0=SQRT((S(1)*S(1)+S(2)*S(2)-S(1)*S(2))*AUX1
  35. C PPw& +(S(3)*S(3)+S(4)*S(4)-S(3)*S(4))*AUX2)
  36. SEQ=(S(1)*S(1)+S(2)*S(2)-S(1)*S(2))*AUX1
  37. & +(S(3)*S(3)+S(4)*S(4)-S(3)*S(4))*AUX2
  38. ENDIF
  39.  
  40. C---------ELEMENTS COQUES EPAISSES MFR=5--------------------------------
  41.  
  42. ELSE IF(MFR.EQ.5) THEN
  43.  
  44. SOM=0.D0
  45. DO I=3,5
  46. SOM=SOM+S(I)*S(I)
  47. ENDDO
  48. C PPw VONMIS0=SQRT(S(1)*S(1)+S(2)*S(2)-S(1)*S(2)+3*SOM)
  49. SEQ=S(1)*S(1)+S(2)*S(2)-S(1)*S(2)+3*SOM
  50.  
  51. C---------ELEMENTS POUTRES MFR=7----------------------------------------
  52.  
  53. ELSE IF(MFR.EQ.7) THEN
  54.  
  55. C PPw VONMIS0=SQRT(S(1)*S(1)+S(4)*S(4)+S(5)*S(5)+S(6)*S(6))
  56. SEQ=S(1)*S(1)+S(4)*S(4)+S(5)*S(5)+S(6)*S(6)
  57.  
  58. ENDIF
  59. C+PPw
  60. VONMIS0=SQRT(ABS(SEQ))
  61. C+PPw
  62. END
  63.  
  64.  
  65.  
  66.  
  67.  

© Cast3M 2003 - Tous droits réservés.
Mentions légales