Télécharger ddfour.procedur

Retour à la liste

Numérotation des lignes :

  1. * DDFOUR PROCEDUR TTMF3 96/06/04 21:15:12 2157
  2. 'DEBPROC' DDFOUR GEO*'MAILLAGE' CHARM*'MAILLAGE' NHARM*'ENTIER'
  3. SOL*'CHPOINT' RHARM*'FLOTTANT' ORIG/'POINT' LIS*'LOGIQUE'
  4. ITYPE/'ENTIER' ;
  5. *********************************************************************
  6. * ANALYSE HARMONIQUE DU POTENTIEL VECTEUR
  7. * GEO MAILLAGE SUPPORT SOLUTION GENERALE *
  8. * CHARM ARC DE CERCLE SUR LEQUEL ON A LE POTENTIEL *
  9. * NHARM NOMBRE D HARMONIQUES A CALCULER *
  10. * SOL SOLUTION EN POTENTIEL *
  11. * RHARM RAYON DE NORMALISATION *
  12. * ORIG ORIGINE CERCLE D ANALYSE *
  13. * LIS LOGIQUE FAUX SI PAS DE LISSAGE POLYNOMIAL *
  14. * ITYPE ORDRE DE SYMETRIE POUR L EXPANSION POLY SI LIS=VRAI*
  15. *********************************************************************
  16. 1P = CHARM 'POINT' 'INITIAL' ;2P= CHARM 'POINT' 'FINAL' ;
  17. 'SI' ( 'EXISTE' ORIG ) ;
  18. RCIRC = 'NORM' ( 1P 'MOINS' ORIG ) ;
  19. 'SINON' ;
  20. RCIRC = 'NORM' 1P ;ORIG = 0. 0. ;
  21. 'FINS' ;
  22. X1 Y1 = 'COOR' ( 1P 'MOINS' ORIG );X2 Y2 = 'COOR' 2P ;
  23. PT = 2P 'PROJETER' (( Y1 * -1. ) X1 ) 'DROITE' ORIG 1P ;
  24. H = 'NORM' ( 2P 'MOINS' PT ) ;
  25. 'SI' ( (X1 + X2 ) < 1.E-3 ) ;
  26. NUM = 2 ;
  27. 'MESS' 'SOLUTION DONNEE SUR 180 DEGRES ' ;
  28. 'SINON' ;
  29. ANG = 'ATG' H ( 'NORM' ( PT 'MOINS' ORIG)) ;
  30. NUM = ( 'ENTI' ( 360.1 / ANG )) ;
  31. 'MESS' 'SOLUTION DONNEE SUR ' ANG ' DEGRES ' ;
  32. 'FINSI' ;
  33. 'SI' ( 'NON' ( LIS )) ;
  34. SCIRC = INT_COMP GEO SOL CHARM ;
  35. SCIRC = SCIRC 'NOMC' 'T' ;
  36. 'SINON' ;
  37. * lissage polynomial
  38. * CHLI = PROI POLY CHARM GEO SOL ITYPE 'PLAN' ;
  39. CHLI = 'PROI' POLY GEO CHARM SOL ITYPE 'PLAN' ;
  40. SCIRC = ('EXCO' CHLI 'A' ) 'NOMC' 'T' ;
  41. 'FINS' ;
  42. RRAP = RHARM / RCIRC ;
  43. CQTT RTOT = F_S2PI CHARM SCIRC NUM ORIG ;
  44. EQTT = 'EVOL' 'ROUG' 'CHPO' CQTT 'T' RTOT ;
  45. AAAA= 'EXTR' EQTT 'ABSC' ;
  46. * QTOT= 'EXTR' EQTT 'ORDO' ;
  47. * LIST QTOT ;
  48. * NI1= (DIME AAAA ) - 1 ;
  49. * IVAL = VALE ELEM ;
  50. * OPTION ELEM SEG2 ;I1 = 0. 0. ; I2 = 360. 0. ;
  51. * RTOT = DROI NI1 I1 I2 ;
  52. OBA_RMA = 'MODEL' RTOT 'MECANIQUE' 'COQ2';
  53. * OBA_RMA= 'AFFECTE' RTOT ('MODELE' 'STANDARD') COQ2;
  54. PERIOD= 'MAXI' AAAA ;
  55. MULT= 360. / PERIOD ;
  56. ANGVEC= AAAA * MULT ;
  57. * MESS ' PERIODE ' PERIOD ;
  58. N=0;
  59. 'MESS' ' ANALYSE CONFORME A POISSON HARMON ';
  60. 'MESS' ' HARMONIQUE POTENTIEL CHAMP ';
  61. 'MESS' ' NOMBRE DE SEGMENTS DE DR ' NI1 ;
  62. 'MESS' ' DIME DE ANGVEC ' ( DIME ANGVEC ) ;
  63. CHPA= 'MANU' 'CHPO' RTOT 1 'SCAL' ANGVEC ;
  64. * CQTT = MANU CHPO RTOT 1 'POT' QTOT ;
  65. 'REPETER' BLOCIT NHARM ;
  66. N= N + 1 ;
  67. RAN= RRAP ** N ;
  68. ANGLEN= CHPA * N ;
  69. CNX= 'COS' ANGLEN ;
  70. SNX= 'SIN' ANGLEN ;
  71. EVREL= CQTT * CNX ;EVIMA= CQTT * SNX ;
  72. CEVREL =EVREL 'CHAN' 'CHAM' OBA_RMA 'RIGIDITE' ;
  73. CEVIMA = EVIMA 'CHAN' 'CHAM' OBA_RMA 'RIGIDITE' ;
  74. AK= (( 'INTG' OBA_RMA CEVREL ) / PERIOD ) * 2. * RAN ;
  75. BK= (( 'INTG' OBA_RMA CEVIMA ) / PERIOD ) * 2. * RAN ;
  76. CK= (( AK * AK ) + ( BK * BK ) ) ** .5 ;
  77. FF= ( 'FLOT' N) / RHARM ;
  78. AKK= AK * FF ;
  79. BKK= BK * FF ;
  80. CKK= CK * FF ;
  81. 'SI' (N 'EGA' 1 );AKN= 'PROG' AKK ; 'SINON' ;
  82. AKN = AKN 'ET' ('PROG' AKK ) ; 'FINSI' ;
  83. 'SAUTER' 1 'LIGNE' ;
  84. * MESS N AK BK CK AKK BKK CKK ;
  85. 'MESS' N AK BK CK ;
  86. 'MESS' N AKK BKK CKK ;
  87. 'FIN' BLOCIT ;
  88. * OPTION ELEM IVAL ;
  89. 'FINPROC' AKN ANGVEC ;
  90.  
  91.  

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