Télécharger la plaquette Cast3M


fr en

Télécharger epthun.eso

Retour à la liste

Numérotation des lignes : 

epthun

  1. C EPTHUN    SOURCE    CB215821  26/08/24    21:16:30     12622          
  2.       SUBROUTINE EPTHUN(MFR,MELE,VALMAT,NSTRS,TEMP,THIF,THM,
  3.      1            THSU,EPAIST,TXR,XLOC,XGLOB,ROTS,RES,KERRE)
  4. *******************************************************************
  5. *  APPELE PAR EPTHP :DEFORMATIONS THERMIQUES POUR LES MATERIAUX
  6. *  UNIDIRECTIONNELS
  7. *
  8. *   ENTREES :
  9. *    MFR    = NUMERO DE LA FORMULATION
  10. *    MELE   = NUMERO DE L'ELEMENT FINI
  11. *    VALMAT = TABLEAU DE MATERIAU
  12. *    NSTRS  = NOMBRE DE COMPOSANTES DES DEFORMATIONS
  13. *    TEMP,THIF,THM,THSU   = TEMPERATURES
  14. *    EPAIST =  EPAISSEUR
  15. *    TXR,XLOC,XGLOB,ROTS =  TABLEAUX DE TRAVAIL
  16. *   SORTIES :
  17. *    RES =  TABLEAU DE DEFORMATIONS
  18. *    KERRE  = INDICATEUR D'ERREUR
  19. *
  20. *   AM1 JANV 95
  21. ***********************************************************************
  22.       IMPLICIT INTEGER(I-N)
  23.       IMPLICIT REAL*8(A-H,O-Z)
  24. *
  25.  
  26. -INC PPARAM
  27. -INC CCOPTIO
  28. -INC CCHAMP
  29. *
  30.       DIMENSION XLOC(3,3),XGLOB(3,3)
  31.       DIMENSION ROTS(NSTRS,*)
  32.       DIMENSION TXR(IDIM,*),RES(NSTRS)
  33.       DIMENSION VALMAT(*)
  34. *
  35.       DIMENSION EPS1(6),EPS2(6)
  36. *
  37.       KERRE=0
  38. *
  39. *     TEST SUR LA FORMULATION
  40. *
  41. *--------------------------------------------------------------
  42. *                      CAS MASSIF ET POREUX
  43. *--------------------------------------------------------------
  44.       IF(MFR.EQ.1.OR.MFR.EQ.33) THEN
  45. *
  46.            CALL ZERO (XGLOB,IDIM,IDIM)
  47.  
  48. *
  49. *     CONTRAINTES PLANES
  50. *            OU
  51. *     DEFORMATIONS PLANES OU AXI OU FOURIER
  52. *
  53.              IF(IFOUR.EQ.-2.OR.IFOUR.EQ.-1.OR.IFOUR.EQ.0
  54.      1        .OR.IFOUR.EQ.1)THEN
  55. *
  56.               ALPH= VALMAT(1)
  57.               XLOC(1,1)=VALMAT(2)
  58.               XLOC(2,1)=VALMAT(3)
  59.               XLOC(1,2)=-XLOC(2,1)
  60.               XLOC(2,2)=XLOC(1,1)
  61. C    DEFINITION DES AXES ORTHO./AXES GLOBAUX
  62.              DO 4312  K=1,IDIM
  63.                DO 4311 J=1,IDIM
  64.                  DO 4310 I=1,IDIM
  65.                    XGLOB(J,K)=TXR(J,I)*XLOC(I,K)+XGLOB(J,K)
  66.  4310       CONTINUE
  67.  4311          CONTINUE
  68.  4312        CONTINUE
  69.               CC=XGLOB(1,1)*XGLOB(1,1)
  70.               SS=XGLOB(2,1)*XGLOB(2,1)
  71.               CS=XGLOB(1,1)*XGLOB(2,1)
  72. *
  73. *  EN AXES D'ORTHOTROPIE
  74. *
  75.                EPS1(1)=ALPH
  76. *
  77. *  CHGT D'AXES
  78. *
  79.                EPS2(1)=CC*EPS1(1)
  80.                EPS2(2)=SS*EPS1(1)
  81.                EPS2(3)=0.D0
  82.                EPS2(4)=(CS*EPS1(1))*2.D0
  83. *
  84.                RES(1)=EPS2(1)*TEMP
  85.                RES(2)=EPS2(2)*TEMP
  86.                RES(3)=0.D0
  87.                RES(4)=EPS2(4)*TEMP
  88.                IF(IFOUR.EQ.1) THEN
  89.                  RES(5)=0.D0
  90.                  RES(6)=0.D0
  91.                ENDIF
  92. *
  93. *  TRIDIMENSIONNEL
  94. *
  95.         ELSE
  96.               ALPH= VALMAT(1)
  97.               XLOC(1,1)=VALMAT(2)
  98.               XLOC(2,1)=VALMAT(3)
  99.               XLOC(3,1)=VALMAT(4)
  100.               XLOC(1,2)=VALMAT(5)
  101.               XLOC(2,2)=VALMAT(6)
  102.               XLOC(3,2)=VALMAT(7)
  103.               CALL CROSS2(XLOC(1,1),XLOC(1,2),XLOC(1,3),IRR)
  104. C    DEFINITION DES AXES ORTHO./AXES GLOBAUX
  105.              DO 4314  K=1,IDIM
  106.                DO 4313 J=1,IDIM
  107.                  DO 2330 I=1,IDIM
  108.                    XGLOB(J,K)=TXR(J,I)*XLOC(I,K)+XGLOB(J,K)
  109.  2330       CONTINUE
  110.  4313          CONTINUE
  111.  4314        CONTINUE
  112. *
  113. *         MATRICE DE TRANSFORMATION
  114. *
  115.            DO 4315 IC=1,3
  116.            DO 2331 IL=1,3
  117.            ROTS(IL,IC)=XGLOB(IL,IC)*XGLOB(IL,IC)
  118.  2331       CONTINUE
  119.  4315      CONTINUE
  120. C
  121.           DO 2332 IL=1,3
  122.            ROTS(IL,4)=XGLOB(IL,1)*XGLOB(IL,2)
  123.            ROTS(IL,5)=XGLOB(IL,2)*XGLOB(IL,3)
  124.            ROTS(IL,6)=XGLOB(IL,1)*XGLOB(IL,3)
  125.  2332       CONTINUE
  126. C
  127.            DO 2333 IC=1,3
  128.            ROTS(4,IC)=2.D0*XGLOB(1,IC)*XGLOB(2,IC)
  129.            ROTS(5,IC)=2.D0*XGLOB(2,IC)*XGLOB(3,IC)
  130.            ROTS(6,IC)=2.D0*XGLOB(1,IC)*XGLOB(3,IC)
  131.  2333       CONTINUE
  132. C
  133.           DO 4316 IL=4,6
  134.           IL1=IL-3
  135.           IL2=IL1+1
  136.           IF(IL2.GT.3)IL2=IL2-3
  137.            DO 2334 IC=4,6
  138.             IC1=IC-3
  139.             IC2=IC1+1
  140.             IF(IC2.GT.3)IC2=IC2-3
  141.             ROTS(IL,IC)=XGLOB(IL1,IC1)*XGLOB(IL2,IC2)+
  142.      .                       XGLOB(IL1,IC2)*XGLOB(IL2,IC1)
  143.  2334        CONTINUE
  144.  4316     CONTINUE
  145.             DO 2335 IC=1,6
  146.             AA=ROTS(6,IC)
  147.             ROTS(6,IC)=ROTS(5,IC)
  148.             ROTS(5,IC)=AA
  149.  2335        CONTINUE
  150.             DO 2336 IL=1,6
  151.             AA=ROTS(IL,6)
  152.             ROTS(IL,6)=ROTS(IL,5)
  153.             ROTS(IL,5)=AA
  154.  2336        CONTINUE
  155. *
  156. *       EN AXES D'ORTHOTROPIE
  157. *
  158.                EPS1(1)=ALPH
  159.               DO 2340 I=2,6
  160.               EPS1(I)=0.D0
  161.  2340         CONTINUE
  162. *
  163. *  CHGT D'AXES
  164. *
  165.              DO 4317 I=1,6
  166.                  EPS2(I)=0.D0
  167.                 DO 2350 J=1,6
  168.                   EPS2(I)= EPS2(I) +ROTS(I,J)*EPS1(J)
  169.  2350         CONTINUE
  170.  4317        CONTINUE
  171. *
  172.             RES(1)=EPS2(1)*TEMP
  173.             RES(2)=EPS2(2)*TEMP
  174.             RES(3)=EPS2(3)*TEMP
  175.             RES(4)=EPS2(4)*TEMP
  176.             RES(5)=EPS2(5)*TEMP
  177.             RES(6)=EPS2(6)*TEMP
  178.         ENDIF
  179. *____________________________________________________________________
  180. *
  181. *     SECTEUR DE CALCUL POUR LES ELEMENTS COQUES MINCES SAUF COQ2
  182. *____________________________________________________________________
  183. *
  184.        ELSE IF(MELE.EQ.27.OR.MELE.EQ.28.OR.
  185.      1                 MELE.EQ.49)THEN
  186.                ALPH=VALMAT(1)
  187.                V1X= VALMAT(2)
  188.                V1Y= VALMAT(3)
  189.                CC=V1X*V1X
  190.                SS=V1Y*V1Y
  191.                CS=V1X*V1Y
  192.                THMO=(THSU+THIF+4.D0*THM)/6.D0
  193.                THDI=THSU-THIF
  194.                FACM=THMO
  195.                FACF=THDI/EPAIST
  196. *
  197. *  EN AXES D'ORTHOTROPIE
  198. *
  199.                EPS1(1)=ALPH
  200.                EPS1(2)=0.D0
  201.                EPS1(3)=0.D0
  202. *
  203. *  CHGT D'AXES
  204. *
  205.                 EPS2(1)=CC*EPS1(1)
  206.                 EPS2(2)=SS*EPS1(1)
  207.                 EPS2(3)=CS*EPS1(1)*2.D0
  208. *
  209.                IF(IFOUR.GT.0) THEN
  210.                 RES(1)=EPS2(1)*FACM
  211.                 RES(2)=EPS2(2)*FACM
  212.                 RES(3)=EPS2(3)*FACM
  213.                 RES(4)=(EPS2(1)*FACF)
  214. *2.D0
  215.                 RES(5)=(EPS2(2)*FACF)
  216. *2.D0
  217.                 RES(6)=(EPS2(3)*FACF)
  218. *2.D0
  219.                ELSE IF(IFOUR.LE.0) THEN
  220. *
  221. *                 PAS GRAND SENS DANS CE CAS
  222. *
  223.                  KERRE=19
  224.                  GO TO 99
  225.                ENDIF
  226. C____________________________________________________________________
  227. C
  228. C              SECTEUR DE CALCUL POUR LES ELEMENTS COQ2
  229. C____________________________________________________________________
  230. C
  231.       ELSEIF(MELE.EQ.44)THEN
  232.                ALPH=VALMAT(1)
  233.                V1X =VALMAT(2)
  234.                V1Y =VALMAT(3)
  235.                CC=V1X*V1X
  236.                SS=V1Y*V1Y
  237.                CS=V1X*V1Y
  238.                THMO=(THSU+THIF+4.D0*THM)/6.D0
  239.                THDI=THSU-THIF
  240.                FACM=THMO
  241.                FACF=THDI/EPAIST
  242. *
  243. *  EN AXES D'ORTHOTROPIE
  244. *
  245.                EPS1(1)=ALPH
  246.                EPS1(2)=0.D0
  247.                EPS1(3)=0.D0
  248. *
  249. *  CHGT D'AXES
  250. *
  251.                EPS2(1)=CC*EPS1(1)
  252.                EPS2(2)=SS*EPS1(1)
  253.                EPS2(3)=CS*EPS1(1)*2.D0
  254. *
  255.                IF(IFOUR.GT.0) THEN
  256.                 RES(1)=EPS2(1)*FACM
  257.                 RES(2)=EPS2(2)*FACM
  258.                 RES(3)=EPS2(3)*FACM
  259.                 RES(4)=EPS2(1)*FACF
  260. *2.D0
  261.                 RES(5)=EPS2(2)*FACF
  262. *2.D0
  263.                 RES(6)=EPS2(3)*FACF
  264. *2.D0
  265.                ELSE IF(IFOUR.LE.0) THEN
  266.                 RES(1)=EPS2(1)*FACM
  267.                 RES(2)=EPS2(2)*FACM
  268.                 RES(3)=EPS2(1)*FACF
  269.                 RES(4)=EPS2(2)*FACF
  270. *2.D0
  271.                ELSE
  272.                  KERRE=19
  273.                  GO TO 99
  274.                ENDIF
  275.       ELSE
  276.         KERRE=19
  277.       ENDIF
  278. *
  279.  99   CONTINUE
  280.       RETURN
  281.       END
  282.  
  283.  
  284.  
  285.  
  286.  
  287.  
  288.  
  289.  
  290.  

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