1 : * FOR_CONT PROCEDUR BAZE 97/01/30 21:15:04 2468 2 : 'DEBPROC' FOR_CONT GEO*'MAILLAGE' SOL1*'CHPOINT' COURI*'FLOTTANT' ; 3 : ************************************************************************ 4 : * MAGNETOSTATIQUE 2D POTENTIEL VECTEUR 5 : * GEO maillage 6 : * SOL1 SOLTUTION EN POTENTIEL DEFINIE AU MOINS SUR BLOC 7 : * calcul des forces par integrale de contour 8 : * COURI densite de courant sur geo ( constante ) 9 : * LOG1 logique valant vrai si pb axisym ( defaut pb plan 10 : * SORTIE 11 : * RESULTANTE DE LA FORCE AU CDG DE GEO 12 : * MOMENT PAR RAPPORT A L ORIGINE 13 : * EN PLAN FORCES PAR UNITE DE LONGUEUR 14 : * EN AXISYMETRIQUE FORCES / RADIAN (attention signification ) 15 : ************************************************************************ 16 : CCONT = 'CONTOUR' ('ORIENTER' GEO ) ; 17 : OBSEG = CCONT 'MODL' 'THERMIQUE' 'ISOTROPE' ; 18 : AA = 'CHAN' 'CHAM' SOL1 OBSEG 'GRAVITE'; 19 : NBSEG = CCONT 'NBEL' ; 20 : IEL = 0 ;SOMX = 0. ; SOMY = 0. ;SMM = 0. ; 21 : * 22 : 'REPETER' BOUC NBSEG ; 23 : IEL = IEL + 1 ; SEGC= CCONT 'ELEM' IEL ; 24 : I1 = SEGC 'POINT' 'INITIAL' ; I2 = SEGC 'POINT' 'FINAL' ; 25 : X1 = 'COOR' 1 I1 ; X2 = 'COOR' 1 I2 ; RX = ( X1 + X2) / 2. ; 26 : Y1 = 'COOR' 2 I1 ; Y2 = 'COOR' 2 I2 ; RY = ( Y1 + Y2) / 2. ; 27 : DX = X2 - X1 ; DY = Y2 - Y1 ; 28 : AME = 'EXTR' AA 'T' 1 IEL 1 ; 29 : SOMX = SOMX + ( AME * DX) ; 30 : SOMY = SOMY + ( AME * DY) ; 31 : MOMM = (RX * DY) + (RY * DX ) ; 32 : SMM= SMM + ( AME * MOMM ) ; 33 : 'FIN' BOUC ; 34 : * 35 : FXX = COURI * SOMY ; FYY = -1. * COURI * SOMX ; 36 : MOMT = -1. * COURI * SMM ; 37 : CDG1 = 'BARY' CCONT ; 38 : RFORC = 'MANU' 'CHPO' CDG1 2 'FX' FXX 'FY' FYY ; 39 : 'FINPROC' RFORC SMM ; 40 :
© Cast3M 2003 - All rights reserved.
Disclaimer