Télécharger la plaquette Cast3M


fr en

Télécharger deconv3D.procedur

Retour à la liste

Numérotation des lignes : 

  1. * DECONV3D  PROCEDUR  WANG      16/01/25    21:15:00     8800           
  2. 'DEBPROC' DECONV3D COUCHE*'TABLE' PC*'POINT' FOND_SOL*'MAILLAGE'
  3.           MOD_SOL*MMODEL MAT_SOL*MCHAML DIR*'MOT' GAMMAO*'EVOLUTIO'
  4.           F1*'FLOTTANT'F2*'FLOTTANT' FC/'FLOTTANT' P_GAMMA/'TABLE' ;
  5.  
  6. ************************************************************************
  7. *
  8. *                            D E C O N V 3 D
  9. *                           -----------------
  10. *
  11. * FONCTION:
  12. * ---------
  13. *
  14. * Cette procedure permet d'effectuer des calculs sismiques
  15. * d'interaction sol-structure (ISS) en 3D par la methode des
  16. * elements finis.
  17. *
  18. * PHRASE D'APPEL (EN GIBIANE):
  19. * ----------------------------
  20. *
  21. *   TABRESU =  DECONV3D COUCHE PC FOND_SOL MOD_SOL MAT_SOL DIR
  22. *                       GAMMAO F1 F2 (FC) (P_GAMMA)
  23. *
  24. * OPERANDES:
  25. * ----------
  26. *
  27. *   COUCHE 'TABLE' : table de donnees
  28. *     COUCHE.'BLK' : Si existe, deconvolution pour la methode de Bielak
  29. *     COUCHE.'BLK'.'EP_H' : epaisseur horizontale de la zone intermediare
  30. *                           verticale
  31. *     COUCHE.'BLK'.'EP_V' : epaisseur veriticle de la zone intermediare
  32. *                           horizontale
  33. *     COUCHE.I.'indice' : donnees pour la ieme couche du sol
  34. *
  35. *       indice 'FRONTIERE'       'MAILLAGE' : frontiere verticale
  36. *       indice 'MASSE_VOLUMIQUE' 'FLOTTANT' : masse volumique
  37. *       indice 'POISSON'         'FLOTTANT' : coefficient de Poisson
  38. *       indice 'YOUNG'           'FLOTTANT' : module d'Young
  39. *       indice 'AMORTISSEMENT    'FLOTTANT' : amortissement reduit
  40. *
  41. *   PC       'MAILLAGE' : point de reference situe au milieu du millage
  42. *                         de solou sur l'axe de symetrie si calcul sur
  43. *                         la moitie et le quart du systeme
  44. *   FOND_SOL 'MAILLAGE' : frontiere horizontale inferieure du sol
  45. *   MOD_SOL  'MMODEL'   : modele du sol
  46. *   MAT_SOL  'MCHAML'   : materiau du sol
  47. *   DIR      'MOT'      : direction de l'acceleration sismique GAMMAO
  48. *                        'UX' pour la direction horizontale X
  49. *                        'UY' pour la direction horizontale Y
  50. *                        'UZ' pour la direction verticale Z
  51. *   GAMMAO   'EVOLUTIO' : Signal sismique (accelerogramme) d'entree pour
  52. *                         la deconvolution, defini par defaut a la surface
  53. *                         libre du sol, sinon, renseigner la table P_GAMMA
  54. *   F1, F2   'FLOTTANT' : frequences sur lesquelles l'amortissement
  55. *                       : reduit est ajuste suivant le modele de
  56. *                         RAYLEIGH
  57. *   FC       'FLOTTANT' : frequence de coupure pour la deconvolution,
  58. *                         par defaut  FC = 50 Hz
  59. *   P_GAMMA  'TABLE'    : Description des accelerogrammes d'entree et de sorti
  60. *                         pour la deconvolution (facultatif)
  61. *     indice 'ENTREE'   :'TABLE' : description de la nature de l'accelerogramme
  62. *                         GAMMAO autre que sur la surface libre ou description
  63. *                         de l'onde incidente inclinée SH dans le cas d'un sol
  64. *                         demi-espace homogène.
  65. *                         (facultatif)
  66. *      sous-indice 'NATURE' : nature du point de controle :
  67. *                            'MOT' INSIDE : dans le sol
  68. *                            'MOT' OUTCROP : outcrop du bedrock
  69. *                            'MOT' SH : Onde incidente SH inclinée dans le cas 
  70. *                                  de la methode Bielak 
  71. *                                  DIR = UX implique onde SH dans le plan YZ, 
  72. *                                  DIR = UY implique onde SH dans le plan XZ,
  73. *
  74. *      sous-indice 'ANGLE' : Angle d'incidence de l'onde SH si 'NATURE' = SH
  75. *                            
  76. *      sous-indice 'CONTROLE' :'MAILLAGE': 
  77. *                               si 'NATURE' = INSIDE, point de controle sur la 
  78. *                               frontiere verticale 
  79. *                               si 'NATURE' = 'SH', point de controle à la 
  80. *                               surface du sol pour lequel on impose 
  81. *                               l'accélérogramme du champ libre GAMMAO
  82. *
  83. *     indice 'I'        :'MAILLAGE': Ieme (i = 1, 2, 3,...) points sur la
  84. *                         frontiere verticale pour lesquels on desire sortir
  85. *                         l'accelerogramme en champ libre (resultats de
  86. *                         deconvolution)
  87. * RESULTATS:
  88. * ----------
  89. *   TABRESU 'TABLE' :  table qui contient les resultats du calcul
  90. *
  91. *     indice 'CHAR'  'CHARGEMENT' : excitation sur la frontiere du sol
  92. *     indice 'AMOR'  'RIGIDITE'   : frontiere absorbante
  93. *     indice 'DEFO'  'EVOLUTIO'   : deformation maximale en fonction
  94. *                                   de la profondeur
  95. *     indice 'ACCE'  'TABLE'      : table qui contient les accelerations
  96. *                                   obtenues par la deconvolution
  97. *        sous-indice  I             acceleration du point I defini dans la
  98. *                                   table P_GAMMA
  99. *        sous-indice 'OUTCROP'      acceleration du outcrop de bedrock
  100. *
  101. *     indice 'PAS'   'FLOTTANT'
  102. *     indice 'FCDYN' 'FLOTTANT'   : pas de temps et frequence de coupure
  103. *                                   a utiliser pour le calcul de
  104. *                                   l'interaction sol-structure
  105. *                                   a l'aide de la procedure DYNAMIC
  106. *
  107. * AUTEUR, DATE DE CREATION ET DE MODIFICATION:
  108. * --------------------------------------------
  109. *
  110. *     F. WANG   20/05/2009 : Creation.
  111. *     F. WANG   07/10/2010 : Mise a zero des vitesses initiales
  112. *                            au sens moyen sur les 70 premiers pas
  113. *                            eliminant ainsi la rotation du bloc de sol,
  114. *                            Changement de format du chargement reduisant
  115. *                            largement le temps de calcul dynamique.
  116. *     F. WANG   01/02/2012 : Correction d'un bug MAT_S -> MAT_SOL,
  117. *                            Elimination des accents dans les
  118. *                            commentaires pour eviter des problemes
  119. *     F. WANG   20/02/2013 : Ajout des options concernant le signal
  120. *                            d'entree : INSIDE et OUTCROP
  121. *                            Ajout du signal de sortie outcrop de bedrock       
  122. *     F. WANG   01/02/2014 : Implantation de la methode de Bielak
  123. *                            permettant d'alleger le calcul
  124. *     A. FRAU   01/07/2015 : Correction Bug
  125. *
  126. *     F. WANG   01/10/2015 : Implantation d'une nouvelle option: Onde SH
  127. *                            inclinée dans l'option Bielak
  128. *                           
  129. ************************************************************************
  130. 'SAUT' 1 'LIGNE' ;
  131.  
  132.  SOL  = 'EXTR' MOD_SOL 'MAILLAGE';
  133.  
  134. 'SI' (EXIS COUCHE BLK) ;
  135.  BIELAK = VRAI ;
  136.  NB_COU = ('DIME' COUCHE) - 1 ;
  137. 'MESS' 'DECONVOLUTION POUR LA METHODE DE BIELAK' ;
  138. 'SINON' ;
  139.  BIELAK = FAUX ;
  140.  NB_COU = 'DIME' COUCHE ;
  141. 'FINSI' ;
  142.  
  143. *-----------DEBUT DE MODIFICATION 01/10/2015. WANG)--------------------
  144.  
  145.  SH = FAUX ;
  146.  
  147. 'SI' (EXIS P_GAMMA) ;
  148.   'SI' (EXIS P_GAMMA ENTREE) ;
  149.    'SI' (EGA P_GAMMA.ENTREE.NATURE 'SH') ;
  150.     ANGLE  = P_GAMMA.ENTREE.ANGLE ;
  151.     SH = VRAI ;
  152.       SI (EGA DIR 'UX') ;
  153.      'MESS' 'ONDE SH SE PROPAGEANT DANS LE PLAN YZ' ;
  154.      'MESS' 'ANGLE D"INCIDENCE = ' ANGLE 'DEGREES'  ;
  155.       FINSI ;
  156.       SI (EGA DIR 'UY') ;
  157.      'MESS' 'ONDE SH SE PROPAGEANT DANS LE PLAN XZ' ;
  158.      'MESS' 'ANGLE D"INCIDENCE =' ANGLE 'DEGREES';
  159.       FINSI ;
  160.       SI (NB_COU > 1);
  161.      'SAUT' 1 'LIGNE' ;
  162.      'MESS' '***** OPTION NON PREVUE POUR UN SOL MULTICOUCHE *****';
  163.      'MESS' '***** ARRET DE LA PROCEDURE DECONV3D *****' ;
  164.      'SAUT' 1 'LIGNE' ;
  165.       QUIT DECONV3D ;
  166.       FINSI ;
  167.       SI (NON BIELAK) ;
  168.      'SAUT' 1 'LIGNE' ;
  169.      'MESS' '***** OPTION PREVUE SEULEMENT DANS OPTION BIELAK *****';
  170.      'MESS' '***** ARRET DE LA PROCEDURE DECONV3D *****' ;
  171.      'SAUT' 1 'LIGNE' ;
  172.       QUIT DECONV3D ;
  173.       FINSI ;    
  174.     'FINSI' ;
  175.   'FINSI' ;
  176. 'FINSI' ;
  177.  
  178. *-----------FIN DE MODIFICATION 01/10/2015 F. WANG)--------------------
  179.  
  180. *   Verification du maillage de la frontiere du sol :
  181. *   Types d'elements de frontiere : QUA4 + TRI3 ou QUA8 + TRI6
  182.  
  183.  FRSOL  =  FOND_SOL;
  184.  
  185.  I = 1;
  186. 'REPE' BCL0 NB_COU;
  187.  FRSOL = FRSOL 'ET' COUCHE.I.'FRONTIERE';
  188.  I = I + 1;
  189. 'FIN' BCL0;
  190.  
  191.  LTYPF = 'ELEM' FRSOL 'TYPE';
  192.  
  193. 'SI' ((DIME LTYPF) > 2);
  194. 'SAUT' 2 'LIGNE' ;
  195. 'MESS' '***** ARRET DE LA PROCEDURE DECONV3D *****';
  196. 'MESS' '***** FRONTIERE : MAILLAGE INCORRECT *****';
  197. 'SAUT' 2 'LIGNE' ;
  198. 'QUIT' DECONV3D;
  199. 'FINSI';
  200.  
  201. 'SI' ((DIME LTYPF) 'EGA' 2);
  202.  TYPF1 = 'EXTR' LTYPF 1;
  203.  TYPF2 = 'EXTR' LTYPF 2;
  204.   'SI' ((('EGA' TYPF1 'QUA4') ET ('EGA' TYPF2 'TRI3')) OU
  205.         (('EGA' TYPF1 'TRI3') ET ('EGA' TYPF2 'QUA4')));
  206.    TY1 = 'CUB8';
  207.    TY2 = 'PRI6';
  208.   'SINON';
  209.     'SI' ((('EGA' TYPF1 'QUA8') ET ('EGA' TYPF2 'TRI6')) OU
  210.           (('EGA' TYPF1 'TRI6') ET ('EGA' TYPF2 'QUA8')));
  211.      TY1 = 'CU20';
  212.      TY2 = 'PR15';
  213.     'SINON';
  214.     'SAUT' 2 'LIGNE' ;
  215.     'MESS' '***** ARRET DE LA PROCEDURE DECONV3D     *****';
  216.     'MESS' '***** FRONTIERE : ELEMENTS INCOMPATIBLES *****';
  217.     'SAUT' 2 'LIGNE' ;
  218.     'QUIT' DECONV3D;
  219.     'FINSI';
  220.   'FINSI';
  221. 'FINSI';
  222.  
  223. 'SI' ((DIME LTYPF) 'EGA' 1);
  224.  TYPF1 = 'EXTR' LTYPF 1;
  225.   'SI' (('EGA' TYPF1 'QUA4') OU ('EGA' TYPF1 'TRI3'));
  226.    TY1 = 'CUB8';
  227.    TY2 = 'PRI6';
  228.   'SINON';
  229.     'SI' (('EGA' TYPF1 'QUA8') OU ('EGA' TYPF1 'TRI6'));
  230.      TY1 = 'CU20';
  231.      TY2 = 'PR15';
  232.     'SINON';
  233.     'SAUT' 2 'LIGNE' ;
  234.     'MESS' '***** ARRET DE LA PROCEDURE DECONV3D     *****';
  235.     'MESS' '***** FRONTIERE : ELEMENTS INCOMPATIBLES *****';
  236.     'SAUT' 2 'LIGNE' ;
  237.     'QUIT' DECONV3D;
  238.     'FINSI';
  239.   'FINSI';
  240. 'FINSI';
  241.  
  242. 'OPTI' 'ELEM' TY1 ;
  243.  
  244. 'SAUT' 5 'LIGNE'                                                    ;
  245. 'MESS' '**********************************************************' ;
  246. 'MESS' '*                                                        *' ;
  247. 'MESS' '*    CALCUL DU CHAMP SISMIQUE INCIDENT PAR LA METHODE    *' ;
  248. 'MESS' '*                   DE LA DECONVOLUTION                  *' ;
  249. 'MESS' '*                                                        *' ;
  250. 'MESS' '**********************************************************' ;
  251. 'SAUT' 2 'LIGNE'                                                    ;
  252. 'MESS' 'CARACTERISTIQUES MECANIQUES DU SOL'                         ;
  253.  
  254.  I = 0                                                              ;
  255. 'REPE' BOUCL0 NB_COU                                                ;
  256.  I = I + 1                                                          ;
  257.  E_I  = COUCHE.I.'YOUNG'                                            ;
  258.  NU_I = COUCHE.I.'POISSON'                                          ;
  259.  COUCHE.I.'MODULE_G'  = E_I / ( 2. * ( 1. + NU_I ) )                ;
  260.  CS_I = ( COUCHE.I.'MODULE_G' / COUCHE.I.'MASSE_VOLUMIQUE' ) ** .5  ;
  261.  CP_I = ((2.0 * (1.0 - NU_I))/(1.0 - (2.0 * NU_I))) ** 0.5          ;
  262.  CP_I = CP_I * CS_I                                                 ;
  263.  FR_I = COUCHE.I.'FRONTIERE'                                        ;
  264.  Z_I   = 'COOR' 3 FR_I                                              ;
  265.  Z_MAX = 'MAXI' Z_I                                                 ;
  266.  Z_MIN = 'MINI' Z_I                                                 ;
  267.  COUCHE.I.'EPAISSEUR' = Z_MAX - Z_MIN                               ;
  268.  
  269. 'SAUT' 'LIGNE'                                                      ;
  270. 'MESS' 'COUCHE NUMERO ' I                                           ;
  271. 'SAUT' 'LIGNE'                                                      ;
  272. 'MESS' 'EPAISSEUR DE LA COUCHE   H  ' COUCHE.I.'EPAISSEUR'          ;
  273. 'MESS' 'MODULE D"YOUNG           E  ' COUCHE.I.'YOUNG'              ;
  274. 'MESS' 'MASSE VOLUMIQUE          RHO' COUCHE.I.'MASSE_VOLUMIQUE'    ;
  275. 'MESS' 'COEFFICIENT DE POISSON   NU ' COUCHE.I.'POISSON'            ;
  276. 'MESS' 'VITESSE DES ONDES S      CS ' CS_I                          ;
  277. 'MESS' 'VITESSE DES ONDES P      CP ' CP_I                          ;
  278. 'MESS' 'AMORTISSEMENT EQUIVALENT KSI' COUCHE.I.'AMORTISSEMENT'      ;
  279. 'FIN' BOUCL0                                                        ;
  280.  
  281.  
  282.  
  283. ************* Construction d'une couche de sol enveloppe ***************
  284. *
  285. *       La couche enveloppe est construite a l'exterieur de la
  286. *       frontiere du sol (en utilisant le point de reference PC)
  287. *
  288. ************************************************************************
  289.  
  290. 'SAUT' 3 'LIGNE'                                            ;
  291. 'MESS' 'CONSTRUCTION D"UNE ENVELOPPE DE SOL'                ;
  292.  A_I   = 'MESU' FR_I 'SURF'                                 ;
  293.  NE_I  = 'NBEL' FR_I                                        ;
  294.  DR    = (A_I/NE_I) ** 0.5                                  ;
  295.     'SI' BIELAK                                             ;
  296.     'SI' ('EXIS' COUCHE.BLK EP_H)                           ;
  297.      DR = COUCHE.BLK.EP_H                                   ;
  298.     'FINSI'                                                 ;
  299.     'FINSI'                                                 ;
  300.  XMAX = 'MAXI' 'ABS' (('COOR' 1 FR_I) - ('COOR' 1 PC))      ;
  301.  YMAX = 'MAXI' 'ABS' (('COOR' 2 FR_I) - ('COOR' 2 PC))      ;
  302.  R_FR = ((XMAX**2) + (YMAX**2)) ** 0.5                      ;
  303.  RAPP = (R_FR + DR) / R_FR                                  ;
  304.  
  305. 'TITR' 'ENVELOPPE DE SOL POUR LA DECONVOLUTION (EN ROSE)'   ;
  306.  
  307.  COLONN = 'TABLE'  ;
  308.  COLMOD = 'TABLE'  ;
  309.  COLMAT = 'TABLE'  ;
  310.  COLMAS = 'TABLE'  ;
  311.  COLRIG = 'TABLE'  ;
  312.  COLBLO = 'TABLE'  ;
  313.  COLAMO = 'TABLE'  ;
  314.  
  315. 'SI' ( 'NON' ( 'EXIS' F1 ) )                    ;
  316.  F1    = 1.0                                    ;
  317.  F2    = 25.0                                   ;
  318. 'FINSI'                                         ;
  319. 'MESS' 'AMORTISSEMENT DE RAYLEIGH POUR LE SOL'  ;
  320. 'MESS' 'F1 = ' F1 'Hz      F2 = ' F2 'Hz'       ;
  321.  ALPHA =  PI * 4 * F1 * F2 / ( F1 + F2 )        ;
  322.  BETA  =  1 / ( PI * ( F1 + F2 ) )              ;
  323.  
  324.   I = 0                                                             ;
  325.  'REPE' BOUCL1 NB_COU                                               ;
  326.   I = I + 1                                                         ;
  327.   KSI      =  COUCHE.I.'AMORTISSEMENT'                              ;
  328.   FR_I     =  COUCHE.I.'FRONTIERE'                                  ;
  329.   FR_I2    = (FR_I 'HOMO' RAPP PC) 'AFFI' (1./RAPP)
  330.               PC ('PLUS' PC (0 0 -1))                               ;
  331.   COLONN.I =  FR_I 'VOLU' 1 FR_I2                                   ;
  332.   COLMOD.I = 'MODE' ( COLONN.I ) 'MECANIQUE' 'ELASTIQUE'
  333.              'ISOTROPE' TY1 TY2                                     ;
  334.   COLMAT.I = 'MATE' ( COLMOD.I ) 'YOUNG' ( COUCHE.I.'YOUNG' )
  335.              'NU'   ( COUCHE.I.'POISSON' )
  336.              'RHO'  ( COUCHE.I.'MASSE_VOLUMIQUE' )                  ;
  337.   COLMAS.I = 'MASS' ( COLMOD.I ) ( COLMAT.I )                       ;
  338.   COLRIG.I = 'RIGI' ( COLMOD.I ) ( COLMAT.I )                       ;
  339.  'SI' (('EGA' DIR 'UX') OU ('EGA' DIR 'UY'))                        ;
  340.   COLBLO.I = 'BLOQ' 'UZ' COLONN.I                                   ;
  341.  'SINON'                                                            ;
  342.  'SI' ('EGA' DIR 'UZ')                                              ;
  343.   COLBLO.I = 'BLOQ' 'UX' 'UY' COLONN.I                              ;
  344.  'FINSI'                                                            ;
  345.  'FINSI'                                                            ;
  346.   COLAMO.I =  KSI * (( ALPHA * ( COLMAS.I ) ) 'ET'
  347.              ( BETA * ( COLRIG.I 'ET' COLBLO.I )))                  ;
  348.  
  349.        'SI' ( I 'EGA' 1 )           ;
  350.         FRON = COUCHE.I.'FRONTIERE' ;
  351.  
  352.         COL1 = COLONN.I             ;
  353.         MOD1 = COLMOD.I             ;
  354.         MAT1 = COLMAT.I             ;
  355.         MAS1 = COLMAS.I             ;
  356.         RIG1 = COLRIG.I             ;
  357.         BLO1 = COLBLO.I             ;
  358.         AMO1 = COLAMO.I             ;
  359.   'SINON'                                    ;
  360.    FRON = FRON 'ET' ( COUCHE.I.'FRONTIERE' ) ;
  361.    COL1 = COL1 'ET' ( COLONN.I )   ;
  362.    MOD1 = MOD1 'ET' ( COLMOD.I )   ;
  363.    MAT1 = MAT1 'ET' ( COLMAT.I )   ;
  364.    MAS1 = MAS1 'ET' ( COLMAS.I )   ;
  365.    RIG1 = RIG1 'ET' ( COLRIG.I )   ;
  366.    BLO1 = BLO1 'ET' ( COLBLO.I )   ;
  367.    AMO1 = AMO1 'ET' ( COLAMO.I )   ;
  368.   'FINSI'                          ;
  369.   'FIN' BOUCL1                     ;
  370.  
  371. 'ELIM' COL1 0.01                       ;
  372.  COL1 = COL1 'COUL' 'ROSE'             ;
  373.  RIGT = RIG1 'ET' BLO1                 ;
  374. 'SI' ( 'EXIS' SOL )                    ;
  375. 'TRAC' ( SOL 'ET' COL1 ) 'FACE' 'NCLK' ;
  376. 'SINON'                                ;
  377. 'TRAC'  COL1  'FACE' 'NCLK'            ;
  378. 'FINSI'                                ;
  379.  
  380. *   Amortisseurs a la base de l'enveloppe
  381.  
  382. 'SAUT' 'LIGNE'                                                  ;
  383. 'MESS' 'AMORTISSEURS A LA BASE DE L ENVELOPPE'                  ;
  384.  
  385.  ZBA = 'MINI' ('COOR' 3 COLONN.NB_COU)                          ;
  386.  EBA = 'ENVE' COLONN.NB_COU                                     ;
  387.  PBA =  EBA 'POIN' 'PLAN' (0 0 ZBA) (1 0 ZBA)
  388.         (0 1 ZBA) 0.001                                         ;
  389.  SBAS = EBA 'ELEM' 'APPUYE' 'STRICTEMENT' PBA                   ;
  390.  BAS_COL = SBAS                                                 ;
  391.  AMO_BAS = AMOR COLMOD.NB_COU SBAS COLMAT.NB_COU                ;
  392.  AMOT    = AMO_BAS 'ET' AMO1                                    ;
  393.  
  394. *  Point surface libre et point base du sol
  395.  
  396.  ZSUR  = 'MAXI' ('COOR' 3 (COUCHE. 1 . 'FRONTIERE'))            ;
  397.  PSUR  = (COUCHE. 1 . 'FRONTIERE') 'POIN' 'PLAN' (0 0 ZSUR)
  398.          (1 0 ZSUR) (0 1 ZSUR) 0.001                            ;
  399.  P_SUR =  PSUR 'POINT' 'INITIAL' ;                              ;
  400.  
  401.  ZBAS  = 'MINI' ('COOR' 3 (COUCHE.NB_COU.'FRONTIERE'))          ;
  402.  PBAS  = (COUCHE.NB_COU.'FRONTIERE') 'POIN' 'PLAN' (0 0 ZBAS)
  403.          (1 0 ZBAS) (0 1 ZBAS) 0.001                            ;
  404.  P_BAS =  PBAS 'POINT' 'INITIAL'                                ;
  405.  
  406. *   Renumerotation des noeuds sur la frontiere :  TASSER "FRON"
  407.  
  409.  NB_NOEU = 'NBNO' FRON                                          ;
  410. 'MESS' 'NOMBRE DE NOEUDS SUR LA FRONTIERE :' NB_NOEU            ;
  411.  
  412. *   Acces aux noeuds sur la frontiere
  413.  
  414.  NUM_SUR  = 'NOEU' P_SUR  ;
  415.  NUM_BAS  = 'NOEU' P_BAS  ;
  416.  
  417.  POI_FRO = 'TABLE'         ;
  418.  I = 0                     ;
  419. 'REPE' BOUCL2 NB_NOEU      ;
  420.  I = I + 1                 ;
  421.  POI_FRO.I = 'NOEU' I      ;
  422. 'FIN'  BOUCL2              ;
  423.  
  424.  
  425. *******************  Reponse dynamique de L'enveloppe  *****************
  426. *
  427. *  Pour eviter des problemes de precision lies a l'interpolation du
  428. *  chargement produit par cette procedure, le pas de temps utilise ici
  429. *  pour la deconvolution sera egalement utilise pour le calcul d'ISS.
  430. *
  431. *  Le pas de temps est egale a 1/2, 1/4, 1/8, ... de celui de GAMMAO
  432. *  pour eviter des erreurs d'interpolation sur le signal d'entree. Il
  433. *  doit etre defini de sorte que la frequence de coupure utilisee dans
  434. *  la procedure DYNAMIC est au moins 4 fois FC.
  435. *
  436. *  Duree de calcul
  437. *
  438. *  sol monocouche  : 3  aller-retours onde S
  439. *  sol multicouche : 5  aller-retours onde S (prise en compte des
  440. *                       reflections sur les interfaces des couches)
  441. *
  442. ************************************************************************
  443.  
  444. *      Frequence de coupure  FC
  445.  
  446. 'SI' ( 'NON' ( 'EXIS' FC ) )               ;
  447.  FC = 50.0                                 ;
  448. 'FINSI'                                    ;
  449. 'SAUT' 'LIGNE'                             ;
  450. 'MESS' 'FREQUENCE DE COUPURE = ' FC 'Hz'   ;
  451.  
  452. 'SAUT' 3 'LIGNE'                                               ;
  453. 'MESS' 'REPONSE DYNAMIQUE DE L ENVELOPPE'                      ;
  454. 'SAUT' 'LIGNE'                                                 ;
  455.  T_GAMO = 'EXTR' GAMMAO 'ABSC'                                 ;
  456.  A_GAMO = 'EXTR' GAMMAO 'ORDO'                                 ;
  457.  NP     = 'DIME' T_GAMO                                        ;
  458.  TEM1   = 'EXTR' T_GAMO 1                                      ;
  459.  TEM2   = 'EXTR' T_GAMO 2                                      ;
  460.  TEMNP  = 'EXTR' T_GAMO NP                                     ;
  461.  PAS_TEM  = TEM2 - TEM1                                        ;
  462.  PAS_2FC  = 1.0 / ( 8 * FC )                                   ;
  463.  
  464.  IDIV = 0                                                      ;
  465. 'REPE' BB1                                                     ;
  466.  ERR = (( PAS_TEM / ( 2**IDIV )) - PAS_2FC ) / PAS_2FC         ;
  467. 'SI' ( ERR < 0.01 )                                            ;
  468. 'QUIT' BB1                                                     ;
  469. 'FINSI'                                                        ;
  470.  IDIV = IDIV + 1                                               ;
  471. 'FIN' BB1                                                      ;
  472.  PAS_TEM = PAS_TEM / (2**IDIV)                                 ;
  473.  
  474. *----------------------------------------------------------------------
  475. *  Aller a la partie (Bielak et SH) (modif du 01/10/2015 par F. WANG)
  476. *
  477. SI (NON (ET BIELAK SH)) ;
  478. *
  479. *----------------------------------------------------------------------   
  480.  
  481. 'SI' ( NB_COU 'EGA' 1 )                                      ;
  482.  G = ( COUCHE. 1 . 'MODULE_G' )                              ;
  483.  CS = ( G / ( COUCHE. 1 . 'MASSE_VOLUMIQUE' )  ) ** 0.5      ;
  484.  DUREE  = 6.0 * ( COUCHE. 1 . 'EPAISSEUR' ) / CS             ;
  485.  
  486. 'SINON'                                                      ;
  487.  DUREE = 0.0                                                 ;
  488.  I = 0                                                       ;
  489. 'REPE' BOUCD NB_COU                                          ;
  490.  I = I + 1                                                   ;
  491.  G = ( COUCHE.I.'MODULE_G' )                                 ;
  492.  CS = ( G / ( COUCHE.I.'MASSE_VOLUMIQUE' )  ) ** 0.5         ;
  493.  DUREE = DUREE + ( 10.0 * ( COUCHE.I.'EPAISSEUR' ) / CS )    ;
  494. 'FIN' BOUCD                                                  ;
  495. 'FINSI'                                                      ;
  496.  
  497. *  Impulsion de force horizontale a la base de l'enveloppe
  498. *  Sa transforme de Fourier doit etre non nulle dans
  499. *  l'intervalle (0, FC). On a DT = 7 *  PAS_TEM
  500.  
  501. 'SAUT'  'LIGNE'                                             ;
  502. 'MESS'  'FORCE APPLIQUEE A LA BASE DE L ENVELOPPE'          ;
  503. 'TITRE' 'FORCE APPLIQUEE A LA BASE DE L ENVELOPPE'          ;
  504.  DT       = 7 * PAS_TEM                                     ;
  505.  NT       = 300                                             ;
  506.  L_TIM    = 'PROG' 0.0 'PAS' DT 'NPAS' NT                   ;
  507.  L_FOR    = 'PROG' 0.  1.E5   0.  -1.E5   297 * 0.          ;
  508.  EV_FOR   = 'EVOL' 'MANU' 'TEMPS(S)' L_TIM 'FORCE(N)' L_FOR ;
  509.  
  510.  DUREE  =  DUREE + ( 4. * DT )                              ;
  511. 'SI' ( DUREE < ( 8. * DT ) )                                ;
  512.  DUREE  =  8. * DT                                          ;
  513. 'FINSI'                                                     ;
  514.  NB_PAS = 'ENTI' ( DUREE / PAS_TEM )                        ;
  515.  DUREE  =  NB_PAS * PAS_TEM                                 ;
  516. 'DESS' EV_FOR 'XBOR' 0.0 DUREE 'NCLK'                       ;
  517. 'SAUT' 'LIGNE'                                              ;
  518. 'MESS' 'CALCUL DYNAMIQUE POUR LA DECONVOLUTION'             ;
  519. 'MESS' 'PAS DE CALCUL = ' PAS_TEM '      DUREE = ' DUREE    ;
  520.  
  521. 'SI' ( 'EGA' DIR 'UX' )                          ;
  522.  FOR_P = 'PRES' 'MASS' COLMOD.NB_COU 1.0 SBAS    ;
  523.  FOR_BAS = 'EXCO' 'FZ' FOR_P 'FX'                ;
  524. 'SINON'                                          ;
  525.  
  526. 'SI' ( 'EGA' DIR 'UY' )                          ;
  527.  FOR_P = 'PRES' 'MASS' COLMOD.NB_COU 1.0 SBAS    ;
  528.  FOR_BAS = 'EXCO' 'FZ' FOR_P 'FY'                ;
  529. 'SINON'                                          ;
  530.  
  531. 'SI' ( 'EGA' DIR 'UZ' )                          ;
  532.  FOR_BAS = 'PRES' 'MASS' COLMOD.NB_COU 1.0 SBAS  ;
  533. 'FINSI'                                          ;
  534. 'FINSI'                                          ;
  535. 'FINSI'                                          ;
  536.  
  537.  CHA_BAS = 'CHAR' FOR_BAS EV_FOR                        ;
  538. 'LIST' FOR_BAS                                          ;
  539.  
  540. CH_DEPI = 'MANU' 'CHPO' 3 COL1  'UX' 0. 'UY' 0. 'UZ' 0. ;
  541. CH_VITI = 'MANU' 'CHPO' 3 COL1  'UX' 0. 'UY' 0. 'UZ' 0. ;
  542.  
  543. LINST ='PROG' TEM1 'PAS' PAS_TEM 'NPAS' (NB_PAS - 1)    ;
  544. FCDYN = 1.0 / (4 * PAS_TEM)                             ;
  545.  
  546. TAB_DYN        ='TABLE'  ;
  547. TAB_DYN.'DEBU' = TEM1    ;
  548. TAB_DYN.'VITE' = CH_VITI ;
  549. TAB_DYN.'DEPL' = CH_DEPI ;
  550. TAB_DYN.'CHAR' = CHA_BAS ;
  551. TAB_DYN.'RIGI' = RIGT    ;
  552. TAB_DYN.'MASS' = MAS1    ;
  553. TAB_DYN.'AMOR' = AMOT    ;
  554. TAB_DYN.'INST' = LINST   ;
  555. TAB_DYN.'FREQ' = FCDYN   ;
  556.  
  557. 'TEMP'                                           ;
  558.  RESU = DYNAMIC TAB_DYN                          ;
  559. 'SAUT' 'LIGNE'                                   ;
  560. 'MESS' 'FIN DE CALCUL DYNAMIQUE SUR L ENVELOPPE' ;
  561. 'MENA'                                           ;
  563. 'TEMP'                                           ;
  564.  
  565.  
  566. *   Deplacements, vitesses, deformations et reactions
  567.  
  568. 'SI' ( 'EGA' DIR 'UX' )         ;
  569.  UMOT = 'MOT' 'UX'              ;
  570.  EMOT = 'MOT' 'GAXZ'            ;
  571.  FMOT = 'MOT' 'FZ'              ;
  572. 'SINON'                         ;
  573.  
  574. 'SI' ( 'EGA' DIR 'UY' )         ;
  575.  UMOT = 'MOT' 'UY'              ;
  576.  EMOT = 'MOT' 'GAYZ'            ;
  577.  FMOT = 'MOT' 'FZ'              ;
  578. 'SINON'                         ;
  579.  
  580. 'SI' ( 'EGA' DIR 'UZ' )         ;
  581.  UMOT  = 'MOT' 'UZ'             ;
  582.  EMOT  = 'MOT' 'EPZZ'           ;
  583.  FMOT1 = 'MOT' 'FX'             ;
  584.  FMOT2 = 'MOT' 'FY'             ;
  585. 'FINSI'                         ;
  586. 'FINSI'                         ;
  587. 'FINSI'                         ;
  588.  
  589. 'SAUT' 3 'LIGNE'                                                    ;
  590. 'MESS' 'FONCTION DE TRANSFERT POUR LES DEPLACEMENTS, LES VITESSES,' ;
  591. 'MESS' 'LES DEFORMATIONS DE CISAILLEMENT ET LES REACTIONS '         ;
  592.  
  593. U_FRON  = 'TABLE' ;
  594. V_FRON  = 'TABLE' ;
  595. E_FRON  = 'TABLE' ;
  596. R_FRON  = 'TABLE' ;
  597. R1_FRON = 'TABLE' ;
  598. R2_FRON = 'TABLE' ;
  599.  
  600.   J = 0                  ;
  601.  'REPE' BOUCL3 NB_NOEU   ;
  602.   J = J + 1              ;
  603.   U_FRON.J  = 'PROG' 0.0 ;
  604.   V_FRON.J  = 'PROG' 0.0 ;
  605.   E_FRON.J  = 'PROG' 0.0 ;
  606.   R_FRON.J  = 'PROG' 0.0 ;
  607.   R1_FRON.J = 'PROG' 0.0 ;
  608.   R2_FRON.J = 'PROG' 0.0 ;
  609.  'FIN' BOUCL3            ;
  610.  
  611.  L_TEM = 'PROG' 0.0                                   ;
  612.  I_TEM = 1                                            ;
  613. 'REPE' BOUCL4  ( NB_PAS - 1 )                         ;
  614.  I_TEM = I_TEM + 1                                    ;
  615.  L_TEM = L_TEM 'ET' ( 'PROG' ( RESU.I_TEM.'TEMP' ) )  ;
  616.  CHD_I = RESU.I_TEM.'DEPL'                            ;
  617.  CHV_I = RESU.I_TEM.'VITE'                            ;
  618.  CHE_I ='CHAN' 'CHPO' MOD1 ( 'EPSI' MOD1 CHD_I MAT1 ) ;
  619.  CHR_I ='REAC' BLO1 CHD_I                             ;
  620.  
  621. 'SI' (('EGA' DIR 'UX') OU ('EGA' DIR 'UY')) ;
  622.  J = 0                                      ;
  623. 'REPE' BOUCL4_1 NB_NOEU                     ;
  624.  J = J + 1                                  ;
  625.  U_FRON.J = U_FRON.J 'ET' ('PROG' ('EXTR' CHD_I POI_FRO.J UMOT)) ;
  626.  V_FRON.J = V_FRON.J 'ET' ('PROG' ('EXTR' CHV_I POI_FRO.J UMOT)) ;
  627.  E_FRON.J = E_FRON.J 'ET' ('PROG' ('EXTR' CHE_I POI_FRO.J EMOT)) ;
  628.  R_FRON.J = R_FRON.J 'ET' ('PROG' ('EXTR' CHR_I POI_FRO.J FMOT)) ;
  629. 'FIN' BOUCL4_1           ;
  630. 'SINON'                  ;
  631.  
  632. 'SI' ('EGA' DIR 'UZ')    ;
  633.  J = 0                   ;
  634. 'REPE' BOUCL4_2 NB_NOEU  ;
  635.  J = J + 1               ;
  636.  U_FRON.J  = U_FRON.J  'ET' ('PROG' ('EXTR' CHD_I POI_FRO.J UMOT))  ;
  637.  V_FRON.J  = V_FRON.J  'ET' ('PROG' ('EXTR' CHV_I POI_FRO.J UMOT))  ;
  638.  E_FRON.J  = E_FRON.J  'ET' ('PROG' ('EXTR' CHE_I POI_FRO.J EMOT))  ;
  639.  R1_FRON.J = R1_FRON.J 'ET' ('PROG' ('EXTR' CHR_I POI_FRO.J FMOT1)) ;
  640.  R2_FRON.J = R2_FRON.J 'ET' ('PROG' ('EXTR' CHR_I POI_FRO.J FMOT2)) ;
  641. 'FIN' BOUCL4_2 ;
  642.  
  643. 'FINSI'        ;
  644. 'FINSI'        ;
  645. 'FIN' BOUCL4   ;
  646.  
  647. *      TFR de l'acceleration en surface GAMMAO
  648.  
  649. 'TITRE' 'TRANSFORMATION DE FOURIER DE GAMMAO' ;
  650.  NP = 'DIME' ( 'EXTR' GAMMAO 'ABSC' )         ;
  651.  XX = ( LOG NP ) / ( LOG 2 )                  ;
  652.  NN = 'ENTI' XX                               ;
  653. 'SI' ( XX > ( 1.0 * NN ) )                    ;
  654.  NN = NN + 1                                  ;
  655. 'FINSI'                                       ;
  656.  
  657.  TF_GAMO = 'TFR'  NN GAMMAO 'MOPH'            ;
  658.  FR_GAMO = 'EXTR' TF_GAMO   'ABSC'            ;
  659.  MO_GAMO = 'EXTR' TF_GAMO 'ORDO' 1            ;
  660.  PH_GAMO = 'EXTR' TF_GAMO 'ORDO' 2            ;
  661.  NFREQ   = 'DIME' FR_GAMO                     ;
  662.  DFREQ   = 'EXTR' FR_GAMO 2                   ;
  663.  FREQN   = 'EXTR' FR_GAMO NFREQ               ;
  664.  MODUN   = 'EXTR' MO_GAMO NFREQ               ;
  665.  
  666. 'SI' ( IDIV NEG 0 )                           ;
  667.  NN = NN + IDIV                               ;
  668. 'REMP' MO_GAMO NFREQ ( MODUN / 2.0 )          ;
  669.  
  670.  NFREQ1  = (2**(NN - 1)) + 1                                ;
  671.  FR_GAMO = 'PROG' 0.0 'PAS' DFREQ 'NPAS' ( NFREQ1 - 1 )     ;
  672. LIST ( DIME FR_GAMO ) ;
  673.  MO_GAMO = MO_GAMO 'ET' ('PROG' (NFREQ1 - NFREQ ) * 0.0)    ;
  674. LIST ( DIME MO_GAMO ) ;
  675.  PH_GAMO = PH_GAMO 'ET' ('PROG' (NFREQ1 - NFREQ ) * 0.0)    ;
  676. LIST ( DIME PH_GAMO ) ;
  677.  TF_GAMO = 'EVOL' 'COMP' 'MOPH' 'FREQ(HZ)' FR_GAMO 'MODULE'
  678.             MO_GAMO 'PHASE' PH_GAMO                         ;
  679. 'FINSI' ;
  680.  
  681. *-------------------- Debut de modif 1 du 20/02/2013 -------------------
  682.  
  683. *      TFR du deplacement de controle
  684.  
  685. 'TITRE' 'TFR DU DEPLACEMENT DE CONTROLE'                ;
  686.  OMEGA2 = ( FR_GAMO * 2. * PI ) ** 2                    ;
  687. 'REMP' OMEGA2 1 1.0                                     ;
  688.  MOD    = MO_GAMO / ( OMEGA2 * ( -1 ) )                 ;
  689. 'REMP' MOD 1 0.0                                        ;
  690.  TF_DEPLO = 'EVOL' 'COMP' 'MOPH' 'FREQ(HZ)' FR_GAMO
  691.             'MODULE' MOD 'PHASE' PH_GAMO                ;
  692.  
  693. *  TFR du deplacement du outcrop de bedrock sous l'impulsion
  694. *  Calcul de vitesse VOUT à partir de l'impulsion à la base F et
  695. *  l'amortissement à la base C (VOUT = F/C, sur une unite de surface)
  696.  
  697.  G  = ( COUCHE.NB_COU.'MODULE_G' )                      ;
  698.  RO = ( COUCHE.NB_COU.'MASSE_VOLUMIQUE' )               ;
  699.  NU = ( COUCHE.NB_COU.'POISSON')                        ;
  700.  CS = ( G / RO ) ** 0.5                                 ;
  701.  CP = ((2.0 * (1.0 - NU))/(1.0 - (2.0 * NU))) ** 0.5    ;
  702.  CP = CP * CS                                           ;
  703.  
  704. 'SI' (('EGA' DIR 'UX') OU ('EGA' DIR 'UY'))             ;
  705.  VOUT =  EV_FOR / (RO * CS)                             ;
  706. 'SINON'                                                 ;
  707. 'SI'( 'EGA' DIR 'UZ')                                   ;
  708.  VOUT = EV_FOR / (RO * CP)                              ;
  709. 'FINSI'                                                 ;
  710. 'FINSI'                                                 ;
  711. 'TITR' 'VITESSE OUTCROP SOUS L IMPULSION'               ;
  712.  VOUT = 'EVOL' 'MANU' L_TEM ('IPOL' L_TEM VOUT)         ;
  713.  TF_VOUT = 'TFR' NN  VOUT 'REIM'                        ;
  714.  
  715. 'TITR' 'TFR DEPLACEMENT OUTCROP SOUS L IMPULSION'       ;
  716.  LF_VOUT = EXTR TF_VOUT ABSC                            ;
  717.  IOMEGA = 2.*PI*LF_VOUT                                 ;
  718. 'REMP' IOMEGA 1 1.0                                     ;
  719.  EVIOMEGA = EVOL COMP 'REIM' LF_VOUT (LF_VOUT*0) IOMEGA ;
  720.  TF_UOUT = TF_VOUT / EVIOMEGA                           ;
  721.  RIMP TF_UOUT                                           ;
  722.  
  723. *     Fonction de transfert
  724.  
  725. H_DEPL  = 'TABLE' ;
  726. H_VITE  = 'TABLE' ;
  727. H_EPSI  = 'TABLE' ;
  728. H_REAC  = 'TABLE' ;
  729. H1_REAC = 'TABLE' ;
  730. H2_REAC = 'TABLE' ;
  731.  
  732. * Point de controle en surface libre (option par defaut)
  733.  
  734. 'TITRE'  'DEPLACEMENT EN SURFACE'                           ;
  735.  UO    = 'EVOL' 'MANU' 'TEMPS'  L_TEM 'DEPL' U_FRON.NUM_SUR ;
  736.  TF_UO = 'TFR' NN UO 'MOPH'                                 ;
  737.  
  738. * Remplacer le debut zero de TF_UOUT par celui de TF_UO
  739.  
  740. UOM  = EXTR (EXTR TF_UO  ORDO 1) 1 ;
  741. UOP  = EXTR (EXTR TF_UO  ORDO 2) 1 ;
  742. LF_OUT = EXTR TF_UOUT ABSC   ;
  743. LM_OUT = EXTR TF_UOUT ORDO 1 ;
  744. LP_OUT = EXTR TF_UOUT ORDO 2 ;
  745. REMP LM_OUT 1 UOM ;
  746. REMP LP_OUT 1 UOP ;
  747. TF_UOUT = 'EVOL' 'COMP' 'MOPH' 'FREQ(HZ)' LF_OUT
  748.           'MODULE' LM_OUT 'PHASE' LP_OUT  ;
  749.  
  750. 'SI' ('EXIS' P_GAMMA ) ;
  751. 'SI' ('EXIS' P_GAMMA 'ENTREE') ;
  752.  
  753. * Option point de contole "INSIDE"
  754.  
  755.    'SI' (EGA P_GAMMA.ENTREE.NATURE 'INSIDE');
  756.     NUM_CTR = 'NOEU' P_GAMMA.ENTREE.CONTROLE;
  757.     UO    = 'EVOL' 'MANU' 'TEMPS'  L_TEM 'DEPL' U_FRON.NUM_CTR ;
  758.     TF_UO = 'TFR' NN UO 'MOPH' ;
  759.    'FINSI' ;
  760.  
  761. * Option point de controle "OUTCROP" de bedrock
  762.  
  763.    'SI' (EGA P_GAMMA.ENTREE.NATURE 'OUTCROP') ;
  764.     TF_UO = TF_UOUT ;
  765.    'FINSI' ;
  766. 'FINSI' ;
  767. 'FINSI' ;
  768.  
  769. *--------------------- Fin de modif 1 du 20/02/2013 --------------------
  770.  
  771. *     Elimination des zeros dans TF_UO
  772.  
  773. LABS   = 'PROG' 0.0  FC    (FC + DFREQ) 9000.0       ;
  774. LFIL1  = 'PROG' 1.0  1.0    0.0         0.0          ;
  775. LFIL2  = 'PROG' 0.0  0.0    0.0         0.0          ;
  776. LFIL3  = 'PROG' 0.0  0.0    1.0         1.0          ;
  777.  
  778. LFIL1  = 'IPOL' FR_GAMO LABS LFIL1                   ;
  779. LFIL2  = 'IPOL' FR_GAMO LABS LFIL2                   ;
  780. LFIL3  = 'IPOL' FR_GAMO LABS LFIL3                   ;
  781.  
  782. EVFIL  = 'EVOL' 'COMP' 'MOPH' FR_GAMO LFIL1 LFIL2    ;
  783. EVADD  = 'EVOL' 'COMP' 'MOPH' FR_GAMO LFIL3 LFIL2    ;
  784.  
  785. TF_UO = TF_UO * EVFIL                                ;
  786. TF_UO = TF_UO + EVADD                                ;
  787.  
  788.  
  789. TF_GAMO  = TF_GAMO  * EVFIL                          ;
  790. TF_DEPLO = TF_DEPLO * EVFIL                          ;
  791.  
  792. 'TITR'  'ACCELERATION FILTREE AU POINT DE CONTROLE'  ;
  793.  GAMMAO   = 'TFRI' TF_GAMO                           ;
  794. 'TITR'  'DEPLACEMENT FILTREE AU POINT DE CONTROLE'   ;
  795.  DEPLO    = 'TFRI' TF_DEPLO                          ;
  796.  
  797. 'TITRE' 'DEPLACEMENT A LA BASE'                             ;
  798.  UB   = 'EVOL' 'MANU' 'TEMPS'  L_TEM 'DEPL' U_FRON.NUM_BAS  ;
  799. 'TITRE' 'VITESSE A LA BASE'                                 ;
  800.  VO   = 'EVOL' 'MANU' 'TEMPS'  L_TEM 'VITE' V_FRON.NUM_BAS  ;
  801. 'TITRE' 'DEFORMATION A LA BASE'                             ;
  802.  EO   = 'EVOL' 'MANU' 'TEMPS'  L_TEM 'EPSI' E_FRON.NUM_BAS  ;
  803. 'TITRE' 'REACTION A LA BASE'                                ;
  804.  
  805. 'SI' (('EGA' DIR 'UX') OU ('EGA' DIR 'UY'))                 ;
  806.  RO   = 'EVOL' 'MANU' 'TEMPS'  L_TEM 'REAC' R_FRON.NUM_BAS  ;
  807. 'SINON'                                                     ;
  808. 'SI' ('EGA' DIR 'UZ')                                       ;
  809.  RO1  = 'EVOL' 'MANU' 'TEMPS'  L_TEM 'REAC' R1_FRON.NUM_BAS ;
  810.  RO2  = 'EVOL' 'MANU' 'TEMPS'  L_TEM 'REAC' R2_FRON.NUM_BAS ;
  811. 'FINSI'                                                     ;
  812. 'FINSI'                                                     ;
  813.  
  814.  J = 0                                                ;
  815. 'REPE' BOUCL5 NB_NOEU                                 ;
  816.  J = J + 1                                            ;
  817. 'TITRE'  'FONCTION DE TRANSFERT  HU    DU NOEUD ' J   ;
  818.  U_J   = 'EVOL' 'MANU' 'TEMPS'  L_TEM 'DEPL' U_FRON.J ;
  819.  TF_UJ = 'TFR' NN  U_J 'MOPH'                         ;
  820.  H_DEPL.J = TF_UJ / TF_UO                             ;
  821.  
  822. 'TITRE'  'FONCTION DE TRANSFERT  HV    DU NOEUD' J    ;
  823.  V_J   = 'EVOL' 'MANU' 'TEMPS'  L_TEM 'VITE' V_FRON.J ;
  824.  TF_VJ = 'TFR' NN  V_J 'MOPH'                         ;
  825.  H_VITE.J = TF_VJ / TF_UO                             ;
  826.  
  827. 'TITRE' 'FONCTION DE TRANSFERT   HE    DU NOEUD' J    ;
  828.  E_J   = 'EVOL' 'MANU' 'TEMPS'  L_TEM 'EPSI' E_FRON.J ;
  829.  TF_EJ = 'TFR' NN  E_J 'MOPH'                         ;
  830.  H_EPSI.J = TF_EJ / TF_UO                             ;
  831.  
  832. 'TITRE' 'FONCTION DE TRANSFERT   HR    DU NOEUD' J    ;
  833.  
  834. 'SI' (('EGA' DIR 'UX') OU ('EGA' DIR 'UY'))           ;
  835.  R_J   = 'EVOL' 'MANU' 'TEMPS'  L_TEM 'REAC' R_FRON.J ;
  836.  TF_RJ = 'TFR' NN  R_J 'MOPH'                         ;
  837.  H_REAC.J = TF_RJ / TF_UO                             ;
  838. 'SINON'                                               ;
  839.  
  840. 'SI' ('EGA' DIR 'UZ')                                   ;
  841.  R1_J   = 'EVOL' 'MANU' 'TEMPS'  L_TEM 'REAC' R1_FRON.J ;
  842.  R2_J   = 'EVOL' 'MANU' 'TEMPS'  L_TEM 'REAC' R2_FRON.J ;
  843.  TF_R1J = 'TFR' NN  R1_J 'MOPH'                         ;
  844.  TF_R2J = 'TFR' NN  R2_J 'MOPH'                         ;
  845.  H1_REAC.J = TF_R1J / TF_UO                             ;
  846.  H2_REAC.J = TF_R2J / TF_UO                             ;
  847.  
  848. 'FINSI'                                                 ;
  849. 'FINSI'                                                 ;
  850.  
  851. 'FIN' BOUCL5  ;
  852. 'MENA'        ;
  854.  
  855. *-------------------- Debut de modif 2 du 20/02/2013 -------------------
  856.  
  857. *     Acceleration apres deconvolution
  858.  
  859. 'SAUT' 2 'LIGNE'                       ;
  860. 'MESS' 'ACCELERATION DU CHAMP LIBRE'   ;
  861. 'SI' ( 'EXIS' P_GAMMA )                ;
  862.    NB_GAM = 'DIME' P_GAMMA             ;
  863.    GAMMA  = 'TABLE'                    ;
  864.  
  865.   'SI' ( 'EXIS' P_GAMMA ENTREE )       ;
  866.    NB_GAM = NB_GAM - 1                 ;
  867.   'FINSI'                              ;
  868.  
  869.    K = 0                               ;
  870.   'REPE' BOUCL6 NB_GAM                 ;
  871.    K = K + 1                           ;
  872.    NUM_P  = 'NOEU' P_GAMMA.K           ;
  873.    HU     =  H_DEPL.NUM_P              ;
  874.  
  875.   'TITRE' 'ACCELERATION DU POINT' K                      ;
  876.    TF_GAM   = HU * TF_GAMO                               ;
  877.    GAMMA.K  = 'TFRI' TF_GAM                              ;
  878.    GAMMA.K  = CHTITR ( GAMMA.K ) 'TEMPS' 'ACCE(M/S*S)'   ;
  879.   'FIN' BOUCL6                                           ;
  880.  
  881.    HU     = TF_UOUT / TF_UO                                    ;
  882.   'TITRE' 'ACCELERATION OUTCROP DE BEDROCK'                    ;
  883.    TF_GAM   = HU * TF_GAMO                                     ;
  884.    GAMMA.OUTCROP  = 'TFRI' TF_GAM                              ;
  885.    GAMMA.OUTCROP  = CHTITR GAMMA.OUTCROP 'TEMPS' 'ACCE(M/S*S)' ;
  886. 'FINSI' ;
  887.  
  888. *--------------------- Fin de modif 2 du 20/02/2013 --------------------
  889.  
  890. 'SAUT' 'LIGNE'                             ;
  891. 'MESS' 'VITESSE, DEFORMATION ET REACTION'  ;
  892.  
  893.  J = 0                                     ;
  894. 'REPE' BOUCL7 NB_NOEU                      ;
  895.  J = J + 1                                 ;
  896.  
  897. *    Vitesse compatible avec GAMMAO
  898.  
  899. 'TITRE' 'VITESSE DU NOEUD' J               ;
  900.  TF_VITE  = ( H_VITE.J ) * TF_DEPLO        ;
  901.  V_FRON.J = 'TFRI' TF_VITE                 ;
  902.  
  903. *    On met la vitesse initiale (moyenne sur NV point) a zero
  904.  
  905.  NV = 1                                              ;
  906.  LT     = 'EXTR' ( V_FRON.J ) 'ABSC'                 ;
  907.  LY     = 'EXTR' ( V_FRON.J ) 'ORDO'                 ;
  908.  
  909.  Y1 = 0                                              ;
  910.  IV = 1                                              ;
  911. 'REPE' BV NV                                         ;
  912.  Y1 = Y1 + ('EXTR' LY IV)                            ;
  913.  IV = IV + 1                                         ;
  914. 'FIN' BV                                             ;
  915.  Y1 = Y1 / NV ;
  916.  
  917.  NY     = 'DIME' LY                                  ;
  918.  LY1    = 'PROG' NY * Y1                             ;
  919.  LY     =  LY - LY1                                  ;
  920.  V_FRON.J = 'EVOL' 'MANU' 'TEMPS' LT 'VITE(M/S)' LY  ;
  921.  
  922. *    Deformation
  923.  
  924. 'TITRE' 'DEFORMATION AU NOEUD' J                    ;
  925.  TF_EPSI  = H_EPSI.J * TF_DEPLO                     ;
  926.  E_FRON.J = 'TFRI' TF_EPSI                          ;
  927.  
  928. *    Reactions des appuis : DIR = UX ou UY, Reaction FZ
  929. *                           DIR = UZ, Reaction Fx et FY
  930.  
  931.  
  932. 'SI' (('EGA' DIR 'UX') OU ('EGA' DIR 'UY'))         ;
  933. 'TITRE' 'REACTION VERTICALE AU NOEUD' J             ;
  934.  TF_REAC  = ( H_REAC.J ) * TF_DEPLO                 ;
  935.  R_FRON.J = 'TFRI' TF_REAC                          ;
  936.  
  937. *    On met la reaction intiale a zero
  938.  
  939.  LT     = 'EXTR' ( R_FRON.J ) 'ABSC'                ;
  940.  LY     = 'EXTR' ( R_FRON.J ) 'ORDO'                ;
  941.  Y1     = 'EXTR' LY 1                               ;
  942.  NY     = 'DIME' LY                                 ;
  943.  LY1    = 'PROG' NY * Y1                            ;
  944.  LY     =  LY - LY1                                 ;
  945.  R_FRON.J = 'EVOL' 'MANU' 'TEMPS' LT 'REAC' LY      ;
  946. 'SINON'                                             ;
  947.  
  948. 'SI' ('EGA' DIR 'UZ')    ;
  949. 'TITRE' 'REACTIONS HORIZONTALES AU NOEUD' J         ;
  950.  TF1_REAC  = ( H1_REAC.J ) * TF_DEPLO               ;
  951.  TF2_REAC  = ( H2_REAC.J ) * TF_DEPLO               ;
  952.  R1_FRON.J = 'TFRI' TF1_REAC                        ;
  953.  R2_FRON.J = 'TFRI' TF2_REAC                        ;
  954.  
  955. *    On met les reactions intiales a zero
  956.  
  957.  LT     = 'EXTR' ( R1_FRON.J ) 'ABSC'               ;
  958.  LY     = 'EXTR' ( R1_FRON.J ) 'ORDO'               ;
  959.  Y1     = 'EXTR' LY 1                               ;
  960.  NY     = 'DIME' LY                                 ;
  961.  LY1    = 'PROG' NY * Y1                            ;
  962.  LY     =  LY - LY1                                 ;
  963.  R1_FRON.J = 'EVOL' 'MANU' 'TEMPS' LT 'REAC' LY     ;
  964.  
  965.  LT     = 'EXTR' ( R2_FRON.J ) 'ABSC'               ;
  966.  LY     = 'EXTR' ( R2_FRON.J ) 'ORDO'               ;
  967.  Y1     = 'EXTR' LY 1                               ;
  968.  NY     = 'DIME' LY                                 ;
  969.  LY1    = 'PROG' NY * Y1                            ;
  970.  LY     =  LY - LY1                                 ;
  971.  R2_FRON.J = 'EVOL' 'MANU' 'TEMPS' LT 'REAC' LY     ;
  972.  
  973. 'FINSI'                                             ;
  974. 'FINSI'                                             ;
  975.  
  976. 'FIN'  BOUCL7                                       ;
  977. ************************************************************************
  978.  
  979. *    Force d'amortissement a la base
  980.  
  981. 'SAUT' 'LIGNE'                                       ;
  982. 'MESS' 'FORCE D"AMORTISSEMENT A LA BASE : LYSMER'    ;
  983. 'TITR' 'CONTRAINTE D"AMORTISSEMENT A LA BASE'        ;
  984.  G  = ( COUCHE.NB_COU.'MODULE_G' )                   ;
  985.  RO = ( COUCHE.NB_COU.'MASSE_VOLUMIQUE' )            ;
  986.  NU = ( COUCHE.NB_COU.'POISSON')                     ;
  987.  CS = ( G / RO ) ** 0.5                              ;
  988.  CP = ((2.0 * (1.0 - NU))/(1.0 - (2.0 * NU))) ** 0.5 ;
  989.  CP = CP * CS                                        ;
  990.  
  991. 'SI' (('EGA' DIR 'UX') OU ('EGA' DIR 'UY'))          ;
  992.  SIG_AMO =  RO * CS * V_FRON.NUM_BAS                 ;
  993.  SIG_AMO =  CHTITR SIG_AMO 'TEMPS(S)' 'SIGMA(PA)'    ;
  994.  
  995. 'SINON'                                              ;
  996. 'SI'('EGA' DIR 'UZ')                                 ;
  997.  SIG_AMO =  RO * CP * V_FRON.NUM_BAS                 ;
  998.  SIG_AMO =  CHTITR SIG_AMO 'TEMPS(S)' 'SIGMA(PA)'    ;
  999. 'FINSI'                                              ;
  1000. 'FINSI'                                              ;
  1001.  
  1002. *    Contrainte a la base
  1003.  
  1004. 'MESS' 'CONTRAINTE A LA BASE'                              ;
  1005. 'MESS'                                                     ;
  1006. 'TITR' 'IMPULSION A LA BASE'                               ;
  1007.  EV_FOR = 'EVOL' 'MANU' L_TEM ( 'IPOL' L_TEM L_TIM L_FOR ) ;
  1008.  
  1009. 'TITR' 'TFR DE L"IMPULSION'                                ;
  1010.  TF_FOR = 'TFR' NN  EV_FOR 'MOPH'                          ;
  1011. 'TITR' 'FONCTION DE TRANSFERT DE L"IMPULSION'              ;
  1012.  HS     = TF_FOR / TF_UO                                   ;
  1013.  TF_SIG =  HS * TF_DEPLO                                   ;
  1014.  
  1015. 'TITR' 'CONTRAINTE A LA BASE'                              ;
  1016.  SIG_H  = ( 'TFRI' TF_SIG  ) - SIG_AMO                     ;
  1017.  
  1018. *   On met la contrainte initiale a zero
  1019.  
  1020. LT     = 'EXTR' SIG_H 'ABSC'                        ;
  1021. LY     = 'EXTR' SIG_H 'ORDO'                        ;
  1022. Y1     = 'EXTR' LY 1                                ;
  1023. NY     = 'DIME' LY                                  ;
  1024. LY1    = 'PROG' NY * Y1                             ;
  1025. LY     =  LY - LY1                                  ;
  1026.  
  1027. SIG_H  = 'EVOL' 'MANU' 'TEMPS(S)' LT 'SIGMA(PA)' LY ;
  1028.  
  1029. 'MENA'                                              ;
  1031. ************************************************************************
  1032. ***              FRONTIERE ABSORBANTE DU SOL PROCHE                  ***
  1033. ************************************************************************
  1034. 'SAUT' 3 'LIGNE'                            ;
  1035. 'MESS' 'FRONTIERE ABSORBANTE DU SOL PROCHE' ;
  1036. 'SAUT' 'LIGNE'                              ;
  1037. 'MESS' 'FRONTIERE VERTICALE DU SOL '        ;
  1038.  
  1039.  I = 0                                                             ;
  1040. 'REPE' BOUCL8 NB_COU                                               ;
  1041.  I = I + 1                                                         ;
  1042. 'SI' ( I 'EGA' 1 )                                                 ;
  1043.  AM_V = 'AMOR' COLMOD.I (COUCHE.I.'FRONTIERE') COLMAT.I            ;
  1044.  
  1045. 'SINON'                                                            ;
  1046.  AM_V = AM_V 'ET'('AMOR' COLMOD.I (COUCHE.I.'FRONTIERE') COLMAT.I) ;
  1047. 'FINSI'                                                            ;
  1048. 'FIN' BOUCL8                                                       ;
  1049.  
  1050. 'MESS' 'FRONTIERE HORIZONTALE DU SOL'                 ;
  1051.  AM_H = 'AMOR' MOD_SOL FOND_SOL MAT_SOL               ;
  1052.  AMOF = AM_V 'ET' AM_H                                ;
  1053.  
  1054. ************************************************************************
  1055. ***         CHARGEMENT SISMIQUE SUR LA FRONTIERE HORIZONTALE         ***
  1056. ************************************************************************
  1057.  
  1058. 'SAUT' 3 'LIGNE'                                          ;
  1059. 'MESS' 'CHARGEMENT SISMIQUE SUR LA FRONTIERE HORIZONTALE' ;
  1060.  
  1061. *    FO : forces nodales equivalentes de la contrainte
  1062.  
  1063. 'SI' ( 'EGA' DIR 'UX' )                                   ;
  1064.  FOR_P = 'PRES' 'MASS' MOD_SOL 1.0 FOND_SOL               ;
  1065.  FOR_H = 'EXCO' 'FZ' FOR_P 'FX'                           ;
  1066. 'SINON'                                                   ;
  1067.  
  1068. 'SI' ( 'EGA' DIR 'UY' )                                   ;
  1069.  FOR_P = 'PRES' 'MASS' MOD_SOL 1.0 FOND_SOL               ;
  1070.  FOR_H = 'EXCO' 'FZ' FOR_P 'FY'                           ;
  1071. 'SINON'                                                   ;
  1072.  
  1073. 'SI' ( 'EGA' DIR 'UZ' )                                   ;
  1074.  FOR_P = 'PRES' 'MASS' MOD_SOL 1.0 FOND_SOL               ;
  1075.  FOR_H =  FOR_P                                           ;
  1076.  
  1077. 'FINSI'                                                   ;
  1078. 'FINSI'                                                   ;
  1079. 'FINSI'                                                   ;
  1080.  
  1081.  CHA_HFO = 'CHAR' FOR_H SIG_H                             ;
  1082.  
  1083. *    Amortissement * vitesse (A*VO)
  1084.  
  1085. 'SI' (NON BIELAK) ;
  1086.  
  1087. 'SI' ( 'EGA' DIR 'UX' )                    ;
  1088.  CH_V = 'MANU' 'CHPO' FOND_SOL 1 'UX' 1.0  ;
  1089. 'SINON'                                    ;
  1090.  
  1091. 'SI' ( 'EGA' DIR 'UY' )                    ;
  1092.  CH_V = 'MANU' 'CHPO' FOND_SOL 1 'UY' 1.0  ;
  1093. 'SINON'                                    ;
  1094.  
  1095. 'SI' ( 'EGA' DIR 'UZ' )                    ;
  1096.  CH_V = 'MANU' 'CHPO' FOND_SOL 1 'UZ' 1.0  ;
  1097. 'FINSI'                                    ;
  1098. 'FINSI'                                    ;
  1099. 'FINSI'                                    ;
  1100.  
  1101.  F_BAS   =  AM_H * CH_V                    ;
  1102.  CHA_HAX = 'CHAR' F_BAS ( V_FRON.NUM_BAS ) ;
  1103. 'MENA'                                     ;
  1105.  
  1106. 'FINSI' ;
  1107.  
  1108. ************************************************************************
  1109. ***         CHARGEMENT SISMIQUE SUR LA FRONTIERE VERTICALE           ***
  1110. ************************************************************************
  1111.  
  1112. *     FO = - Reaction
  1113.  
  1114. 'MESS' 'CHARGEMENT SISMIQUE SUR LA FRONTIERE VERTICALE'             ;
  1115.  
  1116. 'SI' (('EGA' DIR 'UX') OU ('EGA' DIR 'UY'))                         ;
  1117.  J = 0                                                              ;
  1118. 'REPE' BOUCL9_1 NB_NOEU                                             ;
  1119.  J = J + 1                                                          ;
  1120. 'SI' ( J 'EGA'  1 )                                                 ;
  1121.  CHA_VFO =  'CHAR' ( 'FORC' FMOT -1.0 POI_FRO.J ) ( R_FRON.J )      ;
  1122. 'SINON'                                                             ;
  1123.  CHA_VFO = ('CHAR' ( 'FORC' FMOT -1.0 POI_FRO.J ) ( R_FRON.J ))
  1124.             'ET' CHA_VFO                                            ;
  1125. 'FINSI'                                                             ;
  1126. 'FIN' BOUCL9_1                                                      ;
  1127.  
  1128. 'SINON'                                                             ;
  1129. 'SI' ('EGA' DIR 'UZ')                                               ;
  1130.  J = 0                                                              ;
  1131. 'REPE' BOUCL9_2 NB_NOEU                                             ;
  1132.  J = J + 1                                                          ;
  1133. 'SI' ( J 'EGA'  1 )                                                 ;
  1134.  CHA_VFO = ('CHAR' ( 'FORC' FMOT1 -1.0 POI_FRO.J ) ( R1_FRON.J )) ET
  1135.            ('CHAR' ( 'FORC' FMOT2 -1.0 POI_FRO.J ) ( R2_FRON.J ))   ;
  1136. 'SINON'                                                             ;
  1137.  CHA_VFO = ('CHAR' ( 'FORC' FMOT1 -1.0 POI_FRO.J ) ( R1_FRON.J )) ET
  1138.            ('CHAR' ( 'FORC' FMOT2 -1.0 POI_FRO.J ) ( R2_FRON.J )) ET
  1139.              CHA_VFO                                                ;
  1140. 'FINSI'                                                             ;
  1141. 'FIN' BOUCL9_2                                                      ;
  1142.  
  1143. 'FINSI'                                                             ;
  1144. 'FINSI'                                                             ;
  1145.  
  1146. *     Amortissement  * viteese : A*VO
  1147.  
  1148. 'SI' (NON BIELAK) ;
  1149.  
  1150.  V_ABSC = 'EXTR' (V_FRON. 1) 'ABSC'  ;
  1151.  N_PAS  = 'DIME' V_ABSC              ;
  1152.  
  1153. *    DIR = UZ
  1154.  
  1155. 'SI' ( 'EGA' DIR 'UZ' )              ;
  1156.  
  1157.  VISC = 'TABLE'                      ;
  1158.  K = 0                               ;
  1159. 'REPE' BOUCL10 NB_NOEU               ;
  1160.  K = K + 1                           ;
  1161.  EV01 = V_FRON.K                     ;
  1162.  V_FRON.K = 'EXTR' EV01 'ORDO'       ;
  1163.  VISC.K   = 'PROG' N_PAS * 0.0       ;
  1164. 'FIN' BOUCL10                        ;
  1165.  
  1166.  J = 0                                           ;
  1167. 'REPE' BOUCL11 NB_NOEU                           ;
  1168.  J = J + 1                                       ;
  1169.  VEC_U = 'MANU' 'CHPO' POI_FRO.J 1 'UZ' 1.0      ;
  1170.  VEC_F = AM_V * VEC_U                            ;
  1171.  
  1172.  K = 0                                           ;
  1173. 'REPE' BOUC11B NB_NOEU                           ;
  1174.  K = K + 1                                       ;
  1175.  AM_JK = 'EXTR' VEC_F POI_FRO.K 'FZ'             ;
  1176.  VISC.K = VISC.K + ( AM_JK * V_FRON.J )          ;
  1177. 'FIN' BOUC11B                                    ;
  1178. 'FIN' BOUCL11                                    ;
  1179.  
  1180. *     Recombinaison
  1181.  
  1182.  FOR_H = 'FORC' 'FZ' 1.0  POI_FRO.1              ;
  1183.  CHA_VAX =  'CHAR' FOR_H ( 'EVOL' 'MANU'
  1184.              V_ABSC ( VISC. 1 ) )                ;
  1185.  K = 1                                           ;
  1186. 'REPE' BOUCL12 ( NB_NOEU - 1 )                   ;
  1187.  K = K + 1                                       ;
  1188.  FOR_H =  'FORC' 'FZ' 1.0 POI_FRO.K              ;
  1189.  CHA_VAX =( 'CHAR' FOR_H ( 'EVOL' 'MANU'
  1190.              V_ABSC ( VISC.K ) ) ) 'ET' CHA_VAX  ;
  1191. 'FIN' BOUCL12                                    ;
  1192. 'FINSI'                                          ;
  1193.  
  1194. *     Dir = UX ou UY
  1195.  
  1196. 'SI' (( 'EGA' DIR 'UX' ) OU ( 'EGA' DIR 'UY' ))  ;
  1197.  UDIR = DIR                                      ;
  1198.  
  1199.  VISC_X = 'TABLE'                                ;
  1200.  VISC_Y = 'TABLE'                                ;
  1201.  K = 0                                           ;
  1202. 'REPE' BOUCL16 NB_NOEU                           ;
  1203.  K = K + 1                                       ;
  1204.  EV01 = V_FRON.K                                 ;
  1205.  V_FRON.K = 'EXTR' EV01 'ORDO'                   ;
  1206.  VISC_X.K = 'PROG' N_PAS * 0.0                   ;
  1207.  VISC_Y.K = 'PROG' N_PAS * 0.0                   ;
  1208. 'FIN' BOUCL16                                    ;
  1209.  J = 0                                           ;
  1210. 'REPE' BOUCL17 NB_NOEU                           ;
  1211.  J = J + 1                                       ;
  1212.  VEC_U = 'MANU' 'CHPO' POI_FRO.J 1 UDIR 1.0      ;
  1213.  VEC_F = AM_V * VEC_U                            ;
  1214.  K = 0                                           ;
  1215. 'REPE' BOUC17B NB_NOEU                           ;
  1216.  K = K + 1                                       ;
  1217.  AMX_JK = 'EXTR' VEC_F POI_FRO.K 'FX'            ;
  1218.  AMY_JK = 'EXTR' VEC_F POI_FRO.K 'FY'            ;
  1219.  VISC_X.K = VISC_X.K + ( AMX_JK * V_FRON.J )     ;
  1220.  VISC_Y.K = VISC_Y.K + ( AMY_JK * V_FRON.J )     ;
  1221. 'FIN' BOUC17B                                    ;
  1222. 'FIN' BOUCL17                                    ;
  1223.  
  1224. *    Recombinaison
  1225.  
  1226.  FOR_X = 'FORC' 'FX' 1.0 POI_FRO.1               ;
  1227.  FOR_Y = 'FORC' 'FY' 1.0 POI_FRO.1               ;
  1228.  CHA_VAX = 'CHAR' FOR_X ( 'EVOL' 'MANU'
  1229.             V_ABSC ( VISC_X. 1 ) )               ;
  1230.  CHA_VAY = 'CHAR' FOR_Y ( 'EVOL' 'MANU'
  1231.             V_ABSC ( VISC_Y. 1 ) )               ;
  1232.  K = 1                                           ;
  1233. 'REPE' BOUCL18 ( NB_NOEU - 1 )                   ;
  1234.  K = K + 1                                       ;
  1235.  FOR_X =  'FORC' 'FX' 1.0 POI_FRO.K              ;
  1236.  CHA_VAX = CHA_VAX 'ET' ( 'CHAR' FOR_X ( 'EVOL'
  1237.           'MANU' V_ABSC ( VISC_X.K ) ) )         ;
  1238.  FOR_Y =  'FORC' 'FY' 1.0 POI_FRO.K              ;
  1239.  CHA_VAY = CHA_VAY 'ET' ( 'CHAR' FOR_Y ( 'EVOL'
  1240.           'MANU' V_ABSC ( VISC_Y.K ) ) )         ;
  1241. 'FIN' BOUCL18                                    ;
  1242.  CHA_VAX = CHA_VAX 'ET' CHA_VAY                  ;
  1243.  
  1244. 'FINSI'                                          ;
  1245.  
  1246. 'FINSI' ;
  1247.  
  1248. 'SI' BIELAK ;
  1249.  CHA1 = CHA_HFO 'ET' CHA_VFO ;
  1250. 'SINON' ;
  1251.  CHA1 = CHA_HFO 'ET' CHA_HAX 'ET' CHA_VFO 'ET' CHA_VAX ;
  1252. 'FINSI' ;
  1253.  
  1254. * Changement de format du chargement
  1255.  
  1256.  CHA1 = 'CHAN' 'TABL' CHA1 'MECA' ;
  1257.  
  1258. ************************************************************************
  1259. ***               Deformation de cisaillement maximale               ***
  1260. ************************************************************************
  1261. 'SAUT' 2 'LIGNE'                               ;
  1262. 'MESS' 'DEFORMATION MAXIMALE'                  ;
  1263.  PROF = 'PROG' NB_NOEU * 0.0                   ;
  1264.  GMAX = 'PROG' NB_NOEU * 0.0                   ;
  1265.  J = 0                                         ;
  1266. 'REPE' BOUCL19 NB_NOEU                         ;
  1267.  J = J + 1                                     ;
  1268.  ZJ   = 'COOR' 3 ( POI_FRO.J )                 ;
  1269.  CISA = 'EXTR' ( E_FRON.J ) 'ORDO'             ;
  1270.  GX   = 'MAXI' CISA 'ABS'                      ;
  1271. 'REMP' PROF J ZJ                               ;
  1272. 'REMP' GMAX J GX                               ;
  1273. 'FIN' BOUCL19                                  ;
  1274.  GMAXMAX = 'MAXI' GMAX                         ;
  1275.  
  1276. *     Evolution de deformation en profondeur
  1277.  
  1278. 'SI' (('EGA' DIR 'UX') OU ('EGA' DIR 'UY'))                 ;
  1279. 'TITR' 'DEFORMATION DE CISAILLEMENT MAXIMALE'               ;
  1280.  DEFZ  = 'EVOL' 'MANU' 'PROFONDEUR' PROF 'GAM_MAX' GMAX     ;
  1281.  DEFZ  = 'ORDO' 'DECROISSANT' DEFZ                          ;
  1282.  ZFROF = 'EXTR' DEFZ 'ABSC'                                 ;
  1283.  ZGMAX = 'EXTR' DEFZ 'ORDO'                                 ;
  1284.  DEFZ  = 'EVOL' 'MANU' 'GAM_MAX' ZGMAX 'PROFONDEUR' ZFROF   ;
  1285. 'DESS' DEFZ 'NCLK'                                   ;
  1286. 'FINSI'                                                     ;
  1287.  
  1288. 'SI' ( 'EGA' DIR 'UZ' )                                     ;
  1289. 'TITR' 'DEFORMATION VERTICALE MAXIMALE'                     ;
  1290.  DEFZ  = 'EVOL' 'MANU' 'PROFONDEUR' PROF 'EPS_MAX' GMAX     ;
  1291.  DEFZ  = 'ORDO' 'DECROISSANT' DEFZ                          ;
  1292.  ZFROF = 'EXTR' DEFZ 'ABSC'                                 ;
  1293.  ZGMAX = 'EXTR' DEFZ 'ORDO'                                 ;
  1294.  DEFZ  = 'EVOL' 'MANU' 'EPS_MAX' ZGMAX 'PROFONDEUR' ZFROF   ;
  1295. 'DESS' DEFZ 'NCLK'                                   ;
  1296. 'FINSI'                                                     ;
  1297.  
  1298. *-------------DEBUT DE MODIFICATION 01/02/2014 (F. WANG)----------------
  1299. *
  1300. *-----------------------------------------------------------------------
  1301. *  IMPLANTATION DE LA METHODE DE BIELAK :
  1302. *  GENERATION DE LA ZONE INTERMEDIAIRE ET LE CHARGEMENT CORRESPONDANT
  1303. *-----------------------------------------------------------------------
  1304.  
  1305. 'SI' BIELAK ;
  1306.  
  1307. *  Calcul de dépalcement et d'acceleration
  1308.  
  1309.  A_FRON = 'TABLE'                                   ;
  1310.  J = 0                                              ;
  1311. 'REPE' BBLK1 NB_NOEU                                ;
  1312.  J = J + 1                                          ;
  1313. 'TITRE' 'DEPLACEMENT AU NOEUD' J                    ;
  1314.  TF_DEPL  = ( H_DEPL.J ) * TF_DEPLO                 ;
  1315.  U_FRON.J = 'TFRI' TF_DEPL                          ;
  1316.  
  1317.  LT     = 'EXTR' ( U_FRON.J ) 'ABSC'                ;
  1318.  LY     = 'EXTR' ( U_FRON.J ) 'ORDO'                ;
  1319.  Y1     = 'EXTR' LY 1                               ;
  1320.  NY     = 'DIME' LY                                 ;
  1321.  LY1    = 'PROG' NY * Y1                            ;
  1322.  LY     =  LY - LY1                                 ;
  1323.  U_FRON.J = 'EVOL' 'MANU' 'TEMPS' LT 'DEPL' LY      ;
  1324.  
  1325. 'TITR' 'ACCELARATION AU NOEUD' J                    ;
  1326. * A_FRON.J = DERIV1 V_FRON.J                         ;
  1327.  
  1328.  TF_ACCE = ( H_DEPL.J ) * TF_GAMO                   ;
  1329.  A_FRON.J =  'TFRI' TF_ACCE                         ;
  1330.  
  1331.  LT     = 'EXTR' ( A_FRON.J ) 'ABSC'                ;
  1332.  LY     = 'EXTR' ( A_FRON.J ) 'ORDO'                ;
  1333.  Y1     = 'EXTR' LY 1                               ;
  1334.  NY     = 'DIME' LY                                 ;
  1335.  LY1    = 'PROG' NY * Y1                            ;
  1336.  LY     =  LY - LY1                                 ;
  1337.  A_FRON.J = 'EVOL' 'MANU' 'TEMPS' LT 'ACCE' LY      ;
  1338. 'FIN' BBLK1                                         ;
  1339.  
  1340. * RECUPERATION DE LA COLONNE DE DECONVOLUTION
  1341.  
  1342.  BLKV_MAI = COL1 ;
  1343.  BLKV_MOD = MOD1 ;
  1344.  BLKV_MAT = MAT1 ;
  1345.  BLKV_MAS = MAS1 ;
  1346.  BLKV_RIG = RIG1 ;
  1347.  BLKV_AMO = AMO1 ;
  1348.  
  1349. * GENERATION DE LA ZONE HORIZONTALE
  1350.  
  1351.  VEC_V    = 0 0 (-1*COUCHE.BLK.EP_V)                                  ;
  1352.  BLKH_MAI =  VOLU TRAN 1 (FOND_SOL ET BAS_COL) VEC_V                  ;
  1353.  BLKH_MOD = 'MODE' BLKH_MAI 'MECANIQUE' 'ELASTIQUE' 'ISOTROPE'
  1354.              TY1 TY2                                                  ;
  1355.  BLKH_MAT = 'MATE' BLKH_MOD 'YOUNG' (COUCHE.NB_COU.'YOUNG')
  1356.             'NU'  (COUCHE.NB_COU.'POISSON')
  1357.             'RHO' (COUCHE.NB_COU.'MASSE_VOLUMIQUE')                   ;
  1358.  BLKH_MAS = 'MASS' BLKH_MOD BLKH_MAT                                  ;
  1359.  BLKH_RIG = 'RIGI' BLKH_MOD BLKH_MAT                                  ;
  1360.  BLKH_AMO =  KSI * ((ALPHA * BLKH_MAS) 'ET' (BETA * BLKH_RIG))        ;
  1361.  
  1362. * ZONE INTERMEDIAIRE COMPLETE
  1363.  
  1364.  BLK_MAI = BLKV_MAI 'ET' BLKH_MAI ;
  1365.  BLK_MOD = BLKV_MOD 'ET' BLKH_MOD ;
  1366.  BLK_MAT = BLKV_MAT 'ET' BLKH_MAT ;
  1367.  BLK_MAS = BLKV_MAS 'ET' BLKH_MAS ;
  1368.  BLK_RIG = BLKV_RIG 'ET' BLKH_RIG ;
  1369.  BLK_AMO = BLKV_AMO 'ET' BLKH_AMO ;
  1370.  
  1371. * TRACE DE LA ZONE
  1372.  
  1373. 'TITR' 'ZONE INTERNE ET INTERMEDIAIRE'  ;
  1374. 'ELIM' BLK_MAI  0.001                   ;
  1375.  BLK_MAI =  BLK_MAI 'COUL' 'ROUG'       ;
  1376. 'TRAC' (SOL 'ET' BLK_MAI) 'FACE' 'NCLK' ;
  1377.  
  1378. *-----------------------------------------------------------------------
  1379. *    CALCUL DU CHARGEMENT SUR LA ZONE INTERMEDIAIRE
  1380. *-----------------------------------------------------------------------
  1381.    'SI' ('EGA' DIR 'UX') ;
  1382.     FDIR = 'MOT' 'FX'    ;
  1383.    'FINSI'               ;
  1384.    'SI' ('EGA' DIR 'UY') ;
  1385.     FDIR = 'MOT' 'FY'    ;
  1386.    'FINSI'               ;
  1387.    'SI' ('EGA' DIR 'UZ') ;
  1388.     FDIR = 'MOT' 'FZ'    ;
  1389.    'FINSI'               ;
  1390.  
  1391. *  FRONTIERE HORIZONTALE
  1392.  
  1393.  NBFOND = 'NBNO' FOND_SOL                    ;
  1394.  F1     =  1.* NBFOND                        ;
  1395.  FOR_H  = 'FORC' FDIR F1 FOND_SOL            ;
  1396.  CHU_H  = 'CHAR' FOR_H U_FRON.NUM_BAS        ;
  1397.  CHV_H  = 'CHAR' FOR_H V_FRON.NUM_BAS        ;
  1398.  CHA_H  = 'CHAR' FOR_H A_FRON.NUM_BAS        ;
  1399.  
  1400. *  FRONTIERE VERTICALE (SAUF LES NOEUDS A LA BASE
  1401. *  DEJA COMPTES SUR LA FRONTIERE HORIZONTALE)
  1402.  
  1403.     FOR_V  = 'FORC' FDIR 0.0 POI_FRO.1              ;
  1404.     CHU_V  = 'CHAR' FOR_V U_FRON.1                  ;
  1405.     CHV_V  = 'CHAR' FOR_V V_FRON.1                  ;
  1406.     CHA_V  = 'CHAR' FOR_V A_FRON.1                  ;
  1407.  
  1408.     K = 0                                           ;
  1409.    'REPE' BBLK2 NB_NOEU                             ;
  1410.     K = K + 1                                       ;
  1411.       ZK = COOR 3 POI_FRO.K                         ;
  1412.      'SI' ((ABS (ZK - ZBAS)) > 0.001)               ;
  1413.       FOR_V  = 'FORC' FDIR 1.0  POI_FRO.K           ;
  1414.       CHU_V  = CHU_V 'ET' ('CHAR' FOR_V U_FRON.K )  ;
  1415.       CHV_V  = CHV_V 'ET' ('CHAR' FOR_V V_FRON.K )  ;
  1416.       CHA_V  = CHA_V 'ET' ('CHAR' FOR_V A_FRON.K )  ;
  1417.      'FINSI'                                        ;
  1418.    'FIN' BBLK2                                      ;
  1419.  
  1420. *  FRONTIERE HORIZONTALE + VERTICALE
  1421.  
  1422.  CHU = CHU_V 'ET' CHU_H ;
  1423.  CHV = CHV_V 'ET' CHV_H ;
  1424.  CHA = CHA_V 'ET' CHA_H ;
  1425.  
  1426.  CHU = 'CHAN' 'TABL' CHU 'MECA' ;
  1427.  CHV = 'CHAN' 'TABL' CHV 'MECA' ;
  1428.  CHA = 'CHAN' 'TABL' CHA 'MECA' ;
  1429.  
  1430.  TAB_T TAB_U = 'EXTR' CHU 'MECA' 'TABL' ;
  1431.  TAB_T TAB_V = 'EXTR' CHV 'MECA' 'TABL' ;
  1432.  TAB_T TAB_A = 'EXTR' CHA 'MECA' 'TABL' ;
  1433.  
  1434.  LM1 = 'MOTS' FDIR ;
  1435.  LM2 = 'MOTS' DIR ;
  1436.  
  1437.  TAB_F = 'TABLE' ;
  1438.  TAB_FD = 'TABLE' ;
  1439.  TAB_FV = 'TABLE' ;
  1440.  TAB_FA = 'TABLE' ;
  1441.  
  1442.  N_T   = 'DIME' TAB_T ;
  1443.  
  1444.  I = 0 ;
  1445. 'REPE' BBLK3 N_T ;
  1446.  CHU_I = 'EXCO' LM1 (TAB_U.I) LM2 ;
  1447.  CHV_I = 'EXCO' LM1 (TAB_V.I) LM2 ;
  1448.  CHA_I = 'EXCO' LM1 (TAB_A.I) LM2 ;
  1449.  
  1450.  CHU_I = CHU_I * BLK_RIG ;
  1451.  CHV_I = CHV_I * BLK_AMO ;
  1452.  CHA_I = CHA_I * BLK_MAS ;
  1453.  TAB_F.I = CHU_I + CHV_I + CHA_I ;
  1454.  
  1455.  TAB_FD.I = CHU_I ;
  1456.  TAB_FV.I = CHV_I ;
  1457.  TAB_FA.I = CHA_I ;
  1458.  
  1459.  I = I + 1  ;
  1460. 'FIN' BBLK3 ;
  1461.  
  1462. *-----------------------------------------------------------------------
  1463. *-----------------------------------------------------------------------
  1464.  
  1465.  CHA2 = 'CHAR' 'MECA' TAB_T TAB_F ;
  1466.  
  1467. 'FINSI' ;
  1468.  
  1469. *  Sortie des resultats
  1470.  
  1471.  TABRESU = 'TABLE'                ;
  1472.  
  1473. 'SI' BIELAK                       ;
  1474.  TABRESU.'CHAR'    = CHA1 ET CHA2 ;
  1475.  
  1476. * -- pour tester la procédure -----
  1477.  TABRESU.'CHA1'    = CHA1 ;
  1478.  TABRESU.'CHA2'    = CHA2 ;
  1479. *----------------------------------
  1480.  
  1481.  TABRESU.'BLK_MAI' = BLK_MAI      ;
  1482.  TABRESU.'BLK_MOD' = BLK_MOD      ;
  1483.  TABRESU.'BLK_MAT' = BLK_MAT      ;
  1484.  TABRESU.'BLK_MAS' = BLK_MAS      ;
  1485.  TABRESU.'BLK_RIG' = BLK_RIG      ;
  1486.  TABRESU.'BLK_AMO' = BLK_AMO      ;
  1487. 'SINON'                           ;
  1488.  TABRESU.'CHAR'  = CHA1           ;
  1489.  TABRESU.'AMOR'  = AMOF           ;
  1490. 'FINSI'                           ;
  1491.  
  1492. *-----------FIN DE MODIFICATION 01/02/2014 (F. WANG)-------------------
  1493.  
  1494.  TABRESU.'DEFO'  = DEFZ    ;
  1495.  TABRESU.'ACCE'  = GAMMA   ;
  1496.  TABRESU.'PAS'   = PAS_TEM ;
  1497.  TABRESU.'FCDYN' = FCDYN   ;
  1498.  
  1499. *----------------------------------------------------------------------
  1500. 'SINON' ;
  1501. *----------------------------------------------------------------------
  1502. *                         BIELAK ET ONDE SH 
  1503. *                DEBUT DE MODIFICATION 01/10/2015. WANG)
  1504. *----------------------------------------------------------------------
  1505.  
  1506.  SI SH ;
  1507.  MESS 'BIELAK ET ONDE SH' ;
  1508.  
  1509. * RECUPERATION DE L ENVELOPPE DE DECONVOLUTION
  1510.  
  1511.  BLKV_MAI = COL1 ;
  1512.  BLKV_MOD = MOD1 ;
  1513.  BLKV_MAT = MAT1 ;
  1514.  BLKV_MAS = MAS1 ;
  1515.  BLKV_RIG = RIG1 ;
  1516.  BLKV_AMO = AMO1 ;
  1517.  
  1518. * GENERATION DE LA ZONE HORIZONTALE
  1519.  
  1520.  VEC_V    = 0 0 (-1*COUCHE.BLK.EP_V)                                  ;
  1521.  BLKH_MAI = 'VOLU' 'TRAN' 1 (FOND_SOL 'ET' BAS_COL) VEC_V             ;
  1522.  BLKH_MOD = 'MODE' BLKH_MAI 'MECANIQUE' 'ELASTIQUE' 'ISOTROPE'
  1523.              TY1 TY2                                                  ;
  1524.  BLKH_MAT = 'MATE' BLKH_MOD 'YOUNG' (COUCHE.NB_COU.'YOUNG')
  1525.             'NU'  (COUCHE.NB_COU.'POISSON')
  1526.             'RHO' (COUCHE.NB_COU.'MASSE_VOLUMIQUE')                   ;
  1527.  BLKH_MAS = 'MASS' BLKH_MOD BLKH_MAT                                  ;
  1528.  BLKH_RIG = 'RIGI' BLKH_MOD BLKH_MAT                                  ;
  1529.  BLKH_AMO =  KSI * ((ALPHA * BLKH_MAS) 'ET' (BETA * BLKH_RIG))        ;
  1530.  
  1531. * ZONE INTERMEDIAIRE COMPLETE
  1532.  
  1533.  BLK_MAI = BLKV_MAI 'ET' BLKH_MAI ;
  1534.  BLK_MOD = BLKV_MOD 'ET' BLKH_MOD ;
  1535.  BLK_MAT = BLKV_MAT 'ET' BLKH_MAT ;
  1536.  BLK_MAS = BLKV_MAS 'ET' BLKH_MAS ;
  1537.  BLK_RIG = BLKV_RIG 'ET' BLKH_RIG ;
  1538.  BLK_AMO = BLKV_AMO 'ET' BLKH_AMO ;
  1539.  
  1540. * TRACE DE LA ZONE INTERNE ET INTERMEDIAIRE
  1541.  
  1542. 'TITR' 'ZONE INTERNE ET INTERMEDIAIRE'  ;
  1543. 'ELIM' BLK_MAI  0.001                   ;
  1544.  BLK_MAI =  BLK_MAI 'COUL' 'ROUG'       ;
  1545. 'TRAC' (SOL 'ET' BLK_MAI) 'FACE' 'NCLK' ;
  1546.  
  1547. * SURFACE INTERIEURE
  1548.  
  1549. 'TITR' 'SURFACE INTERIEURE' ;
  1550.  BLK_INT = FRON 'ET' FOND_SOL ;
  1551. 'TRAC' BLK_INT CACH ;
  1552.  
  1553. * ACCES AUX NOEUDS 
  1554.  
  1555. 'TASS' BLK_INT ;
  1556.  NB_INT = 'NBNO' BLK_INT ;
  1557. 'MESS' 'NB. DE NOEUDS, SURFACE INTERIEURE' NB_INT ;
  1558.  POI_INT = 'TABLE'         ;
  1559.  J = 0                     ;
  1560. 'REPE' BINT1 NB_INT        ;
  1561.  J = J + 1                 ;
  1562.  POI_INT.J = 'NOEU' J      ;
  1563. 'FIN'  BINT1               ;
  1564.  
  1565. 'TASS' BLK_MAI ;
  1566.  NB_ZON = 'NBNO' BLK_MAI ;
  1567. 'MESS' 'NB. DE NOEUDS, ZONE INTERMEDIAIRE' NB_ZON ;
  1568.  POI_ZON = 'TABLE'         ;
  1569.  J = 0                     ;
  1570. 'REPE' BZON1 NB_ZON        ;
  1571.  J = J + 1                 ;
  1572.  POI_ZON.J = 'NOEU' J      ;
  1573. 'FIN'  BZON1               ;
  1574.  
  1575. * TFR: ACCELERATION, VITESSE ET DEPLACEMENT DU POINT DE CONTROLE 
  1576. * (FC a prendre en compte ...)
  1577.  
  1578.  NP = 'DIME' ( 'EXTR' GAMMAO 'ABSC' )         ;
  1579.  XX = ( 'LOG' NP ) / ( 'LOG' 2 )              ;
  1580.  NN = 'ENTI' XX                               ;
  1581. 'SI' ( XX > ( 1.0 * NN ) )                    ;
  1582.  NN = NN + 1                                  ;
  1583. 'FINSI'                                       ;
  1584.  
  1585. 'TITRE' 'TFR ACCELERATION'                    ;
  1586.  TF_GAMMO = 'TFR'  NN GAMMAO 'MOPH'           ;
  1587.  FR_GAMMO = 'EXTR' TF_GAMMO  'ABSC'           ;
  1588.  MO_GAMMO = 'EXTR' TF_GAMMO 'ORDO' 1          ;
  1589.  PH_GAMMO = 'EXTR' TF_GAMMO 'ORDO' 2          ;
  1590.  NFREQ    = 'DIME' FR_GAMMO                   ;
  1591.  DFREQ    = 'EXTR' FR_GAMMO 2                 ;
  1592.  FREQN    = 'EXTR' FR_GAMMO NFREQ             ;
  1593.  MODUN    = 'EXTR' MO_GAMMO NFREQ             ;
  1594.  
  1595. 'SI' ( IDIV NEG 0 )                           ;
  1596.  NN = NN + IDIV                               ;
  1597. 'REMP' MO_GAMMO NFREQ ( MODUN / 2.0 )         ;
  1598.  
  1599.  NFREQ1  = (2**(NN - 1)) + 1                                  ;
  1600.  FR_GAMMO = 'PROG' 0.0 'PAS' DFREQ 'NPAS' ( NFREQ1 - 1 )      ;
  1601.  LIST ( DIME FR_GAMMO ) ;
  1602.  MO_GAMMO = MO_GAMMO 'ET' ('PROG' (NFREQ1 - NFREQ ) * 0.0)    ;
  1603.  PH_GAMMO = PH_GAMMO 'ET' ('PROG' (NFREQ1 - NFREQ ) * 0.0)    ;
  1604.  TF_GAMMO = 'EVOL' 'COMP' 'MOPH' 'FREQ(HZ)' FR_GAMMO 'MODULE'
  1605.             MO_GAMMO 'PHASE' PH_GAMMO                         ;
  1606. 'FINSI' ;
  1607.  
  1608.  OMEGA = 2*PI*FR_GAMMO ;
  1609. 'TITRE' 'TFR VITESSE' ;
  1610.  TFIOM  = 'EVOL' 'COMP' 'REIM' FR_GAMMO (FR_GAMMO*0) OMEGA ;
  1611. 'RIMP' TFIOM ;
  1612.  TF_VITEO = TF_GAMMO / TFIOM ;
  1613.  
  1614. 'TITRE' 'TFR DEPLACEMENT' ;
  1615.  TF_DEPLO = TF_VITEO / TFIOM ;
  1616.  
  1617. * ACCELERATION, VITESSE ET DEPLACEMENT
  1618.  
  1619.  TABDON = 'TABLE' ;
  1620.  TABDON.PDC    = P_GAMMA.ENTREE.CONTROLE ;
  1621.  TABDON.DIR    = DIR ;
  1622.  TABDON.ANGLE  = ANGLE ;
  1623.  TABDON.VS     = CS_I ;
  1624.  TABDON.KSI    = KSI ;
  1625.  TABDON.F1     = F1 ;
  1626.  TABDON.F2     = F2 ;
  1627.  TABDON.LF     = FR_GAMMO ;
  1628.  
  1629. * SURFACE INTERIEURE
  1630.  
  1631.  A_INT = 'TABLE' ;
  1632.  V_INT = 'TABLE' ;
  1633.  D_INT = 'TABLE' ;
  1634.  
  1635.  J = 0 ;
  1636. 'REPE' BINT2 NB_INT ;
  1637.  J = J + 1 ;
  1638. 'TITR' 'FONCTION DE TRANSFERT BLK_INT, NOEUD' J ;
  1639.  TABDON.PDF = POI_INT.J ;
  1640.  H_INT = SHFDT TABDON ;
  1641.  
  1642. 'TITR' 'ACCELARATION BLK_INT, NOEUD' J          ;
  1643.  A_INT.J = 'TFRI' (H_INT * TF_GAMMO)            ;
  1644.  LT      = 'EXTR' A_INT.J 'ABSC'                ;
  1645.  LY      = 'EXTR' A_INT.J 'ORDO'                ;
  1646.  Y1      = 'EXTR' LY 1                          ;
  1647.  NY      = 'DIME' LY                            ;
  1648.  LY1     = 'PROG' NY * Y1                       ;
  1649.  LY      =  LY - LY1                            ;
  1650.  A_INT.J = 'EVOL' 'MANU' 'TEMPS' LT 'DEPL' LY   ;
  1651.  
  1652. 'TITR' 'VITESSE BLK_INT, NOEUD' J               ;
  1653.  V_INT.J = 'TFRI' (H_INT * TF_VITEO)            ; 
  1654.  LT      = 'EXTR' V_INT.J 'ABSC'                ;
  1655.  LY      = 'EXTR' V_INT.J 'ORDO'                ;
  1656.  Y1      = 'EXTR' LY 1                          ;
  1657.  NY      = 'DIME' LY                            ;
  1658.  LY1     = 'PROG' NY * Y1                       ;
  1659.  LY      =  LY - LY1                            ;
  1660.  V_INT.J = 'EVOL' 'MANU' 'TEMPS' LT 'DEPL' LY   ;
  1661.  
  1662. 'TITR' 'DEPLACEMENT BLK_INT, NOEUD' J           ;
  1663.  D_INT.J = 'TFRI' (H_INT * TF_DEPLO)            ;
  1664.  LT      = 'EXTR' D_INT.J 'ABSC'                ;
  1665.  LY      = 'EXTR' D_INT.J 'ORDO'                ;
  1666.  Y1      = 'EXTR' LY 1                          ;
  1667.  NY      = 'DIME' LY                            ;
  1668.  LY1     = 'PROG' NY * Y1                       ;
  1669.  LY      =  LY - LY1                            ;
  1670.  D_INT.J = 'EVOL' 'MANU' 'TEMPS' LT 'DEPL' LY   ;
  1671.  
  1672. 'DETR' H_INT ;
  1673. 'FIN' BINT2 ;
  1674.  
  1675. * ZONE INTERMEDIAIRE
  1676.  
  1677.  A_ZON = 'TABLE' ;
  1678.  V_ZON = 'TABLE' ;
  1679.  D_ZON = 'TABLE' ;
  1680.  
  1681.  J = 0 ;
  1682. 'REPE' BZON2 NB_ZON ;
  1683.  J = J + 1 ;
  1684. 'TITR' 'FONCTION DE TRANSFERT BLK_MAI, NOEUD' J ;
  1685.  TABDON.PDF = POI_ZON.J ;
  1686.  H_ZON = SHFDT TABDON ;
  1687.  
  1688. 'TITR' 'ACCELARATION BLK_MAI, NOEUD' J          ;
  1689.  A_ZON.J = 'TFRI' (H_ZON * TF_GAMMO)            ;
  1690.  LT      = 'EXTR' A_ZON.J 'ABSC'                ;
  1691.  LY      = 'EXTR' A_ZON.J 'ORDO'                ;
  1692.  Y1      = 'EXTR' LY 1                          ;
  1693.  NY      = 'DIME' LY                            ;
  1694.  LY1     = 'PROG' NY * Y1                       ;
  1695.  LY      =  LY - LY1                            ;
  1696.  A_ZON.J = 'EVOL' 'MANU' 'TEMPS' LT 'DEPL' LY   ;
  1697.  
  1698. 'TITR' 'VITESSE BLK_MAI, NOEUD' J               ;
  1699.  V_ZON.J = 'TFRI' (H_ZON * TF_VITEO)            ; 
  1700.  LT      = 'EXTR' V_ZON.J 'ABSC'                ;
  1701.  LY      = 'EXTR' V_ZON.J 'ORDO'                ;
  1702.  Y1      = 'EXTR' LY 1                          ;
  1703.  NY      = 'DIME' LY                            ;
  1704.  LY1     = 'PROG' NY * Y1                       ;
  1705.  LY      =  LY - LY1                            ;
  1706.  V_ZON.J = 'EVOL' 'MANU' 'TEMPS' LT 'DEPL' LY   ;
  1707.  
  1708. 'TITR' 'DEPLACEMENT BLK_MAI, NOEUD' J           ;
  1709.  D_ZON.J = 'TFRI' (H_ZON * TF_DEPLO)            ; 
  1710.  LT      = 'EXTR' D_ZON.J 'ABSC'                ;
  1711.  LY      = 'EXTR' D_ZON.J 'ORDO'                ;
  1712.  Y1      = 'EXTR' LY 1                          ;
  1713.  NY      = 'DIME' LY                            ;
  1714.  LY1     = 'PROG' NY * Y1                       ;
  1715.  LY      =  LY - LY1                            ;
  1716.  D_ZON.J = 'EVOL' 'MANU' 'TEMPS' LT 'DEPL' LY   ;
  1717.  
  1718. 'DETR' H_ZON ;
  1719. 'FIN' BZON2 ;
  1720.  
  1721. *-----------------------------------------------------------------------
  1722. *  CHARGEMENT SUR LA ZONE INTERNE (REACTION PAR LA ZONE INTERMEDIAIRE) 
  1723. *     ZONE INTERMEDIARE COMPLETE SOUMISE AU MOUVEMENT DU CHAMP LIBRE 
  1724. *-----------------------------------------------------------------------
  1725.  
  1726. 'SI' ('EGA' DIR 'UX') ;
  1727.  FDIR = 'MOT' 'FX'    ;
  1728. 'FINSI'               ;
  1729. 'SI' ('EGA' DIR 'UY') ;
  1730.  FDIR = 'MOT' 'FY'    ;
  1731. 'FINSI'               ;
  1732.  
  1733.  FOR_Z  = 'FORC' FDIR 1.0 POI_ZON.1             ;
  1734.  CHD_Z  = 'CHAR' FOR_Z D_ZON.1                  ;
  1735.  CHV_Z  = 'CHAR' FOR_Z V_ZON.1                  ;
  1736.  CHA_Z  = 'CHAR' FOR_Z A_ZON.1                  ;
  1737.  
  1738.  J = 1                                          ;
  1739. 'REPE' BZON3 (NB_ZON - 1)                       ;
  1740.  J = J + 1                                      ;
  1741.    FOR_Z  = 'FORC' FDIR 1.0  POI_ZON.J          ;
  1742.    CHD_Z  = CHD_Z 'ET' ('CHAR' FOR_Z D_ZON.J )  ;
  1743.    CHV_Z  = CHV_Z 'ET' ('CHAR' FOR_Z V_ZON.J )  ;
  1744.    CHA_Z  = CHA_Z 'ET' ('CHAR' FOR_Z A_ZON.J )  ;
  1745. 'FIN' BZON3                                     ;
  1746.  
  1747.  CHD_Z = 'CHAN' 'TABL' CHD_Z 'MECA' ;
  1748.  CHV_Z = 'CHAN' 'TABL' CHV_Z 'MECA' ;
  1749.  CHA_Z = 'CHAN' 'TABL' CHA_Z 'MECA' ;
  1750.  
  1751.  TAB_T TABD_Z = 'EXTR' CHD_Z 'MECA' 'TABL' ;
  1752.  TAB_T TABV_Z = 'EXTR' CHV_Z 'MECA' 'TABL' ;
  1753.  TAB_T TABA_Z = 'EXTR' CHA_Z 'MECA' 'TABL' ;
  1754.  
  1755.  LM1 = 'MOTS' FDIR ;
  1756.  LM2 = 'MOTS' DIR ;
  1757.  
  1758.  TABF_Z  = 'TABLE' ;
  1759.  N_T     = 'DIME' TAB_T ;
  1760.  
  1761.  I = 0 ;
  1762. 'REPE' BZON4 N_T ;
  1763.  CHD_Z = 'EXCO' LM1 (TABD_Z.I) LM2 ;
  1764.  CHV_Z = 'EXCO' LM1 (TABV_Z.I) LM2 ;
  1765.  CHA_Z = 'EXCO' LM1 (TABA_Z.I) LM2 ;
  1766.  
  1767.  CHD_Z = CHD_Z * BLK_RIG ;
  1768.  CHV_Z = CHV_Z * BLK_AMO ;
  1769.  CHA_Z = CHA_Z * BLK_MAS ;
  1770.  
  1771.  TABF_Z.I  = CHD_Z + CHV_Z + CHA_Z ;
  1772.  I = I + 1  ;
  1773. 'FIN' BZON4 ;
  1774.  
  1775.  TABREA = TABLE ;
  1776.  I = 0 ;
  1777. 'REPE' BINT5 N_T ;
  1778.  TABREA.I = -1*(REDU BLK_INT TABF_Z.I);
  1779.  I = I + 1  ;
  1780. 'FIN' BINT5 ;
  1781.  
  1782.  CHA1 = 'CHAR' 'MECA' TAB_T TABREA ;
  1783.  
  1784. *-----------------------------------------------------------------------
  1785. *                 CHARGEMENT SUR LA ZONE INTERMEDIAIRE
  1786. *
  1787. *  -SURFACE INTERIEURE: SOUMISE AU MOUVEMENT CHAMP LIBRE
  1788. *  -LE RESTE DE LA ZONE INTERMEDIARE Y COMPRIS LES NOEUDS INTERMEDIAIRES
  1789. *   SI ELEMENTS CU20 OU PR15): MOUVEMENT BLOQUE (ON NE MET RIEN)
  1790. *-----------------------------------------------------------------------
  1791.  
  1792.  FOR_I  = 'FORC' FDIR 1.0 POI_INT.1             ;
  1793.  CHD_I  = 'CHAR' FOR_I D_INT.1                  ;
  1794.  CHV_I  = 'CHAR' FOR_I V_INT.1                  ;
  1795.  CHA_I  = 'CHAR' FOR_I A_INT.1                  ;
  1796.  
  1797.  J = 1                                          ;
  1798. 'REPE' BINT3 (NB_INT - 1)                       ;
  1799.  J = J + 1                                      ;
  1800.    FOR_I  = 'FORC' FDIR 1.0  POI_INT.J          ;
  1801.    CHD_I  = CHD_I 'ET' ('CHAR' FOR_I D_INT.J )  ;
  1802.    CHV_I  = CHV_I 'ET' ('CHAR' FOR_I V_INT.J )  ;
  1803.    CHA_I  = CHA_I 'ET' ('CHAR' FOR_I A_INT.J )  ;
  1804. 'FIN' BINT3                                     ;
  1805.  
  1806.  CHD_I = 'CHAN' 'TABL' CHD_I 'MECA' ;
  1807.  CHV_I = 'CHAN' 'TABL' CHV_I 'MECA' ;
  1808.  CHA_I = 'CHAN' 'TABL' CHA_I 'MECA' ;
  1809.  
  1810.  TAB_T TABD_I = 'EXTR' CHD_I 'MECA' 'TABL' ;
  1811.  TAB_T TABV_I = 'EXTR' CHV_I 'MECA' 'TABL' ;
  1812.  TAB_T TABA_I = 'EXTR' CHA_I 'MECA' 'TABL' ;
  1813.  
  1814.  TABF_I = 'TABLE' ;
  1815.  N_T    = 'DIME' TAB_T ;
  1816.  
  1817.  I = 0 ;
  1818. 'REPE' BINT4 N_T ;
  1819.  CHD_I = 'EXCO' LM1 (TABD_I.I) LM2 ;
  1820.  CHV_I = 'EXCO' LM1 (TABV_I.I) LM2 ;
  1821.  CHA_I = 'EXCO' LM1 (TABA_I.I) LM2 ;
  1822.  
  1823.  CHD_I = CHD_I * BLK_RIG ;
  1824.  CHV_I = CHV_I * BLK_AMO ;
  1825.  CHA_I = CHA_I * BLK_MAS ;
  1826.  
  1827.  TABF_I.I  = CHD_I + CHV_I + CHA_I ;
  1828.  I = I + 1  ;
  1829. 'FIN' BINT4 ;
  1830.  
  1831.  CHA2 = 'CHAR' 'MECA' TAB_T TABF_I ;
  1832.  
  1833. *  Sortie des resultats
  1834.  
  1835.  TABRESU = 'TABLE'                ;
  1836.  TABRESU.'CHAR'    = CHA1 ET CHA2 ;
  1837.  TABRESU.'BLK_MAI' = BLK_MAI      ;
  1838.  TABRESU.'BLK_MOD' = BLK_MOD      ;
  1839.  TABRESU.'BLK_MAT' = BLK_MAT      ;
  1840.  TABRESU.'BLK_MAS' = BLK_MAS      ;
  1841.  TABRESU.'BLK_RIG' = BLK_RIG      ;
  1842.  TABRESU.'BLK_AMO' = BLK_AMO      ;
  1843.  
  1844.  TABRESU.'ACCE'  = GAMMA          ;
  1845.  TABRESU.'PAS'   = PAS_TEM        ;
  1846.  TABRESU.'FCDYN' = 1./(4*PAS_TEM) ;
  1847.  
  1848. FINSI ;
  1849. FINSI ;
  1850. *------------FIN DE MODIFICATION 01/10/2015. WANG)---------------------
  1851.  
  1852.  
  1853. 'SAUT' 5 'LIGNE'                                                   ;
  1854. 'MESS' '*********************************************************' ;
  1855. 'MESS' '*                                                       *' ;
  1856. 'MESS' '*        FIN DE LA PROCEDURE DE LA DECONVOLUTION        *' ;
  1857. 'MESS' '*                                                       *' ;
  1858. 'MESS' '*********************************************************' ;
  1859. 'SAUT' 5 'LIGNE'                                                   ;
  1860.  
  1861.  
  1862. 'FINPROC' TABRESU ;
  1863. ************************************************************************
  1864.  
  1865.  

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