* SPPLANC PROCEDUR OF166741 23/02/07 11:58:41 11590 ********************************************************************** * * * UNE PROCEDURE LE CALCUL DES SPECTRES DE PLANCHER * * * ********************************************************************** * ********************************************************************** ********************************************************************** * * * PROCEDURE S P P L A N C * * --------------------------- * * * * AUTEUR F. BEN JEMAA (DMT/SEMT/EMSI) * * * * CALCUL DES SPECTRES DE PLANCHER PAR UNE APPROCHE ANALYTIQUE * * THESE DE DOCTORAT DE 3 EME CYCLE 01/10/1990 - 30/03/1994 * * * ********************************************************************** * * * SPPLANC TAB1 * * * * TAB1 : OBJET DE TYPE TABLE * * * ********************************************************************** * * * TAB1 : objet de type table contenant * * * * * * Indice Type Commentaires * * * * STRUC TABLE Caractéristiques modales * * de la structure support * * * * Indice Type * * * * NMODE ENTIER Nombre de modes * * FREQU TABLE Les frequences modales * * AMORT TABLE Les amortissements modaux* * * * PLANCH TABLE Caractéristiques du * * Plancher étudié (P) * * * * Indice Type * * * * LISFREQ LISTREEL Axe fréquentiel du * * spectre de plancher * * PAR DEFAUT : * * discrétisation par la * * procédure "DISCRFR" * * * * COEFFPL TABLE Les coefficients de * * participation modale * * en P : ( Qn * PHIn ) / Mn* * * * AMORTPL FLOTTANT Amortissement du spectre * * de plancher à calculer * * * * TYPSPPL MOT Type du spectre de * * plancher * * 'DEPL' : déplacement * * 'VITE' : pseudo-vitesse * * 'ACCE' : pseudo-accél. * * * * EXCIT TABLE Données de l'excitation * * ( processus séparable ) * * * * Indice Type * * * * ENVE MOT Type de l'enveloppe * * 'PLATLIN' : * * montée-plat-descente * * ( par défaut : plateau ) * * * * DUREE FLOTTANT Durée du signal * * * * TDEB FLOTTANT Temps où commence * * le plateau ( montée ) * * * * TFIN FLOTTANT Temps où se termine * * le plateau ( début desc.)* * * * DSP EVOLUTION D.S.P. de la fonction * * aléatoire stationnaire * ********************************************************************** ********************************************************************** * *===================================================================== * LES DONNEES DU PROBLEME *===================================================================== * ZPI = PI ; Z2PI = 2.0 * ZPI ; * * ZNMOD = ZTAB. STRUC . NMODE ; TABF = ZTAB. STRUC . FREQU ; TEPS = ZTAB. STRUC . AMORT ; * TX = ZTAB. PLANCH . COEFFPL ; EPSPL = ZTAB. PLANCH . AMORTPL ; ZTYPSP = ZTAB. PLANCH . TYPSPPL ; * ZT = ZTAB. EXCIT . DUREE ; 'SI' ( 'EGA' ZENVE ZPTLINE ) ; ZTDEB = ZTAB. EXCIT . TDEB ; ZTFIN = ZTAB. EXCIT . TFIN ; 'SINON' ; 'MESS' 'L ENVELOPPE EST CONSIDERE PAR DEFAUT PLATEAU' ; ZENVE = ZPTPLAT ; 'FINSI' ; DSPE = ZTAB. EXCIT . DSP ; * *===================================================================== * CARACTERISATION DU MOUVEMENT DU PLANCHER * PAR UN PROCESSUS SEPARABLE *===================================================================== * 'MESS' ' ' ; 'MESS' 'CARACTERISATION DU MOUVEMENT DU PLANCHER' ; 'MESS' ' ' ; 'MESS' 'STRUCTURE SUPPORT : ' ; * * DEFINITION DE L'AXE DES TEMPS : * IMOD = 0 ; ZZ1 = 0.0 ; 'REPETER' BOUCM ZNMOD ; IMOD = IMOD '+' 1 ; ZZ0 = ( TABF.IMOD ) '*' ( TEPS.IMOD ) ; 'SI' ( ZZ0 '>' ZZ1 ) ; ZI0 = IMOD ; ZZ1 = ZZ0 ; 'FINSI' ; 'FIN' BOUCM ; * LTEV = DISCRTM TABF.ZI0 TEPS.ZI0 ZT ; LTEV0 = 'ENLEVER' LTEV 1 ; NTV1 = NTV '-' 1 ; * * CREATION DE L'ENVELOPPE D'EXCITATION * DANS LE CAS OU L'ENVELOPPE EST LINEAAIRE-PLAT-LINEAIRE * TABENV = 'TABLE' ; TABENV.1 = LTEV ; 'SI' ( 'EGA' ZENVE ZPTLINE ) ; ZDT1 = ZTDEB '/' 10.0 ; LENV1 = LT1 '*' ( 1.0 / ZTDEB ) ; ZDT2 = ( ZTFIN '-' ZTDEB ) '/' 10.0 ; ZDT3 = ( ZT '-' ZTFIN ) '/' 10.0 ; LENV3 = LENV3 '*' ( 1.0 '/' ( ZT '-' ZTFIN ) ) ; LT0 = LT1 'ET' LT2 'ET' LT3 ; LENV0 = LENV1 'ET' LENV2 'ET' LENV3 ; TITRE 'ENVELOPPE : TDEBUT =' ZTDEB 'TFIN =' ZTFIN ; DESSIN ENV 'CARR' 'MIMA' ; TABENV.2 = LENV ; 'FINSI' ; * * CALCUL DES COEFFICIENTS D'AJUSTEMENT * DE LA FONCTION DE TRANSFERT DE LA STRUCTURE SUPPORT * * APPEL DE LA PROCEDURE AJUSTE * IMOD = 0 ; LCT0 = 'TABLE' ; LEPST0 = 'TABLE' ; 'REPETER' BOUCM ZNMOD ; IMOD = IMOD '+' 1 ; * OMGI = TABF.IMOD * Z2PI ; LCT0.IMOD = COEFAJ.1 ; LEPST0.IMOD = COEFAJ.2 ; * DESSIN CT0 'CARR' 'MIMA' ; * LCT02 = ( LCT0.IMOD ) * ( LCT0.IMOD ) ; DESSIN CT02 'CARR' 'MIMA' ; * DESSIN EPST0 'CARR' 'MIMA' ; 'MESS' ' ' ; * FIN BOUCM ; * * ***** CALCUL DES FONCTIONS DE TRANSFERT MODALES APPROCHEES ***** * * * DISCRETISATION DE L'AXE DES FREQUENCE * FMIN = 0.0 ; FMAX = 100.0 ; LFREQ = DISCRFR TABF TEPS FMIN FMAX ; LFREQ2 = LFREQ * LFREQ ; * * INTERPOLATION DE LA DSP STATIONNAIRE * SUR LE NOUVEAU AXE DES FREQUENCES * XDSPE = 'EXTRAIRE' DSPE 'ABSC' ; * * ITEM = 0 ; TEM1 = 0.0 ; DSPAPP = EVINIT ; DSPAPP1 = EVINIT ; RVARMAXA = 0.0 ; REPETER BOUC2 NTV1 ; ITEM = ITEM + 1 ; * TEM2 = 'EXTRAIRE' LTEV0 ITEM ; DTEM = ( TEM2 - TEM1 ) / ( 2.0 * ZT ) ; * IMOD = 0 ; HA = ECUN0 ; REPETER BOUC3 ZNMOD ; IMOD = IMOD + 1 ; * LEPSITEM = LEPST0.IMOD ; EPSITEM = 'EXTRAIRE' LEPSITEM ITEM ; LCITEM = LCT0.IMOD ; CITEM = 'EXTRAIRE' LCITEM ITEM ; * FREQI = TABF.IMOD ; FREQI2 = FREQI * FREQI ; LRE = ( LUNF * FREQI2 ) - LFREQ2 ; * RIM = 2.0 * EPSITEM * FREQI ; LIM = LFREQ * RIM ; 'REEL' LRE 'IMAGIN' LIM ; HID = ( ECUN1 / EVCOMP ) * CITEM ; HIA = HID * FREQI2 * ( -1.0 ) ; HA = HA + ( HIA * ( TX.IMOD ) ) ; * FIN BOUC3 ; * ********* MODULE DE LA FONCTION DE TRANSFERT A L'INSTANT T0 ********** * * ******** D.S.P. DE LA REPONSE DE LA STRUCTURE A L'INSTANT T0 ********* * DSPAPP2 = MODHA * MODHA * DSPEC ; * ******* VARTIANCE DE LA REPONSE DE LA STRUCTURE A L'INSTANT T0 ******* * RVAR0 = 'EXTRAIRE' LSOM0 1 ; LVARAPP = LVARAPP 'ET' LSOM0 ; * * LA D.S.P. MOYENNE DE LA REPONSE DE LA STRUCTURE * DSPAPP = DSPAPP '+' ( ( DSPAPP1 '+' DSPAPP2 ) '*' DTEM ) ; * * MAXIMUM DE LA VARIANCE * D.S.P. A L'INSTANT OU LA VARIANCE EST MAXIMUM * LOGIC1 = ( RVAR0 > RVARMAXA ) ; LOGIC2 = (( 'ABS' (( RVAR0 - RVARMAXA ) / RVAR0 )) < 1.0E-2 ) ; SI ( OU LOGIC1 LOGIC2 ) ; TMAX = TEM2 ; RVARMAXA = RVAR0 ; DSPMAX = DSPAPP2 ; FINSI ; * TEM1 = TEM2 ; DSPAPP1 = DSPAPP2 ; * MENAGE ; FIN BOUC2 ; * TITRE 'VARIANCE : NMODE =' ZNMOD ; LVARAPP = 'INSERER' LVARAPP 1 0. ; DESSIN VARAPP 'CARR' 'MIMA' ; * * TITRE 'D.S.P. : NMODE =' ZNMOD ; LABSC = 'EXTRAIRE' DSPAPP 'ABSC' ; LABSC = 'ENLEVER' LABSC 1 ; DESSIN DSPAPP0 'LOGX' 'LOGY' 'CARR' 'MIMA' ; * VARMOY = VARMOY '/' ZT ; LENVPL = ( LVARAPP '/' VARMOY ) '**' 0.50 ; TABENVPL = 'TABLE' ; TABENVPL.1 = LTEV ; TABENVPL.2 = LENVPL ; * DSPPL = DSPAPP ; * *===================================================================== * CALCUL DU SPECTRE DE PLANCHER * AMORTISSEMENT : ZTAB. PLANCH . AMORTPL *===================================================================== * 'MESS' ' ' ; 'MESS' ' ' ; 'MESS' 'AMORTISSEMENT DU SPECTRE : ' EPSPL ; * * AXE DES FREQUENCES DU SPECTRE * 'SI' ( 'EXIST' ( ZTAB. PLANCH ) LISFREQ ) ; 'MESS' 'L AXE FREQUENTIEL DU SPECTRE EST DONNE PAR L UTILISATEUR' ; LFOSC = ZTAB. PLANCH . LISFREQ ; 'SINON' ; 'MESS' ' ' ; 'MESS' 'DISCRETISATION FREQUENTIELLE PAR LA PROCEDURE "DISCRFR" : [0.5HZ 30.0HZ]' ; FMINPL = 0.50 ; FMAXPL = 30.0 ; LFOSC = DISCRFR TABF TEPS FMINPL FMAXPL ; 'FINSI' ; * * IOSC = 0 ; TABFPL = TABF ; TEPSPL = TEPS ; ZNMOD1 = ZNMOD + 1 ; TEPSPL.ZNMOD1 = EPSPL ; 'REPETER' BOUC6 NOSC ; IOSC = IOSC + 1 ; * FREQI = 'EXTRAIRE' LFOSC IOSC ; TABFPL.ZNMOD1 = FREQI ; OMGI = Z2PI * FREQI ; * 'MESS' ' ' ; * * CALCUL DES COEFFICIENTS D'AJUSTEMENT * DE LA FONCTION DE TRANSFERT DE L'OSCILLATEUR AU PLANCHER * * APPEL DE LA PROCEDURE AJUSTE * LCI = COEFAJ.1 ; LEPSI = COEFAJ.2 ; * * DISCRETISATION DU DOMAINE DES FREQUENCES POUR L'INTEGRATION * LFREQ = DISCRFR TABFPL TEPSPL FMIN FMAX ; LFREQ2 = LFREQ '*' LFREQ ; * * INTERPOLATION DE DSPPL SUR LE NOUVEAU AXE FREQUENTIEL * XDSPPL = 'EXTRAIRE' DSPPL 'ABSC' ; * FREQI2 = FREQI '*' FREQI ; LFREQI2 = LUNF '*' FREQI2 ; LRE = LFREQI2 '-' LFREQ2 ; * ITEM = 0 ; TEM1 = 0.0 ; DSPAPP1 = EVINIT ; RVARMAXA = 0.0 ; 'REPETER' BOUC7 NTV1 ; ITEM = ITEM '+' 1 ; * TEM2 = 'EXTRAIRE' LTEV0 ITEM ; DTEM = ( TEM2 '-' TEM1 ) '/' ( 2.0 '*' ZT ) ; CI = 'EXTRAIRE' LCI ITEM ; * LIM = LFREQ '*' RIM ; 'REEL' LRE 'IMAGIN' LIM ; * HDEPL = ( ECUN1 / EVCOMP ) * ( CI / ( Z2PI * Z2PI ) ) ; * * MODULE DE LA FONCTION DE TRANSFERT A L'INSTANT T0 * * ******** D.S.P. DE LA REPONSE DE LA STRUCTURE A L'INSTANT T0 ********* * DSPAPP2 = MODHDEPL * MODHDEPL * DSPPL1 ; * ******* VARTIANCE DE LA REPONSE DE LA STRUCTURE A L'INSTANT T0 ******* * RVAR0 = 'EXTRAIRE' LSOM0 1 ; LVAROSC = LVAROSC 'ET' LSOM0 ; * * LA D.S.P. MOYENNE DE LA REPONSE DE LA STRUCTURE * DSPOSC = ( DSPAPP1 + DSPAPP2 ) * DTEM ; * * MAXIMUM DE LA VARIANCE * D.S.P. A L'INSTANT OU LA VARIANCE EST MAXIMUM * LOGIC1 = ( RVAR0 > RVARMAXA ) ; LOGIC2 = (( 'ABS' (( RVAR0 - RVARMAXA ) / RVAR0 )) < 1.0E-2 ) ; SI ( OU LOGIC1 LOGIC2 ) ; TMAX = TEM2 ; RVARMAXA = RVAR0 ; DSPMAX = DSPAPP2 ; FINSI ; * TEM1 = TEM2 ; DSPAPP1 = DSPAPP2 ; * MENAGE ; FIN BOUC7 ; * ******************** APPEL DE LA PROCEDURE LIVO ********************** ****** CALCUL DE N0, DELTA, TEMPS MODIFIE TMOD ET FACTEUR DE PIC ***** * LVAROSC = 'INSERER' LVAROSC 1 0. ; * * APPEL DE LA PROCEDURE LIVO * TPIC = LIVO VAROSC DSPOSC ; * * * * * FIN BOUC6 ; * TITRE 'SPECTRE N0 : NMODE =' ZNMOD 'AMORT. PL. =' EPSPL ; DESSIN ENBREAPP 'LOGX' 'LOGY' 'CARR' 'MIMA' ; * TITRE 'SPECTRE DELTA : NMODE =' ZNMOD 'AMORT. PL. =' EPSPL ; DESSIN EDELTAPP 'LOGX' 'LOGY' 'CARR' 'MIMA' ; * TITRE 'SPECTRE SIGMA : NMODE =' ZNMOD 'AMORT. PL. =' EPSPL ; DESSIN ESIGMAPP 'LOGX' 'LOGY' 'CARR' 'MIMA' ; * TITRE 'SPECTRE FACTPIC : NMODE =' ZNMOD 'AMORT. PL. =' EPSPL ; DESSIN EFPICAPP 'LOGX' 'LOGY' 'CARR' 'MIMA' ; * LSPPLDEP = LFPICAPP * LSIGMAPP ; * * TRANSFORMATION DU SRO DE PLANCHER CALCULE EN TYPE VOULU * TITRE 'SPECTRE PLANCHER : TYPE_SP =' ZTYPSP 'NMODE =' ZNMOD 'AMORT. PL. =' EPSPL ; 'SI' ( 'EGA' ZTYPSP ZDEPL ) ; ESPPLAPP = ESPPLDEP ; 'SINON' ; 'FINSI' ; DESSIN ESPPLAPP 'LOGX' 'LOGY' 'CARR' 'MIMA' ; * ZTABAPP = 'TABLE' ; ZTABAPP. VARP = VARAPP ; ZTABAPP. DSPP = DSPAPP ; ZTABAPP. NBR0 = ENBREAPP ; ZTABAPP. DELT = EDELTAPP ; ZTABAPP. FPIC = EFPICAPP ; ZTABAPP. SPPL = ESPPLAPP ; 'FINPROC' ZTABAPP ;
© Cast3M 2003 - Tous droits réservés.
Mentions légales