C IDMAT4 SOURCE CHAT 07/02/16 21:15:00 5668 SUBROUTINE IDMAT4 (NUMP1,NUMP2,NUDIR1,NUDIR2,ANG,MELEME,MINTE, . IPVAL,NPG2,MINTE2) IMPLICIT INTEGER(I-N) IMPLICIT REAL*8(A-H,O-Z) * ************************************************************************ * * I D M A T 4 * ----------- * * FONCTION: * --------- * * CALCUL DES COS-DIR. DES AXES ORTH. (OU ANISO.)PAR RAPPORT AU * * REPERE LOCAL D'ELEMENT , A PARTIR D'UNE DEFINITION GLOBALE * * * MODULES UTILISES: * ----------------- * -INC SMCOORD -INC CCOPTIO -INC CCREEL -INC SMELEME -INC SMINTE * * PARAMETRES: (E)=ENTREE (S)=SORTIE (+ = CONTENU DANS UN COMMUN) * ----------- * * NUMP1 (E) NUMERO DU POINT P1 ASSOCIE A NUDIR1 * NUMP2 (E) NUMERO DU POINT P2 ASSOCIE A NUDIR1 * NUDIR1 (E) NUMERO DE DIRECTIVE UTILISEE DANS LA LISTE: * "DIRECTION", "RADIAL" * NUDIR2 (E) NUMERO DE DIRECTIVE UTILISEE DANS LA LISTE: * "PARALLELE", "PERPENDIC.", "INCLINE", * POUR LA DEFINITION DES DIRECTIONS D'ORTHOTROPIE. * ANG (E) ANGLE UTILISE DANS LA DEFINITION DES DIRECTIONS * D ORTHOTROPIE (INCLINE) * MELEME (E) POINTEUR DE "MAILLAGE" A 1 SEUL TYPE D'ELEMENT. * MINTE (E) SEGMENT CONTENANT LES FONCTIONS DE FORME,ET LEURS * DERIVEES * XVAL (S) DIRECTIONS D'ORTHOTROPIE PAR ELEMENT, * ON FOURNIT LES COSINUS DIRECTEURS DE L'AXE 1 * EN BIDIM ,ET LES COSINUS DIRECTEURS DES AXES 1 ET * 2 EN TRIDIM,PAR RAPPORT AU REPERE LOCAL DE L'ELEMENT * * * VARIABLES: *----------- * VGLOB1,VGLOB2,VGLOB3 = COS.DIRECTEURS DES AXES 1, 2 ET 3 D'ORTH. PAR * RAPPOT AU REPERE GLOBAL * TXR = TABLEAU CONTENANT LES COS.DIRECTEURS DES AXES LOCAUX PAR * RAPPORT AU REPERE GLOBAL * * * REMARQUES: * ---------- * * LES DIRECTION P1 ET P2 DEFINISSENT LE PLAN QUI CONTIENT LES AXES 1,2 * NUDIR2 =1 SIGNIFIE QUE L'AXE 1 EST PARALLELE A LA DIRECTION P1, * NUDIR2 =2 SIGNIFIE QUE L'AXE 1 EST PERPENDICULAIRE A LA DIRECTION P1, * NUDIR2 =3 SIGNIFIE QUE L'AXE 1 FAIT AN ANGLE DE (ANG) AVEC LA * DIRECTION P1 * * * AUTEUR, DATE DE CREATION: * ------------------------- * * P. DOWLATYARI OCT. 1990 * * LANGAGE: * -------- * FORTRAN 77 + ESOPE * * ************************************************************************ * * DIMENSION VGLOB1(3),VGLOB2(3),VGLOB3(3),VAUX1(3),VAUX2(3) DIMENSION VPT1(3),VPT2(3) * SEGMENT YVAL REAL*8 VLOC1(IDIM,NPG2,NBELEM),VLOC2(IDIM,NPG2,NBELEM) ENDSEGMENT * SEGMENT WORK REAL*8 XE(3,NBNN),TXR(IDIM,IDIM) ENDSEGMENT * SEGACT MELEME NBELEM=NUM(/2) NBNN=NUM(/1) SEGACT MINTE IF(NUDIR1.EQ.2) SEGACT MINTE2 * * DIRECTIONS FOURNIES * IF (IDIM.EQ.2)THEN * * EXTRAIRE LES COORDONNEES DU POINT NUMP1 * CALL EXCOO1 (NUMP1,VPT1(1),VPT1(2),VPT1(3),REEL1) * IF(NUDIR1.EQ.2) THEN VAUX1(1)=VPT1(1) VAUX1(2)=VPT1(2) GO TO 15 ENDIF * * NORMALISATION * VNORM=DSQRT(VPT1(1)*VPT1(1)+VPT1(2)*VPT1(2)) IF (VNORM.EQ.0.) THEN CALL ERREUR (524) SEGDES,MINTE,MELEME RETURN ENDIF VGLOB1(1)=VPT1(1)/VNORM VGLOB1(2)=VPT1(2)/VNORM VGLOB2(1)=-VGLOB1(2) VGLOB2(2)=VGLOB1(1) * IF(NUDIR2.EQ.2)THEN * * ON EFFECTUE UNE ROTATION DE 90 DEGRE AUTOUR DE L'AXE 3 * VGLOB1(1)=VGLOB2(1) VGLOB1(2)=VGLOB2(2) * ELSEIF(NUDIR2.EQ.3)THEN * * ON EFFECTUE UNE ROTATION DE (ANG) DEGRE AUTOUR DE L'AXE 3 * COSA=DCOS(ANG) SINA=DSIN(ANG) VGLOB1(1)=VGLOB1(1)*COSA+VGLOB2(1)*SINA VGLOB1(2)=VGLOB1(2)*COSA+VGLOB2(2)*SINA ENDIF * ELSEIF(IDIM.EQ.3)THEN * * COORDONNEES DU POINT P1 * CALL EXCOO1 (NUMP1,VPT1(1),VPT1(2),VPT1(3),REEL1) * IF(NUDIR1.NE.2) THEN * * NORMALISATION * VNORM=DSQRT(VPT1(1)*VPT1(1)+VPT1(2)*VPT1(2)+ . VPT1(3)*VPT1(3)) IF (VNORM.EQ.0.) THEN CALL ERREUR (524) SEGDES,MINTE,MELEME RETURN ENDIF VGLOB1(1)=VPT1(1)/VNORM VGLOB1(2)=VPT1(2)/VNORM VGLOB1(3)=VPT1(3)/VNORM ENDIF * * COORDONNEES DU POINT P2 * CALL EXCOO1 (NUMP2,VPT2(1),VPT2(2),VPT2(3),REEL1) * IF(NUDIR1.EQ.2) THEN VGLOB1(1)=VPT2(1)-VPT1(1) VGLOB1(2)=VPT2(2)-VPT1(2) VGLOB1(3)=VPT2(3)-VPT1(3) * * NORMALISATION * VNORM=DSQRT(VGLOB1(1)*VGLOB1(1)+VGLOB1(2)*VGLOB1(2)+ . VGLOB1(3)*VGLOB1(3)) VGLOB1(1)=VGLOB1(1)/VNORM VGLOB1(2)=VGLOB1(2)/VNORM VGLOB1(3)=VGLOB1(3)/VNORM GO TO 15 ENDIF * * PRODUIT VECTORIEL DES DEUX VECTEURS ET NORMALISATION * CALL CROSS2 (VGLOB1(1),VPT2(1),VGLOB3(1),IRR) IF (IRR.EQ.0)THEN * * DIRECTIONS D'ORTHOTROPIE MAL FOURNIES * CALL ERREUR (524) SEGDES,MINTE,MELEME IF(NUDIR1.EQ.2) SEGDES MINTE2 RETURN ENDIF CALL CROSS2 (VGLOB3(1),VGLOB1(1),VGLOB2(1),IRR) * IF(NUDIR2.EQ.2)THEN * * ON EFFECTUE UNE ROTATION DE 90 DEGRE AUTOUR DE L'AXE 3 * REEL1=VGLOB1(1) VGLOB1(1)=VGLOB2(1) VGLOB2(1)=-REEL1 REEL1=VGLOB1(2) VGLOB1(2)=VGLOB2(2) VGLOB2(2)=-REEL1 REEL1=VGLOB1(3) VGLOB1(3)=VGLOB2(3) VGLOB2(3)=-REEL1 * ELSEIF(NUDIR2.EQ.3)THEN * * ON EFFECTUE UNE ROTATION DE (ANG) DEGRE AUTOUR DE L'AXE 3 * COSA=DCOS(ANG) SINA=DSIN(ANG) REEL1=VGLOB1(1) VGLOB1(1)=VGLOB1(1)*COSA+VGLOB2(1)*SINA VGLOB2(1)=-REEL1*SINA+VGLOB2(1)*COSA REEL1=VGLOB1(2) VGLOB1(2)=VGLOB1(2)*COSA+VGLOB2(2)*SINA VGLOB2(2)=-REEL1*SINA+VGLOB2(2)*COSA REEL1=VGLOB1(3) VGLOB1(3)=VGLOB1(3)*COSA+VGLOB2(3)*SINA VGLOB2(3)=-REEL1*SINA+VGLOB2(3)*COSA ENDIF ENDIF * 15 CONTINUE * SEGINI YVAL IPVAL=YVAL SEGINI WORK * * BOUCLE SUR LES ELEMENTS * DO 10 IEL=1,NBELEM * DO 20 NC=1,IDIM DO 20 NV=1,NPG2 VLOC1(NC,NV,IEL)=0.D0 VLOC2(NC,NV,IEL)=0.D0 20 CONTINUE * * COORDONNEES DES NOEUDS * CALL DOXE(XCOOR,IDIM,NBNN,NUM,IEL,XE) * * CALCUL DU REPERE LOCAL DE L'ELEMENT * NBSH=SHPTOT(/2) CALL RLOCAL(XE,SHPTOT,NBSH,NBNN,TXR) if (nbsh.eq.-1) then call erreur(525) return endif IF(IERR.NE.0)THEN SEGSUP,WORK,YVAL SEGDES , MELEME,MINTE IF(NUDIR1.EQ.2) SEGDES MINTE2 RETURN ENDIF * * ON BOUCLE SUR LES POINTS DE GAUSS * DO 80 IGAU=1,NPG2 IF (NUDIR1.EQ.2) THEN DO 25 ICO=1,IDIM VAUX2(ICO)=0.D0 DO 25 ID=1,NBNN VAUX2(ICO)=VAUX2(ICO)+ . MINTE2.SHPTOT(1,ID,IGAU)*XE(ICO,ID) 25 CONTINUE * IF(IDIM.EQ.2) THEN VGLOB1(1)=VAUX1(1)-VAUX2(1) VGLOB1(2)=VAUX1(2)-VAUX2(2) * * NORMALISATION * VNORM=DSQRT(VGLOB1(1)*VGLOB1(1)+VGLOB1(2)*VGLOB1(2)) VGLOB1(1)=VGLOB1(1)/VNORM VGLOB1(2)=VGLOB1(2)/VNORM VGLOB2(1)=-VGLOB1(2) VGLOB2(2)=VGLOB1(1) * IF(NUDIR2.EQ.2)THEN * * ON EFFECTUE UNE ROTATION DE 90 DEGRE AUTOUR DE L'AXE 3 * VGLOB1(1)=VGLOB2(1) VGLOB1(2)=VGLOB2(2) * ELSEIF(NUDIR2.EQ.3)THEN * * ON EFFECTUE UNE ROTATION DE (ANG) DEGRE AUTOUR DE L'AXE 3 * COSA=DCOS(ANG) SINA=DSIN(ANG) VGLOB1(1)=VGLOB1(1)*COSA+VGLOB2(1)*SINA VGLOB1(2)=VGLOB1(2)*COSA+VGLOB2(2)*SINA ENDIF ENDIF * IF(IDIM.EQ.3) THEN VAUX1(1)=VAUX2(1)-VPT2(1) VAUX1(2)=VAUX2(2)-VPT2(2) VAUX1(3)=VAUX2(3)-VPT2(3) PSC=VAUX1(1)*VGLOB1(1)+VAUX1(2)*VGLOB1(2)+VAUX1(3)*VGLOB1(3) VGLOB2(1)=VAUX1(1)-PSC*VGLOB1(1) VGLOB2(2)=VAUX1(2)-PSC*VGLOB1(2) VGLOB2(3)=VAUX1(3)-PSC*VGLOB1(3) * NORMALISATION VNORM=DSQRT(VGLOB2(1)*VGLOB2(1)+VGLOB2(2)*VGLOB2(2)+ . VGLOB2(3)*VGLOB2(3)) VGLOB2(1)=VGLOB2(1)/VNORM VGLOB2(2)=VGLOB2(2)/VNORM VGLOB2(3)=VGLOB2(3)/VNORM * * PRODUIT VECTORIEL DES DEUX VECTEURS ET NORMALISATION * CALL CROSS2 (VGLOB1(1),VGLOB2(1),VGLOB3(1),IRR) * IF (IRR.EQ.0)THEN CALL ERREUR (524) SEGDES,MINTE,MELEME IF(NUDIR1.EQ.2) SEGDES MINTE2 RETURN ENDIF * IF(NUDIR2.EQ.2)THEN * * ON EFFECTUE UNE ROTATION DE 90 DEGRE AUTOUR DE L'AXE 3 * REEL1=VGLOB1(1) VGLOB1(1)=VGLOB2(1) VGLOB2(1)=-REEL1 REEL1=VGLOB1(2) VGLOB1(2)=VGLOB2(2) VGLOB2(2)=-REEL1 REEL1=VGLOB1(3) VGLOB1(3)=VGLOB2(3) VGLOB2(3)=-REEL1 * ELSEIF(NUDIR2.EQ.3)THEN * * ON EFFECTUE UNE ROTATION DE (ANG) DEGRE AUTOUR DE L'AXE 3 * COSA=DCOS(ANG) SINA=DSIN(ANG) REEL1=VGLOB1(1) VGLOB1(1)=VGLOB1(1)*COSA+VGLOB2(1)*SINA VGLOB2(1)=-REEL1*SINA+VGLOB2(1)*COSA REEL1=VGLOB1(2) VGLOB1(2)=VGLOB1(2)*COSA+VGLOB2(2)*SINA VGLOB2(2)=-REEL1*SINA+VGLOB2(2)*COSA REEL1=VGLOB1(3) VGLOB1(3)=VGLOB1(3)*COSA+VGLOB2(3)*SINA VGLOB2(3)=-REEL1*SINA+VGLOB2(3)*COSA ENDIF ENDIF ENDIF * * CALCUL DES COS.DIRECTEURS DU REPERE LOCAL / REPERE GLOBAL * IF(IDIM.EQ.2)THEN VLOC1(1,IGAU,IEL)=TXR(1,1)*VGLOB1(1)+TXR(2,1)*VGLOB1(2) VLOC1(2,IGAU,IEL)=TXR(1,2)*VGLOB1(1)+TXR(2,2)*VGLOB1(2) * ELSEIF(IDIM.EQ.3)THEN DO 40 J=1,3 DO 40 I=1,3 VLOC1(J,IGAU,IEL)=VLOC1(J,IGAU,IEL)+VGLOB1(I)*TXR(I,J) VLOC2(J,IGAU,IEL)=VLOC2(J,IGAU,IEL)+VGLOB2(I)*TXR(I,J) 40 CONTINUE ENDIF * 80 CONTINUE 10 CONTINUE * * DESACTIVATION DES SEGMENTS * SEGSUP,WORK SEGDES , MELEME,MINTE,YVAL IF(NUDIR1.EQ.2) SEGDES MINTE2 * RETURN END