Télécharger la plaquette Cast3M


fr en

Télécharger clendo.eso

Retour à la liste

Numérotation des lignes : 

clendo

  1. C CLENDO    SOURCE    CHAT      05/01/12    22:04:48     5004
  2.       SUBROUTINE CLENDO (DDP,YUNG,XNU,ET,CISA,IENDO,KERRE)
  3.       IMPLICIT INTEGER(I-N)
  4.       IMPLICIT REAL*8(A-H,O-Z)
  5.       DIMENSION ET(*),CISA(*),DDP(6,6)
  6. C
  7. C***********************************************************************
  8. C***************** CALCUL DU MODULE D'ENDOMMAGEMENT ********************
  9. C***********************************************************************
  10. C
  11. C     INITIALISATION
  12. C
  13.       CALL ZDANUL(DDP,36)
  14.       ET1=ET(1)
  15.       ET2=ET(2)
  16.       ET3=ET(3)
  17.       UPUN=1.D0+XNU
  18.       UMUN=1.D0-XNU
  19.       UMDN=1.D0-2.D0*XNU
  20.       EMET1=YUNG-ET1
  21.       EMET2=YUNG-ET2
  22.       EMET3=YUNG-ET3
  23.       COEF=YUNG/(UPUN*UMDN)
  24. C
  25. C***********************************************************************
  26. C************************ PAS DE FISSURATION ***************************
  27. C***********************************************************************
  28. C
  29.       DDP(1,1)=UMUN*COEF
  30.       DDP(1,2)=XNU*COEF
  31.       DDP(1,3)=DDP(1,2)
  32.       DDP(2,1)=DDP(1,2)
  33.       DDP(2,2)=DDP(1,1)
  34.       DDP(2,3)=DDP(1,2)
  35.       DDP(3,1)=DDP(1,2)
  36.       DDP(3,2)=DDP(1,2)
  37.       DDP(3,3)=DDP(1,1)
  38.       DDP(4,4)=YUNG/(2.D0*UPUN)
  39.       DDP(5,5)=DDP(4,4)
  40.       DDP(6,6)=DDP(4,4)
  41.       IF(IENDO.EQ.0) RETURN
  42. C
  43.       GO TO (100,200,300,400,500,600,700),IENDO
  44.       KERRE=468
  45.       RETURN
  46. C
  47. C***********************************************************************
  48. C******************** FISSURATION EN X SEULEMENT ***********************
  49. C***********************************************************************
  50. C
  51.   100 DENO=(EMET1*UMUN)+(ET1*UPUN*UMDN)
  52.       USDENO=1.D0/DENO
  53.       DDP(1,1)=YUNG*ET1*UMUN*USDENO
  54.       DDP(1,2)=XNU*YUNG*ET1*USDENO
  55.       DDP(1,3)=DDP(1,2)
  56.       DDP(2,1)=DDP(1,2)
  57.       DDP(2,2)=YUNG*((EMET1/UPUN)+(ET1*UMUN))*USDENO
  58.       DDP(2,3)=XNU*YUNG*((EMET1/UPUN)+ET1)*USDENO
  59.       DDP(3,1)=DDP(1,3)
  60.       DDP(3,2)=DDP(2,3)
  61.       DDP(3,3)=DDP(2,2)
  62.       DDP(4,4)=CISA(1)
  63.       DDP(6,6)=CISA(3)
  64.       RETURN
  65. C
  66. C***********************************************************************
  67. C******************** FISSURATION EN Y SEULEMENT ***********************
  68. C***********************************************************************
  69. C
  70.   200 DENO=(EMET2*UMUN)+(ET2*UPUN*UMDN)
  71.       USDENO=1.D0/DENO
  72.       DDP(1,1)=YUNG*((EMET2/UPUN)+(ET2*UMUN))*USDENO
  73.       DDP(1,2)=XNU*YUNG*ET2*USDENO
  74.       DDP(1,3)=XNU*YUNG*((EMET2/UPUN)+ET2)*USDENO
  75.       DDP(2,1)=DDP(1,2)
  76.       DDP(2,2)=YUNG*ET2*UMUN*USDENO
  77.       DDP(2,3)=DDP(1,2)
  78.       DDP(3,1)=DDP(1,3)
  79.       DDP(3,2)=DDP(2,3)
  80.       DDP(3,3)=DDP(1,1)
  81.       DDP(4,4)=CISA(1)
  82.       DDP(5,5)=CISA(2)
  83.       RETURN
  84. C
  85. C***********************************************************************
  86. C******************** FISSURATION EN Z SEULEMENT ***********************
  87. C***********************************************************************
  88. C
  89.   400 DENO=(EMET3*UMUN)+(ET3*UPUN*UMDN)
  90.       USDENO=1.D0/DENO
  91.       DDP(1,1)=YUNG*((EMET3/UPUN)+(ET3*UMUN))*USDENO
  92.       DDP(1,2)=XNU*YUNG*((EMET3/UPUN)+ET3)*USDENO
  93.       DDP(1,3)=XNU*YUNG*ET3*USDENO
  94.       DDP(2,1)=DDP(1,2)
  95.       DDP(2,2)=DDP(1,1)
  96.       DDP(2,3)=DDP(1,3)
  97.       DDP(3,1)=DDP(1,3)
  98.       DDP(3,2)=DDP(2,3)
  99.       DDP(3,3)=YUNG*ET3*UMUN*USDENO
  100.       DDP(5,5)=CISA(2)
  101.       DDP(6,6)=CISA(3)
  102.       RETURN
  103. C
  104. C***********************************************************************
  105. C******************** FISSURATION EN XY ENSEMBLE ***********************
  106. C***********************************************************************
  107. C
  108.   300 DENO=(EMET1*EMET2)+(ET1*EMET2*UPUN*UMUN)+(EMET1*ET2*UPUN*UMUN)+
  109.      .     (ET1*ET2*UPUN*UPUN*UMDN)
  110.       USDENO=1.D0/DENO
  111.       DDP(1,1)=YUNG*ET1*(EMET2+(ET2*UPUN*UMUN))*USDENO
  112.       DDP(1,2)=XNU*YUNG*ET1*ET2*UPUN*USDENO
  113.       DDP(1,3)=XNU*YUNG*ET1*(EMET2+(ET2*UPUN))*USDENO
  114.       DDP(2,1)=DDP(1,2)
  115.       DDP(2,2)=YUNG*ET2*(EMET1+(ET1*UPUN*UMUN))*USDENO
  116.       DDP(2,3)=XNU*YUNG*ET2*(EMET1+(ET1*UPUN))*USDENO
  117.       DDP(3,1)=DDP(1,3)
  118.       DDP(3,2)=DDP(2,3)
  119.       DDP(3,3)=YUNG*((EMET1*EMET2)+(ET1*EMET2)+(EMET1*ET2)+
  120.      .       (ET1*ET2*UPUN*UMUN))*USDENO
  121.       DDP(4,4)=CISA(1)
  122.       DDP(5,5)=CISA(2)
  123.       DDP(6,6)=CISA(3)
  124.       RETURN
  125. C
  126. C***********************************************************************
  127. C******************** FISSURATION EN YZ ENSEMBLE ***********************
  128. C***********************************************************************
  129. C
  130.   600 DENO=(EMET2*EMET3)+(ET2*EMET3*UPUN*UMUN)+(EMET2*ET3*UPUN*UMUN)+
  131.      .     (ET2*ET3*UPUN*UPUN*UMDN)
  132.       USDENO=1.D0/DENO
  133.       DDP(1,1)=YUNG*((EMET2*EMET3)+(ET2*EMET3)+(EMET2*ET3)+
  134.      .       (ET2*ET3*UPUN*UMUN))*USDENO
  135.       DDP(1,2)=XNU*YUNG*ET2*(EMET3+(ET3*UPUN))*USDENO
  136.       DDP(1,3)=XNU*YUNG*ET3*(EMET2+(ET2*UPUN))*USDENO
  137.       DDP(2,1)=DDP(1,2)
  138.       DDP(2,2)=YUNG*ET2*(EMET3+(ET3*UPUN*UMUN))*USDENO
  139.       DDP(2,3)=XNU*YUNG*ET2*ET3*UPUN*USDENO
  140.       DDP(3,1)=DDP(1,3)
  141.       DDP(3,2)=DDP(2,3)
  142.       DDP(3,3)=YUNG*ET3*(EMET2+(ET2*UPUN*UMUN))*USDENO
  143.       DDP(4,4)=CISA(1)
  144.       DDP(5,5)=CISA(2)
  145.       DDP(6,6)=CISA(3)
  146.       RETURN
  147. C
  148. C***********************************************************************
  149. C******************** FISSURATION EN ZX ENSEMBLE ***********************
  150. C***********************************************************************
  151. C
  152.   500 DENO=(EMET3*EMET1)+(ET3*EMET1*UPUN*UMUN)+(EMET3*ET1*UPUN*UMUN)+
  153.      .     (ET3*ET1*UPUN*UPUN*UMDN)
  154.       USDENO=1.D0/DENO
  155.       DDP(1,1)=YUNG*ET1*(EMET3+(ET3*UPUN*UMUN))*USDENO
  156.       DDP(1,2)=XNU*YUNG*ET1*(EMET3+(ET3*UPUN))*USDENO
  157.       DDP(1,3)=XNU*YUNG*ET3*ET1*UPUN*USDENO
  158.       DDP(2,1)=DDP(1,2)
  159.       DDP(2,2)=YUNG*((EMET3*EMET1)+(ET3*EMET1)+(EMET3*ET1)+
  160.      .       (ET3*ET1*UPUN*UMUN))*USDENO
  161.       DDP(2,3)=XNU*YUNG*ET3*(EMET1+(ET1*UPUN))*USDENO
  162.       DDP(3,1)=DDP(1,3)
  163.       DDP(3,2)=DDP(2,3)
  164.       DDP(3,3)=YUNG*ET3*(EMET1+(ET1*UPUN*UMUN))*USDENO
  165.       DDP(4,4)=CISA(1)
  166.       DDP(5,5)=CISA(2)
  167.       DDP(6,6)=CISA(3)
  168.       RETURN
  169. C
  170. C***********************************************************************
  171. C******************** FISSURATION EN XYZ ENSEMBLE **********************
  172. C***********************************************************************
  173. C
  174.   700 DENO=(EMET1*EMET2*EMET3)+
  175.      .   (ET1*EMET2*EMET3)+(EMET1*ET2*EMET3)+(EMET1*EMET2*ET3)+
  176.      .   (UPUN*UMUN*((EMET1*ET2*ET3)+(ET1*EMET2*ET3)+(ET1*ET2*EMET3)))+
  177.      .   (UPUN*UPUN*UMDN*(ET1*ET2*ET3))
  178.       USDENO=1.D0/DENO
  179.       DDP(1,1)=((EMET2*EMET3)+(ET2*EMET3)+(EMET2*ET3)+
  180.      .      (ET2*ET3*UPUN*UMUN))*YUNG*ET1*USDENO
  181.       DDP(1,2)=XNU*YUNG*ET1*ET2*(EMET3+(ET3*UPUN))*USDENO
  182.       DDP(1,3)=XNU*YUNG*ET1*ET3*(EMET2+(ET2*UPUN))*USDENO
  183.       DDP(2,1)=DDP(1,2)
  184.       DDP(2,2)=((EMET1*EMET3)+(ET1*EMET3)+(EMET1*ET3)+
  185.      .      (ET1*ET3*UPUN*UMUN))*YUNG*ET2*USDENO
  186.       DDP(2,3)=XNU*YUNG*ET2*ET3*(EMET1+(ET1*UPUN))*USDENO
  187.       DDP(3,1)=DDP(1,3)
  188.       DDP(3,2)=DDP(2,3)
  189.       DDP(3,3)=((EMET1*EMET2)+(ET1*EMET2)+(EMET1*ET2)+
  190.      .      (ET1*ET2*UPUN*UMUN))*YUNG*ET3*USDENO
  191.       DDP(4,4)=CISA(1)
  192.       DDP(5,5)=CISA(2)
  193.       DDP(6,6)=CISA(3)
  194.       RETURN
  195. C
  196.       END
  197.  
  198.  

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