Télécharger la plaquette Cast3M


fr en

Télécharger kelvin.eso

Retour à la liste

Numérotation des lignes : 

kelvin

  1. C KELVIN    SOURCE    PV        17/12/08    21:17:31     9660           
  2. C KELVIN    SOURCE    MOUSTACAS      21/07/15    01:38:12     5004
  3.  
  4.       SUBROUTINE KELVIN (wrk52,wrk53,wrk54,IB,IGAU,NBPGAU)
  5. C
  6. C==============================================================
  7. C      chaine de KELVIN
  8. C==============================================================
  9. C----------------------------------------------------------------------
  10. C-ENTREES:
  11. C-------
  12. C-6        = NBR. DE COMPOSANTES DES CONTR. OU DES DEFORM.
  13. C-SIG0(6)  = CONTR. AU DEBUT DU PAS D'INTEGRATION
  14. C-DEPST(6) = INCREMENT DES DEFORM. TOTALES
  15. C-NVARI        = NBR. DE VARIABLES INTERNES
  16. C-VAR0(NVARI)  = VARIABLES INTERNES AU DEBUT DU PAS D'INTEGRATION
  17. C
  18. C-NCOMAT       = NBR. DE CARACTERISTIQUES MECANIQUES DU MATERIAU
  19. C-IVAL(NCOMAT) = INDICE DES COMPOSANTES DE MATERIAU
  20. C-XMAT(NCOMAT) = CARACTERISTIQUES MECANIQUES DU MATERIAU
  21. C-MFR          = INDICE DE LA FORMULATION MECANIQUE
  22. C-ICARA        = NBR. DE CARACT. GEOMETRIQUES DES ELEMENTS FINIS
  23. C-XCAR(ICARA)  = CARACT. GEOMETRIQUES DES ELEMENTS FINIS
  24. C-IFOUR        = DIMENSION DU MODÈLE 3D = 2
  25. C
  26.  
  27.  
  28.       IMPLICIT INTEGER(I-N)
  29.       IMPLICIT REAL*8 (A-H,O-Z)
  30.  
  31.  
  32. -INC PPARAM
  33. -INC CCOPTIO
  34. -INC CCREEL
  35. -INC DECHE
  36.  
  37.       PARAMETER (UNDEMI=.5D0,UN=1.D0,DEUX=2.D0,TROIS=3.D0)
  38.       PARAMETER (QUATRE=4.D0)
  39.       PARAMETER (HEURE= 24.D0,SECONDE=3600.D0)
  40.       REAL*8 K,K1,K2,K3,TR1,TR2,TR3,ACN
  41.       REAL*8 RMF1,RMF2,RMF3
  42.       REAL*8 JPF1,JPF2,JPF3
  43.       DIMENSION AN(6,6),BN(6,6)
  44.       DIMENSION SR0(6,6),SR1(6,6),SR2(6,6),SR3(6,6)
  45.       DIMENSION EF0(6,6),EF1(6,6),EF2(6,6),EF3(6,6)
  46.       DIMENSION DEFTOT(6)
  47.       DIMENSION SIGF(6)
  48.       NVARI = VAR0(/1)
  49. c      WRITE(*,*) 'NVARI', NVARI
  50.  
  51. C TESTER IFOUR POUR VÉRIFIER SI LE MODELE EST BIEN 3D
  52.       IF (INPLAS.EQ.174) THEN
  53. c       WRITE(*,*) 'INPLAS', INPLAS
  54. c       WRITE(*,*) 'CALCUL FLUAGE CHAINE DE KELVIN'
  55. c       WRITE(*,*) 'ENTRER TAU EN JOURS'
  56.       ENDIF
  57.  
  58. c       WRITE(*,*) 'MATEPL kelvin.eso ', MATEPL
  59. c       write(*,*) ' inplas jnppla ',inplas,jnpla
  60.       IF (IFOUR.NE.2) THEN
  61. c       WRITE(*,*) 'IFOUR', IFOUR
  62. c       WRITE(*,*) 'L OPTION DE CALCUL N EST PAS IMPLEMENTE'
  63. c       WRITE(*,*) 'SIMULATION MASSIF 3D ISOTROPE SEULEMENT'
  64.        KERRE = 99
  65.        RETURN
  66.       ENDIF
  67.  
  68. C----------------------------------------------------------------------
  69. C       CARACTERISTIQUES
  70. C----------------------------------------------------------------------
  71.  
  72.       KERRE=0
  73. C list xmax pour verifier la positon de rho et alpha
  74.       epf1=epstf(1)
  75.       epf2=epstf(2)
  76.       epf3=epstf(3)
  77.       epf4=epstf(4)
  78.       epf5=epstf(5)
  79.       epf6=epstf(6)
  80.  
  81.       YOUN=XMAT(1)
  82.       XNU=XMAT(2)
  83.       RHO=XMAT(3)
  84.       ALP=XMAT(4)     
  85.       YF1=XMAT(5)
  86.       TF1=XMAT(6)
  87.       YF2=XMAT(7)
  88.       TF2=XMAT(8)
  89.       YF3=XMAT(9)
  90.       TF3=XMAT(10)
  91.  
  92.       SSXX = SIG0(1)
  93.       SSYY = SIG0(2)   
  94.       SSZZ = SIG0(3)
  95.       SSXY = SIG0(4)      
  96.       SSXZ = SIG0(5)      
  97.       SSYZ = SIG0(6) 
  98.  
  99.       TINI = TEMP0
  100.       TIFI = TEMPF
  101.       DELT = (TEMPF-TEMP0)
  102.  
  103.  
  104.       EPSFL01 = VAR0(1)     
  105.       EPSFL02 = VAR0(2) 
  106.       EPSFL03 = VAR0(3) 
  107.       EPSFL04 = VAR0(4)
  108.       EPSFL05 = VAR0(5)
  109.       EPSFL06 = VAR0(6)
  110.  
  111.       SIGRE01 = VAR0(7)     
  112.       SIGRE02 = VAR0(8) 
  113.       SIGRE03 = VAR0(9) 
  114.       SIGRE04 = VAR0(10)
  115.       SIGRE05 = VAR0(11)
  116.       SIGRE06 = VAR0(12) 
  117.  
  118.       ELAS01 = VAR0(13)
  119.       ELAS02 = VAR0(14)
  120.       ELAS03 = VAR0(15)
  121.       ELAS04 = VAR0(16)
  122.       ELAS05 = VAR0(17)
  123.       ELAS06 = VAR0(18)
  124.  
  125.  
  126. C-MODULE D YOUNG A LA FIN DE LA RELAXATION = 
  127. C-SOMME DE TOUS LES RESSORT EN PARALLELE
  128.       K = UN/((UN/YOUN)+(UN/YF1)+(UN/YF2)+(UN/YF3))
  129. C-MODULE D YOUNG RESSORT 1      
  130.       K1 = (YOUN*YF1)/(YOUN+YF1)
  131.       K2 = (YOUN*YF2)/(YOUN+YF2)
  132.       K3 = (YOUN*YF3)/(YOUN+YF3)
  133. C -(TEMPF)AU STEME PARALLELE 1      
  134.       TR1 = (YF1*TF1)/(YOUN+YF1)
  135.       TR2 = (YF2*TF2)/(YOUN+YF2)
  136.       TR3 = (YF3*TF3)/(YOUN+YF3)    
  137.  
  138.       ACN = UN/((UN+XNU)*(UN-(DEUX*XNU)))
  139.  
  140. C-CALCUL TRIDIMENSIONEL
  141. C-CREATION DE LA MATRICE DE RIGIDITE ELASTIQUE
  142.  
  143.        DO I=1,6
  144.          DO J=1,6
  145.           AN(I,J)=XZERO
  146.           BN(I,J)=XZERO
  147.           SR0(I,J)=XZERO
  148.           SR1(I,J)=XZERO
  149.           SR2(I,J)=XZERO
  150.           SR3(I,J)=XZERO
  151.           EF0(I,J)=XZERO
  152.           EF1(I,J)=XZERO
  153.           EF2(I,J)=XZERO
  154.           EF3(I,J)=XZERO
  155.          ENDDO
  156.        ENDDO
  157.  
  158.       IF(IFOUR.EQ.2) THEN
  159.        AN(1,1)=1.D0
  160.        AN(1,2)=-XNU
  161.        AN(1,3)=-XNU
  162.        AN(2,1)=-XNU
  163.        AN(2,2)=1.D0
  164.        AN(2,3)=-XNU
  165.        AN(3,1)=-XNU
  166.        AN(3,2)=-XNU
  167.        AN(3,3)=1.D0
  168.        AN(4,4)=1.D0+XNU
  169.        AN(5,5)=1.D0+XNU
  170.        AN(6,6)=1.D0+XNU
  171.  
  172.        BN(1,1)=1.D0-XNU
  173.        BN(1,2)=XNU
  174.        BN(1,3)=XNU
  175.        BN(2,1)=XNU
  176.        BN(2,2)=1.D0-XNU
  177.        BN(2,3)=XNU
  178.        BN(3,1)=XNU
  179.        BN(3,2)=XNU
  180.        BN(3,3)=1.D0-XNU
  181.        BN(4,4)=(1.D0-(2.D0*XNU))/2.D0
  182.        BN(5,5)=(1.D0-(2.D0*XNU))/2.D0
  183.        BN(6,6)=(1.D0-(2.D0*XNU))/2.D0 
  184. C MATRICE RIGIDITÉ SANS YOUNG     
  185.       ENDIF 
  186.  
  187.        DO I=1,6
  188.         DO  J=1,6
  189.          AN(I,J) = AN(I,J)
  190.          BN(I,J) = ACN*BN(I,J)        
  191.         ENDDO
  192.        ENDDO      
  193.  
  194. C===========RELAXATION : CALCUL DES CONTRAINTES=======================
  195. C=====================================================================
  196.        RMF1 = (YOUN-K1)*EXP(-(TEMPF/TR1))
  197.        RMF2 = (YOUN-K2)*EXP(-(TEMPF/TR2))
  198.        RMF3 = (YOUN-K3)*EXP(-(TEMPF/TR3))
  199.  
  200.        DO I=1,6
  201.         DO  J=1,6
  202.          SR0(I,J) = BN(I,J)*K 
  203.          SR1(I,J) = BN(I,J)*RMF1
  204.          SR2(I,J) = BN(I,J)*RMF2
  205.          SR3(I,J) = BN(I,J)*RMF3       
  206.         ENDDO
  207.        ENDDO
  208.  
  209. C======CAS ISOTROPE 3D       
  210. C======RESSORT ELASTICITE INSTANTANNÉE    
  211.        S0XX = SR0(1,1)+SR0(2,1)+SR0(3,1)
  212. C       WRITE(*,*) 'SOXX', S0XX
  213.        S0YY = SR0(1,2)+SR0(2,2)+SR0(3,2)
  214.        S0ZZ = SR0(1,3)+SR0(2,3)+SR0(3,3)
  215.        S0XY = SR0(4,4)
  216.        S0XZ = SR0(5,5)
  217.        S0YZ = SR0(6,6)
  218. C======SYSTEME KV1
  219.        S1XX = SR1(1,1)+SR1(2,1)+SR1(3,1)
  220.        S1YY = SR1(1,2)+SR1(2,2)+SR1(3,2) 
  221.        S1ZZ = SR1(1,3)+SR1(2,3)+SR1(3,3)
  222.        S1XY = SR1(4,4)
  223.        S1XZ = SR1(5,5)
  224.        S1YZ = SR1(6,6) 
  225.  
  226. C======SYSTEME KV2
  227.        S2XX = SR2(1,1)+SR2(2,1)+SR2(3,1)
  228.        S2YY = SR2(1,2)+SR2(2,2)+SR2(3,2) 
  229.        S2ZZ = SR2(1,3)+SR2(2,3)+SR2(3,3)
  230.        S2XY = SR2(4,4)
  231.        S2XZ = SR2(5,5)
  232.        S2YZ = SR2(6,6)     
  233.  
  234. C======SYSTEME KV3
  235.        S3XX = SR3(1,1)+SR3(2,1)+SR3(3,1)
  236.        S3YY = SR3(1,2)+SR3(2,2)+SR3(3,2) 
  237.        S3ZZ = SR3(1,3)+SR3(2,3)+SR3(3,3)
  238.        S3XY = SR3(4,4)
  239.        S3XZ = SR3(5,5)
  240.        S3YZ = SR3(6,6) 
  241.  
  242. C======TOTAL DES CONTRAINTES    
  243.        SIGREF1 = (S0XX+S1XX+S2XX+S3XX)*epf1 
  244.        SIGREF2 = (S0YY+S1YY+S2YY+S3YY)*epf2 
  245.        SIGREF3 = (S0ZZ+S1ZZ+S2ZZ+S3ZZ)*epf3 
  246.        SIGREF4 = (S0XY+S1XY+S2XY+S3XY)*epf4 
  247.        SIGREF5 = (S0XZ+S1XZ+S2XZ+S3XZ)*epf5 
  248.        SIGREF6 = (S0YZ+S1YZ+S2YZ+S3YZ)*epf6
  249.  
  250.        SIGF(1) = SIGREF1
  251.        SIGF(2) = SIGREF2
  252.        SIGF(3) = SIGREF3
  253.        SIGF(4) = SIGREF4
  254.        SIGF(5) = SIGREF5
  255.        SIGF(6) = SIGREF6
  256.  
  257. c       DO I=1,6
  258. c       RELL=SIGF(I)
  259. c       WRITE(6,*) ' I= ',I, 'SIGF = ',SIGF(I), 'RELL', RELL
  260. c       ENDDO       
  261.  
  262. C=============FLUAGE : CALCUL DES DÉFORMATIONS=========================
  263. C======================================================================
  264.  
  265.        JPF1 = (UN/YF1)*(UN-EXP(-(TEMPF/TF1)))
  266.        JPF2 = (UN/YF2)*(UN-EXP(-(TEMPF/TF2)))
  267.        JPF3 = (UN/YF3)*(UN-EXP(-(TEMPF/TF3)))
  268.  
  269.        DO I=1,6
  270.         DO  J=1,6           
  271.          EF0(I,J)=AN(I,J)/YOUN
  272.          EF1(I,J)=AN(I,J)*JPF1
  273.          EF2(I,J)=AN(I,J)*JPF2
  274.          EF3(I,J)=AN(I,J)*JPF3
  275.         ENDDO
  276.        ENDDO   
  277.  
  278. C======RESSORT ELASTIQUE       
  279.        E0XX = EF0(1,1)+EF0(2,1)+EF0(3,1)
  280.        E0YY = EF0(1,2)+EF0(2,2)+EF0(3,2)
  281.        E0ZZ = EF0(1,3)+EF0(2,3)+EF0(3,3) 
  282.        E0XY = EF0(4,4)
  283.        E0XZ = EF0(5,5)
  284.        E0YZ = EF0(6,6)
  285.  
  286. C======SYSTEME KV1
  287.        E1XX = EF1(1,1)+EF1(2,1)+EF1(3,1)
  288.        E1YY = EF1(1,2)+EF1(2,2)+EF1(3,2)
  289.        E1ZZ = EF1(1,3)+EF1(2,3)+EF1(3,3) 
  290.        E1XY = EF1(4,4)
  291.        E1XZ = EF1(5,5)
  292.        E1YZ = EF1(6,6)
  293.  
  294. C======SYSTEME KV2
  295.        E2XX = EF2(1,1)+EF2(2,1)+EF2(3,1)
  296.        E2YY = EF2(1,2)+EF2(2,2)+EF2(3,2)
  297.        E2ZZ = EF2(1,3)+EF2(2,3)+EF2(3,3) 
  298.        E2XY = EF2(4,4)
  299.        E2XZ = EF2(5,5)
  300.        E2YZ = EF2(6,6)    
  301.  
  302. C======SYSTEME KV3
  303.        E3XX = EF3(1,1)+EF3(2,1)+EF3(3,1)
  304.        E3YY = EF3(1,2)+EF3(2,2)+EF3(3,2)
  305.        E3ZZ = EF3(1,3)+EF3(2,3)+EF3(3,3) 
  306.        E3XY = EF3(4,4)
  307.        E3XZ = EF3(5,5)
  308.        E3YZ = EF3(6,6)    
  309.  
  310. C======DÉFORMATIONS DUES AU FLUAGE 
  311.        EPSFLF1 = (E1XX+E2XX+E3XX)*SIGREF1
  312.        EPSFLF2 = (E1YY+E2YY+E3YY)*SIGREF2
  313.        EPSFLF3 = (E1ZZ+E2ZZ+E3ZZ)*SIGREF3
  314.        EPSFLF4 = (E1XY+E2XY+E3XY)*SIGREF4
  315.        EPSFLF5 = (E1XZ+E2XZ+E3XZ)*SIGREF5
  316.        EPSFLF6 = (E1YZ+E2YZ+E3YZ)*SIGREF6
  317.  
  318. c======DEFORMATIONS ELASTIQUES
  319.        ELAS1 = E0XX*SIGREF1
  320.        ELAS2 = E0YY*SIGREF2
  321.        ELAS3 = E0ZZ*SIGREF3
  322.        ELAS4 = E0XY*SIGREF4
  323.        ELAS5 = E0XZ*SIGREF5
  324.        ELAS6 = E0YZ*SIGREF6
  325.  
  326. C======TOTAL DES DÉFORMATION : DEFORMATIONS DUES AU FLUAGE +DEFORMATIONS
  327. C======DUES A L'ELASTICITÉ INSTANTANNÉE DU RESSORT          
  328.        DEF1 = EPSFLF1+ ELAS1
  329.        DEF2 = EPSFLF2+ ELAS2
  330.        DEF3 = EPSFLF3+ ELAS3
  331.        DEF4 = EPSFLF4+ ELAS4 
  332.        DEF5 = EPSFLF5+ ELAS5 
  333.        DEF6 = EPSFLF6+ ELAS6
  334.  
  335.        DEFTOT(1) = DEF1
  336.        DEFTOT(2) = DEF2
  337.        DEFTOT(3) = DEF3
  338.        DEFTOT(4) = DEF4
  339.        DEFTOT(5) = DEF5
  340.        DEFTOT(6) = DEF6
  341.  
  342. c       DO I=1,6
  343. c       IF(EPSTF(I).NE.0.D0) THEN
  344. c          RELL=(DEFTOT(I)-EPSTF(I))/EPSTF(I)
  345. c       ELSE
  346. c          RELL=(DEFTOT(I)-EPSTF(I))
  347. c       ENDIF
  348.  
  349. c       WRITE(6,*) ' I= ',I, 'EPSTF = ',EPSTF(I),
  350. c     &     'DEFTOT = ', DEFTOT(I), 'RELL', RELL
  351. c      ENDDO         
  352.  
  353. C===== VARIABLE INTERNE AU LA FIN DU PAS                       
  354.        VARF(1)= EPSFLF1
  355.        VARF(2)= EPSFLF2
  356.        VARF(3)= EPSFLF3        
  357.        VARF(4)= EPSFLF4
  358.        VARF(5)= EPSFLF5
  359.        VARF(6)= EPSFLF6
  360.        VARF(7)= SIGREF1
  361.        VARF(8)= SIGREF2
  362.        VARF(9)= SIGREF3      
  363.        VARF(10)= SIGREF4      
  364.        VARF(11)= SIGREF5      
  365.        VARF(12)= SIGREF6
  366.        VARF(13)=ELAS1
  367.        VARF(14)=ELAS2
  368.        VARF(15)=ELAS3
  369.        VARF(16)=ELAS4
  370.        VARF(17)=ELAS5
  371.        VARF(18)=ELAS6 
  372.  
  373.       END      
  374.  
  375.  
  376.  
  377.  
  378.  

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