1 : C DYCHPL SOURCE CHAT 05/01/12 23:06:57 5004
2 : SUBROUTINE DYCHPL(XDEP,XVIT,XDPLAS,XRAID,XJEU,XAMO,XSEUIL,XFL,
3 : & IPERM,iannul)
4 : IMPLICIT INTEGER(I-N)
5 : IMPLICIT REAL*8(A-H,O-Z)
6 : *--------------------------------------------------------------------*
7 : * *
8 : * Op{rateur DYNE : algorithme de Fu - de Vogelaere *
9 : * ________________________________________________ *
10 : * *
11 : * Calcul de la force de choc pour un choc {l{mentaire point-plan
12 : * elastoplastique, la force restant toujours positive
13 : * *
14 : * Param}tres: *
15 : * *
16 : * e XDEP Valeur du d{placement. *
17 : * e XVIT Valeur de la vitesse normale
18 : * e XRAID Valeur de la raideur. *
19 : * e XJEU Valeur de la limite de force elastique (positive)
20 : * e XAMO Valeur de l'amortissement
21 : * e XSEUIL Valeur de la force plastique
22 : * e IPERM =1 si liaison permanente *
23 : * es XDPLAS Valeur du deplacement plastique (du pas precedent
24 : * en entree, du pas actuel en sortie)
25 : * s XFL Valeur de la force de choc. *
26 : * *
27 : * *
28 : * Auteur, date de cr{ation: *
29 : * *
30 : * I. POLITOPOULOS 19/09/1994
31 : * *
32 : *--------------------------------------------------------------------*
33 :
34 : XFL = 0.D0
35 : XVAL = XDEP - XJEU - XDPLAS
36 :
37 : IF (IPERM.EQ.1) THEN
38 : IF (iannul.EQ.0) THEN
39 : XE = XSEUIL/XRAID
40 : XFLR = XRAID * XVAL
41 : IF (XFLR.LT.0) XCRIT=-XSEUIL
42 : IF (XFLR.GE.0) XCRIT= XSEUIL
43 :
44 : DXFLP = XFLR - XCRIT
45 : IF (ABS(XFLR).LE.ABS(XCRIT)) THEN
46 : XFL = XFLR
47 : ELSE
48 : XFL = XCRIT
49 : XDPLAS = XDEP - XJEU - XCRIT/XRAID
50 : ENDIF
51 : XFL = - XFL - XAMO*XVIT
52 :
53 : ENDIF
54 : ELSE
55 : IF ((XVAL .GE. 0.D0).AND.(iannul.EQ.0)) THEN
56 : XE = XSEUIL/XRAID
57 : XFLR = - XRAID * XVAL
58 : DXFLP = XSEUIL - (ABS (XFLR))
59 : IF (DXFLP.LE.0.D0) THEN
60 : XDPLAS = XDEP - XJEU - XE
61 : XFLR = - XSEUIL
62 : ENDIF
63 : XFL = XFLR - XAMO*XVIT
64 : IF (XFL.GT.0.D0) XFL = 0.D0
65 : ELSE
66 : XFL = 0.D0
67 : ENDIF
68 : ENDIF
69 :
70 : END
71 :
72 :
73 :
74 :
75 :
© Cast3M 2003 - All rights reserved.
Disclaimer