Télécharger la plaquette Cast3M


fr en

Télécharger mrem.eso

Retour à la liste

Numérotation des lignes : 

mrem

  1. C MREM      SOURCE    CB215821  20/11/25    13:34:31     10792          
  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.  
  84.  

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