Télécharger la plaquette Cast3M


fr en

Télécharger unilater.procedur

Retour à la liste

Numérotation des lignes : 

  1. * UNILATER  PROCEDUR  MB234859  23/12/08    21:15:04     11806          
  2. *  MODIF JUILLET 1987 POUR PERMETTRE DES JEUX DANS LES LES 2 DIRECTIONS
  3. *
  4. *  CETTE PROCEDURE EST APPELER PAR RESOU
  5. *  LES ARGUMENTS SONT : MATAB TABLE CONTENANT LES INFORMATIONS UTILES
  6. *                         C'EST A DIRE RRR BBB STAT SUP RITOU ET LISEA
  7. *                       FFF   FORCE A METTRE AU SECOND MEMBRE
  8. *
  9. *           RRR   MATRICE DE RIGIDITE SANS LES APPUIS
  10. *           BBB   MATRICE DE RIGIDITE DES APPUIS
  11. *           STAT  LOGIQUE INDIQUANT SI C'EST LE PREMIER PASSAGE
  12. *           SUP   SUPER ELEMENT EQUIVALENT
  13. *           RITOU DERNIERE RAIDEUR UTILISEE DANS LA RESOLUTION
  14. *           LISEA DESCRIPTEUR DES APPUIS COURAMMENT ACTIFS
  15. *           FL    CHAMPS DE FORCES LIMITES DE FROTTEMENT
  16. *
  17. 'DEBP' UNILATER MATAB*'TABLE' FFF*'CHPOINT' MCLE/'MOT' FL/'CHPOINT';
  18.  
  19. *  Si il y a lieu,nbnegs est le nombre de multiplicateurs de lagrange contenus da
  20. ** nbnegs = 0;
  21. ** nbnegt = 0;
  22.   mostab ='MOT' 'NOSTAB';
  23.  
  24.   'SI' ('EGA' 'INIB' MCLE);
  25.     'ARGU' TEMP*'RIGIDITE';
  26.      MATAB.7    = TEMP;
  27.     'ARGU' TEMP*'LISTENTI';
  28.      MATAB.6    = TEMP;
  29.      MATAB.'OK' = VRAI ;
  30.   'FINSI';
  31. *
  32. *  en cas de frottement on lit le champs de force limite de frottement, et eventuellement 
  33. *  la partie non symetrique de la matrice de frottement
  34. *  et la matrice de transfert du multiplicateur de contact au multiplicateur de frottement
  35. *  pour avoir les bonnes reactions
  36. *
  37.   LOG_RF = FAUX;
  38.   LOG_FL='EGA' ('TYPE' FL) 'CHPOINT' ;
  39.   'SI' LOG_FL;
  40.     'ARGU' RFROT/'RIGIDITE';
  41.     LOG_RF='EGA' ('TYPE' RFROT) 'RIGIDITE' ;
  42.     'SI' LOG_RF;
  43.       'ARGU' RTRSF*'RIGIDITE';
  44.       'SI' ('VIDE' RFROT); LOG_RF = FAUX;
  45. **    'MESS' 'rfrot vide dans unilater';
  46.     'FINSI';
  47.     'FINSI';
  48.   'FINSI';
  49.    LOG_ITER = 'EGA' ('VALE' 'RESO') 'ITER';
  50.    RRR =  MATAB.1 ;
  51.    BBB =  MATAB.2 ;
  52.  
  53.   'SI' LOG_FL ;
  54.      FADZ  ='EXCO' FL 'FADH' 'NOID' 'LX' ;
  55.      FL    ='EXCO' FL 'LX'   'NOID' 'LX' ;
  56.      chp2  = BBB * FADZ ;
  57.      FFF = FFF + chp2;
  58.      'DETR' chp2  ;
  59. *1     'SI' LOG_RF;
  60. *1     'SINON';
  61.       chp_z = BBB * FL  ;
  62.       FFF   = FFF + chp_z  ;
  63.       'DETR' chp_z ;
  64. *1     'FINSI';
  65.   'FINSI';
  66.  
  67.    STAT = MATAB.3 ;
  68.   'SI' STAT ;
  69.     'SI' MATAB.'VRAISUP';
  70.        SUP  = MATAB.4 ;
  71.     'SINON';
  72.        RISUP= MATAB.4 ;
  73.     'FINSI';
  74.   'SINON'    ;
  75.  
  76. ************************************************************************
  77. *   CALCUL DU SUPER ELEMENT (SI BESOIN)
  78. ************************************************************************
  79.      MMM ='EXTR' MATAB.1 'MAILLAGE';
  80.      MMB ='EXTR' BBB     'MAILLAGE';
  81. **     MMU = MMM 'ELEM' 'APPU' 'LARG' MMB;
  82.      NTOT='NBNO' MMM ;
  83.      NAPP='NBNO' MMB;
  84. **   MESS ' unilater napp ntot' ' ' NAPP ' ' NTOT ;
  85. **   MESS ' comparaison' ' ' (NAPP ** 2.) ' ' (2 * 6  * NTOT ** 1.333333333) ;
  86. *  Le super element est suppose etre en n^2 alors que la matrice est en n ^4/3
  87.     'SI' (((NAPP ** 2.) > ((30 *NTOT) ** 1.333333)) 'OU' MATAB.'NSYM' 'OU' LOG_ITER);
  88. **     MESS ' unilater pas de super element' napp ntot;
  89.        RISUP          = MATAB.1 ;
  90.        MATAB.4        = MATAB.1 ;
  91.        MATAB.'VRAISUP'= FAUX    ;
  92.  
  93.     'SINON';
  94. **     mess ' on calcule le super element';
  95. *      je remets l'option NOMU car il y a un problème de normalisation avec les mults
  96. *      A voir si elle ne pose pas des problèmes de surcontraintes.
  97.        SUP            ='SUPER' 'RIGI' MATAB.1  BBB ;
  98.        MATAB.4        = SUP  ;
  99.        MATAB.'VRAISUP'= VRAI ;
  100.     'FINSI';
  101.   'FINSI' ;
  102. **
  103. ** sauver le souci courant pour le reinjecter au dernier resou
  104.    souc_co = 'VALE' 'SOUCI';
  105. **
  106.    ISUP   = MATAB.'VRAISUP' ;
  107.    MATAB.3= VRAI;
  108.  
  109. ************************************************************************
  110. * AU CAS OU IL N'Y AURAIT PAS DE DEPIMP
  111. ************************************************************************
  112.    FA     = FFF + ('DEPI' BBB 0.);
  113.    DEPTOT = FA 'EXCO' 'FLX' 'FLX' ;
  114.   'SI' ISUP ;
  115.      RISUP= 'EXTRAI' SUP   'RIGI' ;
  116.      MASUP= 'EXTRAI' RISUP 'MAIL' ;
  117.      DEPTOT = DEPTOT - (DEPTOT REDU MASUP);
  118.      F =   'SUPER'  'CHAR' SUP  FA ;
  119. *    f_lx = f 'EXCO' 'FLX' 'NOID';
  120. *    f_lx_m = f_lx 'EXTR' 'MAIL';
  121. *   'SI' ('NEG' ('TYPE' F_LX_m) 'ENTIER');
  122. **     nbnegs = f_lx_m  nbno;
  123. *   'FINSI';
  124.   'SINON';
  125.      RISUP = MATAB.1 * 1;
  126.      F=FA 'ENLE' 'FLX';
  127.   'FINSI';
  128. * 'SI' ('NON' MATAB.'VRAISUP');
  129. **   nbnegs = 0;
  130. * 'FINSI';
  131.    FDEPTO=F + DEPTOT ;
  132.  
  133. ************************************************************************
  134. *    INITIALISATION
  135. ************************************************************************
  136.    conva1 = vrai;
  137.    'SI' ('NON' STAT);
  138.       MATAB.'OK' = VRAI;
  139.    'FINSI';
  140.    'SI' (STAT 'ET' MATAB.'OK' 'ET' ('NON' LOG_RF)) ;
  141.       RITOU  = MATAB.5   ;
  142.       LISEA  = MATAB.6   ;
  143.       FDEPTOT=FDEPTO;
  144.    'SINON'   ;
  145.      'SI' ('EXISTE' MATAB 7) ;
  146.         BLOINI = MATAB.7;
  147.     'DETR' RITOU;
  148.         RITOU  = RISUP 'ET' MATAB.7 ;
  149.         LISEA  = MATAB.6 'ET' ('LECT') ;
  150.         FDEPTOT=FDEPTO;
  151.      'SINON';
  152.         BLOINI = BBB ;
  153.     'DETR' RITOU;
  154.         RITOU  = RISUP 'ET' BLOINI ;
  155. *       INITIALISER LISEA A UNE VALEUR IMPOSSIBLE
  156.         LISEA  ='LECT' -1 ;
  157.         FDEPTOT= F ;
  158.         conva1 = faux;
  159.      'FINSI';
  160.    'FINSI'   ;
  161. *
  162.     MATAB.'OK' = VRAI;
  163.     matab.5    = 1   ;
  164.     matab.6    = 1   ;
  165.  
  166.     LITOT  ='LECT' ;
  167.     BLOTOT = BBB   ;
  168.  
  169.     nbe = 0;
  170.     MAFREIN = 'VIDE' 'MAILLAGE' / 'POI1';
  171. ************************************************************************
  172. *   RECHERCHE DES BLOQUAGES ACTIFS
  173. ************************************************************************
  174.     OPTION ='MOT' 'RAPIDE';
  175.    'SI' ('EXIS' FL) ;
  176.       OPTIO2='MOT' 'INIT';
  177.    'SINON';
  178.       OPTIO2='MOT' 'CONT';
  179.    'FINSI';
  180.  
  181.     I     = 0    ;
  182.     NCONV = VRAI ;
  183.  
  184.   'SI' LOG_RF ;
  185.        BLOTOT = BLOTOT et rfrot;
  186.   'FINSI';
  187.   'REPETER'   BOUCL1 ;
  188.      I = I + 1 ;
  189. * oublier les soucis des iterations precedentes
  190.      'SOUC' souc_co;
  191.      'SI' LOG_ITER;
  192.         DE = 'RESO' 'SOUCI' 'NOID'  'NOUNIL' RITOU FDEPTOT;
  193.      'SINON';
  194.      MAFREIN*'MAILLAGE'  NBE*'ENTIER' DE = 'RESO' 'SOUCI' 'NOID' 'ENSE' 'NOUNIL' MOSTAB RITOU FDEPTOT ;
  195. **   nbnegt = 'DIAG' ritou;
  196.      'FINSI';
  197.  
  198.     'SI' LOG_FL;
  199.        BLOTO LISEN LOG1 = BLOTOT  'EXCITER' DEPTOT  DE MAFREIN OPTION FL ;
  200.     'SINON';
  201.        BLOTO LISEN LOG1 = BLOTOT  'EXCITER' DEPTOT  DE MAFREIN OPTION ;
  202.     'FINSI';
  203.     'SI' ('NON' ('VERI' DE));
  204.        DE = 'VIDE' 'CHPOINT'/'DIFFUS';
  205.       'MESS' 'mise a zero des deplacements';
  206.       'QUIT' boucl1;
  207.     'FINSI';
  208. **   list lisen;
  209.     'SI' ((I 'EGA' 5) 'ET' ('EGA' OPTION 'RAPIDE'));
  210.        OPTION ='MOT' 'MOYEN';
  211.     'FINSI';
  212.  
  213.     'SI' (('NEG' 'INIB' MCLE) 'OU' (I 'NEG' 1));
  214.       'SI' ((LISEN 'EGA' LISEA) 'OU' LOG1) ;
  215.          'SI' conva1     ;
  216.            'SI' (NBE 'EGA' 0) ;
  217.               NCONV=  FAUX ;
  218.              'QUIT' BOUCL1;
  219.            'SINON';
  220.               'SI' ('NEG' OPTION 'RAPIDE') ;
  221.                 M1 = 'EXTR' (RISUP ET BLOTO) 'MAIL';
  222.                 M2 = 'EXTR' RITOU  'MAIL';
  223.                 M3 = M1 'DIFF' M2;
  224.                'SI' (('NBEL' M3) 'EGA' 0) ;
  225.                   NCONV=  FAUX ;
  226.                  'QUIT' BOUCL1 ;
  227.                'FINSI';
  228.              'FINSI';
  229.            'FINSI';
  230.          'FINSI';
  231.        'FINSI';
  232.          conva1 = vrai;
  233.  
  234.        'SI' ((LISEN 'DANS' LITOT) 'OU' (('MULT' &BOUCL1 100) 'ET' ('NEG' OPTION 'LENT'))) ;
  235.          'SI' ('EGA' OPTION 'LENT');
  236.            'SI' ('EGA' OPTIO2 'CONT');
  237.               OPTIO2= 'MOT' 'INIT';
  238.               OPTION= 'MOT' 'RAPIDE';
  239.               BLOTO=BBB;
  240.  
  241.               LITOT= 'LECT' ;
  242.               LISEN=LISEN * -3;
  243.              'MESS' ' *** UNILATER strategie alternative a la' ' ' I ' eme iteration';
  244.               mostab = 'MOT' 'STAB';
  245.               I=0;
  246.            'SINON';
  247.              'QUIT' BOUCL1;
  248.            'FINSI';
  249.          'SINON';
  250.             OPTION = 'MOT' 'LENT';
  251.             LITOT = 'LECT';
  252.            'MESS' ' *** UNILATER strategie lente a la' ' ' I ' eme iteration';
  253.          'FINSI';
  254.        'FINSI';
  255.     'FINSI';
  256.      LITAN=LITOT;
  257.      LITOT=LITOT 'ET' LISEN;
  258.     'DETR' LITAN;
  259.     'DETR' DE   ;
  260. **'SI' (I 'NEG' 1);
  261.     'DETR' RITOU;
  262. **'FINSI';
  263.  
  264.      RITOU  = RISUP 'ET' BLOTO ;
  265.      BLOT   = BLOTO ;
  266.      LISEA  = LISEN ;
  267.      FDEPTOT= FDEPTO;
  268.   'FIN' BOUCL1;
  269.  
  270. * CFBUG = 'CHAINE' ' *** UNILATER convergence a l''iteration ';
  271. * 'MESS'   CFBUG ' ' I;
  272.   'SI' (NBE 'NEG' 0) ;
  273.      MATAB.'OK' = FAUX;
  274.   'FINSI';
  275.  
  276.    MATAB.5  = RITOU ;
  277.    MATAB.6  = LISEA ;
  278.    MATAB.7  = BLOTO ;
  279. ** MATAB.13 = nbnegt - nbnegs;
  280.  
  281.   'SI'  NCONV  ;
  282.  
  283.      NCONV     =FAUX;
  284.      MATAB.'OK'=FAUX;
  285.   'FINSI';
  286.  
  287.   'SI'  NCONV  ;
  288.  
  289.     'MESS'  ' *** UNILATER il n y a pas de solution au systeme ';
  290.      RITOU = RISUP 'ET' BBB ;
  291.      LISEA = 'LECT' 0;
  292.      MATAB.5 =  RITOU ;
  293.      MATAB.6 =  LISEA ;
  294.      MATAB.7 = BLOTO;
  295.       'SOUCI' 2345;
  296.  
  297.   'SINON';
  298. *    CFBUG = 'CHAINE' ' *** UNILATER convergence a l''iteration ';
  299.     'SI' ('NON' MATAB.'OK');
  300.        CFBUG = 'CHAINE' ' *** UNILATER non convergence a l''iteration ';
  301.       'MESS' CFBUG I ' ' nbe;
  302.        MATAB.13 = 1119 ;
  303.       'SOUCI' 1119;
  304.     'FINSI';
  305.      MATAB.'NITER' = I;
  306.  
  307. ************************************************************************
  308. *    CALCUL SUR TOUTE LA STRUCTURE
  309. ************************************************************************
  310.     'SI' ISUP;
  311.        DE1 ='SUPER' 'DEPL' SUP DE FFF;
  312.        DE2 = DE1 'ET' (DE 'ENLE' 'LX') 'NOER';
  313.       'DETR' DE1;
  314.       'DETR' DE ;
  315.       'DETR' FA ;
  316.     'SINON';
  317.        DE2 = DE ;
  318.     'FINSI' ;
  319.  
  320.     'SI' LOG_FL;
  321.        DE3 = DE2 - FADZ ;
  322.        'DETR' DE2;
  323.        DE2 = DE3;
  324.      'SI' LOG_RF;
  325. *  transfert de la reaction de contact sur la reaction de frottement
  326.       DE3 = RTRSF * DE2 ;
  327.       DE3 = DE3 CHAN ATTRIBUT NATURE DIFFUS;
  328. ******   DE3 = DE3 'REDU' (BBB 'EXTR' 'MAIL' 'MULT');
  329.          DE2 = DE2 +  DE3;
  330.          'DETR' DE3;
  331.       DE3 = DE2;
  332.      'FINSI';
  333. *1   'SINON';
  334.        FL  = FL 'REDU' (BBB 'EXTR' 'MAIL' 'MULT');
  335.        DE2 = DE3 - FL;
  336.       'DETR' DE3 ;
  337. *1   'FINSI';
  338.     'FINSI';
  339.   'FINSI' ;
  340.  
  341.   'DETR' F  ;
  342.   'DETR' DEPTOT ;
  343.   'DETR' FDEPTOT ;
  344.  
  345.   'SI' ( 'EXISTE' MATAB 50) ;
  346.     DE3 = 'MREM'  DE2  matab.50 matab.10 matab.11;
  347.    'RESPRO' DE3;
  348.   'SINON' ;
  349. *    mess ' dans unilater de2 ';list de2;
  350.     'RESPRO'  DE2  ;
  351.   'FINSI' ;
  352.  
  353.   'SI' ('EGA' 'INIB' MCLE);
  354.     'SI' LOG_FL ;
  355.       'RESP' 'bid' 'bid' 'bid' 'bid';
  356.     'SINON';
  357.       'RESP' 'bid' 'bid' 'bid'      ;
  358.     'FINSI';
  359.   'FINSI';
  360. 'FINP'  ;
  361.  
  362.  
  363.  

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