Télécharger elftem.eso

Retour à la liste

Numérotation des lignes :

elftem
  1. C ELFTEM SOURCE CHAT 05/01/12 23:32:44 5004
  2. SUBROUTINE ELFTEM
  3. C
  4. C =====================================================================
  5. C
  6. C
  7. C SOLUT1 = ELFE
  8. C
  9. C / 'TEMPS'POUTRE
  10. C STRUC1 ATTAC1 TEMPS DT1 CHARG1 ( M1
  11. C GREEN ( NFOIS NN ) ;
  12. C
  13. C SOLUT1 : OBJET DE TYPE SOLUTION DYNAMIQUE
  14. C STRUC1 : OBJET DE TYPE STRUCTURE, ENSEMBLE DES ELEMENTS DE POUTRE
  15. C ATTAC1 : OBJET DE TYPE ATTACHE, LIAISONS ENTRE LES ELEMENTS
  16. C TEMPS : FLOTTANT, VALEUR DU TEMPS DE CALCUL
  17. C DT1 : FLOTTANT, VALEUR DU PAS DE TEMPS
  18. C CHARG1 : OBJET DE TYPE CHARGEMENT
  19. C GREEN : OBJET EVOLUTION CONTENANT LES FONCTIONS DE GREEN
  20. C FACULTATIVEMENT
  21. C M1 : ENTIER, INDICE DE TRONCATURE , SI 0 PAS DE TRONCATURE
  22. C NFOI : MOT CLE
  23. C NN : CHPOINT RESULTAT DANS L'OBJET SOLUTION TOUT LES NN PAS
  24. C
  25. C VERSION : JUIN 86
  26. C PROGRAMMEUR : GUILBAUD
  27. C MODIFIE PAR P.MANIGOT EN MARS 1988
  28. C MODIFIE PAR X.VACELET EN JUIN 1989
  29. C
  30. C =====================================================================
  31. C
  32. IMPLICIT INTEGER(I-N)
  33. IMPLICIT REAL*8(A-H,O-Z)
  34.  
  35. -INC PPARAM
  36. -INC CCOPTIO
  37. -INC SMCOORD
  38. CHARACTER*4 MFOIS(1)
  39. CHARACTER*40 CHA1
  40. DATA MFOIS/'NFOI'/
  41. C TAILLE DES BLOCS (NOMBRE DE REELS)
  42. DATA LBLO/50000/
  43. C PRECISION
  44. DATA EPS/1.D-11/
  45. M=0
  46. LFOIS=1
  47. C
  48. C 1 - LECTURE DES DONNES
  49. C
  50. CALL LIROBJ('STRUCTUR',MSTRUC,1,IRETOU)
  51. IF (IERR.NE.0) RETURN
  52. CALL LIROBJ ('ATTACHE ',MATTAC,1,IRETOU)
  53. IF (IERR.NE.0) RETURN
  54. C
  55. CALL LIRREE (TEMPS,1,IRETOU)
  56. IF (IERR.NE.0) RETURN
  57. IF (TEMPS.LE.0) THEN
  58. C *** VALEUR DE TEMPS NEGATIVE OU NULLE
  59. CALL ERREUR(406)
  60. ENDIF
  61. CALL LIRREE (DELTAT,1,IRETOU)
  62. IF (IERR.NE.0) RETURN
  63. IF (DELTAT.LE.0.D0) THEN
  64. C *** VALEUR DE PAS DE TEMPS NEGATIVE OU NULLE
  65. CALL ERREUR(406)
  66. ENDIF
  67. IF (TEMPS.LT.DELTAT) THEN
  68. TAMP=TEMPS
  69. TEMPS=DELTAT
  70. DELTAT=TAMP
  71. END IF
  72. N=NINT(TEMPS/DELTAT)
  73. C
  74. CALL LIRENT (M,0,IRETOU)
  75. IF (IRETOU.NE.0 .AND. M.LE.0) THEN
  76. C *** INDICE DE TRONCATURE NEGATIF OU NUL
  77. CALL ERREUR(407)
  78. ENDIF
  79. C
  80. CALL LIROBJ ('CHARGEME',MCHARG,1,IRETOU)
  81. C
  82. CALL LIROBJ ('EVOLUTIO',KGREEN,1,IRETOU)
  83. IF (IERR .NE. 0) RETURN
  84. C
  85. NFOIS=1
  86. CALL LIRMOT(MFOIS,LFOIS,IMOT,0)
  87. IF (IMOT.NE.0) THEN
  88. CALL LIRENT (NFOIS,1,IRETOU)
  89. END IF
  90. SEGACT MCOORD
  91. C
  92. C
  93. C 2 - TRAITEMENT DE L'OBJET STRUCTURE
  94. C
  95. CALL ELFELE(MSTRUC,KNREFE)
  96. IF (IERR.NE.0) RETURN
  97. C
  98. C
  99. C 3 - TRAITEMENT DE L'OBJET SOLUTION GREEN, CALCUL DES MATRICES ANBN
  100. C
  101. CALL ELFGR0 (KGREEN,DELTAT,N, KGREE2)
  102. IF (IERR.NE.0) RETURN
  103. MAXBLO=LBLO
  104. CALL ELFGR2(KGREE2,DELTAT,M,N,MAXBLO,KNREFE,KANBN,KDNCN)
  105. IF (IERR.NE.0) RETURN
  106. C
  107. C
  108. C 4 - TRAITEMENT DE L'OBJET ATTACHE, CALCUL DES MATRICES DE LIAISON
  109. C
  110. CALL ELFTRI(MATTAC,KNREFE,MATTA1)
  111. IF (IERR.NE.0) RETURN
  112. CALL ELFTR2(MATTA1,KNREFE)
  113. IF (IERR.NE.0) RETURN
  114. CALL ELFLIA(MATTA1,KANBN,KNREFE,KOMALI,EPS)
  115. IF (IERR.NE.0) RETURN
  116. C
  117. C
  118. C 5 - TRAITEMENT DU CHARGEMENT, CALCUL DES SECONDS-MEMBRES
  119. C
  120. CALL ELFCHA(MCHARG,KOMALI,KNREFE,KACHAR)
  121. IF (IERR.NE.0) RETURN
  122. C
  123. C
  124. C 6 - RESOLUTION ET CREATION DE L'OBJET SOLUTION
  125. C
  126. CALL ELFRES(KANBN,KDNCN,KNREFE,KOMALI,KACHAR,DELTAT,M,N,NFOIS,
  127. *KSOLUT)
  128. C
  129. IF (KGREEN .NE. KGREE2) CALL DTEVOL (KGREE2,1)
  130. C
  131. CALL ECROBJ('SOLUTION',KSOLUT)
  132. RETURN
  133. C
  134. END
  135.  
  136.  

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