Télécharger la plaquette Cast3M


fr en

Télécharger hookp.eso

Retour à la liste

Numérotation des lignes : 

hookp

  1. C HOOKP     SOURCE    PV        11/02/18    21:15:00     6875
  2. C
  3.       SUBROUTINE HOOKP(VAR,VARD,NVARI,XMAT,NCOMAT,MFR,NSTRS,DINV)
  4. C
  5. C----------------------------------------------------------------------------
  6. C Objet: Calcul de la variation de la matrice de Hooke endommagee inverse
  7. C        dans le cas de la viscoplasticite avec endommagement
  8. C        anisotrope(materiau 142) en regime isotherme
  9. C----------------------------------------------------------------------------
  10. C
  11. C----------------------------------------------------------------------------
  12. C Entree: NSTRS nombre de composantes des contraintes ou
  13. C           des deformations
  14. C         NVARI nombre de variables internes pilotant les
  15. C           equations du modele
  16. C         VAR(NVARI) tableau des variables internes
  17. C         VARD(NVARI) tableau des variations des variables internes
  18. C         XMAT(NCOMAT) tableau des parametres scalaires du materiau
  19. C           a une temperature T donnee
  20. C           dont XMAT(1) module d'Young et XMAT(2) coeff de Poisson
  21. C         MFR indice de la formulation mecanique(seulement massif
  22. C           pour les materiaux endommageables)
  23. C--------------------------------------------------------------------------------
  24. C
  25. C--------------------------------------------------------------------------------
  26. C Sortie: DINV(NSTRS,NSTRS) variation de l'inverse de la matrice de Hooke endommagee
  27. C _______________________________________________________________________________
  28. C
  29.        IMPLICIT INTEGER(I-N)
  30.        IMPLICIT REAL*8(A-H,O-Z)
  31.  
  32. -INC PPARAM
  33. -INC CCOPTIO
  34.        DIMENSION VAR(*),VARD(*),XMAT(*),DINV(NSTRS,NSTRS)
  35. C
  36. C --- Termes d'endommagement ---
  37. C
  38.       D11=VAR(2)
  39.       D22=VAR(3)
  40.       D33=VAR(4)
  41.       D12=VAR(5)
  42.       D13=VAR(6)
  43.       D23=VAR(7)
  44.  
  45.       DE1 = 1.D0-D11
  46.       DE2 = 1.D0-D22
  47.       DE3 = 1.D0-D33
  48.       IF (D11.GE.1.D0) DE1 = 1.D-6
  49.       IF (D22.GE.1.D0) DE2 = 1.D-6
  50.       IF (D33.GE.1.D0) DE3 = 1.D-6
  51. C      write(6,*) DE1,DE2,DE3
  52. C
  53. C --- Coefficients de la matrive de Hooke inversée ---
  54.  
  55.       YUNG=XMAT(1)
  56.       XNU=XMAT(2)
  57.       COEFF1=(1.D0+XNU)/YUNG
  58.       COEFF2=XNU/YUNG
  59.       COEFF3=1.D0/YUNG
  60. C
  61. C --- Calcul du déterminant de I - D ---
  62.  
  63.       det=(DE1)*(DE2)*(DE3)
  64.      &   - 2.D0*D12*D13*D23-D12*D12*(DE3)
  65.      &   - D23*D23*(DE1)-D13*D13*(DE2)
  66.  
  67.       det2=det*det
  68.       det4=det2*det2
  69. C
  70. C --- Calcul des dérivées partielles de det par rapport aux termes du tenseur d'endommagement
  71.  
  72.       detp11=D23*D23-(DE2)*(DE3)
  73.       detp22=D13*D13-(DE1)*(DE3)
  74.       detp33=D12*D12-(DE1)*(DE2)
  75.       detp12=-2.D0*(D23*D13+(DE3)*D12)
  76.       detp13=-2.D0*(D12*D23+(DE2)*D13)
  77.       detp23=-2.D0*(D12*D13+(DE1)*D23)
  78. C
  79. C --- Calcul des termes de la matrice de l'inverse de I - D (multipliés par det)
  80.  
  81.       A11=(DE2)*(DE3)-D23*D23
  82.       A22=(DE1)*(DE3)-D13*D13
  83.       A33=(DE2)*(DE1)-D12*D12
  84.       A12=D12*(DE3)+D23*D13
  85.       A13=D13*(DE2)+D23*D12
  86.       A23=D23*(DE1)+D13*D12
  87. C
  88. C --- Calcul des termes de la variation de la matrice de Hooke inverse
  89.  
  90.       DINV(1,1)=COEFF3*(vard(2)*(-2.D0*detp11*det*A11*A11)
  91.      &      +vard(3)*(2.D0*(D33-1.D0)*A11*det2-2.D0*detp22*det*A11*A11)
  92.      &      +vard(4)*(2.D0*(D22-1.D0)*A11*det2-2.D0*detp33*det*A11*A11)
  93.      &      +vard(5)*(-2.D0*detp12*det*A11*A11)
  94.      &      +vard(6)*(-2.D0*detp13*det*A11*A11)
  95.      &      +vard(7)*(-4.D0*D23*A11*det2-2.D0*detp23*det*A11*A11))/det4
  96.  
  97.       DINV(1,2)=(COEFF1*(vard(2)*(-2.D0*detp11*det*A12*A12)
  98.      &      +vard(3)*(-2.D0*detp22*det*A12*A12)
  99.      &      +vard(4)*(-2.D0*D12*A12*det2-2.D0*detp33*det*A12*A12)
  100.      &      +vard(5)*(2.D0*(DE3)*A12*det2-2.D0*detp12*det*A12*A12)
  101.      &      +vard(6)*(2.D0*D23*A12*det2-2.D0*detp13*det*A12*A12)
  102.      &      +vard(7)*(2.D0*D13*A12*det2-2.D0*detp23*det*A12*A12))
  103.      &      -COEFF2*(vard(2)*((D33-1.D0)*A11*det2
  104.      &                       -2.D0*detp11*det*A11*A22)
  105.      &      +vard(3)*((D33-1.D0)*A22*det2-2.D0*detp22*det*A11*A22)
  106.      &      +vard(4)*(((D22-1.D0)*A22+(D11-1.D0)*A11)*det2
  107.      &                -2.D0*detp33*det*A11*A22)
  108.      &      +vard(5)*(-2.D0*detp12*det*A11*A22)
  109.      &      +vard(6)*(-2.D0*D13*A11*det2-2.D0*detp13*det*A11*A22)
  110.      &      +vard(7)*(-2.D0*D23*A22*det2-2.D0*detp23*det*A11*A22)))/det4
  111.  
  112.       DINV(1,3)=(COEFF1*(vard(2)*(-2.D0*detp11*det*A13*A13)
  113.      &      +vard(3)*(-2.D0*D13*A13*det2-2.D0*detp22*det*A13*A13)
  114.      &      +vard(4)*(-2.D0*detp33*det*A13*A13)
  115.      &      +vard(5)*(2.D0*D23*A13*det2-2.D0*detp12*det*A13*A13)
  116.      &      +vard(6)*(2.D0*(DE2)*A13*det2-2.D0*detp13*det*A13*A13)
  117.      &      +vard(7)*(2.D0*D12*A13*det2-2.D0*detp23*det*A13*A13))
  118.      &      -COEFF2*(vard(2)*((D22-1.D0)*A11*det2
  119.      &                       -2.D0*detp11*det*A11*A33)
  120.      &      +vard(3)*(((D33-1.D0)*A33+(D11-1.D0)*A11)*det2
  121.      &                -2.D0*detp22*det*A11*A33)
  122.      &      +vard(4)*((D22-1.D0)*A33*det2-2.D0*detp33*det*A11*A33)
  123.      &      +vard(5)*(-2.D0*D12*A11*det2-2.D0*detp12*det*A11*A33)
  124.      &      +vard(6)*(-2.D0*detp13*det*A11*A33)
  125.      &      +vard(7)*(-2.D0*D23*A33*det2-2.D0*detp23*det*A11*A33)))/det4
  126.  
  127.       DINV(1,4)=(2.D0*COEFF3/det4)*(vard(2)*(-2.D0*detp11*det*A11*A12)
  128.      &      +vard(3)*((D33-1.D0)*A12*det2-2.D0*detp22*det*A11*A12)
  129.      &      +vard(4)*(((D22-1.D0)*A12-D12*A11)*det2
  130.      &                -2.D0*detp33*det*A11*A12)
  131.      &      +vard(5)*((DE3)*A11*det2-2.D0*detp12*det*A11*A12)
  132.      &      +vard(6)*(D23*A11*det2-2.D0*detp13*det*A11*A12)
  133.      &      +vard(7)*((-2.D0*D23*A12+D13*A11)*det2
  134.      &                -2.D0*detp23*det*A11*A12))
  135.  
  136.       IF (IFOUR.GT.0) THEN
  137.       DINV(1,5)=(2.D0*COEFF3/det4)*(vard(2)*(-2.D0*detp11*det*A11*A13)
  138.      &      +vard(3)*(((D33-1.D0)*A13-D13*A11)*det2
  139.      &                -2.D0*detp22*det*A11*A13)
  140.      &      +vard(4)*((D22-1.D0)*A13*det2-2.D0*detp33*det*A11*A13)
  141.      &      +vard(5)*(D23*A11*det2-2.D0*detp12*det*A11*A13)
  142.      &      +vard(6)*((DE2)*A11*det2-2.D0*detp13*det*A11*A13)
  143.      &      +vard(7)*((-2.D0*D23*A13+D12*A11)*det2
  144.      &                -2.D0*detp23*det*A11*A13))
  145.  
  146.       DINV(1,6)=(COEFF1*(vard(2)*(-2.D0*detp11*det*A12*A13)
  147.      &      +vard(3)*(-D13*A12*det2-2.D0*detp22*det*A12*A13)
  148.      &      +vard(4)*(-D12*A13*det2-2.D0*detp33*det*A12*A13)
  149.      &      +vard(5)*(((DE3)*A13+D23*A12)*det2
  150.      &                -2.D0*detp12*det*A12*A13)
  151.      &      +vard(6)*((D23*A13+(DE2)*A12)*det2
  152.      &                -2.D0*detp13*det*A12*A13)
  153.      &      +vard(7)*((D13*A13+D12*A12)*det2-2.D0*detp23*det*A12*A13))
  154.      &      -COEFF2*(vard(2)*(-D23*A11*det2-2.D0*detp11*det*A11*A23)
  155.      &      +vard(3)*((D33-1.D0)*A23*det2-2.D0*detp22*det*A11*A23)
  156.      &      +vard(4)*((D22-1.D0)*A23*det2-2.D0*detp33*det*A11*A23)
  157.      &      +vard(5)*(D13*A11*det2-2.D0*detp12*det*A11*A23)
  158.      &      +vard(6)*(D12*A11*det2-2.D0*detp13*det*A11*A23)
  159.      &      +vard(7)*(((DE1)*A11-2.D0*D23*A23)*det2
  160.      &                -2.D0*detp23*det*A11*A23)))*2.D0/det4
  161.       ENDIF
  162.  
  163.       DINV(2,2)=COEFF3*(vard(2)*(2.D0*(D33-1.D0)*A22*det2
  164.      &                  -2.D0*detp11*det*A22*A22)
  165.      &      +vard(3)*(-2.D0*detp22*det*A22*A22)
  166.      &      +vard(4)*(2.D0*(D11-1.D0)*A22*det2-2.D0*detp33*det*A22*A22)
  167.      &      +vard(5)*(-2.D0*detp12*det*A22*A22)
  168.      &      +vard(6)*(-4.D0*D13*A22*det2-2.D0*detp13*det*A22*A22))
  169.      &      +vard(7)*(-2.D0*detp23*det*A22*A22)/det4
  170.  
  171.       DINV(2,3)=(COEFF1*(vard(2)*(-2.D0*D23*A23*det2
  172.      &                            -2.D0*detp11*det*A23*A23)
  173.      &      +vard(3)*(-2.D0*detp22*det*A23*A23)
  174.      &      +vard(4)*(-2.D0*detp33*det*A23*A23)
  175.      &      +vard(5)*(2.D0*D13*A23*det2-2.D0*detp12*det*A23*A23)
  176.      &      +vard(6)*(2.D0*D12*A23*det2-2.D0*detp13*det*A23*A23)
  177.      &      +vard(7)*(2.D0*(DE1)*A23*det2-2.D0*detp23*det*A23*A23))
  178.      &      -COEFF2*(vard(2)*(((D33-1.D0)*A33+(D22-1.D0)*A22)*det2
  179.      &                       -2.D0*detp11*det*A22*A33)
  180.      &      +vard(3)*((D11-1.D0)*A22*det2-2.D0*detp22*det*A22*A33)
  181.      &      +vard(4)*((D11-1.D0)*A33*det2-2.D0*detp33*det*A22*A33)
  182.      &      +vard(5)*(-2.D0*D12*A22*det2-2.D0*detp12*det*A22*A33)
  183.      &      +vard(6)*(-2.D0*D13*A33*det2-2.D0*detp13*det*A22*A33)
  184.      &      +vard(7)*(-2.D0*detp23*det*A22*A33)))/det4
  185.  
  186.       DINV(2,4)=(vard(2)*((D33-1.D0)*A12*det2-2.D0*detp11*det*A22*A12)
  187.      &      +vard(3)*(-2.D0*detp22*det*A22*A12)
  188.      &      +vard(4)*(((D11-1.D0)*A12-D12*A22)*det2
  189.      &                -2.D0*detp33*det*A22*A12)
  190.      &      +vard(5)*((DE3)*A22*det2-2.D0*detp12*det*A22*A12)
  191.      &      +vard(6)*((-2.D0*D13*A12+D23*A22)*det2
  192.      &                -2.D0*detp13*det*A22*A12)
  193.      &      +vard(7)*(D13*A22*det2-2.D0*detp23*det*A22*A12))
  194.      &      *(2.D0*COEFF3/det4)
  195.  
  196.       IF (IFOUR.GT.0) THEN
  197.       DINV(2,5)=(COEFF1*(vard(2)*(-D23*A12*det2
  198.      &                            -2.D0*detp11*det*A12*A23)
  199.      &      +vard(3)*(-2.D0*detp22*det*A12*A23)
  200.      &      +vard(4)*(-D12*A23*det2-2.D0*detp33*det*A12*A23)
  201.      &      +vard(5)*(((DE3)*A23+D13*A12)*det2
  202.      &                -2.D0*detp12*det*A12*A23)
  203.      &      +vard(6)*((D23*A23+D12*A12)*det2-2.D0*detp13*det*A12*A23)
  204.      &      +vard(7)*((D13*A23+(DE1)*A12)*det2
  205.      &                -2.D0*detp23*det*A12*A23))
  206.      &      -COEFF2*(vard(2)*((D33-1.D0)*A13*det2
  207.      &                        -2.D0*detp11*det*A22*A13)
  208.      &      +vard(3)*(-D13*A22*det2-2.D0*detp22*det*A22*A13)
  209.      &      +vard(4)*((D11-1.D0)*A13*det2-2.D0*detp33*det*A22*A13)
  210.      &      +vard(5)*(D23*A22*det2-2.D0*detp12*det*A22*A13)
  211.      &      +vard(6)*(((DE2)*A22-2.D0*D13*A13)*det2
  212.      &                -2.D0*detp13*det*A22*A13)
  213.      &      +vard(7)*(D12*A22*det2-2.D0*detp23*det*A22*A13)))*2.D0/det4
  214.  
  215.       DINV(2,6)=2.D0*COEFF3*(vard(2)*(((D33-1.D0)*A23-D23*A22)*det2
  216.      &                    -2.D0*detp11*det*A22*A23)
  217.      &      +vard(3)*(-2.D0*detp22*det*A22*A23)
  218.      &      +vard(4)*((D11-1.D0)*A23*det2-2.D0*detp33*det*A22*A23)
  219.      &      +vard(5)*(D13*A22*det2-2.D0*detp12*det*A22*A23)
  220.      &      +vard(6)*((-2.D0*D13*A23+D12*A22)*det2
  221.      &                -2.D0*detp13*det*A22*A23)
  222.      &      +vard(7)*((DE1)*A22*det2-2.D0*detp23*det*A22*A23))/det4
  223.       ENDIF
  224.  
  225.       DINV(3,3)=COEFF3*(vard(2)*(2.D0*(D22-1.D0)*A33*det2
  226.      &                  -2.D0*detp11*det*A33*A33)
  227.      &      +vard(3)*(2.D0*(D11-1.D0)*A33*det2-2.D0*detp22*det*A33*A33)
  228.      &      +vard(4)*(-2.D0*detp33*det*A33*A33)
  229.      &      +vard(5)*(-4.D0*D12*A33*det2-2.D0*detp12*det*A33*A33)
  230.      &      +vard(6)*(-2.D0*detp13*det*A33*A33))
  231.      &      +vard(7)*(-2.D0*detp23*det*A33*A33)/det4
  232.  
  233.       DINV(3,4)=(COEFF1*(vard(2)*(-D23*A13*det2-2.D0*detp11*det*A13*A23)
  234.      &      +vard(3)*(-D13*A23*det2-2.D0*detp22*det*A13*A23)
  235.      &      +vard(4)*(-2.D0*detp33*det*A13*A23)
  236.      &      +vard(5)*((D23*A23+D13*A13)*det2-2.D0*detp12*det*A13*A23)
  237.      &      +vard(6)*(((DE2)*A23+D12*A13)*det2
  238.      &                -2.D0*detp13*det*A13*A23)
  239.      &      +vard(7)*((D12*A23+(DE1)*A13)*det2
  240.      &                -2.D0*detp23*det*A13*A23))
  241.      &      -COEFF2*(vard(2)*((D22-1.D0)*A12*det2
  242.      &                        -2.D0*detp11*det*A33*A12)
  243.      &      +vard(3)*((D11-1.D0)*A12*det2-2.D0*detp22*det*A33*A12)
  244.      &      +vard(4)*(-D12*A33*det2-2.D0*detp33*det*A33*A12)
  245.      &      +vard(5)*(((DE3)*A33-2.D0*D12*A12)*det2
  246.      &                -2.D0*detp12*det*A33*A12)
  247.      &      +vard(6)*(D23*A33*det2-2.D0*detp13*det*A33*A12)
  248.      &      +vard(7)*(D13*A33*det2-2.D0*detp23*det*A33*A12)))*2.D0/det4
  249.  
  250.       IF (IFOUR.GT.0) THEN
  251.       DINV(3,5)=2.D0*COEFF3*(vard(2)*((D22-1.D0)*A13*det2
  252.      &                      -2.D0*detp11*det*A33*A13)
  253.      &      +vard(3)*(((D11-1.D0)*A13-D13*A33)*det2
  254.      &                -2.D0*detp22*det*A33*A13)
  255.      &      +vard(4)*(-2.D0*detp33*det*A33*A13)
  256.      &      +vard(5)*((-2.D0*D12*A13+D23*A33)*det2
  257.      &                -2.D0*detp12*det*A33*A13)
  258.      &      +vard(6)*((DE2)*A33*det2-2.D0*detp13*det*A33*A13)
  259.      &      +vard(7)*(D12*A33*det2-2.D0*detp23*det*A33*A13))/det4
  260.  
  261.       DINV(3,6)=2.D0*COEFF3*(vard(2)*(((D22-1.D0)*A23-D23*A33)*det2
  262.      &                                -2.D0*detp11*det*A33*A23)
  263.      &      +vard(3)*((D11-1.D0)*A23*det2-2.D0*detp22*det*A33*A23)
  264.      &      +vard(4)*(-2.D0*detp33*det*A33*A23)
  265.      &      +vard(5)*((-2.D0*D12*A23+D13*A33)*det2
  266.      &                -2.D0*detp12*det*A33*A23)
  267.      &      +vard(6)*(D12*A33*det2-2.D0*detp13*det*A33*A23)
  268.      &      +vard(7)*((DE1)*A33*det2
  269.      &                -2.D0*detp23*det*A33*A23))/det4
  270.       ENDIF
  271.  
  272.       DINV(4,4)=(2.D0/det4)*(COEFF1*(vard(2)*((D33-1.D0)*A11*det2
  273.      &                         -2.D0*detp11*det*(A12*A12+A11*A22))
  274.      &      +vard(3)*((D33-1.D0)*A22*det2
  275.      &                -2.D0*detp22*det*(A12*A12+A11*A22))
  276.      &      +vard(4)*((-2.D0*D12*A12+(D22-1.D0)*A22+(D11-1.D0)*A11)
  277.      &                *det2-2.D0*detp33*det*(A12*A12+A11*A22))
  278.      &      +vard(5)*(2.D0*(DE3)*A12*det2
  279.      &                -2.D0*detp12*det*(A12*A12+A11*A22))
  280.      &      +vard(6)*((2.D0*D23*A12-2.D0*D13*A11)*det2
  281.      &                -2.D0*detp13*det*(A12*A12+A11*A22))
  282.      &      +vard(7)*((2.D0*D13*A12-2.D0*D23*A22)*det2
  283.      &                -2.D0*detp23*det*(A12*A12+A11*A22)))
  284.      &      -2.D0*COEFF2*(vard(2)*(-2.D0*detp11*det*A12*A12)
  285.      &      +vard(3)*(-2.D0*detp22*det*A12*A12)
  286.      &      +vard(4)*(-2.D0*D12*A12*det2-2.D0*detp33*det*A12*A12)
  287.      &      +vard(5)*(2.D0*(DE3)*A12*det2-2.D0*detp12*det*A12*A12)
  288.      &      +vard(6)*(2.D0*D23*A12*det2-2.D0*detp13*det*A12*A12)
  289.      &      +vard(7)*(2.D0*D13*A12*det2-2.D0*detp23*det*A12*A12)))
  290.  
  291.       IF (IFOUR.GT.0) THEN
  292.       DINV(4,5)=(2.D0/det4)*(COEFF1*(vard(2)*(-D23*A11*det2
  293.      &                -2.D0*detp11*det*(A13*A12+A11*A23))
  294.      &      +vard(3)*((-D13*A12+(D33-1.D0)*A23)*det2
  295.      &                -2.D0*detp22*det*(A13*A12+A11*A23))
  296.      &      +vard(4)*((-D12*A13+(D22-1.D0)*A23)*det2
  297.      &                -2.D0*detp33*det*(A13*A12+A11*A23))
  298.      &      +vard(5)*((D23*A12+(DE3)*A13+D13*A11)*det2
  299.      &                -2.D0*detp12*det*(A13*A12+A11*A23))
  300.      &      +vard(6)*(((DE2)*A12+D23*A13+D12*A11)*det2
  301.      &                -2.D0*detp13*det*(A13*A12+A11*A23))
  302.      &      +vard(7)*((D12*A12+D13*A13-2.D0*D23*A23+(DE1)*A11)
  303.      &                *det2-2.D0*detp23*det*(A13*A12+A11*A23)))
  304.      &      -2.D0*COEFF2*(vard(2)*(-2.D0*detp11*det*A12*A13)
  305.      &      +vard(3)*(-D13*A12*det2-2.D0*detp22*det*A12*A13)
  306.      &      +vard(4)*(-D12*A13*det2-2.D0*detp33*det*A12*A13)
  307.      &      +vard(5)*(((DE3)*A13+D23*A12)*det2
  308.      &       -2.D0*detp12*det*A12*A13) +vard(6)*((D23*A13+(DE2)
  309.      &          *A12)*det2-2.D0*detp13*det*A12*A13)
  310.      &      +vard(7)*((D13*A13+D12*A12)*det2-2.D0*detp23*det*A12*A13)))
  311.  
  312.       DINV(4,6)=(2.D0/det4)*(COEFF1*(vard(2)*((-D23*A12
  313.      &          +(D33-1.D0)*A13)*det2-2.D0*detp11*det*(A23*A12+A22*A13))
  314.      &      +vard(3)*(-D13*A22*det2-2.D0*detp22*det*(A23*A12+A22*A13))
  315.      &      +vard(4)*((-D12*A23+(D11-1.D0)*A13)*det2
  316.      &                -2.D0*detp33*det*(A23*A12+A22*A13))
  317.      &      +vard(5)*((D13*A12+(DE3)*A23+D23*A22)*det2
  318.      &                -2.D0*detp12*det*(A23*A12+A22*A13))
  319.      &      +vard(6)*((D12*A12+D23*A23-2.D0*D13*A13+(DE2)*A22)
  320.      &                *det2-2.D0*detp13*det*(A23*A12+A22*A13))
  321.      &      +vard(7)*(((DE1)*A12+D13*A23+D12*A22)*det2
  322.      &                -2.D0*detp23*det*(A23*A12+A22*A13)))
  323.      &      -2.D0*COEFF2*(vard(2)*(-D23*A12*det2-2.D0*detp11*det*A12
  324.      &         *A23) + vard(3)*(-2.D0*detp22*det*A12*A23)
  325.      &      +vard(4)*(-D12*A23*det2-2.D0*detp33*det*A12*A23)
  326.      &      +vard(5)*(((DE3)*A23+D13*A12)*det2
  327.      &                -2.D0*detp12*det*A12*A23)
  328.      &      +vard(6)*((D23*A23+D12*A12)*det2-2.D0*detp13*det*A12*A23)
  329.      &      +vard(7)*((D13*A23+(DE1)*A12)*det2
  330.      &                -2.D0*detp23*det*A12*A23)))
  331.  
  332.       DINV(5,5)=(2.D0/det4)*(COEFF1*(vard(2)*((D22-1.D0)*A11*det2
  333.      &                         -2.D0*detp11*det*(A13*A13+A11*A33))
  334.      &      +vard(3)*((-2.D0*D13*A13+(D33-1.D0)*A33+(D11-1.D0)*A11)*det2
  335.      &                -2.D0*detp22*det*(A13*A13+A11*A33))
  336.      &      +vard(4)*((D22-1.D0)*A33*det2
  337.      &                -2.D0*detp33*det*(A13*A13+A11*A33))
  338.      &      +vard(5)*((2.D0*D23*A13-2.D0*D12*A11)*det2
  339.      &                -2.D0*detp12*det*(A13*A13+A11*A33))
  340.      &      +vard(6)*(2.D0*(DE2)*A13*det2
  341.      &                -2.D0*detp13*det*(A13*A13+A11*A33))
  342.      &      +vard(7)*((2.D0*D12*A13-2.D0*D23*A33)*det2
  343.      &                -2.D0*detp23*det*(A13*A13+A11*A33)))
  344.      &      -2.D0*COEFF2*(vard(2)*(-2.D0*detp11*det*A13*A13)
  345.      &      +vard(3)*(-2.D0*D13*A13*det2-2.D0*detp22*det*A13*A13)
  346.      &      +vard(4)*(-2.D0*detp33*det*A13*A13)
  347.      &      +vard(5)*(2.D0*D23*A13*det2-2.D0*detp12*det*A13*A13)
  348.      &      +vard(6)*(2.D0*(DE2)*A13*det2-2.D0*detp13*det*A13*A13)
  349.      &      +vard(7)*(2.D0*D12*A13*det2-2.D0*detp23*det*A13*A13)))
  350.  
  351.       DINV(5,6)=(2.D0/det4)*(COEFF1*(vard(2)*((-D23*A13
  352.      &          +(D22-1.D0)*A12)*det2-2.D0*detp11*det*(A23*A13+A33*A12))
  353.      &      +vard(3)*((-D13*A23+(D11-1.D0)*A12)*det2
  354.      &                -2.D0*detp22*det*(A23*A13+A33*A12))
  355.      &      +vard(4)*(-D12*A33*det2-2.D0*detp33*det*(A23*A13+A33*A12))
  356.      &      +vard(5)*((D13*A13+D23*A23-2.D0*D12*A12+(DE3)*A33)*det2
  357.      &                -2.D0*detp12*det*(A23*A13+A33*A12))
  358.      &      +vard(6)*((D12*A13+(DE2)*A23+D23*A33)*det2
  359.      &                -2.D0*detp13*det*(A23*A13+A33*A12))
  360.      &      +vard(7)*(((DE1)*A13+D12*A23+D13*A33)*det2
  361.      &                -2.D0*detp23*det*(A23*A13+A33*A12)))
  362.      &      -2.D0*COEFF2*(vard(2)*(-D23*A13*det2-2.D0*detp11*det*A13
  363.      &       *A23) +vard(3)*(-D13*A23*det2-2.D0*detp22*det*A12*A23)
  364.      &      +vard(4)*(-2.D0*detp33*det*A13*A23)
  365.      &      +vard(5)*((D23*A23+D13*A13)*det2-2.D0*detp12*det*A13*A23)
  366.      &      +vard(6)*(((DE2)*A23+D12*A13)*det2
  367.      &                -2.D0*detp13*det*A13*A23)
  368.      &      +vard(7)*((D12*A23+(DE1)*A13)*det2
  369.      &                -2.D0*detp23*det*A13*A23)))
  370.  
  371.       DINV(6,6)=(2.D0/det4)*(COEFF1*(vard(2)*((-2.D0*D23*A23
  372.      &           +(D33-1.D0)*A33+(D22-1.D0)*A22)*det2
  373.      &           -2.D0*detp11*det*(A23*A23+A22*A33))
  374.      &      +vard(3)*((D11-1.D0)*A22*det2
  375.      &                -2.D0*detp22*det*(A23*A23+A22*A33))
  376.      &      +vard(4)*((D11-1.D0)*A33*det2
  377.      &                -2.D0*detp33*det*(A23*A23+A22*A33))
  378.      &      +vard(5)*((2.D0*D13*A23-2.D0*D12*A22)*det2
  379.      &                -2.D0*detp12*det*(A23*A23+A22*A33))
  380.      &      +vard(6)*((2.D0*D12*A23-2.D0*D13*A33)*det2
  381.      &                -2.D0*detp13*det*(A23*A23+A22*A33))
  382.      &      +vard(7)*(2.D0*(DE1)*A23*det2
  383.      &                -2.D0*detp23*det*(A23*A23+A22*A33)))
  384.      &      -2.D0*COEFF2*(vard(2)*(-2.D0*D23*A23*det2-2.D0
  385.      &        *detp11*det*A23*A23) +vard(3)*(-2.D0*detp22
  386.      &        *det*A23*A23) +vard(4)*(-2.D0*detp33*det*A23*A23)
  387.      &      +vard(5)*(2.D0*D13*A23*det2-2.D0*detp12*det*A23*A23)
  388.      &      +vard(6)*(2.D0*D12*A23*det2-2.D0*detp13*det*A23*A23)
  389.      &      +vard(7)*(2.D0*(DE1)*A23*det2
  390.      &                -2.D0*detp23*det*A23*A23)))
  391.       ENDIF
  392.  
  393.       DINV(2,1)=DINV(1,2)
  394.  
  395.       DINV(3,1)=DINV(1,3)
  396.       DINV(3,2)=DINV(2,3)
  397.  
  398.       DINV(4,1)=DINV(1,4)
  399.       DINV(4,2)=DINV(2,4)
  400.       DINV(4,3)=DINV(3,4)
  401.  
  402.       IF (IFOUR.GT.0) THEN
  403.       DINV(5,1)=DINV(1,5)
  404.       DINV(5,2)=DINV(2,5)
  405.       DINV(5,3)=DINV(3,5)
  406.       DINV(5,4)=DINV(4,5)
  407.  
  408.       DINV(6,1)=DINV(1,6)
  409.       DINV(6,2)=DINV(2,6)
  410.       DINV(6,3)=DINV(3,6)
  411.       DINV(6,4)=DINV(4,6)
  412.       DINV(6,5)=DINV(5,6)
  413.       ENDIF
  414.  
  415.       RETURN
  416.       END
  417.  
  418.  
  419.  
  420.  
  421.  
  422.  

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