Télécharger ksigma.eso

Retour à la liste

Numérotation des lignes :

  1. C KSIGMA SOURCE PV 20/03/30 21:20:39 10567
  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.  
  29. -INC PPARAM
  30. -INC CCOPTIO
  31. -INC SMCOORD
  32. C
  33. segact mcoord
  34. IPCHA1=0
  35. IPCHA2=0
  36. IPCHE1=0
  37. IPCHE2=0
  38. IFLAM=0
  39. C
  40. CALL LIRMOT('FLAM',1,IFLAM,0)
  41. IF (IERR.NE.0) RETURN
  42. C
  43. CALL LIROBJ('MMODEL',IPMODL,1,IRT1)
  44. IF (IERR.NE.0) RETURN
  45. CALL ACTOBJ('MMODEL ',IPMODL,1)
  46. C
  47. CALL LIROBJ('MCHAML',IPIN,1,IRT1)
  48. IF (IERR.NE.0) RETURN
  49. CALL ACTOBJ('MCHAML ',IPIN,1)
  50. CALL REDUAF(IPIN,IPMODL,IPCHA1,0,IR,KER)
  51. IF(IR .NE. 1) CALL ERREUR(KER)
  52. IF(IERR .NE. 0) RETURN
  53. C
  54. CALL LIROBJ('MCHAML',IPIN,0,IRT1)
  55. IF (IERR.NE.0) RETURN
  56. IPCHA2=0
  57. IF (IRT1 .EQ. 1) THEN
  58. CALL ACTOBJ('MCHAML ',IPIN,1)
  59. CALL REDUAF(IPIN,IPMODL,IPCHA2,0,IR,KER)
  60. IF(IR .NE. 1) CALL ERREUR(KER)
  61. IF(IERR .NE. 0) RETURN
  62. ENDIF
  63. C
  64. CALL RNGCHA(IPCHA1,IPCHA2,'CONTRAINTES',
  65. 1 'CARACTERISTIQUES',IPCHE1,IPCHE2)
  66. C
  67. IF(IPCHE1.EQ.0)THEN
  68. MOTERR(1:16)='CONTRAINTES'
  69. CALL ERREUR(565)
  70. RETURN
  71. ENDIF
  72. C
  73. CALL KSIGMP(IPMODL,IPCHE1,IPCHE2,IFLAM,IPRIGG)
  74. C
  75. IF (IERR.EQ.0.AND.IPRIGG.NE.0) CALL ECROBJ('RIGIDITE',IPRIGG)
  76. C
  77. RETURN
  78. END
  79.  
  80.  
  81.  
  82.  
  83.  

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