Télécharger la plaquette Cast3M


fr en

Télécharger unilater.procedur

Retour à la liste

Numérotation des lignes : 

  1. * UNILATER  PROCEDUR  MB234859  24/11/20    21:15:23     12080          
  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. *  test reporte dans supri
  88.     'SI' (((NAPP ** 2.) > ((30 *NTOT) ** 1.333333)) 'OU' MATAB.'NSYM' 'OU' LOG_ITER 'ET' FAUX);
  89. **     MESS ' unilater pas de super element' napp ntot;
  90.        RISUP          = MATAB.1 ;
  91.        MATAB.4        = MATAB.1 ;
  92.        MATAB.'VRAISUP'= FAUX    ;
  93.  
  94.     'SINON';
  95. **     mess ' on calcule le super element';
  96. *      je remets l'option NOMU car il y a un problème de normalisation avec les mults
  97. *      A voir si elle ne pose pas des problèmes de surcontraintes.
  98.        SUP            ='SUPER' 'RIGI' MATAB.1  BBB ;
  99.        MATAB.4        = SUP  ;
  100.        MATAB.'VRAISUP'= VRAI ;
  101.     'FINSI';
  102.   'FINSI' ;
  103. **
  104. ** sauver le souci courant pour le reinjecter au dernier resou
  105.    souc_co = 'VALE' 'SOUCI';
  106. **
  107.    ISUP   = MATAB.'VRAISUP' ;
  108.    MATAB.3= VRAI;
  109.  
  110. ************************************************************************
  111. * AU CAS OU IL N'Y AURAIT PAS DE DEPIMP
  112. ************************************************************************
  113.    FA     = FFF + ('DEPI' BBB 0.);
  114.    DEPTOT = FA 'EXCO' 'FLX' 'FLX' ;
  115.   'SI' ISUP ;
  116.      RISUP= 'EXTRAI' SUP   'RIGI' ;
  117.      MASUP= 'EXTRAI' RISUP 'MAIL'   ;
  118.      F =   'SUPER'  'CHAR' SUP  FA ;
  119.      F=F  'ENLE' 'FLX';
  120. **  DEPTOT = DEPTOT - (DEPTOT REDU MASUP);
  121.   'SINON';
  122.      RISUP = MATAB.1 * 1;
  123.      F=FA 'ENLE' 'FLX';
  124.   'FINSI';
  125.  
  126.    FDEPTO=F + DEPTOT ;
  127.  
  128. ************************************************************************
  129. *    INITIALISATION
  130. ************************************************************************
  131.    conva1 = vrai;
  132.    'SI' ('NON' STAT);
  133.       MATAB.'OK' = VRAI;
  134.    'FINSI';
  135.    'SI' (STAT 'ET' MATAB.'OK' 'ET' ('NON' LOG_RF)) ;
  136.       RITOU  = MATAB.5   ;
  137.       LISEA  = MATAB.6   ;
  138.       FDEPTOT=FDEPTO;
  139.    'SINON'   ;
  140.      'SI' ('EXISTE' MATAB 7) ;
  141.         BLOINI = MATAB.7;
  142.     'DETR' RITOU;
  143.         RITOU  = RISUP 'ET' MATAB.7 ;
  144.         LISEA  = MATAB.6 'ET' ('LECT') ;
  145.         FDEPTOT=FDEPTO;
  146.      'SINON';
  147.         BLOINI = BBB ;
  148.     'DETR' RITOU;
  149.         RITOU  = RISUP 'ET' BLOINI ;
  150. *       INITIALISER LISEA A UNE VALEUR IMPOSSIBLE
  151.         LISEA  ='LECT' -1 ;
  152.         FDEPTOT= F ;
  153.         conva1 = faux;
  154.      'FINSI';
  155.    'FINSI'   ;
  156. *
  157.     MATAB.'OK' = VRAI;
  158.     matab.5    = 1   ;
  159.     matab.6    = 1   ;
  160.  
  161.     LITOT  ='LECT' ;
  162.     BLOTOT = BBB   ;
  163.  
  164.     nbe = 0;
  165.     MAFREIN = 'VIDE' 'MAILLAGE' / 'POI1';
  166. ************************************************************************
  167. *   RECHERCHE DES BLOQUAGES ACTIFS
  168. ************************************************************************
  169.     OPTION ='MOT' 'RAPIDE';
  170.    'SI' ('EXIS' FL) ;
  171.       OPTIO2='MOT' 'INIT';
  172.    'SINON';
  173.       OPTIO2='MOT' 'CONT';
  174.    'FINSI';
  175.  
  176.     I     = 0    ;
  177.     NCONV = VRAI ;
  178.  
  179.   'SI' LOG_RF ;
  180.        BLOTOT = BLOTOT et rfrot;
  181.   'FINSI';
  182.   'REPETER'   BOUCL1 ;
  183.     I = I + 1 ;
  184. * oublier les soucis des iterations precedentes
  185.     'SOUC' souc_co;
  186.     'SI' LOG_ITER;
  187.       NBE = 0;
  188.       CHFREIN = 'VIDE' CHPOINT;
  189.       DE = 'RESO' 'SOUCI' 'NOID'  'NOUNIL' RITOU FDEPTOT;
  190.     'SINON';
  191.       NBE*'ENTIER' DE*CHPOINT CHFREIN*CHPOINT  = 'RESO' 'SOUCI' 'NOID' 'ENSE' 'NOUNIL' MOSTAB RITOU FDEPTOT ;
  192.     'FINSI';
  193.     maens = (chFrein 'EXCO' 'LX' 'LX') 'EXTR' 'MAILLAGE';
  194.     'SI' ('>' ('NBNO' maens) 0) ;
  195.       fens = 'MANU' 'CHPO' maens 1 'FLX' 1.;
  196.       NBF*'ENTIER' CHFREIN*CHPOINT BID*CHPOINT  = 'RESO' 'SOUCI' 'NOID' 'ENSE' 'NOUNIL' MOSTAB RITOU FENS ;
  197.       mess 'unilater correction mode d ensemble iteration' ' ' I ;
  198.     'FINSI';
  199.  
  200.     'SI' LOG_FL;
  201.       BLOTO LISEN LOG1 = BLOTOT 'EXCITER' DEPTOT DE CHFREIN OPTION FL ;
  202.     'SINON';
  203.       BLOTO LISEN LOG1 = BLOTOT 'EXCITER' DEPTOT DE CHFREIN OPTION ;
  204.     'FINSI';
  205. **  list lisen;
  206.     'SI' ('NON' ('VERI' DE));
  207.       DE = 'VIDE' 'CHPOINT'/'DIFFUS';
  208.       'MESS' 'mise a zero des deplacements';
  209.       'QUIT' boucl1;
  210.     'FINSI';
  211. *
  212.     'SI' ((I 'EGA' 5) 'ET' ('EGA' OPTION 'RAPIDE'));
  213.        OPTION ='MOT' 'MOYEN';
  214.     'FINSI';
  215. *
  216.     'SI' (('NEG' 'INIB' MCLE) 'OU' (I 'NEG' 1));
  217.       'SI' ((LISEN 'EGA' LISEA) 'OU' LOG1) ;
  218.          'SI' conva1     ;
  219.            'SI' (NBE 'EGA' 0) ;
  220.               NCONV=  FAUX ;
  221.              'QUIT' BOUCL1;
  222.            'SINON';
  223.               'SI' ('NEG' OPTION 'RAPIDE') ;
  224.                 M1 = 'EXTR' (RISUP ET BLOTO) 'MAIL';
  225.                 M2 = 'EXTR' RITOU  'MAIL';
  226.                 M3 = M1 'DIFF' M2;
  227.                'SI' (('NBEL' M3) 'EGA' 0) ;
  228.                   NCONV=  FAUX ;
  229.                  'QUIT' BOUCL1 ;
  230.                'FINSI';
  231.              'FINSI';
  232.            'FINSI';
  233.          'FINSI';
  234.        'FINSI';
  235.          conva1 = vrai;
  236.  
  237.        'SI' ((LISEN 'DANS' LITOT) 'OU' (('MULT' &BOUCL1 100) 'ET' ('NEG' OPTION 'LENT'))) ;
  238.          'SI' ('EGA' OPTION 'LENT');
  239.            'SI' ('EGA' OPTIO2 'CONT');
  240.               OPTIO2= 'MOT' 'INIT';
  241.               OPTION= 'MOT' 'RAPIDE';
  242.               BLOTO=BBB;
  243.  
  244.               LITOT= 'LECT' ;
  245.               LISEN=LISEN * -3;
  246.              'MESS' ' *** UNILATER strategie alternative a la' ' ' I ' eme iteration';
  247.               mostab = 'MOT' 'STAB';
  248.               I=0;
  249.            'SINON';
  250.              'QUIT' BOUCL1;
  251.            'FINSI';
  252.          'SINON';
  253.             OPTION = 'MOT' 'LENT';
  254.             LITOT = 'LECT';
  255.            'MESS' ' *** UNILATER strategie lente a la' ' ' I ' eme iteration';
  256.          'FINSI';
  257.        'FINSI';
  258.     'FINSI';
  259.      LITAN=LITOT;
  260.      LITOT=LITOT 'ET' LISEN;
  261.     'DETR' LITAN;
  262.     'DETR' DE   ;
  263. **'SI' (I 'NEG' 1);
  264.     'DETR' RITOU;
  265. **'FINSI';
  266.  
  267.      RITOU  = RISUP 'ET' BLOTO ;
  268.      BLOT   = BLOTO ;
  269.      LISEA  = LISEN ;
  270.      FDEPTOT= FDEPTO;
  271.   'FIN' BOUCL1;
  272.  
  273. * CFBUG = 'CHAINE' ' *** UNILATER convergence a l''iteration ';
  274. * 'MESS'   CFBUG ' ' I;
  275.   'SI' (NBE 'NEG' 0) ;
  276.      MATAB.'OK' = FAUX;
  277.   'FINSI';
  278.  
  279.    MATAB.5  = RITOU ;
  280.    MATAB.6  = LISEA ;
  281.    MATAB.7  = BLOTO ;
  282. ** MATAB.13 = nbnegt - nbnegs;
  283.  
  284.   'SI'  NCONV  ;
  285.  
  286.      NCONV     =FAUX;
  287.      MATAB.'OK'=FAUX;
  288.   'FINSI';
  289.  
  290.   'SI'  NCONV  ;
  291.  
  292.     'MESS'  ' *** UNILATER il n y a pas de solution au systeme ';
  293.      RITOU = RISUP 'ET' BBB ;
  294.      LISEA = 'LECT' 0;
  295.      MATAB.5 =  RITOU ;
  296.      MATAB.6 =  LISEA ;
  297.      MATAB.7 = BLOTO;
  298.       'SOUCI' 2345;
  299.  
  300.   'SINON';
  301. *    CFBUG = 'CHAINE' ' *** UNILATER convergence a l''iteration ';
  302.     'SI' ('NON' MATAB.'OK');
  303.        CFBUG = 'CHAINE' ' *** UNILATER non convergence a l''iteration ';
  304.       'MESS' CFBUG ' ' I ' ' nbe;
  305.        MATAB.13 = 1119 ;
  306.       'SOUCI' 1119;
  307.     'FINSI';
  308.      MATAB.'NITER' = I;
  309.  
  310. ************************************************************************
  311. *    CALCUL SUR TOUTE LA STRUCTURE
  312. ************************************************************************
  313.     'SI' ISUP;
  314.        DE1 ='SUPER' 'DEPL' SUP DE FFF;
  315.        DE2 = DE1 'ET' (DE 'ENLE' 'LX') 'NOER';
  316.       'DETR' DE1;
  317.       'DETR' DE ;
  318.       'DETR' FA ;
  319.     'SINON';
  320.        DE2 = DE ;
  321.     'FINSI' ;
  322.  
  323.     'SI' LOG_FL;
  324.        DE3 = DE2 - FADZ ;
  325.        'DETR' DE2;
  326.        DE2 = DE3;
  327.      'SI' LOG_RF;
  328. *  transfert de la reaction de contact sur la reaction de frottement
  329.       DE3 = RTRSF * DE2 ;
  330.       DE3 = DE3 CHAN ATTRIBUT NATURE DIFFUS;
  331. ******   DE3 = DE3 'REDU' (BBB 'EXTR' 'MAIL' 'MULT');
  332.          DE2 = DE2 +  DE3;
  333.          'DETR' DE3;
  334.       DE3 = DE2;
  335.      'FINSI';
  336. *1   'SINON';
  337.        FL  = FL 'REDU' (BBB 'EXTR' 'MAIL' 'MULT');
  338.        DE2 = DE3 - FL;
  339.       'DETR' DE3 ;
  340. *1   'FINSI';
  341.     'FINSI';
  342.   'FINSI' ;
  343.  
  344.   'DETR' F  ;
  345.   'DETR' DEPTOT ;
  346.   'DETR' FDEPTOT ;
  347.  
  348.   'SI' ( 'EXISTE' MATAB 50) ;
  349.     DE3 = 'MREM'  DE2  matab.50 matab.10 matab.11;
  350.    'RESPRO' DE3;
  351.   'SINON' ;
  352. *    mess ' dans unilater de2 ';list de2;
  353.     'RESPRO'  DE2  ;
  354.   'FINSI' ;
  355.  
  356.   'SI' ('EGA' 'INIB' MCLE);
  357.     'SI' LOG_FL ;
  358.       'RESP' 'bid' 'bid' 'bid' 'bid';
  359.     'SINON';
  360.       'RESP' 'bid' 'bid' 'bid'      ;
  361.     'FINSI';
  362.   'FINSI';
  363. 'FINP'  ;
  364.  
  365.  
  366.  
  367.  
  368.  
  369.  
  370.  
  371.  
  372.  

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