Télécharger rotplaq.dgibi

Retour à la liste

Numérotation des lignes :

  1. * fichier : rotplaq.dgibi
  2. ************************************************************************
  3. ************************************************************************
  4. OPTI DIME 3 ELEM TRI3 ;
  5. P0 = 0. 0. 0. ;
  6. P1 = 1. 0. 0. ;
  7. P2 = 0. 0. 1.;
  8. P3 = 0. 1. 0. ;
  9. L1 = CERC 5 P1 P0 P3 ;
  10. OEIL = 1000. 1000. 1000.;
  11. L2 = L1 TOUR 90. P0 P2 ;
  12. L3 = L2 TOUR 90. P0 P2 ;
  13. L4 = L3 TOUR 90. P0 P2 ;
  14. P4 = -1. 0. 0.;
  15. L01 = L1 ET L2;
  16. L02 = P4 DROIT 11 P1;
  17. L0 = L01 ET L02;
  18. ELIM 0.001 L0;
  19. ELIM 0.001 L01;
  20. ELIM 0.001 L02;
  21. S1 = SURF L0 'PLANE';
  22. ELIM 0.001 L0 S1;
  23. ELIM 0.001 L01 S1;
  24. ELIM 0.001 L02 S1;
  25. S2 = S1 TOUR 180. P0 P2;
  26. S3 = S1 ET S2;
  27. S3 = ORIEN S3 'POINT' P2;
  28. S3 = VERSENS S3;
  29. LB = L1 ET L2 ET L3 ET L4;
  30. ELIM 0.001 LB S3;
  31. *OPTI SAUV 'rotplaq.bin';
  32. *SAUV S3;
  33. *OPTI DONN 5;
  34. MOD1 = 'MODELE' S1 'MAGNETODYNAMIQUE'
  35. 'POTENTIEL_VECTEUR' 'ISOTROPE' 'ROT3';
  36. MAT1 = 'MATE' MOD1 'ETA' 1.E-4 'PERM' 1. 'EPAI' 1.;
  37. M11 = 'MUTU' MOD1 MAT1 S1;
  38. * CAS ROTATION *
  39. M21 = 'MUTU' MOD1 MAT1 S2;
  40. M1 = M11 ET M21;
  41. R1 = 'RESI' MOD1 MAT1 ;
  42. *******************************
  43. * CONSTITUTION DU SYSTEME *****
  44. *******************************
  45. * DONNEE DU CHAMP INUCTEUR : 100T/s ****
  46. B1Z = MANU CHPO S1 1 SCAL 0.;
  47. B2Z = MANU CHPO S1 1 SCAL 100.;
  48. DTB=1.;
  49. DBZDT = (B2Z-B1Z)/DTB;
  50. X1 = (COOR 1 S1) ;
  51. Y1 = (COOR 2 S1) ;
  52. X2 = X1**2;
  53. Y2 = Y1**2;
  54. RAY = (X2+Y2)**0.5 ;
  55. DAPHIDT = 0.5*RAY*DBZDT;
  56. * ON PASSE EN CARTESIEN *
  57. ANGL1 = 'ATG' Y1 X1 ;
  58. PIS4 = 'ATG' 1. 1. ;
  59. PIS2 = 2.*PIS4;
  60. ANGL2 = ANGL1 + PIS2 ;
  61. COSA1 = 'COS' ANGL2 ;
  62. SINA1 = 'SIN' ANGL2 ;
  63. DAXDT = DAPHIDT*COSA1;
  64. DAXDT = EXCO DAXDT SCAL AX;
  65. DAXDT = CHAN DAXDT ATTRIBUT NATURE DISCRET;
  66. DAYDT = DAPHIDT*SINA1;
  67. DAYDT = EXCO DAYDT SCAL AY;
  68. DAYDT = CHAN DAYDT ATTRIBUT NATURE DISCRET;
  69. DAZDT = MANU CHPO S1 1 AZ 0.;
  70. DAZDT = CHAN DAZDT ATTRIBUT NATURE DISCRET;
  71. DADT = DAXDT ET DAYDT ET DAZDT ;
  72. VADT = VECT DADT AX AY AZ 0.005 ROUGE;
  73. DADTN=CNEQ MOD1 DADT;
  74. * DONNEE DES OPERATEURS ********
  75. TETA=0.5;
  76. DT = 1.E-5;
  77. TEMPS=0.;
  78. OPE1 = M1 ET (R1*(TETA*DT)) ;
  79. OPE2 = BLOQ 'FC' L01;
  80. OPE3 = OPE1 ET OPE2;
  81. SOU1 = M1 ET (R1*((TETA-1.)*DT)) ;
  82. SOU2 = SOU1 ;
  83. *SOU2 = SOU1 ET OPE2 ;
  84. SOL1 = MANU CHPO S1 1 FC 0. ;
  85. NDT = 5;
  86. NSORT = NDT ;
  87. REPE BOUC NDT;
  88. IDT = &BOUC ;
  89. TEMPS = TEMPS+DT;
  90. SMB11 = SOU2*SOL1;
  91. SMB12 = DT*DADTN;
  92. SMB1 = SMB11 - SMB12;
  93. SOL2 = RESO OPE3 SMB1;
  94. SOL0 = OPE3*SOL2 ;
  95. SOL1 = SOL2 ;
  96. OUT = (IDT - (NSORT*(IDT/NSORT)));
  97. FIN BOUC;
  98. MESS ' ******** INSTANT FINAL *******';
  99. CURR = DECO MOD1 MAT1 SOL1;
  100. CURX = EXCO 'FC,X' CURR ;
  101. CURX = CHAN 'CHPO' CURX MOD1;
  102. CRX2 = CURX**2;
  103. CURY = EXCO 'FC,Y' CURR ;
  104. CURY = CHAN 'CHPO' CURY MOD1;
  105. CRY2 = CURY**2;
  106. CURZ = EXCO 'FC,Z' CURR ;
  107. CURZ = CHAN 'CHPO' CURZ MOD1;
  108. CRZ2 = CURZ**2;
  109. CURP = CURX ET CURY ET CURZ ;
  110. CURN = (CRX2+CRY2+CRZ2)**0.5;
  111. JMAX = MAXI CURN;
  112. LIST JMAX;
  113. JREF = 17290.;
  114. CURP = CURP/JMAX;
  115. CURV = VECT CURP 'FC,X' 'FC,Y' 'FC,Z' 0.2 ROUGE;
  116. ERRJ = ABS(JREF-JMAX)/JREF;
  117. LIST ERRJ;
  118. SI (ERRJ > 0.02);
  119. ERREUR 5;
  120. FINSI;
  121. *OPTI DONN 5;
  122. FIN;
  123.  
  124.  
  125.  
  126.  
  127.  

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