Télécharger mrem.eso

Retour à la liste

Numérotation des lignes :

  1. C MREM SOURCE PV 19/12/09 21:15:09 10422
  2. SUBROUTINE MREM
  3. ************************************************************************
  4. * remontee de la sotution complete apres resolution a partir
  5. * d une matrice condensee par CMCT ( hors de resou )
  6. *
  7. * Syntaxe :
  8. * chpo3 = MREM chpo1 (rig1 et rig2) chpo2 ;
  9. *
  10. * chpo1 solution reduite sur les ddl non elimines
  11. * rig1 rigidites initiale (hors dependances )
  12. * rig2 rigidites de dependances
  13. *
  14. * chpo3 solution complete en deplacements et LX
  15. *
  16. *************************************************************************
  17. IMPLICIT INTEGER(I-N)
  18. IMPLICIT REAL*8(A-H,O-Z)
  19.  
  20. -INC SMRIGID
  21.  
  22. -INC PPARAM
  23. -INC CCOPTIO
  24. -INC SMCHPOI
  25. -INC SMLCHPO
  26. segment idemem(0)
  27. segment ideme0(idemem(/1),30)
  28. segment ideme1(idemem(/1),30)
  29.  
  30. noen=1
  31. ipt8=0
  32. segini idemem
  33. CALL LIROBJ('CHPOINT',mchpoi,1,IRETOU)
  34. IF(IERR.NE.0) GO TO 5000
  35. idemem(**)=mchpoi
  36. segini ideme0,ideme1
  37. CALL LIROBJ('LISTCHPO',mlchpo,1,IRETOU)
  38. IF(IERR.NE.0) GO TO 5000
  39. CALL LIROBJ('LISTCHPO',mlchp1,1,IRETOU)
  40. IF(IERR.NE.0) GO TO 5000
  41. segact mlchpo,mlchp1
  42. if=mlchpo.ichpoi(/1)
  43. if (if.ne.mlchp1.ichpoi(/1)) call erreur(5)
  44. do 1000 i=1,if
  45. ideme0(1,i)=mlchpo.ichpoi(i)
  46. ideme1(1,i)=mlchp1.ichpoi(i)
  47. 1000 continue
  48.  
  49. CALL LIROBJ('RIGIDITE',mrigid,1,IRETOU)
  50. IF(IERR.NE.0) GO TO 5000
  51. do 2010 ifois=1,30
  52. segact mrigid
  53. mrigid=jrsup
  54. if (mrigid.eq.0) goto 2011
  55. segact mrigid
  56. isouci=1
  57. iverif=0
  58. call resour(idemem,ideme0,ideme1,mrigid,if,noen,ipt8,
  59. > isouci,iverif)
  60. if=if-1
  61. 2010 continue
  62. 2011 continue
  63. if (if.ne.0) call erreur(5)
  64. iret=idemem(1)
  65. call ecrobj('CHPOINT',iret)
  66.  
  67. 5000 continue
  68. return
  69. end
  70.  
  71.  
  72.  
  73.  
  74.  
  75.  
  76.  
  77.  
  78.  
  79.  
  80.  
  81.  
  82.  
  83.  

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