Télécharger dync02.dgibi

Retour à la liste

Numérotation des lignes :

  1. ************************************************************************
  2. *
  3. * Calcul d'un rotor de type Jeffcott avec contact frottant
  4. * avec la methode HBM (DYNC)
  5. *
  6. * ___ ^ Z
  7. * _|_ kc,mu |__
  8. * _______ jeu |_ \_
  9. * | k,c | | | \ \
  10. * |----------| m |----------| +--|--|-----> Y
  11. * | |_______| R R+jeu
  12. * ___ jeu
  13. * _|_ kc,mu
  14. *
  15. * ^ Z
  16. * |
  17. * +----> X
  18. *
  19. *
  20. * Ref : [Xie et al., MSSP, 2015]
  21. * Auteur : PBZ CEA-ENSTA, 2020-07-07
  22. *
  23. ************************************************************************
  24.  
  25.  
  26. *********************************************************************
  27. * OPTIONS GENERALES
  28. *********************************************************************
  29.  
  30. OPTI DIME 3 ELEM SEG2 EPSI LINEAIRE;
  31. OPTION 'TRAC' 'PSC' 'EPTR' 5 'POTR' 'HELVETICA_16';
  32. OPTION 'EPSI' 'LINEAIRE';
  33.  
  34.  
  35. *********************************************************************
  36. * PARAMETRES
  37. *********************************************************************
  38.  
  39. * nombre de modes a calculer
  40. NMODE = 2;
  41. * nombre de harmoniques a calculer et NFFT
  42. NHBM = 7;
  43. NFFT = 256;
  44.  
  45. * mu de frottement
  46. * mu = 0.05;
  47. mu = 0.2;
  48.  
  49.  
  50. *********************************************************************
  51. * DEFINITION DU MODELE
  52. *********************************************************************
  53.  
  54. * rotor data
  55. m = 1. ;
  56. k = 100.;
  57. c = 5. ;
  58. jeu = 0.105;
  59. R = 20.*jeu ;
  60. ebal = 0.1;
  61. * Fbal = ebal * m * (Omega**2) ; avec Omega = vitesse de rotation en rad/s
  62.  
  63. * contact parameters
  64. kchoc= 25. * k;
  65. mu_s = mu;
  66. mu_d = mu;
  67. si (mu_s < mu_d); erreur 21; finsi;
  68. Kt = 10.*kchoc;
  69. Ct = 0.5*(m*Kt)**0.5;
  70.  
  71. * deduced parameter
  72. w = (k/m)**0.5;
  73. wHz = w / (2.*pi);
  74. wchoc = (kchoc/m)**0.5;
  75. xi = 0.5 * c / ((kchoc*m)**0.5);
  76. mess 'w=' w '(' wHz ')Hz < wchoc=' wchoc;
  77. mess 'xi=' xi;
  78.  
  79. * critical timestep (Diff Centrees, no damping)
  80. wadhe = ((k+Kt)/m)**0.5;
  81. dtc = 2. / wadhe;
  82.  
  83. * !!! on adimensionne apres !!!
  84.  
  85.  
  86. ************************************************************************
  87. * BASE MODALE
  88. ************************************************************************
  89.  
  90. * point physique
  91. p1 = 0. 0. 0.;
  92.  
  93. * base modale = 2 modes
  94. palfa1 = 0. 0. 0.;
  95. phi1 = MANU 'CHPO' p1 3 'UX' 0. 'UY' 1. 'UZ' 0. 'NATURE' 'DIFFUS';
  96. palfa2 = 0. 0. 0.;
  97. phi2 = MANU 'CHPO' p1 3 'UX' 0. 'UY' 0. 'UZ' 1. 'NATURE' 'DIFFUS';
  98.  
  99. TBAS = TABL 'BASE_MODALE';
  100. TBAS . 'MODES' = TABL 'BASE_DE_MODES';
  101. TBAS . 'MODES' . 1 = TABL 'MODES';
  102. TBAS . 'MODES' . 1 . 'POINT_REPERE' = palfa1;
  103. TBAS . 'MODES' . 1 . 'FREQUENCE' = wHz/wchoc;
  104. TBAS . 'MODES' . 1 . 'MASSE_GENERALISEE' = 1.;
  105. TBAS . 'MODES' . 1 . 'DEFORMEE_MODALE' = phi1;
  106. TBAS . 'MODES' . 2 = TABL 'MODES';
  107. TBAS . 'MODES' . 2 . 'POINT_REPERE' = palfa2;
  108. TBAS . 'MODES' . 2 . 'FREQUENCE' = wHz/wchoc;
  109. TBAS . 'MODES' . 2 . 'MASSE_GENERALISEE' = 1.;
  110. TBAS . 'MODES' . 2 . 'DEFORMEE_MODALE' = phi2;
  111.  
  112.  
  113. ************************************************************************
  114. * AMORTISSEMENT
  115. ************************************************************************
  116.  
  117. TAMOR = TABL 'AMORTISSEMENT';
  118. Ctot = MANU 'RIGIDITE' 'TYPE' 'AMORTISSEMENT'
  119. (palfa1 et palfa2) (MOTS 'ALFA') (prog (c / (m*wchoc)));
  120. TAMOR . 'AMORTISSEMENT' = Ctot ;
  121.  
  122.  
  123. ************************************************************************
  124. * LIAISONS
  125. ************************************************************************
  126.  
  127. * contact-frottant point_cercle_frottement
  128. TL1 = TABL 'LIAISON_ELEMENTAIRE';
  129. TL1 . 'TYPE_LIAISON' = mot 'POINT_CERCLE_FROTTEMENT';
  130. TL1 . 'SUPPORT' = p1;
  131. TL1 . 'NORMALE' = (1. 0. 0.) ;
  132. TL1 . 'EXCENTRATION' = (0. 0. 0.) ;
  133. TL1 . 'RAIDEUR' = 1.;
  134. TL1 . 'RAYON' = 1. ;
  135. TL1 . 'COEFFICIENT_GLISSEMENT' = mu_d;
  136. TL1 . 'COEFFICIENT_ADHERENCE' = mu_s;
  137. TL1 . 'RAIDEUR_TANGENTIELLE' = -1.*Kt;
  138. * TL1 . 'AMORTISSEMENT_TANGENTIEL'= Ct;
  139. TL1 . 'AMORTISSEMENT_TANGENTIEL'= 1.E-4;
  140. TL1 . 'VITESSE_ENTRAINEMENT' = 0. ;
  141. * TL1 . 'VITESSE_ENTRAINEMENT' = -1.*R*Omega;
  142.  
  143. * stockage
  144. TLB = TABL 'LIAISON_B';
  145. TLB . 1 = TL1;
  146. TLIA = TABL 'LIAISON';
  147. TLIA . 'LIAISON_B' = TLB;
  148.  
  149.  
  150. ************************************************************************
  151. * CHARGEMENT
  152. ************************************************************************
  153.  
  154. Fbal = ebal*m/(m*jeu) ;
  155. Fchpo_Y = MANU 'CHPO' 1 palfa1 'FALF' Fbal;
  156. Fchpo_Z = MANU 'CHPO' 1 palfa2 'FALF' Fbal;
  157. *+++ BALOURD == 1 pour ajouter w**2 a Fbal
  158.  
  159.  
  160. *********************************************************************
  161. * CALCUL NON LINEAIRE
  162. *********************************************************************
  163.  
  164. * chargement spatial : GFC1
  165. TAB_CHAR = TABLE 'CHARGEMENT';
  166. TAB_CHAR . 1 = Fchpo_Y;
  167. TAB_CHAR . -1 = Fchpo_Z;
  168. *Pour ajouter w**2 dans la definition du chargement
  169. TAB_CHAR . 'BALOURD' = 1;
  170.  
  171. * Approximation pour la solution initiale à fréquence fixe
  172. VEC_INIT = TABLE 'INITIAL';
  173. VEC_INIT . 'FREQUENCE' = 0.5/wchoc ;
  174. * VEC_INIT . 3 = 0.1;
  175.  
  176. * Paramètres numériques pour la continuation
  177. PAR_NUM = TABLE 'PARAMETRES_NUMERIQUES';
  178. PAR_NUM . 'TYPE' = MOT 'FORC';
  179. PAR_NUM . 'VAL_FIN' = 60./wchoc;
  180. PAR_NUM . 'DS0' = 0.01;
  181. * PAR_NUM . 'DSMAX' = 0.1;
  182. PAR_NUM . 'DSMAX' = 0.025;
  183. PAR_NUM . 'DSMIN' = 1.E-4;
  184. PAR_NUM . 'ITERMOY' = 3.4;
  185. PAR_NUM . 'ITERMAX' = 20;
  186. PAR_NUM . 'ANGLE_MIN' = 0.;
  187. PAR_NUM . 'ANGLE_MAX' = 40.;
  188. PAR_NUM . 'ISENS' = 1.;
  189. PAR_NUM . 'NBPAS' = 2000;
  190. * PAR_NUM . 'NBPAS' = 200;
  191. PAR_NUM . 'TOLERANCE' = 1.E-7;
  192. * PAR_NUM . 'CALC_JAC' = VRAI;
  193. PAR_NUM . 'CALC_JAC' = FAUX;
  194.  
  195. *----------------------------------------------------------------------*
  196. * Appel a l'operateur
  197. * opti impi 2;
  198. opti impi 1;
  199. TCON = DYNC TBAS TAB_CHAR TLIA TAMOR VEC_INIT PAR_NUM
  200. NHBM NFFT;
  201. opti impi 0;
  202. *----------------------------------------------------------------------*
  203.  
  204.  
  205. ************************************************************************
  206. * POST-TRAITEMENT
  207. ************************************************************************
  208.  
  209.  
  210. *----------- RECUP ------------
  211. * frequence
  212. FREQ = TCON . 'REPONSE' . 'FREQUENCE';
  213. * MAX de la NORME2 des COEFFICIENTS de FOURIER pour chaque freq
  214. AMPS = TCON . 'REPONSE' . 'NORME_DEPLACEMENT';
  215. * stabilite
  216. STAB = TCON . 'REPONSE' . 'STABILITE';
  217.  
  218. * On cree les evolutions correspondants
  219. XY = EVOL MANU '\w' FREQ '|Q|' AMPS;
  220.  
  221. TABSTB = TABLE;
  222. TABSTB . 'LIGNE_VARIABLE' = TABL ENTI;
  223. TABSTB . 'LIGNE_VARIABLE' . 1 = STAB;
  224.  
  225. DESS XY TABSTB 'TITR' 'Jeffcott : courbe de reponse';
  226. DESS XY TABSTB XBOR 0.14 0.16 'TITR' 'Jeffcott : courbe de reponse';
  227.  
  228.  
  229. * On sauvegarde les exposants de FLOQUET et on leur assigne une couleur
  230. TEXP_R = TCON . 'REPONSE' . 'EXPOSANT_REEL';
  231. TEXP_I = TCON . 'REPONSE' . 'EXPOSANT_IMAGINAIRE';
  232. exp_R = VIDE 'EVOLUTIO';
  233. exp_I = VIDE 'EVOLUTIO';
  234. coco = MOTS 'ROUG' 'BLEU' 'VERT' 'AZUR' 'ORAN' 'VIOL' 'TURQ' ;
  235. repe bexp;
  236. si (non (exis TEXP_R &bexp)); quit bexp ; finsi;
  237. coco_i = EXTR coco &bexp;
  238. exp_R = exp_R ET (evol coco_i 'MANU' '\w' FREQ '\m_{R}' TEXP_R . &bexp);
  239. exp_I = exp_I ET (evol coco_i 'MANU' '\w' FREQ '\m_{I}' TEXP_I . &bexp);
  240. fin bexp;
  241.  
  242. DESS exp_R 'TITR' 'Jeffcott : exposants de Floquet';
  243. DESS exp_I 'TITR' 'Jeffcott : exposants de Floquet';
  244.  
  245. * evolution de l'harmonique 1 des modes en y et en z
  246. j = 1;
  247. y1 = TCON . REPONSE . COEFFICIENTS . j . palfa1;
  248. z1 = TCON . REPONSE . COEFFICIENTS . j . palfa2;
  249. ev1 = (EVOL 'AZUR' 'MANU' 'LEGE' 'mode 1 (y)' '\w' FREQ 'Q^{j=1}' y1)
  250. et (EVOL 'ORAN' 'MANU' 'LEGE' 'mode 2 (Z)' '\w' FREQ 'Q^{j=1}' Z1);
  251. DESS ev1 'TITR' 'Jeffcott : Reponse de l harmonique 1 des 2 modes' LEGE SO;
  252.  
  253.  
  254. * trace des bifurcation
  255. NBIF = DIME TCON . 'BIFURCATION' . 'TYPE';
  256. prLP_Q = PROG; prLP_w = PROG;
  257. prBP_Q = PROG; prBP_w = PROG;
  258. prPD_Q = PROG; prPD_w = PROG;
  259. prNS_Q = PROG; prNS_w = PROG;
  260. REPE BBIF NBIF;
  261. typbif = EXTR TCON . 'BIFURCATION' . 'TYPE' &BBIF;
  262. Qbif = EXTR TCON . 'BIFURCATION' . 'NORME_DEPLACEMENT' &BBIF;
  263. wbif = EXTR TCON . 'BIFURCATION' . 'FREQUENCE' &BBIF;
  264. SI (EGA typbif 'LP'); prLP_Q = prLP_Q et Qbif; prLP_w = prLP_w et wbif; FINSI;
  265. SI (EGA typbif 'BP'); prBP_Q = prBP_Q et Qbif; prBP_w = prBP_w et wbif; FINSI;
  266. SI (EGA typbif 'PD'); prPD_Q = prPD_Q et Qbif; prPD_w = prPD_w et wbif; FINSI;
  267. SI (EGA typbif 'NS'); prNS_Q = prNS_Q et Qbif; prNS_w = prNS_w et wbif; FINSI;
  268. FIN BBIF;
  269. evLP = EVOL MANU LEGE 'LP' '\w' prLP_w '|Q|' prLP_Q;
  270. evBP = EVOL MANU LEGE 'BP' '\w' prBP_w '|Q|' prBP_Q;
  271. evPD = EVOL MANU LEGE 'PD' '\w' prPD_w '|Q|' prPD_Q;
  272. evNS = EVOL MANU LEGE 'NS' '\w' prNS_w '|Q|' prNS_Q;
  273. * evBIF = evLP et evBP et evPD et evNS;
  274. * TABSTB . 2 = MOT 'NOLI MARQ TRIR';
  275. * TABSTB . 3 = MOT 'NOLI MARQ TRID';
  276. * TABSTB . 4 = MOT 'NOLI MARQ CARR';
  277. * TABSTB . 5 = MOT 'NOLI MARQ LOSA';
  278. evBIF = evLP ;
  279. TABSTB . 2 = MOT 'NOLI MARQ TRIU';
  280.  
  281.  
  282. DESS (XY ET evBIF )TABSTB 'TITR' 'Jeffcott : courbe de reponse + bifurcations';
  283.  
  284.  
  285. ************************************************************************
  286. * TEST DE BON FONCTIONNEMENT
  287. ************************************************************************
  288. * amplitude max et w corespondant
  289. * nombre, type de localisation des bifurcations
  290.  
  291.  
  292.  
  293. ************************************************************************
  294. * PERFORMANCES
  295. ************************************************************************
  296. TEMP IMPR MAXI CPU;
  297.  
  298. * opti donn 5 trac X;
  299. FIN ;
  300.  
  301.  
  302.  
  303.  

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