Télécharger d2vlfb.eso

Retour à la liste

Numérotation des lignes :

d2vlfb
  1. C D2VLFB SOURCE BP208322 20/09/18 21:15:18 10718
  2. C DEVLFB SOURCE LAVARENN 96/10/30 21:18:15 2349
  3.  
  4. SUBROUTINE D2VLFB(Q1,Q2,FTOTA,NA1,IPALB,IPLIB,XPALB,XVALB,NLIAB,
  5. & XPHILB,JPLIB,NPLB,IDIMB,FTOTB,FTOTBA,XPTB,PDT,T,
  6. & NPAS,IBASB,IPLSB,INMSB,IORSB,NSB,NPLSB,NA2,IND,
  7. & FEXPSM,NPC1,IERRD,FTEST2,XABSCI,XORDON,NIP,
  8. & FEXB,RIGIDE,IAROTA,XCHPFB,
  9. & KTOTXA,KTOTVA,KTOTXB,KTOTVB,GETJAC)
  10.  
  11. *--------------------------------------------------------------------*
  12. * Operateur DYN* : Calcul des forces de choc base B *
  13. * et projection vers la base A *
  14. *--------------------------------------------------------------------*
  15. * *
  16. * Parametres: *
  17. * *
  18. * e Q1(.,.) Vecteur des deplacements generalises. *
  19. * e Q2(.,.) Vecteur des vitesses generalisees. *
  20. * es FTOTA Forces exterieures totalisees sur la base A. *
  21. * es KTOTXA Matrice tangente,X des efforts non-lineaires (base A) *
  22. * es KTOTVA Matrice tangente,V des efforts non-lineaires (base A) *
  23. * es KTOTXB Matrice tangente,X des efforts non-lineaires (base B) *
  24. * es KTOTVB Matrice tangente,V des efforts non-lineaires (base B) *
  25. * e NA1 Nombre total d'inconnues en base A. *
  26. * e IPALB Renseigne sur le type de la liaison. *
  27. * e IPLIB Tableau contenant les numeros "DYNE" de la liaison. *
  28. * e XPALB Tableau contenant les parametres de la liaison. *
  29. * es XVALB Tableau contenant les variables internes des liaisons *
  30. * e NLIAB Nombre de liaisons sur la base B. *
  31. * e XPHILB Tableau des vecteurs propres aux points de liaisons. *
  32. * e JPLIB Tableau contenant les numeros "GIBI" des liaisons. *
  33. * e NPLB Nombre total de points intervenant dans les liaisons. *
  34. * e IDIMB Nombre de directions. *
  35. * e IND Indice du pas. *
  36. * e XABSCI Tableau contenant les abscisses de la loi plastique *
  37. * pour la liaison point-point- ... -plastique *
  38. * e XORDON Tableau contenant les ordonnees de la loi plastique *
  39. * pour la liaison point-point- ... -plastique *
  40. * e NIP Nb de pts dans l'evolution de la loi de comportement *
  41. * e RIGIDE Vrai si corps rigide, faux sinon *
  42. * e USEQ2 .TRUE. si on doit utiliser Q2 *
  43. * e GETJAC .TRUE. si on doit calculer la jacobienne KTOTXA,KTOTVA *
  44. * *
  45. *--------------------------------------------------------------------*
  46.  
  47. IMPLICIT INTEGER(I-N)
  48. IMPLICIT REAL*8(A-H,O-Z)
  49.  
  50. INTEGER IPALB(NLIAB,*),IPLIB(NLIAB,*),JPLIB(*)
  51. INTEGER IBASB(*),IPLSB(*),INMSB(*),IORSB(*),IAROTA(*)
  52. REAL*8 XPALB(NLIAB,*),Q1(NA1,*),Q2(NA1,*),FTOTA(NA1,*)
  53. REAL*8 XVALB(NLIAB,4,*),XPHILB(NSB,NPLSB,NA2,*),XPTB(NPLB,2,*)
  54. REAL*8 FTOTB(NPLB,*),FTOTBA(*),FEXPSM(NPLB,NPC1,2,*)
  55. REAL*8 XABSCI(NLIAB,*),XORDON(NLIAB,*),FEXB(NPLB,2,*)
  56. REAL*8 XCHPFB(2,NLIAB,4,NPLB)
  57. REAL*8 FTEST2(NPLB,6)
  58. LOGICAL RIGIDE
  59. REAL*8 KTOTXA(NA1,NA1), KTOTVA(NA1,NA1)
  60. REAL*8 KTOTXB(NPLB,IDIMB,IDIMB),KTOTVB(NPLB,IDIMB,IDIMB)
  61. c REAL*8 KTOTXBA(NA1,NA1), KTOTVBA(NA1,NA1)
  62. LOGICAL GETJAC
  63.  
  64. *--------------------------------------------------------------------*
  65. * Initialisations a 0
  66. *--------------------------------------------------------------------*
  67. *
  68. DO ID = 1,IDIMB
  69. DO IP = 1,NPLB
  70. FTOTB(IP,ID) = 0.D0
  71. ENDDO
  72. ENDDO
  73.  
  74. IF (GETJAC) THEN
  75. DO ij=1,IDIMB
  76. DO ii=1,IDIMB
  77. DO ip=1,NPLB
  78. KTOTXB(ip,ii,ij)=0.D0
  79. KTOTVB(ip,ii,ij)=0.D0
  80. ENDDO
  81. ENDDO
  82. ENDDO
  83. ENDIF
  84.  
  85. *--------------------------------------------------------------------*
  86. * Recombinaison des deplacements aux points de choc
  87. *--------------------------------------------------------------------*
  88. *
  89. IF (IDIMB.EQ.6) THEN
  90. IDIM=3
  91. ELSE
  92. IDIM=2
  93. ENDIF
  94.  
  95. *old CALL DEVRCO(Q1,NA1,XPTB,NPLB,XPHILB,NSB,NPLSB,NA2,IDIMB,IBASB,
  96. CALL DEVRCO(Q1,Q2,NA1,XPTB,NPLB,XPHILB,NSB,NPLSB,NA2,IDIMB,IBASB,
  97. & IPLSB,INMSB,IORSB,IND,IAROTA)
  98. * --> XPTB(:,1)=x_n XPTB(:,2,:)=\dot{q}_n
  99.  
  100. *--------------------------------------------------------------------*
  101. * Calcul des forces de choc sur base B
  102. *--------------------------------------------------------------------*
  103.  
  104. CALL D2VLB1(FTOTB,XPTB,IPALB,IPLIB,JPLIB,XPALB,XVALB,NLIAB,
  105. & NPLB,IDIMB,PDT,NPAS,IND,FEXPSM,NPC1,IERRD,
  106. & FTEST2,XABSCI,XORDON,NIP,XCHPFB,
  107. & KTOTXB,KTOTVB,GETJAC)
  108.  
  109. IF (IERRD.NE.0) RETURN
  110.  
  111. *--------------------------------------------------------------------*
  112. * Calcul des moments dans le cas des modes de rotation rigide
  113. *--------------------------------------------------------------------*
  114.  
  115. IF (RIGIDE) THEN
  116. CALL DEVMOM(FTOTB,Q1,FEXB,XPHILB,IAROTA,IBASB,IPLSB,INMSB,
  117. & NA2,NA1,NSB,NPLSB,NPLB,IND,IDIM)
  118. ENDIF
  119.  
  120. *--------------------------------------------------------------------*
  121. * Projection des forces base B sur base A
  122. * + eventuellement des jacobiennes
  123. *--------------------------------------------------------------------*
  124.  
  125. IF (GETJAC) THEN
  126. c CALL DEVPR2(XPHILB,FTOTB,FTOTBA,KTOTXB,KTOTVB,KTOTXBA,KTOTVBA,
  127. c & IBASB,INMSB,IPLSB,IORSB,NSB,NPLSB,NA2,IDIMB,NPLB,NA1)
  128. c c * Ajout des matrices projetees sur base A
  129. c DO J=1,NA1
  130. c DO I=1,NA1
  131. c KTOTXA(I,J) = KTOTXA(I,J) + KTOTXBA(I,J)
  132. c KTOTVA(I,J) = KTOTVA(I,J) + KTOTVBA(I,J)
  133. c ENDDO
  134. c ENDDO
  135. cbp : commente car fait directement dans DEVPRO2 -> permet de supprimer KTOTXBA et KTOTVBA
  136. cbp : rem : on pourrait aussi le faire pour FTOTA et FTOTBA
  137. CALL DEVPR2(XPHILB,FTOTB,FTOTBA,KTOTXB,KTOTVB,KTOTXA,KTOTVA,
  138. & IBASB,INMSB,IPLSB,IORSB,NSB,NPLSB,NA2,IDIMB,NPLB,NA1)
  139. ELSE
  140. CALL DEVPRO(XPHILB,FTOTB,FTOTBA,IBASB,INMSB,IPLSB,IORSB,NSB,
  141. & NPLSB,NA2,IDIMB,NPLB,NA1)
  142. ENDIF
  143.  
  144. *--------------------------------------------------------------------*
  145. * Ajout des forces projetees aux forces exterieures sur base A
  146. *--------------------------------------------------------------------*
  147.  
  148. DO I = 1,NA1
  149. FTOTA(I,IND) = FTOTA(I,IND) + FTOTBA(I)
  150. ENDDO
  151.  
  152. END
  153.  
  154.  
  155.  
  156.  
  157.  
  158.  
  159.  
  160.  
  161.  
  162.  
  163.  

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