Considering Second-Order Analysis in Dynamic Analyses in RFEM 6 and RSTAB 9
Technical Article
For the ultimate limit state design, EN 1998‑1, Sections 2.2.2 and 4.4.2.2 require a calculation considering the second‑order theory (P‑Δ effect). This effect may be neglected only if the interstory drift sensitivity coefficient θ is less than 0.1.
Sensitivity Factor θ
The sensitivity factor θ is defined as follows [1]: $$\mathrm\theta\;=\;\frac{\displaystyle{\mathrm P}_\mathrm{tot}\;\cdot\;{\mathrm d}_\mathrm r}{{\mathrm V}_\mathrm{tot}\;\cdot\;\mathrm h}\;$$ where
- θ = interstory drift sensitivity coefficient
- P_{tot} = the total gravity load at and above the story considered in the seismic design situation
- d_{r} = design interstory drift, evaluated as the difference of the average lateral displacements d_{S} at the top and bottom of the story under consideration; for this, the displacement is determined using the linear design response spectrum with q = 1.0
- V_{tot} = total seismic story shear determined using the linear design response spectrum
- h = interstory height
The second-order effects may be taken into account approximately by a factor equal to 1 / (1 − θ), if 0.1 < θ ≤ 0.2. For θ > 0.2, it is necessary to consider the geometric stiffness matrix when calculating eigenvalues and multi-modal response spectrum analysis.
The sensitivity factor can also be calculated in RFEM 6 and RSTAB 9. Find detailed information in the following technical article: KB 001866, "Determination of Sensitivity Coefficient to Investigate the Need for Second-Order Analysis in Dynamic Analyses".
Geometric Stiffness Matrix
For dynamic analyses, the iterative calculations for the nonlinear determination of second-order theory are not suitable. The problem can be linearized, and it is sufficient to use the geometric stiffness matrix based on axial loads to consider the second-order theory. For this, it is assumed that the vertical loads do not change due to horizontal effects, and the deformations are small compared to the building dimensions [2]. The loads to be considered should correspond to those for the seismic design situations in accordance with EN 1990 Section 6.4.3.4 [3]:
${\mathrm{E}}_{\mathrm{d}}=\mathrm{\Sigma}{\mathrm{G}}_{\mathrm{k},\mathrm{j}}\mathrm{\Sigma}{\mathrm{\psi}}_{2,\mathrm{i}}{\mathrm{Q}}_{\mathrm{k},\mathrm{i}}\left(2\right)$
- E_{d} = design value of effects
- G_{k,j} = characteristic value of the permanent action j
- Q_{k,i} = characteristic value of the variable action i
- Ψ_{2,i} = factor for the quasi‑permanent values of variable actions i
Axial tensile forces increase the stiffness, for example, in a prestressed cable. Compression forces reduce the stiffness and can lead to a singularity in the stiffness matrix. The geometric stiffness K_{g} is not dependent on the mechanical properties of the structure, but only on the member length L and axial force N. To illustrate the basic problem, a simple example of a cantilever is displayed in Figure 01. The individual mass points of the cantilever represent the individual stories of a building. The building is subjected to a dynamic analysis considering the second-order theory. The axial forces N_{i} on the individual stories i = 1...n result from the vertical forces in the seismic design situation. The story height is defined by hi.
The geometric stiffness matrix K_{g} can be derived from the static equilibrium conditions:
$\left[\begin{array}{c}{\mathrm{F}}_{\mathrm{i}}\\ {\mathrm{F}}_{\mathrm{i}+1}\end{array}\right]=\underset{{\mathbf{K}}_{\mathbf{g}}}{\underset{\u23df}{\frac{{\mathrm{N}}_{\mathrm{i}}}{{\mathrm{h}}_{\mathrm{i}}}\left[\begin{array}{rc}1.0& -1.0\\ -1.0& 1.0\end{array}\right]}}\left[\begin{array}{c}{\mathrm{u}}_{\mathrm{i}}\\ {\mathrm{u}}_{\mathrm{i}+1}\end{array}\right]\left(3\right)$
For the purpose of simplification, only the degrees of freedom of the horizontal displacement are displayed here. The derivation shown is based on the overturning moment approach due to the linear displacement application. This is a simplification for the bending element, and accurate assumption for the truss element. More precise determination of the geometric stiffness matrix for bending beams can be obtained by using the cubic displacement approach or the analytical solution of the differential equation of the bending line. More information and derivations are provided by Werkle [4]. The geometric stiffness matrix K_{g} is added to the system stiffness matrix K, and thus the modified stiffness matrix K_{mod} is obtained:
K_{mod} = K + K_{g} (4)
In the case of compression normal forces, this consequently leads to stiffness reduction.
Example: Natural Frequencies and Multi-Modal Response Spectrum Analysis Considering Second-Order Theory
The following shows how the geometric stiffness matrix can be considered in RFEM 6 and RSTAB 9. The cantilever shown in Figure 01 is used as an example. The cantilever consists of five concentrated mass points. Here, 4,000 kg act in the global X-direction in each case.
The cross‑section is IPE 300 made of S 235 material with I_{y} = 8,356 ∙ 10^{-5} m^{4} and E = 2.1 ∙ 10^{11} N/m^{2}. To be able to consider the geometric stiffness matrix in the dynamic analysis, a design situation of the seismic/mass combination is initially defined in the RFEM main program. Thus, the load combination (CO1) used subsequently is automatically created by the combinatorics integrated in the program.
The Modal Analysis add-on allows you to determine mode shapes and effective modal masses of a structure. It is possible to select the initial states in order to consider the stiffness modifications based on load cases and load combinations. Two modal analysis load cases are defined. In LC4, CO1 is imported in order to consider the geometric stiffness matrix and thus the second-order analysis. For comparison, LC3 is defined, which does not include any stiffness modifications.
The following table includes the determined natural frequencies f [Hz], natural periods T [sec], and the acceleration values S_{a} [m/s^{2}] based on the response spectrum, with and without the geometric stiffness matrix K_{g} resulting from the axial forces of CO1.
The multi-modal response spectrum analysis uses natural frequencies to determine the acceleration values from the defined response spectrum. Based on these acceleration values, the internal forces and moments of the response spectrum analysis are determined. The graphic display of a user‑defined response spectrum is shown in Figure 06, and the acceleration values S_{a} [m/s^{2}] determined from the response spectrum for each eigenvalue are listed in the table above.
In order to ensure the correct assignment of the modified frequencies, the desired modal analysis must be used as a basis when creating a response spectrum load case.
In the case of compression axial forces, the consideration of the geometric stiffness matrix leads to the reduction of the natural frequency, and can cause lower acceleration values S_{a}, as in our example. Sole modification of the natural frequencies is not enough to consider the second-order theory. In fact, this may actually lead to smaller results, which can thus be incorrect. It is very important to also use the modified stiffness matrix for the determination of internal forces and deformations. In the response spectrum analysis, the modified stiffness from the modal analysis is automatically used for determining the results of the response spectrum analysis. The deformations, internal forces and moments, and support reactions determined in the response spectrum analysis, with and without the geometric stiffness matrix, are shown in Image 08.
The consideration of the geometric stiffness matrix leads to larger deformations and internal forces. However, the resulting support loads are slightly smaller when considering the geometric stiffness matrix.
Author
Thomas Eichner, M.Sc.
Product Engineering & Customer Support
Mr. Eichner is responsible for the development of products for dynamic analysis and provides technical support for our customers.
Keywords
Second-Order Analysis Dynamics Dynamic Analysis Dynamic Calculations Response Spectrum Analysis Sensitivity Coefficient P-Δ Effect
Reference
Links
- RFEM 6 / RSTAB 9 | Dynamic Analysis | Online Manual
- KB 001866, "Determination of Sensitivity Coefficient to Investigate the Need for Second-Order Analysis in Dynamic Analyses"
- Add-ons for Dynamic Analysis in RFEM 6 / RSTAB 9
Write Comment...
Write Comment...
Contact Us
Do you have further questions or need advice? Contact us via phone, email, chat, or forum, or search the FAQ page, available 24/7.
New
This article deals with the calculation of the sensitivity coefficient of the interstory drift to investigate the relevance of the second-order analysis in dynamic analyses.
Associated Products