Télécharger la plaquette Cast3M


fr en

Télécharger coq2mas.procedur

Retour à la liste

Numérotation des lignes : 

  1. * COQ2MAS   PROCEDUR  PASCAL    14/01/23    21:15:02     7911           
  2. *************************************************************************
  3. *
  4. *
  5. *  Génère un maillage volumique à partir d'un modèle 
  6. *            de coques multicouches
  7. *
  8. **************************************************************************
  9. DEBPROC COQ2MAS MOD1*MMODEL MAT1*MCHAML MOTCLE/MOT TAB1*TABLE; 
  10. *
  11. OPTI ELEM CUB8;
  12. *
  13. ZONMOD = extr MOD1 ZONE;
  14. *
  15. FLAGDEFO = FAUX;
  16. FLAGCONT = FAUX;
  17. FLAGVARI = FAUX;
  18. FLAGRELA = FAUX;
  19. *
  20. SI (NEG (TYPE TAB1) TABLE);
  21.   TAB1 = TABLE;
  22. FINSI;
  23. *
  24. * On regarde si la table TAB1 existe sinon on la crée
  25. * On la complète le cas échéant
  26. *
  27. SI (EXISTE TAB1);
  28. *
  29.     SI (EXISTE TAB1 'DEPLACEMENTS');
  30.       FLAGDEFO = VRAI;
  31.       NDEFO = DIME (TAB1.'DEPLACEMENTS');
  32.      SINON;
  33.       FLAGDEFO = FAUX;
  34.       NDEFO = 0;
  35.     FINSI;
  36. *
  37.     SI (EXISTE TAB1 'CONTRAINTES');
  38.       FLAGCONT = VRAI;
  39.       NCONT = DIME (TAB1.'CONTRAINTES');
  40.      SINON;
  41.       FLAGCONT = FAUX;
  42.       NCONT = 0;
  43.     FINSI;
  44. *
  45. *
  46.     SI (EXISTE TAB1 'VARIABLES_INTERNES');
  47.       FLAGVARI = VRAI;
  48.       NVARI = DIME (TAB1.'VARIABLES_INTERNES');
  49.      SINON;
  50.       FLAGVARI = FAUX;
  51.       NVARI = 0;
  52.     FINSI;
  53. *
  54. *
  55. * Si on veut creer les relations pour les sections au droit des noeuds
  56. *
  57.     BLCOQ3D = MOT 'BLCOQ3D';
  58.     SI (EXISTE TAB1 'RELATION_3D');
  59.       SI (EGA (TAB1.'RELATION_3D') VRAI);
  60.         FLAGRELA = VRAI;
  61.       SINON; 
  62.         FLAGRELA = FAUX;      
  63.       FINSI;
  64.     SINON;
  65.       FLAGRELA = FAUX;
  66.     FINSI;
  67. *
  68.     SI (NON (EXISTE TAB1 'MODELE'));
  69.       TAB1. 'MODELE' = TABLE;
  70.     FINSI; 
  71. *
  72.     SI (NON (EXISTE TAB1 'MATERIAU'));
  73.       TAB1. 'MATERIAU' = TABLE;
  74.     FINSI; 
  75. *
  76.     SI (NON (EXISTE TAB1 'MAILLAGE_FIBRE_MOYENNE'));
  77.       TAB1. 'MAILLAGE_FIBRE_MOYENNE' = TABLE;
  78.     FINSI; 
  79. *
  80.     SI (NON (EXISTE TAB1 'MAILLAGE_FIBRE_INFERIEURE'));
  81.       TAB1. 'MAILLAGE_FIBRE_INFERIEURE' = TABLE;
  82.     FINSI; 
  83. *
  84.     SI (NON (EXISTE TAB1 'MAILLAGE_FIBRE_SUPERIEURE'));
  85.       TAB1. 'MAILLAGE_FIBRE_SUPERIEURE' = TABLE;
  86.     FINSI;    
  87. *
  88.     SI (NON (EXISTE TAB1 'MAILLAGE_VOLUMIQUE'));
  89.       TAB1. 'MAILLAGE_VOLUMIQUE' = TABLE;
  90.     FINSI; 
  91. *
  92.     SI (NON (EXISTE TAB1 'DEPLACEMENTS_FIBRE_MOYENNE'));
  93.       TAB1. 'DEPLACEMENTS_FIBRE_MOYENNE' = TABLE;
  94.     FINSI; 
  95. *
  96.     SI (NON (EXISTE TAB1 'VARI_FIBRE_MOYENNE'));
  97.       TAB1. 'VARI_FIBRE_MOYENNE' = TABLE;
  98.     FINSI;   
  99. *
  100.     SI (NON (EXISTE TAB1 'CONTRAINTES_FIBRE_MOYENNE'));
  101.       TAB1. 'CONTRAINTES_FIBRE_MOYENNE' = TABLE;
  102.     FINSI;  
  103. *
  104.     SI (NON (EXISTE TAB1 'CONTRAINTES_FIBRE_SUPERIEURE'));
  105.       TAB1. 'CONTRAINTES_FIBRE_SUPERIEURE' = TABLE;
  106.     FINSI;   
  107. *
  108.     SI (NON (EXISTE TAB1 'CONTRAINTES_FIBRE_INFERIEURE'));
  109.       TAB1. 'CONTRAINTES_FIBRE_INFERIEURE' = TABLE;
  110.     FINSI;   
  111. *
  112.     SI (NON (EXISTE TAB1 'CONTRAINTES_VOLUMIQUE'));
  113.       TAB1. 'CONTRAINTES_VOLUMIQUE' = TABLE;
  114.     FINSI;                      
  115. *
  116.     SI (NON (EXISTE TAB1 'DEPLACEMENTS_FIBRE_INFERIEURE'));
  117.       TAB1. 'DEPLACEMENTS_FIBRE_INFERIEURE' = TABLE;
  118.     FINSI; 
  119. *
  120.     SI (NON (EXISTE TAB1 'DEPLACEMENTS_FIBRE_SUPERIEURE'));
  121.       TAB1. 'DEPLACEMENTS_FIBRE_SUPERIEURE' = TABLE;
  122.     FINSI;    
  123. *
  124.     SI (NON (EXISTE TAB1 'DEPLACEMENTS_VOLUMIQUE'));
  125.       TAB1. 'DEPLACEMENTS_VOLUMIQUE' = TABLE;
  126.     FINSI;       
  127. *
  128.     SI (NON (EXISTE TAB1 'NORMALE'));
  129.       TAB1. 'NORMALE' = TABLE;
  130.     FINSI;       
  131. *
  132. SINON;
  133. *
  134.  TAB1 = TABLE;
  135.  TAB1. 'MODELE' = TABLE;
  136.  TAB1. 'MATERIAU' = TABLE;
  137.  TAB1. 'MAILLAGE_FIBRE_MOYENNE' = TABLE; 
  138.  TAB1. 'MAILLAGE_FIBRE_INFERIEURE' = TABLE; 
  139.  TAB1. 'MAILLAGE_FIBRE_SUPERIEURE' = TABLE; 
  140.  TAB1. 'MAILLAGE_VOLUMIQUE' = TABLE; 
  141.  TAB1. 'DEPLACEMENTS_FIBRE_MOYENNE' = TABLE; 
  142.  TAB1. 'DEPLACEMENTS_FIBRE_INFERIEURE' = TABLE; 
  143.  TAB1. 'DEPLACEMENTS_FIBRE_SUPERIEURE' = TABLE; 
  144.  TAB1. 'DEPLACEMENTS_VOLUMIQUE' = TABLE; 
  145.  TAB1. 'NORMALE' = TABLE;
  146.  TAB1. 'VARI_FIBRE_MOYENNE' = TABLE;
  147.  TAB1. 'CONTRAINTES_FIBRE_MOYENNE' = TABLE;
  148.  TAB1. 'CONTRAINTES_FIBRE_INFERIEURE' = TABLE;
  149.  TAB1. 'CONTRAINTES_FIBRE_SUPERIEURE' = TABLE;
  150.  TAB1. 'CONTRAINTES_VOLUMIQUE' = TABLE;
  151. *
  152. FINSI; 
  153. *
  154. SI FLAGDEFO;
  155.  REPETER LAB3 NDEFO;
  156.    i3 = (INDEX (TAB1. 'DEPLACEMENTS')) . &lab3 ; 
  157.    TAB1. 'DEPLACEMENTS_FIBRE_MOYENNE' . i3 = TABLE;
  158.    TAB1. 'DEPLACEMENTS_FIBRE_INFERIEURE' . i3 = TABLE; 
  159.    TAB1. 'DEPLACEMENTS_FIBRE_SUPERIEURE' . i3 = TABLE; 
  160.    TAB1. 'DEPLACEMENTS_VOLUMIQUE' . i3 = TABLE; 
  161.  FIN LAB3;
  162. FINSI;
  163. SI FLAGCONT;
  164.  REPETER LAB3 NCONT;
  165.    i3 = (INDEX (TAB1. 'CONTRAINTES')) . &lab3 ; 
  166.   TAB1. 'CONTRAINTES_FIBRE_MOYENNE' . i3 = TABLE;  
  167.   TAB1. 'CONTRAINTES_FIBRE_INFERIEURE' . i3 = TABLE;
  168.   TAB1. 'CONTRAINTES_FIBRE_SUPERIEURE' . i3 = TABLE;  
  169.   TAB1. 'CONTRAINTES_VOLUMIQUE' . i3 = TABLE;  
  170.  FIN LAB3;
  171. FINSI;
  172. SI FLAGVARI;
  173.  REPETER LAB3 NVARI;
  174.    i3 = (INDEX (TAB1. 'VARIABLES_INTERNES')) . &lab3 ; 
  175.    TAB1. 'VARI_FIBRE_MOYENNE' . i3 = TABLE;    
  176.  FIN LAB3;
  177. FINSI;
  178. *
  179. * Fin de la lecture de la table d'entrée
  180. *
  181. * Génération du maillage
  182. *
  183.  
  184. NBZONE = ENTIER ((DIME ZONMOD)/2.);
  185. *
  186. * On boucle sur les zones
  187. *
  188. REPETER LAB1 NBZONE;
  189. *
  190.   I1= &LAB1;
  191. *
  192. * On ne crée les maillages que s'ils n'existent pas
  193. *
  194.   SI (NON (EXISTE  TAB1.'MAILLAGE_FIBRE_MOYENNE'  I1));
  195. *
  196.    MESHCOUI =  ZONMOD.(2*I1);
  197.    MODCOUI = ZONMOD.((2*I1) - 1);
  198.    MATCOUI = REDU MAT1 MODCOUI;
  199. *
  200.    VVJ = VSUR MODCOUI NORM;
  201.    VVJP = EXCO (CHAN VVJ CHPO MODCOUI)
  202.         (mots 'VX' 'VY' 'VZ') (mots 'UX' 'UY' 'UZ') NATURE DIFFUS;
  203.    EPAP = CHAN (EXCO MATCOUI 'EPAI' SCAL) CHPO MODCOUI;        
  204.    SI (EXISTE MATCOUI 'EXCE');        
  205.      EXCP = CHAN (EXCO MATCOUI 'EXCE' SCAL) CHPO MODCOUI;
  206.    SINON;
  207.      EXCP = 0.*EPAP;
  208.    FINSI;
  209.  
  210.    MESHMOY = MESHCOUI PLUS (excp*vvjp);
  211.    MESHINF = (MESHCOUI PLUS ((excp - (0.5*epap))*vvjp)); 
  212.    MESHSUP = (MESHCOUI PLUS ((excp + (0.5*epap))*vvjp));
  213.    MESHVOL  = VOLU 1 MESHSUP MESHINF; 
  214. *           
  215.    TAB1. 'MODELE'  . I1 = MODCOUI;
  216.    TAB1. 'MATERIAU' . I1 = MATCOUI; 
  217.    TAB1. 'NORMALE' . I1 = VVJP;
  218.    TAB1. 'MAILLAGE_FIBRE_MOYENNE' . I1 = MESHMOY;
  219.    TAB1. 'MAILLAGE_FIBRE_INFERIEURE' . I1 = MESHINF;
  220.    TAB1. 'MAILLAGE_FIBRE_SUPERIEURE' . I1 = MESHSUP;
  221.    TAB1. 'MAILLAGE_VOLUMIQUE'. I1 = MESHVOL;
  222. *
  223.   SI FLAGRELA;
  224.     NNN1 = NBNO MESHCOUI;   
  225.     REPETER LAB100 NNN1;
  226.        I100 = &LAB100;
  227.        MESHINF2 = CHAN POI1 MESHINF; 
  228.        MESHSUP2 = CHAN POI1 MESHSUP; 
  229.        MESHCOU2 = CHAN POI1 MESHCOUI; 
  230.        PSUPJ = MESHINF2 POIN I100;
  231.        PINFJ = MESHSUP2 POIN I100;
  232.        PCOQJ = MESHCOU2 POIN I100;
  233.        BLCOQJ = RELA CORI DEPL ROTA PCOQJ (PSUPJ ET PINFJ);
  234.        SI (EGA (TYPE BLCOQ3D) MOT);
  235.          BLCOQ3D = BLCOQJ;
  236.        SINON;
  237.          BLCOQ3D = BLCOQ3D ET BLCOQJ;
  238.        FINSI; 
  239.     FIN LAB100;
  240.   FINSI;
  241. *
  242.    SI (I1 EGA 1);
  243.     TAB1. 'MODELE'  . 'TOTAL' = MODCOUI;
  244.     TAB1. 'MATERIAU' . 'TOTAL' = MATCOUI; 
  245. *    TAB1. 'NORMALE' . 'TOTAL' = VVJP;
  246.     TAB1. 'MAILLAGE_FIBRE_MOYENNE' . 'TOTAL' = MESHMOY;
  247.     TAB1. 'MAILLAGE_FIBRE_INFERIEURE' . 'TOTAL' = MESHINF;
  248.     TAB1. 'MAILLAGE_FIBRE_SUPERIEURE' . 'TOTAL' = MESHSUP;
  249.     TAB1. 'MAILLAGE_VOLUMIQUE'. 'TOTAL' = MESHVOL;
  250.    SINON;
  251.     TAB1. 'MODELE'  . 'TOTAL' = 
  252.      TAB1. 'MODELE'  . 'TOTAL'  et MODCOUI;
  253.     TAB1. 'MATERIAU' . 'TOTAL' = 
  254.      TAB1. 'MATERIAU' . 'TOTAL' et MATCOUI; 
  255. *    TAB1. 'NORMALE' . 'TOTAL' = 
  256. *    TAB1. 'NORMALE' . 'TOTAL' et VVJP;
  257.     TAB1. 'MAILLAGE_FIBRE_MOYENNE' . 'TOTAL' = 
  258.       TAB1. 'MAILLAGE_FIBRE_MOYENNE' . 'TOTAL'  
  259.        et MESHMOY;
  260.     TAB1. 'MAILLAGE_FIBRE_INFERIEURE' . 'TOTAL' = 
  261.      TAB1. 'MAILLAGE_FIBRE_INFERIEURE' . 'TOTAL'  
  262.      et MESHSUP;
  263.     TAB1. 'MAILLAGE_FIBRE_SUPERIEURE' . 'TOTAL' = 
  264.      TAB1. 'MAILLAGE_FIBRE_SUPERIEURE' . 'TOTAL'  
  265.      et MESHINF;
  266.     TAB1. 'MAILLAGE_VOLUMIQUE'. 'TOTAL' = 
  267.      TAB1. 'MAILLAGE_VOLUMIQUE'. 'TOTAL' et MESHVOL;
  268.    FINSI;
  269. *
  270.   FINSI;
  271. *
  272. * Transport des contraintes
  273. *
  274.   SI FLAGCONT;
  275.    MESHORI = (EXTR (TAB1 . 'MODELE' . I1) MAILLAGE) ;
  276.    MODCOUI = TAB1 . 'MODELE' . I1;
  277.    MATCOUI = TAB1 . 'MATERIAU' . I1;   
  278.    REPETER LAB21 NCONT;
  279.      i21 = (INDEX (TAB1.'CONTRAINTES')). &LAB21; 
  280.      CONJ = REDU ((TAB1.'CONTRAINTES'). i21) MODCOUI;
  281.      CONJM = CALP CONJ MODCOUI MATCOUI MOYE;
  282.      CONJS = CALP CONJ MODCOUI MATCOUI SUPE;
  283.      CONJI = CALP CONJ MODCOUI MATCOUI INFE;
  284.      CONJMB = CHAN CHPO CONJM MODCOUI;
  285.      CONJSB = CHAN CHPO CONJS MODCOUI;
  286.      CONJIB = CHAN CHPO CONJI MODCOUI;  
  287. *        
  288.      CONJM3D = DEDU CONJMB TRAN MESHORI
  289.          (TAB1. 'MAILLAGE_FIBRE_MOYENNE' . I1);
  290.      TAB1. 'CONTRAINTES_FIBRE_MOYENNE' . i21 . I1 = CONJM3D;
  291.      SI (I1 EGA 1);
  292.        TAB1. 'CONTRAINTES_FIBRE_MOYENNE' . i21 . 'TOTAL' 
  293.          = CONJM3D;
  294.      SINON;
  295.        TAB1. 'CONTRAINTES_FIBRE_MOYENNE' . i21 . 'TOTAL' = 
  296.          (TAB1. 'CONTRAINTES_FIBRE_MOYENNE' . i21 . 'TOTAL') 
  297.          ET CONJM3D;    
  298.      FINSI;  
  299.      CONJS3D = DEDU CONJSB TRAN MESHORI
  300.          (TAB1. 'MAILLAGE_FIBRE_SUPERIEURE' . I1);
  301.      TAB1. 'CONTRAINTES_FIBRE_SUPERIEURE' . i21 . I1 = CONJS3D;
  302.      SI (I1 EGA 1);
  303.        TAB1. 'CONTRAINTES_FIBRE_SUPERIEURE' . i21 . 'TOTAL' 
  304.          = CONJS3D;
  305.      SINON;
  306.        TAB1. 'CONTRAINTES_FIBRE_SUPERIEURE' . i21 . 'TOTAL' = 
  307.          (TAB1. 'CONTRAINTES_FIBRE_SUPERIEURE' . i21 . 'TOTAL') 
  308.          ET CONJS3D;    
  309.      FINSI;  
  310.      CONJI3D = DEDU CONJIB TRAN MESHORI
  311.          (TAB1. 'MAILLAGE_FIBRE_INFERIEURE' . I1);
  312.      TAB1. 'CONTRAINTES_FIBRE_INFERIEURE' . i21 . I1 = CONJi3D;
  313.      SI (I1 EGA 1);
  314.        TAB1. 'CONTRAINTES_FIBRE_INFERIEURE' . i21 . 'TOTAL' 
  315.          = CONJi3D;
  316.      SINON;
  317.        TAB1. 'CONTRAINTES_FIBRE_INFERIEURE' . i21 . 'TOTAL' = 
  318.          (TAB1. 'CONTRAINTES_FIBRE_INFERIEURE' . i21 . 'TOTAL') 
  319.          ET CONJI3D;    
  320.      FINSI;  
  321. *
  322.      TAB1. 'CONTRAINTES_VOLUMIQUE' . i21 . I1 
  323.           = CONJS3D ET CONJI3D;
  324.      SI (I1 EGA 1);
  325.        TAB1. 'CONTRAINTES_VOLUMIQUE' . i21 . 'TOTAL' 
  326.          = CONJS3D ET CONJI3D;
  327.      SINON;
  328.        TAB1. 'CONTRAINTES_VOLUMIQUE' . i21. 'TOTAL' = 
  329.          (TAB1. 'CONTRAINTES_VOLUMIQUE' . i21 . 'TOTAL') 
  330.          ET CONJS3D ET CONJI3D;    
  331.      FINSI;          
  332.    FIN LAB21;       
  333.   FINSI;
  334. *
  335. *
  336. * Transport des variables internes
  337. *
  338.   SI FLAGVARI;
  339.    MESHORI = (EXTR (TAB1 . 'MODELE' . I1) MAILLAGE) ;
  340.    MODCOUI = TAB1 . 'MODELE' . I1;
  341.    MATCOUI = TAB1 . 'MATERIAU' . I1;      
  342.    REPETER LAB22 NVARI;
  343.      i22 = (INDEX (TAB1.'VARIABLES_INTERNES')). &LAB22; 
  344.      VARJ = REDU ((TAB1.'VARIABLES_INTERNES'). i22) MODCOUI;
  345.      VARJB = CHAN CHPO VARJ MODCOUI;
  346.      VARJ3D = DEDU VARJB TRAN MESHORI
  347.          (TAB1. 'MAILLAGE_FIBRE_MOYENNE' . I1);
  348.      TAB1. 'VARI_FIBRE_MOYENNE' . i22 . I1 = VARJ3D;
  349.      SI (I1 EGA 1);
  350.        TAB1. 'VARI_FIBRE_MOYENNE' . i22 . 'TOTAL' 
  351.          = VARJ3D;
  352.      SINON;
  353.        TAB1. 'VARI_FIBRE_MOYENNE' . i22 . 'TOTAL' = 
  354.          (TAB1. 'VARI_FIBRE_MOYENNE' . i22 . 'TOTAL') 
  355.          ET VARJ3D;    
  356.      FINSI;  
  357.    FIN LAB22;       
  358.   FINSI;
  359. *
  360. * Calcul des deformees sur le 3D
  361. *
  362.   SI FLAGDEFO;
  363. *
  364.      VECOP = TAB1. 'NORMALE' . I1;
  365.      MESHORI = (EXTR (TAB1 . 'MODELE' . I1) MAILLAGE) ;
  366.      MODCOUI = TAB1 . 'MODELE' . I1;
  367.      MATCOUI = TAB1 . 'MATERIAU' . I1; 
  368.      EPAP = CHAN (EXCO MATCOUI EPAI SCAL) CHPO MODCOUI;     
  369.      SI (EXISTE MATCOUI 'EXCE');        
  370.        EXCP = CHAN (EXCO MATCOUI 'EXCE' SCAL) CHPO MODCOUI;
  371.      SINON;
  372.        EXCP = 0.*EPAP;
  373.      FINSI;       
  374.      VECOPM= (DEDU (EXCP* VECOP) TRAN MESHORI
  375.             (TAB1. 'MAILLAGE_FIBRE_MOYENNE' . I1));  
  376.      VECOPI= (DEDU ((excp - (0.5*epap))*VECOP)
  377.               TRAN MESHORI
  378.             (TAB1. 'MAILLAGE_FIBRE_INFERIEURE' . I1));  
  379.      VECOPS= (DEDU ((excp + (0.5*epap))*VECOP)
  380.               TRAN MESHORI
  381.             (TAB1. 'MAILLAGE_FIBRE_SUPERIEURE' . I1)); 
  382. *
  383.    REPETER LAB2 NDEFO;
  384.      i2 = (INDEX (TAB1.'DEPLACEMENTS')). &LAB2; 
  385.      DEPJ = REDU ((TAB1.'DEPLACEMENTS'). i2) MESHORI;
  386. *
  387.      DEPT = EXCO DEPJ (mots UX UY UZ) (mots UX UY UZ)
  388.              NATURE DIFFUS;
  389.      DEPR = EXCO DEPJ (mots RX RY RZ) (mots RX RY RZ)
  390.              NATURE DIFFUS;
  391. *
  392.      DEPTM DEPRM = DEDU DEPT DEPR TRAN MESHORI
  393.             (TAB1. 'MAILLAGE_FIBRE_MOYENNE' . I1);
  394.      DEPMOY =  DEPTM +
  395.         ((EXCO DEPRM RY SCAL)*(EXCO VECOPM UZ UX))
  396.       - ((EXCO DEPRM RZ SCAL)*(EXCO VECOPM UY UX))
  397.       + ((EXCO DEPRM RZ SCAL)*(EXCO VECOPM UX UY))
  398.       - ((EXCO DEPRM RX SCAL)*(EXCO VECOPM UZ UY));
  399.      DEPMOY =  DEPMOY +
  400.         ((EXCO DEPRM RX SCAL)*(EXCO VECOPM UY UZ))
  401.       - ((EXCO DEPRM RY SCAL)*(EXCO VECOPM UX UZ));
  402.      DEPMOY = CHAN ATTRIBUT DEPMOY NATURE   DIFFUS;
  403.      TAB1. 'DEPLACEMENTS_FIBRE_MOYENNE' . i2 . I1 = DEPMOY;
  404.      SI (I1 EGA 1);
  405.        TAB1. 'DEPLACEMENTS_FIBRE_MOYENNE' . i2 . 'TOTAL' 
  406.          = DEPMOY;
  407.      SINON;
  408.        TAB1. 'DEPLACEMENTS_FIBRE_MOYENNE' . i2. 'TOTAL' = 
  409.          (TAB1. 'DEPLACEMENTS_FIBRE_MOYENNE' . i2 . 'TOTAL') 
  410.          ET DEPMOY;    
  411.      FINSI;
  412. *
  413. *
  414.      DEPTI DEPRI = DEDU DEPT DEPR TRAN MESHORI
  415.             (TAB1. 'MAILLAGE_FIBRE_INFERIEURE' . I1);
  416.      DEPINF =  DEPTI +
  417.         ((EXCO DEPRI RY SCAL)*(EXCO VECOPI UZ UX))
  418.       - ((EXCO DEPRI RZ SCAL)*(EXCO VECOPI UY UX))
  419.       + ((EXCO DEPRI RZ SCAL)*(EXCO VECOPI UX UY))
  420.       - ((EXCO DEPRI RX SCAL)*(EXCO VECOPI UZ UY));
  421.      DEPINF =  DEPINF +
  422.         ((EXCO DEPRI RX SCAL)*(EXCO VECOPI UY UZ))
  423.       - ((EXCO DEPRI RY SCAL)*(EXCO VECOPI UX UZ));
  424.      DEPINF = CHAN ATTRIBUT DEPINF NATURE   DIFFUS;
  425.      TAB1. 'DEPLACEMENTS_FIBRE_INFERIEURE' . i2 . I1 = DEPINF;
  426.      SI (I1 EGA 1);
  427.        TAB1. 'DEPLACEMENTS_FIBRE_INFERIEURE' . i2 . 'TOTAL' 
  428.          = DEPINF;
  429.      SINON;
  430.        TAB1. 'DEPLACEMENTS_FIBRE_INFERIEURE' . i2. 'TOTAL' = 
  431.          (TAB1. 'DEPLACEMENTS_FIBRE_INFERIEURE' . i2 . 'TOTAL') 
  432.          ET DEPINF;    
  433.      FINSI;
  434. *
  435. *
  436.      DEPTS DEPRS = DEDU DEPT DEPR TRAN MESHORI
  437.             (TAB1. 'MAILLAGE_FIBRE_SUPERIEURE' . I1);
  438.      DEPSUP =  DEPTS +
  439.         ((EXCO DEPRS RY SCAL)*(EXCO VECOPS UZ UX))
  440.       - ((EXCO DEPRS RZ SCAL)*(EXCO VECOPS UY UX))
  441.       + ((EXCO DEPRS RZ SCAL)*(EXCO VECOPS UX UY))
  442.       - ((EXCO DEPRS RX SCAL)*(EXCO VECOPS UZ UY));
  443.      DEPSUP =  DEPSUP +
  444.         ((EXCO DEPRS RX SCAL)*(EXCO VECOPS UY UZ))
  445.       - ((EXCO DEPRS RY SCAL)*(EXCO VECOPS UX UZ));
  446.      DEPSUP = CHAN ATTRIBUT DEPSUP NATURE   DIFFUS;
  447.      TAB1. 'DEPLACEMENTS_FIBRE_SUPERIEURE' . i2 . I1 = DEPSUP;
  448.      SI (I1 EGA 1);
  449.        TAB1. 'DEPLACEMENTS_FIBRE_SUPERIEURE' . i2 . 'TOTAL' 
  450.          = DEPSUP;
  451.      SINON;
  452.        TAB1. 'DEPLACEMENTS_FIBRE_SUPERIEURE' . i2. 'TOTAL' = 
  453.          (TAB1. 'DEPLACEMENTS_FIBRE_SUPERIEURE' . i2 . 'TOTAL') 
  454.          ET DEPSUP;    
  455.      FINSI;
  456. *
  457.      TAB1. 'DEPLACEMENTS_VOLUMIQUE' . i2 . I1 
  458.           = DEPSUP ET DEPINF;
  459.      SI (I1 EGA 1);
  460.        TAB1. 'DEPLACEMENTS_VOLUMIQUE' . i2 . 'TOTAL' 
  461.          = DEPSUP ET DEPINF;
  462.      SINON;
  463.        TAB1. 'DEPLACEMENTS_VOLUMIQUE' . i2. 'TOTAL' = 
  464.          (TAB1. 'DEPLACEMENTS_VOLUMIQUE' . i2 . 'TOTAL') 
  465.          ET DEPSUP ET DEPINF;    
  466.      FINSI;
  467. *
  468.    FIN LAB2; 
  469. *
  470.   FINSI;
  471. *
  472. * Fin de la boucle sur la sous-zone
  473. *
  474. FIN LAB1;
  475. *
  476. MESH3D = TAB1. 'MAILLAGE_VOLUMIQUE'. 'TOTAL';
  477. *
  478. SI FLAGRELA;
  479.      TAB1.'RELATION_3D' = BLCOQ3D;
  480. FINSI;
  481. *
  482. 'FINPROC' MESH3D;
  483. *
  484.  

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