Télécharger ksigma.eso

Retour à la liste

Numérotation des lignes :

  1. C KSIGMA SOURCE CB215821 16/12/05 21:39:54 9237
  2. SUBROUTINE KSIGMA
  3. C_______________________________________________________________________
  4. C
  5. C OPERATEUR KSIGMA
  6. C
  7. C Construction de la matrice de raideur geometrique @ partir d'un
  8. C MCHAML de contraintes
  9. C
  10. C Nouvelle syntaxe:
  11. C _________________
  12. C
  13. C RIG1 =KSIG MOD1 SIG1 ( CAR1 ) 'FLAM';
  14. C
  15. C MOD1=OBJET DE TYPE MMODEL
  16. C SIG1=MCHAML de CONTRAINTES
  17. C CAR1=MCHAML de CARACTERISTIQUES
  18. C
  19. C RIG1=Matrice KSIGMA de sous type RIGIDITE
  20. C
  21. C Passage aux nouveaux MCHAMLs par I.Monnier le 14.06.90
  22. C
  23. C_______________________________________________________________________
  24. C
  25. IMPLICIT INTEGER(I-N)
  26. IMPLICIT REAL*8(A-H,O-Z)
  27. C
  28. -INC CCOPTIO
  29. C
  30. IPCHA1=0
  31. IPCHA2=0
  32. IPCHE1=0
  33. IPCHE2=0
  34. IFLAM=0
  35. C
  36. CALL LIRMOT('FLAM',1,IFLAM,0)
  37. IF (IERR.NE.0) RETURN
  38. C
  39. CALL LIROBJ('MMODEL',IPMODL,1,IRT1)
  40. IF (IERR.NE.0) RETURN
  41. C
  42. CALL LIROBJ('MCHAML',IPIN,1,IRT1)
  43. IF (IERR.NE.0) RETURN
  44. CALL REDUAF(IPIN,IPMODL,IPCHA1,0,IR,KER)
  45. IF(IR .NE. 1) CALL ERREUR(KER)
  46. IF(IERR .NE. 0) RETURN
  47. C
  48. CALL LIROBJ('MCHAML',IPIN,0,IRT1)
  49. IF (IERR.NE.0) RETURN
  50. IPCHA2=0
  51. IF (IRT1 .EQ. 1) THEN
  52. CALL REDUAF(IPIN,IPMODL,IPCHA2,0,IR,KER)
  53. IF(IR .NE. 1) CALL ERREUR(KER)
  54. IF(IERR .NE. 0) RETURN
  55. ENDIF
  56. C
  57. CALL RNGCHA(IPCHA1,IPCHA2,'CONTRAINTES',
  58. 1 'CARACTERISTIQUES',IPCHE1,IPCHE2)
  59. C
  60. IF(IPCHE1.EQ.0)THEN
  61. MOTERR(1:16)='CONTRAINTES'
  62. CALL ERREUR(565)
  63. RETURN
  64. ENDIF
  65. C
  66. CALL KSIGMP(IPMODL,IPCHE1,IPCHE2,IFLAM,IPRIGG)
  67. C
  68. IF (IERR.EQ.0.AND.IPRIGG.NE.0) CALL ECROBJ('RIGIDITE',IPRIGG)
  69. C
  70. RETURN
  71. END
  72.  
  73.  
  74.  

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