Télécharger la plaquette Cast3M


fr en

Télécharger sicfsic.dgibi

Retour à la liste

Numérotation des lignes : 

  1. * fichier :  sicfsic.dgibi
  2. ************************************************************************
  3. ************************************************************************
  4. complet = faux ;
  5. * POUR UN CALCUL COMPLET METTRE À VRAI
  6. *=======================================================================
  7. *
  8. *      CAS TEST DE VALIDATION DES LOIS DE COMPORTEMENT ONERA  
  9. *          SCALAIRE ET PSEUDO-TENSORIEL POUR LE COMPOSITE 
  10. *              TISSE SICf/SIC. 
  11. *        Paramètres pris en compte correspondant au tissé
  12. *        Modèle décrit dans la note SEMT/LM2S/RT/05-034/A. 
  13. *
  14. *=======================================================================
  15. *
  16. *    1 ELEMENT CUB 8 : COMPARAISON RESULTATS AVec ceux obtenu par 
  17. *    fichier MATLAB de l'ONERA  
  18. *    Traction pour le cas Complet = faux 
  19. *    Cycles de traction compressionspour Complet = Vrai 
  20. *=======================================================================
  21. GRAPH = N; 
  22. opti echo 1 dime 3 elem cub8 ; 
  23. OPTI TRAC PSC ;    
  24. *  Maillage   
  25. long = 1.0 ;
  26. larg = 1.0 ;
  27. epai = 1.0 ;
  28. p1  = 0.0 0.0 0.0 ;
  29. p2  = p1 PLUS ( 0.0 larg 0.0) ;
  30. p3  = p2 PLUS (0.0  0.0 long) ;
  31. p4  = p1 PLUS (0.0  0.0 long) ;
  32. p11 = p1 PLUS (epai 0.0  0.0 ) ;
  33. p21 = p2 PLUS (epai 0.0  0.0 ) ;
  34. p31 = p3 PLUS (epai 0.0  0.0 ) ;
  35. p41 = p4 PLUS (epai 0.0  0.0 ) ;
  36. n1 = 1 ; 
  37. p1p2 = DROIT N1 p1 p2 ;
  38. p4p3 = DROIT N1 p4 p3 ;
  39. sbas = p1p2 REGL N1 p4p3 ;
  40. p11p21 = DROIT N1 p11 p21 ;
  41. p41p31 = DROIT N1 p41 p31 ;
  42. shaut = p11p21 REGL N1 p41p31 ;
  43. cube = sbas VOLU N1 shaut ;
  44. vol1 = cube ; 
  45. *
  46. * Cas où on veut atteindre 180 puis 235 MPa pour fibre
  47. * selon la direction de la traction 
  48. * Différent numéros de Cachou selon les angles
  49. *
  50. * tenten = vrai pour le modèle pseudo tensoriel faux sinon
  51. TENTEN = FAUX; 
  52. * On fait une boucle ou on teste successivemenbt le modèle 
  53. * Scalaire et tensoriel ; 
  54. REPETER BOUMO 2 ; 
  55. SI (&BOUMO EGA 2) ; 
  56. TENTEN = VRAI ; 
  57. FINSI ; 
  58. *
  59. *  On tire à 25 ° selon x 
  60. TETE = 22.5 ; 
  61. SI TENTEN ; 
  62. SI COMPLET ; 
  63. L_depl=PROG 0.0  PAS 0.01 0.587 PAS -0.01 -1.13 PAS 0.01 0.65
  64. PAS -0.01 -1.13 PAS 0.01 0.725 PAS -0.01 -1.13 PAS 0.01 0.849 
  65. PAS -0.02 -1.13 PAS 0.02 1.023 PAS -0.02 -1.13 PAS 0.02 1.261
  66. PAS -0.02 -1.13 PAS 0.02 1.611 PAS -0.02 -1.13 PAS 0.02 2.12 ; 
  67. SINON ; 
  68. L_depl=PROG 0.0  PAS 0.01 0.587 ;
  69. FINSI ; 
  70. SINON ; 
  71. SI COMPLET ; 
  72. L_depl=PROG 0.0  PAS 0.01 0.60 PAS -0.01 -1.14 PAS 0.01 0.70
  73. PAS -0.01 -1.14 PAS 0.01 0.806 PAS -0.01 -1.14 PAS 0.01 0.926 
  74. PAS -0.02 -1.14 PAS 0.02 1.067 PAS -0.02 -1.14 PAS 0.02 1.224
  75. PAS -0.02 -1.14 PAS 0.02 1.402 PAS -0.02 -1.14 PAS 0.02 1.613 
  76. PAS -0.02 -1.15 PAS 0.02 1.857 PAS -0.02 -1.16 PAS 0.02 2.141 ; 
  77. SINON ; 
  78. L_depl=PROG 0.0  PAS 0.01 0.60 ;
  79. FINSI ; 
  80. FINSI ; 
  81.  
  82. DIDI = DIME l_DEPL ; 
  83. L_temps = PROG 0.0 PAS 1. (DIDI - 1.) ;
  84. L_TEMPI  = PROG (DIME L_TEMPS) * 1.0  ;
  85. netat = DIME L_temps ;
  86. * 1.4.1 Blocages mécaniques et deplacement axial impose
  87. uz3 = 1.e-3;
  88. *
  89. clim2 = BLOQ P1 'DEPL' ET (BLOQ P2  'UZ')
  90. ET (BLOQ P4  'UY');
  91. clim3 = BLOQ shaut 'UX' ;
  92. clim1 = BLOQ sbas 'UX' ;
  93. depi3 = DEPI clim3 uz3 ;
  94. charmeca = CHAR 'DIMP' (EVOL MANU L_temps L_depl) depi3 ;
  95. vec1 = 'VECT' depi3 1.E3 'FLX' 'FLY' 'FLZ' 'ROUG' ;
  96. titr ' chargement ' ;
  97. trace vec1 vol1 ;
  98. cl0 = clim1 et clim2 et clim3 ;
  99. * Definition du modele et affectation des proprietes de materiau--------
  100. SI TENTEN  ; 
  101. mo_util = MODE vol1 'MECANIQUE' 'ELASTIQUE' ORTHOTROPE 
  102. ENDOMMAGEMENT SICTENS ;
  103. SINON ; 
  104. mo_util = MODE vol1 'MECANIQUE' 'ELASTIQUE' ORTHOTROPE 
  105. ENDOMMAGEMENT SICSCAL ;
  106. FINSI ; 
  107. **** MATERIAU : YOUNG ET ALPHA FONCTION DE T ********************** 
  108. ***** X,Y = directions d'orthotropie dans le plan; Z= epaisseur
  109. SICYGX = 230.E9 ;
  110. SICYGY = 230.E9 ;
  111. SICYGZ = 80.E9 ;
  112. SICNUXY = 0.173 ;
  113. SICNUXZ = 0.26 ;
  114. SICNUYZ = 0.26 ;
  115. SICGXY = 85.E9 ;
  116. SICGXZ = 40.E9 ;
  117. SICGYZ = 40.E9 ;
  118. SICALPX=  5.E-6 ;
  119. SICALPY = 5.E-6;
  120. SICALPZ = 5.0E-6 ;
  121. *
  122. COCO = COS TETE ; 
  123. SISI = SIN TETE ;    
  124. *
  125. SI TENTEN ; 
  126. VH1N = 1.0 ;VH1HP = 0.7 ;VH1P = 0.45 ;
  127. VH2N = 1.0 ;VH2HP = 0.7 ;VH2P = 0.45 ;
  128. VH3N = 1.0 ;VH3P = 0.7 ;
  129. VH4N = 1.0 ;VH4HP = 0.7 ;VH4P = 1.2 ;
  130. VY01 = 173.2;VY02=  173.2;VY03 = 173.2;
  131. VY04 = 100.0 ; VY05 = 100.0 ; 
  132. VYC1 = 1870.8;VYC2=  1870.8; VYC3 = 1870.8;
  133. VYC4 = 3464.1 ; VYC5 = 3464.1; 
  134. VDC1 = 4. ; VDC2 = 4. ; VDC3 = 4. ;VDC4 = 4. ; VDC5 = 4. ; 
  135. VPY1 = 1. ;VPY2 = 1. ;VPY3 = 1. ;VPY4 = 1.2 ;VPY5 = 1.2 ; 
  136. VB1 = 1. ; VB2 = 1. ; VB3 = 1. ; 
  137. VDTAL = 0. ; 
  138. VTER0 = 0. ; 
  139. VSIF1 = 0.0003 ; VSIF2 = 0.0003 ; VSIF3 = 0.0003 ; 
  140. VSIF4 = 0.0005 ; VSIF5 = 0.0005 ; 
  141. VAIF1 = 0.5 ; VAIF2 = 0.5 ; VAIF3 = 0.5 ;
  142. VAIF4 = 1. ; VAIF5 = 1. ;
  143. VETA1 = 0.1 ; VETA2 = 0.1 ; VETA3 = 0.0 ; 
  144. VETA4 = 0.1 ; VETA5 = 0.1 ;            
  145. ma_util1  = MATE mo_util DIRE (COCO SISI 0.)((SISI)(-1.*COCO) 0.)
  146. YG1 SICYGX YG2 SICYGY YG3 SICYGZ NU12 SICNUXY NU13 SICNUXZ 
  147. NU23 SICNUYZ G12  SICGXY  G13  SICGXZ  G23 SICGYZ ALP1 SICALPX ALP2 
  148. SICALPY ALP3 SICALPZ 'H1N' VH1N 'H1HP' VH1HP 'H1P' VH1P 'H2N' VH2N       
  149. 'H2HP' VH2HP 'H2P' VH2P 'H3N' VH3N 'H3P' VH3P  'H4N' VH4N 'H4HP' 
  150. VH4HP 'H4P' VH4P 'Y01' VY01 'Y02' VY02 'Y03' VY03 'Y04' VY04 
  151. 'Y05' VY05 'YC1' VYC1 'YC2' VYC2 'YC3' VYC3 'YC4' VYC4 'YC5' VYC5; 
  152. ma_util2  = MATE mo_util  DIRE (COCO SISI 0.)((SISI)(-1.*COCO) 0.)
  153. 'DC1' VDC1 'DC2' VDC2 'DC3' VDC3 'DC4' VDC4 'DC5' VDC5 
  154.  'PY1' VPY1 'PY2' VPY2 'PY3' VPY3 'PY4' VPY4 'PY5' VPY5 
  155. 'B1' VB1 'B2' VB2 'B3' VB3  'DTAL' VDTAL 'TER0' VTER0
  156. 'SIF1' VSIF1 'SIF2' VSIF2 'SIF3' VSIF3 'SIF4' VSIF4 'SIF5' 
  157. VSIF5 'AIF1' VAIF1 'AIF2' VAIF2 'AIF3' VAIF3 'AIF4' VAIF4 'AIF5' 
  158. VAIF5 ;
  159. ma_util3 = MATE mo_util  DIRE (COCO SISI 0.)((SISI)(-1.*COCO) 0.)
  160. 'ETA1' VETA1 'ETA2' VETA2 'ETA3' VETA3 'ETA4' VETA4 
  161. 'ETA5' VETA5 ; 
  162. ma_util = ma_util1 et ma_util2 et ma_util3;
  163. *
  164. * Essais sans les paramètres 
  165. ma_util  = MATE mo_util DIRE (COCO SISI 0.)((SISI)(-1.*COCO) 0.)
  166. YG1 SICYGX YG2 SICYGY YG3 SICYGZ NU12 SICNUXY NU13 SICNUXZ 
  167. NU23 SICNUYZ G12  SICGXY  G13  SICGXZ  G23 SICGYZ ALP1 SICALPX ALP2 
  168. SICALPY ALP3 SICALPZ ;
  169.  
  170. SINON ; 
  171. VH1N = 1. ;
  172. VH1HP = 0.7 ; 
  173. VH1P = 0.45 ;
  174. VH2N = 1. ;
  175. VH2HP = 0.7 ; 
  176. VH2P = 0.45 ;
  177. VH3N = 1. ; 
  178. VH3P = 0.7 ;
  179. VY0N1 =173.205 ; VY0N2 =173.205 ; VY0N3 =173.205 ; 
  180. VYCN1 =1870.83 ; VYCN2 =1870.83; VYCN3 =1870.83 ; 
  181. VY0T1 =31.623 ; VY0T2 =31.623 ;VY0T3 =173.205 ;
  182. VYCT1 =1870.83 ; VYCT2 =1870.83; VYCT3 =1870.83 ;           
  183. VDCT1 = 4.0 ; VDCT2 = 4.0 ; VDCT3 = 4.0 ; 
  184. VDCN1 = 4.0 ; VDCN2 = 4.0 ; VDCN3 = 4.0 ;  
  185. VPN1 =1. ; VPN2 = 1. ; VPN3 = 1. ;
  186. VPT1 =1.2 ; VPT2 = 1.2 ; VPT3 = 1. ;     
  187. VB =1. ;  VDTAL = 0. ; VTER0 =0. ; 
  188. VSIF1 = 0.0003  ;   VSIF2 = 0.0003  ;   VSIF3 = 0.0003  ;  
  189. VAIF1 = 0.5 ; VAIF2 = 0.5 ;VAIF3 = 0.5 ;  
  190. VETA1 =0.1 ; VETA2 = 0.1 ; VETA3 = 0. ;    
  191. **********************************************************************
  192. ma_util1  = MATE mo_util DIRE (COCO SISI 0.)((SISI)(-1.*COCO) 0.)
  193. YG1 SICYGX YG2 SICYGY YG3 SICYGZ NU12 SICNUXY NU13 SICNUXZ 
  194. NU23 SICNUYZ G12  SICGXY  G13  SICGXZ  G23 SICGYZ ALP1 SICALPX ALP2 
  195. SICALPY ALP3 SICALPZ 'H1N' VH1N 'H1HP' VH1HP 'H1P' VH1P 'H2N' VH2N       
  196. 'H2HP' VH2HP 'H2P' VH2P 'H3N' VH3N 'H3P' VH3P  'Y0N1' VY0N1 
  197. 'Y0N2' VY0N2 'Y0N3' VY0N3 'YCN1' VYCN1 'YCN2' VYCN2 'YCN3' VYCN3 ; 
  198. ma_util2  = MATE mo_util  DIRE (COCO SISI 0.)((SISI)(-1.*COCO) 0.)
  199. 'Y0T1' VY0T1 'Y0T2' VY0T2 'Y0T3' VY0T3 'YCT1' VYCT1 'YCT2' VYCT2 
  200. 'YCT3' VYCT3 'DCT1' VDCT1 'DCT2' VDCT2 'DCT3' VDCT3
  201. 'DCN1' VDCN1 'DCN2' VDCN2 'DCN3' VDCN3 PN1 VPN1 PN2 VPN2 PN3 VPN3 
  202. PT1 VPT1 PT2 VPT2 PT3 VPT3 B VB 'DTAL' VDTAL 'TER0' VTER0
  203. 'SIF1' VSIF1 'SIF2' VSIF2 'SIF3' VSIF3 'AIF1' VAIF1 'AIF2' VAIF2 
  204. 'AIF3' VAIF3 'ETA1' VETA1 'ETA2' VETA2  'ETA3' VETA3 ; 
  205. ma_util = ma_util1 et ma_util2 ;
  206. * Essais sans les paramètres 
  207. ma_util  = MATE mo_util DIRE (COCO SISI 0.)((SISI)(-1.*COCO) 0.)
  208. YG1 SICYGX YG2 SICYGY YG3 SICYGZ NU12 SICNUXY NU13 SICNUXZ 
  209. NU23 SICNUYZ G12  SICGXY  G13  SICGXZ  G23 SICGYZ ALP1 SICALPX ALP2 
  210. SICALPY ALP3 SICALPZ ;
  211. FINSI ;  
  212. * 
  213. * Resolution par PASAPAS------------------------------------------------
  214. *
  215. TABU = TABLE ;
  216. TABU.'MODELE' = mo_util ;
  217. TABU.'CARACTERISTIQUES' = ma_util ;
  218. TT = MANU 'CHPO' VOL1 1 'T'  20. ;
  219. TTI =CHAR 'TIMP' (EVOL MANU L_temps L_TEMPI) TT ;
  220. TABU.'BLOCAGES_MECANIQUES' = cl0 ;
  221. TABU.'CHARGEMENT' = charmeca ET TTI;
  222. TABU.'TEMPS_CALCULES' = L_TEMPS ;
  223. PASAPAS TABU ;
  224. mama = (dime (TABU. 'VARIABLES_INTERNES') ) - 1 ;
  225. * 
  226. * Test sur la valeur finale de SMXX 
  227. TECON = CHAN GRAVITE (TABU .CONTRAINTES . MAMA)
  228. (TABU . MODELE ) ; 
  229. VALTEST = MAXI (exco TECON 'SMXX') ;
  230. TOLERE = 1.E6 ; 
  231. SI COMPLET ; 
  232.   SI TENTEN  ;
  233.      VALVIZ = 190.E6 ; TOLERE = 4.E6 ; 
  234.     SINON ; 
  235.      VALVIZ = 210.E6 ;
  236.   FINSI ;  
  237. SINON ; 
  238.  VALVIZ = 120.E6 ; 
  239. FINSI ; 
  240. BORNI = VALVIZ - TOLERE ; 
  241. BORNS = VALVIZ + TOLERE ; 
  242. SI ((VALTEST < BORNI) OU (VALTEST > BORNS) ) ;
  243. ERRE 5 ; 
  244. FINSI ; 
  245. SI (NEG GRAPH 'N') ; 
  246. SI (NON TENTEN) ; 
  247. Esma1 = prog 0. ; 
  248. Esma2 = prog 0. ; 
  249. Esma3 = prog 0. ; 
  250. Esma4 = prog 0. ;
  251. Esma5 = prog 0. ; 
  252. Esma6 = prog 0. ; 
  253. eEMA1 =  prog 0. ; 
  254. *
  255. OPTI ECHO 1 ; 
  256. evt = prog 0. et L_temps ; 
  257. REPETER BOUBOU (mama + 1); 
  258. i = &BOUBOU - 1 ; 
  259. TADE = CHAN GRAVITE (TABU .DEFORMATIONS . I)
  260. (TABU . MODELE ) ; 
  261. VAZY = CHAN GRAVITE (TABU .VARIABLES_INTERNES . I)
  262. (TABU . MODELE ) ; 
  263. eEMA1 = eEMA1 et (prog (maxi
  264. (exco  TADE 'EPXX'))) ;
  265. esma1 = esma1 et (prog (maxi
  266. (exco VAZY 'DOM1'))) ;
  267. esma2 = esma2 et (prog (maxi
  268. (exco VAZY 'DOM2'))) ;
  269. esma3 = esma3 et (prog (maxi
  270. (exco VAZY 'DOM3'))) ;
  271. ESMA4 = ESMA4 ET (PROG (MAXI
  272. (EXCO VAZY 'NUA1'))) ; 
  273. esma5 = esma5 et (prog (maxi
  274. (exco VAZY 'NUA2'))) ;
  275. esma6 = esma6 et (prog (maxi
  276. (exco VAZY 'NUA3'))) ;
  277. fin boubou ; 
  278. evma1 = evol manu 'EPXX' EEMA1 'DOM1' ESMA1 'COUL' BLEU;
  279. evma2 = evol manu 'EPXX' EEMA1  'DOM2' ESMA2 'COUL' VERT;
  280. evma3 = evol manu 'EPXX' EEMA1  'DOM3' ESMA3 'COUL' ROSE;
  281. evma4 = evol manu 'EPXX' EEMA1  'NUA1' ESMA4 'COUL' JAUN;
  282. evma5 = evol manu 'EPXX' EEMA1  'NUA2' ESMA5 'COUL' ROUG;
  283. evma6 = evol manu 'EPXX' EEMA1 'NUA3' ESMA6 'COUL' TURQ;
  284. DESS (EVMA1 ET EVMA2 ET EVMA3 )
  285.  'TITR' 'dommage 1:bleu 2:vert 3:rose ' ;
  286. DESS (EVMA4 ET EVMA5 ET EVMA6 )
  287.  'TITR' 'COEFF NU 1:JAUN 2:ROUGE 3:TURQ ' ;
  288. *
  289. SINON ; 
  290. * pARAMÈTReS D'ENDOMMAGEMENT 
  291. EdOM1 = prog 0. ; 
  292. EDOM2 = prog 0. ; 
  293. EdOM3= prog 0. ; 
  294. EdOM4 = prog 0. ; 
  295. EdOM5 = prog 0. ; 
  296. EEMA1 = PROG 0. ;
  297. EsNU1 = prog 0. ;
  298. EsNU2 = prog 0. ; 
  299. EsNU3 = prog 0. ;
  300. EsNU4 = prog 0. ;
  301. EsNU5 = prog 0. ;  
  302. evt = prog 0. et L_temps ; 
  303. REPETER BOUBOU (mama + 1); 
  304. i = &BOUBOU - 1 ;
  305. TADE = CHAN GRAVITE (TABU .DEFORMATIONS . I)
  306. (TABU . MODELE ) ; 
  307. VAZY = CHAN GRAVITE (TABU .VARIABLES_INTERNES . I)
  308. (TABU . MODELE ) ; 
  309. eEMA1 = eEMA1 et (prog (maxi
  310. (exco TADE 'EPXX'))) ;
  311.  
  312. EdOm1 = EdOm1 et (prog (maxi
  313. (exco VAZY 'DOM1'))) ;
  314. EdOm2 = EDOM2 et (prog (maxi
  315. (exco VAZY 'DOM2'))) ;
  316. EDOM3 = eDOM3 et (prog (maxi
  317. (exco VAZY 'DOM3'))) ;
  318. EDOM4 = EDOM4 ET (PROG (MAXI
  319. (EXCO VAZY 'DOM4'))) ; 
  320. EDOM5 = EDOM5 et (prog (maxi
  321. (exco VAZY 'DOM5'))) ;
  322. ESNU1 = ESNU1 et (prog (maxi
  323. (exco VAZY 'NUA1'))) ;
  324. ESNU2 = ESNU2 et (prog (maxi
  325. (exco VAZY 'NUA2'))) ;
  326. ESNU3 = ESNU3 et (prog (maxi
  327. (exco VAZY 'NUA3'))) ;
  328. ESNU4 = ESNU4 et (prog (maxi
  329. (exco VAZY 'NUA4'))) ;
  330. ESNU5 = ESNU5 et (prog (maxi
  331. (exco VAZY 'NUA5'))) ;
  332. fin boubou ; 
  333. EVDM1 = EVOL MANU 'EXX' EEMA1 'DOM1' EDOM1 'COUL' BLEU ;
  334. EVDM2 = EVOL MANU 'EXX' EEMA1 'DOM2' EDOM2 'COUL' VERT;
  335. EVDM3 = EVOL MANU 'EXX' EEMA1 'DOM3' EDOM3 'COUL' ROSE ;
  336. EVDM4 = EVOL MANU 'EXX' EEMA1 'DOM4' EDOM4 'COUL' JAUN ;
  337. EVDM5 = EVOL MANU 'EXX' EEMA1 'DOM5' EDOM5 'COUL' ROUG ;
  338. EVNU1 = EVOL MANU 'EXX' EEMA1 'NUA1' ESNU1 'COUL' BLEU ;
  339. EVNU2 = EVOL MANU 'EXX' EEMA1 'NUA2' ESNU2 'COUL' VERT ;
  340. EVNU3 = EVOL MANU 'EXX' EEMA1 'NUA3' ESNU3 'COUL' ROSE ;
  341. EVNU4 = EVOL MANU 'EXX' EEMA1 'NUA4' ESNU4 'COUL' JAUN ;
  342. EVNU5 = EVOL MANU 'EXX' EEMA1 'NUA5' ESNU5 'COUL' ROUG ;
  343.  
  344. DESS (EVDM1 ET EVDM2 ET EVDM3 ET EVDM4 ET EVDM5)
  345. TITR 'DOM1 : B 2: V 3:ROs 4 : JAUN 5: R ' ; 
  346. *CHA11=CHAINE THEURI'T'  TOTO ; 
  347. DESS (EVNU1 ET EVNU2 ET EVNU3 ET EVNU4 ET EVNU5)
  348. TITR 'NU1 : B 2: V 3:ROs 4 : JAUN 5: R ' ; 
  349. FINSI ; 
  350. EsFI1 = prog 0. ; 
  351. EsFI2 = prog 0. ; 
  352. EsFI3 = prog 0. ; 
  353. EsFI4 = prog 0. ; 
  354. EsFI5 = prog 0. ; 
  355. EsFI6 = prog 0. ; 
  356. OPTI ECHO 1 ; 
  357. evt = prog 0. et L_temps ; 
  358. REPETER BOUBOU (mama + 1); 
  359. i = &BOUBOU - 1 ; 
  360. TECON = CHAN GRAVITE (TABU .CONTRAINTES . I)
  361. (TABU . MODELE ) ; 
  362. esFI1 = esFI1 et (prog (maxi
  363. (exco TECON 'SMXX'))) ;
  364. esFI2 = esFI2 et (prog (maxi
  365. (exco TECON 'SMYY'))) ;
  366. esFI3 = esFI3 et (prog (maxi
  367. (exco TECON 'SMZZ'))) ;
  368. ESFI4 = ESFI4 ET (PROG (MAXI
  369. (EXCO TECON 'SMXY'))) ; 
  370. esFI5 = esFI5 et (prog (maxi
  371. (exco TECON 'SMXZ'))) ;
  372. esFI6 = esFI6 et (prog (maxi
  373. (exco TECON 'SMYZ'))) ;
  374. fin boubou ; 
  375.  
  376. EEMA1 = prog 0. ; 
  377. EEMA2 = prog 0. ; 
  378. EEMA3 = prog 0. ; 
  379. EEMA4 = prog 0. ; 
  380. EEMA5 = prog 0. ; 
  381. EEMA6 = prog 0. ; 
  382. OPTI ECHO 1 ; 
  383. evt = prog 0. et L_temps ; 
  384. REPETER BOUBOU (mama + 1); 
  385. i = &BOUBOU - 1 ; 
  386. TADE = CHAN GRAVITE (TABU .DEFORMATIONS . I)
  387. (TABU . MODELE ) ; 
  388. eEMA1 = eEMA1 et (prog (maxi
  389. (exco TADE 'EPXX'))) ;
  390. eEMA2 = eEMA2 et (prog (maxi
  391. (exco TADE 'EPYY'))) ;
  392. eEMA3 = eEMA3 et (prog (maxi
  393. (exco TADE 'EPZZ'))) ;
  394. EEMA4 = EEMA4 ET (PROG (MAXI
  395. (EXCO TADE 'GAYZ'))) ; 
  396. eEMA5 = eEMA5 et (prog (maxi
  397. (exco TADE 'GAXZ'))) ;
  398. eEMA6 = eEMA6 et (prog (maxi
  399. (exco TADE 'GAXY'))) ;
  400. fin boubou ; 
  401. *
  402. EEFI1 = prog 0. ; 
  403. EEFI2 = prog 0. ; 
  404. EEFI3 = prog 0. ; 
  405. EEFI4 = prog 0. ; 
  406. EEFI5 = prog 0. ; 
  407. EEFI6 = prog 0. ; 
  408. * Courbes contraintes defo 
  409. *
  410. evFI1 = evol manu 'EXX' EEMA1 'SMXX' ESFI1  'COUL' BLEU;
  411. evFI2 = evol manu  'EYY' EEMA2  'SMYY' ESFI2  'COUL' VERT;
  412. evFI3 = evol manu 'EZZ' EEMA3  'SMZZ' ESFI3 'COUL' ROSE;
  413. eVFI4 = evol manu  'EYZ' EEMA4  'SMYZ' ESFI6  'COUL' JAUN;
  414. eVFI5 = evol manu  'EXZ' EEMA5  'SMXZ' ESFI5  'COUL' ROUG;
  415. eVFI6 = evol manu 'EXY' EEMA6  'SMXY' ESFI4 'COUL' TURQ;
  416. *
  417. *
  418. * Courbes contraintes defo 
  419. *
  420. evFI1 = evol manu 'EXX' EEMA1 'SMXX' ESFI1  'COUL' BLEU;
  421. evFI2 = evol manu  'EYY' EEMA2  'SMXX' ESFI1  'COUL' VERT;
  422. evFI3 = evol manu 'EZZ' EEMA3  'SMXX' ESFI1 'COUL' ROSE;
  423. eVFI4 = evol manu  'EYZ' EEMA4  'SMXX' ESFI1  'COUL' JAUN;
  424. eVFI5 = evol manu  'EXZ' EEMA5  'SMXX' ESFI1  'COUL' ROUG;
  425. eVFI6 = evol manu 'EXY' EEMA6  'SMXX' ESFI1 'COUL' TURQ;
  426. *
  427. DESS (EVFI1 ET EVFI2 ET EVFI3 ET EVFI4 ET EVFI5 ET EVFI6)
  428.  'TITR' 'Courbe ONERA 1:B 2:V 3:Ros 4:J 5:RG 6:Tur' ; 
  429. *
  430. EsIN1 = prog 0. ; 
  431. EsIN2 = prog 0. ; 
  432. EsIN3 = prog 0. ; 
  433. EsIN4 = prog 0. ;
  434. EsIN5 = prog 0. ; 
  435. EsIN6 = prog 0. ; 
  436. OPTI ECHO 1 ; 
  437. evt = prog 0. et L_temps ; 
  438. REPETER BOUBOU (mama + 1); 
  439. i = &BOUBOU - 1 ; 
  440. VAZY = CHAN GRAVITE (TABU .VARIABLES_INTERNES . I)
  441. (TABU . MODELE ) ; 
  442. esIN1 = esIN1 et (prog (maxi
  443. (exco VAZY 'EPI1'))) ;
  444. esIN2 = esIN2 et (prog (maxi
  445. (exco VAZY 'EPI2'))) ;
  446. esIN3 = esIN3 et (prog (maxi
  447. (exco VAZY 'EPI3'))) ;
  448. ESIN4 = ESIN4 ET (PROG (MAXI
  449. (EXCO VAZY 'EPI4'))) ; 
  450. esIN5 = esIN5 et (prog (maxi
  451. (exco VAZY 'EPI5'))) ;
  452. esIN6 = esIN6 et (prog (maxi
  453. (exco VAZY 'EPI6'))) ;
  454. fin boubou ; 
  455. evIN1 = evol manu 'EXX' EEMA1 'EPR1' ESIN1 'COUL' BLEU;
  456. evIN2 = evol manu 'EXX' EEMA1 'EPR2' ESIN2 'COUL' VERT;
  457. evIN3 = evol manu 'EXX' EEMA1 'EPR3' ESIN3 'COUL' ROSE;
  458. evIN4 = evol manu 'EXX' EEMA1'EPR4' ESIN4 'COUL' JAUN;
  459. evIN5 = evol manu 'EXX' EEMA1 'EPR5' ESIN5 'COUL' ROUG;
  460. evIN6 = evol manu 'EXX' EEMA1 'EPR6' ESIN6 'COUL' TURQ;
  461. DESS (EVIN1 ET EVIN2 ET EVIN3 ET EVIN4 ET EVIN5 ET EVIN6)
  462.  'TITR' 'DEFO INEL 1:B 2:V 3:Ros 4:J 5:RG 6:Tur' ; 
  463. *
  464. FINSI ; 
  465. FIN BOUMO ; 
  466. FIN ;  
  467.  
  468.  
  469.  
  470.  
  471.  
  472.  
  473.  
  474.  
  475.  

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