Télécharger la plaquette Cast3M


fr en

Télécharger adtuy.eso

Retour à la liste

Numérotation des lignes : 

adtuy

  1. C ADTUY     SOURCE    CB215821  26/03/06    21:15:05     12485          
  2.       SUBROUTINE adtuy (NEF,IPMAIL,IPINTE,IMATE,IVAMAT,NMATR,
  3.      &                   IPMATR,NLIGR)
  4.       IMPLICIT INTEGER(I-N)
  5.       IMPLICIT REAL*8(A-H,O-Z)
  6.  
  7. -INC PPARAM
  8. -INC CCOPTIO
  9. -INC SMCHAML
  10. -INC SMCOORD
  11. -INC SMELEME
  12. -INC SMINTE
  13. -INC SMRIGID
  14.  
  15. -INC TMPTVAL
  16.  
  17.       SEGMENT,MMAT1
  18.         REAL*8 CEL(NBNN,NBNN),XE(3,NBNN)
  19.         REAL*8 SHP(6,NBNN)
  20.       ENDSEGMENT
  21.  
  22. C  1 - INITIALISATIONS ET VERIFICATIONS
  23. C ======================================
  24.       MELEME = IPMAIL
  25. c*      SEGACT,MELEME
  26.       NBNN   = NUM(/1)
  27.       NBELEM = NUM(/2)
  28. C =====
  29.       MINTE = IPINTE
  30. c*      SEGACT,MINTE
  31.       NBPGAU = MINTE.POIGAU(/1)
  32.  
  33. C =====
  34.       MPTVAL = IVAMAT
  35. c*      SEGACT,MPTVAL
  36. C =====
  37.       XMATRI = IPMATR
  38.  
  39. C =====
  40. C Initialisation des segments de travail
  41. C =====
  42.       IF (IFOMOD.EQ.1) THEN
  43.          NDIM = 3
  44.       ELSE
  45.          NDIM = IDIM
  46.       ENDIF
  47.       SEGINI,MMAT1
  48.  
  49. C  2 - BOUCLE SUR LES ELEMENTS DU MAILLAGE ELEMENTAIRE IPMAIL
  50. C ============================================================
  51.       DO IEL = 1, NBELEM
  52. *
  53. *       MISE A ZERO DU TABLEAU CEL
  54. *
  55.         CALL ZERO(CEL,NBNN,NBNN)
  56. *
  57. *       COORDONNEES DES NOEUDS DE L'ELEMENT IEL DANS LE REPERE GLOBAL
  58. *
  59.         CALL DOXE(XCOOR,IDIM,NBNN,NUM,IEL,XE)
  60.  
  61.         DO IGAU = 1, NBPGAU
  62.  
  63. C         CB215821 : Cas particulier en 1D pour le Jacobien
  64. C                    il se simplifie entre GRAD(Ni) et l'élément d'integration dL
  65.           XSCAL = poigau(igau)
  66.  
  67. *         Recuperation des valeurs des composantes (ordre dans tadve1.eso)
  68.           DO i = 1, NMATR
  69.             MELVAL = IVAL(i)
  70.             ibmn   = MIN(iel ,VELCHE(/2))
  71.             igmn   = MIN(igau,VELCHE(/1))
  72.             XSCAL  = XSCAL*VELCHE(igmn,ibmn)
  73.           ENDDO
  74.  
  75.           DO i=1,nbnn
  76.             cz = shptot(1,i,igau) * XSCAL
  77.             DO j=1,nbnn
  78.               cel(i,j)=cel(i,j) +cz*shptot(2,j,igau)
  79.             ENDDO
  80.           ENDDO
  81.         ENDDO
  82.  
  83.         CALL rempms(cel,nbnn,re(1,1,iel))
  84.       ENDDO
  85.  
  86.       SEGSUP,MMAT1
  87.       RETURN
  88.       END
  89.  
  90.  

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