Télécharger la plaquette Cast3M


fr en

Télécharger hhoidc.eso

Retour à la liste

Numérotation des lignes : 

hhoidc

  1. C HHOIDC    SOURCE    OF166741  24/05/06    21:15:10     11082          
  2. C HHOidc    SOURCE    FANDEUR                                           
  3.  
  4. C== FORMULATION HHO ====================================================
  5.  
  6.       SUBROUTINE HHOidc (IPMODL,IPNOCA)
  7.  
  8.       IMPLICIT INTEGER(I-N)
  9.       IMPLICIT REAL*8(A-H,O-Z)
  10.  
  11. -INC PPARAM
  12. -INC CCOPTIO
  13. -INC CCHHOPA
  14.  
  15. -INC SMMODEL
  16.  
  17.       imodel = IPMODL
  18. c*      segact,imodel*MOD   <- Actif en E/S
  19.  
  20.       MELE = imodel.nefmod
  21.       IF (MELE .NE. HHO_NUM_ELEMENT) RETURN
  22.  
  23.       mfr0 = nummfr(mele)
  24.       if (mfr0 .ne. HHO_MFR_ELEMENT) then
  25.         write(ioimp,*) 'HHOIDC : MFR0 incorrect !'
  26.         call erreur(5)
  27.         return
  28.       end if
  29. c-dbg        CALL HHONOB(IPMODL, nobHHO, iret)
  30. c-dbg        IF (nobHHO.LE.0)THEN
  31. c-dbg          write(ioimp,*) 'HHOIDC : IPMODL incorrect (not HHO)'
  32. c-dbg          CALL ERREUR(5)
  33. c-dbg          RETURN
  34. c-dbg        END IF
  35.  
  36. C*=====================================================================*
  37. C*= FORMULATION HHO 'MECANIQUE'                                       =*
  38. C*= On ajoute aux composantes materiaux                               =*
  39. C*=   - le coefficient de stabilisation        'STAB' (cellule)       =*
  40. C*=   - le poids de chaque point d'integration 'PIHO' (PG & cellule)  =*
  41. C*=   - la matrice B pour la formulation HHO   'BHHO' (PG & cellule)  =*
  42. C*=   - la matrice de stabilisation            'SHHO' (cellule)       =*
  43. C*=====================================================================*
  44.       IF (imodel.FORMOD(1).EQ.'MECANIQUE') THEN
  45.  
  46.         IF (IPNOCA.EQ.0) IPNOCA = imodel.lnomid(6)
  47.  
  48.         nomid = IPNOCA
  49.         SEGACT,nomid
  50.         nbrob0 = nomid.lesobl(/2)
  51.         nbrfa0 = nomid.lesfac(/2)
  52.         nbrobl = nbrob0 + 4
  53.         nbrfac = nbrfa0
  54. C* Ou facultative ?        nbrobl = nbrob0
  55. C* Ou facultative ?        nbrfac = nbrfac + 4
  56.         SEGADJ,nomid
  57.         nomid.LESOBL(nbrob0+1) = 'STAB'
  58.         nomid.LESOBL(nbrob0+2) = 'PIHO'
  59.         nomid.LESOBL(nbrob0+3) = 'BHHO'
  60.         nomid.LESOBL(nbrob0+4) = 'SHHO'
  61. C* Ou facultative ?        nomid.LESFAC(nbrfa0+1) = 'STAB'
  62. C* Ou facultative ?        nomid.LESFAC(nbrfa0+2) = 'PIHO'
  63. C* Ou facultative ?        nomid.LESFAC(nbrfa0+3) = 'BHHO'
  64. C* Ou facultative ?        nomid.LESFAC(nbrfa0+4) = 'SHHO'
  65. c*        SEGDES,nomid   <- Actif en S
  66.         IPNOCA = nomid
  67.  
  68. C*=====================================================================*
  69. C*= Autres formulations HHO que 'MECANIQUE'                           =*
  70. C*=====================================================================*
  71.       ELSE
  72.       END IF
  73.  
  74. c*      RETURN
  75.       END
  76.  
  77.  
  78.  

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