Télécharger la plaquette Cast3M


fr en

Télécharger dychpl2.eso

Retour à la liste

Numérotation des lignes : 

dychpl2

  1. C DYCHPL2   SOURCE    BP208322  19/02/25    21:15:49     10120          
  2. C DYCHPL2    SOURCE    LAVARENN  96/08/05    21:19:09     2237
  3.       SUBROUTINE DYCHPL2(XDEP,XVIT,XDPLAS,XRAID,XJEU,XAMO,XSEUIL
  4.      &           ,XABSCI,XORDON,NIP,NLIAB,I,XFL,IPERM,iannul)
  5.       IMPLICIT INTEGER(I-N)
  6.       IMPLICIT REAL*8(A-H,O-Z)
  7. *--------------------------------------------------------------------*
  8. *                                                                    *
  9. *     Operateur DYNE : algorithme de Fu - de Vogelaere               *
  10. *     ________________________________________________               *
  11. *                                                                    *
  12. *     Calcul de la force de choc pour un choc elementaire point-plan
  13. *     elastoplastique, la force restant toujours positive
  14. *                                                                    *
  15. *     Parametres:                                                    *
  16. *                                                                    *
  17. * e   XDEP    Valeur du d{placement.                                 *
  18. * e   XVIT    Valeur de la vitesse normale
  19. * e   XJEU    Valeur de la limite de force elastique (positive)
  20. * e   XAMO    Valeur de l'amortissement
  21. * e   XABSCI  Abscisses loi elastique                                *
  22. * e   XORDON  Ordonnees loi elastique                                *
  23. * e   XSEUIL  Valeur de la force plastique
  24. * e   IPERM   =2 si liaison permanente                               *
  25. * es  XDPLAS  Valeur du deplacement plastique (du pas precedent
  26. *              en entree, du pas actuel en sortie)
  27. *  s  XFL     Valeur de la force de choc.                            *
  28. *                                                                    *
  29. *                                                                    *
  30. *     Auteur, date de cr{ation:                                      *
  31. *                                                                    *
  32. *     Bruno FONTAINE, le 26/02/98                                    *
  33. *                                                                    *
  34. *--------------------------------------------------------------------*
  35.       REAL*8 XABSCI(NLIAB,*),XORDON(NLIAB,*)
  36.  
  37.       XFL = 0.D0
  38.       XVAL = XDEP - XJEU - XDPLAS
  39.       IF (IPERM.EQ.2) THEN
  40. cbp         IF (iannul.EQ.0) THEN
  41.         CALL LIRANG(XSEUIL,XORDON,XABSCI,I,NLIAB,NRG,XPENTE,NIP)
  42.           XE = XABSCI(I,NRG) + XPENTE*(XSEUIL-XORDON(I,NRG))
  43.         CALL LIRANG(XVAL,XABSCI,XORDON,I,NLIAB,NRG,XPENTE,NIP)
  44.         XFLR = XORDON(I,NRG) + XPENTE*(XVAL-XABSCI(I,NRG))
  45.           IF (XFLR.LT.0) XCRIT=-XSEUIL
  46.           IF (XFLR.GE.0) XCRIT= XSEUIL
  47.  
  48.           DXFLP =  XFLR - XCRIT
  49.             IF (ABS(XFLR).LE.ABS(XCRIT)) THEN
  50.               XFL = XFLR
  51.             ELSE
  52.               XFL = XCRIT
  53.               XDPLAS = XDEP - XJEU - XCRIT/XRAID
  54.             ENDIF
  55.           XFL = - XFL - XAMO*XVIT
  56. cbp         ENDIF
  57.       ELSE
  58.         IF ((XVAL .GE. 0.D0).AND.(iannul.EQ.0)) THEN
  59.         CALL LIRANG(XSEUIL,XORDON,XABSCI,I,NLIAB,NRG,XPENTE,NIP)
  60.           XE = XABSCI(I,NRG) + XPENTE*(XSEUIL-XORDON(I,NRG))
  61.         CALL LIRANG(XVAL,XABSCI,XORDON,I,NLIAB,NRG,XPENTE,NIP)
  62.         XFLR = - (XORDON(I,NRG) + XPENTE*(XVAL-XABSCI(I,NRG)))
  63.           DXFLP =  XSEUIL - (ABS (XFLR))
  64.             IF (DXFLP.LE.0.D0) THEN
  65.               XDPLAS = XDEP - XJEU - XE
  66.               XFLR = - XSEUIL
  67.             ENDIF
  68.           XFL = XFLR - XAMO*XVIT
  69.           IF (XFL.GT.0.D0) XFL = 0.D0
  70.         ELSE
  71.           XFL = 0.D0
  72.         ENDIF
  73.       ENDIF
  74.  
  75.       END
  76.  
  77.  
  78.  
  79.  
  80.  
  81.  
  82.  
  83.  
  84.  

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