Télécharger sissi.dgibi

Retour à la liste

Numérotation des lignes :

  1. * fichier : sissi.dgibi
  2. ************************************************************************
  3. ************************************************************************
  4.  
  5. *******************************************
  6. * Test sissi.dgibi: jeux de données *
  7. * --------------------------------- *
  8. *******************************************
  9.  
  10. * SI GRAPH = N PAS DE GRAPHIQUE AFFICHE
  11. * SINON SI GRAPH DIFFERENT DE N TOUS
  12. * LES GRAPHIQUES SONT AFFICHES
  13.  
  14. GRAPH = 'N' ;
  15.  
  16. SAUT PAGE;
  17. SI (NEG GRAPH 'N') ;
  18. OPTI ECHO 1 ;
  19. OPTI TRAC PSC ;
  20. SINO ;
  21. OPTI ECHO 0 ;
  22. FINSI ;
  23.  
  24. SAUT PAGE;
  25. *******************************************
  26. * *
  27. * Mots-clés : Vibrations, calcul modal, *
  28. * poutre, seisme *
  29. * *
  30. * SISSI *
  31. * *
  32. * CAS TEST DU 91/06/19 PROVENANCE : PLAF *
  33. * CAS TEST DU 91/06/18 PROVENANCE : PLAF *
  34. * *
  35. * test de la procedure SISSIB *
  36. * *
  37. * 1 poutre encastree, 20 elements finis, *
  38. * IY different de IZ, un spectre *
  39. * d'oscillateurs en ACCE direction du *
  40. * seisme : X *
  41. * *
  42. *******************************************
  43.  
  44. TEMPS ;
  45. OPTI DIME 3 ELEM 'SEG2' ;
  46. OPTI EPSI 'LINEAIRE';
  47. *
  48. C1 = 0. 0. 0. ;
  49. C2 = 0. 0. 5. ;
  50. C3 = 0. 0. 10. ;
  51. L1 = DROITE 10 C1 C2 ;
  52. L2 = DROITE 10 C2 C3 ;
  53. LI = L1 ET L2 ;
  54.  
  55.  
  56. SI (NEG GRAPH 'N');
  57. TRAC 'QUAL' LI;
  58. FINSI;
  59.  
  60. MOD1 = MODE LI MECANIQUE ELASTIQUE POUT ;
  61.  
  62. CH_MAT = MATE MOD1 YOUNG 2.E11 NU 0.3 RHO 7800. ;
  63. CH_CAR = CARA MOD1 SECT 0.25 INRY 0.006 INRZ 0.004 TORS 0.01 VECT ( 0. 1. 0. ) ;
  64. CH_MAT=CH_MAT et CH_CAR;
  65.  
  66. RIG1 = RIGI CH_MAT MOD1 ;
  67. ENC1 = BLOQ C1 DEPL ROTA ;
  68. RIGFI= RIG1 ET ENC1 ;
  69.  
  70. MAS1 = MASS CH_MAT MOD1 ;
  71.  
  72. * Calcul des premiers modes
  73.  
  74. L_FREQ = PROG 3.58 4.38 24.09 25.15 ;
  75. MODE_POU = VIBR PROCHE L_FREQ RIGFI MAS1 ;
  76.  
  77. * Calcul des contraintes modales
  78. MODE_POU = SIGSOL MOD1 CH_MAT MODE_POU ;
  79.  
  80. * Calcul des reactions modales
  81. MODE_POU = REAC ENC1 MODE_POU ;
  82.  
  83. ITAB2 = MODE_POU . 'MODES' ;
  84. NB_MODE = ( DIME ITAB2 ) - 2 ;
  85.  
  86. * Definition du spectre acceleration
  87. LIS_TEMP = PROG 0. PAS 1.E-2 0.50 ;
  88. LIS_ACCE = PROG 0. 0.5 1. 0.5 0. -0.5 -1. -0.5 0. 0.5 1. 0.5 0. -0.5 -1. -0.5 0. 0.5 1. 0.5 0. -0.5 -1. -0.5 0. 0.5 1. 0.5 0. -0.5 -1. -0.5 0. 0.5 1. 0.5 0. -0.35 -0.7 -0.35 0. 0.25 0.5 0.25 0. -0.15 -0.3 -0.15 0. 0.1 0. ;
  89. LIS_CHAR = EVOL MANU LIS_TEMP LIS_ACCE ;
  90.  
  91. SI (NEG GRAPH 'N');
  92. TITR 'Chargement en Acceleration';
  93. DESS LIS_CHAR;
  94. FINSI;
  95.  
  96. LIS_FREQ = PROG 1. PAS 2. 400. ;
  97. LIS_AMOR = PROG 2. 5. 10. 15. ;
  98. LIS_AMOR = LIS_AMOR * 1.e-2 ;
  99. SEISME_X = SPO LIS_CHAR 'AMOR' LIS_AMOR 'FREQ' LIS_FREQ 'ACCE' ;
  100.  
  101. BAS_AMOR = PROG 4. 6. 7. 9. ;
  102. BAS_AMOR = BAS_AMOR * 1.e-2 ;
  103.  
  104. TAB1 = TABL ;
  105.  
  106. TAB1 . 'STRUCTURE' = MODE_POU ;
  107. TAB1 . 'AMORTISSEMENT' = BAS_AMOR ;
  108.  
  109. TAB3 = 'TABLE' 'EXCITATION' ;
  110. TAB1 . 'EXCITATION' = TAB3 ;
  111.  
  112. TAB3 . 1 = 'TABLE' ;
  113. TAB3 . 1 . 'DIRECTION' = 'X' ;
  114. TAB3 . 1 . 'SPECTRE' = SEISME_X ;
  115. TAB3 . 1 . 'AMORTISSEMENT' = LIS_AMOR ;
  116.  
  117. TAB1 . 'RECOMBINAISON_MODES' = 'SRSS';
  118. TAB1 . 'RECOMBINAISON_DIRECTIONS' = 'QUADRATIQUE' ;
  119.  
  120. TAB4 = 'TABLE' 'SORTIES' ;
  121. TAB1 . 'SORTIES' = TAB4 ;
  122. TAB4 . 'DOMAINE' = MOD1 ;
  123. TAB4 . 'DEPLACEMENTS' = VRAI ;
  124. TAB4 . 'CONTRAINTES' = VRAI ;
  125. TAB4 . 'ACCELERATIONS' = VRAI ;
  126.  
  127. TAB2 = SISSIB TAB1 ;
  128.  
  129. MESS ' Chpoint de deplacement du point C3 issu de SISSIB :' ;
  130. MESS ' ' ;
  131. CHP_DEP = TAB2 . 'X' . 'DEPLACEMENTS' ;
  132. LIST ( REDU CHP_DEP C3 ) ;
  133. MESS ' Chpoint d acceleration du point C3 issu de SISSIB :' ;
  134. MESS ' ' ;
  135. CHP_ACC = TAB2 . 'X' . 'ACCELERATIONS' ;
  136. LIST ( REDU CHP_ACC C3 ) ;
  137.  
  138. CHE_CONT = TAB2 . 'X' . 'CONTRAINTES' ;
  139. MAX_CONT = MAXI CHE_CONT ;
  140. MESS ' Valeur maximale des contraintes issue de SISSIB =' MAX_CONT ;
  141. MESS ' ' ;
  142.  
  143. *-----------------------------------------*
  144. * *
  145. * Verification de la procedure *
  146. * *
  147. *-----------------------------------------*
  148. * Verification du deplacement
  149. * Verification de l'acceleration
  150. * Verification de la contrainte maximale
  151.  
  152. * ----- calcul de S(N,X,B) -----
  153. TSEIS = TABL ;
  154. NB_AMOR = 'DIME' LIS_AMOR ;
  155. I_MODE = 0 ;
  156. REPETER BOUC1 NB_MODE ;
  157. I_MODE = I_MODE + 1 ;
  158. ITAB3 = ITAB2 . I_MODE ;
  159. F_N = ITAB3 . 'FREQUENCE' ;
  160. LOG_F_N = LOG F_N ;
  161. B_N = EXTRAIRE BAS_AMOR I_MODE ;
  162.  
  163. P_SPEC = 'PROG' ;
  164. I_AMOR = 0 ;
  165. REPETER BOUC10 NB_AMOR ;
  166. I_AMOR = I_AMOR '+' 1 ;
  167. EVOLS1 = 'EXTR' SEISME_X 'COUR' I_AMOR ;
  168. LISABS1 = 'EXTR' EVOLS1 'ABSC' ;
  169. LISORD1 = 'EXTR' EVOLS1 'ORDO' ;
  170. LOG_ABS1 = 'LOG' LISABS1 ;
  171. LOG_ORD1 = 'LOG' LISORD1 ;
  172. S1 = IPOL LOG_F_N LOG_ABS1 LOG_ORD1 ;
  173. S_1 = 'EXP' S1 ;
  174.  
  175. P_SPEC = P_SPEC ET ( 'PROG' S_1 ) ;
  176. FIN BOUC10 ;
  177.  
  178. S_N = 'IPOL' B_N LIS_AMOR P_SPEC ;
  179.  
  180. TSEIS . I_MODE = S_N ;
  181. FIN BOUC1 ;
  182.  
  183. I_MODE = 0 ;
  184. ZTRON = 'MANU' 'CHPO' LI 1 'UX' 1. ;
  185. *
  186. REPETER BOUC2 NB_MODE ;
  187. I_MODE = I_MODE + 1 ;
  188. S_I = TSEIS . I_MODE ;
  189. ITAB3 = ITAB2 . I_MODE ;
  190. F_I = ITAB3 .'FREQUENCE' ;
  191. D_I = ITAB3 .'DEFORMEE_MODALE' ;
  192. C_I = ITAB3 .'CONTRAINTE_MODALE' ;
  193. ITAB4=ITAB3 .'DEPLACEMENTS_GENERALISES';
  194. Q_I = ITAB4 . 1 ;
  195. M_I = ITAB3 . 'MASSE_GENERALISEE' ;
  196. W_I = 2.0 * PI * F_I ;
  197. QSM = Q_I / M_I ;
  198. COEF_1 = QSM / ( W_I * W_I ) ;
  199. DEPL_I = COEF_1 * S_I * D_I ;
  200. ACCE_I = QSM * S_I * D_I ;
  201. ZTRON = ZTRON - ( QSM * D_I ) ;
  202. CONT_I = COEF_1 * S_I * C_I ;
  203. 'SI' ( I_MODE 'EGA' 1 ) ;
  204. DEPL_T = DEPL_I ** 2 ;
  205. ACCE_T = ACCE_I ** 2 ;
  206. CONT_T = CONT_I ** 2 ;
  207. 'SINON' ;
  208. DEPL_T = DEPL_T + ( DEPL_I ** 2 ) ;
  209. ACCE_T = ACCE_T + ( ACCE_I ** 2 ) ;
  210. CONT_T = CONT_T + ( CONT_I ** 2 ) ;
  211. 'FINSI' ;
  212. FIN BOUC2 ;
  213.  
  214. LISORD1 = 'EXTR' SEISME_X 'ORDO' 1 ;
  215. NVAL_1 = 'DIME' LISORD1 ;
  216. GAMMA0 = 'EXTR' LISORD1 NVAL_1 ;
  217. ACCE_T = ACCE_T + ( ( GAMMA0 * GAMMA0 ) * ( ZTRON ** 2 ) );
  218.  
  219. DEPL_T = DEPL_T ** 0.5 ;
  220. ACCE_T = ACCE_T ** 0.5 ;
  221. CONT_T = CONT_T ** 0.5 ;
  222.  
  223. D_SOM_X = 'EXTR' DEPL_T C3 'UX' ;
  224. A_SOM_X = 'EXTR' ACCE_T C3 'UX' ;
  225. C_MAX = 'MAXI' CONT_T ;
  226.  
  227.  
  228. DEP = EXTR CHP_DEP C3 'UX' ;
  229. ACC = EXTR CHP_ACC C3 'UX' ;
  230.  
  231.  
  232. MESS ' deplacement du point C3 :' ;
  233. LIST ( 'REDU' DEPL_T C3 ) ;
  234. SAUTER 2 LIGNES ;
  235. MESS ' acceleration du point C3 :' ;
  236. LIST ( 'REDU' ACCE_T C3 ) ;
  237. SAUTER 2 LIGNES ;
  238. MESS ' contrainte maximale =' C_MAX ;
  239. SAUTER 2 LIGNES ;
  240.  
  241. TEMPS ;
  242.  
  243. * Code de bon fonctionnement
  244. REF_DEP = D_SOM_X ;
  245. REF_ACC = A_SOM_X ;
  246. REF_CON = C_MAX ;
  247. RES1 = 100 * (ABS ( ( DEP - REF_DEP ) / REF_DEP ));
  248. RES2 = 100 * (ABS ( ( ACC - REF_ACC ) / REF_ACC ));
  249. RES3 = 100 * (ABS ( ( MAX_CONT - REF_CON ) / REF_CON ));
  250.  
  251. SAUTER 2 LIGNES ;
  252.  
  253. MESS 'deplacement theorique :' D_SOM_X 'm';
  254. MESS 'deplacement calculee :' DEP 'm';
  255. MESS ' Soit un ecart de : ' RES1 '%';
  256. SAUTER 1 LIGNES ;
  257.  
  258. MESS 'acceleration theorique:' A_SOM_X 'm';
  259. MESS 'acceleration calculee :' ACC 'm';
  260. MESS ' Soit un ecart de : ' RES2 '%';
  261. SAUTER 1 LIGNES ;
  262. SI ( RES1 <EG 1. ) ;
  263. ERRE 0 ;
  264. SINON ;
  265. ERRE 5 ;
  266. FINSI ;
  267. SI ( RES2 <EG 1. ) ;
  268. ERRE 0 ;
  269. SINON ;
  270. ERRE 5 ;
  271. FINSI ;
  272. SI ( RES3 <EG 1. ) ;
  273. ERRE 0 ;
  274. SINON ;
  275. ERRE 5 ;
  276. FINSI ;
  277.  
  278. FIN ;
  279.  
  280.  
  281.  
  282.  
  283.  
  284.  
  285.  
  286.  
  287.  

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