Procedures | Cast3M

   1 : * PMPB      PROCEDUR  BAZE      98/11/20    21:22:25     3365           
   2 : DEBPROC PMPB TSIG*chpoint  LIG*maillage  grapl*logique
   3 :  nnc*entier  nne*entier ;
   4 : *************************************************************************
   5 : *   Analyse des contraintes sur le segment d appui  LIG  ( SEG2 )suivant
   6 : *                     specification des codes  
   7 : *   ENTREES :    TSIG chpoint des contraintes projetees sur LIG
   8 : *                        (TSIG = PROI LIG  SIG ) 
   9 : *                LIG  segment d'appui
  10 : *                GRAPL Logique   si VRAI --> traces des decompositions  
  11 : *   SORTIES :
  12 : *       TABLES contenant pour chacune des composantes des contraintes
  13 : *              caractéristiques sur les extrémites de segments d'appui
  14 : *     
  15 : ************************************************************************
  16 : P1 = lig point init ;
  17 : P2 = lig point fina ;
  18 : Listcomp = extr Tsig comp ;
  19 : nn = lig nbno ;
  20 : idi = vale DIME ;
  21 : si (ega idi 2);
  22 : x y = coor lig ;
  23 : z = x * 0. ;
  24 : sinon ;
  25 : x y z =   coor lig   ;
  26 : z = z - ( coor 3 P1) ;
  27 : finsi ;
  28 : x = x - ( coor 1 P1);
  29 : y = y - ( coor 2 P1);
  30 : xabs  = ((x * x) + ( y * y) + ( z * z)) ** .5 ;
  31 : ep = maxi xabs ;
  32 : SAUT  LIGNE ;
  33 : MESS '        ETAT ' NNE ' COUPE' NNC ': CONTRAINTES LINEARISEES ' ;
  34 : SAUT  LIGNE ;
  35 : nbcomp = dime Listcomp  ;
  36 : compm  = Table ;
  37 : compb  = Table ;
  38 : compmb = Table ;
  39 :   MESS '          ----------------------------------------------'
  40 :   '----------' ;
  41 :   MESS '          ! MEMBRANE  ! MEM+FLEX ! MEM+FLEX ! TOTALES '
  42 :   ' ! TOTALES  !' ;
  43 :   MESS '          !           !    INT.  !    EXT.  !   INT. '
  44 :   '  !   EXT.   !' ;
  45 :  TIRET= 
  46 : '-----------------------------------------------------------------';
  47 :    MESS TIRET; 
  48 : *
  49 : REPETER BOU1 NBCOMP ;
  50 :  CCC   = extr Listcomp &bou1 ;
  51 :  ev1 = evol roug chpo TSIG CCC lig ;
  52 :  rabs = extr ev1 absc ; 
  53 : *  ep =  maxi rabs  ;                          
  54 : *
  55 : *------------------------   contrainte moyenne membrane
  56 : *
  57 :    pm =(extr (somm ev1 ) 1 ) / ep  ;
  58 :   compm.&bou1 = manu chpo lig 1 SCAL  pm ;                                    
  59 :   
  60 : *     
  61 : *-----------------------    calcul du moment de flexion
  62 : * 
  63 :    xl = rabs - ( prog nn * (ep / 2 ))  ;
  64 :    z = manu chpo lig 1 SCAL xl ;
  65 :    sigzz= ((exco TSIG CCC) * z ) nomc SCAL ;
  66 :    evm = evol bleu chpo sigzz SCAL lig  ;
  67 :    mm =  extr (somm evm ) 1 ;                                              
  68 : *-----------------------------   contrainte   flexion
  69 :    pb    = mm * 6. / ( ep * ep) ;                                       
  70 :    a =( pb * 2. )  / ep  ;
  71 :    compb.&bou1 = ((xabs * a) - pb ) ;                       
  72 : * 
  73 :  compmb.&bou1 = compm.&bou1 + compb.&bou1 ;
  74 :  MEMBR = extr compm.&bou1 SCAL  P1 ;
  75 :  MEFLRI =(extr compmb.&bou1  SCAL  P1 )  ; 
  76 :  MEFLRE =(extr compmb.&bou1  SCAL  P2 )  ; 
  77 :  FLEXI = (extr  compb.&bou1  SCAL  P1 )  ;
  78 :  FLEXE = (extr  compb.&bou1  SCAL  P2 )  ;  
  79 :  TOTIR = EXTR TSIG CCC  P1 ;
  80 :  TOTER = EXTR TSIG CCC  P2 ;
  81 : 
  82 : 
  83 :  si GRAPL  ;
  84 : *
  85 :  'TITRE' 'Seg_app' nnc 'Etat' nne  CCC  '  Membrane ' membr  
  86 :          'Flexion max ' (MEFLRI - membr)   ;  
  87 :  evpm  = evol vert chpo  compm.&bou1    SCAL lig  ;
  88 :  evpmb = evol turq chpo  compmb.&bou1   SCAL lig ;
  89 :  evsignl = ev1 - evpmb ;
  90 :  evt =  ev1 et evpm et evpmb  et evsignl ;
  91 : *  
  92 :    TUB = TABLE ;TIB = TABLE ;
  93 :    TUB.1 = 'MARQ ETOI'; TIB.1 =  'Totale ' ;  
  94 :    TUB.2 = 'MARQ CROI'; TIB.2 =  'Pm ' ;
  95 :    TUB.3 = 'MARQ TRIA'; TIB.3 =  'Pm+Pb ' ;
  96 :    TUB.4 = 'MARQ CARR'; TIB.4 =  'NonLin ' ;
  97 :    TUB.'TITRE' = TIB ;   
  98 : *
  99 :    'DESS' evt lege tub xbor 0. ep    AXES ;
 100 :  finsi ;
 101 : 
 102 : AA = CHAINE FORMAT '(F9.2)' ' ! ' CCC   '   ! 'MEMBR ' !' MEFLRI
 103 : ' !' MEFLRE ' !' TOTIR ' !' TOTER ' !';
 104 :     MESS AA;
 105 :     MESS TIRET; 
 106 :    SI ( EGA &BOU1 1) ;
 107 :  P_MEMB  = PROG  MEMBR ;
 108 :  P_FLEXI = PROG  FLEXI ;    
 109 :  P_MPFI   = PROG  MEFLRI ;
 110 :  P_MPFE   = PROG  MEFLRE;
 111 :  P_FLEXE = PROG  FLEXE ;
 112 :  P_TOTI  = PROG  TOTIR  ;
 113 :  P_TOTE  = PROG  TOTER  ;
 114 :    SINON ;
 115 :  P_MEMB  = P_MEMB  et (PROG  MEMBR  ) ;  
 116 :  P_FLEXI = P_FLEXI et (PROG  FLEXI  ) ;   
 117 :  P_MPFI  = P_MPFI  et (PROG  MEFLRI ) ;
 118 :  P_MPFE  = P_MPFE  et (PROG  MEFLRE ) ;
 119 :  P_FLEXE = P_FLEXE et (PROG  FLEXE  ) ;
 120 :  P_TOTI  = P_TOTI  et (PROG  TOTIR  ) ;
 121 :  P_TOTE  = P_TOTE  et (PROG  TOTER  ) ;  
 122 :    FINSI ;
 123 :               
 124 : FIN BOU1 ;
 125 : *
 126 :  TABV = TABLE ;
 127 :  TABV . 1 = P_FLEXI  ;
 128 :  TABV . 2 = P_MPFI ;
 129 :  TABV . 3 = P_TOTI  ;
 130 : *
 131 :  TABV . 4 = P_FLEXE  ;
 132 :  TABV . 5 = P_MPFE ;
 133 :  TABV . 6 = P_TOTE ;
 134 : *
 135 : * Contraintes équivalentes - Critère de Von Mises  en peau 
 136 : * -----------------------------------------------
 137 : *                          - Critère de Tresca en peau 
 138 : * -------------------------------------------- 
 139 : *
 140 :  tabana = table ;
 141 :  tabana.1 = P_MEMB ;
 142 :  tabana.2 = P_MPFI;
 143 :  tabana.3 = P_MPFE ; 
 144 :  tabana.4 = P_TOTI ;
 145 :  tabana.5 = P_TOTE ;  
 146 :     REPETER BTB 5 ;
 147 : prcour = tabana.&btb ;    
 148 : S1 = EXTR prcour 1 ;S2 = EXTR prcour 2 ;S3 = EXTR prcour 3 ;
 149 : S4 = EXTR prcour 4 ;S5= 0. ;S6 = 0. ;
 150 : SI ( EGA NBCOMP 6 ) ;
 151 : S5 = EXTR prcour 5 ;S6 = EXTR prcour 6 ;
 152 : FINSI ;
 153 : *
 154 : TRMM VMMM = @RCCMTRV S1 S2 S3 S4 S5 S6 ;
 155 :    si ( ega &btb 1 ) ;
 156 :    bbb = prog TRMM ;
 157 :    VVV = prog VMMM ;
 158 :    sinon ;
 159 :    bbb = bbb et (prog TRMM ) ;
 160 :    VVV = VVV et(prog VMMM ) ; 
 161 :    finsi ;   
 162 :     'FIN' BTB ;
 163 : TRESMI = 'EXTR' BBB 1 ; MEMISI  = 'EXTR' VVV 1 ;
 164 : TREMFI = 'EXTR' BBB 2 ; BEMISI  = 'EXTR' VVV 2 ;
 165 : TREMFE = 'EXTR' BBB 3 ; BEMISE  = 'EXTR' VVV 3 ;
 166 : TTREI  = 'EXTR' BBB 4 ; VMV1    = 'EXTR' VVV 4 ;
 167 : TTREE  = 'EXTR' BBB 5 ; VMV2    = 'EXTR' VVV 5 ;
 168 :     MESS TIRET;  
 169 : AA = CHAINE FORMAT '(F9.2)' ' ! ' TRESCA ' ! 'TRESMI ' !' TREMFI
 170 : ' !' TREMFE  ' !' TTREI  ' !' TTREE  ' !';
 171 :     MESS AA;
 172 :     MESS TIRET;      
 173 : AA = CHAINE FORMAT '(F9.2)' ' ! ' VMISES ' ! 'MEMISI ' !' BEMISI
 174 : ' !'BEMISE  ' !' VMV1 ' !' VMV2 ' !';
 175 :     MESS AA;
 176 :     MESS TIRET;
 177 : 'FINPROC'  TABV;
 178 :