Télécharger la plaquette Cast3M


fr en

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

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