Télécharger la plaquette Cast3M


fr en

Télécharger frigi2.eso

Retour à la liste

Numérotation des lignes : 

frigi2

  1. C FRIGI2    SOURCE    JK148537  26/08/26    21:15:13     12627          
  2.       SUBROUTINE FRIGI2(MELE,IPMAIL,MINTE,NBPTEL,
  3.      1     IVAMAT,IVACAR,NMATT,NCARR,
  4.      2     CRIGI,CMASS)
  5. ***********************************************************************
  6. *     CALCUL DES COMPOSANTES DE LA RIGIDITE (HOOK) ELASTIQUE
  7. *     CALCUL DES COMPOSANTES DE LA MATRICE (HOOK) DE MASSE
  8. *     .... AU SIGNE PRES
  9. *     BOUCLE SUR LES SS-ZONES DU MODELE DE SECTION
  10. **********************************************************************
  11. * ENTREES :
  12. *
  13. * MELE    = NUMERO  ELEMENT FINI
  14. * IPMAIL = POINTEUR DU MAILLAGE
  15. * NBPTEL =NOMBRE DE POINTS PAR ELEMENT
  16. * IVAMAT =POINTEUR SUR UN SEGMENT MPTVAL DE MATERIAU
  17. * IVACAR =POINTEUR SUR UN SEGMENT MPTVAL DE CARACT. GEOMETRIQUES
  18. * NVARI  =NOMBRE DE COMPOSANTES DE VARIABLES INTERNES
  19. * NMATT  =NOMBRE DE COMPOSNATES DE PROPRIETES DE MATERIAU
  20. * NCARR  =NOMBRE DE COMPOSNATES DE CARACTERISTIQUES GEOMETRIQUES
  21. *
  22. * SORTIES :
  23. * CRIGI(12) RIGIDITE SUR LA FIBRE MOYENNE
  24. * CMASS(12) MASSE SUR LA FIBRE MOYENNE
  25. *
  26. ************************************************************************
  27. *      Pierre Pegon (ISPRA) Juillet/Aout 1993
  28. ***********************************************************************
  29.       IMPLICIT INTEGER(I-N)
  30.       IMPLICIT REAL*8(A-H,O-Z)
  31.  
  32. -INC PPARAM
  33. -INC CCOPTIO
  34. -INC CCHAMP
  35.  
  36. -INC SMCHAML
  37. -INC SMELEME
  38. -INC SMCOORD
  39. -INC SMMODEL
  40. -INC SMINTE
  41.  
  42. -INC TMPTVAL
  43.  
  44.       SEGMENT WRK0
  45.        REAL*8 XMAT(NCXMAT),XCAR(NCXCAR)
  46.       ENDSEGMENT
  47. *
  48.       SEGMENT WRK2
  49.         REAL*8 XE(3,NBBB),SHP(6,NBBB)
  50.       ENDSEGMENT
  51. *
  52.       DIMENSION CRIGI(12),CMASS(12)
  53.  
  54.       MFR =NUMMFR(MELE)
  55.       MELEME=IPMAIL
  56.       NBNN=NUM(/1)
  57.       NBELEM=NUM(/2)
  58. *
  59. *    SEGMENT D'INTEGRATION
  60. *
  61.       SEGACT,MINTE
  62. *
  63. *    INITIALISATION DES SEGMENTS DE TRAVAIL
  64. *
  65.       NCXMAT=NMATT
  66.       NCXCAR=NCARR
  67.       NBBB=NBNN
  68.       SEGINI WRK0,WRK2
  69. *
  70. *       BOUCLE SUR LES ELEMENTS
  71. *
  72.       SEGACT,MCOORD
  73.       DO 1000 IB=1,NBELEM
  74. *
  75. *       ON CHERCHE  LES COORDONNEES DES NOEUDS DE L ELEMENT IB
  76. *
  77.          CALL DOXE(XCOOR,IDIM,NBNN,NUM,IB,XE)
  78. *
  79. *       BOUCLE SUR LES POINTS DE GAUSS
  80. *
  81.          DO 1100 IGAU=1,NBPTEL
  82. *
  83. *       ON CHERCHE LA POSITION DU POINT DE LA SECTION (X->Y) (Y->Z)
  84. *
  85.            YY=0.D0
  86.            ZZ=0.D0
  87.            DO IE1=1,NBNN
  88.              CGAUSS=SHPTOT(1,IE1,IGAU)
  89.              YY=YY+XE(1,IE1)*CGAUSS
  90.              ZZ=ZZ+XE(2,IE1)*CGAUSS
  91.            END DO
  92.            YY2=YY*YY
  93.            ZZ2=ZZ*ZZ
  94. *
  95. *       ON REMPLIT LES SHP ET ON CALCUL LE JACOBIEN
  96. *
  97.            DO IE1=1,6
  98.              DO IE2=1,NBNN
  99.                SHP(IE1,IE2)=SHPTOT(IE1,IE2,IGAU)
  100.              END DO
  101.            END DO
  102. C PPf     CALL JACOBI(XE,SHP,2,NBNN,DJAC)
  103. *
  104. *     ON RECUPERE LES CONSTANTES DU MATERIAU
  105. *
  106.            MPTVAL=IVAMAT
  107.            DO IC=1,NMATT
  108.                  MELVAL=IVAL(IC)
  109. c jk148537 on n'utilise que xmat(1 à 3)
  110. c                 IF(IC.LE.3)THEN
  111. c                   IIC=IC
  112. c                 ELSEIF(IC.LT.(NMATT-4))THEN
  113. c                   IIC=IC+3
  114. c                 ELSEIF(IC.LE.(NMATT-2))THEN
  115. c                   IIC=5+IC-NMATT+2
  116. c                 ELSE
  117. c                   IIC=IC
  118. c                 ENDIF
  119.                   iic =ic
  120.  
  121.                  IF(MELVAL.NE.0)THEN
  122.                     IF(TYVAL(IC)(1:8).NE.'POINTEUR')THEN
  123.                        IBMN=MIN(IB,VELCHE(/2))
  124.                       IGMN=MIN(IGAU,VELCHE(/1))
  125.                        XMAT(IIC)=VELCHE(IGMN,IBMN)
  126.                     ELSE
  127.                       IBMN=MIN(IB,IELCHE(/2))
  128.                      IGMN=MIN(IGAU,IELCHE(/1))
  129.                       XMAT(IIC)=IELCHE(IGMN,IBMN)
  130.                     ENDIF
  131.                  ELSE
  132.                     XMAT(IIC)=0.D0
  133.                     IF(TYVAL(IC)(1:8).EQ.'POINTEUR') THEN
  134.                        XMAT(IIC)=0.D0
  135.                     END IF
  136.                  ENDIF
  137.            END DO
  138. *
  139. *     ON RECUPERE LES CARACTERISTIQUES GEOMETRIQUES
  140. *
  141.            MPTVAL=IVACAR
  142.            DO IC=1,NCARR
  143.               MELVAL=IVAL(IC)
  144.               IF(MELVAL.NE.0)THEN
  145.                     IF(TYVAL(IC)(1:8).NE.'POINTEUR')THEN
  146.                        IBMN=MIN(IB,VELCHE(/2))
  147.                        IGMN=MIN(IGAU,VELCHE(/1))
  148.                         XCAR(IC)=VELCHE(IGMN,IBMN)
  149.                       ELSE
  150.                        IBMN=MIN(IB,IELCHE(/2))
  151.                         IGMN=MIN(IGAU,IELCHE(/1))
  152.                         XCAR(IC)=IELCHE(IGMN,IBMN)
  153.                     ENDIF
  154.                  ELSE
  155.                     XCAR(IC)=0.D0
  156.                     IF(TYVAL(IC)(1:8).EQ.'POINTEUR') THEN
  157.                        XCAR(IC)=0.D0
  158.                     END IF
  159.                  ENDIF
  160. *
  161. *          write(6,*) 'frigi2 caracteristiques ic xcar ',ic,xcar(ic)
  162.            END DO
  163.  
  164. C+PPf
  165. C
  166. C     TRAITEMENT PARTICULIER DES ELEMENTS SEGS(166) ET POJS(167)
  167. C
  168.            IF(MELE.EQ.167)THEN
  169.              DJAC=XCAR(3)
  170.            ELSEIF(MELE.EQ.166)THEN
  171. C+DC on utilise le cas joi3
  172.              CALL JACOBI(XE,SHP,86,NBNN,DJAC)
  173.              DJAC= DJAC*XCAR(3)
  174.            ELSE
  175.              CALL JACOBI(XE,SHP,2,NBNN,DJAC)
  176.            ENDIF
  177. C+PPf
  178. *
  179. *     CONTRIBUTION A CRIGI
  180. *
  181.            PGAUSS=POIGAU(IGAU)*ABS(DJAC)
  182. *
  183.            YOUNG=XMAT(1)
  184.            GAMMA=XMAT(1)/2.D0/(1.D0+XMAT(2))
  185.            ALPH1=XCAR(1)
  186.            ALPH2=XCAR(2)
  187.            CRIGI( 1)=CRIGI( 1)+YOUNG*PGAUSS
  188.            CRIGI( 2)=CRIGI( 2)+YOUNG*YY*PGAUSS
  189.            CRIGI( 3)=CRIGI( 3)+YOUNG*ZZ*PGAUSS
  190.            CRIGI( 4)=CRIGI( 4)+YOUNG*YY2*PGAUSS
  191.            CRIGI( 5)=CRIGI( 5)+YOUNG*YY*ZZ*PGAUSS
  192.            CRIGI( 6)=CRIGI( 6)+YOUNG*ZZ2*PGAUSS
  193. *
  194.            CRIGI( 7)=CRIGI( 7)+ALPH2*GAMMA*PGAUSS
  195.            CRIGI( 8)=CRIGI( 8)+ALPH1*GAMMA*PGAUSS
  196.            CRIGI( 9)=CRIGI( 9)+ALPH2*GAMMA*YY*PGAUSS
  197.            CRIGI(10)=CRIGI(10)+ALPH1*GAMMA*ZZ*PGAUSS
  198.            CRIGI(11)=CRIGI(11)+ALPH2*GAMMA*YY2*PGAUSS
  199.            CRIGI(12)=CRIGI(12)+ALPH1*GAMMA*ZZ2*PGAUSS
  200. *
  201. *     CONTRIBUTION A CMASS
  202. *
  203.          RHO=XMAT(3)         
  204.  
  205. C
  206. C           RHO=XMAT(NMATT)
  207. C
  208.            CMASS( 1)=CMASS( 1)+RHO*PGAUSS
  209.            CMASS( 2)=CMASS( 2)+RHO*YY*PGAUSS
  210.            CMASS( 3)=CMASS( 3)+RHO*ZZ*PGAUSS
  211.            CMASS( 4)=CMASS( 4)+RHO*YY2*PGAUSS
  212.            CMASS( 5)=CMASS( 5)+RHO*YY*ZZ*PGAUSS
  213.            CMASS( 6)=CMASS( 6)+RHO*ZZ2*PGAUSS
  214. *
  215.            CMASS( 7)=CMASS( 7)+RHO*PGAUSS
  216.            CMASS( 8)=CMASS( 8)+RHO*PGAUSS
  217.            CMASS( 9)=CMASS( 9)+RHO*YY*PGAUSS
  218.            CMASS(10)=CMASS(10)+RHO*ZZ*PGAUSS
  219.            CMASS(11)=CMASS(11)+RHO*YY2*PGAUSS
  220.            CMASS(12)=CMASS(12)+RHO*ZZ2*PGAUSS
  221. C
  222. C    FIN DE LA BOUCLE SUR LES POINTS DE GAUSS
  223. C
  224. 1100     CONTINUE
  225. C
  226. C     FIN DE LA BOUCLE SUR LES ELEMENTS
  227. C
  228. 1000  CONTINUE
  229. *
  230.       SEGSUP WRK0,WRK2
  231. *
  232.       RETURN
  233.       END
  234.  
  235.  
  236.  
  237.  

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