Télécharger la plaquette Cast3M


fr en

Télécharger cmaxot.eso

Retour à la liste

Numérotation des lignes : 

cmaxot

  1. C CMAXOT    SOURCE    PV        17/12/08    21:15:56     9660           
  2.       SUBROUTINE CMAXOT(wrk52,wrk53,wrk54,WR12,IB,IGAU,NCHAIN,DT,NWA,
  3.      &                  EPSFLU)
  4. *
  5. *
  6. *==============================================================
  7. *      modele maxott
  8. *==============================================================
  9. *
  10. * entrees
  11. *
  12. *   WRK12
  13. *        EMi(2,NWA(i+1))    i=0 a 8
  14. *                        = EVOLUTION DU MODULE D'YOUNG DE LA BRANCHE i
  15. *        SMi(NSTRS)         i=1 a 8
  16. *                         = VARIABLES INTERNES AU DEBUT DU PAS D'INTEGRATION
  17. *
  18. *        IVALMA(NCXMAT)    = INDEX CARACTERISTIQUES MECANIQUES DU MATERIAU
  19. *
  20. *
  21. *   IB = NUMERO DE L'ELEMENT
  22. *   IGAU = NUMERO DU POINT DE GAUSS
  23. *   NCHAIN = NOMBRE DE BRANCHE DU MODELE DE MAXWELL
  24. *   DT = INCREMENT DE TEMPS
  25. *   NWA = DIMENSION DES EVOLUTIONS EMi
  26. *   EPSFLU = DEFORMATION DE FLUAGE
  27. *
  28. * sorties
  29. *
  30. *        SIGF(NSTRS)= CONTR. A LA FIN DU PAS D'INTEGRATION
  31. *        DEFP(NSTRS)= INCREMENT DES DEFORM. PLASTIQUES A LA FIN DU PAS
  32. *               D'INTEGRATION
  33. *        EPINF(NSTRS) = DEFORMATION INELASTIQUE A LA FIN DU PAS D'INTEGRATION
  34. *
  35. *==============================================================
  36. *
  37.       IMPLICIT INTEGER(I-N)
  38.       IMPLICIT REAL*8(A-H,O-Z)
  39.  
  40. -INC PPARAM
  41. -INC CCOPTIO
  42. -INC DECHE
  43. *
  44. *
  45.       SEGMENT WR12
  46.         REAL*8  EM0(2,NWA(1)),EM1(2,NWA(2)),EM2(2,NWA(3))
  47.         REAL*8  EM3(2,NWA(4)),EM4(2,NWA(5)),EM5(2,NWA(6))
  48.         REAL*8  EM6(2,NWA(7)),EM7(2,NWA(8)),EM8(2,NWA(9))
  49.         REAL*8  SM0(NSTRS),SM1(NSTRS),SM2(NSTRS),SM3(NSTRS)
  50.         REAL*8  SM4(NSTRS),SM5(NSTRS),SM6(NSTRS),SM7(NSTRS)
  51.         REAL*8  SM8(NSTRS)
  52.       ENDSEGMENT
  53. *
  54. *
  55.       DIMENSION NWA(9)
  56.       DIMENSION EPSFLU(8)
  57.       DIMENSION EPSOTT(8)
  58.       DIMENSION DEPFIS(8)
  59.       DIMENSION DEPSF(8)
  60.       DIMENSION SFLU(8)
  61.       DIMENSION SIGM(8)
  62.       REAL*8 ETILDE,COE
  63.       INTEGER IND
  64.  
  65.       NCXMAT = XMAT(/1)
  66.       NSTRS = SIG0(/1)
  67.       NVARI = VAR0(/1)
  68.       ICARA=XCARB(/1)
  69.       ANU=XMAT(2)
  70.       UNANU=1.D0-ANU
  71.       TPS1=TEMP0
  72.       TPS2=TEMP0+DT
  73.       SE1=0.D0
  74.       SE2=0.D0
  75.       ETILDE = 0.D0
  76.  
  77.       IF (IFOUR.EQ.2.AND.MFR.EQ.1) THEN
  78.         IND = 21
  79.       ELSE IF (IFOUR.EQ.-1.OR.IFOUR.EQ.-3
  80.      &         .OR.IFOUR.EQ.0.OR.IFOUR.EQ.1) THEN
  81.         IND = 16
  82.       ELSE
  83.         IND = 13
  84.       ENDIF
  85.  
  86. *======================================================================
  87. *        Formulation TRIDIMENSIONNELLE MASSIVE
  88. *        =====================================
  89. *        Formulation DEFORMATIONS PLANES /AXISYMETRIQUE
  90. *        ==============================================
  91. *----------------------------------------------------------------------
  92. *
  93.       IF ((IFOUR.EQ.2.AND.MFR.EQ.1).OR.
  94.      &    (IFOMOD.EQ.-1.AND.IFOUR.NE.-2).OR.
  95.      &       (IFOMOD.EQ.0.OR.IFOMOD.EQ.1))  THEN
  96.         CONS=1.D0/((1.D0+ANU)*(1.D0-2.D0*ANU))
  97.         DO M=0,NCHAIN-1
  98. *
  99. *----------------------------------------------------------------------
  100. * module d'YOUNG apparent
  101. *----------------------------------------------------------------------
  102.  
  103.           CALL CFICMO(TPS1,TPS2,M,WR12,wrk52,wrk53,NWA,EFIC,TR,E1,E2)
  104.  
  105.           ETILDE = ETILDE + EFIC
  106.  
  107.           SE1=SE1+E1
  108.           SE2=SE2+E2
  109.  
  110.           IF (E1.EQ.0.) THEN
  111.             COE = 1.
  112.           ELSE
  113.             COE = E2 / E1
  114.           ENDIF
  115. *
  116. *----------------------------------------------------------------------
  117. * deformations de fluage
  118. *----------------------------------------------------------------------
  119. *
  120.           FLUA=EXP(-TR*DT)
  121.           DO II=1,NSTRS
  122.             IF (M.EQ.0) THEN
  123.               SFLU(II)=SM0(II)*(COE*FLUA-1.0D0)
  124.             ELSE IF (M.EQ.1) THEN
  125.               SFLU(II)=SFLU(II)+SM1(II)*(COE*FLUA-1.0D0)
  126.             ELSE IF (M.EQ.2) THEN
  127.               SFLU(II)=SFLU(II)+SM2(II)*(COE*FLUA-1.0D0)
  128.             ELSE IF (M.EQ.3) THEN
  129.               SFLU(II)=SFLU(II)+SM3(II)*(COE*FLUA-1.0D0)
  130.             ELSE IF (M.EQ.4) THEN
  131.               SFLU(II)=SFLU(II)+SM4(II)*(COE*FLUA-1.0D0)
  132.             ELSE IF (M.EQ.5) THEN
  133.               SFLU(II)=SFLU(II)+SM5(II)*(COE*FLUA-1.0D0)
  134.             ELSE IF (M.EQ.6) THEN
  135.               SFLU(II)=SFLU(II)+SM6(II)*(COE*FLUA-1.0D0)
  136.             ELSE IF (M.EQ.7) THEN
  137.               SFLU(II)=SFLU(II)+SM7(II)*(COE*FLUA-1.0D0)
  138.             ELSE IF (M.EQ.8) THEN
  139.               SFLU(II)=SFLU(II)+SM8(II)*(COE*FLUA-1.0D0)
  140.             ENDIF
  141.           END DO
  142.         END DO
  143. *
  144.         EPSFLU(1)=(SFLU(1) - ANU*SFLU(2) - ANU*SFLU(3)) / ETILDE
  145.         EPSFLU(2)=(SFLU(2) - ANU*SFLU(3) - ANU*SFLU(1)) / ETILDE
  146.         EPSFLU(3)=(SFLU(3) - ANU*SFLU(1) - ANU*SFLU(2)) / ETILDE
  147.         DO II=4,NSTRS
  148.           EPSFLU(II)=(SFLU(II)*(1.0D0+ANU))/ETILDE
  149.         END DO
  150. *
  151. *----------------------------------------------------------------------
  152. * stockage du module d'elasticite
  153. *----------------------------------------------------------------------
  154. *
  155.         VARF(IND) = ETILDE
  156.         XMAT(1) = ETILDE
  157. *
  158. *----------------------------------------------------------------------
  159. * deformation corrigee pour Ottosen
  160. *----------------------------------------------------------------------
  161. *
  162.         DO I=1,NSTRS
  163.           EPSOTT(I) = DEPST(I) + EPSFLU(I)
  164.         END DO
  165. *
  166. *----------------------------------------------------------------------
  167. * appel a Ottosen
  168. *----------------------------------------------------------------------
  169.  
  170. *
  171.         CALL OTTOSE(INPLAS,SIG0,NSTRS,EPSOTT,VAR0,XMAT,IVALMA,NCXMAT,
  172.      &       XCARB,ICARA,NVARI,SIGF,VARF,DEPFIS,MFR,KERRE,IB,IGAU)
  173. *
  174. *----------------------------------------------------------------------
  175. * deformation corrigee pour Maxwell
  176. *----------------------------------------------------------------------
  177. *
  178.         DO II=1,NSTRS
  179.           DEPSF(II) = DEPST(II) - DEPFIS(II)
  180.         END DO
  181. *
  182. *----------------------------------------------------------------------
  183. * modele de Maxwell
  184. *----------------------------------------------------------------------
  185. *
  186.         DO M=0,NCHAIN-1
  187.           CALL CFICMO(TPS1,TPS2,M,WR12,wrk52,wrk53,NWA,EFIC,TR,E1,E2)
  188.           IF (E1.EQ.0.) THEN
  189.             COE = 1.
  190.           ELSE
  191.             COE = E2 / E1
  192.           ENDIF
  193.           FLUA=EXP(-TR*DT)
  194.           IF (M.EQ.0) THEN
  195.             SM0(1)=SM0(1)*COE*FLUA
  196.      1           +CONS*EFIC*((UNANU)*DEPSF(1)+ANU*(DEPSF(2)+DEPSF(3)))
  197.             SM0(2)=SM0(2)*COE*FLUA
  198.      1           +CONS*EFIC*((UNANU)*DEPSF(2)+ANU*(DEPSF(1)+DEPSF(3)))
  199.             SM0(3)=SM0(3)*COE*FLUA
  200.      1           +CONS*EFIC*((UNANU)*DEPSF(3)+ANU*(DEPSF(1)+DEPSF(2)))
  201.             DO II = 4,NSTRS
  202.               SM0(II)=SM0(II)*COE*FLUA
  203.      1               +CONS*EFIC*(1.D0-2.D0*ANU)*DEPSF(II)
  204.             END DO
  205.             DO II=1,NSTRS
  206.               SIGM(II)=SM0(II)
  207.             END DO
  208.           ELSE IF (M.EQ.1) THEN
  209.             SM1(1)=SM1(1)*COE*FLUA
  210.      1           +CONS*EFIC*((UNANU)*DEPSF(1)+ANU*(DEPSF(2)+DEPSF(3)))
  211.             SM1(2)=SM1(2)*COE*FLUA
  212.      1           +CONS*EFIC*((UNANU)*DEPSF(2)+ANU*(DEPSF(1)+DEPSF(3)))
  213.             SM1(3)=SM1(3)*COE*FLUA
  214.      1           +CONS*EFIC*((UNANU)*DEPSF(3)+ANU*(DEPSF(1)+DEPSF(2)))
  215.             DO II = 4,NSTRS
  216.               SM1(II)=SM1(II)*COE*FLUA
  217.      1               +CONS*EFIC*(1.D0-2.D0*ANU)*DEPSF(II)
  218.             END DO
  219.             DO II=1,NSTRS
  220.               SIGM(II)=SIGM(II)+SM1(II)
  221.             END DO
  222.           ELSE IF (M.EQ.2) THEN
  223.             SM2(1)=SM2(1)*COE*FLUA
  224.      1           +CONS*EFIC*((UNANU)*DEPSF(1)+ANU*(DEPSF(2)+DEPSF(3)))
  225.             SM2(2)=SM2(2)*COE*FLUA
  226.      1           +CONS*EFIC*((UNANU)*DEPSF(2)+ANU*(DEPSF(1)+DEPSF(3)))
  227.             SM2(3)=SM2(3)*COE*FLUA
  228.      1           +CONS*EFIC*((UNANU)*DEPSF(3)+ANU*(DEPSF(1)+DEPSF(2)))
  229.             DO II = 4,NSTRS
  230.               SM2(II)=SM2(II)*COE*FLUA
  231.      1               +CONS*EFIC*(1.D0-2.D0*ANU)*DEPSF(II)
  232.             END DO
  233.             DO II=1,NSTRS
  234.               SIGM(II)=SIGM(II)+SM2(II)
  235.             END DO
  236.           ELSE IF (M.EQ.3) THEN
  237.             SM3(1)=SM3(1)*COE*FLUA
  238.      1           +CONS*EFIC*((UNANU)*DEPSF(1)+ANU*(DEPSF(2)+DEPSF(3)))
  239.             SM3(2)=SM3(2)*COE*FLUA
  240.      1           +CONS*EFIC*((UNANU)*DEPSF(2)+ANU*(DEPSF(1)+DEPSF(3)))
  241.             SM3(3)=SM3(3)*COE*FLUA
  242.      1           +CONS*EFIC*((UNANU)*DEPSF(3)+ANU*(DEPSF(1)+DEPSF(2)))
  243.             DO II = 4,NSTRS
  244.               SM3(II)=SM3(II)*COE*FLUA
  245.      1               +CONS*EFIC*(1.D0-2.D0*ANU)*DEPSF(II)
  246.             END DO
  247.             DO II=1,NSTRS
  248.               SIGM(II)=SIGM(II)+SM3(II)
  249.             END DO
  250.           ELSE IF (M.EQ.4) THEN
  251.             SM4(1)=SM4(1)*COE*FLUA
  252.      1           +CONS*EFIC*((UNANU)*DEPSF(1)+ANU*(DEPSF(2)+DEPSF(3)))
  253.             SM4(2)=SM4(2)*COE*FLUA
  254.      1           +CONS*EFIC*((UNANU)*DEPSF(2)+ANU*(DEPSF(1)+DEPSF(3)))
  255.             SM4(3)=SM4(3)*COE*FLUA
  256.      1           +CONS*EFIC*((UNANU)*DEPSF(3)+ANU*(DEPSF(1)+DEPSF(2)))
  257.             DO II = 4,NSTRS
  258.               SM4(II)=SM4(II)*COE*FLUA
  259.      1               +CONS*EFIC*(1.D0-2.D0*ANU)*DEPSF(II)
  260.             END DO
  261.             DO II=1,NSTRS
  262.               SIGM(II)=SIGM(II)+SM4(II)
  263.             END DO
  264.           ELSE IF (M.EQ.5) THEN
  265.             SM5(1)=SM5(1)*COE*FLUA
  266.      1           +CONS*EFIC*((UNANU)*DEPSF(1)+ANU*(DEPSF(2)+DEPSF(3)))
  267.             SM5(2)=SM5(2)*COE*FLUA
  268.      1           +CONS*EFIC*((UNANU)*DEPSF(2)+ANU*(DEPSF(1)+DEPSF(3)))
  269.             SM5(3)=SM5(3)*COE*FLUA
  270.      1           +CONS*EFIC*((UNANU)*DEPSF(3)+ANU*(DEPSF(1)+DEPSF(2)))
  271.             DO II = 4,NSTRS
  272.               SM5(II)=SM5(II)*COE*FLUA
  273.      1               +CONS*EFIC*(1.D0-2.D0*ANU)*DEPSF(II)
  274.             END DO
  275.             DO II=1,NSTRS
  276.               SIGM(II)=SIGM(II)+SM5(II)
  277.             END DO
  278.           ELSE IF (M.EQ.6) THEN
  279.             SM6(1)=SM6(1)*COE*FLUA
  280.      1           +CONS*EFIC*((UNANU)*DEPSF(1)+ANU*(DEPSF(2)+DEPSF(3)))
  281.             SM6(2)=SM6(2)*COE*FLUA
  282.      1           +CONS*EFIC*((UNANU)*DEPSF(2)+ANU*(DEPSF(1)+DEPSF(3)))
  283.             SM6(3)=SM6(3)*COE*FLUA
  284.      1           +CONS*EFIC*((UNANU)*DEPSF(3)+ANU*(DEPSF(1)+DEPSF(2)))
  285.             DO II = 4,NSTRS
  286.               SM6(II)=SM6(II)*COE*FLUA
  287.      1               +CONS*EFIC*(1.D0-2.D0*ANU)*DEPSF(II)
  288.             END DO
  289.             DO II=1,NSTRS
  290.               SIGM(II)=SIGM(II)+SM6(II)
  291.             END DO
  292.           ELSE IF (M.EQ.7) THEN
  293.             SM7(1)=SM7(1)*COE*FLUA
  294.      1           +CONS*EFIC*((UNANU)*DEPSF(1)+ANU*(DEPSF(2)+DEPSF(3)))
  295.             SM7(2)=SM7(2)*COE*FLUA
  296.      1           +CONS*EFIC*((UNANU)*DEPSF(2)+ANU*(DEPSF(1)+DEPSF(3)))
  297.             SM7(3)=SM7(3)*COE*FLUA
  298.      1           +CONS*EFIC*((UNANU)*DEPSF(3)+ANU*(DEPSF(1)+DEPSF(2)))
  299.             DO II = 4,NSTRS
  300.               SM7(II)=SM7(II)*COE*FLUA
  301.      1               +CONS*EFIC*(1.D0-2.D0*ANU)*DEPSF(II)
  302.             END DO
  303.             DO II=1,NSTRS
  304.               SIGM(II)=SIGM(II)+SM7(II)
  305.             END DO
  306.           ELSE IF (M.EQ.8) THEN
  307.             SM8(1)=SM8(1)*COE*FLUA
  308.      1           +CONS*EFIC*((UNANU)*DEPSF(1)+ANU*(DEPSF(2)+DEPSF(3)))
  309.             SM8(2)=SM8(2)*COE*FLUA
  310.      1           +CONS*EFIC*((UNANU)*DEPSF(2)+ANU*(DEPSF(1)+DEPSF(3)))
  311.             SM8(3)=SM8(3)*COE*FLUA
  312.      1           +CONS*EFIC*((UNANU)*DEPSF(3)+ANU*(DEPSF(1)+DEPSF(2)))
  313.             DO II = 4,NSTRS
  314.               SM8(II)=SM8(II)*COE*FLUA
  315.      1               +CONS*EFIC*(1.D0-2.D0*ANU)*DEPSF(II)
  316.             END DO
  317.             DO II=1,NSTRS
  318.               SIGM(II)=SIGM(II)+SM8(II)
  319.             END DO
  320.           ENDIF
  321.         END DO
  322.         DEFP(1)=DEPSF(1)-((SIGM(1)-ANU*(SIGM(2)+SIGM(3)))/SE2
  323.      1                  -(SIG0(1)-ANU*(SIG0(2)+SIG0(3)))/SE1)
  324.         DEFP(2)=DEPSF(2)-((SIGM(2)-ANU*(SIGM(1)+SIGM(3)))/SE2
  325.      1                  -(SIG0(2)-ANU*(SIG0(1)+SIG0(3)))/SE1)
  326.         DEFP(3)=DEPSF(3)-((SIGM(3)-ANU*(SIGM(2)+SIGM(1)))/SE2
  327.      1                  -(SIG0(3)-ANU*(SIG0(2)+SIG0(1)))/SE1)
  328.         DO II = 4,NSTRS
  329.           DEFP(II)=DEPSF(II)-(1+ANU)*(SIGM(II)/SE2-SIG0(II)/SE1)
  330.         END DO
  331.         DO I=1,NSTRS
  332.           EPINF(I)=EPIN0(I)+DEFP(I)+DEPFIS(I)
  333.         END DO
  334. *
  335. *======================================================================
  336. *        Formulation CONTRAINTES PLANES
  337. *        ==============================
  338. *----------------------------------------------------------------------
  339. *
  340.       ELSE IF (IFOMOD.EQ.-1.AND.IFOUR.EQ.-2) THEN
  341.         CONS=1.D0/(1.D0-(ANU*ANU))
  342.         DO M=0,NCHAIN-1
  343. *
  344. *----------------------------------------------------------------------
  345. * module d'YOUNG apparent
  346. *----------------------------------------------------------------------
  347.  
  348.           CALL CFICMO(TPS1,TPS2,M,WR12,wrk52,wrk53,NWA,EFIC,TR,E1,E2)
  349.           ETILDE = ETILDE + EFIC
  350.           SE1=SE1+E1
  351.           SE2=SE2+E2
  352.           IF (E1.EQ.0.) THEN
  353.             COE = 1.
  354.           ELSE
  355.             COE = E2 / E1
  356.           ENDIF
  357. *
  358. *----------------------------------------------------------------------
  359. * deformations de fluage
  360. *----------------------------------------------------------------------
  361. *
  362.           FLUA=EXP(-TR*DT)
  363.           DO II=1,NSTRS
  364.             IF (M.EQ.0) THEN
  365.               SFLU(II)=SM0(II)*(COE*FLUA-1.0D0)
  366.             ELSE IF (M.EQ.1) THEN
  367.               SFLU(II)=SFLU(II)+SM1(II)*(COE*FLUA-1.0D0)
  368.             ELSE IF (M.EQ.2) THEN
  369.               SFLU(II)=SFLU(II)+SM2(II)*(COE*FLUA-1.0D0)
  370.             ELSE IF (M.EQ.3) THEN
  371.               SFLU(II)=SFLU(II)+SM3(II)*(COE*FLUA-1.0D0)
  372.             ELSE IF (M.EQ.4) THEN
  373.               SFLU(II)=SFLU(II)+SM4(II)*(COE*FLUA-1.0D0)
  374.             ELSE IF (M.EQ.5) THEN
  375.               SFLU(II)=SFLU(II)+SM5(II)*(COE*FLUA-1.0D0)
  376.             ELSE IF (M.EQ.6) THEN
  377.               SFLU(II)=SFLU(II)+SM6(II)*(COE*FLUA-1.0D0)
  378.             ELSE IF (M.EQ.7) THEN
  379.               SFLU(II)=SFLU(II)+SM7(II)*(COE*FLUA-1.0D0)
  380.             ELSE IF (M.EQ.8) THEN
  381.               SFLU(II)=SFLU(II)+SM8(II)*(COE*FLUA-1.0D0)
  382.             ENDIF
  383.           END DO
  384.         END DO
  385. *
  386.         EPSFLU(1)=(SFLU(1) - ANU*SFLU(2) - ANU*SFLU(3)) / ETILDE
  387.         EPSFLU(2)=(SFLU(2) - ANU*SFLU(3) - ANU*SFLU(1)) / ETILDE
  388.         EPSFLU(3)=(SFLU(3) - ANU*SFLU(1) - ANU*SFLU(2)) / ETILDE
  389.         EPSFLU(4)=(1.0D0+ANU) * SFLU(4) / ETILDE
  390. *
  391. *----------------------------------------------------------------------
  392. * stockage du module d'elasticite
  393. *----------------------------------------------------------------------
  394. *
  395.         VARF(IND) = ETILDE
  396.         XMAT(1) = ETILDE
  397. *
  398. *----------------------------------------------------------------------
  399. * deformation corrigee pour Ottosen
  400. *----------------------------------------------------------------------
  401. *
  402.         DO I=1,NSTRS
  403.           EPSOTT(I) = DEPST(I) + EPSFLU(I)
  404.         END DO
  405. *
  406. *----------------------------------------------------------------------
  407. * appel a Ottosen
  408. *----------------------------------------------------------------------
  409. *
  410.  
  411.         CALL OTTOSE(INPLAS,SIG0,NSTRS,EPSOTT,VAR0,XMAT,IVALMA,NCXMAT,
  412.      &       XCARB,ICARA,NVARI,SIGF,VARF,DEPFIS,MFR,KERRE,IB,IGAU)
  413. *
  414. *----------------------------------------------------------------------
  415. * deformation corrigee pour Maxwell
  416. *----------------------------------------------------------------------
  417. *
  418.         DO II=1,NSTRS
  419.           DEPSF(II) = DEPST(II) - DEPFIS(II)
  420.         END DO
  421. *----------------------------------------------------------------------
  422. * modele de Maxwell
  423. *----------------------------------------------------------------------
  424.         DO M=0,NCHAIN-1
  425.           CALL CFICMO(TPS1,TPS2,M,WR12,wrk52,wrk53,NWA,EFIC,TR,E1,E2)
  426.           IF (E1.EQ.0.) THEN
  427.             COE = 1.
  428.           ELSE
  429.             COE = E2 / E1
  430.           ENDIF
  431.           FLUA=EXP(-TR*DT)
  432.           IF (M.EQ.0) THEN
  433.             SM0(1)=SM0(1)*COE*FLUA
  434.      1           +CONS*EFIC*(DEPSF(1)+ANU*DEPSF(2))
  435.             SM0(2)=SM0(2)*COE*FLUA
  436.      1           +CONS*EFIC*(DEPSF(2)+ANU*DEPSF(1))
  437.             SM0(3)=0.D0
  438.             SM0(4)=SM0(4)*COE*FLUA+CONS*EFIC*(UNANU)/2.D0*DEPSF(4)
  439.             DO II=1,NSTRS
  440.               SIGM(II)=SM0(II)
  441.             END DO
  442.           ELSE IF (M.EQ.1) THEN
  443.             SM1(1)=SM1(1)*COE*FLUA
  444.      1           +CONS*EFIC*(DEPSF(1)+ANU*DEPSF(2))
  445.             SM1(2)=SM1(2)*COE*FLUA
  446.      1           +CONS*EFIC*(DEPSF(2)+ANU*DEPSF(1))
  447.             SM1(3)=0.D0
  448.             SM1(4)=SM1(4)*COE*FLUA+CONS*EFIC*(UNANU)/2.D0*DEPSF(4)
  449.             DO II=1,NSTRS
  450.               SIGM(II)=SIGM(II)+SM1(II)
  451.             END DO
  452.           ELSE IF (M.EQ.2) THEN
  453.             SM2(1)=SM2(1)*COE*FLUA
  454.      1           +CONS*EFIC*(DEPSF(1)+ANU*DEPSF(2))
  455.             SM2(2)=SM2(2)*COE*FLUA
  456.      1           +CONS*EFIC*(DEPSF(2)+ANU*DEPSF(1))
  457.             SM2(3)=0.D0
  458.             SM2(4)=SM2(4)*COE*FLUA+CONS*EFIC*(UNANU)/2.D0*DEPSF(4)
  459.             DO II=1,NSTRS
  460.               SIGM(II)=SIGM(II)+SM2(II)
  461.             END DO
  462.           ELSE IF (M.EQ.3) THEN
  463.             SM3(1)=SM3(1)*COE*FLUA
  464.      1           +CONS*EFIC*(DEPSF(1)+ANU*DEPSF(2))
  465.             SM3(2)=SM3(2)*COE*FLUA
  466.      1           +CONS*EFIC*(DEPSF(2)+ANU*DEPSF(1))
  467.             SM3(3)=0.D0
  468.             SM3(4)=SM3(4)*COE*FLUA+CONS*EFIC*(UNANU)/2.D0*DEPSF(4)
  469.             DO II=1,NSTRS
  470.               SIGM(II)=SIGM(II)+SM3(II)
  471.             END DO
  472.           ELSE IF (M.EQ.4) THEN
  473.             SM4(1)=SM4(1)*COE*FLUA
  474.      1           +CONS*EFIC*(DEPSF(1)+ANU*DEPSF(2))
  475.             SM4(2)=SM4(2)*COE*FLUA
  476.      1           +CONS*EFIC*(DEPSF(2)+ANU*DEPSF(1))
  477.             SM4(3)=0.D0
  478.             SM4(4)=SM4(4)*COE*FLUA+CONS*EFIC*(UNANU)/2.D0*DEPSF(4)
  479.             DO II=1,NSTRS
  480.               SIGM(II)=SIGM(II)+SM4(II)
  481.             END DO
  482.           ELSE IF (M.EQ.5) THEN
  483.             SM5(1)=SM5(1)*COE*FLUA
  484.      1           +CONS*EFIC*(DEPSF(1)+ANU*DEPSF(2))
  485.             SM5(2)=SM5(2)*COE*FLUA
  486.      1           +CONS*EFIC*(DEPSF(2)+ANU*DEPSF(1))
  487.             SM5(3)=0.D0
  488.             SM5(4)=SM5(4)*COE*FLUA+CONS*EFIC*(UNANU)/2.D0*DEPSF(4)
  489.             DO II=1,NSTRS
  490.               SIGM(II)=SIGM(II)+SM5(II)
  491.             END DO
  492.           ELSE IF (M.EQ.6) THEN
  493.             SM6(1)=SM6(1)*COE*FLUA
  494.      1           +CONS*EFIC*(DEPSF(1)+ANU*DEPSF(2))
  495.             SM6(2)=SM6(2)*COE*FLUA
  496.      1           +CONS*EFIC*(DEPSF(2)+ANU*DEPSF(1))
  497.             SM6(3)=0.D0
  498.             SM6(4)=SM6(4)*COE*FLUA+CONS*EFIC*(UNANU)/2.D0*DEPSF(4)
  499.             DO II=1,NSTRS
  500.               SIGM(II)=SIGM(II)+SM6(II)
  501.             END DO
  502.           ELSE IF (M.EQ.7) THEN
  503.             SM7(1)=SM7(1)*COE*FLUA
  504.      1           +CONS*EFIC*(DEPSF(1)+ANU*DEPSF(2))
  505.             SM7(2)=SM7(2)*COE*FLUA
  506.      1           +CONS*EFIC*(DEPSF(2)+ANU*DEPSF(1))
  507.             SM7(3)=0.D0
  508.             SM7(4)=SM7(4)*COE*FLUA+CONS*EFIC*(UNANU)/2.D0*DEPSF(4)
  509.             DO II=1,NSTRS
  510.               SIGM(II)=SIGM(II)+SM7(II)
  511.             END DO
  512.           ELSE IF (M.EQ.8) THEN
  513.             SM8(1)=SM8(1)*COE*FLUA
  514.      1           +CONS*EFIC*(DEPSF(1)+ANU*DEPSF(2))
  515.             SM8(2)=SM8(2)*COE*FLUA
  516.      1           +CONS*EFIC*(DEPSF(2)+ANU*DEPSF(1))
  517.             SM8(3)=0.D0
  518.             SM8(4)=SM8(4)*COE*FLUA+CONS*EFIC*(UNANU)/2.D0*DEPSF(4)
  519.             DO II=1,NSTRS
  520.               SIGM(II)=SIGM(II)+SM8(II)
  521.             END DO
  522.           END IF
  523.         END DO
  524.         DEFP(1)=DEPSF(1)-((SIGM(1)-ANU*(SIGM(2)+SIGM(3)))/SE2
  525.      1                  -(SIG0(1)-ANU*(SIG0(2)+SIG0(3)))/SE1)
  526.         DEFP(2)=DEPSF(2)-((SIGM(2)-ANU*(SIGM(1)+SIGM(3)))/SE2
  527.      1                  -(SIG0(2)-ANU*(SIG0(1)+SIG0(3)))/SE1)
  528.         DEFP(3)=DEPSF(3)-((SIGM(3)-ANU*(SIGM(2)+SIGM(1)))/SE2
  529.      1                  -(SIG0(3)-ANU*(SIG0(2)+SIG0(1)))/SE1)
  530.         DEFP(4)=DEPSF(4)-(1+ANU)*(SIGM(4)/SE2-SIG0(4)/SE1)
  531.         DO I=1,NSTRS
  532.           EPINF(I)=EPIN0(I)+DEFP(I)+DEPFIS(I)
  533.         END DO
  534.       ENDIF
  535.  
  536. *======================================================================
  537. *          QUELQUES IMPRESSIONS
  538. *======================================================================
  539. *
  540.       IF (IIMPI.EQ.103) THEN
  541.         WRITE(*,66770) IB,IGAU
  542.         WRITE(*,*) dt,tps2
  543.         WRITE(*,66775) XMAT(1)
  544.         WRITE(*,66775) (XMAT(I),I=1,NCXMAT)
  545.         WRITE(*,66776) (XCARB(I),I=1,ICARA)
  546.         WRITE(*,66773) (VAR0(I),I=1,NVARI)
  547.         WRITE(*,66774) (DEPST(I),I=1,NSTRS)
  548.         WRITE(*,66777) (EPSFLU(I),I=1,NSTRS)
  549.         WRITE(*,66778) (EPSOTT(I),I=1,NSTRS)
  550.         WRITE(*,66779) (DEPFIS(I),I=1,NSTRS)
  551.         WRITE(*,66771) (DEPSF(I),I=1,NSTRS)
  552.         WRITE(*,66794) (DEFP(I),I=1,NSTRS)
  553.         WRITE(*,66784) (SIGM(I),I=1,NSTRS)
  554.         WRITE(*,66772) (SIG0(I),I=1,NSTRS)
  555.         WRITE(*,66783) (SIGF(I),I=1,NSTRS)
  556.         WRITE(*,66795) (EPIN0(I),I=1,NSTRS)
  557.         WRITE(*,66782) (EPINF(I),I=1,NSTRS)
  558.       ENDIF
  559. 66770 FORMAT(////////2X,'ELEMENT  ',I6,2X,'POINT  ',I3//)
  560. 66771 FORMAT(2X,'  DEPSF '/(6(1X,1PE12.5)))
  561. 66772 FORMAT(2X,'  SIG0 '/(6(1X,1PE12.5)))
  562. 66773 FORMAT(2X,'  VAR0 '/(6(1X,1PE12.5)))
  563. 66774 FORMAT(2X,'  DEPST '/(6(1X,1PE12.5)))
  564. 66775 FORMAT(2X,'  XMAT  '/(6(1X,1PE12.5)))
  565. 66776 FORMAT(2X,'  XCARB  '/(6(1X,1PE12.5)))
  566. 66777 FORMAT(2X,'  EPSFLU  '/(6(1X,1PE12.5)))
  567. 66778 FORMAT(2X,'  EPSOTT  '/(6(1X,1PE12.5)))
  568. 66779 FORMAT(2X,'  DEPFIS  '/(6(1X,1PE12.5)))
  569. 66782 FORMAT(2X,'  EPINF  '/(6(1X,1PE12.5)))
  570. 66783 FORMAT(2X,'  SIGF  '/(6(1X,1PE12.5)))
  571. 66784 FORMAT(2X,'  SIGM  '/(6(1X,1PE12.5)))
  572. 66794 FORMAT(2X,'  DEFP  '/(6(1X,1PE12.5)))
  573. 66795 FORMAT(2X,'  EPIN0  '/(6(1X,1PE12.5)))
  574.  
  575. *----------------------------------------------------------------------
  576.       RETURN
  577.       END
  578.  
  579.  
  580.  
  581.  
  582.  
  583.  
  584.  
  585.  
  586.  
  587.  
  588.  
  589.  
  590.  
  591.  
  592.  
  593.  
  594.  
  595.  
  596.  
  597.  
  598.  
  599.  
  600.  
  601.  
  602.  
  603.  
  604.  
  605.  
  606.  
  607.  
  608.  
  609.  
  610.  

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