Télécharger la plaquette Cast3M


fr en

Télécharger CORF1.dgibi

Retour à la liste

Numérotation des lignes : 

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

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