Télécharger dyna5.dgibi

Retour à la liste

Numérotation des lignes :

  1. * fichier : dyna5.dgibi
  2. ************************************************************************
  3. ************************************************************************
  4. * Test Dyna5.dgibi: Jeux de données *
  5. * --------------------------------- *
  6. * *
  7. * CAS TEST DU 01/07/92 PROVENANCE : PHIL
  8. * MODIFICATION D. COMBESCURE 30/09/2005
  9. OPTI ECHO 1;
  10. *GRAPH='O';
  11. GRAPH='N';
  12. SAUT PAGE ;
  13. *******************************************************
  14. *
  15. * DYNA5 - Test de sous-structuration par 3 méthodes
  16. *
  17. * La fréquence propre d'une poutre est égale à 65 Hz.
  18. * On se propose de retrouver cette fréquence par trois
  19. * méthodes de sous-structuration.
  20. *
  21. * Pour cela on décompose la poutre en deux morceaux
  22. * de meme longueur.
  23. *
  24. *
  25. * Q1 poutre 1 Q2
  26. * .--------------..------------.
  27. * Q3 poutre 2 Q4
  28. *
  29. *
  30. * Les trois méthodes de sous-structuration correspon-
  31. * dent aux conditions aux limites de l'extrémité
  32. * commune Q2 et Q3 des deux morceaux de poutres :
  33. * * Q2 libre - Q3 libre
  34. * * Q2 bloque - Q3 bloque
  35. * * Q2 bloque - Q3 libre
  36. *
  37. * A ces trois méthodes correspondent trois manières de
  38. * décrire les liaisons entre sous-structures
  39. * (opérateur RELA).
  40. *
  41. *******************************************************
  42. *
  43. TITRE 'DYNA5 - METHODES DE SOUS-STRUCTURATION' ;
  44. OPTIO DIME 3 ELEM SEG2 MODE TRID IMPI 0;
  45. *OPTION TRAC BENS;
  46. SOL=TABLE;
  47. TEMPS ;
  48. *
  49. *______________________________________________________
  50. * PREMIERE METHODE LIBRE-LIBRE
  51. *______________________________________________________
  52. *
  53. *---- Définition de la géométrie de la 1ère poutre ----
  54. *
  55. Q1 = 0. 0. 0.; Q2 =2. 0. 0.;
  56. LIG1= D 12 Q1 Q2;
  57. OEIL1 = 5 5 5 ;
  58. * TRACE OEIL1 LIG1 ;
  59. *
  60. *---- Type de calcul, matériau et caractéristiques ----
  61. *---------------- de la 1ère poutre -------------------
  62. *
  63. MO1= MODE LIG1 MECANIQUE ELASTIQUE POUT ;
  64.  
  65. MATPL1=MATE MO1 YOUNG 2.E11 NU 0.3 RHO 7800;
  66. CARPL1= 'CARA' MO1 SECT 0.1 INRY 2.08E-7 INRZ 3.33E-8
  67. TORS 2.41E-7;
  68. MATPL1=MATPL1 ET CARPL1;
  69. *
  70. *------ Rigidité, masse et conditions aux limites -----
  71. *----------------- de la 1ère poutre ------------------
  72. *
  73. RIGPL1 =RIGI MATPL1 MO1 ;
  74. MASPLA1=MASS MATPL1 MO1 ;
  75. BLOQ1=BLOQ DEPL ROTA Q1 ;
  76. RIGPLA1=RIGPL1 ET BLOQ1 ;
  77. *
  78. *---- Recherche des modes propres de la 1ère poutre ---
  79. *
  80. PR=PROG 0.4 ;
  81. MODPLA1=VIBR PROC PR RIGPLA1 MASPLA1 'SOLU';
  82. *
  83. *---- Définition de la géométrie de la 2ème poutre ----
  84. *
  85. Q3 =2. 0. 0.; Q4 =4. 0. 0.;
  86. LIG2= D 12 Q3 Q4;
  87. * TRAC OEIL1 LIG2 ;
  88. *
  89. POUTRE=LIG1 ET LIG2 ;
  90. SI(NEG GRAPH 'N');
  91. TRAC OEIL1 POUTRE;
  92. FINSI;
  93. *
  94. *---- Type de calcul, matériau et caractéristiques ----
  95. *---------------- de la 2ème poutre -------------------
  96. *
  97. MO2= MODE LIG2 MECANIQUE ELASTIQUE POUT ;
  98. MATPL2=MATE MO2 YOUNG 2.E11 NU 0.3 RHO 7800;
  99. CARPL2= 'CARA' MO2 SECT 0.1 INRY 2.08E-7 INRZ 3.33E-8
  100. TORS 2.41E-7;
  101. MATPL2=MATPL2 ET CARPL2;
  102. *
  103. *----- Rigidité, masse et conditions aux limites ------
  104. *---------------- de la 2ème poutre -------------------
  105. *
  106. RIGPL2 =RIGI MATPL2 MO2 ;
  107. MASPLA2=MASS MATPL2 MO2 ;
  108. BLOQ4=BLOQ DEPL ROTA Q4;
  109. RIGPLA2=RIGPL2 ET BLOQ4 ;
  110. *
  111. *---- Recherche des modes propres de la 2ème poutre ---
  112. *
  113. MODPLA2=VIBR PROC PR RIGPLA2 MASPLA2 'SOLU';
  114. *
  115. * Création des sous-structures associées aux 2 poutres
  116. *
  117. SPLA1=STRU RIGPLA1 MASPLA1 ;
  118. SPLA2=STRU RIGPLA2 MASPLA2 ;
  119. ELM2=ELST Q2 SPLA1;
  120. ELM3=ELST Q3 SPLA2;
  121. PR2=PROG 1. ;PR3=PROG -1. ;
  122. *
  123. *--- Description des liaisons entre sous-structures ---
  124. *
  125. LIUX=RELA ELM2 UX PR2 ELM3 UX PR3;
  126. LIUY=RELA ELM2 UY PR2 ELM3 UY PR3;
  127. LIUZ=RELA ELM2 UZ PR2 ELM3 UZ PR3;
  128. LIRX=RELA ELM2 RX PR2 ELM3 RX PR3;
  129. LIRY=RELA ELM2 RY PR2 ELM3 RY PR3;
  130. LIRZ=RELA ELM2 RZ PR2 ELM3 RZ PR3;
  131. ENCL=LIUY ET LIRX ET LIRY ET LIRZ ET LIUX ET LIUZ;
  132. *
  133. * -Création des bases modales associées aux 2 poutres -
  134. *
  135. SOL1 = SOLS ENCL SPLA1;
  136. SOL2 = SOLS ENCL SPLA2;
  137. BASE1=BASE SPLA1 MODPLA1 ENCL SOL1;
  138. BASE2=BASE SPLA2 MODPLA2 ENCL SOL2;
  139. *
  140. *------------------ Base modale totale ----------------
  141. *---------------- correspondant à la réunion des ------
  142. *------------------- deux bases élémentaires ----------
  143. *
  144. BASET= BASE1 ET BASE2 ;
  145. RIGMO = RIGI BASET ;
  146. MASMO = MASS BASET ;
  147. *
  148. *------------ Recherche de la 1ère fréquence ----------
  149. *------------ propre de la structure totale -----------
  150. *
  151. SOL.1= VIBR PROC (PROG 0.65) MASMO RIGMO 'SOLU';
  152. FREQ1=TIRE (SOL.1) FREQ RANG 1;
  153. *
  154. SOLL = SYNT (SOL.1) BASET;
  155. SOLB = TRADUIRE SOLL;
  156. DEFMOD1 = SOLB. 1 . DEFORMEE_MODALE;
  157. SI (EGA GRAPH 'O');
  158. TRAC (DEFO (LIG1 ET LIG2) DEFMOD1);
  159. FINSI;
  160. *
  161. *
  162. *______________________________________________________
  163. * DEUXIEME METHODE BLOQUE-BLOQUE
  164. *______________________________________________________
  165. *
  166. *
  167. *-------------------- Conditions aux limites ----------
  168. *--------------- sur l'extrémité de la 1ère poutre ----
  169. *
  170. BLOQUX2=BLOQ UX Q2;
  171. BLOQUY2=BLOQ UY Q2;
  172. BLOQUZ2=BLOQ UZ Q2;
  173. BLOQRX2=BLOQ RX Q2;
  174. BLOQRY2=BLOQ RY Q2;
  175. BLOQRZ2=BLOQ RZ Q2;
  176. *
  177. *------------ Rigidité de la 1ère poutre --------------
  178. *
  179. RIGPLA1=RIGPL1 ET BLOQ1
  180. ET BLOQUX2 ET BLOQUY2 ET BLOQUZ2 ET BLOQRX2 ET
  181. BLOQRY2 ET BLOQRZ2;
  182. *
  183. *---- Recherche des modes propres de la 1ère poutre ---
  184. *
  185. PR1=PROG 2.6 ;
  186. MODPLA1=VIBR PROC PR1 RIGPLA1 MASPLA1 'SOLU';
  187. *
  188. *---------------- Conditions aux limites --------------
  189. *------------ sur l'extrémité de la 2ème poutre -------
  190. *
  191. *- Le point Q2 de la 1ère poutre et le point Q3 de la -
  192. *-------- 2ème poutre ont les memes coordonnées -------
  193. *
  194. BLOQUX3=BLOQ UX Q3;
  195. BLOQUY3=BLOQ UY Q3;
  196. BLOQUZ3=BLOQ UZ Q3;
  197. BLOQRX3=BLOQ RX Q3;
  198. BLOQRY3=BLOQ RY Q3;
  199. BLOQRZ3=BLOQ RZ Q3;
  200. *
  201. *------------ Rigidité de la 2ème poutre --------------
  202. *
  203. RIGPLA2=RIGPL2 ET BLOQ4
  204. ET BLOQUX3 ET BLOQUY3 ET BLOQUZ3 ET BLOQRX3 ET
  205. BLOQRY3 ET BLOQRZ3;
  206. *
  207. *---- Recherche des modes propres de la 2ème poutre ---
  208. *
  209. MODPLA2=VIBR PROC PR1 RIGPLA2 MASPLA2 'SOLU';
  210. *
  211. * Création des sous structures associées aux 2 poutres
  212. *
  213. SPLA1=STRU RIGPLA1 MASPLA1 ;
  214. SPLA2=STRU RIGPLA2 MASPLA2 ;
  215. CLUX2=CLST SPLA1 BLOQUX2;
  216. CLUY2=CLST SPLA1 BLOQUY2;
  217. CLUZ2=CLST SPLA1 BLOQUZ2;
  218. CLRX2=CLST SPLA1 BLOQRX2;
  219. CLRY2=CLST SPLA1 BLOQRY2;
  220. CLRZ2=CLST SPLA1 BLOQRZ2;
  221. CLUX3=CLST SPLA2 BLOQUX3;
  222. CLUY3=CLST SPLA2 BLOQUY3;
  223. CLUZ3=CLST SPLA2 BLOQUZ3;
  224. CLRX3=CLST SPLA2 BLOQRX3;
  225. CLRY3=CLST SPLA2 BLOQRY3;
  226. CLRZ3=CLST SPLA2 BLOQRZ3;
  227. *
  228. *--- Description des liaisons entre sous-structures ---
  229. *
  230. PR2=PROG 1.;PR3=PROG 1.;
  231. LIUX=RELA CLUX2 LX PR2 CLUX3 LX PR3;
  232. LIUY=RELA CLUY2 LX PR2 CLUY3 LX PR3;
  233. LIUZ=RELA CLUZ2 LX PR2 CLUZ3 LX PR3;
  234. LIRX=RELA CLRX2 LX PR2 CLRX3 LX PR3;
  235. LIRY=RELA CLRY2 LX PR2 CLRY3 LX PR3;
  236. LIRZ=RELA CLRZ2 LX PR2 CLRZ3 LX PR3;
  237. ENCL=LIUY ET LIRX ET LIRY ET LIRZ ET LIUX ET LIUZ;
  238. *
  239. *--------- Création des bases modales élémentaires ----
  240. *------------- associées aux deux poutres -------------
  241. *
  242. SOL1 = SOLS ENCL SPLA1;
  243. SOL2 = SOLS ENCL SPLA2;
  244. BASE1=BASE SPLA1 MODPLA1 ENCL SOL1;
  245. BASE2=BASE SPLA2 MODPLA2 ENCL SOL2;
  246. *
  247. *------------------ Base modale totale ----------------
  248. *------------- correspondant à la réunion des ---------
  249. *---------------- deux bases élémentaires -------------
  250. *
  251. BASET= BASE1 ET BASE2 ;
  252. RIGMO = RIGI BASET ;
  253. MASMO = MASS BASET ;
  254. *
  255. *------------- Recherche de la 1ère fréquence ---------
  256. *-------------- propre de la structure totale ---------
  257. *
  258. SOL.2= VIBR PROC (PROG 0.65) MASMO RIGMO 'SOLU';
  259. FREQ2=TIRE SOL.2 FREQ RANG 1;
  260. *
  261. SOLL = SYNT (SOL.2) BASET;
  262. SOLB = TRADUIRE SOLL;
  263. DEFMOD2 = SOLB. 1 . DEFORMEE_MODALE;
  264. SI (EGA GRAPH 'O');
  265. TRAC (DEFO (LIG1 ET LIG2) DEFMOD2);
  266. FINSI;
  267. *
  268. *
  269. *______________________________________________________
  270. * TROISIEME METHODE LIBRE-BLOQUE
  271. *______________________________________________________
  272. *
  273. *
  274. *------------ Rigidité et conditions aux limites ------
  275. *----------------- de la 1ère poutre ------------------
  276. *
  277. RIGPLA1=RIGPL1 ET BLOQ1
  278. ET BLOQUX2 ET BLOQUY2 ET BLOQUZ2 ET BLOQRX2 ET BLOQRY2
  279. ET BLOQRZ2;
  280. *
  281. *---- Recherche des modes propres de la 1ère poutre ---
  282. *
  283. PR1=PROG 2.6;
  284. MODPLA1=VIBR PROC PR1 RIGPLA1 MASPLA1 'SOLU';
  285. *
  286. *-------- Rigidité et conditions aux limites ----------
  287. *--------------- de la 2ème poutre --------------------
  288. *
  289. RIGPLA2=RIGPL2 ET BLOQ4;
  290. *
  291. *---- Recherche des modes propres de la 2ème poutre ---
  292. *
  293. PR2=PROG 0.41 ;
  294. MODPLA2=VIBR PROC PR2 RIGPLA2 MASPLA2 'SOLU';
  295. *
  296. *------------ Définition des sous-structures ----------
  297. *
  298. SPLA2=STRU RIGPLA2 MASPLA2 ;
  299. CL3=ELST SPLA2 Q3;
  300. *
  301. *---- Définition des liaisons entre sous-structures ---
  302. *
  303. PR=PROG 1;
  304. REL1=RELA CL3 UX 1 CLUX2 LX PR ;
  305. REL2=RELA CL3 UY 1 CLUY2 LX PR ;
  306. REL3=RELA CL3 UZ 1 CLUZ2 LX PR ;
  307. REL4=RELA CL3 RX 1 CLRX2 LX PR ;
  308. REL5=RELA CL3 RY 1 CLRY2 LX PR ;
  309. REL6=RELA CL3 RZ 1 CLRZ2 LX PR ;
  310. ENCL=REL1 ET REL2 ET REL3 ET REL4 ET REL5 ET REL6;
  311. *
  312. *--------- Création des bases modales élèmentaires ----
  313. *
  314. *SOL1 = SOLS ENCL SPLA1;
  315. *SOL2 = SOLS ENCL SPLA2;
  316. BASE1=BASE SPLA1 MODPLA1 ENCL;
  317. BASE2=BASE SPLA2 MODPLA2 ENCL;
  318. *
  319. *------------------ Base modale totale ----------------
  320. *------------- correspondant à la réunion des ---------
  321. *----------------- deux bases élémentaires ------------
  322. *
  323. BASET= BASE1 ET BASE2 ;
  324. RIGMO = RIGI BASET ;
  325. MASMO = MASS BASET ;
  326. *
  327. *---------- Recherche de la 1ère fréquence ------------
  328. *---------- propre de la structure totale -------------
  329. *
  330. SOL.3= VIBR PROC (PROG 0.65) MASMO RIGMO 'SOLU';
  331. FREQ3=TIRE SOL.3 FREQ RANG 1;
  332. *______________________________________________________
  333. * RESULTATS
  334. *______________________________________________________
  335. SAUT PAGE ; SAUT 4 LIGNE ;
  336. MESS 'ON DOIT TROUVER UNE FREQUENCE DE 0.65 HZ';
  337. SAUT 2 LIGNE ;
  338. MESS 'PAR LA PREMIERE METHODE ON TROUVE : ' FREQ1;
  339. MESS 'PAR LA DEUXIEME METHODE ON TROUVE : ' FREQ2;
  340. MESS 'PAR LA TROISIEME METHODE ON TROUVE : ' FREQ3;
  341. TEMPS ;
  342. * CODE FONCTIONNEMENT
  343. SOLL = SYNT (sol.3) BASET;
  344. SOLB = TRADUIRE SOLL;
  345. DEFMOD3 = SOLB. 1 . DEFORMEE_MODALE;
  346. SI (EGA GRAPH 'O');
  347. TRAC (DEFO (LIG1 ET LIG2) DEFMOD3);
  348. FINSI;
  349. *
  350. *
  351. REF=0.65;
  352. RESI1=ABS((FREQ1-REF)/REF);
  353. RESI2=ABS((FREQ2-REF)/REF);
  354. RESI3=ABS((FREQ3-REF)/REF);
  355. SI((RESI1 <EG 5E-2) ET (RESI2 <EG 5E-2) ET
  356. (RESI3 <EG 5E-2));
  357. ERRE 0;
  358. SINO;
  359. ERRE 5;
  360. FINSI;
  361. FIN;
  362. *
  363. *
  364.  
  365.  
  366.  
  367.  
  368.  
  369.  
  370.  
  371.  
  372.  
  373.  
  374.  
  375.  
  376.  
  377.  
  378.  
  379.  
  380.  
  381.  
  382.  
  383.  
  384.  
  385.  
  386.  

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