Télécharger kp.eso

Retour à la liste

Numérotation des lignes :

  1. C KP SOURCE MB234859 16/12/14 21:15:18 9252
  2. SUBROUTINE KP
  3. C_______________________________________________________________________
  4. C
  5. C OPERATEUR KP
  6. C
  7. C Construction de la matrice des pressions associees à la linearisation
  8. C des actions de pression sur la géometrie deformée.
  9. C
  10. C
  11. C Nouvelle syntaxe:
  12. C _________________
  13. *
  14. * 1ere option
  15. C
  16. C RIG1 =KP MODL1 PRES1 ('FLAM') ('ASYM');
  17. C
  18. C MODL1 = objet modele de type MMODEL
  19. C
  20. C PRES1=MCHAML ou CHPOINT de PRESSION
  21. C
  22. C 'FLAM' = MOT-CLE N{CESSAIRE SI L'ON VEUT UTILISER LA MATRICE
  23. C POUR FAIRE UNE ANALYSE DE FLAMBAGE.
  24. C
  25. C RIG1=Matrice KP de sous type RIGIDITE
  26. *
  27. * 2ème option
  28. *
  29. * RIG1=KP MODL1 RG (VEC1) ('ASYM')
  30. *
  31. * MOD1 : MODELE DE CALCUL, TYPE MMODEL
  32. * RG : COEFFICIENT MULTIPLICATIF DE RIGIDITE
  33. * VEC1 : vecteur (type point) qui definit le sens de
  34. * la "pesanTeur"(facultatif). S'il n'est pas
  35. * fourni vec1 est consideré suivant la normale
  36. * à chaque élément.
  37. * ASYM : mot clé pour avoir la matrice dissymetrique
  38. *
  39. * la deuxième version CALCULE LA MATRICE DE RAIDEUR LIEE A LA
  40. * VARIATION DE PRESSION DUE AU MOUVEMENT VERTICAL DANS UN CHAMP DE
  41. * PESANTEUR (on peut l'utiliser dans d'autres cas aussi, p.ex. plaque
  42. * sur sol elastique)
  43. C
  44. C ces matrices sont en general non symetriques, mais leur combinaison peut
  45. c etre symetrique (p.ex. en cas de pression hydrostatique). C'est pourquoi
  46. c nous voulons avoir la possibilité de construire les matrices symetriques
  47. c ou disymetriques.
  48. C
  49. C
  50. C CODE COMBESCURE JANV 87
  51. C
  52. C Passage aux nouveaux CHAMELEMs par P.DOWLATYARI le 5/4/91
  53. C
  54. C Ajout de la 2ème option I. Politopoulos Aout 1995
  55. C_______________________________________________________________________
  56. C
  57. IMPLICIT INTEGER(I-N)
  58. IMPLICIT REAL*8(A-H,O-Z)
  59. C
  60. -INC CCOPTIO
  61. C
  62. C
  63. C
  64. IPMODL=0
  65. IPCHP=0
  66. IPRIG=0
  67. IFLAM=0
  68. ICHA=0
  69. IDISS=0
  70.  
  71. C
  72. CALL LIRMOT('FLAM',1,IFLAM,0)
  73. IF (IERR.NE.0) CALL ERREUR(7)
  74. C
  75. CALL LIRMOT('ASYM',1,IDISS,0)
  76. IF (IERR.NE.0) CALL ERREUR(7)
  77.  
  78. C
  79. CALL LIROBJ('MMODEL',IPMODL,1,IRT1)
  80. IF (IERR.NE.0) RETURN
  81.  
  82. CALL LIROBJ('MCHAML',IPIN,0,IRT2)
  83. IF (IERR.NE.0) RETURN
  84. IF(IRT2.EQ.1) THEN
  85. CALL REDUAF(IPIN,IPMODL,IPCHP,0,IR,KER)
  86. IF(IR .NE. 1) CALL ERREUR(KER)
  87. IF(IERR .NE. 0) RETURN
  88. ICHA=1
  89. ELSE
  90. CALL LIROBJ('CHPOINT ',IPCHP,0,IRT3)
  91. IF (IERR.NE.0) RETURN
  92. IF(IRT3.EQ.0) GO TO 10
  93. ENDIF
  94. C
  95. CALL KPRES(IPMODL,IPCHP,ICHA,ICONV,IFLAM,IDISS,IPRIG,IRET)
  96. GOTO 20
  97. C
  98. *** matrice associee à un gradient de pression imposée
  99. 10 CONTINUE
  100. CALL LIRREE(XRG,1,IRET)
  101. IF (IRET.EQ.0) GOTO 666
  102. CALL LIROBJ('POINT',IPOIN1,0,IRET2)
  103.  
  104. CALL KGRAV1 (IPMODL,XRG,IPOIN1,IPRIG,IRET,IDISS,IFLAM)
  105. *
  106. * ecriture de la rigidite
  107.  
  108. 20 CONTINUE
  109. IF (IRET.EQ.1) THEN
  110. CALL ECROBJ('RIGIDITE',IPRIG)
  111. ENDIF
  112.  
  113. 666 CONTINUE
  114. RETURN
  115. END
  116.  
  117.  
  118.  
  119.  
  120.  
  121.  
  122.  
  123.  

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