Télécharger th_non_boucle.dgibi

Retour à la liste

Numérotation des lignes :

  1. * fichier : th_non_boucle.dgibi
  2. DEBP REEV_THE TB*TABLE TBM1/TABLE ;
  3. mess 'REEV_THE' ;
  4. FINP ;
  5.  
  6. *--------------------------------------------------------*
  7.  
  8. DEBP REEV_MEC TB*TABLE TBM1/TABLE ;
  9. mess 'REEV_MEC' ;
  10. * SI (EXISTE TB 'TT') ;
  11. * TTTTT = TB. 'TT' . 'T_PREC' ;
  12. ** SINON ;
  13. * TTTTT = 0. ;
  14. * FINSI ;
  15. estim = tb.'ESTIMATION';wtab=tb.'WTABLE';
  16. TTTTT= estim.'TEMPS';
  17. CONF0 = FORME ;
  18. TB_REV = TB.'TABREP' ;
  19.  
  20. S1 = TB_REV . 'S1' ;
  21. S2 = TB_REV . 'S2' ;
  22. * FORME TB.'DEPI' ;
  23. FORM estim . 'DEPLACEMENTS';
  24. VAL_H = (0.02*TTTTT) + 0.1 ;
  25. SI (TTTTT > 10.) ;
  26. VAL_H = 0.3 ;
  27. FINSI ;
  28. mess TEMP0 '-->' VAL_H ;
  29.  
  30. wtab . 'MAT_TOT' = TB_REV.'MATT1' ET TB_REV.'MATT2' ET
  31. TB_REV.'MATM1' ET TB_REV.'MATM2' ;
  32. wtab . 'MAT_TOT' = wtab . 'MAT_TOT' ET
  33. (mater (TB_REV.'MOT3') 'K' VAL_H 'C' 0. 'RHO' 0.) ;
  34. wtab.'MAT_COND'='REDU' wtab.'MAT_TOT' wtab.'MOD_COND';
  35. FORME CONF0 ;
  36. FINP ;
  37.  
  38. *--------------------------------------------------------*
  39. *--------------------------------------------------------*
  40.  
  41.  
  42. opti dime 2 elem qua8 ;
  43. dens 1 ;
  44.  
  45. p0 = 0. 0. ; p1 = 1. 0. ;
  46. d1 = d 20 p0 p1;
  47. s1 = d1 tran 4 (0. 1.);
  48. clt1 = cote 4 s1 ;
  49. CLM1 = cote 2 s1 ;
  50.  
  51. p2 = 1.3 0. ; p3 = 1.6 0. ;
  52. d2 = d 7 p2 p3;
  53. s2 = d2 tran 4 (0. 1.);
  54. clt2 = cote 2 s2 ;
  55. CLM2 = COTE 4 s2 ;
  56. *trace (s1 et s2);
  57. coterac1 = cote 2 s1 ;
  58. coterac2 = cote 4 s2 ;
  59.  
  60. L1 = coterac1 ;
  61. L2 = INVE coterac2 ;
  62. II = 0 ;
  63. REPETER BOUC01 (NBELEM L1 ) ;
  64. II = II + 1 ;
  65. rac_int = raccord 1.0 (L1 elem II) (L2 ELEM II);
  66. SI ('EGA' II 1) ;
  67. RACTOT = RAC_INT ;
  68. SINON ;
  69. RACTOT = RACTOT ET RAC_INT ;
  70. FINSI ;
  71. FIN BOUC01 ;
  72.  
  73. ractot = coterac1 regler 7 (inve coterac2) ;
  74. mom1 = modele s1 mecanique elastique cons meca ;
  75. mom2 = modele s2 mecanique elastique cons meca ;
  76. matm1 = mater mom1 'YOUN' 1E10 'NU' 0. 'ALPHA' 1E-4 ;
  77. matm2 = mater mom2 'YOUN' 1E10 'NU' 0. 'ALPHA' 0. ;
  78. bl1 = (BLOQUE 'DEPL' P0) ET (BLOQUE 'DEPL' P3) ;
  79. bl2 = (BLOQUE 'UY' D1) ET (BLOQUE 'UY' D2) ;
  80. BL3 = BLOQ UX CLM1 ;
  81. BL4 = BLOQ UX CLM2 ;
  82.  
  83.  
  84. mot1 = modele s1 thermique isotrope cons ther ;
  85. mot2 = modele s2 thermique isotrope cons ther ;
  86. mot3 = modele RACTOT thermique isotrope cons ther ;
  87. matt1 = mater mot1 'K' .1 'RHO' 0. 'C' 0. ;
  88. matt2 = mater mot2 'K' .1 'RHO' 0. 'C' 0. ;
  89. matt3 = mater mot3 'K' .1 'RHO' 0. 'C' 0. ;
  90.  
  91. CL_THER1 = BLOQUE CLT2 'T' ;
  92. DI_PTH1 = DEPI CL_THER1 1500.;
  93. CL_THER2 = BLOQUE CLT1 'T' ;
  94. DI_PTH2 = DEPI CL_THER2 300. ;
  95. CHARG1 = CHARGE 'TIMP' DI_PTH1 (EVOL MANU (PROG 0. 10. 100.)
  96. (PROG 0.2 1. 1.)) ;
  97. CHARG2 = CHARGE 'TIMP' DI_PTH2 (EVOL MANU (PROG 0. 10. 100.)
  98. (PROG 1. 1. 1. )) ;
  99.  
  100. TAB1 = TABLE ;
  101. TAB1.'CHARGEMENT' = CHARG1 ET CHARG2 ;
  102. TAB1.'TEMPS_CALCULES' = PROG 0. PAS 4. 12. ;
  103. TAB1.'MODELE' = mot1 et mot2 et mot3 et mom1 et mom2 ;
  104. TAB1.'CARACTERISTIQUES' = MATT1 ET MATT2 ET MATT3
  105. ET MATM1 ET MATM2 ;
  106. TAB1.'BLOCAGES_MECANIQUES' = BL1 ET BL2 ET BL3 ET BL4 ;
  107. TAB1.'BLOCAGES_THERMIQUES' = CL_THER1 ET CL_THER2 ;
  108. TAB1.'TEMPERATURES' = TABLE ;
  109. TAB1.'TEMPERATURES' . 0 = MANU CHPO (S1 ET S2) 1 'T' 300. ;
  110. TAB1.'PROCEDURE_REEV_THE' = VRAI ;
  111. TAB1.'PROCEDURE_REEV_MEC' = VRAI ;
  112. TAB1.'CONVERGENCE_MEC_THE' = FAUX ;
  113. TAB1.'RELAXATION_THETA' = 1. ;
  114.  
  115. TABREP = TABLE ;
  116. TABREP . 'S1' = S1 ;
  117. TABREP . 'S2' = S2 ;
  118.  
  119. TABREP . 'MATT1' = MATT1 ; TABREP . 'MATT2' = MATT2 ;
  120. TABREP . 'MATM1' = MATM1 ; TABREP . 'MATM2' = MATM2 ;
  121. TABREP . 'MOT3' = MOT3 ;
  122.  
  123. TAB1 . 'TABREP' = TABREP ;
  124.  
  125. * TAB2 = COPIER TAB1 ;
  126. pasapas tab1;
  127.  
  128. VAL1 = TABLE ;
  129. VAL2 = TABLE ;
  130. II = -1 ;
  131. REPETER BOUCT (dime TAB1.temperatures);
  132. II = II + 1 ;
  133. CHTER = tab1.temperatures . II ;
  134. VAL1.II = (maxi (redu chter S1)) ;
  135. VAL2.II = (mini (redu chter S2)) ;
  136. mess 'T maxi S1: ' (maxi (redu chter S1))
  137. 'T mini S2: ' (mini (redu chter S2)) ;
  138. FIN BOUCT ;
  139.  
  140. NIVERR = FAUX ;
  141. SI ((VAL1.0 'NEG' 300. 1.E-2) OU (VAL2.0 'NEG' 300. 1.E-2)) ;
  142. NIVERR = VRAI ;
  143. FINSI ;
  144. SI ((VAL1.1 'NEG' 600.0 1.E-2) OU (VAL2.1 'NEG' 690.0 1.E-2)) ;
  145. NIVERR = VRAI ;
  146. FINSI ;
  147. SI ((VAL1.2 'NEG' 954.55 1.E-1) OU (VAL2.2 'NEG' 1063.64 1.E-1)) ;
  148. NIVERR = VRAI ;
  149. FINSI ;
  150. SI ((VAL1.3 'NEG' 1147.83 1.E-1) OU (VAL2.3 'NEG' 1245.65 1.E-1)) ;
  151. NIVERR = VRAI ;
  152. FINSI ;
  153.  
  154. SI NIVERR ;
  155. ERREUR 5 ;
  156. FINSI ;
  157.  
  158. FIN ;
  159.  
  160.  
  161.  
  162.  
  163.  
  164.  
  165.  
  166.  
  167.  

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