Télécharger princi.eso

Retour à la liste

Numérotation des lignes :

princi
  1. C PRINCI SOURCE CB215821 19/07/31 21:16:41 10277
  2. SUBROUTINE PRINCI
  3. C=======================================================================
  4. C
  5. C CALCUL DE CHAMP DE CONTRAINTES PRINCIPALES
  6. C
  7. C
  8. C CHAM2 = PRINCI CHAM1 (CAR1) MODL (MOTCL);
  9. C
  10. C MOTCL = 'SUP ' OU 'INF ' OU 'MOYE' POUR LES COQUES
  11. C OU 'TRID' POUR LES MASSIFS
  12. C CAR1 = objet de type MCHAML de sous type CARACTERISTIQUES
  13. C CHAM1 = objet de type MCHAML de sous type CONTRAINTES
  14. C ou DEFORMATIONS
  15. C MODL = objet de type MMODEL
  16. C CHAM2 = objet de type MCHAML de sous type CONTRAINTES
  17. C PRINCIPALES
  18. C
  19. C Passage au nouveau Chamelem par S.RAMAHANDRY le 21/09/90
  20. C
  21. C=======================================================================
  22. IMPLICIT INTEGER(I-N)
  23. IMPLICIT REAL*8(A-H,O-Z)
  24.  
  25. -INC PPARAM
  26. -INC CCOPTIO
  27. CHARACTER*4 MMM
  28. C
  29. C
  30. C IPMODL MODELE MMODEL
  31. C IPCHE1 MCHAML CONTRAINTES ou DEFORMATIONS
  32. C IPCHE2 MCHAML CARACTERISTIQUES
  33. C IPSTRS MCHAML CONTRAINTES PRINCIPALES
  34. C
  35. IPMODL=0
  36. IRETOU=0
  37. IPCHE1=0
  38. IPCHE2=0
  39. IPMODL=0
  40. IPSTRS=0
  41. IRETOU=0
  42. KER =0
  43. IR =0
  44. MMM =' '
  45. C
  46. C LECTURE D'UN MOT CLEF
  47. C
  48. CALL LIRCHA(MMM,0,IRETOU)
  49. IF(IRETOU.EQ.0) MMM='MOYE'
  50. C
  51. C LECTURE D'UN MODEL
  52. C
  53. CALL LIROBJ('MMODEL ',IPMODL,1,IRETOU)
  54. CALL ACTOBJ('MMODEL ',IPMODL,1)
  55. IF(IERR.NE.0) RETURN
  56.  
  57. C LECTURE D'UN PREMIER MCHAML (CONTRAINTES ou DEFORMATIONS)
  58. CALL LIROBJ('MCHAML ',IPIN,1,IRETOU)
  59. CALL ACTOBJ('MCHAML ',IPIN,1)
  60. IF(IERR.NE.0) RETURN
  61. CALL REDUAF(IPIN,IPMODL,IPCHE1,0,ir,ker)
  62. IF (ir.NE.1) CALL erreur(ker)
  63. IF (IERR.NE.0) RETURN
  64. C
  65. C LECTURE D'UN DEUXIEME MCHAML (CARACTERISTIQUES)
  66. CALL LIROBJ('MCHAML ',IPIN,0,IRETOU)
  67. IF(IERR.NE.0) RETURN
  68. IF(IRETOU .EQ. 1)THEN
  69. CALL ACTOBJ('MCHAML ',IPIN,1)
  70.  
  71. CALL REDUAF(IPIN,IPMODL,IPCHE2,0,ir,ker)
  72. IF (ir.NE.1) CALL erreur(ker)
  73. IF (IERR.NE.0) RETURN
  74. ENDIF
  75. C
  76. C APPEL A PRINPO
  77. C ==============
  78. CALL PRINPO(IPCHE1,MMM,IPCHE2,IPMODL,IPSTRS,IRET)
  79. IF(IRET.NE.0 .AND. IERR.EQ.0) THEN
  80. CALL ACTOBJ('MCHAML ',IPSTRS,1)
  81. CALL ECROBJ('MCHAML ',IPSTRS)
  82. ENDIF
  83.  
  84. END
  85.  
  86.  
  87.  

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