1 : C LISPRI SOURCE CHAT 05/01/13 01:23:38 5004 2 : C RAIDEUR LINESPRING 3 : SUBROUTINE LISPRI(XE,CAR,D,WORK,MELE,REL,I70,I195,I157,I158) 4 : C======================================================================= 5 : C 6 : C CALCULE LA MATRICE DE RIGIDITE DU LINE SPRING 7 : C 8 : C ENTREES 9 : C XE (3,4)=COORDONNES DE LA POUTRE LINE SPRING 10 : C EPA1 =EPAISSEUR NOEUDS 1 ,4 CAR(1) 11 : C FISS1=PROFONDEUR NOEUDS 1 ,4 CAR(2) 12 : C V1(3)=VECTEUR ORIENTANT LES NOEUDS 1 ,4 CAR(3,4,5) 13 : C EPA2 =EPAISSEUR NOEUDS 1 ,4 CAR(6) 14 : C FISS2=PROFONDEUR NOEUDS 1 ,4 CAR(7) 15 : C V2(3)=VECTEUR ORIENTANT LES NOEUDS 1 ,4 CAR(8,9,10) 16 : C D(2,2) =MATRICE DE HOOKE 17 : C WORK(90)=TABLEAU DE TRAVAIL 18 : C MELE = 30 OU 50 NUMERO DE L ELEMENT LINESPRING 19 : C SORTIES 20 : C REL(24,24)=MATRICE DE RIGIDITE 21 : C I70,I195,I157,I158 MESSAGE D ERREURS 22 : C 23 : C CODE EBERSOLT OCTOBRE 84 24 : C 25 : C======================================================================= 26 : IMPLICIT INTEGER(I-N) 27 : IMPLICIT REAL*8(A-H,O-Z) 28 : DIMENSION XE(3,*),CAR(*),D(2,*),WORK(*),REL(24,*) 29 : DIMENSION V1(3),V2(3) 30 : EPA1 =CAR(1) 31 : FIS1 =CAR(2) 32 : V1(1)=CAR(3) 33 : V1(2)=CAR(4) 34 : V1(3)=CAR(5) 35 : EPA2 =CAR(6) 36 : FIS2 =CAR(7) 37 : V2(1)=CAR(8) 38 : V2(2)=CAR(9) 39 : V2(3)=CAR(10) 40 : C 41 : IF(MELE.EQ.30) THEN 41 : CALL LISPR2(XE,EPA1,FIS1,V1,EPA2,FIS2,V2,D, 42 : 1 WORK(1),WORK(10),REL,I70,I195,I157,I158) 43 : ENDIF 44 : C 45 : IF(MELE.EQ.50) THEN 45 : CALL LISPR3(XE,EPA1,FIS1,V1,EPA2,FIS2,V2,D, 46 : 1 WORK(1),WORK(10),REL,I70,I195,I157,I158) 47 : ENDIF 48 : C 49 : RETURN 50 : END 51 :
© Cast3M 2003 - All rights reserved.
Disclaimer