Télécharger dyfor1.eso

Retour à la liste

Numérotation des lignes :

  1. C DYFOR1 SOURCE CHAT 05/01/12 23:09:00 5004
  2. C DYFORC SOURCE LAVARENN 96/11/05 21:22:41 2357
  3. SUBROUTINE DYFOR1(XPTB,XVALB,IPALB,IPLIB,XPALB,NLIAB,NPLB,IND
  4. &,IND2,PDTS2,I,iannul,IESC,IROLE,XCHPFB)
  5. IMPLICIT INTEGER(I-N)
  6. IMPLICIT REAL*8(A-H,O-Z)
  7. *--------------------------------------------------------------------*
  8. * *
  9. * Opérateur DYNE : algorithme de Fu - de Vogelaere *
  10. * ________________________________________________ *
  11. * *
  12. * Calcul des forces de choc sur base B pour les liaisons de *
  13. * type ligne_cercle *
  14. * *
  15. * Paramètres: *
  16. * *
  17. * e XPTB Tableau des déplacements des points *
  18. * e IPALB Renseigne sur la liaison. *
  19. * e IPLIB Tableau contenant les numéros "DYNE" de la liaison. *
  20. * e XPALB Tableau contenant les paramètres de la liaison. *
  21. * es XVALB Tableau contenant les variables internes de liaisons. *
  22. * es XCHPFB Tableau contenant les valeurs des futurs chpoints *
  23. * e NLIAB Nombre de liaisons sur la base B. *
  24. * e NPLB Nombre total de points intervenant dans les liaisons. *
  25. * e IND Indice du pas. *
  26. * e I numéro de la liaison. *
  27. * e IESC numéro DYNE du noeud esclave *
  28. * e IROLE indique qui est le maitre et l'esclave *
  29. * = 0 premiere passe , =1 deuxième passe *
  30. * *
  31. * *
  32. * Auteur, date de création: **
  33. * *
  34. * Auteur, date de création: *
  35. * *
  36. * Ibrahim PINTO, 05/97 ,liaisons ligne_cercle *
  37. *--------------------------------------------------------------------*
  38. *
  39. INTEGER IPALB(NLIAB,*),IPLIB(NLIAB,*)
  40. INTEGER ICAND(2)
  41. REAL*8 XPALB(NLIAB,*),XPTB(NPLB,4,*),XVALB(NLIAB,4,*)
  42. REAL*8 XXXN(3),XCHPFB(2,NLIAB,4,*),XTE(3),XNET
  43.  
  44.  
  45.  
  46. *
  47. XEPSI = 0.0001
  48. IGP = IPALB(I,2)
  49. IDIM = IPALB(I,3)
  50. * ILOCAL = IPALB(I,23)
  51. ITYP = IPALB(I,1)
  52.  
  53. XPOID = 0.D0
  54.  
  55.  
  56. XRAY=XPALB(I,2)
  57. *
  58. IF (ITYP.EQ.37 .OR. ITYP.EQ.39) THEN
  59. ID1 = 6
  60. ELSE
  61. ID1 = 7
  62. ENDIF
  63. *
  64. IF (IROLE.EQ.0) THEN
  65. KMAI = 0
  66. IMAI = ID1 +4*IDIM
  67. IBUT = IMAI + IPALB(I,21)*IDIM
  68. KBUT = IPALB(I,21)
  69. LMAI = ID1 + (IPALB(I,21)+IPALB(I,22)+4)*IDIM
  70. LBUT = LMAI +IPALB(I,21)*IDIM
  71. JMAI = 26 + IPALB(I,21)+IPALB(I,22)
  72. JBUT = JMAI + IPALB(I,21)
  73. MBUT = LBUT +IPALB(I,22)*IDIM +
  74. & IPALB(I,21)
  75. NBUT = MBUT +IPALB(I,21)+IPALB(I,22)
  76. ELSE
  77. KMAI = IPALB(I,21)
  78. IBUT = ID1 + 4*IDIM
  79. IMAI = IBUT + IPALB(I,21)*IDIM
  80. KBUT = 0
  81. LBUT = ID1 + (IPALB(I,21)+IPALB(I,22)+4)*IDIM
  82. LMAI = LBUT +IPALB(I,21)*IDIM
  83. JBUT = 26 + 2*(IPALB(I,21)+IPALB(I,22))
  84. JMAI = JBUT +IPALB(I,21)
  85. MBUT = LMAI +IPALB(I,22)*IDIM
  86. NBUT = MBUT +IPALB(I,21)+IPALB(I,22)
  87. ENDIF
  88. ID2 = ID1 + IDIM
  89. ID4 = ID1 +3*IDIM
  90.  
  91. CALL DYCAN1(IPALB,IPLIB,XPALB,XPTB,IND,I,NLIAB,NPLB,
  92. & XXXN,XNET,XTE,XPOID,ICAND,IESC,IROLE)
  93.  
  94.  
  95. XVITN =0.D0
  96. IDCAN1 = IMAI + (ICAND(1)-1)*IDIM
  97. IDCAN2 = IMAI + (ICAND(2)-1)*IDIM
  98.  
  99.  
  100. DO 80 ID=1,IDIM
  101. *
  102. * Déplacement du point de contact au pas courant
  103.  
  104.  
  105. x1=(XPTB(IPLIB(I,KMAI+ICAND(2)),IND,ID))*(1-XPOID)
  106.  
  107. x2=(XPTB(IPLIB(I,KMAI+ICAND(1)),IND,ID))*XPOID
  108.  
  109. xde2=x1+x2
  110.  
  111.  
  112. XVALB(I,IND,3+ID) = XVALB(I,IND,3+ID)+
  113. &XPALB(I,IDCAN1+ID)*XPOID +XPALB(I,IDCAN2+ID)*(1-XPOID)
  114.  
  115.  
  116. * Déplacement du point de contact au pas précédent
  117. XDM2 =(XPTB(IPLIB(I,KMAI+ICAND(2)),IND2,ID))*(1-XPOID)
  118. & + (XPTB(IPLIB(I,KMAI+ICAND(1)),IND2,ID))*XPOID
  119.  
  120.  
  121. *-----Déplacement de la butee(point de contact sur le cercle)
  122. *-----on a neglige les rotations des butees
  123. XDB = (XPTB(IPLIB(I,KBUT+IESC),IND,ID)
  124. & -XPTB(IPLIB(I,KBUT+IESC),IND2,ID))
  125.  
  126.  
  127.  
  128. * Vitesse de glissement
  129. XPALB(I,ID2+ID)=(XDB+XDM2-XDE2)/PDTS2
  130. * Vitesse normale
  131. XVITN = XVITN + ((XDB+XDM2-XDE2)*XXXN(ID)/PDTS2)
  132. 80 CONTINUE
  133.  
  134.  
  135. DO 82 ID=1,IDIM
  136. * Vitesse tangentielle
  137. XPALB(I,ID2+ID) = XPALB(I,ID2+ID) - (XVITN*XXXN(ID))
  138. 82 CONTINUE
  139. * Calcul de la force de choc
  140. XPALB(I,1) = XPALB(I,MBUT+IESC)
  141.  
  142. IF (ITYP.EQ.37 .OR. ITYP.EQ.39) THEN
  143. CALL DYCHE5(XNET,XTE,XXXN,IDIM,IGP,XPALB,NLIAB,I
  144. & ,IPALB,XFN,XFT,XPUS,iannul)
  145. ELSE
  146.  
  147. XPALB(I,7) = XPALB(I,NBUT+IESC)
  148. CALL DYCHA5(XNET,XTE,XXXN,XVITN,IDIM,IGP,XPALB,NLIAB,
  149. &I,IPALB,XFN,XFT,XPUS,iannul)
  150.  
  151.  
  152. ENDIF
  153.  
  154.  
  155. XVALB(I,IND,3) = XVALB(I,IND,3) + XVITN
  156. IPALB(I,2) = IGP
  157.  
  158.  
  159. IF (IGP.EQ.1) THEN
  160. PS =0.D0
  161. DO 84 ID=1,IDIM
  162. PS=PS+(XPALB(I,ID2+ID)*XPALB(I,ID2+ID))
  163. 84 CONTINUE
  164. XVITT = SQRT(PS)
  165. ELSE
  166. XVITT = 0.D0
  167. ENDIF
  168. XVALB(I,IND,11) = XVALB(I,IND,11)+XVITT
  169.  
  170.  
  171.  
  172.  
  173.  
  174. IF (XNET.GE.XRAY) THEN
  175.  
  176. DO 90 ID=1,IDIM
  177. XPALB(I,ID+LMAI+(ICAND(1)-1)*IDIM) =XPALB(I,ID+LMAI+
  178. &(ICAND(1)-1)*IDIM)-XPOID * (XFN*XXXN(ID) + XPALB(I,ID4+ID))
  179. XPALB(I,ID+LMAI+(ICAND(2)-1)*IDIM) =XPALB(I,ID+LMAI+(ICAND
  180. &(2)-1)*IDIM)- ( 1.-XPOID )*(XFN*XXXN(ID) + XPALB(I,ID4+ID))
  181. XPALB(I,ID+LBUT+(IESC-1)*IDIM) =XPALB(I,ID+LBUT+
  182. &(IESC-1)*IDIM) + XFN * XXXN(ID) + XPALB(I,ID4+ID)
  183. 90 CONTINUE
  184. XCHPFB(1,I,IND,IPLIB(I,KBUT+IESC))=
  185. &XCHPFB(1,I,IND,IPLIB(I,KBUT+IESC))+XFN
  186.  
  187.  
  188.  
  189.  
  190. XCHPFB(2,I,IND,IPLIB(I,KBUT+IESC))=
  191. &XCHPFB(2,I,IND,IPLIB(I,KBUT+IESC))+ABS(XFT)
  192.  
  193. XCHPFB(1,I,IND,IPLIB(I,KMAI+ICAND(1)))=
  194. &XCHPFB(1,I,IND,IPLIB(I,KMAI+ICAND(1)))+XPOID *XFN
  195. XCHPFB(2,I,IND,IPLIB(I,KMAI+ICAND(1)))=
  196. &XCHPFB(2,I,IND,IPLIB(I,KMAI+ICAND(1)))+XPOID *ABS(XFT)
  197. XCHPFB(1,I,IND,IPLIB(I,KMAI+ICAND(2)))=
  198. &XCHPFB(1,I,IND,IPLIB(I,KMAI+ICAND(2)))+( 1.-XPOID )*XFN
  199. XCHPFB(2,I,IND,IPLIB(I,KMAI+ICAND(2)))=
  200. &XCHPFB(2,I,IND,IPLIB(I,KMAI+ICAND(2)))+( 1.-XPOID )*ABS(XFT)
  201. IF (XPOID.GT.XEPSI) THEN
  202. IPALB(I,JMAI+ICAND(1))=1
  203. ENDIF
  204. IF ((1.0-XPOID).GT.XEPSI) THEN
  205. IPALB(I,JMAI+ICAND(2))=1
  206. ENDIF
  207. IPALB(I,JBUT+IESC) = 1
  208. XVALB(I,IND,1) = XVALB(I,IND,1)+XFN
  209. XVALB(I,IND,10) = XVALB(I,IND,10)+ ABS(XFT)
  210. XVALB(I,IND,12) = XVALB(I,IND,12)+ XPUS
  211. ENDIF
  212. *
  213.  
  214.  
  215. END
  216.  
  217.  
  218.  
  219.  
  220.  

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