Télécharger dyna13.dgibi

Retour à la liste

Numérotation des lignes :

  1. * fichier : dyna13.dgibi
  2. ************************************************************************
  3. ************************************************************************
  4. **************************************
  5. *
  6. * TEST POUR LA SOUS-STRUCTURATION
  7. *
  8. * UTILISATION DE BASE (OSCAR...)
  9. *
  10. * ASSEMBLAGE DE 2 PLAQUES
  11. *
  12. * D. COMBESCURE 30/09/2005
  13. **************************************
  14. *
  15. *GRAPH = 'O';
  16. GRAPH = 'N';
  17. *
  18. opti dime 2 elem qua4 mode plan cont;
  19. OPTI EPSI LINEAIRE;
  20. *
  21. nmode = 1;
  22. *
  23. nele = 1;
  24. q0 = 0. 0. ; vz = 0. 1;
  25. q1 = 5. 0. ; q2 = 10. 0. ;
  26. q1b = q1 plus (0. 0.);
  27. q3 = q0 plus vz;
  28. q4 = q1 plus vz;
  29. q5 = q2 plus vz;
  30. q4b = q4 plus (0. 0.);
  31. lig1 = d (5*nele) q0 q1 ;
  32. sur1 = tran nele lig1 vz;
  33. ligi0 = d nele q0 q3;
  34. ligi1 = d nele q1 q4;
  35. elim 0.0001 (lig1 et sur1 et ligi1 et ligi0);
  36.  
  37. lig2 = d (5*nele) q1b q2 ;
  38. sur2 = coul (tran nele lig2 vz) rouge;
  39. ligi1b = d nele q1b q4b;
  40. ligi2 = d nele q2 q5;
  41. elim 0.0001 (lig2 et sur2 et ligi2 et ligi1b);
  42. *
  43. si (ega graph 'O');
  44. trac (sur1 et sur2);
  45. finsi;
  46. *
  47. * Plaque 1
  48. *
  49. MO1= MODE SUR1 MECANIQUE ELASTIQUE QUA4 ;
  50. MATE1=MATE MO1 YOUNG 2.E5 NU 0.3 RHO 7.800D-3;
  51. RIGPL1 =RIGI MATE1 MO1 ;
  52. MASPLA1=MASS MATE1 MO1 ;
  53. BLOQ1=(BLOQ DEPL ligi0);
  54. RIGPLA1 = RIGPL1 et BLOQ1 ;
  55. SPLA1=STRU RIGPLA1 MASPLA1;
  56. *ELM1=ELST (q1) SPLA1;
  57. *ELM12=ELST (q4) SPLA1;
  58. Meshi1 = ligi1;
  59. BLOQ1L=BLOQ DEPL Meshi1;
  60. *
  61. * Plaque 2
  62. *
  63. MO2= MODE SUR2 MECANIQUE ELASTIQUE QUA4 ;
  64. MATE2=MATE MO2 YOUNG 2.E5 NU 0.3 RHO 7.80d-3;
  65. RIGPL2 =RIGI MATE2 MO2 ;
  66. MASPLA2=MASS MATE2 MO2 ;
  67. BLOQ2=BLOQ DEPL ligi2;
  68. Meshi2 = ligi1b;
  69. BLOQ2L=BLOQ DEPL Meshi2;
  70. RIGPLA2 = RIGPL2 et BLOQ2 ;
  71. * Liaison
  72. LIUX = RELA ux Meshi1 - ux Meshi2 ;
  73. *LIUX2=RELA UX q4 - UX q4b;
  74. LIUY=RELA UY Meshi1 - UY Meshi2;
  75. *LIUY2=RELA UY q4 - UY q4b;
  76. ENCL = (LIUX et LIUY);
  77. nmodet = 4;
  78. MODPLAT=VIBR PROC (prog 0.) (lect nmodet)
  79. (RIGPLA1 et RIGPLA2 et ENCL)
  80. (MASPLA1 et MASPLA2);
  81. MODPLAT = traduire MODPLAT;
  82. *
  83. si (ega graph 'O');
  84. repeter lab1 nmodet;
  85. list (MODPLAT. &lab1 . FREQUENCE);
  86. trac (defo (sur1 et sur2) (MODPLAT. &lab1 . DEFORMEE_MODALE))
  87. titre 'Structure complete';
  88. fin lab1;
  89. finsi;
  90. *
  91. *
  92. *=======================================================
  93. *
  94. *
  95. * Test de l ancienne méthode CAST3M
  96. *
  97. *
  98. *
  99. *======================================================
  100. *
  101. * Mode libre libre
  102. *=======================================================
  103. MODPLA1=VIBR PROC (prog 0.) (lect nmode)
  104. RIGPLA1 MASPLA1;
  105. * (RIGPLA1 et BLOQ1L) MASPLA1;
  106. MODPLA2=VIBR PROC (prog 0.) (lect nmode)
  107. RIGPLA2 MASPLA2;
  108. * (RIGPLA2 et BLOQ2L) MASPLA2;
  109. *
  110. MODPLA1T = traduire MODPLA1;
  111. MODPLA2T = traduire MODPLA2;
  112. si (ega graph 'O');
  113. trac (defo sur1 (MODPLA1T.1 . DEFORMEE_MODALE))
  114. titre 'Modes propres (CL SS libres)';
  115. trac (defo sur2 (MODPLA2T.1 . DEFORMEE_MODALE))
  116. titre 'Modes propres (CL SS libres)';
  117. finsi;
  118. *
  119. SPLA1=STRU RIGPLA1 MASPLA1 ;
  120. SPLA2=STRU RIGPLA2 MASPLA2 ;
  121. nlab1 = nbno Meshi1;
  122. REPETER LAB1 nlab1;
  123. i = &lab1;
  124. ELM1=ELST (Meshi1 poin i) SPLA1;
  125. ELM2=ELST (Meshi2 poin i) SPLA2;
  126. SI (EGA i 1);
  127. ENCL2 = (RELA ELM1 'UY' 1. ELM2 'UY' -1. )
  128. et (RELA ELM1 'UX' 1. ELM2 'UX' -1. );
  129. SINON;
  130. ENCL2 = ENCL2 et (RELA ELM1 'UY' 1. ELM2 'UY' -1. )
  131. et (RELA ELM1 'UX' 1. ELM2 'UX' -1. );
  132. FINSI;
  133. FIN LAB1;
  134. SOL1 = SOLS ENCL2 SPLA1;
  135. SOL1D = traduire SOL1;
  136. SOL2 = SOLS ENCL2 SPLA2;
  137. SOL2D = traduire SOL2;
  138. si (ega graph 'O');
  139. repeter lab4 ((dime sol1d) - 2);
  140. trac (defo (sur1 et sur2)
  141. ((sol1d. &lab4 .DEFORMEE_MODALE) et
  142. (sol2d. &lab4 .DEFORMEE_MODALE)))
  143. titre 'Modes statiques (CL SS libres)';
  144. fin lab4;
  145. finsi;
  146. BASE1=BASE SPLA1 MODPLA1 ENCL2 SOL1;
  147. BASE2=BASE SPLA2 MODPLA2 ENCL2 SOL2;
  148. *
  149. BASET= BASE1 ET BASE2 ;
  150. RIGMO2 = RIGI BASET ;
  151. MASMO2 = MASS BASET ;
  152. SOL2= VIBR PROC (PROG 0.65) (lect 4) RIGMO2 MASMO2;
  153. soll = synt sol2 BASET;
  154. solb = traduire soll;
  155. *
  156. si (ega graph 'O');
  157. repeter lab2 ((dime solb) - 2);
  158. dp1_t1 = solb. &lab2 . DEFORMEE_MODALE;
  159. frq1 = solb. &lab2 . FREQUENCE;
  160. list frq1;
  161. trac (defo (SUR1 et SUR2) dp1_t1)
  162. titre 'Modes libres';
  163. fin lab2;
  164. finsi;
  165. *
  166. dep11 = MODPLA1T. 1 . DEFORMEE_MODALE;
  167. dep21 = MODPLA2T. 1 . DEFORMEE_MODALE;
  168. depsta11 = SOL1D. 1 . DEFORMEE_MODALE;
  169. depsta12 = SOL1D. 2 . DEFORMEE_MODALE;
  170. depsta13 = SOL1D. 3 . DEFORMEE_MODALE;
  171. depsta14 = SOL1D. 4 . DEFORMEE_MODALE;
  172. *
  173. *===========================================
  174. * Mode bloqué-bloqué
  175. *===========================================
  176. BLOQUX1 = TABLE;
  177. BLOQUX2 = TABLE;
  178. BLOQUY1 = TABLE;
  179. BLOQUY2 = TABLE;
  180. nlab1 = nbno Meshi1;
  181. REPETER LAB1 nlab1;
  182. i = &lab1;
  183. BLOQUX1.i = BLOQ (Meshi1 poin i) UX;
  184. BLOQUX2.i = BLOQ (Meshi2 poin i) UX;
  185. BLOQUY1.i = BLOQ (Meshi1 poin i) UY;
  186. BLOQUY2.i = BLOQ (Meshi2 poin i) UY;
  187. si (i ega 1);
  188. BLOQ1L = (BLOQUX1.i) et (BLOQUY1.i);
  189. BLOQ2L = (BLOQUX2.i) et (BLOQUY2.i);
  190. sinon;
  191. BLOQ1L = BLOQ1L et (BLOQUX1.i) et (BLOQUY1.i);
  192. BLOQ2L = BLOQ2L et (BLOQUX2.i) et (BLOQUY2.i);
  193. finsi;
  194. FIN LAB1;
  195. *
  196. MODPLA1=VIBR PROC (prog 0.) (lect nmode)
  197. (RIGPLA1 et BLOQ1L) MASPLA1;
  198. * (RIGPLA1) MASPLA1;
  199. MODPLA2=VIBR PROC (prog 0.) (lect nmode)
  200. (RIGPLA2 et BLOQ2L) MASPLA2;
  201. * (RIGPLA2) MASPLA2;
  202. *
  203. si (ega graph 'O');
  204. repeter lab3 nmode;
  205. MODPLA1T = traduire MODPLA1;
  206. trac (defo sur1 (MODPLA1T.&lab3 . DEFORMEE_MODALE))
  207. titre 'Modes propres (CL SS bloques)';
  208. MODPLA2T = traduire MODPLA2;
  209. trac (defo sur2 (MODPLA2T.&lab3 . DEFORMEE_MODALE))
  210. titre 'Modes propres (CL SS bloques)';
  211. fin lab3;
  212. finsi;
  213. *
  214. SPLA1=STRU (RIGPLA1 et BLOQ1L) MASPLA1 ;
  215. SPLA2=STRU (RIGPLA2 et BLOQ2L) MASPLA2 ;
  216. *
  217. nlab1 = nbno Meshi1;
  218. REPETER LAB1 nlab1;
  219. i = &lab1;
  220. CLUX1=CLST SPLA1 (BLOQUX1.i);
  221. CLUX2=CLST SPLA2 (BLOQUX2.i);
  222. CLUY1=CLST SPLA1 (BLOQUY1.i);
  223. CLUY2=CLST SPLA2 (BLOQUY2.i);
  224. SI (EGA i 1);
  225. ENCL3 = (RELA CLUX1 'LX' 1. CLUX2 'LX' 1. )
  226. et (RELA CLUY1 'LX' 1. CLUY2 'LX' 1. );
  227. SINON;
  228. ENCL3 = ENCL3 et (RELA CLUX1 'LX' 1. CLUX2 'LX' 1. )
  229. et (RELA CLUY1 'LX' 1. CLUY2 'LX' 1. );
  230. FINSI;
  231. FIN LAB1;
  232. *
  233. SOL1 = SOLS ENCL3 SPLA1;
  234. SOL1D = traduire SOL1;
  235. si (ega graph 'O');
  236. repeter lab4 ((dime sol1d) - 2);
  237. trac (defo (sur1 et sur2) (sol1d. &lab4 .DEFORMEE_MODALE))
  238. titre 'Modes statiques (CL SS bloques)';
  239. fin lab4;
  240. finsi;
  241. SOL2 = SOLS ENCL3 SPLA2;
  242. SOL2D = traduire SOL2;
  243. *
  244. BASE1=BASE SPLA1 ENCL3 MODPLA1 SOL1 ;
  245. BASE2=BASE SPLA2 ENCL3 MODPLA2 SOL2 ;
  246. *
  247. defsta11 = SOL1D. 1 . DEFORMEE_MODALE ;
  248. *
  249. BASET= BASE1 ET BASE2 ;
  250. RIGMO3 = RIGI BASET ;
  251. MASMO3 = MASS BASET ;
  252. SOL3= VIBR PROC (PROG 0.65) (lect 4) RIGMO3 MASMO3;
  253. soll = synt sol3 BASET;
  254. solb = traduire soll;
  255. si (ega graph 'O');
  256. repeter lab2 ((dime solb) - 2);
  257. dp1_t1 = solb. &lab2 . DEFORMEE_MODALE;
  258. frq1 = solb. &lab2 . FREQUENCE;
  259. list frq1;
  260. trac (defo (SUR1 et SUR2) dp1_t1)
  261. titre 'Modes propres bloques';
  262. fin lab2;
  263. finsi;
  264. *
  265. FIN;
  266. *
  267.  
  268.  
  269.  
  270.  
  271.  
  272.  
  273.  
  274.  
  275.  
  276.  
  277.  
  278.  
  279.  

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