Télécharger motan.eso

Retour à la liste

Numérotation des lignes :

  1. C MOTAN SOURCE CB215821 16/12/05 21:40:12 9237
  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. -INC CCOPTIO
  30. *
  31. * Lecture d'un modele
  32. *
  33. CALL LIROBJ('MMODEL',IPMODL,1,IRTM)
  34. IF (IERR.NE.0) RETURN
  35. *
  36. * Lecture du MCHAML de CONTRAINTES
  37. *
  38. CALL LIROBJ('MCHAML',IPIN,1,IRT1)
  39. IF (IERR.NE.0) RETURN
  40. CALL REDUAF(IPIN,IPMODL,IPCHE1,0,IR,KER)
  41. IF(IR .NE. 1) CALL ERREUR(KER)
  42. IF(IERR .NE. 0) RETURN
  43. *
  44. * Lecture du MCHAML de VARIABLES INTERNES
  45. *
  46. CALL LIROBJ('MCHAML',IPIN,1,IRT2)
  47. IF (IERR.NE.0) RETURN
  48. CALL REDUAF(IPIN,IPMODL,IPCHE2,0,IR,KER)
  49. IF(IR .NE. 1) CALL ERREUR(KER)
  50. IF(IERR .NE. 0) RETURN
  51. *
  52. * Lecture du MCHAML de CARACTERISTIQUES
  53. *
  54. CALL LIROBJ('MCHAML',IPIN,1,IRT3)
  55. IF (IERR.NE.0) RETURN
  56. CALL REDUAF(IPIN,IPMODL,IPCHE3,0,IR,KER)
  57. IF(IR .NE. 1) CALL ERREUR(KER)
  58. IF(IERR .NE. 0) RETURN
  59. *
  60. * Lecture eventuelle de la precision des iteration internes
  61. *
  62. CALL LIRREE(PRECIS,0,IRETPR)
  63. IF (IERR.NE.0) RETURN
  64. IF (IRETPR.EQ.0) PRECIS = 0.9
  65. *
  66. CALL MOTANA(IPMODL,IPCHE1,IPCHE2,IPCHE3,PRECIS,IPSCAL,IRET)
  67. *
  68. * Ecriture de la matrice de raideur tangente
  69. *
  70. IF (IRET.EQ.1) CALL ECROBJ('MCHAML',IPSCAL)
  71. *
  72. RETURN
  73. END
  74.  
  75.  
  76.  

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