Télécharger princi.eso

Retour à la liste

Numérotation des lignes :

  1. C PRINCI SOURCE CB215821 16/11/03 21:15:13 9155
  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. -INC CCOPTIO
  25. CHARACTER*4 MMM
  26. C
  27. C
  28. C IPMODL MODELE MMODEL
  29. C IPCHE1 MCHAML CONTRAINTES ou DEFORMATIONS
  30. C IPCHE2 MCHAML CARACTERISTIQUES
  31. C IPSTRS MCHAML CONTRAINTES PRINCIPALES
  32. C
  33. IPMODL=0
  34. IRETOU=0
  35. IPCHE1=0
  36. IPCHE2=0
  37. IPMODL=0
  38. IPSTRS=0
  39. IRETOU=0
  40. KER =0
  41. IR =0
  42. MMM =' '
  43. C
  44. C LECTURE D'UN MOT CLEF
  45. C
  46. CALL LIRCHA(MMM,0,IRETOU)
  47. IF(IRETOU.EQ.0) MMM='MOYE'
  48. C
  49. C LECTURE D'UN MODEL
  50. C
  51. CALL LIROBJ('MMODEL',IPMODL,1,IRETOU)
  52. IF(IERR.NE.0) RETURN
  53.  
  54. C LECTURE D'UN PREMIER MCHAML (CONTRAINTES ou DEFORMATIONS)
  55. CALL LIROBJ('MCHAML',IPIN,1,IRETOU)
  56. IF(IERR.NE.0) RETURN
  57. CALL REDUAF(IPIN,IPMODL,IPCHE1,0,ir,ker)
  58. IF (ir.NE.1) CALL erreur(ker)
  59. IF (IERR.NE.0) RETURN
  60. C
  61. C LECTURE D'UN DEUXIEME MCHAML (CARACTERISTIQUES)
  62. CALL LIROBJ('MCHAML',IPIN,0,IRETOU)
  63. IF(IERR.NE.0) RETURN
  64. IF(IRETOU .EQ. 1)THEN
  65. CALL REDUAF(IPIN,IPMODL,IPCHE2,0,ir,ker)
  66. IF (ir.NE.1) CALL erreur(ker)
  67. IF (IERR.NE.0) RETURN
  68. ENDIF
  69. C
  70. C APPEL A PRINPO
  71. C ==============
  72. CALL PRINPO(IPCHE1,MMM,IPCHE2,IPMODL,IPSTRS,IRET)
  73. IF(IRET.NE.0 .AND. IERR.EQ.0) CALL ECROBJ('MCHAML',IPSTRS)
  74. C
  75. RETURN
  76. END
  77.  
  78.  
  79.  

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