* @RCCMCO2 PROCEDUR CHAT 03/09/23 21:15:03 4697 NSEG*ENTIER SIG1*MCHAML MOD1*MMODEL EPAI*FLOTTANT * ---------------------------------------------------------------------- * * CETTE PROCEDURE PERMET D'EFFECTUER UNE COUPE DANS UN MODELE 2D OU 3D * POUR OBTENIR LES CONTRAINTES ELEMENTAIRES A CHAQUE POINT ET DANS * LE REPERE LOCAL DE CELLE-CI. * CETTE PROCEDURE PERMET EGALEMENT LE TRACE DE CHAQUE CONTRAINTE * ELEMENTAIRE ET DE SA DECOMPOSITION LE LONG DE CETTE COUPE PAR * APPEL A LA PROCEDURE DECOMP * * ---------------------------------------------------------------------- * * LES DONNEES SONT : * ----> le numero d'etat * ----> le numero de la coupe * ----> les 2 points extremes de la coupe * ----> le nombre de sous-divisions de la coupe * ----> le champ de contraintes aux elements (points de GAUSS) * ----> la geometrie * ----> l'epaisseur du trait de la coupe lors du trace * ----> le Logique GRAPL permettant ou non les traces * ----> le mot NOM (optionnel) permettant ou non la sortie d'un * fichier pour NOCLIB * ----> Pour le 3D PP3 DIRECTION (type point) definissant * l axe OY du repere ( OX etant le segment d appui ) en 2D * il ne servira a rien ( mais obligatoire tout de meme * EN SORTIE : * * ----> le champ de contraintes dans la coupe * ----> les contraintes elementaires x distances a chaque point * par rapport au milieu de la coupe * ----> les contraintes 'TOT..' aux extremites de la coupe * ----> les evolutions la decomposition de chaque contrainte elementaire * la coupe si GRAPL est VRAI * * ----> eventuellement un fichier pour NOCLIB * ( voir NOTA ci-dessous ) * * %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% * * NOTA : Les contraintes du fichier NOCLIB sont dans le repere General * * * %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% * ---------------------------------------------------------------------- * PAR Jean-Marc PIETTE * --------------------------------------------------------------------- 'OPTION' 'ELEM' 'SEG2' ; * 'SI' ( 'EGA' IDIM 2 ) ; 'OPTION' 'ELEM' 'QUA4' ; 'SI' GRAPL ; PTD = PTC 'TOURNER' -90. PP1 ; PTE = PTC 'TOURNER' 90. PP1 ; MAILCOU = L1M L2M L3M L4M DALLER 'PLAN' ; MAILCOU = MAILCOU 'COULEUR' 'ROSE' ; 'FINSI' ; SINON ; * ---------------------------------------------------------------------- * Le calcul est donc tridimensionnel * ---------------------------------------------------------------------- OPTION ELEMENT CUB8 ; V1 = ( PP2 MOINS PP1 ) ; 'SI' GRAPL ; V2 = PP3 ; VH1 = VH1 * EPAI / ( NORME VH1 ) ; MAILCOU = VO1 ET VO2 ; MAILCOU = MAILCOU COULEUR ROSE ; 'FINSI' ; CHAMRT = RTENS SIG1 MOD1 V1 V2 ; FINSI ; *----------------------------------------------------------------------- * projection du champ de contraintes sur le segment d appui *----------------------------------------------------------------------- 'SI' ('EGA' ( NSEG+1) NPRO ) ; SI ECRI ; LIST CHCOU ; FINSI ; * ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ * eventuel fichier pour NOCLIB *~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ SI ( NEG NOM 'BLANC' ) ; OPTION 'IMPR' NOM ; MESSAGE ' COORDONNEES GLOBALES ----> NB DE POINTS = ' NBNDC ; REPETER BOUCLN NBNDC ; I = &BOUCLN ; VECTD = P MOINS PP1 ; LONGD = NORME VECTD ; *~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ * ecriture des coordonnees *~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ MESSAGE LONGD ' 0. 0. ' ; FIN BOUCLN ; * ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ * ecriture des contraintes elementaires pour NOCLIB *~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ TTT1 = CHAIN ' COUPE No ' NNC ' ' ; repeter bib NCC; fin bib ; MESSAGE TTT1 ; REPETER BOUCLC ( NBNDC ) ; I = &BOUCLC ; repeter bib (NCC - 1); fin bib ; MESSAGE TTT2 ; FIN BOUCLC ; 'OPTION' 'IMPR' VA_IMP ; FINSI ; 'SI' GRAPL ; * TITRE 'Contraintes segment appui no' nnc 'Etat' nne ; TTTT = TABLE ; TABLEG = TABLE ; TTTT.'TITRE' = TABLEG ; 'REPETER' BIB NCC ; TABLEG.&BIB = CCC ; 'SI' ( 'EGA' &BIB 1) ; 'SINON' ; EVS = EVS ET 'FINSI' ; 'FIN' BIB ; * 'DESS' EVS 'LEGE' TTTT AXES; 'FINSI' ; *---------------------------------------------------------------------- * decomposition en membrane flexion etc --- des contraintes *---------------------------------------------------------------------- 'SINON' ; MESS ' POINTS DU SEGMENT D APPUI ' NNC ' HORS DU DOMAINE ' ; TABV= 'ARRET' ; 'FINSI' ; FINPROC TABV ;
© Cast3M 2003 - Tous droits réservés.
Mentions légales