Télécharger la plaquette Cast3M


fr en

Télécharger ddfds.eso

Retour à la liste

Numérotation des lignes : 

ddfds

  1. C DDFDS     SOURCE    CB215821  26/10/26    21:15:24     12682          
  2.       SUBROUTINE DDFDS(IC,ALFADV,ALFAD1,ALFAD2,YUNG,XNU,S1,S2,DDF,KERRE)
  3. C
  4.       IMPLICIT INTEGER(I-N)
  5.       IMPLICIT REAL*8(A-H,O-Z)
  6.  
  7. -INC PPARAM
  8. -INC CCOPTIO
  9. C
  10.       DIMENSION S1(6),S2(6),DDF(6),WW1(3)
  11. C
  12. C INITIALISATIONS
  13. C
  14.       X=0.D0
  15.       ALF=0.D0
  16.       DO 1 I=1,6
  17.       DDF(I)=0.D0
  18.     1 CONTINUE
  19. C
  20.       IF(XNU.NE.0.5D0.AND.XNU.NE.-1.D0) GO TO 2
  21.       WRITE(IOIMP,10) XNU
  22.       KERRE=640
  23.       RETURN
  24. C
  25.    2  GO TO(100,200,3,400,500),IC
  26.       WRITE(IOIMP,11) IC
  27.       KERRE=640
  28.       RETURN
  29. C
  30. C CRITERE DE LA POROSITE (1)
  31. C
  32.   100 ALF=-0.3333333333333333333333333D0
  33.       GO TO 4
  34. C
  35. C CRITERE DE DRUCKER DUCTILE (2)
  36. C
  37.   200 ALF=ALFADV
  38.       GO TO 4
  39. C
  40. C CRITERE DE DRUCKER FRAGILE (4)
  41. C
  42.   400 ALF=ALFAD1
  43.       GO TO 4
  44. C
  45. C CRITERE DE DRUCKER FRAGILE ECROUI. (5)
  46. C
  47.   500 ALF=ALFAD2
  48. C
  49.     4 X=(ALF*YUNG)/(1.D0-2.D0*XNU)
  50. C
  51.       DDF(1)=X
  52.       DDF(2)=X
  53.       DDF(3)=X
  54. C
  55.       IF(IC.EQ.1) RETURN
  56. C
  57.     3 SEQ2=AVM(S1,S1)
  58.       SEQ=SQRT(SEQ2)
  59.       SMOY=(S1(1)+S1(2)+S1(3))/3.D0
  60.       IF(SEQ.GT.0.D0) GO TO 5
  61.       SEQ2=AVM(S2,S2)
  62.       SEQ=SQRT(SEQ2)
  63.       SMOY=(S2(1)+S2(2)+S2(3))/3.D0
  64.       IF(SEQ.GT.0.D0) GO TO 8
  65.       WW1(1)=S2(1)
  66.       WW1(2)=S2(2)
  67.       WW1(3)=S2(4)
  68.       CALL DIAGOD(WW1)
  69.       ANGLE=WW1(3)*0.01745329252D0
  70.       CO=COS(ANGLE)
  71.       SII=SIN(ANGLE)
  72.       CC=CO*CO
  73.       SS=SII*SII
  74.       CS=CO*SII
  75.       UTIL=YUNG*ALF/((1.D0+XNU)*(1.D0-2.D0*XNU))
  76.       UTI=1.D0-XNU
  77.       WW1(1)=((UTI+2.D0*XNU)*UTIL)
  78.       WW1(2)=WW1(1)
  79.       WW1(3)=0.D0
  80.       CALL RETOUR(WW1,DDF,CC,SS,CS,1)
  81.       DDF(3)=((UTI+2.D0*XNU)*UTIL)
  82.       DDF(5)=0.D0
  83.       DDF(6)=0.D0
  84.       RETURN
  85. C
  86.     5 Y=(YUNG*1.5D0)/(SEQ*(1.D0+XNU))
  87.       DDF(1)=DDF(1)+Y*(S1(1)-SMOY)
  88.       DDF(2)=DDF(2)+Y*(S1(2)-SMOY)
  89.       DDF(3)=DDF(3)+Y*(S1(3)-SMOY)
  90.       DDF(4)=Y*S1(4)
  91.       DDF(5)=Y*S1(5)
  92.       DDF(6)=Y*S1(6)
  93.       GO TO 99
  94.     8 Y=(YUNG*1.5D0)/(SEQ*(1.D0+XNU))
  95.       DDF(1)=DDF(1)+Y*(S2(1)-SMOY)
  96.       DDF(2)=DDF(2)+Y*(S2(2)-SMOY)
  97.       DDF(3)=DDF(3)+Y*(S2(3)-SMOY)
  98.       DDF(4)=Y*S2(4)
  99.       DDF(5)=Y*S2(5)
  100.       DDF(6)=Y*S2(6)
  101. C
  102.    10 FORMAT(1X,'ERREUR DANS DDFDS XNU =',1PD12.5)
  103.    11 FORMAT(1X,'ERREUR DANS DDFDS DANS LA VALEUR DE IC  =',I4)
  104. C
  105.    99 RETURN
  106.       END
  107.  
  108.  
  109.  

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