Télécharger iss2D_z.dgibi

Retour à la liste

Numérotation des lignes :

  1. * fichier : iss2D_z.dgibi
  2. ************************************************************************
  3. ************************************************************************
  4. GRAPH='Y';
  5. SAUT PAGE ;
  6. ************************************************************************
  7. * REPONSE SISMIQUE DU SOL EN ABSENCE DE STRUTURE
  8. ************************************************************************
  9. *
  10. * DESCRIPTION DU PROBLEME
  11. * ------------------------
  12. * IL S'AGIT D'UN PROBLEME D'INTERACTION SOL-STRUCTURE.
  13. * EN ABSENCE DE STRUCTURE, IL N'Y PAS D'INTERACTION. ON DOIT RETROUVER
  14. * L'ACCELEROGRAMME IMPOSE A LA SURFACE DU SOL VIA LE PROCESSUS DE
  15. * DECONVOLUTION (PROCEDURES DECONV OU DECONV3D) ET CONVOLUTION
  16. * (PROCEDURE DYNAMIC OU PASAPAS).
  17. * LE CYLINDRE DANS CE JEU DE DONNEES REPRESENTE LE SOL PROCHE TANDIS
  18. * QUE LE SOL LOINTAIN QUI S'ETEND VERS L'INFINI EST REPRESENTE PAR UNE
  19. * FRONTIERE ABSORBANTE COMPOSEE D'AMORTISSEURS VISQUEUX.
  20. ************************************************************************
  21. OPTION TRAC PSC ;
  22. OPTION ECHO 0 ;
  23. *-----------------------------------------------------------------------
  24. * PROCEDURE DE CALCUL DE LA DERIVEE PREMIERE
  25. *-----------------------------------------------------------------------
  26. 'DEBPROC' DERIV1 EV1 ;
  27. X = 'EXTR' EV1 'ABSC' ;
  28. Y = 'EXTR' EV1 'ORDO' ;
  29. NPOIN = 'DIME' X ;
  30. H = ( 'EXTR' X 2 ) - ( 'EXTR' X 1 ) ;
  31.  
  32. X1 = ( 'PROG' 0.0 0.0 ) 'ET' Y ;
  33. X3 = Y 'ET' ( 'PROG' 0.0 0.0 ) ;
  34. DX = ( X3 - X1 ) / ( H * 2.0 ) ;
  35. Y1 = 'ENLE' ( 'ENLE' DX ( NPOIN + 2 ) ) 1 ;
  36. L1 = ( 'PROG' 0.0 ) 'ET' ( 'PROG' ( NPOIN - 2 ) * 1.0 )
  37. 'ET' ( 'PROG' 0.0 ) ;
  38. Y2 = Y1 * L1 ;
  39. EV2 = 'EVOL' 'MANU' X Y2 ;
  40. 'FINPROC' EV2 ;
  41. *-----------------------------------------------------------------------
  42. * SOL SANS STRUCTURE
  43. *-----------------------------------------------------------------------
  44. OPTI DIME 2 MODE FOUR 0 ELEM QUA8 COUL BLEU ;
  45.  
  46. R = 2. / (PI**0.5) ;
  47. NE = 2 ;
  48. NR = 3 ;
  49.  
  50. P0 = 0. 0. ;
  51. P4 = (NR*R) 0. ;
  52.  
  53. SUR = COUL VERT (DROIT (NE*NR) P0 P4) ;
  54. SOL = TRAN (NE*NR) SUR (0 (-1*NR*R)) ;
  55. BOR = COTE 2 SOL ;
  56. FON = COTE 3 SOL ;
  57. AXE = COTE 4 SOL ;
  58.  
  59. ELIM SOL 0.001 ;
  60. TITR 'MAILLAGE SOL (AXISYMETRIQUE)' ;
  61. SI (NEG GRAPH 'N');
  62. TRAC SOL QUAL NCLK;
  63. FINSI;
  64. *-----------------------------------------------------------------------
  65. * MODELE DE SOL
  66. *-----------------------------------------------------------------------
  67. NU1 = 0.3 ;
  68. G1 = 100.E7 ;
  69. E1 = 2*(1 + NU1)*G1 ;
  70.  
  71. MOD_S = MODE SOL MECANIQUE ELASTIQUE ISOTROPE QUA8 ;
  72. MAT_S = MATE MOD_S YOUN E1 NU NU1 RHO 2000 ;
  73. MAS_S = MASS MOD_S MAT_S ;
  74. RIG_S = RIGI MOD_S MAT_S ;
  75. *-----------------------------------------------------------------------
  76. * AMORTISSEMENT DE TYPE RAYLEIGH
  77. *-----------------------------------------------------------------------
  78. F1 = 5.0 ;
  79. F2 = 25.0 ;
  80. KSI_S = 0.05 ;
  81. ALPHA = 4.0 * PI * F1 * F2 / (F1 + F2) ;
  82. BETA = 1.0 / (PI * (F1 + F2)) ;
  83. AMO_S = KSI_S * ((ALPHA * MAS_S) ET (BETA * RIG_S)) ;
  84. *-----------------------------------------------------------------------
  85. * SIGNAL SISMIQUE (SELON SPECTRE PS92 SITE S1)
  86. *-----------------------------------------------------------------------
  87. LFR = PROG 0.1 0.5 1.0 2.5 5. 33. 50. ;
  88. LSP = PROG 0.1 0.5 1.0 2.5 2.5 1.0 1.0 ;
  89. SP0 = EVOL MANU 'FREQ(HZ)' LFR 'ACCE(M/S*S)' LSP ;
  90.  
  91. TAB = TABLE ;
  92. TAB.'MOTIT' = 'SPECTRE PS92 S1 ' ;
  93. TAB.'SEISME'= TABLE ;
  94. TAB.'SEISME'.'SPECTRE' = SP0 ;
  95. TAB.'SEISME'.'AMORT' = 0.05 ;
  96. TAB.'SEISME'.'TYPSP' = 'ACCE' ;
  97. TAB.'SIGNAL' = TABLE ;
  98. TAB.'SIGNAL'.'ENVE' = 'PLATLIN' ;
  99. TAB.'SIGNAL'.'NP' = 8 ;
  100. TAB.'SIGNAL'.'DUREE' = 2.56 ;
  101. TAB.'SIGNAL'.'TDEBUT' = 0.75 ;
  102. TAB.'SIGNAL'.'TFIN' = 1.5 ;
  103. TAB.'NBITER' = 5 ;
  104. TAB.'NBSIGN' = 1 ;
  105. TAB.'NALEAT' = 3 ;
  106. TAB.'FRCOUP' = 49.9 ;
  107. TAB.'OPTSORT' = 'SPECTRE' ;
  108.  
  109. TABSIG = SIGNSYNT FABR TAB ;
  110.  
  111. * MOYENNE ZERO, PLAGE INITIALE ZERO
  112.  
  113. LT0 = EXTR TABSIG.1 ABSC ;
  114. LA0 = EXTR TABSIG.1 ORDO ;
  115. DT = (EXTR LT0 2) - (EXTR LT0 1) ;
  116. NP0 = DIME LT0 ;
  117.  
  118. LT1 = PROG 0. PAS DT 2.3 ;
  119. LA1 = IPOL LT1 LT0 LA0 ;
  120. NP1 = DIME LT1 ;
  121. LA1 = LA1 - (PROG NP1 * ((SOMM LA1) / NP1)) ;
  122.  
  123. LA1 = (PROG 20 * 0.) ET LA1 ;
  124. LT1 = PROG 0. PAS DT NPAS (19 + NP1) ;
  125.  
  126. ACC1 = EVOL MANU 'TEMPS(S)' LT1 'ACCE(M/S*S)' LA1 ;
  127. SI (NEG GRAPH 'N');
  128. DESS ACC1 MIMA NCLK;
  129. FINSI ;
  130. *-----------------------------------------------------------------------
  131. * DECONVOLUTION ET FRONTIERE ABSORBANTE
  132. *-----------------------------------------------------------------------
  133. TAB = TABLE ;
  134. TAB.1 = TABLE ;
  135. TAB.1 .'FRONTIERE' = BOR ;
  136. TAB.1 .'MASSE_VOLUMIQUE' = 2000. ;
  137. TAB.1 .'POISSON' = 0.3 ;
  138. TAB.1 .'YOUNG' = E1 ;
  139. TAB.1 .'AMORTISSEMENT' = 0.05 ;
  140.  
  141. TYP_F = 'LYSMER' ;
  142. FC = 50.0 ;
  143. DIR = 'VERT' ;
  144.  
  145. TABS = DECONV TAB FON MOD_S DIR ACC1 F1 F2 FC TYP_F ;
  146. *-----------------------------------------------------------------------
  147. * TABLE POUR LA PROCEDURE 'DYNAMIC'
  148. *-----------------------------------------------------------------------
  149. CH_DEPI = MANU CHPO 3 SOL UR 0.0 UT 0.0 UZ 0.0 ;
  150. CH_VITI = MANU CHPO 3 SOL UR 0.0 UT 0.0 UZ 0.0 ;
  151.  
  152. TAB_DYN = TABLE ;
  153. TAB_DYN.'DEPL' = CH_DEPI ;
  154. TAB_DYN.'VITE' = CH_VITI ;
  155. TAB_DYN.'RIGI' = RIG_S ;
  156. TAB_DYN.'MASS' = MAS_S ;
  157. TAB_DYN.'AMOR' = AMO_S ET TABS.'AMOR' ;
  158. TAB_DYN.'CHAR' = TABS.'CHAR' ;
  159. TAB_DYN.'FREQ' = TABS.'FCDYN' ;
  160. DT = TABS.'PAS' ;
  161. NB_PAS = ENTI (2.51 / DT) ;
  162. TAB_DYN.'DEBU' = 0.0 ;
  163. TAB_DYN.'INST' = PROG 0.0 PAS DT NPAS (NB_PAS - 1) ;
  164.  
  165. TDYNA = DYNAMIC TAB_DYN ;
  166. *-----------------------------------------------------------------------
  167. * POST-TRAITEMENT
  168. *-----------------------------------------------------------------------
  169. POI = TABLE ;
  170. POI.1 = P0 ;
  171. POI.2 = P4 ;
  172.  
  173. NP = DIME POI ;
  174.  
  175. L_TEM = PROG NB_PAS * 0.0 ;
  176. L_D = TABLE ;
  177. L_V = TABLE ;
  178.  
  179. I = 1 ;
  180. REPE B1 NP ;
  181. L_V.I = PROG NB_PAS * 0.0 ;
  182.  
  183. I = I + 1 ;
  184. FIN B1 ;
  185.  
  186. I_TEM = 0 ;
  187. REPETER B2 NB_PAS ;
  188. I_TEM = I_TEM + 1 ;
  189. TEM_I = TDYNA.I_TEM.'TEMP' ;
  190. REMP L_TEM I_TEM TEM_I ;
  191. CHV_I = TDYNA.I_TEM.'VITE' ;
  192.  
  193. I = 1 ;
  194. REPE B2_1 NP ;
  195. REMP L_V.I I_TEM (EXTR CHV_I POI.I UZ) ;
  196. I = I + 1 ;
  197. FIN B2_1 ;
  198.  
  199. FIN B2 ;
  200.  
  201. * REPONSE EN ACCELERATION
  202.  
  203. LT1 = EXTR ACC1 ABSC ;
  204. EV_VIT = TABLE ;
  205. EV_ACC = TABLE ;
  206.  
  207. I = 1 ;
  208. REPE B3 NP ;
  209. EV_VIT.I = EVOL MANU L_TEM L_V.I ;
  210. EV_ACC.I = DERIV1 EV_VIT.I ;
  211. LT2 = EXTR EV_ACC.I ABSC ;
  212. LA2 = EXTR EV_ACC.I ORDO ;
  213. LA2 = IPOL LT1 LT2 LA2 ;
  214. EV_ACC.I = EVOL MANU 'TEMPS(S)' LT1 'ACCE(M/S*S)' LA2 ;
  215. I = I + 1 ;
  216. FIN B3 ;
  217.  
  218. * COMPARAISON AVEC L'ACCELEROGRAMME IMPOSE
  219.  
  220. SI (NEG GRAPH 'N');
  221. DESS ((COUL ROUG EV_ACC.1) ET (COUL VERT EV_ACC.2) ET ACC1) NCLK
  222. TITR 'ACCELEROGRAMMES IMPOSE (BLEU) ET CALCULES (ROUGE ET VERT)' MIMA ;
  223. FINSI ;
  224.  
  225. ERR1 = MAXI ABS (EXTR (EV_ACC.1 - ACC1) ORDO) ;
  226. ERR2 = MAXI ABS (EXTR (EV_ACC.2 - ACC1) ORDO) ;
  227.  
  228. SI ((MAXI (PROG ERR1 ERR2)) > 0.01);
  229. ERRE 5;
  230. FINSI;
  231.  
  232. FIN ;
  233.  
  234.  
  235.  
  236.  
  237.  
  238.  
  239.  
  240.  
  241.  
  242.  

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