Télécharger la plaquette Cast3M


fr en

Télécharger motan.eso

Retour à la liste

Numérotation des lignes : 

motan

  1. C MOTAN     SOURCE    PV        20/04/10    12:45:04     10579          
  2.       SUBROUTINE MOTAN
  3. *_______________________________________________________________________
  4. *
  5. *           Operateur permettant de calculer le module tangent
  6. *
  7. *   SYNTAXE :
  8. *   ---------
  9. *
  10. *    ETAN=MOTAN MOD1 SI1 VA1 ( CA1 ) (PRECIS) ;
  11. *
  12. *        PRECIS  Flottant (par defaut 1.D-3)
  13. *        MOD1    Modele de calcul, type MMODEL
  14. *        SI1     Champ par element de contraintes,type MCHAML
  15. *        VA1     Champ par element de variables internes,type MCHAML
  16. *        CA1     Champ par element de caracteristiques (materielle et/ou
  17. *                geometriques),type MCHAML
  18. *        ETAN    Champ par element de type scalaire contenant les
  19. *                modules tangent aux points des elements,type MCHAML
  20. *
  21. *    CODE L.EBERSOLT NOVEMBRE 86
  22. *
  23. *    passage aux nouveaux CHAMELEMs par JM CAMPENON LE 06/91
  24. *_______________________________________________________________________
  25. *
  26.       IMPLICIT INTEGER(I-N)
  27.       IMPLICIT REAL*8(A-H,O-Z)
  28. *
  29.  
  30. -INC PPARAM
  31. -INC CCOPTIO
  32. *
  33. *     Lecture d'un modele
  34. *
  35.       CALL LIROBJ('MMODEL  ',IPMODL,1,IRTM)
  36.       CALL LIROBJ('MMODEL  ',IPMODL,1,IRETOU)
  37.       IF (IERR.NE.0) RETURN
  38. *
  39. *     Lecture du MCHAML de CONTRAINTES
  40. *
  41.       CALL LIROBJ('MCHAML  ',IPIN,1,IRT1)
  42.       CALL ACTOBJ('MCHAML  ',IPIN,1)
  43.       IF (IERR.NE.0) RETURN
  44.       CALL REDUAF(IPIN,IPMODL,IPCHE1,0,IR,KER)
  45.       IF(IR   .NE. 1) CALL ERREUR(KER)
  46.       IF(IERR .NE. 0) RETURN
  47. *
  48. *     Lecture du MCHAML de VARIABLES INTERNES
  49. *
  50.       CALL LIROBJ('MCHAML',IPIN,1,IRT2)
  51.       CALL ACTOBJ('MCHAML',IPIN,1)
  52.       IF (IERR.NE.0) RETURN
  53.       CALL REDUAF(IPIN,IPMODL,IPCHE2,0,IR,KER)
  54.       IF(IR   .NE. 1) CALL ERREUR(KER)
  55.       IF(IERR .NE. 0) RETURN
  56. *
  57. *     Lecture du MCHAML de CARACTERISTIQUES
  58. *
  59.       CALL LIROBJ('MCHAML',IPIN,1,IRT3)
  60.       CALL ACTOBJ('MCHAML',IPIN,1)
  61.       IF (IERR.NE.0) RETURN
  62.       CALL REDUAF(IPIN,IPMODL,IPCHE3,0,IR,KER)
  63.       IF(IR   .NE. 1) CALL ERREUR(KER)
  64.       IF(IERR .NE. 0) RETURN
  65. *
  66. *     Lecture eventuelle de la precision des iteration internes
  67. *
  68.       CALL LIRREE(PRECIS,0,IRETPR)
  69.       IF (IERR.NE.0) RETURN
  70.       IF (IRETPR.EQ.0) PRECIS = 0.9
  71. *
  72.       CALL MOTANA(IPMODL,IPCHE1,IPCHE2,IPCHE3,PRECIS,IPSCAL,IRET)
  73. *
  74. *     Ecriture de la matrice de raideur tangente
  75. *
  76.       IF (IRET.EQ.1) THEN
  77.         CALL ACTOBJ('MCHAML  ',IPSCAL,1)
  78.         CALL ECROBJ('MCHAML  ',IPSCAL)
  79.       ENDIF
  80.  
  81.       END
  82.  
  83.  
  84.  
  85.  

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