Salome HOME
CCAR: gestion des unites d'INCLUDE memorisees
[tools/eficas.git] / Exemples / ex22 / titi.comm
1
2
3 DEBUT();
4 % LECTURE DES PARAMETRES 
5 %
6 % LECTURE DES TEMPS 
7 %CYCLE 1
8 % LECTURE DU MAILLAGE ET DEFINITIONS DES PROPRIETES ELEMENTS FINIS
9
10
11 PRE_GIBI();
12 %
13
14
15 MA1=LIRE_MAILLAGE();
16 %
17
18
19 &MA1=DEFI_GROUP(CREA_GROUP_NO:(GROUP_MA:(FIXATI,APPUI,LIAS_1,LIAS_2,
20                                           EBGUIDE,MAINTIEN,
21                                           MAIN_SUP,EMBO_SUP,
22                                           TUB_GUID,GRILLES,
23                                           EMBO_INF,PIEDS,
24                                           CRAYONS),)
25                 MAILLAGE:MA1);
26 %
27 % AFFECTATION DU MODELE
28 %
29
30
31 MO1=AFFE_MODELE(AFFE:(GROUP_MA:CRAYONS
32                       PHENOMENE:'MECANIQUE'
33                       MODELISATION:'POU_D_E')
34                 AFFE:(GROUP_MA:TUB_GUID
35                       PHENOMENE:'MECANIQUE'
36                       MODELISATION:'POU_D_E')
37                 AFFE:(GROUP_MA:(EMBO_SUP,EMBO_INF)
38                       PHENOMENE:'MECANIQUE'
39                       MODELISATION:'DKT')
40                 AFFE:(GROUP_MA:PIEDS
41                       PHENOMENE:'MECANIQUE'
42                       MODELISATION:'POU_D_E')
43                 AFFE:(GROUP_MA:GRILLES
44                       PHENOMENE:'MECANIQUE'
45                       MODELISATION:'POU_D_E')
46                 AFFE:(GROUP_MA:LI_GR_CR
47                       PHENOMENE:'MECANIQUE'
48                       MODELISATION:'DIS_TR')
49                 AFFE:(GROUP_MA:LI_GR_GU
50                       PHENOMENE:'MECANIQUE'
51                       MODELISATION:'DIS_TR')
52                 AFFE:(GROUP_MA:MAINTIEN
53                       PHENOMENE:'MECANIQUE'
54                       MODELISATION:'POU_D_E')
55                 MAILLAGE:MA1);
56 %
57 % AFFECTATION DES CARACTERISTIQUES GEOMETRIQUES ELEMENTAIRES
58 %
59
60
61 CARA1=AFFE_CARA_ELEM(MODELE:MO1
62                      COQUE:(ANGL_REP:90.0
63                             GROUP_MA:EBOINF
64                             EPAIS:2.0)
65                      COQUE:(ANGL_REP:90.0
66                             GROUP_MA:EBOSUP
67                             EPAIS:3.0)
68                      DISCRET:(VALE:(1.0,2.0,3.0,4.0,5.0,6.0)
69                               REPERE:'LOCAL'
70                               CARA:'K_TR_D_L'
71                               GROUP_MA:ELASTI)
72                      DISCRET:(VALE:(1.0,2.0,3.0,4.0,5.0,6.0)
73                               REPERE:'LOCAL'
74                               CARA:'K_TR_D_L'
75                               GROUP_MA:RIGID)
76                      ORIENTATION:(VALE:(1.0,0.0,0.0)
77                                   CARA:'VECT_Y'
78                                   GROUP_MA:(GRIL_I,GRIL_E))
79                      ORIENTATION:(VALE:(0.0,1.0,0.0)
80                                   CARA:'VECT_Y'
81                                   GROUP_MA:(MAINTI_Y))
82                      ORIENTATION:(VALE:(0.0,0.0,1.0)
83                                   CARA:'VECT_Y'
84                                   GROUP_MA:(MAINTI_Z))
85                      ORIENTATION:(VALE:(0.0,1.0,0.0)
86                                   CARA:'VECT_Y'
87                                   GROUP_MA:(PI_PLE,PI_TRO))
88                      ORIENTATION:(VALE:(1.0,0.0,0.0)
89                                   CARA:'VECT_Y'
90                                   GROUP_MA:(ELASTI,RIGID))
91                      POUTRE:(SECTION:'CERCLE'
92                              GROUP_MA:CRAY
93                              VALE:(3.0,3.0)
94                              CARA:('R','EP'))
95                      POUTRE:(SECTION:'CERCLE'
96                              GROUP_MA:GUID
97                              VALE:(3.0,3.0)
98                              CARA:('R','EP'))
99                      POUTRE:(SECTION:'CERCLE'
100                              GROUP_MA:RETRE
101                              VALE:(3.0,3.0)
102                              CARA:('R','EP'))
103                      POUTRE:(SECTION:'RECTANGLE'
104                              GROUP_MA:(GRIL_I)
105                              VALE:(3.0,3.0)
106                              CARA:('HY','HZ'))
107                      POUTRE:(SECTION:'RECTANGLE'
108                              GROUP_MA:(GRIL_E)
109                              VALE:(3.0,3.0)
110                              CARA:('HY','HZ'))
111                      POUTRE:(SECTION:'RECTANGLE'
112                              GROUP_MA:(MAINTI_Y,MAINTI_Z)
113                              VALE:(3.0,3.0)
114                              CARA:('HY','HZ'))
115                      POUTRE:(SECTION:'GENERALE'
116                              GROUP_MA:PI_PLE
117                              VALE:(1.0,2.0,3.0,4.0)
118                              CARA:('A','IZ','IY','JX'))
119                      POUTRE:(SECTION:'GENERALE'
120                              GROUP_MA:PI_TRO
121                              VALE:(1.0,2.0,3.0,4.0)
122                              CARA:('A','IZ','IY','JX')));
123 %
124 % DEFINITION DES CHARGES ET CONDITIONS LIMITES
125 % DEFINITION DES TEMPERATURES NODALES EVOLUTIVES
126 %
127
128
129 F_TP1_1=DEFI_FONCTION(NOM_RESU:'TEMP'
130                       PROL_GAUCHE:'CONSTANT'
131                       NOM_PARA:'X'
132                       VALE:(0.0,1.0,5.0,4.0)
133                       PROL_DROIT:'CONSTANT');
134
135 F_TP2_1=DEFI_FONCTION(NOM_RESU:'TEMP'
136                       PROL_GAUCHE:'CONSTANT'
137                       NOM_PARA:'X'
138                       VALE:(0.0,1.0,5.0,4.0)
139                       PROL_DROIT:'CONSTANT');
140
141 F_TP3_1=DEFI_FONCTION(NOM_RESU:'TEMP'
142                       PROL_GAUCHE:'CONSTANT'
143                       NOM_PARA:'X'
144                       VALE:(0.0,1.0,5.0,4.0)
145                       PROL_DROIT:'CONSTANT');
146
147 CHTEM1_1=AFFE_CHAM_NO(GRANDEUR:'TEMP_F'
148                       AFFE:(FONCTION:F_TP1_1
149                             GROUP_NO:
150                             
151                                        (MAINTIEN,PIEDS,
152                                        TUB_GUID,GRILLES,
153                                        CRAYONS),
154                             NOM_CMP:'TEMP')
155                       AFFE:(FONCTION:(F_TP1_1,F_TP1_1,F_TP1_1)
156                             GROUP_NO:(EMBO_SUP,EMBO_INF)
157                             NOM_CMP:('TEMP_INF','TEMP','TEMP_SUP'))
158                       MAILLAGE:MA1);
159
160 CHTEM2_1=AFFE_CHAM_NO(GRANDEUR:'TEMP_F'
161                       AFFE:(FONCTION:F_TP2_1
162                             GROUP_NO:
163                             
164                                        (MAINTIEN,PIEDS,
165                                        TUB_GUID,GRILLES,
166                                        CRAYONS),
167                             NOM_CMP:'TEMP')
168                       AFFE:(FONCTION:(F_TP2_1,F_TP2_1,F_TP2_1)
169                             GROUP_NO:(EMBO_SUP,EMBO_INF)
170                             NOM_CMP:('TEMP_INF','TEMP','TEMP_SUP'))
171                       MAILLAGE:MA1);
172
173 CHTEM3_1=AFFE_CHAM_NO(GRANDEUR:'TEMP_F'
174                       AFFE:(FONCTION:F_TP3_1
175                             GROUP_NO:(MAINTIEN,PIEDS,TUB_GUID,GRILLES)
176                             NOM_CMP:'TEMP')
177                       AFFE:(FONCTION:(F_TP3_1,F_TP3_1,F_TP3_1)
178                             GROUP_NO:(EMBO_SUP,EMBO_INF)
179                             NOM_CMP:('TEMP_INF','TEMP','TEMP_SUP'))
180                       AFFE:(FONCTION:F_TP3_1
181                             GROUP_NO:CRAYONS
182                             NOM_CMP:'TEMP')
183                       MAILLAGE:MA1);
184
185 CHTH_1=CREA_RESU(NOM_CHAM:'TEMP'
186                  CHAM_GD:(INST:0.0
187                           CHAM_NO:CHTEM1_1)
188                  CHAM_GD:(INST:1000.0
189                           CHAM_NO:CHTEM1_1)
190                  CHAM_GD:(INST:2000.0
191                           CHAM_NO:CHTEM2_1)
192                  CHAM_GD:(INST:3000.0
193                           CHAM_NO:CHTEM3_1)
194                  CHAM_GD:(INST:4000.0
195                           CHAM_NO:CHTEM3_1)
196                  CHAM_GD:(INST:5000.0
197                           CHAM_NO:CHTEM2_1)
198                  CHAM_GD:(INST:6000.0
199                           CHAM_NO:CHTEM1_1)
200                  CHAM_GD:(INST:7000.0
201                           CHAM_NO:CHTEM1_1)
202                  TYPE_RESU:'EVOL_THER');
203 % CHARGEMENT: CHARGES THERMO-MECANIQUES
204 % ----------------------------------------
205
206
207 THERM_1=AFFE_CHAR_MECA(MODELE:MO1
208                        TEMP_CALCULEE:CHTH_1);
209 % DEFINITION DES CONDITIONS LIMITES
210 %            DES CONDITIONS DE CONTACT
211 %
212 % CHARGEMENT: PESANTEUR
213 %
214
215
216 PESANT_1=AFFE_CHAR_MECA(MODELE:MO1
217                         PESANTEUR:(9.8100000000000005,-1.0,0.0,0.0));
218 %
219 % CONDITIONS AUX LIMITES
220 %
221
222
223 CL_PER_1=AFFE_CHAR_MECA(MODELE:MO1
224                         DDL_IMPO:(DZ:0.0
225                                   DX:0.0
226                                   DY:0.0
227                                   GROUP_NO:FIXATI
228                                   DRZ:0.0
229                                   DRX:0.0
230                                   DRY:0.0)
231                         DDL_IMPO:(DX:0.0
232                                   GROUP_NO:APPUI)
233                         DDL_IMPO:(DZ:0.0
234                                   DY:0.0
235                                   GROUP_NO:EBGUIDE
236                                   DRZ:0.0
237                                   DRX:0.0
238                                   DRY:0.0)
239                         DDL_IMPO:(GROUP_NO:(LIAS_1,LIAS_2)
240                                   DRX:0.0)
241                         LIAISON_GROUP:(COEF_MULT_1:1.0
242                                        DDL_1:'DX'
243                                        COEF_MULT_2:-1.0
244                                        COEF_IMPO:0.0
245                                        GROUP_NO_1:LIAS_1
246                                        GROUP_NO_2:LIAS_2
247                                        DDL_2:'DX')
248                         LIAISON_GROUP:(COEF_MULT_1:1.0
249                                        DDL_1:'DY'
250                                        COEF_MULT_2:-1.0
251                                        COEF_IMPO:0.0
252                                        GROUP_NO_1:LIAS_1
253                                        GROUP_NO_2:LIAS_2
254                                        DDL_2:'DY')
255                         LIAISON_GROUP:(COEF_MULT_1:1.0
256                                        DDL_1:'DZ'
257                                        COEF_MULT_2:-1.0
258                                        COEF_IMPO:0.0
259                                        GROUP_NO_1:LIAS_1
260                                        GROUP_NO_2:LIAS_2
261                                        DDL_2:'DZ')
262                         LIAISON_GROUP:(COEF_MULT_1:1.0
263                                        DDL_1:'DRY'
264                                        COEF_MULT_2:-1.0
265                                        COEF_IMPO:0.0
266                                        GROUP_NO_1:LIAS_1
267                                        GROUP_NO_2:LIAS_2
268                                        DDL_2:'DRY')
269                         LIAISON_GROUP:(COEF_MULT_1:1.0
270                                        DDL_1:'DRZ'
271                                        COEF_MULT_2:-1.0
272                                        COEF_IMPO:0.0
273                                        GROUP_NO_1:LIAS_1
274                                        GROUP_NO_2:LIAS_2
275                                        DDL_2:'DRZ'));
276 %
277
278
279 CL_CTC=AFFE_CHAR_MECA(INFO:1
280                       MODELE:MO1
281                       CONTACT:(GROUP_MA_1:EBOINF
282                                VECT_NORM_2:(-1.0,0.0,0.0)
283                                APPARIEMENT:'NODAL'
284                                GROUP_MA_2:MCRBAS));
285 %                            
286 % DEPLACEMENT DU SYSTEME DE MAINTIEN
287 %
288
289
290 CLDI_F1=DEFI_FONCTION(PROL_GAUCHE:'CONSTANT'
291                       NOM_PARA:'INST'
292                       VALE:
293                       (0.0,0.0,1000.0,
294                              0.10000000000000001,2000.0,
295                              0.20000000000000001,3000.0,
296                              0.20000000000000001,4000.0,
297                              0.20000000000000001,5000.0,
298                              0.20000000000000001,6000.0,
299                              0.10000000000000001,7000.0,
300                              0.10000000000000001),
301                       PROL_DROIT:'CONSTANT');
302
303 CL_DI_1=AFFE_CHAR_MECA(MODELE:MO1
304                        DDL_IMPO:(DX:-1.0
305                                  GROUP_NO:MAIN_SUP));
306 % DEFINITION DES FORCES EXTERIEURES
307 %
308 % FONCTION MULTIPLICATIVE DE LA FORCE HYDRODYNAMIQUE NORMEE
309 %
310
311
312 HYDR_F1=DEFI_FONCTION(PROL_GAUCHE:'CONSTANT'
313                       NOM_PARA:'INST'
314                       VALE:
315                       (0.0,0.0,1000.0,
316                              0.10000000000000001,2000.0,
317                              0.20000000000000001,3000.0,
318                              0.20000000000000001,4000.0,
319                              0.20000000000000001,5000.0,
320                              0.20000000000000001,6000.0,
321                              0.10000000000000001,7000.0,
322                              0.10000000000000001),
323                       PROL_DROIT:'CONSTANT');
324 %
325 % PROFIL DES FORCES HYDRODYNAMIQUES TUBES GUIDES A FROID
326 %
327
328
329 HYFRTU_1=DEFI_FONCTION(PROL_GAUCHE:'CONSTANT'
330                        NOM_PARA:'X'
331                        VALE:(0.0,1.0,5.0,2.0)
332                        PROL_DROIT:'CONSTANT');
333 %
334 % PROFIL DES FORCES HYDRODYNAMIQUES CRAYON A FROID
335 %
336
337
338 HYFRCR_1=DEFI_FONCTION(PROL_GAUCHE:'CONSTANT'
339                        NOM_PARA:'X'
340                        VALE:(0.0,1.0,5.0,2.0)
341                        PROL_DROIT:'CONSTANT');
342 %
343 % CHARGEMENT : FORCES HYDRODYNAMIQUES 
344 %
345
346
347 HYDR_1=AFFE_CHAR_MECA(FORCE_POUTRE:(GROUP_MA:(GRIL10)
348                                     FX:10.0)
349                       FORCE_POUTRE:(GROUP_MA:(GRIL1)
350                                     FX:10.0)
351                       FORCE_POUTRE:(GROUP_MA:(GRIL2)
352                                     FX:10.0)
353                       FORCE_POUTRE:(GROUP_MA:(GRIL3)
354                                     FX:10.0)
355                       FORCE_POUTRE:(GROUP_MA:(GRIL4)
356                                     FX:10.0)
357                       FORCE_POUTRE:(GROUP_MA:(GRIL5)
358                                     FX:10.0)
359                       FORCE_POUTRE:(GROUP_MA:(GRIL6)
360                                     FX:10.0)
361                       FORCE_POUTRE:(GROUP_MA:(GRIL7)
362                                     FX:10.0)
363                       FORCE_POUTRE:(GROUP_MA:(GRIL8)
364                                     FX:10.0)
365                       FORCE_POUTRE:(GROUP_MA:(GRIL9)
366                                     FX:10.0)
367                       MODELE:MO1
368                       FORCE_COQUE:(FX:10.0
369                                    GROUP_MA:EBOSUP)
370                       FORCE_COQUE:(FX:10.0
371                                    GROUP_MA:EBOINF));
372
373 TOTO_1=DEFI_FONCTION(PROL_GAUCHE:'CONSTANT'
374                      NOM_PARA:'X'
375                      VALE:(0.0,10.0,5.0,10.0)
376                      PROL_DROIT:'CONSTANT');
377 %                         
378
379
380 FOHYDR_1=AFFE_CHAR_MECA_F(FORCE_POUTRE:(GROUP_MA:(GUID)
381                                         FX:TOTO_1)
382                           FORCE_POUTRE:(GROUP_MA:CRAY
383                                         FX:TOTO_1)
384                           MODELE:MO1);
385 %
386 % PROFIL DES FORCES D'ARCHIMEDE TUBE GUIDE 
387 %
388
389
390 ARFRTU_1=DEFI_FONCTION(PROL_GAUCHE:'CONSTANT'
391                        NOM_PARA:'X'
392                        VALE:(0.0,10.0,5.0,10.0)
393                        PROL_DROIT:'CONSTANT');
394 %
395 % PROFIL DES FORCES D'ARCHIMEDE CRAYON 
396 %
397
398
399 ARFRCR_1=DEFI_FONCTION(PROL_GAUCHE:'CONSTANT'
400                        NOM_PARA:'X'
401                        VALE:(0.0,10.0,5.0,10.0)
402                        PROL_DROIT:'CONSTANT');
403 %
404 % FONCTION MULTIPLICATIVE DE LA FORCE D'ARCHIMEDE
405 %
406
407
408 ARCH_F1=DEFI_FONCTION(PROL_GAUCHE:'CONSTANT'
409                       NOM_PARA:'INST'
410                       VALE:
411                       (0.0,0.0,1000.0,
412                              0.10000000000000001,2000.0,
413                              0.20000000000000001,3000.0,
414                              0.20000000000000001,4000.0,
415                              0.20000000000000001,5000.0,
416                              0.20000000000000001,6000.0,
417                              0.10000000000000001,7000.0,
418                              0.10000000000000001),
419                       PROL_DROIT:'CONSTANT');
420 %
421 %
422 % CHARGEMENT : FORCES D'ARCHIMEDE 
423 % ----------------------------------------
424 %
425
426
427 ARCH_1=AFFE_CHAR_MECA(FORCE_POUTRE:(GROUP_MA:MAINTIEN
428                                     FX:10.0)
429                       FORCE_POUTRE:(GROUP_MA:(GRIL10)
430                                     FX:10.0)
431                       FORCE_POUTRE:(GROUP_MA:(GRIL1)
432                                     FX:10.0)
433                       FORCE_POUTRE:(GROUP_MA:(GRIL2)
434                                     FX:10.0)
435                       FORCE_POUTRE:(GROUP_MA:(GRIL3)
436                                     FX:10.0)
437                       FORCE_POUTRE:(GROUP_MA:(GRIL4)
438                                     FX:10.0)
439                       FORCE_POUTRE:(GROUP_MA:(GRIL5)
440                                     FX:10.0)
441                       FORCE_POUTRE:(GROUP_MA:(GRIL6)
442                                     FX:10.0)
443                       FORCE_POUTRE:(GROUP_MA:(GRIL7)
444                                     FX:10.0)
445                       FORCE_POUTRE:(GROUP_MA:(GRIL8)
446                                     FX:10.0)
447                       FORCE_POUTRE:(GROUP_MA:(GRIL9)
448                                     FX:10.0)
449                       MODELE:MO1
450                       FORCE_COQUE:(FX:10.0
451                                    GROUP_MA:EBOSUP)
452                       FORCE_COQUE:(FX:10.0
453                                    GROUP_MA:EBOINF));
454
455 FOARCH_1=AFFE_CHAR_MECA_F(FORCE_POUTRE:(GROUP_MA:GUID
456                                         FX:ARCH_F1)
457                           FORCE_POUTRE:(GROUP_MA:CRAY
458                                         FX:ARCH_F1)
459                           MODELE:MO1);
460 %
461
462 % FORCE DE SERRAGE DE LA LIAISON GRILLE - CRAYON 
463 %
464 %
465 %
466 %
467 %
468 % DEFINITION DE LA LISTE D'INSTANTS
469 %
470
471
472 L_FRIG=DEFI_LIST_REEL(INTERVALLE:(JUSQU_A:1000.0
473                                   NOMBRE:100)
474                       DEBUT:0.0);
475 %
476 % INRTERPRETATION DES FONCTIONS
477 %
478
479
480 F1=DEFI_FONCTION(PROL_GAUCHE:'CONSTANT'
481                  NOM_PARA:'INST'
482                  VALE:(0.0,0.0,4.3200E4,0.0)
483                  PROL_DROIT:'CONSTANT');
484
485 F2=DEFI_FONCTION(PROL_GAUCHE:'CONSTANT'
486                  NOM_PARA:'INST'
487                  VALE:(0.0,0.0,4.3200E4,0.0)
488                  PROL_DROIT:'CONSTANT');
489
490 F3=DEFI_FONCTION(PROL_GAUCHE:'CONSTANT'
491                  NOM_PARA:'INST'
492                  VALE:(0.0,0.0,4.3200E4,0.0)
493                  PROL_DROIT:'CONSTANT');
494
495 F4=DEFI_FONCTION(PROL_GAUCHE:'CONSTANT'
496                  NOM_PARA:'INST'
497                  VALE:(0.0,0.0,4.3200E4,0.0)
498                  PROL_DROIT:'CONSTANT');
499
500 F5=DEFI_FONCTION(PROL_GAUCHE:'CONSTANT'
501                  NOM_PARA:'INST'
502                  VALE:(0.0,0.0,4.3200E4,0.0)
503                  PROL_DROIT:'CONSTANT');
504
505 F6=DEFI_FONCTION(PROL_GAUCHE:'CONSTANT'
506                  NOM_PARA:'INST'
507                  VALE:(0.0,0.0,4.3200E4,0.0)
508                  PROL_DROIT:'CONSTANT');
509
510 F7=DEFI_FONCTION(PROL_GAUCHE:'CONSTANT'
511                  NOM_PARA:'INST'
512                  VALE:(0.0,0.0,4.3200E4,0.0)
513                  PROL_DROIT:'CONSTANT');
514
515 F8=DEFI_FONCTION(PROL_GAUCHE:'CONSTANT'
516                  NOM_PARA:'INST'
517                  VALE:(0.0,0.0,4.3200E4,0.0)
518                  PROL_DROIT:'CONSTANT');
519
520 F9=DEFI_FONCTION(PROL_GAUCHE:'CONSTANT'
521                  NOM_PARA:'INST'
522                  VALE:(0.0,0.0,4.3200E4,0.0)
523                  PROL_DROIT:'CONSTANT');
524
525 F10=DEFI_FONCTION(PROL_GAUCHE:'CONSTANT'
526                   NOM_PARA:'INST'
527                   VALE:(0.0,0.0,4.3200E4,0.0)
528                   PROL_DROIT:'CONSTANT');
529
530 FRIG1=CALC_FONC_INTERP(FONCTION:F1
531                        PROL_DROIT:'CONSTANT'
532                        INTERPOL:'LIN'
533                        PROL_GAUCHE:'CONSTANT'
534                        LIST_PARA:L_FRIG);
535
536 FRIG2=CALC_FONC_INTERP(FONCTION:F2
537                        PROL_DROIT:'CONSTANT'
538                        INTERPOL:'LIN'
539                        PROL_GAUCHE:'CONSTANT'
540                        LIST_PARA:L_FRIG);
541
542 FRIG3=CALC_FONC_INTERP(FONCTION:F3
543                        PROL_DROIT:'CONSTANT'
544                        INTERPOL:'LIN'
545                        PROL_GAUCHE:'CONSTANT'
546                        LIST_PARA:L_FRIG);
547
548 FRIG4=CALC_FONC_INTERP(FONCTION:F4
549                        PROL_DROIT:'CONSTANT'
550                        INTERPOL:'LIN'
551                        PROL_GAUCHE:'CONSTANT'
552                        LIST_PARA:L_FRIG);
553
554 FRIG5=CALC_FONC_INTERP(FONCTION:F5
555                        PROL_DROIT:'CONSTANT'
556                        INTERPOL:'LIN'
557                        PROL_GAUCHE:'CONSTANT'
558                        LIST_PARA:L_FRIG);
559
560 FRIG6=CALC_FONC_INTERP(FONCTION:F6
561                        PROL_DROIT:'CONSTANT'
562                        INTERPOL:'LIN'
563                        PROL_GAUCHE:'CONSTANT'
564                        LIST_PARA:L_FRIG);
565
566 FRIG7=CALC_FONC_INTERP(FONCTION:F7
567                        PROL_DROIT:'CONSTANT'
568                        INTERPOL:'LIN'
569                        PROL_GAUCHE:'CONSTANT'
570                        LIST_PARA:L_FRIG);
571
572 FRIG8=CALC_FONC_INTERP(FONCTION:F8
573                        PROL_DROIT:'CONSTANT'
574                        INTERPOL:'LIN'
575                        PROL_GAUCHE:'CONSTANT'
576                        LIST_PARA:L_FRIG);
577
578 FRIG9=CALC_FONC_INTERP(FONCTION:F9
579                        PROL_DROIT:'CONSTANT'
580                        INTERPOL:'LIN'
581                        PROL_GAUCHE:'CONSTANT'
582                        LIST_PARA:L_FRIG);
583
584 FRIG10=CALC_FONC_INTERP(FONCTION:F10
585                         PROL_DROIT:'CONSTANT'
586                         INTERPOL:'LIN'
587                         PROL_GAUCHE:'CONSTANT'
588                         LIST_PARA:L_FRIG);
589 %
590 % DEFINITION DU FLUX NEUTRONIQUE RADIAL
591
592
593 FLUX_F1=DEFI_FONCTION(PROL_GAUCHE:'CONSTANT'
594                       NOM_PARA:'INST'
595                       VALE:
596                       (0.0,0.0,4.3200E4,0.0,
597                              9.74520E5,4.53280E5,2.638800E6,
598                              2.077850E6,4.087800E6,3.512790E6,
599                              5.428440E6,4.856770E6,7.245000E6,
600                              6.676270E6,8.853480E6,8.279770E6,
601                              1.0460900E7,9.872610E6,1.2202900E7,
602                              1.1586900E7,1.3944600E7,1.3295200E7,
603                              1.4143700E7,1.3490200E7,1.5563200E7,
604                              1.4875300E7,1.7236100E7,1.6506900E7,
605                              1.7948200E7,1.7199200E7,1.8190100E7,
606                              1.7433400E7,1.8193700E7,1.7435200E7,
607                              1.8197300E7,1.7435200E7,1.8200900E7,
608                              1.7435200E7),
609                       PROL_DROIT:'CONSTANT');
610 % DEFINITION DU FLUX AXIAL POUR LE CYCLE 1
611
612
613 FLUXAX1=LIRE_FONCTION(PROL_GAUCHE:'CONSTANT'
614                       UNITE:60
615                       PROL_DROIT:'CONSTANT');
616
617 FLY_1_1=DEFI_FONCTION(PROL_GAUCHE:'CONSTANT'
618                       NOM_PARA:'Y'
619                       VALE:
620                       (-0.1008,
621                              3.8032834757835002E-07,
622                              -0.088200000000000001,
623                              3.8083689458688999E-07,
624                              -0.075600000000000001,
625                              3.8296082621083001E-07,-0.063,
626                              3.858400997151E-07,-0.0504,
627                              3.8946100427349999E-07,-0.0378,
628                              3.9259455128205E-07,-0.0252,
629                              3.9258457977207998E-07,-0.0126,
630                              3.9294106125356001E-07,0.0,
631                              3.9399554843304999E-07,0.0126,
632                              3.9293233618234E-07,0.0252,
633                              3.9256712962963002E-07,0.0378,
634                              3.9256712962963002E-07,0.0504,
635                              3.8942610398860001E-07,0.063,
636                              3.8579522792023001E-07,
637                              0.075600000000000001,
638                              3.8290972222222002E-07,
639                              0.088200000000000001,
640                              3.8077706552707E-07,0.1008,
641                              3.8026103988603998E-07),
642                       PROL_DROIT:'CONSTANT');
643
644 FLY_1_2=DEFI_FONCTION(PROL_GAUCHE:'CONSTANT'
645                       NOM_PARA:'Y'
646                       VALE:
647                       (-0.1008,
648                              4.3595566239316001E-07,
649                              -0.088200000000000001,
650                              4.3866542022792E-07,
651                              -0.075600000000000001,
652                              4.4331463675214E-07,-0.063,
653                              4.4911556267806002E-07,-0.0504,
654                              4.5750658831909E-07,-0.0378,
655                              4.6647222222221997E-07,-0.0252,
656                              4.6002938034188002E-07,-0.0126,
657                              4.6004309116809001E-07,0.0,
658                              4.6686235754985997E-07,0.0126,
659                              4.6003311965811999E-07,0.0252,
660                              4.6000943732193999E-07,0.0378,
661                              4.6643981481481001E-07,0.0504,
662                              4.5746420940171002E-07,0.063,
663                              4.4906445868945999E-07,
664                              0.075600000000000001,
665                              4.4325356125356E-07,
666                              0.088200000000000001,
667                              4.3859686609687001E-07,0.1008,
668                              4.3587838319088001E-07),
669                       PROL_DROIT:'CONSTANT');
670
671 FLY_1_3=DEFI_FONCTION(PROL_GAUCHE:'CONSTANT'
672                       NOM_PARA:'Y'
673                       VALE:
674                       (-0.1008,
675                              4.9431891025640996E-07,
676                              -0.088200000000000001,
677                              5.0005128205127998E-07,
678                              -0.075600000000000001,
679                              5.1337197293447004E-07,-0.063,
680                              5.2915313390313001E-07,-0.0504,
681                              5.3788069800570004E-07,-0.0378,
682                              5.3504629629630004E-07,-0.0252,
683                              5.3221189458688998E-07,-0.0126,
684                              5.3136431623932E-07,0.0,
685                              5.3135870726495999E-07,0.0126,
686                              5.3135309829059998E-07,0.0252,
687                              5.3218696581197002E-07,0.0378,
688                              5.3500827991452996E-07,0.0504,
689                              5.3782959401709E-07,0.063,
690                              5.2909205840456002E-07,
691                              0.075600000000000001,
692                              5.1330217236466999E-07,
693                              0.088200000000000001,
694                              4.9997150997151001E-07,0.1008,
695                              4.9423165954416005E-07),
696                       PROL_DROIT:'CONSTANT');
697
698 FLY_1_4=DEFI_FONCTION(PROL_GAUCHE:'CONSTANT'
699                       NOM_PARA:'Y'
700                       VALE:
701                       (-0.1008,
702                              5.5489957264956999E-07,
703                              -0.088200000000000001,
704                              5.6477884615384998E-07,
705                              -0.075600000000000001,
706                              5.9007781339030998E-07,-0.063,
707                              5.9726415598291003E-07,-0.0504,
708                              6.0445049857550001E-07,-0.0378,
709                              5.9840776353275996E-07,-0.0252,
710                              5.8440651709402002E-07,-0.0126,
711                              5.8263906695157003E-07,0.0,
712                              5.9153240740741004E-07,0.0126,
713                              5.8262660256409999E-07,0.0252,
714                              5.8438034188033999E-07,0.0378,
715                              5.9836663105413E-07,0.0504,
716                              6.0439316239316004E-07,0.063,
717                              5.9719497863248E-07,
718                              0.075600000000000001,
719                              5.8999679487179001E-07,
720                              0.088200000000000001,
721                              5.6469159544159996E-07,0.1008,
722                              5.5480110398860005E-07),
723                       PROL_DROIT:'CONSTANT');
724
725 FLY_1_5=DEFI_FONCTION(PROL_GAUCHE:'CONSTANT'
726                       NOM_PARA:'Y'
727                       VALE:
728                       (-0.1008,
729                              6.1267200854701004E-07,
730                              -0.088200000000000001,
731                              6.2962357549858005E-07,
732                              -0.075600000000000001,
733                              6.5655413105413E-07,-0.063,
734                              6.6164583333332999E-07,-0.0504,
735                              6.5190242165241998E-07,-0.0378,
736                              6.5212678062678005E-07,-0.0252,
737                              6.3746616809116999E-07,-0.0126,
738                              6.3547435897436004E-07,0.0,
739                              6.4509686609687001E-07,0.0126,
740                              6.3546064814814995E-07,0.0252,
741                              6.3743625356125E-07,0.0378,
742                              6.5207941595442004E-07,0.0504,
743                              6.5183885327634997E-07,0.063,
744                              6.6156980056979999E-07,
745                              0.075600000000000001,
746                              6.5646563390312997E-07,
747                              0.088200000000000001,
748                              6.2952386039885999E-07,0.1008,
749                              6.1256356837606998E-07),
750                       PROL_DROIT:'CONSTANT');
751
752 FLY_1_6=DEFI_FONCTION(PROL_GAUCHE:'CONSTANT'
753                       NOM_PARA:'Y'
754                       VALE:
755                       (-0.1008,
756                              6.5755128205127996E-07,
757                              -0.088200000000000001,
758                              6.8336253561253996E-07,
759                              -0.075600000000000001,
760                              6.9012072649573003E-07,-0.063,
761                              6.9687891737892E-07,-0.0504,
762                              6.9362571225071001E-07,-0.0378,
763                              6.9074768518518997E-07,-0.0252,
764                              6.8786965811965998E-07,-0.0126,
765                              6.8586413817664004E-07,0.0,
766                              6.8585603632478996E-07,0.0126,
767                              6.8584793447293003E-07,0.0252,
768                              6.8783600427350001E-07,0.0378,
769                              6.9069720441594996E-07,0.0504,
770                              6.9355840455840002E-07,0.063,
771                              6.9679789886040002E-07,
772                              0.075600000000000001,
773                              6.9002412749288002E-07,
774                              0.088200000000000001,
775                              6.8325035612536003E-07,0.1008,
776                              6.5743162393161998E-07),
777                       PROL_DROIT:'CONSTANT');
778
779 FLY_1_7=DEFI_FONCTION(PROL_GAUCHE:'CONSTANT'
780                       NOM_PARA:'Y'
781                       VALE:
782                       (-0.1008,
783                              6.9684526353276003E-07,
784                              -0.088200000000000001,
785                              7.1406606125356005E-07,
786                              -0.075600000000000001,
787                              7.3236502849003E-07,-0.063,
788                              7.2076442307691997E-07,-0.0504,
789                              7.1793500712250995E-07,-0.0378,
790                              7.2835648148147998E-07,-0.0252,
791                              7.1444747150996998E-07,-0.0126,
792                              7.1130893874644005E-07,0.0,
793                              7.2107727920227996E-07,0.0126,
794                              7.1129398148148005E-07,0.0252,
795                              7.1441257122506995E-07,0.0378,
796                              7.2830288461537998E-07,0.0504,
797                              7.1786396011395999E-07,0.063,
798                              7.2067841880341996E-07,
799                              0.075600000000000001,
800                              7.3226157407406997E-07,
801                              0.088200000000000001,
802                              7.1395138888889005E-07,0.1008,
803                              6.9671812678063005E-07),
804                       PROL_DROIT:'CONSTANT');
805
806 FLY_1_8=DEFI_FONCTION(PROL_GAUCHE:'CONSTANT'
807                       NOM_PARA:'Y'
808                       VALE:
809                       (-0.1008,
810                              7.3679362535612996E-07,
811                              -0.088200000000000001,
812                              7.5423130341880005E-07,
813                              -0.075600000000000001,
814                              7.7215633903134002E-07,-0.063,
815                              7.5865740740741005E-07,-0.0504,
816                              7.5547774216523998E-07,-0.0378,
817                              7.6662464387464E-07,-0.0252,
818                              7.5088461538461995E-07,-0.0126,
819                              7.4313924501423998E-07,0.0,
820                              7.486323005698E-07,0.0126,
821                              7.4312054843304996E-07,0.0252,
822                              7.5084722222222E-07,0.0378,
823                              7.6656730769230998E-07,0.0504,
824                              7.5540170940170999E-07,0.063,
825                              7.5856267806267997E-07,
826                              0.075600000000000001,
827                              7.7204415954415998E-07,
828                              0.088200000000000001,
829                              7.5410790598291005E-07,0.1008,
830                              7.3665776353275997E-07),
831                       PROL_DROIT:'CONSTANT');
832
833 FLY_1_9=DEFI_FONCTION(PROL_GAUCHE:'CONSTANT'
834                       NOM_PARA:'Y'
835                       VALE:
836                       (-0.1008,
837                              7.7370940170939995E-07,
838                              -0.088200000000000001,
839                              8.0144764957264998E-07,
840                              -0.075600000000000001,
841                              8.0374358974358997E-07,-0.063,
842                              8.0603952991452995E-07,-0.0504,
843                              8.0238372507123005E-07,-0.0378,
844                              7.9936672008546997E-07,-0.0252,
845                              7.9634971509971995E-07,-0.0126,
846                              7.8317361111111004E-07,0.0,
847                              7.8316426282051005E-07,0.0126,
848                              7.8315491452990996E-07,0.0252,
849                              7.9631107549857995E-07,0.0378,
850                              7.9930689102564003E-07,0.0504,
851                              8.0230270655270996E-07,0.063,
852                              8.0593856837607005E-07,
853                              0.075600000000000001,
854                              8.0362580128205003E-07,
855                              0.088200000000000001,
856                              8.0131303418803E-07,0.1008,
857                              7.7356356837607E-07),
858                       PROL_DROIT:'CONSTANT');
859
860 FLY_1_10=DEFI_FONCTION(PROL_GAUCHE:'CONSTANT'
861                        NOM_PARA:'Y'
862                        VALE:
863                        (-0.1008,
864                               7.9424198717948999E-07,
865                               -0.088200000000000001,
866                               8.1225925925926001E-07,
867                               -0.075600000000000001,
868                               8.3075516381766E-07,-0.063,
869                               8.1540651709402001E-07,-0.0504,
870                               8.1131071937321997E-07,-0.0378,
871                               8.2297489316238996E-07,-0.0252,
872                               8.0577403846154003E-07,-0.0126,
873                               7.9717236467236002E-07,0.0,
874                               8.0284615384615001E-07,0.0126,
875                               7.9715117521368004E-07,0.0252,
876                               8.0573165954416005E-07,0.0378,
877                               8.2291257122506995E-07,0.0504,
878                               8.1122845441594998E-07,0.063,
879                               8.1530306267805998E-07,
880                               0.075600000000000001,
881                               8.3063051994301999E-07,
882                               0.088200000000000001,
883                               8.1211841168090999E-07,0.1008,
884                               7.9408742877492998E-07),
885                        PROL_DROIT:'CONSTANT');
886
887 FLY_1_11=DEFI_FONCTION(PROL_GAUCHE:'CONSTANT'
888                        NOM_PARA:'Y'
889                        VALE:
890                        (-0.1008,
891                               8.1647845441594995E-07,
892                               -0.088200000000000001,
893                               8.3510648148147996E-07,
894                               -0.075600000000000001,
895                               8.5487873931623998E-07,-0.063,
896                               8.3969088319087998E-07,-0.0504,
897                               8.3502920227919995E-07,-0.0378,
898                               8.4651014957265E-07,-0.0252,
899                               8.2973931623932002E-07,-0.0126,
900                               8.2549269943020005E-07,0.0,
901                               8.3637909544160002E-07,0.0126,
902                               8.2547275641026002E-07,0.0252,
903                               8.2969319800569997E-07,0.0378,
904                               8.4644408831909002E-07,0.0504,
905                               8.3493945868945997E-07,0.063,
906                               8.3957995014244996E-07,
907                               0.075600000000000001,
908                               8.5474786324785998E-07,
909                               0.088200000000000001,
910                               8.3495690883191003E-07,0.1008,
911                               8.1631641737892005E-07),
912                        PROL_DROIT:'CONSTANT');
913
914 FLY_1_12=DEFI_FONCTION(PROL_GAUCHE:'CONSTANT'
915                        NOM_PARA:'Y'
916                        VALE:
917                        (-0.1008,
918                               8.3948272792023003E-07,
919                               -0.088200000000000001,
920                               8.7002546296295996E-07,
921                               -0.075600000000000001,
922                               8.7608564814814997E-07,-0.063,
923                               8.8214583333333004E-07,-0.0504,
924                               8.7589245014244996E-07,-0.0378,
925                               8.7128872863248004E-07,-0.0252,
926                               8.6668500712251001E-07,-0.0126,
927                               8.6321616809117001E-07,0.0,
928                               8.6320432692307996E-07,0.0126,
929                               8.6319248575499001E-07,0.0252,
930                               8.6663764245014005E-07,0.0378,
931                               8.7121705840456005E-07,0.0504,
932                               8.7579647435896998E-07,0.063,
933                               8.8202617521368001E-07,
934                               0.075600000000000001,
935                               8.7594667022791999E-07,
936                               0.088200000000000001,
937                               8.6986716524217003E-07,0.1008,
938                               8.3931196581196996E-07),
939                        PROL_DROIT:'CONSTANT');
940
941 FLY_1_13=DEFI_FONCTION(PROL_GAUCHE:'CONSTANT'
942                        NOM_PARA:'Y'
943                        VALE:
944                        (-0.1008,
945                               8.5508315527066002E-07,
946                               -0.088200000000000001,
947                               8.7546616809116998E-07,
948                               -0.075600000000000001,
949                               9.0944533475783005E-07,-0.063,
950                               9.1292663817663998E-07,-0.0504,
951                               8.9653721509972004E-07,-0.0378,
952                               8.9547275641025999E-07,-0.0252,
953                               8.7405644586895003E-07,-0.0126,
954                               8.7004789886040001E-07,0.0,
955                               8.8226549145299002E-07,0.0126,
956                               8.7002546296295996E-07,0.0252,
957                               8.7400534188033999E-07,0.0378,
958                               8.9539672364672004E-07,0.0504,
959                               8.9643376068376001E-07,0.063,
960                               9.1279825498575999E-07,
961                               0.075600000000000001,
962                               9.0929576210826002E-07,
963                               0.088200000000000001,
964                               8.7530288461538001E-07,0.1008,
965                               8.5490491452990995E-07),
966                        PROL_DROIT:'CONSTANT');
967
968 FLY_1_14=DEFI_FONCTION(PROL_GAUCHE:'CONSTANT'
969                        NOM_PARA:'Y'
970                        VALE:
971                        (-0.1008,
972                               8.5798237179486999E-07,
973                               -0.088200000000000001,
974                               8.6880893874644002E-07,
975                               -0.075600000000000001,
976                               9.0299750712251001E-07,-0.063,
977                               9.0948646723646997E-07,-0.0504,
978                               9.1597542735043003E-07,-0.0378,
979                               9.0500302706552995E-07,-0.0252,
980                               8.8214084757835001E-07,-0.0126,
981                               8.7779825498574995E-07,0.0,
982                               8.8993732193732E-07,0.0126,
983                               8.7777207977207998E-07,0.0252,
984                               8.8208974358973997E-07,0.0378,
985                               9.0492325498575998E-07,0.0504,
986                               9.1586574074073996E-07,0.063,
987                               9.0935434472933995E-07,
988                               0.075600000000000001,
989                               9.0284294871795E-07,
990                               0.088200000000000001,
991                               8.6864066951567002E-07,0.1008,
992                               8.5779665242165003E-07),
993                        PROL_DROIT:'CONSTANT');
994
995 FLY_1_15=DEFI_FONCTION(PROL_GAUCHE:'CONSTANT'
996                        NOM_PARA:'Y'
997                        VALE:
998                        (-0.1008,
999                               8.5773557692308004E-07,
1000                               -0.088200000000000001,
1001                               8.6180893874644003E-07,
1002                               -0.075600000000000001,
1003                               8.7862838319088002E-07,-0.063,
1004                               8.9924074074074004E-07,-0.0504,
1005                               9.0872489316239005E-07,-0.0378,
1006                               9.0174545940170996E-07,-0.0252,
1007                               8.9476602564102997E-07,-0.0126,
1008                               8.9122364672365005E-07,0.0,
1009                               8.9120993589743996E-07,0.0126,
1010                               8.9119622507122997E-07,0.0252,
1011                               8.9471118233617996E-07,0.0378,
1012                               9.0166132478632003E-07,0.0504,
1013                               9.0861146723646996E-07,0.063,
1014                               8.9910363247862999E-07,
1015                               0.075600000000000001,
1016                               8.7847257834757996E-07,
1017                               0.088200000000000001,
1018                               8.6163693019942995E-07,0.1008,
1019                               8.5754362535613004E-07),
1020                        PROL_DROIT:'CONSTANT');
1021
1022 FLY_1_16=DEFI_FONCTION(PROL_GAUCHE:'CONSTANT'
1023                        NOM_PARA:'Y'
1024                        VALE:
1025                        (-0.1008,
1026                               8.5903935185184999E-07,
1027                               -0.088200000000000001,
1028                               8.5690046296295998E-07,
1029                               -0.075600000000000001,
1030                               8.5831018518518999E-07,-0.063,
1031                               8.6167307692307999E-07,-0.0504,
1032                               8.7125320512820998E-07,-0.0378,
1033                               8.8563586182335997E-07,-0.0252,
1034                               8.7089298433048003E-07,-0.0126,
1035                               8.6842877492877004E-07,0.0,
1036                               8.7942111823362E-07,0.0126,
1037                               8.6840135327634996E-07,0.0252,
1038                               8.7083814102563997E-07,0.0378,
1039                               8.8554985754985996E-07,0.0504,
1040                               8.7113977920228004E-07,0.063,
1041                               8.6153846153846001E-07,
1042                               0.075600000000000001,
1043                               8.5815313390312996E-07,
1044                               0.088200000000000001,
1045                               8.5672471509971999E-07,0.1008,
1046                               8.5884116809116995E-07),
1047                        PROL_DROIT:'CONSTANT');
1048
1049 FLY_1_17=DEFI_FONCTION(PROL_GAUCHE:'CONSTANT'
1050                        NOM_PARA:'Y'
1051                        VALE:
1052                        (-0.1008,
1053                               8.6388675213674999E-07,
1054                               -0.088200000000000001,
1055                               8.5563782051281999E-07,
1056                               -0.075600000000000001,
1057                               8.5087019230768996E-07,-0.063,
1058                               8.4755965099714995E-07,-0.0504,
1059                               8.4761698717949003E-07,-0.0378,
1060                               8.5133636039885996E-07,-0.0252,
1061                               8.4844088319088001E-07,-0.0126,
1062                               8.4633689458689002E-07,0.0,
1063                               8.4644658119657999E-07,0.0126,
1064                               8.4630947293447004E-07,0.0252,
1065                               8.4838479344729004E-07,0.0378,
1066                               8.512516025641E-07,0.0504,
1067                               8.4750480769231E-07,0.063,
1068                               8.4742378917379002E-07,
1069                               0.075600000000000001,
1070                               8.5071064814815003E-07,
1071                               0.088200000000000001,
1072                               8.5545584045584E-07,0.1008,
1073                               8.6368108974359E-07),
1074                        PROL_DROIT:'CONSTANT');
1075
1076 FLUXRD1=DEFI_NAPPE(FONCTION:(FLY_1_1,FLY_1_2,FLY_1_3,FLY_1_4,FLY_1_5,
1077                               FLY_1_6,FLY_1_7,FLY_1_8,FLY_1_9,
1078                               FLY_1_10,FLY_1_11,FLY_1_12,
1079                               FLY_1_13,FLY_1_14,FLY_1_15,
1080                               FLY_1_16,FLY_1_17),
1081                    PARA:
1082                    (-0.1008,
1083                           -0.088200000000000001,
1084                           -0.075600000000000001,-0.063,-0.0504,
1085                           -0.0378,-0.0252,-0.0126,0.0,0.0126,
1086                           0.0252,0.0378,0.0504,0.063,
1087                           0.075600000000000001,
1088                           0.088200000000000001,0.1008),
1089                    PROL_GAUCHE:'CONSTANT'
1090                    NOM_PARA:'Z'
1091                    PROL_DROIT:'CONSTANT');
1092 %
1093 % DEFINITION DES MATERIAUX ET AFFECTATIONS
1094 %
1095 %           **** DEFINITION DES MATERIAUX SANS FLUX NEUTRONIQUE ****
1096 %
1097 %
1098 % FAMILLE ASSEMBLAGE : NUMERO_FAM
1099 %
1100 % EMBOUT SUPERIEUR
1101
1102
1103 E_ES=DEFI_FONCTION(PROL_GAUCHE:'CONSTANT'
1104                    NOM_PARA:'TEMP'
1105                    VALE:(100.0,2.0,300.0,3.0)
1106                    PROL_DROIT:'CONSTANT');
1107
1108 NU_ES=DEFI_FONCTION(PROL_GAUCHE:'CONSTANT'
1109                     NOM_PARA:'TEMP'
1110                     VALE:(100.0,2.0,300.0,3.0)
1111                     PROL_DROIT:'CONSTANT');
1112
1113 AL_ES=DEFI_FONCTION(PROL_GAUCHE:'CONSTANT'
1114                     NOM_PARA:'TEMP'
1115                     VALE:(100.0,2.0,300.0,3.0)
1116                     PROL_DROIT:'CONSTANT');
1117
1118 MT_ES=DEFI_MATERIAU(ELAS_FO:(NU:NU_ES
1119                              ALPHA:AL_ES
1120                              RHO:7800.0
1121                              E:E_ES));
1122 % EMBOUT INFERIEUR
1123
1124
1125 E_EI=DEFI_FONCTION(PROL_GAUCHE:'CONSTANT'
1126                    NOM_PARA:'TEMP'
1127                    VALE:(100.0,2.0,300.0,3.0)
1128                    PROL_DROIT:'CONSTANT');
1129
1130 NU_EI=DEFI_FONCTION(PROL_GAUCHE:'CONSTANT'
1131                     NOM_PARA:'TEMP'
1132                     VALE:(100.0,2.0,300.0,3.0)
1133                     PROL_DROIT:'CONSTANT');
1134
1135 AL_EI=DEFI_FONCTION(PROL_GAUCHE:'CONSTANT'
1136                     NOM_PARA:'TEMP'
1137                     VALE:(100.0,2.0,300.0,3.0)
1138                     PROL_DROIT:'CONSTANT');
1139
1140 MT_EI=DEFI_MATERIAU(ELAS_FO:(NU:NU_EI
1141                              ALPHA:AL_EI
1142                              RHO:12.0
1143                              E:E_EI));
1144
1145 MT_RIG=DEFI_MATERIAU(ELAS:(NU:0.0
1146                            ALPHA:0.0
1147                            RHO:0.0
1148                            E:1.0));
1149 % FONCTION COUPLE EN FONCTION DE DRZ POUR LA LIAISON GRILLE/CRAYON
1150
1151
1152 FOTRAC=DEFI_FONCTION(PROL_GAUCHE:'LINEAIRE'
1153                      NOM_PARA:'DRZ'
1154                      VALE:(1.0,1.0,0.0,0.0,0.0,0.0,0.0,0.0)
1155                      PROL_DROIT:'LINEAIRE');
1156
1157 MT_RE1=DEFI_MATERIAU(DIS_CONTACT:(EFFO_N_INIT:-35.0
1158                                   RELA_MZ:FOTRAC
1159                                   RIGI_N_FO:FRIG1
1160                                   COULOMB:1.0)
1161                      ELAS:(NU:0.0
1162                            ALPHA:0.0
1163                            RHO:0.0
1164                            E:1.0));
1165
1166 MT_RE2=DEFI_MATERIAU(DIS_CONTACT:(EFFO_N_INIT:-35.0
1167                                   RELA_MZ:FOTRAC
1168                                   RIGI_N_FO:FRIG2
1169                                   COULOMB:1.0)
1170                      ELAS:(NU:0.0
1171                            ALPHA:0.0
1172                            RHO:0.0
1173                            E:1.0));
1174
1175 MT_RE3=DEFI_MATERIAU(DIS_CONTACT:(EFFO_N_INIT:-35.0
1176                                   RELA_MZ:FOTRAC
1177                                   RIGI_N_FO:FRIG3
1178                                   COULOMB:1.0)
1179                      ELAS:(NU:0.0
1180                            ALPHA:0.0
1181                            RHO:0.0
1182                            E:1.0));
1183
1184 MT_RE4=DEFI_MATERIAU(DIS_CONTACT:(EFFO_N_INIT:-35.0
1185                                   RELA_MZ:FOTRAC
1186                                   RIGI_N_FO:FRIG4
1187                                   COULOMB:1.0)
1188                      ELAS:(NU:0.0
1189                            ALPHA:0.0
1190                            RHO:0.0
1191                            E:1.0));
1192
1193 MT_RE5=DEFI_MATERIAU(DIS_CONTACT:(EFFO_N_INIT:-35.0
1194                                   RELA_MZ:FOTRAC
1195                                   RIGI_N_FO:FRIG5
1196                                   COULOMB:1.0)
1197                      ELAS:(NU:0.0
1198                            ALPHA:0.0
1199                            RHO:0.0
1200                            E:1.0));
1201
1202 MT_RE6=DEFI_MATERIAU(DIS_CONTACT:(EFFO_N_INIT:-35.0
1203                                   RELA_MZ:FOTRAC
1204                                   RIGI_N_FO:FRIG6
1205                                   COULOMB:1.0)
1206                      ELAS:(NU:0.0
1207                            ALPHA:0.0
1208                            RHO:0.0
1209                            E:1.0));
1210
1211 MT_RE7=DEFI_MATERIAU(DIS_CONTACT:(EFFO_N_INIT:-35.0
1212                                   RELA_MZ:FOTRAC
1213                                   RIGI_N_FO:FRIG7
1214                                   COULOMB:1.0)
1215                      ELAS:(NU:0.0
1216                            ALPHA:0.0
1217                            RHO:0.0
1218                            E:1.0));
1219
1220 MT_RE8=DEFI_MATERIAU(DIS_CONTACT:(EFFO_N_INIT:-35.0
1221                                   RELA_MZ:FOTRAC
1222                                   RIGI_N_FO:FRIG8
1223                                   COULOMB:1.0)
1224                      ELAS:(NU:0.0
1225                            ALPHA:0.0
1226                            RHO:0.0
1227                            E:1.0));
1228
1229 MT_RE9=DEFI_MATERIAU(DIS_CONTACT:(EFFO_N_INIT:-35.0
1230                                   RELA_MZ:FOTRAC
1231                                   RIGI_N_FO:FRIG9
1232                                   COULOMB:1.0)
1233                      ELAS:(NU:0.0
1234                            ALPHA:0.0
1235                            RHO:0.0
1236                            E:1.0));
1237
1238 MT_RE10=DEFI_MATERIAU(DIS_CONTACT:(EFFO_N_INIT:-35.0
1239                                    RELA_MZ:FOTRAC
1240                                    RIGI_N_FO:FRIG10
1241                                    COULOMB:1.0)
1242                       ELAS:(NU:0.0
1243                             ALPHA:0.0
1244                             RHO:0.0
1245                             E:1.0));
1246 % SYSTEME DE MAINTIEN
1247
1248
1249 E_MA=DEFI_FONCTION(PROL_GAUCHE:'CONSTANT'
1250                    NOM_PARA:'TEMP'
1251                    VALE:(0.0,1.0,5.0,2.0)
1252                    PROL_DROIT:'CONSTANT');
1253
1254 NU_MA=DEFI_FONCTION(PROL_GAUCHE:'CONSTANT'
1255                     NOM_PARA:'TEMP'
1256                     VALE:(0.0,1.0,5.0,2.0)
1257                     PROL_DROIT:'CONSTANT');
1258
1259 AL_MA=DEFI_CONSTANTE(VALE:0.29999999999999999);
1260
1261 SY_MA=DEFI_FONCTION(PROL_GAUCHE:'CONSTANT'
1262                     NOM_PARA:'TEMP'
1263                     VALE:(0.0,1.0,5.0,2.0)
1264                     PROL_DROIT:'CONSTANT');
1265
1266 DS_MA=DEFI_CONSTANTE(VALE:0.29999999999999999);
1267
1268 NP_MA=DEFI_FONCTION(PROL_GAUCHE:'CONSTANT'
1269                     NOM_PARA:'TEMP'
1270                     VALE:(0.0,1.0,5.0,2.0)
1271                     PROL_DROIT:'CONSTANT');
1272
1273 MEY_MA=DEFI_FONCTION(PROL_GAUCHE:'CONSTANT'
1274                      NOM_PARA:'TEMP'
1275                      VALE:(0.0,1.0,5.0,2.0)
1276                      PROL_DROIT:'CONSTANT');
1277
1278 MPY_MA=DEFI_FONCTION(PROL_GAUCHE:'CONSTANT'
1279                      NOM_PARA:'TEMP'
1280                      VALE:(0.0,1.0,5.0,2.0)
1281                      PROL_DROIT:'CONSTANT');
1282
1283 MEZ_MA=DEFI_FONCTION(PROL_GAUCHE:'CONSTANT'
1284                      NOM_PARA:'TEMP'
1285                      VALE:(0.0,1.0,5.0,2.0)
1286                      PROL_DROIT:'CONSTANT');
1287
1288 MPZ_MA=DEFI_FONCTION(PROL_GAUCHE:'CONSTANT'
1289                      NOM_PARA:'TEMP'
1290                      VALE:(0.0,1.0,5.0,2.0)
1291                      PROL_DROIT:'CONSTANT');
1292
1293 CAY_MA=DEFI_CONSTANTE(VALE:0.29999999999999999);
1294
1295 CBY_MA=DEFI_CONSTANTE(VALE:0.29999999999999999);
1296
1297 CAZ_MA=DEFI_CONSTANTE(VALE:0.29999999999999999);
1298
1299 CBZ_MA=DEFI_CONSTANTE(VALE:0.29999999999999999);
1300
1301 MPX_MA=DEFI_FONCTION(PROL_GAUCHE:'CONSTANT'
1302                      NOM_PARA:'TEMP'
1303                      VALE:(0.0,1.0,5.0,2.0)
1304                      PROL_DROIT:'CONSTANT');
1305
1306 MT_MA=DEFI_MATERIAU(VMIS_POUTRE_FO:(MEZ:MEZ_MA
1307                                     MPX:MPX_MA
1308                                     MPZ:MPZ_MA
1309                                     CBY:CBY_MA
1310                                     NP:NP_MA
1311                                     CBZ:CBZ_MA
1312                                     MEY:MEY_MA
1313                                     CAY:CAY_MA
1314                                     CAZ:CAZ_MA
1315                                     MPY:MPY_MA)
1316                     ECRO_LINE_FO:(D_SIGM_EPSI:DS_MA
1317                                   SY:SY_MA)
1318                     ELAS_FO:(NU:NU_MA
1319                              ALPHA:AL_MA
1320                              RHO:7800.0
1321                              E:E_MA));
1322 %
1323 %           **** DEFINITION DES MATERIAUX AVEC FLUX NEUTRONIQUE ****
1324 %
1325 %
1326 % GRILLE
1327
1328
1329 E_GRE=DEFI_FONCTION(PROL_GAUCHE:'CONSTANT'
1330                     NOM_PARA:'TEMP'
1331                     VALE:(0.0,1.0,5.0,2.0)
1332                     PROL_DROIT:'CONSTANT');
1333
1334 E_GRM=DEFI_FONCTION(PROL_GAUCHE:'CONSTANT'
1335                     NOM_PARA:'TEMP'
1336                     VALE:(0.0,1.0,5.0,2.0)
1337                     PROL_DROIT:'CONSTANT');
1338
1339 NU_GRE=DEFI_FONCTION(PROL_GAUCHE:'CONSTANT'
1340                      NOM_PARA:'TEMP'
1341                      VALE:(0.0,1.0,5.0,2.0)
1342                      PROL_DROIT:'CONSTANT');
1343
1344 NU_GRM=DEFI_FONCTION(PROL_GAUCHE:'CONSTANT'
1345                      NOM_PARA:'TEMP'
1346                      VALE:(0.0,1.0,5.0,2.0)
1347                      PROL_DROIT:'CONSTANT');
1348
1349 AL_GRE=DEFI_FONCTION(PROL_GAUCHE:'CONSTANT'
1350                      NOM_PARA:'TEMP'
1351                      VALE:(0.0,1.0,5.0,2.0)
1352                      PROL_DROIT:'CONSTANT');
1353
1354 AL_GRM=DEFI_FONCTION(PROL_GAUCHE:'CONSTANT'
1355                      NOM_PARA:'TEMP'
1356                      VALE:(0.0,1.0,5.0,2.0)
1357                      PROL_DROIT:'CONSTANT');
1358 % TUBE GUIDE
1359
1360
1361 E_TU=DEFI_FONCTION(PROL_GAUCHE:'CONSTANT'
1362                    NOM_PARA:'TEMP'
1363                    VALE:(0.0,1.0,5.0,2.0)
1364                    PROL_DROIT:'CONSTANT');
1365
1366 NU_TU=DEFI_FONCTION(PROL_GAUCHE:'CONSTANT'
1367                     NOM_PARA:'TEMP'
1368                     VALE:(0.0,1.0,5.0,2.0)
1369                     PROL_DROIT:'CONSTANT');
1370
1371 AL_TU=DEFI_FONCTION(PROL_GAUCHE:'CONSTANT'
1372                     NOM_PARA:'TEMP'
1373                     VALE:(0.0,1.0,5.0,2.0)
1374                     PROL_DROIT:'CONSTANT');
1375 % CRAYON
1376
1377
1378 E_CR=DEFI_FONCTION(PROL_GAUCHE:'CONSTANT'
1379                    NOM_PARA:'TEMP'
1380                    VALE:(0.0,1.0,5.0,2.0)
1381                    PROL_DROIT:'CONSTANT');
1382
1383 NU_CR=DEFI_FONCTION(PROL_GAUCHE:'CONSTANT'
1384                     NOM_PARA:'TEMP'
1385                     VALE:(0.0,1.0,5.0,2.0)
1386                     PROL_DROIT:'CONSTANT');
1387
1388 AL_CR=DEFI_FONCTION(PROL_GAUCHE:'CONSTANT'
1389                     NOM_PARA:'TEMP'
1390                     VALE:(0.0,1.0,5.0,2.0)
1391                     PROL_DROIT:'CONSTANT');
1392 % GRANDISSEMENT DES GRILLES
1393 % GRILLES EXTREMES 1 ET 10
1394
1395
1396 MT_GRE1=DEFI_MATERIAU(GRAN_IRRA:(FONC_MULT:FLUX_F1
1397                                  FLUX_L:FLUXAX1
1398                                  FLUX_TN:FLUXRD1
1399                                  S:1.0
1400                                  A:0.0)
1401                       LEMAITRE:(N:0.0
1402                                 UN_SUR_M:0.0
1403                                 QSR_K:1.0
1404                                 UN_SUR_K:0.0)
1405                       ELAS_FO:(NU:NU_GRE
1406                                ALPHA:AL_GRE
1407                                RHO:7800.0
1408                                E:E_GRE));
1409 % GRILLES DE MELANGE       
1410
1411
1412 MT_GRM1=DEFI_MATERIAU(GRAN_IRRA:(FONC_MULT:FLUX_F1
1413                                  FLUX_L:FLUXAX1
1414                                  FLUX_TN:FLUXRD1
1415                                  S:1.0
1416                                  A:0.0)
1417                       LEMAITRE:(N:0.0
1418                                 UN_SUR_M:0.0
1419                                 QSR_K:2.0
1420                                 UN_SUR_K:0.0)
1421                       ELAS_FO:(NU:NU_GRM
1422                                ALPHA:AL_GRM
1423                                RHO:7800.0
1424                                E:E_GRM));
1425 % DEFINITION DES MATERIAUX VISCOPLASTIQUES ( TUBES ET CRAYONS )
1426
1427
1428 MT_TU1=DEFI_MATERIAU(GRAN_IRRA:(FONC_MULT:FLUX_F1
1429                                 FLUX_L:FLUXAX1
1430                                 FLUX_TN:FLUXRD1
1431                                 S:1.0
1432                                 A:0.0)
1433                      LEMAITRE:(N:0.0
1434                                UN_SUR_M:0.0
1435                                QSR_K:2.0
1436                                UN_SUR_K:0.0)
1437                      ELAS_FO:(NU:NU_TU
1438                               ALPHA:AL_TU
1439                               RHO:7800.0
1440                               E:E_TU));
1441
1442 MT_CR1=DEFI_MATERIAU(GRAN_IRRA:(FONC_MULT:FLUX_F1
1443                                 FLUX_L:FLUXAX1
1444                                 FLUX_TN:FLUXRD1
1445                                 S:1.0
1446                                 A:0.0)
1447                      LEMAITRE:(N:0.0
1448                                UN_SUR_M:0.0
1449                                QSR_K:2.0
1450                                UN_SUR_K:0.0)
1451                      ELAS_FO:(NU:NU_CR
1452                               ALPHA:AL_CR
1453                               RHO:7800.0
1454                               E:E_CR));
1455 %
1456 %           AFFECTATION DES MATERIAUX A CHAUD - CYCLE  
1457 %
1458 %
1459
1460
1461 A_MAT_1=AFFE_MATERIAU(AFFE:(TEMP_REF:100.0
1462                             MATER:MT_MA
1463                             GROUP_MA:(MAINTI_Y,MAINTI_Z))
1464                       AFFE:(TEMP_REF:100.0
1465                             MATER:MT_GRM1
1466                             GROUP_MA:GRIL_I)
1467                       AFFE:(TEMP_REF:100.0
1468                             MATER:MT_GRE1
1469                             GROUP_MA:GRIL_E)
1470                       AFFE:(TEMP_REF:100.0
1471                             MATER:MT_CR1
1472                             GROUP_MA:CRAY)
1473                       AFFE:(TEMP_REF:100.0
1474                             MATER:MT_TU1
1475                             GROUP_MA:GUID)
1476                       AFFE:(TEMP_REF:100.0
1477                             MATER:MT_EI
1478                             GROUP_MA:EBOINF)
1479                       AFFE:(TEMP_REF:100.0
1480                             MATER:MT_EI
1481                             GROUP_MA:(PI_PLE,PI_TRO))
1482                       AFFE:(TEMP_REF:100.0
1483                             MATER:MT_ES
1484                             GROUP_MA:EBOSUP)
1485                       AFFE:(TEMP_REF:100.0
1486                             MATER:MT_RIG
1487                             GROUP_MA:RIGID)
1488                       AFFE:(TEMP_REF:100.0
1489                             MATER:MT_RE1
1490                             GROUP_MA:ELA1)
1491                       AFFE:(TEMP_REF:100.0
1492                             MATER:MT_RE2
1493                             GROUP_MA:ELA2)
1494                       AFFE:(TEMP_REF:100.0
1495                             MATER:MT_RE3
1496                             GROUP_MA:ELA3)
1497                       AFFE:(TEMP_REF:100.0
1498                             MATER:MT_RE4
1499                             GROUP_MA:ELA4)
1500                       AFFE:(TEMP_REF:100.0
1501                             MATER:MT_RE5
1502                             GROUP_MA:ELA5)
1503                       AFFE:(TEMP_REF:100.0
1504                             MATER:MT_RE6
1505                             GROUP_MA:ELA6)
1506                       AFFE:(TEMP_REF:100.0
1507                             MATER:MT_RE7
1508                             GROUP_MA:ELA7)
1509                       AFFE:(TEMP_REF:100.0
1510                             MATER:MT_RE8
1511                             GROUP_MA:ELA8)
1512                       AFFE:(TEMP_REF:100.0
1513                             MATER:MT_RE9
1514                             GROUP_MA:ELA9)
1515                       AFFE:(TEMP_REF:100.0
1516                             MATER:MT_RE10
1517                             GROUP_MA:ELA10)
1518                       MAILLAGE:MA1);
1519 %
1520
1521
1522 L_INST_1=DEFI_LIST_REEL(INTERVALLE:(JUSQU_A:0.5
1523                                     NOMBRE:1)
1524                         INTERVALLE:(JUSQU_A:1000.0
1525                                     NOMBRE:1)
1526                         INTERVALLE:(JUSQU_A:2000.0
1527                                     NOMBRE:1)
1528                         INTERVALLE:(JUSQU_A:3000.0
1529                                     NOMBRE:5)
1530                         INTERVALLE:(JUSQU_A:4000.0
1531                                     NOMBRE:1)
1532                         INTERVALLE:(JUSQU_A:5000.0
1533                                     NOMBRE:1)
1534                         INTERVALLE:(JUSQU_A:6000.0
1535                                     NOMBRE:1)
1536                         INTERVALLE:(JUSQU_A:7000.0
1537                                     NOMBRE:1)
1538                         INTERVALLE:(JUSQU_A:8000.0
1539                                     NOMBRE:1)
1540                         DEBUT:0.0);
1541
1542 L_ARCH_1=DEFI_LIST_REEL(INTERVALLE:(JUSQU_A:0.5
1543                                     NOMBRE:1)
1544                         INTERVALLE:(JUSQU_A:1000.0
1545                                     NOMBRE:1)
1546                         INTERVALLE:(JUSQU_A:2000.0
1547                                     NOMBRE:1)
1548                         INTERVALLE:(JUSQU_A:3000.0
1549                                     NOMBRE:5)
1550                         INTERVALLE:(JUSQU_A:4000.0
1551                                     NOMBRE:1)
1552                         INTERVALLE:(JUSQU_A:5000.0
1553                                     NOMBRE:1)
1554                         INTERVALLE:(JUSQU_A:6000.0
1555                                     NOMBRE:1)
1556                         INTERVALLE:(JUSQU_A:7000.0
1557                                     NOMBRE:1)
1558                         INTERVALLE:(JUSQU_A:8000.0
1559                                     NOMBRE:1)
1560                         DEBUT:0.0);
1561 %
1562 % CALCUL PHASE I A V DU CYCLE 1
1563 %
1564
1565
1566 RESU=STAT_NON_LINE(CHAM_MATER:A_MAT_1
1567                    MODELE:MO1
1568                    ARCHIVAGE:(LIST_INST:L_ARCH_1)
1569                    CONVERGENCE:(ITER_GLOB_MAXI:50
1570                                 ITER_INTE_MAXI:500
1571                                 RESI_GLOB_RELA:1.E-4
1572                                 ARRET:'OUI')
1573                    COMP_INCR:(RELATION:'VMIS_POU_LINE'
1574                               GROUP_MA:MAINTIEN)
1575                    COMP_INCR:(RELATION:'ASSE_COMBU'
1576                               GROUP_MA:(CRAYONS,TUB_GUID))
1577                    COMP_INCR:(RELATION:'DIS_CONTACT'
1578                               GROUP_MA:LI_GR_CR)
1579                    COMP_INCR:(COQUE_NCOU:1
1580                               RELATION:'ELAS'
1581                               GROUP_MA:
1582                               
1583                                          (EMBO_INF,EMBO_SUP,
1584                                          PIEDS,GRILLES,
1585                                          LI_GR_GU),)
1586                    INCREMENT:(NUME_INST_FIN:20
1587                               LIST_INST:L_INST_1)
1588                    CARA_ELEM:CARA1
1589                    RECH_LINEAIRE:(ITER_LINE_MAXI:3)
1590                    NEWTON:(REAC_ITER:1
1591                            MATRICE:'TANGENTE')
1592                    SOLVEUR:(METHODE:'MULT_FRONT'
1593                             RENUM:'MD')
1594                    EXCIT:(FONC_MULT:ARCH_F1
1595                           CHARGE:ARCH_1)
1596                    EXCIT:(FONC_MULT:ARCH_F1
1597                           CHARGE:FOARCH_1)
1598                    EXCIT:(FONC_MULT:HYDR_F1
1599                           CHARGE:HYDR_1)
1600                    EXCIT:(FONC_MULT:HYDR_F1
1601                           CHARGE:FOHYDR_1)
1602                    EXCIT:(FONC_MULT:CLDI_F1
1603                           CHARGE:CL_DI_1)
1604                    EXCIT:(CHARGE:CL_PER_1)
1605                    EXCIT:(CHARGE:CL_CTC)
1606                    EXCIT:(CHARGE:PESANT_1)
1607                    EXCIT:(CHARGE:THERM_1)
1608                    PARM_THETA:0.5);
1609
1610 FIN();