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

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