Télécharger for_cont.procedur

Retour à la liste

Numérotation des lignes :

  1. * FOR_CONT PROCEDUR FANDEUR 14/10/10 21:15:14 8178
  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 'MODE' '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.  
  41.  
  42.  

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