Earthquake Simulation Software: Engineering Tools

The Role of Simulation in Earthquake Engineering

Earthquake simulation software enables engineers, researchers, and risk analysts to answer questions that cannot be answered through physical experimentation alone: How will an unbuilt structure perform in a once-in-500-years earthquake? What is the expected annual loss to a portfolio of 10,000 buildings distributed across a seismically active region? What ground motions should be used to design a dam in a region with incomplete historical records? The answers to these questions emerge from computational models that integrate seismological source models, wave propagation physics, soil-structure interaction, and structural dynamics into a unified analytical framework.

Finite Element Structural Analysis

At the building scale, nonlinear finite element analysis (FEA) is the primary tool for earthquake response simulation. Software packages including OpenSees (open-source, UC Berkeley), PERFORM-3D (Computers and Structures), SAP2000 Nonlinear, ABAQUS, and LS-DYNA represent the structural system as a mesh of elements with specified material constitutive laws. For earthquake analysis, elements must be capable of capturing nonlinear hysteretic behavior — the yielding, degrading, and re-stiffening that structural members undergo during cyclic inelastic deformation. The accuracy of simulation results depends critically on the calibration of material models against experimental test data from structural component tests.

Ground Motion Simulation

Probabilistic Seismic Hazard Analysis (PSHA)A method for quantifying earthquake hazard that considers all possible earthquake sources, magnitudes, and ground motion levels, expressing results as probability of exceeding specific shaking levels. requires ground motion prediction equations (GMPEs) to estimate shaking at sites, but complex site-specific assessments benefit from physics-based ground motion simulation. Stochastic simulation methods (stochastic finite fault models) represent the earthquake source as a sum of many small sub-faults radiating stochastic ground motion consistent with seismological scaling relations. Deterministic methods (finite difference and spectral element codes like Hercules, SW4, and Specfem3D) solve the three-dimensional elastic wave equation on computational grids, capturing the effects of basin geometry, velocity structure, and Fault LineThe trace of a fault on the Earth's surface, visible as a line or zone of broken rock. Active fault lines are mapped by geologists to assess earthquake hazard for nearby communities. geometry on ground motion amplification and focusing. The vulnerability functionsA mathematical function describing the probability of various damage states for a specific building type given a level of ground shaking. Essential for loss estimation models. used in regional loss estimation are derived from simulated ground motions applied to large suites of structural models.

Regional Loss Estimation: HAZUS and OpenQuake

For regional Loss EstimationThe process of predicting economic losses and casualties from a potential earthquake scenario. FEMA's HAZUS software is the standard loss estimation tool in the United States., simulation shifts from individual structures to large portfolios. FEMA's HAZUS software performs loss estimation by combining ShakeMapA USGS product that displays the distribution of ground shaking intensity after an earthquake. Combines seismograph data, ground motion models, and 'Did You Feel It?' reports. or Probabilistic Seismic Hazard Analysis (PSHA)A method for quantifying earthquake hazard that considers all possible earthquake sources, magnitudes, and ground motion levels, expressing results as probability of exceeding specific shaking levels. ground motion with a building inventory database and vulnerability functionsA mathematical function describing the probability of various damage states for a specific building type given a level of ground shaking. Essential for loss estimation models. that relate shaking intensity to damage state probabilities. Output includes estimates of casualties, displaced households, dollar losses, and infrastructure damage by census tract or county. HAZUS uses fragility functions derived from engineering judgment, laboratory testing, and post-earthquake building survey data. Fragility functions express the conditional probability of reaching or exceeding a specific damage state (none, slight, moderate, extensive, complete) given a specified shaking intensity.

The GEM OpenQuake Platform

OpenQuake Engine, the GEM Foundation's open-source loss estimation platform, extends HAZUS-style approaches with more flexible input options, probabilistic source characterization, and open exposure databases (GED4ALL global exposure model). OpenQuake can perform scenario-based loss calculations (for specified earthquake scenarios), classical PSHA-based probabilistic loss calculations, and event-set calculations using stochastic event catalogs. The platform supports Seismic Risk AssessmentThe process of evaluating earthquake hazard, building vulnerability, and potential losses for a specific area or structure. Combines hazard maps, building inventory, and damage models. at national to global scales and has been used to produce the Global Seismic Risk Map published by GEM.

Soil-Structure Interaction Simulation

The interaction between a building foundation and the surrounding soil modifies both the ground motion input to the structure and the structural response. Soil-structure interaction (SSI) simulation uses finite element or boundary element methods to model the soil domain surrounding the foundation, capturing effects like foundation embedment, soil compliance, and kinematic interaction (the inability of a rigid foundation to perfectly follow the free-field ground motion). SSI effects can be beneficial — reducing structural demands through foundation rocking — or detrimental — amplifying shaking in soft soils through Soil Amplification (Site Effect)The increase in shaking intensity caused by soft soil or sediment layers amplifying seismic waves. Structures built on soft soil can experience 2-10 times stronger shaking than those on bedrock. and Structural ResonanceThe amplification of building motion when earthquake wave frequency matches the building's natural frequency. Low-rise buildings resonate with high-frequency waves; tall buildings with low-frequency. effects.

Shake Table Physical Testing and Model Calibration

Simulation software must be validated against experimental data to be trustworthy. Shake table testing subjects physical structural models to earthquake ground motion recordings, generating measured response data that validates analytical models. Major shake table facilities include the NEES/NHERI E-Defense facility in Japan (the world's largest 3D shake table, capable of testing full-scale structures) and the facilities at UC San Diego, Lehigh University, and the University of Buffalo. Blind prediction contests, where analysts submit structural response predictions before the shake table test is run, rigorously evaluate simulation accuracy and reveal systematic model biases.

Monte Carlo Methods in Earthquake Risk

Regional Loss EstimationThe process of predicting economic losses and casualties from a potential earthquake scenario. FEMA's HAZUS software is the standard loss estimation tool in the United States. for insurance and financial risk management employs Monte Carlo simulation to characterize loss distributions. A stochastic event catalog is generated containing thousands of simulated years of earthquake activity, sampling from the Probabilistic Seismic Hazard Analysis (PSHA)A method for quantifying earthquake hazard that considers all possible earthquake sources, magnitudes, and ground motion levels, expressing results as probability of exceeding specific shaking levels. source model and Gutenberg-Richter LawA statistical law describing the relationship between earthquake frequency and magnitude: for each unit increase in magnitude, earthquakes become about 10 times less frequent. magnitude-frequency relationships. Each catalog event is combined with a ground motion simulation and a building damage model to produce a loss for that event. Aggregating over many thousands of events and catalog years produces exceedance probability curves for annual, aggregate, and occurrence losses — the basis for catastrophe insurance pricing and Probable Maximum Loss (PML)An estimate of the maximum loss an insurance portfolio or property is likely to experience from a single earthquake event. A key metric for insurers and reinsurers. calculations used in property financing.

High-Performance Computing Requirements

The most detailed earthquake simulations — physics-based ground motion simulation at regional scale with spatial resolution sufficient to capture basin effects — require high-performance computing (HPC) resources. Simulations of the 1994 Northridge and 1906 San Francisco earthquakes using Specfem3D have run on tens of thousands of processor cores simultaneously. As computational resources grow, simulation resolution improves, enabling direct coupling of seismic source physics, wave propagation, soil response, and structural analysis in a single integrated simulation — a long-term research goal that would dramatically improve earthquake risk assessment accuracy.

Summary

Earthquake simulation software spans scales from individual structural elements to entire building inventories across regional loss models. The Probabilistic Seismic Hazard Analysis (PSHA)A method for quantifying earthquake hazard that considers all possible earthquake sources, magnitudes, and ground motion levels, expressing results as probability of exceeding specific shaking levels. framework provides the probabilistic ground motion input; finite element codes simulate structural response; OpenQuake and HAZUS aggregate vulnerability functionsA mathematical function describing the probability of various damage states for a specific building type given a level of ground shaking. Essential for loss estimation models. with exposure to produce Loss EstimationThe process of predicting economic losses and casualties from a potential earthquake scenario. FEMA's HAZUS software is the standard loss estimation tool in the United States. products. These tools collectively enable quantitative earthquake risk management that would be impossible through purely empirical or judgment-based approaches.