Vibtoolbox is a set of command line tools for analysis of vibration and acoustic test data. The interface design used for most of the tools in the package is the standard UNIX filter model. The tools take input on standard input and transform the input in some fashion, and then send the transformed input data to standard output. The tool set also defines ASCII and native binary database formats for storage of vibration and acoustic laboratory and flight test data.
Vibtoolbox was started in 2003 by Greg Rudd, out of the need to have a set of tools to handle vibration and acoustic test data in a general and portable way. An important goal of the project was that the code for the analysis tools would be self contained and easy to port and maintain on different UNIX operating systems.
The programs in this package are designed to work with SI absolute units, US absolute units, and US gravitational units. The programs do not work with the US engineering system of units that require a additional factor to be included in the relation between force and mass. Any consistent units may be used in the programs. This means that all parameters must be expressed in the same units within the program. Examples of consistent units that work in these programs are:
------------------------------------------------------ System Force Length Mass Time ------------------------------------------------------ SI kilonewton meter kilogram second SI newton millimeter kilogram second US pound-force inches lbf-s^2/in second US pound-force foot slug second US poundal foot pound second
The vibtoolbox utilities are used like other standard UNIX command line utilities. Standard output from one program can be passed to standard input of an another program via UNIX pipes. Also the programs process input and output data using UNIX I/O redirection. The programs in the package use a subset of the gnuplot data format for the STDIN and STDOUT streams. The utilities handle commend line options using the getopt utility.
For example lets see how we would use the filters to analysis some flight test time history data.
Assume for the example we have the flight data stored in the vibtoolbox ASCII database format.
The ASCII data can then be packed to native binary format using:
packdb flights channels < DATABASE > DATABASE.BIN Where: flights is a list of flight conditions to process channels is a list of channels to process DATABASE is the input ASCII database DATABASE.BIN is the output binary database
Note flight conditions and channels are input using a colon operator to
for example: 1,2,10:15,6 would expand to 1,2,10,11,12,13,14,15,6
To read a time record from the binary database use:
breaddb conditions channels < DATABASE.BIN > TimeHistory Were TimeHistory contains the extracted time history in gnuplot format.
Calculate basic properties of the time history record using:
checkth < timehistory Sample output is: #BLOCK Flight = 23 Channel = 12104 Samples per Second = 5998.784 Standard Deviation Sample Rate = 0.167 Maximum signal level = 266.659 Minimum signal level = -266.667 Signal duration = 60.012 Number of samples = 360000 Signal standard deviation = 5.932e+01 Signal mean = 2.257e-01 Crest factor = 4.496e+00
Calculate autospectral density of the time history using:
asd PointsPerSegment NumberOfAverages < TimeHistory > AutoSpectralDensity
Calculate autospectral density without writing out the time history to a file using:
cat DATABASE.BIN | breaddb conditions channels | asd PointsPerSegment NumberOfAverages > AutoSpectralDensity
The vibtoolbox tools can be used with other UNIX utility's on the command line. For example to add calculation of RMS signal level to the above calculation use:
cat DATABASE.BIN | breaddb conditions channels | asd PointsPerSegment NumberOfAverages | tee ASD.OUT | rms > RMS.OUT autospectral density is stored in file ASD.OUT using the UNIX tee command.
This program calculates random response spectra given autospectral density input. The response is calculated referenced to inertial space coordinates. USAGE: absrs [OPTIONS] STDIN STDOUT OPTIONS: -h prints help message -q number = quality factor default quality factor = 10 STDIN: The input file is a free-format ASCII text file with two columns. Frequency (Hz) and autospectral density (eu^2/hz) eu is engineering unit. There may be more then one data set in the input file. Each data set is separated by 2 empty lines. STDOUT: The output file is a free-format ASCII text file with two columns. Frequency (hz) and RMS Response Amplitude (eu) eu is engineering unit. There may be more then one data set in the output file. Each data set is separated by 2 empty lines. This file can be processed by gnuplot.
This program calculates airspeed parameters and atmosphere properties. USAGE: airspeed [OPTIONS] STDIN STDOUT OPTIONS: -h prints help message STDIN: To calculate airspeeds given Mach use: "Mach" altitude(ft) Mach To calculate airspeeds given keas use: "keas" altitude(ft) keas To calculate airspeeds given ktas use: "ktas" altitude(ft) ktas To calculate airspeeds given kcas use: "kcas" altitude(ft) kcas To calculate properties at altitude use: "prop" altitude(ft) STDOUT: Given altitude in feet and one of the airspeeds(Mach, KEAS, KCAS, KTAS), the program will return. Altitude(ft) keas kcas ktas Mach Dynamic Pressure(psf) Also given altitude in feet the program can return the following atmosphere properties. Alt(ft) Density(slug/ft^3) Temp(degree R) Press(psf) Speed of Sound(fps) viscosity(slug/ft/sec)
This program calculates autospectral density given input time history. Usage: asd [OPTIONS] "Points Per Segment" averages STDIN STDOUT OPTIONS -m = remove mean from data -a = calculate analysis parameters given the total number of analysis points and the sample rate. Returns a table of points per segment, number of segments to average, frequency resolution, and standard error. -h = print help message Points Per Segment = number points per analysis segment averages = number of segments to average STDIN: The input file is a free-format ASCII text file with two columns: Time (sec) and amplitude (eu) eu is engineering unit. There may be more then one data set in the input file. Each data set in separated by 2 empty lines. STDOUT: The output file is a free-format ASCII text file with two columns: Frequency (hz) and autospectra density (eu^2/hz) eu is engineering unit. There may be more then one data set in the output file. Each data set in separated by 2 empty lines. This file can be processed by gnuplot.
This program transforms acceleration PSD(eu^2/Hz) to velocity PSD(eu^2/Hz). eu = engineering units USAGE: asdvsd [OPTION] STDIN STDOUT OPTIONS -h prints help message -g = factor to normalize acceleration by gravity default = 386.1 STDIN: The input file is a free-format ASCII text file with two columns: Frequency(Hz) and PSD(eu^2/Hz) EU is engineering unit. There may be more then one data set in the input file each data set in separated by 2 blank lines STDOUT: The output file is a free-format ASCII text file with two columns: Frequency(Hz) and PSD(eu^2/Hz) eu is engineering unit. There may be more then one data set in the output file. Each data set in separated by 2 blank lines. This file can be processed by gnuplot.
This program calculates autocorrelation of signal. USAGE: autocorr [OPTION] STDIN STDOUT OPTIONS: -h prints help message -l number of lag points STDIN: The input file is a free-format ASCII text file with two columns: Time and amplitude(EU) EU is engineering unit. There may be more then one data set in the input file each data set in separated by 2 blank lines STDOUT: The output file is a free-format ASCII text file with two columns: lag Rxx There may be more then one data set in the output file. Each data set in separated by 2 blank lines. This file can be processed by gnuplot.
Calculates average of a set of vibration spectra. USAGE: avespec [OPTIONS] STDIN STDOUT OPTIONS: -h prints help message STDIN: The input file is a free-format ASCII text file with 2 columns. Frequency (hz) and amplitude spectra (eu) or autospectral density (eu^2/hz) There may be more then one data set in the input file each data set in separated by 2 blank lines STDOUT: The output file is a free-format ASCII text file with two columns: Frequency (hz) and amplitude spectra (eu) or autospectral density (eu^2/hz) This file can be processed by gnuplot.
Reads time history from binary database and transforms to autospectral density. USAGE: bdb2psd [OPTIONS] CONDITIONS CHANNELS POINTS AVERAGES STDIN STDOUT OPTIONS -s = number = segment input signal into power of 2 size segments number = segment size (power of 2) -m = remove mean from time signal. -a = calculate analysis parameters given the total number of analysis points and the sample rate. returns a table of points per segment, number of segments to average, frequency resolution, and standard error. -h print help message CONDITIONS = list of conditions to parse from database CHANNELS = list of channels to parse from database Note conditions and channels are input using a colon operator to for example: 1,2,10:15,6 would expand to 1,2,10,11,12,13,14,15,6 POINTS = number of points in analysis segment AVERAGES = number of autospectral densities to average STDIN: database file binary format. GOOD/BAD scale x0 dx np condition channel y ... GOOD/BAD = mark data as good or bad scale = scale factor to apply to data x0 = starting time(sec) dx = delta time(sec) np = number of amplitudes condition = condition(integer number) channel = channel(integer number) y = amplitude(eu) eu is engineering unit. STDOUT: The output file is a free-format ASCII text file with two columns: Frequency (hz) and autospectral density (eu^2/hz) eu is engineering unit. There may be more then one data set in the output file. Each data set is separated by 2 blank lines. This file can be processed by gnuplot.
This program calculates natural frequencies of beams. USAGE: beamfn STDOUT Interactive input STDOUT: The output file is a free-format ASCII text file: Natural Frequency(Hz) = fn
Calculates conditions and channels that define spectra maximum envelope from binary database. USAGE: bmaxfc CONDITIONS CHANNELS STDIN STDOUT CONDITIONS = list of conditions to parse from database CHANNELS = list of channel numbers to parse from database Note conditions and channels are input using a colon operator to for example: 1,2,10:15,6 would expand to 1,2,10,11,12,13,14,15,6 STDIN: The database file is a native binary format. GOOD/BAD scale x0 dx np condition channel y ... GOOD/BAD = mark data as good or bad scale = scale factor to apply to data x0 = starting frequency dx = delta frequency np = number of amplitudes condition = condition(integer number) channel = channel(integer number) y = amplitude(eu) or autospectral density (eu^2/hz) eu is engineering unit. STDOUT: The output file is a free-format ASCII text file with three columns: Frequency(hz) and envelope condition number and envelope channel number There may be more then one data set in the output file. Each data set is separated by 2 blank lines. This file can be processed by gnuplot. STOP
Read data from binary database to create setup file for cross spectral density or coherence analysis programs csd and cohere. USAGE: breadcd [OPTIONS] FILE STDIN STDOUT OPTIONS: -h = help message FILE is a file with three columns that defines flights and channels to calculate in format: Condition and Forced Channel and Response Channel STDIN: The database file in a native binary format. GOOD/BAD scale x0 dx np condition channel y ... GOOD/BAD = mark data as good or bad scale = scale factor to apply to data x0 = starting time(sec) dx = delta time(sec) np = number of time points condition = condition(integer number) channel = channel(integer number) y = amplitude(eu) eu is engineering unit. STDOUT: The output file is a free-format ASCII text file with three columns: time (sec) and amplitude (eu) and amplitude (eu) eu is engineering unit. There may be more then one data set in the output file. Each data set is separated by 2 blank lines. This file can be processed by gnuplot.
Read time history or spectra from binary database. USAGE: breaddb [OPTIONS] CONDITIONS CHANNELS STDIN STDOUT OPTIONS: -h = help message CONDITIONS = list of flight conditions to parse from database CHANNELS = list of channel numbers to parse from database Note conditions and channels are input using a colon operator to for example: 1,2,10:15,6 would expand to 1,2,10,11,12,13,14,15,6 STDIN: The database file in native binary format. GOOD/BAD scale x0 dx np condition channel y ... GOOD/BAD = mark data as good or bad scale = scale factor to apply to data x0 = starting time(sec) or frequency dx = delta time(sec) or frequency np = number of amplitudes condition = condition(integer number) channel = channel(integer number) y = amplitude(eu) or spectra amplitude STDOUT: The output file is a free-format ASCII text file with two columns: time(sec) or frequency and amplitude(eu) or spectra amplitude (eu) or autospectral density (eu^2/hz) eu is engineering unit. There may be more then one data set in the output file. Each data set is separated by 2 blank lines. This file can be processed by gnuplot.
Calculate breakpoints on log-log coordinates. USAGE: breakpoint [OPTIONS] STDIN STDOUT OPTIONS: -h prints help message STDIN: The input file is a free-format ASCII text file with two columns: Frequency (hz) and amplitude (eu) or autospectral density (eu^2/hz) eu is engineering unit. There may be more then one data set in the input file. Each data set is separated by 2 empty lines STDOUT: The output file is a free-format ASCII text file with two columns: Frequency (hz) and amplitude (eu) or autospectral density (eu^2/hz) eu is engineering unit. There may be more then one data set in the output file. Each data set in separated by 2 empty lines. This file can be processed by gnuplot.
This program calculates damping of detected peaks in amplitude spectra or autospectral density. USAGE: calq [OPTIONS] STDIN STDOUT OPTIONS: -h prints help message STDIN: The input file is a free-format ASCII text file with two columns: Frequency (hz) and amplitude (eu) or autospectral density (eu^2/hz) eu is engineering unit. There may be more then one data set in the input file each data set in separated by 2 blank lines STDOUT: The output file is a free-format ASCII text file with two columns: Frequency(hz) and Quality factor There may be more then one data set in the output file. Each data set in separated by 2 blank lines. This file can be processed by gnuplot.
Differentiation of time history data. USAGE: cdiff [OPTION] SCALE STDIN STDOUT SCALE = scale factor for calculated result. Use 1/386.1 to scale acceleration in inches/seconds^2 calculated from velocity input time history in inches/seconds to units of g's. OPTIONS: -h prints help message STDIN: The input file is a free-format ASCII text file with two columns: Time(seconds) and Amplitude(EU) EU is engineering unit. There may be more then one data set in the input file each data set is separated by 2 empty lines Output File: The output file is a free-format ASCII text file with two columns: Time(seconds) and Amplitude(EU) EU is engineering unit. There may be more then one data set in the output file. Each data set in separated by 2 empty lines. This file can be processed by gnuplot.
This program calculates fast Fourier transform given complex time history input. The input and output data are in the same engineering units. Usage: cfft [OPTIONS] STDIN STDOUT OPTIONS: -i = real and imag output -h = print help message STDIN: The input file is a free-format ASCII text file with three columns: Frequency(hz) and Amplitude (eu) and Amplitude (eu) There may be more then one data set in the input file. Each data set in separated by 2 empty lines. STDOUT: The output file is a free-format ASCII text file with three columns: Frequency(hz) and [ amplitude (eu) and phase (radian) ] or [ real and imag ] There may be more then one data set in the output file. Each data set in separated by 2 empty lines. This file can be processed by gnuplot.
This program checks database array size matches array size parameter. USAGE: checkdb STDIN STDOUT STDIN: The database file is a free-format ASCII text file. GOOD/BAD scale x0 dx np flight channel y ... GOOD/BAD = mark data as good or bad scale = scale factor to apply to data x0 = starting time(sec) or frequency dx = delta time(sec) or frequency np = number of amplitudes flight = flight condition(integer number) channel = accelerometer channel(integer number) y = amplitude STDOUT: The output file is a ASSCI file. number point parameter = number number of points in array = number or for bad lines number point parameter = number number of points in array = number error
This program calculates Samples per Second, Standard Deviation Sample Rate, Number of samples, Signal duration, Maximum signal level, Minimum signal level, Signal mean, Signal standard deviation, Signal skewness, Signal Kurtosis, Signal excess Kurtosis, and Crest factor. This program also checks if the time history is ascending. USAGE: checkth [OPTION] STDIN STDOUT OPTION -h prints help message STDIN: The input file is a free-format ASCII text file with two columns: time(sec) and amplitude There may be more then one data set in the input file each data set is separated by 2 empty lines STDOUT: The outfile file is a free-format ASCII text file: Samples per Second Standard Deviation Sample Rate Number of samples Signal duration Maximum signal level Minimum signal level Signal mean Signal standard deviation Signal skewness Signal Kurtosis Signal excess Kurtosis Crest factor There may be more then one data set in the output file. Each data set is separated by 2 empty lines.
This program clips the amplitudes of a time history at a maximum level. USAGE: clipth [OPTIONS] STDIN STDOUT OPTIONS: -h prints help message -f number = data clipping level STDIN: The input file is a free-format ASCII text file with two columns: time (sec) and amplitude (eu) There may be more then one data set in the input file. Each data set is separated by 2 empty lines. STDOUT: The output file is a free-format ASCII text file with two columns: time (sec) and amplitude (eu) There may be more then one data set in the output file. Each data set in separated by 2 empty lines. This file can be processed by gnuplot.
This program calculates ordinary coherence for two input channels. Usage: cohere [OPTIONS] "Points per Segment" "Number of Averages" STDIN STDOUT OPTIONS: -m = remove mean from data -a = calculate analysis parameters -h = print help message Points per Segment = number points per analysis segment Number of Averages = number of segments to average STDIN: The input file is a free-format ASCII text file with three columns. The input file can be created with program breadcd: Frequency(hz), amplitude channel 1, amplitude channel 2 There may be more then one data set in the input file. Each data set in separated by 2 empty lines. STDOUT: The output file is a free-format ASCII text file with two columns: Frequency(hz) and coherence There may be more then one data set in the output file. Each data set in separated by 2 empty lines. This file can be processed by gnuplot.
This program calculates cross spectral density for two input channels. Usage: csd [OPTIONS] "Points per Segment" "Number of Averages" STDIN STDOUT OPTIONS: -m = remove mean from data -a = calculate analysis parameters -h = print help message Points per Segment = number points per analysis segment Number of Averages = number of segments to average STDIN: The input file is a free-format ASCII text file with three columns. The input file can be created with program breadcd: time, amplitude channel 1, amplitude channel 2 There may be more then one data set in the input file. Each data set in separated by 2 empty lines. STDOUT: The output file is a free-format ASCII text file with four columns: Frequency (hz) and Real part and Imag part and magnitude There may be more then one data set in the output file. Each data set in separated by 2 empty lines. This file can be processed by gnuplot.
This program calculate cubic spline curve fit to time series data. USAGE: cspline [OPTIONS] STDIN STDOUT OPTIONS: -h prints help message -t number = delta time increment default delta time increment = 0.01 seconds STDIN: The input file is a free-format ASCII text file with two columns: time (sec) and amplitude (eu) eu is engineering unit. There may be more then one data set in the input file. Each data set is separated by 2 empty lines STDOUT: The output file is a free-format ASCII text file with two columns: time (sec) and amplitude (eu) eu is engineering unit. There may be more then one data set in the output file. Each data set in separated by 2 empty lines. This file can be processed by gnuplot.
This program calculates cumulative rms for autospectral density sampled at even delta frequency increments. USAGE: cumrms [OPTIONS] STDIN STDOUT OPTIONS: -h prints help message STDIN: The input file is a free-format ASCII text file with two columns: Frequency (hz) and autospectral density (eu^2/hz) eu is engineering unit. There may be more then one data set in the input file. Each data set in separated by 2 blank lines STDOUT: The output file is a free-format ASCII text file with two columns: Frequency (hz) and cumulative RMS Amplitude (eu) eu is engineering unit. There may be more then one data set in the output file. Each data set in separated by 2 blank lines. This file can be processed by gnuplot.
This program combines sine dwell, sine sweep, and random vibration environments into a final damage equivalent test level. The final equivalent test level can be a sine dwell, sine sweep, or random test type.
Calculates dBA level for octave or third octave band sound pressure level data. USAGE: dba [OPTIONS] STDIN STDOUT OPTIONS: -h prints help message STDIN: The input file is a free-format ASCII text file with two columns: Frequency (hz) and SPL (dB) dB is referenced to 0.00002 N/M^2 There may be more then one data set in the input file. Each data set in separated by 2 blank lines STDOUT: The output file is a free-format ASCII text file Sound pressure (dBA) There may be more then one data set in the output file. Each data set is separated by 2 blank lines.
This program applies a dB scale factor to input data file. USAGE: dbfactor [OPTIONS] TYPE FACTOR STDIN STDOUT OPTIONS -h prints help message -m apply - sign to factor TYPE = spl to apply dB factor SPL spectra TYPE = sine to factor sine(g's) spectra by factor dB TYPE = psd to factor psd(g's) spectra by factor dB FACTOR = in dB to apply to spectra STDIN: The input file is a free-format ASCII text file with two columns: Frequency(hz) and autospectral density (eu^2/hz) or sine amplitude (g's) or SPL eu is engineering unit. There may be more then one data set in the input file each data set is separated by 2 empty lines STDOUT: The output file is a free-format ASCII text file with two columns: Frequency(hz) and autospectral density (eu^2/hz) or sine amplitude (g's) or SPL eu is engineering unit. There may be more then one data set in the output file. Each data set in separated by 2 empty lines. This file can be processed by gnuplot.
This program decimates a time history by getting every other data point. Note you should low pass data to avoid aliasing. USAGE: desample [OPTIONS] STDIN STDOUT OPTIONS -h prints help message STDIN: The input file is a free-format ASCII text file with two columns: Time(seconds) and amplitude (eu) eu is engineering unit. There may be more then one data set in the input file each data set is separated by 2 empty lines STDOUT: The output file is a free-format ASCII text file with two columns: Time(seconds) and amplitude (eu) eu is engineering unit. There may be more then one data set in the output file. Each data set in separated by 2 empty lines. This file can be processed by gnuplot.
This program interpolates data dropouts in time history or psd data. USAGE: dropout [OPTIONS] STDIN STDOUT OPTIONS: -h prints help message -p autospectral density file type (default file type is time history) -f number = data level indicating dropout large negative number for time history (default = -1000) small positive number for autospectral density (default = 1.0E-18) STDIN: The input file is a free-format ASCII text file with two columns: Frequency(hz) or time(sec) and amplitude (eu) or amplitude (eu^2/hz) eu is engineering unit. There may be more then one data set in the input file each data set is separated by 2 empty lines STDOUT: The output file is a free-format ASCII text file with two columns: Frequency(hz) or time(sec) and amplitude (eu) or amplitude (eu^2/hz) eu is engineering unit. There may be more then one data set in the output file. Each data set in separated by 2 empty lines. This file can be processed by gnuplot.
This program calculates a endurance vibration envelope from a ASCII database that contains flight test data. USAGE: enddb [OPTIONS] conditions channels MISSIONS STDIN STDOUT OPTIONS: -h = print help message STDIN: The database file is a free-format ASCII text file. GOOD/BAD scale_factor f0 df np flight channel psd... GOOD/BAD = good or bad channel scale = scale factor for psd fo = starting frequency(hz) df = delta frequency increment(hz) condition = flight condition(integer number) channel = accelerometer channel(integer number) psd ... = autospectral density amplitudes for np analysis frequencies MISSIONS: Free-format ASCII text file in format. hours condition 1 condition 2 ... STDOUT: The output file is a free-format ASCII text file with two columns: Frequency(hz) and autospectral density(eu^2/hz) eu is engineering unit.
This program calculates endurance vibration envelope with constant Q smoothing from ASCII database that contains the flight test data. USAGE: %s [OPTION] flights channels missionfile STDIN STDOUT\n", prog); [OPTIONS] -h = print help message -b = bandwidth increment in delta octave default = 3 for third octave smoothing -t = test duration in hours default = 5 hours -d = damping exponet default = 2 -s = slope of fatigue curve default = 5.68 STDIN: The database file is a free-format ASCII text file. GOOD/BAD f0 df np flight channel psd... GOOD/BAD = good or bad channel scale = scale factor for psd fo = starting frequency(hz) df = delta frequency increment(hz) flight = flight condition(integer number) channel = accelerometer channel(integer number) psd ... = PSD amplitudes for np analysis frequencies Input mission file: The input file is a free-format ASCII text file. hours flight1 flight2 ... STDOUT: The output file is a free-format ASCII text file with two columns: Frequency(Hz) and PSD(EU^2/Hz) EU is engineering unit.
This program calculates a envelope of vibration spectra. USAGE: envel [OPTIONS] STDIN STDOUT OPTIONS: -h prints help message STDIN: The input file is a free-format ASCII text file. Frequency(hz) and amplitude (eu) or autospectral density (eu^2/hz) eu is engineering unit. STDOUT: The output file is a free-format ASCII text file with two columns: Frequency(hz) and amplitude (eu) or autospectral density (eu^2/hz) eu is engineering unit.
This program calculates every spectral density given input time history. Usage: esd [OPTIONS] "Points per Segment" averages STDIN STDOUT OPTIONS: -m = remove mean from data -a = calculate analysis parameters -h = print help message Point per Segment = number points per analysis segment averages = number of segments to average STDIN: The input file is a free-format ASCII text file with two columns: Time (seconds) and Amplitude (eu) eu is engineering unit There may be more then one data set in the input file each data set in separated by 2 empty lines. STDOUT: The output file is a free-format ASCII text file with two columns: Frequency (hz) and energy spectral density (eu^2/hz) eu is engineering unit. There may be more then one data set in the output file. Each data set in separated by 2 empty lines. This file can be processed by gnuplot.
This program calculates level exceedances of a time history or a rainflow cycle count. USAGE: exceed [OPTIONS] STDIN STDOUT OPTIONS: -h prints help message -b number = number of bins -x use xrange for histogram (default = yrange) STDIN: The input file is a free-format ASCII text file with 2 columns: Time (sec) or rainflow range and amplitude (eu) or rain flow mean eu = engineering unit There may be more then one data set in the input file each data set in separated by 2 blank lines STDOUT: The output file is a free-format ASCII text file with two columns: bin center point and bin count There may be more then one data set in the output file. Each data set in separated by 2 blank lines. This file can be processed by gnuplot.
This program calculates fast Fourier transform given time history input time history. The input and output data are in the same engineering units. Usage: fft [OPTIONS] STDIN STDOUT OPTIONS: -i = real and imag output -h = print help message STDIN: The input file is a free-format ASCII text file with two columns: Time(sec) Amplitude There may be more then one data set in the input file each data set in separated by 2 empty lines. STDOUT: The output file is a free-format ASCII text file with three columns: Frequency(Hz) "Amplitude and Phase or Real and Imag output" There may be more then one data set in the output file. Each data set in separated by 2 empty lines. This file can be processed by gnuplot.
This program does fatigue life calculation given a rainflow cycle count. The fatigue strength exponent and fatigue strength coefficient are input from the command line. USAGE: fstress [OPTIONS] "Fatigue Strength Coefficient" "Fatigue Strength Exponent" STDIN STDOUT OPTIONS: -h prints help message STDIN: The input file is a free-format ASCII text file with two columns: Range Count and Mean Count There may be more then one data set in the input file each data set is separated by 2 empty lines STDOUT: life in blocks There may be more then one data set in the output file. Each data set in separated by 2 empty lines. This file can be processed by gnuplot.
This program calculates transfer function of single degree of freedom system. USAGE: gain STDOUT Input to this program is Interactive. STDOUT: The output file is a free-format ASCII text file with two columns: Frequency(hz) and Gain
This program generates a sinusoidal signal. USAGE: gsine td sr freq amp [OPTION] STDOUT td = time duration sr = sample rate freq = frequency amp = amplitude OPTIONS -h prints help message STDOUT: The output file is a free-format ASCII text file with two columns: Time(sec) Amplitude(eu) eu is engineering unit.
This program generates a white noise signal. USAGE: gwhite td sr mean stddev [OPTION] STDOUT td = time duration sr = sample rate mean = mean of signal stddev = standard deviation of signal OPTIONS -h prints help message STDOUT: The output file is a free-format ASCII text file with two columns: Time(sec) Amplitude(eu) eu is engineering unit.
This program applies a Hanning window to spectra. Usage: hann [OPTIONS] STDIN STDOUT -h = print help message STDIN: The input file is a free-format ASCII text file with two columns: Frequency (hz) and autospectral density (eu^2/hz) There may be more then one data set in the input file each data set in separated by 2 empty lines. STDOUT: The output file is a free-format ASCII text file with two columns: Frequency (hz) and autospectral density (eu^2/hz) eu is engineering unit. There may be more then one data set in the output file. Each data set in separated by 2 empty lines. This file can be processed by gnuplot.
This program calculates histograms. Calculates a normal distribution function for comparison by using the -p command line switch. USAGE: histogram [OPTIONS] STDIN STDOUT OPTIONS: -h prints help message -b number = number of bins -p calculate density -x use xrange for histogram(default = yrange) STDIN: The input file is a free-format ASCII text file with 2 columns: Time (sec) or rainflow range and amplitude (eu) or rain flow mean eu = engineering unit There may be more then one data set in the input file each data set in separated by 2 blank lines STDOUT: The output file is a free-format ASCII text file with two or three columns: bin center point and ( bin count and normal distribution for -p switch ) There may be more then one data set in the output file. Each data set in separated by 2 blank lines. This file can be processed by gnuplot.
This program calculates butterworth high-pass filter of time history. USAGE: hpf [OPTIONS] "Frequency Cutoff" "Sample Rate" STDIN STDOUT OPTIONS: -h prints help message Frequency Cutoff = high-pass frequency Sample Rate = sample rate for time history STDIN: The input file is a free-format ASCII text file with two columns: Time (sec) and Amplitude (eu) eu = engineering unit There may be more then one data set in the input file each data set is separated by 2 empty lines STDOUT: The output file is a free-format ASCII text file with two columns: Time (sec) and Amplitude (eu) There may be more then one data set in the output file. Each data set in separated by 2 empty lines. This file can be processed by gnuplot.
This program calculates inverse fast Fourier transform. The input and output data are in the same engineering units. Usage: ifft [OPTIONS] STDIN STDOUT OPTIONS: -h = print help message STDIN: The input file is a free-format ASCII text file with three columns: Frequency(Hz) Real_Magitude Imag_Magitude There may be more then one data set in the input file each data set in separated by 2 empty lines. STDOUT: The output file is a free-format ASCII text file with three columns: Time(sec) Real_Magitude Imag_Magitude There may be more then one data set in the output file. Each data set in separated by 2 empty lines. This file can be processed by gnuplot.
Calculate log-log plot parameters for amplitude spectra plotted as straight lines on log-log coordinates. Slopes are calculated in dB/octave. dB is defined as 20*log10(A2/A1) for amplitude spectra. USAGE: ilpa [OPTIONS] STDIN STDOUT OPTIONS: -h prints help message STDIN: To calculate Amplitude 2 input: 0 slope(dB/octave) freq1 amp1 freq2 To calculate Frequency 2 input: 1 slope(dB/octave) freq1 amp1 amp2 To calculate Frequency 1 input: 2 slope(dB/octave) amp1 freq2 amp2 To calculate Amplitude 1 input: 3 slope(dB/octave) freq1 freq2 amp2 To calculate slope(dB/octave) input: 4 freq1 amp1 freq2 amp2 STDOUT: The output file is a free-format ASCII text file with five columns: Freq 1 and Amplitude 1 and Freq 2 and Amplitude 2 and Slope (dB/Octave)
USAGE: ilpp [OPTIONS] STDIN STDOUT Calculate log-log plot parameters for power spectra plotted as straight lines on log-log coordinates. Slopes are calculated in dB/octave. dB is defined as 10*log10(P2/P1) for power spectra. OPTIONS: -h prints help message STDIN: To calculate Amplitude 2 input: 0 slope(dB/octave) freq1 amp1 freq2 To calculate Frequency 2 input: 1 slope(dB/octave) freq1 amp1 amp2 To calculate Frequency 1 input: 2 slope(dB/octave) amp1 freq2 amp2 To calculate Amplitude 1 input: 3 slope(dB/octave) freq1 freq2 amp2 To calculate slope(dB/octave) input: 4 freq1 amp1 freq2 amp2 STDOUT: The output file is a free-format ASCII text file with five columns: Freq 1 and Amplitude 1 and Freq 2 and Amplitude 2 and Slope (dB/Octave)
This program interpolates autospectral density on log-log coordinates to integer number of delta autospectral density input spectra points. Note this program does not change the frequency resolution of the autospectral density, it just re-samples at a difference frequency spacing. USAGE: int [OPTIONS] STDIN STDOUT OPTIONS -h prints help message -f number = delta interpolation increment default delta interpolation increment = 1 hertz -l number = lower freq extrapolation limit -u number = upper freq extrapolation limit STDIN: The input file is a free-format ASCII text file with two columns: Frequency(hz) and autospectral density (eu^2/hz) eu is engineering unit. There may be more then one data set in the input file each data set is separated by 2 empty lines STDOUT: The output file is a free-format ASCII text file with two columns: Frequency(hz) and autospectral density (eu^2/hz) eu is engineering unit. Each data set in separated by 2 empty lines. This file can be processed by gnuplot.
This program interpolates amplitude data plotted on log-log coordinates to third octave band frequencies. Amplitude data can be sine(g's), psd(g^2/hz), or mobility. USAGE: intampthird [OPTIONS] STDIN STDOUT OPTIONS: -h prints help message STDIN: The input file is a free-format ASCII text file with two columns: Frequency(hz) and amplitude (eu) or amplitude (eu) or spectra (eu^2/hz) or mobility There may be more then one data set in the input file. Each data set in separated by 2 blank lines. STDOUT: The output file is a free-format ASCII text file with two columns: Frequency(hz) and amplitude (eu) or amplitude (eu) or spectra (eu^2/hz) or mobility There may be more then one data set in the output file. Each data set in separated by 2 blank lines. This file can be processed by gnuplot.
This program does integration of time history data. USAGE: integ [OPTION] SCALE STDIN STDOUT SCALE = scale factor for time history data. Use 386.1 to scale accel in g's to inches/seconds^2 before integration to get velocity in in/second OPTIONS: -h prints help message STDIN: The input file is a free-format ASCII text file with two columns: Time(seconds) and Amplitude(EU) EU is engineering unit. There may be more then one data set in the input file each data set is separated by 2 empty lines STDOUT: The output file is a free-format ASCII text file with two columns: Time(seconds) and Amplitude(EU) EU is engineering unit. There may be more then one data set in the output file. Each data set in separated by 2 empty lines. This file can be processed by gnuplot.
This program interpolates octave band SPL data to standard octave band frequencies USAGE: intspl2octave [OPTIONS] STDIN STDOUT OPTIONS -h prints help message -l number = lower freq extrapolation limit -u number = upper freq extrapolation limit STDIN: The input file is a free-format ASCII text file with two columns: Frequency(hz) and SPL There may be more then one data set in the input file each data set in separated by 2 blank lines STDOUT: The output file is a free-format ASCII text file with two columns: Frequency(hz) and SPL There may be more then one data set in the output file. Each data set in separated by 2 blank lines. This file can be processed by gnuplot.
USAGE: intspl2old [OPTIONS] STDIN STDOUT This program interpolates octave band sound pressure level data to data to old sound pressure level octave band frequencies OPTIONS -l number = lower freq extrapolation limit -u number = upper freq extrapolation limit STDIN: The input file is a free-format ASCII text file with two columns: Frequency(hz) and SPL There may be more then one data set in the input file each data set in separated by 2 blank lines STDOUT: The output file is a free-format ASCII text file with two columns: Frequency(hz) and SPL There may be more then one data set in the output file. Each data set in separated by 2 blank lines. This file can be processed by gnuplot.
This program interpolates sound pressure level data to third octave band frequencies. USAGE: intspl2third [OPTIONS] STDIN STDOUT OPTIONS -h prints help message -l number = lower freq extrapolation limit -u number = upper freq extrapolation limit STDIN: The input file is a free-format ASCII text file with two columns: Frequency(hz) and SPL There may be more then one data set in the input file each data set in separated by 2 blank lines STDOUT: The output file is a free-format ASCII text file with two columns: Frequency(hz) and SPL There may be more then one data set in the output file. Each data set in separated by 2 blank lines. This file can be processed by gnuplot.
Calculate Joint acceptance for reverberent noise on plate with simple supports. USAGE: jarev [OPTIONS] STDIN STDOUT STDIN: The output file is a free-format ASCII text file with four columns: Altitude(ft) and mode number and length x and length y OPTIONS -h prints help message STDOUT: The output file is a free-format ASCII text file with two columns: Frequency(hz) and Joint Acceptance
Calculate Longitudinal Joint acceptance upstream of protuberance for simple support beam per Wyle report 71-10 Equation 39. USAGE: joint [OPTIONS] STDIN STDOUT OPTIONS -h prints help message STDIN: Calculate given Mach use: "Mach" altitude(ft) Mach mode length_streamwise Calculate given keas use: "keas" altitude(ft) keas mode length_streamwise Calculate given ktas use: "ktas" altitude(ft) ktas mode length_streamwise Calculate given kcas use: "kcas" altitude(ft) kcas mode length_streamwise Notes: First entry is a string that defines the type of airspeed information to entry. mode is number of half waves STDOUT: The output file is a free-format ASCII text file with two columns: Frequency(hz) and Joint Acceptance
Calculation of Joint acceptance for progressive wave on simply supported beam. USAGE: %s [OPTIONS] STDIN STDOUT\n", prog); OPTIONS -h prints help message STDIN: The input file is a free-format ASCII text file with 3 columns: Altitude(ft) \"Axial mode half waves\" Lenght_X STDOUT: The output file is a free-format ASCII text file with two columns: Frequency(Hz) and Joint Acceptance There may be more then one data set in the output file. Each data set in separated by 2 blank lines. This file can be processed by gnuplot.
This program applies a 6 order Butterworth low pass filter to a time history. USAGE: lpf [OPTIONS] "Low Pass Frequency" "Time HistorySample Rate" STDIN STDOUT OPTIONS: -h prints help message STDIN: The input file is a free-format ASCII text file with two columns: Time (sec) and Amplitude (eu) eu = engineering units There may be more then one data set in the input file each data set is separated by 2 empty lines STDOUT: The output file is a free-format ASCII text file with two columns: Time (sec) and Amplitude (eu) eu = engineering units There may be more then one data set in the output file. Each data set in separated by 2 empty lines. This file can be processed by gnuplot.
Remove linear trend from time history. USAGE: ltrend [OPTIONS] STDIN STDOUT OPTIONS: -h prints help message STDIN: The input file is a free-format ASCII text file with two columns: time(sec) and amplitude(eu) eu is engineering unit. There may be more then one data set in the input file each data set is separated by 2 empty lines STDOUT: The output file is a free-format ASCII text file with two columns: time(sec) and amplitude(eu) eu is engineering unit. There may be more then one data set in the output file. Each data set in separated by 2 empty lines. This file can be processed by gnuplot.
This program calculates a composite condition and channel envelope from a ASCII database that contains the test data. The condition and channel envelopes are the ones that define the maximum envelopes of the composite plotted vibration data. USAGE: maxfc [OPTIONS] conditions channels STDIN STDOUT OPTIONS: -h = help message STDIN: The database file is a free-format ASCII text file. GOOD/BAD scale f0 df np condition channel psd... GOOD/BAD = good or bad channel scale = scale factor for psd fo = starting frequency(hz) df = delta frequency increment(hz) condition = condition (integer number) channel = accelerometer channel(integer number) psd ... = autospectral density amplitudes for np analysis frequencies Note conditions and channels are input using a colon operator to for example: 1,2,10:15,6 would expand to 1,2,10,11,12,13,14,15,6 STDOUT: The output file is a free-format ASCII text file with three columns: Frequency(hz) and Condition and Channel There may be more then one data set in the output file. Each data set in separated by 2 blank lines. This file can be processed by gnuplot.
This program calculates a endurance vibration envelope. USAGE: nocal [OPTIONS] STDIN STDOUT OPTIONS: -h prints help message STDIN: The input file is a free-format ASCII text file with n data blocks as: segment duration in Hours Frequency(hz) and autospectral density(g^2/hz) . . . . . . Each data block in separated by 2 blank lines. STDOUT: The output file is a free-format ASCII text file with two columns: Frequency(hz) and autospectral density (eu^2/hz) eu is engineering unit. This file can be processed by gnuplot.
This program normalizes time history to start at zero time. Default is to subtract offset from time values. USAGE: normth [OPTIONS] STDIN STDOUT OPTIONS: -t calculate time values t = 0+i*dt -h prints help message STDIN: The input file is a free-format ASCII text file with two columns: Time(sec) and Amplitude (eu) eu = engineering unit There may be more then one data set in the input file each data set is separated by 2 empty lines STDOUT: The output file is a free-format ASCII text file with two columns: Time(sec) and Amplitude (eu) eu = engineering unit There may be more then one data set in the output file. Each data set in separated by 2 empty lines. This file can be processed by gnuplot.
This program calculates log normal tolerance limits. USAGE: ntl [OPTIONS] STDIN STDOUT OPTIONS: -h prints help message -b is beta the percentile limit, default = 95% possible values are 90%, 95%, and 99% -g is gama the confidence coefficient, default = 0.50 possible values are 50%, 75%, and 90% STDIN: The input file is a free-format ASCII text file. Frequency(hz) and autospectral density (eu^2/hz) eu is engineering unit. STDOUT: The output file is a free-format ASCII text file with two columns: Frequency(hz) and autospectral density (eu^2/hz) eu is engineering unit.
This program calculates overall sound pressure level. USAGE: oaspl [OPTIONS] STDIN STDOUT OPTIONS: -h prints help message STDIN: The input file is a free-format ASCII text file with two columns: Frequency(hz) and SPL(dB) There may be more then one data set in the input file each data set is separated by 2 empty lines STDOUT: OASPL There may be more then one data set in the output file. Each data set in separated by 2 empty lines. This file can be processed by gnuplot.
This program converts octave band sound pressure levels to autospecral density. USAGE: oct2psd [OPTIONS] STDIN STDOUT OPTIONS: -h prints help message STDIN: The input file is a free-format ASCII text file with two columns: Frequency(hz) and octave band sound pressure levels There may be more then one data set in the input file each data set is separated by 2 empty lines STDOUT: The output file is a free-format ASCII text file with two columns: Frequency(hz) and autospectral density (eu^2/hz) eu is engineering unit. There may be more then one data set in the output file. Each data set in separated by 2 empty lines. This file can be processed by gnuplot.
This program multiplies complex data in real - imag format by i*omega. Usage: omathdiff [OPTIONS] STDIN STDOUT OPTIONS: -h = print help message STDIN: The input file is a free-format ASCII text file with three columns: Frequency(Hz) real imag There may be more then one data set in the input file each data set in separated by 2 empty lines. STDOUT: The output file is a free-format ASCII text file with three columns: Frequency(Hz) real imag There may be more then one data set in the output file. Each data set in separated by 2 empty lines. This file can be processed by gnuplot.
This program packs a ASCII database into a binary database. USAGE: packdb conditions channels STDIN STDOUT conditions = list of flight conditions to pack to binary database channels = list of channel numbers to pack to binary database flight or channel ranges are separated by colon's for example: 100,101:105,110 expands to: 100 101 105 106 107 108 109 110 STDIN: The database file is a free-format ASCII text file. GOOD/BAD scale x0 dx np flight channel y ... GOOD/BAD = mark data as good or bad scale = scale factor to apply to data x0 = starting time(sec) or frequency dx = delta time(sec) or frequency np = number of amplitudes flight = flight condition(integer number) channel = accelerometer channel(integer number) y = amplitude STDOUT: The output database file is a binary file. GOOD/BAD scale x0 dx np flight channel y ... GOOD/BAD = mark data as good or bad scale = scale factor to apply to data x0 = starting time(sec) or frequency dx = delta time(sec) or frequency np = number of amplitudes flight = flight condition(integer number) channel = accelerometer channel(integer number) y = amplitude
This program locates minimum and maximum peaks of time histories USAGE: peaks [OPTIONS] STDIN STDOUT OPTIONS: -h prints help message STDIN: The input file is a free-format ASCII text file with two columns: Time (sec) and amplitude (eu) eu = engineering unit There may be more then one data set in the input file each data set is separated by 2 empty lines STDOUT: The output file is a free-format ASCII text file with two columns: Time (sec) and amplitude (eu) eu = engineering unit There may be more then one data set in the output file. Each data set in separated by 2 empty lines. This file can be processed by gnuplot.
This program calculates the minimum, average, and maximum envelope of a set of vibration spectra. USAGE: envel [OPTION] STDIN STDOUT OPTIONS: -h prints help message STDIN: The input file is a free-format ASCII text file. Frequency(Hz) and PSD(EU^2/Hz) EU is engineering unit. STDOUT: The output file is a free-format ASCII text file with four columns: Frequency(Hz) minimum PSD(EU^2/Hz) average PSD(EU^2/Hz) maximum PSD(EU^2/Hz) EU is engineering unit.
Calculation of modal mass for plate with simple supports. USAGE: %s [OPTIONS] STDIN STDOUT\n", prog); OPTIONS -h prints help message STDIN: The input file is a free-format ASCII text file with 3 columns: "Mass per unit area" Lenght_X Lenght_Y STDOUT: The output file is a free-format ASCII text file: Modal Mass
Convert autospectral density of pressures to sound pressure response. USAGE: psd2spr [OPTIONS] STDIN STDOUT OPTIONS: -h prints help message -q damping Q value default damping Q = 10 -b pressure coupling factor default coupling factor Beta = 2 STDIN: The input file is a free-format ASCII text file with two columns: Frequency(hz) and autospectral density (psi^2/hz) There may be more then one data set in the input file each data set is separated by 2 empty lines STDOUT: The output file is a free-format ASCII text file with two columns: Frequency(hz) and Sound Pressure Response (psi) There may be more then one data set in the output file. Each data set in separated by 2 empty lines. This file can be processed by gnuplot.
This program calculates autospectral density response for a single degree of freedom system relative to a inertial reference frame fixed to earth given a autospectral density base motion input autospectral density spectra. USAGE: psdresponse [OPTIONS] STDIN STDOUT OPTIONS: -h prints help message -f number = natural frequency default natural frequency = 100 hertz -q number = system gain default system gain = 10 STDIN: The input file is a free-format ASCII text file with two columns: Frequency(hz) and autospectral density (eu^2/hz) eu = engineering unit. There may be more then one data set in the input file, each data set in separated by 2 empty lines. STDOUT: The output file is a free-format ASCII text file with two columns: Frequency(hz) and autospectral density (eu^2/hz) eu is engineering unit. There may be more then one data set in the output file. Each data set in separated by 2 empty lines. This file can be processed by gnuplot.
This program does rainflow cycle counting. Input file must be arranged to begin with the maximum peak for this program. Use the program rfsetup to preprocess the data for the rainflow count. The rfsetup program will arrange the data to begin with the maximum peak. USAGE: rainflow [OPTIONS] STDIN STDOUT OPTION -h prints help message STDIN: The input file is a free-format ASCII text file with two columns: Time (sec) and amplitude (eu) eu = engineering unit There may be more then one data set in the input file each data set is separated by 2 empty lines STDOUT: The output file is a free-format ASCII text file with two columns: Range and Mean There may be more then one data set in the output file. Each data set is separated by 2 empty lines.
Usage: ranpan "interactive input" STDOUT This program calculates random panel response to acoustic pressure Reference: R. H. Lyon and R. G. DeJong, "Theory and application of statistical energy analysis," Butterworth-Heinemann, Boston, 1995 STDOUT: The output file is a free-format ASCII text file with two columns: Frequency (hz) and autospectral density (g^2/hz)
Read time history or autospectral density data from ASCII database. USAGE: readdb FLIGHTS CHANNELS STDIN STDOUT FLIGHTS = list of flight conditions to parse from database CHANNELS = list of channel numbers to parse from database Note conditions and channels are input using a colon operator to for example: 1,2,10:15,6 would expand to 1,2,10,11,12,13,14,15,6 Database format: GOOD/BAD scale x0 dx np flight channel y ... STDIN: The database file is a free-format ASCII text file. GOOD/BAD scale x0 dx np flight channel y ... GOOD/BAD = mark data as good or bad scale = scale factor to apply to data x0 = starting time(sec) or frequency dx = delta time(sec) or frequency np = number of amplitudes flight = flight condition(integer number) channel = accelerometer channel(integer number) y = amplitude STDOUT: The output file is a free-format ASCII text file with two columns: time(sec or frequency) and amplitude(eu) eu is engineering unit. There may be more then one data set in the output file. Each data set in separated by 2 blank lines. This file can be processed by gnuplot.
This program reads sound pressure data from ASCII database USAGE: readspl [OPTIONS] BAND FLIGHTS CHANNELS STDIN STDOUT OPTIONS: -h = help BAND = Data type: 1 = octave band,0 = old octave band, 3 = third octave band Flights = list of conditions to parse from database Channels = list of channel numbers to parse from database Note conditions and channels are input using a colon operator to for example: 1,2,10:15,6 would expand to 1,2,10,11,12,13,14,15,6 STDIN: The database file is a free-format ASCII text file. GOOD/BAD dbscale x0 dx np flight channel y ... GOOD/BAD = mark data as good or bad dbscale = dB factor to apply to data x0 = starting frequency(hz) np = number of frequencies flight = flight condition(integer number) channel = channel(integer number) y = SPL amplitude STDOUT: The output file is a free-format ASCII text file with two columns: Frequency(hz) and SPL(dB) There may be more then one data set in the output file. Each data set in separated by 2 blank lines. This file can be processed by gnuplot.
This program calculates autospectral density response spectra for a single degree of freedom system relative to a inertial reference frame fixed to the base of the system given a autospectral density base motion input autospectral density spectra. USAGE: relpsdresponse [OPTIONS] STDIN STDOUT OPTIONS: -h prints help message -f number = natural frequency default natural frequency = 100 hertz -q number = system gain default system gain = 10 STDIN: The input file is a free-format ASCII text file with two columns: Frequency(hz) and autospectral density (eu^2/hz) eu = engineering unit. There may be more then one data set in the input file, each data set in separated by 2 empty lines. Output File: The output file is a free-format ASCII text file with two columns: Frequency(hz) and autospectral density (eu^2/hz) eu is engineering unit. There may be more then one data set in the output file. Each data set in separated by 2 empty lines. This file can be processed by gnuplot.
This program calculates relative random response spectra given autospectral density base input spectra. Response is calculated relative to a axis fixed to the base on the system. USAGE: relrs [OPTIONS] STDIN STDOUT OPTIONS: -q number = damping default damping = 10 STDIN: The input file is a free-format ASCII text file with two columns: Frequency(hz) and autospectral density(eu^2/hz) eu is engineering unit. There may be more then one data set in the input file each data set in separated by 2 blank lines STDOUT: The output file is a free-format ASCII text file with two columns: Frequency(hz) and Amplitude(eu) eu is engineering unit. There may be more then one data set in the output file. Each data set in separated by 2 blank lines. This file can be processed by gnuplot.
This program arranges time history for simplified rainflow counting by rearranging so that time history begins with maximum peak. Data should be preconditioned with program "peaks" to remove data points between peaks. This data can be processed with program rainflow to do the cycle counting. USAGE: rfsetup [OPTIONS] STDIN STDOUT OPTIONS: -h prints help message STDIN: The input file is a free-format ASCII text file with two columns: Time (sec) and amplitude (eu) eu = engineering unit There may be more then one data set in the input file each data set is separated by 2 empty lines STDOUT: The output file is a free-format ASCII text file with two columns: Time (sec) and amplitude (eu) eu = engineering unit There may be more then one data set in the output file. Each data set in separated by 2 empty lines. This file can be processed by gnuplot.
This program removes mean from time history. Usage: rmean [OPTIONS] STDIN STDOUT OPTIONS: -h = print help message STDIN: The input file is a free-format ASCII text file with two columns: Time(sec) and Amplitude(eu) eu is engineering unit. There may be more then one data set in the input file each data set in separated by 2 empty lines. STDOUT: The output file is a free-format ASCII text file with two columns: Time(sec) and Amplitude(eu) eu is engineering unit. There may be more then one data set in the output file. Each data set in separated by 2 empty lines. This file can be processed by gnuplot.
This program calculates root mean square of a autospectral spectral density sampled at small even delta frequency increments using trapezoidal integration. USAGE: rms [OPTIONS] STDIN STDOUT OPTIONS -h prints help message STDIN: The input file is a free-format ASCII text file with two columns: Frequency(hz) and autospectral density(eu^2/hz) eu = engineering unit. There may be more then one data set in the input file each data set in separated by 2 blank lines STDOUT: The output file is a free-format ASCII text file: RMS(eu) eu = engineering unit. There may be more then one data set in the output file. Each data set in separated by 2 blank lines.
Calculate RMS in octave bands from input autospectral density. USAGE: rmsoctave [OPTIONS] STDIN STDOUT OPTIONS: -h prints help message STDIN: The input file is a free-format ASCII text file with two columns: Frequency(hz) and autospectral density(eu^2/hz) eu is engineering unit. There may be more then one data set in the input file each data set in separated by 2 blank lines STDOUT: The output file is a free-format ASCII text file with two columns: Frequency(hz) and RMS(eu) eu is engineering unit. There may be more then one data set in the output file. Each data set in separated by 2 blank lines. This file can be processed by gnuplot.
This program calculates RMS for time history. USAGE: rmsth [OPTIONS] STDIN STDOUT OPTIONS: -h prints help message STDIN: The input file is a free-format ASCII text file with two columns: Frequency(hz) and amplitude(eu) eu is engineering unit. There may be more then one data set in the input file each data set in separated by 2 blank lines STDOUT: The output file is a free-format ASCII text file RMS There may be more then one data set in the output file. Each data set in separated by 2 blank lines. This file can be processed by gnuplot.
This program calculates rms level in third octave bands. USAGE: rmsthird [OPTIONS] STDIN STDOUT OPTIONS: -h prints help message -s scale RMS(G's) to RMS(M^2/sec) for comparison to MIL-STD-1472 criteria STDIN: The input file is a free-format ASCII text file with two columns: Frequency(hz) and autospectral density(eu^2/hz) eu is engineering unit. There may be more then one data set in the input file each data set in separated by 2 blank lines STDOUT: The output file is a free-format ASCII text file with two columns: Frequency(hz) and RMS(eu) eu is engineering unit. There may be more then one data set in the output file. Each data set in separated by 2 blank lines. This file can be processed by gnuplot.
This program estimates standard deviation of stress response in structure USAGE: rstress STDOUT Interactive input STDOUT: The output file is a free-format ASCII text file: Standard Deviation Stress = fn
This program scales autospectral density data file scaled in measurement units to engineering units USAGE: scale [OPTIONS] STDIN STDOUT OPTIONS: -h prints help message STDIN: The input file is a free-format ASCII text file with a header line and two columns: Header line: X1, X2, dX, Y1, Y2, dY X1 = x axis lower limit X2 = x axis upper limit dX = x axis length in measurement units Y1 = y axis lower limit Y2 = y axis upper limit dY = y axis length in measurement units Data: Frequency(MU) and autospectral density(MU) MU = measurement units There may be more then one data set in the input file each data set in separated by 2 blank lines STDOUT: The output file is a free-format ASCII text file with two columns: Frequency(hz) and autospectral density(eu^2/hz) eu is engineering unit. There may be more then one data set in the output file. Each data set in separated by 2 blank lines. This file can be processed by gnuplot.
This program scales SPL data file scaled in measurement units to engineering units USAGE: scaledb [OPTIONS] STDIN STDOUT OPTIONS: -h prints help message STDIN: The input file is a free-format ASCII text file with a header line and two columns: Header line: X1, X2, dX, Y1, Y2, dY X1 = x axis lower limit X2 = x axis upper limit dX = x axis length in measurement units Y1 = y axis lower limit Y2 = y axis upper limit dY = y axis length in measurement units Data: Frequency(MU) and autospectral density(MU) MU = measurement units There may be more then one data set in the input file each data set in separated by 2 blank lines STDOUT: The output file is a free-format ASCII text file with two columns: Frequency(hz) and SPL(dB) eu is engineering unit. There may be more then one data set in the output file. Each data set in separated by 2 blank lines. This file can be processed by gnuplot.
This program applies a scale factor to input data file. USAGE: sfactor [OPTIONS] factor STDIN STDOUT OPTIONS: -h prints help message STDIN: The input file is a free-format ASCII text file with two columns: Frequency(hz) and autospectral density(eu^2/hz) eu is engineering unit. There may be more then one data set in the input file each data set is separated by 2 empty lines STDOUT: The output file is a free-format ASCII text file with two columns: Frequency(hz) and autospectral density(eu^2/hz) eu is engineering unit. There may be more then one data set in the output file. Each data set in separated by 2 empty lines. This file can be processed by gnuplot.
This program calculates separated flow power spectral density based on data from Wyle laboratories report 70-10 Equation 10. USAGE: sflow10 [OPTIONS] STDIN STDOUT OPTIONS: -h prints help message STDIN: Input is one of the following options. "mach" altitude(ft) mach "distance from leading edge" "RMS pressure / dynamic pressure" "keas" altitude(ft) keas "distance from leading edge" "RMS pressure / dynamic pressure" "ktas" altitude(ft) ktas "distance from leading edge" "RMS pressure / dynamic pressure" "kcas" altitude(ft) kcas "distance from leading edge" "RMS pressure / dynamic pressure" STDOUT: The output file is a free-format ASCII text file with two columns: Frequency(hz) and autospectral density(psi^2/hz) There may be more then one data set in the output file. Each data set in separated by 2 blank lines. This file can be processed by gnuplot.
This program calculates attached flow power spectral density based on data from Wyle laboratories report 70-10. USAGE: sflow10 [OPTIONS] STDIN STDOUT OPTIONS: -h prints help message STDIN: Input is one of the following options. "mach" altitude(ft) mach "distance from leading edge" "keas" altitude(ft) keas "distance from leading edge" "ktas" altitude(ft) ktas "distance from leading edge" "kcas" altitude(ft) kcas "distance from leading edge" STDOUT: The output file is a free-format ASCII text file with two columns: Frequency(hz) and autospectral density(psi^2/hz) There may be more then one data set in the output file. Each data set in separated by 2 blank lines. This file can be processed by gnuplot.
This program calculates separated flow power spectral density based on data from Wyle laboratories report 70-10 Equation 43. USAGE: sflow10 [OPTIONS] STDIN STDOUT OPTIONS: -h prints help message STDIN: Input is one of the following options. "mach" altitude(ft) mach "diameter of protuberance" "RMS pressure / dynamic pressure" "keas" altitude(ft) keas "diameter of protuberance" "RMS pressure / dynamic pressure" "ktas" altitude(ft) ktas "diameter of protuberance" "RMS pressure / dynamic pressure" "kcas" altitude(ft) kcas "diameter of protuberance" "RMS pressure / dynamic pressure" STDOUT: The output file is a free-format ASCII text file with two columns: Frequency(hz) and autospectral density(psi^2/hz) There may be more then one data set in the output file. Each data set in separated by 2 blank lines. This file can be processed by gnuplot.
This program generates a half sine or terminal sawtooth shock pulse time history. half sine pulse is default USAGE: shockth [OPTIONS] STDOUT -h prints help message -s terminal sawtooth pulse: default is a half sine pulse -d number = pulse duration default duration = 0.011 seconds -g number = pulse level default level = 15 g
Calculates SIL level for octave octave band sound pressure level data. USAGE: sil [OPTIONS] STDIN STDOUT OPTIONS: -c = common octane bands 848, 1696, and 3392 hertz -t = three band method 500, 1000, and 2000 hertz default is 500, 1000, 2000, and 4000 hertz bands -h prints help message Calculates SIL level. STDIN: The input file is a free-format ASCII text file with two columns: Frequency(hz) and SPL eu is engineering unit. There may be more then one data set in the input file each data set in separated by 2 blank lines STDOUT: The output file is a free-format ASCII text file SIL There may be more then one data set in the output file. Each data set in separated by 2 blank lines.
This program applies constant Q smoothing to autospectral density. USAGE: smooth [OPTIONS] STDIN STDOUT OPTIONS: -h prints help message -b fractional octave: default third = 3 6 would be sixth octave ... STDIN: The input file is a free-format ASCII text file with two columns: Frequency(hz) and autospectral density(eu^2/hz) eu is engineering unit. There may be more then one data set in the input file each data set in separated by 2 blank lines STDOUT: The output file is a free-format ASCII text file with two columns: Frequency(hz) and autospectral density(eu^2/hz) eu is engineering unit. There may be more then one data set in the output file. Each data set in separated by 2 blank lines. This file can be processed by gnuplot.
This program calculates spl level in octave bands from pressure autospectral density data. USAGE: sploctave [OPTIONS] STDIN STDOUT OPTIONS: -h prints help message STDIN: The input file is a free-format ASCII text file with two columns: Frequency(hz) and autospectral density(psi^2/hz) There may be more then one data set in the input file each data set in separated by 2 blank lines STDOUT: The output file is a free-format ASCII text file with two columns: Frequency(hz) and spl in octave bands reference 2.9e-9 psi There may be more then one data set in the output file. Each data set in separated by 2 blank lines. This file can be processed by gnuplot.
This program calculates spl level in third octave bands from pressure autospectral density data. USAGE: splthird [OPTIONS] STDIN STDOUT OPTIONS: -h prints help message STDIN: The input file is a free-format ASCII text file with two columns: Frequency(hz) and autospectral density(psi^2/hz) There may be more then one data set in the input file each data set in separated by 2 blank lines Output File: The output file is a free-format ASCII text file with two columns: Frequency(hz) and spl in third octave band reference 2.9e-9 psi There may be more then one data set in the output file. Each data set in separated by 2 blank lines. This file can be processed by gnuplot.
This program scales random vibration test spectra for test duration and total mission profile duration USAGE: interactive
This program calculates shock response spectra. The calculated acceleration is referenced to a frame of fixed to earth. USAGE: srsa [OPTIONS] STDIN STDOUT OPTIONS: -q = damping default damping of Q = 10 -f = starting frequency default = 10 hz -d = delta frequency default = 1/20 octave increment STDIN: The input file is a free-format ASCII text file with two columns: Frequency(hz) and Amplitude (eu) eu = engineering unit There may be more then one data set in the input file each data set in separated by 2 empty lines. STDOUT: The output file is a free-format ASCII text file with two columns: Frequency(hz) and Amplitude (eu) eu = engineering unit There may be more then one data set in the output file. Each data set in separated by 2 empty lines. This file can be processed by gnuplot.
This program calculates equivalent static acceleration shock response spectra. Reference: Kelly, R.D. and Richman, G., \"Principles and Techniques of Shock Data Analysis\", The Shock and Vibration Information Center, pages 145-146 USAGE: srse [OPTIONS] STDIN STDOUT OPTIONS: -h = help -r = calculate relative deflection response(inches) -g = factor to normalize acceleration by gravity default = 386.1 -q = damping default damping of Q = 10 -f = starting frequency default = 10 Hz -d = octave increment for analysis frequencies default = 1/24 octave increment enter increment as X to get 1/X frequency spacing STDIN: The input file is a free-format ASCII text file with two columns: time(sec) and amplitude(g's) There may be more then one data set in the input file Each data set in separated by 2 empty lines. STDOUT: The output file is a free-format ASCII text file with two columns: Frequency(Hz) and Response(g's) There may be more then one data set in the output file. Each data set in separated by 2 empty lines. This file can be processed by gnuplot.
This program calculates shock response spectra given a force shock pulse input. USAGE: srse [OPTIONS] STDIN STDOUT OPTIONS: -h = help -m = replace default mass default mass = unit mass -q = damping default damping of Q = 10 -g = factor to normalize acceleration by gravity default = 386.1 -f = starting frequency default = 10 Hz -d = octave increment for analysis frequencies default = 1/24 octave increment enter increment as X to get 1/X frequency spacing STDIN: The input file is a free-format ASCII text file with two columns: time(sec) and amplitude(force) There may be more then one data set in the input file Each data set in separated by 2 empty lines. STDOUT: The output file is a free-format ASCII text file with 4 columns: Frequency(Hz) Displacement Velocity Acceleration There may be more then one data set in the output file. Each data set in separated by 2 empty lines. This file can be processed by gnuplot.
This program converts sound spectra level (referenced to 2.9e-9 psi squared) to power spectra density (psi^2/hz). USAGE: ssl2psd [OPTIONS] STDIN STDOUT OPTIONS: -h prints help message STDIN: The input file is a free-format ASCII text file with two columns: Frequency(hz) and SSL(dB) eu is engineering unit. There may be more then one data set in the input file each data set is separated by 2 empty lines STDOUT: The output file is a free-format ASCII text file with two columns: Frequency(hz) and autospectral density(psi^2/hz) eu is engineering unit. There may be more then one data set in the output file. Each data set in separated by 2 empty lines. This file can be processed by gnuplot.
This program converts third octave SPL data octave bands SPL. USAGE: third2octave [OPTIONS] STDIN STDOUT OPTIONS: -h prints help message STDIN: The input file is a free-format ASCII text file with two columns: Frequency(hz) and third octave SPL There may be more then one data set in the input file each data set in separated by 2 blank lines STDOUT: The output file is a free-format ASCII text file with two columns: Frequency(hz) and octave SPL There may be more then one data set in the output file. Each data set in separated by 2 blank lines. This file can be processed by gnuplot.
This program converts a third octave SPL to a autospectral density USAGE: third2psd [OPTIONS] STDIN STDOUT OPTIONS: -h prints help message STDIN: The input file is a free-format ASCII text file with two columns: Frequency(hz) and third octave SPL There may be more then one data set in the input file each data set is separated by 2 empty lines STDOUT: The output file is a free-format ASCII text file with two columns: Frequency(hz) and autospectral density(eu^2/hz) eu is engineering unit. There may be more then one data set in the output file. Each data set in separated by 2 empty lines. This file can be processed by gnuplot.
This program calculates running RMS level for input time history. USAGE: rmsth tr to [OPTION] STDIN STDOUT tr = lenght of RMS analysis time slice in seconds to = RMS time slice overlap in seconds OPTIONS: -h prints help message STDIN: The input file is a free-format ASCII text file with two columns: Frequency(Hz) and amplitude(EU) EU is engineering unit. There may be more then one data set in the input file each data set in separated by 2 blank lines STDOUT: The output file is a free-format ASCII text file with two columns: time RMS There may be more then one data set in the output file. Each data set in separated by 2 blank lines. This file can be processed by gnuplot.
This program splits a time history into segments. USAGE: thsplit [OPTIONS] STDIN STDOUT OPTIONS: -h prints help message -s input segment size default = 2048 STDIN: The input file is a free-format ASCII text file with two columns: time(sec) and amplitude (eu) eu = engineering unit There may be more then one data set in the input file each data set is separated by 2 empty lines STDOUT: The output file is a free-format ASCII text file with two columns: time(sec) and amplitude (eu) eu = engineering unit There may be more then one data set in the output file. Each data set in separated by 2 empty lines. This file can be processed by gnuplot.
Time reverse time series on time axis. Can be used with programs lpf and hpf to do phase correction for filtering operation. USAGE: timer [OPTIONS] STDIN STDOUT OPTIONS: -h prints help message STDIN: The input file is a free-format ASCII text file with two columns: time(sec) and amplitude (eu) eu = engineering unit There may be more then one data set in the input file each data set is separated by 2 empty lines STDOUT: The output file is a free-format ASCII text file with two columns: time(sec) and amplitude (eu) eu = engineering unit There may be more then one data set in the output file. Each data set in separated by 2 empty lines. This file can be processed by gnuplot.
This program truncates autospectral density data file to power of 2. USAGE: tp2 [OPTIONS] STDIN STDOUT OPTIONS -h prints help message STDIN: The input file is a free-format ASCII text file with two columns: Frequency(hz) and amplitude (eu) eu = engineering unit There may be more then one data set in the input file each data set in separated by 2 blank lines STDOUT: The output file is a free-format ASCII text file: Frequency(hz) and amplitude (eu) eu = engineering unit There may be more then one data set in the output file. Each data set in separated by 2 blank lines.
This program calculates transfer function of time history for two input channels. Usage: transfun [OPTIONS] Points per Segment" averages STDIN STDOUT Points per Segment = number points per analysis segment averages = number of segments to average OPTIONS: -m = remove mean from data -a = calculate analysis parameters -h = print help message STDIN: The input file is a free-format ASCII text file with three columns: Frequency(hz) and Amplitude(eu) and Amplitude(eu) eu = engineering unit There may be more then one data set in the input file each data set in separated by 2 empty lines. STDOUT: The output file is a free-format ASCII text file with two columns: Frequency(hz) and Gain There may be more then one data set in the output file. Each data set in separated by 2 empty lines. This file can be processed by gnuplot.
This program calculates sound transmission loss for a panel per "Frank Fahy, "sound and Structural vibration radiation, transmission and response," equation 4.38b. USAGE: interactive
This program selects time slice from time history file. USAGE: tslice [OPTIONS] "starting time" "ending time" STDIN STDOUT OPTIONS -h prints help message STDIN: The input file is a free-format ASCII text file with two columns: Time(sec) and Amplitude(eu) eu is engineering unit. There may be more then one data set in the input file each data set is separated by 2 empty lines STDOUT: The output file is a free-format ASCII text file with two columns: Time(sec) and Amplitude(eu) eu is engineering unit. There may be more then one data set in the input file There may be more then one data set in the output file. Each data set in separated by 2 empty lines. This file can be processed by gnuplot.
Unpacks binary time history or autospectral density database to ASCII database. USAGE: unpackdb STDIN STDOUT STDIN: The database file is a binary database file. GOOD/BAD scale x0 dx np flight channel y ... GOOD/BAD = mark data as good or bad scale = scale factor to apply to data x0 = starting time(sec) or frequency dx = delta time(sec) or frequency np = number of amplitudes flight = flight condition(integer number) channel = accelerometer channel(integer number) y = amplitude STDOUT: The database file is a ASCII database file. GOOD/BAD scale x0 dx np flight channel y ... GOOD/BAD = mark data as good or bad scale = scale factor to apply to data x0 = starting time(sec) or frequency dx = delta time(sec) or frequency np = number of amplitudes flight = flight condition(integer number) channel = accelerometer channel(integer number) y = amplitude
List descriptions of programs in vibtoolbox package. USAGE: interactive
This program transforms velocity autospectral density to displacement autospectral density USAGE: vsddsd [OPTIONS] STDIN STDOUT OPTIONS: -h prints help message STDIN: The input file is a free-format ASCII text file with two columns: Frequency(hz) and autospectral density(ips^2/hz) There may be more then one data set in the input file each data set in separated by 2 blank lines STDOUT: The output file is a free-format ASCII text file with two columns: Frequency(hz) and autospectral density(inch^2/hz) Calculated displacement is 0 to peak eu is engineering unit. There may be more then one data set in the output file. Each data set in separated by 2 blank lines. This file can be processed by gnuplot.
This program calculates zero crossings with positive slope from autospectral density input spectra. USAGE: zerocross [OPTIONS] STDIN STDOUT OPTIONS: -h prints help message STDIN: The input file is a free-format ASCII text file with two columns: Frequency(hz) and autospectral density(eu^2/hz) eu is engineering unit. There may be more then one data set in the input file each data set in separated by 2 blank lines STDOUT: The output file is a free-format ASCII text file. zero crossing frequency in hz There may be more then one data set in the output file. Each data set in separated by 2 blank lines. This file can be processed by gnuplot.
This program adds an array of zeros to each end of a PSD or time history The input and output data are in the same engineering units. Usage: zerofill [OPTIONS] pointsSegment averages STDIN STDOUT OPTIONS: -z = number zero fill points to add -h = print help message STDIN: The input file is a free-format ASCII text file with two columns: Frequency(Hz) or seconds and PSD(EU^2/hz) or g's There may be more then one data set in the input file each data set in separated by 2 empty lines. STDOUT: The output file is a free-format ASCII text file with two columns: Frequency(Hz) or seconds and PSD(EU^2/hz) or g's EU is engineering unit. There may be more then one data set in the output file. Each data set in separated by 2 empty lines. This file can be processed by gnuplot.
This program converts a output file in gnuplot format to vibtoolbox ASSCI database file format. USAGE: %s [OPTION] STDIN STDOUT\n",prog); OPTIONS: -h prints help message -i use slice index number in place of channel number STDIN: The input file is a free-format ASCII text file with two columns: Frequency(Hz) or Time and power spectral density(eu**2/Hz) or amplitude(eu) eu is engineering unit. There may be more then one data set in the input file each data set in separated by 2 blank lines STDOUT: The output file is a free-format ASCII database file