Télécharger la plaquette Cast3M


fr en

Télécharger dy_devo4.dgibi

Retour à la liste

Numérotation des lignes : 

  1. * fichier :  dy_devo4.dgibi
  2. ************************************************************************
  3. ************************************************************************
  4. *         Test Dy_devo4.dgibi: Jeux de donnees        *
  5. *         ------------------------------------        *
  6. *                                                     *
  7. * si GRAPH = N, les graphiques ne sont pas affiches                             
  8. * si GRAPH different de N,  tous les graphiques 
  9. * sont affiches                    
  10. *                                                                               
  11. GRAPH = 'N' ;  d_dyne = faux ;    
  12. * test comparaison pasapas et dyn                                                       
  13. *                                                                               
  14. SAUT PAGE ;                                                                     
  15. *                                                                               
  16. *-----------------------------------------------------*          
  17. *                    DY_DEVO4                         *          
  18. *-----------------------------------------------------*          
  19. *                                                     *          
  20. *    Tests de l'op{rateur DYNE option DE_VOGELAERE    *          
  21. *    _____________________________________________    *          
  22. *                                                     *          
  23. *                                                     *          
  24. * 1 poutre encastree, 10 elements finis,              * 
  25. * masse negligeable                                   *
  26. * 1 masse concentree e son extremite                  *          
  27. * Test PROFIL_PROFIL_EXTERIEUR avec des carres        *          
  28. * Comparaison avec la liaison POINT_PLAN              *          
  29. * NB_M: nombre de segments sur le cote du carre mobile*          
  30. * NB_F: nombre de segments sur le cote du carre fixe  *          
  31. *                                                     *          
  32. *-----------------------------------------------------*        
  33.  
  34. OPTION ECHO 1 ;
  35. TEMPS ;
  36. *
  37. NB_M = 3  ;
  38. NB_F = 1  ;
  39. R_PPL = 1.E+05 ;
  40. R_PP  = R_PPL  ;
  41. EXPO = 0.5 ;
  42. *
  43. OPTION ELEM SEG2 DIME 3 ;
  44. *
  45. *     Definition geometrie
  46. *
  47. P1 = 0. 0. 0. ;
  48. P2 = 0. 0. 0.9 ;
  49. LIG1 = P1 D 10 P2 ;
  50. *
  51. *     Definition materiau
  52. *
  53. MODL1 = MODE LIG1 MECANIQUE ELASTIQUE POUT;
  54.  
  55. MA1 = MATE MODL1 YOUNG 2.1E11 NU 0.3 RHO 78.  ;
  56. CA1 = 'CARA' MODL1 SECT 1.1E-6  INRY 8.33E-10  
  57.              INRZ 8.33E-10 TORS 5.0E-7  VECT (1. 0. 1. );
  58. *
  59. RIG1 = RIGI MODL1 (MA1 ET CA1) ;
  60. MAST = (MASSE MODL1  (MA1 ET CA1)) ET 
  61.        ( MASSE DEPL 2.0 P2 ) ;
  62. *
  63. *     Conditions aux limites
  64. *
  65. ENC1 = BLOQUE DEPL ROTA P1 ;
  66. ENC2 = BLOQUE 'UY' LIG1 ;
  67. RIGT = RIG1 ET ENC1 ET ENC2 ;
  68. *
  69. *     Calcul des modes
  70. *
  71. TBAS = VIBR PROCHE ( PROG 1.0  )  RIGT MAST ;
  72. *
  73. TB1 = TBAS.'MODES';
  74.  
  75.  
  76.     ma1 = manu poi1 tb1 . 1 . point_repere ;
  77.     mo1 = mode ma1 mecanique elastique modal ;
  78.     ca1 = mate mo1 freq tb1 . 1 . frequence 
  79.   mass tb1 . 1 . masse_generalisee defo tb1 . 1 . deformee_modale
  80.   amor 0. ;
  81.  
  82. *
  83. RIGI1 = RIGI TB1 ;
  84. *
  85. *     Amortissement
  86. *
  87. AMOR1 = AMOR TBAS ( PROG 0.0 ) ;
  88. amor2 = amor mo1 ca1 ;
  89. *
  90. TA = TABLE 'AMORTISSEMENT' ;
  91. TA.'AMORTISSEMENT' = AMOR1 ;
  92. *
  93. *     Deplacement initial
  94. *
  95. FE2 = FORCE FX 1. P2 ;
  96. XE2 = RESOU RIGT  FE2   ;
  97. XMAX = MAXI XE2 AVEC (MOTS UX ) ;
  98. FE1 = FE2 * (0.139E-2 / XMAX );
  99. FGE1 = PJBA TBAS FE1 ;
  100. ALO = RESO RIGI1 FGE1 ;
  101. *
  102. TINIT = TABLE 'INITIAL' ;
  103. TINIT.'DEPLACEMENT' = ALO ;
  104. *
  105. *---  Definition des liaisons
  106. *
  107. TLIA = TABLE 'LIAISON' ;
  108. TLB = TABLE 'LIAISON_B' ;
  109. TLIA.'LIAISON_B' = TLB ;
  110. NORM1 = 0 0 1 ;
  111. *
  112. *   Liaison PROFIL_PROFIL_EXTERIEUR
  113. *
  114. *   profil mobile
  115. P1_1 =   -0.5    0.0   0.9 ;
  116. P1_2 =    0.0   -0.5   0.9 ;
  117. P1_3 =    0.5    0.0   0.9 ;
  118. P1_4 =    0.0    0.5   0.9 ;
  119. MAIL_1 = P1_1 D NB_M P1_2 D NB_M P1_3 D NB_M P1_4 D 
  120.          NB_M P1_1 ;
  121. *
  122. *   profil fixe
  123. P2_1 =    0.501   1.0    0.9 ;
  124. P2_2 =    0.501  -1.0    0.9 ;
  125. P2_3 =    0.601  -1.0    0.9 ;
  126. P2_4 =    0.601   1.0    0.9 ;
  127. MAIL_2 = P2_1 D NB_F P2_2 D NB_F P2_3 D NB_F P2_4 D 
  128.          NB_F P2_1 ;
  129. *
  130. SI ( NEG GRAPH 'N' ) ;
  131.    TRAC ( 0 0 1000 ) ( MAIL_1 ET MAIL_2 ) ;
  132. FINSI ;
  133. *
  134. TL1  =  TABLE  'LIAISON_ELEMENTAIRE'    ;
  135. TL1.'TYPE_LIAISON'     = MOT 'PROFIL_PROFIL_EXTERIEUR';
  136. TL1.'SUPPORT'          =      P2     ;
  137. TL1.'PROFIL_FIXE'      =    MAIL_2   ;
  138. TL1.'PROFIL_MOBILE'    =    MAIL_1   ;
  139. TL1.'NORMALE'          =    NORM1    ;
  140. TL1.'RAIDEUR'          =    R_PP     ;
  141. TL1.'EXPOSANT_RAIDEUR' =    EXPO     ;
  142. *
  143. TLB . 1  = TL1  ;
  144. mappe1 = manu poi1 p2 ;
  145. moppe1 = modele mappe1 liaison PROFIL_PROFIL EXTERIEUR ;
  146.  tuppe1 = table ; tuppe1 . moppe1 = vrai ;
  147. cappe1 = mate moppe1 'NORM' norm1 'RAID' r_pp 'ERAI' expo
  148.      'PFIX' mail_2 'PMOB' mail_1 'SORT' tuppe1 ;
  149.  
  150. *
  151. *---  Definition des resultats en sortie
  152. *
  153. TSORT = TABLE 'SORTIE' ;
  154. TSOR1 = TABLE 'SORTIE' ;
  155. TSORT . 'VARIABLE'    = TSOR1 ;
  156. TSOR1 . 'DEPLACEMENT' =  VRAI ;
  157. TSOR1 . 'VITESSE'     =  VRAI ;
  158. *
  159. TSOR2 = TABLE 'SORTIE' ;
  160. TSORT.'LIAISON_B' = TSOR2 ;
  161. TSOR2.TL1  = VRAI ;
  162. *
  163. *     Temps
  164. *
  165. PDT  = 0.1E-3 ;
  166. NP1  =  5000  ;
  167. NINS =     5  ;
  168. *
  169. ev1 = evol manu (prog 0. 1000.) (prog 1. 1.) ;
  170. cha1 = char force ev1 FGE1 ;
  171. *
  172. SI d_dyne ;
  173. tabdyn = table 'PASAPAS' ;
  174. tabdyn . modele = mo1 et moppe1 ;
  175. tabdyn . caracteristiques = ca1 et cappe1 ;
  176. tabdyn . sortie = table ;
  177. tabdyn . sortie . moppe1 = vrai ;
  178. tabdyn . sortie . vitesses = faux ;
  179. tabdyn . pas_de_temps = pdt ;
  180. tabdyn . nombre_pas = np1 ;
  181. tabdyn . pas_de_sortie = nins ;
  182. tabdyn . temps = table ;
  183. tabdyn . temps . 0 = 0. ;
  184. tabdyn . deplacements = table ;
  185. tabdyn . deplacements . 0 =  ALO ;
  186.   dyne de_vogelaere tabdyn ;
  187. SINON ;
  188. * pour pasapas 
  189. * etape precontrainte
  190. t_prec = table ;
  191. t_prec . modele = mo1 ;
  192. t_prec . caracteristiques = ca1 ;
  193. t_prec . chargement = cha1 ;
  194. t_prec . temps_calcules  = prog 0. 1000. ;
  195. t_prec . mes_sauvegardes = table ;
  196. t_prec . mes_sauvegardes . defto = vrai ;
  197. t_prec . hypothese_deformations = 'LINEAIRE' ;
  198.   pasapas t_prec ;
  199.  
  200. *lacher
  201. tabdyn = table 'PASAPAS' ;
  202. tabdyn . modele = mo1 et moppe1 ;
  203. tabdyn . caracteristiques = ca1 et cappe1 ;
  204. *tabdyn . modele = mo1 ;
  205. *tabdyn . caracteristiques = ca1 ;
  206. tabdyn . temps = table ;
  207. tabdyn . temps . 0 = 0. ;
  208. tabdyn . deplacements = table ;
  209. tabdyn . deplacements . 0 =  t_prec . deplacements . 1 ;
  210. tabdyn . contraintes = table ;
  211. tabdyn . contraintes . 0 = t_prec . contraintes . 1 ;
  212. tabdyn . deformations = table;
  213. tabdyn . deformations . 0 = t_prec . deformations . 1;
  214. tabdyn . dynamique = vrai ;
  215. tabdyn . temps_calcules = prog 0 pas (4*nins*pdt) (np1*pdt);
  216. tabdyn . hypothese_deformations = 'LINEAIRE' ;
  217.  pasapas tabdyn ;
  218.  
  219. FINSI ;
  220. *------  Op{rateur DYNE  ------
  221. *
  222. TRESU1 = DYNE 'DE_VOGELAERE' TBAS TA TINIT TSORT TLIA 
  223.          NP1 PDT NINS ;
  224.  
  225.  nsu1 = dime tresu1 . temps_de_sortie ;
  226.  ovi1 = prog ;
  227.  repeter bsu1 nsu1 ;
  228.    ovi1 = ovi1 et 
  229.   (prog (extr tresu1 . &bsu1 . vitesse alfa tb1 . 1 . point_repere)) ;
  230.  fin bsu1 ;
  231.  evvi1 = evol manu 'temps' TRESU1 . 'TEMPS_DE_SORTIE'
  232.       'vitesse' ovi1 ;
  233. *
  234. *     Trac{ des courbes
  235. *
  236. SI ( NEG GRAPH 'N' ) ;
  237.   si d_dyne ;
  238.  l_t2 = prog ;
  239.  repeter bs (dime tabdyn . temps) ;
  240.   si (ega &bs 1) ; iterer bs ; finsi ;
  241.   l_t2 = l_t2 et (prog tabdyn . temps . (&bs - 1)) ;
  242.  fin bs ;
  243.  evfn1 = evol manu l_t2 tabdyn . force_de_choc ;
  244.  evux1 = evol manu l_t2 tabdyn . ux ;
  245.  TITRE '*****  LIAISON PROFIL_PROFIL_EXTERIEUR  *****';
  246.  DESS (( EVOL MANU 'TEMPS' TRESU1 . 'TEMPS_DE_SORTIE'
  247.       'FN'  TRESU1 . TL1 . 'FORCE_DE_CHOC') et (evfn1 coul bleu)) 
  248.                   'MIMA' ;
  249. *
  250.  DESS (( EVOL MANU 'TEMPS' TRESU1 . 'TEMPS_DE_SORTIE'
  251.       'UX'   TRESU1 . TL1 . 'UX' ) et (evux1 coul bleu)) 'MIMA';
  252. *
  253.   sinon ;
  254. *
  255.   l_t20 = prog tabdyn . temps . 1 ;
  256.   tt1 = 'EXTR' tabdyn . variables_internes . 1 'SLIA' 2 1 1 ;
  257.    lutt1 = dime tt1 . force_de_choc ;
  258.   l_fdc20 = extr tt1 . force_de_choc lutt1;
  259.   l_ux20 = extr tt1 . ux lutt1 ;
  260.  l_vi20 = prog 
  261.       (extr tabdyn . vitesses . 1 alfa tb1 . 1 . point_repere);
  262.   indi1 = 1 ;
  263.  
  264. repeter  bs ((dime tabdyn . temps) - 2) ;
  265.  indi1 = indi1 + 1 ;
  266.  ttu1 = 'EXTR' tabdyn . variables_internes . indi1 'SLIA' 2 1 1 ;
  267.  l_ux20 = l_ux20 et ttu1 . ux ;
  268.  l_fdc20 = l_fdc20 et ttu1 . force_de_choc ;
  269.   l_vi20 = l_vi20 et 
  270.   (prog (extr tabdyn . vitesses . indi1 alfa tb1 . 1 . point_repere));
  271.  l_t20 = l_t20 et (prog tabdyn . temps . indi1) ;
  272. fin bs ;
  273. evfdc20 = evol manu 'temps' l_t20 'force' l_fdc20 ;
  274. evux20 = evol manu 'temps' l_t20 'dep_pl' l_ux20 ;
  275. evit20 = evol manu 'temps' l_t20 'vitesse' l_vi20 ;
  276. evfdcr1 = evol manu 'temps' TRESU1 . 'TEMPS_DE_SORTIE'
  277.             'FN' TRESU1 . TL1 . 'FORCE_DE_CHOC' ;
  278. evuxr1 = EVOL MANU 'TEMPS' TRESU1 . 'TEMPS_DE_SORTIE'
  279.       'UX'   TRESU1 . TL1 . 'UX' ;
  280.  titre ' vitesse ' ;
  281.  DESS (evvi1 et (evit20 coul bleu)) 'MIMA' ;
  282.  titre 'force de choc' ;
  283.  DESS (( EVOL MANU 'TEMPS' TRESU1 . 'TEMPS_DE_SORTIE'
  284.       'FN'  TRESU1 . TL1 . 'FORCE_DE_CHOC') et (evfdc20 coul bleu)) 
  285.                   'MIMA' ;
  286.  titre ' ux liaison ' ;
  287.  DESS (( EVOL MANU 'TEMPS' TRESU1 . 'TEMPS_DE_SORTIE'
  288.       'UX'   TRESU1 . TL1 . 'UX' ) et (evux20 coul bleu)) 'MIMA';
  289. *
  290.   finsi ;
  291. FINSI ;
  292.  
  293.  
  294. * mess (extr tresu1 . temps_de_sortie 801) tabdyn . temps . 200 ;
  295. D_FC_1 = EXTRAIRE ( TRESU1 . TL1 . 'FORCE_DE_CHOC' ) 
  296.          801 ;
  297. D_DE_1 = EXTRAIRE ( TRESU1 . TL1 . 'UX' ) 801 ;
  298.  
  299. tt200 = 'EXTR' tabdyn . variables_internes . 200 'SLIA' 2 1 1 ;
  300. D_FC_2 =EXTRAIRE (tt200 . 'FORCE_DE_CHOC') 1;
  301. D_DE_2 =EXTRAIRE (tt200 . 'UX') 1;
  302.  
  303. *
  304. SAUT PAGE ;
  305. MESS '-----> liaison PROFIL_PROFIL_EXTERIEUR' ;
  306. MESS '           valeur de la force de choc :' D_FC_1 D_FC_2;
  307. MESS '                valeur du deplacement :' D_DE_1 D_DE_2;
  308. *
  309. R_FC = ABS ( ( D_FC_1 - D_FC_2 ) / D_FC_2 ) ;
  310. R_DE = ABS ( ( D_DE_1 - D_DE_2 ) / D_DE_2 ) ;
  311. R_DE0 = ABS((D_DE_1 - 1.29833E-03) / 1.29833E-03) ;
  312.  
  313. SI ( ( R_FC <EG (1.4497 * 1.2) ) ET ( R_DE <EG  (0.15738 *1.2) ) 
  314.      ET (R_DE0 < 1.e-2) ) ;
  315.    ERRE 0 ;
  316. SINON ;
  317.    ERRE 5 ;
  318. FINSI ;
  319.  
  320.  
  321. *   Liaison POINT_PLAN
  322. *
  323. NORM1 = 1  0  0 ;
  324. TL1  =  TABLE  'LIAISON_ELEMENTAIRE'    ;
  325. TL1.'TYPE_LIAISON'     =  MOT 'POINT_PLAN' ;
  326. TL1.'SUPPORT'          =      P2     ;
  327. TL1.'NORMALE'          =    NORM1    ;
  328. TL1.'RAIDEUR'          =    R_PPL    ;
  329. TL1.'JEU'              =    0.001    ;
  330. *
  331. TLB . 1  = TL1  ;
  332. TSOR2.TL1  = VRAI ;
  333. mapp2 = manu poi1 p2 ;
  334. mopp2 = modele mapp2 liaison POINT_PLAN ;
  335.  tuppe2 = table ; tuppe2 . mopp2 = vrai ;
  336. capp2 = mate mopp2 'NORM' norm1 'RAID' r_ppl 'JEU' 0.001 'SORT' tuppe2;
  337. *
  338. tabdyn2 = table 'PASAPAS' ;
  339. tabdyn2 . modele = mo1 et mopp2 ;
  340. tabdyn2 . caracteristiques = ca1 et capp2 ;
  341. tabdyn2 . temps = table ;
  342. tabdyn2 . temps . 0 = 0. ;
  343.  
  344. SI d_dyne ;
  345. tabdyn2 . deplacements = table ;
  346. tabdyn2 . deplacements . 0 =  ALO ;
  347. tabdyn2 . sortie = table ;
  348. tabdyn2 . sortie . mopp2 = vrai ;
  349. tabdyn2 . sortie . vitesses = faux ;
  350. tabdyn2 . pas_de_temps = pdt ;
  351. tabdyn2 . nombre_pas = np1 ;
  352. tabdyn2 . pas_de_sortie = nins ;
  353.   dyne de_vogelaere tabdyn2 ;
  354.  
  355. SINON ;
  356. tabdyn2 . deplacements = table ;
  357. tabdyn2 . deplacements . 0 =  t_prec . deplacements . 1 ;
  358. tabdyn2 . contraintes = table ;
  359. tabdyn2 . contraintes . 0 = t_prec . contraintes . 1 ;
  360. tabdyn2 . deformations = table;
  361. tabdyn2 . deformations . 0 = t_prec . deformations . 1;
  362. tabdyn2 . dynamique = vrai ;
  363. tabdyn2 . temps_calcules = prog 0 pas (4*nins*pdt) (np1*pdt);
  364. *tabdyn2 . temps_calcules = prog 0 pas (pdt) .0005;
  365. *tabdyn2 . temps_sauves = prog 0 pas (nins*pdt) (np1*pdt) ;
  366. *tabdyn2 . temps_sauves = prog 0 pas (nins*pdt) .01 ;
  367. tabdyn2 . hypothese_deformations = 'LINEAIRE' ;
  368.   pasapas tabdyn2 ; 
  369.  
  370. FINSI ;
  371.  
  372. *np1 = 5 ; nins = 1;
  373. TRESU2 = DYNE 'DE_VOGELAERE' TBAS TA TINIT TSORT TLIA 
  374.          NP1 PDT NINS ;
  375.  
  376.  nsu2 = dime tresu2 . temps_de_sortie ;
  377.  ovi2 = prog ;
  378.  repeter bsu2 nsu2 ;
  379.    ovi2 = ovi2 et 
  380.   (prog (extr tresu2 . &bsu2 . vitesse alfa tb1 . 1 . point_repere)) ;
  381.  fin bsu2 ;
  382.  evvi2 = evol manu 'temps' TRESU2 . 'TEMPS_DE_SORTIE'
  383.       'vitesse' ovi2 ;
  384. *
  385. *     Trace des courbes
  386. *
  387. SI ( NEG GRAPH 'N' ) ;
  388.   si d_dyne ;
  389. l_t22 = prog ;
  390. repeter bs (dime tabdyn2 . temps) ;
  391.  si (&bs ega 1) ; iterer bs ; finsi ;
  392.  l_t22 = l_t22 et (prog tabdyn2 . temps . (&bs - 1)) ;
  393. fin bs ;
  394.  evfn2 = evol manu l_t22 tabdyn2 . force_de_choc;
  395.  evux2 = evol manu l_t22 tabdyn2 . ux ;
  396.    TITRE '*****  LIAISON POINT_PLAN  *****' ;
  397.    DESS (( EVOL MANU 'TEMPS' TRESU2 . 'TEMPS_DE_SORTIE'
  398.        'FN'    TRESU2 . TL1 . 'FORCE_DE_CHOC') et (evfn2 coul
  399. vert)) 'MIMA' ;
  400. *
  401.    DESS (( EVOL MANU 'TEMPS' TRESU2 . 'TEMPS_DE_SORTIE'
  402.        'UX'    TRESU2 . TL1 . 'UX' ) et (evux2 coul vert)) 'MIMA' ;
  403.   sinon ;
  404. *
  405.   l_t202 = prog tabdyn2 . temps . 1 ;
  406.   tt2 = 'EXTR' tabdyn2 . variables_internes . 1 'SLIA' 2 1 1 ;
  407.    lutt2 = dime tt2 . force_de_choc ;
  408.   l_fdc202 = extr tt2 . force_de_choc lutt2;
  409.   l_ux202 = extr tt2 . ux lutt2 ;
  410.  l_vi202 = prog 
  411.       (extr tabdyn2 . vitesses . 1 alfa tb1 . 1 . point_repere);
  412.   indi1 = 1 ;
  413.  
  414. repeter  bs ((dime tabdyn2 . temps) - 2) ;
  415.  indi1 = indi1 + 1 ;
  416.  ttu2 = 'EXTR' tabdyn2 . variables_internes . indi1 'SLIA' 2 1 1 ;
  417.  l_ux202 = l_ux202 et ttu2 . ux ;
  418.  l_fdc202 = l_fdc202 et ttu2 . force_de_choc ;
  419.   l_vi202 = l_vi202 et 
  420.   (prog (extr tabdyn2 . vitesses .indi1 alfa tb1 . 1 . point_repere)) ;
  421.  l_t202 = l_t202 et (prog tabdyn2 . temps . indi1) ;
  422. fin bs ;
  423. evfdc202 = evol manu 'temps' l_t202 'force' l_fdc202 ;
  424. evux202 = evol manu 'temps' l_t202 'dep_pl' l_ux202 ;
  425. evit202 = evol manu 'temps' l_t202 'vitesse' l_vi202 ;
  426. evfdcr2 = evol manu 'temps' TRESU2 . 'TEMPS_DE_SORTIE'
  427.             'FN' TRESU2 . TL1 . 'FORCE_DE_CHOC' ;
  428. evuxr2 = EVOL MANU 'TEMPS' TRESU2 . 'TEMPS_DE_SORTIE'
  429.       'UX'   TRESU2 . TL1 . 'UX' ;
  430.  titre ' vitesse ' ;
  431.  DESS (evvi2 et (evit202 coul bleu)) 'MIMA' ;
  432.  titre 'force de choc' ;
  433.  DESS (( EVOL MANU 'TEMPS' TRESU2 . 'TEMPS_DE_SORTIE'
  434.       'FN'  TRESU2 . TL1 . 'FORCE_DE_CHOC') et (evfdc202 coul bleu)) 
  435.                   'MIMA' ;
  436.  titre ' ux liaison ' ;
  437.  DESS (( EVOL MANU 'TEMPS' TRESU2 . 'TEMPS_DE_SORTIE'
  438.       'UX'   TRESU2 . TL1 . 'UX' ) et (evux202 coul bleu)) 'MIMA';
  439. *
  440.   finsi ;
  441. FINSI ;
  442. *
  443. TEMPS ;
  444. *
  445. *     Code de bon fonctionnement
  446. *
  447. D_FC_20 =EXTRAIRE ( TRESU2 . TL1 . 'FORCE_DE_CHOC' ) 801;
  448. D_DE_20 =EXTRAIRE ( TRESU2 . TL1 . 'UX' ) 801 ;
  449. tt2200 = 'EXTR' tabdyn2 . variables_internes . 200 'SLIA' 2 1 1 ;
  450. D_FC_21 =EXTRAIRE (tt2200 . 'FORCE_DE_CHOC') 1;
  451. D_DE_21 =EXTRAIRE (tt2200 . 'UX') 1;
  452.  
  453. SAUT 2 LIGNES ;
  454. MESS '-----> liaison POINT_PLAN ' ;
  455. MESS '           valeur de la force de choc :' D_FC_20 D_FC_21 ;
  456. MESS '                valeur du deplacement :' D_DE_20 D_DE_21 ;
  457. SAUT 3 LIGNES ;
  458. R_FC = ABS ( ( D_FC_20 - D_FC_21 ) / D_FC_20 ) ;
  459. R_DE = ABS ( ( D_DE_20 - D_DE_21 ) / D_DE_20 ) ;
  460. R_DE0 = ABS ((D_DE_20 - 1.29833E-03) / 1.29833E-03) ;
  461. SI ( ( R_FC &lt;EG (0.59178 *1.2)) ET ( R_DE &lt;EG (0.13598 * 1.2)) 
  462.     ET (R_DE0 < 1.e-2) ) ;
  463.    ERRE 0 ;
  464. SINON ;
  465.    ERRE 5 ;
  466. FINSI ;
  467. *
  468.  
  469. FIN ;
  470.  
  471.  
  472.  
  473.  
  474.  
  475.  
  476.  
  477.  
  478.  
  479.  
  480.  
  481.  
  482.  
  483.  
  484.  
  485.  
  486.  
  487.  
  488.  
  489.  

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