Response Spectrum Analysis

Response Spectrum Analysis (RSA) is a technique used to estimate the maximum response of a structure for a transient event. Maximum displacement, stresses, and/or forces may be determined in this manner.

The technique combines response spectra for a specified dynamic loading with results of a normal modes analysis. The time-history of the responses are not available.

Response spectra describes the maximum response versus natural frequency of a 1-DOF system for a specified dynamic loading. They are employed to calculate the maximum modal response for each structural mode. These modal maxima may then be combined using various methods, such as the Absolute Sum (ABS) method or the Complete Quadratic Combination (CQC) method, to obtain an estimate of the peak structural response.

RSA is a simple and computationally inexpensive method to provide an approximation of peak response, compared to conventional transient analysis. The major computational effort is to obtain a sufficient number of normal modes in order to represent the entire frequency range of input excitation and resulting response. Response spectra are usually provided by design specifications; given these, peak responses under various dynamic excitations can be quickly calculated. Therefore, it is widely used as a design tool in areas such as seismic analysis of buildings.

Governing Equations

Normal Modes Analysis

The equilibrium equation for a structure performing free vibration appears as the eigenvalue problem:

( K λ M ) A = 0

Where,
K
Stiffness matrix of the structure.
M
Mass matrix.
Damping is neglected.

The solution of the eigenvalue problem yields n eigenvalues λ i , where n is the number of degrees of freedom. The vector A is the eigenvector corresponding to the eigenvalue.

The eigenvalue problem is solved using the Lanczos or the AMSES method. Not all eigenvalues are required and only a small number of the lowest eigenvalues are normally calculated. The results of eigenvalue analysis are the fundamentals of response spectrum analysis.

Response spectrum analysis can be performed together with normal modes analysis in a single run, or eigenvalue analysis with Lanczos solver can be performed first to save eigenvalues and eigenvectors by using EIGVSAVE, which can be retrieved later by using EIGVRETRIEVE for response spectrum analysis.

Modal Combination

It is assumed each individual mode behaves like a single degree-of-freedom system. The transient response at a degree of freedom is:

u k = i A ik ψ i χ MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaebbnrfifHhDYfgasaacH8srps0l bbf9q8WrFfeuY=Hhbbf9v8qqaqFr0xc9pk0xbba9q8WqFfea0=yr0R Yxir=Jbba9q8aq0=yq=He9q8qqQ8frFve9Fve9Ff0dmeaabaqaciGa caGaaeqabaqaaeaadaaakeaacaWG1bWaaSbaaSqaaiaadUgaaeqaaO Gaeyypa0ZaaabuaeaacaWGbbWaaSbaaSqaaiaadMgacaWGRbaabeaa kiabeI8a5naaBaaaleaacaWGPbaabeaakiabeE8aJbWcbaGaamyAaa qab0GaeyyeIuoaaaa@3F74@

Where,
A
Eigenvector
ψ
Modal participation factor
χ
Response spectrum

For loading due to base acceleration, the modal participation factor can be expressed as:

ψ i = A i T M T

Where,
A
Eigenvector
M
Mass matrix
T
Rigid body motion due to excitation

In ABS modal combination, the peak response is estimated by:

u k = i | A i k | | ψ i χ |

In CQC modal combination, the peak response is estimated by:

u k = m n v m ρ m n v n

Where,
v m
m t h MathType@MTEF@5@5@+= feaahqart1ev3aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamyBamaaCa aaleqabaGaamiDaiaadIgaaaaaaa@38F9@ mode's contribution to u k MathType@MTEF@5@5@+= feaahqart1ev3aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamyDamaaBa aaleaacaWGRbaabeaaaaa@380A@ ; equal to A m k ψ k χ MathType@MTEF@5@5@+= feaahqart1ev3aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamyqamaaBa aaleaacaWGTbGaam4AaaqabaGccqaHipqEdaWgaaWcbaGaam4Aaaqa baGccqaHhpWyaaa@3D7D@ .
ρ m n
Cross-modal coefficient

The cross modal coefficient ρ m n between modes m and n is calculated as:

ρ m n = 8 ξ m ξ n ( ξ m + r n m ξ n ) r n m 1.5 ( 1 r n m 2 ) 2 + 4 ξ m ξ n r n m ( 1 + r n m 2 ) + 4 ( ξ m 2 + ξ n 2 ) r n m 2

Where,
r n m = λ n λ m
Ratio of eigenvalues of the modes
ξ m  and  ξ n
Modal damping values of the two modes

In SRSS modal combination, the peak response is estimated by:

u k = i ( A i k ψ i χ ) 2

The SRSS method is less conservative than ABS method. It is more accurate when the modes are well separated.

The NRL method combines ABS and SRSS methods. It adds the maximum modal response by ABS method and the rest of the modes by SRSS method. The peak response is estimated by:

u k =| A ik || ψ i χ |+ ji ( A jk ψ j χ ) 2

Directional Combination

In order to estimate peak response due to dynamic excitations in different directions, the peak response in each direction must be combined to obtain total peak response. Methods such as ALG (algebraic) and SRSS (square root of sum of squares) can be used.

Input Specification

Subcase Definition

An RSA subcase may be explicitly identified by setting ANALYSIS=RSPEC, but it is also implicitly chosen for any subcase containing the RSPEC data selector (when the ANALYSIS entry is not present).

The following data selectors are recognized for an RSA subcase definition.
METHOD
References an eigenvalue extraction Bulk Data Entry definition (EIGRL). Only METHOD(STRUCTURE) is supported. This reference is required.
RSPEC
References an RSPEC Bulk Data Entry where the combination rules, excitation DOF, and the input spectra are identified. This reference is required.
SDAMPING
References damping table Bulk Data Entries (TABDMP1 or TABDMP2) to specify modal damping. This reference is required.
SPC
References single point constraint Bulk Data Entries (SPCADD, SPC and SPC1). For RSA analysis, these entries define the base degrees of freedom where excitation is applied.
MPC
References multi-point constraint Bulk Data Entries (MPCADD or MPC).
STATSUB(PRELOAD)
Pre-loading is supported for a Response Spectrum analysis subcase. STATSUB(PRELOAD) can be used to identify the subcase used to apply the preloading. The eigenvalues are augmented with the pre-loading effect coming from the pre-loading subcase.

Bulk Data

Bulk Data Entries which have particular significance for RSA include:
RSPEC
Specifies combination rules, excitation DOF, and references the input spectra.
DTI,SPECSEL
Defines response spectra.
EIGRL
Defines parameters for eigenvalue extraction.
PARAM, LFREQ and PARAM, HFREQ
Defines the range of modes used in modal combinations.
TABDMP1
Specifies modal damping as a function of frequency.
TABDMP2
Specifies modal damping as a function of a range of mode indices.
SPC, SPC1, and SPCADD
Specifies base where excitation is applied and other constraints.

Example: Input

SUBCASE 100
	RSPEC = 2
	SPC = 5
	SDAMPING = 12
	METHOD = 24
$
BEGIN BULK
$
PARAM, LFREQ, 0.1
PARAM, HFREQ, 1000.
EIGRL, 24, 0.0, 1000.
RSPEC, 2, ABS, CQC, 0.1
, 99, 2.0, 1.0, 0.0, 0.0
DTI, SPECSEL, 99, , A, 2, 0., 3, 0.02,
, 4, 0.04, ENDREC
TABDMP1, 12, …
TABLED1, 2
+,…
TABLED1, 3
+,…
TABLED1, 4
+,…
ENDDATA
$

Output

Results of interest from RSA include maximum displacement, stress, strain and force. These are requested via the I/O Options Entry DISPLACEMENT, STRESS / ELSTRESS, STRAIN and FORCE / ELFORCE, respectively. More details on the supported output formats for the results can be found in Results Output by OptiStruct.

For shell elements, corner stresses are available in the H3D, PCH and OP2 file formats, while corner strains are available in the OP2 file format. For more information on the location of element outputs, refer to STRESS and STRAIN in the Reference Guide.

For bar and beam elements defined using PBARL and PBEAML respectively in RSA, von Mises Stress output is available in the .h3d file format.