Télécharger symplaq.dgibi

Retour à la liste

Numérotation des lignes :

  1. * fichier : symplaq.dgibi
  2. ************************************************************************
  3. ************************************************************************
  4. OPTI DIME 3 ELEM TRI3 ;
  5. *OPTI IMPR 'symplaq.out';
  6. P0 = 0. 0. 0. ;
  7. P1 = 1. 0. 0. ;
  8. P2 = 0. 0. 1.;
  9. P3 = 0. 1. 0. ;
  10. L1 = CERC 5 P1 P0 P3 ;
  11. OEIL = 1000. 1000. 1000.;
  12. L2 = L1 TOUR 90. P0 P2 ;
  13. L3 = L2 TOUR 90. P0 P2 ;
  14. L4 = L3 TOUR 90. P0 P2 ;
  15. P4 = -1. 0. 0.;
  16. L01 = L1 ET L2;
  17. L02 = P4 DROIT 11 P1;
  18. L0 = L01 ET L02;
  19. ELIM 0.001 L0;
  20. ELIM 0.001 L01;
  21. ELIM 0.001 L02;
  22. S1 = SURF L0 'PLANE';
  23. ELIM 0.001 L0 S1;
  24. ELIM 0.001 L01 S1;
  25. ELIM 0.001 L02 S1;
  26. S2 = S1 SYME 'DROIT' P4 P1;
  27. S3 = S1 ET S2;
  28. S3 = ORIEN S3 'POINT' P2;
  29. S3 = VERSENS S3;
  30. LB = L1 ET L2 ET L3 ET L4;
  31. ELIM 0.001 LB S3;
  32. *OPTI SAUV 'symplaq.bin';
  33. *SAUV S3;
  34. *OPTI DONN 5;
  35. MOD1 = 'MODELE' S1 'MAGNETODYNAMIQUE'
  36. 'POTENTIEL_VECTEUR' 'ISOTROPE' 'ROT3';
  37. MAT1 = 'MATE' MOD1 'ETA' 1.E-4 'PERM' 1. 'EPAI' 1.;
  38. M11 = 'MUTU' MOD1 MAT1 S1;
  39. * CAS SYMETRIE : CHANGEMENT D'ORIENTATION *
  40. M21 = 'MUTU' MOD1 MAT1 S2;
  41. M21 = -1.*M21;
  42. M1 = M11 ET M21;
  43. R1 = 'RESI' MOD1 MAT1 ;
  44. *******************************
  45. * CONSTITUTION DU SYSTEME *****
  46. *******************************
  47. * DONNEE DU CHAMP INUCTEUR : 100T/s ****
  48. B1Z = MANU CHPO S1 1 SCAL 0.;
  49. B2Z = MANU CHPO S1 1 SCAL 100.;
  50. DTB=1.;
  51. DBZDT = (B2Z-B1Z)/DTB;
  52. X1 = (COOR 1 S1) ;
  53. Y1 = (COOR 2 S1) ;
  54. X2 = X1**2;
  55. Y2 = Y1**2;
  56. RAY = (X2+Y2)**0.5 ;
  57. DAPHIDT = 0.5*RAY*DBZDT;
  58. * ON PASSE EN CARTESIEN *
  59. ANGL1 = 'ATG' Y1 X1 ;
  60. PIS4 = 'ATG' 1. 1. ;
  61. PIS2 = 2.*PIS4;
  62. ANGL2 = ANGL1 + PIS2 ;
  63. COSA1 = 'COS' ANGL2 ;
  64. SINA1 = 'SIN' ANGL2 ;
  65. DAXDT = DAPHIDT*COSA1;
  66. DAXDT = EXCO DAXDT SCAL AX;
  67. DAXDT = CHAN DAXDT ATTRIBUT NATURE DISCRET;
  68. DAYDT = DAPHIDT*SINA1;
  69. DAYDT = EXCO DAYDT SCAL AY;
  70. DAYDT = CHAN DAYDT ATTRIBUT NATURE DISCRET;
  71. DAZDT = MANU CHPO S1 1 AZ 0.;
  72. DAZDT = CHAN DAZDT ATTRIBUT NATURE DISCRET;
  73. DADT = DAXDT ET DAYDT ET DAZDT ;
  74. VADT = VECT DADT AX AY AZ 0.005 ROUGE;
  75. DADTN=CNEQ MOD1 DADT;
  76. * DONNEE DES OPERATEURS ********
  77. TETA=0.5;
  78. DT = 1.E-5;
  79. TEMPS=0.;
  80. OPE1 = M1 ET (R1*(TETA*DT)) ;
  81. OPE2 = BLOQ 'FC' L01;
  82. OPE3 = OPE1 ET OPE2;
  83. SOU1 = M1 ET (R1*((TETA-1.)*DT)) ;
  84. SOU2 = SOU1 ;
  85. *SOU2 = SOU1 ET OPE2 ;
  86. SOL1 = MANU CHPO S1 1 FC 0. ;
  87. NDT = 5;
  88. NSORT = NDT ;
  89. REPE BOUC NDT;
  90. IDT = &BOUC ;
  91. TEMPS = TEMPS+DT;
  92. SMB11 = SOU2*SOL1;
  93. SMB12 = DT*DADTN;
  94. SMB1 = SMB11 - SMB12;
  95. SOL2 = RESO OPE3 SMB1;
  96. SOL0 = OPE3*SOL2 ;
  97. SOL1 = SOL2 ;
  98. FIN BOUC;
  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