Documentation

ACS SASSI

The latest ACS SASSI Version 3 SSI capabilities incorporate many advanced algorithms and specialized features. In comparison with the standard SASSI methodology, the ACS SASSI incorporates many additional SSI capabilities and specialized features, in addition to its much faster computational speed:

(i) Generation of three-component input acceleration time histories compatible with a given design ground response spectrum with or without time-varying correlation between the components. The user has also the option to generate acceleration histories using the complex Fourier phasing of selected acceleration records (called “seed records” in the new ASCE 4-16 standard). The software provides baseline correction and computes PSD and peak ground accelerations, velocities and displacements to be used by the analyst to check the US NRC SRP 3.7.1 requirements for the simulated ground accelerations.

(ii) Evaluation of the seismic motion incoherency and wave passage effects. ACS SASSI
Version 3 includes state-of-the-art modeling including both isotropic (radial) and anisotropic (directional) incoherency models. Both stochastic and deterministic incoherent SSI approaches could be employed for simple stick models with rigid basemats. These incoherent SSI approaches were validated by EPRI (Short et. Al, 2007) for stick models with rigid basemats, and accepted by US NRC (ISG-01, May 2008) for application to the new NPP seismic analysis. ACS SASSI includes six incoherent SSI approaches, namely, two simplified deterministic approaches that are the AS and SRSS approaches benchmarked by EPRI (Short et al., 2007), three other alternate deterministic approaches, and a rigorous stochastic simulation approach that is called “Simulation Mean” approach included in the 2006-2007 EPRI validation studies. There are seven plane-wave incoherency models that can be used: the Luco-Wong model, 1986 (theoretical, not validated model), five Abrahamson models (empirical, isotropic or anisotropic, based on the statistical dense array records) and user-defined incoherency models. The Abrahamson models include the coherency models published in 1993, 2005 (all sites, surface foundations), 2006 (all sites, embedded foundations), 2007a (rock sites, all foundations), 2007b (soil sites, surface foundations) and user-defined coherency models. The recent ACS SASSI versions include directional or anisotropic Abrahamson coherency models in addition to the isotropic or radial Abrahamson coherency models included in earlier releases and also used in the 2007 EPRI studies.

The user-defined coherence functions are useful for particular sites for which more detailed seismological information is available, or for sensitivity studies including the effect of soil layering inclination or nonuniformity in horizontal plane. They are based on specific-site response data, eventually coming from the 2D nonlinear site response analysis using the equivalent-linear iterative approach for selected soil profile slices. The slice directions should in principle correspond to the two principal orthogonal directions of the soil layer slopes, namely for the maximum and minimum soil layer slopes. Thus, the user-defined coherence functions are usually different in the two selected orthogonal horizontal directions. For the incoherent SSI analysis using refined FE models with elastic foundations, we highly recommend the use of the stochastic simulation approach that includes no intrusion in the SSI system dynamics. The AS and SRSS deterministic approaches are simplified incoherent SSI approaches that should be applied only to rigid basemat stick models, as validated by the 2007 EPRI studies (Short et al., 2007). The use of AS or SRSS to elastic foundation FE models might not be necessarily appropriate since could produce considerably biased results, mostly overly conservative, but sometime unconservative, especially for the vertical direction for which the foundation basemats are much more flexible, as stated also in the new ASCE 4-16 standard.

The SRSS approach is more difficult to apply since it has no clear convergence criteria for the required number of the incoherent spatial modes to be considered for the incoherent SSI analysis. For flexible foundations, the number of required incoherent spatial modes required to reconstruct the free-field coherency matrix could be very large, in order of several tens or even hundreds of modes on a case-by-case basis. This makes SRSS impractical and even dangerous for FE elastic foundation problems.

(iii) Nonlinear hysteretic soil behavior capability is included in the ACS SASSI main software for seismic SSI analysis using the Seed-Idriss iterative equivalent linear procedure for both the global (due to wave propagation in free-field) and the local soil nonlinearity (due to SSI effects). The nonlinear soil site response is included in SOIL module based on SHAKE methodology. For 2D and 3D SSI models, the local soil nonlinear behavior can be included using near-field soil PLANE or SOLID elements (defined in the HOUSE model input). The ACS SASSI code uses a fast SSI reanalysis (or restart) solutions for the equivalent-linear soil SSI iterations that takes advantage of the already computed soil impedance matrix available from the SSI initiation run. This feature reduces the run time per SSI iteration by a factor of 2 to 5 times depending of the foundation embedment size. For nonlinear soil SSI analyses performed in batch mode, the simultaneous X-Y-Z input effects can be considered at each SSI iteration using the COMB_XYZ_STRAIN auxiliary program.

(iv) Nonlinear hysteretic structure behavior capability is included in Option NON. The iterative equivalent-linearization can be applied for modeling the cracking and post-cracking reinforced concrete wall behavior using shell elements, or for modeling the rubber-bearing hysteretic materials for the seismic base-isolation using nonlinear spring elements. The nonlinear springs can be also applied for modeling the local pile-soil interface nonlinear effects due to pile slipping in the vertical direction that can be an important aspect during intense seismic motions.

(v) The nonuniform soil motion, or multiple seismic input motion option, includes the capability to consider variable amplitude seismic input motions. The nonuniform motion input is applicable to continuous foundations assuming that the free-field motion complex amplitude varies in the horizontal plane after specific frequency dependent spectral patterns. These patterns are described by the user using complex amplification factors (or relative transfer functions) computed with respect to the reference amplitude motion. The nonuniform motion assumption could be combined with motion incoherency and wave passage to create more realistic seismic environments.

(vi) Seven interpolation schemes for the complex responses are implemented for complex responses. The newer interpolation scheme that uses bi-cubic splines is recommended for complex FE models under incoherent seismic inputs (for which the number of SSI frequencies should be usually larger than 200). For such cases, when number of frequencies is sufficiently large, the bi-cubic spline interpolation provides most accurate results for incoherent analysis since it does not create any spurious peaks or valleys (option = 6). The bi-cubic spline interpolation should be applied only if the number of SSI frequencies is sufficiently large, so that spectral peaks are not clipped by the smooth spline interpolation. Different interpolation techniques could perform differently on a case-by-case basis, especially for highly complex FE models with coupled responses, especially for incoherent motions. The various interpolation options that are available in the code provide the structural analyst a set of powerful tools for identifying and avoiding the occurrences of spurious spectral peaks in the computed transfer functions of structural motions and stresses. The first six options (options 0-5) were implemented in the original SASSI 1982 scheme that uses a non-overlapping moving window, the university SASSI2000 scheme that uses a weighted average moving window, and four new interpolation schemes including two nonoverlapping window schemes with different shifts and two average overlapping moving window schemes with different numbers of sliding windows. To check the interpolation accuracy, convenient comparative plots of the computed TFs versus the interpolated TFs can be easily obtained using the GUI graphics.

(vii) The new Fast Flexible Volume (FFV) method provides accurate and numerically efficient SSI analysis solutions for deeply embedded structures (DES) such as small modular reactors (SMRs). The FFV method, in addition to the interaction nodes defined at the outer surface of excavation volume, includes interaction nodes defined by internal node layers within excavation volume. Using the INTGEN command, the user can automatically generate the interaction nodes for the FFV method or other methods, as FV and FI-FSIN (SM) or FI-EVBN (MSM). The FFV method speeds the SSI analysis of deeply embedded structures by up to tens of times faster than the traditional, reference FV method.

(vii) Automatic selection of additional SSI calculation frequencies that are required to improve the accuracy of the interpolated TF that is applicable to both the node acceleration TFs (ATF) and the element stress TFs (STF). This feature that can be implemented using UI commands and macros, is an important practical capability, especially for larger size FE model applications, because it saves a lot of labor effort and ensures high quality of the SSI analysis.

(viii) Visualization of complex TF variation patterns within the entire structural model for selected, SSI calculation frequencies. The complex TF patterns are visualized on the structure using colored vector plot animations including all three-directional components (red for X, green for Y and blue for Z). The TF amplitude is given by vector length, and the TF phase is given by vector orientation. This capability is extremely useful for checking the correctness of the FE modeling and understanding the structural dynamic behavior.

(ix) Computation and visualization of the amplitude TF or spectral accelerations for a selected damping value at a given SSI calculation frequency for the entire SSI model using either structural deformed shape or bubble plots. The deformed shape plots are animated structural plots with a controlled movie frame speed, so that they can be also viewed as static plots. For selected resonant frequencies, the spectral amplitudes or the ZPA values can be plotted as a deformed shape plot.

(x) Computation and visualization of structural acceleration and relative displacement time histories using structural deformed shape plots. The deformed shape plots can be static structural plots for selected times, or maximum values, or structural animations of the SSI response variation in time during the earthquake action.

(xi) Computation and contour plotting of the average nodal seismic stresses (for all six components in global coordinates) in the entire structure, or for selected parts of the structure based on the element center stresses for the SHELL and SOLID elements. Both maximum and time-varying values of nodal stresses are computed and available for plotting. The approximation is based on the assumption that element center and node stresses are equal (no shape function extrapolation is included). For sufficiently refined finite element models this approximation appears reasonable. Contour stress plotting can be either static maximum values or animated time-varying values at selected time frames (automatic frame selection is included). Maximum element center stresses values are also available in a convenient text file format.

(xii) Computation and contour plotting of seismic soil pressure on foundation walls using near-field SOLID elements. The nodal pressure is computed based on averaging of adjacent element center pressures. Both maximum and time-varying values of nodal seismic pressures are computed and available for plotting. The analyst can also automatically combine the seismic soil pressures with the static soil bearing pressures and then, plot the resultant soil pressure of foundation walls and mat. Contour seismic soil pressure plotting can be either static contour plots of maximum values or animated contour plots of time-varying values at selected time frames (an automatic frame selection capability is included).

(xiii) Post SSI analysis calculations for superposition of the co-directional SSI effects in terms of acceleration, displacement of stress time-histories and in-structure response spectra using UI commands and macros. For time histories both the algebraic summation available. For in-structure response spectra (ISRS) post-processing i) the weighted linear combination and ii) the square-root of sum of square (SRSS) for the superposition of the X, Y and Z co-directional effects are implemented. The analyst can also compute the envelope, broaden and average ISRS from multiple spectral curves. Post-processing calculations can be also used for computing the maximum structural stresses, forces and moments, and/or the maximum seismic soil pressure on walls and mat with or without including the soil static bearing pressure component.

(xiv) The most recent upgrades, starting with IKTR9, includes a fast-postprocessing option for SSI response time histories using compressed binary databases for compuations and visualizations, and including a new thick shell element to model thick walls and floors in nuclear structures. New visualization of the shell element forces and moments per element is included. The new thick shell element is accurate for both thick and thin shell elements based on MindlinReissner theory. The latest upgrade IKTR10 also includes a coordinate transformation capability for the shell elements and a principal strain computation capability for the upper and lower surfaces of shell elements, applicable to linear type shell elements. Details on the application of the new coordinate transformation are provided in Section 6 for SHELL and TSHELL elements. Also, the IKTR10 upgrade includes new aux program called Post_BinDB for post-processing SSI response time-histories stored in the binary databases for acceleration, displacement and element stress/force SSI time-history responses. It should be noted that this Post_BinDB program partially duplicates and has less capabilities than the ACS SASSI UI commands for performing post-processing of the SSI response binary databases. It might be still an attractive alternate for some users that prefer running a DOS program instead of using the UI commands.The ACS SASSI NQA Version 3 has been tested, verified, documented and released under the Ghiocel Predictive Technologies Nuclear Quality Assurance Program which is in compliance with the requirements of 10 CFR50 Appendix B, 10 CFR21, ASME NQA-1 2008 and Addenda Subpart 2.7 2009. The ACS SASSI NQA version comes with a complete set of software documentations that were developed under the quality assurance requirements of the GP Technologies NQA-1 Level Program. The ACS SASSI NQA version documentation includes the user and verification manuals and the V&V computer files for a large set of various seismic V&V problems, including shallow, embedded and buried foundations, rigid and flexible foundations, piles, subjected to various different seismic environments, different surface and body seismic waves, motion incoherency and directional wave passage along an arbitrary horizontal direction, multiple support excitations for isolated foundations, linear or nonlinear SSI analysis. The current ACS SASSI NQA Version 3 IKTR10 includes a set of 57 SSI verification and validation (V&V) problems, many of these including several subproblems. In these V&V problems, the computed SSI results using ACS SASSI are compared against benchmark results based on published analytical solutions or computed using other validated with computer programs, including SHAKE91, SASSI2000 and ANSYS. Each SSI verification problem tests a different capability of the ACS SASSI NQA code. The total number of the V&V computer input files and output files for all the SSI verification problems of the ACS SASSI NQA version is several thousands of files that require about few hundreds of MB of the hard drive space.