Télécharger la plaquette Cast3M


fr en

Télécharger ellp12.eso

Retour à la liste

Numérotation des lignes : 

ellp12

  1. C ELLP12    SOURCE    CHAT      05/01/12    23:35:16     5004
  2.         SUBROUTINE ELLP12(CORRES,XCL,FLAG,NUMERO,MASS,NMAS,
  3.      *                    NP,NP24,NNT,ZA1,ZSM)
  4. C
  5.       IMPLICIT INTEGER(I-N)
  6.         IMPLICIT REAL*8     (A-H,O-Y)
  7.         IMPLICIT COMPLEX*16 (Z)
  8. C
  9. C CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
  10. C
  11. C   OPERATEUR ELFE LAPLACE POUTRE
  12. C
  13. C   PROGRAMME APPELE PAR ELLP11 : CALCUL LA DEUXIEME PARTIE DE LA
  14. C   MATRICE ZA1 ET DE SECOND MEMBRE ZSM
  15. C
  16. C
  17. C
  18. C   PARAMETRES :
  19. C
  20. C   ENTREE :
  21. C
  22. C   CORRES : TABLEAU D'ENTIER DES NUMEROS DES NOEUDS REELS
  23. C            POUR CHAQUE NOEUD FICTIF
  24. C   XCL    : TABLEAU DES VALEURS DES CONDITIONS AUX LIMITES
  25. C   FLAG   : TABLEAU DE POINTEURS SUR XCL
  26. C   NUMERO : TABLEAU DE TRANSFORMATION NUMERO GLOBAL <--> NUMERO LOCAL
  27. C
  28. C   NP     : NOMBRE DE POUTRES
  29. C   NNT    : NOMBRE DE NOEUDS REELS
  30. C
  31. C
  32. C   SORTIES :
  33. C
  34. C   ZA1    : TABLEAU COMPLEXE REPRESENTANT LA MATRICE DE "RIGIDITE"
  35. C   ZSM    : VECTEUR COMPLEXE SECOND MEMBRE
  36. C
  37. C
  38. C   AUTEUR : SAINT-DIZIER
  39. C   DATE   : 22 DECEMBRE 1989 (VERSION DU 03 AVRIL 1990)
  40. C
  41. C CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
  42. C
  43.         COMPLEX*16 ZA1(NP24,*),ZSM(*)
  44.         REAL*8     XCL(12,*)
  45.         INTEGER    FLAG(*),CORRES(*),NUMERO(*),MASS(NNT,4)
  46. C
  47. C **********************************************************************
  48. C         REMPLISSAGE DE LA 2EME PARTIE DE ZA1 (INDEPENDANT DE W)
  49. C         + 2EME PARTIE DU SECOND MEMBRE
  50. C **********************************************************************
  51. C
  52.         ZI = 1.D0
  53. C
  54.         IZA1 = 12 * NP + 1
  55.         NP2  = 2 * NP
  56. C
  57. C
  58. C ---------------------  CONDITION DE NOEUDS IDENTIQUES
  59. C
  60.         DO 10 I = 1 , NNT
  61. C
  62.            ITEST = 0
  63.            J = 0
  64. C
  65. 11         J = J + 1
  66. C
  67.            IF ((CORRES(J).EQ.NUMERO(I)).AND.(ITEST.EQ.0)) THEN
  68.               J0 = J
  69.               ITEST = 1
  70.            ELSE IF ((CORRES(J).EQ.NUMERO(I)).AND.(ITEST.EQ.1)) THEN
  71.               ZA1 (IZA1   , 12*(J0-1)+1 ) = ZI
  72.               ZA1 (IZA1+1 , 12*(J0-1)+2 ) = ZI
  73.               ZA1 (IZA1+2 , 12*(J0-1)+3 ) = ZI
  74.               ZA1 (IZA1+3 , 12*(J0-1)+4 ) = ZI
  75.               ZA1 (IZA1+4 , 12*(J0-1)+5 ) = ZI
  76.               ZA1 (IZA1+5 , 12*(J0-1)+6 ) = ZI
  77. C
  78.               ZA1 (IZA1   , 12*(J-1)+1 ) = -ZI
  79.               ZA1 (IZA1+1 , 12*(J-1)+2 ) = -ZI
  80.               ZA1 (IZA1+2 , 12*(J-1)+3 ) = -ZI
  81.               ZA1 (IZA1+3 , 12*(J-1)+4 ) = -ZI
  82.               ZA1 (IZA1+4 , 12*(J-1)+5 ) = -ZI
  83.               ZA1 (IZA1+5 , 12*(J-1)+6 ) = -ZI
  84.               IZA1 = IZA1 + 6
  85. C
  86.               ENDIF
  87. C
  88.               IF (J.NE.NP2) GOTO 11
  89. C
  90. 10      CONTINUE
  91. C
  92. C
  93.         DO 100 I = 1 , NNT
  94. C
  95. C ---------------------  CONDITION AUX LIMITES SUR LE DEPLACEMENT
  96. C
  97.             IF (FLAG((I-1)*12+1).EQ.I) THEN
  98. C
  99.             J = 0
  100.             ISTOP = 0
  101. 20          J = J + 1
  102.             IF ((CORRES(J).EQ.NUMERO(I)).AND.(ISTOP.EQ.0)) THEN
  103. C
  104.                 ZA1 (IZA1   , 12*(J-1)+1 ) = ZI
  105.                 ZSM ( IZA1  ) = CMPLX(XCL ( 1,I))
  106. C
  107.                 IZA1 = IZA1 + 1
  108.                 ISTOP = 1
  109. C
  110.               ENDIF
  111. C
  112.               IF (J.NE.NP2) GO TO 20
  113. C
  114.             END IF
  115. C
  116. C
  117.             IF (FLAG((I-1)*12+2).EQ.I) THEN
  118. C
  119.             J = 0
  120.             ISTOP = 0
  121. 21          J = J + 1
  122.             IF ((CORRES(J).EQ.NUMERO(I)).AND.(ISTOP.EQ.0)) THEN
  123. C
  124.                 ZA1 (IZA1   , 12*(J-1)+2 ) = ZI
  125.                 ZSM ( IZA1  ) = CMPLX(XCL ( 2,I))
  126. C
  127.                 IZA1 = IZA1 + 1
  128.                 ISTOP = 1
  129. C
  130.               ENDIF
  131. C
  132.               IF (J.NE.NP2) GO TO 21
  133. C
  134.             END IF
  135. C
  136. C
  137.             IF (FLAG((I-1)*12+3).EQ.I) THEN
  138. C
  139.             J = 0
  140.             ISTOP = 0
  141. 22          J = J + 1
  142.             IF ((CORRES(J).EQ.NUMERO(I)).AND.(ISTOP.EQ.0)) THEN
  143. C
  144.                 ZA1 (IZA1   , 12*(J-1)+3 ) = ZI
  145.                 ZSM ( IZA1  ) = CMPLX(XCL ( 3,I))
  146. C
  147.                 IZA1 = IZA1 + 1
  148.                 ISTOP = 1
  149. C
  150.               ENDIF
  151. C
  152.               IF (J.NE.NP2) GO TO 22
  153. C
  154.             END IF
  155. C
  156. C ---------------------  CONDITION AUX LIMITES SUR LA ROTATION
  157. C
  158.             IF (FLAG((I-1)*12+4).EQ.I) THEN
  159. C
  160.               J = 0
  161.               ISTOP = 0
  162. 30            J = J + 1
  163.               IF ((CORRES(J).EQ.NUMERO(I)).AND.(ISTOP.EQ.0)) THEN
  164. C
  165.                 ZA1 (IZA1   , 12*(J-1)+4 ) = ZI
  166.                 ZSM ( IZA1  ) = CMPLX(XCL ( 4,I))
  167. C
  168.                 IZA1 = IZA1 + 1
  169.                 ISTOP = 1
  170. C
  171.               ENDIF
  172. C
  173.               IF (J.NE.NP2) GO TO 30
  174. C
  175.             END IF
  176. C
  177. C
  178.             IF (FLAG((I-1)*12+5).EQ.I) THEN
  179. C
  180.               J = 0
  181.               ISTOP = 0
  182. 31            J = J + 1
  183.               IF ((CORRES(J).EQ.NUMERO(I)).AND.(ISTOP.EQ.0)) THEN
  184. C
  185.                 ZA1 (IZA1   , 12*(J-1)+5 ) = ZI
  186.                 ZSM ( IZA1  ) = CMPLX(XCL ( 5,I))
  187. C
  188.                 IZA1 = IZA1 + 1
  189.                 ISTOP = 1
  190. C
  191.               ENDIF
  192. C
  193.               IF (J.NE.NP2) GO TO 31
  194. C
  195.             END IF
  196. C
  197. C
  198.             IF (FLAG((I-1)*12+6).EQ.I) THEN
  199. C
  200.               J = 0
  201.               ISTOP = 0
  202. 32            J = J + 1
  203.               IF ((CORRES(J).EQ.NUMERO(I)).AND.(ISTOP.EQ.0)) THEN
  204. C
  205.                 ZA1 (IZA1   , 12*(J-1)+6 ) = ZI
  206.                 ZSM ( IZA1  ) = CMPLX(XCL ( 6,I))
  207. C
  208.                 IZA1 = IZA1 + 1
  209.                 ISTOP = 1
  210. C
  211.               ENDIF
  212. C
  213.               IF (J.NE.NP2) GO TO 32
  214. C
  215.             END IF
  216. C
  217. C ---------------------  CONDITION AUX LIMITES SUR LA FORCE
  218. C
  219.             IF (FLAG((I-1)*12+1).EQ.0) THEN
  220. C
  221.               DO 40 J = 1 , 2*NP
  222.                   IF (CORRES(J).EQ.NUMERO(I)) THEN
  223. C
  224.                     ZA1 (IZA1   , 12*(J-1)+7 ) = ZI
  225. C
  226.                   END IF
  227. 40            CONTINUE
  228. C
  229.               ZSM ( IZA1  ) =  CMPLX(XCL(7,I))
  230. C
  231.               IF (NMAS.GT.0) THEN
  232.                 DO 41 I1 = 1 , NMAS
  233.                   IF (CORRES(MASS(I1,1)).EQ.NUMERO(I)) THEN
  234.                     MASS(I1,4) = IZA1
  235.                   END IF
  236. 41              CONTINUE
  237.               END IF
  238. C
  239.               IZA1 = IZA1 + 1
  240. C
  241.             END IF
  242. C
  243. C
  244.             IF (FLAG((I-1)*12+2).EQ.0) THEN
  245. C
  246.               DO 42 J = 1 , 2*NP
  247.                   IF (CORRES(J).EQ.NUMERO(I)) THEN
  248.                     ZA1 (IZA1   , 12*(J-1)+8 ) = ZI
  249. C
  250.                   END IF
  251. 42            CONTINUE
  252. C
  253.               ZSM ( IZA1  ) =  CMPLX(XCL(8,I))
  254. C
  255.               IF (NMAS.GT.0) THEN
  256.                 DO 43 I1 = 1 , NMAS
  257.                   IF (CORRES(MASS(I1,1)).EQ.NUMERO(I)) THEN
  258.                     MASS(I1,4) = IZA1-1
  259.                   END IF
  260. 43              CONTINUE
  261.               END IF
  262. C
  263.               IZA1 = IZA1 + 1
  264. C
  265.             END IF
  266. C
  267. C
  268.             IF (FLAG((I-1)*12+3).EQ.0) THEN
  269. C
  270.               DO 44 J = 1 , 2*NP
  271.                   IF (CORRES(J).EQ.NUMERO(I)) THEN
  272.                     ZA1 (IZA1   , 12*(J-1)+9 ) = ZI
  273. C
  274.                   END IF
  275. 44            CONTINUE
  276. C
  277.               ZSM ( IZA1  ) =  CMPLX(XCL(9,I))
  278. C
  279.               IF (NMAS.GT.0) THEN
  280.                 DO 45 I1 = 1 , NMAS
  281.                   IF (CORRES(MASS(I1,1)).EQ.NUMERO(I)) THEN
  282.                     MASS(I1,4) = IZA1-2
  283.                   END IF
  284. 45              CONTINUE
  285.               END IF
  286. C
  287.               IZA1 = IZA1 + 1
  288. C
  289.             END IF
  290. C
  291. C ---------------------  CONDITION AUX LIMITES SUR LE MOMENT
  292. C
  293.             IF (FLAG((I-1)*12+4).EQ.0) THEN
  294. C
  295.               DO 50 J = 1 , 2*NP
  296.                   IF (CORRES(J).EQ.NUMERO(I)) THEN
  297.                     ZA1 (IZA1   , 12*(J-1)+10) = ZI
  298.                   END IF
  299. 50            CONTINUE
  300. C
  301.               ZSM ( IZA1  ) =  CMPLX(XCL(10,I))
  302. C
  303.               IF (NMAS.GT.0) THEN
  304.                 DO 51 I1 = 1 , NMAS
  305.                   IF (CORRES(MASS(I1,1)).EQ.NUMERO(I)) THEN
  306.                     MASS(I1,4) = IZA1-3
  307.                   END IF
  308. 51              CONTINUE
  309.               END IF
  310. C
  311.               IZA1 = IZA1 + 1
  312. C
  313.             END IF
  314. C
  315. C
  316.             IF (FLAG((I-1)*12+5).EQ.0) THEN
  317. C
  318.               DO 52 J = 1 , 2*NP
  319.                   IF (CORRES(J).EQ.NUMERO(I)) THEN
  320.                     ZA1 (IZA1   , 12*(J-1)+11) = ZI
  321.                   END IF
  322. 52            CONTINUE
  323. C
  324.               ZSM ( IZA1  ) =  CMPLX(XCL(11,I))
  325. C
  326.               IF (NMAS.GT.0) THEN
  327.                 DO 53 I1 = 1 , NMAS
  328.                   IF (CORRES(MASS(I1,1)).EQ.NUMERO(I)) THEN
  329.                     MASS(I1,4) = IZA1-4
  330.                   END IF
  331. 53              CONTINUE
  332.               END IF
  333. C
  334.               IZA1 = IZA1 + 1
  335. C
  336.             END IF
  337. C
  338. C
  339.             IF (FLAG((I-1)*12+6).EQ.0) THEN
  340. C
  341.               DO 54 J = 1 , 2*NP
  342.                   IF (CORRES(J).EQ.NUMERO(I)) THEN
  343.                     ZA1 (IZA1   , 12*(J-1)+12) = ZI
  344.                   END IF
  345. 54            CONTINUE
  346. C
  347.               ZSM ( IZA1  ) =  CMPLX(XCL(12,I))
  348. C
  349.               IF (NMAS.GT.0) THEN
  350.                 DO 55 I1 = 1 , NMAS
  351.                   IF (CORRES(MASS(I1,1)).EQ.NUMERO(I)) THEN
  352.                     MASS(I1,4) = IZA1-5
  353.                   END IF
  354. 55              CONTINUE
  355.               END IF
  356. C
  357.               IZA1 = IZA1 + 1
  358. C
  359.             END IF
  360. C
  361. C
  362. 100     CONTINUE
  363. C
  364.         END
  365.  
  366.  

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