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. *........... MISE EN COMMENTAIRE
  93. *........... POUR PREPARATION SUPPRESSION
  94. *........... OBJETS SOLUTIONS AND CO ...
  95. * *=======================================================
  96. * *
  97. * *
  98. * * Test de l ancienne méthode CAST3M
  99. * *
  100. * *
  101. * *
  102. * *======================================================
  103. * *
  104. * * Mode libre libre
  105. * *=======================================================
  106. * MODPLA1=VIBR PROC (prog 0.) (lect nmode)
  107. * RIGPLA1 MASPLA1;
  108. * * (RIGPLA1 et BLOQ1L) MASPLA1;
  109. * MODPLA2=VIBR PROC (prog 0.) (lect nmode)
  110. * RIGPLA2 MASPLA2;
  111. * * (RIGPLA2 et BLOQ2L) MASPLA2;
  112. * *
  113. * MODPLA1T = traduire MODPLA1;
  114. * MODPLA2T = traduire MODPLA2;
  115. * si (ega graph 'O');
  116. * trac (defo sur1 (MODPLA1T.1 . DEFORMEE_MODALE))
  117. * titre 'Modes propres (CL SS libres)';
  118. * trac (defo sur2 (MODPLA2T.1 . DEFORMEE_MODALE))
  119. * titre 'Modes propres (CL SS libres)';
  120. * finsi;
  121. * *
  122. * SPLA1=STRU RIGPLA1 MASPLA1 ;
  123. * SPLA2=STRU RIGPLA2 MASPLA2 ;
  124. * nlab1 = nbno Meshi1;
  125. * REPETER LAB1 nlab1;
  126. * i = &lab1;
  127. * ELM1=ELST (Meshi1 poin i) SPLA1;
  128. * ELM2=ELST (Meshi2 poin i) SPLA2;
  129. * SI (EGA i 1);
  130. * ENCL2 = (RELA ELM1 'UY' 1. ELM2 'UY' -1. )
  131. * et (RELA ELM1 'UX' 1. ELM2 'UX' -1. );
  132. * SINON;
  133. * ENCL2 = ENCL2 et (RELA ELM1 'UY' 1. ELM2 'UY' -1. )
  134. * et (RELA ELM1 'UX' 1. ELM2 'UX' -1. );
  135. * FINSI;
  136. * FIN LAB1;
  137. * SOL1 = SOLS ENCL2 SPLA1;
  138. * SOL1D = traduire SOL1;
  139. * SOL2 = SOLS ENCL2 SPLA2;
  140. * SOL2D = traduire SOL2;
  141. * si (ega graph 'O');
  142. * repeter lab4 ((dime sol1d) - 2);
  143. * trac (defo (sur1 et sur2)
  144. * ((sol1d. &lab4 .DEFORMEE_MODALE) et
  145. * (sol2d. &lab4 .DEFORMEE_MODALE)))
  146. * titre 'Modes statiques (CL SS libres)';
  147. * fin lab4;
  148. * finsi;
  149. * BASE1=BASE SPLA1 MODPLA1 ENCL2 SOL1;
  150. * BASE2=BASE SPLA2 MODPLA2 ENCL2 SOL2;
  151. * *
  152. * BASET= BASE1 ET BASE2 ;
  153. * RIGMO2 = RIGI BASET ;
  154. * MASMO2 = MASS BASET ;
  155. * SOL2= VIBR PROC (PROG 0.65) (lect 4) RIGMO2 MASMO2;
  156. * soll = synt sol2 BASET;
  157. * solb = traduire soll;
  158. * *
  159. * si (ega graph 'O');
  160. * repeter lab2 ((dime solb) - 2);
  161. * dp1_t1 = solb. &lab2 . DEFORMEE_MODALE;
  162. * frq1 = solb. &lab2 . FREQUENCE;
  163. * list frq1;
  164. * trac (defo (SUR1 et SUR2) dp1_t1)
  165. * titre 'Modes libres';
  166. * fin lab2;
  167. * finsi;
  168. * *
  169. * dep11 = MODPLA1T. 1 . DEFORMEE_MODALE;
  170. * dep21 = MODPLA2T. 1 . DEFORMEE_MODALE;
  171. * depsta11 = SOL1D. 1 . DEFORMEE_MODALE;
  172. * depsta12 = SOL1D. 2 . DEFORMEE_MODALE;
  173. * depsta13 = SOL1D. 3 . DEFORMEE_MODALE;
  174. * depsta14 = SOL1D. 4 . DEFORMEE_MODALE;
  175. * *
  176. * *===========================================
  177. * * Mode bloqué-bloqué
  178. * *===========================================
  179. * BLOQUX1 = TABLE;
  180. * BLOQUX2 = TABLE;
  181. * BLOQUY1 = TABLE;
  182. * BLOQUY2 = TABLE;
  183. * nlab1 = nbno Meshi1;
  184. * REPETER LAB1 nlab1;
  185. * i = &lab1;
  186. * BLOQUX1.i = BLOQ (Meshi1 poin i) UX;
  187. * BLOQUX2.i = BLOQ (Meshi2 poin i) UX;
  188. * BLOQUY1.i = BLOQ (Meshi1 poin i) UY;
  189. * BLOQUY2.i = BLOQ (Meshi2 poin i) UY;
  190. * si (i ega 1);
  191. * BLOQ1L = (BLOQUX1.i) et (BLOQUY1.i);
  192. * BLOQ2L = (BLOQUX2.i) et (BLOQUY2.i);
  193. * sinon;
  194. * BLOQ1L = BLOQ1L et (BLOQUX1.i) et (BLOQUY1.i);
  195. * BLOQ2L = BLOQ2L et (BLOQUX2.i) et (BLOQUY2.i);
  196. * finsi;
  197. * FIN LAB1;
  198. * *
  199. * MODPLA1=VIBR PROC (prog 0.) (lect nmode)
  200. * (RIGPLA1 et BLOQ1L) MASPLA1;
  201. * * (RIGPLA1) MASPLA1;
  202. * MODPLA2=VIBR PROC (prog 0.) (lect nmode)
  203. * (RIGPLA2 et BLOQ2L) MASPLA2;
  204. * * (RIGPLA2) MASPLA2;
  205. * *
  206. * si (ega graph 'O');
  207. * repeter lab3 nmode;
  208. * MODPLA1T = traduire MODPLA1;
  209. * trac (defo sur1 (MODPLA1T.&lab3 . DEFORMEE_MODALE))
  210. * titre 'Modes propres (CL SS bloques)';
  211. * MODPLA2T = traduire MODPLA2;
  212. * trac (defo sur2 (MODPLA2T.&lab3 . DEFORMEE_MODALE))
  213. * titre 'Modes propres (CL SS bloques)';
  214. * fin lab3;
  215. * finsi;
  216. * *
  217. * SPLA1=STRU (RIGPLA1 et BLOQ1L) MASPLA1 ;
  218. * SPLA2=STRU (RIGPLA2 et BLOQ2L) MASPLA2 ;
  219. * *
  220. * nlab1 = nbno Meshi1;
  221. * REPETER LAB1 nlab1;
  222. * i = &lab1;
  223. * CLUX1=CLST SPLA1 (BLOQUX1.i);
  224. * CLUX2=CLST SPLA2 (BLOQUX2.i);
  225. * CLUY1=CLST SPLA1 (BLOQUY1.i);
  226. * CLUY2=CLST SPLA2 (BLOQUY2.i);
  227. * SI (EGA i 1);
  228. * ENCL3 = (RELA CLUX1 'LX' 1. CLUX2 'LX' 1. )
  229. * et (RELA CLUY1 'LX' 1. CLUY2 'LX' 1. );
  230. * SINON;
  231. * ENCL3 = ENCL3 et (RELA CLUX1 'LX' 1. CLUX2 'LX' 1. )
  232. * et (RELA CLUY1 'LX' 1. CLUY2 'LX' 1. );
  233. * FINSI;
  234. * FIN LAB1;
  235. * *
  236. * SOL1 = SOLS ENCL3 SPLA1;
  237. * SOL1D = traduire SOL1;
  238. * si (ega graph 'O');
  239. * repeter lab4 ((dime sol1d) - 2);
  240. * trac (defo (sur1 et sur2) (sol1d. &lab4 .DEFORMEE_MODALE))
  241. * titre 'Modes statiques (CL SS bloques)';
  242. * fin lab4;
  243. * finsi;
  244. * SOL2 = SOLS ENCL3 SPLA2;
  245. * SOL2D = traduire SOL2;
  246. * *
  247. * BASE1=BASE SPLA1 ENCL3 MODPLA1 SOL1 ;
  248. * BASE2=BASE SPLA2 ENCL3 MODPLA2 SOL2 ;
  249. * *
  250. * defsta11 = SOL1D. 1 . DEFORMEE_MODALE ;
  251. * *
  252. * BASET= BASE1 ET BASE2 ;
  253. * RIGMO3 = RIGI BASET ;
  254. * MASMO3 = MASS BASET ;
  255. * SOL3= VIBR PROC (PROG 0.65) (lect 4) RIGMO3 MASMO3;
  256. * soll = synt sol3 BASET;
  257. * solb = traduire soll;
  258. * si (ega graph 'O');
  259. * repeter lab2 ((dime solb) - 2);
  260. * dp1_t1 = solb. &lab2 . DEFORMEE_MODALE;
  261. * frq1 = solb. &lab2 . FREQUENCE;
  262. * list frq1;
  263. * trac (defo (SUR1 et SUR2) dp1_t1)
  264. * titre 'Modes propres bloques';
  265. * fin lab2;
  266. * finsi;
  267. * *
  268. *........... FIN DE MISE EN COMMENTAIRE
  269. *........... POUR PREPARATION SUPPRESSION
  270. *........... OBJETS SOLUTIONS AND CO ...
  271. FIN;
  272. *
  273.  
  274.  
  275.  
  276.  
  277.  
  278.  
  279.  
  280.  
  281.  
  282.  
  283.  
  284.  
  285.  
  286.  
  287.  

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