Download @ortho.procedur

Back to the list

   1 : * @ORTHO    PROCEDUR  FANDEUR   09/03/31    21:15:00     6315           
   2 : *-------------------------------------------------
   3 : ******          PROCEDURE @ORTHO            ******
   4 : *-------------------------------------------------
   5 : *
   6 : *    CETTE PROCEDURE A ETE MISE GRACIEUSEMENT
   7 : *   A DISPOSOTION DE LA COMMUNAUTE  CASTEM2000
   8 : *    PAR M. LIBEYRE ( CEA/DSM/DRFC )
   9 : *
  10 : *     TEL : ( 33 1 ) 42 25 46 03
  11 : *
  12 : *  PROCEDURE MISE A JOUR PAR OF LE 09/03/2009
  13 : *-------------------------------------------------
  14 : DEBPROC @ORTHO1 MAIL1*'MAILLAGE' LIG1*'MAILLAGE'
  15 :                 CH1*'CHPOINT' CH2*'CHPOINT'
  16 :                 ALPHA*'FLOTTANT' LL*'LISTREEL'
  17 :                 TYPEMOD/'MOT' TYPEELT/'MOT' PTGENE/'POINT' ;
  18 : 
  19 : 'SI' ('<' ('DIMENSION' LL) 13) ;
  20 :   'ERREUR' 'La liste des caracteristiques doit avoir 13 valeurs' ;
  21 : 'FINSI' ;
  22 : 
  23 : DIMGEO = 'VALEUR' 'DIME' ;  MODGEO = 'VALEUR' 'MODE' ;
  24 : 'SI' ('NON' ('EXISTE' TYPEMOD)) ;
  25 :   TYPEMOD = MODGEO ;
  26 : 'FINSI' ;
  27 : 'SI' ('NON' ('EXISTE' TYPEELT)) ;
  28 :   TYPEELT = 'MOT' 'COQ3' ;
  29 : 'FINSI' ;
  30 : Mot_DPGE = 'TEXTE' ' ' ;  Pt_DPGE = 'TEXTE' ' ' ;
  31 : 
  32 : 'SI' (('NEG' TYPEMOD 'COEP') 'ET' ('NEG' TYPEMOD 'COMI') 'ET'
  33 :       ('NEG' TYPEMOD MODGEO)) ;
  34 :   'ERREUR' ('CHAINE' 'Mot cle ' TYPEMOD ' non reconnu') ;
  35 : 'FINSI' ;
  36 : 'SI' ('EGA' TYPEMOD 'COMI') ;
  37 :   'SI' ('NON' ('EXISTE' ('MOTS' 'COQ3' 'DKT ') TYPEELT)) ;
  38 :     'ERREUR' 'Element de COque MInce a support TRI3 inconnu' ;
  39 :   'FINSI' ;
  40 : 'FINSI' ;
  41 : 'SI' ('EGA' TYPEMOD 'COEP') ;
  42 :   'SI' ('NEG' DIMGEO 3) ;
  43 :     'ERREUR' 'Modelisation COEP incompatible avec la DIMEnsion' ;
  44 :   'FINSI' ;
  45 : 'FINSI' ;
  46 : 
  47 : 'SI' ('EGA' DIMGEO 1) ;
  48 :   'ERREUR' 'Cas DIME=1 non prevu actuellement' ;
  49 : 'FINSI' ;
  50 : 'SI' ('EGA' DIMGEO 3) ;
  51 :   'SI' ('EGA' TYPEMOD 'TRID') ;
  52 :     lm_geo = 'MOTS' 'TET4' 'PRI6' 'PYR5' 'CUB8'
  53 :                     'TE10' 'PR15' 'PY14' 'CU20' ;
  54 :     lm_elt = lm_geo ;
  55 :   'FINSI' ;
  56 :   'SI' ('EGA' TYPEMOD 'COMI') ;
  57 :     lm_geo = 'MOTS' 'TRI3' ;
  58 :     lm_elt = 'MOTS' TYPEELT ;
  59 :   'FINSI' ;
  60 :   'SI' ('EGA' TYPEMOD 'COEP') ;
  61 :     lm_geo = 'MOTS' 'TRI3' 'QUA4' 'TRI6' 'QUA8' ;
  62 :     lm_elt = 'MOTS' 'DST ' 'COQ4' 'COQ6' 'COQ8' ;
  63 :   'FINSI' ;
  64 : 'FINSI' ;
  65 : 'SI' ('EGA' DIMGEO 2) ;
  66 :   'SI' ('EGA' MODGEO 'PLANGENE') ;
  67 :     'SI' ('NON' ('EXISTE' PTGENE)) ;
  68 :       'ERREUR' 'Le point support de la DPGE n est pas defini' ;
  69 :     'FINSI' ;
  70 :     Mot_DPGE = 'MOT' 'DPGE' ;  Pt_DPGE = PTGENE ;
  71 :   'FINSI' ;
  72 :   lm_geo = 'MOTS' 'TRI3' 'QUA4' 'TRI6' 'QUA8' ;
  73 :   lm_elt = lm_geo ;
  74 :   'SI' ('EGA' TYPEMOD 'COMI') ;
  75 :     lm_geo = 'MOTS' 'SEG2' ;
  76 :     lm_elt = 'MOTS' 'COQ2' ;
  77 :   'FINSI' ;
  78 : 'FINSI' ;
  79 : 
  80 : lm_mai = 'ELEM' MAIL1 'TYPE' ;
  81 : 'REPETER' Boucle ('DIMENSION' lm_mai) ;
  82 :   str_z = 'EXTRAIRE' lm_mai &Boucle ;
  83 :   'SI' ('NON' ('EXISTE' lm_geo str_z)) ;
  84 :     'ERREUR' 'Element geometrique incompatible avec la modelisation' ;
  85 :   'FINSI' ;
  86 : 'FIN' Boucle ;
  87 : 
  88 : i_z = 0 ;
  89 : 'REPETER' Boucle ('DIMENSION' lm_geo) ;
  90 :   str_z = 'EXTRAIRE' lm_geo &Boucle ;
  91 :   'SI' ('EXISTE' lm_mai str_z) ;
  92 :     mail_z = MAIL1 'ELEM' str_z ;
  93 :     str_elt = 'EXTRAIRE' lm_elt &Boucle ;
  94 :     modl_z = 'MODL' mail_z 'MECANIQUE' 'ELASTIQUE' 'ORTHOTROPE'
  95 :                            str_elt 'CONS' 'MECA' Mot_DPGE Pt_DPGE ;
  96 :     i_z = i_z + 1 ;
  97 :     'SI' ('EGA' i_z 1) ;
  98 :       MODL1 = modl_z ;
  99 :     'SINON' ;
 100 :       MODL1 = MODL1 'ET' modl_z ;
 101 :     'FINSI' ;
 102 :   'FINSI' ;
 103 : 'FIN' Boucle ;
 104 : 
 105 : GEO = 'CHANGER' 'POI1' MAIL1 ;
 106 : NBPT1 = 'NBNO' GEO ;
 107 : 
 108 : L1 = 'EXTRAIRE' CH1 'COMP' ;
 109 : L2 = 'EXTRAIRE' CH2 'COMP' ;
 110 : 
 111 : 'SI' (DIMGEO 'EGA' 3) ;
 112 :   L11 = 'EXTRAIRE' 1 L1 ;
 113 :   L12 = 'EXTRAIRE' 2 L1 ;
 114 :   L13 = 'EXTRAIRE' 3 L1 ;
 115 :   L21 = 'EXTRAIRE' 1 L2 ;
 116 :   L22 = 'EXTRAIRE' 2 L2 ;
 117 :   L23 = 'EXTRAIRE' 3 L2 ;
 118 :   LCOMP_K  = 'MOTS' 'KX  ' 'KY  ' 'KZ  ' ;
 119 :   LCOMP_V1 = 'MOTS' 'V1X ' 'V1Y ' 'V1Z ' ;
 120 :   LCOMP_V2 = 'MOTS' 'V2X ' 'V2Y ' 'V2Z ' ;
 121 :   'REPETER' Boucle NBPT1 ;
 122 :     POINMAIL = GEO 'POINT' &Boucle ;
 123 :     POINCOUR = LIG1 'POINT' 'PROCHE' POINMAIL ;
 124 :     CH1MAIL = 'MANU' 'CHPO' POINMAIL 3
 125 :                             'KX  ' ('EXTRAIRE' CH1 L11 POINCOUR)
 126 :                             'KY  ' ('EXTRAIRE' CH1 L12 POINCOUR)
 127 :                             'KZ  ' ('EXTRAIRE' CH1 L13 POINCOUR)
 128 :                             'NATURE' 'DIFFUS' ;
 129 :     CH2MAIL = 'MANU' 'CHPO' POINMAIL 3
 130 :                             'KX  ' ('EXTRAIRE' CH2 L21 POINCOUR)
 131 :                             'KY  ' ('EXTRAIRE' CH2 L22 POINCOUR)
 132 :                             'KZ  ' ('EXTRAIRE' CH2 L23 POINCOUR)
 133 :                             'NATURE' 'DIFFUS' ;
 134 :     'SI' ('EGA' &Boucle 1) ;
 135 :       CH1GEO = CH1MAIL ;
 136 :       CH2GEO = CH2MAIL ;
 137 :     'SINON' ;
 138 :       CH1GEO = CH1GEO 'ET' CH1MAIL ;
 139 :       CH2GEO = CH2GEO 'ET' CH2MAIL ;
 140 :     'FINSI' ;
 141 :   'FIN' Boucle ;
 142 : 'FINSI' ;
 143 : 'SI' (DIMGEO 'EGA' 2) ;
 144 :   L11 = 'EXTRAIRE' 1 L1 ;
 145 :   L12 = 'EXTRAIRE' 2 L1 ;
 146 :   L21 = 'EXTRAIRE' 1 L2 ;
 147 :   L22 = 'EXTRAIRE' 2 L2 ;
 148 :   LCOMP_K  = 'MOTS' 'KX  ' 'KY  ' ;
 149 :   LCOMP_V1 = 'MOTS' 'V1X ' 'V1Y ' ;
 150 :   LCOMP_V2 = 'MOTS' 'V2X ' 'V2Y ' ;
 151 :   'REPETER' Boucle NBPT1 ;
 152 :     POINMAIL = GEO 'POINT' &Boucle ;
 153 :     POINCOUR = LIG1 'POINT' 'PROCHE' POINMAIL ;
 154 :     CH1MAIL = 'MANU' 'CHPO' POINMAIL 2
 155 :                             'KX  ' ('EXTRAIRE' CH1 L11 POINCOUR)
 156 :                             'KY  ' ('EXTRAIRE' CH1 L12 POINCOUR)
 157 :                             'NATURE' 'DIFFUS' ;
 158 :     CH2MAIL = 'MANU' 'CHPO' POINMAIL 2
 159 :                             'KX  ' ('EXTRAIRE' CH2 L21 POINCOUR)
 160 :                             'KY  ' ('EXTRAIRE' CH2 L22 POINCOUR)
 161 :                             'NATURE' 'DIFFUS' ;
 162 :     'SI' ('EGA' &Boucle 1) ;
 163 :       CH1GEO = CH1MAIL ;
 164 :       CH2GEO = CH2MAIL ;
 165 :     'SINON' ;
 166 :       CH1GEO = CH1GEO 'ET' CH1MAIL ;
 167 :       CH2GEO = CH2GEO 'ET' CH2MAIL ;
 168 :     'FINSI' ;
 169 :   'FIN' Boucle ;
 170 : 'FINSI' ;
 171 : 
 172 : CHP1 = (('COS' ALPHA)*CH1GEO) + (('SIN' ALPHA)*CH2GEO) ;
 173 : CHP2 = (((-1.*('SIN' ALPHA))*CH1GEO)) + (('COS' ALPHA)*CH2GEO) ;
 174 : 
 175 : CHPV1 = 'NOMC' LCOMP_K LCOMP_V1 CHP1 ;
 176 : CHMV1 = 'CHANGER' 'CHAM' CHPV1 MODL1 'RIGIDITE' 'CARACTERISTIQUES' ;
 177 : CHPV2 = 'NOMC' LCOMP_K LCOMP_V2 CHP2 ;
 178 : CHMV2 = 'CHANGER' 'CHAM' CHPV2 MODL1 'RIGIDITE' 'CARACTERISTIQUES' ;
 179 : 
 180 : YOUNG1S = 'EXTRAIRE' LL  1 ;
 181 : YOUNG2S = 'EXTRAIRE' LL  2 ;
 182 : YOUNG3S = 'EXTRAIRE' LL  3 ;
 183 : NU1S    = 'EXTRAIRE' LL  4 ;
 184 : NU2S    = 'EXTRAIRE' LL  5 ;
 185 : NU3S    = 'EXTRAIRE' LL  6 ;
 186 : CIS1S   = 'EXTRAIRE' LL  7 ;
 187 : CIS2S   = 'EXTRAIRE' LL  8 ;
 188 : CIS3S   = 'EXTRAIRE' LL  9 ;
 189 : ALPH1S  = 'EXTRAIRE' LL 10 ;
 190 : ALPH2S  = 'EXTRAIRE' LL 11 ;
 191 : ALPH3S  = 'EXTRAIRE' LL 12 ;
 192 : RHOS    = 'EXTRAIRE' LL 13 ;
 193 : 
 194 : YOUNG1 = 'MANU' 'CHML' MODL1 'YG1 ' YOUNG1S
 195 :                        'TYPE' 'CARACTERISTIQUES' 'RIGIDITE' ;
 196 : YOUNG2 = 'MANU' 'CHML' MODL1 'YG2 ' YOUNG2S
 197 :                        'TYPE' 'CARACTERISTIQUES' 'RIGIDITE' ;
 198 : YOUNG3 = 'MANU' 'CHML' MODL1 'YG3'  YOUNG3S
 199 :                        'TYPE' 'CARACTERISTIQUES' 'RIGIDITE' ;
 200 : NU1    = 'MANU' 'CHML' MODL1 'NU12' NU1S
 201 :                        'TYPE' 'CARACTERISTIQUES' 'RIGIDITE' ;
 202 : NU2    = 'MANU' 'CHML' MODL1 'NU23' NU2S
 203 :                        'TYPE' 'CARACTERISTIQUES' 'RIGIDITE' ;
 204 : NU3    = 'MANU' 'CHML' MODL1 'NU13' NU3S
 205 :                        'TYPE' 'CARACTERISTIQUES' 'RIGIDITE' ;
 206 : CIS1   = 'MANU' 'CHML' MODL1 'G12 ' CIS1S
 207 :                        'TYPE' 'CARACTERISTIQUES' 'RIGIDITE' ;
 208 : CIS2   = 'MANU' 'CHML' MODL1 'G23 ' CIS2S
 209 :                        'TYPE' 'CARACTERISTIQUES' 'RIGIDITE' ;
 210 : CIS3   = 'MANU' 'CHML' MODL1 'G13 ' CIS3S
 211 :                        'TYPE' 'CARACTERISTIQUES' 'RIGIDITE' ;
 212 : ALPH1  = 'MANU' 'CHML' MODL1 'ALP1' ALPH1S
 213 :                        'TYPE' 'CARACTERISTIQUES' 'RIGIDITE' ;
 214 : ALPH2  = 'MANU' 'CHML' MODL1 'ALP2' ALPH2S
 215 :                        'TYPE' 'CARACTERISTIQUES' 'RIGIDITE' ;
 216 : ALPH3  = 'MANU' 'CHML' MODL1 'ALP3' ALPH3S
 217 :                        'TYPE' 'CARACTERISTIQUES' 'RIGIDITE' ;
 218 : RHO    = 'MANU' 'CHML' MODL1 'RHO ' RHOS
 219 :                        'TYPE' 'CARACTERISTIQUES' 'RIGIDITE' ;
 220 : 
 221 : 'SI' ('EGA' TYPEMOD 'COMI') ;
 222 :   MAT11 = YOUNG1 'ET' YOUNG2 'ET' NU1 ;
 223 :   MAT22 = CIS1 ;
 224 :   'SI' ('EGA' DIMGEO 3) ;
 225 :     MAT33 = ('EXCO' ('MOTS' 'V1X ' 'V1Y ') CHMV1) 'ET' ;
 226 :     MAT33 = ('EXCO' ('MOTS' 'V2X ' 'V2Y ') CHMV2) ;
 227 :   'SINON' ;
 228 :     MAT33 = CHMV1 'ET' CHMV2 ;
 229 :   'FINSI' ;
 230 :   MAT44 = ALPH1 'ET' ALPH2 'ET' RHO ;
 231 : 'FINSI' ;
 232 : 'SI' ('EGA' TYPEMOD 'COEP') ;
 233 :   MAT11 = YOUNG1 'ET' YOUNG2 'ET' NU1 ;
 234 :   MAT22 = CIS1 'ET' CIS2 'ET' CIS3 ;
 235 :   'SI' ('EGA' DIMGEO 3) ;
 236 :     MAT33 = ('EXCO' ('MOTS' 'V1X ' 'V1Y ') CHMV1) 'ET' ;
 237 :     MAT33 = ('EXCO' ('MOTS' 'V2X ' 'V2Y ') CHMV2) ;
 238 :   'SINON' ;
 239 :     MAT33 = CHMV1 'ET' CHMV2 ;
 240 :   'FINSI' ;
 241 :   MAT44 = ALPH1 'ET' ALPH2 'ET' RHO ;
 242 : 'FINSI' ;
 243 : 'SI' ('EGA' TYPEMOD 'PLANCONT') ;
 244 :   MAT11 = YOUNG1 'ET' YOUNG2 'ET' YOUNG3 'ET' NU1 'ET' NU2 'ET' NU3 ;
 245 :   MAT22 = CIS1 ;
 246 :   MAT33 = CHMV1 'ET' CHMV2 ;
 247 :   MAT44 = ALPH1 'ET' ALPH2 'ET' RHO ;
 248 : 'FINSI' ;
 249 : 'SI' (('EGA' TYPEMOD 'PLANDEFO') 'OU' ('EGA' TYPEMOD 'PLANGENE') 'OU'
 250 :       ('EGA' TYPEMOD 'AXIS')) ;
 251 :   MAT11 = YOUNG1 'ET' YOUNG2 'ET' YOUNG3 'ET' NU1 'ET' NU2 'ET' NU3 ;
 252 :   MAT22 = CIS1 ;
 253 :   MAT33 = CHMV1 'ET' CHMV2 ;
 254 :   MAT44 = ALPH1 'ET' ALPH2 'ET' ALPH3 'ET' RHO ;
 255 : 'FINSI' ;
 256 : 'SI' (('EGA' TYPEMOD 'FOUR') 'OU' ('EGA' TYPEMOD 'TRID')) ;
 257 :   MAT11 = YOUNG1 'ET' YOUNG2 'ET' YOUNG3 'ET' NU1 'ET' NU2 'ET' NU3 ;
 258 :   MAT22 = CIS1 'ET' CIS2 'ET' CIS3 ;
 259 :   MAT33 = CHMV1 'ET' CHMV2 ;
 260 :   MAT44 = ALPH1 'ET' ALPH2 'ET' ALPH3 'ET' RHO ;
 261 : 'FINSI' ;
 262 : 
 263 : MATGEO = MAT11 'ET' MAT22 'ET' MAT33 'ET' MAT44 ;
 264 : 
 265 : 'FINPROC' MATGEO MODL1 ;
 266 : 
 267 :  

© Cast3M 2003 - All rights reserved.
Disclaimer