Abstract — Earthquake Resistant Building Design
Civil and Environmental Engineering. GEM implemented the scientific methodology developed by UCLA [BDL+15] and integrated it into the QGIS framework as part of the OpenQuake IRMT plugin. This computational tool was developed to quantify the effectiveness of specific resilience-building actions (preparedness and response) that would enhance the process of recovery following an earthquake. Building performance is assessed using metrics that are explicitly linked to community seismic resilience such as functionality, occupiability and repairability. A probabilistic approach to mapping the fragility function parameters for the loss-based damage states used in OpenQuake and HAZUS (e.g. slight, moderate, extensive, complete) to the recovery-based damage states (e.g. functional loss, uninhabitable, irreparable) that are needed to assess seismic resilience is applied. A discrete state-based stochastic simulation model is used to quantify recovery trajectories for damaged buildings. Deterministic recovery paths conditioned on the post-earthquake damage state are described by discrete recovery states and the time spent in each state. The uncertainty in the recovery trajectory is captured by generating multiple realizations of the recovery path (Monte Carlo Simulation) accounting for the uncertainty in the time spent in each state. The function that describes community-level (or neighborhood-level) recovery is obtained by aggregating the recovery curves for the individual buildings after accounting for the variation in shaking and damage at each site and various exogenous and endogenous factors that affect the pace of recovery.
Analytical & Numerical Modeling Methodology
This tool provides an effective way to engage key stakeholders from government, industry and other groups, who are concerned with enhancing post-disaster recovery of communities. For example, focus groups with practitioners and emergency managers can be held to demonstrate its capabilities and explore their questions and concerns with regards to post-earthquake recovery. By allowing them to explore various trajectories for specific earthquake scenarios, the simulator will assist policy-makers, municipal governments and planners in understanding the possible interventions and tradeoffs associated with post-earthquake recovery. The simulation model can also be used to investigate hypotheses derived from findings and conclusions from previous earthquakes. Hypothesis examples include the following: designing for higher building performance levels prevents outmigration from the city, and âIIresidential red zonesâII or temporary building moratoriums are preferable to patchwork reconstruction.
Finite Element Results & Code Compliance
The recovery framework previously described was incorporated into this plugin, resulting in a dynamic and user-friendly software tool for generating both building-by-building and community level post-earthquake recovery predictions. The OpenQuake IRMT allows users to: develop composite indicators (or indices) to measure social vulnerability and/or disaster resilience; integrate these indices with physical risk estimates from OpenQuake or other software platforms; and In the following steps, a brief description of the basic workflow to develop an end-to-end recovery prediction is presented, where the main features and capabilities of the tool are highlighted.
Scope & Practical Applications
- High-rise building and spatial frame analysis under seismic & wind loading
- Finite element continuum modeling in ABAQUS for non-linear material plasticity
- Classical lamination theory (CLT) and failure prediction for composite structures
- Spatial hazard mapping, landslide susceptibility, and catchment terrain analysis in QGIS
- Detailed bar bending schedules (BBS), structural drawings, and calculation booklets
3D Frame Analysis
Spatial modeling in ETABS & STAAD.Pro
Non-Linear FEA
Constitutive damage models in ABAQUS
QGIS Geomatics
DEM terrain modeling & catchment maps
Stiffness Solver
Direct matrix assembly in MATLAB
Code Compliance
IS 456, IS 1893, ACI 318, Eurocodes
Stress Envelopes
Critical axial, shear & moment vectors
Engineering Software & Technologies
Built using industry-standard engineering analysis platforms, numerical code solvers, and geomatics tools.
| Software Platform | Analytical Function & Role | Engineering Tier |
|---|---|---|
| ETABS v21 | Multi-story building analysis, response spectrum, P-Delta, wind & seismic lateral loads | Structural Core |
| STAAD.Pro Connect | 3D space frame analysis, steel connection design, foundation reaction extraction | Steel & Concrete |
| ABAQUS Standard | Non-linear continuum mechanics, concrete damaged plasticity, composite lamina failure | Non-Linear FEA |
| QGIS Desktop | Digital Elevation Model (DEM) hydrological mapping, slope stability raster overlays | Spatial Geomatics |
| MATLAB R2024b | Global stiffness matrix assembly, modal eigenvalue extraction, composite ABD matrix calculation | Mathematical Solver |
Analytical Engineering Workflow
Systematic integration from CAD architectural geometry and GIS terrain mapping to finite element simulation and code validation.
Geometry & BIM
(Concrete/Steel/FRP)
Spatial Analysis
Non-Linear Mesh
& Foundation Coupling
& Detailing
Verified Structural Diagrams & FEA Stress Contours
Advantages of the Earthquake Resistant Building Design Project
Why this project delivers superior academic grades and provides high technical evaluation standing.
Multi-Software Cross-Validation
Integrates ETABS, STAAD.Pro, and ABAQUS modeling to verify results across both linear and non-linear finite element solvers.
Rigorous Code Compliance
Fully compliant with IS 456, IS 1893 (Seismic Criteria), IS 800 (Steel Design), ACI 318, and Eurocode provisions.
Non-Linear Material Plasticity
Captures crack propagation, yield surfaces, and post-peak softening behavior using ABAQUS Concrete Damaged Plasticity.
Composite Lamination Analysis
Evaluates Classical Lamination Theory (CLT), ABD matrices, and Hashin/Tsai-Wu failure envelopes for high-strength FRP laminates.
Spatial Geomatics Coupling
Combines high-resolution QGIS Digital Elevation Models with slope stability analysis for accurate terrain hazard zoning.
Deterministic Matrix Solvers
Custom MATLAB scripts assemble global stiffness matrices [K] and extract natural frequencies without relying purely on black-box GUIs.
Second-Order P-Delta Verification
Accurately assesses geometric non-linearity and sway amplification induced by gravity and lateral wind/seismic combinations.
Detailed Structural Detailing
Delivers complete Bar Bending Schedules (BBS), ductile detailing drawings per IS 13920, and steel connection drawings.
Direct Industry Relevance
Builds competencies demanded by top structural design consultancies, EPC contractors, and infrastructure developers.
Complete Turnkey Package
Includes full model files (.EDB, .STD, .CAE), MATLAB code, calculation sheets, IEEE-format report, and viva defense presentation.
Sample MATLAB Code, FEA Output & Member Schedules
Representative excerpts illustrating the finite element stiffness assembly, command solver output, and structural capacity envelopes.
[*] Initializing ETABS API Engine v21 & STAAD.Pro Analytical Space Frame... [✓] Dead Load (DL), Live Load (LL), Wind (WL), and Seismic (EQX/EQY) Loaded ──────────────────────────────────────────────────────────────────────── [MODAL ANALYSIS] Mode 1 (Translational X): T1 = 0.842 sec (Mass Part: 81.4%) [MODAL ANALYSIS] Mode 2 (Translational Y): T2 = 0.791 sec (Mass Part: 82.8%) [MODAL ANALYSIS] Mode 3 (Torsional Theta): T3 = 0.512 sec (T/Tx = 0.608 < 0.90 OK) [DRIFT CHECK] Max Inter-Story Drift: Δ = 0.0018 < 0.0040 (IS 1893 Clause 7.11.1 PASSED) [P-DELTA EFFECT] Maximum Stability Index θ: 0.041 < 0.10 (Secondary Moments Negligible) ──────────────────────────────────────────────────────────────────────── [✓] Design Envelope Complete: Zero Overstressed Members across 1.5(DL + LL) & 1.2(DL + LL ± EL)
[*] Generating critical member capacity envelope across factored combinations...
| Member Tag | Section Designation | Axial Pu (kN) | Bending Mz (kNm) | Shear Vy (kN) | Capacity Ratio | IS Code Status |
|---|---|---|---|---|---|---|
| Col C1-G1 | 450 x 450 mm M40 | 2,410.8 | 184.2 | 84.5 | 0.742 | PASSED |
| Col C2-G1 | 500 x 500 mm M40 | 3,180.4 | 218.6 | 112.4 | 0.814 | PASSED |
| Beam B101 | 300 x 600 mm M35 | 42.5 | 264.8 | 168.2 | 0.688 | PASSED |
| Beam B102 | 300 x 600 mm M35 | 38.1 | 295.4 | 182.1 | 0.762 | PASSED |
[✓] Structural Analysis Complete: Zero overstressed members. Design capacity factor > 1.0 throughout.
Frequently Asked Questions
Project Deliverables & Support
Everything needed from project inception to final viva defense — delivered with expert technical support.
Model Files
Input Decks
Scripts
Report
Slides
Drawings
Sheets
Q&A
Support
Extension