Télécharger model2.eso

Retour à la liste

Numérotation des lignes :

  1. C MODEL2 SOURCE KICH 18/12/05 21:15:09 10024
  2.  
  3. C=======================================================================
  4. C= Objet MODELE - Cas de la formulation 'MECANIQUE' =
  5. C= =
  6. C= Definition des noms de modeles de MATERIAU =
  7. C= et des noms de types d'ELEMENTS FINIS =
  8. C=======================================================================
  9.  
  10. SUBROUTINE MODEL2 (LESPRO,NBPRO,LESTEF,NBTEF,N1MAX,N2MAX)
  11.  
  12. IMPLICIT INTEGER(I-N)
  13. IMPLICIT REAL*8 (A-H,O-Z)
  14.  
  15. -INC CCOPTIO
  16.  
  17. CHARACTER*16 MOMODL(200)
  18. CHARACTER*16 LESPRO(*)
  19. CHARACTER*4 LESTEF(*)
  20.  
  21. C 1 - Noms des comportements LINEAIRES
  22. C --------------------------------------
  23. CALL MODLIN(MOMODL,NMOD)
  24. DO i=1,NMOD
  25. LESPRO(i)=MOMODL(i)
  26. ENDDO
  27. NBPRO=NMOD
  28.  
  29. C 2 - Noms des MATERIAUx pour le comportement LINEAIRE
  30. C ------------------------------------------------------
  31. CALL MODELA(MOMODL,NMOD)
  32. DO i=NBPRO+1,NBPRO+NMOD
  33. LESPRO(i)=MOMODL(i-NBPRO)
  34. ENDDO
  35. NBPRO=NBPRO+NMOD
  36.  
  37. C 3 - Noms des comportements NON LINEAIRES
  38. C ------------------------------------------
  39. CALL MODNLI(MOMODL,NMOD)
  40. DO i=NBPRO+1,NBPRO+NMOD
  41. LESPRO(i)=MOMODL(i-NBPRO)
  42. ENDDO
  43. NBPRO=NBPRO+NMOD
  44.  
  45. C 4 - Noms des modeles PLASTIQUES
  46. C ---------------------------------
  47. CALL MODPLA(MOMODL,NMOD)
  48. DO i=NBPRO+1,NBPRO+NMOD
  49. LESPRO(i)=MOMODL(i-NBPRO)
  50. ENDDO
  51. NBPRO=NBPRO+NMOD
  52.  
  53. C 5 - Noms des modeles de FLUAGE
  54. C --------------------------------
  55. CALL MODFLU(MOMODL,NMOD)
  56. DO i=NBPRO+1,NBPRO+NMOD
  57. LESPRO(i)=MOMODL(i-NBPRO)
  58. ENDDO
  59. NBPRO=NBPRO+NMOD
  60.  
  61. C 6 - Noms des modeles VISCOPLASTIQUES
  62. C --------------------------------------
  63. CALL MODVIS(MOMODL,NMOD)
  64. DO i=NBPRO+1,NBPRO+NMOD
  65. LESPRO(i)=MOMODL(i-NBPRO)
  66. ENDDO
  67. NBPRO=NBPRO+NMOD
  68.  
  69. C 7 - Noms des modeles ENDOMMAGEABLES
  70. C -------------------------------------
  71. CALL MODEND(MOMODL,NMOD)
  72. DO i=NBPRO+1,NBPRO+NMOD
  73. LESPRO(i)=MOMODL(i-NBPRO)
  74. ENDDO
  75. NBPRO=NBPRO+NMOD
  76.  
  77. C 8 - Noms des modeles PLASTIQUES ENDOMMAGEABLES
  78. C ------------------------------------------------
  79. CALL MODPLE(MOMODL,NMOD)
  80. DO i=NBPRO+1,NBPRO+NMOD
  81. LESPRO(i)=MOMODL(i-NBPRO)
  82. ENDDO
  83. NBPRO=NBPRO+NMOD
  84.  
  85. C 9 - Noms des modeles ELASTIQUES NON LINEAIRES
  86. C -----------------------------------------------
  87. CALL MODENL(MOMODL,NMOD)
  88. DO i=NBPRO+1,NBPRO+NMOD
  89. LESPRO(i)=MOMODL(i-NBPRO)
  90. ENDDO
  91. NBPRO=NBPRO+NMOD
  92.  
  93. C 10 - Noms des modeles VISCOPLASTIQUES EXTERNES
  94. C ------------------------------------------------
  95. CALL MODVIX(MOMODL,NMOD)
  96. DO i=NBPRO+1,NBPRO+NMOD
  97. LESPRO(i)=MOMODL(i-NBPRO)
  98. ENDDO
  99. NBPRO=NBPRO+NMOD
  100.  
  101. C 12 - Noms des modeles LIAISONS
  102. C ------------------------------------------------
  103. CALL MODLIA(MOMODL,NMOD)
  104. DO i=NBPRO+1,NBPRO+NMOD
  105. LESPRO(i)=MOMODL(i-NBPRO)
  106. ENDDO
  107. NBPRO=NBPRO+NMOD
  108.  
  109. C 13 - Noms des modeles IMPEDANCE
  110. C ------------------------------------------------
  111. CALL MODIMP(MOMODL,NMOD)
  112. DO i=NBPRO+1,NBPRO+NMOD
  113. LESPRO(i)=MOMODL(i-NBPRO)
  114. ENDDO
  115. NBPRO=NBPRO+NMOD
  116.  
  117. C CRItere Pression
  118. LESPRO(NBPRO + 1) = 'CRIP'
  119. NBPRO = NBPRO + 1
  120.  
  121. IF (NBPRO.GT.N1MAX) GOTO 99
  122.  
  123. C 11 - Noms de types d'ELEMENTS FINIS
  124. C -------------------------------------
  125. C DIMENSION 1 :
  126. C ---------------
  127. IF (IDIM.EQ.1) THEN
  128. NBTEF=2
  129. IF (NBTEF.GT.N2MAX) GOTO 99
  130. LESTEF(1)='M1D2'
  131. LESTEF(2)='M1D3'
  132. C DIMENSION 2 :
  133. C ---------------
  134. C Les EF sont definis suivant la valeur de IFOUR
  135. ELSE IF (IDIM.EQ.2) THEN
  136. IF (IFOUR.EQ.0.OR.IFOUR.EQ.1) THEN
  137. NBTEF=19
  138. IF (NBTEF.GT.N2MAX) GOTO 99
  139. LESTEF(1)='TRI3'
  140. LESTEF(2)='TRI6'
  141. LESTEF(3)='QUA4'
  142. LESTEF(4)='QUA8'
  143. LESTEF(5)='COQ2'
  144. LESTEF(6)='BARR'
  145. LESTEF(7)='TRIH'
  146. LESTEF(8)='CERC'
  147. LESTEF(9)='JOI2'
  148. LESTEF(10)='JOI3'
  149. LESTEF(11)='POLY'
  150. LESTEF(12)='BAR3'
  151. LESTEF(13)='QUAH'
  152. LESTEF(14)='Q4RI'
  153. LESTEF(15)='Q8RI'
  154. LESTEF(16)='XQ4R'
  155. LESTEF(17)='ZCO2'
  156. LESTEF(18)='COS2'
  157. LESTEF(19)='COA2'
  158. IF (IFOUR.EQ.0) THEN
  159. C Les elements incompressibles ne sont pas utilisables pour le calcul
  160. C en serie de FOURIER
  161. NBTEF=22
  162. IF (NBTEF.GT.N2MAX) GOTO 99
  163. LESTEF(18)='ICT3'
  164. LESTEF(19)='ICQ4'
  165. LESTEF(20)='ICT6'
  166. LESTEF(21)='ICQ8'
  167. LESTEF(22)='BBAR'
  168. ENDIF
  169. ELSE
  170. NBTEF=37
  171. IF (NBTEF.GT.N2MAX) GOTO 99
  172. LESTEF(1)='TRI3'
  173. LESTEF(2)='TRI6'
  174. LESTEF(3)='QUA4'
  175. LESTEF(4)='QUA8'
  176. LESTEF(5)='COQ2'
  177. LESTEF(6)='BARR'
  178. LESTEF(7)='TRIH'
  179. LESTEF(8)='JOI2'
  180. LESTEF(9)='POI1'
  181. LESTEF(10)='JOI3'
  182. LESTEF(11)='ICT3'
  183. LESTEF(12)='ICQ4'
  184. LESTEF(13)='ICT6'
  185. LESTEF(14)='ICQ8'
  186. LESTEF(15)='POLY'
  187. LESTEF(16)='BAR3'
  188. LESTEF(17)='QUAH'
  189. LESTEF(18)='TIMO'
  190. LESTEF(19)='POUT'
  191. LESTEF(20)='TRIS'
  192. LESTEF(21)='QUAS'
  193. LESTEF(22)='POIS'
  194. LESTEF(23)='TRH6'
  195. LESTEF(24)='SEGS'
  196. LESTEF(25)='POJS'
  197. LESTEF(26)='Q4RI'
  198. LESTEF(27)='Q8RI'
  199. LESTEF(28)='JCT3'
  200. LESTEF(29)='JCI4'
  201. LESTEF(30)='JGI2'
  202. LESTEF(31)='CIFL'
  203. LESTEF(32)='XQ4R'
  204. LESTEF(33)='JOI1'
  205. LESTEF(34)='ZCO2'
  206. LESTEF(35)='COS2'
  207. LESTEF(36)='COA2'
  208. LESTEF(37)='BBAR'
  209. ENDIF
  210. C DIMENSION 3 :
  211. C ---------------
  212. ELSE IF (IDIM.EQ.3) THEN
  213. NBTEF=54
  214. IF (NBTEF.GT.N2MAX) GOTO 99
  215. LESTEF( 1)='CUB8'
  216. LESTEF( 2)='CU20'
  217. LESTEF( 3)='PRI6'
  218. LESTEF( 4)='PR15'
  219. LESTEF( 5)='TET4'
  220. LESTEF( 6)='TE10'
  221. LESTEF( 7)='PYR5'
  222. LESTEF( 8)='PY13'
  223. LESTEF( 9)='COQ3'
  224. LESTEF(10)='COQ4'
  225. LESTEF(11)='COQ6'
  226. LESTEF(12)='COQ8'
  227. LESTEF(13)='DKT '
  228. LESTEF(14)='BARR'
  229. LESTEF(15)='POUT'
  230. LESTEF(16)='TUYA'
  231. LESTEF(17)='LISP'
  232. LESTEF(18)='TUFI'
  233. LESTEF(19)='DST '
  234. LESTEF(20)='JOI4'
  235. LESTEF(21)='TUYO'
  236. LESTEF(22)='TIMO'
  237. LESTEF(23)='TRIS'
  238. LESTEF(24)='QUAS'
  239. LESTEF(25)='POIS'
  240. LESTEF(26)='LISM'
  241. LESTEF(27)='JOT3'
  242. LESTEF(28)='BAR3'
  243. LESTEF(29)='BAEX'
  244. LESTEF(30)='LIA2'
  245. LESTEF(31)='CUBH'
  246. LESTEF(32)='POJS'
  247. LESTEF(33)='Q4RI'
  248. LESTEF(34)='Q8RI'
  249. LESTEF(35)='JGT3'
  250. LESTEF(36)='JGI4'
  251. LESTEF(37)='SEGS'
  252. LESTEF(38)='SHB8'
  253. LESTEF(39)='POI1'
  254. LESTEF(40)='XC8R'
  255. LESTEF(41)='JOI1'
  256. LESTEF(42)='ZCO3'
  257. LESTEF(43)='ZCO4'
  258. LESTEF(44)='COS2'
  259. LESTEF(45)='COA2'
  260. LESTEF(46)='ICC8'
  261. LESTEF(47)='ICT4'
  262. LESTEF(48)='ICP6'
  263. LESTEF(49)='IC20'
  264. LESTEF(50)='IC10'
  265. LESTEF(51)='IC15'
  266. LESTEF(52)='BBAR'
  267. LESTEF(53)='ICY5'
  268. LESTEF(54)='IC13'
  269. ENDIF
  270. RETURN
  271.  
  272. C 12 - Traitement des ERREURS
  273. C -----------------------------
  274. 99 CALL ERREUR(5)
  275. RETURN
  276.  
  277. END
  278.  
  279.  
  280.  
  281.  
  282.  
  283.  
  284.  
  285.  
  286.  
  287.  
  288.  
  289.  
  290.  
  291.  
  292.  
  293.  
  294.  

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