Télécharger soudage18.dgibi

Retour à la liste

Numérotation des lignes :

  1. *
  2. graph = faux ;
  3. complet = faux ;
  4.  
  5. * donnees geometriques qu il faudrait sauver
  6. opti dime 3 elem cub8 trac x;
  7. p_ori = 0. 0. 0. ; d0 = 1.d-2 ;
  8.  
  9. e_x = 1. 0. 0. ; e_y = 0. 1. 0. ; e_z = 0. 0. 1. ;
  10. l0_x = 5.d-3 ;l_x = 1.d-2 ; l_y = 10*l_x ; l_zh = 3*l_x ; l_zb = 2*l_x ;
  11. l_cz = 0.5*l_x ; l_cy = 1.05*l_x ;
  12. dens d0 ; n_x = 1 ;
  13.  
  14. *******************
  15. debproc prosym vut1*maillage ;
  16. * surface symetrie
  17. chu1 = coor 2 (manu chml (enve vut1) scal 1.) ;
  18. ssu1 = elem chu1 inferieur (d0*1.e-2) ;
  19. chu2 = coor 2 (manu chml ssu1 scal 1.) ;
  20. ssu2 = elem chu2 superieur (d0*1.e-1) ;
  21. ssu3 = diff ssu1 ssu2 ;
  22. finproc ssu3 ;
  23. *******************
  24.  
  25. **
  26. *si d_ther ;
  27. *opti sauv 'soudage18_lin_a.sauv' ;
  28. **
  29.  
  30. * plaque
  31. pc1m = p_ori moins (l0_x*e_x) ;
  32. pp1m = pc1m moins (l_zb*e_z) ;
  33. pc2m = pc1m plus (l_cy*e_y) plus (l_cz*e_z) ;
  34. pp2m = pp1m plus (l_cy*e_y) ;
  35. pp3m = pp1m plus (l_y*e_y);
  36. pp4m = pp3m plus ((l_zh + l_zb)*e_z) ;
  37. pc3m = pc2m plus (l_cz*e_z) ;
  38. pp5m = pp3m plus ((l_zb + l_cz)*e_z) ;
  39. pp6m = pp5m plus (l_cz*e_z) ;
  40. pp7m = pp2m plus ((l_zh + l_zb)*e_z) ;
  41.  
  42.  
  43. l1c = d pc1m pp1m ; nncz1 = nbel l1c ;
  44. l2c = d pp1m pp2m ; nnc2 = nbel l2c ;
  45. l3c = d nncz1 pp2m pc2m ;
  46. l4c = d nnc2 pc2m pc1m ;
  47. s1c = dall l1c l2c l3c l4c ;
  48.  
  49. l11 = d pp2m pp3m ;
  50. l12 = d nncz1 pp3m pp5m ; l13 = d pp5m pc2m ;
  51. s11 = dall l11 l12 l13 (inve l3c) ;
  52. l21 = d pc3m pc2m ; nncz2 = nbel l21 ;
  53. l23 = d nncz2 pp5m pp6m ; l24 = d pp6m pc3m ;
  54. s22 = dall l21 (inve l13) l23 l24 ;
  55. l31 = d pp7m pc3m ; nncz3 = nbel l31 ;
  56. l33 = d nncz3 pp6m pp4m ; l34 = d pp4m pp7m ;
  57. s33 = dall l31 (inve l24) l33 l34 ;
  58.  
  59. *trac (s1c et s11 et s22 et s33) ;
  60. sp0 = s1c et s11 et s22 et s33 ;
  61. vp0 = sp0 volu n_x trans (l_x*e_x) ;
  62. vpp1 = (face 2 vp0) volu n_x trans (l_x*e_x) ;
  63. vpm1 = sp0 volu n_x trans (-1*l_x*e_x) ;
  64. vv1 = vpm1 et vp0 et vpp1 ;
  65. *trac vv1 cach ;
  66.  
  67. * capteurs thermiques
  68. ppthm1 = pp1m plus (3*l_x*e_y) ;
  69. ppthm2 = ppthm1 plus ((l_zh + l_zb)*e_z) ;
  70.  
  71. *apport
  72. pa1m pa2m pa3m la4c la21 = pc1m pc2m pc3m l4c l21 plus p_ori ;
  73. pa4m = pa1m plus (l_x*e_z) ;
  74. la22 = d pa3m pa4m ; la23 = d pa4m pa1m ;
  75. sa0 = surf plan ((inve la21) et la22 et la23 et (inve la4c)) ;
  76. va0 = sa0 volu n_x trans (l_x*e_x) ; va0 = va0 coul rouge ;
  77.  
  78. *trac (vv1 et va0) cach ;
  79.  
  80. ********************
  81. T_INIT = 20. ;
  82. T_CONV = T_INIT ;
  83. T_FUSION = 1400. ;
  84. t_dep = 1600. ;
  85. ********************
  86.  
  87. mothep = mode vv1 thermique conduction ;
  88. motha0 = mode va0 thermique conduction ;
  89.  
  90. * materiau
  91. TTMP = TABL ;
  92. TTMP. 1 = PROG 20 200 400 600 800 1000 1200 1400 1500 3000 ;
  93. TTMP. 2 = PROG 8.0e3 7.93e3 7.84e3 7.75e3 7.65e3 7.55e3
  94. 7.45e3 7.35e3 7.30e3 7.30e3 ;
  95. EVRHO = EVOL MANU 'T' (TTMP. 1) 'RHO' (TTMP. 2) ;
  96.  
  97. * Parametres thermiques
  98. TTMP = TABL ;
  99. TTMP. 1 = PROG 20 100 200 300 400 500 600 700 800 900 1000
  100. 1200 1400 3000 ;
  101. TTMP. 2 = PROG 0.0147e3 0.0158e3 0.0172e3 0.0186e3 0.02e3
  102. 0.0211e3 0.0222e3 0.0232e3 0.0241e3 0.0248e3 0.0255e3
  103. 0.0269e3 0.0283e3 0.0283e3 ;
  104. EVK = EVOL MANU 'T' (TTMP. 1) 'K' (TTMP. 2);
  105.  
  106. TTMP = TABL ;
  107. TTMP. 1 = PROG 20 100 200 300 400 600 800 1000 1200 1500 3000 ;
  108. TTMP. 2 = PROG 450 490 525 545 560 580 625 660 670 690 690 ;
  109. EVCP = EVOL MANU 'T' (TTMP. 1) 'C' (TTMP. 2);
  110.  
  111. SI GRAPH ;
  112. DESS EVRHO TITR 'Masse Volumique';
  113. DESS EVK TITR 'Conductivite Thermique';
  114. DESS EVCP TITR 'Chaleur Specifique';
  115. FINS;
  116.  
  117. * coefficient de convection avec l'exterieur
  118. he1 = 20. ;
  119.  
  120.  
  121. * MATERIAUX MECANIQUE ET THERMIQUE
  122.  
  123. CATHEP = MATE MOTHEP 'K' EVK 'RHO' EVRHO 'C' EVCP ;
  124.  
  125. cathea0 = mate motha0 'K' EVK 'RHO' EVRHO 'C' EVCP ;
  126.  
  127. **********************************
  128. * procede
  129.  
  130. vpu1 = vv1 ;
  131. mothu1 = mothep ;
  132. cathu1 = cathep ;
  133.  
  134. symvpu1 = prosym vpu1 ;
  135. symva0 = prosym va0 ;
  136. * trac (symvpu1 et symva0) ;
  137.  
  138. * surface convection
  139.  
  140. elim (d0*1.d-3) vpu1 va0 ;
  141. sconvu1 = enve vpu1 ;
  142. sdep1 = inte sconvu1 (enve va0) ;
  143. sconvu1 = diff sconvu1 sdep1 ;
  144. sconvu1 = diff sconvu1 symvpu1 ;
  145. *trac sconvv1 ;
  146.  
  147. * convection
  148. ctconvv1 = manu chpo sconvu1 T 20. ;
  149. ev_conv = evol manu 't (s)' (prog 0. 1.e5) 'Te (oC)' (prog 1. 1.) ;
  150. cha_conv = char teco ev_conv ctconvv1 ;
  151.  
  152. moconvu1 = mode sconvu1 thermique convection ;
  153. caconvu1 = mate moconvu1 'H' he1 ;
  154.  
  155. * temperatures initiale
  156. mpp1_u = diff (changer poi1 vpu1) (changer poi1 sdep1) ;
  157. cptempi = manu chpo 1 mpp1_u T t_init nature diffuse ;
  158.  
  159. ca0tempi = manu chpo 1 (changer poi1 va0) T t_dep nature diffuse ;
  160.  
  161. si graph ;
  162. trac (cptempi et ca0tempi) (vpu1 et va0)
  163. titre 'temperatures initiales' nclk;
  164. finsi ;
  165. * deplacements initiaux
  166. ma0p1_u = diff (changer poi1 va0) (changer poi1 sdep1) ;
  167. ca0depi = manu chpo 3 ma0p1_u ux 0. uy 0. uz 0. nature diffuse ;
  168.  
  169.  
  170.  
  171. tdep = table 'PASAPAS' ; tdep.'PROCESSEURS' = MOT 'AUTOMATIQUE';
  172. tdep . modele = mothu1 et motha0 et moconvu1 ;
  173. tdep . caracteristiques = cathu1 et cathea0 et caconvu1 ;
  174. tdep . chargement = cha_conv ;
  175. tdep . procedure_thermique = DUPONT ;
  176.  
  177. tdep . temperatures = table ;
  178. tdep . temperatures . 0 = cptempi et ca0tempi ;
  179.  
  180. si complet ;
  181. tdep . temps_calcules = prog 0. pas 1. 10. pas 10. 100. pas 1.e2 1.e3;
  182. tdep . temps_sauves = prog 0. 1. 5. 10. 20. 50. 100. 200. 500. 1.e3 ;
  183. sinon ;
  184. tdep . temps_sauves = prog 0. 1. ;
  185. finsi ;
  186.  
  187. pasapas tdep ;
  188.  
  189. ****************************************
  190. mailtra1 = enve (extr (mothu1 et motha0) mail) ;
  191. ppthm1 = point mailtra1 proc ppthm1 ; op1 = prog ;
  192. ppthm2 = point mailtra1 proc ppthm2 ; op2 = prog ;
  193. abt = prog ;
  194. * temperature
  195. l_temp = prog 100. 200. 300. 400. 500. 600. 700. 800. 900. 1000.
  196. 1100. 1200. 1300. 1400. ;
  197. si graph ;
  198. repeter bdes (dime tdep . temps) ;
  199. ind1 = &bdes - 1 ;
  200. abt = abt et (prog tdep . temps . ind1) ;
  201. op1 = op1 et (prog (extr tdep . temperatures . ind1 t ppthm1)) ;
  202. op2 = op2 et (prog (extr tdep . temperatures . ind1 t ppthm2)) ;
  203.  
  204. titre 's18_lin_a / temperature / date = ' tdep . temps . ind1 ;
  205. *trac tdep . temperatures . ind1 mailtra1 (arete mailtra1) l_temp ;
  206. trac tdep . temperatures . ind1 mailtra1 l_temp nclk ;
  207. fin bdes ;
  208.  
  209. titre 's18_lin_a / thermocouples ' ;
  210. evcapt = (evol manu abt 't(s)' op1 'T(°C)' coul jaune)
  211. et (evol manu abt 't(s)' op2 'T(°C)' coul rouge) ;
  212. dess evcapt ;
  213. finsi ;
  214. **
  215. *d_ther = faux ;
  216. *sauv tdep;
  217. *finsi ;
  218. **
  219.  
  220. ********************
  221. TAREF = -273.15 ;
  222. ********************
  223. *******************
  224. debpro prpmec tu1*table ;
  225. *calcul mecanique
  226.  
  227. tu1 . blocages_mecaniques = (bloq ux pp3mm1) et
  228. (bloq uz (pp3mm1 et pp3mp2)) et (bloq (symvpu1 et symva0) uy);
  229. tu1 . chargement = ch_temp ;
  230. tu1 . talpha_reference = TAREF ;
  231. tu1 . temps_calcules = tdep . temps_sauves ;
  232.  
  233. pasapas tu1 ;
  234.  
  235. opsy1 = prog ; abst = prog ;
  236.  
  237. repeter beliso (dime tu1 . temps) ;
  238. ind2 = &beliso - 1 ;
  239. si graph ;
  240. titre 's18_lin_a / T /' tu1 . temps . ind2 ' s';
  241. trac tdep . temperatures . ind2 mailtra1 nclk;
  242. titre (et tu1 . etiquette ' / vmises /') tu1 . temps . ind2 ' s';
  243. trac tu1 . modele tu1 . contraintes . ind2 ;
  244. finsi ;
  245. chpu5 = changer chpo tu1 . modele (exco tu1 . contraintes . ind2 smyy);
  246. opsy1 = opsy1 et (prog ((extr chpu5 smyy ppmem1)/1.e6)) ;
  247. abst = abst et (prog tu1 . temps . ind2) ;
  248. fin beliso ;
  249.  
  250. tu1 . ecsmy1 = evol manu abst 't(s)' opsy1 ' smyy (MPa)' ;
  251. si graph ;
  252. titre 's18_lin_a / synthese SMYY' ;
  253. dess tu1 . ecsmy1 nclk;
  254. finsi ;
  255. finproc ;
  256. *******************
  257.  
  258. *******************
  259. *Procédure pour calculer YOUN à partir d'une évolution
  260. 'DEBP' CALYO EV1*EVOLUTION;
  261. X1 = EXTR EV1 'ABSC' ;
  262. Y1 = EXTR EV1 'ORDO' ;
  263. YOU1 = '/' (EXTR Y1 2) (EXTR X1 2) ;
  264. 'FINP' YOU1 ;
  265. *******************
  266.  
  267. **
  268. *si (non d_ther) ;
  269. *opti rest 'soudage18_lin_a.sauv' ;
  270. *rest ;
  271. **
  272.  
  273. *************
  274. * donnees materiau
  275. TTMP = TABL ;
  276. TTMP. 1 = PROG 20. 200. 300. 400. 500. 600.
  277. 700. 800. 900. 1000. 1500. 3000. ;
  278. TTMP. 2 = PROG 197e9 184e9 176e9 168e9
  279. 160e9 151e9 142e9 130e9 108e9 81e9 1e9 1e9 ;
  280. EVYOUN = EVOL MANU 'T' (TTMP. 1) 'YOUN' (TTMP. 2) ;
  281. MAXIT = MAXI (TTMP. 1) ;
  282.  
  283.  
  284. TTMP = TABL ;
  285. TTMP. 1 = PROG 25 100 200 400 500 600 700 800
  286. 900 1000 1100 1200 1300 1390 1420 1450 3000 ;
  287. TTMP. 2 = PROG 0.296 0.298 0.304 0.315 0.32
  288. 0.323 0.326 0.33 0.336 0.339 0.346 0.349 0.353 0.353 0.353 0.353 0.353 ;
  289. DIMNU = DIME TTMP. 1 ;
  290. VALNU = TTMP. 2 ;
  291. EVNU = EVOL MANU 'T' (TTMP. 1) 'NU' (TTMP. 2) ;
  292.  
  293.  
  294. TTMP = TABL ;
  295. TTMP. 1 = PROG 20 200 300 400 500 600 700 800 900 1000 1500 3000 ;
  296. TTMP. 2 = PROG 287e6 198e6 172e6 157e6 152e6 145e6 136e6
  297. 127e6 115e6 79e6 5e6 1e6 ;
  298. EVSIGY = EVOL MANU 'T' (TTMP. 1) 'R0' (TTMP. 2) ;
  299.  
  300. TTMP = TABL ;
  301. TTMP. 1 = PROG 0 100 200 300 400 500 600 700 800 900 1000 1100
  302. 1200 1420 1460 3000 ;
  303. TTMP. 2 = PROG -3.17e-4 0.0013 0.003 0.0048 0.0066 0.0086 0.0106
  304. 0.0126 0.0147 0.0168 0.019 0.0212 0.0235 0.029 0.0298 0.0298 ;
  305. EVDEFT = EVOL MANU 'T' (TTMP. 1) 'EPTH' (TTMP. 2) ;
  306. LALPHA = TTMP. 2 / (TTMP. 1 - TAREF) ;
  307. EVALPHA = EVOL MANU 'T' (TTMP. 1) 'ALPH' LALPHA ;
  308.  
  309. TTMP1 = TABL ;
  310. TTMP1. 1 = PROG 20 200 300 400 500 600 700 800 900 1000 1500 3000 ;
  311. TTMP2 = TABL ;
  312. TTMP2. 1 = PROG 0. 0.01 0.05 0.2 1. 2. ;
  313. TTMP2. 2 = PROG 0 28e6 132e6 274e6 913e6 1369e6 ;
  314. TTMP2. 3 = PROG 0 33e6 136e6 272e6 802e6 1203e6 ;
  315. TTMP2. 4 = PROG 0 29e6 135e6 279e6 728e6 1092e6 ;
  316. TTMP2. 5 = PROG 0 29e6 137e6 278e6 442e6 663e6 ;
  317. TTMP2. 6 = PROG 0 27e6 130e6 267e6 398e6 597e6 ;
  318. TTMP2. 7 = PROG 0 25e6 117e6 242e6 355e6 532e6 ;
  319. TTMP2. 8 = PROG 0 25e6 76e6 204e6 264e6 396e6 ;
  320. TTMP2. 9 = PROG 0 25e6 72e6 167e6 223e6 334e6 ;
  321. TTMP2. 10 = PROG 0 22e6 46e6 54e6 65e6 97e6 ;
  322. TTMP2. 11 = PROG 0 3e6 18e6 21e6 26e6 39e6 ;
  323. TTMP2. 12 = PROG 0 2e6 3e6 4e6 5e6 6e6 ;
  324. TTMP2. 13 = PROG 0 1e6 2e6 3e6 4e6 5e6 ;
  325.  
  326. EVTOT = VIDE 'EVOLUTION';
  327. LTEMPER = VIDE 'LISTREEL' ;
  328. LYOU2 = VIDE 'LISTREEL' ;
  329. LYOU3 = VIDE 'LISTREEL' ;
  330. REPE ITEMPER (DIME TTMP1. 1);
  331. J=&ITEMPER;
  332. TEMPER = EXTR (TTMP1 . 1) J;
  333. LTEMPER= LTEMPER ET TEMPER ;
  334.  
  335. * Interpolation du module de YOUNG et SIGY a la temperature TEMPER
  336. YOUNIPOL = IPOL TEMPER EVYOUN;
  337. SIGYIPOL = IPOL TEMPER EVSIGY;
  338. LYOU2 = LYOU2 ET YOUNIPOL ;
  339.  
  340. * Calcul de EPSIY a partir de YOUNIPOL et SIGYIPOL
  341. EPSIY = SIGYIPOL / YOUNIPOL;
  342.  
  343. * Construction de la courbe de traction Complete pour Cast3M
  344. LESPI = (PROG 0.) ET ((TTMP2. 1) + EPSIY );
  345. LSIGM = (PROG 0.) ET ((TTMP2. (J+1)) + SIGYIPOL);
  346. EVTRAC= EVOL MANU 'EPS' LESPI 'CONT' LSIGM ;
  347. EVTOT = EVTOT ET EVTRAC;
  348. LYOU3 = LYOU3 ET (CALYO EVTRAC) ;
  349.  
  350. * Construction du NUAGE qui associe une EVOLUTION a une TEMPERATURE
  351. SI (J EGA 1);
  352. NUATRAC = NUAGE 'COMP' 'T ' TEMPER
  353. 'COMP' 'TRAC' EVTRAC ;
  354. SINO;
  355. NUATRAC = NUATRAC ET (NUAGE 'COMP' 'T ' TEMPER
  356. 'COMP' 'TRAC' EVTRAC ) ;
  357. FINS;
  358. FIN ITEMPER;
  359.  
  360. EVYOUN = EVOL MANU 'T' LTEMPER 'YOUN' LYOU3 ;
  361. evlimp = EVOL MANU 'T' LTEMPER 'LIMP' (LYOU3*2.D-3) ;
  362.  
  363. SI GRAPH ;
  364. DESS EVYOUN TITR 't_initModule de YOUNG';
  365. DESS EVNU TITR 'Coefficient de POISSON';
  366. DESS EVSIGY TITR 'Limite Elastique';
  367. DESS EVALPHA TITR 'Coefficient de dilatation Thermique';
  368. DESS EVTOT TITR 'Courbe de Traction a differentes temperatures';
  369. dess evlimp TITR 'Limite 1er invariant ' ;
  370. FINS;
  371. *************
  372.  
  373. * point d appui
  374. l1cm1 pp3mm1 = l1c pp3m moins (l_x*e_x) ;
  375. l1cp2 pp3mp2 = l1c pp3m plus (2*l_x*e_x) ;
  376. * symetrie
  377. symvpu1 = prosym vpu1 ;
  378. symva0 = prosym va0 ;
  379. * trac (symvpu1 et symva0) ;
  380. * chargement
  381. ch_temp = char t (tdep . temps) (tdep . temperatures) ;
  382.  
  383. **
  384. momecp = mode vpu1 'MECANIQUE' 'ELASTIQUE' 'ISOTROPE'
  385. 'PLASTIQUE' 'ISOTROPE' ;
  386. momca0 = mode va0 'MECANIQUE' 'ELASTIQUE' 'ISOTROPE'
  387. 'PLASTIQUE' 'ISOTROPE' ;
  388. *
  389. momecb = mode (vpu1 et va0) 'MECANIQUE' 'ELASTIQUE' 'ISOTROPE'
  390. 'PLASTIQUE' 'ISOTROPE' 'BBAR' ;
  391. *
  392. momecl = mode (vpu1 et va0) 'MECANIQUE' 'ELASTIQUE' 'ISOTROPE'
  393. 'PLASTIQUE' 'ISOTROPE' 'CRIP' ;
  394. *
  395. momect =mode (vpu1 et va0) 'MECANIQUE' 'ELASTIQUE' 'ISOTROPE'
  396. 'PLASTIQUE' 'ISOTROPE' 'CRIP' 'BBAR' ;
  397.  
  398. **
  399. camecp = MATE MOMECP 'YOUN' EVYOUN 'NU' EVNU 'RHO' EVRHO 'ALPH' EVALPHA
  400. 'TRAC' NUATRAC ;
  401. *
  402. camca0 = MATE MOMCA0 'YOUN' EVYOUN 'NU' EVNU 'RHO' EVRHO 'ALPH' EVALPHA
  403. 'TRAC' NUATRAC ;
  404. *
  405. camecb = MATE MOMECB 'YOUN' EVYOUN 'NU' EVNU 'RHO' EVRHO 'ALPH' EVALPHA
  406. 'TRAC' NUATRAC ;
  407. *
  408. camecl = MATE MOMECL 'YOUN' EVYOUN 'NU' EVNU 'RHO' EVRHO 'ALPH' EVALPHA
  409. 'TRAC' NUATRAC 'LIMP' evlimp ;
  410. *
  411. camect = MATE MOMECT 'YOUN' EVYOUN 'NU' EVNU 'RHO' EVRHO 'ALPH' EVALPHA
  412. 'TRAC' NUATRAC 'LIMP' evlimp ;
  413. **
  414. mailtra1 = vpu1 et va0 ;
  415. *trac mailtra1 ;
  416. mm125 = elem 125 mailtra1 ;
  417. myu1 = elem mailtra1 appuye largement va0 ;
  418. myu1 = diff myu1 va0 ;
  419.  
  420.  
  421. pp3mm1 = point mailtra1 proc pp3mm1 ;
  422. pp3mp2 = point mailtra1 proc pp3mp2 ;
  423. ppmem1 = point mailtra1 proc ppthm1 ;
  424.  
  425. tmec0 = table 'PASAPAS' ;
  426. tmec0 . modele = momecp et momca0 ;
  427. tmec0 . caracteristiques = camecp et camca0 ;
  428. tmec0 . etiquette = 's18_lin_a / standard' ;
  429. *tu1 = tmec0 ;
  430. prpmec tmec0 ;
  431.  
  432. tmecb = table 'PASAPAS' ;
  433. tmecb . modele = momecb ;
  434. tmecb . caracteristiques = camecb ;
  435. tmecb . etiquette = 's18_lin_a / bbar' ;
  436. prpmec tmecb ;
  437.  
  438. tmecl = table 'PASAPAS' ;
  439. tmecl . modele = momecl ;
  440. tmecl . caracteristiques = camecl ;
  441. tmecl . etiquette = 's18_lin_a / critp' ;
  442. prpmec tmecl ;
  443.  
  444. tmect = table 'PASAPAS' ;
  445. tmect . modele = momect ;
  446. tmect . caracteristiques = camect ;
  447. tmect . etiquette = 's18_lin_a / critp et bbar' ;
  448. prpmec tmect ;
  449.  
  450. tsymb = table ; tsymb . titre = table ;
  451. tsymb . 1 = 'MARQ PLUS' ; tsymb . titre . 1 = 'stand' ;
  452. tsymb . 2 = 'MARQ CARR' ; tsymb . titre . 2 = 'bbar' ;
  453. tsymb . 3 = 'MARQ ETOI' ; tsymb . titre . 3 = 'critp' ;
  454. tsymb . 4 = 'MARQ LOSA' ; tsymb . titre . 4 = 'bbar et critp' ;
  455.  
  456. si graph ;
  457. titre ' s18_lin_a / capteurs contraites ';
  458. dess (tmec0 . ecsmy1 et tmecb . ecsmy1 et tmecl . ecsmy1
  459. et tmect . ecsmy1) tsymb 'LEGE' ;
  460. finsi ;
  461.  
  462. **
  463. *finsi ;
  464. **
  465.  
  466. *test octobre 2018
  467. or0 = extr (extr tmec0 . ecsmy1 ordo) 2 ;
  468. e0 = abs ((or0 - 7.20114E+01) / 7.20114E+01) ;
  469. orb = extr (extr tmecb . ecsmy1 ordo) 2 ;
  470. eb = abs ((orb - 1.38370E+01) / 1.38370E+01) ;
  471. orl = extr (extr tmecl . ecsmy1 ordo) 2 ;
  472. el = abs ((orl - 1.02034E+01) / 1.02034E+01) ;
  473.  
  474. si ((e0 > 1.e-3) ou (eb > 1.e-3) ou (el > 1.e-3)) ;
  475. erre 5 ;
  476. sinon ;
  477. erre 0 ;
  478. finsi ;
  479.  
  480. fin ;
  481.  
  482.  
  483.  
  484.  
  485.  

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