Télécharger idendo.eso

Retour à la liste

Numérotation des lignes :

  1. C IDENDO SOURCE BR232186 16/12/05 21:39:41 9238
  2. SUBROUTINE IDENDO (MFR,MOOBL,MOFAC,IPLAC,IRET)
  3. *--------------------------------------------------------------------*
  4. * Noms de composantes de materiaux endommageables *
  5. *--------------------------------------------------------------------*
  6. * *
  7. * ENTREES: *
  8. * ------- *
  9. * MFR num{ro de formulation *
  10. * MOOBL pointeur sur les mots obligatoires *
  11. * MOFAC POINTEUR SUR LES MOTS FACULTATIFS *
  12. * IPLAC numero du materiau dans MODEND *
  13. * *
  14. * SORTIE : *
  15. * ------- *
  16. * *
  17. * IRET : 0 SI PROBLEME *
  18. * *
  19. *--------------------------------------------------------------------*
  20. *
  21. IMPLICIT INTEGER(I-N)
  22. -INC SMLMOTS
  23. -INC CCOPTIO
  24. *
  25. JGN=8
  26. IRET=1
  27. *
  28. *--------------------------------------------------------------------*
  29. * *
  30. * T Y P E E N D O M M A G E A B L E *
  31. * *
  32. *--------------------------------------------------------------------*
  33. IF (IPLAC.EQ.1) THEN
  34. *
  35. * MODELE ENDOMMAGEMENT SCALAIRE DE MAZARS
  36. *
  37. MLMOTS=MOOBL
  38. JGA=MOTS(/2)
  39. JGM=JGA+6
  40. SEGADJ MLMOTS
  41. MOOBL=MLMOTS
  42. MOTS(JGA+1)= 'KTR0'
  43. MOTS(JGA+2)= 'ACOM'
  44. MOTS(JGA+3)= 'BCOM'
  45. MOTS(JGA+4)= 'ATRA'
  46. MOTS(JGA+5)= 'BTRA'
  47. MOTS(JGA+6)= 'BETA'
  48. ELSE IF (IPLAC.EQ.2) THEN
  49. *
  50. * MODELE ENDOMMAGEMENT UNILATERAL LA BORDERIE
  51. *
  52. MLMOTS=MOOBL
  53. JGA=MOTS(/2)
  54. JGM=JGA+9
  55. SEGADJ MLMOTS
  56. MOOBL=MLMOTS
  57. MOTS(JGA+1)= 'YS1 '
  58. MOTS(JGA+2)= 'YS2 '
  59. MOTS(JGA+3)= 'A1 '
  60. MOTS(JGA+4)= 'A2 '
  61. MOTS(JGA+5)= 'B1 '
  62. MOTS(JGA+6)= 'B2 '
  63. MOTS(JGA+7)= 'BET1'
  64. MOTS(JGA+8)= 'BET2'
  65. MOTS(JGA+9)= 'SIGF'
  66. ELSE IF (IPLAC.EQ.3) THEN
  67. *
  68. * MODELE ROTATING_CRACK
  69. *
  70. MLMOTS=MOOBL
  71. JGA=MOTS(/2)
  72. JGM=JGA+2
  73. SEGADJ MLMOTS
  74. MOOBL=MLMOTS
  75. MOTS(JGA+1)= 'EPCR'
  76. MOTS(JGA+2)= 'MUP '
  77. ELSE IF (IPLAC.EQ.4) THEN
  78. *
  79. * MODELE SIC_SIC
  80. *
  81. MLMOTS=MOOBL
  82. JGA=MOTS(/2)
  83. JGM=JGA+15
  84. SEGADJ MLMOTS
  85. MOOBL=MLMOTS
  86. MOTS(JGA+1)= 'G1DC'
  87. MOTS(JGA+2)= 'G1Y0'
  88. MOTS(JGA+3)= 'G1YC'
  89. MOTS(JGA+4)= 'G1P'
  90. MOTS(JGA+5)= 'G2DC'
  91. MOTS(JGA+6)= 'G2Y0'
  92. MOTS(JGA+7)= 'G2YC'
  93. MOTS(JGA+8)= 'G2P'
  94. MOTS(JGA+9)= 'G3DC'
  95. MOTS(JGA+10)= 'G3Y0'
  96. MOTS(JGA+11)= 'G3YC'
  97. MOTS(JGA+12)= 'G3P'
  98. MOTS(JGA+13)= 'H23'
  99. MOTS(JGA+14)= 'H13'
  100. MOTS(JGA+15)= 'H12'
  101. MLMOTS=MOFAC
  102. JGA=MOTS(/2)
  103. JGM=JGA+3
  104. SEGADJ MLMOTS
  105. MOFAC=MLMOTS
  106. MOTS(JGA+1)= 'A12'
  107. MOTS(JGA+2)= 'A23'
  108. MOTS(JGA+3)= 'A13'
  109. *
  110. ELSE IF (IPLAC.EQ.5 .OR. IPLAC .EQ. 7) THEN
  111. *
  112. * MODELE ENDOMMAGEMENT MICROPLANS
  113. *
  114. MLMOTS=MOOBL
  115. JGA=MOTS(/2)
  116. JGM=JGA+5
  117. SEGADJ MLMOTS
  118. MOOBL=MLMOTS
  119. MOTS(JGA+1)= 'KTR0'
  120. MOTS(JGA+2)= 'BTRA'
  121. MOTS(JGA+3)= 'FPC '
  122. MOTS(JGA+4)= 'MP '
  123. MOTS(JGA+5)= 'ALF1'
  124. *
  125. ELSE IF (IPLAC.EQ.6) THEN
  126. *
  127. * MODELE ENDOMMAGEMENT UNILATERAL LA BORDERIE
  128. * + VISCOSITE J.F. DUBE
  129. *
  130. MLMOTS=MOOBL
  131. JGA=MOTS(/2)
  132. JGM=JGA+13
  133. SEGADJ MLMOTS
  134. MOOBL=MLMOTS
  135. MOTS(JGA+1)= 'YS1 '
  136. MOTS(JGA+2)= 'YS2 '
  137. MOTS(JGA+3)= 'A1 '
  138. MOTS(JGA+4)= 'A2 '
  139. MOTS(JGA+5)= 'B1 '
  140. MOTS(JGA+6)= 'B2 '
  141. MOTS(JGA+7)= 'BET1'
  142. MOTS(JGA+8)= 'BET2'
  143. MOTS(JGA+9)= 'SIGF'
  144. MOTS(JGA+10)= 'M1 '
  145. MOTS(JGA+11)= 'M2 '
  146. MOTS(JGA+12)= 'ALP1'
  147. MOTS(JGA+13)= 'ALP2'
  148. * BCN
  149. ELSE IF (IPLAC.EQ.8) THEN
  150. *
  151. * Modified Von mises model
  152. *
  153. MLMOTS=MOOBL
  154. JGA=MOTS(/2)
  155. JGM=JGA+5
  156. SEGADJ MLMOTS
  157. MOOBL=MLMOTS
  158. MOTS(JGA+1)= 'Y0 '
  159. MOTS(JGA+2)= 'B1 '
  160. MOTS(JGA+3)= 'B2 '
  161. MOTS(JGA+4)= 'RATI'
  162. MOTS(JGA+5)= 'LOI '
  163. * BCN
  164.  
  165. * RDj
  166. ELSE IF (IPLAC.EQ.9) THEN
  167. *
  168. * MODELE ENDOMMAGEMENT SCALAIRE AVEC DEFORMATION ANELASTIQUE
  169. * CRUSH_CRACK
  170. MLMOTS=MOOBL
  171. JGA=MOTS(/2)
  172. JGM=JGA+9
  173. SEGADJ MLMOTS
  174. MOOBL=MLMOTS
  175. MOTS(JGA+1)= 'KCO0'
  176. MOTS(JGA+2)= 'KTR0'
  177. MOTS(JGA+3)= 'ACRU'
  178. MOTS(JGA+4)= 'BCRU'
  179. MOTS(JGA+5)= 'ACRA'
  180. MOTS(JGA+6)= 'BCRA'
  181. MOTS(JGA+7)= 'ECOM'
  182. MOTS(JGA+8)= 'ECOB'
  183. MOTS(JGA+9)= 'ETRA'
  184. * RDj
  185.  
  186. * CCd
  187. ELSE IF (IPLAC.EQ.10) THEN
  188. *
  189. *+++++++++++++++++++++++++++++++++++
  190. * MODELE DAMAGE_TENSION&COMPRESSION = DAMAGE_TC
  191. *+++++++++++++++++++++++++++++++++++
  192. *
  193. * COMPULSORY CARACTERISTICS
  194. *
  195. MLMOTS=MOOBL
  196. JGA=MOTS(/2)
  197. JGM=JGA+15
  198. SEGADJ MLMOTS
  199. MOOBL=MLMOTS
  200. MOTS(JGA+1)='HLEN'
  201. MOTS(JGA+2)='GVAL'
  202. MOTS(JGA+3)='FTUL'
  203. MOTS(JGA+4)='REDC'
  204. MOTS(JGA+5)='FC01'
  205. MOTS(JGA+6)='RT45'
  206. MOTS(JGA+7)='FCU1'
  207. MOTS(JGA+8)='STRU'
  208. MOTS(JGA+9)='EXTP'
  209. MOTS(JGA+10)='STRP'
  210. MOTS(JGA+11)='EXT1'
  211. MOTS(JGA+12)='STR1'
  212. MOTS(JGA+13)='EXT2'
  213. MOTS(JGA+14)='STR2'
  214. MOTS(JGA+15)='NCRI'
  215. *
  216. * FACULTATIVE CARACTERISTICS
  217. *
  218. MLMOTS=MOFAC
  219. JGA=MOTS(/2)
  220. JGM=JGA+10
  221. SEGADJ MLMOTS
  222. MOFAC=MLMOTS
  223. C WARNING !!!!!!!!!! RHO et ALPH doivent etre a la fin! voir defmat!!!
  224. DO ie1=1,JGA
  225. MOTS(JGM+1-ie1)=MOTS(JGA+1-ie1)
  226. ENDDO
  227. C
  228.  
  229. MOTS(1)='TANB'
  230. MOTS(2)='BPLC'
  231. MOTS(3)='TR0T'
  232. MOTS(4)='TR0C'
  233. MOTS(5)='TRUC'
  234. MOTS(6)='AVAT'
  235. MOTS(7)='BVAT'
  236. MOTS(8)='AVAC'
  237. MOTS(9)='BVAC'
  238. MOTS(10)='DVOL'
  239. *
  240. * CCd
  241. *
  242. ELSE IF (IPLAC.EQ.11) THEN
  243. *
  244. * MODELE SICSCAL
  245. *
  246. MLMOTS=MOFAC
  247. JGA=MOTS(/2)
  248. JGM=JGA+44
  249. SEGADJ MLMOTS
  250. MOFAC=MLMOTS
  251. MOTS(JGA+1)= 'H1N'
  252. MOTS(JGA+2)= 'H1HP'
  253. MOTS(JGA+3)= 'H1P'
  254. MOTS(JGA+4)= 'H2N'
  255. MOTS(JGA+5)= 'H2HP'
  256. MOTS(JGA+6)= 'H2P'
  257. MOTS(JGA+7)= 'H3N'
  258. MOTS(JGA+8)= 'H3P'
  259. MOTS(JGA+9)= 'Y0N1'
  260. MOTS(JGA+10)='Y0N2'
  261. MOTS(JGA+11)='Y0N3'
  262. MOTS(JGA+12)= 'YCN1'
  263. MOTS(JGA+13)= 'YCN2'
  264. MOTS(JGA+14)= 'YCN3'
  265. MOTS(JGA+15)='Y0T1'
  266. MOTS(JGA+16)='Y0T2'
  267. MOTS(JGA+17)='Y0T3'
  268. MOTS(JGA+18)='YCT1'
  269. MOTS(JGA+19)='YCT2'
  270. MOTS(JGA+20)='YCT3'
  271. MOTS(JGA+21)='DCT1'
  272. MOTS(JGA+22)='DCT2'
  273. MOTS(JGA+23)='DCT3'
  274. MOTS(JGA+24)='DCN1'
  275. MOTS(JGA+25)='DCN2'
  276. MOTS(JGA+26)='DCN3'
  277. MOTS(JGA+27)='PN1'
  278. MOTS(JGA+28)='PN2'
  279. MOTS(JGA+29)='PN3'
  280. MOTS(JGA+30)='PT1'
  281. MOTS(JGA+31)='PT2'
  282. MOTS(JGA+32)='PT3'
  283. MOTS(JGA+33)='B'
  284. MOTS(JGA+34)='DTAL'
  285. MOTS(JGA+35)='TER0'
  286. MOTS(JGA+36)='SIF1'
  287. MOTS(JGA+37)='SIF2'
  288. MOTS(JGA+38)='SIF3'
  289. MOTS(JGA+39)='AIF1'
  290. MOTS(JGA+40)='AIF2'
  291. MOTS(JGA+41)='AIF3'
  292. MOTS(JGA+42)='ETA1'
  293. MOTS(JGA+43)='ETA2'
  294. MOTS(JGA+44)='ETA3'
  295. *
  296. ELSE IF (IPLAC.EQ.12) THEN
  297. *
  298. * MODELE SICTENS
  299. *
  300. MLMOTS=MOFAC
  301. JGA=MOTS(/2)
  302. JGM=JGA+51
  303. SEGADJ MLMOTS
  304. MOFAC=MLMOTS
  305. MOTS(JGA+1)= 'H1N'
  306. MOTS(JGA+2)= 'H1HP'
  307. MOTS(JGA+3)= 'H1P'
  308. MOTS(JGA+4)= 'H2N'
  309. MOTS(JGA+5)= 'H2HP'
  310. MOTS(JGA+6)= 'H2P'
  311. MOTS(JGA+7)= 'H3N'
  312. MOTS(JGA+8)= 'H3P'
  313. MOTS(JGA+9)= 'H4N'
  314. MOTS(JGA+10)= 'H4HP'
  315. MOTS(JGA+11)= 'H4P'
  316. MOTS(JGA+12)= 'Y01'
  317. MOTS(JGA+13)='Y02'
  318. MOTS(JGA+14)='Y03'
  319. MOTS(JGA+15)='Y04'
  320. MOTS(JGA+16)='Y05'
  321. MOTS(JGA+17)= 'YC1'
  322. MOTS(JGA+18)= 'YC2'
  323. MOTS(JGA+19)= 'YC3'
  324. MOTS(JGA+20)= 'YC4'
  325. MOTS(JGA+21)= 'YC5'
  326. MOTS(JGA+22)='DC1'
  327. MOTS(JGA+23)='DC2'
  328. MOTS(JGA+24)='DC3'
  329. MOTS(JGA+25)='DC4'
  330. MOTS(JGA+26)='DC5'
  331. MOTS(JGA+27)='PY1'
  332. MOTS(JGA+28)='PY2'
  333. MOTS(JGA+29)='PY3'
  334. MOTS(JGA+30)='PY4'
  335. MOTS(JGA+31)='PY5'
  336. MOTS(JGA+32)='B1'
  337. MOTS(JGA+33)='B2'
  338. MOTS(JGA+34)='B3'
  339. MOTS(JGA+35)='DTAL'
  340. MOTS(JGA+36)='TER0'
  341. MOTS(JGA+37)='SIF1'
  342. MOTS(JGA+38)='SIF2'
  343. MOTS(JGA+39)='SIF3'
  344. MOTS(JGA+40)='SIF4'
  345. MOTS(JGA+41)='SIF5'
  346. MOTS(JGA+42)='AIF1'
  347. MOTS(JGA+43)='AIF2'
  348. MOTS(JGA+44)='AIF3'
  349. MOTS(JGA+45)='AIF4'
  350. MOTS(JGA+46)='AIF5'
  351. MOTS(JGA+47)='ETA1'
  352. MOTS(JGA+48)='ETA2'
  353. MOTS(JGA+49)='ETA3'
  354. MOTS(JGA+50)='ETA4'
  355. MOTS(JGA+51)='ETA5'
  356. *
  357. ELSE IF (IPLAC.EQ.13) THEN
  358. *
  359. * MODELE DESMORAT
  360. *
  361. MLMOTS=MOOBL
  362. JGA=MOTS(/2)
  363. JGM=JGA+7
  364. SEGADJ MLMOTS
  365. MOOBL=MLMOTS
  366. MOTS(JGA+1)='K0'
  367. MOTS(JGA+2)='A'
  368. MOTS(JGA+3)='a'
  369. MOTS(JGA+4)='etaC'
  370. MOTS(JGA+5)='etaT'
  371. MOTS(JGA+6)='k'
  372. MOTS(JGA+7)='Dc'
  373. *
  374. ELSE IF (IPLAC.EQ.14) THEN
  375. *
  376. * MODELE ENDOMMAGEMENT FATIGUE 'FATSIN'
  377. *
  378. MLMOTS=MOOBL
  379. JGA=MOTS(/2)
  380. JGM=JGA+8
  381. SEGADJ MLMOTS
  382. MOOBL=MLMOTS
  383. MOTS(JGA+1)= 'KTR0'
  384. MOTS(JGA+2)= 'BETA'
  385. MOTS(JGA+3)= 'LOI '
  386. * pour la loi "L2R"
  387. MOTS(JGA+4)= 'C '
  388. MOTS(JGA+5)= 'ALFA'
  389. * pour la loi "L3R"
  390. MOTS(JGA+6)= 'ALF1'
  391. MOTS(JGA+7)= 'ALF2'
  392. MOTS(JGA+8)= 'ALF3'
  393. ELSE IF (IPLAC.EQ.15) THEN
  394. *
  395. * MODELE ENDOMMAGEMENT 'RICRAG'
  396. *
  397. MLMOTS=MOOBL
  398. JGA=MOTS(/2)
  399. JGM=JGA+5
  400. SEGADJ MLMOTS
  401. MOOBL=MLMOTS
  402. MOTS(JGA+1)='FT '
  403. MOTS(JGA+2)='ALIN'
  404. MOTS(JGA+3)='GAM1'
  405. MOTS(JGA+4)='A1 '
  406. MOTS(JGA+5)='ALDI'
  407.  
  408. C +BR
  409. ELSE IF (IPLAC.EQ.16) THEN
  410. *
  411. * MODELE ENDOMMAGEMENT 'GLRC_DM'
  412. *
  413. MLMOTS=MOOBL
  414. JGA=MOTS(/2)
  415. JGM=JGA+7
  416. SEGADJMLMOTS
  417. MOOBL=MLMOTS
  418. MOTS(JGA+1)='YOUF'
  419. MOTS(JGA+2)='NUF '
  420. MOTS(JGA+3)='GAMT'
  421. MOTS(JGA+4)='GAMC'
  422. MOTS(JGA+5)='GAMF'
  423. MOTS(JGA+6)='SEUI'
  424. MOTS(JGA+7)='ALF '
  425.  
  426. ELSE IF (IPLAC.EQ.17) THEN
  427. *
  428. * MODELE ENDOMMAGEMENT 'RICBET'
  429. *
  430. MLMOTS=MOOBL
  431. JGA=MOTS(/2)
  432. JGM=JGA+13
  433. SEGADJ MLMOTS
  434. MOOBL=MLMOTS
  435. MOTS(JGA+ 1)= 'FT '
  436. MOTS(JGA+ 2)= 'ALDI'
  437. MOTS(JGA+ 3)= 'GAM1'
  438. MOTS(JGA+ 4)= 'A1 '
  439. MOTS(JGA+ 5)= 'SREF'
  440. MOTS(JGA+ 6)= 'AF '
  441. MOTS(JGA+ 7)= 'BF '
  442. MOTS(JGA+ 8)= 'AG '
  443. MOTS(JGA+ 9)= 'BG '
  444. MOTS(JGA+10)= 'AC '
  445. MOTS(JGA+11)= 'BC '
  446. MOTS(JGA+12)= 'SIGU'
  447. MOTS(JGA+13)= 'FC '
  448.  
  449. ELSE IF (IPLAC.EQ.18) THEN
  450. *
  451. * MODELE ENDOMMAGEMENT 'RICCOQ'
  452. *
  453. MLMOTS=MOOBL
  454. JGA=MOTS(/2)
  455. JGM=JGA+4
  456. SEGADJ MLMOTS
  457. MOOBL=MLMOTS
  458. MOTS(JGA+1)='FT '
  459. MOTS(JGA+2)='EPUT'
  460. MOTS(JGA+3)='FC '
  461. MOTS(JGA+4)='EPUC'
  462.  
  463. ELSE IF (IPLAC.EQ.19) THEN
  464. *
  465. * MODELE ENDOMMAGEMENT 'CONCYC'
  466. *
  467. MLMOTS=MOOBL
  468. JGA=MOTS(/2)
  469. JGM=JGA+6
  470. SEGADJMLMOTS
  471. MOOBL=MLMOTS
  472. MOTS(JGA+1)='NEND'
  473. MOTS(JGA+2)='SIGT'
  474. MOTS(JGA+3)='ATRA'
  475. MOTS(JGA+4)='BTRA'
  476. MOTS(JGA+5)='QP '
  477. MOTS(JGA+6)='CF '
  478.  
  479. C -BR
  480. ELSE IF (IPLAC.EQ.20) THEN
  481. *
  482. * MODELE ENDOMMAGEMENT 'PBEFEM'
  483. *
  484. MLMOTS=MOOBL
  485. JGA=MOTS(/2)
  486. JGM=JGA+14
  487. SEGADJ MLMOTS
  488. MOOBL=MLMOTS
  489. MOTS(JGA+1) ='ALP1'
  490. MOTS(JGA+2) ='BET1'
  491. MOTS(JGA+3) ='AD1 '
  492. MOTS(JGA+4) ='AD2 '
  493. MOTS(JGA+5) ='AD3 '
  494. MOTS(JGA+6) ='AD4 '
  495. MOTS(JGA+7) ='FT '
  496. MOTS(JGA+8) ='XZ0 '
  497. MOTS(JGA+9) ='XNX '
  498. MOTS(JGA+10) ='XNY '
  499. MOTS(JGA+11) ='UCRT'
  500. MOTS(JGA+12) ='IND1'
  501. MOTS(JGA+13) ='XIE'
  502. MOTS(JGA+14) ='YIE'
  503.  
  504. ELSE
  505. *
  506. * SI PROBLEME , ON SORT AVEC IRET A 0
  507. *
  508. IRET=0
  509. CALL ERREUR(261)
  510. ENDIF
  511. RETURN
  512. END
  513.  
  514.  
  515.  
  516.  
  517.  
  518.  
  519.  
  520.  
  521.  
  522.  
  523.  
  524.  
  525.  
  526.  
  527.  
  528.  

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