Télécharger th_boucle.dgibi

Retour à la liste

Numérotation des lignes :

  1. * fichier : th_boucle.dgibi
  2. DEBP REEV_THE TB*TABLE IND*ENTIER ;
  3. mess 'REEV_THE' ;
  4. FINP ;
  5.  
  6. *--------------------------------------------------------*
  7.  
  8. DEBP REEV_MEC TB*TABLE IND*ENTIER ;
  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';
  16. TTTTT= estim.'TEMPS';
  17. wtab=tb.'WTABLE';
  18. CONF0 = FORME ;
  19. TB_REV = TB.'TABREP' ;
  20.  
  21. S1 = TB_REV . 'S1' ;
  22. S2 = TB_REV . 'S2' ;
  23. * FORME TB.'DEPI' ;
  24. FORME estim . 'DEPLACEMENTS';
  25. VAL_H = (0.02*TTTTT) + 0.1 ;
  26. SI (TTTTT > 10.) ;
  27. VAL_H = 0.3 ;
  28. FINSI ;
  29. mess 'TEMP0 -->' VAL_H ;
  30.  
  31. wtab . 'MAT_TOT' = TB_REV.'MATT1' ET TB_REV.'MATT2' ET
  32. TB_REV.'MATM1' ET TB_REV.'MATM2' ;
  33. wtab . 'MAT_TOT' = wtab . 'MAT_TOT' ET
  34. (mater (TB_REV.'MOT3') 'K' VAL_H 'C' 0. 'RHO' 0.) ;
  35. wtab.'MAT_COND'='REDU' wtab.'MAT_TOT' WTAB.'MOD_COND';
  36. FORME CONF0 ;
  37. FINP ;
  38.  
  39. *--------------------------------------------------------*
  40. *--------------------------------------------------------*
  41.  
  42.  
  43. opti dime 2 elem qua8 ;
  44. dens 1 ;
  45. * ajout de option epsilon lineaire pour la precision des test!
  46. OPTION epsilon lineaire;
  47. p0 = 0. 0. ; p1 = 1. 0. ;
  48. d1 = d 20 p0 p1;
  49. s1 = d1 tran 4 (0. 1.);
  50. clt1 = cote 4 s1 ;
  51. CLM1 = cote 2 s1 ;
  52.  
  53. p2 = 1.3 0. ; p3 = 1.6 0. ;
  54. d2 = d 7 p2 p3;
  55. s2 = d2 tran 4 (0. 1.);
  56. clt2 = cote 2 s2 ;
  57. CLM2 = COTE 4 s2 ;
  58. *trace (s1 et s2);
  59. coterac1 = cote 2 s1 ;
  60. coterac2 = cote 4 s2 ;
  61.  
  62. L1 = coterac1 ;
  63. L2 = INVE coterac2 ;
  64. II = 0 ;
  65. REPETER BOUC01 (NBELEM L1 ) ;
  66. II = II + 1 ;
  67. rac_int = raccord 1.0 (L1 elem II) (L2 ELEM II);
  68. SI ('EGA' II 1) ;
  69. RACTOT = RAC_INT ;
  70. SINON ;
  71. RACTOT = RACTOT ET RAC_INT ;
  72. FINSI ;
  73. FIN BOUC01 ;
  74.  
  75. ractot = coterac1 regler 7 (inve coterac2) ;
  76. mom1 = modele s1 mecanique elastique cons meca ;
  77. mom2 = modele s2 mecanique elastique cons meca ;
  78. matm1 = mater mom1 'YOUN' 1E10 'NU' 0. 'ALPHA' 1E-4 ;
  79. matm2 = mater mom2 'YOUN' 1E10 'NU' 0. 'ALPHA' 0. ;
  80. bl1 = (BLOQUE 'DEPL' P0) ET (BLOQUE 'DEPL' P3) ;
  81. bl2 = (BLOQUE 'UY' D1) ET (BLOQUE 'UY' D2) ;
  82. BL3 = 'BLOQ' UX CLM1 ;
  83. BL4 = 'BLOQ' UX CLM2 ;
  84.  
  85.  
  86. mot1 = modele s1 thermique isotrope cons ther ;
  87. mot2 = modele s2 thermique isotrope cons ther ;
  88. mot3 = modele RACTOT thermique isotrope cons ther ;
  89. matt1 = mater mot1 'K' .1 'RHO' 0. 'C' 0. ;
  90. matt2 = mater mot2 'K' .1 'RHO' 0. 'C' 0. ;
  91. matt3 = mater mot3 'K' .1 'RHO' 0. 'C' 0. ;
  92.  
  93. CL_THER1 = BLOQUE CLT2 'T' ;
  94. DI_PTH1 = DEPI CL_THER1 1500.;
  95. CL_THER2 = BLOQUE CLT1 'T' ;
  96. DI_PTH2 = DEPI CL_THER2 300. ;
  97. CHARG1 = CHARGE 'TIMP' DI_PTH1 (EVOL MANU (PROG 0. 10. 100.)
  98. (PROG 0.2 1. 1.)) ;
  99. CHARG2 = CHARGE 'TIMP' DI_PTH2 (EVOL MANU (PROG 0. 10. 100.)
  100. (PROG 1. 1. 1. )) ;
  101.  
  102. TAB1 = TABLE ;
  103. TAB1.'CHARGEMENT' = CHARG1 ET CHARG2 ;
  104. TAB1.'TEMPS_CALCULES' = PROG 0. PAS 4. 12. ;
  105. TAB1.'MODELE' = mot1 et mot2 et mot3 et mom1 et mom2 ;
  106. TAB1.'CARACTERISTIQUES' = MATT1 ET MATT2 ET MATT3
  107. ET MATM1 ET MATM2 ;
  108. TAB1.'BLOCAGES_MECANIQUES' = BL1 ET BL2 ET BL3 ET BL4 ;
  109. TAB1.'BLOCAGES_THERMIQUES' = CL_THER1 ET CL_THER2 ;
  110. TAB1.'TEMPERATURES' = TABLE ;
  111. TAB1.'TEMPERATURES' . 0 = MANU CHPO (S1 ET S2) 1 'T' 300. ;
  112. TAB1.'PROCEDURE_REAC_THER' = VRAI ;
  113. TAB1.'PROCEDURE_REAC_MECA' = VRAI ;
  114. TAB1.'CONVERGENCE_MEC_THE' = VRAI ;
  115. TAB1.'RELAXATION_THETA' = 1. ;
  116. TAB1.'NB_BOTH' = 5 ;
  117. TAB1.'CRITERE_COHERENCE' = 1.E-8 ;
  118. TAB1.'CONVERGENCE_FORCEE' = FAUX ;
  119. * TAB1.'MEC_THER' = FAUX ;
  120.  
  121. TABREP = TABLE ;
  122. TABREP . 'S1' = S1 ;
  123. TABREP . 'S2' = S2 ;
  124.  
  125. TABREP . 'MATT1' = MATT1 ; TABREP . 'MATT2' = MATT2 ;
  126. TABREP . 'MATM1' = MATM1 ; TABREP . 'MATM2' = MATM2 ;
  127. TABREP . 'MOT3' = MOT3 ;
  128.  
  129. TAB1 . 'TABREP' = TABREP ;
  130. *opti impi 1754;
  131. TAB1.'PROCESSEURS' = 'AUTOMATIQUE';
  132. tab1 =pasapas tab1;
  133.  
  134. VAL1 = TABLE ;
  135. VAL2 = TABLE ;
  136. II = -1 ;
  137. REPETER BOUCT (dime TAB1.temperatures);
  138. II = II + 1 ;
  139. CHTER = tab1.temperatures . II ;
  140. VAL1.II = (maxi (redu chter S1)) ;
  141. VAL2.II = (mini (redu chter S2)) ;
  142. mess 'T maxi S1: ' (maxi (redu chter S1))
  143. 'T mini S2: ' (mini (redu chter S2)) ;
  144. FIN BOUCT ;
  145.  
  146.  
  147. NIVERR = FAUX ;
  148. SI ((VAL1.0 'NEG' 300. 1.E-2) OU (VAL2.0 'NEG' 300. 1.E-2)) ;
  149. NIVERR = VRAI ;mess ' erreur 1' VAL1.0 VAL2.0 ;
  150. FINSI ;
  151. SI ((VAL1.1 'NEG' 627.27 1.E-2) OU (VAL2.1 'NEG' 681.82 1.E-2)) ;
  152. NIVERR = VRAI ;mess ' erreur 2' VAL1.1 VAL2.1 ;
  153. FINSI ;
  154. SI ((VAL1.2 'NEG' 978.26 1.E-2) OU (VAL2.2 'NEG' 1056.52 1.E-2)) ;
  155. NIVERR = VRAI ;mess ' erreur 3' VAL1.2 VAL2.2 ;
  156. FINSI ;
  157. SI ((VAL1.3 'NEG' 1157.14 1.E-2) OU (VAL2.3 'NEG' 1242.86 1.E-2)) ;
  158. NIVERR = VRAI ;mess ' erreur 4' VAL1.3 VAL2.3 ;
  159. FINSI ;
  160.  
  161. * je calcule les déplacements à la main sur mon exemple
  162.  
  163. tem1 = redu s1 (tab1.'TEMPERATURES' . 3) ;
  164. tem1 = tem1 - 300. ;
  165.  
  166. rig1 = (rigi matm1 mom1) et ('BLOQ' 'DEPL' P0)
  167. et BL3 et (BLOQUE UY D1);
  168. sisi = theta tem1 mom1 matm1 ;
  169. fequ = bsigma mom1 sisi ;
  170. chf1 = fequ ;
  171.  
  172. dep2 = ENLE (reso rig1 chf1) 'LX' ;
  173. dep1 = redu s1 (tab1 . deplacements . 3);
  174.  
  175. ddep = dep2 - dep1 ;
  176.  
  177. 'SI' ((MAXI (ABS DDEP)) > 1.E-8) ;
  178. NIVERR = VRAI ; mess ' erreur 5' ddep;
  179. 'FINSI' ;
  180.  
  181. SI NIVERR ;
  182. ERREUR 5 ;
  183. FINSI ;
  184.  
  185. FIN ;
  186.  
  187.  
  188.  
  189.  
  190.  
  191.  
  192.  
  193.  
  194.  
  195.  
  196.  
  197.  
  198.  
  199.  
  200.  

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