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 :
© Cast3M 2003 - All rights reserved.
Disclaimer