ML18029A932

Flora:Program for Direct Generation of Floor Response Spectra,Users Manual,Version 3A.

ML18029A932
Person / Time
Site:	San Onofre
Issue date:	04/15/1984
From:	SOUTHERN CALIFORNIA EDISON CO.
To:
Shared Package
ML13324A582	List: ML18029A931 ML18029A932 ML18029A933 ML18029A934
References
PROC-840415, NUDOCS 8504020158
Download: ML18029A932 (40)

Text

FLORA Secondary Spectra Generation From Known Primary Spectra User's Manual Version: 3A April 15, 1984 Approved By: Date:

s >sson Nanager 850 402 850329 AD C< 0 500020'i58 PDR L

P

l DEFINITIONS of VARIABLES REFERENCED in the MANUAL VARIABLE DESCRIPTION REFERENCED NST*RT Program restart option. Use with the narr,ow-band option to save time I

NPSD The number of output power spectral density functions to be calculated.

PROPS Variable representing a file name assigned by the user. PROPS refers to the file in which are stored the modal properties of the structure.

MAXI Maximum number of iterations to be used when computing the input spectra's power spectral density function.

SPCNAM Variables representing file names PSDNAM assigned by the user. The files OUTFILE are created by FLORA dur ing execution and contain the output spectra~ the output power spectr a and a complete execution listing.

NDG Number of input spectra dampings.

NM Max number of input spectrum points.

KWT Input response spectra scale type

( fr equency/per iod) switch.

IFSC Output spectra scale type switch.

NMP Number of power spectral density def inition points.

OMAX The frequency above which no er ror checking is performed when iterating on the power spectrum. The PSD typi cal ly conver ges faster at its lo~er fr equen cy points an d s 1 ower as the frequency incr eases. Execut ion goes up disproportionately with increasing OMAX.

(vi)

FMAX The highest frequency at which the PSD is nonzer o. It should be equal to or greater than the highest signif i cant input spectrum frequency.

It must be higher than OMAX~ a good value is 15% higher.

The PSD frequen cy sca l e.

ZMAX The PSD gener at ion a 1gor i thm requir es that the PSD ordinate goes to zero at a 'high frequency. This value should be much higher than ZM*X.

NWDG The number of frequency points in the referenced input response spectrum.

NDP The number of damping ratios for which a response spectrum will be generated at a single output degree of freedom.

Ply 1

P

FLORA Version 3A User' Manual NOTE: Two files are. required to run FLORA. The first file contains this input data and must reside in the user's file space and be named 'FLORIN'. The second file contains the structure's modal properties and output response location information. It is described in Section II of this manual. If the wide band option is chosen FLORA creates two output files, one containing the entire execution'output listing and one containing all generated output spectra. If the narrow band option is used the same two are created plus a.-third containing the input and all output PSD's.

I. HEADING LINE (l 9AO)

Columns 1-76 : TITLE Any desired heading II. EXECUTION MODE AND RESTART OPTION (4I5, 3x, 7al, 3x, 7al, 3x, 7al, 3x, 7al)

Columns 1-5 IM Execution Mode EQ.O Narrow Band Option EQ.l Wide Band Option 6-10 NSTART Restart Cycle (Only for Narrow Band)

EQ.O Normal Execution EQ.l Compute PSDF using MAXI iterations and stop EQ.2 Iterate MAXI times on previously computed PSDF and stop EQ.3 Iterate (or not if MAXI = 0) on previously computed PSDF and complete response calculations

~ i 0

11-15 MAXI (Only for Narrow Band) Maximum allowable number of iterations to be performed in computing the PSDF (or reiterating on the PSDF for NSTART.GT.l) 16-20 NPSD (Only for Narrow Band) Number of output frequency points for which an output power spectral density function is calculated.

Output PSDs use the same frequency scale as the input PSD. See Appendix I, "Notes About Input" 24-30 SPCNAM Seven-character filename in which are stored the frequency, period, acceleration and standard deviation on acceleration for each output response spectrum (4E16.8) 34-40 PSDNAM (Only for Narrow Band) Seven-character filename in which are stored the frequency and PSD for the input PSD and each output PSD in order of their calculation (2E16.8) 44-50 OUTFILE Seven-character filename in which the entire output listing is stored 54-60 PROPS Seven-character filename from which the formatted primary structure modal data are to be read NOTES: 1) For NSTART = 2, input only Sections I and II.

For NSTART = 3, input Sections I, II, IX and X.

2) The four file names must be left justified and may be up to seven characters in length.

3) PROPS contains the primary system modal data stored as given in this pseudo code:

OPEN (4,FILE= 'PRIMARY'ORM='FORMATTED',

1 ACCESS='SEQUENTIAL')

WRITE(4,1000) nf, neqpt WRITE(4,1001)(px(i,l), px(i,2), px(i,3), i=1, nf)

WRITE(4,1002)(freq(i), i=1, nf)

~

Pg

~

0

00 100 i=1, nf 100 WRITE(4, 1002) (xk( i, j), j= 1, neqpt) 1000 FORMAT(2I5) 1001 FORMAT(3E16.8) 1002 FORMAT(4E16.8)

CLOSE(4)

ENO where nf number of modes used (must be

'greater than one) neqpt number. of attachment degrees of freedom (DOFs) for which secondary spectra are sought in the FLORA run reading these properties px(nf,3) = x, y II z mass-normalized mass participation factors for each mode, T

Mr P" t freq(nf) = frequencies in CPS

" xk(nf,neqpt) = nf mass normalized mode shape elements for each mode at each of ttteneqpt attachment degrees of freedom

'If the structure's modal properties are calculated using Impell's program EDSGAP, the reformatting program GAP must be run before FLORA. GAP will create the file PROPS. See the GAP user's manual as Appendix II to this manual.

I

• A word of caution. None of the elements of the xk array can be zero. A zero will cause a divide by zero error. Before running the program change any zero elements to some very small number relative to the other elements. If the element is small it will have a negligible effect on the results.

3

e III. PRIMARY SPECTRUM INPUT CONTROL (5I5, F10.3)

Co lumns 1-5 NDG Number of dampings for which the input spectrum is specified. Leave blank if the "built-in" NRC Regulatory Guide 1.60 spectra are used.

I 6-10 IDIR Direction of input primary spectra application*

EQ.l - X Direction EQ.2 - Y Direction EQ.3 - Z Direction

• One input excitation direction per execution of FLORA.

11-15; NM Maximum number of frequency points in the input spectrum. Leave blank if RG-1.60 spectra are used 16-20 KWT Input spectrum scale code NE.l Frequency vs. Acceleration EQ.1 Period vs. Acceleration 21-25 KHV RG-1.60 Spectrum Direction NE.l Horizontal Spectrum EQ.l Vertical Spectrum 26-35 ZPA Zero period acceleration for RG-1.60 spectra IV. SECONDARY SPECTRA COMPUTATION CONTROL (2I5$ 3F10.4)

Columns 1-5 NFP Number of frequency points at which secondary spectral ordinates are to be generated 4

~p

~ J ~

~

e 6-10 IFSC Scale Generation:

LT.1 Spectral frequency points are generated evenly distributed on a logarithmic scale Eg.l 'Spectral frequency points are generated evenly distributed on a linear scale GT.1 Frequency points are specified by the user 11-20 Lowest Frequency (CPS) 21-30 W2 Highest Frequency (CPS) 31-40 FRAg Exceedance expectation of response. .

If zero, the mean of the peak response is calculated. See Appendix I, "Notes About Input" for more detailed information V. POWER SPECTRAL DENSITY FUNCTION CONTROL (I5$ 6F10.40 I5)

For the wide-band option insert one blank line. For the narrow-band option input the following:

Co 1 umns 1-5 NMP Number of frequency points for which the PSDF ordinates are generated EQ.O, NMP = 50 6-15 OMAX Maximum frequency'pon which the error check is performed EQ.O, OMAX = 33 cps 16-25 FMAX Maximum frequency in the generated PSDF for which the PSD ordinate is nonzero Eg.0, FMAX = 40 cps

E ~

4

26-35 ZMAX Zero PSDF Frequency (CPS)

EQ.O, ZMAX = 10.* FMAX 36-45 XCEE Damping ratio for which the PSDF is computed EQ.O, XCEE = 0.02 46-45 TOL Allowable error in percent used to check convergence of PSDF points at frequencies less than OMAX EQ.O, TOL ~ 5%

56-65 TAU Duration in seconds of original input time history used in generating the input acceleration response spectra EQ.O, TAU = 10. seconds 66-70 IIN Flag for printing the generated PSDF NE.O Print EQ.O Do Not Print NOTE: FLORA assumes the PSDF values are specified at NMP+2 frequency points, Wi, =. 0, 1, ...NMP+1, with PSDF =

i 0 at WO = 0.0 and at WNMP+1 = ZMAX. FLORA generates NMP PSDF values at the frequencies W;, i =

1, 2, ...., NMP evenly distributed between WO = 0.0 and WNMP = FMAX on a linear scale. Tolerance checks are made only at points with frequencies below OMAX.

VI. OUTPUT FREQUENCY POINTS AT WHICH OUTPUT PSDs ARE CALCULATED (SI10)

Enter NPSD values. If NPSD of Section II is zero insert one blank line.

1-10 First output response frequency point with respect to which an output PSD is calculated. See Appendix I, "Notes About Input" 11-20 Second point, etc

0 VII. MODAL DAMPING FACTORS OF PRIMARY STRUCTURE Input as many cards as required to specify the structure's modal damping factors in the format (8F10.0).

Note: Each mode has its own damping ratio. The admissible range for damping ratio is between 0.0 (no damping) and 1.0 (100K of critical damping). If the modal damping ratio for a particular mode is left blank, then the modal damping factor of the previous mode is used for the analysis. Thus, if the modal damping factors are the same for all modes, only the first modal damping value need be specified, with appropriate number of blanks.

VIII. PRIMARY INPUT SPECTRA If NDG'Eg.O in the primary spectrum input control card (Section III), skip this section. Otherwise, input the following cards; (1 A5i - 0 L; Columns 1-80 HED Heading for the spectral table to be printed with output Repeat the following set of lines for each input spectra (NDG set of lines).

S ectrum Control Line (F10.0, I10, F10.0) - One Line Columns 1-10 DGJ Damping factor for the input spectrum 11-20 NWGJ Number of frequency points defining this spectrum. This number must not exceed NM specified in Section III 21-30 SFJ Scale factor for the input spectral values. Default to 1.0

~SC d (2E16.8)

As many lines as required to specify NWGJ frequency (or period) - acceleration pairs, one pair to a line. Lines must be consistent with scale code KWT specified in Section III.

~y The first number of each pair is the frequency (or period) and the second the corresponding spectral acceleration. Lines must be in ascending order of frequency (or period).

Fre uenc S ecification. Lines (8F10.0)

For IFSC = 0 or 1, ignore this section. If IFSC (Columns 6-10,Section IV) is greater than 1, provide as many lines as required to specify NFP (Columns 1-5,Section IV) frequencies (in cps) according to the format (8F10.0).

IX. SECONDARY SPECTRA GENERATION CONTROL (3I5, F10.4)

Input one line as defined here followed by as many lines as required by Section X for each location and degree of freedom for which secondary spectra will be generated. This information must be consistent with the PROPS file information of Section II.

Columns 1-5  : NODE Structure node number where secondary spectra is to be generated 6-10 IDOF Degree of freedom of above node num-ber for secondary spectra generation*

EQ.l X-Direction EQ.2 Y-Direction EQ.3 Z-Direction

• Both in-line and cross-direction spectra may be generated using either the wide-band or narrow-band option.

11-15 NDP Number of damping values for which the secondary spectra are to be generated at the above node number and degree of freedom 16-25 EN Attached Equipment mass. If the mass of the equipment attached to the location of interest is included, equipment-structure interaction will be considered.

Inputting a zero mass will neglect this interaction and result in conventional floor spectra

4~

C'

X. SPECTRAL DAMPING LINE (F10.0)

As many lines as required to specify NDP (Columns 11-15,Section IX) damping values expressed'as ratio to critical damping, according to the format (F10.0).

XI. DATA FILE INPUT/OUTPUT In all cases, the input data file (containing this input information) is the file named 'FLORIN'.

A) EDSGAP SUPPLIED MODAL PROPERTIES INPUT: Tapes 3, 7, 8 and 9 of the EDSGAP "NSTART = 2" run are used to define the modal properties of the primary structure. The program 'GAP'ust be run to tran-scribe these tapes into FLORA's format. See Appendix II for the GAP User's Manual.

B) OTHER SOURCES OF MODAL PROPERTIES INPUT: The notes to the file PROPS of Section II may be followed to prepare an input data file where the input data file, 'FLORIN', and the primary structure property file, PROPS, exist in the user's work file area.

OUTPUT: The output is stored as defined by SPCNAM, OUTFILE, and PSDNAM, described in Section II.

For runs with NSTART greater than 0, and for most error termination when IM = 0, the input data along with the intermediate calculations are stored in a file created by FLORA and called 'RESTART'.

XII. CYBER SYSTEM FLORA ACCESS FLORA can be accessed and executed on the IMPELL in-house computer as follows:

GET FLORA/UN=IMPLIB.

LL FLORA.

J l

XIII. VAX A SYSTEN ACCESS An executable copy of FLORA can be copied into any local file space by typing:

COPY "USERl:)MILLER]FLORA.EXE"

• followed by a carriage'eturn. Assuming 'FLORIN'nd the PROPS file exist in the local file space, run FLORA by:

RUN FLORA Be sure to er ase all FLORA related files when finished, disk space is always critically low.

I APPENDIX I: NOTES ABOUT INPUT Section Variable Notes II IM 1) The narrow band option calculates a power spectral density function consistent with the input response spectra, and uses this PSDF to -compute each output response spectra.

2) The wide band option assumes that the PSDF of the input response spectra is constant over the entire frequency range, and uses approximate equations for output response.

NSTART Used only for the narrow band option. For runs with NSTART greater than zero, and for most error induced terminations, the control information and the PSD are written to files

'RESTART'nd PSDNAM respectively.

MAXI A range of 10 to 25 should be sufficient for most applications.

NPSD Since FLORA perturbs the structure's modal properties differently at each output frequency point to account for interaction, different output PSDs will result when calculated with respect to different output frequencies. If the attached equipment mass is small, the PSD differences will also be small.

NDG One spectra for each damping value.

IV FRAQ The specification of FRAg allows a complete probabilistic description of the response.

Specifying a value greater than zero but less than 1.0 will produce an output spectrum where each ordinate can be expected to be exceeded FRAg X 100% of the time. That is, the spectrum will represent responses for that uniform level of exceedance.

~ ~t

~

(

k

e Section Variable Notes If FRAg is 0.0, FLORA will produce a response spectrum corresponding to other "non-probabilistic" spectra generation methods. A deterministic, time history method produces such a spectrum. Each ordin'ate will represent the average peak response with an unknown exceedance level.

Different ordinates will have different exceedance probabilities. The probability density .function of the peak response depends on many parameters, typically resulting in the mean not necessarily being the 50X exceedance level.

IV Must be greater than zero for the logarithmic scale.

NMP The output response spectra is sensitive to

- the number of PSDF points. The default value of 50 is recomended as a starting value.

Use the nstart = 1 option to select the number of PSD points giving a relatively smooth PSD in the frequency range of interest.

OMAX The PSDF is calculated iteratively and for a given PSDF, the error tends to increase with the frequency of the PSDF point. ONX defines the highest frequency for which TOL must be satisfied.

Though the error of the PSDF for a frequency higher than OMAX is not checked, gust having the higher frequency terms there improves the accuracy of the PSDF at lower frequencies.

ZNX Use the default value. This is the frequency where the PSD ordinate returns to zero. This frequency should be very far from the frequency range of interest so as not to introduce numerical errors.

XCEE Theoretically, if all the input response spectra are derived from a single time history, the choice of XCEE, which is the damping ratio for .which the PSDF is computed,

Section Variable Notes is arbitrary. Since the input response spectra are sometimes artificial it's best to use a damping the same as the average modal damping of the structure.

TOL This is the percent error you will tolerate for any PSDF point equal to or below OMAX.

TAV This is the duration in seconds of the stationary portion of the original excitation time-history., The random vibration based formulation requires this. Since we rarely know the Eg duration use the default value

- Al.3-

a f U

,

'Appendix II Nodal Data Reformatting Program GAP Coupling EDSGAP to FLORA

- A2.1-

Interface Program for fDSGAP to FLORA User's Manual Version: lA April 1, 1984 Approved by: Date: /

>vssson anager /

.- A2.2-

0 J

GAP - User's Manual Version Log Version Date Prepared by Approved by Sumary of Changes

<Pa < (3~~) ~4/(e~ r, na 5f crW

- A2.3-

Introduction GAP reformats structural modal properties generated by EDSGAP into a format acceptable by the program FLORA. The tapes created by a NSTART 2 run of EDSGAP necessary for a FLORA run are "TAPE3",

"TAPET " "TAPES" and "TAPES". Infonsatson in these $ s read t

an d h en wr itt t en to a file with a name referenced by tables e variable PROPS. i To successfully run GAP, those four tapes must be indirect access files in the user's file space with actual filenames, "TAPE3," etc. A file with a name referenced by PROPS will be created and saved during the GAP execution. The data described here must be entered into a file named "GAPIN". This file must reside in the user's file space. An echo print of the input data is written to the file "GAPOUT"; this file is created by the FORTRAN77.

First Line {A7, 2X, 2I5)

Columns Variable Entry 1 - 7 PROPS Name of file to contain results generated by this program. This name also needs to be entered in columns 54-60 on the second line of FLORA's input data file.

s 10 - 15 NF Number of modes to use in the FLORA analysis. Less than or equal to the number generated by EDSGAP but must be greater than one.

16 - 20 NEgPT Number of output response degrees of freedom; i.e. the number of spectra generation locations.

Second Line (2I5) *(HEgPT lines required) 1 -2 NODE Node number 6 - 10 IDOF direction of output spectrum calculation

• Use same order asSection IX if results at multiple DOF's are to be calculated in a single run.

- A2.4-

{i -I A

~,