Background

psuedo Spectral Accelerations (pSA) values are calculated in a cython script called rspectra.

It was found that the RotD values for pSA intensity measures were smaller than expected. Upon investigation it was determined that the spectral displacements were being returned instead of the pSA values.

This did not effect components such as the basic components and geom as they were calculated and returned separately.

Changes made

In order to get back the correct pSA values, the variable u was multiplied buy the square of the circular frequency (ω) before being returned. As the maximum pSA value for each period can be determined from this data, this return value has been removed from the script.

The maximum pSA values are now calculated outside of the rspectra script. Any future issues with pSA values should be reflected here.

Further changes in the same pull request include name changes to clarify that the RotD component can be taken of any intensity measure that takes an acceleration waveform as an input.

Verification

This change was tested using a copy of a Hossack broadband binary that was available to the tester.

Because the method of calculating the geom component was changed this is also included in the below table

While only one station is presented, this held true for all stations in the resulting intensity measure files.

SourcestationcomponentpSA_0.02pSA_0.05pSA_0.1pSA_0.2pSA_0.3pSA_0.4pSA_0.5pSA_0.75pSA_1.0pSA_2.0pSA_3.0pSA_4.0pSA_5.0pSA_7.5pSA_10.0
Original rspectraWRAVgeom0.060133366752730.0891157453224540.1251331811138550.1210053035665890.1445127780306970.1038017833252560.0805236637693360.0671675791841180.0285985481439780.010283515311410.0045877555574690.0019933425131910.0011020700040630.0004431698024140.000250627974779
Original rspectraWRAV

rotd50

6.2335278555554E-075.72149302246633E-062.96559252638543E-050.0001490286142830.0003446190640280.000420198400630.0005354451913590.0009502979216580.0007554943391190.0010709226960950.0010410795034840.0008494427413230.0006978942406380.0005750259839320.000597973586903
Original rspectra with ω2WRAVgeom0.060133366752730.0891157453224540.1251331811138550.1210053035665890.1445127780306970.1038017833252560.0805236637693360.0671675791841180.0285985481439780.010283515311410.0045877555574690.0019933425131910.0011020700040630.0004431698024140.000250627974779
Original rspectra with ω2WRAVrotd500.0615224539575030.0903501963445370.1170769002010060.1470853467412020.1511668369345520.1036797996047790.0845541154748250.0666955701329520.029825721017470.0105695833546050.0045666857108750.0020959159545590.0011020704110230.0004035753942580.000236070509801
New rspectraWRAVgeom0.060133366752730.0891157453224540.1251331811138550.1210053035665890.1445127780306970.1038017833252560.0805236637693360.0671675791841180.0285985481439780.010283515311410.0045877555574690.0019933425131910.0011020700040630.0004431698024140.000250627974779
New rspectraWRAVrotd500.0615224539575030.0903501963445370.1170769002010060.1470853467412020.1511668369345520.1036797996047790.0845541154748250.0666955701329520.029825721017470.0105695833546050.0045666857108750.0020959159545590.0011020704110230.0004035753942580.000236070509801
  • No labels