Télécharger la plaquette Cast3M


fr en

Télécharger ohno1.dgibi

Retour à la liste

Numérotation des lignes : 

  1. * fichier :  ohno1.dgibi
  2. ************************************************************************
  3. ************************************************************************
  4. * 
  5. complet = faux;
  6.  
  7. *
  8. *************************************************
  9. *  pour calcul complet mettre complet à : vrai;
  10. ***************************************************
  11. **************************************************
  12. *                                                *  
  13. *  TEST DE VALIDATION D'UNE LOI DE COMPORTEMENT  *  
  14. *  --------------------------------------------  * 
  15. *              DE MATERIAU                       * 
  16. *              -----------                       *   
  17. *                                                *  
  18. *  LOI DE COMPORTEMENT VISCOPLASTIQUE DE:        * 
  19. *   OHNO ( CHABOCHE MODIFIE OHNO )               * 
  20. *                                                *   
  21. *  TEST POUR DES ELEMENTS MASSIFS                *  
  22. *  ------------------------------                *  
  23. *                                                *  
  24. *  MAILLAGE:                                     *  
  25. *   UNE BARRE DE SECTION CARREE                  * 
  26. *   LONGUEUR L=.5 M                              *  
  27. *   LARGEUR  l=.01 M                             * 
  28. *                                                *  
  29. *  CHARGEMENT:                                   * 
  30. *   ESSAI DE TRACTION                            * 
  31. *   DEPLACEMENTS IMPOSES                         * 
  32. *                                                * 
  33. **************************************************  
  34. *
  35. *
  36. 'DEBPROC' VISOHNO ta1*'TABLE' i0*'ENTIER' 
  37.                   sigx01*'FLOTTANT' eps0*'FLOTTANT'  
  38.                   R01*'FLOTTANT' 
  39.                   epsx01*'FLOTTANT' 
  40.                   Xx11*'FLOTTANT' 
  41.                   Xx21*'FLOTTANT' ;   
  42. *  
  43. *  Résolution des lois viscoplastiques de OHNO dans 
  44. *  le cas unidirectionnel  
  45. *  
  46. mat0 = ta1 .'CARACTERISTIQUES' ;
  47. mod0 = ta1 .'MODELE' ;
  48. *q00 = 'MAXI' ( ta1 .'QQ' ) ;
  49. q00 = 0. ; 
  50. *
  51. youn0 = 'MAXI' ( 'EXCO' mat0 'YOUN' ) ; 
  52. n0 = 'MAXI' ( 'EXCO' mat0 'N' ) ;
  53. kk0 = 'MAXI' ( 'EXCO' mat0 'KK' ) ;
  54. k00 = 'MAXI' ( 'EXCO' mat0 'K0' ) ;
  55. alfk0 = 'MAXI' ( 'EXCO' mat0 'ALFK0' ) ;
  56. alfr0 = 'MAXI' ( 'EXCO' mat0 'ALFR0' ) ;
  57. alf0 = 'MAXI' ( 'EXCO' mat0 'ALF' ) ;
  58. a10 = 'MAXI' ( 'EXCO' mat0 'A1' ) ;
  59. c10 = 'MAXI' ( 'EXCO' mat0 'C1' ) ;
  60. bet10 = 'MAXI' ( 'EXCO' mat0 'BET1' ) ;
  61. r10 = 'MAXI' ( 'EXCO' mat0 'R1' ) ;
  62. a20 = 'MAXI' ( 'EXCO' mat0 'A2' ) ;
  63. c20 = 'MAXI' ( 'EXCO' mat0 'C2' ) ;
  64. bet20 = 'MAXI' ( 'EXCO' mat0 'BET2' ) ;
  65. r20 = 'MAXI' ( 'EXCO' mat0 'R2' ) ; 
  66. phi0 = 'MAXI' ( 'EXCO' mat0 'PHI' ) ;
  67. b0 = 'MAXI' ( 'EXCO' mat0 'B' ) ;
  68. gama0 = 'MAXI' ('EXCO' mat0 'GAMA' ) ;
  69. m0 = 'MAXI' ( 'EXCO' mat0 'M' ) ;
  70. qmax0 = 'MAXI ' ( 'EXCO' mat0 'QMAX' ) ;
  71. qsta0 = 'MAXI' ( 'EXCO' mat0 'QSTA' ) ;
  72. mu0 = 'MAXI' ( 'EXCO' mat0 'MU' ) ;
  73. eta0 = 'MAXI' ( 'EXCO' mat0 'ETA' ) ;
  74. exp10 = 'MAXI' ( 'EXCO' mat0 'EXP1' ) ;
  75. exp20 = 'MAXI' ( 'EXCO' mat0 'EXP2' ) ;
  76. * 
  77. depl0 = ta1 .'DEPLACEMENTS'.( i0 - 1 ) ; 
  78. temp1 = ta1 .'TEMPS'.i0 ; 
  79. temp0 = ta1 .'TEMPS'.( i0 - 1 ) ; 
  80. temp01 = temp0 ; 
  81. eps_old = 'MAXI' ( 'EXCO' ( 'EPSI' mod0 depl0 'LINE' ) 'EPXX' ) ;  
  82. *  
  83. niter0 = 50 ; 
  84. dt0 = ( temp1 - temp0 ) / niter0 ; 
  85. sigx00 = sigx01 ; 
  86. R00 = R01 ; 
  87. epsx00 = epsx01 ; 
  88. Xx10 = Xx11 ; 
  89. Xx20 = Xx21 ;  
  90. *  
  91. *  Incrément de contrainte test 
  92. *  
  93. spt0 = youn0 * ( eps0 - eps_old ) / ( temp1 - temp0 ) ;  
  94. *mess ' sigx00 R00 epsx00 Xx10 Xx20=' sigx00 R00 epsx00 Xx10 Xx20 ; 
  95. *  
  96. 'REPETER' bloc0 niter0 ; 
  97. *  
  98. *  Initialisation du calcul pour un pas de temps   
  99. *  
  100. temp01 = temp01 + dt0 ; 
  101. depsx9 = 0. ; 
  102. dXx8 = 0. ; 
  103. dXx9 = 0. ; 
  104. dr9 = 0. ;  
  105. sigx0 = sigx00 ; 
  106. R0 = R00 ; 
  107. epsx0 = epsx00 ; 
  108. Xx1 = Xx10 ; 
  109. Xx2 = Xx20 ; 
  110. *  
  111. 'REPETER' bloc1 2 ; 
  112.   i1 = &bloc1 ; 
  113. * 
  114. *  Memoire en déformation inélastique  
  115. *  
  116.   qq = eta0 * epsx0 ; 
  117.   Q0 = qmax0 - 
  118. ( ( qmax0 - q00 ) * ( 'EXP' ( -2. * mu0 * eta0 * epsx0 ) ) ) ;  
  119.   YX0 = ( 1. - eta0 ) * epsx0 ; 
  120.   QR = Q0 - ( qsta0 * ( 1. - ( ( ( qmax0 - Q0 ) / qmax0 ) ** 2. ) ) ) ;  
  121.   RR0 = alfr0 * R0 ; 
  122.   KR0 = k00 + ( alfk0 * R0 ) ;    
  123. *  mess ' qq, Q0, YX0, QR, RR0, KR0=' qq Q0 YX0 QR RR0 KR0 ; 
  124. *  
  125. *  Critère  
  126. *   
  127.   Xx0 = Xx1 + Xx2 ; 
  128.   F0 = sigx0 - ( 3. * Xx0 / 2. ) - RR0 - kk0 ;    
  129. * 
  130. *  Ecrouissages  
  131. *  
  132.   'SI' ( F0 '>' 0. ) ;  
  133.     depsx0 = ( F0 / KR0 ) ** ( n0 ) ; 
  134.     depsx1 = alf0 * ( ( F0 / KR0 ) ** ( n0 + 1. ) ) ; 
  135.     depsx1 = 'EXP' ( depsx1 ) ; 
  136.     depsx0 = depsx0 * depsx1 * dt0 ;  
  137.   'SINON' ; 
  138.     depsx0 = 0. ; 
  139.   'FINSI' ; 
  140. *  mess ' depsx0 F0 = ' depsx0 F0 ; 
  141.   dr0 = b0 * ( Q0 - R0 ) * depsx0 ; 
  142.   dr1 = gama0 * ( ( 'ABS' ( QR - R0 ) ) ** ( m0 - 1. ) ) ; 
  143.   dr1 = dr1 * ( QR - R0 ) ;  
  144.   dr0 = dr0 + ( dr1 * dt0 ) ;  
  145.   phi = phi0 + ( ( 1. - phi0 ) * ( 'EXP' ( -1. * b0 * epsx0 ) ) ) ; 
  146.   X1lim = a10 / phi ;
  147.   X2lim = a20 / phi ; 
  148.   dXx1 = 2. * a10 * c10 * depsx0 / 3. ; 
  149.   dXx1 = dXx1 - 
  150. ( c10 * Xx1 * depsx0 * phi * 
  151. (( 3. * Xx1 / ( 2. * X1lim )) ** ( exp10 ))) ; 
  152.   dXx1 = dXx1 - 
  153. ( ( ( 3. * Xx1 / ( 2. * bet10 ) ) ** ( r10 - 1. ) ) * Xx1 * dt0 ) ; 
  154.   dXx2 = 2. * a20 * c20 * depsx0 / 3. ;
  155.   dXx2 = dXx2 - 
  156. ( c20 * Xx2 * depsx0 * phi * 
  157. (( 3. * Xx2 / ( 2. * X2lim )) ** ( exp20 ))) ; 
  158.   dXx2 = dXx2 - 
  159. ( ( ( 3. * Xx2 / ( 2. * bet20 ) ) ** ( r20 - 1. ) ) * Xx2 * dt0 ) ; 
  160. *  mess ' dr0 phi dXx1 dXx2 i1 =' dr0 phi dXx1 dXx2 i1 ;  
  161. *  
  162. *  Somme ( Runge-Kutta ) 
  163. *  
  164.   dr0 = ( dr0 + dr9 ) / i1 ;  
  165.   depsx0 = ( depsx0 + depsx9 ) / i1 ; 
  166.   dXx1 = ( dXx1 + dXx8 ) / i1 ;  
  167.   dXx2 = ( dXx2 + dXx9 ) / i1 ; 
  168.   dr9 = dr0 ; 
  169.   depsx9 = depsx0 ; 
  170.   dXx8 = dXx1 ; 
  171.   dXx9 = dXx2 ; 
  172. *  mess ' dr0 depsx0 dXx1 dXx2 i1 =' dr0 depsx0 dXx1 dXx2 i1 ;
  173. *  
  174. *  Calcul des grandeurs actuelles 
  175. *  
  176.   R0 = R00 + dr0 ;   
  177.   epsx0 = epsx00 + depsx0 ; 
  178.   Xx1 = Xx10 + dXx1 ; 
  179.   Xx2 = Xx20 + dXx2 ;  
  180.   sigx0 = sigx00 + ( dt0 * spt0 ) 
  181.           - ( youn0 * depsx0 ) ;  
  182. *  mess ' R00 epsx00 Xx10 Xx20 sigx00 =' R00 epsx00 Xx10 Xx20 sigx00 ; 
  183. *  mess ' R0 epsx0 Xx1 Xx2 sigx0 =' R0 epsx0 Xx1 Xx2 sigx0 ;  
  184. *  
  185. 'FIN' bloc1 ; 
  186. *  
  187. *  Mise à jour des valeurs  
  188. *  
  189. R00 = R0 ; 
  190. epsx00 = epsx0 ; 
  191. Xx10 = Xx1 ; 
  192. Xx20 = Xx2 ; 
  193. sigx00 = sigx0 ; 
  194. *  
  195. 'FIN' bloc0 ; 
  196. * 
  197. 'MESS' ' Temps final ' temp01 temp1 ; 
  198. * 
  199. 'FINP' sigx00 R00 
  200.        epsx0 Xx10 Xx20 ;  
  201. *
  202. *--------------------------------------------------------------- 
  203. opti dime 3 elem cub8 ; 
  204. *  
  205. *  
  206. *  Maillage 
  207. *  
  208. l1 = .5 ; 
  209. l2 = .01 ; 
  210. l3 = .01 ;  
  211. n1 = 50 ; 
  212. n2 = 1 ; 
  213. n3 = 1 ; 
  214. p1 = 0. 0. 0. ; 
  215. p2 = l1 0. 0. ; 
  216. p3 = l1 l2 0. ; 
  217. p4 = 0. l2 0. ; 
  218. p5 = 0. 0. l3 ; 
  219. p6 = l1 0. l3 ; 
  220. p7 = l1 l2 l3 ; 
  221. p8 = 0. l2 l3 ; 
  222. d1 = d p1 n3 p5 ; 
  223. d2 = d p5 n2 p8 ; 
  224. d3 = d p8 n3 p4 ; 
  225. d4 = d p4 n2 p1 ; 
  226. sur1 = 'DALL' d1 d2 d3 d4 ; 
  227. d5 = d p2 n3 p6 ; 
  228. d6 = d p6 n2 p7 ;  
  229. d7 = d p7 n3 p3 ; 
  230. d8 = d p3 n2 p2 ; 
  231. sur2 = 'DALL' d5 d6 d7 d8 ; 
  232. vol1 = sur1 'VOLU' n1 sur2 ;    
  233. *  
  234. *    Modele de calcul 
  235. * 
  236. mod0 = 'MODE' vol1 'MECANIQUE' 'ELASTIQUE' 'ISOTROPE'   
  237.        'VISCOPLASTIQUE' 'OHNO' ; 
  238. mat0 = 'MATE' mod0 'YOUN' 7.34E10 'NU' 0.33   
  239.        'RHO' 7.8E3 'ALPHA' 1.E-5 'N' 24  
  240.        'KK' 10.E6 'K0' 116.E6 'ALFK' 1.5 'ALFR' 0.35 'ALF' 2.E6   
  241.        'A1' 67.5E6 'C1' 1300 'BET1' 4807.E6 'R1' 4 
  242.        'A2' 80.E6 'C2' 45 'BET2' 58480.E6 'R2' 4 'PHI' 1.   
  243.        'B' 12 'GAMA' 2.E-7 'M' 2 
  244.        'QMAX' 455.E6 'QSTA' 200.E6 'MU' 19 'ETA' 0.06    
  245.        'EXP1' 2. 'EXP2' 2. ; 
  246. *q00 = 'MANU' 'CHML' vol1 'QQ' 30.E6 ;  
  247. *  
  248. *    Conditions aux limites 
  249. *  
  250. cl1 = 'BLOQ' sur1 'UX' ; 
  251. cl2 = 'BLOQ' sur2 'UX' ; 
  252. cl3 = 'BLOQ' p1 'UX' 'UY' 'UZ' ;  
  253. cl4 = 'BLOQ' p5 'UX' 'UY' ; 
  254. cl0 = cl1 'ET' cl2 'ET' cl3 'ET' cl4 ;  
  255. *  
  256. *    Chargement 
  257. *  
  258. dep1 = 'DEPI' cl2 .02 ; 
  259. ev0 = 'EVOL' 'MANU' temps ( 'PROG' 0. 10000. ) y ( 'PROG' 0. 1. ) ; 
  260. cha0 = 'CHAR' 'DIMP' dep1 ev0 ; 
  261. *  
  262. *  Temps du calcul  
  263. *  
  264. dt0 = 50. ; 
  265. si complet;
  266. tfin0 = 1000. ;
  267. sinon;
  268. tfin0 = 300.;
  269. finsi;  
  270. *   
  271. *   
  272. *  Resolution par PASAPAS  
  273. *---------------------------------------
  274. *    
  275. ta1 = 'TABLE' ; 
  276. *ta1 .'QQ' = 'TABLE' ; 
  277. ta1 .'MODELE' = mod0 ; 
  278. ta1 .'CARACTERISTIQUES' = mat0 ; 
  279. ta1 .'BLOCAGES_MECANIQUES' = cl0 ; 
  280. ta1 .'CHARGEMENT' = cha0 ; 
  281. *ta1 .'QQ'. 0 = q00 ;  
  282. ta1 .'TEMPS_CALCULES' = 'PROG' 0. 'PAS' dt0 tfin0 ; 
  283. ta1 .'HYPOTHESE_DEFORMATIONS' = 'LINEAIRE' ;
  284. *  
  285. TMASAU=table;
  286. ta1 . 'MES_SAUVEGARDES'=TMASAU;
  287. TMASAU .'DEFTO'=VRAI;
  288. TMASAU .'DEFIN'=VRAI;
  289. PASAPAS ta1 ; 
  290. *  
  291. mat0 = ta1 .'CARACTERISTIQUES' ;  
  292. mod0 = ta1 .'MODELE' ; 
  293. *q00 = ta1 .'QQ' ; 
  294. q00 = 0. ; 
  295. * 
  296. n0    = 'MAXI' ( 'EXCO' mat0 'N'    ) ; 
  297. kk0   = 'MAXI' ( 'EXCO' mat0 'KK'   ) ; 
  298. k00   = 'MAXI' ( 'EXCO' mat0 'K0'   ) ; 
  299. alfk0 = 'MAXI' ( 'EXCO' mat0 'ALFK' ) ; 
  300. alfr0 = 'MAXI' ( 'EXCO' mat0 'ALFR' ) ; 
  301. alf0  = 'MAXI' ( 'EXCO' mat0 'ALF'  ) ; 
  302. a10   = 'MAXI' ( 'EXCO' mat0 'A1'   ) ; 
  303. c10   = 'MAXI' ( 'EXCO' mat0 'C1'   ) ; 
  304. bet10 = 'MAXI' ( 'EXCO' mat0 'BET1' ) ; 
  305. r10   = 'MAXI' ( 'EXCO' mat0 'R1'   ) ; 
  306. a20   = 'MAXI' ( 'EXCO' mat0 'A2'   ) ; 
  307. c20   = 'MAXI' ( 'EXCO' mat0 'C2'   ) ; 
  308. bet20 = 'MAXI' ( 'EXCO' mat0 'BET2' ) ; 
  309. r20   = 'MAXI' ( 'EXCO' mat0 'R2'   ) ; 
  310. phi0  = 'MAXI' ( 'EXCO' mat0 'PHI'  ) ; 
  311. b0    = 'MAXI' ( 'EXCO' mat0 'B'    ) ; 
  312. gama0 = 'MAXI' ( 'EXCO' mat0 'GAMA' ) ; 
  313. m0    = 'MAXI' ( 'EXCO' mat0 'M'    ) ; 
  314. qmax0 = 'MAXI' ( 'EXCO' mat0 'QMAX' ) ; 
  315. qsta0 = 'MAXI' ( 'EXCO' mat0 'QSTA' ) ; 
  316. mu0   = 'MAXI' ( 'EXCO' mat0 'MU'   ) ; 
  317. eta0  = 'MAXI' ( 'EXCO' mat0 'ETA'  ) ; 
  318. exp10 = 'MAXI' ( 'EXCO' mat0 'EXP1' ) ;
  319. exp20 = 'MAXI' ( 'EXCO' mat0 'EXP2' ) ; 
  320. *  
  321. dim0 = 'DIME' ta1 .'TEMPS' ; 
  322. *
  323. sig_th0 = 0. ; 
  324. r_th0 = 0. ; 
  325. eps_th0 = 0. ; 
  326. xx_th10 = 0. ; 
  327. xx_th20 = 0. ; 
  328. sig_th1 = 'PROG' 0. ; 
  329. eps_th1 = 'PROG' 0. ;
  330. si ( non complet);
  331. sig_th1 = 'PROG' 0. 
  332.  1.46800E+07  2.93600E+07  4.40400E+07  5.87196E+07  7.30534E+07
  333.   8.20293E+07  9.01555E+07  9.82258E+07  1.06163E+08
  334.   1.13855E+08  1.21148E+08  1.27849E+08  1.33737E+08  1.38604E+08 
  335.  1.42342E+08  1.45023E+08  1.46884E+08  1.48214E+08  1.49240E+08
  336.   1.50108E+08;
  337. eps_th1 = 'PROG' 0.
  338.  2.62535E-34  5.81971E-19  7.73859E-13  6.03725E-09  4.72149E-06 
  339.  8.24351E-05  1.71724E-04  2.61774E-04  3.53633E-04
  340.   4.48845E-04  5.49487E-04  6.58189E-04  7.77970E-04  9.11658E-04 
  341.  1.06073E-03  1.22421E-03  1.39885E-03  1.58074E-03  1.76676E-03
  342.   1.95494E-03 ;
  343. finsi;
  344. *  
  345. *  
  346. 'REPETER' bloc0 ( dim0 - 1 ) ; 
  347.   i0 = &bloc0 ; 
  348.   depl0 = ta1 .'DEPLACEMENTS'.i0 ; 
  349.   sigm0 = ta1 .'CONTRAINTES'.i0 ; 
  350.   var0 = ta1 .'VARIABLES_INTERNES'.i0 ; 
  351.   def0 = ta1 .'DEFORMATIONS_INELASTIQUES'.i0 ; 
  352.   temp0 = ta1 .'TEMPS'.i0 ; 
  353.   eps0 = 'MAXI' ( 'EXCO' ( 'EPSI' mod0 depl0 'LINE' ) 'EPXX' ) ;  
  354. *  
  355. *  Contraintes et déformations inélastiques équivalentes 
  356. *  
  357.   vm0 = 'MAXI' ( 'VMIS' sigm0 mod0 ) ; 
  358.   epsp0 = 'MAXI' ( 'EXCO' var0 'EPSE' ) ; 
  359. *  
  360. *  Déformation totale et inélastique suivant x 
  361. *  
  362.   epst0 = 'MAXI' ( 'EXCO' ( 'EPSI' depl0 mod0 'LINE' ) 'EPXX' ) ; 
  363.   epsx0 = 'MAXI' ( 'EXCO' def0 'EIXX' ) ; 
  364. *  
  365. *  Contraintes suivant x 
  366. *  
  367.   sigx0 = 'MAXI' ( 'EXCO' sigm0 'SMXX' ) ; 
  368. *  
  369. *  Calcul theorique  
  370. *  
  371. si complet;
  372. sig_th00 r_th00 eps_th00 xx_th1 xx_th2 = VISOHNO ta1 i0 sig_th0 
  373. eps0 r_th0 eps_th0 xx_th10 xx_th20 ;   
  374. *  
  375.   sig_th0 = sig_th00 ; 
  376.   r_th0 = r_th00 ; 
  377.   eps_th0 = eps_th00 ; 
  378.   xx_th10 = xx_th1 ; 
  379.   xx_th20 = xx_th2 ; 
  380.  sig_th1 = sig_th1 'ET' ( 'PROG' sig_th0 ) ; 
  381.  eps_th1 = eps_th1 'ET' ( 'PROG' eps_th0 ) ; 
  382. sinon;
  383.   sig_th00 = extr  sig_th1 ( i0 + 1);
  384.   eps_th00 = extr eps_th1 ( i0 + 1);
  385. finsi;
  386.  
  387. *  
  388. *  Calcul d'erreur  
  389. *  
  390.   'SI' ( epsx0 '<' ( 1.E-2 * eps0) ) ; 
  391.     err_e0 = 0. ; 
  392.   'SINON' ; 
  393.     err_e0 = ( 'ABS' ( eps_th00 - epsx0 ) ) / epsx0 ;  
  394.     err_e0 = err_e0 * 100. ; 
  395.   'FINSI' ; 
  396.   'SI' ( sigx0 'EGA' 0. ) ; 
  397.     err_s0 = 0. ; 
  398.   'SINON' ; 
  399.     err_s0 = ( 'ABS' ( sig_th00 - sigx0 ) ) / sigx0 ; 
  400.     err_s0 = err_s0 * 100. ; 
  401.   'FINSI' ; 
  402. *   
  403.   'SI' ( (err_e0 '>' 1.E-1) 'OU' (err_s0 '>' 1.E-1 ) ) ; 
  404.      prob0 = vrai ; 
  405.   'SINON' ; 
  406.      prob0 = faux ; 
  407.   'FINSI' ; 
  408.   'MESS' ' Fin du calcul de restitution au temps:' temp0 ; 
  409. *  
  410. 'FIN' bloc0 ;   
  411. *  
  412. 'SI' ( prob0 'EGA' vrai ) ; 
  413.    'ERRE' 5 ; 
  414. 'SINON' ; 
  415.    'ERRE' 0 ; 
  416. 'FINSI' ; 
  417. *
  418. 'FIN' ;   
  419.  
  420.  
  421.  
  422.  
  423.  
  424.  
  425.  
  426.  
  427.  
  428.  
  429.  
  430.  
  431.  
  432.  
  433.  
  434.  
  435.  

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