Télécharger la plaquette Cast3M


fr en

Télécharger dycham2.eso

Retour à la liste

Numérotation des lignes : 

dycham2

  1. C DYCHAM2   SOURCE    BP208322  19/02/25    21:15:47     10120          
  2. C DYCHAM2    SOURCE    LAVARENN  96/08/05    21:17:56     2237
  3.       SUBROUTINE DYCHAM2(XDEP,XVIT,XJEU,XAMO,XABSCI,XORDON,NIP,
  4.      &            NLIAB,I,XFL,IPERM,iannul)
  5.       IMPLICIT INTEGER(I-N)
  6.       IMPLICIT REAL*8(A-H,O-Z)
  7. *--------------------------------------------------------------------*
  8. *                                                                    *
  9. *     Operateur DYNE                                                 *
  10. *     ________________________________________________               *
  11. *                                                                    *
  12. *     Calcul de la force de choc pour un choc elementaire.           *
  13. *                Cas d'une LOI_DE_COMPORTEMENT elastique             *
  14. *     Parametres:                                                    *
  15. *                                                                    *
  16. * e   XDEP    Valeur du d{placement.                                 *
  17. * e   XVIT    Valeur de la vitesse.                                  *
  18. * e   XJEU    Valeur du jeu.                                         *
  19. * e   XAMO    Valeur de l'amortissement.                             *
  20. * e   IPERM   =1 si liaison permanente                               *
  21. * e   XABSCI  Abscisses loi elastique                                *
  22. * e   XORDON  Ordonnees loi elastique                                *
  23. *                                                                    *
  24. *  s  XFL     Valeur de la force de choc.                            *
  25. *                                                                    *
  26. *                                                                    *
  27. *     Auteur, date de creation:                                      *
  28. *                                                                    *
  29. *     Bruno FONTAINE, le 26/02/98                                    *
  30. *                                                                    *
  31. *--------------------------------------------------------------------*
  32.       REAL*8 XABSCI(NLIAB,*),XORDON(NLIAB,*)
  33.       XFL = 0.D0
  34. *
  35. cbp      if (iannul.eq.0) then
  36.  
  37.       IF (IPERM.EQ.2) THEN
  38. cbp        XVAL = abs(XDEP)
  39.         XVAL = abs(XDEP-XJEU)
  40.         CALL LIRANG(XVAL,XABSCI,XORDON,I,NLIAB,NRG,XPENTE,NIP)
  41.         XFLA = (XORDON(I,NRG) + XPENTE*(XVAL-XABSCI(I,NRG)))
  42. cbp        IF (XDEP.LT.0D0) XFLA = -XFLA
  43.         IF (XDEP.LT.XJEU) XFLA = -XFLA
  44.         XFL = -XFLA - XAMO * XVIT
  45.       ELSE
  46. cbp         IF (XDEP.GE.XJEU) THEN
  47.         IF ((XDEP.GE.XJEU).AND.(iannul.EQ.0)) THEN
  48.             XVAL = XDEP - XJEU
  49.             CALL LIRANG(XVAL,XABSCI,XORDON,I,NLIAB,NRG,XPENTE,NIP)
  50.             XFL = - (XORDON(I,NRG) + XPENTE*(XVAL-XABSCI(I,NRG)))
  51.      &            - XAMO * XVIT
  52.             IF (XFL.GT.0.D0) XFL = 0.D0
  53.         ENDIF
  54.       ENDIF
  55.  
  56. cbp      endif
  57. *
  58. *
  59.       END
  60.  
  61.  
  62.  
  63.  
  64.  
  65.  
  66.  

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