Télécharger la plaquette Cast3M


fr en

Télécharger unilater.procedur

Retour à la liste

Numérotation des lignes : 

  1. * UNILATER  PROCEDUR  PV090527  24/10/24    21:15:01     12047          
  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.         DE = 'RESO' 'SOUCI' 'NOID'  'NOUNIL' RITOU FDEPTOT;
  188.      'SINON';
  189.      MAFREIN*'MAILLAGE'  NBE*'ENTIER' DE = 'RESO' 'SOUCI' 'NOID' 'ENSE' 'NOUNIL' MOSTAB RITOU FDEPTOT ;
  190. **   nbnegt = 'DIAG' ritou;
  191.      'FINSI';
  192.  
  193.     'SI' LOG_FL;
  194.        BLOTO LISEN LOG1 = BLOTOT  'EXCITER' DEPTOT  DE MAFREIN OPTION FL ;
  195.     'SINON';
  196.        BLOTO LISEN LOG1 = BLOTOT  'EXCITER' DEPTOT  DE MAFREIN OPTION ;
  197.     'FINSI';
  198.     'SI' ('NON' ('VERI' DE));
  199.        DE = 'VIDE' 'CHPOINT'/'DIFFUS';
  200.       'MESS' 'mise a zero des deplacements';
  201.       'QUIT' boucl1;
  202.     'FINSI';
  203. **   list lisen;
  204.     'SI' ((I 'EGA' 5) 'ET' ('EGA' OPTION 'RAPIDE'));
  205.        OPTION ='MOT' 'MOYEN';
  206.     'FINSI';
  207.  
  208.     'SI' (('NEG' 'INIB' MCLE) 'OU' (I 'NEG' 1));
  209.       'SI' ((LISEN 'EGA' LISEA) 'OU' LOG1) ;
  210.          'SI' conva1     ;
  211.            'SI' (NBE 'EGA' 0) ;
  212.               NCONV=  FAUX ;
  213.              'QUIT' BOUCL1;
  214.            'SINON';
  215.               'SI' ('NEG' OPTION 'RAPIDE') ;
  216.                 M1 = 'EXTR' (RISUP ET BLOTO) 'MAIL';
  217.                 M2 = 'EXTR' RITOU  'MAIL';
  218.                 M3 = M1 'DIFF' M2;
  219.                'SI' (('NBEL' M3) 'EGA' 0) ;
  220.                   NCONV=  FAUX ;
  221.                  'QUIT' BOUCL1 ;
  222.                'FINSI';
  223.              'FINSI';
  224.            'FINSI';
  225.          'FINSI';
  226.        'FINSI';
  227.          conva1 = vrai;
  228.  
  229.        'SI' ((LISEN 'DANS' LITOT) 'OU' (('MULT' &BOUCL1 100) 'ET' ('NEG' OPTION 'LENT'))) ;
  230.          'SI' ('EGA' OPTION 'LENT');
  231.            'SI' ('EGA' OPTIO2 'CONT');
  232.               OPTIO2= 'MOT' 'INIT';
  233.               OPTION= 'MOT' 'RAPIDE';
  234.               BLOTO=BBB;
  235.  
  236.               LITOT= 'LECT' ;
  237.               LISEN=LISEN * -3;
  238.              'MESS' ' *** UNILATER strategie alternative a la' ' ' I ' eme iteration';
  239.               mostab = 'MOT' 'STAB';
  240.               I=0;
  241.            'SINON';
  242.              'QUIT' BOUCL1;
  243.            'FINSI';
  244.          'SINON';
  245.             OPTION = 'MOT' 'LENT';
  246.             LITOT = 'LECT';
  247.            'MESS' ' *** UNILATER strategie lente a la' ' ' I ' eme iteration';
  248.          'FINSI';
  249.        'FINSI';
  250.     'FINSI';
  251.      LITAN=LITOT;
  252.      LITOT=LITOT 'ET' LISEN;
  253.     'DETR' LITAN;
  254.     'DETR' DE   ;
  255. **'SI' (I 'NEG' 1);
  256.     'DETR' RITOU;
  257. **'FINSI';
  258.  
  259.      RITOU  = RISUP 'ET' BLOTO ;
  260.      BLOT   = BLOTO ;
  261.      LISEA  = LISEN ;
  262.      FDEPTOT= FDEPTO;
  263.   'FIN' BOUCL1;
  264.  
  265. * CFBUG = 'CHAINE' ' *** UNILATER convergence a l''iteration ';
  266. * 'MESS'   CFBUG ' ' I;
  267.   'SI' (NBE 'NEG' 0) ;
  268.      MATAB.'OK' = FAUX;
  269.   'FINSI';
  270.  
  271.    MATAB.5  = RITOU ;
  272.    MATAB.6  = LISEA ;
  273.    MATAB.7  = BLOTO ;
  274. ** MATAB.13 = nbnegt - nbnegs;
  275.  
  276.   'SI'  NCONV  ;
  277.  
  278.      NCONV     =FAUX;
  279.      MATAB.'OK'=FAUX;
  280.   'FINSI';
  281.  
  282.   'SI'  NCONV  ;
  283.  
  284.     'MESS'  ' *** UNILATER il n y a pas de solution au systeme ';
  285.      RITOU = RISUP 'ET' BBB ;
  286.      LISEA = 'LECT' 0;
  287.      MATAB.5 =  RITOU ;
  288.      MATAB.6 =  LISEA ;
  289.      MATAB.7 = BLOTO;
  290.       'SOUCI' 2345;
  291.  
  292.   'SINON';
  293. *    CFBUG = 'CHAINE' ' *** UNILATER convergence a l''iteration ';
  294.     'SI' ('NON' MATAB.'OK');
  295.        CFBUG = 'CHAINE' ' *** UNILATER non convergence a l''iteration ';
  296.       'MESS' CFBUG I ' ' nbe;
  297.        MATAB.13 = 1119 ;
  298.       'SOUCI' 1119;
  299.     'FINSI';
  300.      MATAB.'NITER' = I;
  301.  
  302. ************************************************************************
  303. *    CALCUL SUR TOUTE LA STRUCTURE
  304. ************************************************************************
  305.     'SI' ISUP;
  306.        DE1 ='SUPER' 'DEPL' SUP DE FFF;
  307.        DE2 = DE1 'ET' (DE 'ENLE' 'LX') 'NOER';
  308.       'DETR' DE1;
  309.       'DETR' DE ;
  310.       'DETR' FA ;
  311.     'SINON';
  312.        DE2 = DE ;
  313.     'FINSI' ;
  314.  
  315.     'SI' LOG_FL;
  316.        DE3 = DE2 - FADZ ;
  317.        'DETR' DE2;
  318.        DE2 = DE3;
  319.      'SI' LOG_RF;
  320. *  transfert de la reaction de contact sur la reaction de frottement
  321.       DE3 = RTRSF * DE2 ;
  322.       DE3 = DE3 CHAN ATTRIBUT NATURE DIFFUS;
  323. ******   DE3 = DE3 'REDU' (BBB 'EXTR' 'MAIL' 'MULT');
  324.          DE2 = DE2 +  DE3;
  325.          'DETR' DE3;
  326.       DE3 = DE2;
  327.      'FINSI';
  328. *1   'SINON';
  329.        FL  = FL 'REDU' (BBB 'EXTR' 'MAIL' 'MULT');
  330.        DE2 = DE3 - FL;
  331.       'DETR' DE3 ;
  332. *1   'FINSI';
  333.     'FINSI';
  334.   'FINSI' ;
  335.  
  336.   'DETR' F  ;
  337.   'DETR' DEPTOT ;
  338.   'DETR' FDEPTOT ;
  339.  
  340.   'SI' ( 'EXISTE' MATAB 50) ;
  341.     DE3 = 'MREM'  DE2  matab.50 matab.10 matab.11;
  342.    'RESPRO' DE3;
  343.   'SINON' ;
  344. *    mess ' dans unilater de2 ';list de2;
  345.     'RESPRO'  DE2  ;
  346.   'FINSI' ;
  347.  
  348.   'SI' ('EGA' 'INIB' MCLE);
  349.     'SI' LOG_FL ;
  350.       'RESP' 'bid' 'bid' 'bid' 'bid';
  351.     'SINON';
  352.       'RESP' 'bid' 'bid' 'bid'      ;
  353.     'FINSI';
  354.   'FINSI';
  355. 'FINP'  ;
  356.  
  357.  
  358.  
  359.  
  360.  
  361.  

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