Télécharger la plaquette Cast3M


fr en

Télécharger incoms.eso

Retour à la liste

Numérotation des lignes : 

incoms

  1. C INCOMS    SOURCE    GOUNAND   21/06/02    21:16:24     11022          
  2.       SUBROUTINE INCOMS(MYCOMS,CGEOM2,LERF,IMPR,IRET)
  3.       IMPLICIT REAL*8 (A-H,O-Z)
  4.       IMPLICIT INTEGER (I-N)
  5. C***********************************************************************
  6. C NOM         : INCOMS
  7. C PROJET      : Noyau linéaire NLIN
  8. C DESCRIPTION : Initialise le segment contenant les informations sur
  9. C               les lois de comportement
  10. C
  11. C LANGAGE     : ESOPE
  12. C AUTEUR      : Stéphane GOUNAND (CEA/DRN/DMT/SEMT/LTMF)
  13. C               mél : gounand@semt2.smts.cea.fr
  14. C***********************************************************************
  15. C APPELES          :
  16. C APPELES (E/S)    :
  17. C APPELE PAR       :
  18. C***********************************************************************
  19. C ENTREES            : -
  20. C SORTIES            : MYCOMS
  21. C CODE RETOUR (IRET) : = 0 si tout s'est bien passé
  22. C***********************************************************************
  23. C VERSION    : v1, 10/05/04, version initiale
  24. C HISTORIQUE : v1, 10/05/04, création
  25. C HISTORIQUE :
  26. C HISTORIQUE :
  27. C***********************************************************************
  28. C Prière de PRENDRE LE TEMPS de compléter les commentaires
  29. C en cas de modification de ce sous-programme afin de faciliter
  30. C la maintenance !
  31. C***********************************************************************
  32.  
  33. -INC PPARAM
  34. -INC CCOPTIO
  35. -INC SMELEME
  36.       POINTEUR CGEOM2.MELEME
  37. -INC SMLMOTS
  38.       POINTEUR LNCOM.MLMOTS
  39. -INC TNLIN      
  40. *-INC SLCOMP
  41.       POINTEUR MYCOMS.COMPS
  42.       POINTEUR COCOUR.COMP
  43. *
  44.       INTEGER IMPR,IRET
  45.       CHARACTER*8 BLANC
  46.       CHARACTER*8 NCOM
  47.       INTEGER IDIM0,IDIM1,IDIM2,IDIM3,IDIM4
  48.       CHARACTER*1 CDIM0,CDIM1,CDIM2,CDIM3,CDIM4,CMETH,CIH
  49.       CHARACTER*1 CID
  50. *
  51. *
  52. * Executable statements
  53. *
  54.       IF (IMPR.GT.1) WRITE(IOIMP,*) 'Entrée dans incoms'
  55.       BLANC='        '
  57.       SEGINI MYCOMS
  58. *
  59. * Détermination de la dimension de l'espace de dérivation
  60. *
  61.       IF (LERF.NE.0) THEN
  62.          CALL DIMESH(CGEOM2,IDMESH,IMPR,IRET)
  63.          IF (IRET.NE.0) GOTO 9999
  64.          IESDER=IDMESH
  65.       ELSE
  66.          IESDER=IDIM
  67.       ENDIF
  68. *
  69. * Loi de comportement RIEN
  70. *
  71.       NCOM=BLANC
  72.       NCOM='RIEN'
  73.       NCOCOF=0
  74.       SEGINI COCOUR
  75.       COCOUR.NOMCOM=NCOM
  76.       COCOUR.LTREF=.FALSE.
  77.       SEGDES COCOUR
  78.       MYCOMS.LISCOM(**)=COCOUR
  79. *
  80. * Loi de comportement IDEN
  81. *
  82.       NCOM=BLANC
  83.       NCOM='IDEN'
  84.       NCOCOF=1
  85.       SEGINI COCOUR
  86.       COCOUR.NOMCOM=NCOM
  87.       COCOUR.DERCOF(1)=0
  88.       COCOUR.LTREF=.FALSE.
  89.       SEGDES COCOUR
  90.       MYCOMS.LISCOM(**)=COCOUR
  91. *
  92. * Loi de comportement rays = epsi * sigma * T^3
  93. *
  94.       NCOM=BLANC
  95.       NCOM='RAYS'
  96.       NCOCOF=3
  97.       SEGINI COCOUR
  98.       COCOUR.NOMCOM=NCOM
  99.       COCOUR.DERCOF(1)=0
  100.       COCOUR.DERCOF(2)=0
  101.       COCOUR.DERCOF(3)=0
  102.       COCOUR.LTREF=.FALSE.
  103.       SEGDES COCOUR
  104.       MYCOMS.LISCOM(**)=COCOUR
  105. *
  106. * Loi de comportement MUR
  107. *
  108.       NCOM=BLANC
  109.       NCOM='MUR'
  110.       NCOCOF=4
  111.       SEGINI COCOUR
  112.       COCOUR.NOMCOM=NCOM
  113.       COCOUR.DERCOF(1)=0
  114.       COCOUR.DERCOF(2)=0
  115.       COCOUR.DERCOF(3)=0
  116.       COCOUR.DERCOF(4)=0
  117.       COCOUR.LTREF=.FALSE.
  118.       SEGDES COCOUR
  119.       MYCOMS.LISCOM(**)=COCOUR
  120. *
  121. * Loi de comportement SUTH
  122. *
  123.       NCOM=BLANC
  124.       NCOM='SUTH'
  125.       NCOCOF=3
  126.       SEGINI COCOUR
  127.       COCOUR.NOMCOM=NCOM
  128.       COCOUR.DERCOF(1)=0
  129.       COCOUR.DERCOF(2)=0
  130.       COCOUR.DERCOF(3)=0
  131.       COCOUR.LTREF=.FALSE.
  132.       SEGDES COCOUR
  133.       MYCOMS.LISCOM(**)=COCOUR
  134. *
  135. * Loi de comportement DDX{1..IESDER}
  136. *
  137.       DO IID=1,IESDER
  138.          CALL INT2CH(IID,CID,IMPR,IRET)
  139.          IF (IRET.NE.0) GOTO 9999
  140.          NCOM=BLANC
  141.          NCOM='D/DX'
  142.          NCOM(5:5)=CID
  143.          NCOCOF=1
  144.          SEGINI COCOUR
  145.          COCOUR.NOMCOM=NCOM
  146.          COCOUR.DERCOF(1)=1
  147.          COCOUR.LTREF=.FALSE.
  148.          SEGDES COCOUR
  149.          MYCOMS.LISCOM(**)=COCOUR
  150.       ENDDO
  151. *
  152. * Loi de comportement DIV
  153. *
  154.       NCOM=BLANC
  155.       NCOM='DIV'
  156.       NCOCOF=IESDER
  157.       SEGINI COCOUR
  158.       COCOUR.NOMCOM=NCOM
  159.       DO IIDIM=1,IESDER
  160.          COCOUR.DERCOF(IIDIM)=1
  161.       ENDDO
  162.       COCOUR.LTREF=.FALSE.
  163.       SEGDES COCOUR
  164.       MYCOMS.LISCOM(**)=COCOUR
  165. *
  166. * Loi de comportement MAXI
  167. *
  168.       NCOM=BLANC
  169.       NCOM='MAXI'
  170.       NCOCOF=1
  171.       SEGINI COCOUR
  172.       COCOUR.NOMCOM=NCOM
  173.       COCOUR.DERCOF(1)=0
  174.       COCOUR.LTREF=.FALSE.
  175.       SEGDES COCOUR
  176.       MYCOMS.LISCOM(**)=COCOUR
  177. *
  178. * Loi de comportement TAU
  179. *
  180.       NCOM=BLANC
  181.       NCOM='TAU'
  182.       NCOCOF=5
  183.       SEGINI COCOUR
  184.       COCOUR.NOMCOM=NCOM
  185.       COCOUR.DERCOF(1)=0
  186.       COCOUR.DERCOF(2)=0
  187.       COCOUR.DERCOF(3)=0
  188.       COCOUR.DERCOF(4)=0
  189.       COCOUR.DERCOF(5)=0
  190.       COCOUR.LTREF=.FALSE.
  191.       SEGDES COCOUR
  192.       MYCOMS.LISCOM(**)=COCOUR
  193. *
  194. * Loi de comportement SIGMA
  195. *
  196.       NCOM=BLANC
  197.       NCOM='SIGMA'
  198.       NCOCOF=5
  199.       SEGINI COCOUR
  200.       COCOUR.NOMCOM=NCOM
  201.       COCOUR.DERCOF(1)=0
  202.       COCOUR.DERCOF(2)=0
  203.       COCOUR.DERCOF(3)=0
  204.       COCOUR.DERCOF(4)=0
  205.       COCOUR.DERCOF(5)=0
  206.       COCOUR.LTREF=.FALSE.
  207.       SEGDES COCOUR
  208.       MYCOMS.LISCOM(**)=COCOUR
  209. *
  210. * Loi de comportement TAILDIRE
  211. *
  212.       NCOM=BLANC
  213.       NCOM='TAILDIRE'
  214.       NCOCOF=IDIM
  215.       SEGINI COCOUR
  216.       COCOUR.NOMCOM=NCOM
  217.       DO IIDIM=1,IDIM
  218.          COCOUR.DERCOF(IIDIM)=0
  219.       ENDDO
  220.       COCOUR.LTREF=.TRUE.
  221.       SEGDES COCOUR
  222.       MYCOMS.LISCOM(**)=COCOUR
  223. *
  224. * Loi de comportement VOLORI
  225. *
  226.       NCOM=BLANC
  227.       NCOM='VOLORI'
  228.       NCOCOF=0
  229.       SEGINI COCOUR
  230.       COCOUR.NOMCOM=NCOM
  231.       COCOUR.LTREF=.TRUE.
  232.       SEGDES COCOUR
  233.       MYCOMS.LISCOM(**)=COCOUR
  234. *
  235. * Loi de comportement MADSMID
  236. * max |d|/min|d| où d=det J
  237. * avec un signe moins si d change de signe
  238. *
  239.       NCOM=BLANC
  240.       NCOM='MADSMID'
  241.       NCOCOF=0
  242.       SEGINI COCOUR
  243.       COCOUR.NOMCOM=NCOM
  244.       COCOUR.LTREF=.TRUE.
  245.       SEGDES COCOUR
  246.       MYCOMS.LISCOM(**)=COCOUR
  247. *
  248. * Loi de comportement AHUF : Huang
  249. *                     AHPF : Huang préconditionné
  250. *
  251.       DO ICOPO=1,2
  252.          IF (ICOPO.EQ.1) THEN
  253.             CIH='U'
  254.          ELSE
  255.             CIH='P'
  256.          ENDIF
  257.          NARG=2
  258.          NCOM=BLANC
  259.          NCOM='AH?F'
  260.          NCOM(3:3)=CIH
  261.          NCOCOF=((IDIM*(IDIM+1))/2)+NARG
  262.          SEGINI COCOUR
  263.          COCOUR.NOMCOM=NCOM
  264. * Plus clair, mais inutile
  265. *      DO ICOCOF=1,NCOCOF
  266. *         COCOUR.DERCOF(ICOCOF)=0
  267. *      ENDDO
  268.          COCOUR.LTREF=.TRUE.
  269.          SEGDES COCOUR
  270.          MYCOMS.LISCOM(**)=COCOUR
  271. *
  272. * Lois de comportement AHUR{1..idim}*{0..iesder}
  273. *
  274.          DO IDIM1=1,IDIM
  275.             CALL INT2CH(IDIM1,CDIM1,IMPR,IRET)
  276.             IF (IRET.NE.0) GOTO 9999
  277.             DO IDIM2=0,IESDER
  278.                CALL INT2CH(IDIM2,CDIM2,IMPR,IRET)
  279.                IF (IRET.NE.0) GOTO 9999
  280.                NCOM=BLANC
  281.                NCOM='AH?R'
  282.                NCOM(3:3)=CIH
  283.                NCOM(5:5)=CDIM1
  284.                NCOM(6:6)=CDIM2
  285.                NCOCOF=((IDIM*(IDIM+1))/2)+NARG
  286.                SEGINI COCOUR
  287.                COCOUR.NOMCOM=NCOM
  288.                COCOUR.LTREF=.TRUE.
  289.                SEGDES COCOUR
  290.                MYCOMS.LISCOM(**)=COCOUR
  291.             ENDDO
  292.          ENDDO
  293. *
  294. * Lois de comportement AHUJ{1..idim}*4
  295. *
  296.          DO IDIM1=1,IDIM
  297.             CALL INT2CH(IDIM1,CDIM1,IMPR,IRET)
  298.             IF (IRET.NE.0) GOTO 9999
  299.             DO IDIM2=0,IESDER
  300.                CALL INT2CH(IDIM2,CDIM2,IMPR,IRET)
  301.                IF (IRET.NE.0) GOTO 9999
  302.                DO IDIM3=1,IDIM
  303.                   CALL INT2CH(IDIM3,CDIM3,IMPR,IRET)
  304.                   IF (IRET.NE.0) GOTO 9999
  305.                   DO IDIM4=0,IESDER
  306.                      CALL INT2CH(IDIM4,CDIM4,IMPR,IRET)
  307.                      IF (IRET.NE.0) GOTO 9999
  308.                      NCOM=BLANC
  309.                      NCOM='AH?J'
  310.                      NCOM(3:3)=CIH
  311.                      NCOM(5:5)=CDIM1
  312.                      NCOM(6:6)=CDIM2
  313.                      NCOM(7:7)=CDIM3
  314.                      NCOM(8:8)=CDIM4
  315.                      NCOCOF=((IDIM*(IDIM+1))/2)+NARG
  316.                      SEGINI COCOUR
  317.                      COCOUR.NOMCOM=NCOM
  318.                      COCOUR.LTREF=.TRUE.
  319.                      SEGDES COCOUR
  320.                      MYCOMS.LISCOM(**)=COCOUR
  321.                   ENDDO
  322.                ENDDO
  323.             ENDDO
  324.          ENDDO
  325.       ENDDO
  326. *
  327. * Loi de comportement Qualité de maillage
  328. *
  329.       DO IQ=1,2
  330.          NCOM=BLANC
  331.          IF (IQ.EQ.1) THEN
  332.             NCOM='QEQU'
  333.          ELSE
  334.             NCOM='QALI'
  335.          ENDIF
  336.          NCOCOF=((IDIM*(IDIM+1))/2)
  337.          SEGINI COCOUR
  338.          COCOUR.NOMCOM=NCOM
  339.          COCOUR.LTREF=.TRUE.
  340.          SEGDES COCOUR
  341.          MYCOMS.LISCOM(**)=COCOUR
  342.       ENDDO
  343. *
  344. * Lois de comportement MUSTAB{1,2,3}{1..idim}
  345. * (viscosité numérique)
  346. *
  347.       DO IMETH=1,3
  348.          CALL INT2CH(IMETH,CMETH,IMPR,IRET)
  349.          IF (IRET.NE.0) GOTO 9999
  350.          DO IDIM1=0,IDIM
  351.             CALL INT2CH(IDIM1,CDIM1,IMPR,IRET)
  352.             IF (IRET.NE.0) GOTO 9999
  353.             NCOM=BLANC
  354.             NCOM='MUSTAB'
  355.             NCOM(7:7)=CMETH
  356.             NCOM(8:8)=CDIM1
  357.             NCOCOF=IDIM+3
  358.             SEGINI COCOUR
  359.             COCOUR.NOMCOM=NCOM
  360.             COCOUR.LTREF=.TRUE.
  361.             SEGDES COCOUR
  362.             MYCOMS.LISCOM(**)=COCOUR
  363.          ENDDO
  364.       ENDDO
  365. *
  366. * Loi de comportement TSUF : Tension de surface
  367. *
  368.       NCOM=BLANC
  369.       NCOM='TSUF'
  370.       NCOCOF=1
  371.       SEGINI COCOUR
  372.       COCOUR.NOMCOM=NCOM
  373.       COCOUR.LTREF=.TRUE.
  374.       SEGDES COCOUR
  375.       MYCOMS.LISCOM(**)=COCOUR
  376. *
  377. * Lois de comportement TSUR{1..idim}*{0..iesder}
  378. *
  379.       DO IDIM1=0,IDIM
  380.          CALL INT2CH(IDIM1,CDIM1,IMPR,IRET)
  381.          IF (IRET.NE.0) GOTO 9999
  382.          DO IDIM2=0,IESDER
  383.             CALL INT2CH(IDIM2,CDIM2,IMPR,IRET)
  384.             IF (IRET.NE.0) GOTO 9999
  385.             NCOM=BLANC
  386.             NCOM='TSUR'
  387.             NCOM(5:5)=CDIM1
  388.             NCOM(6:6)=CDIM2
  389.             NCOCOF=1
  390.             SEGINI COCOUR
  391.             COCOUR.NOMCOM=NCOM
  392.             COCOUR.LTREF=.TRUE.
  393.             SEGDES COCOUR
  394.             MYCOMS.LISCOM(**)=COCOUR
  395.          ENDDO
  396.       ENDDO
  397. *
  398. * Lois de comportement TSUJ{1..idim}*4
  399. *
  400.       DO IDIM1=1,IDIM
  401.          CALL INT2CH(IDIM1,CDIM1,IMPR,IRET)
  402.          IF (IRET.NE.0) GOTO 9999
  403.          DO IDIM2=0,IESDER
  404.             CALL INT2CH(IDIM2,CDIM2,IMPR,IRET)
  405.             IF (IRET.NE.0) GOTO 9999
  406.             DO IDIM3=1,IDIM
  407.                CALL INT2CH(IDIM3,CDIM3,IMPR,IRET)
  408.                IF (IRET.NE.0) GOTO 9999
  409.                DO IDIM4=0,IESDER
  410.                   CALL INT2CH(IDIM4,CDIM4,IMPR,IRET)
  411.                   IF (IRET.NE.0) GOTO 9999
  412.                   NCOM=BLANC
  413.                   NCOM='TSUJ'
  414.                   NCOM(5:5)=CDIM1
  415.                   NCOM(6:6)=CDIM2
  416.                   NCOM(7:7)=CDIM3
  417.                   NCOM(8:8)=CDIM4
  418.                   NCOCOF=1
  419.                   SEGINI COCOUR
  420.                   COCOUR.NOMCOM=NCOM
  421.                   COCOUR.LTREF=.TRUE.
  422.                   SEGDES COCOUR
  423.                   MYCOMS.LISCOM(**)=COCOUR
  424.                ENDDO
  425.             ENDDO
  426.          ENDDO
  427.       ENDDO
  428. *
  429. * Lois de comportement TSU{1..4}{1..idim}*4
  430. *
  431.       DO IDIM0=1,4
  432.          CALL INT2CH(IDIM0,CDIM0,IMPR,IRET)
  433.          IF (IRET.NE.0) GOTO 9999
  434.          DO IDIM1=1,IDIM
  435.             CALL INT2CH(IDIM1,CDIM1,IMPR,IRET)
  436.             IF (IRET.NE.0) GOTO 9999
  437.             DO IDIM2=0,IESDER
  438.                CALL INT2CH(IDIM2,CDIM2,IMPR,IRET)
  439.                IF (IRET.NE.0) GOTO 9999
  440.                DO IDIM3=1,IDIM
  441.                   CALL INT2CH(IDIM3,CDIM3,IMPR,IRET)
  442.                   IF (IRET.NE.0) GOTO 9999
  443.                   DO IDIM4=0,IESDER
  444.                      CALL INT2CH(IDIM4,CDIM4,IMPR,IRET)
  445.                      IF (IRET.NE.0) GOTO 9999
  446.                      NCOM=BLANC
  447.                      NCOM='TSU'
  448.                      NCOM(4:4)=CDIM0
  449.                      NCOM(5:5)=CDIM1
  450.                      NCOM(6:6)=CDIM2
  451.                      NCOM(7:7)=CDIM3
  452.                      NCOM(8:8)=CDIM4
  453.                      NCOCOF=1
  454.                      SEGINI COCOUR
  455.                      COCOUR.NOMCOM=NCOM
  456.                      COCOUR.LTREF=.TRUE.
  457.                      SEGDES COCOUR
  458.                      MYCOMS.LISCOM(**)=COCOUR
  459.                   ENDDO
  460.                ENDDO
  461.             ENDDO
  462.          ENDDO
  463.       ENDDO
  464. *
  465. * Lois de comportement VNOR{1..idim}
  466. *
  467.       DO IDIM1=1,IDIM
  468.          CALL INT2CH(IDIM1,CDIM1,IMPR,IRET)
  469.          IF (IRET.NE.0) GOTO 9999
  470.          NCOM=BLANC
  471.          NCOM='VNOR'
  472.          NCOM(5:5)=CDIM1
  473.          NCOCOF=0
  474.          SEGINI COCOUR
  475.          COCOUR.NOMCOM=NCOM
  476.          COCOUR.LTREF=.TRUE.
  477.          SEGDES COCOUR
  478.          MYCOMS.LISCOM(**)=COCOUR
  479.       ENDDO
  480. *
  481. * Lois de comportement VNOJ{1..idim}*3
  482. *
  483.       DO IDIM1=1,IDIM
  484.          CALL INT2CH(IDIM1,CDIM1,IMPR,IRET)
  485.          IF (IRET.NE.0) GOTO 9999
  486.          DO IDIM2=0,IDIM
  487.             CALL INT2CH(IDIM2,CDIM2,IMPR,IRET)
  488.             IF (IRET.NE.0) GOTO 9999
  489.             DO IDIM3=1,IDIM
  490.                CALL INT2CH(IDIM3,CDIM3,IMPR,IRET)
  491.                IF (IRET.NE.0) GOTO 9999
  492.                NCOM=BLANC
  493.                NCOM='VNOJ'
  494.                NCOM(5:5)=CDIM1
  495.                NCOM(6:6)=CDIM2
  496.                NCOM(7:7)=CDIM3
  497.                NCOCOF=0
  498.                SEGINI COCOUR
  499.                COCOUR.NOMCOM=NCOM
  500.                COCOUR.LTREF=.TRUE.
  501.                SEGDES COCOUR
  502.                MYCOMS.LISCOM(**)=COCOUR
  503.             ENDDO
  504.          ENDDO
  505.       ENDDO
  506. *
  507. * Lois de comportement IMET{1..idim}{idim1..idim}
  508. *
  509.       DO IDIM1=1,IDIM
  510.          CALL INT2CH(IDIM1,CDIM1,IMPR,IRET)
  511.          IF (IRET.NE.0) GOTO 9999
  512.          DO IDIM2=IDIM1,IDIM
  513.             CALL INT2CH(IDIM2,CDIM2,IMPR,IRET)
  514.             IF (IRET.NE.0) GOTO 9999
  515.             NCOM=BLANC
  516.             NCOM='IMET'
  517.             NCOM(5:5)=CDIM1
  518.             NCOM(6:6)=CDIM2
  519.             NCOCOF=0
  520.             SEGINI COCOUR
  521.             COCOUR.NOMCOM=NCOM
  522.             COCOUR.LTREF=.TRUE.
  523.             SEGDES COCOUR
  524.             MYCOMS.LISCOM(**)=COCOUR
  525.          ENDDO
  526.       ENDDO
  527. *
  528. *     Loi de comportement X1 : première coordonnée des noeuds
  529. *     Uniquement si mode AXI car sinon JPc non définie (cf. nlin.eso)      
  530. *
  531.       IF (IFOMOD.NE.-1) THEN
  532.          NCOM=BLANC
  533.          NCOM='X1'
  534.          NCOCOF=0
  535.          SEGINI COCOUR
  536.          COCOUR.NOMCOM=NCOM
  537.          COCOUR.LTREF=.TRUE.
  538.          SEGDES COCOUR
  539.          MYCOMS.LISCOM(**)=COCOUR
  540.       ENDIF
  541. *
  542. * Impression finale
  543. *
  544.       IF (IMPR.GT.1) THEN
  545.          SEGPRT,MYCOMS
  546.          NBCOMP=MYCOMS.LISCOM(/1)
  547.          DO IBCOMP=1,NBCOMP
  548.             COCOUR=MYCOMS.LISCOM(IBCOMP)
  549.             SEGPRT,COCOUR
  550.          ENDDO
  551.       ENDIF
  552.       SEGDES MYCOMS
  553. *
  554. * Normal termination
  555. *
  556.       IRET=0
  557.       RETURN
  558. *
  559. * Format handling
  560. *
  561. *
  562. * Error handling
  563. *
  564.  9999 CONTINUE
  565.       IRET=1
  566.       WRITE(IOIMP,*) 'An error was detected in subroutine incoms'
  567.       RETURN
  568. *
  569. * End of subroutine INCOMS
  570. *
  571.       END
  572.  
  573.  
  574.  
  575.  
  576.  
  577.  
  578.  
  579.  
  580.  
  581.  
  582.  

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