C CFLUO3D SOURCE OF166741 25/11/04 21:15:19 12349 c A.Sellier nov 2023 SUBROUTINE cfluo3d(WRK52,WRK53,WRK54,MWRKXE,WR14, # nbnb,idimb,teta1b,teta2b,nvarib,nstrsb,ifourb,dtb,trefb) C WRK52,53,54 segments déclarés dans le common DECHE C XMAT(NCOMAT) = COMPOSANTES DE MATERIAU C IVAL(NCOMAT) = INDICE DES COMPOSANTES DE MATERIAU C NCOMAT = NOMBRE DE COMPOSANTES DE MATERIAU C XCAR(ICARA) = CARACTERISTIQUES C MFR = NUMERO DE LA FORMULATION DE L'ELEMENT FINI C = 1 MASSIF C = 3 COQUE MINCE ( COQ2 , COQ3 ET DKT ) C = 5 COQUE EPAISSE ( COQ6 , COQ8 ) C = 7 POUTRE C = 9 COQUE MINCE AVEC CISAILLEMENT TRANSVERSE ( COQ4 ET DST ) C = 11 LIQUIDE C = 13 TUYAU C = 15 LINESPRING C = 17 TUYAU FISSURE C = 19 RACCORD MASSIF C = 21 RACCORD COQUE C = 23 SURFACE LIBRE C = 25 MEMBRANE C = 27 UNIAXIALE C = 29 THERMIQUE C = 31 INCOMPRESSIBLES C = 33 POREUX C = 35 JOINT C = 37 HOMOGENEISE C = 39 TUYO C = 41 TUYAU ACOUSTIQUE PURE C = 43 RACCORD TUYAU FLUIDE c* DDAUX = MATRICE DE HOOKE ELASTIQUE c* NSTRS = NBRE DE COMPOSANTES DES DEFORMATIONS c* CMATE = NOM DU MATERIAU c* VALMAT= TABLEAU DE CARACTERISTIQUES DU MATERIAU c* VALCAR= TABLEAU DE CARACTERISTIQUES GEOMETRIQUES c* N2EL = NBRE D ELEMENTS DANS SEGMENT DE HOOKE c* N2PTEL= NBRE DE POINTS DANS SEGMENT DE HOOKE c* MFR = NUMERO DE LA FORMULATION c* IFOU = type de formulation c* IB = NUMERO DE L ELEMENT COURANT c* IGAU = NUMERO DU POINT COURANT c* EPAIST= EPAISSEUR c* NBPGAU= NBRE DE POINTS DE GAUSS c* MELE = NUMERO DE L ELEMENT FINI c* NPINT = NBRE DE POINTS D INTEGRATION c* NBGMAT= NBRE DE POINTS DANS SEGMENT DE CARACTERISTIQUES c* NELMAT= NBRE D ELEMENTS DANS SEGMENT DE CARACTERISTIQUES c* SECT = SECTION c* LHOOK = TAILLE DE LA MATRICE DE HOOKE c* TXR,XLOC,XGLOB,D1HOOK,ROTHOO,DDHOMU,CRIGI = TABLEAUX UTILISES c* POUR LE CALCUL DE LA MATRICE DE HOOKE c* matrice de Hoocke DDHOOK(LHOOK,LHOOK) C----------------------------------------------------------------------- C VARIABLES PASSEES PAR LES COMMONS COPTIO , ECOU ET NECOU C C IFOUR INDICE DU TYPE DE PROBLEME C -3 DEFORMATIONS PLANES GENERALISEES C -2 CONTRAINTES PLANES C -1 DEFORMATIONS PLANES C 0 AXISYMETRIQUE C 1 SERIE DE FOURIER C 2 TRIDIMENSIONNEL C ITYP TYPE DE FORMULATION MECANIQUE C --------------- ATTENTION --------------- C IL EST ACTIF APRES L APPEL DE VISAVI C ----------------------------------------- C ITYP=1 CAS DES ELEMENTS MASSIFS C ITYP=2 CAS DES COQUES C ITYP=3 CAS DES MEMBRANES C ITYP=4 CAS DES CABLES ET DES BARRES C ITYP=5 CAS QUELCONQUE C ITYP=6 CAS DES CONTRAINTES PLANES C ITYP=7 CAS DES COQUES A NU=0. OU CONTRAINTES PLANES C ITYP=8 CAS DES MEMBRANES A NU=0. OU CONTRAINTES PLANES C ITYP=9 CAS DES COQUES EPAISSES C ITYP=10 CAS DES JOINTS C ITYP=11 CAS DES POUTRES C ITYP=12 CAS DES TUYAUX C ITYP=13 CAS DES COQUES AVEC CISAILLEMENT TRANSVERSE C C ISTEP flag utilise pour separer les etapes dans un calcul non local C ISTEP=0 -----> calcul local C ISTEP=1 -----> calcul non local etape 1 on calcule les seuils C ISTEP=2 -----> calcul non local etape 2 on continue le calcul C a partir des seuils moyennes C ISTEP=3 ----->sous iteration de Helmholtz C ISTEP est défini par DEFINI.ESO C C EPIN0 déformations inélastiques initiales C EPINF déformations inélastiques finales C epst0 déformation totale initiale C EPSTF déformation totale finale C C----------------------------------------------------------------------- C SORTIES C SIGF(NSTRS) = CONTRAINTES FINALES C VARF(NVARI) = VARIABLES INTERNES FINALES C DEFP = DEFORMATIONS PLASTIQUES C KERRE = 0 TOUT OK C----------------------------------------------------------------------- IMPLICIT INTEGER(I-N) IMPLICIT REAL*8(A-H,O-Z) c ajout sellier 2022 -INC PPARAM c fin ajout sellier 2022 -INC CCOPTIO -INC DECHE c segment de coordonnees des noeuds de l element integer nbnb,insb,idimb SEGMENT MWRKXE REAL*8 XE(3,nbnb) ENDSEGMENT * sellier/millard segment pour les Helmholtz 04 04 2020 integer NBVIA SEGMENT WR14 INTEGER INLVIA(NBVIA) ENDSEGMENT c----- nombre de noeuds max et table locale des coordonnees ------------ -INC HNBRNEU c rappel du contenu de HNBRNEU C tableau local des coordonnees des noeuds de lEF C integer NBNMAX3D C nbr de noeuds maxi par element C parameter (NBNMAX3D=20) C tableau des coordonnees des noeuds C real*8 XE3D(3,NBNMAX3D) c----------------------------------------------------------------------- c-------------- Parametres pour le non locale -------------------------- -INC HNBRHEL c NB_HELM: Nombre de variable d etat principales c NB_PARA_PAR_HELM: Nombre de parametres par variable d etat c NB_VARI_PAR_HELM: Nombre de variable internes par variables d etat c NB_PARA_HELM: nombre de parametre total resultant c NB_VARI_HELM: nombre de variables internes totales resultantes c INLVIA3D(NB_HELM):copie du pointeur sur les numeros des vari traitee par helmholtz c NBVIA3D: recuperation du nombre de Helmholtz actives c declaration des tableaux pour helmholtz -INC HDECHEL c----------------Parametres pour les renforts -------------------------- c-INC HNBRREN c NB_RENF c NB_PARA_PAR_RENF c NB_PARA_RENF c NB_VARI_PAR_RENF c NB_VARI_RENF c------------- Nombre de parametres pour le modele de fibre ------------ c-INC HNBRF3D c pas de fibre pour l instant dans FLUISO3D c----------------------------------------------------------------------- c-- transfert des coordonnees des noeuds ------------------------------- c recuperation de la dimension hors plan en cas de calcul 2D c boucle sur les noeuds if(nbnb.le.NBNMAX3D) then NBNB3D=nbnb if(idimb.le.3) then idimb3d=idimb else print*,'pb table de coord des noeuds ds cfluo3d' kerre=1 return end if else print*,'Element avec + de noeuds que capacite declaree dans' print*,'HNBRNEU.INC' kerre=1 return end if c chargement des coordonnees des noeuds do insb=1,NBNB3D do j=1,idimb3d c print*,'xe(',j,insb,')=',xe(j,insb) XE3D(j,insb)=xe(j,insb) end do if (idimb.lt.3) then do j=idimb+1,3 XE3D(j,insb)=0.d0 end do end if end do c on double les neouds si le probleme est 2D if ((ifourb.eq.0).or.(ifourb.eq.-1)) then if(mfr.ne.33) then c position de la dimension en non poreux dimension3=XMAT3D(16) else c position en poreux print*,'dans cfluo3d ajouter dimension3' kerre=1 return end if NBNB3D=2*nbnb do insb=nbnb+1,NBNB3D do j=1,2 c print*,'xe(',j,insb,')=',xe(j,insb) XE3D(j,insb)=XE3D(j,(insb-nbnb)) end do XE3D(3,insb)=dimension3 end do end if c------------- test du nombre de variables de Helmholtz ---------------- c transfert dans un tableau de fluendo3d dimensionne au max des c facultative helmholtz de IDMATR (actuellement 10) if(NBVIA.GT.NB_HELM) then print*,'Pb de dimensionnement avec les variables de Helmholtz' print*,'dans cflui3d:',NBVIA,' > ',NB_HELM print*,'Verifier NB_HELM et nbr de non locales' print*,'declarees nombre_helmholtz.h (NBRHEL.INC)' else c transfert vers les variables locales if (NBVIA.gt.0) then NBVIA3D=NBVIA do i=1,NB_HELM if(i.le.NBVIA3D) then INLVIA3D(i)=INLVIA(i) else INLVIA3D(i)=0 end if end do else NBVIA3D=NBVIA do i=1,NB_HELM INLVIA3D(i)=0 end do end if end if c------------ fin de recuperation des numeros des variables Helmholtz -- c------------- appel au modele ---------------------------------------- call FLUO3D(XMAT,NMAT,SIG0,SIGF,DEPST,NSTRSb,VAR0,VARF, # NVARIb,TETA1b,TETA2b,DTb,KERRE,MFR,IFOURB,ISTEP,EPST0, # EPSTF,TREFB,NBVIA3d,INLVIA3d,NBNMAX3D,NBNB3D,IDIMB3D,XE3D) return end c***********************************************************************