Télécharger ella12.eso

Retour à la liste

Numérotation des lignes :

  1. C ELLA12 SOURCE CHAT 05/01/12 23:34:04 5004
  2. SUBROUTINE ELLA12(CORRES ,XCL ,FLAG ,NUMERO, CARACT ,COOR,
  3. & MASS ,NMAS ,IRAILO ,VALRAI, NRAIDE ,
  4. & NP,NP28,NNT,ZA1,ZSM,ZS,XPI,IMP)
  5. C
  6. C OPERATEUR ELFE LAPLACE ACOU
  7. C
  8. C PROGRAMME APPELE PAR ELLA11 : CALCUL LA DEUXIEME PARTIE DE LA
  9. C MATRICE ZA1 ET DE SECOND MEMBRE ZSM
  10. C
  11. C PARAMETRES :
  12. C
  13. C ENTREE :
  14. C
  15. C CORRES : TABLEAU D'ENTIER DES NUMEROS DES NOEUDS REELS
  16. C POUR CHAQUE NOEUD FICTIF
  17. C XCL : TABLEAU DES VALEURS DES CONDITIONS AUX LIMITES
  18. C FLAG : TABLEAU DE POINTEURS SUR XCL
  19. C NUMERO : TABLEAU DE TRANSFORMATION NUMERO GLOBAL <--> NUMERO LOCAL
  20. C CARACT : TABLEAU DES CARACTERISTIQUES DES POUTRES
  21. C COOR : TABLEAU DES COORDONNEES DES NOEUDS FICTIFS
  22. C
  23. C NP : NOMBRE DE POUTRES
  24. C NNT : NOMBRE DE NOEUDS REELS
  25. C ZS : VALEUR DE S = S0 + I W
  26. C XPI : VALEUR PRECISE DE PI
  27. C
  28. C
  29. C SORTIES :
  30. C
  31. C ZA1 : TABLEAU COMPLEXE REPRESENTANT LA MATRICE DE "RIGIDITE"
  32. C ZSM : VECTEUR COMPLEXE SECOND MEMBRE
  33. C
  34. C
  35. C AUTEURS : SAINT-DIZIER ET GORCY
  36.  
  37. C DATE : 29 JANVIER 1991
  38. C MODIFICATION : 26 JUIN 1991
  39. C
  40. C
  41. IMPLICIT INTEGER(I-N)
  42. COMPLEX*16 ZA1,ZSM,ZB,ZI,ZS
  43. REAL*8 XCL,CARACT,COOR,VALRAI,XPI,X1,Y1,Z1,X2,Y2,Z2,
  44. & XI1,XI2,XI3,XL,CSEF,CROF
  45. INTEGER FLAG,CORRES,NUMERO,MASS,IRAILO,J0,JCOLO
  46. C
  47. DIMENSION ZA1(NP28,*),ZSM(*)
  48. DIMENSION XCL(17,*),CARACT(10,*),COOR(3,*),VALRAI(6,*),
  49. & FLAG(*),CORRES(*),NUMERO(*),MASS(4,*),IRAILO(4,*)
  50. C
  51. C REMPLISSAGE DE LA 2EME PARTIE DE ZA1 (INDEPENDANT DE W)
  52. C + 2EME PARTIE DU SECOND MEMBRE
  53. C
  54. J0=0
  55. JCOLO=0
  56. ZI = CMPLX(1.D0,0.D0)
  57. C
  58. IZA1 = 14 * NP + 1
  59. NP2 = 2 * NP
  60. C
  61. C --------------------- CONDITION DE NOEUDS IDENTIQUES
  62. C (CONTINUITE DU DEPLACEMENT ET DE LA PRESSION)
  63. C
  64. DO 10 I = 1 , NNT
  65. C
  66. ITEST = 0
  67. J = 0
  68. C
  69. 11 J = J + 1
  70. C
  71. IF ((CORRES(J).EQ.NUMERO(I)).AND.(ITEST.EQ.0)) THEN
  72. J0 = J
  73. ITEST = 1
  74. ELSE IF ((CORRES(J).EQ.NUMERO(I)).AND.(ITEST.EQ.1)) THEN
  75. C
  76. ZA1 (IZA1 , 14*(J0-1)+1 ) = -ZI
  77. ZA1 (IZA1+1 , 14*(J0-1)+2 ) = -ZI
  78. ZA1 (IZA1+2 , 14*(J0-1)+3 ) = -ZI
  79. ZA1 (IZA1+3 , 14*(J0-1)+4 ) = -ZI
  80. ZA1 (IZA1+4 , 14*(J0-1)+5 ) = -ZI
  81. ZA1 (IZA1+5 , 14*(J0-1)+6 ) = -ZI
  82. ZA1 (IZA1+6 , 14*(J0-1)+13) = -ZI
  83. C
  84. ZA1 (IZA1 , 14*(J-1)+1 ) = ZI
  85. ZA1 (IZA1+1 , 14*(J-1)+2 ) = ZI
  86. ZA1 (IZA1+2 , 14*(J-1)+3 ) = ZI
  87. ZA1 (IZA1+3 , 14*(J-1)+4 ) = ZI
  88. ZA1 (IZA1+4 , 14*(J-1)+5 ) = ZI
  89. ZA1 (IZA1+5 , 14*(J-1)+6 ) = ZI
  90. ZA1 (IZA1+6 , 14*(J-1)+13) = ZI
  91. C
  92. C INTRODUCTION DES SOUPLESSES DANS LA MATRICE ZA1
  93. C
  94. IF (NRAIDE.GT.0) THEN
  95. C
  96. DO 15 IM1 = 1 , NRAIDE
  97. C
  98. IF (IRAILO(4,IM1).EQ.2) THEN
  99. C
  100. C COUDE
  101. C
  102. JCOLO=(14*(IRAILO(2,IM1)-1) ) + 7
  103. ELSE IF (IRAILO(4,IM1).EQ.3) THEN
  104. C
  105. C TE
  106. C
  107. JCOLO=(14*(IRAILO(1,IM1)-1) ) + 7
  108. ELSE
  109. IF (IMP.NE.0) THEN
  110. WRITE(IMP,*)
  111. ENDIF
  112. 600 FORMAT(1X,'IL Y A UNE ERREUR')
  113. END IF
  114. C
  115. IF ( (IRAILO(2,IM1).EQ.J0).OR.(IRAILO(1,IM1).EQ.J) ) THEN
  116. C
  117. ZA1 (IZA1 , JCOLO ) = CMPLX( (-1.D0)**IRAILO(4,IM1)/
  118. & VALRAI(1,IM1) )
  119. ZA1 (IZA1+1, JCOLO + 1 ) = CMPLX( (-1.D0)**IRAILO(4,IM1)/
  120. & VALRAI(2,IM1) )
  121. ZA1 (IZA1+2, JCOLO + 2 ) = CMPLX( (-1.D0)**IRAILO(4,IM1)/
  122. & VALRAI(3,IM1) )
  123. ZA1 (IZA1+3, JCOLO + 3 ) = CMPLX( (-1.D0)**IRAILO(4,IM1)/
  124. & VALRAI(4,IM1) )
  125. ZA1 (IZA1+4, JCOLO + 4 ) = CMPLX( (-1.D0)**IRAILO(4,IM1)/
  126. & VALRAI(5,IM1) )
  127. ZA1 (IZA1+5, JCOLO + 5 ) = CMPLX( (-1.D0)**IRAILO(4,IM1)/
  128. & VALRAI(6,IM1) )
  129.  
  130. END IF
  131. C
  132. 15 CONTINUE
  133. END IF
  134. C
  135. ZSM ( IZA1+6 ) = CMPLX(XCL ( 13,I))
  136. C
  137. IZA1 = IZA1 + 7
  138. C
  139. END IF
  140. C
  141. IF (J.NE.NP2) GOTO 11
  142. C
  143. 10 CONTINUE
  144. C
  145. DO 100 I = 1 , NNT
  146. C
  147. C --------------------- CONDITION AUX LIMITES SUR LE DEPLACEMENT
  148. C
  149. IF (FLAG((I-1)*17+1).EQ.I) THEN
  150. C
  151. J = 0
  152. ISTOP = 0
  153. 20 J = J + 1
  154. IF ((CORRES(J).EQ.NUMERO(I)).AND.(ISTOP.EQ.0)) THEN
  155. C
  156. ZA1 (IZA1 , 14*(J-1)+1 ) = ZI
  157. ZSM ( IZA1 ) = CMPLX(XCL ( 1,I))
  158. C
  159. IZA1 = IZA1 + 1
  160. ISTOP = 1
  161. C
  162. END IF
  163. C
  164. IF (J.NE.NP2) GO TO 20
  165. C
  166. END IF
  167. C
  168. IF (FLAG((I-1)*17+2).EQ.I) THEN
  169. C
  170. J = 0
  171. ISTOP = 0
  172. 21 J = J + 1
  173. IF ((CORRES(J).EQ.NUMERO(I)).AND.(ISTOP.EQ.0)) THEN
  174. C
  175. ZA1 (IZA1 , 14*(J-1)+2 ) = ZI
  176. ZSM ( IZA1 ) = CMPLX(XCL ( 2,I))
  177. C
  178. IZA1 = IZA1 + 1
  179. ISTOP = 1
  180. C
  181. END IF
  182. C
  183. IF (J.NE.NP2) GO TO 21
  184. C
  185. END IF
  186. C
  187. C
  188. IF (FLAG((I-1)*17+3).EQ.I) THEN
  189. C
  190. J = 0
  191. ISTOP = 0
  192. 22 J = J + 1
  193. IF ((CORRES(J).EQ.NUMERO(I)).AND.(ISTOP.EQ.0)) THEN
  194. C
  195. ZA1 (IZA1 , 14*(J-1)+3 ) = ZI
  196. ZSM ( IZA1 ) = CMPLX(XCL ( 3,I))
  197. C
  198. IZA1 = IZA1 + 1
  199. ISTOP = 1
  200. C
  201. END IF
  202. C
  203. IF (J.NE.NP2) GO TO 22
  204. C
  205. END IF
  206. C
  207. C --------------------- CONDITION AUX LIMITES SUR LA ROTATION
  208. C
  209. IF (FLAG((I-1)*17+4).EQ.I) THEN
  210. C
  211. J = 0
  212. ISTOP = 0
  213. 30 J = J + 1
  214. IF ((CORRES(J).EQ.NUMERO(I)).AND.(ISTOP.EQ.0)) THEN
  215. C
  216. ZA1 (IZA1 , 14*(J-1)+4 ) = ZI
  217. ZSM ( IZA1 ) = CMPLX(XCL ( 4,I))
  218. C
  219. IZA1 = IZA1 + 1
  220. ISTOP = 1
  221. C
  222. END IF
  223. C
  224. IF (J.NE.NP2) GO TO 30
  225. C
  226. END IF
  227. C
  228. IF (FLAG((I-1)*17+5).EQ.I) THEN
  229. C
  230. J = 0
  231. ISTOP = 0
  232. 31 J = J + 1
  233. IF ((CORRES(J).EQ.NUMERO(I)).AND.(ISTOP.EQ.0)) THEN
  234. C
  235. ZA1 (IZA1 , 14*(J-1)+5 ) = ZI
  236. ZSM ( IZA1 ) = CMPLX(XCL ( 5,I))
  237. C
  238. IZA1 = IZA1 + 1
  239. ISTOP = 1
  240. C
  241. END IF
  242. C
  243. IF (J.NE.NP2) GO TO 31
  244. C
  245. END IF
  246. C
  247. IF (FLAG((I-1)*17+6).EQ.I) THEN
  248. C
  249. J = 0
  250. ISTOP = 0
  251. 32 J = J + 1
  252. IF ((CORRES(J).EQ.NUMERO(I)).AND.(ISTOP.EQ.0)) THEN
  253. C
  254. ZA1 (IZA1 , 14*(J-1)+6 ) = ZI
  255. ZSM ( IZA1 ) = CMPLX(XCL ( 6,I))
  256. C
  257. IZA1 = IZA1 + 1
  258. ISTOP = 1
  259. C
  260. END IF
  261. C
  262. IF (J.NE.NP2) GO TO 32
  263. C
  264. END IF
  265. C
  266. C --------------------- CONDITION AUX LIMITES SUR LA FORCE
  267. C
  268. IF (FLAG((I-1)*17+1).EQ.0) THEN
  269. C
  270. DO 40 J = 1 , 2*NP
  271. IF (CORRES(J).EQ.NUMERO(I)) THEN
  272. C
  273. INP = INT((J+1)/2)
  274. CSEF = XPI * CARACT(4,INP) * CARACT(4,INP)
  275. N1 = 2*INP-1
  276. N2 = 2*INP
  277. X1 = COOR (1,N1)
  278. Y1 = COOR (2,N1)
  279. Z1 = COOR (3,N1)
  280. X2 = COOR (1,N2)
  281. Y2 = COOR (2,N2)
  282. Z2 = COOR (3,N2)
  283. XL = SQRT((X2-X1)**2 + (Y2-Y1)**2 + (Z2-Z1)**2)
  284. C
  285. XI1 = (X2-X1)/XL
  286. IF (J.EQ.N2) THEN
  287. XI1 = -XI1
  288. END IF
  289. C
  290. ZA1 (IZA1 , 14*(J-1)+7 ) = ZI
  291. ZA1 (IZA1 , 14*(J-1)+13) = -CMPLX( CSEF*XI1 )
  292. C
  293. END IF
  294. 40 CONTINUE
  295. C
  296. ZSM ( IZA1 ) = -CMPLX(XCL(7,I))
  297. C
  298. IF (NMAS.GT.0) THEN
  299. DO 41 I1 = 1 , NMAS
  300. IF (CORRES(MASS(1,I1)).EQ.NUMERO(I)) THEN
  301. MASS(4,I1) = IZA1
  302. END IF
  303. 41 CONTINUE
  304. END IF
  305. C
  306. IZA1 = IZA1 + 1
  307. C
  308. END IF
  309. C
  310. C
  311. IF (FLAG((I-1)*17+2).EQ.0) THEN
  312. C
  313. DO 42 J = 1 , 2*NP
  314. IF (CORRES(J).EQ.NUMERO(I)) THEN
  315. C
  316. INP = INT((J+1)/2)
  317. CSEF = XPI * CARACT(4,INP) * CARACT(4,INP)
  318. N1 = 2*INP-1
  319. N2 = 2*INP
  320. X1 = COOR (1,N1)
  321. Y1 = COOR (2,N1)
  322. Z1 = COOR (3,N1)
  323. X2 = COOR (1,N2)
  324. Y2 = COOR (2,N2)
  325. Z2 = COOR (3,N2)
  326. XL = SQRT((X2-X1)**2 + (Y2-Y1)**2 + (Z2-Z1)**2)
  327. C
  328. XI2 = (Y2-Y1)/XL
  329. IF (J.EQ.N2) THEN
  330. XI2 = -XI2
  331. END IF
  332. C
  333. ZA1 (IZA1 , 14*(J-1)+8 ) = ZI
  334. ZA1 (IZA1 , 14*(J-1)+13) = -CMPLX( CSEF*XI2 )
  335. C
  336. END IF
  337. 42 CONTINUE
  338. C
  339. ZSM ( IZA1 ) = -CMPLX(XCL(8,I))
  340. C
  341. IF (NMAS.GT.0) THEN
  342. DO 43 I1 = 1 , NMAS
  343. IF (CORRES(MASS(1,I1)).EQ.NUMERO(I)) THEN
  344. MASS(4,I1) = IZA1-1
  345. END IF
  346. 43 CONTINUE
  347. END IF
  348. C
  349. IZA1 = IZA1 + 1
  350. C
  351. END IF
  352. C
  353. C
  354. IF (FLAG((I-1)*17+3).EQ.0) THEN
  355. C
  356. DO 44 J = 1 , 2*NP
  357. IF (CORRES(J).EQ.NUMERO(I)) THEN
  358. C
  359. INP = INT((J+1)/2)
  360. CSEF = XPI * CARACT(4,INP) * CARACT(4,INP)
  361. N1 = 2*INP-1
  362. N2 = 2*INP
  363. X1 = COOR (1,N1)
  364. Y1 = COOR (2,N1)
  365. Z1 = COOR (3,N1)
  366. X2 = COOR (1,N2)
  367. Y2 = COOR (2,N2)
  368. Z2 = COOR (3,N2)
  369. XL = SQRT((X2-X1)**2 + (Y2-Y1)**2 + (Z2-Z1)**2)
  370. C
  371. XI3 = (Z2-Z1)/XL
  372. IF (J.EQ.N2) THEN
  373. XI3 = -XI3
  374. END IF
  375. C
  376. ZA1 (IZA1 , 14*(J-1)+9 ) = ZI
  377. ZA1 (IZA1 , 14*(J-1)+13) = -CMPLX( CSEF*XI3 )
  378. C
  379. END IF
  380. 44 CONTINUE
  381. C
  382. ZSM ( IZA1 ) = -CMPLX(XCL(9,I))
  383. C
  384. IF (NMAS.GT.0) THEN
  385. DO 45 I1 = 1 , NMAS
  386. IF (CORRES(MASS(1,I1)).EQ.NUMERO(I)) THEN
  387. MASS(4,I1) = IZA1-2
  388. END IF
  389. 45 CONTINUE
  390. END IF
  391. C
  392. IZA1 = IZA1 + 1
  393. C
  394. END IF
  395. C
  396. C --------------------- CONDITION AUX LIMITES SUR LE MOMENT
  397. C
  398. IF (FLAG((I-1)*17+4).EQ.0) THEN
  399. C
  400. DO 50 J = 1 , 2*NP
  401. IF (CORRES(J).EQ.NUMERO(I)) THEN
  402. ZA1 (IZA1 , 14*(J-1)+10) = ZI
  403. END IF
  404. 50 CONTINUE
  405. C
  406. ZSM ( IZA1 ) = -CMPLX(XCL(10,I))
  407. C
  408. IF (NMAS.GT.0) THEN
  409. DO 51 I1 = 1 , NMAS
  410. IF (CORRES(MASS(1,I1)).EQ.NUMERO(I)) THEN
  411. MASS(4,I1) = IZA1-3
  412. END IF
  413. 51 CONTINUE
  414. END IF
  415. C
  416. IZA1 = IZA1 + 1
  417. C
  418. END IF
  419. C
  420. C
  421. IF (FLAG((I-1)*17+5).EQ.0) THEN
  422. C
  423. DO 52 J = 1 , 2*NP
  424. IF (CORRES(J).EQ.NUMERO(I)) THEN
  425. ZA1 (IZA1 , 14*(J-1)+11) = ZI
  426. END IF
  427. 52 CONTINUE
  428. C
  429. ZSM ( IZA1 ) = -CMPLX(XCL(11,I))
  430. C
  431. IF (NMAS.GT.0) THEN
  432. DO 53 I1 = 1 , NMAS
  433. IF (CORRES(MASS(1,I1)).EQ.NUMERO(I)) THEN
  434. MASS(4,I1) = IZA1-4
  435. END IF
  436. 53 CONTINUE
  437. END IF
  438. C
  439. IZA1 = IZA1 + 1
  440. C
  441. END IF
  442. C
  443. C
  444. IF (FLAG((I-1)*17+6).EQ.0) THEN
  445. C
  446. DO 54 J = 1 , 2*NP
  447. IF (CORRES(J).EQ.NUMERO(I)) THEN
  448. ZA1 (IZA1 , 14*(J-1)+12) = ZI
  449. END IF
  450. 54 CONTINUE
  451. C
  452. ZSM ( IZA1 ) = -CMPLX(XCL(12,I))
  453. C
  454. IF (NMAS.GT.0) THEN
  455. DO 55 I1 = 1 , NMAS
  456. IF (CORRES(MASS(1,I1)).EQ.NUMERO(I)) THEN
  457. MASS(4,I1) = IZA1-5
  458. END IF
  459. 55 CONTINUE
  460. END IF
  461. C
  462. IZA1 = IZA1 + 1
  463. C
  464. END IF
  465. C
  466. C--------------------- CONDITION AUX LIMITES SUR LE DEBIT
  467. C
  468. IF (FLAG((I-1)*17+15).EQ.0) THEN
  469. C
  470. DO 60 J = 1 , 2*NP
  471. IF (CORRES(J).EQ.NUMERO(I)) THEN
  472. C
  473. INP = INT((J+1)/2)
  474. CSEF = XPI * CARACT(4,INP) * CARACT(4,INP)
  475. CROF = CARACT(9,INP)
  476. N1 = 2*INP-1
  477. N2 = 2*INP
  478. X1 = COOR (1,N1)
  479. Y1 = COOR (2,N1)
  480. Z1 = COOR (3,N1)
  481. X2 = COOR (1,N2)
  482. Y2 = COOR (2,N2)
  483. Z2 = COOR (3,N2)
  484. XL = SQRT((X2-X1)**2 + (Y2-Y1)**2 + (Z2-Z1)**2)
  485. C
  486. XI1 = (X2-X1)/XL
  487. XI2 = (Y2-Y1)/XL
  488. XI3 = (Z2-Z1)/XL
  489. IF (J.EQ.N2) THEN
  490. XI1 = -XI1
  491. XI2 = -XI2
  492. XI3 = -XI3
  493. END IF
  494. C
  495. ZA1 (IZA1 , 14*(J-1)+14) = ZI
  496. IF (J.EQ.N2) ZA1 (IZA1 , 14*(J-1)+14) = -ZI
  497. ZA1 (IZA1 , 14*(J-1)+1) = -ZS * ( CROF*CSEF*XI1 )
  498. ZA1 (IZA1 , 14*(J-1)+2) = -ZS * ( CROF*CSEF*XI2 )
  499. ZA1 (IZA1 , 14*(J-1)+3) = -ZS * ( CROF*CSEF*XI3 )
  500. C
  501. END IF
  502. C
  503. 60 CONTINUE
  504. C
  505. ZSM (IZA1 ) = CMPLX(XCL(14,I))
  506. C
  507. IZA1 =IZA1 + 1
  508. END IF
  509. C
  510. C--------------------- CONDITION AUX LIMITES D'IMPEDANCE ACOUSTIQUE
  511. C
  512. IF (FLAG((I-1)*17+15).EQ.I) THEN
  513. C
  514. ISTOP = 0
  515. J = 0
  516. C
  517. 70 J = J + 1
  518. IF ((CORRES(J).EQ.NUMERO(I)).AND.(ISTOP.EQ.0)) THEN
  519.  
  520. C
  521. INP = INT((J+1)/2)
  522. CSEF = XPI * CARACT(4,INP) * CARACT(4,INP)
  523. CROF = CARACT(9,INP)
  524. N1 = 2*INP-1
  525. N2 = 2*INP
  526. X1 = COOR (1,N1)
  527. Y1 = COOR (2,N1)
  528. Z1 = COOR (3,N1)
  529. X2 = COOR (1,N2)
  530. Y2 = COOR (2,N2)
  531. Z2 = COOR (3,N2)
  532. XL = SQRT((X2-X1)**2 + (Y2-Y1)**2 + (Z2-Z1)**2)
  533. C
  534. XI1 = (X2-X1)/XL
  535. XI2 = (Y2-Y1)/XL
  536. XI3 = (Z2-Z1)/XL
  537. C
  538. ZA1 (IZA1 , 14*(J-1)+13 ) = CMPLX(XCL (15,I))
  539. C
  540. ZB = CMPLX(XCL (16,I))
  541. ZA1 (IZA1 , 14*(J-1)+14) = ZB
  542. ZA1 (IZA1 , 14*(J-1)+1) = -ZS * ( CROF*CSEF*XI1 )* ZB
  543. ZA1 (IZA1 , 14*(J-1)+2) = -ZS * ( CROF*CSEF*XI2 )* ZB
  544. ZA1 (IZA1 , 14*(J-1)+3) = -ZS * ( CROF*CSEF*XI3 )* ZB
  545. C
  546. ZSM ( IZA1 ) = CMPLX(XCL (17,I))
  547. C
  548. IZA1 = IZA1+1
  549. ISTOP = 1
  550. END IF
  551. C
  552. IF (J.NE.NP2) GO TO 70
  553. C
  554. END IF
  555. C
  556. 100 CONTINUE
  557. C
  558. RETURN
  559.  
  560. END
  561.  
  562.  
  563.  

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